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     

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.     

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.     

Use of the L2C signal for inversions of GPS radio occultation data in the neutral atmosphere

Results from processing FORMOSAT-3/COSMIC radio occultations (RO) with the new GPS L2C signal acquired both in phase locked loop (PLL) and open loop (OL) modes are presented. Analysis of L2P, L2C, and L1CA signals acquired in PLL mode shows that in the presence of strong ionospheric scintillation not only L2P tracking, but also L1CA tracking often fails, while L2C tracking is most stable. The use of L2C improves current RO processing in the neutral atmosphere mainly by increasing the number of processed occultations (due to significant reduction in the number of L2 tracking failures) and marginally by a reduction in noise in statistics. The latter is due to the combination of reduced L2C noise (compared to L2P) and increased L1CA noise in those occultations where L2P would have failed. This result suggests application of OL tracking for L1CA and L2C signals throughout an entire occultation to optimally acquire RO data. Two methods of concurrent processing of L1CA and L2C RO signals are considered. Based on testing of individual occultations, these methods allow: (1) reduction in uncertainty of bending angles retrieved by wave optics in the lower troposphere and (2) reduction in small-scale residual errors of the ionospheric correction in the stratosphere.     

Comparative evaluation of Indian remote sensing satellite and landsat thematic mapper data for geological and geomorphological applications

In order to evaluate the potentials of IRS‐1A Linear Imaging Self‐scanning Sensor (LISS‐I) data for geological and geomorphological applications and also to compare the IRS‐1A LISS‐I data with Landsat Thematic Mapper (TM) data, a study has been attempted for parts of Uttar Pradesh and Madhya Pradesh in Northern India. The first four spectral bands of Landsat TM sensor data which are similar and close to IRS‐1A LISS‐I senor have been utilised for the comparative evaluation. Various techniques employed for both the data set to derive the required geology and geomorphology related information include (i) band combination (ii) spectral response analysis (iii) principal component analysis (iv) supervised classification techniques and (v) visual observation of various outputs generated by the above methods. The Optimum Index Factor (OIF) method adopted for selecting suitable band combinations showed similar OIF rankings for IRS‐1A LISS‐I data and Landsat TM data. It has been visually observed that the band combination 1, 3 & 4 offers relatively better feature display. The spectral responses derived for various major geologic rock units such as Deccan Trap, Vindhyan Formation, Bundelkhand Granite and for a few landcovers such as surface water bodies and black soil show striking similarity in pattern for both LISS‐I and TM. The Principal Component (PC) analysis of both data sets suggested that the total scene brightness tends to dominate in the first PC. The percentage information contributed by PCs 1&2 as also by PCs 1,2 & 3 in both the LISS‐I and TM are comparable. It was observed from the classified image generated by performing supervised classification with a maximum likelihood algorithm that major geomorphic landforms were clearly distinguishable. Thus the qualitative and quantitative evaluation of both IRS‐1A LISS‐I and Landsat TM data showed that significant similarities exist between them. The study also revealed that IRS‐1A LISS‐I data can be effectively used for deriving geology and geomorphology related details.     

Mapping spiny aster infestations with QuickBird imagery

QuickBird satellite imagery acquired in June 2003 and September 2004 was evaluated for detecting the noxious weed spiny aster [Leucosyris spinosa (Benth.) Greene] on a south Texas, USA rangeland area. A subset of each of the satellite images representing a diversity of cover types was extracted and used as a study site. The satellite imagery had a spatial resolution of 2.8 m and contained 11-bit data. Unsupervised and supervised classification techniques were used to classify false colour composite (green, red, and near-infrared bands) images of the study site. Imagery acquired in June was superior to that obtained in September for distinguishing spiny aster infestations. This was attributed to differences in spiny aster phenology between the two dates. An unsupervised classification of the June image showed that spiny aster had producer's and user's accuracies of 90% and 93.1%, respectively, whereas a supervised classification of the June image had producer's and user's accuracies of 90% and 81.8%, respectively. These results indicate that high resolution satellite imagery coupled with image analysis techniques can be used successfully for detecting spiny aster infestations on rangelands.     

A global inventory of coral reef stressors based on satellite observed nighttime lights

In this article, we present a satellite-based approach to gather information about the threat to coral reefs worldwide. Three chosen reef stressors – development, gas flaring and heavily lit fishing boat activity – are analysed using nighttime lights data derived from the Defense Meteorological Satellite Program (DMSP) produced at the National Oceanic & Atmospheric Administration, National Geophysical Data Center (NOAA/NGDC). Nighttime lights represent a direct threat to coral reef ecosystems and are an excellent proxy measure for associated human-caused stressors. A lights proximity index (LPI) is calculated, measuring the distance of coral reef sites to each of the stressors and incorporating the stressor's intensity. Colourized maps visualize the results on a global scale. Area rankings clarify the effects of artificial night lighting on coral reefs on a regional scale. The results should be very useful for reef managers and for state administrations to implement coral reef conservation projects and for the scientific world to conduct further research.     

Generation and study of satellite gravity over Gujarat,India and their possible correlation with earthquake occurrences

High-resolution satellite gravity data have been generated and utilized to infer subsurface geological structures in the area of devastating earthquake that struck the Bhuj region in Gujarat on 26 January 2001. Latitudinal gravity profiles have been generated in the Bhuj, Anjar and IBF regions across the epicentres (23.5° N, 69.8° E/M_w 7.0 in 2001; 23.2° N, 70° E/M_w 7.0 in 1956; 24.2° N, 69.2° E/M_w 7.8 in 1819). Substantial differences in gravity anomaly patterns as high as 37 mGal could be observed existing near the epicentre regions. These gravitational differences might have caused due to the plate tectonic processes and due to the changes in densities of different lithospheric zones/sedimentary layers. Temporal variations of the satellite-derived gravity and their probable relations with already occurred major earthquakes in this region have been studied. Hence we conclude that drastic changes in gravity anomalies can be considered as a precursor for occurrences of substantially large earthquakes.