Calibration of uncertain flood inundation models using remotely sensed water levels

A traditional method of validating the performance of a flood model when remotely sensed data of the flood extent are available is to compare the predicted flood extent to that observed. The performance measure employed often uses areal pattern-matching to assess the degree to which the two extents overlap. Recently, remote sensing of flood extents using synthetic aperture radar (SAR) and airborne scanning laser altimetry (LIDAR) has made more straightforward the synoptic measurement of water surface elevations along flood waterlines, and this has emphasised the possibility of using alternative performance measures based on height. This paper considers the advantages that can accrue from using a performance measure based on waterline elevations rather than one based on areal patterns of wet and dry pixels. The two measures were compared for their ability to estimate flood inundation uncertainty maps from a set of model runs carried out to span the acceptable model parameter range in a GLUE-based analysis. A 1 in 5-year flood on the Thames in 1992 was used as a test event. As is typical for UK floods, only a single SAR image of observed flood extent was available for model calibration and validation. A simple implementation of a two-dimensional flood model (LISFLOOD-FP) was used to generate model flood extents for comparison with that observed. The performance measure based on height differences of corresponding points along the observed and modelled waterlines was found to be significantly more sensitive to the channel friction parameter than the measure based on areal patterns of flood extent. The former was able to restrict the parameter range of acceptable model runs and hence reduce the number of runs necessary to generate an inundation uncertainty map. A result of this was that there was less uncertainty in the final flood risk map. The uncertainty analysis included the effects of uncertainties in the observed flood extent as well as in model parameters. The height-based measure was found to be more sensitive when increased heighting accuracy was achieved by requiring that observed waterline heights varied slowly along the reach. The technique allows for the decomposition of the reach into sections, with different effective channel friction parameters used in different sections, which in this case resulted in lower r.m.s. height differences between observed and modelled waterlines than those achieved by runs using a single friction parameter for the whole reach. However, a validation of the modelled inundation uncertainty using the calibration event showed a significant difference between the uncertainty map and the observed flood extent. While this was true for both measures, the difference was especially significant for the height-based one. This is likely to be due to the conceptually simple flood inundation model and the coarse application resolution employed in this case. The increased sensitivity of the height-based measure may lead to an increased onus being placed on the model developer in the production of a valid model

Automatic near real-time selection of flood water levels from high resolution Synthetic Aperture Radar images for assimilation into hydraulic models: a case study

Flood extents caused by fluvial floods in urban and rural areas may be predicted by hydraulic models. Assimilation may be used to correct the model state and improve the estimates of the model parameters or external forcing. One common observation assimilated is the water level at various points along the modelled reach. Distributed water levels may be estimated indirectly along the flood extents in Synthetic Aperture Radar (SAR) images by intersecting the extents with the floodplain topography. It is necessary to select a subset of levels for assimilation because adjacent levels along the flood extent will be strongly correlated. A method for selecting such a subset automatically and in near real-time is described, which would allow the SAR water levels to be used in a forecasting model. The method first selects candidate waterline points in flooded rural areas having low slope. The waterline levels and positions are corrected for the effects of double reflections between the water surface and emergent vegetation at the flood edge. Waterline points are also selected in flooded urban areas away from radar shadow and layover caused by buildings, with levels similar to those in adjacent rural areas. The resulting points are thinned to reduce spatial autocorrelation using a top-down clustering approach. The method was developed using a TerraSAR-X image from a particular case study involving urban and rural flooding. The waterline points extracted proved to be spatially uncorrelated, with levels reasonably similar to those determined manually from aerial photographs, and in good agreement with those of nearby gauges.

Modelling the hydraulics of the Carlisle 2005 flood event

The performance of a 2D numerical model of flood hydraulics is tested for a major event in Carlisle, UK, in 2005. This event is associated with a unique data set, with GPS surveyed wrack lines and flood extent surveyed 3 weeks after the flood. The Simple Finite Volume (SFV) model is used to solve the 2D Saint-Venant equations over an unstructured mesh of 30000 elements representing channel and floodplain, and allowing detailed hydraulics of flow around bridge piers and other influential features to be represented. The SFV model is also used to corroborate flows recorded for the event at two gauging stations. Calibration of Manning's n is performed with a two stage strategy, with channel values determined by calibration of the gauging station models, and floodplain values determined by optimising the fit between model results and observed water levels and flood extent for the 2005 event. RMS error for the calibrated model compared with surveyed water levels is ~±0.4m, the same order of magnitude as the estimated error in the survey data. The study demonstrates the ability of unstructured mesh hydraulic models to represent important hydraulic processes across a range of scales, with potential applications to flood risk management.

