Simulating improved combinations tillage-rotation under dryland conditions.

Crop simulation models allow analyzing various tillage-rotation combinations and exploring management scenarios. This study was conducted to test the DSSAT (Decision Support System for Agrotechnology Transfer) modelling system in rainfed semiarid central Spain. The focus is on the combined effect of tillage system and winter cereal-based rotations (cereal/legume/fallow) on the crop yield and soil quality. The observed data come from a 16-year field experiment. The CERES and CROPGRO models, included in DSSAT v4.5, were used to simulate crop growth and yield, and DSSAT- CENTURY was used in the soil organic carbon (SOC) and soil nitrogen (SN) simulations. Genetic coefficients were calibrated using part of the observed data. Field observations showed that barley grain yield was lower for continuous cereal (BB) than for vetch (VB) and fallow (FB) rotations for both tillage systems. The CERES-Barley model also reflected this trend. The model predicted higher yield in the conventional tillage (CT) than in the no tillage (NT) probably due to the higher nitrogen availability in the CT, shown in the simulations. The SOC and SN in the top layer only, were higher in NT than in CT, and decreased with depth in both simulated and observed values. These results suggest that CT-VB and CT-FB were the best combinations for the dry land conditions studied. However, CT presented lower SN and SOC content than NT. This study shows how models can be a useful tool for assessing and predicting crop growth and yield, under different management systems and under specific edapho-climatic conditions. Additional key words: CENTURY model; CERES-Barley; crop simulation models; DSSAT; sequential simula- tion; soil organic carbon.

The role of cover crops in irrigated systems: water balance, nitrate leaching and soil mineral nitrogen accumulation

Soil salinity and salt leaching are a risk for sustainable agricultural production in many irrigated areas. This study was conducted over 3.5 years to determine how replacing the usual winter fallow with a cover crop (CC) affects soil salt accumulation and salt leaching in irrigated systems. Treatments studied during the period between summer crops were: barley (Hordeum vulgare L.), vetch (Vicia villosa L.) and fallow. Soil water content was monitored daily to a depth of 1.3 m and used with the numerical model WAVE to calculate drainage. Electrical conductivity (EC) was measured in soil solutions periodically, and in the soil saturated paste extracts before sowing CC and maize. Salt leaching was calculated multiplying drainage by total dissolved salts in the soil solution, and use to obtain a salt balance. Total salt leaching over the four winter fallow periods was 26 Mg ha−1, whereas less than 18 Mg ha−1 in the presence of a CC. Periods of salt gain occurred more often in the CC than in the fallow. By the end of the experiment, net salt losses occurred in all treatments, owing to occasional periods of heavy rainfall. The CC were more prone than the fallow to reduce soil salt accumulation during the early growth stages of the subsequent cash crop.

Nitrogen deposition alters nitrogen cycling and reduces soil carbon content in low-productivity semiarid Mediterranean ecosystems.

Anthropogenic N deposition poses a threat to European Mediterranean ecosystems. We combined data from an extant N deposition gradient (4.3–7.3 kg N ha−1 yr−1) from semiarid areas of Spain and a field experiment in central Spain to evaluate N deposition effects on soil fertility, function and cyanobacteria community. Soil organic N did not increase along the extant gradient. Nitrogen fixation decreased along existing and experimental N deposition gradients, a result possibly related to compositional shifts in soil cyanobacteria community. Net ammonification and nitrification (which dominated N-mineralization) were reduced and increased, respectively, by N fertilization, suggesting alterations in the N cycle. Soil organic C content, C:N ratios and the activity of β-glucosidase decreased along the extant gradient in most locations. Our results suggest that semiarid soils in low-productivity sites are unable to store additional N inputs, and that are also unable to mitigate increasing C emissions when experiencing increased N deposition.

Seasonal evolution of soil respiration in a mixed forestry plantation of walnuts (Juglans × intermedia Carr.) and alders (Alnus cordata (Loisel.) Duby) in Domaine de Restinclières (France)

The ecological intensification of crops is proposed as a solution to the growing demand of agricultural and forest resources, in opposition to intensive monocultures. The introduction of mixed cultures as mixtures between nitrogen fixing species and non nitrogen fixing species intended to increase crop yield as a result of an improvement of the available nitrogen and phosphorus in soil. Relationship between crops have received little attention despite the wide range of advantages that confers species diversity to these systems, such as increased productivity, resilience to disruption and ecological sustainability. Forests and forestry plantations can develop an important role in storing carbon in their tissues, especially in wood which become into durable product. A simplifying parameter to analyze the amount allocated carbon by plantation is the TBCA (total belowground carbon allocation), whereby, for short periods and mature plantations, is admitted as the subtraction between soil carbon efflux and litterfall. Soil respiration depends on a wide range of factors, such as soil temperature and soil water content, soil fertility, presence and type of vegetation, among others. The studied orchard is a mixed forestry plantation of hybrid walnuts(Juglans × intermedia Carr.) for wood and alders (Alnus cordata (Loisel.) Duby.), a nitrogen fixing specie through the actinomycete Frankia alni ((Woronin, 1866) Von Tubeuf 1895). The study area is sited at Restinclières, a green area near Montpellier (South of France). In the present work, soil respiration varied greatly throughout the year, mainly influenced by soil temperature. Soil water content did not significantly influence the response of soil respiration as it was constant during the measurement period and under no water stress conditions. Distance between nearest walnut and measurement was also a highly influential factor in soil respiration. Generally there was a decreasing trend in soil respiration when the distance to the nearest tree increased. It was also analyzed the response of soil respiration according to alder presence and fertilizer management (50 kg N·ha-1·año-1 from 1999 to 2010). None of these treatments significantly influenced soil respiration, although previous studies noticed an inhibition in rates of soil respiration under fertilized conditions and high rates of available nitrogen. However, treatments without fertilization and without alder presence obtained higher respiration rates in those cases with significant differences. The lack of significant differences between treatments may be due to the high coefficient of variation experienced by soil respiration measurements. Finally an asynchronous fluctuation was observed between soil respiration and litterfall during senescence period. This is possibly due to the slowdown in the emission of exudates by roots during senescence period, which are largely related to microbial activity.

