Modeling the responses of the middle latitude ionosphere to solar flares

In this paper, we investigate the solar flare effects of the ionosphere at middle latitude with a one-dimensional ionosphere theoretical model. The measurements of solar irradiance from the SOHO/Solar EUV Monitor (SEM) and GOES satellites have been used to construct a simple time-dependent solar flare spectrum model, which serves as the irradiance spectrum during solar flares. The model calculations show that the ionospheric responses to solar flares are largely related to the solar zenith angle. During the daytime most of the relative increases in electron density occur at an altitude lower than 300 km, with a peak at about 115 km, whereas around sunrise and sunset the strongest ionospheric responses occur at much higher altitudes (e.g. 210 km for a summer flare). The ionospheric responses to flares in equinox and winter show an obvious asymmetry to local midday with a relative increase in total electron content (TEC) in the morning larger than that in the afternoon. The flare-induced TEC enhancement increases slowly around sunrise and reaches a peak at about 60 min after the flare onset.     

Enhancement of OI 630.0 nm emission at mid-latitudes during an intense magnetic storm

All-sky camera (ASC), Global Positioning System (GPS), and ionosonde measurements were used to investigate the upper atmospheric variations at mid-latitude during the strong geomagnetic storm on October 29–31, 2003. An arc-shaped 630.0 nm emission was observed in the northern sky on all-sky images taken at Mt. Bohyun (36.2°N, 128.9°E, GMLAT=29°N) in Korea during 17:48–8:58 UT in the main phase of the geomagnetic storm on October 29. The NmF₂ and hmF₂ from the ionosonde show strong disturbances at that time. The vertical profiles of electron densities, calculated by the ionospheric tomographic method using ground-based GPS slant total electron contents measurements, show the largest value at ∼440 km height at 18:30 UT on October 29 when the enhancements of OI 630.0 nm emission were observed. The arc-shaped red emission observed during the main phase of the magnetic storm is likely a low-latitude red aurora due to its short duration of ∼1 h. The result implies that the plasmapause was at L=1.4–1.6 during the geomagnetic storm. The fact that the arc did not follow a constant L-value appears to suggest that neutral precipitation and a traveling ionospheric disturbance could also be the cause of the arc.     

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.     

The effects of thunderstorm-generated atmospheric gravity waves on mid-latitude F-region drifts

Superposed epoch analysis (SEA) was used to examine ionospheric drift velocities measured by a digital ionosonde located at the mid-latitude station Bundoora (145.1°E, 37.7°S geographic), near Melbourne. The control times for the SEA were the times of cloud-to-ground (CG) lightning strokes measured from August 2003 to August 2004 by the World Wide Lightning Location Network (WWLLN). Statistically, regions of concentrated lightning activity migrated from west to east across Bundoora, and the stroke frequency was higher the day prior the activity reached the station, and lower on the day after it passed to the east. For the SEA, CG strokes were separated into four directional quadrants centred on north, south, east and west. No SEA results are shown for the south quadrant due to the relatively low detection frequency of strokes across the Southern Ocean (6% of all events). The strongest downward vertical perturbations in F-region drifts, ?4.5 m s^?1, were found for lightning located towards the west during ?30 to ?16 h (i.e., the afternoon prior the activity passed near the station at t=0 h). The downward perturbation decreased in amplitude to ?1.5 m s^?1 for lightning located towards the north during ?6–+6 h, and was weakest (?0.7 m s^?1) for lightning located towards the east during +16–+28 h (i.e., the next afternoon). There were directionally consistent perturbations in the drift azimuths associated with the lightning located in their respective quadrants; lightning located to the west of the station caused eastward azimuth enhancements, northward lightning caused southward enhancements, and eastward lightning caused westward enhancements. Velocity magnitudes and fluctuations tended to increase during the passage of lightning. The observed responses were stronger when the SEA was performed with data selected using time windows of <2 min on either side of each lightning stroke. However, they persisted at longer time scales and were strong when thunderstorm onsets (instead of lightning times) were used as controls. Our results can be explained by thunderstorm-generated atmospheric gravity waves (AGWs) which subsequently gave rise to medium-scale travelling ionospheric disturbances (MSTIDs), with the lightning strokes acting merely as a proxy for this coupling. The prevailing thermospheric winds were flowing from east to west across the study region, and may have acted as a directional ‘filter’ for the MSTIDs, allowing waves generated in the west quadrant to reach the station and preventing those generated in other quadrants. Displacement of the MSTIDs in the direction anti-parallel to mean neutral wind flow has been observed by (Waldock, J.A., Jones, T.B., 1986. HF Doppler observations of medium-scale travelling ionospheric disturbances at mid-latitudes. Journal of atmospheric and terrestrial physics 48(3), 245–260).     

