Login| Sign Up| Help| Contact|

Patent Searching and Data


Title:
SYSTEMS AND METHODS FOR IMPROVED ATMOSPHERIC MONITORING AND GPS POSITIONING UTILIZING GNSS TOMOGRAPHIC REFRACTIVITY
Document Type and Number:
WIPO Patent Application WO/2019/094255
Kind Code:
A1
Abstract:
The disclosed technology relates to systems and methods for determining three-dimensional atmospheric and ionospheric density using refraction of electromagnetic waves. A method is provided for receiving, at a processing system, and from a plurality of Global Navigation Satellite Systems (GNSS) stations, navigation data corresponding to computed positions of the plurality of GNSS stations. The method can further include determining, based at least in part on received navigation data and received GNSS transmitter information, ionosphere and atmosphere refractivity corresponding to intersections of two or more GNSS signals. The method can include calculating, based on the determined 3D density states, data fields of a model representing the three-3D density states. The method can include transmitting position adjustment data to calibrate a navigation position of at least one of the plurality of the GNSS stations based at least in part on the calculated data fields of the model.

Inventors:
MACDONALD ALEXANDER (US)
PLATZER PETER (US)
Application Number:
PCT/US2018/058760
Publication Date:
May 16, 2019
Filing Date:
November 01, 2018
Export Citation:
Click for automatic bibliography generation   Help
Assignee:
SPIRE GLOBAL INC (US)
International Classes:
G01S19/14; G01S1/04; G01W1/00
Foreign References:
US20140253375A12014-09-11
US20060080038A12006-04-13
Other References:
SOLHEIM ET AL.: "Propagation delays induced in GPS signals by dry air, water vapor, hydrometeors, and other particulates", JOURNAL OF GEOPHYSICAL RESEARCH, vol. 104, 1 April 1999 (1999-04-01), Retrieved from the Internet [retrieved on 20190107]
Attorney, Agent or Firm:
LI, Dan (US)
Download PDF:
Claims:
Claims

What is claimed is:

1. A method comprising:

receiving, at a processing system, and from a plurality of Global Navigation Satellite Systems (GNSS) stations, navigation data corresponding to computed positions of the plurality of GNSS stations;

receiving, at the processing system, GNSS transmitter information comprising one or more of:

GNSS satellite positional information;

GNSS clock information; and

GNSS transmission frequency information;

determining, based at least in part on the received navigation data and the received GNSS transmitter information, ionosphere and atmosphere refractivity associated with navigation data; determining, based at least in part on the refractivity, three-dimensional (3D) density states of the atmosphere and ionosphere corresponding to intersections of two or more GNSS signals;

calculating, by the processing system, and based on the determined 3D density states, data fields of a model representing the three-3D density states; and

transmitting position adjustment data to calibrate a navigation position of at least one of the plurality of the GNSS stations based at least in part on the calculated data fields of the model.

2. The method of claim 1, wherein determining the ionosphere and atmosphere refractivity comprises computing density delay components for dry air, water, and ions.

3. The method of claim 1, wherein the plurality of GNSS stations include one or more of: a CubeSat, a mobile GPS receiver, a mobile computing device, a cell phone, a receiver of an airplane, and a smartphone.

4. The method of claim 1, wherein at least a portion of the plurality of GNSS stations comprise one or more local processors, and wherein the computed positions are determined by the corresponding local processors.

5. The method of claim 1, wherein the GNSS transmitter information is received from one or more of:

a central command and control (CCC) system;

one or more GNSS stations; and

a precise positioning service.

6. The method of claim 1, wherein the navigation data and the GNSS transmitter information corresponds to two or more GNSS satellites, and wherein determining the 3D density states comprises determining intersection angles of the two or more GNSS signals.

7. The method of claim 1, wherein determining the 3D density states comprises determining dynamic changes associated with the atmosphere and ionosphere.

8. The method of claim 1, wherein determining the 3D density states comprises predicting future density states associated with the atmosphere and ionosphere.

9. The method of claim 1, wherein the data fields of the model are refined by feedback.

10. The method of claim 1, wherein determining the ionosphere and atmosphere refractivity comprises determining the 3D density states by minimizing a difference in least-squares between integrated refraction of density over one-dimensional rays to determine an optimal three- dimensional structure.

11. A communications system comprising:

a central command and control system configured for communication with at least one satellite associated with a constellation, and by at least one ground station;

at least one memory for storing data and computer-executable instructions, the at least one memory comprising a resource database configured for storing knowledge data

corresponding to a plurality of components associated with the satellite constellation

communications system, wherein the plurality of components comprises the at least one satellite associated with the constellation;

at least one processor in communication with the at least one memory, wherein the at least one processor is further configured to execute the computer-executable instructions to cause the system to: receive, at a processing system, and from a plurality of Global Navigation Satellite Systems (GNSS) stations, navigation data corresponding to computed positions of the plurality of GNSS stations;

receive, at the processing system, GNSS transmitter information comprising one or more of:

GNSS satellite positional information;

GNSS clock information; and

GNSS transmission frequency information;

determine, based at least in part on the received navigation data and the received GNSS transmitter information, ionosphere and atmosphere refractivity associated with navigation data;

determine, based at least in part on the refractivity, three-dimensional (3D) density states of the atmosphere and ionosphere corresponding to intersections of two or more GNSS signals;

calculate, by the processing system, and based on the determined 3D density states, data fields of a model representing the three-3D density states; and

transmit position adjustment data to calibrate a navigation position of at least one of the plurality of the GNSS stations based at least in part on the calculated data fields of the model.

12. The system of claim 1 1, wherein the ionosphere and atmosphere refractivity is determined by computing density delay components for dry air, water, and ions.

13. The system of claim 1 1, wherein at least a portion of the plurality of GNSS stations comprise one or more local processors, and wherein the computed positions are determined by the corresponding local processors.

14. The system of claim 1 1, wherein the GNSS transmitter information is received from one or more of:

a central command and control (CCC) system;

one or more GNSS stations; and

a precise positioning service.

15. The system of claim 11, wherein the navigation data and the GNSS transmitter information corresponds to two or more GNSS satellites, and wherein determining the 3D density states comprises determining intersection angles of the two or more GNSS signals.

16. The system of claim 11, wherein determining the 3D density states comprises determining dynamic changes associated with the atmosphere and ionosphere.

17. The system of claim 11, wherein the at least one processor is further configured to execute the computer-executable instructions to cause the system to determine at least one atmospheric state based at least in part on the calculated data fields of the model.

18. The system of claim 11, wherein the data fields of the model are refined by feedback.

19. The system of claim 11, wherein determining the ionosphere and atmosphere refractivity comprises determining the 3D density states by minimizing a difference in least-squares between integrated refraction of density over one-dimensional rays to determine an optimal three- dimensional structure.

20. The system of claim 11, wherein the plurality of GNSS stations include one or more of: a CubeSat, a mobile GPS receiver, a mobile computing device, a cell phone, a receiver of an airplane, and a smartphone.

Description:
SYSTEMS AND METHODS FOR IMPROVED ATMOSPHERIC MONITORING AND GPS POSITIONING UTILIZING GNSS TOMOGRAPHIC REFRACTIVITY

FIELD

[0001] The disclosed technology relates to systems and methods for determining three- dimensional atmospheric and ionospheric density using refraction of electromagnetic waves. Improved atmospheric and ionosphere three-dimensional density fields are used to constrain errors in GNSS navigation devices, allowing more accurate estimate of location. Additionally, iterative improvements between three-dimensional refraction fields and GNSS locations provides improved atmospheric and space weather predictions, and improved GNSS location capabilities.

BACKGROUND