Nitrogen uptake dynamics, yield and quality as influenced by nitrogen fertilization in "Piel de sapo" melon.

The need to reduce nitrogen (N) fertilizer pollution strengthens the importance of improving the utilization efficiency of applied N to crops. This requires knowledge of crop N uptake characteristics and how fertilization management affects it. A three-year field experiment was conducted from May to September in central Spain to investigate the influence of different N rates, which ranged from 11 to 393 kg ha-1, applied through drip irrigation, on the dynamics of N uptake, nitrogen use efficiency (NUE), fruit yield and quality of a ?Piel de sapo? melon crop (Cucumis melo L. cv. Sancho). Both N concentration and N content increased in different plant parts with the N rate. Leaves had the highest N concentration, which declined by 40-50% from 34-41 days after transplanting (DAT), while the highest N uptake rate was observed from 30-35 to 70-80 DAT, coinciding with fruit development. In each year, NUE declined with increasing N rate. With N fertilizer applications close to the optimum N rate of 90-100 kg ha-1, the fruits removed approximately 60 kg N ha-1, and the amount of N in the crop residue was about 80 kg N ha-1; this serves to replenish the organic nutrient pool in the soil and may be used by subsequent crops following mineralization.

Role of maize stover incorporation on nitrogen oxide emissions in a non-irrigated Mediterranean barley field.

Aims Agricultural soils in semiarid Mediterranean areas are characterized by low organic matter contents and low fertility levels. Application of crop residues and/or manures as amendments is a cost-effective and sustainable alternative to overcome this problem. However, these management practices may induce important changes in the nitrogen oxide emissions from these agroecosystems, with additional impacts on carbon dioxide emissions. In this context, a field experiment was carried out with a barley (Hordeum vulgare L.) crop under Mediterranean conditions to evaluate the effect of combining maize (Zea mays L.) residues and N fertilizer inputs (organic and/or mineral) on these emissions. Methods Crop yield and N uptake, soil mineral N concentrations, dissolved organic carbon (DOC), denitrification capacity, N2O, NO and CO2 fluxes were measured during the growing season. Results The incorporation of maize stover increased N2O emissions during the experimental period by c. 105 %. Conversely, NO emissions were significantly reduced in the plots amended with crop residues. The partial substitution of urea by pig slurry reduced net N2O emissions by 46 and 39 %, with and without the incorporation of crop residues respectively. Net emissions of NO were reduced 38 and 17 % for the same treatments. Molar DOC:NO 3 − ratio was found to be a robust predictor of N2O and NO fluxes. Conclusions The main effect of the interaction between crop residue and N fertilizer application occurred in the medium term (4–6 month after application), enhancing N2O emissions and decreasing NO emissions as consequence of residue incorporation. The substitution of urea by pig slurry can be considered a good management strategy since N2O and NO emissions were reduced by the use of the organic residue.

Evaluation of cover crops based on characteristics of nitrogen use, and aerial and root growth

