GGTm-Ts: A global grid model of weighted mean temperature (Tm) based on surface temperature (Ts) with two modes

With the development of Global Navigation Satellite System (GNSS), the detection of precipitable water vapor (PWV) using the GNSS atmospheric sounding technique becomes a research interest in GNSS meteorology. In the conversion of zenith tropospheric delay (ZTD) to PWV, the weighted mean temperature (Tm) plays a crucial role. Generally, the Tm estimated by the linear regression models based on surface temperature (Ts) cannot meet the requirement for global use, and the accuracy of Tm derived from the empirical models is limited. In this study, a new Tm model, named GGTm-Ts model, was developed using the global geodetic observing system (GGOS) atmosphere Tm data and European Centre for Medium-Range Weather Forecasts (ECMWF) data from 2011 to 2015. Resting upon a global 2.5°*2° grid of coefficients of Tm-Ts linear function, the new model can provide Tm at any site in two modes, one for the case with measured Ts provided, i.e., the accurate mode, the other for the case that Ts provided by a subroutine, i.e., the normal mode. The performance of GGTm-Ts model was assessed against the Bevis formula, GPT2w and GPT2wh model using different data sources in 2016-the GGOS atmosphere and radiosonde data. The results show that the GGTm-Ts model in accurate mode achieves best performance with an improvement of 46.9 %/15.3 %, 37.8 %/19.5 % and 34.4 %/14.2 % over other three models in the GGOS atmosphere/radiosonde comparison. For the normal mode, the GGTm-Ts model outperforms the GPT2w model and achieves equivalence results with the GPT2wh model. Moreover, the impact of Tm on GNSS-PWV was analyzed to validate the performance of the GGTm-Ts model.     

高纬度地区GPT2w模型的适应性分析

GPT2w模型是现有精度最高的天顶对流层模型,但是应用在高纬度地区时存在较大误差.为更好地保障卫星导航定位系统在高纬度地区的高精度应用,评定了GPT2w模型在高纬度地区的精度,获取天顶对流层湿延迟、干延迟和总延迟,探讨了GPT2w模型改正对精密单点定位的影响.试验结果表明: GPT2w模型在高纬度地区的精度为厘米级,优于其在中低纬度地区的精度;南北极地区天顶对流层呈现明显季节变化特征和区域一致性特征,夏季天顶对流层总延迟高于冬季,北极地区天顶对流层湿延迟明显高于南极地区,北极地区天顶对流层随季节的变化幅度大于南极地区.PPP试验结果表明,GPT2w模型能够有效改善定位精度,适应高纬度地区的高精度定位需求.      

High temporal resolution global PWV dataset of 2005–2016 by using a neural network approach to determine the mean temperature of the atmosphere

Atmospheric water vapour plays an important role in phenomena related to the global hydrologic cycle and climate change. However, the rapid temporal–spatial variation in global tropospheric water vapour has not been well investigated due to a lack of long-term, high-temporal-resolution precipitable water vapour (PWV). Accordingly, this study generates an hourly PWV dataset for 272 ground-based International Global Navigation Satellite System (GNSS) Service (IGS) stations over the period of 2005–2016 using the zenith troposphere delay (ZTD) derived from global-scale GNSS observation. The root mean square (RMS) of the hourly ZTD obtained from the IGS tropospheric product is approximately 4 mm. A fifth-generation reanalysis dataset of the European Centre for Medium-range Weather Forecasting (ECMWF ERA5) is used to obtain hourly surface temperature (T) and pressure (P), which are first validated with GNSS synoptic station data and radiosonde data, respectively. Then, T and P are used to calculate the water vapour-weighted atmospheric mean temperature (Tm) and zenith hydrostatic delay (ZHD), respectively. T and P at the GNSS stations are obtained via an interpolation in the horizontal and vertical directions using the grid-based ERA5 reanalysis dataset. Here, Tm is calculated using a neural network model, whereas ZHD is obtained using an empirical Saastamoinen model. The RMS values of T and P at the collocated 693 radiosonde stations are 1.6 K and 3.1 hPa, respectively. Therefore, the theoretical error of PWV caused by the errors in ZTD, T and P is on the order of approximately 2.1 mm. A practical comparison experiment is performed using 97 collocated radiosonde stations and 23 GNSS stations equipped with meteorological sensors. The RMS and bias of the hourly PWV dataset are 2.87/?0.16 and 2.45/0.55 mm, respectively, when compared with radiosonde and GNSS stations equipped with meteorological sensors. Additionally, preliminary analysis of the hourly PWV dataset during the EI Niño event of 2014–2016 further indicates the capability of monitoring the daily changes in atmospheric water vapour. This finding is interesting and significant for further climate research.     