[0002] Global Navigation Satellite Systems (GNSS) are satellite systems that can include the United States' Global Positioning System (GPS), Russia's GLONASS, and the European Union's GALILEO. Such systems are used primarily in navigation applications, for example, to determine receiver positioning. GPS has been in operation since 1978 (globally since 1994) and was initially developed to provide precise positioning for military purposes. Today, GNSS systems are used for many civilian and military applications including navigation, surveying, time referencing, geo-fencing, weather data, self-driving cars, etc.

[0003] As the GNSS electromagnetic signals traverse the intervening distances from satellites in orbit to various receivers, the signals can be refracted due to the negatively charged ionosphere, moisture in the atmosphere, and the mass of the dry atmosphere. Such refraction can create errors in positional accuracy of GPS trilateration. Previous GNSS stations utilize mapping functions to determine the appropriate dry atmosphere delay for signals arriving at different angles above the horizon. These mapping functions generally use an estimate, such as surface pressure, to calculate the amount of mass traversed for any angle. In addition, the error introduced by the ionosphere can be relatively large when the satellite is near the observer's horizon, when the vernal equinox is near, and/or when sunspot activity is severe. Error can also be introduced by non-uniform moisture content in the atmosphere. The corresponding positional errors can vary as a function of weather, magnetic activity, location, time of day, direction of observation, temperature, atmospheric pressure, etc. [0004] A need exists for improved systems and methods to address such challenges.

BRIEF SUMMARY

[0005] Some, or all, of the above needs may be addressed by certain embodiments and implementations disclosed herein.

[0006] Certain implementations of the disclosed technology may include a method of receiving, at a processing system, and from a plurality of Global Navigation Satellite Systems (GNSS) stations, navigation data corresponding to computed positions of the plurality of GNSS stations. The method can include receiving, at the processing system, GNSS transmitter information including one or more of: GNSS satellite positional information; GNSS clock information; and GNSS transmission frequency information. The method can further include: determining, based at least in part on the received navigation data and the received GNSS transmitter information, ionosphere and atmosphere refractivity associated with navigation data; determining, based at least in part on the refractivity, three-dimensional (3D) density states of the atmosphere and ionosphere corresponding to intersections of two or more GNSS signals; calculating, by the processing system, and based on the determined 3D density states, data fields of a model representing the three-3D density states; and transmitting position adjustment data to calibrate a navigation position of at least one of the plurality of the GNSS stations based at least in part on the calculated data fields of the model.

[0007] According to another example implementation, a satellite constellation communications system is provided, which may include a central command and control system configured for communication with at least one satellite associated with a constellation, and by at least one ground station. The system may include at least one memory for storing data and computer-executable instructions, the at least one memory including a resource database configured for storing knowledge data corresponding to a plurality of components associated with the satellite constellation communications system, wherein the plurality of components include the at least one satellite associated with the constellation. The system may include at least one processor in communication with the at least one memory, wherein the at least one processor is further configured to execute the computer-executable instructions to cause the system to: receive, at the processing system, GNSS transmitter information including one or more of: GNSS satellite positional information; GNSS clock information; and GNSS transmission frequency information. The at least one processor is further configured to: determine, based at least in part on the received navigation data and the received GNSS transmitter information, ionosphere and atmosphere refractivity associated with navigation data; determine, based at least in part on the refractivity, three-dimensional (3D) density states of the atmosphere and ionosphere corresponding to intersections of two or more GNSS signals; calculate, based on the determined 3D density states, data fields of a model representing the three-3D density states; and transmit position adjustment data to calibrate a navigation position of at least one of the plurality of the GNSS stations based at least in part on the calculated data fields of the model.

[0008] Other implementations, features, and aspects of the disclosed technology are described in detail herein and are considered a part of the claimed disclosed technology. Other implementations, features, and aspects can be understood with reference to the following detailed description, accompanying drawings, and claims.

BRIEF DESCRIPTION OF THE FIGURES

Reference will now be made to the accompanying figures and flow diagrams, which are not necessarily drawn to scale, and wherein:

[0009] FIG. 1 depicts a block diagram representation of an example GNSS system 100 according to an example implementation of the disclosed technology.

[0010] FIG. 2 depicts an example GNSS system 200 (not to scale) relative to earth 201 with associated orbital (226), airborne (218), oceanic (222), and earthbound (214, 216, 220) GNSS stations 104.

[0011] FIG. 3 A shows an initial surface wind and pressure "true" field.

[0012] FIG. 3B shows calculated surface wind and pressure field from the Observing System Simulation Experiments (OSSEs) processes after a few hundred iterations.

[0013] FIG. 4A shows an initial wind and pressure "true" field at 5km elevation.

[0014] FIG. 4B shows calculated 5km elevation wind and pressure field from the OSSEs processes after a few hundred iterations. [0015] FIG. 5 depicts a block diagram of an illustrative computing device 500 according to an example implementation.

[0016] FIG. 6 is flow diagram of a method 600, according to an example implementation of the disclosed technology.

[0017] FIG. 7 is flow diagram of another method 700, according to an example implementation of the disclosed technology

DETAILED DESCRIPTION

[0018] Some implementations of the disclosed technology will be described more fully hereinafter with reference to the accompanying drawings. This disclosed technology may however, be embodied in many different forms and should not be construed as limited to the implementations set forth herein. Some of the components illustrated in the accompanying figures are shown for illustration purposes only, and may or may not be drawn to scale. In the following detailed description, numerous specific details are set forth by way of examples to provide a thorough understanding of the relevant teachings. However, it should be apparent to those skilled in the art that the present teachings may be practiced without such details. In other instances, well known methods, procedures, components, and/or circuitry have been described at a relatively high-level, without detail, to avoid unnecessarily obscuring aspects of the present teachings.

[0019] The disclosed technology generally relates to systems and methods for determining static and/or dynamic three-dimensional density states of ionospheric electrons and atmospheric mass (hereinafter, referred to "3D density states") using received GNSS electromagnetic signals that traverse through these mediums. Certain implementations extend and improve earlier work entitled "Diagnosis of Three-Dimensional Water Vapor Using a GPS Network," by MacDonald et al., Monthly Weather Review, Vol. 130, pp. 386-397, 2002, incorporated herein by reference as presented in full.

[0020] Prior GNSS systems utilize mapping functions to approximate an appropriate delay for GNSS signals arriving at GNSS stations from different angles above the horizon. Such mapping functions generally use an estimate, such as surface pressure, to calculate the amount of mass traversed for a given angle, and the density is assumed to be horizontally uniform at all levels around the receiver, which is a false simplifying assumption leading to a previously acceptable degree of error in device location.

[0021] Certain example implementations of the disclosed technology measure the density delay of the GNSS signals as caused by the ionosphere and/or atmosphere, recognizing that the mass is not horizontally uniform. Certain example implementations utilize the density delay to determine the non-uniformity, leading to a better estimate of the 3D density state. [0022] Certain implementations of the disclosed technology may utilize new measurement, analysis, and/or feedback techniques to improve the accuracy/resolution of the determined 3D density states. Certain implementation may be further configured to provide additional information, such as momentum in three directions, temperature, and pressure. For example, a significant increase in the number of (and improved utilization of) measurements may provide enhanced information about the 3D density states, which can provide improvements in weather prediction. The disclosed technology can apply to very short-range time-scale events (such as tornados and flash floods), medium-range time-scale events (such as snowstorms and hurricanes), and/or long-range time-scale events (such as droughts and heat waves).

[0023] In certain example implementations, the improved information about the 3D density states may be utilized to correct for refractive delays to reduce errors in positional accuracy of navigational systems, as will be further discussed below.

