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.     

Giant radio pulses from the Crab pulsar

Individual giant radio pulses (GRPs) from the Crab pulsar last only a few microseconds. However, during that time they rank among the brightest objects in the radio sky reaching peak flux densities of up to 1500 Jy even at high radio frequencies. Our observations show that GRPs can be found in all phases of ordinary radio emission including the two high frequency components (HFCs) visible only between 5 and 9 GHz [Moffett, D.A., Hankins, T.H. Multifrequency radio observations of the Crab pulsar. Astrophys. J. 468, 779–783, 1996]. This leads us to believe that there is no difference in the emission mechanism of the main pulse (MP), inter pulse (IP) and HFCs. High resolution dynamic spectra from our recent observations of giant pulses with the Effelsberg telescope at a center frequency of 8.35 GHz show distinct spectral maxima within our observational bandwidth of 500 MHz for individual pulses. Their narrow band components appear to be brighter at higher frequencies (8.6 GHz) than at lower ones (8.1 GHz). Moreover, there is an evidence for spectral evolution within and between those structures. High frequency features occur earlier than low frequency ones. Strong plasma turbulence might be a feasible mechanism for the creation of the high energy densities of ∼6.7 × 10⁴ erg cm⁻³ and brightness temperatures of ∼10³¹ K.     

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.     

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).     

The coronal loop controversy: TRACE analysis

The temperature distribution along coronal loops provides an important clue for solving the coronal heating problem. Recent analysis, however, has produced conflicting results. Here, we analyze in detail one component of this analysis – the effect of background subtraction on the temperature of loops observed with the Transition Region and Coronal Explorer (TRACE). Specifically, we selected 10 coronal loops that were visible in the TRACE 171 Å and 195 Å passbands. We chose between 20 and 30 pixel along each loop and background pixels to correspond with the loop pixels. Temperature analysis was done in three different ways: (1) standard TRACE analysis of the loop pixels with no background subtraction; (2) constant background subtraction for each TRACE image; (3) pixel pair background subtraction. Each method produced a temperature estimate for the selected pixels. We find that a flat line is an excellent fit to the temperature results – the analysis indicates that the temperature of the loop is uniform along the length visible by TRACE. However, if we select random pixels and plot the temperature results in the same way, these pixels indicate that the temperature of this “structure” is also uniform. We conclude therefore, that in the cases considered here, the image ratio analysis does not produce a physically meaningful value of plasma temperature; in addition, background subtraction makes no significant difference to the temperatures results.     

The total ozone and UV solar radiation over Stara Zagora,Bulgaria

The results from direct ground-based solar UV irradiance measurements and the total ozone content (TOC) over Stara Zagora (42° 25′N, 25° 37′E), Bulgaria are presented. During the period 1999–2003 the TOC data show seasonal variations, typical for the middle latitudes – maximum in the spring and minimum in the autumn. The comparison between TOC ground-based data and Global Ozone Monitoring Experiment (GOME) satellite-borne ones shows a seasonal dependence of the differences between them.A strong negative relationship between the total ozone and the 305 nm wavelength irradiance was found. The dependence between the two variables is significant (r = −0.62 ± 0.18) at 98% confidence level.The direct sun UV doses for some specific biological effects (erythema and eyes) are obtained. The estimation of the radiation amplification factor RAF shows that the ozone reduction by 1% increases the erythemal dose by 2.3%. The eye-damaging doses are more influenced by the TOC changes and in this case RAF = −2.7%.The amount of these biological doses depended on the solar altitude over the horizon. This dependence was not so strong when the total ozone content in the atmosphere was lower.     

Applicationof remote sensing techniques in geological studies of various mineralized belts in India. (National paper on IGCP — 143)

Geological interpretation, analysis and evaluation of LANDSAT and airborne MSS data, panchromatic and colour infrred photographs and spectral reflectance data for a number of test sites in different mineralized belts in India have led to certain significant results. Lithological discrimination on digitally processed enhanced MSS data products has helped in further subdivisions of major rock groups. Lineament pattern analyses reveals three prominent sets of fractures in Indian Precambrian Shield with predominant ENE-WSW megalineament probably representing earliest deep crustal fractures and subsequent NNW-SSE and NE-SW fractures, associated with polyphase tectonic movements. Correlation of structures and mineralization indicate structural control for most of the mineralized belts with intersecting locii type concentration of ore deposits in some cases. In a few cases control of mineralization is lithology and stratigraphy. Spectral groundtruth and laboratory studies indicate that major rock types have characteristic spectral patterns (signatures) which would be useful for lithological mapping by automatic classification techniques. The work also has led to the improvement of enhanced MSS data products for geological studies.     

Spectral analysis of a two-ribbon microflare

