Articles | Volume 20, issue 14
Research article
31 Jul 2023
Research article |  | 31 Jul 2023

Influence of manganese cycling on alkalinity in the redox stratified water column of Chesapeake Bay

Aubin Thibault de Chanvalon, George W. Luther, Emily R. Estes, Jennifer Necker, Bradley M. Tebo, Jianzhong Su, and Wei-Jun Cai

The alkalinity dynamics in coastal environments play a crucial role in controlling the global burial of carbonate minerals and the ocean's capacity to sequester anthropogenic CO2. This study presents results from high vertical resolution profiles obtained during two summers in the temperate Chesapeake Bay estuary, enabling detailed investigation of carbonate dynamics over salinity and redox gradients, along with measurement of the speciation of most redox-sensitive elements. Under oxygen-rich conditions, carbonate dissolution, primary production and aerobic respiration explain the evolution of total alkalinity (TA) versus dissolved inorganic carbon (DIC), once adjusted for fresh and oceanic water mixing. A significant flooding event in 2018 promoted carbonate dissolution. In oxygen-depleted waters, we observed a previously unreported 2.4 mol increase in DIC per 1 mol of TA production, which was consistent over the 2 years. Stoichiometric changes suggest that MnO2 reduction followed by Mn carbonate precipitation is responsible for this characteristic carbonate signature, likely produced in sediment pore water and then transferred to the water column along with other by-products of anoxic respiration at the onset of summer. Our findings highlight the critical role of Mn in alkalinity dynamics in the Chesapeake Bay and potentially other river-dominated environments where it can limit H2S oxidation to SO42- and promote sulfur burial.

1 Introduction

About 30 % of anthropogenic CO2 emissions are rapidly trapped by dissolution in the ocean as dissolved inorganic carbon (DIC), which is dominated by bicarbonate ions (HCO3-; Friedlingstein et al., 2019). At the century timescale, atmosphere–ocean exchanges result in oceanic HCO3- enrichment not associated with a cationic enrichment, as measured by total alkalinity (TA). This is in contrast to silicate or carbonate weathering, which preponderates at timescales of thousands to millions of years and increases TA simultaneously with DIC (Urey, 1952). This cationic deficiency corresponds with an excess of proton release that can be tracked by an indicator such as TA minus DIC (Xue and Cai, 2020) and informs about the reduction of the CO2 buffering capacity of seawater (Gustafsson et al., 2019). The proton excess is only balanced in the deep ocean by the production of Ca2+ from carbonate dissolution, a process named chemical carbonate compensation (Boudreau et al., 2018). However, in shallow waters, that accounts for two-thirds of global buried carbonate (Smith and Mackenzie, 2016); carbonate precipitation largely predominates over dissolution, and other localised processes may constrain carbonate dynamics (Borges et al., 2006; Lohrenz et al., 2010). For example, some calcifying species may slow down their carbonate precipitation in case of pH decrease (so-called biological carbonate compensation; Boudreau et al., 2018). As another example, in estuaries, the seasonality of river flow, temperature and continental erosion modulates CaCO3 dissolution (e.g. Su et al., 2020b), atmospheric CO2 exchanges (e.g. Borges et al., 2018), respiratory activity (e.g. Abril et al., 2003) and transfer of carbonate particles (Meybeck, 1987; Middelburg et al., 2020).

The global trends of human migration towards littoral areas and global warming favour eutrophication and a decrease in oxygen levels in coastal water (Breitburg et al., 2018; Rabalais et al., 2014). Sometimes driven by recurrent natural process (e.g. Gupta et al., 2021), exceptional events (e.g. Hulot et al., 2023), stratification reinforced by global warming (e.g. Meire et al., 2013) or anthropogenic nutrient loading (e.g. Hagy et al., 2004; Carstensen et al., 2018), coastal deoxygenation enhanced the possible build-up of anoxic conditions to build up as a permanent or seasonal feature. In anoxic environments, the drivers of the carbonate cycle change as the terminal electron acceptor for organic matter remineralisation shifts from O2 to Mn oxides, Fe oxides and then sulfate. Notably, the consumption of sulfate, followed by sulfur burial, produces an significant alkalinity flux toward coastal waters, accounting for approximately 10 % of the global alkalinity input to the ocean (Middelburg et al., 2020). Fe and Mn oxides play a critical role in sulfur burial, as they prevent SO42- regeneration and lead to the production of S0 or FeS instead (Findlay et al., 2014; Avetisyan et al., 2021), leaving a specific fingerprint in the TA and DIC effluxes.

While TA controls the CO2 buffering capacity of the ocean, riverine input of carbonate to the ocean is poorly constrained (Middelburg et al., 2020), and only rare publications take into account the estuarine transformations of the carbonate species (e.g. Su et al., 2020a; Abril et al., 2003) furthermore in the context of oxygen depletion (e.g. Abril et al., 2004). To better constrain the carbonate cycle in oxygen-depleted estuaries, we sampled a stratified water column in the Chesapeake Bay 11 times over two campaigns with a high vertical resolution (down to 10 cm). This protocol allows a precise description of carbonate dynamics over a redox gradient along with the measurement of the speciation of most redox-sensitive elements. Such sampling illustrates carbonate chemistry on transitioning from oxygenated waters to waters devoid of oxygen, as usually only encountered in sediments or in anoxic lakes or seas (e.g. the Black Sea). The original observed changes of alkalinity versus dissolved inorganic carbon changes are interpreted based on typical geochemical reactions occurring along the redox gradient.

2 Materials and methods

2.1 Sampling

During two sampling campaigns from 3 to 9 August 2017 and 28 July to 3 August 2018, 11 profile casts were conducted in a single station in the Chesapeake Bay with a water depth of 25 m (station 858; 3858.54 N, 07622.22 W). The Susquehanna River is the main tributary of the bay representing on average two-thirds of the fresh water input (Zhang et al., 2015). Despite the similar season, the two campaigns occurred at very different river flow, with about 850 m3 s−1 in 2017 versus 8500 m3 s−1 in 2018 due to release of flood waters from the Conowingo Dam. The August 2018 condition corresponds to flooding which occurs on average every 3.5 years (return period of 3.5 years, USGS survey).

Each CTD (conductivity, temperature, depth) cast was performed during low or high tide slack periods. An oxygen sensor (Clark electrode, SBE Inc.; detection limit of 1 µM) and fluorescence sensor (ECO-FL Fluorometer, WET Labs) were part of the CTD Rosette to take measurements during sampling. Also, a submersible all-plastic pump profiler was attached with the pump near the sensor orifices, allowing measurement and sampling at a resolution of a few centimetres over 25 m water depth. Water was pumped to the deck within 1 min, and water passed through a flow through voltammetry system measuring O2, Mn(II), Fe(II), organically complexed Fe(III), FeS clusters, H2S and polysulfides continuously (Hudson et al., 2019). When redox interfaces were identified, samples were filtered through an acetate cartridge filter (pore size 0.45 µm) for pH and inorganic carbon parameters, which were processed onboard within a few hours after sampling in order to conserve chemical speciation. The pump profiler system was cleaned with deionised water (18 ) onboard the deck of the ship after deployment. No coating effects were observed with the pump system.

2.2 Discrete measurements

For each sample, all redox species were determined in the through flow voltammetry system using cyclic voltammetry with a 100 µm diameter Au/Hg amalgam PEEK microelectrode prepared according to Luther et al. (2008) connected to a DLK-60 electrochemical analyser from Analytical Instrument Systems Inc. The detection limit of this method is 0.2 µM for sulfide and polysulfides. Discrete samples for the determination of NO2-, Fe and Mn species were filtered through nylon luer-lock syringe filters (Millipore, 0.20 µm). Iron was measured on both bulk and filtered samples using the ferrozine method (Stookey, 1970): after HCl acidification (for Fe(II)) and an optional reduction step (for Fe(III)+Fe(II)) with hydroxylamine hydrochloric (final concentration 0.7 M for 1 h), ammonium acetate (final concentration 0.5 M) and ferrozine (final concentration 1 mM). Absorbance at 562 nm was read with a diode array spectrophotometer (Hewlett Packard 8452B). Limit of detection is 100 nM for Fe(II) and Fe(III) with a 1 cm cell. Shipboard nitrite determination was performed using the method of Grasshoff (1983). To 25 mL of sample, 0.5 mL of 58 mM sulfanilamide in 10 % v/v HCl and 0.5 mL of a 4 mM N(1-naphthyl)ethylene diamine hydrochloride solution were added. Samples with added reagents were shaken and left to sit for 15 min, followed by UV–Vis analysis at 540 nm using a 10 cm cell to increase detection limits. Calibration curves were constructed using sodium nitrite. Limit of detection is 10 nM for NO2-.

Dissolved manganese was determined by displacement of a Cd(II)-porphyrin complex with Mn(II) to form the Mn(II)-porphyrin complex (Ishii et al., 1982). Mn(III) species were identified based on slower reactivity with the Cd complex (Madison et al., 2011) as modified in Thibault de Chanvalon and Luther (2019). Alternatively, Mn(III) species were identified after HCl treatment (down to pH=1.5) followed by filtration in order to flocculate and eliminate the dissolved manganese bound to humic material by filtration (Oldham et al., 2017b). Limit of detection is 50 nM for Mn(II) and Mn(III) in a 1 cm cell. MnOx was measured on 20 mL samples of suspended material retained on 0.2 µm filters by the Leucoberbelin blue (LBB) method (Jones et al., 2019). Four millilitres of a reagent solution ([LBB]=78µM, [acetic acid]=14 mM) react with the filter, and the absorbance is read at 624 nm. KMnO4 was used to calibrate the LBB method, which allows the calculation of the electron equivalents obtained from particulate MnOx. Results are given as MnO2 equivalent, with a limit of detection of 0.1 µM and an uncertainty below 5 %.