La caracterización de los cultivos cubierta (cover crops) puede permitir comparar la idoneidad de diferentes especies para proporcionar servicios ecológicos como el control de la erosión, el reciclado de nutrientes o la producción de forrajes. En este trabajo se estudiaron bajo condiciones de campo diferentes técnicas para caracterizar el dosel vegetal con objeto de establecer una metodología para medir y comparar las arquitecturas de los cultivos cubierta más comunes. Se estableció un ensayo de campo en Madrid (España central) para determinar la relación entre el índice de área foliar (LAI) y la cobertura del suelo (GC) para un cultivo de gramínea, uno de leguminosa y uno de crucífera. Para ello se sembraron doce parcelas con cebada (Hordeum vulgare L.), veza (Vicia sativa L.), y colza (Brassica napus L.). En 10 fechas de muestreo se midieron el LAI (con estimaciones directas y del LAI-2000), la fracción interceptada de la radiación fotosintéticamente activa (FIPAR) y la GC. Un experimento de campo de dos años (Octubre-Abril) se estableció en la misma localización para evaluar diferentes especies (Hordeum vulgare L., Secale cereale L., x Triticosecale Whim, Sinapis alba L., Vicia sativa L.) y cultivares (20) en relación con su idoneidad para ser usadas como cultivos cubierta. La GC se monitorizó mediante análisis de imágenes digitales con 21 y 22 muestreos, y la biomasa se midió 8 y 10 veces, respectivamente para cada año. Un modelo de Gompertz caracterizó la cobertura del suelo hasta el decaimiento observado tras las heladas, mientras que la biomasa se ajustó a ecuaciones de Gompertz, logísticas y lineales-exponenciales. Al final del experimento se determinaron el C, el N y el contenido en fibra (neutrodetergente, ácidodetergente y lignina), así como el N fijado por las leguminosas. Se aplicó el análisis de decisión multicriterio (MCDA) con objeto de obtener un ranking de especies y cultivares de acuerdo con su idoneidad para actuar como cultivos cubierta en cuatro modalidades diferentes: cultivo de cobertura, cultivo captura, abono verde y forraje. Las asociaciones de cultivos leguminosas con no leguminosas pueden afectar al crecimiento radicular y a la absorción de N de ambos componentes de la mezcla. El conocimiento de cómo los sistemas radiculares específicos afectan al crecimiento individual de las especies es útil para entender las interacciones en las asociaciones, así como para planificar estrategias de cultivos cubierta. En un tercer ensayo se combinaron estudios en rhizotrones con extracción de raíces e identificación de especies por microscopía, así como con estudios de crecimiento, absorción de N y 15N en capas profundas del suelo. Las interacciones entre raíces en su crecimiento y en el aprovisionamiento de N se estudiaron para dos de los cultivares mejor valorados en el estudio previo: uno de cebada (Hordeum vulgare L. cv. Hispanic) y otro de veza (Vicia sativa L. cv. Aitana). Se añadió N en dosis de 0 (N0), 50 (N1) y 150 (N2) kg N ha-1. Como resultados del primer estudio, se ajustaron correctamente modelos lineales y cuadráticos a la relación entre la GC y el LAI para todos los cultivos, pero en la gramínea alcanzaron una meseta para un LAI>4. Antes de alcanzar la cobertura total, la pendiente de la relación lineal entre ambas variables se situó en un rango entre 0.025 y 0.030. Las lecturas del LAI-2000 estuvieron correlacionadas linealmente con el LAI, aunque con tendencia a la sobreestimación. Las correcciones basadas en el efecto de aglutinación redujeron el error cuadrático medio del LAI estimado por el LAI-2000 desde 1.2 hasta 0.5 para la crucífera y la leguminosa, no siendo efectivas para la cebada. Esto determinó que para los siguientes estudios se midieran únicamente la GC y la biomasa. En el segundo experimento, las gramíneas alcanzaron la mayor cobertura del suelo (83-99%) y la mayor biomasa (1226-1928 g m-2) al final del mismo. Con la mayor relación C/N (27-39) y contenido en fibra digestible (53-60%) y la menor calidad de residuo (~68%). La mostaza presentó elevadas GC, biomasa y absorción de N en el año más templado en similitud con las gramíneas, aunque escasa calidad como forraje en ambos años. La veza presentó la menor absorción de N (2.4-0.7 g N m-2) debido a la fijación de N (9.8-1.6 g N m-2) y escasa acumulación de N. El tiempo térmico hasta alcanzar el 30% de GC constituyó un buen indicador de especies de rápida cubrición. La cuantificación de las variables permitió hallar variabilidad entre las especies y proporcionó información para posteriores decisiones sobre la selección y manejo de los cultivos cubierta. La agregación de dichas variables a través de funciones de utilidad permitió confeccionar rankings de especies y cultivares para cada uso. Las gramíneas fueron las más indicadas para los usos de cultivo de cobertura, cultivo captura y forraje, mientras que las vezas fueron las mejor como abono verde. La mostaza alcanzó altos valores como cultivo de cobertura y captura en el primer año, pero el segundo decayó debido a su pobre actuación en los inviernos fríos. Hispanic fue el mejor cultivar de cebada como cultivo de cobertura y captura, mientras que Albacete como forraje. El triticale Titania alcanzó la posición más alta como cultiva de cobertura, captura y forraje. Las vezas Aitana y BGE014897 mostraron buenas aptitudes como abono verde y cultivo captura. El MCDA permitió la comparación entre especies y cultivares proporcionando información relevante para la selección y manejo de cultivos cubierta. En el estudio en rhizotrones tanto la mezcla de especies como la cebada alcanzaron mayor intensidad de raíces (RI) y profundidad (RD) que la veza, con valores alrededor de 150 cruces m-1 y 1.4 m respectivamente, comparados con 50 cruces m-1 y 0.9 m para la veza. En las capas más profundas del suelo, la asociación de cultivos mostró valores de RI ligeramente mayores que la cebada en monocultivo. La cebada y la asociación obtuvieron mayores valores de densidad de raíces (RLD) (200-600 m m-3) que la veza (25-130) entre 0.8 y 1.2 m de profundidad. Los niveles de N no mostraron efectos claros en RI, RD ó RLD, sin embargo, el incremento de N favoreció la proliferación de raíces de veza en la asociación en capas profundas del suelo, con un ratio cebada/veza situado entre 25 a N0 y 5 a N2. La absorción de N de la cebada se incrementó en la asociación a expensas de la veza (de ~100 a 200 mg planta-1). Las raíces de cebada en la asociación absorbieron también más nitrógeno marcado de las capas profundas del suelo (0.6 mg 15N planta-1) que en el monocultivo (0.3 mg 15N planta-1). ABSTRACT Cover crop characterization may allow comparing the suitability of different species to provide ecological services such as erosion control, nutrient recycling or fodder production. Different techniques to characterize plant canopy were studied under field conditions in order to establish a methodology for measuring and comparing cover crops canopies. A field trial was established in Madrid (central Spain) to determine the relationship between leaf area index (LAI) and ground cover (GC) in a grass, a legume and a crucifer crop. Twelve plots were sown with either barley (Hordeum vulgare L.), vetch (Vicia sativa L.), or rape (Brassica napus L.). On 10 sampling dates the LAI (both direct and LAI-2000 estimations), fraction intercepted of photosynthetically active radiation (FIPAR) and GC were measured. A two-year field experiment (October-April) was established in the same location to evaluate different species (Hordeum vulgare L., Secale cereale L., x Triticosecale Whim, Sinapis alba L., Vicia sativa L.) and cultivars (20) according to their suitability to be used as cover crops. GC was monitored through digital image analysis with 21 and 22 samples, and biomass measured 8 and 10 times, respectively for each season. A Gompertz model characterized ground cover until the decay observed after frosts, while biomass was fitted to Gompertz, logistic and linear-exponential equations. At the end of the experiment C, N, and fiber (neutral detergent, acid and lignin) contents, and the N fixed by the legumes were determined. Multicriteria decision analysis (MCDA) was applied in order to rank the species and cultivars according to their suitability to perform as cover crops in four different modalities: cover crop, catch crop, green manure and fodder. Intercropping legumes and non-legumes may affect the root growth and N uptake of both components in the mixture. The knowledge of how specific root systems affect the growth of the individual species is useful for understanding the interactions in intercrops as well as for planning cover cropping strategies. In a third trial rhizotron studies were combined with root extraction and species identification by microscopy and with studies of growth, N uptake and 15N uptake from deeper soil layers. The root interactions of root growth and N foraging were studied for two of the best ranked cultivars in the previous study: a barley (Hordeum vulgare L. cv. Hispanic) and a vetch (Vicia sativa L. cv. Aitana). N was added at 0 (N0), 50 (N1) and 150 (N2) kg N ha-1. As a result, linear and quadratic models fitted to the relationship between the GC and LAI for all of the crops, but they reached a plateau in the grass when the LAI > 4. Before reaching full cover, the slope of the linear relationship between both variables was within the range of 0.025 to 0.030. The LAI-2000 readings were linearly correlated with the LAI but they tended to overestimation. Corrections based on the clumping effect reduced the root mean square error of the estimated LAI from the LAI-2000 readings from 1.2 to less than 0.50 for the crucifer and the legume, but were not effective for barley. This determined that in the following studies only the GC and biomass were measured. In the second experiment, the grasses reached the highest ground cover (83- 99%) and biomass (1226-1928 g/m2) at the end of the experiment. The grasses had the highest C/N ratio (27-39) and dietary fiber (53-60%) and the lowest residue quality (~68%). The mustard presented high GC, biomass and N uptake in the warmer year with similarity to grasses, but low fodder capability in both years. The vetch presented the lowest N uptake (2.4-0.7 g N/m2) due to N fixation (9.8-1.6 g N/m2) and low biomass accumulation. The thermal time until reaching 30% ground cover was a good indicator of early coverage species. Variable quantification allowed finding variability among the species and provided information for further decisions involving cover crops selection and management. Aggregation of these variables through utility functions allowed ranking species and cultivars for each usage. Grasses were the most suitable for the cover crop, catch crop and fodder uses, while the vetches were the best as green manures. The mustard attained high ranks as cover and catch crop the first season, but the second decayed due to low performance in cold winters. Hispanic was the most suitable barley cultivar as cover and catch crop, and Albacete as fodder. The triticale Titania attained the highest rank as cover and catch crop and fodder. Vetches Aitana and BGE014897 showed good aptitudes as green manures and catch crops. MCDA allowed comparison among species and cultivars and might provide relevant information for cover crops selection and management. In the rhizotron study the intercrop and the barley attained slightly higher root intensity (RI) and root depth (RD) than the vetch, with values around 150 crosses m-1 and 1.4 m respectively, compared to 50 crosses m-1 and 0.9 m for the vetch. At deep soil layers, intercropping showed slightly larger RI values compared to the sole cropped barley. The barley and the intercropping had larger root length density (RLD) values (200-600 m m-3) than the vetch (25-130) at 0.8-1.2 m depth. The topsoil N supply did not show a clear effect on the RI, RD or RLD; however increasing topsoil N favored the proliferation of vetch roots in the intercropping at deep soil layers, with the barley/vetch root ratio ranging from 25 at N0 to 5 at N2. The N uptake of the barley was enhanced in the intercropping at the expense of the vetch (from ~100 mg plant-1 to 200). The intercropped barley roots took up more labeled nitrogen (0.6 mg 15N plant-1) than the sole-cropped barley roots (0.3 mg 15N plant-1) from deep layers.

