Practical simulation of GNSS signals in the presence of ionospheric scintillaton

We propose a low complexity technique to generate amplitude correlated time-series with Nakagami-m distribution and phase correlated Gaussian-distributed time-series, which is useful in the simulation of ionospheric scintillation effects during the transmission of GNSS signals. The method requires only the knowledge of parameters S4 (scintillation index) and σΦ (phase standard deviation) besides the definition of models for the amplitude and phase power spectra. The Zhang algorithm is used to produce Nakagami-distributed signals from a set of Gaussian autoregressive processes.

An enhancement of the ionospheric sporadic‐E layer in response to negative polarity cloud‐to‐ground lightning

Lightning data, collected using a Boltek Storm Tracker system installed at Chilton, UK, were used to investigate the mean response of the ionospheric sporadic-E layer to lightning strokes in a superposed epoch study. The lightning detector can discriminate between positive and negative lightning strokes and between cloud-to-ground ( CG) and inter-cloud ( IC) lightning. Superposed epoch studies carried out separately using these subsets of lightning strokes as trigger events have revealed that the dominant cause of the observed ionospheric enhancement in the Es layer is negative cloud-to-ground lightning.

The location of lightning affecting the ionospheric sporadic-E layer as evidence for multiple enhancement mechanisms

We present a study of the geographic location of lightning affecting the ionospheric sporadic-E (Es) layer over the ionospheric monitoring station at Chilton, UK. Data from the UK Met Office's Arrival Time Difference (ATD) lightning detection system were used to locate lightning strokes in the vicinity of the ionospheric monitoring station. A superposed epoch study of this data has previously revealed an enhancement in the Es layer caused by lightning within 200km of Chilton. In the current paper, we use the same data to investigate the location of the lightning strokes which have the largest effect on the Es layer above Chilton. We find that there are several locations where the effect of lightning on the ionosphere is most significant statistically, each producing different ionospheric responses. We interpret this as evidence that there is more than one mechanism combining to produce the previously observed enhancement in the ionosphere.

Lightning-induced intensification of the ionospheric sporadic E layer

A connection between thunderstorms and the ionosphere has been hypothesized since the mid-1920s(1). Several mechanisms have been proposed to explain this connection(2-7), and evidence from modelling(8) as well as various types of measurements(9-14) demonstrate that lightning can interact with the lower ionosphere. It has been proposed, on the basis of a few observed events(15), that the ionospheric 'sporadic E' layer - transient, localized patches of relatively high electron density in the mid-ionosphere E layer, which significantly affect radio-wave propagation - can be modulated by thunderstorms, but a more formal statistical analysis is still needed. Here we identify a statistically significant intensification and descent in altitude of the mid-latitude sporadic E layer directly above thunderstorms. Because no ionospheric response to low-pressure systems without lightning is detected, we conclude that this localized intensification of the sporadic E layer can be attributed to lightning. We suggest that the co-location of lightning and ionospheric enhancement can be explained by either vertically propagating gravity waves that transfer energy from the site of lightning into the ionosphere, or vertical electrical discharge, or by a combination of these two mechanisms.

Timing and localization of ionospheric signatures associated with substorm expansion phase onset

In this paper, we present case studies of the optical and magnetic signatures of the characteristics of the first minute of substorm expansion phase onset observed in the ionosphere. We find that for two isolated substorms, the onset of magnetic pulsations in the 24–96 s period wavelet band are colocated in time and space with the formation and development of small-scale optical undulations along the most equatorward preexisting auroral arc prior to auroral breakup. These undulations undergo an inverse spatial cascade into vortices prior to the release of the westward traveling surge. We also present a case study of a multiple activation substorm, whereby discrete onsets of ULF wave power above a predetermined quiet time threshold are shown to be associated with specific optical intensifications and brightenings. Moreover, in the multiple activation substorm event, we show that neither the formation of the small-scale undulations nor the formation of similar structures along a north–south aligned arc is sufficient to produce auroral breakup associated with expansion phase onset. It is only ∼10 min after these two disparate activation regions initiate that auroral breakup and the subsequent formation of a westward traveling surge occur. We discuss the implications of these results in terms of the triggering mechanisms likely to be occurring during these specific events.

