Numerical simulation of groundwater flow patterns using flux as upper boundary

Since the 1960s, most of the studies on groundwater flow systems by analytical and numerical modelling have been based on given‐head upper boundaries. The disadvantage of the given‐head approach is that the recharge into and discharge from a basin vary with changes in hydraulic conductivity and/or basin geometry. Consequently, flow patterns simulated with given‐head boundaries but with different hydraulic conductivities and/or basin geometry may not reflect the effects of these variables. We conducted, therefore, numerical simulations of groundwater flow in theoretical drainage basins using flux as the upper boundary and realistically positioned fluid‐potential sinks while changing the infiltration intensity, hydraulic conductivities, and geometric configuration of the basin. The simulated results demonstrate that these variables are dominant factors controlling the flow pattern in a laterally closed drainage basin. The ratio of infiltration intensity to hydraulic conductivity (R_ic) has been shown to be an integrated pattern‐parameter in a basin with a given geometric configuration and possible fluid‐potential‐sink distribution. Successively, the changes in flow patterns induced by stepwise reductions in R_ic are identical, regardless of whether the reductions are due to a decrease in infiltration intensity or an increase in hydraulic conductivity. The calculated examples show five sequential flow patterns containing (i) only local, (ii) local–intermediate, (iii) local–intermediate–regional, (iv) local–regional, and (v) just regional flow systems. The R_ic was found to determine also whether a particular sink is active or not as a site of discharge. Flux upper boundary is preferable for numerical simulation when discussing the flow patterns affected by a change of infiltration, the hydraulic conductivity, or the geometry of a basin. Copyright © 2012 John Wiley & Sons, Ltd.     

Stochastic analysis of flux and head moments in a heterogeneous aquifer system

This study investigates the behavior of flux and head in a strongly heterogeneous three-dimensional aquifer system. The analyses relied on data from 520 slug tests together with 38,000 one-foot core intervals lithological data from the site of the General Separations Area in central Savannah River Site, South Carolina, USA. The skewness in the hydraulic conductivity histograms supported the geologic information for the top two aquifers, but revealed stronger clay content, than was reported for the bottom aquifer. The log-normal distribution model described adequately the hydraulic conductivity measurements for all three aquifers although, other distributions described equally well the bottom aquifer measurements. No apparent anisotropy on the horizontal plane was found for the three aquifers, but ratios of horizontal to vertical correlation lengths between 33 and 75 indicated a strong stratification at the site. Three-dimensional Monte Carlo stochastic simulations utilized a grid with larger elements than the support volume of measurements, but of sub-REV (representative elementary volume) dimensions. This necessitated, on one hand, the use of upscaled hydraulic conductivity expressions, but on the other hand did not allow for the use of anisotropic effective hydraulic conductivity expressions (Sarris and Paleologos in J Stoch Environ Res Risk Assess 18: 188–197, 2004). Flux mean and standard deviations components were evaluated on three vertical cross-sections. The mean and variance of the horizontal flux component normal to a no-flow boundary tended to zero at approximately two to three integral scales from that boundary. Close to a prescribed head boundary both the mean and variance of the horizontal flux component normal to the boundary increased from a stable value attained at a distance of about five integral scales from that boundary. The velocity field 〈q_x〉 was found to be mildly anisotropic in the top two aquifers, becoming highly anisotropic in the bottom aquifer; 〈q_y〉 was anisotropic in all three aquifers with directions of high continuity normal to those of the 〈q_x〉 field; finally, 〈q_z〉 was highly anisotropic in all three aquifers, with higher continuity along the east–west direction. The mean head field was found to be continuous, despite the high heterogeneity of the underlying hydraulic conductivity field. Directions of high continuity were in alignment with field boundaries and mean flow direction. Conditioning did not influence significantly the expected value of the flux terms, with more pronounced being the effect on the standard deviation of the flux vector components. Conditioning reduced the standard deviations of the horizontal flux components by as much as 50% in the bottom aquifer. Variability in the head cross-sections was affected only marginally, with an average 10% reduction in the respective standard deviation. Finally, the location of the conditioning data did not appear to have a significant effect on the surrounding area, with uniform reduction in standard deviations.     

Numerical method of moments for solute transport in a nonstationary flow field conditioned on hydraulic conductivity and head measurements

