Weather data assimilation support device, weather data assimilation support method, program, and weather data assimilation support system

Abstract

This weather data assimilation assistance device (1) comprises: a wet delay amount calculation unit (12) for calculating a wet delay amount in a satellite line-of-sight direction on the basis of observation amount data indicating an observation amount of radio waves from a positioning satellite as measured by a receiver that receives radio waves from the positioning satellite, and atmospheric delay amount data indicating an atmospheric delay amount in a satellite line-of-sight direction calculated on the basis of a correction value for the satellite clock and / or satellite orbit of the positioning satellite, and outputting data indicating the wet delay amount; and a wet delay amount variance calculation unit (13) for calculating a variance of the wet delay amount, on the basis of the wet delay amount and the covariance matrix of the correction value for the satellite clock and / or satellite orbit of the positioning satellite, and outputting data indicating the variance of the wet delay amount.

Description

[Technical Field]

[0001] The present disclosure relates to a meteorological data assimilation support device, a meteorological data assimilation support method, a program, and a meteorological data assimilation support system. [Background technology]

[0002] The atmosphere is one of the error sources that causes delays in signals from positioning satellites, such as the Global Positioning System (GPS). By utilizing this characteristic, the atmospheric state can be determined by receiving signals from positioning satellites on the ground and calculating the delay caused by the atmosphere (hereinafter referred to as atmospheric delay). By subtracting the delay in dry conditions (hereinafter referred to as dry delay) from the atmospheric delay, the delay caused by water vapor in the atmosphere (hereinafter referred to as wet delay) can be calculated. This makes it possible to estimate the amount of water vapor in the atmosphere, and assimilating the amount of water vapor in the atmosphere with other meteorological data is expected to improve the accuracy of numerical weather forecasts.

[0003] The estimated atmospheric delay on the signal path connecting the positioning satellite and the receiver contains errors caused by factors other than the atmosphere, such as satellite clock and orbit correction errors. Therefore, a method has been proposed to cancel out errors caused by factors other than the atmosphere by taking the difference between the atmospheric delay estimates from multiple receivers.

[0004] For example, Patent Document 1 discloses an atmospheric-related quantity derivation device that acquires a first phase difference based on the phase of radio waves transmitted from a transmitting station and received by a receiver at a first location, and a second phase difference based on the phase of radio waves transmitted from the transmitting station and received by a receiver at a second location, calculates the difference between the first phase difference and the second phase difference, and derives a relative related quantity, which is the difference between a first atmospheric-related quantity between the transmitting station and the first location and a second atmospheric-related quantity between the transmitting station and the second location, based on the calculated difference. [Prior art documents] [Patent documents]

[0005] [Patent Document 1] Japanese Patent Application Laid-Open No. 2017-207459 Summary of the Invention [Problem to be solved by the invention]

[0006] The technology described in Patent Document 1 is not effective when applied to a wide area when the distance between receivers is long, when there are no other receivers in the vicinity, etc. If the error cannot be canceled out, when assimilating with other meteorological data, it is necessary to correctly estimate the degree of error contained in the estimated value of the wet delay, and if the error is estimated to be large, to reduce the assimilation weight and perform assimilation by assigning a weight according to the degree of error.

[0007] The present disclosure has been made in consideration of the above-mentioned problems, and aims to make it possible to more accurately estimate the degree of error contained in the estimated value of the amount of wet delay. [Means for solving the problem]

[0008] To achieve the above object, a meteorological data assimilation support device according to the present disclosure includes a wet delay calculation unit and a wet delay variance calculation unit. The wet delay calculation unit calculates a wet delay in the satellite line of sight based on observation data indicating observation amounts of radio waves from positioning satellites measured by a receiver that receives radio waves from the positioning satellites, and atmospheric delay data indicating atmospheric delay in the satellite line of sight calculated based on correction values ​​for at least one of the satellite clock and the satellite orbit of the positioning satellite, and outputs data indicating the wet delay. The wet delay variance calculation unit calculates a variance of the wet delay based on a covariance matrix of the wet delay and the correction values ​​for at least one of the satellite clock and the satellite orbit of the positioning satellite, and outputs data indicating the variance of the wet delay. [Effects of the Invention]

[0009] According to the present disclosure, by using the covariance matrix of the correction values ​​of at least one of the satellite clock and satellite orbit of the positioning satellite to calculate the variance of the wet delay amount in the satellite line of sight, and providing the wet delay amount in the satellite line of sight and its variance as wet delay amount data, it becomes possible to more accurately estimate the amount of error contained in the estimated value of the wet delay amount. [Brief explanation of the drawings]

[0010] [Figure 1] FIG. 1 is a block diagram showing an example of the functional configuration of a meteorological data assimilation support device according to a first embodiment. [Figure 2] 1 is a flowchart showing a meteorological data assimilation support process according to the first embodiment. [Figure 3] A block diagram showing an example of the functional configuration of a meteorological data assimilation support device according to a second embodiment. [Figure 4] FIG. 10 is a block diagram showing a configuration example of a meteorological data assimilation support system according to a modification of the second embodiment. [Figure 5] A block diagram showing an example of the functional configuration of a meteorological data assimilation support device according to a third embodiment. [Figure 6] A block diagram showing an example of the functional configuration of a meteorological data assimilation support device according to a fourth embodiment. [Figure 7] FIG. 1 is a diagram showing an example of the hardware configuration of a meteorological data assimilation support device according to first to fourth embodiments. [Figure 8] FIG. 1 is a block diagram showing an example of the functional configuration of a meteorological data assimilation support device according to a modification of the first embodiment. [Figure 9] FIG. 10 is a block diagram showing an example of the functional configuration of a meteorological data assimilation support device according to a modification of the second embodiment. [Figure 10] FIG. 10 is a block diagram showing an example of the functional configuration of a meteorological data assimilation support device according to a modification of the third embodiment. [Figure 11] FIG. 10 is a block diagram showing an example of the functional configuration of a meteorological data assimilation support device according to a modification of the fourth embodiment. DETAILED DESCRIPTION OF THE INVENTION

[0011] A meteorological data assimilation support device, meteorological data assimilation support method, program, and meteorological data assimilation support system according to the present embodiment will be described in detail below with reference to the drawings. Note that identical or corresponding parts in the drawings are designated by the same reference numerals. In this embodiment, an example will be described in which the atmospheric delay in the satellite line of sight is calculated based on the observed amount of radio waves from a positioning satellite and correction information for the satellite clock and satellite orbit, and the moist delay in the satellite line of sight and the variance of the moist delay are calculated based on the atmospheric delay in the satellite line of sight and the correction information for the satellite clock and satellite orbit. Hereinafter, a positioning satellite will be simply referred to as a satellite.

