Classical globally reflected gravity field determination in modern locally oriented multiscale framework

The purpose of this paper is the canonical connection of classical global gravity field determination following the concept of Stokes (Trans Camb Philos Soc 8:672–712, 1849), Bruns (Die Figur der Erde, Publikation Königl. Preussisch. Geodätisches Institut, P. Stankiewicz Buchdruckerei, Berlin, 1878), and Neumann (Vorlesungen über die Theorie des Potentials und der Kugelfunktionen. Teubner, Leipzig, pp 135–154, 1887) on the one hand and modern locally oriented multiscale computation by use of adaptive locally supported wavelets on the other hand. The essential tools are regularization methods of the Green, Neumann, and Stokes integral representations. The multiscale approximation is guaranteed simply as linear difference scheme by use of Green, Neumann, and Stokes wavelets. As an application, gravity anomalies caused by plumes are investigated for the Hawaiian and Iceland areas.     

On the topographic bias in geoid determination by the external gravity field

The topographic bias is defined as the error/bias committed by continuing the external gravity field inside the topographic masses by a harmonic function. We study the topographic bias given by a digital terrain model defined by a spherical template, and we show that the topographic bias is given only by the potential of an inner-zone cap, and it equals the bias of the Bouguer shell, independent of the size of the cap. Then we study the effect on the real Earth by decomposing its topography into a template, and we show also in this case that the topographic bias is that of the Bouguer shell, independent of the shape of the terrain. Finally, we show that the topographic potential of the terrain at the geoid can be determined to any precision by a Taylor expansion outside the Earth’s surface. The last statement is demonstrated by a Taylor expansion to fourth order.     

Long term shoreline oscillation and changes of Cauvery delta coastline inferred from satellite imageries

Coastal zone is highly volatile ecosystem which is always in adjustments. Loss of shore line will cause severe impact on human life and as well as their properties. Remote sensing is a reliable technique to study the historical shoreline changes. Therefore in this paper long term shoreline oscillations of Cauvery delta shorelines at Poompuhar, Tharangambadi and Nagapattinam were studied using satellite imageries and the same was physically observed at the above three locations with the help of reference pillars and compared mutually. It was observed that the shoreline at Poompuhar is under accretion at the rate of 1.79m/ year and other shoreline stretches at Tharangambadi and Nagapattinam were under erosion at 0.4888m/ year and 0.4985m/ year respectively. It was also observed that the remote sensing study qualitatively matches with the physical observation for all the three coastal stretches of the study area.     

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.     

Flood inundation modeling using MIKE FLOOD and remote sensing data

A coupled 1D-2D hydrodynamic model, MIKE FLOOD was used to simulate the flood inundation extent and flooding depth in the delta region of Mahanadi River basin in India. Initially, the 1D model MIKE 11 was calibrated using river water level and discharge data of various gauging sites for the monsoon period (June to October) of the year 2002. Subsequently, the calibrated set up was validated using both discharge and water level data for the same period of the year 2001. The performance of calibration and validation results of MIKE 11 were evaluated using different performance indices. A bathymetry of the study area with a spatial resolution of 90m was prepared from SRTM DEM and provided as an input to the 2D model, MIKE 21. MIKE 11 and MIKE 21 models were then coupled using lateral links to form the MIKE FLOOD model set up for simulating the two dimensional flood inundations in the study area. Flood inundation is simulated for the year 2001 and the maximum flood inundation extent simulated by the model was compared with the corresponding actual inundated area obtained from IRS-1D WiFS image.     

On the precision and accuracy of IGS orbits

