This work documents the first version of the U.S. Department of Energy (DOE) new EnergyExascale Earth System Model (E3SMv1). We focus on the standard resolution of the fully coupled physical model designed to address DOE mission-relevant water cycle questions. Its components include atmosphere and land (110-km grid spacing), ocean and sea ice (60 km in the midlatitudes and 30 km at the equator and poles), and river transport (55 km) models. This base configuration will also serve as a foundation for additional configurations exploring higher horizontal resolution as well as augmented capabilities in the form of biogeochemistry and cryosphere configurations. The performance of E3SMv1 is evaluated by means of a standard set of Coupled Model Intercomparison Project Phase 6 (CMIP6) Diagnosis, Evaluation, and Characterization of Klima simulations consisting of a long preindustrial control, historical simulations (ensembles of fully coupled and prescribed SSTs) as well as idealized CO 2 forcing simulations. The model performs well overall with biases typical of other CMIP-class models, although the simulated Atlantic Meridional Overturning Circulation is weaker than many CMIP-class models. While the E3SMv1 historical ensemble captures the bulk of the observed warming between preindustrial (1850) and present day, the trajectory of the warming diverges from observations in the Key Points: • This work documents E3SMv1, the first version of the U.S. DOE Energy Exascale Earth System Model • The performance of E3SMv1 is documented with a set of standard CMIP6 DECK and historical simulations comprising nearly 3,000 years • E3SMv1 has a high equilibrium climate sensitivity (5.3 K) and strong aerosol-related effective radiative forcing (-1.65 W/m 2 ) Correspondence to: Chris Golaz, golaz1@llnl.gov Citation: Golaz, J.-C., Caldwell, P. M., Van Roekel, L. P., Petersen, M. R., Tang, Q., Wolfe, J. D., et al. (2019). The DOE E3SM coupled model version 1: Overview and evaluation at standard resolution. second half of the twentieth century with a period of delayed warming followed by an excessive warming trend. Using a two-layer energy balance model, we attribute this divergence to the model's strong aerosol-related effective radiative forcing (ERF ari+aci = −1.65 W/m 2 ) and high equilibrium climate sensitivity (ECS = 5.3 K). Plain Language Summary The U.S. Department of Energy funded the development of a new state-of-the-art Earth system model for research and applications relevant to its mission. The Energy Exascale Earth System Model version 1 (E3SMv1) consists of five interacting components for the global atmosphere, land surface, ocean, sea ice, and rivers. Three of these components (ocean, sea ice, and river) are new and have not been coupled into an Earth system model previously. The atmosphere and land surface components were created by extending existing components part of the Community Earth System Model, Version 1. E3SMv1's capabilities are demonstrated by performing a set of standardized simulation experiments described by...

The Trilinos Project is an effort to facilitate the design, development, integration and ongoing support of mathematical software libraries within an object-oriented framework for the solution of large-scale, complex multi-physics engineering and scientific problems. Trilinos addresses two fundamental issues of developing software for these problems: (i) Providing a streamlined process and set of tools for development of new algorithmic implementations and (ii) promoting interoperability of independently developed software. Trilinos uses a two-level software structure designed around collections of packages. A Trilinos package is an integral unit usually developed by a small team of experts in a particular algorithms area such as algebraic preconditioners, nonlinear solvers, etc. Packages exist underneath the Trilinos top level, which provides a common look-and-feel, including configuration, documentation, licensing, and bug-tracking.Here we present the overall Trilinos design, describing our use of abstract interfaces and default concrete implementations. We discuss the services that Trilinos provides to a prospective package and how these services are used by various packages. We also illustrate how packages can be combined to rapidly develop new algorithms. Finally, we discuss how Trilinos facilitates highquality software engineering practices that are increasingly required from simulation software.

