Daily actual evapotranspiration estimation of different land use types based on SEBAL model in the agro-pastoral ecotone of northwest China.

Evapotranspiration (ET) plays a crucial role in hydrological and energy cycles, as well as in the assessments of water resources and irrigation demands. On a regional scale, particularly in the agro-pastoral ecotone, clarification of the distribution of surface ET and its influencing factors is critical for the rational use of water resources, restoration of the ecological environment, and protection of ecological water sources. The SEBAL model was used to invert the regional ET based on Landsat8 images in the agro-pastoral ecotone of northwest China. The results were indirectly verified by monitoring data from meteorological stations. The correlation between ET and surface parameters was analyzed. Thus, the main factors that affect the surface ET were identified. The results show that the SEBAL model determines an accurate inversion, with a correlation coefficient of 0.81 and an average root mean square error of 0.9 mm/d, which is highly suitable for research on water resources. The correlation coefficients of normalized vegetation index, surface temperature, land surface albedo, net radiation flux with daily ET were 0.5830, 0.8425, 0.3428 and 0.9111, respectively. The normalized vegetation index and the net radiation flux positively correlated with the daily ET, while the surface temperature and land surface albedo negatively correlated with the daily ET. The correlation from strong to weak is the net radiation flux > surface temperature > normalized vegetation index > surface albedo. In terms of spatial distribution, the daily ET of water was the highest, followed by woodland, wetland, cropland, built-up land, shrub land, grassland and bare land. However, the SEBAL model overestimates the inversion of daily ET of built-up land.

Introduction

Evapotranspiration (ET) is an important part of the global water cycle [1]. It is the channel that transforms surface and atmospheric water, which directly affects the spatial distribution of global precipitation and vegetation [2,3]. Studies have shown that > 60% of the precipitation that reaches the surface through condensation returns to the atmosphere through the process of ET, and this proportion can be as high as 90% owing to a reduction in precipitation and a higher ET in arid areas [4-6]. Solar short-wave radiation energy returns to the atmosphere in the form of latent heat during the process of ET on the surface, thus, forming an energy cycle process. Therefore, the surface ET plays an important role in the global water and energy cycles [7]. Simultaneously, the actual surface ET is a key parameter for research on environmental changes [8], hydrological processes [9], drought trends [10,11], and the potential evaluation of agricultural development [12,13]. Therefore, the study of surface ET is highly important scientifically to ensure the rational use of water resources and the protection of ecosystems in the study area [14,15].

Owing to the complexity of the surface, current ET observation methods include lysimeters [16], eddy correlator large aperture scintillation flux meters [17] and other methods. Such ET monitoring methods can obtain the ET more accurately at a specific location, but simultaneously, it is time-consuming, labor-intensive and expensive, and only a single point of ET data can be obtained. Therefore, the acquisition of regional surface ET data that is highly accurate over a continuous time-scale is always the key to solving difficult issues of research in the fields of meteorological science, hydrology, agriculture and geography. The development of remote sensing technology provides a new platform and way to obtain the actual surface ET, and it solves the problem of a lack of measured data owing to geographical reasons. In recent years, researchers have made substantial progress in utilizing remote sensing ET models for research, which successively established the surface energy balance index (SEBI) model [18], simplified surface energy balance index (S-SEBI) model [19], surface energy balance algorithm for land (SEBAL) model [20] and surface energy balance system (SEBS) model [21]. Among them, the SEBAL model benefits from a clearer mechanism, fewer model input parameters, simple data acquisition, high inversion accuracy, and a wide application range. It has become one of the most commonly used remote sensing methods to determine the inversion of surface ET [22].

However, the accuracy of SEBAL model varies in the estimation of ET for different products and environmental conditions. Liu and Hu investigated the effects of land use/cover change and climate change on wetland ET using SEBAL, and the results showed that the average relative error of regional ET estimated by the SEBAL model was 9.01% [23]. Rahimzadegan et al. studied the efficiency of SEBAL algorithm in estimating the ET of pistachio crops. The results of this study indicated that the determination coefficient and Root Mean Square Error of the estimated actual ET for pistachio plants was 0.8 and 2.5 mm, respectively [24]. The agro-pastoral ecotone of northwest China receives little rainfall, which is distributed unevenly, and has a strong surface ET and an extremely fragile ecological environment. A shortage of water is the most serious problem in this area. In addition, the research on the influence factors of ET mostly focuses on meteorological factors [25,26], and there is a lack of in-depth research on the influence of surface parameters on ET.

The objective of this study was to explore the applicability of SEBAL model and analyze the influencing factors in a wind and sandy grass beach area of northwest China. Therefore, in this study, (1) Landsat8 imagery and meteorological data were used as the source of data to estimate daily ET based on the SEBAL model; (2) the accuracy of estimation of daily ET was verified and used to analyze its spatiotemporal distribution, and (3) the influence of surface parameters on the daily ET was studied, providing an analysis of the response of daily ET to different land use types. These conclusions will enhance our understanding of the factors that affect ET and thus, provide a basis for the sustainable development of water resources and the restoration of the ecological environment in arid and semiarid regions in northwest China.

Materials and methods

Study area

The study area (107°58′-109°15′E, 37°31′-38°49′) is located at the junction of Yulin City and Ordos City, China (Fig 1) in a sandy and grassy beach area southeast of Mu Us Sandy Land. It has a typical farming-pastoral transition zone with a composite underlying surface of bare land-grassland-farmland in China. The study area has a typical continental monsoon climate. The northwest is arid and lacks water, which results in infrequent precipitation. Consequently, the ecological environment is fragile. The land use types are primarily grassland and bare land. The southeast is low-lying, with rivers that pass through, and the land use types are primarily grassland and cropland. The average annual temperature is approximately 7°C and the average annual precipitation is between 350–450 mm, which decreases from southeast to northwest in the study area. The interannual precipitation is highly uneven with a large seasonal variation. The rainy season is from July to September, and it comprises more than 50% of the rain that falls during a whole year. The potential ET is relatively large, and the annual average ET is 2,000–3,000 mm in the study area. The sparse precipitation and strong ET activities are the main reasons for the drought in study area.

Fig 1

Location of the study area in agro-pastoral ecotone of northwest China.

Data sources

