Abstract. Barium is widely used as a proxy for dissolved silicon and particulateorganic carbon fluxes in seawater. However, these proxy applications arelimited by insufficient knowledge of the dissolved distribution of Ba([Ba]). For example, there is significant spatial variability in thebarium–silicon relationship, and ocean chemistry may influence sedimentaryBa preservation. To help address these issues, we developed 4095 models forpredicting [Ba] using Gaussian process regression machine learning. Thesemodels were trained to predict [Ba] from standard oceanographic observationsusing GEOTRACES data from the Arctic, Atlantic, Pacific, and Southernoceans. Trained models were then validated by comparing predictions againstwithheld [Ba] data from the Indian Ocean. We find that a model trained usingdepth, temperature, and salinity, as well as dissolved dioxygen, phosphate,nitrate, and silicate, can accurately predict [Ba] in the Indian Ocean with amean absolute percentage deviation of 6.0 %. We use this model tosimulate [Ba] on a global basis using these same seven predictors in theWorld Ocean Atlas. The resulting [Ba] distribution constrains the Ba budgetof the ocean to 122(±7) × 1012 mol and revealsoceanographically consistent variability in the barium–silicon relationship. We then calculate the saturation state of seawater with respect to barite. This calculation reveals systematic spatial and vertical variations in marine barite saturation and shows that the ocean below 1000 m is at equilibrium with respect tobarite. We describe a number of possible applications for our model outputs, ranging from use in mechanistic biogeochemical models to paleoproxy calibration. Ourapproach demonstrates the utility of machine learning in accurately simulatingthe distributions of tracers in the sea and provides a framework that couldbe extended to other trace elements. Our model, the data used in training and validation, and global outputs are available in Horner and Mete (2023, https://doi.org/10.26008/1912/bco-dmo.885506.2).
more »
« less
A spatially and vertically resolved global grid of dissolved barium concentrations in seawater determined using Gaussian Process Regression machine learning
We present a spatially and vertically resolved global grid of dissolved barium concentrations ([Ba]) in seawater determined using Gaussian Process Regression machine learning. This model was trained using 4,345 quality-controlled GEOTRACES data from the Arctic, Atlantic, Pacific, and Southern Oceans. Model output was validated by assessing the accuracy of [Ba] simulations in the Indian Ocean, noting that none of the Indian Ocean data were seen by the model during training. We identify a model that can accurate predict [Ba] in the Indian Ocean using seven features: depth, temperature, salinity, as well as dissolved dioxygen, phosphate, nitrate, and silicate concentrations. This model achieves a mean absolute percentage error of 6.0 %, which we assume represents the generalization error. This model was used to simulate [Ba] on a global basis using predictor data from the World Ocean Atlas 2018. The global model of [Ba] is on a 1°x 1° grid with 102 depth levels from 0 to 5,500 m. The dissolved [Ba] output was then used to simulate dissolved Ba* (barium-star), which is the difference between 'observed' and [Ba] predicted from co-located [Si]. Lastly, [Ba] data were combined with temperature, salinity, and pressure data from the World Ocean Atlas to calculate the saturation state of seawater with respect to barite. The model reveals that the volume-weighted mean oceanic [Ba] and and saturation state are 89 nmol/kg and 0.82, respectively. These results imply that the total marine Ba inventory is 122(±7) ×10¹² mol and that the ocean below 1,000 m is at barite equilibrium.
more »
« less
- NSF-PAR ID:
- 10462075
- Publisher / Repository:
- Biological and Chemical Oceanography Data Management Office (BCO-DMO)
- Date Published:
- Edition / Version:
- 2
- Subject(s) / Keyword(s):
- ["barium","barite","machine learning"]
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract Processes controlling dissolved barium (dBa) were investigated along the GEOTRACES GA03 North Atlantic and GP16 Eastern Tropical Pacific transects, which traversed similar physical and biogeochemical provinces. Dissolved Ba concentrations are lowest in surface waters (∼35–50 nmol kg−1) and increase to 70–80 and 140–150 nmol kg−1in deep waters of the Atlantic and Pacific transects, respectively. Using water mass mixing models, we estimate conservative mixing that accounts for most of dBa variability in both transects. To examine nonconservative processes, particulate excess Ba (pBaxs) formation and dissolution rates were tracked by normalizing particulate excess230Th activities. Th‐normalized pBaxsfluxes, with barite as the likely phase, have subsurface maxima in the top 1,000 m (∼100–200 μmol m−2 year−1average) in both basins. Barite precipitation depletes dBa within oxygen minimum zones from concentrations predicted by water mass mixing, whereas inputs from continental margins, particle dissolution in the water column, and benthic diffusive flux raise dBa above predications. Average pBaxsburial efficiencies along GA03 and GP16 are ∼37% and 17%–100%, respectively, and do not seem to be predicated on barite saturation indices in the overlying water column. Using published values, we reevaluate the global freshwater dBa river input as 6.6 ± 3.9 Gmol year−1. Estuarine mixing processes may add another 3–13 Gmol year−1. Dissolved Ba inputs from broad shallow continental margins, previously unaccounted for in global marine summaries, are substantial (∼17 Gmol year−1), exceeding terrestrial freshwater inputs. Revising river and shelf dBa inputs may help bring the marine Ba isotope budget more into balance.
-
more » « less
Abstract Recent studies have explored the sensitivity of global ocean model simulations to the treatment of riverine freshwater and the representation of estuarine processes via an estuary box model applied within Community Earth System Model (CESM). This study builds on these efforts by assessing the model skill score relative to a new salinity climatology. The new climatology averages the original observational data of the World Ocean Database directly onto the CESM ocean component tracer grid cells without spatial interpolation, smoothing, or other gap‐filling techniques to mitigate coastal ocean salinity bias present in the World Ocean Atlas. The mean square error for coastal upper ocean salinity relative to climatology is reduced by up to 14%, and the mean square error of near‐surface salinity stratification is reduced by up to 28% near major river mouths in the simulations with improved treatments of river runoff. The improvement in upper ocean bulk salinity is attributed primarily to focusing runoff as point sources thereby avoiding the artificial horizontal spreading of the control run and to applying a locally varying instead of a global constant reference salinity for riverine virtual salt fluxes. The improvements in near‐surface salinity stratification are primarily attributed to adding parameterized estuarine mixing with the estuary box model. Salinity and salinity stratification skill improvements are achieved not just near large rivers but also along the global coast and skill improvements extend far offshore. Despite these improvements, many other sources of model‐climatology mismatch in coastal salinity and stratification remain and merit further attention.
-
Keynote points • Thermal expansion from a warming ocean and land ice melt are the main causes of the accelerating global rise in the mean sea level. • Global warming is also affecting many circulation systems. The Atlantic meridional overturning circulation has already weakened and will most likely continue to do so in the future. The impacts of ocean circulation changes include a regional rise in sea levels, changes in the nutrient distribution and carbon uptake of the ocean and feedbacks with the atmosphere, such as altering the distribution of precipitation. • More than 90 per cent of the heat from global warming is stored in the global ocean. Oceans have exhibited robust warming since the 1950s from the surface to a depth of 2,000 m. The proportion of ocean heat content has more than doubled since the 1990s compared with long-term trends. Ocean warming can be seen in most of the global ocean, with a few regions exhibiting long-term cooling. • The ocean shows a marked pattern of salinity changes in multidecadal observations, with surface and subsurface patterns providing clear evidence of a water cycle amplification over the ocean. That is manifested in enhanced salinities in the near-surface, high-salinity subtropical regions and freshening in the low-salinity regions such as the West Pacific Warm Pool and the poles. • An increase in atmospheric CO2 levels, and a subsequent increase in carbon in the oceans, has changed the chemistry of the oceans to include changes to pH and aragonite saturation. A more carbon-enriched marine environment, especially when coupled with other environmental stressors, has been demonstrated through field studies and experiments to have negative impacts on a wide range of organisms, in particular those that form calcium carbonate shells, and alter biodiversity and ecosystem structure. • Decades of oxygen observations allow for robust trend analyses. Long-term measurements have shown decreases in dissolved oxygen concentrations for most ocean regions and the expansion of oxygen-depleted zones. A temperature-driven solubility decrease is responsible for most near-surface oxygen loss, though oxygen decrease is not limited to the upper ocean and is present throughout the water column in many areas. • Total sea ice extent has been declining rapidly in the Arctic, but trends are insignificant in the Antarctic. In the Arctic, the summer trends are most striking in the Pacific sector of the Arctic Ocean, while, in the Antarctic, the summer trends show increases in the Weddell Sea and decreases in the West Antarctic sector of the Southern Ocean. Variations in sea ice extent result from changes in wind and ocean currents.more » « less
-
Abstract To overcome challenges with observing ocean heat content (OHC) over the entire ocean, we propose a novel approach that exploits the abundance of satellite data, including data from modern satellite geomagnetic surveys such as Swarm. The method considers a novel combination of conventional in situ (temperature and pressure) as well as satellite (altimetry and gravimetry) data with estimates of ocean electrical conductance (depth-integrated conductivity), which can potentially be obtained from magnetic observations (by satellite, land, seafloor, ocean, and airborne magnetometers). To demonstrate the potential benefit of the proposed method, we sample model output of an ocean state estimate to reflect existing observations and train a machine learning algorithm [Generalized Additive Model (GAM)] on these samples. We then calculate OHC everywhere using information potentially derivable from various global satellite coverage—including magnetic observations—to gauge the GAM’s goodness of fit on a global scale. Inclusion of in situ observations of OHC in the upper 2000 m from Argo-like floats and conductance data each reduce the root-mean-square error by an order of magnitude. Retraining the GAM with recent ship-based hydrographic data attains a smaller RMSE in polar oceans than training the GAM only once on all available historical ship-based hydrographic data; the opposite is true elsewhere. The GAM more accurately calculates OHC anomalies throughout the water column than below 2000 m and can detect global OHC anomalies over multiyear time scales, even when considering hypothetical measurement errors. Our method could complement existing methods and its accuracy could be improved through careful ship-based campaign planning. Significance Statement The purpose of this manuscript is to demonstrate the potential for practical implementation of a remote monitoring method for ocean heat content (OHC) anomalies. To do this, we sample data from a reanalysis product primarily because of the dearth of observations below 2000 m depth that can be used for validation and the fact that full-depth-integrated electrical seawater conductivity data products derived from satellite magnetometry are not yet available. We evaluate multiple factors related to the accuracy of OHC anomaly estimation and find that, even with hypothetical measurement errors, our method can be used to monitor OHC anomalies on multiyear time scales.more » « less