Space weather in the thermospheric–ionospheric domain over the Brazilian region: Climatology of ionospheric plasma bubbles in the subequatorial and low-latitude region

The climatology of ionospheric plasma bubbles is studied here by means of a comparison of the frequency of occurrence of the spread-F/plasma bubble events over the South American region using the images from two OI 630 nm imager systems located at the subequatorial station São João do Cariri—CA (7.4 S, 36.5 W, 20 S dip) and the low-latitude station Cachoeira Paulista—CP (22.5 S, 45 W, 33 S dip) in Brazil during the years of 2004 and 2005. The results are discussed in the light of current theory and geomagnetic parameters of the two observation stations.     

The coherence between the IMF and high-latitude ionospheric flows: The dayside magnetosphere–ionosphere low-pass filter

This study seeks to establish a new system characteristic describing dayside convective flows in the coupled magnetosphere–ionosphere: the low-pass filter function through which interplanetary magnetic field (IMF) fluctuations are processed as they are communicated from the magnetopause to the high-latitude ionosphere near local noon. In doing so, this study confirms that variations in the ionospheric flows at high-latitudes near local noon are well correlated with variations in the IMF orientation and magnitude on short timescales. We construct the filter function by comparing time series of the ionospheric equivalent flows at a fixed location at magnetic local noon and 80° latitude with time series of the IMF. The coherence spectra of these two parameters—averaged over 330 h of comparison—indicate that there is a low-pass cutoff in the ionospheric response to IMF driving at a periods shorter than 20 min (frequencies higher than 0.8 mHz). When there is sufficient power in the IMF fluctuations, this cutoff is relatively sharp—the coherence drops by roughly a factor of three between the periods 32 and 21 min (0.5 and 0.8 mHz). The results also show that on average the coherence between the east–west component of the equivalent flows and IMF B_y tends to be less than the coherence between the north–south component of the equivalent flows and IMF B_z.     

On the connection between the atmospheric electric field measured at the surface and the ionospheric electric field in the Central Antarctica

Regular measurements of the atmospheric electric field made at Vostok Station (φ=78.45°S; λ=106.87°E, elevation 3500 m) in Antarctica demonstrate that extremely intense electric fields (1000–5000 V/m) can be observed during snow storms. Usually the measured value of the atmospheric electric field at Vostok is about 100–250 V/m during periods with “fair weather” conditions. Actual relation between near-surface electric fields and ionospheric electric fields remain to be a controversial problem. Some people claimed that these intense electric fields produced by snowstorms or appearing before strong earthquakes can re-distribute electric potential in the ionosphere at the heights up to 300 km. We investigated interrelation between the atmospheric and ionospheric electric fields by both experimental and theoretical methods. Our conclusion is that increased near-surface atmospheric electric fields do not contribute notably to distribution of ionospheric electric potential.     

Ionospheric and geomagnetic disturbances during the 2005 Sumatran earthquakes

This paper investigates the ionospheric and geomagnetic responses during the 28 March 2005 and 14 May 2005 Sumatran earthquakes using GPS and magnetometer stations located in the near zone of the epicenters. These events occurred during low solar and geomagnetic activity. TEC oscillations with periods of 5–10 min were observed about 10–24 min after the earthquakes and have horizontal propagation velocities of 922–1259 m/s. Ionospheric disturbances were observed at GPS stations located to the northeast of the epicenters, while no significant disturbances were seen relatively east and south of the epicenters. The magnetic field measurements show rapid fluctuations of 4–5 s shortly after the earthquake, followed by a Pc5 pulsation of 4.8 min about 11 min after the event. The correlation between the ionospheric and geomagnetic responses shows a good agreement in the period and time lag of the peak disturbance arrival, i.e. about 11–13 min after the earthquake.     