[0024] Certain example implementations of the disclosed technology may utilize two-way feedback to improve knowledge of the 3D density states. In certain example implementations, information identifying a transmitter's location, together with the determined refraction-related delays of the received signals may reduce inaccuracies of the 3D density states. Therefore, as the 3D density states are better known, such information can create better estimates of refractive effects for receivers in a large, proximate geographic region, which then may constrain the location error and further improve the estimate of the 3D density states. This continual iteration between 3D density states and location information may enable very high-resolution weather prediction models to extend improved state knowledge into the future, resulting in better estimates for future geo-location determination. According to an example implementation of the disclosed technology, continual feedback and iteration may be utilized to improve weather prediction and geo-location.

[0025] FIG. 1 depicts a block diagram representation of an example GNSS system 100 according to an example implementation of the disclosed technology. The system 100 can include a plurality of GNSS transmitters 102 deployed into orbit. For example, the GNSS transmitters 102 may be installed on GNSS satellite constellations, such as United States' Global Positioning System (GPS), Russia's GLONASS, the European Union's GALILEO systems, etc. In certain example implementations, the system 100 includes a plurality of GNSS stations 104 configured with GNSS stations. The GNSS stations 104, for example, may be mounted on ships, fixed terrestrial monuments, airplanes, vehicles, CubeSats, mobile GPS receivers, mobile computing devices, cell phones, smartphones, etc. GNSS signals originating from the transmitters 102 may traverse the ionosphere and atmosphere and may be detected and processed by the receivers 104. In certain example implementations, at least a portion of the receivers 104 may be configured for two-way communication with a data exchange system 106.

[0026] Since the effective refractive indices of the ionosphere and the atmosphere are nonuniform, and dynamically changing as a function of local electron density and water vapor content, the various GNSS signals undergo different refraction-related delays, depending on the relative positions of the transmitters, receivers, and intervening media through which the signals traverse.

[0027] In accordance with certain example implementations of the disclosed technology, the various calculations may be done in the data exchange system 106, the command and control center (CCC) 107, and/or the processing system(s) 108. For example, the calculation of refraction-related delay may utilize highly accurate position and clock information 111 for the GNSS satellite constellation 102. Such position and clock information 111, for example, may be available from services provided by Precision Point Position companies. In certain example implementations, the CCC 107, and/or the processing system(s) 108 may also compute updated analysis of the atmosphere and ionosphere that is communicated back to receivers in the GNSS stations 104 to improve their location.

[0028] In certain example implementations, the data exchange system 106 may be in communication with a command and control center (CCC) system 107, for example, to receive information about the GNSS transmitters 102, such as position, clock information, broadcast frequencies, etc. According to an example implementation of the disclosed technology, the data exchange system 106 may be in communication with one or more processing systems 108 configured with special-purpose software that instructs one or more processor(s) of the system 108 to analyze and/or assimilate the data received by the CCC system 107 and the GNSS stations 104 via the data exchange system 106. [0029] In certain example implementations, the data exchange system 106 may receive information about the GNSS transmitters 102 (such as position, clock information, and/or broadcast frequencies) from one or more of the GNSS stations 104, as appropriate.

[0030] In certain example implementations, the analysis/assimilation carried out by the processing system 108 may be used to produce one or more models 110 representing the 3D density states of the ionosphere and/or atmosphere. According to an example implementation of the disclosed technology, the determined 3D density states may be continuously updated (via the model 110) as new information becomes available.

[0031] In accordance with an example implementation of the disclosed technology, the data exchange system 106 and/or the processing system 108 may be utilized to provide feedback from the model 110 to at least a portion of the receivers 104, for example, for calibration, derived positional corrections, refined local weather events, predictions, etc. In certain example implementations, the system 100 may be utilized to continuously refine the information represented by the model 110. In certain example implementations, the feedback sent back to the receivers 104 may be processed/utilized by the receivers 104 to reduce positional errors, for example, due to GNSS signal delay errors caused by the corresponding variations in the ionosphere and atmosphere.

[0032] As indicated by the double arrows between the receivers 104, the data exchange system 106, the CCC system 107, and the processing system(s) 108 data may be passed in both directions. For example, delay and/or location computed at the GNSS stations 104 may be transmitted via the data exchange system 106 to the processing system 108, which may be used in conjunction with GNSS transmitter 102 positional, clock, and/or frequency information to update the model(s) 110. In certain example implementations, three and four-dimensional atmospheric states (including predictions) may be fed back from processing system 108 to the GNSS stations 104 to improve their location accuracy.

[0033] The continual iterative process as disclosed herein may be utilized to improve the weather and ionospheric models 110 by using "density delay distances" (explained below) to better estimate the initial states (i.e. three-dimensional densities).

[0034] The densities, together with suitable assumptions (such as a hydrostatic atmosphere) enables derivation of the complete dynamic state of the atmosphere, including winds, temperature and pressure. Accordingly, the processing system 108 may use these observations, and other data sources such as GNSS radio occultation soundings and other global weather data, to refine the atmospheric state. The initial state of the atmosphere and ionosphere can be used together with a neutral atmosphere and space weather prediction model to project data forward a short time (e.g. 5 to 60 minutes), generating a new, slightly different state. In accordance with an example implementation of the disclosed technology, the estimate of future state may be fed back to the GNSS stations 104, allowing more accurate constraints on errors caused by dry atmosphere and volume electron content, resulting in better location estimates. Such a system could bring average errors to less than a centimeter, which makes the data being fed back to the central processor more accurate.

[0035] FIG. 2 depicts an example GNSS system 200 (not to scale) relative to earth 201 with associated orbital (226), airborne (218), oceanic (222), and earthbound (214, 216, 220) GNSS stations 104. The GNSS system 200 may include all or a portion of the GNSS transmitters 102, GNSS stations 104, data exchange system 106, CCC 107, processing system 108, and model 110 as discussed above with respect to the GNSS system 100 of FIG. 1. The associated GNSS signals 230, 232, 234 from the corresponding GNSS transmitters 202, 204, and 206 are illustrated in FIG. 2 as having localized refraction "bends" in the rays as higher refractive index media in the atmosphere is encountered, however, this is simplified for illustration purposes only, as the actual refraction may be less pronounced, and/or distributed over longer distances.

[0036] In accordance with an example implementation of the disclosed technology, the data exchange system 106 may include and/or may be in communication with a network 228, such as the internet, etc. According to certain example implementations, the network 228 may include and/or be in communication with one or more cell towers 225 configured for communication with the mobile GNSS stations 214, terrestrial vehicular GNSS stations 220, etc. In certain example implementations, the fixed GNSS stations 216 may be in communication with the network 228.

[0037] In an example implementation, the transmitters 102 may be associated with constellation satellites, such as United States' Global Positioning System (GPS), Russia's GLONASS, the European Union's GALILEO systems, etc. In an example implementation, a first GNSS satellite transmitter 202 may transmit a first GNSS signal 230, which may traverse through the ionosphere 208 (having non-uniform electron density) and the atmosphere 210 (having non-uniform dry mass 212), each of which can vary by time and relative positioning with respect to the signal 230. Accordingly, the first transmitted GNSS signal 230 may be received by one or more GNSS stations 104, which can include mobile GNSS stations 214, fixed GNSS stations 216, airborne GNSS stations 218, terrestrial vehicular GNSS stations 220, shipping GNSS stations 222 at sea 224, CubeSat radio occultation receivers 226, etc.

[0038] With continued reference to FIG. 2, and in an example implementation of the disclosed technology, additional GNSS satellite transmitters (such as a second GNSS satellite transmitter 204, and/or a third GNSS satellite transmitter 206, etc.) may also transmit their corresponding GNSS signals (such as a second GNSS signal 232 from the second GNSS transmitter 204, a third GNSS signal 234 from the third GNSS transmitter 206, etc.), which may be received by one or more of the GNSS stations 104.

