Nonlinear Interactions between Climate and Atmospheric Carbon Dioxide Drivers of Terrestrial and Marine Carbon Cycle Changes
Quantifying interactions between global biogeochemical cycles and the Earth system is important for predicting future atmospheric composition and informing energy policy. We applied a feedback analysis framework to three sets of Historical (1850–2005), Representative Concentration Pathway 8.5 (2006–2100), and its extension (2101–2300) simulations from the Community Earth System Model version 1.0 (CESM1(BGC)) to quantify drivers of terrestrial and ocean responses of carbon uptake. In the biogeochemically coupled simulation (BGC), the effects of CO2fertilization and nitrogen deposition influenced marine and terrestrial carbon cycling. In the radiatively coupled simulation (RAD), the effects of rising temperature and circulation changes due to radiative forcing from CO2, other greenhouse gases, and aerosols were the sole drivers of carbon cycle changes. In the third, fully coupled simulation (FC), both the biogeochemical and radiative coupling effects acted simultaneously. We found that climate–carbon sensitivities derived from RAD simulations produced a net ocean carbon storage climate sensitivity that was weaker and a net land carbon storage climate sensitivity that was stronger than those diagnosed from the FC and BGC simulations. For the ocean, this nonlinearity was associated with warming-induced weakening of ocean circulation and mixing that limited exchange of dissolved inorganic carbon between surface and deeper water masses. For the land, this nonlinearity was associated with strong gains in gross primary production in the FC simulation, driven by enhancements in the hydrological cycle and increased nutrient availability. We developed and applied a nonlinearity metric to rank model responses and driver variables. The climate–carbon cycle feedback gain at 2300 was 42% higher when estimated from climate–carbon sensitivities derived from the difference between FC and BGC than when derived from RAD. We re-analyzed other CMIP5 model results to quantify the effects of such nonlinearities on their projected climate–carbon cycle feedback gains.