Triangulated spherical splines for geopotential reconstruction

We present an alternate mathematical technique than contemporary spherical harmonics to approximate the geopotential based on triangulated spherical spline functions, which are smooth piecewise spherical harmonic polynomials over spherical triangulations. The new method is capable of multi-spatial resolution modeling and could thus enhance spatial resolutions for regional gravity field inversion using data from space gravimetry missions such as CHAMP, GRACE or GOCE. First, we propose to use the minimal energy spherical spline interpolation to find a good approximation of the geopotential at the orbital altitude of the satellite. Then we explain how to solve Laplace’s equation on the Earth’s exterior to compute a spherical spline to approximate the geopotential at the Earth’s surface. We propose a domain decomposition technique, which can compute an approximation of the minimal energy spherical spline interpolation on the orbital altitude and a multiple star technique to compute the spherical spline approximation by the collocation method. We prove that the spherical spline constructed by means of the domain decomposition technique converges to the minimal energy spline interpolation. We also prove that the modeled spline geopotential is continuous from the satellite altitude down to the Earth’s surface. We have implemented the two computational algorithms and applied them in a numerical experiment using simulated CHAMP geopotential observations computed at satellite altitude (450 km) assuming EGM96 (n _max = 90) is the truth model. We then validate our approach by comparing the computed geopotential values using the resulting spherical spline model down to the Earth’s surface, with the truth EGM96 values over several study regions. Our numerical evidence demonstrates that the algorithms produce a viable alternative of regional gravity field solution potentially exploiting the full accuracy of data from space gravimetry missions. The major advantage of our method is that it allows us to compute the geopotential over the regions of interest as well as enhancing the spatial resolution commensurable with the characteristics of satellite coverage, which could not be done using a global spherical harmonic representation. The results in this paper are based on the research supported by the National Science Foundation under the grant no. 0327577.     

A comparison of three image-object methods for the multiscale analysis of landscape structure

Within the conceptual framework of Complex Systems, we discuss the importance and challenges in extracting and linking multiscale objects from high-resolution remote sensing imagery to improve the monitoring, modeling and management of complex landscapes. In particular, we emphasize that remote sensing data are a particular case of the modifiable areal unit problem (MAUP) and describe how image-objects provide a way to reduce this problem. We then hypothesize that multiscale analysis should be guided by the intrinsic scale of the dominant landscape objects composing a scene and describe three different multiscale image-processing techniques with the potential to achieve this. Each of these techniques, i.e., Fractal Net Evolution Approach (FNEA), Linear Scale-Space and Blob-Feature Detection (SS), and Multiscale Object-Specific Analysis (MOSA), facilitates the multiscale pattern analysis, exploration and hierarchical linking of image-objects based on methods that derive spatially explicit multiscale contextual information from a single resolution of remote sensing imagery. We then outline the weaknesses and strengths of each technique and provide strategies for their improvement.     

Land damage assessment - a case study

The paper discusses the land damage assessment and change detection analysis with reference to a mineral bearing zone in Manjhi, Manuni and Churan valleys. The area is located in environmentally sensitive and fragile region of Himalaya and constitutes of nearly 400 small-scale mines of slate, which were operative since last one hundred years and are stopped by court of law since 1995 on account of environment deterioration. The status of land degradation has been studied using IRS-1B satellite data of 1988,1992 and 1995. The geo-coded data on 1:50,000 scale has been interpreted and an increase in land degradation status was noticed. Finally, the management strategy for arresting the further land damage in a broader perspective is suggested.     

Characterization of Compass M-1 signals

An analysis of observations from China’s first medium earth orbit satellite Compass M-1 is presented, with main focus on the first orbit and clock solution for this satellite. The orbit is computed from laser ranging measurements. Based on this orbit solution, the apparent clock offset is estimated using measurements from two GNSS receivers, which allow Compass tracking. The analysis of the clock solutions reveals unexpectedly high dynamics in the pseudorange and carrier-phase observations. Furthermore, carrier-to-noise density ratio, pseudorange noise, and multipath are analyzed and compared to GPS and GIOVE. The results of the clock analysis motivate further research on the signals of the geostationary satellites of the Compass constellation.     

Hyperspectral remote sensing of evaporate minerals and associated sediments in Lake Magadi area, Kenya

