Satellite monitoring of anomalous effects in the ionosphere probably related to strong earthquakes

The paper presents simultaneously a satellite monitoring of plasma density disturbances and quasi-static electric fields on the one hand, and a ground earthquake monitoring, on the other. The INTERCOSMOS-BULGARIA-1300 satellite was launched on 7th August 1981. It had a perigee of 825 km, an apogee of 906 km and orbit inclination of 81.2°. For satellite’s orbits over sources of earthquakes M ⩾ 5.5 (seismic data of the World Data Center – Denver, Colorado), in the time period 12.08.1981–30.12.1981, some ionospheric effects are investigated. These effects can be considered eventually as possible earthquake precursors. The paper discusses specific anomalous effects observed in the ionosphere, which cannot be explained by factors of solar-magnetospheric origin. They could be associated with the earthquake growing processes in the lithosphere and troposphere zones above the source. In conclusion the results of ionosphere anomalous effects monitoring are proposed.     

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.     

Spectral filter for signal identification in the kHz SLR measurements of the fast spinning satellite Ajisai

The high repetition rate satellite laser ranging (SLR) measurements to the fast spinning satellites contain a frequency signal caused by the rotational motion of the corner cube reflector (CCR) array. The spectral filter, developed here, is based on the Lomb algorithm, and is tested with the simulated and the observed high repetition rate SLR data of the geodetic satellite Ajisai (spin period ∼2 s). The filter allows for the noise elimination from the SLR data, and for identification of the returns from the single CCRs of the array – even for the low return rates. Applying the spectral filter to the simulated SLR data increases the S/N ratio by a factor 40–45% for all return rates. Filtering out the noise from the observed data strengthens the frequency signal by factor of ∼25 for the low return rates, which significantly helps to determine the spin phase of the satellite. The spectral filter is applied to the Graz SLR data and the spin rates of Ajisai are determined by two different methods: the frequency analysis and the phase determination of the spinning retroreflector array.The analysis of more than 8 years of the Graz SLR measurements indicates an exponential spin rate trend: f = 0.67034 exp(−0.0148542 Y) [Hz], RMS = 0.085 mHz, where Y is the year since launch. The highly accurate spin rate information demonstrates periodic changes related to the precession of the orbital plane of Ajisai, as it determines the amount of energy received by the satellite from the Sun. The rate of deceleration of Ajisai is not constant: the half life period of the satellite’s spin oscillates around 46.7 years with an amplitude of about 5 years.     

Use of the ASTER satellite images for evaluation of structural changes in the Popocatépetl volcano related to microseismicity

Mexico is one of the most volcanically active regions in North America. Volcanic activity in central Mexico is associated with the subduction of the Cocos and Rivera plates beneath the North American plate. Periods of enhanced microseismic activity, associated with the volcanic activity of the Popocatépetl volcano are compared with periods, during which the microseismic activity was low. We detected systematical changes in the number of lineaments, associated with the microseismic activity due to lineament analysis of a temporal sequence of high resolution satellite images of the Popocatépetl volcano, provided by the ASTER/VNIR instrument. The Lineament Extraction and Stripes Statistic Analysis (LESSA) software package was used for the lineament extraction. In the future it would allow develop a methodology for detection of possible elevation of pressure in volcano edifice.     

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.     

On the mid-latitude ionospheric storm association with intense geomagnetic storms

The bulk association between ionospheric storms and geomagnetic storms has been studied. Hemispheric features of seasonal variation of ionospheric storms in the mid-latitude were also investigated. 188 intense geomagnetic storms (Dst ≤ 100 nT) that occurred during solar cycles 22 and 23 were considered, of which 143 were observed to be identified with an ionospheric storm. Individual ionospheric storms were identified as maximum deviations of the F2 layer peak electron density from quiet time values. Only ionospheric storms that could clearly be associated with the peak of a geomagnetic storm were considered. Data from two mid-latitude ionosonde stations; one in the northern hemisphere (i.e. Moscow) and the other in the southern hemisphere (Grahamstown) were used to study ionospheric conditions at the time of the individual geomagnetic storms. Results show hemispheric and latitudinal differences in the intensity and nature of ionospheric storms association with different types of geomagnetic storms. These results are significant for our present understanding of the mechanisms which drive the changes in electron density during different types of ionospheric storms.     