The DIC samples were preserved in 250 mL borosilicate glass bottles with 50 µL saturated HgCl2 solution. The total alkalinity (TA) samples were not poisoned to prevent HgS precipitation and H+ release in anoxic and low-salinity waters (Cai et al., 2017). Then TA was analysed by Gran titration in an open-cell setting (AS-ALK2, Apollo SciTech) within 24 h of collection (Cai et al., 2010). The DIC samples were measured by a nondispersive infrared analyser (AS-C3, Apollo SciTech) within a week (Huang et al., 2012). The precision for DIC and TA was about 0.1 %. Both DIC and TA measurements were calibrated against certified reference materials (CRMs Batch 163 and Batch 173 provided by Andrew Dickson of the Scripps Institution of Oceanography). The pH samples were measured onboard at 25 C within 1 h of collection using an Orion ROSS glass electrode and calibrated with NIST standard buffers. The pCO2, calcite saturation and TA were calculated from measured DIC and pH via the CO2SYS program using Cai and Wang (1998) constants. The measured TA was found to be highly correlated to the calculated TA (r2=0.995 and 0.998, slope=0.995 and 1.017 for the 2017 and 2018 campaign respectively), and their difference was always below 45 µM with an average of 7.5 µM for 2017 and of 22.2 µM for 2018. These results suggest a low contribution of non-carbonate species (e.g. nitrite, ammonium or organic matter; Cotovicz et al., 2016), and measured TA was used for the interpretation.

2.3 Models of biogeochemical process on TA and DIC

2.3.1 Identification of biogeochemical processes from scatter plots

This section identifies the required conditions to interpret a scatter plot of two species in terms of chemical reaction stoichiometry: the “reaction-driven” approximation. In case of turbulent diffusion mixing (sometimes called eddy diffusion) in only one direction (no lateral input), at steady state, on a portion of space where occurs one chemical reaction, the changes of concentration, C and D, of two species can be described by Eq. (1) (see Sect. A1 for more details):

(1) d C d x = α C α D d D d x + G ,

with αC and αD the respective stoichiometric coefficient, x the direction in the space and G an unknown constant. Here, C and D are concentrations of any elements produced or consumed during the reaction and can be TA, DIC or a linear combination with salinity. When G is null, Eq. (2) is valid and the data points would be distributed along a straight line in a C versus D scatter plot.

(2) d C d D = α C α D

Theoretically, Eq. (2) is valid in at least three cases. First, it is valid if the contribution from one endmember is negligible compared to the timescale of reactions, which corresponds also to a system with only one endmember. Second, it is valid when the slope of ΔC/ΔD between the endmembers equals αC/αD, which is achieved, in particular, if one endmember was previously generated from the other endmember by a chemical reaction with similar stoichiometry. And third, it is valid when the two endmembers have similar concentrations (see Sect. A2), as in this study where the excess of TA (TAex) and the excess of DIC (DICex) are calculated by linear combination with salinity to be equal to zero for the upstream and downstream endmembers. Thus, if only one reaction occurs in between the endmembers, the slope ΔTAex/ΔDICex would correspond to the stoichiometry of αTA/αDIC. The oceanic endmember for both August campaigns was the one proposed by Su et al. (2020a). The oceanic endmember varies mainly with the season (Cai et al., 2020), and a maximal change of 50 µM results in 5 % uncertainty on the slope of the mixing line. For the upstream estuary, larger variations exist mainly due to changes of river discharge and water temperature (Meybeck, 2003; Joesoef et al., 2017), and a one-off endmember has to be determined by fit with the in situ measurements at the lowest measured salinity. Between the two campaigns, the upstream endmember changed by 77 µM, generating 5 % of change on the slope (see Fig. A1). The upstream endmember is not a river endmember (Su et al., 2020a) but corresponds to a salinity above 1.5, avoiding any interpretation for biological activity in the fresh water part of the estuary (Meybeck et al., 1988). However, it corresponds to a larger water mass pool, less sensitive to short term changes, with a residence time being higher than 240 d in the Chesapeake upstream part (Du and Shen, 2016), and thus is more likely to satisfy the condition of steady state.

The relative uncertainty on ΔTAex/ΔDICex is equal to the sum of the relative uncertainty of ΔTAex and ΔDICex. Posing ΔTA as the change of TA measured, ΔS as the change of salinity and sml_TA as the slope of the mixing line for TA, we have ΔTAex=ΔTA-sml_TA×ΔS. Uncertainty on ΔTA and ΔS are negligible to the relative uncertainty of the slope of the mixing line. Posing δ(x) as the uncertainty on x, we get Eq. (3) that describes the fact that the uncertainty is much lower on ΔDICex than on DICex because most of the error associated with the calculation of the endmember is cancelled when calculating the difference of DICex on two points with close salinity:

(3) δ ( Δ TAex / Δ DICex ) Δ TAex / Δ DICex = δ ( Δ TAex ) Δ TAex + δ ( Δ DICex ) Δ DICex = δ ( sml_TA ) sml_TA + δ ( sml_DIC ) sml_DIC = 0.1 .

However, in a stratified water column, not just one but several successive reactions occur, limiting the validity of Eq. (1) to each reactional stratum. The general case is not straightforward to solve, but in the particular case where the C versus D plot represents a straight line between two endmembers with different concentrations, the previous discussion of Eq. (2) indicates (second case) that one endmember would have been previously generated from the second by a chemical reaction with similar stoichiometry. Thus, the depths corresponding to the straight line define a reactional stratum characterised by a constant αC/αD and delimitated by two local endmembers maintained in steady state by chemical reactions with similar stoichiometry than the one that produced them, i.e. ΔC/ΔD=αC/αD. The local endmembers should have been produced before the steady-state achievement by a reaction of similar stoichiometry, but the reaction could have been faster than the observed one or could have occurred in a different place, including in the sediment.

Therefore, in a system defined between only two endmembers, away from atmospheric exchanges, in case of turbulent diffusion mixing, at steady state and with negligible lateral mixing, the reaction-driven approximation allows us to interpret linear variations of TAex versus DICex as a sum of biogeochemical reactions spread all over the water column that can be broken into several reactional stratum. In each stratum, if the local ΔTAex/ΔDICex ratio is constant, it corresponds to the apparent stoichiometry of a combination of the biogeochemical reactions occurring in this stratum. In the case of multiple simultaneous reactions in the same stratum, by posing αCi as the stoichiometry of the ith reaction concerning the reagent C and vi the reaction rate of the ith reaction, we obtain


To maintain a global reaction rate independent of the species, we have

v=ivi,(4)so     αC=1viαCivi.

Equation (4) indicates that the apparent stoichiometry in a given stratum corresponds to the sum of the stoichiometric coefficients of each reaction weighted by the relative rate of each reaction. Therefore, to estimate the relative rate of each reaction to the observed local changes of TAex, DICex, and AOU (apparent oxygen uptake) or H2S, a linear combination of reactions is calculated. This combination has to fit three equations (one for each parameter), which allows a maximum of three reactions to be used to solve the system. A limited number of reactions is selected as candidates based on the discussion (see Table 1 and Sects. 4.2 and 4.3). Then, the system is solved with the minimum possible reactions, and the weighted coefficients, vi/v, are calculated.

Table 1Summary of the main net reactions occurring among the different zones of a redox gradient with reactants starting in equilibrium with the atmosphere, adapted from Soetaert et al. (2007). The calculations assume the Redfield ratio of the organic matter, i.e. γN=0.156 and γP=0.0094. Redox zones describe the different conditions required for completing the reaction.

Download Print Version | Download XLSX

2.3.2 TA changes indicated by the charge transfer during a reaction

The simplest way to calculate the TA changes induced by an individual reaction is to look at charge transfer induced by the stoichiometry of the given reaction. Indeed, the total alkalinity (TA) corresponds to the quantity of acid added to titrate a solution down to pH 4.5 (Dickson, 1981). It can be described by Eq. (5), with the example of HCl as an acid and B any titrated base.

(5) HCl + B - BH + Cl -

Assuming a complete reaction, the quantity of acid added is equal to the negative charges initially present in the sample consumed plus the positive charges added to the species initially present in the sample. Thus, the total alkalinity corresponds to the loss of negative charges (or gain of positive charges) for species initially present in the sample produced by the pH change from the initial pH = pHini to pH=4.5. Writing zipH, the charges held by species i, initially present in the sample at a given pH, we get Eq. (6):

(6) TA = i z i pH = 4.5 - i z i pHini .

However, the electroneutrality of water induces

izipHini=0(7)so     TA=izipH=4.5.

Equation (7) demonstrates that the total alkalinity is simply the sum of charges that each species present in the sample would have at pH 4.5. From Eq. (7), one can easily deduce the changes of alkalinity from any reaction stoichiometry as soon as the charges of all species at pH=4.5 are known. For example, for a 0.0020 M NaOH dissolution, its TA is 0.0020 M as at pH=4.5, TA=([Na+]+[H+]-[OH-])=0.002+10-4.5-10-14+4.5=0.0020 M. For natural waters, most of the time, the only acid/base species charged at pH=4.5 are H2PO4- and NH4+. In that respect, whatever the initial pH and the acid–base equilibrium of species in the sample is, the sum of phosphate species will count negatively, and the sum of ammonium species will count positively. Strictly speaking, at pH=4.5, acid species with pKa between 2.5 and 6.5, such as F and NO2-, would be only partially titrated, and the charge equals their concentration multiplied by a correction factor of (1+10pKa-4.5)-1, but this correction can be neglected to a first approximation. Equation (7) is equivalent to those published in Soetaert et al. (2007) or Wolf-Gladrow et al. (2007), whose Eq. (32) can be refined considering that


However, Eq. (7) is more concise and more general. For example, in suboxic water, specific species such as polysulfides (as HS82-, Rickard and Luther, 2007), and in highly productive environments, carboxylic groups from DOC, can be easily added as soon as the bearing charges at pH=4.5 are known.

3 Results

3.1 Water column stratification

High-resolution profiles of salinity and temperature plotted in Fig. 1 show the stratified water column at station 858. Carbonate and redox chemistry are plotted against salinity in Fig. 2. These data result from the 11 CTD casts performed over 1 week during each campaign and correspond to depths ranging from 0.7 to 25 m. Plots against depth generate noisier profiles and are shown in Fig. A2, while plots against salinity follow the water masses. Despite the overall lower salinity due to a near-7-fold river flow increase in 2018 than 2017 (5800 m s−1 in 2018 versus 850 m3 s−1 in 2017), similar zonation of the water column occurred.

Figure 1Superimposed salinity (in white) and temperature (in black) profiles over the 11 casts for each campaign.


Figure 2Superimposed carbonate and redox chemistry profiles of 11 casts done in August 2017 and in August 2018 at station 858 of the Chesapeake Bay.


A surface layer sampled only in 2018 is visible in the top 3 m depth (Fig. 1). It presents highly variable temperature and oxygen concentration and oversaturation of pCO2 (Fig. 2), indicating emission to the atmosphere. Below, at 3 m depth, a subsurface layer (named primary production zone or PP in Fig. 2) is characterised by a high amount of O2 (about or above 100 % saturation), high pH (about 8; 8.11 ± 0.07, n=13 in 2017 and 7.94 ± 0.08, n=14 in 2018) and high day-to-day temperature variation (above 1 C between different days). The layer presents relatively low pCO2 (505 ±75µatm, n=13 in 2017 and 770 ± 130 µatm, n=14 in 2018), with minimal values at 110 µatm in 2017 and 205 µatm in 2018 which are below the atmospheric pCO2 of 407 µatm (Chen et al., 2020). This signature corresponds to primary production (PP). Fluorescence (not shown) correlates with pH as expected for primary production (pH=Fluo(mV)×13+7.14, r2=0.8 in 2017 and r2=0.9 in 2018). DICex and TAex reach their minimal value in the surface layer. While DICex minimum is similar between the two campaigns (78 ± 17 µM, n=13 in 2017 and 76 ± 10 µM, n=14 in 2018); TAex minimum is much lower in 2017 (7 ± 3 µM, n=12) than in 2018 (48 ± 7 µM, n=14).

Below, with increasing salinity, an important increase of pCO2 accompanying the decrease of O2, pH and temperature is visible. A relatively invariable low oxygen zone (called ILO in Fig. 2) is here defined by the salinity invariance of O2 concentrations at about 24 ± 3 µM (n=15) in 2017 and 105 ± 5 µM (n=16) in 2018. Other species are also relatively stable for this depth, such as pCO2 at about 3250 ± 150 µatm (n=10) in 2017 and 2590 ± 100 µatm (n=16) in 2018 and pH at about 7.28 ± 0.01 (n=13) in 2017 and 7.37 ± 0.02 (n=16) in 2018. Deeper, where the oxygen is not detectable (≲1µM) and in absence of free sulfide, the so-called anoxic zone corresponds to a pH minimum at 7.24 ± 0.01 (n=40), similar in 2017 and 2018, that generates a pCO2 maximum. The deepest layer is a sulfidic layer, [H2S]=11.2± 2.8 µM (n=36), in which the pH seems quite stable at 7.33 ± 0.01 (n=29) and 7.32 ± 0.01 (n=9) in 2017 and 2018 respectively. The Ca2+ concentrations observed by Su et al. (2021) and during the 2018 cruise (data not shown) vary linearly with salinity (calcium excess stays below 200 µM or 10 % of total Ca). Assuming similar behaviour in 2017, calculations show that the whole water column (except four samples from the PP zone) is under saturated (0.36<Ωcal<1; mean=0.68) with respect to calcite in 2018, while undersaturation is only valid below S=10 in 2017. The main changes between the two campaigns correspond to a greater oxygen penetration in 2018, especially visible in the ILO zone. Additionally, a surface layer (with salinity below 3) is visible above the primary production zone in 2018 but is related to the more superficial sampling in 2018.

3.2 Intermediate redox species

The development of the anoxic zone during summer (Su et al., 2021) and the regularity of this development over the years (Sholkovitz et al., 1992; Trouwborst et al., 2006; Lewis et al., 2007; Cai et al., 2017; Oldham et al., 2017a) requires the presence of species able to rapidly oxidise the H2S mixing upward and to reduce the O2 mixing downward. The three main redox couples known to play this role, NO3-/NO2-, MnOx/Mn2+ and Fe3+/Fe2+, are described in Fig. 3 by the superimposition of all cast results against salinity. Four representative casts are plotted in Fig. A3. The primary production zone is depleted in dissolved Mn(II) and Mn(III) and shows an average value for solid MnOx (0.7 ± 0.1 µM, n=4) and NO2- (2.1 ± 0.1 µM, n=18). In 2018, iron speciation was investigated in the primary production zone, showing an important iron pool (1.6 ± 0.1 µM, n=9) dominated by solid Fe(II) (95 % ± 2 %). Below, in the ILO zone, NO2- reaches a maximum plateau at 5.3 ± 0.2 µM, n=11 in 2017 that is not visible in 2018; MnOx dominates the Mn pool in the ILO zone with a concentration of 1.8 ± 0.1 µM, n=14 in 2017 and 1.3 ± 0.2 µM, and n=3 in 2018. Fe(II) represents only 49 %± 13 %, n=4 of the iron pool (for a total iron concentration of 2.0 ± 0.4 µM, n=4) in 2017, while it represents 87 % ± 5 %, n=4 of the iron pool in 2018 (for a total concentration of 1.2 ± 0.2 µM, n=4). Just below oxygen depletion, in the anoxic zone, MnOx reaches a maximum (2.0 ± 0.1 µM, n=25 in 2017 and 2.9 ± 0.4 µM, n=12 in 2018). During both campaigns, the sulfidic layer is characterised by the absence of NO2- and MnOx, while dissolved manganese concentration increases up to 7.7 ± 0.8 µM, n=15, and the iron pool increases to 5.0 ± 0.7 µM, n=13. No Mn(III) was detected with the porphyrin kinetics method (Thibault de Chanvalon and Luther, 2019), but about 30 % of the total dissolved manganese flocculated after acidification down to pH 1.5 when analysed in 2017, indicating the existence of Mn(III).

Figure 3Mn, Fe and nitrite profiles built by superimposition of 11 casts for each campaign.


4 Discussion

In terms of salinity, station 858 of the Chesapeake Bay shows very similar water column zonation between summer 2017 and summer 2018, despite a 10-fold difference in freshwater discharge rates (Figs. 1 and 2). Major features are first, a surface layer characterised by intense atmospheric exchange that was only sampled in 2018. Below, at about 3 m depth, there is a subsurface layer associated with high primary production (PP zone) with high pH of about 8, oversaturation of dioxygen and low CO2 partial pressure (down to 110 µatm). Below, a low oxygen layer with invariant concentration of most species surveyed (the invariant low oxygen (ILO) zone) is characterised by significant nitrite accumulation in 2017 (Fig. 2), probably due to oxidation of NH4+ diffusing upward and/or produced by in situ remineralisation. This feature is not visible in 2018 probably because the higher O2 concentration in 2018 accelerates nitrite oxidation into nitrate and prevents any significant accumulation. Below, the anoxic zone, with neither O2 nor H2S detectable, is characterised by an increase of MnOx concentration and a pH minimum. This MnOx maximum can be explained by the upward diffusing Mn2+ that is biologically oxidised by the downward diffusing O2, though at low, undetectable, concentration (Clement et al., 2009). Additionally, Mn2+ could be oxidised by the nitrite or the nitrate (not measured) diffusing downward (thermodynamically favourable, Luther, 2010). Compared to the ILO zone, the anoxic MnOx maximum corresponds to an increase of 0.2 µM in 2017, while it is much more marked in 2018 with an increase of 1.6 µM (Fig. 3). This difference could come from a faster Mn2+ oxidation produced by the steeper oxygen gradient above the anoxic zone and by the thinness of the anoxic layer in 2018. Finally, in the deeper sulfidic layer, the MnOx disappearance corresponds to the Mn2+ increases (Fig. 3) according to the reduction of settling MnOx by H2S. The concentration increase of the manganese pool and of the iron pool with depth in the anoxic and sulfidic layers probably results from important sedimentary efflux.

