Data assimilation: From photon counts to Earth System forecasts

Data assimilation – the set of techniques whereby information from observing systems and models is combined optimally – is rapidly becoming prominent in endeavours to exploit Earth Observation for Earth sciences, including climate prediction. This paper explains the broad principles of data assimilation, outlining different approaches (optimal interpolation, three-dimensional and four-dimensional variational methods, the Kalman Filter), together with the approximations that are often necessary to make them practicable. After pointing out a variety of benefits of data assimilation, the paper then outlines some practical applications of the exploitation of Earth Observation by data assimilation in the areas of operational oceanography, chemical weather forecasting and carbon cycle modelling. Finally, some challenges for the future are noted.

Simulating metals and mine discharges in river basins using a new integrated catchment model for metals: pollution impacts and restoration strategies in the Aries-Mures river system in Transylvania, Romania

The INtegrated CAtchment (INCA) model has been developed to simulate the impact of mine discharges on river systems. The model accounts for the key kinetic chemical processes operating as well as the dilution, mixing and redistribution of pollutants in rivers downstream of mine discharges or acid rock drainage sites. The model is dynamic and simulates the day-to-day behaviour of hydrology and eight metals (cadmium, mercury, copper, zinc, lead, arsenic, manganese and chromium) as well as cyanide and ammonia. The model is semi-distributed and can simulate catchments, sub-catchment and in-stream river behaviour. The model has been applied to the Roia Montan Mine in Transylvania, Romania, and used to assess the impacts of old mine adits on the local catchments as well as on the downstream Aries and Mures river system. The question of mine restoration is investigated and a set of clean-up scenarios investigated. It is shown that the planned restoration will generate a much improved water quality from the mine and also alleviate the metal pollution of the river system.

Electrochemical determination of plant-derived leuco-indigo after chemical reduction by glucose

Electrochemical determination of redox active dye species is demonstrated in indigo samples contaminated with high levels of organic and inorganic impurities. The use of a hydrodynamic electrode system based on a vibrating probe (250 Hz, 200 mu m lateral amplitude) allows time-independent diffusion controlled signals to be enhanced and reliable concentration data to be obtained under steady state conditions at relatively fast scan rates up to 4 V s-1In this work the indigo content of a complex plant-derived indigo sample (dye content typically 30%) is determined after indigo is reduced by addition of glucose in aqueous 0.2 M NaOH. The soluble leuco-indigo is measured by its oxidation response at a vibrating electrode. The vibrating electrode, which consisted of a laterally vibrating 500 mu m diameter gold disc, is calibrated with Fe(CN)(6) 3-/4- in 0.1 M KCl and employed for indigo determination at 55, 65, and 75 C in 0.2 M NaOH. Determinations of the indigo content of 25 different samples of plant-derived indigo are compared with those obtained by conventional spectrophotometry. This comparison suggests a significant improvement by the electrochemical method, which appears to be less sensitive to impurities.

Phytohormones and plant-herbivore-pathogen interactions: integrating the molecular with the ecological

Current research into indirect phytopathogen–herbivore interactions (i.e., interactions mediated by the host plant) is carried out in two largely independent directions: ecological/mechanistic and molecular. We investigate the origin of these approaches and their strengths and weaknesses. Ecological studies have determined the effect of herbivores and phytopathogens on their host plants and are often correlative: the need for long-term manipulative experiments is pressing. Molecular/cellular studies have concentrated on the role of signaling pathways for systemic induced resistance, mainly involving salicylic acid and jasmonic acid, and more recently the cross-talk between these pathways. This cross-talk demonstrates how interactions between signaling mechanisms and phytohormones could mediate plant–herbivore–pathogen interactions. A bridge between these approaches may be provided by field studies using chemical induction of defense, or investigating whole-organism mechanisms of interactions among the three species. To determine the role of phytohormones in induced resistance in the field, researchers must combine ecological and molecular methods. We discuss how these methods can be integrated and present the concept of “kaleidoscopic defense.” Our recent molecular-level investigations of interactions between the herbivore Gastrophysa viridula and the rust fungus Uromyces rumicis on Rumex obtusifolius, which were well studied at the mechanistic and ecological levels, illustrate the difficulty in combining these different approaches. We suggest that the choice of the right study system (possibly wild relatives of model species) is important, and that molecular studies must consider the environmental conditions under which experiments are performed. The generalization of molecular predictions to ecologically realistic settings will be facilitated by “middle-ground studies” concentrating on the outcomes of the interactions.

