Nonlinear model of short-scale electrodynamics in the auroral ionosphere

The optical detection of auroral subarcs a few tens of m wide as well as the direct observation of shears several m/s per m over km to sub km scales by rocket instrumentation both indicate that violent and highly localized electrodynamics can occur at times in the auroral ionosphere over scales 100 m or less in width. These observations as well as the detection of unstable ion-acoustic waves observed by incoherent radars along the geomagnetic field lines has motivated us to develop a detailed time-dependent two-dimensional model of short-scale auroral electrodynamics that uses current continuity, Ohms law, and 8-moment transport equations for the ions and electrons in the presence of large ambient electric fields to describe wide auroral arcs with sharp edges in response to sharp cut-offs in precipitation (even though it may be possible to describe thin arcs and ultra-thin arcs with our model, we have left such a study for future work). We present the essential elements of this new model and illustrate the models usefulness with a sample run for which the ambient electric field is 100 mV/m away from the arc and for which electron precipitation cuts off over a region 100 m wide. The sample run demonstrates that parallel current densities of the order of several hundred A m^-2 can be triggered in these circumstances, together with shears several m/s per m in magnitude and parallel electric fields of the order of 0.1 mV/m around 130 km altitude. It also illustrates that the local ionospheric properties like densities, temperature and composition can strongly be affected by the violent localized electrodynamics and vice-versa.     

Upward acceleration of magnetospheric ions by oscillatory centrifugal force

The paper addresses the issue of upward acceleration of ions along geomagnetic field lines. It has been shown that ion acceleration by electric field oscillations (formerly known as magnetic moment “pumping” or MMP) may be treated as a centrifugal acceleration mechanism. More precisely, the case in point is oscillatory centrifugal acceleration; this brings up the question on comparing the MMP with the centrifugal acceleration caused by the quasi-static magnetospheric convection field. It has been found that at high geomagnetic latitudes, the oscillatory centrifugal force is weaker or stronger than the centrifugal force of magnetospheric convection if the ratio of the electric field oscillation amplitude to the mean field is correspondingly lower or higher than \(\sqrt 2 \). Analysis of data from measurements and calculations of magnetospheric electric fields suggests that, contrary to current opinion, the oscillatory centrifugal force may be comparable to the centrifugal force of magnetospheric convection and even exceed it when strong global Pc5 pulsations are excited in the magnetosphere.     

Outflowing ionospheric oxygen-ion motion in a reconfigurating magnetosphere

During substorms, large-scale changes of the topology of the Earths magnetosphere following the variation of the characteristics of the interplanetary medium are accompanied by the induction of the electric field. In this study a model of a time-dependent magnetosphere is constructed and the large-scale features of the induced electric field are described. Local-time sectors with upward or downward field-aligned component and with intense perpendicular component of the electric field are distinguished. The electric-field structure implies the existence of outflow regions particularly effective in ion energization. With the vector potential adopted in the study, the region from which the most energized ions originate is defined by the local-time sector near 2100 MLT and latitude zone near 71° MLAT. The motion of ionospheric oxygen ions of energy 0.3–3 keV is investigated during a 5-min reconfiguration event when the tail-like magneto-spheric field relaxes to the dipole-like field. As the characteristics of plasma in the regions near the equatorial plane affect the substorm evolution, the energy, pitch angle, and the magnetic moment of ions in these regions are analyzed. These quantities depend on the initial energy and pitch angle of the ion and on the magnetic and electric field it encounters on its way. With the vector potential adopted, the energy attained in the equatorial regions can reach hundreds of keV. Three regimes of magnetic-moment changes are identified: adiabatic, oscillating, and monotonous, depending on the ion initial energy and pitch angle and on the magnetic- and electric-field spatial and temporal scales. The implications for the global substorm dynamics are discussed.     

Numerical modeling of the light ion trough and heat balance of the topside ionosphere in quiet geomagnetic conditions

