Next Article in Journal
The Influence of Hydrometeorological Conditions on Changes in Littoral and Riparian Vegetation of a Meromictic Lake in the Last Half-Century
Next Article in Special Issue
Parameterization and Application of Stanghellini Model for Estimating Greenhouse Cucumber Transpiration
Previous Article in Journal
Documenting a Century of Coastline Change along Central California and Associated Challenges: From the Qualitative to the Quantitative
Previous Article in Special Issue
Measuring and Modelling Soil Evaporation in an Irrigated Olive Orchard to Improve Water Management
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Maize Evapotranspiration Estimation Using Penman-Monteith Equation and Modeling the Bulk Canopy Resistance

Colegio de Postgraduados, Carretera México-Texcoco Km. 36.5, Montecillo, Texcoco 56230, Mexico
*
Author to whom correspondence should be addressed.
Water 2019, 11(12), 2650; https://doi.org/10.3390/w11122650
Submission received: 5 November 2019 / Revised: 4 December 2019 / Accepted: 11 December 2019 / Published: 15 December 2019
(This article belongs to the Special Issue Evapotranspiration and Plant Irrigation Strategies)

Abstract

:
Some techniques, such as the Katerji and Perrier approach, estimate the bulk canopy resistance (rc) as a function of meteorological variables and then calculate the hourly evapotranspiration using the Penman–Monteith equation, so that traditional crop coefficients are not needed. As far as we know, there are no published studies regarding using this method for a maize crop. The objective of this study was to calibrate and validate the canopy resistance for an irrigated continuous maize crop in the Midwestern United States (US). In addition, we determined the effect of derivation year, bowen ratio, and the extent of canopy. In this study we derive empirical coefficients necessary to estimate rc for maize, five years (2001–2005) were considered. A split-sample approach was taken, in which each year’s data was taken as a potential calibration data set and validation was accomplished while using the other four years of data. We grouped the data by green leaf area index (GLAI) and the Bowen ratio (β) by parsing the data into a 3 × 3 grouping: LAI (≥2, ≥3, and ≥4) and |β| (≤0.1, ≤0.2, and ≤0.3). The best fit data indicated reasonably good results for all nine groupings, so that the calibration coefficients derived for the conditions LAI ≥ 2 and |β| ≤ 0.3 were taken in light of the longer span associated with LAI ≥ 2 and the larger number of hours. For the calibrations in this subgroup, the results indicate that the annual empirical coefficients for rc are nearly the same and equally effective, regardless of the year used for calibration. Our validation included all the daytime hours regardless of β. Thus, it was concluded that the calibration at our site was independent of the derivation year. Knowledge of the Bowen ratio was useful in calibration, but accurate ET estimates (validation) can be obtained without knowledge of the Bowen ratio. Validation resulted in hourly ET estimates for irrigated maize that explained 83% to 86% of the variation in measured ET with an accuracy of ± 0.2 mm.

1. Introduction

