A statistical mechanical approach for the computation of the climatic response to general forcings

The climate belongs to the class of non-equilibrium forced and dissipative systems, for which most results of quasi-equilibrium statistical mechanics, including the fluctuation-dissipation theorem, do not apply. In this paper we show for the first time how the Ruelle linear response theory, developed for studying rigorously the impact of perturbations on general observables of non-equilibrium statistical mechanical systems, can be applied with great success to analyze the climatic response to general forcings. The crucial value of the Ruelle theory lies in the fact that it allows to compute the response of the system in terms of expectation values of explicit and computable functions of the phase space averaged over the invariant measure of the unperturbed state. We choose as test bed a classical version of the Lorenz 96 model, which, in spite of its simplicity, has a well-recognized prototypical value as it is a spatially extended one-dimensional model and presents the basic ingredients, such as dissipation, advection and the presence of an external forcing, of the actual atmosphere. We recapitulate the main aspects of the general response theory and propose some new general results. We then analyze the frequency dependence of the response of both local and global observables to perturbations having localized as well as global spatial patterns. We derive analytically several properties of the corresponding susceptibilities, such as asymptotic behavior, validity of Kramers-Kronig relations, and sum rules, whose main ingredient is the causality principle. We show that all the coefficients of the leading asymptotic expansions as well as the integral constraints can be written as linear function of parameters that describe the unperturbed properties of the system, such as its average energy. Some newly obtained empirical closure equations for such parameters allow to define such properties as an explicit function of the unperturbed forcing parameter alone for a general class of chaotic Lorenz 96 models. We then verify the theoretical predictions from the outputs of the simulations up to a high degree of precision. The theory is used to explain differences in the response of local and global observables, to define the intensive properties of the system, which do not depend on the spatial resolution of the Lorenz 96 model, and to generalize the concept of climate sensitivity to all time scales. We also show how to reconstruct the linear Green function, which maps perturbations of general time patterns into changes in the expectation value of the considered observable for finite as well as infinite time. Finally, we propose a simple yet general methodology to study general Climate Change problems on virtually any time scale by resorting to only well selected simulations, and by taking full advantage of ensemble methods. The specific case of globally averaged surface temperature response to a general pattern of change of the CO2 concentration is discussed. We believe that the proposed approach may constitute a mathematically rigorous and practically very effective way to approach the problem of climate sensitivity, climate prediction, and climate change from a radically new perspective.

Geodetic mass balance of Azarova glacier, Kodar Mountains, eastern Siberia, and its links to observed and projected climatic change

The Kodar Mountains in eastern Siberia accommodate 30 small, cold-based glaciers with a combined surface area of about 19 km2. Very little is known about these glaciers, with the first survey conducted in the late 1950s. In this paper, we use terrestrial photogrammetry to calculate changes in surface area, elevation, volume and geodetic mass balance of the Azarova Glacier between 1979 and 2007 and relate these to meteorological data from nearby Chara weather station (1938-2007). The glacier surface area declined by 20±6.9% and surface lowered on average by 20±1.8 m (mean thinning: 0.71 m a-1) resulting in a strongly negative cumulative and average mass balance of -18±1.6 m w.e. and -640±60 mm w.e.a-1 respectively. The July-August air temperature increased at a rate of 0.036oC a-1 between 1979 and 2007 and the 1980-2007 period was, on average, around 1oC warmer than 1938-1979. The regional climate projections for A2 and B2 CO2 emission scenarios developed using PRECIS regional climate model indicate that summer temperatures will increase in 2071–2100 by 2.6-4.7°C and 4.9-6.2°C respectively in comparison with 1961–1990. The annual total of solid precipitation will increase by 20% under B2 scenario but decline by 3% under A2 scenario. The length of the ablation season will extend from July–August to June-September. The Azarova Glacier exhibits high sensitivity to climatic warming due to its low elevation, exposure to comparatively high summer temperatures, and the absence of a compensating impact of cold season precipitation. Further summer warming and decline of solid precipitation projected under the A2 scenario will force Azarova to retreat further while impacts of an increase in solid precipitation projected under the B2 scenario require further investigation.