In this study, the agro-pastoral ecotone of northwest China was used as the study area, and the Landsat8 images that cover the study area were numbered 128,033 and 128,034. They are all freely available at the Geospatial Data Cloud (http://www.gscloud.cn). Since the data is affected by cloud coverage, and parts of it are missing, it is impossible to obtain the image data for each period of the same year. To ensure that the data images can cover the study area, this study screened out 2017-1-2, 2016-4-21, 2017-5-26, 2017-7-13, 2016-9-28, 2016-12-17 and six other Landsat8 original remote sensing images. The selected data covered the four seasons of the year and are suitable for the study of the spatiotemporal distribution of surface ET. The Landsat8 remote sensing data requires preprocessing, such as radiometric calibration, atmospheric correction and image mosaics, before it can be substituted into the model. ENVI software (L3Harris Geospatial, Boulder, CA, USA) was used to process the images in this research.

Principle of the SEBAL model

In 1995, Bastiaanssen proposed the SEBAL inversion model, which was verified and optimized in 1998. In recent years, the surface ET model has been effectively developed and improved. The SEBAL model has been utilized for its advantages and widely used throughout the world. It has produced good inversion results and has become a high-precision regional ET remote sensing inversion model. The research foundation of the SEBAL model is the energy balance equation, which is composed of four parts, including net radiance, soil heat, sensible heat and latent heat fluxes, which form a closed-loop energy transmission. The formula is shown as follows: Where R is the net solar radiation flux (W/m2); G is the soil heat flux (W/m2); H is the turbulent sensible heat flux (W/m2); λET is the turbulent latent heat flux (W/m2), and λ is the coefficient of latent heat of vaporization of water, which is a function of temperature.

Net radiant flux (R)

The net surface radiant flux is expressed per unit area, which includes the difference between total energy received from the solar short-wave radiation and atmospheric long-wave radiation on the surface and the long-wave radiation emitted by the surface, reflected solar short-wave radiation and atmospheric long-wave radiation [27]. The process of calculation is as follows: Where R is the short-wave radiation energy that the sun reaches the surface (W/m2); R is the long-wave radiation from the atmosphere to the surface (W/m2); R is the ground upgoing long-wave radiation (W/m2); α is the surface albedo; ε is the surface emissivity; G is the solar constant, with a value of 1,367 W/m2; θ is the zenith angle of sun; d is the distance factor between the sun and the earth; τ is the atmospheric transmittance; DAY is the number of days of the image in the current year; ε is the specific emissivity of the atmosphere; σ is a constant with a value of 5.67×10−8 W/(m2·K4); T represents the temperature at the ground reference altitude, and LST represents the surface temperature.

Soil heat flux (Gn)

The soil heat flux is the difference between surface and soil temperatures that enables the transmission of energy. A small part of the energy lost in the surface soil during this process is primarily affected by topography and land use types. The soil heat flux generally occupies a relatively small proportion in the energy balance equation. Factoring in the greater difficulty in its direct calculation, Bastiaanssen used calculations of net radiant flux, vegetation index, land surface temperature and land surface reflectance in the SEBAL model [28], and the calculations are shown as follows: Where NDVI is the normalized vegetation index.

Sensible heat flux (H)

Sensible heat flux is the temperature difference between the air and surface temperatures, and the energy that flows into the air from the surface is in the form of vertical transmission. The sensible heat flux is primarily related to the ground temperature difference, density, aerodynamic impedance and other parameters, and it is calculated is as follows: Where C is the specific heat capacity of the atmosphere; ρ is the density of the atmosphere; r is the aerodynamic impedance (s·m-1), and d is the ground temperature difference (K). r and d are unknown parameters. The specific process for calculation has been previously described [29].

Latent heat flux (λET)

Latent heat flux is the heat exchange per unit area at constant temperature, which is the most important parameter in the inversion study of surface ET. It is the energy absorbed during the vaporization of water. Based on the surface energy balance formula in the remote sensing inversion ET model, the latent heat flux is obtained as follows: λET represents the instantaneous latent heat flux (W/m2).

Daily ET

The surface radiation flux is substantially affected by factors, such as surface temperature, wind speed, and surface coverage, and cannot be obtained directly. The difference in ratio of latent heat flux to the net radiant and soil heat fluxes in a day always remain basically stable, i.e., the theory of constant evaporation ratio and the instantaneous ET can be calculated by the theory of constant evaporation ratio, which can be extended to the daily ET [30,31]. The formula is calculated as follows: Where EF represents the instantaneous evaporation ratio; EF24 represents the entire day evaporation ratio; R represents the daily net radiant flux (W/m2), and G represents the daily soil heat flux (W/m2).

Results

Verification of the accuracy of daily ET

The verification and evaluation of the inversion accuracy of the surface ET model is the focus and difficulty of the research on ET. To verify the reliability of the daily ET retrieved by the SEBAL model, the water surface evaporation monitored by four weather stations in the study area was selected for relative verification. However, the evaporation on the water surface that was measured was significantly higher than the actual evaporation on the surface, and the water surface evaporator of the weather station has two types that include E-601 and 20cm-caliber. There are differences in the installation plan, structure and observational methods of the two instruments. There are large deviations in the verification of data when the data of weather stations is used directly. Therefore, to more accurately explore the inversion accuracy of the SEBAL model and maintain the continuity and consistency of the data at each station, it is necessary to explore the conversion coefficients of the two surface evaporators and unify the data. E-601 and 20cm-caliber meteorological stations, such as those at the Hengshan, Jingbian, Shenmu, Dingbian, Yanchi and Etuoke Banner Stations in China were compared in the study and surrounding areas. The monitoring data were analyzed and indicated that the average conversion coefficient of the two surface evaporators was 0.63, and the determination coefficient R2 was 0.95, indicating that the conversion coefficients of the E-601 type and the 20cm-caliber small evaporator are highly reliable, respectively, and can meet the requirements of the two types of surface evaporators.

The monitoring data of 20cm-caliber small evaporator was converted into the E-601 monitoring data using the conversion coefficients of the two surface evaporators. This enabled the use of data to verify the SEBAL model inversion result. Since the meteorological sites were small meteorological stations that were established by research teams in Wushenqi and Henan in China, they did not have evaporation dishes to obtain the data on water surface evaporation. The only data involved in the validation were the monitoring data from the meteorological stations in Hengshan and Jingbian, China. Fig 2 is a scatter diagram of the water surface evaporation of the meteorological station and the inversion results. This figure indicates that the actual surface ET based on the SEBAL model is highly correlated with the surface water surface ET. The coefficient of determination was 0.81, and the average root mean square error was 0.9 mm/d. These values indicate that the SEBAL model can be used to retrieve the surface ET in agro-pastoral ecotone of northwest China with a high degree of accuracy.

Fig 2

Validation scatter plot of the SEBAL model inversion results and meteorological data.

Spatiotemporal distribution of the daily ET

In this study, the daily ET of the study area was obtained based on the theory of constant ET ratio. Table 1 shows the data on daily ET in different periods in the study area. An examination of the distribution revealed a pattern in the daily ET in the northwest agro-pastoral ecotone area. The pattern included July (5.17 mm/d)>May (3.51 mm/d)>September (3.31 mm/d)>April (2.02 mm/d)>December (0.72 mm/d)>January (0.56 mm/d). The possible reason is related to the local climate change during the four seasons. The temperature and precipitation were the lowest during the winter in January and December, and insufficient hydrothermal conditions limited the ET process. In April and May, the air temperature began to increase; the vegetation gradually turned green, and the crops were just beginning to be cultivated and sprouting. These factors indicate that the ET gradually increased. Since July is in the summer in this study area, the monthly average temperature reaches the highest value of the year, and this is the rainy season. The vegetation grows vigorously, and there is sufficient soil moisture. The hydrothermal conditions were the highest of the year. Therefore, the study area had its highest daily ET in the summer. In September, the vegetation in the sand of study area began to senesce; the temperature decreased; the rainfall was significantly reduced; the hydrothermal conditions were all lower than those in the summer, and the daily ET was less than in the summer. In summary, hydrothermal conditions are the main factor that affect the changes in evapotranspiration during the year.

Table 1

Daily average evapotranspiration in the study area.

Date	Min. (mm/d)	Max. (mm/d)	Daily average ET (mm/d)	SD
2017-01-02	0.02	3.44	0.56	0.3477
2016-04-21	0.15	12.58	2.02	1.4135
2017-05-26	0.35	15.13	3.51	2.1273
2017-07-13	0.51	14.56	5.17	2.0476
2016-09-28	0.22	8.07	3.31	1.1839
2016-12-17	0.07	3.36	0.72	0.3363

ET, evapotranspiration; Min., minimum; Max, maximum; SD, standard deviation.

Fig 3 shows the spatial distribution of the daily ET in different periods in the study area. The evapotranspiration varied with regions within the study area. It was higher overall in the southeast and lower in the northwest, and such a pattern of distribution corresponds to the terrain of the study area. There are rivers that pass through the southeast, and there are many areas of crops, woodland and grassland. The type of land use in northwest China primarily included denuded dunes or bare land. Thus, there is less vegetation coverage, and the soil is dry and lacks water. From a local perspective, the daily ET of different land use types differs significantly.

Fig 3

Spatiotemporal distribution of daily evapotranspiration in the agro-pastoral ecotone of northwest China.

Discussion

A correlation analysis of ET and meteorological factors

In this study, the annual daily data of the Hengshan station in 2016 were selected, including four meteorological factors: ET (mm/day), air temperature (TEM, °C), relative humidity (RHU, %), and wind speed (WIN, m/s). A correlation analysis between meteorological factors and the surface ET was performed using SPSS, and the strength of correlation between the meteorological factors and ET was evaluated using the correlation coefficient R2. Fig 4 shows the linear relationship between ET and the meteorological factors. ET significantly correlated with the air temperature with an R2 of 0.6274. The wind speed weakly correlated with an R2 of 0.3307. The RHU had a negative R2 of 0.3374. In areas with higher air humidity, the surface moisture faces more resistance when entering the air, so the ET was relatively small.

Fig 4

The linear relationship between evapotranspiration and meteorological factors.

Sensitivity analysis of model parameters

To better understand the extent to which each model parameter affects the model inversion results, this study performed a sensitivity analysis of the model parameters to derive the amount of change in the final ET owing to changes in each parameter. The sensitivity analysis will be performed using the one-parameter variable method, i.e., the interaction of the parameters will not be considered. Only the parameters to be analyzed will be treated for variation. Since hydrothermal conditions have a large influence on ET, this study explored the sensitivity of atmospheric water vapor (wv) content, hot spot surface temperature (LSTh) and cold pixels to ET. The sensitivity index I was used to reflect the sensitivity of the evapotranspiration model to water vapor content. Larger values of I indicated that the model is highly sensitivity to this parameter. The formula is calculated as follows [32]: Where ΔX represents the amount of change in the parameter, and ΔET represents the amount of change in the model output results owing to the change in the parameter. The sensitivity indices were classified into four classes by the magnitude of the value of I as shown in Table 2.

Table 2

Sensitivity classification.

Grade	Sensitivity index	Sensitivity classification
I	0.0<|I|<0.05	poor
II	0.05<|I|<0.20	fair
III	0.20<|I|<1.00	moderate
IV	|I|>1.00	excellent

The inversion accuracy of wv and LSTh were considered to be approximately 0.5 g/cm2 and 2°C, respectively. Δwv = 0.5 g/cm2 and ΔLST = 2°C were established. Tables 3 and 4 show the mean values of the model outputs, as well as the sensitivity indices.

Table 3

Evapotranspiration changes and sensitivity indices from the changes in water vapor.

Daily ET	Average value (mm/d)	Average difference (mm/d)	|I|	Sensitivity classification
wv
wv+1.5 g/cm2	2.73	-2.44	0.83	moderate
wv+1.0 g/cm2	3.24	-1.93	0.98	moderate
wv+0.5 g/cm2	4.04	-1.13	1.15	excellent
wv-0.5 g/cm2	6.05	0.88	0.90	moderate
wv-1.0 g/cm2	6.73	1.56	0.79	moderate
wv-1.5g/cm2	7.24	2.07	0.70	moderate

ET, evapotranspiration; wv, atmospheric water vapor.

Table 4

Changes in evapotranspiration and sensitivity indices from the changes in hot surface temperature.

Daily ET	Average value (mm/d)	Average difference (mm/d)	|I|	Sensitivity classification
LSTh
LSTh+6°C	5.54	0.37	0.56	moderate
LSTh+4°C	5.42	0.25	0.57	moderate
LSTh+2°C	5.30	0.13	0.59	moderate
LSTh-2°C	5.09	-0.08	0.36	moderate
LSTh-4°C	5.05	-0.13	0.30	moderate
LSTh-6°C	4.97	-0.20	0.30	moderate

ET, evapotranspiration; LSTh, hot surface temperature.

Table 3 shows the daily changes in ET and the sensitivity indices owing to changes in the wv. Increasing the wv by 0.5 g/cm2 caused a daily change of 1.13 mm/d in the ET, and a decrease of 0.5 g/cm2 wv caused a daily change of 0.88 mm/d in the ET. The sensitivity indices for increasing and decreasing the wv by Δ, 2Δ, and 3Δ were all > 0.2, indicating that the output of the SEBAL model-based inversion of daily ET is sensitive to the wv. The main reason is that the magnitude of wv directly affects the atmospheric transmittance, which causes a change in the net radiation flux. The effect on daily ET is weakened by the large changes in wv, which is owing to the large number of parameters in the SEBAL model. The changes in wv affect other parameters, thus, weakening the effect on daily ET.

Table 4 shows the daily changes in the ET and sensitivity indices owing to changes in the LSTh. A 2°C increase in LSTh caused a daily change of 0.15 mm/d in the ET. A 2°C decrease in LSTh caused an ET change of 0.08 mm/d. The sensitivity indices for increasing and decreasing the LSTh by Δ, 2Δ, and 3Δ were all > 0.2, indicating that the inversion results are sensitive to the LSTh based on the SEBAL model inversion of daily ET. The model output results on the LSTh are less sensitive compared with those of the wv. The main reason for this is that the LSTh is used to calculate the sensible heat flux, and the result of the sensible heat flux requires calculations that are continuously iterative. Multiple iterations weaken the influence of the LSTh. Thus, the change of a single parameter of the LSTh does not have a large influence on the daily ET.

Correlation analysis of the daily ET and surface parameters

An analysis of the spatiotemporal changes of ET in the study area indicated that the underlying surface characteristics, such as hydrothermal conditions and land use types, are important factors that affect energy transmission, and their interactions in concert affect the spatial distribution pattern of the surface ET [22]. To analyze the correlation between surface parameters and ET in more detail, this study obtained the spatial distribution map of the NDVI, LST, surface albedo and net radiant flux (Fig 5) and their correlation with the daily ET (Fig 6). The figures indicate terms of spatial distribution. For example, the high-value areas of NDVI are primarily concentrated in the southern part of the study area and distributed in strips and points. This pattern is consistent with the spatial distribution of daily ET. The spatial distribution of surface temperature is opposite to that of the daily ET. The high values are distributed in the west and north of the study area on primarily bare and sandy lands. The high value areas of surface albedo are primarily distributed in the west and north, similar to the distribution of surface temperature. The main reason for this is that crops are grown in the south and east, and the northwest is primarily bare and sandy lands. The surface is not very rough, and there is strong reflection from solar energy. The pattern of spatial distribution of the net radiation flux and daily evaporation remains highly consistent. The correlation coefficients between the NDVI, LST, surface albedo, net radiant flux and daily ET were 0.5830, 0.8425, 0.3428, and 0.9111, respectively. Among them, the NDVI and net radiant flux positively correlated with the daily ET, and the LST and surface albedo negatively correlated with the daily ET. Overall, the correlation between daily ET and the underlying surface parameters from strong to weak is the net surface radiation > LST > NDVI > surface albedo. The correlation of surface parameters and daily ET as shown in this study has also been widely reported from many studies [6,33,34].

Fig 5

Spatial distribution of the normalized vegetation index, surface temperature, surface albedo and the net radiation flux in the agro-pastoral ecotone of northwest China.

Fig 6

Correlation between the daily evapotranspiration and surface parameters in the agro-pastoral ecotone of northwest China.

Analysis of the daily ET with different land use types

Fig 7 shows the spatial distribution of land use types and daily ET in the agro-pastoral ecotone of northwest China. The amount of evapotranspiration varies from one land use type to another owing to differences in the physicochemical properties of the subsurface. The spatial distribution of the larger areas of ET was significantly similar to that of water and cropland. The spatial distribution of the smaller areas matches that of the bare land, and the daily ET of the built-up land did not differ significantly from that of the other surrounding land use types. Table 5 shows the data of daily ET values for different land use types. The daily ET of water and woodland were the highest in study area, with average daily ET values of 9.90 mm/d and 8.69 mm/d, respectively, which is because water provides sufficient water resources. Woodland have a dual role of water connotation and ET and can provide good moisture conditions for ET, which should be relatively high. This result is consistent with a study by Jeremiah et al. [35]. The average daily ET of wetlands and croplands was 7.70 mm/d and 7.59 mm/d, respectively, which are relatively high categories of evapotranspiration. The combination of plant emissions in the results and evaporation from the soil and water surfaces results in high evapotranspiration. The croplands are predominantly dryland, and crop irrigation provides sufficient water for ET during warmer temperatures. The daily ET of built-up land, shrub land, grassland, and bare land were 6.39 mm/d, 5.74 mm/d, 4.92 mm/d and 3.78 mm/d, respectively. Among them, shrub lands are dominated by plantation woodland that are young and have a low density and few species of trees. Therefore, its ET is relatively low. Grassland and bare land comprise approximately 67.21% of the study area, which is the dominant land use type in the study area. The grasslands are primarily desert and sandy vegetation, which are sparsely distributed. Owing to the unique mosaic substratum characteristics of the region, grassland and bare lands are distributed in intervals. Transpiration is low. Therefore, it has a small value for its daily ET. There is only a small area of built-up land, and the surrounding land use types are mostly cropland. As a result, the daily ET data for the built-up land are higher overall than the daily ET of grassland and shrub land. In summary, water and woodland are the two land use types that had the highest daily ET, indicating that water resources and vegetation cover are the main factors that determined the daily ET in the study area. Wetlands have a higher ET owing to the dual effect of plant ET and evaporation from soil and water surfaces, and cropland has a higher ET owing to irrigation and other integrated human management measures. Shrub land, grassland, and bare land have a lower daily ET owing to the sparse vegetation and insufficient soil moisture. The daily ET of built-up land is clearly overestimated. The exact reasons for this merit further analysis and research.

Fig 7

Spatial distribution of (a) land use types and (b) daily evapotranspiration in agro-pastoral ecotone of northwest China.

Table 5

Daily evapotranspiration data for different land use types.

Land use types	Area/km2	Percentage of total area/%	Daily average ET (mm/d)
Cropland	1982.86	13.92	7.59
Woodland	9.40	0.07	8.69
Grassland	9571.70	67.21	4.92
Shrub land	312.56	2.19	5.74
Wetland	14.88	0.10	7.70
Water	55.89	0.39	9.90
Built-up land	126.57	0.89	6.39
Bare land	2167.58	15.22	3.78

ET, evapotranspiration.

The verification of the ET model is an important part of the remote sensing inversion of ET. The use of remote sensing to estimate the surface flux in areas < 1 km2 may have major problems, while the remote sensing inversion of large-scale areas is problematic owing to the problems of inadequate field observation data and errors in the observation instrument itself [14]. To more effectively vary the accuracy of the ET model, this study used the data of ET from the water surface of the weather station for relative verification and considered the problem of inconsistent water surface evaporators. The conversion coefficients of the E-601 type and 20cm-caliber small evaporator in the study area were explored using multiple sets of meteorological data, and the errors caused by the inconsistent evaporators were reduced. The conversion coefficients that were obtained are basically consistent with those in China provided by Mao Yue [36]. When the SEBAL model was used to retrieve the ET, the absolute error < 1 mm/d, indicating that the SEBAL model is highly applicable to retrieve the ET in agro-pastoral ecotone of northwest China. Similar results can also be found in other studies of different areas [6,37]. There is no doubt that the improvement in inversion accuracy of the SEBAL model merits celebration. However, surface parameters are important ones that affect the accuracy of the ET model. Norman et al. delineated the error possible from the ground temperature, i.e., with a speed of 5 m/s for a canopy that is 10 m high and an error of the ground temperature difference that exceeds 1 K, the error of H can reach 87 W/m2 [38]. A study by Li et al. [22] found that the parameters with the strongest correlation between ET and surface characteristic parameters were surface temperature and net radiant flux, which are consistent with the results of this study. The uncertainty of current radiation temperature inversion is usually 1~3 K [39]. While the remote sensing ET model is widely used, it is also necessary to realize that there is still substantial uncertainty in the inversion of the surface parameters in the remote sensing ET model. Therefore, the future research of remote sensing ET models should ensure that each step of the inversion result and each surface parameter can be verified, so that the inversion result is accurate. The analysis of daily ET of different land use types found that the daily ET of built-up land was higher than that of the grasslands and shrub land. This may be owing to two reasons. First, the drought and water shortage in the study area can be owing to the sparseness of grasslands and shrub lands and their lower daily ET. The built-up land is mostly surrounded by cropland, which results in inaccurate data on the daily ET of built-up land. Secondly, the SEBAL model overestimates the inversion of the ET of built-up land, which is consistent with the results of Jin et al. [40].

Conclusions

The SEBAL model was used to retrieve the daily surface ET based on Landsat8 remote sensing images in the northwestern windy sand and grassland area. This data was combined with that of the water surface evaporation from a weather station to verify the accuracy of the inversion. In addition, the correlation between the ET and surface parameters was analyzed to clarify the impact factor of surface ET. This study produced the following conclusions:

An analysis of the monitoring data of the E-601 type and 20cm-caliber small evaporator in the same period indicated that the average conversion coefficient of the two surface evaporators was 0.63, and the correlation coefficient R2 was 0.95, indicating that the conversion coefficients of E-601 type and the 20cm-caliber small evaporator are highly reliable and can meet the conversion needs of the monitoring data of the two surface evaporators.

The actual surface ET retrieved based on the SEBAL model correlated highly with the surface water ET. The correlation coefficient was 0.81, and the average root mean square error was 0.9 mm/d, indicating that the SEBAL model that is used to retrieve surface ET is highly accurate in the agro-pastoral ecotone of northwest China.

The correlation coefficients between NDVI, LST, surface albedo and net radiant flux and daily ET were 0.5830, 0.8425, 0.3428, and 0.9111, respectively. Among them, the NDVI and net radiant flux positively correlated with the daily ET. The LST and surface albedo negatively correlated with the daily ET. The overall correlation between the daily ET and underlying surface parameters from strong to weak were the net radiant flux > LST > NDVI > surface albedo.

The amount of daily ET varied substantially in terms of spatial distribution. The daily ET of water was the highest, followed by woodland, wetland, cropland, built-up land, and shrubland. The area of grassland and bare land were the dominant land use types in the study area, but they had the lowest value of daily ET. The inversion of daily ET of built-up land based on the SEBAL model was overestimated.

Shape files of image coverage of four meteorological station.

(ZIP)

Click here for additional data file.

ET inversion image based on SEBAL model.

(ZIP)

Click here for additional data file.

Land use types in the study area.

(ZIP)

Click here for additional data file.

29 Oct 2021

A rebuttal letter that responds to each point raised by the academic editor and reviewer(s). You should upload this letter as a separate file labeled 'Response to Reviewers'.

A marked-up copy of your manuscript that highlights changes made to the original version. You should upload this as a separate file labeled 'Revised Manuscript with Track Changes'.

An unmarked version of your revised paper without tracked changes. You should upload this as a separate file labeled 'Manuscript'.

If you would like to make changes to your financial disclosure, please include your updated statement in your cover letter. Guidelines for resubmitting your figure files are available below the reviewer comments at the end of this letter.

If applicable, we recommend that you deposit your laboratory protocols in protocols.io to enhance the reproducibility of your results. Protocols.io assigns your protocol its own identifier (DOI) so that it can be cited independently in the future. For instructions see: https://journals.plos.org/plosone/s/submission-guidelines#loc-laboratory-protocols. Additionally, PLOS ONE offers an option for publishing peer-reviewed Lab Protocol articles, which describe protocols hosted on protocols.io. Read more information on sharing protocols at https://plos.org/protocols?utm_medium=editorial-email&utm_source=authorletters&utm_campaign=protocols.

We look forward to receiving your revised manuscript.

Kind regards,

Débora Regina Roberti, Ph.D.

Academic Editor

PLOS ONE

Journal Requirements:

When submitting your revision, we need you to address these additional requirements.

1. Please ensure that your manuscript meets PLOS ONE's style requirements, including those for file naming. The PLOS ONE style templates can be found at

https://journals.plos.org/plosone/s/file?id=wjVg/PLOSOne_formatting_sample_main_body.pdf and

https://journals.plos.org/plosone/s/file?id=ba62/PLOSOne_formatting_sample_title_authors_affiliations.pdf

2. Thank you for stating the following financial disclosure:

This research was funded by Technology Innovation Center for Land Engineering and Human Settlements, Shaanxi Land Engineering Construction Group Co., Ltd and Xian Jiaotong University (2021WHZ0088, 2021WHZ0091) and Scientific Research Item of Shaanxi Provincial Land Engineering Built-up Group (DJNY2021-33).

Please state what role the funders took in the study.  If the funders had no role, please state: "The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript."

If this statement is not correct you must amend it as needed.

Please include this amended Role of Funder statement in your cover letter; we will change the online submission form on your behalf.

3. Thank you for stating the following in the Acknowledgments / Funding Section of your manuscript:

This research was funded by Technology Innovation Center for Land Engineering and Human Settlements, Shaanxi Land Engineering Construction Group Co., Ltd and Xian Jiaotong University (2021WHZ0088, 2021WHZ0091) and Scientific Research Item of Shaanxi Provincial Land Engineering Built-up Group (DJNY2021-33).

Please note that funding information should not appear in the Acknowledgments section or other areas of your manuscript. We will only publish funding information present in the Funding Statement section of the online submission form.

Please remove any funding-related text from the manuscript and let us know how you would like to update your Funding Statement. Currently, your Funding Statement reads as follows:

This research was funded by Technology Innovation Center for Land Engineering and Human Settlements, Shaanxi Land Engineering Construction Group Co., Ltd and Xian Jiaotong University (2021WHZ0088, 2021WHZ0091) and Scientific Research Item of Shaanxi Provincial Land Engineering Built-up Group (DJNY2021-33).

Please include your amended statements within your cover letter; we will change the online submission form on your behalf.

4. We note that Figure 1, 3, 4 and 6 in your submission contain [map/satellite] images which may be copyrighted. All PLOS content is published under the Creative Commons Attribution License (CC BY 4.0), which means that the manuscript, images, and Supporting Information files will be freely available online, and any third party is permitted to access, download, copy, distribute, and use these materials in any way, even commercially, with proper attribution. For these reasons, we cannot publish previously copyrighted maps or satellite images created using proprietary data, such as Google software (Google Maps, Street View, and Earth). For more information, see our copyright guidelines: http://journals.plos.org/plosone/s/licenses-and-copyright.

We require you to either (1) present written permission from the copyright holder to publish these figures specifically under the CC BY 4.0 license, or (2) remove the figures from your submission:

a. You may seek permission from the original copyright holder of Figure 1, 3, 4 and 6 to publish the content specifically under the CC BY 4.0 license.

We recommend that you contact the original copyright holder with the Content Permission Form (http://journals.plos.org/plosone/s/file?id=7c09/content-permission-form.pdf) and the following text:

“I request permission for the open-access journal PLOS ONE to publish XXX under the Creative Commons Attribution License (CCAL) CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). Please be aware that this license allows unrestricted use and distribution, even commercially, by third parties. Please reply and provide explicit written permission to publish XXX under a CC BY license and complete the attached form.”

Please upload the completed Content Permission Form or other proof of granted permissions as an "Other" file with your submission.

In the figure caption of the copyrighted figure, please include the following text: “Reprinted from [ref] under a CC BY license, with permission from [name of publisher], original copyright [original copyright year].”

b. If you are unable to obtain permission from the original copyright holder to publish these figures under the CC BY 4.0 license or if the copyright holder’s requirements are incompatible with the CC BY 4.0 license, please either i) remove the figure or ii) supply a replacement figure that complies with the CC BY 4.0 license. Please check copyright information on all replacement figures and update the figure caption with source information. If applicable, please specify in the figure caption text when a figure is similar but not identical to the original image and is therefore for illustrative purposes only.

