Anisotropic statistical scaling of vadose zone hydraulic property estimates near Maricopa, Arizona

A. Guadagnini, Shlomo P Neuman, Marcel Schaap, M. Riva

Research output: Contribution to journalArticle

17 Citations (Scopus)

Abstract

Fluid flow and mass transport in the vadose zone are strongly influenced by spatial variability in soil hydraulic properties. It has become common to characterize spatial variability geostatistically and to solve corresponding flow and/or transport problems stochastically. The typical (though not only) approach is to treat log saturated hydraulic conductivity, Y=log10 Ks, and perhaps other medium properties as statistically homogeneous, isotropic, or anisotropic multivariate-Gaussian random fields with unique variance and autocorrelation scale. A growing body of literature suggests that Y as well as many other variables may possess heavy-tailed, non-Gaussian distributions, and/or scale-dependent statistical parameters representing a multiscale, hierarchical structure. Elsewhere, we have demonstrated that these important phenomena are difficult to detect, and are not fully tractable, with standard geostatistical methods. They are however detectable and tractable with a novel geostatistical method of analysis which treats such variables as samples from sub-Gaussian random fields subordinated to truncated fractional Brownian motion (tfBm) or truncated fractional Gaussian noise (tfGn). Each such field is a mixture of Gaussian components having random variances. The purpose of this paper is to explore the extent to which hydraulic parameters of unsaturated soils exhibit heavy-tailed distributions and statistical scaling of the above type. As a test bed we have selected an experimental site near Maricopa, Arizona, at which soil texture data have been collected in boreholes and trenches to a depth of 15 m over an area of 3600 m2. Since hydraulic parameters are difficult to measure to such depths we estimate them using a neural network pedotransfer model. Input data include percent sand, silt, and clay. Outputs include estimates of saturated and residual volumetric water content, saturated hydraulic conductivity, and (due to their wide applicability) parameters of the van Genuchten-Mualem constitutive model of unsaturated soil property variations with capillary pressure head. We find indeed that vertical and horizontal increments of our hydraulic parameter estimates exhibit heavy-tailed frequency distributions and statistical scaling which conform closely to our novel geostatistical framework. Its generality and applicability to the Maricopa site has far reaching implications vis-à-vis the geostatistical characterization of vadose zone properties and the stochastic modeling of flow and transport in unsaturated media at other sites.

Original languageEnglish (US)
Pages (from-to)8463-8479
Number of pages17
JournalWater Resources Research
Volume49
Issue number12
DOIs
StatePublished - Dec 2013

Fingerprint

hydraulic property
vadose zone
hydraulics
hydraulic conductivity
unsaturated medium
Brownian motion
capillary pressure
mass transport
soil texture
autocorrelation
fluid flow
trench
soil property
silt
soil
borehole
water content
parameter
clay
sand

ASJC Scopus subject areas

  • Water Science and Technology

Cite this

Anisotropic statistical scaling of vadose zone hydraulic property estimates near Maricopa, Arizona. / Guadagnini, A.; Neuman, Shlomo P; Schaap, Marcel; Riva, M.

In: Water Resources Research, Vol. 49, No. 12, 12.2013, p. 8463-8479.

Research output: Contribution to journalArticle

@article{1f1502d0878d48fb817dded8faa7a561,
title = "Anisotropic statistical scaling of vadose zone hydraulic property estimates near Maricopa, Arizona",
abstract = "Fluid flow and mass transport in the vadose zone are strongly influenced by spatial variability in soil hydraulic properties. It has become common to characterize spatial variability geostatistically and to solve corresponding flow and/or transport problems stochastically. The typical (though not only) approach is to treat log saturated hydraulic conductivity, Y=log10 Ks, and perhaps other medium properties as statistically homogeneous, isotropic, or anisotropic multivariate-Gaussian random fields with unique variance and autocorrelation scale. A growing body of literature suggests that Y as well as many other variables may possess heavy-tailed, non-Gaussian distributions, and/or scale-dependent statistical parameters representing a multiscale, hierarchical structure. Elsewhere, we have demonstrated that these important phenomena are difficult to detect, and are not fully tractable, with standard geostatistical methods. They are however detectable and tractable with a novel geostatistical method of analysis which treats such variables as samples from sub-Gaussian random fields subordinated to truncated fractional Brownian motion (tfBm) or truncated fractional Gaussian noise (tfGn). Each such field is a mixture of Gaussian components having random variances. The purpose of this paper is to explore the extent to which hydraulic parameters of unsaturated soils exhibit heavy-tailed distributions and statistical scaling of the above type. As a test bed we have selected an experimental site near Maricopa, Arizona, at which soil texture data have been collected in boreholes and trenches to a depth of 15 m over an area of 3600 m2. Since hydraulic parameters are difficult to measure to such depths we estimate them using a neural network pedotransfer model. Input data include percent sand, silt, and clay. Outputs include estimates of saturated and residual volumetric water content, saturated hydraulic conductivity, and (due to their wide applicability) parameters of the van Genuchten-Mualem constitutive model of unsaturated soil property variations with capillary pressure head. We find indeed that vertical and horizontal increments of our hydraulic parameter estimates exhibit heavy-tailed frequency distributions and statistical scaling which conform closely to our novel geostatistical framework. Its generality and applicability to the Maricopa site has far reaching implications vis-{\`a}-vis the geostatistical characterization of vadose zone properties and the stochastic modeling of flow and transport in unsaturated media at other sites.",
author = "A. Guadagnini and Neuman, {Shlomo P} and Marcel Schaap and M. Riva",
year = "2013",
month = "12",
doi = "10.1002/2013WR014286",
language = "English (US)",
volume = "49",
pages = "8463--8479",
journal = "Water Resources Research",
issn = "0043-1397",
publisher = "American Geophysical Union",
number = "12",

}

