Next Article in Journal
Attribution Analysis of Runoff Variation in the Second Songhua River Based on the Non-Steady Budyko Framework
Next Article in Special Issue
Characterization of Evapotranspiration in the Orange River Basin of South Africa-Lesotho with Climate and MODIS Data
Previous Article in Journal
Monitoring of Plant Cultivation Conditions Using Ground Measurements and Satellite Products
Previous Article in Special Issue
Correlation between Ground Measurements and UAV Sensed Vegetation Indices for Yield Prediction of Common Bean Grown under Different Irrigation Treatments and Sowing Periods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of Actual Evapotranspiration Using Satellite-Based Surface Energy Balance Derived from Landsat Imagery in Northern Thailand

by
Teerawat Suwanlertcharoen
1,
Thodsapol Chaturabul
2,
Thanaporn Supriyasilp
3,* and
Kobkiat Pongput
1
1
Department of Water Resources Engineering, Faculty of Engineering, Kasetsart University, Bangkhen Campus, Bangkok 10900, Thailand
2
Department of Civil and Environmental Engineering, Faculty of Science and Engineering, Kasetsart University, Chalermphrakiat Sakon Nakhon Province Campus, Sakon Nakhon 47000, Thailand
3
Department of Civil Engineering, Faculty of Engineering, Chiang Mai University, Chiang Mai 50200, Thailand
*
Author to whom correspondence should be addressed.
Water 2023, 15(3), 450; https://doi.org/10.3390/w15030450
Submission received: 19 December 2022 / Revised: 10 January 2023 / Accepted: 18 January 2023 / Published: 22 January 2023

Abstract

:
In this study, satellite-based measures of surface energy balance and the mapping evapotranspiration at high resolution with internalized calibration (METRIC) from Landsat imagery were used to estimate the spatiotemporal distribution of actual evapotranspiration (ETa) in northern Thailand, constituting a procedure that has rarely been performed in southeast Asia. Subsequently, we compared the ETa obtained from METRIC with that calculated using the FAO-56 dual-crop coefficient method via the SIMDualKc software and found a strong correlation. An assessment of the accuracy of all the sample plots revealed the R2, Root-Mean-Square Error (RMSE), and mean absolute error (MAE) values to be 0.830, 0.730, and 0.575 mm d−1, respectively. Differences in the cumulative ETa values derived from SIMDualKc and METRIC ranged in magnitude from 0.93–3.57% for rice and 3.08–7.99% for longan. The ETa values for forestland and waterbodies were higher than those for agricultural areas and areas with other forms of land use. The spatiotemporal distribution of the seasonal ETa during the dry season was consistent with the climate, vegetation, and anthropogenic activity. Thus, our results indicate that METRIC is a reliable tool for estimating ETa for water resource management under different environmental conditions and improving water use efficiency over large areas.

1. Introduction

Southeast Asia is one of the world’s most vital agricultural regions, accounting for 15% of food production in Asia and 7.7% globally [1]. Thailand, located in the south-eastern region of the continent of Asia, covers an area of 513,120 km2, with 43% of the country devoted to agricultural land [2]. Thailand is the foremost longan-producing country globally and produced approximately 0.98 million metric tons of longan from 2015–2017 [3].
At present, many parts of the world are affected by water scarcity, including Thailand. Higher food demands and population growth have directly affected water use in the agricultural and other sectors. Additionally, climate change has increased water scarcity and drought, with potential impacts including increased evapotranspiration and decreased water availability [4]. Consequently, the demand for water in agriculture and for irrigation water use have increased.
Longan orchards occupy a considerable portion of the agricultural land area in northern Thailand. Longan trees are susceptible to drought during the flowering and early fruit development stages, which generally occur during the dry season from the end of October to mid-May [5]. Premium fruit can only be produced under intensive water and nutrient management. In addition to longan, rice is commonly cultivated in this area, and also requires a large amount of water during the dry season, especially in flooded paddies. The Upper Ping River Basin, which is one of the headwaters of the Chao Phraya River Basin, is located in northern Thailand. To understand crop water consumption over large areas and achieve effective water resource management, it is necessary to estimate the spatiotemporal distribution of actual evapotranspiration (ETa) in agricultural areas, especially in longan orchards and paddy fields.
Evapotranspiration (ET) is an important component of water balance associated with the hydrological cycle and biological processes [6] and represents the total water loss from the Earth’s surface, including water evaporated from the soil surface and that has transpired through plants. Therefore, the accurate estimation of the actual ET (ETa) rate is crucial for understanding changes in water availability and sustainably managing water resources.
According to Allen et al. [7], ET data can be obtained using various measurement techniques, such as lysimeters, the Bowen ratio, the eddy covariance, sap flow, root zone soil water balance, and even satellite-based remote sensing, as well as direct modelling. Among these techniques, direct measurement techniques necessitate extensive experimental care and are susceptible to biases based on reported results [7]. Satellite imagery-based techniques can be used to estimate ET using the energy balance and can potentially be applied over large areas to identify regions characterised by ET reductions due to water stress [7]. The surface energy balance algorithm for land (SEBAL) [8] and mapping evapotranspiration at high resolution with internalized calibration (METRIC) [9] models are two widely used satellite-based surface energy balance techniques in which thermal infrared sensors are utilised to estimate ETa. The METRIC model, which is a variant of the image-processing tool SEBAL, is used to estimate ETa from the energy balance residual at the Earth’s surface. Moreover, METRIC has been enhanced to provide better integration with ground weather station data used to compute reference ET and has been applied to satellite images obtained using thermal infrared sensors such as the moderate resolution imaging spectroradiometer (MODIS) and Landsat [10]. Remote-sensing-based energy balance models have major limitations, such as time gaps, because remote-sensing images can only be obtained at a specific location on a periodic basis. Therefore, the effects of evaporation caused by precipitation events that occur between these time gaps may be overlooked [7]. Another technique that uses remote-sensing image data involves satellite-based vegetation indices (VIs). In this method, VIs obtained from image data are used to estimate the crop coefficient (Kc) and ‘basal’ crop coefficient (Kcb) by determining the relationship between the VI and these coefficients. A typical VI is the normalised difference vegetation index (NDVI). The ‘basal’ Kcb condition is a condition under which the soil surface is dry enough to reduce evaporation from the soil surface to relatively low levels compared with transpiration. However, transpiration continues to occur [7]. A reliable soil water balance model would be able to estimate soil water behaviour over time, thereby helping to overcome the majority of the difficulties in estimating ET from observations of soil water via direct modelling [7]. The advantages and disadvantages of these ET measurement techniques have previously been discussed by Allen et al. [7].
In Thailand, several approaches have been used to estimate ET with satellite-based remote sensing at spatial resolutions from regional to global scales. Zheng et al. [11] estimated ET with satellite-based Earth observation datasets as the primary inputs to derive daily ET in north-eastern Thailand. Sriwongsitanon et al. [12] validated ET values from seven global remote-sensing products used for water resource management in Thailand. Furthermore, other researchers have applied satellite-based surface energy balance using the SEBAL model to estimate ETa from MODIS and Landsat-7 images in the Chao Phraya River Basin [13] and from Landsat-5 images in the Lam Ta Kong Basin [14]. However, to date, no studies using high-resolution satellite-based images and employing the METRIC model, which can provide better integration with data from ground weather stations than SEBAL for ETa estimation at the field scale, have been conducted in Thailand. Therefore, it is of interest to study the accuracy of ETa estimation using METRIC and Landsat-8 images in this region.
To our knowledge, few researchers have studied the accuracy of satellite-based ETa values obtained using the soil water balance method. However, Paço et al. [15] analysed and compared the ET and crop coefficients of olive orchards obtained using the SIMDualKc and METRIC algorithms. They reported that the ETa values obtained from the SIMDualKc and METRIC algorithms were similar, and that the crop coefficients derived from each model showed comparable patterns throughout the year. Thus, although the SIMDualKc and METRIC algorithms represent different approaches, their results appear to be comparable and complementary with respect to spatiotemporal scale. Furthermore, Tasumi [16] compared the ETa value obtained from METRIC to that estimated independently using the FAO-56 single-crop approach. However, ETa estimation based on surface energy balance using the METRIC model with a 30 m spatial resolution has rarely been performed in Southeast Asia, and the dual-crop coefficient method and SIMDualKc model have seldom been utilised for the evaluation of ETa, especially in Thailand.
Herein, we focus on the estimation of ET over a large area for water management purposes using two approaches. First, we combine the METRIC algorithm and NDVI to overcome the time gap limitation and use Landsat images, which currently represent the highest resolution (a spatial resolution of 30 m) over a large area with a revisit time of 16 d. Thermal Landsat imagery can be used to evaluate vegetation phenology in the growing period and observe changes in the water stress of plants [17]. In the second approach, we apply direct soil water balance modelling using SIMDualKc. The SIMDualKc model employs the FAO-56 dual-crop coefficient method to estimate the crop ET. Then, we compare the results obtained with the two approaches.
The objectives of this study are as follows: (1) to estimate ETa during each crop growth period using the METRIC model and Landsat-8 satellite data; (2) to evaluate the ETa obtained from the METRIC model relative to the ETa calculated from the FAO-56 dual-crop coefficient method; and (3) to study the spatiotemporal distribution of the ETa for various land-use types in order to enhance water resource management in northern Thailand.
The ETa obtained from the satellite-based surface energy balance data derived from Landsat imagery can be applied in water management schemes based on appropriate technology at field to regional scales, thereby covering a large area. These data are beneficial with respect to monitoring crop water use and stress. Furthermore, ETa estimates based on the METRIC model can be used for water resource management planning and the localisation of water to meet irrigation water requirements in the dry season, especially in cultivation areas on both sides of the Ping River in the Chiang Mai and Lamphun provinces. The accurate estimation of the spatiotemporal distribution of ETa in the Upper Ping Basin will allow for the establishment of better irrigation strategies, enhanced water use efficiency, and increased crop production yields.