[0012] (Embodiment 1) As shown in Figure 1, the meteorological data assimilation support system 100 includes a receiver i that receives radio waves from a satellite l, an atmospheric delay calculation device 2 that calculates the atmospheric delay in the satellite line of sight based on the radio wave observations from satellite l measured by receiver i and correction information for the satellite clock and satellite orbit of satellite l, and a meteorological data assimilation support device 1 that calculates the moist delay in the satellite line of sight based on the atmospheric delay in the satellite line of sight calculated by atmospheric delay calculation device 2, calculates the variance of the moist delay in the satellite line of sight based on the covariance matrix of the moist delay in the satellite line of sight and the correction values ​​for the satellite clock and satellite orbit, and outputs moist delay data including data indicating the moist delay in the satellite line of sight and data indicating the variance of the moist delay in the satellite line of sight. Although one receiver i is shown as a representative example, multiple receivers i may be used.

[0013] Receiver i receives radio waves from satellite l, measures the pseudorange and carrier phase, and transmits observation data indicating the observation amounts to atmospheric delay calculation device 2. Note that if only one of the observation amounts, pseudorange and carrier phase, can be obtained, only one of them may be used.

[0014] The atmospheric delay calculation device 2 includes an observation quantity acquisition unit 21 that acquires observation quantity data from the receiver i, a correction information acquisition unit 22 that acquires correction information for the satellite clock and satellite orbit, and an atmospheric delay calculation unit 23 that calculates the atmospheric delay in the satellite line of sight based on the observation quantity data and the correction information for the satellite clock and satellite orbit. Hereinafter, the atmospheric delay in the satellite line of sight will be simply referred to as the atmospheric delay. The correction information for the satellite clock and satellite orbit acquired by the correction information acquisition unit 22 includes correction values ​​for the satellite clock and satellite orbit and the covariance matrix of the correction values. The correction information acquisition unit 22 sends the correction values ​​for the satellite clock and satellite orbit to the atmospheric delay calculation unit 23.

[0015] A reference for the technology of calculating atmospheric delay is, for example, "GPS Meteorology: Research on the construction of a GPS water vapor information system and its application to meteorology, geodesy, and hydrology" by Koji Yoshinori et al., Journal of the Geodetic Society, Vol. 55, No. 1 (2009).

[0016] Here, an example will be described in which the atmospheric delay calculation unit 23 calculates the atmospheric delay using the observation model described in the above-mentioned reference. i l is the observed carrier phase of the radio wave emitted from satellite l measured by receiver i. On the right side, the atmospheric delay dT i l Includes:

[0017] [Number 1] Φ i l =ρ i l +N i l λ+dT i l -dI i l +c(δ i -δ l )+ε i

[0018] ρ i l is the distance between the satellite and the receiver, N i lis the carrier phase ambiguity, λ is the wavelength of the radio wave, and dT i l is the atmospheric delay, dI i l is the ionospheric delay, c is the speed of light in a vacuum, δ l is the satellite clock, δ i is the receiver clock, ε i is the residual. Atmospheric delay dT i l can be replaced by the model shown in the following equation 2.

[0019] [Number 2] dT l i =m(θ l )·[ZTD i +cotθ l (G ni cosφ l +G ei sinφ l )]

[0020] The atmospheric delay calculation unit 23 calculates the zenith atmospheric delay ZTD i and the linear gradients of atmospheric delay in the north-south and east-west directions, respectively, G ni and G ei These three are estimated as unknown parameters. ni and G ei is not estimated and can be omitted. l ) is a mapping function that is the ratio of the zenith atmospheric delay to the atmospheric delay. θ l is the satellite elevation angle, φ l is the azimuth angle of the satellite, which is calculated based on the approximate position of the receiver and the approximate position of the satellite. l and φ l is calculated, for example, from the broadcast ephemeris delivered by the satellite.

[0021] As shown in the following equation (3), the estimated values ​​of these unknown parameters are calculated using, for example, a Kalman filter, and the a posteriori residual dz i l Using this, we can obtain an estimate of the atmospheric delay. ni and G eiIf no estimate is made, it will also be excluded here.

[0022]

number

[0023] In order to distinguish Equation 3 from Equation 2, a hat is added to indicate that it is an estimated value. i l Since the atmospheric delay estimates are corrected by adding a dash, the corrected values ​​are indicated by the dash. i l includes correction errors of the satellite clock and the satellite orbit. The atmospheric delay calculation unit 23 calculates the a posteriori residual dz i l Calculate.

[0024] The posterior residual is dz i l , can be expressed by the following formula 4. In formula 4, the estimated value by the Kalman filter is indicated by a hat. i l and δ l are marked with a dash to indicate that they have been corrected using satellite clock and orbit corrections.

[0025]

number

[0026] The satellite clock is δ l The satellite orbit is the distance ρ between the satellite and the receiver. i l This corresponds to the satellite position for calculating the carrier phase ambiguity N i l , ionospheric delay dI i l , and the receiver clock δ i At least N of il , δ i is an unknown parameter, and is estimated using a Kalman filter together with an estimated atmospheric delay. i l The receiver position for calculating dI can be either known in advance or an unknown parameter. i l is sometimes given as correction information for ionospheric delay, sometimes estimated, and sometimes eliminated using two-frequency observations, as described below.

[0027] The atmospheric delay calculation unit 23 transmits atmospheric delay data indicating the calculated estimated value of the atmospheric delay to the meteorological data assimilation support device 1. In addition, the correction information acquisition unit 22 transmits to the meteorological data assimilation support device 1 the covariance matrix of the correction values ​​for the satellite clock and satellite orbit included in the acquired correction information for the satellite clock and satellite orbit.

[0028] The meteorological data assimilation support device 1 includes a pressure information acquisition unit 11 that acquires pressure information indicating the pressure at the receiver i, a wet delay calculation unit 12 that calculates the wet delay in the satellite line of sight, and a wet delay variance calculation unit 13 that calculates the variance of the wet delay in the satellite line of sight based on the wet delay in the satellite line of sight and the covariance matrix of the correction values ​​for the satellite clock and satellite orbit. Hereinafter, the wet delay in the satellite line of sight will be simply referred to as the wet delay.

