Uses of airborne scanning laser altimetry for river flood modelling and coastal geomorphology

Two ongoing projects at ESSC that involve the development of new techniques for extracting information from airborne LiDAR data and combining this information with environmental models will be discussed. The first project in conjunction with Bristol University is aiming to improve 2-D river flood flow models by using remote sensing to provide distributed data for model calibration and validation. Airborne LiDAR can provide such models with a dense and accurate floodplain topography together with vegetation heights for parameterisation of model friction. The vegetation height data can be used to specify a friction factor at each node of a model’s finite element mesh. A LiDAR range image segmenter has been developed which converts a LiDAR image into separate raster maps of surface topography and vegetation height for use in the model. Satellite and airborne SAR data have been used to measure flood extent remotely in order to validate the modelled flood extent. Methods have also been developed for improving the models by decomposing the model’s finite element mesh to reflect floodplain features such as hedges and trees having different frictional properties to their surroundings. Originally developed for rural floodplains, the segmenter is currently being extended to provide DEMs and friction parameter maps for urban floods, by fusing the LiDAR data with digital map data. The second project is concerned with the extraction of tidal channel networks from LiDAR. These networks are important features of the inter-tidal zone, and play a key role in tidal propagation and in the evolution of salt-marshes and tidal flats. The study of their morphology is currently an active area of research, and a number of theories related to networks have been developed which require validation using dense and extensive observations of network forms and cross-sections. The conventional method of measuring networks is cumbersome and subjective, involving manual digitisation of aerial photographs in conjunction with field measurement of channel depths and widths for selected parts of the network. A semi-automatic technique has been developed to extract networks from LiDAR data of the inter-tidal zone. A multi-level knowledge-based approach has been implemented, whereby low level algorithms first extract channel fragments based mainly on image properties then a high level processing stage improves the network using domain knowledge. The approach adopted at low level uses multi-scale edge detection to detect channel edges, then associates adjacent anti-parallel edges together to form channels. The higher level processing includes a channel repair mechanism.

The relationship between multi-species catch and effort: among fishery comparisons

The control of fishing mortality via fishing effort remains fundamental to most fisheries management strategies even at the local community or co-management level. Decisions to support such strategies require knowledge of the underlying response of the catch to changes in effort. Even under adaptive management strategies, imprecise knowledge of the response is likely to help accelerate the adaptive learning process. Data and institutional capacity requirements to employ multi-species biomass dynamics and age-structured models invariably render their use impractical particularly in less developed regions of the world. Surplus production models fitted to catch and effort data aggregated across all species offer viable alternatives. The current paper seeks models of this type that best describe the multi-species catch–effort responses in floodplain-rivers, lakes and reservoirs and reef-based fisheries based upon among fishery comparisons, building on earlier work. Three alternative surplus production models were fitted to estimates of catch per unit area (CPUA) and fisher density for 258 fisheries in Africa, Asia and South America. In all cases examined, the best or equal best fitting model was the Fox type, explaining up to 90% of the variation in CPUA. For lake and reservoir fisheries in Africa and Asia, the Schaefer and an asymptotic model fitted equally well. The Fox model estimates of fisher density (fishers km−2) at maximum yield (iMY) for floodplain-rivers, African lakes and reservoirs and reef-based fisheries are 13.7 (95% CI [11.8, 16.4]); 27.8 (95% CI [17.5, 66.7]) and 643 (95% CI [459,1075]), respectively and compare well with earlier estimates. Corresponding estimates of maximum yield are also given. The significantly higher value of iMY for reef-based fisheries compared to estimates for rivers and lakes reflects the use of a different measure of fisher density based upon human population size estimates. The models predict that maximum yield is achieved at a higher fishing intensity in Asian lakes compared to those in Africa. This may reflect the common practice in Asia of stocking lakes to augment natural recruitment. Because of the equilibrium assumptions underlying the models, all the estimates of maximum yield and corresponding levels of effort should be treated with caution.

Evaluating a new LISFLOOD-FP formulation with data from the summer 2007 floods in Tewkesbury, UK

Two-dimensional flood inundation modelling is a widely used tool to aid flood risk management. In urban areas, the model spatial resolution required to represent flows through a typical street network often results in an impractical computational cost at the city scale. This paper presents the calibration and evaluation of a recently developed formulation of the LISFLOOD-FP model, which is more computationally efficient at these resolutions. Aerial photography was available for model evaluation on 3 days from the 24 to the 31 of July. The new formulation was benchmarked against the original version of the model at 20 and 40 m resolutions, demonstrating equally accurate simulation, given the evaluation data but at a 67 times faster computation time. The July event was then simulated at the 2 m resolution of the available airborne LiDAR DEM. This resulted in more accurate simulation of the floodplain drying dynamics compared with the coarse resolution models, although maximum inundation levels were simulated equally well at all resolutions tested.

