The lowest model layer is assumed to be the surface layer and apply the Monin-Obukhov similarity profile relationship to obtain the surface stress and sensible and latent heat fluxes. The formulation was based on Miyakoda and Sirutis (1986) [141] and has been modified by P. Long [124] [125] in the very stable and very unstable situations.
GFS Noah LSM Model are largely responsible for the quality of model forecasts produced for near-surface weather parameters, such as 2-meter air temperature ( \(T_{2m}\)) and surface skin temperature ( \(LST\)). \(LST\) is derived from the surface energy budget, and is particularly important to remote sensing and data assimilation. How precise these two parameters can be simulated by the model strongly depends on how accurate the surface heat fluxes are parameterized, particularly the surface sensible heat flux ( \(SH\)). The surface thermal roughness length is a key parameter to determine \(SH\). Previous GFS version do not distinguish between the roughness length for heat and momentum. The aerodynamic roughness \(Z_{0m}\) is used for wind, while the thermal roughness \(Z_{0t}\) is used for heat and water vapor. In the surface layer, the GFS applies the MO similarity profile scheme with modified stability functions (Miyakoda and Sirutis 1986 [141]; Long 1984, 1986 [124] [125]). Momentum and thermal roughness lengths are necessary to estimate the surface fluxes from the atmospheric surface layer similarity theory. In the current operational GFS, the momentum roughtness length \(Z_{0m}\) is specified according to the fixed vegetation types but has no seasonal variation.
In May 2011, the new vegetation-dependent formulations of thermal roughness formulation ( Zheng et al. (2012) [214]) was implemented to deal with the cold \(LST\) bias over the arid western continental United States (CONUS) during daytime. The thermal roughness length \(Z_{0H}\) is derived by a seasonlly varying formulation dependent on the seasonal cycle of green vegetation fraction. In this \(Z_{0H}\) formulation, a key parameter known as \(C_{zil}\) is specified according to a dependence on canopy height.
The NCEP GFS global prediction model has experienced a longstanding problem of severe cold bias in the \(T_{2m}\) forecasts over land in the late afternoon and nighttime during moist seasons. This cold bias is closely associated with the nocturnal stable boundary layer and is accompanied by a corresponding warm air temperature bias in the first model level above the ground. In 2017, Zheng et al. (2017) [215] identified the bias and introduced a stability parameter constraint \((z/L)_{lim}\) to prevent the land-atmosphere system from fully decoupling:
\[ (z/L)_{lim}=\frac{ln(\frac{z}{z_{0M}})}{2\alpha(1-\frac{z_{0M}}{z})} \]
Here \(z\) is the height, \(L\) is the Obukhov length, \(z_{0M}\) is the momentum roughness length, and \(\alpha = 5\).
The pertinent features of the GFS stable surface layer parameterization scheme are described in the appendix of Zheng et al. (2017) [215].
Recently, Han et al.(2021) [84] and Han et al.(2022) [82] introduced a few new features: