Modeling errors in upward continuation for INS gravity compensation

Accurate upward continuation of gravity anomalies supports future precision, free-inertial navigation systems, since the latter cannot by themselves sense the gravitational field and thus require appropriate gravity compensation. This compensation is in the form of horizontal gravity components. An analysis of the model errors in upward continuation using derivatives of the standard Pizzetti integral solution (spherical approximation) shows that discretization of the data and truncation of the integral are the major sources of error in the predicted horizontal components of the gravity disturbance. The irregular shape of the data boundary, even the relatively rough topography of a simulated mountainous region, has only secondary effect, except when the data resolution is very high (small discretization error). Other errors due to spherical approximation are even less important. The analysis excluded all measurement errors in the gravity anomaly data in order to quantify just the model errors. Based on a consistent gravity field/topographic surface simulation, upward continuation errors in the derivatives of the Pizzetti integral to mean altitudes of about 3,000 and 1,500 m above the mean surface ranged from less than 1 mGal (standard deviation) to less than 2 mGal (standard deviation), respectively, in the case of 2 arcmin data resolution. Least-squares collocation performs better than this, but may require significantly greater computational resources.     

Small-world characteristics on transportation networks: a perspective from network autocorrelation

Studies on small-world networks have received intensive interdisciplinary attention during the past several years. It is well-known among researchers that a small-world network is often characterized by high connectivity and clustering, but so far there exist few effective approaches to evaluate small-world properties, especially for spatial networks. This paper proposes a method to examine the small-world properties of spatial networks from the perspective of network autocorrelation. Two network autocorrelation statistics, Moran’s I and Getis–Ord’s G, are used to monitor the structural properties of networks in a process of “rewiring” networks from a regular to a random network. We discovered that Moran’s I and Getis–Ord’s G tend to converge and have relatively low values when properties of small-world networks emerge. Three transportation networks at the national, metropolitan, and intra-city levels are analyzed using this approach. It is found that spatial networks at these three scales possess small-world properties when the correlation lag distances reach certain thresholds, implying that the manifestation of small-world phenomena result from the interplay between the network structure and the dynamics taking place on the network.      

Correcting superconducting gravity time-series using rainfall modelling at the Vienna and Membach stations and application to Earth tide analysis

We demonstrate the possibility to improve the signal-to-noise ratio of superconducting gravity time-series by correcting for local hydrological effects. Short-term atmospheric events associated with heavy rain induce step-like gravity signals that deteriorate the frequency spectrum estimates. Based on 4D modeling constrained by high temporal resolution rain gauge data, rainfall admittances for the Vienna and Membach superconducting gravity stations are calculated. This allows routine correction for Newtonian rain water effects, which reduces the standard deviation of residuals after tidal parameter adjustment by 10%. It also improves the correction of steps of instrumental origin when they coincide with step-like water mass signals.     

Universal time from VLBI single-baseline observations during CONT08

The IVS Intensive sessions are single-baseline, 1-h VLBI sessions carried out everyday in order to determine Universal Time (UT1). We investigate different possibilities to improve the results of such sessions. We do this investigation by extracting 2-h single-baseline sessions from the CONT08 data set. These are analysed like normal Intensives, and the results are compared to the results of the analysis of the full CONT08 data set. We find that tropospheric asymmetry is the major error source for the single-baseline sessions. It is possible to improve the accuracy of the estimated UT1 either by using accurate a priori tropospheric gradients or by estimating gradients in the data analysis.     

Hyperspectral Radiometry to Quantify the Grades of Iron Ores of Noamundi and Joda Mines,Eastern India

This paper present the results of a preliminary study to assess the potential of the visible, NIR and SWIR energy of the EMR in differentiating iron ores of different grades in a rapid manner using hyperspectral radiometry. Using different iron ore samples from Noamundi and Joda mines, Jharkhand and Orissa, states of India, certain spectro-radiometric measurements and geochemical analysis were carried out and the results have been presented. It was observed that the primary spectral characteristics of these iron ores lie in the 850 to 900 nm and 650–750 nm regions. The spectral parameters for each curve used for studying the iron ores are: (i) the slopes of the spectral curve in 685–725 nm region; (ii) position of the peak with respect to wavelength in 730–750 nm region and (iii) radius of curvature of the absorption trough in the 850–900 nm region. Comparison of these spectral parameters and the geochemistry of the samples indicates that the position of the peak of the curve in 730–750 nm region shifts towards longer wavelength with increasing iron oxide content, while the slope of the curvature in the 685–725 nm region has a strong negative correlation with the iron oxide content of the samples. Similarly, a strong negative correlation is observed between the radius of curvature of the 850–900 nm absorption trough and the iron oxide content. Such strong correlations indicate that hyperspectral radiometry in the visible and NIR regions can give a better estimate and quantification of the grades of iron ores. This study has demonstrated that generation of empirical models using hyperspectral radiometric techniques is helpful to quantify the grade of iron ores with limited geochemical analysis.     