Holocene palaeoenvironmental history and the impact of prehistoric salt production in the Seille Valley, Eastern France

The Seille Valley in eastern France was home to one of Europe’s largest Iron Age salt industries. Sedimentology, palynology and geochronology have been integrated within ongoing archaeological investigations to reconstruct the Holocene palaeoenvironmental history of the Seille Valley and to elucidate the human–environment relationship of salt production. A sedimentary model of the valley has been constructed from a borehole survey of the floodplain and pollen analyses have been undertaken to reconstruct the vegetation history. Alluvial records have been successfully dated using optically stimulated luminescence and radiocarbon techniques, thereby providing a robust chronological framework. The results have provided an insight into the development of favourable conditions for salt production and there is evidence in the sedimentary record to suggest that salt production may have taken place during the mid-to-late Bronze Age. The latter has yet to be identified in the archaeological record and targeted excavation is therefore underway to test this finding. The development of the Iron Age industry had a major impact on the hydrological regime of the valley and its sedimentological history, with evidence for accelerated alluviation arising from floodplain erosion at salt production sites and modification of the local fluvial regime due to briquetage accumulation on the floodplain. This research provides an important insight into the environmental implications of early industrial activities, in addition to advancing knowledge about the Holocene palaeoenvironmental and social history of this previously poorly studied region of France.

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.

Faunal evidence of a Holocene pluvial phase in Southern Arabia with remarks on the morphological variability of Helenina anderseni

Although foraminifera have been found living in inland saline lakes isolated from the sea, this phenomenon has rarely been recognized in the fossil record. This study documents the occurrence of benthic foraminifera in Holocene lake sediments located nearly 500 km inland from the Red Sea, in the Al-Mundafan region of southern Saudi Arabia. The lake formed during a regional pluvial period, 10,500–6000 yr BP. The presence of foraminifera and brackish charophytes in the studied section represent an interval when the lake was slightly brackish due to high evaporation. The studied sediments yielded a bispecific benthic foraminiferal fauna comprised of Helenina anderseni and Trichohyalus aguayoi, as well as the brackish charophyte genus Lamprothamnium. The benthic foraminifera are species characteristic of mangrove swamps, salt marshes, and lagoons, which are environments currently widespread along the Red Sea coasts. Because the Al Mundafan area was never connected to the sea during the Quaternary, wading birds must have been the vector that transported the foraminifera to the paleolake.

Scheduling satellite-based SAR acquisition for sequential assimilation of water level observations into flood modelling

Satellite-based Synthetic Aperture Radar (SAR) has proved useful for obtaining information on flood extent, which, when intersected with a Digital Elevation Model (DEM) of the floodplain, provides water level observations that can be assimilated into a hydrodynamic model to decrease forecast uncertainty. With an increasing number of operational satellites with SAR capability, information on the relationship between satellite first visit and revisit times and forecast performance is required to optimise the operational scheduling of satellite imagery. By using an Ensemble Transform Kalman Filter (ETKF) and a synthetic analysis with the 2D hydrodynamic model LISFLOOD-FP based on a real flooding case affecting an urban area (summer 2007,Tewkesbury, Southwest UK), we evaluate the sensitivity of the forecast performance to visit parameters. We emulate a generic hydrologic-hydrodynamic modelling cascade by imposing a bias and spatiotemporal correlations to the inflow error ensemble into the hydrodynamic domain. First, in agreement with previous research, estimation and correction for this bias leads to a clear improvement in keeping the forecast on track. Second, imagery obtained early in the flood is shown to have a large influence on forecast statistics. Revisit interval is most influential for early observations. The results are promising for the future of remote sensing-based water level observations for real-time flood forecasting in complex scenarios.

Marine Isotope Stage 9 environments of fluvial deposits at Hackney, north London, UK

Middle Pleistocene deposits at Hackney, north London comprise a thick unit of organic sands and silts occupying a channel near the confluence of the River Thames in south-eastern England and its left-bank tributary the River Lea. They represent a short time interval, perhaps no more than a few years, within a late Middle Pleistocene interglacial. The organic sediments are overlain by unfossiliferous sands and gravels indicating deposition on the floodplain of a braided river under cool or cold climatic conditions. The fossil plant, insect, mollusc and vertebrate remains from the interglacial deposits all indicate climatic conditions with summers warmer than the present in SE England, and winters with a similar thermal climate. The biostratigraphic evidence suggests that the time period represented by the organic unit is part of MIS 9, although the geochronological evidence for such an age is inconclusive. The palaeontological evidence strongly suggests that this temperate stage was warmer than the succeeding temperate stage MIS 7 or the Holocene, and approaching the Ipswichian (MISs 5e) in its warmth. The multidisciplinary description of the Hackney deposits is one of the first to reconstruct terrestrial conditions in Marine Isotope Stage 9 in Western Europe.

Reducing inconsistencies in point observations of maximum flood inundation level

Flood simulation models and hazard maps are only as good as the underlying data against which they are calibrated and tested. However, extreme flood events are by definition rare, so the observational data of flood inundation extent are limited in both quality and quantity. The relative importance of these observational uncertainties has increased now that computing power and accurate lidar scans make it possible to run high-resolution 2D models to simulate floods in urban areas. However, the value of these simulations is limited by the uncertainty in the true extent of the flood. This paper addresses that challenge by analyzing a point dataset of maximum water extent from a flood event on the River Eden at Carlisle, United Kingdom, in January 2005. The observation dataset is based on a collection of wrack and water marks from two postevent surveys. A smoothing algorithm for identifying, quantifying, and reducing localized inconsistencies in the dataset is proposed and evaluated showing positive results. The proposed smoothing algorithm can be applied in order to improve flood inundation modeling assessment and the determination of risk zones on the floodplain.

Holocene alluvial environments at Barking, Lower Thames Valley (London, UK)

Deposit modelling based on archived borehole logs supplemented by a small number of dedicated boreholes is used to reconstruct the main boundary surfaces and the thickness of the main sediment units within the succession of Holocene alluvial deposits underlying the floodplain in the Barking Reach of the Lower Thames Valley. The basis of the modelling exercise is discussed and the models are used to assess the significance of floodplain relief in determining patterns of sedimentation. This evidence is combined with the results of biostratigraphical and geochronological investigations to reconstruct the environmental conditions associated with each successive stage of floodplain aggradation. The two main factors affecting the history and spatial pattern of Holocene sedimentation are shown to be the regional behaviour of relative sea level and the pattern of relief on the surface of the sub-alluvial, Late Devensian Shepperton Gravel. As is generally the case in the Lower Thames Valley, three main stratigraphic units are recognised, the Lower Alluvium, a peat bed broadly equivalent to the Tilbury III peat of Devoy (1979) and an Upper Alluvium. There is no evidence to suggest that the floodplain was substantially re-shaped by erosion during the Holocene. Instead, the relief inherited from the Shepperton Gravel surface was gradually buried either by the accumulation of peat or by deposition of fine-grained sediment from suspension in standing or slow-moving water. The palaeoenvironmental record from Barking confirms important details of the Holocene record observed elsewhere in the Lower Thames Valley, including the presence of Taxus in the valley-floor fen carr woodland between about 5000 and 4000 cal BP, and the subsequent growth of Ulmus on the peat surface.

Geophysical characterization of the complex dynamics of groundwater and seawater exchange in a highly stressed aquifer system linked to a coastal lagoon (SE Spain)

Anthropogenic pressure influences the two-way interactions between shallow aquifers and coastal lagoons. Aquifer overexploitation may lead to seawater intrusion, and aquifer recharge from rainfall plus irrigation may, in turn, increase the groundwater discharge into the lagoon. We analyse the evolution, since the 1950s up to the present, of the interactions between the Campo de Cartagena Quaternary aquifer and the Mar Menor coastal lagoon (SE Spain). This is a very heterogeneous and anisotropic detrital aquifer, where aquifer–lagoon interface has a very irregular geometry. Using electrical resistivity tomography, we clearly identified the freshwater–saltwater transition zone and detected areas affected by seawater intrusion. Severity of the intrusion was spatially variable and significantly related to the density of irrigation wells in 1950s–1960s, suggesting the role of groundwater overexploitation. We distinguish two different mechanisms by which water from the sea invades the land: (a) horizontal advance of the interface due to a wide exploitation area and (b) vertical rise (upconing) caused by local intensive pumping. In general, shallow parts of the geophysical profiles show higher electrical resistivity associated with freshwater mainly coming from irrigation return flows, with water resources mostly from deep confined aquifers and imported from Tagus river, 400 km north. This indicates a likely reversal of the former seawater intrusion process.