Irrigation is the largest use of water resources in many agricultural regions of the world. The challenges of planning for water use and water resources management in a changing climate include the increased stress on the world’s fresh water reserves due to increased use of water for food production, contamination of rivers, lakes, and aquifers, and a lack of support for major new water projects. There is potential need to alleviate the stress on the world’s fresh water system with human population increasing and the need for more food production; conserving water through irrigation scheduling provides that potential [1].
The method, frequency, and duration of irrigation strongly affect crop yields and farm productivity [2]. An efficient use of water resources through irrigation scheduling requires an accurate estimation of crop water use to realize a potential savings.
Evapotranspiration (ET) estimates are required for determining timely irrigation needs when scheduling irrigation with the checkbook approach [3].
ET is one of the main components of the water cycle [4]. ET plays an important role in several disciplines, such as hydrology, agriculture, climatology, and meteorology [5].
The traditional approach to estimate crop ET is to first estimate reference ET from a standard surface and to apply an empirical crop coefficient, such as those that were presented by Doorenbos and Pruitt [6] and Wright [7,8]. The Penman–Monteith (PM) equation was accepted as the standard for reference ET by the FAO in 1990 along with a set of procedures for calculating the various parameters for short and tall crops [9,10]. The PM computes the latent heat flux (λET) as:
λ E T = Δ ( R n G ) + ρ a c p ( e s e a ) r a Δ + γ ( 1 + r c r a )
where λ E T is the latent heat of vaporization ( M J m 2 t ) , Δ is the slope of the saturation vapor-pressure curve ( k P a ° C ) , Rn is the net radiation at the crop surface ( M J m 2 t ) , G is the soil heat flux at the soil surface ( M J m 2 t ) , ρa is the mean air density ( k g m 3 ) , cp is the specific heat of air at constant pressure ( M J k g ° C ) , es is the saturated vapor pressure at the air temperature ( k P a ) , ea is the vapor pressure of air ( k P a ) , esea is the vapor pressure deficit ( k P a ) , γ is the psychrometric constant ( k P a ° C ) , ra is aerodynamic resistance ( t m ) , and rc is bulk canopy resistance ( t m ) , hereafter referred to as canopy resistance.
The PM equation, strictly speaking, applies to crops with full cover and it has been shown to be unreliable for ET estimates under partial cover [11,12]. Shuttleworth and Wallace provided a physically-based approach for modeling ET for the entire range of canopy cover [13]. However, the Shuttleworth and Wallace model requires a number of model parameters that are not commonly available [14].
When a sensitivity analysis of ET estimated for crops of different heights and in contrasting water status while using the PM equation was conducted, it was concluded that it is very important to correctly represent aerodynamic resistance for short crops under well-watered conditions. Canopy resistance plays a secondary role for short crops, accounting for only 10% to 25% of the variation in ET [15]. This could explain the good performance of simple models of canopy resistance to estimate reference grass ET while using the PM equation [9]. However, it was concluded that vapor pressure deficit and canopy resistance are critical in estimating ET for tall crops and for water stressed crops [15].
Research on ET has provided sound theoretical knowledge and practical applications that have been validated through field measurements. However, the transfer of theoretical advances into field practice remains far short of its potential [16,17]. Verstraeten et al. [18] provided an extensive survey on different methods to estimate ecosystem ET and soil moisture content and summarize the ET assessment techniques (from sap flow to remote sensing). Current operational approaches to estimate ET for an agricultural crop require crop coefficients [6,7,8]. Worldwide research provides crop coefficients for most existing crops as a function of phenological stage. However, variations in the accuracy of crop ET estimates due to the use of fixed crop coefficient values have been identified [19,20,21]. Differences as large as 40% between the crop coefficient derived ET values that were reported by Allen et al. [9] and ET values that were experimentally obtained have been found [22].
Several corrections have been proposed for standard crop coefficients to account for varying crop roughness values as functions of crop height, wind speed, and relative humidity [9,23,24]. Other studies proposed using a functional relationship to generate crop coefficient values that consider crop development stage and canopy density as functions of leaf area index and soil water availability [20,25]. Several researchers have developed theoretical expressions for crop coefficients based on the PM formulation to explain crop coefficient variations, and demonstrated that crop coefficients are a function of canopy resistance, aerodynamic resistance, and meteorological conditions [17,26,27]. Shuttleworth provides theoretical analyses that support the use of the PM formulation to make a direct estimation of crop water requirements, and emphasizes the need for effective surface resistances that are equivalent to those which exists for crop coefficients, which is a major hurdle in using the PM-based approach [27].
The direct estimation of ET from the PM formulation for agricultural crops when the canopy and aerodynamic resistances (rc and ra) are known avoids the uncertainty associated with the use of crop coefficients [27]. The single layer or “big leaf” model assumes that latent heat lost by a crop is controlled by the entire canopy resistance to water vapor loss to the atmosphere [28]. For vegetative surfaces, rc is the combined resistance of all leaves and soil surfaces and the resistance to vapor loss inside the canopy [27]. Meteorological and biophysical variables influenced the resistance. It has been argued that rc has an indirect aerodynamic component, but it is often assumed that rc mainly represents stomatal response for full canopies [29].
Bulk canopy and aerodynamic resistances have been parameterized for a reference crop, most often defined as an actively growing grass canopy of uniform height under non-limiting water conditions, but not for other crops [9]. Alves et al. proposed an approach, where the aerodynamic resistance is computed from the top of the canopy (h) to the reference height (z), where wind speed (u) and humidity are measured [30]. This expression estimates the aerodynamic resistance from the plane of momentum sink to the reference height and from the plane of momentum sink to the canopy height. Afterwards, the difference of these two values is calculated. While assuming neutral conditions, the wind speed at the top of the canopy and the wind profile under neutral stability conditions, friction velocity, is obtained [30]. An equation for the aerodynamic resistance, ra, is resolved while using estimates of zero plane displacement height (d) [31] and roughness length (zo) based on crop height (h), [32]:
r a = ln ( z d h d ) ln ( z d z o ) k 2 u = ln ( z 0.75 h 0.25 h ) ln ( z 0.75 h 0.13 h ) k 2 u
where k is von Karman’s constant, unitless. For stable and unstable cases, stability functions are needed to calculate the aerodynamic resistance [33,34]. The height h can be estimated while using the Richards plant growth model [35]:
h = a h ( 1 + e b h c h c D ) 1 d h
where D is day of year and ah, bh, ch, and dh are annually derived from non-linear regression analysis.
Using Equation (2) to estimate aerodynamic resistance implies that the canopy resistance has two components: one corresponding mainly to stomatal resistance and the other to leaf boundary layer and turbulent transfer inside the canopy [17]. The use of this approach has the advantage in that it does not require any assumptions about the roughness lengths for heat and vapor transfer that may have no precise physical meaning [30].
Crop canopy resistance values are typically determined in research studies by inverting the PM equation [30], or by scaling the leaf resistance to canopy level by weighting the average stomatal resistance of the leaves by the leaf area index (LAI) [9,32]. Methods that determine the resistance through alternative means have been derived. For example, Jarvis considered that stomatal conductance (the reciprocal of resistance) is a function of environmental variables (as obtained from controlled environment studies) and their effects multiplicative [36]. Estimates using this approach have yielded good canopy resistance approximations explaining 70% of the variation; however, its use is year and site specific [37]. Alves and Pereira argue that rc is dependent on environmental variables, so that using purely physiological approaches might not be suitable [37]. Hatfield and Allen found consistent ET estimations of cotton, grain sorghum, and grass forage, while using an empirical canopy resistance [38]. Katerji and Perrier (1983) proposed an estimate of canopy resistance as a function of an equilibrium resistance (r*) and aerodynamic resistance (ra), which has been used to estimate ET of alfalfa, rice, grass, lettuce, sweet sorghum, grain sorghum, soybean, tomato, and wheat [22,37,39]. Katerji and Rana found that treating canopy resistance as a constant might be acceptable in short crops, but not in taller crops, like sunflower, soybean, and sorghum [22].
A study concluded that the Katerji and Perrier approach considers a physiological component that responds to the main environmental variables in a similar way to the multiplicative method of Jarvis [36], and also considers an aerodynamic component [8]. According to Katerji and Rana [22], the work of Katerji and Perrier [39] pointed out the failure of many researchers in recognizing vapor transfer inside the canopy is a component of canopy resistance and it is dependent on wind speed. Katerji and Perrier [39] recommend a simpler alternative to the Jarvis [36] and Shuttleworth [27] methods. This approach simulates the canopy resistance by the following empirical relationship:
r c r a = a r * r a + b
where a and b are parameters to be experimentally determined by calibration when bulk canopy resistance (rc), aerodynamic resistance (ra), and equilibrium resistance (r*) are known. We will refer to the evapotranspiration that was estimated from the use of rc (Equation (4)) in Equation (1) as the Katerji Perrier evapotranspiration approach (ETKP). The equilibrium resistance (r*) is the resistance when only only the radiative energy streams determine evapotranspiration, i.e., λ E T e = Δ ( R n G ) Δ + γ . Substituting this value for λET in Equation (1), rc can be found through a relation between rc/ra and r*/ra, where the radiative controlled r* is calculated as [17]:
r * = Δ + γ Δ ρ a c p ( e s e a ) γ ( R n G )
The relation between rc/ra and r*/ra depends on the phenological stage of the vegetation and soil water status [29]. The generation of the coefficients a and b for use in Equation (4) requires an estimate of crop specific bulk canopy resistance. Inverting the PM formulation and solving for rc we obtain:
r c = r a Δ γ ( R n G ) λ E T r a Δ + γ γ + ρ c p γ ( e s e a ) λ E T
It is suggested that an approach that estimates rc as a function of meteorological conditions, as in Equation (4), could offer a way to overcome the limitations of applying the PM method to estimate ET [29]. The approach has been applied to estimate ET from grass, alfalfa, wheat, tomato, and rice [37]; it has received some criticism due to the need for calibrating for each location. However, it was found that the soybean coefficients a and b that were calibrated in Italy worked well for another soybean crop that was grown in France [40]. For short crops, the use of Equation (4) gave little improvement as compared to the constant rc approach of FAO [41]. In 2011, it was shown that this semi-empirical approach of Equation (4) for rc is capable of producing unbiased estimates of ET, while other approaches are biased too low for short crops and too high for tall crops [42]. In 2013, eight models were tested of the Jarvis type for rc in soybean, where multiple variables, both micrometeorological (net radiation: Rn; air temperature: Ta; relative humidity: RH; wind speed at 3 m; U3; vapor pressure deficit: VPD; solar zenith angle; and aerodynamic resistance: ra) and plant variables (Leaf Area Index: LAI; crop height: h; and, leaf stomatal resistance rs) were used in different combinations. The models were able to explain 66 to 93% of the variability in the measured ET that was obtained from a Bowen Ratio energy balance [43].
In 2017, the Katerji and Perrier approach was validated for rice and potato in china, the study concluded that this approach is a better alternative for estimating crop ET than the traditional FAO Penman–Monteith approach (PM) crop coefficient used [44].
Another consideration is that the formulation for ra is only for neutral conditions. Previous studies have calibrated the method (a and b in Equation (4)) for near normal conditions. For example, in a study from Lecina et al., the calibration data was restricted to conditions of −0.5 ≤ β ≤ 0.5, while Alves et al. used only data with −0.3 ≤ β ≤ 0.3, and determined near neutral conditions while using the Richardson number [30,41].
Based on our literature review, we did not find any results, where the approach of Katerji and Perrier for estimating hourly evapotranspiration was calibrated for a maize crop. Thus, this study was undertaken. The objective of this study was to analyze the feasibility of using this approach to estimate the bulk canopy resistance for maize, a tall crop, so that the PM formulation can be directly used to estimate the hourly ET from weather data without the use of a crop coefficient. We will refer to the use of Equations (1)–(6) as the Penman Monteith Katerji and Perrier (ETKP) approach. By deriving two fixed coefficients, a and b, to estimate the canopy resistance for a maize crop, the ETKP from Equation (1) can provide surface ET from a canopy through the growing season on an operational basis while using weather station data and the canopy height (h). The extent to which the ETKP coefficients, a and b, change from year to year for maize was determined and the accuracy with which the measured ET can be estimated while using the ETKP under varying leaf area index (LAI) and under different ranks of Bowen ratio (β).

