Seepage boundary conditions are commonly used in groundwater simulations to allow groundwater to discharge at the upper surface of the model when groundwater head exceeds atmospheric pressure. However, the extent and transient behavior of the seepage zone is often unknown a priori and is difficult to predict. A mathematical description of the boundary condition is straightforward, such that head is equivalent to elevation only when groundwater flow indicates a seepage condition, which is a mixed conditional Dirichlet and Neumann boundary condition. This standard representation of the boundary condition has been successfully implemented and applied in a real-world context by most groundwater models. However, it is rarely reported that convergence is only guaranteed when both the efflux and zero pressure conditions are simultaneously satisfied, often requiring unnecessarily small timestep sizes, which results in low computational efficiency. This study suggests a continuous differentiable equation as an alternative to model the seepage boundary. The new formulation is derived by analogy to the first-order exchange equation, which is commonly used to represent the interactions between surface water and groundwater flow in integrated hydrologic simulations. The results of this study suggest that mixed Dirichlet and Neumann boundary conditions can be effectively converted into a Robin boundary condition, which is a head-dependent flux condition that incorporates appropriate physical considerations. This new approach has the potential to significantly improve the accuracy and efficiency of groundwater flow simulations and can help to advance the understanding of subsurface hydrology.

This study proposes the utilization of municipal well records as an alternative dataset for large-scale heterogeneity characterization of hydraulic conductivity () and specific storage () using hydraulic tomography (HT). To investigate the performance of HT and the feasibility of utilizing municipal well records, a three-dimensional aquifer/aquitard system is constructed and synthetic groundwater flow and solute transport experiments are conducted to generate data for inverse modeling and validation of results. In particular, we simultaneously calibrate four groundwater models with varying parameterization complexity using five datasets consisting of different time durations and periods. Calibration and validation results are qualitatively and quantitatively assessed to evaluate the performance of investigated models. The estimated and tomograms from different model cases are also validated through the simulation of independently conducted pumping tests and conservative solute transport. Our study reveals that: 1) the HT analysis of municipal well records is feasible and yields reliable heterogeneous and distributions where drawdown records are available; 2) accurate geological information is of critical importance when data density is low and should be incorporated for geostatistical inversions; 3) the estimated and tomograms from the geostatistical model with geological information are capable in providing robust predictions of both groundwater flow and solute transport. Overall, this synthetic study provides a general framework for large-scale heterogeneity characterization using HT through the interpretation of municipal well records, and provides guidance for applying this concept to field problems.