The unconditional stochastic studies on groundwater flow and solute transport in a nonstationary conductivity field show that the standard deviations of the hydraulic head and solute flux are very large in comparison with their mean values (Zhang et al. in Water Resour Res 36:2107–2120, 2000; Wu et al. in J Hydrol 275:208–228, 2003; Hu et al. in Adv Water Resour 26:513–531, 2003). In this study, we develop a numerical method of moments conditioning on measurements of hydraulic conductivity and head to reduce the variances of the head and the solute flux. A Lagrangian perturbation method is applied to develop the framework for solute transport in a nonstationary flow field. Since analytically derived moments equations are too complicated to solve analytically, a numerical finite difference method is implemented to obtain the solutions. Instead of using an unconditional conductivity field as an input to calculate groundwater velocity, we combine a geostatistical method and a method of moment for flow to conditionally simulate the distributions of head and velocity based on the measurements of hydraulic conductivity and head at some points. The developed theory is applied in several case studies to investigate the influences of the measurements of hydraulic conductivity and/or the hydraulic head on the variances of the predictive head and the solute flux in nonstationary flow fields. The study results show that the conditional calculation will significantly reduce the head variance. Since the hydraulic head measurement points are treated as the interior boundary (Dirichlet boundary) conditions, conditioning on both the hydraulic conductivity and the head measurements is much better than conditioning only on conductivity measurements for reduction of head variance. However, for solute flux, variance reduction by the conditional study is not so significant.     

Near‐saturated surface soil hydraulic properties under different land uses in the St Denis National Wildlife Area,Saskatchewan, Canada

Surface soil hydraulic properties are key factors controlling the partition of rainfall and snowmelt into runoff and soil water storage, and their knowledge is needed for sound land management. The objective of this study was to evaluate the effects of three land uses (native grass, brome grass and cultivated) on surface soil hydraulic properties under near‐saturated conditions at the St Denis National Wildlife Area, Saskatchewan, Canada. For each land use, water infiltration rates were measured using double‐ring and tension infiltrometers at ?0·3, ?0·7, ?1·5 and ?2·2 kPa pressure heads. Macroporosity and unsaturated hydraulic properties of the surface soil were estimated. Mean field‐saturated hydraulic conductivity (K_fs), unsaturated hydraulic conductivity at ?0·3 kPa pressure head, inverse capillary length scale (α) and water‐conducting macroporosity were compared for different land uses. These parameters of the native grass and brome grass sites were significantly (p < 0·1) higher than that of the cultivated sites. At the ?0·3 kPa pressure head, hydraulic conductivity of grasslands was two to three times greater than that of cultivated lands. Values of α were about two times and values of K_fs about four times greater in grasslands than in cultivated fields. Water‐conducting macroporosity of grasslands and cultivated fields were 0·04% and 0·01% of the total soil volume, respectively. Over 90% of the total water flux at ?0·06 kPa pressure head was transmitted through pores > 1·36 × 10^?4 m in diameter in the three land uses. Land use modified near‐saturated hydraulic properties of surface soil and consequently may alter the water balance of the area by changing the amount of surface runoff and soil water storage. Copyright © 2004 John Wiley & Sons, Ltd.     

Synergistic effects of water repellency and macropore flow on the hydraulic conductivity of a burned forest soil,south‐east Australia

Research shows that water repellency is a key hydraulic property that results in reduced infiltration rates in burned soils. However, more work is required in order to link the hydrological behaviour of water repellent soils to observed runoff responses at the plot and hillslope scale. This study used 5 M ethanol and water in disc infiltrometers to quantify the role of macropore flow and water repellency on spatial and temporal infiltration patterns in a burned soil at plot (<10 m²) scale in a wet eucalypt forest in south‐east Australia. In the first summer and winter after wildfire, an average of 70% and 60%, respectively, of the plot area was water repellent and did not contribute to infiltration. Macropores (r > 0·5 mm), comprising just 5·5% of the soil volume, contributed to 70% and 95%, respectively, of the field‐saturated and ponded hydraulic conductivity (K_p). Because flow occurred almost entirely via macropores in non‐repellent areas, this meant that less than 2·5% of the soil surface effectively contributed to infiltration. The hydraulic conductivity increased by a factor of up to 2·5 as the hydraulic head increased from 0 to 5 mm. Due to the synergistic effect of macropore flow and water repellency, the coefficient of variation (CV) in K_p was three times higher in the water‐repellent soil (CV = 175%) than under the simulated non‐repellent conditions (CV = 66%). The high spatial variability in K_p would act to reduce the effective infiltration rate during runoff generation at plot scale. Ponding, which tend to increase with increasing scale, activates flow through macropores and would raise the effective infiltration rates at larger scales. Field experiments designed to provide representative measurements of infiltration after fire in these systems must therefore consider both the inherent variability in hydraulic conductivity and the variability in infiltration caused by interactions between surface runoff and hydraulic conductivity. Copyright © 2010 John Wiley & Sons, Ltd.     