[0039] In accordance with an example implementation of the disclosed technology, the system 200 may include a CubeSat GNSS radio occultation receiver 226 configured for use with s GPS Radio Occultation (GNSS-RO) system. GNSS-RO is a technique for measuring properties of the Earth's atmosphere from space, and it can be used for weather forecasts, climate modeling, and space weather prediction, as described in U.S. Patent Publication No. US2015/0192696, the contents or which are incorporated herein by reference as if presented in full.

[0040] A GNSS-RO receiver 226, for example, can receive a radio signal 236 transmitted by a GNSS satellite (such as depicted by the second GPS satellite transmitter 204) in medium Earth Orbit. Depending on the density of the ionosphere and atmosphere through which the transmitted signal 236 travels, the transmitted signal 236 can be "bent" towards the earth 201 and slowed down to travel a longer path both spatially and temporally through the ionosphere/atmosphere, and the bending angle can be calculated based on the delay in the signal's arrival. The bending angle and density-induced delay can be used to calculate the ionospheric/atmospheric density along the signal's path, which in turn can be used to calculate temperature, pressure, humidity and electron density.

[0041] With continued reference to FIG. 2, and in an example implementation of the disclosed technology, the combination of GNSS received rays 230, 232, 234 crossing in different somewhat orthogonal directions, allows tomographic recovery of the 3D density states, as discussed above. In certain example implementations, the received, refractivity signals are assimilated by the processing system 108 to determine the atmospheric mass and ionospheric electron density.

[0042] Traditionally, tomography has been utilized as a technique to determine the three- dimensional structure of an object by measuring attenuation of multiple penetrating waves. For example, in three-dimensional medical tomography, an electromagnetic wave (typically an x- ray) is transmitted through a biological object, and total signal attenuation along each ray is measured. When the media is probed from many different angles, a two or three-dimensional portrayal of the interior structure (e.g. internal bones and organs) can be mathematically derived.

[0043] Certain example implementations of the disclosed technology utilize a new form of tomography in which the refraction of the signal is utilized (rather than attenuation) to determine the internal structure of the medium.

[0044] Refraction is defined as a change in direction of wave propagation due to a change in its associated transmission media. For example, as an electromagnetic wave propagates through media having different refractive indices, the phase velocity of the wave is changed but its frequency remains constant. The refractive index of a media is defined as:

[0045] n =—,

V

[0046] where c is the speed of light, and v is the velocity of the electromagnetic wave through the media. Thus, in a vacuum, the refractive index is n = 1 since v = c. The index of refraction for air at pressure of one atmosphere is 1.000293, with an associated electromagnetic wave velocity through air expressed as:

[0047] v air

i r = — 1.00 -0—293 299,700— s ,

[0048] or about 90 km/s slower than the speed of light in a vacuum. When electromagnetic waves traverse from a first media having a first refractive index ni and are incident at a first angle θι with respect to normal incidence at an interface of a second media of refractive index τΐ2, the emerging wave is bent travels through the second media at a second angle Θ2 with respect to normal according to Snell's Law:

[0050] It is common to define another index, N, referred to as refractivity and defined as:

[0051] N = 10 6 (n - 1).

[0052] The refractivity N can be expressed as a sum of the various contributing components, including the dry atmosphere density, water vapor density, liquid water content, and ionosphere electron density. For GNSS satellite electromagnetic waves, the refractivity may be expressed according to the expressions:

[0053] N = dry + wet + liquid + ion,

[0054] with the corresponding terms expressed as:

[0055] N = 222.76 p a + 1.72xl0 6 (^) + lAxW ~3 W L + 4.03xl0 7 ^, [0056] where:

[0057] p a = atmospheric density in kg/m 3 ; [0058] Pv= water vapor density in kg/m 3 ; [0059] T = atmospheric temperature in Kelvin [0060] de = electron density (ions per m 3 ); [0061] f = transmitter frequency in Hz; and [0062] WL = liquid water droplet density in kg/m 3 .

[0063] In accordance with an example implementation of the disclosed technology, the refraction of GNSS signals caused by the "dry" and "ion" terms may be utilized to recover the 3D density states for those substances. In an example implementation, the "liquid" term due to water droplets are negligible at electromagnetic wavelengths used by GNSS and can be eliminated from the calculation. According to an example implementation of the disclosed technology, the "wet" (water vapor) term can be subtracted from the integrated refractivity N to estimate the contribution of the other two remaining terms: "dry" and "ion," for which the "ion" term may be separated using two different L-band signals, as will be explained below.

[0064] Accordingly, and in an example implementation of the disclosed technology, the effects of intervening media on the GNSS electromagnetic waves to determine the integrated refractivity N may be characterized by accurately measuring the distance that a media-modified wave is reduced in a period. For example, as the wave velocity is slowed down and/or bent by the ionosphere and atmosphere, it does not reach the same point that it would reach in a vacuum in a given amount of time. The delayed distance difference between where the wave would be in a vacuum, versus where it actually is, can be inferred accurately by utilizing clocks in the receivers. This delayed distance is herein referred to as the "density delay distance," abbreviated as D3. In certain example implementations, the D3 may be represented as a time delay. As defined herein, the term "density delay information" may include time delay and/or distance delay due to refractive effects of the atmosphere and ionosphere.

[0065] In accordance with an example implementation of the disclosed technology, the D3 may be calculated in real-time using known locations and clock corrections for the GNSS satellites. In another example implementation, information about the locations and clock corrections for the GNSS satellites may be obtained, for example, by via a command and control center (CCC) 107 (or via a service that provides satellite position and clock information 111 as previously described) and such information may be communicated to and processed by the processing system 108 for more precise satellite location and/or clock correction information.

[0066] The GPS system utilizes at least 24 satellites in a constellation having an orbital altitude of 20, 180 km. Each satellite in the GPS constellation broadcasts signals over the LI (1.57542 GHZ), L2 (1.2276 GHZ), and L5 (1.17645 GHz) frequencies. The transmitted signals are encoded with the corresponding satellite's orbital position and time of transmission, with synchronization of the constellation maintained by atomic clocks. A receiver compares the time of broadcast encoded in the transmission of three or more different satellites, allowing the receiver's position to be calculated using trilateration based on the corresponding time-of-flight delays of each received signal. In accordance with certain example implementations, the system (100 and or 200) may augment the GNSS satellite position and clock information 111 obtained from a precise positioning service with the three-dimensional density fields of atmospheric mass, moisture and electron density to achieve a more accurate position. In certain example implementations, this calculation may be done at the central processor 108, and may or may not be transmitted back to any particular GNSS station having a receiver.

[0067] A consequence of the dispersive nature of the ionosphere is that the apparent time delay for a higher frequency carrier wave is less than it is for a lower frequency wave. That means that LI (1.57542 GHZ) is not affected as much as L2 (1.2276 GHZ), and L2 is not affected as much as L5 (1.17645 GHz). In accordance with an example implementation of the disclosed technology, the separation between the LI, L2, and/or L5 frequencies may be utilized to facilitate an estimation of the ionospheric group delay to determine the "ion" term, as discussed above. In certain example implementations, multiple-frequency receivers may track the carrier frequencies to determine the total electron content, which has a frequency dependence given by:

At*c*f 2

[0068] TEC =

40.3

[0069] and where:

[0070] TEC = the total electron content per square meter; [0071] /= transmitter frequency in Hz; [0072] At = the time delay; and [0073] c = the speed of light in m/s.

