Performance of GLASS and MODIS Satellite Albedo Products in Diagnosing Albedo Variations during Di ﬀ erent Time Scales and Special Weather Conditions in the Tibetan Plateau

: Surface albedo is a crucial parameter in accurately and quantitatively estimating energy and water budget on the Tibetan Plateau (TP) and is also one of the largest radiative uncertainties in land surface modelling attempts. Based on an 8-year ground-based observation of the surface albedo over typical alpine meadows at Maqu and Maduo sites in the eastern TP, the performance of surface albedo products of Global LAnd Surface Satellite (GLASS) and Moderate Resolution Imaging Spectroradiometer (MODIS) in describing albedo variations at daily, 8-day, seasonal timescales, and during di ﬀ erent special weather conditions were analyzed. Compared with the ground-based observation in Maqu, the 8-day albedo products from GLASS and MCD43B3 present maximum negative biases of − 0.030 and − 0.027 at Maqu, respectively. The black-sky albedo (BSA) of GLASS product coincides well with the ground-based observation in Maduo, with root mean square error (RMSE) of 0.092 and correlation coe ﬃ cient (R) of 0.833, whereas that of MCD43B3 had an RMSE of 0.072 and R of 0.752. However, they are underestimated when the albedo is greater than 0.4. At the seasonal timescale, the BSA of GLASS and MCD43B3 underestimated the ground-based observation of Maqu by 0.015 in summer, while their white-sky albedo (WSA) are slightly overestimated and closer to the ground-based observation. In daily timescale, the response of surface albedo to soil moisture is di ﬀ erent in semihumid and semiarid areas in summer. For both sites, the blue-sky-albedo of MCD43A3 has better agreement with the ground-based observation than GLASS and MCD43B3, as it improves the temporal resolution and calculates the albedo by weighting multiple observations within 16 days to be closer to the actual surface. However, even MCD43A3 could not capture the slowdown processes of albedo changes resulted by small snowfall processes or the snow aging due to cloud cover and inversion algorithms. 2-day period and then, comparing the results on a sunny day. RMSE of 0.088 and 0.092 and R of 0.830 and 0.833, respectively. And the WSA and BSA of MCD43B3 have the RMSE of 0.069 and 0.072 and R of 0.743 and 0.752, respectively. The results demonstrate that GLASS albedo is comparable to MCD43B3 albedo in capturing the variation of surface albedo. The temporal and spatial resolutions between them are the same, and the error caused by the scale di ﬀ erence is small. However, the 8-day average albedo of remote sensing products have negative biases and reduced accuracy in characterizing high values of surface albedo.

the products in describing the changes of albedo responses to rapid snow cover changes or different weather conditions such as rainfall and snowfall events are still left to be evaluated.
This article is an advanced study which pays more attention on how about the products reflect the albedo changes in different time scales, particularly on how they performed under different weather situation. A recent study illustrated that the rapid change of albedo caused by the snowfall and the quick snow melting on the TP in early spring and winter can be very important for 3-7 days' short-term climate prediction [31]. Thus, this is the first step of our work to see how the satellite albedo reflect the snowfall-and melting-induced albedo changes, and then, in the second step, a better time varying albedo products will be used to update the lower boundary condition to actively perform its influence on the atmospheric circulation. In this study, we further assessed the MODIS daily albedo product to investigate the albedo changes under special weather conditions and within an 8-day timescale and its relationship with air temperature and 5 cm soil moisture, as well as the main causes of surface albedo product bias were analyzed in detail. The study is outlined as follows. Section 2 describes the remote sensing albedo products, in situ measurements, and methods used in this study. Section 3 provides a performance assessment of the remote sensing albedo products. Section 4 presents a discussion on analysis results and methods. Section 5 concludes with a summary of the findings.

GLASS Product
Global LAnd Surface Satellite (GLASS) is a land surface product developed by Beijing Normal University. GLASS surface albedo product has stable, high-quality data output, and is produced in two steps. The first step is the traditional remote sensing inversion, which is to extract the surface albedo from the remote sensing observation data by inverting the surface-atmospheric radiation transmission model. Then, a variety of primary albedo products are fused and filled by filtering methods to form high-quality fused products. GLASS establishes advanced high-level direct estimation algorithms including Angular Bin and statistics-based temporal filter (STF) with multisource remote sensing observations [37]. The Angular Bin algorithm establishes a linear regression equation between surface directional reflectance and shortwave broadband albedo, specifically the shortwave white-sky albedo and black-sky albedo corresponding to the solar angle at local noon. The STF algorithm regards different intermediate products as samples of the "true" surface albedo with an inversion error and time discrepancy [38]. The fourth edition product named GLASS02A06 is used in this study, which extends to black-sky albedo (BSA), white-sky albedo (WSA) at noon local time, and quality control information in the visible (VIS) and near-infrared (NIR) bands, except for short waves ( Table 1). The data can be downloaded for free at http://glass-product.bnu.edu.cn:8088/.

MODIS Products
Moderate Resolution Imaging Spectroradiometer (MODIS) is a key instrument of the NASA Earth Observation System (EOS) aboard the Terra and Aqua satellites with high spatial and temporal resolution [39,40]. The albedo product number of sinusoidal projections is MCD43, which provides the BSA and WSA albedos for three broadbands and seven spectral bands ( Table 1). The inversed algorithm of MCD43B3 and MCD43A3 products applies the linear kernel-based BRDF model driven by the RossThick-LiSparse-Reciprocal, based on the multiangle observation data accumulated by the change of satellite-surface-sun relative position over a 16-day window [41,42].
MCD43B3 Version 5 (V005) data set is a level-3 gridded product, providing a standard suite of albedo at 1 km resolution. It is widely used in various research fields [19][20][21][22]32,33]. MCD43A3 (V006) albedo with a spatial resolution of 500 m is the latest version of albedo product, weighting observations to estimate BRDF/Albedo on the ninth day of the 16-day period and improving the temporal resolution from 8-day of V005 to 1-day. Its quality information is provided in the corresponding MCD43A2 data file. The data was downloaded from the https://ladsweb.modaps.eosdis.nasa.gov/search/order/.

