Predicting floodplain inundation: raster‐based modelling versus the finite‐element approach

We compare two approaches to modelling floodplain inundation: a raster‐based approach, which uses a relatively simple process representation, with channel flows being resolved separately from the floodplain using either a kinematic or diffusive wave approximation, and a finite‐element hydraulic model aiming to solve the full two‐dimensional shallow‐water equations. A flood event on a short (c. 4 km) reach of the upper River Thames in the UK is simulated, the models being validated against inundation extent as determined from satellite synthetic aperture radar (SAR) imagery. The unconstrained friction parameters are found through a calibration procedure, where a measure of fit between predicted and observed shorelines is maximized. The raster and finite‐element models offer similar levels of performance, both classifying approximately 84% of the model domain correctly, compared with 65% for a simple planar prediction of water surface elevation. Further discrimination between models is not possible given the errors in the validation data. The simple raster‐based model is shown to have considerable advantages in terms of producing a straightforward calibration process, and being robust with respect to channel specification. Copyright © 2001 John Wiley & Sons, Ltd.     

The use of remotely sensed land cover to derive floodplain friction coefficients for flood inundation modelling

Remotely sensed land cover was used to generate spatially‐distributed friction coefficients for use in a two‐dimensional model of flood inundation. Such models are at the forefront of research into the prediction of river flooding. Standard practice, however, is to use single (static) friction coefficients on both the channel and floodplain, which are varied in a calibration procedure to provide a “best fit” to a known inundation extent. Spatially‐distributed friction provides a physically grounded estimate of friction that does not require fitting to a known inundation extent, but which can be fitted if desired. Remote sensing offers the opportunity to map these friction coefficients relatively straightforwardly and for low cost. Inundation was predicted using the LISFLOOD‐FP model for a reach on the River Nene, UK. Friction coefficients were produced from land cover predicted from Landsat TM imagery using both ML and fuzzy c‐means classifiction. The elevetion data used were from combined contour and differential global positioning system (GPS) elevation data. Predicted inundation using spatially‐distributed and static friction were compared. Spatially‐distributed friction had the greatest effect on the timing of flood inundation, but a small effect on predicted inundation extent. The results indicate that spatially‐distributed friction should be considered where the timing of initial flooding (e.g. for early warning) is important. Copyright © 2007 John Wiley & Sons, Ltd.     

Modelling floodplain flows using a two-dimensional finite element model

A prototype two-dimensional finite element flow model for depth-averaged free surface flows was developed for floodplain environments. Limited refinement of the model's physical representation was undertaken and the enhanced scheme applied to an 11 km river channel/floodplain reach in the U.K. Preliminary model results indicate that this modelling approach can be used to identify dynamic variations in the flow field parameters over length scales of the order of 10-100 m. Potentially, such data have the ability to permit detailed analysis of short-term floodplain sedimentary dynamics.     

Characterisation of flood inundated areas and delineation of poor drainage soil using ERS-1 SAR imagery

Recent years have been marked by a continuous availability of spatial SAR data since the launch of the European remote sensing satellite (ERS-1) in 1991. Consequently, remote sensing techniques now offer an opportunity to map flood inundation fields caused by river overflow or waterlogging in environments characterized by frequent cloud cover. Indeed, inundation fields can clearly be seen on ERS-1 SAR images taken during flooding periods. However, such an identification can be constrained by the similarity in behaviour between water surfaces and other features of the landscape such as extended asphalt areas, permanent water bodies and less illuminated slopes. For consistent flood inundation extent mapping a more robust approach is required. This is provided by a conceptual flood inundation index that is physically sound in relation to radar imaging. Moreover, this index has proved to be useful for highlighting soils located within inundation fields and having significantly different internal drainage. The results achieved in the framework of the research must be seen in the context of intensive use of remote sensing data to support decision methods for sustainable management of land and water resources. Such decision support methods could be provided by river hydraulic models aimed at assessing environmental effects of inundation floods and at early flood warning systems. © 1997 John Wiley & Sons, Ltd.     