2. Materials and Methods

2.1. Study Area

The research area was divided into two sections: the study area, which was the entire area used to study the spatiotemporal distribution of the ETa across various land-use types, and the testing area, which was used to compare the ETa obtained from the METRIC model and the FAO-56 dual-crop coefficient method (Figure 1).
The study area (~12,900 km2 located in the plains of the upper northern region of Thailand along the Ping River (Figure 1b)), includes a variety of land utilisation areas, especially agricultural zones in irrigated areas, such as government, private, and community irrigation systems, and non-irrigated agricultural areas in the Chiang Mai and Lamphun provinces. Most of the agricultural portions consist of longan orchards and paddy fields. The general climate of this area is affected by southwest and northeast monsoons, and the annual climate is influenced by depression from the South China Sea between July and September, causing heavy rain from May to October. The mean annual temperature is 26.1 °C. The highest and lowest temperatures are 30 °C (in April) and 21.3 °C (in December), respectively. The average annual precipitation is 993.2 mm [18]. Types of land use in the study area include forestland, orchards, urban and built-up land, paddy fields, miscellaneous land, field crops, perennial trees, waterbodies, and horticulture, accounting for 65.88%, 10.20%, 7.02%, 6.95%, 4.49%, 2.61%, 1.10%, 1.08%, and 0.67% of the entire area, respectively.
To evaluate the ETa obtained from the METRIC model, the FAO-56 dual-crop coefficient method was used. The testing area used for ETa evaluation is part of the Muang Luang Irrigation Project located in the Chom Thong District, Chiang Mai Province (Figure 1c). The Mueang Luang Weir is located upstream of this area. Paddy fields and longan orchards cover ~ 68% and 26% of the total area, respectively, and the Tae Sam Nga diversion weir, which splits into three irrigation canals, occupies the northern boundary. The testing area covers 2.12 km2 and is located between the Na Sai and Na Ngu ditches [19].

2.2. Data

2.2.1. Satellite Imagery

In this study, Landsat-8 Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS) images with a 30 m spatial resolution were collected from level-1 products with terrain precision correction (L1TP) using U.S. Geological Survey (USGS) Earth Explorer. In total, 41 Landsat-8 OLI/TIRS images (Path 131/Row 47) were collected under cloud-free conditions throughout the crop growth period in the dry season from October 2017 to May 2020 as shown in Table 1 (12, 14, and 15 images for 2018, 2019, and 2020, respectively).

2.2.2. Meteorological Data

Meteorological data for this study were obtained from the closest weather station in order to compute the standardised alfalfa reference ET equation (ETr) and provide the inputs required for METRIC computations. The weather station records meteorological data every 10 min, including precipitation, air temperature, wind speed, relative humidity, and atmospheric pressure. Data were collected at the Lamphun station using the Automatic Weather System (AWS) of the Thailand Meteorological Department. The only climate data not available at this station were solar radiation data. Therefore, solar radiation data were derived from the air temperature difference using the Hargreaves–Samani radiation formula, which is an option available for calculating solar radiation on the REF-ET software. These meteorological data, except for precipitation, were also used to compute the daily standardised grass reference ET (ETo) for the Muang Luang testing area. Precipitation data for this area were collected at Chom Thong station. The ETo and precipitation data were then used as inputs for SIMDualKc. Precipitation data were obtained from different stations. For the study area (red square in Figure 1a), precipitation data from Lamphun station were chosen for METRIC computations. For the testing area, precipitation data from Chom Thong station were used.

2.2.3. Geographic Information System (GIS) Data

The GIS data required for this study consisted of land-use, irrigated area, and digital elevation model (DEM) data. Land-use data were obtained from the Land Development Department, whereas data regarding the irrigated area were obtained from the Royal Irrigation Department. For the DEM, Shuttle Radar Topography Mission (SRTM) data with 30 m spatial resolution were used, which were obtained from the USGS Earth Resources Observation and Science (EROS) Centre. The DEM was then used to calculate the slope and aspect of the area. The slope and aspect were later combined with the latitude, declination, hour angle, and solar incidence to obtain solar angles, which were then utilised to calculate the incoming shortwave radiation (RS↓) and net radiation (Rn), respectively.

2.2.4. Soil Sampling

Soil samples were taken from five plots to obtain soil texture information. The main soil type in most plots (A2_Rice, A5_Rice, and A6_Longan) was sandy loam, while the soil types of the A1_Rice and A4_Longan plots were loam and sandy clay loam, respectively. These soil properties were necessary for the soil water balance calculation and estimation of ETa using the FAO-56 dual-crop coefficient method. The sampling locations and soil properties are summarised in Table 2.

2.2.5. Water Use and Water Measurement in the Testing Area

According to the data obtained by Supriyasilp and Pongput [20], water from various sources is used in the dry season in the testing area, including rainfall, irrigation water supplied from the upstream Muang Luang Irrigation Scheme, irrigation water supplied from two pumping stations on the Ping River, and ponds. The proportions of the total amounts of water derived from each source are 11%, 15%, 57%, 16%, and 1%, respectively. The flow data observed in the dry season were measured from 21 March to 15 May 2019. The agricultural water supply for this testing area is jointly managed by the water management committee and farmers.
Descriptions and sources of the datasets used in this study are summarised in Table 3.

2.3. Software

Four software products were used in this study. Their capabilities are presented in Table 4.

2.4. Overview of Methods

An outline of the methods, procedures, and subsequent workflow used in this study is presented in Figure 2.

2.4.1. METRIC Model

