Assessing Nitrogen Use Efficiency and Nitrogen Loss in a Forage-Based System Using a Modeling Approach

In intensive agriculture, N supply often exceeds crop requirements, even in nitrate vulnerable zones (NVZ). In farmland, the N surplus gives rise to NO3 ́ leaching and consequent groundwater pollution. The present study aimed at proposing measures to reduce N leaching and hence improve N efficiency in a buffalo livestock farm located in the NVZ of Latina plain (Central Italy). The farm was cultivated with forage crops in a double annual crop rotation: Italian ryegrass (Lolium multiflorum Lam.) in winter and silage corn (Zea mays L.) in summer. Mineral and organic fertilizers were supplied to both crops. The annual N budget and soil solution NO3-N concentrations were evaluated using a modeling approach. The performance of the WinEPIC model in simulating the response of the NO3-N concentration in percolation to the N application rate was assessed and validated by field measurements of the NO3-N concentration in the soil solution. Three scenarios were proposed to identify the best practice to minimize the environmental impact of N application without significant yield loss. Also, recommendations of best practices in N fertilization and animal manure spreading were given. This study thus provides useful preliminary information for decision-making in agriculture/environmental policies.


Introduction
In the majority of European Member States, the agricultural sector contributes over 50% of the total nitrate (NO 3 ´) leaching to surface and groundwater [1].NO 3 ´leaching from agricultural soils is a complex process closely related to local environmental factors such as soil characteristics and climatic variables [2,3], and to farm management practices in intensive agriculture.Here, inappropriate application of manure, fertilizer, and sewage sludge [4], or irrigation and planting patterns [5], are liable to affect nitrate leaching in both livestock and arable farms [6], as well as leading to erosion and surface runoff.As in Knudsen et al. [7] and in Kros et al. [8], the proportion of soil N loss ranged between 30% and 60%, mainly leached below the root zone.NO 3 ´leaching and groundwater pollution are significant public concerns worldwide [9] and are often linked to livestock farms [10,11].As reported in Colombo et al. [12], more than 85% of the European agricultural land exceeds the threshold limit of NO 3 ´levels (25 mg¨L ´1) due to N leaching, and approximately 22% of the groundwater collected in these areas overtakes the maximum admissible NO 3 ´concentration (50 mg¨L ´1) recommended by the World Health Organization for drinking water [13,14].Managing this pollution is challenging because it profoundly affects the nitrogen (N) cycle [15] and has strong environmental impact due to its direct influence on drinking water, eutrophication and its indirect contribution to atmospheric emissions with ammonia or nitrogen oxides [5,16,17].
Thus, in the last decades the European Union (EU) has implemented political strategies aiming at reducing and controlling NO 3 ´contamination of groundwater caused by N losses from agricultural land [12,18].Among these, the Nitrates Directive [19] aims to reduce water contamination caused by NO 3 ´.It legally restricts farm manure application to 170 kg N ha ´1 [20], addressing farmers to get them to comply with the water quality protection measures.The directive also promotes the management action programs imposed in nitrate vulnerable zones (NVZs) [3,21] and it is strictly linked to other EU policies focused on air and water quality, climate change and agriculture.NO 3 ´pollution is a relevant issue for intensive livestock farms where liquid and solid manure as organic effluents are used [22][23][24][25].NO 3 ´leaching is particularly serious in Mediterranean semi-arid countries [23,26], where substantial N losses by leaching were found to occur during the rainy period in winter and after irrigation practice in summer.To date, the relationships between soil NO 3 ´accumulation and groundwater NO 3 ćoncentrations in agricultural soils are still not entirely understood [1].Over the last decade, numerical models linked to N cycling in soil water and plants have received increasing attention as useful tools to predict the risk of NO 3 -N contamination in surface and groundwater and to define environmentally and economically suitable agricultural systems [27][28][29][30][31]. Models, when properly validated, provide a fast and cost-effective way of estimating NO 3 -N leaching under different agricultural management practices.Another frequently used approach to quantify N use efficiency at various scales is the N balance [32][33][34][35], widely recognized as a powerful tool employed by policy-makers for raising farmers' awareness of their current management practices [20,34] and to support sustainable farming planning [17,36] at field [37], farm-gate [38] or territorial scale [39].
In an intensive buffalo livestock farm located in the NVZ of Latina plain (Central Italy) cultivated with forage crops in a double annual crop rotation-Italian ryegrass (Lolium multiflorum Lam.) in winter and silage corn (Zea mays L.) in summer-a study was carried out addressing the following issues: (1) monitoring of NO 3 ´concentrations in the soil solution, (2) evaluation of the EPIC model against the data in terms of its ability to simulate the process of NO 3 -N leaching loss in field conditions, and (3) proposal of alternative fertilizing management strategies as possible measures to reduce N leaching and to improve the N use efficiency in the Mediterranean conditions.

