Variance-based global sensitivity analysis (e.g., the Sobol' sensitivity index) can be used to identify the important parameters over the entire parameter space. However, one often cannot afford the computational costs of sampling-based approaches in combination with expensive high-fidelity forward models. Reduced-order models (ROM) can substantially accelerate calculation of these sensitivities. However, it is usually difficult to determine what type of ROM should be used and how accurately the ROM represents the high-fidelity model (HFM) results. In this paper, we propose to concurrently use multiple ROMs as a way to assess the robustness of the model-reduction method. Two sets of HFM simulations are needed, one set for building ROMs and the other for validating ROMs. Our goal is to keep the total number of HFM simulations to a minimum. Ideally some of the HFM simulations in the first set can be shared by different ROMs. Based on validation results, the ROMs can be combined with different schemes. We demonstrate that we can achieve the goal by using four different ROMs and still considerably save computational time compared to using traditional HFM simulation for calculating sensitivity indices. We apply the approach to an example problem of a large-scale geological carbon dioxide storage system, in which the objective is to calculate a sensitivity index to identify important parameters. For this problem, the locally best ROM provides better estimates than the weighted average from all ROMs.

High-resolution predictions of land surface hydrological dynamics are desirable for improved investigations of regional- and watershed-scale processes. Direct deterministic simulations of fine-resolution land surface variables present many challenges, including high computational cost. We therefore propose the use of reduced-order modeling techniques to facilitate emulation of fine-resolution simulations. We use an emulator, Gaussian process regression, to approximate fine-resolution four-dimensional soil moisture fields predicted using a three-dimensional surface-subsurface hydrological simulator (PFLOTRAN). A dimension-reduction technique known as “proper orthogonal decomposition” is further used to improve the efficiency of the resulting reduced-order model (ROM). The ROM reduces simulation computational demand to negligible levels compared to the underlying fine-resolution model. In addition, the ROM that we constructed is equipped with an uncertainty estimate, allowing modelers to construct a ROM consistent with uncertainty in the measured data. The ROM is also capable of constructing statistically equivalent analogs that can be used in uncertainty and sensitivity analyses. We apply the technique to four polygonal tundra sites near Barrow, Alaska that are part of the Department of Energy’s Next-Generation Ecosystem Experiments (NGEE)-Arctic project. The ROM is trained for each site using simulated soil moisture from 1998-2000 and validated using the simulated data for 2002 and 2006. The average relative RMSEs of the ROMs are under 1%.

The topography, and the biotic and abiotic parameters are typically upscaled to make watershed-scale hydrologic-biogeochemical models computationally tractable. However, upscaling procedure can produce biases when nonlinear interactions between different processes are not fully captured at coarse resolutions. Here we applied the Proper Orthogonal Decomposition Mapping Method (PODMM) to downscale the field solutions from a coarse (7 km) resolution grid to a fine (220 m) resolution grid. PODMM trains a reduced-order model (ROM) with coarse-resolution and fine-resolution solutions, here obtained using PAWS+CLM, a quasi-3-D watershed processes model that has been validated for many temperate watersheds. Subsequent fine-resolution solutions were approximated based only on coarse-resolution solutions and the ROM. The approximation errors were efficiently quantified using an error estimator. By jointly estimating correlated variables and temporally varying the ROM parameters, we further reduced the approximation errors by up to 20%. We also improved the method's robustness by constructing multiple ROMs using different set of variables, and selecting the best approximation based on the error estimator. The ROMs produced accurate downscaling of soil moisture, latent heat flux, and net primary production with O(1000) reduction in computational cost. The subgrid distributions were also nearly indistinguishable from the ones obtained using the fine-resolution model. Compared to coarse-resolution solutions, biases in upscaled ROM solutions were reduced by up to 80%. This method has the potential to help address the long-standing spatial scaling problem in hydrology and enable long-time integration, parameter estimation, and stochastic uncertainty analysis while accurately representing the heterogeneities.