The METRIC model was used to estimate the ETa from remote-sensing data including satellite images at the time of image acquisition. The resulting ETa value is called instantaneous ET. As part of the surface energy balance method, energy consumed by the ET process was calculated as the residual from the surface energy balance equation [9]:
L E = R n G H ,
where LE is the latent heat flux consumed by ET (W·m−2), Rn is the net radiation (W·m−2), G is the soil heat flux conducted into the ground (W·m−2), and H is the sensible heat flux transferred to the atmosphere (W·m−2).
The parameter Rn was calculated from the surface reflectance and temperature at the time of satellite image acquisition by summing incoming and outgoing shortwave and longwave radiation, including solar and thermal radiation [9]:
R n = R S α R S + R L R L ( 1 ε o ) R L
where Rn is the net radiation (W·m−2), RS↓ is the incoming shortwave radiation (W·m−2), α is the surface albedo (dimensionless), RL↓ is the incoming longwave radiation (W·m−2), RL↑ is the outgoing longwave radiation (W·m−2), and εo is the surface thermal emissivity (dimensionless).
Soil heat flux (G) is the magnitude of the soil energy stored or convected into the ground and can be computed as follows [21]:
G R n = 0.05 + 0.18 e 0.521 L A I ,
L A I 0.5 ;
G R n = 180   ( T s 273.16 ) R n + 0.084 ,
L A I < 0.5 ,
where LAI is the leaf area index (dimensionless) and Ts is the land surface temperature (Kelvin).
Sensible heat flux transferred to the atmosphere (H) was calculated from the aerodynamic function-based heat transfer. The METRIC model was estimated based on a calibration using inverse modelling under extreme conditions (CIMEC) [22]. This method involves the selection of pixels with near-extreme conditions (hot and cold anchor pixels) on whose basis ETa can be estimated and assigned. The sensible heat flux was calculated as follows:
H = ρ a i r C p d T r a h
where ρair is the air density (kg·m−3), Cp is the specific heat of air under constant pressure (1004 ·J·kg−1·K−1), dT is the near-surface temperature difference between the two near-surface heights z1 (0.1 m) and z2 (2 m; Kelvin), and rah is the aerodynamic resistance to heat transport (s·m−1).
Note that dT was solved using linear regression between dT and Ts of hot and cold anchor pixels selected from different areas under two extreme conditions. Cold and hot conditions were selected for the full vegetation canopy area and empty agricultural fields without vegetation cover, respectively.
Based on the calculation of the LE process as the residual from the surface energy balance equation, LE was calculated as the ETa during the time of satellite image acquisition:
E T i n s t = 3600 × L E λ ρ w ,
where ETinst is the instantaneous ETa at the moment of satellite image acquisition (mm·h−1), 3600 is the factor for the conversion from seconds to hours, ρw is the density of water (100 kg·m−3), and λ is the latent heat of vaporisation (J·kg−1), which is calculated as λ = [ 2.51 0.00236   ( T s 273.15 ) ] × 10 6 .
In this step, the daily ETa (ET24) of each pixel was calculated as follows:
E T 24 = E T i n s t E T r * E T r 24 ,
where ETinst is the instantaneous ETa at the time of the satellite image acquisition (mm·h−1); ETr is the alfalfa reference ET at the moment of satellite image acquisition, which is calculated at the weather station (mm·h−1); ETr24 is the daily alfalfa reference ET (mm·day−1); and ET24 is the daily ETa (mm·day−1).
The ETinst/ETr ratio is the reference ET fraction (ETrF) and is the same as the well-known Kc [9].
The ETrF at the date of the image acquisition was interpolated with the cubic spline interpolation method to obtain the time series of the daily ETrF. Subsequently, the daily ETa (ET24) was calculated by multiplying the daily ETrF by ETr.
The cumulative ETa can be calculated for any period, such as a month or season, using the time series of the daily ETa (ET24) and the following equation:
E T p e r i o d = i = m n E T r F i × E T r 24 i ,
where ETperiod is the cumulative ETa for the period starting on day m and ending on day n, E T r F i is the reference evapotranspiration fraction for day i, and ETr24i is the daily ETa for day i.

2.4.2. Reference Evapotranspiration Calculation

The standardised Penman–Monteith equation for ETref has been proposed by the American Society of Civil Engineers (ASCE) for both the grass reference (ETo) and alfalfa reference (ETr) according to the format adopted in FAO-56 [6,17]
E T r e f = 0.408 Δ ( R n G ) + γ C n ( T + 273 ) u 2 ( e s e a ) Δ + γ ( 1 + C d u 2 ) ,
where ETref is the reference ET (mm·day−1 or mm·h−1), R n is the net radiation (MJ m−2 d−1 or MJ m−2 h−1), G is the soil heat flux transferred to the ground (MJ m−2 d−1 or MJ m−2 h−1), T is the mean daily or hourly air temperature (°C), u2 is the mean daily or hourly wind speed at a height of 2 m (m·s−1), es is the saturation vapour pressure (kPa), ea is the actual vapour pressure (kPa), ∆ is the slope vapour pressure curve (kPa °C−1), γ is the psychrometric constant (kPa °C−1), and Cn and Cd are coefficients that change depending on the calculation time step (hourly or daily) and reference vegetation type (ETo for grass or ETr for alfalfa).

2.4.3. Dual-Crop Coefficient Approach

The dual-crop coefficient approach is based on the standard method for estimating the crop ET (ETc), which is calculated from ETo and the crop coefficient (Kc). The FAO-56 dual-crop coefficient approach is based on crop transpiration and soil evaporation components [6,23]. The basal crop coefficient (Kcb) and soil evaporation coefficient (Ke) are considered in the dual-crop coefficient approach. In this study, ETc values for paddy rice and longan in the testing area were determined using SIMDualKc software. SIMDualKc is a soil water balance model that estimates ETc using the dual-crop coefficient approach [15,24]. Consequently, the ETc value was calculated as follows:
E T c = ( K s * K c b + K e ) × E T o ,
where Kcb is the basal coefficient, Ke is the soil evaporation coefficient, and Ks is the water stress coefficient. Ks ranges from 0 and 1 depending on the soil water content in the root zone. The Kcb of rice varies between 0.15 and 1.15, as suggested by Allen et al. [6] and Pereira et al. [23], and the Kcb of fruit trees varies between 0.75–1.15, as suggested by Allen et al. [6]. Date of each growing stage was obtained from the crop phenology identified by the NDVI. For the computation of Ke, soil data and soil evaporable layer data are required to compute the soil water balance in the root zone and the surface soil layer. As suggested by Allen et al. [6], the SIMDualKc model provides values for the total available soil water, soil water content at field capacity, wilting point, and bulk density based on soil texture to assist users in selecting appropriate values for these parameters.

2.4.4. Estimation of Actual Evapotranspiration Based on the METRIC Model

Calculation of the Reference ET

In this step, weather data from the AWS at the Lamphun station were quality-assessed and quality-controlled before being used for the ET calculation. The weather data were then used as input for the REF-ET software to obtain the ETref equation at the time of satellite image acquisition, with hourly (for ETr) and daily (for ETr and ETo) time steps.

METRIC Model Processing

The METRIC model was applied to estimate the ETa from Landsat-8 images taken between October 2017 and May 2020, which covers the growth period in the dry season. Cloud masking was the first process applied to the satellite images in each scene. Hot pixels, representing empty agricultural fields without vegetation cover, and cold pixels, representing areas with full vegetation canopy, were selected for the study area. METRIC corresponds to the satellite-based surface energy balance. Thus, the latent heat flux (LE) was computed by subtracting the soil heat flux (G) and sensible heat flux (H) from the net radiation (Rn). The LE was computed for each pixel at the instant the satellite overpassed the monitored location and converted to ETinst [9]. The ETrF was then calculated as the ratio of ETinst to ETr at the time the satellite overpassed the monitored location and assumed to be equivalent to the 24 h average. Finally, ET24 was calculated by multiplying ETrF by ETr using the ETref equation.

Filling Image Gaps and Interpolation of the Daily Image

ETrF and NDVI were used to fill data gaps due to cloud cover, provide missing values, and generate the daily ETrF image between two dates of satellite image captures using the cubic spline interpolation method. Time series of daily ETrF images was used to calculate the daily ETa, and the NDVI was applied to extract crop phenology data and identify crop development stages, including initial, developing, middle, and late stages.

Spatiotemporal Distribution of ETa for Various Land Use Types

Spatiotemporal distributions of ETa for various land-use types were generated by multiplying the daily image of ETrF by ETr for each day. The output was ETa, which was then accumulated over the period from October to May, covering the dry seasons in 2018, 2019, and 2020.

