David L. Freyberg

Bio: David L. Freyberg is an academic researcher from Stanford University. The author has contributed to research in topics: Aquifer & Groundwater. The author has an hindex of 23, co-authored 36 publications receiving 2241 citations.

A natural gradient experiment on solute transport in a sand aquifer: 2. Spatial moments and the advection and dispersion of nonreactive tracers

Abstract: The three-dimensional movement of a tracer plume containing bromide and chloride is investigated using the data base from a large-scale natural gradient field experiment on groundwater solute transport. The analysis focuses on the zeroth-, first-, and second-order spatial moments of the concentration distribution. These moments define integrated measures of the dissolved mass, mean solute velocity, and dispersion of the plume. Moments are estimated from the point observations using quadrature approximations tailored to the density of the sampling network. The estimators appear to be robust, with acceptable sampling variability. Estimates of the mass in solution for both bromide and chloride demonstrate that the tracers behaved conservatively, as expected. Analysis of the first-order moment estimates indicates that the experimental tracer plumes traveled along identical trajectories. The horizontal trajectory is linear and aligned with the hydraulic gradient. The vertical trajectory is curvilinear, concave upward. The total vertical displacement is small, however, so that the vertical component of the mean solute velocity vector is negligible. The estimated mean solute velocity is identical for both tracers (0.091 m/day) and is spatially and temporally uniform for the first 647 days of travel time. After 647 days of transport, the plume apparently encountered a relatively large-scale heterogeneity in the velocity field, leading to a distinct vertical layering, and slowing the rate of advance of the center of mass of the plume as a whole. The estimated horizontal components of the covariance tensor evolve over time in a manner consistent with the qualitative shape changes observed from plots of the concentration data. The major principal axis, initially aligned roughly perpendicular to the hydraulic gradient, rotates smoothly over time until it is nearly aligned with the mean solute velocity vector, as the plume itself elongates and orients its long axis with the direction of movement. Plots of the components of the covariance tensor as functions of time show evidence of what is commonly called “scale-dependent” dispersion: the rate of growth of the covariance over time is not linear. The theoretical results of G. Dagan (1984) calibrate well to the estimated covariance data for the first 647 days of transport. The calibrated values of the parameters of the hydrualic conductivity distribution closely match independently measured values from the site. The asymptotic longitudinal dispersivity obtained from the calibration is 0.49 m, although asymptotic conditions were apparently not reached. The estimated covariance terms for the last sampling session, 1038 days after injection, are inconsistent with the earlier data and with the Dagan model, particularly for the transverse and off-diagonal components. This behavior is probably attributable to the observed large-scale heterogeneity in the velocity field.

A natural gradient experiment on solute transport in a sand aquifer: 1. Approach and overview of plume movement

Abstract: A large-scale field experiment on natural gradient transport of solutes in groundwater has been conducted at a site in Borden, Ontario. Well-defined initial conditions were achieved by the pulse injection of 12 m3 of a uniform solution containing known masses of two inorganic tracers (chloride and bromide) and five halogenated organic chemicals (bromoform, carbon tetrachloride, tetrachloroethylene, 1,2-dichlorobenzene, and hexachloroethane). A dense, three-dimensional array of over 5000 sampling points was installed throughout the zone traversed by the solutes. Over 19,900 samples have been collected over a 3-year period. The tracers followed a linear horizontal trajectory at an approximately constant velocity, both of which compare well with expectations based on water table contours and estimates of hydraulic head gradient, porosity, and hydraulic conductivity. The vertical displacement over the duration of the experiment was small. Spreading was much more pronounced in the horizontal longitudinal than in the horizontal transverse direction; vertical spreading was very small. The organic solutes were retarded in mobility, as expected.

Use of sedimentological information for geometric simulation of natural porous media structure

Abstract: A geometric simulation method was used to develop a three-dimensional, highly detailed synthetic representation of point bar sediments in the Wabash River system. Geometric simulation methods, in comparison to well-known second-order stochastic methods, offer the advantage of being more closely related to depositional processes, which are often similarly conceptualized (i.e., described in terms of shapes of discrete bed forms, trends in grain size, and spatial relationships of defined geologic facies). Multiple scales of geometric variation were defined within a sedimentologically prescribed framework, and shapes of discrete geometric elements were established at each scale. The selected shapes were based on published field studies including sedimentological bed form studies and trench studies in active point bar sediments. The parameterization of the shapes allowed for random variability of the shape descriptors; discrete shapes were then generated and assimilated by computer. Hydraulic conductivity values were assigned to the discrete elements based on reports of observed variations in grain size and field measurements of hydraulic conductivity. The synthetic model, referred to as a numerical aquifer, is being used as the basis for extensive numerical experimentation to study the relationship between natural spatial structure and subsurface flow and transport.

