New Eco-Friendly Polymeric-Coated Urea Fertilizers Enhanced Crop Yield in Wheat

Presently, there is a growing interest in developing new controlled-release fertilizers based on ecological raw materials. The present study aims to compare the efficacy of two new ureic-based controlled-release fertilizers formulated with water-soluble polymeric coatings enriched with humic acids or seaweed extracts. To this end, an experimental approach was designed under controlled greenhouse conditions by carrying out its subsequent field scaling. Different physiological parameters and crop yield were measured by comparing the new fertilizers with another non polymeric-coated fertilizer, ammonium nitrate, and an untreated ‘Control’. As a result, on the microscale the fertilizer enriched with humic acids favored a better global response in the photosynthetic parameters and nutritional status of wheat plants. A significant 1.2-fold increase in grain weight yield and grain number was obtained with the humic acid polymeric fertilizer versus that enriched with seaweed extracts; and also, in average, higher in respect to the uncoated one. At the field level, similar results were confirmed by lowering N doses by 20% when applying the humic acid polymeric-coated produce compared to ammonium nitrate. Our results showed that the new humic acid polymeric fertilizer facilitated crop management and reduced the environmental impact generated by N losses, which are usually produced by traditional fertilizers.


Introduction
According to the Food and Agriculture Organization of the United Nations (FAO), wheat is the world's largest cultivated crop per hectare and the third largest cereal to be produced [1]. In the European Union, 144.5 million tons were harvested in 2016 and production is estimated to increase by 3.5% each year. In fact, the world's production in 2017 was expected to come to 744.5 million tons, an increase that comes close to 1000% since 1990/1991. Current cereal production demand and gradual soil impoverishment mean that it is increasingly necessary to apply more fertilizers, mainly nitrogenous ones [2]. High quantities of nitrogen (N) per hectare need to be applied to soil to produce optimum wheat grain yields [3]. N-organic mineralization in soil is a slow process that requires the action of soil microorganisms and must also be given the necessary environmental conditions [4,5]. Plants absorb N in the form of exchangeable ammonium (NH 4 + ) or nitrates (NO 3 -), which are highly mobile

Experimental Design
A comparison of the effectiveness of the different polymeric-coated formulations was made under greenhouse conditions and then scaled to field essays. Experiment 1 (Triticum aestivum). Wheat was grown in a greenhouse at ambient temperature and humidity at the facilities of the Valencian Institute of Agrarian Research (IVIA) (Moncada, Spain) from autumn 2014 to spring 2015. Sowing was carried out in pots (22 cm high × 16 cm Ø) placed in watertight trays, sized 54 × 39 × 9 cm (6 pots/tray) with three repetitions per treatment on 7 November. The pots were filled with a non-fertilized soil from a fallow area close to the greenhouse in IVIA's grounds. Culture was developed with extra lighting for a 12:12 photoperiod and irrigated with distilled water. High temperatures were reduced and partly controlled by means of automatic systems for shading, ventilation, and water cooling. Temperature and humidity were measured every hour with a digital thermo-hygrometer. The average temperatures during the whole culture cycle, from November to May were 18.8 ± 7.2 • C (39 • C max. and 4.9 • C min.), and average relative humidity 55.3 ± 12.8% (88.0% max. and 20.4% min.). Experiment 2 (Triticum spelta). Three separated grids of 400 m 2 , divided into 16 individual surfaces of 25 m 2 each, were designed at the field level in a plot located in Teruel (Spain) at GPS coordinates 40 • 22 17.4" N, 1 • 05 58.9" W. Four treatments, including an untreated Control, were placed by quadruplicate. Sowing was carried out mechanically in the winter at a dose of 250,000 seeds/ha. Culture was surface-irrigated on a bi-weekly basis with well water.

