Estimating the mean and standard deviation of concentration distributions using censored samples

Environmental data sets of pollutant concentrations in air, water, and soil frequently include unquantified sample values reported only as being below the analytical method detection limit. These values, referred to as censored values, should be considered in the estimation of distribution parameters as each represents some value of pollutant concentration between zero and the detection limit. Most of the currently accepted methods for estimating the population parameters of environmental data sets containing censored values rely upon the assumption of an underlying normal (or transformed normal) distribution. This assumption can result in unacceptable levels of error in parameter estimation due to the unbounded left tail of the normal distribution. With the beta distribution, which is bounded by the same range of a distribution of concentrations, $\rm\lbrack0\le x\le1\rbrack,$ parameter estimation errors resulting from improper distribution bounds are avoided. This work developed a method that uses the beta distribution to estimate population parameters from censored environmental data sets and evaluated its performance in comparison to currently accepted methods that rely upon an underlying normal (or transformed normal) distribution. Data sets were generated assuming typical values encountered in environmental pollutant evaluation for mean, standard deviation, and number of variates. For each set of model values, data sets were generated assuming that the data was distributed either normally, lognormally, or according to a beta distribution. For varying levels of censoring, two established methods of parameter estimation, regression on normal ordered statistics, and regression on lognormal ordered statistics, were used to estimate the known mean and standard deviation of each data set. The method developed for this study, employing a beta distribution assumption, was also used to estimate parameters and the relative accuracy of all three methods were compared. For data sets of all three distribution types, and for censoring levels up to 50%, the performance of the new method equaled, if not exceeded, the performance of the two established methods. Because of its robustness in parameter estimation regardless of distribution type or censoring level, the method employing the beta distribution should be considered for full development in estimating parameters for censored environmental data sets. ^

Polycyclic aromatic compounds (PAHs and oxygenated PAHs) and trace metals in fish species from Ghana (West Africa): Bioaccumulation and health risk assessment

We report the concentrations of 28 PAHs, 15 oxygenated PAHs (OPAHs) and 11 trace metals/metalloids (As, Cd, Co, Cr, Cu, Fe, Mn, Ni, Pb, Se, and Zn) in muscle and gut + gill tissues of demersal fishes (Drapane africana, Cynoglossus senegalensis and Pomadasys peroteti) from three locations along the coast of the Gulf of Guinea (Ghana). The concentrations of ∑ 28PAHs in muscle tissues averaged 192 ng g− 1 dw (range: 71–481 ng g− 1 dw) and were not statistically different between locations. The concentrations of ∑ 28 PAHs were higher in guts + gills than in muscles. The PAH composition pattern was dominated by low molecular weight compounds (naphthalene, alkyl-naphthalenes and phenanthrene). All fish tissues had benzo[a]pyrene concentrations lower than the EU limit for food safety. Excess cancer risk from consumption of some fish was higher than the guideline value of 1 × 10− 6. The concentrations of ∑ 15 OPAHs in fish muscles averaged 422 ng g− 1 dw (range: 28–1715 ng g− 1dw). The ∑ 15 OPAHs/∑ 16 US-EPA PAHs concentration ratio was > 1 in 68% of the fish muscles and 100% of guts + gills. The log-transformed concentrations of PAHs and OPAHs in muscles, guts + gills were significantly (p < 0.05) correlated with their octanol–water partitioning coefficients, strongly suggesting that equilibrium partitioning from water/sediment into fish tissue was the main mechanism of bioaccumulation. The trace metal concentrations in the fish tissues were in the medium range when compared to fish from other parts of the world. The concentrations of some trace metals (Cd, Cu, Fe, Mn, Zn) were higher in guts + gills than in muscle tissues. The target hazard quotients for metals were < 1 and did not indicate a danger to the local population. We conclude that the health risk arising from the consumption of the studied fish (due to their PAHs and trace metals content) is minimal.

Schmiede ─ Heiligtum ─ Wassermühle. Cham-Hagendorn (Kanton Zug) in römischer Zeit. Ausgrabungen 1944/45 und 2003/04

