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 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
  3. ; ; ; ; ; ; ; ; ; ; 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
  4. ; ; ; ; (, 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
  5. ; ; ; (, Journal of Physical Oceanography)
    Abstract The age of seawater refers to the amount of time that has elapsed since that water encountered the surface. This age measures the ventilation rate of the ocean, and the spatial distribution of age can be influenced by multiple processes, such as overturning circulation, ocean mixing, and air–sea exchange. In this work, we aim to gain new quantitative insights about how the ocean’s age tracer distribution reflects the strength of the meridional overturning circulation and diapycnal diffusivity. We propose an integral constraint that relates the age tracer flow across an isopycnal surface to the geometry of the surface. With the integral constraint, a relationship between the globally averaged effective diapycnal diffusivity and the meridional overturning strength at an arbitrary density level can be inferred from the age tracer concentration near that level. The theory is tested in a set of idealized single-basin simulations. A key insight from this study is that the age difference between regions of upwelling and downwelling, rather than any single absolute age value, is the best indicator of overturning strength. The framework has also been adapted to estimate the strength of abyssal overturning circulation in the modern North Pacific, and we demonstrate that the age field provides an estimate of the circulation strength consistent with previous studies. This framework could potentially constrain ocean circulation and mixing rates from age-like realistic tracers (e.g., radiocarbon) in both past and present climates. Significance StatementThe age of seawater—the local mean time since local water from different pathways was last at the surface—is a valuable indicator of ocean circulation and the transport time scale of heat and carbon. We introduce a novel constraint that relates total age flow across a density surface to its geometry, which provides new insights into constraining ocean circulation and mixing rates from age-like realistic tracers (e.g., radiocarbon). 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email