[0074] In accordance with an example implementation of the disclosed technology, the availability of TEC along a large number of intersecting paths allow a tomography calculation of the three-dimensional electron density of the ionosphere. The availability of the three- dimensional electron content of the ionosphere allow GNSS station receivers that do not have the two frequencies needed to estimate TEC to use the ray path through the ionosphere to estimate the D3 due to TEC. Thus, certain example implementations of the disclosed technology may enable inexpensive receivers in GNSS stations to have significantly higher accuracy than would be otherwise possible. [0075] In one example implementation, the TEC may be calculated via a processor associated with the local receiver 104. In this implementation, the precise satellite location and clock corrections may be obtained via the data exchange system 106. In certain example implementations, the local receiver 104 may transmit data to a processing system 108 via the data exchange system 106, where corrections are available to do the D3 calculations. In certain example implementations, data flow may be optimized for the particular communication path, associated processors, etc.

[0076] As previously indicated, prior GNSS systems may utilize so-called "mapping functions" to approximate delays for GNSS signals arriving at GNSS stations from different angles above the horizon. Such mapping functions make a simplifying assumption that the density is horizontally uniform at all levels around the receiver, which can lead to errors in the computed location of the device.

[0077] Consider a scenario where GNSS signals are detected by a surface receiver. For example, a signal transmitted by a GNSS satellite at the zenith (directly above the receiver) could have a D3 of approximately 2.1 meters. A signal transmitted by a GNSS satellite at 45 degrees (with respect to zenith) from the west would have a D3 of about 2.97 meters. If the air temperature were uniform, a signal transmitted by a GNSS satellite at -45 degrees (with respect to zenith) from the east would have a similar D3 of about 2.97 meters. However, if the air to the east of the receiver is much colder, and denser, it could make the corresponding D3 about 3.02 meters, or 5 cm longer than the D3 measured from the west due to the integrated density difference.

[0078] Certain example implementations of the disclosed technology may utilize very accurate receivers that can measure the time within the carrier phase signal to determine D3 to within about a centimeter of root mean square accuracy. For GNSS stations that have the atmosphere and ionosphere as intervening media, if the angle between the in-situ receiver and the GNSS satellite is greater and 15 degrees above the horizon, the bending effect can be ignored. For signals coming in a low angle, the bending of the wave creates a longer path length, which are included in the D3 calculation.

[0079] Certain example implementations of the disclosed technology utilize iterative ray tracing of GNSS signals to determine D3 with extremely high accuracy. Certain example implementations may improve navigational positioning accuracy by measuring the intervening media and/or making higher accuracy estimates of the associated density.

[0080] Experiments confirm that GNSS signals may be utilized to determine the three- dimensional density of dry air in the atmosphere. Furthermore, water vapor in the lower atmosphere can be determined by using minimization techniques. Certain example implementations of the disclosed technology may calculate a density delay due to electrons via differential delay between the LI and L2 band frequencies used by GNSS constellations, as discussed above. Certain example implementations of the disclosed technology derive all three (dry, wet, ion) fields using variational assimilation techniques. In certain example implementations of the disclosed technology, the atmospheric dynamic fields may be calculated by removing the D3 effect of electron density before the assimilation. Briefly described, the assimilation system arrives as a three-dimensional field by minimizing the difference between the D3 observations and the iterated three-dimensional fields using numerical variational techniques. Such minimization, variational assimilation, and numerical techniques are discussed in "Atmospheric Modeling, Data Assimilation and Predictability," Kalney 2003, Cambridge University Press, the contents of which are incorporated herein by reference, as if presented in full.

[0081] Certain example implementations of the disclosed technology utilize large numbers of receivers 104. In one example implementation of the disclosed technology, 5-10 receivers 104 are utilized to improve the accuracy of the 3D density states. In another example implementation, 11-50 receivers 104 are utilized. In another example implementation, 51-100 receivers 104 are utilized. In yet another example implementation, 101-1000 receivers 104 are utilized.

[0082] As an illustration of large numbers of observations, consider a 1000 km by 1000 km domain over populated advanced countries like the US or Europe. Such a domain could have about 200 Radio Occultation observations, a couple of hundred aircraft observations and a couple of thousand surface based observations. The scale of hourly weather that affects GNSS and prediction can be estimated as about 100 km by 100 km. Such a domain would have a couple hundred observations, which would allow horizontal and vertical characterization of the weather phenomena of importance. According to an example implementation of the disclosed technology, the more receivers 104 that are utilized in the determination of the 3D density states, the more accurate the model becomes.

[0083] As previously discussed, certain example implementations of the disclosed technology may share certain similarities with medical tomography, in which an x-ray generation device rotates around a patient such that the x-rays traverse through the patient's tissue in many different directions, and from which a tomogram may be reconstructed from the detected transmissions.

[0084] In accordance with certain example implementations of the disclosed technology, different receivers may be utilized to receive the different signals arriving from different directions by different satellites to perform the tomography. For example, in the technique of Radio Occultation, a GNSS satellite signal 236 may be transmitted by a GNSS satellite 204 and received by another satellite 226 after a dominantly horizontal traverse through the atmosphere. In certain example implementations of the disclosed technology, surface receivers such as ships 222, automobiles 220, fixed ground based devices 216, mobile devices 214, unmanned aerial vehicles and other aircraft 218, etc., may receive predominantly vertical signals (such as may be received from an overhead GNSS satellite 202). Thus, the horizontal Radio Occultation signals and the vertical in-situ signals may be approximately orthogonal, and may provide the necessary delay information for constructing the 3D density states. Furthermore, location information may be transmitted from at least a portion of the receivers 104 to the processing system 108 for additional data to further reduce errors. The relative position of the different receivers 104 with respect to the GNSS transmitters 102 may create multiple crossings from different directions for horizontal domains within each level of the atmosphere.

[0085] In accordance with certain example implementations of the disclosed technology, GNSS signals from two or more GNSS transmitters 102 may cross or intersect at an angle a 238 and the intersection may define a location, volume, and/or region in the intervening space (i.e., in the ionosphere 208 and/or atmosphere 210) between the associated GNSS transmitters 102 and GNSS stations 104. According to an example implementation of the disclosed technology, the density state of the location, volume, and/or region (at and/or surrounding) the intersection may contribute to the same or similar density delay for the associated GNSS signals penetrating the same location, volume, and/or region. While an orthogonal intersection angle a 238 (i.e., 90 degrees) may prove ideal for precisely determining locations of ionospheric/atmospheric content associated with the density delay contributions, certain example implementations of the disclosed technology may utilize weightings to account for intersecting angles a 238 that are less than 90 degrees. In one example implementation of the disclosed technology, a trigonometric function, such as sin(a) may be utilized as weighting so that certain non-orthogonal crossing signals may still be utilized by the processing system 108, and incorporated into the model 110.

[0086] In accordance with certain example implementations of the disclosed technology, and as determined by developmental testing and experimentation, tomography using the techniques discussed herein may be sufficiently accurate when approximately 20 to 40 percent (or greater) of the rays cross at angles greater than 30 degrees.

