Abstract. Climate simulation codes, such as the Community Earth System Model (CESM), are especially complex and continually evolving. Their ongoing state of development requires frequent software verification in the form of quality assurance to both preserve the quality of the code and instill model confidence. To formalize and simplify this previously subjective and computationally expensive aspect of the verification process, we have developed a new tool for evaluating climate consistency. Because an ensemble of simulations allows us to gauge the natural variability of the model's climate, our new tool uses an ensemble approach for consistency testing. In particular, an ensemble of CESM climate runs is created, from which we obtain a statistical distribution that can be used to determine whether a new climate run is statistically distinguishable from the original ensemble. The CESM ensemble consistency test, referred to as CESM-ECT, is objective in nature and accessible to CESM developers and users. The tool has proven its utility in detecting errors in software and hardware environments and providing rapid feedback to model developers.

No abstract

A 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 28 (;222 km) grid with a refined patch of 0.258 (;28 km) resolution centered at the equator. Climatology statistics from these simulations are compared to globally uniform runs of 28 and 0.258. 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. Using CAM version 5 (CAM5) physics significantly improves cloud scaling at different grid resolutions. Additional precipitation at the equator in the highresolution mesh results in collocated zonally anomalous divergence in both var-res simulations, although this feature is much weaker in CAM5 than CAM4. The equilibrium solution at each grid spacing within the var-res simulations captures the majority of the resolution signal of the corresponding globally uniform grids. The var-res simulation exhibits good performance with respect to wave propagation, including equatorial regions where waves pass through grid transitions. In addition, the increased frequency of high-precipitation events in the refined 0.258 area within the var-res simulations matches that observed in the global 0.258 simulations.

We evaluate the Community ocean Vertical Mixing project version of the K‐profile parameterization (KPP) for modeling upper ocean turbulent mixing. For this purpose, one‐dimensional KPP simulations are compared across a suite of oceanographically relevant regimes against horizontally averaged large eddy simulations (LESs). We find the standard configuration of KPP consistent with LES across many forcing regimes, supporting its physical basis. Our evaluation also motivates recommendations for KPP best practices within ocean circulation models and identifies areas where further research is warranted. The original treatment of KPP recommends the matching of interior diffusivities and their gradients to the KPP‐predicted values computed in the ocean surface boundary layer (OSBL). However, we find that difficulties in representing derivatives of rapidly changing diffusivities near the base of the OSBL can lead to loss of simulation fidelity. To mitigate this difficulty, we propose and evaluate two computationally simpler approaches: (1) match to the internal predicted diffusivity alone and (2) set the KPP diffusivity to 0 at the OSBL base. We find the KPP entrainment buoyancy flux to be sensitive to vertical grid resolution and details of how to diagnose the KPP boundary layer depth. We modify the KPP turbulent shear velocity parameterization to reduce resolution dependence. Additionally, an examination of LES vertical turbulent scalar flux budgets shows that the KPP‐parameterized nonlocal tracer flux is incomplete due to the assumption that it solely redistributes the surface tracer flux. This result motivates further studies of the nonlocal flux parameterization.