Determining hydraulic properties of a loam soil by alternative infiltrometer techniques

Testing infiltrometer techniques to determine soil hydraulic properties is necessary for specific soils. For a loam soil, the water retention and hydraulic conductivity predicted by the BEST (Beerkan Estimation of Soil Transfer parameters) procedure of soil hydraulic characterization was compared with data collected by more standard laboratory and field techniques. Six infiltrometer techniques were also compared in terms of saturated soil hydraulic conductivity, K_s. BEST yielded water retention values statistically similar to those obtained in the laboratory and K_s values practically coinciding with those determined in the field with the pressure infiltrometer (PI). The unsaturated soil hydraulic conductivity measured with the tension infiltrometer (TI) was reproduced satisfactorily by BEST only close to saturation. BEST, the PI, one‐potential experiments with both the TI and the mini disk infiltrometer (MDI), the simplified falling head (SFH) technique and the bottomless bucket (BB) method yielded statistically similar estimates of K_s, differing at the most by a factor of three. Smaller values were obtained with longer and more soil‐disturbing infiltration runs. Any of the tested infiltration techniques appears usable to obtain the order of magnitude of K_s at the field site, but the BEST, BB and PI data appear more appropriate to characterize the soil at some stage during a rainfall event. Additional investigations on both similar and different soils would allow development of more general procedures to apply infiltrometer techniques for soil hydraulic characterization. Copyright © 2015 John Wiley & Sons, Ltd.     

Alternative analysis of transient infiltration experiment to estimate soil water repellency

The repellency index (RI) defined as the adjusted ratio between soil‐ethanol, S_e, and soil‐water, S_w, sorptivities estimated from minidisk infiltrometer experiments has been used instead of the widely used water drop penetration time and molarity of ethanol drop tests to assess soil water repellency. However, sorptivity calculated by the usual early‐time infiltration equation may be overestimated as the effects of gravity and lateral capillary are neglected. With the aim to establish the best applicative procedure to assess RI, different approaches to estimate S_e and S_w were compared that make use of both the early‐time infiltration equation (namely, the 1 min, S1, and the short‐time linearization approaches), and the two‐term axisymmetric infiltration equation, valid for early to intermediate times (namely, the cumulative linearization and differentiated linearization approaches). The dataset included 85 minidisk infiltrometer tests conducted in three sites in Italy and Spain under different vegetation habitats (forest of Pinus pinaster and Pinus halepensis, burned pine forest, and annual grasses), soil horizons (organic and mineral), postfire treatments, and initial soil water contents. The S1 approach was inapplicable in 42% of experiments as water infiltration did not start in the first minute. The short‐time linearization approach yielded a systematic overestimation of S_e and S_w that resulted in an overestimation of RI by a factor of 1.57 and 1.23 as compared with the cumulative linearization and differentiated linearization approaches. A new repellency index, RI_s, was proposed as the ratio between the slopes of the linearized data for the wettable and hydrophobic stages obtained by a single water infiltration test. For the experimental conditions considered, RI_s was significantly correlated with RI and WDPT. Compared with RI, RI_s includes information on both soil sorptivity and hydraulic conductivity and, therefore, it can be considered more physically linked to the hydrological processes affected by soil water repellency.     

Impact of land use on the hydraulic properties of the topsoil in a small French catchment

