Exploring strategies for coupled 4D-Var data assimilation using an idealised atmosphere-ocean model

Operational forecasting centres are currently developing data assimilation systems for coupled atmosphere-ocean models. Strongly coupled assimilation, in which a single assimilation system is applied to a coupled model, presents significant technical and scientific challenges. Hence weakly coupled assimilation systems are being developed as a first step, in which the coupled model is used to compare the current state estimate with observations, but corrections to the atmosphere and ocean initial conditions are then calculated independently. In this paper we provide a comprehensive description of the different coupled assimilation methodologies in the context of four dimensional variational assimilation (4D-Var) and use an idealised framework to assess the expected benefits of moving towards coupled data assimilation. We implement an incremental 4D-Var system within an idealised single column atmosphere-ocean model. The system has the capability to run both strongly and weakly coupled assimilations as well as uncoupled atmosphere or ocean only assimilations, thus allowing a systematic comparison of the different strategies for treating the coupled data assimilation problem. We present results from a series of identical twin experiments devised to investigate the behaviour and sensitivities of the different approaches. Overall, our study demonstrates the potential benefits that may be expected from coupled data assimilation. When compared to uncoupled initialisation, coupled assimilation is able to produce more balanced initial analysis fields, thus reducing initialisation shock and its impact on the subsequent forecast. Single observation experiments demonstrate how coupled assimilation systems are able to pass information between the atmosphere and ocean and therefore use near-surface data to greater effect. We show that much of this benefit may also be gained from a weakly coupled assimilation system, but that this can be sensitive to the parameters used in the assimilation.

Characterization and identification of mixed-metal phosphates in soils: the application, of Raman spectroscopy

The development of protocols for the identification of metal phosphates in phosphate-treated, metal-contaminated soils is a necessary yet problematical step in the validation of remediation schemes involving immobilization of metals as phosphate phases. The potential for Raman spectroscopy to be applied to the identification of these phosphates in soils has yet to be fully explored. With this in mind, a range of synthetic mixed-metal hydroxylapatites has been characterized and added to soils at known concentrations for analysis using both bulk X-ray powder diffraction (XRD) and Raman spectroscopy. Mixed-metal hydroxylapatites in the binary series Ca-Cd, Ca-Pb, Ca-Sr and Cd-Pb synthesized in the presence of acetate and carbonate ions, were characterized using a range of analytical techniques including XRD, analytical scanning electron microscopy (SEM), infrared spectroscopy (IR), inductively coupled plasma-atomic emission spectrometry (ICP-AES) and Raman spectroscopy. Only the Ca-Cd series displays complete solid solution, although under the synthesis conditions of this study the Cd-5(PO4)(3)OH end member could not be synthesized as a pure phase. Within the Ca-Cd series the cell parameters, IR active modes and Raman active bands vary linearly as a function of Cd content. X-ray diffraction and extended X-ray absorption fine structure spectroscopy (EXAFS) suggest that the Cd is distributed across both the Ca(1) and Ca(2) sites, even at low Cd concentrations. In order to explore the likely detection limits for mixed-metal phosphates in soils for XRD and Raman spectroscopy, soils doped with mixed-metal hydroxylapatites at concentrations of 5, 1 and 0.5 wt.% were then studied. X-ray diffraction could not confirm unambiguously the presence or identity of mixed-metal phosphates in soils at concentrations below 5 wt.%. Raman spectroscopy proved a far more sensitive method for the identification of mixed-metal hydroxylapatites in soils, which could positively identify the presence of such phases in soils at all the dopant concentrations used in this study. Moreover, Raman spectroscopy could also provide an accurate assessment of the degree of chemical substitution in the hydroxylapatites even when present in soils at concentrations as low as 0.1%.

A dendrochemical study of Pinus sylvestris from Siljansfors Experimental Forest, central Sweden