2. Materials and Methods

2.1. Study Area

The study was conducted while using data from 2001 through 2005 from a production field in the University of Nebraska Carbon Sequestration Program (CSP) [45]. The production field (41°09′54.2″ N, 96°28′35.9″ W, and 361 m), which is located in Eastern Nebraska at the University of Nebraska Agricultural Research and Development Center near Ithaca, NE, is 47 ha in size, has a center pivot irrigation system, and is maintained as a no-till continuous maize production system. The study site provides sufficient upwind fetch over a uniform cover to adequately measure mass and energy fluxes while using an eddy covariance system. We refer to the measurement of ET from this system as ETEC. The soil is a deep silty clay loam that consists of four soil series: Yutan, Tomek, Filbert, and Filmore. Prior to initiating the study, the site had a 10-year history of maize-soybean rotation under no-till management. In 2001, the site was uniformly tilled (to incorporate fertilizers) by disking prior to initiating the study and it remained under no-till through 2005. Seed was planted below the crop residue from previous years and the standard best management practices were followed. The amount of N fertilizer applied was adjusted in the spring before planting to augment nitrate already in the soil, according to recommended guidelines [46]. Maize Pioneer 33P67 (Bt hybrid) was planted in 2001 and 2002, Pioneer 33B51 in 2003, and Pioneer 33B53 in 2004, all at approximately 84,000 plants ha−1; DeKalb 63–75 was planted in 2005 at approximately 82,000 plants ha−1.

2.2. Field Measurements

Water was applied in the field to maintain a minimum of 50% available soil moisture in the root depth zone. The soil water conditions in the root zone were continuously monitored at four depths (0.10, 0.25, 0.5, and 1.0 m) while using Theta probes (Delta-T Devices, Cambridge, UK) at three locations within the field.
From the CSP data, we used daylight (7:30–18:30 h) growing season hourly averages of maize ET flux from eddy covariance measurements, net radiation, soil heat flux, sensible heat flux, air temperature, relative humidity, wind speed, and incoming solar radiant flux (hours during which precipitation/irrigation occurred were excluded from the study). The hourly Bowen ratios were calculated as the ratio of the sensible heat flux to latent heat flux. The measurements began around planting time in 2001 (May 14) and they have continued to the present day [47]. Net radiation (Rn) and soil heat flux (G) were available for maize at the site but they were not used in the calibration and validation, since these measurements are typically not provided from local weather stations. Thus, the methodology would be consistent with the available data to future users of the calibration coefficients. These variables were estimated in our analysis, as follows:
R n = ( 1 α ) R s + R n l
G = 0.1 R n
where α is the albedo of the crop (0.23), Rs is the incoming short wave radiation, and Rnl is the net long wave radiation that we estimate according to Allen [9], and G is the soil heat flux. The albedo for the maize crop was assumed to be 0.23; this value is also assumed for a well-irrigated short grass or a tall alfalfa crop when the reference ET is estimated [9,23].
The leaf area index (LAI) was calculated from periodic harvesting of one meter row lengths of crop measured with the LI-3100C Area Meter (Li-Cor, Lincoln, NE, USA). Samples were harvested at six plots every 10 to 14 days until the crops reached reproductive stage and again before harvest. The total and green LAI (GLAI) were calculated based upon the sorting of leaf material into green and brown components. Maize crop height was concurrently measured with measurements of leaf area. The height measurements were used in a non-linear regression to derive annual values of ah, bh, ch, and dh for use in prescribing daily height estimates (Equation (3)).

