Soil moisture estimation over vegetated terrains using multitemporal remote sensing data

A new method for retrieving soil moisture content over vegetated fields, employing multitemporal radar and optical images, is presented. It is based on the integration of the temporal series of radar data within an inversion scheme and on the correction of the vegetation effects. The retrieval algorithm uses the Bayesian maximum posterior probability and assumes the existence of a relation among the soil conditions at the different times of the series. The correction of the vegetation effects models the variation, with respect to the initial time of the series, of the component of the backscattering coefficient due to the soil characteristics as function of the variations of the measured backscattering coefficient and of the biomass. The method is tested on the data acquired throughout the SMEX02 experiment. The results show that measured and estimated soil moistures are fairly well correlated and that the performances of multitemporal retrieval algorithm are better than those obtained by employing one radar acquisition, especially in terms of capability to detect soil moisture changes. Although the approach to correct the vegetation effects on radar observations needs to be further assessed on different sets of data, this finding demonstrates that the proposed method has a potential to improve the quality of the soil moisture retrievals.     

Satellite radar and optical remote sensing for earthquake damage detection: results from different case studies

In case of a seismic event, a fast and draft damage map of the hit urban areas can be very useful, in particular when the epicentre of the earthquake is located in remote regions, or the main communication systems are damaged. Our aim is to analyse the capability of remote sensing techniques for damage detection in urban areas and to explore the combined use of radar (SAR) and optical satellite data. Two case studies have been proposed: Izmit (1999; Turkey) and Bam (2003; Iran). Both areas have been affected by strong earthquakes causing heavy and extended damage in the urban settlements close to the epicentre. Different procedures for damage assessment have been successfully tested, either to perform a pixel by pixel classification or to assess damage within homogeneous extended areas. We have compared change detection capabilities of different features extracted from optical and radar data, and analysed the potential of combining measurements at different frequency ranges. Regarding the Izmit case, SAR features alone have reached 70% of correct classification of damaged areas and 5 m panchromatic optical images have given 82%; the fusion of SAR and optical data raised up to 89% of correct pixel‐to‐pixel classification. The same procedures applied to the Bam test case achieved about 61% of correct classification from SAR alone, 70% from optical data, while data fusion reached 76%. The results of the correlation between satellite remote sensing and ground surveys data have been presented by comparing remotely change detection features averaged within homogeneous blocks of buildings with ground survey data.     

Uplift and subsidence due to the 26 December 2004 Indonesian earthquake detected by SAR data

The Indonesian earthquake took place on 26 December 2004 at 00:58 GMT (moment magnitude 9.3) in the Indian Ocean, offshore the west coast of Sumatra, at a depth of about 30 km. This earthquake is one of the largest of the past 100 years, comparable only with those in Chile (1960) and Alaska (1964). The earthquake originated in the subduction zone of the Indian and Burma plates, moving at a relative velocity of 6 cm/year. The aftershocks were distributed along a plate boundary of about 1000–1300 km between Sumatra and the Andaman Islands. Some hours after the earthquake a destructive tsunami followed and hit the coastlines of the surrounding regions, causing widespread destruction in Indonesia, India, Thailand and Sri Lanka. The European Space Agency (ESA) made available a data package composed of European Remote Sensing Satellite Synthetic Aperture Radar (ERS‐SAR) and Environment Satellite Advanced SAR (ENVISAT‐ASAR) data covering the affected area, acquired before (four acquisition dates) and after (five acquisition dates) the earthquake. A total of 26 frames were analysed. We used this dataset to evaluate the effects of the earthquake and tsunami on the human settlements and the physiographic conditions along the coast. The proposed method is based on a visual comparison between pre‐ and post‐seismic SAR intensity images, and on an analysis of their correlation coefficients. No complex data were made available by the ESA to exploit phase coherence. Analysis of pre‐ and post‐earthquake SAR backscattering showed wide uplift areas between the Andaman Islands and Simeulue Island, and large modifications of the coastline of Sumatra. Subsiding areas were detected along the southeast coast of Andaman up to the west coast of Nicobar Island. Tidal effects were filtered out of the SAR images to identify the consequences of the earthquake. Global Positioning System (GPS) stations in the Andaman provided results confirming the surface displacement pattern detected by SAR. The analysis enabled us to draw a boundary line separating the uplift and subsidence.     