A numerical model of the ionospheric signatures of time-varying magnetic reconnection: III. Quasi-instantaneous convection responses in the Cowley-Lockwood paradigm

Using a numerical implementation of the Cowley and Lockwood (1992) model of flow excitation in the magnetosphere–ionosphere (MI) system, we show that both an expanding (on a _12-min timescale) and a quasiinstantaneous response in ionospheric convection to the onset of magnetopause reconnection can be accommodated by the Cowley–Lockwood conceptual framework. This model has a key feature of time dependence, necessarily considering the history of the coupled MI system. We show that a residual flow, driven by prior magnetopause reconnection, can produce a quasi-instantaneous global ionospheric convection response; perturbations from an equilibrium state may also be present from tail reconnection, which will superpose constructively to give a similar effect. On the other hand, when the MI system is relatively free of pre-existing flow, we can most clearly see the expanding nature of the response. As the open-closed field line boundary will frequently be in motion from such prior reconnection (both at the dayside magnetopause and in the cross-tail current sheet), it is expected that there will usually be some level of combined response to dayside reconnection.

Modeling the observed proton aurora and ionospheric convection responses to changes in the IMF clock angle: 2. Persistence of ionospheric convection

We apply a numerical model of time-dependent ionospheric convection to two directly driven reconnection pulses during a 15-min interval of southward IMF on 26 November 2000. The model requires an input magnetopause reconnection rate variation, which is here derived from the observed variation in the upstream IMF clock angle, q. The reconnection rate is mapped to an ionospheric merging gap, the MLT extent of which is inferred from the Doppler-shifted Lyman-a emission on newly opened field lines, as observed by the FUV instrument on the IMAGE spacecraft. The model is used to reproduce a variety of features observed during this event: SuperDARN observations of the ionospheric convection pattern and transpolar voltage; FUV observations of the growth of patches of newly opened flux; FUVand in situ observations of the location of the Open-Closed field line Boundary (OCB) and a cusp ion step. We adopt a clock angle dependence of the magnetopause reconnection electric field, mapped to the ionosphere, of the form Enosin4(q/2) and estimate the peak value, Eno, by matching observed and modeled variations of both the latitude, LOCB, of the dayside OCB (as inferred from the equatorward edge of cusp proton emissions seen by FUV) and the transpolar voltage FPC (as derived using the mapped potential technique from SuperDARN HF radar data). This analysis also yields the time constant tOCB with which the open-closed boundary relaxes back toward its equilibrium configuration. For the case studied here, we find tOCB = 9.7 ± 1.3 min, consistent with previous inferences from the observed response of ionospheric flow to southward turnings of the IMF. The analysis confirms quantitatively the concepts of ionospheric flow excitation on which the model is based and explains some otherwise anomalous features of the cusp precipitation morphology.

Modeling the observed proton aurora and ionospheric convection responses to changes in the IMF clock angle: 1. Persistence of cusp proton aurora

We employ a numerical model of cusp ion precipitation and proton aurora emission to fit variations of the peak Doppler-shifted Lyman-a intensity observed on 26 November 2000 by the SI-12 channel of the FUV instrument on the IMAGE satellite. The major features of this event appeared in response to two brief swings of the interplanetary magnetic field (IMF) toward a southward orientation. We reproduce the observed spatial distributions of this emission on newly opened field lines by combining the proton emission model with a model of the response of ionospheric convection. The simulations are based on the observed variations of the solar wind proton temperature and concentration and the interplanetary magnetic field clock angle. They also allow for the efficiency, sampling rate, integration time and spatial resolution of the FUV instrument. The good match (correlation coefficient 0.91, significant at the 98% level) between observed and modeled variations confirms the time constant (about 4 min) for the rise and decay of the proton emissions predicted by the model for southward IMF conditions. The implications for the detection of pulsed magnetopause reconnection using proton aurora are discussed for a range of interplanetary conditions.

A numerical model of the ionospheric signatures of time-varying magnetic reconnection: II. Measuring expansions in the ionospheric flow response