[0029] The atmospheric pressure information acquisition unit 11 sends the acquired atmospheric pressure information to the wet delay amount calculation unit 12 and the wet delay amount variance calculation unit 13. The atmospheric pressure information acquisition unit 11 may acquire the atmospheric pressure information based on the latitude and date data of the receiver i, or may acquire the information from an atmospheric pressure sensor provided in the receiver i.

[0030] The wet delay calculation unit 12 calculates the wet delay based on the atmospheric delay data received from the atmospheric delay calculation unit 23 and the atmospheric pressure information received from the atmospheric pressure information acquisition unit 11. There are several methods for calculating the wet delay, such as a method that uses only atmospheric pressure information or a method that uses a numerical weather model, but in this embodiment, an example will be described in which a method that uses only atmospheric pressure information is adopted. For example, a reference for a calculation method that uses only atmospheric pressure information is Pratap Misra and Per Enge, "Detailed Explanation of GPS Basic Concepts, Positioning Processing, Signals and Receivers," 1st Edition, translated by the GPS Study Group of the Japan Institute of Navigation.

[0031] When the wet delay calculation unit 12 calculates the wet delay using the calculation method described in the above-mentioned reference, first, the atmospheric pressure Pi at the receiver i indicated by the atmospheric pressure information is used, and the zenith dry delay ZHDmodel, which is the dry delay in the zenith direction, is calculated by the following formula 5: i (P i ) to get lat i is the latitude of receiver i, H i is the antenna height above sea level.

[0032] [Number 5] ZHDmodel i (P i )=0.002277(1+0.0026cos2lat i +0.00028H i )P i

[0033] As an example, a mapping function m(θ l ) is multiplied to obtain the drying delay in the line of sight dThydrostatic i l θ l is the elevation angle of the satellite. Hereinafter, the drying delay amount in the line of sight direction will be simply referred to as the drying delay amount.

[0034] [Number 6] m(θ l )=1 / (sinθ l +0.00143 / tanθ l +0.0445)

[0035] [Number 7] dThydrostatic i l =m(θ l )·ZHDmodel i (P i )

[0036] Many models have been proposed for the zenith dry delay amount and mapping function, and the above examples are the most representative ones, but are not limited to these.

[0037] The wet delay amount calculation unit 12 calculates the wet delay amount by subtracting the dry delay amount from the atmospheric delay amount. The wet delay amount calculation unit 12 outputs data indicating the calculated wet delay amount as wet delay amount data. The wet delay amount calculation unit 12 also sends data indicating the calculated wet delay amount to the wet delay amount variance calculation unit 13.

[0038] The wet delay amount variance calculation unit 13 calculates the variance of the wet delay amount based on the atmospheric pressure information received from the atmospheric pressure information acquisition unit 11 and the data indicating the wet delay amount received from the wet delay amount calculation unit 12, using the covariance matrix of the correction values ​​of the satellite clock and satellite orbit contained in the correction information received from the correction information acquisition unit 22 as an index representing the correction error of the satellite clock and satellite orbit.

[0039] First, the amount of wet retardation can be calculated, for example, by the following formula (8).

[0040]

number

[0041] Therefore, the wet delay dTwet i l As an example, the variance of can be obtained by the following formula (9).

[0042]

number

[0043] ZWD i is an estimate of the zenith wet delay. The zenith wet delay is calculated by subtracting the zenith dry delay from the estimated zenith atmospheric delay. There are several methods for calculating the zenith dry delay, such as using only atmospheric pressure information or using a numerical weather model. When only atmospheric pressure information is used, the zenith wet delay can be obtained using the following formula:

[0044] [Number 10] ZWD i =ZTD i -ZWDmodel i (P i )

[0045] σ ZWDi is the standard deviation of the zenith wet delay, and is calculated, for example, by the following formula (11).

[0046]

number

[0047] σ ZTDi is the estimated zenith atmospheric delay ZTD i The standard deviation of σ is used, and for example, the value obtained when the atmospheric delay calculation unit 23 estimates it using a Kalman filter. ZHDmodeli is the standard deviation of the model error of the zenith dry delay, and for example, a value proportional to the amount of cumulonimbus clouds above receiver i is used. σ PR is the error standard deviation of the atmospheric pressure observation. σ ZTDi , σ ZHDmodeli and σ PR Any or all of the can be omitted.

[0048] σ Gni and σ Gei are the standard deviations of the estimated first-order gradients of atmospheric delay in the north-south and east-west directions, respectively. σ Gni and σ Gei Either or both of the σ can be omitted. ZTDi , σ Gni and σ Geiis included in the data indicating the amount of wet delay, and is included in the amount of atmospheric delay data.

[0049] The second term on the right side of Equation 9 is the correction error of the satellite clock and satellite orbit. This term is the covariance matrix P clk,orb l and the line of sight vector los i l It can be obtained by the following formula 12 using

[0050]

number

[0051] The wet delay amount variance calculation unit 13 outputs data indicating the variance of the calculated wet delay amount as wet delay amount data. The wet delay amount calculation unit 12 and the wet delay amount variance calculation unit 13 may, for example, display the wet delay amount data on a screen, transmit it to a user terminal, or transmit it to a meteorological data assimilation device. When transmitted to a meteorological data assimilation device, the meteorological data assimilation device assimilates meteorological data using the received wet delay amount data. By obtaining the wet delay amount data, the user or the meteorological data assimilation device can more accurately estimate the degree of error contained in the estimated value of the wet delay amount.

[0052] Here, the flow of the meteorological data assimilation support process executed by the meteorological data assimilation support device 1 will be described with reference to Figure 2. The meteorological data assimilation support process shown in Figure 2 starts, for example, when an instruction to generate wet delay data is input to the meteorological data assimilation support device 1. The wet delay calculation unit 12 calculates the zenith dry delay using the atmospheric pressure at receiver i indicated by the atmospheric pressure information acquired by the atmospheric pressure information acquisition unit 11 (step S11). The wet delay calculation unit 12 multiplies the zenith dry delay by a mapping function to calculate the dry delay, which is the dry delay in the line of sight direction (step S12). The wet delay calculation unit 12 calculates the wet delay by subtracting the dry delay from the atmospheric delay (step S13).

[0053] The wet delay amount calculation unit 12 outputs wet delay amount data including the calculated wet delay amount (step S14). In addition, the wet delay amount calculation unit 12 sends data indicating the calculated wet delay amount to the wet delay amount variance calculation unit 13.

