The accuracy of ionogram-derived N(h) profiles

At the Millstone Hill station the Incoherent Scatter Radar (ISR) and a Digisonde 256 are simultaneously operating. Some characteristic true heights determined by both instruments are compared with each other, possible reasons for observed difference are indicated.     

Kalman filter based estimation of differential hardware biases with triangular interpolation technique for IRNSS

Ionosphere delay correction is the main error correction to the computation of single frequency user position using satellite navigation. However ionosphere delay consists of not only delay but also frequency dependent differential hardware biases from satellite and receiver ends. For ionosphere point of view, Indian Regional Navigation Satellite System (IRNSS) service area comes in equatorial anomaly region. It is a unique satellite navigation system which operates at L5 and S frequencies and consists of Geostationary Earth Orbit (GEO) and Geo Synchronous Orbit (GSO) satellite constellation. With IRNSS measurements availability, there is a good opportunity to estimate and analyse differential hardware biases with GEO/GSO combination and with equatorial ionosphere variation. In this paper, Kalman filter based estimation with triangular interpolation technique is used to estimate differential hardware biases for all IRNSS satellites and reference receivers at L5 frequency. The standard deviation of the 15?days of daily estimation of satellite differential hardware biases is in the range of 0.32 to 1.17 TECU for all IRNSS satellites. Similarly, the standard deviation of the 15?days of daily estimation varies up to 2.85 and 6.0 TECU for receiver differential hardware biases during calm and stormy period respectively. The ionosphere delay computed using estimated differential hardware biases is compared with Global Ionosphere Map (GIM) data. A rigorous analysis is carried out to study the error in the estimation in terms of input data noise level, satellite constellation and effect of latitude. Our result reveals that over IRNSS service area, there is an exponential increase in the error in the estimation of receiver differential hardware biases with respect to latitude.     

Nonlinear filtering for sequential spacecraft attitude estimation with real data: Cubature Kalman Filter,Unscented Kalman Filter and Extended Kalman Filter

This article compares the attitude estimated by nonlinear estimator Cubature Kalman Filter with results obtained by the Extended Kalman Filter and Unscented Kalman Filter. Currently these estimators are the subject of great interest in attitude estimation problems, however, mostly the Extended Kalman Filter has been applied to real problems of this nature. In order to evaluate the behavior of the Extended Kalman Filter, Unscented Kalman Filter and Cubature Kalman Filter algorithms when submitted to realistic situations, this paper uses real data of sensors on-board the CBERS-2 remote sensing satellite (China Brazil Earth Resources Satellite). It is observed that, for the case studied in this article, the filters are very competitive and present advantages and disadvantages that should be dealt with according to the requirements of the problem.     

GRACE accelerometer calibration by high precision non-gravitational force modeling

The precise modeling and knowledge of non-gravitational forces acting on satellites is of big interest to many scientific tasks and missions. Since 2002, the twin GRACE satellites have measured these forces in a low Earth orbit with highly precise accelerometers, for about 15?years. Besides the significance for the GRACE mission, these measurement data allow the evaluation of modeling approaches and the improvement of force models. Unfortunately, before any scientific usage, the accelerometer measurements need to be calibrated, namely scale factor and bias have to be regularly estimated.In this study we demonstrate an accelerometer calibration approach, solely based on high precision non-gravitational force modeling without any use of empirically or stochastically estimated parameters, using our in-house developed satellite simulation tool XHPS. The aim of this work is twofold, first we use the accelerometer data and the residuals resulting from the calibration to quantitatively analyze and validate different non-gravitational force model approaches. In a second step, we compare the calibration results to three different calibration methods from different authors, based on gravity field recovery, GPS-based precise orbit determination, and based on modeled accelerations.We consider atmospheric drag forces and winds, as well as radiation forces due to solar radiation pressure, albedo, Earth infrared and thermal radiation (TRP) of the satellite itself. For TRP, we investigate different transient temperature calculation approaches for the satellite surfaces with absorbed power from the aforementioned radiation sources. A detailed finite element model of the satellite is utilized for every force, considering orientation, material properties and shadowing conditions for each element.For cross-track and radial direction, which are mainly affected by the radiative forces, our calibration residuals are quite small when drag is not super dominant (1–3?$\text{nm}/{\text{s}}^2$ for total accelerations around $\pm$50?$\text{nm}/{\text{s}}^2$). For these directions the calibration seems to perform better than the other compared methods, where some bigger differences were found. For the drag dominated along-track direction it is vice versa, here our method is not sensitive enough because the difference between modeled and measured drag is bigger (e.g. residuals around 10?$\text{nm}/{\text{s}}^2$ for total accelerations around $\pm$70?$\text{nm}/{\text{s}}^2$ for low solar activity). In along-track direction the orbit determination based methods are more sensitive and produce more reliable results. Results for the complete GRACE mission time span from 2003 to 2017 are shown, covering different seasonal environmental conditions.     

On the detection of anomalous seismo-ionospheric behavior in the presence of space weather stimuli for large earthquakes

Anomalous behavior of ionospheric total electron content (TEC) prior to earthquake has been observed in many studies. Evidence of such seismo-ionospheric coupling effects suggests that it is plausible to rely on TEC signatures for early earthquake warning. However, the detection of pre-earthquake TEC anomalies (PETA) has not been adopted in practice due to two pertinent issues. Firstly, the effects of space weather activity can affect TEC levels and cause anomalous behavior in the TEC. Usually arbitrary thresholds are set for space weather indices to eliminate TEC anomaly due to space weather effects. Secondly, the choice regarding moving time-window length used to characterise background variation of TEC within the statistical envelope approach has an effect on detection of PETA. While the rule-of-thumb in selecting the moving window length is to have a time window capable of capturing background variability and short-term fluctuations, the length of the time window used in the literature varies with little justification. In this study, a critical examination is conducted on the statistical envelope approach and in particular, to eliminate the effect of space weather activity without the use of arbitrary space indices to detect PETA. A two-part PETA identification procedure is proposed, consisting of wavelet analyses isolating non-earthquake TEC contributions, followed by the statistical envelope method using a moving window length standardized based on observed periodicities and statistical implications is suggested. The approach is tested on a database of 30 large earthquakes (M?≥?7.0). The proposed procedure shows that PETA can be detected prior to earthquakes at higher confidence levels without the need to separately check for space weather activity. More importantly, the procedure was able to detect PETA for studies where it was previously reported that PETA could not be detected or detected convincingly.     