[0087] In accordance with certain example implementations of the disclosed technology, the ionosphere 208 three-dimensional electron density structure may be determined in a similar way as discussed above with respect to the air density, and may have an additional advantage in that certain satellite 226 constellations may orbit within the ionosphere 208 itself. Specifically, satellites 226 in Low Earth Orbit are within line of sight of many GNSS satellites 102. In certain example implementations of the disclosed technology, such satellites 226 may be utilized measure the density delay distance (D3) from the GNSS satellites 102 at high temporal rates. The availability of such GNSS signals 236 that traverse near the surface of the earth 201 may enable the system 200 to use tomography to determine the difference between the part of the ionosphere 208 that is lower than the receiver satellite 226, and the part that is above the satellite 226. Since the satellites 226 are near the peak of ionospheric electron density, this is very helpful for the mapping of the electrons in the ionosphere 208. The rays emanating from all directions within the desired three-dimensional substance create a naturally orthogonal set. Again, discussed above, a the predominantly horizontal signals (such as those received via radio occultation satellite 226) combined with predominantly vertical signals (such as may be received from an overhead GNSS satellite 202) may be utilized for the determination of the 3D ionospheric electron density.

[0088] In accordance with certain example implementations of the disclosed technology, a feedback loop may be utilized between the GNSS stations 104 and the processing system 108 For example, the data exchange system 106 may be utilized to provide a two-way communication between at least a portion of the receivers 104 and the processing system 108, which may allow the GNSS stations 104 to have more accurate positioning without the simplifying assumption discussed above. In an example implementation, accurate mass density data can be fed back into GNSS stations 104 to provide enhanced location accuracy.

[0089] Certain example implementations of the disclosed technology may provide the technical effects of preventing rare but large location errors. For example, a projected safety standard for a driverless vehicle is that it should not stray outside of a one or two-meter locational accuracy more than once in a hundred years. Certain example implementations of the disclosed technology may help reduce and/or eliminate such location errors.

[0090] Numerical testing was utilized to validate the processes describe herein. For example, simulated GNSS signals received at different receivers were utilized to determine that the density states and atmospheric fields could be recovered. The numerical tests used current estimates of location accuracy and realistic errors. The tests were conducted for a domain consisting of the US 48 contiguous states, and areas immediately around that domain. Only dry atmosphere data was used for the initial validation tests, as it was assumed that the electron density component could be removeable via L-band differential frequency response, as described above.

[0091] During the validation testing, certain elements, factors, and specifications of the system 200 were identified for refinement and/or weighting. For example, it became apparent during testing that the ability to provide highly accurately density states depends on the measurement of nearly orthogonal GNSS signals in the constituent small volumes. In accordance with certain example implementations of the disclosed technology, the orthogonality of the received GNSS signals may be determined, with associated calculations weighted according to the crossing angles. In certain example implementations of the disclosed technology, weighting may also be used to account for the vertical and/or horizontal offset between received GNSS signals, as such offsets may define a larger constituent volume, and hence, less resolution in the measurements.

[0092] The validation testing also provided information about horizontal and temporal resolution that the system 200 may provide for a given density of receivers 104. For example, previous weather models typically assimilate data for 12 hours to create an initial state. However, it was determined that the systems and methods described herein could support more frequent data, and may be utilized to create a new state every hour, which may in turn enable hourly weather predictions. There are many cases in the Earth's atmosphere, mainly related to thunderstorms and convection, where the hourly changes in atmospheric density, water vapor and associated weather are quite rapid. In the last decade, weather models and analysis have shown that the use of hourly observations greatly improve prediction of these important phenomena. Since the density inhomogeneities can be large (several percent), the associated errors in location could also be large, especially for GNSS satellites at low angles. In certain example implementations, the three-dimensional density fields may significantly improve location accuracy in convective areas.

[0093] Additionally, the validation testing was utilized to provide information regarding the preferred use of aircraft 218 receiver data at heights to eliminate errors due to water vapor variability. For example, it was determined that data from aircraft at higher levels of the atmosphere (e.g. typical flight levels of 10 km) could be preferable since almost all water vapor is concentrated below 10 km. Additionally, it was determined that data taken during aircraft take-off and landing could be utilized in specifying the full atmospheric state. The validation testing confirmed that highly detailed weather prediction could be achieved via such observations.

[0094] In accordance with certain example implementations of the disclosed technology, a three-dimensional density field of the dry atmosphere, together with assumptions of various complexity may be made to lead to the full momentum and thermodynamic structure of the atmosphere and ionosphere. According to an example implementation of the disclosed technology, hydrostatic and balance assumptions may be utilized to create the full state.

[0095] Certain example implementations of the disclosed technology may utilize the methodology of Observing System Simulation Experiments (OSSEs), such as described in "Observing System Simulation Experiments at the National Centers for Environmental Prediction," by Masutani, M., et al. (2010), J. Geophys. Res., 115, D07101, doi: 10.1029/2009JD012528, the contents of which are incorporated herein by reference as if presented in full. According to an example implementation of the disclosed technology, OSSE may utilize analysis and statistics techniques to estimate the atmospheric conditions using the system 200 described herein. Certain example implementations of the disclosed technology may use OSSE techniques to determine the extent to which the system 200 could improve weather forecasts and/or associated GNSS location data.

[0096] According to an example implementation of the disclosed technology, an OSSE process may start with a hypothetical "true" atmosphere over the experimental domain. In example embodiment, a low-pressure system may be defined that extends from the surface of the earth 201 into the stratosphere. In certain example implementations, it may be assumed that the dynamic fields, the associated directional momentum components, the density, and the temperature may vary over the same volumes. Conservative estimates may be made regarding the number of GNSS signals that could be received from surface, aircraft and satellite receivers 104.

[0097] In accordance with certain example implementations of the disclosed technology, the "true" atmosphere may be used to estimate the D3 for the incoming GNSS signals, based only on the dry air density. In an example implementation, and for the "first guess," a constant temperature and pressure may be utilized on each of the horizontal pressure levels, for example, from the earth surface (at approximately 1 bar) to the top level in the stratosphere (at approximately 5 milli-bar). In certain example implementations, the input fields of D3 from all available GNSS stations 104, together with the initial first guess may be iterated using a "relaxation" search to arrive at an approximate solution. In certain example implementations, empirically determined search coefficients (such as the D3 increment) may provide a robust system that can converge toward a good approximation of the solution representing the true atmosphere.

[0098] In one example experiment conducted, about 53,000 received D3 distances were utilized, from which about 21,000 were received from land, 22,000 from aircraft, and about 9,000 from satellite Radio Occultation. The associated calculated D3 densities from such measurements may be considered realistic for about a one-hour period over the 48-state domain. An average error of 2 cm was used for the receivers, using a random number generator to make realistic simulations of noise, which is a very conservative choice for errors, since it is believed average errors can be reduced significantly below the 2-cm level. [0099] FIG. 3 A shows an initial surface wind and pressure "true" field. FIG. 3B shows calculated surface wind and pressure field from the Observing System Simulation Experiments (OSSEs) processes after a few hundred iterations. FIG. 4 A shows an initial wind and pressure "true" field at 5km elevation. FIG. 4B shows calculated 5km elevation wind and pressure field from the Observing System Simulation Experiments (OSSEs) processes after a few hundred iterations.

[00100] FIGs. 3A and 3B, and FIGs. 4A and 4B, for example, show results of the OSSE processes after a few hundred iterations, which took a few hours on a small PC. For example, FIG. 3A shows the initial surface wind and pressure "true" field, and FIG. 3B shows the field derived from the simulated D3 observations. The pressure and wind fields show a deep low in central US, with strong winds around it (no surface friction included). The low in the true field, as shown in FIG. 3A, is closely approximated by the low in the derived field, as shown in FIG. 3B. The random error inserted can be seen in the wind field away from the low, but the important structure is well duplicated at the surface by the experiments. FIGs. 4A and 4B show corresponding the pressure and wind fields at 5 km high in the atmosphere. For example, FIG. 4A shows the initial wind and pressure "true" field at 5km, and FIG. 4B shows the field at 5km derived from the simulated D3 observations. As demonstrated by these experimental simulation, the desired fields may be well recovered for the central US low using the techniques described herein, with random error perturbations added. The results of the OSSE experiments show that air density (and derived fields of wind, temperature and pressure) may be recovered from appropriately distributed D3 measurements. Certain example implementations may also be used for the recovery of three-dimensional electron density via a similar process using L-band differential measurements.

