Improvement of reference crop evapotranspiration estimates using limited data for the Brazilian Cerrado

ABSTRACT Evapotranspiration (ET) is a key component of the hydrological cycle. Therefore, adequately estimating it is crucial to improving water resource planning and management.One of the most affordable methods of estimating ET is first to estimate reference crop evapotranspiration (ETo) and later associate it to crop and soil coefficients. The FAO Kc-ETo approach can be used only when ETo is computed with the FAO Penman-Monteith equation. However, low data availability may restict the equations used to estimate ETo. In this study, we assess and calibrate common methods used to estimate ETo under such conditions of limited data availability. Based on the annual calibration, the Makkink (NSE = 0.85) outperformed the Priestley-Taylor (NSE = 0.73), Hargreaves-Samani (NSE = 0.56), and Penman-Monteith temperature approach (NSE = 0.58). The seasonal calibration of parameters showed no significant improvement to the methods assessed (ΔNSE ≤ 0.01), except for the Priestley-Taylor (ΔNSE = 0.06). The performance of temperature-based equations was particularly limited due to the performance of the equation adopted to estimate global solar radiation. Thus, improving the representation of global solar radiation for limited data availability can also play a key role in improving ETo prediction.


Introduction
In Brazil, disputes over water use seem to be more widespread in the northeastern and southeastern regions, where rainfall is highly variable.This being the case, special attention should be given to the savannah biome (Cerrado), which represents ~24 % of Brazil's territory, and is the most important and the last agricultural frontier in Brazil and one of the few areas with the potential to increase food production sustainably (Klink, 2014;Wendt et al., 2015).
In this context, crop evapotranspiration (ETc) is a crucial component of the hydrological cycle and water balance in the soil.Thus, obtaining more accurate evapotranspiration (ET) estimates is crucial to the development of more efficient water source management, irrigation practices, and trustworthy climate studies.Because of the high costs associated with obtaining direct ET observations, it is commonly estimated by physics-based or empirical equations.ET can be estimated from meteorological data in the case of crops by first calculating a reference crop evapotranspiration (ETo) as defined in the Food and Agriculture Organization Irrigation and Drainage paper 56 (FAO56, Allen et al., 1998).Crop evapotranspiration results from adjusting ETo to water-stress conditions and specific crop coefficients.
The Penman-Monteith (PM) equation parametrized in FAO56 (Allen et al., 1998) has long been used as the standard method for estimating daily ETo, which requires data on solar radiation (or sunshine hours), air temperature, relative humidity, and wind speed.
However, many regions worldwide lack one or more meteorological inputs required by the PM equation (Zanetti et al., 2019).
Several empirical equations have been developed to provide ETo estimates under conditions of limited meteorological data availability, such as the Hargreaves and Samani (1985) (HS), Priestley and Taylor (1972) (PT), and Makkink (1957) (MK) equations.These equations were developed considering specific climatic conditions and may require calibration to adequately estimate ETo in other locations and under different conditions (Gavilán et al., 2006).
Given the importance of both providing more accurate estimates for environmental and hydrological studies, and developing sustainable irrigation management practices, the objectives of the present study were to calibrate new coefficients for the PT, MK, HS, and the PM temperature approach (PMT) methods, on both the annual and seasonal scales for the entire Cerrado region.Although solar radiation data is less commonly available, temperature-radiation-based methods were included in this study to compare with temperature-based methods that might have errors from estimating solar radiation prior to the ETo calculation.

Study area
The Cerrado is Brazilian's second-largest biome, covering 23 % of the territory and extending into 11 states.The geographic coordinates are from 60°28'22" W to 41°16'48" W and from 2°19'48" S to 24°40'55" S (Figure 1), and its altitude ranges from 0 m to 2854 m above sea level (mean altitude: 396 m).Tropical grasslands and savannah characterize vegetation in the Cerrado.The region's climate is predominantly classified as tropical with dry winter (Alvares et al., 2013).Thus, the region is characterized by two welldefined seasons: a wet season from Oct to Apr, when ~85 % to 90 % of the annual rainfall occurs, and a dry season from May to Sept.

