Floodplain environmental change during the Younger Dryas and Holocene in Northwest Europe: Insights from the lower Kennet Valley, south central England

Many lowland rivers across northwest Europe exhibit broadly similar behavioural responses to glacial-interglacial transitions and landscape development. Difficulties exist in assessing these, largely because the evidence from many rivers remains limited and fragmentary. Here we address this issue in the context of the river Kennet, a tributary of the Thames, since c. 13,000 cal BP. Some similarities with other rivers are present, suggesting that regional climatic shifts are important controls. The Kennet differs from the regional pattern in a number of ways. The rate of response to sudden climatic change, particularly at the start of the Holocene and also mid-Holocene forest clearance, appears very high. This may reflect abrupt shifts between two catchment scale hydrological states arising from contemporary climates, land use change and geology. Stadial hydrology is dominated by nival regimes, with limited winter infiltration and high spring and summer runoff. Under an interglacial climate, infiltration is more significant. The probable absence of permafrost in the catchment means that a lag between the two states due to its gradual decay is unlikely. Palaeoecology, supported by radiocarbon dates, suggests that, at the very start of the Holocene, a dramatic episode of fine sediment deposition across most of the valley floor occurred, lasting 500-1000 years. A phase of peat accumulation followed as mineral sediment supply declined. A further shift led to tufa deposition, initially in small pools, then across the whole floodplain area, with the river flowing through channels cut in tufa and experiencing repeated avulsion. Major floods, leaving large gravel bars that still form positive relief features on the floodplain, followed mid-Holocene floodplain stability. Prehistoric deforestation is likely to be the cause of this flooding, inducing a major environmental shift with significantly increased surface runoff. Since the Bronze Age, predominantly fine sediments were deposited along the valley with apparently stable channels and vertical floodplain accretion associated with soil erosion and less catastrophic flooding. The Kennet demonstrates that, while a general pattern of river behaviour over time, within a region, may be identifiable, individual rivers are likely to diverge from this. Consequently, it is essential to understand catchment controls, particularly the relative significance of surface and subsurface hydrology. (c) 2005 Elsevier B.V. All rights reserved.

Rice versus fish revisited: On the integrated management of floodplain resources in Bangladesh

Disproportionately little attention has been paid to the dry season trade-off between rice and (inland capture) fish production on the floodplains of Bangladesh, compared to the same trade-off during the flood season. As the rural economy grows increasingly dominated by dry-season irrigated rice production, and floodplain land and water come under ever-increasing pressure during the dry winter months, there is an urgent need to focus attention on these dry months that are so critical to the survival and propagation of the floodplain resident fish, and to the poor people that depend on these fish for their livelihood. This article examines three important dry-season natural resource constraints to floodplain livelihoods in Bangladesh, and finds a common factor at the heart of all three: rice cultivation on lands at low and very low elevations. The article articulates the system interlinkages that bind these constraints and the long-run trend towards irrigated rice cropping on lower-lying lands, and suggests a management approach based on locally tailored strategies to arrest this trend. Apart from its direct relevance to the floodplains of Bangladesh, which support more than 100 million people, these lessons have relevance for river floodplain systems elsewhere in the developing world, notably the Mekong Delta.

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.

Multi‐method global sensitivity analysis (MMGSA) for modelling floodplain hydrological processes

When studying hydrological processes with a numerical model, global sensitivity analysis (GSA) is essential if one is to understand the impact of model parameters and model formulation on results. However, different definitions of sensitivity can lead to a difference in the ranking of importance of the different model factors. Here we combine a fuzzy performance function with different methods of calculating global sensitivity to perform a multi-method global sensitivity analysis (MMGSA). We use an application of a finite element subsurface flow model (ESTEL-2D) on a flood inundation event on a floodplain of the River Severn to illustrate this new methodology. We demonstrate the utility of the method for model understanding and show how the prediction of state variables, such as Darcian velocity vectors, can be affected by such a MMGSA. This paper is a first attempt to use GSA with a numerically intensive hydrological model.

Multi-method global sensitivity analysis (MMGSA) for modelling floodplain hydrological processes

When studying hydrological processes with a numerical model, global sensitivity analysis (GSA) is essential if one is to understand the impact of model parameters and model formulation on results. However, different definitions of sensitivity can lead to a difference in the ranking of importance of the different model factors. Here we combine a fuzzy performance function with different methods of calculating global sensitivity to perform a multi-method global sensitivity analysis (MMGSA). We use an application of a finite element subsurface flow model (ESTEL-2D) on a flood inundation event on a floodplain of the River Severn to illustrate this new methodology. We demonstrate the utility of the method for model understanding and show how the prediction of state variables, such as Darcian velocity vectors, can be affected by such a MMGSA. This paper is a first attempt to use GSA with a numerically intensive hydrological model

