Wetlands Cool Land Surface Temperature in Tropical Regions but Warm in Boreal Regions

: Wetlands play a critical role in global hydrological and biogeochemical cycles. Regulating the regional climate is one of the most important ecosystem services of natural wetlands. However, the impact of wetlands on local temperature on the global scale and the attribution is still unclear. This study utilizes the satellite-based products (land surface temperature (LST), albedo, and evapotranspiration (ET)) to evaluate the difference in LST between wetlands and their adjacent landcover types and the possible drivers. Here we show that on average for the whole year, wetlands have a cooling effect in tropical regions, but have a warming effect in boreal regions. The impacts of wetlands on LST show great seasonality in the boreal regions; i.e., the wetlands have a warming effect in winter but a cooling effect in summer. The difference in albedo and ET between wetlands and the other landcover types only interprets 30% of temporal variation of the difference in LST. Due to the large water storage in wetlands, the ground heat ﬂux (G) may interpret the rest of the impact, absorbing energy in summer and releasing energy in winter in wetlands, which has often been neglected in previous studies. Our results indicate that it is critical to comprehensively consider the effects of wetland restoration in different regions to realize potential climatic beneﬁts in the future.


Introduction
Wetlands cover only 10% (13.0 million km 2 ) of the world's land surface area [1], but exert critical impacts on global hydrological and biogeochemical cycles, as well as biodiversity conservation [2][3][4][5]. Due to climate change and the intensifying human activities, however, both satellite images and wetland inventories suggest a decline in the area of natural wetlands around the world [6][7][8]. For instance, the Global Wetland Outlook 2018 reported an approximately 35% decrease of inland and marine/coastal wetlands between 1970 and 2015 [4]. In turn, some essential ecosystem services of natural wetlands such as water purification, flood control, and climate regulation could be jeopardized. As the threat of climate change has grown more severe, the feedbacks between climate and wetlands could be intensified and cause unexpected consequences. Therefore, identifying and understanding the ecosystem services of natural wetlands is critical to sustainable development goals [9]. Regulating the regional climate is one of the most important ecosystem services of natural wetlands [10]. According to the radiation budget equation on an ecosystem scale, the net radiation, R net , i.e., the balance between the inputs and outputs of shortwave (K in and K out ) and longwave radiation (L in and L out ), is approximately balanced by the energy that is transferred out of the ecosystem by non-radiative processes including sensible heat flux (H), latent heat flux (LE), ground heat flux (G), and the change in chemical energy stored in the ecosystem (∆S) [11]: R net = K in − K out + L in − L out = H + LE + G + ∆S With dark water surface and dense vegetation, wetlands tend to have a lower albedo and thus receive more net radiation in contrast to adjacent drylands (the incoming solar radiation is expected to be very comparable for all grids due to proximity) [12]. On the other hand, the existence of vegetation and water surface in wetland grids could result in significant differences in the transformation of available energy into latent and sensible heat fluxes (i.e., LE and H, respectively) as well as surface roughness against ambient other landscapes. If the ground heat flux G is always negligible, the wetland will have a cooler surface if more net radiation is dissipated as LE, compared with other landscapes (∆S is only <5% of net shortwave radiation [13]); otherwise, the wetland will have a warmer surface. Most station-based studies suggested a surface cooling effect for wetlands, probably due to the increased latent heat when comparing the surface temperature in wetlands and other adjacent landscapes [14,15]. For example, Bai et al. [16] indicated that soil temperature is higher in an area surrounded with less wetland. Zheng et al. [17] found that the land surface temperature (LST) in the Minjiang River estuary increased during the wetland deterioration period (1993-2013), indicating a cooling effect of wetlands without deterioration [17]. Given that, in summer and autumn, more heat could be moved away from the wetland surface due to the large heat capacity of the water body and turbulent mixing capability of water/vapor fluxes, the ground heat flux G in wetlands is substantial and plays an important role on LST in these two seasons. Yet the ground heat flux G has been neglected in most existing studies about the impacts of land-use change on LST [18,19], and also in the evaluation of the impacts of wetlands on local/regional climate. Moreover, the attribution of the climatic impacts of wetlands on the global scale is still scarce and needs to be investigated.
In this study, we use a satellite-derived LST product during the period 2001-2017 to verify the impacts of wetlands on climate. The differences of albedo and evapotranspiration (ET) between wetlands and adjacent other landscapes (including forest, shrubland, and open land) are adopted to investigate the radiative and nonradiative impacts of wetlands, respectively. The primary objectives of this investigation are to (1) characterize the annual and seasonal climatic impacts of global natural wetlands, and (2) interpret the net climatic impacts of wetlands from the radiative (i.e., albedo) and nonradiative (i.e., H, LE, and G) biophysical aspects.