[00101] FIG. 5 depicts a block diagram of an illustrative computing device 500 according to an example implementation. Certain aspects of FIG. 5 may be embodied in a transmitter, satellite, ground station, receiver, data exchange system, central command and control system, and/or associated processing system (for example, the satellites 202, 204, 206, transmitters 102, GNSS stations 104, CCC 107, data exchange system 106, processing system 108, etc., as shown in FIGs. 1 and/or 2). According to one example implementation, the term "computing device," as used herein, may be a CPU, processor, or conceptualized as a CPU (for example, the CPU 502 of FIG. 5). In this example implementation, the computing device (CPU) may be coupled, connected, and/or in communication with one or more peripheral devices, such as display. In another example implementation, the term computing device, as used herein, may refer to a processor and associated components, for example, that may be installed and reside on satellites in a constellation (such one or more satellites 202, 204, 206 as shown in FIG. 2).

[00102] In certain example embodiments, a satellite may have various sensors connected to and in communication with a CPU (such as the CPU 502 of FIG. 5). For example, the sensors can include one or more of the following: a frequency specific monitor such as UV (Ultraviolet) and IR (infrared), a sensor for remote detection of surface temperature, a spectrometer (UV and/or IR, for example), an accelerometer, a camera or vision system for still and video capture, a gravimetric sensor, a radar, and a radio transceiver, among other possibilities. In certain example implementations, a satellite may embody a computing device that includes memory (such as memory 518 depicted in FIG. 5) to store data collected by one or more sensors.

[00103] In an example implementation, the computing device may output content to its local display and may transmit and receive messages via the antenna interface 510, the network connection interface 512, telephony subsystem 532, etc. In example implementation, the computing device may output content to an external display device (e.g., over Wi-Fi) such as a TV or an external computing system. It will be understood that the computing device 500 is provided for example purposes only and does not limit the scope of the various implementations of the communication systems and methods.

[00104] The computing device 500 of FIG. 5 includes a central processing unit (CPU) 502, where computer instructions are processed. Certain example implementations can include a display interface 504 that acts as a communication interface and provides functions for rendering video, graphics, images, and texts on the display. In certain example implementations of the disclosed technology, the display interface 504 may be directly connected to a local display, such as a touch-screen display associated with a mobile computing device. In another example implementation, the display interface 504 may be configured to provide content (for example, data, images, and other information as previously discussed) for an external/remote display that is not necessarily physically connected to the computing device 500. For example, a desktop monitor may be utilized for mirroring graphics and other information that is presented on a mobile computing device. In certain example implementations, the display interface 504 may wirelessly communicate, for example, via a Wi-Fi channel or other available network connection interface 512 to an external/remote display.

[00105] In an example implementation, the network connection interface 512 may be configured as a communication interface and may provide functions for rendering video, graphics, images, text, other information, or any combination thereof on the display. In one example, the computing device 500 may include a communication interface that may include one or more of: a serial port, a parallel port, a general-purpose input and output (GPIO) port, a universal serial bus (USB), a micro-USB port, a high definition multimedia (HDMI) port, a video port, an audio port, a Bluetooth™ port, a near-field communication (NFC) port, another like communication interface, or any combination thereof.

[00106] According to an example implementation of the disclosed technology, the computing device 500 may include a keyboard interface 506 that provides a communication interface to a keyboard. In one example implementation, the computing device 500 may include a pointing device interface 508 for connecting to a presence-sensitive input interface. According to certain example implementations of the disclosed technology, the pointing device interface 508 may provide a communication interface to various devices such as a touch screen, a depth camera, etc.

[00107] The computing device 500 may be configured to use an input device via one or more of input/output interfaces (for example, the keyboard interface 506, the display interface 504, the pointing device interface 508, the antenna interface 510, the network connection interface 512, camera interface 514, sound interface 516, etc.,) to allow capturing information into the computing device 500. The input device may include a mouse, a trackball, a directional pad, a track pad, a touch-verified track pad, a presence-sensitive track pad, a presence-sensitive display, a scroll wheel, a digital camera, a digital video camera, a web camera, a microphone, a sensor, a smartcard, and the like. Additionally, the input device may be integrated with the computing device 500 or may be a separate device. For example, the input device may be an accelerometer, a magnetometer, a digital camera, a microphone, and/or an optical sensor.

[00108] Certain example implementations of the computing device 500 may include an antenna interface 510 in communication with an antenna. For example, certain satellites (such as the satellites as shown in FIGs. 1-2) may include one or more antennas with high gain capabilities. Certain example implementations of the antenna interface 510 can include one or more of: a receiver, analog-to-digital converter, sampler, buffers, memory, and memory. Certain example implementations can include a network connection interface 512 that provides a communication interface to a network. In certain implementations, a camera interface 514 may act as a communication interface to provide functions for capturing digital images from a camera. In certain implementations, a sound interface 516 may be provided as a communication interface for converting sound into electrical signals using a microphone and for converting electrical signals into sound using a speaker. According to example implementations, a random- access memory (RAM) 518 is provided, where computer instructions and data may be stored in a volatile memory device for processing by the CPU 502.

[00109] According to an example implementation, the computing device 500 includes a readonly memory (ROM) 520 where invariant low-level system code or data for basic system functions such as basic input and output (I/O), startup, or reception of keystrokes from a keyboard are stored in a non-volatile memory device. According to an example implementation, the computing device 500 includes a storage medium 522 or other suitable type of memory (e.g. such as RAM, ROM, programmable read-only memory (PROM), erasable programmable readonly memory (EPROM), electrically erasable programmable read-only memory (EEPROM), magnetic disks, optical disks, floppy disks, hard disks, removable cartridges, flash drives), where the files include an operating system 524, application programs 526 and content files 528 are stored. In accordance with certain example implementations of the disclosed technology, the application programs 526 can include one or more of programs to dynamically adjust the various parameters, such as one or more of: data packet size; data chunk size; data coding rate; number of re-transmits on uplink; number of tolerated timeouts; maximum transmission time; data compression parameters; bit rate factors; recovery factors; resource burden factors; batching; latency control; advance filtering; etc.

[00110] According to an example implementation, the computing device 500 includes a power source 530 that provides an appropriate alternating current (AC) or direct current (DC) to power components. According to an example implementation, the computing device 500 can include a telephony subsystem 532 that allows the device 500 to transmit and receive sound over a telephone network. The constituent devices and the CPU 502 communicate with each other over a bus 534. [00111] In accordance with an example implementation, the CPU 502 has appropriate structure to be a computer processor. In one arrangement, the computer CPU 502 may include more than one processing unit. The RAM 518 interfaces with the computer bus 534 to provide quick RAM storage to the CPU 502 during the execution of software programs such as the operating system application programs, and device drivers. More specifically, the CPU 502 loads computer-executable process steps from the storage medium 522 or other media into a field of the RAM 518 to execute software programs. Content may be stored in the RAM 518, where the content may be accessed by the computer CPU 502 during execution. In one example configuration, the device 500 includes at least 128 MB of RAM, and 256 MB of flash memory.

