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.     

Stationary spatially modulated fourier transform spectro-radiometer

Journal of the Indian Society of Remote Sensing - A stationary, compact, spatially modulated Fourier Transform spectro-radiometer based on triangular, common path Sagnac interferometer has been...     

Change detection in Sal forest in Dehradun Forest Division using remote sensing and geographical information system

The review of study site have revealed the change in vegetation cover of Sal Dense to Sal Medium and Sal Open in 6 forest Mosaics owing to biotic and abiotic conditions prevailing in the specific areas. Analysis carried out using thematic map derived from aerial photograph of 1976 and satellite data of IRS 1C LISS III False Colour Composite (FCC) of March 1999 revealed the cause for change in forest density classes. Deforestation, encroachment and agriculture have been identified as the underlying causes, which have affected some specific locations to a marked extent. There has been a progressive and remarkable change among vegetation classes from 1976 to 1999. It is evident from forest type and density map that Sal density has significantly reduced from Sal Dense 65.61 % in 1976 to Sal Dense 11.12% in the year 1999 followed by Sal Open 11.18 % and Sal Medium 18.24 %. The overall change has been estimated to be 42.11% of the total forested area.     

Remote Sensing Study of Sittampundi Anorthosite Complex,India

In the present study, The Landsat 7 ETM satellite data was collected for the Sittampundi anorthosites complex and digital image analysis was carried out. The anorthositic rocks available at Sittampundi complex is considered as an equivalent of lunar highland rocks. Hence, a remote sensing study comprises of image analysis and spectral profile analysis was carried out. The satellite data was digitally processed and generated various outputs like band combinations, color composites, stretched outputs, and PCA. The suitable processed outputs were identified for delineating the anorthosite complex. The diagnostic absorption features of reflectance spectra are the sensitive indicators of mineralogy and chemical composition of rocks, which are interest to the planetary scientists. The spectral profile of Landsat ETM plotted for pure and mixed anorthosite pixels and compared with the field and lab reflectance spectra. The percentages of image spectra vary from 30% to 60% for Sittampundi anorthosite. The spectral bands 2, 4 and 6 have low reflectance and bands 3 and 5 have high reflectance. The spectral range of bands 2,3,4,5 and 6 are 525 nm–605 nm, 630 nm–690 nm, 750 nm–900 nm, 1550 nm–1750 nm and 10400 nm–12500 nm respectively. The field spectral curve has weak absorptions at 650 nm and 1000 nm due to the iron transition absorption and low ca- pyroxene respectively available in the anorthosite, matching with the image spectra. However, hyperspectal image with narrow bandwidth could be more useful in selecting the suitable spectrum for remotely mapping the anorthosite region, as equivalent test site for lunar highland region.     

Artificial Neural Network (ANN) Based Inversion of Benthic Substrate Bottom Type and Bathymetry in Optically Shallow Waters - Initial Model Results

Ocean-colour remote sensing in optically shallow waters is influenced by contribution from the water column depth as well as by the substrate type. Therefore, it is required to include the contribution from the water column and substrate bottom type for bathymetry estimation. In this report we demonstrate the use of Artificial Neural Network (ANN) based approach to spectrally distinguish various benthic bottom types and estimate depth of substrate bottom simultaneously in optically shallow waters. We have used in-water radiative transfer simulation modeling to generate simulated top-of-the-water column reflectance the four major benthic bottom types viz. sea grass, coral sand, green algae and red algae using Hydrolight simulation model. The simulated remote sensing reflectance, for the four benthic bottom types having benthic bottom depth up to 30 m were generated for moderately clear waters. A multi-layer perceptron (MLP) type neural network was trained using the simulated data. ANN based approach was used for classification of the benthic bottom type and simultaneous inversion of bathymetry. Simulated data was inverted to yield benthic bottom type classification with an accuracy of ~98% for the four benthic substrate types and the substrate depth were estimated with an error of 0% for sea grass, 1% for coral sand and 1–3% for green and red algae up to 25 m, whereas for substrate bottom deeper than 25 m depth the classification errors increased by 2–5% for three substrate bottom types except sea grass bottom type. The initial results are promising which needs validation using the in-situ measured remote sensing reflectance spectra for implementing further on satellite data.     

Simulation study of a follow-on gravity mission to GRACE

