### Abstract

We analyze the statistical scaling properties of vertical and horizontal increments in soil and sediment texture data measured to a depth of 15m over an area of 3600m^{2} in a vadose zone near Maricopa, Arizona. The data include sand, silt and clay fractions, their principal components and logit transforms. The statistical scaling properties we uncover are difficult to detect with standard geostatistical methods. They include (a) pronounced intermittency (rough, irregular spatial variability) and antipersistence (tendency for large and small values to alternate rapidly); (b) symmetric, non-Gaussian frequency distributions characterized by heavy tails that often decay with separation distance or lag; (c) nonlinear power-law scaling of sample structure functions (statistical moments of absolute increments) in a midrange of lags, with breakdown in such scaling at small and large lags; (d) extended power-law scaling (linear relations between log structure functions of successive orders) at all lags; (e) nonlinear scaling of power-law exponent with order of sample structure function; and (f) pronounced anisotropy in these behaviors. Similar statistical scaling behaviors are known to be exhibited by a wide variety of earth, environmental and other variables (including ecological, biological, physical, astrophysical and financial). The literature has traditionally interpreted this to imply that the variables are multifractal, which however explains neither the observed breakdown in power-law scaling at small and large lags nor extended power-law scaling. We offer an alternative interpretation that is simpler and consistent with all the above phenomena. Our interpretation views the data as samples from stationary, anisotropic sub-Gaussian random fields subordinated to truncated fractional Brownian motion (tfBm) or truncated fractional Gaussian noise (tfGn). Such sub-Gaussian fields are mixtures of Gaussian fields with random variances. Truncation of monofractal fBm (which is non-stationary) and fGn entails filtering out components below the measurement or resolution scale of the data and above the scale of their sampling domain. Our novel interpretation of the data allows us to obtain maximum likelihood estimates of all parameters characterizing the underlying truncated sub-Gaussian fields. These parameters in turn make it possible to downscale or upscale all statistical moments to situations entailing smaller or larger measurement or resolution and sampling scales, respectively. They also allow one to perform conditional or unconditional Monte Carlo simulations of random field realizations corresponding to these scales.

Original language | English (US) |
---|---|

Pages (from-to) | 217-227 |

Number of pages | 11 |

Journal | Geoderma |

Volume | 214-215 |

DOIs | |

State | Published - Feb 2014 |

### Fingerprint

### Keywords

- Anisotropy
- Antipersistent spatial variability
- Random fields
- Scaling
- Soil and sediment texture
- Truncated fractional Brownian motion

### ASJC Scopus subject areas

- Soil Science

### Cite this

*Geoderma*,

*214-215*, 217-227. https://doi.org/10.1016/j.geoderma.2013.09.008

**Anisotropic statistical scaling of soil and sediment texture in a stratified deep vadose zone near Maricopa, Arizona.** / Guadagnini, A.; Neuman, Shlomo P; Schaap, Marcel; Riva, M.

Research output: Contribution to journal › Article

*Geoderma*, vol. 214-215, pp. 217-227. https://doi.org/10.1016/j.geoderma.2013.09.008

}

TY - JOUR

T1 - Anisotropic statistical scaling of soil and sediment texture in a stratified deep vadose zone near Maricopa, Arizona

AU - Guadagnini, A.

AU - Neuman, Shlomo P

AU - Schaap, Marcel

AU - Riva, M.

PY - 2014/2

Y1 - 2014/2

