Modelling the responses of marine ecosystems to increasing levels of CO2

Report from Advances in Marine Ecosystem Modelling Workshop February 12-14 2007 Plymouth UK

Severe tissue damage in Atlantic cod larvae under increasing ocean acidification

Ocean acidification, caused by increasing atmospheric concentrations of CO2 (refs 1-3), is one of the most critical anthropogenic threats to marine life. Changes in seawater carbonate chemistry have the potential to disturb calcification, acid-base regulation, blood circulation and respiration, as well as the nervous system of marine organisms, leading to long-term effects such as reduced growth rates and reproduction(4,5). In teleost fishes, early life-history stages are particularly vulnerable as they lack specialized internal pH regulatory mechanisms(6,7). So far, impacts of relevant CO2 concentrations on larval fish have been found in behaviour(8,9) and otolith size(10,11), mainly in tropical, non-commercial species. Here we show detrimental effects of ocean acidification on the development of a mass-spawning fish species of high. commercial importance. We reared Atlantic cod larvae at three levels of CO2, (1) present day, (2) end of next century and (3) an extreme, coastal upwelling scenario, in a long-term (2; months) mesocosm experiment. Exposure to CO2 resulted in severe to lethal tissue damage in many internal organs, with the degree of damage increasing with CO2 concentration. As larval survival is the bottleneck to recruitment, ocean acidification has the potential to act as an additional source of natural mortality, affecting populations of already exploited fish stocks.

Climate change impacts on sea–air fluxes of CO2 in three Arctic seas: a sensitivity study using Earth observation

We applied coincident Earth observation data collected during 2008 and 2009 from multiple sensors (RA2, AATSR and MERIS, mounted on the European Space Agency satellite Envisat) to characterise environmental conditions and integrated sea-air fluxes of CO2 in three Arctic seas (Greenland, Barents, Kara). We assessed net CO2 sink sensitivity due to changes in temperature, salinity and sea ice duration arising from future climate scenarios. During the study period the Greenland and Barents seas were net sinks for atmospheric CO2, with integrated sea-air fluxes of -36 +/- 14 and -11 +/- 5 Tg C yr(-1), respectively, and the Kara Sea was a weak net CO2 source with an integrated sea-air flux of +2.2 +/- 1.4 TgC yr(-1). The combined integrated CO2 sea-air flux from all three was -45 +/- 18 TgC yr(-1). In a sensitivity analysis we varied temperature, salinity and sea ice duration. Variations in temperature and salinity led to modification of the transfer velocity, solubility and partial pressure of CO2 taking into account the resultant variations in alkalinity and dissolved organic carbon (DOC). Our results showed that warming had a strong positive effect on the annual integrated sea-air flux of CO2 (i.e. reducing the sink), freshening had a strong negative effect and reduced sea ice duration had a small but measurable positive effect. In the climate change scenario examined, the effects of warming in just over a decade of climate change up to 2020 outweighed the combined effects of freshening and reduced sea ice duration. Collectively these effects gave an integrated sea-air flux change of +4.0 TgC in the Greenland Sea, +6.0 Tg C in the Barents Sea and +1.7 Tg C in the Kara Sea, reducing the Greenland and Barents sinks by 11% and 53 %, respectively, and increasing the weak Kara Sea source by 81 %. Overall, the regional integrated flux changed by +11.7 Tg C, which is a 26% reduction in the regional sink. In terms of CO2 sink strength, we conclude that the Barents Sea is the most susceptible of the three regions to the climate changes examined. Our results imply that the region will cease to be a net CO2 sink in the 2050s.

Are jellyfish increasing in response to ocean acidification?

Ocean acidification may negatively affect calcifying plankton, opening ecological space for non-calcifying species. Recently, a study of climate-forcing of jellyfish reported the first analysis suggesting that there were more jellyfish (generally considered a noncalcifying group) when conditions were more acidic (lower pH) from one area within the North Sea. We examine this suggestion for a number of areas in the North Sea and beyond in the Northeast Atlantic using coelenterate records from the Continuous Plankton Recorder and pH data from the International Council for the Exploration of the Sea for the period 1946-2003. We could find no significant relationships between jellyfish abundance and acidic conditions in any of the regions investigated. We conclude that the role of pH in structuring zooplankton communities in the North Sea and further afield at present is tenuous.

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.

The effect of the North Atlantic Subpolar Front as a boundary in pelagic biogeography decreases with increasing depth and organism size

Broad-scale patterns in the distribution of deep-sea pelagic species and communities are poorly known. An important question is whether biogeographic boundaries identified from surface features are important in the deep mesopelagic and bathypelagic. We present community analyses of discrete-depth samples of mesozooplankton and micronekton to full-ocean depth collected in the area where the Mid-Atlantic Ridge is crossed by the Subpolar Front. The results show that the distributional discontinuity associated with the front, which is strong near the surface, decreases with increasing depth. Both the frontal separation near the surface and the community convergence at increasing depths were clearer for mesozooplankton than for micronekton.

The effect of the North Atlantic Subpolar Front as a boundary in pelagic biogeography decreases with increasing depth and organism size