The influence of pomegranate by-product and punicalagins on selected groups of human intestinal microbiota.

We have examined the gut bacterial metabolism of pomegranate by-product (POMx) and major pomegranate polyphenols, punicalagins, using pH-controlled, stirred, batch culture fermentation systems reflective of the distal region of the human large intestine. Incubation of POMx or punicalagins with faecal bacteria resulted in formation of the dibenzopyranone-type urolithins. The time course profile confirmed the tetrahydroxylated urolithin D as the first product of microbial transformation, followed by compounds with decreasing number of phenolic hydroxy groups: the trihydroxy analogue urolithin C and dihydroxylated urolithin A. POMx exposure enhanced the growth of total bacteria, Bifidobacterium spp. and Lactobacillus spp., without influencing the Clostridium coccoides–Eubacterium rectale group and the C. histolyticum group. In addition, POMx increased concentrations of short chain fatty acids (SCFA) viz. acetate, propionate and butyrate in the fermentation medium. Punicalagins did not affect the growth of bacteria or production of SCFA. The results suggest that POMx oligomers, composed of gallic acid, ellagic acid and glucose units, may account for the enhanced growth of probiotic bacteria.

Assessment of projected changes in upland environments using simple climatic indices

Climate controls upland habitats, soils and their associated ecosystem services; therefore, understanding possible changes in upland climatic conditions can provide a rapid assessment of climatic vulnerability over the next century. We used 3 different climatic indices that were optimised to fit the upland area classified by the EU as a Severely Disadvantaged Area (SDA) 1961–1990. Upland areas within the SDA covered all altitudinal ranges, whereas the maximum altitude of lowland areas outside of the SDA was ca. 300 m. In general, the climatic index based on the ratio between annual accumulated temperature (as a measure of growing season length) and annual precipitation predicted 96% of the SDA mapped area, which was slightly better than those indices based on annual or seasonal water deficit. Overall, all climatic indices showed that upland environments were exposed to some degree of change by 2071–2100 under UKCIP02 climate projections for high and low emissions scenarios. The projected area declined by 13 to 51% across 3 indices for the low emissions scenario and by 24 to 84% for the high emissions scenario. Mean altitude of the upland area increased by +11 to +86 m for the low scenario and +21 to +178 m for the high scenario. Low altitude areas in eastern and southern Great Britain were most vulnerable to change. These projected climatic changes are likely to affect upland habitat composition, long-term soil carbon storage and wider ecosystem service provision, although it is not yet possible to determine the rate at which this might occur.

Climatic and oceanic associations with daily rainfall extremes over southern Africa

Changes in climate variability and, in particular, changes in extreme climate events are likely to be of far more significance for environmentally vulnerable regions than changes in the mean state. It is generally accepted that sea-surface temperatures (SSTs) play an important role in modulating rainfall variability. Consequently, SSTs can be prescribed in global and regional climate modelling in order to study the physical mechanisms behind rainfall and its extremes. Using a satellite-based daily rainfall historical data set, this paper describes the main patterns of rainfall variability over southern Africa, identifies the dates when extreme rainfall occurs within these patterns, and shows the effect of resolution in trying to identify the location and intensity of SST anomalies associated with these extremes in the Atlantic and southwest Indian Ocean. Derived from a Principal Component Analysis (PCA), the results also suggest that, for the spatial pattern accounting for the highest amount of variability, extremes extracted at a higher spatial resolution do give a clearer indication regarding the location and intensity of anomalous SST regions. As the amount of variability explained by each spatial pattern defined by the PCA decreases, it would appear that extremes extracted at a lower resolution give a clearer indication of anomalous SST regions.

The influence of carbon surface oxygen groups on Dubinin–Astakhov equation parameters calculated from CO2 adsorption isotherm