2.4.5. Evaluation of the Actual Evapotranspiration

In this study, the Muang Luang Irrigation Project was used as the field site for the testing area. Within this testing area, soil samples were collected, and soil properties were used as inputs for the SIMDualKc software, in addition to climate, crop, and irrigation data. For the crop data, each crop’s planting date and crop growth period were required. This information was obtained from the crop phenology identified by the NDVI. The ETa values obtained from METRIC and SIMDualKc were used for the evaluation. Finally, the accuracy of the ETa values obtained from the METRIC and FAO-56 dual-crop coefficient methods was assessed using statistical analysis methods and the following quality measures: coefficient of determination (R2), Root-Mean-Square Error (RMSE), and Mean Absolute Error (MAE).

3. Results and Discussion

3.1. Evaluation of the Actual Evapotranspiration

3.1.1. Analysis of the Crop Growth Period in the Testing Area during the Dry Season

Analysis of Landsat-8 satellite imagery collected between 1 November 2018 and 31 May 2019 showed that areas in which rice was cultivated accounted for 0.53 km2 of the testing area during this period. The rice-planting periods in each area differed due to the joint allocation of water. Farmers generally started planting rice in the southern part of the study area, followed by the northern part, with the period from the beginning of rice cultivation to harvest being ~120 days in length. The amplitudes of the NDVI values of each sampling plot were analysed during the growth period in the testing area. Typically, the peak NDVI occurs in the middle of the rice season. Thus, the start date of rice planting can be specified by counting back ~90 days from the date with the highest NDVI value, as shown in the NDVI graph of each area (Figure 3). In addition, the NDVI curve corresponds to Kcb, which was used at each stage of plant growth; this indicated that Kcb was the lowest in the initial stage, increased in the development stage, and peaked in the mid-season stage. The NDVI did not vary significantly for areas covered with longan as a perennial tree, with values ranging from ~0.4–0.6. Supriyasilp et al. [19] interviewed longan farmers and reported that the flowering period in this area ranges from February to May (dry season).
The farmers working on the land in the testing area first started cultivating rice in the A5_Rice sample plot, followed by the plots A2_Rice and A1_Rice (Figure 3). The highest NDVI values for these plots were observed on 26 November 2018, 13 January 2019, and 19 April 2019, respectively. Using these data, we identified the crops’ planting dates and the crop growth periods of each sample plot. The Kcb and specific date of each plant growth stage were used as inputs for the crop section of the SIMDualKc. For the Kcb of rice, the values recommended by Allen et al. [6] and Pereira et al. [23] were used. A recommended value for the Kcb of longan was not available. Therefore, the Kcb value of fruit trees recommended by Allen et al. [6] was used. The SIMDualKc model was then utilised to estimate an ETa value for comparison with the ETa output of the METRIC model based on Landsat-8 satellite data.

3.1.2. Comparison of Actual Evapotranspiration Obtained from the METRIC Model and FAO-56 Method

The spatial distribution of the daily ETa was obtained from the METRIC model based on Landsat-8 images with a resolution of 30 m in the testing area covering the dry season (1 October 2018 to 31 May 2019; Figure 4). The daily ETa exhibited a significant spatial distribution and temporal changes across various land-use types and sample plots. The spatial distribution and temporal changes in ETa during the dry season showed that the ETa started to increase in February and reached its peak in April, corresponding to the crop growth period in irrigated agricultural areas. The spatial distribution and temporal changes in the ETr showed the same patterns as those of the ETa. Figure 4 shows the ETa estimate obtained with METRIC, indicating that the ETa was highest on 19 April 2019. Further, an average ETa > 7 mm/day was observed in the A1 and A2 areas when rice was in its mid-season stage, and in the A4 and A6 areas, which contain longan orchards.
The spatiotemporal distribution obtained for the testing area is consistent with that reported by Supriyasilp et al. [19]: longan is mainly planted in A4 and A6. The amount of water required by longan trees is very high during March and April, i.e., the middle of the longan season. Farmers near the Ping River pump water from the river and directly supply it to the longan area.
Five sample plots in the testing area, including rice and longan plots, were used to evaluate the ETa values obtained with the METRIC model relative to those obtained with the FAO-56 dual-crop coefficient method. Correlations between the daily ETa from METRIC and SIMDualKc for each crop and all sample plots at the satellite acquisition date were evaluated (Figure 5). During the dry season in 2019, 34 observations (N = 34) were used to analyse the correlation between the ETa from the METRIC model and that obtained from the SIMDualKc model (Figure 5). The data plots of the linear regression indicate that the ETa estimates obtained from the METRIC model were strongly correlated with those obtained from SIMDualKc. The ETa values based on METRIC and SIMDualKc were in good agreement. For the rice sample plots, the R2, RMSE, and MAE values were 0.827, 0.827, and 0.711 mm d−1, respectively. For the longan sample plots, the R2, RMSE, and MAE values were 0.847, 0.631, and 0.455 mm d−1, respectively. A comparison of the R2, RMSE, and MAE values between the rice and longan plots showed that the accuracy obtained for the longan plots was higher than that of the rice plots. The overall R2, RMSE, and MAE values of all sample plots (including rice and longan) were 0.830, 0.730, and 0.575 mm d−1, respectively.
The cumulative ETa values derived from METRIC and SIMDualKc for each sample plot covering the dry season in 2019 were used to assess the accuracy of the seasonal ETa (Table 5). Generally, rice was ready to be harvested 120 days after planting. Furthermore, the planting dates of the rice plots in the A2_Rice and A5_Rice plots did not significantly differ (Table 5). Both rice plots were planted in October and harvested from late January to mid-February. The A1_Rice plot was planted later than the other two plots. It was planted in mid-January and harvested in mid-May, with peak growth in April. We also found that ETr was generally higher during March and May than in other months. Therefore, the cumulative ETa of the rice plot in A1_Rice in the dry season was higher than that of the rice plots in A2_Rice and A5_Rice.
The difference between the cumulative ETa derived from SIMDualKc and that from the METRIC model ranged in magnitude from 0.93–3.57% for the rice and 3.08–7.99% for the longan sample plots. In all plots, the cumulative ETa derived from METRIC was lower (4.28%) than that from SIMDualKc.
Tasumi [16] estimated the ET of the western Urmia Lake Basin using METRIC with Landsat satellite images from 2014–2016 for irrigated apple orchards, grape fields, and non-irrigated bare soil. The ET estimates obtained from METRIC (ETMETRIC) and the FAO-56 single-crop coefficient method (ETFAO56) at the satellite image acquisition dates for the irrigated agricultural plots were in good agreement. The ETFAO56 and ETMETRIC differences for apple-planted, grape-planted, and bare soils were indicated by the RMSE, MAE, and Mean Bias Error (MBE), with values of 0.73, 0.84, and 0.68 mm d−1 (positive numbers represent an overestimation); 0.57, 0.72, and 0.59 mm d−1; and 0.20, 0.16, and 0.36 mm d−1, respectively. Overall, the ET estimate from ETMETRIC was greater than that related to ETFAO56. The error range indicated that wet apple fields had a high relative accuracy, while dry bare soil surfaces had a low relative accuracy. Similarly, Paço et al. [15] compared the ET at the dates of image acquisition derived from the METRIC model (ETMETRIC) and FAO-56 dual-crop coefficient method with that based on the SIMDualKc model (ETsim). They reported that ETMETRIC and ETsim had similar patterns. The overall MBE, defined as (ETMETRICETsim)/ETsim), between the ET derived from METRIC and that from SIMDualKc was 18% (the ET from METRIC was greater than the ET from SIMDualKc), corresponding to a mean absolute difference of 0.4 mm/d. Our results agree with previous reports [15,16], apart from the fact that the ET derived from METRIC was smaller than the ET derived from SIMDualKc. It is possible that ETa was underestimated due to a lack of vegetation cover during the initial stage of crop growth [25,26], especially in the rice plots used in this study. In addition, METRIC uses a constant ETrF to estimate the daily ETa; as a result, water stress conditions may be underestimated. [7]. Furthermore, an underestimation of ETa might have occurred due to an underestimation of the daily average Rn computed by empirical equations when upscaling the instantaneous ETa to daily ETa [26]. In addition, satellite pixels over small, cultivated fields may cause vegetation types and densities to mix, resulting in mixed surface temperature signals [7]. Therefore, it is difficult to interpret ET retrievals.