Improving the management of water and nitrogen resources in cropping systems by using simulation models

Los modelos de simulación de cultivos permiten analizar varias combinaciones de laboreo-rotación y explorar escenarios de manejo. El modelo DSSAT fue evaluado bajo condiciones de secano en un experimento de campo de 16 años en la semiárida España central. Se evaluó el efecto del sistema de laboreo y las rotaciones basadas en cereales de invierno, en el rendimiento del cultivo y la calidad del suelo. Los modelos CERES y CROPGRO se utilizaron para simular el crecimiento y rendimiento del cultivo, mientras que el modelo DSSAT CENTURY se utilizó en las simulaciones de SOC y SN. Tanto las observaciones de campo como las simulaciones con CERES-Barley, mostraron que el rendimiento en grano de la cebada era mas bajo para el cereal continuo (BB) que para las rotaciones de veza (VB) y barbecho (FB) en ambos sistemas de laboreo. El modelo predijo más nitrógeno disponible en el laboreo convencional (CT) que en el no laboreo (NT) conduciendo a un mayor rendimiento en el CT. El SOC y el SN en la capa superficial del suelo, fueron mayores en NT que en CT, y disminuyeron con la profundidad en los valores tanto observados como simulados. Las mejores combinaciones para las condiciones de secano estudiadas fueron CT-VB y CT-FB, pero CT presentó menor contenido en SN y SOC que NT. El efecto beneficioso del NT en SOC y SN bajo condiciones Mediterráneas semiáridas puede ser identificado por observaciones de campo y por simulaciones de modelos de cultivos. La simulación del balance de agua en sistemas de cultivo es una herramienta útil para estudiar como el agua puede ser utilizado eficientemente. La comparación del balance de agua de DSSAT , con una simple aproximación “tipping bucket”, con el modelo WAVE más mecanicista, el cual integra la ecuación de Richard , es un potente método para valorar el funcionamiento del modelo. Los parámetros de suelo fueron calibrados usando el método de optimización global Simulated Annealing (SA). Un lisímetro continuo de pesada en suelo desnudo suministró los valores observados de drenaje y evapotranspiración (ET) mientras que el contenido de agua en el suelo (SW) fue suministrado por sensores de capacitancia. Ambos modelos funcionaron bien después de la optimización de los parámetros de suelo con SA, simulando el balance de agua en el suelo para el período de calibración. Para el período de validación, los modelos optimizados predijeron bien el contenido de agua en el suelo y la evaporación del suelo a lo largo del tiempo. Sin embargo, el drenaje fue predicho mejor con WAVE que con DSSAT, el cual presentó mayores errores en los valores acumulados. Esto podría ser debido a la naturaleza mecanicista de WAVE frente a la naturaleza más funcional de DSSAT. Los buenos resultados de WAVE indican que, después de la calibración, este puede ser utilizado como "benchmark" para otros modelos para periodos en los que no haya medidas de campo del drenaje. El funcionamiento de DSSAT-CENTURY en la simulación de SOC y N depende fuertemente del proceso de inicialización. Se propuso como método alternativo (Met.2) la inicialización de las fracciones de SOC a partir de medidas de mineralización aparente del suelo (Napmin). El Met.2 se comparó con el método de inicialización de Basso et al. (2011) (Met.1), aplicando ambos métodos a un experimento de campo de 4 años en un área en regadío de España central. Nmin y Napmin fueron sobreestimados con el Met.1, ya que la fracción estable obtenida (SOC3) en las capas superficiales del suelo fue más baja que con Met.2. El N lixiviado simulado fue similar en los dos métodos, con buenos resultados en los tratamientos de barbecho y cebada. El Met.1 subestimó el SOC en la capa superficial del suelo cuando se comparó con una serie observada de 12 años. El crecimiento y rendimiento del cultivo fueron adecuadamente simulados con ambos métodos, pero el N en la parte aérea de la planta y en el grano fueron sobreestimados con el Met.1. Los resultados variaron significativamente con las fracciones iniciales de SOC, resaltando la importancia del método de inicialización. El Met.2 ofrece una alternativa para la inicialización del modelo CENTURY, mejorando la simulación de procesos de N en el suelo. La continua emergencia de nuevas variedades de híbridos modernos de maíz limita la aplicación de modelos de simulación de cultivos, ya que estos nuevos híbridos necesitan ser calibrados en el campo para ser adecuados para su uso en los modelos. El desarrollo de relaciones basadas en la duración del ciclo, simplificaría los requerimientos de calibración facilitando la rápida incorporación de nuevos cultivares en DSSAT. Seis híbridos de maiz (FAO 300 hasta FAO 700) fueron cultivados en un experimento de campo de dos años en un área semiárida de regadío en España central. Los coeficientes genéticos fueron obtenidos secuencialmente, comenzando con los parámetros de desarrollo fenológico (P1, P2, P5 and PHINT), seguido de los parámetros de crecimiento del cultivo (G2 and G3). Se continuó el procedimiento hasta que la salida de las simulaciones estuvo en concordancia con las observaciones fenológicas de campo. Después de la calibración, los parámetros simulados se ajustaron bien a los parámetros observados, con bajos RMSE en todos los casos. Los P1 y P5 calibrados, incrementaron con la duración del ciclo. P1 fue una función lineal del tiempo térmico (TT) desde emergencia hasta floración y P5 estuvo linealmente relacionada con el TT desde floración a madurez. No hubo diferencias significativas en PHINT entre híbridos de FAO-500 a 700 , ya que tuvieron un número de hojas similar. Como los coeficientes fenológicos estuvieron directamente relacionados con la duración del ciclo, sería posible desarrollar rangos y correlaciones que permitan estimar dichos coeficientes a partir de la clasificación del ciclo. ABSTRACT Crop simulation models allow analyzing various tillage-rotation combinations and exploring management scenarios. DSSAT model was tested under rainfed conditions in a 16-year field experiment in semiarid central Spain. The effect of tillage system and winter cereal-based rotations on the crop yield and soil quality was evaluated. The CERES and CROPGRO models were used to simulate crop growth and yield, while the DSSAT CENTURY was used in the SOC and SN simulations. Both field observations and CERES-Barley simulations, showed that barley grain yield was lower for continuous cereal (BB) than for vetch (VB) and fallow (FB) rotations for both tillage systems. The model predicted higher nitrogen availability in the conventional tillage (CT) than in the no tillage (NT) leading to a higher yield in the CT. The SOC and SN in the top layer, were higher in NT than in CT, and decreased with depth in both simulated and observed values. The best combinations for the dry land conditions studied were CT-VB and CT-FB, but CT presented lower SN and SOC content than NT. The beneficial effect of NT on SOC and SN under semiarid Mediterranean conditions can be identified by field observations and by crop model simulations. The simulation of the water balance in cropping systems is a useful tool to study how water can be used efficiently. The comparison of DSSAT soil water balance, with a simpler “tipping bucket” approach, with the more mechanistic WAVE model, which integrates Richard’s equation, is a powerful method to assess model performance. The soil parameters were calibrated by using the Simulated Annealing (SA) global optimizing method. A continuous weighing lysimeter in a bare fallow provided the observed values of drainage and evapotranspiration (ET) while soil water content (SW) was supplied by capacitance sensors. Both models performed well after optimizing soil parameters with SA, simulating the soil water balance components for the calibrated period. For the validation period, the optimized models predicted well soil water content and soil evaporation over time. However, drainage was predicted better by WAVE than by DSSAT, which presented larger errors in the cumulative values. That could be due to the mechanistic nature of WAVE against the more functional nature of DSSAT. The good results from WAVE indicate that, after calibration, it could be used as benchmark for other models for periods when no drainage field measurements are available. The performance of DSSAT-CENTURY when simulating SOC and N strongly depends on the initialization process. Initialization of the SOC pools from apparent soil N mineralization (Napmin) measurements was proposed as alternative method (Met.2). Method 2 was compared to the Basso et al. (2011) initialization method (Met.1), by applying both methods to a 4-year field experiment in a irrigated area of central Spain. Nmin and Napmin were overestimated by Met.1, since the obtained stable pool (SOC3) in the upper layers was lower than from Met.2. Simulated N leaching was similar for both methods, with good results in fallow and barley treatments. Method 1 underestimated topsoil SOC when compared with a 12-year observed serial. Crop growth and yield were properly simulated by both methods, but N in shoots and grain were overestimated by Met.1. Results varied significantly with the initial SOC pools, highlighting the importance of the initialization procedure. Method 2 offers an alternative to initialize the CENTURY model, enhancing the simulation of soil N processes. The continuous emergence of new varieties of modern maize hybrids limits the application of crop simulation models, since these new hybrids should be calibrated in the field to be suitable for model use. The development of relationships based on the cycle duration, would simplify the calibration requirements facilitating the rapid incorporation of new cultivars into DSSAT. Six maize hybrids (FAO 300 through FAO 700) were grown in a 2-year field experiment in a semiarid irrigated area of central Spain. Genetic coefficients were obtained sequentially, starting with the phenological development parameters (P1, P2, P5 and PHINT), followed by the crop growth parameters (G2 and G3). The procedure was continued until the simulated outputs were in good agreement with the field phenological observations. After calibration, simulated parameters matched observed parameters well, with low RMSE in most cases. The calibrated P1 and P5 increased with the duration of the cycle. P1 was a linear function of the thermal time (TT) from emergence to silking and P5 was linearly related with the TT from silking to maturity . There were no significant differences in PHINT between hybrids from FAO-500 to 700 , as they had similar leaf number. Since phenological coefficients were directly related with the cycle duration, it would be possible to develop ranges and correlations which allow to estimate such coefficients from the cycle classification.