Remote sensing of river stage using the cross‐sectional inundation area‐river stage relationship (IARSR) constructed from digital elevation model data

Remote sensing of discharge and river stage from space provides us with a promising alternative approach to monitor watersheds, no matter if they are ungauged, poorly gauged, or fully gauged. One approach is to estimate river stage from satellite measured inundation area based on the inundation area – river stage relationship (IARSR). However, this approach is not easy to implement because of a lack of data for constructing the IARSR. In this study, an innovative and robust approach to construct the IARSR from digital elevation model (DEM) data was developed and tested. It was shown that the constructed IARSR from DEM data could be used to retrieve water level or river stage from satellite‐measured inundation area. To reduce the uncertainty in the estimated inundation area, a dual‐thresholding method was proposed. The first threshold is the lower limit of pixel value for classifying water body pixels with a relatively high‐level certainty. The second threshold is the upper limit of pixel value for classifying potentially flooded pixels. All pixels with values between the first threshold and the second threshold and adjacent to the classified water body pixels may be partially flooded. A linear interpolation method was used to estimate the wetted area of each partially flooded pixel. In applying the constructed IARSR to the estimated inundation areas from 11 Landsat TM images, 11 water levels were obtained. The root mean square error (RMSE) of the estimated water levels compared with the observed water levels at the US Geological Survey (USGS) gauging station on the Trinity River at Liberty in Liberty County, Texas, is about 0.38 m. Copyright © 2012 John Wiley & Sons, Ltd.     

Estimation of uncertainty propagation in flood inundation mapping using a 1‐D hydraulic model

A need for more accurate flood inundation maps has recently arisen because of the increasing frequency and extremity of flood events. The accuracy of flood inundation maps is determined by the uncertainty propagated from all of the variables involved in the overall process of flood inundation modelling. Despite our advanced understanding of flood progression, it is impossible to eliminate the uncertainty because of the constraints involving cost, time, knowledge, and technology. Nevertheless, uncertainty analysis in flood inundation mapping can provide useful information for flood risk management. The twin objectives of this study were firstly to estimate the propagated uncertainty rates of key variables in flood inundation mapping by using the first‐order approximation method and secondly to evaluate the relative sensitivities of the model variables by using the Hornberger–Spear–Young (HSY) method. Monte Carlo simulations using the Hydrologic Engineering Center's River Analysis System and triangle‐based interpolation were performed to investigate the uncertainty arising from discharge, topography, and Manning's n in the East Fork of the White River near Seymour, Indiana, and in Strouds Creek in Orange County, North Carolina. We found that the uncertainty of a single variable is propagated differently to the flood inundation area depending on the effects of other variables in the overall process. The uncertainty was linearly/nonlinearly propagated corresponding to valley shapes of the reaches. In addition, the HSY sensitivity analysis revealed the topography of Seymour reach and the discharge of Strouds Creek to be major contributors to the change of flood inundation area. Copyright © 2014 John Wiley & Sons, Ltd.     

A comparison of continuous mass‐lumped finite elements with finite differences for 3‐D wave propagation

The finite‐difference method on rectangular meshes is widely used for time‐domain modelling of the wave equation. It is relatively easy to implement high‐order spatial discretization schemes and parallelization. Also, the method is computationally efficient. However, the use of finite elements on tetrahedral unstructured meshes is more accurate in complex geometries near sharp interfaces. We compared the standard eighth‐order finite‐difference method to fourth‐order continuous mass‐lumped finite elements in terms of accuracy and computational cost. The results show that, for simple models like a cube with constant density and velocity, the finite‐difference method outperforms the finite‐element method by at least an order of magnitude. Outside the application area of rectangular meshes, i.e., for a model with interior complexity and topography well described by tetrahedra, however, finite‐element methods are about two orders of magnitude faster than finite‐difference methods, for a given accuracy.     

Relative role and extent of marine and groundwater inundation on a dune‐dominated barrier island under sea‐level rise scenarios

