Modeling of phosphate flux induced by flood resuspension on a macrotidal estuarine mudflat (Seine, France)

Coastal marine sediments can be either major scrubbers or eutrophication contributors to surface waters. Standard methods for direct measurement of nutrient fluxes at the sediment-water interface do not consider hydrodynamic forcing although several ex-situ studies suggest that sediment resuspension can dramatically increase dissolved fluxes. We provide a new model to quantify dissolved phosphate (PO43−) resuspension flux (JR) based on physical representation of its identified components: diffusion stimulation by exposure of deeper sediment layer with higher PO43− concentration in the porewater (JD), pore water mixing with overlying water (JM) and net adsorption/desorption from suspended sediments (JK). This approach was applied to field data from a Seine intertidal mudflat periodically submitted to millimetric erosion. On a tidal scale, the model output reveals a JR of 272.3 ± 360.0 μmol m−2 h−1 (± 52% from parameter uncertainty), well above flux calculated by application of Fick's first law (0.15 ± 0.85 μmol m−2 h−1) or by ex situ core incubation (40.8 μmol m−2 h−1). Iron bound phosphorus within suboxic layers buffers PO43− concentrations in superficial sediments leading to negligible contributions of JD and JM to total fluxes. Conversely, JK appears to be the main exchange pathway, even though improvements in turbidity measurement would allow this term to be defined more precisely. Correction required to enhance and control model robustness are described. These results show the importance of considering the dissolved PO43− resuspension flux in dynamic environments.