A comparative study of the thermal behaviour of length-fast chalcedony, length-slow chalcedony (quartzine) and moganite

The thermal behaviour of silica rocks upon heat treatment is dependent on the constituent minerals and petrographic texture types. These constituents can be shown to be mainly quartz in the form of two types of chalcedony (Length-fast (LF) chalcedony and Length-slow (LS) chalcedony, the latter also being termed quartzine) and moganite. Even though the thermal behaviour of LF-chalcedony is well understood, major uncertainties persist concerning the high-temperature behaviour of LS-chalcedony and moganite. We present here a comparative study of these three constituents of common silica rocks. Our results show that the chemical reaction is the same in all three, Si–OH + HO–Si → Si–O–Si + H₂O, but that the reaction kinetics and activation temperatures are very different. LS-chalcedony begins to react from 200 °C upwards, that is at temperatures 50 °C below the ones observed in LF-chalcedony, and shows the fastest reaction kinetics of this ‘water’ loss. Chemically bound water (SiOH) in moganite is more stable at high temperatures and no specific activation temperature is necessary for triggering the temperature-induced ‘water’ loss. Moganite is also found to act as a stabilizer in silica rocks preventing them from temperature-induced fracturing. These findings have implications for the study of potential heat treatment temperatures of silica rocks (in industry and heritage studies), but they also shed light on the different structures of SiO₂ minerals and the role of OH impurities therein.     

Formation and Evolution of Intermediate Mass Black Hole X-Ray Binaries

The evolution of young (≲ 10 Myr) star clusters with a density exceeding about 10⁵ star pc⁻³ are strongly affected by physical stellar collisions during their early lifetime. In such environments the same star may participate in several tens to hundreds of collisions ultimately leading to the collapse of the star to a black hole of intermediate mass. At later time, the black hole may acquire a companion star by tidal capture or by dynamical – three-body – capture. When the captured star evolves it starts to fill its Roche-lobe and transfers mass to its accompanying black hole. This then leads to a bright phase of X-ray emission, which lasts for the remaining main-sequence lifetime of the donor. If the star captured by the intermediate mass black hole is relatively low mass ≲ 2 M⊙) the binary will also be visible as a bright source in gravitational waves. Based on empirical models we argue that, for as long as the donor remains on the main sequence, the source will be ultraluminous L_x >rsim 10⁴⁰ ergs^-1 for about a week every few month. When the donor star is more massive >15 M⊙, or evolved off the main sequence the bright time is longer, but the total accretion phase lasts much shorter.     

Density variations in the thickened crust as a function of pressure, temperature, and composition

Constraints on density as a function of pressure, temperature, and composition are crucial to understand isostatic movements during geodynamic processes. Here, we provide a systematic series of density diagrams extracted from thermodynamic calculations for a variety of crustal compositions within a wide P–T range. We quantify systematic density changes in collisional settings for relevant compositional variations and attempt to simplify the density–composition relationships. Rock densities depend strongly on pressure, temperature, and composition. Densities at some selected pressure–temperature conditions increase linearly with increasing Al₂O₃ as well as MgO/FeO contents in pelitic rocks. Al- and Fe-rich pelites yield the highest densities, which is mostly due to the formation of garnet but also depends on other minerals and changes of reactions. The effect of loading on densities is investigated, and we show that for deep burial, a meta-pelite rich in Fe and Mg yields much larger density changes than a dry basalt and that the burial of such a rock with a composition close to typical lower crust may result in significant negative buoyancy. Metamorphism of hydrous lower crust due to pressurization and heating thus leads to densification of thickened lower crust, while heating of dry crust leads to a decrease in density. Hence, water-loaded isostatic subsidence due to metamorphism of water-saturated lower crust is substantial and increases with the thickness and depth of the reacting layer, while dry compositions show much less or only transient densification and subsidence. The density change due to thermal expansion, an extensively used concept in geodynamic models, predicts uplift under the same P–T conditions and is an order of magnitude smaller than the density variation calculated from petrologically consistent diagrams.     

Platinum solubility and partitioning in a felsic melt-vapor-brine assemblage