Polar mesosphere summer echoes (PMSE) studied at Bragg wavelengths of 2.8 m, 67 cm,and 16 cm

We present observations of radar volume reflectivities under conditions of polar mesosphere summer echoes (PMSE) at three frequencies, i.e., 53.5, 224, and 930 MHz corresponding to Bragg wavelengths of 2.8, 0.67, and 0.16 m. These measurements were made with the ALWIN radar in Andenes and the EISCAT VHF and UHF radars in Tromsø. Contributions to the signal at 930 MHz by incoherent scatter are used to estimate electron number densities and their gradient at PMSE altitudes, and spectral width measurements of Doppler spectra recorded at 224 MHz are used to estimate the turbulent energy dissipation rate. We further derive a theoretical expression for the radar volume reflectivity for the case of turbulent scatter aided by a large Schmidt number (i.e., the current standard theory of PMSE) and show that our observations quantitatively agree with this theory if Schmidt numbers between 2500 and 5000 are assumed. We then show that these Schmidt numbers correspond to ice particles with radii in the range 20–30 nm which should frequently occur in the polar summer mesopause region. In addition, we show that for the short period when PMSE was observed at UHF frequencies the volume reflectivity is proportional to a factor determined by the turbulent energy dissipation rate, electron number density, and the electron number density gradient in agreement with theory. We consider our findings as strong support that PMSE at all considered frequencies is indeed created by turbulent scatter in the presence of a large Schmidt number. We finally highlight that ultimate proof of this concept will require the direct measurement of ice particle sizes in a PMSE environment probed by radars covering frequencies between 50 MHz and 1 GHz.     

The wave-like structures of ionospheric perturbation associated with Sumatra earthquake of 26 December 2004, as revealed from VLF observation in Japan of NWC signals

The observations of subionospheric VLF waves from the Australian VLF transmitter NWC (frequency=19.8 kHz) at the Japanese receiving stations Chofu, Chiba and Kochi have been utilized to identify a possible precursor of ionospheric perturbations to the huge Sumatra earthquake of 26 December 2004. The VLF amplitude data at Japanese stations have indicated the depression in amplitude and also the enhancement in nighttime amplitude fluctuation before the earthquake. The nighttime fluctuation is composed of wave-like structures, and the wavelet analysis and cross-correlation analyses have been performed for those fluctuations. A significant enhancement in the fluctuation spectra in the period 20–30 min to ∼100 min (the frequency range of atmospheric gravity waves) is observed only before the earthquake. Then, the wave-like structures tend to propagate from the NWC–Kochi path to NWC–Chiba path with the time delay of ∼2 h, and so the wave propagation speed is estimated as ∼20 m/s. This finding might be important when we think of lithosphere–ionosphere coupling mechanism.     

Prompt derivation of TEC from GEONET data for space weather monitoring in Japan

Continuous monitoring of ionospheric conditions is essential to monitoring and forecasting space weather. The worldwide use of global navigation satellite systems like the Gobal Positioning System (GPS) makes it possible to continuously monitor the total electron content (TEC) of the ionosphere and plasmasphere up to a height of about 20,000 km. We have developed a system for deriving the TEC from GEONET data rapidly and we use the TEC distribution over Japan in the daily operations of the Space Weather Forecast Center at NICT (RWC Tokyo of ISES). Using instrumental biases from a few days before enables us to drastically shorten the processing time for deriving TEC. The latest TEC values (with a delay of about 1 h) are obtained every 3 h, and most of the values are within 2 TEC units of the actual TEC. We have found our system for deriving TEC rapidly to be useful for continuously monitoring the progress of ionospheric storms under any ionospheric conditions, even those under which the usual ionosonde observations are unable to obtain F-region profiles.     

