Estimated errors in a global gravity wave climatology from GPS radio occultation temperature profiles

In a previous paper by Schmidt et al. (2008), from CHAllenging Minisatellite Payload (CHAMP) Global Positioning System (GPS) radio occultation data, a comparison was made between a Gaussian filter applied to the “complete” temperature profile and to its “separate” tropospheric and stratospheric height intervals, for gravity wave analyses. It was found that the separate filtering method considerably reduces a wave activity artificial enhancement near the tropopause, presumably due to the isolation process of the wave component. We now propose a simple approach to estimate the uncertainty in the calculation of the mean specific wave potential energy content, due exclusively to the filtering process of vertical temperature profiles, independently of the experimental origin of the data. The approach is developed through a statistical simulation, built up from the superposition of synthetic wave perturbations. These are adjusted by a recent gravity wave (GW) climatology and temperature profiles from reanalyses. A systematic overestimation of the mean specific wave potential energy content is detected and its variability with latitude, altitude, season and averaging height interval is highlighted.     

A comparative and numerical study of effects of gravity waves in small miss-distance and miss-time GPS radio occultation temperature profiles

The Global Positioning System (GPS) Radio Occultation (RO) technique has global coverage and is capable of generating high vertical resolution temperature profiles of the upper troposphere and lower stratosphere with sub-Kelvin accuracy and long-term stability, regardless of weather conditions. In this work, we take advantage of the anomalously high density of occultation events at the eastern side of the highest Andes Mountains during the initial mission months of COSMIC (Constellation Observing System for Meteorology, Ionosphere and Climate). This region is well-known for its high wave activity. We choose to study two pairs of GPS RO, both containing two occultations that occurred close in time and space. One pair shows significant differences between both temperature profiles. Numerical simulations with a mesoscale model were performed, in order to understand this discrepancy. It is attributed to the presence of a horizontal inhomogeneous structure caused by gravity waves.     

Higher order ionospheric propagation effects on GPS radio occultation signals

With the increasing number of remote sensing satellites using the GPS radio occultation technique for atmospheric sounding, the estimation of higher order ionospheric effects and their mitigation have become relevant and important. Due to long ionospheric limb paths, GPS signals are strongly affected by ionospheric refraction during radio occultation. Standard dual-frequency GPS measurements may be used to estimate the first order term of the refractive index. However, non-linear terms such as the second and third order ionospheric terms and ray path bending effects are not considered in occultation measurements so far. Analysing selected CHAMP–GPS occultation events different higher order ionospheric terms are estimated and their effects on dual-frequency range estimation and total electron content (TEC) estimation are discussed. We have found that the separation between the GPS L1 and L2 ray paths exceeds the kilometer level during occultation for a vertical TEC level of more than 160 TEC units. Corresponding errors in the GPS dual-frequency range estimation and TEC estimation are found to exceed the meter and 10 TEC units level, respectively.     

Temperature comparison between CHAMP radio occultation and TIMED/SABER measurements in the lower stratosphere

Global Positioning System (GPS) receiver on the CHAllenging Mini-satellite Payload (CHAMP) and the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) instrument, one of four on board the TIMED satellite, provide middle atmosphere temperature profiles by Radio Occultation (RO) and limb viewing infrared emission measurements, respectively. These temperature profiles retrieved by two different techniques in the stratosphere are compared with each other using more than 1300 correlative profiles in March, September and December 2005. The over-all mean differences averaged over 15 and 35 km are approximately −2 K and standard deviation is less than 3 K. Below 20 km of altitude, relatively small mean temperature differences ∼1 K are observed in wide latitudinal range except for June (during the SABER nighttime observation). In the middle to low latitudes, between 30°S and 30°N, the temperature difference increases with height from ∼0–1 K at 15 km, to ∼−4 K at 35 km of altitude. Large temperature differences about −4 to −6 K are observed between 60°S and 30°N and 31–35 km of altitude for all months and between 0° and 30°N below 16 km during June (nighttime).     

GRAS radio occultation on-board of Metop

