Showing papers in "Annales Geophysicae in 2017"

Exceptionally strong summer-like zonal wind reversal in the upper mesosphere during winter 2015/16

Abstract: . The 2015/16 Northern Hemisphere winter season was marked by peculiarities in the circulation pattern in the high-latitude mesopause region. Wind measurements from the Andenes (69° N, 13° E) meteor radar show westward winds below 84 km and eastward winds above. This wind pattern in the zonal wind was observable between the end of December 2015 and the end of January 2016, i.e., conditions that are typical for the summer were found during winter. Additional meteor radar measurements at midlatitude stations did not show such a zonal wind reversal but indicate, together with the polar latitude stations, a reversal of the horizontal temperature gradient. This is confirmed by global satellite measurements. Therefore, it is plausible that the polar latitude summer-like zonal wind reversal in December–January is in accordance with the reversed horizontal temperature gradient assuming a thermal wind balance between mid- and polar latitudes. The reversed horizontal temperature gradient itself is induced by stationary planetary waves at lower and midlatitudes in the mesosphere, leading to a weakening of the residual circulation above the European sector.

Modelling geomagnetically induced currents in midlatitude Central Europe using a thin-sheet approach

Abstract: Geomagnetically induced currents (GICs) in power systems, which can lead to transformer damage over the short and the long term, are a result of space weather events and geomagnetic variations. For a long time, only high-latitude areas were considered to be at risk from these currents, but recent studies show that considerable GICs also appear in midlatitude and equatorial countries. In this paper, we present initial results from a GIC model using a thin-sheet approach with detailed surface and subsurface conductivity models to compute the induced geoelectric field. The results are compared to measurements of direct currents in a transformer neutral and show very good agreement for short-period variations such as geomagnetic storms. Long-period signals such as quiet-day diurnal variations are not represented accurately, and we examine the cause of this misfit. The modelling of GICs from regionally varying geoelectric fields is discussed and shown to be an important factor contributing to overall model accuracy. We demonstrate that the Austrian power grid is susceptible to large GICs in the range of tens of amperes, particularly from strong geomagnetic variations in the east–west direction.

A comparison of 11-year mesospheric and lower thermospheric winds determined by meteor and MF radar at 69 ° N

Abstract: . The Andenes Meteor Radar (MR) and the Saura Medium Frequency (MF) Radar are located in northern Norway (69° N, 16° E) and operate continuously to provide wind measurements of the mesosphere and lower thermosphere (MLT) region. We compare the two systems to find potential biases between the radars and combine the data from both systems to enhance altitudinal coverage between 60 and 110 km. The systems have altitudinal overlap between 78 and 100 km at which we compare winds and tides on the basis of hourly winds with 2 km altitude bins. Our results indicate reasonable agreement for the zonal and meridional wind components between 78 and 92 km. An exception to this is the altitude range below 84 km during the summer, at which the correlation decreases. We also compare semidiurnal and diurnal tides according to their amplitudes and phases with good agreement below 90 km for the diurnal and below 96 km for the semidiurnal tides. Based on these findings we have taken the MR data as a reference. By comparing the MF and MR winds within the overlapping region, we have empirically estimated correction factors to be applied to the MF winds. Existing gaps in that data set will be filled with weighted MF data. This weighting is done due to underestimated wind values of the MF compared to the MR, and the resulting correction factors fit to a polynomial function of second degree within the overlapping area. We are therefore able to construct a consistent and homogenous wind from approximately 60 to 110 km.

Climatology and modeling of ionospheric scintillations and irregularity zonal drifts at the equatorial anomaly crest region

Abstract: . In this study the climatology of ionospheric scintillations and the zonal drift velocities of scintillation-producing irregularities are depicted for a station located under the southern crest of the equatorial ionization anomaly. Then, the α − μ ionospheric fading model is used for the first- and second-order statistical characterization of amplitude scintillations. In the statistical analyzes, data are used from single-frequency GPS receivers acquired during ∼ 17 years (September 1997–November 2014) at Cachoeira Paulista (22.4° S; 45.0° W), Brazil. The results reveal that the nocturnal occurrence of scintillations follows the seasonal distribution of plasma bubble irregularities observed in the longitudinal sector of eastern South America. In addition to the solar cycle dependence, the results suggest that the occurrence climatology of scintillations is also modulated by the secular variation in the dip latitude of Cachoeira Paulista, since the maximum occurrence of scintillations during the peak of solar cycle 24 was ∼ 20 % lower than that observed during the maximum of solar cycle 23. The dynamics of the irregularities throughout a solar cycle, as investigated from the estimates of the mean zonal drift velocities, presented a good correlation with the EUV and F10.7 cm solar fluxes. Meanwhile, the seasonal behavior showed that the magnitude of the zonal drift velocities is larger during the December solstice months than during the equinoxes. In terms of modeling, the results for the α − μ distribution fit quite well with the experimental data and with the temporal characteristics of fading events independently of the solar activity level.

The response of the H geocorona between 3 and 8 R e to geomagnetic disturbances studied using TWINS stereo Lyman- α data

