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.

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.

Biomass and essential oil yields of cornmint (Mentha arvensis L.) grown in the irrigation area of Río Dulce, Santiago del Estero, Argentina.

The aim of this study is to determine the yield and composition of the essential oil of cornmint (Mentha arvensis L.) grown in the irrigation area of Santiago del Estero, Argentina. Field tests were carried out under irrigation conditions, harvesting when 70% flowering was reached (in the summer and at the end of the winter seasons). Essential oil yields were 2% in the first cut and 1.6% in the second cut, respectively, the major constituents of the essential oil being menthol, menthone, isomenthone and menthofuran. In both cases, a high concentration of menthol was obtained, although during the winter the content decreased, increasing the concentration of menthofuran. It is concluded that during the summer a higher yield and better quality of essential oil are produced.

Growth dynamics and yield of melon as influenced by nitrogen fertilizer

Nitrogen (N) is an important nutrient for melon (Cucumis melo L.) production. However there is scanty information about the amount necessary to maintain an appropriate balance between growth and yield. Melon vegetative organs must develop sufficiently to intercept light and accumulate water and nutrients but it is also important to obtain a large reproductive-vegetative dry weight ratio to maximize the fruit yield. We evaluated the influence of different N amounts on the growth, production of dry matter and fruit yield of a melon ‘Piel de sapo’ type. A three-year field experiment was carried out from May to September. Melons were subjected to an irrigation depth of 100% crop evapotranspiration and to 11 N fertilization rates, ranging 11 to 393 kg ha –1 in the three years. The dry matter production of leaves and stems increased as the N amount increased. The dry matter of the whole plant was affected similarly, while the fruit dry matter decreased as the N amount was increased above 112, 93 and 95 kg ha –1 , in 2005, 2006 and 2007, respectively. The maximum Leaf Area Index (LAI), 3.1, was obtained at 393 kg ha –1 of N. The lowest N supply reduced the fruit yield by 21%, while the highest increased the vegetative growth, LAI and Leaf Area Duration (LAD), but reduced yield by 24% relative to the N93 treatment. Excessive applications of N increase vegetative growth at the expense of reproductive growth. For this melon type, rates about 90-100 kg ha –1 of N are sufficient for adequate plant growth, development and maximum production. To obtain fruit yield close to the maximum, the leaf N concentration at the end of the crop cycle should be higher than 19.5 g kg –1

Summer deficit irrigation strategies in hedgerow olive orchard cv. ‘Arbequina’: Effect on fruit characteristics and yield.

Maximum production in hedgerow olive orchards is likely not achieved with maximum evapotranspiration over the long-term. Thus, regulated deficit irrigation (RDI) should be considered as a management option. Four irrigation treatments were evaluated during the summer when olive is most drought resistant. Control (CON) was irrigated to maintain the rootzone close to field capacity. Severe water deficit was applied by irrigating 30% CON from end of fruit drop to end July (DI-J) and from end July until beginning of oil synthesis (DI-A). Less severe water deficit was applied during July and August (DI-JA) by irrigating 50% CON. Flowering, fruiting, abscission, fruit development, fresh and dry weight of fruits, and oil production were evaluated. There were not significant differences in number of buds initiated, number of fruits per inflorescence and fruit drop. Oil production was significantly different between irrigation treatments in all experimental years. CON produced more oil and fruit with higher oil% than DI-A and DI-JA. Oil production of DI-J was not significantly reduced compared to CON and oil% was greater. DI-J was the most effective RDI strategy; with 16% less applied water relative to CON average loss in oil production of 8% was not significantly different to CON. While DI-JA saved most water (27%), oil production was reduced by 15%. Greatest loss in oil production (21%) was observed in DI-A with water saving of 22%.

Effects of irrigation and time of application on the yield and quality of cv. Cabernet sauvignon in a warm climate

Mediterranean climate is characterized by hot summer, high evapotranspiration rates, and scarce precipitations (400 mm per year) during grapevine cycle. These extremely dry conditions affect vineyard productivity and sustainability. Supplementary irrigation is a needed practice in order to maintain yield and quality. Almost all Spanish grape growing regions are characterized by these within this context, especially in the center region, where this study was performed. The main objective of this work was to study the influence of irrigation on yield and quality. For this aim, we applied different levels of irrigation (mm of water applied) during different stages of growth and berry maturity. Four experimental treatments were applied considering the amount of water and the moment of the application: T1: Water irrigation (420 mm) applied from bloom to maturity. T2: Corresponded to the traditional irrigation scheduling, from preveraison to maturity (154 mm). T3: Water irrigation from bloom to preveraison, and water deficit from veraison to maturity (312 mm). T4: Irrigation applied from preveraison to maturity (230 mm) Experimental vineyard, cv. Cabernet Sauvignon, was located in a commercial vineyard (Bodegas Licinia S.L.) in the hot region of Morata de Tajuña (Madrid). The trial was performed during 2010 and 2011 seasons. Our results showed that yield increased from 2010 to 2011 in the treatments with a higher amount of water appli ed, T1 and T3 (24 and 10 % of yield increase respectively). This was mainly due to an increase in bud fertility (nº of bunches per shoot). Furthermore, sugar content was higher in T3 (27.3 ºBrix), followed by T2 (27 ºBrix). By contrast, T4 (irrigation from veraison) presented the lowest solid soluble concentration and the highest acidity. These results suggest that grapevine has an intrinsic capacity to adapt to its environment. However, this adaptation capacity should be evaluated considering the sensibility of quality parameters during the maturity period (acidity, pH, aroma, color...) and its impact on yield. Here, we demonstrated that a higher amount of water irrigation applied was no linked to a negative effect on quality.