In Situ Measurements
The Zoige Plateau Wetland Ecosystem Research Station (Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou, Gansu, China), established in 2008, is located in the source region of the Yellow River in the eastern TP ( Figure 1). To select representative measurement sites, the two steps are followed. First, using the camera carried by unmanned drone aircraft (to obtain overall consistency), satellite remote sensing observations (to maintain uniformity within a relative large area), and the Digital Elevation Model (DEM) data (to explore the complexity of the surrounding topography), a relatively representative and uniform region was selected. Second, three sets of instruments including radiation components, automatic weather stations, and eddy covariance were installed to test the data variance and consistency within this region, and then, the best site was selected. According to the footprint model in Maqu, 90% turbulent heat fluxes extent is about 200 m, which also shows the rationality of the site selection [43]. Wang et al. [44] examined the radiation balance and the response of albedo to environmental factors above an alpine meadow and alpine wetland surfaces in the eastern TP, using 2014 data collected at the Zoige Plateau Wetlands Ecosystem Research Stations. The standard deviation of elevation (σ) is a measure of the variability in elevation within the window. The DEM used to calculate the topographic complexity indices was derived from the NASADEM data released in February 2020, with a spatial resolution of 30 m. The σ of 1, 2, and 3 km away from the ground measurement fields of Maqu and Maduo are 2.33, 5.68, 6.70 m and 40.94, 61.39, 71.28 m, respectively. This also proves that Maqu site is very flat, and the Maduo area is hilly (Figure 1d,f). In addition, the vegetation cover types of the ground measurement sites are consistent with that of the surroundings, obtained using MODIS land cover products with a spatial resolution of 0.05 • from

Grassland Observation Site in Maqu
The Maqu site (33°53′ N, 102°08′ E, 3423 m.a.s.l.) is typically characterized by an alpine meadow south of the Maqu County (Figure 1a,c). The land surface is uniformly covered by 0.2 m meadow in the summer and about 0.1 m in the winter. Its soil type is silt clay loam consisting of 29.8% sand, 66.7% silt, and 3.5% clay in the top 40 cm [43]. The annual mean temperature from 1967 to 2010 is 1.7 °C and annual mean precipitation is 600 mm but there is little snow frequency in winter [44]. The radiation sensor at Maqu site is a CNR1 net radiometer (Kipp & Zonen, Delft, The Netherlands), including two CM3 pyranometers and two CG3 pyrgeometers, which records data every 30 min [45]. The parameters of the instrument installed at this site are listed in Table 2. Solar radiation is measured with two CM3 pyranometers, one for measuring incoming solar radiation from the sky and the other facing downward for measuring the reflected solar radiation. From these two pyranometers, albedo, the ratio of reflected and incoming shortwave radiation, can also be determined. The pyrgeometers for measuring far infrared radiation are not involved, as they do not measure shortwave solar radiation. This unitless value ranges between 0 and 1. Checking this can be used as a tool for quality assurance of data. All instruments are calibrated annually to ensure their reliability. The CNR1 should be calibrated using the calibration factor added to each pyranometer or by running a higher standard parallel with it over a 2-day period and then, comparing the results on a sunny day.

Grassland Observation Site in Maqu
The Maqu site (33 • 53 N, 102 • 08 E, 3423 m.a.s.l.) is typically characterized by an alpine meadow south of the Maqu County (Figure 1a,c). The land surface is uniformly covered by 0.2 m meadow in the summer and about 0.1 m in the winter. Its soil type is silt clay loam consisting of 29.8% sand, 66.7% silt, and 3.5% clay in the top 40 cm [43]. The annual mean temperature from 1967 to 2010 is 1.7 • C and annual mean precipitation is 600 mm but there is little snow frequency in winter [44]. The radiation sensor at Maqu site is a CNR1 net radiometer (Kipp & Zonen, Delft, The Netherlands), including two CM3 pyranometers and two CG3 pyrgeometers, which records data every 30 min [45]. The parameters of the instrument installed at this site are listed in Table 2. Solar radiation is measured with two CM3 pyranometers, one for measuring incoming solar radiation from the sky and the other facing downward for measuring the reflected solar radiation. From these two pyranometers, albedo, the ratio of reflected and incoming shortwave radiation, can also be determined. The pyrgeometers for measuring far infrared radiation are not involved, as they do not measure shortwave solar radiation. This unitless value ranges between 0 and 1. Checking this can be used as a tool for quality assurance of data. All instruments are calibrated annually to ensure their reliability. The CNR1 should be calibrated using the calibration factor added to each pyranometer or by running a higher standard parallel with it over a 2-day period and then, comparing the results on a sunny day. ) is located in a grassland about 1.7 km away from the northwest side of Ngoring Lake. There are no fences in the ground measurement field of Maduo, and the vegetation is the same as the surrounding area, with a height of about 0.05-1.10 m ( Table 2). The landscape is homogeneous and wide open with soil dominated by alpine meadow soil, which is rich in nitrogen, phosphorus, and potassium. A few hills with relative height of 50 to 100 m are distributed about 1 km away from the northwest. The instrument for measuring radiation in the site is CNR1 net radiometer (Kipp & Zonen, Delft, The Netherlands) before June 2014, and then CNR4 net radiometer (Kipp & Zonen, Delft, The Netherlands) is used. The CNR1 instrument is the same as that installed at Maqu site. The CNR4 is a four-component net radiometer, consisting of a pyranometer pair, one facing upward and the other facing downward, for measuring shortwave solar radiation [46]. The measured values of these two pyranometers are used to determine the albedo. The output is expressed in W/m 2 . Unlike the CNR1, four probes in the CNR4 have different sensitivity values, which make each measurement from four sensors more accurate. The calibration of albedo from CNR4 is similar to CNR1. Radiation data are recorded every 30 min by the pyranometers installed at a height of 1.5 m at the site.

