Current and Future Applications of 3-D Global Earth-Ionosphere Models Based on the Full-Vector Maxwell’s Equations FDTD Method

Advances in computing technologies in recent decades have provided a means of generating and performing highly sophisticated computational simulations of electromagnetic phenomena. In particular, just after the turn of the twenty-first century, improvements to computing infrastructures provided for the first time the opportunity to conduct advanced, high-resolution three-dimensional full-vector Maxwell’s equations investigations of electromagnetic propagation throughout the global Earth-ionosphere spherical volume. These models, based on the finite-difference time-domain (FDTD) method, are capable of including such details as the Earth’s topography and bathymetry, as well as arbitrary horizontal/vertical geometrical and electrical inhomogeneities and anisotropies of the ionosphere, lithosphere, and oceans. Studies at this level of detail simply are not achievable using analytical methods. The goal of this paper is to provide an historical overview and future prospectus of global FDTD computational research for both natural and man-made electromagnetic phenomena around the world. Current and future applications of global FDTD models relating to lightning sources and radiation, Schumann resonances, hypothesized earthquake precursors, remote sensing, and space weather are discussed.

电离层不规则结构漂移的GPS测量及其初步结果

本文阐述了利用GPS接收机台阵测量到的闪烁和TEC变化率ROT快速起伏图样估计F层不规则结构漂移的原理和方法,并利用实测数据估计了静日和暴时电离层不规则结构的水平漂移速度.短间距台网和超短间距台链观测实例的计算结果表明,暴时武汉地区引起TEC快速起伏的电离层不规则结构沿纬圈向西漂移,21∶30至03∶00 LT,西向漂移速度在约40 m/s至130 m/s的范围内变化;在桂林地区,磁静日午夜前后引起L波段电波闪烁的电离层不规则结构沿纬圈向东漂移,漂移速度从约70 m/s下降到约55 m/s,磁扰日午夜前不规则结构向西漂移,速度从约150 m/s下降到约50 m/s,午夜后转为向东漂移,速度从约25 m/s上升到约65 m/s.文中还提出了由单站多卫星观测估计F层不规则结构漂移的设想.实例分析与计算结果表明,利用单站多卫星观测估计电离层不规则结构漂移是一种合理可行的方法.     

Global plasmaspheric TEC and its relative contribution to GPS TEC

The plasmaspheric electron content is directly estimated from the global positioning system (GPS) data onboard JASON-1 Satellite for the first time. Similarly, the ground-based GPS total electron content (TEC) is estimated using about 1000 GPS receivers distributed around the globe. The relative contribution of the plasmaspheric electron content to the ground-based GPS TEC is then estimated globally using these two independent simultaneous measurements; namely ground-based GPS TEC and JASON-1 GPS TEC. Results presented here include data from 3 months of different solar cycle conditions (October 2003, May 2005, and December 2006). The global comparison between the two independent measurements was performed by dividing the data into three different regions; equatorial, mid- and high-latitude regions. This division is essential as the GPS raypaths traverse different distances through the plasmasphere at different latitudes. The raypath length through the plasmasphere decreases as latitude increases. The relative contribution of the plasmaspheric electron content exhibits a diurnal variation that depends on latitude with minimum contribution (10%) during daytime and maximum (up to 60%) at night. The contribution is also maximum at the equatorial region where the GPS raypath traverses a long distance through the plasmasphere compared to its length in mid- and high-latitude regions. Finally, the solar cycle variation of plasmaspheric contribution is also reported globally.     

Response of the ionosphere–thermosphere system to magnetospheric processes

The magnetosphere–ionosphere–thermosphere system at high latitudes is strongly coupled via electric fields, particle precipitation, plasma and neutral outflows, and field-aligned currents. Although the climatology of the coupled system is fairly well established, our understanding of the variability of the disturbed state (weather) is rudimentary. This variability is associated with magnetic storms and substorms, nonlinear processes that operate over a range of spatial scales, time delays, and feedback mechanisms between the different domains. The variability and resultant structure of the ionosphere can appear in the form of propagating plasma patches and polar wind jets, pulsing ion and neutral polar winds, auroral and boundary blobs, and ionization channels associated with polar cap arcs, discrete auroral arcs, and storm-enhanced densities (SEDs). The variability and structure of the thermosphere can appear in the form of propagating atmospheric holes, neutral gas fountains, neutral density patches, and transient neutral jets. In addition, during periods of enhanced plasma convection, the neutral winds can become supersonic in relatively narrow regions of the polar cap. The spatial structure in the ionosphere–thermosphere system not only affects the local environment, but the cumulative effect of multiple structures may affect the global circulation and energy balance. A focused topical review of recent results in our modeling the variability and structure of the high-latitude ionosphere–thermosphere system is presented. This review was given at the Greenland Space Science Symposium (May 2007).     

Altitude variation of the plasmapause signature in the main ionospheric trough

