Modeling the effect of hypoxia on macrobenthos production in the lower Rappahannock River, Chesapeake Bay, USA

Hypoxia in Chesapeake Bay has substantially increased in recent decades, with detrimental effects on macrobenthic production; the production of these fauna link energy transfer from primary consumers to epibenthic and demersal predators. As such, the development of accurate predictive models that determine the impact of hypoxia on macrobenthic production is important. A continuous-time, biomass-based model was developed for the lower Rappahannock River, a Bay tributary prone to seasonal hypoxia. Phytoplankton, zooplankton, and macrobenthic state variables were modeled, with a focus on quantitatively constraining the effect of hypoxia on macrobenthic biomass. This was accomplished through regression with Z': a sigmoidal function between macrobenthic biomass and dissolved oxygen concentration, derived using macrobenthic data collected from the Rappahannock River during the summers of 2007 and 2008, and applied to compute hypoxia-induced mortality as a rate process. The model was verified using independent monitoring data collected by the Chesapeake Bay Program. Simulations showed that macrobenthic biomass was strongly linked to dissolved oxygen concentrations, with fluctuations in biomass related to the duration and severity of hypoxia. Our model demonstrated that hypoxia negatively affected macrobenthic biomass, as longer durations of hypoxia and greater hypoxic severity resulted in an increasing loss in biomass. This exercise represents an important contribution to modeling anthropogenically impacted coastal ecosystems, by providing an empirically constrained relationship between hypoxia and macrobenthic biomass, and applying that empirical relationship in a mechanistic model to quantify the effect of the severity, duration, and frequency of hypoxia on benthic biomass dynamics.