Weak constraints in four-dimensional variational data assimilation

The formulation of four-dimensional variational data assimilation allows the incorporation of constraints into the cost function which need only be weakly satisfied. In this paper we investigate the value of imposing conservation properties as weak constraints. Using the example of the two-body problem of celestial mechanics we compare weak constraints based on conservation laws with a constraint on the background state.We show how the imposition of conservation-based weak constraints changes the nature of the gradient equation. Assimilation experiments demonstrate how this can add extra information to the assimilation process, even when the underlying numerical model is conserving.

Effect of ENSO phase on large-scale snow water equivalent distribution in a GCM

Understanding links between the El Nino-Southern Oscillation (ENSO) and snow would be useful for seasonal forecasting, but also for understanding natural variability and interpreting climate change predictions. Here, a 545-year run of the general circulation model HadCM3, with prescribed external forcings and fixed greenhouse gas concentrations, is used to explore the impact of ENSO on snow water equivalent (SWE) anomalies. In North America, positive ENSO events reduce the mean SWE and skew the distribution towards lower values, and vice versa during negative ENSO events. This is associated with a dipole SWE anomaly structure, with anomalies of opposite sign centered in western Canada and the central United States. In Eurasia, warm episodes lead to a more positively skewed distribution and the mean SWE is raised. Again, the opposite effect is seen during cold episodes. In Eurasia the largest anomalies are concentrated in the Himalayas. These correlations with February SWE distribution are seen to exist from the previous June-July-August (JJA) ENSO index onwards, and are weakly detected in 50-year subsections of the control run, but only a shifted North American response can be detected in the anaylsis of 40 years of ERA40 reanalysis data. The ENSO signal in SWE from the long run could still contribute to regional predictions although it would be a weak indicator only

Water vapor variability in the tropics and its links to dynamics and precipitation

The distribution and variability of water vapor and its links with radiative cooling and latent heating via precipitation are crucial to understanding feedbacks and processes operating within the climate system. Column-integrated water vapor (CWV) and additional variables from the European Centre for Medium-Range Weather Forecasts (ECMWF) 40-year reanalysis (ERA40) are utilized to quantify the spatial and temporal variability in tropical water vapor over the period 1979–2001. The moisture variability is partitioned between dynamical and thermodynamic influences and compared with variations in precipitation provided by the Climate Prediction Center Merged Analysis of Precipitation (CMAP) and the Global Precipitation Climatology Project (GPCP). The spatial distribution of CWV is strongly determined by thermodynamic constraints. Spatial variability in CWV is dominated by changes in the large-scale dynamics, in particular associated with the El Niño–Southern Oscillation (ENSO). Trends in CWV are also dominated by dynamics rather than thermodynamics over the period considered. However, increases in CWV associated with changes in temperature are significant over the equatorial east Pacific when analyzing interannual variability and over the north and northwest Pacific when analyzing trends. Significant positive trends in CWV tend to predominate over the oceans while negative trends in CWV are found over equatorial Africa and Brazil. Links between changes in CWV and vertical motion fields are identified over these regions and also the equatorial Atlantic. However, trends in precipitation are generally incoherent and show little association with the CWV trends. This may in part reflect the inadequacies of the precipitation data sets and reanalysis products when analyzing decadal variability. Though the dynamic component of CWV is a major factor in determining precipitation variability in the tropics, in some regions/seasons the thermodynamic component cancels its effect on precipitation variability.

Summertime precipitation variability over europe and its links to atmospheric dynamics and evaporation

Gridded monthly precipitation data for 1979-2006 from the Global Precipitation Climatology Project are used to investigate interannual summer precipitation variability over Europe and its links to regional atmospheric circulation and evaporation. The first empirical orthogonal function (EOF) mode of European precipitation, explaining 17.2%-22.8% of its total variance, is stable during the summer season and is associated with the North Atlantic Oscillation. The spatialtemporal structure of the second EOF mode is less stable and shows monthtomonth variations during the summer season. This mode is linked to the Scandinavian teleconnection pattern. Analysis of links between leading EOF modes of regional precipitation and evaporation has revealed a significant link between precipitation and evaporation from the European land surface, thus, indicating an important role of the local processes in summertime precipitation variability over Europe. Weaker, but statistically significant links have been found for evaporation from the surface of the Mediterranean and Baltic Seas. Finally, in contrast to winter, no significant links have been revealed between European precipitation and evaporation in the North Atlantic during the summer season.