Use of fused airborne scanning laser altimetry and digital map data for urban flood modelling

Flood modelling of urban areas is still at an early stage, partly because until recently topographic data of sufficiently high resolution and accuracy have been lacking in urban areas. However, Digital Surface Models (DSMs) generated from airborne scanning laser altimetry (LiDAR) having sub-metre spatial resolution have now become available, and these are able to represent the complexities of urban topography. The paper describes the development of a LiDAR post-processor for urban flood modelling based on the fusion of LiDAR and digital map data. The map data are used in conjunction with LiDAR data to identify different object types in urban areas, though pattern recognition techniques are also employed. Post-processing produces a Digital Terrain Model (DTM) for use as model bathymetry, and also a friction parameter map for use in estimating spatially-distributed friction coefficients. In vegetated areas, friction is estimated from LiDAR-derived vegetation height, and (unlike most vegetation removal software) the method copes with short vegetation less than ~1m high, which may occupy a substantial fraction of even an urban floodplain. The DTM and friction parameter map may also be used to help to generate an unstructured mesh of a vegetated urban floodplain for use by a 2D finite element model. The mesh is decomposed to reflect floodplain features having different frictional properties to their surroundings, including urban features such as buildings and roads as well as taller vegetation features such as trees and hedges. This allows a more accurate estimation of local friction. The method produces a substantial node density due to the small dimensions of many urban features.

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.

Improving flood inundation models using remotely sensed data

The paper discusses the wide variety of ways in which remotely sensed data are being utilized in river flood inundation modeling. Model parameterization is being aided using airborne LiDAR data to provide topography of the floodplain for use as model bathymetry, and vegetation heights in the floodplain for use in estimating floodplain friction factors. Model calibration and validation are being aided by comparing the flood extent observed in SAR images with the extent predicted by the model. The recent extension of this to the observation of urban flooding using high resolution TerraSAR-X data is described. Possible future research directions are considered.

Trace fossils and pseudofossils from the Wealden strata (non-marine Lower Cretaceous) of southern England

The trace fossils of the Wealden (non-marine Lower Cretaceous) of southern England are described. Sixteen invertebrate ichnotaxa include Agrichnium fimbriatus, Beaconites antarcticus, B. barretti, Cochlichnus anguineus, Diplichnites triassicus, Diplocraterion parallelum, Lockeia siliquaria, L. serialis, Monocraterion cf. tentaculum, Palaeophycus striatus, P. tubularis, Planolites montanus, Protovirgularia rugosa, Rhizocorallium isp., Scoyenia cf. gracilis, Unisulcus minutus, insect and root traces. Tetrapod tracks and trackways include tridactyl Iguanodontipus burreyi and other ornithopods, theropod, and tetradactyl sauropod (or possibly ankylosaur), together with extensive dinosaur tramplings. Coprolites are referred to two broad types: spiral, with or without included fish scales (attributable to sharks), and elongate and irregular (possibly produced by reptiles). A skinprint and two types of pseudofossil are also included. Five environmental associations are recognised: (1) lacustrine/lagoonal; (2) brackish incursions (flooding events) into the lacustrine/lagoonal environment; (3) a marginal lacustrine association with fluvial input; (4) a fluvial (lacustrine delta) association; (5) floodplain sediments (seasonal wetlands). These associations are assigned to the fluvial-lacustrine Scoyenia Ichnofacies and the incursions to Glossifungites lchnofacies. (c) 2005 Elsevier Ltd. All rights reserved.

Case study of the use of remotely sensed data for modeling flood inundation on the river Severn, UK.

A methodology for using remotely sensed data to both generate and evaluate a hydraulic model of floodplain inundation is presented for a rural case study in the United Kingdom: Upton-upon-Severn. Remotely sensed data have been processed and assembled to provide an excellent test data set for both model construction and validation. In order to assess the usefulness of the data and the issues encountered in their use, two models for floodplain inundation were constructed: one based on an industry standard one-dimensional approach and the other based on a simple two-dimensional approach. The results and their implications for the future use of remotely sensed data for predicting flood inundation are discussed. Key conclusions for the use of remotely sensed data are that care must be taken to integrate different data sources for both model construction and validation and that improvements in ground height data shift the focus in terms of model uncertainties to other sources such as boundary conditions. The differences between the two models are found to be of minor significance.

