Managing fragile island freshwater resources requires identifying pumping strategies that trade off the financial cost of groundwater supply against controlling the seawater intrusion (SWI) associated with aquifer pumping. In this work, these tradeoffs are investigated through a sensitivity analysis conducted in the context of an optimization formulation of the groundwater management problem, which aims at minimizing the groundwater supply operation cost associated with groundwater pumping and desalination treatment, subject to constraints on SWI control, as quantified by the water table drawdown over the well (∆s), the reduction in freshwater volume (∆FV) in the aquifer, or the salt mass increase (∆SM) in the aquifer. This study focuses on a simplified two‐dimensional model of the San Salvador Island aquifer (Bahamas). Pumping strategies are characterized by the distance of the pumping system from the shoreline (WL), the abstraction screen depth (D) and overall pumping rate (Q), constituting the decision variables of the optimization problem. We investigate the impacts of pumping strategies on the operation cost, ∆s, ∆FV and ∆SM. Findings indicate increasing D or decreasing WL reduces ∆s, ∆FV and ∆SM, thus preserving the aquifer hydrogeologic stability, but also leads to extracting saltier groundwater, thus increasing the water treatment requirements, which have a strong impact on the overall groundwater supply cost. From a financial perspective, groundwater abstraction near the island center and at shallow depths seems the most convenient strategy. However, the analysis of the optimization constraints reveals that strategies where the pumping system approaches the island center tend to cause more severe SWI, highlighting the need to trade off groundwater supply cost against SWI control.

Based on the results of the test among the joint roughness coefficient (JRC) of rock fracture, mechanical aperture, and hydraulic aperture proposed by Barton, this paper deduces and proposes a permeability coefficient formula of single fracture stress‐seepage coupling considering microroughness by the introduction of effect variables considering the microparticle size and structural morphology of facture surface. Quasi‐sandstone fracture of different particle size is made by the laboratory test, and the respective modification is made on the coupled shear‐seepage test system of JAW‐600 rock. Under this condition, the laboratory test of stress‐seepage coupling of fracture of different particle size is carried out. The test results show that, for the different particle‐sized fracture surface of the same JRC, the permeability coefficient is different, which means the smaller particle size, the smaller permeability coefficient, and the larger particle size, the larger permeability coefficient; with the increase of cranny hydraulic pressure, the permeability coefficient increases exponentially, and under the same cranny hydraulic pressure, there is relation of power function between the permeability coefficient and normal stress. Meanwhile, according to the theoretical formula, the microroughness coefficient of the fractures with different particle size is obtained by the calculation, and its accuracy and validity are verified by experiments. The theoretical verification values are in good agreement with the measured values.

Sea-level rise due to global warming and over-exploitation of groundwater resources in coastal areas will induce seawater intrusion into inland groundwater which is leading to the migration of the transition zone, and is affecting the security of regional social economy and water resources. Taking the blue economic zone of Shandong Peninsula as the research background, selecting the coastal area of Longkou and Zhaoyuan as the study area, this paper firstly depicts the groundwater flow field pattern in the research area, and then carries out the quantitative analysis for the migration patterns of the transition zone under changing groundwater levels. The results demonstrate that the width of sea-land transitional zone in the research area is about 1.5–4.5 km. The groundwater level is higher, the migration speed of the sea-land transitional zone is smaller. The results are of significance to study the migration of the sea-land transitional zone in the blue economic zone of Shandong peninsula.

