基于X射线掩星探测反演行星大气密度

X射线掩星是一种常见的天文现象，基于X射线掩星探测的大气密度反演是一种涉及学科交叉的新方法，其通过处理高能X射线天体辐射源的掩星观测数据实现大气密度的反演，基本原理为X射线在大气中传播时，X射线光子被大气中的原子（包括分子中的原子）吸收和散射，从而导致X射线强度发生衰减，根据衰减后X射线信号的强度反演对应的密度廓线。本文根据X射线掩星探测的应用需求，论证了基于X射线掩星实现大气密度反演的新方法，重点介绍了光变曲线拟合和能谱拟合两种地球中高层大气密度反演算法，分析了X射线掩星探测反演大气密度的研究进展和研究方法，对基于X射线掩星反演大气密度的优点进行分析和讨论，进而对X射线掩星探测的应用场景进行展望。结果表明，作为一种新型中高层大气密度测量手段，X射线掩星探测可对中高层大气密度实现有效探测，弥补了目前中高层大气密度实测数据的不足。     

无线电掩星和激光雷达观测结果比较

简要阐述了激光雷达和ＧＰＳ掩星探测大气的基本原理，并利用１９９６年１０月１４日１６２５ ＵＴ时 Ｂ本 ＴＭＵ钠和瑞利散射激光雷Ｂ（３９．２３°Ｅ、 ３５．３７°Ｎ）资＃，反演得到 ３０－５５ ｋｍ高度范围内大气密度和温度的观测结果，与此同时通过激光雷达上空的 ＧＰＳ掩星观测到的掩星切点处的观测结果进行了比较．结果表明：两种方法获得测出的密度廓线符合得很好；温度廓线的变化趋势也基本一致，但在某些高度上有较大的偏离．文中最后对比较结果作了简要分析，讨论了重力波对大气密度和温度廓线的影响．     

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

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

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

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

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

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

利用COSMIC掩星数据监测电离层的异常变化

分析了COSMIC掩星数据反演电子密度的方法,利用实例研究反演方法的特点,并采用ISR非相干散射雷达获取的电子密度数据进行验证,进而反演了长三角区域SHAO（IGS）站上空在日全食和太阳风暴期间的电子密度廓线图. 通过与平静日期间电离层电子密度进行比较,发现日全食及太阳风暴导致电离层发生的异常变化,从而提出COSMIC掩星数据反演电子密度在监测电离层变化时所具有的优势.      

瑞利激光雷达反演中层大气温度算法

瑞利激光雷达是探测中层大气温度的重要设备,具有高时空分辨率且不存在探测盲区的优点.为充分利用激光雷达探测到的原始数据,改进了传统的Chanin-Haunchecorne方法,采用由均匀搜索生成温度初值的方式反演大气温度.位于北京延庆（40.3°N,116.2°E）的瑞利激光雷达具有589nm和532nm两个通道,将589nm通道用于计算,532nm通道作为参考.在150~250K区间等间隔选取多个温度初值,利用589nm通道反演60~70km高度范围内的大气温度廓线和大气密度廓线,利用参考密度廓线选取准确的温度初值,反演得到准确的大气温度廓线.将589nm通道和532nm通道反演的温度廓线相比较,发现二者具有较高的一致性.改进后的方法有效利用了信噪比较差的光信号,使589nm通道的温度探测上限从60km提高到70km.      

改进的GPS掩星滑动频谱方法

在大气多路径和噪声的条件下,高斯白噪声会造成滑动频谱方法获得的弯曲角与真值之间的偏离,无法获得较好的反演结果。为此,文章提出了一种改进的滑动频谱方法,即利用信号的振幅和谱能量信息对滑动频谱方法进行修正,削弱了信号中噪声的影响,与真值较为接近。分别用改进前后两种方法对中华卫星三号计划（Constellation Observing System for Meteorology, Ionosphere and Climate,COSMIC）的掩星进行反演,并将其折射率计算结果和通过全谱反演方法获得的折射率一起,与欧洲中期天气预报中心（European Centre for Medium Range Weather Forecasts,ECMWF）的分析场资料进行了统计比较。结果表明：改进的滑动频谱方法删除了信号中的部分噪声,减少了系统偏差;与全谱反演方法进行比较,发现两者具有相当的反演精度。     

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

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

