Next Article in Journal
Seasonal, Weekly, and Diurnal Black Carbon in Moscow Megacity Background under Impact of Urban and Regional Sources
Next Article in Special Issue
An Estimation of Precipitation Retention Time Using Depth Metres in the Northern Basin of Lake Biwa
Previous Article in Journal
Improving Combustion Technology for Cooking Activities for Pollutant Emission Reduction and Carbon Neutrality
Previous Article in Special Issue
Spatiotemporal Changes of sc-PDSI and Its Dynamic Drivers in Yellow River Basin
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evapotranspiration Variations of the Minjiang River Basin in Southeastern China from 2000 to 2019

1
National Demonstration Center for Experimental Geography Education, Fujian Normal University, Fuzhou 350007, China
2
Institute of Geography, Fujian Normal University, Fuzhou 350007, China
3
Fujian Provincial Engineering Research Center for Monitoring and Assessing Terrestrial Disasters, Fujian Normal University, Fuzhou 350007, China
4
School of Urban and Environmental Science, Huaiyin Normal University, Huai’an 223300, China
*
Authors to whom correspondence should be addressed.
Atmosphere 2022, 13(4), 562; https://doi.org/10.3390/atmos13040562
Submission received: 30 December 2021 / Revised: 25 March 2022 / Accepted: 28 March 2022 / Published: 30 March 2022
(This article belongs to the Special Issue The Water Cycle and Climate Change)

Abstract

:
Evapotranspiration is one of the critical processes in the terrestrial hydrological cycle, and the assessment of evapotranspiration is essential for understanding the regional hydrological cycle. In this study, the Minjiang River Basin, a typical watershed in the humid subtropical climate zone, is selected as the study region. The Penman-Monteith equation and the dual crop coefficient method are used to calculate the actual evapotranspiration (ETa) at seven meteorological stations within the study basin. Meanwhile, the applicability of the Global Land Data Assimilation System-Noah (GLDAS-Noah) ETa data in the Minjiang River Basin is evaluated based on stations P-M equation results, then to analyze the changes of the ETa in the Minjiang River Basin from 2000 to 2019. The results show that the GLDAS-Noah ETa data are well applicable in the Minjiang River Basin (R2 > 0.9 and NSE > 0.8). The ETa in the basin shows an increasing trend since 2000, and the increasing rate is 3.60 mm·yr−1 (p < 0.01). The seasonal variation results show that ETa tends to increase in winter and spring, with increasing rates of 1.10 mm·yr−1 (p < 0.01) and 2.60 mm·yr−1 (p < 0.01), respectively, while the ETa did not change significantly in summer and autumn. Annual air temperature has the largest effect on annual ETa (59.6%), followed by precipitation at 33.9%. ETa increased in spring was mainly influenced by increasing temperatures (89.4%) in the Minjiang River Basin from 2000 to 2019. The research results are of great benefit to further improve the understanding of ETa variations in the basin under global warming.

1. Introduction

Evapotranspiration bridges land-atmosphere interactions and plays an essential role in regional hydrological cycle and water balance [1,2]. Changes in actual evapotranspiration (ETa) are closely related to meteorological conditions (such as air temperature, precipitation, and wind speed) and surface characteristics (including soil moisture, crop type and growth status). These factors significantly change under climate change and human activities, thereby affecting ETa [3,4]. ETa is an essential factor affecting runoff, and its accurate estimation at the basin scale is crucial for further understanding hydrological cycle processes and water resource management, and even for promoting the sustainable development of water resources [5].
ETa estimation methods mainly include traditional and remote sensing methods. Traditional methods for estimating ETa include the water balance method [3], the Bowen ratio-energy balance method [6], Penman-Monteith (P-M) equation [7], and the method proposed by Dalton (1802) reflecting the relationship between the surface evaporation rate and impact factors (atmospheric physical properties and physiological characteristics of vegetation) [8]. These methods have several advantages, i.e., simple structure, high accuracy, and relatively good applicability. However, they only perform well at the single-point and station scales [9], while the errors are relatively large at the regional scale. In the last 30 years, with the development of remote sensing technology with a high spatio-temporal resolution, it has been possible to quantitatively estimate the regional ETa [10], with the advantages of all-weather, real-time and better accuracy. However, this technology is highly influenced by clouds, thus requiring validation by ground-based stations. In addition, land surface modeling and data assimilation, such as, the Global Land Data Assimilation System (GLDAS) model and reanalysis data (such as the ERA5 dataset), can also provide ETa data. However, these estimates needs to be validated based on observations before application. Therefore, combining traditional methods with remote sensing methods, land surface process models and reanalysis data can improve the estimation accuracy of ETa.
Under the background of global warming, ETa changes show significant spatio-temporal differences in China [11]. ETa in the northeastern part of the Songhua River Basin [12], the North China Plain [13] and the coastal region of southeastern China [14] have shown increasing trends. In the western, northern, and southeastern parts of the Loess Plateau, ETa tended to decrease [15]. Tibetan Plateau averaged terrestrial ETa increased significantly (1.87 mm yr−1, p < 0.001) from 1982 to 2016 and is due primarily to precipitation increased [16]. The ETa in Southwest China had tended to decrease from 1960 to 2013, mainly influenced by the shorter than solar radiation and lower wind speed [17]. The ETa in arid Northwest China showed a decreasing trend from 1958 to 1993 and an increasing trend from 1994 to 2010, mainly closely related to the near-surface wind speed variation [18]. Liu et al. (2011) [19] revealed that the potential evapotranspiration (PET) in the basins of Southeast China showed a decreasing trend between 1960 and2007, which was most sensitive to changes in the maximum daily temperature. The regional differences in ETa variability reflect regional hydrothermal conditions in the warming context and regional surface characteristics. Therefore, exploring ETa variations can contribute to understanding the characteristics of regional hydrological changes in the warming context.
As the largest river of Fujian Province, the Minjiang River Basin (MRB) is of great importance to the economic, social, and ecological development of Fujian. In the past 20 years, notable warming in the MRB [20], and the frequency and intensity of successive autumn-winter-spring meteorological droughts in the basin have tended to increase [21]. The increase in the temperature and ETa is likely one of the critical reasons for droughts in the MRB. However, current research results lack attention to ETa variations in the MRB. There are only eight meteorological stations in the MRB, which makes it difficult to assess ETa variations at the basin scale. Due to the complex and diverse topographic conditions in the basin, its middle and upper reaches usually experience cloudy and foggy weather, which seriously affects the accuracy of remote sensing products.
In this study, we first analyze whether GLADS-Noah ETa can produce results that are compatible with station P-M results in the MRB. Meanwhile, we use the GLDAS-Noah ETa data to describe the spatio-temporal patterns of ETa. Then, we explore the influencing factors of ETa spatio-temporal variability in the MRB.

