Multi-year behaviour of the midnight OH1 temperature according to observations at Zvenigorod over 2000–2016

Using spectral measurements of the hydroxyl airglow at the Zvenigorod station (56° N, 37° E), Moscow region, over 2000–2016, we obtained the long-term set of data comprising 1822 midnight values of the OH¹ temperature in the mesopause region. These data revealed a 17-year series of its mean annual values, as well as amplitudes and phases of the first two harmonics of its annual variation. The obtained parameters were analyzed to determine statistically relevant characteristics of their long-term variations. As a result, we found that the long-term behaviour of the mean annual OH¹ temperature features a small negative linear trend (?0.07 ± 0.03 K/year) over the addressed period. Besides, its dependence on solar activity is shown to be 4.1 ± 0.5 K/100 SFU. Regarding the long-term behaviour of the mean annual OH¹ temperature, we revealed the existence of two oscillations with 3-year (the amplitude being 1.3 ± 0.2 K) and 4.1-year (the amplitude being 0.6 ± 0.2 K) periods. We obtained empirical relations describing year-to-year variations in the amplitudes and phases of the annual and semi-annual harmonics.     

Relationship of solar flare accelerated particles to solar energetic particles (SEPs) observed in the interplanetary medium

Observations of hard X-ray (HXR)/γ-ray continuum and γ-ray lines produced by energetic electrons and ions, respectively, colliding with the solar atmosphere, have shown that large solar flares can accelerate ions up to many GeV and electrons up to hundreds of MeV. Solar energetic particles (SEPs) are observed by spacecraft near 1 AU and by ground-based instrumentation to extend up to similar energies as in large SEP events, but it appears that a different acceleration process, one associated with fast coronal mass ejections is responsible. Much weaker SEP events are observed that are generally rich in electrons, ³He, and heavy elements. The energetic particles in these events appear to be similar to those accelerated in flares. The Ramaty high energy solar spectroscopic imager (RHESSI) mission provides high-resolution spectroscopy and imaging of flare HXRs and γ-rays. Such observations can provide information on the location, energy spectra, and composition of the flare accelerated energetic particles at the Sun. Here, preliminary comparisons of the RHESSI observations with observations of both energetic electron and ion near 1 AU are reviewed, and the implications for the particle acceleration and escape processes are discussed.     

Multi precursors analysis associated with the powerful Ecuador (MW = 7.8) earthquake of 16 April 2016 using Swarm satellites data in conjunction with other multi-platform satellite and ground data

After DEMETER satellite mission (2004–2010), the launch of the Swarm satellites (Alpha (A), Bravo (B) and Charlie (C)) has created a new opportunity in the study of earthquake ionospheric precursors. Nowadays, there is no doubt that multi precursors analysis is a necessary phase to better understand the LAIC (Lithosphere Atmosphere Ionosphere Coupling) mechanism before large earthquakes. In this study, using absolute scalar magnetometer, vector field magnetometer and electric field instrument on board Swarm satellites, GPS (Global Positioning System) measurements, MODIS-Aqua satellite and ECMWF (European Centre for Medium-Range Weather Forecasts) data, the variations of the electron density and temperature, magnetic field, TEC (Total Electron Content), LST (Land Surface Temperature), AOD (Aerosol Optical Depth) and SKT (SKin Temperature) have been surveyed to find the potential seismic anomalies around the strong Ecuador (M_w = 7.8) earthquake of 16 April 2016. The four solar and geomagnetic indices: F10.7, D_st, K_p and a_p were investigated to distinguish whether the preliminary detected anomalies might be associated with the solar-geomagnetic activities instead of the seismo-ionospheric anomalies. The Swarm satellites (A, B and C) data analysis indicate the anomalies in time series of electron density variations on 7, 11 and 12 days before the event; the unusual variations in time series of electron temperature on 8 days preceding the earthquake; the analysis of the magnetic field scalar and vectors data show the considerable anomalies 52, 48, 23, 16, 11, 9 and 7 days before the main shock. A striking anomaly is detected in TEC variations on 1 day before earthquake at 9:00 UTC. The analysis of MODIS-Aqua night-time images shows that LST increase unusually on 11 days prior to main shock. In addition, the AOD variations obtained from MODIS measurements reach the maximum value on 10 days before the earthquake. The SKT around epicentral region presents anomalous higher value about 40 days before the earthquake. It should be noted that the different lead times of the observed anomalies could be acknowledged based on a reasonable LAIC earthquake mechanism. Our results emphasize that the Swarm satellites measurements play an undeniable role in progress the studies of the ionospheric precursors.     