Pleistocene to present evaporitic lacustrine sediments in Lake Magadi, East African Rift Valley, Kenya were studied and mapped using spectral remote sensing methods. This approach incorporated surface mineral mapping using space-borne hyperspectral Hyperion imagery together with laboratory analysis, including visible, near-infrared diffuse reflectance spectroscopy (VNIR) measurements and X-ray diffraction for selected rock and soil samples of the study area. The spectral signatures of Magadiite and Kenyaite, which have not been previously reported, were established and the spectral signatures of trona, chert series, volcanic tuff and the High Magadi bed were also analyzed.Image processing techniques, MNF (Minimum Noise Fraction) and MTMF (Mixture Tuned Matched Filtering) using a stratified approach (image analysis with and without the lake area), were used to enhance the mapping of evaporates. High Magadi beds, chert series and volcanic tuff were identified from the Hyperion image with an overall mapping accuracy of 84.3%. Even though, the spatial distribution of evaporites and sediments in Lake Magadi area change in response to climate variations, the mineralogy of this area has not been mapped recently. The results of this study shows the usefulness of the hypersspectral remote sensing to map the surface geology of this kind of environment and to locate promising sites for industrial open-pit trona mining in a qualitative and quantitative manner.     

The impact of accelerometry on CHAMP orbit determination

 The contribution of the STAR accelerometer to the CHAMP orbit precision is evaluated and quantified by means of the following results: orbital fit to the satellite laser ranging (SLR) observations, GPS reduced-dynamic vs SLR dynamic orbit comparisons, and comparison of the measured to the modeled non-gravitational accelerations (atmospheric drag in particular). In each of the four test periods in 2001, five CHAMP arcs of 2 days' length were analyzed. The mean RMS-of-fit of the SLR observations of the orbits computed with STAR data or the non-gravitational force model were 11 and 24 cm, respectively. If the accelerometer calibration parameters are not known at least at the few percent level, the SLR orbit fit deteriorates. This was tested by applying a 10% error to the along-track scale factor of the accelerometer, which increased the SLR RMS-of-fit on average to 17 cm. Reference orbits were computed employing the reduced-dynamic technique with GPS tracking data. This technique yields the most accurate orbit positions thanks to the estimation of a large number of empirical accelerations, which compensate for dynamic modeling errors. Comparison of the SLR orbits, computed with STAR data or the non-gravitational force model, to the GPS-based orbits showed that the SLR orbits employing accelerometer observations are twice as accurate. Finally, comparison of measured to modeled accelerations showed that the level of geomagnetic activity is highly correlated with the atmospheric drag model error, and that the largest errors occur around the geomagnetic poles. Received: 7 May 2002 / Accepted: 18 November 2002 Correspondence to: S. Bruinsma Acknowledgments. The TIGCM results were obtained from the CEDAR database. This study was supported by the Centre National d'Etudes Spatiales (CNES). The referees are thanked for their helpful remarks and suggestions.     

Groebner-basis solution of the three-dimensional resection problem (P4P)

The three-dimensional (3-D) resection problem is usually solved by first obtaining the distances connecting the unknown point P{X,Y,Z} to the known points P_i{X_i,Y_i,Z_i}i=1,2,3 through the solution of the three nonlinear Grunert equations and then using the obtained distances to determine the position {X,Y,Z} and the 3-D orientation parameters {,, }. Starting from the work of the German J. A. Grunert (1841), the Grunert equations have been solved in several substitutional steps and the desire as evidenced by several publications has been to reduce these number of steps. Similarly, the 3-D ranging step for position determination which follows the distance determination step involves the solution of three nonlinear ranging (`Bogenschnitt') equations solved in several substitution steps. It is illustrated how the algebraic technique of Groebner basis solves explicitly the nonlinear Grunert distance equations and the nonlinear 3-D ranging (`Bogenschnitt') equations in a single step once the equations have been converted into algebraic (polynomial) form. In particular, the algebraic tool of the Groebner basis provides symbolic solutions to the problem of 3-D resection. The various forward and backward substitution steps inherent in the classical closed-form solutions of the problem are avoided. Similar to the Gauss elimination technique in linear systems of equations, the Groebner basis eliminates several variables in a multivariate system of nonlinear equations in such a manner that the end product normally consists of a univariate polynomial whose roots can be determined by existing programs e.g. by using the roots command in Matlab.Acknowledgments.The first author wishes to acknowledge the support of JSPS (Japan Society of Promotion of Science) for the financial support that enabled the completion of the write-up of the paper at Kyoto University, Japan. The author is further grateful for the warm welcome and the good working atmosphere provided by his hosts Professors S. Takemoto and Y. Fukuda of the Department of Geophysics, Graduate School of Science, Kyoto University, Japan.     

A formula for computing the gravity disturbance from the second radial derivative of the disturbing potential

 A formula for computing the gravity disturbance and gravity anomaly from the second radial derivative of the disturbing potential is derived in detail using the basic differential equation with spherical approximation in physical geodesy and the modified Poisson integral formula. The derived integral in the space domain, expressed by a spherical geometric quantity, is then converted to a convolution form in the local planar rectangular coordinate system tangent to the geoid at the computing point, and the corresponding spectral formulae of 1-D FFT and 2-D FFT are presented for numerical computation. Received: 27 December 2000 / Accepted: 3 September 2001