Spatiotemporal Dynamics of Land Surface Albedo and Its Inﬂuencing Factors in the Qilian Mountains, Northeastern Tibetan Plateau

: Land surface albedo directly determines the distribution of radiant energy between the surface and the atmosphere, and it is a key parameter affecting the energy balance on the land surface. However, the spatiotemporal dynamics of land surface albedo and associated inﬂuencing factors in the Qilian Mountains (QM) have been rarely reported. By using the long-time series data products of MODIS shortwave albedo, normalized difference vegetation index (NDVI), and snow cover with a spatial resolution of 0.05 ◦ from 2001 to 2020, this paper analyzes the temporal and spatial variations of land surface albedo in the QM over the past 20 years and its inﬂuencing factors. The analysis results show that the multi-year average surface albedo in the QM has obvious differences in spatial distribution: it increases with the altitude, and it is high in the west (at the west of 98 ◦ E) and low in the east. Meanwhile, the surface albedo has different distribution characteristics in different seasons: the spatial distribution of surface albedo is similar in spring and autumn; the areas with a high surface albedo in summer are signiﬁcantly fewer than those in other seasons. Besides, in the past 20 years, the annual average surface albedo has shown a weak growth trend in the QM, with a change rate of 5 × 10 − 3 /10a, and the minimum and maximum values were reached in 2001 and 2019, respectively. In addition, the annual variation of the surface albedo in the QM showed a “U” shape, with the largest variation in January and the smallest variation in August. The annual variation of surface albedo is signiﬁcantly positively correlated with snow cover (r = 0.96) and signiﬁcantly negatively correlated with NDVI (r = − 0.91). Moreover, the interannual variation of the surface albedo in the QM is closely related to land surface cover and is greatly affected by snow cover. Spatially, the annual variation of surface albedo in most areas of the QM is dominated by the change of snow cover, and the increase of surface albedo in the middle area is consistent with the increase of snow cover, while the decrease of albedo in the edge area is related to the improvement of vegetation cover. The results of this study provide a scientiﬁc basis for studying the climate and environmental changes caused by changes in the surface of the QM and making ecological environment restoration strategies.


Introduction
Since the Industrial Revolution, global warming has become increasingly serious, and multiple indicators of the Earth's climate system such as the atmosphere, ocean, land, and cryosphere have shown the impact of human activities (IPCC). The surface albedo is the ratio of the reflected solar radiation on the surface to the incident solar radiation, which glaciers in mountainous areas below 4000 m have completely disappeared. Compared with the 1960s, the area of glaciers in the QM has decreased by 20.5% ± 6.04%, and the loss rate between 2007 and 2015 is as high as 5.82%/10a [29]. In addition, with the intensification of human activities, the areas of industrial and mining, housing, transportation, and other land uses have increased [27]. Furthermore, the extensive development of mineral resources, overgrazing of local grasslands, irregular operation of tourism facilities, and overloaded groundwater use has seriously damaged the ecological environment of the QM and posed serious challenges to the sustainable development of the society and economy [31].
Affected by various driving factors, surface albedo exhibits obvious temporal variability and spatial heterogeneity. To obtain credible and accurate albedo data, satellite remote sensing is indispensable. Compared with measured and simulated data, satellite remote sensing has a wider coverage, and fewer natural limitations and uncertainties. Thus, it has become an efficient method to acquire continuous albedo [32]. There is good consistency between the satellite-retrieved annual mean and field measurements, when the surface albedo varies seasonally, forest-covered areas are better matched than nonforested areas [33]. In satellite datasets, all latitude assessments show good agreement in summer, while albedo assessments have significant uncertainty in winter, especially at high latitudes [34]. The MODIS V006 albedo product has an improved temporal resolution and accuracy, and the daily MODIS V006 product can reproduce the dynamics of albedo well [32,35,36]. In recent years, some scholars have studied the albedo of glaciers in the QM. It was found that the dust on the surface of the glacier significantly reduced the albedo [18], and there was a significant negative correlation between the diurnal variation of glacial albedo and air temperature [37]. Another study found a significant positive correlation between glacier area change and annual mean albedo [38]. However, the spatiotemporal variation characteristics of surface albedo and its influencing factors in the QM are less studied. Therefore, this paper uses the MODIS MCD43C3 V006 product provided by NASA (National Aeronautics and Space Administration, NASA, Washington, DC, USA), combined with NDVI and snow cover data, to analyze the spatiotemporal distribution of surface albedo and its influencing factors in the QM. The results of the study will help to understand the feedback mechanism of surface albedo and reveal the relation between climate change and human activities, thus providing an important basis for ecological barrier construction and sustainable development in the QM.

