Approaches for comparative evaluation of raster GIS-based landslide susceptibility zonation maps

Evaluation of maps generated from different conceptual models or data processing approaches at spatial level has importance in many geoenvironmental applications. This paper addresses the spatial comparison of different landslide susceptibility zonation (LSZ) raster maps of the same area derived from various procedures.     

Using spectral information from the NIR water absorption features for the retrieval of canopy water content

Canopy water content (CWC) is important for mapping and monitoring the condition of the terrestrial ecosystem. Spectral information related to the water absorption features at 970 nm and 1200 nm offers possibilities for deriving information on CWC. In this study, we compare the use of derivative spectra, spectral indices and continuum removal techniques for these regions. Hyperspectral reflectance data representing a range of canopies were simulated using the combined PROSPECT + SAILH model. Best results in estimating CWC were obtained by using spectral derivatives at the slopes of the 970 nm and 1200 nm water absorption features. Real data from two different test sites were analysed. Spectral information at both test sites was obtained with an ASD FieldSpec spectrometer, whereas at the second site HyMap airborne imaging spectrometer data were also acquired. Best results were obtained for the derivative spectra. In order to avoid the potential influence of atmospheric water vapour absorption bands the derivative of the reflectance on the right slope of the canopy water absorption feature at 970 nm can best be used for estimating CWC.     

Application of remote sensing technology for Impact Assessment of Watershed Development Programme

The current study was taken up to investigate the utility of remote sensing and GIS tools for evaluation of Integrated Wasteland Development Programme (IWDP) implemented during 1997–2001 in Katangidda Nala watershed, Chincholi taluk, Gulbarga district, Karnataka. The study was carried out using IRS 1C, LISS III data of December 11, 1997 (pre-treatment) and November 15, 2002 (post-treatment) covering the watershed to assess the changes in land use / land cover and biomass that have changed over a period of five years (1997–2002). The images were classified into different land use/land cover categories using supervised classification by maximum likelihood algorithm. They were also classified into different biomass levels using Normalized Difference Vegetation Index (NDVI) approach. The results indicated that the area under agriculture crops and forest land were increased by 671 ha (5.7%) and 1,414 ha (11.94%) respectively. This is due to the fact that parts of wastelands and fallow lands were brought into cultivation. This increase in the area may be attributed to better utilization of surface and ground waters, adoption of soil and water conservation practices and changes in cropping pattern. The area under waste lands and fallow lands decreased by 1,667 ha (14.07%) and 467 ha (3.94%), respectively. The vegetation vigour of the area was classified into three classes using NDVI. Substantial increase in the area under high and low biomass levels was observed (502 ha and 19 ha respectively). The benefit-cost analysis indicates that the use of remote sensing and GIS was 2.2 times cheaper than the conventional methods. Thus, the repetitive coverage of the satellite data provides an excellent opportunity to monitor the land resources and evaluate the land cover changes through comparison of images for the watershed at different periods.     

Distribution analysis of salt affected soils under canal and non-canal command area in a part of Etah district,U. P., using remote sensing technique

A World Bank-aided project on sodic land reclamation in Uttar Pradesh is being executed by U.P. Bhumi Sudhar Nigam, Lucknow, and Remote Sensing Applications Centre, U.P., Lucknow has the responsibility of sodic land mapping for the execution of land reclamation programme at the cadastral level. Sodic lands are mainly concentrated in the Gangetic alluvial plains but the problem of sodicity is particularly acute in the canal-irrigated areas. A study of the distribution pattern of sodic lands in canal and noncanal command areas in a reclamation site (covering 60 villages out of which sodic lands were mapped in 51 villages) of Etah district in Uttar Pradesh, indicates that 18.39 per cent area of the canal command villages was barren sodic which was 13.41 per cent of the total geographical area of the site (15417 ha), however, 11.69 per cent area was recorded to be barren sodic in the non-canal command villages which was only 3.16 per cent of the geographical area of the site. The results of soil chemical analysis indicate that barren sodic lands of canal command area are saline-sodic with higher concentration of soluble salts (pH₂ >8.5, EC₂ >4 dSm⁻¹), however, those of non-canal command area are sodic (pH₂ >8.5, EC₂ <4 dSm⁻¹). The post-monsoon ground level in the canal-irrigated areas was in the critical and semicritical zone (< 3.0 mbgl) whereas it was well below the semi-critical zone in the non-canal command area, which indicates that the high ground water level is a major factor to higher the area under sodicity.     

Urban cadastral mapping using very high resolution remote sensing data

With growing urban expanses, one of the pre-requisites for effective governance is Urban Information Systems (UIS) with content down to individual properties (and individuals). The basic input i.e., a map, in UIS should show individual property boundaries showing the plan outline of all structures existing within, at a scale of 1:1000 and larger with sub-metre to centimeters planimetric and geometric accuracy. With very high resolution remote sensing data of the order of 1m available in hand, it is possible to prepare maps with high resolution spatial content. The present exercise demonstrates a method of preparing a geometrically and planimetrically accurate urban cadastral map on very large scale for a small area of about 5 sq km. IKONOS merged data with 1m resolution is used for the purpose. Mapping was done in conjunction with on-site measurements and sketches. Guides are used to maintain shape symmetry and accuracy of buildings and other features. Working out cost of mapping per unit area is another objective in the present exercise. For want of fully or semi-automatic methods of information extraction from very high resolution remote sensing data, it is imperative that mapping should be carried out in conjunction with some on-site measurements wherever necessary.     

