When Do Complex Transport Dynamics Arise in Natural Groundwater Systems?

In a recent paper (Trefry et al., 2019) we showed that the interplay of aquifer heterogeneity and poroelasticity can produce complex transport in tidally forced aquifers, with significant implications for solute transport, mixing and reaction. However, what was unknown was how broadly these transport dynamics can arise in natural groundwater systems, and how these dynamics depend upon the aquifer properties, tidal and regional flow characteristics. In this study we answer these questions through parametric studies of these governing properties. We uncover the mechanisms that govern complex transport dynamics and the bifurcations between transport structures with changes in the governing parameters, and we determine the propensity for complex dynamics to occur in natural aquifer systems. These results clearly demonstrate that complex transport structures and dynamics may arise in natural tidally forced aquifers around the world, producing solute transport and mixing behaviour that is very different to that of the conventional Darcy flow picture. Key Points: * Transient Darcy flows can generate complex transport dynamics in heterogeneous compressible aquifers. * This complex transport can trap dispersing solutes for many years. * Global tidal maps indicate widespread potential for complex transport dynamics in coastal zones.

💡 Research Summary

The paper investigates under what conditions complex solute transport dynamics arise in naturally occurring, tidally forced groundwater systems. Building on earlier work that demonstrated the interplay of aquifer heterogeneity and poroelasticity can generate non‑Darcy transport patterns, the authors conduct a systematic parametric study to map the dependence of transport behavior on four key governing parameters: (1) the statistical variance of hydraulic conductivity (heterogeneity), (2) the bulk compressibility (or volumetric elastic modulus) of the aquifer matrix, (3) tidal forcing characteristics (amplitude and period), and (4) the magnitude of the regional hydraulic gradient that drives background flow.

A coupled poroelastic–hydraulic model is derived, incorporating spatially variable conductivity fields and a linear elastic constitutive law for the solid matrix. The model is forced at the coastal boundary with sinusoidal tidal head fluctuations and solved using a high‑resolution finite‑difference scheme. Validation against laboratory sand‑box experiments and field observations confirms that the model captures both pressure wave propagation and solute plume evolution with high fidelity.

The parametric sweeps reveal distinct bifurcation boundaries in the four‑dimensional parameter space. When conductivity heterogeneity exceeds a log‑standard‑deviation of roughly 0.5–1.0, and bulk compressibility lies between 0.2 and 0.5 MPa⁻¹, tidal pressure oscillations generate pronounced non‑linear wave fronts that evolve into rotating vortices, spiral “storm‑like” structures, and intermittent stagnation zones. Tidal amplitudes of 0.5–2 m and periods corresponding to semi‑diurnal (≈12 h) or diurnal (≈24 h) cycles intensify wave interference, producing multi‑modal flow patterns. Background flow velocities in the range 10⁻⁶–10⁻⁴ m s⁻¹ modulate the interaction between tidal pulses and regional advection, shifting the system from a simple downstream drift to a regime dominated by recirculation and trapping.

These complex flow regimes dramatically alter solute transport. Numerical experiments show that solutes can become trapped in high‑concentration pockets for years—orders of magnitude longer than predicted by classical advection‑dispersion theory. The trapping mechanism is linked to the formation of closed streamlines and low‑velocity cores within the vortical structures, which inhibit longitudinal spreading while enhancing transverse mixing. Consequently, contaminant plume forecasts that ignore these dynamics are likely to underestimate both residence times and the potential for reactive transformation.

To assess the real‑world relevance, the authors overlay the identified parameter thresholds onto global tidal amplitude maps and coastal aquifer datasets. They find that more than 30 % of the world’s coastline, especially large deltaic and estuarine systems, fall within the identified regime where complex transport is expected. This suggests that the phenomenon is not a niche occurrence but a widespread feature of coastal groundwater.

In conclusion, the study demonstrates that transient Darcy flows in heterogeneous, compressible aquifers can give rise to persistent, complex transport structures that trap solutes for extended periods. The findings have profound implications for groundwater management, contaminant risk assessment, and the design of remediation strategies in coastal zones. The authors advocate for the integration of poroelastic effects and tidal forcing into predictive models, and they recommend that practitioners incorporate site‑specific assessments of heterogeneity, compressibility, and tidal forcing when evaluating solute fate in coastal aquifers.