3.2. Spatiotemporal Distribution of the Actual Evapotranspiration for Various Land Use Types

3.2.1. Temporal Variation of the Actual Evapotranspiration

The ET variation and trends in northern Thailand are affected by rainfall and weather conditions. The average rainfall and humidity are lower during the dry season. Each year, the air temperature decreases to its minimum in December and increases to its maximum in April. The average ETr from the AWS in the study area depends on the weather conditions. The average daily ETa was obtained from the METRIC model for nine land-use types between October and May during the dry seasons in 2018–2020 (Figure 6). The average ETr decreased from October–December, with a minimum of 2.88 mm. It gradually increased in January, started to rapidly rise in February, and reached its peak in April (6.37 mm). The ETa variation—mainly influenced by ETr—for waterbodies was relatively consistent with the average trend in ETr. The ETa variation in forestland between October and February was consistent with the trend of ETr, with the lowest ETa in December (4.10 mm), when ETa was greater than ETr. From March–May, the ETa/ETr ratio decreased, and ETa was lower than ETr. The ETa of miscellaneous land (mainly pastures in the study area) was related to the ETa of natural vegetation, such as forests, but was smaller than that of forests. Perennial trees naturally grow without water management during the dry season and had slightly different monthly ETa variations ranging from 2.55–3.26 mm. Field crops were primarily cultivated in non-irrigated areas. Field crops were planted during the rainy season and harvested between November and December. Their ETa variations ranged from 3.07–3.68 mm, with slightly differing monthly ETa values. Paddy fields in the study area were planted and harvested differently depending on the type of water management in the irrigated area. The ETa for paddy fields was relatively consistent with the average trend in ETr but lower than ETr. The lowest and highest ETa values were observed in December (1.67 mm) and April (2.78 mm), respectively. The ETa variation of horticulture between October and January was consistent with the trend in ETr. From February–May, the ratio of the ETa to the ETr of horticulture decreased, with the ETa being lower than the ETr. Orchards showed an ETa pattern similar to that of perennial trees, with slightly different monthly ETa variations. The lowest and highest ETa estimates were obtained in November (2.6 mm) and May (3.13 mm), respectively. Finally, urban and built-up land showed the lowest ETa values, with slightly different monthly ETa variations. The average accumulated seasonal ETa values for forestland, waterbodies, miscellaneous land, horticulture, field crops, perennial trees, orchards, paddy fields, and urban and built-up land over the dry seasons of 2018–2020 were 1,145.15, 1,067.65, 945.48, 935.97, 807.45, 664.76, 655.79, 547.29, and 511.27 mm, respectively (Table 6).
Our results show that the average daily ETa of waterbodies was similar to the ETr. Additionally, the ETa values of forestland and waterbodies were higher than those of irrigated and non-irrigated agricultural areas, including rice, longan, and other land-use types, while the ETa of urban areas was the lowest. These results are consistent with a previous report [27] which showed that forests and wetlands have higher average ETa values than other land-use types, and that developed areas have the lowest average ETa. Thus, ETa appears to be determined by land use type, land cover, and water availability.
The spatiotemporal variation in ETa depended on several factors controlling the ET rate, including weather, vegetation, soil, waterbodies, and management [17]. The ET trends of each land-use type were associated with their biophysical properties and factors controlling the ET. For example, forest areas are resistant to water stress in the regular dry season because the water storage in the soil’s top 3 m layer is adequate for maintaining a stable ET rate [28]. Furthermore, the peak rate of ET was influenced by the height of vegetation. Forests have greater stomatal control than agricultural plants, which leads to higher ET values than other land-use types [17]. Similarly, urban and built-up land represent developed areas with lower ETa values due to lower soil moisture availability [27].

3.2.2. Spatial Distribution of the Actual Evapotranspiration

The spatial distribution of the seasonal ETa was calculated from accumulated time series of daily ETa in the study area during the dry seasons of 2018, 2019, and 2020 (Figure 7). The characteristics of the ETa distribution were consistent with variations in the climate, vegetation, and anthropogenic activities such as water resource management in the study area. The spatial distribution and variation of ET can be used as indicators of vegetation variation and water consumption [29]. Additionally, ETa estimated from the METRIC model can capture crop water stress conditions [16].
Within the study area, area A is an irrigated region belonging to two irrigation projects: Mae Tang and Mae Faek-Mae Ngat Somboon Chon (Figure 7). Area A is mainly covered by paddy fields; the ETa in the cultivated area in the dry season is shown as the blue-shaded area in Figure 7. Areas B and D belong to the Mae Kuang Udom Thara irrigation project and are examples of upstream and downstream areas within the same irrigation project. Most rice-cultivation areas with adequate water allocation are located upstream of the irrigation project area (area B) and are shaded blue for all three years. Area D consists of primarily urban and built-up land in Lamphun Province. The paddy fields in area D were not cultivated during the dry season (shown as orange-shaded areas). In area C, which contains longan orchards, blue shading can be observed throughout the whole area across all three years. As mentioned before, longan is an economically valuable fruit; thus, sufficient water is required for its maintenance.

4. Conclusions

In this study, satellite-based surface energy balance data derived from Landsat-8 imagery and the METRIC model were applied to estimate the spatiotemporal distribution of ETa at a field scale, and to compare the ETa obtained from the METRIC model with that derived from the FAO-56 dual-crop coefficient method using SIMDualKc.
The daily ETa exhibited substantial spatial variation and temporal changes across various land-use types and sample plots in the testing area throughout the dry season. The accuracy of the ETa values obtained from the METRIC model and FAO-56 dual-crop coefficient method was determined based on statistical analyses performed on the satellite acquisition dates and the cumulative ETa of all sample plots. The data plots indicated that the ETa estimates derived from the METRIC model were similar to those derived from SIMDualKc, and we observed a strong linear correlation between the estimates calculated using the two methods. Overall, the accuracy assessment for all sample plots on the satellite acquisition date yielded R2, RMSE, and MAE values of 0.830, 0.730, and 0.575 mm d−1, respectively, and the ETa accuracy obtained for the longan plots was higher than that for the rice plots. The differences between the cumulative ETa derived from SIMDualKc and that derived from the METRIC model ranged in magnitude from 0.93%–3.57% and 3.08%–7.99% for the rice and longan sample plots, respectively. For all plots, the cumulative ETa derived from METRIC was lower than that derived from SIMDualKc. Analysis of the ETa variation for various land-use types showed that the average daily ETa for waterbodies is similar to the ETr. The ETa values obtained for forestland and waterbodies were higher than those for irrigated and non-irrigated agricultural areas, including rice, longan, and other land-use types, while the ETa for urban areas was the lowest. Overall, the observed spatial distributions and temporal changes of the seasonal ETa during the dry seasons of 2018, 2019, and 2020 correlated with variations in the climate, vegetation, and anthropogenic activities such as water management in the study area each year.
The results of this study indicate that the METRIC model is an efficient tool for estimating ETa during crop growth in the dry season, as well as the spatiotemporal distribution of the ETa for various land-use types at a field scale, in order to facilitate water resource management. Furthermore, the METRIC model can be applied to estimate the spatiotemporal distribution of ETa at field to regional scales, making it suitable for application to water resource management under various environmental conditions and for the improvement of the water use efficiency over a large area. In the future, ETa studies should be expanded to other areas in Southeast Asia, as this region is an essential source of agricultural production.

Author Contributions

