Nothing Special   »   [go: up one dir, main page]
Next Article in Journal
The Multiple Aperture SAR Interferometry (MAI) Technique for the Detection of Large Ground Displacement Dynamics: An Overview
Next Article in Special Issue
Direct and Lagged Effects of Spring Phenology on Net Primary Productivity in the Alpine Grasslands on the Tibetan Plateau
Previous Article in Journal
Use of Automated Change Detection and VGI Sources for Identifying and Validating Urban Land Use Change
Previous Article in Special Issue
Radar-Derived Internal Structure and Basal Roughness Characterization along a Traverse from Zhongshan Station to Dome A, East Antarctica
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Remote Sensing
Volume 12
Issue 7
10.3390/rs12071188
remotesensing-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatiotemporal Variability of Land Surface Albedo over the Tibet Plateau from 2001 to 2019

by
Xingwen Lin
1,
Jianguang Wen
2,3,*,
Qinhuo Liu
2,3,
Dongqin You
2,
Shengbiao Wu
2,3,
Dalei Hao
2,3,
Qing Xiao
2,3,
Zhaoyang Zhang
1 and
Zhenzhen Zhang
1
1
College of Geography and Environmental Sciences, Zhejiang Normal University, Hangzhou 321004, China
2
State Key Laboratory of Remote Sensing Science, Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100101, China
3
University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(7), 1188; https://doi.org/10.3390/rs12071188
Submission received: 13 February 2020 / Revised: 27 March 2020 / Accepted: 2 April 2020 / Published: 7 April 2020
(This article belongs to the Special Issue Remote Sensing of Climatic and Environmental Changes over the Antarctic, Arctic, and the Qinghai-Tibet Plateau)
Browse Figures
Graphical abstract
"> Figure 1
<p>The land cover types (<b>A</b>) and the DEM (<b>B</b>) of the Tibetan Plateau (TP); The distributions of these three types of sites following the change of land cover change (<b>C</b>) and topographic slope and aspect (<b>D</b>). The subplot in <a href="#remotesensing-12-01188-f001" class="html-fig">Figure 1</a>A shows the percentage of the area covered by different land cover types to the TP. The hollowed triangle symbols in <a href="#remotesensing-12-01188-f001" class="html-fig">Figure 1</a>B show the distribution of selected sites. The filled triangle symbols in <a href="#remotesensing-12-01188-f001" class="html-fig">Figure 1</a>B show the location of Mt. Qogir K2 (in the northwest) and Mt. Everest (in the south). The names of counties (black color) and mountainous area (blue color) are also listed in <a href="#remotesensing-12-01188-f001" class="html-fig">Figure 1</a>B. GRA = grasslands; Shb = shrublands; BC = broadleaf croplands; SAV = savannas; EBF = evergreen broadleaf forests; DBF = deciduous broadleaf forests; ENF = evergreen needleleaf forests; and UNV = unvegetated.</p> "> Figure 2
<p>The topographic features in TP. (<b>A</b>) The slope of the TP, (<b>B</b>) details of a subarea in (<b>A</b>), (<b>C</b>) the slope of a selected site and its surroundings; (<b>D</b>) the aspect of the TP, (<b>E</b>) details of a subarea in (<b>D</b>), (<b>F</b>) the aspect of a selected site and its surroundings.</p> "> Figure 3
<p>Evolution of the land surface albedo from 2001 to 2019: (<b>A</b>) the annual albedo; (<b>B</b>) the albedo variation showed by the Sen’s slope over the entire TP, (<b>C</b>) the significant land surface albedo variation displayed by the Sen’s slope, and (<b>D</b>) the Z test results.</p> "> Figure 4
<p>The albedo SD over the entire TP for the four seasons (<b>A</b>). Spring season; (<b>B</b>). Summer season; (<b>C</b>). Autumn season; (<b>D</b>). Winter season.</p> "> Figure 5
<p>The evolution of the monthly averaged albedo (<b>A</b>), land surface temperature (LST), normalized difference vegetation index (NDVI) (<b>B</b>), and the fraction of snow cover anomaly (<b>C</b>).</p> "> Figure 6
<p>The decadal variability of yearly land surface albedo (<b>A</b>), NDVI (<b>B</b>), LST (<b>C</b>), and the fraction of snow cover (<b>D</b>) for the mountainous terrain over the years 2001 to 2019.</p> "> Figure 7
<p>Surface albedo variation following the increasing of the slope at the grassland surface over mountain terrain over the years 2001 to 2019. (<b>A</b>) were located at the surface with the north aspect; (<b>B</b>) were located at the surface with the south aspect.</p> "> Figure 8
<p>The albedo variation over different land cover types at the sites, with the mean slope of 35 degrees from 2001 to 2019. (<b>A</b>). The sites at the south-facing slope. (<b>B</b>) The sites at the north-facing slope.</p> "> Figure 9
<p>The evaluation of the percentage of land cover changes, spanning the years 2001 to 2019. WAR = Water bodies, MF = Mixed Forest, SAV = Savannas, WET = permanent Wetlands, CRO = Croplands, UB = Urban and Built-up lands, CM = Cropland/Natural Vegetation Mosaics, BAR = Barren and soil, UN = Unclassified.</p> ">
Review Reports Versions Notes

Abstract

:
As an essential climate variable (ECV), land surface albedo plays an important role in the Earth surface radiation budget and regional or global climate change. The Tibetan Plateau (TP) is a sensitive environment to climate change, and understanding its albedo seasonal and inter-annual variations is thus important to help capture the climate change rules. In this paper, we analyzed the large-scale spatial patterns, temporal trends, and seasonal variability of land surface albedo overall the TP, based on the moderate resolution imaging spectroradiometer (MODIS) MCD43 albedo products from 2001 to 2019. Specifically, we assessed the correlations between the albedo anomaly and the anomalies of normalized difference vegetation index (NDVI), the fraction of snow cover (snow cover), and land surface temperature (LST). The results show that there are larger albedo variations distributed in the mountainous terrain of the TP. Approximately 10.06% of the land surface is identified to have been influenced by the significant albedo variation from the year 2001 to 2019. The yearly averaged albedo was decreased significantly at a rate of 0.0007 (Sen’s slope) over the TP. Additionally, the yearly average snow cover was decreased at a rate of 0.0756. However, the yearly average NDVI and LST were increased with slopes of 0.0004 and 0.0253 over the TP, respectively. The relative radiative forcing (RRF) caused by the land cover change (LCC) is larger than that caused by gradual albedo variation in steady land cover types. Overall, the RRF due to gradual albedo variation varied from 0.0005 to 0.0170 W/m2, and the RRF due to LCC variation varied from 0.0037 to 0.0243 W/m2 during the years 2001 to 2019. The positive RRF caused by gradual albedo variation or the LCC can strengthen the warming effects in the TP. The impact of the gradual albedo variations occurring in the steady land cover types was very low between 2001 and 2019 because the time series was short, and it therefore cannot be neglected when examining radiative forcing for a long time series regarding climate change.

Graphical Abstract">

Graphical Abstract

1. Introduction

