Digital signal processing and numerical analysis for radar in geophysical applications

Numerical solutions for signal processing are described in this work as a contribution to study of echo detection methods for ionospheric sounder design. The ionospheric sounder is a high frequency radar for geophysical applications. The main detection approach has been done by implementing the spread-spectrum techniques using coding methods to improve the radar’s range resolution by transmitting low power. Digital signal processing has been performed and the numerical methods were checked. An algorithm was proposed and its computational complexity was calculated.     

Systematic examination of the geomagnetic storm sudden commencement using multi resolution analysis

Geomagnetic storm sudden commencements (SSC) contain a wealth of information which is useful in many applications. It is important to point out that the SSCs used in this study are sampled at the rate of one sample per second in order get use of such high resolution data. In this paper, two studies are made on the geomagnetic SSC and an SSC onset automatic detection algorithm is introduced. The first study is about finding the relationship between the SSC rise time and its amplitude. Where it is found that there is a positive correlation between the amplitude and the amplitude gradient which is the amplitude divided by rise time. The second study is checking the spectrum of the SSC, starting from its onset until the end of the SSC rise time. This check had proved that the SSC contains both low and high frequency regions. This led us to introduce a new term, namely the SSC variation rate (VR). This VR is defined as the maximum rate of change of the field in the higher-frequency region of the SSC. These two studies were the guide to build an SSC automatic detector of one sample per second data using multi resolution analysis (MRA) of the discrete wavelet transform (DWT). The data set contains 134 SSCs with different VRs that were collected from the Circum-pan Pacific Magnetometer Network (CPMN). It is found that the standard deviation of the detection error is 41 s and that the average error is 9 s. From the calculated error distribution function, it is found that the detection error is within the range of −1.5 to 3 min. The detection process, as will be shown in the article, takes 70 s for one station and 3 min if the decision is related to the detection(s) of other stations. These results demonstrate the superiority of the proposed algorithm over other algorithms in which the detection error ranges between −8 and 5 min and the detection process takes 2–10 min. In addition to being faster and more accurate than the other algorithms, the proposed algorithm is the first algorithm that automatically detects the SSC onset times from high-resolution data unlike previous studies that focused on determining the SSC times automatically using one-minute resolution data.     

Spatial heterogeneity in the radiogenic activity of the lunar interior: Inferences from CHACE and LLRI on Chandrayaan-1

In the past, clues on the potential radiogenic activity of the lunar interior have been obtained from the isotopic composition of noble gases like Argon. Excess Argon (40) relative to Argon (36), as compared to the solar wind composition, is generally ascribed to the radiogenic activity of the lunar interior. Almost all the previous estimates were based on, ‘on-the-spot’ measurements from the landing sites. Relative concentration of the isotopes of ⁴⁰Ar and ³⁶Ar along a meridian by the Chandra’s Altitudinal Composition Explorer (CHACE) experiment, on the Moon Impact Probe (MIP) of India’s first mission to Moon, has independently yielded clues on the possible spatial heterogeneity in the radiogenic activity of the lunar interior in addition to providing indicative ‘antiquity’ of the lunar surface along the ground track over the near side of the moon. These results are shown to broadly corroborate the independent topography measurements by the Lunar Laser Ranging Instrument (LLRI) in the main orbiter Chandrayaan-1. The unique combination of these experiments provided high spatial resolution data while indicating the possible close linkages between the lunar interior and the lunar ambience.     

Concept of wireless sensor network for future in-situ exploration of lunar ice using wireless impedance sensor

It is known that a wireless sensor network uses some sort of sensors to detect a physical quantity of interest, in general. The wireless sensor network is a potential tool for exploring the difficult-to-access area on the earth and the concept may be extended to space applications in future. Recently, lunar water has been detected by a few lunar missions using remote sensing techniques. The lunar water is expected to be in the form of ice at very low temperatures of permanently dark regions on the moon. To support the remote observations and also to find out potential ice bearing sites on the moon, in-situ measurement of the lunar ice is essential. However, a rover may not be able to reach the permanently shadowed regions due to terrain irregularity. One possibility to access such areas is to use a wireless sensor network on the lunar surface.     

VHF radio response of the near Earth space during solar activity growth in October, 2003

The analysis of observations of very high frequency radio noise intensity at the middle latitude on a frequency f = 500 MHz from 14th till 26th of October, 2003 is presented. These data are compared with the solar radio bursts in the range of frequencies 1–14 MHz registered by RAD2 receiver of the WAVES device installed on board the WIND spacecraft.     

Modelling and observing Jovian electron propagation times in the inner heliosphere

The propagation of Jovian electrons in interplanetary space was modelled by solving the relevant transport equation numerically through the use of stochastic differential equations. This approach allows us to calculate, for the first time, the propagation time of Jovian electrons from the Jovian magnetosphere to Earth. Using observed quiet-time increases of electron intensities at Earth, we also derive values for this quantity. Comparing the modelled and observed propagation times we can gauge the magnitude of the transport parameters sufficiently to place a limit on the 6 MeV Jovian electron flux reaching Earth. We also investigate how the modelled propagation time, and corresponding Jovian electron flux, varies with the well-known ∼13 month periodicity in the magnetic connectivity of Earth and Jupiter. The results show that the Jovian electron intensity varies by a factor of ∼10 during this cycle of magnetic connectivity.     

