River Seepage Conductance in Large-Scale Regional Studies

Surface and groundwater interactions are essential for managing water resources and legal water rights. Prior research has shown that numerical models often use large grid sizes, which limit the precision of flow exchange simulations. Established methods rely on empirical coefficients calibrated to match observed data, but these coefficients lack a strong physical basis. This gap motivated the search for an analytical approach that retains accuracy even with coarse grids. No prior work had resolved how to incorporate physical parameters into leakance coefficients without extensive recalibration. That uncertainty drove the development of a new method to estimate the leakance coefficient using physical principles. This approach aims to reduce reliance on arbitrary calibration while maintaining model efficiency. Understanding the physical basis of leakance is crucial for improving large-scale groundwater models.

The goal of this study is to provide an analytical framework for estimating the leakance coefficient in river-aquifer interactions. The researchers propose an alternative to the commonly used empirical leakance coefficient, which is typically determined through calibration. The specific problem addressed is the lack of a physically grounded method for calculating seepage discharge in large-scale models. This approach seeks to maintain accuracy even with coarse grid sizes. The motivation stems from the need to reduce model recalibration when grid resolution changes. The study focuses on translating physical parameters into a leakance coefficient. It aims to simplify model adjustments without sacrificing precision. This work is intended to improve the reliability of groundwater simulations at regional scales.

The researchers employed an analytical solution to derive an empirical leakance coefficient. They used a Cauchy boundary condition to model the exchange between rivers and aquifers. The method assumes a water-table aquifer treated as a single layer in the model. Physical principles were used to relate the empirical coefficient to exact conductance values. The approach incorporates factors like wetted perimeter and river penetration. Anisotropy and clogging layers were included without increasing computational cost. The method avoids recalibration by linking the coefficient to physical parameters. The solution is designed for use in models like MODFLOW with minimal grid refinement.

The study found that an analytical solution can replace the empirical leakance coefficient in large-scale models. The derived coefficient matches the accuracy of fine-grid simulations while using coarse grids. The method allows for physical interpretation of the leakance coefficient. Normalized wetted perimeter and river penetration were key variables in the calculation. The model maintained accuracy even when grid size changed. The presence of a clogging layer was easily incorporated into the framework. Anisotropy was accounted for without increasing computational demand. The approach eliminates the need for recalibration when grid resolution is altered.

The authors propose that the analytical leakance coefficient provides a physically based alternative to traditional calibration methods. This method supports accurate simulations even with coarse grid sizes. The study suggests that the approach improves model reliability in large-scale groundwater studies. The researchers propose that the inclusion of physical parameters enhances model transparency. The findings suggest that the method is computationally efficient and adaptable. The approach allows for easier modification of the river coefficient in models like MODFLOW. The study concludes that the method reduces the need for recalibration when grid sizes are adjusted. The authors suggest that this framework can be applied in regional groundwater modeling with minimal effort.