Dataset to the article 'Redox-zoning in high-energy subterranean estuaries as a function of storm floods, temperatures, seasonal groundwater recharge and morphodynamics'

This repository contains model input files, and pre- and post-processing scripts for the research article:

Janek Greskowiak, Stephan L. Seibert, Vincent E.A. Post, Gudrun Massmann (2023), Redox-zoning in high-energy subterranean estuaries as a function of storm floods, temperatures, seasonal groundwater recharge and morphodynamics, Estuarine, Coastal and Shelf Science, https://doi.org/10.1016/j.ecss.2023.108418

Data:

The folder Laserascannerdata_top_and_hydr_heads contains five cross-shore topography profiles of a high energy beach on Spiekeroog Island, Germany, including a python-script that assembles them to a daily times-series over one year and that calculates the tide-averaged hydraulic heads in the intertidal zone as detailed in Greskowiak and Massmann (2021), The impact of morphodynamics and storm floods on pore water flow and transport in the subterranean estuary, Hydrological Processes, 35:e14050, https://doi.org/10.1002/hyp.14050

Model input:

The file Dynamod_1_redox_basecase.zip contains a folder with the model input textfiles (SEAWAT and PHT3D) for the case that considers all dynamic influence factors investigated in the paper, i.e., seasonal meteoric groundwater recharge, stormfloods, temperature-dependence of reaction rates and morphodynamics.

The files pht3d_ph.dat and pht3d_datab.dat are PHT3D model specific files defining the reacants and reactions, and have to be copied into the model folder before PHT3D is started.

Run models:

First the SEAWAT model needs to be run with the name file 'Dynamod_1_redox_basecase.nam'. This creates the mt3d link file 'mt3d_link.ftl'. After that, PHT3D has to be started with the name file 'pht3d.nam'

Note that running the models will generate 22 Gybte output.

Pre-processing:

All model input files were generated with the python-script Dynamod_1_redox_basebase.py using Flopy, a python-based groundwater modelling user-interface:
https://www.usgs.gov/software/flopy-python-package-creating-running-and-post-processing-modflow-based-models

Post-processing:

For the figures 3,4 and S1, S2 in the paper, the corresponding python scripts are provided in this database. Note that with respect to Figures 5 and S3, only the results for the case with all dynamic influence factors (Figure S3_h) are being plotted.

Other model output:
Animation_A1_Dynamod_redox.mp4 is a video animation showing the dynamic modelled salinity and temperature distribution, and redox zoning in the subterranean estuary all dynamic influence factors investigated in the paper, i.e., seasonal meteoric groundwater recharge, stormfloods, temperature-dependence of reaction rates and morphodynamics.