We present the results of a systematic study of the influence of carbon surface oxidation on Dubinin–Astakhov isotherm parameters obtained from the fitting of CO2 adsorption data. Using GCMC simulations of adsorption on realistic VPC models differing in porosity and containing the most frequently occurring carbon surface functionalities (carboxyls, hydroxyls and carbonyls) and their mixtures, it is concluded that the maximum adsorption calculated from the DA model is not strongly affected by the presence of oxygen groups. Unfortunately, the same cannot be said of the remaining two parameters of this model i.e. the heterogeneity parameter (n) and the characteristic energy of adsorption (E0). Since from the latter the pore diameters of carbons are usually calculated, by inverse-type relationships, it is concluded that they are questionable for carbons containing surface oxides, especially carboxyls.

Secondary uranium mineralization in southern Finland and its relationship to recent glacial events

Uranium series dating has been carried out on secondary uranyl silicate minerals formed during sub-glacial and post-glacial weathering of Proterozoic uraninite ores in south west Finland. The samples were obtained from two sites adjacent to the Salpauselkä III ice marginal formation and cover a range of depths, from the surface to more than 60 m. Measured ages fall into three distinct groups, 70–100 ka, 28–36 ka and < 2500 yr. The youngest set is associated with surface exposures and the crystals display clear evidence of re-working. The most likely trigger for uranium release at depths below the surface weathering zone is intrusion of oxidising glacial melt water. The latter is often characterised by very high discharge rates along channels, which close once the overpressure generated at the ice margin is released. There is excellent correspondence between the two Finnish sites and published data for similar deposits over a large area of southern and central Sweden. None of the seventy samples analysed gave a U–Th age between 40 and 70 ka; a second hiatus is apparent at 20 ka, coinciding with the Last Glacial Maximum. Thus, the process responsible for uranyl silicate formation was halted for significant periods, owing to a change in geochemical conditions or the hydrogeological regime. These data support the presence of interstadial conditions during the Early and Middle Weichselian since in the absence of major climatic perturbations the uranium phases at depth are stable. When viewed in conjunction with proxy data from mammoth remains it would appear that the region was ice-free prior to the Last Glacial Maximum.

Nickel(II) and copper(II) complexes of Schiff base ligands containing N4O2 and N4S2 donors with pyrrole terminal binding groups: synthesis, characterization, X-ray structures, DFT and electrochemical studies

A series of hexadentate ligands, H2Lm (m = 1−4), [1H-pyrrol-2-ylmethylene]{2-[2-(2-{[1H-pyrrol-2-ylmethylene]amino}phenoxy)ethoxy]phenyl}amine (H2L1), [1H-pyrrol-2-ylmethylene]{2-[4-(2-{[1H-pyrrol-2-ylmethylene]amino}phenoxy)butoxy]phenyl}amine (H2L2), [1H-pyrrol-2-ylmethylene][2-({2-[(2-{[1H-pyrrol-2-ylmethylene]amino}phenyl)thio]ethyl}thio)phenyl]amine (H2L3) and [1H-pyrrol-2-ylmethylene][2-({4-[(2-{[1H-pyrrol-2-lmethylene]amino}phenyl)thio]butyl}thio) phenyl]amine (H2L4) were prepared by condensation reaction of pyrrol-2-carboxaldehyde with {2-[2-(2-aminophenoxy)ethoxy]phenyl}amine, {2-[4-(2-aminophenoxy)butoxy]phenyl}amine, [2-({2-[(2-aminophenyl)thio]ethyl}thio)phenyl]amine and [2-({4-[(2-aminophenyl)thio]butyl}thio)phenyl]amine respectively. Reaction of these ligands with nickel(II) and copper(II) acetate gave complexes of the form MLm (m = 1−4), and the synthesized ligands and their complexes have been characterized by a variety of physico-chemical techniques. The solid and solution states investigations show that the complexes are neutral. The molecular structures of NiL3 and CuL2, which have been determined by single crystal X-ray diffraction, indicate that the NiL3 complex has a distorted octahedral coordination environment around the metal while the CuL2 complex has a seesaw coordination geometry. DFT calculations were used to analyse the electronic structure and simulation of the electronic absorption spectrum of the CuL2 complex using TDDFT gives results that are consistent with the measured spectroscopic behavior of the complex. Cyclic voltammetry indicates that all copper complexes are electrochemically inactive but the nickel complexes with softer thioethers are more easily oxidized than their oxygen analogs.