The following resources for replacing copyrighted map figures may be helpful:

USGS National Map Viewer (public domain): http://viewer.nationalmap.gov/viewer/

The Gateway to Astronaut Photography of Earth (public domain): http://eol.jsc.nasa.gov/sseop/clickmap/

Maps at the CIA (public domain): https://www.cia.gov/library/publications/the-world-factbook/index.html and https://www.cia.gov/library/publications/cia-maps-publications/index.html

NASA Earth Observatory (public domain): http://earthobservatory.nasa.gov/

Landsat: http://landsat.visibleearth.nasa.gov/

USGS EROS (Earth Resources Observatory and Science (EROS) Center) (public domain): http://eros.usgs.gov/#

Natural Earth (public domain): http://www.naturalearthdata.com/

[Note: HTML markup is below. Please do not edit.]

Reviewers' comments:

Reviewer's Responses to Questions

Comments to the Author

1. Is the manuscript technically sound, and do the data support the conclusions?

The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented.

Reviewer #1: Partly

Reviewer #2: Yes

**********

2. Has the statistical analysis been performed appropriately and rigorously?

Reviewer #1: No

Reviewer #2: Yes

**********

3. Have the authors made all data underlying the findings in their manuscript fully available?

The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified.

Reviewer #1: No

Reviewer #2: Yes

