嫦娥一号和二号卫星微波探测仪数据的交叉定标与一致性

嫦娥一号（CE-1）和嫦娥二号（CE-2）卫星分别于2007年和2010年发射，其上搭载着相同的微波探测仪（CELMS）.由于在轨定标和探测时间不同，两星的CELMS数据存在差异.为了提高数据的一致性，基于月球极区永久阴影区（PSR），利用两步定标法对两星的CELMS数据进行交叉定标和时间校正.通过对不同纬度的亮温进行分析可知:两星的CELMS数据在3GHz，19.35GHz和37GHz通道的一致性良好；在中低纬度区域，CE-2卫星7.8GHz通道的亮温较CE-1卫星高，随着纬度的升高，差异逐渐变小并趋于一致.为避免经度跨度大对亮温的影响，选取云海盆地对交叉定标结果进行定量分析，结果表明两星在该区域的一致性良好，平均误差小于0.2K.由此可知，交叉定标和时间校正明显提高了两次观测数据的一致性，减小了由测量误差和时间变化引起的差异.      

Diviner红外温度与嫦娥卫星微波探测仪月表亮温的比较分析

提出使用高分辨率的Diviner辐射计红外温度作为月壤热传导的边界输入条件,使得月壤物理温度廓线的模拟更接近月表真实情况. 利用Diviner物理温度廓线模拟的微波亮温与嫦娥二号卫星微波探测仪（CELMS2）的实测亮温匹配,来分析CELMS2数据的有效性和稳定性,为月壤微波特性的反演奠定基础. 以哥白尼坑区域为例,对比分析了Diviner辐射计红外通道温度与CELMS1实测微波亮温的变化趋势. 选择月球虹湾地区内一位置点（45.1°-45.5°N,33.1°-33.5°W）,分析了CELMS2在该位置点内5个月球时刻的亮温数据. 通过对实测亮温与模拟亮温在经纬度和太阳高度角的匹配分析得到,在高频微波通道内两者的相关性达到96%以上,表明嫦娥二号微波探测仪亮温数据具有较高的稳定性,可以为月壤特性反演研究提供可靠数据来源.      

基于TechDemoSat-1卫星的GPS反射信号海面高度反演

针对星载GPS反射信号（GPS-R）海面测高的误差问题，基于星载GPS-R实测数据进行星载海面测高模型和误差修正模型的研究，并验证其有效性。利用TechDemoSat-1（TDS-1）数据，使用时延多普勒图（DDM）海面高度反演技术，着重分析了星载GPS-R海面高度反演中的各类误差，并建立了相应的误差模型。对星载GPS-R海面高度反演模型进行优化，采用DTU15全球平均海面模型、DTU全球海潮模型验证反演精度。结果证明:优化后反演模型得到的全球海面高度反演结果的平均绝对误差（MAD）为6.05 m，精度提高了约29%，有效提高了海面高度反演的精度。研究成果对于推广星载GNSS反射信号（GNSS-R）的海面测高应用具有一定的意义。      

观测路径和光谱分辨率对吸收性气溶胶指数的影响

吸收性气溶胶指数（AAI）在监测污染物方面有较好的应用,其反演结果受边界层高度、云高、相对湿度以及仪器指标和观测几何路径等因素影响.利用大气辐射传输模型MODTRAN,对辐射传输方程中多次散射计算的不同近似方法和不同观测几何角度下的AAI进行模拟,研究这些因素对AAI结果的作用,并分析仪器光谱分辨率对AAI探测结果的影响.结果表明,折衷考虑计算效率与精度的情况下,应选择8流近似的离散坐标法进行带有多次散射计算的辐射方程解析.AAI在不同浓度和不同类型气溶胶下随观测角度变化的趋势相同:相对方位角<120°时,AAI误差在太阳天顶角和卫星方位角均为40°～60°时最大;相对方位角在120°～180°时误差均较小;光谱分辨率对AAI反演结果无明显影响.      

AIRS观测资料研究全球平流层重力波特性