Using high-resolution Hα, CaII 8542 Å and FeI 6302.5 Å Stokes spectral data obtained simultaneously with THEMIS in 2002 September, we have analyzed the spectra and the characteristics of a two-ribbon microflare (MF). The hard X-ray emission provides evidence of non-thermal particle acceleration in the microflare. The two-ribbons are located on either sides of the magnetic polarity inversion line. The non-thermal characteristics mainly appeared at the outer edges of the flare ribbons. It indicates that the instantaneous magnetic reconnection and the particle acceleration mainly took place at the outer edges of the flare ribbons. Using the Hα and CaII 8542 Å line profiles and the non-LTE calculation, we obtain the semi-empirical atmospheric model for the bright kernel of the MF. The result indicates that the temperature enhancement in the chromosphere is about 2000–2500 K.     

The hard X-ray emission from the Vela pulsar wind nebula

We present the results of a preliminary spectral analysis performed on the BeppoSAX and XMM observations of the Vela plerion. The broad energy range covered by the instruments on board the two observatories allows an evaluation of the spectral parameters of the high energy emission model and provides an indication on the morphology of the source emission above 10 keV. We confirm the softening of the PWN spectrum (3–10 keV band) at distances greater than 4′ from the pulsar and estimate the diameter of the high energy (>10 keV) emission region to be on the order of 25′–30′.     

Hard X-ray component in Sco X-1 spectrum: Synchrotron emission from a nano quasar

Sco X-1 is a low mass X-ray binary system and with the recent observations of a resolved radio jet, the source has been included in the list of galactic microquasars. The observed spectral data in the 2–20 keV energy band fits a thermal emission. Above 20 keV, a hard tail has been reported on occasions. During our continuing balloon borne X-ray survey in the 20–200 keV region using high sensitivity Large Area Scintillation counter Experiment, Sco X-1 was observed on two different occasions. Even though the total X-ray luminosity of the source was different, the spectral nature of the source did not show any variation. The presence of hard X-ray flux is unmistakable. We present the spectral and temporal data in the hard X-ray band and discuss the results in terms of geometrical characteristics of X-ray source and its observed temporal properties. We note that the jet activity is similar to the microquasars, however, the absence of the large magnitude abrupt changes in X-ray light curve compared to GRS1915 + 105 suggest that the quasar-like behaviour is at a nano scale.     

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.     

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.     

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.     

Detection of meteors and sub-relativistic dust grains by the fluorescence detectors of ultra high energy cosmic rays

Fluorescence detectors of ultra high energy cosmic rays (UHECR) allow to record not only the extensive air showers, initiated by the UHECR particles, but also to detect light, produced by meteors and by the fast dust grains. It is shown that the fluorescence detector operated at the mountain site can register signals from meteors with kinetic energy threshold of about 25 J (meteor mass ∼ 5 × 10⁻⁶ g, velocity ∼ 3 × 10⁶ cm/s). The same detector might be used for recording of the dust grains of smaller mass (of about 10⁻¹⁰ g) but with velocity 10⁹ cm/s, close to the light velocity (sub-relativistic dust grains). The light signal from a sub-relativistic dust grain is expected in much shorter time scale (∼0.001 s), in comparison with the meteor signal (∼0.1–1 s), and much longer than duration of the UHECR signals (tens of μs). The fluorescence detector capable to register various phenomena: from meteors to UHECR – should have a variable pixel and selecting system integration time. A study of the new phenomenon of sub-relativistic grains will help to understand the mechanism of particle and dust grain acceleration in astrophysical objects (in SN explosions, for example).     

Comparison of vertical E × B drift velocities and ground-based magnetometer observations of DELTA H in the low latitude under geomagnetically disturbed conditions

In the present work, we analyzed the daytime vertical E × B drift velocities obtained from Jicamarca Unattended Long-term Ionosphere Atmosphere (JULIA) radar and ΔH component of geomagnetic field measured as the difference between the magnitudes of the horizontal (H) components between two magnetometers deployed at two different locations Jicamarca, and Piura in Peru for 22 geomagnetically disturbed events in which either SC has occurred or Dst_max < ?50 nT during the period 2006–2011. The ΔH component of geomagnetic field is measured as the differences in the magnitudes of horizontal H component between magnetometer placed directly on the magnetic equator and one displaced 6–9° away. It will provide a direct measure of the daytime electrojet current, due to the eastward electric field. This will in turn gives the magnitude of vertical E × B drift velocity in the F region. A positive correlation exists between peak values of daytime vertical E × B drift velocity and peak value of ΔH for the three consecutive days of the events. It was observed that 45% of the events have daytime vertical E × B drift velocity peak in the magnitude range 10–20 m/s and 20–30 m/s and 20% have peak ΔH in the magnitude range 50–60 nT and 80–90 nT. It was observed that the time of occurrence of the peak value of both the vertical E × B drift velocity and the ΔH have a maximum (40%) probability in the same time range 11:00–13:00 LT. We also investigated the correlation between E × B drift velocity and D_st index and the correlation between delta H and D_st index. A strong positive correlation is found between E × B drift and D_st index as well as between delta H and D_st Index. Three different techniques of data analysis – linear, polynomial (order 2), and polynomial (order 3) regression analysis were considered. The regression parameters in all the three cases were calculated using the Least Square Method (LSM), using the daytime vertical E × B drift velocity and ΔH. A formula was developed which indicates the relationship between daytime vertical E × B drift velocity and ΔH, for the disturbed periods. The E × B drift velocity was then evaluated using the formulae thus found for the three regression analysis and validated for the ‘disturbed periods’ of 3 selected events. The E × B drift velocities estimated by the three regression analysis have a fairly good agreement with JULIA radar observed values under different seasons and solar activity conditions. Root Mean Square (RMS) errors calculated for each case suggest that polynomial (order 3) regression analysis provides a better agreement with the observations from among the three.     

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.     