Phosphorus dynamics in the ditch system of a restored peat wetland

The restoration of wetlands as bird habitats often involves the maintenance of a fluctuating water regime by careful, localised ditch water management using pumps and sluices. However, there is evidence in the literature to Suggest that alternate flood/drainage cycles can accelerate nutrient cycling and transport within the soil and, therefore, pose a threat to water quality through the process of eutrophication. This study focused on the dynamics and losses of soil P in a recently re-wetted, eutrophic fen peat developed on alluvium ill South west England. During the 2-year Study (2001 and 2002), soil water tensiometry revealed that the field water table (fluctuating annually between +20 and 60 cm relative to ground level) was extensively influenced across most of the 8.4 ha field site by the management of the adjacent ditch water levels. This conservation-led, prescribed water balance was facilitated by the high hydraulic conductivity (1.1 x 10(-s) ms(-1)) of the lower (70-140 cm), degraded layer of peat. However, only during a 7-day period of water table drawdown by intermittent pump drainage, approximately 45 g ha(-1) of dissolved reactive P (DRP) entered the pumped ditch from the field via this degraded layer. Summer rainfall events >35 mm d(-1) also coincided with significant peaks ill ditch water P concentration (up to 200 mu g L-1 DRP). Even larger peaks (Up to 700 mu g L-1 DRP) Occurred With the annual onset of autumn reflooding. These episodic P loss events pose a serious potential threat to biological water quality. (C) 2009 Elsevier B.V. All rights reserved.

Recent habitat degradation in karstic Lake Uluabat, western Turkey: A coupled limnological-palaeolimnological approach

The Ramsar site of Lake Uluabat, western Turkey, suffers from eutrophication, urban and industrial pollution and water abstraction, and its water levels are managed artificially. Here we combine monitoring and palaeolimnological. techniques to investigate spatial and temporal limnological variability and ecosystem impact, using an ostracod and mollusc survey to strengthen interpretation of the fossil record. A combination of low invertebrate Biological Monitoring Working Party scores (<10) and the dominance of eutrophic diatoms in the modern lake confirms its poor ecological status. Palaeolimnological analysis of recent (last >200 yr) changes in organic and carbonate content, diatoms, stable isotopes, ostracods and molluscs in a lake sediment core (UL20A) indicates a 20th century trend towards increased sediment accumulation rates and eutrophication which was probably initiated by deforestation and agriculture. The most marked ecological shift occurs in the early 1960s, however. A subtle rise in diatom-inferred total phosphorus and an inferred reduction in submerged aquatic macrophyte cover accompanies a major increase in sediment accumulation rate. An associated marked shift in ostracod stable isotope data indicative of reduced seasonality and a change in hydrological input suggests major impact from artificial water management practices, all of which appears to have culminated in the sustained loss of submerged macrophytes since 2000. The study indicates it is vital to take both land-use and water management practices into account in devising restoration strategies. in a wider context, the results have important implications for the conservation of shallow karstic lakes, the functioning of which is still poorly understood. (c) 2008 Elsevier Ltd. All rights reserved.

Calibrating flood inundation models: methods for assessing model performance

Satellite observed data for flood events have been used to calibrate and validate flood inundation models, providing valuable information on the spatial extent of the flood. Improvements in the resolution of this satellite imagery have enabled indirect remote sensing of water levels by using an underlying LiDAR DEM to extract the water surface elevation at the flood margin. Further to comparison of the spatial extent, this now allows for direct comparison between modelled and observed water surface elevations. Using a 12.5m ERS-1 image of a flood event in 2006 on the River Dee, North Wales, UK, both of these data types are extracted and each assessed for their value in the calibration of flood inundation models. A LiDAR guided snake algorithm is used to extract an outline of the flood from the satellite image. From the extracted outline a binary grid of wet / dry cells is created at the same resolution as the model, using this the spatial extent of the modelled and observed flood can be compared using a measure of fit between the two binary patterns of flooding. Water heights are extracted using points at intervals of approximately 100m along the extracted outline, and the students T-test is used to compare modelled and observed water surface elevations. A LISFLOOD-FP model of the catchment is set up using LiDAR topographic data resampled to the 12.5m resolution of the satellite image, and calibration of the friction parameter in the model is undertaken using each of the two approaches. Comparison between the two approaches highlights the sensitivity of the spatial measure of fit to uncertainty in the observed data and the potential drawbacks of using the spatial extent when parts of the flood are contained by the topography.

