Search

Version 1 of the Community Earth SystemModel, in the configurationwhere its full carbon cycle is enabled, is introduced and documented. In this configuration, the terrestrial biogeochemical model, which includes carbon-nitrogen dynamics and is present in earlier model versions, is coupled to an ocean biogeochemical model and atmospheric CO2 tracers. The authors provide a description of the model, detail how preindustrial-control and twentieth-century experiments were initialized and forced, and examine the behavior of the carbon cycle in those experiments. They examine how sea- and land-to-air CO2 fluxes contribute to the increase of atmospheric CO2 in the twentieth century, analyze how atmospheric CO2 and its surface fluxes vary on interannual time scales, including how they respond to ENSO, and describe the seasonal cycle of atmospheric CO2 and its surface fluxes. While the model broadly reproduces observed aspects of the carbon cycle, there are several notable biases, including having too large of an increase in atmospheric CO2 over the twentieth century and too small of a seasonal cycle of atmospheric CO2 in the Northern Hemisphere. The biases are related to a weak response of the carbon cycle to climatic variations on interannual and seasonal time scales and to twentieth-century anthropogenic forcings, including rising CO2, land-use change, and atmospheric deposition of nitrogen.

The authors compare Community Earth System Model results to marine observations for the 1990s and examine climate change impacts on biogeochemistry at the end of the twenty-first century under two future scenarios (Representative Concentration Pathways RCP4.5 and RCP8.5). Late-twentieth-century seasonally varying mixed layer depths are generally within 10 m of observations, with a Southern Ocean shallow bias. Surface nutrient and chlorophyll concentrations exhibit positive biases at low latitudes and negative biases at high latitudes. The volume of the oxygen minimum zones is overestimated.The impacts of climate change on biogeochemistry have similar spatial patterns under RCP4.5 and RCP8.5, but perturbation magnitudes are larger under RCP8.5. Increasing stratification leads to weaker nutrient entrainment and decreased primary and export production (>30% over large areas). The global-scale decreases in primary and export production scale linearly with the increases in mean sea surface temperature. There are production increases in the high nitrate, low chlorophyll (HNLC) regions, driven by lateral iron inputs from adjacent areas. The increased HNLC export partially compensates for the reductions in non-HNLC waters (~25% offset). Stabilizing greenhouse gas emissions and climate by the end of this century (as in RCP4.5) will minimize the changes to nutrient cycling and primary production in the oceans. In contrast, continued increasing emission of CO2 (as in RCP8.5) will lead to reduced productivity and significant modifications to ocean circulation and biogeochemistry by the end of this century, with more drastic changes beyond the year 2100 as the climate continues to rapidly warm.

Ocean carbon uptake and storage simulated by the Community Earth System Model, version 1–Biogeochemistry [CESM1(BGC)], is described and compared to observations. Fully coupled and ocean-ice configurations are examined; both capture many aspects of the spatial structure and seasonality of surface carbon fields. Nearly ubiquitous negative biases in surface alkalinity result from the prescribed carbonate dissolution profile. The modeled sea–air CO2 fluxes match observationally based estimates over much of the ocean; significant deviations appear in the Southern Ocean. Surface ocean pCO2 is biased high in the subantarctic and low in the sea ice zone. Formation of the water masses dominating anthropogenic CO2 (Cant) uptake in the Southern Hemisphere is weak in the model, leading to significant negative biases in Cant and chlorofluorocarbon (CFC) storage at intermediate depths. Column inventories of Cant appear too high, by contrast, in the North Atlantic. In spite of the positive bias, this marks an improvement over prior versions of the model, which underestimated North Atlantic uptake. The change in behavior is attributable to a new parameterization of density-driven overflows. CESM1(BGC) provides a relatively robust representation of the ocean–carbon cycle response to climate variability. Statistical metrics of modeled interannual variability in sea–air CO2 fluxes compare reasonably well to observationally based estimates. The carbon cycle response to key modes of climate variability is basically similar in the coupled and forced ocean-ice models; however, the two differ in regional detail and in the strength of teleconnections.

Changes in atmospheric CO2 variability during the twenty-first century may provide insight about ecosystem responses to climate change and have implications for the design of carbon monitoring programs. This paper describes changes in the three-dimensional structure of atmospheric CO2 for several representative concentration pathways (RCPs 4.5 and 8.5) using the Community Earth System Model–Biogeochemistry (CESM1-BGC). CO2 simulated for the historical period was first compared to surface, aircraft, and column observations. In a second step, the evolution of spatial and temporal gradients during the twenty-first century was examined. The mean annual cycle in atmospheric CO2 was underestimated for the historical period throughout the Northern Hemisphere, suggesting that the growing season net flux in the Community Land Model (the land component of CESM) was too weak. Consistent with weak summer drawdown in Northern Hemisphere high latitudes, simulated CO2 showed correspondingly weak north–south and vertical gradients during the summer. In the simulations of the twenty-first century, CESM predicted increases in the mean annual cycle of atmospheric CO2 and larger horizontal gradients. Not only did the mean north–south gradient increase due to fossil fuel emissions, but east–west contrasts in CO2 also strengthened because of changing patterns in fossil fuel emissions and terrestrial carbon exchange. In the RCP8.5 simulation, where CO2 increased to 1150 ppm by 2100, the CESM predicted increases in interannual variability in the Northern Hemisphere midlatitudes of up to 60% relative to present variability for time series filtered with a 2–10-yr bandpass. Such an increase in variability may impact detection of changing surface fluxes from atmospheric observations.