Gravity wave mixing effects on the OH*-layer

Based on an advanced numerical model for excited hydroxyl (OH*) we simulate the effects of gravity waves (GWs) on the OH*-layer in the upper mesosphere. The OH* model takes into account (1) production by the reaction of atomic hydrogen (H) with ozone (O₃), (2) deactivation by atomic oxygen (O), molecular oxygen (O₂), and molecular nitrogen (N₂), (3) spontaneous emission, and (4) loss due to chemical reaction with O. This OH* model is part of a chemistry-transport model (CTM) which is driven by the high-resolution dynamics from the KMCM (Kühlungsborn Mechanistic general Circulation Model) which simulates mid-frequency GWs and their effects on the mean flow in the MLT explicitly. We find that the maximum number density and the height of the OH*-layer peak are strongly determined by the distribution of atomic oxygen and by the temperature. As a results, there are two ways how GWs influence the OH*-layer: (1) through the instantaneous modulation by O and T on short time scales (a few hours), and (2) through vertical mixing of O (days to weeks). The instantaneous variations of the OH*-layer peak altitude due to GWs amount to 5–10 km. Such variations would introduce significant biases in the GW parameters derived from airglow when assuming a constant pressure level of the emission height. Performing a sensitivity experiment we find that on average, the vertical mixing by GWs moves the OH*-layer down by ~2 to 7 km and increases its number density by more than 50%. This effect is strongest at middle and high latitudes during winter where secondary GWs generated in the stratopause region account for large GW amplitudes.     

Validation of GOCE densities and evaluation of thermosphere models

Atmospheric densities from ESA’s GOCE satellite at a mean altitude of 270 km are validated by comparison with predictions from the near real time model HASDM along the GOCE orbit in the time frame 1 November 2009 through 31 May 2012. Except for a scale factor of 1.29, which is due to different aerodynamic models being used in HASDM and GOCE, the agreement is at the 3% (standard deviation) level when comparing daily averages. The models NRLMSISE-00, JB2008 and DTM2012 are compared with the GOCE data. They match at the 10% level, but significant latitude-dependent errors as well as errors with semiannual periodicity are detected. Using the 0.1 Hz sampled data leads to much larger differences locally, and this dataset can be used presently to analyze variations down to scales as small as 150 km.     

利用静电双探针对辉光放电等离子体诊断分析

阐述了利用Langmuir双探针诊断等离子体参数的基本原理,对圆柱腔中的空气辉光放电等离子体进行了双探针诊断,通过实验测得伏安特性曲线,计算得出等离子体密度和电子温度。分析了密度参数变化与放电电流大小之间的关系。最后将诊断结果与准光腔系统进行比较,二者虽略有差别,但差别不大,且变化趋势一致。验证了采用本双探针系统诊断辉光放电等离子体密度的可靠性。     

Secular variations of OH nightglow emission and of the OH intensity-weighted temperature induced by gravity–wave forcing in the MLT region

Secular variations of OH airglow (8,3) band and of the intensity-weighted temperature induced by a gravity wave packet at three latitudes in the northern hemisphere in the Mesosphere/Lower Thermosphere region are investigated using a 2-D, time-dependent, fully nonlinear OH chemistry–dynamics model. A net-cycled average survives due to wave transience and dissipation. The integrated OH (8,3) volume emission rates show a 20 percent increase induced by the wave packet. The corresponding intensity-weighted temperature shows an initial increase up to 2 percent then a gradual decrease for the remainder of the simulation time by the same wave packet. These secular variations could be mistaken as long-period or short-period waves in the airglow observations. Therefore, care must be taken when analyzing the data from observations.     

Evolution of Dst index,cosmic rays and global temperature during solar cycles 20–23

