Calculation of X-ray polarization during flares

The polarization of X-radiation emitted by electrons which are accelerated during flares is investigated in a non-thermal model. Monte Carlo calculations combining analytically treated multiple-scattering and random large-angle scattering processes yield the energy and angular distributions of flare electrons penetrating the upper solar atmosphere. The X-ray polarization derived from these distributions is computed as a function of energy, observation angle and column density traversed by the electrons. Initially monoenergetic electrons as well as electrons with a power-law spectrum are considered. The degree of polarization of the total X-radiation obtained by summing over all layers of the atmosphere is compared with observations.     

内重力波传播的3维传输函数模式研究

在考虑背景风场及大气耗散的条件下,建立了3维内重力波传输函数数值模式.分析了300 km高度3维传输函数在频率波数域的特性,并以近地面单位脉冲点源为激发源,得到了内重力波在3维空间中的时空分布.讨论了不同时空尺度地面方波源激发的内重力波在电离层高度的能量分布特征.结果表明,(1)对内重力波而言,背景大气相当于一个带通滤波器,只有波动周期和波长分别在15～30 min和200～400 km之间的重力波扰动最容易上传到300km高度;(2)在背景风场的作用下水平面上以同心圆扩散的波阵面以及垂直方向上成漏斗状的波阵面发生了变形,并且逆风方向比顺风方向更有利于声重力波由对流层向电离层高度传播;(3)300km高度对时间尺度和空间尺度分别在20～30 min和150～250 km之间的地面方波源响应的总能量最强.      

Atmospheric dynamics on Venus and Mars

New results from Pioneer Orbiter observations indicate a continued vortex organization of the cloud level atmosphere in either hemisphere, centered over respective poles. Significant changes in the magnitude of the cloud level zonal circulation over a period of several years have been detected. A strong signature of the solar tidal circulation has been detected in the atmospheric circulation with the lowest speeds occurring in equatorial latitudes about 20° upstream of the sub-solar point. Finally, a solar-locked persistent spatial structure has been discovered in the variance of the ultraviolet brightness measured from brightness normalized images of Venus. Vega balloons (drifting at about 53 km altitude near 7°N and 7°S latitudes) have also provided some unique observations of atmospheric circulation, significant among them being the strong vertical motions, the zonality of their drift speeds as well as a significant temperature difference between the two balloons. The temperature difference which amounts to 6.5°K on average is currently being interpreted as a temperature variation with longitude or time.

Diagnostic modelling efforts towards simulating the atmospheric circulation on Venus are continuing and have provided some clues about the processes that maintain them but have not yet been successful in explaining the superrotation of the atmosphere.

Knowledge of the Martian atmospheric dynamics on the other hand is still limited by lack of adequate observations. Numerical modelling of the Martian atmosphere continues to provide most of the information about the atmospheric circulation. The situation regarding the paucity of observations should improve with the completion of the proposed Mars Observer mission. The low circular polar orbit planned provides an excellent opportunity to study the Martian atmosphere.     

Enhancement of High Energy Electron Fluxes and Variation of Atmospheric Electric Field in the Antarctic Region ormalsize

High-energy electron precipitation in the high latitude regions enhances the ionization of the atmosphere, and subsequently increases the atmospheric conductivities and the vertical electric field of the atmosphere near the ground as well. The High-Energy Electron Flux (HEEF) data measured by the Fengyun-3 meteorological satellite are analyzed together with the data of near-surface atmospheric vertical electric field measured at the Russian Vostok Station. Three HEEF enhancements are identified and it is shown that when the HEEF increases to a certain level, the local atmospheric vertical electric field near the ground can increase substantially than usual. The response time of the electric field to HEEF enhancement is about 3.7 to 4 days.      

Photochemical response of neutral and ionized middle atmosphere composition to the strong solar proton event of October 1989

One of the strongest solar proton events (SPE) occurred in October 1989. Its forcing of the middle atmosphere chemistry including ionized components in the D-region is examined. The ionization rate, and ozone, NO and OH density temporal and spatial (vertical) deviations induced by the SPE, calculated by a 1-D time-dependent photochemical model separately for daytime and nighttime (not shown here), are used in a 1-D model of the lower ionosphere to calculate the response of ionized components to combined forcing by ionization rate and neutral chemical composition disturbances. The radio wave absorption caused by electron density disturbances after the SPE is calculated and compared with observations. The computed ozone values are compared with observations, as well.     