The GRAS radio occultation instrument is flying on Metop-A and belongs to the EPS (EUMETSAT Polar System). GRAS observes GPS satellites in occultation. Within this work, validation of GRAS closed-loop bending angle data against co-located ECMWF profiles extracted from model fields and occultations from the COSMIC constellation of radio occultation instruments is shown. Results confirm the high data quality and robustness, where GRAS shows lower bending angle noise against ECMWF than COSMIC and in terms of occultations per day, one GRAS ≈ two COSMIC satellites. This is partly due to the operational setup of EPS. For the investigation we focus on two observation periods where updates in the ECMWF (March 2009) and COSMIC processing (October 2009) have improved the statistics further. Bending angles biases agree to within 0.5% against ECMWF and to within 0.1% against COSMIC after the updates for altitudes between 8 and 40 km. In addition, we also analyze the impact of the Metop orbit processing on the derived GRAS bending angle data, where different GPS and Metop orbit solutions are analyzed. Results show that a batch based orbit processing would improve in particular the bending angle bias behavior at higher altitudes. Requirements for the operational processing of GRAS data are briefly outlined, options to ease the use of other positioning system satellites in the near future are discussed. A simplified analysis on the observation of several of these systems, e.g. GPS and Galileo, from one platform shows that about 16% of occultations are found within 300 km, ±3 h, thus providing similar information. A constellation of 2 GRAS like instruments would have only about 10% close-by.     

A theoretical assessment of the feasibility of potential Lunar Reconnaissance Orbiter radio occultation observations of the lunar ionosphere

The existence of a “dense” lunar ionosphere has been controversial for decades. Positive ions produced from the lunar surface and exosphere are inferred to have densities that are $?$${10}^6$ – ${10}^7$ m^?3 near the surface and smaller at higher altitudes, yet electron densities derived from radio occultation measurements occasionally exceed these values by orders of magnitude. For example, about 4% of the single-spacecraft radio occultation measurements from Kaguya/SELENE were consistent with peak electron densities of $~$$3\times {10}^8$ m^?3. Space plasmas should be neutral on macroscopic scales, so this represents a substantial discrepancy. Aditional observations of electron densities in the lunar ionosphere are critical to resolving this longstanding paradox. Here we theoretically assess whether radio occultation observations using two-way coherent S-band radio signals from the Lunar Reconnaissance Orbiter (LRO) spacecraft could provide useful measurements of electron densities in the lunar ionosphere. We predict the uncertainty in a single LRO radio occultation measurement of electron density to be $~$$3\times {10}^8$ m^?3, comparable to occasional observations by Kaguya/SELENE of a dense lunar ionosphere. Thus an individual profile from LRO is unlikely to reliably detect the lunar ionosphere; however, averages of multiple ($~$10) LRO profiles acquired under similar geophysical and viewing conditions should be able to make reliable detections. An observing rate of six ingress occultations per day ($~$2000 per year) could be achieved with minimal impact on current LRO operations. This rate compares favorably with the 378 observations reported from the single-spacecraft experiment on Kaguya/SELENE between November 2007 and June 2009. The large number of observations possible for LRO would be sufficient to permit wide-ranging investigations of spatial and temporal variations in the poorly understood lunar ionosphere. These findings strengthen efforts to conduct such observations with LRO.     

Measurements of the upper troposphere and lower stratosphere during tropical cyclones using the GPS radio occultation technique

Water vapour transport to the upper troposphere and lower stratosphere by deep convective storms affects the radiation balance of the atmosphere and has been proposed as an important component of climate change. The aim of the work presented here is to understand if the GPS radio occultation technique is useful for characterization of this process. Our assessment addresses the question if severe storms leave a significant signature in radio occultation profiles in the upper troposphere/lower stratosphere. Radio occultation data from the Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) were analyzed, focusing on two particular tropical cyclones with completely different characteristics, the hurricane Bertha, which formed in the Atlantic Basin during July 2008 and reached a maximum intensity of Category 3, and the typhoon Hondo, which formed in the south Indian Ocean during 2008 reaching a maximum intensity of Category 4. The result is positive, suggesting that the bending angle of a GPS radio occultation signal contains interesting information on the atmosphere around the tropopause, but not any information regarding the water vapour. The maximum percentage anomaly of bending angle between 14 and 18 km of altitude during tropical cyclones is typically larger than the annual mean by 5–15% and it can reach 20% for extreme cases. The results are discussed in connection to the GPS radio occultation receiver which will be part of the Atomic Clock Ensemble in Space (ACES) payload on the International Space Station.     

A gravity wave analysis near to the Andes Range from GPS radio occultation data and mesoscale numerical simulations: Two case studies