Site Description and Data Collection
In the framework of a local research project (Nitrate Vulnerability of Latium Region-VULNRELA), a monitoring study about NO 3 ´percolation was performed in a private buffalo livestock farm located in Sabaudia, Latina Province plain, Central Italy.The farm, covering an area of about 17 hectares, is located in a coastal plain with mean altitude of 9 m a.s.l., drained in the 1930s from wetlands.Outcropping lithologies are mixed yellow sands, tuffs and lacustrine clays, and a multi-level aquifer system is present; the plain is included in a NVZ for the presence of intensive agriculture and livestock farming (Figure 1).
The mean temperatures during winter and summer are 10.5 and 20.5 ˝C, respectively, the mean annual rainfall is about 800 mm, and the rainy season-with about 60% of the total annual rainfall occurring-is November-March.Temperate winter and dry summer are typical of the Mediterranean climate (Figure 2).For the purposes of the research, time series of climatic data (2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015) were acquired from a weather station managed by ARSIAL (Regional Agency for Extension Services of Latium), about 7 km from the experimental site, which provided records of rainfall, relative humidity, air temperatures, wind speed and solar radiation.The mean temperatures during winter and summer are 10.5 and 20.5 °C, respectively, the mean annual rainfall is about 800 mm, and the rainy season-with about 60% of the total annual rainfall occurring-is November-March.Temperate winter and dry summer are typical of the Mediterranean climate (Figure 2).For the purposes of the research, time series of climatic data (2004-2015) were acquired from a weather station managed by ARSIAL (Regional Agency for Extension Services of Latium), about 7 km from the experimental site, which provided records of rainfall, relative humidity, air temperatures, wind speed and solar radiation.N fertilization in the study area often exceeds the limit of 170 kg ha −1 fixed for NVZs by the EU,  The mean temperatures during winter and summer are 10.5 and 20.5 °C, respectively, the mean annual rainfall is about 800 mm, and the rainy season-with about 60% of the total annual rainfall occurring-is November-March.Temperate winter and dry summer are typical of the Mediterranean climate (Figure 2).For the purposes of the research, time series of climatic data (2004-2015) were acquired from a weather station managed by ARSIAL (Regional Agency for Extension Services of Latium), about 7 km from the experimental site, which provided records of rainfall, relative humidity, air temperatures, wind speed and solar radiation.N fertilization in the study area often exceeds the limit of 170 kg ha −1 fixed for NVZs by the EU, and gross irrigation supplies can reach a total of 600 mm for each summer cultivation cycle.The frequent and large use of fertilizers, and the intensive irrigation combined with a shallow groundwater table create a favorable environment for the leaching of nutrients to subsoil layers and groundwater.
The soil of the experimental site, classified as Vertic Stagnosols [40], was analyzed for physical, hydrological, and chemical properties according to the standard methods [41,42].In Table 1, the main soil physical and chemical characteristics are reported.N fertilization in the study area often exceeds the limit of 170 kg¨ha ´1 fixed for NVZs by the EU, and gross irrigation supplies can reach a total of 600 mm for each summer cultivation cycle.The frequent and large use of fertilizers, and the intensive irrigation combined with a shallow groundwater table create a favorable environment for the leaching of nutrients to subsoil layers and groundwater.
The soil of the experimental site, classified as Vertic Stagnosols [40], was analyzed for physical, hydrological, and chemical properties according to the standard methods [41,42].In Table 1, the main soil physical and chemical characteristics are reported.At present in the farm there are 430 heads of buffalo.The monitored field is cultivated with forage crops in a double crop annual rotation: Italian ryegrass (Lolium multiflorum Lam.) grown from October to May and silage corn (Zea mays L.) grown from June to September (Table 2).Organic fertilizers are provided by livestock as buffalo slurry and manure and are applied in combination with Nitrophoska® (12% N) prior to sowing each crop.Urea 46% is provided as a top dressing.Considering an average N content of the liquid and solid manure of 1.96 g kg ´1 [43], the amount of organic N applied yearly to the field is about 157 kg¨ha ´1.During the silage corn cycle (dry season), sprinkler irrigation is provided weekly at an average rate of 50 mm.All the information was provided by the farmer.The monitored field was equipped with three micro-lysimeters, made up by a PVC pipe with a porous ceramic cup at one end.These were inserted in the soil to collect soil solution down to 30, 60 and 90 cm depths, respectively.By means of a vacuum pump, a depression is created inside the lysimeters, establishing a pressure gradient from soil to the inner lysimeter that allows the soil solution to enter the pipe through the porous ceramic cup.Samples of this solution were collected every about two weeks during a two-year period, between November 2013 and November 2015, transferred in plastic vials and stored at ´18 ˝C to prevent any N transformation before analysis.
At a later stage these samples were defrosted and analyzed for NO 3 -N concentration by means of a spectrophotometer with an automatic procedure: nitrates are reduced into nitrites that form a colored complex read at 550 nm.Due to technical problems occurred during summer 2015, the measurements are considered as reliable until May 2015, and hence later measurements were not reported.