Data
from 2000 to 2020 featuring 82 automatic weather stations (AWS) and 50 conventional weather stations (Figure 1) were used to assess reference crop evapotranspiration equations for the Cerrado.Data were acquired from the Brazilian National Institute of Meteorology (INMET) Meteorological Database.The INMET's weather stations were installed and are maintained according to the World Meteorological Organization guide to meteorological instruments and observation methods (WMO, 2008).Although there is no information pertaining to the stations' management of the grass surface, we adopted these as reference sites.The stations are listed in Appendix I together with their geographic coordinates, elevation, and number of observations during the study period.
INMET's CWS records daily data consisted of daily mean, maximum and minimum air temperature (T avg , T max , and T min , °C); average relative humidity (RH, %), wind speed at 10 m above ground (Ws, m s -1 ); and sunshine hours (SH, h).SH and Ws were converted into global solar radiation (Rs, MJ m -2 d -1 ) and wind speed at 2 m above ground (u 2 , m s -1 ), respectively, adopting the FAO56 methodology (Allen et al., 1998).The AWS records hourly data for the same variables but instead of sunshine hours, records solar radiation directly.Although hourly computations of ETo provide more accurate predictions, hourly records from the AWS dataset were summarised into daily values to render CWS and AWS equivalent, e.g., the minimum and maximum daily values from hourly temperatures were denominated as T min and T max .Temperature, relative humidity, and wind speed were averaged, while solar radiation was summed over daily periods.
The time series for all weather stations were carefully screened to eliminate flawed daily observations, e.g., infinite values, negative wind speed values, T min higher than T max , SH or Rs above the maximum physically possible, and RH outside the range of 5 to 100 %.The Rs time series was compared to the clear-sky solar radiation time series according to Allen (1996), and periods with shifts concerning clear-sky solar radiation were removed.Data were also visually screened for homogeneity (Althoff et al., 2020;Xavier et al., 2016).Weather stations with less than ten years of observations were excluded from the analysis.

Reference crop evapotranspiration
This section describes all ETo equations considered in this study.The Penman-Monteith equation using daily input was selected in this study as the reference ETo equation for the evaluation of ETo methods which can be adopted under conditions of limited data availability.The Penman-Monteith equation has been shown to be the most accurate equation for estimating the ETo measured by lysimeters in regions with the same climate classification as the predominant one in the Cerrado (Barros et al., 2009;Mendonça et al., 2003).Thus, we first describe the Penman-Monteith equation (following section), which requires that all meteorological variables described in the previous section (Data) are available.Next, we describe the equations used to estimate ETo under limited data availability, including temperature-based and temperature-radiation-based methods.These equations are presented in a generic form, i.e., with coefficients instead of the originally proposed values.

FAO Penman-Monteith
From a diverse number of models developed in the literature for estimating ETo, the Penman-Monteith (PM) standardized by FAO 56 (Allen et al., 1998) has been widely used as a reference for calibrating a diverse number of methods (Althoff et al., 2018;Feng et al., 2017;Ferreira et al., 2019;Shiri et al., 2014).Penman-Monteith's ETo is calculated as follows: the mean saturation vapor pressure for a day, kPa; e a = actual vapor pressure, kPa; Δ = slope of saturation vapor pressure curve at air temperature (kPa °C-1 ); and, γ = psychrometric constant (kPa °C-1 ).Rn, e s , e a, Δ, and γ were calculated according to FAO56 (Allen et al., 1998), while G can be assumed to equal zero as the soil heat flux beneath the reference grass is very limited on a daily time scale.

Penman-Monteith temperature approach
Since the acquisition of all meteorological data needed to calculate PM can be expensive, it is a widespread practice to work with reduced data.To cope with fewer variables, several alternatives were proposed by Allen et al. (1998).In the Penman-Monteith temperature approach (PMT), we assume that only the temperature daily range (T max and T min ) is available.Alternative methods were used to estimate missing variables, i.e., relative humidity, solar radiation, and wind speed.When relative humidity is unavailable, Allen et al. (1998) recommend computing actual vapor pressure, e a , considering dew point temperature (T dew ) to be equal to T min (Eq.2).Although T min can exceed T dew in semi-arid regions, the opposite can also be true in humid regions (Paredes and Pereira, 2019).The high seasonality in tropical regions may also affect this relationship between wet (Oct-Apr) and dry seasons (May-Sept).It is common to assume T min = T dew , however, a temperature correction factor (a T ) can be used to improve the model.Once ea is calculated using T max , T min , and RH (see Allen et al., 1998), the equations below can be used to derive a T : (2) where a T is the dew temperature correction coefficient (°C).
In the absence of sunshine hours or solar radiation, Allen et al. (1998) suggest that Rs can be estimated using the equation proposed by Hargreaves and Samani (1985), which expresses Rs as a function of temperature range: where k Rs is an empirical radiation adjustment coefficient (°C -0.5 ) and Ra the extraterrestrial solar radiation (MJ m -2 d -1 ).Ra is calculated considering the station latitude and day of the year (Allen et al., 1998).A k Rs equal to 0.16 is suggested for inland sites.
When wind speed data is not available, Allen et al. (1998) recommend using the world average wind speed value of 2.0 m s -1 .An average wind speed for the Cerrado and its wet and dry seasons were also evaluated instead of using the world average for model improvement.Finally, the Penman-Monteith temperature approach (PMT) can be calculated using Eq. ( 1).

