skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Evaluation of the accuracy of an offline seasonally-varying matrix transport model for simulating ideal age

Abstract

Newton-Krylov solvers for ocean tracers have the potential to greatly decrease the computational costs of spinning up deep-ocean tracers, which can take several thousand model years to reach equilibrium with surface processes. One version of the algorithm uses offline tracer transport matrices to simulate an annual cycle of tracer concentrations and applies Newton’s method to find concentrations that are periodic in time. Here we present the impact of time-averaging the transport matrices on the equilibrium values of an ideal-age tracer. We compared annually-averaged, monthly-averaged, and 5-day-averaged transport matrices to an online simulation using the ocean component of the Community Earth System Model (CESM) with a nominal horizontal resolution of 1° × 1° and 60 vertical levels. We found that increasing the time resolution of the offline transport model reduced a low age bias from 12% for the annually-averaged transport matrices, to 4% for the monthly-averaged transport matrices, and to less than 2% for the transport matrices constructed from 5-day averages. The largest differences were in areas with strong seasonal changes in the circulation, such as the Northern Indian Ocean. As a result, for many applications the relatively small bias obtained using the offline model makes the offline approach attractive becausemore » it uses significantly less computer resources and is simpler to set up and run.« less

Authors:
 [1];  [1];  [2];  [3]
  1. Univ. of California, Irvine, CA (United States)
  2. National Center for Atmospheric Research, Boulder, CO (United States)
  3. Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
Publication Date:
Research Org.:
Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA)
OSTI Identifier:
1275512
Alternate Identifier(s):
OSTI ID: 1398012
Report Number(s):
SAND-2016-6880J
Journal ID: ISSN 1463-5003; PII: S1463500316300750
Grant/Contract Number:
AC04-94AL85000; SC0012550
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Ocean Modelling
Additional Journal Information:
Journal Volume: 105; Journal Issue: C; Journal ID: ISSN 1463-5003
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
58 GEOSCIENCES; 54 ENVIRONMENTAL SCIENCES; implicit solver; ideal age; Newton–Krylov; offline; transport matrix; global ocean modeling

Citation Formats

