The Southern California Bight (SCB), an eastern boundary upwelling system, is impacted by global warming, acidification and deoxygetation, and receives anthropogenic nutrients from a coastal population of 20 million people. We describe the configuration, forcing, and validation of a realistic, submesoscale resolving ocean model as a tool to investigate coastal eutrophication. This modeling system represents an important achievement because it strikes a balance of capturing the forcing by U.S. Pacific Coast-wide phenomena, while representing the bathymetric features and submesoscale circulation that affect the vertical and horizontal transport of nutrients from natural and human sources. Moreover, the model allows to run simulations at timescales that approach the interannual frequencies of ocean variability, making the grand challenge of disentangling natural variability, climate change, and local anthropogenic forcing a tractable task in the near-term. The model simulation is evaluated against a broad suite of observational data throughout the SCB, showing realistic depiction of mean state and its variability with remote sensing and in situ physical-biogeochemical measurements of state variables and biogeochemical rates. The simulation reproduces the main structure of the seasonal upwelling front, the mean current patterns, the dispersion of plumes, as well as their seasonal variability. It reproduces the mean distributions of key biogeochemical and ecosystem properties. Biogeochemical rates reproduced by the model, such as primary productivity and nitrification, are also consistent with measured rates. Results of this validation exercise demonstrate the utility of fine-scale resolution modeling in support of management decisions on local anthropogenic nutrient discharges to coastal zones.

Continental margins play an essential role in global ocean biogeochemistry and the carbon cycle; however, global assessments of this role remain highly uncertain. This uncertainty arises from large variability over a broad range of temporal and spatial scales of the processes that characterize these environments. High-resolution simulations with ocean biogeochemical models have emerged as essential tools to advance biogeochemical assessments at regional scales. Here, we examine the processes and balances for carbon, oxygen, and nitrogen cycles along the U.S. West Coast in an 11-year hindcast simulation with a submesoscale-permitting oceanic circulation-biogeochemical model. We highlight the importance of biogeochemical cycles on the continental shelf, and their connection to the broader regional context encompassing the California Current System. On the shelf, coastal and wind stress curl upwelling drive a vigorous overturning circulation that supports biogeochemical rates and fluxes that are approximately twice as large as offshore. Exchanges with the proximate sediments, submesoscale shelf currents, bottom boundary layer transport, and intensified cross-shelf export of shelf-produced materials impact coastal and open-ocean balances. While regional variability prevents extrapolation of our results to global margins, our approach provides a powerful tool to identify the dominant dynamics in different shelf setting and quantify their large-scale consequences.

Earth System Models project a decline of dissolved oxygen in the oceans under warming climate. Observational studies suggest that the ratio of O2 inventory to ocean heat content (O2-OHC) is several fold larger than can be explained by solubility alone, but the ratio remains poorly understood. In this work, models of different complexity are used to understand the factors controlling the O2-OHC ratio during deep convection, with a focus on the Labrador Sea, a site of deep water formation in the North Atlantic Ocean. A simple one-dimensional convective adjustment model suggests two limit case scenarios. When the near-surface oxygen level is dominated by the entrainment of subsurface water, surface buoyancy forcing, air-sea gas exchange coefficient and vertical structure of sea water together affect the O2-OHC ratio. In contrast, vertical gradients of temperature and oxygen become important when the surface oxygen flux dominates. The former describes the O2-OHC ratio of individual convective event in agreement with model simulations of deep convection. The latter captures the O2-OHC ratio of interannual variability, where the pre-conditioning of interior ocean gradients dominates. The relative vertical gradients of temperature and oxygen, which in turn depend on the lateral transport and regional biological productivity, control the year-to-year variations of O2-OHC ratio. These theoretical predictions are tested against the output of a three-dimensional regional circulation and biogeochemistry model which captures the observed large-scale distribution of the O2-OHC ratio, and agrees broadly with the prediction by the simpler model.