The results from the numerical calculations of the global distribution of topside ionospheric parameters such as H⁺ ions and ion and electron temperatures up to 1500 km height are presented for equinoctial conditions at solar minimum. Calculations are carried out using the Global Self-consistent Model of Thermosphere, Ionosphere and Protonosphere (GSM TIP) developed in WD IZMIRAN, and using a new calculation block for electric fields due to dynamo and of magnetospheric origin. A comparison of two sets of calculations of magnetospheric convection electric field for a given potential difference is carried out, one through polar caps and other through field aligned currents of first zone. It is shown that the distribution of the electric potential obtained through field aligned currents of first zone is more self-consistent than that through polar caps. The light ion trough in H⁺ ions is deeper and occupies larger region for the potential difference through polar cap. For a given potential difference through field aligned current, at 1500 km, the maximum ion temperature is 150 K higher, minimum ion temperature is 200 K lower and maximum electron temperature is 100 K higher than those obtained for the same potential difference through polar caps. It is concluded that for modeling the electric field of magnetospheric origin, it is necessary to use the potential difference through field aligned current of first zone instead of through polar caps.     

Equatorial spread-F (ESF) and vertical winds

The Equatorial Spread-F (ESF) phenomenon is recorded in ionograms as a hierarchy of plasma instabilities in the F-layer of the equatorial ionosphere. The ESF is characterized by irregularities in the plasma (electron and ion) density and electric field distributions perpendicular to the Earth’s magnetic field. Large scale irregularities are generated by a primary plasma instability that develops in electric fields and plasma densities. Other secondary instabilities then develop and generate irregularities at several scale sizes that often produce a plasma ‘hole’ or ‘bubble’ that rises up with high E×B velocities. The ESF/plasma bubble phenomenon has been studied extensively with experimental techniques and modeling, which revealed important features. In the bottom side F-layer, near sunset, when the vertical density gradient steepens as the layer is supported by the horizontal (North–South) Earth’s magnetic field lines against the omnipresent Earth’s gravitational acceleration (g), the plasma conditions can give rise to Rayleigh–Taylor (RT) type instability. But the observed day to day variability of the ESF occurrence suggested that other agencies may also be involved in generating the instability. Sekar and Raghavarao (1987) with linear theory, and Raghavarao, Sekar and Suhasini (1992), with non-linear numerical modeling, suggested that vertical downward (upward) winds in the ambient gas have the potential to cause (inhibit) the ESF/bubble phenomenon. The presence of downward winds near the equator was reported earlier. In this paper, we show evidence for the presence of downward winds collocated with irregularities in electric fields and plasma densities as revealed by an unique combination of highly accurate measurements with instruments onboard the DE-2 satellite. The observations reported here are also consistent with the notion that the build-up of the equatorial ionization anomaly (EIA) prior to local sunset is important for the ESF instability.     

The equatorial electrojet

The typical quiet day variations of the equatorial electrojet (EEJ) current intensity with time of the day, season, sunspot number, and geomagnetic latitude are presented in terms of the corresponding variations of H which is the deviation of the horizontal component (H) of the geomagnetic field from its steady nighttime level. The observed height structure of the current density in the EEJ as measured in rocket flights is presented, along with the theoretically computed structure. Theoretical model results on the polarization electric fields and east-west currents as generated by the local interactions of height-varying winds in the EEJ show large height gradients and reversals for both currents and electric fields; experimental evidence for the reality of such height structures is also shown. The characteristics of the counter-electrojet events are presented and the possible causative mechanisms are discussed critically.Some typical experimental results are presented on the electric field changes in the EEJ which result from its sensitive response to electrodynamic disturbances in the magnetosphere and the auroral-polar latitude ionosphere during geomagnetic substorms and storms; and their implications are discussed. Possibilities for utilizing the EEJ as a very useful medium for important scientific studies on the larger space domain of ionosphere-magnetosphere system, on plasma waves, and on the earth's conductivity are emphasized.     

Estimation of the value of the electric field in the upper ionosphere before an earthquake

Summary Plasma parameters were measured by the Intercosmos 24 satellite at altitudes between 2300 km and 2500 km about 40 hours prior to the strong shock of the Iranian earthquake of 6 July 1990. Anomalous concentrations of light ions and electrons were observed mainly above the active seismic zone up to two days before the shock (Boková et al., 1993, 1994). This paper presents a model which, with the use of the anomalous electric field, is able to describe the observed changes of light ion concentrations in the plasmasphere above regions of imminent seismic events. The value of the intensity of the anomalous electric field, required to explain the observed changes, is estimated. For the seismic event analyzed intensity of the anomalous electric field is comparable with the value of the ambipolar electric field, but of opposite direction.     

Interplanetary control of thermospheric densities during large magnetic storms