**********

4. Is the manuscript presented in an intelligible fashion and written in standard English?

PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.

Reviewer #1: No

Reviewer #2: No

**********

5. Review Comments to the Author

Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)

Reviewer #1: This paper was used SEBAL model to estimate ET by Landsat imagery. The topic is quite interesting and could be stimulating for the scientific community after the major revision.

1. As a general comment, I can say that the paper/work is still in an initial phase of implementation, and deeper analysis is needs.

2. The novelty, originality, and scientific contribution of this study in its respective field is not well presented. I think authors should clarify why the study is important, i.e., what is the contribution of this study inside of its branch?

3. Can Authors explain the verification section in more detail?

- How many ground stations were used for validation? (Line 172 mentioned 4. But 8 station were mentioned between line 180-181)

- What is temporal resolution of station data?

- What is error of the data verification? How the author prevents the errors of data validation?

4. E-601 and 20cm-caliber have proven “meet the requirements” (line 185), why only 20cm-caliber data was selected as the validation data?

5. Please please clarify the fig 2 in more detail. Author mentioned 4 stations, investigated 6 days, why were 12 “point” used in the competition?

6. The results section is not clear and need to be improved.

- The solid explanation should be given instead of rather generic explanations are provided (“Spatiotemporal distribution of daily ET” section).

