Separating Weight and Likelihood Tail Effects in Bayesian Neural Groundwater Inversion
Abstract
Groundwater inversion is commonly built on Gaussian assumptions for both network weights and output errors. The assumptions are convenient, but it is unclear which one should be relaxed when hydraulic conductivity is heavy tailed. Here, Gaussian and q-Gaussian distributions were assigned independently to variational weights and the conditional output likelihood in a 2×2 Bayesian neural design. Hydraulic conductivity, porosity, and dispersion were inferred jointly from sparse head and transport observations. Four q regimes (q=1.0, 1.2, 1.5, and 1.8) were evaluated with five independently generated data sets and matched training seeds, giving 130 held-out cases per regime and 500 posterior draws per case. At q=1.5, q-Gaussian weights with a Gaussian likelihood reduced log-conductivity CRPS by 2.85\% (95\% paired CI 0.82--4.88) and calibration error by 3.72 percentage points. Conductivity RMSE was unchanged, while 95\% interval width increased by 9.86\%. At q=1.8, assigning q-Gaussian tails to both components produced no stable accuracy gain and increased interval width by 138.46\%. The Gaussian-generated porosity and dispersion controls showed negligible accuracy changes. Moderate weight-side tails can thus improve distributional prediction, whereas an overly heavy output likelihood trades sharpness for coverage. Separating these two roles provides a clearer basis for uncertainty modeling in groundwater inversion.
Related
- Input-data perturbation in groundwater level prediction: a controlled experiment on accuracy–reliability trade-offs among data augmentation and Bayesian-type uncertainty estimation methods
- Bayesian Joint Velocity and Impedance Inversion via Diffusion Models Conditioned on Common Image Gathers
- Bayesian inversion of multilayer CO2 migration from seismic plume observations using a graph-based finite-rate invasion-percolation model
- Structure-aware Bayesian full waveform inversion via diffusion-regularized Stein variational gradient descent
- Bayesian Geophysical Inversion using Stein Mixture Variational Inference for Multimodal or Complex Solutions
BibTeX
@misc{shen2026separating,
title = {{Separating Weight and Likelihood Tail Effects in Bayesian Neural Groundwater Inversion}},
author = {Ke-Ming Shen},
howpublished = {arXiv},
year = {2026},
month = nov,
doi = {10.48550/arxiv.2611.03669},
eprint = {2611.03669},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2611.03669},
note = {CC BY-NC-ND 4.0},
}
RIS
TY - UNPB
TI - Separating Weight and Likelihood Tail Effects in Bayesian Neural Groundwater Inversion
AU - Ke-Ming Shen
PY - 2026
DA - 2026/11/01
PB - arXiv
DO - 10.48550/arxiv.2611.03669
UR - https://arxiv.org/abs/2611.03669
AB - Groundwater inversion is commonly built on Gaussian assumptions for both network weights and output errors. The assumptions are convenient, but it is unclear which one should be relaxed when hydraulic conductivity is heavy tailed. Here, Gaussian and q-Gaussian distributions were assigned independently to variational weights and the conditional output likelihood in a 2×2 Bayesian neural design. Hydraulic conductivity, porosity, and dispersion were inferred jointly from sparse head and transport observations. Four q regimes (q=1.0, 1.2, 1.5, and 1.8) were evaluated with five independently generated data sets and matched training seeds, giving 130 held-out cases per regime and 500 posterior draws per case. At q=1.5, q-Gaussian weights with a Gaussian likelihood reduced log-conductivity CRPS by 2.85\% (95\% paired CI 0.82--4.88) and calibration error by 3.72 percentage points. Conductivity RMSE was unchanged, while 95\% interval width increased by 9.86\%. At q=1.8, assigning q-Gaussian tails to both components produced no stable accuracy gain and increased interval width by 138.46\%. The Gaussian-generated porosity and dispersion controls showed negligible accuracy changes. Moderate weight-side tails can thus improve distributional prediction, whereas an overly heavy output likelihood trades sharpness for coverage. Separating these two roles provides a clearer basis for uncertainty modeling in groundwater inversion.
ER -