模拟三维重力波传播过程的并行数值模式

从三维大气运动的基本控制方程出发, 基于MPI消息传递接口和区域分解的思想, 建立了模拟三维可压缩大气中重力波传播过程的并行数值模式. 在对垂直方向进行区域分解的基础上, 针对跳点网格的特点, 以垂直速度主格点线为子区域的下边界, 状态变量主格点线为子区域的上边界进行区域分解. 利用MPI消息传递接口来传递计算各子区域中变量时所要用到的相邻子区域上相应变量, 从而顺利地完成对整个计算域的并行计算. 根据线性重力波理论, 通过模拟小振幅重力波的传播过程对所建模式进行了验证. 结果表明, 模式可以很好地模拟小振幅重力波在三维空间中的传播过程, 模拟的重力波振幅随着传播高度的增加以指数形式增长, 对能量传播路径以及有效扰动位能和扰动动能的模拟结果均与理论预测的结果吻合很好. 另外, 随着进程个数的增加, 完成相同的计算量所用的计算时间也显著减少. 这表明本文建立的并行数值模式不但能够很好的模拟重力波的传播过程, 而且能够有效地节约计算时间.      

大气重力波引起的突发钠层的模拟研究

应用钠原子和中性大气分子的质量连续性方程来模拟突发纳层（SSL），垂直风场采用接近实际大气重力波的正弦行波模式，结果较好地反映了SSL的形成过程。SSL的形成时刻在5－15min之间并可持续到30min之后，形成高度大约在90－100km之间，峰宽为0.5-2km之间，这些都与实际观测SSL的特点相符，同时还进一步地研究了当重力波参数（主要指垂直波长和周期）、风速以及常态钠层半宽度发生变化时SSL的变化趋势。     

Enhancement of High Energy Electron Fluxes and Variation of Atmospheric Electric Field in the Antarctic Region

High-energy electron precipitation in the high latitude regions enhances the ionization of the atmosphere,and subsequently increases the atmospheric conductivities and the vertical electric field of the atmosphere near the ground as well.The High-Energy Electron Flux(HEEF) data measured by the Fengyun-3 meteorological satellite are analyzed together with the data of nearsurface atmospheric vertical electric field measured at the Russian Vostok Station.Three HEEF enhancements are identified and it is shown that when the HEEF increases to a certain level,the local atmospheric vertical electric field near the ground can increase substantially than usual.The response time of the electric field to HEEF enhancement is about 3.7 to 4 days.     

E region electric fields at the dip equator and anomalous conductivity effects

Zonal and vertical electric fields were estimated at E region heights in the Brazilian sector. Zonal electric fields are obtained from the vertical electric fields based on their relation through the Hall-to-Pedersen ionospheric conductivities ratio. The technique for obtaining the vertical electric field is based on its proportionality to the Doppler velocities of type 2 irregularities as detected by coherent radars. The 50 MHz backscatter coherent (RESCO) radar was used to estimate the Doppler velocities of the type 2 irregularities embedded in the equatorial electrojet. A magnetic field-line integrated conductivity model was developed to provide the conductivities. It considers a multi-species ionosphere and a multi-species neutral atmosphere, and uses the IRI 2007, the MISIS 2000 and the IGRF 10 models as input parameters for ionosphere, neutral atmosphere and Earth’s magnetic field, respectively. The ion-neutral collision frequencies of all the species are combined through the momentum transfer collision frequency equation, and different percentages of electron-neutral collisions were artificially included for studying the implication of such increase in the zonal electric field, which resulted ranging from 0.13 to 0.49 mV/m between the 8 and 18 h (LT), under quiet magnetic conditions.     

Analysis of NASA’s DSN Venus Express radio occultation data for year 2014

The Venus Express Radio Science Experiment (VeRa) was part of the scientific payload of the Venus Express (VEX) spacecraft and was targeted at the investigation of Venus’ atmosphere, surface, and gravity field as well as the interplanetary medium. This paper describes the methods and the required calibrations applied to VEX-VeRa raw radio occultation data used to retrieve vertical profiles of Venus’ ionosphere and neutral atmosphere. In this work we perform an independent analysis of a set of 25 VEX, single-frequency (X-band), occultations carried out in 2014, recorded in open-loop at the NASA Deep Space Network. Our temperature, pressure and electron density vertical profiles are in agreement with previous studies available in the literature. Furthermore, our analysis shows that Venus’ ionosphere is more influenced by the day/night condition than the latitude variations, while the neutral atmosphere experiences the opposite. Our scientific interpretation of these results is based on two major responsible effects: Venus’ high thermal inertia and the zonal winds. Their presence within Venus’ neutral atmosphere determine why in these regions a latitude dependence is predominant on the day/night condition. On the contrary, at higher altitudes the two aforementioned effects are less important or null, and Venus’ ionosphere shows higher electron density peaks in the probed day-time occultations, regardless of the latitude.     