Data
In this study, the LST products and landcover maps derived from MODIS (Moderate Resolution Imaging Spectroradiometer) with a spatial resolution of 1 × 1 km were utilized to investigate the annual and seasonal climatic impacts of global natural wetlands. Specifically, the difference in LST between wetland pixels and other landcover types (∆LST) was estimated for each 0.5 × 0.5 • grid cell. Following that, the albedo and ET data from MODIS were analyzed to interpret the ∆LST between wetland pixels and other landscapes, from radiative and nonradiative aspects, respectively. More details for these data and relevant analyses are provided below.

Land Use and Landcover Data
The MODIS land use and landcover map (MCD12Q1) in 2001, with an overall accuracy of 75% from the cross-validation accuracy assessment [20], were used to distinguish wetland, forest, shrubland, and open land pixels. This dataset encompasses 17 landcover classes proposed by the International Geosphere-Biosphere Program (IGBP), of which class 11, class 1-5, class 6-9, and class 10, 12, 14 represent permanent wetlands, forest, and open lands, respectively. Open lands include cropland, grassland, or a mix of cropland and grassland. The absolute area and the fraction of the land surface area of the 17 classes are summarized by the statistics in Table 1. According to the IGBP classification, only permanently inundated lands with 30-60% water cover and >10% vegetated cover were identified as wetlands, such that in the following analysis global natural wetlands covered only 0.73% (1.1 million km 2 ) of the land surface area (see Table 1), which is much less than those in other studies where the seasonally inundated wetlands are included [1]. Note that the area of permanent wetlands in 2001 was the lowest during the period 2001-2017. We used the landcover map of 2001 as the minimum area of wetlands. We also tested that our results were robust with landcover maps of other years.

. LST Data
The LST is the radiative skin temperature of the land surface, as measured in the direction of the remote sensor, and LST is a mixture of vegetation and bare soil temperatures [21]. The LST data were obtained from the MODIS Earth Observing System Terra instruments, which are sun-synchronous and image the entire Earth every 1-2 days. Thus, a daytime (local solar time~10:30 a.m. from Terra) and a nighttime (~10:30 p.m. from Terra) temperature observation at a 1 × 1 km resolution could be provided every 8 days. This product has been widely used at more than 42 sites in Greenland, China, Europe, and Africa and is proven to be of high quality [22], with errors less than 1 • C [23]. In our analysis, the daytime and nighttime LST were averaged to obtain the daily mean values during the period 2001-2017.

Albedo Data
The MODIS black-sky albedo (or directional hemispherical reflectance) (MCD43D51), a concept of albedo in the absence of a diffuse component with shortwave broadband (0.3-5.0 µm) of solar radiation, was utilized because it can describe the energy difference between different landcover types. In MODIS, black sky albedo products are available from 2001 to 2017 with a spatial resolution of 30 s (1 km) and a time resolution of 8-day intervals. The uncertainty of albedo data is up to 0.02 (10%) [24], and has been extensively validated in typical landcover types [25].

ET Data
The MODIS ET data (MOD16A2), 8-day composite data available from 2001 to 2017, is utilized here because it has a spatial resolution of 500 m, which is much finer compared to other ET products. The data was derived from the daily meteorological reanalysis products and 8-day remotely sensed dynamics of vegetation characteristics; thus, it was independent of the MODIS LST product. The mean absolute errors between ET estimates in this dataset and observations from ground-based flux measurements are 0.33 kg m −3 day −1 [26]. It should be noted that the MOD16 evaporation algorithm has the best performance for forests while it has mismatches for grasslands [27].

Air Temperature, Radiation, and Elevation Data
The CERES-EBAF shortwave fluxes at 1 • resolution in the all-sky scene from the Clouds and the Earth's Radiant Energy System (CERES) project website (https://ceres. larc.nasa.gov/data/, accessed on 10 May 2019) were used to calculate energy from solar shortwave radiation. We also used the monthly air temperature (Tair) data during the period 2001-2017 from the Climate Research Unit (CRU) time-series v4.02 climate data [28] to compare with the LST with a spatial resolution of 0.5 × 0.5 • . Furthermore, to minimize the impact of elevation on ∆LST, the SRTM (Shuttle Radar Topography Mission) elevation data at 1 arc-second (30 m) resolution were used to remove grid cells with elevation difference between wetlands, forests, shrublands, and open lands larger than 100 m.