Satellite evidence of harmful algal blooms and related oceanographic features in the Bohai Sea during autumn 1998

Harmful algal blooms (HABs) are truly global marine phenomena of increasing significance. Some HAB occurrences are different to observe because of their high spatial and temporal variability and their advection, once formed, by surface currents. A serious HAB occurred in the Bohai Sea during autumn 1998, causing the largest fisheries economic loss. The present study analyzes the formation, distribution, and advection of HAB using satellite SeaWiFS ocean color data and other oceanographic data. The results show that the bloom originated in the western coastal waters of the Bohai Sea in early September, and developed southeastward when sea surface temperature (SST) increased to 25–26 °C. The bloom with a high Chl-a concentration (6.5 mg m⁻³) in center portion covered an area of 60 × 65 km². At the end of September, the bloom decayed when SST decreased to 22–23 °C. The HAB may have been initiated by a combination of the river discharge nutrients in the west coastal waters and the increase of SST; afterwards it may have been transported eastward by the local circulation that was enhanced by northwesterly winds in late September and early October.     

Temporal profiles of solar energetic particle events from SOHO/EPHIN data

Temporal profiles of energetic ions and electrons observed at 1 AU during solar energetic particle events are mainly determined by particle injection features, the observer location relative to the source region at the Sun, and the interplanetary space plasma and field conditions during particle transport. In this work, temporal profiles of 18 solar energetic particle events have been analyzed by fitting a pulse function to them in order to find a set of parameters which can be used to characterize the events.     

The lower ionosphere response to its disturbances by powerful radio waves

The paper presents data from some campaigns at Sura heating facility in 2011–1016. The experiments on probing of the artificial disturbed region of the lower ionosphere were carried out at two observation sites. One of them was located near Vasil’sursk 1 km from Sura facility (56.1°N; 46.1°E) and the other site was located at the Observatory (55.85°N; 48.8°E) of Kazan State University, 170 km to the East. Investigation of the features of the disturbed region of the lower ionosphere based on its diagnostics by the methods of the vertical sounding and oblique backscattering is the main goal of this paper. Ionosphere disturbance was fulfilled by the effect of the powerful radio wave of the ordinary or extraordinary polarization emitted by transmitters of the Sura facility with effective radiated power ERP = 50–120 MW at the frequency of 4.3, 4.7 and 5.6 MHz. Pumping waves were emitted with period from 30 s to 15 min. The disturbed region of the ionosphere in Vasil’sursk was probed by the vertical sounding technique using the partial reflexion radar at the frequency of 2.95 and 4.7 MHz. For the oblique sounding of the disturbed region the modified ionosonde Cyclon-M, operating at ten frequencies from 2.01 to 6.51 MHz was used at the Observatory site. On many heating sessions simultaneous variations of the probing partial reflection signals in Vasil’sursk and backscattered signals in Observatory were observed at the height at 40–100 km below the reflection height of the pumping wave. These observations were correlated with the pumping periods of the Sura facility. Possible mechanisms of the appearance of the disturbance in the lower ionosphere and its effect on the probing radio waves are discussed.     

Dependence of E-sporadic layer response on solar and geomagnetic activity variations from its ion composition