地形产生的山地波及其传播过程模拟研究

气流经过地形产生的山地波是大气重力波研究的重要类别之一.从大气运动的控制方程组出发,建立模拟地形产生的山地波及其传播过程的二维数值模式.利用水平背景风场、地形和垂直速度之间的关系,在模式中引入垂直速度扰动作为地形产生山地波的激发源.通过模拟该激发源引起的扰动即山地波在大气中的传播过程,再现了山地波的产生、传播及充分发展过程.通过分析水平波长、垂直波长、位温扰动、流线,在空间尺度上描述了山地波的产生、传播及充分发展的过程.在气流经过地形产生的山地波的传播过程中,其水平波长λ_x的范围为2.5~5km,垂直波长λ_z约为2.5km.这些结果与利用山地波线性理论计算的垂直波长一致,从而验证了本模式能够模拟地形产生的山地波及其传播过程,为深入了解山地波的产生过程及其在中高层大气中的传播机制和效应奠定了基础.      

Interpretation of data on laser sensing of aerosol atmosphere from space

A single-frequency lidar, using aerosol scattering as an informative component, is the simplest and reliable facility for remote sensing of the atmosphere. The information on vertical distribution of atmospheric aerosol which can be obtained using such a lidar is necessary for investigating the physics of atmospheric processes and forecast of optical state of the atmosphere. At the same time, the interpretation of data on single-frequency sounding is associated with some difficulties of fundamental character, mainly due to insufficient software of the experiment. Under such conditions the problems of optimal processing of lidar returns aiming at extracting the useful information on aerosol are of great importance, especially if one takes into account the hindering effects of atmospheric background and optical noises. This paper presents a statistical approach to this problem, and the possibilities of single-frequency sounding from space are analyzed.     

太阳大气中非线性波的数字计算

本文用多方过程方程式封闭一维带引力场的流体力学方程组,用不同的γ值计算了不同类型的扰动在HSRA太阳大气模型中向上传播和演化。说明在太阳大气中向上传播的小扰动能够发展成为激波;激波的形成高度基本上不受热力学过程的影响,而与扰动强度关系很大,愈小的扰动传播愈高;但由于太阳大气中湍流场的存在,完全有可能使这样小的扰动在未形成激波前就被破坏掉。这就证实了一些人认为的对流层机械波并不是一个对色球和日冕的有效加热机制的说法。      

Thermal structure of the Martian atmosphere retrieved from the IR spectrometry in the 15 μm CO2 band: input to MIRA

This paper describes one of the sources of the data concerning the thermal structure of the Martian atmosphere, based on the thermal IR spectrometry method. It allows to investigate the Martian atmosphere below 55 km by retrieving the temperature profiles from the 15 μm CO₂ band. This approach enables to reach the vertical resolution of several kilometers and the temperature accuracy of several Kelvins. An aerosol abundance, which influences the temperature profile, is obtained from the continuum of the same spectrum parallel with the temperature profile and is taken into account in the temperature retrieval procedure in a self consistent way. Although this method has the limited vertical resolution, it possesses a significant advantage: the thermal IR spectrometry allows to monitor the temperature profiles with a good coverage both in space and local time. The Planetary Fourier spectrometer on board of Mars Express has the spectral range from 250 to 8000 cm⁻¹ and a high spectral resolution of about 2 cm⁻¹. Vertical temperature profiles retrieval is one of the main scientific goals of the experiment. The important data are expected to be obtained on the vertical thermal structure of the atmosphere, and its dependence on latitude, longitude, season, local time, clouds and dust loadings. These results should give a significant input in the future MIRA, being included in the Chapter “Structure of the atmosphere from the surface to 100 km”.     

Sterilisation properties of the Mars surface and atmospheric environment.

