Snow depth on the interior of Tibetan Plateau (TP) in state-of-the-art reanalysis products is almost an order of magnitude higher than observed. This huge bias stems primarily from excessive snowfall, but inappropriate process representation of shallow snow also causes excessive snow depth and snow cover. This study investigated the issue with respect to the parameterization of fresh snow albedo. The characteristics of TP snowfall were investigated using ground truth data. Snow in the interior of the TP is usually only some centimeters in depth. The albedo of fresh snow depends on snow depth, and is frequently less than 0.4. Such low albedo values contrast with the high values (~0.8) used in the existing snow schemes of land surface models. The SNICAR radiative transfer model can reproduce the observations that fresh shallow snow has a low albedo value, based on which a fresh snow albedo scheme was derived in this study. Finally, the impact of the fresh snow albedo on snow ablation was examined at 45 meteorological stations on TP using the land surface model Noah-MP which incorporated the new scheme. Allowing albedo to change with snow depth can produce quite realistic snow depths compared with observations. In contrast, the typically assumed fresh snow albedo of 0.82 leads to too large snow depths in the snow ablation period averaged across 45 stations. The shallow snow transparency impact on snow ablation is therefore particularly important in the TP interior, where snow is rather thin and radiation is strong.
The average elevation of the Tibetan Plateau (TP), the “roof of the world,” is more than 4000 m above the sea level. Strong uplifting of the TP since the Neogene has produced a natural environment distinctive from other regions of the world. Characteristics include large seasonal and very large diurnal variations of temperature, intense solar radiation, numerous rivers and lakes, extensive glacier and permafrost, and abundant biodiversity (Zhou et al. 2000). The TP has experienced warming at a higher rate than the global mean especially during the 1980s–2000s (Liu and Chen 2000), and is expected to warm faster in the future climates as well (IPCC 2013). This makes the cryosphere and ecosystem of the TP particularly vulnerable, but until now the mechanisms of glacier retreat, permafrost degradation, lake expansion and/or shrinkage, vegetation succession along with the feedbacks to the atmosphere are far from understood (Wu et al. 2016; Sun et al. 2018; Yang et al. 2018). More comprehensive station observations are expected to deepen our understanding of the land–atmosphere interactions on TP, where meteorological stations are very sparse.
In this unique region, snow plays a crucial role in the water and energy cycles, and can have profound feedbacks to climate (Zhang 2005). Snow modifies the ground energy balance by its high albedo and also provides strong insulation due to its low thermal conductivity. The spatial and temporal variation of TP snow cover at various scales affects both the local weather status and the Asian summer monsoon (Blanford 1884; Xiao and Duan 2016). Snowmelt in spring also serves as an important water supply to some large river systems on the TP (Han et al. 2019). Snow is therefore a key component in the Tibetan climate system.
The characteristics of snow on the TP is quite different from that in the circum-Arctic region (Wang et al. 2018). Due to the low latitude of the TP, the solar radiation is strong throughout the year, and snow ablates quickly on clear and dry winter and spring days. The snow is rather heterogeneous across the TP, thick snow from tens of centimeters to several meters on the mountains and shallow snow of only a few centimeters across the interior (Qin et al. 2006). In addition, where the wind is strong, blowing snow both leads to patchy snowpack and bare ground, and reduces snow accumulation by sublimation (Shao et al. 2017).
State-of-the-art reanalysis datasets can produce 10 times the observed snow depth in the TP interior (Orsolini et al. 2019). The huge bias stems primarily from too much snowfall in the reanalysis model (Orsolini et al. 2019). However, various snow processes that might be important may also be missing or misrepresented. These include blowing snow sublimation (Orsolini et al. 2019), snow patch effects, and dust deposition (Qian et al. 2011). Current snow schemes used in land surface and reanalysis models have emphasized development of increasing number of snow layers, dynamic evolution of snow density and thermal conductivity (Wang et al. 2016). These are important processes to simulate snow depth and the insulation effect on the ground (Wang et al. 2016). Another key component of the snow scheme is the parameterization of snow surface albedo, which greatly determines the ground surface energy balance. Although many snow schemes assume fresh snow albedo as high as ~0.8, Giddings and LaChapelle (1961) and McKay and Gray (1981) showed that the combined albedo of shallow snow over soil can be far lower than 0.8, and only achieves the value of 0.8 after accumulation of several centimeters of snow depth.
Snow albedo is described by various models and schemes of varying complexity. Table 1 summarizes the current fresh snow albedo schemes used by state-of-the-art land surface models. Almost all the land surface models adopt simplified snow albedo schemes, which assume a high fixed value as the fresh snow albedo (0.78 for HTESSEL, 0.8 for CoLM, 0.84 for CLASS, and 0.8 for BATS in Noah-MP, etc.) and use a snow aging factor to describe the decrease of snow albedo with time. The CLASS scheme in Noah-MP includes the snow depth in fresh snow albedo calculation only when snowfall amount is less than 1 mm and sets a lower bound of snow albedo as 0.55 (see details given in the appendix). The CLM uses the Snow and Ice Aerosol Radiation model (SNICAR) to calculate snow albedo. This model calculates the radiation transfer in snow with the consideration of snow grain size and snow depth (Flanner and Zender 2006), and outputs the spectral reflectance and broadband albedo. SNICAR is excellent for investigating the impact of light-absorbing impurities on snow aging and has been widely applied in remote sensing and point-scale albedo simulations (Singh et al. 2018; He et al. 2018; Lawrence et al. 2011; Flanner et al. 2007; Flanner and Zender 2005). However, SNICAR is computationally expensive when it comes to high-resolution or long time series of climate simulations. Hence, a simplified version of SNICAR is preferable for regional land surface modeling.
These listed fixed high fresh snow albedo parameters (Table 1) are designed for thick snow, but the snow cover over most of TP is quite shallow which can come and go several times a year and thus it is usually difficult to explicitly distinguish the accumulation and ablation seasons. Snow of only a few centimeters can be somewhat transparent, so the albedo of fresh shallow snow is prone to be affected by the ground beneath. While SNICAR can describe the effect of bare ground on the snow albedo, other schemes with fixed fresh snow albedo values (~0.78–0.9) run the risk of being too high for the shallow snow. Thus, an important question here is how the albedo of fresh shallow snow differs from that of fresh thick snow, and how the albedo difference affects the snow ablation on the interior of the TP.
The aim of this study is to investigate the characteristics of fresh snow albedo on the TP interior and its effect on snow ablation. We used observations of upward/downward shortwave radiation and snow depth at four experimental stations (section 2) to investigate the characteristics of snow depth and albedo (section 3). We demonstrated that the SNICAR model can reproduce the observed variability of snow albedo in shallow snow cases and then developed a new fresh snow albedo scheme which was fitted with the SNICAR simulations (section 4). Finally, we tested the new scheme in Noah-MP and performed validation with the snow depth observed at 45 meteorological stations on the TP (section 5).
2. Data and models
a. Observations at experimental stations
The experimental station observations were used to analyze the characteristics of snow depth and albedo in the TP. Ground truth from station-observed snow accumulation and ablation with simultaneous radiation and snow depth data on the TP is quite few. We used all data that are accessible, which means that the observation periods are not the same across the stations. These stations include ANNI-AWS (1 October 2002–14 August 2004 hourly), D105 (1 October 2002–28 November 2004, hourly) obtained through CEOP (Coordinated Enhanced Observing Period) Asia–Australia monsoon project (CAMP) on the TP (Koike 2004; Ma et al. 2005), Tanggula Mountain Station (1 January 2012–31 December 2012, half-hourly) (Zhao et al. 2017; Zhao and Sheng 2019), and Dadongshu-Yakou Station (17 December 2013–31 August 2016 hourly) (Che et al. 2019; Li et al. 2019). Their locations are shown in Fig. 1. The altitudes of the four stations are higher than 4000 m, with a cold and dry climate that is typical on the TP. Among them, three stations (ANNI-AWS, D105, Tanggula) are located in the interior of the TP and one (Dadongshu-Yakou) is located in the Qilian Mountains, Northeast TP. The mountain station experiences more snowfall events than the three interior stations, but they share similar albedo features when snow is shallow.
The measurements are described in detail in the corresponding references. The four-component radiation sensor at Dadongshu-Yakou and Tanggula station is a CNR4. The shortwave radiation sensor at D105 and ANNI station is a CM21 and the longwave radiation sensor is a precision infrared radiometer. Radiation data were recorded every 30 min or 1 h by the radiometers installed at a height of 1 m at the Dadongshu-Yakou station, 2 m at the Tanggula station, and at 1.28 m for upward shortwave/longwave radiation and 1.58 m for downward shortwave/longwave radiation at both the D105 and ANNI station.
All the measurements have been quality controlled as far as possible. Snow depth can be difficult to reliably ascertain, especially in harsh and rugged terrain and therefore we made cross checks of data quality. Observed meteorological and snow depth data with good quality are expected to meet all of the following requirements: 1) the upward longwave radiation must be less than 315.6 W m−2, that is, the ground surface temperature cannot be above 0°C when snow exists; 2) air temperature must be less than 5°C during snowfall, since snowfall can happen on the TP when air temperature is above 0°C; 3) measured snow depth should be greater than 1 cm, which is the measurement uncertainty of snow depth; and 4) observed snow depth increases when snowfall is recorded–this requirement was not applied to Tanggula station because precipitation events were not recorded during all the measurement period. Observed snow albedo was calculated as the ratio of upward to downward shortwave radiation following Che et al. (2019). To further eliminate the impact of solar angle on albedo, we took averages between 1000 and 1400 local time of albedo, snow depth, air temperature, and precipitation, as the “daily” values for the four experimental sites in the statistical analysis.
b. Observations at routine meteorological stations
Snow depth is measured at routine meteorological stations of China Meteorological Administration (CMA) on the TP. The station data were used to validate the impact of fresh snow albedo on the snow ablation simulation.
Daily mean datasets of air temperature, precipitation, relative humidity, pressure and snow depth observed during 2012–14 were used. The observed precipitation was prone to be underestimated due to the wind. We therefore performed bias correction according to the methodology described in Ye et al. (2004). The daily snow depth was recorded at 0800 local time of each day. We further quality controlled the data to eliminate stations with questionable data. A station was removed if there were frequent missing data or there were discrepancies in the time occurrences between the observed precipitation and snow depth. After quality control, there are 45 meteorological stations available (Fig. 1).
The forcing data to drive the land surface model Noah-MP (see section 2c) were made available through fusing the daily mean of the station data with the diurnal cycle of the China Meteorological Forcing Dataset (CMFD). The CMFD at 0.1° spatial resolution and at 3 h temporal resolution was developed by the Institute of Tibetan Plateau Research, Chinese Academy of Sciences (Yang et al. 2010; He et al. 2020). CMFD includes precipitation, downward shortwave radiation, downward longwave radiation, 2-m air temperature, specific humidity, wind speed, and surface pressure. Missing air temperatures in the meteorological observation were replaced with the CMFD grid’s air temperature and the air temperature lapse rate (6.5° km−1) was used to correct the impact of altitude difference between CMFD grids and stations. Missing pressures were replaced with the observed pressure of the previous day. Missing wind speed and relative humidities were replaced with the grid values of CMFD. The meteorological specific humidity was then calculated with the corrected air temperature, pressure and relative humidity. Downward shortwave and longwave radiation data were not provided by CMA and they were directly extracted from CMFD.
We used two models: SNICAR was used to demonstrate its capability for reproducing the observed characteristics of fresh snow albedo, and the Noah-MP land surface model was used to verify the impact of fresh snow albedo on the TP snow ablation.
SNICAR is a physically based radiative transfer model for snow, which incorporates a two-stream radiative transfer solution from Toon et al. (1989). Following the adapted calculation methodology from Wiscombe (1980), optical properties from Mie theory are calculated for each snow layer and spectral band (Flanner and Zender 2006; Flanner et al. 2007). With the inclusion of extinction optical depth, SNICAR can consider the impact of snow depth on albedo. SNICAR simulates solar absorption both in the snow layers and the underlying substrate. Snow aging is calculated with the ice effective grain size, which conserves the specific surface area (Grenfell and Warren 1999). The change in effective grain size is caused both by the dry and wet snow processes, including liquid-water-induced metamorphism, refreezing of liquid water, the addition of freshly fallen snow, and dry snow metamorphism. SNICAR can reproduce some typically observed phenomena like rapid snow aging from the combination of warm snow, large temperature gradient, and low density.
Noah-MP is a unique multiparameterization land surface model that allows the user to combine different schemes for simulating individual land processes (Niu et al. 2011). The snow model in Noah-MP is of up to three layers, depending on the snow depth. Snow layer temperature, snow height, snow water equivalent, and snow cover fraction have to meet the requirement of energy and mass balance. Snow layer division and combination, liquid water storage and melt/refreeze capability, and snow interception by vegetation are all described in the model.
There are two options for snow albedo scheme in Noah-MP: CLASS and BATS. The equations of the two snow albedo schemes are listed in the appendix. Noah-MP calculates albedo for both the visible and near-infrared wave bands under both direct and diffuse radiation. The snow albedo scheme of CLASS for the two wave bands under both direct and diffuse radiation is the same. The fresh snow albedo of CLASS is 0.84 if the new fallen snow mass >1 mm. Below the critical value the fresh snow albedo is an average of old snow albedo and 0.84, weighted by the old and new snow mass. Snow ages as an exponential function with time and the minimum value of snow albedo is 0.55. In BATS, the fresh snow albedo for the visible band is 0.95 and the near-infrared band is 0.65. The snow albedo under direct radiation is a linear function of the snow albedo under diffuse radiation, and snow aging is a function of time, ground temperature and snow mass. Neither of the schemes strictly describes the impact of snow transparency on fresh snow albedo and neither of the schemes includes the calculation of snow grain size.
3. Observed results
a. Observed snow depth and fresh snow albedo
The histograms of snow depth for the four experimental stations shown in Fig. 1 are presented in Fig. 2. The fraction of winter period (DJF) with snow depth >1 cm is labeled, and a few days with snow cover for D105 station occurred from March to June. The snow is quite shallow, except at Dadongshu-Yakou, which is in the Qilian Mountains. The percentage of days with snow depth <3 cm during each station’s research period exceeds 50% for ANNI-AWS, 100% for D105%, and 87% for Tanggula. Even for Dadongshu-Yakou, days with snow depth <3 cm still amount to 30% of the total snow-cover days. These observed snow depths are in agreement with Qin et al. (2006), who noted that the snow on the interior of the TP is quite shallow, usually less than 3 cm.
Within the hourly or half-hourly observations averaged over 1000–1400 local time, a snowfall case (a day with snowfall) should meet the following requirements: 1) the averaged air temperature is <5°C, 2) the precipitation amount >0, and 3) the averaged snow depth increases compared with the previous day. To investigate the albedo characteristics of fresh snow on shallow snowpack or bare ground, we show in Fig. 3 the surface albedo change after all snowfall events for the four experimental stations where surface albedo was <0.3 before snowfall. The snowfall events are represented by spikes. Current snow schemes set a fixed fresh snow albedo of 0.7–0.9, which is in agreement with the black spikes plotted in Fig. 3 describing an albedo increase to more than 0.7. However, most of the changes in albedo are much less, ranging from 0.1 to 0.4, as indicated by the colored spikes. In total across the four stations, 78 of the 94 albedo changes were gradual increases in albedo rather than an immediate jump to 0.7 or more. For ANNI-AWS, only 1 of the 5 cases shows a sharp increase of albedo with the other 4 cases presenting gradual increases. For D105, there is only one snowfall case, which shows a rather small increase (<0.1) of fresh snow albedo. For Tanggula, only 2 of the 25 snowfall cases show a sharp albedo increase and all the remaining are with a gradual change. For the 63 snowfall cases at Dadongshu-Yakou station, 13 cases are with a sharp albedo increase and 50 cases are with a gradual change. In general, surface albedo gradual increases after snowfall on the TP rather than undergoing an abrupt change to high albedo.
The distribution of fresh snow albedo is shown in Fig. 4. Over 30% of snowfall cases produce a fresh snow albedo below 0.5. Clearly, the fresh snow albedo on the TP is far below the default value (~0.8) used in most of the current snow schemes.
b. Observed relationship between fresh snow albedo and snow depth
In this section we explored the factors that affect the fresh snow albedo. In addition to the snow depth mentioned in Giddings and LaChapelle (1961), air temperature and water vapor also drive the growth of snow grain (Rasmus 2005) and affect albedo. We therefore looked into the correlation coefficients between fresh snow albedo and each of air temperature, relative humidity and snow depth, which are shown in Table 2. The correlation coefficient at ANNI-AWS, and D105 is not defined, as the sample size is too small (see Fig. 2). Significant correlations at both remaining stations only occur between fresh snow albedo and snow depth. For Dadongshu-Yakou and Tanggula, the correlation coefficient is 0.61 and 0.79, respectively. As might be expected, relative humidity is a significant factor at the high snow-accumulation Dadongshu-Yakou site, and surface temperature at the low snow-accumulation Tanggula station. Snow depth is therefore the dominant factor determining fresh snow albedo at both high and low accumulation sites, and the deeper the snow depth (<4 cm), the higher the fresh snow albedo.
Figure 5 shows the boxplot of fresh snow albedo varying with respect to snow depth. Again, the snowfall cases with snow depth >1 cm are very limited for ANNI-AWS and D105, and thus they are not shown in Fig. 5. The boxplot for each station shows the range of fresh snow albedo within the individual snow-depth categories. When the snow depth is less than 4 cm, the fresh snow albedo of Dadongshu-Yakou in more than half of the snowfall cases remains below 0.4, although it occasionally reaches 0.8 or more; when the snow is thicker than 4 cm, the albedo in more than half of the snowfall cases is higher than 0.7 and occasionally drops below 0.5. The fresh snow albedo at Tanggula shows an obvious increase with respect to snow depth when snow depth is less than about 5 cm. The results of the TP stations are consistent with the much earlier finding in the Arctic, that several centimeters of snow are required for the albedo of freshly fallen snow to reach an asymptotic high value (Giddings and LaChapelle 1961).
The existence of high albedo values in thin snow cases or low albedo values in thick snow cases indicates that second-order factors can also affect the fresh snow albedo. For example, increased cloud fraction will change the solar spectral distribution (i.e., after passing through clouds, the energy ratio of visible band to near-infrared band increases) and therefore enhance snow surface albedo. However, although important, investigation on these second-order factors is beyond the scope of this study and is expected to be addressed in future.
4. Simulated snow albedo–snow depth relationship
First, SNICAR is demonstrated to be capable of reproducing the observed dynamics at the four experimental stations. Then, the SNICAR-simulated albedo values are used to support the development of a simple fresh snow albedo scheme.
a. SNICAR-simulated dependence of fresh snow albedo on snow depth and grain size
Snow grain size is an important snowpack variable influencing the albedo is. It is the effective snow grain radius or diameter if the ice crystal is nonspherical (Grenfell and Warren 1999) that is represented by a collection of spheres with the same volume-to-area ratio. Considering snow albedo is greatly affected by the snow grain size as well as snow depth, we run SNICAR with a group of snow depth values ranging from 0.1 to 20 cm and snow grain diameters ranging from 20 to 1000 μm. The results are shown by the symbols in Fig. 6. The observed increase of fresh snow albedo with respect to snow depth shown in Fig. 5 can be approximately reproduced by SNICAR.
b. Development of a fresh snow albedo scheme
Weather forecasting and climate modeling usually require time-saving parameterization schemes. We therefore made efforts to simplify SNICAR and derived a parameterization which is an approximate relationship between fresh snow albedo, snow grain size and snow depth (Fig. 6).
The scheme is a heuristic one, which consists of two terms: one is the back scattering when solar radiation penetrates the snow, and the other is the forward scattering for the reflected solar radiation from the ground surface. The scheme is given below:
where albsnow is the fresh snow albedo, z is the snow depth (cm), k is the extinction coefficient (cm−1) that describes the absorption and scattering capability of the snowpack.
Considering two extreme conditions—1) albsnow = alb∞ if z approaches infinite snow depth, where alb∞ is the albedo of fresh semi-infinite snow, and 2) albsnow = albground if z = 0, where albground is the background albedo (i.e., the albedo of the old snow or the snow-free ground before snowfall)—we can determine β1 = alb∞, β2 = albground. Therefore,
alb∞ and the extinction coefficient k are dependent on snow grain size. Based on the SNICAR simulated albedo values, alb∞ and k are parameterized as follows:
where Dsnw is the snow grain size in μm. Parameters a, b, c, and d are the regression coefficients and were fitted by the SNICAR simulations. The fitted values are presented in Table 3.
Other processes like blowing snow and dust deposition could also decrease the snow albedo (Rasmus 2005), and there is space to work further on the current scheme.
5. Impact of depth-dependent fresh snow albedo on snow ablation
We tested the impact of fresh snow albedo on the Tibetan snow processes using the Noah-MP (Niu et al. 2011). The new fresh snow albedo scheme derived in section 4b considers the impact of snow grain size, but the original snow schemes in Noah-MP (BATS and CLASS) do not provide this parameter. To demonstrate the impact of fresh snow albedo on snow ablation, we introduced a new set of snow albedo scheme in this study, which consists of the derived new fresh snow albedo scheme, the snow grain size scheme from Anderson (1976) as function of dewpoint temperature and the snow aging scheme from Baker et al. (1990) and Verseghy (1991). The details of snow grain size and snow aging schemes used in this study are listed in the appendix.
The model was run with this introduced snow albedo scheme which includes snow aging and grain size calculation. We also made a simulation with a fixed fresh snow albedo. The simulations were conducted at the 45 TP meteorological stations as shown in Fig. 1 over the period of 2012–14. Because of uncertainties in measuring snowfall and snow depth, the averaged snow depth across the 45 stations was used to validate the simulations. Since the observed precipitation does not distinguish snowfall and rainfall, a precipitation partition scheme based on air temperature, relative humidity and elevation (Ding et al. 2014, 2017) was used in this study.
The new fresh snow albedo scheme exerts a considerable impact on the snow depth on the TP. The simulated time series of daily mean snow depth from 2012 to 2014 was validated with station-averaged observation in Fig. 7. The simulation with the new fresh snow albedo fitted quite well with the observation in both the 2012–13 snow season when there are many snowfall events and the 2013–14 snow season when there are fewer snowfall events. In contrast, the simulation with the fixed fresh snow albedo of 0.82 overestimated snow during the snow ablation period in both the 2012–13 snow season and the 2013–14 snow season based on the averaged snow depth across 45 stations. The new fresh snow albedo is clearly a controlling factor in the snow ablation simulation.
Nevertheless, we found that the difference between the simulations in winter (DJF) is much less than the difference in spring (MAM). We analyzed the energy budget process in winter and in spring (not shown). The results showed that snow sublimation instead of melting dominates the snow ablation in winter. The absorbed solar radiation is used to warm the snow surface but the surface temperature still cannot reach the melting point, and thus more solar radiation absorbed by the surface in the new scheme is used to increase sensible heat flux and upward longwave radiation. Therefore, the snow depth between the two simulations is similar to each other in winter. However, in spring, the snow surface temperature increases to the melting point and thus the more solar radiation absorbed by the surface in the new scheme is mainly used for snow melting, causing large difference in snow depth between the two simulations.
We also compared the simulated results using the two original snow albedo schemes in Noah-MP with the observed snow depth, as presented in Fig. 8. Although the snow aging process is more complex in BATS, the simulated snow depth is very similar to the one simulated with the fixed fresh snow albedo in Fig. 7. The CLASS scheme outperforms BATS in the simulation because CLASS calculates the fresh snow albedo as the weighted average of the old snow albedo at the previous step and 0.84 when the new snowfall mass is <1 mm. This again confirms that the ground visibility impacts fresh snow albedo and impacts much the snow ablation on TP. However, since CLASS sets a lower bound of snow albedo as 0.55, it still significantly overestimates the snow depth during the snow ablation period based on the averaged results across 45 stations.
6. Conclusions and discussion
This study investigated the characteristics of fresh snow albedo in the TP interior and its effect on the snow ablation. According to experimental station data, snow is quite thin on the TP interior and the albedo of freshly fallen snow cannot reach 0.7 unless the snow depth reaches around 4 cm. This relationship can be reproduced by the physics-based radiative transfer model SNICAR. However, this fact is not well supported in parameterizations of snow albedo used in most of current land surface models, in which a fixed high value for fresh snow albedo is usually assumed regardless of the snow depth. Therefore, a new scheme was developed to account for the relationship between fresh snow albedo, snow grain size and snow depth based on SNICAR simulations.
The impact of the fresh snow albedo on the snow ablation was then investigated with the land surface model Noah-MP over 45 TP meteorological stations, using either a fixed high value of fresh snow albedo or the new fresh snow albedo scheme, or each of two original snow albedo schemes in Noah-MP. The results are validated with station observed daily snow depth. The modeling revealed that the fresh snow albedo that considers the impact of snow depth is a controlling factor in the snow ablation simulation. The simulated TP station snow depth corresponded well with the observation when the new fresh snow albedo scheme was applied. The simulation with the fixed fresh snow albedo, however, simulated too much snow in the snow ablation period based on the averaged snow depth across 45 TP stations. The simulations with both the original snow albedo schemes in Noah-MP are poorer than the one with the new fresh snow albedo scheme. This confirmed the importance of considering snow depth impact on fresh snow albedo.
Other factors in addition to snow depth and grain size can affect fresh snow albedo, too. Cold weather favors high albedo snow (Rasmus 2005; Armstrong and Brun 2008), whereas warm and moist air tends to lower the albedo due to the larger size of snow grain (Rasmus 2005). Clouds also produce effects on snow albedo by modifying the spectrum of solar radiation (i.e., cloud absorbs less visible energy than near-infrared energy), which causes snow broadband albedo higher (Salby 1996). Impurities such as black carbon and dust can greatly reduce the snow albedo in the marginal area of the TP. All these factors likely contribute to scatter in the observed relationship between fresh snow albedo and snow depth (see Fig. 5), and they should be accounted for in follow-on fresh snow albedo scheme development.
In summary, this study revealed the fresh snow albedo as a key factor in snow ablation on the interior of the TP, where the snow is rather thin and the radiation is strong. The absence of the dependence of fresh snow albedo on snow depth in the current snow schemes could potentially cause large overestimation of snow depth in state-of-the-art reanalysis datasets. Implementation of the relationship between fresh snow albedo and snow depth in reanalysis systems can be helpful to understand the overestimation of snow depth on the Tibetan Plateau. This has implications both for analyzing snow cover and water resources under climate change scenarios, and for reconstructing past environments on TP in paleoclimate simulations.
This work was supported by National Key Research and Development Program of China (Grant 2017YFA0603703), National Science Foundation of China (Grant 91537210), Frontier Science Project of Chinese Academy of Sciences (Grant QYZDY-SSW-DQC011-03), and Project supported by the National Science Foundation for Young Scientists of China (Grant 41901075).
Data availability statement: The observed data under CAMP and of Dadongshu-Yakou station is available at http://data.tpdc.ac.cn/en/data/8028b944-daaa-4511-8769-965612652c49/, the website of National Tibetan Plateau Data Center. The observed data of Tanggula station is provided by the Cryosphere Research Station on Qinghai-Xizang Plateau, Chinese Academy of Sciences and is available at http://www.crs.ac.cn. The daily meteorological observations is provided by the National Meteorological Information Center of Chinese Meteorological Administration (CMA).
For diffuse radiation, the snow albedo albdiffuse is calculated by Eq. (A1), with the fresh snow albedo for visible band being 0.95 and for near-infrared band being 0.65:
For the direct radiation, the snow albedo albdirect is calculated by Eq. (A2):
where Fzenith is a factor to account for solar zenith angle impact, calculated as Eq. (A3):
where TFRZ is the freezing temperature set as 273.16 K in Noah-MP, TG is the ground temperature (K), and swet and swet−1 are the snow water equivalents of the current time step and the previous time step. Parameter swec is the critical value of new snow water equivalent to fully cover the old snow, assumed to be 1 mm. Parameter dt is the time step (s) in Noah-MP.
CLASS assumes an identical snow albedo for both direct and diffuse radiation. When there is snowfall and the snow depth surpasses 1 cm, the fresh snow albedo is 0.84; otherwise, the fresh snow albedo is calculated as the weighted average of 0.84 and the albedo of the underground before snowfall, as in Eq. (A11):
where alb is the snow albedo, is the snow albedo before snowfall, swec = 1 mm is the critical value of new snow water equivalent that fully covers the old snow, Qsnow is the snowfall rate (mm s−1), and dt is the modeling time step (s).
The snow albedo is assumed to be no less than 0.55. The snow aging is expressed by an exponential function of the modeling time step dt,
where is the snow albedo before snowfall at the last time step.
c. Introduced snow grain size and snow aging schemes in this study
This study introduced a new snow albedo scheme, which is a combination of a derived fresh snow albedo scheme described in detail in section 4b and the snow grain size scheme from Anderson (1976) and the snow aging scheme from Baker et al. (1990) and Verseghy (1991). Note that BATS and CLASS are not used in the new scheme because they do not calculate snow grain size which is needed in the new scheme.
The snow grain diameter is calculated as follows:
where Dsn is the snow particle diameter in μm, ρsn is snow density in kg m−3, and Tw is the dewpoint temperature in °C.
Snow aging is calculated by
where alb is the snow albedo, albsoil is the soil albedo, and dt is the modeling time step (s). Parameters τf and τa are scaling parameters that are set to 0.24 and 0.008, respectively, according to Sun et al. (1999).