The addition reaction between silylene and ethyne: further isotope studies, pressure dependence studies, and quantum chemical calculations

Time-resolved kinetic studies of the reaction of dideutero-silylene, SiD2, generated by laser flash photolysis of phenylsilane-d(3), have been carried out to obtain rate constants for its bimolecular reaction with C2H2. The reaction was studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas, at five temperatures in the range 297-600 K. The second-order rate constants obtained by extrapolation to the high-pressure limits at each temperature fitted the Arrhenius equation log(k(infinity)/cm(3) molecule(-1) s(-1)) = (-10.05 +/- 0.05) + (3.43 +/- 0.36 kJ mol(-1))/RT ln 10. The rate constants were used to obtain a comprehensive set of isotope effects by comparison with earlier obtained rate constants for the reactions of SiH2 with C2H2 and C2D2. Additionally, pressure-dependent rate constants for the reaction of SiH2 with C2H2 in the presence of He (1-100 Tort) were obtained at 300, 399, and 613 K. Quantum chemical (ab initio) calculations of the SiC2H4 reaction system at the G3 level support the initial formation of silirene, which rapidly isomerizes to ethynylsilane as the major pathway. Reversible formation of vinylsilylene is also an important process. The calculations also indicate the involvement of several other intermediates, not previously suggested in the mechanism. RRKM calculations are in semiquantitative agreement with the pressure dependences and isotope effects suggested by the ab initio calculations, but residual discrepancies suggest the possible involvement of the minor reaction channel, SiH2 + C2H2 - SWPO + C2H4. The results are compared and contrasted with previous studies of this reaction system.

A self-repairing, supramolecular polymer system: healability as a consequence of donor-acceptor pi-pi stacking interactions

A novel supramolecular polymer system, in which the terminal pyrenyl groups of a polyamide intercalate into the chain-folds of a polyimide via electronically-complementary pi-pi stacking, shows both enhanced mechanical properties relative to those of its individual components and facile healing characteristics as a result of the thermoreversibility of non-covalent interactions.

Cross-resistance following artificial selection for increased defense against parasitoids in Drosophila melanogaster

An increase in resistance to one natural enemy may result in no correlated change, a positive correlated change, or a negative correlated change in the ability of the host or prey to resist other natural enemies. The type of specificity is important in understanding the evolutionary response to natural enemies and was studied here in a Drosaphila-parasitoid system. Drosophila melanogaster lines selected for increased larval resistance to the endoparasitoid wasps Asobara tabida or Leptopilina boulardi were exposed to attack by A. tabida, L. boulardi and Leptopilina heterotama at 15 degrees C, 20 degrees C, and 25 degrees C. In general, encapsulation ability increased with temperature, with the exception of the lines selected against L. boulardi, which showed the opposite trend. Lines selected against L, boulardi showed large increases in resistance against all three parasitoid species, and showed similar levels of defense against A. tabida to the lines selected against that parasitoid. In contrast, lines selected against A. tabida showed a large increase in resistance to A. tabida and generally to L. heterotoma, but displayed only a small change in their ability to survive attack by L. boulardi. Such asymmetries in correlated responses to selection for increased resistance to natural enemies may influence host-parasitoid community structure.

PEGylated amyloid peptide nanocontainer delivery and release system