Vegetation biomass change of the Bosoh Peninsula: Impacted by the volcano fumes from the Miyakejima

Since the 2001 eruption, the volcanic fumes from the Oyama of Miyakejima Island, located approximately 200 km southwest from the mainland of Japan, have affected the vegetation health and biomass of the Tokyo metropolitan area. In this study, we evaluated the potential for measuring/mapping relative forest damage/recovery in forests in the Bosoh Peninsula, Chiba Prefecture, Japan, using LANDSAT Enhanced Thematic Mapper data and ASTER Level-1B data. The simple ratio which is derived from the near-infrared and red reflective response was found to correlate well with ground-based measurements of forest damage caused by continuous SO₂ contained in the volcanic fumes. Accumulated data by the synchronized field campaign observation of the Japanese cedar field spectra and by satellite analysis indicates the following interesting features: (1) regional differences of SO₂ levels in the forestal area can be estimated by the simple spectrum ratio of the near-infrared and red bands of multispectral satellite data; (2) time-series images using the simple ratios indicated the location of still damaged and recovered forests between the 2001 eruption and 2003; (3) if the vegetations under study are homogeneous single-species with a similar leaf area index and a similar age, the simple ratio was useful to evaluate vegetation damage or recovery among individual forests throughout many of the forests in the Bosoh Peninsula; (4) the damage level of the Tokyo Bay side was found to be consistently higher than the Pacific Ocean side of the peninsula. This obvious difference is likely caused, not only by SO₂ from Miyakejima Island, but also by long term and chronic discharged gas from heavy industrial complexes.     

A search for the signature of microquasars in the cosmic ray iron spectrum measured by TIGER

In a recent paper Heinz and Sunyaev suggest that relativistic jets observed in microquasars might result in narrow features in the energy spectra of heavy cosmic rays with ≈1 to ≈10 GeV/nuc. They further argue that such features might be observable if there has been one or more microquasars nearby within the last few million years. We report preliminary results of a search for evidence of such features using data from a 32-day balloon flight of the Trans-Iron Galactic Element Recorder (TIGER). Although this flight took place near solar maximum, calculations of the broadening effects of solar modulation indicate that a narrow feature of sufficient intensity should still be observable. An energy spectrum for iron with high statistical significance has been derived from ≈100,000 Fe events in the energy range from about 2.5 to 10 GeV/nuc. Although our preliminary results do not reveal any obvious features, we will discuss the possibility of observing such features with TIGER and other instruments.     

Study of >100 GeV electrons with BETS detector using a long duration balloon flight in Antarctica

We have observed cosmic-ray electrons from 10 to 1000 GeV by a long duration balloon flight using Polar Patrol Balloon (PPB) in Antarctica. The observation was carried out for 13 days at an altitude of 35 km in January 2004. The detector is an imaging calorimeter composed of scintillating-fiber belts and plastic scintillators inserted between lead plates. The geometrical factor of detector is about 600 cm²sr and the total thickness of lead absorber is 9 radiation lengths. The performance of the detector has been confirmed by the CERN-SPS beam test and also investigated by Monte-Carlo simulations. New telemetry system using a commercial satellite of iridium, power supply by solar batteries, and automatic level control using CPU have successfully been developed and operated during the flight. We have collected 5.7 × 10³ events over 100 GeV including nearly 100 candidates of primary electrons.     

VLBI observations to the APOD satellite

