[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
Distribution Characteristics and Genesis Mechanism of Ground Fissures in Three Northern Counties of the North China Plain
Previous Article in Journal
Evaluation of Environmental Impacts in Legume Crops: A Case Study of PGI White Bean Production in Southern Europe
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Sustainability
Volume 16
Issue 18
10.3390/su16188025
sustainability-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Tracking Evapotranspiration Patterns on the Yinchuan Plain with Multispectral Remote Sensing

by
Junzhen Meng
1,
Xiaoquan Yang
1,
Zhiping Li
2,*,
Guizhang Zhao
1,
Peipei He
1,
Yabing Xuan
1 and
Yunfei Wang
1
1
College of Surveying and Geo-Informatics, North China University of Water Resources and Electric Power, Zhengzhou 450046, China
2
Henan Vocational College of Water Conservancy and Environment, The Education Department Henan Province, Zhengzhou 450046, China
*
Author to whom correspondence should be addressed.
Sustainability 2024, 16(18), 8025; https://doi.org/10.3390/su16188025
Submission received: 7 August 2024 / Revised: 10 September 2024 / Accepted: 11 September 2024 / Published: 13 September 2024
(This article belongs to the Section Sustainable Water Management)
Browse Figures
Figure 1
<p>Location of the meteorological stations and Landsat images that were used to illustrate the land cover conditions in the Yinchuan Plain area.</p> "> Figure 2
<p>Comparison between meteorological stations with large ET and geeSEBAL ET. Diamond-shaped points represent outliers lying outside the 150% inter-quartile range.</p> "> Figure 3
<p>Comparison between ET<sub>p</sub> and geeSEBAL ET (<b>a</b>–<b>d</b>). Compared with large-scale evapotranspiration, R<sup>2</sup> has significantly improved, indicating that the model correlation is influenced by external factors.</p> "> Figure 4
<p>Comparison between small ET and geeSEBAL ET ((<b>a</b>–<b>f</b>) presents different meteorological stations in the Yinchuan Plain).</p> "> Figure 5
<p>Seasonal ET<sub>a</sub> changes in the Yinchuan Plain. (<b>a</b>) Spring ET distribution; (<b>b</b>) Summer ET distribution; (<b>c</b>) Autum ET distribution; (<b>d</b>) Winter ET distribution.</p> "> Figure 6
<p>Trends in ET<sub>a</sub> on the Yinchuan Plain from 1987 to 2020.</p> "> Figure 7
<p>Area of Yinchuan Plain land use types.</p> "> Figure 8
<p>ET<sub>a</sub> in different subsurface types.</p> "> Figure 9
<p>Comparison of remote sensing imagery (<b>a</b>), land use classification (<b>b</b>), and ET imagery (<b>c</b>) on the Yinchuan Plain.</p> "> Figure 10
<p>Impervious areas misclassified in some intersecting land types. water bodies are identified as impervious areas (red color).</p> "> Figure 11
<p>ET<sub>a</sub> contribution of different subsurface types.</p> "> Figure 12
<p>geeSEBAL with the batch image estimation mode.</p> ">
Review Reports Versions Notes

Abstract

:
Evapotranspiration (ET) is a critical component of the hydrological cycle, and it has a decisive impact on the ecosystem balance in arid and semi-arid regions. The Yinchuan Plain, located in the Gobi of Northwest China, has a strong surface ET, which has a significant impact on the regional water resource cycle. However, there is a current lack of high-resolution evapotranspiration datasets and a substantial amount of time is required for long-time series remote sensing evapotranspiration estimation. In order to assess the ET pattern in this region, we obtained the actual ET (ETa) of the Yinchuan Plain between 1987 and 2020 using the Google Earth Engine (GEE) platform. Specifically, we used Landsat TM+/OLI remote sensing imagery and the GEE Surface Energy Balance Model (geeSEBAL) to analyze the spatial distribution pattern of ET over different seasons. We then reproduced the interannual variation in ET from 1987 to 2020, and statistically analyzed the distribution patterns and contributions of ET with regard to different land use types. The results show that (1) the daily ETa of the Yinchuan Plain is the highest in the central lake wetland area in spring, with a maximum value of 4.32 mm day−1; in summer, it is concentrated around the croplands and water bodies, with a maximum value of 6.90 mm day−1; in autumn and winter, it is mainly concentrated around the water bodies and impervious areas, with maximum values of 3.93 and 1.56 mm day−1, respectively. (2) From 1987 to 2020, the ET of the Yinchuan Plain showed an obvious upward and downward trend in some areas with significant land use changes, but the overall ET of the region remained relatively stable without dramatic fluctuations. (3) The ETa values for different land use types in the Yinchuan Plain region are ranked as follows: water body > cultivated land > impervious > grassland > bare land. Our results showed that geeSEBAL is highly applicable in the Yinchuan Plain area. It allows for the accurate and detailed inversion of ET and has great potential for evaluating long-term ET in data-scarce areas due to its low meteorological sensitivity, which facilitates the study of the regional hydrological cycle and water governance.

1. Introduction

Evapotranspiration (ET) is a complex physical process involving factors including vegetation transpiration, soil, and river and lake evaporation, playing an important role in the hydrological cycle and energy exchange in ecosystems such as the biosphere and the atmosphere [1]. Simultaneously, it is crucial for water resource management, improving water security, and achieving the allocation of water for agricultural crops. Studies have shown that 70% of the precipitation on the earth’s surface is returned to the atmosphere through ET [2]. Especially in agricultural ecosystems of arid and semi-arid zones, ET is an important factor in maintaining water and energy balance. In these areas, the evaporative losses account for more than 80% of the total water consumption [3]. The surface heat flux and water balance are critical factors in achieving water budget closure at the plot scale. Simulating their variations using remote sensing models is an essential method used for accurately quantifying the hydrological components of regional ecosystems [4,5], thereby contributing to a more comprehensive and precise assessment of local evapotranspiration. Remote sensing technology has been used to estimate surface evapotranspiration (ET) since the 1970s [6], effectively monitoring and quantifying water exchange processes [7]. In addition to measuring the daily surface ET using spatial and subsurface parameters, as well as the vegetation and soil information obtained from remote sensing [8,9,10,11], existing remote sensing data can also be utilized to predict ET variations; this approach helps mitigate issues related to the scarcity of observational data [12]. Models for estimating ET based on remote sensing can be mainly classified into three categories: statistical empirical methods, energy residual methods, and full remote sensing models.
Statistical empirical methods are easily affected by vegetation cover and weather conditions [13,14], and the full remote sensing model can achieve good inversion results in some arid areas. However, the inversion of ET in croplands or forest needs to be further researched and validated [15]. The residual method is used to invert the net radiative flux, soil heat flux and sensible heat flux by remote sensing technology to calculate the latent heat flux residual according to the surface energy balance equation, so as to derive the surface ET [16]. The residual method can be generally categorized into two types: single-layer models and double-layer models. The single-layer model does not distinguish the fluxes between soil and vegetation but take that as a unit for calculating the heat flux, and is commonly used to estimate the surface ET, which is strongly influenced by the land surface temperature (Ts). In cases of high vegetation cover, single-layer models exhibit greater accuracy in estimating land surface evapotranspiration. Additionally, the aerodynamic resistance parameters in the model can be determined based on surface characteristics and conventional meteorological observations. Furthermore, this single-layer model simplifies land surface processes and requires fewer parameters, meaning that it has broader applicability. Widely used single-layer models include SEBI (Surface Energy Balance Index) [17], SEBAL (Surface Energy Balance Algorithm for Land) [18,19,20,21,22], SEBS (Surface Energy Balance System) [23,24,25], and Simplified Surface Energy Balance (SSEBop) [26]. These models have different sensitivities to Ts and meteorological inputs. The SEBS model, which is an extension of SEBI, is more suitable for estimating the actual evapotranspiration over heterogeneous surfaces. The SSEBop requires high-quality daily meteorological aggregation data inputs [27,28,29], which has a large impact on the applicability of ET estimation, especially in areas lacking ground meteorological data. On the other hand, the SEBAL model is suitable for estimating the actual evapotranspiration over homogeneous vegetation cover, as it has a lower sensitivity to the meteorological data input and a higher sensitivity to Ts [30,31]. Additionally, the SEBAL model does not incorporate air temperature into the calculation of sensible heat, allowing for fewer parameters and internal adjustments, thus avoiding empirical adjustments to residual aerodynamic resistance. This characteristic enhances its stability in long-term crop evapotranspiration estimation compared to other single-layer models, meaning that it has demonstrated good applicability and reliability in estimating crop ET in the Gansu Plain in China [32], Shahrekord Plain, Chaharmahal and Bakhtiari Province, Iran [33], and Habra Plain in western Algeria [34]. The double-layer models proposed by Shuttleworth [35,36,37] in 1988 not only consider the energy exchange at the surface, but also incorporate vegetation characteristics and soil heat conduction. Compared to single-layer models, this kind of model can more accurately describe the surface fluxes in sparse vegetation areas, integrate geometric optics and consider the effects of variations in sensor observation angles, and avoid empirical corrections for residual impedance. However, the double-layer model is not extensively utilized because its impedance calculations employ empirical methods, making them more sophisticated and vulnerable to geo-graphical limitations.
The Yinchuan Plain, one of the three major alluvial plains of the Yellow River (including the Yinchuan Plain, Hetao Plain, and North China Plain), is characterized by its flat terrain, fertile soil, and convenient irrigation conditions, which provide significant advantages for local agricultural development. The hydrological cycle is a crucial component of agricultural production. However, the underdeveloped economic conditions and the absence of comprehensive water resource monitoring infrastructure contribute to a lack of high-precision hydrological observation data. Moreover, few studies have applied remote sensing methods to invert ET using high-resolution data series. Meanwhile, due to the complexity of traditional remote sensing methods, surface ET inversions have been conducted in 5-t10-year intervals and without spatial and temporal continuity. To compensate for this issue, we utilized a new tool called geeSEBAL that was developed based on the GEE platform and SEBAL for the estimation of surface ET, using Landsat5(TM+) and Landsat8(OLI) remote sensing images. The state-of-the-art meteorological reanalysis ERA5-Land data were used to obtain the multi-year actual evapotranspiration (ETa) in the Yinchuan Plain area from 1987 to 2020. The main objective was to analyze the changes in the ET of the Yinchuan Plain area over different seasons, years and land types.

2. Materials and Methods

2.1. Study Area

The Yinchuan Plain is situated in the central part of the Ningxia Hui Autonomous Region, flanked by the Yellow River on both sides. Bounded by the Ordos Plateau to the east, the Helan Mountains to the west, the Loess Plateau to the south and Shizuishan to the north, the study area lies between approximately 37°29′ to 38°53′ N and 105°49′ to 106°53′ E (Figure 1). Although the plain is situated in a temperate arid zone, the annual average precipitation reaches merely 200 mm. The annual water discharge of the Yellow River passing through the region reaches an average of 300 × 108 m³, making it an important industrial and agricultural hub in the middle and upper reaches of the Yellow River, and an important area for the national economic development of the Ningxia Hui Autonomous Region and surrounding areas. The annual evaporation of the Yinchuan Plain is as high as 1600 mm, which has a significant impact on the ecology of the Yinchuan Plain. This area relies heavily on the river runoff supply, which can become an important restricting factor the development of local industries and agriculture [38].

2.2. Data Sources

Since the Yinchuan Plain covers a wide area of 7790 km2, it requires high-quality remote sensing images. In this study, some 46 high-quality Landsat5(TM+) and Landsat8(OLI) images from between 1989 and 2020 were selected to cover the whole study area. These images were employed to estimate the ETa from the surface parameters. The ERA5 Land Reanalysis dataset was used as the meteorological data input for geeSEBAL [39,40]. This dataset contains hourly data on air temperature at 2 m, the dew point temperature, and the eastward and northward wind speed at 10 m. The relative humidity estimation follows Shuttleworth [41], and more detailed information on data input can be found in [42].
The meteorological data were collected from Ningxia Meteorological Bureau (Table 1) and used to estimate the reference ET based on the FAO–Penman–Monteith method [43]. Due to the lack of flux tower observations, the large and small ET data from meteorological stations were used as geeSEBAL ET validation data. The large evaporation pan was made from a 60 cm wide fiberglass plastic drum buried in the ground, and the evapotranspiration was close to the surface ETa. The small evaporator (20 cm in diameter) was made from galvanized iron or other alloys. In dry areas or dry seasons, the outer wall of the pan had higher temperatures given its small size, which makes the observed values significantly larger than the real water surface evaporation. Therefore, the obtained values cannot accurately represent the real free water surface evaporation. Likewise, estimating and tracking the actual land surface evaporation from these values become challenging, but they have certain reference value for understanding the changes and trends in water surface evaporation over time [44]. More detail about the evaporator used can be found in the Specification for Surface Meteorological Observations [45].
The China Land Cover Dataset (CLCD), used to calculate the ET for different subsurfaces, is provided by Xin Huang from Wuhan University [46], who produced the annual CLCD based on 335,709 Landsat images from Google Earth Engine. This dataset contains annual land cover information for China from 1990 to 2020. The author elicited spatiotemporal features based on all available Landsat data on GEE using the Random Forest Classifier. The author also proposed a post-processing method including spatiotemporal filtering and logical inference to further improve the consistency of the CLCD. Finally, the overall accuracy of CLCD reaches 80% based on 5463 visual decoding samples, meeting the accuracy requirements for this experiment.

2.3. Methods

2.3.1. The geeSEBAL Algorithm

SEBAL is a commonly used large-scale remote sensing ET assessment model that uses the energy of longwave and shortwave radiation from the sun and the atmosphere to estimate soil and surface-layer heating and ET. It divides the net surface radiative fluxes into three components: soil heat flux, sensible heat flux, and latent heat flux. The energy balance equation is as follows:
λ E = R n − H − G
where λ E is the latent heat flux (W/m²), R n is the surface net radiation flux (W/m²) (Equation (2)), H is the sensible heat flux (W/m²) and G is the soil heat flux (W/m²) (Equation (3)).
Net surface radiation (netradiation) is the main source of energy for the process of heat flux transfer and exchange in surface water. Soil heat flux is the state of heat exchange between the surface and deeper soil layers.
R n = ( 1 − α ) R S ↓ + R L ↓ − R L ↑ − ( 1 − ε 0 ) R L ↓
G R n = α ( T s − 273.15 ) ( 0.0038 α + 0.0074 α 2 ) ( 1 − 0.98 N D V I 4 )
where α is the surface albedo, calculated according to Tasumi and Ke [38,47], T s is the surface temperature, R S ↓ is the incoming shortwave radiation, R L ↓ is the incoming longwave radiation, R L ↑ is the outgoing longwave radiation, and ε 0 is the surface thermal emissivity.
The daily ET is computed by upscaling instantaneous λ E to E T a ( ∧ is the intermediate variable) (Equation (4)).
E T a 24 h = ∧ R n 24 h λ
where E T a 24 h is the daily ET (mm); R n 24 h is the daily average net radiation (W/m2); ∧ is the evaporative fraction; and λ is the latent heat of water evaporation (J/Kg).
The evaporative fraction ∧ is the intermediate variable.
∧ = λ E R n − G
where ∧ is the evaporative fraction; λ E is the latent heat flux (W/m²); R n is the net radiation (w/m²); and G is the soil heat flux (w/m²).
The SEBAL algorithm used in this paper was implemented in the GEE platform using the JavaScript language [42]. The GEE platform provides well-established Landsat datasets and meteorological data such as ERA5. ERA5 is the fifth generation of the ECMWF (European Centre for Medium-Range Weather Forecasts) atmospheric reanalysis dataset of the global climate from January 1950 until now. This dataset provides hourly estimates of a large number of atmospheric, terrestrial and oceanic climate variables, which allows geeSEBAL to handle ET on global to regional scales, especially in areas where ground meteorological data are limited. geeSEBAL provides long-term and comprehensive reference ET values. The operation process of this program can be found in the references [42].
It is worth noting that in the process of calculating the actual vapor pressure and relative humidity in the geeSEBAL model, the specific humidity is a crucial variable, but it is not archived directly in the ERA5 dataset in the GEE platform. However, the ERA5 contains the near-surface (2 m from the surface) temperature, dew point temperature, and surface pressure, from which one can calculate the specific humidity at 2 m. The specific humidity at 2 m q s a t was calculated according to Equations (6) and (7) (for more information, refer to [48]). In these equations, the dew point temperature T and surface pressure p (which is approximately equal to the pressure at 2 m) are approximated at the height of 2 m.
q s a t = R d r y R v a p e s a t ( T ) p − ( 1 − R d r y R v a p ) e s a t ( T )
where R d r y and R v a p are the gas constants for dry air and water vapor. The constants are defined in [42]. The saturation vapor pressure is expressed with Teten’s formula:
e s a t ( T ) = a 1 exp a 2 ( T − T 0 T − a 3 )
with the parameters set according to Buck (1981) [49] for saturation over water ( a 1 = 611.21 Pa, a 2 = 17.502 and a 3 = 32.19 K) and according to the AERKi formula of Alduchov and Eskridge [50] for saturation over ice ( a 1 = 611.21 Pa, a 2 = 22.587 and a 3 = −0.7 K), with T 0 = 273.16 K.

2.3.2. Hot and Cold Endmembers: Automated Calibration

There are two important assumptions in the SEBAL model: (i) there is a linear relationship between the geothermal temperature difference dT and the surface temperature Ts (Equation (8)); (ii) there are both cold and hot image elements in the study area, with the H of the cold and hot image elements being 0 and ( R n − G ), respectively.
d T = a + b ⋅ T s
where a and b are the empirical regression coefficients.
The first assumption serves to simplify the problem of solving the linear coefficients of the surface radiant temperature T s ; the second assumption, which introduces two cases of ET extremes (0 for the cold image element H, 0 for the hot image element LE), will thus solve for the coefficients of the linear relationship described above, and thus for H over the entire region.
geeSEBAL supports a simplified automated statistical algorithm to select the hot and cold endmembers, which is based on the version of the Calibration using Inverse Modeling at Extreme Conditions (CIMEC) algorithm used in METRIC [51]. The NDVI and T are usually used as the cold and hot endmember candidates, for example, in a semi-arid climate. Here, the cold endmember is selected within a set of candidates with the highest NDVI (5%) and the lowest Ts (20%) percentiles. Conversely, the hot endmember is selected based on the lowest NDVI (10%) and the highest Ts (20%) percentiles [51]. In this study, we performed the selection of hot and cold end members based on the NDVI and Ts of each image according to the thresholds provided in Table 2 and Table 3.

2.3.3. Statistical Analysis

RMSD (Root Mean Square Deviation) and the bias value were used to evaluate the consistency between the ET estimated by geSEBAL and the ET measured by meteorological stations (Equations (9) and (10)). R2 is used as the correlation analysis between ETa and the potential evapotranspiration by PM model to comprehensively evaluate the applicability of geeSEBAL in the Yinchuan Plain.
R M S D = ∑ i = 1 n ( Y i − Y − i ) 2 n
B i a s = ∑ i = 0 n ( Y i − Y − i ) n

2.3.4. Mann–Kendall Trend Analysis

The Mann–Kendall (MK) nonparametric test was used in this study to reveal the trends in the ET changes in the Yinchuan Plain area. The MK method has no requirements regarding the distribution of the sample data, which works well for not normally distributed hydrological data [52]. The MK method defines the statistic S as Equation (11):
S = ∑ j = 1 n − 1 ∑ k = j + 1 n s ( Y k − Y j )
where Y 1 , Y 2 …, Y n are the time series variables and n is the length of the time series.
s ( Y k − Y j ) = + 1 ,   Y k > Y j 0 ,   Y k = Y j − 1 ,   Y k < Y j
where Y j and Y k are the corresponding measurements for the years j and k, respectively, and k > j.
Z = S − 1 V ( S ) ,   S > 0 0 ,   S = 0 S + 1 V ( S ) ,   S < 0
where Z is a normally distributed statistic and V ( S ) is the variance. At a given α confidence level, Z > 0 indicates an upward trend and Z < 0 indicates a downward trend. The statistics Z ≥ 1.96 and Z ≥ 2.58 mean significance at 95% and 99% confidence levels, respectively.

3. Results

3.1. geeSEBAL ET Validation

Four meteorological sites with longer time series and large ET values were selected for validation (2000–2020) (Figure 2). Meanwhile, to ensure the accuracy of the long ET time series, meteorological stations with relatively complete and small ET values from 1990 to 2010 were selected as geeSABAL validation sites (Figure 3).
The comparison of geeSEBAL ET and the observed values (Figure 2) showed moderate-to-high accuracy in the Yinchuan plain area. The average bias was estimated as −0.11 mm day−1 close to 0, ranging from −0.87 mm day−1 to 0.80 mm day−1. The average RMSD was measured as 1.12 mm day−1 without a clear bias pattern, ranging from 0.83 to 1.47 mm day−1, which demonstrates a wide range of applicability.
The comparison of geeSEBAL values with the large evapotranspiration values at some of the meteorological stations, such as meteorological stations 53615 and 53519 distributed along the Yellow River, showed lower R2 values. Station 53615 had the highest RMSD and bias at 1.47 mm day−1 and −0.87 mm day−1, respectively, and the lowest R2 at 0.12. To further investigate the reasons for the low R2 values, we calculated the potential evapotranspiration using the PM model based on local meteorological data (Figure 3) and the result showed a significant improvement in R2. Additionally, the R2 values for large evapotranspiration and the potential evapotranspiration from the geeSEBAL results are consistent with stations 53518, 53519, and 53614, and station 53615 exhibited the lowest R2 for both large and potential evapotranspiration. This indicates that these results are not due to the geeSEBAL model itself, but rather because the ET calculated is influenced by the meteorological or image conditions affecting the actual ET. On the one hand, the Landsat remote sensing images were scarce in dry seasons. On the other hand, the limited number of years with actual large ET values resulted in large time spans between the data. Although an R2 was not as high as expected in four meteorological stations when comparing large ET values, the R2 showed a higher correlation for most of the stations with small ET values, reaching up to 0.82 (Figure 4). As for stations 53614 and 53615, there was a good consistency between the large and small ETs. Although a weak correlation was detected between the geeSEBAL and observed ET values, the estimated values were still within a reasonable range, with an RMSD of 1.47 mm day−1.
Due to the large span of the Yinchuan Plain and its complex geomorphology, the surface ET results for the remote sensing inversion are uncertain. This paper analyzes the inversion results from the geeSEBAL and PM based on quantitative aspects. The surface ET obtained by inversion based on the geeSEBAL model was analyzed in correlation with the values obtained from the potential ET(ETp) calculated by the P-M model (Figure 3). Only in site 53619 was a weak correlation with an R2 of 0.43 detected. The rest of the stations showed strong correlations, with the R2 ranging from 0.62 to 0.79. In summary, the surface ET inversion accuracy on the Yinchuan Plain using the geeSEBAL model was relatively high, indicating that the model has high applicability in arid and semi-arid areas.
Sustainability 16 08025 g003
Figure 3. Comparison between ETp and geeSEBAL ET (a–d). Compared with large-scale evapotranspiration, R2 has significantly improved, indicating that the model correlation is influenced by external factors.
Figure 3. Comparison between ETp and geeSEBAL ET (a–d). Compared with large-scale evapotranspiration, R2 has significantly improved, indicating that the model correlation is influenced by external factors.
Sustainability 16 08025 g003
Sustainability 16 08025 g004
Figure 4. Comparison between small ET and geeSEBAL ET ((a–f) presents different meteorological stations in the Yinchuan Plain).
Figure 4. Comparison between small ET and geeSEBAL ET ((a–f) presents different meteorological stations in the Yinchuan Plain).
Sustainability 16 08025 g004

3.2. Changing Patterns of Evapotranspiration in Yinchuan Plain

3.2.1. Changing Trend of ET

The geeSEBAL was driven by ERA5 Land reanalysis using the optimized percentiles and the standard percentiles recommended by Allen et al. [48]. Table 2 shows the dates of the images in four seasons in the Yinchuan Plain. From the spatial distribution of the spring ET in Figure 5a, it can be seen that the northern and southern regions have lower ET values in the range of 0–3 mm day−1. The ET is greater at water bodies and around Yinchuan City, with a maximum value of 4.32 mm day−1. This is mainly due to the existence of broad lakes and wetlands around Yinchuan City, which is known as the “International Wetland City”. There are five national wetland parks around this city, with the Yinchuan Plain wetland covering more than 80% of the total area and acting as the main source of spring ET in this region.
Figure 5b shows the spatial distribution of ET on the Yinchuan Plain in summer. The ET values from bare land and impervious areas are relatively lower, ranging between 0 and 4 mm day−1. These lower ET values are mainly due to the presence of impervious areas. From the ET spatial distribution, it can be seen that water bodies and arable lands are the main sources of ET in summer, with the highest value being 6.90 mm day−1.
Figure 5c shows that the ET in the fall is more uniformly distributed. In this case, water bodies have the highest ET, at 3.93 mm day−1. Leaving crop lands fallow has resulted in lower ET values in the range of 1–3 mm day−1. Figure 5d shows that the highest ET of 1.56 mm day−1 occurred in winter. Although this result is affected by the quality of the image, it can be seen that the ET in water bodies and impervious areas is generally higher than that in cultivated land and bare land.
In summary, the seasonal ET variation on the Yinchuan Plain is concentrated in agricultural fields, wetlands, rivers, and lakes during spring and summer, and primarily in water bodies like rivers and lakes during autumn and winter. This ET distribution pattern indicates that local agricultural and water resources are relatively abundant. The local government should focus on protecting and developing the ecological environment with regard to these two aspects to promote the sustainable development of key industries.
To obtain long-term ET trends in the Yinchuan Plain area, a total of 22 representative images were selected. These images, acquired from June to August between 1989 and 2020, offer a complete coverage of the study area (Table 3). The trend in the surface ET was obtained using MK trend analysis (Figure 6).
The overall ET changes on the Yinchuan Plain are mostly concentrated in the western section, impervious areas and along the Yellow River, with confidence intervals exceeding 99%. ET changes in the vegetated areas (arable land and grassland) are not obvious, and the majority of the confidence intervals are lower than 95%. The ET in the western region and some areas along the Yellow River showed an increasing trend. However, in several major residential areas (Shizuishan City, Yinchuan City, Wuzhong City), it showed a decreasing trend. Our results showed that the areas with significant ET alterations have undergone considerable land use changes. Land use changes from bare lands to croplands, wetlands and water bodies resulted in a decreasing ET trend. At the same time, land use conversion from cropland to impervious areas resulted in an increasing ET trend. The trend in ET reflects the land use changes to a certain extent. In order to further explore the distribution of ET in the Yinchuan Plain, the contribution of different land uses to ET values is discussed in detail in Section 3.2.2.

3.2.2. Contribution of Soil Evaporation and Transpiration to ET Change

The variation in evapotranspiration (ET) is influenced by meteorology, the vegetation cover type, and other factors. Additionally, there are significant disparities in the physical and chemical properties of different land use types on the subsurface. These disparities result in diverse responses on the surface to solar radiation. Furthermore, there can be substantial differences in the surface radiation and energy distribution within similar subsurface areas, leading to the non-uniform distribution of ET across the land surface. In order to investigate the ETa changes on different subsurfaces of the Yinchuan Plain, the land surface ET from 1990 to 2020 was used along with the CLCD land use dataset. The land use types on the Yinchuan Plain could be grouped into five main categories: cropland, grassland, water bodies, bare land and impervious land (Figure 7).
As shown in Figure 8, the daily ET of different land use types is in the following order: water bodies > cropland > impervious > grassland > bare land. The ET of different types showed obvious stratified characteristics and highly consistent trends. The fluctuation amplitude of each line reflects the response of different subsurface types to climatic factors. Among them, the ET of water bodies fluctuates between 5.03 and 7.62 mm day−1, being the most sensitive to the influence of meteorological factors such as temperature, wind speed, and surface radiation, etc. Cultivated land and grassland, as the main land use types on the Yinchuan Plain, show opposite trends, with the ET of cultivated land ranging from 4.16 to 6.36 mm day−1, and the ET of grassland ranging from 2.07 to 3.78 mm day−1. This is mainly due to low transpiration in grasslands. On the Yinchuan Plain, rice paddy is the main agricultural land use type where water is more abundant during the growing season in summer. Increased sunlight hours lead to increased evapotranspiration (ET) values, which are twice as large as those in grasslands. Bare lands and impervious areas are the two land use types with high change rates on the Yinchuan Plain. The conversion of bare lands into built-up and impervious areas has altered ET in this region. In this manner, the amount of ET in impervious areas ranges from 3.28 to 5.27 mm day−1, while it ranges from 0.81 to 2.67 mm day−1 in bare lands.
It should be noted that an anomalous phenomenon was observed, where the ET for impervious surfaces exceeded that of grasslands; to clarify the specific reason for this phenomenon in the study, we compared remote sensing imagery, land use classification data, and ET imagery for the region (Figure 9). Locally, instances exist where the ET of grasslands surpasses that of impervious surfaces, particularly in the eastern irrigation areas where grassland types are more prevalent. However, on a broader scale, the primary distribution of grasslands on the Yinchuan Plain differs from that in the eastern irrigation areas. Most grasslands are situated in the western arid region and are quite sparse. As depicted in Figure 3c, the ET from these grasslands is generally lower compared to the ET from impervious surfaces in the central irrigation area. Additionally, the land use types within impervious areas are heterogeneous, encompassing various vegetation types, wetlands, and water bodies. As illustrated in Figure 10, the CLCD dataset lacks finer classifications, which is another significant reason why the overall ET from impervious surfaces exceeds that from local grasslands.
Human activities often have significant impacts on natural systems and may even determine the evolution of regional microclimates and ecosystems in arid and semi-arid regions. For example, human activities have dominated the evolution of hydro-ecological processes in the lower Tarim River for centuries [53]. The most apparent manifestation of human activities in the Yinchuan Plain region could be regarded as land use conversion. In order to more intuitively respond to the distribution pattern of ET on the Yinchuan Plain and under different land use types, the contribution of ET ( σ ) of different types (Figure 9) was calculated using Equation (14).
σ C = E T C E T C + E T G + E T W + E T B + E T I
where σ C is the contribution of cropland to ET, and E T C − E T I is the evaportranpiration of different land use types.
As Figure 11 showed, before 1996, σ was in the order of cropland > grassland > bare land > water bodies > impervious areas. After 1996, the order changed to cropland > grassland > impervious areas > water body > bare land. Among the different vegetated land use types, croplands contributed between 68% and 78% to the total ET, followed by grassland contributing between 13% and 21%. Areas with sparse vegetation including water bodies, bare lands and impervious areas accounted for less than 10% of the total ET. In summary, croplands are the main source of evapotranspiration on the Yinchuan Plain. This implies that the impact of human activities has become one of the most important factors dominating evapotranspiration in the region. The primary focal points for hydrological cycle conservation and water governance in the future appear to be cropland and vegetation management.

4. Discussion

4.1. The Effects on ET Trends

In this paper, the spatiotemporal variation pattern of evapotranspiration on the Yinchuan Plain from 1987 to 2020 was studied from two aspects: one is the distribution of ET in different seasons, and the other is the trend in ET changes in long-term series. Comparing the summer ET distribution in Figure 5b with the summer variation trend in Figure 8, the maximum difference in ET variation among different land types can reach 5.46 mm day−1, while the maximum difference in ET variation among the same land type is only 2.59 mm day−1. This indicates that the variation in surface ET is greatly influenced by different underlying surfaces, and that human activities are one of the important factors affecting regional ET variation. Studies have shown that meteorological factors such as temperature and rainfall also have an important impact on ET. Xiaomei Jin et al. [54] analyzed the impact of the monthly average temperature and rainfall on ET on the Yinchuan Plain using actual data, and found an R2 of 0.82 between ET and temperature, while it was only 0.38 between ET and rainfall. At the same time, meteorological factors such as the sunshine hours, wind speed, and relative humidity also affect surface evapotranspiration to varying degrees [55,56]. For example, Wang et al. [57] found that a decrease in the local sunshine hours in the Yangtze River basin would reduce the energy source of surface evapotranspiration, thereby weakening the process of surface evapotranspiration. In summary, it can be seen that the factors affecting surface ET changes are complex and diverse. However, land use types cannot fully explain the changes in ET, and the influence of climate factors is also an important factor.
Generally speaking, climate factors are a prerequisite for surface energy supply, determining the maximum actual surface evapotranspiration. However, in arid and semi-arid environments, vegetation canopy serves as an interceptor of energy flux and hence determines the range of surface evapotranspiration [58]. As the main vegetation and land use type on the Yinchuan Plain, cropland accounts for about 60% of the total area. In Figure 1, the vegetation cover of the Yinchuan Plain during the growing season (June to August) is shown. In Figure 6, it can be seen that in the time series studied, the evapotranspiration changed only significantly in a small number of areas, such as bare land conversions into farmlands, wetlands, or water bodies. Most areas did not show significant changes in evapotranspiration within the study period, such as farmland. This is consistent with the findings of Xiaomei Jin et al. in their study of the Yinchuan Plain from 2001 to 2014. This indicates that the trend in ET changes in the Yinchuan Plain area is mainly dominated by changes in cropland and land cover types, with meteorological factors being secondary influences. Through studying the changes in ETa, it is possible to better understand the water consumption of vegetation and land, adjust the allocation of agricultural water resources in a timely manner (such as irrigation) and improve the efficiency of freshwater resource utilization (avoiding inefficient water use or irrigation water shortages). Meanwhile, the attribution analysis of surface evapotranspiration changes on the Yinchuan Plain will be the focus of subsequent research, clarifying the direction of water resource management and sustainable development in the region.

4.2. The Uncertainty Analysis of geeSEBAL

Because of the impact of observational conditions, the validation data gathered from weather stations were not continuous, and there were missing measurements for some study dates and years. Therefore, the large evapotranspiration values from four stations with a relatively complete dataset and the small evapotranspiration values from six stations were selected for validation. Based on our findings, the estimation results of the geeSEBAL were within a reasonable margin of error. Even if there were some stations with weak correlations, these would have been due to the closure mode of the surface energy [42]. According to Laipelt [42], the main challenge faced by the geeSEBAL program is the selection of end members and the definition of a spatial scope. The geeSEBAL model is highly sensitive to the spatial scope, and the calibration function of cold and hot end members can significantly improve the accuracy of ET estimation [59]. This study upgraded the estimation method of ET from single-image to batch estimation in order to achieve high consistency in ET inversion over a spatial range and time span (Figure 12). Therefore, the same calibration end members were set for multiple images on the same date. This affected the estimation accuracy in some areas to some extent, but the overall ET estimation showed satisfactory results (Figure 2). It can also be seen that geeSEBAL has a great potential for estimating ET over small areas, and adjusting the calibration end members based on the actual situation to achieve a higher estimation accuracy.
Due to the scarcity of high-quality remote sensing images in Google Earth Engine and the lack of observation data from meteorological stations, it was difficult to extrapolate ET to both monthly and annual time scales. Therefore, the ET trend analysis for long time series only considered the estimated ET representing dates as an example. By the same token, we conducted no further analysis of the sensitivity and contribution of meteorological factors to ET. This can be used in future research to fully comprehend the pattern of changes in evapotranspiration on the Yinchuan Plain.
Unlike previous studies on ET on the Yinchuan Plain by Wang Zhuoyue et al. [60], the GEE online dataset in this study can be used to quickly obtain continuous high-spatiotemporal-resolution ET values around the world. While the accuracy of meteorological datasets can impact the geeSEBAL model, this problem can be mitigated by using the automatic calibration option of endpoints or by restricting the study area. This ensures that the ET distribution maps generated fulfill the necessary criteria. The ET values for different land use types can also be quickly calculated based on the CLCD, avoiding the cumbersome image processing required in the traditional remote sensing techniques to estimate evapotranspiration. This greatly reduces the inversion time for regional evapotranspiration, and fully leverages the advantages of remote sensing technology in regional water resource management.

5. Conclusions

This study was conducted on the Yinchuan Plain, the irrigation area of the Yellow River Basin. The open-source SEBAL framework was used to verify the applicability of the model to the Yinchuan Plain through the Application Programming Interface (API) to estimate ET and to evaluate its distribution pattern across different underlying surfaces.
In spring, the daily ET of the Yinchuan Plain was found to be mainly concentrated in the central lake wetland area, with a maximum value of 4.32 mm day−1. In summer, the daily ET was found to be concentrated in the cultivated land and water areas, with a maximum value of 6.90 mm day−1. The ET values were highest in autumn and winter in water bodies and impervious areas, with maximum values of 3.93 mm day−1 and 1.56 mm day−1, respectively.
From 1987 to 2020, the ET of the Yinchuan Plain showed an obvious upward and downward trend in some areas with significant land use changes, but the overall ET of the region remained relatively stable without dramatic fluctuations.
The evapotranspiration of different underlying land use types on the Yinchuan Plain region is as follows: water body > cultivated land > impervious areas > grassland > bare land. The evapotranspiration of water bodies ranges from 5.03 to 7.62 mm day−1, that of cropland ranges from 4.16 to 6.36 mm day−1, that of impervious areas ranges from 3.28 to 5.27 mm day−1, that of grasslands ranges from 2.07 to 3.78 mm day−1, and that of bare land ranges from 0.81 to 2.67 mm day−1. Thid study shows that the distribution pattern of regional evapotranspiration is not only influenced by climate conditions, but it is also closely related to human activities.

Author Contributions

Conceptualization, J.M. and G.Z.; Methodology, J.M. and X.Y.; Software, P.H., Y.X. and Y.W.; Validation, X.Y.; Formal analysis, X.Y.; Data curation, X.Y. and G.Z.; Writing – original draft, X.Y.; Writing – review & editing, J.M.; Supervision, Z.L.; Funding acquisition, J.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, Zhiping Li (Grant number 41901285 and 41972261) and Training Plan for Young Backbone Teachers in Higher Education Institutions in Henan Province, Junzhen Meng (Grant number 2023GGJS073).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The data are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Priestley, C.H.B.; Taylor, R.H. On the assessment of surface heat flux and evaporation using large-scale parameters. Mon. Weather Rev. 1972, 100, 81–92. [Google Scholar] [CrossRef]
  2. Resenberg, N.J.; Blad, B.L.; Verma, S.B. Microclimate: The Biological Environment of Plants; John Wiley & Sons: New York, NY, USA, 1983. [Google Scholar]
  3. Chen, H.; Yang, D.W.; Lv, H.F. Comparison of temporal extrapolation methods for evapotranspiration over variant underlying croplands. Trans. Chin. Soc. Agric. Eng. 2013, 29, 73–81. [Google Scholar]
  4. Kahan, D.S.; Xue, Y.; Allen, S.J. The impact of vegetation and soil parameters in simulations of surface energy and water balance in the semi-arid Sahel: A case study using SEBEX and HAPEX-Sahel data. J. Hydrol. 2006, 320, 238–259. [Google Scholar] [CrossRef]
  5. Sridhar, V.; Wedin, D.A. Hydrological behaviour of grasslands of the Sandhills of Nebraska: Water and energy-balance assessment from measurements, treatments, and modelling. Ecohydrol. Ecosyst. Land Water Process Interact. Ecohydrogeomorphol. 2009, 2, 195–212. [Google Scholar] [CrossRef]
  6. Brown, K.W.; Rosenberg, N.J. A resistance model to predict evapotranspiration and its application to a sugar beet field 1. Agron. J. 1973, 65, 341–347. [Google Scholar] [CrossRef]
  7. Paniconi, C.; Putti, M. Physically based modeling in catchment hydrology at 50 Survey and outlook. Water Resour. Res. 2015, 51, 7090–7129. [Google Scholar] [CrossRef]
  8. Gao, Y.C.; Long, D. Progress in Models for Evapotranspiration Estimation Using Remotely Sensed Data. Natl. Remote Sens. Bull. 2008, 12, 515–528. [Google Scholar]
  9. Li, X.Y.; Yu, D.Y. Progress on Evapotranspiration Estimation Methods and Driving Forces in Arid and Semiarid Regions. Arid Zone Res. 2020, 37, 26–36. [Google Scholar]
  10. Zhang, R.H.; Sun, X.M.; Liu, J.Y. Regional variations in crop transpiration and soil water utilization by quantitative remote sensing inversion. Sci. Sin. Terrae 2001, 11, 959–968. [Google Scholar]
  11. Zhao, T.W.; Zhu, W.B.; Pei, L. Remote Sensing Estimation of Terrestrial Evapotranspiration and Analysis of Its Temporal-spatial Distribution Characteristics over the Three-River Headwater Region. Remote Sens. Technol. Appl. 2022, 37, 137–147. [Google Scholar]
  12. Jaksa, W.T.; Sridhar, V.; Huntington, J.L.; Khanal, M. Evaluation of the complementary relationship using Noah Land Surface Model and North American Regional Reanalysis (NARR) data to estimate evapotranspiration in semiarid ecosystems. J. Hydrometeorol. 2013, 14, 345–359. [Google Scholar] [CrossRef]
  13. Carlson, T.N.; Capehart, W.J.; Gillies, R.R. A new look at the simplified method for remote sensing of daily evapotranspiration. Remote Sens. Environ. 1995, 54, 161–167. [Google Scholar] [CrossRef]
  14. Seguin, B.; Itier, B. Using midday surface temperature to estimate daily evaporation from satellite thermal IR data. Int. J. Remote Sens. 1983, 4, 371–383. [Google Scholar] [CrossRef]
  15. Zhang, R.H.; Sun, X.M.; Zhu, Z.L. Differential thermal inertia-based full remote sensing information model of surface evaporation and its validation in Shapotou area, Gansu, China. Sci. Sin. Terrae 2002, 12, 1041–1050. [Google Scholar]
  16. Engman, E.T. Recent advances in remote sensing in hydrology. Rev. Geophys. 1995, 33, 967–975. [Google Scholar] [CrossRef]
  17. Menenti, M.; Choudhury, B.J. Parameterization of land surface evaporation by means of location dependent potential evaporation and surface temperature range. In Proceedings of the International Symposium, Yokohama, Japan, 13–16 July 1993. [Google Scholar]
  18. Bastiaanssen, W.G.M.; Menenti, M.; Feddes, R.A. A remote sensing surface energy balance algorithm for land (SEBAL). 1. Formulation. J. Hydrol. 1998, 212, 198–212. [Google Scholar] [CrossRef]
  19. Bastiaanssen, W.G.M.; Pelgrum, H.; Wang, J.; Al, E. A remote sensing surface energy balance algorithm for land (SEBAL). Part 2 Validation. J. Hydrol. 1998, 212, 213–229. [Google Scholar] [CrossRef]
  20. Jin, K.L.; Hao, L. Evapotranspiration estimation in the Jiangsu-Zhejiang-Shanghai Area based on remote sensing data and SEBAL model. Remote Sens. Nat. Resour. 2020, 32, 204–212. [Google Scholar]
  21. Gong, Z.N.; Lu, L.; Jin, D.D.; Ackermann, K. Remote sensing estimation of evapotranspiration and ecological water demand in Zhalong wetland under land use/cover change. Acta Ecol. Sin. 2021, 41, 3572–3587. [Google Scholar]
  22. Li, X.L.; Yang, L.X.; Xue, X.X. Analysis of evapotranspiration pattern by SEBAL model during the growing season in the agro-pastoral ecotone in Northwest China. Acta Ecol. Sin. 2020, 40, 2175–2185. [Google Scholar]
  23. Su, Z.; Schmugge, T.; Kustas, W.P. An Evaluation of Two Models for Estimation of the Roughness Height for Heat Transfer between the Land Surface and the Atmosphere. J. Appl. Meteorol. Clim. 2001, 40, 1933–1951. [Google Scholar] [CrossRef]
  24. Zhang, D.; Wei, Z.; Feng, T.; Chen, Z.; Wang, F.; Huang, K. Analysis of evapotranspiration estimation and its spatial-temporal characteristics: Taking Zhanghe irrigation district as an example. Bull. Surv. Mapp. 2022, 12, 57–63. [Google Scholar]
  25. Su, Z. The Surface Energy Balance System (SEBS) for estimation of turbulent heat fluxes. Hydrol. Earth Syst. Sci. 2002, 6, 85–100. [Google Scholar] [CrossRef]
  26. Senay, G.B.; Bohms, S.; Singh, R.K. Operational Evapotranspiration Mapping Using Remote Sensing and Weather Datasets a New Parameterization for the SSEB Approach. Jawra J. Am. Water Resour. Assoc. 2013, 49, 577–591. [Google Scholar] [CrossRef]
  27. Bhattarai, N.; Shaw, S.B.; Quackenbush, L.J. Evaluating five remote sensing based single-source surface energy balance models for estimating daily evapotranspiration in a humid subtropical climate. Int. J. Appl. Earth Obs. 2016, 49, 75–86. [Google Scholar] [CrossRef]
  28. Tasumi, M.; Trezza, R.; Allen, R.G.; Al, E. Operational aspects of satellite-based energy balance models for irrigated crops in the semi-arid US. Irrig. Drain. Syst. 2005, 19, 355–379. [Google Scholar] [CrossRef]
  29. Allen, R.; Irmak, A.; Trezza, R.; Al, E. Satellite-based ET estimation in agriculture using SEBAL and METRIC. Hydrol. Process. 2011, 25, 4011–4027. [Google Scholar] [CrossRef]
  30. Laipelt, L.; Ruhoff, A.L.; Fleischmann, A.S. Assessment of an automated calibration of the SEBAL algorithm to estimate dry-season surface-energy partitioning in a forest–savanna transition in Brazil. Remote Sens. 2020, 12, 1108. [Google Scholar] [CrossRef]
  31. Long, D.; Singh, V.P.; Li, Z.L. How sensitive is SEBAL to changes in input variables, domain size and satellite sensor? J. Geophys. Res. 2011, 116, D21. [Google Scholar] [CrossRef]
  32. Wei, G.; Cao, J.; Xie, H.; Xie, H.; Yang, Y.; Wu, C.; Cui, Y.; Luo, Y. Spatial-temporal variation in paddy evapotranspiration in subtropical climate regions based on the SEBAL model: A case study of the Ganfu Plain irrigation system, southern China. Remote Sens. 2022, 14, 1201. [Google Scholar] [CrossRef]
  33. Huang, D.; Wang, J.; Khayatnezhad, M. Estimation of actual evapotranspiration using soil moisture balance and remote sensing. Iran. J. Sci. Technol. Trans. Civ. Eng. 2021, 45, 2779–2786. [Google Scholar] [CrossRef]
  34. Sahnoun, F.; Abderrahmane, H.; Kaddour, M.; Abdelkader, K.; Mohamed, B.; Castro, T.A.H.D. Application of SEBAL and Ts/VI trapezoid models for estimating actual evapotranspiration in the Algerian semi-arid environment to improve agricultural water management. Rev. Bras. De Meteorol. 2021, 36, 219–236. [Google Scholar] [CrossRef]
  35. Shuttleworth, W.J.; Gurney, R.J.; Hsu, A.Y. FIFE The variation in energy partition at surface flux sites. IAHS Publ. 1989, 186, 523–534. [Google Scholar]
  36. Shuttleworth, W.J. Macrohydrology—The new challenge for process hydrology. J. Hydrol. 1988, 100, 31–56. [Google Scholar] [CrossRef]
  37. Shuttleworth, W.J.; Gurney, R.J. The theoretical relationship between foliage temperature and canopy resistance in sparse crops. Q. J. R. Meteorol. Soc. 1990, 116, 497–519. [Google Scholar] [CrossRef]
  38. Ke, Y.; Im, J.; Park, S.; Al, E. Downscaling of MODIS one kilometer evapotranspiration using Landsat-8 data and machine learning approaches. Remote Sens. 2016, 8, 215. [Google Scholar] [CrossRef]
  39. Muñoz-Sabater, J.; Dutra, E.; Agustí-Panareda, A.; Al, E. ERA5-Land: A state-of-the-art global reanalysis dataset for land applications. Earth Syst. Sci. Data 2021, 13, 4349–4383. [Google Scholar] [CrossRef]
  40. Hersbach, H.; Bell, B.; Berrisford, P.; Al, E. The ERA5 global reanalysis. Q. J. R. Meteorol. Soc. 2020, 146, 1999–2049. [Google Scholar] [CrossRef]
  41. Shuttleworth, W.J. Terrestrial Hydrometeorology; John Wiley and Sons: West Sussex, UK, 2012. [Google Scholar]
  42. Laipelt, L.; Kayser, R.H.B.; Fleischmann, A.S. Long-term monitoring of evapotranspiration using the SEBAL algorithm and Google Earth Engine cloud computing. Isprs J. Photogramm. Remote Sens. 2021, 178, 81–96. [Google Scholar] [CrossRef]
  43. Allen, R.G.; Pereira, L.S.; Raes, D.; Al, E. FAO Irrigation and drainage paper No. 56. Rome Food Agric. Organ. United Nations 1998, 56, e156. [Google Scholar]
  44. Ren, G.; Guo, J. Change in pan evaporation and the influential factors over China 1956-2000. J. Nat. Resour. 2006, 21, 31–44. [Google Scholar]
  45. GB/T 35230-2017; Specification for Surface Meteorological Observations—Evaporation. China Meteorological Administration: Beijing, China, 2017.
  46. Yang, J.; Huang, X. 30 m annual land cover and its dynamics in China from 1990 to 2019. Earth Syst. Sci. Data Discuss. 2021, 2021, 1–29. [Google Scholar]
  47. Tasumi, M.; Allen, R.G.; Trezza, R. At-surface reflectance and albedo from satellite for operational calculation of land surface energy balance. J. Hydrol. Eng. 2008, 13, 51–63. [Google Scholar] [CrossRef]
  48. IFS Documentation CY41R2; IFS Documentation CY41R2—Part VII: ECMWF Wave Model. European Centre for Medium-Range Weather Forecasts: Reading, UK, 2016.
  49. Buck, A.L. New equations for computing vapor pressure and enhancement factor. J. Appl. Meteorol. 1981, 20, 1527–1532. [Google Scholar] [CrossRef]
  50. Alduchov, O.A.; Eskridge, R.E. Improved Magnus form approximation of saturation vapor pressure. J. Appl. Meteorol. Clim. 1996, 35, 601–609. [Google Scholar] [CrossRef]
  51. Allen, R.G.; Burnett, B.; Kramber, W.; Al, E. Automated calibration of the metric-landsat evapotranspiration process. Jawra J. Am. Water Resour. Assoc. 2013, 49, 563–576. [Google Scholar] [CrossRef]
  52. Sung, J.H.; Chung, E.S.; Kim, Y.; Al, E. Meteorological hazard assessment based on trends and abrupt changes in rainfall characteristics on the Korean peninsula. Theor. Appl. Clim. 2017, 127, 305–326. [Google Scholar] [CrossRef]
  53. Hao, X.; Chen, Y.; Xu, C. Impacts of climate change and human activities on the surface runoff in the Tarim River Basin over the last fifty years. Water Resour. Manag. 2008, 22, 1159–1171. [Google Scholar] [CrossRef]
  54. Jin, X.; Zhu, X.; Xue, Y. Satellite-based analysis of regional evapotranspiration trends in a semi-arid area. Int. J. Remote Sens. 2019, 40, 3267–3288. [Google Scholar] [CrossRef]
  55. Yang, L.; Feng, Q.; Zhu, M.; Wang, L.; Alizadeh, M.R.; Adamowski, J.F.; Wen, X.; Yin, Z. Variation in actual evapotranspiration and its ties to climate change and vegetation dynamics in northwest China. J. Hydrol. 2022, 607, 127533. [Google Scholar] [CrossRef]
  56. Ning, T.; Li, Z.; Feng, Q.; Li, Z.; Qin, Y. Attribution of growing season evapotranspiration variability considering snowmelt and vegetation changes in the arid alpine basins. Hydrol. Earth Syst. Sci. 2021, 25, 3455–3469. [Google Scholar] [CrossRef]
  57. Wang, Y.; Liu, B.; Su, B.; Zhai, J.; Gemmer, M. Trends of Calculated and Simulated Actual Evaporation in the Yangtze River Basin. J Clim. 2011, 24, 4494–4507. [Google Scholar] [CrossRef]
  58. Qiu, G.Y. Water and Energy: Evapotranspiration, Thermal Environment, and Energy Balance; Science Press: Beijing, China, 2014. [Google Scholar]
  59. Tang, R.; Li, Z.L.; Chen, K.S.; Jia, Y.; Li, C.; Sun, X. Spatial-scale effect on the SEBAL model for evapotranspiration estimation using remote sensing data. Agric. For. Meteorol. 2013, 174, 28–42. [Google Scholar] [CrossRef]
  60. Wang, Z.Y.; Kong, J.L.; Li, Y.; Zhang, Z.Y.; Liu, H.H.; Jiang, Y.Z.; Zhong, Y.L.; Zhang, J.Y. An analysis of spatio-temporal characteristics and influencing factors of surface evapotranspiration in the Yinchuan Plain based on MOD16 data. Hydrogeol. Eng. Geol. 2021, 48, 53–61. [Google Scholar]
Sustainability 16 08025 g001
Figure 1. Location of the meteorological stations and Landsat images that were used to illustrate the land cover conditions in the Yinchuan Plain area.
Figure 1. Location of the meteorological stations and Landsat images that were used to illustrate the land cover conditions in the Yinchuan Plain area.
Sustainability 16 08025 g001
Sustainability 16 08025 g002
Figure 2. Comparison between meteorological stations with large ET and geeSEBAL ET. Diamond-shaped points represent outliers lying outside the 150% inter-quartile range.
Figure 2. Comparison between meteorological stations with large ET and geeSEBAL ET. Diamond-shaped points represent outliers lying outside the 150% inter-quartile range.
Sustainability 16 08025 g002
Sustainability 16 08025 g005
Figure 5. Seasonal ETa changes in the Yinchuan Plain. (a) Spring ET distribution; (b) Summer ET distribution; (c) Autum ET distribution; (d) Winter ET distribution.
Figure 5. Seasonal ETa changes in the Yinchuan Plain. (a) Spring ET distribution; (b) Summer ET distribution; (c) Autum ET distribution; (d) Winter ET distribution.
Sustainability 16 08025 g005
Sustainability 16 08025 g006
Figure 6. Trends in ETa on the Yinchuan Plain from 1987 to 2020.
Figure 6. Trends in ETa on the Yinchuan Plain from 1987 to 2020.
Sustainability 16 08025 g006
Sustainability 16 08025 g007
Figure 7. Area of Yinchuan Plain land use types.
Figure 7. Area of Yinchuan Plain land use types.
Sustainability 16 08025 g007
Sustainability 16 08025 g008
Figure 8. ETa in different subsurface types.
Figure 8. ETa in different subsurface types.
Sustainability 16 08025 g008
Sustainability 16 08025 g009
Figure 9. Comparison of remote sensing imagery (a), land use classification (b), and ET imagery (c) on the Yinchuan Plain.
Figure 9. Comparison of remote sensing imagery (a), land use classification (b), and ET imagery (c) on the Yinchuan Plain.
Sustainability 16 08025 g009
Sustainability 16 08025 g010
Figure 10. Impervious areas misclassified in some intersecting land types. water bodies are identified as impervious areas (red color).
Figure 10. Impervious areas misclassified in some intersecting land types. water bodies are identified as impervious areas (red color).
Sustainability 16 08025 g010
Sustainability 16 08025 g011
Figure 11. ETa contribution of different subsurface types.
Figure 11. ETa contribution of different subsurface types.
Sustainability 16 08025 g011
Sustainability 16 08025 g012
Figure 12. geeSEBAL with the batch image estimation mode.
Figure 12. geeSEBAL with the batch image estimation mode.
Sustainability 16 08025 g012
Table 1. Ningxia Meteorological Bureau Observed Data. Here, 53518, 53519, 53614 and 53615 include a large evaporator (L) and small evaporator (S), both of which are 24 h ET values.
Table 1. Ningxia Meteorological Bureau Observed Data. Here, 53518, 53519, 53614 and 53615 include a large evaporator (L) and small evaporator (S), both of which are 24 h ET values.
Meteorological StationsLatitudeLongitudeAltitude(m)Measured InstrumentsMeasured Time Series
535183904106351126.9L1991–2018
S1991–2014
535193921106761092.2L1991–2018
S1991–2014
536103856106341106.0S1991–2014
536143847106211111.6L1991–2018
S1991–2014
536153881106701101.0L1991–2018
S1991–2014
536183822106201118.4S1991–2014
536193812106301115.7S1991–2014
Table 2. Seasonal images and end member criterion.
Table 2. Seasonal images and end member criterion.
SeasonImage DataNDVIcoldTscoldNDVIhotTshot
SpringLC08_129034_202003065%20%10%10%
SummerLC08_129033_202008135%20%10%20%
AutumLC08_129033_202010165%20%10%20%
WinterLC08_129033_201401175%10%1%10%
Table 3. Remote sensing images and end member criterion.
Table 3. Remote sensing images and end member criterion.
Remote Sensing ImageNDVIcoldTscoldNDVIhotTshot
LT05_129033_198709205%20%10%20%
LT05_129033_198908245%20%10%20%
LT05_129033_199108305%20%10%1%
LT05_129033_199208165%20%10%20%
LT05_129033_199306165%20%10%20%
LT05_129033_199408225%20%10%20%
LT05_129033_199608115%20%10%1%
LT05_129033_199908205%20%10%1%
LT05_129033_200108095%20%10%1%
LT05_129033_200206255%20%10%1%
LT05_129033_200609085%20%10%1%
LT05_129033_200707095%10%1%10%
LC08_129034_201407285%20%10%20%
LC08_129033_201709065%20%10%20%
LC08_129033_201808245%10%1%10%
LC08_129034_201908115%20%10%20%
LC08_129033_202008135%10%1%10%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Meng, J.; Yang, X.; Li, Z.; Zhao, G.; He, P.; Xuan, Y.; Wang, Y. Tracking Evapotranspiration Patterns on the Yinchuan Plain with Multispectral Remote Sensing. Sustainability 2024, 16, 8025. https://doi.org/10.3390/su16188025

AMA Style

Meng J, Yang X, Li Z, Zhao G, He P, Xuan Y, Wang Y. Tracking Evapotranspiration Patterns on the Yinchuan Plain with Multispectral Remote Sensing. Sustainability. 2024; 16(18):8025. https://doi.org/10.3390/su16188025

Chicago/Turabian Style

Meng, Junzhen, Xiaoquan Yang, Zhiping Li, Guizhang Zhao, Peipei He, Yabing Xuan, and Yunfei Wang. 2024. "Tracking Evapotranspiration Patterns on the Yinchuan Plain with Multispectral Remote Sensing" Sustainability 16, no. 18: 8025. https://doi.org/10.3390/su16188025

APA Style

Meng, J., Yang, X., Li, Z., Zhao, G., He, P., Xuan, Y., & Wang, Y. (2024). Tracking Evapotranspiration Patterns on the Yinchuan Plain with Multispectral Remote Sensing. Sustainability, 16(18), 8025. https://doi.org/10.3390/su16188025

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