A numerical model approximating extreme energetic electron events involved in the physical and chemical processes of the middle atmosphere

Relativistic electrons (with energies >150 keV) which originate in the outer radiation belt and detected by the Russian ‘Meteor’ series of satellites have been correlated with the atmospheric total ozone data compiled by almost 90 stations located around the world within the latitude zone 40°–70°N. In more than 60% of the stations examined we have detected a clear decrease of the ozone 3–5 days after the electron flux excess. A numerical model has been applied to approximate this effect based on relativistic electron initiated nitric oxides creation in the upper mesosphere with subsequent atmospheric transport (both vertical and horizontal) towards the upper stratosphere. A first attempt of local and temporal prediction of ozone depletion because of energetic electrons impact in the middle atmosphere has been illustrated.     

New evidence of solar variation in temperature proxies from Northern Fennoscandia

A new summer temperature proxy was built for northern Fennoscandia in AD 1000–2004 using parameters of tree growth from a large region, extending from the Swedish Scandes to the Kola Peninsula. It was found that century-scale (55–140 year) cyclicity is present in this series during the entire time interval. This periodicity is highly significant and has a bi-modal structure, i.e. it consists of two oscillation modes, 55–100 year and 100–140 year variations. A comparison of the century-long variation in the northern Fennoscandian temperature proxy with the corresponding variations in Wolf numbers and concentration of cosmogenic ¹⁰Be in glacial ice shows that a probable cause of this periodicity is the modulation of regional climate by the secular solar cycle of Gleissberg. This is in line with the results obtained previously for a more limited part of the region (Finnish Lapland: 68–70° N, 20–30° E). Thus the reality of a link between long-term changes in solar activity and climate in Fennoscandia has been confirmed. Possible mechanisms of solar influence on the lower troposphere are discussed.     

The response of sporadic E-layer to the total solar eclipse of July 22, 2009 over the equatorial ionization anomaly region of the Indian zone

The digital ionosonde located in Bhopal (23.2°N, 77.2°E), India has been used to investigate the responses of the Es layer in the equatorial ionization anomaly (EIA) crest to the total solar eclipse (TSE) of July 22, 2009. Results show the presence of intense Es layer during and after the eclipse period. The gravity waves induced by the solar eclipse propagated upward in the Es layer and produced the periodic disturbance. The results of the wavelet analysis display the presence of dominant oscillation of about 24–32, 16–20 and 8 min. The appearance of intense sporadic-E concomitantly with the signatures of gravity wave suggests that the wind shear introduced by the solar eclipse induced gravity wave might be the plausible mechanism behind the intensification of Es-layer ionization.     

Behavior of the atomic oxygen 557.7 nm atmospheric emission in the solar cycle 23

We present the results of nightglow observation of the atomic oxygen 557.7 nm line emission in the solar cycle 23. We use the experimental data obtained at Geophysical observatory near Irkutsk (52°N, 103°E), Russia, for the 1997–2006 period. The 557.7 nm emission observations data are compared with atmospheric and solar parameters. We note a difference in correlation coefficients between the 557.7 nm emission intensity and the solar activity indices in different phases of the solar cycle. Airglow observation results are compared with the observational data obtained by other authors.     

Solar energetic particle events during the rise phases of solar cycles 23 and 24

We present a comparative study of the properties of coronal mass ejections (CMEs) and flares associated with the solar energetic particle (SEP) events in the rising phases of solar cycles (SC) 23 (1996–1998) (22 events) and 24 (2009–2011) (20 events), which are associated with type II radio bursts. Based on the SEP intensity, we divided the events into three categories, i.e. weak (intensity < 1 pfu), minor (1 pfu < intensity < 10 pfu) and major (intensity ? 10 pfu) events. We used the GOES data for the minor and major SEP events and SOHO/ERNE data for the weak SEP event. We examine the correlation of SEP intensity with flare size and CME properties. We find that most of the major SEP events are associated with halo or partial halo CMEs originating close to the sun center and western-hemisphere. The fraction of halo CMEs in SC 24 is larger than the SC 23. For the minor SEP events one event in SC23 and one event in SC24 have widths < 120° and all other events are associated with halo or partial halo CMEs as in the case of major SEP events. In case of weak SEP events, majority (more than 60%) of events are associated with CME width < 120°. For both the SC the average CMEs speeds are similar. For major SEP events, average CME speeds are higher in comparison to minor and weak events. The SEP event intensity and GOES X-ray flare size are poorly correlated. During the rise phase of solar cycle 23 and 24, we find north–south asymmetry in the SEP event source locations: in cycle 23 most sources are located in the south, whereas during cycle 24 most sources are located in the north. This result is consistent with the asymmetry found with sunspot area and intense flares.     

