Longitudinal variation of the equatorial ionosphere: Modeling and experimental results

We describe a new version of the Parameterized Regional Ionospheric Model (PARIM) which has been modified to include the longitudinal dependences. This model has been reconstructed using multidimensional Fourier series. To validate PARIM results, the South America maps of critical frequencies for the E (foE) and F (foF2) regions were compared with the values calculated by Sheffield Plasmasphere-Ionosphere Model (SUPIM) and IRI representations. PARIM presents very good results, the general characteristics of both regions, mainly the presence of the equatorial ionization anomaly, were well reproduced for equinoctial conditions of solar minimum and maximum. The values of foF2 and hmF2 recorded over Jicamarca (12°S; 77°W; dip lat. 1°N; mag. declination 0.3°) and sites of the conjugate point equatorial experiment (COPEX) campaign Boa Vista (2.8°N; 60.7°W; dip lat. 11.4°; mag. declination −13.1°), Cachimbo (9.5°S; 54.8°W; dip lat. −1.8°; mag. declination −15.5°), and Campo Grande (20.4°S; 54.6°W; dip lat. −11.1°; mag. declination −14.0°) have been used in this work. foF2 calculated by PARIM show good agreement with the observations, except during morning over Boa Vista and midnight-morning over Campo Grande. Some discrepancies were also found for the F-region peak height (hmF2) near the geomagnetic equator during times of F₃ layer occurrences. IRI has underestimated both foF2 and hmF2 over equatorial and low latitude sectors during evening-nighttimes, except for Jicamarca where foF2 values were overestimated.     

Universal scaling laws for fully-developed magnetic field turbulence near and far upstream of the Earth’s bow shock

We analyze the multifractal scaling of the modulus of the interplanetary magnetic field near and far upstream of the Earth’s bow shock, measured by Cluster and ACE, respectively, from 1 to 3 February 2002. The maximum order of the structure function is carefully estimated for each time series using two different techniques, to ensure the validity of our high-order statistics. The first technique consists of plotting the integrand of the pth order structure function, and the second technique is a quantitative method which relies on the power-law scaling of the extreme events. We compare the scaling exponents computed from the structure functions of magnetic field differences with the predictions obtained by the She–Lévêque model of intermittency in anisotropic magnetohydrodynamic turbulence. Our results show a good agreement between the model and the observations near and far upstream of the Earth’s bow shock, rendering support for the modelling of universal scaling laws based on the Kolmogorov phenomenology in the presence of sheet-like dissipative structures.     

The effect of geocenter motion on Jason-2 orbits and the mean sea level

We compute a series of Jason-2 GPS and SLR/DORIS-based orbits using ITRF2005 and the std0905 standards ( Lemoine et al., 2010). Our GPS and SLR/DORIS orbit data sets span a period of 2 years from cycle 3 (July 2008) to cycle 74 (July 2010). We extract the Jason-2 orbit frame translational parameters per cycle by the means of a Helmert transformation between a set of reference orbits and a set of test orbits. We compare the annual terms of these time-series to the annual terms of two different geocenter motion models where biases and trends have been removed. Subsequently, we include the annual terms of the modeled geocenter motion as a degree-1 loading displacement correction to the GPS and SLR/DORIS tracking network of the POD process. Although the annual geocenter motion correction would reflect a stationary signal in time, under ideal conditions, the whole geocenter motion is a non-stationary process that includes secular trends. Our results suggest that our GSFC Jason-2 GPS-based orbits are closely tied to the center of mass (CM) of the Earth consistent with our current force modeling, whereas GSFC’s SLR/DORIS-based orbits are tied to the origin of ITRF2005, which is the center of figure (CF) for sub-secular scales. We quantify the GPS and SLR/DORIS orbit centering and how this impacts the orbit radial error over the globe, which is assimilated into mean sea level (MSL) error, from the omission of the annual term of the geocenter correction. We find that for the SLR/DORIS std0905 orbits, currently used by the oceanographic community, only the negligence of the annual term of the geocenter motion correction results in a – 4.67 ± 3.40 mm error in the Z-component of the orbit frame which creates 1.06 ± 2.66 mm of systematic error in the MSL estimates, mainly due to the uneven distribution of the oceans between the North and South hemisphere.     

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.     

