Laura Giustarini

A Change Detection Approach to Flood Mapping in Urban Areas Using TerraSAR-X

Very high resolution synthetic aperture radar (SAR) sensors represent an alternative to aerial photography for delineating floods in built-up environments where flood risk is highest. However, even with currently available SAR image resolutions of 3 m and higher, signal returns from man-made structures hamper the accurate mapping of flooded areas. Enhanced image processing algorithms and a better exploitation of image archives are required to facilitate the use of microwave remote-sensing data for monitoring flood dynamics in urban areas. In this paper, a hybrid methodology combining backscatter thresholding, region growing, and change detection (CD) is introduced as an approach enabling the automated, objective, and reliable flood extent extraction from very high resolution urban SAR images. The method is based on the calibration of a statistical distribution of “open water” backscatter values from images of floods. Images acquired during dry conditions enable the identification of areas that are not “visible” to the sensor (i.e., regions affected by “shadow”) and that systematically behave as specular reflectors (e.g., smooth tarmac, permanent water bodies). CD with respect to a reference image thereby reduces overdetection of inundated areas. A case study of the July 2007 Severn River flood (UK) observed by airborne photography and the very high resolution SAR sensor on board TerraSAR-X highlights advantages and limitations of the method. Even though the proposed fully automated SAR-based flood-mapping technique overcomes some limitations of previous methods, further technological and methodological improvements are necessary for SAR-based flood detection in urban areas to match the mapping capability of high-quality aerial photography.

Detection of flooded urban areas in high resolution Synthetic Aperture Radar images using double scattering

Abstract Flooding is a particular hazard in urban areas worldwide due to the increased risks to life and property in these regions. Synthetic Aperture Radar (SAR) sensors are often used to image flooding because of their all-weather day–night capability, and now possess sufficient resolution to image urban flooding. The flood extents extracted from the images may be used for flood relief management and improved urban flood inundation modelling. A difficulty with using SAR for urban flood detection is that, due to its side-looking nature, substantial areas of urban ground surface may not be visible to the SAR due to radar layover and shadow caused by buildings and taller vegetation. This paper investigates whether urban flooding can be detected in layover regions (where flooding may not normally be apparent) using double scattering between the (possibly flooded) ground surface and the walls of adjacent buildings. The method estimates double scattering strengths using a SAR image in conjunction with a high resolution LiDAR (Light Detection and Ranging) height map of the urban area. A SAR simulator is applied to the LiDAR data to generate maps of layover and shadow, and estimate the positions of double scattering curves in the SAR image. Observations of double scattering strengths were compared to the predictions from an electromagnetic scattering model, for both the case of a single image containing flooding, and a change detection case in which the flooded image was compared to an un-flooded image of the same area acquired with the same radar parameters. The method proved successful in detecting double scattering due to flooding in the single-image case, for which flooded double scattering curves were detected with 100% classification accuracy (albeit using a small sample set) and un-flooded curves with 91% classification accuracy. The same measures of success were achieved using change detection between flooded and un-flooded images. Depending on the particular flooding situation, the method could lead to improved detection of flooding in urban areas.

Probabilistic Flood Mapping Using Synthetic Aperture Radar Data

Probabilistic flood mapping offers flood managers, decision makers, insurance agencies, and humanitarian relief organizations a useful characterization of uncertainty in flood mapping delineation. Probabilistic flood maps are also of high interest for data assimilation into numerical models. The direct assimilation of probabilistic flood maps into hydrodynamic models would be beneficial because it would eliminate the intermediate step of having to extract water levels first. This paper introduces a probabilistic flood mapping procedure based on synthetic aperture radar (SAR) data. Given a SAR image of backscatter values, we construct a total histogram of backscatter values and decompose this histogram into probability distribution functions of backscatter values associated with flooded (open water) and non-flooded pixels, respectively. These distributions are then used to estimate, for each pixel, its probability of being flooded. The new approach improves on binary SAR-based flood mapping procedures, which do not inform on the uncertainty in the pixel state. The proposed approach is tested using four SAR images from two floodplains, i.e., the Severn River (U.K.) and the Red River (U.S.). In all four test cases, reliability diagrams, with error values ranging from 0.04 to 0.23, indicate a good agreement between the SAR-derived probabilistic flood map and an independently available validation map, which is obtained from aerial photography.