2.3. Estimating Evapotranspiration and Canopy Resistance

For each hour (7:30–18:30 Central Standard Time) during the growing season the values of ra, r*, and λ E T were used to derive the apparent canopy resistance, rc. These data were sub-divided with attention to the GLAI and Bowen ratio ranks (β). For calibration, the values of rc/ra and r*/ra were treated as the dependent and independent values, respectively, in a linear regression. This resulted in coefficients for a and b for each year and for the defined conditions. The estimated values that exceeded ± 3σest were considered to be outliers and eliminated from the analysis, where σest is the standard error of the estimate.
The use of the ETKP approach supposes that a full canopy is present and is contributing to the aerial and canopy resistances. Thus, data when GLAI < 2 were not included in our study. The data were grouped according to GLAI ≥ 2, 3, and 4 and |β| ≤ 0.1, 0.2, and 0.3 for each year from which the calibration coefficients a and b were derived. This provided a 3 × 3 set of coefficients for each of the years 2001 to 2005. In this manner, we determined the calibration effects of different ranks of Bowen ratio (β) and GLAI and their consistency from year to year.
The calibration coefficients a and b were applied in the validation process to the entire growing season data (GLAI ≥ 2) for all years, with the exception of the calibration year, regardless of the Bowen ratio value.
Hourly rc values were estimated (Equation (4)) for daylight conditions (7:30–18:30 h); nocturnal hours and hours during which precipitation occurred were excluded from the study. Some studies recommend nocturnal times be excluded, because the vapor pressure deficit and the available energy are relatively small during these times and, thus, can generate negative or undefined values for the equilibrium and canopy resistances (r* and rc, respectively) [29].
Estimated canopy resistance values (rc from Equation (4)) and aerodynamic resistance, ra, (estimated while using the Perrier model in Equation (2)) were applied in the PM equation (Equation (1)) to calculate λET (Wm−2), following the recommended methods [9,48] while giving a direct approach (ETKP) to estimating ET without crop coefficients. The volume of data was around 3500 values of measured hourly ET for the five years studied.

2.4. Data Analyses

The data analysis included a calibration for each year of the Katerji and Perrier approach for each stratification i.e., each combination of GLAI thresholds (GLAI ≥ 2, 3, and 4) with each classification of Bowen ratios (|β| ≤ 0.1, 0.2, and 0.3)]. The means and standard deviations of a, b, and σ were also calculated for the data within each stratification.
A split-sample approach was used in calibration and validation [49]. Coefficients a and b were derived for each year’s data and used in a validation by estimating the hourly ET (ETKP) for the remaining years (calibration year excluded). The performance was documented in terms of variance explained (r2), the standard error of estimate (σest), and bias that was present in the estimates (slope of the ET comparison curve). The calibration coefficients that produce the best overall validation (ETKP compared to ETEC) are assumed to be the ‘best’ values to use with a maize crop.

3. Results

In Table 1, the derived coefficients a and b for estimating rc (Equation (4)) and the standard error of estimate for the calibrations (σest) are shown based on the combinations of GLAI and Bowen ratio (β) for each year (hours were included on the basis of GLAI thresholds (GLAI ≥ 2, 3, and 4) and Bowen ratio classes (|β| ≤ 0.1, 0.2, and 0.3)).
Figure 1 represents a typical calibration result where coefficients a and b in Equation (4) are estimated by linear regression. In this case, the solution is for the conditions GLAI ≥ 2 and|β| ≤ 0.2, which correspond to the top middle part of Table 1. The linear relationship between rc/ra and r*/ra was statistically tested for all of the cases shown in Table 1 and a statistically significant relationship was found at a 95.0% confidence level for all of the studied cases.
The standard deviation (σ) of the slope (a) across the years is a measure of the consistency of the calibration from year to year. In examining Table 1, we can see that σ for the slope for any GLAI threshold increases from |β| ≤ 0.3 to |β| ≤ 0.1. For example, in the case of GLAI ≥ 2, the σ is 0.04, 0.03, and 0.07 for|β| ≤ 0.3, 0.2, and 0.1, respectively. Contributing to the inconsistency from year to year is the fewer number of hours in the more restrictive conditions of |β| ≤ 0.1 (regardless of LAI), which result in a smaller range in r*/ra. For example, the mode of r*/ra for |β| ≤ 0.1, 0.2, and 0.3 at LAI ≥ 2 in 2001 is 7.81, 9.17, and 22.94, respectively.
While the magnitude of σest decreases as β gets closer to zero (i.e., values of 2.52, 2.37, and 2.14 for GLAI ≥ 2 and |β| ≤ 0.3, 0.2, and 0.1. respectively), the standard deviation (σ) of the σest values indicate higher consistency from year to year at |β| ≤ 0.3 with a value of 0.30.
Figure 2 further demonstrates the consistency among the calibration results, which shows the linear relationship between rc/ra and r*/ra. It was decided to choose the coefficients from the data group |β| ≤ 0.3 and GLAI ≥ 2 to be used in the validation of the ETKP, because this case provides a larger and more consistent dataset to validate the approach.
Table 2 shows the results of the split-sample validation for GLAI ≥ 2 and |β| ≤ 0.3 for all of the years. The ET estimations using the ETKP approach explained approximately 85% of the variance in the measured ETEC, with a standard error of estimate about 0.2 mm for all the cases studied. This was attained without employing the β value for each hour.
The coefficients for 2003 (a = 0.74 and b = 2.35) produced an ET validation slope of 1.0 and the highest r2, but showed little difference in validation when compared to the other coefficient pairs. The slope of the relationship was nearly one in all cases, which indicated no bias in the estimator, demonstrating the accuracy with which ETEC can be estimated while using the ETKP under varying LAI conditions. As an example, the 2003 calibration coefficients were applied to all daylight growing season hours, regardless of the β value. It was excluded the 2003 data and those hours in which precipitation occurred. (Figure 3).