Nitrogen Balance and Nitrogen Use Efficiency (NUE) Calculation
A simplified N balance in the soil-plant system was computed following the approach proposed by Constantin et al. [44], Pieri et al. [45], and Di Bene et al. [46].The method accounts for the input and output components, calculated as the difference between the N entering and leaving the soil.The balance is an indicator of either surplus (when input > output), indicating N accumulation in the soil, or deficit (when output > input), showing a net N depletion from the soil; the result is expressed as kg N ha ´1 year ´1.To consider the N surplus as an indicator of possible N leaching is a simplification; however, the available measured data did not allow us to take into account the other components of N loss.
The N balance was calculated as follows: where N balance refers to total N variation; A is N from atmospheric deposition; F is N applied as mineral and organic fertilizers; R is N input from organic sources (plant roots); U is the N amount taken up by the total above-ground crop biomass.
A was estimated by multiplying the average N concentration in dry and wet atmospheric depositions for the yearly rainfall amount.According to Grignani et al. [33] and Pieri et al. [45], the N concentration in dry and wet atmospheric depositions was, on average, 2 mg¨L ´1 for Central Italy.F was calculated by multiplying the amount of fertilizers applied for their respective N concentrations (Table 2).R was estimated as 15% of the total above-ground biomass of each crop [47].U was calculated by multiplying the total above-ground biomass of each crop by its N concentration (Table 2).
In order to assess the nitrogen use efficiency (NUE) of the Italian ryegrass-silage corn rotation, the partial nutrient balance proposed by Dobermann et al. [48] was used: This is the simplest form of NUE: a value close to 1 suggests that soil fertility will be sustained at a steady state, while values well below 1 suggest avoidable nutrient losses.Also values >1 must be regarded as unsustainable, because more nutrients are removed with the harvested crop than applied by fertilizer and/or manure.This is a simplification of a complex process, since nutrient removal by erosion, leaching and gaseous emissions are not considered.Nevertheless, this index provides qualitative information about the N available to the crops and the environmental impact in terms of N surplus.The obtained NUE allows farmers to make comparisons among crops in the same environment and time, and to decide which crop or which form of agronomic management is advisable.