Climate change and sea‐level rise will have severe impacts on coastal water resources around the world. However, whereas the influence of marine inundation is well documented in the literature, the impact of groundwater inundation on coastal communities is not well known. Here, core analysis, groundwater monitoring, and ground penetrating radar are utilized to assess the groundwater regime of the surficial aquifer on Bogue Banks Barrier Island (USA). Then, geospatial techniques are used to assess the relative roles and extents of groundwater and marine inundation on the dune‐dominated barrier island under sea‐level rise scenarios of 0.2, 0.5, and 1.0 m above current conditions by 2100. Additionally, the effects of rising water tables on onsite wastewater treatment systems (OWTS) are modelled using the projected sea‐level rise scenarios. The results indicate that the surficial aquifer comprising fine to medium sands responds quickly to precipitation. Water‐level measurements reveal varying thicknesses of the vadose zone (>3 to 0 m) and several groundwater mounds with radial flow patterns. Results from projected sea‐level rise scenarios suggest that owing to aquifer properties and morphology of the island, groundwater inundation may occur at the same rate as marine inundation. Furthermore, the area inundated by groundwater may be as significant as that affected by marine inundation. The results also show that the proportion of land in the study area where OWTS may be perpetually compromised by rising water tables under worst case scenarios may range from ~43 to ~54% over an 86‐year‐period. Copyright © 2014 John Wiley & Sons, Ltd.     

Comparing 2D capabilities of HEC-RAS and LISFLOOD-FP on complex topography

ABSTRACT

This study evaluates and compares two-dimensional (2D) numerical models of different complexities by testing them on a floodplain inundation event that occurred on the Secchia River (Italy). We test 2D capabilities of LISFLOOD-FP and HEC-RAS (5.0.3), implemented using various grid sizes (25–100 m) based on 1-m DEM resolution. As expected, the best results were shown by the higher-resolution grids (25 m) for both models, which is justified by the complex terrain of the area. However, the coarser resolution simulations (50 and 100 m) performed virtually identically compared to the high-resolution simulations. Nevertheless, the spatial distribution of flood characteristics varies: the 50 and 100 m results of LISFLOOD-FP and HEC-RAS misestimated flood extent and water depth in selected control areas (built-up zones). We suggest that the specific terrain of the area can cause ambiguities in large-scale modelling, while providing plausible results in terms of the overall model performance.     

A 2D finite volume model for bebris flow and its application to events occurred in the Eastern Pyrenees

FLATModel is a 2D finite volume code that contains several original approaches to improve debris-flow simulation. Firstly, FLATModel incorporates a ＂stop-and-go＂ technique in each cell to allow continuous collapses and remobilizations of the debris-flow mass. Secondly, flow velocity and consequently yield stress is directly associated with the type of rheology to improve boundary accuracy. Thirdly, a simple approach for entrainment is also included in the model to analyse the effect of basal erosion of debris flows. FLATMODEL was tested at several events that occurred in the Eastern Pyrenees and simulation results indicated that the model can represent rather well the different characteristics observed in the field.     

Satellite remote sensing of river inundation area,stage, and discharge: a review

The growing availability of multi-temporal satellite data has increased opportunities for monitoring large rivers from space. A variety of passive and active sensors operating in the visible and microwave range are currently operating, or planned, which can estimate inundation area and delineate flood boundaries. Radar altimeters show great promise for directly measuring stage variation in large rivers. It also appears to be possible to obtain estimates of river discharge from space, using ground measurements and satellite data to construct empirical curves that relate water surface area to discharge. Extrapolation of these curves to ungauged sites may be possible for the special case of braided rivers. Where clouds, trees and floating vegetation do not obscure the water surface, high-resolution visible/infrared sensors provide good delineation of inundated areas. Synthetic aperture radar (SAR) sensors can penetrate clouds and can also detect standing water through emergent aquatic plants and forest canopies. However, multiple frequencies and polarizations are required for optimal discrimination of various inundated vegetation cover types. Existing single-polarization, fixed-frequency SARs are not sufficient for mapping inundation area in all riverine environments. In the absence of a space-borne multi-parameter SAR, a synergistic approach using single-frequency, fixed-polarization SAR and visible/infrared data will provide the best results over densely vegetated river floodplains. © 1997 John Wiley & Sons, Ltd.     