Calculation of 8-Day Average Albedo of GLASS and MCD43B3
From the data, unreasonable albedos measuring <0 and >1 in the ground-based observations were removed. As the GLASS albedo is calculated at noontime, the ground-based daily mean albedo is the hourly average from 12:00 p.m. to 2:00 p.m. local time in Maqu and Maduo. The albedos of the two pixels closest to the Maqu and Maduo sites in the GLASS and MCD43B3 were averaged, respectively. The substantial aggregation method of GLASS albedo product is used to calculate the arithmetic mean of the valid albedo data selected before and after 8 days (total duration of 16 days) of the specific time. In other words, when comparing the surface albedo averaged over 8 days, the daily albedo on the 9th day is the composite result of the daily surface albedo from day 1 to 16. For an accurate comparison and analysis, the ground-based observations were processed into the 8-day average data and were consistent with the temporal resolution of the GLASS product. Considering the effect of snow on albedo, if there were snow cover during most days of the 16-day period, the albedo at noon is greater than 0.4 and only snow-covered pixels were used for the 8-day albedo retrieval [47].
The surface albedo products of GLASS and MCD43B3 directly supply the BSA and WSA with 1 km spatial resolution, 8-day temporal resolution, and sinusoidal projection, whereas the synthesis periods between GLASS and MCD43B3 are inconsistent. To reduce uncertainty from different data sources, the inversion data during the 16-day period of the sequential accumulation periods from the 9th day of GLASS is compared with that from the 1st day of MCD43B3 and the inversion data from the 17th day of GLASS is compared with that from the 9th day of MCD43B3, and so on.
The surface albedo under actual atmospheric conditions should be weighted based on the fraction of the diffuse radiation to the total radiation and aerosol optical depth. However, it is difficult to obtain the aerosol optical depths at Maqu and Maduo sites. Also, there is no such information in GLASS and MCD43B3 albedo products. Hence, their BSA and WSA albedos were directly evaluated by observations in this study.

The Blue-Sky-Albedo of MCD43A3
First, it needs to calculate the proportion of diffuse illumination, according to the solar zenith angle (SZA) at noon provided by BRDF_Albedo_LocalSolorNoon in the MCD43A2 (V006) product. Then, the linear weighting method is used to compute the blue-sky-albedo (actual surface albedo) using the BSA and WSA albedos of the shortwave provided by the MCD43A3 product (Equations (1) and (2)) [48,49].
where r is the proportion of diffuse illumination, θ is the SZA, α GA is the actual surface albedo, α BSA is the BSA, and α WSA is the WSA. The spatial representativeness of the selected ground observation site is firstly judged. For the ground-based observation data sets with sufficient spatial representation, the accuracy of the remote sensing products to be verified is usually obtained by direct verification. The uniform grasslands at Maqu and Maduo sites perfectly meet this requirement. The albedos of the two pixels closest to the Maqu and Maduo sites in MCD43A3 were averaged.

Interpretation of Statistical Metrics Used
The RMSE, bias, and R had been regularly used as the standard statistical metrics to measure product performance for many years. They are calculated for the data set as in which X and y stand for the surface albedo, x observed and x estimate are the ground-based and satellite-based albedos, and n is the number of data set.

Evaluation of 8-Day Average Albedo
The annual magnitude of ground-based albedo in Maqu area is smaller than that of Maduo, because there are relatively more snowfall events in Maduo ( Figure 2a). Nevertheless, on the inversions beginning on the 49th day of 2014 and the 9th day of 2015, the ground-based albedo rose to 0.67 and 0.65, respectively. Particularly, the albedo at 12:00-2:00 p.m. on 17 February 2014, reached 0.83-0.88. Shi and Liang [50] found that the albedo increase over TP is highly correlated with the increase in snow respectively. GLASS used the STF algorithm that performs time series smoothing and filters fusion of surface albedo, which ignores sharp increase [51]. GLASS's BSA albedo exhibited the highest RMSE of 0.055 and smallest R of 0.327 from January 2009 to June 2017, while the WSA and BSA of MCD43B3 obtained data closer to the ground-based observation with RMSE of 0.055 and R of 0.503 (Table 3 and Figure 3a). Compared to the ground-based observation at Maqu, GLASS and MCD43B3 albedo products have the maximum negative biases of -0.030 and -0.027, especially in spring (Figure 4a). Less snow frequency at Maqu in winter, combined with strong solar radiation, causes the snow to melt quickly. Therefore, the ground-based albedo in 8-day greater than 0.4 is not very large, resulting in the bias of albedo products.
Remote Sens. 2020, 12, x FOR PEER REVIEW 8 of 20 the 9th day of 2015, respectively. GLASS used the STF algorithm that performs time series smoothing and filters fusion of surface albedo, which ignores sharp increase [51]. GLASS's BSA albedo exhibited the highest RMSE of 0.055 and smallest R of 0.327 from January 2009 to June 2017, while the WSA and BSA of MCD43B3 obtained data closer to the ground-based observation with RMSE of 0.055 and R of 0.503 (Table 3 and Figure 3a). Compared to the ground-based observation at Maqu, GLASS and MCD43B3 albedo products have the maximum negative biases of -0.030 and -0.027, especially in spring ( Figure 4a). Less snow frequency at Maqu in winter, combined with strong solar radiation, causes the snow to melt quickly. Therefore, the ground-based albedo in 8-day greater than 0.4 is not very large, resulting in the bias of albedo products.

RMSE Bias RMSE Bias
Unlike Maqu site, the ground-based albedo at the Maduo site has more obvious seasonal cycle, ranging from 0.16 in summer to 0.76 in winter (Figure 2b), and the high values appearing in time series may be closely associated with snow. The appearance of snow in the inversions from the 9th to 25th day of 2012, 1st to 49th day of 2014, and 1st to 9th day of 2015 quickly increased the albedo to around 0.6. Particularly, the 16-day albedo from the 9th day of 2012 and 313rd day of 2013 reached 0.76. GLASS and MCD43B3 albedo significantly increased during these periods, but GLASS is underestimated. At Maduo site, the annual negative bias, RMSE, and R between the ground-based albedo and WSA and BSA of GLASS are 0.019-0.073, 0.026-0.128, and 0.514-0.922, while those of MCD43B3 are 0.015-0.076, 0.040-0.137, and 0.387-0.964, respectively (Figures 3b and 4b). The WSA and BSA of GLASS product coincide well with the ground-based observation in Maduo, with the Remote Sens. 2020, 12, 2456 9 of 20 RMSE of 0.088 and 0.092 and R of 0.830 and 0.833, respectively. And the WSA and BSA of MCD43B3 have the RMSE of 0.069 and 0.072 and R of 0.743 and 0.752, respectively. The results demonstrate that GLASS albedo is comparable to MCD43B3 albedo in capturing the variation of surface albedo. The temporal and spatial resolutions between them are the same, and the error caused by the scale difference is small. However, the 8-day average albedo of remote sensing products have negative biases and reduced accuracy in characterizing high values of surface albedo.      Unlike Maqu site, the ground-based albedo at the Maduo site has more obvious seasonal cycle, ranging from 0.16 in summer to 0.76 in winter (Figure 2b), and the high values appearing in time series may be closely associated with snow. The appearance of snow in the inversions from the 9th to 25th day of 2012, 1st to 49th day of 2014, and 1st to 9th day of 2015 quickly increased the albedo to around 0.6. Particularly, the 16-day albedo from the 9th day of 2012 and 313rd day of 2013 reached 0.76. GLASS and MCD43B3 albedo significantly increased during these periods, but GLASS is underestimated. At Maduo site, the annual negative bias, RMSE, and R between the ground-based albedo and WSA and BSA of GLASS are 0.019-0.073, 0.026-0.128, and 0.514-0.922, while those of MCD43B3 are 0.015-0.076, 0.040-0.137, and 0.387-0.964, respectively (Figures 3b and 4b). The WSA and BSA of GLASS product coincide well with the ground-based observation in Maduo, with the RMSE of 0.088 and 0.092 and R of 0.830 and 0.833, respectively. And the WSA and BSA of MCD43B3 have the RMSE of 0.069 and 0.072 and R of 0.743 and 0.752, respectively. The results demonstrate that GLASS albedo is comparable to MCD43B3 albedo in capturing the variation of surface albedo. The temporal and spatial resolutions between them are the same, and the error caused by the scale difference is small. However, the 8-day average albedo of remote sensing products have negative biases and reduced accuracy in characterizing high values of surface albedo.

