ABSTRACT:

Coastal aquifers commonly exhibit heterogeneity in hydraulic conductivity, a critical factor controlling groundwater extraction. Most existing analytical solutions for groundwater withdrawal in coastal regions assume homogeneous or vertically stratified aquifers, neglecting anthropogenic modifications such as land reclamation and subsurface construction that introduce horizontal heterogeneity. This HydroShare resource contains the numerical models and MATLAB codes developed to support the study of “Maximum Sustainable Pumping Rate in Coastal Aquifers with Landward Horizontal Heterogeneity: A Two-Zone Semi-Analytical Solution”. The archived materials include numerical groundwater flow and solute transport models (case setting: L = 400 m, x_w = 300 m, K1 = 2 m/d, K2 = 1 m/d, αL = 10.0 m, and αT = 0.50 m) and a MATLAB script for calculating the maximum sustainable pumping rate. These materials are provided to facilitate reproducibility of the numerical experiments and semi-analytical calculations presented in the study and to support further investigation of groundwater management in heterogeneous coastal aquifer systems.

ABSTRACT:

This HydroShare resource contains the analytical and numerical modeling files developed to investigate the maximum pumping rate in a horizontally heterogeneous coastal aquifer affected by land reclamation. The aquifer is represented as a two-zone unconfined system consisting of the original land and a reclaimed area with contrasting hydraulic conductivities. Inland and seaward constant-head boundaries and uniform areal recharge are considered. The resource includes MATLAB source code implementing the analytical solution based on discharge-potential theory, superposition, the Ghyben–Herzberg relation, and the Fourier cosine transform. The MATLAB scripts also determine the limiting pumping rate using a bisection procedure in which the pumping-induced stagnation point coincides with the toe of the seawater wedge. Numerical model files developed with SEAWAT-2000 are provided for variable-density groundwater-flow and solute-transport simulations used to validate the analytical results. The resource additionally includes model parameters, simulation scenarios, numerical outputs, and post-processing scripts used to compare analytical and numerical estimates. These files support reproduction of the principal results and facilitate further investigation of the effects of reclamation length, hydraulic-conductivity contrast, well location, inland head conditions, and precipitation recharge on groundwater pumping limits in coastal aquifers.