- There are many surface parameters affect to ET, such as wind speed, air humidity…. Why four factor, including NDVI, LST, land surface albedo, net radiation flux, were selected?

- NDVI, LST, land surface albedo, net radiation flux are inputs the SEBAL model, so it's a problem of measuring similarities.

Reviewer #2: The paper is a straightforward analysis for invert the regional ET based on Landsat8 images and compare with observed ET and explore the mechanisms of surface parameters on daily ET. I recommend publication after the following minor revisions.

(1) The latent heat flux should be calculated by P-M equation

(2) The author should clarify the distinction between daily ET and latent heat flux

(3) It is interesting to explore the mechanism of how hydrothermal conditions effect the energy transmission.

(4) Please explain why the results of albedo v.s. ET shows a larges bias

**********

6. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #1: No

Reviewer #2: No

[NOTE: If reviewer comments were submitted as an attachment file, they will be attached to this email and accessible via the submission site. Please log into your account, locate the manuscript record, and check for the action link "View Attachments". If this link does not appear, there are no attachment files.]

While revising your submission, please upload your figure files to the Preflight Analysis and Conversion Engine (PACE) digital diagnostic tool, https://pacev2.apexcovantage.com/. PACE helps ensure that figures meet PLOS requirements. To use PACE, you must first register as a user. Registration is free. Then, login and navigate to the UPLOAD tab, where you will find detailed instructions on how to use the tool. If you encounter any issues or have any questions when using PACE, please email PLOS at figures@plos.org. Please note that Supporting Information files do not need this step.