Methods
All datasets are in the WGS84 geographic coordinate system. First, we split the world into regular 0.5 × 0.5 • grid cells. In each grid cell, we calculated the mean LST, albedo and ET of all 1 km or 500 m pixels of each landcover type. Then, the ∆LST, ∆albedo and ∆ET between wetlands and other adjacent landcover types were derived for each 0.5 × 0.5 • grid. The land covers in the same 0.5 × 0.5 • grid cell were regarded as adjacent landscapes. We defined the 0.5 × 0.5 • grid cells with wetland fraction larger than 1% in the land use and landcover map as "wetlands" grids, and grids with a wetland fraction smaller than 1% were masked out. Then we removed grid cells with elevation difference between wetlands, forests, shrublands, and open lands larger than 100 m. Finally, only 2335, 2336, and 1599 sample grids were obtained for the comparison of ∆LST between wetlands and forests, shrublands, and open lands, respectively. Globally, permanent wetlands are mainly located in northwest Canada, West Siberia, Amazon, Congo Basin, and the Malay Archipelago ( Figure 1). The global mean annual air temperature, LST, and albedo of these permanent wetlands are 9.7 ± 14.4 • C, 7.8 ± 14.6 • C, and 0.2 ± 0.1, respectively (Figures S1 and S2). Wetlands in different regions may have different spatial patterns of ∆LST, which can be governed by different underlying mechanisms. The typical regions were selected to characterize a class/type of underlying mechanisms for corresponding spatial patterns of ∆LST ( Figure 1). Comparing the difference in surface energy fluxes between wetlands and other adjacent landcover types, Equation (1) can be reformulated as follows: where ∆K and ∆L are the differences of shortwave and longwave radiation, respectively; ∆R net is the difference in net radiation; ∆H, ∆LE, and ∆G are the differences in sensible, latent, and ground heat flux respectively; and ∆S' is the difference in the change of chemical energy storage. In the calculation of the difference in surface energy fluxes between wetlands and the other adjacent landcover types, ∆K (the difference of shortwave radiation) was calculated by multiplying ∆albedo and the energy flux of solar shortwave radiation, and ∆LE (the difference of latent heat flux) was calculated by multiplying ∆ET and latent heat of vaporization at 2.45 MJ/kg. The sum of effects of ∆albedo and ∆ET is noted as ∆SEF, equal to the sum of ∆K and ∆LE, which are two major energy fluxes that impact ∆LST. For the effects of other energy fluxes (∆L, ∆H, ∆G, and ∆S') on ∆LST, we discussed these effects in Section 4.2. To attribute the effect of ∆albedo and ∆LE to the temporal variation of ∆LST, we performed a dominance analysis based on the multiple linear regression model [29].
The dominance analysis has been widely used for quantifying the relative importance of predictors based on the decomposition of R 2 [30]. In our study, the dominance analysis was conducted with the "relaimpo" package in R to derive the relative importance of ∆albedo and ∆ET for the temporal variation of monthly ∆LST. Figure 2 illustrates the spatial patterns of mean annual ∆LST between wetlands and three ambient landcover types, denoted by ∆LST w f , ∆LST w s , and ∆LST w o , respectively. As shown in Figure 2, although the global mean values of ∆LST w f , ∆LST w s , and ∆LST w o were relatively small (0.0 ± 0.6, −0.1 ± 0.6, and −0.1 ± 0.7 • C, respectively) compared to the difference between LST and air temperature in wetlands (∆LST w a , −2.1 ± 1.3 • C), regional and latitudinal differences were found in the ∆LST w f , ∆LST w s , and ∆LST w o . In boreal regions (>60 • N, Northeast Canada and West Siberia), more than 55% of the regions (except the south part of West Siberia) had positive ∆LST (with mean values of 0.0 ± 0.2, 0.1 ± 0.2, and 0.1 ± 0.2 • C for ∆LST w f , ∆LST w s , and ∆LST w o , respectively), indicating that the LST of wetlands was higher than those in forests, shrublands, and open lands. In tropical regions (30 • S-30 • N, Amazon, Congo Basin, and the Malay Archipelago), more than 70% of the region (most in the Amazon Basin and Malaysia) had negative ∆LST (−0.3 ± 0.9 • C), indicating that the LST of wetlands was lower than those in the three other landcover types. Therefore, the general spatial/latitudinal patterns of ∆LST indicated a warming effect of natural wetlands in boreal regions but a cooling effect in the low-latitude tropical regions.

