Catchment-scale modelling of flow and nutrient transport in the Chalk unsaturated zone

The unsaturated zone exerts a major control on the delivery of nutrients to Chalk streams, yet flow and transport processes in this complex, dual-porosity medium have remained controversial. A major challenge arises in characterising these processes, both at the detailed mechanistic level and at an appropriate level for inclusion within catchment-scale models for nutrient management. The lowland catchment research (LOCAR) programme in the UK has provided a unique set of comprehensively instrumented groundwater-dominated catchments. Of these, the Pang and Lambourn, tributaries of the Thames near Reading, have been a particular focus for research into subsurface processes and surface water-groundwater interactions. Data from LOCAR and other sources, along with a new dual permeability numerical model of the Chalk, have been used to explore the relative roles of matrix and fracture flow within the unsaturated zone and resolve conflicting hypotheses of response. From the improved understanding gained through these explorations, a parsimonious conceptualisation of the general response of flow and transport within the Chalk unsaturated zone was formulated. This paper summarises the modelling and data findings of these explorations, and describes the integration of the new simplified unsaturated zone representation with a catchment-scale model of nutrients (INCA), resulting in a new model for catchment-scale flow and transport within Chalk systems: INCA-Chalk. This model is applied to the Lambourn, and results, including hindcast and forecast simulations, are presented. These clearly illustrate the decadal time-scales that need to be considered in the context of nutrient management and the EU Water Framework Directive. (C) 2007 Elsevier B.V. All rights reserved.

Vinasse application to sugar cane fields. Effect on the unsaturated zone and groundwater at Valle del Cauca (Colombia)

Extensive application of vinasse, a subproduct from sugar cane plantations for bioethanol production, is currently taking place as a source of nutrients that forms part of agricultural management in different agroclimatic regions. Liquid vinasse composition is characterised by high variability of organic compounds and major ions, acid pH (4.7), high TDS concentration (117,416–599,400 mg L− 1) and elevated EC (14,350–64,099 μS cm− 1). A large-scale sugar cane field application is taking place in Valle del Cauca (Colombia), where monitoring of soil, unsaturated zone and the aquifer underneath has been made since 2006 to evaluate possible impacts on three experimental plots. For this assessment, monitoring wells and piezometers were installed to determine groundwater flow and water samples were collected for chemical analysis. In the unsaturated zone, tensiometers were installed at different depths to determine flow patterns, while suction lysimeters were used for water sample chemical determinations. The findings show that in the sandy loam plot (Hacienda Real), the unsaturated zone is characterised by low water retention, showing a high transport capacity, while the other two plots of silty composition presented temporal saturation due to La Niña event (2010–2011). The strong La Niña effect on aquifer recharge which would dilute the infiltrated water during the monitoring period and, on the other hand dissolution of possible precipitated salts bringing them back into solution may occur. A slight increase in the concentration of major ions was observed in groundwater (~ 5% of TDS), which can be attributed to a combination of factors: vinasse dilution produced by water input and hydrochemical processes along with nutrient removal produced by sugar cane uptake. This fact may make the aquifer vulnerable to contamination.

Stratigraphy of the unsaturated zone and uppermost part of the Snake River Plain aquifer at Test Area North, Idaho National Engineering Laboratory, Idaho /

Mode of access: Internet.

Zum Transport von partikelgebundenen Schadstoffen in der ungesättigten Zone