2. Study Area, Dataset and Methods

2.1. Study Area

The MRB is located at 25°23′ N–28°19′ N, 116°23′ E–119°43′ E (Figure 1), with an area of 6.1 × 104 km2, a total length of 562 km and an average annual runoff of 940 mm. The MRB belongs to the subtropical monsoon climate zone, with an average annual rainfall of 1600–1800 mm and annual temperature 16–20 °C [20]. The MRB has a forestland area of 4.15 × 104 km2, and a cultivated land area of 0.92 × 104 km2.

2.2. Data Sources

2.2.1. Meteorological Data

Meteorological data used in this study come from the National Meteorological Information Center. In the study area, the daily data from eight meteorological stations (Table 1) during January 2000–December 2019 were collected, including the average, maximum and minimum temperature, precipitation, relative humidity, and hours of sunshine. The ETa of the eight meteorological stations can be calculated based on these daily observation data.

2.2.2. GLDAS-Noah Data

The GLDAS version 2.1 is forced by a combination of satellite-based and ground-based observation products, and this system uses advanced land surface modeling and data assimilation techniques to generate optimal fields of land surface states and fluxes [22]. The land surface model parameters include four categories: elevation, soils, vegetation type and extent, and leaf area index [22]. Therefore, the study did not account for changes in land cover. Currently, the GLDAS data are used to drive four land surface models including Noah, catchment and community land surface models and the variable infiltration capacity (VIC) land surface model. In this research, we select the monthly ETa and precipitation data with a spatial resolution of 0.25° × 0.25° from the GLDAS-Noah model during 2000–2019. Figure 2 shows that the basin average of monthly GLDAS-Noah precipitation is quite close to average of eight stations observations, and the determination coefficient (R2) is 0.9, indicating that the monthly GLDAS-Noah precipitation products are effective in simulating the precipitation in the MRB.

2.3. Method

2.3.1. Station ETa Calculation

In this study, ETa is calculated in two steps. Firstly, the PET of station is estimated by the P-M equation. Then, ETa at station is calculated by combining the station PET with the dual crop coefficient.
ET a = P E T × K c
where ETa is actual evapotranspiration, PET is potential evapotranspiration, and Kc is crop coefficient.
The P-M equation (Equation (2)) used to calculate PET in this study. Previous study indicated that P-M calculate PET are quite consistent with the measured values [23].
E T 0   = 0.408 Δ ( R n G ) + r 900 T + 273 U 2 ( e s e a ) Δ + r ( 1 + 0.34 U 2 )
where ET0 indicates the PET (mm·day−1), Rn the net radiation (MJ·m−2·day−1), G the soil heat flux density (MJ·m−2·day−1), T the 2-m air temperature (°C), U2 the 2-m wind speed (m·s−1), es the saturated water vapor pressure (kPa), ea the actual vapor pressure (kPa), esea the saturated vapor pressure deficit (kPa), ∆ the slope of vapor pressure curve (kPa·°C−1), γ the psychrometric constant (kPa·°C−1).
The method to calculated Kc was proposed by Wright in 1982 [24]. The previous study shows that the NDVI-based Kc values agree with the measured Kc values [25]. The Kc was calculated from NDVI and applied in this study.
K c = K s × K c b + K e
K c b = 1.07 × [ 1 ( NDVI max NDVI NDVI max NDVI min ) 0.84 0.54 ]
K e = β × ( 1 f c )  
f c = 1.3514 × NDVI 0.2811
where Kcb the basic crop coefficient reflecting the crop evapotranspiration, and Ke the evapotranspiration coefficient reflecting the soil surface evaporation. Ks represents the water stress coefficient reflecting the effect on the crop evapotranspiration when the soil water content in the root zone is insufficient (Ks = 1 means water is not a limiting factor) [26]. In this study, the value of Ks is taken as 1 because the MRB is a subtropical monsoon zone, so there is sufficient water. NDVI indicates the vegetation index, which is from NOAA Global Inventory Monitoring and Modeling System, where NDVImax and NDVImin are the maximum and minimum values of monthly NDVI from 2000 to 2014, respectively. fc denotes the effective area ratio of vegetation covering the soil surface, β an empirical coefficient depending on the Kc values at the early and middle stages of crop growth. Based on previous research results [26], β is taken as 0.25 in this study.