Abstract: . Circumterrestrial Lyman-α column brightness observations from 3–8 Earth radii (Re) have been used to study temporal density variations in the exospheric neutral hydrogen as response to geomagnetic disturbances of different strength, i.e., Dst peak values between −26 and −147 nT. The data used were measured by the two Lyman-α detectors (LAD1/2) onboard both TWINS satellites between the solar minimum of 2008 and near the solar maximum of 2013. The solar Lyman-α flux at 121.6 nm is resonantly scattered near line center by exospheric H atoms and measured by the TWINS LADs. Along a line of sight (LOS), the scattered LOS-column intensity is proportional to the LOS H column density, assuming optically thin conditions above 3 Re. In the case of the eight analyzed geomagnetic storms we found a significant increase in the exospheric Lyman-α flux between 9 and 23 % (equal to the same increase in H column density ΔnH) compared to the undisturbed case short before the storm event. Even weak geomagnetic storms (e.g., Dst peak values ≥ −41 nT) under solar minimum conditions show increases up to 23 % of the exospheric H densities. The strong H density increase in the observed outer exosphere is also a sign of an enhanced H escape flux during storms. For the majority of the storms we found an average time shift of about 11 h between the time when the first significant dynamic solar wind pressure peak (pSW) hits the Earth and the time when the exospheric Lyman-α flux variation reaches its maximum. The results show that the (relative) exospheric density reaction of ΔnH have a tendency to decrease with increasing peak values of Dst index or the Kp index daily sum. Nevertheless, a simple linear correlation between ΔnH and these two geomagnetic indices does not seem to exist. In contrast, when recovering from the peak back to the undisturbed case, the Kp index daily sum and the ΔnH essentially show the same temporal recovery.

THEMIS satellite observations of hot flow anomalies at Earth's bow shock

Abstract: . Hot flow anomalies (HFAs) at Earth's bow shock were identified in Time History of Events and Macroscale Interactions During Substorms (THEMIS) satellite data from 2007 to 2009. The events were classified as young or mature and also as regular or spontaneous hot flow anomalies (SHFAs). The dataset has 17 young SHFAs, 49 mature SHFAs, 15 young HFAs, and 55 mature HFAs. They span a wide range of magnetic local times (MLTs) from approximately 7 to 16.5 MLT. The largest ratio of solar wind to HFA core density occurred near dusk and at larger distances from the bow shock. In this study, HFAs and SHFAs were observed up to 6.3 RE and 6.1 RE (Earth radii), respectively, upstream from the model bow shock. HFA–SHFA occurrence decreases with distance upstream from the bow shock. HFAs of the highest event core ion temperatures were not seen at the flanks. The ratio of HFA ion temperature increase to HFA electron temperature increase is highest around 12 MLT and slightly duskward. For SHFAs, (Tihfa∕Tisw)/(Tehfa∕Tesw) generally increased with distance from the bow shock. Both mature and young HFAs are more prevalent when there is an approximately radial interplanetary magnetic field. HFAs occur most preferentially for solar wind speeds from 550 to 600 km s−1. The correlation coefficient between the HFA increase in thermal energy density from solar wind values and the decrease in kinetic energy density from solar wind values is 0.62. SHFAs and HFAs do not show major differences in this study.

Atmospheric loss from the dayside open polar region and its dependence on geomagnetic activity: implications for atmospheric escape on evolutionary timescales

Abstract: . We have investigated the total O+ escape rate from the dayside open polar region and its dependence on geomagnetic activity, specifically Kp. Two different escape routes of magnetospheric plasma into the solar wind, the plasma mantle, and the high-latitude dayside magnetosheath have been investigated separately. The flux of O+ in the plasma mantle is sufficiently fast to subsequently escape further down the magnetotail passing the neutral point, and it is nearly 3 times larger than that in the dayside magnetosheath. The contribution from the plasma mantle route is estimated as ∼ 3. 9 × 1024exp(0. 45 Kp) [s−1] with a 1 to 2 order of magnitude range for a given geomagnetic activity condition. The extrapolation of this result, including escape via the dayside magnetosheath, indicates an average O+ escape of 3 × 1026 s−1 for the most extreme geomagnetic storms. Assuming that the range is mainly caused by the solar EUV level, which was also larger in the past, the average O+ escape could have reached 1027–28 s−1 a few billion years ago. Integration over time suggests a total oxygen escape from ancient times until the present roughly equal to the atmospheric oxygen content today.

Longitudinal variation of equatorial electrojet and the occurrence of its counter electrojet