Overview of the Study Area
The QM are located on the border between the northeastern Qinghai Province and the western Gansu Province in China, and it is also located at the intersection of three major plateaus, i.e., Qinghai-Tibet, Mengxin, and Loess. It consists of a number of parallel mountains and wide valleys from northwest to southeast (Figure 1a). Because of the high altitude of the overall mountain range, the peaks over 4000 m are covered with snow all year round. Meanwhile, there are 3306 modern glaciers [39], which are the source of the Shiyang River, Heihe River, Shule River, and other rivers [40]. As the altitude decreases from northwest to southeast, the land cover types are bare land, grassland, and cultivated land in sequence [41]. Bare land is mainly in the west of 98 • E, and grassland is mainly in the east of 98 • E (Figure 1b). The annual average temperature in the QM is below 4 • C, higher in barren/rocky areas, lower in snow/glacier areas [42], and the temperature gradually decreases with the increase of elevation. Different from air temperature, precipitation is not only affected by altitude, aspect, and slope but also by latitude, longitude, and atmospheric circulation; it shows a decreasing trend from east to west [28]. Remote Sens. 2022, 14, x FOR PEER REVIEW 4 of 17

Dataset and Preprocessing
This study used three sets of remote sensing products and reanalysis datasets (Table  1), including MODIS surface albedo, snow cover, normalized difference vegetation index (NDVI) products, and National Centers for Environmental Prediction (NCEP) downward solar radiation reanalysis data. The albedo data and NDVI data used in this paper are obtained from the 2001-2020 MODIS MCD43C3 and MOD13C2 data provided by NASA (https://ladsweb.modaps.eosdis.nasa.gov (accessed on 1 November 2020)). The HEG tool was adopted to extract the black sky albedo (BSA) and the white sky albedo (WSA) in the albedo data. BSA represents the albedo under complete direct solar radiation, and WSA

Dataset and Preprocessing
This study used three sets of remote sensing products and reanalysis datasets (Table 1), including MODIS surface albedo, snow cover, normalized difference vegetation index (NDVI) products, and National Centers for Environmental Prediction (NCEP) downward solar radiation reanalysis data. The albedo data and NDVI data used in this paper are obtained from the 2001-2020 MODIS MCD43C3 and MOD13C2 data provided by NASA (https://ladsweb.modaps. eosdis.nasa.gov (accessed on 1 November 2020)). The HEG tool was adopted to extract the black sky albedo (BSA) and the white sky albedo (WSA) in the albedo data. BSA represents the albedo under complete direct solar radiation, and WSA represents the albedo under complete diffusion of solar radiation. Then, MATLAB software was used to calculate the daily surface albedo, and the daily data were spatially aggregated into monthly, seasonal, and annual data. For NDVI data, considering that the areas with NDVI < 0.1 are bare soil and sparse vegetation areas, the values of NDVI > 0.1 were screened [47]. However, since this study took the MODIS MOD10CM product provided by the EOS/MODIS data center (http://nsidc.org/NASA/MODIS (accessed on 29 November 2020)) in the United States as snow cover data, a value between 0 and 100 was used to represent the snow cover rate of each pixel.

NCEP Reanalysis Products
To obtain the true surface albedo, the diffuse skylight ratio needs to be calculated from the downward solar flux reanalysis data provided by NCEP [46,48]. Meanwhile, the calculation of the regional average albedo requires surface downward radiation [34]. The data can be downloaded from the NCEP data website (https://psl.noaa.gov/data/ gridded/data.ncep.reanalysis.html (accessed on 5 December 2020)). Because the spatial resolution of the NCEP reanalysis data is inconsistent with that of the albedo data, this paper used the nearest-neighbor interpolation to interpolate all NCEP reanalysis data to the same resolution of 0.05 • as MODIS data.

Vegetation Coverage Data
The vegetation cover uses the global land cover data GlobeLand30 developed by China in 2020 with a spatial resolution of 30 m (http://www.webmap.cn/mapDataAction. do?method=globalLandCover (accessed on 3 November 2020)). The data covers the land range of 80 degrees north-south latitude, including 10 surface coverage types: cultivated land, forest, grassland, shrub land, wetland, water body, tundra, artificial surface, bare land and, glacier and snow.

Calculation of Surface Albedo
The true albedo (blue sky albedo) is approximately equal to the weighted combination of the black sky albedo and the white sky albedo [48], and the surface albedo can be used the following formula calculates: where BSA and WSA are the black sky albedo and white sky albedo, f di f is the diffuse skylight ratio, and DD v , DD n , BD v , and BD n are the visible diffuse downward solar flux, near IR diffuse downward solar flux, visible beam downward solar flux and near IR beam downward solar flux, respectively.

Calculation of Regional Average Albedo
The monthly average surface albedo of an area is calculated according to Equation (3) [34]. For Equation (3), we need to consider the downward radiation.
where α is the spatially aggregated shortwave albedo. For pixel i, A i is the area of the pixel, F i d is the surface downward radiation under all-sky condition, BSA i and WSA i are the BSA and WSA, respectively.

Trend Analysis
In this paper, the univariate linear regression method is used to estimate the interannual change rates of surface albedo, NDVI and snow cover in the QM in the past 20 years. The calculation formula of the change rate θ slope is as follows: where n is the total number of years in the study period, and α i is the mean value of the variable in year i. θ slope > 0 means that the change of the variable in n years is an increasing trend; on the contrary, θ slope < 0 means that the variable is in a decreasing trend. Then, at the confidence level of 0.05, the F-test method is used to test the significance of the change trend of each pixel.

Multi-Year Average Characteristics of Surface Albedo
First, the albedo data of the QM from 2001 to 2020 was used to calculate the monthly average data, and then the multi-year average and spring (March-May), summer (June-August), autumn (September-November), and winter (December-February of the following year) surface albedo was synthesized. The multi-year average surface albedo represents the overall albedo situation of the QM. According to the results shown in Figure 2a, the multi-year average surface albedo of the QM is about 0.25, showing obvious differences in the spatial distribution. The whole surface albedo increases with altitude, and the surface albedo at northwest is high and the southeast is low. The areas with a high-value of surface albedo are mainly distributed in the Daxue Mountain, Shulenan Mountain, Danghenan Mountain, Tergun Daban Mountain, Hark Mountain, and Lenglongling Mountain, etc., showing a northwest-southeast trend similar to the trend of the mountains. Meanwhile, the land cover is mainly desert, bare rock, glacier, and snow, and the surface albedo value is between 0.4 and 0.6. The areas with a low value of surface albedo are mainly distributed at Qinghai Lake, Har Lake, Laji Mountain, and Heihe River Basin, with values between 0.05 and 0.1. The land cover is dominated by lakes, rivers, and high vegetation coverage. Since the lake is similar to a black body and has a strong ability to absorb short waves, Qinghai Lake is the area with a minimum surface albedo.
where n is the total number of years in the study period, and is the mean value of t variable in year .
> 0 means that the change of the variable in n years is an increa ing trend; on the contrary, < 0 means that the variable is in a decreasing trend. The at the confidence level of 0.05, the F-test method is used to test the significance of t change trend of each pixel.

Multi-Year Average Characteristics of Surface Albedo
First, the albedo data of the QM from 2001 to 2020 was used to calculate the month average data, and then the multi-year average and spring (March-May), summer (Jun August), autumn (September-November), and winter (December-February of the follo ing year) surface albedo was synthesized. The multi-year average surface albedo repr sents the overall albedo situation of the QM. According to the results shown in Figure 2 the multi-year average surface albedo of the QM is about 0.25, showing obvious diffe ences in the spatial distribution. The whole surface albedo increases with altitude, and t surface albedo at northwest is high and the southeast is low. The areas with a high-val of surface albedo are mainly distributed in the Daxue Mountain, Shulenan Mounta Danghenan Mountain, Tergun Daban Mountain, Hark Mountain, and Lenglongli Mountain, etc., showing a northwest-southeast trend similar to the trend of the mou tains. Meanwhile, the land cover is mainly desert, bare rock, glacier, and snow, and t surface albedo value is between 0.4 and 0.6. The areas with a low value of surface albe are mainly distributed at Qinghai Lake, Har Lake, Laji Mountain, and Heihe River Bas with values between 0.05 and 0.1. The land cover is dominated by lakes, rivers, and hi vegetation coverage. Since the lake is similar to a black body and has a strong ability absorb short waves, Qinghai Lake is the area with a minimum surface albedo.
Then, the standard deviation of the surface albedo was used to represent the degr of abnormality and deviation of the surface albedo from the average value. It can be se from Figure 2b that the surface albedo sensitive areas in the QM are mainly distributed the alpine belts of the Tergun Daban Mountain, Hark Mountain, Qaidam Mountain, a Zhongwunong Mountain. Also, the surface albedo has been relatively stable over the pa 20 years in the Heihe River Basin, Shiyang River Basin, Laji Mountain, the vicinity of t Qinghai Lake, and the west of Qaidam Mountain with low elevations. Then, the standard deviation of the surface albedo was used to represent the degree of abnormality and deviation of the surface albedo from the average value. It can be seen from Figure 2b that the surface albedo sensitive areas in the QM are mainly distributed in the alpine belts of the Tergun Daban Mountain, Hark Mountain, Qaidam Mountain, and Zhongwunong Mountain. Also, the surface albedo has been relatively stable over the past 20 years in the Heihe River Basin, Shiyang River Basin, Laji Mountain, the vicinity of the Qinghai Lake, and the west of Qaidam Mountain with low elevations.
The spatial distribution of the multi-year average surface albedo in the QM in different seasons was further analyzed. As shown in Figure 3, the surface albedo in the QM is high in the west and low in the east in all four seasons, but there are differences in different seasons.
The average values of surface albedo in spring, summer, autumn, and winter in the QM are 0.25, 0.18, 0.24, and 0.30, respectively, and the variation curve is in a single-valley "V" shape. It is the smallest in summer and the largest in winter. The spatial distribution of surface albedo in spring and autumn is similar, but the surface albedo is slightly higher in spring. As shown in Figure 3b, the spatial distribution of surface albedo in summer is different from other seasons. The number of areas with a high surface albedo is obviously smaller than that in other seasons. In the Shulenan Mountain and Tergun Daban Mountain, the overall spatial distribution of surface albedo is high in the west and low in the east, which is consistent with the spatial distribution in other seasons. The areas with a high surface albedo in winter are similar to those in spring and autumn. They are concentrated in Daxue Mountain, Shulenan Mountain, Danghenan Mountain, and Tergun Daban Mountain in the west and Lenglongling Mountain in the east. However, the areas with a high surface albedo in winter cover a wider range. The spatial distribution of the multi-year average surface albedo in the QM in different seasons was further analyzed. As shown in Figure 3, the surface albedo in the QM is high in the west and low in the east in all four seasons, but there are differences in different seasons. The average values of surface albedo in spring, summer, autumn, and winter in the QM are 0.25, 0.18, 0.24, and 0.30, respectively, and the variation curve is in a singlevalley "V" shape. It is the smallest in summer and the largest in winter. The spatial distribution of surface albedo in spring and autumn is similar, but the surface albedo is slightly higher in spring. As shown in Figure 3b, the spatial distribution of surface albedo in summer is different from other seasons. The number of areas with a high surface albedo is obviously smaller than that in other seasons. In the Shulenan Mountain and Tergun Daban Mountain, the overall spatial distribution of surface albedo is high in the west and low in the east, which is consistent with the spatial distribution in other seasons. The areas with a high surface albedo in winter are similar to those in spring and autumn. They are concentrated in Daxue Mountain, Shulenan Mountain, Danghenan Mountain, and Tergun Daban Mountain in the west and Lenglongling Mountain in the east. However, the areas with a high surface albedo in winter cover a wider range.  Figure 4, the annual variation of the average surface albedo for many years exhibits a "U" shaped, with the "left valley slope" indicating the variation from January to July, the "valley bottom" indicating the variation from July and August, and the "right valley slope" indicating the variation from September to December. The surface albedo was the largest in January and the smallest in August, with a value of 0.31 and 0.17, respectively. It can be seen from the variation trend line that the change slopes from January to August and August to December are respectively −0.02 and 0.03. The surface albedo shows a slow downward trend from January to August and reaches the minimum value in August, and the albedo in July and August is not much different; then, it shows a strong upward trend from September to November and then slowly rises, reaching the maximum in January of the following year.  Figure 4, the annual variation of the average surface albedo for many years exhibits a "U" shaped, with the "left valley slope" indicating the variation from January to July, the "valley bottom" indicating the variation from July and August, and the "right valley slope" indicating the variation from September to December. The surface albedo was the largest in January and the smallest in August, with a value of 0.31 and 0.17, respectively. It can be seen from the variation trend line that the change slopes from January to August and August to December are respectively −0.02 and 0.03. The surface albedo shows a slow downward trend from January to August and reaches the minimum value in August, and the albedo in July and August is not much different; then, it shows The influence of the factors (i.e., NDVI and snow cover) on the annual change of sur face albedo in the QM were further analyzed. Figure 5a shows that the annual change trend of snow cover and surface albedo is consistent, and there is a significant positive correlation between the two at the 0.01 level (both sides), with a correlation coefficient o 0.96. When the snow cover gradually decreases from January to August, the surface al bedo also decreases slowly; when the snow cover increases rapidly from August to De cember, the surface albedo also increases rapidly. The annual change of NDVI shows a trend of increasing first and then decreasing, which is opposite to the variation trend o surface albedo (Figure 5b). Through Pearson correlation analysis, it was found that the annual change of surface albedo is significantly correlated with that of NDVI at the 0.01 level (both sides), with a correlation coefficient of −0.91. In autumn, the vegetation begin to turn yellow and new snowfall occurs, and the surface albedo gradually increases and peaks in January of the following year. The above analysis results show that the annua change of surface albedo is positively related to surface cover and negatively related to NDVI.

Characteristics of Interannual Variation of Surface Albedo in QM
The annual average value of surface albedo, snow cover, and NDVI in the QM from 2001 to 2020 was calculated to obtain the interannual variation trend ( Figure 6). It can be seen that the annual average surface albedo in the QM shows a slight upward trend, with a change rate of 5.0 × 10 −3 /10a. Since the surface albedo reached the minimum and maxi The influence of the factors (i.e., NDVI and snow cover) on the annual change of surface albedo in the QM were further analyzed. Figure 5a shows that the annual change trend of snow cover and surface albedo is consistent, and there is a significant positive correlation between the two at the 0.01 level (both sides), with a correlation coefficient of 0.96. When the snow cover gradually decreases from January to August, the surface albedo also decreases slowly; when the snow cover increases rapidly from August to December, the surface albedo also increases rapidly. The annual change of NDVI shows a trend of increasing first and then decreasing, which is opposite to the variation trend of surface albedo (Figure 5b). Through Pearson correlation analysis, it was found that the annual change of surface albedo is significantly correlated with that of NDVI at the 0.01 level (both sides), with a correlation coefficient of −0.91. In autumn, the vegetation begins to turn yellow and new snowfall occurs, and the surface albedo gradually increases and peaks in January of the following year. The above analysis results show that the annual change of surface albedo is positively related to surface cover and negatively related to NDVI. The influence of the factors (i.e., NDVI and snow cover) on the annual change of surface albedo in the QM were further analyzed. Figure 5a shows that the annual change trend of snow cover and surface albedo is consistent, and there is a significant positive correlation between the two at the 0.01 level (both sides), with a correlation coefficient of 0.96. When the snow cover gradually decreases from January to August, the surface albedo also decreases slowly; when the snow cover increases rapidly from August to December, the surface albedo also increases rapidly. The annual change of NDVI shows a trend of increasing first and then decreasing, which is opposite to the variation trend of surface albedo (Figure 5b). Through Pearson correlation analysis, it was found that the annual change of surface albedo is significantly correlated with that of NDVI at the 0.01 level (both sides), with a correlation coefficient of −0.91. In autumn, the vegetation begins to turn yellow and new snowfall occurs, and the surface albedo gradually increases and peaks in January of the following year. The above analysis results show that the annual change of surface albedo is positively related to surface cover and negatively related to NDVI.

Characteristics of Interannual Variation of Surface Albedo in QM
The annual average value of surface albedo, snow cover, and NDVI in the QM from 2001 to 2020 was calculated to obtain the interannual variation trend ( Figure 6). It can be seen that the annual average surface albedo in the QM shows a slight upward trend, with a change rate of 5.

Characteristics of Interannual Variation of Surface Albedo in QM
The annual average value of surface albedo, snow cover, and NDVI in the QM from 2001 to 2020 was calculated to obtain the interannual variation trend ( Figure 6). It can be seen that the annual average surface albedo in the QM shows a slight upward trend, with a change rate of 5.0 × 10 −3 /10a. Since the surface albedo reached the minimum and maximum respectively in 2001 and 2019, the annual average surface albedo of the  (Figure 6c). This result is consistent with the observation of Zhang et al. (2021) in analyzing the variation trend of NDVI in the QM during the growing season. Through Pearson correlation analysis, it was found that there is a weak positive correlation between the surface albedo and NDVI, with a correlation coefficient of 0.12. When the linear trend of the two is eliminated, the surface albedo and NDVI are negatively correlated (r = −0.02). In addition, the spatial resolution of the surface albedo and NDVI data is still low, which may not accurately describe the interannual fluctuations of the surface conditions in the QM, resulting in a weak positive correlation between the two in this region [5]. Since the interannual variation cannot fully reflect the relationship between the surface albedo and NDVI, it is necessary to further explore their relationship by using high-precision remote sensing data.

Spatial Variation Trend of Surface Albedo in QM
The univariate linear regression method was used to calculate the interannual variation trends of the annual average surface albedo, NDVI, and snow cover in the QM. Meanwhile, the F-test method was employed to test the significance of the variation trend of each pixel at a confidence level of 0.05. As shown in Figure 7, the interannual variation of the average annual surface albedo in the QM shows obvious spatial heterogeneity, and the areas with a decreasing annual surface albedo are mainly distributed in the edge of the QM, such as the Qinghai Nanshan, the Shule River Basin, and the Heihe River Basin, accounting for about 39.6% of the total area. Most regions show an increasing trend, and the areas with increased annual surface albedo are mainly distributed in Hark Mountain, Shulenan Mountain, Tulainan Mountain, Tulai Mountain, Datong Mountain and Daban Mountain, accounting for about 60.4% of the total area. Besides, the insignificant areas account for 85.7% of the total area, while the significant areas only account for 14.3% of the total area. The areas with significantly increased annual surface albedo are located in the Tulai Mountain, Datong Mountain and Daban Mountain, accounting for 8.2% of the total area, and areas with significantly decreased annual surface albedo are located in the Qinghainan Mountain, Shule River Basin, and Heihe River Basin, accounting for 6.1% of  [49]. The comparison of the interannual changes of surface albedo and snow cover indicates that the two change trends are highly consistent. Through Pearson correlation analysis, it was found that there is a significant positive correlation between the surface albedo and snow cover rate, with a correlation coefficient of 0.95. The indicates that the interannual variation of the surface albedo in the QM is largely affected by the snow cover changes.  (Figure 6c). This result is consistent with the observation of Zhang et al. (2021) in analyzing the variation trend of NDVI in the QM during the growing season. Through Pearson correlation analysis, it was found that there is a weak positive correlation between the surface albedo and NDVI, with a correlation coefficient of 0.12. When the linear trend of the two is eliminated, the surface albedo and NDVI are negatively correlated (r = −0.02). In addition, the spatial resolution of the surface albedo and NDVI data is still low, which may not accurately describe the interannual fluctuations of the surface conditions in the QM, resulting in a weak positive correlation between the two in this region [5]. Since the interannual variation cannot fully reflect the relationship between the surface albedo and NDVI, it is necessary to further explore their relationship by using high-precision remote sensing data.

Spatial Variation Trend of Surface Albedo in QM
The univariate linear regression method was used to calculate the interannual variation trends of the annual average surface albedo, NDVI, and snow cover in the QM. Meanwhile, the F-test method was employed to test the significance of the variation trend of each pixel at a confidence level of 0.05. As shown in Figure 7, the interannual variation of the average annual surface albedo in the QM shows obvious spatial heterogeneity, and the areas with a decreasing annual surface albedo are mainly distributed in the edge of the QM, such as the Qinghai Nanshan, the Shule River Basin, and the Heihe River Basin, accounting for about 39.

Spatial Variation Trend of Surface Albedo in QM
The univariate linear regression method was used to calculate the interannual variation trends of the annual average surface albedo, NDVI, and snow cover in the QM. Meanwhile, the F-test method was employed to test the significance of the variation trend of each pixel at a confidence level of 0.05. As shown in Figure 7, the interannual variation of the average annual surface albedo in the QM shows obvious spatial heterogeneity, and the areas with a decreasing annual surface albedo are mainly distributed in the edge of the QM, such as the Qinghai Nanshan, the Shule River Basin, and the Heihe River Basin, accounting for about 39

Analysis of Influencing Factors of Surface Albedo at a Spatial Scale
Section 4.3 discusses the influencing factors for the regional average interannual variation of the surface albedo. This part focuses on the analysis of the spatially influencing factors of the surface albedo variation. The surface albedo of the QM is not only sensitive to changes in vegetation and snow cover, but also has a differentiation law with changes in terrain such as altitude. To analyze the distribution characteristics of surface albedo and the influencing factors in QM with different altitude gradients, the surface albedo and the

Analysis of Influencing Factors of Surface Albedo at a Spatial Scale
Section 4.3 discusses the influencing factors for the regional average interannual variation of the surface albedo. This part focuses on the analysis of the spatially influencing factors of the surface albedo variation. The surface albedo of the QM is not only sensitive to changes in vegetation and snow cover, but also has a differentiation law with changes in terrain such as altitude. To analyze the distribution characteristics of surface albedo and the influencing factors in QM with different altitude gradients, the surface albedo and the influencing factors were divided according to the altitude of QM at an interval of 500 m (Figure 8). With the increase of altitude, the vertical distribution of surface albedo, snow cover, and vegetation cover in the QM changes significantly. A total of 81.8% of the QM are located at the altitude of 3000-4500 m. Overall, the surface albedo and snow cover increase with the increase of altitude; the NDVI shows a trend of first increase and then decrease with the increase of altitude; the snow cover dominates the spatial variation of surface albedo. Further research found that the vertical changes of the surface albedo and its influencing factors in different seasons are consistent with the average vertical changes in many years. There are significant positive correlations between albedo and snow cover and elevation. Therefore, altitude is another important influencing factor of the surface albedo variation in the QM. It is worth noting that when the altitude exceeds 4500 m, the increment rate of surface albedo in winter is not as obvious as that in other seasons. Compared with Figure 9b, it is found that the snow cover rate in winter also decreases in high altitude areas, and it is even smaller than that in other seasons. This may be the reason why the surface albedo in winter is lower than that in spring in the area above 5000 m above sea level.
Remote Sens. 2022, 14, x FOR PEER REVIEW influencing factors were divided according to the altitude of QM at an interval o (Figure 8). With the increase of altitude, the vertical distribution of surface albed cover, and vegetation cover in the QM changes significantly. A total of 81.8% of are located at the altitude of 3000-4500 m. Overall, the surface albedo and snow increase with the increase of altitude; the NDVI shows a trend of first increase a decrease with the increase of altitude; the snow cover dominates the spatial vari surface albedo. Further research found that the vertical changes of the surface albe its influencing factors in different seasons are consistent with the average vertical c in many years. There are significant positive correlations between albedo and snow and elevation. Therefore, altitude is another important influencing factor of the albedo variation in the QM. It is worth noting that when the altitude exceeds 4500 increment rate of surface albedo in winter is not as obvious as that in other season pared with Figure 9b, it is found that the snow cover rate in winter also decreases altitude areas, and it is even smaller than that in other seasons. This may be the why the surface albedo in winter is lower than that in spring in the area above above sea level.  The spatial distribution of the interannual variation trend of snow cover has a good correspondence with the surface albedo. Specifically, the significant increase in snow cover is distributed in the south of Har Lake, the south of Qinghai Lake, Tulai Mountain, Datong Mountain, and Daban Mountain (as shown in Figure 10a), and the surface albedo also shows an increasing trend in these regions. This indicates that snow cover is an influencing factor of the significant increase in surface albedo in the central QM. Comparing the interannual variation trend of surface albedo and NDVI (Figure 10b), it is found that the spatial distribution of surface albedo and the interannual variation trend of NDVI also has a good correspondence: the areas with improved vegetation coverage is almost the same as those with reduced surface albedo. This suggests that the changes in the vegetation cover in these regions is also an important factor of surface albedo changes. The spatial distribution of the interannual variation trend of snow cover has a good correspondence with the surface albedo. Specifically, the significant increase in snow cover is distributed in the south of Har Lake, the south of Qinghai Lake, Tulai Mountain, Datong Mountain, and Daban Mountain (as shown in Figure 10a), and the surface albedo also shows an increasing trend in these regions. This indicates that snow cover is an influencing factor of the significant increase in surface albedo in the central QM. Comparing the interannual variation trend of surface albedo and NDVI (Figure 10b), it is found that the spatial distribution of surface albedo and the interannual variation trend of NDVI also has a good correspondence: the areas with improved vegetation coverage is almost the same as those with reduced surface albedo. This suggests that the changes in the vegetation cover in these regions is also an important factor of surface albedo changes. The distribution of spatial correlation coefficients between surface albedo and influencing factors (Figure 11) shows that snow cover and vegetation have significant feedback on the changes of surface albedo. Surface albedo is significantly positively correlated with snow cover in most parts of the QM but negatively correlated with NDVI ( Figure 11). Specifically, the surface albedo and snow cover rate are significantly positively correlated (with a correlation coefficient higher than 0.8) in the areas including Datong Mountain, Tulai Mountain, Hark Mountain, Zhongwunong Mountain, and the west of Tergun The spatial distribution of the interannual variation trend of snow cover has a good correspondence with the surface albedo. Specifically, the significant increase in snow cover is distributed in the south of Har Lake, the south of Qinghai Lake, Tulai Mountain, Datong Mountain, and Daban Mountain (as shown in Figure 10a), and the surface albedo also shows an increasing trend in these regions. This indicates that snow cover is an influencing factor of the significant increase in surface albedo in the central QM. Comparing the interannual variation trend of surface albedo and NDVI (Figure 10b), it is found that the spatial distribution of surface albedo and the interannual variation trend of NDVI also has a good correspondence: the areas with improved vegetation coverage is almost the same as those with reduced surface albedo. This suggests that the changes in the vegetation cover in these regions is also an important factor of surface albedo changes. The distribution of spatial correlation coefficients between surface albedo and influencing factors ( Figure 11) shows that snow cover and vegetation have significant feedback on the changes of surface albedo. Surface albedo is significantly positively correlated with snow cover in most parts of the QM but negatively correlated with NDVI ( Figure 11). Specifically, the surface albedo and snow cover rate are significantly positively correlated (with a correlation coefficient higher than 0.8) in the areas including Datong Mountain, Tulai Mountain, Hark Mountain, Zhongwunong Mountain, and the west of Tergun Daban Mountain (Figure 11a), accounting for about 72% of the total area. The areas showing a negative correlation between surface albedo and NDVI are distributed in the Shule Lake Basin, Heihe River Basin, Qinghainan Mountain, etc. (Figure 11b), which correspond to the areas with significantly improved vegetation coverage. Besides, there are sporadic positive correlations near the Qaidam Basin and Qinghai Lake, and these areas may have a combined effect of vegetation and snow cover on surface albedo changes [6]. Combining The distribution of spatial correlation coefficients between surface albedo and influencing factors ( Figure 11) shows that snow cover and vegetation have significant feedback on the changes of surface albedo. Surface albedo is significantly positively correlated with snow cover in most parts of the QM but negatively correlated with NDVI ( Figure 11). Specifically, the surface albedo and snow cover rate are significantly positively correlated (with a correlation coefficient higher than 0.8) in the areas including Datong Mountain, Tulai Mountain, Hark Mountain, Zhongwunong Mountain, and the west of Tergun Daban Mountain (Figure 11a), accounting for about 72% of the total area. The areas showing a negative correlation between surface albedo and NDVI are distributed in the Shule Lake Basin, Heihe River Basin, Qinghainan Mountain, etc. (Figure 11b), which correspond to the areas with significantly improved vegetation coverage. Besides, there are sporadic positive correlations near the Qaidam Basin and Qinghai Lake, and these areas may have a combined effect of vegetation and snow cover on surface albedo changes [6]. Combining  Figures 7b and 10b, it is found that an increase in NDVI in the low vegetation coverage area in the northern QM contributes to a decrease in albedo. Figures 7b and 10b, it is found that an increase in NDVI in the low vegetation coverage area in the northern QM contributes to a decrease in albedo.

Conclusions and Discussion
In this paper, the temporal and spatial dynamic distribution of surface albedo and its influencing factors in the QM from 2001 to 2020 are analyzed by using the MODIS MCD43C3 V006 product, combined with snow cover and NDVI data. The following conclusions are drawn. the QM showed a slight upward trend, with a change rate of 5.0 × 10 −3 /10a. The fluctuation of the surface albedo in the study period was obviously more significant than that in the previous 10 years. Snow cover showed a slight increase during this period, and there was a significant positive correlation between surface albedo and snow cover rate. NDVI showed a significant upward trend, indicating that the overall vegetation was improving. Before removing the linear trend of the surface albedo and NDVI, there was a weak positive correlation between the two; when the linear trend was removed, a negative correlation was found between the two. (4) The annual and interannual variations of the surface albedo in the QM are closely related to the surface cover. The annual variation of surface albedo is "U" shaped, with the largest variation in January and the smallest variation in August, which is positively correlated with snow cover and negatively correlated with NDVI. As for interannual variation, the increase of the regional average surface albedo is significantly related to the increase of snow cover. In terms of spatial distribution, the interannual variation of the surface albedo in most areas of the QM is mainly affected by the change of snow cover. The improvement of vegetation cover in marginal areas is the main factor of the significant decrease of surface albedo in these areas.
The MODIS surface albedo maintains a good consistency with the ground observations in most time periods, and the inversion accuracy is high [33,50]. Zhang et al. (2021) found that MODIS remote sensing data has good applicability in the QM, and it can be used to study the temporal and spatial changes of snow albedo. In the plateau area, the accuracy of MODIS surface albedo inversion results is not high enough, but after manual

Conclusions and Discussion
In this paper, the temporal and spatial dynamic distribution of surface albedo and its influencing factors in the QM from 2001 to 2020 are analyzed by using the MODIS MCD43C3 V006 product, combined with snow cover and NDVI data. The following conclusions are drawn. the QM showed a slight upward trend, with a change rate of 5.0 × 10 −3 /10a. The fluctuation of the surface albedo in the study period was obviously more significant than that in the previous 10 years. Snow cover showed a slight increase during this period, and there was a significant positive correlation between surface albedo and snow cover rate. NDVI showed a significant upward trend, indicating that the overall vegetation was improving. Before removing the linear trend of the surface albedo and NDVI, there was a weak positive correlation between the two; when the linear trend was removed, a negative correlation was found between the two. (4) The annual and interannual variations of the surface albedo in the QM are closely related to the surface cover. The annual variation of surface albedo is "U" shaped, with the largest variation in January and the smallest variation in August, which is positively correlated with snow cover and negatively correlated with NDVI. As for interannual variation, the increase of the regional average surface albedo is significantly related to the increase of snow cover. In terms of spatial distribution, the interannual variation of the surface albedo in most areas of the QM is mainly affected by the change of snow cover. The improvement of vegetation cover in marginal areas is the main factor of the significant decrease of surface albedo in these areas.
The MODIS surface albedo maintains a good consistency with the ground observations in most time periods, and the inversion accuracy is high [33,50]. Zhang et al. (2021) found that MODIS remote sensing data has good applicability in the QM, and it can be used to study the temporal and spatial changes of snow albedo. In the plateau area, the accuracy of MODIS surface albedo inversion results is not high enough, but after manual analysis and elimination of errors, the root mean square error of the surface observation results can be reduced to 0.00329, which can fully meet the accuracy requirements of climate and land surface process models [51]. In addition, long-time-series surface albedo datasets remove the uncertainty of trend analysis greatly, which is crucial for evaluating surface albedo trends [8].
The Earth's surface comprehensively manifests various ground features and is complex and changeable, which makes the surface albedo vary greatly in space and time. NDVI is one of the most commonly used vegetation indices, which can accurately reflect characteristics such as vegetation cover density and growth [41,47]. When NDVI gradually increases from a small value, the surface albedo decreases rapidly. Then, as the NDVI continues to increase, the surface albedo decreases slowly. When NDVI > 0.5, the surface albedo almost does not change with the change of NDVI [52]. In addition, ice and snow are highly positively correlated with surface albedo. Usually, the albedo of old snow is 0.7-0.8, and that of new snow is 0.8-0.9 [19]. As time passes, the new snow will be granulated, the grains will be deformed and continuously densified, the particle size will increase, and the pollution will deteriorate. The relationship between the albedo with ice and snow states is as follows: fresh snow > old snow > grain snow > glacier ice > contaminated glacier ice [5,52]. In recent years, under the warming and humidification in the plateau region, the significant changes in the surface cover have caused corresponding changes in the surface albedo, and the surface albedo will react to climate change, thus affecting the regional climate system substantially [2,5,8].
The change of surface albedo in different regions is closely related to the terrain and surface coverage [53]. This study mainly focuses on the relationship between surface albedo and its influencing factors (snow cover, vegetation, and altitude). However, different vegetation cover types make different contributions to surface albedo [54], and further investigation is needed to more accurately monitor vegetation dynamics and understand the impact of surface albedo on ecosystems. In recent years, the local ecological environment damage in the QM has attracted much attention. In response to the ecological problems in the QM, the state has successively adopted a number of remediation measures to strengthen ecological environment protection and ecological restoration in the QM. The rectification and supervision of human activities such as mineral resources development, tourism activities, grazing, and water and soil resources in the QM have achieved initial results. From 2001 to 2020, the vegetation in the QM showed significant improvement, mainly distributed in the Shule River Basin, the Heihe River Basin, and the relatively low-altitude mountains in the Qinghai Lake Basin. Meanwhile, the snow cover increased significantly in the relatively high-altitude Qinghainan Mountain, Datong Mountain, Daban Mountain, and other areas. The improvement of vegetation coverage will cause the surface albedo to decrease, and through the negative feedback of the surface albedo, it will further promote the vegetation, so that the albedo will continue to decrease. The increase of ice and snow coverage will cause the surface albedo to increase, and positive feedback will be formed between ice and snow, albedo and temperature. Based on this, the area of ice and snow will be further expanded, and the albedo will continue to increase. Besides, other factors, such as soil moisture and slope aspect, also affect the surface albedo, and future studies will investigate these factors with high-resolution data. Through the understanding of the surface albedo feedback mechanism, the relationship between climate change and human activities can be deeply understood, thereby providing a scientific basis for the win-win green development of ecological livelihoods in the QM.