The determination of potential difference by the joint application of measured and synthetical gravity data: a case study in Hungary

In an elementary approach every geometrical height difference between the staff points of a levelling line should have a corresponding average g value for the determination of potential difference in the Earth’s gravity field. In practice this condition requires as many gravity data as the number of staff points if linear variation of g is assumed between them. Because of the expensive fieldwork, the necessary data should be supplied from different sources. This study proposes an alternative solution, which is proved at a test bed located in the Mecsek Mountains, Southwest Hungary, where a detailed gravity survey, as dense as the staff point density (~1 point/34 m), is available along a 4.3-km-long levelling line. In the first part of the paper the effect of point density of gravity data on the accuracy of potential difference is investigated. The average g value is simply derived from two neighbouring g measurements along the levelling line, which are incrementally decimated in the consecutive turns of processing. The results show that the error of the potential difference between the endpoints of the line exceeds 0.1 mm in terms of length unit if the sampling distance is greater than 2 km. Thereafter, a suitable method for the densification of the decimated g measurements is provided. It is based on forward gravity modelling utilising a high-resolution digital terrain model, the normal gravity and the complete Bouguer anomalies. The test shows that the error is only in the order of 10⁻³mm even if the sampling distance of g measurements is 4 km. As a component of the error sources of levelling, the ambiguity of the levelled height difference which is the Euclidean distance between the inclined equipotential surfaces is also investigated. Although its effect accumulated along the test line is almost zero, it reaches 0.15 mm in a 1-km-long intermediate section of the line.     

Phase center modeling for LEO GPS receiver antennas and its impact on precise orbit determination

Most satellites in a low-Earth orbit (LEO) with demanding requirements on precise orbit determination (POD) are equipped with on-board receivers to collect the observations from Global Navigation Satellite systems (GNSS), such as the Global Positioning System (GPS). Limiting factors for LEO POD are nowadays mainly encountered with the modeling of the carrier phase observations, where a precise knowledge of the phase center location of the GNSS antennas is a prerequisite for high-precision orbit analyses. Since 5 November 2006 (GPS week 1400), absolute instead of relative values for the phase center location of GNSS receiver and transmitter antennas are adopted in the processing standards of the International GNSS Service (IGS). The absolute phase center modeling is based on robot calibrations for a number of terrestrial receiver antennas, whereas compatible antenna models were subsequently derived for the remaining terrestrial receiver antennas by conversion (from relative corrections), and for the GNSS transmitter antennas by estimation. However, consistent receiver antenna models for space missions such as GRACE and TerraSAR-X, which are equipped with non-geodetic receiver antennas, are only available since a short time from robot calibrations. We use GPS data of the aforementioned LEOs of the year 2007 together with the absolute antenna modeling to assess the presently achieved accuracy from state-of-the-art reduced-dynamic LEO POD strategies for absolute and relative navigation. Near-field multipath and cross-talk with active GPS occultation antennas turn out to be important and significant sources for systematic carrier phase measurement errors that are encountered in the actual spacecraft environments. We assess different methodologies for the in-flight determination of empirical phase pattern corrections for LEO receiver antennas and discuss their impact on POD. By means of independent K-band measurements, we show that zero-difference GRACE orbits can be significantly improved from about 10 to 6 mm K-band standard deviation when taking empirical phase corrections into account, and assess the impact of the corrections on precise baseline estimates and further applications such as gravity field recovery from kinematic LEO positions.     

Using gravity and topography-implied anomalies to assess data requirements for precise geoid computation

Many regions around the world require improved gravimetric data bases to support very accurate geoid modeling for the modernization of height systems using GPS. We present a simple yet effective method to assess gravity data requirements, particularly the necessary resolution, for a desired precision in geoid computation. The approach is based on simulating high-resolution gravimetry using a topography-correlated model that is adjusted to be consistent with an existing network of gravity data. Analysis of these adjusted, simulated data through Stokes’s integral indicates where existing gravity data must be supplemented by new surveys in order to achieve an acceptable level of omission error in the geoid undulation. The simulated model can equally be used to analyze commission error, as well as model error and data inconsistencies to a limited extent. The proposed method is applied to South Korea and shows clearly where existing gravity data are too scarce for precise geoid computation.     

A hierarchical model for estimating methane emission from wetlands using MODIS data and ARIMA modeling

A three-step hierarchical Semi Automated Empirical Methane Emission Model (SEMEM) has been used to estimate methane emission from wetlands and waterlogged areas in India using Moderate Resolution Imagine Spectroradiometer (MODIS) sensor data onboard Terra satellite. Wetland Surface Temperature (WST), methane emission fluxes and wetland extent have been incorporated as parameters in order to model the methane emission. Analysis of monthly MODIS data covering the whole of India from November 2004 to April 2006 was carried out and monthly methane emissions have been estimated. Interpolation techniques were adopted to fill the data gaps due to cloudy conditions during the monsoon period. AutoRegressive Integrated Moving Average (ARIMA) model has been fitted to estimate the emitted methane for the months of May 2006 to August 2006 using SPSS software.     

HPABM: A Hierarchical Parallel Simulation Framework for Spatially-explicit Agent-based Models

A Hierarchical Parallel simulation framework for spatially-explicit Agent-Based Models (HPABM) is developed to enable computationally intensive agent-based models for the investigation of large-scale geospatial problems. HPABM allows for the utilization of high-performance and parallel computing resources to address computational challenges in agent-based models. Within HPABM, an agent-based model is decomposed into a set of sub-models that function as computational units for parallel computing. Each sub-model is comprised of a sub-set of agents and their spatially-explicit environments. Sub-models are aggregated into a group of super-models that represent computing tasks. HPABM based on the design of super- and sub-models leads to the loose coupling of agent-based models and underlying parallel computing architectures. The utility of HPABM in enabling the development of parallel agent-based models was examined in a case study. Results of computational experiments indicate that HPABM is scalable for developing large-scale agent-based models and, thus, demonstrates efficient support for enhancing the capability of agent-based modeling for large-scale geospatial simulation.     