A novel method for the quantitative assessment of the ionosphere effect on high accuracy GNSS applications,which require ambiguity resolution

Real time kinematic, or RTK, is a high-accuracy GPS relative positioning technique, which allows to measure positions in real time with an accuracy usually better than 1 decimeter. Ionospheric small-scale variability can strongly degrade RTK accuracy. In this paper, we present a method allowing to assess in a direct quantitative way the influence of the ionospheric activity on RTK accuracy. We apply this method to two different ionospheric situations: a day where strong travelling ionospheric disturbances (TIDs) were detected (December 24, 2004) and a day where a severe geomagnetic storm was observed (November 20, 2003). We show that on a 4 km baseline, strong TIDs have the same influence as the ionospheric variability induced by a geomagnetic storm on RTK accuracy: in both cases errors of more than 1.5 m are observed.     

Ionospheric activity and possible connection with seismicity: Contribution from the analysis of long time series of GNSS signals

The modifications of some atmospheric physical properties prior to a high magnitude earthquake were debated in the frame of the Lithosphere Atmosphere Ionosphere Coupling (LAIC) model. In this work, among the variety of involved phenomena, the ionisation of air at the ionospheric levels triggered by the leaking of gases from the Earth’s crust was investigated through the analysis of GNSS (Global Navigation Satellite System) signals. In particular, the authors analysed a 5 year (2008–2012) long series of GNSS based ionospheric TEC to produce maps over an area surrounding the epicentre of the L’Aquila (Italy, M_w = 6.3) earthquake of April 6th, 2009. The series was used to detect and quantify amplitude and duration of episodes of ionospheric disturbances by a statistical approach and to discriminate local and global effects on the ionosphere comparing these series with TEC values provided by the analysis of GNSS data from international permanent trackers distributed over a wider region. The study found that during this time interval only three statistically meaningful episodes of ionospheric disturbances were observed. One of them, occurring during the night of 16th of March 2009, anticipated the main shock by 3 weeks and could be connected with the strong earthquake of 6th of April. The other two significant episodes were detected within periods that were not close to the main seismic events and are more likely due to various and global reasons.     

Detection of large scale TIDs associated with the dayside cusp using SuperDARN data

Variations in the dayside ionosphere parameters as a result of a large-scale acoustic gravity wave (LS AGW) were studied for the 17 February 1998 substorm using the super dual auroral radar network (SuperDARN) measurements. This event was characterised by a sharp rise in the AE index with a maximum of ～900 nT. The source of the disturbance responsible for the LS AGW appears to have been located within the plasma convection throat and in the dayside cusp region. The location of the source was obtained from studies of a number of datasets including high-latitude convection maps, data from 4 DMSP satellites and networks of ground-based magnetometers. The propagation of the LS AGWs caused quasi-periodic variations in the skip distance (with an amplitude up to 220–260 km) of the ground backscatter measured by up to 6 SuperDARN radars, including Goose Bay and Kapuskasing, resulting in two large-scale travelling ionospheric disturbances (LS TIDs). The LS TIDs had wave periods of 1.5 and 2 h, a velocity of ～400 m/s for both, and wavelengths of 2200 and 2900 km, respectively. These quasi-periodic variations were also present in the peak electron density and height of the F2 layer measured by the Goose Bay ionosonde. The numerical simulation of the inverse problem show good agreement between Goose Bay radar and Goose Bay ionosonde measurements. But these LS TIDs would be difficult to deduce from the ground based ionospheric station data alone, because h_mF2 variations were 10–40 km only and f_OF2 variations between 10% and 20%. The results demonstrate how important SuperDARN radars can be, and that this is a more powerful technique than routine ground-based sounding for studies of weak quasi-periodic variations in the dayside subauroral ionosphere related to LS AGW.     

Forsterite to wadsleyite phase transformation under shear stress and consequences for the Earth's mantle transition zone

