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-correlated to investigate seasonal variations in groundwater-level response times, based on periods of 3-month duration. The time lags (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 lag 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.     

Estimation of water saturation dependence of dispersion in unsaturated porous media: experiments and modelling analysis

In the dispersion theory, a linear relationship has been verified between the coefficient of hydrodynamic dispersion and water velocity, both in saturated and in unsaturated porous media. But for unsaturated soils the variability of flow directions and microscopic velocities can be larger than in saturated soils because of the lower degree of water saturation. This leads to an increased dispersion. Therefore, relationships between water content and relative water velocity fluctuations and water content together with the coefficient of dispersivity in unsaturated porous media respectively have been investigated systematically by displacement experiments in glass beads and coarse-textured sandy soil columns. The breakthrough curves (BTCs) of chloride showed that an increase of solute mixing with a decrease of water content was caused by an increase of flow velocity fluctuations for different pathways. In order to explain the observed tailing effect in unsaturated flow, two mathematical models were used to fit theoretically derived nonlinear functions of water content dependent dispersivities for both porous media. The close agreement between the observed and computed results suggests that the theoretical model of hydrodynamic dispersion can be extended to transport in unsaturated porous media, providing that BTCs of the effluent water are used to estimate representative dispersivity parameters of soils.     

Predicting Vertical LNAPL Distribution in the Subsurface under the Fluctuating Water Table Effect

The present study proposes a methodology for predicting the vertical light nonaqueous-phase liquids (LNAPLs) distribution within an aquifer by considering the influence of water table fluctuations. The LNAPL distribution is predicted by combining (1) information on air/LNAPL and LNAPL/water interface elevations with (2) the initial elevation of the water table without LNAPL effect. Data used in the present study were collected during groundwater monitoring undertaken over a period of 4 months at a LNAPL-impacted observation well. In this study, the water table fluctuations raised the free LNAPL in the subsurface to an elevation of 206.63 m, while the lowest elevation was 205.70 m, forming a thickness of 0.93 m of LNAPL-impacted soil. Results show that the apparent LNAPL thickness in the observation well is found to be three times greater than the actual free LNAPL thickness in soil; a finding that agrees with previous studies reporting that apparent LNAPL thickness in observation wells typically exceeds the free LNAPL thickness within soil by a factor estimated to range between 2 and 10. The present study provides insights concerning the transient variation of LNAPL distribution within the subsurface and highlights the capability of the proposed methodology to mathematically predict the actual LNAPL thickness in the subsurface, without the need to conduct laborious field tests. Practitioners can use the proposed methodology to determine by how much the water table should be lowered, through pumping, to isolate the LNAPL-impacted soil within the unsaturated zone, which can then be subjected to in situ vadose zone remedial treatment.     

Quantification of vertical water fluxes in the vadose zone using particle‐size distribution and pedology‐based approaches to model soil heterogeneities

One‐dimensional flow simulations were conducted at four locations of the shallow alluvial aquifer of the upper Rhine River (at the Erstein polder) to quantify the time‐dependent moisture distribution, the water flux and the water volume infiltrated in the unsaturated zone as a function of soil heterogeneities during a five‐day‐long flooding event. Three methods of estimating the hydraulic parameters of soil in the vadose zone were tested. They are based on the following: (1) experimental data, (2) soil particle‐size distribution and (3) pedology information on soils. Water fluxes calculated from modelling approaches 2 and 3 were compared with those of the experiment‐based values and the effect of these differences on the arrival time and velocity of water at the water table were analysed. Major differences in water fluxes were found among the methods of estimating the hydrodynamic parameters. At the Terrace location, the groundwater recharge predicted using soil data from methods 1 and 2 are approximately 4500 and 2400 mm, respectively. Flow simulations using soil data and the experiment‐based method show the highest velocities of infiltrating water at the soil surface and largest volume of groundwater infiltration but result in the lowest centres of the moisture content mass. The results obtained using soil data based on the pedological method are similar to those calculated using soil parameters based on the particle‐size distribution of extracted soil samples. Water pressure profiles calculated on Terrace and Channel location, 3 and 7 days after the inundation event agreed reasonably well with those observed when using hydrodynamic parameters from the experiment‐based method. However, the flow model using the pedology‐based parameters largely underestimates the time needed to achieve hydrostatic conditions of the soil water profile once water flooding at the soil surface stops. This can be mainly attributed to the low values of estimated van Genuchten parameter α. Copyright © 2012 John Wiley & Sons, Ltd.     

Water uptake by saltcedar (Tamarix ramosissima) in a desert riparian forest: responses to intra‐annual water table fluctuation

There is considerable interest in naturalizing flow regime on managed rivers to slow the spread of saltcedar (Tamarix ramosissima) invasion in southwestern USA or to preserve riparian forests dominated by saltcedar and other species in northwestern China. However, little is known about the responses of established saltcedar in water sources to frequent intra‐annual fluctuation of water table resulting from this new, more dynamic flow regime. This study investigates how saltcedar at a riparian site in the middle reaches of the Heihe River, northwest China, responds in water sources use to intra‐annual water table fluctuations. Stable oxygen isotope was employed to determine accurate depth at which saltcedar obtains its water supply, and soil moisture monitoring was used to determine sources of plant‐available soil water. We found that the primary zone of water uptake by saltcedar were stable at 25–60 cm depth, but the water sources used by saltcedar switched between groundwater and soil moisture with the water table fluctuations. Saltcedar derived its water from groundwater when water table was at depth less than 60 cm but switched to soil moisture at 25–60 cm depth when water table declined. It is supposed that the well‐developed clay layer at 60–80 cm depth constrained lateral roots of saltcedar to the soil layers above 60 cm, while the fine‐textured soils at this site, which were periodically resaturated by rising groundwater before the stored soil moisture had become depleted, provided an important water reservoir for saltcedar when groundwater dropped below the primary zone of fine roots. The root distribution of saltcedar may also be related to local groundwater history. The quick decline in water table in the early 1980s when the riparian saltcedar had established may strand its roots in the shallow unsaturated zone. We suggested that raising the water table periodically instead of maintaining it invariably above the rooting depth could sustain desired facultative phreatophytes while maximizing water deliveries. Copyright © 2015 John Wiley & Sons, Ltd.     

The role of surface and sub-surface runoff processes in controlling cation export from a wetland watershed

Water and cation budgets were calculated for two sub-basins within a small low relief watershed in South-Central Ontario during a period of ephemeral runoff which was initiated by spring snow melt. The hydrology of one (upland) sub-basin was strongly influenced by seasonal fluctuations in the level of regional ground water. Saturated contributing areas formed in low lying regions adjacent to the stream channel where the water table rose to the surface, and stream discharge was a mixture of ground water and saturation overland flow. In the second sub-basin a wetland provided a large and spatially less variable saturated contributing area. Clay soils underlying the wetland resulted in a shallow perched water table, poorly drained and highly organic soils, and greatly reduced inputs of regional ground water. Stream discharge was largely the result of surface runoff from the wetland and adjacent areas of saturated soil.Inter-basin variations in water export were by far greater than variations in stream chemistry. As a result, inter-basin variations in cation export strongly reflected variations in water export over the time interval in which the majority of a given ion was lost from the watershed. Spatial differences in water export were least at the onset of runoff when basin saturation was greatest and overland flow made large contributions to the discharge from both sub-basins. Potassium and hydrogen had high concentrations at this time which caused these ions to show only small spatial differences in export. With decreases in the areal extent of soil saturation, and increases in the storage capacity of the wetland, the hydrologic contrast between sub-basins increased. Greater water loss from the upland area resulted from a greater discharge of regional ground water, and a more rapid expansion of the saturated contributing areas during storm events. Calcium, magnesium, and sodium concentrations increased steadily during the first 3 weeks of runoff, so that the peak export of these cations occurred later in the runoff period at times of higher concentration, but lower and spatially more variable discharges. Consequently, spatial differences in the loss of these ions was great and favoured the upland sub-basin, since the majority of export occurred when the hydrologic contrast between sub-basins was largest.     

Groundwater Recharge at Five Representative Sites in the Hebei Plain,China

Accurate estimates of groundwater recharge are essential for effective management of groundwater, especially when supplies are limited such as in many arid and semiarid areas. In the Hebei Plain, China, water shortage is increasingly restricting socioeconomic development, especially for agriculture, which heavily relies on groundwater. Human activities have greatly changed groundwater recharge there during the past several decades. To obtain better estimates of recharge in the plain, five representative sites were selected to investigate the effects of irrigation and water table depth on groundwater recharge. At each site, a one‐dimensional unsaturated flow model (Hydrus‐1D) was calibrated using field data of climate, soil moisture, and groundwater levels. A sensitivity analysis of evapotranspirative fluxes and various soil hydraulic parameters confirmed that fine‐textured surface soils generally generate less recharge. Model calculations showed that recharge on average is about 175 mm/year in the piedmont plain to the west, and 133 mm/year in both the central alluvial and lacustrine plains and the coastal plain to the east. Temporal and spatial variations in the recharge processes were significant in response to rainfall and irrigation. Peak time‐lags between infiltration (rainfall plus irrigation) and recharge were 18 to 35 days in the piedmont plain and 3 to 5 days in the central alluvial and lacustrine plains, but only 1 or 2 days in the coastal plain. This implies that different time‐lags corresponding to different water table depths must be considered when estimating or modeling groundwater recharge.     

Evaluation of the soil water balance in an alluvial flood plain with a shallow groundwater table

Abstract

Many of the hydrological and ecological functions of alluvial flood plains within watersheds depend on the water flow exchanges between the vadoze soil zone and the shallow groundwater. The water balance of the soil in the flood plain is investigated, in order to evaluate the main hydrological processes that underlie the temporal dynamics of soil moisture and groundwater levels. The soil moisture and the groundwater level in the flood plain were monitored continuously for a three-year period. These data were integrated with the results derived from applying a physically-based numerical model which simulated the variably-saturated vertical water flow in the soil. The analysis indicated that the simultaneous processes of lateral groundwater flow and the vertical recharge from the unsaturated zone caused the observed water table fluctuations. The importance of these flows in determining the rises in the water table varied, depending on soil moisture and groundwater depth before precipitation. The monitoring period included two hydrological years (September 2009–September 2011). About 13% of the precipitation vertically recharged the groundwater in the first year and about 50% in the second. The difference in the two recharge coefficients was in part due to the lower groundwater levels in the recharge season of the first hydrological year, compared to those observed in the second. In the latter year, the shallow groundwater increased the soil moisture in the unsaturated zone due to capillary rise, and so the mean hydraulic conductivity of the unsaturated soil was high. This moisture state of soil favoured a more efficient conversion of infiltrated precipitation into vertical groundwater recharge. The results show that groundwater dynamics in the flood plain are an important source of temporal variability in soil moisture and vertical recharge processes, and this variability must be properly taken into account when the water balance is investigated in shallow groundwater environments.

Citation Pirastru, M. and Niedda, M., 2013. Evaluation of the soil water balance in an alluvial flood plain with a shallow groundwater table. Hydrological Sciences Journal, 58 (4), 898–911.     

Measurements of the downslope flow of water in a soil

Measurements have been made of unsaturated and saturated lateral soil water flow on a convex hill-slope with a good soil cover and impermeable bedrock during natural rainstorms. The hydraulics of flow are examined in detail with particular reference to the role of breaks in vertical permeability, the change from saturated to unsaturated flow and the velocity of flow. In this instance, after rainfall slope flow is dominated by vertical unsaturated movement towards the profile base. Preceding upslope moisture gradients result in the growth of a zone of soil saturation upwards from the slope base. Slope discharge, through the B and B/C horizons, is related to the form of the saturation zone, within which flow is lateral, according to Darcy's law. The time required for vertical percolation and the low hydraulic conductivity of the lower soil horizons result in a hillslope hydrograph which is delayed and attenuated and cannot be regarded as stormflow. During drainage the saturation zone contracts and is replaced by a lateral unsaturated flow system at the profile base which supplies discharge from the B/C horizon for up to 42 days without further recharge. It is concluded that, in general, either distinct soil horizons or impermeable bedrock are essential for the initiation of lateral flow. Saturated flow is likely to dominate hillslope hydrographs through non-capillary pore spaces but these may be integrated to the point where Darcy's law still holds. Although lateral soil water flow must be a widespread phenomenon, it is unlikely to provide storm runoff to the stream unless saturated conditions are generated within the organic horizons for flow within the lower soil horizons is dominated by non-Darcian flow through non-capillary spaces in the soil.     

Assessment of evaporation and water fluxes in a column of dry saline soil subject to different water table levels

Dry saline soils are common in the arid and hyper‐arid basins located in the Chilean Altiplano, where evaporation from shallow groundwater is typically the major component of the water balance. Thus, a good understanding of evaporation processes is necessary for improving water resource planning and management in these regions. In this study, we conducted laboratory experiments with a natural saline soil column to estimate evaporation rates and assess the liquid and water vapor fluxes under different water table levels. Water content, electrical conductivity and temperature at different depths were utilized to assess the liquid and water vapor fluxes in the soil column. We observed movement of water that dissolves salts from the soil and transports them to areas in the column where they accumulate. Isothermal liquid flux was predominant, while thermal and isothermal liquid and thermal water vapor fluxes were negligible, except for deep water table levels where isothermal and thermal water vapor fluxes had similar magnitude but opposite directions. Differences observed in total fluxes for all water table levels were due to different upward and downward fluxes, which depend on changes in water content and temperature within the soil profile. Both the vapor flux magnitude and direction were found to be very sensitive to the choice of empirical parameters used in flux quantification, such as tortuosity and the enhancement factor for local temperature gradients in the air phase within the column. Copyright © 2013 John Wiley & Sons, Ltd.     

A hydropedological approach to quantifying groundwater recharge in various glacial settings of the mid‐continental USA

Data from six monitoring stations were combined with a soil‐water dynamics model (HYDRUS 1D) to achieve physically‐based estimates of shallow water‐table recharge in representative hydropedological settings of the glaciated midwestern U.S.A. Calibration involved inverse modeling that yielded optimized hydraulic parameters. Root mean square errors for modeled versus measured soil moisture contents were generally within 3% for all soil layers at the six study sites. The optimized models also accurately simulated recharge values that corresponded to observed water‐table fluctuations. Optimized parameter values were consistent with estimates from a pedotransfer function, lab analyses, and field experiments. Forward modeling indicated that shallow water‐table recharge in mid‐continent glacial settings is approximately 35% of precipitation, but interannual and monthly variability is significant. Soil parent materials and horizon characteristics influence recharge primarily through their control on K_s with clay‐rich till parent materials producing values as low as 16% and coarse‐grained outwash parent materials producing values as high as 58% of precipitation. During the three‐year study period, distinct seasonality of recharge was observed with most recharge occurring in the winter (seasonal mean of all sites was 66% of precipitation) and lesser but interannually stable amounts in the spring (44%), summer (13%), and autumn (16%). This research underscores the importance of incorporating pedological information into models of soil‐water dynamics and groundwater recharge. © 2015 The Authors. Hydrological Processes published by John Wiley & Sons Ltd.     

Water and solute movement in a coarse-textured water-repellent field soil

Unstable water flow in water-repellent unsaturated soils can significantly affect the processes of infiltration and soil water redistribution. A field experiment was carried out to study the effect of water-repellency on water and bromide movement in a coarse-textured soil in the southwestern part of The Netherlands. The field data were analyzed using a relatively simple numerical model based on the standard Richards' equation for unsaturated water flow and the Fickian-based convection-dispersion equation for solute transport. Water-repellency was accounted for by multiplying the water content and the unsaturated hydraulic conductivity of the soil with F, a factor equal to the volumetric fraction of soil occupied by preferential flow paths resulting from the unstable flow process. The good comparison of simulated and measured bromide concentrations suggests that the model provides a viable method for simulating unstable water flow in water-repellent soils.     

Water fluxes at a fluctuating water table and groundwater contributions to wheat water use in the lower Yellow River flood plain,China

Capillary upflow from and deep percolation to a water table may be important in crop water supply in irrigated areas of the lower Yellow River flood plain, north China. These fluxes at the water table and the variations of the capillary upflow in relation to crop evapotranspiration need to be investigated to quantify the effect of a water table on soil water balance and to improve agricultural water management. A large weighing lysimeter was used to determine daily crop evapotranspiration, daily capillary upflow from and daily percolation to a fluctuating water table during a rotation period with wheat growing in a dry season and maize in a rainy season. The water table depth varied in the range 0·7–2·3 m during the maize growth period and 1·6–2·4 m during the wheat growth period. Experimental results showed that the capillary upflow and the percolation were significant components of the soil water balance. Three distinctly different phases for the water fluxes at the water table were observed through the rotation period: water downward period, the period of no or small water fluxes, and water upward period. It implied that the temporal pattern of these water fluxes at the water table was intimately associated with the temporal distribution of rainfall through the rotation period. An empirical equation was determined to estimate the capillary upflow in relation to wheat evapotranspiration and root zone soil water content for local irrigation scheduling. Coupled with the FAO‐Penman–Monteith equation, the equation offers a fast and low cost solution to assess the effect of capillary upflow from a water table on wheat water use. Copyright © 2007 John Wiley & Sons, Ltd.     

Effects of Airflow Induced by Rainfall on Shallow Groundwater Table Fluctuations

An investigation of groundwater table fluctuations induced by rainfall should consider interactions between the liquid and gas phases in soils. In this study, a water‐air two‐phase flow model was initially verified by simulating an infiltration experiment. It was then employed to model the interactions between liquid and gas phases regarding actions of airflow on the groundwater table and the fluctuations of the phreatic level and water level in the well induced by rainfall. The effects of airflo7w caused by rainfall on phreatic level fluctuations were also studied quantitatively by comparing the results obtained using the proposed model with those obtained from a water single‐phase flow model. The simulation results show that in addition to actual recharge, compressed airflow in unsaturated zones causes the phreatic level to increase, but the rise in the phreatic level is lower than that in the pore‐air pressure head in unsaturated zones due to the mitigation of capillary fringe. The existence of airflow enhances the phreatic level rise during and after rainfall. In addition, the water level in the well, pushed by the phreatic level fluctuations, varies similarly to the phreatic level, but it experiences somewhat delayed and slightly attenuated. The Lisse effect precisely reflects the phreatic level fluctuations before actual recharge. Furthermore, the fluctuations in the phreatic level and water level in the well and the contributions of airflow to phreatic level fluctuations are affected by many factors: rain intensity, initial moisture, overlying aquitard, groundwater table depths, and screen depths of the well.     

Recharge processes during snowmelt: An isotopic and hydrometric investigation

Results from hydrometric and isotopic investigations of unsaturated flow during snowmelt are presented for a hillslope underlain by well-sorted sands. Passage of melt and rainwater through the vadose zone was detected from temporal changes in soil water ²H concentrations obtained from sequential soil cores. Bypassing flow was indicated during the initial snowmelt phase, but was confined to the near-surface zone. Recharge below this zone was via translatory flow, as meltwater inputs displaced premelt soil water. Estimates of premelt water fluxes indicate that up to 19 per cent of the premelt soil water may have been immobile. Average water particle velocities during snowmelt ranged from 6.2 × 10^?7 to 1.1 × 10^?6 ms^?1, suggesting that direct groundwater recharge by meltwater during snowmelt was confined to areas where the premelt water table was within 1 m of the ground surface. Soil water ²H signatures showed a rapid response to isotopically-heavy rain-on-snow inputs late in the melt. In addition, spatial variations in soil moisture content at a given depth induced a pronounced lateral component to the predominantly vertical transport of water. Both factors may complicate isotopic profiles in the vadose zone, and should be considered when employing environmental isotopes to infer recharge processes during snowmelt.     

Water transport and tracer mixing in volcanic ash soils at a tropical hillslope: A wet layered sloping sponge

Andosol soils formed in volcanic ash provide key hydrological services in montane environments. To unravel the subsurface water transport and tracer mixing in these soils we conducted a detailed characterization of soil properties and analyzed a 3-year data set of sub-hourly hydrometric and weekly stable isotope data collected at three locations along a steep hillslope. A weakly developed (52–61 cm depth), highly organic andic (Ah) horizon overlaying a mineral (C) horizon was identified, both showing relatively similar properties and subsurface flow dynamics along the hillslope. Soil moisture observations in the Ah horizon showed a fast responding (few hours) “rooted” layer to a depth of 15 cm, overlying a “perched” layer that remained near saturated year-round. The formation of the latter results from the high organic matter (33–42%) and clay (29–31%) content of the Ah horizon and an abrupt hydraulic conductivity reduction in this layer with respect to the rooted layer above. Isotopic signatures revealed that water resides within this soil horizon for short periods, both at the rooted (2 weeks) and perched (4 weeks) layer. A fast soil moisture reaction during rainfall events was also observed in the C horizon, with response times similar to those in the rooted layer. These results indicate that despite the perched layer, which helps sustain the water storage of the soil, a fast vertical mobilization of water through the entire soil profile occurs during rainfall events. The latter being the result of the fast transmissivity of hydraulic potentials through the porous matrix of the Andosols, as evidenced by the exponential shape of the water retention curves of the subsequent horizons. These findings demonstrate that the hydrological behavior of volcanic ash soils resembles that of a “layered sponge,” in which vertical flow paths dominate.     

Changes in water table level influence solute transport in uniform porous media

Changes in the water table level result in variable water saturation and variable hydrological fluxes at the interface between the unsaturated and saturated zone. This may influence the transport and fate of contaminants in the subsurface. The objective of this study was to examine the impact of a decreasing and an increasing water table on solute transport. We conducted tracer experiments at downward flow conditions in laboratory columns filled with two different uniform porous media under static and transient flow conditions either increasing or decreasing the water table. Tracer breakthrough curves were simulated using a mobile–immobile transport model. The resulting transport parameters were compared to identify dominant transport processes. Changes in the water table level affected dispersivities and mobile water fractions depending on the direction of water table movement and the grain size of the porous media. In fine glass beads, the water flow velocity was similar to the decline rate of the water table, and the mobile water fraction was decreased compared with steady‐state saturated conditions. However, immobile water was negligible. In coarse glass beads, water flow was faster because of fingered flow in the unsaturated part, and the mobile water fraction was smaller than in the fine material. Here, a rising water table led to an even smaller mobile water fraction and increased solute spreading because of diffusive interaction with immobile water. We conclude that changes of the water table need to be considered to correctly simulate transport in the subsurface at the transition of the unsaturated–saturated zone. Copyright © 2014 John Wiley & Sons, Ltd.     

Prediction of the arrival of peak nitrate concentrations at the water table at the regional scale in Great Britain

A simple process‐based approach to predict regional‐scale loading of nitrate at the water table was implemented in a GIS for Great Britain. This links a nitrate input function, unsaturated zone thickness, and lithologically dependent rate of nitrate unsaturated zone travel to estimate arrival time of nitrate at the water table. The nitrate input function is the loading at the base of the soil and has been validated using unsaturated zone porewater profiles. The unsaturated zone thickness uses groundwater levels based on regional‐scale observations infilled by interpolated river base levels. Estimates of the rate of unsaturated zone travel are attributed from regional‐scale hydrogeological mapping. The results indicate that peak nitrate loading may have already arrived at the water table for many aquifers, but that it has not where the unsaturated zone is relatively thick There are contrasting outcomes for the two main aquifers which have similar unsaturated zone velocities, the predominantly low relief Permo‐Triassic sandstones, and the Chalk, which forms significant topographic features. For about 60% of the Chalk, the peak input has not yet reached the water table and will continue to arrive over the next 60 years. The methodology is readily transferable and provides a robust method for estimating peak arrival time for any diffuse conservative pollutant where an input function can be defined at a regional scale and requires only depth to groundwater and a hydrogeological classification. The methodology is extendable in that if additional information is available this can easily be incorporated into the model scheme. British Geology Survey © NERC 2011. Hydrological Process © 2011 John Wiley & Sons, Ltd