A micellar nanocontainer delivery and release system is designed on the basis of a peptide-polymer conjugate. The hybrid molecules self-assemble into micelles comprising a modified amyloid peptide core surrounded by a PEG corona. The modified amyloid peptide previously studied in our group forms helical ribbons based on a beta-sheet motif and contains beta-amino acids that are excluded from the beta-sheet structure, thus being potentially useful as fibrillization inhibitors. In the model peptide-PEG hybrid system studied, enzymatic degradation using alpha-chymotrypsin leads to selective cleavage close to the PEG-peptide linkage, break up of the micelles, and release of peptides in unassociated form. The release of monomeric peptide is useful because aggregation of the released peptide into beta-sheet amyloid fibrils is not observed. This concept has considerable potential in the targeted delivery of peptides for therapeutic applications.

Chemical and aerosol characterisation of the troposphere over West Africa during the monsoon period as part of AMMA

During June, July and August 2006 five aircraft took part in a campaign over West Africa to observe the aerosol content and chemical composition of the troposphere and lower stratosphere as part of the African Monsoon Multidisciplinary Analysis (AMMA) project. These are the first such measurements in this region during the monsoon period. In addition to providing an overview of the tropospheric composition, this paper provides a description of the measurement strategy (flights performed, instrumental payloads, wing-tip to wing-tip comparisons) and points to some of the important findings discussed in more detail in other papers in this special issue. The ozone data exhibits an "S" shaped vertical profile which appears to result from significant losses in the lower troposphere due to rapid deposition to forested areas and photochemical destruction in the moist monsoon air, and convective uplift of ozone-poor air to the upper troposphere. This profile is disturbed, particularly in the south of the region, by the intrusions in the lower and middle troposphere of air from the southern hemisphere impacted by biomass burning. Comparisons with longer term data sets suggest the impact of these intrusions on West Africa in 2006 was greater than in other recent wet seasons. There is evidence for net photochemical production of ozone in these biomass burning plumes as well as in urban plumes, in particular that from Lagos, convective outflow in the upper troposphere and in boundary layer air affected by nitrogen oxide emissions from recently wetted soils. This latter effect, along with enhanced deposition to the forested areas, contributes to a latitudinal gradient of ozone in the lower troposphere. Biogenic volatile organic compounds are also important in defining the composition both for the boundary layer and upper tropospheric convective outflow. Mineral dust was found to be the most abundant and ubiquitous aerosol type in the atmosphere over Western Africa. Data collected within AMMA indicate that injection of dust to altitudes favourable for long-range transport (i.e. in the upper Sahelian planetary boundary layer) can occur behind the leading edge of mesoscale convective system (MCS) cold-pools. Research within AMMA also provides the first estimates of secondary organic aerosols across the West African Sahel and have shown that organic mass loadings vary between 0 and 2 μg m−3 with a median concentration of 1.07 μg m−3. The vertical distribution of nucleation mode particle concentrations reveals that significant and fairly strong particle formation events did occur for a considerable fraction of measurement time above 8 km (and only there). Very low concentrations were observed in general in the fresh outflow of active MCSs, likely as the result of efficient wet removal of aerosol particles due to heavy precipitation inside the convective cells of the MCSs. This wet removal initially affects all particle size ranges as clearly shown by all measurements in the vicinity of MCSs.

The gas-phase reaction between silylene and 2-butyne: kinetics, isotope studies, pressure dependence studies and quantum chemical calculations