4.1 Validity of the reaction-driven approximation

The simultaneous and high-resolution sampling of multiple carbonate parameters and redox species gives us the rare opportunity to investigate in detail the interaction between carbonate species and redox-sensitive elements. In the Chesapeake Bay, main changes of DIC and TA can be explained by mixing between upstream and oceanic endmembers (Appendix, Fig. A1). This “endmember-driven” interpretation leads to the calculation of an excess of DIC and TA, DICex and TAex, relative to the mixing line, as shown in Fig. 4a. At station 858, the steep gradient observed, for example the pH and pCO2 gradients in the PP zone; the O2 and NO2- gradients above the anoxic zone; and the Mn, Fe and H2S gradients at depth, suggest that ongoing in situ processes control the changes of concentrations and dominate the time-dependent endmember variability or the mixing with an unknown third endmember. Additionally, the TAex versus DICex plot (Fig. 4a) shows steep changes of direction, such as in the case of straight lines between local endmembers maintained by ongoing reactions, while a preponderance of mixing would produce more progressive changes. Finally, at depth, Fig. 4a shows a similar slope for both years studied (ΔTAex/ΔDICex=2.4) rather than a similar TAex and DICex concentration, which reinforces the validity of the reaction-driven approximation.

Figure 4Description of the TAex/DICex/AOU or H2S system. Panels (a), (b) and (d) show samples measured along with interpretative slopes discussed in the main text. Panel (c) shows the theoretical slope of TAex/DICex from combination of reactions presented in Table 1.


Assuming (1) that mixing is efficiently described by vertical turbulent diffusion mixing; (2) that the measured concentrations correspond to a steady state – no changes observed over the 1 week sampling; (3) that no additional endmember contributes significantly to the excess calculation, in particular that the samples are isolated from atmospheric exchanges; and (4) that lateral mixing is negligible, which is equivalent to the lateral invariance of the system – as in the stratified water column of station 858, the reaction-driven approximation (Sect. 2.3.1) permits interpretation of the concentration changes as a linear combination of the stoichiometry of several chemical reactions (Eq. 5). The weighted coefficients of each reaction are equal to the rate of each reaction relative to the sum of the rate of all occurring reactions, vi/v (Eq. 4). Accordingly, this interpretation does not identify reactions with minor impact on the carbonate cycle (Eq. 7) or reactions that cancelled each other later during the journey. For example, ΔTAex introduced by PP is frequently cancelled by a similar amount or even higher ΔTAex by AR. The starting point corresponds to the upstream (S=1.5) endmember at DICex=0 and TAex=0. Figure 4a shows a slight DICex enrichment at TAex=0, which reflects a bias from the endmember calculation.

4.2 Identification of major reactions in oxygenated water column using ΔTAex/ΔDICex/ΔAOU signature

In oxygenated water, the ΔTAex/ΔDICex analysis can be supported and strengthened by looking also at the ΔAOU/ΔDICex ratio (Fig. 4b) and establishing a ΔTAex/ΔDICex/ΔAOU signature for each water layer. Then, a linear combination of up to three reactions can be fit to the observed ΔTAex/ΔDICex/ΔAOU signature. The main candidate reactions to be combined are aerobic respiration (AR), primary production (whose overall mass balance equation is here summarised as negative AR), carbonate dissolution (CD) and nitrification of ammonium (Nit) as presented in Table 1 and in Fig. 4c. In 2017, the whole oxygenated zone, including the PP zone and the ILO zone, is characterised by a ΔTAex/ΔDICex/ΔAOU=0.1/1/1 (Fig. 4a and b), which corresponds to the occurrence of only net aerobic respiration (AR) (theoretical values are 0.15/1/1, see Table 1 and Fig. 4c). Note (a) that “net aerobic respiration” indicates that primary production is possible at a significant rate but slower than AR and (b) that this ΔTAex/ΔDICex/ΔAOU signature indicates weak nitrification following respiration. In case of full nitrification, the theoretical slopes should be ΔTAex/ΔDICex/ΔAOU=-0.167/1/1.31, as proposed by Zeebe and Wolf-Gladrow (2001). A combination of 1 AR and 0.025 Nit (nitrification of 16 % of the produced NH3) improves the modelled value to 0.1/1/1.05. The relatively slow nitrification can be explained by slow kinetics of NH4+ oxidation, with a half-life time estimated between a few days in estuaries (Horrigan et al., 1990) to multiple years in coastal environments (Heiss and Fulweiler, 2016). As a comparison, other NH4+ fates, such as adsorption, lead to an ammonium half-life time of about a few minutes (Alshameri et al., 2018) to a few hours (Raaphorst and Malschaert, 1996) depending on the concentration of fine particles. Additionally, algae are known to use NH4+ as an N source (Raven et al., 1992), and NH4+ can be directly assimilated by heterotrophic organisms.

In 2018, fresh water masses brought by the exceptional flood modified the carbonate system equilibrium. First, a low-salinity layer with pCO2 at 1540 µatm overlays the primary production layer (Fig. 2), preventing the uptake of atmospheric CO2 by primary production as was observed in 2016 (Chen et al., 2020). Just below the surface layer, in the PP zone, the lock down of atmospheric exchange by the low-salinity layer produces supersaturation of trapped O2 (Fig. 2). In Fig. 4a and b, this process translates into a vertical distribution at DICex=40µM with ΔTAex/ΔDICex/ΔAOU=1.37/0/-1. This original signature can be modelled by the combination of simultaneous carbonate dissolution (CD), the water column being undersaturated, and PP, no important turbidity was visible as modelled by Cerco et al. (2013), in equal proportion (second line in Table 2); the carbonate dissolution buffers the DIC consumption by the PP. Note that the ratio between ΔTAex/ΔAOU implicates an important nitrate assimilation superior or equal to the amount of N required for the PP, as modelised by negative nitrification in Table 2.

Table 2Linear combination of reactions from Table 1 that fit the observations (see text for details, H2O molecules are omitted).

Download Print Version

In 2018, below the PP zone down to DICex=240µM, the beginning of the ILO zone, the TAex increases significantly with a ΔTAex/ΔDICex/ΔAOU=0.67/1/0.67 (Fig. 4a and b), incompatible with AR. To explain this signature, a contribution of CD superimposed on AR seems most likely. A linear combination fitting leads to 0.5 CD for 1 AR (third line in Table 2) and results in ΔTAex/ΔDICex/ΔAOU=0.77/1/0.67. An explanation of its occurrence solely in 2018 could be the increase in carbonate-rich suspended material at high flow conditions (Su et al., 2021). Excess of Ca2+ compared to the mixing line with oceanic end member (not shown) indicates that up to 200 µM of Ca2+ is produced in the oxygenated layer. Deeper, in the ILO zone, ΔTAex/ΔDICex/ΔAOU=0.2/1/1.25 (Fig. 4a and b), which results mainly from AR (0.15/1/1) with the possible addition of CD and Nit, the exact signature being fitted for 0.54 CD and 0.46 Nit for 1 AR (fourth line in Table 2), as a continuum of AR and CD relative rates from the overlaying layer. This important nitrification is also in good agreement with the nitrification observed in the PP zone during this campaign, the lack of nitrite build-up in the ILO zone and the relatively high oxygen concentration (at 105 µM) in the ILO zone in 2018 which is able to sustain nitrification. Overall, our results show that the higher river flow of 2018 increases carbonate dissolution for the top 5–10 m water depth that is superimposed on primary production or aerobic respiration. Additionally, in 2017 the ΔTAex/ΔDICex/ΔAOU system indicates weak nitrification, while in 2018 significant nitrification in the ILO and PP zones is suggested by the reaction-driven approximation. The role of nitrification in explaining TAex depletion is only hypothetical since no direct measurements of NH4+ and NO3- were performed. In particular, TAex depletion is particularly intense during high flow, high suspended particles season and could be produced by NH4+ adsorption to the particles rather than by nitrification.

4.3 Identification of major reactions in the anoxic water column using the ΔTAex/ΔDICex/ΔH2S signature

In the absence of oxygen, the ΔTAex/ΔDICex/ΔH2S will be used for rates calculation. In the anoxic zone, the signature is similar in summer 2017 and 2018 at ΔTAex/ΔDICex/ΔH2S=2.4/1/0. The reaction-driven approximation of this signature is more difficult than in oxidised water because more reactions are known to occur simultaneously in the absence of oxygen. However, to fit with the reaction-driven approximation, it is not necessary to describe each reaction step but only the overall changes concerning the journey of a water mass over different redox conditions resulting in a combined result or reaction. This approach has recently been proposed for FeS burial by Hiscock and Millero (2006), Rassmann et al. (2020) or Su et al. (2020b). A scenario combining sulfate reduction (SR) is particularly attractive since SR represents the main carbon remineralisation pathway in absence of oxygen. However, a combination of SR with CD would result in a ΔTAex/ΔDICex between 1.15 and 2 (see Table 1 or Fig. 4c) and fails to reach the ΔTAex/ΔDICex of 2.4. Moreover, SR alone underestimates the importance of the H2S oxidation pathway that can consume all the alkalinity produced during SR. For example, SR follow by oxygenated oxidation results in a ΔTAex/ΔDICex/ΔAOU signature equal to AR only. In the Chesapeake Bay, H2S oxidation is critical since no H2S is measurable in the anoxic zone, while the gradient at the sediment–water interface indicates high H2S sedimentary efflux (Fig. 2).