1944/1945 wurde in Cham-Hagendorn eine Wassermühle ausgegraben, die dank ihrer aussergewöhnlich guten Holzerhaltung seit langem einen prominenten Platz in der Forschung einnimmt. 2003 und 2004 konnte die Kantonsarchäologie Zug den Platz erneut archäologisch untersuchen. Dabei wurden nicht nur weitere Reste der Wassermühle, sondern auch Spuren älterer und jüngerer Anlagen geborgen: eine ältere und eine jüngere Schmiedewerkstatt (Horizont 1a/Horizont 3) sowie ein zweiphasiges Heiligtum (Horizonte 1a/1b). All diese Anlagen lassen sich nun in das in den neuen Grabungen erkannte stratigraphische Gerüst einhängen (s. Beil. 2). Dank der Holzerhaltung können die meisten Phasen dendrochronologisch datiert werden (s. Abb. 4.1/1a): Horizont 1a mit Schlagdaten zwischen 162(?)/173 und 200 n. Chr., Horizont 1b um 215/218 n. Chr. und Horizont 2 um 231 n. Chr. Ferner konnten in den neuen Grabungen Proben für mikromorphologische und archäobotanische Untersuchungen entnommen werden (Kap. 2.2; 3.11). In der vorliegenden Publikation werden der Befund und die Baustrukturen vorgelegt, (Kap. 2), desgleichen sämtliche stratifizierten Funde und eine umfassende Auswahl der 1944/1945 geborgenen Funde (Kap. 3). Dank anpassender Fragmente, sog. Passscherben, lassen sich diese zum Teil nachträglich in die Schichtenabfolge einbinden. Die mikromorphologischen und die archäobotanischen Untersuchungen (Kap. 2.2; 3.11) zeigen, dass der Fundplatz in römischer Zeit inmitten einer stark vom Wald und dem Fluss Lorze geprägten Landschaft lag. In unmittelbarer Nähe können weder eine Siedlung noch einzelne Wohnbauten gelegen haben. Die demnach nur gewerblich und sakral genutzten Anlagen standen an einem Bach, der vermutlich mit jenem Bach identisch ist, der noch heute das Groppenmoos entwässert und bei Cham-Hagendorn in die Lorze mündet (s. Abb. 2.4/1). Der antike Bach führte wiederholt Hochwasser ─ insgesamt sind fünf grössere Überschwemmungsphasen auszumachen (Kap. 2.2; 2.4). Wohl anlässlich eines Seehochstandes durch ein Überschwappen der Lorze in den Bach ausgelöst, müssen diese Überschwemmungen eine enorme Gewalt entwickelt haben, der die einzelnen Anlagen zum Opfer fielen. Wie die Untersuchung der Siedlungslandschaft römischer Zeit rund um den Zugersee wahrscheinlich macht (Kap. 6 mit Abb. 6.2/2), dürften die Anlagen von Cham-Hagendorn zu einer in Cham-Heiligkreuz vermuteten Villa gehören, einem von fünf grösseren Landgütern in diesem Gebiet. Hinweise auf Vorgängeranlagen fehlen, mit denen die vereinzelten Funde des 1. Jh. n. Chr. (Kap. 4.5) in Verbindung gebracht werden könnten. Diese dürften eher von einer der Überschwemmungen bachaufwärts weggerissen und nach Cham-Hagendorn eingeschwemmt worden sein. Die Nutzung des Fundplatzes (Horizont 1a; s. Beil. 6) setzte um 170 n. Chr. mit einer Schmiedewerkstatt ein (Kap. 2.5.1). Der Fundanfall, insbesondere die Schmiedeschlacken (Kap. 3.9) belegen, dass hier nur hin und wieder Geräte hergestellt und repariert wurden (Kap. 5.2). Diese Werkstatt war vermutlich schon aufgelassen und dem Verfall preisgegeben, als man 200 n. Chr. (Kap. 4.2.4) auf einer Insel zwischen dem Bach und einem Lorzearm ein Heiligtum errichtete (Kap. 5.3). Beleg für den sakralen Status dieser Insel ist in erster Linie mindestens ein eigens gepflanzter Pfirsichbaum, nachgewiesen mit Pollen, einem Holz und über 400 Pfirsichsteinen (Kap. 3.11). Die im Bach verlaufende Grenze zwischen dem sakralen Platz und der profanen Umgebung markierte man zusätzlich mit einer Pfahlreihe (Kap. 2.5.3). In diese war ein schmaler Langbau integriert (Kap. 2.5.2), der an die oft an Temenosmauern antiker Heiligtümer angebauten Portiken erinnert und wohl auch die gleiche Funktion wie diese gehabt hatte, nämlich das Aufbewahren von Weihegaben und Kultgerät (Kap. 5.3). Das reiche Fundmaterial, das sich in den Schichten der ersten Überschwemmung fand (s. Abb. 5./5), die um 205/210 n. Chr. dieses Heiligtum zerstört hatte, insbesondere die zahlreiche Keramik (Kap. 3.2.4), und die zum Teil auffallend wertvollen Kleinfunde (Kap. 3.3.3), dürften zum grössten Teil einst in diesem Langbau untergebracht gewesen sein. Ein als Glockenklöppel interpretiertes, stratifiziertes Objekt spricht dafür, dass die fünf grossen, 1944/1945 als Stapel aufgefundenen Eisenglocken vielleicht auch dem Heiligtum zuzuweisen sind (Kap. 3.4). In diesen Kontext passen zudem die überdurchschnittlich häufig kalzinierten Tierknochen (Kap. 3.10). Nach der Überschwemmung befestigte man für 215 n. Chr. (Kap. 4.2.4) das unterspülte Bachufer mit einer Uferverbauung (Kap. 2.6.1). Mit dem Bau eines weiteren, im Bach stehenden Langbaus (Kap. 2.6.2) stellte man 218 n. Chr. das Heiligtum auf der Insel in ähnlicher Form wieder her (Horizont 1b; s. Beil. 7). Von der Pfahlreihe, die wiederum die sakrale Insel von der profanen Umgebung abgrenzte, blieben indes nur wenige Pfähle erhalten. Dennoch ist der sakrale Charakter der Anlage gesichert. Ausser dem immer noch blühenden Pfirsichbaum ist es ein vor dem Langbau aufgestelltes Ensemble von mindestens 23 Terrakottafigurinen (s. Abb. 3.6/1), elf Veneres, zehn Matres, einem Jugendlichen in Kapuzenmantel und einem kindlichen Risus (Kap. 3.6; s. auch Kap. 2.6.3). In den Sedimenten der zweiten Überschwemmung, der diese Anlage um 225/230 n. Chr. zum Opfer gefallen war, fanden sich wiederum zahlreiche Keramikgefässe (Kap. 3.2.4) und zum Teil wertvolle Kleinfunde wie eine Glasperle mit Goldfolie (Kap. 3.8.2) und eine Fibel aus Silber (Kap. 3.3.3), die wohl ursprünglich im Langbau untergebracht waren (Kap. 5.3.2 mit Abb. 5/7). Weitere Funde mit sicherem oder möglichem sakralem Charakter finden sich unter den 1944/1945 geborgenen Funden (s. Abb. 5/8), etwa ein silberner Fingerring mit Merkurinschrift, ein silberner Lunula-Anhänger, eine silberne Kasserolle (Kap. 3.3.3), eine Glasflasche mit Schlangenfadenauflage (Kap. 3.8.2) und einige Bergkristalle (Kap. 3.8.4). Im Bereich der Terrakotten kamen ferner mehrere Münzen (Kap. 3.7) zum Vorschein, die vielleicht dort niedergelegt worden waren. Nach der zweiten Überschwemmung errichtete man um 231 n. Chr. am Bach eine Wassermühle (Horizont 2; Kap. 2.7; Beil. 8; Abb. 2.7/49). Ob das Heiligtum auf der Insel wieder aufgebaut oder aufgelassen wurde, muss mangels Hinweisen offen bleiben. Für den abgehobenen Zuflusskanal der Wassermühle verwendete man mehrere stehen gebliebene Pfähle der vorangegangenen Anlagen der Horizonte 1a und 1b. Obwohl die Wassermühle den 28 jährlichen Überschwemmungshorizonten (Kap. 2.2) und den Funden (Kap. 4.3.2; 4.4.4; 45) zufolge nur bis um 260 n. Chr., während gut einer Generation, bestand, musste sie mindestens zweimal erneuert werden – nachgewiesen sind drei Wasserräder, drei Mühlsteinpaare und vermutlich drei Podeste, auf denen jeweils das Mahlwerk ruhte. Grund für diese Umbauten war wohl der weiche, instabile Untergrund, der zu Verschiebungen geführt hatte, so dass das Zusammenspiel von Wellbaum bzw. Sternnabe und Übersetzungsrad nicht mehr funktionierte und das ganze System zerbrach. Die Analyse von Pollen aus dem Gehhorizont hat als Mahlgut Getreide vom Weizentyp nachgewiesen (Kap. 3.11.4). Das Abzeichen eines Benefiziariers (Kap. 3.3.2 mit Abb. 3.3/23,B71) könnte dafür sprechen, dass das verarbeitete Getreide zumindest zum Teil für das römische Militär bestimmt war (s. auch Kap. 6.2.3). Ein im Horizont 2 gefundener Schreibgriffel und weitere stili sowie eine Waage für das Wägen bis zu 35-40 kg schweren Waren aus dem Fundbestand von 1944/1945 könnten davon zeugen, dass das Getreide zu wägen und zu registrieren war (Kap. 3.4.2). Kurz nach 260 n. Chr. fiel die Wassermühle einem weiteren Hochwasser zum Opfer. Für den folgenden Horizont 3 (Beil. 9) brachte man einen Kiesboden ein und errichtete ein kleines Gebäude (Kap. 2.8). Hier war wohl wiederum eine Schmiede untergebracht, wie die zahlreichen Kalottenschlacken belegen (Kap. 3.9), die im Umfeld der kleinen Baus zum Vorschein kamen. Aufgrund der Funde (Kap. 4.4.4; 4.5) kann diese Werkstatt nur kurze Zeit bestanden haben, höchstens bis um 270 n. Chr., bevor sie einem weiteren Hochwasser zum Opfer fiel. Von der jüngsten Anlage, die wohl noch in römische Zeit datiert (Horizont 4; Beil. 10), war lediglich eine Konstruktion aus grossen Steinplatten zu fassen (Kap. 2.9.1). Wozu sie diente, muss offen bleiben. Auch der geringe Fundanfall spricht dafür, dass die Nutzung des Platzes, zumindest für die römische Zeit, allmählich ein Ende fand (Kap. 4.5). Zu den jüngsten Strukturen gehören mehrere Gruben (Kap. 2.9.2), die vielleicht der Lehmentnahme dienten. Mangels Funden bleibt ihre Datierung indes ungewiss. Insbesondere wissen wir nicht, ob sie noch in römische Zeit datieren oder jünger sind. Spätestens mit der fünften Überschwemmung, die zur endgültigen Verlandung führte und wohl schon in die frühe Neuzeit zu setzen ist, wurde der Platz aufgelassen und erst mit dem Bau der bestehenden Fensterfabrik Baumgartner wieder besetzt.