Estimating discharge from inundation imagery

Remote sensing from space-borne platforms is often seen as an appealing method of monitoring components of the hydrological cycle, including river discharge, due to its spatial coverage. However, data from these platforms is often less than ideal because the geophysical properties of interest are rarely measured directly and the measurements that are taken can be subject to significant errors. This study assimilated water levels derived from a TerraSAR-X synthetic aperture radar image and digital aerial photography with simulations from a two dimensional hydraulic model to estimate discharge, inundation extent, depths and velocities at the confluence of the rivers Severn and Avon, UK. An ensemble Kalman filter was used to assimilate spot heights water levels derived by intersecting shorelines from the imagery with a digital elevation model. Discharge was estimated from the ensemble of simulations using state augmentation and then compared with gauge data. Assimilating the real data reduced the error between analyzed mean water levels and levels from three gauging stations to less than 0.3 m, which is less than typically found in post event water marks data from the field at these scales. Measurement bias was evident, but the method still provided a means of improving estimates of discharge for high flows where gauge data are unavailable or of poor quality. Posterior estimates of discharge had standard deviations between 63.3 m3s-1 and 52.7 m3s-1, which were below 15% of the gauged flows along the reach. Therefore, assuming a roughness uncertainty of 0.03-0.05 and no model structural errors discharge could be estimated by the EnKF with accuracy similar to that arguably expected from gauging stations during flood events. Quality control prior to assimilation, where measurements were rejected for being in areas of high topographic slope or close to tall vegetation and trees, was found to be essential. The study demonstrates the potential, but also the significant limitations of currently available imagery to reduce discharge uncertainty in un-gauged or poorly gauged basins when combined with model simulations in a data assimilation framework.

The impact of uncertainty in satellite data on the assessment of flood inundation models

The performance of flood inundation models is often assessed using satellite observed data; however these data have inherent uncertainty. In this study we assess the impact of this uncertainty when calibrating a flood inundation model (LISFLOOD-FP) for a flood event in December 2006 on the River Dee, North Wales, UK. The flood extent is delineated from an ERS-2 SAR image of the event using an active contour model (snake), and water levels at the flood margin calculated through intersection of the shoreline vector with LiDAR topographic data. Gauged water levels are used to create a reference water surface slope for comparison with the satellite-derived water levels. Residuals between the satellite observed data points and those from the reference line are spatially clustered into groups of similar values. We show that model calibration achieved using pattern matching of observed and predicted flood extent is negatively influenced by this spatial dependency in the data. By contrast, model calibration using water elevations produces realistic calibrated optimum friction parameters even when spatial dependency is present. To test the impact of removing spatial dependency a new method of evaluating flood inundation model performance is developed by using multiple random subsamples of the water surface elevation data points. By testing for spatial dependency using Moran’s I, multiple subsamples of water elevations that have no significant spatial dependency are selected. The model is then calibrated against these data and the results averaged. This gives a near identical result to calibration using spatially dependent data, but has the advantage of being a statistically robust assessment of model performance in which we can have more confidence. Moreover, by using the variations found in the subsamples of the observed data it is possible to assess the effects of observational uncertainty on the assessment of flooding risk.

Characterization of a highly biodiverse floodplain meadow using hyperspectral remote sensing within a plant functional trait framework