4. Discussion

Some techniques estimate the bulk canopy resistance (rc) as a function of meteorological variables and then calculate the hourly evapotranspiration while using the Penman–Monteith equation so that traditional crop coefficients are not needed. The objective of this study was to calibrate and validate the canopy resistance for an irrigated continuous maize crop. In addition, the effect of derivation year, bowen ratio, and the extent of canopy was determined. The empirical coefficients that were necessary to estimate rc for maize in each of the five years (2001–2005) were derived. A split-sample approach was taken in which each year’s data was taken as a potential calibration data set and validation was accomplished while using the other four years of data.
The data was grouped by green leaf area index (GLAI) and the Bowen ratio (β) by parsing the data into a 3 × 3 grouping: LAI (≥2, ≥3, and ≥4) and |β| (≤0.1, ≤0.2, and ≤0.3). The best fit data indicated reasonably good results for all nine groupings, so the calibration coefficients that were derived for the conditions LAI ≥ 2 and |β| ≤ 0.3 were taken in light of the longer span that was associated with LAI ≥ 2 and the larger number of hours. The results indicate the annual empirical coefficients for rc are nearly the same and equally effective, regardless of the year used for calibration, for the calibrations in this subgroup. Our validation included all daytime hours, regardless of β. Thus, it was concluded that the calibration at our site was independent of the derivation year. This is consistent with the findings by Rana and Katerji [40], proposing that calibration shows an interannual stability. Knowledge of the Bowen ratio was useful in calibration, but accurate evapotranspiration estimation can be obtained without knowing it. Validation resulted in hourly ET estimation for irrigated maize that explained 83% to 86% of the variation in measured ET with an accuracy of ± 0.2 mm.
A limitation of the ETKP approach is that the method worked reasonably well for well watered conditions and with LAI bigger than 2. It is necessary to investigate an approach that could be used for LAI smaller than 2; in such conditions, the soil is not full covered and the soil evaporation is important for the ET estimation. The canopy height is needed, but it is not a difficult value to observe in an operational sense. Future work should include additional comparisons of the calibration coefficients across sites with different climates and across hybrids of maize before a common use of this method. Additionally, it is necessary to validate whether the approach can be used or not under rainfed conditions.
The main advantage of the proposed method is that uncertainties that were introduced when estimating the crop ET using a crop coefficient will be avoided; however, it is necessary to first calibrate the method before using it. An advantage of the present study with respect to other studies is the volume of data used, five years of information were used in the validation of the ETKP approach, under different LAI conditions and Bowen ratio values. Our study showed that the calibration was independent of the derivation year.

5. Conclusions

It was concluded that hourly estimates of ET for maize in daylight hours can be made without systematic bias over multiple years. The validation results indicate that the calibration year does not affect the performance of the technique. Note that the Bowen ratio was not used in the validation, so, although the calibration should be conducted with the knowledge of the Bowen ratio, the application does not require knowing it.
The results that are presented in this study show that it is possible to use the ETKP formulation to obtain an unbiased estimate of maize ET when the GLAI is 2 or above and for the Bowen ratio ranks analyzed (|β| ≤ 0.3).
The results indicate that the annual empirical coefficients to estimate the canopy resistance are nearly the same and equally effective, regardless of the year used for calibration. The validation of the ETKP method resulted in hourly ET estimates for irrigated maize that explained 83% to 86% of the variation in measured ET with an accuracy of ±0.2 mm. This study concludes that reliable maize crop evapotranspiration can be made without the use of crop coefficients.

Author Contributions

H.F.-M. conceived the idea of the research; he was in charge of the data review and collection, as well as the data analysis and the writing of the article; N.M.-M. contributed to the analysis, writing and formatting of the article.

Funding

This research received no external funding.

Acknowledgments