Land surface albedo, i.e., the fraction of incident solar radiation reflected by the land surface, is a key climate-regulating variable of the Earth’s surface radiation budget [1,2]. Land surface albedo is an essential climate variable and its variations have significant impacts on climate change, in particular for the snow-covered land surfaces at mid-latitudes and high latitudes in the Northern Hemisphere [3,4,5]. According to the report of the Intergovernmental Panel on Climate Change (IPCC), the cooling effects from surface albedo following historical land-use changes are up to −0.15 Wm−2, with an uncertainty of approximately 67% [6]. Deforestation in tropical areas has been found to increase surface albedo and contribute to Earth warming through increased solar heating of the land, which is offset by strong evapotranspiration cooling through feedbacks with clouds and precipitation [1,4,7,8,9]. A full understanding of season and inter-annual variations and longer-term trends of albedo is essential for scientific research on climate change, hydrological processes, and weather forecasting [10,11,12].
Currently, the Moderate Resolution Imaging Spectroradiometer (MODIS) [13,14,15], Global Land Surface Satellite (GLASS) [16,17], GlobAlbedo [18], and QA4ECV products [19] provide long-term albedo with global coverage. Among these products, the MODIS albedo product is the most widely and efficiently used for inter-annual albedo variation analysis due to its high quality with an uncertainty of 0.02–0.05 at both snow-free and snow-cover surface [20,21,22]. Investigations have been made over the last three decades to show that albedo variation is correlated with the ecological or meteorology parameters [23], such as snow cover [24], surface temperature [25], soil moisture [26], and land cover change [6,27]. For example, land surface albedo generally showed seasonal variation over the Northern Hemisphere, with a decreasing trend in July and an increasing trend in January [23]. The opposite trend existed between surface albedo and land surface temperature (LST). This phenomenon can be illustrated by the example that albedo decreased at a rate of 0.029 ± 0.025, while the LST increased at a rate of 0.049 + 0.038 °C decade−1 during the summers of the year from 2001 to 2016 at the Queen Elizabeth Islands glaciers located in Arctic Canada [28]. Land cover change can lead to a significant change in albedo, for example, a decrease of 0.31% in albedo is reported due to the land cover change in France from 2001 to 2013 [5]. These conclusions demonstrate that the land surface albedo has strong within-year seasonal characteristics [29,30], and is being affected by land cover change [31].
The Tibetan Plateau (TP) is very vulnerable to climate change, thus it is a key area for global climate change research [32,33,34,35,36]. An important feature of the TP is its complex topography, which has a great influence on the variation of surface albedo [37]. Although albedo variation over the TP has been analyzed over the last three decades, the underlying factors (e.g., winter snowfall, snowmelt, land cover change, and land surface temperature) influencing the spatiotemporal characterization of the TP albedo over rugged terrain are not explicit [38,39,40] and are still a challenge for regional climate change research [23,30]. Particularly, the rugged terrain of the TP increases the difficulty for land surface albedo field observation and has resulted in limited in situ albedo observations [30,33,41,42]. Multiple reanalysis albedo datasets, such as NCEP/NCAR, MERRA2, and JCDAS, can support the spatial-temporal continuous distribution of albedo variation [43,44,45], however, they are too coarse (including the spatial and temporal scales) to capture the heterogeneous land surface characteristics over the TP, which is undergoing significant changes [30].
The main objective of this paper is to characterize the spatial-temporal variation of land surface albedo over the TP, based on the 500 m daily MODIS albedo product (MCD43A3), and to analyze the relationships of the anomalies between albedo and the ecological and meteorological parameter over rugged terrain, in order to understand the climate and ecological environment change of the TP. Meanwhile, we attempt to quantify the rate of albedo and radiative forcing change to explore the feedback mechanism caused by the gradual albedo variation for different land covers over the mountainous terrain of the whole TP. The organization of this paper is as follows. Firstly, the TP and the datasets used in this paper are introduced. Secondly, we explained the statistical methods for the dataset analysis. Then, the results about the albedo variation and the radiative forcing shifts are analyzed and discussed. Finally, we summarize the conclusions for this paper.

2. Study Area and Datasets

2.1. Study Area

The TP extends from subtropical to mid-latitude regions, spanning over 25° of longitude and 10° of latitude and covering 11 biomes (based on the MODIS land cover product, MCD12Q1) (Figure 1A). Known as the ‘Third Pole’, the TP has an average elevation of approximately 4000 m, with the highest elevation of more than 8000 m located in the west (e.g., Mt. Qogir K2) and south-east (e.g., Mt. Everest). The lowest elevations are distributed over south Nepal (Figure 1B).
Forty-one sites, marked with hollowed triangles in Figure 1B, were selected randomly to demonstrate albedo variation following topographic slopes, aspects, and land cover types. These sites were divided into three types, including the sites with different topographic aspects but the same land covers and topographic slopes; the sites with different topographic slopes but the same land covers and topographic aspects; and the sites with different land covers but similar topographic slopes and aspects. Detailed information about the three types of selected sites is also displayed in Figure 1C,D.
Figure 2 has shown the TP topographic features, where Figure 2A,D are the slope and aspect for the entire TP, with details of two small regions presented in subfigures (B and C for slope, E and F for aspect). The larger slopes of the TP are distributed over the northwest Karakoram Range and the southeast Himalayan Mountains, with more than 45% of pixels marked as slopes larger than 30° (Figure 2A). The counties (e.g., Naqu, Zhiduo, and Ulan) in the middle of the TP appear as gentle slopes, with the mean slope smaller than 15°. While, for the aspect, more variations are revealed in local scales, shown as dispersed and fusing textures in Figure 2D. Such textures are formed by a wide range of aspect variation concentrated in a local region as Figure 2E,F, but suppressed in the entire TP area. Topographic aspect information (Figure 2D) indicated that more than 70% of the slopes face toward the sun, with the aspects distributed between 135° and 225°. The mean surface temperatures bias for the coldest and warmest months is more than 30 °C; for example, the Himalayan’s mean surface temperatures bias of the coldest and warmest months are approximately −25 °C and 10 °C, respectively [35,46,47]. Snowfall appears at altitudes above 3000 m, especially for the mountainous regions in the Karakoram Range and Himalayan Mountains.

2.2. Datasets and Pre-processing

Both remote sensing products and reanalysis datasets were utilized in this study (Table 1), including (a) MODIS albedo, LST, the fraction of snow cover, and normalized difference vegetation index (NDVI) products; (b) the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis data; and (c) the SRTM DEM products, which were used to derive topographic slope and aspect.

2.2.1. Remote Sensing Products

The daily MCD43A3 C6 was selected from the year 2001 to 2019, which is a 500 m spatial resolution albedo product and has a daily temporal resolution. To explain the ecological controlling factors of albedo variation, daily daytime MODIS land surface temperature (MOD11A1 LST, 1 km resolution, 1 K accuracy [48,49]), 16-day MODIS NDVI product (MOD13A1, 500 m resolution [50]), yearly MODIS land cover products (MCD12Q1, 500 m resolution [51]), and the daily fraction of snow cover product (MOD10A1, 500 m resolution [52]) were prepared in this paper. The preprocessing of these products included image reprojections, resampling, and removal of invalid data points. All MODIS products were re-projected to Geographic Lat/Lon projection, WGS84 datum to match with the SRTM DEM files. To analyze the albedo variation and the relationship with the NDVI products, the daily albedo products were compiled into 16-day composites corresponding to the NDVI products and resampled to 1000 m spatial resolution using the nearest-neighbor resample method. The invalid data of remote sensing products were removed by using the quality assessment (QA) sub-datasets. The 90 m SRTM DEM products were used to calculate the slope and aspect. Moreover, the calculated slope and aspect were resampled to a 1000 m resolution by using the nearest-neighbor resample method [53].

2.2.2. ERA-Interim Reanalysis Products

The ERA-Interim is produced by the European Center for Medium-Range Weather Forecasts (ECMWF), which can be downloaded for free from the ECMWF Public Datasets web interface or the Meteorological Archival and Retrieval System (MARS) [54]. The daily ERA-Interim product was selected, spanning the years 2001 to 2019. It can provide the TOA inclined solar radiation and downward solar radiation of the TP at a spatial scale of 0.125°. These two datasets were interpolated to a 1000 m scale to match with the pixel scale of MODIS products.

3. Statistical Methods

3.1. Trend Analysis

The albedo variation was characterized based on three factors, including (1) seasonal and annual averaged albedo and the standard deviation (SD) over the TP. The months corresponding to the four seasons are spring (MAM, Mar-Apr-May), summer (JJA, Jun-Jul-Aug), autumn (SON, Sep-Oct-Nov), and winter (DJF, Dec-Jan-Feb). (2) The correlations between the albedo and ecological and meteorological parameters, such as the NDVI, snow cover, and LST; (3) decadal trend analysis over the 19 years for the TP using the Mann–Kendall (MK) and Sen’s slope (SS) estimator nonparametric tests, which were calculated as follows:
S S = M e d i a n ( x j − x i j − i ) j > i
S = ∑ i = 1 n − 1 ∑ k = i + 1 n sgn ( x k − x i )
Z = { S − 1 V a r ( S ) S > 0 0 S = 0 S + 1 V a r ( S ) S < 0 }
sgn ( x j − x i ) = { 1 i f ( x j − x i ) > 0 0 i f ( x j − x i ) = 0 − 1 i f ( x j − x i ) < 0 }
where x j and x i are the albedo values in the j t h and i t h years, respectively. A positive value of SS indicates an increasing trend, otherwise it indicates a decreasing trend. Test statistics S and Z, shown in Equations (2) and (3), were used to determine significant variation by coupling with the two-tails of Z test. The null hypothesis (H0) of the MK test, which is “no apparent trend” for the significance level of 0.05, was used; the H0 was rejected and a significant trend was detected at the situation of the absolute value |Zs| > Z(1 − α/2) (α is the confidence level, usually set to 0.01 or 0.05).

3.2. Radiative Forcing Calculation

