1. Introduction
The Mediterranean basin is a recognized hot spot for climate change for the next decades [
1,
2], with modifications of rainfall amount and pattern and temperature increase, and where extreme events are expected to severely affect agricultural sector and food security [
3]. As reported by Duveiller et al. [
4], in Southern Europe future impacts of climate change on agriculture can be generalized by a decline in both productivity and suitability. Moreover, in Mediterranean regions, characterized by high interannual and seasonal rainfall variability (wet and cool periods from autumn to spring, and long dry periods in summer) nowadays, one of the most important issues caused by intensive agricultural farming systems is the reduction of soil organic carbon (SOC), with a possible worsening in the perspective of climate change, with major side effects on soil functioning [
5]. Several authors confirm that multiple forms of physical, chemical, and biological degradation affect the Mediterranean soils [
6,
7].
The EU Soil Strategy for 2030 provides the framework towards protecting and restoring soils, and ensuring that they are used sustainably, through a Proposal for a new Soil Health Law [
8]. The proposal provides a harmonized definition of soil health, puts in place a comprehensive and coherent monitoring framework, and fosters sustainable soil management and remediation of contaminated sites. In this context, organic and conservation agriculture should be possible solutions to achieve a more sustainable agricultural production, which includes ensuring and maintaining productive capacity for today and the future, and increasing productivity without harming the environment and natural resources. Sustainable agricultural practices such as crop rotation, cover crops, use of compost, and organic fertilizers can reduce the external inputs (e.g., pesticides, fertilizers, and herbicides) with the effect of increasing crop yield stability and biodiversity in the rhizosphere over time [
9,
10]. Such management practices help in maintaining the soil functions - food, feed, fiber, and fuel production; water regulation, purification, storage, transformation; carbon sequestration and climate regulation; habitat for functional and intrinsic biodiversity; nutrient cycling and provision - and water quality and quantity [
11]. All these functions should be maintained in an integrated, holistic way so that one function is not maximized at the cost of another [
12,
13].
While current climate change is certain and measurable, thanks to the direct observations and the long-term past data series comparison, future climate projections include uncertainties [
14]. Nevertheless, most of the studies over the Mediterranean basin indicated that the observed temperature and precipitation trends are expected to worsen in the next future [
15]. Future climate will exhibit an increased frequency of extreme events with maximum temperature exceeding 40°C which will represent the normal conditions in the next future [
16]. Indicators are useful tools for interpreting and summarizing the complexity of the impact of alternative scenarios of practices [
17]. Thus, appropriate prediction tools are required to characterize the vulnerability of agricultural systems in a future changing climate, and to implement the best practices.
Deterministic crop growth modelling has proven to be a major tool for analyzing the impacts of climate change on agricultural production. The Environmental Policy Integrated Climate (EPIC) agroecosystem model is extensively applied at field-scale and tested in many pedo-climatic conditions [
18]. It simulates crop production as a function of weather, soil conditions, and management practices [
19,
20]. EPIC model v.0810 [
21] was selected because it has been widely and successfully used for assessing the effects of management on crop productivity, soil water balance, and soil C and N dynamics in a range of environments and agricultural systems, including the United States [
22,
23], Argentina [
24], and Europe [
18,
25,
26]. The EPIC model was used also in a long-term organic vegetable field experiment, to evaluate the performance of agro-ecological practices as adaptation and mitigation measures to cope with climate change in Southern Italy [
27].
The present study aimed to evaluate the long-term (30-years) agro-environmental sustainability of a typical Mediterranean cropping system using a modelling approach under future climate change scenarios. The agricultural management practices were assessed by means of several indicators. Within the H2020 FATIMA project, measured data from the Italian field trail and the EPIC model were used to assess the long-term agro-environmental impacts and sustainability of two different nitrogen fertilization treatments on crop yields, Water Use Efficiency (WUE), Nitrogen Use Efficiency (NUE), SOC stock change, soil bulk density change, and soil N cycle (nitrate leaching and N2O emissions), under future climate change scenarios. The tested treatments, applied in a processing tomato/durum wheat rotation, were: i) conservative N fertilization methods, based on the adoption of compost, leguminous cover crops (fava bean), and poultry manure; and ii) synthetic N fertilization method, based on mineral N fertilizer - used as control.
2. Materials and Methods
2.1. Study Area
The study area is located in the Tarquinia coastal plain (Viterbo Province, Latium, Central Italy), 7 km NW of Tarquinia city, 2.7 km from seashore - 42° 69’ N and 11° 69’ E, at an average altitude of 25 m above sea level, with 3% mean slope (
Figure 1). The area is an intensive agricultural land, characterized by the cultivation of rainfed winter cereals, often in rotation with irrigated summer crops. The area lies within a Nitrate Vulnerable Zone (NVZ), where the excess in N input applied to the soil - coming from livestock and fertilizers - directly contributed to the groundwater pollution. The Nitrates Directive (Directive 91/676/EEC - Council of the European Communities Council, 1991) legally restricted annual N application to 170 kg ha
-1 as the maximum.
In autumn 2015, a field experiment was set-up within a 20 ha private farm, cultivated with a durum wheat (
Triticum durum Desf. var. Iride)/processing tomato (
Solanum lycopersicum L. var. Vulcano) rotation since 2005. The field experiment continued until June 2017. The soil of the experimental site, with a clay loam texture, was classified as Calcaric Cambic Phaeozems according to WRB [
28]. The experiment aimed to evaluate the sustainability of the durum wheat/processing tomato production system in the long-term. Two different field plots (A and B) were set up to test the effects of conservative N fertilization methods (CONS) – based on the use of compost, cover crops (for irrigated tomato) and poultry manure (for rainfed wheat) - in comparison with synthetic N management (SYN) on crop yield and some selected environmental quality indicators. In plot A, the rotation started with processing tomato as main crop, while in plot B the rotation started with durum wheat. To improve soil fertility and NUE, and to reduce the potential nitrate (NO
3-) leaching, the CONS method in tomato included the cultivation of fava bean (
Vicia faba L.), in the autumn-winter period as cover crop, incorporated as green manure before the tomato transplanting together with an organic amendment, i.e. compost derived by vegetal local agro-forestry residues. For durum wheat, the CONS method included the application of poultry manure as organic N fertilizer.
Table 1 and
Table 2 show an overview about nitrogen fertilization treatment for durum wheat and tomato, respectively.
2.2. Evaluation Procedure
The overall assessment of the long-term (30 years) agro-environmental sustainability and soil health of the two N fertilization treatments tested in the study area under future climate change scenarios was obtained in two steps. First, a set of agro-environmental indicators was considered, defining also some thresholds to obtain three classes of sustainability, low, medium and high (
Table 3). In the second step the predicted values of the agro-environmental indicators were assigned to the corresponding class and the long-term overall evaluation was assessed, considering climate change scenarios.
Where:
- –
Yield (Mg ha-1): grain for wheat (Mg ha dry matter weight) and fruits for tomato (Mg ha of fresh weight).
- –
Nitrogen Recovery Efficiency (NRE) (%): expressed as the partial factor productivity, is the ratio between crop yield and N applied with fertilizers [
30].
- –
WUE (Kg mm): ratio between the crop yield and the evapotranspiration (ET) during the growing season.
- –
NUE (kg kg
-1): expressed as the partial nutrient balance, is the simplest form of nutrient RE. It is calculated as the ratio between N content in grain or fruits and the N applied by fertilizers [
30]. 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 [
31].
- –
Cumulative nitrate loss by leaching (Kg ha): amount of NO3- lost from the soil (below the rooting depth) for the whole cropping season/year, expressed as the relative variation (%) to annual N input source.
- –
SOC stock change (kg ha): relative variation between final and initial values.
- –
Soil bulk density change (g cm3): relative variation between final and initial values.
The considered threshold values were based on the results of the FATIMA project – a 2-years field experiment - and on the historical experience and knowledge of the farmer involved in the experiment. For the NRE (%) and NUE (kg kg
-1) the threshold values were based on the study by Dobermann et al. [
30].
The NUE index value for the whole crop rotation indicates the percentage amount of N absorbed by crop yield. WUE is crucial for Mediterranean areas and represents an important indicator when evaluating long-term sustainability under climate change scenarios, particularly for rainfed crops - such as durum wheat - depending entirely on rainfall. The sustainability threshold for tomato lays within a very small range, due to irrigation and to the possibility for the farmer to tune the doses very finely according to the plant needs. For nitrate leaching the thresholds were not defined as absolute values, but as relative variation (%) in comparison with N inputs as fertilizer, both under SYN and under CONS treatments. The N given to the soil as NH4+ can be quickly converted to NO3− by nitrifier microorganisms. This anion is not absorbed by the soil, thus is easily released into the soil liquid phase, possibly moving to the leaching water flow. Sandy soils are particularly sensitive to nitrate leaching towards groundwater, due to their higher permeability. SOC is considered one of the most important indicators of soil quality and soil health, strictly linked to most of the ecosystem services provided by soil, such as nutrient and water cycles regulation, buffer capacity, biodiversity, and GHG emission regulation. Any reduction of SOC in Mediterranean conditions must be considered as negative and not desirable. Finally, bulk density change is an important soil quality indicator in terms of physical quality. Soil compaction due to a very intensive management could result in a reduction of water infiltration, an increase in the energy required to plough the soil, crusting and erosion.
2.3. EPIC Model
EPIC is a field-scale and a daily time step process-based model that has been developed to assess the impacts of soil management on biophysical and biogeochemical processes [
32], such as plant growth and development, soil water balance, C and nutrient cycling, soil erosion, and GHG emissions. The model, developed and maintained by researchers at Blackland Research and Extension Center, Texas A&M Agrilife Research (USA), was designed originally to explore the impacts of soil erosion on crop productivity [
19]. Afterwards, it was refined including additional sub-models to predict water quality and the response of crops to atmospheric CO
2 [
33]. EPIC and the derived models have been applied extensively to a variety of soils and cropping systems worldwide [
34]. As reported in Parton et al. [
35], EPIC has eight major components - modules on weather generation, crop growth, soil water dynamics, erosion, nutrient and carbon cycling, soil temperature, tillage, and soil-crop management - and operates on a continuous basis using a daily time step performing short- and long-term predictions. Simulated processes include the effects of tillage, fertilizer and irrigation on crop yield and soil agro-environmental quality (surface residue, soil bulk density, and biogeochemical cycles) in the considered crop rotation and cropping system.
Information about the cropping system management (such as tillage, irrigation volumes, fertilizers supply, and scheduling of operations), soil and weather data, and crop growth data, such as plants density and crop growing period, is mandatory to run the model. As reported by Folberth et al. [
36], in addition to plant growth and yield formation, EPIC estimates a wide range of environmental externalities, such as wind and water erosion rates, turnover and partitioning of soil organic carbon, N and P, evapotranspiration, fluxes of selected gases, and soil hydrological processes. Depending on N and lignin content, crop residues including roots are split in two litter compartments: metabolic and structural. From them, as a function of soil temperature and moisture, C is allocated in three compartments: microbial biomass, slow humus and passive humus, which are different in size, function and turnover times [
32]. Furthermore, the model accounts for the effects of change in CO
2 concentration and vapor pressure deficit on radiation-use efficiency, leaf resistance, and transpiration of crops to estimate the increase of plant growth and WUE [
37]. In this study, EPIC model v.0810 was used [
21].
Model calibration is the process of adjusting influential model parameters within their reasonable ranges to obtain realistic model results, consistent with the available observed data such as crop yields, soil nutrient content, soil carbon, soil water content, water infiltration rate and flow and water quality [
38]. In the study area the calibration procedure was performed using only the first year’s data.
Validation process consists in the assessment of the accuracy of the model predictions, comparing the results to additional and independent observed data. For the validation process we used the second-year field measures. Coefficient of determination (R2), slope and intercept of the linear regression, and correlation coefficient (r) between observed and simulated values were used to measure the model performance. The differences between model outputs from simulations at different steps were examined by analyzing the mean values, for better quantifying the effects of calibration on simulated crop yield and for understanding the uncertainties associated with the calibration procedures.
2.4. Future Climate Scenarios
In the present study, three different climate scenarios were used to run the EPIC model for long-term assessment. The climate scenarios were obtained from the MARS-AGRI4CAST website (
http://agri4cast.jrc.ec.europa.eu/DataPortal/Index.aspx?o=d), where present and future climate scenarios (two time projections, TPs) were generated by General Circulation Models (GCMs) from a consolidated daily weather dataset with a grid of 25x25 km. The GCM were: (1) METO-HC (METO); (2) DMI-HIRHAM5-ECHAM5 (ECHAM); and (3) ETHZ-CLM-HadCM3Q0 (ETHZ) and both present and the corresponding future climate for each of them were obtained. In terms of annual surface air temperature, the ECHAM future simulation is the coldest (15.7°C), while METO and ETHZ future simulation are the warmest (16.2°C). The temperature and rainfall in the GCMs for the two TPs are reported in
Table 4.
As regards annual rainfall, METO and ETHZ showed a similar precipitation pattern based on an increase in precipitation in comparison with the baseline. Conversely, the ECHAM future climate scenario showed a reduction in precipitation regime with respect to the baseline, and markedly different patterns than under the others. An annual increase of mean temperature compared to the corresponding baseline by 0.8, 0.6, and 1.2°C was predicted with METO, ECHAM and ETHZ, respectively. An increase in rainfall was observed with METO (52.2 mm, +17.4%) and ETHZ (14.4 mm, +4.2%), while a slight reduction in rainfall was predicted with ECHAM (-2.0 mm, 0.6%). Hence, each GCMs climate was run for two TPs, baseline and future climate, for 30 years. The TPs chosen were: (i) “2000” for the baseline, representing mean climate change for the period 1985–2015; and (ii) “2030” for climate change predictions, representing mean climate change for the period 2015-2044. Atmospheric CO2 concentrations for the considered periods were 400 ppm for baseline and 450 ppm for climate change. For all the baseline and climate change simulations, the predicted yield trend and the SOC stocks, mineral N, and bulk density changes were considered. Within the same experimental plot, all simulations were performed for each baseline and the corresponding future climate change scenario. In all the simulations the same values of soil parameters were considered as data input, in order to calculate the percentage of variation both for each baseline and future climate change scenario (soil parameter change = [(final value – initial value)/initial value]). Similarly, to compare the effects of the three future climate projections considered (METO, ECHAM, ETHZ), the relative variation of soil parameters between each climate change scenario and the corresponding baseline used as control was computed.
3. Results and Discussion
3.1. Agronomic Indicators
Figure 2 shows the changes in durum wheat and tomato average yields for the SYN and CONS treatment and the three GCMs in a climate change scenario with respect to the corresponding baseline. Under the METO climate change scenario, the average yield of durum wheat and processing tomato increased for each treatment except for tomato SYN, where a slight decrease was observed (-2%). The highest increase was obtained for durum wheat SYN (14%), while the CONS treatment for both crops showed an increase of 8%. This behavior is likely due to the positive effect of increased CO
2, combined with the increase of both rainfall and temperature observed under METO GCM. Under the ECHAM climate change scenario, which considers no significant changes in rainfall, the average yield of both crops increased under CONS treatment, while under SYN treatment only tomato yield increased, reaching a steady-state for durum-wheat. Finally, under the ETHZ climate change scenario, the highest increase was obtained for tomato CONS (+10%), followed by durum-wheat SYN (+3%).
Changes in NUE in the 30-yrs period are reported in
Figure 3. NUE of durum wheat and tomato crops benefit of climate change under ETHZ both for SYN and CONS treatments. Nevertheless, the highest NUE was obtained for durum wheat CONS (+23%) under ECHAM, while the lowest one was observed in METO for durum wheat under CONS treatment and processing tomato under SYN treatment (-3%).
Figure 4 shows the durum wheat and the tomato relative changes in WUE in the 30-yrs period for the SYN and CONS treatments, under the three GCMs in climate change scenarios with respect to the corresponding baseline. All the climate change scenarios showed higher performance than the relative baseline, except for ETHZ under tomato in SYN treatment.
In
Figure 5 the relative variations of crop yield and WUE for tomato SYN and CONS treatments under relative variations of temperature (
Figure 5a) and rainfall (
Figure 5b) in future climate change with respect to the baseline are reported. The effect of the foreseen increase in temperature on tomato yield is always positive (up to 10%) for CONS and is not significant for SYN (
Figure 5a). The change in rainfall, both positive (as in ETHZ) and negative (as in METO), is associated with an increase in yield in CONS, particularly with the intermediate rain increase. In SYN the yield is slightly influenced by the change in rain, and the lower increases in yield are associated with the higher rate of rainfall change. WUE is positively influenced by the most pronounced change in temperature and rainfall in CONS, while WUE is less affected by climate change in SYN. As a general consideration, yield and WUE in SYN are less influenced by climate change and are more stable. Yield and WUE in CONS are positively influenced by climate change, to a varying extent following the change in rainfall and temperature.
The higher stability of yield under climate change in SYN is confirmed also for durum wheat (
Figure 6).
3.2. Environmental Indicators
The effect of rainfall change on nitrate leaching in the tomato-wheat rotation is clearly different in SYN and CONS. In CONS both the increase in temperature and the decrease in rainfall reduced nitrate leaching, while the opposite can be observed for SYN (
Figure 7).
The bulk density in the durum wheat/tomato cropping system increased in the two simulated N management treatment (SYN and CONS) under current and future climate change scenarios, by about 0.02 g cm−3 on average, in the 30-yrs period.
Several studies have shown that SOC is affected by the management [
39,
40,
41,
42]. The SOC stock in the durum wheat/tomato cropping system decreased in the two simulated N management treatment (SYN and CONS) under current and future climate change scenarios, by about 0.62 Mg ha
−1 yr
−1, on average, in the 30-yrs period.
Figure 8 shows the relative SOC stock changes for the SYN and CONS treatments, under the three GCMs in climate change scenarios with respect to the corresponding baseline. Under the METO and ETHZ climate change scenarios, the relative SOC stock changes were always negative in both SYN and CONS treatments, ranging from -11% to -3%, on average. On the contrary, under ECHAM the SOC stock for the SYN and CONS treatments increased by 4%.
Changes in nitrate leaching in the 30-yrs period was reported in
Figure 9. Nitrate leaching decreased under CONS treatment in all the three-climate change scenario, while in the SYN treatment only in METO scenario a decrease with respect to the relative baseline is predicted.
The predicted cumulated N2O emissions were mainly influenced by N management treatment (SYN or CONS). The cumulated N2O emissions were null under SYN and were slightly positive under CONS (0.07 kg ha−1 on average), ranging from 0.004 to 0.009 kg ha−1. Future climate change scenarios only affected N2O emissions under CONS treatment. The higher emissions were observed under ECHAM climate change scenario, while similar values were detected under METO and ETHZ. Finally, the predicted cumulated N2O emissions were negligible or slightly positive. These low emissions can be attributed to the reduced N fertilization rates applied for durum wheat/processing tomato cropping systems, since the experimental field lies within a NVZ.
3.3. Overall Evaluation of Agro-Environmental Sustainability
In
Table 5 the results of the long-term simulations obtained by running the EPIC model for 30 consecutive years under three different climate change scenarios (METO, ECHAM and ETHZ) are reported. The table summarizes the overall agro-environmental sustainability assessment of the nitrogen management strategies adopted in the study area of the FATIMA project for durum wheat/processing tomato rotation. The overall evaluation compares agronomic and environmental indicators under SYN and CONS treatments. The agronomic indicators were evaluated for each treatment and crop, while the environmental indicators were grouped by cropping systems and treatments. In the table, the values of the indicators averaged on the three climate change scenarios were reported. For both agronomic and environmental indicators, a qualitative evaluation was performed according to three classes, and the comparison between CONS and SYN was computed as a relative variation. Each value is represented with the corresponding color of the evaluation class, thus fitting the data in the qualitative class they belong to.
In the context of future climate change scenario, characterized by higher annual surface air temperature and by differently distributed rainfall over the year -e.g. concentrated in few months, the trend of the agronomic and environmental indicators reported in the figure (i.e. reduction of some performances and improvement of some other aspects) as well as their performances must be read together and considered stable in the framework of the total agro-environmental sustainability assessment of cropping systems, in the mid- long-term (30-years on average).
In the long-term predictions, crop yield under future climate change showed a different trend for durum wheat and tomato in SYN and CONS treatments. As regards rainfed durum wheat, in SYN the yield was stable and remained within the medium class range (
Table 3), while a decrease in grain yield was observed in CONS treatment moving to low class (-34% in comparison with SYN). Aiming to increase durum wheat yield under future climate change scenarios in CONS treatment, some variations in farm management practices - such as supplemental irrigation and/or use of wheat varieties more resistant to drought - can be feasible solutions to address the problem. Conversely, in irrigated tomato both treatments maintained the yield in the best class range (+15% in CONS with respect to SYN). As regards RE, values observed for both crops and treatments are in the low-range class (below 50%). This indicates the strong relationship between crop yield and the capability of crops to use the N applied by fertilization. The lower NRE values observed under CONS treatment with respect to SYN (-47%) can be linked to the different patterns of N release in the two treatments (slow-release organic N under CONS, and faster release mineral N under SYN). In both cases, plant N uptake is linked to soil water availability in the critical crop growth stages. Considering the NUE, the higher value observed under SYN (>100) means that more N is removed with the harvested crop than applied by fertilizer. This situation is equivalent to “soil mining” of N, since soil N stock is used. On the other hand, the lower NUE observed under CONS is favorable, demonstrating a better efficiency of the crops in N uptake. In this case, organic N fertilizers become a positive factor, because they slowly release N for crops, contributing to increase soil organic matter into the soil. WUE of durum wheat (rainfed) and tomato (irrigated) crops were higher under CONS (high class range) than under SYN treatment (low class range). Therefore, incorporating compost and cover crop as green manure for tomato and using poultry manure for durum wheat showed a very positive effect in increasing the water retention of soil over time. Looking at the other indicators, the value of cumulative NO
3- loss by leaching in durum wheat/tomato cropping system is in the low class range in SYN and in the medium class range under CONS, showing a reduction of 38% in NO
3- percolation in CONS in comparison with SYN. Since the nitrates tend to accumulate in groundwater, these findings are particularly significant, especially considering the nitrate leaching at a wider geographical scale. The SOC stock change showed a greater decrease under SYN than under CONS. This is consistent with the behavior of the soil bulk density, where a lower value - favorable in terms of soil quality - was observed under CONS in comparison with SYN.
Several authors found that conservative agriculture had a positive impact on soil characteristics [
43,
44,
45]. Francaviglia et al. [
46] showed that longer crop rotations (3–5 years) and the introduction of legumes resulted in higher increases in SOC contents (18%) in Mediterranean sites. Williams et al. [
47] showed the impact of soil management on soil health. Organic matter inputs, such as on-farm compost, crop residue recycling, manure, or other organic fertilizers can improve soil fertility and SOC sequestration, under various climates and cropping systems [
48,
49]. In our case, since a greater SOC stock and a lower bulk density are desirable conditions for the objective of increasing the long-term environmental sustainability and the soil health, more conservative management strategies such as compost application (different types and rates), minimum or no-tillage, and agroecological service crops might be suggested. Anyway, it should be kept in mind that the conversion from one fertilization strategy to others considered more “environmental-friendly” - e.g., from the mineral fertilization to the organic one - not always is an assurance of higher sustainability [
50].
4. Conclusions
In the study area - which represents a typical Mediterranean cropping system - the agro-environmental sustainability of two different fertilization strategies was evaluated in the long-term using a modelling approach under future climate change scenarios. Soil fertility and crop productivity were affected by the management, since CONS treatment shows higher SOC stock and WUE compared to SYN. The rainfall influenced crop yield. The overall evaluation of the alternative fertilization strategy proposed is strongly dependent from both the environmental and the productive aspects, and should take into account the local applicability of the option and its profitability for the farmer. In terms of productivity, i.e. relative yield change in CONS in comparison with SYN, the effect is positive for tomato and negative for wheat in all the climate scenarios. Given the highest profitability of tomato compared to wheat, the proposed change is considered a feasible strategy under climate change scenario, and could be sustainable in the long-term. Considering the environmental indicators, SOC stock change and nitrate leaching, the effect of CONS in the foreseen climate change is strongly positive. Hence, despite some weakness of the strategy, i.e. type and rate of organic fertilizers and selection of cover crop, the proposed management represents a good option for the farmer, for maintaining the soil health and for protecting the environment.
Author Contributions
Conceptualization, R.N. and R.F.; methodology, R.N. and R.F.; validation, C.D.B. and C.P.; data curation, C.D.B., S.V., C.P. and R.F.; writing—original draft preparation, C.P.; writing—review and editing, S.V., R.N., C.D.B and R.F.; project administration, R.N.; funding acquisition, R.N. All the authors have read and agreed to the published version of the manuscript.
Funding
The research was conducted within the FATIMA (Farming Tools for external nutrient Inputs and Water Management) project (
http://fatima-h2020.eu/), funded by the European Union´s Horizon 2020 research and innovation programme (Grant Agreement no. 633945).
Data Availability Statement
Acknowledgments
The authors wish to thank Dr. Bruno Pennelli and Dr. Melania Migliore for their valuable help in field sampling and lab analyses. Special thanks are dedicated to the farmer Mr. Vincenzo Fava for making his land available for experimental set-up and data collection.
Conflicts of Interest
The authors declare no conflict of interest.
References
- IPCC. Climate Change 2013: The Physical Science Basis.; Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Stocker, T.F., Qin, D., Plattner, G.K., Tignor, M., Allen, S.K., Boschung, J., Nauels, A., Xia, Y., Bex, V., Midgley, P.M., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2013; 1535p. [Google Scholar]
- Cos, J.; Doblas-Reyes, F.; Jury, M.; Marcos, R.; Bretonnière, P.A.; Samsó, M. The Mediterranean climate change hotspot in the CMIP5 and CMIP6 projections. Earth Syst. Dynam. 2022, 13, 321–340. [Google Scholar] [CrossRef]
- Ventrella, D.; Charfeddine, M.; Moriondo, M.; Rinaldi, M.; Bindi, M. Agronomic adaptation strategies under climate change for winter durum wheat and tomato in Southern Italy: irrigation and nitrogen fertilization. Reg. Environ. Chang. 2012, 12, 407–419. [Google Scholar] [CrossRef]
- Duveiller, G.; Donatelli, M.; Fumagalli, D.; Zucchini, A.; Nelson, R.; Baruth, B. A dataset of future daily weather data for crop modelling over Europe derived from climate change scenarios. Theor. Appl. Climatol. 2017, 127, 573–585. [Google Scholar] [CrossRef]
- Álvaro-Fuentes, J.; Paustian, K. Potential soil carbon sequestration in a semiarid Mediterranean agroecosystem under climate change: quantifying management and climate effects. Plant Soil 2011, 338, 261–272. [Google Scholar] [CrossRef]
- Zdruli, P.; Jones, R.J.A.; Montanarella, L. Organic Matter in the Soils of Southern Europe. European Soil Bureau Technical Report, EUR 21083 EN, Office for Official Publications of the European Communities, Luxembourg, 2004; 16.
- Ferreira, C.S.S.; Seifollahi-Aghmiuni, S.; Destouni, G.; Ghajarnia, N.; Kalantari, Z. Soil degradation in the European Mediterranean region: Processes, status and consequences. Sci. Total Environ. 2022, 805, 150106. [Google Scholar] [CrossRef]
- COM 416 (2023), EU Directorate General for Environment, “Proposal for a Directive on Soil Monitoring and Resilience”. Available online: https://environment.ec.europa.eu/publications/proposal-directive-soil-monitoring-and-resilience_en (accessed on 28 September 2023).
- Farina, R.; Marchetti, A.; Francaviglia, R.; Napoli, R.; Di Bene, C. Modeling regional soil C stocks and CO2 emissions under Mediterranean cropping systems and soil types. Agric. Ecosyst. Environ. 2017, 238, 128–141. [Google Scholar] [CrossRef]
- Quintarelli, V.; Radicetti, E.; Allevato, E.; Stazi, S.R.; Haider, G.; Abideen, Z.; Bibi, S.; Jamal, A.; Mancinelli, R. Cover Crops for Sustainable Cropping Systems: A Review. Agriculture 2022, 12, 2076. [Google Scholar] [CrossRef]
- Coyle, C.; Creamer, R.E.; Schulte, R.P.O.; O’Sullivan, L.; Jordan, P. A functional land management conceptual framework under soil drainage and land use scenarios. Environ. Sci. Policy 2016, 56, 39–48. [Google Scholar] [CrossRef]
- Schulte, R.P.O.; Bampa, F.; Bardy, M.; Coyle, C.; Creamer, R.E.; Fealy, R.; Gardi, C.; Ghaley, B.B.; Jordan, P.; Laudon, H.; O'Donoghue, C.; Ó'hUallacháin, D.; O'Sullivan, L.; Rutgers, M.; Six, J.; Toth, G.L.; Vrebos, D. Making the most of our land: Managing soil functions from local to continental scale. Front. Environ. Sci. 2015, 3, 81. [Google Scholar] [CrossRef]
- Ghaley, B.B.; Rusu, T.; Sandén, T.; Spiegel, H.; Menta, C.; Visioli, G.; O’Sullivan, L.; Trinsoutrot Gattin, I.; Delgado, A.; Liebig, M.A.; Vrebos, D.; Szegi, T.; Michéli, E.; Cacovean, H.; Bugge Henriksen, C. Assessment of benefits of conservation agriculture on soil functions in arable production systems in Europe. Sustainability 2018, 10, 794. [Google Scholar] [CrossRef]
- Brilli, L.; Moriondo, M.; Ferrise, R.; Dibari, C.; Bindi, M. Climate change and Mediterranean crops: 2003 and 2012, two possible examples of the near future. Agrochimica 2014, 58, 1–14. [Google Scholar]
- IPCC (2007) Climate Change 2007: Synthesis report. Contribution of Working Group I, II and III to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. [Core Writing Team, Pachauri, R.K and Reisinger, A. (eds.)]. IPCC, Geneva, Switzerland, 2007; 104 pp.
- Battisti, D.S.; Naylor, R.L. Historical warnings of future food insecurity with unprecedented seasonal heat. Science 2009, 323, 240–244. [Google Scholar] [CrossRef] [PubMed]
- Bockstaller, C.; Guichard, L.; Makowski, D.; Aveline, A.; Girardin, P.; Plantureux, S. Agri-environmental indicators to assess cropping and farming systems. A review. Agron Sustain Dev 2008, 28, 139–149. [Google Scholar] [CrossRef]
- Farina, R.; Seddaiu, G.; Orsini, R.; Steglich, E.; Roggero, P.P.; Francaviglia, R. Soil carbon dynamics and crop productivity as influenced by climate change in a rainfed cereal system under contrasting tillage using EPIC. Soil Tillage Res. 2011, 112, 36–46. [Google Scholar] [CrossRef]
- Williams, J.R.; Jones, C.A.; Dyke, P.T. A modeling approach to determining the relationship between erosion and soil productivity. Trans. ASAE 1984, 27, 129–144. [Google Scholar] [CrossRef]
- Williams, J.R. The EPIC model. In: Singh, V.P. (Ed.), Computer Models of Watershed Hydrology. Water Res. Publ., Littleton, CO, 1995; pp. 909–1000.
- Gerik, T.; Williams, J.; Francis, L.; Greiner, J.; Magre, M.; Meinardus, A.; Steglich, E.; Taylor, R. Environmental Policy Integrated Climate Model - User’s Manual Version 0810, Blackland Research and Extension Center, Texas A&M AgriLife, Temple, USA. 2015. Available online: https://epicapex.tamu.edu/manuals-and-publications/ (accessed on 21 September 2023).
- Causarano, H.J.; Shaw, J.N.; Franzluebbers, A.J.; Reeves, D.W.; Raper, R.L.; Balkcom, K.S.; Norfleet, M.L.; Izaurralde, R.C. Simulating field-scale soil organic carbon dynamics using EPIC. Soil Sci. Soc. Am. J. 2007, 71, 1174–1185. [Google Scholar] [CrossRef]
- Causarano, H.J.; Doraiswamy, P.C.; McCarty, G.W.; Hatfield, J.L.; Milak, S.; Stern, A.J. EPIC modeling of soil organic carbon sequestration in croplands of Iowa. J. Environ. Qual. 2008, 37, 1345–1353. [Google Scholar] [CrossRef]
- Apezteguia, H.P.; Izaurralde, R.C.; Sereno, R. Simulation study of soil organic matter dynamics as affected by land use and agricultural practices in semiarid Córdoba, Argentina. Soil Tillage Res. 2009, 102, 101–108. [Google Scholar] [CrossRef]
- Billen, N.; Roder, C.; Gaiser, T.; Stahr, K. Carbon sequestration in soils of SW-Germany as affected by agricultural management-calibration of the EPIC model for regional simulations. Ecol. Model. 2009, 220, 71–80. [Google Scholar] [CrossRef]
- Rinaldi, M.; De Luca, D. Application of EPIC model to assess climate change impact on sorghum in southern Italy. Ital. J. Agron. 2012, 7, 74–85. [Google Scholar] [CrossRef]
- Di Bene, C.; Diacono, M.; Montemurro, F.; Testani, E.; Farina, R. EPIC model simulation to assess effective agro-ecological practices for climate change mitigation and adaptation in organic vegetable system. Agron Sustain Dev. 2022, 42, 7. [Google Scholar] [CrossRef]
- IUSS Working Group WRB. World reference base for soil resources 2014. Update 2015. International soil classification system for naming soils and creating legends for soil maps. In World Soil Resources Report No. 106; Food and Agriculture Organization of the United Nations: Rome, 2015; Available online: https://www.fao.org/3/i3794en/I3794en.pdf (accessed on 15 September 2023).
- Di Bene, C.; Vanino, S.; Nino, P.; Migliore, M.; Anzano, E.; Farina, R.; Pennelli, B.; Fabiani, S.; D’Urso, G.; Marchetti, A.; Piccini, C.; De Michele, C.; Canali, S.; Tittarelli, F.; Napoli, R. Coupling remote sensing and modeling approach for optimizing input management in a typical Mediterranean cropping system. In Proceedings of the XX AIAM Congress and XLVI SIA Congress, Milan, Italy, 12-14 September 2017. [Google Scholar] [CrossRef]
- Dobermann, A. Nutrient use efficiency - Measurement and management. In Fertilizer Best Management Practices; IFA International Workshop on Fertilizer Best Management Practices (FBMPs): Brussels, Belgium, 2017; pp. 1–28. [Google Scholar]
- Piccini, C.; Di Bene, C.; Farina, R.; Pennelli, B.; Napoli, R. Assessing Nitrogen Use Efficiency and Nitrogen Loss in a Forage-Based System Using a Modeling Approach. Agronomy 2016, 6, 23. [Google Scholar] [CrossRef]
- Izaurralde, R.C.; Williams, J.R.; McGill, W.B.; Rosenberg, N.J.; Jakas, M.C.Q. Simulating soil C dynamics with EPIC. Model description and testing against long-term data. Ecol. Model. 2006, 192, 362–384. [Google Scholar] [CrossRef]
- Gassman, P.W.; Williams, J.R.; Benson, V.W.; Izaurralde, R.C.; Hauck, L.; Jones, C.A.; Atwood, J.D.; Kiniry, J.; Flowers, J.D. Historical development and applications of the EPIC and APEX models. In Working paper 05-WP 397; Center for Agricultural and Rural Development, Iowa State University: Ames, 2005. [Google Scholar]
- He, X.; Izaurralde, R.C.; Vanotti, M.B.; Williams, J.R.; Thomson, A.M. Simulating long-term and residual effects of nitrogen fertilization on corn yields, soil carbon sequestration, and soil nitrogen dynamics. J. Environ. Qual. 2006, 35, 1608–1619. [Google Scholar] [CrossRef] [PubMed]
- Parton, W.J.; Schime, D.S.; Ojima, D.S.; Cole, C.V. A general model for soil organic matter dynamics: sensitivity to litter chemistry, texture and management. In Quantitative Modeling of Soil Forming Processes; Bryant, R.B., Arnold, R.W., Eds.; SSSA Special Publication: SSSA, Madison, WI; Volume 39, pp. 147–167.
- Folberth, C.; Elliott, J.; Müller, C.; Balkovic, J.; Chryssanthacopoulos, J.; Izaurralde, R.C.; Jones, C.D.; Khabarov, N.; Liu, W.; Reddy, A.; Schmid, E.; Skalský, R.; Yang, H.; Arneth, A.; Ciais, P.; Deryng, D.; Lawrence, P.J.; Olin, S.; Pugh, T.A.M.; Ruane, A.C.; Wang, X. Uncertainties in global crop model frameworks: effects of cultivar distribution, crop management and soil handling on crop yield estimates. Biogeosciences Discuss 2016. [Google Scholar] [CrossRef]
- Stockle, C.O.; Williams, J.R.; Jones, C.A.; Rosenberg, N.J. A method for estimating the direct and climatic effects of rising atmospheric carbon dioxide on growth and yield of crops. I. Modification of the EPIC model for climate change analysis. Agric. Syst. 1992, 38, 225–238. [Google Scholar] [CrossRef]
- Kersebaum, K.C.; Boote, K.J.; Jorgenson, J.S.; Nendel, C.; Bindi, M.; Frühauf, C.; Gaiser, T.; Hoogenboom, G.; Kollas, C.; Olesen, J.E.; Rotter, R.P.; Ruget, F.; Thorburn, P.J.; Trnka, M.; Wegehenkel, M. Analysis and classification of data sets for calibration and validation of agro-ecosystem models. Environ. Model. Softw. 2015, 72, 402–417. [Google Scholar] [CrossRef]
- Spargo, J.T.; Alley, M.M.; Follett, R.F.; Wallace, J.V. Soil carbon sequestration with continuous no-till management of grain cropping systems in the Virginia coastal plain. Soil Tillage Res 2008, 100, 133–140. [Google Scholar] [CrossRef]
- Alijani, K.; Bahrani, M.J.; Kazemeini, S.A. Short-term responses of soil and wheat yield to tillage, corn residue management, and nitrogen fertilization. Soil Tillage Res 2012, 124, 78–82. [Google Scholar] [CrossRef]
- Khorami, S.S.; Kazemeini, S.A.; Afzalinia, S.; Gathala, M.K. Changes in soil properties and productivity under different tillage practices and wheat genotypes: A short-term study in Iran. Sustainability 2018, 10, 3273. [Google Scholar] [CrossRef]
- Liu, Z.; Cao, S.; Sun, Z.; Wang, H.; Qu, S.; Lei, N.; He, J.; Dong, Q. Tillage effects on soil properties and crop yield after land reclamation. Sci. Rep. 2021, 11, 4611. [Google Scholar] [CrossRef] [PubMed]
- Troccoli, A.; Maddaluno, C.; Mucci, M.; Russo, M.; Rinaldi, M. Is it appropriate to support the farmers for adopting conservation agriculture? Economic and environmental impact assessment. Ital. J. Agron. 2015, 10, 169–177. [Google Scholar] [CrossRef]
- Bai, Z.; Caspari, T.; Ruiperez-Gonzalez, M.; Batjes, N.; Mäder, P.; Bünemann, E.K. Critical review of soil quality indicators. iSQAPER Deliverable 3.2, 2016; Ref. Ares(2018)14240.
- Vanino, S.; Di Bene, C.; Piccini, C.; Fila, G.; Pennelli, B.; Zornoza, R.; Sanchez-Navarro, V.; Álvaro-Fuentes, J.; Hüppi, R.; Six, J.; Farina, R. A comprehensive assessment of diversified cropping systems on agro-environmental sustainability in three Mediterranean long-term field experiments. Eur J Agron 2022, 140, 126598. [Google Scholar] [CrossRef]
- Francaviglia, R.; Álvaro-Fuentes, J.; Di Bene, C.; Gai, L.; Regina, K.; Turtola, E. Diversified arable cropping systems and management schemes in selected European regions have positive effects on soil organic carbon content. Agriculture 2019, 9, 261. [Google Scholar] [CrossRef]
- Williams, H.; Colombi, T.; Keller, T. The influence of soil management on soil health: An on-farm study in southern Sweden. Geoderma 2020, 360, 114010. [Google Scholar] [CrossRef]
- Arunrat, N.; Pumijumnong, N.; Hatano, R. Predicting local-scale impact of climate change on rice yield and soil organic carbon sequestration: a case study in Roi Et Province, Northeast Thailand. Agric Syst 2018, 164, 58–70. [Google Scholar] [CrossRef]
- Farina, R.; Testani, E.; Campanelli, G.; Leteo, F.; Napoli, R.; Canali, S.; Tittarelli, F. Potential carbon sequestration in a Mediterranean organic vegetable cropping system. A model approach for evaluating the effects of compost and Agro-ecological Service Crops (ASCs). Agric Syst 2018, 162, 239–248. [Google Scholar] [CrossRef]
- Fabiani, S.; Vanino, S.; Napoli, R.; Nino, P. Water energy food nexus approach for sustainability assessment at farm level: An experience from an intensive agricultural area in central Italy. Environ. Sci. Policy 2020, 104, 1–12. [Google Scholar] [CrossRef]
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).