During the main phase of large magnetic storms significant energy can be deposited in the ionosphere but produce no commensurate magnetic perturbations on the ground. Consequently, models designed to predict and specify thermospheric energy budgets based on ground magnetic data are negatively impacted. To quantify these effects we compare thermospheric densities predicted by the MSIS model with those inferred from accelerometer measurements by the Gravity Recovery and Climate Experiment (GRACE) satellites during two magnetic storm periods in 2004. Although predictions and measurements are in substantial agreement during quiet times, the model significantly underpredicts densities during storms. Also, the model's maxima occur several hours after observed stormtime peaks. We show that polar cap potentials and magnetospheric electric fields derived from interplanetary parameters measured by the Advanced Composition Explorer satellite are roughly proportional to neutral densities observed by GRACE with lead times of ∼4 h. Finally, ion drift meter data from Defense Meteorological Satellite Program spacecraft suggest that unpredicted positive and negative spikes found in high latitude accelerometer data reflect encounters with strong head and tail thermospheric winds driven by anti-sunward convecting plasma.     

The effect of velocity veering on sand transport in a shallow sea

This study aims at comparing and contrasting two different models for sand transport by currents in a shallow sea to illustrate the effect of velocity veering. The first model uses the Bailard-type formulation, which allows calculation of erosion/deposition rates at a fixed location on the sea floor via the divergence of horizontal sediment fluxes. The second model is a semi-analytical 2.5-dimensional model, which takes into account the time lag between erosion and deposition events and the velocity veering within the sediment-laden (nepheloid) layer caused by the Coriolis force. The velocity veering implies that the direction of the sediment flux is generally different from the direction of the surface flow. The latter model was designed for rapid, semi-analytical computations of sediment transport, using flow fields from 2-DH numerical models. The two models use a matching set of parameters to provide identical values for the bottom stress and suspended sediment load for a uniform steady current at any given surface velocity. The two models were compared in a range of sand grain sizes 50–500 m and current speeds up to 1 m s^–1 for an idealised square region (100 × 100 km) of a shelf sea of constant depth. The erosion/deposition patterns and suspension load were examined in three settings: (1) uniform steady flow, (2) straight jet, (3) meandering jet. It was found that both the rates and, in particular, the spatial distribution of the areas of erosion/deposition differ significantly between the models in cases (2) and (3). This difference can be attributed to additional flux divergence due to velocity veering. A comparison of model results with field data, collected at Long Island Shelf, supports the relevance of Coriolis-induced veering of currents on the direction of the sediment flux.Responsible Editor: Jens Kappenberg     

Disturbance of the ionospheric <Emphasis Type="Italic">D</Emphasis> region by the electric current of the atmospheric-ionospheric electric circuit

The mechanism by which electron and ion densities change in the ionospheric D region due to the electric current flowing in the atmospheric-ionospheric electric circuit is studied. The current disturbance in this circuit exists over the regions of increased seismic, meteorological, and thunderstorm activity. In the framework of the model, the influence of the electron and ion transportation under the action of the electric field on the formation of a disturbance in the D region and heating of the plasma electron component by the field are considered. The calculation results show that the densities of electrons and ions can change by an order of magnitude at an increase in the current density up to ∼(10⁻⁹–10⁻⁸) A m⁻², the sign of the disturbance depending on the current direction.     

Large-scale E-region electric field structure due to gravity wave winds

Oscillatory neutral wind motions, such as those of atmospheric gravity waves which propagate through the E-region of the ionosphere, appear to produce local electric fields in the source region. Although the net effect of these oscillatory fields vanishes outside the source region, the local fields themselves are not shorted out along the magnetic field lines, as is usually assumed. We present in situ measurements of neutral winds and the correlated electric fields, and show that local electric fields of the order of a few mV/m can be sustained by the neutral wind motions.     

Equatorial ionization anomaly and thermospheric meridional winds during two major storms over Brazilian low latitudes

The equatorial ionosphere responses over Brazil to two intense magnetic storms that occurred during 2001 are investigated. The equatorial ionization anomaly (EIA) and variations in the zonal electric field and meridional winds at different storms phases are studied using data collected by digisondes and GPS receivers. The difference between the F layer peak density (foF2) at an equatorial and a low latitude sites was used to quantify the EIA; while the difference between the true heights (h_F) at the equatorial and an off-equatorial site was used to calculate the magnetic meridional winds. The vertical drift was calculated as dh_F/dt. The results show prompt penetration electric fields causing unusual early morning development of the EIA, and disturbed dynamo electric field producing significant modification in the F region parameters. Variations to different degrees in the vertical drift, the thermospheric meridional winds and the EIA developments were observed depending on the storm phases.     

