Modeling marine ecosystem responses to global climate change: Where are we now and where should we be going?

Modeling of global climate change is moving from global circulation model (GCM)-type projections with coupled biogeochemical models to projections of ecological responses, including food web and upper trophic levels. Marine and coastal ecosystems are highly susceptible to the impacts of global climate change and also produce significant ecosystem services. The effects of global climate change on coastal and marine ecosystems involve a much wider array of effects than the usual temperature, sea level rise, and precipitation. This paper is an overview for a collection of 12 papers that examined various aspects of global climate change on marine ecosystems and comprise this special issue. We summarized the major features of the models and analyses in the papers to determine general patterns. A wide range of ecosystems were simulated using a diverse set of modeling approaches. Models were either 3-dimensional or used a few spatial boxes, and responses to global climate change were mostly expressed as changes from a baseline condition. Three issues were identified from the across-model comparison: (a) lack of standardization of climate change scenarios, (b) the prevalence of site-specific and even unique models for upper trophic levels, and (c) emphasis on hypothesis evaluation versus forecasting. We discuss why these issues are important as global climate change assessment continues to progress up the food chain, and, when possible, offer some initial steps for going forward.

Role of zooplankton dynamics for Southern Ocean phytoplankton biomass and global biogeochemical cycles

Global ocean biogeochemistry models currently employed in climate change projections use highly simplified representations of pelagic food webs. These food webs do not necessarily include critical pathways by which ecosystems interact with ocean biogeochemistry and climate. Here we present a global biogeochemical model which incorporates ecosystem dynamics based on the representation of ten plankton functional types (PFTs); six types of phytoplankton, three types of zooplankton, and heterotrophic bacteria. We improved the representation of zooplankton dynamics in our model through (a) the explicit inclusion of large, slow-growing zooplankton, and (b) the introduction of trophic cascades among the three zooplankton types. We use the model to quantitatively assess the relative roles of iron vs. grazing in determining phytoplankton biomass in the Southern Ocean High Nutrient Low Chlorophyll (HNLC) region during summer. When model simulations do not represent crustacean macrozooplankton grazing, they systematically overestimate Southern Ocean chlorophyll biomass during the summer, even when there was no iron deposition from dust. When model simulations included the developments of the zooplankton component, the simulation of phytoplankton biomass improved and the high chlorophyll summer bias in the Southern Ocean HNLC region largely disappeared. Our model results suggest that the observed low phytoplankton biomass in the Southern Ocean during summer is primarily explained by the dynamics of the Southern Ocean zooplankton community rather than iron limitation. This result has implications for the representation of global biogeochemical cycles in models as zooplankton faecal pellets sink rapidly and partly control the carbon export to the intermediate and deep ocean.

Multiple environmental controls on phytoplankton growth strategies determine adaptive responses of the N : P ratio

The controls on the 'Redfield' N:P stoichiometry of marine phytoplankton and hence the N:P ratio of the deep ocean remain incompletely understood. Here, we use a model for phytoplankton ecophysiology and growth, based on functional traits and resource-allocation trade-offs, to show how environmental filtering, biotic interactions, and element cycling in a global ecosystem model determine phytoplankton biogeography, growth strategies and macromolecular composition. Emergent growth strategies capture major observed patterns in marine biomes. Using a new synthesis of experimental RNA and protein measurements to constrain per-ribosome translation rates, we determine a spatially variable lower limit on adaptive rRNA:protein allocation and hence on the relationship between the largest cellular P and N pools. Comparison with the lowest observed phytoplankton N:P ratios and N:P export fluxes in the Southern Ocean suggests that additional contributions from phospholipid and phosphorus storage compounds play a fundamental role in determining the marine biogeochemical cycling of these elements.

Projecting marine fish production and catch potential in Bangladesh in the 21st century under long-term environmental change and management scenarios

The fisheries sector is crucial to the Bangladeshi economy and wellbeing, accounting for 4.4% of national Gross Domestic Product (GDP) and 22.8% of agriculture sector production, and supplying ca.60% of the national animal protein intake. Fish is vital to the 16 million Bangladeshis living near the coast, a number that has doubled since the 1980s. Here we develop and apply tools to project the long term productive capacity of Bangladesh marine fisheries under climate and fisheries management scenarios, based on downscaling a global climate model, using associated river flow and nutrient loading estimates, projecting high resolution changes in physical and biochemical ocean properties, and eventually projecting fish production and catch potential under different fishing mortality targets. We place particular interest on Hilsa shad (Tenualosa ilisha), which accounts for ca.11% of total catches, and Bombay duck (Harpadon nehereus), a low price fish that is the second highest catch in Bangladesh and is highly consumed by low income communities. It is concluded that the impacts of climate change, under greenhouse emissions scenario A1B, are likely to reduce the potential fish production in the Bangladesh Exclusive Economic Zone (EEZ) by less than 10%. However, these impacts are larger for the two target species. Under sustainable management practices we expect Hilsa shad catches to show a minor decline in potential catch by 2030 but a significant (25%) decline by 2060. However, if overexploitation is allowed catches are projected to fall much further, by almost 95% by 2060, compared to the Business as Usual scenario for the start of the 21st century. For Bombay duck, potential catches by 2060 under sustainable scenarios will produce a decline of less than 20% compared to current catches. The results demonstrate that management can mitigate or exacerbate the effects of climate change on ecosystem productivity.

