Estimating Total Electron Content Using 1,000+ GPS Receivers

A computer program uses data from more than 1,000 Global Positioning System (GPS) receivers in an Internet-accessible global network to generate daily estimates of the global distribution of vertical total electron content (VTEC) of the ionosphere. This program supersedes an older program capable of processing readings from only about 200 GPS receivers. This program downloads the data via the Internet, then processes the data in three stages. In the first stage, raw data from a global subnetwork of about 200 receivers are preprocessed, station by station, in a Kalman-filter-based least-squares estimation scheme that estimates satellite and receiver differential biases for these receivers and for satellites. In the second stage, an observation equation that incorporates the results from the first stage and the raw data from the remaining 800 receivers is solved to obtain the differential biases for these receivers. The only remaining error sources for which an account cannot be given are multipath and receiver noise contributions. The third stage is a postprocessing stage in which all the processed data are combined and used to generate new data products, including receiver differential biases and global and regional VTEC maps and animations

Multi-Cone Model for Estimating GPS Ionospheric Delays

The multi-cone model is a computational model for estimating ionospheric delays of Global Positioning System (GPS) signals. It is a direct descendant of the conical-domain model. A primary motivation for the development of this model is the need to find alternatives for modeling slant delays at low latitudes, where ionospheric behavior poses an acute challenge for GPS signal-delay estimates based upon the thin-shell model of the ionosphere

Generating high precision ionospheric ground-truth measurements

A method, apparatus and article of manufacture provide ionospheric ground-truth measurements for use in a wide-area augmentation system (WAAS). Ionospheric pseudorange/code and carrier phase data as primary observables is received by a WAAS receiver. A polynomial fit is performed on the phase data that is examined to identify any cycle slips in the phase data. The phase data is then leveled. Satellite and receiver biases are obtained and applied to the leveled phase data to obtain unbiased phase-leveled ionospheric measurements that are used in a WAAS system. In addition, one of several measurements may be selected and data is output that provides information on the quality of the measurements that are used to determine corrective messages as part of the WAAS system

Effects Of The Ionosphere On Ground-Based Detection Of The Global 21 CM Signal From The Cosmic Dawn And The Dark Ages

Detection of the global HI 21 cm signal from Cosmic Dawn and Epoch of Reionization is the key science driver for several ongoing ground-based and future ground/space-based experiments. The crucial spectral features in the global 21 cm signal (turning points) occur at low radio frequencies <100 MHz. In addition to the human-generated RFI, Earth's ionosphere drastically corrupts low-frequency radio observations from the ground. In this paper, we examine the effects of time-varying ionospheric refraction, absorption and thermal emission at these low radio frequencies and their combined effect on any ground-based global 21 cm experiment. It should be noted that this is the first study of the effect of a dynamic ionosphere on global 21 cm experiments. The fluctuations in the ionosphere are influenced by solar activity with flicker noise characteristics. The same characteristics are reflected in the ionospheric corruption to any radio signal passing through the ionosphere. As a result, any ground based observations of the faint global 21 cm signal are corrupted by flicker noise (or "$1/f$" noise, where "$f$" is the dynamical frequency) which scales as $\nu^{-2}$ (where $\nu$ is the frequency of observation) in the presence of a bright galactic foreground ($\propto \nu^{-s}$, where $s$ is radio spectral index). Hence, the calibration of the ionosphere for any such experiment is critical. Any attempt to calibrate the ionospheric effects will be subject to the inaccuracies in the current ionospheric measurements using GPS ionospheric measurements, riometer measurements, ionospheric soundings, etc. Even considering an optimistic improvement in the accuracy of GPS-TEC (Total Electron Content) measurements, we conclude that the detection of the global 21 cm signal below 100 MHz is best done from above the Earth's atmosphere in orbit of the Moon.Comment: 15 pages, 11 figures. Submitted to The Astrophysical Journal. This is an updated version after addressing the comments from the refere

Real-time detection of tsunami ionospheric disturbances with a stand-alone GNSS receiver. A preliminary feasibility demonstration

It is well known that tsunamis can produce gravity waves that propagate up to the ionosphere generating disturbed electron densities in the E and F regions. These ionospheric disturbances can be studied in detail using ionospheric total electron content (TEC) measurements collected by continuously operating ground-based receivers from the Global Navigation Satellite Systems (GNSS). Here, we present results using a new approach, named VARION (Variometric Approach for Real-Time Ionosphere Observation), and estimate slant TEC (sTEC) variations in a real-time scenario. Using the VARION algorithm we compute TEC variations at 56 GPS receivers in Hawaii as induced by the 2012 Haida Gwaii tsunami event. We observe TEC perturbations with amplitudes of up to 0.25 TEC units and traveling ionospheric perturbations (TIDs) moving away from the earthquake epicenter at an approximate speed of 316 m/s. We perform a wavelet analysis to analyze localized variations of power in the TEC time series and we find perturbation periods consistent with a tsunami typical deep ocean period. Finally, we present comparisons with the real-time tsunami MOST (Method of Splitting Tsunami) model produced by the NOAA Center for Tsunami Research and we observe variations in TEC that correlate in time and space with the tsunami waves

