James R. Overfelt

Aquaplanet Experiments Using CAM’s Variable-Resolution Dynamical Core

AbstractA variable-resolution option has been added within the spectral element (SE) dynamical core of the U.S. Department of Energy (DOE)–NCAR Community Atmosphere Model (CAM). CAM-SE allows for static refinement via conforming quadrilateral meshes on the cubed sphere. This paper investigates the effect of mesh refinement in a climate model by running variable-resolution (var-res) simulations on an aquaplanet. The variable-resolution grid is a 2° (~222 km) grid with a refined patch of 0.25° (~28 km) resolution centered at the equator. Climatology statistics from these simulations are compared to globally uniform runs of 2° and 0.25°.A significant resolution dependence exists when using the CAM version 4 (CAM4) subgrid physical parameterization package across scales. Global cloud fraction decreases and equatorial precipitation increases with finer horizontal resolution, resulting in drastically different climates between the uniform grid runs and a physics-induced grid imprinting in the var-res simulation...

Elastic wave propagation in variable media using a discontinuous Galerkin method.

Motivated by the needs of seismic inversion and building on our prior experience for fluid-dynamics systems, we present a high-order discontinuous Galerkin (DG) Runge-Kutta method applied to isotropic, linearized elasto-dynamics. Unlike other DG methods recently presented in the literature, our method allows for inhomogeneous material variations within each element that enables representation of realistic earth models - a feature critical for future use in seismic inversion. Likewise, our method supports curved elements and hybrid meshes that include both simplicial and nonsimplicial elements. We demonstrate the capabilities of this method through a series of numerical experiments including hybrid mesh discretizations of the Marmousi2 model as well as a modified Marmousi2 model with a oscillatory ocean bottom that is exactly captured by our discretization. A discontinuous Galerkin method for solving the equations of linear isotropic elasticity has been presented. The formulation is designed to accommodate variation of media parameters within elements, curved elements and unstructured heterogeneous meshes. We have demonstrated that each of these important features of the formulation can produce results that are significantly different from formulations that do not possess these capabilities suggesting that each of these capabilities may be important for effective full waveform inversion of elastic medium.

CAM-SE–CSLAM: Consistent Coupling of a Conservative Semi-Lagrangian Finite-Volume Method with Spectral Element Dynamics

AbstractAn algorithm to consistently couple a conservative semi-Lagrangian finite-volume transport scheme with a spectral element (SE) dynamical core is presented. The semi-Lagrangian finite-volume scheme is the Conservative Semi-Lagrangian Multitracer (CSLAM), and the SE dynamical core is the National Center for Atmospheric Research (NCAR)’s Community Atmosphere Model–Spectral Elements (CAM-SE). The primary motivation for coupling CSLAM with CAM-SE is to accelerate tracer transport for multitracer applications. The coupling algorithm result is an inherently mass-conservative, shape-preserving, and consistent (for a constant mixing ratio, the CSLAM solution reduces to the SE solution for air mass) transport that is efficient and accurate. This is achieved by first deriving formulas for diagnosing SE airmass flux through the CSLAM control volume faces. Thereafter, the upstream Lagrangian CSLAM areas are iteratively perturbed to match the diagnosed SE airmass flux, resulting in an equivalent upstream Lagran...

Risk assessment of climate systems for national security.

Climate change, through drought, flooding, storms, heat waves, and melting Arctic ice, affects the production and flow of resource within and among geographical regions. The interactions among governments, populations, and sectors of the economy require integrated assessment based on risk, through uncertainty quantification (UQ). This project evaluated the capabilities with Sandia National Laboratories to perform such integrated analyses, as they relate to (inter)national security. The combining of the UQ results from climate models with hydrological and economic/infrastructure impact modeling appears to offer the best capability for national security risk assessments.