K波段微波高光谱辐射计的晴空定标方法

K波段微波高光谱辐射计应用于大气探测方向,能够减少地表至对流层内的大气廓线误差和提高廓线的垂直分辨率.针对平均辐射温度在宽大气频谱和不同天空仰角下存在浮动的问题,结合所统计的北京地区往年廓线数据和晴空条件下的Liebe大气吸收模型（1993）,提出了一种微波高光谱下改进的晴空定标方法.相比于传统的晴空定标方法,该方法利用先验的气象数据对定标过程进行优化计算,有利于提高平均辐射温度的预测精度.误差分析结果表明当K波段微波高光谱辐射计采用改进的晴空定标方法时,能够有效减少平均辐射温度对其定标误差的影响.      

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.     

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.     

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.     

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.     

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.     

A comprehensive study on middle atmospheric thermal structure over a tropic and sub-tropic stations

The general characteristics of middle atmospheric thermal structure have been studied by making use of the Rayleigh lidar data collected over the period of about four years (1998–2001). Here, the data has been used from two different stations in the Indian sub-continent in tropics (Gadanki; 13.5°N, 79.2°E) and in sub-tropics (Mt. Abu; 24.5°N, 72.7°E). The observed monthly mean temperature profiles are compared with different model atmospheres (CIRA-86 and MSISE-90). We observed, the mean temperature profiles have closer agreement with MSISE-90 than CIRA-86. The temperature profiles measured by lidar and HALOE satellite overpass nearby lidar site are generally in agreement with each other. The systematic and statistical errors in deriving temperature are found to be uniform for both the stations, as 1 K at 50 km, 3 K at 60 km and 10 K at 70 km. The special features of mesospheric temperature inversion (MTI) and double stratopause structure (DBS) are also addressed for both the stations.     

Seasonal variations in vertical distribution of meteor decay time as observed from meteor radars at 8.5°N and 80°N

The decay times of meteor radar echoes have been used for decades to investigate characteristics of the mesosphere and lower thermosphere (MLT) region. As the meteor echo decay time depends on background atmospheric parameters, in the present communication, we examine the seasonal variation of the vertical distributions of underdense meteor echo decay times with respect to echo strength. Observations from two similar radars located at two distinct geographical locations, Thumba (8.5°N, 77°E) and Eureka (80°N, 85.8°W) were used for the present study. Here, the radar received signal power is categorized into strong and weak echoes and vertical profiles of their decay times are constructed. It has been noticed that the monthly mean decay time vertical profile turning altitude (i.e., inflection point) varies in the range of 80–87?km of altitude depending on latitude. The turning altitude is observed at relatively lower heights in the winter than in summer at both the latitudes. The present analysis shows that the meteor decay time below the mean turning altitude follows a decreasing trend with decreasing altitude, which is quite distinct to the behaviour of ambipolar diffusion. It is also observed that there is a difference in mean decay time of strong and weak echoes below 90?km of altitude, which is very prominently seen at lower altitudes. This difference shows a seasonal pattern at high latitude, but does not show any seasonal variation at low latitude. The present results are discussed in light of current understanding of the meteor decay time.     

Proposed reference models for ozone

Ozone reference models are proposed here similar to the Keating and Young 1985 models which were prepared for the new COSPAR International Reference Atmosphere. This paper updates tables provided in the Keating and Young ozone model, giving improved monthly zonal mean total column ozone in 10° latitude increments, improved monthly zonal mean ozone volume mixing ratios (ppmv) from 20 to 0.003 mb in 10° latitude increments, and conversion tables providing ozone vertical structure in other units. Also, a new table is provided giving ozone vertical structure as a function of altitude (from 25 to 80 km), latitude, and month. The models are based on measurements from six contemporary satellite instruments.