Comparison of two measurement fusion methods forKalman-filter-based multisensor data fusion

Currently there exist two commonly used measurement fusion methods for Kalman-filter-based multisensor data fusion. The first (Method I) simply merges the multisensor data through the observation vector of the Kalman filter, whereas the second (Method II) combines the multisensor data based on a minimum-mean-square-error criterion. This paper, based on an analysis of the fused state estimate covariances of the two measurement fusion methods, shows that the two measurement fusion methods are functionally equivalent if the sensors used for data fusion, with different and independent noise characteristics, have identical measurement matrices. Also presented are simulation results on state estimation using the two measurement fusion methods, followed by the analysis of the computational advantages of each method     

Huygens Probe Aerosol Collector Pyrolyser Experiment

ACP's main objective is the chemical analysis of the aerosols in Titan's atmosphere. For this purpose, it will sample the aerosols during descent and prepare the collected matter (by evaporation, pyrolysis and gas products transfer) for analysis by the Huygens Gas Chromatograph Mass Spectrometer (GCMS). A sampling system is required for sampling the aerosols in the 135'32 km and 22'17 km altitude regions of Titan's atmosphere. A pump unit is used to force the gas flow through a filter. In its sampling position, the filter front face extends a few mm beyond the inlet tube. The oven is a pyrolysis furnace where a heating element can heat the filter and hence the sampled aerosols to 250 °C or 600 °C. The oven contains the filter, which has a thimble-like shape (height 28 mm). For transferring effluent gas and pyrolysis products to GCMS, the carrier gas is a labeled nitrogen ¹⁵N₂, to avoid unwanted secondary reactions with Titan's atmospheric nitrogen. Aeraulic tests under cold temperature conditions were conducted by using a cold gas test system developed by ONERA. The objective of the test was to demonstrate the functional ability of the instrument during the descent of the probe and to understand its thermal behavior, that is to test the performance of all its components, pump unit and mechanisms. In order to validate ACP's scientific performance, pyrolysis tests were conducted at LISA on solid phase material synthesized from experimental simulation. The chromatogram obtained by GCMS analysis shows many organic compounds. Some GC peaks appear clearly from the total mass spectra, with specific ions well identified thanks to the very high sensitivity of the mass spectrometer. The program selected for calibrating the flight model is directly linked to the GCMS calibration plan. In order not to pollute the two flight models with products of solid samples such as tholins, we excluded any direct pyrolysis tests through the ACP oven during the first phase of the calibration. Post probe descent simulation of flight results are planned, using the much representative GCMS and ACP spare models. This revised version was published online in August 2006 with corrections to the Cover Date.     

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.     

Vertical motions in Thailand after the 2004 Sumatra–Andaman Earthquake from GPS observations and its geophysical modelling

Following previous findings from ongoing GPS research in Thailand since 2004 we continue to exploit the GPS technique to monitor and model land motions induced by the Sumatra–Andaman Earthquake. Our latest results show that up to the end of 2010, Thailand has been co-seismically displaced and is subsequently undergoing a post-seismic horizontal deformation with total displacements (co-seismic plus post-seismic) ranging from 10.5 to 74.7 cm. We observed the largest horizontal displacements in the southern part of Thailand and moderate and small displacements in the central and northern parts. In addition to horizontal displacements throughout Thailand, continuous GPS measurements show that large parts of Thailand are subsiding at rates up to 1 cm/yr. It is the first time that such vertical post-seismic deformations at large distances (650–1500 km away from the Earthquake’s epicentre) have been recorded. We have investigated the physical processes leading to the observed subsidence. While after-slip on the subduction interface induces negligible or even slightly positive vertical motions, relaxation in the asthenosphere is associated with a sizable subsidence. Predictions from a 3D finite element model feature an asthenosphere with an effective viscosity of the order of 3 * 10¹⁸ Pas, fit the horizontal post-seismic data and the observed subsidence well. This model is then used to predict the subsidence over the whole seismic cycle. The subsidence should go on with a diminishing rate through the next two decades and its final magnitude should not exceed 10 cm in the Bangkok area.     

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.     

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.