The aim of this paper is to consider an influence of solar and geomagnetic activity level variations on frequency and stochastic parameters of mid-latitude ionosphere sporadic-E layer. Critical frequency of sporadic-E layer, foEs, relative excess of sporadic ionization over monthly median values, foEs − foEm/foEm, and probability of Es layer appearance, PEs, are considered. It has been found that sporadic-E layer parameters response to solar and geomagnetic activity level variations can be both positive (foEs and PEs values are increase) and negative (foEs and PEs values are decrease). In particular, sporadic-E layers response to solar and geomagnetic activity variations are different depending upon layer intensity. It is suggested that revealed differences may be associated with dissimilarity of the layers ion composition (high-intensive layers are composed from metallic ions and low intensive composed from molecular ions).     

Atmospheric pressure variations at extratropical latitudes associated with Forbush decreases of galactic cosmic rays

Changes of troposphere pressure associated with short-time variations of galactic cosmic rays (GCRs) taking place in the Northern hemisphere’s cold months (October–March) were analyzed for the period 1980–2006, NCEP/NCAR reanalysis data being used. Noticeable pressure variations during Forbush decreases of GCRs were revealed at extratropical latitudes of both hemispheres. The maxima of pressure increase were observed on the 3rd–4th days after the event onsets over Northern Europe and the European part of Russia in the Northern hemisphere, as well as on the 4th–5th days over the eastern part of the South Atlantic opposite Queen Maud Land and over the d’Urville Sea in the Southern Ocean. According to the weather chart analysis, the observed pressure growth, as a rule, results from the weakening of cyclones and intensification of anticyclone development in these areas. The presented results suggest that cosmic ray variations may influence the evolution of extratropical baric systems and play an important role in solar-terrestrial relationships.     

Experimental and numerical simulation of super-pressure balloon apex section: Mechanical behavior in realistic flight conditions

In order to predict the flight parameters and to improve the life time of long duration super-pressure balloon, a research program on modelization and experimental simulation of the balloon envelope mechanical behavior is carried out.The balloon is a 10 m diameter type made with multilayer polymeric films and tapes. A facility was developed to measure, via a stereo-correlation system, the 3D displacement and the in plane solicitation strain of a 1.5 m diameter balloon envelope part in realistic flight conditions, i.e. pressure, temperature and loading at the sample boundaries. A time dependant non-linear Maxwell model of the polymeric material behavior was identified from uniaxial creep and relaxation tests and implemented in a Finite Elements code, simulating the sample tested in the facility. The Poisson ratio of the transparent and supple balloon film has been measured with an image correlation system.Experimental results are obtained both at room and in cold conditions (−60 °C) for various values of differential pressure.Vertical displacement and in plane 2D strain of apex part deduced from the numerical modeling are compared to experimental results.     

Solar rhythms in the characteristics of the Arctic frontal zone in the North Atlantic

Long-term changes of the Arctic frontal zone characteristics near the south-eastern coasts of Greenland were considered, the NCEP/NCAR reanalysis data being used. It was found that in the cold half of the year the temperature gradients in the layer 1000–500 hPa in the region under study reveal strong ∼10-yr and ∼22-yr periodicities that seem to be related to solar activity cycles. The results obtained suggest the influence of solar activity and cosmic ray variations on the structure of the temperature field of the troposphere resulting in the changes of the temperature contrasts in the Arctic frontal zone that, in turn, may affect the intensity of cyclogenesis at middle latitudes. The detected effects seem to indicate an important part of frontal zones in the mechanism of solar activity and cosmic ray variation influence on the development of extratropical baric systems. It is suggested that the variations of the temperature gradients revealed in the Arctic frontal zone are due to the radiative forcing of cloudiness changes which may be associated with geomagnetic activity and cosmic ray variations.     

An analysis of solar radio burst events on December 1, 2004

We present the analysis of the radio observations of December 1, 2004 from 07:00 UT to 07:40 UT in the 1.100–1.340 GHz band by Solar Broadband Radio Dynamic Spectrometer (SBRS) in Huairou Station. There are three groups of radio fine structures during the impulsive phase of this flare denoted by N1, Z2, and Z3. N1 has several emission lines with mixed fast and slow frequency drift rate which may reflect the conditions of flare loop and fast flows out from reconnection site; Z2 and Z3 are zebra patterns. The radio observations combined with hard X-ray and other observations show that the fine structures are connected with energetic particles. The information about magnetic field and energetic particle during the burst are also estimated based on our model.     