The solubility and partitioning of Pt in a S-free vapor - brine - rhyolite melt - Pt metal assemblage has been quantified at 800 °C, f_O2=NNO and pressures of 100 and 140 MPa. Vapor and brine were sampled at run conditions by trapping these phases as glass-hosted fluid inclusions as the melt cooled through the glass transition temperature. The vapor and brine were in equilibrium with the melt at the time of trapping and, thus, represent fluids which were sampled at the termination of each experimental run. The microthermometrically determined salinities of vapor and brine are ∼2 and ∼63 wt.% NaCl eq. and ∼9 and ∼43 wt.% NaCl eq. at 100 and 140 MPa, respectively. Platinum solubilities in vapor, brine and glass (i.e., quenched melt) were quantified by using laser ablation - inductively coupled plasma - mass spectrometry (LA-ICP-MS). Equilibrium is discussed with reference to the major and trace element concentrations of glass-hosted fluid inclusions as well as the silicate melt over run times that varied from 110 to 377 h at 140 MPa and 159 to 564 h at 100 MPa. Platinum solubility values (±1σ) in H₂O-saturated felsic melt are 0.28 ± 0.13 μg/g and 0.38 ± 0.06 μg/g at 140 and 100 MPa, respectively. Platinum solubility values () at 140 and 100 MPa, respectively, in aqueous vapor are 0.91 ± 0.29 μg/g and 0.37 ± 0.17 μg/g and in are brine 16 ± 10 μg/g and 3.3 ± 1.0 μg/g. The measured solubility data were used to calculate Nernst-type partition coefficients for Pt between vapor/melt, brine/melt and vapor/brine. The partition coefficient values () for vapor/melt, brine/melt and vapor/brine at 140 MPa are 2.9 ± 1.0, 67 ± 27, and 0.13 ± 0.05 and at 100 MPa are 1.0 ± 0.2, 6.8 ± 2.4, and 0.15 ± 0.05. The partitioning data were used to model the Pt-scavenging capacity of vapor and brine during the crystallization-driven degassing (i.e., second boiling) of a felsic silicate melt over a depth range (i.e., 3-6 km) consistent with the evolution of magmatic-hydrothermal ore deposits. Model calculations suggest that aqueous vapor and brine can scavenge sufficient quantities of Pt, and by analogy other platinum group elements (PGE), to produce economically important PGE-rich magmatic-hydrothermal ore deposits in Earth’s upper continental crust.     

Thermal history analysis by integrated modelling of apatite fission track and vitrinite reflectance data: application to an inverted basin (Buller Coalfield, New Zealand)

Integrated analysis and modelling of apatite fission track with vitrinite reflectance (VR) data allows the timing, magnitude and pattern of Palaeogene subsidence and Neogene inversion to be established for an uplifted and largely denuded basin: the Buller Coalfield, New Zealand. At the time of maximum subsidence in the late Oligocene, the basin consisted of an extensional half graben, bounded to the west by the Kongahu Fault Zone (KFZ), with up to 6 km of upper Eocene to Oligocene section adjacent to it; currently, only a few tens of metres of basal coal measures on basement are preserved on top of a range 800–1000 m above sea level. Integrated modelling of the VR and fission track data show that the deepest parts of the basin were inverted during two Miocene compressional phases (24–19 Ma and 13–8 Ma), and are consistent with a further phase of inversion during the Quaternary that formed the present topography. Palinspastic restoration of the three phases of inversion shows that the basin was not inverted in a simple way: most of the rock uplift/denudation adjacent to the KFZ occurred during the early Miocene phase, and at the same time burial occurred in the south-eastern part of the basin (maximum temperatures were experienced at different times at different places in the basin); during the middle to late Miocene there was broad uplift in the central and eastern parts of the coalfield. Because the timing and magnitude of uplift have been derived from the zone of inversion, they can be compared independently with the timing of unconformity development and rapid subsidence in the adjacent foredeeps, particularly the Westport Trough. For the middle to late Miocene phase of inversion, we show that during the first 1–2 million years of compression, the uplift within the coalfield also involved the margins of the Westport Trough, contributing to unconformity development; subsequently, uplift continued on the inversion structure but the margins of the Westport Trough subsided rapidly. This is explained by a model of stick slip behaviour on the boundary faults, especially for the KFZ. When compression started the fault zone has locked and uplift extends into the basin, whereas subsequently the fault zone unlocks, and the inversion structure overrides the basin margin, thereby loading it and causing subsidence.     

Mineshaft imaging using surface and crosshole 3D electrical resistivity tomography: A case history from the East Pennine Coalfield, UK

Hidden mineshafts located in urban areas are a significant problem across much of the industrialized world. Electrical resistivity tomography (ERT) is a technique that can detect and characterize hidden mine entries by exploiting resistivity contrasts between the shaft and surrounding materials, resulting from either compositional or structural differences. A case study is presented in which both surface and crosshole 3D ERT surveys are used to image a hidden backfilled mineshaft at a built environment site, situated on Carboniferous Lower Coal Measures strata in the UK.Backfilled shafts generally present the greatest challenge for detection using geophysical methods, as contrasts between the fill and bedrock are typically low compared to air or water-filled conditions. Nevertheless, the shaft in this case was identified by both the surface and crosshole 3D surveys. The shaft appeared as a strongly resistive anomaly relative to background materials, which we interpreted as resulting from the disturbed fabric of the fill materials rather than any significant compositional differences. The study highlighted the respective strengths and weaknesses of the surface and crosshole ERT methodologies for this type of problem. The surface survey, which covered a non-rectangular area to accommodate irregular boundaries and other physical obstructions, provided a relatively rapid means of investigating the study site. However, this method had a limited depth of investigation and was constrained in its coverage by the locations of buildings. By contrast, the 3D crosshole method was able to image the shaft to the level of the deepest borehole electrodes. Although crosshole ERT is too expensive to be used for large-scale mineshaft surveys, this study clearly demonstrates its suitability for targeted investigations where surface methods cannot be deployed, such as scanning beneath surface structures or in situations where it is essential for resolution to be maintained with depth.     