Hyalotekite, (Ba,Pb,K)(4)(Ca,Y)(2)Si-8(B,Be)(2)(Si,B)(2)O28F, a Tectosilicate Related to Scapolite: New Structure Refinement, Phase Transitions and a Short-Range Ordered 3B Superstructure

Hyalotekite, a framework silicate of composition (Ba,Pb,K)(4)(Ca,Y)(2)Si-8(B,Be)(2) (Si,B)(2)O28F, is found in relatively high-temperature(greater than or equal to 500 degrees C) Mn skarns at Langban, Sweden, and peralkaline pegmatites at Dara-i-Pioz, Tajikistan. A new paragenesis at Dara-i-Pioz is pegmatite consisting of the Ba borosilicates leucosphenite and tienshanite, as well as caesium kupletskite, aegirine, pyrochlore, microcline and quartz. Hyalotekite has been partially replaced by barylite and danburite. This hyalotekite contains 1.29-1.78 wt.% Y2O3, equivalent to 0.172-0.238 Y pfu or 8-11% Y on the Ca site; its Pb/(Pb+Ba) ratio ranges 0.36-0.44. Electron microprobe F contents of Langban and Dara-i-Pioz hyalotekite range 1.04-1.45 wt.%, consistent with full occupancy of the F site. A new refinement of the structure factor data used in the original structural determination of a Langban hyalotekite resulted in a structural formula, (Pb1.96Ba1.86K0.18)Ca-2(B1.76Be0.24)(Si1.56B0.44)Si8O28F, consistent with chemical data and all cations with positive-definite thermal parameters, although with a slight excess of positive charge (+57.14 as opposed to the ideal +57.00). An unusual feature of the hyalotekite framework is that 4 of 28 oxygens are non-bridging; by merging these 4 oxygens into two, the framework topology of scapolite is obtained. The triclinic symmetry of hyalotekite observed at room temperature is obtained from a hypothetical tetragonal parent structure via a sequence of displacive phase transitions. Some of these transitions are associated with cation ordering, either Pb-Ba ordering in the large cation sites, or B-Be and Si-B ordering on tetrahedral sites. Others are largely displacive but affect the coordination of the large cations (Pb, Ba, K, Ca). High-resolution electron microscopy suggests that the undulatory extinction characteristic of hyalotekite is due to a fine mosaic microstructure. This suggests that at least one of these transitions occurs in nature during cooling, and that it is first order with a large volume change. A diffuse superstructure observed by electron diffraction implies the existence of a further stage of short-range cation ordering which probably involves both (Pb,K)-Ba and (BeSi,BB)-BSi.