The gravity recovery and climate experiment (GRACE) has been providing monthly estimates of the Earth’s time-variable gravity field since its launch in March 2002. The GRACE gravity estimates are used to study temporal mass variations on global and regional scales, which are largely caused by a redistribution of water mass in the Earth system. The accuracy of the GRACE gravity fields are primarily limited by the satellite-to-satellite range-rate measurement noise, accelerometer errors, attitude errors, orbit errors, and temporal aliasing caused by un-modeled high-frequency variations in the gravity signal. Recent work by Ball Aerospace & Technologies Corp., Boulder, CO has resulted in the successful development of an interferometric laser ranging system to specifically address the limitations of the K-band microwave ranging system that provides the satellite-to-satellite measurements for the GRACE mission. Full numerical simulations are performed for several possible configurations of a GRACE Follow-On (GFO) mission to determine if a future satellite gravity recovery mission equipped with a laser ranging system will provide better estimates of time-variable gravity, thus benefiting many areas of Earth systems research. The laser ranging system improves the range-rate measurement precision to ~0.6 nm/s as compared to ~0.2 μm/s for the GRACE K-band microwave ranging instrument. Four different mission scenarios are simulated to investigate the effect of the better instrument at two different altitudes. The first pair of simulated missions is flown at GRACE altitude (~480 km) assuming on-board accelerometers with the same noise characteristics as those currently used for GRACE. The second pair of missions is flown at an altitude of ~250 km which requires a drag-free system to prevent satellite re-entry. In addition to allowing a lower satellite altitude, the drag-free system also reduces the errors associated with the accelerometer. All simulated mission scenarios assume a two satellite co-orbiting pair similar to GRACE in a near-polar, near-circular orbit. A method for local time variable gravity recovery through mass concentration blocks (mascons) is used to form simulated gravity estimates for Greenland and the Amazon region for three GFO configurations and GRACE. Simulation results show that the increased precision of the laser does not improve gravity estimation when flown with on-board accelerometers at the same altitude and spacecraft separation as GRACE, even when time-varying background models are not included. This study also shows that only modest improvement is realized for the best-case scenario (laser, low-altitude, drag-free) as compared to GRACE due to temporal aliasing errors. These errors are caused by high-frequency variations in the hydrology signal and imperfections in the atmospheric, oceanographic, and tidal models which are used to remove unwanted signal. This work concludes that applying the updated technologies alone will not immediately advance the accuracy of the gravity estimates. If the scientific objectives of a GFO mission require more accurate gravity estimates, then future work should focus on improvements in the geophysical models, and ways in which the mission design or data processing could reduce the effects of temporal aliasing.     

Evaluation of Almond (Amygdalus Scopartia) Habitat Using Multi-criteria Decision Making Method

Bolhasan Forest region with an area of 5,725?ha is located on east north of Dezful County, Iran. The region belongs to natural forests of Dezful. Considering the area is mainly covered by valuable species of Amygdalus Scopartia, its sustainable exploitation and development as well as restoration enjoys great importance. Study ahead aims at selection of suitable habitats for under studied species using Analytical Hierarchy Process (AHP). Therewith, the required thematic maps were imported in to GIS Software and final suitability map was prepared. The results indicated that around 2,119?ha (37%) out of all study area has high suitability for habitat of Amygdalus Scopartia. In the meanwhile, 1,603?ha [equal to 28%] is categorized as good suitability class and 2,003?ha [35%] has poor suitability.     

A bilinear approximation of the surface relief in terrain correction computations

A new method is presented for the computation of the gravitational attraction of topographic masses when their height information is given on a regular grid. It is shown that the representation of the terrain relief by means of a bilinear surface not only offers a serious alternative to the polyhedra modeling, but also approaches even more smoothly the continuous reality. Inserting a bilinear approximation into the known scheme of deriving closed analytical expressions for the potential and its first-order derivatives for an arbitrarily shaped polyhedron leads to a one-dimensional integration with – apparently – no analytical solution. However, due to the high degree of smoothness of the integrand function, the numerical computation of this integral is very efficient. Numerical tests using synthetic data and a densely sampled digital terrain model in the Bavarian Alps prove that the new method is comparable to or even faster than a terrain modeling using polyhedra.     

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.