Spatial Patterns of Seasonal ∆LST
The spatial patterns of ∆LST (i.e., ∆LST w f , ∆LST w s , and ∆LST w o ) were found to have significant differences as the seasons shifted. Because of the similarity embedded in spatial patterns in spring  Figure S3 and the statistics of ∆LST in four typical regions is shown in Table S1. In tropical regions, about 70% of wetland grids had negative ∆LST w s , and ∆LST w o (i.e., −0.2 ± 0.6 and −0.1 ± 0.6 • C) but positive ∆LST w f (0.3 ± 0.6 • C) across all seasons. In summer, most of the wetlands in the boreal latitude regions had negative ∆LST, indicating that the LST in wetlands was lower than those in the other three landcover types; i.e., wetlands had a cooling effect. In winter, about 70% of wetlands in the high-latitude regions (mostly in Canada) had positive ∆LST w s (0.2 ± 0.5 • C) and ∆LST w o (0.1 ± 0.5 • C) but negative ∆LST w f (−0.2 ± 0.7 • C), indicating that LST in wetlands was higher than those in shrublands and open lands but lower than those in forests; i.e., wetlands had a warming effect over shrublands and open lands but a cooling effect over forests. The winter (and autumn) ∆LST between wetlands and other landcover types in the boreal region dominated the signal of mean annual ∆LST.