Seasonal Distribution of 8-Day Albedo
Surface albedo is a dynamic parameter influenced by the vegetation phenology, soil moisture, soil color, surface roughness, SZA, cloud, snow cover, etc. [52][53][54][55]. Figure 5 shows the relationship in the 8-day albedo between satellite-based albedo and ground-based observation in different seasons.

Seasonal Distribution of 8-Day Albedo
Surface albedo is a dynamic parameter influenced by the vegetation phenology, soil moisture, soil color, surface roughness, SZA, cloud, snow cover, etc. [52][53][54][55]. Figure 5 shows the relationship in the 8-day albedo between satellite-based albedo and ground-based observation in different seasons. At Maqu site, the ground-based seasonal variation of surface albedo is in the order of winter >spring >autumn >summer with means of 0.25, 0.22, 0.18, and 0.19, respectively. The ground-based albedos vary with the largest range in winter (0.20-0.67) and the smallest in summer (0.17-0.22), while the satellite-based albedo are 0.18-0.67 and 0.15-0.23, respectively. Massive albedo data have good agreement with the ground-based observation of Maqu. In summer, the BSA of GLASS and MCD43B3 underestimated a value of about 0.015, while their WSA are slightly overestimated and closer to the ground-based observation. GLASS and MCD43B3 products provide both WSA-diffuse reflectance and BSA-direct reflectance, i.e., surface albedo on sunny and cloudy days. In spring and summer, apparent diurnal variations of clouds and precipitation cells were observed on the TP, and the amplitude of diurnal cloud activity fluctuates on intraseasonal time scales [56]. Generally, many convective-type clouds in Maqu develop quickly after 11:00 a.m. local time and exist until mid-night in summer. The transit time of MODIS in the area is approximately 1:00 p.m. local time, and its monitoring is significantly affected by clouds. And if the atmospheric conditions are hazy, the atmospheric correction algorithm may introduce error in the high-level products of GLASS and MODIS [17]. A few albedos of GLASS and MCD43B3 in winter have overestimation, the maximum bias of GLASS albedo occurred in January and February of about 0.13, while that of MCD43B3 is 0.14 in January. The atmospheric correction is performed by individual aerosol estimates for each data point in the algorithm, while a static model is currently employed in the land-based analysis operational system [57]. Furthermore, the SZA on the TP in winter is large. The accuracy of the MCD43B3 albedo decreases significantly when the SZA exceeds 70 • [20,58].
Remote Sens. 2020, 12, x FOR PEER REVIEW 10 of 20 operational system [57]. Furthermore, the SZA on the TP in winter is large. The accuracy of the MCD43B3 albedo decreases significantly when the SZA exceeds 70° [20,58].  (Figure 6a,c). However, the WSA and BSA of GLASS are less than the mean ground-based observation by 0.031 and 0.042, respectively, and those of MCD43B3 are 0.046 and 0.053, respectively. The fluctuation of surface albedo in winter is intuitively clear at Maduo site, and the satellite-based albedos have higher dispersion. The 8-day average albedo of ground-based observation in winter ranges from 0.21 to 0.76, may be mainly caused by (1) its higher elevations than Maqu, (2) more frequent precipitation, (3) long duration of snow cover, and (4) significant vegetation changes. The maximum biases of GLASS and MCD43B3 albedos from the ground-based observation in winter are 0.155 and 0.171. The "dark-object" method in MCD43B3 product is used in the current operational atmospheric correction algorithms for estimating aerosol optical depth, which fails over nonvegetated regions that often have much larger albedo values, such as snow, ice, and bare soils [59]. It also displays that the STF algorithm in GLASS has some shortcomings, such as the fact that filtering fusion is likely to filter out some sharply changed points, which correspond to the actual situation [60].  (Figure 6a,c). However, the WSA and BSA of GLASS are less than the mean ground-based observation by 0.031 and 0.042, respectively, and those of MCD43B3 are 0.046 and 0.053, respectively. The fluctuation of surface albedo in winter is intuitively clear at Maduo site, and the satellite-based albedos have higher dispersion. The 8-day average albedo of ground-based observation in winter ranges from 0.21 to 0.76, may be mainly caused by (1) its higher elevations than Maqu, (2) more frequent precipitation, (3) long duration of snow cover, and (4) significant vegetation changes. The maximum biases of GLASS and MCD43B3 albedos from the ground-based observation in winter are 0.155 and 0.171. The "dark-object" method in MCD43B3 product is used in the current operational atmospheric correction algorithms for estimating aerosol optical depth, which fails over nonvegetated regions that often have much larger albedo values, such as snow, ice, and bare soils [59]. It also displays that the STF algorithm in GLASS has some shortcomings, such as the fact that filtering fusion is likely to filter out some sharply changed points, which correspond to the actual situation [60].

