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.

Influence of nitrogen deficiency on senescence and the amounts of RNA and proteins in wheat leaves

Senescence-associated coordination in amounts of enzymes localized in different cellular compartments were determined in attached leaves of young wheat (Triticum aestivum L. cv. Arina) plants. Senescence was initiated at the time of full leaf elongation based on declines in total RNA and soluble protein. Removal of N from the growth medium just at the time of full leaf elongation enhanced the rate of senescence. Sustained declines in the amount of ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco, EC 4.1.1.39), and a marked decrease in the rbcS transcripts, just after full leaf elongation indicated that Rubisco synthesis/degradation was very sensitive to the onset of senescence. Rubisco activase amount also declined during senescence but the proportion of rca transcript relative to the total poly A RNA pool increased 3-fold during senescence. Thus, continued synthesis of activase may be required to maintain functional Rubisco throughout senescence. N stress led to declines in the amount of proteins located in the chloroplast, the peroxisome and the cytosol. Transcripts of the Clp protease subunits also declined in response to N stress, indicating that Clp is not a senescence-specific protease. In contrast to the other proteins, mitochondrial NADH-glutamate dehydrogenase (EC 1.4.1.2) was relatively stable during senescence and was not affected by N stress. During natural senescence with adequate plant nitrate supply the amount of nitrite reductase (EC 1.7.7.1) increased, and those of glutamine synthetase (EC 1.4.7.1) and glutamate synthase (EC 6.3.1.2) were stable. These results indicated that N assimilatory capacity can continue or even increase during senescence if the substrate supply is maintained. Differential stabilities of proteins, even within the same cellular compartment, indicate that proteolytic activity during senescence must be highly regulated.

Nitrogen isotopic composition and chemical composition of altered oceanic crust from ODP Legs 129 and 185

Knowledge of the subduction input flux of nitrogen (N) in altered oceanic crust (AOC) is critical in any attempt to mass-balance N across arc-trench systems on a global or individual-margin basis. We have employed sealed-tube, carrier-gas-based methods to examine the N concentrations and isotopic compositions of AOC. Analyses of 53 AOC samples recovered on DSDP/ODP legs from the North and South Pacific, the North Atlantic, and the Antarctic oceans (with larger numbers of samples from Site 801 outboard of the Mariana trench and Site 1149 outboard of the Izu trench), and 14 composites for the AOC sections at Site 801, give N concentrations of 1.3 to 18.2 ppm and d15N_air of -11.6? to +8.3?, indicating significant N enrichment probably during the early stages of hydrothermal alteration of the oceanic basalts. The N-d15N modeling for samples from Sites 801 and 1149 (n=39) shows that the secondary N may come from (1) the sedimentary N in the intercalated sediments and possibly overlying sediments via fluid-sediment/rock interaction, and (2) degassed mantle N2 in seawater via alteration-related abiotic reduction processes. For all Site 801 samples, weak correlation of N and K2O contents indicates that the siting of N in potassic alteration phases strongly depends on N availability and is possibly influenced by highly heterogeneous temperature and redox conditions during hydrothermal alteration. The upper 470-m AOC recovered by ODP Legs 129 and 185 delivers approximately 800 kg/km N annually into the Mariana margin. If the remaining less-altered oceanic crust (assuming 6.5 km, mostly dikes and gabbros) has MORB-like N of 1.5 ppm, the entire oceanic crust transfers 5100 kg/km N annually into that trench. This N input flux is twice as large as the annual N input of 2500 kg/km in seafloor sediments subducting into the same margin, demonstrating that the N input in oceanic crust, and its isotopic consequences, must be considered in any assessment of convergent margin N flux.

Mass of particulate carbon and nitrogen in suspended matter as a function of volume filtered in several oceanic areas

We address two issues in the determination of particulate carbon and nitrogen in suspended matter of aquatic environments. One is the adsorption of dissolved organic matter on filters, leading to overestimate particulate matter. The second is the material loss during filtration due to fragile algal cells breaking up. Examples from both laboratory cultures and natural samples are presented. We recommend using stacked filters in order to estimate thefirst and filtering different volumes of water in order to evaluate the second.

Stable carbon and nitrogen isotope composition of organic matter in sediments of the Oman Margin