The hydraulic properties of the topsoil control the partition of rainfall into infiltration and runoff at the soil surface. They must be characterized for distributed hydrological modelling. This study presents the results of a field campaign documenting topsoil hydraulic properties in a small French suburban catchment (7 km²) located near Lyon, France. Two types of infiltration tests were performed: single ring infiltration tests under positive head and tension‐disk infiltration using a mini‐disk. Both categories were processed using the BEST—Beerkan Estimation of Soil Transfer parameters—method to derive parameters describing the retention and hydraulic conductivity curves. Dry bulk density and particle size data were also sampled. Almost all the topsoils were found to belong to the sandy loam soil class. No significant differences in hydraulic properties were found in terms of pedologic units, but the results showed a high impact of land use on these properties. The lowest dry bulk density values were obtained in forested soils with the highest organic matter content. Permanent pasture soils showed intermediate values, whereas the highest values were encountered in cultivated lands. For saturated hydraulic conductivity, the highest values were found in broad‐leaved forests and small woods. The complementary use of tension‐disk and positive head infiltration tests highlighted a sharp increase of hydraulic conductivity between near saturation and saturated conditions, attributed to macroporosity effect. The ratio of median saturated hydraulic conductivity to median hydraulic conductivity at a pressure of − 20 mm of water was about 50. The study suggests that soil texture, such as used in most pedo‐transfer functions, might not be sufficient to properly map the variability of soil hydraulic properties. Land use information should be considered in the parameterizations of topsoil within hydrological models to better represent in situ conditions, as illustrated in the paper. Copyright © 2010 John Wiley & Sons, Ltd.     

A Correction on Coastal Heads for Groundwater Flow Models

We introduce a simple correction to coastal heads for constant‐density groundwater flow models that contain a coastal boundary, based on previous analytical solutions for interface flow. The results demonstrate that accurate discharge to the sea in confined aquifers can be obtained by direct application of Darcy's law (for constant‐density flow) if the coastal heads are corrected to ((α + 1)/α)h_s ? B/2α, in which h_s is the mean sea level above the aquifer base, B is the aquifer thickness, and α is the density factor. For unconfined aquifers, the coastal head should be assigned the value . The accuracy of using these corrections is demonstrated by consistency between constant‐density Darcy's solution and variable‐density flow numerical simulations. The errors introduced by adopting two previous approaches (i.e., no correction and using the equivalent fresh water head at the middle position of the aquifer to represent the hydraulic head at the coastal boundary) are evaluated. Sensitivity analysis shows that errors in discharge to the sea could be larger than 100% for typical coastal aquifer parameter ranges. The location of observation wells relative to the toe is a key factor controlling the estimation error, as it determines the relative aquifer length of constant‐density flow relative to variable‐density flow. The coastal head correction method introduced in this study facilitates the rapid and accurate estimation of the fresh water flux from a given hydraulic head measurement and allows for an improved representation of the coastal boundary condition in regional constant‐density groundwater flow models.     

Simple Hydraulic Conductivity Estimation by the Kalman Filtered Double Constraint Method

This paper presents the Kalman Filtered Double Constraint Method (DCM‐KF) as a technique to estimate the hydraulic conductivities in the grid blocks of a groundwater flow model. The DCM is based on two forward runs with the same initial grid block conductivities, but with alternating flux‐head conditions specified on parts of the boundary and the wells. These two runs are defined as: (1) the flux run, with specified fluxes (recharge and well abstractions), and (2) the head run, with specified heads (measured in piezometers). Conductivities are then estimated as the initial conductivities multiplied by the fluxes obtained from the flux run and divided by the fluxes obtained from the head run. The DCM is easy to implement in combination with existing models (e.g., MODFLOW). Sufficiently accurate conductivities are obtained after a few iterations. Because of errors in the specified head‐flux couples, repeated estimation under varying hydrological conditions results in different conductivities. A time‐independent estimate of the conductivities and their inaccuracy can be obtained by a simple linear KF with modest computational requirements. For the Kleine Nete catchment, Belgium, the DCM‐KF yields sufficiently accurate calibrated conductivities. The method also results in distinguishing regions where the head‐flux observations influence the calibration from areas where it is not able to influence the hydraulic conductivity.     

Screening‐Level Sensitivity Analysis for the Design of Pump‐and‐Treat Systems