We have studied the phase transformation of forsterite to wadsleyite under shear stress at the Earth's transition zone pressure and temperature conditions. Two-step experiments were performed using a multi-anvil press. First, we hot pressed iron-free forsterite at 6 or 11 GPa and 1100 °C. Then we deformed a slab of this starting material using a direct simple shear assembly at 16 GPa and 1400 °C for 1, 15, 35, 40, or 60 min. Both the starting material and the deformed samples were characterized using optical and scanning electron microscopy including measurements of crystal preferred orientations (CPO) by electron back scattered diffraction (EBSD), transmission electron microscopy (TEM), and Fourier transform infrared spectroscopy (FTIR). The phase transformation occurs very rapidly, in less than 1 min, and metastable forsterite relics are not observed after deformation. The grain size of wadsleyite is slightly smaller than the forsterite starting material. The water contents obtained from FTIR analyses in forsterite and wadsleyite are 65–124 wt ppm H₂O and 114–736 wt ppm H₂O, respectively, which are well below water solubility at similar conditions in the presence of free water. Wadsleyite aggregates display weak CPO patterns with [1 0 0] axes concentrated at low angle to the shear direction, [0 1 0] axes perpendicular to the shear plane and nearly random [0 0 1] axes. Only a few dislocations were observed in wadsleyite with TEM. This observation is consistent with the assumption that most dislocations formed during the initial high-stress stages of these stress-relaxation experiments, were consumed in the phase transformation, probably enhancing the transformation rate. CPO patterns vary as a function of the water content: with increasing water content the density of [1 0 0] axes parallel to the shear direction decreases, and the density of [0 0 1] axes increases. Viscoplastic self-consistent modeling of CPO evolution using previously reported glide systems for wadsleyite, i.e., [1 0 0]{0 k l} and 1/2 〈1 1 1〉{1 0 1}, cannot reproduce the measured CPO, unless the [0 0 1](0 1 0) system, for which dislocations have not been observed by TEM, is also activated. In addition, wadsleyite grain growth suggests the participation of diffusion-assisted processes in deformation. Calculated anisotropies for P and S-waves using measured CPO are always below 1%. This very low anisotropy is due to both the low finite strain achieved in the experiments, which leads to weak wadsleyite CPO, and to the diluting effect of added majorite. The present experiments emphasize the importance of stress, grain size evolution and water content in the forsterite to wadsleyite phase transformation and subsequent deformation in the transition zone.     

First 10 min intervals of Pi2 onset at geosynchronous altitudes during the expansion of energetic ion regions in the nighttime sector

We examined the temporal variations of the geomagnetic field and energetic ions at geosynchronous altitudes associated with substorms during the nighttime using a superposed epoch analysis timed by Pi2 onset. We focused on the first 10 min intervals of Pi2 onset and on subsequent intervals to study the substorm expansion. We conclude that the first 10 min interval of Pi2 onset is a transitional state of the substorm dominated by MHD processes associated with earthward flow and its bifurcation. Intervals of field line variations following the first 10 min were well organized by dipolarization (substorm current wedge) due to the reduced cross-tail current. We also show that energetic ion regions localized in the local time sector from 2000 to 0000 LT in the first 10 min intervals of Pi2 onset expanded to the post-midnight sector, reaching 0400 LT within 20 min after Pi2 onset. We conclude that the expansion of the energetic plasma regions can be attributed to the inflation of the inner magnetosphere during dipolarization.     

Spectacular low- and mid-latitude electrical fields and neutral winds during a superstorm

In November 2004, a major magnetic storm occurred, a lengthy portion of which was recorded by the Upper Atmospheric Radar Chain. On the 9th and 10th, the Jicamarca Radar detected the highest magnitude penetrating electric fields (±3 mV/m) and vertical drifts (±120 m/s) ever seen at this premiere facility. These large and variable drifts were highly correlated with the interplanetary magnetic and electric fields and created a double F layer on the dayside and unusual TEC behavior throughout the low-latitude zone. These solar wind-induced drifts both suppressed and generated irregularities at the magnetic equator at different times. Large-scale thermospheric disturbances were generated by high-latitude heating and tracked through the middle- to low-latitude zones where both parallel and perpendicular plasma drifts created major ionospheric changes. The auroral oval was located at a magnetic L shell of about three for many hours.     