利用2012—2014年1月和7月AIRS（Atmospheric Infrared Sounder）第79通道的观测数据，分析了平流层重力波活动强弱的全球分布以及重力波发生频率的全球分布；分析了重力波活动随纬度和经度的变化特征，给出了重力波活动在全球范围内的热点区域及其活动强度；对比了白天与夜间的重力波活动强度及发生频率.研究表明重力波活动强度呈现出随纬度变化的特征，在低纬度地区（0°—30°），冬季半球重力波活动强度低，夏季半球重力波活动高；在中高纬度地区，冬季半球重力波活动强度高而夏季半球重力波活动强度低.在1月，全球重力波活动有4个突出的热点区域，分别为50°N附近欧洲大陆与大西洋交接地带、北美洲与大西洋交接地，20°S附近南美洲与大西洋交接地区、非洲与印度洋交接地区.在7月，重力波活动突出的地方为巴塔哥尼亚至南极半岛地区，以及50°S和75°E附近的印度洋区域.重力波活动强度在夜间大于白天，但是夜间的强重力波活动区域小于白天.      

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

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

HXMT卫星空间环境监测器及初步观测结果

硬X射线调制望远镜（HXMT）卫星是中国首个专门进行天文探测的空间科学实验卫星,运行于高度约550km、倾角约43°的低地球轨道.星载空间环境监测器为星上科学任务开展提供背景辐射实测资料.该监测器采用固体探测器望远镜系统和扇形阵列全新组合设计,可获取轨道空间高能质子和高能电子能谱、方向综合动态结果,给出更为全面的粒子辐射分布图像.初步探测结果显示,卫星运行轨道遭遇的带电粒子辐射集中分布在经度80°W-20°E,纬度0°-40°S的南大西洋异常区,粒子辐射在该区域表现出不同程度的方向差异分布,高能电子方向差异分布显著强于高能质子.2017年9月空间环境扰动期间,爆发的太阳质子事件并未对该轨道粒子辐射产生影响,而地磁活动导致该轨道穿越经度120°W-60°E,纬度40°-43°N的北美上空和经度60°-120°E,纬度43°-40°S的澳大利亚西南区域时遭遇增强粒子辐射影响,增强的粒子辐射表现出极强的方向分布.      

基于Fabry-Perot干涉仪气辉观测数据反演中性风速方法比较

利用半径法、完整傅里叶级数描述法和非线性回归法分别对安装在河北省境内国家天文台兴隆园区（40.2°N,117.4°E）的Fabry-Perot干涉仪（FPI）获得的地基中高层大气气辉观测数据进行风速反演.反演结果表明,半径法、完整傅里叶级数描述法和非线性回归法均能很好地反演风速,在误差范围内反演结果存在很好的一致性.相对于半径法和完整傅里叶级数描述法,非线性回归方法更完整地利用了观测数据信息,反演结果可信度较高.      

基于QZSS L1信号的岸基GNSS-R码延迟海面测高

为了探究GNSS L1波段中信号不同调制方式的GNSS-R码相位延迟海面的测高性能，开展了双天线岸基GNSS-R海面测高实验，收集了采样率为40 MHz的原始中频数据。利用自主开发的GNSS-R测高软件接收机和反演软件对实验数据进行了处理分析，同时获得了基于QZSS L1C/A码和L1C码的GNSS-R海面测高结果，并分别与岸基同步观测的雷达高度计的测高值进行对比，以对反演结果进行精度评定。实验结果表明:基于QZSS L1C/A码和L1C码的海面测高精度最优可达0.63 m和0.4 m，L1C码的延迟GNSS-R测高精度明显高于L1C/A码。此外，GNSS-R测高精度会随着卫星高度角的增加而有所增加。      

基于COSMIC掩星精密定轨数据的等离子体层电子含量研究