Bei einer Risikoabschätzung bezüglich einer Gefährdung des Schutzgutes Grundwasser müssen alle relevanten Transportpfade, auf denen Schadstoffe durch die Bodenzone bis ins Grundwasser verlagert werden, identifiziert und quantifiziert werden. Die Verlagerung von Schadstoffen gebunden an mobile Partikel im Sickerwasser wird dabei oft vernachlässigt. In dieser Arbeit wurden sowohl experimentelle Untersuchungen zum Partikeltransport in der Bodenzone als auch Szenarienmodellierungen hinsichtlich der Wechselwirkung Partikel/Schadstoff durchgeführt. Die experimentellen ungesättigten Säulenversuche wurden unter naturnahen stationären und instationären hydraulischen und hydrochemischen Bedingungen durchgeführt. Dabei wurde der Einfluss der Parameter Durchmesser Bodenmatrix, Partikelgröße, Beregnungsintensität, Oberflächenspannung und Hydrochemie auf den Transport von natürlichen und synthetischen Partikeln untersucht. Des Weiteren wurden Untersuchungen zur partikelgebundenen Verlagerung von Phenanthren durchgeführt. In einer numerischen Szenarienmodellierung mit dem Modell SMART wurde untersucht, unter welchen Randbedingungen der Transport von Partikeln gleichzeitig zu signifikanten partikelgebundenen Schadstoffkonzentrationen im Grundwasser führt. Dabei wurden die Parameter Lithologie Partikel/Boden, Hydrophobizität Schadstoff, Partikelkonzentration, Partikeldurchmesser sowie Körnung Bodenmatrix variiert. Die Ergebnisse dieser Arbeit zeigen, dass der partikelgebundene Schadstofftransportpfad in der ungesättigten Bodenzone in verschiedenen Szenarien den Anteil mobiler Schadstoffe, die mit dem Sickerwasser ins Grundwasser gelangen, signifikant erhöht. Auf Basis der experimentellen und theoretischen Untersuchungen wurde ein zweistufiges Bewertungsschema entwickelt, das bereits im Vorfeld einer Risikoabschätzung als Entscheidungshilfe hinsichtlich der Relevanz einer Mobilisierung, eines Transports und des Rückhalts von partikelgebundenen Schadstoffen in der ungesättigten Zone dient.

Measurements of a water potential and water content in unsaturated crystalline rock

A water desaturation zone develops around a tunnel in water-saturated rock when the evaporative water loss at the rock surface is larger than the water flow from the surrounding saturated region of restricted permeability. We describe the methods with which such water desaturation processes in rock materials can be quantified. The water retention characteristic theta(psi) of crystalline rock samples was determined with a pressure membrane apparatus. The negative water potential, identical to the capillary pressure, psi, below the tensiometric range (psi < -0.1 MPa) can be measured with thermocouple psychrometers (TP), and the volumetric water contents, theta, by means of time domain reflectometry (TDR). These standard methods were adapted for measuring the water status in a macroscopically unfissured granodiorite with a total porosity of approximately 0.01. The measured water retention curve of granodiorite samples from the Grimsel test site (central Switzerland) exhibits a shape which is typical for bimodal pore size distributions. The measured bimodality is probably an artifact of a large surface ratio of solid/voids. The thermocouples were installed without a metallic screen using the cavity drilled into the granodiorite as a measuring chamber. The water potentials observed in a cylindrical granodiorite monolith ranged between -0.1 and -3.0 MPa; those near the wall in a ventilated tunnel between -0.1 and -2.2 MPa. Two types of three-rod TDR Probes were used, one as a depth probe inserted into the rock, the other as a surface probe using three copper stripes attached to the surface for detecting water content changes in the rock-to-air boundary. The TDR signal was smoothed with a low-pass filter, and the signal length determined based on the first derivative of the trace. Despite the low porosity of crystalline rock these standard methods are applicable to describe the unsaturated zone in solid rock and may also be used in other consolidated materials such as concrete.

Understanding the Aquifier : 3D Visualisation and Groundwater Time-Series Animation in Lockyer Valley

The Lockyer Valley, southeast Queensland, hosts intensive irrigated agriculture using groundwater from over 5000 alluvial bores. A current project is considering introduction of PRW (purified recycled water) to augment groundwater supplies. To assess this, a valley-wide MODFLOW simulation model is being developed plus a new unsaturated zone flow model. To underpin these models and provide a realistic understanding of the aquifer framework a 3D visualisation model has been developed using Groundwater Visualisation System (GVS) software produced at QUT.