Generalising these observations, recent efforts to build an alkalinity budget on the global scale (Hu and Cai, 2011; Middelburg et al., 2020) highlight that the alkalinity produced by anaerobic respiration corresponds to the uncharged species produced, mostly in solid or gaseous phases. Indeed, the alkalinity changes produced during a natural reaction equal the “charge transfer” from species having some charge at pH=4.5, such as NO3- and SO42-, to a species that would lose its charges at pH=4.5, mainly HCO3-, that is not counted in the alkalinity calculation (see Eq. 12). Although correct, this approach tends to neglect the roles of Fe and Mn oxides (Middelburg et al., 2020) since their transformation from (oxyhydr)oxides into sulfur or carbonate solid species does not involve any charge transfer. When looking in detail at these processes, the metal oxides are critical since they are the main H2S oxidation pathway that does not regenerate H2SO4 but rather produces S0, Sn2- or FeS (Findlay et al., 2014; Avetisyan et al., 2021), which limits alkalinity consumption.

To build a pool of candidate reactions for the fitting, first, dissolved species at too low a concentration (e.g. Mn2+aq, Fe2+aq) to be a net reagent to affect the carbon cycle at steady state are not taken into account. These species are usually recycled rapidly and hold the role of catalyser or electron shuttle between other redox species and did not reach 10 µM during the campaigns (Fig. 3). Second, many minerals are expected to be at low concentration or thermodynamically not favoured, and their associated reactions are neglected (e.g. iron phosphate, ferrous or manganous oxide, metal sulfur clusters, MnS, FeCO3, adsorption processes, reverse weathering). Therefore, only aqueous species with important stock concentrations (that can exceed 0.1 mM in anoxic water) are taken into account, i.e. SO42-, Ca2+, H2S and NH4+ together with gaseous (N2, CO2) and main solid phases (FeS2, FeS, S0, MnCO3, FeOOH, MnO2). Third, the combination of many carbon remineralisation reactions with re-oxidation reactions or the net result produces a net chemical equation equal to another remineralisation reaction. As an example, the chemical equation of SR followed by H2S oxidation with oxygen is equal to the equation of aerobic respiration: the proposed model confounds both pathways because the resulting chemical changes are similar.

Table 1 lists the resulting combined reactions, and the calculated ΔTAex/ΔDICex slopes are represented in Fig. 4c. For the anoxic zone, the nitrate respiration can be associated with N2 production (NR-N2) or with NH3 production (NR-NH3), and the manganese oxide respiration produces carbonate precipitation (MnR-MnC). For the sulfidic zone, SR can occur alone, producing a build-up of H2S. However, at a certain point, H2S gets significantly oxidised either by MnO2 which produces MnCO3 and S0 (SR-SMnC) or by FeOOH, producing FeS and S0 (SR-SFeS) and ultimately FeS2 (SR-FeS2). Direct respiration of FeOOH is also taken into account, but as the only final Fe product in the model is FeS or FeS2, it has to be accompanied by some SR (FeSR-FeS).

Figure 4c demonstrates that the slope of ΔTAex/ΔDICex=2.4 can be obtained for any reaction in combination with MnR-MnC. Combinations without MnR-MnC, however, lead to a negative SR, whose overall equation could be interpreted as a possible small participation of anoxygenic phototrophic (purple) bacteria (Findlay et al., 2015, 2017), but are not considered further, as the amount of ΔTAex involved would be tiny. Therefore, in the absence of nitrate, oxygen and H2S, only a combination of MnR-MnC with SR-SMnC (producing S0, fifth line in Table 2), SR-FeS (producing FeS, sixth line in Table 2) or CD (releasing Ca2+, seventh line in Table 2) gives the particularly high ΔTAex/ΔDICex of 2.4. S0 was not measured during our campaign, but it has been previously reported in this water column (Findlay et al., 2014), and the S0 produced by SR-SMnC can react with FeS to form FeS2. All the three identified combinations require a critical role of MnO2. Since Fig. 3 does not indicate any clear reaction for iron while the steep gradient of MnOx in proximity to the H2S rich layer suggests reaction between MnOx and H2S, the combination with SR-SMnC is the most likely (fifth line in Table 2). Deeper, the vertical gradient of sulfide suggests that part of the H2S came by diffusion from the sediment's porewater (Fig. 2). The assumptions required for the reaction-driven approximation are still valid as soon as steady state is maintained by ongoing reactions, even if one of the local endmembers has not been sampled, since it is probably located in the sediment. In the presence of sulfide, the ΔTAex/ΔDICex/ΔH2S signature is 2.4/1/1.2 in 2017 and 2.4/1/3.2 in 2018 (Fig. 4a and d) and can be explained by the same combination of reactions without complete oxidation of H2S from SR to take into account the build-up of H2S (Table 2).

4.4 Comparison with other studies

The high observed ratio of ΔTAex/ΔDICex=2.4 seems very specific to the Chesapeake Bay. Moreover, the reaction-driven interpretation can be applied to other published datasets for which the ΔTA/ΔDIC/ΔH2S system can be calculated (Table 3). In the water column, most of the available datasets are not suitable for the reaction-driven approximation, since either they focus on surface water where DIC and TA are strongly impacted by atmospheric exchange or the water masses change too fast to consider that reactions dominate over water mixing. However, porewater measured in the Gulf of Mexico has ΔTAex/ΔDICex/ΔH2S=1.15/1/0.53 (Hu et al., 2010), as expected when sulfate reduction is associated with H2S accumulation (SR reaction; ΔTAex/ΔDICex/ΔH2S=1.15/1/0.5; Table 1). Similarly, Hiscock and Millero (2006) report ΔTAex/ΔDICex/ΔH2S=1.3/1/0.5 in the western Black Sea close to the SR signature. In the Baltic Sea sediment, ΔTAex/ΔDICex/ΔH2S=1.3/1/0.07 was reported (Lukawska-Matuszewska, 2016), which is close to the expected signature in the case of important H2S consumption by Fe oxides and consequent precipitation as pyrite (SR-FeS2 reaction; ΔTAex/ΔDICex/ΔH2S=1.2/1/0). In Rhône river prodelta sediments, the reported ΔTAex/ΔDICex/ΔH2S is 1/1/0, which can be related to the 1/1/0 signature of FeSR-FeS (Table 1) that is expected in iron-rich sediments, with high sedimentation rates preventing FeS2 formation in the pore water (Rassmann et al., 2020). In permeable, carbonate-rich sediments, the reported signature of ΔTAex/ΔDICex/ΔH2S in Hawaii sands is 0.86/1/ND (Drupp et al., 2016). The lack of salinity and oxygen datasets prevents further model fits, but the ΔTAex/ΔDICex is below 2.4. For oxygen-depleted data from the whole Chesapeake Bay described in Su et al. (2020a), the signatures are ΔTAex/ΔDICex/ΔAOU=0.2/1/1 in the presence of oxygen (typical signature of AR) and 0.8/1/ND in the absence of oxygen that could correspond to NR-N2 or a combination of CD and Nit. Overall, this bibliographic survey highlights the effectiveness of the reaction-driven approximation to identify preponderant reactions controlling the carbon cycle and puts into perspective the originality of the ΔTAex/ΔDICex/ΔAOU signature of 2.4/1/0 observed in the Chesapeake Bay.

Table 3Overview of the ΔTA/ΔDIC/ΔH2S signature observed in different environments.

 Lack of data on oxygen concentration prevent any interpretation.  Important air–water exchange prevents any reaction-driven approximation.

Download Print Version | Download XLSX

4.5 Local budget

While the reaction-driven approximation indicates a dominant role of the SR-SMnC reaction, this possibility needs to be validated by looking at the saturation state of rhodochrosite (the main MnCO3 mineral) and looking at the mass budget between MnOx consumed and TAex produced. The rhodochrosite saturation (Luo and Millero, 2003) is always below 0.3 in our samples, which stands against the occurrence of in situ SR-SMnC reaction. When inspecting the mass budget, the 88 to 155 µM of MnO2 required to produce the 100 µM TAex increase observed (Fig. 4a) is 1 order of magnitude higher than the observed MnOx or Mn2+ concentration (Fig. 3). This mass budget discrepancy cannot be solved by invoking suspended material, since the 88 µM of MnO2 would require a suspended material concentration of about 4.4 g L−1 (assuming an average concentration of 20 µmol g−1 of Mn), which is again 1 or 2 orders of magnitude higher than the 0.01–0.1 g L−1 usually found in the Chesapeake Bay (Cerco et al., 2013). However, a fast settling rate could satisfy and explain the discrepancy between water and solid concentration. But another process dephasing aqueous from solid reaction products is also possible at station 858, since the dissolved phase could have moved up, rather than the particles settling down. In this case, the SR-SMnC reaction was not happening only in the water column of the Chesapeake Bay, and part of the TAex and DICex pool could have been produced in the sediment during the previous year (e.g. Aller, 2014), then diffused out of the sediment simultaneously with other reduced elements as the summer begins. Indeed, many previous studies at station 858 (e.g. Sholkovitz et al., 1992) explained the seasonality of anoxia with an upward move of the redox front from the sediments to bottom waters during the start of summer. Important sedimentary efflux of H2S, Feaq and Mnaq were still visible during both of our August campaigns. Therefore, the 100 µM TA pool does not fit with the ambient Mn2+ or MnO2 in the water column but rather with the MnCO3 deposited in the sediment. The sedimentary solid Mn stock of the Chesapeake Bay is particularly important, up to 70 µmol g−1 at station 858 (Sinex and Helz, 1981) compared to an average value of 15 µmol g−1 for the upper continental crust (Rudnick and Gao, 2003). Indeed semi-enclosed basins are known to concentrate manganese at the deeper sediment (Thamdrup and Dalsgaard, 2000; Lenstra et al., 2020). Recent investigations at a close station (ET 5.1; 3848.36 N; 7554.66 W) in the Chesapeake Bay (Lenstra et al., 2021) report that about 60 % of the surface sedimentary Mn pool is MnOx (acid ascorbic extractable) and 25 % is Mn carbonate (1 M HCl extraction). Assuming a porosity of 0.8 and a bulk solid density of 2.6, the sedimentary pool corresponds to 35 mM of manganese, which largely exceeds the 88 µM required to produce the 100 µM TAex increase. Therefore, the Chesapeake Bay sediment is particularly rich in manganese and could host important SR-SMnC reactions in the surface pore water, whose soluble products diffuse up to the water column during summer and could bear with them the high ΔTAex/ΔDICex signature observed.