A Critical Review of Data on Field-Scale Dispersion in Aquifers

Abstract: found that field-scale dispersivities are several orders of.magnitude greater than lab-scale values for the same material; it is generally agreed that this difference is a reflection of the influence of natural heterogeneities which produce irregular flow patterns at the field scale. Consequently, laboratory measurements of dispersivity cannot be used to predict field values of dispersivity. Instead field-scale tracer tests are sometimes conducted to estimate dispersivity at a particular site. Early efforts to document the scale dependence of dispersivity [Lallemand-Barres and Peaudecerf, 1978; Anderson, 1979; Pickens and Grisak, 1981; Beims, 1983; Neretnieks, 1985] were based on field values of dispersivity reported in the literature and the test scales associated with those values. These studies were useful in that they indeed documented field evidence of the scale effect, but they were lacking in that they did not assess the reliability of the data presented. Because we felt that the data would be more

A natural gradient experiment on solute transport in a sand aquifer: Spatial variability of hydraulic conductivity and its role in the dispersion process

Abstract: The spatial variability of hydraulic conductivity at the site of a long-term tracer test performed in the Borden aquifer was examined in great detail by conducting permeability measurements on a series of cores taken along two cross sections, one along and the other transverse to the mean flow direction. Along the two cross sections, a regular-spaced grid of hydraulic conductivity data with 0.05 m vertical and 1.0 m horizontal spatial discretization revealed that the aquifer is comprised of numerous thin, discontinuous lenses of contrasting hydraulic conductivity. Estimation of the three-dimensional covariance structure of the aquifer from the log-transformed data indicates that an exponential covariance model with a variance equal to 0.29, an isotropic horizontal correlation length equal to about 2.8 m, and a vertical correlation length equal to 0.12 m is representative. A value for the longitudinal macrodispersivity calculated from these statistical parameters using three-dimensional stochastic transport theory developed by L. W. Gelhar and C. L. Axness (1983) is about 0.6 m. For the vertically averaged case, the two-dimensional theory developed by G. Dagan (1982, 1984) yields a longitudinal djspersivity equal to 0.45 m. Use of the estimated statistical parameters describing the ln (K) variability in Dagan's transient equations closely predicted the observed longitudinal and horizontal transverse spread of the tracer with time. Weak vertical and horizontal dispersion that is controlled essentially by local-scale dispersion was obtained from the analysis. Because the dispersion predicted independently from the statistical description of the Borden aquifer is consistent with the spread of the injected tracer, it is felt that the theory holds promise for providing meaningful estimates of effective transport parameters in other complex-structured aquifers.

A review of global terrestrial evapotranspiration: observation, modeling, climatology, and climatic variability

Abstract: [1] This review surveys the basic theories, observational methods, satellite algorithms, and land surface models for terrestrial evapotranspiration, E (or λE, i.e., latent heat flux), including a long-term variability and trends perspective. The basic theories used to estimate E are the Monin-Obukhov similarity theory (MOST), the Bowen ratio method, and the Penman-Monteith equation. The latter two theoretical expressions combine MOST with surface energy balance. Estimates of E can differ substantially between these three approaches because of their use of different input data. Surface and satellite-based measurement systems can provide accurate estimates of diurnal, daily, and annual variability of E. But their estimation of longer time variability is largely not established. A reasonable estimate of E as a global mean can be obtained from a surface water budget method, but its regional distribution is still rather uncertain. Current land surface models provide widely different ratios of the transpiration by vegetation to total E. This source of uncertainty therefore limits the capability of models to provide the sensitivities of E to precipitation deficits and land cover change.

Stochastic subsurface hydrology from theory to applications

Abstract: Research on stochastic analysis of subsurface flow has developed rapidly in the last decade, but applications of this approach have been very limited. The purpose of this paper is to illustrate how currently available techniques and results can be used to answer important questions about the large-scale behavior of naturally heterogeneous aquifers. Perturbation-based spectral theory, which presumes local statistical homogeneity, provides generic theoretical results for the head variance, effective conductivity tensor, and macrodispersivity tensor in a field. These results emphasize the key role of the variance and spatial correlation scales of the log hydraulic conductivity field. Field information of variances and correlation scales of natural materials is summarized. The validity of some of the generic stochastic results is evaluated through comparisons with Monte Carlo simulations and field observations. A specific field application example is developed to illustrate how the stochastic results are used to estimate large-scale parameters and determine the reliability of three-dimensional numerical simulations. Using typical log conductivity covariance parameters, the effective hydraulic conductivity tensor, and the macrodispersivity tensor are estimated. The calculated head variance, based on the simulated mean hydraulic gradient, is used as a measure of the adequacy of the calculation of the steady state flow model. Discussion emphasizes limitations and extensions of this approach, and ongoing field evaluations of the results.