Lower ionosphere model of the polar regions during summer noontime

Measurements of various parameters in both Arctic and Antarctic mesosphere are available in the literature. These have been first examined to identify the similarities and differences between them in the two polar regions. Then an attempt has been made to reproduce them using an ion-chemical scheme in which both positive and negative ion processes have been included. Temperature and nitric oxide density are found to play a crucial role in the formation of these features. Finally a theoretical model of ion composition for Antarctic summer noontime condition is given.     

A global weighted mean temperature model based on empirical orthogonal function analysis

A global empirical orthogonal function (EOF) model of the tropospheric weighted mean temperature called GEOFM_Tm was developed using high-precision Global Geodetic Observing System (GGOS) Atmosphere T_m data during the years 2008–2014. Due to the quick convergence of EOF decomposition, it is possible to use the first four EOF series, which consists base functions U_k and associated coefficients P_k, to represent 99.99% of the overall variance of the original data sets and its spatial-temporal variations. Results show that U₁ displays a prominent latitude distribution profile with positive peaks located at low latitude region. U₂ manifests an asymmetric pattern that positive values occurred over 30° in the Northern Hemisphere, and negative values were observed at other regions. U₃ and U₄ displayed significant anomalies in Tibet and North America, respectively. Annual variation is the major component of the first and second associated coefficients P₁ and P₂, whereas P₃ and P₄ mainly reflects both annual and semi-annual variation components. Furthermore, the performance of constructed GEOFM_Tm was validated by comparison with GTm_III and GTm_N with different kinds of data including GGOS Atmosphere T_m data in 2015 and radiosonde data from Integrated Global Radiosonde Archive (IGRA) in 2014. Generally speaking, GEOFM_Tm can achieve the same accuracy and reliability as GTm_III and GTm_N models in a global scale, even has improved in the Antarctic and Greenland regions. The MAE and RMS of GEOFM_Tm tend to be 2.49?K and 3.14?K with respect to GGOS T_m data, respectively; and 3.38?K and 4.23?K with respect to IGRA sounding data, respectively. In addition, those three models have higher precision at low latitude than middle and high latitude regions. The magnitude of T_m remains at the range of 220–300?K, presented a high correlation with geographic latitude. In the Northern Hemisphere, there was a significant enhancement at high latitude region reaching 270?K during summer. GEOFM_Tm is capable to represent the spatiotemporal variations of T_m, with the high accuracy and reliability in a global scale, therefore, will be of great significance to the real-time GNSS water vapor inversion and climate studies.     

两极地区局部点质量模型构建方法

针对两极地区位系数模型计算扰动引力时的奇异问题和点质量模型方程解不稳定问题,分别引入扰动引力去奇异算法和基于换极法的两极地区局部点质量模型构造方法.数值实验结果表明,去奇异算法能有效解决靠近极点处扰动引力计算过程中出现的奇异性问题,基于换极法的点质量模型在两极地区方程结构稳定,精度与位系数模型相当.本文方法可为极区飞行器扰动引力快速赋值提供新的思路.      

Radiative and dynamical impacts of Arctic and Antarctic ozone holes: General circulation model experiments