Modeling of Upper Atmospheric Responses to Acoustic-Gravity Waves Generated by Earthquakes and Tsunamis

Recent studies have shown that upper atmospheric observations can be used to examine the properties of acoustic and gravity waves (AGWs) induced by natural hazards (NHs), including earthquakes and tsunamis (e.g., Komjathy et al., Radio Sci., 51, 2016). In addition to rapid processing, analysis, and retrieval of the AGW signals in data, the need remains to investigate a broad parameter space of atmospheric and ionospheric state observables for the robust constraint of coupled and nonlinear processes. Here, we present several earthquake/tsunami-atmosphere-ionosphere case studies that demonstrate the possibility to detect AGWs and constrain the characteristics of their sources. Direct numerical simulations of the triggering and wave dynamical processes, from Earth\u27s interior to the exobase, are carried out based on coupled forward seismic wave and tsunami propagation models and our state-of-the-art nonlinear neutral atmosphere and ionosphere models MAGIC and GEMINI (Zettergren and Snively, JGR, 120, 2015)

The 2009 Samoa and 2010 Chile Tsunamis as Observed in the Ionosphere using GPS Total Electron Content

Ground‐based Global Positioning System (GPS) measurements of ionospheric total electron content (TEC) show variations consistent with atmospheric internal gravity waves caused by ocean tsunamis following two recent seismic events: the Samoa earthquake of 29 September 2009 and the Chile earthquake of 27 February 2010. Both earthquakes produced ocean tsunamis that were destructive to coastal communities near the epicenters, and both were observed in tidal gauge and buoy measurements throughout the Pacific Ocean. We observe fluctuations in TEC correlated in time, space, and wave properties with these tsunamis using the Jet Propulsion Laboratory’s Global Ionospheric Mapping software. These TEC measurements were band‐pass filtered to remove ionospheric TEC variations with wavelengths and periods outside the typical range for tsunamis. Observable variations in TEC appear correlated with the tsunamis in some locations (Hawaii and Japan), but not in others (Southern California or near the epicenters). Where variations are observed, the typical amplitude tends to be ∼0.1–0.2 TEC units for these events, on the order of ∼1% of the background TEC value. These observations are compared to estimates of expected tsunami‐driven TEC variations produced by Embry Riddle Aeronautical University’s Spectral Full Wave Model, an atmosphere‐ionosphere coupled model, and are found to be in good agreement. Significant TEC variations are not always seen when a tsunami is present, but in these two events the regions where a strong ocean tsunami was observed coincided with clear TEC observations, while a lack of clear TEC observations coincided with smaller sea surface height amplitudes. There exists the potential to apply these detection techniques to real‐time GPS TEC data, providing estimates of tsunami speed and amplitude that may be useful for early warning systems

Processing Satellite Data for Slant Total Electron Content Measurements

A method, system, and apparatus provide the ability to estimate ionospheric observables using space-borne observations. Space-borne global positioning system (GPS) data of ionospheric delay are obtained from a satellite. The space-borne GPS data are combined with ground-based GPS observations. The combination is utilized in a model to estimate a global three-dimensional (3D) electron density field

GLAS Spacecraft Pointing Study

Science requirements for the GLAS mission demand that the laser altimeter be pointed to within 50 m of the location of the previous repeat ground track. The satellite will be flown in a repeat orbit of 182 days. Operationally, the required pointing information will be determined on the ground using the nominal ground track, to which pointing is desired, and the current propagated orbit of the satellite as inputs to the roll computation algorithm developed by CCAR. The roll profile will be used to generate a set of fit coefficients which can be uploaded on a daily basis and used by the on-board attitude control system. In addition, an algorithm has been developed for computation of the associated command quaternions which will be necessary when pointing at targets of opportunity. It may be desirable in the future to perform the roll calculation in an autonomous real-time mode on-board the spacecraft. GPS can provide near real-time tracking of the satellite, and the nominal ground track can be stored in the on-board computer. It will be necessary to choose the spacing of this nominal ground track to meet storage requirements in the on-board environment. Several methods for generating the roll profile from a sparse reference ground track are presented