The asymmetrical features in electron density during extreme solar minimum

The variations of plasma density in topside ionosphere during 23rd/24th solar cycle minimum attract more attentions in recently years. In this analysis, we use the data of electron density (Ne) from DEMETER (Detection of Electromagnetic Emissions Transmitted from Earthquake Regions) satellite at the altitude of 660–710 km to investigate the solstitial and equinoctial asymmetry under geomagnetic coordinate system at LT (local time) 1030 and 2230 during 2005–2010, especially in solar minimum years of 2008–2009. The results reveal that ΔNe (December–June) is always positive over Southern Hemisphere and negative over northern part whatever at LT 1030 or 2230, only at 0–10°N the winter anomaly occurs with ΔNe (December–June) > 0, and its amplitude becomes smaller with the declining of solar flux from 2005 to 2009. The ΔNe between September and March is completely negative during 2005–2008, but in 2009, it turns to be positive at latitudes of 20°S–40°N at LT 1030 and 10°S–20°N at LT 2230. Furthermore, the solstitial and equinoctial asymmetry index (AI) are calculated and studied respectively, which all depends on local time, latitude and longitude. The notable differences occur at higher latitudes in solar minimum year of 2009 with those in 2005–2008. The equinoctial AI at LT 2230 is quite consistent with the variational trend of solar flux with the lowest absolute AI occurring in 2009, the extreme solar minimum, but the solstitial AI exhibits abnormal enhancement during 2008 and 2009 with bigger AI than those in 2005–2007. Compared with the neutral compositions at 500 km altitude, it illustrates that [O/N₂] and [O] play some roles in daytime and nighttime asymmetry of Ne at topside ionosphere.     

Equatorial Ionospheric Anomaly (EIA) and comparison with IRI model during descending phase of solar activity (2005–2009)

The ionospheric variability at equatorial and low latitude region is known to be extreme as compared to mid latitude region. In this study the ionospheric total electron content (TEC), is derived by analyzing dual frequency Global Positioning System (GPS) data recorded at two stations separated by 325 km near the Indian equatorial anomaly region, Varanasi (Geog latitude 25°, 16^/ N, longitude 82°, 59^/ E, Geomagnetic latitude 16°, 08^/ N) and Kanpur (Geog latitude 26°, 18^/ N, longitude 80°, 12^/ E, Geomagnetic latitude 17°, 18^/ N). Specifically, we studied monthly, seasonal and annual variations as well as solar and geomagnetic effects on the equatorial ionospheric anomaly (EIA) during the descending phase of solar activity from 2005 to 2009. It is found that the maximum TEC (EIA) near equatorial anomaly crest yield their maximum values during the equinox months and their minimum values during the summer. Using monthly averaged peak magnitude of TEC, a clear semi-annual variation is seen with two maxima occurring in both spring and autumn. Results also showed the presence of winter anomaly or seasonal anomaly in the EIA crest throughout the period 2005–2009 only except during the deep solar minimum year 2007–2008. The correlation analysis indicate that the variation of EIA crest is more affected by solar activity compared to geomagnetic activity with maximum dependence on the solar EUV flux, which is attributed to direct link of EUV flux on the formation of ionosphere and main agent of the ionization. The statistical mean occurrence of EIA crest in TEC during the year from 2005 to 2009 is found to around 12:54 LT hour and at 21.12° N geographic latitude. The crest of EIA shifts towards lower latitudes and the rate of shift of the crest latitude during this period is found to be 0.87° N/per year. The comparison between IRI models with observation during this period has been made and comparison is poor with increasing solar activity with maximum difference during the year 2005.     

Changes in total electron content (TEC) during the annular solar eclipse of 15 January 2010

