Evolution of metamorphic volatiles during exhumation of microdiamond-bearing granulites in the Western Gneiss Region, Norway

Fluid inclusions in garnet, kyanite and quartz from microdiamond-bearing granulites in the Western Gneiss Region, Norway, document a conspicuous fluid evolution as the rocks were exhumed following Caledonian high- and ultrahigh-pressure (HP–UHP) metamorphism. The most important of the various fluid mixtures and daughter minerals in these rocks are: (N₂ + CO₂ + magnesian calcite), (N₂ + CO₂ + CH₄ + graphite + magnesian calcite), (N₂ + CH₄), (N₂ + CH₄ + H₂O), (CO₂) and (H₂O + NaCl + CaCl₂ + nahcolite). Rutile also occurs in the N₂ + CO₂ inclusions as a product of titanium diffusion from the garnet host into the fluid inclusions. Volatiles composed of N₂ + CO₂ + magnesian calcite characterise the ambient metamorphic environment between HP–UHP (peak) and early retrograde metamorphism. During progressive decompression, the mole fraction of N₂ increased in the fluid mixtures; as amphibolite-facies conditions were reached, CH₄ and later, H₂O, appeared in the fluids, concomitant with the disappearance of CO₂ and magnesian calcite. Graphite is ubiquitous in the host lithologies and fluid inclusions. Thermodynamic modelling of the metamorphic volatiles in a graphite-buffered C-O-H system demonstrates that the observed metamorphic volatile evolution was attainable only if the f _O2 increased from c. −3.5 (±0.3) to −0.8 (±0.3) log units relative to the FMQ oxygen buffer. External introduction of oxidising aqueous solutions along a system of interconnected ductile shear zones adequately explains the dramatic increase in the f _O2. The oxidising fluids introduced during exhumation were likely derived from dehydration of oceanic crust and continental sediments previously subducted during an extended period of continental collision in conjunction with the Caledonian orogeny. Received: 15 December 1997 / Accepted: 25 May 1998     

Chlorine in basalts from Iceland

Chlorine has a higher solubility in basaltic magmas than other volatiles. A theoretical degassing model predicts that less than 10% of the chlorine originally present in the magma is lost to the atmosphere during surface degassing. This prediction is born out by strong correlation between chlorine and nonvolatile elements in samples from recent volcanic eruptions in Iceland. It is concluded that the chlorine content of subaerial basaltic lavas is proportional to the chlorine content of the magmas. Assuming a roughly constant ratio between chlorine and water in the magmatic gas phase and assuming further that this ratio approximates the $\text{Cl}\text{H}{}_2\text{O}$ ratio in sea water and sediments it is possible to assign each basaltic magma a fixed amount of water. The figures thus obtained (0.3–1.2% H₂O) are identical with previously assumed water contents of basaltic magmas.The implied water contents of the melts are discussed in terms of derivation from hydrous mantle phases. It is concluded that such phases can supply only a fraction of the assumed water content, the rest is derived from an additional volatile source in the mantle. The mantle source giving rise to alkali basalts in Iceland has a lower content of volatiles than the tholeiite source.     

Do paleoclimate proxies agree? A test comparing 19 late Holocene climate and sea‐ice reconstructions from Icelandic marine and lake sediments

Geochemical, mineralogical and biological indicators preserved in sediments are widely used to reconstruct past climate change, but proxies differ in the degree to which their utility as climate indicators has been validated via laboratory experiments, modern spatial calibrations, or down‐core comparisons with instrumental climate data. Multi‐proxy studies provide another means of evaluating interpretations of proxies. This paper presents a multi‐proxy assessment comparing 19 sub‐centennially resolved late Holocene proxy records, covering the period 300–1900 AD, from seven Icelandic marine and lacustrine core sites. We employ simple statistical comparisons between proxy reconstructions to evaluate their correlations over time and, ultimately, their utility as proxies for regional climate. Proxies examined include oxygen isotopic composition of benthic and planktonic foraminifera, abundance of the sea‐ice biomarker IP₂₅, allochthonous quartz in marine sediments (a proxy for drift ice around Iceland), marine carbonate abundance, total organic carbon concentration, chironomid assemblages, lacustrine biogenic silica and carbon/nitrogen ratios in lake sediments. Most of the examined proxy records, including temperature and sea‐ice proxies, correlate strongly with each other over multi‐centennial timescales, and thus do appear to record changes in regional climate. Copyright © 2011 John Wiley & Sons, Ltd.     