Stable isotopic compositions of carbon and nitrogen and organic carbon content of sediments ranging from the Pliocene to the Pleistocene-Holocene in age from the Oman Margin (ODP Sites 724 and 725) are reported. In general, the organic carbon content is greater than 2% at Site 724. Prior to the Pleistocene-Holocene at this site, sediments with higher content of organic matter were deposited owing to favorable preservation conditions and/or higher productivity. In the Pleistocene, lower amounts of organic matter have been preserved; this material generally has more enriched nitrogen isotopic compositions. This may indicate intensification of the Oxygen Minimum Zone and denitrification with the onset of the Pleistocene. A correlation of carbon isotope content of these sediments with oxygen isotope stages at Site 724 indicates an enrichment in 13C during glacial events. Based on the stable isotope evidence of both carbon and nitrogen, there does not appear to be major input of terrigenous-derived allochthonous material in this marine environment. The timing and extent of monsoon winds on the productivity of this region are not evident, but require further studies for collaborative interpretation of small-scale features in the isotopic and carbon content of this environment.

A laboratory guide to the study of qualitative analysis, based upon the application of the theory of electrolytic dissociation and the law of mass action,

Mode of access: Internet.

Effects of nitrogen supply on xylem cytokinin delivery, transpiration and leaf expansion of pea genotypes differing in xylem-cytokinin concentration

The rms2 and rms4 pea ( Pisum sativum L.) branching mutants have higher and lower xylem-cytokinin concentration, respectively, relative to wild type (WT) plants. These genotypes were grown at two levels of nitrogen (N) supply for 18 - 20 d to determine whether or not xylem-cytokinin concentration (X-CK) or delivery altered the transpiration and leaf growth responses to N deprivation. Xylem sap was collected by pressurising de-topped root systems. As sap-flow rate increased, X-CK declined in WT and rms2, but did not change in rms4. When grown at 5.0 mM N, X-CKs of rms2 and rms4 were 36% higher and 6-fold lower, respectively, than WT at sap-flow rates equivalent to whole-plant transpiration. Photoperiod cytokinin (CK) delivery rates ( the product of transpiration and X-CK) decreased more than 6-fold in rms4. Growth of plants at 0.5 mM N had negligible (< 10%) effects on transpiration rates expressed on a leaf area basis in WT and rms4, but decreased transpiration rates of rms2. The low-N treatment decreased leaf expansion by 20 - 25% and expanding leaflet N concentration by 15%. These changes were similar in all genotypes. At sap-flow rates equivalent to whole-plant transpiration, the low N treatment decreased X-CK in rms2 but had no discernible effect in WT and rms4. Since the low N treatment decreased transpiration of all genotypes, photoperiod CK delivery rates also decreased in all genotypes. The similar leaf growth response of all genotypes to N deprivation despite differences in both absolute and relative X-CKs and deliveries suggests that shoot N status is more important in regulating leaf expansion than xylem-supplied cytokinins. The decreased X-CK and transpiration rate of rms2 following N deprivation suggests that changes in xylem-supplied CKs may modify water use.

Decomposition of nitrogen-15 labeled hoop pine harvest residues in subtropical Australia

Information on decomposition of harvest residues may assist in the maintenance of soil fertility in second rotation (2R) hoop pine plantations (Araucaria cunninghamii Aiton ex A. Cunn.) of subtropical Australia. The experiment was undertaken to determine the dynamics of residue decomposition and fate of residue-derived N. We used N-15-labeled hoop pine foliage, branch, and stem material in microplots, over a 30-mo period following harvesting. We examined the decomposition of each component both singly and combined, and used C-13 cross-polarization and magic-angle spinning nuclear magnetic resonance (C-13 CPMAS NMR) to chart C transformations in decomposing foliage. Residue-derived N-15 was immobilized in the 0- to 5-cm soil layer, with approximately 40% N-15 recovery in the soil from the combined residues by the end of the 30-mo period. Total recovery of N-15 in residues and soil varied between 60 and 80% for the combined-residue microplots, with 20 to 40% of the residue N-15 apparently lost. When residues were combined within microplots the rate of foliage decomposition decreased by 30% while the rate of branch and stem decomposition increased by 50 and 40% compared with rates for these components when decomposed separately. Residue decomposition studies should include a combined-residue treatment. Based on C-15 CPMAS NMR spectra for decomposing foliage, we obtained good correlations for methoxyl C, aryl C, carbohydrate C and phenolic C with residue mass, N-15 enrichment, and total N. The ratio of carbohydrate C to methoxyl C may be useful as an indicator of harvest residue decomposition in hoop pine plantations.