Water Harvesting Structure Positioning by Using Geo-Visualization Concept and Prioritization of Mini-Watersheds Through Morphometric Analysis in the Lower Tapi Basin

Geo-visualization concept has been used for positioning water harvesting structures in Varekhadi watershed consisting of 26 mini watersheds, falling in Lower Tapi Basin (LTB), Surat district, Gujarat state. For prioritization of the mini watersheds, morphometric analysis was utilized by using the linear parameters such as bifurcation ratio (R_b), drainage density (D_d), stream frequency (F_u), texture ratio (T), length of overland flow (L_o) and the shape parameter such as form factor (R_f), shape factor (B_s), elongation ratio (R_e), compactness constant (C_c) and circularity ratio (R_c). The different prioritization ranks were assigned after evaluation of the compound factor. 3 Dimensional (3D) Elevation Model (DEM) from Shuttle Radar Topography Mission (SRTM) and DEM from topo contour were analyzed in ArcScene 9.1 and the fly tool was utilized for the Geo-visualization of Varekhadi mini watersheds as per the priority ranks. Combining this with soil map and slope map, the best feasibility of positioning check dams in mini-watershed no. 1, 5 and 24 has been proposed, after validation of the sites.     

Geostatistical approaches to refinement of digital elevation data

Data refinement refers to the processes by which a dataset’s resolution, in particular, the spatial one, is refined, and is thus synonymous to spatial downscaling. Spatial resolution indicates measurement scale and can be seen as an index for regular data support. As a type of change of scale, data refinement is useful for many scenarios where spatial scales of existing data, desired analyses, or specific applications need to be made commensurate and refined. As spatial data are related to certain data support, they can be conceived of as support-specific realizations of random fields, suggesting that multivariate geostatistics should be explored for refining datasets from their coarser-resolution versions to the finer-resolution ones. In this paper, geostatistical methods for downscaling are described, and were implemented using GTOPO30 data and sampled Shuttle Radar Topography Mission data at a site in northwest China, with the latter’s majority grid cells used as surrogate reference data. It was found that proper structural modeling is important for achieving increased accuracy in data refinement; here, structural modeling can be done through proper decomposition of elevation fields into trends and residuals and thereafter. It was confirmed that effects of semantic differences on data refinement can be reduced through properly estimating and incorporating biases in local means.     

Spatial-Correlation Based Persistent Scatterer Interferometric Study for Ground Deformation

Interferometric Synthetic Aperture Radar (InSAR), nowadays, is a precise technique for monitoring and detecting ground deformation at a millimetric level over large areas using multi-temporal SAR images. Persistent Scatterer Interferometric SAR (PSInSAR), an advanced version of InSAR, is an effective tool for measuring ground deformation using temporally stable reference points or persistent scatterers. We have applied both PSInSAR and Small Baseline Subset (SBAS) methods, based on the spatial correlation of interferometric phase, to estimate the ground deformation and time-series analysis. In this study, we select Las Vegas, Nevada, USA as our test area to detect the ground deformation along satellite line-of-sight (LOS) during November 1992–September 2000 using 44 C-band SAR images of the European Remote Sensing (ERS-1 and ERS-2) satellites. We observe the ground displacement rate of Las Vegas is in the range of ?19 to 8 mm/year in the same period. We also cross-compare PSInSAR and SBAS using mean LOS velocity and time-series. The comparison shows a correlation coefficient of 0.9467 in the case of mean LOS velocity. Along this study, we validate the ground deformation results from the satellite with the ground water depth of Las Vegas using time-series analysis, and the InSAR measurements show similar patterns with ground water data.     

Ionospheric and receiver DCB-constrained multi-GNSS single-frequency PPP integrated with MEMS inertial measurements

Single-frequency precise point positioning (SF-PPP) is a potential precise positioning technique due to the advantages of the high accuracy in positioning after convergence and the low cost in operation. However, there are still challenges limiting its applications at present, such as the long convergence time, the low reliability, and the poor satellite availability and continuity in kinematic applications. In recent years, the achievements in the dual-frequency PPP have confirmed that its performance can be significantly enhanced by employing the slant ionospheric delay and receiver differential code bias (DCB) constraint model, and the multi-constellation Global Navigation Satellite Systems (GNSS) data. Accordingly, we introduce the slant ionospheric delay and receiver DCB constraint model, and the multi-GNSS data in SF-PPP modular together. In order to further overcome the drawbacks of SF-PPP in terms of reliability, continuity, and accuracy in the signal easily blocking environments, the inertial measurements are also adopted in this paper. Finally, we form a new approach to tightly integrate the multi-GNSS single-frequency observations and inertial measurements together to ameliorate the performance of the ionospheric delay and receiver DCB-constrained SF-PPP. In such model, the inter-system bias between each two GNSS systems, the inter-frequency bias between each two GLONASS frequencies, the hardware errors of the inertial sensors, the slant ionospheric delays of each user-satellite pair, and the receiver DCB are estimated together with other parameters in a unique Kalman filter. To demonstrate its performance, the multi-GNSS and low-cost inertial data from a land-borne experiment are analyzed. The results indicate that visible positioning improvements in terms of accuracy, continuity, and reliability can be achieved in both open-sky and complex conditions while using the proposed model in this study compared to the conventional GPS SF-PPP.