TEC variations over the Mediterranean before and during the strong earthquake (M = 6.5) of 12th October 2013 in Crete,Greece

In this paper, the total electron content (TEC) data from eight global positioning system (GPS) stations of the EUREF network, provided by IONOLAB (Turkey), were analyzed using discrete Fourier analysis to investigate the TEC variations over the Mediterranean before and during the strong earthquake of 12th October 2013, which occurred west of Crete, Greece. In accordance with the results of similar analyses in the area, the main conclusions of this study are the following: (a) TEC oscillations in a broad range of frequencies occur randomly over an area of several hundred km from the earthquake and (b) high frequency oscillations (f ⩾ 0.0003 Hz, periods T ⩽ 60 m) may point to the location of the earthquake with questionable accuracy. The fractal characteristics of the frequency distribution may point to the locus of the earthquake with higher accuracy. We conclude that the lithosphere–atmosphere–ionosphere coupling (LAIC) mechanism through acoustic or gravity waves could explain this phenomenology.     

Seasonal and diurnal variation of electron and iron concentrations at the meteor heights above Arecibo

We report observations of seasonal and local time variation of the averaged electron and iron concentrations, as well as simultaneous measurements of the two species, above the Arecibo Observatory (18.35°N, 66.75°N), Puerto Rico. The average Fe profile between 21:00 and 24:00 LT has a single peak at about 85 km with the exception of the summer when an additional peak exists at about 95 km. The higher Fe peak in the summer is correlated with higher electron concentrations in this season. The three nights of simultaneous measurements of electron and iron concentrations show that narrow layers of Fe and electrons are well correlated. Comparison of the climatological and simultaneous Fe and electron data suggests that recombination of Fe⁺ plays an important role in determining the Fe profile in the upper part of the Fe layer. Above 93 km, the Fe concentration appears to increase after sunset if the electron concentration exceeds about 4000 electrons cm⁻³. The average rate of Fe production is about 0.1 atom cm⁻³ s⁻¹ for all seasons at 100 km in the early evening hours. A chemical model reveals that the concentration of Fe⁺ must be 50–80% of the total ionization over Arecibo for typical equinox conditions to explain the observed rate of Fe production. These high relative Fe⁺ concentrations are consistent with in situ observations that Fe⁺ is usually the dominant ion in sporadic E layers in the nighttime lower E region. This suggests that the source of Fe⁺ is provided by sporadic E layers descending over Arecibo after sunset. The Fe density between 80 and 85 km decreases during the night, for all seasons. This is attributed to the formation of stable molecular Fe species, such as FeOH, due to the increase in O₃ and decrease in atomic O and H during the night at these altitudes.     

Storm time phase space density radial profiles of energetic electrons for small and large K values: SCATHA results

Measurements of the radial gradient of the phase space density (PSD) at constant first and second invariants provide a test of whether storm time electron acceleration processes are dominated by inward radial transport of electrons or whether other processes must be considered. We used the detailed energetic electron angular and spectral data from the SCATHA satellite to determine the evolution of electron PSD radial profiles through a moderate magnetic storm (2 May 1986; D_st=−95 nT). We compared the changes that occur in the profiles in the storm recovery period to the pre-storm profiles for first invariant values M=200–2500 MeV/G and for K values of 0.055–0.65 Re√G over the L* range of ∼5.2–7.3 Re. The PSD radial profiles showed a range of features from peaked in L* at small K and M during the pre-storm period to those in the late storm recovery phase, that were flat, had negative slopes for small K and had peaks in the L*=5.2–6.5 range for intermediate to large K. There were significant differences in the radial profiles for small and intermediate K values at constant M. For example, during the recovery period the PSD profiles were flat or decreasing for K=0.06 and M>1200 MeV/G, while they were peaked near L*=5.75 for K>0.2. These results imply that radial diffusion is a reasonable explanation of the near-equatorial post-storm PSD enhancements for L*>5.2 for this storm but that either significant electron pitch angle transport, losses, and/or acceleration of off-equatorial mirroring electrons by waves play an important role in the evolution off-equator PSD profiles during the storm recovery.