Broad-scale patterns in the distribution of deep-sea pelagic species and communities are poorly known. An important question is whether biogeographic boundaries identified from surface features are important in the deep mesopelagic and bathypelagic. We present community analyses of discrete-depth samples of mesozooplankton and micronekton to full-ocean depth collected in the area where the Mid-Atlantic Ridge is crossed by the Subpolar Front. The results show that the distributional discontinuity associated with the front, which is strong near the surface, decreases with increasing depth. Both the frontal separation near the surface and the community convergence at increasing depths were clearer for mesozooplankton than for micronekton.

Making modelling count - increasing the contribution of shelf-seas community and ecosystem models to policy development and management

Marine legislation is becoming more complex and marine ecosystem-based management is specified in national and regional legislative frameworks. Shelf-seas community and ecosystem models (hereafter termed ecosystem models) are central to the delivery of ecosystem-based management, but there is limited uptake and use of model products by decision makers in Europe and the UK in comparison with other countries. In this study, the challenges to the uptake and use of ecosystem models in support of marine environmental management are assessed using the UK capability as an example. The UK has a broad capability in marine ecosystem modelling, with at least 14 different models that support management, but few examples exist of ecosystem modelling that underpin policy or management decisions. To improve understanding of policy and management issues that can be addressed using ecosystem models, a workshop was convened that brought together advisors, assessors, biologists, social scientists, economists, modellers, statisticians, policy makers, and funders. Some policy requirements were identified that can be addressed without further model development including: attribution of environmental change to underlying drivers, integration of models and observations to develop more efficient monitoring programmes, assessment of indicator performance for different management goals, and the costs and benefit of legislation. Multi-model ensembles are being developed in cases where many models exist, but model structures are very diverse making a standardised approach of combining outputs a significant challenge, and there is a need for new methodologies for describing, analysing, and visualising uncertainties. A stronger link to social and economic systems is needed to increase the range of policy-related questions that can be addressed. It is also important to improve communication between policy and modelling communities so that there is a shared understanding of the strengths and limitations of ecosystem models.

Difference equation approach to two-thermocouple sensor characterisation in constant velocity flow environments

Thermocouples are one of the most popular devices for temperature measurement due to their robustness, ease of manufacture and installation, and low cost. However, when used in certain harsh environments, for example, in combustion systems and engine exhausts, large wire diameters are required, and consequently the measurement bandwidth is reduced. This article discusses a software compensation technique to address the loss of high frequency fluctuations based on measurements from two thermocouples. In particular, a difference equation sDEd approach is proposed and compared with existing methods both in simulation and on experimental test rig data with constant flow velocity. It is found that the DE algorithm, combined with the use of generalized total least squares for parameter identification, provides better performance in terms of time constant estimation without any a priori assumption on the time constant ratios of the thermocouples.

A note on the prediction of liquid hold-up with the stratified roll wave regime for gas/liquidco-current flow in horizontal pipes.

This note presents a simple model for prediction of liquid hold-up in two-phase horizontal pipe flow for the stratified roll wave (St+RW) flow regime. Liquid hold-up data for horizontal two-phase pipe flow [1, 2, 3, 4, 5 and 6] exhibit a steady increase with liquid velocity and a more dramatic fall with increasing gas rate as shown by Hand et al. [7 and 8] for example. In addition the liquid hold-up is reported to show an additional variation with pipe diameter. Generally, if the initial liquid rate for the no-gas flow condition gives a liquid height below the pipe centre line, the flow patterns pass successively through the stratified (St), stratified ripple (St+R), stratified roll wave, film plus droplet (F+D) and finally the annular (A+D, A+RW, A+BTS) regimes as the gas rate is increased. Hand et al. [7 and 8] have given a detailed description of this progression in flow regime development and definitions of the patterns involved. Despite the fact that there are over one hundred models which have been developed to predict liquid hold-up, none have been shown to be universally useful, while only a handful have proven to be applicable to specific flow regimes [9, 10, 11 and 12]. One of the most intractable regimes to predict has been the stratified roll wave pattern where the liquid hold-up shows the most dramatic change with gas flow rate. It has been suggested that the momentum balance-type models, which give both hold-up and pressure drop prediction, can predict universally for all flow regimes but particularly in the case of the difficult stratified roll wave pattern. Donnelly [1] recently demonstrated that the momentum balance models experienced some difficulties in the prediction of this regime. Without going into lengthy details, these models differ in the assumed friction factor or shear stress on the surfaces within the pipe particularly at the liquid–gas interface. The Baker–Jardine model [13] when tested against the 0.0454 m i.d. data of Nguyen [2] exhibited a wide scatter for both liquid hold-up and pressure drop as shown in Fig. 1. The Andritsos–Hanratty model [14] gave better prediction of pressure drop but a wide scatter for liquid hold-up estimation (cf. Fig. 2) when tested against the 0.0935 m i.d. data of Hand [5]. The Spedding–Hand model [15], shown in Fig. 3 against the data of Hand [5], gave improved performance but was still unsatisfactory with the prediction of hold-up for stratified-type flows. The MARS model of Grolman [6] gave better prediction of hold-up (cf. Fig. 4) but deterioration in the estimation of pressure drop when tested against the data of Nguyen [2]. Thus no method is available that will accurately predict liquid hold-up across the whole range of flow patterns but particularly for the stratified plus roll wavy regime. The position is particularly unfortunate since the stratified-type regimes are perhaps the most predominant pattern found in multiphase lines.