In groundwater pumping optimization (GPO), offline-trained data-driven surrogates can be used to replace numerical-intensive simulators in order to save computing time. The traditional offline training approach involves building surrogates prior to optimization, fitting training datasets that cover the input space uniformly or randomly, which can prove inefficient due to the potential oversampling of low-gradient areas and under-sampling of high-gradient areas. This study proposes an offline machine-learning (ML) algorithm that ranks candidate training points by scoring them based on their distance to the closest training point and on the local gradient of the surrogate estimate and then choosing the highest-rank point. This method is applied to develop surrogates for solving a two-objective GPO problem formulated on a three-dimensional (3D) island aquifer, using hydrogeological conditions representative of San Salvador Island, Bahamas. The objectives are to minimise the supply cost (f_OC) resulting from groundwater pumping and desalination and maximise fresh groundwater supply (Q_p), subject to constraints on seawater intrusion (SWI) control expressed in terms of aquifer drawdown Δs at pumping locations and aquifer salt mass increase ΔSM. Gaussian Process (GP) is the technique applied to construct surrogates of objectives and constraints, alongside the estimation of uncertainties. Using GP models, it is possible to estimate the probability of “Pareto optimality” for each pumping scheme by Monte Carlo simulation. Pareto optimal pumping schemes (POPS) are then characterized by a probability of occurrence, which can be verified by numerical simulation. The GP training strategy's effectiveness in generating POPS is compared to traditional training approaches, showing that such a strategy can efficiently identify reliable POPS.

Sea-level rise (SLR) and regional precipitation pattern change cause island subsurface freshwater, typically shaped like a thin lens, to be at higher risk of contamination from seawater intrusion (SWI). Installing a cutoff wall is considered a feasible strategy for protecting coastal fresh groundwater from SWI. However, the performance of the cutoff wall in managing freshwater lens (FWL) development and mitigating SWI into island aquifers under SLR and aquifer recharge (RCH) fluctuations remains inadequately quantified. This study investigates how water table elevation (WTE), FWL depth, thickness, and SWI extent, measured by aquifer salt mass and freshwater volume, in an island aquifer equipped with cutoff walls, respond to SLR and RCH fluctuations. It focuses on a two-dimensional, variable-density island groundwater simulation model based on hydrogeological conditions of San Salvador Island, Bahamas. The results demonstrate that RCH critically influences cutoff wall effectiveness for FWL development and SWI mitigation, with higher RCH amplifying gains in WTE, FWL metrics, freshwater storage, and aquifer salt removal, but this influence diminishes with wall depth increasing. SLR elevates WTE in a stable manner associated with its magnitude but negligibly affects the cutoff wall performance in FWL enhancement and SWI mitigation. Under simultaneous SLR and RCH fluctuations, SLR can offset the WTE reduction caused by reduced RCH, but the joint effects of SLR and RCH on FWL metrics, freshwater storage and aquifer salt removal align with their individual impacts. Moreover, cutoff walls are more efficient in low-RCH settings, yielding greater relative improvements in FWL development and SWI mitigation per unit wall depth increase.

Subsurface freshwater in oceanic islands is typically shaped like a thin lens due to limited land area and recharge, often the primary freshwater source for local communities and highly vulnerable to seawater intrusion (SWI). Freshwater injection (FI) is considered as a feasible strategy for mitigating SWI in coastal aquifers. However, its transient effectiveness for freshwater lens (FWL) development and SWI mitigation in island aquifers and how the design parameters like FI depth, intensity, duration and injectant concentration affect its performance remain poorly understood. To address this, this study employs a two-dimensional, variable-density island groundwater model to simulate the transient responses of FWL development and SWI mitigation to various FI patterns. Five indicators are developed for comprehensive evaluation, including (1) freshwater recovery efficiency (FRE), and the relative changes in (2) average water table elevation (WTE), (3) FWL depth, (4) FWL volume, and (5) total aquifer salt mass. Results reveal FI universally raises average WTE, expands FWL dimensions, and promotes aquifer desalinization. Injection intensity is the primary driver of WTE rises and salt mass reduction, with higher intensities consistently yielding greater WTE rises and salt mass reductions. Deeper injection within the mixing zone increases FWL depth, but reduces the net gain in FWL volume. Moreover, early-stage FI is highly efficient for expanding FWL volume, often yielding FRE values above 100%, but FRE converges toward zero over time as the system moves toward a new hydrodynamic equilibrium, returning diminishing marginal benefits for long-term FI.