skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Fast Spin‐Up of Geochemical Tracers in Ocean Circulation and Climate Models
Abstract Ocean geochemical tracers such as radiocarbon, protactinium and thorium isotopes, and noble gases are widely used to constrain a range of physical and biogeochemical processes in the ocean. However, their routine simulation in global ocean circulation and climate models is hindered by the computational expense of integrating them to a steady state. Here, a new approach to this long‐standing “spin‐up” problem is introduced to efficiently compute equilibrium distributions of such tracers in seasonally‐forced models. Based on “Anderson Acceleration,” a sequence acceleration technique developed in the 1960s to solve nonlinear integral equations, the new method is entirely “black box” and offers significant speed‐up over conventional direct time integration. Moreover, it requires no preconditioning, ensures tracer conservation and is fully consistent with the numerical time‐stepping scheme of the underlying model. It thus circumvents some of the drawbacks of other schemes such as matrix‐free Newton Krylov that have been proposed to address this problem. An implementation specifically tailored for the batch HPC systems on which ocean and climate models are typically run is described, and the method illustrated by applying it to a variety of geochemical tracer problems. The new method, which provides speed‐ups by over an order of magnitude, should make simulations of such tracers more feasible and enable their inclusion in climate change assessments such as IPCC.  more » « less
Award ID(s):
1924215
PAR ID:
10395288
Author(s) / Creator(s):
 
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Advances in Modeling Earth Systems
Volume:
15
Issue:
2
ISSN:
1942-2466
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; (, Journal of Hydrometeorology)
    null (Ed.)
    Abstract Sampling intervals of precipitation geochemistry measurements are often coarser than those required by fine-scale hydrometeorological models. This study presents a statistical method to temporally downscale geochemical tracer signals in precipitation so that they can be used in high-resolution, tracer-enabled applications. In this method, we separated the deterministic component of the time series and the remaining daily stochastic component, which was approximated by a conditional multivariate Gaussian distribution. Specifically, statistics of the stochastic component could be explained from coarser data using a newly identified power-law decay function, which relates data aggregation intervals to changes in tracer concentration variance and correlations with precipitation amounts. These statistics were used within a copula framework to generate synthetic tracer values from the deterministic and stochastic time series components based on daily precipitation amounts. The method was evaluated at 27 sites located worldwide using daily precipitation isotope ratios, which were aggregated in time to provide low resolution testing datasets with known daily values. At each site, the downscaling method was applied on weekly, biweekly and monthly aggregated series to yield an ensemble of daily tracer realizations. Daily tracer concentrations downscaled from a biweekly series had average (+/- standard deviation) absolute errors of 1.69‰ (1.61‰) for δ 2 H and 0.23‰ (0.24‰) for δ 18 O relative to observations. The results suggest coarsely sampled precipitation tracers can be accurately downscaled to daily values. This method may be extended to other geochemical tracers in order to generate downscaled datasets needed to drive complex, fine-scale models of hydrometeorological processes. 
    more » « less
  2. ; (, Journal of Geophysical Research: Oceans)
    Abstract Marine net community production (NCP), a metric of ecosystem functionality, is often estimated as the residual term in a mass balance equation that aims to describe upper ocean variations in the time series of a chemical tracer. The advent of biogeochemical (BGC) Argo profiling floats equipped with nitrate, pH, and oxygen sensors has enabled such NCP estimation across vast ocean regions. Floats typically drift at 1,000 m depth between profiling from ∼2,000 m to the surface every 10 days, resulting in quasi‐Lagrangian time series that can reflect different upper ocean water masses over time. However, limited information about real‐time horizontal tracer gradients often leads to lateral processes being omitted during tracer budget closure, which can bias the residual‐term NCP estimates. To determine the potential magnitude of such biases, we developed a method to quantify and adjust for the impact of lateral float movement across horizontal tracer gradients using dissolved inorganic carbon (DIC) as our case study. We evaluated the method by extracting artificial float profiles from a depth‐resolved observation‐based DIC product to generate an artificial DIC time series. We then estimated NCP before and after accounting for horizontal gradient effects and compared the results to NCP estimates from an artificial DIC time series extracted at a fixed location along the float trajectory. Testing 10 biogeographical domains with moderate to substantial horizontal DIC gradients, our method significantly improved the precision (by ∼50 to ∼80%) and accuracy (by ∼10 to ∼100%) of regional NCP estimates. This method can be applied to other tracers with multi‐month‐long residence times. 
    more » « less
  3. ; ; ; (, Journal of Geophysical Research: Oceans)
    Abstract Oceanic transient tracers, such as chlorofluorocarbons (CFCs) and sulfur‐hexafluoride (SF6), trace the propagation of intermediate‐to‐abyssal water masses in the ocean interior. Their temporal and spatial sparsity, however, has limited their utility in quantifying the global ocean circulation and its decadal variability. TheTime‐Correction Method(TCM) presented here is a new approach to leverage the available CFCs and SF6observations to solve for the Green's functions (GFs) describing the steady‐state transport from the surface to the ocean interior. From the GFs, we reconstruct global tracer concentrations (and associated uncertainties) in the ocean interior at annual resolution (1940–2021). The spatial resolution includes 50 neutral density levels that span the water column along World Ocean Circulation Experiment/Global Ocean Ship‐Based Hydrographic Investigations Program lines. The reconstructed tracer concentrations return a global view of CFCs and SF6spreading into new regions of the interior ocean, such as the deep north‐western Pacific. For example, they capture the southward spreading and equatorial recirculation of distinct North Atlantic Deep Water components, and the spreading of CFC‐rich Antarctic Bottom Water out of the Southern Ocean and into the North Pacific, East Indian, and West Atlantic. The reconstructed tracer concentrations fit the data in most locations (∼75%), indicating that a steady‐state circulation holds for the most part. Discrepancies between the reconstructed and observed concentrations offer insight into ventilation rate changes on decadal timescales. As an example, we infer decadal changes in Subantartic Mode Water (SAMW) and find an increase in SAMW ventilation from 1992 to 2014, highlighting the skill of the TCM in leveraging the sparse tracer observations. 
    more » « less
  4. ; ; ; ; ; ; ; ; ; ; et al (, Geoscientific Model Development)
    null (Ed.)
    Abstract. The University of Victoria Earth System Climate Model (UVic ESCM) ofintermediate complexity has been a useful tool in recent assessments oflong-term climate changes, including both paleo-climate modelling anduncertainty assessments of future warming. Since the last official releaseof the UVic ESCM 2.9 and the two official updates during the last decade,considerable model development has taken place among multiple researchgroups. The new version 2.10 of the University of Victoria Earth SystemClimate Model presented here will be part of the sixth phaseof the Coupled Model Intercomparison Project (CMIP6). More precisely it willbe used in the intercomparison of Earth system models of intermediatecomplexity (EMIC), such as the C4MIP, the Carbon Dioxide Removal and ZeroEmissions Commitment model intercomparison projects (CDR-MIP and ZECMIP,respectively). It now brings together and combines multiple modeldevelopments and new components that have come about since the lastofficial release of the model. The main additions to the base model are(i) an improved biogeochemistry module for the ocean, (ii) a vertically resolvedsoil model including dynamic hydrology and soil carbon processes, and (iii) arepresentation of permafrost carbon. To set the foundation of its use, wehere describe the UVic ESCM 2.10 and evaluate results from transienthistorical simulations against observational data. We find that the UVicESCM 2.10 is capable of reproducing changes in historical temperature andcarbon fluxes well. The spatial distribution of many ocean tracers,including temperature, salinity, phosphate and nitrate, also agree well withobserved tracer profiles. The good performance in the ocean tracers isconnected to an improved representation of ocean physical properties. Forthe moment, the main biases that remain are a vegetation carbon density thatis too high in the tropics, a higher than observed change in the ocean heatcontent (OHC) and an oxygen utilization in the Southern Ocean that is too low.All of these biases will be addressed in the next updates to the model. 
    more » « less
  5. ; ; ; ; (, Geophysical Research Letters)
    Abstract Stirring of water by mesoscale currents (“eddies”) leads to large‐scale transport of many important oceanic properties (“tracers”). These eddy‐induced transports can be related to the large‐scale tracer gradients, using the concept of turbulent diffusion. The concept is widely used to describe these transports in the real ocean and to represent them in climate models. This study focuses on the inherent complexity of the corresponding coefficient tensor (“K‐tensor”) and its components, defined here in all its spatio‐temporal complexity. Results demonstrate that this comprehensiveK‐tensor is space‐, time‐, direction‐ and tracer‐dependent. Using numerical simulations with both idealized and comprehensive models of the Atlantic circulation, we show that these properties lead to upgradient eddy fluxes and the potential importance of all tensor components. The uncovered complexity of the eddy transports calls for reconsideration of how they are estimated in practice, included in the general circulation models and theoretically interpreted. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email