Pump‐and‐treat systems can prevent the migration of groundwater contaminants and candidate systems are typically evaluated with groundwater models. Such models should be rigorously assessed to determine predictive capabilities and numerous tools and techniques for model assessment are available. While various assessment methodologies (e.g., model calibration, uncertainty analysis, and Bayesian inference) are well‐established for groundwater modeling, this paper calls attention to an alternative assessment technique known as screening‐level sensitivity analysis (SLSA). SLSA can quickly quantify first‐order (i.e., main effects) measures of parameter influence in connection with various model outputs. Subsequent comparisons of parameter influence with respect to calibration vs. prediction outputs can suggest gaps in model structure and/or data. Thus, while SLSA has received little attention in the context of groundwater modeling and remedial system design, it can nonetheless serve as a useful and computationally efficient tool for preliminary model assessment. To illustrate the use of SLSA in the context of designing groundwater remediation systems, four SLSA techniques were applied to a hypothetical, yet realistic, pump‐and‐treat case study to determine the relative influence of six hydraulic conductivity parameters. Considered methods were: Taguchi design‐of‐experiments (TDOE); Monte Carlo statistical independence (MCSI) tests; average composite scaled sensitivities (ACSS); and elementary effects sensitivity analysis (EESA). In terms of performance, the various methods identified the same parameters as being the most influential for a given simulation output. Furthermore, results indicate that the background hydraulic conductivity is important for predicting system performance, but calibration outputs are insensitive to this parameter (K_BK). The observed insensitivity is attributed to a nonphysical specified‐head boundary condition used in the model formulation which effectively “staples” head values located within the conductivity zone. Thus, potential strategies for improving model predictive capabilities include additional data collection targeting the K_BK parameter and/or revision of model structure to reduce the influence of the specified head boundary.     

First approximations of soil moisture retention curves using the filter‐paper method

Relationships between gravimetric soil moisture content (w) and matric potential (ϕ), and between volumetric soil moisture content (θ_v) and pressure head (h) were approximated for the unsaturated zone on Long Island, New York. Soil samples were collected from two sites using a hand auger. The soil moisture content was determined using the filter‐paper (w_f) and gravimetric (w) methods, respectively. The w_f was then used in an empirical equation to estimate ϕ_m. Each set of ϕ_m and w was combined with a straight‐line empirical model to obtain a w(ϕ_m) relationship. Soil ϕ_m was converted to h, and w to the volumetric moisture content θ_v, in order to produce a θ_v(h) curve. Graphical and statistical comparison showed that the resulting θ_v(h) curves fell within one order of magnitude of similar curves generated by a more sophisticated non‐linear model developed previously. The simplicity and low cost of the filter‐paper approach described in this study recommends it for preliminary studies of hydraulic properties in the unsaturated zone. Copyright © 2000 John Wiley & Sons, Ltd.     

Effects of the roots of Cynodon dactylon and Schefflera heptaphylla on water infiltration rate and soil hydraulic conductivity

Water infiltration rate and hydraulic conductivity in vegetated soil are two vital hydrological parameters for agriculturists to determine availability of soil moisture for assessing crop growths and yields, and also for engineers to carry out stability calculations of vegetated slopes. However, any effects of roots on these two parameters are not well‐understood. This study aims to quantify the effects of a grass species, Cynodon dactylon, and a tree species, Schefflera heptaphylla, on infiltration rate and hydraulic conductivity in relation to their root characteristics and suction responses. The two selected species are commonly used for ecological restoration and rehabilitation in many parts of the world and South China, respectively. A series of in‐situ double‐ring infiltration tests was conducted during a wet summer, while the responses of soil suction were monitored by tensiometers. When compared to bare soil, the vegetated soil has lower infiltration rate and hydraulic conductivity. This results in at least 50% higher suction retained in the vegetated soil. It is revealed that the effects of root‐water uptake by the selected species on suction were insignificant because of the small evapotranspiration (<0.2 mm) when the tests were conducted under the wet climate. There appears to have no significant difference (less than 10%) of infiltration rates, hydraulic conductivity and suction retained between the grass‐covered and the tree‐covered soil. However, the grass and tree species having deeper root depth and greater Root Area Index (RAI) retained higher suction. Copyright © 2015 John Wiley & Sons, Ltd.     