Abstract: . We examined the longitudinal variability of the equatorial electrojet (EEJ) and the occurrence of its counter electrojet (CEJ) using the available records of the horizontal component H of the geomagnetic field simultaneously recorded in the year 2009 (mean annual sunspot number Rz = 3.1) along the magnetic equator in the South American, African, and Philippine sectors. Our results indicate that the EEJ undergoes variability from one longitudinal representative station to another, with the strongest EEJ of about 192.5 nT at the South American axis at Huancayo and a minimum peak of 40.7 nT at Ilorin in western Africa. Obtained longitudinal inequality in the EEJ was explicable in terms of the effects of local winds, dynamics of migratory tides, propagating diurnal tide, and meridional winds. The African stations of Ilorin and Addis Ababa registered the greatest % of CEJ occurrence. Huancayo in South America, with the strongest electrojet strength, was found to have the least occurrence of the CEJ. It is suggested that activities that support strong EEJ inhibits the occurrence of the CEJ. Percentage of occurrence of the CEJ varied with seasons across the longitudes. The order of seasonal variation of morning occurrence does not tally with the evening occurrence order at any station. A semiannual equinoctial maximum in percentage of morning occurrence of the CEJ was obtained at Huancayo and Addis Ababa. Only Addis Ababa recorded equal equinoctial maxima in percentage of evening occurrence of the CEJ. The seasonal distribution of the occurrences of the CEJ at different time regimes implies a seasonal variability of causative mechanisms responsible for the occurrence of the CEJ.

Two-stream instabilities from the lower-hybrid frequency to the electron cyclotron frequency: application to the front of quasi-perpendicular shocks

Abstract: . Quasi-perpendicular supercritical shocks are characterized by the presence of a magnetic foot due to the accumulation of a fraction of the incoming ions that is reflected by the shock front. There, three different plasma populations coexist (incoming ion core, reflected ion beam, electrons) and can excite various two-stream instabilities (TSIs) owing to their relative drifts. These instabilities represent local sources of turbulence with a wide frequency range extending from the lower hybrid to the electron cyclotron. Their linear features are analyzed by means of both a dispersion study and numerical PIC simulations. Three main types of TSI and correspondingly excited waves are identified: i. Oblique whistlers due to the (so-called fast ) relative drift between reflected ions/electrons; the waves propagate toward upstream away from the shock front at a strongly oblique angle (θ ∼ 50°) to the ambient magnetic field Bo, have frequencies a few times the lower hybrid, and have wavelengths a fraction of the ion inertia length c∕ωpi. ii. Quasi-perpendicular whistlers due to the (so-called slow ) relative drift between incoming ions/electrons; the waves propagate toward the shock ramp at an angle θ a few degrees off 90°, have frequencies around the lower hybrid, and have wavelengths several times the electron inertia length c∕ωpe. iii. Extended Bernstein waves which also propagate in the quasi-perpendicular domain, yet are due to the (so-called fast ) relative drift between reflected ions/electrons; the instability is an extension of the electron cyclotron drift instability (normally strictly perpendicular and electrostatic) and produces waves with a magnetic component which have frequencies close to the electron cyclotron as well as wavelengths close to the electron gyroradius and which propagate toward upstream. Present results are compared with previous works in order to stress some features not previously analyzed and to define a more synthetic view of these TSIs.

Near real-time estimation of ionosphere vertical total electron content from GNSS satellites using B-splines in a Kalman filter

Abstract: . Although the number of terrestrial global navigation satellite system (GNSS) receivers supported by the International GNSS Service (IGS) is rapidly growing, the worldwide rather inhomogeneously distributed observation sites do not allow the generation of high-resolution global ionosphere products. Conversely, with the regionally enormous increase in highly precise GNSS data, the demands on (near) real-time ionosphere products, necessary in many applications such as navigation, are growing very fast. Consequently, many analysis centers accepted the responsibility of generating such products. In this regard, the primary objective of our work is to develop a near real-time processing framework for the estimation of the vertical total electron content (VTEC) of the ionosphere using proper models that are capable of a global representation adapted to the real data distribution. The global VTEC representation developed in this work is based on a series expansion in terms of compactly supported B-spline functions, which allow for an appropriate handling of the heterogeneous data distribution, including data gaps. The corresponding series coefficients and additional parameters such as differential code biases of the GNSS satellites and receivers constitute the set of unknown parameters. The Kalman filter (KF), as a popular recursive estimator, allows processing of the data immediately after acquisition and paves the way of sequential (near) real-time estimation of the unknown parameters. To exploit the advantages of the chosen data representation and the estimation procedure, the B-spline model is incorporated into the KF under the consideration of necessary constraints. Based on a preprocessing strategy, the developed approach utilizes hourly batches of GPS and GLONASS observations provided by the IGS data centers with a latency of 1 h in its current realization. Two methods for validation of the results are performed, namely the self consistency analysis and a comparison with Jason-2 altimetry data. The highly promising validation results allow the conclusion that under the investigated conditions our derived near real-time product is of the same accuracy level as the so-called final post-processed products provided by the IGS with a latency of several days or even weeks.

Tail reconnection in the global magnetospheric context: Vlasiator first results