CsRAV1: a candidate gene for improving lignocellulosic biomass yield

Fast-growing tree species of Populus spp.,Salix spp. and Eucalyptus spp. are cultivated to produce wood in a short time. Poplars are cultivated with cycles of 15-18 years to obtain saw timber and peeler logs, but when grown as short -rotation coppice(SRC) to produce biomass, planting density increases and rotation is considerably reduced (3-5 years). In this regard, research efforts are focused in the identification of traits and loci that allow the generation of improved SRC biomass-yielding genotypes. Biomass yield is a highly complex trait as it is the combined outcome of many other complex traits, each under separate polygenic control. Among profitable biomass yield-related traits are the amount of sylleptic branching and the length of winter dormancy. In poplar and in a few other Salicaceae species some lateral buds grow out sylleptically, the same season in which they form without the need of an intervening rest period. Sylleptic branching in poplar increases branch number, leaf area and general growth of the tree in its early years, and is a reasonable predictor of coppice yield. On the other hand, the length of winter dormancy determines the extent of the growth period. Our group has characterized the RAV1 gene of Castanea sativa (CsRAV1), encoding a transcription factor of the subfamily RAV (Related to ABI3/VP1). CsRAV1 expression shows a marked seasonal pattern, being higher in autumn and winter both in stems and buds. We generated transgenic lines of the hybrid clone Populus tremulax P. alba INRA 717 1B4 constitutively expressing CsRAV 1. These CsRAV1-expressing poplars develop sylleptic branches only a few weeks after potting. In addition to the sylleptic branching phenotype, these trees show phenological features that could give rise to an extended growth period. We are currently assessing the phenotype and behavior of these transgenic trees in a field trial, and ultimately, we will evaluate the impact on lignocellulosic biomass quality and production.

Short-term effects of four tillage practices on soil physical properties, soil water potential, and maize yield

The area cultivated using conservation tillage has recently increased in central Spain. However, soil compaction and water retention with conservation tillage still remains a genuine concern for landowners in this region be- cause of its potential effect on the crop growth and yield. The aim of this research is to determine the short- term influences of four tillage treatments on soil physical properties. In the experiment, bulk density, cone index, soil water potential, soil temperature and maize (Zea mays L.) productivity have been measured. A field experiment was established in spring of 2013 on a loamy soil. The experiment compared four tillage methods (zero tillage, ZT; reservoir tillage, RT; minimum tillage, MT; and conventional tillage, CT). Soil bulk density and soil cone index were measured during maize growing season and at harvesting time. Furthermore, the soil water potential was monitored by using a wireless sensors network with sensors at 20 and 40 cm depths. Also, soil temperatures were registered at depths of 5 and 12 cm. Results indicated that there were significant differ- ences between soil bulk density and cone index of ZT method and those of RT, MT, and CT, during the growing season; although, this difference was not significant at the time of harvesting in some soil layers. Overall, in most soil layers, tillage practice affected bulk density and cone index in the order: ZT N RT N MT N CT. Regardless oftheentireobservationperiod,results exhibited that soils under ZT and RT treatments usually resulted in higher water potential and lower soil temperature than the other two treatments at both soil depths. In addition, clear differences in maize grain yield were observed between ZT and CT treatments, with a grain yield (up to 15.4%) increase with the CT treatment. On the other hand, no significant differences among (RT, MT, and CT) on maizeyieldwerefound.Inconclusion,the impact of soil compaction increase and soil temperature decrease,pro- duced by ZT treatment is a potential reason for maize yield reduction in this tillage method. We found that RT could be certainly a viable option for farmers incentral Spain,particularly when switching to conservation tillage from conventional tillage. This technique showed a moderate and positive effect on soil physical properties and increased maize yields compared to ZT and MT, and provides an opportunity to stabilize maize yields compared to CT.