Effect of Snowfall on Daily Albedo
Snow has a high albedo and can modify the ground energy balance [35]. Li et al. [31] have pointed out that the subseasonal variability of TP snow cover is closely related to the subsequent East Asian atmospheric circulation at medium-range time scales (approximately 3-8 days later) during wintertime. However, various snow processes that might be important may be missing or misrepresented. The daily ground-based albedo and satellite-based albedo are used to study the effect of snow. As shown in Figure 7a, continuous snowfall occurred at Maqu site from 1 to 17 January 2015. The maximum ground-based albedo remained 0.21 from 1 to 4 January, while it sharply increased to 0.75, 0.72, and 0.69, respectively, when the snowfall increased to 2.7 mm from 5 to 7 January. Meanwhile, the temperatures at noon were less than 0 °C and the daily average temperatures were below -10 °C. When the snowfall is small, the ground-based albedo decreases accordingly, about 0.24-0.29 on 15-17 January. The albedo of MCD43A3 was 0.19 from January 1 to 5, suddenly increased to 0.66 on the 6th, and continued until the 10th, and then, it decreased to 0.19. This indicates the superiority of the MCD43A3 albedo product. For the 8-day average albedo, GLASS and MCD43B3 product were less than the ground-based albedo at Maqu site, and the effects of snowfall on albedo were not presented (Figures 4a and 7a). MCD43A3 albedo is more accurate than GLASS and MCD43B3 in the magnitude of albedo during the snowfall periods, embodying its enhanced capability to separate most snow and clouds [61]. However, it missed capturing some short or small snowfall processes such as the slowdown processes of albedo changes from 11 to 16 January. GLASS and MCD43B3 products applied the bidirectional reflection model to calculate reflectance assuming a constant surface in a certain period, leading to the missing or misrepresentation of rapid changes in albedo over TP [60]. The MCD43A3 product is accumulated based on the ground data of 8-day before and after the center day, and then, it assigns different weights according to the distance of the observation date, thereby producing the albedo of the marked date through the weighting of multiple days.

Effect of Snowfall on Daily Albedo
Snow has a high albedo and can modify the ground energy balance [35]. Li et al. [31] have pointed out that the subseasonal variability of TP snow cover is closely related to the subsequent East Asian atmospheric circulation at medium-range time scales (approximately 3-8 days later) during wintertime. However, various snow processes that might be important may be missing or misrepresented. The daily ground-based albedo and satellite-based albedo are used to study the effect of snow. As shown in Figure 7a, continuous snowfall occurred at Maqu site from 1 to 17 January 2015. The maximum ground-based albedo remained 0.21 from 1 to 4 January, while it sharply increased to 0.75, 0.72, and 0.69, respectively, when the snowfall increased to 2.7 mm from 5 to 7 January. Meanwhile, the temperatures at noon were less than 0 • C and the daily average temperatures were below -10 • C. When the snowfall is small, the ground-based albedo decreases accordingly, about 0.24-0.29 on 15-17 January. The albedo of MCD43A3 was 0.19 from January 1 to 5, suddenly increased to 0.66 on the 6th, and continued until the 10th, and then, it decreased to 0.19. This indicates the superiority of the MCD43A3 albedo product. For the 8-day average albedo, GLASS and MCD43B3 product were less than the ground-based albedo at Maqu site, and the effects of snowfall on albedo were not presented (Figures 4a and 7a). MCD43A3 albedo is more accurate than GLASS and MCD43B3 in the magnitude of albedo during the snowfall periods, embodying its enhanced capability to separate most snow and clouds [61]. However, it missed capturing some short or small snowfall processes such as the slowdown processes of albedo changes from 11 to 16 January. GLASS and MCD43B3 products applied the bidirectional reflection model to calculate reflectance assuming a constant surface in a certain period, leading to the missing or misrepresentation of rapid changes in albedo over TP [60]. The MCD43A3 product is accumulated based on the ground data of 8-day before and after the center day, and then, it assigns different weights according to the distance of the observation date, thereby producing the albedo of the marked date through the weighting of multiple days.  Figure 7b shows the snowfall in Maduo from 2 to 18 December 2016. On 2 December, the groundbased albedo at noon reached 0.8 with 2.2 mm snowfall. The daily average albedo rapidly rose from 0.21 to 0.75, while the temperature dropped from −8.8 to −11.4 °C. When the snow coverage reached 100%, the ground-based albedo varied from 0.55 to 0.68. As the temperature increased, the snow began to melt and the ground-based albedo was only around 0.23 from December 11 to 18. The bluesky-albedo of MCD43A3 had a similar pattern, which was 0.79 on 2 December, and then maintained at 0.21. The WSA and BSA of GLASS reduced from 0.51 to 0.21, highlighting a capture of the slowdown processes of snow albedo. However, the WSA and BSA of MCD43B3 almost presented an opposite trend during this period, increasing from 0.19 to 0.27, because the algorithm included in it is mainly directed at surfaces with plants and the accuracy of highly reflective snow is not sufficient [60]. When there was no snowfall after 11 December, the daily average albedo based on the ground was 0.22-0.32. The biases of albedo products from the ground measurements were small, and MCD43A3 albedo performed relatively well with a bias of 0.02. During the cold season in the absence of snow, the surface is dominated by sparsely dried grass and naked soil, inducing little changes in the albedo. Snow cover makes the surface albedo reach the maximum value of the whole year with albedos of approximately 0.9 for fresh snow and 0.7 for old snow [52]. Snow albedo on the TP depends on the both snow depth and snow density. And the dry-cold snow cover on the TP is characterized by small density and irregular distribution [53]. Although the MODIS sensor has an absolute advantage in snow monitoring, it does not detect snow with depths of less than 0.5 cm [62]. The patchy snow and reduced snow accumulation by sublimation on TP increase the bias of GLASS and MODIS product [35]. In addition, the snow cover in GLASS and MCD43B3 albedo is determined by the presence of snow for most of the time in 16 days, leading to the removal of snow data in the short term. The MCD43A3 product uses the snow status weighted to the ninth day instead of the majority snow observations from the 16-day period. These results are reasonable as the MCD43A3 algorithm  Figure 7b shows the snowfall in Maduo from 2 to 18 December 2016. On 2 December, the ground-based albedo at noon reached 0.8 with 2.2 mm snowfall. The daily average albedo rapidly rose from 0.21 to 0.75, while the temperature dropped from −8.8 to −11.4 • C. When the snow coverage reached 100%, the ground-based albedo varied from 0.55 to 0.68. As the temperature increased, the snow began to melt and the ground-based albedo was only around 0.23 from 11 to 18 December. The blue-sky-albedo of MCD43A3 had a similar pattern, which was 0.79 on 2 December, and then maintained at 0.21. The WSA and BSA of GLASS reduced from 0.51 to 0.21, highlighting a capture of the slowdown processes of snow albedo. However, the WSA and BSA of MCD43B3 almost presented an opposite trend during this period, increasing from 0.19 to 0.27, because the algorithm included in it is mainly directed at surfaces with plants and the accuracy of highly reflective snow is not sufficient [60]. When there was no snowfall after 11 December, the daily average albedo based on the ground was 0.22-0.32. The biases of albedo products from the ground measurements were small, and MCD43A3 albedo performed relatively well with a bias of 0.02. During the cold season in the absence of snow, the surface is dominated by sparsely dried grass and naked soil, inducing little changes in the albedo. Snow cover makes the surface albedo reach the maximum value of the whole year with albedos of approximately 0.9 for fresh snow and 0.7 for old snow [52]. Snow albedo on the TP depends on the both snow depth and snow density. And the dry-cold snow cover on the TP is characterized by small density and irregular distribution [53]. Although the MODIS sensor has an absolute advantage in snow monitoring, it does not detect snow with depths of less than 0.5 cm [62]. The patchy snow and reduced snow accumulation by sublimation on TP increase the bias of GLASS and MODIS product [35]. In addition, the snow cover in GLASS and MCD43B3 albedo is determined by the presence of snow for most of the time in 16 days, leading to the removal of snow data in the short term. The MCD43A3 product uses the snow status weighted to the ninth day instead of the majority snow observations from the 16-day period. These results are reasonable as the MCD43A3 algorithm emphasizes the contribution of the single day of interest more strongly, which enables it to capture subtle changes in surface albedo [15].