Decoupling of soil nutrient cycles as a function of aridity in global drylands

The biogeochemical cycles of carbon (C), nitrogen (N) and phosphorus (P) are interlinked by primary production, respiration and decomposition in terrestrial ecosystems. It has been suggested that the C, N and P cycles could become uncoupled under rapid climate change because of the different degrees of control exerted on the supply of these elements by biological and geochemical processes. Climatic controls on biogeochemical cycles are particularly relevant in arid, semi-arid and dry sub-humid ecosystems (drylands) because their biological activity is mainly driven by water availability. The increase in aridity predicted for the twenty-first century in many drylands worldwide may therefore threaten the balance between these cycles, differentially affecting the availability of essential nutrients. Here we evaluate how aridity affects the balance between C, N and P in soils collected from 224 dryland sites from all continents except Antarctica. We find a negative effect of aridity on the concentration of soil organic C and total N, but a positive effect on the concentration of inorganic P. Aridity is negatively related to plant cover, which may favour the dominance of physical processes such as rock weathering, a major source of P to ecosystems, over biological processes that provide more C and N, such as litter decomposition. Our findings suggest that any predicted increase in aridity with climate change will probably reduce the concentrations of N and C in global drylands, but increase that of P. These changes would uncouple the C, N and P cycles in drylands and could negatively affect the provision of key services provided by these ecosystems.