More accurate,calibrated bootstrap confidence intervals for estimating the correlation between two time series

Estimation of Pearson’s correlation coefficient between two time series, in the evaluation of the influences of one time-dependent variable on another, is an often used statistical method in climate sciences. Data properties common to climate time series, namely non-normal distributional shape, serial correlation, and small data sizes, call for advanced, robust methods to estimate accurate confidence intervals to support the correlation point estimate. Bootstrap confidence intervals are estimated in the Fortran 90 program PearsonT (Mudelsee, Math Geol 35(6):651–665, 2003), where the main intention is to obtain accurate confidence intervals for correlation coefficients between two time series by taking the serial dependence of the data-generating process into account. However, Monte Carlo experiments show that the coverage accuracy of the confidence intervals for smaller data sizes can be substantially improved. In the present paper, the existing program is adapted into a new version, called PearsonT3, by calibrating the confidence interval to increase the coverage accuracy. Calibration is a bootstrap resampling technique that performs a second bootstrap loop (it resamples from the bootstrap resamples). It offers, like the non-calibrated bootstrap confidence intervals, robustness against the data distribution. Pairwise moving block bootstrap resampling is used to preserve the serial dependence of both time series. The calibration is applied to standard error-based bootstrap Student’s $t$ confidence intervals. The performance of the calibrated confidence interval is examined with Monte Carlo simulations and compared with the performance of confidence intervals without calibration. The coverage accuracy is evidently better for the calibrated confidence intervals where the coverage error is acceptably small already (i.e., within a few percentage points) for data sizes as small as 20.     

Bayesian modelling the retreat of the Irish Sea Ice Stream

We present an 8000‐year history spanning 650 km of ice margin retreat for the largest marine‐terminating ice stream draining the former British–Irish Ice Sheet. Bayesian modelling of the geochronological data shows the ISIS expanded 34.0–25.3 ka, accelerating into the Celtic Sea to reach maximum limits 25.3–24.5 ka before a collapse with rapid marginal retreat to the northern Irish Sea Basin (ISB). This retreat was rapid and driven by climatic warming, sea‐level rise, mega‐tidal amplitudes and reactivation of meridional circulation in the North Atlantic. The retreat, though rapid, is uneven, with the stepped retreat pattern possibly a function of the passage of the ice stream between normal and adverse ice bed gradients and changing ice stream geometry. Initially, wide calving margins and adverse slopes encouraged rapid retreat (～550 m a^?1) that slowed (～100 m a^?1) at the topographic constriction and bathymetric high between southern Ireland and Wales before rates increased (～200 m a^?1) across adverse bed slopes and wider and deeper basin configuration in the northern ISB. These data point to the importance of the ice bed slope and lateral extent in predicting the longer‐term (>1000 a) patterns and rates of ice‐marginal retreat during phases of rapid collapse, which has implications for the modelling of projected rapid retreat of present‐day ice streams. Copyright © 2013 John Wiley & Sons, Ltd.     

Subsurface lateral flow generation in aspen and conifer‐dominated hillslopes of a first order catchment in northern Utah