Inductively coupled plasma mass spectrometry (ICP-MS) was used to investigate changes in trace element concentration in two high resolution sequences of tree rings from central Sweden. Individual annual growth increments from 18002002 to 1930-2002 were sampled from two Scots pine (Pinus sylvestris) trees from the Siljansfors Experimental Forest. The aims of the study were: to test the viability of conventional solution induction ICP-MS as a technique for investigating the multi-elemental chemistry of long tree ring sequences at annual resolution, and, to test this specifically with a view to detecting changes in elemental concentrations of Swedish tree rings contemporary with the major (and relatively proximal) Icelandic eruption of Askja (1875). It was found that despite a time consuming sample preparation process, it was possible to use conventional ICP-MS for multi-elemental analysis of a long sequence of tree rings at annual resolution. Although promising data were produced, no truly conclusive concentration anomaly could be detected in the sequence to indicate the impact of the Askja eruption on environmental chemistry. Overall findings underlined the complexity of the tree/environment interaction and the cautious approach to data interpretation essential for any dendrochemical study. (c) 2006 Elsevier Ltd. All rights reserved.

Record of natural and anthropogenic changes in reef environments (Barbados West Indies) using laser ablation ICP-MS and sclerochronology on coral cores

Coral growth rate can be affected by environmental parameters such as seawater temperature, depth, and light intensity. The natural reef environment is also disturbed by human influences such as anthropogenic pollutants, which in Barbados are released close to the reefs. Here we describe a relatively new method of assessing the history of pollution and explain how these effects have influenced the coral communities off the west coast of Barbados. We evaluate the relative impact of both anthropogenic pollutants and natural stresses. Sclerochronology documents framework and skeletal growth rate and records pollution history (recorded as reduced growth) for a suite of sampled Montastraea annularis coral cores. X-radiography shows annual growth band patterns of the corals extending back over several decades and indicates significantly lower growth rate in polluted sites. Results using laser-ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) on the whole sample (aragonite, organic matter, trapped particulate matter, etc.), have shown contrasting concentrations of the trace elements (Cu, Sn, Zn, and Pb) between corals at different locations and within a single coral. Deepwater corals 7 km apart, record different levels of Pb and Sn, suggesting that a current transported the metal pollution in the water. In addition, the 1995 hurricanes are associated with anomalous values for Sn and Cu from most sites. These are believed to result from dispersion of nearshore polluted water. We compared the concentrations of trace elements in the coral growth of particular years to those in the relevant contemporaneous seawater. Mean values for the concentration factor in the coral, relative to the water, ranged from 10 for Cu and Ni to 2.4 and 0.7 for Cd and Zn, respectively. Although the uncertainties are large (60-80%), the coral record enabled us to demonstrate the possibility of calculating a history of seawater pollution for these elements from the 1940s to 1997. Our values were much higher than those obtained from analysis of carefully cleaned coral aragonite; they demonstrate the incorporation of more contamination including that from particulate material as well as dissolved metals.

Diorganoruthenium complexes incorporating noninnocent [C6H2(CH2ER2)2-3,5]22-(E = N, P) Bis-Pincer bridging ligands: synthesis, spectroelectrochemistry, and DFT studies

