Analysis of water-level response to rainfall and implications for recharge pathways in the Chalk aquifer, SE England

Water table response to rainfall was investigated at six sites in the Upper, Middle and Lower Chalk of southern England. Daily time series of rainfall and borehole water level were cross-corretated to investigate seasonal variations in groundwater-level response times, based on periods of 3-month duration. The time tags (in days) yielding significant correlations were compared with the average unsaturated zone thickness during each 3-month period. In general, for cases when the unsaturated zone was greater than 18 m thick, the time tag for a significant water-level response increased rapidly once the depth to the water table exceeded a critical value, which varied from site to site. For shallower water tables, a linear relationship between the depth to the water table and the water-level response time was evident. The observed variations in response time can only be partially accounted for using a diffusive model for propagation through the unsaturated matrix, suggesting that some fissure flow was occurring. The majority of rapid responses were observed during the winter/spring recharge period, when the unsaturated zone is thinnest and the unsaturated zone moisture content is highest, and were more likely to occur when the rainfall intensity exceeded 5 mm/day. At some sites, a very rapid response within 24 h of rainfall was observed in addition to the longer term responses even when the unsaturated zone was up to 64 m thick. This response was generally associated with the autumn period. The results of the cross-correlation analysis provide statistical support for the presence of fissure flow and for the contribution of multiple pathways through the unsaturated zone to groundwater recharge. (c) 2006 Elsevier B.V. All rights reserved.

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.

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.

Pediastrum species as potential indicators of lake-level change in tropical South America

We explored the potential for using Pediastrum (Meyen), a genus of green alga commonly found in palaeoecological studies, as a proxy for lake-level change in tropical South America. The study site, Laguna La Gaiba (LLG) (17°45′S, 57°40′W), is a broad, shallow lake located along the course of the Paraguay River in the Pantanal, a 135,000-km2 tropical wetland located mostly in western Brazil, but extending into eastern Bolivia. Fourteen surface sediment samples were taken from LLG across a range of lake depths (2-5.2 m) and analyzed for Pediastrum. We found seven species, of which P. musteri (Tell et Mataloni), P. argentiniense (Bourr. et Tell), and P. cf. angulosum (Ehrenb.) ex Menegh. were identified as potential indicators of lake level. Results of the modern dataset were applied to 31 fossil Pediastrum assemblages spanning the early Holocene (12.0 kyr BP) to present to infer past lake level changes qualitatively. Early Holocene (12.0-9.8 kyr BP) assemblages do not show a clear signal, though abundance of P. simplex (Meyen) suggests relatively high lake levels. Absence of P. musteri, characteristic of deep, open water, and abundance of macrophyte-associated taxa indicate lake levels were lowest from 9.8 to 3.0 kyr BP. A shift to wetter conditions began at 4.4 kyr BP, indicated by the appearance of P. musteri, though inferred lake levels did not reach modern values until 1.4 kyr BP. The Pediastrum-inferred mid-Holocene lowstand is consistent with lower precipitation, previously inferred using pollen from this site, and is also in agreement with evidence for widespread drought in the South American tropics during the middle Holocene. An inference for steadily increasing lake level from 4.4 kyr BP to present is consistent with diatom-inferred water level rise at Lake Titicaca, and demonstrates coherence with the broad pattern of increasing monsoon strength from the late Holocene until present in tropical South America.

Simulating the Holocene lake-level record of Lake Bysjon, southern Sweden

Lake Bysjön, southern Sweden, has experienced major lake-level lowerings during the Holocene, with one interval about 900014C yr B.P. when water level dropped ca. 7 m and the lake became closed. These changes were not solely due to known changes in radiation budgets or seasonal temperatures. Simulations with a lake-catchment model indicate that, given the actual changes in radiation and temperatures, all the observed lake-level lowerings (including the major lowering at 900014C yr B.P.) could have occurred in response to precipitation changes of <75 mm/yr when winter temperatures were warmer than today. In these circumstances, the reduction of runoff into the lake caused by increased evapotranspiration during the late winter and spring, combined with relatively small changes in precipitation, was sufficient for the lake to become closed. When winter temperatures were colder than today, the reduction in winter runoff related to reduced precipitation was only very slight and insufficient to lower the lake below threshold. In such circumstances, changes in outflow were sufficient to compensate for the combined changes in precipitation and runoff, and lake level therefore remained unchanged.