Development and parameterization of an infiltration model accounting for water depth and rainfall intensity

Experimental work has clearly shown that the effective hydraulic conductivity (K_e) or effective infiltration rate (f_e) on the local scale of a plot cannot be considered as constant but are dependent on water depth and rainfall intensity because non‐random microtopography‐related variations in hydraulic conductivity occur. Rainfall–runoff models generally do not account for this: models assume that excess water is uniformly spread over the soil surface and within‐plot variations are neglected. In the present study, we propose a model that is based on the concepts of microtopography‐related water depth‐dependent infiltration and partial contributing area. Expressions for the plot scale K_e and f_e were developed that depend on rainfall intensity and runon from upslope (and thus on water depth). To calibrate and validate the model, steady state infiltration experiments were conducted on maize fields on silt loam soils in Belgium, with different stages and combinations of rainfall intensity and inflow, simulating rainfall and runon. Water depth–discharge and depth–inundation relationships were established and used to estimate the effect of inundation on K_e. Although inflow‐only experiments were found to be unsuitable for calibration, the model was successfully calibrated and validated with the rainfall simulation data and combined rainfall–runon data (R²: 0.43–0.91). Calibrated and validated with steady state infiltration experiments, the model was combined with the Green–Ampt infiltration equation and can be applied within a two‐dimensional distributed rainfall–runoff model. The effect of water depth–dependency and rainfall intensity on infiltration was illustrated for a hillslope. Copyright © 2012 John Wiley & Sons, Ltd.     

Estimating field‐saturated soil hydraulic conductivity by a simplified Beerkan infiltration experiment

Field‐saturated soil hydraulic conductivity, K_fs, is highly variable. Therefore, interpreting and simulating hydrological processes, such as rainfall excess generation, need a large number of K_fs data even at the plot scale. Simple and reasonably rapid experiments should be carried out in the field. In this investigation, a simple infiltration experiment with a ring inserted shortly into the soil and the estimation of the so‐called α* parameter allowed to obtain an approximate measurement of K_fs. The theoretical approach was tested with reference to 149 sampling points established on Burundian soils. The estimated K_fs with the value of first approximation of α* for most agricultural field soils (α* = 0.012 mm^?1) differed by a practically negligible maximum factor of two from the saturated conductivity obtained by the complete Beerkan Estimation of Soil Transfer parameters (BEST) procedure for soil hydraulic characterization. The measured infiltration curve contained the necessary information to obtain a site‐specific prediction of α*. The empirically derived α* relationship gave similar results for K_fs (mean = 0.085 mm s^?1; coefficient of variation (CV) = 71%) to those obtained with BEST (mean = 0.086 mm s^?1; CV = 67%), and it was also successfully tested with reference to a few Sicilian sampling points, since it yielded a mean and a CV of K_fs (0.0094 mm s^?1 and 102%, respectively) close to the values obtained with BEST (mean = 0.0092 mm s^?1; CV = 113%). The developed method appears attractive due to the extreme simplicity of the experiment. Copyright © 2012 John Wiley & Sons, Ltd.     

Analytical Interpretation of Slug Test in a Vertical Cutoff Wall

An analysis method for slug tests performed in a partially penetrating well within a vertical cutoff wall is presented. A steady‐state shape factor for evaluating hydraulic conductivity of the material within the wall was derived by applying the method of images to the previously developed analytical solution of Zlotnik et al. (2010) for an infinite aquifer. Two distinct boundary conditions were considered: constant‐head boundary for the case of direct contact between the wall and the aquifer, and no‐flux boundary representing an impermeable filter cake on the sides of the wall. The constant‐head and no‐flux boundary conditions yield significantly higher and lower shape factors, respectively, than those for the infinite aquifer. Consequently the conventional line‐fitting method for slug test analysis would yield an inaccurate estimate of the hydraulic conductivity of a vertical cutoff wall.     

Sensitivity analysis of hydraulic conductivity and Manning's n parameters lead to new method to scale effective hydraulic conductivity across model resolutions