Fertilizer Treatments
Different N fertilizers developed by Fertinagro Biotech S.L. (Teruel, Spain) were tested and their efficacy was compared. As these fertilizers are patented, the exact composition is not herein presented.
To study the influence of the different coating compositions and thicknesses, the following fertilizers were tested in Experiment 1: (1) DURAMON ® (Fertinagro Biotech S.L., Teruel, Spain), composed of urea, including a urease inhibitor (monocarbamida dihidrogenosulfate-MCDHS) with no coating (ES 2 204 307 patent); (2) a new controlled-released urea fertilizer based on DURAMON ® technology, but also 3% lignosulfonate-coated with humic acids (hereafter CRF A ); (3) the same as (2), but 5% lignosulfonate-coated with seaweed extracts (CRF B ). The three formulates had the same N composition (24-0-0), but a different coverture percentage. Fertilizers were applied in wheat at doses of 150 kg ha −1 (nitrogen fertilizer units-NFU) as a basal dressing for maximum yields based on theoretical extraction by crops. Experiment 2: Based on the physiological responses observed in the greenhouse experiment, the best CRF was selected and applied to the field at 100% doses and with a reduction to 80%. Both doses were compared with ammonium nitrosulfate (NSA, 26-0-0) (Fertinagro Biotech S.L., Teruel, Spain) as the traditional fertilizer. Maximum doses of 80 kg ha −1 were applied in the phenological state of tillering (27 April), based on the historical average yields obtained in the area. In both experiments, the same repetitions with untreated plants were included (CONTROL).

Soil Fertility Characterization
Several soil properties were measured to characterize soil fertility in both experiments. pH and EC were determined in a 1/5 (w/v) aqueous soil extract by shaking for 2 h, followed by centrifugation at 26916 g for 15 min and filtration. pH was measured by a pH meter (Crison mod.2001, Barcelona, Spain) and EC with a Conductivity meter (Crison micro CM2200, Barcelona, Spain). Total and organic soil C (SOC) and total N (N) were determined by combustion gas chromatography in a Flash EA 1112 Thermo Finnigan (Franklin, MA, USA) elemental analyzer after eliminating carbonate by acid digestion with HCl. The total nutrient contents (P, K, Ca, Mg, Cu, Fe, K, Mg, Mn, and Zn) were extracted by aqua regia digestion (3:1, v/v, HCl/HNO 3 ) and determined by ICP-AES (Thermo Elemental Iris Intrepid II XDL, Franklin, MA, USA). Analysis showed that both cultures grew on N-poor soil (Table 1). Table 1. Fertility of the soil used in the experimental analysis from the first 15 cm of soil surface. Data on total nitrogen (N), total carbon (C) and organic carbon (CO) and other macro-and micronutrients are shown. Values are means ± SD (n = 5) at the beginning of the experiment.

Plant Growth
Photosynthetically active flag leaves (PAFL) were characterized in the phenological state of panicles swelling (booting stage) by fresh weight (g) and foliar surface (cm 2 ) using a LI-3100C area meter (LI-COR ® , Lincoln, Nebraska, USA). Some plant material was weighed before being dried at 65 • C until a constant mass was obtained to determine dry weight (g). Differences in the total dry weight, length (cm), primary stem length (cm), and tillers number were determined at the end of the culture.

Leaf Greenness and Effective Quantum Yield of Photosystem II
Leaf greenness was measured in the booting stage using an SPAD-502 Chlorophyll meter (Konica-Minolta, Osaka, Japan). The effective quantum yield of photosystem II electron transport (ΦPSII), which represents the electron transport efficiency between photosystems within light-adapted leaves, was checked with a leaf fluorometer (Fluorpen FP100, Photos System Instrument, Drásov, Czech Republic). Both parameters were measured in a minimum of 25 PAFL.

Gas Exchange Analysis
Gas exchange measurements were taken at noon in five plants per treatment using a portable infrared gas analyzer LCpro-SD, equipped with a PLU5 LED light unit (ADC BioScientific Ltd., Hoddesdon, UK). The selected flag leaves in wheat (booting stage) of Experiment 1 were analyzed in a leaf chamber (6.25 cm 2 ) to determine the following parameters: stomatal conductance (gs) (expressed as mmol m −2 s −1 ), net photosynthetic rate (A) (µmol m −2 s −1 ), transpiration (E) (mol m −2 s −1 ), and intercellular CO 2 concentration (Ci) (µmol mol −1 ) under ambient CO 2 , temperature, and relative humidity conditions. They were recorded by programming increasing photosynthetically active radiations (PAR) of 348, 522, 696, 870, 1218, and 1566 µmol m −2 s −1 . Water-use efficiency (WUE) and intrinsic WUE were calculated as the ratio between A/gs and A/E, expressed as µmol (CO 2 assimilated) mol −1 (H 2 O transpired). per treatment and culture were collected and kept at −80 • C until they were biochemically analyzed. The compositions in macro-(N, P, K, Ca, Mg and S) and micronutrients (Fe, Cu, Mn, Zn, B, and Mo) were determined by Inductively Coupled Plasma Optical Emission Spectrometry (ICP-OES). N content was estimated by an N-Pen N 100 apparatus (Photon System Instruments, Drásov, Czech Republic).

