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.

The global distribution of phytoplankton size spectrum and size classes from their light-absorption spectra derived from satellite data

The absorption spectra of phytoplankton in the visible domain hold implicit information on the phytoplankton community structure. Here we use this information to retrieve quantitative information on phytoplankton size structure by developing a novel method to compute the exponent of an assumed power-law for their particle-size spectrum. This quantity, in combination with total chlorophyll-a concentration, can be used to estimate the fractional concentration of chlorophyll in any arbitrarily-defined size class of phytoplankton. We further define and derive expressions for two distinct measures of cell size of mixed. populations, namely, the average spherical diameter of a bio-optically equivalent homogeneous population of cells of equal size, and the average equivalent spherical diameter of a population of cells that follow a power-law particle-size distribution. The method relies on measurements of two quantities of a phytoplankton sample: the concentration of chlorophyll-a, which is an operational index of phytoplankton biomass, and the total absorption coefficient of phytoplankton in the red peak of visible spectrum at 676 nm. A sensitivity analysis confirms that the relative errors in the estimates of the exponent of particle size spectra are reasonably low. The exponents of phytoplankton size spectra, estimated for a large set of in situ data from a variety of oceanic environments (similar to 2400 samples), are within a reasonable range; and the estimated fractions of chlorophyll in pico-, nano- and micro-phytoplankton are generally consistent with those obtained by an independent, indirect method based on diagnostic pigments determined using high-performance liquid chromatography. The estimates of cell size for in situ samples dominated by different phytoplankton types (diatoms, prymnesiophytes, Prochlorococcus, other cyanobacteria and green algae) yield nominal sizes consistent with the taxonomic classification. To estimate the same quantities from satellite-derived ocean-colour data, we combine our method with algorithms for obtaining inherent optical properties from remote sensing. The spatial distribution of the size-spectrum exponent and the chlorophyll fractions of pico-, nano- and micro-phytoplankton estimated from satellite remote sensing are in agreement with the current understanding of the biogeography of phytoplankton functional types in the global oceans. This study contributes to our understanding of the distribution and time evolution of phytoplankton size structure in the global oceans.

Global sea-to-air flux climatology for bromoform, dibromomethane and methyl iodide