The APOD (Atmospheric density detection and Precise Orbit Determination) is the first LEO (Low Earth Orbit) satellite in orbit co-located with a dual-frequency GNSS (GPS/BD) receiver, an SLR reflector, and a VLBI X/S dual band beacon. From the overlap statistics between consecutive solution arcs and the independent validation by SLR measurements, the orbit position deviation was below 10?cm before the on-board GNSS receiver got partially operational. In this paper, the focus is on the VLBI observations to the LEO satellite from multiple geodetic VLBI radio telescopes, since this is the first implementation of a dedicated VLBI transmitter in low Earth orbit. The practical problems of tracking a fast moving spacecraft with current VLBI ground infrastructure were solved and strong interferometric fringes were obtained by cross-correlation of APOD carrier and DOR (Differential One-way Ranging) signals. The precision in X-band time delay derived from 0.1?s integration time of the correlator output is on the level of 0.1?ns. The APOD observations demonstrate encouraging prospects of co-location of multiple space geodetic techniques in space, as a first prototype.     

Response of equatorial,low- and mid-latitude F-region in the American sector during the intense geomagnetic storm on 24–25 October 2011

In this paper, we present and discuss the response of the ionospheric F-region in the American sector during the intense geomagnetic storm which occurred on 24–25 October 2011. In this investigation ionospheric sounding data obtained of 23, 24, 25, and 26 October 2011 at Puerto Rico (United States), Jicamarca (Peru), Palmas, São José dos Campos (Brazil), and Port Stanley, are presented. Also, the GPS observations obtained at 12 stations in the equatorial, low-, mid- and high-mid-latitude regions in the American sector are presented. During the fast decrease of Dst (about ∼54 nT/h between 23:00 and 01:00 UT) on the night of 24–25 October (main phase), there is a prompt penetration of electric field of magnetospheric origin resulting an unusual uplifting of the F region at equatorial stations. On the night of 24–25 October 2011 (recovery phase) equatorial, low- and mid-latitude stations show h′F variations much larger than the average variations possibly associated with traveling ionospheric disturbances (TIDs) caused by Joule heating at high latitudes. The foF2 variations at mid-latitude stations and the GPS-VTEC observations at mid- and low-latitude stations show a positive ionospheric storm on the night of 24–25 October, possibly due to changes in the large-scale wind circulation. The foF2 observations at mid-latitude station and the GPS-VTEC observations at mid- and high-mid-latitude stations show a negative ionospheric storm on the night of 24–25 October, probably associated with an increase in the density of molecular nitrogen. During the daytime on 25 October, the variations in foF2 at mid-latitude stations show large negative ionospheric storm, possibly due to changes in the O/N2 ratio. On the night of 24–25, ionospheric plasma bubbles (equatorial irregularities that extended to the low- and mid-latitude regions) are observed at equatorial, low- and mid-latitude stations. Also, on the night of 25–26, ionospheric plasma bubbles are observed at equatorial and low-latitude regions.     

Possible DC electric field in the ionosphere related to seismicity

This paper presents the method for calculation of DC electric field in the atmosphere and the ionosphere generated by model distribution of external electric current in the lower atmosphere. Appearance of such current is associated with enhancement of seismic activity that is accompanied by emanation of soil gases into the atmosphere. These gases transfer positive and negative charged aerosols. Atmospheric convection of charged aerosols forms external electric current, which works as a source of conductivity current in the atmosphere–ionosphere electric circuit. It is shown that DC electric field generated in the ionosphere by this current reaches up to 10 mV/m, while the long-term vertical electric field disturbances excited near the Earth surface do not exceed 100 V/m. Such limitation of the near-ground field is caused by the formation of potential barrier for charged particles at the Earth surface in a process of their transport from soil to atmosphere.     

Observed flux density enhancement at submillimeter wavelengths during an X-class flare