Global maps of potential wave energy per unit mass, recently performed with the Global Positioning System (GPS) Radio Occultation (RO) technique and different satellite missions (CHAMP and SAC-C since 2001, GRACE and COSMIC since 2006) revealed in Argentina, at the eastern side of the highest Andes Mountains, a considerable wave activity (WA) in comparison with other extra-tropical regions. The main gravity wave (GW) sources in this natural laboratory are deep convection (mainly during late Spring and Summer), topographic forcing and geostrophic adjustment.     

Global ionospheric scintillations revealed by GPS radio occultation data with FY3C satellite before midnight during the March 2015 storm

Global observations of S4 amplitude scintillation index by the GPS Occultation Sounder (GNOS) on FengYun-3 C (FY3C) satellite reveal global dynamic patterns of a strong pre-midnight scintillations in F-region of the ionosphere during the St. Patrick’s Day geomagnetic super storm of 17–19 March 2015. The observed strong scintillations mainly occurred in the low latitudes, caused by equatorial plasma bubbles. During the main storm phase (March 17), the scintillations were first triggered in the New Zealand sector near 160°E longitudes, extending beyond 40°S dip latitude. They were also enhanced in the Indian sector, but significantly suppressed in East Asia near 120°E longitude and in Africa around 30°E longitude. During the initial recovery phase (March 18–19), the global scintillations were seldom observed in GNOS data. During the later recovery phase (after March 19), the scintillations recovered to the pre-storm level in Indian, African, and American sectors, but not in East Asian and any of Pacific sectors. These results closely correlate with observations of the density depletion structures by the Communication/Navigation Outage Forecasting System (C/NOFS) satellite, and ground-based instruments. Such consistency indicates reliability of our scintillation sensing approach even in a case-by-case comparison study. The prompt penetration electric field and disturbance dynamo electric field are suggested as the main factors that control the enhancement and inhibition of the scintillations during the storm, respectively.     

An introduction of mountain-based GPS radio occultation experiments in China

The mountain-based GPS radio occultation is a novel approach to lower atmospheric profiling. The experiments of the mountain-based GPS radio occultation were conducted on the top of Mt. Yaogu (29.38°N, 113.68°E, ∼1240 m) on December 17, 2003, and on the top of Mt. Jiugong (29.39°N, 114.65°E, ∼1550 m) on July 24, 2004. Based on these observation data, the scientific data processing software has been developed and is used to retrieve successfully the atmospheric refractivity profiles. The validation experiment was performed on the top of Mt. Wuling (40.60°N, 117.48°E, ∼2118 m) during August 1–29, 2005. Collocated automatic weather station and the radiosondes nearby were operated simultaneously for the comparison campaign. Results show that the radio occultation technique obtained about 40 profiles every day with the receiver antenna pointing to the south. Comparisons show that the refractivity measured by occultation agree well with those by the radiosondes, but not well with those by the automatic weather station due to their much different geographic locations of measurements. Results of these experiments suggest that the mountain-based GPS radio occultation is an economic reliable novel technique monitoring temporal and spatial variations of local lower atmospheric environments.     

Global monitoring of tropospheric water vapor with GPS radio occultation aboard CHAMP

The global positioning system radio occultation (GPS RO) technique provides a powerful tool for atmospheric sounding which requires no calibration, is not affected by clouds, aerosols or precipitation, and provides an almost uniform global coverage. The paper deals with application of GPS RO measurements from CHAllenging Minisatellite Payload (CHAMP) for the retrieval of tropospheric water vapor profiles. CHAMP RO data are available since 2001 with up to 200 high resolution atmospheric profiles per day. We introduce a new direct method for water vapor retrieval from GPS RO data. Additionally, a 1Dvar algorithm is used for this purpose. The so derived CHAMP water vapor profiles are validated with radiosonde data on a global scale. Here, both methods come to statistically comparable results revealing a negative bias of less than 0.1 g/kg and a standard deviation of less than 1 g/kg specific humidity in the mid troposphere. Potentials of CHAMP RO retrievals for monitoring the mean tropospheric water vapor distribution on a global scale are presented.     

The impact of spherical symmetry assumption on radio occultation data inversion in the ionosphere: An assessment study