8 Dec 2021

Response to reviewers’ comments

We would like to thank the academic editor and reviewers for their constructive and insightful comments and suggestions. We have revised the manuscript accordingly. To facilitate the reviewers to see the changes we made in response to their comments, the changes are marked in blue in the revised manuscript. Below we provide a point to point response to all the comments.

Academic Editor

Question 1: Please ensure that your manuscript meets PLOS ONE's style requirements, including those for file naming. The PLOS ONE style templates can be found at https://journals.plos.org/plosone/s/file?id=wjVg/PLOSOne_formatting_sample_main_body.pdf and https://journals.plos.org/plosone/s/file?id=ba62/PLOSOne_formatting_ sample_title_authors_affiliations.pdf

Response: Thanks for the suggestion. We have revised the whole article according to PLOS ONE style templates.

Question 2: Thank you for stating the following financial disclosure:

This research was funded by Technology Innovation Center for Land Engineering and Human Settlements, Shaanxi Land Engineering Construction Group Co., Ltd and Xian Jiaotong University (2021WHZ0088, 2021WHZ0091) and Scientific Research Item of Shaanxi Provincial Land Engineering Built-up Group (DJNY2021-33).

Please state what role the funders took in the study. If the funders had no role, please state: "The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript."

If this statement is not correct you must amend it as needed.

Please include this amended Role of Funder statement in your cover letter; we will change the online submission form on your behalf.

Response: Thanks for the suggestion. We have updated the Funding Statement and Competing Interests Statement in cover letter.

Funding: This paper was supported by the Fund for Less Developed Regions of the National Natural Science Foundation of China (No. 42167039), Technology Innovation Center for Land Engineering and Human Settlements, Shaanxi Land Engineering Construction Group Co., Ltd and Xian Jiaotong University (2021WHZ0088, 2021WHZ0091) and Scientific Research Item of Shaanxi Provincial Land Engineering Built-up Group (DJNY2021-33). The funder provided support in the form of salaries for authors [Yang LY, Li JF, Sun ZH, Liu JB], but did not have any additional role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. The specific roles of these authors are articulated in the ‘author contributions’ section.

Competing interests: Institute of Land Engineering and Technology, Shaanxi Provincial Land Engineering Construction Group Co., Ltd and Shaanxi Provincial Land Engineering Construction Group Co., Ltd do not alter our adherence to PLOS ONE policies on sharing data and materials. The authors have declared that no competing interests exist.

Question 3: Thank you for stating the following in the Acknowledgments / Funding Section of your manuscript:

This research was funded by Technology Innovation Center for Land Engineering and Human Settlements, Shaanxi Land Engineering Construction Group Co., Ltd and Xian Jiaotong University (2021WHZ0088, 2021WHZ0091) and Scientific Research Item of Shaanxi Provincial Land Engineering Built-up Group (DJNY2021-33).

Please note that funding information should not appear in the Acknowledgments section or other areas of your manuscript. We will only publish funding information present in the Funding Statement section of the online submission form.

Please remove any funding-related text from the manuscript and let us know how you would like to update your Funding Statement. Currently, your Funding Statement reads as follows:

This research was funded by Technology Innovation Center for Land Engineering and Human Settlements, Shaanxi Land Engineering Construction Group Co., Ltd and Xian Jiaotong University (2021WHZ0088, 2021WHZ0091) and Scientific Research Item of Shaanxi Provincial Land Engineering Built-up Group (DJNY2021-33).

Please include your amended statements within your cover letter; we will change the online submission form on your behalf.

Response: Thanks for the suggestion. We have updated the Funding Statement and Competing Interests Statement in cover letter. And remove any funding-related text from the manuscript.