Radiative forcing (RF) is defined as the difference in net irradiance between the incoming energy on Earth and the energy radiated back to space [55] and denotes the externally imposed perturbation in the radiative energy budget. RF describes the imbalance in the planet’s radiation budget caused by human interventions. A positive RF illustrates a warming trend to the climatic system, whereas negative radiative forcing leads to a cooling effect [5,31,56]. The RF is expressed with the following equation [55,56,57].
R F T O A = − R s T a Δ α s
where RFTOA (in W/m2) is the instantaneous RF at the top of the atmosphere (TOA), RS is the downward solar radiation at the Earth’s surface, T a is atmospheric transmittance, which is expressed as the fraction of the radiation reflected from the surface that reaches the TOA. The Δαs is the shift in surface albedo, which can be calculated by multiplying the SS and the length of time series. Considering that T a can be expressed as the ratio of RS and the solar incident radiation at the TOA level (RTOA), the equation can be rewritten as follows [9]:
R F T O A = − R s R s R T O A Δ α s
In this study, the RS and the RTOA can be downloaded from the ERA-interim daily downward surface solar radiation and the TOA incident solar radiation datasets. As the significant albedo trend did not affect the entire study area [5], the radiative forcing was weighed against the total land surface of the TP. In this paper, the area that was significantly influenced by albedo variation was selected to be weighed against the whole TP area to calculate the relative radiative forcing over rugged terrain, which can be expressed as follows:
R R F = R F T O A * A a f f e c t e d / A T o t a l
where RRF is the relative radiative forcing (in W/m2), Aaffected is the affected area, which passed the Z test with the two-tailed value of 0.05, and ATotal is the land surface area of the whole TP.

4. Results

4.1. Spatial and Seasonal Variation of Land Surface Albedo

Figure 3 shows the land surface albedo spatial variations over the TP from 2001 to 2019. The average annual albedo over the TP exhibited the larger albedo distributed at the north-west of the TP with the mean albedo larger than 0.25, and the smaller albedo at the south-east of the TP with a mean albedo lower than 0.10 (Figure 3A). The yearly average albedo in the whole TP was approximately 0.21 (Figure 3A). The biggest land surface albedo was distributed in the mountainous areas with permanent glaciers covered, such as the Kaqing, Ruoguo, Midui, and Gongzha glaciers at the Nyenchen Tanglha Mountains, and the Kyagar and Derenmang glaciers at the Hindu Kush and Karakoram mountains (Figure 3A). The lowest albedo, which is smaller than 0.05, was located in the Hengduanshan mountainous range, which is covered by dense forests during all four seasons, in the south-east of the TP.
According to Figure 3B, the albedo shows an obvious variation, spanning the years 2001 to 2019 over the whole TP. However, land surface albedo increased at a rate larger than 0.5% at the western and northern central Himalayan mountainous range, as the area was suffering from the rapid retreat of debris-covered glaciers during the observed periods [25,39,58]. The albedo at the Naqu and Zhiduo counties also showed a larger increasing trend. These places are distributed in the zonal grasslands and unvegetated surfaces (Figure 1A), and have been undergoing rapid land cover change during the years 2001 to 2019. The surface albedo showed a decreasing trend at the Lhasa and Linzhi cities, with a rate of −0.00708 yr−1 (Figure 3B). Notably, significant albedo variation also occurred at these places, especially at the edge of the mountainous area, such as the Kunlun Mountains and the eastern Nyenchen Tanglha mountains (Figure 3C,D). However, the surface albedo did not show a significant variation during 2001 to 2019 at the middle TP and the southern central Himalayan mountainous range (Figure 3C,D), as in these places, the land cover types or the permanent glaciers have been stable or have been retreating at slower rates than those at the western Himalayas [39].
Figure 4 displays the spatial distributions of the seasonal and annual standard deviations of albedo over the TP from 2001 to 2019. The seasonal albedo varied with the elevation, land surface heterogeneity, and land cover changes, especially over the snow-covered land surface with high elevation in complex terrain. From the four subplots of Figure 4, we can see that the larger SD of albedo was at the mountainous regions during the spring, summer, autumn, and winter (such as the Hindu Kush and Karakoram mountains, western Himalayan mountains, Hengduanshan mountains, and the middle Kunlun Mountains). Notably, the Qaidam basin also had a larger SD of albedo in spring, autumn, and winter, with the SD larger than 0.10, especially in winter (Figure 4D), as this location was extremely continental, with long, cold, dry winters and short summers. A small SD of albedo occurred at the southern Hengduanshan and the south-eastern Himalayan with values smaller than 0.02 during these four seasons (Figure 4). These places were covered by dense forests (evergreen broadleaf forest (EBF), deciduous broadleaf forest (DBF), and mixed forest (MF)). The land covers were changed slowly in these places, spanning the years 2001 to 2019. The seasonal land cover change has a larger influence on albedo variation, with the SD larger than 0.1 at a larger slope surface (Figure 4A,C,D). In spring and winter, the SD of albedo peaks over the Karakoram Range, the Hengduanshan mountainous range, and the western Himalaya range, with standard deviation, was larger than 0.14 (Figure 4A,D). In autumn, there was a larger SD of albedo at the Kunlun mountainous range, the Danggula mountainous range, and the Qaidam basin, with values larger than 0.12 (Figure 4C). In summer, the surface albedo does not display large variation over the entire TP, with a low SD of 0.02, except for the Danggula mountainous range (Figure 4B).

4.2. Albedo Anomaly and the Interrelation with Ecological-Meteorological Parameters

The temporal variation illustrated a potential seasonal dependency of the correlation between the albedo anomalies and the LST, NDVI, and the fraction of snow cover anomalies over the whole TP, respectively (Figure 5A–C). The albedo and LST anomalies showed opposite variations from 2001 to 2019, with a negative correlation coefficient (Pearson’s coefficient, r) of −0.5787 (Figure 5A). The albedo anomalies and NDVI anomalies showed a weak negative correlation, with an r-value of −0.2139 (Figure 5B). However, the albedo anomalies showed a high positive correlation, with the fraction of snow covers from 2001 to 2019 having an r-value as high as 0.6086 (Figure 5C). Generally, the maximum albedo anomalies were found in the spring of 2008, with values of 0.048 and 0.050, which were as the same as the fractions of snow cover and LST anomalies. The fraction of snow cover showed a positive anomaly from the years 2018 to 2019; this means that the fraction of snow cover in these periods has little more increase than in other years. Additionally, this can be illustrated by the albedo anomaly in the winter of 2018 and 2019. Figure 5B also showed the larger NDVI variation was delayed to the summer of 2009. This is because the larger fraction of monthly snow cover in the spring and winter of 2008 may lead to an increase of monthly averaged albedo and then result in a decrease of LST. However, once the summer is coming, the snow melts rapidly with the increase of LST, and this leads to sufficient water during the growing season, which is the benefit to the vegetation growth, and this then resulted in the larger NDVI variation in the year 2009.
Decadal variations of the albedo, NDVI, LST and the fraction of snow cover were shown at the surfaces, with a mean slope larger than 15° (Figure 6). Generally, the yearly albedo was distributed in the range of 0.21 to 0.25. The yearly albedo has a slightly decreasing trend, spanning 2001 to 2019 with a Sen’s slope of −0.0007 and a p-value of 0.0021 (Figure 6A), which means that there is a significant albedo variation during these years over the TP. The yearly fraction of snow cover also showed a significant decreasing trend from 2001 to 2019, with a slope of −0.0756 and a p-value of 0.0078 (Figure 6D). Additionally, the decadal variability of yearly NDVI showed a significant increase trend during these years, with a Sen’s slope of 0.0004 and a p-value of 0.0051 (Figure 6B). It is distributed in the range of 0.11 to 0.135, with a larger value of 0.134 in 2015 and a small value of 0.119 in 2008 (Figure 6B). The yearly LST was distributed in the range of 8 to 13 °C, with a larger LST of 12.83 °C in 2015, and a small LST of 8.23 °C in 2005. The yearly LST showed a significant increasing trend over the whole TP, with a slope of −0.0253 and a p-value of 0.0078 (Figure 6C).

4.3. Albedo Variation Following the Change of Mean Slope, Aspect, and Land Covers

Figure 7 shows the long time series of land surface albedo at the sites with different mean slopes in grasslands (shown as hollowed triangles in Figure 1). It shows that the largest albedo was in winter when the land surface is snow for almost the whole season, with a value larger than 0.5 (Figure 7). A heavier snow cover situation happened over the slope far away from the sun (called the north slope) than for the slope facing toward the sun (called the south slope), particularly in the spring and winter. For the pixels on the south slope, the grassland surface showed a larger monthly albedo value of 0.24 at the sites, with a mean slope of 8.2 degrees and a smaller albedo value of 0.15 at the sites, with a mean slope of 42 degrees. The albedo recorded at the sites distributed in the south aspects also showed a similar variation compared with the north-facing surface (Figure 7B). Following the increase of the mean slope, the snow-free land surface albedo shows a decreasing trend, whether there is a change of mean topographic aspects or not.
Figure 8 explains the evolution of the land surface albedo at the sites, with a slope of 35 degrees and covered by different land cover types. Generally, the land surface albedo showed the seasonal variation from the years 2001 to 2019, with larger albedo occurring in winter and smaller ones occurring in summer (Figure 8A). The albedos recorded at the sites covered by the grassland (GRA) and unvegetated (UNV) land were larger than those at the sites covered by other land covers, with the mean albedo larger than 0.15 at the north slopes (Figure 8). The north-facing slope suffered from more snowfall than that at the south-facing slope during the winter and spring seasons. Additionally, the snow persists longer at the north-facing slope than at the south-facing slope. This phenomenon can be detected indirectly from the frequency of the larger albedo that occurred (e.g., larger than 0.4) at the grasslands or the savannas (Figure 8A). Except for the snow-covered periods, the land surface albedo does not show differences at the sites located in the opposite topographic aspects. Both the results of Figure 7 and Figure 8 demonstrated that the topography slope and land cover types have a large influence on land surface albedo over the mountainous ranges.

