Modeling of Groundwater Pollution by Former Seawater in Submarine Groundwater Discharge Driven by Saltwater Intrusion

Coastal aquifers are usually vulnerable to contamination by saltwater intrusion. The degree of contamination depends on the amount of former seawater intruding the aquifer as submarine groundwater discharge (SGD). A three-dimensional numerical SEAWAT model was developed to provide insights to the responses of saltwater intrusion and amount of seawater portion of SGD, to groundwater withdrawal from four wells at different spatial locations. The results showed that saltwater wedge encroachment varies with a withdrawal rate and well location. Increasing withdrawal rates from 0.165 $m^3$/s to 0.53166 $m^3$/s resulted in a noticeable increase in wedge encroachment into the aquifer from 589 m to 1319 m, respectively, regardless of the distance from the coastline. However, higher withdrawals from the wells closer to the coastline caused the hydraulic head near the seaside to drop below the terrestrial head, causing the saltwater wedge to be pushed back toward the sea. Simulations also showed that a coastal well might act as a hydraulic barrier that prevents the brackish zone from moving further inland, which is critically important in terms of groundwater management. Seawater contributed from 5% to above 33% in SGD. Although the withdrawal rates from Well 3 were 10 to 20 times larger than the Well 2 rates, the seawater contribution associated with Well 3 was about the same. This is attributed to the horizontal vicinity of Well 2 to the seacoast boundary compared to Well 3 even though Well 3 is deeper than Well 2, demonstrating that it is the horizontal location that governs the amount of abstracted former seawater.