Radiative and dynamical impacts of Arctic and Antarctic ozone holes on the general circulation are investigated with the aid of a general circulation model developed at Kyushu University. The model includes a simplified ozone photochemistry interactively coupled with radiation and dynamics. Resultant temperature structure consisting of a cooling in the polar lower stratosphere and a warming in the polar upper stratosphere brings about the intensification of the polar night jet. The cooling is caused by the decrease of solar ultraviolet heating due to the ozone depletion, while the warming is caused by adiabatic heating due to the enhancement of downward motion.     

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

利用2012-2014年1月和7月AIRS(Atmospheric Infrared Sounder)第79通道的观测数据,分析了平流层重力波活动强弱的全球分布以及重力波发生频率的全球分布;分析了重力波活动随纬度和经度的变化特征,给出了重力波活动在全球范围内的热点区域及其活动强度;对比了白天与夜间的重力波活动强度及发生...     

Satellite altimeter measurements of the geoid in sea ice zones

The Seasat radar altimeter provided surface height measurements to a precision better than 10cm over the open ocean, The data have been used to produce maps of the ocean geoid which reveal details of sub-surface topography such as sea mounts, ocean trenches and mid-ocean ridges /1/. In areas of the ocean covered by sea ice, however, the quasi-specular ice returns which occurred were incorrectly handled by the on-board processor. This resulted in a significant decrease in the precision of the surface height estimates. Consequently, researchers have generally eliminated data from regions where sea ice is suspected to have been present, including large areas of the Antarctic ocean. We have developed a technique for significantly improving the height measurements over such areas permitting the mapping of the geoid in such regions. The short wavelength RMS deviation of elevation measurements from collinear passes over such areas has been reduced from 1.032m to 0.632m. The application of the technique to ERS-1 altimeter data will be particularly important, since coverage of a substantial area of the Arctic ocean (up to 82° latitude) will be possible for the first time     

Influence of lunar topography on simulated surface temperature

The surface temperature of the Moon is one of the essential parameters for the lunar exploration, especially to evaluate the Moon thermophysical features. The distribution of the temperature is heavily influenced by the Moon topography, which, however, is rarely studied in the state-of-art surface temperature models. Therefore, this paper takes the Moon topography into account to improve the surface temperature model, Racca model. The main parameters, such as slopes along the longitude and latitude directions, are estimated with the topography data from Chang’E-1 satellite and the Horn algorithm. Then the effective solar illumination model is then constructed with the slopes and the relative position to the subsolar point. Finally, the temperature distribution over the Moon surface is obtained with the effective illumination model and the improved Racca model. The results indicate that the distribution of the temperature is very sensitive to the fluctuation of the Moon surface. The change of the surface temperature is up to 150 K in some places compared to the result without considering the topography. In addition, the variation of the surface temperature increases with the distance from the subsolar point and the elevation, along both latitude and longitude directions. Furthermore, the simulated surface temperature coincides well with the brightness temperature in 37 GHz observed by the microwave sounder onboard Chang’E-2 satellite. The corresponded emissivity map not only eliminates the influence of the topography, but also hints the inherent properties of the lunar regolith just below the surface. Last but not the least, the distribution of the permanently shadowed regions (PSRs) in the lunar pole area is also evaluated with the simulated surface temperature result.     

Stratosphere-mesosphere coupling during stratospheric sudden warming events

In this review article we summarize recent results in the coupling of the stratosphere–mesosphere during stratospheric sudden warming (SSW) events. We focus on the role of planetary and gravity waves in driving the middle atmosphere circulation and illustrate the stratosphere–mesosphere coupling during undisturbed wintertime circulation, during an SSW event, and after an SSW event during the formation of an elevated stratopause using simulations of past Arctic and Antarctic winters from the Specified Dynamics version of the Whole Atmosphere Community Climate Model (SD-WACCM). We illustrate the transition of the polar stratopause from being a gravity wave driven phenomena to a planetary wave driven phenomena during SSW events and its subsequent reestablishment and control by gravity waves. We also examine the synoptic structure of the stratosphere, mesosphere, and lower thermosphere using SD-WACCM data fields that show the structure of the vortex during specific dynamical events in both hemispheres. We illustrate the longitudinal asymmetry in the thermal structure in the stratosphere and mesosphere driven by differences in circulation over the polar cap regions during an SSW event. We complement this analysis of the middle atmosphere circulation with a classification of both the Arctic and Antarctic winters since 1979 into major, minor, elevated stratopause or quiet winters based on the level of disturbance using the Modern Era-Retrospective Analysis for Research and Applications (MERRA) reanalysis data. From the MERRA data we find that the combined occurrences of both major and minor warmings in the Arctic have remained constant over the past three decades while we find a minor increase in their occurrences in the Antarctic.     