[0054] The wet delay amount variance calculation unit 13 calculates the variance of the wet delay amount based on the atmospheric pressure information received from the atmospheric pressure information acquisition unit 11, the data indicating the wet delay amount received from the wet delay amount calculation unit 12, and the covariance matrix of the correction values ​​of the satellite clock and satellite orbit received from the correction information acquisition unit 22 of the atmospheric delay amount calculation device 2 (step S15). The wet delay amount variance calculation unit 13 outputs wet delay amount data including the variance of the calculated wet delay amount (step S16), and ends the process.

[0055] According to the meteorological data assimilation support device 1 of embodiment 1, the variance of the wet delay amount is calculated using the covariance matrix of the correction values ​​of the satellite clock and satellite orbit, and the wet delay amount and its variance are provided as wet delay amount data, making it possible to more accurately estimate the amount of error contained in the estimated value of the wet delay amount.

[0056] (Embodiment 2) In the second embodiment, the meteorological data assimilation support device 1 includes the configuration of the atmospheric delay calculation device 2 of the first embodiment, and the meteorological data assimilation support device 1 generates correction information for the satellite clock and satellite orbit internally.

[0057] A meteorological data assimilation support system 200 according to the second embodiment is shown in Fig. 3. Similar to the meteorological data assimilation support system 100, the meteorological data assimilation support system 200 includes a receiver i, an atmospheric delay calculation device 2, and a meteorological data assimilation support device 1. The functions of the receiver i are similar to those of the meteorological data assimilation support system 100.

[0058] The atmospheric delay calculation device 2 includes an observation quantity acquisition unit 21 and an atmospheric delay calculation unit 23. The meteorological data assimilation support device 1 includes a pressure information acquisition unit 11, a wet delay calculation unit 12, a wet delay variance calculation unit 13, and in addition, a satellite clock and orbit correction information generation unit 14 that generates correction information for the satellite clock and satellite orbit from the observation quantity data acquired by the observation quantity acquisition unit 21.

[0059] The observation quantity acquisition unit 21 sends the observation quantity data to the atmospheric delay calculation unit 23 and the satellite clock and orbit correction information generation unit 14 .

[0060] The satellite clock and orbit correction information generator 14 generates satellite clock and orbit correction information, including satellite clock and orbit correction values ​​and the covariance matrix of the correction values, based on the observation data acquired by the observation data acquisition unit. For example, the satellite clock and orbit correction information generator 14 treats the observation data, which indicates the observation values ​​measured by receivers located around the world receiving radio waves from each satellite, as time-series data and processes it using a Kalman filter. An example of using a Kalman filter is A. Rovira-Garcia et al., "Fast Precise Point Positioning: A System to Provide1 Correction for Single and Multi-frequency Navigation," Navigation, 63(3), 2016.

[0061] Here, we will explain an example in which the satellite clock and orbit correction information generator 14 generates a covariance matrix using the Kalman filter described in the above-mentioned reference. In the above-mentioned Kalman filter, in addition to the clock and orbit errors (correction information) of each satellite, the zenith tropospheric delay for each receiver, tropospheric tilt, receiver clock error, and carrier phase bias for each receiver and satellite (ionospheric-free coupling of two-frequency ambiguity) are simultaneously estimated as state quantities. The satellite clock and orbit correction information generator 14 extracts terms related to the clock and orbit errors of each satellite from the covariance matrix of these state quantities held by the Kalman filter.

[0062] As an example, an observation equation using the observation amount measured by receiver i receiving radio waves from satellite l is shown in Equation 13. In the above-mentioned Kalman filter, simultaneous observation equations are set up for the number of combinations of receivers and satellites that have acquired the observation amount, and the state amount is estimated.

[0063]

number

[0064] In the example of Equation 13, two frequency signals s1 and s2 are used to eliminate the influence of ionospheric delays, and ionospheric free coupling is used. f is the signal frequency, Φ is the carrier phase observation, and P is the pseudorange observation. x i is the known coordinate value of receiver i, x l is the reference position of satellite l calculated from the broadcast ephemeris, the ultra-fast ephemeris, the fast ephemeris, etc., and δx l is the error of the reference position relative to the true position. i l is the line of sight vector. l is the clock error of receiver i, δclk l is the clock error of satellite l. dT i l is the atmospheric delay, which can be replaced by the model shown in Equation 2 above. i l is the carrier phase bias, ε is the observation error, and IF (Ionosphere Free) refers to ionosphere free.

[0065] The covariance matrix of the satellite clock and satellite orbit correction value of satellite l is the covariance matrix of the total state quantity, l and δx l The satellite clock and orbit correction information generation unit 14 sends the generated satellite clock and satellite orbit correction values ​​to the atmospheric delay calculation unit 23 of the atmospheric delay calculation device 2, and sends the covariance matrix of the generated satellite clock and satellite orbit correction values ​​to the wet delay variance calculation unit 13. Other functional configurations are the same as those in the first embodiment.

[0066] The atmospheric delay calculation device 2 may be configured to include a satellite clock and orbit correction information generation unit 14. An example of the functional configuration of the atmospheric delay calculation device 2 and the meteorological data assimilation support device 1 of the meteorological data assimilation support system 200' according to this modification is shown in FIG. 4. Furthermore, the moist delay calculation unit 12 does not necessarily generate moist delay data for all receivers i located around the globe. In other words, the receiver i whose observation data is used by the satellite clock and orbit correction information generation unit 14 is not necessarily the same as the receiver i whose observation data is used by the atmospheric delay calculation unit 23.

[0067] The meteorological data assimilation support device 1 according to the second embodiment calculates the variance of the wet delay using the covariance matrix of the correction values ​​for the satellite clock and satellite orbit, and provides the wet delay and its variance as wet delay data, thereby enabling a more accurate estimation of the degree of error contained in the estimated value of the wet delay. Furthermore, by internally generating the correction information for the satellite clock and satellite orbit, the meteorological data assimilation support device 1 can be implemented even when the correction information for the satellite clock and satellite orbit cannot be obtained externally.

[0068] (Embodiment 3) In embodiment 3, in addition to the configuration of embodiment 1, the amount of cumulonimbus clouds above receiver i is determined, and the error variance of the dry delay amount is calculated to be used in calculating the variance of the wet delay amount depending on the amount of cumulonimbus clouds above receiver i.

[0069] A meteorological data assimilation support system 300 according to the third embodiment is shown in Fig. 5. Similar to the meteorological data assimilation support system 100, the meteorological data assimilation support system 300 includes a receiver i, an atmospheric delay calculation device 2, and a meteorological data assimilation support device 1. The functions of the receiver i and the atmospheric delay calculation device 2 are similar to those of the meteorological data assimilation support system 100.