In order to explore the precision and accuracy of International GNSS Service (IGS) orbits, we difference geocentric satellite positions midway between successive daily Final orbits for the period starting 5 November 2006, when the IGS switched its method of antenna calibration, through 31 December 2007. This yields a time series of orbit repeatabilities analogous to the classical geodetic test for position determinations. If we compare our average positional discontinuities to the official IGS accuracy codes, root-sum-squared (RSS) for each pair of days, we find the discontinuities are not well correlated with the predicted performance values. If instead the IGS weighted root-mean-square (WRMS) values from the Final combination long-arc analyses are taken as the measure of IGS accuracy, we find the position differences and long-arc values are correlated, but the long-arc values are exaggerated, particularly around eclipses, despite the fact that our day-boundary position differences apply to a single epoch each day and the long-arc analyses consider variations over a week. Our method is not well suited to probe the extent to which systematic effects dominate over random orbit errors, as indicated by satellite laser ranging residuals, but eclipsing satellites often display the most problematic behavior. A better metric than the current IGS orbit accuracy codes would probably be one based on the orbit discontinuities between successive days.     

The IGS VTEC maps: a reliable source of ionospheric information since 1998

The International GNSS Service (IGS) Working Group on Ionosphere was created in 1998. Since then, the Scientific community behind IGS, in particular CODE, ESA, JPL and UPC, have been continuosly contributing to reliable IGS combined vertical total electron content (VTEC) maps in both rapid and final schedules. The details on how these products are being generated, performance numbers, proposed improvement as far as VTEC evolution trends during near one Solar Cycle, are summarized in this paper. The confirmation of (1) the good performance of the IGS combined VTEC maps, and (2) the characteristic VTEC variability periods, are two main results of this work.     

CPLat: first operational experimental processing center for SIRGAS in Argentina

The SIRGAS permanent GPS network which is in fact the IGS network densification for the American continent, consists today of more than 200 stations covering the continent and islands. It is currently processed by the IGS RNAAC SIR centre at Deutsches Geodätisches Forschungsinstitut producing weekly free solutions relying on IGS final orbits and EOP that contribute to the ITRF through IGS. By August 2006, the SIRGAS Working Group I had accepted five proposals for experimental processing centers within the region that would collaborate with IGS RNAAC SIR. One of them, Centro de Procesamiento La Plata (CPLat) in Argentina, began processing 60 stations on October 2006. By January 2007 CPLat reached operational capability, delivering weekly free solution SINEX files, with an internal consistency of 1.5 mm average for the horizontal components, and 3 mm in the vertical. Comparisons with IGS global and IGS RNAAC SIR weekly solutions were taken as external consistency indications, showing average RMS residuals of 1.8, 2.4 and 5 mm for the north, east, and vertical component, respectively. Analysis and comparison of adjusted solution time series from CPLat and other processing centers has proved to be highly valuable for solution QC, namely detection and identification of station anomalous behavior or modelling problems. These procedures will ensure the maintenance of the performance specifications for CPLat solutions. Action is being taken in order to guarantee the continuity of this effort beyond the experimental phase.     

The IGS-combined station coordinates, earth rotation parameters and apparent geocenter

The International GNSS Service (IGS) routinely generates a number of weekly, daily and sub-daily products. Station coordinates and velocities, earth rotation parameters (ERPs) and apparent geocenter are among these products generated weekly by the IGS Reference Frame Coordinator. They have been determined since 1999 by combining independent estimates from at least seven IGS Analysis Centers (ACs). Two Global Network Associate Analysis Centers (GNAACs) also provide independent combinations using the same AC weekly solutions and they are currently used to quality control the IGS combination. The combined solutions are aligned to an IGS realization (IGS05) of the ITRF2005 using a carefully selected set of the IGS Reference Frame (RF) stations (nominally 132). During the combination process, the contributing solutions are compared and outliers are removed to ensure a high level of consistency of the estimated parameters. The ACs and the weekly combined solution are consistent at the 1–2 and 3–4 mm levels for the horizontal and vertical components. Similarly, the excess Length of Day (LOD), the pole positions and pole rates are consistent at the 10μs, 0.03–0.05 mas and 0.10–0.20 mas/day levels, respectively. The consistency of the apparent geocenter estimate is about 5 mm in the X and Y components and 10 mm in the Z component. Comparison of the IGS-combined ERP estimates with the IERS Bulletin A suggests a small bias of the order of ?0.04 mas and + 0.05 mas (both ±0.05 mas) in the x and y components.