The solar eclipse of 15 January 2010 was an annular eclipse of the Sun with a maximum magnitude of 0.96 at 1.62°N, 69.29°E. To study the effect of this solar eclipse on the ionosphere the GPS data recorded at three different Indian stations Varanasi (Geographic latitude 25°, 16′N, longitude 82°, 59′E), Hyderabad (Geographic latitude 17°, 20′N, longitude 78°, 30′E) and Bengaluru (Geographic latitude 12°, 58′N, longitude 77°, 33′E) have been used to retrieve ionospheric total electron content (TEC). The ionospheric response to this rare event has been studied in terms of GPS-derived TEC observed at all the three Indian stations. A significant reduction in TEC reflected by all PRNs at all the three stations has been observed. The magnitude of the reduction in VTEC compared to quiet mean VTEC depends on latitude as well as longitude. The amount of reduction observed from different satellites (PRN) is different and depends on the location of the satellite from the solar eclipse path.     

The solar eclipse and its associated ionospheric TEC variations over Indian stations on January 15, 2010

An annular solar eclipse occurred over the Indian subcontinent during the afternoon hours of January 15, 2010. This event was unique in the sense that solar activity was minimum and the eclipse period coincides with the peak ionization time at the Indian equatorial and low latitudes. The number of GPS receivers situated along the path of solar eclipse were used to investigate the response of total electron content (TEC) under the influence of this solar eclipse. These GPS receivers are part of the Indian Satellite Based Augmentation System (SBAS) named as ‘GAGAN’ (GPS Aided Geo Augmented Navigation) program. The eight GPS stations located over the wide range of longitudes allows us to differentiate between the various factors induced due to solar eclipse over the equatorial and low latitude ionosphere. The effect of the eclipse was detected in diurnal variations of TEC at all the stations along the eclipse path. The solar eclipse has altered the ionospheric behavior along its path by inducing atmospheric gravity waves, localized counter-electrojet and attenuation of solar radiation intensity. These three factors primarily control the production, loss and transport of plasma over the equatorial and low latitudes. The localized counter-electrojet had inhibited the equatorial ionization anomaly (EIA) in the longitude belt of 72°E–85°E. Thus, there was a negative deviation of the order of 20–40% at the equatorial anomaly stations lying in this ‘inhibited EIA region’. The negative deviation of only 10–20% is observed for the stations lying outside the ‘inhibited EIA region’. The pre-eclipse effect in the form of early morning enhancement of TEC associated with atmospheric gravity waves was also observed during this solar eclipse. More clear and distinctive spatial and temporal variations of TEC were detected along the individual satellite passes. It is also observed that TEC starts responding to the eclipse after 30 min from start of eclipse and the delay of the maximum TEC deviation from normal trend with respect to the maximum phase of the eclipse was close to one hour in the solar eclipse path.     

F2 layer response to geomagnetic disturbances in Eastern Asia under the low solar activity

In this paper, the peculiarities of ionospheric response to geomagnetic disturbances observed at the decay and minimum of solar activity (SA) in the period 2004–2007 are investigated with respect to different geomagnetic conditions. Data from ionospheric stations and results of total electron content (TEC) measurements made at the network of GPS ground-based receivers located within the latitude–longitude sector (20–70°N, 90–160°Е) are used in this study. Three groups of anomalous ionospheric response to geomagnetic disturbances have been observed during low solar activity. At daytime, the large-scale traveling ionospheric disturbances (LSTIDs) could generally be related to the main phase of magnetic storm. Quasi-two-days wavelike disturbances (WLDs) have been also observed in the main phase independent of the geomagnetic storm intensity. Sharp electron density oscillations of short duration (OSD) occurred in the response to the onset of both main and recovery phases of the magnetic storm in the daytime at middle latitudes. A numerical model for ionosphere–plasmasphere coupling was used to interpret the occurrence of LS TIDs. Results showed that the LSTIDs might be associated with the unexpected lifting of F2 layer to the region with the lower recombination rate by reinforced meridional winds that produces the increase of the electron density in the F2 layer maximum.     

Differential emission measure distributions in X-ray solar flares