Trends in aerosol optical depth in the Russian Arctic and their links with synoptic climatology

Temporal and spatial variability of aerosol optical depth (AOD) are examined using observations of direct solar radiation in the Eurasian Arctic for 1940-1990. AOD is estimated using empirical methods for 14 stations located between 66.2 degrees N and 80.6 degrees N, from the Kara Sea to the Chukchi Sea. While AOD exhibits a well-known springtime maximum and summertime minimum at all stations, atmospheric turbidity is higher in spring in the western (Kara-Laptev) part of the Eurasian Arctic. Between June and August, the eastern (East Siberian-Chukchi) sector experiences higher transparency than the western part. A statistically significant positive trend in AOD was observed in the Kara-Laptev sector between the late 1950s and the early 1930s predominantly in spring when pollution-derived aerosol dominates the Arctic atmosphere but not in the eastern sector. Although all stations are remote, those with positive trends are located closer to the anthropogenic sources of air pollution. By contrast, a widespread decline in AOD was observed between 1982 and 1990 in the eastern Arctic in spring but was limited to two sites in the western Arctic. These results suggest that the post-1982 decline in anthropogenic emissions in Europe and the former Soviet Union has had a limited effect on aerosol load in the Arctic. The post-1982 negative trends in AOD in summer, when marine aerosol is present in the atmosphere, were more common in the west. The relationships between AOD and atmospheric circulation are examined using a synoptic climatology approach. In spring, AOD depends primarily on the strength and direction of air flow. Thus strong westerly and northerly flows result in low AOD values in the East Siberian-Chukchi sector. By contrast, strong southerly flow associated with the passage of depressions results in high A OD in the Kara-Laptev sector and trajectory analysis points to the contribution of industrial regions of the sub-Arctic. In summer, low pressure gradient or anticyclonic conditions result in high atmospheric turbidity. The frequency of this weather type has declined significantly since the early 1980s in the Kara-Laptev sector, which partly explains the decline in summer AOD values. (c) 2004 Elsevier B.V. All rights reserved.

How to save a bad element with weak boundary conditions

The P-1-P-1 finite element pair is known to allow the existence of spurious pressure (surface elevation) modes for the shallow water equations and to be unstable for mixed formulations. We show that this behavior is strongly influenced by the strong or the weak enforcement of the impermeability boundary conditions. A numerical analysis of the Stommel model is performed for both P-1-P-1 and P-1(NC)-P-1 mixed formulations. Steady and transient test cases are considered. We observe that the P-1-P-1 element exhibits stable discrete solutions with weak boundary conditions or with fully unstructured meshes. (c) 2005 Elsevier Ltd. All rights reserved.

Formation of dityrosine cross-links during breadmaking

To establish its significance during commercial breadmaking, dityrosine formation was quantified in flours and doughs of six commercial wheat types at various stages of the Chorleywood Bread Process. Dityrosine was formed mainly during mixing and baking, at the levels of nmol/g dry weight. Good breadmaking flours tended to exhibit higher dityrosine content in the final bread than low quality ones, but no relationship was found for dityrosine as a proportion of flour protein content, indicating that the latter was still a dominant factor in the analysis. There was no correlation between gluten yield of the six wheat types and their typical dityrosine concentrations, suggesting that dityrosine crosslinks were not a determinant factor for gluten formation. Ascorbic acid was found to inhibit dityrosine formation during mixing and proving, and have no significant effect on dityrosine in the final bread. Hydrogen peroxide promoted dityrosine formation, which suggests a radical mechanism involving endogenous peroxidases might be the responsible for dityrosine formation during breadmaking.

Multitrophic links between arbuscular mycorrhizal fungi and insect parasitoids