Sensitivity of a hydraulic model to changes in channel erosion during extreme flooding

Recent research into flood modelling has primarily concentrated on the simulation of inundation flow without considering the influences of channel morphology. River channels are often represented by a simplified geometry that is implicitly assumed to remain unchanged during flood simulations. However, field evidence demonstrates that significant morphological changes can occur during floods to mobilise the boundary sediments. Despite this, the effect of channel morphology on model results has been largely unexplored. To address this issue, the impact of channel cross-section geometry and channel long-profile variability on flood dynamics is examined using an ensemble of a 1D-2D hydraulic model (LISFLOOD-FP) of the 1:2102 year recurrence interval floods in Cockermouth, UK, within an uncertainty framework. A series of hypothetical scenarios of channel morphology were constructed based on a simple velocity based model of critical entrainment. A Monte-Carlo simulation framework was used to quantify the effects of channel morphology together with variations in the channel and floodplain roughness coefficients, grain size characteristics, and critical shear stress on measures of flood inundation. The results showed that the bed elevation modifications generated by the simplistic equations reflected a good approximation of the observed patterns of spatial erosion despite its overestimation of erosion depths. The effect of uncertainty on channel long-profile variability only affected the local flood dynamics and did not significantly affect the friction sensitivity and flood inundation mapping. The results imply that hydraulic models generally do not need to account for within event morphodynamic changes of the type and magnitude modelled, as these have a negligible impact that is smaller than other uncertainties, e.g. boundary conditions. Instead morphodynamic change needs to happen over a series of events to become large enough to change the hydrodynamics of floods in supply limited gravel-bed rivers like the one used in this research.

Satellite-supported flood forecasting in river networks: a real case study

Satellite-based (e.g., Synthetic Aperture Radar [SAR]) water level observations (WLOs) of the floodplain can be sequentially assimilated into a hydrodynamic model to decrease forecast uncertainty. This has the potential to keep the forecast on track, so providing an Earth Observation (EO) based flood forecast system. However, the operational applicability of such a system for floods developed over river networks requires further testing. One of the promising techniques for assimilation in this field is the family of ensemble Kalman (EnKF) filters. These filters use a limited-size ensemble representation of the forecast error covariance matrix. This representation tends to develop spurious correlations as the forecast-assimilation cycle proceeds, which is a further complication for dealing with floods in either urban areas or river junctions in rural environments. Here we evaluate the assimilation of WLOs obtained from a sequence of real SAR overpasses (the X-band COSMO-Skymed constellation) in a case study. We show that a direct application of a global Ensemble Transform Kalman Filter (ETKF) suffers from filter divergence caused by spurious correlations. However, a spatially-based filter localization provides a substantial moderation in the development of the forecast error covariance matrix, directly improving the forecast and also making it possible to further benefit from a simultaneous online inflow error estimation and correction. Additionally, we propose and evaluate a novel along-network metric for filter localization, which is physically-meaningful for the flood over a network problem. Using this metric, we further evaluate the simultaneous estimation of channel friction and spatially-variable channel bathymetry, for which the filter seems able to converge simultaneously to sensible values. Results also indicate that friction is a second order effect in flood inundation models applied to gradually varied flow in large rivers. The study is not conclusive regarding whether in an operational situation the simultaneous estimation of friction and bathymetry helps the current forecast. Overall, the results indicate the feasibility of stand-alone EO-based operational flood forecasting.

A palaeoenvironmental context for terminal upper Palaeolithic and Mesolithic activity in the Colne Valley: offsite records contemporary with occupation at three ways Wharf, Uxbridge

Multi-proxy analyses from floodplain deposits in the Colne Valley, southern England, have provided a palaeoenvironmental context for the immediately adjacent Terminal Upper Palaeolithic and Early Mesolithic site of Three Ways Wharf. These deposits show the transition from an open cool environment to fully developed heterogeneous floodplain vegetation during the Early Mesolithic. Several distinct phases of burning are shown to have occurred that are chronologically contemporary with the local archaeological record. The floodplain itself is shown to have supported a number of rare Urwaldrelikt insect species implying human manipulation of the floodplain at this time must have been limited or episodic. By the Late Mesolithic a reed-sedge swamp had developed across much of the floodplain, within which repeated burning of the in situ vegetation took place. This indicates deliberate land management practices utilising fire, comparable with findings from other floodplain sequences in southern Britain. With similar sedimentary sequences known to exist across the Colne Valley, often closely associated with contemporary archaeology, the potential for placing the archaeological record within a spatially explicit palaeoenvironmental context is great.