Abstract: . The key dynamics of the magnetotail have been researched for decades and have been associated with either three-dimensional (3-D) plasma instabilities and/or magnetic reconnection. We apply a global hybrid-Vlasov code, Vlasiator, to simulate reconnection self-consistently in the ion kinetic scales in the noon–midnight meridional plane, including both dayside and nightside reconnection regions within the same simulation box. Our simulation represents a numerical experiment, which turns off the 3-D instabilities but models ion-scale reconnection physically accurately in 2-D. We demonstrate that many known tail dynamics are present in the simulation without a full description of 3-D instabilities or without the detailed description of the electrons. While multiple reconnection sites can coexist in the plasma sheet, one reconnection point can start a global reconfiguration process, in which magnetic field lines become detached and a plasmoid is released. As the simulation run features temporally steady solar wind input, this global reconfiguration is not associated with sudden changes in the solar wind. Further, we show that lobe density variations originating from dayside reconnection may play an important role in stabilising tail reconnection.

Climatology of thermospheric neutral winds over Oukaïmeden Observatory in Morocco

Abstract: . In order to explore coupling between the thermosphere and ionosphere and to address the lack of data relating to thermospheric neutral winds and temperatures over the African sector, a new system of instruments was installed at the Oukaimeden Observatory located in the high Atlas Mountains, 75 km south of Marrakesh, Morocco (31.206° N, 7.866° W, 22.84° N magnetic). In this work we present the first multi-year results of the climatology of meridional and zonal winds obtained during the period from January 2014 to February 2016, including observations from 648 nights. The measurements are obtained using an imaging Fabry–Perot interferometer, which measures the 630.0 nm emissions caused by dissociative recombination of O2+. The basic climatology of the winds is as expected, showing zonal winds that are strongly eastward in the early evening just after sunset with a speed of 50 to 100 m s−1 decreasing in magnitude, and reversing directions in the local summer months, towards sunrise. The meridional winds are slightly poleward in the early evening during the local winter, before reversing directions around 21:00 LT. In the local summer months, the meridional winds are equatorward for the entire night, reaching a maximum equatorward speed of 75 m s−1. We compare the observed climatologies of neutral winds to that provided by the recently updated Horizontal Wind Model (HWM14) in order to validate that model's predictions of the thermospheric wind patterns over the eastern portion of Africa. The model captures much of the features in the observational climatologies. The most notable exception is for the zonal winds during local summer, when the maximum eastward wind in the observations occurs approximately 4 h later than seen in the model results.

New parameterized model for GPS water vapor tomography

Abstract: . Water vapor is the basic parameter used to describe atmospheric conditions. It is rarely contained in the atmosphere during the water cycle, but it is the most active element in rapid space–time changes. Measuring and monitoring the distribution and quantity of water vapor is a necessary task. GPS tomography is a powerful means of providing high spatiotemporal resolution of water vapor density. In this paper, a spatial structure model of a humidity field is constructed using voxel nodes, and new parameterizations for acquiring data about water vapor in the troposphere via GPS are proposed based on inverse distance weighted (IDW) interpolation. Unlike the density of water vapor that is constant within a voxel, the density at a certain point is determined by IDW interpolation. This algorithm avoids the use of horizontal constraints to smooth voxels that are not crossed by satellite rays. A prime number decomposition (PND) access order scheme is introduced to minimize correlation between slant wet delay (SWD) observations. Four experimental schemes for GPS tomography are carried out in dry weather from 2 to 8 August 2015 and rainy days from 9 to 15 August 2015. Using 14 days of data from the Hong Kong Satellite Positioning Reference Station Network (SatRef), the results indicate that water vapor density derived from 4-node methods is more robust than that derived from that of 8 nodes or 12 nodes, or that derived from constant refractivity schemes and the new method has better performance under stable weather conditions than unstable weather (e.g., rainy days). The results also indicate that an excessive number of interpolations in each layer reduce accuracy. However, the accuracy of the tomography results is gradually reduced with increases in altitude below 7000 m. Moreover, in the case of altitudes between 7000 m and the upper boundary layer, the accuracy can be improved by a boundary constraint.

Statistical study of auroral omega bands

Abstract: . The presence of very few statistical studies on auroral omega bands motivated us to test-use a semi-automatic method for identifying large-scale undulations of the diffuse aurora boundary and to investigate their occurrence. Five identical all-sky cameras with overlapping fields of view provided data for 438 auroral omega-like structures over Fennoscandian Lapland from 1996 to 2007. The results from this set of omega band events agree remarkably well with previous observations of omega band occurrence in magnetic local time (MLT), lifetime, location between the region 1 and 2 field-aligned currents, as well as current density estimates. The average peak emission height of omega forms corresponds to the estimated precipitation energies of a few keV, which experienced no significant change during the events. Analysis of both local and global magnetic indices demonstrates that omega bands are observed during substorm expansion and recovery phases that are more intense than average substorm expansion and recovery phases in the same region. The omega occurrence with respect to the substorm expansion and recovery phases is in a very good agreement with an earlier observed distribution of fast earthward flows in the plasma sheet during expansion and recovery phases. These findings support the theory that omegas are produced by fast earthward flows and auroral streamers, despite the rarity of good conjugate observations.

MONITOR Ionospheric Network: two case studies on scintillation and electron content variability

