# Evaluation of Two Irrigation Scheduling Methods and Nitrogen Rates on Corn Production in Alabama

Authors: Jose Franco Da Cunha Leme Filho; Brenda V. Ortiz; Kipling S. Balkcom; Damianos Damianidis; Thorsten J. Knappenberger; Mark Dougherty
Journal: International Journal of Agronomy
 (2020)

---

## Abstract

Regulations on nutrient application amounts and environmental impacts of fertilizers are promoting advances in agricultural management strategies to optimize irrigation application and N fertilization in corn. Previous studies have found a relationship between irrigation application, available water in the soil, and N fertilizer uptake. The objective of this study was to evaluate interactions between two irrigation scheduling methods and four N rate applications (0-control, 202, 269, and 336 kg ha−1) on grain yield, aboveground biomass, plant N concentration, N uptake, and nitrogen use efficiency in corn. The study was conducted at the Tennessee Valley Research and Extension Center (TVREC) during two growing seasons (2014 and 2015). The irrigation scheduling methods consisted of (i) the pan evaporation method, which is based on managing the crop’s estimated evapotranspiration (ET) using pan evaporation values and the crop’s consumptive water use and (ii) the sensor-based irrigation scheduling method based on soil matric potential values recorded by soil moisture tension sensors installed in the field. Irrigation amounts from both irrigation scheduling methods indicated that less water was applied with the sensor-based method. The different amounts of irrigation applied associated with the two irrigation scheduling methods did not impact grain yield, aboveground biomass, and NUE. In general, NUEs values decreased with increased N rates, which means that additional N fertilizer added to the soil was not converted into grain yield or/and adsorbed by plants; therefore, more N remained in the soil, increasing the risk for environmental problems.


---

## Body

## 1. Introduction

The challenge of making farming profitable and minimizing the environmental effects of agriculture has increased the adoption of new technologies and practices to increase crop production efficiency. Regarding corn production, two major factors affecting yield are soil water availability and nitrogen (N) fertilization; therefore, irrigation and N application are important management factors for producers. Irrigation is a relevant tool in Alabama and this practice has been increasing among farmers, even though the state enjoys average rainfall amounts of about 1422 mm per year [1]. However, the amount and distribution of rainfall are highly unpredictable because they can vary drastically depending on the year. Rainfed corn production is risky with large variations in yield from one year to the next. In Alabama, even though the number of irrigated harvested acres of cropland increased from 1997 to 2012, only approximately 44 thousand hectares out of the 870 thousand hectares of harvested cropland were irrigated [2]. Lack of water from rainfall and/or irrigation can impact the morpho-physiology of corn in terms of cellular and whole-plant effects [3, 4]. Thus, it can result in substantial yield reduction and reasserting the importance of irrigation for consistent corn yields. Meanwhile, N fertilizer is also crucial for corn production; Qiu et al. [5] stated that N is the most limiting nutrient for agricultural production; and N fertilizer is an important component to maximize the yield of most nonlegume crops. Woli et al. [6] evaluated data of corn hybrids from 1960 to 2000 and determined that grain yield increased by 65% and total plant biomass by 45% with agronomic optimum N input.Several approaches, such as improved irrigation technologies that include variable rate irrigation and more efficient irrigation scheduling methods, can be adopted by farmers for more effective use of limited water supplies to maximize water use efficiency (WUE) [7, 8]. The aim of an efficient irrigation scheduling program is to optimize the WUE, replenishing the water deficit within the root zone while minimizing N leaching below this depth [9]. Currently, multiple irrigation scheduling methods have been developed to help producers maintaining adequate soil water content levels in the root zone. Technological advances have improved the devices that continuously monitor changes in soil moisture status. These devices are able to monitor soil moisture during and after irrigation; also they control the amount of water applied. Phene and Howell [10] confirmed that irrigation systems can be controlled accurately with sensors that monitor the soil matric potential within the root zone. There are, however, simpler methods to conduct irrigation scheduling. The pan evaporation method using the data from weather stations can be the simplest, cheapest, and most practical meteorological method to measure local atmospheric evaporation demand [11]. Additionally, values obtained using the pan evaporation method are important references for water resource assessment and monitoring evaporative climate change [12]. It is recognized that the adoption of appropriate irrigation scheduling practices could increase or maintain yields and maximize profit for producers. Also, the water savings can reduce the risk of overirrigation and environmental impacts and consequently improving agricultural sustainability [13].Nitrogen fertilization indicates that the highest yield increasing effect for corn on different soils occurs when the three macroelements (NPK) are considered [14, 15]. Many physiological processes associated with corn development are enhanced with N fertilization [16]. Variation in N application affects crop growth and yield components like potential kernel set [17]. Numerous studies have shown that N fertilization is highly correlated with corn growth; thus, lower N supply can promote the reduction of plant height, leaf area index (LAI), crop photosynthetic rate, radiation use efficiency, and plant N uptake [18–20].The integration of effective water and fertilizer management strategies is essential for increasing crop production and maximizing the economic net return while sustainability is maintained [16, 21–23]. Management strategies for water and N can minimize crop production costs and environmental impacts while maintaining crop performance and economic returns [24]. Efficient N utilization to obtain high yields requires adequate water supply for the crop. For instance, a reduction in plants transpiration as a result of water stress can cause decrease in N uptake [25], including reduced development of steam and leaf cells [26]. On the other hand, a surplus of irrigation and N fertilizer can increase residual soil NO3-N and it may leach or denitrify [27], reducing the efficiency of applied N. Studies have demonstrated that crop growth and yield response to N fertilization varies under different water management conditions. However, there are few studies addressing interactions between N management strategies and the soil water status in Alabama. Therefore, studies to understand N uptake dynamics and determine the best N management practices under two different irrigation scheduling methods are important to help Alabama producers optimize corn yields and avoid adverse environmental impacts. The objectives of this study were (1) to determine the effects of two irrigation scheduling methods, pan evaporation, and sensor-based scheduling, on N supply on corn N uptake, (2) to evaluate the interaction of irrigation scheduling methods and N on Nitrogen Use Efficiency (NUE), and (3) to assess the impact of different N rates on corn yield and NUE.


## 2. Materials and Methods

### 2.1. Site Description

The field experiment was conducted at Tennessee Valley Research and Extension Center located in Belle Mina, Alabama (34°39′24″N 86°52′45″W, 183 m above mean sea level) during the 2014 and 2015 growing seasons. The dominant soil series at the research site is a Decatur silty clay loam (fine, kaolinitic, thermic Rhodic Paleudults) with 0% to 10% slope [28]; rainfall conditions during both growing seasons (April–August) were below the 30-year average (1971–2000) according to the experimental station data. Total rainfall during the 2014 growing season was 11 mm below average and 19 mm below average in 2015. However, in 2014, rainfall distribution per month with respect to the historic average (1971–2000) was higher in April (+39 mm) and June (+42 mm), but rainfall in May (−50 mm), July (−13 mm), and August (−29 mm) was below average. In 2015, above normal rainfall with respect to the 30-year average was observed in April (+21 mm) and May (+6 mm) contrasting with rainfall in June (−21 mm), July (−1 mm), and August (−24 mm) that was below average.


### 2.2. Experimental Design

A split-plot design with three replications was implemented for this study in which irrigation regime was the main plot and N fertilizer rates were the subplots. Two irrigation scheduling methods (pan evaporation method and sensor-based irrigation method), and four N rates (0-control, 202, 269, and 336 kg ha−1) were tested for plant N uptake, NUE, and yield differences. Plots consisted of eight rows 11.8 m long and 5.3 m wide with a row spacing of 76.2 cm and the seeding rate was 89,000 plants ha−1. Urea ammonium nitrate (UAN 32-0-0) was used as the fertilizer source, which was split applied with 1/3 of the total applied at planting and the remainder in a side-dress application at the V6 growth stage. Each eight-row plot was divided in half, 4 rows received 202, 269, or 336 kg N ha−1 depending on the randomization, and the other adjacent 4 rows did not receive any N which allowed for NUE estimation. Plots received 224.1 kg ha−1 in a preplant application of P and K. All the cultural practices were performed according to Alabama Cooperative Extension System (ACES) recommendation to maximize corn yield. Each plot was independently irrigated with four overhead sprinkler nozzles located in each corner of the plot.According to ACES, the recommended N rate for irrigated corn in the state is 224.1 kg ha−1 in order to achieve 11299 kg ha−1 yield. If the yield goal is greater than 12555 kg ha−1, the N recommendation is to apply 1 kg of N per 45.3 kg of anticipated corn yield [29]. The Pioneer 1690YHR corn hybrid was planted on April 4th in 2014 and April 24th in 2015 in single rows and managed with conservation tillage and rye (Secale cereale) cover crop. In both growing seasons, the cover crop was planted on 11/21/13 and 11/07/14, respectively. It received 33.6 kg ha−1 of N fertilizer (UAN 32-0-0) in both seasons to maximize biomass production.


### 2.3. Water Management Treatments

Two irrigation scheduling methods tested in our study differed from each other in the way to estimate the amount of water stored in the soil during the growing season. The pan evaporation method was based on calculating the crop’s estimated evapotranspiration (ET) using the pan evaporation values and the crop’s consumptive water use. On the other hand, the sensor-based irrigation relied on measurements of soil matric potential provided by sensors in order to apply the correct amount of irrigation at the right time. The water use efficiency (WUE) of the two irrigation scheduling methods was included in our evaluation to establish differences in the total amount of water applied in the plots.


#### 2.3.1. Pan Evaporation Irrigation Scheduling