4.4. The Radiative Forcing Shifts Due to the Albedo Variation over the Topographic Terrain

The albedo shifts (delta albedo) of the black-sky albedo (BSA) and white-sky albedo (WSA) are used to calculate the radiative forcing (RF) and the relative radiative forcing (RRF), based on Equations (6) and (7). Additionally, the RF and RRF were explored in the areas, which were found to be rarely influenced by land cover changes (LCC), to only assess the contribution of gradual albedo variations from 2001 to 2019 (see Table 2). The nine types of land cover were selected in this scenario, based on the MCD12Q1 products type 2 dataset (IGBP classification). Furthermore, the barren soil, urban and built surface were combined with the non-vegetated surface in order to characterize the radiative forcing variation affected by albedo shift on the vegetated surface and non-vegetated surface. Generally, the mean daily surface solar radiation downwards (SSRD) was 232.25 W/m2. The mean summer daily TOA incident solar radiation (TIRS) was 467.36 W/m2, recorded at the reanalysis ERA-interim datasets over the TP from 2001 to 2019. Table 2 showed that a larger negative albedo shift occurred at the GRA, with a BSA shift of −0.0040 and a WSA shift of −0.0038. A smaller negative albedo shift occurred at unvegetated (UNV) surfaces, with a BSA and WSA shift of approximately −0.0019. The negative BSA and WSA shift showed the positive radiative forcing for overall land cover type, with a larger radiative forcing of 0.4639 (by using the BSA shift) and 0.4408 W/m2 (by using the WSA shift) occurring at the GRA, and a small radiative forcing of 0.2140 W/m2 for the UNV land surface. In terms of the RF shift, the BSA and WSA variations have a weaker impact on a non-vegetated surface than on a vegetated-covered surface (Table 2).
The surface albedo affecting the RF was weighted against the surface to evaluate the real impacts on radiative forcing. The pixels were masked as having sensible effects on the albedo variation at the rugged terrain from 2001 to 2019, which passed the check of the Z test (as the two-sided p-value was <0.05). Generally, 10.06% of the pixels were detected as having sensible effects by albedo variation over the whole TP. Approximately 4.61% of the pixels covered by UNV at the TP were detected to suffer from the significant albedo variation from 2001 to 2019, which occupied the largest significantly changed area of the TP compared with the pixels covered by other land cover types. On the contrary, only 0.21% of the pixels covered by the EBF were detected to have been significantly influenced by albedo variation, which occupied the smallest changed area compared with that at other land cover types (Table 2). Considering that the albedo variation must be weighted against the significantly influenced surface of the mountainous terrain to evaluate the real impacts on radiative forcing, the RRF was calculated, as shown in Table 2. The larger RRF due to BSA and WSA trends occurring in GRA were 0.0170 and 0.0162 W/m2, respectively. The smaller RRF occurred at the ENF surface, with values of 0.0006 and 0.0005 W/m2 for BSA and WSA, respectively. Indeed, the impact of gradual albedo trends occurring in the steady land cover types on the RF was weak compared to that due to land cover change.

5. Discussion

The land cover change (LCC) can partly explain the albedo, LST, NDVI, and radiative forcing variation during these years. To analyze the LCC, the land cover types which were recorded in the MCD12Q1 products type 2 dataset (IGBP classification) were re-divided into 13 types, by combining the EBF and the DBF into the broadleaf forest (BF), combining the ENF and the DNF into the needle leaf forest (NF), combining the open shrubland and the closed shrubland into Shb, and combining the woody savanna and savanna into SAV. Generally, the percentage of the land surface covered by the 13 land cover types showed slight variations over the whole TP (Figure 9). The percentage of barren and soil (BAR) to the whole TP was decreased significantly from the year 2001 to 2019, with an SS of −0.0877 and a p-value smaller than 0.05. In contrast, the area of Shb, Water, and SAV were increased significantly during the observed periods (Figure 9). However, the areas covered by other land cover types (e.g., the GRA, BF, NF, MF, CRO, WET, CM, UN, and UB) displayed insignificant increasing/decreasing trends from the year 2001 to 2019.
Statistics results also demonstrated that about 4.36% of the land surfaces covered by barren and soil were affected by LCC, with about 1.74% changing to savannas, 0.47% changing to shrublands, 0.29% changing to water bodies, and 1.85% changed to other land cover types (Table 2). The mutual change of LCC from a barren and soil zone to a vegetated-covered area resulted in yearly albedo having a decreasing trend during these years, based on the fact that the vegetated-covered land surface (e.g., the shrublands and savannas) has a smaller albedo than that at the barren and soil in the TP. The average albedo of barren and soil is approximately 0.23, while the albedo of the vegetated-covered area is approximately 0.13–0.19 in the TP (Table 3). The LCC led to a radiative forcing variation ranging from 0.7953 W/m2 to 5.0712 W/m2, depending on the LCC type. Using the WSA, the LCC led to RF variation ranging from 0.9220 W/m2 to 5.3017 W/m2, for Shb and water bodies, respectively. The positive radiative forcing was able to enhance the warming effect. Furthermore, the decreasing trend (Figure 6D) of snow cover also increases the absorption of solar shortwave radiation and then enhances atmospheric warming [59]. The warming trend leads to more snow melting and snow-line ascension, which forms positive feedback to climate change over the TP [36] and results in the melting of the glacial deposits at a much faster rate [59] following a significant increase of water-covered area (Table 3). The LCCs, such as the snow cover changes, which include snow melting or snow falling during winter or summer seasons over the whole mountain region, directly and quickly affected the surface albedo [28]. These phenomena reduce the albedo and increase the greenness of the land surface and hence trap more solar heat flux. The workflows formed a positive feedback mechanism for snow melting in the mountain area over the TP [60,61].
The variation of albedo caused by LCC from barren and soil to vegetation results in the higher sensitivity of albedo to NDVI, because the changes of canopy cover will rapidly change the contributions of the vegetation and barren soil albedo to the overall surface albedo. Over the whole TP, the significant LCC from barren soil to the vegetated-covered surface (Figure 9) has a good consistency with the results in other studies, such as the work of Chen [62], Peng [63] and Song [64]. These researchers showed slightly increasing trends in the forest-covered area over China and India from 2000 to 2015, which can support the results above [62,63,64,65]. The p-value of 0.05 over the south-east area of the Himalayas emphasizes the significant albedo change (Figure 3C,D), which was affected by land cover change due to afforestation over the whole Tibetan Plateau over the last 30 years [46].
The albedo variation caused by the land cover changes at the regional scale may lead to large radiative forcing and result in regional climate change. Surface radiative forcing (RF) can be used to measure and evaluate the contributing factors (e.g., greenhouse gases, land cover changes, and the albedo) to climate change [66,67]. Generally, the land cover change provides a greater contribution to the radiative forcing than do other factors [5,68,69].

6. Conclusions

