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.
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- Award ID(s):
- 1739705
- PAR ID:
- 10445738
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Atmospheres
- Volume:
- 124
- Issue:
- 23
- ISSN:
- 2169-897X
- Page Range / eLocation ID:
- p. 12631-12654
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract -
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
-
Abstract Accurate prediction of snow water equivalent (SWE) can be valuable for water resource managers. Recently, deep learning methods such as long short-term memory (LSTM) have exhibited high accuracy in simulating hydrologic variables and can integrate lagged observations to improve prediction, but their benefits were not clear for SWE simulations. Here we tested an LSTM network with data integration (DI) for SWE in the western United States to integrate 30-day-lagged or 7-day-lagged observations of either SWE or satellite-observed snow cover fraction (SCF) to improve future predictions. SCF proved beneficial only for shallow-snow sites during snowmelt, while lagged SWE integration significantly improved prediction accuracy for both shallow- and deep-snow sites. The median Nash–Sutcliffe model efficiency coefficient (NSE) in temporal testing improved from 0.92 to 0.97 with 30-day-lagged SWE integration, and root-mean-square error (RMSE) and the difference between estimated and observed peak SWE values
d maxwere reduced by 41% and 57%, respectively. DI effectively mitigated accumulated model and forcing errors that would otherwise be persistent. Moreover, by applying DI to different observations (30-day-lagged, 7-day-lagged), we revealed the spatial distribution of errors with different persistent lengths. For example, integrating 30-day-lagged SWE was ineffective for ephemeral snow sites in the southwestern United States, but significantly reduced monthly-scale biases for regions with stable seasonal snowpack such as high-elevation sites in California. These biases are likely attributable to large interannual variability in snowfall or site-specific snow redistribution patterns that can accumulate to impactful levels over time for nonephemeral sites. These results set up benchmark levels and provide guidance for future model improvement strategies. -
Abstract Water stored in mountain snowpacks (i.e., snow water equivalent, SWE) represents an important but poorly characterized component of the terrestrial water cycle. The Western United States snow reanalysis (WUS–SR) dataset is novel in its combination of spatial resolution (~500 m), spatial extent (31°–49° N; 102°–125° W), and temporal continuity (daily over 1985–2021). WUS–SR is generated using a Bayesian framework with model-based snow estimates updated through the assimilation of cloud-free Landsat fractional snow-covered area observations. Over the WUS, the peak SWE verification with independent
in situ measurements show correlation coefficient, mean difference (MD), and root mean squared difference (RMSD) of 0.77, −0.15 m, and 0.28 m, respectively. The effects of forest cover and Landsat image availability on peak SWE are assessed. WUS–SR peak SWE is well correlated (ranging from 0.75 to 0.91) against independent lidar-derived SWE taken near April 1st, with MD <0.15 m and RMSD <0.38 m. The dataset is useful for characterizing WUS mountain snow storage, and ultimately for improving snow-derived water resources management. -
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.