We have studied conditions in interplanetary space, which can have an influence on galactic cosmic ray (CR) and climate change. In this connection the solar wind and interplanetary magnetic field parameters and cosmic ray variations have been compared with geomagnetic activity represented by the equatorial Dst index from the beginning 1965 to the end of 2012. Dst index is commonly used as the solar wind–magnetosphere–ionosphere interaction characteristic. The important drivers in interplanetary medium which have effect on cosmic rays as CMEs (coronal mass ejections) and CIRs (corotating interaction regions) undergo very strong changes during their propagation to the Earth. Because of this CMEs, coronal holes and the solar spot numbers (SSN) do not adequately reflect peculiarities concerned with the solar wind arrival to 1 AU. Therefore, the geomagnetic indices have some inestimable advantage as continuous series other the irregular solar wind measurements. We have compared the yearly average variations of Dst index and the solar wind parameters with cosmic ray data from Moscow, Climax, and Haleakala neutron monitors during the solar cycles 20–23. The descending phases of these solar cycles (CSs) had the long-lasting solar wind high speed streams occurred frequently and were the primary contributors to the recurrent Dst variations. They also had effects on cosmic rays variations. We show that long-term Dst variations in these solar cycles were correlated with the cosmic ray count rate and can be used for study of CR variations. Global temperature variations in connection with evolution of Dst index and CR variations is discussed.     

Comparison of GPS-TEC measurements with NeQuick2 and IRI model predictions in the low latitude East African region during varying solar activity period (1998 and 2008–2015)

This paper examines the performances of NeQuick2, the latest available IRI-2016, IRI-2012 and IRI-2007 models in describing the monthly and seasonal mean total electron content (TEC) over the East African region. This is to gain insight into the success of the various model types and versions at characterizing the ionosphere within the equatorial ionization anomaly. TEC derived from five Global Positioning System (GPS) receivers installed at Addis Ababa (ADD, 5.33°N, 111.99°E Geog.), Asab (ASAB, 8.67°N, 116.44°E Geog.), Ambo (ABOO, 5.43°N, 111.05°E Geog.), Nairobi (RCMN, ?4.48°N, 108.46°E Geog.) and Nazret (NAZR, 4.78°N, 112.43°E Geog.), are compared with the corresponding values computed using those models during varying solar activity period (1998 and 2008–2015). We found that different models describe the equatorial and anomaly region ionosphere best depending on solar cycle, season and geomagnetic activity levels. Our results show that IRI-2016 is the best model (compared to others in terms of discrepancy range) in estimating the monthly mean GPS-TEC at NAZR, ADD and RCMN stations except at ADD during 2008 and 2012. It is also found that IRI-2012 is the best model in estimating the monthly mean TEC at ABOO station in 2014. IRI show better agreement with observations during June solstice for all the years studied at ADD except in 2012 where NeQuick2 better performs. At NAZR, NeQuick2 better performs in estimating seasonal mean GPS-TEC during 2011, while IRI models are best during 2008–2009. Both NeQuick2 and IRI models underestimate measured TEC for all the seasons at ADD in 2010 but overestimate at NAZR in 2009 and RCMN in 2008. The periodic variations of experimental and modeled TEC have been compared with solar and geomagnetic indices at ABOO and ASAB in 2014 and results indicate that the F10.7 and sunspot number as indices of solar activity seriously affects the TEC variations with periods of 16–32?days followed by the geomagnetic activity on shorter timescales (roughly periods of less than 16?days). In this case, NeQuick2 derived TEC shows better agreement with a long term period variations of GPS-TEC, while IRI-2016 and IRI-2007 show better agreement with observations during short term periodic variations. This indicates that the dependence of NeQuick2 derived TEC on F10.7 is seasonal. Hence, we suggest that representation of geomagnetic activity indices is required for better performance over the low latitude region.     

Comparison of time-periodic variations in temperature and wind from meteorological rockets and satellites

