Incorporating temporal and spatial variability of salt-marsh foraminifera into sea-level reconstructions

Foraminifera from salt-marsh environments have been used extensively in quantitative relative sea-level reconstructions due to their strong relationship with tidal level. However, the influence of temporal and spatial variability of salt-marsh foraminifera on quantitative reconstructions remains unconstrained. Here, we conducted a monitoring study of foraminifera from four intertidal monitoring stations in New Jersey from high marsh environments over three years that included several extreme weather (temperature, precipitation, and storm surge) events. We sampled four replicates from each station seasonally (four times per year) for a total of 188 samples. The dead foraminiferal assemblages were separated into four site-specific assemblages. After accounting for systematic trends in changes in foraminifera over time among stations, the distribution of foraminiferal assemblages across monitoring stations explained ~87% of the remaining variation, while ~13% can be explained by temporal and/or spatial variability among the replicate samples. We applied a Bayesian transfer function to estimate the elevation of the four monitoring stations. All samples from each station predicted an elevation estimate within a 95% uncertainty interval consistent with the observed elevation of that station. Combining samples into replicate- and seasonal-aggregate datasets decreased elevation estimate uncertainty, with the greatest decrease in aggregate datasets from Fall and Winter. Information about the temporal and spatial variability of modern foraminiferal distributions was formally incorporated into the Bayesian transfer function through informative foraminifera variability priors and was applied to a Common Era relative sea-level record in New Jersey. The average difference in paleomarsh elevation estimates and uncertainties using an informative vs uninformative prior was minimal (