Analysis of the positive ionospheric response to a moderate geomagnetic storm using a global numerical model

Current theories of F-layer storms are discussed using numerical simulations with the Upper Atmosphere Model, a global self-consistent, time dependent numerical model of the thermosphere-ionosphere-plasmasphere-magnetosphere system including electrodynamical coupling effects. A case study of a moderate geomagnetic storm at low solar activity during the northern winter solstice exemplifies the complex storm phenomena. The study focuses on positive ionospheric storm effects in relation to thermospheric disturbances in general and thermospheric composition changes in particular. It investigates the dynamical effects of both neutral meridional winds and electric fields caused by the disturbance dynamo effect. The penetration of short-time electric fields of magnetospheric origin during storm intensification phases is shown for the first time in this model study. Comparisons of the calculated thermospheric composition changes with satellite observations of AE-C and ESRO-4 during storm time show a good agreement. The empirical MSISE90 model, however, is less consistent with the simulations. It does not show the equatorward propagation of the disturbances and predicts that they have a gentler latitudinal gradient. Both theoretical and experimental data reveal that although the ratio of [O]/[N₂] at high latitudes decreases significantly during the magnetic storm compared with the quiet time level, at mid to low latitudes it does not increase (at fixed altitudes) above the quiet reference level. Meanwhile, the ionospheric storm is positive there. We conclude that the positive phase of the ionospheric storm is mainly due to uplifting of ionospheric F2-region plasma at mid latitudes and its equatorward movement at low latitudes along geomagnetic field lines caused by large-scale neutral wind circulation and the passage of travelling atmospheric disturbances (TADs). The calculated zonal electric field disturbances also help to create the positive ionospheric disturbances both at middle and low latitudes. Minor contributions arise from the general density enhancement of all constituents during geomagnetic storms, which favours ion production processes above ion losses at fixed height under day-light conditions.     

二维电离层发电机理论模式及其初步应用

在地理坐标系下推导出二维电离层发电机理论方程,采用逐线迭代法求解得到全球二维电离层发电机电流函数,进而得到电离层发电机电流和电场.模式中使用的电导率是根据外部经验模式给出的背景大气和电离层参数,采用理论公式计算得出;输入的中性风场和磁场分别由HWM93和IGRF2000模型给出,该电离层发电机理论模式很好地给出了全球Sq电流形态及电离层E层发电机电场的基本特征.利用该模式研究了外部模式风场以及地磁场随高度的变化对模拟结果的影响,发现在90~180 km高度上,风场随高度变化对电流影响较大,而地磁场影响较小;重点模拟研究了地磁平静时期,Sq电流涡旋中心位置和总电流强度的变化规律,初步研究发现,电流中心位置在地理纬度±30°附近,不同的地方时电流随地磁纬度线平行移动,且南北半球两个电流涡中心电流强度之和变化不大.分析发现这种规律与发电机高度上的磁场总强度及地磁倾角的全球分布有很好的相关性.     

Three-dimensional Simulations of the Mean Air Transport During the 1997 Forest Fires in Kalimantan, Indonesia Using a Mesoscale Numerical Model

— This paper describes the meteorological processes responsible for the mean transport of air pollutants during the ENSO-related forest fires in Kalimantan, Indonesia from 00 UTC 21 September to 00 UTC 25 September, 1997. The Fifth Generation of the Pennsylvania State University-National Center for Atmospheric Research (PSU-NCAR) Mesoscale Model (MM5) is used to simulate three-dimensional winds at 6-hourly intervals. A nonhydrostatic version of the model is run using two nested grids with horizontal resolutions of 45 km and 15 km. From the simulated wind fields, the backward and forward trajectories of the air parcel are investigated using the Vis5D model.¶The results indicate that the large-scale subsidence over Indonesia, the southwest monsoon low-level flows (2–8 m s^?1), and the shallow planetary boundary layer height (400–800 m) play a key role in the transport of air pollutants from Kalimantan to Malaysia, Singapore and Brunei.