Generation and validation of DEM using SAR interferometry and differential GPS supported by multispectral optical data

Digital elevation model (DEM) and the derived terrain parameters e.g. contour, slope, aspect, drainage pattern, etc are required for natural resources management, infrastructure planning and disaster management. The present paper aims at generating DEM from ERS tandem pair using interferometric technique supported by differential GPS measurements (DGPS) and multispectral optical data. Validation of DEM has been carried out by DGPS measurements. Ground Control Points (GCP) established by DGPS measurements have been used to georeference the IRS-1D optical data that has finally been co-registered with SAR amplitude image. Optical data, co-registered with ERS - I SAR data has helped in locating the GCP’s and check points, precisely, for refinement of DEM and its validation.     

Mapping of suspended sediments using site specific seasonal algorithms

This article reports a preliminary work in which two site specific seasonal algorithms have been proposed for estimating the suspended sediments concentration (SSC) from the digital numbers recorded on Indian Remote sensing Satellite, IRS-P4 Ocean Colour Monitor (OCM) sensor. For estimation of SSC, the proposed algorithms utilize dark pixel deduction atmospheric correction technique. The computations are performed with respect to north east monsoon phase situations of Palk Strait coastal stretch. The algorithms performance was satisfactory during the north east monsoon period. Although the results obtained cannot be generalized, we suggest that the authority of proposed algorithms can be extended to other seasons with the addition of more temporal experimental validation data sets and with numeric constants adjusted to present existing conditions. (As this area was severely affected by Tsunami, it may have dissimilar conditions at present).     

A study on the polarimetric properties of various features using SIR-C data

Spaceborne Imaging Radar (SIR-C) data acquired over Gujarat, India in 1994 were processed and analysed using differnet techniques applicable to polarimetric SAR data such as polarization signatures, polarization index, decomposition of the signal and polarization phase difference and limited groundtruth data. It has been observed that multi-frequency polarimetric data enhances the potential of retrieving geo-physical parameters. The polarization signatures are found to vary with the nature of the target. Target decomposition of the returned signal will be useful for the classification of various features. Polarization Phase Difference (PPD) gives good information about the vegetation parameters.     

The GeoForschungsZentrum Potsdam/Groupe de Recherche de Gèodésie Spatiale satellite-only and combined gravity field models: EIGEN-GL04S1 and EIGEN-GL04C

The recent improvements in the Gravity Recovery And Climate Experiment (GRACE) tracking data processing at GeoForschungsZentrum Potsdam (GFZ) and Groupe de Recherche de Géodésie Spatiale (GRGS) Toulouse, the availability of newer surface gravity data sets in the Arctic, Antarctica and North-America, and the availability of a new mean sea surface height model from altimetry processing at GFZ gave rise to the generation of two new global gravity field models. The first, EIGEN-GL04S1, a satellite-only model complete to degree and order 150 in terms of spherical harmonics, was derived by combination of the latest GFZ Potsdam GRACE-only (EIGEN-GRACE04S) and GRGS Toulouse GRACE/LAGEOS (EIGEN-GL04S) mean field solutions. The second, EIGEN-GL04S1 was combined with surface gravity data from altimetry over the oceans and gravimetry over the continents to derive a new high-resolution global gravity field model called EIGEN-GL04C. This model is complete to degree and order 360 and thus resolves geoid and gravity anomalies at half- wavelengths of 55 km at the equator. A degree-dependent combination method has been applied in order to preserve the high accuracy from the GRACE satellite data in the lower frequency band of the geopotential and to form a smooth transition to the high-frequency information coming from the surface data. Compared to pre-CHAMP global high-resolution models, the accuracy was improved at a spatial resolution of 200 km (half-wavelength) by one order of magnitude to 3 cm in terms of geoid heights. The accuracy of this model (i.e. the commission error) at its full spatial resolution is estimated to be 15 cm. The model shows a reduced artificial meridional striping and an increased correlation of EIGEN-GL04C-derived geostrophic meridional currents with World Ocean Atlas 2001 (WOA01) data. These improvements have led to select EIGEN-GL04C for JASON-1 satellite altimeter data reprocessing. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Methodology and use of tensor invariants for satellite gravity gradiometry

Although its use is widespread in several other scientific disciplines, the theory of tensor invariants is only marginally adopted in gravity field modeling. We aim to close this gap by developing and applying the invariants approach for geopotential recovery. Gravitational tensor invariants are deduced from products of second-order derivatives of the gravitational potential. The benefit of the method presented arises from its independence of the gradiometer instrument’s orientation in space. Thus, we refrain from the classical methods for satellite gravity gradiometry analysis, i.e., in terms of individual gravity gradients, in favor of the alternative invariants approach. The invariants approach requires a tailored processing strategy. Firstly, the non-linear functionals with regard to the potential series expansion in spherical harmonics necessitates the linearization and iterative solution of the resulting least-squares problem. From the computational point of view, efficient linearization by means of perturbation theory has been adopted. It only requires the computation of reference gravity gradients. Secondly, the deduced pseudo-observations are composed of all the gravitational tensor elements, all of which require a comparable level of accuracy. Additionally, implementation of the invariants method for large data sets is a challenging task. We show the fundamentals of tensor invariants theory adapted to satellite gradiometry. With regard to the GOCE (Gravity field and steady-state Ocean Circulation Explorer) satellite gradiometry mission, we demonstrate that the iterative parameter estimation process converges within only two iterations. Additionally, for the GOCE configuration, we show the invariants approach to be insensitive to the synthesis of unobserved gravity gradients.