Conceptualisation, T.S. (Teerawat Suwanlertcharoen), T.S. (Thanaporn Supriyasilp) and K.P.; Formal analysis, T.S. (Teerawat Suwanlertcharoen); Investigation, T.S. (Teerawat Suwanlertcharoen); Methodology, T.S. (Teerawat Suwanlertcharoen); Supervision, T.C., T.S. (Thanaporn Supriyasilp) and K.P.; Visualisation, T.S. (Teerawat Suwanlertcharoen); Writing—original draft, T.S. (Teerawat Suwanlertcharoen); Writing—review and editing, T.S. (Thanaporn Supriyasilp). All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data will be provided upon request.

Acknowledgments

We thank the Thai Meteorological Department for providing the climate and rainfall data used in this study. We also thank the Land Development Department for providing the land-use data used in this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Takeshima, H.; Joshi, P.K. Overview of the Agricultural Modernization in Southeast Asia; IFPRI Discussion Paper 1819; International Food Policy Research Institute (IFPRI): Washington, DC, USA, 2019. [Google Scholar] [CrossRef]
  2. Association of Southeast Asian Nations. ASEAN Statistical Yearbook 2020. 2021. Available online: https://aseandse.org/wp-content/uploads/2021/02/ASYB_2020.pdf (accessed on 30 October 2021).
  3. Altendorf, S. Minor Tropical Fruits (Mainstreaming a Niche Market). Food Outlook 2018, 8, 67–75. Available online: http://www.fao.org/fileadmin/templates/est/COMM_MARKETS_MONITORING/Tropical_Fruits/Documents/Minor_Tropical_Fruits_FoodOutlook_1_2018.pdf (accessed on 21 March 2022).
  4. Mancosu, N.; Snyder, R.L.; Kyriakakis, G.; Spano, D. Water scarcity and future challenges for food production. Water 2015, 7, 975–992. [Google Scholar] [CrossRef]
  5. Spreer, W.; Schulze, K.; Ongprasert, S.; Wiriya-Alongkorn, W.; Müller, J. Mango and longan production in northern Thailand: The role of water saving irrigation and water stress monitoring. In Sustainable Land Use and Rural Development in Southeast Asia: Innovations and Policies for Mountainous Areas; Fröhlich, H.L., Schreinemachers, P., Stahr, K., Clemens, G., Eds.; Springer Environmental Science and Engineering; Springer: Berlin, Germany, 2013; pp. 215–228. [Google Scholar] [CrossRef] [Green Version]
  6. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop evapotranspiration-Guidelines for computing crop water requirements-FAO Irrigation and drainage paper 56; FAO: Rome, Italy, 1998; Volume 300, p. D05109.
  7. Allen, R.G.; Pereira, L.S.; Howell, T.A.; Jensen, M.E. Evapotranspiration information reporting: I. Factors governing measurement accuracy. Agric. Water Manag 2011, 98, 899–920. [Google Scholar] [CrossRef] [Green Version]
  8. Bastiaanssen, W.G.M.; Menenti, M.; Feddes, R.A.; Holtslag, A.A.M. A remote sensing Surface Energy Balance Algorithm for Land (SEBAL). J. Hydrol. 1998, 212–213, 198–212. [Google Scholar] [CrossRef]
  9. Allen, R.G.; Tasumi, M.; Trezza, R. Satellite-based energy balance for Mapping Evapotranspiration with Internalized Calibration (METRIC)—Model. J. Irrig. Drain. Eng. 2007, 133, 380–394. [Google Scholar] [CrossRef]
  10. Allen, R.G.; Irmak, A.; Trezza, R.; Hendrickx, J.M.H.; Bastiaanssen, W.; Kjaersgaard, J. Satellite-based ET estimation in agriculture using SEBAL and METRIC. Hydrol. Process 2011, 25, 4011–4027. [Google Scholar] [CrossRef]
  11. Zheng, C.; Jia, L.; Hu, G.; Lu, J. Earth observations-based evapotranspiration in northeastern Thailand. Remote Sens. 2019, 11, 138. [Google Scholar] [CrossRef] [Green Version]
  12. Sriwongsitanon, N.; Suwawong, T.; Thianpopirug, S.; Williams, J.; Jia, J.; Bastiaanssen, W. Validation of seven global remotely sensed ET products across Thailand using water balance measurements and land use classifications. J. Hydrol. Reg. Stud. 2020, 30, 100709. [Google Scholar] [CrossRef]
  13. Kosa, P.; Pongput, K.; Nontananandh, S.; Khobklag, P. Spatial and temporal distribution of irrigation requirements in the Chao Phraya River Basin, Thailand. Kasetsart J. Nat. Sci. 2008, 42, 758–766. [Google Scholar]
  14. Maskong, H.; Kosa, P.; Jothityangkoon, C. Estimation of evapotranspiration in Lam Ta Kong Basin using Surface Energy Balance Algorithm for Land (SEBAL) Model. In Proceedings of the THA 2015 International Conference on Climate Change and Water & Environment Management in Monsoon Asia, Bangkok, Thailand, 28–30 January 2015. [Google Scholar]
  15. Paço, T.A.; Pôças, I.; Cunha, M.; Silvestre, J.C.; Santos, F.L.; Paredes, P.; Pereira, L.S. Evapotranspiration and crop coefficients for a super intensive olive orchard. An Application of SIMDualKc and METRIC models using ground and satellite observations. J. Hydrol. 2014, 519, 2067–2080. [Google Scholar] [CrossRef] [Green Version]
  16. Tasumi, M. Estimating evapotranspiration using METRIC Model and Landsat data for better understandings of regional hydrology in the Western Urmia lake basin. Agric. Water Manag 2019, 226, 105805. [Google Scholar] [CrossRef]
  17. Jensen, M.E.; Allen, R.G. Evaporation, Evapotranspiration, and Irrigation Water Requirements, 2nd ed.; American Society of Civil Engineers: Reston, VA, USA, 2016. [Google Scholar]
  18. HAII. Data Warehouse System Development Project of 25 Basin and Flood and Drought Modeling: Ping Basin. [in Thai]. 2012. Available online: https://tiwrm.hii.or.th/web/attachments/25basins/06-ping.pdf (accessed on 31 October 2021).
  19. Supriyasilp, T.; Pongput, K.; Boonyanupong, S.; Suwanlertcharoen, T. Enhanced water management for Muang Fai irrigation systems through remote sensing and SWOT analysis. Water Resour Manag. 2021, 35, 263–277. [Google Scholar] [CrossRef]
  20. Supriyasilp, T.; Pongput, K. Water Balance Model as a tool for building participation and joint water use promotion at the river basin level. Irrig. Drain 2021, 70, 254–268. [Google Scholar] [CrossRef]
  21. Tasumi, M. Progress in Operational Estimation of Regional Evapotranspiration Using Satellite Imagery. PhD Dissertation, University of Idaho, Moscow, Russia, 2003. [Google Scholar]
  22. Allen, R.G.; Burnett, B.; Kramber, W.; Huntington, J.; Kjaersgaard, J.; Kilic, A.; Kelly, C.; Trezza, R. Automated calibration of the METRIC-Landsat Evapotranspiration Process. J. Am. Water Resour. Assoc. 2013, 49, 563–576. [Google Scholar] [CrossRef]
  23. Pereira, L.S.; Paredes, P.; López-Urrea, R.; Hunsaker, D.J.; Mota, M.; Mohammadi Shad, Z. Standard single and basal crop coefficients for vegetable crops, an Update of FAO56 Crop Water Requirements Approach. Agric. Water Manag. 2021, 243, 106196. [Google Scholar] [CrossRef]
  24. Rosa, R.D.; Paredes, P.; Rodrigues, G.C.; Fernando, R.M.; Alves, I.; Pereira, L.S.; Allen, R.G. Implementing the dual crop coefficient approach in interactive software: 2. Model Testing. Agricultural Water Manag. 2012, 103, 62–77. [Google Scholar] [CrossRef]
  25. Mkhwanazi, M.; Chávez, J.L.; Andales, A. SEBAL-a: A remote sensing ET algorithm that accounts for advection with limited data. Part I: Development and Validation. Remote Sens. 2015, 7, 15046–15067. [Google Scholar] [CrossRef] [Green Version]
  26. Xue, J.; Bali, K.M.; Light, S.; Hessels, T.; Kisekka, I. 2020. Evaluation of remote sensing-based evapotranspiration models against surface renewal in almonds, tomatoes and maize. Agric. Water Manag. 2020, 238, 106228. [Google Scholar] [CrossRef]
  27. Liu, W.; Hong, Y.; Khan, S.I.; Huang, M.; Vieux, M.; Caliskan, S.; Grout, T. Actual evapotranspiration estimation for different land use and land cover in urban regions using Landsat 5 data. J. Appl. Remote Sens. 2010, 4, 041873. [Google Scholar] [CrossRef]
  28. Khand, K.; Numata, I.; Kjaersgaard, J.; Vourlitis, G.L. Dry season evapotranspiration dynamics over human-impacted landscapes in the Southern Amazon using the Landsat-based METRIC Model. Remote Sens. 2017, 9, 706. [Google Scholar] [CrossRef] [Green Version]
  29. Liu, J.; Xiong, Y.; Tian, J.; Tan, Z. Spatiotemporal changes in evapotranspiration from an overexploited water resources basin in arid northern China and their implications for ecosystem management. Sustainability 2019, 11, 445. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Study and testing areas: (a) general location of the study area in Thailand; (b) location of the study areas in the Chiang Mai and Lamphun provinces, which are situated in the plains of the upper northern region of Thailand along the Ping River; (c) testing area, which forms part of the Muang Luang Irrigation Project located in the Chom Thong District, Chiang Mai Province, downstream of the Muang Luang Weir.