Diurnal variations in middle-atmospheric water vapor by ground-based microwave radiometry

Abstract. In this paper, we compare the diurnal variations in middle-atmospheric water vapor as measured by two ground-based microwave radiometers in the Alpine region near Bern, Switzerland. The observational data set is also compared to data from the chemistry–climate model WACCM. Due to the small diurnal variations of usually less than 1%, averages over extended time periods are required. Therefore, two time periods of five months each, December to April and June to October, were taken for the comparison. The diurnal variations from the observational data agree well with each other in amplitude and phase. The linear correlation coefficients range from 0.8 in the upper stratosphere to 0.5 in the upper mesosphere. The observed diurnal variability is significant at all pressure levels within the sensitivity of the instruments. Comparing our observations with WACCM, we find that the agreement of the phase of the diurnal cycle between observations and model is better from December to April than from June to October. The amplitudes of the diurnal variations for both time periods increase with altitude in WACCM, but remain approximately constant at 0.05 ppm in the observations. The WACCM data are used to separate the processes that lead to diurnal variations in middle-atmospheric water vapor above Bern. The dominating processes were found to be meridional advection below 0.1 hPa, vertical advection between 0.1 and 0.02 hPa and (photo-)chemistry above 0.02 hPa. The contribution of zonal advection is small. The highest diurnal variations in water vapor as seen in the WACCM data are found in the mesopause region during the time period from June to October with diurnal amplitudes of 0.2 ppm (approximately 5% in relative units).