Annual textural banding in Holocene estuarine silts, Severn Estuary Levels (SW Britain): patterns, cause and implications

A high-resolution textural study has been made of laminated and banded estuarine silts exposed intertidally at representative localities and horizons in the Holocene deposits of the Severn Estuary Levels. The laminae, on a submillimetre to millimetre scale, are sharp-based, graded couplets formed of a lower silty part overlain by a finer-textured clayey element. The centimetre- to decimetre-scale banding is formed of laminae in alternating, gradually intergrading sets of relatively coarse and relative fine-grained examples. At outcrop in the field, the banding is recognizable because the coarse sets prove to be recessive to varying degrees under the influence of weathering and current action. Independent evidence at two localities points toward an annual origin for the banding; at a third it arose during part of what appears to have been a relatively short period. Quantified physical arguments suggest that the textural banding is a response of suspended fine sediment to marked seasonal changes in sea temperature and windiness. The banded silts occur in four distinct stratigraphical contexts and record high deposition rates (order 0.01-0.1 m/yr). Because physical factors determine their textures, the silts potentially afford insights in all contexts into aspects of changing Holocene climatic conditions. In one context, the thickness of the bands points to high (order 0.01-0.1 m/yr) but comparatively short-lived (order 10s-100s yrs) rates of relative water-level rise. In the others, however, the banding has no implications for sea-level behaviour, and simply records gross environmental disequilibrium, for example, the recovery of mudflats/marshes after an erosional episode. Similarly, because on account of their rapid accumulation the banded silts preserve animal and human tracks and trackways especially well, they provide an archive of animal and human behaviour in the area during the Holocene.

Hydrological controls on soil redox dynamics in a peat-based, restored wetland

Increasing areas of altered wetland are being restored by re-flooding the soil. Evidence in the literature indicates that this practice can induce the redox-mediated release of soil nutrients, thereby increasing the risk of diffuse water pollution. However, for the sake of improving wedand management decisions, there is a need for more detailed studies of the underlying relationship between the hydrological and redox dynamics that explain this risk; this is particularly the case in agricultural peatlands that are commonly targeted for the creation of lowland wet grassland. A 12-month field study was conducted to evaluate the relationship between hydrological fluctuations and soil redox potential (Eh) in a nutrient-rich peat field (32 g N kg(-1) and 1100 mg P kg(-1) in the surface 0-30 cm soil) that had been restored as lowland wet grassland from intensive arable production. Field tensiometers were installed at the 30-, 60- and 90-cm soil depths, and Pt electrodes at the 10-, 30-, 60- and 90-cm depths, for daily logging of soil water tension and Eh, respectively. The values for soil water tension displayed a strong negative relationship (P < 0.001) with monthly dip well observations of water table height. Calculations of soil water potential from the logged tension values were used, therefore, to provide a detailed profile of field water level and, together with precipitation data, explained some of the variation in Eh. For example, during the summer, alternating periods of aerobism (Eh > 330 mV) in the surface, 0-10 cm layer of peat coincided with intense precipitation events. Redox potential throughout the 30-100 cm profile also fluctuated seasonally; indeed, at all depths Eh displayed a strong, negative relationship (P < 0.001) with water table height over the 12-month study period. However, Eh throughout the 30-100 cm profile remained relatively low (< 230 mV), indicating permanently reduced conditions that are associated with denitrification and reductive dissolution of Fe-bound P. The implications of these processes in the N- and P-rich peat for wetland plant diversity and water quality are discussed. (c) 2006 Elsevier B.V. All rights reserved.

The sediment trapping efficiency of the macro-tidal Daly Estuary, tropical Australia