The effects of arbuscular mycorrhizal colonization of Leucanthemum vulgare on parasitism of a leaf-mining insect was studied in a field and a laboratory experiment. In the field, parasitism of Chromatomyia syngenesiae by Diglyphus isaea was lower on mycorrhizal plants, compared with plants where the association was reduced. A laboratory experiment, in which L. vulgare was inoculated with three species of AM fungi, showed that the effects on parasitism rates were mycorrhizal species-dependent. Some fungal combinations increased parasitism, some decreased it, while others had no effect. It is concluded that the most likely cause of these differences is plant size, with parasitoid searching efficiency being reduced on the larger plants, resulting from certain mycorrhizal species combinations. However, a mycorrhizal effect on herbivore-produced plant volatiles cannot be ruled out.

Predicting functional gene links from phylogenetic-statistical analyses of whole genomes

An important element of the developing field of proteomics is to understand protein-protein interactions and other functional links amongst genes. Across-species correlation methods for detecting functional links work on the premise that functionally linked proteins will tend to show a common pattern of presence and absence across a range of genomes. We describe a maximum likelihood statistical model for predicting functional gene linkages. The method detects independent instances of the correlated gain or loss of pairs of proteins on phylogenetic trees, reducing the high rates of false positives observed in conventional across-species methods that do not explicitly incorporate a phylogeny. We show, in a dataset of 10,551 protein pairs, that the phylogenetic method improves by up to 35% on across-species analyses at identifying known functionally linked proteins. The method shows that protein pairs with at least two to three correlated events of gain or loss are almost certainly functionally linked. Contingent evolution, in which one gene's presence or absence depends upon the presence of another, can also be detected phylogenetically, and may identify genes whose functional significance depends upon its interaction with other genes. Incorporating phylogenetic information improves the prediction of functional linkages. The improvement derives from having a lower rate of false positives and from detecting trends that across-species analyses miss. Phylogenetic methods can easily be incorporated into the screening of large-scale bioinformatics datasets to identify sets of protein links and to characterise gene networks.

Wheat archive links long-term population dynamics of pathogens to air pollution

We used the PCR to study the presence of two plant pathogens in archived wheat samples from a long-term experiment started in 1843. The data were used to construct a unique 160-yr time-series of the abundance of Phaeosphaeria nodorum and Mycosphaerella graminicola, two important pathogens of wheat. During the period since 1970, the relative abundance of DNA of these two pathogens in the samples has reflected the relative importance of the two wheat diseases they cause in U.K. disease surveys. Unexpectedly, changes in the ratio of the pathogens over the 160-yr period were very strongly correlated with changes in atmospheric pollution, as measured by SO2 emissions. This finding suggests that long-term, economically important, changes in pathogen populations can be influenced by anthropogenically induced environmental changes.

A new cyano-bridged one-dimensional (GdFeIII)-Fe-III coordination polymer with o-phenanthroline as the blocking ligand: synthesis, structure, and magnetic properties

Synthesis, structural characterization, and magnetic properties of a new cyano-bridged one-dimensional iron (III)-gadolinium (III) compound, trans-[Gd(o-phen)(2)(H2O)(2)(mu-CN)(2)Fe(CN)(4)], - 2no-phen (o-phen = 1,10-phenanthroline), have been described. The compound crystallizes in the triclinic P (1) over bar space group with the following unit cell parameters: a = 10.538(14) angstrom, b = 12.004(14) angstrom, c = 20.61(2) angstrom, alpha = 92.41(1)degrees, beta = 92.76(1)degrees, gamma = 11 2.72(1)degrees, and Z = 2. In this complex, each gadolinium (III) is coordinated to two nitrile nitrogens of the CN groups coming from two different ferricyanides, the mutually trans cyanides of each of which links another different Gd-III to create -NC-Fe(CN)(4)-CN-Gd-NC- type 1-D chain structure. The one-dimensional chains are self-assembled in two-dimensions via weak C-H center dot center dot center dot N hydrogen bonds. Both the variable-temperature (2-300 K, 0.01 T and 0.8 T) and variable-field (0-50 000 Gauss, 2 K) magnetic measurements reveal the existence of very weak interaction in this molecule. The temperature dependence of the susceptibilities has been analyzed using a model for a chain of alternating classic (7/2) and quantum (1/2) spins. (c) 2005 Elsevier B.V. All rights reserved.