The Bowen ratio (β) is used to quantify heat transfer from the land surface into the air, which is becoming a hot topic in research on the biogeophysical effects of land use and cover changes. The Three-River Headwaters (TRH), as a sensitive and fragile region, was selected as the study area. The β for 2001–2018 was estimated from the evapotranspiration product (ETMOD16) of MODIS and the net radiation of the land surface through the albedo from GLASS. The ETMOD16 data were evaluated against the observation data (ETOBS) at two alpine grassland flux towers obtained from ChinaFLUX. The interannual trend of the β was analyzed by multiple linear regression (MLR) and structure model (SEM) with the multiple factors of precipitation, temperature, humidity, albedo, and normalized difference vegetation index (NDVI, MOD09Q1). The results show that the ETMOD16 values were significantly correlated with ETOBS, with a correlation coefficient above 0.70 (P < 0.01) for the two sites. In 2001–2018, the regional mean β was 2.52 ± 0.77 for the whole grassland, and its spatial distribution gradually increased from the eastern to western region. The interannual β showed a downward trend with a slope of –0.025 and a multiple regression coefficient (R2) of 0.21 (P = 0.056). Most of the variability (51%) in the interannual β can be explained by the linear regression of the above multiple factors, and the temperature plays a dominant role for the whole region. The SEM analysis further shows that an increasing NDVI results in a decreasing albedo with a path coefficient of –0.57, because the albedo was negatively correlated with NDVI (R2 = 0.52, P < 0.01), which indicates a negative and indirect effect on β from vegetation restoration. An obvious warming climate was found to prompt more evapotranspiration, and restoring vegetation makes the land surface receive more radiation, which both resulted in a decreasing trend in the annual β. This study revealed the biogeophysical mechanisms of vegetation restoration under a changing climate, and demonstrated the Bowen ratio can be applied as an indicator of climate-regulating functions in ecosystem assessments.
1 Introduction
Terrestrial ecosystems are both driven by and exert impacts on climate (Bonan, 2008, 2017). Land use, degradation, or restoration will change the physical properties of the ecosystem, such as albedo, leaf area index (LAI), and roughness length, which further alters the exchanges of energy and water fluxes between the ecosystem and the atmosphere, with consequent impacts on the local and regional climate. The effects of land-use changes have been widely studied (Betts, 2000; Cao and Li, 2000; Pielke et al., 2002; Feddema et al., 2005). To the best of our knowledge, however, the effects of land degradation and restoration on energy exchange have rarely been investigated in previous studies.
Grassland degradation and restoration have occurred on the Qinghai-Tibet Plateau, so it is an ideal location for determining the effects of land-use changes. Three-River Headwaters (TRH), located in the hinterland of the Qinghai-Tibet Plateau, is known as the water tower for the whole of China, and even for all of Southern Asia, because it includes the headwaters of the Yangtze, Yellow, and Lancang Rivers (Mei, 2000). In TRH, the Alpine grassland has been suffering from degradation since the 1970s, and the degraded area was 36.12% of the total grassland and continued to increase by a rate of 3.87% between 1990 and 2004 (Liu et al., 2008). Since 2005, many ecological measures have been implemented to restore the grassland and its degradation has been controlled to some extent (Shao et al., 2013; Shao et al., 2017; Huang et al., 2018). Meanwhile, the evapotranspiration of the grasslands increased significantly across the Qinghai-Tibet Plateau from 1982 to 2014 (Cui et al., 2019), which means that the latent heat has increased. Our previous analysis for the region showed a significantly warming and wetting climatic trend in the interpolated meteorological data in 1990–2012 (Yang et al., 2019) and an insignificant greening trend in the remote sensing-based normalized difference vegetation index (NDVI) since 2000 to 2015 (Wang et al., 2019). These findings raised the question of how the albedo, net radiation, and Bowen ratio will change with the greening vegetation and increasing latent heat under the background of the impacts of grassland degradation and restoration. The answer to this question will help us to better understand the processes and mechanisms of land degradation and restoration as two kinds of land-use change.
The energy exchanges were mainly determined through three methods, including on-ground observations, process-based modeling, and remote sensing-based estimations. The first method estimates the energy exchanges based on the Bowen ratio/energy balance (BREB), which requires the measurements of temperature and humidity at two levels, net radiation and soil heat flux (Dugas et al., 1991; Zhang et al., 2006; Zhang et al., 2017; Wang et al., 2018). For the second method, the process-based model describes the processes of the water cycle, energy exchange, carbon cycle, and numerous other biogeophysical and biogeochemical processes to capture the various land–atmosphere interactions (Bastiaanssen, 2000; Zhao and Li, 2015; Hao et al., 2016). The third method relies on remote sensing to retrieve some variables in BREB, such as land surface temperature and leaf area index (Liu et al., 2003; Geli et al., 2010). As remote sensing can provide continuous spatial–temporal information on vegetation and surface temperature, the remote sensing-based method has been widely applied. For example, the evapotranspiration products (MOD16) of the moderate resolution imaging spectroradiometer (MODIS) could be classified for the third method; these products have a 1 km spatial resolution and an 8-day time step from 2000 to the present, and are considered to have high accuracy (Mu et al., 2007; Mu et al., 2011).
Therefore, this study applied the newest evapotranspiration product (MOD16 C006) to analyze the spatial–temporal changes in the Bowen ratio and its impacting factors from climate and vegetation in the Three-River Headwaters region. We aimed to gain an understanding of the effects of ecosystem degradation and restoration on the local climate. With developments in ecosystem quality assessment (Wang et al., 2019), the second aim of this study is to provide a methodological reference for the monitoring and assessment of the climate-regulating function of the terrestrial ecosystem.
2 Data and methods
2.1 Overview of the study area
Three-River Headwaters (TRH) is located in Qinghai Province on the Qinghai-Tibet Plateau, with geographic co-ordinates from 31°39′ to 36°12′N and from 89°45′ to 102°23′E, and an altitude from 3335 to 6564 m (Fig. 1). The TRH region serves as the headwaters of the Yangtze River (YTR), Yellow River (YLR), and Lancang River (LCR), with the respective basin areas covering 46%, 28%, and 10% of the total TRH area (Shao and Fan, 2012). It has a typical highland continental climate, with an annual mean temperature from –18.96 to 7.18 °C and annual total precipitation from 208.58 to 845.26 mm according to the meteorological data from 2000 to 2018. The dominant vegetation of the study site is alpine grassland.
2.2 Data and processing
2.2.1 Evapotranspiration data
In this study, the land surface Bowen ratio was calculated from the land surface energy balance based on the remote sensing-based evapotranspiration (ET) data. The ET data were obtained from the MODIS product of MOD16A2 C006 with an 8-day temporal resolution and a 500 m spatial resolution for the period from 2001 to 2018 ( https://earthdata.nasa.gov/). The algorithm of the MOD16 product is based on the Penman–Monteith equation (Mu et al., 2011):
where λET is latent heat, s is the slope of the relationship between saturation vapor pressure and temperature (Pa °C–1); Rn is net radiation; ρ is air density and Cp is constant pressure specific heat (J kg–1 °C–1); esat is saturated water pressure and ea is actual water vapor pressure; ra is air resistance; γ is the humidity constant; and rs is canopy resistance. The algorithm was improved by considering vegetation transpiration, vegetation canopy wet surface evaporation, and soil evaporation (Mu et al., 2011):
where, λ is Latent heat transfer coefficient; ET is surface evapotranspiration; ETWET_C is evapotranspiration from the wetted canopy, ETtrans is evapotranspiration from the canopy, and ETsoil is evapotranspiration from the soil surface.
2.2.2 Net surface radiation data
The net surface radiation was estimated from the solar shortwave radiation and the albedo data was derived from remote sensing, which was developed in the GLOPEMCEVSA model (Wang, 2007):
where, Rn is the net surface radiation; Rns is the shortwave net radiation; Rnl is the long-wave net radiation; S ↓ is the downward short-wave radiation; r is the surface albedo; ε is the surface emissivity; L ↓ is the downward long-wave radiation; Ts is the surface temperature; σ is Stephen Boltzmann's constant; ea is actual water vapor pressure; n is actual sunshine hours; and N is reasonable sunshine hours. For detailed calculations, see Wang (2007).
2.2.3 Meteorological data
The meteorological data were obtained from the China Meteorological Science Data Sharing Network ( http://date.cma.cn/), including daily precipitation, maximum and min-minimum temperatures, relative humidity, wind speed, and hours of sunshine observed at meteorological stations in the study area and surrounding area. After a quality check and temporal interpolation of the missing data for a given time series in a given year (to provide a missing values rate of no more than 10%), the meteorological grid data (MeteoGrid) were interpolated through ANUSPLIN software (Hutchinson, 1991, 1998a, 1998b) with the digital elevation model data from SRTM of 90 m spatial resolution (Zhu et al., 2013) as a covariate. The interpolated MeteoGrid data have a spatial resolution of 250 m and a time step of eight days, to match with the MODIS data. The interpolated data were validated against the observed temperature and precipitation data from the seven eddy covariance towers. The data were significantly correlated with the observations and could explain 67% and 94% of the temporal variability in the observed precipitation and temperature data, respectively (Wang et al., 2017).
2.2.4 Observed latent heat data used for validation
The observed latent heat flux data from the two sites of ChinaFlux were used to validate the evapotranspiration product of MODIS (MOD16). The two sites were the alpine shrub at Haibei (101°20′E, 37°40′N, and 3190 m above sea level) and the alpine meadow at Dangxiong (91°03′E, 30°29′N, and 4300 m above sea level). The data were processed by a consistent method for all observations at the ChinaFLUX sites (Yu et al., 2006; Yu et al., 2008; Yu et al., 2014).
2.2.5 The data of vegetation index and surface albedo
The normalized difference vegetation index (NDVI) and land surface albedo data were both from MODIS for the period from 2001 to 2018. The NDVI data were calculated from the band reflectance data product of MODIS (MOD09Q1) of 250 m spatial resolution, and the spikes in the time series were removed through the Savitzky-Golay (SG) filtering algorithm in Timesat software (Jönsson and Eklundh, 2004). The albedo data were from the GLASS dataset with 1km spatial resolution that is based on a direct-estimation algorithm (Liu et al., 2013a; Liu et al., 2013b; Qu et al., 2014), and they were spatially resampled to a spatial resolution of 250 m.
2.2.6 Digital Elevation Model (DEM)
The DEM data were derived from the 90 m spatial resolution SRTM (Shuttle Radar Topography Mission) dataset produced jointly by the National Aeronautics and Space Administration (NASA) and the National Institute of Mapping (NIMA) of the Ministry of Defense, and they were resampled to a spatial interpolation of 250 m.
2.3 Calculation of the Bowen ratio and analysis methods
2.3.1 Calculation of the Bowen ratio
The Bowen ratio (Bowen, 1926) is defined as the ratio of sensible heat to latent heat (Equation 6). In this study, it is calculated from the net surface radiation and the evapotranspiration-based latent heat according to the surface energy balance (Equation 7).
where, β is Bowen ratio; Rn is the net radiation flux estimated from the remote sensing-based albedo; H is the sensible heat flux estimated as the difference between the Rn and the latent heat flux; and λET is estimated from the evapotranspiration of MOD16 because the soil heat flux G could be ignored due to the small variability of the intra-annual fluctuations.
2.3.2 Analysis method on impacts from multiple factors
Multiple linear regression analysis was applied to quantify the impacts from climatic and vegetation changes on the Bowen ratio through the following equation:
β =b0 +bP PAT +bTTAM +bH RHAM +bA AAM +bN NAM (8) where, β is the Bowen ratio; PAT, TAM, RHAM, AAM and NAM are the annual total precipitation, the annual mean temperature, annual mean relative humidity, annual mean Albedo and annual mean normalized difference vegetation index, respectively. The terms bP, bT, bH, bA and bN are the normalized regression coefficients which characterize the sensitivity of the Bowen ratio to the changes of the PAT, TAM, RHAM, AAM and NAM, respectively, and b0 is the intercept. The multiple correlation coefficient R2 of the regression equation and the significance level were applied to determine the performance of the fitted regression equation. The data were normalized in order to compare the relative contributions of the different independent variables.
2.3.3 Structural Equation Model (SEM)
Structural Equation Model (SEM) is a statistical method for analyzing the direct and indirect effects of independent variables on a dependent variable based on the correlation coefficient or covariance matrix between the variables; the main methods commonly used in SEM include path analysis, validation factor analysis, latent variable structural model, growth curve model, etc. (Hua and Cheng, 1999; Shipley, 2000; McDonald and Ho, 2002). Here, the path analysis to resolve the effects of TAM, PAT, RHAM, AAM, and NAM on the Bowen ratio was performed using the Lavaan package (Rosseel, 2012) in the R language, which does not require potential variables. Firstly, the correlation analysis is performed by the GGally package ggcorr function of the R language to determine the correlation between the independent variables, and then the following correlations (Equations 9 to 13) were considered to analyze the direct and indirect effects of the independent variables on the Bowen ratio by the path analysis function in the lavaan package.
3 Results and analysis
3.1 Changes of climate and vegetation
The interannual changes from 2000 to 2018 of the annual total precipitation (PAT), the annual mean temperature (TAM) and NDVI (NAM) for the whole region and the three sub-basins in Three-River Headwaters were shown in Fig. 2. According to the multi-year averaged values, the annual total precipitation was 484.19 ± 130.43 mm, the annual mean temperature was –3.33 ± 2.39 °C, the annual mean humidity was 56.27% ± 4.66%, the annual mean NDVI was 0.19 ± 0.11, and the annual mean albedo was 0.22 ± 0.04 for the whole region.
During the period 2000–2018, the climatic changes showed insignificant wetting but significant warming trends, with increasing rates of 29.7 mm per decade (P = 0.20) for the PAT, and 1.0 °C per decade (P < 0.01) for the TAM, respectively. The vegetation was greening by 0.6% per decade (P = 0.05), while the AAM was decreasing by 1% per decade (P = 0.04).
The climatic changes and vegetation responses varied among the three basins of the Yellow River, Yangtze River and Lancang River. The trends of precipitation were insignificant, with slopes of 14.02 mm, 45.05 mm and 0.37 mm per decade and significance levels (P) of 0.09, 0.41, and 0.62, respectively, for the three basins. However, the temperatures showed significant warming trends, with slopes of 1.02 °C, 1.14 °C and 0.98 °C per decade and the P values were less than 0.001 for each of the three basins. These results indicated that the whole region experienced a warming climate, and the Yellow River Basin had a warming and wetting climate.
Fig. 2
Interannual changes of the annual total precipitation (PAT), the annual mean temperature (TAM) and NDVI (NAM) of the whole region (TRH) and the three sub-basins of Yellow River (YLR), Yangtze River (YTR) and Lancang River (LCR) in the Three-River Headwaters, Qinghai, China, from 2000 to 2018.
![img-z5-1_305.jpg](ContentImages/Journals/jore/12/3/j.issn.1674-764x.2021.03.001/graphic/WebImages/img-z5-1_305.jpg)
The trends of the NDVI were also insignificant, but they were increasing at rates of 0.5%, 0.8% and 0.7% per decade for the three basins. In contrast, the albedo showed a decreasing trend in each basin, with the rates of –1.3%, –1.2% and –1.1% per decade, which meant that more solar radiation would be absorbed by the land surface according to Equation 4 along with the greening of vegetation in that region.
3.2 Evapotranspiration validation
The data were sampled for the 3×3 pixel areas centered around the eddy covariance towers in the MOD16A2 products (ETMOD16), and the mean of the data was compared with the observations (ETOBS) for the eddy covariance towers at Haibei and Dangxiong grasslands. The results showed that the ETMOD was significantly correlated with the ETOBS, and the multiple correlation coefficients were 0.61 and 0.45 for the two sites respectively, each with a P value of less than 0.01; which indicated the MOD16 product could explain 61% and 45% of the seasonal changes in the observations and could capture the seasonal changes in the ET of the alpine grass reasonably well, although some underestimations were found as shown by the RMSE of 6.44 and 8.19 mm m–2 (8 day)–1 at Haibei and Dangxiong, repectively (Fig. 3).
3.3 Bowen ratio and its spatio-temporal changes
3.3.1 Energy exchanges and Bowen ratio
Figure 4 shows the Bowen ratio, latent heat, net radiation, and sensible heat in the Three-River Headwaters region for the period from 2001 to 2018. The net surface radiation was 3304.73 ± 342.79 MJ m–2 yr–1 with the latent heat flux from the MOD16A2 C006 of 980.67 ± 150.71 MJ m–2 yr–1. According to surface energy balance (Equation 7), if the soil heat flux was ignored, the sensible heat flux was 2327.18 ± 334.09 MJ m–2 yr–1 and the resulting Bowen ratio would be 2.52 ± 0.77, which means that the sensible heat was 2.52 times the latent heat in the Three-River Headwaters in the study period from 2001 to 2018.
The spatial distribution of the Bowen ratio shows it is gradually increasing from the eastern to western region (Fig. 4a). The highest Bowen ratio was found in western Tanggulashan, with an average value of 3.13 and a coefficient of variation of 20.37%; while the southeastern Gander has the smallest Bowen, with an average value of 2.00 and a smaller spatial variation with a coefficient of variation of 21.8%.
Fig. 3
The evapotranspiration (ETMOD16) products of MODIS (MOD16A2) were validated against the observations (ETOBS) through the linear regression analysis (a, c) and the seasonal change comparisons (b, d) on the alpine grasslands at Haibei Station (a, b) and Dangxiong Station (c, d) on the Qinghai-Tibet Plateau, China.
![img-z6-1_305.jpg](ContentImages/Journals/jore/12/3/j.issn.1674-764x.2021.03.001/graphic/WebImages/img-z6-1_305.jpg)
Fig. 4
Distribution patterns of multi-year averages of Bowen ratio (a), Latent heat (b), Net Radiation (c), and Sensible Heat (d) of the Three-River Headwaters from 2001 to 2018.
![img-z6-3_305.jpg](ContentImages/Journals/jore/12/3/j.issn.1674-764x.2021.03.001/graphic/WebImages/img-z6-3_305.jpg)
Among the three basins, the Yellow River has the lowest Bowen ratio, with net radiation of 3271.30±302.37 MJ m–2 yr–1, latent heat based on MOD16 products of 1051.41 ± 126.15 MJ m–2 yr–1, and sensible heat according to the surface energy balance of 2218.06 ± 279.52 MJ m–2 yr–1. Therefore, its Bowen ratio was 2.21 ± 0.35, that is, about 67.80% of the net surface radiation was used to heat the air, while only 32.14% of that was used for vegetation transpiration and soil evaporation.
The net radiation from the Yangtze River source was 3375.44 ± 290.51 MJ m–2 yr–1; the latent heat was 945.92 ± 172.92 MJ m–2 yr–1, accounting for 28.28% of the net radiation; the sensible heat was 2424.32 ± 284.51 MJ m–2 yr–1, accounting for 71.82% of the net radiation; and the Bowen ratio was 2.68 ± 0.30. The net radiation from the Lancang River source was 3168.90 ± 426.34 MJ m–2 yr–1, the latent heat was 945.92 ± 172.92 MJ m–2 yr–1 (29.85% of the net radiation), the sensible heat was 2227.56 ± 365.25 MJ m–2 yr–1 (70.29% of the net radiation), and the Bowen ratio was 2.55 ± 0.35. The above results of the three basins showed that nearly one-third of the energy is used for the evapotranspiration and that more than two-thirds is used for sensible heat in the whole area, and Yangtze River has the highest Bowen ratio, followed by Lancang River, while the Yellow River has the lowest ratio.
3.3.2 Inter-annual trends of Bowen ratio
According to the linear regression analysis, the Bowen ratio of the whole TRH was significantly decreasing at a rate of –2.5% yr–1 (R2 = 0.21, P = 0.056) in the period from 2001 to 2018. This trend meant that the amount of net surface radiation allocated as sensible heat was decreasing while the amount allocated as latent heat was increasing during the study period, with a confidence level of 94.4% (Fig. 5a). The Bowen ratio of the three basins showed decreasing trends, at a rate of 0.03 per year for the Yellow River basin and 0.02 per year for both the Yangtze and Lancang River basins (Fig. 5b).
The spatial pattern of the Bowen ratio trend showed an increasing trend in 13.57% and a decreasing trend in 72.75% of the area of the whole region (Fig. 6). The increasing trend occurred over the southwestern and central regions, especially in the southeastern Tanggulashan and the western Zaduo, while the decreasing trend occurred over the northeastern region, including Xinghai, Tongde, Zeku and Henan.
3.4 Analysis of the impact factors on the Bowen ratio
3.4.1 Multiple linear regression analysis
Multivariate linear regression analysis was applied to explore the impacts on the Bowen ratio from the climatic factors (PAT, TAM, RHAM), vegetation (NAM) and albedo (AAM), and the results based on the regional mean values are shown in Table 1. With the confidence level of 84%, the five factors can explain 51% of the interannual variation of the Bowen ratio for the whole region, and 57%, 48%, and 46% for the Yellow River Basin, Yangtze River Basin, and Lancang River Basin, respectively, according to the correlation coefficients (R2). Additionally, the TAM, RHAM, NAM each has a negative effect on the Bowen ratio, while the effects of PAT and AAM varied among the basins according to the standardized regression coefficients. The temperature (TAM) had the greatest effect on the Bowen ratio at the whole regional scale, with a standardized regression coefficient of –0.43, followed by NAM (–0.27), RHAM (–0.14), PAT (–0.08) and AAM (–0.05).
Table 1
Multiple linear normalized regression coefficients for the three sub-basins in the Three-River Headwaters and the whole region
![img-z8-8_305.gif](ContentImages/Journals/jore/12/3/j.issn.1674-764x.2021.03.001/graphic/WebImages/img-z8-8_305.gif)
Fig. 7
The trends of Bowen ratio attributed to climate (PAT, TAM and RHAM), vegetation greenness (NAM) and AAM through the multiple linear regression with its multiple correlation coefficient (a), significance level (b), and the regression coefficients multiplied by the trends of the dependent variables (c,d).
![img-z8-9_305.jpg](ContentImages/Journals/jore/12/3/j.issn.1674-764x.2021.03.001/graphic/WebImages/img-z8-9_305.jpg)
The multiple linear regression was performed at the pixel scale to explore the spatial patterns of the effects of the five factors on the Bowen ratio (Fig. 7). Over an area of 86% of the whole region, these five factors can explain more than 50% of the variability in the Bowen ratio (Table 2), with the exception of the upstream and valleys of the rivers where the regression with the five factors was insignificant and other unknown factors probably dominate in driving the changes of the Bowen ratio (Fig. 7a, b, Fig. 8).
Among the five factors, the changes in annual mean temperature contributed more than 90% of the variability in the Bowen ratio with the specific contribution rates of 92%, 95% and 118% in the Yellow River, Yangtze River and Lancang River, respectively, and 90% for the whole region. This was followed by the NDVI, with contribution rates of 10.2%, 9.7% and 11.7% for the three basins, respectively. However, the contribution of the humidity was contrasted between the two basins of Lancang River (–31.6%) and Yellow River (–0.16%). The contribution of precipitation was less than 10% and varied among the three basins from –1.47% (Lancang), to 0.41% (Yellow) and 8.47% (Yangtze).
Fig. 8
The standard regression coefficients (a) and their contributions (b) on the Bowen ratio from the annual mean temperature (TAM), annual total precipitation (PAT), annual mean humidity (RHAM), annual mean albedo (AAM) and annual mean normalized difference vegetation index (NAM) in the Three-River Headwaters in the period from 2001 to 2018.
![img-z9-1_305.jpg](ContentImages/Journals/jore/12/3/j.issn.1674-764x.2021.03.001/graphic/WebImages/img-z9-1_305.jpg)
The effects of the three types of factors on the spatial pattern varied (Fig. 7b, c). Over the southwestern region, precipitation had the largest influence, but the greatest influence was vegetation status (expressed by NDVI) over the central region, and temperature over the northeastern region.
The Bowen ratio trend was dominated by the increasing air temperature, with a contribution above 90%, which would be attributed to its sensitivity to air temperature and the warming trend in TAM. The Bowen ratio was most sensitive to the interannual variability in the TAM time series (–43%), followed by NAM (–27%) and RHAM (–14%), but was less sensitive to the changes in PAT and AAM (Table 1). On the other hand, the region experienced a significant warming change at the rate of 1.02 °C per decade (P < 0.01) with an insignificant change in PAT, and the vegetation was greening at the rate of 0.6% per decade (P = 0.05) along with its decreasing AAM at the rate of 1% per decade (P = 0.04). Therefore, the Bowen ratio showed a decreasing trend induced by temperature (–31.7% per decade), NDVI (–4.5% per decade) and albedo (1.1% per decade) (Table 2). 3.4.2 Analysis by the Structural Equation Model (SEM) The five key variables were correlated with each other (Fig. 10a). First of all, temperature was significantly and positively correlated with PAT (R = 0.55, P < 0.001). The NAM was correlated with TAM (0.71), PAT (0.43) and RHAM (0.42). Those correlations among the variables influenced the interannual trend in the Bowen ratio (β). Therefore, the SEM was applied to determine the direct and indirect effects from the above five variables on the Bowen ratio and the results are shown in Fig. 10b and Table 3.
Table 2
The 2001–2018 interannual trends of Bowen ratio and related climate variables in the Three-River Headwaters
![img-z9-9_305.gif](ContentImages/Journals/jore/12/3/j.issn.1674-764x.2021.03.001/graphic/WebImages/img-z9-9_305.gif)
The Bowen ratio was negatively influenced by TAM (–0.53) and RHAM (–0.30), while positively influenced by AAM (0.43). An effect of albedo was not revealed by the MLR, but more information was given by SEM. The AAM was affected by the NAM (–0.57) and RHAM (0.53), and it then played a positive and direct role in influencing the Bowen ratio. The effect of TAM was too strong to be revealed by both methods. PAT directly affected RHAM (0.59) and TAM (0.55), and indirectly affected the Bowen ratio. Furthermore, TAM strongly affected vegetation changes with a path coefficient of 0.72 with NAM. In addition, vegetation changed its AAM, that is, the increasing NAM would result in a decreasing AAM according to the path coefficient of –0.57, and NAM was significantly correlated with AAM (R2 = 0.52, P < 0.001), according to the linear regression (Fig. 11). These results meant that the warming climate drove vegetation restoration and both affected the interannual variability in the Bowen ratio by altering the albedo of the insignificant vegetation greening.
Fig. 9
Effects of different factors on the interannual variation of the Bowen ratio in the Three-River Headwaters: (a) PAT, (b) TAM, (c) RHAM, (d) AAM, and (e) NAM.
![img-z10-1_305.jpg](ContentImages/Journals/jore/12/3/j.issn.1674-764x.2021.03.001/graphic/WebImages/img-z10-1_305.jpg)
Table 3
Fitness test results of SEM
![img-z10-4_305.gif](ContentImages/Journals/jore/12/3/j.issn.1674-764x.2021.03.001/graphic/WebImages/img-z10-4_305.gif)
Fig. 10.
The correlation matrix (a) and the structural equation model (b) for the variables of annual total precipitation (PAT), the annual mean temperature (TAM), annual mean relative humidity (RHAM), and the vegetation factors of the annual mean Albedo (AAM) and annual mean normalized difference vegetation index (NAM), and the Bowen ratio (β). Note: In subplot (b), the width of a line indicates the correlation and the green and red colors indicate a positive and a negative effect, respectively. The labeled values are the normalized path coefficients.
![img-z11-1_305.jpg](ContentImages/Journals/jore/12/3/j.issn.1674-764x.2021.03.001/graphic/WebImages/img-z11-1_305.jpg)
4 Discussion
4.1 Uncertainty in the Bowen ratio calculation
In this study, the Bowen ratio was calculated from the remote sensing-based evapotranspiration (ET) data from the MODIS product (MOD16 V006). Therefore, the uncertainty mainly came from the MOD16 data. Here the product (ETMODIS) was validated against the ET (ETOBS) which was estimated from the latent heat data observed on the eddy covariance towers in the alpine grassland. As the validation showed, the MOD16 data were significantly well-correlated with the observations. However, the MOD16 could only explain 45% and 61% of the variability in the observations in the two alpine grassland sites of Dangxiong and Haibei on the Qinghai-Tibet Plateau, respectively. The uncertainties, however, probably came from the observations since the eddy covariance technique is considered to be the standard method for determining energy and substance fluxes (Tian et al., 2006). Despite advances in the development and improvement of instrumentation, energy balance closure is a problem which persists (Li et al., 2004). On the other hand, there are some uncertainties in the MOD16 products which arise from its algorithm and the input data; so the algorithm should be further improved, especially the soil water balance estimation which was not considered in the most recent ET algorithm (Yan et al., 2012). In the future, more flux sites should be applied to improve the ET algorithm by considering surface energy balance through the optimization of the parameters.
4.2 Factors affecting the inter-annual variation of the Bowen ratio
As the ratio of sensible heat (H) to latent heat (λET) in Equation 6, the Bowen ratio is affected by climate, vegetation and soil properties (Monteith and Unsworth, 2013; Tang et al., 2014). Previous studies have reported that the Bowen ratio was affected by soil moisture through the water supply for surface evaporation (Wang et al., 2003; Zhao et al., 2008), by surface roughness through changes in the aerodynamic impedance (Pitman, 2003), and by the freeze-thaw cycle (Ge et al., 2016). However, in this study, temperature was found to be the dominant factor along with humidity and albedo as the secondary factors, which was also found in a sub-tropical forest (Tang et al., 2014). Temperature plays a dominant role because it is one of important factors controlling the many physiological and chemical reaction processes in ecosystems and it acts as a regulator and allocator of energy, water and nutrients in ecosystems (Yu and Xu, 2009; Yao and Zhang, 2015). Humidity has a strong influence on vegetation growth and production, as well as biological diversity in arid and semi-arid regions (Ye et al., 2012). Vegetation factors, here represented by NDVI, seem have no effect on the Bowen ratio, and in fact, through altering albedo to affect the Bowen ratio, such interactions were well quantified by the SEM analysis (Fig. 9b).
4.3 Importance and aspects for further study
Land use and cover change, such as planting, grazing, deforestation and afforestation, was considered as one of the important drivers of climate by altering surface energy balance (Turner et al., 1995; Li, 1996; Pitman, 2003). Alpine grasslands have experienced degradation since the early 1970s and restoration since 2005 (Liu et al., 2008; Shao et al., 2013; Shao et al., 2017; Huang et al., 2018), which has probably resulted in the changes of surface energy balance, known as the biogeophysical effect (Zhao et al., 2019). On the interannual scale, the Bowen ratio showed a decreasing trend for the warming climate changes in the background of global warming, and a decreasing albedo due to vegetation restoration as shown by the increasing NDVI (Fig. 2). A decreasing Bowen ratio, that is, an increasing latent heat or a decreasing sensible heat, will inhibit the warming effect of the surface on the atmosphere. In this sense, the Bowen ratio could be applied to assess ecosystem quality as an indicator of the regulating function of the ecosystem.
5 Conclusions
In this paper, the Bowen ratio was calculated from the evapotranspiration products of MODIS (MOD16) and its interannual trends and affecting factors were analyzed in the Three-River Headwaters region of China over the period of 2001–2018. The Bowen ratio was found to have a trend of decreasing by 25% per decade, which was mainly attributed to a warming climate, changing humidity and greening vegetation. Specifically, a changing climate dominated by global warming, superimposed on slightly increasing precipitation, will promote vegetation growth and a decreasing surface albedo, leading to an increasing net surface radiation; but at the same time, climate warming dominates and directly contributes to increased surface evapotranspiration, leading to a decrease in the Bowen ratio. However, the validation of the data from MOD16 indicated that its algorithm has some room for improvement. It is essential to develop a new algorithm and improve the accuracy of latent and sensible heat estimations by the fusion of the ecosystem process model and observations from more eddy covariance towers. This study expands our understanding of the biogeophyscial effects of grassland restoration as one type of land use and cover change, and provides a knowledge foundation for using the Bowen ratio as an indicator of the climate regulating function of ecosystem assessments in future.