We analyse the 30 October, 2004, X1.2/SF solar event that occurred in AR 10691 (N13 W18) at around 11:44 UT. Observations at 212 and 405 GHz of the Solar Submillimeter Telescope (SST), with high time resolution (5 ms), show an intense impulsive burst followed by a long-lasting thermal phase. EUV images from the Extreme Ultraviolet Imaging Telescope (SOHO/EIT) are used to identify the possible emitting sources. Data from the Radio Solar Telescope Network (RSTN) complement our spectral observations below 15 GHz. During the impulsive phase the turnover frequency is above 15.4 GHz. The long-lasting phase is analysed in terms of thermal emission and compared with GOES observations. From the ratio between the two GOES soft X-ray bands, we derive the temperature and emission measure, which is used to estimate the free-free submillimeter flux density. Good temporal agreement is found between the estimated and observed profiles, however the former is larger than the latter.     

Automatic lineament extraction in a heavily vegetated region using Landsat Enhanced Thematic Mapper (ETM+) imagery

Lineament extraction from satellite remotely sensed data has been one of the widely used applications of remote sensing in geology. In fact, recent advances in digital image processing allow such lineament extraction to be accomplished in semi-automatic to fully automatic approaches. However, satellite remotely sensed data acquired in heavily vegetated regions such as tropical rainforest, are vulnerable to higher inherent noise levels attributed to the resultant effects of scattering by clouds and adjacency effects of highly inhomogeneous vegetation cover within the pixel dimension. In this study, we examined the effects of noise levels to lineament extraction using a fully automatic approach, consisting of a combination of edge-line detection algorithms. Ancillary information from a digitized topographic map and image classification was used to discriminate between cultural and natural lineaments from the extracted lineaments. Adapting the combination of edge detection and a line-linking algorithm, we have found the optimal parameters for automatic lineament extraction of such complex areas using Enhanced Thematic Mapper (ETM+) data. A noise level of 30% is the maximum threshold before artifacts are generated. It is therefore concluded that the combination of edge-based and line-linking digital image processing operations with the priori local optimal parameters is crucial in lineament feature extraction in heavily vegetated regions.     

Single event upset investigations on the “Flying Laptop” satellite mission

The radiation effects in electronic parts are called single-event effects, which are deemed to be critical for space missions. This paper presents the Single Event Upsets that were observed in an onboard memory device of the Low Earth Orbit “Flying Laptop” satellite mission during its in-orbit operation. The Single Event Upsets were carefully mapped on the satellite orbital space itself and their root causes were investigated together with their rates of occurrence. Subsequently, the events were traced to show several root cause sources such as (i) trapped energetic protons leaking to low altitudes within the South Atlantic Anomaly, (ii) Solar Energetic Particles emitted by an impulsive event on 10 September 2017, and (iii) Galactic Cosmic Rays. A profound analysis was carried out on the observed flight data, and its corresponding results are actually in agreement with the standard energetic particle models. The presented results provide another important insight on the Single Event Upsets for future Low Earth Orbit satellite missions.     

Ionospheric response to the storm-time disturbance of 29 May, 2010

The Ionospheric F2-layer peak parameters response to a magnetic storm had been investigated over Ilorin, Nigeria (Lat. 8:53°N, Long. 4.5°E, dip angle, −2.96°), Jicamarca, Peru (11.95°S, 76.87°W, dip angle, 0.8°) and Hermanus, South Africa (34.42°S, 19.22°E, dip angle, −60.77°), using percentage enhancement/depletion values. Our results showed an enhancement in NmF2 at all of these stations. Averagely, pre-noon and post-noon peaks are highest at Ilorin during quiet time. The similar pattern observed for quiet condition between Ilorin and Jicamarca was due to their latitudinal positions. For disturbed NmF2 condition, Jicamarca and Ilorin recorded higher peaks at nighttime than during the daytime for the storms main phase, and the reverse over Hermanus. The nighttime and daytime increases were observed respectively at Ilorin and Hermanus during the recovery period. The hmF2 variation recorded higher enhancement at Jicamarca during the daytime and at Hermanus at nighttime during the main phase. During the recovery phase, the highest enhancement was recorded during the daytime at Jicamarca, and over Hermanus at nighttime. These observations find their explanation in the magnetospheric current, solar wind and E × B drift.     

Evidence for a solar coronal thick-target hard X-ray source observed by RHESSI