Hargreaves-Samani
The Hargreaves-Samani (HS) model (Eq.5) is a temperature-based method also recommended by Allen et al. (1998) for occasions when solar radiation, relative humidity, and wind speed data are missing.The HS equation has shown reliable results with global validity (Allen et al., 1998) and is calculated as follows: where 0.0135 is a factor in the conversion of units, offset an empirical coefficient, and l is the latent heat of vaporization, commonly taken to be 2.45 MJ kg -1 .Hargreaves and Samani (1985) suggest the offset = 17.8.

Priestley-Taylor
The Priestley-Taylor (PT) is a temperature-radiation model (Eq.6) that was originally proposed for estimating evaporation from saturated land and open water sites in advection-free scenarios.As the method does not consider vapor deficit and aerodynamic influence on evapotranspiration, it does not require wind speed and relative humidity data.The PT equation can be described as follows: where α = empirical coefficient, and α is recommended to be equal to 1.26 by Priestley and Taylor (1972).

Makkink
The Makkink (MK) model is also a temperatureradiation-based method (Eq.7), whose coefficients have already been parameterized by many authors.However, its original form is written as follows: where c 1 and c 2 = empirical coefficients.The values suggested for the empirical coefficients c 1 and c 2 are 0.61 and -0.12 (Makkink, 1957), respectively.

Model improvement
The meteorological dataset was divided into two subsets: one set ranging from 2000 to 2014 (approximately 70 %), termed the calibration set, and a second from 2015 to 2020 (approximately 30 %), termed the validation set.The equations were assessed for the validation set, considering both the original empirical coefficients and coefficients calibrated for the calibration set.This study differs from previous studies in that it considers the whole Cerrado biome and uses a single coefficient parametrization for its entirety.The empirical coefficients were calibrated for the Cerrado using the calibration set considering the (i) whole period (annual) and (ii) wet and dry seasons (seasonal).The empirical coefficients were calibrated by minimizing the sum of squares with the Levenberg-Marquardt algorithm (Moré, 1978).The performance of the equations using their original and calibrated coefficients (annual and seasonal) was evaluated for the validation set considering the (i) entire validation set, (ii) monthly periods, and (iii) spatially (station-wise).Figure 2 presents a flowchart of the study structure.

Evaluation criteria
The statistical indices adopted as performance criteria for the ETo equations were the percentage bias (PBIAS), the absolute mean error (MAE), the root mean square error (RMSE), and the Nash and Sutcliffe (1970) efficiency (NSE): where X i and Y i denote ETo values estimated by Penman-Monteith and remaining methods, respectively.X and Y represent the corresponding mean ETo values; i refers to the i th value of ETo and n is the number of data observations.Percentage bias indicates if a model generally under or overestimates PM-ETo.Mean absolute error indicates the average error magnitude, while RMSE assigns more weight to larger errors.Nash-Sutcliffe efficiency ranges from -∞ to 1, where 1 indicates a perfect model and negative values are negligible and indicate that the model is as good as the average of the observations (Schaefli and Gupta, 2007).Uncalibrated models are used here as benchmarks to evaluate if there is any improvement obtained by calibrating the models' coefficients or if the model must, nevertheless, be considered unacceptable (NSE < 0).

