Impact of sampling volume on the probability density function of steady state concentration

Abstract

[1] In recent years, statistical theory has been used to compute the ensemble mean and variance of solute concentration in aquifer formations with second-order stationary velocity fields. The merit of accurately estimating the mean and variance of concentration, however, remains unclear without knowing the shape of the probability density function (pdf). In a setup where a conservative solute is continuously injected into a domain, the concentration is bounded between zero and the concentration value in the injected solution. At small travel distances close to the fringe of the plume, an observation point may fall into the plume or outside, so that the statistical concentration distribution clusters at the two limiting values. Obviously, this results in non-Gaussian pdf's of concentration. With increasing travel distance, the lateral plume boundaries are smoothed, resulting in increased probability of intermediate concentrations. Likewise, averaging the concentration in a larger sampling volume, as typically done in field measurements, leads to higher probabilities of intermediate concentrations. We present semianalytical results of concentration pdf's for measurements with point-like or larger support volumes based on stochastic theory applied to stationary media. To this end, we employ a reversed auxiliary transport problem, in which we use analytical expressions for first and second central spatial lateral moments with an assumed Gaussian pdf for the uncertainty of the first lateral moment and Gauss-like shapes in individual cross sections. The resulting concentration pdf can be reasonably fitted by beta distributions. The results are compared to Monte Carlo simulations of flow and steady state transport in 3-D heterogeneous domains. In both methods the shape of the concentration pdf changes with distance to the contaminant source: Near the source, the distribution is multimodal, whereas it becomes a unimodal beta distribution far away from the contaminant source. The semianalytical and empirical pdf's differ slightly, which we contribute to the numerical artifacts in the Monte Carlo simulations but also to hard assumptions made in the semianalytical approach. Our results imply that geostatistical techniques for interpolation and other statistical inferences based on Gaussian distributions, such as kriging and cokriging, may be feasible only far away from the contaminant source. For calculations near the source, the beta-like distribution of concentration should be accounted for.

Concentration statistics of solute transport for the near field zone of an estuary

“…In order to obtain the probability of exceeding these limit values one must have the estimate of the concentration probability density function. For demonstration purposes we use the first two concentration moments and select the beta distribution for concentration pdf here, as previously used in other environmental flows ( Chatwin et al, 1995;Schwede et al, 2008 ).…”

Kriging groundwater solute concentrations using flow coordinates and nonstationary covariance functions

“…Many studies show that the solute concentrations pdf can be modeled as a scaled beta distribution (Fiorotto and Caroni, 2002;Caroni and Fiorotto, 2005;Bellin and Tonina, 2007;Schwede et al, 2008;Cirpka et al, 2011a). Schwede and Cirpka (2010b) assumed the kriging estimate and estimation variance are the first two moments of a beta distribution.…”

Simulating an exclusion zone for vapour intrusion of TCE from groundwater into indoor air

“…The second strategy begins with the statistical hypothesis that the indoor air TCE concentrations for all Monte Carlo realizations C c TCE IA i;j;l with l ={1,2,…,n realizations } for a given monitoring interval (t i −1 ,t i ] and the jth house are log-normally distributed. In contrast, Schwede et al (2008) suggest average concentrations within a large sampling volume follow a beta distribution. However, their analysis is specific to a groundwater plume whereas this work focuses on that of house located above the water table.…”