X-ray spectrometer RESIK has observed spectra in the four wavelength bands from 3.3 Å to 6.1 Å. This spectral range contains many emission lines of H- and He-like ions for Si, S, Ar and K. These lines are formed in plasma of coronal temperatures (T > 3 MK). Analysis of their intensities allows studying differential emission measure distributions (DEM) in temperature range roughly between 3 MK and 30 MK. The aim of present study was to check whether any relationship exists between the character of DEM distribution, the event phase and the X-ray flare class. To do this we have calculated and analyzed the DEM distributions for a set of flares belonging to different GOES classes from the range B5.6–X1. The DEM distributions have been calculated using “Withbroe–Sylwester” multiplicative, maximum likelihood iterative algorithm. As the input data we have used absolute fluxes observed by RESIK in several spectral bands (lines + continuum). Respective emission functions have been calculated using the CHIANTI v. 5.2 atomic data package.     

Equatorial total electron content (TEC) at low and high solar activity

Total electron content (TEC) derived from ionosonde data recorded at the station of Korhogo (Lat = 9.33°N, Long = 5.43°W, Dip = 0.67°S) are compared to the International Reference Ionosphere (IRI) model predicted TEC for high (1999) and low (1994) solar activity conditions. The results show that the model represents the diurnal variation of the TEC as well as a solar activity and seasonal dependence. This variation is closer to that of the ionosonde-inferred TEC at high solar activity. However, at low solar activity the IRI overestimates the ionosonde-inferred TEC. The relative deviation ΔTEC is more prominent in the equinoctial seasons during nighttime hours where it is as high as 70%. At daytime hours, the relative deviation is estimated to 0–30%.     

Correlation between the cosmic noise absorption calculated from the SARINET data and the energetic particles measured by MEPED: Simultaneous observations over SAMA region

The cosmic noise absorption is presented in terms of two-dimensional images obtained from the imaging riometers operated at the Southern Space Observatory (geographic coordinate: 29.4° S, 53.1° W), in São Martinho da Serra, Brazil, Concepcion (geographic coordinate: 36.5° S, 73.0° W) and Punta Arenas (geographic coordinate: 53.0° S, 70.5° W) in Chile, which belong to the South American Riometer Network and are located at the central and periphery regions of the South American Magnetic Anomaly. Correlations are performed between the maximum cosmic noise absorption observed at these stations and the energetic electron flux in two energy channels (>30 and >300 keV) and the proton flux in three energy channels (80–240, 800–2500 and >6900 keV) as measured by the Medium Energy Proton and Electron Detector, during a moderate geomagnetic storm that occurred on September 3, 2008. The results show high correlations between the cosmic noise absorption detected at São Martinho da Serra and the flux of protons with energy between 80 and 240 keV, and the flux of electrons with energies higher than 300 keV, while an additional ionization at Concepcion was correlated with electrons of energies higher than 30 keV. The cosmic noise absorption detected at Punta Arenas was probably caused by the increase of the protons flux with energy between 80 and 240 keV.     

Occurrence rate of SAR arcs during the 23rd solar activity cycle

The occurrence rate of SAR arcs during 1997–2007 has been analyzed based on the photometric observations at the Yakutsk meridian (Maimaga station, corrected geomagnetic coordinates: 57°N, 200°E). SAR arcs appeared in 114 cases (∼500 h) during ∼370 nights of observations (∼3170 h). The occurrence frequency of SAR arcs increases to 27% during the growth phase of solar activity and has a clearly defined maximum at a decline of cycle 23. The SAR arc registration probability corresponds to the variations in geomagnetic activity in this solar cycle. The dates, intervals of UT, and geomagnetic latitudes of SAR arc observations at the Yakutsk meridian are presented.     

Ozone mini-hole observation over the Balkan Peninsula in March 2005

Areas with dimensions of 1000–3000 km in which the total ozone content (TOC) decreases fast are called ozone mini-holes. They are generated mainly dynamically in two ways, either by poor-ozone air mass transport from the tropics to higher latitudes by planetary wave activity or, they are connected with strong adiabatic uplifting of the tropopause height. An ozone mini-hole, generated by the second mechanism, was observed over the Balkan Peninsula on 19/21 March 2005. In the middle of March, the polar vortex was strongly disturbed by Rossby waves, reaching up to the lower stratosphere. Warming episodes over a geographical area, covering the Barents Sea and the Polar Sea north from Central Siberia, displaced a polar vortex fragment extremely southwards. However, the vorticity was weak and the stratospheric temperatures did not reach low values, providing conditions for ozone chemical destruction via heterogenic reactions. At the same time, a Rossby wave ridge was located below the European polar fragment. In the period from 13 to 19 March, the thermal tropopause over Sofia was uplifted almost by 3 km. Ozone distributions observed by the SCIAMACHY instrument on 18–21 March show a fast TOC decrease westwards from Ireland, which was moving eastwards during the next days, increasing the area in which the ozone content decreased. On 20/21 March low ozone content was observed above the Stara Zagora (42°N, 25°E) ground-based station by means of the GASCOD instrument, using DOAS technique. The TOMS Earth probe instrument detected 237 DU over Sofia. This is a record low March value from the beginning of the TOMS instrument measurements in 1978. In March/April the ozone distribution was characterized by its mean annual maximum of 360 DU at 42°N.     