Although the Meteorological Rocket Network operated by or in cooperation with the United States has decreased from fourteen to nine stations in the past five years, there have been many observations accumulated in the ten years since CIRA 1972 was prepared with data up to 1969. The mean, annual and semiannual variations of temperature and wind are presented and special attention is directed to the polar semiannual wave. The results are compared with the Oxford SCR-PMR five-year data set, the CDC-SCR seven-year data, and CIRA 1972 with respect to both temperature and zonal winds, as far as presently available. The agreement among the data sets is very good.     

The influence of radiative cooling and turbulence on the heat budget of the thermosphere

Numerical models of the thermal budget of the Earth's upper atmosphere in the height range of 90–500km are developed. The main sources and sinks of energy including infra-red radiative cooling by vibrational-rotational bands of NO, CO₂, OH and O₃ as well as heating and cooling arising from dissipation of turbulent energy and eddy heat transport are taken into account. The calculated temperature and density height profiles are in good agreement with the respective profiles from CIRA 72 and Jacchia 1977 models. It is shown for the models considered that IR-radiative cooling by CO₂ and NO in the 15μ and 5.3μ bands, not eddy turbulence provides the major loss of heat from 90 to 180km.     

Description of the mean behaviour of ionospheric plasma temperatures

A status report on the empirical modeling of ionospheric electron and ion temperatures is given with special emphasis on the models used in the International Reference Ionosphere (IRI).Electron temperature models have now reached a state where reliable prediction of the mean altitudinal, latitudinal and diurnal variations is possible. These models are largely based on satellite measurements, but comparisons with incoherent scatter radar measurements have shown excellent agreement. Variations with season and magnetic and solar activity seem to be small and are not yet included consistently in these models.Similar to the electron temperature, the ion temperature shows the largest variations with altitude, latitude and local time. But due to the larger mass, these variations are smoother and more steady in the case of the ions and therefore easier to model. Nevertheless, very few ion temperature models exist. The IRI model takes advantage of the observed concurrence of the ion temperature with the neutral temperature at low altitudes and with the electron temperature at high altitudes.     

Comparison of model and experimental ionospheric parameters at high latitudes

Modern use and study of the auroral region needs to attract a wider class of models for describing conditions of radio wave propagation in the ionosphere. In this paper the possibilities of the International Reference Ionosphere model, well-proven and widespread in the mid-latitudes, are investigated in the high latitude zone. Model and measured values of the critical frequency foF2 for two mid-latitude stations (Juliusruh and Goosebay) and four high-latitude ones (Loparskaya, Sodankyla, Sondrestrom, Thule) are compared. Deviations of medians, variations from day to day and solar activity trends seemed to be similar for both areas. This similarity is irrespective of the RZ12 index. Special attention is paid to the TEC parameter and its determination using 6 versions of models, a new version of the model IRI2010 (IRI-Plas) among them. It is shown that the IRI-Plas model significantly improves the definition of TEC in contrast to the versions of IRI2007 and the new model NeQuick. The use of the median of the experimental equivalent slab thickness, together with the current values of the TEC, increases by a factor of two the agreement between calculated and measured foF2 values as compared with the variations from day to day. This allows foF2 to be defined in near-real time.     

利用Aura卫星资料计算全球中层大气背景风场

利用2005年第二代地球观测卫星系统(EOS)的Aura卫星上MLS观测的压强、温度、密度等数据,推算出全球中层大气的平均背景风场,分析了中层大气风场随时间和高度变化的特点.与武汉流星雷达及澳大利亚Adelaide台站观测的比较结果显示,用Aura数据推算出来的风场与实际观测比较符合.与HWM-93模式的比较显示,Aura风场随时间和空间变化的总体趋势与HWM-93基本吻合.特别是在80km以下的高度范围内Aura数据与HWM-93数据符合得比较好;在80 km以上的高度,Aura所算得的风值与HWM-93风值的差别逐渐增大,Aura风值普遍比HWM-93的要大.      

Ozone reference models for CIRA