Performance evaluation of IRI-2007 at equatorial latitudes and its Matlab version for GNSS applications

International Reference Ionosphere (IRI) model is the widely used empirical model for ionospheric predictions, especially TEC which is an important parameter for radio navigation and communication. The Fortran based IRI-2007 does not support real-time interactive visualization and debugging. Therefore, the source code is converted into Matlab and is validated for the purposes of this study. This facilitates easy representation of results and for near real-time implementation of IRI in the applications including spacecraft launching, now casting, pseudolite based navigation systems etc. In addition, the vertical delay results over the equatorial region derived from IRI and GPS data of three IGS stations namely Libreville (Garbon, Africa), Brasilia (Brazil, South America) and Hyderabad (India, Asia) are compared. As the IRI model does not account for plasmasphere TEC, the vertical delays are underestimated compared to vertical delays of GPS signals. Therefore, the model should be modified accordingly for precise TEC estimation.     

Active debris removal of multiple priority targets

Today’s space debris environment shows major concentrations of objects within distinct orbital regions for nearly all size regimes. The most critical region is found at orbital altitudes near 800 km with high declinations. Within this region many satellites are operated in so called sun-synchronous orbits (SSO). Among those, there are Earth observation, communication and weather satellites. Due to the orbital geometry in SSO, head-on encounters with relative velocities of about 15 km/s are most probable and would thus result in highly energetic collisions, which are often referred to as catastrophic collisions, leading to the complete fragmentation of the participating objects. So called feedback collisions can then be triggered by the newly generated fragments, thus leading to a further population increase in the affected orbital region. This effect is known as the Kessler syndrome.     

Height of shock formation in the solar corona inferred from observations of type II radio bursts and coronal mass ejections

Employing coronagraphic and EUV observations close to the solar surface made by the Solar Terrestrial Relations Observatory (STEREO) mission, we determined the heliocentric distance of coronal mass ejections (CMEs) at the starting time of associated metric type II bursts. We used the wave diameter and leading edge methods and measured the CME heights for a set of 32 metric type II bursts from solar cycle 24. We minimized the projection effects by making the measurements from a view that is roughly orthogonal to the direction of the ejection. We also chose image frames close to the onset times of the type II bursts, so no extrapolation was necessary. We found that the CMEs were located in the heliocentric distance range from 1.20 to 1.93 solar radii (Rs), with mean and median values of 1.43 and 1.38 Rs, respectively. We conclusively find that the shock formation can occur at heights substantially below 1.5 Rs. In a few cases, the CME height at type II onset was close to 2 Rs. In these cases, the starting frequency of the type II bursts was very low, in the range 25–40 MHz, which confirms that the shock can also form at larger heights. The starting frequencies of metric type II bursts have a weak correlation with the measured CME/shock heights and are consistent with the rapid decline of density with height in the inner corona.     

Subsidence rate monitoring of Aghajari oil field based on Differential SAR Interferometry

Land subsidence, due to natural or anthropogenic processes, causes significant costs in both economic and structural aspects. That part of subsidence observed most is the result of human activities, which relates to underground exploitation. Since the gradual surface deformation is a consequence of hydrocarbon reservoirs extraction, the process of displacement monitoring is amongst the petroleum industry priorities. Nowadays, Differential SAR Interferometry, in which satellite images are utilized for elevation change detection and analysis – in a millimetre scale, has proved to be a more real-time and cost-effective technology in contrast to the traditional surveying method. In this study, surface displacements in Aghajari oil field, i.e. one of the most industrious Iranian hydrocarbon sites, are being examined using radar observations. As in a number of interferograms, the production wells inspection reveals that surface deformation signals develop likely due to extraction in a period of several months. In other words, different subsidence or uplift rates and deformation styles occur locally depending on the geological conditions and excavation rates in place.