The assumptions detailed in Sect. 4.1 permit the reaction-driven approach to be reconsidered in regard to this sediment efflux scenario. Indeed, the sediment efflux does not need to be considered as an additional endmember; that would violate the third assumption, since its salinity and the pore water concentrations results from the upstream and oceanic endmembers superimposed to geochemical reactions. In Sect. 2.3.1, we point out that Eq. (1) was valid in between each local endmember and that the straight lines on the ΔTAex/ΔDICex plot between them indicate that the local endmembers are maintained in steady state by a chemical reaction with a similar stoichiometry. These results indicate that, even if most of the MnCO3 was produced when the local endmembers were localised in the sediment, (a) their migration does not alter the TAex/DICex signature, and (b) the chemical reaction that produced them is still ongoing at a sufficient rate to maintain a steady state characterised by the steep changes of direction observed in Fig. 4a.

5 Conclusion

The reaction-driven approximation proposed in this study is a powerful interpretative framework that can identify the major reactions controlling the carbonate cycle. In the Chesapeake Bay, similar redox stratification can support varying intensity of carbonate dissolution, absent as in 2017 or important as in 2018. In 2018, the reaction-driven approximation also suggests an important role of the N cycle to consume alkalinity, while in 2017, nitrification was limited to 16 % of the nitrogen mineralised. The summer anoxia observed in the Chesapeake Bay is characterised by an exceptionally high ΔTAex/ΔDICex of 2.4 which has never been reported in anoxic water columns or sediment pore waters. The reaction-driven approximation suggests that it comes from sulfate reduction with almost complete hydrogen sulfide oxidation by MnOx followed by MnCO3 precipitation. This interpretation is supported by the important manganese dynamics observed in this and previous papers (Oldham et al., 2017a; Sholkovitz et al., 1992; Trouwborst et al., 2006). However, most of the reactions would have occurred in the sediment from which the components of the ΔTAex/ΔDICex of 2.4 diffused into the bottom water with the redox front during the set-up of summer water column anoxia. Our study demonstrates that the manganese cycle can have a strong impact on alkalinity, as it prevents H2S oxidation to SO42- and favours sulfur burial.

Appendix A

A1 Demonstration of Eq. (1)

First, Eq. (A1) describes the behaviour of a solute, in the case of turbulent diffusion mixing, with Ds the effective diffusion coefficient, superimposed on a chemical reaction of rate v as described in Eq. (A2).


Assuming steady state, i.e. dC/dt=0, and applying Eq. (1) to two different species sharing the same reactions at rate v on each point and αD≠0, we can express the rate of the reaction as

(A3) - v α C D s = d 2 C d x 2 = α C α D d 2 D d x 2 .

Equation (A3) is also true for any linear combination of a solute with a conservative element such as the salinity, S, because αS=0, hence (with k any real)


By integration of Eq. (A3) with respect to x on a portion of space where αC/αD is constant, it becomes Eq. (1).

A2 Validity of Eq. (2) when the two endmembers have similar concentrations

From Eq. (A3), and assuming that the rate v can be described by a polynomial function of x, with ai, any real coefficient, we obtain


That can be solved analytically and gives the following equation:


Here, ΔC is the difference of concentration between the two endmembers. With similar reasoning on D and with ΔC=ΔD=0, we find Eq. (1) with G=0.

A3 Additional figures

Figure A1Variation of total alkalinity (TA) during oceanic and river mixing. Dashed lines represent the theoretical DIC and TA if only mixing occurs.


Figure A2Superimposed carbonate and redox chemistry profiles over 11 casts against depth.


Figure A3Examples of redox and metal chemistry profiles. Note that the Fe(III)s and Fe(III)aq scales are zoomed in by 3 and the O2 scale is zoomed out by 100.


Data availability

The data used in this paper are available on request to the corresponding author.

Author contributions

ATdC, ERE, JN, BMT and JS performed the data analysis. ATdC and JS processed the data. ATdC, GWL, JS and WJC interpreted the results. GWL, BMT and WJC got the funding. ATdC wrote the paper with contributions from all authors.

Competing interests

The contact author has declared that none of the authors has any competing interests.


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


We gratefully acknowledge the support of the captain and crew of R/V Hugh R. Sharp.

Financial support

This research has been supported by the Division of Ocean Sciences (grant nos. 1558738, 1558692 and 1756815).

Review statement

This paper was edited by Jack Middelburg and reviewed by three anonymous referees.


Abril, G., Etcheber, H., Delille, B., Frankignoulle, M., and Borges, A. V.: Carbonate dissolution in the turbid and eutrophic Loire estuary, Mar. Ecol. Prog. Ser., 259, 129–138, 2003. 

Abril, G., Commarieu, M.-V., Maro, D., Fontugne, M., Guérin, F., and Etcheber, H.: A massive dissolved inorganic carbon release at spring tide in a highly turbid estuary: DIC at spring tide in a estuary, Geophys. Res. Lett., 31, L09316,, 2004. 

Aller, R. C.: 8.11 – Sedimentary Diagenesis, Depositional Environments, and Benthic Fluxes, in: Treatise on Geochemistry, 2nd edn., edited by: Turekian, H. D. H. K., Elsevier, Oxford, 293–334, ISBN 9780128182345, 2014. 

Alshameri, A., He, H., Zhu, J., Xi, Y., Zhu, R., Ma, L., and Tao, Q.: Adsorption of ammonium by different natural clay minerals: Characterization, kinetics and adsorption isotherms, Appl. Clay Sci., 159, 83–93,, 2018. 

Avetisyan, K., Zweig, I., Luther, G. W., and Kamyshny, A.: Kinetics and mechanism of polysulfides and elemental sulfur formation by a reaction between hydrogen sulfide and δ-MnO2, Geochim. Cosmochim. Ac., 313, 21–37,, 2021. 

Borges, A. V., Schiettecatte, L.-S., Abril, G., Delille, B., and Gazeau, F.: Carbon dioxide in European coastal waters, Estuar. Coast. Shelf S., 70, 375–387,, 2006. 

Borges, A. V., Abril, G., and Bouillon, S.: Carbon dynamics and CO2 and CH4 outgassing in the Mekong delta, Biogeosciences, 15, 1093–1114,, 2018. 

Boudreau, B. P., Middelburg, J. J., and Luo, Y.: The role of calcification in carbonate compensation, Nat. Geosci., 11, 894–900,, 2018. 

Breitburg, D., Levin, L. A., Oschlies, A., Grégoire, M., Chavez, F. P., Conley, D. J., Garçon, V., Gilbert, D., Gutiérrez, D., Isensee, K., Jacinto, G. S., Limburg, K. E., Montes, I., Naqvi, S. W. A., Pitcher, G. C., Rabalais, N. N., Roman, M. R., Rose, K. A., Seibel, B. A., Telszewski, M., Yasuhara, M., and Zhang, J.: Declining oxygen in the global ocean and coastal waters, Science, 359, eaam7240,, 2018. 

Cai, W.-J. and Wang, Y.: The chemistry, fluxes, and sources of carbon dioxide in the estuarine waters of the Satilla and Altamaha Rivers, Georgia, Limnol. Oceanogr., 43, 657–668,, 1998. 

Cai, W.-J., Hu, X., Huang, W.-J., Jiang, L.-Q., Wang, Y., Peng, T.-H., and Zhang, X.: Alkalinity distribution in the western North Atlantic Ocean margins, J. Geophys. Res., 115, C08014,, 2010. 

Cai, W.-J., Huang, W.-J., Luther, G. W., Pierrot, D., Li, M., Testa, J., Xue, M., Joesoef, A., Mann, R., Brodeur, J., Xu, Y.-Y., Chen, B., Hussain, N., Waldbusser, G. G., Cornwell, J., and Kemp, W. M.: Redox reactions and weak buffering capacity lead to acidification in the Chesapeake Bay, Nat. Commun., 8, 369,, 2017. 

