Evaluation of uncertainty in gravity wave potential energy calculations through GPS radio occultation measurements

The application of the Global Positioning System (GPS) radio occultation (RO) method to the atmosphere enables the determination of height profiles of temperature, among other variables. From these measurements, gravity wave activity is usually quantified by calculating the potential energy through the integration of the ratio of perturbation and background temperatures between two given altitudes in each profile. The uncertainty in the estimation of wave activity depends on the systematic biases and random errors of the measured temperature, but also on additional factors like the selected vertical integration layer and the separation method between background and perturbation temperatures. In this study, the contributions of different parameters and variables to the uncertainty in the calculation of gravity wave potential energy in the lower stratosphere are investigated and quantified. In particular, a Monte Carlo method is used to evaluate the uncertainty that results from different GPS RO temperature error distributions. In addition, our analysis shows that RO data above 30 km height becomes dubious for gravity waves potential energy calculations.     

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.     

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.     

Characterisation of residual ionospheric errors in bending angles using GNSS RO end-to-end simulations

Global Navigation Satellite System (GNSS) radio occultation (RO) is an innovative meteorological remote sensing technique for measuring atmospheric parameters such as refractivity, temperature, water vapour and pressure for the improvement of numerical weather prediction (NWP) and global climate monitoring (GCM). GNSS RO has many unique characteristics including global coverage, long-term stability of observations, as well as high accuracy and high vertical resolution of the derived atmospheric profiles. One of the main error sources in GNSS RO observations that significantly affect the accuracy of the derived atmospheric parameters in the stratosphere is the ionospheric error. In order to mitigate the effect of this error, the linear ionospheric correction approach for dual-frequency GNSS RO observations is commonly used. However, the residual ionospheric errors (RIEs) can be still significant, especially when large ionospheric disturbances occur and prevail such as during the periods of active space weather. In this study, the RIEs were investigated under different local time, propagation direction and solar activity conditions and their effects on RO bending angles are characterised using end-to-end simulations. A three-step simulation study was designed to investigate the characteristics of the RIEs through comparing the bending angles with and without the effects of the RIEs. This research forms an important step forward in improving the accuracy of the atmospheric profiles derived from the GNSS RO technique.     

Atmosphere sounding by GPS radio occultation: First results from TanDEM-X and comparison with TerraSAR-X

On 21 June 2010 the TerraSAR-X satellite was joined by the TanDEM-X satellite. A Global Positioning System (GPS) radio occultation (RO) experiment using the twin satellites has been carried out to estimate the precision of GPS atmospheric soundings. For the Day Of Year (DOY) 330–336, 2011, we analyze phase and amplitude data recorded by GPS receivers separated by a few hundred meters in a low earth orbit and derive collocated atmospheric refractivity profiles. In the altitude range 10–20 km the standard deviation between TerraSAR-X and TanDEM-X refractivity does not exceed 0.15%. The standard deviation is rapidly increasing for lower and higher altitudes; close to the surface and at an altitude of 30 km the standard deviation reaches 0.8% and 0.5%, respectively. Systematic deviations between TerraSAR-X and TanDEM-X refractivity in the considered altitude range (0–30 km) are negligible. The results confirm the anticipated high precision of the GPS RO technique. However, the difference in the retrieved refractivity in the lower troposphere for different Open Loop (OL) signal tracking parameters, altered onboard TanDEM-X for DOY 49–55, 2012, calls for an in depth analysis. At the moment we can not exclude that a potential bias in the OL Doppler model introduces a bias in our retrieved refractivity at altitudes <8$<8$ km.     

利用COSMIC数据分析全球对流层顶温度和高度的变化特性

利用COSMIC卫星无线电掩星观测的温度数据分析了全球范围的对流层顶参量的变化特性.COSMIC多颗卫星连续观测的温度数据有很好的全球覆盖和很好的垂直分辨率,是分析全球范围对流层顶变化的重要探测数据.经过数据筛选,得到日平均温度的全年分布数据.采用大气温度递减率判据得到对流层顼温度和对流层顶高度的全球分布和变化.结果表明,全球范国内的对流层顶温度和对流层顶高度分布的基本结果与利用无线电探空资料和欧洲气象局再分析数据得到的结果基本一致;对流层项温度和对流层顶高度呈现明显的季节变化,赤道附近和中纬度地区的对流层顶高度变化与高纬地区的对流层顶高度变化出现反相:对流层顶温度和高度的变化呈现明显的南北半球不对称性.     

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).     

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.     