The assessment of a global marine ecosystem model on the basis of emergent properties and ecosystem function: a case study with ERSEM

Ecosystem models are often assessed using quantitative metrics of absolute ecosystem state, but these model-data comparisons are disproportionately vulnerable to discrepancies in the location of important circulation features. An alternative method is to demonstrate the models capacity to represent ecosystem function; the emergence of a coherent natural relationship in a simulation indicates that the model may have an appropriate representation of the ecosystem functions that lead to the emergent relationship. Furthermore, as emergent properties are large-scale properties of the system, model validation with emergent properties is possible even when there is very little or no appropriate data for the region under study, or when the hydrodynamic component of the model differs significantly from that observed in nature at the same location and time. A selection of published meta-analyses are used to establish the validity of a complex marine ecosystem model and to demonstrate the power of validation with emergent properties. These relationships include the phytoplankton community structure, the ratio of carbon to chlorophyll in phytoplankton and particulate organic matter, the ratio of particulate organic carbon to particulate organic nitrogen and the stoichiometric balance of the ecosystem. These metrics can also inform aspects of the marine ecosystem model not available from traditional quantitative and qualitative methods. For instance, these emergent properties can be used to validate the design decisions of the model, such as the range of phytoplankton functional types and their behaviour, the stoichiometric flexibility with regards to each nutrient, and the choice of fixed or variable carbon to nitrogen ratios.

Modelling the effects of climate change on the distribution and production of marine fishes: accounting for trophic interactions in a dynamic bioclimate envelope model

Climate change has already altered the distribution of marine fishes. Future predictions of fish distributions and catches based on bioclimate envelope models are available, but to date they have not considered interspecific interactions. We address this by combining the species-based Dynamic Bioclimate Envelope Model (DBEM) with a size-based trophic model. The new approach provides spatially and temporally resolved predictions of changes in species' size, abundance and catch potential that account for the effects of ecological interactions. Predicted latitudinal shifts are, on average, reduced by 20% when species interactions are incorporated, compared to DBEM predictions, with pelagic species showing the greatest reductions. Goodness-of-fit of biomass data from fish stock assessments in the North Atlantic between 1991 and 2003 is improved slightly by including species interactions. The differences between predictions from the two models may be relatively modest because, at the North Atlantic basin scale, (i) predators and competitors may respond to climate change together; (ii) existing parameterization of the DBEM might implicitly incorporate trophic interactions; and/or (iii) trophic interactions might not be the main driver of responses to climate. Future analyses using ecologically explicit models and data will improve understanding of the effects of inter-specific interactions on responses to climate change, and better inform managers about plausible ecological and fishery consequences of a changing environment.

Generalized ocean color inversion model for retrieving marine inherent optical properties

Ocean color measured from satellites provides daily, global estimates of marine inherent optical properties (IOPs). Semi-analytical algorithms (SAAs) provide one mechanism for inverting the color of the water observed by the satellite into IOPs. While numerous SAAs exist, most are similarly constructed and few are appropriately parameterized for all water masses for all seasons. To initiate community-wide discussion of these limitations, NASA organized two workshops that deconstructed SAAs to identify similarities and uniqueness and to progress toward consensus on a unified SAA. This effort resulted in the development of the generalized IOP (GIOP) model software that allows for the construction of different SAAs at runtime by selection from an assortment of model parameterizations. As such, GIOP permits isolation and evaluation of specific modeling assumptions, construction of SAAs, development of regionally tuned SAAs, and execution of ensemble inversion modeling. Working groups associated with the workshops proposed a preliminary default configuration for GIOP (GIOP-DC), with alternative model parameterizations and features defined for subsequent evaluation. In this paper, we: (1) describe the theoretical basis of GIOP; (2) present GIOP-DC and verify its comparable performance to other popular SAAs using both in situ and synthetic data sets; and, (3) quantify the sensitivities of their output to their parameterization. We use the latter to develop a hierarchical sensitivity of SAAs to various model parameterizations, to identify components of SAAs that merit focus in future research, and to provide material for discussion on algorithm uncertainties and future emsemble applications.

Technical Note: Ensuring consistent, global measurements of short-lived halocarbon gases in the ocean and atmosphere