[0070] The meteorological data assimilation support device 1 includes a pressure information acquisition unit 11, a wet delay calculation unit 12, a wet delay variance calculation unit 13, and a cumulonimbus determination unit 15 that determines the amount of cumulonimbus clouds above the receiver i.

[0071] The cumulonimbus cloud determination unit 15 determines the amount of cumulonimbus clouds above the receiver i using inputs such as weather radar images, which are images from a weather radar, and temperature and atmospheric pressure dynamics information indicating temperature and atmospheric pressure perturbations in other weather models. The cumulonimbus cloud determination unit 15 calculates an error variance of the dry delay amount, which increases as the amount of cumulonimbus clouds above the receiver i increases. For example, the cumulonimbus cloud determination unit 15 calculates an error variance of the dry delay amount that is directly proportional to the amount of cumulonimbus clouds above the receiver i. The cumulonimbus cloud determination unit 15 sends data indicating the calculated error variance of the dry delay amount to the wet delay amount variance calculation unit 13. The cumulonimbus cloud determination unit 15 may also use as input the variance of atmospheric pressure indicated by the atmospheric pressure information acquired by the atmospheric pressure information acquisition unit 11 to determine the amount of cumulonimbus clouds above the receiver i.

[0072] The wet delay amount variance calculation unit 13 calculates the variance of the wet delay amount based on the atmospheric pressure information received from the atmospheric pressure information acquisition unit 11, data indicating the wet delay amount received from the wet delay amount calculation unit 12, the covariance matrix of the satellite clock and satellite orbit correction values ​​received from the satellite clock and orbit correction information generation unit 14, and data indicating the error variance of the dry delay amount received from the cumulonimbus cloud determination unit 15. Other functional configurations are the same as those in the first embodiment.

[0073] According to the meteorological data assimilation support device 1 of the third embodiment, the variance of the wet delay is calculated using the covariance matrix of the correction values ​​of the satellite clock and satellite orbit, and the wet delay and its variance are provided as wet delay data, thereby making it possible to more accurately estimate the degree of error contained in the estimated value of the wet delay. Furthermore, by taking into account the amount of cumulonimbus clouds above the receiver i, the estimation accuracy of the wet delay variance calculated by the wet delay variance calculation unit 13 can be improved.

[0074] (Fourth embodiment) In embodiment 4, in addition to the configuration of embodiment 1, the amount of cumulonimbus clouds above receiver i is determined, and depending on the result of determining the amount of cumulonimbus clouds above receiver i, the method for calculating the amount of wet delay is switched between a method that uses only atmospheric pressure information and a method that uses a numerical weather model.

[0075] A meteorological data assimilation support system 400 according to the fourth embodiment is shown in Fig. 6. Similar to the meteorological data assimilation support systems 100 and 300, the meteorological data assimilation support system 400 includes a receiver i, an atmospheric delay calculation device 2, and a meteorological data assimilation support device 1. The functions of the receiver i and the atmospheric delay calculation device 2 are similar to those of the meteorological data assimilation support systems 100 and 300.

[0076] The meteorological data assimilation support device 1 includes a pressure information acquisition unit 11, a wet delay amount calculation unit 12, a wet delay amount variance calculation unit 13, a cumulonimbus cloud determination unit 15, and a numerical meteorological data acquisition unit 16 that acquires numerical meteorological data related to temperature and pressure.

[0077] The numerical weather data acquisition unit 16 sends the acquired numerical weather data related to the temperature and atmospheric pressure to the wet delay amount calculation unit 12. The wet delay amount calculation unit 12 can switch between a method of using only atmospheric pressure information and a method of using a numerical weather model based on numerical weather data as a method of calculating the dry delay amount for calculating the wet delay amount. When calculating the dry delay amount using a method of using a numerical weather model, the wet delay amount calculation unit 12 calculates the dry delay amount using a pre-stored numerical weather model based on the numerical weather data received from the numerical weather data acquisition unit 16.

[0078] The cumulonimbus cloud determination unit 15 determines the amount of cumulonimbus clouds above the receiver i using, as input, weather radar images, which are images from a weather radar, and temperature and atmospheric pressure dynamics information indicating temperature and atmospheric pressure perturbations in a weather model. The cumulonimbus cloud determination unit 15 sends an instruction to the wet delay amount calculation unit 12 to switch between calculating the dry delay amount using only atmospheric pressure information and calculating the dry delay amount using a numerical weather model, depending on the amount of cumulonimbus clouds above the receiver i. For example, if hydrostatic equilibrium is established for the amount of cumulonimbus clouds above the receiver i, the cumulonimbus cloud determination unit 15 sends an instruction to the wet delay amount calculation unit 12 to calculate the dry delay amount using only atmospheric pressure information. On the other hand, if hydrostatic equilibrium is not established for the amount of cumulonimbus clouds above the receiver i, the cumulonimbus cloud determination unit 15 sends an instruction to the wet delay amount calculation unit 12 to calculate the dry delay amount using a numerical weather model. Other functional configurations are the same as those in the first and second embodiments.

[0079] According to the meteorological data assimilation support device 1 of the third embodiment, the variance of the wet delay is calculated using the covariance matrix of the correction values ​​of the satellite clock and satellite orbit, and the wet delay and its variance are provided as wet delay data, thereby making it possible to more accurately estimate the degree of error contained in the estimated value of the wet delay. Also, by switching the calculation method of the dry delay calculated by the wet delay calculation unit 12 to a more appropriate method depending on the amount of cumulonimbus clouds above the receiver i, the estimation accuracy of the wet delay calculated by the wet delay variance calculation unit 13 can be improved.

[0080] The hardware configuration of the meteorological data assimilation support device 1 will be described with reference to Fig. 7. As shown in Fig. 7, the meteorological data assimilation support device 1 includes a temporary storage unit 111, a storage unit 112, a calculation unit 113, an input unit 114, a transmission / reception unit 115, and a display unit 116. The temporary storage unit 111, the storage unit 112, the input unit 114, the transmission / reception unit 115, and the display unit 116 are all connected to the calculation unit 113 via a BUS.

[0081] The calculation unit 113 is, for example, a CPU (Central Processing Unit). The calculation unit 113 executes the processes of the wet delay calculation unit 12, the wet delay variance calculation unit 13, the satellite clock and orbit correction information generation unit 14, and the cumulonimbus cloud determination unit 15 of the meteorological data assimilation support device 1 in accordance with the control program stored in the storage unit 112.

[0082] The temporary storage unit 111 is, for example, a RAM (Random-Access Memory). The temporary storage unit 111 loads the control program stored in the storage unit 112 and is used as a work area for the calculation unit 113.