Spatial Patterns of Seasonal ∆albedo and ∆ET
To illustrate the radiative and nonradiative impacts of wetlands on LST, the spatial patterns of seasonal ∆albedo and ∆ET between wetlands and the other three landcover types were further investigated (Figures 4 and 5). Figure 4 shows the spatial patterns of ∆albedo of JJA and DJF between wetlands and forests (∆albedo w f , Figure 4a Figure S4. In tropical regions, ∆albedo (i.e., ∆albedo w f , ∆albedo w s and ∆albedo w o were all negative (i.e., −0.02 ± 0.02, −0.03 ± 0.03, and −0.03 ± 0.02, respectively) throughout the year, indicating that wetlands had a smaller albedo and thus reflected less solar radiation than the other three landcover types. However, in the boreal regions, ∆albedo had different signals between summer and winter. In summer, the mean of ∆albedo was close to zero (i.e., 0.00 ± 0.02, −0.01 ± 0.02, and −0.01 ± 0.02) for ∆albedo w f , ∆albedo w s and ∆albedo w o , respectively), while in winter, the ∆albedo values were all positive (i.e., 0.12 ± 0.1, 0.02 ± 0.04, and 0.08 ± 0.07 for ∆albedo w f , ∆albedo w s and ∆albedo w o respectively) (Figure 4b,d,f). This indicates that wetlands have a higher albedo in winter due to snow cover [18,31] and thus reflect more solar radiation than the other three landcover types.     Figure S5. Similar to ∆albedo, in tropical regions, ∆ET had an almost consistent pattern throughout the year, while in the boreal regions ∆ET had distinct patterns between summer and winter. Specifically, in tropical regions, ∆ET had similar values (i.e., 0.00 ± 0.20, 0.00 ± 0.25, and 0.01 ± 0.28 mm/day for ∆ET w f , ∆ET w s and ∆ET w o , respectively) and had little seasonality. In boreal regions, although ∆ET had negative values throughout the year; the absolute values of ∆ET decreased as the season shifted from summer (0.00 ± 0.18 mm/day) to winter (0.01 ± 0.02 mm/day).

The Cooling or Warming Effect of Wetlands
The global pattern of ∆LST between wetlands and other landcover types provides an important insight that the overall behavior of wetlands throughout the year has a warming effect in boreal regions but a cooling effect in the tropical regions. Many local scale studies located in different regions confirm the global spatial patterns in ∆LST found in our results. For example, in tropical regions, the similar cooling effect (−1.7 • C) of wetlands is observed in the Minjiang River (25 • N-26 • N, 119 • E-120 • E) by comparing the satellite-based LST over wetlands and other adjacent landscapes from 1993-2013 [32]. In the Kilombero catchment (7 • S-10 • S, 34 • E-37 • E), from 2000 to 2017, the shrinking of wetlands increased the MODIS LST of the corresponding grid cell by 2-3 • C [33]. In temperate regions, in a three-year (2015-2017) observation of the eddy flux tower, it was found that the annual aerodynamic surface temperature of wetlands was up to 0.2 • C warmer than that of adjacent croplands in the Sacramento-San Joaquin Delta, California, USA [14]. In boreal regions, the annual warming effects of wetlands were seldom reported in previous case studies, probably due to the warming or cooling effect of wetlands being only observed in the growing season or summer. Seasonally, we find that ∆LST indicates a cooling effect of wetlands in summer but a warming effect in winter, which could be supported by several reported cases. In Northern Wisconsin, USA, using 2.4m air temperature observations of the eddy flux tower, Turner et al. indicated that the wetlands with more sedge meadow and forests (i.e., sheltered fens) had a stronger cooling effect than wetlands with less vegetation (i.e., open fens) due to the lower latent heat flux of the latter during the growing season [34]. The warming effect in winter was confirmed in a three-year (2015-2017) observation of eddy flux tower data which showed that the LST of wetlands was up to 0.73 • C higher than that of adjacent croplands in the Sacramento-San Joaquin Delta, California, USA [14]. Our satellite-based methods found the ∆LST ranging from -1.08 to 0.65 • C, indicating large spatial heterogeneity of ∆LST around the world. This warns that site-level observation of ∆LST may under-or overrepresent the large-scale impact of wetlands on local LST.

Mechanisms of the Warming or Cooling Effects of Wetlands
The spatial patterns of ∆LST (i.e., LST w f , LST w s , and LST w o ) can be interpreted by the surface energy budget as Equation (1) [15]. ∆S in the ecosystem is the chemical energy produced through photosynthesis and the temperature increase of the plant biomass [11]. The energy of ∆S is generally less than 5% of the net shortwave radiation [15]. In our study, the difference of ∆S between wetlands and other landcover types could be estimated by difference in gross primary production (∆GPP) as chemical energy fixed through photosynthesis. The ∆GPP between adjacent different biomass was less than 2 g C m −2 day −1 [35], equal to~80 kJ m −2 day −1 (0.9 W m −2 ). Thus, the difference in ∆S between wetlands and other landcover types was less than 1% of energy derived from ∆ET and ∆albedo, and was neglected in our study. By neglecting the ecosystem chemical energy variation ∆S' in Equation (2), Figure 6 shows the difference in surface energy budget between wetlands and other adjacent landcover types for summer and winter. Generally, ∆H was highly positively correlated with ∆LST [11], and ∆K and ∆LE were generally regarded as the two major energy fluxes to explain ∆H.
To be more intuitive, the contributions of ∆albedo and ∆ET to ∆LST in the entire seasonal cycle as well as in two seasons (i.e., summer and winter) are identified in Figure 7. In tropical regions, the contributions of the ∆albedo and ∆ET indicated a slight difference between summer (Figure 7d-f) and winter (Figure 7g-i) which explains the little seasonality of ∆LST; whereas in boreal regions, in winter, the contribution of ∆albedo was much higher than that of ∆ET (Figure 7e,f). In summer, the contributions ∆albedo and ∆ET to ∆LST were comparable (Figure 7c,d). Throughout the multiple years, 30% of temporal variation in ∆LST can be attributed to ∆albedo and ∆ET together, with their contributions being respectively 20% and 10% (Figure 7a-c). This suggests that in wetlands, other factors (e.g., G) may play a more important part in ∆LST. To better interpret the contributions of ∆albedo and ∆ET to ∆LST in different regions, Figure 8 compares the seasonal cycle of mean values of ∆LST and involved surface energy fluxes for typical regions I, II, III, and IV. In West Siberia (Region III), the seasonal cycle of ∆LST matched well with that of ∆SEF (i.e., the energy fluxes derived from ∆albedo and ∆ET) (Figure 8c,g,h), suggesting that the ∆LST between wetlands and other landcover types can be well interpreted by the corresponding ∆albedo and ∆ET in typical region III. In Canada, however, ∆LST started to increase in May and June and reached its peak in October and November, which was approximately 4~6 months later than the peak of ∆SEF (around June). A similar seasonal cycle of ∆LST was confirmed in three-year in situ observations in three different wetlands [14]. As ∆L can be calculated by multiplying the difference of the fourth power of wetlands' and other landcover types' LST by the Stefan-Boltzmann constant 5.67 × 10 −8 W m −2 K −4 , ∆L is up to 2.5 W/m 2 when the LST is 20 • C in the tropical regions and ∆LST is 0.7 • C. This is much smaller than ∆SEF so that ∆L is negligible. Here, we hypothesize that the time lag between the peaks of ∆LST and ∆SEF can be attributed to the seasonality of ground heat flux (G), as wetlands can store more energy in summer and then release it in winter due to larger water storage and the larger specific heat capacity of water (4200 J/Kg/ • C) compared to soils (~2500 J/Kg/ • C) [36]. Region I and II (northern and southern Canada) have a lot of wetlands including large lakes (e.g., the Great Lakes) that store plentiful water [37]. By in situ meteorological measurements, G of wetlands (14 W m −2 ) is two times higher than G of adjacent pasture during the growing season in the Třeboň Biosphere Reserve in southern Bohemia (48 • 49 -49 • 20 N, 14 • 39 -15 • 00 E) [38]. Global G has been modeled such that in the boreal regions, the highest G is found in summer at 19 W/m 2 , while the lowest G is in winter at −17 W/m 2 [39]. This confirms that the ground absorbs energy in summer and releases it in winter. However, the lack of large-scale ground flux G data hinders the interpretation of global patterns of ∆LST in different seasons, which needs further efforts. Besides, in tropical regions (IV), with their dense vegetation cover, the landcover maps have difficulty in distinguishing wetlands under forests (e.g., swamps) and tropical forests [40]. Consequently, ∆LST calculated in tropical regions (IV) based on satellite products in this study has a large degree of uncertainty, and ∆LST data in tropical regions (IV) should be interpreted with caution. Our study investigates the warming or cooling effects of wetlands and provides a cohesive picture of the difference of summer and winter surface heat budgets between wetlands and other landcover types ( Figure 6). Generally, wetlands had a larger absolute value of LE in summer than other landcover types while receiving less net shortwave radiation in winter due to more snow cover. The wetlands absorbed more energy (G) down to the ground than other landcover types in summer while releasing more energy in winter.

Conclusions
In this study, using satellite-based datasets, we show the global patterns of annual and seasonal ∆LST between wetlands and other landcover types, as well as ∆albedo and ∆ET. Our results show that in the tropical regions, the annual LST of wetlands are lower than those in forest, shrubland and open land; i.e., in these regions wetlands have a cooling effect while in the boreal regions, the wetlands have a warming effect. While ∆LST in the tropical regions has little seasonality, in the boreal regions, wetlands have a cooling effect in summer but a warming effect in winter. Only 30% of the spatial patterns of ∆LST and their seasonal variations can be attributed to the ∆albedo and ∆ET in west Siberia. The remaining part of ∆LST that cannot be interpreted by ∆albedo and ∆ET can be attributed to the ground heat flux, G, for instance in Canada, which can also interpret the time lag of the peaks of ∆LST and involved surface energy fluxes. Nowadays, there is an increasing interest in wetland restoration for mitigating climate change. Wetland restoration could greatly reduce nitrogen pollution [41], and methane and CO 2 emission [42], as well as change the local biophysical processes. This study highlights the potential climatic benefits of widespread wetland restoration in the future. As wetlands in different regions have different climatic effects (warming or cooling local temperature), there should be careful consideration of climatic benefits of wetland restoration in different regions.

Supplementary Materials:
The following are available online at https://www.mdpi.com/article/ 10.3390/rs13081439/s1, Table S1: Statistics of ∆LST, ∆albedo, and ∆ET in four regions; Figure S1: Spatial patterns of mean annual air temperature and LST of wetlands; Figure S2: Spatial patterns of mean annual albedo of wetlands; Figure S3: Differences in seasonal LST between wetlands, forests, and open lands; Figure S4: Differences in seasonal albedo between wetlands, forests, and open lands; Figure S5: Differences in seasonal ET between wetlands, forests, and open lands.
Author Contributions: Data analysis and interpretation was conducted primarily by Y.W.; writingoriginal draft preparation, Y.W.; writing-review and editing, Y.X., M.F. and S.P. This study was designed by S.P. All authors have read and agreed to the published version of the manuscript.