Abstract:Coupled physical–biogeochemical models can fill the spatial and temporal gap in ocean carbon observations. Challenges of applying a coupled physical–biogeochemical model in the regional ocean include the reasonable prescription of carbon model boundary conditions, lack of in situ observations, and the oversimplification of certain biogeochemical processes. In this study, we applied a coupled physical–biogeochemical model (Regional Ocean Modelling System, ROMS) to the Gulf of Mexico (GoM) and achieved an unprecedented 20-year high-resolution (5 km, 1/22∘ ) hindcast covering the period of 2000 to 2019. The biogeochemical model incorporated the dynamics of dissolved organic carbon (DOC) pools and the formation and dissolution of carbonate minerals. The biogeochemical boundaries were interpolated from NCAR's CESM2-WACCM-FV2 solution after evaluating the performance of 17 GCMs in the GoM waters. Model outputs included carbon system variables of wide interest, such as p CO 2 , pH, aragonite saturation state ( Ω Arag ), calcite saturation state ( Ω Calc ), CO 2 air–sea flux, and carbon burial rate. The model's robustness is evaluated via extensive model–data comparison against buoys, remote-sensing-based machine learning (ML) products, and ship-based measurements. A reassessment of air–sea CO 2 flux with previous modeling and observational studies gives us confidence that our model provides a robust and updated CO 2 flux estimation, and NGoM is a stronger carbon sink than previously reported. Model results reveal that the GoM water has been experiencing a ∼ 0.0016 yr −1 decrease in surface pH over the past 2 decades, accompanied by a ∼ 1.66 μ atm yr −1 increase in sea surface p CO 2 . The air–sea CO 2 exchange estimation confirms in accordance with several previous models and ocean surface p CO 2 observations that the river-dominated northern GoM (NGoM) is a substantial carbon sink, and the open GoM is a carbon source during summer and a carbon sink for the rest of the year. Sensitivity experiments are conducted to evaluate the impacts of river inputs and the global ocean via model boundaries. The NGoM carbon system is directly modified by the enormous carbon inputs ( ∼ 15.5 Tg C yr −1 DIC and ∼ 2.3 Tg C yr −1 DOC) from the Mississippi–Atchafalaya River System (MARS). Additionally, nutrient-stimulated biological activities create a ∼ 105 times higher particulate organic matter burial rate in NGoM sediment than in the case without river-delivered nutrients. The carbon system condition of the open ocean is driven by inputs from the Caribbean Sea via the Yucatan Channel and is affected more by thermal effects than biological factors.

A Numerical reassessment of the Gulf of Mexico carbon system in connection with the Mississippi River and global ocean

Marine biogeochemical cycling and oceanic CO&lt;sub&gt;2&lt;/sub&gt; uptake simulated by the NUIST Earth System Model version 3

Assessment of Global Ocean Biogeochemistry Models for Ocean Carbon Sink Estimates in RECCAP2 and Recommendations for Future Studies

Marine Biogeochemical Cycling and Oceanic CO2 Uptake Simulated by the NUIST Earth System Model Version 3 (NESM V3)

Consistency and Challenges in the Ocean Carbon Sink Estimate for the Global Carbon Budget

A regional hindcast model simulating ecosystem dynamics, inorganic carbon chemistry and ocean acidification in the Gulf of Alaska

Mississippi River Chemistry Impacts on the Interannual Variability of Aragonite Saturation State in the Northern Gulf of Mexico

An Assessment of CO2 Storage and Sea‐Air Fluxes for the Atlantic Ocean and Mediterranean Sea Between 1985 and 2018

Early Diagenesis in the Hypoxic and Acidified Zone of the Northern Gulf of Mexico: Is Organic Matter Recycling in Sediments Disconnected From the Water Column?

Seasonal Mixing and Biological Controls of the Carbonate System in a River-Dominated Continental Shelf Subject to Eutrophication and Hypoxia in the Northern Gulf of Mexico

The Role of Benthic Fluxes in Acidifying the Bottom Waters in the Northern Gulf of Mexico Hypoxic Zone Based on an Updated Water Column Biogeochemical‐Seabed Diagenetic and Sediment Transport Model

An Earth System Model of Intermediate Complexity: Simulation of the Role of Ocean Mixing Parameterizations and Climate Change in Estimated Uptake for Natural and Bomb Radiocarbon and Anthropogenic Co2

Evaluation of CMIP models with IOMB: Rates of contemporary ocean carbon uptake linked with vertical temperature gradients and transport to the ocean interior

Microbial metabolite fluxes in a model marine anoxic ecosystem

Learning about the Ocean Carbon Cycle from Observational Constraints and Model Simulations of Multiple Tracers

Spatial and temporal variation in surface nitrate and phosphate in the Northern Gulf of Mexico over 35 years

Marine biogeochemical cycling and oceanic CO 2 uptake simulated by the NUIST Climate System Model : NUIST-CSM-2 . 0 . 1

Gulf of Mexico Hypoxia: Exploring Increasing Sensitivity to Nitrogen Loads

Variability in the Global Ocean Carbon Sink From 1959 to 2020 by Correcting Models With Observations

Dissolved inorganic carbon and alkalinity fluxes from coastal marine sediments: model estimates for different shelf environments and sensitivity to global change

Enhanced CO 2 uptake of the coastal ocean is dominated by biological carbon fixation