Time-resolved kinetic studies of the reactions of silylene, SiH2, and dideutero-silylene, SiD2, generated by laser. ash photolysis of phenylsilane and phenylsilane-d(3), respectively, have been carried out to obtain rate coefficients for their bimolecular reactions with 2-butyne, CH3C CCH3. The reactions were studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas at five temperatures in the range 294-612 K. The second-order rate coefficients, obtained by extrapolation to the high pressure limits at each temperature, fitted the Arrhenius equations where the error limits are single standard deviations: log(k(H)(infinity)/cm(3) molecule(-1) s(-1)) = (-9.67 +/- 0.04) + (1.71 +/- 0.33) kJ mol(-1)/RTln10 log(k(D)(infinity)/cm(3) molecule(-1) s(-1)) = (-9.65 +/- 0.01) + (1.92 +/- 0.13) kJ mol(-1)/RTln10 Additionally, pressure-dependent rate coefficients for the reaction of SiH2 with 2-butyne in the presence of He (1-100 Torr) were obtained at 301, 429 and 613 K. Quantum chemical (ab initio) calculations of the SiC4H8 reaction system at the G3 level support the formation of 2,3-dimethylsilirene [cyclo-SiH2C(CH3)=C(CH3)-] as the sole end product. However, reversible formation of 2,3-dimethylvinylsilylene [CH3CH=C(CH3)SiH] is also an important process. The calculations also indicate the probable involvement of several other intermediates, and possible products. RRKM calculations are in reasonable agreement with the pressure dependences at an enthalpy value for 2,3-dimethylsilirene fairly close to that suggested by the ab initio calculations. The experimental isotope effects deviate significantly from those predicted by RRKM theory. The differences can be explained by an isotopic scrambling mechanism, involving H - D exchange between the hydrogens of the methyl groups and the D-atoms in the ring in 2,3-dimethylsilirene-1,1-d(2). A detailed mechanism involving several intermediate species, which is consistent with the G3 energy surface, is proposed to account for this.

Coordination behavior of pyridylmethylthioether system with cupric chloride and cupric bromide: CS bond cleavage and crystal structures

The coordination behavior of pyridylmethylthioether type of organic moieties having N2S2 donor set [L-1=1,2-bis(2-pyridylmethylthio)ethane, L-2 = 1,3-bis(2-pyridylmethyl-thio)propane and L-3 = 1,4-bis(2-pyridylmethylthio)butane] with copper(II) chloride and copper(II) bromide have been studied in different chemical environments. Copper(II) chloride assisted C-S bond cleavage of the organic moieties leading to the formation of copper(II) picolinate derivatives, whereas, under similar experimental conditions, no C-S bond cleavage was observed in the reaction with copper(II) bromide. The resulted copper(II) complexes isolated from the different mediums have been characterized by spectroscopic and X-ray crystallographic tools.

Kinetic multi-layer model of gas-particle interactions in aerosols and clouds (KM-GAP): linking condensation, evaporation and chemical reactions of organics, oxidants and water

We present a novel kinetic multi-layer model for gas-particle interactions in aerosols and clouds (KM-GAP) that treats explicitly all steps of mass transport and chemical reaction of semi-volatile species partitioning between gas phase, particle surface and particle bulk. KM-GAP is based on the PRA model framework (Pöschl-Rudich-Ammann, 2007), and it includes gas phase diffusion, reversible adsorption, surface reactions, bulk diffusion and reaction, as well as condensation, evaporation and heat transfer. The size change of atmospheric particles and the temporal evolution and spatial profile of the concentration of individual chemical species can be modelled along with gas uptake and accommodation coefficients. Depending on the complexity of the investigated system, unlimited numbers of semi-volatile species, chemical reactions, and physical processes can be treated, and the model shall help to bridge gaps in the understanding and quantification of multiphase chemistry and microphysics in atmo- spheric aerosols and clouds. In this study we demonstrate how KM-GAP can be used to analyze, interpret and design experimental investigations of changes in particle size and chemical composition in response to condensation, evaporation, and chemical reaction. For the condensational growth of water droplets, our kinetic model results provide a direct link between laboratory observations and molecular dynamic simulations, confirming that the accommodation coefficient of water at 270 K is close to unity. Literature data on the evaporation of dioctyl phthalate as a function of particle size and time can be reproduced, and the model results suggest that changes in the experimental conditions like aerosol particle concentration and chamber geometry may influence the evaporation kinetics and can be optimized for eðcient probing of specific physical effects and parameters. With regard to oxidative aging of organic aerosol particles, we illustrate how the formation and evaporation of volatile reaction products like nonanal can cause a decrease in the size of oleic acid particles exposed to ozone.