The projection of the plasmapause magnetic-field lines to low altitudes, where the light-ion chemistry is dominated by O⁺, tends to occur near the minimum electron density in the main (midlatitude) electron density trough at night. With increasing altitude in the trough, where H⁺ emerges as the dominant ion on the low-latitude boundary, we have found cases where the plasmapause field lines are located on the sharp low-latitude side of the trough as expected if this topside ionosphere H⁺ distribution varies in step with the plasmapause gradient in the distant plasmasphere. These conclusions are based on near-equatorial crossings of the plasmapause (corresponding to the steep gradient in the dominant species H⁺) by the Explorer-45 satellite as determined from electric-field measurements by Maynard and Cauffman in the early 1970s and ISIS-2 ionospheric topside-sounder measurements. The former data have now been converted to digital form and made available at http://nssdcftp.gsfc.nasa.gov. The latter provide samples of nearly coincident observations of ionospheric main trough crossings near the same magnetic-field lines of the Explorer 45-determined equatorial plasmapause. The ISIS-2 vertical electron density profiles are used to infer where the F-region transitions from an O⁺ to a H⁺ dominated plasma through the main trough boundaries.     

On the spectrum of mid-latitude sporadic-E irregularities

We discuss the creation of mid-latitude sporadic-E plasma irregularities (with length-scales smaller than sporadic layer thickness) by the neutral atmosphere turbulence. Using fluid equations, the relation between plasma density fluctuations and the velocity field of neutrals is derived. After a brief discussion of the relevant neutral turbulence, the analytical expression for the power spectrum of plasma irregularities is obtained. This expression allows a power-law type of experimental irregularity spectra (the spectral index depends on sporadic-E characteristics) and possible departures in detail of the irregularity spectra from the power-law form to be explained. In addition, it allows us to make estimates of length-scales at which such departures must occur.     

Geomagnetic control of polar mesosphere summer echoes

Using observations with the ALOMAR SOUSY radar near Andenes (69.3°N, 16.0°E) from 1994 until 1997 polar mesosphere summer echoes (PMSE) have been investigated in dependence on geomagnetic K indices derived at the Auroral Observatory Tromsø (69.66°N, 18.94°E). During night-time and morning hours a significant correlation between the signal-to-noise ratio (SNR) of the radar results and the geomagnetic K indices could be detected with a maximum correlation near midnight. The correlation becomes markedly smaller in the afternoon and early evening hours with a minimum near 17 UT. This diurnal variation is in reasonable agreement with riometer absorption at Ivalo (68.55°N, 27.28°E) and can be explained by the diurnal variation of ionization due to precipitating high energetic particles. Therefore, a part of the diurnal PMSE variation is caused by this particle precipitation. The variability of the solar EUV variation, however, has no significant influence on the PMSE during the observation period.     

Letter to the editor Electric field fluctuations (25–35 min) in the midnight dip equatorial ionosphere

Measurements with a HF Doppler sounder at Kodaikanal (10.2°N, 77.5°E, geomagnetic latitude 0.8°N) showed conspicuous quasi-periodic fluctuations (period 25/35 min) in F region vertical plasma drift, V_z in the interval 0047/0210 IST on the night of 23/24 December, 1991 (A_p = 14, K_p < 4^–). The fluctuations in F region vertical drift are found to be coherent with variations in B_z (north-south) component of interplanetary magnetic field (IMF), in geomagnetic H/X components at high-mid latitude locations both in the sunlit and dark hemispheres and near the dayside dip equator, suggestive of DP2 origin. But the polarity of the electric field fluctuations at the midnight dip equator (eastward) is the same as the dayside equator inferred from magnetic variations, contrary to what is expected of equatorial DP2. The origin of the coherent occurrence of equatorial electric field fluctuations in the DP2 range of the same sign in the day and night hemispheres is unclear and merits further investigations.     

The generation of non aspect sensitive plasma density irregularities by field aligned drifts in the lower ionosphere

A theory of the generation of plasma density irregularities with virtually no aspect sensitivity, in the lower ionosphere at high latitudes, by electron drifts aligned with the geomagnetic field, is presented. The theory is developed through fluid equations in which the destabilising mechanism involves positive feedback from electron collisional heating. When field aligned electron drift speeds exceed a few km s^–1, this effect destabilises waves with wavelengths in excess of a few tens of metres in the lower E-region, where collisional effects are sufficiently large. Furthermore, the threshold conditions are almost independent of the wave propagation direction and the unstable waves propagate at speeds well below the ion acoustic speed. The role that this new instability may play in recent radar backscatter observations of short scale irregularities propagating in directions close to that of the geomagnetic field, in the lower E-region is also considered.     

Effects of substorms on the stormtime ring current index Dst

There has been some discussion in recent times regarding whether or not substorm expansive phase activity plays any role of importance in the formation of the stormtime ring current. I explore this question using the K_p index as a proxy for substorm expansive phase activity and the D_st index as a proxy for symmetric ring current strength. I find that increases in D_st are mildly related to the strength of substorm expansive phase activity during the development of the storm main phase. More surprisingly, I find that the strength of D_st during the storm recovery phase is positively correlated with the strength of substorm expansive phase activity. This result has an important bearing on the question of how much the Dst index reflects activity other than that of the stormtime symmetric ring current strength for which it is supposed to be a proxy.