An open water mechanism of seawater intrusion in the coastal Yucatan aquifer

The Yucatan state, Mexico, encloses a large karstic aquifer, which is confined near the coast. There are places where this confinement is fractured creating springs that discharge to wetlands or shallow lagoons, which are important because of the freshwater ecosystems that develop around them (locally known as Petenes), and because they provide a thermohaline equilibrium in these systems. Previous studies reported that during spring low tide, the spring becomes a sink, reversing the flow in the Peten channel connecting the lagoon to the spring. Potentially, this inversion of the flow can cause the intrusion of brackish – and sometimes hypersaline - water from the lagoon into the aquifer, which can also have ecological implications. The flow reversal was observed in two sampling campaigns and it occurred a total of 16 times. In the last sampling campaign, it was confirmed that the spring becomes a sink, with measured water velocities of about 0.37 m/s at the channel end close to the spring. Preliminary results suggest that the flow reversal is controlled by the discharge from the aquifer, the confinement of the aquifer, the geometry of the lagoon and the inlet channel, as well as the sea tides. A simple 1D analytical solution was used to explain this phenomenon, which describes the hydrodynamic interaction between the lagoon and the aquifer by solving the heat equation, which can be used to simulate flow in both systems, a confined aquifer and a highly frictional lagoon. Therefore, the same solution was used to simulate the gradients observed between the confined aquifer and the lagoon. The solution was implemented as a Python library using Fourier series, and has the advantage that it can be used with more general boundary conditions and finite length systems. Finally, we present coefficients that can be used to determine under which conditions the flow reverses at the lagoon.