Funding: This paper was supported by the Fund for Less Developed Regions of the National Natural Science Foundation of China (No. 42167039), Technology Innovation Center for Land Engineering and Human Settlements, Shaanxi Land Engineering Construction Group Co., Ltd and Xian Jiaotong University (2021WHZ0088, 2021WHZ0091) and Scientific Research Item of Shaanxi Provincial Land Engineering Built-up Group (DJNY2021-33). The funder provided support in the form of salaries for authors [Yang LY, Li JF, Sun ZH, Liu JB], but did not have any additional role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. The specific roles of these authors are articulated in the ‘author contributions’ section.

Competing interests: Institute of Land Engineering and Technology, Shaanxi Provincial Land Engineering Construction Group Co., Ltd and Shaanxi Provincial Land Engineering Construction Group Co., Ltd do not alter our adherence to PLOS ONE policies on sharing data and materials. The authors have declared that no competing interests exist.

Question 4: We note that Figure 1, 3, 4 and 6 in your submission contain [map/satellite] images which may be copyrighted. All PLOS content is published under the Creative Commons Attribution License (CC BY 4.0), which means that the manuscript, images, and Supporting Information files will be freely available online, and any third party is permitted to access, download, copy, distribute, and use these materials in any way, even commercially, with proper attribution. For these reasons, we cannot publish previously copyrighted maps or satellite images created using proprietary data, such as Google software (Google Maps, Street View, and Earth). For more information, see our copyright guidelines: http://journals.plos.org/plosone/s/licenses-and-copyright.