Contribution of geophysical surveys to groundwater modelling of a porous aquifer in semiarid Niger: An overview

Subsurface geophysical surveys were carried out using a large range of methods in an unconfined sandstone aquifer in semiarid south-western Niger for improving both the conceptual model of water flow through the unsaturated zone and the parameterization of numerical a groundwater model of the aquifer. Methods included: electromagnetic mapping, electrical resistivity tomography (ERT), resistivity logging, time domain electromagnetic sounding (TDEM), and magnetic resonance sounding (MRS). Analyses of electrical conductivities, complemented by geochemical measurements, allowed us to identify preferential pathways for infiltration and drainage beneath gullies and alluvial fans. The mean water content estimated by MRS (13%) was used for computing the regional groundwater recharge from long-term change in the water table. The ranges in permeability and water content obtained with MRS allowed a reduction of the degree of freedom of aquifer parameters used in groundwater modelling.

Estimation of soil hydraulic properties and their uncertainty: comparison between laboratory and field experiment

Often the soil hydraulic parameters are obtained by the inversion of measured data (e.g. soil moisture, pressure head, and cumulative infiltration, etc.). However, the inverse problem in unsaturated zone is ill-posed due to various reasons, and hence the parameters become non-unique. The presence of multiple soil layers brings the additional complexities in the inverse modelling. The generalized likelihood uncertainty estimate (GLUE) is a useful approach to estimate the parameters and their uncertainty when dealing with soil moisture dynamics which is a highly non-linear problem. Because the estimated parameters depend on the modelling scale, inverse modelling carried out on laboratory data and field data may provide independent estimates. The objective of this paper is to compare the parameters and their uncertainty estimated through experiments in the laboratory and in the field and to assess which of the soil hydraulic parameters are independent of the experiment. The first two layers in the field site are characterized by Loamy sand and Loamy. The mean soil moisture and pressure head at three depths are measured with an interval of half hour for a period of 1 week using the evaporation method for the laboratory experiment, whereas soil moisture at three different depths (60, 110, and 200 cm) is measured with an interval of 1 h for 2 years for the field experiment. A one-dimensional soil moisture model on the basis of the finite difference method was used. The calibration and validation are approximately for 1 year each. The model performance was found to be good with root mean square error (RMSE) varying from 2 to 4 cm(3) cm(-3). It is found from the two experiments that mean and uncertainty in the saturated soil moisture (theta(s)) and shape parameter (n) of van Genuchten equations are similar for both the soil types. Copyright (C) 2010 John Wiley & Sons, Ltd.

Implications of using On-Farm Flood Flow Capture to recharge groundwater and mitigate flood risks along the Kings River, CA