The dinuclear complex [(tpy)Ru-II(PCP-PCP)Ru-II(tPY)]Cl-2 (bridging PCP-PCP = 3,3',5,5'-tetrakis(diphenylphosphinomethyl)biphenyl, [C6H2(CH2PPh2)(2)-3,5](2)(2-)) was prepared via a transcyclometalation reaction of the bis-pincer ligand [PC(H)P-PC(H)P] and the Ru(II) precursor [Ru(NCN)(tpy)]Cl (NCN = [C6H3(CH2NMe2)(2)-2,6](-)) followed by a reaction with 2,2':6',2 ''-terpyridine (tpy). Electrochemical and spectroscopic properties of [(tpy)Ru-II(PCP-PCP)Ru-II(tPY)]Cl-2 are compared with those of the closely related [(tpy)Ru-II(NCN-NCN)Ru-II(tpy)](PF6)(2) (NCN-NCN = [C6H2(CH2- NMe2)(2)-3,5](2)(2-)) obtained by two-electron reduction of [(tpy)Ru-III(NCN-NCN)Ru-III(tpy)](PF6)(4). The molecular structure of the latter complex has been determined by single-crystal X-ray structure determination. One-electron reduction of [(tpy)Ru-III(NCN-NCN)Ru-III(tpy)](PF6)(4) and one-electron oxidation of [(tpy)Ru-II(PCP-PCP)RUII(tpy)]Cl-2 yielded the mixed-valence species [(tpy)Ru-III(NCN-NCN)RUII(tpy)](3+) and [(tpy)Ru-III(PCP-PCP)RUII(tpy)](3+), respectively. The comproportionation equilibrium constants K-c (900 and 748 for [(tpy)Ru-III(NCN-NCN)Ru-III(tpy)](4+) and [(tpy)Ru-II(PCP-PCP)RUII(tpy)](2+), respectively) determined from cyclic voltammetric data reveal comparable stability of the [Ru-III-Ru-II] state of both complexes. Spectroelectrochemical measurements and near-infrared (NIR) spectroscopy were employed to further characterize the different redox states with special focus on the mixed-valence species and their NIR bands. Analysis of these bands in the framework of Hush theory indicates that the mixed-valence complexes [(tpy)Ru-III(PCP-PCP)RUII(tpy)](3+) and [(tpy)Ru-III(NCN-NCN)RUII(tpy)](3+) belong to strongly coupled borderline Class II/Class III and intrinsically coupled Class III systems, respectively. Preliminary DFT calculations suggest that extensive delocalization of the spin density over the metal centers and the bridging ligand exists. TD-DFT calculations then suggested a substantial MLCT character of the NIR electronic transitions. The results obtained in this study point to a decreased metal-metal electronic interaction accommodated by the double-cyclometalated bis-pincer bridge when strong sigma-donor NMe2 groups are replaced by weak sigma-donor, pi-acceptor PPh2 groups

Dielectric relaxation around the nematic-isotropic transition of liquid crystalline polymers

This study is concerned with a series of acrylate based side-chain liquid crystalline (LC) polymers. Previous studies have shown that these LC polymers have a preference for parallel or perpendicular alignment with respect to the polymer chain which depends on the length of the coupling chain joining the mesogenic unit to the polymer backbone. On the other hand, the dielectric relaxation of these side-chain LC polymers shows a strong relaxation associated to the mesogenic unit dynamics. For samples with parallel alignment, it was found that the dielectric relaxation of the nematic is weaker and broader than the relaxation of the isotropic. By contrast, for samples with perpendicular alignment, the isotropic to nematic transition reduces the broadening the relaxation and increases the relaxation strength. These two features are more evident for samples with short coupling units for which the dielectric relaxation observed appears to be strongly coupled with the backbone dynamics.

Cloning, expression, and characterization of human cytosolic aminopeptidase P: a single manganese(II)-dependent enzyme

The mammalian bradykinin-degrading enzyme aminopeptidase P (AP-P; E. C. 3.4.11.9) is a metal-dependent enzyme and is a member of the peptidase clan MG. AP-P exists as membrane-bound and cytosolic forms, which represent distinct gene products. A partially truncated clone encoding the cytosolic form was obtained from a human pancreatic cDNA library and the 5' region containing the initiating Met was obtained by 5' rapid accumulation of cDNA ends (RACE). The open reading frame encodes a protein of 623 amino acids with a calculated molecular mass of 69,886 Da. The full-length cDNA with a C-terminal hexahistidine tag was expressed in Escherichia coli and COS-1 cells and migrated on SDS-PAGE with a molecular mass of 71 kDa. The expressed cytosolic AP-P hydrolyzed the X-Pro bond of bradykinin and substance P but did not hydrolyze Gly-Pro-hydroxyPro. Hydrolysis of bradykinin was inhibited by 1,10-phenanthroline and by the specific inhibitor of the membrane-bound form of mammalian AP-P, apstatin. Inductively coupled plasma atomic emission spectroscopy of AP-P expressed in E. coli revealed the presence of 1 mol of manganese/mol of protein and insignificant amounts of cobalt, iron, and zinc. The enzymatic activity of AP-P was promoted in the presence of Mn(II), and this activation was increased further by the addition of glutathione. The only other metal ion to cause slight activation of the enzyme was Co(II), with Ca(II), Cu(II), Mg(II), Ni(II), and Zn(II) all being inhibitory. Removal of the metal ion from the protein was achieved by treatment with 1,10-phenanthroline. The metal-free enzyme was reactivated by the addition of Mn(II) and, partially, by Fe(II). Neither Co(II) nor Zn(II) reactivated the metal-free enzyme. On the basis of these data we propose that human cytosolic AP-P is a single metal ion-dependent enzyme and that manganese is most likely the metal ion used in vivo.