The Tibetan Plateau (TP) is very vulnerable to climate change and plays an important role in regional weather and climate research in East and South Asia [46,70,71]. The land surface albedo, LST, NDVI, fraction of snow cover, and the radiative forcing over the mountainous terrain of the TP were selected to understand the regional radiation transfer and climate change.
(1) The mountainous terrain was more sensitive to albedo variation than the places with gentle slopes, which can be identified as the significant albedo anomaly occurred at the mountainous terrain, with larger topographic slopes over the years 2001 to 2019 (Figure 3C,D). The land surface albedo and albedo anomaly exhibited the spatial and seasonal variability at the Tibetan Plateau. The larger yearly albedo was distributed in the southwest of the Hindu Kush, the western and central Himalayas, and the central Kunlun mountainous range, where many permanent glaciers are distributed and where they were suffering from rapid retreating during the observed time period (Figure 3).
(2) The topographic slope has a larger influence on the surface albedo than the topographic aspect, with the decreasing trend of albedo following the increase of mean slope. However, the topographic aspects did not show any obvious differences in the albedo variation when comparing the albedo at the south-facing slope to the albedo at the north-facing slope.
(3) The yearly albedo and the fraction of snow cover have a significant decreasing trend from 2001 to 2019. However, the yearly LST and NDVI showed an increasing trend during this period (Figure 6). The slow albedo variations impact on RF. The RF due to LCC was shown to be larger than that due to gradual albedo variation in steady land cover types. Overall, the positive RF due to gradual albedo variation in steady land cover types, which were detected to not land cover changes during the years from 2001 to 2019, can reach up to 0.0005~0.4408 W/m2. The RF due to LCC can reach up to 5.0712 W/m2 during these years. Weighed against the surface area of the study area, approximately 10.06% of the area suffered from significant albedo variation over the whole TP due to gradual albedo variation. Nevertheless, the magnitude of RRF due to gradual albedo variation decreased by 90% than that the magnitude of RF, and still has the warming effects from the year 2001 to 2019. The LCC also led to a positive RF over the TP, since the land cover was primarily changed from barren and soil to savanna, shrubland, and water bodies, which have smaller albedos than the barren and soil. The RRF due to LCC varied from about 0.0037 to 0.0243 W/m2, when using the BSA shift, and 0.0043 to 0.0239 W/m2 using the WSA shift.
In conclusion, the albedo variation needs to be taken into account in radiative forcing research for long duration climate change studies, whatever the causes of the albedo trends. The investigation of the NDVI, albedo, snow cover, and LST gradual variation reveals that the gradual decreases in snow cover had a positive impact on RF via the decreasing of the albedo, which might lead to warming effects. However, the gradual increase in the NDVI showed a slight increase in forest cover and gave rise to additional evapotranspiration and surface roughness, which feeds back to climate cooling. The gradual increases of LST reveal that the decrease of snow cover has a more significant influence than other parameters on climate change. However, the mechanisms linking the surface relative parameters and climate change are very complex at temperate latitudes, and therefore more relative parameters (e.g., evapotranspiration and surface roughness) need to be measured and analyzed to support climate change research.

Author Contributions

Conceptualization, X.L. and J.W.; Methodology, X.L. and D.Y.; Software, S.W. and D.H.; Formal Analysis, X.L. and J.W.; Investigation, Z.Z. (Zhaoyang Zhang) and X.L.; Resources, Z.Z. (Zhenzhen Zhang); Data Curation, D.H.; Writing-Original Draft Preparation, X.L. and J.W.; Writing-Review & Editing, X.L., J.W. and D.Y.; Visualization, X.L. and J.W.; Supervision, Q.L. and Q.X.; Project Administration, X.L.; Funding Acquisition, X.L. and D.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was jointly supported by the Chinese Natural Science Foundation Project (NO.: 41901290), “Strategic Priority Research Program” of the Chinese Academy of Sciences (Grant No. XDA15012102), Open Fund of State Key Laboratory of Remote Sensing Science (Grant No. OFSLRSS201907).

Acknowledgments

