Improving snow albedo modeling in E3SM land model (version 2.0) and assessing its impacts on snow and surface fluxes over the Tibetan Plateau

Abstract

Abstract. With the highest albedo of the land surface, snow plays a vital role in Earth’s surface energy and water cycles. Snow albedo is greatly affected by snow grain properties (e.g., size and shape) and light absorbing particles (LAPs) such as black carbon (BC) and dust. The mixing state of LAPs in snow also has large impacts on LAP-induced snow albedo reduction and surface radiative forcing (RF). However, most land surface models assume that snow grain shape is spherical and LAPs are externally mixed with the snow grains. This study improves the snow radiative transfer model in the land model (ELM v2.0) of the Energy Exascale Earth System Model version 2.0 (E3SM v2.0) by considering non-spherical snow grain shapes (i.e., spheroid, hexagonal plate and Koch snowflake) and internal mixing of dust-snow and systematically evaluates the impacts on surface energy and water balances over the Tibetan Plateau (TP). A series of ELM simulations with different treatments of snow grain shape, mixing state of BC-snow and dust-snow, and sub-grid topographic effects (TOP) on solar radiation are performed. Compared with two remote sensing snow products derived from the Moderate Resolution Imaging Spectroradiometer, the control ELM simulation with the default settings of spherical snow grain shape, internal mixing of BC-snow, external mixing of dust-snow and without TOP can capture the overall snow distribution reasonably. The estimated LAP-induced RF ranging from 0 to 21.9 W/m2 with the area-weighted average value of 1.3 W/m2 is comparable to reported values. Focusing on the snow-related processes and surface energy and water balances, Koch snowflake shape, among other non-spherical shapes, shows the largest difference from spherical shape in spring. The impacts of the mixing state of LAP-snow are smaller than the shape effects and depend on snow grain shape. Compared to external mixing, internal mixing of LAP-snow can lead to larger snow albedo reduction and snowmelt, which further affect surface energy and water cycles. Compared to the control simulation, the individual contributions of non-spherical snow shape, mixing state of LAP-snow, and local topography to the change of snow and surface fluxes have different signs and magnitudes, and their combined effects may be negative or positive due to complex and non-linear interactions among the factors. Overall, the changes of net solar radiation in spring due to individual and combined effects range from -28.6 to 16.9 W/m2 and -29.7 to 12.2 W/m2, respectively. This study advances understanding of the role of snow grain shape and mixing state of LAP-snow in land surface processes and offers guidance for improving snow simulations and RF estimates in Earth system models under climate change.

Wildfire impacts on western United States snowpacks

“…Earth system models have shown notable sensitivity to changes in snow-albedo (Letcher and Minder, 2015). However, radiative transfer models still miss the mark on incorporating LAP in snow radiative forcing (Hao et al, 2022). Gleason and Nolin (2016) developed an albedo decay function that accounts for BC.…”

“…Ultimately, we need to bring together multiscale models to incorporate these changes into models spanning the atmosphere to bedrock to fully address this issue (Siirila-Woodburn et al, 2021). If net shortwave and topography are the biggest drivers of changes in SWE (Maina and Siirila-Woodburn, 2020), then incorporating canopy changes and snow albedo is imperative to reduce uncertainty in land surface (Letcher and Minder, 2015;Hao et al, 2022). Additional physically-based model studies examining the direct impacts of wildfire-induced landscape changes on snowpack dynamics are needed to supplement growing observational studies.…”

Evaluation of snow processes over the Western United States in E3SM land model

“…First, the original SNICAR model has been replaced by a hybrid model (SNICAR-AD) of SNICAR and delta-Eddington adding-doubling radiative transfer solver, which corrects the snow albedo bias for large SZAs and can better represent the shortwave radiative properties of snow (Dang et al, 2019). Second, compared to only external mixing in CLM4.5, both external mixing and internal mixing of hydrophilic BC-snow and dust-snow are now represented in ELM (Hao et al, 2022a;Wang et al, 2020). Third, the direct and diffuse irradiance under different atmospheric profiles and their dependence on SZA are included (Hao et al, 2022a).…”

“…Second, compared to only external mixing in CLM4.5, both external mixing and internal mixing of hydrophilic BC-snow and dust-snow are now represented in ELM (Hao et al, 2022a;Wang et al, 2020). Third, the direct and diffuse irradiance under different atmospheric profiles and their dependence on SZA are included (Hao et al, 2022a).…”

“…The Energy Exascale Earth System Model (E3SM) Land Model (ELM) (Leung et al, 2020) includes a multi-layer snow scheme to simulate the prognostic snow processes such as snow accumulation, snow interception, snow compaction, and snow melt. Recently, the snow albedo model in ELM was improved to include new radiative transfer solvers with improved accuracy (Dang et al, 2019), add non-spherical snow grain shape (Hao et al, 2022a), account for the internal mixing of light-absorbing particles (LAPs) with snow (Böttcher et al, 2014); Hao et al, 2022a), and incorporate new parameterizations to account for the subgrid topographic effects on solar radiation (Hao et al, 2021;Hao et al, 2022b) (see Section 2.1 for details). With these enhancements and improvements, ELM may skillfully simulate snow dynamics at a regional scale (e.g., WUS).…”