等离子体层是日地环境重要的组成部分.本文利用COSMIC掩星精密定轨数据经处理后得到的podTec文件获取等离子体层电子含量（PEC）对等离子体层进行研究.将podTec数据进行处理后获得的PEC（pod-PEC）和IRI-Plas经验模型提供的PEC （IRI-PEC）进行对比，发现pod-PEC与IRI-PEC符合得较好.在低（0°—20°）、中（20°—50°）、高（50°—90°）修正地磁纬度带下，分析了COSMIC在太阳活动极大年（2014年）3，6，9和12月的pod-PEC，得到如下结论:PEC随着纬度升高而逐渐减少，且3，9月PEC在中低纬关于磁赤道的南北对称性较好，6月北半球各纬度带的PEC明显高于南半球同一纬度带的值，而12月情况则完全相反，南半球中纬的PEC甚至会等于北半球低纬的PEC值；PEC在白天高而晚上低，高纬地区的PEC昼夜变化不明显；PEC具有明显的季节性.对于北半球，一年中PEC最大值出现在春季，冬秋季次之，夏季最低，具有明显的年度异常现象.      

D-region ionospheric disturbances associated with the Extremely Severe Cyclone Fani over North Indian Ocean as observed from two tropical VLF stations

The D-region ionospheric disturbances due to the tropical cyclone Fani over the Indian Ocean have been analysed using Very Low Frequency (VLF) radio communication signals from three transmitters (VTX, NWC and JJI) received at two low latitude stations (Kolkata-CUB and Cooch Behar-CHB). The cyclone Fani formed from a depression on 26th April, 2019 over the Bay of Bengal (Northeastern part of the Indian Ocean) and turned into an extremely severe cyclone with maximum 1-min sustained winds of 250 km/h on 2 May, 2019 which made landfall on 3 May, 2019. Out of six propagation paths, five propagation paths, except the JJI-CHB which was far away from the cyclone track, showed strong perturbations beyond $3\sigma$ level compared to unperturbed signals. Consistent good correlations of VLF signal perturbations with the wind speed and cyclone pressure have been seen for both the receiving stations. Computations of radio signal perturbations at CUB and CHB using the Long Wave Propagation Capability (LWPC) code revealed a Gaussian perturbation in the D-region ionosphere. Analysis of atmospheric temperature at different layers from the NASA’s TIMED satellite revealed a cooling effect near the tropopause and warming effects near the stratopause and upper mesosphere regions on 3 May, 2019. This study shows that the cyclone Fani perturbed the whole atmosphere, from troposphere to ionosphere and the VLF waves responded to the disturbances in the conductivity profiles of the lower ionosphere.     

Characterization of the low latitude plasma density irregularities observed using C/NOFS and SCINDA data

Complex electrodynamic processes over the low latitude region often result in post sunset plasma density irregularities which degrade satellite communication and navigation. In order to forecast the density irregularities, their occurrence time, duration and location need to be quantified. Data from the Communication/Navigation Outage Forecasting System (C/NOFS) satellite was used to characterize the low latitude ion density irregularities from 2011 to 2013. This was supported by ground based data from the SCIntillation Network Decision Aid (SCINDA) receivers at Makerere (Geographic coordinate 32.6°E, 0.3°N, and dip latitude ?9.3°N) and Nairobi (Geographic coordinate 36.8°E, ?1.3°N, and dip latitude ?10.8°N). The results show that irregularities in ion density have a daily pattern with peaks from 20:00 to 24:00 Local Time (LT). Scintillation activity at L band and VHF over East Africa peaked in 2011 and 2012 from 20:00 to 24:00 LT, though in many cases scintillation at VHF persisted longer than that at L band. A longitudinal pattern in ion density irregularity occurrence was observed with peaks over 135–180°E and 270–300°E. The likelihood of ion density irregularity occurrence decreased with increasing altitude. Analysis of C/NOFS zonal ion drift velocities showed that the largest nighttime and daytime drifts were in 270–300°E and 300–330°E longitude regions respectively. Zonal irregularity drift velocities over East Africa were for the first time estimated from L-band scintillation indices. The results show that the velocity of plasma density irregularities in 2011 and 2012 varied daily, and hourly in the range of 50–150 m s^?1. The zonal drift velocity estimates from the L-band scintillation indices had good positive correlation with the zonal drift velocities derived from VHF receivers by the spaced receiver technique.     

