1. Introduction
Nitrogen (N) is an essential macronutrient for all organisms and N-containing fertilizers are the most widely applied fertilizers in agriculture [
1]. Among N-based fertilizers, urea is the most commonly used [
2] due to its high N content (46%), its relatively low cost, and its ease of use [
3]. However, the efficiency of urea application is low due to its high hydrolysis rate, which generates more mineral N than the crop can assimilate in the first stage. Some of this N can be immobilized in the soil, but the excess is lost through ammonia (NH
3) volatilization or nitrate (NO
3−) leaching [
4]. In order to avoid NH
4+ losses in soil, there are several techniques that have proven their effectiveness. It is possible to hinder the nitrification process (either through anaerobic conditions or through the use of nitrification inhibitors) to avoid excess NO
3− in the soil and its potential leaching. To avoid losses due to NH
3 volatilization, the tactic of incorporating NH
4+ into the soil is generally used, either by burying the fertilizer with mechanical labor, by introducing it with irrigation or taking advantage of rain, or by injecting the fertilizer directly into the soil. In this way, the release into the atmosphere of spontaneously generated gaseous NH
3 is inhibited and its conversion to NH
4+ is improved. Finally, the fixation of NH
4+ in the soil can be enhanced by increasing the ion exchange capacity of the soil, which is related to the clay proportions in the soil, or by increasing the amount of organic matter [
5]. However, these processes, except when in artificial soils, are slow and expensive.
The global average losses of NH
3 from urea fertilizers are estimated to be close to 14% (range of 10–19%) [
6], but they can reach up to 40% of applied urea-N in warm and tropical regions [
7]. This not only implies economic losses for farmers, but also generates a negative impact on the environment and human health due to the eutrophication generated by the leaching of nitrates, the greenhouse gas emissions, and the increase in air fine particles caused by ammonia volatilization [
8].
One way to reduce N losses consists of adjusting the N concentration in soil to the crop demand, which can be achieved by dividing the N application into several dressings during the crop development [
9]. As some authors have pointed out, split applications of N can improve crop yields and N use efficiency parameters, such as the agronomic efficiency (AE
N) and the N recovery efficiency (RE
N) [
10]. However, split applications have some downsides, such as an increase in the amount and difficulty of field work, which typically leads to an increase in the overall crop costs. Another proposed method for reducing N loses is the use of urease inhibitors. Urea hydrolysis is facilitated by urease, a nickel-dependent enzyme that catalyzes the hydrolysis of urea to one mole of carbon dioxide (CO
2) and two moles of ammonia, which can oxidize to nitrate. As a key enzyme for the global N cycle, this hydrolase is widely distributed in nature and is released naturally by soil microorganisms [
11]. Urease inhibitors reduce urea hydrolysis, adjusting the release of N mineral forms to match crop demand, therefore increasing the crop nitrogen use efficiency (NUE). However, their effect on crop yield is variable [
12]. The reduction in NH
3 losses caused by the urease inhibitors can range from 0% to 94% of the NH
3 lost when using urea alone. Crop yield can also vary depending on the crop management and conditions, but overall, Silva et al. [
13] calculated in a meta-analysis that urease inhibitors can increase crop yield by up to a 5%. In the last decade, climate change has become a growing threat to crop development, with more frequent harsh climate conditions leading to a crop yield reduction. In this context, water shortages and droughts will become an increasingly common and severe problem in warm countries. Therefore, it is of great interest to study methods to reduce the impact of water shortage on crops, making the effects derived from the interaction between N nutrition and water availability especially relevant [
14].
Globally, maize and wheat are the second and third most produced crops, respectively. In 2019, 1148 Mt of maize and 766 Mt of wheat were produced, which respectively represented 38.5% and 25.7% of the world’s total cereal production [
15]. In the 2020–2029 period, global cereal production is projected to expand by 375 Mt, reaching 3054 Mt in 2029, mainly driven by higher yields. In Spain, wheat and maize are two of the most important crops and, therefore, increasing the economic and environmental benefits in these systems is key to increasing their medium- and long-term sustainability.
This experiment was based on the concept that the better the coupling between the mineral nitrogen in the soil and the crop demand, the better the nitrogen use efficiency. We expected to prove that the application of urease inhibitor could produce an increase in the nitrogen use efficiency, reducing, at the same time, machinery use. Therefore, the main goal of this experiment was to quantify the effect the application of a urease inhibitor, joined to the traditional urea fertilizer with two different fertilizer application schedules and irrigation rates, had on nitrogen use efficiency. The impact of these factors was quantified based on the yield response, the grain quality, and the soil mineral nitrogen and enzymatic activity.
2. Materials and Methods
2.1. Field Experiment
This field experiment was conducted from April 2018 to July 2019 at the “La Canaleja” field station of the National Institute of Agricultural and Food Research and Technology (INIA). The field station is located in Alcalá de Henares, (Madrid province), Spain (40°32′ N, 3°20′ W, 600 masl). The climate is Mediterranean semi-arid, with an annual rainfall of 389 mm (186–547 mm), mainly occurring in autumn and spring and almost negligible in summer, and a mean annual temperature of 13.5 °C. The soil was classified as a
Typic Calcixerept by the soil taxonomy classification [
16] and as a
Calcic Cambisol by the Food and Agriculture Organization (FAO) classification [
17]. Further information can be found in
Table 1.
The experiment included a maize and wheat crop rotation, which consisted of maize sown on 14 May and harvested on 8 October, followed by wheat sown on 14 February and harvested on 10 July. The maize (plant density of 80,000 plants ha−1) and wheat (sowing rate of 230 kg ha−1) were grown in an irrigated field with a 16 × 16 m2 permanent sprinkle system. The P and K fertilization consisted of 50 kg P2O5 and 50 kg K2O, before maize sowing, following the traditional dose in the region, in order to avoid a phosphorus (P) or potassium (K) deficit. The soil management consisted of conventional tillage plowing to 20 cm depth two weeks before maize or wheat sowing, followed by disc plowing for sowing bed preparation. All crop residues were left over the surface until plowing.
The maize fertilizer rate was 170 kg N ha
−1, based on the optimum observed for the region in previous experiments [
18,
19]. The five fertilization treatments consisted of: C (control without N fertilization), U
2d (170 kg N ha
−1 of urea split into two dressings at 4–6 and 8 fully expanded leaves), U
1d (170 kg N ha
−1 of urea applied at 4–6 fully expanded leaves), IN
2d (170 kg N ha
−1 of urea with urease inhibitor split into two dressings at 4–6 and 8 fully expanded leaves), and IN
1d (170 kg N ha
−1 of urea with urease inhibitor applied at 4–6 fully expanded leaves). The first fertilizer dressing (85 or 170 kg N ha
−1 depending on the treatment) was applied on 18 June and the second on 17 July. The five treatments were duplicated based on two different irrigation rates. Each treatment had four replications randomly distributed in forty 8 × 8 m
2 plots included in a 1.25 ha maize field. The spring wheat crop was sown over the previous maize treatments and received no fertilization for a better appreciation of the residual effect of urease inhibitor, application schedule, and irrigation intensity on yield and grain quality.
The irrigation schedule and doses were estimated from the daily values of crop evapotranspiration (ETc) and rainfall. This was calculated as ETc = Kc × ET
0, where ET
0 was the reference evapotranspiration, calculated by the FAO Penman-Monteith model [
20] using daily local data, and Kc was daily obtained by the Martinez-Cob [
21] method. Total water input was divided into plots with an optimal irrigation intensity equivalent to the ETc demand (100%) and plots with suboptimal irrigation (75%) in order to analyze the effect of water deficit. The total water irrigated during the maize crop was 700 mm (in the 100% treatments) and 525 mm (in the 75% treatments). As the applied water presented an average of 12.3 mg N-NO
3− L
−1 over the course of the cropping season, the total amount of nitrate applied via irrigation was 86.0 and 64.5 kg N ha
−1 in the 100% and 75% treatments, respectively. The spring wheat was also irrigated in order to supply the difference between the potential evapotranspiration and the actual precipitation. The total water irrigated was 145 mm (in the 100% treatments) and 110 mm (in the 75% treatments). In this case, the N incorporated with the irrigation water was 17.5 and 13.1 kg N ha
−1 in the 100% and 75% treatments, respectively.
2.2. Crop Analysis: Yield, Grain Quality, N Content, and Nutritional Status
At maize harvest, a 5 m stripe in the central row was harvested by hand and separated into plant components (grain vs. rest of aboveground biomass), and a subsample of each component was oven-dried (65 °C) and weighed. A subsample of each plant component was also used to determine total N concentration by the Dumas combustion method (LECO FP–428 analyzer, St. Joseph, MI, USA). For each plot, the N content of each crop component was calculated by multiplying its dry biomass by its N concentration and adding both to obtain the total crop N content. At wheat harvest, the grain biomass was registered by sampling a 5 m stripe in a central row by hand. A subsample was oven-dried (65 °C), weighed, and analyzed for N concentration and total N content.
The crop nutritional status of the maize was evaluated with a Dualex
® Scientific (Force-A, Orsay, France) chlorophyll meter, a leaf clip sensor that measures chlorophyll content (Chl) as the difference between the light transmitted at the red and infrared wavelengths. The device also measured flavonols (Flav) and anthocyanins (Anth), two substances which are highly correlated with the crop stress and the nutritional balance index (NBI) [
22,
23]. Readings were carried out over the maize crop cycle from 12 June (before first dressing) to 3 October (before harvest). On each sampling date, 20 measurements were taken from the uppermost fully developed leaf (or in the ear leaf once flowering) of 20 representative plants in the central rows of each plot, following the sampling method defined for maize by Gabriel et al. [
24]. The representative value of each plot was the twenty-measurement average.
The crop nutritional status of the wheat was evaluated with RapidScan® (Holland Scientific Inc., Lincoln, NE, USA). In this case, the measurements consisted of the proximal measurement of the normalized difference vegetation index (NDVI, structural and growth index) and the normalized difference red edge (NDRE, nitrogen content index), integrating measurements every 0.5 s in the central 7 m row per plot. These measurements were taken three times: 9 April, 3 May (flowering), and 30 May.
2.3. Soil Inorganic N Content (Nmin) and Urease Activity
The soil nitrate and ammonium content were determined at four dates: before the maize crop sowing (11 May), at flowering (6 August), after maize harvest (11 October), and after wheat harvest (14 July). A combined sample from three soil cores was taken from each plot by 0.2 m intervals with an Eijkelkamp
® helicoidal auger (Eijkelkamp Agrisearch Equipment, Geisbeek, the Netherlands). The samples were taken at a depth of up to 0.6 m. The samples were placed in a plastic box and firmly closed immediately, transported, and refrigerated (4–6 °C). Within the five consecutive days, a soil subsample of each box was extracted with 2M KCl (~5 g of soil: 50 mL of KCl), centrifuged, decanted, and a subsample of the supernatant volume was stored in a freezer until later analysis (FIAstar™ 5000, FOSS Analytical AB, Höganäs, Sweden). The nitrate concentration was determined by the Griess-Ilosvay method [
25], and ammonium was measured using the method of Solorzano [
26]. Soil N
min was calculated for each layer and plot. The soil urease activity was also measured in the soil samples obtained at maize flowering and harvest, quantifying the residual effect of the inhibitor in the soil. The methodology followed was that suggested by Kandeler and Gerber [
27].
2.4. N Efficiency Parameters
The components of N use efficiency that were calculated for each crop included: agronomic efficiency (AE
N), the N recovery efficiency (RE
N), nitrogen use efficiency (NUE), and nitrogen surplus (N
surplus) based on Ayuso, Gabriel and Quemada [
28]. The AE
N (Equation (1)) refers to the kg of crop yield increase obtained per kg of N applied (i.e., the ratio of the difference between the grain yield of a treatment (Grain_yield
Treatment) and the average grain yield of the control (Grain_yield
Control) to the N applied in the specific treatment from fertilizer (N
Fertilizer) and irrigation water (N
irrigation)). The RE
N (Equation (2)) refers to the kg of crop N uptake per kg of N applied, calculated as the ratio of the difference between the crop N uptake by a treatment (N
uptake-Treatment) and the average crop N uptake by the control (N
uptake-Control), to the N fertilizer applied. The NUE (Equation (3)) was calculated as the N exported from the system as grain (N
uptake-grain) divided by the N input (N from fertilizer and irrigation). Finally, the N
surplus (Equation (4)) was calculated as the difference between the N inputs from the fertilizer, the irrigation water, and initial mineral N (N
min_ini) and the N outputs (as crop N uptake (N
uptake) and final soil mineral N (N
min_final).
2.5. Statistical Analysis
Statistical analysis was carried out using SPSS for Windows (v. 23.0). A Shapiro-Wilk test was carried out in order to assess data homogeneity. All multifactor analyses were performed using a GLM or Kruskal-Wallis test. Post hoc multiple comparisons of means were carried out using Duncan’s or Games-Howell´s test as appropriate. Only a p value of < 0.10 or lower was considered statistically significant.
3. Results
3.1. Maize Crop Analysis
Over the maize crop cycle, there were differences in maize nutritional status between fertilized and non-fertilized treatments and between optimally (100%) and sub-optimally (75%) irrigated plots, but differences were not so clear between the different fertilizer treatments (
Figure 1). Nevertheless, some tendencies require further analysis. For instance, both treatments with split dressings presented some incipient deficiencies (lower NBI and Chl indexes and higher Flav and Anth) during the period between the first and the second dressing. However, while the U
2d treatment recovered after the second dressing, the IN
2d needed more time. This effect was more evident in the NBI index. Moreover, the Anth were more related to water stress and the Flav to nitrogen.
At maize harvest, there were differences in maize height, grain yield, and straw biomass between the control and the fertilized treatments (
Figure 2). Looking at the fertilized treatments, there were no final differences between both regular urea treatments (U
1d and U
2d) and IN
1d, but the three presented larger growth values than the IN
2d. These results (joined to the maize nitrogen status observed with sensors) indicated that the urease inhibitor delayed N availability for the crop. On the one hand, when the fertilizer with the inhibitor was applied at early crop stages, the N mineralization was coupled with the crop demand. However, when the inhibitor application was split, the N availability was decoupled.
Irrigation intensity was the most significant factor (
Table 2) in terms of yield, aerial biomass, and plant height, with a consistent reduction in the values obtained for each of these variables under sub-optimal irrigation conditions (75%). However, the detrimental effect of sub-optimal irrigation was attenuated by using urease inhibitors, although this relationship was only significantly established for plant height (
Figure 3).
The number of fertilizer dressings also had an important impact on the data obtained. Regardless of fertilization type, one dressing treatments showed the largest grain yield, aerial biomass, plant height, N concentration, and N uptake in grain, although these results were only significant compared to the non-fertilized control. However, on the other hand, N concentration and uptake in the rest of the aerial biomass was higher in the two dressings treatments, although these differences were not significant (
Table 3).
When the use of urease inhibitors was studied alone, data showed low and a statistically non-significant impact of this factor on the data obtained. However, significant differences were observed in the joint effect of inhibitors with other factors. As seen previously, the detrimental effect caused by the reduction in irrigation intensity to sub-optimal levels (75%) can, in turn, significantly reduce the yield, biomass, and N concentration of the crop. However, this effect was slightly attenuated in our experiment after the use of urease inhibitors, as this effect was statistically significant for plant height. A similar case was observed regarding the combination of urease inhibitors with the different number of dressings. In this case, a light and non-significant increase was observed in the N uptake in grain when all the fertilizer was applied in one dressing. On the other hand, when the fertilizer was divided in two dressings, the N uptake in grain was reduced with the use of urease inhibitors (
Table 3).
It is interesting to analyze the result obtained between the interaction of the three factors (inhibitor, irrigation, and dressings) and the N concentration in grain. It was observed that when no inhibitor and only one dressing were applied, a non-significant increase in the concentration of N in the grain was obtained in the sub-optimal irrigation conditions (75%). This difference disappeared when the crop was fertilized with two dressings. However, the opposite effect was observed when urease inhibitors were applied, obtaining non-significantly higher N concentrations in the grain of the optimally irrigated maize (100%) when one dressing was applied, while, when two dressings were applied, the N concentration in grain was not significantly higher in corn with suboptimal irrigation (75%). The results, therefore, showed that urease inhibitors could increase the N concentration in grain when applied together with two dressings under sub-optimal irrigation conditions. In all cases, the absence of fertilization (C) showed the worst results (
Figure 3).
Grain N concentration in IN1d,100%, IN2d,75%, and all urease alone (U) treatments, except for U1d,75%, showed higher values than both controls. Nevertheless, while the differences were generally not significant, both IN1d,100% and IN2d,75% presented the highest N concentrations in grain among all the treatments studied. On the other hand, the optimally irrigated IN1d, U2d, and U1d showed a tendency to increase Nuptake in grain when compared to the other treatments, especially the controls. Regarding Nuptake and N concentration in aerial biomass, the differences between treatments were not statistically significant, and, therefore, it was not possible to obtain any remarkable data from them.
3.2. Wheat Crop Analysis
As the wheat crop was not fertilized, it showed a great dependency on the residual nitrogen content left in the soil by the previous crop, maize. In this case, the greater effect observed was the N availability in the previous maize sub-optimally irrigated treatments. These treatments presented larger wheat development, observed in a larger NDVI, and a better nitrogen status, observed in a larger NDRE (
Figure 4). This also corresponded with the N content and absorption observed at harvest in the grain. Sub-optimal irrigation (75%) was associated with an increase in N concentration (
p < 0.01) and N uptake (
p < 0.05) in grain. However, the irrigation effect was not so clear in the plant development. Urease inhibitor use and the number of dressings showed no statistically significant differences between treatments, although wheat treated with the urease inhibitor appeared to have higher grain yield than the urea-alone fertilized counterpart. The number of dressing seemed to have very little impact on the data obtained, although the two dressing treatments showed a slightly and non-statistically significant increase in the results obtained compared to the one dressing treatments (
Figure 5).
Among all the treatments studied, IN
2d seemed to have the best performance for all the studied factors, followed by U
2d, reinforcing the previous maize nitrogen residual effect (
Figure 5). Within the fertilized treatments, U
1d and IN
1d showed the lowest performance, being the no fertilized control, as the treatment with the lowest results, even though in all cases, these differences were not significant.
3.3. Soil Inorganic N Content
Initial soil inorganic N content in the 0–60 cm layer was slightly high (110 kg N ha
−1), even though rapeseed was sown in order to homogenize and extract the residual N
min excess. The entire field was homogeneous and no differences were found between treatments. At flowering, the N
min changed from one treatment to the other and both the number of dressings and the use of inhibitor induced an effect. On the one hand, when the entire N was applied in one dressing, the use of urease inhibitors in the IN
1d treatment allowed a larger amount of N
min in the upper layer at flowering, reducing the N
min at the lower, and reducing the leaching risk (
Figure 6). Comparatively, U
2d increased the N
min in the upper layer greatly, but a considerable amount was leached down along the profile, increasing the nitrate leaching risk. Finally, the IN
2d drastically reduced the amount of N
min available at flowering along the profile, which reduced the nitrate leaching risk, but also reduced the N availability for the maize to values similar to those of the non-fertilized control. In general, sub-optimal irrigation produced larger accumulation of N
min in the upper layer, mainly due to the reduction in drainage and a reduction in the plant uptake.
After maize harvest, the soil was depleted of Nmin. Both U treatments presented 62 kg N ha−1 at the upper 60 cm, between 34 and 44 kg N ha−1 in the IN treatments, and 28 kg N ha−1 in the control. The amount of NH4+ was similar in all treatments and depths (around 7 kg N ha−1 per treatment and depth), but NO3− changed depending on the treatment and at different depths. On the surface, all the fertilized treatments presented a similar amount of NO3− (12 kg N ha−1 on average), larger than the control treatment (4 kg N ha−1). Moreover, the fertilized treatments applied in two dressings tended to increase the residual NO3− with respect to that applied in one dressing at the surface (on average, 13.5 vs 10.3 kg N ha−1, respectively). At a depth of 30 cm, differences increased between the treatments. In general, treatments fertilized with regular urea presented larger residual NO3− than the IN treatments (with an average difference of 7.7 N ha −1), and treatments that applied the fertilizer in two dressings presented larger values than those seen in the one dressing application (with an average difference of 4.1 N ha−1). Differences with respect to the control treatments were reduced for the IN treatments, but not for the U treatments. Finally, in the deeper layer, this tendency persisted and the IN treatments and the control presented similar results (2.3 kg N ha−1 on average), much lower than the U treatments (15.7 kg N ha−1 on average). Without plants in the field, this NO3− in the lower part of the soil profile was prone to be leached via successive rainfall events, but it could also be indicative of previous leaching during the maize cropping season. It was also remarkable how inhibitors reduced the proportion of Nmin as NO3− with respect to the NH4+ (approaching the natural values observed in unfertilized soil), when compared with the regular urea fertilized treatment, reducing the risk of NO3− leaching. Finally, the same tendencies observed after maize harvest were observed after wheat harvest. In this case, the differences were even smaller and all soils reached similar values of depletion as those observed in the unfertilized control plots. There was no effect from irrigation, the use of inhibitors, or the fertilizer timing.
3.4. Soil Urease Activity
The urease activity data showed greater enzymatic activity for treatments with sub-optimal irrigation. Regarding the treatments with poor irrigation, the treatments without inhibitors (U
1d, U
2d and C) showed a similar behavior over time with high urease activities on the surface, which tended to decrease around 30 cm in depth, and a general increase at 50 cm depth (
Figure 7). Treatments with urease inhibitors, on the other hand, presented other behaviors, with the IN
2d treatment showing a constant urease activity increase while descending in the soil profile. The IN
1d treatment showed relatively low activity in the entire soil profile during flowering, while at harvest this activity increased at the 30 cm depth.
Regarding the optimal irrigated treatments, soil urease activity at flowering showed differences between treatments. On the one hand, IN treatments presented lower urease activity (as expected) in the upper layer (at 0–20 cm depth in both IN1d and IN2d treatments, and at 0–5 cm depth for the IN2d treatment only). This effect showed that the inhibition power was still working even after the 49 days between the first dressing and the flowering sampling. On the other hand, there was also an increase in the upper layer urease activity in the U treatments (larger for U1d) with respect to the unfertilized control. This effect suggested an increase in the microorganisms decomposing urea due to the larger urea availability. This tendency was not so clear in the deeper layers. In this case, the treatments where urea was applied in one dressing (U1d and IN1d) presented some tendency to increase urease activity, likely due to urea leaching down without the inhibitor. The variability was very high at these layers, and thus the conclusions were not so clear. However, the reduction in the irrigation presented the most important effect. All treatments, including the unfertilized control, presented higher values for the urease activity along the soil profile, indicating that under these drier conditions, soil microorganisms were more active.
Soil urease activity at harvest was very similar, and values moved on the same range as at flowering. Between 20 and 60 cm, the values were very similar to those obtained at flowering and with the same variability, making conclusions unclear. However, on the surface, some differences appeared. In this case, both the U treatments demonstrated larger urease activity than the control. Moreover, the IN2d treatment showed a similar tendency, with a fast recovery of soil microorganism populations. However, the IN1d treatment still presented lower urease activity than the unfertilized control.
3.5. N Efficiency Parameters
Regarding the N efficiency parameters, irrigation was the factor with the greatest effect on the AE
N, RE
N, and NUE parameters, being able to reduce up to 50% of the efficiency of the use of N in crops developed without urease inhibitors, while N
surplus could be multiplied by three (
Table 4). However, the results indicated that the detrimental effect on N efficiency observed in the sub-optimally irrigated treatments could be ameliorated significantly through the use of urease inhibitors, achieving reductions in AE
N between 6–27% depending on the treatment. In the case of the IN
2d treatment, the use of urease inhibitors increased RE
N by more than 100%. At the same time, N
surplus was reduced as a consequence of the use of urease inhibitors in the IN
2d treatment. However, this was not the case with the IN
1d treatment, where N
surplus increased with insufficient irrigation (75%), even though this increase was lower for the IN
1d treatment than for any of the U treatments.
In general, when the irrigation effect was not considered, the IN1d treatment presented with a high efficiency, which was almost equal to the efficiency obtained by the U2d treatment. The efficiency obtained by the U1d treatment was slightly lower than it was for U2d and IN1d. Finally, the IN2d treatment presented with the lowest of values, mostly because of its lower yield and total N uptake.
3.6. Economic Efficiency Parameters
There were three main differences in the economic inputs between treatments: the amount of fertilizer, the price of the fertilizer, and the number of applications. The amount of fertilizer was the same (370 kg of urea ha
−1) in the four fertilized treatments, while it was zero in the control. With respect to the price of the fertilizer, it was EUR 27.50 for the sack of 40 kg of regular urea and EUR 29.50 for the urea with inhibitor. Finally, each fertilizer application had a cost of EUR 5 ha
−1 (considering only fuel consumption and tractor/machinery amortization). The average maize grain price (not including the recent global increase) in Spain was EUR 0.174 kg
−1, standardized to 14% of humidity [
29]. Equally, the average wheat grain price was EUR 0.185 kg
−1, standardized to 11% of humidity [
29]. Considering all these values, the final benefit for each treatment was calculated. As a result, the IN
1d treatment had the most interesting relative benefit with EUR 2714.45 ha
−1 (3099.3 if the wheat period was included), but without statistical differences with U
2d or U
1d (
Table 5). The IN
2d treatment presented a lower economic benefit, mostly because of the combination of the high cost and the lower incomes, but without statistical differences with respect U
1d. The unfertilized control presented the lowest benefits, even when the relative cost was EUR 0 ha
−1, due to the very low yield. Moreover, the IN
1d presented larger benefits when irrigation water was a limitation, increasing the differences between U
1d and U
2d. Under these conditions, IN
2d also increased in benefits, reaching the margins obtained by the U treatments.
5. Conclusions
A single application of urease inhibitor (IN1d) coupled with the conventional urea can help to reduce the nitrate leaching risk both during the maize period (even when compared to the two dressing treatments) and after harvest. In addition, this improvement was achieved together with an increase in economic benefit, even when compared with the application of the same amount of regular urea split into two dressings. Moreover, the N efficiency parameters showed that it was not only the most efficient treatment studied, but also that the amount of fertilizer applied could be reduced in a larger proportion than in the others (as the larger Nsurplus suggests), and the number of dressings could be reduced from two to one. However, when the urea with the inhibitor was split in two dressings, the Nmin availability was uncoupled with the maize demand, reducing maize yield, and N use efficiency.
The inhibitor effect seemed to persist more than 100 days in the soil when the application rate was high, but it also seemed that the inhibitor presented a reduced movement capacity along the profile, reducing the risk of ground water contamination with these molecules. Finally, under low water availability systems, the benefits of applying urease inhibitors increased with respect to the application of regular urea, making this technique a very promising strategy for adaptation to climate change under arid and semiarid regions.
All the results observed in this study require further confirmation from additional experiments from different fields and laboratories to test not only the possible different behaviors, but also the different aspects that we were not able to demonstrate in this one.