Growth, Yield, and Cereal Grain Composition
Once grain ripening had been completed at 138 days in greenhouses after emergence, the remaining plants per culture were harvested and characterized in growth and grain yield terms. The growth parameters of the total dry weight of aerial parts, primary stem length, tillers number, ears number, and ear length and weight were measured. Yield was determined by measuring the total dry grain weight, one-hundred grain weight and grain number. Nitrogen Use Efficiency (NUE) for each fertilizer treatment was calculated as agronomic efficiency according to [37]: NUE kg kg− 1 = Grain yield of the fertilized area − Grain yield of the unfertilized area Quantity of N applied as N fertilizer (1) At the field level and 90 days after applying fertilizers, the biomass of aerial parts, ears weight, and grain weight was studied. The Harvest Index (HI) was calculated as the grain weight/biomass of aerial parts. Different quality parameters were also measured in the grain. A representative composite sample was prepared separately for each treatment, pooling fractions of plant material for each replication. Subsequently, each composite mixture was ground and analyzed based on food quality analysis methods (Comission Regulation EC Nº 152/2009 of 27 January): humidity (gravimetric by drying in an oven at 130 • C), ashes (gravimetric by incineration at 550 • C), lipids (extraction without hydrolisis in Soxtec Avanti-Foss), protein (Kjeldahl method using Foss automatic distillation equipment, Foss, Hillerød, Denmark), crude fiber (gravimetric), and total carbohydrates (volumetric using Luff Schoorl reagent). Analysis were carried out by the Valencia's Agrifood Laboratory (Burjassot, Spain).

Statistics
The differences between fertilizers treatments were tested by analysis of variance (ANOVA) at 95% confidence. Prior to the ANOVA, the data requirements of normality and homogeneity of variances were checked according to Levene's and Shapiro-Wilk tests. When the null ANOVA hypothesis was rejected, post hoc comparisons were made to establish the possible statistical differences among the different treatments applied using the Fisher's LSD test. The statistical Statgraphics Centurion XV, version 15.2.05 software program (Statpoint Technologies, Inc., Warrenton, VA, USA) was used to perform the analysis.

Plant Growth, Leaf Greenness and Effective Quantum Yield of Photosystem II
No significant differences were observed for the various treatments performed for chlorophyll content and ΦPSII, measured by nondestructive techniques in the phenological state of panicles swelling (booting stage), although they were significant compared to the CONTROL (Table 2). Similar results were obtained when studying PAFL fresh weight content, dry weight and area. The total fresh weight of aerial parts was 1.4-and 1.7-fold significantly higher for the plants fertilized with DURAMON ® compared to CRF A and CRF B , but no differences were observed for dry weight. The responses of fertilizer treatments on plant growth, foliar area, and root development are shown in Figure 1. Table 2. Effects of fertilizer treatments CRF A , CRF B, and DURAMON ® on photosynthetic parameters (effective quantum yield of photosystem II-ΦPSII, leaf greenness, nitrogen content, total fresh weight of aerial parts (g), dry weight of aerial parts (%), photosynthetically active flag leaves-PAFL fresh weight, PAFL area, and leaf area index-LAI) in Triticum aestivum leaves compared to the Control in the phenological state of panicles swelling (booting state). Values are means ± SD (n = 18 for ΦPSII, leaf greenness, and N content; n = 8 for the other growth parameters).