太阳高能粒子事件上升时间统计研究

选取1997-2006年共66个较大的缓变型太阳高能粒子(SEP)事件, 分析了不同条件下太阳高能粒子通量廓线上升时间与源区日面经向分布之间的相关关系, 研究了日冕物质抛射(CME)和耀斑在SEP上升阶段的作用特点.统计结果表明,大SEP事件的源区主要分布在太阳西半球, 特别是磁足点东西两侧45°范围内; 在高速太阳风条件下, 低能通道的通量上升时间与日面相对经度有较好的相关性,即离磁足点越远, 上升时间越长,而高能通道相关性则不明显; 全晕状CME产生的SEP事件对应的上升时间与源区位置没有明显的相关性, 而部分晕状CME伴随的SEP事件则与二次拟合曲线符合很好.分析表明,在缓变型SEP事件的通量上升阶段, 耀斑加速过程起着重要作用,这在部分晕状CME伴随的SEP事件中尤为显著.      

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.     

基于陆态网的区域电离层TEC空间变化特性分析

为研究中国陆态网区域电离层TEC在空间小尺度、高分辨率情况下的变化特性及适用精度范围,利用陆态网260个GNSS连续运行观测站数据,解算并生成2016-2017年731天陆态网区域电离层RIM格网,并进行精度验证.在同一RIM格网中,分别在经度和纬度方向上对间隔不同经纬度的TEC格网点作差分析.结果表明:陆态网区域内经度方向上TEC最大变化率和平均变化率分别为0.30TECU·（°）-1和0.11TECU·（°）-1;经度间隔1°时,TEC差值小于2TECU,且随着经度间隔的增大,其TEC差值也随之增大,并表现出一定的半年和周年变化规律;纬度方向上TEC最大变化率和平均变化率分别为1.7TECU·（°）-1和0.46TECU·（°）-1;陆态网区域内电离层TEC随纬度减小而增大,纬度间隔1°时,99.4%的TEC差值小于4TECU,且随着纬度间隔的增大,其TEC差值也随之增大,并表现出一定的半年和周年变化规律;间隔相同情况下,纬度方向上TEC的变化比经度方向上大.      

TIMED/SABER与AURA/MLS临近空间探测温度数据比较

利用AURA/MLS数据（V4.2）和TIMED/SABER数据（V2.0）对20~92km高度的大气温度进行比较分析,计算AURA/MLS数据与TIMED/SABER数据的温度绝对偏差,并分析平均温度偏差在不同季节中随经度、纬度和高度的变化特征.结果表明:20~80km高度的平均温度偏差在±6K以内,相对偏差在3%以内;80~90km高度平均温度偏差减小至-10K以下,相对偏差在9%以内.中低纬度地区平均温度偏差廓线的变化趋势一致,从20km高度的-3K左右的负偏差逐渐增加,在45~50km高度的平流层顶处有较明显的3K左右的正偏差峰值.平均温度偏差随纬度变化明显,随经度变化很小.研究结果可为卫星数据的应用提供参考依据.      

Evaluation of the improvement of IRI-2016 over IRI-2012 at the India low-latitude region during the ascending phase of cycle 24