Abstract: . The ESA MONITOR network is composed of high-frequency-sampling global navigation satellite systems (GNSS) receivers deployed mainly at low and high latitudes to study ionosphere variability and jointly with global GNSS data and ionospheric processing software in support of the GNSS and its satellite-based augmentation systems (SBAS) like the European EGNOS. In a recent phase of the project, the network was merged with the CNES/ASECNA network and new receivers were added to complement the latter in the western African sector. This paper summarizes MONITOR, presenting two case studies on scintillations (using almost 2 years of data measurements). The first case occurred during the major St. Patrick's Day geomagnetic storm in 2015. The second case study was performed in the last phase of the project, which was supported by ESA EGNOS Project Office, when we paid special attention to extreme events that might degrade the system performance of the European EGNOS.

Diurnal variations in integrated water vapor derived from a GPS ground network in the Volga–Ural region of Russia

Abstract: . In this article, we present estimates of diurnal and semidiurnal harmonics of variations in integrated water vapor content (IWV) according to data from 16 GPS stations in the Volga–Ural region of Russia during 2013–2015. Amplitudes of diurnal harmonics are maximal in summer and reach values from 0.37 to 1.01 mm. Time at the maximum of diurnal harmonic is typically in the period from 14:00 to 17:00. Semidiurnal harmonics have the largest amplitudes in spring and autumn, but they do not exceed 0.19 mm. A comparison of the diurnal cycle from GPS data and ERA-Interim reanalysis has revealed significant differences in the phase. It is established that, as a result of evaporation from the underlying surface and convective lifting of moist air, the summer diurnal variations in IWV and surface density of water vapor are in antiphase. The diurnal cycle of IWV is determined by surface air temperature to be 88 % in summer and less than at 35 % in other seasons. It is noted that maximal amplitudes of diurnal harmonic of IWV are observed at stations located on the windward side of mountains.

An improved troposphere tomographic approach considering the signals coming from the side face of the tomographic area

Abstract: . The spatio-temporal distribution of atmospheric water vapour information plays a crucial role in the establishment of modern numerical weather forecast models and description of the different weather variations. A troposphere tomographic method has been proposed considering the signal rays penetrating from the side of the area of interest to solve the problem of the low utilisation rate of global navigation satellite system (GNSS) observations. Given the method above needs the establishment of a unit scale factor model using the radiosonde data at only one location in the research area, an improved approach is proposed by considering the reasonability of modelling data and the diversity of the modelling parameters for building a more accurate unit scale factor model. The new established model is established using grid point data derived from the European Centre for Medium-Range Weather Forecasts (ECMWF) and evenly distributed in the tomographic area, which can enhance the number of calculated initial water vapour density values with high accuracy. We validated the improved method with respect to the previous methods, as well as the result from a radiosonde using data from 12 stations from the Hong Kong Satellite Positioning Reference Station Network. The obtained result shows that the number of initial values estimated by the new model is increased by 6.83 %, while the internal and external accuracies are 0.08 and 0.24 g m−3, respectively. Integrated water vapour (IWV) and water vapour density profile comparisons show that the improved method is superior to previous studies in terms of RMS, MAE, and bias, which suggests higher accuracy and reliability.

Large-scale gravity wave perturbations in the mesopause region above Northern Hemisphere midlatitudes during autumnal equinox: a joint study by the USU Na lidar and Whole Atmosphere Community Climate Model

Abstract: . To investigate gravity wave (GW) perturbations in the midlatitude mesopause region during boreal equinox, 433 h of continuous Na lidar full diurnal cycle temperature measurements in September between 2011 and 2015 are utilized to derive the monthly profiles of GW-induced temperature variance, T′2, and the potential energy density (PED). Operating at Utah State University (42° N, 112° W), these lidar measurements reveal severe GW dissipation near 90 km, where both parameters drop to their minima (∼ 20 K2 and ∼ 50 m2 s−2, respectively). The study also shows that GWs with periods of 3–5 h dominate the midlatitude mesopause region during the summer–winter transition. To derive the precise temperature perturbations a new tide removal algorithm suitable for all ground-based observations is developed to de-trend the lidar temperature measurements and to isolate GW-induced perturbations. It removes the tidal perturbations completely and provides the most accurate GW perturbations for the ground-based observations. This algorithm is validated by comparing the true GW perturbations in the latest mesoscale-resolving Whole Atmosphere Community Climate Model (WACCM) with those derived from the WACCM local outputs by applying this newly developed tidal removal algorithm.

Seismo-ionospheric anomalies in ionospheric TEC and plasma density before the 17 July 2006 M 7.7 south of Java earthquake