Localized climate change scenarios of mean temperature and precipitation over Switzerland

There is a growing need of the climate change impact modeling and adaptation community to have more localized climate change scenario information available over complex topography such as in Switzerland. A gridded dataset of expected future climate change signals for seasonal averages of daily mean temperature and precipitation in Switzerland is presented. The basic scenarios are taken from the CH2011 initiative. In CH2011, a Bayesian framework was applied to obtain probabilistic scenarios for three regions within Switzerland. Here, the results for two additional Alpine sub-regions are presented. The regional estimates have then been downscaled onto a regular latitude-longitude grid with a resolution of 0.02° or roughly 2 km. The downscaling procedure is based on the spatial structure of the climate change signals as simulated by the underlying regional climate models and relies on a Kriging with external drift using height as auxiliary predictor. The considered emission scenarios are A1B, A2 and the mitigation scenario RCP3PD. The new dataset shows an expected warming of about 1 to 6 °C until the end of the 21st century, strongly depending on the scenario and the lead time. Owing to a large vertical gradient, the warming is about 1 °C stronger in the Alps than in the Swiss lowlands. In case of precipitation, the projection uncertainty is large and in most seasons precipitation can increase or decrease. In summer a distinct decrease of precipitation can be found, again strongly depending on the emission scenario.     

Estimation and propagation of volcanic source parameter uncertainty in an ash transport and dispersal model: application to the Eyjafjallajokull plume of 14–16 April 2010

Data on source conditions for the 14 April 2010 paroxysmal phase of the Eyjafjallaj?kull eruption, Iceland, have been used as inputs to a trajectory-based eruption column model, bent. This model has in turn been adapted to generate output suitable as input to the volcanic ash transport and dispersal model, puff, which was used to propagate the paroxysmal ash cloud toward and over Europe over the following days. Some of the source parameters, specifically vent radius, vent source velocity, mean grain size of ejecta, and standard deviation of ejecta grain size have been assigned probability distributions based on our lack of knowledge of exact conditions at the source. These probability distributions for the input variables have been sampled in a Monte Carlo fashion using a technique that yields what we herein call the polynomial chaos quadrature weighted estimate (PCQWE) of output parameters from the ash transport and dispersal model. The advantage of PCQWE over Monte Carlo is that since it intelligently samples the input parameter space, fewer model runs are needed to yield estimates of moments and probabilities for the output variables. At each of these sample points for the input variables, a model run is performed. Output moments and probabilities are then computed by properly summing the weighted values of the output parameters of interest. Use of a computational eruption column model coupled with known weather conditions as given by radiosonde data gathered near the vent allows us to estimate that initial mass eruption rate on 14 April 2010 may have been as high as 10⁸?kg/s and was almost certainly above 10⁷?kg/s. This estimate is consistent with the probabilistic envelope computed by PCQWE for the downwind plume. The results furthermore show that statistical moments and probabilities can be computed in a reasonable time by using 9⁴?=?6,561 PCQWE model runs as opposed to millions of model runs that might be required by standard Monte Carlo techniques. The output mean ash cloud height plus three standard deviations??encompassing c. 99.7?% of the probability mass??compares well with four-dimensional ash cloud position as retrieved from Meteosat-9 SEVIRI data for 16 April 2010 as the ash cloud drifted over north-central Europe. Finally, the ability to compute statistical moments and probabilities may allow for the better separation of science and decision-making, by making it possible for scientists to better focus on error reduction and decision makers to focus on ??drawing the line?? for risk assessment.     

Isothermal compression of an eclogite from the Western Gneiss Region (Norway)

In the Western Gneiss Region in Norway, mafic eclogites form lenses within granitoid orthogneiss and contain the best record of the pressure and temperature evolution of this ultrahigh-pressure (UHP) terrane. Their exhumation from the UHP conditions has been extensively studied, but their prograde evolution has been rarely quantified although it represents a key constraint for the tectonic history of this area. This study focused on a well-preserved phengite-bearing eclogite sample from the Nordfjord region. The sample was investigated using phase-equilibrium modelling, trace-element analyses of garnet, trace- and major-element thermobarometry and quartz-in-garnet barometry by Raman spectroscopy. Inclusions in garnet core point to crystallization conditions in the amphibolite facies at 510–600°C and 11–16 kbar, whereas chemical zoning in garnet suggests growth during isothermal compression up to the peak pressure of 28 kbar at 600°C, followed by near-isobaric heating to 660–680°C. Near-isothermal decompression to 10–14 kbar is recorded in fine-grained clinopyroxene–amphibole–plagioclase symplectites. The absence of a temperature increase during compression seems incompatible with the classic view of crystallization along a geothermal gradient in a subduction zone and may question the tectonic significance of eclogite facies metamorphism. Two end-member tectonic scenarios are proposed to explain such an isothermal compression: Either (1) the mafic rocks were originally at depth within the lower crust and were consecutively buried along the isothermal portion of the subducting slab or (2) the mafic rocks recorded up to 14 kbar of tectonic overpressure at constant depth and temperature during the collisional stage of the orogeny.