We study a solar flare hard X-ray (HXR) source observed by the Reuven Ramaty high energy solar spectroscopic imager (RHESSI) in which the HXR emission is almost entirely in a coronal loop so dense as to be collisionally thick at electron energies up to ∼45−60 keV. This contrasts with most events previously reported in which the HXR emission is primarily from the loop footpoints in the collisionally dense chromosphere. In particular, we show that the high loop column densities inferred from the GOES and RHESSI soft X-ray emission measure and the volume of the flare loop are consistent with the coronal thick-target interpretation of the HXR images and spectra. The high column densities observed already at the very beginning of the impulsive phase are explained by chromospheric evaporation during a preflare which, as Nobeyama 17 GHz radio images reveal, took place in the same set of nested loops as the main flare.     

Cyclone regeneration in the North Atlantic intensified by energetic solar proton events

Atmospheric effects of energetic solar proton events (SPE) were studied in the North Atlantic region, for particle energies above 90 MeV, using NCEP/NCER reanalysis data and weather charts. A significant lowering of the pressure levels in the troposphere accompanied by an increase of the cyclonic vorticity was found near the south-eastern coast of Greenland on days following the event onsets. According to the weather charts, the detected effects are caused by the re-deepening (the regeneration) of well developed cyclones that seems to be intensified during the SPE under study. A joint analysis of the pressure and temperature variations showed a noticeable decrease of the temperature in the rear of the deepening cyclones that may be due to the cold advection increase. The results obtained suggest the influence of energetic SPE on the cyclone development as well as the importance of the frontal zone situated near the Greenland coast for this influence. The physical mechanism may involve the increase of cold advection due to changes in the temperature gradients in this region, resulting from radiative forcing and/or latent heat release related to variations of cloudiness.     

The temperature map of the loop and coronal sources of an occulted solar flare

An occulted solar flare occurred at about 06:07 UT on 2002, November 2. The RHESSI X-ray images show two separate parts. The lower part consists of a complete loop and the upper part a coronal source which well extends above the solar limb. The loop source shrank for about 3 min with a speed of ∼24 km s⁻¹ during the early impulsive phase and then expanded at ∼7 km s⁻¹, while the coronal source presented an upward motion at about 6 km s⁻¹. We obtained the temperature map of the loop source from RHESSI image spectrum. The temperature of the loop increases with altitude, indicating that the reconnection X-point of this flare is located above the loop source. However, the apparent coronal source is the top of another independent large-scale loop.     

The plerionic supernova remnant G21.5-0.9: In and out

The absence of a supernova remnant (SNR) shell surrounding the Crab and other plerions (pulsar wind nebulae) has been a mystery for three decades. G21.5-0.9 is a particularly intriguing plerionic SNR in which the central powering engine is not yet detected. Early CHANDRA observations revealed a faint extended X-ray halo which was suggested to be associated with the SNR shell; however its spectrum was non-thermal, unlike what is expected from an SNR shell. On the other hand, a plerionic origin to the halo is problematic since the X-ray plerion would be larger than the radio plerion. We present here our analysis of an integrated 245 ks of archival CHANDRA data acquired with the High-Resolution Camera (HRC) and 520 ks acquired with the Advanced CCD Imaging Spectrometer (ACIS). This study provides the deepest and highest resolution images obtained to date. The resulting images reveal for the first time: (1) a limb-brightened morphology in the eastern section of the halo, and (2) a rich structure in the inner (40″-radius) bright plerion including wisps and a double-lobed morphology with an axis of symmetry running in the northwest–southeast direction. Our spatially resolved spectroscopic study of the ACIS-I data indicates that the photon index steepens with increasing distance from the central point source out to a radius of 40″ then becomes constant at ∼2.4 in the X-ray halo (for a column density N_H = 2.2 × 10²² cm⁻²). No line emission was found from the eastern limb; however marginal evidence for line emission in the halo’s northern knots was found. This study illustrates the need for deep CHANDRA observations to reveal the missing SNR material in Crab-like plerions.