SEP events and multi-spacecraft observations: Constraints on theory

Fifteen solar energetic particle (SEP) events have been analyzed using proton flux data recorded by the Helios 1, Helios 2, and IMP 8 spacecraft in the energy range ∼4–40 MeV during 1974–1982. For each of the events at least two of the spacecraft have their nominal magnetic footpoint within 20° in heliocentric longitude from each other. The SEP events are sub-grouped as a function of their heliocentric longitudinal separation and heliocentric radial distance from the SEP associated solar flare and several case studies are presented in this paper. Main results concerning their usage in estimating the SEP radial dependence are given. Moreover, we investigate the behavior of the third not connected spacecraft in order to study the dependence of the proton flux as a function of flare location. It is found that the contribution of the longitudinal gradient in determining variations in the SEP proton flux is particularly relevant for spacecraft having their magnetic connection footpoint separated from the flare between 30° and 50°.     

Flaring loop parameters estimated from solar decimeter type U-like and type J-like fine structures

This work presents the analysis of five fine structures in the solar radio emission, observed between June 2000 and October 2001 by the Brazilian Solar Spectroscope (BSS), in the decimeter frequency band of 950–2500 MHz. Based on their morphological characteristics identified in the dynamic spectra, the fine structures had been classified as type U-like or type J-like bursts. Such emissions are variants of the type III bursts. They support the hypothesis of generation by plasma emission mechanism, from interaction of electron beams accelerated during solar flares, propagating along closed magnetic structures, within the trapped plasma of the solar corona. The spectral and temporal characteristics of the five fine structures had been obtained from the dynamic spectra and the parameters of the agent and the emitting source have been determined, assuming both fundamental and harmonic emissions. The analysis revealed the flux density of the structures is less than 20–80 s.f.u. For assumption of harmonic emission, the interval of values for the source parameters estimated are: the loop size is (0.3–5.1) × 10¹⁰ cm; the electron beam velocity is in the range of 0.16–0.53 c; the temperature of coronal loop top is of the order of (0.25–1.55) × 10⁷ K; and the low limit for the magnetic field is of 7–26 G. These results are in agreement with previous determinations reported in the literature.     

Effective non-vertical and apparent cutoff rigidities for a cosmic ray latitude survey from Antarctica to Italy in minimum of solar activity

In this paper we will report the results of the computation of cutoff rigidities of vertical and non-vertical incident cosmic ray particles. Non-vertical effective cutoff rigidities have been computed by tracing particle trajectories through the “real” geomagnetic magnetic field comprising the International Geomagnetic Reference Field model (IGRF95, IAGA Division 5 Working Group 8, 1996: Sabaka, T.J., Langel, R.A., Baldwin, R.T., Conrad, J.A. The geomagnetic field, 1900–1995, including the large scale fields from magnetospheric sources and NASA candidate models for the 1995 IGRF revision. J. Geomag. Geoelect. 49, 157–206, 1997.) and the Tsyganenko [Tsyganenko, N.A. A magnetospheric magnetic field model with a warped tail current sheet. Planet. Space Sci. 37, 5–20, 1989.] magnetosphere model. The computation have been done for the backward route (from Antarctica to Italy) of the Italian Antarctic ship survey 1996–1997, for geographic points corresponding to the daily average coordinates of the ship; for zenith angles 15°, 30°, 45° and 60°, and azimuth angles from 0° to 360° in steps of 45°. By means of the obtained non-vertical cutoffs the apparent cutoff rigidities have been calculated. The information on integral multiplicities of secondary neutrons detected by the neutron monitor in dependence of the zenith angle of incoming primary cosmic ray particles have also been used. This information is based on the theoretical calculations of meson-nuclear cascades of primary protons with different rigidities arriving to the Earth’s atmosphere at the zenith angles of 0°, 15°, 30°, 45°, 60° and 75°. The difference between the computed apparent and vertical cutoff rigidities reaches ∼1 GV at rigidities >7–8 GV. At rigidities of 10–16 GV, the difference between the apparent and vertical cutoff rigidities is larger than that obtained earlier by Clem et al. [Clem, J.M., Bieber, J.W., Duldig, M., Evenson, P., Hall, D., Humble, J.E. Contribution of obliquely incident particles to neutron monitor counting rate. J. Geophys. Res. 102, 26919–26926, 1997.] and Dorman et al. [Dorman, L.I., Villoresi, G., Iucci, N., Parisi, M., Tyasto, M.I., Danilova, O.A., Ptitsyna, N.G. Cosmic ray survey to Antarctica and coupling functions for neutron component near solar minimum (1996–1997), 3. Geomagnetic effects and coupling functions. J. Geophys. Res. 105, 21047–21056, 2000.].     