[00112] The storage medium 522 itself may include a number of physical drive units, such as a redundant array of independent disks (RAID), a floppy disk drive, a flash memory, a USB flash drive, an external hard disk drive, thumb drive, pen drive, key drive, a High-Density Digital Versatile Disc (HD-DVD) optical disc drive, an internal hard disk drive, a Blu-Ray optical disc drive, or a Holographic Digital Data Storage (HDDS) optical disc drive, an external mini-dual inline memory module (DIMM) synchronous dynamic random access memory (SDRAM), or an external micro-DIMM SDRAM. Such computer readable storage media allow the device 500 to access computer-executable process steps, application programs and the like, stored on removable and non-removable memory media, to off-load data from the device 500 or to upload data onto the device 500. A computer program product, such as one utilizing a communication system may be tangibly embodied in storage medium 522, which may comprise a machine- readable storage medium.

[00113] According to one example implementation, the terms computing device or mobile computing device, as used herein, may be a central processing unit (CPU), controller or processor, or may be conceptualized as a CPU, controller or processor (for example, the CPU processor 502 of FIG. 5). In yet other instances, a computing device may be a CPU, controller or processor combined with one or more additional hardware components. In certain example implementations, the computing device operating as a CPU, controller or processor may be operatively coupled with one or more peripheral devices, such as a navigation system. In another example implementation, the term computing device, as used herein, may refer to a satellite processor. [00114] In accordance with an example implementation of the disclosed technology, the various available resources may be selected and utilized to affect and/or improve communications efficiency. For example, an effective bit-rate can be influenced by many factors including but not limited to noise, re-transmission of data due to errors, transmit power, error correction encoding/decoding, etc. Certain example implementations of the disclosed technology may provide a flexible, tunable optimization process to enhance the efficiency of the communication and control of the various resources.

[00115] FIG. 6 is flow diagram of a method 600, according to an example implementation of the disclosed technology. In block 602, the method 600 includes receiving, at a processing system, and from a plurality of Global Navigation Satellite Systems (GNSS) stations, navigation data corresponding to computed positions of the plurality of GNSS stations. In block 604, the method 600 includes receiving, at the processing system, GNSS transmitter information including one or more of: GNSS satellite positional information; GNSS clock information; and GNSS transmission frequency information. In block 606, the method 600 includes determining, based at least in part on the received navigation data and the received GNSS transmitter information, ionosphere and atmosphere refractivity associated with navigation data. In block 608, the method 600 includes determining, based at least in part on the refractivity, three- dimensional (3D) density states of the atmosphere and ionosphere corresponding to intersections of two or more GNSS signals. In block 610, the method 600 includes calculating, by the processing system, and based on the determined 3D density states, data fields of a model representing the three-3D density states. In block 612, the method 600 includes transmitting position adjustment data to calibrate a navigation position of at least one of the plurality of the GNSS stations based at least in part on the calculated data fields of the model.

[00116] FIG. 7 is flow diagram of a method 700, according to an example implementation of the disclosed technology. In block 702, the method 700 includes receiving, at a processing system, and from a plurality of Global Navigation Satellite Systems (GNSS) stations, navigation data corresponding to computed positions of the plurality of GNSS stations. In block 704, the method 700 includes receiving, at the processing system, GNSS transmitter information including one or more of: GNSS satellite positional information; GNSS clock information; and GNSS transmission frequency information. In block 706, the method 700 includes determining, based at least in part on the received navigation data and the received GNSS transmitter information, ionosphere and atmosphere refractivity associated with navigation data. In block 708, the method 700 includes determining, based at least in part on the refractivity, three- dimensional (3D) density states of the atmosphere and ionosphere corresponding to intersections of two or more GNSS signals. In block 710, the method 700 includes calculating, by the processing system, and based on the determined 3D density states, data fields of a model representing the three-3D density states. In block 712, the method 700 includes determining at least one atmospheric state based at least in part on the calculated data fields of the model. In block 714, the method 700 includes transmitting data representing the at least one atmospheric state.

[00117] The method 600 and/or the method 700 may include additional steps or elements as discussed below.

[00118] In certain example implementations, determining the ionosphere and atmosphere refractivity can include computing density delay components for dry air, water, and ions.

[00119] In certain example implementations, variational assimilation may be utilized to compute the density delay components.

[00120] In certain example implementations, at least a portion of the plurality of GNSS stations can include one or more local processors, and the computed positions may be determined by the corresponding local processors.

[00121] According to an example implementation of the disclosed technology, the GNSS transmitter information may be received from a central command and control (CCC) system. In certain example implementations, the GNSS transmitter information may be received from one or more GNSS stations. In some example implementations, the GNSS transmitter information may be received from a service provider.

[00122] In certain example implementations, the navigation data and the GNSS transmitter information may correspond to two or more GNSS satellites.

[00123] According to an example implementation of the disclosed technology, determining the 3D density states include determining intersection angles of the two or more GNSS signals. [00124] According to an example implementation of the disclosed technology, updating the model can include applying intersection angle weightings, such as trigonometric functions of the corresponding intersection angles.

[00125] In certain example implementations, determining the 3D density states can include determining dynamic changes associated with the atmosphere and ionosphere.

[00126] According to an example implementation of the disclosed technology, determining the 3D density states can include predicting future density states associated with the atmosphere and ionosphere.

[00127] According to an example implementation of the disclosed technology, the data fields of the model may be refined by feedback.

[00128] In certain example implementations, determining the ionosphere and atmosphere refractivity can include numerically solving the three-dimensional state using the mathematical technique of minimization. For example, an optimal three-dimensional density may be derived by minimizing a least square difference between the integrated refraction effects of density over one-dimensional rays to find the three-dimensional structure that is optimal.

[00129] In certain example implementations, three-dimensional density fields may be recovered using refractive effects detected in a large number of receivers. In certain example implementations, a sufficient number of rays from the GNSS satellites to the receivers may be accounted for to resolve the field(s). In certain example implementations, the rays may have significant numbers crossing at substantially different angles. In other words, in certain example implementations, the rays may not traverse a single predominant direction or path.

[00130] In accordance with an example implementation of the disclosed technology, Radio Occultation data may be communicated to a central processor and combined with the other in- situ data (such as from ships, aircraft and land surface receivers) to provide field resolution using the additional ray paths.

[00131] In the preceding description, numerous specific details have been set forth. However, it is to be understood that implementations of the disclosed technology may be practiced without these specific details. In other instances, well-known methods, structures, and techniques have not been shown in detail in order not to obscure an understanding of this description. References to "one implementation," "an implementation," "example implementation," "various implementations," etc., indicate that the implementation(s) of the disclosed technology so described may include a feature, structure, or characteristic, but not every implementation necessarily includes the particular feature, structure, or characteristic. Further, repeated use of the phrase "in one implementation" does not necessarily refer to the same implementation, although it may.

[00132] As used herein, unless otherwise specified the use of the ordinal adjectives "first," "second," "third," etc., to describe a common object, merely indicate that different instances of like objects are being referred to, and are not intended to imply that the objects so described must be in a sequence, either temporally, spatially, in ranking, or in any other manner.

[00133] While certain implementations of the disclosed technology have been described regarding what is presently considered to be the most practical and various implementations, it is to be understood that the disclosed technology is not to be limited to the disclosed implementations, but on the contrary, is intended to cover various modifications and equivalent arrangements included within the scope of the appended claims. Although specific terms are employed herein, they are used in a generic and descriptive sense only and not for purposes of limitation.

[00134] This written description uses examples to disclose certain implementations of the disclosed technology, including the best mode, and to enable any person skilled in the art to practice certain implementations of the disclosed technology, including making and using any devices or systems and performing any incorporated methods. The patentable scope of certain implementations of the disclosed technology is defined in the claims, and may include other examples that occur to those skilled in the art. Such other examples are intended to be within the scope of the claims if they have structural elements that do not differ from the literal language of the claims, or if they include equivalent structural elements with insubstantial differences from the literal language of the claims.