The radiative and chemical conditions at the surface and in the lower Martian atmosphere are computed at various latitudes and seasons combining a 2D photochemical model and radiation simulations. In most situations, the solar UV B and C radiations reach the surface however, suspended dust and, in polar cases, ozone can constitute an effective UV shield. The daytime and night time concentrations of the sterilizing oxidants: OH, H2O2 and O3 are determined, as well as the concentration of the substances which could influence the metabolism of microorganisms. The possible habitats of a remaining Mar's life as well as the possibilities of contamination by resistant earth life forms will be described.     

Absorbers seen near the Venus cloud tops from pioneer Venus

Spin-scan images from the Pioneer Venus Orbiter UV Spectrometer and the Cloud Photopolarimeter provide a set of planetary contrast measurements in the wavelength range 1990A to 3650A and phase angles from 33°–130°. The planet is darkest at the point where the UVS line of sight penetrates perpendicular to the cloud tops: thus the absorbing material responsible must be deep in the atmosphere. Sulfur dioxide absorption can explain the amount of contrast seen between 2000A and 3200A. At the longer wavelengths, the persistence of contrast requires another absorber which is deeper in the atmosphere and strongly associated with the location of the SO₂. Part of the observed contrast is due to the high-lying haze discovered from Pioneer Venus polarimetry. The correlation between planetary contrast and polarization does not support large scale clearing or major vertical motions of the cloud tops as the sole cause of the observed contrast. However, a scheme in which absorbers subject to photochemical destruction are mixed upward into the cloud top region provides a consistent explanation for the origin of these markings.     

Model for induced ionization by galactic cosmic rays in the Earth atmosphere and ionosphere

The ionization profiles produced by galactic cosmic rays in the Earth atmosphere and ionosphere are obtained on the basis of Monte Carlo simulations. Cascade processes in the atmosphere are simulated using CORSIKA 6.52 code with FLUKA 2006 and QGSJET II hadronic interaction subroutines. Proton induced showers are considered using a realistic atmospheric model (US Standard Atmosphere). The energy deposit from different components is taken into account, namely electromagnetic, hadron and muon components. The curvature of the atmosphere is considered in the computer code. On the basis of the computed ionization yield function the ion pair production rate in the atmosphere is obtained for different conditions and locations. The model is applicable to the entire atmosphere, from ground level to upper atmosphere. Several applications of the obtained results are discussed. The Monte Carlo simulation model considers nuclear interactions below the altitude of 35 km. It is compared with analytical–numerical electron production rate model. The latter model which takes into account the electromagnetic interactions above altitudes of 35 km has two main regions of application: above 50 km (thin target model) and between 35 and 50 km (intermediate target model). A good agreement between the CORSIKA results and analytical–numerical model results is found above altitude of 35 km.     

Importance of electric field measurement over low latitudes at stratospheric heights by balloons

The vertical field in the stratosphere around 35 km is predominantly of atmospheric origin whereas the horizontal electric field at these altitude is mainly of ionospheric origin. The electrical coupling between ionosphere and atmosphere is not known for low latitudes. Balloon borne electric field measurements are planned from Hyderabad, India (geographic latitude 17.5° N) to understand this coupling. Measurement of stratospheric electric fields are also important from the point of view of the sun-weather relationship. It si suggested that the balloon borne electric field measurements are important to understand the electrodynamics of the middle atmosphere.     

Wavelet modelling in support of IRI

The slant total electron content (STEC) of the ionosphere is defined as the integral of the electron density along the ray-path of the signal between the transmitter and the receiver. So-called geometry free GPS measurements provide information on the electron density, which is basically a four-dimensional function depending on spatial position and time. Since ground-based measurements are not very sensitive to the vertical structure within the atmosphere, the ionosphere is often represented by a spherical layer, where all electrons are concentrated. Then the STEC is transformed into the vertical total electron content (VTEC), which is a three-dimensional function depending on longitude, latitude and time.In our approach, we decompose an ionospheric function, i.e. the electron density or the VTEC, into a reference part computed from a given model like the International Reference Ionosphere (IRI) and an unknown correction term expanded in a multi-dimensional series in terms of localizing base functions. The corresponding series coefficients are calculable from GPS measurements applying parameter estimation procedures. Since the GPS receivers are located rather unbalanced, finer structures are modelable just in regions with a sufficient number of observation sites. Due to the localizing feature of B-spline functions we apply a tensor product spline expansion to model the correction term regionally. Furthermore, the multi-resolution representation derived from wavelet analysis allows monitoring the ionosphere at different resolutions levels. We demonstrate the advantages of this procedure by representing a simulated VTEC data set over South America.