Transient variations of vertical total electron content over some African stations from 2002 to 2012

This paper presents the vertical total electron content vTEC variations for three African stations, located at mid-low and equatorial latitudes, and operating since more than 10 years. The vTEC of the middle latitude GPS station in Alexandria, Egypt (31.2167°N; 29.9667°E, geographic) is compared to the vTEC of two others GPS stations: the first one in Rabat/Morocco (33.9981°N; 353.1457°E, geographic), and the second in Libreville/Gabon (0.3539°N; 9.6721°E, geographic). Our results discussed the diurnal, seasonal, and solar cycle dependences of vTEC at the local ionospheric conditions, during different phases of solar cycle in the light of the classification of Legrand and Simon. The vTEC over Alexandria exhibits the well-known equinoctial asymmetry which changes with the phases of the solar cycle; the spring vTEC is larger than that of autumn during the maximum, decreasing and minimum phases of solar cycle 23. During the increasing phase of solar cycle 24, it is the contrary. The diurnal variation of the vTEC presents multiple maxima during the equinox from 2005 to 2008 and during the summer solstice from 2006 to 2012. A nighttime vTEC enhancement and winter anomaly are also observed. During the deep solar minimum (2006–2009) the diurnal variation of the vTEC observed over Alexandria is similar to the diurnal variation observed during quiet magnetic period at equatorial latitudes. We observed also that the amplitude of vTEC at Libreville is larger than the amplitude of vTEC observed at Alexandria and Rabat, indeed Libreville is near the southern crest of the Equatorial Ionization anomaly. Finally, the correlation coefficient between vTEC and the sunspot number Rz is high and changes with solar cycle phases.     

The variation of equatorial spread-F occurrences observed by ionosondes at Thailand longitude sector

The equatorial spread-F (ESF) is a phenomenon of ionopheric irregularities which are mainly generated by the generalized Rayleigh–Taylor (R–T) instability mechanism in conjunction with the other physical mechanisms, originated at the bottom side of the F-layer in the equatorial region after sunset. It degrades the quality of signals that propagate through these irregularities, especially in the navigation satellite system, which requires the high integrity signals. In this work, we analyze the ESF statistics obtained from the FM/CW ionosonde stations over Thailand longitude sector. One is at Chumphon (10.72°N, 99.37°E, dip latitude 3.0°), located near the geomagnetic equator, and the other station is located at Chiangmai (18.76°N, 98.93°E, dip latitude 12.7°). Both stations are as part of the South-East Asia Low Latitude Ionospheric Network (SEALION) project. The ionograms are obtained at every 15 min from September 2004 to August 2005, which has the monthly mean of solar 10.7 cm flux (F10.7) from ∼80 to ∼110. In addition, we compare the diurnal patterns between the ESF occurrences and the variation of virtual height of the F-layer bottom side (h’F) of these two stations. The results show that the ESF occurrences at Chumphon stations are higher than Chiangmai station in all seasons. The high ESF occurrences of both stations mostly occur in equinoctial months corresponded with the rapid rising of the monthly mean h’F in the post-sunset. However, some inconsistent results are still observed, implying the role of other factors such as gravity waves and planetary waves to ESF occurrences.