We assessed the potential for using optical functional types as effective markers to monitor changes in vegetation in floodplain meadows associated with changes in their local environment. Floodplain meadows are challenging ecosystems for monitoring and conservation because of their highly biodiverse nature. Our aim was to understand and explain spectral differences among key members of floodplain meadows and also characterize differences with respect to functional traits. The study was conducted on a typical floodplain meadow in UK (MG4-type, mesotrophic grassland type 4, according to British National Vegetation Classification). We compared two approaches to characterize floodplain communities using field spectroscopy. The first approach was sub-community based, in which we collected spectral signatures for species groupings indicating two distinct eco-hydrological conditions (dry and wet soil indicator species). The other approach was “species-specific”, in which we focused on the spectral reflectance of three key species found on the meadow. One herb species is a typical member of the MG4 floodplain meadow community, while the other two species, sedge and rush, represent wetland vegetation. We also monitored vegetation biophysical and functional properties as well as soil nutrients and ground water levels. We found that the vegetation classes representing meadow sub-communities could not be spectrally distinguished from each other, whereas the individual herb species was found to have a distinctly different spectral signature from the sedge and rush species. The spectral differences between these three species could be explained by their observed differences in plant biophysical parameters, as corroborated through radiative transfer model simulations. These parameters, such as leaf area index, leaf dry matter content, leaf water content, and specific leaf area, along with other functional parameters, such as maximum carboxylation capacity and leaf nitrogen content, also helped explain the species’ differences in functional dynamics. Groundwater level and soil nitrogen availability, which are important factors governing plant nutrient status, were also found to be significantly different for the herb/wetland species’ locations. The study concludes that spectrally distinguishable species, typical for a highly biodiverse site such as a floodplain meadow, could potentially be used as target species to monitor vegetation dynamics under changing environmental conditions.

Phosphorus flux from wetland ditch sediments

The accumulation of phosphorus (P) in the bottom sediment of field drainage ditches poses a threat to the ecology both of the ditch water and downstream water courses. We investigated the amounts, forms and internal loading of sediment-bound P along two drainage ditches that regulate water levels in a basin fen (~ 200 ha) supporting a mixture of restored wetland and drained agricultural fields. Water levels in the Lady's Drove Rhyne are currently managed to enhance the biodiversity of the wetland (Catcott Lows Reserve — an area formerly cultivated for arable crop production); whereas, the East Ditch is managed to drain adjoining land that remains under arable and livestock production. Laboratory-based chemical fractionation schemes were used to characterise the forms and potential mobility of the sediment-bound P, whilst pore-water equilibrators were employed in situ to evaluate the diffusive flux of P through the sediment–water column, and to characterise the corresponding redox conditions. Along both ditches, sediment pore-water profiles indicated conditions ranging from weakly to very reducing conditions with increasing depth, and net fluxes of P from the sediment to overlying water. P flux values ranged from 0.33 to 1.30 mg m− 2 day− 1. Both the degree of P saturation (DPS) of the sediment and NaOH extractable (Fe/Al-bound) P correlated significantly (P < 0.05) with P flux. Both in the wetland and agricultural ditches, by far the highest values for P flux were recorded at sites closest to points of drainage water entry from the corresponding, adjoining land. Although the P flux data were obtained from only a single sampling event, this study highlights the contribution of historical as well as ongoing agricultural land use on the sustained elevated P status of ditch sediments in lowland catchments.

Hydrology and climatology at Laguna La Gaiba, lowland Bolivia: complex responses to climatic forcings over the last 25,000 years

Diatom, geochemical and isotopic data provide a record of environmental change in Laguna La Gaiba, lowland Bolivia (17°450S, 57°350W), over the last ca. 25 000 years. High-resolution diatom analysis around the Last Glacial–Interglacial Transition provides new insights into this period of change. The full and late glacial lake was generally quite shallow, but with evidence of periodic flooding. At about 13 100 cal a BP, just before the start of the Younger Dryas chronozone, the diatoms indicate shallower water conditions, but there is a marked change at about 12 200 cal a BP indicating the onset of a period of high variability, with rising water levels punctuated by periodic drying. From ca. 11 800 to 10 000 cal a BP, stable, deeper water conditions persisted. There is evidence for drying in the early to middle Holocene, but not as pronounced as that reported from elsewhere in the southern hemisphere tropics of South America. This was followed by the onset of wetter conditions in the late Holocene consistent with insolation forcing. Conditions very similar to present were established about 2100 cal a BP. A complex response to both insolation forcing and millennial-scale events originating in the North Atlantic is noted.

Vibrational energy levels of water