Distribution of calcium (Ca) and magnesium (Mg) in the leaves of Brassica rapa under varying exogenous Ca and Mg supply

Background and Aims Leafy vegetable Brassica crops are an important source of dietary calcium (Ca) and magnesium (Mg) and represent potential targets for increasing leaf Ca and Mg concentrations through agronomy or breeding. Although the internal distribution of Ca and Mg within leaves affects the accumulation of these elements, such data are not available for Brassica. The aim of this study was to characterize the internal distribution of Ca and Mg in the leaves of a vegetable Brassica and to determine the effects of altered exogenous Ca and Mg supply on this distribution. Methods Brassica rapa ssp. trilocularis ‘R-o-18’ was grown at four different Ca:Mg treatments for 21 d in a controlled environment. Concentrations of Ca and Mg were determined in fully expanded leaves using inductively coupled plasma-mass spectrometry (ICP-MS). Internal distributions of Ca and Mg were determined in transverse leaf sections at the base and apex of leaves using energy-dispersive X-ray spectroscopy (EDS) with cryo-scanning electron microscopy (cryo-SEM). Key Results Leaf Ca and Mg concentrations were greatest in palisade and spongy mesophyll cells, respectively, although this was dependent on exogenous supply. Calcium accumulation in palisade mesophyll cells was enhanced slightly under high Mg supply; in contrast, Mg accumulation in spongy mesophyll cells was not affected by Ca supply. Conclusions The results are consistent with Arabidopsis thaliana and other Brassicaceae, providing phenotypic evidence that conserved mechanisms regulate leaf Ca and Mg distribution at a cellular scale. The future study of Arabidopsis gene orthologues in mutants of this reference B. rapa genotype will improve our understanding of Ca and Mg homeostasis in plants and may provide a model-to-crop translation pathway for targeted breeding.

Asymmetric forcing from stratospheric aerosols impacts Sahelian rainfall

The Sahelian drought of the 1970s–1990s was one of the largest humanitarian disasters of the past 50 years, causing up to 250,000 deaths and creating 10 million refugees1. It has been attributed to natural variability2–5, overgrazing6 and the impact of industrial emissions of sulphur dioxide7,8. Each mechanism can influence the Atlantic sea surface temperature gradient, which is strongly coupled to Sahelian precipitation2,3. We suggest that sporadic volcanic eruptions in the Northern Hemisphere also strongly influence this gradient and cause Sahelian drought. Using de-trended observations from 1900 to 2010, we show that three of the four driest Sahelian summers were preceded by substantial Northern Hemisphere volcanic eruptions. We use a state-ofthe- art coupled global atmosphere–ocean model to simulate both episodic volcanic eruptions and geoengineering by continuous deliberate injection into the stratosphere. In either case, large asymmetric stratospheric aerosol loadings concentrated in the Northern Hemisphere are a harbinger of Sahelian drought whereas those concentrated in the Southern Hemisphere induce a greening of the Sahel. Further studies of the detailed regional impacts on the Sahel and other vulnerable areas are required to inform policymakers in developing careful consensual global governance before any practical solar radiation management geoengineering scheme is implemented.

Mode-specific chemisorption of CH4 on Pt{110}-(1 x 2) explored by first-principles molecular dynamics

The chemisorption of CH4 on Pt{110}-(1 x 2) has been studied by vibrational analysis of the reaction pathway defined by the potential energy surface and, in time reversal, by first-principles molecular dynamics simulations of CH4 associative desorption, with the electronic structure treated explicitly using density functional theory. We find that the symmetric stretch vibration ν1 is strongly coupled to the reaction coordinate; our results therefore provide a firm theoretical basis for recently reported state-resolved reactivity measurements, which show that excitation of the ν1 normal mode is the most efficient way to enhance the reaction probability

Mode-specificity and transition state-specific energy redistribution in the chemisorption of CH4 on Ni{100}