N2 - We analyze the statistical scaling properties of vertical and horizontal increments in soil and sediment texture data measured to a depth of 15m over an area of 3600m2 in a vadose zone near Maricopa, Arizona. The data include sand, silt and clay fractions, their principal components and logit transforms. The statistical scaling properties we uncover are difficult to detect with standard geostatistical methods. They include (a) pronounced intermittency (rough, irregular spatial variability) and antipersistence (tendency for large and small values to alternate rapidly); (b) symmetric, non-Gaussian frequency distributions characterized by heavy tails that often decay with separation distance or lag; (c) nonlinear power-law scaling of sample structure functions (statistical moments of absolute increments) in a midrange of lags, with breakdown in such scaling at small and large lags; (d) extended power-law scaling (linear relations between log structure functions of successive orders) at all lags; (e) nonlinear scaling of power-law exponent with order of sample structure function; and (f) pronounced anisotropy in these behaviors. Similar statistical scaling behaviors are known to be exhibited by a wide variety of earth, environmental and other variables (including ecological, biological, physical, astrophysical and financial). The literature has traditionally interpreted this to imply that the variables are multifractal, which however explains neither the observed breakdown in power-law scaling at small and large lags nor extended power-law scaling. We offer an alternative interpretation that is simpler and consistent with all the above phenomena. Our interpretation views the data as samples from stationary, anisotropic sub-Gaussian random fields subordinated to truncated fractional Brownian motion (tfBm) or truncated fractional Gaussian noise (tfGn). Such sub-Gaussian fields are mixtures of Gaussian fields with random variances. Truncation of monofractal fBm (which is non-stationary) and fGn entails filtering out components below the measurement or resolution scale of the data and above the scale of their sampling domain. Our novel interpretation of the data allows us to obtain maximum likelihood estimates of all parameters characterizing the underlying truncated sub-Gaussian fields. These parameters in turn make it possible to downscale or upscale all statistical moments to situations entailing smaller or larger measurement or resolution and sampling scales, respectively. They also allow one to perform conditional or unconditional Monte Carlo simulations of random field realizations corresponding to these scales.

AB - We analyze the statistical scaling properties of vertical and horizontal increments in soil and sediment texture data measured to a depth of 15m over an area of 3600m2 in a vadose zone near Maricopa, Arizona. The data include sand, silt and clay fractions, their principal components and logit transforms. The statistical scaling properties we uncover are difficult to detect with standard geostatistical methods. They include (a) pronounced intermittency (rough, irregular spatial variability) and antipersistence (tendency for large and small values to alternate rapidly); (b) symmetric, non-Gaussian frequency distributions characterized by heavy tails that often decay with separation distance or lag; (c) nonlinear power-law scaling of sample structure functions (statistical moments of absolute increments) in a midrange of lags, with breakdown in such scaling at small and large lags; (d) extended power-law scaling (linear relations between log structure functions of successive orders) at all lags; (e) nonlinear scaling of power-law exponent with order of sample structure function; and (f) pronounced anisotropy in these behaviors. Similar statistical scaling behaviors are known to be exhibited by a wide variety of earth, environmental and other variables (including ecological, biological, physical, astrophysical and financial). The literature has traditionally interpreted this to imply that the variables are multifractal, which however explains neither the observed breakdown in power-law scaling at small and large lags nor extended power-law scaling. We offer an alternative interpretation that is simpler and consistent with all the above phenomena. Our interpretation views the data as samples from stationary, anisotropic sub-Gaussian random fields subordinated to truncated fractional Brownian motion (tfBm) or truncated fractional Gaussian noise (tfGn). Such sub-Gaussian fields are mixtures of Gaussian fields with random variances. Truncation of monofractal fBm (which is non-stationary) and fGn entails filtering out components below the measurement or resolution scale of the data and above the scale of their sampling domain. Our novel interpretation of the data allows us to obtain maximum likelihood estimates of all parameters characterizing the underlying truncated sub-Gaussian fields. These parameters in turn make it possible to downscale or upscale all statistical moments to situations entailing smaller or larger measurement or resolution and sampling scales, respectively. They also allow one to perform conditional or unconditional Monte Carlo simulations of random field realizations corresponding to these scales.

KW - Anisotropy

KW - Antipersistent spatial variability

KW - Random fields

KW - Scaling

KW - Soil and sediment texture

KW - Truncated fractional Brownian motion

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

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

U2 - 10.1016/j.geoderma.2013.09.008

DO - 10.1016/j.geoderma.2013.09.008

M3 - Article

AN - SCOPUS:84888201320

VL - 214-215

SP - 217

EP - 227

JO - Geoderma

JF - Geoderma

SN - 0016-7061

ER -