Field studies were carried out on the water and sediment dynamics in the tropical, macro-tidal, Daly Estuary. The estuary is shallow, very-turbid, about 100 km long, and the entrance is funnel-shape. In the wet, high flow season, normal tidal ranges can be suppressed in the estuary, depending on inflow rates, and freshwater becomes dominant up to the mouth. At that time a fraction of the fine sediment load is exported offshore as a bottom-tagging nepheloid layer after the sediment falls out of suspension of the thin, near-surface, river plume. The remaining fraction and the riverine coarse sediment form a large sediment bar 10 km long, up to 6 m in height and extending across the whole width of the channel near the mouth. This bar, as well as shoals in the estuary, partially pond the mid- to upper-estuary. This bar builds up from the deposition of riverine sediment during a wet season with high runoff and can raise mean water level by up to 2 m in the upper estuary in the low flow season. This ponding effect takes about three successive dry years to disappear by the sediment forming the bar being redistributed all over the estuary by tidal pumping of fine and coarse sediment in the dry season, which is the low flow season. The swift reversal of the tidal currents from ebb to flood results in macro-turbulence that lasts about 20 min. Bed load transport is preferentially landward and occurs only for water currents greater than 0.6 m s(-1). This high value of the threshold velocity suggests that the sand may be cemented by the mud. The Daly Estuary thus is a leaky sediment trap with an efficiency varying both seasonally and inter-annually. (c) 2006 Elsevier Ltd. All rights reserved.

Sediment fingerprinting as an environmental forensics tool explaining cyanobacteria blooms in lakes

Cyanobacteria (blue-green algae) blooms in water bodies present serious public health issues with attendant economic and ecological impacts. Llyn Tegid (Lake Bala) is an important conservation and amenity asset within Snowdonia National Park, Wales which since the mid-1990s has experienced multiple toxic cyanobacteria blooms threatening the ecology and tourism-dependent local economy. Multiple working hypotheses explain the emergence of this problem, including climate change, land management linked to increased nutrient flux, hydromorphological alterations or changing trophic structure - any of which may operate individually or cumulatively to impair lake function. This paper reports the findings of a sedimentfingerprinting study using dated lake cores to explore the linkages between catchment and lake management practices and the emergence of the algal blooms problem. Since 1900 AD lake bed sedimentation rates have varied from 0.06 to 1.07 g cm−2 yr−1, with a pronounced acceleration since the early 1980s. Geochemical analysis revealed increases in the concentrations of total phosphorus (TP), calcium and heavy metals such as zinc and lead consistent with eutrophication and a rising pollution burden, particularly since the late 1970s. An uncertainty-inclusive sedimentfingerprinting approach was used to apportion the relative fluxes from the major catchment land cover types of improved pasture, rough grazing, forestry and channel banks. This showed improved pasture and channel banks are the dominant diffuse sources of sediment in the catchment, though forestry sources were important historically. Conversion of rough grazing to improved grassland, coupled with intensified land management and year-round livestock grazing, is concluded to provide the principal source of rising TP levels. Lake Habitat Survey and particle size analysis of lake cores demonstrate the hydromorphological impact of the River Dee Regulation Scheme, which controls water level and periodically diverts flow into Llyn Tegid from the adjacent Afon Tryweryn catchment. This hydromorphological impact has also been most pronounced since the late 1970s. It is concluded that an integrated approach combining land management to reduce agricultural runoff allied to improved water level regulation enabling recovery of littoral macrophytes offers the greatest chance halting the on-going cyanobacteria issue in Llyn Tegid.

Technical Note: The normal quantile transformation and its application in a flood forecasting system.

The Normal Quantile Transform (NQT) has been used in many hydrological and meteorological applications in order to make the Cumulated Distribution Function (CDF) of the observed, simulated and forecast river discharge, water level or precipitation data Gaussian. It is also the heart of the meta-Gaussian model for assessing the total predictive uncertainty of the Hydrological Uncertainty Processor (HUP) developed by Krzysztofowicz. In the field of geo-statistics this transformation is better known as the Normal-Score Transform. In this paper some possible problems caused by small sample sizes when applying the NQT in flood forecasting systems will be discussed and a novel way to solve the problem will be outlined by combining extreme value analysis and non-parametric regression methods. The method will be illustrated by examples of hydrological stream-flow forecasts.

Problems with binary pattern measures for flood model evaluation