Figure 1. Study and testing areas: (a) general location of the study area in Thailand; (b) location of the study areas in the Chiang Mai and Lamphun provinces, which are situated in the plains of the upper northern region of Thailand along the Ping River; (c) testing area, which forms part of the Muang Luang Irrigation Project located in the Chom Thong District, Chiang Mai Province, downstream of the Muang Luang Weir.
Water 15 00450 g001
Figure 2. Frameworks and workflow for estimating actual evapotranspiration (ETa) during each crop growth period using the METRIC model and Landsat-8 satellite data; evaluating the ETa obtained from the METRIC model relative to the ETa calculated using the FAO-56 dual-crop coefficient method; and studying the spatiotemporal distribution of the ETa for various land-use types to improve water resource management in northern Thailand.
Figure 2. Frameworks and workflow for estimating actual evapotranspiration (ETa) during each crop growth period using the METRIC model and Landsat-8 satellite data; evaluating the ETa obtained from the METRIC model relative to the ETa calculated using the FAO-56 dual-crop coefficient method; and studying the spatiotemporal distribution of the ETa for various land-use types to improve water resource management in northern Thailand.
Water 15 00450 g002
Figure 3. Crop phenology for plots in the testing area located in the Chom Thong District, Chiang Mai Province, downstream of the Muang Luang Weir in Thailand based on the Normalised Difference Vegetation Index (NDVI) in the dry season.
Figure 3. Crop phenology for plots in the testing area located in the Chom Thong District, Chiang Mai Province, downstream of the Muang Luang Weir in Thailand based on the Normalised Difference Vegetation Index (NDVI) in the dry season.
Water 15 00450 g003
Figure 4. Daily actual evapotranspiration (ETa) estimates for the testing area in the Chom Thong District, Chiang Mai Province, downstream of the Muang Luang Weir, Thailand, based on the mapping evapotranspiration at high resolution with internalized calibration (METRIC) model.
Figure 4. Daily actual evapotranspiration (ETa) estimates for the testing area in the Chom Thong District, Chiang Mai Province, downstream of the Muang Luang Weir, Thailand, based on the mapping evapotranspiration at high resolution with internalized calibration (METRIC) model.
Water 15 00450 g004
Figure 5. Correlation between the actual evapotranspiration (ETa) values derived from the mapping evapotranspiration at high resolution with internalized calibration (METRIC) model and SIMDualKc at the satellite acquisition date for all sample plots from the testing area in the Chom Thong District, Chiang Mai Province, downstream of the Muang Luang Weir, Thailand.
Figure 5. Correlation between the actual evapotranspiration (ETa) values derived from the mapping evapotranspiration at high resolution with internalized calibration (METRIC) model and SIMDualKc at the satellite acquisition date for all sample plots from the testing area in the Chom Thong District, Chiang Mai Province, downstream of the Muang Luang Weir, Thailand.
Water 15 00450 g005
Figure 6. Average daily actual evapotranspiration (ETa) and standardised alfalfa reference ET equation (ETr) for nine land-use types in the study area in the Chiang Mai and Lamphun provinces, which are located on the plains of the upper northern region of Thailand along the Ping River, during the dry seasons of 2018–2020.
Figure 6. Average daily actual evapotranspiration (ETa) and standardised alfalfa reference ET equation (ETr) for nine land-use types in the study area in the Chiang Mai and Lamphun provinces, which are located on the plains of the upper northern region of Thailand along the Ping River, during the dry seasons of 2018–2020.
Water 15 00450 g006
Figure 7. Spatial distribution of the seasonal actual evapotranspiration (ETa) in the study area located in the Chiang Mai and Lamphun provinces, which is located on the plains of the upper northern region of Thailand along the Ping River, during the dry seasons of 2018–2020.
Figure 7. Spatial distribution of the seasonal actual evapotranspiration (ETa) in the study area located in the Chiang Mai and Lamphun provinces, which is located on the plains of the upper northern region of Thailand along the Ping River, during the dry seasons of 2018–2020.
Water 15 00450 g007
Table 1. Landsat-8 OLI/TIRS images (Path 131/Row 47) for the study area.
Table 1. Landsat-8 OLI/TIRS images (Path 131/Row 47) for the study area.
Dry Season in Year 2018Dry Season in Year 2019Dry Season in Year 2020
DateDay of Year (DOY)Time of
Acquisition (hh:mm:ss) [UTC+7]
DateDay of Year (DOY)Time of
Acquisition (hh:mm:ss) [UTC+7]
DateDay of Year (DOY)Time of
Acquisition (hh:mm:ss) [UTC+7]
22 Oct 201729510:49:149 Oct 201828210:48:4412 Oct 201928510:49:19
9 Dec 201734310:49:0425 Oct 201829810:48:4928 Oct 201930110:49:19
25 Dec 201735910:49:0726 Nov 201833010:48:5113 Nov 201931710:49:17
10 Jan 20181010:49:0128 Dec 201836210:48:4829 Nov 201933310:49:15
26 Jan 20182610:48:5213 Jan 20191310:48:4615 Dec 201934910:49:13
11 Feb 20184210:48:4529 Jan 20192910:48:4231 Dec 201936510:49:09
27 Feb 20185810:48:3914 Feb 20194510:48:4016 Jan 20201610:49:06
15 Mar 20187410:48:302 Mar 20196110:48:361 Feb 20203210:49:01
31 Mar 20189010:48:2318 Mar 20197710:48:3117 Feb 20204810:48:57
16 Apr 201810610:48:143 Apr 20199310:48:284 Mar 20206410:48:52
18 May 201813810:47:5419 Apr 201910910:48:2220 Mar 20208010:48:45
3 Jun 201815410:47:435 May 201912510:48:235 Apr 20209610:48:35
21 May 201914110:48:3321 Apr 202011210:48:29
6 Jun 201915710:48:427 May 202012810:48:20
23 May 202014410:48:22
Table 2. Summary of soil properties for each plot in testing area.
Table 2. Summary of soil properties for each plot in testing area.
Plot NameSoil TextureSand (%)Silt (%)Clay (%)Rock (%)Location
A1_RiceLoam49.333.7177.8718°24′48.1″ N, 98°40′44.1″ E
A2_RiceSandy Loam64191718.2418°24′16.2″ N, 98°41′01.1″ E
A4_LonganSandy Clay Loam45.82529.215.418°23′48.2″ N, 98°41′24.3″ E
A5_RiceSandy Loam6422.813.20.1418°23′51.4″ N, 98°41′11.0″ E
A6_LonganSandy Loam67.920.611.50.7718°23′17.4″ N, 98°41′05.4″ E
Table 3. Datasets used in this study.
Table 3. Datasets used in this study.
DatasetDescriptionSource
Satellite imageryLandsat-8 Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS) images with 30 m spatial resolution were collected from level-1 products with terrain precision correction (L1TP)USGS Earth Explorer (https://earthexplorer.usgs.gov/) (accessed on 3 January 2013)
Meteorological dataThe Automatic Weather System (AWS) from the Lamphun station records meteorological data every 10 min, including the precipitation, air temperature, wind speed, relative humidity, and atmospheric pressure. In addition, precipitation data were collected at the Chom Thong station.Thailand Meteorological Department
Digital Elevation Model (DEM)Shuttle Radar Topography Mission (SRTM) data with 30 m spatial resolutionUSGS Earth Resources Observation and Science (EROS) Centre
Land UseGIS data regarding land-use classification at 1:25,000 scaleLand Development Department, Thailand
Irrigated AreaGIS data of irrigated area at 1:50,000 scaleRoyal Irrigation Department, Thailand
Soil SamplingSoil properties for each plot in testing areaSoil samples were taken from sample plots in testing area
Water Use and Water MeasurementWater supplied from various sources and the flow data observed in the dry season were measured from 21 March to 15 May 2019 in testing areaAccording to Supriyasilp and Pongput [20]
Table 4. Software used in this study.
Table 4. Software used in this study.
SoftwareCapabilities
ArcGIS 10.3.1Image processing, spatial analysis, spatial statistics, mapping, visualisation, and exploration for estimation of satellite-based ETa.
SIMDualKcEstimation of ETa using the FAO-56 dual-crop coefficient method.
REF-ETCalculation of the reference ET.
DATimeSFill data gaps and generate the time series of daily ETrF and NDVI image using the cubic spline interpolation method.
Table 5. Cumulative actual evapotranspiration (ETa) values derived from the mapping evapotranspiration at high resolution with internalized calibration (METRIC) model and SIMDualKc for each sample plot in the testing area in the Chom Thong District, Chiang Mai Province, downstream of the Muang Luang Weir in Thailand, covering the dry season in 2019.
Table 5. Cumulative actual evapotranspiration (ETa) values derived from the mapping evapotranspiration at high resolution with internalized calibration (METRIC) model and SIMDualKc for each sample plot in the testing area in the Chom Thong District, Chiang Mai Province, downstream of the Muang Luang Weir in Thailand, covering the dry season in 2019.
Plot NameSeasonal DurationCumulative ETa (SIMDualKc)Cumulative ETa (METRIC)Percentage Difference (%)
A1_Rice11/01/2019 to 11/05/2019549.96530.683.57
A2_Rice17/10/2018 to 14/02/2019328.14318.852.87
A5_Rice02/10/2018 to 29/01/2019313.68310.780.93
A4_Longan01/11/2018 to 31/05/2019817.99755.147.99
A6_Longan01/11/2018 to 31/05/2019872.40845.933.08
Sum2882.172761.384.28
Table 6. Accumulated monthly and seasonal ETa values for various land-use types and ETr in the study area during the dry seasons of 2018, 2019, and 2020.
Table 6. Accumulated monthly and seasonal ETa values for various land-use types and ETr in the study area during the dry seasons of 2018, 2019, and 2020.
DateETa (mm)ETr
(mm)
Paddy FieldField CropMiscellaneous LandUrban and Built-Up LandPerennial TreeOrchardHorticultureForestlandWaterbody
2017/10/01 to 2017/10/3182.16100.69110.9459.4588.6880.01107.76132.3690.09110.09
2017/11/01 to 2017/11/3058.7682.0194.4750.5667.5262.1095.66113.7466.1493.51
2017/12/01 to 2017/12/3155.0685.9998.1862.8776.2074.1797.97117.8386.4887.71
2018/01/01 to 2018/01/3157.9394.99114.7065.3586.5084.99111.73141.24120.4994.61
2018/02/01 to 2018/02/2860.3781.20105.6858.7875.2175.29104.76130.49125.54115.45
2018/03/01 to 2018/03/3158.1566.96102.3644.7059.6659.46103.86128.14134.18170.43
2018/04/01 to 2018/04/3055.8665.24101.3943.6564.1661.37100.17135.63136.11170.84
2018/05/01 to 2018/05/3183.67112.58137.4973.51122.52105.04128.70171.64157.21165.49
Sum511.96689.66865.22458.86640.44602.44850.611071.07916.241008.13
2018/10/01 to 2018/10/3169.31106.17110.4046.6281.1866.75104.69135.4082.59113.48
2018/11/01 to 2018/11/3073.95120.14118.2873.29110.6387.86110.60138.60108.1998.85
2018/12/01 to 2018/12/3159.0897.38103.0566.0085.4278.7599.61118.51107.2987.05
2019/01/01 to 2019/01/3154.9097.24109.9756.1475.0076.20109.45132.00104.1494.61
2019/02/01 to 2019/02/2875.8593.41116.0864.7573.2682.40115.02136.79139.07131.77
2019/03/01 to 2019/03/3185.9785.84122.9261.8664.4078.04126.09150.08180.29187.3
2019/04/01 to 2019/04/30130.36112.32158.1198.47102.09115.45157.23192.03243.42204.21
2019/05/01 to 2019/05/3162.7277.75112.5954.6273.6774.14115.71142.46135.50195.71
Sum612.12790.25951.39521.76665.64659.61938.391145.861100.511112.98
2019/10/01 to 2019/10/3162.87135.57135.6558.8599.8186.35132.89162.87109.69122.6
2019/11/01 to 2019/11/3030.8996.7898.6733.0259.6253.43100.11123.5380.28102.56
2019/12/01 to 2019/12/3140.72111.45115.9153.5975.1176.43113.05144.65107.2893.48
2020/01/01 to 2020/01/3164.52112.37122.1473.6781.7489.80118.96141.63144.87103.02
2020/02/01 to 2020/02/2973.67109.83123.2478.2379.8490.00123.06140.18165.74135.31
2020/03/01 to 2020/03/31102.99146.18163.38104.62113.16118.19168.52186.89236.65202.05
2020/04/01 to 2020/04/3063.6098.44114.3169.9872.1079.62116.91139.14173.38198.4
2020/05/01 to 2020/05/3178.55131.83146.5381.23106.80111.51145.43179.62168.32198.77
Sum517.80942.441,019.84553.19688.18705.331018.921218.521186.201156.19
Average seasonal of ETa
(2018–2020)
547.29807.45945.48511.27664.76655.79935.971145.151067.651092.43
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

Suwanlertcharoen, T.; Chaturabul, T.; Supriyasilp, T.; Pongput, K. Estimation of Actual Evapotranspiration Using Satellite-Based Surface Energy Balance Derived from Landsat Imagery in Northern Thailand. Water 2023, 15, 450. https://doi.org/10.3390/w15030450

AMA Style

Suwanlertcharoen T, Chaturabul T, Supriyasilp T, Pongput K. Estimation of Actual Evapotranspiration Using Satellite-Based Surface Energy Balance Derived from Landsat Imagery in Northern Thailand. Water. 2023; 15(3):450. https://doi.org/10.3390/w15030450

Chicago/Turabian Style

Suwanlertcharoen, Teerawat, Thodsapol Chaturabul, Thanaporn Supriyasilp, and Kobkiat Pongput. 2023. "Estimation of Actual Evapotranspiration Using Satellite-Based Surface Energy Balance Derived from Landsat Imagery in Northern Thailand" Water 15, no. 3: 450. https://doi.org/10.3390/w15030450

APA Style

Suwanlertcharoen, T., Chaturabul, T., Supriyasilp, T., & Pongput, K. (2023). Estimation of Actual Evapotranspiration Using Satellite-Based Surface Energy Balance Derived from Landsat Imagery in Northern Thailand. Water, 15(3), 450. https://doi.org/10.3390/w15030450

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