The sensitivity of estuarine aragonite saturation state and pH to the carbonate chemistry of a freshet-dominated river

Ocean acidification threatens to reduce pH and aragonite saturation state (ΩA) in estuaries, potentially damaging their ecosystems. However, the impact of highly variable river total alkalinity (TA) and dissolved inorganic carbon (DIC) on pH and ΩA in these estuaries is unknown. We assess the sensitivity of estuarine surface pH and ΩA to river TA and DIC using a coupled biogeochemical model of the Strait of Georgia on the Canadian Pacific coast and place the results in the context of global rivers. The productive Strait of Georgia estuary has a large, seasonally variable freshwater input from the glacially fed, undammed Fraser River. Analyzing TA observations from this river plume and pH from the river mouth, we find that the Fraser is moderately alkaline (TA 500–1000 µmol kg−1) but relatively DIC-rich. Model results show that estuarine pH and ΩA are sensitive to freshwater DIC and TA, but do not vary in synchrony except at high DIC : TA. The asynchrony occurs because increased freshwater TA is associated with increased DIC, which contributes to an increased estuarine DIC : TA and reduces pH, while the resulting higher carbonate ion concentration causes an increase in estuarine ΩA. When freshwater DIC : TA increases (beyond  ∼  1.1), the shifting chemistry causes a paucity of the carbonate ion that overwhelms the simple dilution/enhancement effect. At this high DIC : TA ratio, estuarine sensitivity to river chemistry increases overall. Furthermore, this increased sensitivity extends to reduced flow regimes that are expected in future. Modulating these negative impacts is the seasonal productivity in the estuary which draws down DIC and reduces the sensitivity of estuarine pH to increasing DIC during the summer season

Concentrations and cycling of DMS, DMSP, and DMSO in coastal and offshore waters of the Subarctic Pacific during summer, 2010-2011

Author Posting. © American Geophysical Union, 2017. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 122 (2017): 3269–3286, doi:10.1002/2016JC012465.Concentrations of dimethylsulfide (DMS), measured in the Subarctic Pacific during summer 2010 and 2011, ranged from ∼1 to 40 nM, while dissolved dimethylsulfoxide (DMSO) concentrations (range 13-23 nM) exceeded those of dissolved dimethyl sulfoniopropionate (DMSP) (range 1.3–8.8 nM). Particulate DMSP dominated the reduced sulfur pool, reaching maximum concentrations of 100 nM. Coastal and off shore waters exhibited similar overall DMS concentration ranges, but sea-air DMS fluxes were lower in the oceanic waters due to lower wind speeds. Surface DMS concentrations showed statistically significant correlations with various hydrographic variables including the upwelling intensity (r2 = 0.52, p < 0.001) and the Chlorophyll a/mixed layer depth ratio (r2 = 0.52, p < 0.001), but these relationships provided little predictive power at small scales. Stable isotope tracer experiments indicated that the DMSP cleavage pathway always exceeded the DMSO reduction pathway as a DMS source, leading to at least 85% more DMS production in each experiment. Gross DMS production rates were positively correlated with the upwelling intensity, while net rates of DMS production were significantly correlated to surface water DMS concentrations. This latter result suggests that our measurements captured dominant processes driving surface DMS accumulation across a coastal-oceanic gradient.Natural Sciences and Engineering Research Council of Canada, from the Peter Wall Institute for Advanced Studies2017-10-2

Estimating marine carbon uptake in the northeast Pacific using a neural network approach

The global ocean takes up nearly a quarter of anthropogenic CO2 emissions annually, but the variability in this uptake at regional scales remains poorly understood. Here we use a neural network approach to interpolate sparse observations, creating a monthly gridded seawater partial pressure of CO2 (pCO2) data product from January 1998 to December 2019, at 1/12∘ × 1/12∘ spatial resolution, in the northeast Pacific open ocean, a net sink region. The data product (ANN-NEP; NCEI Accession 0277836) was created from pCO2 observations within the 2021 version of the Surface Ocean CO2 Atlas (SOCAT) and a range of predictor variables acting as proxies for processes affecting pCO2 to create nonlinear relationships to interpolate observations at a spatial resolution 4 times greater than leading global products and with better overall performance. In moving to a higher resolution, we show that the internal division of training data is the most important parameter for reducing overfitting. Using our pCO2 product, wind speed, and atmospheric CO2, we evaluate air–sea CO2 flux variability. On sub-decadal to decadal timescales, we find that the upwelling strength of the subpolar Alaskan Gyre, driven by large-scale atmospheric forcing, acts as the primary control on air–sea CO2 flux variability (r2=0.93, p&lt;0.01). In the northern part of our study region, divergence from atmospheric CO2 is enhanced by increased local wind stress curl, enhancing upwelling and entrainment of naturally CO2-rich subsurface waters, leading to decade-long intervals of strong winter outgassing. During recent Pacific marine heat waves from 2013 on, we find enhanced atmospheric CO2 uptake (by as much as 45 %) due to limited wintertime entrainment. Our product estimates long-term surface ocean pCO2 increase at a rate below the atmospheric trend (1.4 ± 0.1 µatm yr−1) with the slowest increase in the center of the subpolar gyre where there is strong interaction with subsurface waters. This mismatch suggests the northeast Pacific Ocean sink for atmospheric CO2 may be increasing.</p