Global scale water stages from SAR imagery (GlosSARi) to support global flood forecasting

Despite the success of studies attempting to integrate remotely sensed data and flood modelling and the need to provide near-real time data routinely on a global scale as well as setting up online data archives, there is to date a lack of spatially and temporally distributed hydraulic parameters to support ongoing efforts in modelling. Therefore, the objective of this project is to provide a global evaluation and benchmark data set of floodplain water stages with uncertainties and assimilation in a large scale flood model using space-borne radar imagery. An algorithm is developed for automated retrieval of water stages with uncertainties from a sequence of radar imagery and data are assimilated in a flood model using the Tewkesbury 2007 flood event as a feasibility study. The retrieval method that we employ is based on possibility theory which is an extension of fuzzy sets and that encompasses probability theory. In our case we first attempt to identify main sources of uncertainty in the retrieval of water stages from radar imagery for which we define physically meaningful ranges of parameter values. Possibilities of values are then computed for each parameter using a triangular ‘membership’ function. This procedure allows the computation of possible values of water stages at maximum flood extents along a river at many different locations. At a later stage in the project these data are then used in assimilation, calibration or validation of a flood model. The application is subsequently extended to a global scale using wide swath radar imagery and a simple global flood forecasting model thereby providing improved river discharge estimates to update the latter.

Improving flood models using LiDAR: progress and problems

Airborne scanning laser altimetry (LiDAR) is an important new data source for river flood modelling. LiDAR can give dense and accurate DTMs of floodplains for use as model bathymetry. Spatial resolutions of 0.5m or less are possible, with a height accuracy of 0.15m. LiDAR gives a Digital Surface Model (DSM), so vegetation removal software (e.g. TERRASCAN) must be used to obtain a DTM. An example used to illustrate the current state of the art will be the LiDAR data provided by the EA, which has been processed by their in-house software to convert the raw data to a ground DTM and separate vegetation height map. Their method distinguishes trees from buildings on the basis of object size. EA data products include the DTM with or without buildings removed, a vegetation height map, a DTM with bridges removed, etc. Most vegetation removal software ignores short vegetation less than say 1m high. We have attempted to extend vegetation height measurement to short vegetation using local height texture. Typically most of a floodplain may be covered in such vegetation. The idea is to assign friction coefficients depending on local vegetation height, so that friction is spatially varying. This obviates the need to calibrate a global floodplain friction coefficient. It’s not clear at present if the method is useful, but it’s worth testing further. The LiDAR DTM is usually determined by looking for local minima in the raw data, then interpolating between these to form a space-filling height surface. This is a low pass filtering operation, in which objects of high spatial frequency such as buildings, river embankments and walls may be incorrectly classed as vegetation. The problem is particularly acute in urban areas. A solution may be to apply pattern recognition techniques to LiDAR height data fused with other data types such as LiDAR intensity or multispectral CASI data. We are attempting to use digital map data (Mastermap structured topography data) to help to distinguish buildings from trees, and roads from areas of short vegetation. The problems involved in doing this will be discussed. A related problem of how best to merge historic river cross-section data with a LiDAR DTM will also be considered. LiDAR data may also be used to help generate a finite element mesh. In rural area we have decomposed a floodplain mesh according to taller vegetation features such as hedges and trees, so that e.g. hedge elements can be assigned higher friction coefficients than those in adjacent fields. We are attempting to extend this approach to urban area, so that the mesh is decomposed in the vicinity of buildings, roads, etc as well as trees and hedges. A dominant points algorithm is used to identify points of high curvature on a building or road, which act as initial nodes in the meshing process. A difficulty is that the resulting mesh may contain a very large number of nodes. However, the mesh generated may be useful to allow a high resolution FE model to act as a benchmark for a more practical lower resolution model. A further problem discussed will be how best to exploit data redundancy due to the high resolution of the LiDAR compared to that of a typical flood model. Problems occur if features have dimensions smaller than the model cell size e.g. for a 5m-wide embankment within a raster grid model with 15m cell size, the maximum height of the embankment locally could be assigned to each cell covering the embankment. But how could a 5m-wide ditch be represented? Again, this redundancy has been exploited to improve wetting/drying algorithms using the sub-grid-scale LiDAR heights within finite elements at the waterline.