Subsurface physical barriers have been recognized as effective in mitigating seawater intrusion in coastal aquifers, although mainly 2D (cross‐sectional) barrier effects have been considered. In this study, impermeable barriers with finite shore‐parallel lengths are investigated through 3D numerical simulation, thereby extending previous analyses. Two scenarios are considered: (a) barrier‐only and (b) barrier‐well systems; and three available barrier types are analyzed and compared: (1) subsurface dam, (2) cutoff wall, and (3) fully penetrating barrier. Barrier location, length, and height are investigated, and barrier effectiveness is evaluated from seawater volumes, seawater wedge toe positions, and maximum safe pumping rates. In the barrier‐only system, a better performance in preventing seawater intrusion was achieved by cutoff walls rather than subsurface dams. Finite‐length subsurface dams may slightly enhance seawater extent along parts of the coastline that are beyond the dam's length. Cutoff walls performed best when located at relatively small distances from the coast in the barrier‐only system, whereas with a well at 450 m from the shoreline, the subsurface dam located at a critical distance from the sea (i.e., 300 m in the current study) performed optimally (from the tested cases) and was superior to cutoff walls in terms of the maximum safe pumping rate. A fully penetrating barrier outperformed cutoff walls and subsurface dams, as expected. Our investigation indicates that subsurface barrier design should consider the effect of the shore‐parallel length, because barrier benefits may otherwise be significantly overestimated.

We propose a new concept that has the potential to mitigate seawater intrusion and increase the fresh groundwater storage of oceanic islands by creating a less permeable slice along the shoreline. We present a proof‐of‐concept study to examine its effectiveness through analytical and experimental studies. Analytical expressions for calculating the freshwater‐seawater interface location, water table elevation, fresh groundwater volume, and groundwater travel time are presented for both barrier and circular islands, which are found dependent on three different scenarios of interface locations. The analytical solution of the interface location in a barrier island is verified through sand‐tank experiments. Sensitivity analyses based on a simplified conceptual model of St. George Island in Florida, USA, indicate that the fresh groundwater volume monotonically increases with the decrease in the hydraulic conductivity of the coastal less permeable hydrogeologic unit. On the other hand, the increase of the coastal less permeable unit extent leads to an increased fresh groundwater volume. However, when the interface tip is on the aquifer bed of the coastal less permeable unit, a further increase of the less permeable unit extent only slightly increases the fresh groundwater volume, since the interface does not change any more and only the water table is elevated. We demonstrate here that the concept proposed has the potential in increasing the fresh groundwater storage of oceanic islands. Analytical expressions presented can improve our understanding of seawater intrusion in a dual‐unit oceanic island.

Sloping coastal aquifers in reality are ubiquitous and well documented. Steady state sharp‐interface analytical solutions for describing seawater intrusion in sloping confined and unconfined coastal aquifers are developed based on the Dupuit‐Forchheimer approximation. Specifically, analytical solutions based on the constant‐flux inland boundary condition are derived by solving the discharge equation for the interface zone with the continuity conditions of the head and flux applied at the interface between the freshwater zone and the interface zone. Analytical solutions for the constant‐head inland boundary are then obtained by developing the relationship between the inland freshwater flux and hydraulic head and combining this relationship with the solutions of the constant‐flux inland boundary. It is found that for the constant‐flux inland boundary, the shape of the saltwater interface is independent of the geometry of the bottom confining layer for both aquifer types, despite that the geometry of the bottom confining layer determines the location of the interface tip. This is attributed to that the hydraulic head at the interface is identical to that of the coastal boundary, so the shape of the bed below the interface is irrelevant to the interface position. Moreover, developed analytical solutions with an empirical factor on the density factor are in good agreement with the results of variable‐density flow numerical modeling. Analytical solutions developed in this study provide a powerful tool for assessment of seawater intrusion in sloping coastal aquifers as well as in coastal aquifers with a known freshwater flux but an arbitrary geometry of the bottom confining layer.