Parameters CRF
Different letters in the same row indicate significant statistical differences (Fisher's LSD test, P < 0.05).
Agronomy 2020, 10, x FOR PEER REVIEW 6 of 15 The responses of fertilizer treatments on plant growth, foliar area, and root development are shown in Figure 1.

Gas Exchange Analysis
Significant increases in gs, A, and E were noted in the plants treated with the different fertilizers as increasing PAR levels were applied from 348 to 1566 µmol m −2 s −1 (Figure 2A-C). Conversely, Ci showed a decreasing tendency in all the applied treatments ( Figure 2D).
between CRFs and DURAMON ® . After the maximum PAR application of 1566 µmol m −2 s −1 , the levels of A did not significantly differ for both CRFs, but gs and Ci were 1.3-and 1.1-fold significantly higher for CRFA than for CRFB. At the same PAR, the A levels were 1.3-fold significantly higher for CRFs than for DURAMON ® . The gs and A levels for the CONTROL plants were significantly lower at all the studied PAR compared to those of the different treatments. However, the E and Ci levels for the CONTROL plants were only significantly different for DURAMON ® , as was Ci with CRFA. A/gs significantly differed when globally comparing CRFs with DURAMON ® and the CONTROL ( Figure 3A), but the A/E levels were significantly higher for CRFA than for CRFB, DURAMON ® and the CONTROL, which also significantly differed from one another. The maximum A/E levels were produced within the 870 to 1218 µmol m −2 s −1 range ( Figure 3B). The A/E levels in the CRFA leaves were 1.3-and 1.8-fold significantly higher than CRFB and DURAMON ® at a PAR of 870 µmol m −2 s −1 .  The levels of gs, A and Ci for CRF A were higher than those found for CRF B and DURAMON ® but were lower for E. Significant differences were found in all the studied gas exchange parameters between CRFs and DURAMON ® . After the maximum PAR application of 1566 µmol m −2 s −1 , the levels of A did not significantly differ for both CRFs, but gs and Ci were 1.3-and 1.1-fold significantly higher for CRF A than for CRF B . At the same PAR, the A levels were 1.3-fold significantly higher for CRFs than for DURAMON ® . The gs and A levels for the CONTROL plants were significantly lower at all the studied PAR compared to those of the different treatments. However, the E and Ci levels for the CONTROL plants were only significantly different for DURAMON ® , as was Ci with CRF A . A/gs significantly differed when globally comparing CRFs with DURAMON ® and the CONTROL ( Figure 3A), but the A/E levels were significantly higher for CRF A than for CRF B , DURAMON ® and the CONTROL, which also significantly differed from one another. The maximum A/E levels were produced within the 870 to 1218 µmol m −2 s −1 range ( Figure 3B). The A/E levels in the CRF A leaves were 1.3-and 1.8-fold significantly higher than CRF B and DURAMON ® at a PAR of 870 µmol m −2 s −1 . Agronomy 2020, 10, x FOR PEER REVIEW 8 of 15

Foliar Nutrient Content
No significant differences appeared in the PAFL macronutrient concentrations of N, K, Ca, Mg, and S at the beginning of ears formation ( Figure 4A). The plants fertilized with CRFA presented 1.1and 1.2-fold higher foliar N average levels than CRFB and DURAMON ® . The foliar P concentrations in the plants fertilized with CRFA were 1.2-fold significantly higher than for CRFB and DURAMON ® . On average, the plants fertilized with CRFB presented slightly higher contents of K, Ca, Mg, and S. Regarding micronutrient foliar contents, the plants fertilized with CRFA presented 1.7-and 1.5-fold higher Fe levels than CRFB and DURAMON ® , respectively ( Figure 4B). Cu and Zn contents were 1.4and 1.2-fold significantly higher in the plants fertilized with CRFs than in those fertilized with DURAMON ® , respectively. Foliar Mn concentrations did not differ significantly among the distinct fertilizer treatments and the CONTROL levels came close to the critical thresholds. The foliar B levels were significantly higher in the plants treated with CRFB compared to those treated with DURAMON ® , but the quantitative CONTROL levels came close to CRFB. Mo in PAFL content was < 2 mg kg -1 DW (dry weight) in all the treatments and the CONTROL.
No clear correspondence was obtained for the macro-and micronutrient concentrations quantified at the foliar level compared to those found in roots ( Figure 4C, 4D). Quantitatively, the root N concentrations were around half of those obtained at the foliar level. The remaining macronutrient contents were slightly lower in roots, except for Ca. Micronutrients almost doubled in roots compared to foliar content but were 30-fold higher in Fe content. No significant differences were obtained for the three compared treatments in the macro-and micronutrient contents at the root level, except for Ca and Mg, which were higher in the plants treated with CRFs compared to DURAMON ® . Values represent means ± SE (n = 6). Different letters for each treatment (same color) indicate statistically significant differences (ANOVA, P < 0.05).

Foliar Nutrient Content
No significant differences appeared in the PAFL macronutrient concentrations of N, K, Ca, Mg, and S at the beginning of ears formation ( Figure 4A). The plants fertilized with CRF A presented 1.1and 1.2-fold higher foliar N average levels than CRF B and DURAMON ® . The foliar P concentrations in the plants fertilized with CRF A were 1.2-fold significantly higher than for CRF B and DURAMON ® . On average, the plants fertilized with CRF B presented slightly higher contents of K, Ca, Mg, and S. Regarding micronutrient foliar contents, the plants fertilized with CRF A presented 1.7-and 1.5-fold higher Fe levels than CRF B and DURAMON ® , respectively ( Figure 4B). Cu and Zn contents were 1.4-and 1.2-fold significantly higher in the plants fertilized with CRFs than in those fertilized with DURAMON ® , respectively. Foliar Mn concentrations did not differ significantly among the distinct fertilizer treatments and the CONTROL levels came close to the critical thresholds. The foliar B levels were significantly higher in the plants treated with CRF B compared to those treated with DURAMON ® , but the quantitative CONTROL levels came close to CRF B . Mo in PAFL content was < 2 mg kg -1 DW (dry weight) in all the treatments and the CONTROL.
No clear correspondence was obtained for the macro-and micronutrient concentrations quantified at the foliar level compared to those found in roots ( Figure 4C, 4D). Quantitatively, the root N concentrations were around half of those obtained at the foliar level. The remaining macronutrient contents were slightly lower in roots, except for Ca. Micronutrients almost doubled in roots compared to foliar content but were 30-fold higher in Fe content. No significant differences were obtained for the three compared treatments in the macro-and micronutrient contents at the root level, except for Ca and Mg, which were higher in the plants treated with CRFs compared to DURAMON ® .

Growth, Yield, and Cereal Grain Composition
The results about the measured growth and yield parameters in greenhouses are shown in Table 3. No significant differences were found in the different treatments for the measured growth parameters. Significant differences were observed in the CONTROL with the CRFA-and DURAMON ® -treated plants for the total dry weight of aerial parts, and with DURAMON ® for ear weight and number. Regarding the yield parameters, total dry grain weight was significantly higher for CRFA than for CRFB. No significant differences were observed in grain weight among treatments, but differences were significant compared to the CONTROL. Total grain number was 1.2-fold significantly higher for CRFA than for CRFB but was not significant compared to DURAMON ® . At the field level, no significant differences were found in the various applied treatments when comparing growth and grain yield parameters (Table 4). However, NUE was 28 and 38% higher for CRFA 80%, compared to CRFA 100% and NSA; these differences were statistically significant. No significant differences were observed among treatments and the CONTROL for the analyzed grain parameters. As a result, on average the values were 1.6% CRFA ash, 12.9% humidity, 1.7% lipids, 11.4% protein, 2.5% crude fiber, and 65.8% total carbohydrates. No significant differences were obtained comparing the studied quality parameters between CRFA and the rest of treatments and CONTROL.

Growth, Yield, and Cereal Grain Composition
The results about the measured growth and yield parameters in greenhouses are shown in Table 3. No significant differences were found in the different treatments for the measured growth parameters. Significant differences were observed in the CONTROL with the CRF A -and DURAMON ® -treated plants for the total dry weight of aerial parts, and with DURAMON ® for ear weight and number. Regarding the yield parameters, total dry grain weight was significantly higher for CRF A than for CRF B . No significant differences were observed in grain weight among treatments, but differences were significant compared to the CONTROL. Total grain number was 1.2-fold significantly higher for CRF A than for CRF B but was not significant compared to DURAMON ® . At the field level, no significant differences were found in the various applied treatments when comparing growth and grain yield parameters (Table 4). However, NUE was 28 and 38% higher for CRF A 80%, compared to CRF A 100% and NSA; these differences were statistically significant. No significant differences were observed among treatments and the CONTROL for the analyzed grain parameters. As a result, on average the values were 1.6% CRF A ash, 12.9% humidity, 1.7% lipids, 11.4% protein, 2.5% crude fiber, and 65.8% total carbohydrates. No significant differences were obtained comparing the studied quality parameters between CRF A and the rest of treatments and CONTROL.

Discussion
Nutrient release contained in fertilizers depends on many factors, including environmental conditions, crop management, and the chemical composition of fertilizers [2,38,39]. Nowadays, the slowing down of nutrient emissions in fertilization is a challenge that has already been overcome [11,18,24]. However, high production costs and contamination linked to synthetic fertilizers and waste materials, together with increasingly restrictive environmental policies, have forced new ecological materials to be sought to allow sustainable fertilization [40,41]. The use of water-soluble synthetic products or natural polymers based on lignin has been indicated as an alternative to these problems because they can be obtained in large quantities and at cheap prices from the waste generated in the paper industry, from wood and other sources [42,43].
Research conducted with CRFs has shown that their efficiency is generally higher to that of traditional fertilizers and SRFs [9]. In fact, SRFs are more sensitive to high temperature and sandy soils [44]. Nevertheless, most research works conducted to date with CRFs have focused mainly on crops with a high added value, such as horticultural, ornamental or wood products, and have obtained different results [45][46][47][48][49][50][51][52]. It is important to point out that the main challenge of CRFs application to crops is to successfully provide the amount of nutrients that plants need and in a fractional manner. Moreover, there are also the important advantages that CRFs offer 'per se' in both crop management and the environment. In fact, CRF applications are usually unique, which means savings in crop handling from avoiding successive top-dressing fertilizer applications. Finally, nutrient doses are usually lower in CRFs than those applied with traditional fertilizers, and N losses by evaporation or leaching consistently lower.
In two experiments, this research compares the effectiveness of two lignosulfonate-based polymer-coated urea fertilizers: an analogous non-coated urea, and ammonium nitrosulfate as a traditional fertilizer. Based on the experimental design, the CRF with the best behavior was selected based on the responses of wheat to different physiological and yield parameters on the microscale. In a second stage, the selected CRF was compared with ammonium nitrosulfate in the field. In physiological terms, significant differences were found in growth, chlorophyll content and ΦPSII among CRF A , CRF B and DURAMON ® compared with the untreated plants in early crop development stages. Lower values of E and higher of Ci were detected in plants treated with DURAMON ® , as compared to the CRFs treatments. Intrinsic efficiency in the water use of CRF A was significantly higher compared to CRF B and DURAMON ® . To explain these results, the enhanced effects of CRFs were produced by a combination of the individual effects of lignosulfonates and biostimulants on nutrient supply. In fact, lignosulfonates or sulfonated lignin have a variety of functional groups that provide unique colloidal properties and act as chelating agents [53]. The humic substances contained in CRF A may promote plant development by stimulating root and shoot growth as they can enhance nutrient use efficiency by facilitating the assimilation of macro-and microelements [54,55]. Seaweed extracts can enhance chlorophyll and carotenoid contents in plant shoots, root thickness, and biomass [56]. However, the effectiveness of marketed biostimulants depends very much on their origin because their composition and proportions usually change [57].
Better physiological responses in the state of panicles formation suggest that plants would increase yield and biomass parameters at the end of the crop. A significant correlation was also found between higher levels of photosynthesis during grain formation and increased crop yields obtained in wheat [58,59]. In fact, on a microscale, the total yield expressed in dry grain weight was significantly higher in CRF A than in CRF B ; and was, on average, also higher than DURAMON ® . The observed differences were due to a large number of grains harvested by the production of 1.1-fold more spikes in the plants fertilized with CRF A compared to CRF B and DURAMON ® . On average the plants fertilized with DURAMON ® produced larger sized and heavier spikes with more grains. This could be explained by the faster N-release kinetics of DURAMON ® compared to CRFs, as lacking lignosulfonate-polymeric coverage. A bigger N supply in the first crop stages could explain why the plants fertilized with DURAMON ® seemed to be slightly more advanced in their phenological status compared to CRFs. Physiological requirements of wheat are established as 3 kg N Qm −1 grain; therefore, the theoretical yield that should have been obtained at a dose of 150 kg N ha −1 was 5 t ha −1 . Our results showed exceeded yields of 6.4 grain t ha −1 with CRF A , 5.3 grain t ha −1 with CRF B and 5.8 grain t ha −1 with DURAMON ® . Despite DURAMON ® not being a CRF because it lacks polymer-coating, it could be considered an SRF for being formulated with urease inhibitor MCDHS, which is also contained in CRFs. This would explain why the DURAMON ® -treated plants gave yields close to CRFs. In fact, the minor differences between DURAMON ® and CRFs might indicate that DURAMON ® could be also used successfully to maintain N availability for plants over time in wheat. As examples, when using nitrification and urease inhibitors in wheat, maize and barley, it was obtained better crop yields and N 2 O mitigation than SRCFs [2]. Further, better performance for CRFs and those formulated with nitrification inhibitors compared to traditional ones in maize with a reduction in N 2 O emissions up to 21% that did not affect yields [39]. The lower yields obtained with CRF B could be explained by excessive N emission slowdown by having formulated with a 2% thicker polymeric coverage. The best results obtained with CRF A were confirmed at field level when comparing the NUE between the applied fertilizer treatments. It was possible to obtain yields close to those observed with CRF A 100% and NSA, by applying CRF A with a 20% less N content, but significantly increasing the NUE by more than four times. Even though it has been reported that NUE can vary depending on factors like the doses applied or climatic conditions [2], no reductions in grain yield were observed when applying different CRF formulations by reducing N content in a similar proportion [60,61].
The macronutrients analysis showed that wheat plants had a good NPK nutritional status in phenological state at the beginning of ear formation. The N concentrations in the treated plants fell within the N leaf DW 4-6% range, which is considered suitable for obtaining high yields [62], but no statistically significant differences appeared in the applied treatments. On the contrary, P concentrations were 1.2-fold significantly higher in CRF A than in CRF B and DURAMON ® . In all cases, P levels fell within the range considered optimal for good plant development (0.2-0.5% leaf DW). Although the applied N fertilizer was not mixed with P, it is known that soil N applications can stimulate root growth and increase cation exchange capacity to favor Ca uptake, and P uptake indirectly [63]. In cereals, increased yields and improvements in the content of macro-and micronutrients of crops have been achieved in barley, maize, rice, or wheat using SCRFs [15,[64][65][66][67][68][69]. No significant differences were found in the K, Ca, Mg, and S contents in leaves among treatments, and their levels were medium to high. Regarding micronutrients, the plants fertilized with CRF A presented 1.7-and 1.5-fold higher Fe levels than CRF B and DURAMON ® , respectively. However, the Fe levels were optimum for maximum yields (21-200 mg kg −1 leaf DW). Cu and Zn contents were 1.4-and 1.2-fold significantly higher in the plants fertilized with CRFs compared to those fertilized with DURAMON ® , but concentrations were at the lowest levels within the range considered normal for Cu and Zn (5-50 and 20-70 mg kg −1 foliar DW, respectively). The Mn content did not differ significantly for the different fertilizer treatments and presented lower levels (16-200 mg kg −1 leaf DW). Fertilizer treatments did not significantly affect B content as the CONTROL plants had similar levels. The Mo levels were very low and were lower than 2 mg kg −1 leaf DW in all the treatments.

Conclusions
In the present study, we have carried out a comparison of two lignin-coated controlled release fertilizers enriched with humic substances (CRF A ) or seaweed extracts (CRF B ) with a similar non polymeric-coated fertilizer (DURAMON ® ) and with an ammonium nitrosulfate one (NSA). Our results showed that plants performed better when they were fertilized with CRFs coated with humic substances, although the improvement in crop yield was not excessive compared to the seaweed-coated one and that the uncoated one. However, a significant improvement in crop yield and the measured physiological parameters of wheat plants was achieved with Fertinagro's controlled release fertilizers compared to the traditional NSA. Fertilization with these new technified CRFs greatly favored the wheat crop management by making it possible to carry out one single application as a basal dressing, which simplified crop handling. Smaller amounts of N in formulations gave important advantages, such as reduced costs and minimized N losses, which thus avoids the common contamination problems that usually occur when applying traditional fertilizers. We conclude that it is possible to use this technology in extensive cropping, as lignin-based polymers are economically feasible and environmentally friendly.