Preliminary evaluation of the Simplified Higher Order Closure (SHOC) in the NASA Goddard Earth Observing System (GEOS)

Here we evaluate alternative planetary boundary layer (PBL) schemes in the NASA Goddard Earth Observing System (GEOS) model. The model currently calculates turbulent diffusivities with a local Richardson number-based scheme for stable layers, and a K-profile scheme with surface- and cloud-driven plumes based on Lock et al. (2000) for unstable layers. In this study, these are replaced with the TKE-based Simplified Higher Order Closure (SHOC). SHOC combines a 1.5-order turbulence closure with a double-Gaussian PDF serving as a subgrid condensation scheme. The new diffusivity calculations are tested with two new shallow mass flux schemes. First, the University of Washington shallow cumulus scheme of Park and Bretherton (2009), a buoyancy-sorting entraining-detraining bulk plume with overshooting entrainment and a TKE/CIN-based closure. Second, an Eddy Diffusivity Mass Flux (EDMF) scheme based on the multiple plume model of Suselj et al. (2013) with stochastic lateral entrainment. Global climatologies of low level temperature, humidity, and wind are evaluated against MERRA-2 and ERA-Interim reanalyses, and significant reductions in bias are seen with both mass flux schemes. The diurnal cycle over the continental United States is evaluated using MERRA-2 and an in situ network of surface Lidar profilers.