As the calibration and evaluation of flood inundation models are a prerequisite for their successful application, there is a clear need to ensure that the performance measures that quantify how well models match the available observations are fit for purpose. This paper evaluates the binary pattern performance measures that are frequently used to compare flood inundation models with observations of flood extent. This evaluation considers whether these measures are able to calibrate and evaluate model predictions in a credible and consistent way, i.e. identifying the underlying model behaviour for a number of different purposes such as comparing models of floods of different magnitudes or on different catchments. Through theoretical examples, it is shown that the binary pattern measures are not consistent for floods of different sizes, such that for the same vertical error in water level, a model of a flood of large magnitude appears to perform better than a model of a smaller magnitude flood. Further, the commonly used Critical Success Index (usually referred to as F<2 >) is biased in favour of overprediction of the flood extent, and is also biased towards correctly predicting areas of the domain with smaller topographic gradients. Consequently, it is recommended that future studies consider carefully the implications of reporting conclusions using these performance measures. Additionally, future research should consider whether a more robust and consistent analysis could be achieved by using elevation comparison methods instead.

Urban microclimates and urban heat island in Chongqing, China

Chongqing is the largest directly-controlled municipality in China, which is now undergoing a rapid urbanization. The urbanization rate increased from 35.6% in 2000 to 48.3% in 2007, and it is estimated to reach at least 70% by 2020. The question remains open: What are the consequences of such rapid urbanization in Chongqing in terms of urban microclimate? Furthermore, Chongqing is located within the Three Gorges Reservoir (TGR) region and the upper Yangtze River, where the Three Gorges Reservoir (TGR) project started in 1993 and was completed in 2010. As one of the biggest construction projects in the world with a rising water level of 175m and water storage capacity of about 39.3 billion m3, it would be interesting to investigate how such a gigantic project impacts the surrounding micro-environment, especially in Chongqing. Different research approaches are adopted in the study. Our literature review indicates present studies on the urban climate in Chongqing are mainly confined within the historical trend analysis of several weather stations operated by the Chongqing government, little is known about the spatial distribution of urban air temperature and how the local land cover influences the air temperature, especially when there are rivers running through the Chongqing urban area. To contribute to the present knowledge, a series of field measurement campaigns and numerical simulations were carried out. Two complementary types of field measurements are included: fixed weather stations and mobile transverse measurement. Numerical simulations using a house-developed program are able to predict the urban air temperature in Chongqing.

A real-time measurement system for long-life flood monitoring and warning applications

A flood warning system incorporates telemetered rainfall and flow/water level data measured at various locations in the catchment area. Real-time accurate data collection is required for this use, and sensor networks improve the system capabilities. However, existing sensor nodes struggle to satisfy the hydrological requirements in terms of autonomy, sensor hardware compatibility, reliability and long-range communication. We describe the design and development of a real-time measurement system for flood monitoring, and its deployment in a flash-flood prone 650 km2 semiarid watershed in Southern Spain. A developed low-power and long-range communication device, so-called DatalogV1, provides automatic data gathering and reliable transmission. DatalogV1 incorporates self-monitoring for adapting measurement schedules for consumption management and to capture events of interest. Two tests are used to assess the success of the development. The results show an autonomous and robust monitoring system for long-term collection of water level data in many sparse locations during flood events.

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.

Efficient incorporation of channel cross-section geometry uncertainty into regional and global scale flood inundation models

This paper investigates the challenge of representing structural differences in river channel cross-section geometry for regional to global scale river hydraulic models and the effect this can have on simulations of wave dynamics. Classically, channel geometry is defined using data, yet at larger scales the necessary information and model structures do not exist to take this approach. We therefore propose a fundamentally different approach where the structural uncertainty in channel geometry is represented using a simple parameterization, which could then be estimated through calibration or data assimilation. This paper first outlines the development of a computationally efficient numerical scheme to represent generalised channel shapes using a single parameter, which is then validated using a simple straight channel test case and shown to predict wetted perimeter to within 2% for the channels tested. An application to the River Severn, UK is also presented, along with an analysis of model sensitivity to channel shape, depth and friction. The channel shape parameter was shown to improve model simulations of river level, particularly for more physically plausible channel roughness and depth parameter ranges. Calibrating channel Manning’s coefficient in a rectangular channel provided similar water level simulation accuracy in terms of Nash-Sutcliffe efficiency to a model where friction and shape or depth were calibrated. However, the calibrated Manning coefficient in the rectangular channel model was ~2/3 greater than the likely physically realistic value for this reach and this erroneously slowed wave propagation times through the reach by several hours. Therefore, for large scale models applied in data sparse areas, calibrating channel depth and/or shape may be preferable to assuming a rectangular geometry and calibrating friction alone.