Evapotranspiration and components of corn ( Zea mays L.) under micro irrigation systems in a semi-arid environment

Aim of study: This work summarizes the influence of surface drip irrigation (DI) and subsurface drip irrigation (SDI) systems on corn growth indices and actual evapotranspiration (ET c-act ) and its components of plant transpiration (T p ) and soil evaporation (E). Area of study: Karaj, Iran Material and methods: The experimental soil was loamy. The corn ET c-act of each mini-lysimeter was measured based on the water balance method. The E was measured using two mini-lysimeters and T p was estimated from the difference between ET c-act and E. Main results: The resulting data showed that the ET c-act was lower under SDI (384.8 mm) than under DI (423.4 mm). The K cb-m for the corn increased after sowing and peaked during the mid-season stage, with an average value of 0.47, a minimum value of 0.0 and maximum value of 1.52 under DI and 0.53, 0.0 and 1.74 respectively, under SDI. For K e-m , the average, minimum and maximum values were 0.33, 0.20 and 0.58 under DI and 0.23, 0.15 and 0.46 respectively, under SDI. The biomass yield was much higher under SDI (81.90 ton/ha) than under DI (63.21 ton/ha). Less E and more T p occurred under SDI than under DI. SDI achived superior WUE (8.32 kg/m 3 ) compared with DI. Research highlights: SDI was superior to DI based on biomass yield, corn height, stem diameter, and leaf area index which contributed to more favorable soil moisture conditions and low weed incidence; Thus, the SDI system is more productive and would better increase WUE than the DI system. c-act (actual cb cb-adj (adjusted basal crop coefficient by K cb-m (measured basal crop coeffi-cient); K e (soil evaporation coefficient); K e-adj (adjusted soil evaporation coefficient by FAO-56); K e-m (measured evaporation coefficient); LAI (leaf area index); SDI (subsurface drip irrigation); SPAC (soil-plant-atmosphere continuum); T p (plant transpiration); WUE (water use efficiency).


Introduction
Corn (Zea mays L.) is a major cereal crop in Iran, ranking third in cultivated area and production after wheat and rice. Since the susceptibility of corn to drought is one of the production challenges in arid and semi-arid regions, providing enough irrigation water for its growth is required. Gi-ven the importance of this crop and the decreasing availability of agricultural water resources, increasing water use efficiency (WUE) to produce more crops with the available water is highly important to stable agricultural development. Since water is the most limited resource in arid and semi-arid regions, the agricultural sector should produce more food with less water (Zwart & Bastiaanssen, 2004).

Spanish Journal of Agricultural Research
June 2020 • Volume 18 • Issue 2 • e1202 The use of drip irrigation (DI) systems is an effective strategy for increasing water availability in the future (Enciso et al., 2015). The subsurface drip irrigation (SDI) system is the most efficient micro-irrigation methods in arid and semi-arid regions, where the evaporation rate is high during the growing season (Kalfountzos et al., 2007;Sharma et al., 2010). The SDI has greater WUE and saves more water, helping to preserve the nutrients used by crops in comparison to other irrigation methods (Schneider & Howell, 2001;Paul et al., 2013;Panigrahi et al., 2016;Liu et al., 2017;Zhang et al., 2017). The SDI has a larger wetted soil volume than the DI and hence, the volume of available soil for root growth is higher, while the wet radius in the SDI is smaller than that under the DI. In general, under similar irrigation conditions, access to water and nutrients under SDI is increased and root rot and other soil diseases are minimized (Phene & Ruskin, 1995;Kalfountzos et al., 2007).
For the above reason, the SDI has been recommended as a high-efficiency method, reducing water losses through soil surface evaporation and creating more suitable conditions for plant transpiration (Tiwari et al., 2014, Parthasarathi et al., 2017Reddy et al., 2018). A large number of experiments have been conducted to define the main advantages of SDI for several crops, and its superior performance was confirmed in all of them (Seyfi & Rashidi, 2007;Van Donk et al., 2013;Albasha et al., 2015;Biswas et al., 2015;Lamm, 2016). Thus, the SDI can be an alternative to other irrigation methods and could be used to increase the growth of fruits, vegetables and row crops due to the precise application of water and provision of adequate moisture in the root zone (Imtiyaz et al., 2000). Unfortunately, knowledge of irrigation water management at the farm level is poorly developed in Iran and most of the farmers act based on their own experience. Therefore, detailed on-farm information on water and crop could support the design and management of sustainable and beneficial irrigation systems.
The aim of this work was to evaluate the corn growth indices and actual evapotranspiration (ET c-act ), i.e., plant transpiration (T p ) and soil evaporation (E) under DI and SDI systems to understand how the advanced irrigation systems could benefit growers from the point of view of yield and water use.

Experimental site
The experiments were carried out at the Agricultural Engineering Research Institute (35°46ˊ N, 50°55ˊE, 1260 m a.s.l.) during the 2014 and 2015 growing seasons. The area is semi-arid with an average annual precipitation of approximately 279.3 mm. Daily meteorological data (air temperature, relative humidity, wind speed, rainfall, and solar radiation data) were collected from a synoptic meteorology station, 5 km from the field site. Average daily values of meteorological characteristics are shown in Fig. 1. The minimum air temperature during the 2014 and 2015 corn growing seasons ranged from 5.7 to 23.6°C and 9.5 to 25.5°C, respectively. Also, the average maximum air temperature ranged from 20.8 to 41.4°C and 23.2 to 39.7°C, respectively. The daily relative humidity was 11 to 90% and 15 to 87%, respectively. On average, wind speed and solar radiation were higher in 2014 by about 6.4 and 2.7%, respectively (Fig 1).
The experimental soil was loamy with mean volumetric field capacity and permanent wilting point of 22.3 and 9.63%, respectively; the mean soil bulk density was 1.42 g/cm 3 (Table 1).
To monitor water consumption by corn under DI and SDI, eight mini-lysimeters (Dugas & Bland, 1989;Kong et al., 2012) were placed within a corn farm of 18-ha considering having adequate fetch and filled with soil excavated from the study site to resemble the original soil profile conditions. The mini-lysimeters had a diameter of 40 cm and a depth of 70 cm. Every 6 mini-lysimeters were used as three replicates for DI and SDI, respectively. The DI and SDI were equipped with 40 cm emitter apart and discharge of 4 L/h. For the subsurface drip irrigation, the drip-line was buried 30 cm below the soil surface. Inside each mini-lysimeters, three forage corn (Single Cross 704) seeds were planted with 13 cm spacing on 6 th August 2014 and 2015. Water and nutritients were optimally provided for the mini-lysimeters. Each of the treatments received the same amount of water and nutrients through DI or SDI during each growing season. The depth of water applied for each of the two experimental seasons is shown in Fig. 2. Additionally, two unplanted mini-lysimeters were used to measure evaporation from the soil surface under DI and SDI, and were placed near the other mini-lysimeters.

Irrigation management
The required irrigation water depth was estimated daily using the FAO-Penman-Monteith model (Allen et al., 1998;Eqs. (1) and (2)) confirmed for the Karaj region by Dehghanisanij et al. (2004) and corn crop coefficient (K c ) recommended for Karaj by Farshi et al. (1997) where ET c is the crop evapotranspiration; ET o is the reference evapotranspiration (mm/day); R n is the net radiation 3 Evapotranspiration and components of corn (Zea mays L.) under micro irrigation systems  (MJ/m 2 / day); G is the soil heat flux density (MJ/m 2 /day); T is the mean temperature (ºC); U 2 is the wind speed at 2-m height (m/s); ɣ is the psychrometric constant (kPa/ ºC); ∆ is the slope vapor pressure curve (kPa/ºC); e a is the actual vapor pressure (kPa); e s is the saturation vapor pressure (kPa). Fertilizers were applied through the irrigation water from the 3 and 4 leafed stage of corn growth to 45 days before harvesting. The corn received 250 kg/ha ammonium phosphate and 200 kg/ha urea.

Actual evapotranspiration (ET c-act )
Daily corn water consumption or actual evapotranspiration (ET c-act ) of each mini-lysimeter was measured based on the water balance method using Eq. (3) (Allen et al., 1998): where P is the rainfall (mm); I is the irrigation depth (mm); D p is the water loss through drainage from the mini-lysimeter (mm); R is the runoff (mm) which here was zero and Δs is the change in soil water storage in the mini-lysimeter (mm). The change in soil water storage (Δs) was determined using Eq. (4) (Allen et al., 1998): where S t and S t-1 are the available water in the root zone at the beginning and end of the selected period (mm), respectively.

Soil surface evaporation (E) and plant transpiration (T p )
Evaporation from the soil surface (E) was measured using two mini-lysimeters. When these mini-lysimeters were placed inside the soil, their edges were about one cm above the soil surface, and the soil inside the mini-lysimeter was about 1 to 1.5 cm below the edge. The E was estimated from the difference between the amount of intake and drainage water in the mini-lysimeter at each irrigation interval. Plant transpiration (T p ) was estimated from the difference between ET c-act and E, by the following equation (Moran et al., 2009):

Basal crop coefficient (K cb ) and soil evaporation coefficient (K e )
Crop coefficient can be applied as a single crop coefficient which is influenced by evaporation and transpiration together and dual crop coefficient that is expres-sed by soil evaporation coefficient (K e ) and basal crop coefficient (K cb ), separately (Allen et al., 1998). The measured K e and K cb (K cb-m and K e-m ) are defined using Eqs. (6) and (7) K cb-m and K e-m were also compared with estimated values based on FAO-56 (Allen et al., 1998). In this study, the K cb and K e suggested by FAO-56 were adjusted (K cb-adj and K e-adj ) based on the climatic conditions of the study area.

Corn growth indices and water use efficiency (WUE)
Leaf area index (LAI), plant height (H), stem diameter (D) and yield (biomass) were measured during the growing season. LAI was measured with the electronic leaf area-meter, CI-202, seven times during the growing season. WUE (kg/m 3 ) was calculated using Eq. (8) (Sakthivadivel et al., 1999) where yield represents the biomass of corn:

Statistical analysis
Statistical analysis was done to evaluate the influence of different point source irrigation systems on actual evapotranspiration (ET c-act ) and growth indices of corn by using the SAS package. LSD (least significant difference) tests were used to compare and rank the treatment means. Differences were declared significant at p ≤ 0.05.

Variations in actual evapotranspiration (ET c-act ), plant transpiration (T p ) and soil evaporation (E)
The values for ET c-act and T p of corn and E under DI and SDI during the corn growing season are shown in Fig.  3. The ET c-act varied from 2.1 to 9.7 mm/day, depending on the crop growth stage and climatic condition. The daily ET c-act and T p increased rapidly and peaked 48 days after planting at the mid-season stage with average values of 9.7 and 7.6 mm/day under DI and 9.5 and 8.7 mm/day under SDI. The minimum ET c-act values occurred at the initial stage with an average of 2.5 and 2.1 mm/day under DI and SDI. During the first week after planting, the most share of ET c-act was E, because the soil was kept wet for corn seed germination. Transpiration increased faster under the SDI than the DI (Fig. 3) because of moisture distribution close to the root mass and better root development under the SDI at the initial stage.
Significantly lower corn ET c-act and E values were obtained under the SDI than under DI but for both seasons, the T p was higher under SDI than DI (  (2007), Liu et al. (2017) and Zhang et al. (2017). Evaporation accounted for a small proportion of ET c-act and decreased with time and increasing LAI. Transpiration was smaller than E at the initial growth stage (Fig.  3), and the total contribution to ET c-act at that stage was from E because more water was lost through soil surface evaporation and the crop canopy was not fully developed yet as reported in Allen et al. (1998) andValentin et al. (2020). The T p increased to a peak at the mid-season stage (Liu et al., 2002;Kang et al., 2003;Majnooni-Heris et al., 2012). During this stage, E was reduced and the reduction was compensated with higher amounts of water for T p , coinciding with increasing LAI until the crop achieved near or full ground cover ( Table 2).
The average annual T p during the two growing seasons was 221.9 and 250.2 mm, and that for E was 198.0 and 157.6 mm under DI and SDI, respectively. Under SDI, the soil surface usually remained drier than under DI. Accordingly, the total T p was 53.0 and 64.2% of corn ET c-act under DI and SDI. Also, the average annual E was 46.9 and 35.9 % of corn ET c-act under DI and SDI during two growing seasons (Table 2). At the initial stage during which the portion of the bare or low covered soil surface was high and considering that under DI, moisture accumulation occured on the soil surface, the DI showed 27.70% more evaporation (82.9 mm) than SDI (64.9 mm).

Relationship between the ratio of plant transpiration (T p ) and soil evaporation (E) to actual corn evapotranspiration (ET c-act ) and leaf area index (LAI)
Based on the soil-plant-atmosphere continuum (SPAC), the variations in T p , and ET c-act are affected by meteorological data, soil moisture and plant factors (Zhou et al., 2017). The variations in T p /ET c-act and LAI for each growing season are presented in Fig. 4. The data indicated that the ratio was controlled by soil surface cover and LAI. The average T p /ET c-act varied from 0.0 at sowing to 82 and 87% at full growth under DI and SDI. The value of T p /ET c-act increased faster when LAI was smaller than 3.0. Thus, two polynomial functions were computed to show the relation between the ratio of T p /ET c-act and E/ ET c-act with LAI under SDI and DI (Eqs. 8 and 9; Fig. 4), as reported by other authors (e.g., Kang et al., 2003;Majnooni-Heris et al., 2012).  (10) It was presumed, as occurred for the T p /ET c-act ratio, that the E/ETc-act was also affected by LAI and surface soil moisture (Liu et al., 2002). Fig. 5 shows the trends in the ratio of E/ET c-act with LAI under DI and SDI during the growing seasons. The E/ET c-act decreased significantly witth increases in LAI. The relationship between LAI and E/ET c-act based on the experimental data was calculated as follows: From these relationships, the E/ET c-act decreased sharply when LAI was at about 3, which was attributed to lower evaporation due to the development of crop canopy (Fig. 5). Under SDI, the soil surface between crops usually remains dry, so that evaporation loss from the soil surface is low but moisture distribution at the root zone is suitable for crop growth. Thus, the SDI can improve water uptake by reducing soil evaporation (Liu et al., 2002).

Variations in basal crop coefficient (K cb ), soil evaporation (K e ) and LAI
The mean variations for the two years for K cb-m and K e-m compared to K cb-adj , and K e-adj are presented in Fig. 6. The values of K cb-m increased from 0.0 to its peak value in the mid-season stage. The average, minimum For the whole growing season, the values of K cb-m were smaller than the adjusted values based on FAO-56.
Evapotranspiration and components of corn (Zea mays L.) under micro irrigation systems Table 2

. The ratio of soil evaporation (E) and plant transpiration (T p ) to actual corn evapotranspiration (ET c-act ) under surface drip irrigation (DI) and subsurface drip irrigation (SDI) systems during 2014 and 2015 corn growing seasons.
For each year, means within a column bearing the same letter do not differ significantly at the 0.05 level of probility This was attributed to the overestimation of ET o by FAO-Penman-Monteith model for this region (Dehghanisanij et al., 2004). For the different irrigation systems, the peak values of K cb-.m were obtained under SDI at the mid-season stage, and the K e-m under the DI peaked at the initial stage. The K e-m decreased with crop development during the growing season due to the increase in the percentage of the shaded area by the plant canopy. Again, Fig. 6 emphasizes that the evaporation from soil surface was higher than the transpiration from the crop at the initial stage and with an increase in plant shading, evaporation was lower than transpiration during crop development and mid-season stages. Moreover, the K e-m varied temporally during the corn growing season. Thus, the average K e-m peaked at the initial stage, and decreased gradually during the growing season, reaching a minimum value at the mid-season stage ( Table 3). The amounts of K cb-FAO during the initial stage, crop development and mid-season stage of corn growth (K cbini , K cb-dev , and K cb-mid ) were 0.15, 0.15-1.15 and 1.15, respectively (Allen et al., 1998). The values of K cb-dev and K cb-mid changed on the basis of plant height, wind speed and relative humidity in different regions. Therefore, the recommended K cb values were adjusted to 0.15, 0.70, and 1.20 under DI and 0.15, 0.76, and 1.29 under SDI during the initial stage, crop development, and mid-season stages, respectively (Table 3). The maximum value of K cb-adj was 1.30 under SDI at the mid-season stage, which was attributed to the higher T p . The K e-adj varied temporally during the corn growing season. Also, the K e-adj value was higher at the initial stage and gradually decreased, reaching a minimum value of only 0.14 under SDI at the mid-season stage (Fig. 6).
The K cb-adj and K e-adj were higher than K cb-m and K e-m during the initial stage and crop development. Whe-ther for DI or SDI, the K cb-adj values were overestimated when compared to K cb-m for the whole growing season. Comparisons of K cb-m and K e-m vs K cb -adj , and K e -adj values are shown in Fig. 7. For both DI and SDI, the relationships between K cb-m and K e-m as well as K cb -adj , and K e -adj were linear linear. The K cb -adj under DI performed better than that under SDI. The slopes of the linear regression were 0.93 and 0.85 with a coefficients of determination of 0.75 and 0.71 for K cb -adj under DI and SDI, respectively.
The K e-adj under DI and SDI was overestimated as compared to K e-m during the initial and crop development stages and underestimated during the mid-season stage. Also, K e-adj under DI provided a better performance than that under SDI. The slopes of the linear regression were 0.67 and 0.67 with a coefficients of determination of 0.25 and 0.19 for K e -adj under DI and SDI, respectively.
The difference between K cb-m and K e-m and adjusted values by FAO-56 clearly emphasizes the difficulty of applying K cb and K e values across locations due to varying climatic and agricultural management factors like irrigation method and frequency (Katerji & Rana, 2014). As a response to the crop development, the K e was higher than the K cb at the initial stage, and with an increase in LAI and plant shading, the resulting decreases in soil surface evaporation could be used to estimate the variations in K cb and K e as functions of LAI (Fig. 8).
Overall, the relationship between LAI and K cb and K e based on the experimental data was calculated using second-order polynomial equations with high coefficients of determination (R 2 ) as follows: The K cb of corn at the mid-season stage was larger because of the large LAI at this stage and consequently, the smaller soil evaporation, compared to plant transpiration. Similar relationships have been reported for bean and canola (De Medeiros et al., 2001;Majnooni-Heris et al., 2012).

Variations in leaf area index (LAI), plant height (H), stem diameter (D), wet and dry mass yield and water use efficiency (WUE)
The variations in mean LAI, H, D, biomass yield and WUE under DI and SDI are presented in Table 4 and Fig. 9. The LAI increased slowly in the initial stages and more rapidly reached its peak at about 3.4 m 2 /m 2 under DI and 4.6 m 2 /m 2 under SDI. Also, the H increased slowly at the initial stages and more rapidly at the mid-season stage, with the tallest plants (229.2 cm) under SDI at the mid-season stage. Moreover, the D also increased during the growing season and peaked during mid-season with values of 3.2 cm under DI and 3.6 cm under SDI. Between the irrigation systems, the lowest values of LAI, H and D occurred under DI (Table 4 and Fig. 9) as reported in Chuanyan & Zhongren (2007).

Spanish Journal of Agricultural Research
June 2020 • Volume 18 • Issue 2 • e1202 Table 3. Mean values of measured (m) and adjusted (adj) soil evaporation coefficient (K e ) and basal crop coefficient (K cb ) under surface drip irrigation (DI) and subsurface drip irrigation (SDI) systems during 2014 and 2015 corn growing seasons.
For each year, means within a column bearing the same letter do not differ significantly at the 0.05 level of probility  (Table 4). This better performance under SDI could be explained by the conservation of optimal moisture status in the root zone, which favoured water and nutrient uptake by the crop (Zotarelli et al., 2008;Badr et al., 2010).
In summary, the results of this study showed that the lowest ET c-act for corn occurred under SDI, mainly because of reduced evaporation (E) from the soil surface in comparison to the DI system. Consequently, higher T p and lower E rates were observed under the SDI system. The T p /ET c-act started from 0 at sowing and peaked at the mid-season stage when the LAI was also at peak levels.
The DI and SDI systems had different influences on the LAI of corn during the growing season. A better understanding of these components (T p , E, E/ET c-act and T p /ET c-act ) can provide important insights to water saving under irrigated corn production. The highest LAI occurred under the SDI system. Overall, the SDI system reduced soil evaporation loss and increased the efficiency of water consumption; it also produced superior biomass yield. Thus, the SDI system is more productive and would better increase WUE than the DI system.