2.3.2. Evaluation Metrics

In this study, we use a new comprehensive index DISO (distance of indices between simulation and observation) [27], R2 and Nash-Sutcliffe efficiency coefficient (NSE) to evaluate the applicability of the GLDAS-Noah ETa data in the MRB. The DISO can characterize the performance of the simulated data evaluated based on the distance between the simulated and the observed data [27]. The Smaller the DISO value, the higher the accuracy of the GLDAS-Noah ETa data. The expression is as follows (Equation (6)).
DISO = NRMSE 2 + RB 2 + ( R 1 ) 2
where NRMSE, RB and R are the normalized root mean square error, relative bias and Pearson correlation coefficient, respectively, which are calculated as follows.
NRMSE = 1 OB ¯ 1 n i = 1 n ( SI i OB i ) 2  
RB = 1 n i = 1 n ( SI i OB i ) 1 n i = 1 n OB i = i = 1 n ( SI i OB i ) i = 1 n OB i  
R = i = 1 n ( SI i SI ¯ ) ( OB i OB ¯ ) i = 1 n ( SI i SI ¯ ) 2 i = 1 n ( OB i OB ¯ ) 2    
where OB denotes the results from the P-M method at the station, SI the GLDAS-Noah data, n the time series, OB ¯ the mean of OB, and SI ¯ the mean of SI.
The NSE (Equation (10)) is usually used to assess the fit of the simulated values to the observed values, and NSE values range from 0 to 1.
N S E = 1 1 n ( OB i SI i ) 2 1 n ( OB i OB ¯ ) 2
Therefore, if the DISO is close to 0, the R2 and NSE values are close to 1, this is indication that the simulation performance of the GLDAS-Noah model is relatively better.

2.3.3. Partial Least-Squares (PLS) Regression Model