基于STK的GPS/LEO无线电掩星技术仿真研究

GPS/LEO无线电掩星技术反演地球大气参数剖面已经具有较高的精度. 国外开展了多个GPS/LEO掩星项目, 但中国还尚未深入进行相关的实验, 这制约了中国掩星技术的发展. 本文提出基于STK进行GPS/LEO掩星技术研究的方法; 根据GPS/LEO掩星的原理, 推导出掩星事件发生的条件和掩星切点的计算公式; 利用STK对掩星过程进行模拟, 得到掩星数据. 在大气球对称假设和大气模型已知的条件下, 反演得到中性大气折射指数. 通过比较模型和反演数据, 表明反演数据精度较高, 验证了利用STK模拟GPS/LEO掩星实验方法的可行性.      

Comparison of RO tropopause height based on different tropopause determination methods

The tropopause is an important boundary in the Earth’s atmosphere, and has been the subject of close attention from atmosphere and climate researchers. To monitor the global tropopause using radio occultation (RO) data, there are two primary methods, one is the widely used temperature lapse rate method, and the other is the bending angle covariance transform method which is unique to RO data. We use FengYun3-C (FY3C) and Meteorological Operational Satellite Program (MetOp) RO data and European Centre for Medium-Range Weather Forecasts (ECMWF) operational analysis data to determine differences in RO tropopause height calculated by these two methods. We compute biases of the RO lapse rate tropopause height (LRTH) and the RO bending angle tropopause height (BATH) relative to the ECMWF LRTH. The dependences of the tropopause height biases on tropopause height (TPH) retrieval method, latitude, season and RO mission are investigated. The results indicate that BATH show a consistent 0.8–1.2 km positive bias over the tropics and high latitude regions compared with LRTH, however, over 25° to 40° latitude in both hemisphere, BATH results are less stable. Furthermore, the mean bias between BATH and LRTH displayed a different symmetrical characteristic from 2017.12 to 2018.2 (DJF) compared to 2018.6–2018.8 (JJA). However, except for some bias over Antarctica, the mean value of both LRTH and BATH show a similar tropopause variation, indicating the consistency of both methods.     

Case study of inclined sporadic E layers in the Earth’s ionosphere observed by CHAMP/GPS radio occultations: Coupling between the tilted plasma layers and internal waves

We have used the radio occultation (RO) satellite data CHAMP/GPS (Challenging Minisatellite Payload/Global Positioning System) for studying the ionosphere of the Earth. A method for deriving the parameters of ionospheric structures is based upon an analysis of the RO signal variations in the phase path and intensity. This method allows one to estimate the spatial displacement of a plasma layer with respect to the ray perigee, and to determine the layer inclination and height correction values. In this paper, we focus on the case study of inclined sporadic E (E_s) layers in the high-latitude ionosphere based on available CHAMP RO data. Assuming that the internal gravity waves (IGWs) with the phase-fronts parallel to the ionization layer surfaces are responsible for the tilt angles of sporadic plasma layers, we have developed a new technique for determining the parameters of IGWs linked with the inclined E_s structures. A small-scale internal wave may be modulating initially horizontal E_s layer in height and causing a direction of the plasma density gradient to be rotated and aligned with that of the wave propagation vector k. The results of determination of the intrinsic wave frequency and period, vertical and horizontal wavelengths, intrinsic vertical and horizontal phase speeds, and other characteristics of IGWs under study are presented and discussed.     

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.     

The use of ionospheric tomography and elevation masks to reduce the overall error in single-frequency GPS timing applications

Signals from Global Positioning System (GPS) satellites at the horizon or at low elevations are often excluded from a GPS solution because they experience considerable ionospheric delays and multipath effects. Their exclusion can degrade the overall satellite geometry for the calculations, resulting in greater errors; an effect known as the Dilution of Precision (DOP). In contrast, signals from high elevation satellites experience less ionospheric delays and multipath effects. The aim is to find a balance in the choice of elevation mask, to reduce the propagation delays and multipath whilst maintaining good satellite geometry, and to use tomography to correct for the ionosphere and thus improve single-frequency GPS timing accuracy. GPS data, collected from a global network of dual-frequency GPS receivers, have been used to produce four GPS timing solutions, each with a different ionospheric compensation technique. One solution uses a 4D tomographic algorithm, Multi-Instrument Data Analysis System (MIDAS), to compensate for the ionospheric delay. Maps of ionospheric electron density are produced and used to correct the single-frequency pseudorange observations. This method is compared to a dual-frequency solution and two other single-frequency solutions: one does not include any ionospheric compensation and the other uses the broadcast Klobuchar model. Data from the solar maximum year 2002 and October 2003 have been investigated to display results when the ionospheric delays are large and variable. The study focuses on Europe and results are produced for the chosen test site, VILL (Villafranca, Spain). The effects of excluding all of the GPS satellites below various elevation masks, ranging from 5° to 40°, on timing solutions for fixed (static) and mobile (moving) situations are presented. The greatest timing accuracies when using the fixed GPS receiver technique are obtained by using a 40° mask, rather than a 5° mask. The mobile GPS timing solutions are most accurate when satellites at lower elevations continue to be included: using a mask between 10° and 20°. MIDAS offers the most accurate and least variable single-frequency timing solution and accuracies to within 10 ns are achieved for fixed GPS receiver situations. Future improvements are anticipated by combining both GPS and Galileo data towards computing a timing solution.     

Comparison of GPS S/C orbits determined from GPS and SLR tracking data

The last two Global Positioning System (GPS) spacecraft launched on August 30, 1993 and March 10, 1994 carried identical laser retroreflector arrays that allow laser ranging from ground stations. Tracking of these targets is a low priority for most of the ground stations due to the present mandatory support of missions with no alternative tracking. The available data is sparse and concentrated primarily on the first of the two spacecraft, GPS-35 (PRN 5). Despite this, analysis of this data set indicates that there is great potential for engineering and scientific experiments in a synergistic application of Satellite Laser Ranging (SLR) and radiometric GPS techniques.     

基于BDS/GPS和经验模型的区域电离层建模

基于北斗卫星导航系统（BDS）和全球定位系统（GPS）实测电离层穿刺点（IPP）数据，结合国际参考电离层（IRI）经验模型历史数据，提出一种对区域二维电离层总电子含量（TEC）进行高精度建模的方法.针对缺乏穿刺点的区域内短时间电离层建模时精度较低且各时段穿刺点空间分布不同的问题，该方法使用IRI模型在建模区域内均匀添加虚拟穿刺点数据，并根据与实测穿刺点的距离，使用构造的权重计算公式赋予其动态权重值，通过加权最小二乘法进行球谐模型参数解算.与欧洲定轨中心（CODE）发布的全球电离层图（GIM）进行数据比对发现，相对于只使用BDS/GPS实测穿刺点数据的建模方法，利用本文建模方法计算获得的垂直总电子含量（VTEC）值对缺乏实测穿刺点的区域精度有明显的提升.      

Coastal sea level measurements using a single geodetic GPS receiver

This paper presents a method to derive local sea level variations using data from a single geodetic-quality Global Navigation Satellite System (GNSS) receiver using GPS (Global Positioning System) signals. This method is based on multipath theory for specular reflections and the use of Signal-to-Noise Ratio (SNR) data. The technique could be valuable for altimeter calibration and validation. Data from two test sites, a dedicated GPS tide gauge at the Onsala Space Observatory (OSO) in Sweden and the Friday Harbor GPS site of the EarthScope Plate Boundary Observatory (PBO) in USA, are analyzed. The sea level results are compared to independently observed sea level data from nearby and in situ tide gauges. For OSO, the Root-Mean-Square (RMS) agreement is better than 5 cm, while it is in the order of 10 cm for Friday Harbor. The correlation coefficients are better than 0.97 for both sites. For OSO, the SNR-based results are also compared with results from a geodetic analysis of GPS data of a two receivers/antennae tide gauge installation. The SNR-based analysis results in a slightly worse RMS agreement with respect to the independent tide gauge data than the geodetic analysis (4.8 cm and 4.0 cm, respectively). However, it provides results even for rough sea surface conditions when the two receivers/antennae installation no longer records the necessary data for a geodetic analysis.     

Correlation of GPS signal fades due to ionospheric scintillation for aviation applications

Deep and frequent Global Positioning System (GPS) signal fading due to strong ionospheric scintillation is a major concern for GPS-guided aviation in equatorial areas during high solar activity. A GPS aviation receiver may lose carrier tracking lock under deep fading, and a lost channel cannot be used for position calculation until lock is reestablished. Hence, frequent loss of lock due to frequent fading can significantly reduce the availability of GPS aviation. However, the geometric diversity of the satellites can mitigate scintillation impact on GPS aviation depending on the correlation level of deep fades between satellites. This paper proposes a metric to measure the correlation level of two fading channels from the perspective of GPS aviation. Using this metric, the satellite-to-satellite correlation is studied based on real scintillation data. The low satellite-to-satellite correlation shown in this paper envisions notable availability benefit from the geometric diversity of satellites under strong scintillation. In addition, this paper proposes a way to generate correlated fading processes with arbitrary correlation coefficients. Using this correlated fading process model, the availability of Localizer Performance with Vertical guidance (LPV)-200 under severe scintillation scenarios is analyzed. The result emphasizes the importance of a fast reacquisition capability of an aviation receiver after a brief outage, which is not currently mandated by the aviation receiver performance standards.     

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.     

Performance assessment of real-time orbit determination for the Haiyang-2D using onboard BDS-3/GPS observations

The Global Navigation Satellite System (GNSS) receivers equipped on the Haiyang-2D (HY-2D) satellite is capable of tracking the signals of both the third generation of BeiDou satellite navigation System (BDS-3) and the Global Positioning System (GPS), which make it feasible to assess the performance of real-time orbit determination (RTOD) for the HY-2D using onboard GNSS observations. In this study, the achievable accuracy and convergence time of RTOD for the HY-2D using onboard BDS-3 and GPS observations are analyzed. Benefiting from the binary-offset-carrier (BOC) modulation, the BDS-3 C1X signal includes less noise than the GPS C1C signal, which has the same signal frequency and chipping rate. The root mean squares (RMS) of the noises of C1X and C1C code measurements are 0.579 m and 1.636 m, respectively. Thanks to a ten-times higher chipping rate, the code measurements of BDS-3 C5P, GPS C1W and C2W are less noisy. The RMS of code noises of BDS-3 C5P, GPS C1W, and C2W are 0.044 m, 0.386 m, and 0.272 m, respectively. For the HY-2D orbit, the three-dimensional (3D) and radial accuracies can reach 31.8 cm and 7.5 cm with only BDS-3 observations, around 50 % better than the corresponding accuracies with GPS. Better performance of the BDS-3 in RTOD for the HY-2D is attributed to the high quality of its broadcast ephemeris. When random parameters are used to absorb ephemeris errors, substantial improvement is seen in the accuracy of HY-2D orbit with either BDS-3 or GPS. The 3D RMS of HY-2D orbit errors with BDS-3 and GPS are enhanced to 23.1 cm and 33.6 cm, and the RMS of the radial components are improved to 6.1 cm and 13.3 cm, respectively. The convergence time is 41.6 and 75.5 min for the RTOD with BDS-3 and GPS, while it is reduced to 39.2 and 27.4 min after the broadcast ephemeris errors are absorbed by random parameters. Overall, the achievable accuracy of RTOD with BDS-3 reaches decimeter level, which is even better than that with GPS, making real-time navigation using onboard BDS-3 observations a feasible choice for future remote sensing missions.     

An update global model of hmF2 from values estimated from ionosonde and COSMIC/FORMOSAT-3 radio occultation

In this paper, we present our recent work on developing an updated global model of the ionospheric F2 peak height hmF2 parameter by combining data from the Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC/FORMOSAT-3) radio occultation (RO) measurements and from the extended global ionosonde stations. In particular, 10 Chinese ionosonde stations’ data are newly introduced into this study. The modeling technique used is based on a two-layer empirical orthogonal function (EOF) expansion. Global distributions of hmF2 maps calculated using the newly constructed global model and the one provided by the International Reference Ionosphere model (IRI-ITU-R) are compared with the global distributions of hmF2 obtained by the COSMIC RO measurements and quantitative statistical analysis of the differences between the model results and those of the COSMIC RO measurements is made for the low (2008) and high (2012) solar activity years. The obtained average root-mean-square differences (RMSEs) for our model are 27.7 km (11.1%) and 31.0 km (9.8%), respectively for the years 2008 and 2012, whereas those for the IRI-ITU-R model are 39.9 km (16.9%) and 35.0 km (11.6%), respectively. Comparison of the results calculated both by our model and the IRI-ITU-R model with the digisonde observation is also made. The comparisons show that the newly constructed global hmF2 model can reproduce reasonably well the observations and perform better than IRI-ITU-R model.