Cai, W.-J., Xu, Y.-Y., Feely, R. A., Wanninkhof, R., Jönsson, B., Alin, S. R., Barbero, L., Cross, J. N., Azetsu-Scott, K., Fassbender, A. J., Carter, B. R., Jiang, L.-Q., Pepin, P., Chen, B., Hussain, N., Reimer, J. J., Xue, L., Salisbury, J. E., Hernández-Ayón, J. M., Langdon, C., Li, Q., Sutton, A. J., Chen, C.-T. A., and Gledhill, D. K.: Controls on surface water carbonate chemistry along North American ocean margins, Nat. Commun., 11, 2691,, 2020. 

Carstensen, J., Chierici, M., Gustafsson, B. G., and Gustafsson, E.: Long-Term and Seasonal Trends in Estuarine and Coastal Carbonate Systems, Global Biogeochem. Cy., 32, 497–513,, 2018. 

Cerco, C. F., Kim, S.-C., and Noel, M. R.: Management modeling of suspended solids in the Chesapeake Bay, USA, Estuar. Coast. Shelf S., 116, 87–98,, 2013. 

Chen, B., Cai, W.-J., Brodeur, J. R., Hussain, N., Testa, J. M., Ni, W., and Li, Q.: Seasonal and spatial variability in surface pCO2 and air–water CO2 flux in the Chesapeake Bay, Limnol. Oceanogr., 65, 3046–3065,, 2020. 

Clement, B. G., Luther, G. W., and Tebo, B. M.: Rapid, oxygen-dependent microbial Mn(II) oxidation kinetics at sub-micromolar oxygen concentrations in the Black Sea suboxic zone, Geochim. Cosmochim. Ac., 73, 1878–1889,, 2009. 

Cotovicz Jr., L. C., Libardoni, B. G., Brandini, N., Knoppers, B. A., and Abril, G.: Comparisons between real-Time PCO2 measurements with indirect estimates in two contrasting Brazilian estuaries: The eutrophic guanabara bay (RJ) and the oligotrophic sao francisco River estuary (AL), Quim. Nova, 39, 1206–1214,, 2016. 

Dickson, A. G.: An exact definition of total alkalinity and a procedure for the estimation of alkalinity and total inorganic carbon from titration data, Deep-Sea Res., 28, 609–623,, 1981. 

Drupp, P. S., De Carlo, E. H., and Mackenzie, F. T.: Porewater CO2–carbonic acid system chemistry in permeable carbonate reef sands, Mar. Chem., 185, 48–64,, 2016. 

Du, J. and Shen, J.: Water residence time in Chesapeake Bay for 1980–2012, J. Marine Syst., 164, 101–111,, 2016. 

Findlay, A. J., Gartman, A., MacDonald, D. J., Hanson, T. E., Shaw, T. J., and Luther, G. W.: Distribution and size fractionation of elemental sulfur in aqueous environments: The Chesapeake Bay and Mid-Atlantic Ridge, Geochim. Cosmochim. Ac., 142, 334–348, 2014. 

Findlay, A. J., Bennett, A. J., Hanson, T. E., and Luther, G. W.: Light-Dependent Sulfide Oxidation in the Anoxic Zone of the Chesapeake Bay Can Be Explained by Small Populations of Phototrophic Bacteria, Appl. Environ. Microb., 81, 7560–7569,, 2015. 

Findlay, A. J., Di Toro, D. M., and Luther, G. W.: A model of phototrophic sulfide oxidation in a stratified estuary, Limnol. Oceanogr., 62, 1853–1867,, 2017. 

Friedlingstein, P., Jones, M. W., O'Sullivan, M., Andrew, R. M., Hauck, J., Peters, G. P., Peters, W., Pongratz, J., Sitch, S., Le Quéré, C., Bakker, D. C. E., Canadell, J. G., Ciais, P., Jackson, R. B., Anthoni, P., Barbero, L., Bastos, A., Bastrikov, V., Becker, M., Bopp, L., Buitenhuis, E., Chandra, N., Chevallier, F., Chini, L. P., Currie, K. I., Feely, R. A., Gehlen, M., Gilfillan, D., Gkritzalis, T., Goll, D. S., Gruber, N., Gutekunst, S., Harris, I., Haverd, V., Houghton, R. A., Hurtt, G., Ilyina, T., Jain, A. K., Joetzjer, E., Kaplan, J. O., Kato, E., Klein Goldewijk, K., Korsbakken, J. I., Landschützer, P., Lauvset, S. K., Lefèvre, N., Lenton, A., Lienert, S., Lombardozzi, D., Marland, G., McGuire, P. C., Melton, J. R., Metzl, N., Munro, D. R., Nabel, J. E. M. S., Nakaoka, S.-I., Neill, C., Omar, A. M., Ono, T., Peregon, A., Pierrot, D., Poulter, B., Rehder, G., Resplandy, L., Robertson, E., Rödenbeck, C., Séférian, R., Schwinger, J., Smith, N., Tans, P. P., Tian, H., Tilbrook, B., Tubiello, F. N., van der Werf, G. R., Wiltshire, A. J., and Zaehle, S.: Global Carbon Budget 2019, Earth Syst. Sci. Data, 11, 1783–1838,, 2019. 

Grasshoff, K., Kremling, K., and Ehrhardt, M.: Methods of seawater analysis, Third, Completely Revised and Extended Edition, John Wiley & Sons, 1999. 

Gupta, G. V. M., Jyothibabu, R., Ramu, C. V., Reddy, A. Y., Balachandran, K. K., Sudheesh, V., Kumar, S., Chari, N. V. H. K., Bepari, K. F., Marathe, P. H., Reddy, B. B., and Vijayan, A. K.: The world's largest coastal deoxygenation zone is not anthropogenically driven, Environ. Res. Lett., 16, 054009,, 2021. 

Gustafsson, E., Hagens, M., Sun, X., Reed, D. C., Humborg, C., Slomp, C. P., and Gustafsson, B. G.: Sedimentary alkalinity generation and long-term alkalinity development in the Baltic Sea, Biogeosciences, 16, 437–456,, 2019. 

Hagy, J. D., Boynton, W. R., Keefe, C. W., and Wood, K. V.: Hypoxia in Chesapeake Bay, 1950–2001: Long-Term Change in Relation to Nutrient Loading and River Flow, Estuaries, 27, 634–658, 2004. 

Heiss, E. M. and Fulweiler, R. W.: Coastal water column ammonium and nitrite oxidation are decoupled in summer, Estuar. Coast. Shelf S., 178, 110–119,, 2016. 

Hiscock, W. T. and Millero, F. J.: Alkalinity of the anoxic waters in the Western Black Sea, Deep-Sea Res. Pt. II, 53, 1787–1801,, 2006. 

Horrigan, S. G., Montoya, J. P., Nevins, J. L., McCarthy, J. J., Ducklow, H., Goericke, R., and Malone, T.: Nitrogenous nutrient transformations in the spring and fall in the Chesapeake Bay, Estuar. Coast. Shelf S., 30, 369–391,, 1990. 

Hu, X. and Cai, W.-J.: An assessment of ocean margin anaerobic processes on oceanic alkalinity budget, Global Biogeochem. Cy., 25, GB3003,, 2011. 

Hu, X., Cai, W.-J., Wang, Y., Luo, S., and Guo, X.: Pore-water geochemistry of two contrasting brine-charged seep sites in the northern Gulf of Mexico continental slope, Mar. Chem., 118, 99–107,, 2010. 

Huang, W.-J., Wang, Y., and Cai, W.-J.: Assessment of sample storage techniques for total alkalinity and dissolved inorganic carbon in seawater, Limnol. Oceanogr.-Meth., 10, 711–717,, 2012. 

Hudson, J. M., MacDonald, D. J., Estes, E. R., and Luther, G. W.: A durable and inexpensive pump profiler to monitor stratified water columns with high vertical resolution, Talanta, 199, 415–424,, 2019. 

Hulot, V., Metzger, E., Thibault De Chanvalon, A., Mouret, A., Schmidt, S., Deflandre, B., Rigaud, S., Beneteau, E., Savoye, N., Souchu, P., Le Merrer, Y., and Maillet, G. M.: Impact of an exceptional winter flood on benthic oxygen and nutrient fluxes in a temperate macrotidal estuary: Potential consequences on summer deoxygenation, Front. Mar. Sci., 10, 1083377,, 2023. 

Ishii, H., Koh, H., and Satoh, K.: Spectrophotometric determination of manganese utilizing metal ion substitution in the cadmium -α, β,-γ,δ -tetrakis (4-carboxyphenyl) porphine complex, Anal. Chim. Acta, 136, 347–352, 1982. 

Joesoef, A., Kirchman, D. L., Sommerfield, C. K., and Cai, W.-J.: Seasonal variability of the inorganic carbon system in a large coastal plain estuary, Biogeosciences, 14, 4949–4963,, 2017. 

Jones, M. R., Luther, G. W., Mucci, A., and Tebo, B. M.: Concentrations of reactive Mn(III)-L and MnO2 in estuarine and marine waters determined using spectrophotometry and the leuco base, leucoberbelin blue, Talanta, 200, 91–99,, 2019. 

Lenstra, W. K., Séguret, M. J. M., Behrends, T., Groeneveld, R. K., Hermans, M., Witbaard, R., and Slomp, C. P.: Controls on the shuttling of manganese over the northwestern Black Sea shelf and its fate in the euxinic deep basin, Geochim. Cosmochim. Ac., 273, 177–204,, 2020. 