城镇建筑震害SAR遥感探测与评估研究综述

合成孔径雷达(SAR)遥感是城镇建筑震害监测的重要手段之一.随着SAR传感器技术的最新发展,利用高分辨率、多极化SAR图像对城镇建筑地震灾害损毁进行探测和评估成为当前的研究趋势和热点.文中首先分析了建筑物及其震害损毁SAR成像特点,并系统综述了近15年来国内外基于SAR图像进行城镇建筑震害损毁探测与评估技术方法的研究和应用现状,评述了各类方法的优缺点,最后针对当前SAR传感器高分辨率、全极化的发展趋势,对有待于进一步研究的问题和技术发展趋势进行了展望.     

Correlation of dynamic characteristics of a super‐tall building from full‐scale measurements and numerical analysis with various finite element models

The Di Wang Tower located in Shenzhen has 79 storeys and is about 325 m high. Field measurements have been conducted to investigate the dynamic characteristics of the super‐tall building. In parallel with the field measurements, seven finite element models have been established to model the multi‐outrigger‐braced tall building and to analyse the effects of various arrangements of outrigger belts and vertical bracings on the dynamic characteristics and responses of the Di Wang Tower under the design wind load and earthquake action. The distributions of shear forces in vertical structural components along the building height are also presented and discussed. The results from detailed modelling of group shear walls with several types of finite elements are addressed and compared to investigate various modelling assumptions. Finally, the performance of the finite element models is evaluated by correlating the natural frequencies and mode shapes from the numerical analysis with the finite element models and the field measurements. The results generated from this study are expected to be of interest to professionals and researchers involved with the design of tall buildings. Copyright © 2004 John Wiley & Sons, Ltd.     

基于微增模型的海洋可控源电磁法三维非结构化矢量有限元数值模拟

为了分析海洋中含金属矿的储油气砂层的成分、结构对于海洋可控源电磁场的影响特征,引入多相微增等效介质模型,结合非结构化网格实现了海洋可控源电磁三维矢量有限元正演.首先对多相微增模型进行了介绍,并分析不同参数影响下的等效电导率变化特征.而后,详细推导了频率域可控源电磁法非结构化矢量有限元方程.接着,利用预处理的IDR（s）迭代算法求解线性方程组.最后设计典型的海洋地电模型,通过正演计算,验证了本文算法的正确性和有效性,同时研究了海洋储油砂层中的金属矿的含量、孔隙度、含水饱和度等参数对于可控源电磁场响应特征的影响.

     

A comparison of one‐ and two‐dimensional approaches to modelling flood inundation over complex upland floodplains

A much understudied aspect of flood inundation is examined, i.e. upland environments with topographically complex floodplains. Although the presence of high‐resolution topographic data (e.g. lidar) has improved the quality of river flood inundation predictions, the optimum dimensionality of hydraulic models for this purpose has yet to be fully evaluated for situations of both topographic and topological (i.e. the connectivity of floodplain features) complexity. In this paper, we present the comparison of three treatments of upland flood inundation using: (a) a one‐dimensional (1D) model (HEC‐RAS v. 3·1·2) with the domain defined as series of extended cross‐sections; (b) the same 1D model, but with the floodplain defined by a series of storage cells, hydraulically connected to the main river channel and other storage cells on the floodplain according to floodplain topological characteristics; (c) a two‐dimensional (2D) diffusion wave treatment, again with explicit representation of floodplain structural features. The necessary topographic and topological data were derived using lidar and Ordnance Survey Landline data. The three models were tested on a 6 km upland reach of the River Wharfe, UK. The models were assessed by comparison with measured inundation extent. The results showed that both the extended cross‐section and the storage cell 1D modes were conceptually problematic. They also resulted in poorer model predictions, requiring incorrect parameterization of the main river to floodplain flux in order to approach anything like the level of agreement observed when the 2D diffusion wave treatment was assessed. We conclude that a coupled 1D–2D treatment is likely to provide the best modelling approach, with currently available technology, for complex floodplain configurations. Copyright © 2007 John Wiley & Sons, Ltd.     