Data overview
The daily observations from weather stations gathered for the Cerrado from 2000 to 2020 are displayed in Figure 3.The number of observations for each station ranged from 10.1 to 20.5 years.AWSs yield more data from recent periods, while CWSs yield more data from the early 2000s, although some have been decommissioned more recently.To overcome the discontinuity of observation data sets in space and time, satellite data and reanalysis products may be considered a possible alternative, though they might require bias correction or adjustment to the well-watered conditions of reference sites.
The meteorological variables show a clear distinction between wet and dry seasons.For example, relative humidity decreases throughout the dry season until Oct, when the rainy season begins.The dry season is also associated with less cloud cover due to lower relative humidity.This results in a lower spread of solar radiation distribution during these months.In the dry season (winter), minimum daily temperatures are lower, while wind speed is generally higher.
Daily mean air temperature values varied from 4.9 to 34.6 °C (avg = 23.3°C) in the dry season and from 9.2 to 35.6 °C (avg = 24.9°C) in the wet season.Similarly, wind speed ranged from 0 to approximately 8.0 m s -1 in both seasons, with averages equal to 1.15 and 1.29 for the wet and dry seasons, respectively.Average relative humidity for the dry and wet seasons were 59.9 % and 73.6 %, respectively.The average solar radiation was similar for both wet (19.3 MJ m -2 d -1 ) and dry (18.5 MJ m -2 d -1 ) seasons.However, it registered a higher variance in the wet season due to rainfall and the increased cloudiness.
Reference evapotranspiration generally follows the same behavior of maximum daily air temperature and solar radiation.Solar radiation and wind speed generally increase at the end of the dry season.or the mainpart the high solar radiation results fom low relative humidity and cloudiness.Thus, low relative humidity associated with high solar radiation and wind speed produced the highest average ETo in Sept (4.97 mm d -1 ).In contrast, lower temperatures and solar radiation resulted in the lowest average ETo in June (3.30mm d -1 ).
In the wet season, Rs showed the highest linear correlation to ETo (cor = 0.94), followed by RH (-0.75) and T max (0.68).The variables least correlated to ETo in the wet season were u 2 (0.25) and T min (0.09).In contrast, Rs (0.81), RH (-0.67), and T max (0.66) correlation to ETo declined in the dry season, while for u 2 (0.43) and T min (0.46) it increased.

Improved parameters for the Cerrado
The alternative methods for computing ETo were calibrated using the calibration set considering the annual and season periods.The parameters, i.e., empirical coefficients and correction factors, are displayed in Table 1.The dew point temperature correction factor, a T ,   was negative (-0.42) for the annual period, implying that T dew is generally higher than T min .Negative a T values are generally seen in humid regions (Paredes and Pereira, 2019), similar to here, i.e., negative [positive] value calibrated for the wet [dry] season.T dew is commonly higher than T min when daytime absolute humidity measurements exceed nighttime measurements due to evaporative conditions and humidification from daytime evapotranspiration (Allen, 1996).In contrast, despite the dry season's lower relative humidity and higher wind speeds, a T was lower than 1 °C.In both cases, a T was lower than 1 °C because the temperature (T min ) will approach T dew at reference sites provided the wind speed is relatively calm (Allen, 1996).The average wind speed (1.23 m s -1 ) in the Cerrado is lower than the global average, with lower averages seen in the wet season, as discussed in the previous section.The k Rs factor resulted in values very close to that suggested for the inland region.Although higher values were generally observed near the ocean, they were recorded in the dry season (0.165) instead of the wet season (0.155).
The offset used in the Hargreaves-Samani equation was lower for the Cerrado (14.9), with a slight variance between wet and dry seasons.The empirical coefficient of the PT equation was lower in the wet season (1.13) and higher in the dry season (1.30), when advection is likely higher due to increased wind speed and lower relative humidity.The MK equation calibrated a more significant slope (c 1 ) for the dry season (0.76) as opposed to the wet season (0.67), while the intercept (c 2 ) was higher during the wet season (0.28) and lower in the dry season (-0.26).