A numerical model embodying the concepts of the Cowley-Lockwood (Cowley and Lockwood, 1992, 1997) paradigm has been used to produce a simple Cowley– Lockwood type expanding flow pattern and to calculate the resulting change in ion temperature. Cross-correlation, fixed threshold analysis and threshold relative to peak are used to determine the phase speed of the change in convection pattern, in response to a change in applied reconnection. Each of these methods fails to fully recover the expansion of the onset of the convection response that is inherent in the simulations. The results of this study indicate that any expansion of the convection pattern will be best observed in time-series data using a threshold which is a fixed fraction of the peak response. We show that these methods used to determine the expansion velocity can be used to discriminate between the two main models for the convection response to a change in reconnection.

A numerical model of the ionospheric signatures of time-varying magneticreconnection: I. ionospheric convection

This paper presents a numerical model for predicting the evolution of the pattern of ionospheric convection in response to general time-dependent magnetic reconnection at the dayside magnetopause and in the cross-tail current sheet of the geomagnetic tail. The model quantifies the concepts of ionospheric flow excitation by Cowley and Lockwood (1992), assuming a uniform spatial distribution of ionospheric conductivity. The model is demonstrated using an example in which travelling reconnection pulses commence near noon and then move across the dayside magnetopause towards both dawn and dusk. Two such pulses, 8 min apart, are used and each causes the reconnection to be active for 1 min at every MLT that they pass over. This example demonstrates how the convection response to a given change in the interplanetary magnetic field (via the reconnection rate) depends on the previous reconnection history. The causes of this effect are explained. The inherent assumptions and the potential applications of the model are discussed.

Ionospheric ion and electron heating at the poleward boundary of a poleward expanding substorm-disturbed region

We present observations of a poleward propagating substorm-disturbed region which was observed by the European Incoherent SCATter (EISCAT) radar and the Svalbard International Monitor for Auroral Geomagnetic Effects (IMAGE) magnetometers in the postmidnight sector. The expansion of the disturbance was launched by a substorm intensification which started similar to 25 min after the initial onset, and similar to 10 min before the disturbance arrived over Svalbard. In association with the magnetic disturbance, a poleward expanding enduring enhancement in the F region electron temperature was observed, indicative of soft electron precipitation, with a narrow band of enhanced ion temperature straddling its poleward edge, indicative of fast ion flows and ion-neutral collisional heating. This electron temperature boundary was coincident with the poleward propagating electrojet current system detected by the high-latitude IMAGE magnetometer stations and is taken to be a proxy for the observation of a substorm auroral bulge. The electron temperature boundary is inferred to have a width comparable or less than one radar range gate (similar to 60 km transverse to the magnetic field), while the region of high ion temperature was found to be approximately three gates wide, extending approximately two gates (similar to 120 km) poleward of the electron temperature boundary, and approximately one gate (similar to 60 km) equatorward. The two-beam radar line-of-sight velocity data are found to be consistent with the existence of a layer of high-speed flow in the boundary, peaking at values similar to1.5-3 km s(-1), roughly consistent with the ion temperature data. The flow is directed either east or west along the boundary depending on the direction of the flow in the poleward region. We infer that the flow is deflected along and around the boundary of the substorm-disturbed region due to the high conductivity of the latter. Variations in the flow poleward of the boundary produced no discernible magnetic effects on the ground, confirming the low conductivity of the preboundary ionosphere.

Ionospheric measurements of relative coronal brightness during the total solar eclipses of 11 August, 1999 and 9 July, 1945

Swept-frequency (1-10 MHz) ionosonde measurements were made at Helston, Cornwall (50 degrees 06'N, 5 degrees 18'W) during the total solar eclipse on August 11, 1999. Soundings were made every three minutes. We present a method for estimating the percentage of the ionising solar radiation which remains unobscured at any time during the eclipse by comparing the variation of the ionospheric E-layer with the behaviour of the layer during a control day. Application to the ionosonde date for II August, 1999, shows that the flux of solar ionising radiation fell to a minimum of 25 +/- 2% of the value before and after the eclipse. For comparison, the same technique was also applied to measurements made during the total solar eclipse of 9 July, 1945, at Sormjole (63 degrees 68'N, 20 degrees 20'E) and yielded a corresponding minimum of 16 +/- 2%. Therefore the method can detect variations in the fraction of solar emissions that originate from the unobscured corona and chromosphere. We discuss the differences between these two eclipses in terms of the nature of the eclipse, short-term fluctuations, the sunspot cycle and the recently-discovered long-term change in the coronal magnetic field.

