Modeling method of tropospheric delay for arid-areas including sandy, gobi and desert regions based on GNSS

Patent metadata

Field	Value
Publication number	US-12517261-B2
Application number	US-202519083465-A
Country	US
Kind code	B2
Filing date	Mar 19, 2025
Priority date	Mar 22, 2024
Publication date	Jan 6, 2026
Grant date	Jan 6, 2026

Abstract

The provided is a modeling method of tropospheric delay for arid-areas including sandy, gobi and desert regions based on GNSS. The modeling method incudes: obtaining a troposphere zenith hydrostatic delay ZHD and a troposphere zenith wet delay ZWD by utilizing a GPT-3 model, and summing to obtain a troposphere zenith total delay ZTD; inversely calculating the ZTD by utilizing a ground-based GNSS, obtaining a high-precision ZWD inversed by the GNSS by calculating a difference between the ZTD and the troposphere zenith hydrostatic delay ZHD obtained by a Saastamonien model, inversely calculating an atmospheric weighted average temperature T′m by utilizing an ERA5 data set, performing linear correction to obtain a corrected temperature Tmr, and calculating a priori value ZWD0 of the high-precision ZWD; and correcting the ZWD0, and adding the corrected ZWD0 with the ZHD to obtain a corrected value of the high-precision ZTD.

First claim

The invention claimed is:

1 . A modeling method of tropospheric delay for arid-areas including sandy, gobi and desert regions based on GNSS, comprising:
step 1: obtaining an atmospheric pressure P, an atmospheric weighted average temperature Tm, a water vapor pressure lapse rate λm and a water vapor pressure e according to a geographical position and time of a station by utilizing a global pressure and temperature model; obtaining a troposphere zenith hydrostatic delay ZHD and a troposphere zenith wet delay ZWD according to meteorological parameters obtained by the global pressure and temperature model, and obtaining a troposphere zenith total delay ZTD by adding the ZHD with the ZWD; step 2: inversely calculating the ZTD by utilizing a ground-based GNSS, obtaining the troposphere zenith hydrostatic delay ZHD by inputting an air pressure in a fifth generation of atmospheric reanalysis data set of global climate into a Saastamonien model, and obtaining a troposphere zenith wet delay ZWD inversely calculated according to the GNSS by calculating a difference between the ZTD and the ZHD; step 3: inversely calculating the atmospheric weighted average temperature T′m according to a water pressure and a temperature in the fifth generation of atmospheric reanalysis data set of global climate, and obtaining a corrected temperature T m r by linearly correcting the Tm calculated by the global pressure and temperature model in the step 1 utilizing the T′m; step 4: calculating a priori value ZWD0 of the troposphere zenith wet delay utilizing an Askne model according to the T m r obtained in the step 3; step 5: establishing a correction model of troposphere zenith wet delay by utilizing a bidirectional long and short term memory network BiLSTM, taking the priori value ZWD0 of the troposphere zenith wet delay obtained in the step 4 as input data of the correction model, taking the troposphere zenith wet delay ZWD inversely calculated by the GNSS in the step 2 as reference data for a corrected value ZWDm of the troposphere zenith wet delay which is an output result of the correction model, establishing a training data set of the correction model of troposphere zenith wet delay, and setting a loss function to continuously optimize the BiLSTM; and step 6: completing establishment of the troposphere zenith total delay ZTD correction model in the arid region and calculating a corrected troposphere zenith total delay ZTD by summing the corrected value ZWDm of the troposphere zenith wet delay obtained by the correction model of troposphere zenith wet delay in the step 5 and the ZHD calculated by the global pressure and temperature model in the step 1; selecting multiple IGS stations in the arid region, which cover different regions in the northern and southern hemispheres; calculating values of a mean absolute error MAE and a root mean square error RMS of the troposphere zenith total delay ZTD of the IGS stations according to the corrected troposphere zenith total delay ZTD; wherein the global pressure and temperature model is used for: calculating the atmospheric pressure P, the atmospheric weighted average temperature Tm, the water vapor pressure lapse rate λ m and the water vapor pressure e according to the geographical position and the time of the station,

r
⁡
(

P
,

T
⁢
m

,

λ
m

,
e

)

=

A
0

+

A
1

⁢
sin
⁢

(

doy

3
⁢
6
⁢

5
.
2

⁢
5

⁢

2
⁢
π

)

+

B
1

⁢
cos
⁢

(

doy

3
⁢
6
⁢

5
.
2

⁢
5

⁢
2
⁢
π

)

+

A
2

⁢
sin
⁢

(

doy

3
⁢
6
⁢

5
.
2

⁢
5

⁢
4
⁢
π

)

+

B
2

⁢
cos
⁢

(

doy

3
⁢
6
⁢

5
.
2

⁢
5

⁢
4
⁢
π

)

wherein, A 0 is a average value calculated by a least square method, A 1 and B 1 are both annual period parameters, A 2 and B 2 are both sem

Assignees

• Three Gorges Group Ind Development Beijing Co Ltd
• China Three Gorges Corp

Inventors

• Lv Tingyan
• Zhou Jingliang
• Zhang Rui
• Chen Rongmin
• Liu Jiaqi
• Qi Zhuoya

Classifications

• G01S19/40

  Correcting position, velocity or attitude · CPC title

• Y02A90/10

  Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation · CPC title

• G06F2113/08

  Fluids · CPC title

• G06F2119/14

  Force analysis or force optimisation, e.g. static or dynamic forces · CPC title

• G06F2119/08

  Thermal analysis or thermal optimisation · CPC title