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.     

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.     

Testing the three-dimensional IRI-SIRMUP-P mapping of the ionosphere for disturbed periods

This paper describes the three-dimensional (3-D) electron density mapping of the ionosphere given as output by the assimilative IRI-SIRMUP-P (ISP) model for three different geomagnetic storms. Results of the 3-D model are shown by comparing the electron density profiles given by the model with the ones measured at two testing ionospheric stations: Roquetes (40.8°N, 0.5°E), Spain, and San Vito (40.6°N, 17.8°E), Italy. The reference ionospheric stations from which the autoscaled foF2 and M(3000)F2 data as well as the real-time vertical electron density profiles are assimilated by the ISP model are those of El Arenosillo (37.1°N, 353.3°E), Spain, Rome (41.8°N, 12.5°E), and Gibilmanna (37.9°N, 14.0°E), Italy. Overall, the representation of the ionosphere made by the ISP model is better than the climatological representation made by only the IRI-URSI and the IRI-CCIR models. However, there are few cases for which the assimilation of the autoscaled data from the reference stations causes either a strong underestimation or a strong overestimation of the real conditions of the ionosphere, which is in these cases better represented by only the IRI-URSI model. This ISP misrepresentation is mainly due to the fact that the reference ionospheric stations covering the region mapped by the model turn out to be few, especially for disturbed periods when the ionosphere is very variable both in time and in space and hence a larger number of stations would be required. The inclusion of new additional reference ionospheric stations could surely smooth out this concern.     

Estimation of ballistic coefficients of low altitude debris objects from historical two line elements

This paper presents a new method for estimating ballistic coefficients (BCs) of low perigee debris objects from their historical two line elements (TLEs). The method uses the drag perturbation equation of the semi-major axis of the orbit. For an object with perigee altitude below 700 km, the variation in the mean semi-major axis derived from the TLE is mainly caused by the atmospheric drag effect, and therefore is used as the source in the estimation of the ballistic coefficient. The method is tested using the GRACE satellites, and a number of debris objects with external ballistic coefficient values, and agreements of about 10% are achieved.     

Seasonally steady planetary disturbances in the troposphere and stratosphere as seen in 30 years of NCEP reanalysis data

Zonal velocity and temperature daily global reanalysis data of 30 years are used to search seasonally steady planetary disturbances in the middle troposphere (400 hPa) and middle stratosphere (10 hPa). Significant wavenumber 1, 2 and 3 modes are found. Constant phase lines of zonal velocity 1 modes exhibit significant inclination angles with respect to the meridians. The winter hemisphere generally shows a more significant presence of structures. The Northern Hemisphere (NH) exhibits all over the year a larger amount of structures and more intense amplitudes than the Southern Hemisphere (SH). Middle latitudes exhibit the most significant cases and low latitudes the least significant ones. Longitudinally oriented land–sea transitions at ±$\pm$65° and −35° latitudes appear to play a significant role for the presence of steady planetary modes. The stratosphere exhibits a much simpler picture than the troposphere. Large scale structures with respectively NE–SW (NH) and NW–SE (SH) tilts in the observed temperature and zonal velocity constant phase lines recall the quasi-stationary Rossby wave trains that favor the poleward transport of angular momentum.     

Powered-descent trajectory optimization scheme for Mars landing

This paper presents a trajectory optimization scheme for powered-descent phase of Mars landing with considerations of disturbance. Firstly, θ–D$\theta ''–''D$ method is applied to design a suboptimal control law with descent model in the absence of disturbance. Secondly, disturbance is estimated by disturbance observer, and the disturbance estimation is as feedforward compensation. Then, semi-global stability analysis of the composite controller consisting of the nonlinear suboptimal controller and the disturbance feedforward compensation is proposed. Finally, to verify the effectiveness of proposed control scheme, an application including relevant simulations on a Mars landing mission is demonstrated.     

Design and performance of a vacuum-UV simulator for material testing under space conditions

This paper describes the construction and performance of a VUV-simulator that has been designed to study degradation of materials under space conditions. It is part of the Complex Irradiation Facility at DLR in Bremen, Germany, that has been built for testing of material under irradiation in the complete UV-range as well as under proton and electron irradiation. Presently available UV-sources used for material tests do not allow the irradiation with wavelengths smaller than about 115 nm where common Deuterium lamps show an intensity cut-off. The VUV-simulator generates radiation by excitation of a gas-flow with an electron beam. The intensity of the radiation can be varied by manipulating the gas-flow and/or the electron beam.     

Understanding Substorms in the Magnetotail: Early Development and Recent Progress

This is a concise review of physics of the substorm in the magnetotail. It consists of two parts. The first part summarizes historical developments in the early days of the space age (1960--1975) when the basic concepts such as magnetotail and reconnection were established and the leading model of the substorm was introduced. The second part is an overview of the research conducted in recent years (1995--2010) when very significant advances have been achieved in understanding the substorm physics by virtue of several major satellites missions that addressed the magnetotail physics intensively.