Pinatubo erupted during the first decadal survey of ocean biogeochemistry, embedding its climate fingerprint into foundational ocean biogeochemical observations and complicating the interpretation of long-term biogeochemical change. Here, we quantify the influence of the Pinatubo climate perturbation on externally forced decadal and multi-decadal changes in key ocean biogeochemical quantities using a large ensemble simulation of the Community Earth System Model designed to isolate the effects of Pinatubo, which cleanly captures the ocean biogeochemical response to the eruption. We find increased uptake of apparent oxygen utilization and preindustrial carbon over 1993-2003. Nearly 100\% of the forced response in these quantities are attributable to Pinatubo. The eruption caused enhanced ventilation of the North Atlantic, as evidenced by deep ocean chlorofluorocarbon changes that appear 10-15 years after the eruption. Our results help contextualize observed change and contribute to improved constraints on uncertainty in the global carbon budget and ocean deoxygenation.

We use a statistical emulation technique to construct synthetic ensembles of global and regional sea-air carbon dioxide (CO2) flux from four observation-based products over 1985-2014. Much like ensembles of Earth system models that are constructed by perturbing their initial conditions, our synthetic ensemble members exhibit different phasing of internal variability and a common externally forced signal. Our synthetic ensembles illustrate an important role for internal variability in the temporal evolution of global and regional CO2 flux and produce a wide range of possible trends over 1990-1999 and 2000-2009. We assume a specific externally forced signal and calculate the likelihood of the observed trend given the distribution of synthetic trends during these two periods. Over the decade 1990-1999, three of the four observation-based products exhibit small negative trends in globally integrated sea-air CO2 flux (i.e., enhanced ocean CO2 absorption with time) that are highly probable (44-72% chance of occurrence) in their respective synthetic trend distributions. Over the decade 2000-2009, however, three of the four products show large negative trends in globally integrated sea-air CO2 flux that are somewhat improbable (17-19% chance of occurrence). Our synthetic ensembles suggest that the largest observation-based positive trends in global and Southern Ocean CO2 flux over 1990-1999 and the largest negative trends over 2000-2009 are somewhat improbable (<30% chance of occurrence). Our approach provides a new understanding of the role of internal and external processes in driving sea-air CO2 flux variability.

Large volcanic eruptions drive significant climate perturbations through major anomalies in radiative fluxes and the resulting widespread cooling of the surface and upper ocean. Recent studies suggest that these eruptions also drive important variability in air-sea carbon and oxygen fluxes. By simulating the Earth system using two initial-condition large ensembles, with and without the aerosol forcing associated with the Mt. Pinatubo eruption in June 1991, we isolate the impact of this event on ocean physical and biogeochemical properties. The Mt. Pinatubo eruption generated significant anomalies in surface fluxes and the ocean interior inventories of heat, oxygen, and carbon. Pinatubo-driven changes persist for multiple years in the upper ocean and permanently modify the ocean’s heat, oxygen, and carbon inventories. Positive anomalies in oxygen concentrations emerge immediately post-eruption and penetrate into the deep ocean. In contrast, carbon anomalies intensify in the upper ocean over several years post-eruption, and are largely confined to the upper 150 m. In the tropics and northern high latitudes, the change in oxygen is dominated by surface cooling and subsequent ventilation to mid-depths, while the carbon anomaly is associated with solubility changes and eruption-generated ENSO variability. Our results indicate that Pinatubo does not substantially impact oxygen or carbon in the Southern Ocean; forced signals do not emerge from the large internal variability in this region.