skip to main content


Title: Reconciling and Improving Formulations for Thermodynamics and Conservation Principles in Earth System Models (ESMs)
Abstract

This paper provides a comprehensive derivation of the total energy equations for the atmospheric components of Earth System Models (ESMs). The assumptions and approximations made in this derivation are motivated and discussed. In particular, it is emphasized that closing the energy budget is conceptually challenging and hard to achieve in practice without resorting to ad hoc fixers. As a concrete example, the energy budget terms are diagnosed in a realistic climate simulation using a global atmosphere model. The largest total energy errors in this example are spurious dynamical core energy dissipation, thermodynamic inconsistencies (e.g., coupling parameterizations with the host model) and missing processes/terms associated with falling precipitation and evaporation (e.g., enthalpy flux between components). The latter two errors are not, in general, reduced by increasing horizontal resolution. They are due to incomplete thermodynamic and dynamic formulations. Future research directions are proposed to reconcile and improve thermodynamics formulations and conservation principles.

 
more » « less
PAR ID:
10444539
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Advances in Modeling Earth Systems
Volume:
14
Issue:
9
ISSN:
1942-2466
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Analyses of atmospheric heat and moisture budgets serve as an effective tool to study convective characteristics over a region and to provide large‐scale forcing fields for various modeling applications. This paper examines two popular methods for computing large‐scale atmospheric budgets: the conventional budget method (CBM) using objectively gridded analyses based primarily on radiosonde data and the constrained variational analysis (CVA) approach which supplements vertical profiles of atmospheric fields with measurements at the top of the atmosphere and at the surface to conserve mass, water, energy, and momentum. Successful budget computations are dependent on accurate sampling and analyses of the thermodynamic state of the atmosphere and the divergence field associated with convection and the large‐scale circulation that influences it. Utilizing analyses generated from data taken during Dynamics of the Madden‐Julian Oscillation (DYNAMO) field campaign conducted over the central Indian Ocean from October to December 2011, we evaluate the merits of these budget approaches and examine their limitations. While many of the shortcomings of the CBM, in particular effects of sampling errors in sounding data, are effectively minimized with CVA, accurate large‐scale diagnostics in CVA are dependent on reliable background fields and rainfall constraints. For the DYNAMO analyses examined, the operational model fields used as the CVA background state provided wind fields that accurately resolved the vertical structure of convection in the vicinity of Gan Island. However, biases in the model thermodynamic fields were somewhat amplified in CVA resulting in a convective environment much weaker than observed.

     
    more » « less
  2. Abstract Optimal control theory can be a useful tool to identify the best strategies for the management of infectious diseases. In most of the applications to disease control with ordinary differential equations, the objective functional to be optimized is formulated in monetary terms as the sum of intervention costs and the cost associated with the burden of disease. We present alternate formulations that express epidemiological outcomes via health metrics and reframe the problem to include features such as budget constraints and epidemiological targets. These alternate formulations are illustrated with a compartmental cholera model. The alternate formulations permit us to better explore the sensitivity of the optimal control solutions to changes in available budget or the desired epidemiological target. We also discuss some limitations of comprehensive cost assessment in epidemiology. 
    more » « less
  3. ABSTRACT

    A popular numerical method to model the dynamics of a ‘full spectrum’ of cosmic rays (CRs), also applicable to radiation/neutrino hydrodynamics, is to discretize the spectrum at each location/cell as a piecewise power law in ‘bins’ of momentum (or frequency) space. This gives rise to a pair of conserved quantities (e.g. CR number and energy) that are exchanged between cells or bins, which in turn give the update to the normalization and slope of the spectrum in each bin. While these methods can be evolved exactly in momentum-space (e.g. considering injection, absorption, continuous losses/gains), numerical challenges arise dealing with spatial fluxes, if the scattering rates depend on momentum. This has often been treated either by neglecting variation of those rates ‘within the bin,’ or sacrificing conservation – introducing significant errors. Here, we derive a rigorous treatment of these terms, and show that the variation within the bin can be accounted for accurately with a simple set of scalar correction coefficients that can be written entirely in terms of other, explicitly evolved ‘bin-integrated’ quantities. This eliminates the relevant errors without added computational cost, has no effect on the numerical stability of the method, and retains manifest conservation. We derive correction terms both for methods that explicitly integrate flux variables (e.g. two-moment or M1-like) methods, as well as single-moment (advection-diffusion, FLD-like) methods, and approximate corrections valid in various limits.

     
    more » « less
  4. Abstract

    The water budget components of an atmospheric column are precipitation, evaporation, and horizontal water vapor divergence. This study finds that when precipitation from the Global Precipitation Climatology Project (GPCP), evaporation from the SeaFlux product, and water vapor divergence from ERA5 are employed, the degree of budget closure depends strongly on the location and time period. Variations in this error are not random, and this study seeks to better understand these biases as the climate system evolves. Errors are particularly significant over ocean regions in and near the west Pacific warm pool, where there are multiyear budget residuals of roughly 10% of the magnitude of precipitation. Biases in other tropical ocean basins are more seasonal and smaller in magnitude. Time-varying budget errors are strongly linked to variations in convective organization that vary with the large-scale environment; errors correlating with deep organized rain show coefficients of 0.62 and 0.56 in the west Pacific and central Pacific, respectively. Errors are linked to a lesser extent with cloud microphysical structures in the East Indian. Both factors affect rainfall retrievals through well-known precipitation bias mechanisms (beam-filling, convective/stratiform effects). Characteristics of the large-scale environment that produce changes in convective organization and precipitable ice water content are explored.

     
    more » « less
  5. Abstract For the Community Atmosphere Model version 6 (CAM6), an adjustment is needed to conserve dry air mass. This adjustment exposes an inconsistency in how CAM6’s energy budget incorporates water—in CAM6 water in the vapor phase has energy, but condensed phases of water do not. When water vapor condenses, only its latent energy is retained in the model, while its remaining internal, potential, and kinetic energy are lost. A global fixer is used in the default CAM6 model to maintain global energy conservation, but locally the energy tendency associated with water changing phase violates the divergence theorem. This error in energy tendency is intrinsically tied to the water vapor tendency, and reaches its highest values in regions of heavy rainfall, where the error can be as high as 40 W m −2 annually averaged. Several possible changes are outlined within this manuscript that would allow CAM6 to satisfy the divergence theorem locally. These fall into one of two categories: 1) modifying the surface flux to balance the local atmospheric energy tendency and 2) modifying the local atmospheric tendency to balance the surface plus top-of-atmosphere energy fluxes. To gauge which aspects of the simulated climate are most sensitive to this error, the simplest possible change—where condensed water still does not carry energy and a local energy fixer is used in place of the global one—is implemented within CAM6. Comparing this experiment with the default configuration of CAM6 reveals precipitation, particularly its variability, to be highly sensitive to the energy budget formulation. Significance Statement This study examines and explains spurious regional sources and sinks of energy in a widely used climate model. These energy errors result from not tracking energy associated with water after it transitions from the vapor phase to either liquid or ice. Instead, the model used a global fixer to offset the energy tendency related to the energy sources and sinks associated with condensed water species. We replace this global fixer with a local one to examine the model sensitivity to the regional energy error and find a large sensitivity in the simulated hydrologic cycle. This work suggests that the underlying thermodynamic assumptions in the model should be revisited to build confidence in the model-simulated regional-scale water and energy cycles. 
    more » « less