We thank the financial support of the University of Nebraska and the National Council of Science and Technology of Mexico (CONACyT) to make possible this study. We also thank to the Professors from the Natural Resources School of the University of Nebraska-Lincoln.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sivakumar, B. Global climate change and its impacts on water resources planning and management; assessment and challenges. Stoch. Environ. Res. Risk Assess. 2011, 25, 583–600. [Google Scholar] [CrossRef]
  2. Walker, W.R. Chapter 1. The Practice of Irrigation in Guidelines for Designing and Evaluating Surface Irrigation Systems. Paper No. 45. 1989. Available online: http://www.fao.org/docrep/T0231E/T0231E00.htm (accessed on 15 September 2019).
  3. Melvin, S.R.; Yonts, C.D. Irrigation Scheduling: Checkbook Method; Extension Publication EC709; University of Nebraska-Lincoln: Lincoln, NE, USA, 2004. [Google Scholar]
  4. Aguilar, A.; Flores, H.; Crespo, G.; Marin, M.; Campos, I.; Calera, A. Performance Assessment of MOD16 in Evapotranspiration Evaluation in Northwestern Mexico. Water 2018, 10, 901. [Google Scholar] [CrossRef] [Green Version]
  5. Ramoelo, A.; Majozi, N.; Mathieu, R.; Jovanovic, N.; Nickless, A.; Dzikiti, S. Validation of Global Evapotranspiration Product (MOD16) using Flux Tower Data in the African Savanna, South Africa. Remote Sens. 2014, 6, 7406–7423. [Google Scholar] [CrossRef] [Green Version]
  6. Doorenbos, J.; Pruitt, W.O. Guidelines for Predicting Crop Water Requirements; Irrigation and Drainage Paper No. 24 (Review); FAO-ONU: Roma, Italy, 1977; p. 144. [Google Scholar]
  7. Wright, J.L. Crop coefficients for estimates of daily crop evapotranspiration. In Irrigation Scheduling for Water and Energy Conservation in the 80s, Proceedings of the Irrigation Scheduling Conference, Chicago, IL, USA, 14–15 December 1981; ASAE: St. Joseph, MI, USA, 1981; pp. 18–26. [Google Scholar]
  8. Wright, J.L. New evapotranspiration crop coefficients. J. Irrig. Drain. Div. ASCE 1982, 108, 57–74. [Google Scholar]
  9. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration (Guidelines for Computer Crop Water Requirements); FAO Irrig. and Drain. Paper No. 56; Food and Agricultural Organization of the United Nations: Rome, Italy, 1998; Available online: http://www.fao.org/docrep/X0490E/X0490E00.htm (accessed on 15 September 2019).
  10. Allen, R.G.; Pruitt, W.O.; Wright, J.L.; Howell, T.A.; Ventura, F.; Snyder, R.; Itenfisu, D.; Steduto, P.; Berengena, J.; Yrisarry, J.B.; et al. A recommendation on standardized surface resistance for hourly calculation of reference ETo by the FAO56 Penman-Monteith method. Agric. Water Manag. 2005, 81, 1–22. [Google Scholar] [CrossRef]
  11. Lafleur, P.M.; Rouse, W.R. Application of an energy combination model for evaporation from sparse canopies. Agric. For. Meteorol. 1990, 49, 135–153. [Google Scholar] [CrossRef]
  12. Irmak, S.; Mutiibwa, D. On the dynamics of evaporative losses from Penman-Monteith with fixed and variable canopy resistance during partial and complete Maize cover. Trans. ASABE 2009, 52, 1139–1153. [Google Scholar] [CrossRef]
  13. Shuttleworth, W.J.; Wallace, J.S. Evaporation from sparse crops: An energy combination theory. Quart. J. R. Meteorol. Soc. 1985, 111, 839–855. [Google Scholar] [CrossRef]
  14. Farahani, H.J.; Ahuja, L.R. Evapotranspiration modeling of partial canopy/residue-covered fields. Trans. ASAE 1996, 39, 2051–2064. [Google Scholar] [CrossRef]
  15. Rana, G.; Katerji, N. A measurement based sensitivity analysis of the Penman-Monteith actual evapotranspiration model for crops of different height and in contrasting water status. Theor. Appl. Climatol. 1998, 60, 141–149. [Google Scholar] [CrossRef]
  16. Farahani, H.J.; Howell, T.A.; Shuttleworth, W.J.; Bausch, W.C. Evapotranspiration: Progress in measurement and modeling in agriculture. Trans. ASABE 2007, 50, 1627–1638. [Google Scholar] [CrossRef]
  17. Pereira, L.S.; Perrier, A.; Allen, R.G.; Alves, I. Evapotranspiration: Review of concepts and future trends. J. Irrig. Drain. Eng. 1999, 125, 45–51. [Google Scholar] [CrossRef]
  18. Verstraeten, W.W.; Veroustraete, F.; Feyen, J. Assessment of evapotranspiration and soil moisture content across different scales of observation. Sensors 2008, 8, 70–117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Jagtap, S.S.; Jones, J.W. Stability of crop coefficients under different climate and irrigation management practices. Irrig. Sci. 1989, 10, 231–244. [Google Scholar] [CrossRef]
  20. Ritchie, J.T. Model for predicting evaporation from a row crop with incomplete cover. Water Resour. 1972, 8, 1204–1213. [Google Scholar] [CrossRef] [Green Version]
  21. Ritchie, J.T.; Johnson, B.S. Soil and plant factors affecting evaporation. In Irrigation of Agricultural Crops; Stewart, B.A., Nielsen, D.R., Eds.; Michigan State University: East Lansing, MI, USA, 1990; pp. 363–390. [Google Scholar]
  22. Katerji, N.; Rana, G. Modelling evapotranspiration of six irrigated crops under mediterranean climate conditions. Agric. For. Meteorol. 2006, 138, 142–155. [Google Scholar] [CrossRef]
  23. Allen, R.G.; Pruitt, W.O.; Businger, J.A.; Fritschen, L.J.; Jensen, M.E.; Quinn, F.H. Evaporation and transpiration. In ASCE Handbook of Hydrology; ASCE: Reston, VA, USA, 1996; pp. 125–252. [Google Scholar]
  24. Allen, R.G.; Smith, M.; Pruitt, W.O.; Pereira, L.S. Modifications to the FAO crop coefficient approach. In Proceedings of the International Conference Evapotranspiration Irrigation Scheduling, San Antonio, TX, USA, 3–6 November 1996; Camp, C.R., Sadler, E.J., Yoder, R.E., Eds.; ASAE/IA/ICID: St. Joseph, MI, USA, 1996; pp. 124–132. [Google Scholar]
  25. Ritchie, J.T. Evaluating Irrigation Needs for Southeastern U.S.A. 1975, pp. 262–273. Available online: http://agris.fao.org/agris-search/search.do?recordID=US201303058926 (accessed on 15 December 2019).
  26. Hubbard, K.G. Climatic factors that limit daily evapotranspiration in sorghum. Clim. Res. 1992, 2, 73–80. [Google Scholar] [CrossRef]
  27. Shuttleworth, W.J. Towards one-step estimation of crop water requirements. Trans. ASABE 2006, 49, 925–935. [Google Scholar] [CrossRef]
  28. Komatsu, H. A general method of parameterizing the big-leaf model to predict the dry-canopy evaporation rate of individual coniferous forest stands. Hydrol. Process. 2004, 18, 3019–3036. [Google Scholar] [CrossRef]
  29. Perez, P.J.; Lecina, S.; Castellvi, F.; Martínez-Cob, A.; Villalobos, F.J. A simple parameterization of bulk canopy resistance from climatic variables for estimating hourly evapotranspiration. Hydrol. Process. 2006, 20, 515–532. [Google Scholar] [CrossRef]
  30. Alves, I.; Perrier, A.; Pereira, L.S. Aerodynamic and surface resistances of complete cover crops: How good is the “big leaf”? Trans. ASAE 1998, 41, 345–351. [Google Scholar] [CrossRef]
  31. Thom, A.S. Momentum, mass and heat exchange of plant communities. In Vegetation and the Atmosphere; Monteith, J.L., Ed.; Vol. 1 Principles; Academic Press: New York, NY, USA, 1975; pp. 57–109. [Google Scholar]
  32. Monteith, J.L. Principles of Environmental Physics; Edward Arnold: London, UK, 1973. [Google Scholar]
  33. Tolk, J.A.; Howell, T.A.; Steiner, J.L.; Krieg, D.R. Aerodynamic characteristics of corn as determined by energy balance techniques. Agron. J. 1995, 87, 464–473. [Google Scholar] [CrossRef]
  34. Verma, S.B. Aerodynamic resistances to transfers of heat, mass, and momentum. In Estimation of Areal Evapotranspiration; IAHS Pub. No. 177; IAHS: Vancouver, BC, Canada, 1989. [Google Scholar]
  35. Gregorczyk, A. Richards plant growth model. J. Agron. Crop Sci. 1998, 181, 243–247. [Google Scholar] [CrossRef]
  36. Jarvis, P.G. The interpretation of the variations in leaf water potential and stomatal conductance found in canopies in the field. Phil. Trans. R. Soc. Lon. B 1976, 273, 595–610. [Google Scholar] [CrossRef]
  37. Alves, I.; Pereira, L.S. Modeling surface resistance from climatic variables? Agric. Water Manag. 2000, 42, 371–385. [Google Scholar] [CrossRef]
  38. Hatfield, J.L.; Allen, R.G. Evapotranspiration estimates under deficient water supplies. J. Irrig. Drain. Eng. 1996, 122, 301–308. [Google Scholar] [CrossRef]
  39. Katerji, N.; Perrier, A. Modelisation de l’evapotranspiration reelle d’une parcelle de luzerne: Role d’une coefficient culturale. Agronomie 1983, 3, 513–521. [Google Scholar] [CrossRef] [Green Version]
  40. Rana, G.; Katerji, N. Direct and indirect methods to simulate the actual evapotranspiration of irrigated overhead table grape vineyard under Mediterranean conditions. Hydrol. Process. 2008, 22, 181–188. [Google Scholar] [CrossRef]
  41. Lecina, S.; Martínez-Cob, A.; Pérez, P.J.; Villalobos, F.J.; Baselga, J.J. Fixed versus variable bulk canopy resistance for reference evapotranspiration estimation using the Penman-Monteith equation under semiarid conditions. Agric. Water Manag. 2003, 60, 181–198. [Google Scholar] [CrossRef] [Green Version]
  42. Katerji, N.; Rana, G.; Fahed, S. Parameterizing canopy resistance using mechanistic and semi-empirical estimates of hourly evapotranspiration: Critical evaluation for irrigated crops in the Mediterranean. Hydrol. Process. 2011, 25, 117–129. [Google Scholar] [CrossRef]
  43. Irmak, S.; Mutibwa, D.; Payero, J.; Marek, T.; Porter, D. Modeling soybean canopy resistance from micrometeorological and plant variables for estimating evapotranspiration using one-step Penman-Monteith approach. J. Hydrol. 2013, 507, 1–18. [Google Scholar] [CrossRef]
  44. Zhao, P.; Luers, J. Parameterization of Evapotranspiration Estimation for Two Typical East Asian Crops. Atmosphere 2017, 8, 111. [Google Scholar] [CrossRef] [Green Version]
  45. Verma, S.B.; Dobermann, A.; Cassman, K.G.; Walters, D.T.; Knops, J.M.; Arkebauer, T.J.; Suyker, A.E.; Burba, G.G.; Amos, B.; Yang, H.; et al. Annual carbon dioxide exchange in irrigated and rainfed maize-based agroecosystems. Agric. For. Meteorol. 2005, 131, 77–96. [Google Scholar] [CrossRef] [Green Version]
  46. Shapiro, C.A.; Ferguson, R.B.; Hergert, G.W.; Dobermann, A.; Wortmann, C.S. Fertilizer Suggestions for Corn; NebGuide G74-174-A; Cooperative Extension, Institute of Agriculture and Natural Resources, University of Nebraska-Lincoln: Lincoln, NE, USA, 2001. [Google Scholar]
  47. Suyker, A.; Verma, S.B. Evapotranspiration of irrigated and rainfed maize-soybean cropping systems. Agric. For. Meteorol. 2009, 149, 443–452. [Google Scholar] [CrossRef] [Green Version]
  48. Allen, R.G.; Walter, I.A.; Elliott, R.L.; Howell, T.A.; Itenfisu DJensen, M.E.; Snyder, R.A. The ASCE Standardized Reference Evapotranspiration Equation; American Association of Civil Engineers: Reston, VA, USA, 2005; 203p. [Google Scholar]
  49. Abbot, M.B.; Refsgaard, J.C. (Eds.) Distributed Hydrological Modelling; Kluwer Academic Publishers: Norwell, MA, USA, 1996; 336p. [Google Scholar]