Very short-lived halocarbons are significant sources of reactive halogen in the marine boundary layer, and likely in the upper troposphere and lower stratosphere. Quantifying ambient concentrations in the surface ocean and atmosphere is essential for understanding the atmospheric impact of these trace gas fluxes. Despite the body of literature increasing substantially over recent years, calibration issues complicate the comparison of results and limit the utility of building larger-scale databases that would enable further development of the science (e.g. sea-air flux quantification, model validation, etc.). With this in mind, thirty-one scientists from both atmospheric and oceanic halocarbon communities in eight nations gathered in London in February 2008 to discuss the scientific issues and plan an international effort toward developing common calibration scales (http://tinyurl.com/c9cg58). Here, we discuss the outputs from this meeting, suggest the compounds that should be targeted initially, identify opportunities for beginning calibration and comparison efforts, and make recommendations for ways to improve the comparability of previous and future measurements.

Global assessment of ocean carbon export by combining satellite observations and food-web models

The export of organic carbon from the surface ocean by sinking particles is an important, yet highly uncertain, component of the global carbon cycle. Here we introduce a mechanistic assessment of the global ocean carbon export using satellite observations, including determinations of net primary production and the slope of the particle size spectrum, to drive a food-web model that estimates the production of sinking zooplankton feces and algal aggregates comprising the sinking particle flux at the base of the euphotic zone. The synthesis of observations and models reveals fundamentally different and ecologically consistent regional-scale patterns in export and export efficiency not found in previous global carbon export assessments. The model reproduces regional-scale particle export field observations and predicts a climatological mean global carbon export from the euphotic zone of ~6 Pg C yr−1. Global export estimates show small variation (typically < 10%) to factor of 2 changes in model parameter values. The model is also robust to the choices of the satellite data products used and enables interannual changes to be quantified. The present synthesis of observations and models provides a path for quantifying the ocean's biological pump.

Vulnerability of coastal ecosystems to changes in harmful algal bloom distribution in response to climate change: projections based on model analysis

Harmful algal blooms (HABs), those proliferations of algae that can cause fish kills, contaminate seafood with toxins, form unsightly scums, or detrimentally alter ecosystem function have been increasing in frequency, magnitude, and duration worldwide. Here, using a global modeling approach, we show, for three regions of the globe, the potential effects of nutrient loading and climate change for two HAB genera, pelagic Prorocentrum and Karenia, each with differing physiological characteristics for growth. The projections (end of century, 2090-2100) are based on climate change resulting from the A1B scenario of the Intergovernmental Panel on Climate Change Institut Pierre Simon Laplace Climate Model (IPCC, IPSL-CM4), applied in a coupled oceanographic-biogeochemical model, combined with a suite of assumed physiological 'rules' for genera-specific bloom development. Based on these models, an expansion in area and/or number of months annually conducive to development of these HABs along the NW European Shelf-Baltic Sea system and NE Asia was projected for both HAB genera, but no expansion (Prorocentrum spp.), or actual contraction in area and months conducive for blooms (Karenia spp.), was projected in the SE Asian domain. The implications of these projections, especially for Northern Europe, are shifts in vulnerability of coastal systems to HAB events, increased regional HAB impacts to aquaculture, increased risks to human health and ecosystems, and economic consequences of these events due to losses to fisheries and ecosystem services.

Exploiting satellite earth observation to quantify current global oceanic DMS flux and its future climate sensitivity

We used coincident Envisat RA2 and AATSR temperature and wind speed data from 2008/2009 to calculate the global net sea-air flux of dimethyl sulfide (DMS), which we estimate to be 19.6 Tg S a21. Our monthly flux calculations are compared to open ocean eddy correlation measurements of DMS flux from 10 recent cruises, with a root mean square difference of 3.1 lmol m22 day21. In a sensitivity analysis, we varied temperature, salinity, surface wind speed, and aqueous DMS concentration, using fixed global changes as well as CMIP5 model output. The range of DMS flux in future climate scenarios is discussed. The CMIP5 model predicts a reduction in surface wind speed and we estimate that this will decrease the global annual sea-air flux of DMS by 22% over 25 years. Concurrent changes in temperature, salinity, and DMS concentration increase the global flux by much smaller amounts. The net effect of all CMIP5 modelled 25 year predictions was a 19% reduction in global DMS flux. 25 year DMS concentration changes had significant regional effects, some positive (Southern Ocean, North Atlantic, Northwest Pacific) and some negative (isolated regions along the Equator and in the Indian Ocean). Using satellite-detected coverage of coccolithophore blooms, our estimate of their contribution to North Atlantic DMS emissions suggests that the coccolithophores contribute only a small percentage of the North Atlantic annual flux estimate, but may be more important in the summertime and in the northeast Atlantic.