Daily pan evaporation data was reported from the weather station located at the research site (http://www.awis.com). The daily irrigation requirement, ET, was calculated using the following equation: ET = 0.90 × PAN × CC where evapotranspiration, ET (mm day−1), was calculated from 90% pan evaporation, PAN (mm day−1) adjusted for fractional canopy cover (CC) [30, 31]. Canopy cover was determined weekly by measuring the open canopy distance (cm) between rows with a tape measure. Fractional canopy cover was calculated using the following formula: (row width - open canopy distance)/row width. Canopy closure measurements were a critical component used in this study to determine the daily amount of irrigation water. Rainfall was not accounted for in the equation because calculated pan evaporation values on cloudy and rainy days decreased considerably, so this method accounts for rainfall without directly subtracting rainfall from irrigation amounts. Thus, irrigation was applied based on the above equation, replacing water lost by evapotranspiration. Soil moisture sensors were not used as a tool to schedule irrigation with pan evaporation plots. However, some sensors were installed in the pan evaporation plots in order to monitor soil water dynamics.


#### 2.3.2. Sensor-Based Irrigation Scheduling

Irrigation events in the sensor-based irrigation scheduling plots were conducted using information from the soil water tension sensors installed between two corn plants in the experimental area. The smart sensor array system used in our research consisted of a centrally located receiver connected to a laptop computer and multiple sensor nodes installed in the field. Each sensor node consisted of sensors (3 soil moisture sensors–Watermark ®) and thermocouples, a sensor circuit board, and an active transmitter, which transmitted data to the receiver [32]. The smart sensor board acquired sensor values and wirelessly transmitted those values to a centrally located radio frequency receiver, and then via cellphone signal, the data was transmitted to a website. The board of each node was able to read up to three Watermark ® granular resistive-type soil moisture sensors located at three soil depths (15, 30, and 60 cm) enabling a better assessment of the water availability through the soil profile. In terms of water adsorption, these three soil depths corresponded to the most active root zone in corn. A Watermark ® sensor is a granular matrix device that was used to measure soil water tension; therefore the data obtained from the Watermark ® sensors is transferred to the database in Kilopascal units. These soil water tension measuring devices provided a continuous measurement analogous to the force (soil matric potential) necessary for corn plants to extract water from the soil and these values were used to schedule irrigation. Watermark ® sensors have been used to measure soil water status for irrigation management and other purposes for more than two decades [33–36].


#### 2.3.3. Soil Water Dynamics and Irrigation Strategy

Though soil moisture sensors provided hourly data of soil matric potential in the field, the next important step was to set the Management Allowed Depletion (MAD) which is how much water to decrease from the soil before starting irrigation. The MAD represented the level of plant available water (PAW) used by the plant or evaporated before irrigation was applied without exposing the plants to water stress. Several studies have recommended irrigating row crops such as corn or cotton (Gossypium hirsutum L.) when the MAD approaches 50% [37–39]. On the other hand, Irmak et al. [40] used 35% depletion of PAW in a study testing irrigation strategies based on soil matric potential sensors. Therefore, as it was the first time that this type of experiment was conducted in Alabama, a MAD of 35% was selected as a more conservative approach. Once the MAD in volumetric water content units was identified, it was necessary to identify the soil matric potential that corresponded with the specific MAD value. In order to define this MAD, soil water level, it was necessary to estimate first the field capacity (FC), permanent wilting point (PWP), and PAW values, which were important parameters to determine the soil moisture threshold to start irrigation on the sensor-based irrigation plots. Soil texture data from the soil samples collected from the study site was used to calculate a pedotransfer function which was used to estimate soil water retention curves (SWRC) and subsequently the FC and PWP parameters. Saxton et al. [41] stated that experiences have shown that soil texture predominately determines the water-holding characteristics of most agricultural soils.Disturbed soil samples were collected from three depths (15, 30, and 60 cm) in a representative area of the field. The disturbed soil samples were used to run a soil texture analysis, which allowed for estimation of the percentage of sand, silt, and clay for each sample. These data were used to estimate the van Genuchten equation parameters using the retention curve (RECT) computer program. The percent of sand, silt, and clay data were input in the (RETC) computer program; thus this software defined theθr, θs, α, and n values representing the experimental field. All these variables were plugged in the van Genuchten equation. The van Genuchten model has the following form:(1)θh=θr+θs−θr1+αhnm,where θ(h) is the actual soil water content (cm3 cm−3) at the suction h (cm, taken positive for increasing suctions); θr and θs are the residual and saturated soil water contents (cm3 cm−3), respectively; α is a parameter related to the inverse of the air entry suction (cm−1); and m and n are curve shape parameters [42]. Note that m here characterizes the asymmetry using the constraint m = 1 −−1/n; therefore the water retention characteristics defined by this equation only contained four unknown parameters that were θr, θs, α, and n and these unknown parameters were defined using the (RETC) computer program (Table 1). The solution of this equation provided enough data to generate the SWRC representing the soil water content according to the soil matric potential, which was essential to set the irrigation threshold (Figure 1).Table 1
van Genuchten equation parameters from 2014, percentage (sand, silt, and clay), and soil textural class for three soil depths at the TVREC. The parameters from 2014 were used during both growing seasons.
Depth (cm)0–1515–3030–60θrǂ0.08080.09690.1048θs§0.40950.47790.536α¶0.0890.01320.0178N£1.46691.3651.3098m¤0.31830.26740.2365θs − θr¥0.32870.3810.4312% sand16.18109.76% silt51.0543.4535.74% clay32.7746.5554.5Textural classSilty clay loamSilty clayClayǂ Residual soil water content (cm3 cm−3). §Saturated soil water content (cm3 cm−3). ¶Parameter related to the inverse of the air entry suction (cm−1); £n curve shape parameters (empirical parameter). ¤m = (1 − 1/n). ¥θr and θs are the residual and saturated soil water contents (cm3 cm−3).Figure 1
Soil water retention curves created using the van Genuchten equation for each depth (15, 30, and 60 cm) to convert the soil water content data into soil matric potential in 2014 and 2015 at TVREC.Each SWRC was generated in order to estimate the amount of water retained in a soil (expressed as mass or volume water content) under equilibrium at a given soil matric potential. In practical applications, soil matric potential has a negative sign (i.e., more negative soil matric potential values indicate drier soil); however, a positive sign is used in this study as an indicator of soil water tension.Once the SWRC was generated, the FC and PWP were identified. These parameters are the key for the PAW estimation. Richards and Weaver [43] found that a coarse-textured soil can reach FC at 10 kPa and a fine-textured soil can reach FC at 33 kPa. The PWP is usually found at a soil matric potential of 1500 kPa. Therefore, for the purpose of this study, soil matric potential values representing the FC were selected based on soil texture. In most cases, a FC value at 33 kPa was used in this study when the soil textural class in the experimental field was silty clay loam, silty clay, or clay (Tables 1 and 2).Table 2
Field capacity, permanent wilting point, plant available water, water depletion, and 35% depletion of plant available water in water content and soil matric potential units for each soil depth at TVREC in 2014. Parameters from 2014 were used during both growing seasons.
Depth (cm)0–1515–3030–60Field capacity (FC)ǂ0.14810.31150.3478Permanent wilting point (PWP)§0.09210.15170.1806Plant available water (PAW)¶0.0560.15980.1672Water depletion (35%) (cm3/cm3)£0.01960.05590.058535% of PAW (cm3/cm3)¤0.12850.25550.289235 % of PAW (kPa)¥708083ǂ(cm3 cm3), §(cm3 cm3), ¶PAW = FC − PWP, £water depletion (35%) = PAW, ∗0.35–35% was the soil water depletion selected in the study, ¤35% of PAW (cm3 cm3) = PAW − water depletion (35%), and ¥35% of PAW (kPa) = soil matric potential referred to the depletion of 35% of PAW.The FC and PWP were established; thus the plant available water (PAW) was calculated using FC minus PWP (Table2). Once PAW was calculated, the irrigation threshold was set using the MAD value of 35%. The specific soil matric potential (tension) representing the irrigation threshold was equal to 76 kPa, which was the average value of the three soil depths (Table 2). Every time the sensors installed on the sensor-based irrigation plots were approaching 76 kPa, an irrigation event was triggered. A weighted average of the sensor values located at the three different depths was used and it changed according to the corn growth stage. At the beginning of the season, the upper sensor (15 cm depth) had more importance as a parameter to schedule irrigation events because corn roots have not reached the deeper layers of the soil. However, as the season progressed, the average of the three sensor tension values was used. Dates and amounts of irrigation events for both pan evaporation method and sensor-based irrigation plots were registered; also each irrigation event did not exceed 15 to 25 mm in order to avoid runoff. Sensors located in the pan evaporation plots were also monitored; however, the irrigation trigger was not based on sensor values but on the method described previously.


### 2.4. Nitrogen Rate Management Treatments

Three N rates (202, 269, and 336 kg N ha−1) and the control (0 kg N ha−1) were evaluated under the two irrigation methods in order to observe the impact of different N rates and irrigation practices including their interaction on corn yield, N uptake, and NUE. Aboveground biomass samples were collected at harvest to determine N concentration in the grain and stover and then establish the impact of different N rates on NUE.


### 2.5. Field Data Collection and Data Analysis

At maturity, corn plants were harvested from two 1 m row sections of each plot. Samples were harvested from rows 1 and 4 because the middle two rows were harvested for yield data using the plot combine. Harvested plants were separated in ears and stover (stem + leaves); then the stover was weighted at the field and a subsample was collected. The ears and stover subsamples were oven-dried at 60°C for one week to reach a constant weight. Also, the ears were shelled prior to N analyses. Posteriorly, the grain, cob, and stover were ground through a 2 mm sieve and analyzed for total C and N with the LECOR C/N analyzer (LECO Corp., St. Joseph, MI).The two center rows were harvested using a plot combine, and the grain yield was based on the data provided by the combine. The stover yield was the sum of cob and stalk yield. Stalk yield was calculated based on the fresh stalk yield harvested in 2 m row and the proportional oven-dried stalk yield in the subsample. Cob yield was calculated based on total weight of oven-dried cobs after the grain was shelled. Grain yield, aboveground biomass, plant N concentration, and plant N uptake are reported on an oven-dried basis.For each plot, the aboveground N uptake, recovery efficiency (REn), nitrogen agronomic efficiency (AEn), nitrogen internal efficiency (IEn), and nitrogen partial productivity (PFPn) values were calculated using the following equations:(2)aboveground N uptake=grain N concentration×grain yield+Stover N concentration×Stover biomass,REN=UN−U0N rate×100%,AEN=YN−Y0N rate,IEN=YNUN,PFPN=YNN rate,where UN and YN represented the N uptake by aboveground biomass and grain yield in the fertilized treatments, respectively, and U0 and Y0 represented the N uptake by aboveground biomass and grain yield in the control treatments, respectively. The N rate was the total amount of N applied during the corn growing season. The NUE was calculated using four different formulas in order to evaluate the interaction between yield and N uptake as N rate change.


### 2.6. Statistical Analysis

Yield, aboveground biomass, grain N concentration, stover N concentration, grain N uptake, aboveground N uptake, REn, AEn, IEn, and PFPn differences among N rates were analyzed for both years of the study and for each irrigation scenario (pan evaporation irrigation scheduling method and sensor-based irrigation scheduling). The statistical analysis was conducted using the procedure for generalized linear mixed models (PROC GLIMMIX) implemented in SAS 9.1 (SAS for Windows v. 9.1, SAS Institute Inc., Cary, NC). All dependent variables were first analyzed to measure the effect of the irrigation scheduling methods, N rates, and years as fixed effects. It allowed us to measure the magnitude of irrigation method × N rate × years interactions. Preliminary analyses suggested the need for treating years as nonfixed effect in order to observe the effects of each year independently and scheduling methods and N rates as fixed effects. The mean separation between irrigation methods, N rates, years, and all the interactions were obtained by Tukey’s significant difference test (P<0.05).


## 2.1. Site Description

The field experiment was conducted at Tennessee Valley Research and Extension Center located in Belle Mina, Alabama (34°39′24″N 86°52′45″W, 183 m above mean sea level) during the 2014 and 2015 growing seasons. The dominant soil series at the research site is a Decatur silty clay loam (fine, kaolinitic, thermic Rhodic Paleudults) with 0% to 10% slope [28]; rainfall conditions during both growing seasons (April–August) were below the 30-year average (1971–2000) according to the experimental station data. Total rainfall during the 2014 growing season was 11 mm below average and 19 mm below average in 2015. However, in 2014, rainfall distribution per month with respect to the historic average (1971–2000) was higher in April (+39 mm) and June (+42 mm), but rainfall in May (−50 mm), July (−13 mm), and August (−29 mm) was below average. In 2015, above normal rainfall with respect to the 30-year average was observed in April (+21 mm) and May (+6 mm) contrasting with rainfall in June (−21 mm), July (−1 mm), and August (−24 mm) that was below average.


## 2.2. Experimental Design

A split-plot design with three replications was implemented for this study in which irrigation regime was the main plot and N fertilizer rates were the subplots. Two irrigation scheduling methods (pan evaporation method and sensor-based irrigation method), and four N rates (0-control, 202, 269, and 336 kg ha−1) were tested for plant N uptake, NUE, and yield differences. Plots consisted of eight rows 11.8 m long and 5.3 m wide with a row spacing of 76.2 cm and the seeding rate was 89,000 plants ha−1. Urea ammonium nitrate (UAN 32-0-0) was used as the fertilizer source, which was split applied with 1/3 of the total applied at planting and the remainder in a side-dress application at the V6 growth stage. Each eight-row plot was divided in half, 4 rows received 202, 269, or 336 kg N ha−1 depending on the randomization, and the other adjacent 4 rows did not receive any N which allowed for NUE estimation. Plots received 224.1 kg ha−1 in a preplant application of P and K. All the cultural practices were performed according to Alabama Cooperative Extension System (ACES) recommendation to maximize corn yield. Each plot was independently irrigated with four overhead sprinkler nozzles located in each corner of the plot.According to ACES, the recommended N rate for irrigated corn in the state is 224.1 kg ha−1 in order to achieve 11299 kg ha−1 yield. If the yield goal is greater than 12555 kg ha−1, the N recommendation is to apply 1 kg of N per 45.3 kg of anticipated corn yield [29]. The Pioneer 1690YHR corn hybrid was planted on April 4th in 2014 and April 24th in 2015 in single rows and managed with conservation tillage and rye (Secale cereale) cover crop. In both growing seasons, the cover crop was planted on 11/21/13 and 11/07/14, respectively. It received 33.6 kg ha−1 of N fertilizer (UAN 32-0-0) in both seasons to maximize biomass production.


## 2.3. Water Management Treatments

Two irrigation scheduling methods tested in our study differed from each other in the way to estimate the amount of water stored in the soil during the growing season. The pan evaporation method was based on calculating the crop’s estimated evapotranspiration (ET) using the pan evaporation values and the crop’s consumptive water use. On the other hand, the sensor-based irrigation relied on measurements of soil matric potential provided by sensors in order to apply the correct amount of irrigation at the right time. The water use efficiency (WUE) of the two irrigation scheduling methods was included in our evaluation to establish differences in the total amount of water applied in the plots.


### 2.3.1. Pan Evaporation Irrigation Scheduling

Daily pan evaporation data was reported from the weather station located at the research site (http://www.awis.com). The daily irrigation requirement, ET, was calculated using the following equation: ET = 0.90 × PAN × CC where evapotranspiration, ET (mm day−1), was calculated from 90% pan evaporation, PAN (mm day−1) adjusted for fractional canopy cover (CC) [30, 31]. Canopy cover was determined weekly by measuring the open canopy distance (cm) between rows with a tape measure. Fractional canopy cover was calculated using the following formula: (row width - open canopy distance)/row width. Canopy closure measurements were a critical component used in this study to determine the daily amount of irrigation water. Rainfall was not accounted for in the equation because calculated pan evaporation values on cloudy and rainy days decreased considerably, so this method accounts for rainfall without directly subtracting rainfall from irrigation amounts. Thus, irrigation was applied based on the above equation, replacing water lost by evapotranspiration. Soil moisture sensors were not used as a tool to schedule irrigation with pan evaporation plots. However, some sensors were installed in the pan evaporation plots in order to monitor soil water dynamics.


### 2.3.2. Sensor-Based Irrigation Scheduling

Irrigation events in the sensor-based irrigation scheduling plots were conducted using information from the soil water tension sensors installed between two corn plants in the experimental area. The smart sensor array system used in our research consisted of a centrally located receiver connected to a laptop computer and multiple sensor nodes installed in the field. Each sensor node consisted of sensors (3 soil moisture sensors–Watermark ®) and thermocouples, a sensor circuit board, and an active transmitter, which transmitted data to the receiver [32]. The smart sensor board acquired sensor values and wirelessly transmitted those values to a centrally located radio frequency receiver, and then via cellphone signal, the data was transmitted to a website. The board of each node was able to read up to three Watermark ® granular resistive-type soil moisture sensors located at three soil depths (15, 30, and 60 cm) enabling a better assessment of the water availability through the soil profile. In terms of water adsorption, these three soil depths corresponded to the most active root zone in corn. A Watermark ® sensor is a granular matrix device that was used to measure soil water tension; therefore the data obtained from the Watermark ® sensors is transferred to the database in Kilopascal units. These soil water tension measuring devices provided a continuous measurement analogous to the force (soil matric potential) necessary for corn plants to extract water from the soil and these values were used to schedule irrigation. Watermark ® sensors have been used to measure soil water status for irrigation management and other purposes for more than two decades [33–36].


### 2.3.3. Soil Water Dynamics and Irrigation Strategy

Though soil moisture sensors provided hourly data of soil matric potential in the field, the next important step was to set the Management Allowed Depletion (MAD) which is how much water to decrease from the soil before starting irrigation. The MAD represented the level of plant available water (PAW) used by the plant or evaporated before irrigation was applied without exposing the plants to water stress. Several studies have recommended irrigating row crops such as corn or cotton (Gossypium hirsutum L.) when the MAD approaches 50% [37–39]. On the other hand, Irmak et al. [40] used 35% depletion of PAW in a study testing irrigation strategies based on soil matric potential sensors. Therefore, as it was the first time that this type of experiment was conducted in Alabama, a MAD of 35% was selected as a more conservative approach. Once the MAD in volumetric water content units was identified, it was necessary to identify the soil matric potential that corresponded with the specific MAD value. In order to define this MAD, soil water level, it was necessary to estimate first the field capacity (FC), permanent wilting point (PWP), and PAW values, which were important parameters to determine the soil moisture threshold to start irrigation on the sensor-based irrigation plots. Soil texture data from the soil samples collected from the study site was used to calculate a pedotransfer function which was used to estimate soil water retention curves (SWRC) and subsequently the FC and PWP parameters. Saxton et al. [41] stated that experiences have shown that soil texture predominately determines the water-holding characteristics of most agricultural soils.Disturbed soil samples were collected from three depths (15, 30, and 60 cm) in a representative area of the field. The disturbed soil samples were used to run a soil texture analysis, which allowed for estimation of the percentage of sand, silt, and clay for each sample. These data were used to estimate the van Genuchten equation parameters using the retention curve (RECT) computer program. The percent of sand, silt, and clay data were input in the (RETC) computer program; thus this software defined theθr, θs, α, and n values representing the experimental field. All these variables were plugged in the van Genuchten equation. The van Genuchten model has the following form:(1)θh=θr+θs−θr1+αhnm,where θ(h) is the actual soil water content (cm3 cm−3) at the suction h (cm, taken positive for increasing suctions); θr and θs are the residual and saturated soil water contents (cm3 cm−3), respectively; α is a parameter related to the inverse of the air entry suction (cm−1); and m and n are curve shape parameters [42]. Note that m here characterizes the asymmetry using the constraint m = 1 −−1/n; therefore the water retention characteristics defined by this equation only contained four unknown parameters that were θr, θs, α, and n and these unknown parameters were defined using the (RETC) computer program (Table 1). The solution of this equation provided enough data to generate the SWRC representing the soil water content according to the soil matric potential, which was essential to set the irrigation threshold (Figure 1).Table 1
van Genuchten equation parameters from 2014, percentage (sand, silt, and clay), and soil textural class for three soil depths at the TVREC. The parameters from 2014 were used during both growing seasons.
Depth (cm)0–1515–3030–60θrǂ0.08080.09690.1048θs§0.40950.47790.536α¶0.0890.01320.0178N£1.46691.3651.3098m¤0.31830.26740.2365θs − θr¥0.32870.3810.4312% sand16.18109.76% silt51.0543.4535.74% clay32.7746.5554.5Textural classSilty clay loamSilty clayClayǂ Residual soil water content (cm3 cm−3). §Saturated soil water content (cm3 cm−3). ¶Parameter related to the inverse of the air entry suction (cm−1); £n curve shape parameters (empirical parameter). ¤m = (1 − 1/n). ¥θr and θs are the residual and saturated soil water contents (cm3 cm−3).Figure 1
Soil water retention curves created using the van Genuchten equation for each depth (15, 30, and 60 cm) to convert the soil water content data into soil matric potential in 2014 and 2015 at TVREC.Each SWRC was generated in order to estimate the amount of water retained in a soil (expressed as mass or volume water content) under equilibrium at a given soil matric potential. In practical applications, soil matric potential has a negative sign (i.e., more negative soil matric potential values indicate drier soil); however, a positive sign is used in this study as an indicator of soil water tension.Once the SWRC was generated, the FC and PWP were identified. These parameters are the key for the PAW estimation. Richards and Weaver [43] found that a coarse-textured soil can reach FC at 10 kPa and a fine-textured soil can reach FC at 33 kPa. The PWP is usually found at a soil matric potential of 1500 kPa. Therefore, for the purpose of this study, soil matric potential values representing the FC were selected based on soil texture. In most cases, a FC value at 33 kPa was used in this study when the soil textural class in the experimental field was silty clay loam, silty clay, or clay (Tables 1 and 2).Table 2
Field capacity, permanent wilting point, plant available water, water depletion, and 35% depletion of plant available water in water content and soil matric potential units for each soil depth at TVREC in 2014. Parameters from 2014 were used during both growing seasons.
Depth (cm)0–1515–3030–60Field capacity (FC)ǂ0.14810.31150.3478Permanent wilting point (PWP)§0.09210.15170.1806Plant available water (PAW)¶0.0560.15980.1672Water depletion (35%) (cm3/cm3)£0.01960.05590.058535% of PAW (cm3/cm3)¤0.12850.25550.289235 % of PAW (kPa)¥708083ǂ(cm3 cm3), §(cm3 cm3), ¶PAW = FC − PWP, £water depletion (35%) = PAW, ∗0.35–35% was the soil water depletion selected in the study, ¤35% of PAW (cm3 cm3) = PAW − water depletion (35%), and ¥35% of PAW (kPa) = soil matric potential referred to the depletion of 35% of PAW.The FC and PWP were established; thus the plant available water (PAW) was calculated using FC minus PWP (Table2). Once PAW was calculated, the irrigation threshold was set using the MAD value of 35%. The specific soil matric potential (tension) representing the irrigation threshold was equal to 76 kPa, which was the average value of the three soil depths (Table 2). Every time the sensors installed on the sensor-based irrigation plots were approaching 76 kPa, an irrigation event was triggered. A weighted average of the sensor values located at the three different depths was used and it changed according to the corn growth stage. At the beginning of the season, the upper sensor (15 cm depth) had more importance as a parameter to schedule irrigation events because corn roots have not reached the deeper layers of the soil. However, as the season progressed, the average of the three sensor tension values was used. Dates and amounts of irrigation events for both pan evaporation method and sensor-based irrigation plots were registered; also each irrigation event did not exceed 15 to 25 mm in order to avoid runoff. Sensors located in the pan evaporation plots were also monitored; however, the irrigation trigger was not based on sensor values but on the method described previously.


## 2.3.1. Pan Evaporation Irrigation Scheduling

Daily pan evaporation data was reported from the weather station located at the research site (http://www.awis.com). The daily irrigation requirement, ET, was calculated using the following equation: ET = 0.90 × PAN × CC where evapotranspiration, ET (mm day−1), was calculated from 90% pan evaporation, PAN (mm day−1) adjusted for fractional canopy cover (CC) [30, 31]. Canopy cover was determined weekly by measuring the open canopy distance (cm) between rows with a tape measure. Fractional canopy cover was calculated using the following formula: (row width - open canopy distance)/row width. Canopy closure measurements were a critical component used in this study to determine the daily amount of irrigation water. Rainfall was not accounted for in the equation because calculated pan evaporation values on cloudy and rainy days decreased considerably, so this method accounts for rainfall without directly subtracting rainfall from irrigation amounts. Thus, irrigation was applied based on the above equation, replacing water lost by evapotranspiration. Soil moisture sensors were not used as a tool to schedule irrigation with pan evaporation plots. However, some sensors were installed in the pan evaporation plots in order to monitor soil water dynamics.


## 2.3.2. Sensor-Based Irrigation Scheduling

Irrigation events in the sensor-based irrigation scheduling plots were conducted using information from the soil water tension sensors installed between two corn plants in the experimental area. The smart sensor array system used in our research consisted of a centrally located receiver connected to a laptop computer and multiple sensor nodes installed in the field. Each sensor node consisted of sensors (3 soil moisture sensors–Watermark ®) and thermocouples, a sensor circuit board, and an active transmitter, which transmitted data to the receiver [32]. The smart sensor board acquired sensor values and wirelessly transmitted those values to a centrally located radio frequency receiver, and then via cellphone signal, the data was transmitted to a website. The board of each node was able to read up to three Watermark ® granular resistive-type soil moisture sensors located at three soil depths (15, 30, and 60 cm) enabling a better assessment of the water availability through the soil profile. In terms of water adsorption, these three soil depths corresponded to the most active root zone in corn. A Watermark ® sensor is a granular matrix device that was used to measure soil water tension; therefore the data obtained from the Watermark ® sensors is transferred to the database in Kilopascal units. These soil water tension measuring devices provided a continuous measurement analogous to the force (soil matric potential) necessary for corn plants to extract water from the soil and these values were used to schedule irrigation. Watermark ® sensors have been used to measure soil water status for irrigation management and other purposes for more than two decades [33–36].


## 2.3.3. Soil Water Dynamics and Irrigation Strategy

Though soil moisture sensors provided hourly data of soil matric potential in the field, the next important step was to set the Management Allowed Depletion (MAD) which is how much water to decrease from the soil before starting irrigation. The MAD represented the level of plant available water (PAW) used by the plant or evaporated before irrigation was applied without exposing the plants to water stress. Several studies have recommended irrigating row crops such as corn or cotton (Gossypium hirsutum L.) when the MAD approaches 50% [37–39]. On the other hand, Irmak et al. [40] used 35% depletion of PAW in a study testing irrigation strategies based on soil matric potential sensors. Therefore, as it was the first time that this type of experiment was conducted in Alabama, a MAD of 35% was selected as a more conservative approach. Once the MAD in volumetric water content units was identified, it was necessary to identify the soil matric potential that corresponded with the specific MAD value. In order to define this MAD, soil water level, it was necessary to estimate first the field capacity (FC), permanent wilting point (PWP), and PAW values, which were important parameters to determine the soil moisture threshold to start irrigation on the sensor-based irrigation plots. Soil texture data from the soil samples collected from the study site was used to calculate a pedotransfer function which was used to estimate soil water retention curves (SWRC) and subsequently the FC and PWP parameters. Saxton et al. [41] stated that experiences have shown that soil texture predominately determines the water-holding characteristics of most agricultural soils.Disturbed soil samples were collected from three depths (15, 30, and 60 cm) in a representative area of the field. The disturbed soil samples were used to run a soil texture analysis, which allowed for estimation of the percentage of sand, silt, and clay for each sample. These data were used to estimate the van Genuchten equation parameters using the retention curve (RECT) computer program. The percent of sand, silt, and clay data were input in the (RETC) computer program; thus this software defined theθr, θs, α, and n values representing the experimental field. All these variables were plugged in the van Genuchten equation. The van Genuchten model has the following form:(1)θh=θr+θs−θr1+αhnm,where θ(h) is the actual soil water content (cm3 cm−3) at the suction h (cm, taken positive for increasing suctions); θr and θs are the residual and saturated soil water contents (cm3 cm−3), respectively; α is a parameter related to the inverse of the air entry suction (cm−1); and m and n are curve shape parameters [42]. Note that m here characterizes the asymmetry using the constraint m = 1 −−1/n; therefore the water retention characteristics defined by this equation only contained four unknown parameters that were θr, θs, α, and n and these unknown parameters were defined using the (RETC) computer program (Table 1). The solution of this equation provided enough data to generate the SWRC representing the soil water content according to the soil matric potential, which was essential to set the irrigation threshold (Figure 1).Table 1
van Genuchten equation parameters from 2014, percentage (sand, silt, and clay), and soil textural class for three soil depths at the TVREC. The parameters from 2014 were used during both growing seasons.
Depth (cm)0–1515–3030–60θrǂ0.08080.09690.1048θs§0.40950.47790.536α¶0.0890.01320.0178N£1.46691.3651.3098m¤0.31830.26740.2365θs − θr¥0.32870.3810.4312% sand16.18109.76% silt51.0543.4535.74% clay32.7746.5554.5Textural classSilty clay loamSilty clayClayǂ Residual soil water content (cm3 cm−3). §Saturated soil water content (cm3 cm−3). ¶Parameter related to the inverse of the air entry suction (cm−1); £n curve shape parameters (empirical parameter). ¤m = (1 − 1/n). ¥θr and θs are the residual and saturated soil water contents (cm3 cm−3).Figure 1
Soil water retention curves created using the van Genuchten equation for each depth (15, 30, and 60 cm) to convert the soil water content data into soil matric potential in 2014 and 2015 at TVREC.Each SWRC was generated in order to estimate the amount of water retained in a soil (expressed as mass or volume water content) under equilibrium at a given soil matric potential. In practical applications, soil matric potential has a negative sign (i.e., more negative soil matric potential values indicate drier soil); however, a positive sign is used in this study as an indicator of soil water tension.Once the SWRC was generated, the FC and PWP were identified. These parameters are the key for the PAW estimation. Richards and Weaver [43] found that a coarse-textured soil can reach FC at 10 kPa and a fine-textured soil can reach FC at 33 kPa. The PWP is usually found at a soil matric potential of 1500 kPa. Therefore, for the purpose of this study, soil matric potential values representing the FC were selected based on soil texture. In most cases, a FC value at 33 kPa was used in this study when the soil textural class in the experimental field was silty clay loam, silty clay, or clay (Tables 1 and 2).Table 2
Field capacity, permanent wilting point, plant available water, water depletion, and 35% depletion of plant available water in water content and soil matric potential units for each soil depth at TVREC in 2014. Parameters from 2014 were used during both growing seasons.
Depth (cm)0–1515–3030–60Field capacity (FC)ǂ0.14810.31150.3478Permanent wilting point (PWP)§0.09210.15170.1806Plant available water (PAW)¶0.0560.15980.1672Water depletion (35%) (cm3/cm3)£0.01960.05590.058535% of PAW (cm3/cm3)¤0.12850.25550.289235 % of PAW (kPa)¥708083ǂ(cm3 cm3), §(cm3 cm3), ¶PAW = FC − PWP, £water depletion (35%) = PAW, ∗0.35–35% was the soil water depletion selected in the study, ¤35% of PAW (cm3 cm3) = PAW − water depletion (35%), and ¥35% of PAW (kPa) = soil matric potential referred to the depletion of 35% of PAW.The FC and PWP were established; thus the plant available water (PAW) was calculated using FC minus PWP (Table2). Once PAW was calculated, the irrigation threshold was set using the MAD value of 35%. The specific soil matric potential (tension) representing the irrigation threshold was equal to 76 kPa, which was the average value of the three soil depths (Table 2). Every time the sensors installed on the sensor-based irrigation plots were approaching 76 kPa, an irrigation event was triggered. A weighted average of the sensor values located at the three different depths was used and it changed according to the corn growth stage. At the beginning of the season, the upper sensor (15 cm depth) had more importance as a parameter to schedule irrigation events because corn roots have not reached the deeper layers of the soil. However, as the season progressed, the average of the three sensor tension values was used. Dates and amounts of irrigation events for both pan evaporation method and sensor-based irrigation plots were registered; also each irrigation event did not exceed 15 to 25 mm in order to avoid runoff. Sensors located in the pan evaporation plots were also monitored; however, the irrigation trigger was not based on sensor values but on the method described previously.


## 2.4. Nitrogen Rate Management Treatments

Three N rates (202, 269, and 336 kg N ha−1) and the control (0 kg N ha−1) were evaluated under the two irrigation methods in order to observe the impact of different N rates and irrigation practices including their interaction on corn yield, N uptake, and NUE. Aboveground biomass samples were collected at harvest to determine N concentration in the grain and stover and then establish the impact of different N rates on NUE.


## 2.5. Field Data Collection and Data Analysis

At maturity, corn plants were harvested from two 1 m row sections of each plot. Samples were harvested from rows 1 and 4 because the middle two rows were harvested for yield data using the plot combine. Harvested plants were separated in ears and stover (stem + leaves); then the stover was weighted at the field and a subsample was collected. The ears and stover subsamples were oven-dried at 60°C for one week to reach a constant weight. Also, the ears were shelled prior to N analyses. Posteriorly, the grain, cob, and stover were ground through a 2 mm sieve and analyzed for total C and N with the LECOR C/N analyzer (LECO Corp., St. Joseph, MI).The two center rows were harvested using a plot combine, and the grain yield was based on the data provided by the combine. The stover yield was the sum of cob and stalk yield. Stalk yield was calculated based on the fresh stalk yield harvested in 2 m row and the proportional oven-dried stalk yield in the subsample. Cob yield was calculated based on total weight of oven-dried cobs after the grain was shelled. Grain yield, aboveground biomass, plant N concentration, and plant N uptake are reported on an oven-dried basis.For each plot, the aboveground N uptake, recovery efficiency (REn), nitrogen agronomic efficiency (AEn), nitrogen internal efficiency (IEn), and nitrogen partial productivity (PFPn) values were calculated using the following equations:(2)aboveground N uptake=grain N concentration×grain yield+Stover N concentration×Stover biomass,REN=UN−U0N rate×100%,AEN=YN−Y0N rate,IEN=YNUN,PFPN=YNN rate,where UN and YN represented the N uptake by aboveground biomass and grain yield in the fertilized treatments, respectively, and U0 and Y0 represented the N uptake by aboveground biomass and grain yield in the control treatments, respectively. The N rate was the total amount of N applied during the corn growing season. The NUE was calculated using four different formulas in order to evaluate the interaction between yield and N uptake as N rate change.


## 2.6. Statistical Analysis

Yield, aboveground biomass, grain N concentration, stover N concentration, grain N uptake, aboveground N uptake, REn, AEn, IEn, and PFPn differences among N rates were analyzed for both years of the study and for each irrigation scenario (pan evaporation irrigation scheduling method and sensor-based irrigation scheduling). The statistical analysis was conducted using the procedure for generalized linear mixed models (PROC GLIMMIX) implemented in SAS 9.1 (SAS for Windows v. 9.1, SAS Institute Inc., Cary, NC). All dependent variables were first analyzed to measure the effect of the irrigation scheduling methods, N rates, and years as fixed effects. It allowed us to measure the magnitude of irrigation method × N rate × years interactions. Preliminary analyses suggested the need for treating years as nonfixed effect in order to observe the effects of each year independently and scheduling methods and N rates as fixed effects. The mean separation between irrigation methods, N rates, years, and all the interactions were obtained by Tukey’s significant difference test (P<0.05).


## 3. Results and Discussion

### 3.1. Climatic Conditions and Irrigation Performance

Monthly average weather variables for the 2013, 2014, and 2015 growing seasons including the long-term average values (1971–2000) are summarized in Table3. Overall, climatic conditions helped explain differences in irrigation amount and grain yield. Differences in rainfall and temperature influence crop evapotranspiration in a given environment. Runge [44] found that daily temperature and rainfall affect corn yield from 25 days before to 15 days after anthesis. The seasonal average temperature was 22.60°C and 23.50°C in 2014 and 2015, respectively. The seasonal total rainfall in 2014 and 2015, respectively, was 502 mm (2.1% below normal) and 495 mm (3.7% below normal). The total rainfall was relatively similar across the two years and slightly below the (30 years) normal, which was 513 mm. However, rainfall distribution through the growing season will have greater effect on irrigation amount than total rainfall. Thus, if rainfall is concentrated in the period of high-water demand, irrigation requirements will be reduced and if rainfall is concentrated at the beginning and the end of the season, it will also limit irrigation requirements. In this study, rainfall was distributed differently before and after tasseling in 2014 but evenly in 2015. According to Lee (2011), tasseling is the highest corn water demand period. In 2014 and 2015, 41% and 58% of total rainfall and 49% and 51% of the total rainfall were evenly distributed before and after tasseling, respectively.Table 3
Monthly average climatic conditions, during the 2013, 2014, and 2015 growing seasons and long-term average values (30 years) measured at the research site TVREC.
MonthRainfall (mm)Temperature max (°C)Temperature min (°C)Average temperature (°C)2013April482210.315.9May5925.814.516June3031.519.825.2July10630.620.625.2August1530.720.124.6Total/average25928.117.121.42014April14723.31016.9May6127.615.221.2June1513120.325July933119.324.9August5032.719.825.5Total/average50229.116.922.62015April12923.311.917.9May11727.915.621.5June8832.320.626.1July10533.12227.2August5531.219.824.9Total/average49429.617.923.5Historic average (1971–2000)April10822.28.315.3May11126.713.620.1June10930.617.924.3July10632.42026.2August7932.118.925.5Total/average51328.815.722.3


### 3.2. Irrigation Amounts

Average irrigation amounts applied in the sensor plots were 106 mm and 142 mm in 2014 and 2015, respectively. However, the temporal distribution of rainfall between the two growing seasons was not comparable. Rainfall in the 2014 growing season was concentrated after tasseling, 294 mm, and in 2015, the rainfall distribution was similar before and after tasseling, 246 mm and 248 mm, respectively. Thus, considering that after tasseling corresponds to a period of greater water demand in corn, consequently, the 2015 season received much irrigation.Irrigation amounts for the two irrigation scheduling methods indicated that the sensor-based method received less irrigation during both growing seasons. The 2014 growing season received 138 mm of irrigation water on the plots under the pan evaporation irrigation scheduling method and the sensor-based irrigation scheduling method received 122 mm in the plots fertilized with 202 kg N ha−1 and 99 mm in the plots fertilized with 269 and 336 kg N ha−1. In 2015, the pan evaporation plots received 215 mm of irrigation and the sensor-based treatments received 152 mm at 202 and 336 kg N ha−1 plots and 127 mm at 269 kg N ha−1 plots (Table 4).Table 4
Seasonal irrigation on corn under different irrigation scheduling methods and N rates in 2014 and 2015 at TVREC.
`Irrigation amount (mm)20142015Pan evaporation0138215202138215269138215336138215Sensor-based0BNISǂBNIS2021221522699912733699152ǂBased on neighboring plot irrigation scheduling (BNIS).Figures2 and 3 illustrate soil water dynamics in the pan evaporation and sensor-based plots during the 2015 growing season, respectively. Daily soil matric potential data in the figures was collected from sensors installed in the plots fertilized with 269 kg N ha−1. The figures show how irrigation and rainfall affected soil matric potential values. From the pan evaporation method, it was possible to observe that even though a specific irrigation threshold was not set, the soil matric potential was maintained below 76 kPa, at approximately 50 kPa, due to more frequent irrigation events (Figure 2). In the sensor-based approach, the irrigation events were applied when the soil matric potential approached the threshold (76 kPa) (Figure 3). Therefore, since soil matric potential was maintained at lower values in the pan evaporation method, it infers that the soil had higher moisture during the growing season. Also, the sensor-based method depleted more water from the soil for irrigation events scheduled at a higher soil matric potential. After July 22, malfunctions occurred mainly in the sensors located in the 30 and 60 cm depth in this treatment (sensor-based irrigation fertilized with 269 kg N ha−1). However, irrigation practices in the sensor-based irrigation were continued using the shallowest sensor.Figure 2
Daily soil matric potential in the pan evaporation plot fertilized with 269 kg N ha−1 in 2015 at TVREC. The amount of irrigation referred to in the black arrows was value between 15 and 25 mm to avoid runoff.Figure 3
Daily soil matric potential in the sensor-based plot fertilized with 269 kg N ha−1 in 2015 at TVREC. The amount of irrigation referred to in the black arrows was value between 15 and 25 mm to avoid runoff.


### 3.3. Irrigation Scheduling Methods Effects: Pan Evaporation and Sensor-Based Irrigation

Differences between pan evaporation and sensor-based method were not significant with respect to all parameters analyzed in this study, except for REn (Table 5). Therefore, different amounts of irrigation applied using the two irrigation scheduling methods did not result in significant interactions among treatments that include irrigation (I x N, I x Y, and I x N x Y) for examined variables during 2014 and 2015.Table 5
Summary of ANOVA for grain yield, aboveground biomass, N uptake, and NUE with respect to N rates (0, 202, 269, and 336 kg ha-1), irrigation, and years (2014–2015) at TVREC.
Sources of variationGrain yield (Mg ha−1)Aboveground biomass (Mg ha−1)Grain N concentration (g kg−1)Stover N concentration (g kg−1)Grain N uptake (g kg−1)Aboveground N uptake (kg N ha−1)Nitrogen use efficiencyAEnǂREn§IEn¶PFPn#Irrigation (I)Ns£NsNsNsNsNsNsP<0.05NsNsN rate (N)P<0.05∗NsNsNsNsP<0.05P<0.05NsNsP<0.05Year (Y)P<0.05NsP<0.05P<0.05NsP<0.05P<0.05P<0.05P<0.05P<0.05I x NNsNsNsNsNsNsNsNsNsNsY x NP<0.05NsNsNsNsP<0.05P<0.05NsNsP<0.05Y x INsNsNsNsNsNsNsP<0.05NsNsY x I x NNsNsNsNsNsNsNsNsNsNs∗Significant at 0.05 probability level. £Nonsignificant at 0.05 probability level. ǂAEn: nitrogen agronomy efficiency. §REn: recovery efficiency. ¶IEn: nitrogen internal efficiency. #PFPn: nitrogen partial productivity.Grain yield is important indicator of profitability for farmers and there were not differences between irrigation methods when these variables were analyzed (Table5). During both growing seasons, pan evaporation scheduling method on average used 29% more water than sensor-based irrigation method. In 2014, the average reduction of water applied was 30 mm and, in 2015, the average reduction of water applied was 71 mm without statistically significant differences in grain yield or aboveground biomass within each year (Table 4). Based on the 2014 USDA Crop Production Summary (http://www.usda.gov), the nonirrigated corn yield average in the state of Alabama was 9.97 Mg ha−1 and the yield average obtained from the fertilized treatments of both irrigation methods was 13.84 Mg ha−1 and 13.18 Mg ha−1 during the 2014 and 2015 growing seasons, respectively. Glass et al. [45] conducted an evaluation of corn hybrids also at the Tennessee Valley research station and the irrigated corn yield average was 13.99 Mg ha−1.


### 3.4. Effects of N Rate on Corn Grain Yield and Aboveground Biomass

Significant grain yield differences were observed with respect to N rate, year, and interaction of year × N rate (Table5). In contrast, there was no aboveground biomass differences observed with respect to the same factors (Table 5). The lowest grain yield was observed on the control treatment during both years (Table 6). In 2014, the greatest grain yield was 14.23 Mg ha−1 achieved in the plot fertilized with 269 kg ha−1 and, in 2015, the greatest yield was 13.81 Mg ha−1 achieved in the plot fertilized with 336 kg N ha−1 (Table 6). In 2014, grain yield differences between the control and fertilized treatments were statistically significant; however no statistical differences were found between the fertilized plots (Table 6). In other words, the increase of N rate from 202 to 336 kg N ha−1 did not statistically change grain yield. However, the growing season of 2015 did not follow the same trend observed in 2014 and differences between control and fertilized plots and among fertilized plots were presented.Table 6
Grain yield and aboveground biomass differences between years and N rates (0, 202, 269, and 336 kg ha−1) at the TVREC.
N rate (kg ha)Grain yield (Mg ha−1)Aboveground biomass (Mg ha−1)201420152015MeanSEǂMeanSEMeanSE07.62b§0.525.52c0.3612.56b0.7920213.53a0.5212.54b0.3615.79a0.7926914.23a0.5213.22ba0.3614.65ba0.7933613.80a0.5213.81a0.3617.18a0.79ǂStandard error. §Numbers at the same column followed by the same letter are nonsignificant at P<0.05.The aboveground biomass was not significantly affected in 2014; however in 2015, the treatments receiving no N fertilizer presented the lowest value (Table6). In 2015, the aboveground biomass showed to be statistically responsive to the increasing N rates, which is commonly found in the literature [46, 47].Hagedorn et al. [48] mentioned that, in years with heavy rains, the mineral N content in the topsoil decrease by 50–70%; however, in low rainfall seasons, the pronounced mineralization occurs but the mineral N decreases only slightly. Therefore, the low rainfall in 2013 may have contributed to the no statistical differences in the fertilized plots in 2014. The 2013 season received 50% less precipitation than the normal of 30 years normal (Table 3). Furthermore, residual mineral N may have increase because leaching was limited due to low rainfall in 2013 and it could explain the low corn yield responses to N application observed in 2014. Extra N in 2014 could also explain why control plots in 2014 obtained higher yields than 2015.The impact of N rate on grain yield is important because of its effect on production profitability. Therefore, if the N fertilizer investment does not translate into a yield increase, then excess N could cause environmental problems. Some studies in this area suggest that increasing N application above optimal rates will lead to only small increases or even decreases in the grain yield and aboveground biomass. However, it may cause an increase in the energy use and production costs as well as increasing the risk of negative environmental effects such as greenhouse gas emission, soil acidification, and nitrate leaching [49–51].


### 3.5. Variation in Corn N Concentration

Grain and stover N concentrations were different between the two growing seasons with the N concentration values in 2015 greater than the values observed in 2014 (Table7). Overall, grain and stover N concentration increased as the N rate increased (Table 7). During the two growing seasons, grain N concentration ranged from 8.75 g kg−1 to 14.00 g kg −1 and the stover N concentration ranged from 0.61 g kg−1 to 13.23 g kg−1. Overall, grain and stover N concentration increased as the N rate increased (Table 7). Ciampitti and Vyn [52] reported the N concentration range for the new era corn cultivars was 3.0 g kg−1 to 26.8 g kg and 1.2 g kg−1 to 21.1 g kg−1 for grain and stover, respectively. Setiyono et al. [53] also reported large variation of N concentration values in corn plants across several studies; therefore this may explain the oscillation of N concentration in 2014 and 2015.Table 7
Grain N, stover N concentrations and grain, and aboveground N uptake differences between years and N rates (0, 202, 269, and 336 kg ha-1) at TVREC.
N rate (kg ha−1)Grain N concentration (g kg−1)Stover N concentration (g kg−1)Aboveground N uptake (kg N ha−1)201420152014201520142015MeanSEǂMeanSEMeanSEMeanSEMeanSEMeanSE08.75b§0.6611.70c0.210.61a1.688.11b0.87125.43b22.03166.68c39.0320210.53a0.6612.52b0.213.05a1.6810.50ba0.87216.55a22.03321.68b39.0326911.10a0.6613.51a0.211.15a1.6812.16a0.87230.51a22.03351.36ba39.0333611.54a0.6614.00a0.212.45a1.6813.23a0.87241.79a22.03421.64a39.03ǂStandard error. §Numbers at the same column followed by the same letter are nonsignificant at P<0.05.In 2014, yields were higher which indicates that plants were able to use the applied N. Therefore, the N dilution could occur because higher yields and more biomass may reduce an equivalent N concentration in the plant in 2014. Moreover, temperature could be another reason of the increase of N concentration in the corn plants. Wienhold et al. [54] compared N concentration in corn and obtained 75% higher N content values in years when the temperature was above the 30-year average. Therefore, higher temperature in 2015 (Table 3) could explain the increase of N concentration for this growing season.In 2015, grain N concentration differences between the control and fertilized treatments were smaller than the differences observed in the 2014 growing season. Therefore, in 2015 the corn plants uptake more N from the unfertilized plots, showing the ability of corn to uptake nitrogen from the soil N pools even when no N fertilizer is applied (Table7). In 2014, N rate did not affect the N concentration in the stover, but in 2015, stover N concentration increased following increment of nitrogen fertilizer.Aboveground N uptake ranged from 125.43 kg N ha−1 to 421.64 kg N ha−1 during the two growing seasons. Ciampitti and Vyn [49] reported the mean, minimum, and maximum N uptake of new era maize cultivars and the values were 184, 2, and 427 kg N ha−1, respectively. Other numbers reported by Setiyono et al. [53] showed 239, 7, and 471 kg N ha−1 as mean, minimum, and maximum, respectively. Therefore, the maximum and minimum N uptake values obtained in Alabama’s environment in this present study agreed with the range of new era cultivars reported by other authors. Aboveground N uptake was higher in 2015 when compared to 2014. In 2014, the aboveground N uptake between control and fertilized treatments was statistically significant; however considering only the fertilized treatments, the N uptake values increased numerically but these differences were not significant (Table 7). During the season of 2015, the aboveground N uptake increased significantly with the increase of N rate (Table 7). Abbasi et al. [55] reported that aboveground N uptake increased with N rate; therefore in our study, the N uptake increased as the N rate increased but the increment was not significant in all cases.Setiyono et al. [53] and Chua et al. [56] reported that aboveground N uptake has been closely linked with grain and aboveground biomass. However, we observed that, in general, grain yield was higher in 2014 but the overall aboveground N uptake was lower when comparing with 2015. Herrmann and Taube [57] found lower N uptake in higher corn biomass yield and it might be due the N dilution because the N is spread out in the plant tissue.Furthermore, other studies evidenced that aboveground N uptake is not proportionately linked with grain yield. Jokela and Randall [58] presented some data from a three-year study (1982–1984) comparing grain yield and aboveground N uptake and found years with higher yield and lower N uptake or vice versa. Uribelarrea et al. [59] compared four commercial corn hybrids in three years and also found higher yield and lower N uptake; furthermore, some hybrids had a negative N balance removing more N than was applied as fertilizer-N. It also occurred in our study in the 2015 growing season (Table 7).


### 3.6. Nitrogen Use Efficiency

The Nitrogen Use Efficiency (NUE) increase indicates that most N fertilizer applied in the soil was taken up by plants and converted into grain yield. Therefore, the amount of N that remains in the soil is lower decreasing the potential of leaching and can effectively decrease the potential negative impacts of N fertilizer to the environment [55, 60].The AEnand PFPn followed the same trend in 2014 and 2015; higher AEn and PFPn were observed for the treatment fertilized with 202 kg N ha−1 and gradually decreased as N rate increased (Table 8). Also, the lowest NUE value was observed from the plots fertilized with the highest N rate (336 kg N ha-1). According to Ciampitti and Vyn [60], it is well accepted the NUEs values are high at low N rates and decrease with increasing N rates. Furthermore, AEn and PFPn may be more practical parameters because grain yield and N rate are used in the equation; therefore they are easier parameters to be calculated by farmers. The REn and IEn did not present significant differences during both seasons. Also considering that these parameters use aboveground N uptake in the equation, it might be more laborious for the farmers because it is necessary to do the N concentration analysis.Table 8
Nitrogen recovery efficiency (REn), nitrogen agronomic efficiency (AEn), nitrogen internal efficiency (IEn), and nitrogen partial productivity (PFPn) at four N rates (0, 202, 269, and 280 kg N ha−1) in 2014 and 2015 at TVREC.
N rate (kg ha−1)REn (%)ǂAEn (kg kg−1)§IEn (kg kg−1)¶PFPn (kg kg−1)#20142015201420152014201520142015MeanSEǂMeanSEMeanSEMeanSEMeanSEMeanSEMeanSEMeanSE0¤¤¤¤¤¤¤¤62.57a2.3835.04 a2.16¤¤¤¤20236.94a¥4.4498.03a16.7924.96a2.7936.32a0.6763.64a2.3840.38a2.1667.96a2.5962.33a1.8226943.96a4.4480.27a16.7922.7a2.7928.79b0.6762.76a2.3838.16a2.1652.32b2.5949.32b1.8233635.61a4.4485.58a16.7916.78b2.7923.62c0.6758.37a2.3833.94a2.1639.56c2.5941.22c1.82ǂRen: recovery efficiency. §AEn: nitrogen agronomy efficiency. ¶IEn: nitrogen internal efficiency. #PFPn: nitrogen partial productivity. ¥Numbers at the same column followed by the same letter are nonsignificant at P<0.05. ¤Do not apply.Considering the values obtained in 2014 and 2015, REn decreased from 98.03 to 35.61, AEn from 36.32 to 16.78 kg kg−1, IEn from 63.64 to 33.94 kg kg−1, and PFPn from 67.96 to 39.56 kg kg−1 and, in most of the cases, the value decreased due increasing N rate. Ciampitti and Vyn [60] reported that the mean values of REn, AEn, IEn, and PFPn for new era hybrid corn were 44%, 22.9, 55.0, and 66.0 kg kg−1, respectively. In general, all NUE parameters from 2014 were similar or comparable to the mean values of the new era corn hybrid. However, in 2015, REn values were high especially in the pan evaporation treatments and the IEn values were low (Table 8).High REn values in 2015 may be attributed to increased aboveground N uptake; however, control treatments had lower N uptake values compared to fertilized treatments (Table 8). Therefore, N uptake differences between control and fertilized treatments were higher in 2015; consequently, this increased the REn values for 2015. Lower IEn values in 2015 can be explained by the high aboveground N uptake values.The NUEs equations are simple; however, they depend on several parameters, which may or may not be controlled by producers, for example, N rate or weather conditions, respectively. Therefore, there is potential to improve NUE and minimize negative environmental impacts if optimum practices for soil management, agronomy, ecology, and genetics are adopted with sustainable agriculture [52, 61].


## 3.1. Climatic Conditions and Irrigation Performance

Monthly average weather variables for the 2013, 2014, and 2015 growing seasons including the long-term average values (1971–2000) are summarized in Table3. Overall, climatic conditions helped explain differences in irrigation amount and grain yield. Differences in rainfall and temperature influence crop evapotranspiration in a given environment. Runge [44] found that daily temperature and rainfall affect corn yield from 25 days before to 15 days after anthesis. The seasonal average temperature was 22.60°C and 23.50°C in 2014 and 2015, respectively. The seasonal total rainfall in 2014 and 2015, respectively, was 502 mm (2.1% below normal) and 495 mm (3.7% below normal). The total rainfall was relatively similar across the two years and slightly below the (30 years) normal, which was 513 mm. However, rainfall distribution through the growing season will have greater effect on irrigation amount than total rainfall. Thus, if rainfall is concentrated in the period of high-water demand, irrigation requirements will be reduced and if rainfall is concentrated at the beginning and the end of the season, it will also limit irrigation requirements. In this study, rainfall was distributed differently before and after tasseling in 2014 but evenly in 2015. According to Lee (2011), tasseling is the highest corn water demand period. In 2014 and 2015, 41% and 58% of total rainfall and 49% and 51% of the total rainfall were evenly distributed before and after tasseling, respectively.Table 3
Monthly average climatic conditions, during the 2013, 2014, and 2015 growing seasons and long-term average values (30 years) measured at the research site TVREC.
MonthRainfall (mm)Temperature max (°C)Temperature min (°C)Average temperature (°C)2013April482210.315.9May5925.814.516June3031.519.825.2July10630.620.625.2August1530.720.124.6Total/average25928.117.121.42014April14723.31016.9May6127.615.221.2June1513120.325July933119.324.9August5032.719.825.5Total/average50229.116.922.62015April12923.311.917.9May11727.915.621.5June8832.320.626.1July10533.12227.2August5531.219.824.9Total/average49429.617.923.5Historic average (1971–2000)April10822.28.315.3May11126.713.620.1June10930.617.924.3July10632.42026.2August7932.118.925.5Total/average51328.815.722.3


## 3.2. Irrigation Amounts

Average irrigation amounts applied in the sensor plots were 106 mm and 142 mm in 2014 and 2015, respectively. However, the temporal distribution of rainfall between the two growing seasons was not comparable. Rainfall in the 2014 growing season was concentrated after tasseling, 294 mm, and in 2015, the rainfall distribution was similar before and after tasseling, 246 mm and 248 mm, respectively. Thus, considering that after tasseling corresponds to a period of greater water demand in corn, consequently, the 2015 season received much irrigation.Irrigation amounts for the two irrigation scheduling methods indicated that the sensor-based method received less irrigation during both growing seasons. The 2014 growing season received 138 mm of irrigation water on the plots under the pan evaporation irrigation scheduling method and the sensor-based irrigation scheduling method received 122 mm in the plots fertilized with 202 kg N ha−1 and 99 mm in the plots fertilized with 269 and 336 kg N ha−1. In 2015, the pan evaporation plots received 215 mm of irrigation and the sensor-based treatments received 152 mm at 202 and 336 kg N ha−1 plots and 127 mm at 269 kg N ha−1 plots (Table 4).Table 4
Seasonal irrigation on corn under different irrigation scheduling methods and N rates in 2014 and 2015 at TVREC.
`Irrigation amount (mm)20142015Pan evaporation0138215202138215269138215336138215Sensor-based0BNISǂBNIS2021221522699912733699152ǂBased on neighboring plot irrigation scheduling (BNIS).Figures2 and 3 illustrate soil water dynamics in the pan evaporation and sensor-based plots during the 2015 growing season, respectively. Daily soil matric potential data in the figures was collected from sensors installed in the plots fertilized with 269 kg N ha−1. The figures show how irrigation and rainfall affected soil matric potential values. From the pan evaporation method, it was possible to observe that even though a specific irrigation threshold was not set, the soil matric potential was maintained below 76 kPa, at approximately 50 kPa, due to more frequent irrigation events (Figure 2). In the sensor-based approach, the irrigation events were applied when the soil matric potential approached the threshold (76 kPa) (Figure 3). Therefore, since soil matric potential was maintained at lower values in the pan evaporation method, it infers that the soil had higher moisture during the growing season. Also, the sensor-based method depleted more water from the soil for irrigation events scheduled at a higher soil matric potential. After July 22, malfunctions occurred mainly in the sensors located in the 30 and 60 cm depth in this treatment (sensor-based irrigation fertilized with 269 kg N ha−1). However, irrigation practices in the sensor-based irrigation were continued using the shallowest sensor.Figure 2
Daily soil matric potential in the pan evaporation plot fertilized with 269 kg N ha−1 in 2015 at TVREC. The amount of irrigation referred to in the black arrows was value between 15 and 25 mm to avoid runoff.Figure 3
Daily soil matric potential in the sensor-based plot fertilized with 269 kg N ha−1 in 2015 at TVREC. The amount of irrigation referred to in the black arrows was value between 15 and 25 mm to avoid runoff.


## 3.3. Irrigation Scheduling Methods Effects: Pan Evaporation and Sensor-Based Irrigation

Differences between pan evaporation and sensor-based method were not significant with respect to all parameters analyzed in this study, except for REn (Table 5). Therefore, different amounts of irrigation applied using the two irrigation scheduling methods did not result in significant interactions among treatments that include irrigation (I x N, I x Y, and I x N x Y) for examined variables during 2014 and 2015.Table 5
Summary of ANOVA for grain yield, aboveground biomass, N uptake, and NUE with respect to N rates (0, 202, 269, and 336 kg ha-1), irrigation, and years (2014–2015) at TVREC.
Sources of variationGrain yield (Mg ha−1)Aboveground biomass (Mg ha−1)Grain N concentration (g kg−1)Stover N concentration (g kg−1)Grain N uptake (g kg−1)Aboveground N uptake (kg N ha−1)Nitrogen use efficiencyAEnǂREn§IEn¶PFPn#Irrigation (I)Ns£NsNsNsNsNsNsP<0.05NsNsN rate (N)P<0.05∗NsNsNsNsP<0.05P<0.05NsNsP<0.05Year (Y)P<0.05NsP<0.05P<0.05NsP<0.05P<0.05P<0.05P<0.05P<0.05I x NNsNsNsNsNsNsNsNsNsNsY x NP<0.05NsNsNsNsP<0.05P<0.05NsNsP<0.05Y x INsNsNsNsNsNsNsP<0.05NsNsY x I x NNsNsNsNsNsNsNsNsNsNs∗Significant at 0.05 probability level. £Nonsignificant at 0.05 probability level. ǂAEn: nitrogen agronomy efficiency. §REn: recovery efficiency. ¶IEn: nitrogen internal efficiency. #PFPn: nitrogen partial productivity.Grain yield is important indicator of profitability for farmers and there were not differences between irrigation methods when these variables were analyzed (Table5). During both growing seasons, pan evaporation scheduling method on average used 29% more water than sensor-based irrigation method. In 2014, the average reduction of water applied was 30 mm and, in 2015, the average reduction of water applied was 71 mm without statistically significant differences in grain yield or aboveground biomass within each year (Table 4). Based on the 2014 USDA Crop Production Summary (http://www.usda.gov), the nonirrigated corn yield average in the state of Alabama was 9.97 Mg ha−1 and the yield average obtained from the fertilized treatments of both irrigation methods was 13.84 Mg ha−1 and 13.18 Mg ha−1 during the 2014 and 2015 growing seasons, respectively. Glass et al. [45] conducted an evaluation of corn hybrids also at the Tennessee Valley research station and the irrigated corn yield average was 13.99 Mg ha−1.


## 3.4. Effects of N Rate on Corn Grain Yield and Aboveground Biomass

Significant grain yield differences were observed with respect to N rate, year, and interaction of year × N rate (Table5). In contrast, there was no aboveground biomass differences observed with respect to the same factors (Table 5). The lowest grain yield was observed on the control treatment during both years (Table 6). In 2014, the greatest grain yield was 14.23 Mg ha−1 achieved in the plot fertilized with 269 kg ha−1 and, in 2015, the greatest yield was 13.81 Mg ha−1 achieved in the plot fertilized with 336 kg N ha−1 (Table 6). In 2014, grain yield differences between the control and fertilized treatments were statistically significant; however no statistical differences were found between the fertilized plots (Table 6). In other words, the increase of N rate from 202 to 336 kg N ha−1 did not statistically change grain yield. However, the growing season of 2015 did not follow the same trend observed in 2014 and differences between control and fertilized plots and among fertilized plots were presented.Table 6
Grain yield and aboveground biomass differences between years and N rates (0, 202, 269, and 336 kg ha−1) at the TVREC.
N rate (kg ha)Grain yield (Mg ha−1)Aboveground biomass (Mg ha−1)201420152015MeanSEǂMeanSEMeanSE07.62b§0.525.52c0.3612.56b0.7920213.53a0.5212.54b0.3615.79a0.7926914.23a0.5213.22ba0.3614.65ba0.7933613.80a0.5213.81a0.3617.18a0.79ǂStandard error. §Numbers at the same column followed by the same letter are nonsignificant at P<0.05.The aboveground biomass was not significantly affected in 2014; however in 2015, the treatments receiving no N fertilizer presented the lowest value (Table6). In 2015, the aboveground biomass showed to be statistically responsive to the increasing N rates, which is commonly found in the literature [46, 47].Hagedorn et al. [48] mentioned that, in years with heavy rains, the mineral N content in the topsoil decrease by 50–70%; however, in low rainfall seasons, the pronounced mineralization occurs but the mineral N decreases only slightly. Therefore, the low rainfall in 2013 may have contributed to the no statistical differences in the fertilized plots in 2014. The 2013 season received 50% less precipitation than the normal of 30 years normal (Table 3). Furthermore, residual mineral N may have increase because leaching was limited due to low rainfall in 2013 and it could explain the low corn yield responses to N application observed in 2014. Extra N in 2014 could also explain why control plots in 2014 obtained higher yields than 2015.The impact of N rate on grain yield is important because of its effect on production profitability. Therefore, if the N fertilizer investment does not translate into a yield increase, then excess N could cause environmental problems. Some studies in this area suggest that increasing N application above optimal rates will lead to only small increases or even decreases in the grain yield and aboveground biomass. However, it may cause an increase in the energy use and production costs as well as increasing the risk of negative environmental effects such as greenhouse gas emission, soil acidification, and nitrate leaching [49–51].


## 3.5. Variation in Corn N Concentration

Grain and stover N concentrations were different between the two growing seasons with the N concentration values in 2015 greater than the values observed in 2014 (Table7). Overall, grain and stover N concentration increased as the N rate increased (Table 7). During the two growing seasons, grain N concentration ranged from 8.75 g kg−1 to 14.00 g kg −1 and the stover N concentration ranged from 0.61 g kg−1 to 13.23 g kg−1. Overall, grain and stover N concentration increased as the N rate increased (Table 7). Ciampitti and Vyn [52] reported the N concentration range for the new era corn cultivars was 3.0 g kg−1 to 26.8 g kg and 1.2 g kg−1 to 21.1 g kg−1 for grain and stover, respectively. Setiyono et al. [53] also reported large variation of N concentration values in corn plants across several studies; therefore this may explain the oscillation of N concentration in 2014 and 2015.Table 7
Grain N, stover N concentrations and grain, and aboveground N uptake differences between years and N rates (0, 202, 269, and 336 kg ha-1) at TVREC.
N rate (kg ha−1)Grain N concentration (g kg−1)Stover N concentration (g kg−1)Aboveground N uptake (kg N ha−1)201420152014201520142015MeanSEǂMeanSEMeanSEMeanSEMeanSEMeanSE08.75b§0.6611.70c0.210.61a1.688.11b0.87125.43b22.03166.68c39.0320210.53a0.6612.52b0.213.05a1.6810.50ba0.87216.55a22.03321.68b39.0326911.10a0.6613.51a0.211.15a1.6812.16a0.87230.51a22.03351.36ba39.0333611.54a0.6614.00a0.212.45a1.6813.23a0.87241.79a22.03421.64a39.03ǂStandard error. §Numbers at the same column followed by the same letter are nonsignificant at P<0.05.In 2014, yields were higher which indicates that plants were able to use the applied N. Therefore, the N dilution could occur because higher yields and more biomass may reduce an equivalent N concentration in the plant in 2014. Moreover, temperature could be another reason of the increase of N concentration in the corn plants. Wienhold et al. [54] compared N concentration in corn and obtained 75% higher N content values in years when the temperature was above the 30-year average. Therefore, higher temperature in 2015 (Table 3) could explain the increase of N concentration for this growing season.In 2015, grain N concentration differences between the control and fertilized treatments were smaller than the differences observed in the 2014 growing season. Therefore, in 2015 the corn plants uptake more N from the unfertilized plots, showing the ability of corn to uptake nitrogen from the soil N pools even when no N fertilizer is applied (Table7). In 2014, N rate did not affect the N concentration in the stover, but in 2015, stover N concentration increased following increment of nitrogen fertilizer.Aboveground N uptake ranged from 125.43 kg N ha−1 to 421.64 kg N ha−1 during the two growing seasons. Ciampitti and Vyn [49] reported the mean, minimum, and maximum N uptake of new era maize cultivars and the values were 184, 2, and 427 kg N ha−1, respectively. Other numbers reported by Setiyono et al. [53] showed 239, 7, and 471 kg N ha−1 as mean, minimum, and maximum, respectively. Therefore, the maximum and minimum N uptake values obtained in Alabama’s environment in this present study agreed with the range of new era cultivars reported by other authors. Aboveground N uptake was higher in 2015 when compared to 2014. In 2014, the aboveground N uptake between control and fertilized treatments was statistically significant; however considering only the fertilized treatments, the N uptake values increased numerically but these differences were not significant (Table 7). During the season of 2015, the aboveground N uptake increased significantly with the increase of N rate (Table 7). Abbasi et al. [55] reported that aboveground N uptake increased with N rate; therefore in our study, the N uptake increased as the N rate increased but the increment was not significant in all cases.Setiyono et al. [53] and Chua et al. [56] reported that aboveground N uptake has been closely linked with grain and aboveground biomass. However, we observed that, in general, grain yield was higher in 2014 but the overall aboveground N uptake was lower when comparing with 2015. Herrmann and Taube [57] found lower N uptake in higher corn biomass yield and it might be due the N dilution because the N is spread out in the plant tissue.Furthermore, other studies evidenced that aboveground N uptake is not proportionately linked with grain yield. Jokela and Randall [58] presented some data from a three-year study (1982–1984) comparing grain yield and aboveground N uptake and found years with higher yield and lower N uptake or vice versa. Uribelarrea et al. [59] compared four commercial corn hybrids in three years and also found higher yield and lower N uptake; furthermore, some hybrids had a negative N balance removing more N than was applied as fertilizer-N. It also occurred in our study in the 2015 growing season (Table 7).


## 3.6. Nitrogen Use Efficiency

The Nitrogen Use Efficiency (NUE) increase indicates that most N fertilizer applied in the soil was taken up by plants and converted into grain yield. Therefore, the amount of N that remains in the soil is lower decreasing the potential of leaching and can effectively decrease the potential negative impacts of N fertilizer to the environment [55, 60].The AEnand PFPn followed the same trend in 2014 and 2015; higher AEn and PFPn were observed for the treatment fertilized with 202 kg N ha−1 and gradually decreased as N rate increased (Table 8). Also, the lowest NUE value was observed from the plots fertilized with the highest N rate (336 kg N ha-1). According to Ciampitti and Vyn [60], it is well accepted the NUEs values are high at low N rates and decrease with increasing N rates. Furthermore, AEn and PFPn may be more practical parameters because grain yield and N rate are used in the equation; therefore they are easier parameters to be calculated by farmers. The REn and IEn did not present significant differences during both seasons. Also considering that these parameters use aboveground N uptake in the equation, it might be more laborious for the farmers because it is necessary to do the N concentration analysis.Table 8
Nitrogen recovery efficiency (REn), nitrogen agronomic efficiency (AEn), nitrogen internal efficiency (IEn), and nitrogen partial productivity (PFPn) at four N rates (0, 202, 269, and 280 kg N ha−1) in 2014 and 2015 at TVREC.
N rate (kg ha−1)REn (%)ǂAEn (kg kg−1)§IEn (kg kg−1)¶PFPn (kg kg−1)#20142015201420152014201520142015MeanSEǂMeanSEMeanSEMeanSEMeanSEMeanSEMeanSEMeanSE0¤¤¤¤¤¤¤¤62.57a2.3835.04 a2.16¤¤¤¤20236.94a¥4.4498.03a16.7924.96a2.7936.32a0.6763.64a2.3840.38a2.1667.96a2.5962.33a1.8226943.96a4.4480.27a16.7922.7a2.7928.79b0.6762.76a2.3838.16a2.1652.32b2.5949.32b1.8233635.61a4.4485.58a16.7916.78b2.7923.62c0.6758.37a2.3833.94a2.1639.56c2.5941.22c1.82ǂRen: recovery efficiency. §AEn: nitrogen agronomy efficiency. ¶IEn: nitrogen internal efficiency. #PFPn: nitrogen partial productivity. ¥Numbers at the same column followed by the same letter are nonsignificant at P<0.05. ¤Do not apply.Considering the values obtained in 2014 and 2015, REn decreased from 98.03 to 35.61, AEn from 36.32 to 16.78 kg kg−1, IEn from 63.64 to 33.94 kg kg−1, and PFPn from 67.96 to 39.56 kg kg−1 and, in most of the cases, the value decreased due increasing N rate. Ciampitti and Vyn [60] reported that the mean values of REn, AEn, IEn, and PFPn for new era hybrid corn were 44%, 22.9, 55.0, and 66.0 kg kg−1, respectively. In general, all NUE parameters from 2014 were similar or comparable to the mean values of the new era corn hybrid. However, in 2015, REn values were high especially in the pan evaporation treatments and the IEn values were low (Table 8).High REn values in 2015 may be attributed to increased aboveground N uptake; however, control treatments had lower N uptake values compared to fertilized treatments (Table 8). Therefore, N uptake differences between control and fertilized treatments were higher in 2015; consequently, this increased the REn values for 2015. Lower IEn values in 2015 can be explained by the high aboveground N uptake values.The NUEs equations are simple; however, they depend on several parameters, which may or may not be controlled by producers, for example, N rate or weather conditions, respectively. Therefore, there is potential to improve NUE and minimize negative environmental impacts if optimum practices for soil management, agronomy, ecology, and genetics are adopted with sustainable agriculture [52, 61].


## 4. Conclusions

The sensor-based Irrigation scheduling method tends to have a better performance during the 2014 and 2015 relative to the pan evaporation method. Even though less water was applied using the sensor-based irrigation scheduling method, the differences between the methods might be affected by the lack of irrigation that occurred on both seasons due to the failure of the irrigation pump during the period of high crop water demand. However, the results showed that the sensor-based method was a reliable method to monitor soil water status and provide useful data for producers to manage irrigation scheduling. The pan evaporation method might be less accurate in terms of soil moisture content because it assumes uniform soil conditions and is based on estimations of the soil water balance. When comparing both irrigation scheduling methods tested in this experiment, 29% more water was applied on average using the pan evaporation irrigation scheduling method compared to the sensor-based method. Despite these differences, there were no significant grain yields or aboveground biomass differences observed between the two irrigation scheduling methods. Moreover, the lack of water application caused by the problems in the irrigation pump may have affected the yield potential of grain and aboveground biomass; thus no significant differences were presented. Nitrogen rates evaluated in this study affected grain yield, aboveground N uptake, and NUE expressed as AEn and PFPn. In general, NUEs values decreased with an increase in N rates, which means that extra N fertilizer added to the soil was not converted into grain yield or/and uptake by plants, which might increase the risk of future environmental problems.



---

*Source: 8869383-2020-09-15.xml*