Impact of mineral reactions and surface complexation on the transport of dissolved species in a subterranean estuary: Application of a comprehensive reactive transport modeling approach

Subterranean estuaries (STE) are hotspots of biogeochemical reactions. Here, dissolved constituents in waters of terrestrial and marine origin are transformed before they discharge to the coastal oceans. The involved biogeochemical reactions are complex and non-linear, calling for the application of numerical reactive transport modeling (RTM) to improve the process understanding. The aim of this study was to assess the roles of organic matter degradation and coupled secondary mineral reactions for the fate of dissolved species in STEs of sandy beaches. A comprehensive RTM approach was applied for this purpose, accounting for the effects of ion activities, pH, pe, redox reactions, mineral equilibria (calcite, goethite, siderite, iron sulfide, hydroxyapatite and vivianite) as well as surface complexation. Results show that the STE biogeochemistry and associated species fluxes are very sensitive to the assumed reaction network. For example, inorganic carbon and pH were largely controlled by calcite and siderite dynamics, and dissolved Fe²⁺ and HS^- were precipitated as goethite, siderite and/or iron sulfides. Moreover, PO₄^3- concentrations were affected by both the formation of vivianite or hydroxyapatite as well as surface complexation. This work helped to establish the relative importance of some of the major biogeochemical processes in the STE. However, further field studies are needed to understand which processes play a role in real-world STEs, including an exploration of the deep subsurface of STEs. Such field-based observations will improve our conceptual process understanding, which is key to developing well-constrained RTMs.