Propagation of cosmic-ray electrons in the Galaxy

We formulate the global propagation model of cosmic-ray electrons including the source region, which is currently considered to be supernova remnants (SNRs). The model is characterized by the escape rate of electrons from SNRs into the interstellar space. It becomes clear that the energy index of the escape rate influences the high energy side of the interstellar spectrum and makes it possible to explain the observed data up to 2 TeV in the case of source spectral index smaller than 2.2 that is expected from the radio spectrum in SNRs. The escape lifetime of electrons in SNRs is also discussed by using the ratio of the radio flux in two regions: SNRs and the Galaxy. The result shows the mean lifetime in SNRs of ∼10⁴ yr around 1 GeV, which corresponds to the SNR age in the Sedov phase.     

Progress in laboratory performance verification of the R-process Isotope Observer

The source of galactic cosmic rays (GCR) remains unknown, despite many decades of research. It is generally accepted among high-energy astrophysicists that GCR are accelerated by supernova (SN) shocks in the interstellar medium, but the evidence for this remains indirect. However, the fact that the cosmic-ray spectrum continues far past the limit of SN shock acceleration is a serious problem for this model. Most SNae occur in superbubbles – regions of the interstellar medium containing massive stars and SN remnants. Thus, SNae should be accelerating the ejecta of other SNae. The local medium of the superbubble will be enhanced in freshly synthesized R-process material. We are currently developing the R-process Isotope Observer (RIO), which could make the first measurements of the isotopic abundances of the “ultraheavy” GCR (those in the range 32 ⩽ Z ⩽ 42) in order to determine the fractional contribution of freshly synthesized R-process material in GCR through the measurement of several key isotopic ratios. We have recently performed a series of tests of the RIO detector with 400 A MeV ⁸⁴Kr at the HIMAC accelerator in Japan. The results of these tests indicate that a mass resolution of ∼0.3 u can be achieved with RIO. We also discuss the prospects for accommodation of the RIO mission.     

Nonlinear processes in cosmic-ray precursor of strong supernova shock: Maximum energy and average energy spectrum of accelerated particles

The instability in the cosmic-ray precursor of a supernova shock is studied. The level of turbulence in this region determines the maximum energy of accelerated particles. The consideration is not limited by the case of weak turbulence. It is assumed that the Kolmogorov type nonlinear wave interactions together with the ion-neutral collisions restrict the amplitude of random magnetic field. As a result, the maximum energy of accelerated particles strongly depends on the age of a SNR. The average spectrum of cosmic rays injected in the interstellar medium in the course of adiabatic SNR evolution takes the approximate form E⁻² at energies larger than 10–30 GeV/nucleon with the maximum energy that is close to the position of the knee in cosmic-ray spectrum at 4 × 10¹⁵ eV. At an earlier stage of SNR evolution – the ejecta-dominated stage, the particles are accelerated to higher energies and have a rather steep power-law distribution. These results suggest that the knee may mark the transition from the ejecta-dominated to the adiabatic evolution of SNR shocks which accelerate cosmic rays.     

Signatures of magnetic reconnection in solar radio observations?

Magnetic reconnection occurs during eruptive processes (flares, CMEs) in the solar corona. This leads to a change of magnetic connectivity. Nonthermal electrons propagate along the coronal magnetic field thereby exciting dm- and m-wave radio burst emission after acceleration during reconnection or other energy release processes in heights of some Mm to ⩾700 Mm. We summarize the results of some case studies which can be interpreted as radio evidence of magnetic reconnection: under certain conditions, simple spectral structures (pulsation pulses, reverse drift bursts) are formed by simultaneously acting but widely spaced radio sources. Narrowband spikes are emitted as a side-effect during large-scale coronal loop collisions. In dynamic radio spectra, the lower fast mode shock formed in the reconnection outflow appears as type II burst-like but nondrifting emission lane. It has been several times observed at the harmonic mode of the local plasma frequency between 250 and 500 MHz and at heights of ≈200 Mm.