Retrieval of atmospheric and surface parameters from satellite microwave radiometers over the Mediterranean Sea

A procedure to estimate atmospheric and sea surface parameters in the Mediterranean area from satellite microwave radiometric measurements is described. The method is founded on a simulator of brightness temperatures at the top of the atmosphere. The simulator is based on microwave sea emissivity and scattering model functions, derived from the outputs of the SEAWIND software, which implements a two-scale microwave sea surface model and a radiative transfer scheme in a nonscattering atmosphere. The development of the model functions aims to reduce the SEAWIND computational time, still maintaining its sensitivity to the main geophysical variables. Different adaptations of the simulation model have been performed to better reproduce the radiometric data in the region of interest. A comparison between the simulations and the Special Sensor Microwave/Imager (SSM/I) observations acquired throughout year 2000 over the Mediterranean Sea has permitted us to refine the model functions as well as to assess the whole simulation procedure. As for the inversion problem, a regression analysis has been applied to two different synthetic datasets to retrieve integrated precipitable water vapor, liquid water path and wind speed. The first dataset simulates the observations of SSM/I, whilst the second one concerns the Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E). Both have been generated by using the ECMWF atmospheric profiles and the measurements of the SeaWinds scatterometer aboard QuickSCAT. The SSM/I data have been used to carry out a statistical validation of the estimators. AMSR-E observations of a Tramontane-Mistral event, typical of the Mediterranean Sea, have been analyzed to evaluate the benefits of its expanded channel capability.     

Mapping sea surface temperature from aircraft using a multi-angle technique: An experiment over the Orbetello lagoon

Airborne scanner data collected in the thermal band over the Orbetello lagoon on 14 May 1993 were used to draw out a water surface temperature map. The geometric correction method was chosen after a quantitative comparison between the classic polynomial method and a technique based on matching Delaunay triangles. Sensor raw data were calibrated and striping noise was removed. Finally, a new atmospheric correction method was developed. The method consists of extracting, from the superimposed portion of two images, couples of measurements relative to the same resolution cell observed under two different angles. From these measurements we retrieved the atmospheric transmittance and integrated emitted radiance using an iterative algorithm capable of reducing the effects of measurement errors and taking into account the nonlinearity of the derived equations. The variation of surface emissivity with observation angle is also taken into account in the equations. In order to apply the method, the unknown covariance matrix of the parameters to be retrieved was valued. To this aim, we generated many statistical realizations of such atmospheric radiative quantities using the Lowtran 7 code applied to a large number of radiosoundings. The obtained atmospheric parameters were then used to correct the entire set of collected images and to produce a temperature map of the lagoon.     

Intercomparison of inversion algorithms to retrieve rain rate from SSM/I by using an extended validation set over the Mediterranean area

The capability of some inversion algorithms to estimate surface rain rate at the midlatitude basin scale from the Special Sensor Microwave Imager (SSM/I) data is analyzed. For this purpose, an extended database has been derived from coincident SSM/I images and half-hourly rain rate data obtained from a rain gauge network, placed along the Tiber River basin in Central Italy, during nine years (from 1992 to 2000). The database has been divided in a training set, to calibrate the empirical algorithms, and in a validation one, to compare the results of the considered techniques. The proposed retrieval methods are based on both empirical and physical approaches. Among the empirical methods, a regression, an artificial feedforward neural network, and a Bayesian maximum a posteriori (MAP) inversion have been considered. Three algorithms available in the literature are also included as benchmarks. As physical algorithms, the MAP method and the minimum mean square estimator have been used. Moreover, in order to test the behavior of the algorithms with different kinds of precipitation, a classification of rainy events, based on some statistical parameters derived from rain gauge measurements, has been performed. From this classification, an attempt to identify the type of event from radiometric data has been carried out. The purposes of this paper are to determine whether the use of an extended training set, referred to a limited geographical area, can improve the SSM/I skill in rain detection and estimation and, mainly, to confirm the validity of the physical approach adopted in previous works. It will be shown that, among all the estimators, the neural network presents the best performances and that the physical techniques provide results only slightly worse than those given by empirical methods, but with the well-known advantage of an easy application to different geographical zones and different sensors.     