Figure 1. The calibration result for irrigated maize in 2001 with hours for calibration restricted to green leaf area index GLAI ≥ 2 and only hours with |β, Bowen ratio| ≤ 0.2. The variables are canopy resistance (rc), aerodynamic resistance (ra) and equilibrium resistance (r*). The variation of rc/ra explained by r*/ra is 86% and the empirical coefficients are a = 0.73 and b = 2.36. p value is 0.0000.
Figure 1. The calibration result for irrigated maize in 2001 with hours for calibration restricted to green leaf area index GLAI ≥ 2 and only hours with |β, Bowen ratio| ≤ 0.2. The variables are canopy resistance (rc), aerodynamic resistance (ra) and equilibrium resistance (r*). The variation of rc/ra explained by r*/ra is 86% and the empirical coefficients are a = 0.73 and b = 2.36. p value is 0.0000.
Water 11 02650 g001
Figure 2. Regression lines resulting from a comparison of rc/ra to r*/ra from which coefficients a and b were derived for 2001–2005 while using different green leaf area index (GLAI) ranges and bowen ratio (β) combinations. The variables are canopy resistance (rc), aerodynamic resistance (ra) and equilibrium resistance (r*).
Figure 2. Regression lines resulting from a comparison of rc/ra to r*/ra from which coefficients a and b were derived for 2001–2005 while using different green leaf area index (GLAI) ranges and bowen ratio (β) combinations. The variables are canopy resistance (rc), aerodynamic resistance (ra) and equilibrium resistance (r*).
Water 11 02650 g002
Figure 3. Validation of ETKP against ETEC for hours when GLAI ≥ 2 using the 2003 coefficients with slope, r2, and σest of 1.0, 86%, and 0.20 mm, respectively, for all hours, regardless of β value. ETEC is the ET measured by the eddy correlation approach and ETKP is the estimated ET from the Penman Monteith equation with the Katerji and Perrier approach [19] for canopy resistance.
Figure 3. Validation of ETKP against ETEC for hours when GLAI ≥ 2 using the 2003 coefficients with slope, r2, and σest of 1.0, 86%, and 0.20 mm, respectively, for all hours, regardless of β value. ETEC is the ET measured by the eddy correlation approach and ETKP is the estimated ET from the Penman Monteith equation with the Katerji and Perrier approach [19] for canopy resistance.
Water 11 02650 g003
Table 1. All possible calibration combinations where three Bowen ratio ranks and three GLAI (green leaf area index) thresholds were applied to the period 2001–2005. The empirical coefficients a and b are shown along with the standard error of estimate for the best fit equation rc/ra = ar*/ra + b. Other statistics over the five-year period are given at the bottom of each column including the Mean, and the standard deviation (σ).
Table 1. All possible calibration combinations where three Bowen ratio ranks and three GLAI (green leaf area index) thresholds were applied to the period 2001–2005. The empirical coefficients a and b are shown along with the standard error of estimate for the best fit equation rc/ra = ar*/ra + b. Other statistics over the five-year period are given at the bottom of each column including the Mean, and the standard deviation (σ).
YearGLAI ≥|β| ≤ 0.3|β| ≤ 0.2|β| ≤ 0.1
abσestabσestabσest
200120.673.152.430.732.362.240.742.151.95
200220.672.982.260.742.232.060.791.61.78
200320.742.352.240.81.582.070.880.611.85
200420.742.352.750.752.122.710.761.612.41
200520.752.562.910.81.962.770.880.932.69
Mean0.712.682.520.762.052.370.811.382.14
σ0.040.370.300.030.300.350.070.610.40
200130.682.922.190.722.372.030.81.551.81
200230.672.922.230.742.192.040.781.631.76
200330.732.372.250.791.622.080.880.611.85
200430.712.662.690.72.562.620.692.392.24
200530.732.622.660.782.12.650.851.272.61
Mean0.702.702.400.752.172.280.801.492.05
σ0.030.230.250.040.350.320.070.640.36
200140.672.952.240.722.332.070.791.551.84
200240.662.822.10.732.11.930.771.631.73
200340.732.312.240.791.662.050.890.531.85
200440.742.132.530.741.972.420.771.362.01
200540.752.412.760.821.662.710.920.422.55
Mean0.712.522.370.761.942.240.831.102.00
σ0.040.350.270.040.290.320.070.580.33
Table 2. The validation of estimated hourly evapotranspiration from the Katerji Perrier approach (ETKP) against measured evapotranspiration (ETEC) from the eddy correlation system. The symbols indicate canopy resistance (rc), aerodynamic resistance (ra), and equilibrium resistance (r*). The empirical coefficients necessary to estimate rc for maize in each of the five years were derived (rc/ra = ar*/ra + b). A and B and σest (validate) are from the linear regression between ETEC and ETKP. Calibration years and coefficients are repeated here from Table 1 and shown on the left. Validation data includes all hours from 730 to 1830 CST with LAI ≥ 2, regardless of thermal stability. No hours from the calibration year were included in the validation.
Table 2. The validation of estimated hourly evapotranspiration from the Katerji Perrier approach (ETKP) against measured evapotranspiration (ETEC) from the eddy correlation system. The symbols indicate canopy resistance (rc), aerodynamic resistance (ra), and equilibrium resistance (r*). The empirical coefficients necessary to estimate rc for maize in each of the five years were derived (rc/ra = ar*/ra + b). A and B and σest (validate) are from the linear regression between ETEC and ETKP. Calibration years and coefficients are repeated here from Table 1 and shown on the left. Validation data includes all hours from 730 to 1830 CST with LAI ≥ 2, regardless of thermal stability. No hours from the calibration year were included in the validation.
YearCalibration rc/ra = ar*/ra + bValidation ETEC Measured Vs ETKP Estimated
abσestr2σest (validate)AB
(mm)
20010.673.152.430.830.210.970.04
20020.672.982.260.850.200.990.01
20030.742.352.240.860.201.000.01
20040.742.352.750.850.201.000.01
20050.752.562.910.850.201.010.02

Share and Cite

MDPI and ACS Style

Meraz-Maldonado, N.; Flores-Magdaleno, H. Maize Evapotranspiration Estimation Using Penman-Monteith Equation and Modeling the Bulk Canopy Resistance. Water 2019, 11, 2650. https://doi.org/10.3390/w11122650

AMA Style

Meraz-Maldonado N, Flores-Magdaleno H. Maize Evapotranspiration Estimation Using Penman-Monteith Equation and Modeling the Bulk Canopy Resistance. Water. 2019; 11(12):2650. https://doi.org/10.3390/w11122650

Chicago/Turabian Style

Meraz-Maldonado, Nora, and Héctor Flores-Magdaleno. 2019. "Maize Evapotranspiration Estimation Using Penman-Monteith Equation and Modeling the Bulk Canopy Resistance" Water 11, no. 12: 2650. https://doi.org/10.3390/w11122650

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