Effect of Rainfall on Daily Albedo
There were three dominant rainfalls in the Maqu region from 28 July to 13 August 2014, during 28-31 July, 2-8 and 10-13 August. Typical weather conditions with rainfall of more than 20 mm occurred on the 5 and 8 August 2014. As can be seen in Figure 8a, the ground-based albedo ranged from 0.16 to 0.21 during this period. On August 8, the ground-based albedo decreased from 0.20 to 0.17 with 23.0 mm rainfall and high soil moisture content of 0.38 m 3 m −3 . The albedo then rose to 0.20 in the absence of rainfall, and the soil moisture decreased by 5 cm. The relationship among rainfall, soil moisture content, soil color, and albedo were found to be more complicated on the grassland surface. The blue-sky-albedo of MCD43A3 increased twice on 30 July and 5 August, but did not capture the significant decrease in albedo on July 31 and August 8. The ground-based 8-day average albedo was 0.18-0.19, at which time, the BSA of GLASS and MCD43B3 were underestimated and their WSA were overestimated. The cloud activity on the TP is the most active and complicated in the summer afternoon, and its influence on albedo is dominant and non-negligible [56]. Even with some degree of cloud discrimination provided by MODIS, it does not always reliably separate optically thin cirrus clouds [62].
Remote Sens. 2020, 12, x FOR PEER REVIEW 13 of 20 emphasizes the contribution of the single day of interest more strongly, which enables it to capture subtle changes in surface albedo [15].