‘Onion-peeling’ is a very common technique used to invert Radio Occultation (RO) data in the ionosphere. Because of the implicit assumption of spherical symmetry for the electron density (N(e)) distribution in the ionosphere, the standard Onion-peeling algorithm could give erroneous concentration values in the retrieved electron density profile. In particular, this happens when strong horizontal ionospheric electron density gradients are present, like for example in the Equatorial Ionization Anomaly (EIA) region during high solar activity periods. In this work, using simulated RO Total Electron Content (TEC) data computed by means of the NeQuick2 ionospheric electron density model and ideal RO geometries, we tried to formulate and evaluate an asymmetry level index for quasi-horizontal TEC observations. The asymmetry index is based on the electron density variation that a signal may experience along its path (satellite to satellite link) in a RO event and is strictly dependent on the occultation geometry (e.g. azimuth of the occultation plane). A very good correlation has been found between the asymmetry index and errors related to the inversion products, in particular those concerning the peak electron density NmF2 estimate and the Vertical TEC (VTEC) evaluation.     

Ionospheric F2-Layer critical frequency retrieved from COSMIC radio occultation: A statistical comparison with measurements from a meridional ionosonde chain over Southeast Asia

In this paper, we compared the F2-Layer critical frequency (foF2) derived from FORMOSAT-3/COSMIC radio occultation (RO) and ionosondes at Chiang Mai, Chumphon and Kototabang during the years 2008–2015 to evaluate the performance of COSMIC RO over Southeast Asia region. The results show that the time development of foF2 values derived from COSMIC RO generally agrees well with those from ionosonde measurements. However, the differences between the foF2 derived from COSMIC RO and that derived from ionosonde observations display latitudinal dependence. COSMIC RO tends to underestimate foF2 at Chiang Mai and Kototabang, which is near to the north EIA crest and the south one, respectively, while a little overestimate foF2 at Chumphon, which is close to the geomagnetic equator. COSMIC RO agrees best with ionosonde at Chumphon and worst at Chiang Mai. At each ionosonde station, the quality of COSMIC RO data degrades with the increase of solar activity. In addition, at the station Chiang Mai and Kototabang, COSMIC RO performs better in summer than in equinox and winter. Furthermore, the differences in foF2 derived from COSMIC RO and that from ionosonde measurements vary with local time, i.e., the differences in foF2 are generally smaller at night and larger in noontime when equatorial ionization anomaly (EIA) is well developed.     

Validation of FORMOSAT-3/COSMIC radio occultation electron density profiles by incoherent scatter radar data

Low Earth Orbiting satellites carrying a dual frequency GPS receiver onboard offer a unique opportunity to remote sensing of the global ionosphere on a continuous basis. No other profiling technique unifies profiling through the entire F2-layer with global coverage. The FORMOSAT-3/COSMIC data can make a positive impact on the global ionosphere study providing essential information about the height electron density distribution and particularly over regions that are not accessible with ground-based measuring instruments such as ionosondes and GPS dual frequency receivers. Therefore, it is important to verify occultation profiles with other techniques and to obtain experience in the reliability of their derivation. In the given study we present results of comparison of the electron density profiles derived from radio occultation measurements on-board FS-3/COSMIC and from the Kharkov incoherent scatter radar sounding.     

Variability of the upper troposphere and lower stratosphere observed with GPS radio occultation bending angles and temperatures

Recently, Lewis (2009) introduced a new method for the identification of tropopause heights (TPHs) from GPS radio occultation (RO) bending angles (α)$(\alpha )$. The method uses a covariance transform to identify transitions in a ln(α)$\mathit{ln}(\alpha )$ profile. Lewis validates the results with lapse rate tropopause (LRT) heights from one year of FORMOSAT-3/COSMIC data and radiosondes. In this study we apply the new method to the RO data sets from CHAMP/GRACE (2001–2009) and FORMOSAT-3/COSMIC (2006–2009). These results are the basis for TPH trend estimations for the time period between May 2001 and August 2009 (100 months) based on zonal monthly mean GPS RO data from CHAMP (2001–2008), GRACE (since 2006) and FORMOSAT-3/COSMIC (since 2006). Further, we compare the α$\alpha$ based TPH trends with LRT height trends and discuss the differences, which are largest in the subtropical regions (20°–40°) on both the northern and southern hemisphere. A global increase of the TPH between 5 and 9 m/yr is found for both methods and different data sets (CHAMP/GRACE alone and CHAMP/GRACE plus FORMOSAT-3/COSMIC). The results for the TPH trends are linked with bending angle and temperature trends in the upper troposphere and lower stratosphere region. Generally, an upper tropospheric warming (bending angle decrease) and a lower stratospheric cooling (bending angle increase) is noted.     

Location of layered structures in the ionosphere and atmosphere by use of GPS occultation data

A method is introduced to locate the layered structures in the atmosphere and ionosphere based on simultaneous observations of radio wave intensity and phase variations in trans-ionospheric satellite-to-satellite links. The method determines location of a tangent point on the trans-ionospheric ray trajectory where gradient of refractivity is perpendicular to the ray trajectory and influence of a layered structure on radio wave parameters is maximal. An estimate of the location of a layer can be obtained from a combination of the phase and intensity variations. This new technique was applied to measurements provided during FORMOSAT-3 and CHAMP radio occultation (RO) missions. For the considered RO events the location of the inclined plasma layer in the lower ionosphere is found and the electron density distribution is retrieved. The method is checked by measuring the location of the tangent point on the ray trajectory in the neutral gas in the atmosphere. The results showed a fairly good agreement.     

Assessment of precision in ionospheric electron density profiles retrieved by GPS radio occultations

The Constellation Observing System for Meteorology Ionosphere and Climate (COSMIC) is a six satellite radio occultation mission that was launched in April 2006. The close proximity of these satellites during some months after launch provides a unique opportunity to evaluate the precision of Global Positioning System (GPS) radio occultation (RO) retrievals of ionospheric electron density from nearly collocated and simultaneous observations. RO data from 30 consecutive days during July and August 2006 are divided into ten groups in terms of daytime or nighttime and latitude. In all cases, the best precision values (about 1%) are found at the F peak height and they slightly degrade upwards. For all daytime groups, it is seen that electron density profiles above about 120 km height exhibit a substantial improvement in precision. Nighttime groups are rather diverse: in particular, the precision becomes better than 10% above different levels between 120 and 200 km height. Our overall results show that up to 100–200 km (depending on each group), the uncertainty associated with the precision is in the order of the measured electron density values. Even worse, the retrieved values tend sometimes to be negative. Although we cannot rely directly on electron density values at these altitudes, the shape of the profiles could be indicative of some ionospheric features (e.g. waves and sporadic E layers). Above 200 km, the profiles of precision are qualitatively quite independent from daytime or latitude. From all the nearly collocated pairs studied, only 49 exhibited a difference between line of sight angles of both RO at the F peak height larger than 10°. After analyzing them we find no clear indications of a significant representativeness error in electron density profiles due to the spherical assumption above 120 km height. Differences in precision between setting and rising GPS RO may be attributed to the modification of the processing algorithms applied to rising cases during the initial period of the COSMIC mission.     

Secular variations of OH nightglow emission and of the OH intensity-weighted temperature induced by gravity–wave forcing in the MLT region

Secular variations of OH airglow (8,3) band and of the intensity-weighted temperature induced by a gravity wave packet at three latitudes in the northern hemisphere in the Mesosphere/Lower Thermosphere region are investigated using a 2-D, time-dependent, fully nonlinear OH chemistry–dynamics model. A net-cycled average survives due to wave transience and dissipation. The integrated OH (8,3) volume emission rates show a 20 percent increase induced by the wave packet. The corresponding intensity-weighted temperature shows an initial increase up to 2 percent then a gradual decrease for the remainder of the simulation time by the same wave packet. These secular variations could be mistaken as long-period or short-period waves in the airglow observations. Therefore, care must be taken when analyzing the data from observations.     

Comparison of the CHAMP radio occultation data with the Canadian advanced digital ionosonde in the Polar Regions

In this paper we present the results of the comparison of the retrieved electron density profiles of the Ionospheric Radio Occultation (IRO) experiment on board CHAMP (CHAllenging Minisatellite Payload), with the ground ionosonde profiles for the Polar Regions. IRO retrieved electron density profiles from CHAMP are compared with Canadian Advanced Digital Ionosonde (CADI) measurements at two vertical sounding stations well within the Polar Cap, Eureka (geog. 80°13′ N; 86°11′ W) and Resolute Bay (geog. 74°41′ N; 94°54′ W). We compared the ionospheric parameters such as the peak electron density of the F-layer (NmF₂) and the peak height of the F-layer (h_mF₂) for a 3-year period, 2004–2006. CHAMP derived NmF₂ shows reasonable agreement with the ionosonde retrieved NmF₂ for both the stations (0.76 and 0.71 correlation coefficient, for Eureka and Resolute Bay, respectively) whereas the h_mF₂ agreement is not that acceptable (0.25 and 0.37 correlation coefficient, respectively). The h_mF₂ from vertical sounding showed less spread than the CHAMP h_mF₂.