Phytoplankton assemblage characteristics in recurrently fluctuating environments

Annual variations in biogeochemical and physical processes can lead to nutrient variability and seasonal patterns in phytoplankton productivity and assemblage structure. In many coastal systems river inflow and water exchange with the ocean varies seasonally, and alternating periods can arise where the nutrient most limiting to phytoplankton growth switches. Transitions between these alternating periods can be sudden or gradual and this depends on human activities, such as reservoir construction and interbasin water transfers. How such activities might influence phytoplankton assemblages is largely unknown. Here, we employed a multispecies, multi-nutrient model to explore how nutrient loading switching mode might affect characteristics of phytoplankton assemblages. The model is based on the Monod-relationship, predicting an instantaneous growth rate from ambient inorganic nutrient concentrations whereas the limiting nutrient at any given time was determined by Liebig’s Law of the Minimum. Our simulated phytoplankton assemblages self-organized from species rich pools over a 15-year period, and only the surviving species were considered as assemblage members. Using the model, we explored the interactive effects of complementarity level in trait trade-offs within phytoplankton assemblages and the amount of noise in the resource supply concentrations. We found that the effect of shift from a sudden resource supply transition to a gradual one, as observed in systems impacted by watershed development, was dependent on the level of complementarity. In the extremes, phytoplankton species richness and relative overyielding increased when complementarity was lowest, and phytoplankton biomass increased greatly when complementarity was highest. For low-complementarity simulations, the persistence of poorer-performing phytoplankton species of intermediate R*s led to higher richness and relative overyielding. For high-complementarity simulations, the formation of phytoplankton species clusters and niche compression enabled higher biomass accumulation. Our findings suggest that an understanding of factors influencing the emergence of life history traits important to complementarity is necessary to predict the impact of watershed development on phytoplankton productivity and assemblage structure

An NPZ Model with State-Dependent Delay due to Size-Structure in Juvenile Zooplankton

The study of planktonic ecosystems is important as they make up the bottom trophic levels of aquatic food webs. We study a closed Nutrient-Phytoplankton-Zooplankton (NPZ) model that includes size structure in the juvenile zooplankton. The closed nature of the system allows the formulation of a conservation law of biomass that governs the system. The model consists of a system of nonlinear ordinary differential equation coupled to a partial differential equation. We are able to transform this system into a system of delay differential equations where the delay is of threshold type and is state-dependent. The system of delay differential equations can be further transformed into one with fixed delay. Using the different forms of the model we perform a qualitative analysis of the solutions, which includes studying existence and uniqueness, positivity and boundedness, local and global stability, and conditions for extinction. Key parameters that are explored are the total biomass in the system and the maturity level at which the juvenile zooplankton reach maturity. Numerical simulations are also performed to verify our analytical results

Modeling dimethylsulphide production in the upper ocean

Dimethylsulphide (DMS) is produced by upper ocean ecosystems and emitted to the atmosphere, where it may have an important role in climate regulation. Several attempts to quantify the role of DMS in climate change have been undertaken in modeling studies. We examine a model of biogenic DMS production and describe its endogenous dynamics and sensitivities. We extend the model to develop a one-dimensional version that more accurately resolves the important processes of the mixed layer in determining the ecosystem dynamics. Comparisons of the results of the one-dimensional model with an empirical relationship that describes the global distribution of DMS, and also with vertical profiles of DMS in the upper ocean measured at the Bermuda Atlantic Time Series, suggest that the model represents the interaction between the biological and physical processes well on local and global scales. Our analysis of the model confirms its veracity and provides insights into the important processes determining DMS concentration in the oceans

Estuary-enhanced upwelling of marine nutrients fuels coastal productivity in the U.S. Pacific Northwest

© 2014. American Geophysical Union. The Pacific Northwest (PNW) shelf is the most biologically productive region in the California Current System. A coupled physical-biogeochemical model is used to investigate the influence of freshwater inputs on the productivity of PNW shelf waters using realistic hindcasts and model experiments that omit outflow from the Columbia River and Strait of Juan de Fuca (outlet for the Salish Sea estuary). Outflow from the Strait represents a critical source of nitrogen to the PNW shelf-accounting for almost half of the primary productivity on the Vancouver Island shelf, a third of productivity on the Washington shelf, and a fifth of productivity on the Oregon shelf during the upwelling season. The Columbia River has regional effects on the redistribution of phytoplankton, but does not affect PNW productivity as strongly as does the Salish Sea. A regional nutrient budget shows that nitrogen exiting the Strait is almost entirely (98%) of ocean-origin-upwelled into the Strait at depth, mixed into surface waters by tidal mixing, and returned to the coastal ocean. From the standpoint of nitrogen availability in the coastal euphotic zone, the estuarine circulation driven by freshwater inputs to the Salish Sea is more important than the supply of terrigenous nitrogen by rivers. Nitrogen-rich surface waters exiting the Strait follow two primary pathways-to the northwest in the Vancouver Island Coastal Current and southward toward the Washington and Oregon shelves. Nitrogen flux from the Juan de Fuca Strait and Eddy Region to these shelves is comparable to flux from local wind-driven upwelling

Modular System for Shelves and Coasts (MOSSCO v1.0) - a flexible and multi-component framework for coupled coastal ocean ecosystem modelling

Shelf and coastal sea processes extend from the atmosphere through the water column and into the sea bed. These processes are driven by physical, chemical, and biological interactions at local scales, and they are influenced by transport and cross strong spatial gradients. The linkages between domains and many different processes are not adequately described in current model systems. Their limited integration level in part reflects lacking modularity and flexibility; this shortcoming hinders the exchange of data and model components and has historically imposed supremacy of specific physical driver models. We here present the Modular System for Shelves and Coasts (MOSSCO, http://www.mossco.de), a novel domain and process coupling system tailored---but not limited--- to the coupling challenges of and applications in the coastal ocean. MOSSCO builds on the existing coupling technology Earth System Modeling Framework and on the Framework for Aquatic Biogeochemical Models, thereby creating a unique level of modularity in both domain and process coupling; the new framework adds rich metadata, flexible scheduling, configurations that allow several tens of models to be coupled, and tested setups for coastal coupled applications. That way, MOSSCO addresses the technology needs of a growing marine coastal Earth System community that encompasses very different disciplines, numerical tools, and research questions.Comment: 30 pages, 6 figures, submitted to Geoscientific Model Development Discussion

Plankton functional group models – An assessment

This Discussant’s Report provides a summary of the discussions that followed presentation of the approaches and ideas described in Thingstad et al. (this volume). The discussions, which addressed aspects of conceptual understanding and parameterization that are relevant to development of ecosystem models capable of emergent behavior at a range of scales, the benefits of functional group modeling, and some of the limitations of this approach, provide insights that are relevant to setting directions for future research efforts. One important point emerging from the discussions was that reconciling the requirements of simplicity versus complexity with the desire to obtain predictive capability is an important area where biogeochemical and ecosystem models can be improved