The D-region ion composition and electron density response to strong solar proton events (model simulations)

Fluxes of energetic solar protons penetrate deep into the Earth’s polar cap middle atmosphere. Interacting with molecules of the air they cause additional dissociation and ionization, and the formed NO_x, OH_y and ions enter chemical and ion-molecular reactions. Induced changes of the ionospheric D-layer are modeled by a 1D model of lower ionosphere with chemistry, using neutral species concentrations calculated by a 1D photochemical time-dependent model. Changes of the electron and ion densities, and the most important ionospheric parameters are calculated after SPE with the onset on July 14, 2000 and the results are compared with our results obtained previously for the October 19, 1989 SPE. It is shown that not only electron density increases after SPE, but also the amount of clusters. It is found that the magnitude of the ionospheric response depends on season.     

Origin and energetics of the Io plasma torus

The ionization of the gas ejections from the Io satellite into the Jovian magnetosphere by the corotating magnetospheric plasma flow is considered. It is shown that the plasma flow velocity at the Io orbit exceeds the critical velocity at which the anomalous electron ionization of the heavy gas components takes place due to collisionless energy transfer from ionized gas atoms to plasma electrons. The energy, number density and spatial distribution of suprathermal electrons is calculated using the quasilinear theory of newly ionized atoms instability in a moving plasma. Saturation of the plasma density build up in a plasma is considered in terms of the instability quenching by Coulomb collisions.     

Cosmic ray influence on chemical composition of the atmosphere of the earth

Photochemical modeling and satellite data has been used to investigate the response of ozonosphere to solar proton events (SPEs) of the current (23d) maximum of solar activity. First SPE after the minima of solar activity occurred in November 1997. One of the strongest SPE after this occurred in July 2000. It was assumed in photochemical calculations that the ionization caused by solar protons in the atmosphere produced additional amount of the NOX and HOX compounds. Model runs showed strong ozone depletion in the mesosphere after SPE of 14 July 2000 for both polar regions. Corresponding study of ozone variations measured by HALOE instrument placed on board of UARS gave similar picture as the model calculations for North polar region.     

Ionospheric response to solar and magnetospheric protons during January 15–22, 2005: EAGLE whole atmosphere model results

We present an analysis of the ionosphere and thermosphere response to Solar Proton Events (SPE) and magnetospheric proton precipitation in January 2005, which was carried out using the model of the entire atmosphere EAGLE. The ionization rates for the considered period were acquired from the AIMOS (Atmospheric Ionization Module Osnabrück) dataset. For numerical experiments, we applied only the proton-induced ionization rates of that period, while all the other model input parameters, including the electron precipitations, corresponded to the quiet conditions. In January 2005, two major solar proton events with different energy spectra and proton fluxes occurred on January 17 and January 20. Since two geomagnetic storms and several sub-storms took place during the considered period, not only solar protons but also less energetic magnetospheric protons contributed to the calculated ionization rates. Despite the relative transparency of the thermosphere for high-energy protons, an ionospheric response to the SPE and proton precipitation from the magnetotail was obtained in numerical experiments. In the ionospheric E layer, the maximum increase in the electron concentration is localized at high latitudes, and at heights of the ionospheric F2 layer, the positive perturbations were formed in the near-equatorial region. An analysis of the model-derived results showed that changes in the ionospheric F2 layer were caused by a change in the neutral composition of the thermosphere. We found that in the recovery phase after both solar proton events and the enhancement of magnetospheric proton precipitations associated with geomagnetic disturbances, the TEC and electron density in the F region and in topside ionosphere/plasmasphere increase at low- and mid-latitudes due to an enhancement of atomic oxygen concentration. Our results demonstrate an important role of magnetospheric protons in the formation of negative F-region ionospheric storms. According to our results, the topside ionosphere/plasmasphere and bottom-side ionosphere can react to solar and magnetospheric protons both with the same sign of disturbances or in different way. The same statement is true for TEC and foF2 disturbances. Different disturbances of foF2 and TEC at high and low latitudes can be explained by topside electron temperature disturbances.     

Precision measurements of stratospheric ozone profiles by rocket-borne optical ozonesondes

Sunset observations of the upper stratospheric and mesospheric ozone were made at Uchinoura (31.25°N, 131.08°E) with rocket-borne optical ozonesondes, which consist of multi-color solar ultraviolet radiometers and sun tracking devices. Three ozone density profiles were obtained in this study. A comparison with the 30°N zonal and monthly average of the interim reference ozone model shows a variability that our present ozone mixing ratios below −50 km are larger in January and February and smaller in September than those of the model.     

太阳光球磁场特征物理量在质子事件短期预报中的应用

