In Part 2 of this two‐part paper, documentation is provided of key aspects of a version of the AM4.0/LM4.0 atmosphere/land model that will serve as a base for a new set of climate and Earth system models (CM4 and ESM4) under development at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). The quality of the simulation in AMIP (Atmospheric Model Intercomparison Project) mode has been provided in Part 1. Part 2 provides documentation of key components and some sensitivities to choices of model formulation and values of parameters, highlighting the convection parameterization and orographic gravity wave drag. The approach taken to tune the model's clouds to observations is a particular focal point. Care is taken to describe the extent to which aerosol effective forcing and Cess sensitivity have been tuned through the model development process, both of which are relevant to the ability of the model to simulate the evolution of temperatures over the last century when coupled to an ocean model.

This study compares five planetary boundary-layer (PBL) parametrizations in the Weather Research and Forecasting (WRF) numerical model for a single day from the Cooperative Atmosphere-Surface Exchange Study (CASES-99) field program. The five schemes include two first-order closure schemes-the Yonsei University (YSU) PBL and Asymmetric Convective Model version 2 (ACM2), and three turbulent kinetic energy (TKE) closure schemes-the Mellor-Yamada-Janjić (MYJ), quasi-normal scale elimination (QNSE), and Bougeault-Lacarrére (BouLac) PBL. The comparison results reveal that discrepancies among thermodynamic surface variables from different schemes are large at daytime, while the variables converge at nighttime with large deviations from those observed. On the other hand, wind components are more divergent at nighttime with significant biases. Regarding PBL structures, a non-local scheme with the entrainment flux proportional to the surface flux is favourable in unstable conditions. In stable conditions, the local TKE closure schemes show better performance. The sensitivity of simulated variables to surface-layer parametrizations is also investigated to assess relative contributions of the surface-layer parametrizations to typical features of each PBL scheme. In the surface layer, temperature and moisture are more strongly influenced by surface-layer formulations than by PBL mixing algorithms in both convective and stable regimes, while wind speed depends on vertical diffusion formulations in the convective regime. Regarding PBL structures, surface-layer formulations only contribute to near-surface variability and then PBL mean properties, whereas shapes of the profiles are determined by PBL mixing algorithms.

In this two‐part paper, a description is provided of a version of the AM4.0/LM4.0 atmosphere/land model that will serve as a base for a new set of climate and Earth system models (CM4 and ESM4) under development at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). This version, with roughly 100 km horizontal resolution and 33 levels in the vertical, contains an aerosol model that generates aerosol fields from emissions and a “light” chemistry mechanism designed to support the aerosol model but with prescribed ozone. In Part 1, the quality of the simulation in AMIP (Atmospheric Model Intercomparison Project) mode—with prescribed sea surface temperatures (SSTs) and sea‐ice distribution—is described and compared with previous GFDL models and with the CMIP5 archive of AMIP simulations. The model's Cess sensitivity (response in the top‐of‐atmosphere radiative flux to uniform warming of SSTs) and effective radiative forcing are also presented. In Part 2, the model formulation is described more fully and key sensitivities to aspects of the model formulation are discussed, along with the approach to model tuning.

Parameterization of the unresolved vertical transport in the planetary boundary layer (PBL) is one of the key physics algorithms in atmospheric models. This study attempts to represent the subgrid-scale (SGS) turbulent transport in convective boundary layers (CBLs) at gray-zone resolutions by investigating the effects of grid-size dependency in the vertical heat transport parameterization for CBL simulations. The SGS transport profile is parameterized based on the 2013 conceptual derivation by Shin and Hong. First, nonlocal transport via strong updrafts and local transport via the remaining small-scale eddies are separately calculated. Second, the SGS nonlocal transport is formulated by multiplying a grid-size dependency function with the total nonlocal transport profile fit to the large-eddy simulation (LES) output. Finally, the SGS local transport is formulated by multiplying a grid-size dependency function with the total local transport profile, which is calculated using an eddy-diffusivity formula. The new algorithm is evaluated against the LES output and compared with a conventional nonlocal PBL parameterization. For ideal CBL cases, by considering the scale dependency in the parameterized vertical heat transport, improvements over the conventional nonlocal K-profile model appear in mean profiles, resolved and SGS vertical transport profiles with their grid-size dependency, and the energy spectrum. Real-case simulations for convective rolls show that the simulated roll structures are more robust with stronger intensity when the new algorithm is used.

The gray zone of a physical process in numerical models is defined as the range of model resolution in which the process is partly resolved by model dynamics and partly parameterized. In this study, the authors examine the grid-size dependencies of resolved and parameterized vertical transports in convective boundary layers (CBLs) for horizontal grid scales including the gray zone. To assess how stability alters the dependencies on grid size, four CBLs with different surface heating and geostrophic winds are considered. For this purpose, reference data for grid-scale (GS) and subgrid-scale (SGS) fields are constructed for 50–4000-m mesh sizes by filtering 25-m large-eddy simulation (LES) data. As relative importance of shear increases, the ratio of resolved turbulent kinetic energy increases for a given grid spacing. Vertical transports of potential temperature, momentum, and a bottom-up diffusion passive scalar behave in a similar fashion. The effects of stability are related to the horizontal scale of coherent large-eddy structures that change in the different stability. The subgrid-scale vertical transport of heat and the bottom-up scalar are divided into a nonlocal mixing owing to the coherent structures and remaining local mixing. The separate treatment of the nonlocal and local transports shows that the grid-size dependency of the SGS nonlocal flux and its sensitivity to stability predominantly determine the dependency of total (nonlocal plus local) SGS transport.