Challenges for simultaneous nitrification, denitrification, and phosphorus removal in microbial aggregates: mass transfer limitation and nitrous oxide production

The microbial community composition and activity was investigated in aggregates from a lab-scale bioreactor, in which nitrification, denitrification and phosphorus removal occurred simultaneously. The biomass was highly enriched for polyphosphate accumulating organisms facilitating complete removal of phosphorus from the bulk liquid; however, some inorganic nitrogen still remained at the end of the reactor cycle. This was ascribed to incomplete coupling of nitrification and denitrification causing NO3- accumulation. After 2 h of aeration, denitrification was dependent on the activity of nitrifying bacteria facilitating the formation of anoxic zones in the aggregates; hence, denitrification could not occur without simultaneous nitrification towards the end of the reactor cycle. Nitrous oxide was identified as a product of denitrification, when based on stored PHA as carbon source. This observation is of critical importance to the outlook of applying PHA-driven denitrification in activated sludge processes. (c) 2004 Federation of European Microbiological Societies. Published by Elsevier B.V. All rights reserved.

Reporter ion-based mass spectrometry approaches for the detection of non-enzymatic protein modifications in biological samples

Development of mass spectrometry techniques to detect protein oxidation, which contributes to signalling and inflammation, is important. Label-free approaches have the advantage of reduced sample manipulation, but are challenging in complex samples owing to undirected analysis of large data sets using statistical search engines. To identify oxidised proteins in biological samples, we previously developed a targeted approach involving precursor ion scanning for diagnostic MS3 ions from oxidised residues. Here, we tested this approach for other oxidations, and compared it with an alternative approach involving the use of extracted ion chromatograms (XICs) generated from high-resolution MSMS data using very narrow mass windows. This accurate mass XIC data methodology was effective at identifying nitrotyrosine, chlorotyrosine, and oxidative deamination of lysine, and for tyrosine oxidations highlighted more modified peptide species than precursor ion scanning or statistical database searches. Although some false positive peaks still occurred in the XICs, these could be identified by comparative assessment of the peak intensities. The method has the advantage that a number of different modifications can be analysed simultaneously in a single LC-MSMS run. This article is part of a Special Issue entitled: Posttranslational Protein modifications in biology and Medicine. Biological significance: The use of accurate mass extracted product ion chromatograms to detect oxidised peptides could improve the identification of oxidatively damaged proteins in inflammatory conditions. © 2013 Elsevier B.V.

An operational research-based integrated approach for mass evacuation planning of a city

Large-scale disasters are constantly occurring around the world, and in many cases evacuation of regions of city is needed. ‘Operational Research/Management Science’ (OR/MS) has been widely used in emergency planning for over five decades. Warning dissemination, evacuee transportation and shelter management are three ‘Evacuation Support Functions’ (ESF) generic to many hazards. This thesis has adopted a case study approach to illustrate the importance of integrated approach of evacuation planning and particularly the role of OR/MS models. In the warning dissemination phase, uncertainty in the household’s behaviour as ‘warning informants’ has been investigated along with uncertainties in the warning system. An agentbased model (ABM) was developed for ESF-1 with households as agents and ‘warning informants’ behaviour as the agent behaviour. The model was used to study warning dissemination effectiveness under various conditions of the official channel. In the transportation phase, uncertainties in the household’s behaviour such as departure time (a function of ESF-1), means of transport and destination have been. Households could evacuate as pedestrians, using car or evacuation buses. An ABM was developed to study the evacuation performance (measured in evacuation travel time). In this thesis, a holistic approach for planning the public evacuation shelters called ‘Shelter Information Management System’ (SIMS) has been developed. A generic allocation framework of was developed to available shelter capacity to the shelter demand by considering the evacuation travel time. This was formulated using integer programming. In the sheltering phase, the uncertainty in household shelter choices (either nearest/allocated/convenient) has been studied for its impact on allocation policies using sensitivity analyses. Using analyses from the models and detailed examination of household states from ‘warning to safety’, it was found that the three ESFs though sequential in time, however have lot of interdependencies from the perspective of evacuation planning. This thesis has illustrated an OR/MS based integrated approach including and beyond single ESF preparedness. The developed approach will help in understanding the inter-linkages of the three evacuation phases and preparing a multi-agency-based evacuation planning evacuation