Volatile halogenated organic compounds containing bromine and iodine, which are naturally produced in the ocean, are involved in ozone depletion in both the troposphere and stratosphere. Three prominent compounds transporting large amounts of marine halogens into the atmosphere are bromoform (CHBr3), dibromomethane (CH2Br2) and methyl iodide (CH3I). The input of marine halogens to the stratosphere has been estimated from observations and modelling studies using low-resolution oceanic emission scenarios derived from top-down approaches. In order to improve emission inventory estimates, we calculate data-based high resolution global sea-to-air flux estimates of these compounds from surface observations within the HalOcAt (Halocarbons in the Ocean and Atmosphere) database (https://halocat.geomar.de/). Global maps of marine and atmospheric surface concentrations are derived from the data which are divided into coastal, shelf and open ocean regions. Considering physical and biogeochemical characteristics of ocean and atmosphere, the open ocean water and atmosphere data are classified into 21 regions. The available data are interpolated onto a 1 degrees x 1 degrees grid while missing grid values are interpolated with latitudinal and longitudinal dependent regression techniques reflecting the compounds' distributions. With the generated surface concentration climatologies for the ocean and atmosphere, global sea-to-air concentration gradients and sea-to-air fluxes are calculated. Based on these calculations we estimate a total global flux of 1.5/2.5 Gmol Br yr(-1) for CHBr3, 0.78/0.98 Gmol Br yr(-1) for CH2Br2 and 1.24/1.45 Gmol Br yr(-1) for CH3I (robust fit/ordinary least squares regression techniques). Contrary to recent studies, negative fluxes occur in each sea-to-air flux climatology, mainly in the Arctic and Antarctic regions. "Hot spots" for global polybromomethane emissions are located in the equatorial region, whereas methyl iodide emissions are enhanced in the subtropical gyre regions. Inter-annual and seasonal variation is contained within our flux calculations for all three compounds. Compared to earlier studies, our global fluxes are at the lower end of estimates, especially for bromoform. An under-representation of coastal emissions and of extreme events in our estimate might explain the mismatch between our bottom-up emission estimate and top-down approaches.

The Continuous Plankton Recorder: towards a global perspective

The Continuous Plankton Recorder (CPR) began its first routine deployment to collect plankton samples in the North Sea in 1931. The better part of a century after that event, the Recorder is used for sampling plankton widely in many oceans. The CPR is designed to be towed behind ships of opportunity (such as commercial and passenger vessels) and to collect plankton samples along the ship’s path. The samples provide information on the broad scale spatial distributions of larger hard-shelled phytoplankton and robust mesozooplankton organisms (Richardson et al., 2006). Changes in these distributions provide indices for the biological health of the ocean including those required for fisheries and for assessment of climate change impacts.

Unanticipated biological changes and global warming

Evidence of global warming is now unequivocal, and studies suggest that it has started to influence natural systems of the planet, including the oceans. However, in the marine environment, it is well-known that species and ecosystems can also be influenced by natural sources of large-scale hydro-climatological variability. The North Atlantic Oscillation (NAO) was negatively correlated with the mean abundance of one of the subarctic key species Calanus finmarchicus in the North Sea. This correlation was thought to have broken down in 1996, however, the timing has never been tested statistically. The present study revisits this unanticipated change and reveals that the correlation did not break down in 1996 as originally proposed but earlier, at the time of an abrupt ecosystem shift in the North Sea in the 1980s. Furthermore, the analyses demonstrate that the correlation between the NAO and C. finmarchicus abundance is modulated by the thermal regime of the North Sea, which in turn covaries positively with global temperature anomalies. This study thereby provides evidence that global climate change is likely to alter some empirical relationships found in the past between species abundance or the ecosystem state and large-scale natural sources of hydro-climatological variability. A theory is proposed to explain how this might happen. These unanticipated changes, also called ‘surprises’ in climatic research, are a direct consequence of the complexity of both climatic and biological systems. In this period of rapid climate change, it is therefore hazardous to integrate meteo-oceanic indices such as the NAO in models used in the management of living resources, as it has been sometimes attempted in the past.

Global synchrony of an accelerating risein sea surface temperature

The oceans have shown a recent rapid and accelerating rise in temperature with, given the close link between temperature and marine organisms, pronounced effects on ecosystems. Here we describe for the first time a globally synchronous pattern of pulsed short period (�1 year long) emanations of warm sea surface temperature anomalies from tropical seas towards the poles on the shelf/slope with an intensification of the warming after the 1976/1977, 1986/1987 and 1997/1998 El Nin˜os. On the eastern margins of continents the anomalies propagate towards the poles in part by largely baroclinic boundary currents, reinforced by regional atmospheric warming. The processes contributing to the less continuous warm anomalies on western margins are linked to the transfer of warmth from adjacent western boundary currents. These climate induced events show a close parallelism with the timing of ecosystem changes in shelf seas, important for fisheries and ecosystem services, and melting of sea-ice.

Global latitudinal variations in marine copepod diversity and environmental factors

Latitudinal gradients in diversity are among the most striking features in ecology. For terrestrial species, climate (i.e. temperature and precipitation) is believed to exert a strong influence on the geographical distributions of diversity through its effects on energy availability. Here, we provide the first global description of geographical variation in the diversity of marine copepods, a key trophic link between phytoplankton and fish, in relation to environmental variables. We found a polar-tropical difference in copepod diversity in the Northern Hemisphere where diversity peaked at subtropical latitudes. In the Southern Hemisphere, diversity showed a tropical plateau into the temperate regions. This asymmetry around the Equator may be explained by climatic conditions, in particular the influence of the Inter-Tropical Convergence Zone, prevailing mainly in the northern tropical region. Ocean temperature was the most important explanatory factor among all environmental variables tested, accounting for 54 per cent of the variation in diversity. Given the strong positive correlation between diversity and temperature, local copepod diversity, especially in extra-tropical regions, is likely to increase with climate change as their large-scale distributions respond to climate warming.

Going Global: planning the next 80 years of the Continuous Plankton Recorder Survey

Going Global: planning the next 80 years of the Continuous Plankton Recorder Survey. Operated by the Sir Alister Hardy Foundation for Ocean Science (SAHFOS), the Continuous Plankton Recorder (CPR) survey is the world’s largest, sampling 4 ocean basins, and longest running (since 1931) plankton biodiversity monitoring programme. Having sampled enough miles to circumnavigate the globe over 200 times, the CPR database houses over 2.5 million entries, describing the distribution of 500 phytoplankton and zooplankton taxa. Routinely sampling in the Arctic, Atlantic, Pacific and Southern Oceans, the survey analyses 4000 samples yearly. Data collected from these samples are made freely available for bona fide scientific purposes. The CPR survey data is used to generate a better understanding of changes in the plankton and to date some 1000 papers have been published on plankton biodiversity. This year sees the 80th anniversary of the CPR survey and to celebrate and build upon this unique monitoring programme, SAHFOS intends to further develop its global plankton perspective. Work will be extended into the South Atlantic and Indian Ocean and an international partnership with complementary surveys in Australia, Canada, America, Japan and South Africa will be implemented. The Digital Object will describe the CPR survey using compilations made by Plymouth Art College and BBC film footage.

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.

An updated climatology of surface dimethlysulfide concentrations and emission fluxes in the global ocean

The potentially significant role of the biogenic trace gas dimethylsulfide (DMS) in determining the Earth's radiation budget makes it necessary to accurately reproduce seawater DMS distribution and quantify its global flux across the sea/air interface. Following a threefold increase of data (from 15,000 to over 47,000) in the global surface ocean DMS database over the last decade, new global monthly climatologies of surface ocean DMS concentration and sea-to-air emission flux are presented as updates of those constructed 10 years ago. Interpolation/extrapolation techniques were applied to project the discrete concentration data onto a first guess field based on Longhurst's biogeographic provinces. Further objective analysis allowed us to obtain the final monthly maps. The new climatology projects DMS concentrations typically in the range of 1–7 nM, with higher levels occurring in the high latitudes, and with a general trend toward increasing concentration in summer. The increased size and distribution of the observations in the DMS database have produced in the new climatology substantially lower DMS concentrations in the polar latitudes and generally higher DMS concentrations in regions that were severely undersampled 10 years ago, such as the southern Indian Ocean. Using the new DMS concentration climatology in conjunction with state-of-the-art parameterizations for the sea/air gas transfer velocity and climatological wind fields, we estimate that 28.1 (17.6–34.4) Tg of sulfur are transferred from the oceans into the atmosphere annually in the form of DMS. This represents a global emission increase of 17% with respect to the equivalent calculation using the previous climatology. This new DMS climatology represents a valuable tool for atmospheric chemistry, climate, and Earth System models.