Lenstra, W. K., Klomp, R., Molema, F., Behrends, T., and Slomp, C. P.: A sequential extraction procedure for particulate manganese and its application to coastal marine sediments, Chem. Geol., 584, 120538,, 2021. 

Lewis, B. L., Glazer, B. T., Montbriand, P. J., Luther, G. W., Nuzzio, D. B., Deering, T., Ma, S., and Theberge, S.: Short-term and interannual variability of redox-sensitive chemical parameters in hypoxic/anoxic bottom waters of the Chesapeake Bay, Mar. Chem., 105, 296–308,, 2007. 

Lohrenz, S. E., Cai, W.-J., Chen, F., Chen, X., and Tuel, M.: Seasonal variability in air-sea fluxes of CO2 in a river-influenced coastal margin, J. Geophys. Res., 115, C10034,, 2010. 

Lukawska-Matuszewska, K.: Contribution of non-carbonate inorganic and organic alkalinity to total measured alkalinity in pore waters in marine sediments (Gulf of Gdansk, S-E Baltic Sea), Mar. Chem., 186, 211–220,, 2016. 

Luo, Y. and Millero, F. J.: Solubility of rhodochrosite (MnCO3) in NaCl solutions, J. Solut. Chem., 32, 405–416, 2003. 

Luther, G. W., Glazer, B. T., Ma, S., Trouwborst, R. E., Moore, T. S., Metzger, E., Kraiya, C., Waite, T. J., Druschel, G., Sundby, B., Taillefert, M., Nuzzio, D. B., Shank, T. M., Lewis, B. L., and Brendel, P. J.: Use of voltammetric solid-state (micro)electrodes for studying biogeochemical processes: Laboratory measurements to real time measurements with an in situ electrochemical analyzer (ISEA), Mar. Chem., 108, 221–235,, 2008. 

Luther, G. W.: The role of one-and two-electron transfer reactions in forming thermodynamically unstable intermediates as barriers in multi-electron redox reactions, Aquat. Geochem., 16, 395–420, 2010. 

Madison, A. S., Tebo, B. M., and Luther, G. W.: Simultaneous determination of soluble manganese(III), manganese(II) and total manganese in natural (pore)waters, Talanta, 84, 374–381,, 2011. 

Meire, L., Soetaert, K. E. R., and Meysman, F. J. R.: Impact of global change on coastal oxygen dynamics and risk of hypoxia, Biogeosciences, 10, 2633–2653,, 2013. 

Meybeck, M.: Global chemical weathering of surficial rocks estimated from river dissolved loads, Am. J. Sci., 287, 401–428, 1987. 

Meybeck, M.: Global occurrence of major elements in rivers, Treatise Geochem., 5, 605, 2003. 

Meybeck, M., Cauwet, G., Dessery, S., Somville, M., Gouleau, D., and Billen, G.: Nutrients (organic C, P, N, Si) in the eutrophic River Loire (France) and its estuary, Estuar. Coast. Shelf S., 27, 595–624,, 1988. 

Middelburg, J. J., Soetaert, K., and Hagens, M.: Ocean Alkalinity, Buffering and Biogeochemical Processes, Rev. Geophys., 58, e2019RG000681,, 2020. 

Oldham, V. E., Jones, M. R., Tebo, B. M., and Luther, G. W.: Oxidative and reductive processes contributing to manganese cycling at oxic-anoxic interfaces, Mar. Chem., 195, 122–128,, 2017a. 

Oldham, V. E., Miller, M. T., Jensen, L. T., and Luther, G. W.: Revisiting Mn and Fe removal in humic rich estuaries, Geochim. Cosmochim. Ac., 209, 267–283,, 2017b. 

Raaphorst, W. V. and Malschaert, J. F. P.: Ammonium adsorption in superficial North Sea sediments, Cont. Shelf Res., 16, 1415–1435,, 1996. 

Rabalais, N. N., Cai, W.-J., Carstensen, J., Conley, D. J., Fry, B., Hu, X., Quiñones-Rivera, Z., Rosenberg, R., Slomp, C. P., Turner, R. E., Voss, M., Wissel, B., and Zhang, J.: Eutrophication-Driven Deoxygenation in the Coastal Ocean, Oceanography, 27, 172–183, 2014. 

Rassmann, J., Eitel, E. M., Lansard, B., Cathalot, C., Brandily, C., Taillefert, M., and Rabouille, C.: Benthic alkalinity and dissolved inorganic carbon fluxes in the Rhône River prodelta generated by decoupled aerobic and anaerobic processes, Biogeosciences, 17, 13–33,, 2020. 

Raven, J. A., Wollenweber, B., and Handley, L. L.: A comparison of ammonium and nitrate as nitrogen sources for photolithotrophs, New Phytol., 121, 19–32,, 1992. 

Rickard, D. and Luther, G. W.: Chemistry of Iron Sulfides, Chem. Rev., 107, 514–562,, 2007. 

Rudnick, R. L. and Gao, S.: Composition of the continental crust, Treatise Geochem., 3, 1–64, 2003. 

Sholkovitz, E. R., Shaw, T. J., and Schneider, D. L.: The geochemistry of rare earth elements in the seasonally anoxic water column and porewaters of Chesapeake Bay, Geochim. Cosmochim. Ac., 56, 3389–3402, 1992. 

Sinex, S. A. and Helz, G. R.: Regional geochemistry of trace elements in Chesapeake Bay sediments, Environ. Geol., 3, 315–323,, 1981. 

Smith, S. V. and Mackenzie, F. T.: The Role of CaCO3 Reactions in the Contemporary Oceanic CO2 Cycle, Aquat. Geochem., 22, 153–175,, 2016. 

Soetaert, K., Hofmann, A. F., Middelburg, J. J., Meysman, F. J. R., and Greenwood, J.: The effect of biogeochemical processes on pH, Mar. Chem., 105, 30–51,, 2007. 

Stookey, L. L.: Ferrozine – a new spectrophotometric reagent for iron, Anal. Chem., 42, 779–781, 1970. 

Su, J., Cai, W.-J., Brodeur, J., Chen, B., Hussain, N., Yao, Y., Ni, C., Testa, J. M., Li, M., Xie, X., Ni, W., Scaboo, K. M., Xu, Y., Cornwell, J., Gurbisz, C., Owens, M. S., Waldbusser, G. G., Dai, M., and Kemp, W. M.: Chesapeake Bay acidification buffered by spatially decoupled carbonate mineral cycling, Nat. Geosci., 13, 441–447,, 2020a. 

Su, J., Cai, W.-J., Brodeur, J., Hussain, N., Chen, B., Testa, J. M., Scaboo, K. M., Jaisi, D. P., Li, Q., Dai, M., and Cornwell, J.: Source partitioning of oxygen-consuming organic matter in the hypoxic zone of the Chesapeake Bay, Limnol. Oceanogr., 65, 1801–1817,, 2020b. 

Su, J., Cai, W.-J., Testa, J. M., Brodeur, J. R., Chen, B., Scaboo, K. M., Li, M., Shen, C., Dolan, M., Xu, Y., Zhang, Y., and Hussain, N.: Supply-controlled calcium carbonate dissolution decouples the seasonal dissolved oxygen and pH minima in Chesapeake Bay, Limnol. Oceanogr., 66, 3796–3810,, 2021. 

Thamdrup, B. and Dalsgaard, T.: The fate of ammonium in anoxic manganese oxide-rich marine sediment, Geochim. Cosmochim. Ac., 64, 4157–4164,, 2000. 

Thibault de Chanvalon, A. and Luther, G. W.: Mn speciation at nanomolar concentrations with a porphyrin competitive ligand and UV–vis measurements, Talanta, 200, 15–21,, 2019. 

Trouwborst, R. E., Clement, B. G., Tebo, B. M., Glazer, B. T., and Luther, G. W.: Soluble Mn(III) in Suboxic Zones, Science, 313, 1955–1957,, 2006. 

Urey, H. C.: On the Early Chemical History of the Earth and the Origin of Life, P. Natl. Acad. Sci. USA, 38, 351–363,, 1952. 

Wolf-Gladrow, D. A., Zeebe, R. E., Klaas, C., Körtzinger, A., and Dickson, A. G.: Total alkalinity: The explicit conservative expression and its application to biogeochemical processes, Mar. Chem., 106, 287–300,, 2007. 

Xue, L. and Cai, W.-J.: Total alkalinity minus dissolved inorganic carbon as a proxy for deciphering ocean acidification mechanisms, Mar. Chem., 222, 103791,, 2020. 

Zeebe, R. E. and Wolf-Gladrow, D.: CO2 in seawater: equilibrium, kinetics, isotopes, Gulf Professional Publishing, ISBN 0 444 50579 2, 2001.  

Zhang, Q., Brady, D. C., Boynton, W. R., and Ball, W. P.: Long-Term Trends of Nutrients and Sediment from the Nontidal Chesapeake Watershed: An Assessment of Progress by River and Season, J. Am. Water Resour. As., 51, 1534–1555,, 2015. 

Short summary
The intensity of the oceanic trap of CO2 released by anthropogenic activities depends on the alkalinity brought by continental weathering. Between ocean and continent, coastal water and estuaries can limit or favour the alkalinity transfer. This study investigate new interactions between dissolved metals and alkalinity in the oxygen-depleted zone of estuaries.
Final-revised paper