We have investigated methane (CH4) dissociative chemisorption on the Ni{100} surface by first-principles molecular dynamics (MD) simulations. Our results show that this reaction is mode-specific, with the n1 state being the most strongly coupled to efficient energy flow into the reaction coordinate when the molecule reaches the transition state. By performing MD simulations for two different transition state (TS) structures we provide evidence of TS structure-specific energy redistribution in methane chemisorption. Our results are compared with recently reported state-resolved measurement of methane adsorption probability on nickel surfaces, and we find that a strong correlation exists between the highest vibrational efficacy measured on Ni{100} for the n1 state and the calculated highest fractional vibrational energy content in this mode.

Biomineralisation by earthworms: an investigation into the stability and distribution of amorphous calcium carbonate

Background Many biominerals form from amorphous calcium carbonate (ACC), but this phase is highly unstable when synthesised in its pure form inorganically. Several species of earthworm secrete calcium carbonate granules which contain highly stable ACC. We analysed the milky fluid from which granules form and solid granules for amino acid (by liquid chromatography) and functional group (by Fourier transform infrared (FTIR) spectroscopy) compositions. Granule elemental composition was determined using inductively coupled plasma-optical emission spectroscopy (ICP-OES) and electron microprobe analysis (EMPA). Mass of ACC present in solid granules was quantified using FTIR and compared to granule elemental and amino acid compositions. Bulk analysis of granules was of powdered bulk material. Spatially resolved analysis was of thin sections of granules using synchrotron-based μ-FTIR and EMPA electron microprobe analysis. Results The milky fluid from which granules form is amino acid-rich (≤ 136 ± 3 nmol mg−1 (n = 3; ± std dev) per individual amino acid); the CaCO3 phase present is ACC. Even four years after production, granules contain ACC. No correlation exists between mass of ACC present and granule elemental composition. Granule amino acid concentrations correlate well with ACC content (r ≥ 0.7, p ≤ 0.05) consistent with a role for amino acids (or the proteins they make up) in ACC stabilisation. Intra-granule variation in ACC (RSD = 16%) and amino acid concentration (RSD = 22–35%) was high for granules produced by the same earthworm. Maps of ACC distribution produced using synchrotron-based μ-FTIR mapping of granule thin sections and the relative intensity of the ν2: ν4 peak ratio, cluster analysis and component regression using ACC and calcite standards showed similar spatial distributions of likely ACC-rich and calcite-rich areas. We could not identify organic peaks in the μ-FTIR spectra and thus could not determine whether ACC-rich domains also had relatively high amino acid concentrations. No correlation exists between ACC distribution and elemental concentrations determined by EMPA. Conclusions ACC present in earthworm CaCO3 granules is highly stable. Our results suggest a role for amino acids (or proteins) in this stability. We see no evidence for stabilisation of ACC by incorporation of inorganic components.

Improved seasonal prediction of the hot summer of 2003 over Europe through better representation of uncertainty in the land surface

Methods to explicitly represent uncertainties in weather and climate models have been developed and refined over the past decade, and have reduced biases and improved forecast skill when implemented in the atmospheric component of models. These methods have not yet been applied to the land surface component of models. Since the land surface is strongly coupled to the atmospheric state at certain times and in certain places (such as the European summer of 2003), improvements in the representation of land surface uncertainty may potentially lead to improvements in atmospheric forecasts for such events. Here we analyse seasonal retrospective forecasts for 1981–2012 performed with the European Centre for Medium-Range Weather Forecasts’ (ECMWF) coupled ensemble forecast model. We consider two methods of incorporating uncertainty into the land surface model (H-TESSEL): stochastic perturbation of tendencies, and static perturbation of key soil parameters. We find that the perturbed parameter approach considerably improves the forecast of extreme air temperature for summer 2003, through better representation of negative soil moisture anomalies and upward sensible heat flux. Averaged across all the reforecasts the perturbed parameter experiment shows relatively little impact on the mean bias, suggesting perturbations of at least this magnitude can be applied to the land surface without any degradation of model climate. There is also little impact on skill averaged across all reforecasts and some evidence of overdispersion for soil moisture. The stochastic tendency experiments show a large overdispersion for the soil temperature fields, indicating that the perturbation here is too strong. There is also some indication that the forecast of the 2003 warm event is improved for the stochastic experiments, however the improvement is not as large as observed for the perturbed parameter experiment.