[0083] The storage unit 112 is a non-volatile memory such as a flash memory, a hard disk, a DVD-RAM (Digital Versatile Disc - Random Access Memory), a DVD-RW (Digital Versatile Disc - Rewritable), etc. The storage unit 112 stores in advance a program for causing the calculation unit 113 to perform the processing of the meteorological data assimilation support device 1, and also supplies the data stored by this program to the calculation unit 113 in accordance with instructions from the calculation unit 113, and stores the data supplied from the calculation unit 113.

[0084] The input unit 114 is an interface device that connects input devices such as a keyboard and a pointing device to the BUS. Information input by the user is supplied to the calculation unit 113 via the input unit 114. In a configuration in which an instruction to generate wet delay amount data is input to the meteorological data assimilation support device 1, the user inputs the instruction to the input unit 114.

[0085] The transmitter / receiver 115 is a network termination device or a wireless communication device that connects to the network, and a serial interface or a LAN (Local Area Network) interface that connects to them. The transmitter / receiver 115 functions as the atmospheric pressure information acquisition unit 11, the cumulonimbus cloud determination unit 15, and the numerical meteorological data acquisition unit 16. In a configuration in which the wet delay amount calculation unit 12 and the wet delay amount variance calculation unit 13 transmit wet delay amount data to a user terminal, a meteorological data assimilation device, etc., the transmitter / receiver 115 functions as the wet delay amount calculation unit 12 and the wet delay amount variance calculation unit 13.

[0086] The display unit 116 is a display device such as a CRT (Cathode Ray Tube) or LCD (Liquid Crystal Display). For example, the display unit 116 displays an operation screen for a user to input information. In a configuration in which a user directly inputs instructions to the meteorological data assimilation support device 1, the display unit 116 displays a screen for inputting instructions. In a configuration in which the wet delay amount calculation unit 12 and the wet delay amount variance calculation unit 13 display wet delay amount data on a screen, the display unit 116 functions as the wet delay amount calculation unit 12 and the wet delay amount variance calculation unit 13.

[0087] The processing of the atmospheric pressure information acquisition unit 11, the wet delay calculation unit 12, the wet delay variance calculation unit 13, the satellite clock and orbit correction information generation unit 14, the cumulonimbus cloud determination unit 15, and the numerical meteorological data acquisition unit 16 of the meteorological data assimilation support device 1 shown in Figures 1, 3, 4, 5, and 6 is executed by a control program using the temporary memory unit 111, the calculation unit 113, the memory unit 112, the input unit 114, the transmission / reception unit 115, the display unit 116, etc. as resources.

[0088] Furthermore, the above hardware configuration and flowchart are merely examples and can be changed and modified as desired.

[0089] The core processing components of the meteorological data assimilation support device 1, such as the calculation unit 113, temporary storage unit 111, storage unit 112, input unit 114, transmission / reception unit 115, and display unit 116, can be realized using an ordinary computer system rather than a dedicated system. For example, a computer program for executing the above operations may be stored and distributed on a computer-readable recording medium such as a flexible disk, CD-ROM (Compact Disc - Read Only Memory), or DVD-ROM (Digital Versatile Disc - Read Only Memory), and the meteorological data assimilation support device 1 that executes the above processing may be configured by installing the computer program on a computer. Alternatively, the meteorological data assimilation support device 1 may be configured by storing the computer program in a storage device of a server device on a communication network, such as the Internet, and downloading it to an ordinary computer system.

[0090] In addition, when the functions of the meteorological data assimilation support device 1 are realized by sharing the responsibilities between an OS (Operating System) and an application program, or by collaboration between the OS and an application program, only the application program portion may be stored on a recording medium or storage device.

[0091] It is also possible to provide a computer program via a communications network. For example, the computer program may be posted on a bulletin board system (BBS) on the communications network and provided via the communications network. The computer program may then be started and executed under the control of the OS in the same way as any other application program, thereby enabling the above-described processing to be performed.

[0092] In the above embodiments, the meteorological data assimilation support systems 100, 200, 200', 300 and 400 are equipped with the atmospheric delay amount calculation device 2 and the meteorological data assimilation support device 1 separately, but the meteorological data assimilation support device 1 may also be configured to have the functions of the atmospheric delay amount calculation device 2.

[0093] A modified example of the meteorological data assimilation support device 1 having the functions of the atmospheric delay calculation device 2 according to the first embodiment is shown in Fig. 8. In the meteorological data assimilation support system 100' shown in Fig. 8, the meteorological data assimilation support device 1 includes an observation quantity acquisition unit 21, a correction information acquisition unit 22, and an atmospheric delay calculation unit 23 in addition to a pressure information acquisition unit 11, a wet delay calculation unit 12, and a wet delay variance calculation unit 13. Although not shown, the meteorological data assimilation support device 1 may also include a receiver i.

[0094] A modified example of the meteorological data assimilation support device 1 having the functions of the atmospheric delay calculation device 2 according to the second embodiment is shown in Fig. 9. In the meteorological data assimilation support system 200'' shown in Fig. 9, the meteorological data assimilation support device 1 includes an observation quantity acquisition unit 21 and an atmospheric delay calculation unit 23 in addition to a pressure information acquisition unit 11, a wet delay calculation unit 12, a wet delay variance calculation unit 13, and a satellite clock and orbit correction information generation unit 14. Although not shown, the meteorological data assimilation support device 1 may also include a receiver i.

[0095] A modified example of the meteorological data assimilation support device 1 having the functions of the atmospheric delay calculation device 2 according to the third embodiment is shown in Fig. 10. In the meteorological data assimilation support system 300' shown in Fig. 10, the meteorological data assimilation support device 1 includes an observation quantity acquisition unit 21, an atmospheric delay calculation unit 23, and a correction information acquisition unit 22 in addition to a pressure information acquisition unit 11, a wet delay calculation unit 12, a wet delay variance calculation unit 13, and a cumulonimbus cloud determination unit 15. Although not shown, the meteorological data assimilation support device 1 may also include a receiver i.

[0096] Fig. 11 shows a modified example of the meteorological data assimilation support device 1 having the functions of the atmospheric delay calculation device 2 according to the fourth embodiment. In the meteorological data assimilation support system 400' shown in Fig. 11, the meteorological data assimilation support device 1 includes an observation quantity acquisition unit 21, a correction information acquisition unit 22, and an atmospheric delay calculation unit 23 in addition to a pressure information acquisition unit 11, a wet delay calculation unit 12, a wet delay variance calculation unit 13, a cumulonimbus cloud determination unit 15, and a numerical meteorological data acquisition unit 16. Although not shown, the meteorological data assimilation support device 1 may also include a receiver i.

