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.     

Near real-time high-resolution airborne camera,AEROCam, for precision agriculture

Precision agriculture often relies on high-resolution imagery to delineate the variability within a field. Airborne Environmental Research Observational Camera (AEROCam) was designed to meet the needs of agriculture producers, ranchers, and researchers, who require high-resolution imagery in a near real-time environment for rapid decision support. AEROCam was developed and operated through a unique collaboration between several departments at the University of North Dakota, including the Upper Midwest Aerospace Consortium (UMAC), the School of Engineering and Mines, and flight operations at the John D. Odegard School of Aerospace Sciences. AEROCam consists of a Redlake MS4100 area-scan multi-spectral digital camera that features a 1920 × 1080 CCD array (7.4-μm detector) with 8-bit quantization. When operated at ～2 km above ground level, multispectral images with four bands in the visible and near infrared have a ground sample distance of 1 m with a horizontal extent of just over 1.6 km. Depending on the applications, flying at different altitudes can adjust the spatial resolution from 0.25 to 2 m. Rigorous spectral and radiometric calibrations allow AEROCam to be used in a variety of applications, qualitative and quantitative. Equipped with an inertial measurement unit (IMU) system, the images acquired can be geo-referenced automatically and delivered to end users near real time through our Digital Northern Great Plains system (DNGP). The images are also available to zone mapping application for precision farming (ZoneMAP), an online decision support tool for creating management zones from remote sensing imagery and data from other sources. Operational since 2004, AEROCam has flown over 250 sorties and delivered over 150,000 images to the users in the Northern Great Plains region, resulting in numerous applications in precision agriculture and resource management.     

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.     

A two-step ionospheric modeling algorithm considering the impact of GLONASS pseudo-range inter-channel biases

The Global Navigation Satellite System presents a plausible and cost-effective way of computing the total electron content (TEC). But TEC estimated value could be seriously affected by the differential code biases (DCB) of frequency-dependent satellites and receivers. Unlike GPS and other satellite systems, GLONASS adopts a frequency-division multiplexing access mode to distinguish different satellites. This strategy leads to different wavelengths and inter-frequency biases (IFBs) for both pseudo-range and carrier phase observations, whose impacts are rarely considered in ionospheric modeling. We obtained observations from four groups of co-stations to analyze the characteristics of the GLONASS receiver P1P2 pseudo-range IFB with a double-difference method. The results showed that the GLONASS P1P2 pseudo-range IFB remained stable for a period of time and could catch up to several meters, which cannot be absorbed by the receiver DCB during ionospheric modeling. Given the characteristics of the GLONASS P1P2 pseudo-range IFB, we proposed a two-step ionosphere modeling method with the priori IFB information. The experimental analysis showed that the new algorithm can effectively eliminate the adverse effects on ionospheric model and hardware delay parameters estimation in different space environments. During high solar activity period, compared to the traditional GPS + GLONASS modeling algorithm, the absolute average deviation of TEC decreased from 2.17 to 2.07 TECu (TEC unit); simultaneously, the average RMS of GPS satellite DCB decreased from 0.225 to 0.219 ns, and the average deviation of GLONASS satellite DCB decreased from 0.253 to 0.113 ns with a great improvement in over 55%.     

Tropospheric refractivity and zenith path delays from least-squares collocation of meteorological and GNSS data

Precise positioning requires an accurate a priori troposphere model to enhance the solution quality. Several empirical models are available, but they may not properly characterize the state of troposphere, especially in severe weather conditions. Another possible solution is to use regional troposphere models based on real-time or near-real time measurements. In this study, we present the total refractivity and zenith total delay (ZTD) models based on a numerical weather prediction (NWP) model, Global Navigation Satellite System (GNSS) data and ground-based meteorological observations. We reconstruct the total refractivity profiles over the western part of Switzerland and the total refractivity profiles as well as ZTDs over Poland using the least-squares collocation software COMEDIE (Collocation of Meteorological Data for Interpretation and Estimation of Tropospheric Pathdelays) developed at ETH Zürich. In these two case studies, profiles of the total refractivity and ZTDs are calculated from different data sets. For Switzerland, the data set with the best agreement with the reference radiosonde (RS) measurements is the combination of ground-based meteorological observations and GNSS ZTDs. Introducing the horizontal gradients does not improve the vertical interpolation, and results in slightly larger biases and standard deviations. For Poland, the data set based on meteorological parameters from the NWP Weather Research and Forecasting (WRF) model and from a combination of the NWP model and GNSS ZTDs shows the best agreement with the reference RS data. In terms of ZTD, the combined NWP-GNSS observations and GNSS-only data set exhibit the best accuracy with an average bias (from all stations) of 3.7 mm and average standard deviations of 17.0 mm w.r.t. the reference GNSS stations.     

Spherical integral transforms of second-order gravitational tensor components onto third-order gravitational tensor components

New spherical integral formulas between components of the second- and third-order gravitational tensors are formulated in this article. First, we review the nomenclature and basic properties of the second- and third-order gravitational tensors. Initial points of mathematical derivations, i.e., the second- and third-order differential operators defined in the spherical local North-oriented reference frame and the analytical solutions of the gradiometric boundary-value problem, are also summarized. Secondly, we apply the third-order differential operators to the analytical solutions of the gradiometric boundary-value problem which gives 30 new integral formulas transforming (1) vertical-vertical, (2) vertical-horizontal and (3) horizontal-horizontal second-order gravitational tensor components onto their third-order counterparts. Using spherical polar coordinates related sub-integral kernels can efficiently be decomposed into azimuthal and isotropic parts. Both spectral and closed forms of the isotropic kernels are provided and their limits are investigated. Thirdly, numerical experiments are performed to test the consistency of the new integral transforms and to investigate properties of the sub-integral kernels. The new mathematical apparatus is valid for any harmonic potential field and may be exploited, e.g., when gravitational/magnetic second- and third-order tensor components become available in the future. The new integral formulas also extend the well-known Meissl diagram and enrich the theoretical apparatus of geodesy.     

Group delay variations of GPS transmitting and receiving antennas

GPS code pseudorange measurements exhibit group delay variations at the transmitting and the receiving antenna. We calibrated C1 and P2 delay variations with respect to dual-frequency carrier phase observations and obtained nadir-dependent corrections for 32 satellites of the GPS constellation in early 2015 as well as elevation-dependent corrections for 13 receiving antenna models. The combined delay variations reach up to 1.0 m (3.3 ns) in the ionosphere-free linear combination for specific pairs of satellite and receiving antennas. Applying these corrections to the code measurements improves code/carrier single-frequency precise point positioning, ambiguity fixing based on the Melbourne–Wübbena linear combination, and determination of ionospheric total electron content. It also affects fractional cycle biases and differential code biases.