F2 layer characteristics and electrojet strength over an equatorial station

The data presented in this work describes the diurnal and seasonal variation in hmF2, NmF2, and the electrojet current strength over an African equatorial station during a period of low solar activity. The F2 region horizontal magnetic element H revealed that the Solar quiet Sq(H) daily variation rises from early morning period to maximum around local noon and falls to lower values towards evening. The F2 ionospheric current responsible for the magnetic field variations is inferred to build up at the early morning hours, attaining maximum strength around 1200 LT. The Sq variation across the entire months was higher during the daytime than nighttime. This is ascribed to the variability of the ionospheric parameters like conductivity and winds structure in this region. Seasonal daytime electrojet (EEJ) current strength for June solstice, March and September equinoxes, respectively had peak values ranging within 27–35 nT (at 1400 LT) , 30–40 nT (at 1200 LT) and 35–45 nT (at 1500 LT). The different peak periods of the EEJ strength were attributed to the combined effects of the peak electron density and electric field. Lastly, the EEJ strength was observed to be higher during the equinoxes than the solstice period.     

Storm time spatial variations in TEC during moderate geomagnetic storms in extremely low solar activity conditions (2007–2009) over Indian region

The total electron content (TEC) measurements from a network of GPS receivers were analyzed to investigate the storm time spatial response of ionosphere over the Indian longitude sector. The GPS receivers from the GPS Aided Geo Augmented Navigation (GAGAN) network which are uniquely located around the ∼77°E longitude are used in the present study so as to get the complete latitudinal coverage from the magnetic equator to low mid-latitude region. We have selected the most intense storms but of moderate intensity (−100 nT < Dst < −50 nT) which occurred during the unusually extremely low solar activity conditions in 2007–2009. Though the storms are of moderate intensity, their effects on equatorial to low mid-latitude ionosphere are found to be very severe as TEC deviations are more than 100% during all the storms studied. Interesting results in terms of spatial distribution of positive/negative effects during the main/early recovery phase of storms are noticed. The maximum effect was observed at crest region during two storms whereas another two storms had maximum effect near the low mid-latitude region. The associated mechanisms like equatorial electrodynamics and neutral dynamics are segregated and explained using the TIMED/GUVI and EEJ data during these storms. The TEC maps are generated to investigate the storm time development/inhibition of equatorial ionization anomaly (EIA).     

Airglow OH emission height inferred from the OH temperature and meteor trail diffusion coefficient

Simultaneous observations of the airglow OH(6,2) band rotational temperature, TOH, and meteor trail ambipolar diffusion coefficient, D, were carried out at Shigaraki (35°N, 136°E), during PSMOS 2003 Campaign, January 28 to February 8, 2003. The OH emission height was estimated by cross correlation analysis of the TOH and D nocturnal variations. A good correlation between TOH and D was obtained at 85 km of altitude. From the nocturnal variations of TOH and D, it is found that the OH emission peak height varied from 88 km before the midnight to 84 km in the early morning. The height variation could be caused by an atmospheric tidal effect in the emission height.     

On long-term changes of electron density in the upper mesosphere