The Energy Exascale Earth System Model Atmosphere Model version 1, the atmospheric component of the Department of Energy's Energy Exascale Earth System Model is described. The model began as a fork of the well‐known Community Atmosphere Model, but it has evolved in new ways, and coding, performance, resolution, physical processes (primarily cloud and aerosols formulations), testing and development procedures now differ significantly. Vertical resolution was increased (from 30 to 72 layers), and the model top extended to 60 km (~0.1 hPa). A simple ozone photochemistry predicts stratospheric ozone, and the model now supports increased and more realistic variability in the upper troposphere and stratosphere. An optional improved treatment of light‐absorbing particle deposition to snowpack and ice is available, and stronger connections with Earth system biogeochemistry can be used for some science problems. Satellite and ground‐based cloud and aerosol simulators were implemented to facilitate evaluation of clouds, aerosols, and aerosol‐cloud interactions. Higher horizontal and vertical resolution, increased complexity, and more predicted and transported variables have increased the model computational cost and changed the simulations considerably. These changes required development of alternate strategies for tuning and evaluation as it was not feasible to “brute force” tune the high‐resolution configurations, so short‐term hindcasts, perturbed parameter ensemble simulations, and regionally refined simulations provided guidance on tuning and parameterization sensitivity to higher resolution. A brief overview of the model and model climate is provided. Model fidelity has generally improved compared to its predecessors and the CMIP5 generation of climate models.

An adaptively stabilized finite element scheme is proposed for a strongly coupled hydro-mechanical problem in fluid-infiltrating porous solids at finite strain. We first present the derivation of the poromechanics model via mixture theory in large deformation. By exploiting assumed deformation gradient techniques, we develop a numerical procedure capable of simultaneously curing the multiple-locking phenomena related to shear failure, incompressibility imposed by pore fluid, and/or incompressible solid skeleton and produce solutions that satisfy the inf-sup condition. The template-based generic programming and automatic differentiation (AD) techniques used to implement the stabilized model are also highlighted. Finally, numerical examples are given to show the versatility and efficiency of this model. Darcy's velocity [7], by using inf-sup stable finite elements (e.g., Talyor-Hood, Raviart-Thomas finite elements) [9,[11][12][13][14], or by applying stabilization procedures to cure the otherwise unstable finite elements [7,[15][16][17][18]. The displacement-Darcy-velocity coupling scheme is relatively easy to implement. However, computation time will be significantly increased because of the extra degrees of freedom added for the pore fluid velocity. Implementing inf-sup stable displacement-pressure finite element requires only one extra degree of freedom for the pore pressure. However, inf-sup stable finite elements require special meshing, data structure, and preprocessing and postprocessing tools to accommodate the need to have different basis functions for the displacement and pore pressure solutions. This difficulty is accompanied with a significant increase in computation cost, as higher-order mixed finite elements lead to larger systems of equations and require more integration points to perform numerical integration, as pointed out in [14].An equal-order displacement-pore pressure mixed finite element formulation is computationally more efficient and does not require significant modification to data structures (except adding an additional degree of freedom to all nodes), and it is therefore more feasible for both code development and maintenance. The price for using the equal-order mixed finite element method, however, is that the mixed finite element formulation must be stabilized by either adding additional terms or introducing enhanced shape functions to eliminate the spurious mode encountered as a result of a failure to satisfy the inf-sup condition. In recent years, effort has been invested to develop stabilization procedures for poromechanics problem under the small strain assumption, for example, [7,[16][17][18]. However, there are numerous occasions in which porous solids experience significant deformation such that a finite strain formulation becomes essential. Examples include water-saturated soil near critical state [7] and hydrated biological tissue during normal physiological activities [19]. In those situations, a stabilized formulation for large deformation poromechanics inheriting the same ...

LOCA, the Library of Continuation Algorithms, is a software library for performing stability analysis of large-scale applications. LOC A enables the tracking of solution branches as a function of a system parameter, the direct tracking of bifurcation points, and, when linked with the ARPACK library, a linear stability analysis capability. It is designed to be easy to implement around codes that already use Newton's method to converge to steady-state solutions. The algorithms are chosen to work for large problems, such as those that arise from discretizations of partial differential equations, and to run on distributed memory parallel machines. This manual presents LOCA's continuation and bifurcation analysis algorithms, and instructions on how to implement LOCA with an application code. The LOCA code is being made publicly available at www . cs . sandia. gov/loca. 3