Kinetic multi-layer model of gas-particle interactions in aerosols and clouds (KM-GAP): linking condensation, evaporation and chemical reactions of organics, oxidants and water

We present a novel kinetic multi-layer model for gas-particle interactions in aerosols and clouds (KMGAP) that treats explicitly all steps of mass transport and chemical reaction of semi-volatile species partitioning between gas phase, particle surface and particle bulk. KMGAP is based on the PRA model framework (P¨oschl-Rudich- Ammann, 2007), and it includes gas phase diffusion, reversible adsorption, surface reactions, bulk diffusion and reaction, as well as condensation, evaporation and heat transfer. The size change of atmospheric particles and the temporal evolution and spatial profile of the concentration of individual chemical species can be modeled along with gas uptake and accommodation coefficients. Depending on the complexity of the investigated system and the computational constraints, unlimited numbers of semi-volatile species, chemical reactions, and physical processes can be treated, and the model shall help to bridge gaps in the understanding and quantification of multiphase chemistry and microphysics in atmospheric aerosols and clouds. In this study we demonstrate how KM-GAP can be used to analyze, interpret and design experimental investigations of changes in particle size and chemical composition in response to condensation, evaporation, and chemical reaction. For the condensational growth of water droplets, our kinetic model results provide a direct link between laboratory observations and molecular dynamic simulations, confirming that the accommodation coefficient of water at 270K is close to unity (Winkler et al., 2006). Literature data on the evaporation of dioctyl phthalate as a function of particle size and time can be reproduced, and the model results suggest that changes in the experimental conditions like aerosol particle concentration and chamber geometry may influence the evaporation kinetics and can be optimized for efficient probing of specific physical effects and parameters. With regard to oxidative aging of organic aerosol particles, we illustrate how the formation and evaporation of volatile reaction products like nonanal can cause a decrease in the size of oleic acid particles exposed to ozone.

Characterisation of the intracellular-glutamate decarboxylase system (GAD): analysis of its function, transcription and role in the acid resistance of various strains of Listeria monocytogenes

The glutamate decarboxylase (GAD) system is important for the acid resistance of Listeria monocytogenes. We previously showed that under acidic conditions, glutamate (Glt)/γ-aminobutyrate (GABA) antiport is impaired in minimal media but not in rich ones, like brain heart infusion. Here we demonstrate that this behavior is more complex and it is subject to strain and medium variation. Despite the impaired Glt/GABA antiport, cells accumulate intracellular GABA (GABA(i)) as a standard response against acid in any medium, and this occurs in all strains tested. Since these systems can occur independently of one another, we refer to them as the extracellular (GAD(e)) and intracellular (GAD(i)) systems. We show here that GAD(i) contributes to acid resistance since in a ΔgadD1D2 mutant, reduced GABA(i) accumulation coincided with a 3.2-log-unit reduction in survival at pH 3.0 compared to that of wild-type strain LO28. Among 20 different strains, the GAD(i) system was found to remove 23.11% ± 18.87% of the protons removed by the overall GAD system. Furthermore, the GAD(i) system is activated at milder pH values (4.5 to 5.0) than the GAD(e) system (pH 4.0 to 4.5), suggesting that GAD(i) is the more responsive of the two and the first line of defense against acid. Through functional genomics, we found a major role for GadD2 in the function of GAD(i), while that of GadD1 was minor. Furthermore, the transcription of the gad genes in three common reference strains (10403S, LO28, and EGD-e) during an acid challenge correlated well with their relative acid sensitivity. No transcriptional upregulation of the gadT2D2 operon, which is the most important component of the GAD system, was observed, while gadD3 transcription was the highest among all gad genes in all strains. In this study, we present a revised model for the function of the GAD system and highlight the important role of GAD(i) in the acid resistance of L. monocytogenes.