Periodic and aperiodic ozone variations in the middle and upper stratosphere

Umkehr, ozonesonde and satellite observations were used to determine the height/latitude distribution of the amplitude and phase of the periodic components of the variation of the ozone mixing ratio in the middle and upper stratosphere. The amplitude of the first (annual) harmonic is small in the subtropics and increases to a maximum at polar latitudes. It also increases with height in the mid and upper stratosphere to an apparent maximum just below the stratopause. The second (semi-annual) harmonic has an amplitude that is largest in tropical regions and in subpolar regions at a level of about 40 km. There seems to be very little ozone variation above 30 km with dominant periods close to the quasi-biennial period of total ozone observed in the tropics. The percent of the total variance of the ozone mixing ratio accounted for by the first harmonic is larger than 60 percent at all heights from 20° – 60° latitude in both hemispheres (except near 40 km in the Northern Hemisphere). The percent of the total variance accounted for by the second harmonic is maximum at a height of about 40 km in the tropics and at subpolar latitudes where, as mentioned, its amplitude is also largest.The phase of the first harmonic shows a marked transition from a winter/spring maximum below 30 km to a summer maximum at 30 km, changing rapidly to a maximum in winter in both hemispheres. The regions of minimum amplitude of the annual variation and the marked phase shifts with height both indicate the separation by levels of the dominant physical control mechanisms on the periodic changes of the ozone mixing ratio in the middle and upper stratosphere. Changes below 30 km respond primarily to dynamic influences in the lower stratosphere while above 30 km the periodic variations result mainly from photochemical processes. Above 40 km these variations are strongly temperature dependent.     

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.     

Global cloud water distribution derived from special sensor microwave imager/sounder and its comparison with GCM simulation

An SSM/I algorithm is developed for measuring cloud liquid water of raining and non-raining clouds. Cloud ice water corresponding to precipitation sized ice particles is obtained from the SSM/I scattering index which is calibrated against radar- derived ice water content. Smaller ice particles in upper level clouds are detected using scattering index based on SSM/T2 183±3 and 183±1 GHz channels which allow for removing water vapor emission. It is shown that SSM/I derived cloud liquid water distribution agrees well with GCM simulations, particularly in tropical latitudes.     

Global MGS TES data and Mars-GRAM validation

Mars Global Reference Atmospheric Model (Mars-GRAM 2001) is an engineering-level Mars atmosphere model widely used for many Mars mission applications. From 0 to 80 km, it is based on NASA Ames Mars General Circulation Model (MGCM), while above 80 km it is based on University of Michigan Mars Thermospheric General Circulation Model. Mars-GRAM 2001 and MGCM use surface topography from Mars Global Surveyor Mars Orbiter Laser Altimeter (MOLA). Validation studies are described comparing Mars-GRAM with a global summary data set of Mars Global Surveyor Thermal Emission Spectrometer (TES) data. TES averages and standard deviations were assembled from binned TES data which covered surface to ∼40 km, over more than a full Mars year (February 1999–June 2001, just before start of a Mars global dust storm). TES data were binned in 10° × 10° latitude–longitude bins (36 longitude bins, centered at 5°–355°, by 18 latitude bins, centered at −85° to +85°), and 12 seasonal bins (based on 30° increments of L_s angle). Bin averages and standard deviations were assembled at 23 data levels (temperature at 21 pressure levels, plus surface temperature and surface pressure). Two time-of-day bins were used: local time near 2 or 14 h. Two dust optical depth bins were used: infrared optical depth, either less than or greater than 0.25 (which corresponds to visible optical depth less than or greater than about 0.5). For interests in aerocapture and precision entry and landing, comparisons focused on atmospheric density. TES densities versus height were computed from TES temperature versus pressure, using assumptions of perfect gas law and hydrostatics. Mars-GRAM validation studies used density ratio (TES/Mars-GRAM) evaluated at data bin center points in space and time. Observed average TES/Mars-GRAM density ratios were generally 1 ± 0.05, except at high altitudes (15–30 km, depending on season) and high latitudes (>45°N), or at most altitudes in the southern hemisphere at L_s ∼ 90° and 180°. Compared to TES averages for a given latitude and season, TES data had average density standard deviation about the mean of ∼2.5% for all data, or ∼1–4%, depending on time of day and dust optical depth. Average standard deviation of TES/Mars-GRAM density ratio was 8.9% for local time 2 h and 7.1% for local time 14 h. Thus standard deviation of observed TES/Mars-GRAM density ratio, evaluated at matching positions and times, is about three times the standard deviation of TES data about the TES mean value at a given position and season.     

