Theory and generation of conditional, scalable sub-Gaussian random fields

M. Panzeri, M. Riva, A. Guadagnini, Shlomo P Neuman

Research output: Contribution to journalArticle

8 Citations (Scopus)

Abstract

Many earth and environmental (as well as a host of other) variables, Y, and their spatial (or temporal) increments, ΔY, exhibit non-Gaussian statistical scaling. Previously we were able to capture key aspects of such non-Gaussian scaling by treating Y and/or ΔY as sub-Gaussian random fields (or processes). This however left unaddressed the empirical finding that whereas sample frequency distributions of Y tend to display relatively mild non-Gaussian peaks and tails, those of ΔY often reveal peaks that grow sharper and tails that become heavier with decreasing separation distance or lag. Recently we proposed a generalized sub-Gaussian model (GSG) which resolves this apparent inconsistency between the statistical scaling behaviors of observed variables and their increments. We presented an algorithm to generate unconditional random realizations of statistically isotropic or anisotropic GSG functions and illustrated it in two dimensions. Most importantly, we demonstrated the feasibility of estimating all parameters of a GSG model underlying a single realization of Y by analyzing jointly spatial moments of Y data and corresponding increments, ΔY. Here, we extend our GSG model to account for noisy measurements of Y at a discrete set of points in space (or time), present an algorithm to generate conditional realizations of corresponding isotropic or anisotropic random fields, introduce two approximate versions of this algorithm to reduce CPU time, and explore them on one and two-dimensional synthetic test cases.

Original languageEnglish (US)
Pages (from-to)1746-1761
Number of pages16
JournalWater Resources Research
Volume52
Issue number3
DOIs
StatePublished - Mar 1 2016

Fingerprint

distribution
test
parameter

Keywords

  • anisotropic random fields
  • conditional simulation
  • generalized sub-Gaussian model
  • non-Gaussian geostatistics
  • non-Gaussian random fields

ASJC Scopus subject areas

  • Water Science and Technology

Cite this

Theory and generation of conditional, scalable sub-Gaussian random fields. / Panzeri, M.; Riva, M.; Guadagnini, A.; Neuman, Shlomo P.

In: Water Resources Research, Vol. 52, No. 3, 01.03.2016, p. 1746-1761.

Research output: Contribution to journalArticle

Panzeri, M. ; Riva, M. ; Guadagnini, A. ; Neuman, Shlomo P. / Theory and generation of conditional, scalable sub-Gaussian random fields. In: Water Resources Research. 2016 ; Vol. 52, No. 3. pp. 1746-1761.
@article{d374461d92c4448eb365bae67d790e9d,
title = "Theory and generation of conditional, scalable sub-Gaussian random fields",
abstract = "Many earth and environmental (as well as a host of other) variables, Y, and their spatial (or temporal) increments, ΔY, exhibit non-Gaussian statistical scaling. Previously we were able to capture key aspects of such non-Gaussian scaling by treating Y and/or ΔY as sub-Gaussian random fields (or processes). This however left unaddressed the empirical finding that whereas sample frequency distributions of Y tend to display relatively mild non-Gaussian peaks and tails, those of ΔY often reveal peaks that grow sharper and tails that become heavier with decreasing separation distance or lag. Recently we proposed a generalized sub-Gaussian model (GSG) which resolves this apparent inconsistency between the statistical scaling behaviors of observed variables and their increments. We presented an algorithm to generate unconditional random realizations of statistically isotropic or anisotropic GSG functions and illustrated it in two dimensions. Most importantly, we demonstrated the feasibility of estimating all parameters of a GSG model underlying a single realization of Y by analyzing jointly spatial moments of Y data and corresponding increments, ΔY. Here, we extend our GSG model to account for noisy measurements of Y at a discrete set of points in space (or time), present an algorithm to generate conditional realizations of corresponding isotropic or anisotropic random fields, introduce two approximate versions of this algorithm to reduce CPU time, and explore them on one and two-dimensional synthetic test cases.",
keywords = "anisotropic random fields, conditional simulation, generalized sub-Gaussian model, non-Gaussian geostatistics, non-Gaussian random fields",
author = "M. Panzeri and M. Riva and A. Guadagnini and Neuman, {Shlomo P}",
year = "2016",
month = "3",
day = "1",
doi = "10.1002/2015WR018348",
language = "English (US)",
volume = "52",
pages = "1746--1761",
journal = "Water Resources Research",
issn = "0043-1397",
publisher = "American Geophysical Union",
number = "3",

}