Model Simulation
Data about soil, weather and agronomic management were used as input for the EPIC model-WinEPIC0509 version 1.0, a field-scale agroecosystem model that simulates crop production as a function of weather, soil conditions, and production practices employed (e.g., tillage types, tillage frequency, crop rotations).EPIC was designed, in principle, to explore the impacts of soil erosion on crop productivity [49], and evolved with continued refinements to approach carbon and nutrient cycling via submodels; also, additional capacity was introduced to predict water quality and the response of crops to atmospheric CO 2 [50].EPIC and the models which have evolved from it have been applied extensively to a variety of soils and cropping systems worldwide [51].
EPIC's nutrient cycling sub-routine tracks soil N movement through five organic C/N pools and two inorganic pools.The five organic C/N pools consist of two surface litter pools-metabolic and structural litter-and three subsurface pools-microbial biomass, an active soil humus pool (called "slow" humus) and a stable (or "passive") soil humus pool.The inorganic N pools are NO 3 ´and NH 4 + .
N moves into the organic pools through plant residues, roots, green manures or fertilizers, and out of the organic pools through sediment losses or through mineralization to the inorganic ammonium pool, which also receives inputs from ammonium fertilizers.The ammonium pool increases based on mineralization and fertilizer inputs, and decreases via nitrification losses to the NO 3 ´pool or volatilization into the atmosphere as NH 3 .
The NO 3 ´pool in EPIC represents the pool of N that is available for plant uptake.N enters this pool through precipitation deposition, fertilization and nitrification, and exits through runoff, subsurface flow, leaching, denitrification and plant uptake [52].A thorough description of the functions underlying the various transformation processes in the organic C/N cycle can be found in Izaurralde et al. [53].
In the considered farm, the N cycle was modeled for the same period in which climatic data are available (2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015), assuming no changes of agronomic management.In order to compare the modeled NO 3 ´concentration in soil solution for the whole soil profile to the NO 3 ´concentration measured into the lysimeters, only the period from November 1 2013 to September 30 2015 was reported as a result.In WinEPIC outputs, the NO 3 ´concentration is given as the daily concentration in stored water along the soil profile, while we measured the soil solution concentration at three different depths at specific sampling dates (every about two weeks).Rather than comparing these two types of data, we considered it more appropriate to compare their trends, in order to assess whether predicted (modeled) NO 3 ´concentrations were consistent with measured NO 3 ´concentrations, especially when the critical (highest) NO 3 ´concentrations in soil solution occurred.Therefore, at each sampling date the measured cumulative concentration was put in relation to the sum of the WinEPIC daily values of a two-week period preceding the same dates.
In order to evaluate the efficiency of a change in the N fertilization rate, the WinEPIC model was used to carry out a scenario analysis for the 2004-2015 period.The current N fertilizing management was used as a baseline (S 0 ) for simulating three alternative scenarios with different fertilization rates.The first scenario (S 1 ) was run with a reduced N supply considering the high N fertilization in the baseline; in the other scenarios (S 2 and S 3 ) we assumed that the mineral fertilization supplies more readily available N, in comparison with the organic fertilization, when applied during the rainy season or an irrigation period, thus increasing the amount of N leaching.Scenarios are defined as follows: S 0 : Baseline; S 1 : 50% reduction of both organic and mineral fertilization rates, maintaining the same application dates; S 2 : omission of the mineral fertilization in Italian ryegrass cycle only; S 3 : elimination of the mineral fertilization on both Italian ryegrass and silage corn cycles.
These could represent possible options for farmers for more sustainable fertilizing management, mitigating the N leaching risk.

Nitrogen Balance and Nitrogen Use Efficiency
The N balance was obtained as the difference between the mean annual N inputs and the mean annual N outputs.It was calculated separately for each crop and then summed up for the entire cropping system, as rotational crops can take advantage of previous applications of fertilizers and residual plant inputs (Table 3).
The average annual N balance was equal to 224 kg¨ha ´1 with the Italian ryegrass management contributing 65% to this amount.This N surplus is mainly due to the high N application of manure and mineral fertilizations that largely exceeded N removals by plant uptake, and it is the maximum amount allowed by the Nitrates Directive (see Table 2).The NUE index value for the whole crop rotation indicates that, yearly, the above-ground biomass crop absorbed only 49% of the total amount of N, originating from A (atmospheric deposition), F (mineral and organic fertilizers), R (plant root additions).Likewise, Bassanino et al. [34] and Demurtas et al. [31] obtained NUE values ranging between 47% and 50% in similar climatic and agronomic conditions with cattle livestock.In particular, for the considered cropping system, the Italian ryegrass had a much lower efficiency than silage corn in converting N input into crop N biomass.The N surplus may constitute an indicator for the agricultural pressure on the environment and of N loss, which is driven mainly by the N fertilization rate and the soil's physical, chemical and microbiological properties.The N given to the soil as NH 4 + can be quickly converted to NO 3 ´by the nitrifying microorganisms.This anion is not absorbed by the soil, and thus is easily released into the soil liquid phase, possibly moving to the leaching water flow.Sandy soils are particularly sensitive to the leaching of nitrate to groundwater due to their higher permeability.* N balance: total N variation; F: N applied as mineral and organic fertilizers; R: N input by organic addition (plant roots); U: N amount taken up by the above-ground crop biomass; A: N from atmospheric deposition, estimated by multiplying the N concentration for the yearly rainfall amounts.

