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     

Evaluation of pre-CHAMP gravity models `GRIM5-S1' and `GRIM5-C1' with satellite crossover altimetry

 The new GFZ/GRGS gravity field models GRIM5-S1 and GRIM5-C1, currently used as initial models for the CHAMP mission, have been compared with other recent models (JGM 3, EGM 96) for radial orbit accuracy (by means of latitude lumped coefficients) in computations on altimetry satellite orbits. The bases for accuracy judgements are multi-year averages of crossover sea height differences from Geosat and ERS 1/2 missions. This radially sensitive data is fully independent of the data used to develop these gravity models. There is good agreement between the observed differences in all of the world's oceans and projections of the same errors from the scaled covariance matrix of their harmonic geopotential coefficients. It was found that the tentative scale factor of five for the formal standard deviations of the harmonic coefficients of the new GRIM fields is justified, i.e. the accuracy estimates, provided together with the GRIM geopotential coefficients, are realistic. Received: 20 February 2001 / Accepted: 24 October 2001     

Conjoint modeling of residential group preferences: A comparison of the internal validity of hierarchical information integration approaches

 In this paper, two approaches for measuring residential group preferences, based on the method of Hierarchical Information Integration (HII), are compared. In particular, the hypothesis that group-based preference models estimated from integrated HII experiments better predict group preferences than part individual-based group models estimated from classical HII experiments is tested. To that effect, the models' ability to predict group preferences for new residential alternatives is compared in a study of residential preferences of co-ops. Results indicate that integrated HII group experiments indeed result in better predictions of residential preferences.     

GPS network monitors the Western Alps' deformation over a five-year period: 1993–1998

  The Western Alps are among the best studied collisional belts with both detailed structural mapping and also crustal geophysical investigations such as the ECORS and EGT seismic profile. By contrast, the present-day kinematics of the belt is still largely unknown due to small relative motions and the insufficient accuracy of the triangulation data. As a consequence, several tectonic problems still remain to be solved, such as the amount of N–S convergence in the Occidental Alps, the repartition of the deformation between the Alpine tectonic units, and the relation between deformation and rotation across the Alpine arc. In order to address these problems, the GPS ALPES group, made up of French, Swiss and Italian research organizations, has achieved the first large-scale GPS surveys of the Western Alps. More than 60 sites were surveyed in 1993 and 1998 with a minimum observation of 3 days at each site. GPS data processing has been done by three independent teams using different software. The different solutions have horizontal repeatabilities (N–E) of 4–7 mm in 1993 and 2–3 mm in 1998 and compare at the 3–5-mm level in position and 2-mm/yr level in velocity. A comparison of 1993 and 1998 coordinates shows that residual velocities of the GPS marks are generally smaller than 2 mm/yr, precluding a detailed tectonic interpretation of the differential motions. However, these data seem to suggest that the N–S compression of the Western Alps is quite mild (less than 2 mm/yr) compared to the global convergence between the African and Eurasian plate (6 mm/yr). This implies that the shortening must be accomodated elsewhere by the deformation of the Maghrebids and/or by rotations of Mediterranean microplates. Also, E–W velocity components analysis supports the idea that E–W extension exists, as already suggested by recent structural and seismotectonic data interpretation. Received: 27 November 2000 / Accepted: 17 September 2001     

Comparing seven candidate mission configurations for temporal gravity field retrieval through full-scale numerical simulation

The goal of this contribution is to focus on improving the quality of gravity field models in the form of spherical harmonic representation via alternative configuration scenarios applied in future gravimetric satellite missions. We performed full-scale simulations of various mission scenarios within the frame work of the German joint research project “Concepts for future gravity field satellite missions” as part of the Geotechnologies Program, funded by the German Federal Ministry of Education and Research and the German Research Foundation. In contrast to most previous simulation studies including our own previous work, we extended the simulated time span from one to three consecutive months to improve the robustness of the assessed performance. New is that we performed simulations for seven dedicated satellite configurations in addition to the GRACE scenario, serving as a reference baseline. These scenarios include a “GRACE Follow-on” mission (with some modifications to the currently implemented GRACE-FO mission), and an in-line “Bender” mission, in addition to five mission scenarios that include additional cross-track and radial information. Our results clearly confirm the benefit of radial and cross-track measurement information compared to the GRACE along-track observable: the gravity fields recovered from the related alternative mission scenarios are superior in terms of error level and error isotropy. In fact, one of our main findings is that although the noise levels achievable with the particular configurations do vary between the simulated months, their order of performance remains the same. Our findings show also that the advanced pendulums provide the best performance of the investigated single formations, however an accuracy reduced by about 2–4 times in the important long-wavelength part of the spectrum (for spherical harmonic degrees ${<}50$ ), compared to the Bender mission, can be observed. Concerning state-of-the-art mission constraints, in particular the severe restriction of heterodyne lasers on maximum range-rates, only the moderate Pendulum and the Bender-mission are beneficial options, of course in addition to GRACE and GRACE-FO. Furthermore, a Bender-type constellation would result in the most accurate gravity field solution by a factor of about 12 at long wavelengths (up to degree/order 40) and by a factor of about 200 at short wavelengths (up to degree/order 120) compared to the present GRACE solution. Finally, we suggest the Pendulum and the Bender missions as candidate mission configurations depending on the available budget and technological progress.     