Performances of alternative reference evapotranspiration methods
The performance achieved by the alternative ETo models considering the original parametrization and the improved model for the Cerrado is displayed in Figure 4.The original parametrization for HS, PMT, and PT revealed a an overestimated ETo, while MK generally underestimated it.After improving the parameters, all the methods resulted in negligible PBIAS (< 3 %).Before calibration, the MAE was equal to 0.74, 0.75, 0.56, and 0.74 mm d -1 for PMT, HS, PT, and MK, respectively.After calibrating the parameters for the annual period, MAE declined by 17 %, 16 %, 18 %, and 55 % for PMT, HS, PT, and MK, respectively.The MK had the most significant improvement since it showed the greatest bias before calibration.PMT and HS had minor improvements, likely due to limitation of using only temperature data.
The methods showed slight improvement by considering a calibration on a seasonal versus an annual basis.For example, the parameters calibrated on an annual or seasonal basis were very similar for PMT and HS.The temperature correction factor, average wind speed, and k Rs presented minor differences between the seasons (Table 1).Only the PT equation showed significant performance improvement by considering the seasons, with MAE declining up to 31 % in this case.The MK method seemed to have already improved sufficiently for the annual period by considering the slope (c 1 ) and intercept in calibration (c 2 ), while PT relied entirely on the empirical coefficient to scale its estimates in different seasons.
The small or insignificant improvements for the methods on a seasonal basis are, despite the better-adjusted parameters, the result of their lack of information on variables such as wind speed and relative humidity.For instance, other studies also concluded that only minor improvements could be achieved for the HS equation in the dry season for a border region between the Cerrado, Caatinga, and Atlantic rainforest with climate classified as tropical, semi-arid, and humid subtropical for the weather stations studied (Althoff et al., 2019).As wind speed correlation to ETo increases in the dry season, and relative humidity also has a relatively high correlation to ETo, it is difficult for the models to capture the ETo dynamics in this period without them.
Overall, the best performing method after calibration was MK, with NSE = 0.85 for both seasonal and annual parametrizations, followed by PT (seasonal: NSE = 0.79; annual NSE = 0.73), PMT (seasonal: NSE = 0.59; annual: NSE = 0.58), and HS (seasonal: NSE = 0.57; annual: NSE = 0.56).For instance, annual parametrization resulted in an RMSE for MK equal to clustering stations for a region with a similar climate by the Cerrado border did not improve the calibration performance (Althoff et al., 2019).Since the Cerrado climate is predominantly characterized by a tropical savanna climate with strong seasonality, parametrization of coefficients for the entire biome is more practical for stakeholders than station-wise or sub-regional parametrizations.This further reinforces that choosing a single set of coefficients calibrated for the Cerrado is practical and the recommended approach.The scatterplots in Figure 5 reinforce the lower level of sensitivity shown by methods based solely on temperature, where the regression slopes did not show any significant improvement with calibration.This shows that most of the improvements in errors came from reducing bias, but the methods still result in relatively large errors on a daily time scale.Despite both MK and PT being temperature-radiation based, MK was the method that returned the best overall performance, and slope, intercept, and coefficient of determination (R 2 ) closer to the desired values.
However, even when these methods are adopted for water management purposes, e.g., irrigation, the under and/or overestimation in consecutive periods can balance each other out.Thus, to conclude that the methods systematically under or overestimated ETo in specific months, the performance was assessed for each month separately.Figure 6 presents the performance criteria for the annual parametrization, while Figure 7 presents the performance criteria for the seasonal parametrization.Annual parametrization (Figure 0.48 mm d -1 , which is lower by 40 %, 41 %, and 25 % when compared to PMT (0.80 mm d -1 ), HS (0.79 mm d -1 ), and PT (0.63 mm d -1 ), respectively.
Other alternative methods based on machine learning have been reported with good performance in different studies.For instance, Ferreira et al. (2019) used artificial neural networks and a single temperaturebased model for Brazil and achieved performance similar to our results (RMSE = 0.81 mm d -1 ).The authors demonstrated that this was an improvement compared to the HS and PMT equations calibrated for the entire Brazilian territory.However, distributing an equation is still more practical than sharing a machine learning model with stakeholders.Following this idea, using a genetic algorithm, Valle Júnior et al. ( 2020) generated new equations for a Cerrado region under a tropical savanna climate with dry winter similar to ours.However, these equations' performance was still short when compared to the calibrated version of the equations used here with similar input.
Clustering weather stations under similar climates has improved the final performance when modeling ETo for Brazil (Ferreira et al., 2019).However, the clustering algorithm used on a national scale for Brazil grouped stations in regions similar to that delimited by the Brazilian biomes Amazon rainforest, Cerrado, Caatinga, and Atlantic rainforest (Ferreira et al., 2019).This is because, under different climate types, the chances are that a single parametrization would be biased towards the climate type more often seen by weather stations.Additionally, it has also been shown that 6) presents a strong seasonal presence of bias in predictions.For example, temperature-based methods tend to overestimate ETo during the wet season and underestimate it during the dry season.PT shows similar biased behavior.In contrast, MK tends to overestimate ETo from Feb to July, then underestimates it from Aug to Dec.This shows that MK generally overestimates ETo in the period where it is supposedly lower and underestimates it in the period when it is higher.
For parametrization on a seasonal basis, the overall bias across the months was reduced.Interestingly, despite showing sparse performance improvement, the biases for temperature-based methods were better balanced across all months after the seasonal parametrization.PMT and HS presented PBIAS below 5 % for all months except May.This signifies that over more extended periods, e.g., weeks or month-long periods, the estimates made by more straightforward methods can result in better estimates than expected.In contrast, the biases for MK and PT across the months maintained a certain seasonality, resulting in PBIAS values exceeding 5 % in several months.
The alternative ETo methods are also assessed for each station to provide a regional assessment of their performance (Figure 8).Central Cerrado showed a somewhat lower performance.For the temperature-based methods PMT and HS, a few stations resulted in negative NSE, which highlights their inadequacy for predicting ETo.The MK method outperformed the temperature-based methods for all stations and PT in 79 % of the stations.Even when PT outperformed MK, the difference was negligible.Thus, the MK method is recommended regardless of the site within the Cerrado.As for the temperature-based equations, the HS method outperformed the PMT in 72 % of the stations.
Our main recommendation for future research is that a more concentrated effort should be channeled into estimating solar radiation using limited data.Solar radiation is the variable with the highest correlation to ETo, and it is fundamental to any endeavor to improve the estimation process when data is missing.Using the suggested equation (Eq.4), a closer look at solar radiation estimates reveals poor overall performance.For example, even when seasonal calibration was included, Eq. ( 4) resulted in an MAE = 2.88 MJ m -2 d -1 , RMSE = 3.63 MJ m -2 d -1 , and NSE = 0.46.Several alternative methods (Fan et al., 2018) have been proposed and assessed with better outcomes than the traditional formulation proposed by Hargreaves and Samani (1985).Thus, there seems to be room for significant improvement for PMT and HS, especially when their balanced biases across all months are taken into account.