Nitrogen Leaching and WinEPIC Simulation
Based on the field data, WinEPIC model was used to simulate the soil N cycle.Among the various components of this cycle, we focused on the average total soluble N and on the NO 3 ´concentration in the soil profile, parameters somewhat involved in the evaluation of groundwater pollution risk for NVZ.
Average total soluble N in the soil profile estimated by EPIC is 198 kg¨ha ´1, which represents 88% of the N surplus obtained by the calculated balance.This fraction is consistent with the observations by Schröder et al. [54] in The Netherlands: on average the potentially leachable N in sandy soils is 73%, rising up to 90% in arable lands as reported also in Roelsma et al. [30].
Figure 3 shows the results of the WinEPIC simulation of the NO 3 ´concentration in soil water for two complete annual crop cycles, together with the values of measured NO 3 ´concentration in the lysimeters; during the second corn crop no samples were considered.
Based on a visual inspection, except for the initial values of the period, the WinEPIC model was able to successfully simulate the NO 3 ´concentration in the soil profile.The initial values were not well represented, probably because the model was still stabilizing.During the considered period, the observed NO 3 ´concentration was higher than 50 mg¨L ´1 in 48% of the measured data, and ranged from 2 to 302 mg¨L ´1 with an average of 66 mg¨L ´1.WinEPIC gives modeled values ranging from 1 to 342 mg¨L ´1 with an average of 53 mg¨L ´1, representing the high variability of measured data well.Higher values observed during winter were in conjunction with Italian ryegrass fertilization operations and/or high rainfall events, while those observed during summer were related to fertilizer application and irrigation on silage corn.In intensive agriculture and livestock areas, in particular where sandy soils are prevalent, excess of irrigation and fertilization may lead to pollution of groundwater.Similar findings in adjacent areas with intensive buffalo breeding were obtained by Infascelli et al. [55].Based on a visual inspection, except for the initial values of the period, the WinEPIC model was able to successfully simulate the NO3 − concentration in the soil profile.The initial values were not well represented, probably because the model was still stabilizing.During the considered period, the observed NO3 − concentration was higher than 50 mg L −1 in 48% of the measured data, and ranged from 2 to 302 mg L −1 with an average of 66 mg L −1 .WinEPIC gives modeled values ranging from 1 to 342 mg L −1 with an average of 53 mg L −1 , representing the high variability of measured data well.Higher values observed during winter were in conjunction with Italian ryegrass fertilization operations and/or high rainfall events, while those observed during summer were related to fertilizer application and irrigation on silage corn.In intensive agriculture and livestock areas, in particular where sandy soils are prevalent, excess of irrigation and fertilization may lead to pollution of groundwater.Similar findings in adjacent areas with intensive buffalo breeding were obtained by Infascelli et al. [55].
Obviously, evaluating the model performance taking into account only NO3 − concentrations is not sufficient for the assessment of its predictive power [29].Further simulation results such as the other N cycle components and the soil water cycle in relation to the N cycle should be considered [31].

Fertilizer Management Scenarios
From measured results we identified as critical aspects on NO3 − leaching the fertilizing practices adopted, and thus three different scenarios with different fertilizing strategies were simulated for the 2004-2015 period by WinEPIC, assuming current field conditions and management as the baseline (S0).In Table 4, average values (over the 12-yr period) for forage annual yield, total soluble N in the soil profile and N leaching rate, as given by WinEpic, are reported.Obviously, evaluating the model performance taking into account only NO 3 ´concentrations is not sufficient for the assessment of its predictive power [29].Further simulation results such as the other N cycle components and the soil water cycle in relation to the N cycle should be considered [31].

Fertilizer Management Scenarios
From measured results we identified as critical aspects on NO 3 ´leaching the fertilizing practices adopted, and thus three different scenarios with different fertilizing strategies were simulated for the 2004-2015 period by WinEPIC, assuming current field conditions and management as the baseline (S 0 ).In Table 4, average values (over the 12-yr period) for forage annual yield, total soluble N in the soil profile and N leaching rate, as given by WinEpic, are reported.The proposed scenarios allowed us to predict the effect of alternative N fertilization management in terms of its environmental (soluble N and N leaching rate) and agronomic (forage yield) performance.In the 12 yrs simulation, forage yield for the three alternative scenarios did not substantially differ from the baseline value: the predicted reduction, less than 5%, can be considered as negligible.We speculate this is likely due to the soil organic N buffer capacity.
In regards to the N cycle, the choice of the most sustainable scenario should be based on other considerations.Reducing both mineral and organic fertilization by 50% (S 1 ), a decrease of 40% in the average total soluble N was obtained, but only a decrease of 22% in the average annual N leaching rate was obtained.The NUE is increased up to 82% for both crops.Eliminating only the mineral fertilization on Italian ryegrass (S 2 ), the amount of total soluble N is reduced even by 60%, and the leaching rate is reduced by 33%.The NUE is 65% for the whole crop rotation, due to the increase of N efficiency in the Italian ryegrass cycle, while in silage corn the N efficiency was the same as the baseline.Finally, eliminating the mineral fertilization (S 3 ), the average total soluble N declined by 51% while the average annual N leaching rate was reduced by 80%.The NUE is higher than in the other scenarios (89%).However, in this last scenario, the silage corn plants might suffer an N deficit during growth stage.
Based on these considerations, the N mineral fertilization on Italian ryegrass could be eliminated, thus consistently reducing the N leaching rate and increasing the NUE with respect to the baseline, without considerable yield loss.This is consistent with the findings of Demurtas et al. [31] in similar climatic conditions, with the same crop rotation and similar management (cattle slurry and mineral fertilization).Zavattaro et al. [56] highlighted that in intensively corn-based cropping systems in northern Italy, the total supply of fertilizer could be reduced to minimize leaching and other potential losses.Morari and Giupponi [57] reported that reduced fertilization was highly effective in reducing the N impact on groundwater quality compared to high-input management.On the whole, considering the current situation and management of intensive agriculture and livestock areas, optimizing fertilization would result in a reduction of N leaching loss and represent a means of control for non-point source pollution.