Watershed Delineation in Flat Terrain of Thar Desert Region in North West India – A Semi Automated Approach Using DEM

Thar desert spreads in western part of Rajasthan, northern part of Gujarat, and some parts of Punjab and Haryana. The terrain is dominated by slightly sloping plains, broken by some dunes and low barren hills. The area is characterized by low average annual rainfall which is erratic in distribution and intensity. Drought will remain a major hindrance for agricultural production in Thar desert. Due to water stress condition, many watershed based development activities has been adopted by government and non-government organizations for the growth and sustainable development of this region. The need of this hour is preparation of a national level watershed atlas of 1:50,000 scale because majority of thematic maps are being produced presently on same or 1:10,000 scale. The manual delineation of watershed boundary in flat terrain based on topographic map will be time consuming and less accurate in the absence of prominent contour lines. Automated approach for watershed delineation using Digital Elevation Model (DEM) along a suitable algorithm has the advantage because the output is not only less time consuming but also independent from human decisions. Hence, a case study has been carried out in Churu sub-basin part of Indus basin which is located in Thar desert region. Depression less DEM with different spatial resolutions was used as input in hydrology tool of ArcGIS spatial analyst function for characterization of watersheds. The Churu sub-basin has been divided into various numbers of watersheds with an average size of 600 km². These watershed boundaries have been validated with respect to high resolution satellite imageries (IRS P6 LISS IV), Survey of India toposheets, ancillary data and limited field checks.     

Reducing the draconitic errors in GNSS geodetic products

Systematic errors at harmonics of the GPS draconitic year have been found in diverse GPS-derived geodetic products like the geocenter $Z$ -component, station coordinates, $Y$ -pole rate and orbits (i.e. orbit overlaps). The GPS draconitic year is the repeat period of the GPS constellation w.r.t. the Sun which is about 351 days. Different error sources have been proposed which could generate these spurious signals at the draconitic harmonics. In this study, we focus on one of these error sources, namely the radiation pressure orbit modeling deficiencies. For this purpose, three GPS+GLONASS solutions of 8 years (2004–2011) were computed which differ only in the solar radiation pressure (SRP) and satellite attitude models. The models employed in the solutions are: (1) the CODE (5-parameter) radiation pressure model widely used within the International GNSS Service community, (2) the adjustable box-wing model for SRP impacting GPS (and GLONASS) satellites, and (3) the adjustable box-wing model upgraded to use non-nominal yaw attitude, specially for satellites in eclipse seasons. When comparing the first solution with the third one we achieved the following in the GNSS geodetic products. Orbits: the draconitic errors in the orbit overlaps are reduced for the GPS satellites in all the harmonics on average 46, 38 and 57 % for the radial, along-track and cross-track components, while for GLONASS satellites they are mainly reduced in the cross-track component by 39 %. Geocenter $Z$ -component: all the odd draconitic harmonics found when the CODE model is used show a very important reduction (almost disappearing with a 92 % average reduction) with the new radiation pressure models. Earth orientation parameters: the draconitic errors are reduced for the $X$ -pole rate and especially for the $Y$ -pole rate by 24 and 50 % respectively. Station coordinates: all the draconitic harmonics (except the 2nd harmonic in the North component) are reduced in the North, East and Height components, with average reductions of 41, 39 and 35 % respectively. This shows, that part of the draconitic errors currently found in GNSS geodetic products are definitely induced by the CODE radiation pressure orbit modeling deficiencies.