TY - JOUR

T1 - Theory and generation of conditional, scalable sub-Gaussian random fields

AU - Panzeri, M.

AU - Riva, M.

AU - Guadagnini, A.

AU - Neuman, Shlomo P

PY - 2016/3/1

Y1 - 2016/3/1

N2 - Many earth and environmental (as well as a host of other) variables, Y, and their spatial (or temporal) increments, ΔY, exhibit non-Gaussian statistical scaling. Previously we were able to capture key aspects of such non-Gaussian scaling by treating Y and/or ΔY as sub-Gaussian random fields (or processes). This however left unaddressed the empirical finding that whereas sample frequency distributions of Y tend to display relatively mild non-Gaussian peaks and tails, those of ΔY often reveal peaks that grow sharper and tails that become heavier with decreasing separation distance or lag. Recently we proposed a generalized sub-Gaussian model (GSG) which resolves this apparent inconsistency between the statistical scaling behaviors of observed variables and their increments. We presented an algorithm to generate unconditional random realizations of statistically isotropic or anisotropic GSG functions and illustrated it in two dimensions. Most importantly, we demonstrated the feasibility of estimating all parameters of a GSG model underlying a single realization of Y by analyzing jointly spatial moments of Y data and corresponding increments, ΔY. Here, we extend our GSG model to account for noisy measurements of Y at a discrete set of points in space (or time), present an algorithm to generate conditional realizations of corresponding isotropic or anisotropic random fields, introduce two approximate versions of this algorithm to reduce CPU time, and explore them on one and two-dimensional synthetic test cases.

AB - Many earth and environmental (as well as a host of other) variables, Y, and their spatial (or temporal) increments, ΔY, exhibit non-Gaussian statistical scaling. Previously we were able to capture key aspects of such non-Gaussian scaling by treating Y and/or ΔY as sub-Gaussian random fields (or processes). This however left unaddressed the empirical finding that whereas sample frequency distributions of Y tend to display relatively mild non-Gaussian peaks and tails, those of ΔY often reveal peaks that grow sharper and tails that become heavier with decreasing separation distance or lag. Recently we proposed a generalized sub-Gaussian model (GSG) which resolves this apparent inconsistency between the statistical scaling behaviors of observed variables and their increments. We presented an algorithm to generate unconditional random realizations of statistically isotropic or anisotropic GSG functions and illustrated it in two dimensions. Most importantly, we demonstrated the feasibility of estimating all parameters of a GSG model underlying a single realization of Y by analyzing jointly spatial moments of Y data and corresponding increments, ΔY. Here, we extend our GSG model to account for noisy measurements of Y at a discrete set of points in space (or time), present an algorithm to generate conditional realizations of corresponding isotropic or anisotropic random fields, introduce two approximate versions of this algorithm to reduce CPU time, and explore them on one and two-dimensional synthetic test cases.

KW - anisotropic random fields

KW - conditional simulation

KW - generalized sub-Gaussian model

KW - non-Gaussian geostatistics

KW - non-Gaussian random fields

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

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

U2 - 10.1002/2015WR018348

DO - 10.1002/2015WR018348

M3 - Article

AN - SCOPUS:84960331281

VL - 52

SP - 1746

EP - 1761

JO - Water Resources Research

JF - Water Resources Research

SN - 0043-1397

IS - 3

ER -