Marine carbon dioxide removal (mCDR) is gaining interest as a tool to meet global climate goals. Because the response of the ocean–atmosphere system to mCDR takes years to centuries, modeling is required to assess the impact of mCDR on atmospheric CO2reduction. Here, we use a coupled ocean–atmosphere model to quantify the atmospheric CO2reduction in response to a CDR perturbation. We define two metrics to characterize the atmospheric CO2response to both instantaneous ocean alkalinity enhancement (OAE) and direct air capture (DAC): the cumulative additionality (
Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher.
Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?
Some links on this page may take you to non-federal websites. Their policies may differ from this site.
-
Abstract α ) measures the reduction in atmospheric CO2relative to the magnitude of the CDR perturbation, while the relative efficiency (ϵ ) quantifies the cumulative additionality of mCDR relative to that of DAC. For DAC,α is 100% immediately following CDR deployment, but declines to roughly 50% by 100 years post-deployment as the ocean degasses CO2in response to the removal of carbon from the atmosphere. For instantaneous OAE,α is zero initially and reaches a maximum of 40%–90% several years to decades later, depending on regional CO2equilibration rates and ocean circulation processes. The global meanϵ approaches 100% after 40 years, showing that instantaneous OAE is nearly as effective as DAC after several decades. However, there are significant geographic variations, withϵ approaching 100% most rapidly in the low latitudes whileϵ stays well under 100% for decades to centuries near deep and intermediate water formation sites. These metrics provide a quantitative framework for evaluating sequestration timescales and carbon market valuation that can be applied to any mCDR strategy. -
Abstract The ocean's organic carbon export is a key control on atmospheric pCO2and stimulating this export could potentially mitigate climate change. We use a data‐constrained model to calculate the sensitivity of atmospheric pCO2to local changes in export using an adjoint approach. A perpetual enhancement of the biological pump's export by 0.1 PgC/yr could achieve a roughly 1% reduction in pCO2at average sensitivity. The sensitivity varies roughly 5‐fold across different ocean regions and is proportional to the difference between the mean sequestration time
τ seqof regenerated carbon and the response timeτ preof performed carbon, which is the reduction in the preformed carbon inventory per unit increase in local export production. Air‐sea CO2disequilibrium modulates the geographic pattern ofτ pre, causing particularly high sensitivities (2–3 times the global mean) in the Antarctic Divergence region of the Southern Ocean.Free, publicly-accessible full text available June 28, 2025 -
Abstract The ocean is a major carbon sink and takes up 25%–30% of the anthropogenically emitted CO2. A state‐of‐the‐art method to quantify this sink are global ocean biogeochemistry models (GOBMs), but their simulated CO2uptake differs between models and is systematically lower than estimates based on statistical methods using surface ocean
p CO2and interior ocean measurements. Here, we provide an in‐depth evaluation of ocean carbon sink estimates from 1980 to 2018 from a GOBM ensemble. As sources of inter‐model differences and ensemble‐mean biases our study identifies (a) the model setup, such as the length of the spin‐up, the starting date of the simulation, and carbon fluxes from rivers and into sediments, (b) the simulated ocean circulation, such as Atlantic Meridional Overturning Circulation and Southern Ocean mode and intermediate water formation, and (c) the simulated oceanic buffer capacity. Our analysis suggests that a late starting date and biases in the ocean circulation cause a too low anthropogenic CO2uptake across the GOBM ensemble. Surface ocean biogeochemistry biases might also cause simulated anthropogenic fluxes to be too low, but the current setup prevents a robust assessment. For simulations of the ocean carbon sink, we recommend in the short‐term to (a) start simulations at a common date before the industrialization and the associated atmospheric CO2increase, (b) conduct a sufficiently long spin‐up such that the GOBMs reach steady‐state, and (c) provide key metrics for circulation, biogeochemistry, and the land‐ocean interface. In the long‐term, we recommend improving the representation of these metrics in the GOBMs.Free, publicly-accessible full text available March 1, 2025 -
Trawling the seafloor can disturb carbon that took millennia to accumulate, but the fate of that carbon and its impact on climate and ecosystems remains unknown. Using satellite-inferred fishing events and carbon cycle models, we find that 55-60% of trawling-induced aqueous CO2is released to the atmosphere over 7-9 years. Using recent estimates of bottom trawling’s impact on sedimentary carbon, we found that between 1996-2020 trawling could have released, at the global scale, up to 0.34-0.37 Pg CO2yr-1to the atmosphere, and locally altered water pH in some semi-enclosed and heavy trawled seas. Our results suggest that the management of bottom-trawling efforts could be an important climate solution.
-
Abstract Chromophoric dissolved organic matter (CDOM) is an important part of ocean carbon biogeochemistry with relevance to long‐term observations of ocean biology due to its dominant light absorption properties. Thus, understanding the underlying processes controlling CDOM distribution is important for predicting changes in light availability, primary production, and the cycling of biogeochemically important matter. We present a biogeochemical CDOM model for the open ocean with two classes of biological lability and uncertainty estimates derived from 43 ensemble members that provide a range of model parameter variations. Ensemble members were optimized to match global ocean in situ CDOM measurements and independently assessed against satellite CDOM estimates, which showed good agreement in spatial patterns. Based on the ensemble median, we estimate that about 7% of open‐ocean CDOM is of terrestrial origin, but the ensemble range is large (<0.1–26%). CDOM is rapidly removed in the surface ocean (<200 m) due to biological degradation for short‐lived CDOM and photodegradation for long‐lived CDOM, leading to a net flux of CDOM to the surface ocean from the dark ocean. This deep‐water source (ensemble median 0.001 m−1 yr−1) is similar in magnitude to the riverine flux (0.005 m−1 yr−1) into the surface ocean. Though discrepancies between the model and observational data remain, this work serves as a foundational framework for a mechanistic assessment of global CDOM distribution that is independent of satellite data.
-
The ocean holds vast quantities of carbon that it continually exchanges with the atmosphere through the air-sea interface. Because of its enormous size and relatively rapid exchange of carbon with the atmosphere, the ocean controls atmospheric CO 2 concentration and thereby Earth's climate on timescales of tens to thousands of years. This review examines the basic functions of the ocean's carbon cycle, demonstrating that the ocean carbon inventory is determined primarily by the mass of the ocean, by the chemical speciation of CO 2 in seawater, and by the action of the solubility and biological pumps that draw carbon into the ocean's deeper layers, where it can be sequestered for decades to millennia. The ocean also plays a critical role in moderating the impacts of climate change by absorbing an amount of carbon equivalent to about 25% of anthropogenic CO 2 emissions over the past several decades. However, this also leads to ocean acidification and reduces the chemical buffering capacity of the ocean and its future ability to take up CO 2 . This review closes with a look at the uncertain future of the ocean carbon cycle and the scientific challenges that this uncertainty brings.more » « less
-
Abstract This study characterized ocean biological carbon pump metrics in the second iteration of the REgional Carbon Cycle Assessment and Processes (RECCAP2) project. The analysis here focused on comparisons of global and biome‐scale regional patterns in particulate organic carbon (POC) production and sinking flux from the RECCAP2 ocean biogeochemical model ensemble against observational products derived from satellite remote sensing, sediment traps, and geochemical methods. There was generally good model‐data agreement in mean large‐scale spatial patterns, but with substantial spread across the model ensemble and observational products. The global‐integrated, model ensemble‐mean export production, taken as the sinking POC flux at 100 m (6.08 ± 1.17 Pg C yr−1), and export ratio defined as sinking flux divided by net primary production (0.154 ± 0.026) both fell at the lower end of observational estimates. Comparison with observational constraints also suggested that the model ensemble may have underestimated regional biological CO2drawdown and air‐sea CO2flux in high productivity regions. Reasonable model‐data agreement was found for global‐integrated, ensemble‐mean sinking POC flux into the deep ocean at 1,000 m (0.65 ± 0.24 Pg C yr−1) and the transfer efficiency defined as flux at 1,000 m divided by flux at 100 m (0.122 ± 0.041), with both variables exhibiting considerable regional variability. The RECCAP2 analysis presents standard ocean biological carbon pump metrics for assessing biogeochemical model skill, metrics that are crucial for further modeling efforts to resolve remaining uncertainties involving system‐level interactions between ocean physics and biogeochemistry.
-
Abstract Mid-depth North Pacific waters are rich in nutrients and respired carbon accumulated over centuries. The rates and pathways with which these waters exchange with the surface ocean are uncertain, with divergent paradigms of the Pacific overturning: one envisions bottom waters upwelling to 1.5 km depth; the other confines overturning beneath a mid-depth Pacific shadow zone (PSZ) shielded from mean advection. Here global inverse modelling reveals a PSZ where mean ages exceed 1400 years with overturning beneath. The PSZ is supplied primarily by Antarctic and North-Atlantic ventilated waters diffusing from below and from the south. Half of PSZ waters re-surface in the Southern Ocean, a quarter in the subarctic Pacific. The abyssal North Pacific, despite strong overturning, has mean re-surfacing times also exceeding 1400 years because of diffusion into the overlying PSZ. These results imply that diffusive transports – distinct from overturning transports – are a leading control on Pacific nutrient and carbon storage.more » « less
-
Abstract Barystatic sea level rise (SLR) caused by the addition of freshwater to the ocean from melting ice can in principle be recorded by a reduction in seawater salinity, but detection of this signal has been hindered by sparse data coverage and the small trends compared to natural variability. Here, we develop an autoregressive machine learning method to estimate salinity changes in the global ocean from 2001 to 2019 that reduces uncertainties in ocean freshening trends by a factor of four compared to previous estimates. We find that the ocean mass rose by 13,000 ± 3,000 Gt from 2001 to 2019, implying a barystatic SLR of 2.0 ± 0.5 mm/yr. Combined with SLR of 1.3 ± 0.1 mm/yr due to ocean thermal expansion, these results suggest that global mean sea level rose by 3.4 ± 0.6 mm/yr from 2001 to 2019. These results provide an important validation of remote‐sensing measurements of ocean mass changes, global SLR, and global ice budgets.