Recent review study done jointly by 19 experts of 17 institutes shows zero trend of temperature in the upper mesosphere. In the light of this latest development, we have examined the long-term changes in electron density, [e], in this region. The study has been concentrated at 80 km. At this altitude, electrons are mainly produced by the interaction of nitric oxide, NO, by solar Ly-α. Any long-term change in this flux will affect trend of [e]. Considering this flux proportional to 10.7 cm solar flux, analysis of available 10.7 cm solar flux data from 1948 to 2003 has been made. A decreasing trend up to about 1970 and then an increasing trend are found. The over-all increasing trend of Ly-α flux during the past five decades is ∼0.17% per year. This increase also gives a ∼0.17% increasing trend per year in [e]. This non-anthropogenic increase is much less compared to the observed increase in [e] which is reported to be >0.7% per year. The observed increase in [e] of this magnitude will then, predominantly, be due to the anthropogenic effect. In zero trend in temperature, significant change in electron loss coefficient, α_eff, and [NO] are unlikely to take place to cause a significant change in [e]. The increase in [e] > 0.7% per year then can be explained by considering a decreasing trend in [O₂].     

Supernova remnant 1987A: The latest report from the Chandra X-ray Observatory

We continue monitoring supernova remnant (SNR) 1987A with the Chandra X-ray Observatory. As of 2004 January, bright X-ray spots in the northwest and the southwest are now evident in addition to the bright eastern ring. The overall X-ray spectrum, since 2002 December, can be described by a planar shock with an electron temperature of ∼2.1 keV. The soft X-ray flux is now 8 × 10⁻¹³ ergs cm⁻² s⁻¹, which is about five times higher than four years ago. This flux increase rate is consistent with our prediction based on an exponential density distribution along the radius of the SNR between the HII region and the inner ring. We still have no direct evidence of a central point source, and place an upper limit of L_X = 1.3 × 10³⁴ ergs s⁻¹ on the 3–10 keV band X-ray luminosity.     

The role of drift and convection–diffusion mechanisms in small energy global cosmic ray modulation: Application to the hysteresis phenomenon of proton and alpha particle satellite data

The hysteresis effect for small energies of galactic cosmic rays is due to two effects. The first is the same as for neutron monitor energies – the delay of the interplanetary processes responsible for cosmic ray modulation with respect to the initiating solar processes, according to the effective velocity of solar wind and shock waves propagation. Then, the observed cosmic ray intensity is connected to the solar activity variations during many months before the time of cosmic ray measurement. The second is caused by the time delay of small energy cosmic ray diffusion from the boundary of modulation region to the Earth’s orbit. The model describing the connection between solar activity variation and cosmic ray convection–diffusion global modulation for neutron monitor energies is here developed by taking into account also the time-lag of the small energy particle diffusion in the Heliosphere. We use theoretical results on drifts and analytically approximate the dependences of drifts from tilt angle, and take into account the dependence from the sign of primary particles, and from the sign of polar magnetic field (A > 0 or A < 0). The obtained results are applied on proton and alpha-particle satellite data. We analyze satellite 5-min data of proton fluxes with energies >1 MeV, >2 MeV, >5 MeV, >10 MeV, >30 MeV, >50 MeV, >60 MeV, >100 MeV, and in intervals 10–30 MeV, 30–60 MeV, and 60–100 MeV during January 1986–December 1999. We exclude periods with great cosmic ray increases caused by particle acceleration in solar flare events. Then, we determine monthly averaged fluxes, as well as 5-month and 11-month smoothed data. We analyze also satellite 5-min data on alpha-particle fluxes in the energy intervals 60-160 MeV, 160–260 MeV and 330–500 MeV during January 1986–May 2000. We correct observation data for drifts and then compare with what is expected according to the convection–diffusion mechanism. We assume different dimensions of the modulation region (by the time propagation X₀ of solar wind from the Sun to the boundary of modulation region), for X₀ values from 1 to 60 average months, by one-month steps. For each value of X₀ we determine the correlation coefficient between variations of expected and observed cosmic ray intensities (the estimation of cosmic ray intensities values is given in Section 3 by Eq. (9), and the determination of correlation and regression coefficients in Section 3 by Eq. (8)). The dimension of modulation region is determined by the value of X_0 max, for which the correlation coefficient reaches the maximum value. Then the effective radial diffusion coefficient and residual modulation in small energy region can be estimated.     

Regional tropospheric responses to long-term solar activity variations

