skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Implementation and Evaluation of a Unified Turbulence Parameterization Throughout the Canopy and Roughness Sublayer in Noah‐MP Snow Simulations
Abstract The Noah‐MP land surface model (LSM) relies on the Monin‐Obukhov (M‐O) Similarity Theory (MOST) to calculate land‐atmosphere exchanges of water, energy, and momentum fluxes. However, MOST flux‐profile relationships neglect canopy‐induced turbulence in the roughness sublayer (RSL) and parameterize within‐canopy turbulence in an ad hoc manner. We implement a new physics scheme (M‐O‐RSL) into Noah‐MP that explicitly parameterizes turbulence in RSL. We compare Noah‐MP simulations employing the M‐O‐RSL scheme (M‐O‐RSL simulations) and the default M‐O scheme (M‐O simulations) against observations obtained from 647 Snow Telemetry (SNOTEL) stations and two AmeriFlux stations in the western United States. M‐O‐RSL simulations of snow water equivalent (SWE) outperform M‐O simulations over 64% and 69% of SNOTEL sites in terms of root‐mean‐square‐error (RMSE) and correlation, respectively. The largest improvements in skill for M‐O‐RSL occur over closed shrubland sites, and the largest degradations in skill occur over deciduous broadleaf forest sites. Differences between M‐O and M‐O‐RSL simulated snowpack are primarily attributable to differences in aerodynamic conductance for heat underneath the canopy top, which modulates sensible heat flux. Differences between M‐O and M‐O‐RSL within‐canopy and below‐canopy sensible heat fluxes affect the amount of heat transported into snowpack and hence change snowmelt when temperatures are close to or above the melting point. The surface energy budget analysis over two AmeriFlux stations shows that differences between M‐O and M‐O‐RSL simulations can be smaller than other model biases (e.g., surface albedo). We intend for the M‐O‐RSL physics scheme to improve performance and uncertainty estimates in weather and hydrological applications that rely on Noah‐MP.  more » « less
Award ID(s):
1637686
PAR ID:
10360619
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Advances in Modeling Earth Systems
Volume:
13
Issue:
11
ISSN:
1942-2466
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract When compared with differences in snow accumulation predicted by widely used hydrological models, there is a much greater divergence among otherwise “good” models in their simulation of the snow ablation process. Here, we explore differences in the performance of the Variable Infiltration Capacity model (VIC), Noah land surface model with multiparameterization options (Noah-MP), the Catchment model, and the third-generation Simplified Simple Biosphere model (SiB3) in their ability to reproduce observed snow water equivalent (SWE) during the ablation season at 10 Snowpack Telemetry (SNOTEL) stations over 1992–2012. During the ablation period, net radiation generally has stronger correlations with observed melt rates than does air temperature. Average ablation rates tend to be higher (in both model predictions and observations) at stations with a large accumulation of SWE. The differences in the dates of last snow between models and observations range from several days to approximately a month (on average 5.1 days earlier than in observations). If the surface cover in the models is changed from observed vegetation to bare soil in all of the models, only the melt rate of the VIC model increases. The differences in responses of models to canopy removal are directly related to snowpack energy inputs, which are further affected by different algorithms for surface albedo and energy allocation across the models. We also find that the melt rates become higher in VIC and lower in Noah-MP if the shrub/grass present at the observation sites is switched to trees. 
    more » « less
  2. Abstract Land surface processes are vital to the performance of regional climate models in dynamic downscaling application. In this study, we investigate the sensitivity of the simulation by using the weather research and forecasting (WRF) model at 10-km resolution to the land surface schemes over Central Asia. The WRF model was run for 19 summers from 2000 to 2018 configured with four different land surface schemes including CLM4, Noah-MP, Pleim-Xiu and SSiB, hereafter referred as Exp-CLM4, Exp-Noah-MP, Exp-PX and Exp-SSiB respectively. The initial and boundary conditions for the WRF model simulations were provided by the National Centers for Environmental Prediction Final (NCEP-FNL) Operational Global Analysis data. The ERA-Interim reanalysis (ERAI), the GHCN-CAMS and the CRU gridded data were used to comprehensively evaluate the WRF simulations. Compared with the reanalysis and observational data, the WRF model can reasonably reproduce the spatial patterns of summer mean 2-m temperature, precipitation, and large- scale atmospheric circulation. The simulations, however, are sensitive to the option of land surface scheme. The performance of Exp-CLM4 and Exp-SSiB are better than that of Exp-Noah-MP and Exp-PX assessed by Multivariable Integrated Evaluation (MVIE) method. To comprehensively understand the dynamic and physical mechanisms for the WRF model’s sensitivity to land surface schemes, the differences in the surface energy balance between Ave-CLM4-SSiB (the ensemble average of Exp-CLM4 and Exp-SSiB) and Ave-NoanMP-PX (the ensemble average of Exp-Noah-MP and Exp-PX) are analyzed in detail. The results demonstrate that the sensible and latent heat fluxes are respectively lower by 30.42 W·m−2and higher by 14.86 W·m−2in Ave-CLM4-SSiB than that in Ave-NoahMP-PX. As a result, large differences in geopotential height occur over the simulation domain. The simulated wind fields are subsequently influenced by the geostrophic adjustment process, thus the simulations of 2-m temperature, surface skin temperature and precipitation are respectively lower by about 2.08 ℃, 2.23 ℃ and 18.56 mm·month−1in Ave-CLM4-SSiB than that in Ave-NoahMP-PX over Central Asia continent. 
    more » « less
  3. Abstract Accurate precipitation estimates are critical to simulating seasonal snowpack evolution. We conduct and evaluate high‐resolution (4‐km) snowpack simulations over the western United States (WUS) mountains in Water Year 2013 using the Noah with multi‐parameterization (Noah‐MP) land surface model driven by precipitation forcing from convection‐permitting (4‐km) Weather Research and Forecasting (WRF) modeling and four widely used high‐resolution datasets that are derived from statistical interpolation based on in situ measurements. Substantial differences in the precipitation amount among these five datasets, particularly over the western and northern portions of WUS mountains, significantly affect simulated snow water equivalent (SWE) and snow depth (SD) but have relatively limited effects on snow cover fraction (SCF) and surface albedo. WRF generally captures observed precipitation patterns and results in an overall best‐performed SWE and SD in the western and northern portions of WUS mountains, where the statistically interpolated datasets lead to underpredicted precipitation, SWE, and SD. Over the interior WUS mountains, all the datasets consistently underestimate precipitation, causing significant negative biases in SWE and SD, among which the results driven by the WRF precipitation show an average performance. Further analysis reveals systematic positive biases in SCF and surface albedo across the WUS mountains, with similar bias patterns and magnitudes for simulations driven by different precipitation datasets, suggesting an urgent need to improve the Noah‐MP snowpack physics. This study highlights that convection‐permitting modeling with proper configurations can have added values in providing decent precipitation for high‐resolution snowpack simulations over the WUS mountains in a typical ENSO‐neutral year, particularly over observation‐scarce regions. 
    more » « less
  4. The influence of the Unified Noah and Noah-MP land surface models (LSMs) on the evolution of cumulus clouds reaching convective initiation (CI) is assessed using infrared brightness temperatures (BT) from GOES-16. Cloud properties from individual cloud objects are examined using output from high-resolution (500 m horizontal grid spacing) model simulations. Cloud objects are tracked over time and related to observed clouds reaching CI to examine differences in cloud extent, longevity, and growth rate. The results demonstrate that differences in assumed surface properties can lead to large discrepancies in the net surface radiative budget, particularly in the sensible and latent heating components where differences exceed 40 W m−2. These differences lead to changes in the local mesoscale circulation patterns that are more pronounced near the edges of forested and grassland boundaries where lower-level convergence is stronger. Higher sensible heating from the Noah-MP LSM produced growth of CI clouds earlier and with increased longevity, which was closer to the timing and growth observed from GOES-16. The increased cloud growth in the Noah-MP experiment results from stronger and deeper updrafts, which lofts more cloud water into the upper levels of the troposphere. The weaker updrafts from the Noah LSM experiment results in shallower convection after CI is detected due to slower growth rates. The differences in cloud properties and growth are directly related to the land surfaces they develop above and point to the importance of accurately representing land properties and radiative characteristics when simulating convection in numerical weather prediction models. 
    more » « less
  5. Abstract The influence of the Unified Noah and Noah‐MP land surface models (LSMs) on the evolution of cumulus clouds reaching convective initiation (CI) is assessed using infrared brightness temperatures (BT) from GOES‐16. Cloud properties from individual cloud objects are examined using output from high‐resolution (500 m horizontal grid spacing) model simulations. Cloud objects are tracked over time and related to observed clouds reaching CI to examine differences in cloud extent, longevity, and growth rate. The results demonstrate that differences in assumed surface properties can lead to large discrepancies in the net surface radiative budget, particularly in the sensible and latent heating components where differences exceed 40 W m−2. These differences lead to changes in the local mesoscale circulation patterns that are more pronounced near the edges of forested and grassland boundaries where lower‐level convergence is stronger. Higher sensible heating from the Noah‐MP LSM produced growth of CI clouds earlier and with increased longevity, which was closer to the timing and growth observed from GOES‐16. The increased cloud growth in the Noah‐MP experiment results from stronger and deeper updrafts, which lofts more cloud water into the upper levels of the troposphere. The weaker updrafts from the Noah LSM experiment results in shallower convection after CI is detected due to slower growth rates. The differences in cloud properties and growth are directly related to the land surfaces they develop above and point to the importance of accurately representing land properties and radiative characteristics when simulating convection in numerical weather prediction models. 
    more » « less