Keith Lindsay

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.

Climate-carbon cycle feedback analysis: Results from the C4MIP model intercomparison

Eleven coupled climate–carbon cycle models used a common protocol to study the coupling between climate change and the carbon cycle. The models were forced by historical emissions and the Intergovernmental Panel on Climate Change (IPCC) Special Report on Emissions Scenarios (SRES) A2 anthropogenic emissions of CO2 for the 1850–2100 time period. For each model, two simulations were performed in order to isolate the impact of climate change on the land and ocean carbon cycle, and therefore the climate feedback on the atmospheric CO2 concentration growth rate. There was unanimous agreement among the models that future climate change will reduce the efficiency of the earth system to absorb the anthropogenic carbon perturbation. A larger fraction of anthropogenic CO2 will stay airborne if climate change is accounted for. By the end of the twenty-first century, this additional CO2 varied between 20 and 200 ppm for the two extreme models, the majority of the models lying between 50 and 100 ppm. The higher CO2 levels led to an additional climate warming ranging between 0.1° and 1.5°C. All models simulated a negative sensitivity for both the land and the ocean carbon cycle to future climate. However, there was still a large uncertainty on the magnitude of these sensitivities. Eight models attributed most of the changes to the land, while three attributed it to the ocean. Also, a majority of the models located the reduction of land carbon uptake in the Tropics. However, the attribution of the land sensitivity to changes in net primary productivity versus changes in respiration is still subject to debate; no consensus emerged among the models.

The Community Earth System Model: A Framework for Collaborative Research

The Community Earth System Model (CESM) is a flexible and extensible community tool used to investigate a diverse set of Earth system interactions across multiple time and space scales. This global coupled model significantly extends its predecessor, the Community Climate System Model, by incorporating new Earth system simulation capabilities. These comprise the ability to simulate biogeochemical cycles, including those of carbon and nitrogen, a variety of atmospheric chemistry options, the Greenland Ice Sheet, and an atmosphere that extends to the lower thermosphere. These and other new model capabilities are enabling investigations into a wide range of pressing scientific questions, providing new foresight into possible future climates and increasing our collective knowledge about the behavior and interactions of the Earth system. Simulations with numerous configurations of the CESM have been provided to phase 5 of the Coupled Model Intercomparison Project (CMIP5) and are being analyzed by the broad com...

Upper ocean ecosystem dynamics and iron cycling in a global three-dimensional model

A global 3-D marine ecosystem model with several key phytoplankton functional groups, multiple limiting nutrients, explicit iron cycling, and a mineral ballast/organic matter parameterization is run within a global ocean circulation model. The coupled biogeochemistry/ecosystem/circulation (BEC) model reproduces known basin-scale patterns of primary production, biogenic silica production, calcification, chlorophyll, macronutrient and dissolved iron concentrations. The model captures observed High Nitrate, Low Chlorophyll (HNLC) conditions in the Southern Ocean, subarctic and equatorial Pacific. Spatial distributions of nitrogen fixation are in general agreement with field data, with total N-fixation of 54 Tg N. Diazotrophs directly account for a small fraction of primary production (0.54%) but indirectly support 10% of primary production and nearly 8% of particulate organic carbon (POC) export. Diatoms disproportionately contribute to export of POC out of surface waters, but CaCO3 from the coccolithophores is the key driver of POC flux to the deep ocean in the model. An iron source from shallow ocean sediments is found critical in preventing iron limitation in shelf regions, most notably in the Arctic Ocean, but has a relatively localized impact. In contrast, global-scale primary production, export production, and nitrogen fixation are all sensitive to variations in atmospheric mineral dust inputs. The residence time for dissolved iron in surface waters is estimated to be a few years to a decade. Most of the iron utilized by phytoplankton is from subsurface sources supplied by mixing, entrainment, and ocean circulation. However, due to the short residence time of iron in the upper ocean, this subsurface iron pool is critically dependent on continual replenishment from atmospheric dust deposition and, to a lesser extent, lateral transport from shelf regions.A global three-dimensional marine ecosystem model with several key phytoplankton functional groups, multiple limiting nutrients, explicit iron cycling, and a mineral ballast/organic matter parameterization is run within a global ocean circulation model. The coupled biogeochemistry/ecosystem/circulation (BEC) model reproduces known basin-scale patterns of primary and export production, biogenic silica production, calcification, chlorophyll, macronutrient and dissolved iron concentrations. The model captures observed high nitrate, low chlorophyll (HNLC) conditions in the Southern Ocean, subarctic and equatorial Pacific. Spatial distributions of nitrogen fixation are in general agreement with field data, with total N-fixation of 55 Tg N. Diazotrophs directly account for a small fraction of primary production (0.5%) but indirectly support 10% of primary production and 8% of sinking particulate organic carbon (POC) export. Diatoms disproportionately contribute to export of POC out of surface waters, but CaCO3 from the coccolithophores is the key driver of POC flux to the deep ocean in the model. An iron source from shallow ocean sediments is found critical in preventing iron limitation in shelf regions, most notably in the Arctic Ocean, but has a relatively localized impact. In contrast, global-scale primary production, export production, and nitrogen fixation are all sensitive to variations in atmospheric mineral dust inputs. The residence time for dissolved iron in the upper ocean is estimated to be a few years to a decade. Most of the iron utilized by phytoplankton is from subsurface sources supplied by mixing, entrainment, and ocean circulation. However, owing to the short residence time of iron in the upper ocean, this subsurface iron pool is critically dependent on continual replenishment from atmospheric dust deposition and, to a lesser extent, lateral transport from shelf regions.