In this study, we used the PLS regression model to estimate the impacts of climate variables (precipitation, temperature, and wind speed) on the ETa, and determine the contribution of each climate variables.
ET a = a + b × P + c × T + d × W
where a, b, c, and d are regression coefficients; P, T, and W are standardization of precipitation, temperature, and wind speed, respectively. Once the regression coefficients of this equation were defined, we used the function to analysis the relative contribution of each climate variables on ETa [28,29].
{ ρ p = | b | / ( | b | + | c | + | d | ) ρ T = | c | / ( | b | + | c | + | d | ) ρ W = | d | / ( | b | + | c | + | d | )
where ρ P , ρ T , and ρ W are relative contribution of precipitation, temperature, and wind speed on the ETa.

3. Results

3.1. Applicability Assessment of the GLDAS-Noah ETa Data

For the applicability assessment of ETa data, the calculated results from the P-M equation at seven stations are used as observations, and the GLDAS-Noah ETa data in the corresponding grid points at the stations are taken as simulations. The results (Table 2) show that the DISO for each station is less than 0.3, the NSE is greater than or equal to 0.8, and the R2 values are all greater than 0.9. Thus, these results indicated that there is highly consistence between stations P-M result and GLDAS-Noah ETa data in the corresponding grid points of meteorological stations.
Figure 3a shows that the GLDAS-Noah ETa agrees with the station P-M ETa. Although the GLDAS-Noah ETa underestimates the actual evapotranspiration, the Pbias is less than 15%. Meanwhile, the decidability coefficient of monthly GLDAS-Noah ETa and station P-M results reaches 0.97 (Figure 3b). Therefore, the GLDAS-Noah ETa data have good applicability in the MRB and can assess the ETa variation characteristics in this region.

3.2. Spatio-Temporal Variations of the ETa

During 2000–2019, ETa in the MRB shows an increasing trend, and the overall increasing rate is 3.60 mm·yr−1, passing the significance test of 0.01 confidence level (i.e., p < 0.01 as shown in Table 3). This result is consistent with the findings of Huang (2020) [30]. Moreover, there is a significant seasonal variation in the ETa in the MRB. Winter includes December, January and February (DJF), spring for March, April, and May (MAM), summer for June, July and August (JJA), and autumn for September, October and November (SON). As expected, the largest ETa occurs in summer (multi-year average of 342.1 mm), followed by spring and autumn (multi-year average of approximately 200 mm), and the smallest in winter (multi-year average of 92.6 mm). Most of the increasing trends in ETa is occurring during the Spring months (2.60 mm yr−1, p < 0.01). Interesting that a significant amount is happening during winter (1.10 mm yr−1, p < 0.01). However, the ETa variations in summer and autumn are nonsignificant.
The multi-year average ETa in the MRB ranges from 700 mm to 1100 mm, regional average is 846 mm (Figure 4a). The average annual ETa of GLDAS-Noah is reasonable, because of average annual precipitation range of 1600–1800 mm and annual runoff depth is about 900 mm. The annual ETa generally shows increasing trends in the entire MRB during 2000–2019 (Figure 4b), except for slightly decreasing trends in individual grid points which indicate that ETa tends to increases of most region in this basin under the background of climate warming. ETa increasing rates above 4 mm·yr−1 pass the significance test at the 0.05 confidence level, while ETa decreasing trends are nonsignificant. Moreover, the spatial differences in annual ETa increase are obvious, with the most apparent increasing trend in the southwestern and southern parts of the MRB and weaker increasing trends in the central and northern parts.
The spatial change in the ETa seasonal variations in the MRB is obvious. There are generally significant increasing trends of ETa (p < 0.05) in winter and spring throughout the basin. Specifically, the increasing rates in winter range from 1 mm·yr−1 to 2 mm·yr−1, and the increasing rates are most significant in spring, with 2–4 mm·yr−1 in most areas (Figure 4b). The spatial differences in the ETa variations in summer and autumn are noticeable. In summer, there are increasing trends (p < 0.05) from 2000 to 2019 in the southwestern and southeastern regions of the basin, while trends in the central and northern regions are decreasing but nonsignificant (Figure 4). In autumn, the ETa tends to increase (0–1 mm yr−1) in the western and southern parts of this region (Figure 5) and decrease (−1–0 mm yr−1) in northeastern parts. However, none of these variation rates pass the significance test at the 0.05 confidence level.

3.3. Attribution Analysis of ETa Variations in the MRB

Table 4 shows that annual ETa was positively related to air temperature (r = 0.35) and negatively related to precipitation (r = −0.20). Annual wind speed had little effect on annual ETa in the region. ETa in winter was negatively related to precipitation (r = −0.35) and wind speed (r = −0.24). The spring ETa was controlled by air temperature (r = 0.60, p < 0.01). Summer ETa was controlled by precipitation (r = −0.52, p < 0.05) and air temperature (r = 0.45, p < 0.05). Autumn ETa was weakly negatively related to precipitation (r = −0.44, p < 0.1).
To quantify the effects of climatic variables on ETa in the MRB, we conduct PLS regression model. The PLS regression coefficient in Table 5 shows that precipitation, air temperature, and wind speed have seasonal differences. Annual air temperature has the largest effect on annual ETa (59.6%). ETa was mainly influenced by precipitation (64%) and wind speed (35.5%) in winter. ETa increased (1.10 mm yr−1, p < 0.05) in winter was due to increased precipitation and wind speed. In spring, ETa is dominated by air temperature (89.4%). ETa changes in summer and autumn are influenced by precipitation, with a relative contribution is 59.8% and 54.9%, respectively.

4. Discussion

Evapotranspiration is one of critical processes in land-atmosphere interactions. The extrapolation of the ETa estimation at a station to regional scale can result in significant errors. The reason for this is that the differences of underlying surfaces can lead to variations of meteorological elements, which consequently cause regional ETa differences [31]. With the improvement of the Earth system observation capability, especially the continuous advancement of high-precision remote sensing technology [32,33]. Moreover, with the development of land surface models [34], such as the GLDAS-Noah and Famine Early Warning Systems Network Land Data Assimilation System of the Noah-MP land surface model, they can provide the ETa from the land surface. Previous studies have extensively verified the applicability of GLDAS ETa data in China. The monthly GLDAS ETa had high stability and reliability in southwestern China [35]. In the Yellow River Basin, although there is little difference between GLDAS ETa and watershed ETa based on water balance, the overall uncertainty is high [36]. In the Weihe River Basin, the GLDAS ETa data products are fully meeting the needs of evapotranspiration research and have excellent application [37]. In the Yangtze River Basin, comparing the commonly used evapotranspiration data such as GLEAMV3.2a, MOD16, and GLDAS-Noah 2, suggests that the error of the GLDAS-Noah 2 dataset is small [38].
The factors of ETa variations can be divided into two categories: water limitation and energy limitation [2]. Water limitation mainly occurs in arid and semiarid areas, and energy limitation mainly occurs in humid climate areas [9]. ET change is more influenced by increased wind speed than the increased temperature in southern Italy [5]. In Canadian grasslands, wind speed significantly influences the decreasing evapotranspiration, while water vapor pressure difference (VPD) controls the increasing evapotranspiration [6]. In the Heihe River Basin of China, the most significant influencing factor on the daily variation of evapotranspiration is air humidity, followed by wind speed and soil moisture [7]. In the Jinsha River Basin, evapotranspiration is more sensitive to precipitation and temperature, followed by wind speed [8]. MRB is located in a subtropical monsoon zone, with an annual rainfall amount of about 1800 mm and an ETa of approximately 846 mm. The ETa variation in this region depends on energy limitation (Table 4), while a seasonal difference. In the context of global warming, the temperature in the MRB also tends to increase (0.3 °C·decade−1, Table 6) from 2000 to 2019, which is an indicator of the increase in the energy term at the higher temperature. Meanwhile, precipitation also increased (~18–19 mm yr−1) and wind speed insignificantly changed in the MRB. Therefore, combining Table 4, Table 5 and Table 6, the increase in ETa in winter and spring in the MRB is mainly due to the rise in temperature and wind speed. The increase in ETa was insignificant in summer, and it may be due to precipitation increases suppressed the effect of temperature rise on ETa. There was a weakly decreasing trend of ETa in autumn, and although the temperature increased, wind speed did not change in autumn.
In addition, the input parameters of the current GLDAS-Noah model do not include land-use types, which is an important reason that the GLDAS-Noah model underestimates the ETa in the urban region. Therefore, our subsequent research will investigate the influence of vegetation changes on ETa in the MRB.

5. Conclusions

In this study, we evaluated the applicability of the GLDAS-Noah ETa data and then used the GLDAS-Noah products to analyze the ETa variations in the MRB from 2000 to 2019. The main conclusions are as follows.
The accuracy of the GLDAS-Noah ETa data was assessed by using the calculated results from the P-M equation in 2000–2019 at eight stations of the MRB. The assessment results showed that the GLDAS-Noah ETa data had applicability in the MRB. The R2 value was close to 1, the NSE value was larger than 0.8, and the DISO was less than 0.3, which indicated that the GLDAS-Noah ETa data could be used to analyze the ETa variations in the MRB.
Since 2000, there have been significant spatio-temporal differences in the ETa variations in the MRB. During 2000–2019, ETa in the MRB shows an increasing trend, and the overall increasing rate is 3.60 mm·yr−1 (p < 0.01). The ETa showed a significant increase in winter and spring, with increasing rates of 1.10 mm·yr−1 (p < 0.01) and 2.60 mm·yr−1 (p < 0.01), respectively, while the ETa did not change significantly in other seasons. Spatially, the ETa tended to increase in winter and spring in the whole MRB. Therefore, we should pay more attention to the relationship between increased winter and spring ETa and the increased frequency of winter and spring droughts in the MRB, as well as the seasonal contradiction between supply and demand of water resources in the basin.
In the MRB, annual ETa was positively related to air temperature (r = 0.35) and negatively related to precipitation (r = −0.20). There are seasonal differences in precipitation, air temperature, and wind speed that affect ETa. ETa was mainly influenced by precipitation (64%) and wind speed (35.5%) in winter. In spring, ETa is dominated by air temperature, and the relative contribution rate of 89.4%.

Author Contributions

Conceptualization, H.D. and Q.L.; methodology, Y.L., Y.W. and H.D.; validation, X.C., Y.Z. and L.G.; formal analysis, Q.L., Y.C. and M.L.; writing—original draft preparation, Y.L. and Y.W.; writing—review and editing, H.D., Q.L., X.C., Y.Z., L.G., Y.C. and M.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Projects for National Natural Science Foundation of China (41877167), and the National Key Research and Development Program (2018YFE0206400).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets generated and analyzed during the current study are not publicly available but are available from the corresponding author on reasonable request.

Acknowledgments

The air temperature, precipitation, air pressure, wind speed and relative humidity data are supported by the China Meteorological Data Service Center. The authors appreciate the comments and encouragement given by the reviewers, editor, and associate editor.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gan, G.; Liu, Y.; Sun, G. Understanding interactions among climate, water, and vegetation with the Budyko framework. Earth-Sci. Rev. 2021, 212, 103451. [Google Scholar] [CrossRef]
  2. Roderick, M.L.; Farquhar, G.D. The Cause of Decreased Pan Evaporation over the Past 50 Years. Science 2002, 298, 1410–1411. [Google Scholar] [CrossRef] [PubMed]
  3. Beven, K. A sensitivity analysis of the Penman-Monteith actual evapotranspiration estimates. J. Hydrol. 1979, 44, 169–190. [Google Scholar] [CrossRef]
  4. Zeng, R.; Cai, X. Assessing the temporal variance of evapotranspiration considering climate and catchment storage factors. Adv. Water Resour. 2015, 79, 51–60. [Google Scholar] [CrossRef]
  5. McCabe, M.F.; Rodell, M.; Alsdorf, D.E.; Miralles, D.G.; Uijlenhoet, R.; Wagner, W.; Lucieer, A.; Houborg, R.; Verhoest, N.E.C.; Franz, T.E.; et al. The future of Earth observation in hydrology. Hydrol. Earth Syst. Sci. 2017, 21, 3879–3914. [Google Scholar] [CrossRef] [Green Version]
  6. Bowen, I.S. The ratio of heat losses by conductions and by evaporation from any water surface. Phys. Rev. 1926, 27, 779–787. [Google Scholar] [CrossRef] [Green Version]
  7. Shuttleworth, W.J.; Wallace, J.S. Evaporation from sparse crops-an energy combination theory. Q. J. R. Meteorol. Soc. 1985, 111, 839–855. [Google Scholar] [CrossRef]
  8. Zhang, B.; Xu, D.; Liu, Y.; Chen, H. Review of multi-scale evapotranspiration estimation and spatio-temporal scale expansion. Trans. Chin. Soc. Agric. Eng. 2015, 31, 8–16, (Chinese with English Abstract). [Google Scholar]
  9. Liu, Y.; Qiu, G.; Zhang, H.; Yang, Y.; Zhang, Y.; Wang, Q.; Zhao, W.; Jia, L.; Ji, X.; Xiong, Y.; et al. Shifting from homogeneous to heterogeneous surfaces in estimating terrestrial evapotranspiration: Review and perspectives. Sci. China Earth Sci. 2022, 65, 197–214. [Google Scholar] [CrossRef]
  10. Guerschman, J.P.; McVicar, T.R.; Vleeshower, J.; Van Niel, T.G.; Peña-Arancibia, J.L.; Chen, Y. Estimating actual evapotranspiration at field-to-continent scales by calibrating the CMRSET algorithm with MODIS, VIIRS, Landsat and Sentinel-2 data. J. Hydrol. 2021, 605, 127318. [Google Scholar] [CrossRef]
  11. Yin, Y.; Wu, S.; Dai, E. Determining factors in potential evapotranspiration changes over China in the period 1971–2008. Chin. Sci. Bull. 2010, 55, 3329–3337. [Google Scholar] [CrossRef]
  12. Wen, S.; Jiang, T.; Li, X.; Wang, T.; Wang, Y.; Thomas, F. Changes of actual evapotranspiration over the Songhua River Basin from 1961 to 2010. Adv. Clim. Chang. Res. 2014, 10, 79–86, (Chinese with English Abstract). [Google Scholar]
  13. Mo, X.; Chen, X.; Hu, S.; Hu, S.; Liu, S.; Xia, J. Attributing regional trends of evapotranspiration and gross primary productivity with remote sensing: A case study in the North China Plain. Hydrol. Earth Syst. Sci. 2017, 21, 295–310. [Google Scholar] [CrossRef] [Green Version]
  14. Liu, Y.; Zhou, M.; Chen, Z.; Li, S. Sensitivity of the potential evapotranspiration to climate and vegetation in Hanjiang River basin based on S-W Model. Trans. Chin. Soc. Agric. Eng. 2013, 29, 92–100, (Chinese with English Abstract). [Google Scholar]
  15. Zhao, Y.; Zou, X.; Zhang, J.; Cao, L.; Xu, X.; Zhang, K.; Chen, Y. Spatio-temporal variation of reference evapotranspiration and aridity index in the Loess Plateau Region of China, during 1961–2012. Quat. Int. 2014, 349, 196–206. [Google Scholar] [CrossRef]
  16. Ma, N.; Zhang, Y. Increasing Tibetan Plateau terrestrial evapotranspiration primarily driven by precipitation. Agric. For. Meteorol. 2022, 317, 108887. [Google Scholar] [CrossRef]
  17. Zhao, Y.; Zou, X.; Cao, L.; Yao, W.; Fu, J. Spatiotemporal variations of potential evapotranspiration and aridity index in relation to influencing factors over Southwest China during 1960–2013. Theor. Appl. Climatol. 2018, 133, 711–726. [Google Scholar] [CrossRef]
  18. Li, Z.; Chen, Y.; Shen, Y.; Liu, Y.; Zhang, S. Analysis of changing pan evaporation in the arid region of Northwest China. Water Resour. Res. 2013, 49, 2205–2212. [Google Scholar] [CrossRef]
  19. Liu, C.; Zhang, D. Temporal and spatial change analysis of the sensitivity of potential evapotranspiration to meteorological influencing factors in China. Sci. Geogr. Sin. 2011, 66, 579–588, (Chinese with English Abstract). [Google Scholar]
  20. Chen, Y.; Chen, X.-W.; Yin, Y.-X. Characteristics of Runoff Changes in the Minjiang River Basin from 1960 to 2006. J. Nat. Resour. 2011, 26, 1401–1411. [Google Scholar]
  21. Huang, W.; Sui, Y.; Yang, X.; Wei, D.; Song, L. Characteristics and adaptation of seasonal drought in southern China under the background of climate change.III. Spatiotemporal characteristics of seasonal drought in southern China based on the percentage of precipitation anomalies. Chin. J. Appl. Ecol. 2013, 24, 397–406, (Chinese with English Abstract). [Google Scholar]
  22. Rodell, M.; Houser, P.R.; Jambor, U.; Gottschalck, J.; Bosilovich, M. The Global Land Data Assimilation System. Bull. Am. Meteorol. Soc. 2004, 85, 381–394. [Google Scholar] [CrossRef] [Green Version]
  23. Zhao, L.; Ji, X. Quantification of transpiration and evaporation over agricultural field using the FAO-56 dual crop coefficient approach—A case study of the maize field in an oasis in the middlestream of the Heihe River Basin in northwest China. Sci. Agric. Sin. 2010, 43, 4016–4026, (Chinese with English Abstract). [Google Scholar]
  24. Wright, J.L. New evapotranspiration crop coefficients. Proc. Am. Soc. Civ. Eng. J. Irrig. Drain. Div. 1982, 108, 57–74. [Google Scholar] [CrossRef]
  25. Baburao, K.; Ayse, K.; Kenneth, H. Estimating Crop Coefficients Using Remote Sensing-Based Vegetation Index. Remote Sens. 2013, 5, 1588–1602. [Google Scholar] [CrossRef] [Green Version]
  26. Rocha, J.; Perdigo, A.; Melo, R.; Henriques, C. Managing water in agriculture through remote sNSEing applications. In Proceedings of the 30th EARSeL Symposium on Remote Sensing for Science, Education, and Natural and Cultural Heritage, Paris, France, 31 May–3 June 2010; Volume 31, p. 223230. [Google Scholar] [CrossRef]
  27. Hu, Z.; Chen, X.; Zhou, Q.; Chen, D.; Li, J. DISO: A rethink of Taylor diagram. Int. J. Climatol. 2019, 39, 2825–2832. [Google Scholar] [CrossRef]
  28. Hu, Z.; Zhang, Z.; Sang, Y.-F.; Qian, J.; Feng, W.; Chen, X.; Zhou, Q. Temporal and spatial variations in the terrestrial water storage across Central Asia based on multiple satellite datasets and global hydrological models. J. Hydrol. 2021, 596, 126013. [Google Scholar] [CrossRef]
  29. Zhang, S.H.; Liu, S.X.; Mo, X.G.; Shu, C.; Sun, Y.; Zhang, C. Assessing the impact of climate change on potential evapotranspiration in aksu river basin. J. Geogr. Sci. 2011, 21, 609–620. [Google Scholar] [CrossRef]
  30. Huang, Q.; Qiu, Y.; Chen, X.; Sun, X.; Lin, Y. Climate Change Characteristics and Future Trends of Fujian Province during 1980—2017. J. Wuyi Univ. 2020, 39, 27–32, (Chinese with English Abstract). [Google Scholar]
  31. Ramillien, G.; Frappart, F.; Güntner, A.; Ngo-Duc, T.; Cazenave, A.; Laval, K. Time variations of the regional evapotranspiration rate from Gravity Recovery and Climate Experiment (GRACE) satellite gravimetry. Water Resour. Res. 2006, 42, 10403-1–10403-8. [Google Scholar] [CrossRef] [Green Version]
  32. Mu, Q.; Zhao, M.; Running, S.W. Improvements to a MODIS global terrestrial evapotranspiration algorithm. Remote Sens. Environ. 2011, 115, 1781–1800. [Google Scholar] [CrossRef]
  33. Rodell, M.; Famiglietti, J.S.; Chen, J.; Seneviratne, S.I.; Viterbo, P.; Holl, S.; Wilson, C.R. Basin scale estimates of evapotranspiration using GRACE and other observations. Geophys. Res. Lett. 2004, 31, L20504. [Google Scholar] [CrossRef] [Green Version]
  34. Di, L.; Longuevergne, L.; Scanlon, B.R. Uncertainty in evapotranspiration from land surface modeling, remote sensing, and GRACE satellites. Water Resour. Res. 2014, 50, 1131–1151. [Google Scholar] [CrossRef] [Green Version]
  35. Li, H.; Wang, S.; Bai, X.; Tang, H.; Cao, Y.; Wang, M.; Wu, L. Inversion and Spatiotemporal Evolution of Actual Evapotranspiration in Southwest China for the Past 50 Years. Acta Ecol. Sin. 2018, 38, 8835–8848. [Google Scholar] [CrossRef]
  36. Chen, X.; Deng, C. Evaluation and Integration of Multiple Evapotranspiration Products in the Yellow River Basin. Water Resour. Power 2021, 39, 5–8. [Google Scholar]
  37. Wang, Y.; Wang, S.; Song, W.; Yang, S. Application of GLDAS Data to the Potential Evapotranspiration Monitoring in Weihe River Basin. J. Arid Land Resour. Environ. 2013, 27, 53–58. [Google Scholar] [CrossRef]
  38. Zhan, Y.; Zhang, W.; Yan, Y.; Wang, C.; Rong, Y.; Zhu, J.; Lu, H.; Zheng, T. Analysis of actual evapotranspiration evolution and influencing factors in the Yangtze River Basin. Acta Ecol. Sin. 2021, 41, 6924–6935. [Google Scholar] [CrossRef]
Figure 1. Topography and river network of the Minjiang River Basin (MRB).
Figure 1. Topography and river network of the Minjiang River Basin (MRB).
Atmosphere 13 00562 g001
Figure 2. Scatter plot of monthly precipitation from station observations and the GLDAS-Noah in the MRB.
Figure 2. Scatter plot of monthly precipitation from station observations and the GLDAS-Noah in the MRB.
Atmosphere 13 00562 g002
Figure 3. The compared analysis of monthly ETa between average of egiht stations P-M results and GLDAS-Noah data over MRB in 2000–2019 as (a) time series and (b) as a scatter plot.
Figure 3. The compared analysis of monthly ETa between average of egiht stations P-M results and GLDAS-Noah data over MRB in 2000–2019 as (a) time series and (b) as a scatter plot.
Atmosphere 13 00562 g003
Figure 4. Average annual ETa (a) and its trends (b) in the MRB from 2000 to 2019, and grid at 0.25° × 0.25°. The black dots in Figure 3b indicate that the variation rates passing the significance test (M-K trend test) at the 0.05 confidence level.
Figure 4. Average annual ETa (a) and its trends (b) in the MRB from 2000 to 2019, and grid at 0.25° × 0.25°. The black dots in Figure 3b indicate that the variation rates passing the significance test (M-K trend test) at the 0.05 confidence level.
Atmosphere 13 00562 g004
Figure 5. Seasonal trends of ETa in the MRB in (a) winter, (b) spring, (c) summer and (d) autumn from 2000–2019, and grid at 0.25° × 0.25°.
Figure 5. Seasonal trends of ETa in the MRB in (a) winter, (b) spring, (c) summer and (d) autumn from 2000–2019, and grid at 0.25° × 0.25°.
Atmosphere 13 00562 g005
Table 1. Basic information of the eight meteorological stations.
Table 1. Basic information of the eight meteorological stations.
NameLat. (°N)Lon. (°E)ELev. (m)Pre. (mm)Temp. (°C)LUCC
Shaowu27117218185518.3Grass land
Wuyishan27118222189918.7Forest land
Pucheng28118277175618.0Forest land
Jianyang27118197166318.6Cultivated land
Jianou27118155168219.4Forest land
Taining26117343180417.9Forest land
Nanping26118152161120.1Grass land
Note: ‘Elev.’ is elevation; ‘Pre.’ is average annual precipitation in 2000–2019; ‘Temp.’ is average annual precipitation in 2000–2019.
Table 2. Applicability assessment of the GLDAS ETa data.
Table 2. Applicability assessment of the GLDAS ETa data.
NameR2NSEDISO
Shaowu0.930.910.16
Wuyishan0.920.900.16
Pucheng0.930.820.29
Jianyang0.940.930.15
Jianou0.940.850.25
Taining0.930.850.23
Nanping0.920.800.26
Table 3. The annual and seasonal ETa trend (by M-K trend test) in the MRB from 2000 to 2019.
Table 3. The annual and seasonal ETa trend (by M-K trend test) in the MRB from 2000 to 2019.
Mean Value (mm)Trend (mm yr−1)
Annual846.23.60 **
Winter92.61.10 **
Spring207.52.60 **
Summer342.10.19
Autumn202.5−0.10
“**” indicate that the results passing the significance test at the 0.01 confidence level (p < 0.01).
Table 4. The Pearson correlation coefficients of between climate variables and ETa from 2000 to 2019.
Table 4. The Pearson correlation coefficients of between climate variables and ETa from 2000 to 2019.
PrecipitationTemperatureWind Speed
Annual ETa−0.200.35−0.08
Winter ETa−0.350.16−0.24
Spring ETa−0.100.60 ***−0.14
Summer ETa−0.52 **0.45 **0.26
Autumn ETa−0.44 *−0.21−0.10
Note: “***” indicate that the results passing the significance test at the 0.01 confidence level (p < 0.01), “**” denotes that the results passing the significance test at the 0.05 confidence level (p < 0.05), and “*” means that the results passing the significance test at the 0.1 confidence (p < 0.1).
Table 5. PLS regression coefficient and relative contribution rates of climatic variables over MRB in 2000–2019.
Table 5. PLS regression coefficient and relative contribution rates of climatic variables over MRB in 2000–2019.
bcd ρ P ρ T ρ W
Annual−0.200.35−0.0433.9%59.6%6.5%
Winter−0.310.002−0.1764%0.5%35.5%
Spring0.060.61−0.018.6%89.4%2.0%
Summer−0.400.170.1059.8%24.8%15.5%
Autumn−0.53−0.20−0.2354.9%21.0%24.1%
Table 6. Trend rate of ETa, precipitation, air temperature, and wind speed in the MRB from 2000 to 2019.
Table 6. Trend rate of ETa, precipitation, air temperature, and wind speed in the MRB from 2000 to 2019.
Precipitation (mm yr−1)Temperature (°C yr−1)Wind Speed (m yr−1)
Annual18.560.03 *−0.001
Winter0.400.010.005
Spring5.980.010.002
Summer6.390.020.001
Autumn0.370.06 **0.000
Note: “**” denotes that the results passing the significance test at the 0.05 confidence level (p < 0.05), and “*” means that the results passing the significance test at the 0.1 confidence (p < 0.1).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lu, Y.; Wang, Y.; Liu, Q.; Chen, X.; Zhang, Y.; Gao, L.; Chen, Y.; Liu, M.; Deng, H. Evapotranspiration Variations of the Minjiang River Basin in Southeastern China from 2000 to 2019. Atmosphere 2022, 13, 562. https://doi.org/10.3390/atmos13040562

AMA Style

Lu Y, Wang Y, Liu Q, Chen X, Zhang Y, Gao L, Chen Y, Liu M, Deng H. Evapotranspiration Variations of the Minjiang River Basin in Southeastern China from 2000 to 2019. Atmosphere. 2022; 13(4):562. https://doi.org/10.3390/atmos13040562

Chicago/Turabian Style

Lu, Yijin, Yuanyuan Wang, Qun Liu, Xingwei Chen, Yuqing Zhang, Lu Gao, Ying Chen, Meibing Liu, and Haijun Deng. 2022. "Evapotranspiration Variations of the Minjiang River Basin in Southeastern China from 2000 to 2019" Atmosphere 13, no. 4: 562. https://doi.org/10.3390/atmos13040562

APA Style

Lu, Y., Wang, Y., Liu, Q., Chen, X., Zhang, Y., Gao, L., Chen, Y., Liu, M., & Deng, H. (2022). Evapotranspiration Variations of the Minjiang River Basin in Southeastern China from 2000 to 2019. Atmosphere, 13(4), 562. https://doi.org/10.3390/atmos13040562

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