Comparison of Two Data Smoothing Techniques for Vegetation Spectra Derived From EO-1 Hyperion

Hyperspectral data are generally noisier compared to broadband multispectral data because their narrow bandwidth can only capture very little energy that may be overcome by the self-generated noise inside the sensors. It is desirable to smoothen the reflectance spectra. This study was carried out to see the effect of smoothing algorithms - Fast-Fourier Transform (FFT) and Savitzky–Golay (SG) methods on the statistical properties of the vegetation spectra at varying filter sizes. The data used in the study is the reflectance spectra data obtained from Hyperion sensor over an agriculturally dominated area in Modipuram (Uttar Pradesh). The reflectance spectra were extracted for wheat crop at different growth stages. Filter sizes were varied between 3 and 15 with the increment of 2. Paired t-test was carried out between the original and the smoothed data for all the filter sizes in order to see the extent of distortion with changing filter sizes. The study showed that in FFT, beyond filter size 11, the number of locations within the spectra where the smooth spectra were statistically different from its original counterpart increased to 14 and reaches 21 at the filter size 15. However, in SG method, number of statistically different locations were more than those found in the FFT, but the number of locations did not changing drastically. The number of statistically disturbed locations in SG method varied between 16 and 19. The optimum filter size for smoothing the vegetation spectra was found to be 11 in FFT and 9 in SG method.     

Precise ionosphere-free single-frequency GNSS positioning

Ionospheric delays can be efficiently eliminated from single-frequency data using a combination of carrier phases and code ranges. Unfortunately, GPS and GLONASS ranges are relatively noisy which can limit the use of the positioning method. Nevertheless, position standard deviations are in the range of 6–8 cm (horizontal) and 7–9 cm (3d) obtained from diurnal data batches from selected IGS reference stations can be further reduced to 2–3 cm (3d) for weekly smoothed averages. GPS data sets collected in Ghana (Africa) reveal a typical level of 10 cm of deviation that must be anticipated under average conditions. Looking at the future of GNSS, the European Galileo system will, in contrast to GPS, provide the broadband signal E5 that is by far less affected by multipath thus providing rather precise range measurements. Simulated processing runs featuring both high ionospheric and tropospheric delay variations show a 3d position precision of 4 cm even for a data batch as short as just 1 h, whereas GPS L1/Galileo E1 performance is close to 13 cm for the same data set.     

Attitude estimation based on fusion of gyroscopes and single antenna GPS for small UAVs under the influence of vibration

An attitude estimation method is presented for small unmanned aerial vehicles (UAVs) powered by a piston engine which is the major source of vibration. Vibration of the engine significantly degrades the accuracy of the inertial measurement unit, especially for low-cost sensors that are based on micro electro-mechanical system. Therefore, a vibration model for a small UAV is proposed in order to examine the influence of vibration on attitude estimation with different sensors. The model is derived based on spectrum analysis with short-time Fourier transform. The vibration is compared with the drift of the gyroscope in order to examine the impact on attitude estimation. An attitude estimation method that fuses the gyroscopes and single antenna global positioning system (GPS) is proposed to mitigate the influence of engine vibration and gyroscope drift. The quaternion-based extended Kalman filter is implemented to fuse the sensors. This filter fuses the angular rates measured by the gyroscopes and the pseudo-attitude derived from the GPS velocity to estimate the attitude of the UAV. Simulations and experiment results indicate that the proposed method performs well both in short-term and long-term accuracy even though the gyroscopes are affected by drift and vibration noise, while the pseudo-attitude contains severe noise.     

GPS tomography in the polar cap: comparison with ionosondes and in situ spacecraft data

Tomographic 4D reconstructions of ionospheric anomalies appearing in the high-latitude polar cap region are compared with plasma density measurements by digital ionosonde located near the north magnetic pole at Eureka station and with in situ plasma measurements on-board DMSP spacecraft. The moderate magnetic storm of 14–17 October 2002 is taken as an example of a geomagnetic disturbance which generates large-scale ionospheric plasma anomalies at mid-latitudes and in the polar cap region. Comparison of the GPS tomographic reconstructions over Eureka station with the ionosonde measurements of the F layer peak densities indicates that the GPS tomography correctly predicts the time of arrival and passage of the ionospheric tongue of ionization over the magnetic pole area, although the tomographic technique appears to under-estimate the value of F peak plasma density. Comparison with the in situ plasma measurements by the DMSP SSIES instruments shows that the GPS tomography correctly reproduces the large-scale spatial structure of ionospheric anomalies over a wide range of latitudes from mid-latitudes to the high-latitude polar cap region, though the tomographic reconstructions tend to over-estimate the density of the topside ionosphere at 840 km DMSP orbit. This study is essential for understanding the quality and limitations of the tomographic reconstruction techniques, particularly in high-latitude regions where GPS TEC measurements and other ionospheric data sources are limited.