The authors would like to acknowledge the National Aeronautics and Space Administration (NASA) for providing MODIS MCD43A3, MCD12Q1, MOD13A1, MOD11A1 and MOD10A1 products and ASTER STRM products, the European Centre for Medium-Range Weather Forecasts (ECMWF) for offering a platform to download ERA-Interim datasets for TOA incident solar radiation and downward solar radiation data. Special thanks are given to the anonymous reviewers for their fruitful comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bonan, G.B. Forests and Climate Change: Forcings, Feedbacks, and the Climate Benefits of Forests. Nature 2008, 320, 1444–1449. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Kiehl, J.T.; Hack, J.J.; Bonan, G.B.; Boville, B.A.; Briegleb, B.P.; Williamson, D.L.; Rasch, P.J. Description of the Ncar Community Climate Model (ccm3). Available online: https://www.osti.gov/biblio/442361 (accessed on 13 February 2020).
  3. Anderson, R.G.; Canadell, J.G.; Randerson, J.T.; Jackson, R.B.; Hungate, B.A.; Baldocchi, D.D.; Ban-Weiss, G.A.; Bonan, G.B.; Caldeira, K.; Cao, L. Biophysical considerations in forestry for climate protection. Front. Ecol. Environ. 2011, 9, 174–182. [Google Scholar] [CrossRef] [Green Version]
  4. Bala, G.; Caldeira, K.; Wickett, M.; Phillips, T.J.; Lobell, D.B.; Delire, C.; Mirin, A. Combined climate and carbon-cycle effects of large-scale deforestation. Proc. Natl. Acad. Sci. USA 2007, 104, 6550–6555. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Planque, C.; Carrer, D.; Roujean, J.-L. Analysis of MODIS albedo changes over steady woody covers in France during the period of 2001–2013. Remote Sens. Environ. 2017, 191, 13–29. [Google Scholar] [CrossRef]
  6. Bauer, S.E.; Menon, S. Aerosol direct, indirect, semidirect, and surface albedo effects from sector contributions based on the IPCC AR5 emissions for preindustrial and present-day conditions. J. Geophys. Res. Atmos. 2012, 117–130. [Google Scholar] [CrossRef]
  7. Bonan, G.B. Effects of land use on the climate of the United States. J. Clim. Chang. 1997, 37, 449–486. [Google Scholar] [CrossRef]
  8. Pimentel, R.; Aguilar, C.; Herrero, J.; Pérezpalazón, M.J.; Polo, M.J. Comparison between Snow Albedo Obtained from Landsat TM, ETM+ Imagery and the SPOT VEGETATION Albedo Product in a Mediterranean Mountainous Site. Hydrology 2016, 3, 10. [Google Scholar] [CrossRef] [Green Version]
  9. Jiao, T.; Williams, C.A.; Ghimire, B.; Masek, J.; Gao, F.; Schaaf, C. Global Climate Forcing from Albedo Change Caused by Large-scale Deforestation and Reforestation: Quantification and Attribution of Geographic Variation. Clim. Chang. 2017, 142, 463–476. [Google Scholar] [CrossRef]
  10. Cess, R.D. Biosphere-Albedo Feedback and Climate Modeling. J. Atmos. Sci. 1978, 35, 1765–1768. [Google Scholar] [CrossRef]
  11. Van de Kerchove, R.; Lhermitte, S.; Veraverbeke, S.; Goossens, R. Spatio-temporal variability in remotely sensed land surface temperature, and its relationship with physiographic variables in the Russian Altay; Mountains. Int. J. Appl. Earth Obs. 2013, 20, 4–19. [Google Scholar] [CrossRef] [Green Version]
  12. Ollinger, S.V.; Richardson, A.D.; Martin, M.E.; Hollinger, D.Y.; Frolking, S.E.; Reich, P.B.; Plourde, L.C.; Katul, G.G.; Munger, J.W.; Oren, R. Canopy nitrogen, carbon assimilation, and albedo in temperate and boreal forests: Functional relations and potential climate feedbacks. Proc. Natl. Acad. Sci. USA 2008, 105, 19336. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Lucht, W.; Hyman, A.H.; Strahler, A.H.; Barnsley, M.J.; Hobson, P.; Muller, J.P. A comparison of satellite-derived spectral albedos to ground-based broadband albedo measurements modeled to satellite spatial scale for a semidesert landscape. Remote Sens. Environ. 2000, 74, 85–98. [Google Scholar] [CrossRef]
  14. Schaaf, C.B.; Gao, F.; Strahler, A.H.; Lucht, W.; Li, X.W.; Tsang, T.; Strugnell, N.C.; Zhang, X.Y.; Jin, Y.F.; Muller, J.P.; et al. First operational BRDF, albedo nadir reflectance products from MODIS. Remote Sens. Environ. 2002, 83, 135–148. [Google Scholar] [CrossRef] [Green Version]
  15. Strahler, A.H.; Muller, J.; Lucht, W.; Schaaf, C.; Tsang, T.; Gao, F.; Li, X.; Lewis, P.; Barnsley, M.J. MODIS BRDF/albedo product: Algorithm theoretical basis document version 5.0. MODIS Doc. 1999, 23, 42–47. [Google Scholar]
  16. Liu, Q.; Wang, L.Z.; Qu, Y.; Liu, N.F.; Liu, S.H.; Tang, H.R.; Liang, S.L. Preliminary evaluation of the long-term GLASS albedo product. Int. J. Digit. Earth 2013, 6, 69–95. [Google Scholar] [CrossRef]
  17. Liang, S.; Zhao, X.; Liu, S.; Yuan, W.; Cheng, X.; Xiao, Z.; Zhang, X.; Liu, Q.; Cheng, J.; Tang, H.; et al. A long-term Global LAnd Surface Satellite (GLASS) data-set for environmental studies. Int. J. Digit. Earth 2013, 6, 5–33. [Google Scholar] [CrossRef]
  18. Muller, J.-P.; López, G.; Watson, G.; Shane, N.; Kennedy, T.; Yuen, P.; Lewis, P.; Fischer, J.; Guanter, L.; Domench, C. The ESA GlobAlbedo Project for Mapping the Earth’s Land Surface Albedo for 15 Years from European Sensors. Available online: https://www.semanticscholar.org/paper/The-ESA-GlobAlbedo-Project-for-mapping-the-Earth%27s-Muller-Lopez/5c67e348414afec872360c8d1d5063a57a1ab722 (accessed on 13 February 2020).
  19. Folkert Boersma, K.; Muller, J.P. Quality Assurance for Essential Climate Variables. In Proceedings of the Egu General Assembly Conference, Vienna, Austria, 23–28 April 2017. [Google Scholar]
  20. Cescatti, A.; Marcolla, B.; Vannan, S.K.S.; Pan, J.Y.; Román, M.O.; Yang, X.; Ciais, P.; Cook, R.B.; Law, B.E.; Matteucci, G. Intercomparison of MODIS albedo retrievals and in situ measurements across the global FLUXNET network. Remote Sens. Environ. 2012, 121, 323–334. [Google Scholar] [CrossRef] [Green Version]
  21. Mira, M.; Weiss, M.; Baret, F.; Courault, D.; Hagolle, O.; Gallego-Elvira, B.; Olioso, A. The MODIS (collection V006) BRDF/albedo product MCD43D: Temporal course evaluated over agricultural landscape. Remote Sens. Environ. 2015, 170, 216–228. [Google Scholar] [CrossRef] [Green Version]
  22. Wang, Z.; Schaaf, C.B.; Sun, Q.; Shuai, Y.; Román, M.O. Capturing rapid land surface dynamics with Collection V006 MODIS BRDF/NBAR/Albedo (MCD43) products. Remote Sens. Environ. 2018, 207, 50–64. [Google Scholar] [CrossRef]
  23. He, T.; Liang, S.; Song, D.X. Analysis of global land surface albedo climatology and spatial-temporal variation during 1981–2010 from multiple satellite products. J. Geophys. Res. Atmos. 2014, 119, 10281–10298. [Google Scholar] [CrossRef]
  24. Loranty, M.M.; Berner, L.T.; Goetz, S.J.; Jin, Y.; Randerson, J.T. Vegetation controls on northern high latitude snow-albedo feedback: Observations and CMIP5 model simulations. Glob. Chang. Biol. 2014, 20, 594–606. [Google Scholar] [CrossRef] [PubMed]
  25. Shekhar, M.; Chand, H.; Kumar, S.; Srinivasan, K.; Ganju, A. Climate-change studies in the western Himalaya. Ann. Glaciol. 2010, 51, 105–112. [Google Scholar] [CrossRef] [Green Version]
  26. Verheijen, F.G.A.; Jeffery, S.; van der Velde, M.; Penízek, V.; Beland, M.; Bastos, A.C.; Keizer, J.J. Reductions in soil surface albedo as a function of biochar application rate: Implications for global radiative forcing. Environ. Res. Lett. 2013, 8, 4008. [Google Scholar] [CrossRef] [Green Version]
  27. Ghimire, B.; Williams, C.A.; Masek, J.; Feng, G.; Wang, Z.; Schaaf, C.; Tao, H. Global albedo change and radiative cooling from anthropogenic land-cover change, 1700 to 2005 based on MODIS, land-use harmonization, radiative kernels and reanalysis. Geophys. Res. Lett. 2014, 41, 9087–9096. [Google Scholar] [CrossRef]
  28. Mortimer, C.A.; Sharp, M. Spatiotemporal variability of Canadian High Arctic glacier surface albedo from MODIS data, 2001–2016. Cryosphere 2018, 12, 701–720. [Google Scholar] [CrossRef] [Green Version]
  29. Fang, H.; Liang, S.; Kim, H.Y.; Townshend, J.R.; Schaaf, C.L.; Strahler, A.H.; Dickinson, R.E. Developing a spatially continuous 1 km surface albedo data set over North America from Terra MODIS products. J. Geophys. Res. Atmos. 2007, 112, 1–20. [Google Scholar] [CrossRef] [Green Version]
  30. Shi, Q.; Liang, S. Characterizing the surface radiation budget over the Tibetan Plateau with ground-measured, reanalysis, and remote sensing data sets: 2. Spatiotemporal analysis: SRB OF THE TIBETAN PLATEAU, 2. J. Geophys. Res. Atmos. 2013, 118, 8921–8934. [Google Scholar] [CrossRef]
  31. Shine, K.P. Radiative Forcing of Climate Change. Weather 2002, 57, 278–283. [Google Scholar]
  32. Beniston, M. Climatic Change in Mountain Regions: A Review of Possible Impacts. Clim. Chang. 2003, 59, 5–31. [Google Scholar] [CrossRef]
  33. Ma, Y.; Fan, S.; Ishikawa, H.; Tsukamoto, O.; Yao, T.; Koike, T.; Zuo, H.; Hu, Z.; Su, Z. Diurnal and inter-monthly variation of land surface heat fluxes over the central Tibetan Plateau area. Theor. Appl. Climatol. 2005, 80, 259–273. [Google Scholar] [CrossRef]
  34. Shen, M.; Piao, S.; Jeong, S.-J.; Zhou, L.; Zeng, Z.; Ciais, P.; Chen, D.; Huang, M.; Jin, C.-S.; Li, L.Z.X.; et al. Evaporative cooling over the Tibetan Plateau induced by vegetation growth. Proc. Natl. Acad. Sci. USA 2015, 112, 9299–9304. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Tian, L.; Zhang, Y.; Zhu, J. Decreased surface albedo driven by denser vegetation on the Tibetan Plateau. Environ. Res. Lett. 2014, 9, 104001. [Google Scholar] [CrossRef]
  36. Yao, T.; Xue, Y.; Chen, D.; Chen, F.; Thompson, L.; Cui, P.; Koike, T.; Lau, W.K.-M.; Lettenmaier, D.; Mosbrugger, V.; et al. Recent Third Pole’s Rapid Warming Accompanies Cryospheric Melt and Water Cycle Intensification and Interactions between Monsoon and Environment: Multidisciplinary Approach with Observations, Modeling, and Analysis. Bull. Am. Meteorol. Soc. 2019, 100, 423–444. [Google Scholar] [CrossRef]
  37. Wen, J.; Liu, Q.; Xiao, Q.; Liu, Q.; You, D.; Hao, D.; Wu, S.; Lin, X. Characterizing Land Surface Anisotropic Reflectance over Rugged Terrain: A Review of Concepts and Recent Developments. Remote Sens. 2018, 10, 370. [Google Scholar] [CrossRef] [Green Version]
  38. Zheng, L.; Zhao, G.; Dong, J.; Ge, Q.; Tao, J.; Zhang, X.; Qi, Y.; Doughty, R.B.; Xiao, X. Spatial, temporal, and spectral variations in albedo due to vegetation changes in China’s grasslands. ISPRS J. Photogramm. Remote Sens. 2019, 152, 1–12. [Google Scholar] [CrossRef]
  39. Scherler, D.; Bookhagen, B.; Strecker, M.R. Spatially variable response of Himalayan glaciers to climate change affected by debris cover. Nat. Geosci. 2011, 4, 156. [Google Scholar] [CrossRef]
  40. Xie, X. A Climatic Estimation of Surface Albedo in Summer and the Distribution of Albedo over the Qinghai-Xizang Plateau; American Meteorological Society: Boston, MA, USA; pp. 162–169.
  41. Guo, X.; Yang, K.; Long, Z.; Wei, Y.; Li, S.; Zhu, M.; Yao, T.; Chen, Y. Critical Evaluation of Scalar Roughness Length Parametrizations Over a Melting Valley Glacier. Bound. Layer Meteorol. 2011, 139, 307–332. [Google Scholar] [CrossRef]
  42. Yang, W.; Guo, X.; Yao, T.; Yang, K.; Zhao, L.; Li, S.; Zhu, M. Summertime surface energy budget and ablation modeling in the ablation zone of a maritime Tibetan glacier. J. Geophys. Res. Atmos. 2011, 116, 1–11. [Google Scholar] [CrossRef]
  43. Balsamo, G.; Albergel, C.; Beljaars, A.; Boussetta, S.; Brun, E.; Cloke, H.; Dee, D.; Dutra, E.; Munoz-Sabater, J.; Papenberger, F. ERA-Interim/Land: A global land surface reanalysis dataset. Hydrol. Earth Syst. Sci. 2015, 19, 389–407. [Google Scholar] [CrossRef] [Green Version]
  44. Bao, X.; Zhang, F. Evaluation of NCEP-CFSR, NCEP-NCAR, ERA-Interim, and ERA-40 reanalysis datasets against independent sounding observations over the Tibetan Plateau. J. Clim. 2013, 26, 206–214. [Google Scholar] [CrossRef] [Green Version]
  45. Dee, D.P.; Uppala, S.M.; Simmons, A.; Berrisford, P.; Poli, P.; Kobayashi, S.; Andrae, U.; Balmaseda, M.; Balsamo, G.; Bauer, D.P. The ERA-Interim reanalysis: Configuration and performance of the data assimilation system. Q. J. R. Meteorol. Soc. 2011, 137, 553–597. [Google Scholar] [CrossRef]
  46. Cui, X.; Graf, H.F. Recent land cover changes on the Tibetan Plateau: A review. Clim. Chang. 2009, 94, 47–61. [Google Scholar] [CrossRef] [Green Version]
  47. Tanaka, K.; Tamagawa, I.; Ishikawa, H.; Ma, Y.; Hu, Z. Surface energy budget and closure of the eastern Tibetan Plateau during the GAME-Tibet IOP 1998. J. Hydrol. 2003, 283, 169–183. [Google Scholar] [CrossRef]
  48. Coll, C.; Caselles, V.; Galve, J.M.; Valor, E.; Niclòs, R.; Sánchez, J.M.; Rivas, R. Ground measurements for the validation of land surface temperatures derived from AATSR and MODIS data. Remote Sens. Environ. 2005, 97, 288–300. [Google Scholar] [CrossRef]
  49. Wan, Z. New refinements and validation of the MODIS Land-Surface Temperature/Emissivity products. Remote Sens. Environ. 2008, 112, 59–74. [Google Scholar] [CrossRef]
  50. Fensholt, R.; Sandholt, I.; Rasmussen, M.S. Evaluation of MODIS LAI, fAPAR and the relation between fAPAR and NDVI in a semi-arid environment using in situ measurements. Remote Sens. Environ. 2004, 91, 490–507. [Google Scholar] [CrossRef]
  51. Sulla-Menashe, D.; Friedl, M.A. User Guide to Collection 6 MODIS Land Cover (MCD12Q1 and MCD12C1) Product; USGS: Reston, VA, USA, 2018.
  52. Stroeve, J.C.; Box, J.E.; Haran, T. Evaluation of the MODIS (MOD10A1) daily snow albedo product over the Greenland ice sheet. Remote Sens. Environ. 2006, 105, 155–171. [Google Scholar] [CrossRef]
  53. Wu, S.; Wen, J.; You, D.; Zhang, H.; Xiao, Q.; Liu, Q. Algorithms for Calculating Topographic Parameters and Their Uncertainties in Downward Surface Solar Radiation (DSSR) Estimation. IEEE Geosci. Remote Sens. Lett. 2018, 15, 1149–1153. [Google Scholar] [CrossRef]
  54. Berrisford, P.; Dee, D.P.; Fielding, K.; Fuentes, M.; Kållberg, P.W.; Kobayashi, S.; Uppala, S.M. The ERA-interim archive. ERA Rep. 2009, 1, 1–16. [Google Scholar]
  55. Ramaswamy, V.; Boucher, O.; Haigh, J.; Hauglustine, D.; Haywood, J.; Myhre, G.; Nakajima, T.; Shi, G.; Solomon, S. Radiative forcing of climate. Clim. Chang. 2001, 349–416. [Google Scholar]
  56. Muñoz, I.; Campra, P.; Fernández-Alba, A.R. Including CO2-emission equivalence of changes in land surface albedo in life cycle assessment. Methodology and case study on greenhouse agriculture. Int. J. Life Cycle Assess. 2010, 15, 672–681. [Google Scholar] [CrossRef]
  57. Lenton, T.M.; Vaughan, N.E. Physics. The radiative forcing potential of different climate geoengineering options. Atmos. Chem. 2009, 9, 5539–5561. [Google Scholar] [CrossRef] [Green Version]
  58. Li, Z.; He, Y.; An, W.; Song, L.; Zhang, W.; Catto, N.; Wang, Y.; Wang, S.; Liu, H.; Cao, W. Climate and glacier change in southwestern China during the past several decades. Environ. Res. Lett. 2011, 6, 045404. [Google Scholar] [CrossRef]
  59. Rangwala, I.; Miller, J.R.; Ming, X. Warming in the Tibetan Plateau: Possible influences of the changes in surface water vapor. Geophys. Res. Lett. 2009, 36, 295–311. [Google Scholar] [CrossRef] [Green Version]
  60. Dimri, A.P.; Kumar, D.; Choudhary, A.; Maharana, P. Future changes over the Himalayas: Mean temperature. Glob. Planet. Chang. 2018, 162, 212–234. [Google Scholar] [CrossRef]
  61. Betts, R.A. Offset of the potential carbon sink from boreal forestation by decreases in surface albedo. Nature 2000, 408, 187–190. [Google Scholar] [CrossRef]
  62. Chen, C.; Park, T.; Wang, X.; Piao, S.; Xu, B.; Chaturvedi, R.K.; Fuchs, R.; Brovkin, V.; Ciais, P.; Fensholt, R.; et al. China and India lead in greening of the world through land-use management. Nat. Sustain. 2019, 2, 122–129. [Google Scholar] [CrossRef]
  63. Jian, P.; Liu, Z.; Liu, Y.; Wu, J.; Han, Y. Trend analysis of vegetation dynamics in Qinghai–Tibet Plateau using Hurst Exponent. Ecol. Indic. 2012, 14, 28–39. [Google Scholar]
  64. Yi, S.; Long, J.; Wang, H. Vegetation Changes along the Qinghai-Tibet Plateau Engineering Corridor Since 2000 Induced by Climate Change and Human Activities. Remote Sens. 2018, 10, 95. [Google Scholar]
  65. Deng, H.; Pepin, N.C.; Liu, Q.; Chen, Y. Understanding the spatial differences in terrestrial water storage variations in the Tibetan Plateau from 2002 to 2016. Clim. Chang. 2018, 151, 379–393. [Google Scholar] [CrossRef] [Green Version]
  66. Andrews, T.; Forster, P.M.; Gregory, J.M. A surface energy perspective on climate change. J. Clim. 2009, 22, 2557. [Google Scholar] [CrossRef] [Green Version]
  67. Zhai, J.; Liu, R.; Liu, J.; Zhao, G.; Huang, L. Radiative forcing over China due to albedo change caused by land cover change during 1990–2010. Acta Geogr. Sin. 2014, 24, 789–801. [Google Scholar] [CrossRef] [Green Version]
  68. Mahmood, R.; Pielke, R.A., Sr.; Hubbard, K.G.; Niyogi, D.; Bonan, G.; Lawrence, P.; Mcnider, R.; Mcalpine, C.; Etter, A.; Gameda, S. Impacts of land use/land cover change on climate and future research priorities. Bull. Am. Meteorol. Soc. 2010, 91, 37–46. [Google Scholar] [CrossRef]
  69. Pielke, R.A. Land Use and Climate Change. Science 2005, 310, 1625–1626. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Manabe, S.; Broccoli, A.J. Mountains and Arid Climates of Middle Latitudes. Sci. China Earth Sci. 1990, 247, 192–195. [Google Scholar] [CrossRef]
  71. Ye, D.Z.; Wu, G.X. The role of the heat source of the Tibetan Plateau in the general circulation. Meteorol. Atmos. Phys. 1998, 67, 181–198. [Google Scholar] [CrossRef]
