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.     

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.     

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

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

Recent results on the Venus atmosphere from pioneer Venus radio occultations

Radio occultation measurements of the temperature structure of the Venus atmosphere have been obtained during seven occultation “seasons” extending from December 1978 to December 1983. Approximately 123 vertical profiles of temperature from about 40 km to about 85 km altitudes have been derived. Since these measurements cover latitudes from both poles to the equator, they have shown the latitudinal dependence of thermal structure. There is a smooth transition from the troposphere to the mesosphere at latitudes below about 45°, with the tropopause at about 56 km. The troposphere then rises to about 62 km in the “collar cloud” region between about 60° and 80° latitude, where a strong temperature inversion (up to 30 K) is present. In the polar areas, 80°–90°, the mesosphere becomes isothermal and there is no inversion. This latitudinal behavior is related to the persistent circulation pattern, in which a predominantly zonal retrograde motion at latitudes below 45° gradually changes to a circumpolar vortex at the “collar cloud” latitudes. Indeed, the radio occultation data have been used in a cyclostrophic balance model to derive zonal winds in the Venus atmosphere, which showed a mid-latitude (50°–55°) jet with a speed of about 120–140 ms^?1 at about 70 km altitude /1,2/. The observations obtained in 1983 and 1984 have shown that above the tropopause there is considerable temporal variability in the detailed thermal structure, suggesting that the persistent circulation pattern is subject to weather-like variability.     

Observation of wave structures in the upper atmosphere by means of radio holographic analysis of the radio occultation data

A radio holographic approach, developed by Pavelyev (1998), Hocke (1999), Igarashi (2000), is applied to observation of wave phenomena in the upper atmosphere using Global Positioning System — “Microlab-1” satellite (GPS/MET) radio occultation data. In the current state the radio holography approach uses the radar focused synthetic aperture principle to obtain high spatial resolution, and to remove the interference part corresponding to scattering from the upper ionosphere. High spatial resolution and accuracy of the radio halographic method is validated by means of revealing the weak signal reflected from the sea in the GPS/MET radio occultation data. The radio holographic method gives a new possibility to measure directly the vertical gradient of the electron density altitude profile in the D-layer using the radio occultation signal. The results of the application of radio holographic analysis to two GPS/MET occultation events (07 February 1997, No. 0447, 0158), in the D-region of the ionosphere, are discussed. Wave structures in the electron density concentration with a vertical spatial period of 1.4–6 km, and variations in the electron density gradient from ±5·10⁹ to ±8·10⁹ [1/(m³km)], have been retrieved from the D-layer data. The features observed in the vertical electron density profiles may be connected with breaking of gravity waves in the D-layer of the ionosphere.     

FY-3D卫星的北斗掩星分布特征与误差特性

利用2018年1－3月FY-3D卫星的掩星折射率数据,研究了北斗导航卫星系统的掩星分布特点、数据精度以及误差统计特征。北斗导航卫星系统同步静止轨道掩星沿卫星轨道呈弧状分布在南北两极地区,倾斜轨道掩星在东西半球低纬度地区分别形成一小一大两个空洞,中地球轨道掩星则全球均匀分布。北斗掩星折射率数据精度在探测核心区域,即12~32 km范围内,与ERA5再分析资料计算的折射率相比,平均偏差的标准差约为1.5%,在核心区外,标准差从1.5%逐渐增大到6%。静止轨道掩星的平均偏差在高层略大于倾斜轨道和中地球轨道掩星。下降掩星在20 km以上区域的标准差大于上升掩星,20 km以下区域小于上升掩星。高纬地区北斗掩星标准差最小,低纬地区最大,对流层中下层尤其明显。分析结果表明,北斗掩星的数据精度和误差特征与GPS掩星数据相似。      

Characteristics of the seasonal variation of the global tropopause revealed by COSMIC/GPS data

Using the Global Navigation Satellite System (GNSS) radio occultation observations from Formosa Satellite mission-3/Constellation Observing System for Meteorology, Ionosphere, and Climate (FORMOSAT-3/COSMIC) from 2007 to 2012, the climatological characteristics of the global tropopause was studied, with the following features identified. The overall results generally agree with previous studies. The tropopause has an obvious zonal structure, with more zonal characteristics in the Southern Hemisphere than the Northern Hemisphere. The vertical shape of the tropopause is sharp in the tropics and broad in the sub-tropical latitudes, with the sharpest latitudinal gradient in the mid-latitudes of both hemispheres. The global tropopause exists in a large range between 8 km and 17 km (or between 100 hPa and 340 hPa). The highest tropopause is over the South Asian monsoon regions for the entire year. The spatial structure of the tropopause in the polar region is of concentric structure, with an altitude between 7.5 km and 10 km. It is more symmetric in the Antarctic than the Arctic. Differing from other places, the height of the tropopause in the Antarctic is higher in winter as opposed to summer. The tropopause has distinct seasonal variability, especially in polar regions.     

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.     

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.     

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.     

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.     

A global survey of COSMIC ionospheric peak electron density and its height: A comparison with ground-based ionosonde measurements

With a network of ground-based ionosondes distributed around the world, the ionospheric peak electron density and its height measured by FORMOSAT-3/COSMIC satellites in terms of GPS radio occultation technique are extensively examined in this article. It is found that, in spite of the latitude, the mean values of the peak electron density measured by COSMIC satellites are systematically smaller than those observed by ground-based ionosondes. The discrepancy between them is dependent on the latitude, namely, it is small in low and mid-latitudes and large in high-latitude region. Moreover, statistical analysis shows that the slopes of the regression line that is best fitted to the scatter diagram of occultation-retrieved peak electron density (ordinate axis) versus ionosonde-observed peak density (abscissa axis) are universally less than one. This feature is believed to be the result of path average effect of non-uniform distribution of the electron density along the GSP ray during the occultation. A comparison between COSMIC-measured peak height and ionosonde-derived peak height hmF2 indicates that the former is systematically higher than the latter. The difference in the two can be as large as 20% or more in equatorial and low-latitude regions. This result implies that the peak height hmF2 derived from the virtual height through true height analysis based on Titheridge method seems to underestimate the true peak height. The correlation between COSMIC and ionosonde peak electron densities is analyzed and the result reveals that correlation coefficient seems to be dependent on the fluctuation of the occultation-retrieved electron density profile. The correlation will be higher (lower) for the electron density profiles with smaller (larger) fluctuations. This feature suggests that the inhomogeneous distribution of the electron density along the GPS ray path during the occultation plays an important role affecting the correlation between COSMIC and ionosonde measurements.     

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.     

GNSS reflectometry and remote sensing: New objectives and results

The Global Navigation Satellite System (GNSS) has been a very powerful and important contributor to all scientific questions related to precise positioning on Earth’s surface, particularly as a mature technique in geodesy and geosciences. With the development of GNSS as a satellite microwave (L-band) technique, more and wider applications and new potentials are explored and utilized. The versatile and available GNSS signals can image the Earth’s surface environments as a new, highly precise, continuous, all-weather and near-real-time remote sensing tool. The refracted signals from GNSS radio occultation satellites together with ground GNSS observations can provide the high-resolution tropospheric water vapor, temperature and pressure, tropopause parameters and ionospheric total electron content (TEC) and electron density profile as well. The GNSS reflected signals from the ocean and land surface could determine the ocean height, wind speed and wind direction of ocean surface, soil moisture, ice and snow thickness. In this paper, GNSS remote sensing applications in the atmosphere, oceans, land and hydrology are presented as well as new objectives and results discussed.     

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.     

地基行星大气掩星观测资料流程和数据整理

在中国火星探测萤火一号(YH-1)计划中, 包括了地基掩星观测反演火星大气的科研任务. 观测资料整理是反演流程的第一步. 本文描述了地基火星大气掩星观测处理软件系统的观测数据流程和观测资料整理模块,并详细介绍了观测资料整理模块的结构和功能, 其中包括时间系统转换、历表插值、坐标系变换、信号时延改正以及掩星平面建立. 利用行星数据系统公布的火星快车无线电科学数据和由SPICE得到的地球、火星历表以及火星快车的轨道数据, 结合本文的算法, 得到了一些实验结果.      

Radio occultation by Saturn's rings: Observations of structure and particle size with Voyager 1

Voyager 1 radio occultation study of Saturn's rings gives detailed information regarding the rings' radial structure and particle sizes. Structure within the rings is mapped to a radial resolution of few hundred m in the tenuous parts of ring C and the Cassini Division, and few km over ost of ring A. Fine resolution profiles reveal extremely sharp edges, very narrow gaps, and a host of wave phenomena. Particle size distributions obtained from occultation data within several ring regions are roughly consistent with an inverse cube power law with upper size cutoff in the 5 to 10 m radius range.     

Comparative analysis of Venusian ionosphere dual-frequency radio soundings with the satellites Venera-9, 10 and Pioneer-Venus

A comparison of the Venusian ionosphere electron density profiles obtained by the dual-frequency radio occultation method for the satellites Venera-9,10 (1975) and the Pioneer-Venus orbiter (1978–1979) has been carried out. It is shown that the general nature of the profiles, the main maximum heights and electron densities, the ionopause height positions determined with the satellites Venera-9, 10 on the one hand and Pioneer-Venus orbiter on the other hand are in good agreement. This fact testifies to the reliability of the experimental dual-frequency radio occultation data. An attempt to compare the radio occultation data with the direct measurements points out the essential contradictions between them, which need detailed analysis and discussion.     

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.     

基于FY-3C掩星数据偶发E层的研究

利用FY-3C极轨卫星提供的2014年6月至2015年5月的GPS无线电掩星数据,统计分析了全球范围内抽样频率为50 Hz的C/A码SNR扰动情况,进而对偶发E层进行了研究.结果表明:偶发E层在夏季半球中纬地区的扰动强度远远大于冬季半球同一纬度地区的扰动强度,偶发E层在纬度40°附近扰动明显增强;在E层100 km高度附近,Es层在10:00 LT和22:00 LT达到峰值;Es层在夏季半球的出现率明显高于冬季半球;FY-3C卫星的掩星观测结果与COSMIC系统的观测结果较一致,可以利用FY-3C卫星的掩星数据研究电离层偶发E层等的变化.