[0097] In the above description, the second, third and fourth embodiments have been described separately, but any or all of the second, third and fourth embodiments may be combined.

[0098] In the above embodiment, correction information for the satellite clock and the satellite orbit is used, but if performance is not an issue, correction information for either the satellite clock or the satellite orbit may be used.

[0099] In the above embodiment, the meteorological data assimilation support device 1 is equipped with a barometric pressure information acquisition unit 11, but for example, if the wet delay amount calculation unit 12 and the wet delay amount variance calculation unit 13 adopt a method that uses a numerical meteorological model without using barometric pressure information to calculate the wet delay amount and the variance of the wet delay amount, respectively, the meteorological data assimilation support device 1 does not need to be equipped with the barometric pressure information acquisition unit 11.

[0100] It should be noted that the present disclosure allows various embodiments and modifications without departing from the broad spirit and scope of the present disclosure. Furthermore, the above-described embodiments are intended to explain the disclosure and do not limit the scope of the present disclosure. That is, the scope of the present disclosure is defined by the claims, not the embodiments. Various modifications made within the scope of the claims and within the meaning of the disclosure equivalent thereto are considered to be within the scope of the present disclosure. Various aspects of the present disclosure are summarized below as appendices. (Appendix 1) a wet delay amount calculation unit that calculates a wet delay amount in the satellite line of sight based on observation amount data indicating an observation amount of radio waves from the positioning satellite measured by a receiver that receives radio waves from the positioning satellite, and atmospheric delay amount data indicating an atmospheric delay amount in the satellite line of sight calculated based on at least one of a satellite clock and a satellite orbit correction value of the positioning satellite, and outputs data indicating the wet delay amount; a wet delay amount variance calculation unit that calculates the variance of the wet delay amount based on the wet delay amount and a covariance matrix of at least one of a correction value of a satellite clock and a satellite orbit of the positioning satellite, and outputs data indicating the variance of the wet delay amount; A meteorological data assimilation support device comprising: (Appendix 2) a satellite clock and orbit correction information generating unit that generates correction information including a correction value and a covariance matrix of at least one of the satellite clock and the satellite orbit, 2. The meteorological data assimilation support device according to claim 1. (Appendix 3) a pressure information acquisition unit that acquires pressure information indicating the pressure in the receiver; the wet delay amount calculation unit calculates a dry delay amount in the satellite line of sight direction based on the atmospheric pressure information, and calculates the wet delay amount by subtracting the dry delay amount from the atmospheric delay amount indicated by the atmospheric delay amount data. 3. The meteorological data assimilation support device according to claim 1 or 2. (Appendix 4) the atmospheric pressure information acquisition unit acquires the atmospheric pressure information based on latitude and date data of the receiver; 4. The meteorological data assimilation support device according to claim 3. (Appendix 5) The atmospheric pressure information acquisition unit acquires the atmospheric pressure information from an atmospheric pressure sensor provided in the receiver. 4. The meteorological data assimilation support device according to claim 3. (Appendix 6) A cumulonimbus cloud determination unit is further provided to determine the amount of cumulonimbus clouds above the receiver, The cumulonimbus cloud determination unit calculates an error variance of the dry delay amount that increases as the amount of cumulonimbus clouds above the receiver increases, the wet delay amount variance calculation unit calculates the variance of the wet delay amount based on the atmospheric pressure information, the wet delay amount, a covariance matrix of a correction value of at least one of a satellite clock and a satellite orbit of the positioning satellite, and an error variance of the dry delay amount; 6. The meteorological data assimilation support device according to any one of appendices 3 to 5. (Appendix 7) a numerical meteorological data acquisition unit that acquires numerical meteorological data related to temperature and atmospheric pressure; a cumulonimbus cloud determination unit that determines the amount of cumulonimbus clouds above the receiver; Furthermore, The cumulonimbus cloud determination unit sends an instruction to the wet delay amount calculation unit to switch between calculating the dry delay amount by a method using only the atmospheric pressure information and calculating the dry delay amount by a method using a numerical weather model based on the numerical weather data, depending on the amount of cumulonimbus clouds above the receiver; The wet delay amount calculation unit switches between calculating the dry delay amount by a method using only the atmospheric pressure information and calculating the dry delay amount by a method using the numerical meteorological model according to an instruction received from the cumulonimbus cloud determination unit. 7. The meteorological data assimilation support device according to any one of appendices 3 to 6. (Appendix 8) The meteorological data assimilation support system executes calculating a satellite line-of-sight wet delay based on observation data indicating an observation amount of radio waves from the positioning satellite measured by a receiver that receives radio waves from the positioning satellite, and atmospheric delay data indicating an atmospheric delay in the satellite line-of-sight calculated based on at least one of a satellite clock and a satellite orbit correction value of the positioning satellite; outputting data indicative of the amount of wet retardation; calculating a variance of the wet delay amount based on a covariance matrix of the wet delay amount and a correction value of at least one of a satellite clock and a satellite orbit of the positioning satellite; outputting data indicative of the variance of the wet retardation; A meteorological data assimilation support method comprising: (Appendix 9) Computer, a wet delay amount calculation unit that calculates a wet delay amount in the satellite line of sight based on observation amount data that indicates an observation amount of radio waves from the positioning satellite measured by a receiver that receives radio waves from the positioning satellite, and atmospheric delay amount data that indicates an atmospheric delay amount in the satellite line of sight calculated based on at least one of a satellite clock and a satellite orbit correction value of the positioning satellite; and a wet delay amount variance calculation unit that calculates the variance of the wet delay amount based on a covariance matrix of the wet delay amount and at least one of a correction value of a satellite clock and a satellite orbit of the positioning satellite; A program that functions as a (Appendix 10) a receiver for receiving radio waves from a positioning satellite; an atmospheric delay calculation device that calculates an atmospheric delay in the satellite line of sight based on observation data indicating the observation amount of the radio waves from the positioning satellite measured by the receiver, and the observation data and at least one of a satellite clock correction value and a satellite orbit correction value of the positioning satellite; a meteorological data assimilation support device that calculates a satellite line-of-sight wet delay based on the satellite line-of-sight atmospheric delay calculated by the atmospheric delay calculation device, calculates a variance of the satellite line-of-sight wet delay based on the satellite line-of-sight wet delay and a covariance matrix of a correction value for at least one of a satellite clock and a satellite orbit of the positioning satellite, and outputs data including the satellite line-of-sight wet delay and the variance of the satellite line-of-sight wet delay; A meteorological data assimilation support system. (Appendix 11) The atmospheric delay calculation device a satellite clock and orbit correction information generating unit that generates correction information including a correction value and a covariance matrix of at least one of the satellite clock and the satellite orbit; 11. The meteorological data assimilation support system according to claim 10. [Explanation of symbols]

[0101] 1 Meteorological data assimilation support device, 2 Atmospheric delay calculation device, 11 Pressure information acquisition unit, 12 Wet delay calculation unit, 13 Wet delay variance calculation unit, 14 Satellite clock and orbit correction information generation unit, 15 Cumulonimbus cloud determination unit, 16 Numerical meteorological data acquisition unit, 21 Observation amount acquisition unit, 22 Correction information acquisition unit, 23 Atmospheric delay calculation unit, 100,100',200,200',200'',300,300',400,400' Meteorological data assimilation support system, 111 Temporary memory unit, 112 Memory unit, 113 Calculation unit, 114 Input unit, 115 Transmitting / receiving unit, 116 Display unit, i Receiver, l Satellite.

Claims

1. a wet delay amount calculation unit that calculates a wet delay amount in the satellite line of sight based on observation amount data indicating an observation amount of radio waves from the positioning satellite measured by a receiver that receives radio waves from the positioning satellite, and atmospheric delay amount data indicating an atmospheric delay amount in the satellite line of sight calculated based on at least one of a satellite clock and a satellite orbit correction value of the positioning satellite, and outputs data indicating the wet delay amount; a wet delay amount variance calculation unit that calculates the variance of the wet delay amount based on the wet delay amount and a covariance matrix of at least one of a correction value of a satellite clock and a satellite orbit of the positioning satellite, and outputs data indicating the variance of the wet delay amount; A meteorological data assimilation support device comprising:

2. a satellite clock and orbit correction information generating unit that generates correction information including a correction value and a covariance matrix of at least one of the satellite clock and the satellite orbit; The meteorological data assimilation support device according to claim 1 .

3. a pressure information acquisition unit that acquires pressure information indicating the pressure in the receiver; the wet delay amount calculation unit calculates a dry delay amount in the satellite line of sight direction based on the atmospheric pressure information, and calculates the wet delay amount by subtracting the dry delay amount from the atmospheric delay amount indicated by the atmospheric delay amount data.

3. The meteorological data assimilation support device according to claim 1.

4. the atmospheric pressure information acquisition unit acquires the atmospheric pressure information based on latitude and date data of the receiver; The meteorological data assimilation support device according to claim 3.

5. The atmospheric pressure information acquisition unit acquires the atmospheric pressure information from an atmospheric pressure sensor provided in the receiver. The meteorological data assimilation support device according to claim 3.

6. A cumulonimbus cloud determination unit is further provided to determine the amount of cumulonimbus clouds above the receiver, The cumulonimbus cloud determination unit calculates an error variance of the dry delay amount that increases as the amount of cumulonimbus clouds above the receiver increases, the wet delay amount variance calculation unit calculates the variance of the wet delay amount based on the atmospheric pressure information, the wet delay amount, a covariance matrix of a correction value of at least one of a satellite clock and a satellite orbit of the positioning satellite, and an error variance of the dry delay amount; The meteorological data assimilation support device according to claim 3.

7. a numerical meteorological data acquisition unit that acquires numerical meteorological data related to temperature and atmospheric pressure; a cumulonimbus cloud determination unit that determines the amount of cumulonimbus clouds above the receiver; Furthermore, The cumulonimbus cloud determination unit sends an instruction to the wet delay amount calculation unit to switch between calculating the dry delay amount by a method using only the atmospheric pressure information and calculating the dry delay amount by a method using a numerical weather model based on the numerical weather data, depending on the amount of cumulonimbus clouds above the receiver; The wet delay amount calculation unit switches between calculating the dry delay amount by a method using only the atmospheric pressure information and calculating the dry delay amount by a method using the numerical meteorological model according to an instruction received from the cumulonimbus cloud determination unit. The meteorological data assimilation support device according to claim 3.

8. The meteorological data assimilation support system executes calculating a satellite line-of-sight wet delay based on observation data indicating an observation amount of radio waves from the positioning satellite measured by a receiver that receives radio waves from the positioning satellite, and atmospheric delay data indicating an atmospheric delay in the satellite line-of-sight calculated based on at least one of a satellite clock and a satellite orbit correction value of the positioning satellite; outputting data indicative of the amount of wet retardation; calculating a variance of the wet delay amount based on a covariance matrix of the wet delay amount and a correction value of at least one of a satellite clock and a satellite orbit of the positioning satellite; outputting data indicative of the variance of the wet retardation; A meteorological data assimilation support method comprising:

9. Computer, a wet delay amount calculation unit that calculates a wet delay amount in the satellite line of sight based on observation amount data that indicates an observation amount of radio waves from the positioning satellite measured by a receiver that receives radio waves from the positioning satellite, and atmospheric delay amount data that indicates an atmospheric delay amount in the satellite line of sight calculated based on at least one of a satellite clock and a satellite orbit correction value of the positioning satellite; and a wet delay amount variance calculation unit that calculates the variance of the wet delay amount based on a covariance matrix of the wet delay amount and at least one of a correction value of a satellite clock and a satellite orbit of the positioning satellite; A program that functions as a

10. a receiver for receiving radio waves from a positioning satellite; an atmospheric delay calculation device that calculates an atmospheric delay in the satellite line of sight based on observation data indicating the observation amount of the radio waves from the positioning satellite measured by the receiver, and the observation data and at least one of a satellite clock correction value and a satellite orbit correction value of the positioning satellite; a meteorological data assimilation support device that calculates a satellite line-of-sight wet delay based on the satellite line-of-sight atmospheric delay calculated by the atmospheric delay calculation device, calculates a variance of the satellite line-of-sight wet delay based on the satellite line-of-sight wet delay and a covariance matrix of a correction value for at least one of a satellite clock and a satellite orbit of the positioning satellite, and outputs data including the satellite line-of-sight wet delay and the variance of the satellite line-of-sight wet delay; A meteorological data assimilation support system.

11. The atmospheric delay calculation device a satellite clock and orbit correction information generating unit that generates correction information including a correction value and a covariance matrix of at least one of the satellite clock and the satellite orbit; The meteorological data assimilation support system according to claim 10.

Images