Eddy-resolving simulation of plankton ecosystem dynamics in the California Current System

Author Posting. © Elsevier B.V., 2006. This is the author's version of the work. It is posted here by permission of Elsevier B.V. for personal use, not for redistribution. The definitive version was published in Deep Sea Research Part I: Oceanographic Research Papers 53 (2006): 1483-1516, doi:10.1016/j.dsr.2006.06.005.We study the dynamics of the planktonic ecosystem in the coastal upwelling zone within the California Current System using a three-dimensional, eddy-resolving circulation model coupled to an ecosystem/biogeochemistry model. The physical model is based on the Regional Oceanic Modeling System (ROMS), configured at a resolution of 15 km for a domain covering the entire U.S. West Coast, with an embedded child grid covering the central California upwelling region at a resolution of 5 km. The model is forced with monthly mean boundary conditions at the open lateral boundaries as well as at the surface. The ecological/biogeochemical model is nitrogen based, includes single classes for phytoplankton and zooplankton, and considers two detrital pools with different sinking speeds. The model also explicitly simulates a variable chlorophyll-to-carbon ratio. Comparisons of model results with either remote sensing observations (AVHRR, SeaWiFS) or in situ measurements from the CalCOFI program indicate that our model is capable of replicating many of the large-scale, time averaged features of the coastal upwelling system. An exception is the underestimation of the chlorophyll levels in the northern part of the domain, perhaps because of the lack of short-term variations in the forcing from the atmosphere. Another shortcoming is that the modeled thermocline is too diffuse, and that the upward slope of the isolines toward the coast is too small. Detailed time-series comparisons with observations from Monterey Bay reveal similar agreements and discrepancies. We attribute the good agreement between the modeled and observed ecological properties in large part to the accuracy of the physical fields. In turn, many of the discrepancies can be traced back to our use of monthly mean forcing. Analysis of the ecosystem structure and dynamics reveal that the magnitude and pattern of phytoplankton biomass in the nearshore region are determined largely by the balance of growth and zooplankton grazing, while in the offshore region, growth is balanced by mortality. The latter appears to be inconsistent with in situ observations and is a result of our consideration of only one zooplankton size class (mesozooplankton), neglecting the importance of microzooplankton grazing in the offshore region. A comparison of the allocation of nitrogen into the different pools of the ecosystem in the 3-D results with those obtained from a box model configuration of the same ecosystem model reveals that only a few components of the ecosystem reach a local steady-state, i.e. where biological sources and sinks balance each other. The balances for the majority of the components are achieved by local biological source and sink terms balancing the net physical divergence, confirming the importance of the 3-D nature of circulation and mixing in a coastal upwelling system.Most of this work has been made possible by two grants from NASA. Additional support is acknowledged from NSF’s ITR program

A vertical model of particle size distributions and fluxes in the midwater column that includes biological and physical processes - Part I: model formulation

The downward transport of surface particle production constitutes an important mechanism for carbon sequestration by the ocean. Only a small fraction (similar or equal to 10%) of the flux that leaves the euphotic zone reaches 1000 m depth because the particulate organic matter is consumed and transformed in the oceanic midwater column. The depth at which this transformation occurs is crucial to estimate carbon sequestration. Description of the particle flux and remineralization with depth below the euphotic zone has previously been limited to empirical relationships that neglect physical and biological mechanisms. Because several particle properties and functions (settling speed, rates of coagulation and consumption) are related to particle size, measurements of particle size spectra provide an important tool to understand particle dynamics. Most of the mechanisms affecting particle dynamics known for the surface layer should continue in the mesopelagic region. We review the different mechanisms and formulate equations describing changes in particle size distributions throughout the water column as a result of particle sinking, coagulation, disaggregation, and bacterial and zooplankton consumption. The resulting model describes the midwater particle population by its mass distribution for the size range I pm to I cm. Most of the particle data sets have a narrower size range but we suggest that model results are not greatly affected by the lack of data on particles smaller than 200 mum because most of the particle mass is in particles larger than 100 mum, as seen in a full size-range particle spectrum. The combination of this model with a unique particle size spectra data set, obtained during an inter-annual survey at the French JGOFS site in the NW Mediterranean Sea, give an insight into the key processes for particle dynamics in the unknown midwater layers. (C) 2004 Elsevier Ltd. All rights reserved