This paper investigated the performance of the latest International Reference Ionosphere model (IRI-2016) over that of IRI-2012 in predicting total electron content (TEC) at three different stations in the Indian region. The data used were Global Positioning System (GPS) data collected during the ascending phase of solar cycle 24 over three low-latitude stations in India namely; Bangalore (13.02°N Geographic latitude, 77.57°E Geographic longitude), Hyderabad (17.25°N Geographic latitude, 78.30°E Geographic longitude) and Surat (21.16°N Geographic latitude, 72.78°E Geographic longitude). Monthly, the seasonal and annual variability of GPS-TEC have been compared with those derived from International Reference Ionosphere IRI-2016 and IRI-2012 with two different options of topside electron density: NeQuick and IRI01-corr. It is observed that both versions of IRI (i.e., IRI-2012 and IRI-2016) predict the GPS-TEC with some deviations, the latest version of IRI (IRI-2016) predicted the TEC similar to those predicted by IRI-2012 for all the seasons at all stations except for morning hours (0500 LT to 1000?LT). This shows that the effect of the updated version is seen only during morning hours and also that there is no change in TEC values by IRI-2016 from those predicted by IRI-2012 for the rest of the time of the day in the Indian low latitude region. The semiannual variations in the daytime maximum values of TEC are clearly observed from both GPS and model-derived TEC values with two peaks around March-April and September-October months of each year. Further, the Correlation of TEC derived by IRI-2016 and IRI-2012 with EUV and F10.7 shows similar results. This shows that the solar input to the IRI-2016 is similar to IRI 2012. There is no significant difference observed in TEC, bottom-side thickness (B0) and shape (B1) parameter predictions by both the versions of the IRI model. However, a clear improvement is visible in hmF2 and NmF2 predictions by IRI-2016 to that by IRI-2012. The SHU-2015 option of the IRI-2016 gives a better prediction of NmF2 for all the months at low latitude station Ahmedabad compare to AMTB 2013.     

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.     

Response of the extratropical middle atmosphere to the September 2002 major stratospheric sudden warming

The effects of a major stratospheric sudden warming (SSW) at extratropical latitudes have been investigated with wind and temperature observations over a Brazilian station, Cachoeira Paulista (22.7°S, 45°W) during September–October 2002. In response to the warming at polar latitudes a corresponding cooling at tropical and extratropical latitudes is prominent in the stratosphere. A conspicuous signature of latitudinal propagation of a planetary wave of zonal wavenumbers 1 and 2 from polar to low latitude has been observed during the warming period. The polar vortex which split into two parts of different size is found to travel considerably low latitude. Significant air mass mixing between low and high latitudes is caused by planetary wave breaking. The meridional wind exhibits oscillations of period 2–4 days during the warming period in the stratosphere. No wave feature is evident in the mesosphere during the warming period, although a 12–14 day periodicity is observed after 2 weeks of the warming event, indicating close resemblance to the results of other simultaneous investigations carried out from high latitude Antarctic stations. Convective activity over the present extratropical station diminishes remarkably during the warming period. This behavior is possibly due to destabilization and shift of equatorial convective active regions towards the opposite hemisphere in response to changes in the mean meridional circulation in concert with the SSW.     

一种火星大气密度三维解析模型

根据火星大气数据库,建立了以高度和经纬度为输入的火星三维解析大气密度模型;在竖直方向上,采用分层指数模型;在水平方向上,将指数模型中的参考密度和参考高度视为经纬度的多项式函数,并通过最小二乘曲面拟合得到了多项式系数。相比传统的指数模型,三维指数模型不仅反映了火星大气密度随高度的变化规律,而且能够体现出大气密度在水平方向上的变化。与大气数据库相比,三维密度模型具有解析形式,计算速度快,仿真过程中只需存储简单的多项式系数即可实现密度的实时计算和调用,适用于地面实时仿真。     

A statistical analysis of low frequency geomagnetic field pulsations at two Antarctic geomagnetic observatories in the polar cap region

The aim of this study is to investigate the characteristics of low frequency (∼0.5–5 mHz) geomagnetic field fluctuations as recorded at two Antarctic stations within the polar cap: the Italian observatory Mario Zucchelli Station (TNB) and the French–Italian observatory Dome C (DMC) in order to investigate the spatial extension and propagation characteristics of the phenomena observed at very high latitude. The stations have approximately the same geographic latitude, but a very different corrected geomagnetic latitude, being DMC close to the geomagnetic pole and TNB closer to the auroral oval.