Characterizing the impact of model error in hydrologic time series recovery inverse problems
Abstract
Hydrologic models are commonly oversmoothed relative to reality, owing to computational limitations and to the difficulty of obtaining accurate highresolution information. When used in an inversion context, such models may introduce systematic biases which cannot be encapsulated by an unbiased “observation noise” term of the type assumed by standard regularization theory and typical Bayesian formulations. Despite its importance, model error is difficult to encapsulate systematically and is often neglected. In this paper, model error is considered for an important class of inverse problems that includes interpretation of hydraulic transients and contaminant source history inference: reconstruction of a time series that has been convolved against a transfer function (i.e., impulse response) that is only approximately known. Using established harmonic theory along with two results established here regarding triangular Toeplitz matrices, upper and lower error bounds are derived for the effect of systematic model error on time series recovery for both welldetermined and overdetermined inverse problems. It is seen that use of additional measurement locations does not improve expected performance in the face of model error. A Monte Carlo study of a realistic hydraulic reconstruction problem is presented, and the lower error bound is seen informative about expected behavior. Finally, a possiblemore »
 Authors:
 Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
 Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Univ. of Texas, Austin, TX (United States). Computational Hydraulics Group (CHG). Inst. for Computational Engineering and Sciences (ICES)
 Publication Date:
 Research Org.:
 Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
 Sponsoring Org.:
 USDOE
 OSTI Identifier:
 1412889
 Report Number(s):
 LAUR1628825
Journal ID: ISSN 03091708
 Grant/Contract Number:
 AC5206NA25396
 Resource Type:
 Journal Article: Accepted Manuscript
 Journal Name:
 Advances in Water Resources
 Additional Journal Information:
 Journal Name: Advances in Water Resources; Journal ID: ISSN 03091708
 Publisher:
 Elsevier
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS
Citation Formats
Hansen, Scott K., He, Jiachuan, and Vesselinov, Velimir V.. Characterizing the impact of model error in hydrologic time series recovery inverse problems. United States: N. p., 2017.
Web. doi:10.1016/j.advwatres.2017.09.030.
Hansen, Scott K., He, Jiachuan, & Vesselinov, Velimir V.. Characterizing the impact of model error in hydrologic time series recovery inverse problems. United States. doi:10.1016/j.advwatres.2017.09.030.
Hansen, Scott K., He, Jiachuan, and Vesselinov, Velimir V.. 2017.
"Characterizing the impact of model error in hydrologic time series recovery inverse problems". United States.
doi:10.1016/j.advwatres.2017.09.030.
@article{osti_1412889,
title = {Characterizing the impact of model error in hydrologic time series recovery inverse problems},
author = {Hansen, Scott K. and He, Jiachuan and Vesselinov, Velimir V.},
abstractNote = {Hydrologic models are commonly oversmoothed relative to reality, owing to computational limitations and to the difficulty of obtaining accurate highresolution information. When used in an inversion context, such models may introduce systematic biases which cannot be encapsulated by an unbiased “observation noise” term of the type assumed by standard regularization theory and typical Bayesian formulations. Despite its importance, model error is difficult to encapsulate systematically and is often neglected. In this paper, model error is considered for an important class of inverse problems that includes interpretation of hydraulic transients and contaminant source history inference: reconstruction of a time series that has been convolved against a transfer function (i.e., impulse response) that is only approximately known. Using established harmonic theory along with two results established here regarding triangular Toeplitz matrices, upper and lower error bounds are derived for the effect of systematic model error on time series recovery for both welldetermined and overdetermined inverse problems. It is seen that use of additional measurement locations does not improve expected performance in the face of model error. A Monte Carlo study of a realistic hydraulic reconstruction problem is presented, and the lower error bound is seen informative about expected behavior. Finally, a possible diagnostic criterion for blind transfer function characterization is also uncovered.},
doi = {10.1016/j.advwatres.2017.09.030},
journal = {Advances in Water Resources},
number = ,
volume = ,
place = {United States},
year = 2017,
month =
}

This presentation describes some methods for modeling error in geophysical inverse problems.

Interpolation error operator: time series error detection and correction
The interpolation error operator is an effective tool for the detection of digital time series data errors characterized by discontinuities in the function itself or its loworder derivatives. The generalized interpolation operator of multiplicity n simultaneously interpolates n consecutive terms. This operator is effective in correcting short lengths of erroneous data. 1 figure, 2 tables. 
Model reduction of dynamical systems by proper orthogonal decomposition: Error bounds and comparison of methods using snapshots from the solution and the time derivatives [Proper orthogonal decomposition model reduction of dynamical systems: error bounds and comparison of methods using snapshots from the solution and the time derivatives]
In this study, we consider two proper orthogonal decomposition (POD) methods for dimension reduction of dynamical systems. The first method (M1) uses only time snapshots of the solution, while the second method (M2) augments the snapshot set with timederivative snapshots. The goal of the paper is to analyze and compare the approximation errors resulting from the two methods by using error bounds. We derive several new bounds of the error from POD model reduction by each of the two methods. The new error bounds involve a multiplicative factor depending on the time steps between the snapshots. For method M1 themore » 
Inverse Modeling of Hydrologic Parameters Using Surface Flux and Runoff Observations in the Community Land Model
This study demonstrates the possibility of inverting hydrologic parameters using surface flux and runoff observations in version 4 of the Community Land Model (CLM4). Previous studies showed that surface flux and runoff calculations are sensitive to major hydrologic parameters in CLM4 over different watersheds, and illustrated the necessity and possibility of parameter calibration. Two inversion strategies, the deterministic leastsquare fitting and stochastic MarkovChain MonteCarlo (MCMC)  Bayesian inversion approaches, are evaluated by applying them to CLM4 at selected sites. The unknowns to be estimated include surface and subsurface runoff generation parameters and vadose zone soil water parameters. We find thatmore »