Remotesensing 12 01188 g001 550
Figure 1. The land cover types (A) and the DEM (B) of the Tibetan Plateau (TP); The distributions of these three types of sites following the change of land cover change (C) and topographic slope and aspect (D). The subplot in Figure 1A shows the percentage of the area covered by different land cover types to the TP. The hollowed triangle symbols in Figure 1B show the distribution of selected sites. The filled triangle symbols in Figure 1B show the location of Mt. Qogir K2 (in the northwest) and Mt. Everest (in the south). The names of counties (black color) and mountainous area (blue color) are also listed in Figure 1B. GRA = grasslands; Shb = shrublands; BC = broadleaf croplands; SAV = savannas; EBF = evergreen broadleaf forests; DBF = deciduous broadleaf forests; ENF = evergreen needleleaf forests; and UNV = unvegetated.
Figure 1. The land cover types (A) and the DEM (B) of the Tibetan Plateau (TP); The distributions of these three types of sites following the change of land cover change (C) and topographic slope and aspect (D). The subplot in Figure 1A shows the percentage of the area covered by different land cover types to the TP. The hollowed triangle symbols in Figure 1B show the distribution of selected sites. The filled triangle symbols in Figure 1B show the location of Mt. Qogir K2 (in the northwest) and Mt. Everest (in the south). The names of counties (black color) and mountainous area (blue color) are also listed in Figure 1B. GRA = grasslands; Shb = shrublands; BC = broadleaf croplands; SAV = savannas; EBF = evergreen broadleaf forests; DBF = deciduous broadleaf forests; ENF = evergreen needleleaf forests; and UNV = unvegetated.
Remotesensing 12 01188 g001
Remotesensing 12 01188 g002 550
Figure 2. The topographic features in TP. (A) The slope of the TP, (B) details of a subarea in (A), (C) the slope of a selected site and its surroundings; (D) the aspect of the TP, (E) details of a subarea in (D), (F) the aspect of a selected site and its surroundings.
Figure 2. The topographic features in TP. (A) The slope of the TP, (B) details of a subarea in (A), (C) the slope of a selected site and its surroundings; (D) the aspect of the TP, (E) details of a subarea in (D), (F) the aspect of a selected site and its surroundings.
Remotesensing 12 01188 g002
Remotesensing 12 01188 g003 550
Figure 3. Evolution of the land surface albedo from 2001 to 2019: (A) the annual albedo; (B) the albedo variation showed by the Sen’s slope over the entire TP, (C) the significant land surface albedo variation displayed by the Sen’s slope, and (D) the Z test results.
Figure 3. Evolution of the land surface albedo from 2001 to 2019: (A) the annual albedo; (B) the albedo variation showed by the Sen’s slope over the entire TP, (C) the significant land surface albedo variation displayed by the Sen’s slope, and (D) the Z test results.
Remotesensing 12 01188 g003
Remotesensing 12 01188 g004 550
Figure 4. The albedo SD over the entire TP for the four seasons (A). Spring season; (B). Summer season; (C). Autumn season; (D). Winter season.
Figure 4. The albedo SD over the entire TP for the four seasons (A). Spring season; (B). Summer season; (C). Autumn season; (D). Winter season.
Remotesensing 12 01188 g004
Remotesensing 12 01188 g005 550
Figure 5. The evolution of the monthly averaged albedo (A), land surface temperature (LST), normalized difference vegetation index (NDVI) (B), and the fraction of snow cover anomaly (C).
Figure 5. The evolution of the monthly averaged albedo (A), land surface temperature (LST), normalized difference vegetation index (NDVI) (B), and the fraction of snow cover anomaly (C).
Remotesensing 12 01188 g005
Remotesensing 12 01188 g006 550
Figure 6. The decadal variability of yearly land surface albedo (A), NDVI (B), LST (C), and the fraction of snow cover (D) for the mountainous terrain over the years 2001 to 2019.
Figure 6. The decadal variability of yearly land surface albedo (A), NDVI (B), LST (C), and the fraction of snow cover (D) for the mountainous terrain over the years 2001 to 2019.
Remotesensing 12 01188 g006
Remotesensing 12 01188 g007 550
Figure 7. Surface albedo variation following the increasing of the slope at the grassland surface over mountain terrain over the years 2001 to 2019. (A) were located at the surface with the north aspect; (B) were located at the surface with the south aspect.
Figure 7. Surface albedo variation following the increasing of the slope at the grassland surface over mountain terrain over the years 2001 to 2019. (A) were located at the surface with the north aspect; (B) were located at the surface with the south aspect.
Remotesensing 12 01188 g007
Remotesensing 12 01188 g008 550
Figure 8. The albedo variation over different land cover types at the sites, with the mean slope of 35 degrees from 2001 to 2019. (A). The sites at the south-facing slope. (B) The sites at the north-facing slope.
Figure 8. The albedo variation over different land cover types at the sites, with the mean slope of 35 degrees from 2001 to 2019. (A). The sites at the south-facing slope. (B) The sites at the north-facing slope.
Remotesensing 12 01188 g008
Remotesensing 12 01188 g009 550
Figure 9. The evaluation of the percentage of land cover changes, spanning the years 2001 to 2019. WAR = Water bodies, MF = Mixed Forest, SAV = Savannas, WET = permanent Wetlands, CRO = Croplands, UB = Urban and Built-up lands, CM = Cropland/Natural Vegetation Mosaics, BAR = Barren and soil, UN = Unclassified.
Figure 9. The evaluation of the percentage of land cover changes, spanning the years 2001 to 2019. WAR = Water bodies, MF = Mixed Forest, SAV = Savannas, WET = permanent Wetlands, CRO = Croplands, UB = Urban and Built-up lands, CM = Cropland/Natural Vegetation Mosaics, BAR = Barren and soil, UN = Unclassified.
Remotesensing 12 01188 g009
Table
Table 1. The datasets used in this study.
Table 1. The datasets used in this study.
ParameterDatasetSpatial ResolutionTemporal Resolution
AlbedoMCD43A3500 mDaily
Normalized difference vegetation index (NDVI)MOD13A1500 m16-day
Land surface temperature (LST)MOD11A11000 mDaily
Land Cover TypeMCD12Q1500 mYearly
The fraction of snow coverMOD10A1500 mDaily
Top of Atmosphere (TOA) incident solar radiationERA-Interim0.125°Daily
Downward solar radiationERA -Interim0.125°Daily
Table
Table 2. Radiative forcing (RF) and Relative radiative forcing (RRF) induced by the black-sky albedo (BSA) and white-sky albedo (WSA) trend over mountainous terrain at TP from 2001 to 2019.
Table 2. Radiative forcing (RF) and Relative radiative forcing (RRF) induced by the black-sky albedo (BSA) and white-sky albedo (WSA) trend over mountainous terrain at TP from 2001 to 2019.
Land CoverDelta BSADelta WSAAffected AreaRF BSA
(W·m−2)		RF WSA
(W·m−2)		RRF BSA
(W·m−2)		RRF WSA
(W·m−2)
GRA−0.0040−0.00383.67%0.46390.44080.01700.0162
Shb−0.0028−0.00300.24%0.32270.34950.00080.0008
BC−0.0037−0.00380.26%0.42620.43290.00110.0011
SAV−0.0030−0.00310.27%0.34660.35940.00090.0010
EBF−0.0031−0.00300.21%0.36180.34110.00080.0007
DBF−0.0024−0.00250.24%0.27810.28690.00070.0007
ENF−0.0022−0.00200.24%0.24990.22960.00060.0005
DNF−0.0019−0.00220.32%0.21830.25000.00070.0008
UNV−0.0019−0.00204.61%0.21400.23430.00990.0108
Table
Table 3. Surface albedo variation and radiative forcing affected by land cover changes in the TP between 2000 and 2019; BSA* = average BSA, WSA* = average WSA; Area LCC = surface area affected by the land cover change (LCC) from 2001 to 2019; Delta BSA* = the BSA variation affected by the LCC; Delta WSA* = the WSA variation affected by the LCC.
Table 3. Surface albedo variation and radiative forcing affected by land cover changes in the TP between 2000 and 2019; BSA* = average BSA, WSA* = average WSA; Area LCC = surface area affected by the land cover change (LCC) from 2001 to 2019; Delta BSA* = the BSA variation affected by the LCC; Delta WSA* = the WSA variation affected by the LCC.
Land CoverLCCAreaLCCBSA*WSA*Delta BSA*Delta WSA*RF BSARF WSARRF BSARRF WSA
Barren and SoilSAV1.74%0.17160.1815−0.0108−0.01121.24481.29090.02170.0225
Shb0.47%0.13880.1474−0.0069−0.00800.79530.92200.00370.0043
WAT0.29%0.05910.0621−0.0440−0.04605.07125.30170.01470.0154
Other0.44%0.15400.1628−0.0114−0.01121.31391.29090.02430.0239