Abstract: . In this paper, we report significant evidence for preseismic ionospheric anomalies in total electron content (TEC) of the global ionosphere map (GIM) and plasma density appearing on day 2 before the 17 July 2006 M7.7 south of Java earthquake. After distinguishing other anomalies related to the geomagnetic activities, we found a temporal precursor around the epicenter on day 2 before the earthquake (15 July 2006), which agrees well with the spatial variations in latitude–longitude–time (LLT) maps. Meanwhile, the sequences of latitude–time–TEC (LTT) plots reveal that the TECs on epicenter side anomalously decrease and lead to an anomalous asymmetric structure with respect to the magnetic equator in the daytime from day 2 before the earthquake. This anomalous asymmetric structure disappears after the earthquake. To further confirm these anomalies, we studied the plasma data from DEMETER satellite in the earthquake preparation zone (2046.4 km in radius) during the period from day 45 before to day 10 after the earthquake, and also found that the densities of both electron and total ion in the daytime significantly increase on day 2 before the earthquake. Very interestingly, O+ density increases significantly and H+ density decreases, while He+ remains relatively stable. These results indicate that there exists a distinct preseismic signal (preseismic ionospheric anomaly) over the epicenter.

Scintillation measurements at Bahir Dar during the high solar activity phase of solar cycle 24

Abstract: . Small-scale ionospheric disturbances may cause severe radio scintillations of signals transmitted from global navigation satellite systems (GNSSs). Consequently, small-scale plasma irregularities may heavily degrade the performance of current GNSSs such as GPS, GLONASS or Galileo. This paper presents analysis results obtained primarily from two high-rate GNSS receiver stations designed and operated by the German Aerospace Center (DLR) in cooperation with Bahir Dar University (BDU) at 11.6° N, 37.4° E. Both receivers collect raw data sampled at up to 50 Hz, from which characteristic scintillation parameters such as the S4 index are deduced. This paper gives a first overview of the measurement set-up and the observed scintillation events over Bahir Dar in 2015. Both stations are located close to one another and aligned in an east–west, direction which allows us to estimate the zonal drift velocity and spatial dimension of equatorial ionospheric plasma irregularities. Therefore, the lag times of moving electron density irregularities and scintillation patterns are derived by applying cross-correlation analysis to high-rate measurements of the slant total electron content (sTEC) along radio links between a GPS satellite and both receivers and to the associated signal power, respectively. Finally, the drift velocity is derived from the estimated lag time, taking into account the geometric constellation of both receiving antennas and the observed GPS satellites.

Measurements of the total ozone column using a Brewer spectrophotometer and TOMS and OMI satellite instruments over the Southern Space Observatory in Brazil

Abstract: This paper presents 23 years (1992–2014) of quasi-continuous measurements of the total ozone column (TOC) over the Southern Space Observatory (SSO) in Sao Martinho da Serra, Brazil (29.26 • S, 53.48 • and 488 m altitude). The TOC was measured by a Brewer spectrometer, and the results are also compared to daily and monthly observations from the TOMS (Total Ozone Mapping Spectrometer) and OMI (Ozone Monitoring Instrument) satellite instruments. Analyses of the main interannual modes of variability computed using the wavelet transform method were performed. A favorable agreement between the Brewer spec-trophotometer and satellite datasets was found. The seasonal TOC variation is dominated by an annual cycle, with a minimum of approximately 260 DU in April and a maximum of approximately 295 DU in September. The wavelet analysis applied in the SSO TOC anomaly time series revealed that the Quasi-Biennial Oscillation (QBO) modulation was the main mode of interannual variability. The comparison between the SSO TOC anomaly time series with the QBO index revealed that the two are in opposite phases.

Similarity and differences in morphology and mechanisms of the fo F2 and TEC disturbances during the geomagnetic storms on 26–30 September 2011

Abstract: . This study presents an analysis of the ground-based observations and model simulations of ionospheric electron density disturbances at three longitudinal sectors (eastern European, Siberian and American) during geomagnetic storms that occurred on 26–30 September 2011. We use the Global Self-consistent Model of the Thermosphere, Ionosphere and Protonosphere (GSM TIP) to reveal the main mechanisms influencing the storm-time behavior of the total electron content (TEC) and the ionospheric F2 peak critical frequency (foF2) during different phases of geomagnetic storms. During the storm's main phase the long-lasting positive disturbances in TEC and foF2 at sunlit mid-latitudes are mainly explained by the storm-time equatorward neutral wind. The effects of eastward electric field can only explain the positive ionospheric storm in the first few hours of the initial storm phase. During the main phase the ionosphere was more changeable than the plasmasphere. The positive disturbances in the electron content at the plasmaspheric heights (800–20 000 km) at high latitudes can appear simultaneously with the negative disturbances in TEC and foF2. The daytime positive disturbances in foF2 and TEC occurred at middle and low latitudes and at the Equator due to n(O) ∕ n(N2) enhancement during later stage of the main phase and during the recovery phase of the geomagnetic storm. The plasma tube diffusional depletion and negative disturbances in electron and neutral temperature were the main formation mechanisms of the simultaneous formation of the positive disturbances in foF2 and negative disturbances in TEC at low latitudes during the storm's recovery phase.

New results on equatorial thermospheric winds and temperatures from Ethiopia, Africa

