The ocean has absorbed the equivalent of 39% of industrial-age fossil carbon emissions, significantly modulating the growth rate of atmospheric CO 2 and its associated impacts on climate. Despite the importance of the ocean carbon sink to climate, our understanding of the causes of its interannual-to-decadal variability remains limited. This hinders our ability to attribute its past behavior and project its future. A key period of interest is the 1990s, when the ocean carbon sink did not grow as expected. Previous explanations of this behavior have focused on variability internal to the ocean or associated with coupled atmosphere/ocean modes. Here, we use an idealized upper ocean box model to illustrate that two external forcings are sufficient to explain the pattern and magnitude of sink variability since the mid-1980s. First, the global-scale reduction in the decadal-average ocean carbon sink in the 1990s is attributable to the slowed growth rate of atmospheric pCO 2. The acceleration of atmospheric pCO 2 growth after 2001 drove recovery of the sink. Second, the global sea surface temperature response to the 1991 eruption of Mt Pinatubo explains the timing of the global sink within the 1990s. These results are consistent with previous experiments using ocean hindcast models with variable atmospheric pCO 2 and with and without climate variability. The fact that variability in the growth rate of atmospheric pCO 2 directly imprints on the ocean sink implies that there will be an immediate reduction in ocean carbon uptake as atmospheric pCO 2 responds to cuts in anthropogenic emissions. Plain Language Summary Humans have added 440 Pg of fossil fuel carbon to the atmosphere since 1750, driving up the atmospheric CO 2 concentration. But not all of this carbon remains in the atmosphere. The ocean has absorbed 39%, substantially mitigating anthropogenic climate change. Though this "ocean carbon sink" is a critical climate process, our understanding of its mechanisms remains limited. Of great interest is the unexplained slowdown of the ocean carbon sink in the 1990s and a subsequent recovery. In this work, we use a simple globally-averaged model to show that two processes external to the ocean are sufficient to explain the slowing of the ocean carbon sink in the 1990s. First, a reduced rate of accumulation of carbon in the atmosphere after 1989 reduced the atmosphere-ocean gradient that drives the ocean sink. Second, the eruption of Mt Pinatubo led to changes in ocean temperature that modified the timing of the sink from 1991 to 2001. We illustrate that the most important control on the decade-averaged magnitude of the ocean sink is variability in the growth rate of atmospheric CO 2. This implies that as future fossil fuel emission cuts drive reduced growth of atmospheric CO 2 , the ocean sink will immediately slow down.

The oceanic response to recent tropical eruptions is examined in Large Ensemble (LE) experiments from two fully coupled global climate models, the Community Earth System Model (CESM) and the Geophysical Fluid Dynamics Laboratory Earth System Model (ESM2M), each forced by a distinct volcanic forcing dataset. Following the simulated eruptions of Agung, El Chichón, and Pinatubo, the ocean loses heat and gains oxygen and carbon, in general agreement with available observations. In both models, substantial global surface cooling is accompanied by El Niño–like equatorial Pacific surface warming a year after the volcanic forcing peaks. A mechanistic analysis of the CESM and ESM2M responses to Pinatubo identifies remote wind forcing from the western Pacific as a major driver of this El Niño–like response. Following eruption, faster cooling over the Maritime Continent than adjacent oceans suppresses convection and leads to persistent westerly wind anomalies over the western tropical Pacific. These wind anomalies excite equatorial downwelling Kelvin waves and the upwelling of warm subsurface anomalies in the eastern Pacific, promoting the development of El Niño conditions through Bjerknes feedbacks a year after eruption. This El Niño–like response drives further ocean heat loss through enhanced equatorial cloud albedo, and dominates global carbon uptake as upwelling of carbon-rich waters is suppressed in the tropical Pacific. Oxygen uptake occurs primarily at high latitudes, where surface cooling intensifies the ventilation of subtropical thermocline waters. These volcanically forced ocean responses are large enough to contribute to the observed decadal variability in oceanic heat, carbon, and oxygen.

The ocean is the main source of thermal inertia in the climate system. Ocean heat uptake during recent decades has been quantified using ocean temperature measurements. However, these estimates all use the same imperfect ocean dataset and share additional uncertainty due to sparse coverage, especially before 2007. Here, we provide an independent estimate by using measurements of atmospheric oxygen (O2) and carbon dioxide (CO2) – levels of which increase as the ocean warms and releases gases – as a whole ocean thermometer. We show that the ocean gained 1.29 ± 0.79 × 1022 Joules of heat per year between 1991 and 2016, equivalent to a planetary energy imbalance of 0.80 ± 0.49 W watts per square metre of Earth’s surface. We also find that the ocean-warming effect that led to the outgassing of O2 and CO2 can be isolated from the direct effects of anthropogenic emissions and CO2 sinks. Our result – which relies on high-precision O2 atmospheric measurements dating back to 1991 – leverages an integrative Earth system approach and provides much needed independent confirmation of heat uptake estimated from ocean data.