Mountain headwater catchments in the semi‐arid Intermountain West are important sources of surface water because these high elevations receive more precipitation than neighboring lowlands. This study examined subsurface runoff in two hillslopes, one aspen dominated, the other conifer dominated, adjacent to a first order stream in snow‐driven northern Utah. Snow accumulation, soil moisture, trenchflow and streamflow were examined in hillslopes and their adjacent stream. Snow water equivalents (SWEs) were greater under aspen stands compared to conifer, the difference increasing with higher annual precipitation. Semi‐variograms of shallow spatial soil moisture patterns and transects of continuous soil moisture showed no increase in soil moisture downslope, suggesting the absence of subsurface flow in shallow (～12 cm) soil layers of either vegetation type. However, a clear threshold relationship between soil moisture and streamflow indicated hillslope–stream connectivity, deeper within the soil profile. Subsurface flow was detected at ～50 cm depth, which was sustained for longer in the conifer hillslope. Soil profiles under the two vegetation types varied, with deep aspen soils having greater water storage capacity than shallow rocky conifer soils. Though SWEs were less under the conifers, the soil profile had less water storage capacity and produced more subsurface lateral flow during the spring snowmelt. Copyright © 2010 John Wiley & Sons, Ltd.     

An international intercalibration for mercury in seawater

An intercalibration between 32 laboratories in 16 countries was conducted using acidified and ultraviolet-irradiated seawater. One sample was natural seawater but two were spiked with mercuric chloride, 15.4 and 143 ngl^?1, respectively. The results show, for the majority of the laboratories, good accuracy and precision in the recovery of spikes but they also demonstrate that the most serious errors in low-level determinations are systematic.     

A multiphysics and multiscale software environment for modeling astrophysical systems

We present MUSE, a software framework for combining existing computational tools for different astrophysical domains into a single multiphysics, multiscale application. MUSE facilitates the coupling of existing codes written in different languages by providing inter-language tools and by specifying an interface between each module and the framework that represents a balance between generality and computational efficiency. This approach allows scientists to use combinations of codes to solve highly coupled problems without the need to write new codes for other domains or significantly alter their existing codes. MUSE currently incorporates the domains of stellar dynamics, stellar evolution and stellar hydrodynamics for studying generalized stellar systems. We have now reached a “Noah’s Ark” milestone, with (at least) two available numerical solvers for each domain. MUSE can treat multiscale and multiphysics systems in which the time- and size-scales are well separated, like simulating the evolution of planetary systems, small stellar associations, dense stellar clusters, galaxies and galactic nuclei. In this paper we describe three examples calculated using MUSE: the merger of two galaxies, the merger of two evolving stars, and a hybrid N-body simulation. In addition, we demonstrate an implementation of MUSE on a distributed computer which may also include special-purpose hardware, such as GRAPEs or GPUs, to accelerate computations. The current MUSE code base is publicly available as open source at http://muse.li.     

Geology and impact features of Vargeão Dome,southern Brazil

Abstract– Vargeão Dome (southern Brazil) is a circular feature formed in lava flows of the Lower Cretaceous Serra Geral Formation and in sandstones of the Paraná Basin. Even though its impact origin was already proposed in the 1980s, little information about its geological and impact features is available in the literature. The structure has a rim‐rim diameter of approximately 12 km and comprises several ring‐like concentric features with multiple concentric lineaments. The presence of a central uplift is suggested by the occurrence of deformed sandstone strata of the Botucatu and Pirambóia formations. We present the morphological/structural characteristics of Vargeão Dome, characterize the different rock types that occur in its interior, mainly brecciated volcanic rocks (BVR) of the Serra Geral Formation, and discuss the deformation and shock features in the volcanic rocks and in sandstones. These features comprise shatter cones in sandstone and basalt, as well as planar microstructures in quartz. A geochemical comparison of the target rock equivalents from outside the structure with the shocked rocks from its interior shows that both the BVRs and the brecciated sandstone have a composition largely similar to that of the corresponding unshocked lithologies. No traces of meteoritic material have been found so far. The results confirm the impact origin of Vargeão Dome, making it one of the largest among the rare impact craters in basaltic targets known on Earth.