Andrew Yool

Anthropogenic ocean acidification over the twenty-first century and its impact on calcifying organisms

Todays surface ocean is saturated with respect to calcium carbonate, but increasing atmospheric carbon dioxide concentrations are reducing ocean pH and carbonate ion concentrations, and thus the level of calcium carbonate saturation. Experimental evidence suggests that if these trends continue, key marine organisms—such as corals and some plankton—will have difficulty maintaining their external calcium carbonate skeletons. Here we use 13 models of the ocean–carbon cycle to assess calcium carbonate saturation under the IS92a ‘business-as-usual’ scenario for future emissions of anthropogenic carbon dioxide. In our projections, Southern Ocean surface waters will begin to become undersaturated with respect to aragonite, a metastable form of calcium carbonate, by the year 2050. By 2100, this undersaturation could extend throughout the entire Southern Ocean and into the subarctic Pacific Ocean. When live pteropods were exposed to our predicted level of undersaturation during a two-day shipboard experiment, their aragonite shells showed notable dissolution. Our findings indicate that conditions detrimental to high-latitude ecosystems could develop within decades, not centuries as suggested previously.

The significance of nitrification for oceanic new production

The flux of organic material sinking to depth is a major control on the inventory of carbon in the ocean. To first order, the oceanic system is at equilibrium such that what goes down must come up. Because the export flux is difficult to measure directly, it is routinely estimated indirectly by quantifying the amount of phytoplankton growth, or primary production, fuelled by the upward flux of nitrate. To do so it is necessary to take into account other sources of biologically available nitrogen. However, the generation of nitrate by nitrification in surface waters has only recently received attention. Here we perform the first synthesis of open-ocean measurements of the specific rate of surface nitrification and use these to configure a global biogeochemical model to quantify the global role of nitrification. We show that for much of the world ocean a substantial fraction of the nitrate taken up is generated through recent nitrification near the surface. At the global scale, nitrification accounts for about half of the nitrate consumed by growing phytoplankton. A consequence is that many previous attempts to quantify marine carbon export, particularly those based on inappropriate use of the f-ratio (a measure of the efficiency of the ‘biological pump’), are significant overestimates.

Evaluation of ocean carbon cycle models with data-based metrics

New radiocarbon and chlorofluorocarbon-11 data from the World Ocean Circulation Experiment are used to assess a suite of 19 ocean carbon cycle models. We use the distributions and inventories of these tracers as quantitative metrics of model skill and find that only about a quarter of the suite is consistent with the new data-based metrics. This should serve as a warning bell to the larger community that not all is well with current generation of ocean carbon cycle models. At the same time, this highlights the danger in simply using the available models to represent the state-of-the-art modeling without considering the credibility of each model.

Evaluation of ocean model ventilation with CFC-11: comparison of 13 global ocean models

We compared the 13 models participating in the Ocean Carbon Model Intercomparison Project (OCMIP) with regards to their skill in matching observed distributions of CFC-11. This analysis characterizes the abilities of these models to ventilate the ocean on timescales relevant for anthropogenic CO2 uptake. We found a large range in the modeled global inventory (±30%), mainly due to differences in ventilation from the high latitudes. In the Southern Ocean, models differ particularly in the longitudinal distribution of the CFC uptake in the intermediate water, whereas the latitudinal distribution is mainly controlled by the subgrid-scale parameterization. Models with isopycnal diffusion and eddy-induced velocity parameterization produce more realistic intermediate water ventilation. Deep and bottom water ventilation also varies substantially between the models. Models coupled to a sea-ice model systematically provide more realistic AABW formation source region; however these same models also largely overestimate AABW ventilation if no specific parameterization of brine rejection during sea-ice formation is included. In the North Pacific Ocean, all models exhibit a systematic large underestimation of the CFC uptake in the thermocline of the subtropical gyre, while no systematic difference toward the observations is found in the subpolar gyre. In the North Atlantic Ocean, the CFC uptake is globally underestimated in subsurface. In the deep ocean, all but the adjoint model, failed to produce the two recently ventilated branches observed in the North Atlantic Deep Water (NADW). Furthermore, simulated transport in the Deep Western Boundary Current (DWBC) is too sluggish in all but the isopycnal model, where it is too rapid.