Several quartic force fields and a full sextic anharmonic force field for H,O have been determined from high-quality ab initio calculations, the highest at the aug-cc-pVQZ CCSD(T) level of theory. These force fields have been used to determine vibrational excited state band origins up to 15 000 cm - ’ above the zero-point level, using both a perturbation-resonancea pproach and a variational approach. An optimisedq uartic force field hasb eeno btained by least squares refinement of our best ab initio results to fit the observed overtone levels of 5 symmetrically substituted isotopomers of water (Hi60, Hi70, HisO, D,O, and T,O) with an rms error of less than 10 cm-‘, using the perturbation-resonancem odel for the vibrational calculation. Predicatel east squaresr efinement was usedt o provide a loose constraint of the refined force field to the ab initio results. The results obtained prove the viability of the perturbation-resonancem odel for usei n larger molecular systemsa nd also highlight someo f its weaknesse

Processes influencing ozone levels in Alaskan forest fire plumes during long-range transport over the North Atlantic

A case of long-range transport of a biomass burning plume from Alaska to Europe is analyzed using a Lagrangian approach. This plume was sampled several times in the free troposphere over North America, the North Atlantic and Europe by three different aircraft during the IGAC Lagrangian 2K4 experiment which was part of the ICARTT/ITOP measurement intensive in summer 2004. Measurements in the plume showed enhanced values of CO, VOCs and NOy, mainly in form of PAN. Observed O3 levels increased by 17 ppbv over 5 days. A photochemical trajectory model, CiTTyCAT, was used to examine processes responsible for the chemical evolution of the plume. The model was initialized with upwind data and compared with downwind measurements. The influence of high aerosol loading on photolysis rates in the plume was investigated using in situ aerosol measurements in the plume and lidar retrievals of optical depth as input into a photolysis code (Fast-J), run in the model. Significant impacts on photochemistry are found with a decrease of 18% in O3 production and 24% in O3 destruction over 5 days when including aerosols. The plume is found to be chemically active with large O3 increases attributed primarily to PAN decomposition during descent of the plume toward Europe. The predicted O3 changes are very dependent on temperature changes during transport and also on water vapor levels in the lower troposphere which can lead to O3 destruction. Simulation of mixing/dilution was necessary to reproduce observed pollutant levels in the plume. Mixing was simulated using background concentrations from measurements in air masses in close proximity to the plume, and mixing timescales (averaging 6.25 days) were derived from CO changes. Observed and simulated O3/CO correlations in the plume were also compared in order to evaluate the photochemistry in the model. Observed slopes change from negative to positive over 5 days. This change, which can be attributed largely to photochemistry, is well reproduced by multiple model runs even if slope values are slightly underestimated suggesting a small underestimation in modeled photochemical O3 production. The possible impact of this biomass burning plume on O3 levels in the European boundary layer was also examined by running the model for a further 5 days and comparing with data collected at surface sites, such as Jungfraujoch, which showed small O3 increases and elevated CO levels. The model predicts significant changes in O3 over the entire 10 day period due to photochemistry but the signal is largely lost because of the effects of dilution. However, measurements in several other BB plumes over Europe show that O3 impact of Alaskan fires can be potentially significant over Europe.

A review of the potential impacts of climate change on surface water quality

It is now accepted that some human-induced climate change is unavoidable. Potential impacts on water supply have received much attention, but relatively little is known about the concomitant changes in water quality. Projected changes in air temperature and rainfall could affect river flows and, hence, the mobility and dilution of contaminants. Increased water temperatures will affect chemical reaction kinetics and, combined with deteriorations in quality, freshwater ecological status. With increased flows there will be changes in stream power and, hence, sediment loads with the potential to alter the morphology of rivers and the transfer of sediments to lakes, thereby impacting freshwater habitats in both lake and stream systems. This paper reviews such impacts through the lens of UK surface water quality. Widely accepted climate change scenarios suggest more frequent droughts in summer, as well as flash-flooding, leading to uncontrolled discharges from urban areas to receiving water courses and estuaries. Invasion by alien species is highly likely, as is migration of species within the UK adapting to changing temperatures and flow regimes. Lower flows, reduced velocities and, hence, higher water residence times in rivers and lakes will enhance the potential for toxic algal blooms and reduce dissolved oxygen levels. Upland streams could experience increased dissolved organic carbon and colour levels, requiring action at water treatment plants to prevent toxic by-products entering public water supplies. Storms that terminate drought periods will flush nutrients from urban and rural areas or generate acid pulses in acidified upland catchments. Policy responses to climate change, such as the growth of bio-fuels or emission controls, will further impact freshwater quality.