Flood Hazard Mapping Combining Hydrodynamic Modeling and Multi Annual Remote Sensing data

This paper explores a method to combine the time and space continuity of a large-scale inundation model with discontinuous satellite microwave observations, for high-resolution flood hazard mapping. The assumption behind this approach is that hydraulic variables computed from continuous spatially-distributed hydrodynamic modeling and observed as discrete satellite-derived flood extents are correlated in time, so that probabilities can be transferred from the model series to the observations. A prerequisite is, therefore, the existence of a significant correlation between a modeled variable (i.e., flood extent or volume) and the synchronously-observed flood extent. If this is the case, the availability of model simulations over a long time period allows for a robust estimate of non-exceedance probabilities that can be attributed to corresponding synchronously-available satellite observations. The generated flood hazard map has a spatial resolution equal to that of the satellite images, which is higher than that of currently available large scale inundation models. The method was applied on the Severn River (UK), using the outputs of a global inundation model provided by the European Centre for Medium-range Weather Forecasts

Enhanced biomass prediction by assimilating satellite data into a crop growth model

Complex crop growth models (CGM) require a large number of input parameters, which can cause large errors if they are uncertain. Furthermore, they often lack spatial information. The coupling of a CGM with a radiative transfer model offers the possibility to assimilate remote sensing data while taking into account uncertainties in input parameters. A particle filter was used to assimilate satellite data into a CGM coupled with a leaf-canopy radiative transfer model to update biomass simulations of maize. The synthetic experiment set up to test the reliability of the procedure, highlighted the importance of the acquisition time. The real case study with RapidEye observations confirmed these findings. Data assimilation increased the accuracy of biomass predictions in the majority of the six maize fields where biomass validation data was available, with improvements of up to 15%. The smallest and largest errors in biomass prediction after assimilation were 82?kg/ha and 2116?kg/ha, respectively. Furthermore, data assimilation enabled the production of biomass maps showing detailed spatial variability. Data assimilation using a particle filter for biomass estimation was conducted.Proof of concept with synthetic case studies.Multispectral satellite data (visible and near infrared) was found to be suitable for data assimilation.Assimilation of satellite data allowed biomass prediction on a pixel basis.

A Hierarchical Split-Based Approach for Parametric Thresholding of SAR Images: Flood Inundation as a Test Case

Parametric thresholding algorithms applied to synthetic aperture radar (SAR) imagery typically require the estimation of two distribution functions, i.e., one representing the target class and one its background. They are eventually used for selecting the threshold that allows binarizing the image in an optimal way. In this context, one of the main difficulties in parameterizing these functions originates from the fact that the target class often represents only a small fraction of the image. Under such circumstances, the histogram of the image values is often not obviously bimodal and it becomes difficult, if not impossible, to accurately parameterize distribution functions. Here we introduce a hierarchical split-based approach that searches for tiles of variable size allowing the parameterization of the distributions of two classes. The method is integrated into a flood-mapping algorithm in order to evaluate its capacity for parameterizing distribution functions attributed to floodwater and changes caused by floods. We analyzed a data set acquired during a flood event along the Severn River (U.K.) in 2007. It is composed of moderate (ENVISAT-WS) and high (TerraSAR-X)-resolution SAR images. The obtained classification accuracies as well as the similarity of performance levels to a benchmark obtained with an established method based on the manual selection of tiles indicate the validity of the new method.

A user-driven case-based reasoning tool for infilling missing values in daily mean river flow records