Soil moisture content determines the effect of the urease inhibitor NBPT on N2O emissions

Among the mitigation strategies to prevent nitrogen (N) losses from ureic fertilizers, urease inhibitors (UIs) have been demonstrated to promote high N use efficiency by reducing ammonia (NH3) volatilization. In the last few years, some field experiments have also shown its effectiveness in reducing nitrous oxide (N2O) losses from fertilized soils under conditions of low soil moisture. An incubation experiment was carried out with the aim of assessing the main biotic mechanisms behind N2O emissions once that the UIs N-(n-butyl) thiophosphoric triamid (NBPT) and phenil phosphorodiamidate (PPDA) were applied with Urea (U) under different soil moisture conditions (40, 60 and 80 % water-filled pore space, WFPS). In the same study we tried to analyze to what extent soil WFPS regulates the effect of these inhibitors on N2O emissions. The use of PPDA in our study allowed us to compare the effect of NBPT with that of another commercially available urease inhibitor, aiming to see if the results were inhibitor-specific or not. Based on the results from this experiment, a WFPS (i.e. 60 %) was chosen for a second study (i.e. mesocosm experiment) aiming to assess the efficiency of the UIs to indirectly affect N2O emissions through influencing the pool of soil mineral N. The N2O emissions at 40 % WFPS were almost negligible, being significantly lower from all fertilized treatments than that produced at 60 and 80 % WFPS. When compared to U alone, NBPT+U reduced the N2O emissions at 60 % WFPS but had no effect at 80 % WFPS. The application of PPDA significantly increased the emissions with respect to U at 80 % WFPS whereas no significant effect was found at 60 %. At 80 % WFPS, denitrification was the main source of N2O emissions for all treatments. In the mesocosm study, the application of NBPT+U was an effective strategy to reduce N2O emissions (75 % reduction compared to U alone), due to a lower soil ammonium (NH4 +) content induced by the inhibitor. These results suggest that adequate management of the UI NBPT could provide, under certain soil conditions, an opportunity for mitigation of N2O emissions from fertilized soils.