Share and Cite

MDPI and ACS Style

Lin, X.; Wen, J.; Liu, Q.; You, D.; Wu, S.; Hao, D.; Xiao, Q.; Zhang, Z.; Zhang, Z. Spatiotemporal Variability of Land Surface Albedo over the Tibet Plateau from 2001 to 2019. Remote Sens. 2020, 12, 1188. https://doi.org/10.3390/rs12071188

AMA Style

Lin X, Wen J, Liu Q, You D, Wu S, Hao D, Xiao Q, Zhang Z, Zhang Z. Spatiotemporal Variability of Land Surface Albedo over the Tibet Plateau from 2001 to 2019. Remote Sensing. 2020; 12(7):1188. https://doi.org/10.3390/rs12071188

Chicago/Turabian Style

Lin, Xingwen, Jianguang Wen, Qinhuo Liu, Dongqin You, Shengbiao Wu, Dalei Hao, Qing Xiao, Zhaoyang Zhang, and Zhenzhen Zhang. 2020. "Spatiotemporal Variability of Land Surface Albedo over the Tibet Plateau from 2001 to 2019" Remote Sensing 12, no. 7: 1188. https://doi.org/10.3390/rs12071188

APA Style

Lin, X., Wen, J., Liu, Q., You, D., Wu, S., Hao, D., Xiao, Q., Zhang, Z., & Zhang, Z. (2020). Spatiotemporal Variability of Land Surface Albedo over the Tibet Plateau from 2001 to 2019. Remote Sensing, 12(7), 1188. https://doi.org/10.3390/rs12071188

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop