武汉与海口地区对流层顶的变化特征

通过对武汉和海口地区2000-2004年Radiosonde常规观测数据的分析,讨论了当地对流层顶温度与高度的年变化特征,并与热带地区站点的结果进行了比较.各等压面温度的扰动、温度递减率及对流层厚度的年变化特征显示,对流层顶更多地受到了平流层的影响.通过对武汉地区2006年1月11-15日Radiosonde加密观测数据的分析,对武汉地区对流层顶周日变化的特征进行了探讨,结果表明武汉地区对流层顶的温度存在显著的周日变化.     

国际参考电离层模型输出参数峰值高度性能

电离层峰值高度H_mF₂是描述电离层形态的重要参数之一,国际参考电离层模型IRI-2016中融入了大量电离层测高仪和无线电掩星探测数据,用以提升H_mF₂的预测精度.本文利用太阳活动低年（2007—2010年）气象、电离层和气候卫星联合观测系统COSMIC探测数据描述全球范围内COSMIC H_mF₂的三维形态变化,对比分析了IRI-2016与IRI-2012模型的预测结果,同时分析了IRI-2016模型输出H_mF₂的性能.结果表明,IRI模型在中高纬度地区的输出结果高于COSMIC反演结果,而赤道及低纬地区则大都偏低.与IRI-2012模型相比,IRI-2016模型的输出结果在夜间至凌晨时段呈现较为明显的纬向梯度变化且大部分区域输出值偏高,但在白天时段赤道附近区域的输出值大都偏低.上述结果为电离层IRI模型的完善提供了一定参考.     

TIMED卫星探测的全球大气温度分布及其与经验模式的比较

利用TIMED卫星遥感探测的全球温度分布与NRLMSISE-00大气经验模式进行了对比研究.研究表明,在中间层下部以下的高度范围内,经验模式与卫星探测的大气温度分布有很好的一致性.但是比较发现,在中层顶区域,经验模式的计算结果与实测结果有较大的差异.卫星探测表明,在春分季节的低纬地区中层顶区存在稳定的逆温层,但是经验模式不能给出低纬地区春分季节中间层逆温层的分布特征.卫星观测表明在全球范围内中层顶有两个非常不同的高度,一个处于100km附近,另一个处于85km附近,但是经验模式不能给出这一中层顶高度的分布特征.同时在低热层,经验模式计算的温度分布与卫星遥感的探测结果有很大的差异.     

基于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层等的变化.     

利用掩星温度数据推算大气月平均纬向风场

为了弥补中高层大气风场探测数据的不足,对2007年COSMIC全球掩星在20～60km高度内的温度数据进行网格化插值,利用梯度风计算方法,计算得到20～60km高度内的月平均纬向风场,分别与ECMWF再分析数据和HWM07模式进行对比分析.结果表明,利用COSMIC温度数据计算得出的纬向风与ECMWF的纬向风十分接近,而与HWM07结果有一定差异,但总体变化趋势一致.与ECMWF再分析数据相比,利用COSMIC温度数据计算的月平均纬向风场平均偏差为-1.50～-0.08m·^-1,标准差为1.50～11.95m·^-1.与HWM07模式风场相比,利用COSMIC温度计算的月平均纬向风场平均偏差为-0.83～1.21m·^-1,标准差为3.69～11.14m·^-1.     

利用COSMIC数据构建的全球垂直TEC图的验证研究

在日固坐标系(地磁纬度和地方时)下, 累积地方时过去24h的COSMIC(Constellation Observing System for Meteorology Ionosphere andClimate)观测资料, 通过对110$sim$750km高度范围内的电子密度进行数值积分得到各掩星点的垂直TEC值, 进而利用Kriging方法插值产生近实时的全球地方时MAGLat2.5°×2h的COSMIC TEC图. 利用2008年1月1日至2010年6月30日共30个月的COSMIC数据, 逐日构建COSMICTEC图, 将其与全球导航卫星系统服务组织(International GNSS Service,IGS)发布的全球电离层TEC图(Global Ionospheric Maps, GIMs)以及OSTM/JASON-2卫星高度计观测值分别进行比对,证明利用COSMIC掩星资料构建全球电离层垂直TEC图是可行的.     

ROPP反演软件算法及其精度分析

介绍了ROPP反演软件中无线电掩星反演的算法与精度分析. 采用COSMIC卫星2008年1 月1日全天的附加相位数据, 反演得到折射率、温度、压强与湿度等参数, 并与CDAAC 相应结果进行对比. 实验结果表明, 在30km高度以下, 折射率、压强和湿度的相对 误差在2%以内, 温度误差不超过2K.     

对流层顶变化对上对流层/下平流层臭氧分布的影响

上对流层和下平流层(UT／LS)，位于8-25km高度之间，是大气中一个很特殊的区域．大部分的臭氧分布在下平流层，在下平流层臭氧的含量发生一个很小的变化，就会对气候和地面的紫外辐射产生很大的影响．而作为气象参数的对流层顶，是充分混合、缺乏臭氧的上对流层和层结稳定、臭氧丰富的下平流层之间的边界或过渡层，其变化对臭氧总量和分布有直接和明显的影响．本文使用二维模式模拟研究对流层顶变化对臭氧在UT／LS分布的影响．模拟结果表明对流层顶的季节变化对UT/LS的臭氧分布有明显的影响，臭氧的局地变化可以超过10％在冬季北半球中纬度对流层顶高度升高1km时，模式结果表明对臭氧分布的影响比较显著，局地变化可超过6％，但是对臭氧总量的影响较小，变化不超过5DU，小于观测资料统计分析的结果。     

基于COSMIC的平流层重力波参数分析

基于COSMIC卫星观测的2006年12月29日到2008年1月3日30°-40°N纬度内的温度剖面,分别利用垂直滑动窗、双滤波器和单滤波器三种方法计算低平流层重力波的扰动和势能,获得重力波扰动和势能随高度、经度的分布以及多时间尺度变化特性,分析重力波扰动势能与背景温度及风场的变化趋势和特点.比较三种方法得到的结果发现：垂直滑动窗方法只能去除大垂直尺度的背景,无法抑制小尺度的扰动,其得到的结果误差较大;双滤波器法对温度剖面中的大尺度背景和小尺度扰动都能很好地抑制;单滤波器法得到的重力波扰动中基本不包含垂直方向的大尺度背景,但是包含一些小垂直尺度的扰动.因此,对于垂直波长为10km左右的重力波,采用双滤波器法合适;如果需要得到小尺度重力波的变化特性,采用单滤波器法合适.采用双滤波器法无法得到势能随高度的变化,而采用单滤波器法能够给出每月势能随高度的分布.对30°-40°N纬度内的重力波参数进行统计分析得到重力波扰动、势能与背景温度和水平风场的关系.     

急流与低层大气重力波能量的相关性研究

通过分析武汉(30.5ºN, 114ºE)上空2000至2004年的Radiosonde常规观测数据, 对当地对流层(1~10 km)与低平流层(18~25km)的某些惯性重力波特征进行了研究, 发现重力波能量与急流强度变化之间存在显著的相关. 同时分析了海口(20ºN, 114ºE)(2000至2004年)、北京(40ºN, 116ºE)(2001年12月至2003年2 月)的Radiosonde观测数据并进行了比较, 发现较高纬度地区的相关性明显大于低纬度地区. 通过对武汉地区2006年1月5天Radiosonde加密观测数据的分析, 讨论了重力波动能与势能的高度变化, 进一步发现纬向风的垂直剪切与急流中心的分布分别与重力波势能和动能的强弱分布相对应.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.