Abstract: . Measurements of equatorial thermospheric winds, temperatures, and 630 nm relative intensities were obtained using an imaging Fabry–Perot interferometer (FPI), which was recently deployed at Bahir Dar University in Ethiopia (11.6° N, 37.4° E, 3.7° N magnetic). The results obtained in this study cover 6 months (53 nights of useable data) between November 2015 and April 2016. The monthly-averaged values, which include local winter and equinox seasons, show the magnitude of the maximum monthly-averaged zonal wind is typically within the range of 70 to 90 ms−1 and is eastward between 19:00 and 21:00 LT. Compared to prior studies of the equatorial thermospheric wind for this local time period, the magnitude is considerably weaker as compared to the maximum zonal wind speed observed in the Peruvian sector but comparable to Brazilian FPI results. During the early evening, the meridional wind speeds are 30 to 50 ms−1 poleward during the winter months and 10 to 25 ms−1 equatorward in the equinox months. The direction of the poleward wind during the winter months is believed to be mainly caused by the existence of the interhemispheric wind flow from the summer to winter hemispheres. An equatorial wind surge is observed later in the evening and is shifted to later local times during the winter months and to earlier local times during the equinox months. Significant night-to-night variations are also observed in the maximum speed of both zonal and meridional winds. The temperature observations show the midnight temperature maximum (MTM) to be generally present between 00:30 and 02:00 LT. The amplitude of the MTM was ∼ 110 K in January 2016 with values smaller than this in the other months. The local time difference between the appearance of the MTM and a pre-midnight equatorial wind was generally 60 to 180 min. A meridional wind reversal was also observed after the appearance of the MTM (after 02:00 LT). Climatological models, HWM14 and MSIS-00, were compared to the observations and the HWM14 model generally predicted the zonal wind observations well with the exception of higher model values by 25 ms−1 in the winter months. The HWM14 model meridional wind showed generally good agreement with the observations. Finally, the MSIS-00 model overestimated the temperature by 50 to 75 K during the early evening hours of local winter months. Otherwise, the agreement was generally good, although, in line with prior studies, the model failed to reproduce the MTM peak for any of the 6 months compared with the FPI data.

Observations of equatorial ionization anomaly over Africa and Middle East during a year of deep minimum

Abstract: . In this work, we investigated the veracity of an ion continuity equation in controlling equatorial ionization anomaly (EIA) morphology using total electron content (TEC) of 22 GPS receivers and three ground-based magnetometers (Magnetic Data Acquisition System, MAGDAS) over Africa and the Middle East (Africa–Middle East) during the quietest periods. Apart from further confirmation of the roles of equatorial electrojet (EEJ) and integrated equatorial electrojet (IEEJ) in determining hemispheric extent of EIA crest over higher latitudes, we found some additional roles played by thermospheric meridional neutral wind. Interestingly, the simultaneous observations of EIA crests in both hemispheres of Africa–Middle East showed different morphology compared to that reported over Asia. We also observed interesting latitudinal twin EIA crests domiciled at the low latitudes of the Northern Hemisphere. Our results further showed that weak EEJ strength associated with counter electrojet (CEJ) during sunrise hours could also trigger twin EIA crests over higher latitudes.

Urbanization effect on trends in sunshine duration in China

Abstract: . There is an ongoing debate on whether the observed decadal variations in surface solar radiation, known as dimming and brightening periods, are a large-scale or solely local phenomenon. We investigated this issue using long-term sunshine duration records from China, which experienced a rapid increase in urbanization during the past decades. Over the period 1960–2013, 172 pairs of urban and nearby rural stations were analyzed. Urban and rural sunshine duration trends show similar spatial patterns during a dimming phase (1960–1989) and a subsequent period during which trends were leveling off (1990–2013). This indicates that rather than local effects, the trends in sunshine duration are on more of a national or regional scale in China. Nevertheless, in the dimming phase, the declining rate of sunshine duration in rural areas is around two-thirds of that in urban areas. The ratio of rural to urban dimming generally increases from a minimum of 0.39 to a maximum of 0.87 with increasing indices of urbanization calculated based on the year 2013. It reaches a maximum when the urbanization level exceeds 50 %, the urban population exceeds 20 million, or the population density becomes higher than 250 person km−2. After the transition into the leveling-off period, sunshine duration trends are no longer significantly affected by urbanization. Meanwhile, the number of laws and regulations related to air pollution and investment in pollution treatment have been increasing in China.

Studying the variability in the diurnal and seasonal variations in GPS total electron content over Nigeria