Current ozone levels threaten gross primary production and yield of Mediterranean annual pastures and nitrogen modulates the response

Pastures are among the most important ecosystems in Europe considering their biodiversity and dis- tribution area. However, their response to increasing tropospheric ozone (O 3 ) and nitrogen (N) deposi- tion, two of the main drivers of global change, is still uncertain. A new Open-Top Chamber (OTC) experiment was performed in central Spain, aiming to study annual pasture response to O 3 and N in close to natural growing conditions. A mixture of six species of three representative families was sowed in the fi eld. Plants were exposed for 40 days to four O 3 treatments: fi ltered air, non- fi ltered air (NFA) repro- ducing ambient levels and NFA supplemented with 20 and 40 nl l � 1 O 3 . Three N treatments were considered to reach the N integrated doses of “ background ” , þ 20 or þ 40 kg N ha � 1 . Ozone signi fi cantly reduced green and total aboveground biomass (maximum reduction 25%) and increased the senescent biomass (maximum increase 40%). Accordingly, O 3 decreased community Gross Primary Production due to both a global reduction of ecosystem CO 2 exchange and an increase of ecosystem respiration. Nitrogen could partially counterbalance O 3 effects on aboveground biomass when the levels of O 3 were moderate, but at the same time O 3 exposure reduced the fertilization effect of higher N availability. Therefore, O 3 must be considered as a stress factor for annual pastures in the Mediterranean areas.

Temperature and precipitation effects on wheat yield across a European transect: a crop model ensemble analysis using impact response surfaces

This study explored the utility of the impact response surface (IRS) approach for investigating model ensemble crop yield responses under a large range of changes in climate. IRSs of spring and winter wheat Triticum aestivum yields were constructed from a 26-member ensemble of process-based crop simulation models for sites in Finland, Germany and Spain across a latitudinal transect. The sensitivity of modelled yield to systematic increments of changes in temperature (-2 to +9°C) and precipitation (-50 to +50%) was tested by modifying values of baseline (1981 to 2010) daily weather, with CO2 concentration fixed at 360 ppm. The IRS approach offers an effective method of portraying model behaviour under changing climate as well as advantages for analysing, comparing and presenting results from multi-model ensemble simulations. Though individual model behaviour occasionally departed markedly from the average, ensemble median responses across sites and crop varieties indicated that yields decline with higher temperatures and decreased precipitation and increase with higher precipitation. Across the uncertainty ranges defined for the IRSs, yields were more sensitive to temperature than precipitation changes at the Finnish site while sensitivities were mixed at the German and Spanish sites. Precipitation effects diminished under higher temperature changes. While the bivariate and multi-model characteristics of the analysis impose some limits to interpretation, the IRS approach nonetheless provides additional insights into sensitivities to inter-model and inter-annual variability. Taken together, these sensitivities may help to pinpoint processes such as heat stress, vernalisation or drought effects requiring refinement in future model development.

Yield and sugar accumulating capacities of Airén cultivar. A preliminary study

Airén is the most worldwide spread white grape cultivar, high yielding, well adapted to hot, dry conditions, and not very sensitive to fungal diseases. Its largest growing region is La Mancha, where Airén has been traditionally bush trained, spur pruned and grown with no irrigation. However, grape growing has evolved to meet the need for higher yields and harvest mechanization; and modern cultural practices train grape vines to simple multi-wire trellis systems, cane pruned, and usually irrigated. The aim of the present study was to evaluate the yield and sugar accumulating capacities of Airén cultivar with regard to leaf area, and to assess the influence that different yield components have on yield. In 2014, five commercial irrigated vineyards, located in La Mancha, of different ages, and grafted onto different rootstocks were selected for this study. Canopy surface area (SA) was measured at maturity. Berry weight and sugar concentration were measured during ripening on a weekly basis. Yield and yield components were determined at harvest. Values for shoot density ranged 2.3-5.1 shoots/m2; SA, 0.6-1.1 m2/m2; yield, 20-40 t/ha; fertility, 1.1-1.7 bunches/shoot; bunch weight, 450-650 g; berry weight, 2.5-2.9 g; and sugar concentration, 17-21 ºBrix. The number of bunches per shoot was the yield component that had the greatest influence on yield. The number of berries was the main contributing factor to bunch weight. A lineal relationship between SA/yield and sugar concentration was observed, with values of SA/yield ranging from 0.20 to 0.45 m2/kg. A ratio SA/yield of approximately 0.4 m2/kg was needed to reach a value of 20 ºBrix. Hence it would be necessary a SA of 12000 m2/ha, under the conditions of this study, to achieve a 30 t/ha yield, and a sugar concentration of 20 ºBrix. These results are a step forward in the study of the Airén cultivar, being of help for grape growers in the center area of Spain in order to maximize crop yield and sugar accumulation.