Response: Thanks for the suggestion. The remote sensing images used in this study are all freely available at the Geospatial Data Cloud (http://www.gscloud.cn). We have explained the acquisition of Data sources of the revised manuscript.

Reviewer #1 This paper was used SEBAL model to estimate ET by Landsat imagery. The topic is quite interesting and could be stimulating for the scientific community after the major revision.

Questions 1: As a general comment, I can say that the paper/work is still in an initial phase of implementation, and deeper analysis is needs.

Response: We are extremely grateful to reviewer for pointing out this problem. To be more clearly and in accordance with the reviewer concerns, we have added a more detailed discussion of meteorological factors and hydrothermal conditions. In future research, we will deeply analyze the mechanism of remote sensing inversion of evapotranspiration and improve the accuracy of model inversion. It can provide data support for the rational use of water resources, restoration of the ecological environment, and protection of ecological water sources. Thank you again for your support and affirmation of our work.

Questions 2: The novelty, originality, and scientific contribution of this study in its respective field is not well presented. I think authors should clarify why the study is important, i.e., what is the contribution of this study inside of its branch?

Response: Evapotranspiration is an important part of the water cycle. Remote sensing technology provides a new method for the estimation of regional land surface water evapotranspiration. This study uses Landsat data and energy balance models to estimate the evapotranspiration in the agro-pastoral ecotone, which solves the problem of transient remote sensing to estimate all-day evapotranspiration. Through the simulation of the daily evapotranspiration in the study area, the distribution law of evapotranspiration in this area is analyzed, which has certain support and guiding significance for the study of land surface process. In the context of increasingly scarce global water resources, accurate estimation of evapotranspiration is not only of great significance for theoretical research on global climate evolution, environmental issues, and evaluation of water resources, but also for guiding agricultural drainage and irrigation, monitoring agricultural drought, and improving agricultural utilization rate of water resources is also very practical significance.

Questions 3: Can Authors explain the verification section in more detail?

- How many ground stations were used for validation? (Line 172 mentioned 4. But 8 station were mentioned between line 180-181)

- What is temporal resolution of station data?

- What is error of the data verification? How the author prevents the errors of data validation?

Response: Our deepest gratitude goes to you for your careful work and thoughtful suggestions that have helped improve this paper substantially. There are four weather stations in the study area, plus 4 surrounding weather stations, a total of eight stations. Since the meteorological sites were small meteorological stations that were established by research teams in Wushenqi and Henan in China, they did not have evaporation dishes to obtain the data on water surface evaporation. Therefore, E-601 and 20cm-caliber meteorological stations, such as those at the Hengshan, Jingbian, Shenmu, Dingbian, Yanchi and Etuoke Banner Stations in China were used for validation in the study and surrounding areas. The water surface evaporator of the weather station has two types that include E-601 and 20cm-caliber. There are differences in the installation plan, structure and observational methods of the two instruments. The two monitoring results are quite different and cannot be directly used to verify the accuracy of the evapotranspiration inversion results. Therefore, to more accurately explore the inversion accuracy of the SEBAL model and maintain the continuity and consistency of the data at each station, it is necessary to explore the conversion coefficients of the two surface evaporators and unify the data. So as to avoid the limitations brought by the data itself.

Questions 4: E-601 and 20cm-caliber have proven “meet the requirements” (line 185), why only 20cm-caliber data was selected as the validation data?

Response: We are very sorry for the inconvenience and confusion caused by my presentation. At present, most weather stations in China have two types of evaporating pans, E-601 and 20cm-caliber, which leads to inconsistent evapotranspiration data. In order to meet the needs, I explored the conversion coefficients of the two by using the weather stations around the study area. Data unification. And because the E-601 evaporating pan monitoring data is closer to the water surface evapotranspiration, this study converted the 20cm-caliber data into E-601 data, and used E-601 data as the verification data. In order to avoid unnecessary mistakes, I have sought out English-speaking experts with professional backgrounds to polish the language.

Questions 5: Please please clarify the fig 2 in more detail. Author mentioned 4 stations, investigated 6 days, why were 12 “point” used in the competition?

Response: Since the meteorological sites were small meteorological stations that were established by research teams in Wushenqi and Henan in China, they did not have evaporation dishes to obtain the data on water surface evaporation. The data involved in the validation were the monitoring data from the meteorological stations with a total of 12 images in 6 scenes in Hengshan and Jingbian, China. We have added explanations of accuracy verification of daily ET in Section 3.1 of the revised manuscript.

Questions 6: The results section is not clear and need to be improved.

- The solid explanation should be given instead of rather generic explanations are provided (“Spatiotemporal distribution of daily ET” section).

- There are many surface parameters affect to ET, such as wind speed, air humidity…. Why four factor, including NDVI, LST, land surface albedo, net radiation flux, were selected?

- NDVI, LST, land surface albedo, net radiation flux are inputs the SEBAL model, so it's a problem of measuring similarities

Response: The suggestion is taken. Our deepest gratitude goes to you for your careful work and thoughtful suggestions that have helped improve this paper substantially. We re-organized the content of the “Spatiotemporal distribution of daily ET” section, and explained the seasonal variation process and spatial distribution of daily evapotranspiration. In summary, hydrothermal conditions are the main factor that affect the changes in evapotranspiration during the year. From a local perspective, the daily ET of different land use types differs significantly. The characteristics of the underlying surface are important factors that affect the exchange of energy and material on the surface. NDVI, surface temperature, surface albedo and net surface radiation are several important parameters describing the nature of the underlying surface. They interact and jointly affect the spatial distribution pattern of evapotranspiration. Therefore, I chose NDVI, LST, land surface albedo and net radiation flux to explore the influencing factors of evapotranspiration. In order to better clarify the factors affecting evapotranspiration, we added Sensitivity analysis of model parameters, including hot spot surface temperature and water vapor. We have explained the details in Sections 4.2 and 4.3 of the revised manuscript.

Reviewer #2: The paper is a straightforward analysis for invert the regional ET based on Landsat8 images and compare with observed ET and explore the mechanisms of surface parameters on daily ET. I recommend publication after the following minor revisions.

Questions 1: The latent heat flux should be calculated by P-M equation

Response: The P-M formula is a standard formula for calculating the reference crop water requirement. The P-M can be used to calculate the reference crop evapotranspiration at a single point, and sufficient meteorological parameters are required as support. The latent heat flux obtained in this study is regional, and the meteorological data is not abundant, so the SEBAL model calculation is more suitable.

Questions 2: The author should clarify the distinction between daily ET and latent heat flux.

Response: Latent Heat Flux is the heat exchange per unit area under the condition of constant temperature, in watts per square meter (W/m2); daily ET is the amount of water actually evaporated in a day (mm/d), which is the result of dividing the total latent heat flux of the whole day by the latent heat of vaporization (λ) of water. We have added explanations of the latent heat flux and daily ET in Section 2.3.4 and2.3.5 of the revised manuscript.

Questions 3: It is interesting to explore the mechanism of how hydrothermal conditions effect the energy transmission.

Response: Thanks for the suggestion. The suggestion is taken. This research has added Sensitivity analysis of model parameters in the discussion section. The influence of hotspot surface temperature and water vapor content on the model results is discussed, and the importance of hydrothermal conditions to SEBAL is further proved. The results show that the model output results on the LSTh are less sensitive compared with those of the wv. The main reason for this is that the LSTh is used to calculate the sensible heat flux, and the result of the sensible heat flux requires calculations that are continuously iterative. Multiple iterations weaken the influence of the LSTh. Thus, the change of a single parameter of the LSTh does not have a large influence on the daily ET. See the discussion section in the text for detailed analysis results.

Questions 4: Please explain why the results of albedo v.s. ET shows a larges bias

Response: The surface albedo determines the amount of solar energy received at the pixel scale. The higher the surface albedo, the smaller the effective radiation reaching the surface, the less energy used for evapotranspiration, and the smaller the evapotranspiration. Therefore, there is a negative correlation between surface albedo and evapotranspiration, which is consistent with the results reflected in Figure 4. However, ET is affected by both energy and moisture. In arid areas, moisture is the decisive factor limiting the amount of surface evapotranspiration. Water in the study area is not sufficient, which limits the amount of evapotranspiration, which results in a large deviation between the surface albedo and ET.

Submitted filename: Response to Reviewers.docx

Click here for additional data file.

9 Feb 2022

Submitted filename: PONE-D-21-26971R2.docx

Click here for additional data file.

13 Feb 2022

Response to reviewers’ comments

We would like to thank the academic editor and reviewers for their constructive and insightful comments and suggestions. We have revised the manuscript accordingly. To facilitate the reviewers to see the changes we made in response to their comments, the changes are marked in blue in the revised manuscript. Below we provide a point-to-point response to all the comments.

Academic Editor

Question 1: Please review your reference list to ensure that it is complete and correct. If you have cited papers that have been retracted, please include the rationale for doing so in the manuscript text, or remove these references and replace them with relevant current references. Any changes to the reference list should be mentioned in the rebuttal letter that accompanies your revised manuscript. If you need to cite a retracted article, indicate the article’s retracted status in the References list and also include a citation and full reference for the retraction notice.

Response: Our deepest gratitude goes to you for your careful work and thoughtful suggestions that have helped improve this paper substantially. We have revised the reference list to ensure that it is complete and correct. The changes of the reference list are marked in blue in the Revised Manuscript with Track Changes.

Reviewer：Dear authors, I do appreciate your further efforts in producing an advanced version of your paper that takes into account my comments.

The paper has improved. I suggest accepting it, after the authors make the following necessary revision:

Questions 1: Normalized Difference Latent Heat (NDLI) index was proposed and proved as a reference for the assessment and monitoring water availability. you are encouraged to use it for your validation.

Response: Thank you again for your support and affirmation of our work. Thank you very much for your suggestion. Since we have not used the Normalized Difference Latent Heat (NDLI) index in past research, this study has not used it for verification. Thank you for pointing us in new directions. And in future research work, we will introduce the Normalized Difference Latent Heat (NDLI) index to improve our research results. Thank you again for your affirmation and support of this research.

Questions 2: Line 30 : “… coefficients of normalized vegetation indices…” Only NDVI was used for valiation as spectral indices…so used as singular noun.

Response: We are extremely grateful to reviewer for pointing out this problem. We have fixed this question. “… coefficients of normalized vegetation indices…” modified to “… coefficients of normalized vegetation index…”. Thank you again for your support and affirmation of our work.

Questions 3: Please, improve the quality of the figure (as grip, scale…).

Response: Thanks for the suggestion. We have improved the quality of Figures 1, 2, 4, and 6, where Figure 1 has modified the scale, and Figures 2 and 4 have increased the resolution of the figures.

Submitted filename: Response to Reviewers-R2.docx

Click here for additional data file.

24 Feb 2022

Daily actual evapotranspiration estimation of different land use types based on SEBAL model in the agro-pastoral ecotone of northwest China

PONE-D-21-26971R2

Dear Dr. Yang,

We’re pleased to inform you that your manuscript has been judged scientifically suitable for publication and will be formally accepted for publication once it meets all outstanding technical requirements.

Within one week, you’ll receive an e-mail detailing the required amendments. When these have been addressed, you’ll receive a formal acceptance letter and your manuscript will be scheduled for publication.

An invoice for payment will follow shortly after the formal acceptance. To ensure an efficient process, please log into Editorial Manager at http://www.editorialmanager.com/pone/, click the 'Update My Information' link at the top of the page, and double check that your user information is up-to-date. If you have any billing related questions, please contact our Author Billing department directly at authorbilling@plos.org.

If your institution or institutions have a press office, please notify them about your upcoming paper to help maximize its impact. If they’ll be preparing press materials, please inform our press team as soon as possible -- no later than 48 hours after receiving the formal acceptance. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information, please contact onepress@plos.org.

Kind regards,

Débora Regina Roberti, Ph.D.

Academic Editor

PLOS ONE

Additional Editor Comments (optional):

Reviewers' comments:

4 Mar 2022

PONE-D-21-26971R2

Daily actual evapotranspiration estimation of different land use types based on SEBAL model in the agro-pastoral ecotone of northwest China

Dear Dr. Yang:

I'm pleased to inform you that your manuscript has been deemed suitable for publication in PLOS ONE. Congratulations! Your manuscript is now with our production department.

If your institution or institutions have a press office, please let them know about your upcoming paper now to help maximize its impact. If they'll be preparing press materials, please inform our press team within the next 48 hours. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information please contact onepress@plos.org.

If we can help with anything else, please email us at plosone@plos.org.

Thank you for submitting your work to PLOS ONE and supporting open access.

Kind regards,

PLOS ONE Editorial Office Staff

on behalf of

Dr. Débora Regina Roberti

Academic Editor

PLOS ONE