We analyze the statistical scaling properties of vertical and horizontal increments in soil and sediment texture data measured to a depth of 15 meters over an area of 3600 m2 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.

Anisotropic Statistical Scaling of Soil and Sediment Texture in a Stratified Deep Vadose Zone near Maricopa, Arizona.

GUADAGNINI, ALBERTO;RIVA, MONICA
2014-01-01

Abstract

We analyze the statistical scaling properties of vertical and horizontal increments in soil and sediment texture data measured to a depth of 15 meters over an area of 3600 m2 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.
2014
soil and sediment texture; random fields; anisotropy; truncated fractional Brownian motion; scaling; antipersistent spatial variability.
File in questo prodotto:
File Dimensione Formato  
2014 - Guadagnini et al (Geoderma - Maricopa GEODER11469).pdf

Accesso riservato

: Post-Print (DRAFT o Author’s Accepted Manuscript-AAM)
Dimensione 997.09 kB
Formato Adobe PDF
997.09 kB Adobe PDF   Visualizza/Apri
Geoderma - Maricopa GEODER11469_11311-761456_Guadagnini.pdf

accesso aperto

: Post-Print (DRAFT o Author’s Accepted Manuscript-AAM)
Dimensione 789.62 kB
Formato Adobe PDF
789.62 kB Adobe PDF Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11311/761456
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 28
  • ???jsp.display-item.citation.isi??? 26
social impact