The response of four winter wheat (Triticum aestivum L.) varieties to field water deficit: leaf anti-oxidative protection and proteolytic activity

Introduction: Drought is one of the most significant factors that limit plant productivity. Oxidative stress is a secondary event in many unfavorable environmental conditions. Intracellular proteases have a role in the metabolism reorganisation and nutrient remobilization under stress. In order to under stand the relative significance of oxidative stress and proteolysis in the yield reduction under drought, four varieties of Triticum aestivum L. with different field drought resistance were examined. Methods: A two-year field experiment was conducted. Analyses were performed on the upper most leaf of control plants and plants under water deficitat the stages most critical for yield reduction under drought (from jointing till milk ripeness). Leaf water deficit and electrolyte leakage, malondyaldehyde level, activities and isoenzymes of superoxide dismutase, catalase and peroxidase, leaf protein content and proteolytic activity were studied. Yield components were analyzed. Results: A general trend of increasing the membrane in stability and accumulation of lipid hydroperoxides was observed with some differences among varieties, especially under drought. The anti-oxidative enzyme activities were progressively enhanced, as well as the azocaseinolytic activities. The leaf protein content decreased under drought at the last phase. Differences among varieties were observed in the parameters under study. They were compared to yield components` reduction under water deprivation.

New online method for water isotope analysis of speleothem fluid inclusions using laser absorption spectroscopy (WS-CRDS)

A new online method to analyse water isotopes of speleothem fluid inclusions using a wavelength scanned cavity ring down spectroscopy (WS-CRDS) instrument is presented. This novel technique allows us simultaneously to measure hydrogen and oxygen isotopes for a released aliquot of water. To do so, we designed a new simple line that allows the online water extraction and isotope analysis of speleothem samples. The specificity of the method lies in the fact that fluid inclusions release is made on a standard water background, which mainly improves the δ D robustness. To saturate the line, a peristaltic pump continuously injects standard water into the line that is permanently heated to 140 °C and flushed with dry nitrogen gas. This permits instantaneous and complete vaporisation of the standard water, resulting in an artificial water background with well-known δ D and δ18O values. The speleothem sample is placed in a copper tube, attached to the line, and after system stabilisation it is crushed using a simple hydraulic device to liberate speleothem fluid inclusions water. The released water is carried by the nitrogen/standard water gas stream directly to a Picarro L1102-i for isotope determination. To test the accuracy and reproducibility of the line and to measure standard water during speleothem measurements, a syringe injection unit was added to the line. Peak evaluation is done similarly as in gas chromatography to obtain &delta D; and δ18O isotopic compositions of measured water aliquots. Precision is better than 1.5 ‰ for δ D and 0.4 ‰ for δ18O for water measurements for an extended range (−210 to 0 ‰ for δ D and −27 to 0 ‰ for δ18O) primarily dependent on the amount of water released from speleothem fluid inclusions and secondarily on the isotopic composition of the sample. The results show that WS-CRDS technology is suitable for speleothem fluid inclusion measurements and gives results that are comparable to the isotope ratio mass spectrometry (IRMS) technique.

Plant-plant interactions, environmental gradients and plant diversity: A global synthesis of community-level studies

Previous syntheses on the effects of environmental conditions on the outcome of plant-plant interactions summarize results from pairwise studies. However, the upscaling to the community-level of such studies is problematic because of the existence of multiple species assemblages and species-specific responses to both the environmental conditions and the presence of neighbors. We conducted the first global synthesis of community-level studies from harsh environments, which included data from 71 alpine and 137 dryland communities to: (i) test how important are facilitative interactions as a driver of community structure, (ii) evaluate whether we can predict the frequency of positive plant-plant interactions across differing environmental conditions and habitats, and (iii) assess whether thresholds in the response of plant-plant interactions to environmental gradients exists between ``moderate'' and ``extreme'' environments. We also used those community-level studies performed across gradients of at least three points to evaluate how the average environmental conditions, the length of the gradient studied, and the number of points sampled across such gradient affect the form and strength of the facilitation-environmental conditions relationship. Over 25% of the species present were more spatially associated to nurse plants than expected by chance in both alpine and chyland areas, illustrating the high importance of positive plant-plant interactions for the maintenance of plant diversity in these environments. Facilitative interactions were more frequent, and more related to environmental conditions, in alpine than in dryland areas, perhaps because drylands are generally characterized by a larger variety of environmental stress factors and plant functional traits. The frequency of facilitative interactions in alpine communities peaked at 1000 mm of annual rainfall, and globally decreased with elevation. The frequency of positive interactions in dtyland communities decreased globally with water scarcity or temperature annual range. Positive facilitation-drought stress relationships are more likely in shorter regional gradients, but these relationships are obscured in regions with a greater species turnover or with complex environmental gradients. By showing the different climatic drivers and behaviors of plant-plant interactions in dryland and alpine areas, our results will improve predictions regarding the effect of facilitation on the assembly of plant communities and their response to changes in environmental conditions.

Seasonal Variation in the Capacity for Plant Trait Measures to Predict Grassland Carbon and Water Fluxes

There is a need for accurate predictions of ecosystem carbon (C) and water fluxes in field conditions. Previous research has shown that ecosystem properties can be predicted from community abundance-weighted means (CWM) of plant functional traits and measures of trait variability within a community (FDvar). The capacity for traits to predict carbon (C) and water fluxes, and the seasonal dependency of these trait-function relationships has not been fully explored. Here we measured daytime C and water fluxes over four seasons in grasslands of a range of successional ages in southern England. In a model selection procedure, we related these fluxes to environmental covariates and plant biomass measures before adding CWM and FDvar plant trait measures that were scaled up from measures of individual plants grown in greenhouse conditions. Models describing fluxes in periods of low biological activity contained few predictors, which were usually abiotic factors. In more biologically active periods, models contained more predictors, including plant trait measures. Field-based plant biomass measures were generally better predictors of fluxes than CWM and FDvar traits. However, when these measures were used in combination traits accounted for additional variation. Where traits were significant predictors their identity often reflected seasonal vegetation dynamics. These results suggest that database derived trait measures can improve the prediction of ecosystem C and water fluxes. Controlled studies and those involving more detailed flux measurements are required to validate and explore these findings, a worthwhile effort given the potential for using simple vegetation measures to help predict landscape-scale fluxes.

The coupled atmosphere–chemistry–ocean model SOCOL-MPIOM

The newly developed atmosphere–ocean-chemistry-climate model SOCOL-MPIOM is presented by demonstrating the influence of the interactive chemistry module on the climate state and the variability. Therefore, we compare pre-industrial control simulations with (CHEM) and without (NOCHEM) interactive chemistry. In general, the influence of the chemistry on the mean state and the variability is small and mainly restricted to the stratosphere and mesosphere. The largest differences are found for the atmospheric dynamics in the polar regions, with slightly stronger northern and southern winter polar vortices in CHEM. The strengthening of the vortex is related to larger stratospheric temperature gradients, which are attributed to a parametrization of the absorption of ozone and oxygen in the Lyman-alpha, Schumann–Runge, Hartley, and Higgins bands. This effect is parametrized in the version with interactive chemistry only. A second reason for the temperature differences between CHEM and NOCHEM is related to diurnal variations in the ozone concentrations in the higher atmosphere, which are missing in NOCHEM. Furthermore, stratospheric water vapour concentrations differ substantially between the two experiments, but their effect on the temperatures is small. In both setups, the simulated intensity and variability of the northern polar vortex is inside the range of present day observations. Sudden stratospheric warming events are well reproduced in terms of their frequency, but the distribution amongst the winter months is too uniform. Additionally, the performance of SOCOL-MPIOM under changing external forcings is assessed for the period 1600–2000 using an ensemble of simulations driven by a spectral solar forcing reconstruction. The amplitude of the reconstruction is large in comparison to other state-of-the-art reconstructions, providing an upper limit for the importance of the solar signal. In the pre-industrial period (1600–1850) the simulated surface temperature trends are in reasonable agreement with temperature reconstructions, although the multi-decadal variability is more pronounced. This enhanced variability can be attributed to the variability in the solar forcing. The simulated temperature reductions during the Maunder Minimum are in the lowest probability range of the proxy records. During the Dalton Minimum, when also volcanic forcing is an important driver of temperature variations, the agreement is better. In the industrial period from 1850 onward SOCOL-MPIOM overestimates the temperature increase in comparison to observational data sets. Sensitivity simulations show that this overestimation can be attributed to the increasing trend in the solar forcing reconstruction that is used in this study and an additional warming induced by the simulated ozone changes.

ANALYTICAL MODELS OF EXOPLANETARY ATMOSPHERES. I. ATMOSPHERIC DYNAMICS VIA THE SHALLOW WATER SYSTEM

Within the context of exoplanetary atmospheres, we present a comprehensive linear analysis of forced, damped, magnetized shallow water systems, exploring the effects of dimensionality, geometry (Cartesian, pseudo-spherical, and spherical), rotation, magnetic tension, and hydrodynamic and magnetic sources of friction. Across a broad range of conditions, we find that the key governing equation for atmospheres and quantum harmonic oscillators are identical, even when forcing (stellar irradiation), sources of friction (molecular viscosity, Rayleigh drag, and magnetic drag), and magnetic tension are included. The global atmospheric structure is largely controlled by a single key parameter that involves the Rossby and Prandtl numbers. This near-universality breaks down when either molecular viscosity or magnetic drag acts non-uniformly across latitude or a poloidal magnetic field is present, suggesting that these effects will introduce qualitative changes to the familiar chevron-shaped feature witnessed in simulations of atmospheric circulation. We also find that hydrodynamic and magnetic sources of friction have dissimilar phase signatures and affect the flow in fundamentally different ways, implying that using Rayleigh drag to mimic magnetic drag is inaccurate. We exhaustively lay down the theoretical formalism (dispersion relations, governing equations, and time-dependent wave solutions) for a broad suite of models. In all situations, we derive the steady state of an atmosphere, which is relevant to interpreting infrared phase and eclipse maps of exoplanetary atmospheres. We elucidate a pinching effect that confines the atmospheric structure to be near the equator. Our suite of analytical models may be used to develop decisively physical intuition and as a reference point for three-dimensional magnetohydrodynamic simulations of atmospheric circulation.

Liquid–Vapor Interface of Formic Acid Solutions in Salt Water: A Comparison of Macroscopic Surface Tension and Microscopic in Situ X-ray Photoelectron Spectroscopy Measurements

The liquid–vapor interface is difficult to access experimentally but is of interest from a theoretical and applied point of view and has particular importance in atmospheric aerosol chemistry. Here we examine the liquid–vapor interface for mixtures of water, sodium chloride, and formic acid, an abundant chemical in the atmosphere. We compare the results of surface tension and X-ray photoelectron spectroscopy (XPS) measurements over a wide range of formic acid concentrations. Surface tension measurements provide a macroscopic characterization of solutions ranging from 0 to 3 M sodium chloride and from 0 to over 0.5 mole fraction formic acid. Sodium chloride was found to be a weak salting out agent for formic acid with surface excess depending only slightly on salt concentration. In situ XPS provides a complementary molecular level description about the liquid–vapor interface. XPS measurements over an experimental probe depth of 51 Å gave the C 1s to O 1s ratio for both total oxygen and oxygen from water. XPS also provides detailed electronic structure information that is inaccessible by surface tension. Density functional theory calculations were performed to understand the observed shift in C 1s binding energies to lower values with increasing formic acid concentration. Part of the experimental −0.2 eV shift can be assigned to the solution composition changing from predominantly monomers of formic acid to a combination of monomers and dimers; however, the lack of an appropriate reference to calibrate the absolute BE scale at high formic acid mole fraction complicates the interpretation. Our data are consistent with surface tension measurements yielding a significantly more surface sensitive measurement than XPS due to the relatively weak propensity of formic acid for the interface. A simple model allowed us to replicate the XPS results under the assumption that the surface excess was contained in the top four angstroms of solution.

Detection and quantification of Flavobacterium psychrophilum in water and fish tissue samples by quantitative real time PCR.

BACKGROUND Flavobacterium psychrophilum is the agent of Bacterial Cold Water Disease and Rainbow Trout Fry Syndrome, two diseases leading to high mortality. Pathogen detection is mainly carried out using cultures and more rapid and sensitive methods are needed. RESULTS We describe a qPCR technique based on the single copy gene β' DNA-dependent RNA polymerase (rpoC). Its detection limit was 20 gene copies and the quantification limit 103 gene copies per reaction. Tests on spiked spleens with known concentrations of F. psychrophilum (106 to 101 cells per reaction) showed no cross-reactions between the spleen tissue and the primers and probe. Screening of water samples and spleens from symptomless and infected fishes indicated that the pathogen was already present before the outbreaks, but F. psychrophilum was only quantifiable in spleens from diseased fishes. CONCLUSIONS This qPCR can be used as a highly sensitive and specific method to detect F. psychrophilum in different sample types without the need for culturing. qPCR allows a reliable detection and quantification of F. psychrophilum in samples with low pathogen densities. Quantitative data on F. psychrophilum abundance could be useful to investigate risk factors linked to infections and also as early warning system prior to potential devastating outbreak.

A thermodynamic model for di-trioctahedral chlorite from experimental and natural data in the system MgO-FeO-Al2O3-SiO2-H2O. Applications to P-T sections and geothermometry

We present a new thermodynamic activity-composition model for di-trioctahedral chlorite in the system FeO–MgO–Al2O3–SiO2–H2O that is based on the Holland–Powell internally consistent thermodynamic data set. The model is formulated in terms of four linearly independent end-members, which are amesite, clinochlore, daphnite and sudoite. These account for the most important crystal-chemical substitutions in chlorite, the Fe–Mg, Tschermak and di-trioctahedral substitution. The ideal part of end-member activities is modeled with a mixing-on-site formalism, and non-ideality is described by a macroscopic symmetric (regular) formalism. The symmetric interaction parameters were calibrated using a set of 271 published chlorite analyses for which robust independent temperature estimates are available. In addition, adjustment of the standard state thermodynamic properties of sudoite was required to accurately reproduce experimental brackets involving sudoite. This new model was tested by calculating representative P–T sections for metasediments at low temperatures (<400 °C), in particular sudoite and chlorite bearing metapelites from Crete. Comparison between the calculated mineral assemblages and field data shows that the new model is able to predict the coexistence of chlorite and sudoite at low metamorphic temperatures. The predicted lower limit of the chloritoid stability field is also in better agreement with petrological observations. For practical applications to metamorphic and hydrothermal environments, two new semi-empirical chlorite geothermometers named Chl(1) and Chl(2) were calibrated based on the chlorite + quartz + water equilibrium (2 clinochlore + 3 sudoite = 4 amesite + 4 H2O + 7 quartz). The Chl(1) thermometer requires knowledge of the (Fe3+/ΣFe) ratio in chlorite and predicts correct temperatures for a range of redox conditions. The Chl(2) geothermometer which assumes that all iron in chlorite is ferrous has been applied to partially recrystallized detrital chlorite from the Zone houillère in the French Western Alps.

Genetic consequences of habitat fragmentation during a range expansion

We investigate the effect of habitat fragmentation on the genetic diversity of a species experiencing a range expansion. These two evolutionary processes have not been studied yet, at the same time, owing to the difficulties of deriving analytic results for non-equilibrium models. Here we provide a description of their interaction by using extensive spatial and temporal coalescent simulations and we suggest guidelines for a proper genetic sampling to detect fragmentation. To model habitat fragmentation, we simulated a two-dimensional lattice of demes partitioned into groups (patches) by adding barriers to dispersal. After letting a population expand on this grid, we sampled lineages from the lattice at several scales and studied their coalescent history. We find that in order to detect fragmentation, one needs to extensively sample at a local level rather than at a landscape level. This is because the gene genealogy of a scattered sample is less sensitive to the presence of genetic barriers. Considering the effect of temporal changes of fragmentation intensities, we find that at least 10, but often >100, generations are needed to affect local genetic diversity and population structure. This result explains why recent habitat fragmentation does not always lead to detectable signatures in the genetic structure of populations. Finally, as expected, long-distance dispersal increases local genetic diversity and decreases levels of population differentiation, efficiently counteracting the effects of fragmentation.