The ozone reference model which is to be incorporated in the COSPAR International Reference Atmosphere (CIRA) is described and compared with other measurements of the Earth's ozone distribution. Ozone vertical structure models from approximately 25 to 90 km are provided combining data from five recent satellite experiments (Nimbus-7 LIMS, Nimbus-7 SBUV, AE-2 SAGE, Solar Mesosphere Explorer (SME) UVS, and SME IR). The results include the latest improvements in the SBUV algorithms using the most recent estimates of ozone cross sections. Also, the latest refinements in SME algorithms are incorporated. These algorithm improvements have improved agreement between the satellite data sets. Standard deviations are provided of monthly zonal means, and an estimate of the interannual variability is given. The models based on satallite data compare well with the Krueger and Minzner mid-latitude model incorporated into the U. S. Standard Atmosphere which is based on rocket and balloon measurements. Other comparisons are shown with Umkehr and more recent balloon data. Models are also provided of total columnar ozone reflecting recent improved estimates of ozone cross section. Information is provided on semiannual and annual variations. Other systematic variations including estimates of diurnal variations in the mesosphere will be included in the CIRA document.     

Chemical heating rates derived from SCIAMACHY vibrationally excited OH limb emission spectra

The SCIAMACHY instrument on board Envisat is able to measure nearly all vibrational transitions of mesospheric hydroxyl – from the ultraviolet to the near infrared spectral region.

In this paper, we analyze SCIAMACHY limb emission data in the 1000–1750 nm spectral region by means of a new vibrational non-LTE model of OH. Several hydroxyl hotbands are identified. Vibrational non-LTE model calculations applying different collisional relaxation models are able to reproduce the measured spectra. Best agreement between model calculations and measured spectra is obtained, if a combination of multiquantum and step ladder single-quantum relaxation model is applied. Emissions from the OH(v = 9) vibrational state are used to derive chemical heating rates from the SCIAMACHY spectra. Instantaneous heating rates are in the order of 10 K/day.     

Connection between winter mesopause region temperatures and diurnal LF reflection height variations measured at Collm

Scale height, H, estimates are calculated from the decrease/increase of ionospheric virtual reflection heights of low-frequency (LF) radio waves at oblique incidence in suitably defined morning intervals around sunrise during winter months. The day-to-day variations of H qualitatively agree with daily mean temperature variations around 90 km from meteor radar measurements. Since mesospheric long-period temperature variations are generally accepted to be the signature of atmospheric planetary waves, this shows that LF reflection height measurements can be used for monitoring the dynamics of the upper middle atmosphere. The long-term variations of monthly mean H estimates have also been analysed. There is no significant trend, which is in agreement with other measurements of mesopause region temperature trends.     

Solar eclipse induced variations in mesospheric ozone concentrations

A model study of the solar eclipse induced variations in mesospheric ozone concentrations was undertaken. This study includes, in addition to the Chapman reactions, the chemistry involving hydrogen species (H, OH and HO₂) which are found to be important in the destruction of odd oxygen in this altitude region. Coupled time dependent continuity equations are solved for the eclipse duration. The results from the present study are compared with earlier theoretical model and the experimental observations during the 16 February 1980 solar eclipse as well as the results obtained during earlier solar eclipses.     

低轨航天器弹道系数估算及热层大气模型误差分析

利用低轨（LEO）航天器在轨期间两行轨道根数（TLEs）数据,结合经验大气密度模型NRLMSISE00,反演计算得到其在轨期间的弹道系数B’,以31年B’的平均值代替弹道系数真值,分别通过标准球形目标卫星对比以及物理参数基本相同的非球形目标卫星对比,对弹道系数真值进行了检验;利用不同外形目标卫星弹道系数在不同太阳活动周内的变化规律,结合太阳和地磁活动变化,估计经验大气密度模型的误差分布. 结果表明,利用反演弹道系数31年的平均值来代替真值,其在理论值的正常误差范围内;大气密度模型误差在210～526km高度范围内存在相同的变化趋势,且模型误差随高度增加而增大;在短周期内B’变化与太阳活动指数F_10.7存在反相关性;密度模型不能有效模拟2008年出现的大气密度异常低. 以上结果表明,经验大气密度模型结果需要修正,尤其是在太阳活动峰年和谷年,此外,磁暴期间模型误差的修正对卫星定轨和轨道预报等也具有重要意义.      