Introduction to Chinese Meridian Project-Phase II

The Chinese Meridian Project is a ground-based space environment monitoring network, which is constructed in two steps. The first step (Phase I) of the project consists of 15 observation stations located roughly along 120°E longitude and 30°N latitude. The second step (Phase II) of the project will additionally deploy 16 stations to better cover China's territory, and build a stereo monitoring capability to monitor the cause and effect of the space weather chain in the solar terrestrial system. Based on the existing two monitoring chains in Phase I, two more chains will be established along 100°E longitude and 40°N latitude, respectively, forming a double-cross network configuration. After the two-step construction, the whole project will run nearly 300 instruments deployed at 31 stations. Aside from standard instruments, quite a few innovative and powerful instruments will be developed, such as radioheliographs with a very wide frequency band, a 3-station incoherent scattering radar to make a 3D measurement of the ionosphere, and a helium lidar to measure atmosphere density up to an altitude of 1000 km.      

Mesospheric doppler wind measurements from Aura Microwave Limb Sounder (MLS)

This paper describes a microwave limb technique for measuring Doppler wind in the Earth’s mesosphere. The research algorithm has been applied to Aura Microwave Limb Sounder (MLS) 118.75 GHz measurements where the O₂ Zeeman lines are resolved by a digital autocorrelation spectrometer. A precision of ∼17 m/s for the line-of-sight (LOS) wind is achieved at 80–92 km, which corresponds to radiometric noise during 1/6 s integration time. The LOS winds from Aura MLS are mostly in the meridional direction at low- and mid-latitudes with vertical resolution of ∼8 km. This microwave Doppler technique has potential to obtain useful winds down to ∼40 km of the Earth’s atmosphere if measurements from other MLS frequencies (near H₂O, O₃, and CO lines) are used. Initial analyses show that the MLS winds from the 118.75 GHz measurements agree well with the TIDI (Thermosphere Ionosphere Mesosphere Energetics and Dynamics Doppler Interferometer) winds for the perturbations induced by a strong quasi 2-day wave (QTDW) in January 2005. Time series of MLS winds reveal many interesting climatological and planetary wave features, including the diurnal, semidiurnal tides, and the QTDW. Interactions between the tides and the QTDW are clearly evident, indicating possible large tidal structural changes after the QTDW events dissipate.     

Seismogenic ionospheric anomalies associated with the strong Indonesian earthquake occurred on 11 April 2012 (M = 8.5)

Ionospheric perturbations in possible association with a major earthquake (EQ) (M?=?8.5) which occurred in India-Oceania region are investigated by monitoring subionospheric propagation of VLF signals transmitted from the NWC transmitter (F?=?19.8?kHz), Australia to a receiving station at Varanasi (geographic lat. 25.3°N, long 82.99°E), India. The EQ occurred on 11 April 2012 at 08:38:35?h UT (magnitude?≈?8.5, depth?=?10?km, and lat.?=?2.3°N, long.?=?93.0°E). A significant increase of few days before the EQ has been observed by using the VLF nighttime amplitude fluctuation method (fixed frequency transmitter signal). The analysis of total electron contents (TEC) derived from the global positioning system (GPS) at three different stations namely, Hyderabad (latitude 17.38°N, longitude 78.48°E), Singapore (latitude 1.37°N, longitude 103.84°E) and Port Blair (latitude 11.62°N, longitude 92.72°E) due to this EQ has also been presented. Significant perturbation in TEC data (enhancements and depletion) is noted before and after the main shock of the EQ. The possible mechanisms behind these perturbations due to EQ have also been discussed.