Statistical properties of the IMF clock angle in the solar wind with northward and southward interplanetary magnetic field based on ACE observation from 1998 to 2009: Dependence on the temporal scale of the solar wind

Utilizing ACE satellite observations from 1998 to 2009, we performed the elaborate study on the properties of the clock angle θ_CA (arctan(By/Bz) (?90° to 90°) of the interplanetary magnetic field (IMF) in the solar wind at 1?AU. The solar wind with northward IMF (NW-IMF) and southward IMF (SW-IMF) are analyzed, independently. Statistical analysis shows that the solar wind with SW-IMF and NW-IMF has similar properties in general, including their durations, the IMF Bz and By components, and the IMF θ_CA. Then, the solar wind with NW-IMF (SW-IMF) is classified into five different temporal scales according to the duration of the NW-IMF (SW-IMF), i.e., very-short wind of 10–30?min, short-scale wind of 0.5–1?h, moderate-scale wind of 1–3?h, long-scale wind of 3–5?h, and super-long wind >5?h. Our analysis reveals that the IMF θ_CA has a distinct decrease with increase of the temporal scale of the solar wind. Next, the solar wind is classified into two groups, i.e., the high-speed solar wind (>450?km/s) and the low-speed solar wind (<450?km/s). Our analysis indicates that the IMF θ_CA depends highly on the solar wind speed. Statistically, high-speed solar wind tends to have larger IMF θ_CA than low-speed solar wind. The evolutions of the solar wind and IMF with the solar activity are further studied, revealing no clear solar variation of the IMF θ_CA. Finally, we analyze the monthly variation of the IMF θ_CA. Superposed epoch result strongly suggests the seasonal variation of the IMF θ_CA.     

Effects observed in the equatorial and low latitude ionospheric F-region in the Brazilian sector during low solar activity geomagnetic storms and comparison with the COSMIC measurements

The main objective of the present investigation has been to compare the ionospheric parameters (NmF2 and hmF2) observed by two ground-based ionospheric sounders (one at PALMAS- located near the magnetic equator and the other at Sao Jose dos Campos-located in the low-latitude region) in the Brazilian sector with that by the satellite FORMOSAT-3/COSMIC radio occultation (RO) measurements during two geomagnetic storms which occurred in December 2006 and July 2009. It should be pointed out that in spite of increasing the latitude (to 10°) and longitude (to 20°) around the stations; we had very few common observations. It has been observed that both the peak electron density (NmF2) and peak height (hmF2) observed by two different techniques (space-borne COSMIC and ground-based ionosondes) during both the geomagnetic storm events compares fairly well (with high correlation coefficients) at the two stations in the Brazilian sector. It should be pointed out that due to equatorial spread F (ESF) in the first storm (December 2006) and no-reflections from the ionosphere during nighttime in the second storm (July 2009), we had virtually daytime data from the two ionosondes.     

Magnetically suspended MEMS spinning wheel gyro

A concept of a magnetically suspended spinning wheel gyro is described. The gyro is based on a patented, planar magnetic actuator and position sensor configuration with a gyro wheel that is rotated by a multi-phase electromagnetic spin motor. The motor and actuator/sensor configurations, in addition to providing the necessary forces to suspend the spinning wheel, are amenable to fabrication using some of the developing MEMS fabrication technologies, making batch fabrication of the gyro possible. This provide the capability for high production yield in high volume through the use of batch processing, without the need for skilled hand labor, resulting in high yields, high reliability, and low cost of ownership. The gyro concept is described in detail. The high-speed rotation of the wheel, somewhat smaller than a dime, can produce an angular momentum much larger than that of the Coriolis force based MEMS gyros now receiving intensive development, and is expected to provide considerably higher performance than currently available from this class of instrument. Two axes of angular rate information and three axes of acceleration information are provided by the control loops that center the wheel within the case. Thus, two such magnetically suspended spinning wheel gyros can be used to implement a complete IMU, and can also provide redundancy in three of the five degrees of freedom     

Estimation and classification of FM signals using time frequency transforms

This work deals with electronic warfare: we have to process an intercepted radar-emitted signal assumed to be frequency-modulated (FM). A battery of parametric models being available, we face the two main problems: the parameter estimation, given a model, and the classification, given the battery. For the former, we propose the classical maximum likelihood estimator (MLE) and an original numerical scheme to reach it. For the latter, we construct a test based on the estimated variances cross-ratios. Both perform very efficiently with respect to (w.r.t.) the Cramer-Rao lower bound (CRLB) and the rate of correct classification, respectively, for low signal-to-noise ratios (SNRs).     

Analysis of CGSDS file delivery protocol: immediate NAK mode

We analyze the immediate negative acknowledgment (NAK) mode of the consultative committee for space data systems (CCSDS) file delivery protocol for the single-hop file transfer operation. We propose a timer setting rule that minimizes the expected time taken to transfer a file under the constraint that the throughput efficiency is maximized. Then, we derive the expected file delivery time and compare it with that of the deferred NAK mode. The main contribution is the closed-form expression for evaluating the performance metrics.