利用SOHO/MDI全日面纵向磁图, 计算了三个描述太阳活动区磁场复杂性和非势性的特征物理量, 即纵向磁场最大水平梯度Bz, 强梯度中性线长度L, 孤立奇点数目η. 为检验太阳光球磁场特征在质子事件短期预报中是否有效, 采用BP神经网络方法, 建立了基于这三个磁场特征物理量简单的太阳质子事件短期(24h)预报模型. 模 型在对2002年和2003年连续两年的样本检测中, 有很高的准确率(2002年和2003年 分别为90 %, 87.54 %)和较高的 质子事件报准率(2002年和2003年分别为60 %, 75 %),从而为光球磁场特征物理量作为质子事件预报的有效因子提供了依据.      

Solar particle events observed at Mars: dosimetry measurements and model calculations

During the period from March 13, 2002 to mid-September, 2002, six solar particle events (SPE) were observed by the MARIE instrument onboard the Odyssey Spacecraft in Martian Orbit. These events were observed also by the GOES 8 satellite in Earth orbit, and thus represent the first time that the same SPE have been observed at these separate locations. The characteristics of these SPE are examined, given that the active regions of the solar disc from which the event originated can usually be identified. The dose rates at Martian orbit are calculated, both for the galactic and solar components of the ionizing particle radiation environment. The dose rates due to galactic cosmic rays (GCR) agree well with the HZETRN model calculations.     

On the ionization ratio of the local interstellar medium

The abundance ratio of neutral hydrogen to neutral helium, as deduced from interplanetary observations of Lyman-alpha and He 584 A radiation by Mariner 10, is significantly lower than the cosmic abundance ratio of these elements, thus showing that the local interstellar medium (LISM) is partly ionized. It is shown that an important source of ionization of the LISM can be thermal collisions, yielding an ionization degree of about 50% for the hydrogen component.     

高层大气臭氧对于太阳紫外辐射27天波动的响应

本文利用Ｓｏｌａｒ Ｍｅｓｏｓｐｈｅｒｅ Ｅｘｐｌｏｒｅｒ（ＳＭＥ）卫星１９８２年和１９８３年太阳紫外辐射和赤道地区５０－９０ｋｍ臭氧分布的红外观测资料，对臭氧地太阳紫外辐射２７天振荡的响应进行了研究，并且利用大气光化波动模式对其进行了理论计算，计算结果与实测结果基本一致。     

Excitation of flare-induced waves in coronal loops and the effects of radiative cooling

EUV imaging observations from several space missions (SOHO/EIT, TRACE, and SDO/AIA) have revealed a presence of propagating intensity disturbances in solar coronal loops. These disturbances are typically interpreted as slow magnetoacoustic waves. However, recent spectroscopic observations with Hinode/EIS of active region loops revealed that the propagating intensity disturbances are associated with intermittent plasma upflows (or jets) at the footpoints which are presumably generated by magnetic reconnection. For this reason, whether these disturbances are waves or periodic flows is still being studied. This study is aimed at understanding the physical properties of observed disturbances by investigating the excitation of waves by hot plasma injections from below and the evolution of flows and wave propagation along the loop. We expand our previous studies based on isothermal 3D MHD models of an active region to a more realistic model that includes full energy equation accounting for the effects of radiative losses. Computations are initialized with an equilibrium state of a model active region using potential (dipole) magnetic field, gravitationally stratified density and temperature obtained from the polytropic equation of state. We model an impulsive injection of hot plasma into the steady plasma outflow along the loops of different temperatures, warm (～1 MK) and hot (～6 MK). The simulations show that hot jets launched at the coronal base excite slow magnetoacoustic waves that propagate to high altitudes along the loops, while the injected hot flows decelerate rapidly with heights. Our results support that propagating disturbances observed in EUV are mainly the wave features. We also find that the effect of radiative cooling on the damping of slow-mode waves in 1–6 MK coronal loops is small, in agreement with the previous conclusion based on 1D MHD models.     

Physical and model interpretation of HF radio propagation on the St. Petersburg–Longyearbyen (Svalbard) path