The influence of ∼200-year solar activity variations (de Vries cyclicity) on climatic parameters has been analyzed. Analysis of palaeoclimatic data from different regions of the Earth for the last millennium has shown that ∼200-year variations in solar activity give rise to a pronounced climatic response. Owing to a nonlinear character of the processes in the atmosphere–ocean system and the inertia of this system, the climatic response to the global influence of solar activity variations has been found to have a regional character. The regions where the climatic response to long-term solar activity variations is stable and the regions where the climatic response is unstable, both in time and space, have been revealed. It has also been found that a considerable lag of the climatic response and reversal of its sign with respect to the solar signal can occur. Comparison of the obtained results with the simulation predictions of the atmosphere–ocean system response to long-term solar irradiance variations (T > 40 years) has shown that there is a good agreement between experimental and simulation results.     

Detection of long-term trends in the mesosphere/lower thermosphere from ground-based radio propagation measurements

Mesospheric temperature trends can be derived from LF phase-height observations in mid-latitudes supported by ionospheric absorption and ionosonde observations. Analysing the full observation period from 1959 until 2003, a mean yearly temperature trend has been derived with −0.25 K/yr for the height interval from 48 to 82 km. Subdividing the whole observation interval in two parts before and after 1979, the trend is markedly stronger in the second period with −0.38 K/yr compared with −0.20 K/yr in the first part before 1979. These differences can at least partly be explained by a steeper CO₂ increase and ozone decrease in the second interval. The differences in the mesospheric temperature trends are most evidently expressed during winter months and are markedly smaller during summer season. The reason of this seasonal difference is not quite clear; it may be related with detected ozone trends which are clearly stronger during winter months on both hemispheres.     

Global MGS TES data and Mars-GRAM validation

Mars Global Reference Atmospheric Model (Mars-GRAM 2001) is an engineering-level Mars atmosphere model widely used for many Mars mission applications. From 0 to 80 km, it is based on NASA Ames Mars General Circulation Model (MGCM), while above 80 km it is based on University of Michigan Mars Thermospheric General Circulation Model. Mars-GRAM 2001 and MGCM use surface topography from Mars Global Surveyor Mars Orbiter Laser Altimeter (MOLA). Validation studies are described comparing Mars-GRAM with a global summary data set of Mars Global Surveyor Thermal Emission Spectrometer (TES) data. TES averages and standard deviations were assembled from binned TES data which covered surface to ∼40 km, over more than a full Mars year (February 1999–June 2001, just before start of a Mars global dust storm). TES data were binned in 10° × 10° latitude–longitude bins (36 longitude bins, centered at 5°–355°, by 18 latitude bins, centered at −85° to +85°), and 12 seasonal bins (based on 30° increments of L_s angle). Bin averages and standard deviations were assembled at 23 data levels (temperature at 21 pressure levels, plus surface temperature and surface pressure). Two time-of-day bins were used: local time near 2 or 14 h. Two dust optical depth bins were used: infrared optical depth, either less than or greater than 0.25 (which corresponds to visible optical depth less than or greater than about 0.5). For interests in aerocapture and precision entry and landing, comparisons focused on atmospheric density. TES densities versus height were computed from TES temperature versus pressure, using assumptions of perfect gas law and hydrostatics. Mars-GRAM validation studies used density ratio (TES/Mars-GRAM) evaluated at data bin center points in space and time. Observed average TES/Mars-GRAM density ratios were generally 1 ± 0.05, except at high altitudes (15–30 km, depending on season) and high latitudes (>45°N), or at most altitudes in the southern hemisphere at L_s ∼ 90° and 180°. Compared to TES averages for a given latitude and season, TES data had average density standard deviation about the mean of ∼2.5% for all data, or ∼1–4%, depending on time of day and dust optical depth. Average standard deviation of TES/Mars-GRAM density ratio was 8.9% for local time 2 h and 7.1% for local time 14 h. Thus standard deviation of observed TES/Mars-GRAM density ratio, evaluated at matching positions and times, is about three times the standard deviation of TES data about the TES mean value at a given position and season.