Global Patterns and Predictions of Seafloor Biomass Using Random Forests

A comprehensive seafloor biomass and abundance database has been constructed from 24 oceanographic institutions worldwide within the Census of Marine Life (CoML) field projects. The machine-learning algorithm, Random Forests, was employed to model and predict seafloor standing stocks from surface primary production, water-column integrated and export particulate organic matter (POM), seafloor relief, and bottom water properties. The predictive models explain 63% to 88% of stock variance among the major size groups. Individual and composite maps of predicted global seafloor biomass and abundance are generated for bacteria, meiofauna, macrofauna, and megafauna (invertebrates and fishes). Patterns of benthic standing stocks were positive functions of surface primary production and delivery of the particulate organic carbon (POC) flux to the seafloor. At a regional scale, the census maps illustrate that integrated biomass is highest at the poles, on continental margins associated with coastal upwelling and with broad zones associated with equatorial divergence. Lowest values are consistently encountered on the central abyssal plains of major ocean basins The shift of biomass dominance groups with depth is shown to be affected by the decrease in average body size rather than abundance, presumably due to decrease in quantity and quality of food supply. This biomass census and associated maps are vital components of mechanistic deep-sea food web models and global carbon cycling, and as such provide fundamental information that can be incorporated into evidence-based management.

Evaluating global ocean carbon models: The importance of realistic physics

A suite of standard ocean hydrographic and circulation metrics are applied to the equilibrium physical solutions from 13 global carbon models participating in phase 2 of the Ocean Carbon-cycle Model Intercomparison Project (OCMIP-2). Model-data comparisons are presented for sea surface temperature and salinity, seasonal mixed layer depth, meridional heat and freshwater transport, 3-D hydrographic fields, and meridional overturning. Considerable variation exists among the OCMIP-2 simulations, with some of the solutions falling noticeably outside available observational constraints. For some cases, model-model and model-data differences can be related to variations in surface forcing, subgrid-scale parameterizations, and model architecture. These errors in the physical metrics point to significant problems in the underlying model representations of ocean transport and dynamics, problems that directly affect the OCMIP predicted ocean tracer and carbon cycle variables (e.g., air-sea CO2 flux, chlorofluorocarbon and anthropogenic CO2 uptake, and export production). A substantial fraction of the large model-model ranges in OCMIP-2 biogeochemical fields (±25–40%) represents the propagation of known errors in model physics. Therefore the model-model spread likely overstates the uncertainty in our current understanding of the ocean carbon system, particularly for transport-dominated fields such as the historical uptake of anthropogenic CO2. A full error assessment, however, would need to account for additional sources of uncertainty such as more complex biological-chemical-physical interactions, biases arising from poorly resolved or neglected physical processes, and climate change.

What controls primary production in the Arctic Ocean? Results from an intercomparison of five general circulation models with biogeochemistry

As a part of Arctic Ocean Intercomparison Project, results from five coupled physical and biological ocean models were compared for the Arctic domain, defined here as north of 66.6°N. The global and regional (Arctic Ocean (AO)–only) models included in the intercomparison show similar features in terms of the distribution of present-day water column–integrated primary production and are broadly in agreement with in situ and satellite-derived data. However, the physical factors controlling this distribution differ between the models. The intercomparison between models finds substantial variation in the depth of winter mixing, one of the main mechanisms supplying inorganic nutrients over the majority of the AO. Although all models manifest similar level of light limitation owing to general agreement on the ice distribution, the amount of nutrients available for plankton utilization is different between models. Thus the participating models disagree on a fundamental question: which factor, light or nutrients, controls present-day Arctic productivity. These differences between models may not be detrimental in determining present-day AO primary production since both light and nutrient limitation are tightly coupled to the presence of sea ice. Essentially, as long as at least one of the two limiting factors is reproduced correctly, simulated total primary production will be close to that observed. However, if the retreat of Arctic sea ice continues into the future as expected, a decoupling between sea ice and nutrient limitation will occur, and the predictive capabilities of the models may potentially diminish unless more effort is spent on verifying the mechanisms of nutrient supply. Our study once again emphasizes the importance of a realistic representation of ocean physics, in particular vertical mixing, as a necessary foundation for ecosystem modeling and predictions.