Missing data in river flow records represent a loss of information and a serious drawback in water management. In this work, we introduce gapIt, a user-driven case-based reasoning tool for infilling gaps in daily mean river flow records. Given a set of flow time series, gapIt builds a database of artificial gaps for which it computes several flow estimates, to find the best combinations of infilling algorithm and automatically selected donor station(s), according to state-of-the-art performance indicators. We obtained satisfactory results with Nash-Sutcliffe >0.7 for more than half of the ~5000 synthetic gaps of various lengths and positions, randomly created along the available records. gapIt was evaluated on 24 daily river discharge time series recorded in Luxembourg over seven years from 01/01/2007 to 31/12/2013. We also discuss the benefits of coupling this approach with user-expertise for an improved infilling of real data gaps. We propose a user-driven case-based reasoning tool to infill gaps in flow records.With data analytics the tool finds the best combination: algorithm - donor stations.Satisfactory results were obtained on synthetic gaps of different characteristics.The tool bridges a gap between data-driven and user-expertise approaches.

Compilation and Validation of SAR and Optical Data Products for a Complete and Global Map of Inland/Ocean Water Tailored to the Climate Modeling Community

Accurate maps of surface water extent are of paramount importance for water management, satellite data processing and climate modeling. Several maps of water bodies based on remote sensing data have been released during the last decade. Nonetheless, none has a truly (90 ∘ N/90 ∘ S) global coverage while being thoroughly validated. This paper describes a global, spatially-complete (void-free) and accurate mask of inland/ocean water for the 2000–2012 period, built in the framework of the European Space Agency (ESA) Climate Change Initiative (CCI). This map results from the synergistic combination of multiple individual SAR and optical water body and auxiliary datasets. A key aspect of this work is the original and rigorous stratified random sampling designed for the quality assessment of binary classifications where one class is marginally distributed. Input and consolidated products were assessed qualitatively and quantitatively against a reference validation database of 2110 samples spread throughout the globe. Using all samples, overall accuracy was always very high among all products, between 98 % and 100 % . The CCI global map of open water bodies provided the best water class representation (F-score of 89 % ) compared to its constitutive inputs. When focusing on the challenging areas for water bodies’ mapping, such as shorelines, lakes and river banks, all products yielded substantially lower accuracy figures with overall accuracies ranging between 74 % and 89 % . The inland water area of the CCI global map of open water bodies was estimated to be 3.17 million km 2 ± 0.24 million km 2 . The dataset is freely available through the ESA CCI Land Cover viewer.

From SAR-derived flood mapping to water level data assimilation into hydraulic models

This paper describes a fully automatic processing chain that makes use of SAR images for retrieving water stage information to be assimilated into a hydraulic forecasting model. This chain is composed of three steps: flood extent delineation, water stage retrieval and data assimilation of stage information into a hydraulic model. The flood-mapping step is addressed with a fully automatic algorithm, based on image statistics and applicable to all existing SAR datasets. Uncertainty on the flood extent map is represented with an ensemble of flood extent maps, obtained following a bootstrap methodology. Water stage observations are then retrieved by intersecting the flood shoreline with the floodplain topography. The ensemble of flood extent maps allows extracting multiple water levels at any river cross section of the hydraulic model, thereby taking into account the uncertainty associated with the floodmapping step. Finally, data assimilation consists in integrating uncertain observations, i.e. SAR-derived water stages, with uncertain hydraulic model simulations. The proposed processing chain was applied to two case studies. For the test case of June 2008 on the Po River (Italy), only low resolution but freely available satellite data were used. For the January 2011 flood on the Sure River (Luxembourg), higher resolution data were used and obtained at a cost. The results show that with the assimilation of SAR-derived water stages significant improvements can be achieved in the forecasting performance of the hydraulic model.

Flood hazard mapping combining high resolution multi-temporal SAR data and coarse resolution global hydrodynamic modelling

A new method for flood hazard mapping that integrates global flood inundation modeling and microwave remote sensing is presented. It combines the time and space continuity of a global inundation model with the limited revisit time but high spatial resolution of satellite observations. The availability of model simulations over a long time period allows a robust estimate of non-exceedance probabilities that can be attributed to the corresponding satellite observations. The resulting flood hazard map will have a spatial resolution equal to that of the used satellite images, generally higher than that of the global inundation model. This can theoretically be done for any point in the world, allowing the estimation of flood hazard at a global scale, provided that a sufficient number of remote sensing images are available. The method is tested on the Severn River (UK), with a high number of flood events observed by ENVISAT ASAR. The global ECMWF flood inundation model is considered for this study.