Coupling basin- and site-scale inverse models of the Española aquifer

Large-scale models are frequently used to estimate fluxes to small-scale models. The uncertainty associated with these flux estimates, however, is rarely addressed. We present a case study from the Española Basin, northern New Mexico, where we use a basin-scale model coupled with a high-resolution, nested site-scale model. Both models are three-dimensional and are analyzed by codes FEHM and PEST. Using constrained nonlinear optimization, we examine the effect of parameter uncertainty in the basin-scale model on the nonlinear confidence limits of predicted fluxes to the site-scale model. We find that some of the fluxes are very well constrained, while for others there is fairly large uncertainty. Site-scale transport simulation results, however, are relatively insensitive to the estimated uncertainty in the fluxes. We also compare parameter estimates obtained by the basin- and site-scale inverse models. Differences in the model grid resolution (scale of parameter estimation) result in differing delineation of hydrostratigraphic units, so the two models produce different estimates for some units. The effect is similar to the observed scale effect in medium properties owing to differences in tested volume. More important, estimation uncertainty of model parameters is quite different at the two scales. Overall, the basin inverse model resulted in significantly lower estimates of uncertainty, because of the larger calibration dataset available. This suggests that the basin-scale model contributes not only important boundary condition information but also improved parameter identification for some units. Our results demonstrate that caution is warranted when applying parameter estimates inferred from a large-scale model to small-scale simulations, and vice versa.