The Community Earth System Model (CESM) large ensemble project: A community resource for studying climate change in the presence of internal climate variability

AbstractWhile internal climate variability is known to affect climate projections, its influence is often underappreciated and confused with model error. Why? In general, modeling centers contribute a small number of realizations to international climate model assessments [e.g., phase 5 of the Coupled Model Intercomparison Project (CMIP5)]. As a result, model error and internal climate variability are difficult, and at times impossible, to disentangle. In response, the Community Earth System Model (CESM) community designed the CESM Large Ensemble (CESM-LE) with the explicit goal of enabling assessment of climate change in the presence of internal climate variability. All CESM-LE simulations use a single CMIP5 model (CESM with the Community Atmosphere Model, version 5). The core simulations replay the twenty to twenty-first century (1920–2100) 30 times under historical and representative concentration pathway 8.5 external forcing with small initial condition differences. Two companion 1000+-yr-long preindu...

Oceanic sources, sinks, and transport of atmospheric CO2

We synthesize estimates of the contemporary net air-sea CO2 flux on the basis of an inversion of interior ocean carbon observations using a suite of 10 ocean general circulation models (Mikaloff Fletcher et al., 2006, 2007) and compare them to estimates based on a new climatology of the air-sea difference of the partial pressure of CO2 (pCO2) (Takahashi et al., 2008). These two independent flux estimates reveal a consistent description of the regional distribution of annual mean sources and sinks of atmospheric CO2 for the decade of the 1990s and the early 2000s with differences at the regional level of generally less than 0.1 Pg C a−1. This distribution is characterized by outgassing in the tropics, uptake in midlatitudes, and comparatively small fluxes in thehigh latitudes. Both estimates point toward a small (∼ −0.3 Pg C a−1) contemporary CO2 sink in the Southern Ocean (south of 44°S), a result of the near cancellation between a substantial outgassing of natural CO2 and a strong uptake of anthropogenic CO2. A notable exception in the generally good agreement between the two estimates exists within the Southern Ocean: the ocean inversion suggests a relatively uniform uptake, while the pCO2-based estimate suggests strong uptake in the region between 58°S and 44°S, and a source in the region south of 58°S. Globally and for a nominal period between 1995 and 2000, the contemporary net air-sea flux of CO2 is estimated to be −1.7 ± 0.4 Pg C a−1 (inversion) and −1.4 ± 0.7 Pg C a−1 (pCO2-climatology), respectively, consisting of an outgassing flux of river-derived carbon of ∼+0.5 Pg C a−1, and an uptake flux of anthropogenic carbon of −2.2 ± 0.3 Pg C a−1 (inversion) and −1.9 ± 0.7 Pg C a−1 (pCO2-climatology). The two flux estimates also imply a consistent description of the contemporary meridional transport of carbon with southward ocean transport throughout most of the Atlantic basin, and strong equatorward convergence in the Indo-Pacific basins. Both transport estimates suggest a small hemispheric asymmetry with a southward transport of between −0.2 and −0.3 Pg C a−1 across the equator. While the convergence of these two independent estimates is encouraging and suggests that it is now possible to provide relatively tight constraints for the net air-sea CO2 fluxes at the regional basis, both studies are limited by their lack of consideration of long-term changes in the ocean carbon cycle, such as the recent possible stalling in the expected growth of the Southern Ocean carbon sink.

Carbon-concentration and carbon-climate feedbacks in CMIP5 Earth System Models

AbstractThe magnitude and evolution of parameters that characterize feedbacks in the coupled carbon–climate system are compared across nine Earth system models (ESMs). The analysis is based on results from biogeochemically, radiatively, and fully coupled simulations in which CO2 increases at a rate of 1% yr−1. These simulations are part of phase 5 of the Coupled Model Intercomparison Project (CMIP5). The CO2 fluxes between the atmosphere and underlying land and ocean respond to changes in atmospheric CO2 concentration and to changes in temperature and other climate variables. The carbon–concentration and carbon–climate feedback parameters characterize the response of the CO2 flux between the atmosphere and the underlying surface to these changes. Feedback parameters are calculated using two different approaches. The two approaches are equivalent and either may be used to calculate the contribution of the feedback terms to diagnosed cumulative emissions. The contribution of carbon–concentration feedback to...

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.

Inverse estimates of anthropogenic CO2 uptake, transport, and storage by the ocean

deviation of the models weighted by a CFC-based model skill score, which reduces the error range and emphasizes those models that have been shown to reproduce observed tracer concentrations most accurately. The greatest anthropogenic CO2 uptake occurs in the Southern Ocean and in the tropics. The flux estimates imply vigorous northward transport in the Southern Hemisphere, northward cross-equatorial transport, and equatorward transport at high northern latitudes. Compared with forward simulations, we find substantially more uptake in the Southern Ocean, less uptake in the Pacific Ocean, and less global uptake. The large-scale spatial pattern of the estimated flux is generally insensitive to possible biases in the data and the models employed. However, the global uptake scales approximately linearly with changes in the global anthropogenic CO2 inventory. Considerable uncertainties remain in some regions, particularly the Southern Ocean.

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.