Ocean fertilization: a potential means of geoengineering?

The oceans sequester carbon from the atmosphere partly as a result of biological productivity. Over much of the ocean surface, this productivity is limited by essential nutrients and we discuss whether it is likely that sequestration can be enhanced by supplying limiting nutrients. Various methods of supply have been suggested and we discuss the efficacy of each and the potential side effects that may develop as a result. Our conclusion is that these methods have the potential to enhance sequestration but that the current level of knowledge from the observations and modelling carried out to date does not provide a sound foundation on which to make clear predictions or recommendations. For ocean fertilization to become a viable option to sequester CO2, we need more extensive and targeted fieldwork and better mathematical models of ocean biogeochemical processes. Models are needed both to interpret field observations and to make reliable predictions about the side effects of large-scale fertilization. They would also be an essential tool with which to verify that sequestration has effectively taken place. There is considerable urgency to address climate change mitigation and this demands that new fieldwork plans are developed rapidly. In contrast to previous experiments, these must focus on the specific objective which is to assess the possibilities of CO2 sequestration through fertilization.

An Examination of the “Continental shelf pump” in an open ocean general circulation model

In a recent study of the shelf region of the East China Sea, Tsunogai et al. [1999] estimated that a combination of air-sea exchange and biological and physical transport processes could transfer carbon from the shelf region into the open ocean at a rate of 35 g Cm-2 yr-1. Contrasting with the solubility and biological pumps of the open ocean, they described this collective activity as the ‘‘continental shelf pump’’ and suggested that if this pump operated throughout the world’s shelf regions, it could be responsible for ocean uptake of ~1 Gt C yr-1 (~50% current ocean uptake of anthropogenic CO2). In this work a general circulation model (GCM) is used to explore the potential strength of this pump across the world’s shelves. Since the GCM does not represent the continental shelf regions explicitly, a parameterization of the pump has been used. Results of simulations find modeled pump activity very variable between shelf regions, with the East China Sea shelf behaving very similarly to the global average. Storage of pump carbon is particularly high in the Atlantic Ocean and other regions where deep water is formed. A considerable reservoir of pump carbon becomes trapped under the Arctic ice sheet. Simple extrapolations from the results suggest that should shelf regions absorb CO2 at the rate of the East China Sea, the pump would account for a net oceanic uptake of 0.6 Gt C yr-1.

Scaling up experimental ocean acidification and warming research: from individuals to the ecosystem

Understanding long-term, ecosystem-level impacts of climate change is challenging because experimental research frequently focuses on short-term, individual-level impacts in isolation. We address this shortcoming first through an interdisciplinary ensemble of novel experimental techniques to investigate the impacts of 14-month exposure to ocean acidification and warming (OAW) on the physiology, activity, predatory behaviour and susceptibility to predation of an important marine gastropod (Nucella lapillus). We simultaneously estimated the potential impacts of these global drivers on N. lapillus population dynamics and dispersal parameters. We then used these data to parameterize a dynamic bioclimatic envelope model, to investigate the consequences of OAW on the distribution of the species in the wider NE Atlantic region by 2100. The model accounts also for changes in the distribution of resources, suitable habitat and environment simulated by finely resolved biogeochemical models, under three IPCC global emissions scenarios. The experiments showed that temperature had the greatest impact on individual-level responses, while acidification had a similarly important role in the mediation of predatory behaviour and susceptibility to predators. Changes in Nucella predatory behaviour appeared to serve as a strategy to mitigate individual-level impacts of acidification, but the development of this response may be limited in the presence of predators. The model projected significant large-scale changes in the distribution of Nucella by the year 2100 that were exacerbated by rising greenhouse gas emissions. These changes were spatially heterogeneous, as the degree of impact of OAW on the combination of responses considered by the model varied depending on local-environmental conditions and resource availability. Such changes in macro-scale distributions cannot be predicted by investigating individual-level impacts in isolation, or by considering climate stressors separately. Scaling up the results of experimental climate change research requires approaches that account for long-term, multiscale responses to multiple stressors, in an ecosystem context.