Conclusions
In this study, the most popular alternatives to compute reference evapotranspiration were assessed concerning the Cerrado biome.These methods are the Penman-Monteith Temperature approach (PMT), Hargreaves-Samani (HS), Priestley-Taylor (PT), and Makkink (MK), which are particularly useful when only limited data is available.The parametrization of these methods was also improved by considering an annual and seasonal basis for the region.Improved parametrization for the Cerrado reduced the mean absolute error between 16 % to 55 % for these methods for the annual parametrization, except for the PT method.
Overall, the best performing temperature-radiationbased method was the MK, with an NSE of 0.85 for both annual and seasonal parametrization.The best temperature-based method was the PMT, with an NSE of 0.58 and 0.59 for seasonal and annual parametrization, respectively.Despite the seasonal parametrization failing to improve performance criteria, it was seen to reduce systematic biases on a monthly scale, especially in the case of the temperature-based methods.Future research should focus efforts in the main part on improving solar radiation estimates from limited data.Solar radiation showed the highest correlation with reference evapotranspiration, and the traditional method proposed for its estimation performed poorly in this study.This likely constrained the opportunity for temperature-based methods to estimate reference evapotranspiration adequately.

where
Figure 1 -The Cerrado biome in relation to Brazil and the weather stations used in this study.

Figure 3 -
Figure 3 -Daily observations of climatic data in the Cerrado.

Figure 4 -
Figure 4 -Performance achieved by the alternative reference evapotranspiration methods computed in their original form, and with parameters improved for the Cerrado on an annual and seasonal basis.PBIAS = percentual bias; MAE = mean absolute error; RMSE = root mean square error; NSE = Nash-Sutcliffe Efficiency; PMT = Penman-Monteith temperature approach; HS = Hargreaves-Samani; PT = Priestley-Taylor; MK = Makkink.

Figure 8 -
Figure 8 -Performance achieved by the alternative reference evapotranspiration methods for each station in the Cerrado with parameters improved for the Cerrado on a seasonal basis.Red points are stations where NSE scored below 0. NSE = Nash-Sutcliffe Efficiency; PMT = Penman-Monteith temperature approach; HS = Hargreaves-Samani; PT = Priestley-Taylor; MK = Makkink.

Table 1 -
Summary of calibrated parameters.