HF radio wave observations have been carried out with an oblique ionospheric sounding (OIS) method on the radio path from St. Petersburg to Longyearbyen (Svalbard), and experimental ionograms were obtained for December 2001. These ionograms have been analysed to investigate the impact of the main ionospheric trough (MIT) and magnetic disturbances on the signals on this path. The observations during weakly disturbed (K_р = 2) magnetic conditions on 14–15 December 2001 were compared with predictions from ray-tracing through a numerical model of the ionosphere. The ray-tracing computer program synthesizes the OIS ionograms by means of the “shooting method”. This method calculates trajectories of HF radio waves for different values of elevation angle and transmission frequency. There was a variety of calculated trajectories, from which we choose those which reach the receiver, and the selected paths provide a synthesis of the oblique ionograms. To simulate HF radio wave propagation, we apply a three-dimensional distribution of the electron density calculated with the mathematical model of the high-latitude ionosphere developed in the Polar Geophysical Institute (PGI). These numerical simulations permit us to interpret specific peculiarities of the OIS data such as abnormal propagation modes, increased delays of signals, enhanced MOF (maximum observed frequency) values etc. New results of the study are summarised as follows. (1) An unusual feature of the propagation along the path is the change of propagation mechanism during substorms on entering a path midpoint (or 1-hop reflection point) to the MIT. (2) Even weak substorms, having the distinguished intensities, lead to the appearance of different types of irregularities observed by the CUTLASS radar and therefore to the different propagation modes and F2MOF values. (3) The PGI model of the ionosphere was first used for ray-tracing at high latitudes. The model results are basically in a good qualitative agreement with experimental observations. This model provides the satisfactory agreement between the calculated and experimental F2MOF values while not correctly representing the fine structure of the experimental OIS ionograms at night. An agreement between the calculated and experimental data is better for day and evening hours than at night.     

Ozone trends in the mid-latitude stratopause region based on microwave measurements at Lindau (51.66° N, 10.13° E), the ozone reference model,and model calculations

We compared 8 years of ozone measurements taken at Lindau (51.66° N, 10.13° E) at altitudes between 40 and 60 km using the microwave technique with the CIRA ozone reference model that was established 20 years ago (Keating et al., 1990). We observed a remarkable decrease in ozone density in the stratopause region (i.e., an altitude of 50 km), but the decrease in ozone density in the middle mesosphere (i.e., up to 60 km in altitude) is slight. Likewise, we observed only a moderate decrease in the atmospheric region below the stratopause. Other studies have found the strongest ozone decrease at 40 km and a more moderate decrease at 50 km, which is somewhat in contradiction to our results. This decrease in ozone density also strongly depends on the season. Similar results showed model calculations using the GCM COMMA-IAP when considering the increase in methane. In the lower mesosphere/stratopause region, the strongest impact on the concentration of odd oxygen (i.e., O₃ and O) was observed due to a catalytic cycle that destroys odd oxygen, including atomic oxygen and hydrogen radicals. The hydrogen radicals mainly result from an increase in water vapor with the growing anthropogenic release of methane. The finding suggesting that the stratopause region is apparently attacked more strongly by the water vapor increase has been interpreted in terms of the action of this catalytic cycle, which is most effective near the stratopause and amplified by a positive feedback between the ozone column density and the ozone dissociation rate, thereby chemically influencing the ozone density. However, the rising carbon dioxide concentration cools the middle atmosphere, thereby damping the ozone decline by hydrogen radicals.     

Long-term tropospheric variations of ozone content caused by galactic cosmic ray influence

Galactic cosmic ray (GCR) fluxes measured by balloons in the troposphere and stratosphere at several points in Russia, and total ozone (TO) records have been used to detect cosmic signal by linear regression analysis. It was shown that TO response is in phase with decadal variations of GCR in contrast to the assumption about ozone destruction by GCR due to the nitrogen catalytic cycle intensification. 1-D photochemical model was used to understand the situation. The results of calculations show positive ozone response in the troposphere caused by additional production of NO by GCRs.     

电离层Es离子分布数值建模

从离子连续性方程和动量方程出发,比较全面地考虑了太阳光电离、星光电离、地冕电离、星际背景电离和流星离子流等电离源,综合分析Es层的主要动力学过程和光化学过程,以风剪切理论为基础,研究中性成分碰撞、电场作用、金属离子作用和分子离子作用等机制对离子分布的影响,建立了一维时变电离层Es物理模型.对离子产生率、垂直方向的离子速度在高度和时间上的变化进行仿真和计算,得到了在电场及风剪切作用下的离子密度剖面24h内的变化规律.根据建立的模型,以昆明相干散射雷达观测的径向风结果作为模型输入参考,仿真并得到了电离层Es电子密度的时空分布,据此反演出电离层Es临界频率f₀E_s.与昆明站电离层测高仪同时间观测的情况进行对比,结果比较一致,验证了建立的Es层物理模型正确性.      

Periodic variations of ionization rate in the low and middle atmosphere

The cosmic ray ionization source functions which were obtained using a simplified extensive air shower model are used to calculate the eleven year cycle, seasonal and diurnal variations of ionization rate in the low and middle atmosphere. The ionization source function, as a function of the penetrating depth and the energy of cosmic ray particles, is the ionization rate per unit depth for a unit flux of incoming cosmic ray particles with certain energy.The calculation of the eleven year cycle variation of ionization rate in the low and middle atmosphere due to the modulation of galactic cosmic ray intensity by solar activity shows that the amplitude is larger at a higher magnetic latitude and is generally larger at higher altitudes. The relative amplitude of fluctuation of the ionization peak value (at altitudes near 15 km) is up to 45% in the magnetic polar region. The ionization rate, due to the seasonal variation of the atmospheric density, varies from several per cent below the ionization peak to several tens per cent above the peak. This seasonal variation of ionization rate reaches 35% at 70 km. The diurnal variation of atmospheric densities caused by atmospheric tidal oscillation can produce a diurnal variation of the ionization rate to an amplitude of several per cent at altitudes above 40 km. The diurnal oscillation is less than 1% below 35 km.     