Effect of Rainfall on Daily Albedo
There were three dominant rainfalls in the Maqu region from 28 July to 13 August 2014, during 28-31 July, 2-8 and 10-13 August. Typical weather conditions with rainfall of more than 20 mm occurred on the 5 and 8 August 2014. As can be seen in Figure 8a, the ground-based albedo ranged from 0.16 to 0.21 during this period. On August 8, the ground-based albedo decreased from 0.20 to 0.17 with 23.0 mm rainfall and high soil moisture content of 0.38 m 3 m −3 . The albedo then rose to 0.20 in the absence of rainfall, and the soil moisture decreased by 5 cm. The relationship among rainfall, soil moisture content, soil color, and albedo were found to be more complicated on the grassland surface. The blue-sky-albedo of MCD43A3 increased twice on 30 July and 5 August, but did not capture the significant decrease in albedo on July 31 and August 8. The ground-based 8-day average albedo was 0.18-0.19, at which time, the BSA of GLASS and MCD43B3 were underestimated and their WSA were overestimated. The cloud activity on the TP is the most active and complicated in the summer afternoon, and its influence on albedo is dominant and non-negligible [56]. Even with some degree of cloud discrimination provided by MODIS, it does not always reliably separate optically thin cirrus clouds [62]. Five rainfall events from 25 June to 13 July 2012, were analyzed at the Maduo site, as shown in Figure 8b. With a rainfall of 14.2, 7.8, and 9.2 mm, respectively, for 26-28 days, the ground-based albedo decreased from 0.20 to 0.16 and the 5 cm soil moisture increased from 0.19 to 0.26 m 3 m −3 . The reduction in albedo is due to the direct effect of rainfall on soil moisture, which increased the soil moisture content and soil color. These results agreed well with Liu et al. [63] that showed a decrease in albedo when the moisture increased for low soil moisture levels. When rainfall reappeared on July 5 and 10, the 5 cm soil moisture increased again and the ground-based albedo decreased Five rainfall events from 25 June to 13 July 2012, were analyzed at the Maduo site, as shown in Figure 8b. With a rainfall of 14.2, 7.8, and 9.2 mm, respectively, for 26-28 days, the ground-based albedo decreased from 0.20 to 0.16 and the 5 cm soil moisture increased from 0.19 to 0.26 m 3 m −3 . The reduction in albedo is due to the direct effect of rainfall on soil moisture, which increased the soil moisture content and soil color. These results agreed well with Liu et al. [63] that showed a decrease in albedo when the moisture increased for low soil moisture levels. When rainfall reappeared on July 5 and 10, the 5 cm soil moisture increased again and the ground-based albedo decreased simultaneously. GLASS, MCD43B3, and MCD43A3 albedos were significantly smaller than the ground-based observation, especially the reduction of MCD43B3 from 0.16 to 0.12.
The relationship between surface albedo and soil moisture was explored using the daily observation data of ground-based albedo and soil moisture content at the 5 cm depth from June to August at Maqu site in 2014 and Maduo site in 2012. Figure 9 shows the difference in surface albedo response to soil moisture content in semihumid Maqu and semiarid Maduo regions. The goodness of fit of linear regression of surface albedo and soil moisture in Maqu and Maduo sites were 0.0265 and 0.2772, respectively. The grassland grows well from June to August in the Maqu area, with uniform coverage and high height. The color of the underlying grassland does not change when rainfall occurs, and the effect of soil moisture on the surface albedo is very limited. In semiarid areas, changes in soil moisture affect albedo to a certain extent, due to sparse grass and low height. Besides, the cloud cover will partially mask the dependence of albedo on moisture content [63].
Remote Sens. 2020, 12, x FOR PEER REVIEW  14 of 20 simultaneously. GLASS, MCD43B3, and MCD43A3 albedos were significantly smaller than the ground-based observation, especially the reduction of MCD43B3 from 0.16 to 0.12. The relationship between surface albedo and soil moisture was explored using the daily observation data of ground-based albedo and soil moisture content at the 5 cm depth from June to August at Maqu site in 2014 and Maduo site in 2012. Figure 9 shows the difference in surface albedo response to soil moisture content in semihumid Maqu and semiarid Maduo regions. The goodness of fit of linear regression of surface albedo and soil moisture in Maqu and Maduo sites were 0.0265 and 0.2772, respectively. The grassland grows well from June to August in the Maqu area, with uniform coverage and high height. The color of the underlying grassland does not change when rainfall occurs, and the effect of soil moisture on the surface albedo is very limited. In semiarid areas, changes in soil moisture affect albedo to a certain extent, due to sparse grass and low height. Besides, the cloud cover will partially mask the dependence of albedo on moisture content [63]. The validation results show that MCD43A3 albedo has better agreement with the ground-based albedo than the GLASS and MCD43B3 for both sites and can also capture the rapid changes of surface albedo in the case of snowfall and rainfall within 8 days. The main reasons are that MCD43A3 product improves the temporal resolution and is more suitable for reflecting the albedo changes in a short period and it calculates the albedo by weighting multiple observations within 16 days to be closer to the actual surface. Figure 10 shows scatter plots of the 8-day average albedo of ground-based albedo and surface albedo products of GLASS and MCD43B3. At Maqu site, a large amount of the WSA and BSA of GLASS and MCD43B3 were concentrated in the low-value area of the 1:1 line, with small negative bias (Figure 10a,b). Both GLASS and MCD43B3 albedos had individual cases, which were significantly smaller or larger than the ground-based observation of Maqu. The distribution of surface albedo in Maduo site was more discrete. The coincidence of the WSA and BSA of GLASS and MCD43B3 and ground-based observation was relatively high in the low-value range, but they were underestimated when the albedo was greater than 0.4 (Figure 10c,d). The 8-day albedo of this period happened to be averaged based on the number of snow-covered days. When the albedo at Maduo site varied from 0.5 to 0.8, the WSA and BSA of MCD43B3 were 0.4-0.8 and those of GLASS were only 0.2-0.6. The blue-sky-albedo distribution of MCD43A3 was further verified (Figure 11), and the data were concentrated on the 1:1 line when the ground-based albedo were less than 0.3, but the albedos had varying degrees of dispersion when ground-based observations were greater than 0.4 at Maqu and Maduo sites. The validation results show that MCD43A3 albedo has better agreement with the ground-based albedo than the GLASS and MCD43B3 for both sites and can also capture the rapid changes of surface albedo in the case of snowfall and rainfall within 8 days. The main reasons are that MCD43A3 product improves the temporal resolution and is more suitable for reflecting the albedo changes in a short period and it calculates the albedo by weighting multiple observations within 16 days to be closer to the actual surface. Figure 10 shows scatter plots of the 8-day average albedo of ground-based albedo and surface albedo products of GLASS and MCD43B3. At Maqu site, a large amount of the WSA and BSA of GLASS and MCD43B3 were concentrated in the low-value area of the 1:1 line, with small negative bias (Figure 10a,b). Both GLASS and MCD43B3 albedos had individual cases, which were significantly smaller or larger than the ground-based observation of Maqu. The distribution of surface albedo in Maduo site was more discrete. The coincidence of the WSA and BSA of GLASS and MCD43B3 and ground-based observation was relatively high in the low-value range, but they were underestimated when the albedo was greater than 0.4 (Figure 10c,d). The 8-day albedo of this period happened to be averaged based on the number of snow-covered days. When the albedo at Maduo site varied from 0.5 to 0.8, the WSA and BSA of MCD43B3 were 0.4-0.8 and those of GLASS were only 0.2-0.6. The blue-sky-albedo distribution of MCD43A3 was further verified (Figure 11), and the data were  From the above results, it can be concluded that GLASS and MCD43B3 albedo products have certain accuracy at the sites of surface cover type represented by alpine meadows on TP, which can basically reflect the annual and seasonal changes of the surface albedo. GLASS and MODIS products use different inversion algorithms, namely, STF and AMBRALS. The final GLASS product is obtained by fusing its initial products and MCD43 products through the STF algorithm. GLASS albedo change is stable when reflecting a long albedo time sequence. The R with the ground-based albedo of Maduo site reached 0.830 from 2011 to 2016. The accuracy of GLASS product relies on Angular Bin and STF algorithms. The Angular Bin algorithm uses the input data from a single angle, which is suboptimal in capturing the anisotropic reflectance behavior [15,60]. STF algorithm effectively reduces the missing data caused by clouds, but it also causes data distortion in the fusion process. Although GLASS captured the slowdown processes of snow albedo at Maduo site from December 2-10 December 2016, many processes of rapid changes in albedo are still missing.