TY - JOUR

T1 - Anisotropic statistical scaling of vadose zone hydraulic property estimates near Maricopa, Arizona

AU - Guadagnini, A.

AU - Neuman, Shlomo P

AU - Schaap, Marcel

AU - Riva, M.

PY - 2013/12

Y1 - 2013/12

N2 - Fluid flow and mass transport in the vadose zone are strongly influenced by spatial variability in soil hydraulic properties. It has become common to characterize spatial variability geostatistically and to solve corresponding flow and/or transport problems stochastically. The typical (though not only) approach is to treat log saturated hydraulic conductivity, Y=log10 Ks, and perhaps other medium properties as statistically homogeneous, isotropic, or anisotropic multivariate-Gaussian random fields with unique variance and autocorrelation scale. A growing body of literature suggests that Y as well as many other variables may possess heavy-tailed, non-Gaussian distributions, and/or scale-dependent statistical parameters representing a multiscale, hierarchical structure. Elsewhere, we have demonstrated that these important phenomena are difficult to detect, and are not fully tractable, with standard geostatistical methods. They are however detectable and tractable with a novel geostatistical method of analysis which treats such variables as samples from sub-Gaussian random fields subordinated to truncated fractional Brownian motion (tfBm) or truncated fractional Gaussian noise (tfGn). Each such field is a mixture of Gaussian components having random variances. The purpose of this paper is to explore the extent to which hydraulic parameters of unsaturated soils exhibit heavy-tailed distributions and statistical scaling of the above type. As a test bed we have selected an experimental site near Maricopa, Arizona, at which soil texture data have been collected in boreholes and trenches to a depth of 15 m over an area of 3600 m2. Since hydraulic parameters are difficult to measure to such depths we estimate them using a neural network pedotransfer model. Input data include percent sand, silt, and clay. Outputs include estimates of saturated and residual volumetric water content, saturated hydraulic conductivity, and (due to their wide applicability) parameters of the van Genuchten-Mualem constitutive model of unsaturated soil property variations with capillary pressure head. We find indeed that vertical and horizontal increments of our hydraulic parameter estimates exhibit heavy-tailed frequency distributions and statistical scaling which conform closely to our novel geostatistical framework. Its generality and applicability to the Maricopa site has far reaching implications vis-à-vis the geostatistical characterization of vadose zone properties and the stochastic modeling of flow and transport in unsaturated media at other sites.

AB - Fluid flow and mass transport in the vadose zone are strongly influenced by spatial variability in soil hydraulic properties. It has become common to characterize spatial variability geostatistically and to solve corresponding flow and/or transport problems stochastically. The typical (though not only) approach is to treat log saturated hydraulic conductivity, Y=log10 Ks, and perhaps other medium properties as statistically homogeneous, isotropic, or anisotropic multivariate-Gaussian random fields with unique variance and autocorrelation scale. A growing body of literature suggests that Y as well as many other variables may possess heavy-tailed, non-Gaussian distributions, and/or scale-dependent statistical parameters representing a multiscale, hierarchical structure. Elsewhere, we have demonstrated that these important phenomena are difficult to detect, and are not fully tractable, with standard geostatistical methods. They are however detectable and tractable with a novel geostatistical method of analysis which treats such variables as samples from sub-Gaussian random fields subordinated to truncated fractional Brownian motion (tfBm) or truncated fractional Gaussian noise (tfGn). Each such field is a mixture of Gaussian components having random variances. The purpose of this paper is to explore the extent to which hydraulic parameters of unsaturated soils exhibit heavy-tailed distributions and statistical scaling of the above type. As a test bed we have selected an experimental site near Maricopa, Arizona, at which soil texture data have been collected in boreholes and trenches to a depth of 15 m over an area of 3600 m2. Since hydraulic parameters are difficult to measure to such depths we estimate them using a neural network pedotransfer model. Input data include percent sand, silt, and clay. Outputs include estimates of saturated and residual volumetric water content, saturated hydraulic conductivity, and (due to their wide applicability) parameters of the van Genuchten-Mualem constitutive model of unsaturated soil property variations with capillary pressure head. We find indeed that vertical and horizontal increments of our hydraulic parameter estimates exhibit heavy-tailed frequency distributions and statistical scaling which conform closely to our novel geostatistical framework. Its generality and applicability to the Maricopa site has far reaching implications vis-à-vis the geostatistical characterization of vadose zone properties and the stochastic modeling of flow and transport in unsaturated media at other sites.

UR - http://www.scopus.com/inward/record.url?scp=84896711774&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84896711774&partnerID=8YFLogxK

U2 - 10.1002/2013WR014286

DO - 10.1002/2013WR014286

M3 - Article

AN - SCOPUS:84896711774

VL - 49

SP - 8463

EP - 8479

JO - Water Resources Research

JF - Water Resources Research

SN - 0043-1397

IS - 12

ER -