Management of irrigation frequency and nitrogen fertilization to mitigate GHG and NO emissions from drip-fertigated crops

Drip irrigation combined with split application of fertilizer nitrogen (N) dissolved in the irrigation water (i.e. drip fertigation) is commonly considered best management practice for water and nutrient efficiency. As a consequence, its use is becoming widespread. Some of the main factors (water-filled pore space, NH4+ and NO3−) regulating the emissions of greenhouse gases (i.e. N2O, CO2 and CH4) and NO from agroecosystems can easily be manipulated by drip fertigation without yield penalties. In this study, we tested management options to reduce these emissions in a field experiment with a melon (Cucumis melo L.) crop. Treatments included drip irrigation frequency (weekly/daily) and type of N fertilizer (urea/calcium nitrate) applied by fertigation. Crop yield, environmental parameters, soil mineral N concentrations and fluxes of N2O, NO, CH4 and CO2 were measured during 85 days. Fertigation with urea instead of calcium nitrate increased N2O and NO emissions by a factor of 2.4 and 2.9, respectively (P < 0.005). Daily irrigation reduced NO emissions by 42% (P < 0.005) but increased CO2 emissions by 21% (P < 0.05) compared with weekly irrigation. We found no relation between irrigation frequency and N2O emissions. Based on yield-scaled Global Warming Potential as well as NO cumulative emissions, we conclude that weekly fertigation with a NO3−-based fertilizer is the best option to combine agronomic productivity with environmental sustainability. Our study shows that adequate management of drip fertigation, while contributing to the attainment of water and food security, may provide an opportunity for climate change mitigation.

Biogeochemical indicators of elevated nitrogen deposition in semiarid Mediterranean ecosystems