Hydrological modelling is an important tool for research, policy, and management, but uncertainty remains about parameters transferability from field observations made at small scale to models at the catchment scale and larger. This uncertainty compels the need to develop parameter relationships that are translatable across scale. In this study, we compare the changes to modelled processes as resolution is coarsened from 100‐m to 1‐km in a topographically complex, 255‐km² Colorado River headwater catchment. We conducted a sensitivity analysis for hydraulic conductivity (K) and Manning's n parameters across four orders of magnitude. Results showed that K acts as a moderator between surface and subsurface contributions to streamflow, whereas n moderates the duration of high intensity, infiltration‐excess flow. The parametric sensitivity analysis informed development of a new method to scale effective hydraulic conductivity across modelling resolutions in order to compensate for the loss of topographic gradients as resolution is coarsened. A similar mathematical relationship between n and lateral resolution changes was not found, possibly because n is also sensitive to time discretization. This research provides an approach to translate hydraulic conductivity parameters from a calibrated coarse model to higher resolutions where the number of simulations are limited by computational demand.     

Spatial scale effect on the upper soil effective parameters of a distributed hydrological model

Nonlinear dynamics and spatial variability in hydrological systems make the formulation of scaling theories difficult. Therefore, the development of knowledge related to scale effects, scaling techniques, parameterization and linkages of parameters across scales is highly relevant. The main purpose of this work is to analyse the spatial effect of the static storage capacity parameter H_u and the saturated hydraulic conductivity parameter k_s from microscale (sub‐grid level) to mesoscale (grid level) and its implication to the definition of an optimum cell size. These two parameters describe the upper soil water characteristics in the infiltration process conceptualization of the TETIS hydrological model. At microscale, the spatial heterogeneity of H_u and k_s was obtained generating random parameter fields through probability distribution functions and a spatial dependence model with pre‐established correlation lengths. The effective parameters at mesoscale were calculated by solving the inverse problem for each parameter field. Results indicate that the adopted inverse formulation allows transferring the nonlinearity of the system from microscale to the mesoscale via non‐stationary effective parameters. Their values at each cell and time step are in the range of zero to the mean value of the parameter at microscale. The stochastic simulations showed that the variance of the estimated effective parameters decreases when the ratio between mesoscale cell size and correlation length at microscale increases. For a ratio greater than 1, we found cell sizes having the characteristics of a representative elementary area (REA); in such case, the microscale variability pattern did not affect the system response at mesoscale. Copyright © 2011 John Wiley & Sons, Ltd.     

An Integrated Laboratory Method to Measure and Verify Directional Hydraulic Conductivity in Fine‐to‐Medium Sandy Sediments

The constant‐head permeameter test (CHPT) is widely used in sandy samples as a standard method in the laboratory to investigate hydraulic conductivity (K). However, it neither can be used to consistently determine directional hydraulic conductivity (DHC) nor guarantee the comparability of measured K values of samples with different sizes. Therefore, this paper proposes an integrated laboratory method, called modified CHPT (MCHPT), for the efficient determination and verification of consistent DHC values in fine‐to‐medium sandy sediments, based on a new methodological framework. A precise and standardized procedure for preparing the experimental setup of MCHPT was conducted, based on the integrated experimental setup of CHPT and tracer tests. Moreover, a formula was yielded for the time‐optimized sample saturation control. In comparison with grain size‐based methods, the validity of consistent K_h and K_v values determined by MCHPT was convincing.     

A Tube Seepage Meter for In Situ Measurement of Seepage Rate and Groundwater Sampling

We designed and evaluated a “tube seepage meter” for point measurements of vertical seepage rates (q), collecting groundwater samples, and estimating vertical hydraulic conductivity (K) in streambeds. Laboratory testing in artificial streambeds show that seepage rates from the tube seepage meter agreed well with expected values. Results of field testing of the tube seepage meter in a sandy‐bottom stream with a mean seepage rate of about 0.5 m/day agreed well with Darcian estimates (vertical hydraulic conductivity times head gradient) when averaged over multiple measurements. The uncertainties in q and K were evaluated with a Monte Carlo method and are typically 20% and 60%, respectively, for field data, and depend on the magnitude of the hydraulic gradient and the uncertainty in head measurements. The primary advantages of the tube seepage meter are its small footprint, concurrent and colocated assessments of q and K, and that it can also be configured as a self‐purging groundwater‐sampling device.