Conclusions
In a farm located in Sabaudia, Latina Province plain, Central Italy, the monitoring results of NO 3 concentration in soil solution exceeded the maximum allowable NO 3 ´concentration (50 mg¨L ´1), following fertilizer application and irrigation or abundant rainfall.Data simulated by the WinEPIC model were reliable since the predicted data were consistent with the measured data.Thus, WinEPIC has the potential to predict the fate of N added to soil in relation to soluble N leaching loss below the soil profile using parameters derived from field observations.Three alternative scenarios with different fertilizing rates were developed, all of them showing a beneficial effect on N loss reduction and an improvement in the NUE, without any substantial decrease in both Italian ryegrass and silage corn biomass yields.Hence, the model can be used to predict the effect of fertilizing practices on N balance, soil N loss and crop biomass yield in NVZs.
The results of this paper suggest a possible improvement of current fertilizing practices of the area, which might be useful in formulating management strategies for intensive cropped catchment to reduce diffuse pollution from agriculture.
Further field-testing using data from other soils, crops and management is needed to extend the model's application to the whole catchment.

Figure 1 .
Figure 1.Location of the study area in Latina Province (Central Italy).The yellow area outlines the experimental field.

Figure 2 .
Figure 2. Monthly mean maximum and minimum temperature, and monthly total rainfall in Sabaudia area, Latina Province (Central Italy).

Figure 1 .
Figure 1.Location of the study area in Latina Province (Central Italy).The yellow area outlines the experimental field.

Figure 1 .
Figure 1.Location of the study area in Latina Province (Central Italy).The yellow area outlines the experimental field.

Figure 2 .
Figure 2. Monthly mean maximum and minimum temperature, and monthly total rainfall in Sabaudia area, Latina Province (Central Italy).

Figure 2 .
Figure 2. Monthly mean maximum and minimum temperature, and monthly total rainfall in Sabaudia area, Latina Province (Central Italy).

Figure 3
Figure 3 shows the results of the WinEPIC simulation of the NO3 − concentration in soil water for two complete annual crop cycles, together with the values of measured NO3 − concentration in the lysimeters; during the second corn crop no samples were considered.

Figure 3 .
Figure 3. Dynamics of the NO3 − concentration (mg L −1 ) in the soil solution measured and predicted by WinEPIC model during the Italian Ryegrass-silage corn rotation (2013-2015).

Figure 3 .
Figure 3. Dynamics of the NO 3 ´concentration (mg¨L ´1) in the soil solution measured and predicted by WinEPIC model during the Italian Ryegrass-silage corn rotation (2013-2015).

Table 1 .
Main soil physical and chemical characteristics at the water buffalo livestock farm in Sabaudia Municipality (Latina Province Plain).

Table 2 .
Agronomic management of Italian ryegrass (Lolim multiflorum L.) and silage corn (Zea mays L.) at the intensive buffalo livestock farm in Sabaudia.

Table 3 .
Mean annual nitrogen balance (kg¨ha ´1) of Italian ryegrass (Lolium multiflorum Lam.) and silage corn (Zea mays L.) rotation, calculated as the difference of input´output, and nitrogen use efficiency (NUE) evaluation.

Table 4 .
Predicted forage average annual yields, total soluble nitrogen and nitrogen leaching rate, in different scenarios.

Table 4 .
Predicted forage average annual yields, total soluble nitrogen and nitrogen leaching rate, in different scenarios.