Nitrogen (N) deposition has doubled the natural N inputs received by ecosystems through biological N fixation and is currently a global problem that is affecting the Mediterranean regions. We evaluated the existing relationships between increased atmospheric N deposition and biogeochemical indicators related to soil chemical factors and cryptogam species across semiarid central, southern, and eastern Spain. The cryptogam species studied were the biocrust-forming species Pleurochaete squarrosa (moss) and Cladonia foliacea (lichen). Sampling sites were chosen in Quercus coccifera (kermes oak) shrublands and Pinus halepensis (Aleppo pine) forests to cover a range of inorganic N deposition representative of the levels found in the Iberian Peninsula (between 4.4 and 8.1 kg N ha(-1) year(-1)). We extended the ambient N deposition gradient by including experimental plots to which N had been added for 3 years at rates of 10, 20, and 50 kg N ha(-1) year(-1). Overall, N deposition (extant plus simulated) increased soil inorganic N availability and caused soil acidification. Nitrogen deposition increased phosphomonoesterase (PME) enzyme activity and PME/nitrate reductase (NR) ratio in both species, whereas the NR activity was reduced only in the moss. Responses of PME and NR activities were attributed to an induced N to phosphorus imbalance and to N saturation, respectively. When only considering the ambient N deposition, soil organic C and N contents were positively related to N deposition, a response driven by pine forests. The PME/NR ratios of the moss were better predictors of N deposition rates than PME or NR activities alone in shrublands, whereas no correlation between N deposition and the lichen physiology was observed. We conclude that integrative physiological measurements, such as PME/NR ratios, measured on sensitive species such as P. squarrosa, can provide useful data for national-scale biomonitoring programs, whereas soil acidification and soil C and N storage could be useful as additional corroborating ecosystem indicators of chronic N pollution.

The NH3-Abating effect of Slurry Injection as Influenced by Soil Moisture in a Semiarid Arable Soil

Accumulation of large volumes of dilute slurries is considered one of the major problems related to intensive farming (Sommer et al., 2004). In the EU-27, more than half of the total N excretion is applied to croplands due to technical advantages for farmers (e.g. reuse of nutrients). However, the N use efficiency of slurries produced by livestock is low, i.e. only 20-52% of the excreted N is recovered by crops. Much of the remainder can be lost into the atmosphere as ammonia (NH3), nitrous oxide (N2O), dinitrogen (N2) and nitrogen oxides (NOx).

Initializing the DSSAT-Century model: inverse calibration of carbon pools from apparent soil N mineralization

The CENTURY soil organic matter model was adapted for the DSSAT (Decision Support System for Agrotechnology Transfer), modular format in order to better simulate the dynamics of soil organic nutrient processes (Gijsman et al., 2002). The CENTURY model divides the soil organic carbon (SOC) into three hypothetical pools: microbial or active material (SOC1), intermediate (SOC2) and the largely inert and stable material (SOC3) (Jones et al., 2003). At the beginning of the simulation, CENTURY model needs a value of SOC3 per soil layer which can be estimated by the model (based on soil texture and management history) or given as an input. Then, the model assigns about 5% and 95% of the remaining SOC to SOC1 and SOC2, respectively. The model performance when simulating SOC and nitrogen (N) dynamics strongly depends on the initialization process. The common methods (e.g. Basso et al., 2011) to initialize SOC pools deal mostly with carbon (C) mineralization processes and less with N. Dynamics of SOM, SOC, and soil organic N are linked in the CENTURY-DSSAT model through the C/N ratio of decomposing material that determines either mineralization or immobilization of N (Gijsman et al., 2002). The aim of this study was to evaluate an alternative method to initialize the SOC pools in the DSSAT-CENTURY model from apparent soil N mineralization (Napmin) field measurements by using automatic inverse calibration (simulated annealing). The results were compared with the ones obtained by the iterative initialization procedure developed by Basso et al., 2011.

Water and nitrogen footprint in an irrigated crop under mineral and organic fertilization

In order to establish rational nitrogen (N) application and reduce groundwater contamination, a clearer understanding of the N distribution through the growing season and its balance is crucial. Excessive doses of N and/or water applied to fertigated crops involve a substantial risk of aquifer contamination by nitrate; but knowledge of N cycling and availability within the soil could assist in avoiding this excess. In central Spain, the main horticultural fertigated crop is the melon type ?piel de sapo¿ and it is cultivated in vulnerable zones to nitrate pollution (Directive 91/676/CEE). However, until few years ago there were not antecedents related to the optimization of nitrogen fertilization together with irrigation. Water and N footprint are indicators that allow assessing the impact generated by different agricultural practices, so they can be used to improve the management strategies in fertigated crop systems. The water footprint distinguishes between blue water (sources of water applied to the crop, like irrigation and precipitation), green water (water used by the crop and stored in the soil), and it is furthermore possible to quantify the impact of pollution by calculating the grey water, which is defined as the volume of polluted water created from the growing and production of crops. On the other hand, the N footprint considers green N (nitrogen consumed by the crops and stored in the soil), blue N (N available for crop, like N applied with mineral and/or organic fertilizers, N applied with irrigation water and N mineralized during the crop period), whereas grey N is the amount of N-NO3- washed from the soil to the aquifer. All these components are expressed as the ratio between the components of water or N footprint and the yield (m3 t-1 or kg N t-1 respectively). The objetives of this work were to evaluate the impact derivated from the use of different fertilizer practices in a melon crop using water and N footprint.