Symmetry and asymmetry of ionospheric weather at magnetic conjugate points for two midlatitude observatories

Variations of the ionospheric weather W-index for two midlatitude observatories, namely, Grahamstown and Hermanus, and their conjugate counterpart locations in Africa are studied for a period from October 2010 to December 2011. The observatories are located in the longitude sector, which has consistent magnetic equator and geographic equator so that geomagnetic latitudes of the line of force are very close to the corresponding geographic latitudes providing opportunity to ignore the impact of the difference of the gravitational field and the geomagnetic field at the conjugate points on the ionosphere structure and dynamics. The ionosondes of Grahamstown and Hermanus provide data of the critical frequency (foF2), and Global Ionospheric Maps (GIM) provide the total electron content (TECgps) along the magnetic field line up to the conjugate point in the opposite hemisphere. The global model of the ionosphere, International Reference Ionosphere, extended to the plasmasphere altitude of 20,200 km (IRI-Plas) is used to deliver the F2 layer peak parameters from TECgps at the magnetic conjugate area. The evidence is obtained that the electron gas heated by day and cooled by night at the summer hemisphere as compared with the opposite features in the conjugate winter hemisphere testifies on a reversal of plasma fluxes along the magnetic field line by the solar terminator. The ionospheric weather W-index is derived from NmF2 (related with foF2) and TECgps data. It is found that symmetry of W-index behavior in the magnetic conjugate hemispheres is dominant for the equinoxes when plasma movement along the magnetic line of force is imposed on symmetrical background electron density and electron content. Asymmetry of the ionospheric storm effects is observed for solstices when the plasma diffuse down more slowly into the colder winter hemisphere than into the warmer summer hemisphere inducing either plasma increase (positive phase) or decrease (negative phase of W-index) in the ionospheric and plasmaspheric plasma density.     

Total electron content of the ionosphere at two stations in East Africa during the 24–25 October 2011 geomagnetic storm

The equatorial ionosphere has been known to become highly disturbed and thus rendering space-based navigation unreliable during space weather events, such as geomagnetic storms. Modern navigation systems, such as the Global Positioning System (GPS) use radio-wave signals that reflect from or propagate through the ionosphere as a means of determining range or distance. Such systems are vulnerable to effects caused by geomagnetic storms, and their performance can be severely degraded. This paper analyses total electron content (TEC) and the corresponding GPS scintillations using two GPS SCINDA receivers located at Makerere University, Uganda (Lat: 0.3^o N; Lon: 32.5^o E) and at the University of Nairobi, Kenya (Lat: 1.3^o S; Lon: 36.8^o E), both in East Africa. The analysis shows that the scintillations actually correspond to plasma bubbles. The occurrence of plasma bubbles at one station was correlated with those at the other station by using observations from the same satellite. It was noted that some bubbles develop at one station and presumably “die off” before reaching the other station. The paper also discusses the effects of the geomagnetic storm of the 24–25 October 2011 on the ionospheric TEC at the two East African stations. Reductions in the diurnal TEC at the two stations during the period of the storm were observed and the TEC depletions observed during that period showed much deeper depletions than on the non-storm days. The effects during the storm have been attributed to the uplift of the ionospheric plasma, which was then transported away from this region by diffusion along magnetic field lines.     

Quality variation of GPS satellite clocks on-orbit using IGS clock products

Over 60% clocks on board of the GPS satellites are working longer than their designed life. Therefore realizing their stabilities in a long time scales is essential to GPS navigation and positioning plus IGS time scale maintaining. IGS clock products from 2001 to 2010 are used to analyze the GPS satellite clock qualities such as frequency stabilities and clock noise level. We find out that for the clocks of Block IIA satellites the frequency stabilities and clock noise are 10 times worse than that of the Block IIR and IIR-M satellites. Moreover, the linear relationships between frequency stabilities and clock residuals have been deduced with an accuracy of better than 0.02 ns. Specially, it is noticed that the clock of the PRN27 is instable and the relationship between the frequency stability and residuals is at least a quadratic curve. Therefore, we suggested that GPS satellite clocks should be weighted by their quality levels in application, and the observations of the Block IIA should not be used for real-time positioning which required precision better than one meter.     

Numerical modeling of propagation of breaking nonlinear acoustic-gravity waves from the lower to the upper atmosphere

Acoustic-gravity waves (AGWs) observed in the upper atmosphere may be generated near the Earth’s surface due to a variety of meteorological sources. Two-dimensional simulations of vertical propagation and breaking of nonlinear AGWs in the atmosphere are performed. Forcing near the Earth’s surface is used as the AGW source in the model. We use a numerical method based on finite-difference analogues of fundamental conservation laws for solving atmospheric hydrodynamic equations. This approach selects physically correct generalized solutions of the wave hydrodynamic equations. Numerical simulations are performed in a representative region of the Earth’s atmosphere up to altitude 500 km. Vertical profiles of temperature, density, molecular viscosity and heat conductivity were taken from the standard atmosphere model MSIS-90 for January. Calculations were made for different amplitudes and frequencies of lower boundary wave forcing. It is shown that after activating the tropospheric wave forcing, the initial pulse of AGWs may very quickly propagate to altitudes of 100 km and above and relatively slowly dissipate due to molecular viscosity and heat conduction. This may increase the role of transient nonstationary waves in effective energy transport and variations of atmospheric parameters and gas admixtures in a broad altitude range.