The first real-time worldwide ionospheric predictions network: an advance in support of spaceborne experimentation, on-line model validation, and space weather

We report on the first realtime ionospheric predictions network and its capabilities to ingest a global database and forecast F-layer characteristics and "in situ" electron densities along the track of an orbiting spacecraft. A global network of ionosonde stations reported around-the-clock observations of F-region heights and densities, and an on-line library of models provided forecasting capabilities. Each model was tested against the incoming data; relative accuracies were intercompared to determine the best overall fit to the prevailing conditions; and the best-fit model was used to predict ionospheric conditions on an orbit-to-orbit basis for the 12-hour period following a twice-daily model test and validation procedure. It was found that the best-fit model often provided averaged (i.e., climatologically-based) accuracies better than 5% in predicting the heights and critical frequencies of the F-region peaks in the latitudinal domain of the TSS-1R flight path. There was a sharp contrast however, in model-measurement comparisons involving predictions of actual, unaveraged, along-track densities at the 295 km orbital altitude of TSS-1R In this case, extrema in the first-principle models varied by as much as an order of magnitude in density predictions, and the best-fit models were found to disagree with the "in situ" observations of Ne by as much as 140%. The discrepancies are interpreted as a manifestation of difficulties in accurately and self-consistently modeling the external controls of solar and magnetospheric inputs and the spatial and temporal variabilities in electric fields, thermospheric winds, plasmaspheric fluxes, and chemistry.

Ionospheric and geomagnetic responses to changes in IMF B_Z: a superposed epoch study

Superposed epoch studies have been carried out in order to determine the ionospheric response at mid-latitudes to southward turnings of the interplanetary magnetic field (IMF). This is compared with the geomagnetic response, as seen in the indices K p, AE and Dst. The solar wind, IMF and geomagnetic data used were hourly averages from the years 1967–1989 and thus cover a full 22-year cycle in the solar magnetic field. These data were divided into subsets, determined by the magnitudes of the southward turnings and the concomitant increase in solar wind pressure. The superposed epoch studies were carried out using the time of the southward turning as time zero. The response of the mid-latitude ionosphere is studied by looking at the F-layer critical frequencies, f o F2, from hourly soundings by the Slough ionosonde and their deviation from the monthly median values, δf o F2. For the southward turnings with a change in B z of δB z > 11.5 nT accompanied by a solar wind dynamic pressure P exceeding 5 nPa, the F region critical frequency, f o F2, shows a marked decrease, reaching a minimum value about 20 h after the southward turning. This recovers to pre-event values over the subsequent 24 h, on average. The Dst index shows the classic storm-time decrease to about −60 nT. Four days later, the index has still to fully recover and is at about −25 nT. Both the K p and AE indices show rises before the southward turnings, when the IMF is strongly northward but the solar wind dynamic pressure is enhanced. The average AE index does register a clear isolated pulse (averaging 650 nT for 2 h, compared with a background peak level of near 450 nT at these times) showing enhanced energy deposition at high latitudes in substorms but, like K p, remains somewhat enhanced for several days, even after the average IMF has returned to zero after 1 day. This AE background decays away over several days as the Dst index recovers, indicating that there is some contamination of the currents observed at the AE stations by the continuing enhanced equatorial ring current. For data averaged over all seasons, the critical frequencies are depressed at Slough by 1.3 MHz, which is close to the lower decile of the overall distribution of δf o Fl values. Taking 30-day periods around summer and winter solstice, the largest depression is 1.6 and 1.2 MHz, respectively. This seasonal dependence is confirmed by a similar study for a Southern Hemisphere station, Argentine Island, giving peak depressions of 1.8 MHz and 0.5 MHz for summer and winter. For the subset of turnings where δB z > 11.5 nT and P ≤ 5 nPa, the response of the geomagnetic indices is similar but smaller, while the change in δf o F2 has all but disappeared. This confirms that the energy deposited at high latitudes, which leads to the geomagnetic and ionospheric disturbances following a southward turning of the IMF, increases with the energy density (dynamic pressure) of the solar wind flow. The magnitude of all responses are shown to depend on δB z . At Slough, the peak depression always occurs when Slough rotates into the noon sector. The largest ionospheric response is for southward turnings seen between 15–21 UT.