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.     

Geostatistical incorporation of spatial coordinates into supervised classification of hyperspectral data

 This paper presents a methodology to incorporate both hyperspectral properties and spatial coordinates of pixels in maximum likelihood classification. Indicator kriging of ground data is used to estimate, for each pixel, the prior probabilities of occurrence of classes which are then combined with spectral-based probabilities within a Bayesian framework. In the case study (mapping of in-stream habitats), accounting for spatial coordinates increases the overall producer's accuracy from 85.8% to 93.8%, while the Kappa statistic rises from 0.74 to 0.88. Best results are obtained using only indicator kriging-based probabilities, with a stunning overall accuracy of 97.2%. Significant improvements are observed for environmentally important units, such as pools (Kappa: 0.17 to 0.74) and eddy drop zones (Kappa: 0.65 to 0.87). The lack of benefit of using hyperspectral information in the present study can be explained by the dense network of ground observations and the high spatial continuity of field classification which might be spurious. Received: 12 April 2001 / Accepted: 7 September 2001     

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.     

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.     

Land Cover Mapping for Namdapha National Park (Arunachal Pradesh), India Using Harmonized Land Cover Legends

Since last few decades RS-GIS is playing vital role in studying and mapping spatiotemporal responses of land cover, however, as a matter of fact, the mapping outputs largely depend on the expert's/user's preferences because location specific and people specific land cover classification systems are adopted autonomously for image classification in GIS. This may actually lead to an ambiguous definition of a particular land cover type when such different maps are compared at global level. In 1993, FAO and UNEP started efforts for development of a software tool know as LCCS which is a comprehensive standardized tool capable of providing land cover characterization to all possible land cover types in the world regardless of spatial relevance, mapping scale, data collection method etc. Adding to the global efforts of land cover legend harmonization and mapping, this study presents development of harmonized land cover legends for Namdapha National Park located in north-eastern Indian Himalayan region using LCCS and subsequent mapping. The potential of Remote Sensing (RS) and Geographical Information Systems (GIS) in forest/land cover mapping is very well recognized. Therefore, adopting the developed harmonized legends for the study area, land cover mapping was done using RS-GIS approach.     

Multifrequency algorithms for precise point positioning: MAP3

We present the new MAP3 algorithms to perform static precise point positioning (PPP) from multifrequency and multisystem GNSS observations. MAP3 represents a two-step strategy in which the least squares theory is applied twice to estimate smoothed pseudo-distances, initial phase ambiguities, and slant ionospheric delay first, and the absolute receiver position and its clock offset in a second adjustment. Unlike the classic PPP technique, in our new approach, the ionospheric-free linear combination is not used. The combination of signals from different satellite systems is accomplished by taking into account the receiver inter-system bias. MAP3 has been implemented in MATLAB and integrated within a complete PPP software developed on site and named PCube. We test the MAP3 performance numerically and contrast it with other external PPP programs. In general, MAP3 positioning accuracy with low-noise GPS dual-frequency observations is about 2.5 cm in 2-h observation periods, 1 cm in 10 h, and 7 mm after 1 day. This means an improvement in the accuracy in short observation periods of at least 7 mm with respect to the other PPP programs. The MAP3 convergence time is also analyzed and some results obtained from real triple-frequency GPS and GIOVE observations are presented.     

Non-linear motions of Australian geodetic stations induced by non-tidal ocean loading and the passage of tropical cyclones

We investigate daily and sub-daily non-tidal oceanic and atmospheric loading (NTOAL) in the Australian region and put an upper bound on potential site motion examining the effects of tropical cyclone Yasi that crossed the Australian coast in January/February 2011. The dynamic nature of the ocean is important, particularly for northern Australia where the long-term scatter due to daily and sub-daily oceanic changes increases by 20–55 % compared to that estimated using the inverted barometer (IB) assumption. Correcting the daily Global Positioning System (GPS) time series for NTOAL employing either a dynamic ocean model or the IB assumption leads to a reduction of up to 52 % in the weighted scatter of daily coordinate estimates. Differences between the approaches are obscured by seasonal variations in the GPS precision along the northern coast. Two compensating signals during the cyclone require modelling at high spatial and temporal resolution: uplift induced by the atmospheric depression, and subsidence induced by storm surge. The latter dominates ( \(>\) 135 %) the combined net effect that reaches a maximum of 14 mm, and 10 mm near the closest GPS site TOW2. Here, 96 % of the displacement is reached within 15 h due to the rapid transit of cyclones and the quasi-linear nature of the coastline. Consequently, estimating sub-daily NTOAL is necessary to properly account for such a signal that can be 3.5 times larger than its daily-averaged value. We were unable to detect the deformation signal in 2-hourly GPS processing and show that seasonal noise in the Austral summer dominates and precludes GPS detection of the cyclone-related subsidence.