一种全球临近空间大气密度建模方法及应用

基于TIMED/SABER卫星2002—2018年观测的20~100 km大气密度数据,统计获得多年月平均值和标准偏差的全球网格数据。利用网格数据,分析了大气密度的变化特征。以网格数据为基准,计算了USSA76的相对偏差,分析了USSA76相对偏差的分布特征。以网格数据为驱动,将大气密度表征为平均值与大尺度扰动量和小尺度扰动量的加和,大尺度扰动和小尺度扰动分别采用余弦函数和一阶自回归模型表征,初步建立了全球临近空间大气密度模型。通过对比模型仿真值与激光雷达观测值,表明模型仿真值与观测值具有较好的吻合度,验证了建模方法的可行性。利用蒙特卡罗方法可再现给定轨迹上所有可能的大气状态。      

Ozone chemical equilibrium in the extended mesopause under the nighttime conditions

For retrieval of atomic oxygen and atomic hydrogen via ozone observations in the extended mesopause region (～70–100?km) under nighttime conditions, an assumption on photochemical equilibrium of ozone is often used in research. In this work, an assumption on chemical equilibrium of ozone near mesopause region during nighttime is proofed. We examine 3D chemistry-transport model (CTM) annual calculations and determine the ratio between the correct (modeled) distributions of the O₃ density and its equilibrium values depending on the altitude, latitude, and season.The results show that the retrieval of atomic oxygen and atomic hydrogen distributions using an assumption on ozone chemical equilibrium may lead to large errors below ～81–87?km. We give simple and clear semi-empirical criterion for practical utilization of the lower boundary of the area with ozone’s chemical equilibrium near mesopause.     

Mesospheric ozone densities deduced from solar occultation measurements

Quartz-UV occultation measurements by the satellite Interkosmos-16 have been used to calculate ozone densities at altitudes between 50 and 90 km for the period August to October 1976. Below 65 km densities agree well with the Krueger-Minzner-model. Mesopause densities have been studied in some detail. A certain percentage of the profiles show close correlation with the model of Shimazaki and Laird (with a pronounced minimum below the mesopause) while others fit better to the Park and London model (no minimum). This variability of the ozone density may be caused by different processes in the photo-chemistry of ozone. Two possible causes, the temperature dependent rate coefficients and the odd hydrogen processes are discussed in greater detail.     

Long term solar induced variations in total ozone, stratospheric temperatures and the tropopause

The variations of total ozone, stratospheric temperature and tropopause temperature are presented for the past 3 solar cycles for the summer months of the northern hemisphere. Ground-based, 30-year total column ozone series, filtered from its seasonal, QBO, El Nino/Southern Oscillation (ENSO) and trend components are found to be correlated to the 11-year solar cycle. Model calculations with a 2D chemical transport model are consistent with the observations. Mean stratospheric temperature variations, between levels 100 and 10 hPa, show also the same variation, correlated with the observed 11-year solar cycle, and the tropopause temperature increases in the same manner, in response to a warmer stratosphere during solar maxima.     

Model simulations of the impact of the 27-day solar rotation period on stratospheric ozone and temperature

Results are presented from two-year simulations of the effects of short-term solar ultraviolet (UV) variability using the Met. Office coupled chemistry-climate model. The model extends from the ground to 0.1 mbar and contains a complete range of chemical reactions allowing representation of all the main ozone formation and destruction processes in the stratosphere. The simulations were achieved by incorporating a 27-day oscillation in the pre-calculated model photolysis rates. Amplitudes for this signal were determined using solar spectral UV observations from the SOLar STellar Irradiance Comparison Experiment (SOLSTICE) instrument. Two experiments were carried out, one in which the UV variability was included in both the photolysis and radiation schemes and one in which only the photolysis scheme was modified.

The model reproduced several main features of observed correlations between short-term solar UV variability and both ozone and temperature in the tropical upper stratosphere, including the downward propagation of the phase lag and sensitivities of ozone and temperature to solar UV which are similar in magnitude to those observed. In the lower stratosphere, the ozone response to solar UV variability has not been well characterised from observations. Both model runs show a reversal of the propagation of phase lag below 10mb. The model response was found to be different between the two runs indicating that radiatively induced dynamical effects may play a significant role in the ozone response to solar UV variability.