Perspectives in using a remotely sensed dryness index in distributed hydrological models at the river‐basin scale

In a previous study a spatially distributed hydrological model, based on the MIKE SHE code, was constructed and validated for the 375 000 km² Senegal River basin in West Africa. The model was constructed using spatial data on topography, soil types and vegetation characteristics together with time‐series of precipitation from 112 stations in the basin. The model was calibrated and validated based on river discharge data from nine stations in the basin for 11 years. Calibration and validation results suggested that the spatial resolution of the input data in parts of the area was not sufficient for a satisfactory evaluation of the modelling performance. The study further examined the spatial patterns in the model input and output, and it was found that particularly the spatial resolution of the precipitation input had a major impact on the model response. In an attempt to improve the model performance, this study examines a remotely sensed dryness index for its relationship to simulated soil moisture and evaporation for six days in the wet season 1990. The index is derived from observations of surface temperature and vegetation index as measured by the NOAA Advanced Very High Resolution Radiometer (AVHRR) sensor. The correlation results between the index and the simulation results are of mixed quality. A sensitivity analysis, conducted on both estimates, reveals significant uncertainties in both. The study suggests that the remotely sensed dryness index with its current use of NOAA AVHRR data does not offer information that leads to a better calibration or validation of the simulation model in a spatial sense. The method potentially may become more suitable with the use of the upcoming high‐resolution temporal Meteosat Second Generation data. Copyright © 2002 John Wiley & Sons, Ltd.     

Three‐dimensional controlled‐source electromagnetic modelling with a well casing as a grounded source: a hybrid method of moments and finite element scheme

Steel well casings in or near a hydrocarbon reservoir can be used as source electrodes in time‐lapse monitoring using grounded line electromagnetic methods. A requisite component of carrying out such monitoring is the capability to numerically model the electromagnetic response of a set of source electrodes of finite length. We present a modelling algorithm using the finite‐element method for calculating the electromagnetic response of a three‐dimensional conductivity model excited using a vertical steel‐cased borehole as a source. The method is based on a combination of the method of moments and the Coulomb‐gauged primary–secondary potential formulation. Using the method of moments, we obtain the primary field in a half‐space due to an energized vertical steel casing by dividing the casing into a set of segments, each assumed to carry a piecewise constant alternating current density. The primary field is then substituted into the primary–secondary potential finite‐element formulation of the three‐dimensional problem to obtain the secondary field. To validate the algorithm, we compare our numerical results with: (i) the analytical solution for an infinite length casing in a whole space, excited by a line source, and (ii) a three‐layered Earth model without a casing. The agreement between the numerical and analytical solutions demonstrates the effectiveness of our algorithm. As an illustration, we also present the time‐lapse electromagnetic response of a synthetic model representing a gas reservoir undergoing water flooding.     

Delineation of riparian vegetation from Landsat multi‐temporal imagery using PCA

A deficiency in crucial digital data, such as vegetation cover, in remote regions is a challenging issue for water management and planning, especially for areas undergoing rapid development, such as mining in the Pilbara, Western Australia. This is particularly relevant to riparian vegetation, which provides important ecological services and, as such, requires regional protection. The objective of this research was to develop an approach to riparian vegetation mapping at a regional scale using remotely sensed data. The proposed method was based on principal component analysis applied to multi‐temporal Normalized Difference Vegetation Index datasets derived from Landsat TM 5 imagery. To delimit the spatial extent of riparian vegetation, a thresholding method was required and various thresholding algorithms were tested. The accuracy of results was estimated for various Normalized Difference Vegetation Index multi‐temporal datasets using available ground‐truth data. The combination of a 14‐dry‐date dataset and Kittler's thresholding method provided the most accurate delineation of riparian vegetation.