Precipitation retrieval from spaceborne microwave radiometers basedon maximum a a posteriori probability estimation

A retrieval technique for estimating rainfall rate and precipitating cloud parameters from spaceborne multifrequency microwave radiometers is described. The algorithm is based on the maximum a posteriori probability criterion (MAP) applied to a simulated data base of cloud structures and related upward brightness temperatures. The cloud data base is randomly generated by imposing the mean values, the variances, and the correlations among the hydrometeor contents at each layer of the cloud vertical structure, derived from the outputs of a time-dependent microphysical cloud model. The simulated upward brightness temperatures are computed by applying a plane-parallel radiative transfer scheme. Given a multifrequency brightness temperature measurement, the MAP criterion is used to select the most probable cloud structure within the cloud-radiation data base. The algorithm is computationally efficient and has been numerically tested and compared against other methods. Its potential to retrieve rainfall over land has been explored by means of Special Sensor Microwave/Imager measurements for a rainfall event over Central Italy. The comparison of estimated rain rates with available raingauge measurements is also shown     

A physical-statistical approach to match passive microwave retrieval of rainfall to Mediterranean climatology

A physical-statistical approach to simulate cloud structures and their upward radiation over the Mediterranean is described. It aims to construct a synthetic database of microwave passive observations matching the climatological conditions of this geographical region. The synthetic database is conceived to train a Bayesian maximum a posteriori probability inversion scheme to retrieve precipitating cloud parameters from spaceborne microwave radiometric data. The initial microphysical a priori information on vertical profiles of cloud parameters is derived from a mesoscale cloud-resolving model. In order to complement information from cloud models and to match simulations to the conditions of the area of interest, a new approach is proposed. Climatological constraints over the Mediterranean are derived on a monthly basis from available radiosounding profiles, rain-gauge network measurements, and colocated METEOSAT infrared measurements. In order to introduce the actual surface background in the radiative-transfer simulations, a further constraint is represented by the monthly average and variance maps of surface emissivity derived from Special Sensor Microwave Imager (SSM/I) clear-air observations. A validation of the forward model is carried out by comparing a large set of brightness temperatures measured by the SSM/I with the synthetic cloud radiative database to asses its representativeness and range of variability.     

Sensitivity of bistatic scattering to soil moisture and surface roughness of bare soils

The sensitivity of bistatic scattering coefficient σ° to soil moisture content (SMC) and surface roughness was investigated by means of model simulations of the incoherent scattered fields performed with the advanced integral equation model (AIEM) and the second order small perturbation model (SPM). The study was performed by simulating scattering on the whole upper half space, for different values of incident angles. The achieved results, represented as maps of σ° as a function of azimuth and zenith angles, were evaluated by means of a quality index which takes into consideration the effect of roughness on SMC measurement. The sensitivity analysis has pointed out that for measuring SMC a bistatic observation, by itself or combined with the monostatic one, can make appreciable improvements with respect to classical monostatic radar. Appendix A contains the AIEM formulas corrected for several typographical errors present in the specific literature.     

Retrieving atmospheric temperature profiles by microwave radiometryusing a priori information on atmospheric spatial-temporalevolution

A new approach is presented to determine atmospheric temperature profiles by combining measurements coming from different sources and taking into account evolution models derived by conventional meteorological observations. Using a historical database of atmospheric parameters and related microwave brightness temperatures, the authors have developed a data assimilation procedure based on the geostatistical Kriging method and the Kalman filtering suitable for processing satellite radiometric measurements available at each satellite pass, data of a ground-based radiometer, and temperature profiles from radiosondes released at specific times and locations. The Kalman filter technique and the geostatistical Kriging method as well as the principal component analysis have proved very powerful in exploiting climatological a priori information to build spatial and temporal evolution models of the atmospheric temperature field. The use of both historical radiosoundings (RAOBs) and a radiative transfer code allowed the estimation of the statistical parameters that appears in the models themselves (covariance and cross-covariance matrices, observation matrix, etc.). The authors have developed an algorithm, based on a Kalman filter supplemented with a Kriging geostatistical interpolator, that shows a significant improvement of accuracy in vertical profile estimations with respect to the results of a standard Kalman filter when applied to real satellite radiometric data