Middle atmosphere models and comparison with shuttle reentry density data

Several middle atmosphere models will be reviewed, including a new set of models produced by Groves in 1985. The latter models are based on rocket and rawinsonde in situ measurements and satellite remote sounding temperature data. The models are compared with measurements made with instruments on board U. S. Shuttles during their reentry. Very useful atmospheric density data have been obtained in the altitude region from 50 to 80 km. The measurements are unique in that they are made by a vehicle travelling almost horizontally through the atmosphere at a velocity of 6 to 7 km/sec. This results in measurements along a path of approximately 8,000 km in a time interval of about 20 minutes. The results show some unique features.     

Results from the LIMS experiment for the PMP-1 winter 1978/79

The Limb Infrared Monitor of the Stratosphere (LIMS) experiment utilized a 6-channel limb scanning infrared radiometer on the Nimbus 7 spacecraft, for the purpose of determining global distributions of temperature and trace constituents in the middle atmosphere (15–65 km) with high vertical resolution.LIMS observations provide a detailed picture of atmospheric variability during the disturbed winter of 1978/79. During January and February, three large disturbances caused simultaneous temperature changes through most of the global middle atmosphere. The amplitudes and phases of the planetary temperature and height waves show large variations during this time. Their values and derived quantities are larger than those resulting from data with lower vertical resolution. The Eliassen-Palm fluxes calculated for late January are in qualitative agreement with the observed acceleration of the zonal winds.     

Mesospheric OH temperatures: Simultaneous ground-based and SABER OH measurements over Millstone Hill

We present rotational temperature measurements of the mesospheric OH emission layer using a meridional imaging spectrograph at Millstone Hill (42.6°N, 72.5°W). The system is equipped with a state-of-the-art bare-CCD detector and can yield simultaneous quasi-meridional images of the mesospheric OH and O₂ intensity and temperature fields at 87 and 94 km altitude during the course of each night. A cross-validation study of the rotational OH temperature measurements obtained on 61 nights during the autumnal months of 2005–2007 was undertaken with near-simultaneous kinetic temperature measurements made by the SABER instrument aboard the NASA TIMED satellite during overpasses of Millstone Hill. Excellent agreement was obtained between the two datasets with the small differences being attributable to differences in the spatial and temporal averaging inherent between the two datasets.     

Measurement of mesopause temperature from hydroxyl nightglow at Kolhapur (16.8°N, 74.2°E), India

The paper reports the nightglow observations of hydroxyl (8–3), (7–2) and (6–2) Meinel band carried out at a low latitude station Kolhapur (16.8°N, 74.2°E, dip latitude 10.6°N), India during November 2002 to May 2005 with the objective of investigating mesopause dynamics based on derived OH rotational temperature. Overall, 132 nights of quality data were collected using filter-tilting photometer and an all-sky scanning photometer. The mean mesopause temperature observed at Kolhapur is 195 ± 11, 196 ± 9 and 195 ± 7 K from OH (8–3), (7–2) and (6–2) band emissions, respectively, using transition probabilities given by Langhoff et al. [Langhoff, S.R., Werner, H.J., Rosmus, P. Theoretical transition probabilities for the OH Meinel system. Journal of Molecular Spectroscopy 118, 507–529, 1986]. Small wave-like variations (periodicities ∼ few hours) existing over long period variations in derived temperatures are also present. A steady decrease of emission intensities from evening to dawn hours has been observed in approximately 59% of nights. No significant change of nightly mean temperatures has been noted. Furthermore, about 62% of observed nightly mean temperatures lie within one error bar of MSISE-90 model predictions.