Two large hydrologic issues face the Kings Basin, severe and chronic overdraft of about 0.16M ac-ft annually, and flood risks along the Kings River and the downstream San Joaquin River. Since 1983, these floods have caused over $1B in damage in today’s dollars. Capturing flood flows of sufficient volume could help address these two pressing issues which are relevant to many regions of the Central Valley and will only be exacerbated with climate change. However, the Kings River has high variability associated with flow magnitudes which suggests that standard engineering approaches and acquisition of sufficient acreage through purchase and easements to capture and recharge flood waters would not be cost effective. An alternative approach investigated in this study, termed On-Farm Flood Flow Capture, involved leveraging large areas of private farmland to capture flood flows for both direct and in lieu recharge. This study investigated the technical and logistical feasibility of best management practices (BMPs) associated with On-Farm Flood Flow Capture. The investigation was conducted near Helm, CA, about 20 miles west of Fresno, CA. The experimental design identified a coordinated plan to determine infiltration rates for different soil series and different crops; develop a water budget for water applied throughout the program and estimate direct and in lieu recharge; provide a preliminary assessment of potential water quality impacts; assess logistical issues associated with implementation; and provide an economic summary of the program. At check locations, we measured average infiltration rates of 4.2 in/d for all fields and noted that infiltration rates decreased asymptotically over time to about 2 – 2.5 in/d. Rates did not differ significantly between the different crops and soils tested, but were found to be about an order of magnitude higher in one field. At a 2.5 in/d infiltration rate, 100 acres are required to infiltrate 10 CFS of captured flood flows. Water quality of applied flood flows from the Kings River had concentrations of COC (constituents of concern; i.e. nitrate, electrical conductivity or EC, phosphate, ammonium, total dissolved solids or TDS) one order of magnitude or more lower than for pumped groundwater at Terranova Ranch and similarly for a broader survey of regional groundwater. Applied flood flows flushed the root zone and upper vadose zone of nitrate and salts, leading to much lower EC and nitrate concentrations to a depth of 8 feet when compared to fields in which more limited flood flows were applied or for which drip irrigation with groundwater was the sole water source. In demonstrating this technology on the farm, approximately 3,100 ac-ft was diverted, primarily from April through mid-July, with about 70% towards in lieu and 30% towards direct recharge. Substantial flood flow volumes were applied to alfalfa, wine grapes and pistachio fields. A subset of those fields, primarily wine grapes and pistachios, were used primarily to demonstrate direct recharge. For those fields about 50 – 75% of water applied was calculated going to direct recharge. Data from the check studies suggests more flood flows could have been applied and infiltrated, effectively driving up the amount of water towards direct recharge. Costs to capture flood flows for in lieu and direct recharge for this project were low compared to recharge costs for other nearby systems and in comparison to irrigating with groundwater. Moreover, the potentially high flood capture capacity of this project suggests significant flood avoidance costs savings to downstream communities along the Kings and San Joaquin Rivers. Our analyses for Terranova Ranch suggest that allocating 25% or more flood flow water towards in lieu recharge and the rest toward direct recharge will result in an economically sustainable recharge approach paid through savings from reduced groundwater pumping. Two important issues need further consideration. First, these practices are likely to leach legacy salts and nitrates from the unsaturated zone into groundwater. We develop a conceptual model of EC movement through the unsaturated zone and estimated through mass balance calculations that approximately 10 kilograms per square meter of salts will be flushed into the groundwater through displacing 12 cubic meters per square meter of unsaturated zone pore water. This flux would increase groundwater salinity but an equivalent amount of water added subsequently is predicted as needed to return to current groundwater salinity levels. All subsequent flood flow capture and recharge is expected to further decrease groundwater salinity levels. Second, the project identified important farm-scale logistical issues including irrigator training; developing cropping plans to integrate farming and recharge activities; upgrading conveyance; and quantifying results. Regional logistical issues also exist related to conveyance, integration with agricultural management, economics, required acreage and Operation and Maintenance (O&M).

Controlling nitrate pollution: an integrated approach

The introduction of the EU Water Framework Directive requires policy to address non-point source pollution as part of an overall integrated strategy to improve the ecological status of water bodies. In this paper, we combine an economic optimisation framework with a dynamic simulation model of N transport in the Kennet Catchment to link decisions taken at the farm level to reductions in nitrate concentrations in the River Kennet. We examine a variety of policies targeted at reducing fertiliser use and changing the way in which farm land is used. We find that a tax on nitrogen emerges as the best policy both in terms of cost- and environmental effectiveness. Such a policy involves a considerable reduction in fertiliser use, as well as, a restructuring of land-use away from arable towards increased use of set-aside. Budgetary implications of such a radical move towards set-aside would be huge and hence unlikely to be politically palatable given the objective of reducing the EU budgetary allocation to agriculture. Additionally, the current rise in world demand for food may also mitigate calls for increasing the proportion of land taken out of agricultural production. Although the study succeeds in establishing a link between actions on the farm and nitrate concentrations in the stream water, further work is required to explore the effect of the retention of nitrates in the unsaturated zone and groundwater on this link.

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.