Discussion
Both MCD43A3 and MCD43B3 albedos use the AMBRALS algorithm to achieve a linear description of the bidirectional reflectance on the basis of BRDF model. In the kernel-driven semiempirical approach, the BRDF is modeled as a weighted sum of a volume scattering function and a surface scattering function (called kernels) and a constant [64]. This combination has proved to  From the above results, it can be concluded that GLASS and MCD43B3 albedo products have certain accuracy at the sites of surface cover type represented by alpine meadows on TP, which can basically reflect the annual and seasonal changes of the surface albedo. GLASS and MODIS products use different inversion algorithms, namely, STF and AMBRALS. The final GLASS product is obtained by fusing its initial products and MCD43 products through the STF algorithm. GLASS albedo change is stable when reflecting a long albedo time sequence. The R with the ground-based albedo of Maduo site reached 0.830 from 2011 to 2016. The accuracy of GLASS product relies on Angular Bin and STF algorithms. The Angular Bin algorithm uses the input data from a single angle, which is suboptimal in capturing the anisotropic reflectance behavior [15,60]. STF algorithm effectively reduces the missing data caused by clouds, but it also causes data distortion in the fusion process. Although GLASS captured the slowdown processes of snow albedo at Maduo site from December 2-10 December 2016, many processes of rapid changes in albedo are still missing.
Both MCD43A3 and MCD43B3 albedos use the AMBRALS algorithm to achieve a linear description of the bidirectional reflectance on the basis of BRDF model. In the kernel-driven semiempirical approach, the BRDF is modeled as a weighted sum of a volume scattering function and a surface scattering function (called kernels) and a constant [64]. This combination has proved to be very suitable for describing the bidirectional reflectance anisotropy of the land surface for different From the above results, it can be concluded that GLASS and MCD43B3 albedo products have certain accuracy at the sites of surface cover type represented by alpine meadows on TP, which can basically reflect the annual and seasonal changes of the surface albedo. GLASS and MODIS products use different inversion algorithms, namely, STF and AMBRALS. The final GLASS product is obtained by fusing its initial products and MCD43 products through the STF algorithm. GLASS albedo change is stable when reflecting a long albedo time sequence. The R with the ground-based albedo of Maduo site reached 0.830 from 2011 to 2016. The accuracy of GLASS product relies on Angular Bin and STF algorithms. The Angular Bin algorithm uses the input data from a single angle, which is suboptimal in capturing the anisotropic reflectance behavior [15,60]. STF algorithm effectively reduces the missing data caused by clouds, but it also causes data distortion in the fusion process. Although GLASS captured the slowdown processes of snow albedo at Maduo site from 2 December to 10 December 2016, many processes of rapid changes in albedo are still missing.
Both MCD43A3 and MCD43B3 albedos use the AMBRALS algorithm to achieve a linear description of the bidirectional reflectance on the basis of BRDF model. In the kernel-driven semiempirical approach, the BRDF is modeled as a weighted sum of a volume scattering function and a surface scattering function (called kernels) and a constant [64]. This combination has proved to be very suitable for describing the bidirectional reflectance anisotropy of the land surface for different plants [65]. However, the AMBRALS algorithm usually requires long-term sequential ground measurements to calculate albedo, resulting in a low temporal resolution, which cannot reflect rapid changes of albedo during rainfall, snowfall, or rapid plant growth. MCD43A3 product have shown some advantages over MCD43B3 in characterizing the albedo changes during the snowfall and rainfall periods, because it uses L2G-lite surface reflectance as input and uses a database with archetypal BRDF parameters to supplement the observational data and performs a lower quality magnitude inversion when high-quality reflectance is insufficient. Also, the database will be continually updated with the latest full inversion retrieval for each pixel. Although the algorithm of MCD43A3 product has been considerably improved, it still missed capturing the slowdown processes of albedo changes resulted by small snowfall processes or the snow aging, and its accuracy is greatly affected by cloud cover.
Remote sensing data matches nonremote sensing data and input into quantitative remote sensing models to obtain the products through inversion. Extremely complex topography, diverse surface types, and severe weather changes have increased the difficulty of obtaining data on the TP. More snowfalls in spring, strong solar radiation, wind-blowing snow, thin snow cover, and fast snowmelt in the study area are also different from other areas. Davidson and Wang [6] have pointed out that snow cover increases not only the albedo of all cover types but also their sensitivities to changes in SZA. The density of ground observation sites on the TP is rarely sufficient to accurately characterize albedo at a regional scale. MCD43A3 product has better quality at high latitudes by using all available observations during the acquisition period, while the MCD43B3 is limited to only four valid observations per day, so that it cannot capture many subtle dynamic changes in albedo.
Accurately characterizing the diurnal variation of surface albedo is very important for monitoring radiative transfers between soil, vegetation, and atmosphere [66]. Most current climate models rely on field measurements to calculate or deduce surface albedo [67]. The absolute accuracy required ranges from 0.02 to 0.05 [68]. But GLASS and MODIS albedo products are still limited by temporal resolution. For example, in the Maqu region, the morning snowfall has melted in the afternoon due to strong solar radiation at noon. Even if the MODIS sensor transits twice a day, it cannot detect the rapid diurnal changes in snow accumulation and melting. At this time, e.g., Chinese geostationary Fengyun meteorological satellite data could be used as a supplement. In addition, GLASS and MODIS albedos have obvious biases in winter due to the effects of snow and cloud cover. For example, none of them reflected the gradual change of the albedo of Maqu in Figure 7a from 12 to 15 January 2015.
Scale gap between ground-based observation and pixel scale is another factor affecting the accuracy of remote sensing albedo product in direct verification. Certainly, more ground-based measurements or upscaling observational technique should be further performed for advanced evaluation. Improving the temporal resolution of remote sensing albedo product and accurately reflecting the rapid change in albedo are still difficult for future inversion. Despite the limited geographic region of the study area, it still provides a basis for the accuracy verification of surface albedo products and also lays a foundation for the direct application of remote sensing albedo products in the land surface model.

Conclusions
Based on the 8-year ground-based observation of the surface albedo over typical alpine meadows at Maqu and Maduo sites in the eastern TP, the performance of surface albedo products of GLASS and MODIS were assessed to investigate surface albedo changes at daily, 8-day, seasonal timescales, and during different special weather conditions, as well as the main causes of surface albedo product bias were analyzed in detail. The main conclusions are as follows: (1) In an 8-day timescale, the ground-based albedos are concentrated in 0. 16 in semihumid and semiarid areas in summer. For both sites, the blue-sky-albedo of MCD43A3 has better agreement with the ground-based albedo than the GLASS and MCD43B3, as it is more suitable for reflecting the albedo changes in a short period by improving the temporal resolution, and calculates the albedo by weighting multiple observations within 16 days to be closer to the actual surface. However, it still missed capturing the slowdown processes of albedo changes resulted by small snowfall processes or the snow aging due to the cloud cover and inversion algorithms. GLASS and MCD43B3 albedos are frequently small and have the missing or misrepresentation of rapid changes in albedo caused by snowfall and rainfall on the TP.