Bardin, Ann, Primeau, Francois, Lindsay, Keith, and Bradley, Andrew. Evaluation of the accuracy of an offline seasonally-varying matrix transport model for simulating ideal age. United States: N. p., 2016. Web. doi:10.1016/j.ocemod.2016.07.003.
Bardin, Ann, Primeau, Francois, Lindsay, Keith, & Bradley, Andrew. Evaluation of the accuracy of an offline seasonally-varying matrix transport model for simulating ideal age. United States. doi:10.1016/j.ocemod.2016.07.003.
Bardin, Ann, Primeau, Francois, Lindsay, Keith, and Bradley, Andrew. Thu . "Evaluation of the accuracy of an offline seasonally-varying matrix transport model for simulating ideal age". United States. doi:10.1016/j.ocemod.2016.07.003. https://www.osti.gov/servlets/purl/1275512.
@article{osti_1275512,
title = {Evaluation of the accuracy of an offline seasonally-varying matrix transport model for simulating ideal age},
author = {Bardin, Ann and Primeau, Francois and Lindsay, Keith and Bradley, Andrew},
abstractNote = {Newton-Krylov solvers for ocean tracers have the potential to greatly decrease the computational costs of spinning up deep-ocean tracers, which can take several thousand model years to reach equilibrium with surface processes. One version of the algorithm uses offline tracer transport matrices to simulate an annual cycle of tracer concentrations and applies Newton’s method to find concentrations that are periodic in time. Here we present the impact of time-averaging the transport matrices on the equilibrium values of an ideal-age tracer. We compared annually-averaged, monthly-averaged, and 5-day-averaged transport matrices to an online simulation using the ocean component of the Community Earth System Model (CESM) with a nominal horizontal resolution of 1° × 1° and 60 vertical levels. We found that increasing the time resolution of the offline transport model reduced a low age bias from 12% for the annually-averaged transport matrices, to 4% for the monthly-averaged transport matrices, and to less than 2% for the transport matrices constructed from 5-day averages. The largest differences were in areas with strong seasonal changes in the circulation, such as the Northern Indian Ocean. As a result, for many applications the relatively small bias obtained using the offline model makes the offline approach attractive because it uses significantly less computer resources and is simpler to set up and run.},
doi = {10.1016/j.ocemod.2016.07.003},
journal = {Ocean Modelling},
number = C,
volume = 105,
place = {United States},
year = {Thu Jul 21 00:00:00 EDT 2016},
month = {Thu Jul 21 00:00:00 EDT 2016}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 2works
Citation information provided by
Web of Science

Save / Share:
  • Leaf area index is an important input for many climate and carbon models. The widely used leaf area products derived from satellite-observed surface reflectances contain substantial erratic fluctuations in time due to inadequate atmospheric corrections and observational and retrieval uncertainties. These fluctuations are inconsistent with the seasonal dynamics of leaf area, known to be gradual. Their use in process-based terrestrial carbon models corrupts model behavior, making diagnosis of model performance difficult. We propose a data assimilation approach that combines the satellite observations of Moderate Resolution Imaging Spectroradiometer (MODIS) albedo with a dynamical leaf model. Its novelty is that the seasonalmore » cycle of the directly retrieved leaf areas is smooth and consistent with both observations and current understandings of processes controlling leaf area dynamics. The approach optimizes the dynamical model parameters such that the difference between the estimated surface reflectances based on the modeled leaf area and those of satellite observations is minimized. We demonstrate the usefulness and advantage of our new approach at multiple deciduous forest sites in the United States.« less
  • There has been renewed interest in the application of functional models to the transport of nonpoint source pollutants at polypedon and watershed scales because of the ease of their coupling to a geographic information system and to the accepted organizational hierarchy of pedogenetic modeling approaches. However, very little work has been done to evaluate the performance of a functional transient-state model for the transport of a reactive solute over an extensive study period. Subsequently, the functional model TETrans (Trace Element Transport) was evaluated for model performance with boron (B) transport data collected from a meso-scale soil lysimeter column over amore » 1,000-day study period. Because the ability to simulate water flow has been evaluated previously for TETrans, the focus of this evaluation centered around the performance of various functional models of B adsorption used as subroutines within the TETrans model, including the (1) Freundlich, (2) kinetic Freundlich, (3) Langmuir, (4) temperature-dependent Langmuir, and (5) pH-dependent Keren adsorption isotherm equations. Model performance was evaluated with statistical functions, specifically the Average Absolute Prediction Error, the Root Mean Square Error, the Reduced Error Estimate and the Coefficient of Residual Mass, and graphic displays of observed and predicted B concentration profiles. Even though no single adsorption isotherm equation, when coupled to TETrans, could be considered poor in its performance, results indicated that the order of model performance was the pH-dependent Keren equation first, followed by the temperature-dependent Langmuir and kinetic Freundlich equations, the Freundlich equation, and, finally, the Langmuir equation. Overall, the TETrans model was able to simulate the transport of B with deviations because no functional adsorption equation incorporated all the influences of pH, ionic strength, temperature, and kinetic effects into a single equation. The inability to correctly predict one of the measured peaks in B concentration near the soil surface suggests that problems with the timing of the sample collection may have occurred for the shallowest sampling depth.« less
  • Purpose: To evaluate the utility of intraprostatic markers in the treatment verification of prostate cancer radiotherapy. Specific aims were: to compare the effectiveness of offline correction protocols, either using gold markers or bony anatomy; to estimate the potential benefit of online correction protocol's using gold markers; to determine the presence and effect of intrafraction motion. Methods and Materials: Thirty patients with three gold markers inserted had pretreatment and posttreatment images acquired and were treated using an offline correction protocol and gold markers. Retrospectively, an offline protocol was applied using bony anatomy and an online protocol using gold markers. Results: Themore » systematic errors were reduced from 1.3, 1.9, and 2.5 mm to 1.1, 1.1, and 1.5 mm in the right-left (RL), superoinferior (SI), and anteroposterior (AP) directions, respectively, using the offline correction protocol and gold markers instead of bony anatomy. The subsequent decrease in margins was 1.7, 3.3, and 4 mm in the RL, SI, and AP directions, respectively. An offline correction protocol combined with an online correction protocol in the first four fractions reduced random errors further to 0.9, 1.1, and 1.0 mm in the RL, SI, and AP directions, respectively. A daily online protocol reduced all errors to <1 mm. Intrafraction motion had greater impact on the effectiveness of the online protocol than the offline protocols. Conclusions: An offline protocol using gold markers is effective in reducing the systematic error. The value of online protocols is reduced by intrafraction motion.« less
  • Recent model studies have indicated that observed stratospheric ozone decline can have a cooling effect on climate. This study intends to investigate the climate response due to changes in the radiative fluxes caused by prescribed changes in stratospheric as well as tropospheric ozone. For this purpose the authors use a simplified climate model, basically consisting of an energy balance atmosphere model coupled to an advection-diffusion ocean model. The coupled climate model simulates the latitudinal and seasonal variations in zonal mean surface air temperature and the average lower (12-22 km) and higher (22-100 km) stratospheric temperatures. First, the quasi-equilibrium response ofmore » the model to various uniform ozone perturbations is examined. For instance, a uniform 50% reduction in lower stratospheric ozone results in a global average cooling of 3.5{degrees}C in the lower stratosphere with maximum values in the Tropics and of 0.46{degrees}C at the surface with maximum cooling in the polar winter. The latter is largely due to the albedo-temperature feedback, mainly through increases in sea ice. The albedo-temperature feedback is consistently stronger in the case of tropospheric and lower stratospheric ozone perturbations than in the case of, for instance, CO{sub 2} perturbations. This can be attributed mainly to differences in the meridional gradient in tropopause radiative forcing. This study indicates that one must be cautious when using concepts such as global radiative forcing and global climate sensitivity in quantifying climate change. Finally, the transient model response to various ozone trend scenarios indicates that the net effect of tropospheric ozone increases and stratospheric ozone depletions is a slight global average cooling (-0.001 to -0.003 K yr{sup -1}), which offsets by approximately 10% the projected surface warming due to increases in the other greenhouse gases.« less