Abstract: The study of diurnal and seasonal variations in total electron content (TEC) over Nigeria has been prompted by the recent increase in the number of GPS continuously operating reference stations (CORSs) across Nigeria as well as the reduced costs of microcomputing. The GPS data engaged in this study were recorded in the year 2012 at nine stations in Nigeria located between geomagnetic latitudes – 4.33 and 0.72° N. The GPS data were used to derive GPS TEC, which was analysed for diurnal and seasonal variations. The results obtained were used to produce local GPS TEC maps and bar charts. The derived GPS TEC across all the stations demonstrates consistent minimum diurnal variations during the pre-sunrise hours 04:00 to 06:00 LT, increases with sharp gradient during the sunrise period (∼ 07:00 to 09:00 LT), attains postnoon maximum at about 14:00 LT, and then falls to a minimum just before sunset. Generally, daytime variations are found to be greater than nighttime variations, which range between 0 and 5 TECU. The seasonal variation depicts a semi-annual distribution with higher values (∼ 25–30 TECU) around equinoxes and lower values (∼ 20–25 TECU) around solstices. The December Solstice magnitude is slightly higher than the June Solstice magnitude at all stations, while March Equinox magnitude is also slightly higher than September Equinox magnitude at all stations. Thus, the seasonal variation shows an asymmetry in equinoxes and solstices, with the month of October displaying the highest values of GPS TEC across the latitudes.

A statistical study of the motion of pulsating aurora patches: using the THEMIS All-Sky Imager

Abstract: . Previous studies of the motion of patches that comprise patchy pulsating aurora (PPA) have been based on a limited number of events. In this study, we use a large database of PPA patches obtained from the THEMIS all-sky imager at Gillam (66.18° magnetic latitude, 332.78° magnetic longitude) between May 2006 and July 2013 to explore the velocity of the PPA patches. Our results show that PPA patches mainly drift eastward after midnight and westward before midnight. In addition, we found that patch velocities are in the expected range of convection given the magnetic latitude and that the velocities do not seem to depend on auroral electrojet (AE) index. The results suggest that the drifts of auroral patches could be a proxy for the ionospheric convection, and possibly provide a convenient and accurate method to remotely sense the magnetospheric convection.

Comparison between IRI-2012 and GPS-TEC observations over the western Black Sea

Abstract: . The ionosphere is a dynamic layer which generally changes according to radiation emitted by the sun, the movement of the earth around the sun, and sunspot activity. Variations can generally be categorized as regular or irregular variations. Both types of variation have a huge effect on radio wave propagation. In this study, we have focused on the seasonal variation effect, which is one of the regular forms of variation in terms of the ionosphere. We examined the seasonal variation over the ZONG station in Turkey for the year 2014. Our analysis results and IRI-2012 present different ideas about ionospheric activity. According to our analysed results, the standard deviation reached a maximum value in April 2014. However, the maximum standard deviation obtained from IRI-2012 was seen in February 2014. Furthermore, it is clear that IRI-2012 underestimated the VTEC values when compared to our results for all the months analysed. The main source of difference between the two models is the IRI-2012 topside ionospheric representation. IRI-2012 VTEC has been produced as a result of the integration of an electron density profile within altitudinal limits of 60–2000 km. In other words, the main problem with regard to the IRI-2012 VTEC representation is not being situated in the plasmaspheric part of the ionosphere. Therefore we propose that the plasmaspheric part should be taken into account to calculate the correct TEC values in mid-latitude regions, and we note that IRI-2012 does not supply precise TEC values for use in ionospheric studies.

Error estimate of Taylor's frozen-in flow hypothesis in the spectral domain

Abstract: . The quality of Taylor's frozen-in flow hypothesis can be measured by estimating the amount of the fluctuation energy mapped from the streamwise wavenumbers onto the Doppler-shifted frequencies in the spectral domain. For a random sweeping case with a Gaussian variation of the large-scale flow, the mapping quality is expressed by the error function which depends on the mean flow speed, the sweeping velocity, the frequency bin, and the frequency of interest. Both hydrodynamic and magnetohydrodynamic treatments are presented on the error estimate of Taylor's hypothesis with examples from the solar wind measurements.

Spatial and temporal variability in MLT turbulence inferred from in situ and ground-based observations during the WADIS-1 sounding rocket campaign

Abstract: . In summer 2013 the WADIS-1 sounding rocket campaign was conducted at the Andoya Space Center (ACS) in northern Norway (69° N, 16° E). Among other things, it addressed the question of the variability in mesosphere/lower thermosphere (MLT) turbulence, both in time and space. A unique feature of the WADIS project was multi-point turbulence sounding applying different measurement techniques including rocket-borne ionization gauges, VHF MAARSY radar, and VHF EISCAT radar near Tromso. This allowed for horizontal variability to be observed in the turbulence field in the MLT at scales from a few to 100 km. We found that the turbulence dissipation rate, e varied in space in a wavelike manner both horizontally and in the vertical direction. This wavelike modulation reveals the same vertical wavelengths as those seen in gravity waves. We also found that the vertical mean value of radar observations of e agrees reasonably with rocket-borne measurements. In this way defined 〈eradar〉 value reveals clear tidal modulation and results in variation by up to 2 orders of magnitude with periods of 24 h. The 〈eradar〉 value also shows 12 h and shorter (1 to a few hours) modulations resulting in one decade of variation in 〈eradar〉 magnitude. The 24 h modulation appeared to be in phase with tidal change of horizontal wind observed by SAURA-MF radar. Such wavelike and, in particular, tidal modulation of the turbulence dissipation field in the MLT region inferred from our analysis is a new finding of this work.