Tracking the magmatic evolution of island arc volcanism : insights from a high-precision Pb isotope record of Montserrat, Lesser Antilles

The volcanic succession on Montserrat provides an opportunity to examine the magmatic evolution of island arc volcanism over a ∼2.5 Ma period, extending from the andesites of the Silver Hills center, to the currently active Soufrière Hills volcano (February 2010). Here we present high-precision double-spike Pb isotope data, combined with trace element and Sr-Nd isotope data throughout this period of Montserrat's volcanic evolution. We demonstrate that each volcanic center; South Soufrière Hills, Soufrière Hills, Centre Hills and Silver Hills, can be clearly discriminated using trace element and isotopic parameters. Variations in these parameters suggest there have been systematic and episodic changes in the subduction input. The SSH center, in particular, has a greater slab fluid signature, as indicated by low Ce/Pb, but less sediment addition than the other volcanic centers, which have higher Th/Ce. Pb isotope data from Montserrat fall along two trends, the Silver Hills, Centre Hills and Soufrière Hills lie on a general trend of the Lesser Antilles volcanics, whereas SSH volcanics define a separate trend. The Soufrière Hills and SSH volcanic centers were erupted at approximately the same time, but retain distinctive isotopic signatures, suggesting that the SSH magmas have a different source to the other volcanic centers. We hypothesize that this rapid magmatic source change is controlled by the regional transtensional regime, which allowed the SSH magma to be extracted from a shallower source. The Pb isotopes indicate an interplay between subduction derived components and a MORB-like mantle wedge influenced by a Galapagos plume-like source.

Pore formation during dehydration of polycrystalline gypsum observed and quantified in a time-series synchrotron radiation based X-ray micro-tomography experiment

We conducted an in-situ X-ray micro-computed tomography heating experiment at the Advanced Photon Source (USA) to dehydrate an unconfined 2.3 mm diameter cylinder of Volterra Gypsum. We used a purpose-built X-ray transparent furnace to heat the sample to 388 K for a total of 310 min to acquire a three-dimensional time-series tomography dataset comprising nine time steps. The voxel size of 2.2 μm3 proved sufficient to pinpoint reaction initiation and the organization of drainage architecture in space and time. We observed that dehydration commences across a narrow front, which propagates from the margins to the centre of the sample in more than four hours. The advance of this front can be fitted with a square-root function, implying that the initiation of the reaction in the sample can be described as a diffusion process. Novel parallelized computer codes allow quantifying the geometry of the porosity and the drainage architecture from the very large tomographic datasets (20483 voxels) in unprecedented detail. We determined position, volume, shape and orientation of each resolvable pore and tracked these properties over the duration of the experiment. We found that the pore-size distribution follows a power law. Pores tend to be anisotropic but rarely crack-shaped and have a preferred orientation, likely controlled by a pre-existing fabric in the sample. With on-going dehydration, pores coalesce into a single interconnected pore cluster that is connected to the surface of the sample cylinder and provides an effective drainage pathway. Our observations can be summarized in a model in which gypsum is stabilized by thermal expansion stresses and locally increased pore fluid pressures until the dehydration front approaches to within about 100 μm. Then, the internal stresses are released and dehydration happens efficiently, resulting in new pore space. Pressure release, the production of pores and the advance of the front are coupled in a feedback loop.

Micas in experimentally shocked gneiss

Biotites and muscovites from a gneiss have been experimentally shocked between 18 and 70 GPa using powder-propellant guns at NASA Johnson Space Center and at the California Institute of Technology. This study shows that shock in biotite and muscovite can produce homogeneous and devolatilized glasses within microseconds. Shock-deformed micas display fracturing, kinking, and complex extinction patterns over the entire pressure range investigated. However, these deformation features are not a sensitive pressure indicator. Localized melting of micas begins at 33 GPa and goes to completion at 70 GPa. Melted biotite and muscovite are optically opaque, but show extensive microvesiculation and flow when observed with the SEM. Electron diffraction confirms that biotite and muscovite have transformed to a glass. The distribution of vesicles in shock-vitrified mica shows escape of volatiles within the short duration of the shock experiment. Experimentally shocked biotite and muscovite undergo congruent melting. Compositions of the glasses are similar to the unshocked micas except for volatiles (H2O loss and K loss). These unusual glasses derived from mica may be quenched by rapid cooling conditions during the shock experiment. Based on these results, the extremely low H2O content of tektites may be reconciled with a terrestrial origin by impact. Release of volatiles in shock-melted micas affects the melting behavior of coexisting dry silicates during the short duration of the shock experiment. Transportation and escape of volatiles released from shock-melted micas may provide plausible mechanisms for the origin of protoatmospheres on terrestrial planets, hydrothermal activity on phyllosilicate-rich meteorite parent bodies, and fluid entrapment in meteorites.

Mineralogy of chondritic interplanetary dust particles

From a mineralogical survey of approximately 30 chondritic micrometeorites collected from the lower stratosphere and studied in detail using current electron microscopy techniques, it is concluded that these particles represent a unique group of extraterrestrial materials. These micrometeorites differ significantly in form and texture from components of carbonaceous chondrites and contain some mineral assemblages which do not occur in any meteorite class. Electron microscope investigations of chondritic micrometeorites have established that these materials (1) are extraterrestrial in origin, (2) existed in space as small objects, (3) endured minimal alteration by planetary processes since formation, and (4) can suffer minimal pulse heating (<600°C) on entering earth's atmosphere. The probable sources for chondritic interplanetary dust particles (IDPs) are cometary and asteroidal debris and, perhaps to a lesser extent, interstellar regions. These sources have not been conclusively linked to any specific mineralogical subset of IDP, although the chondritic porous (CP) aggregate is considered of likely cometary origin. Chondritic IDPs occur in two predominant mineral assemblages: (1) carbonaceous phases and phyllosilicates and (2) carbonaceous phases and nesosilicates or inosilicates, although particles with both types of silicate assemblages are observed. Olivines, pyroxenes, layer silicates, and carbon-rich phases are the most commonly occurring minerals in many chondritic IDPs. Other phases often observed in variable proportions include sulphides, spinels, metals, metal carbides, carbonates, and minor amounts of sulphates and phosphates. Individual mineral grain sizes range from micrometers (primarily pyroxenes and olivines) to nanometers, with the predominant size for all phases less than 100 nm. Specific mineral characteristics for particular chondritic IDPs provide an indication of processes which may have occurred prior to collection in the earth's stratosphere. For example, pyroxene mineralogy in some chondritic aggregates is consistent with condensation from a vapor phase and, we consider, with condensation in a turbulent solar nebula at relatively low temperatures (<1000°C). Carbonaceous phases present in other CP aggregates have been used to imply low-temperature formation processes such as Fischer-Tropsch synthesis (∼530°C) or carbonization and graphitization (∼315°C). Alteration processes have been implicated in the formation of some layer silicates in CP aggregates and may have involved hydrocryogenic alteration at <0°C. In general, interpretations of transformation processes on submicrometer-size minerals in chondritic IDPs are consistent with formation at a radius equivalent to the asteroid belt or greater during the later stages of solar nebula evolution using currently available models.

Preliminary study on the stability of highly sensitive fiber Bragg grating temperature sensors [对高灵敏光纤光栅温度传感器的稳定性的初探]

As a novel sensitive element and due to its advantages of immunity to electrical interference, distributed measurement, etc., fiber Bragg grating (FBG) has been researched widely. To realize the substitution of high accurate electronic temperature sensors, high sensitive FBG temperature sensors can be made by taking advantage of its characters of being sensitive to both temperature and strain. Although there are reports about high sensitive FBG temperature sensors, however, few about their stability have been done. We manufactured a high sensitive FBG temperature sensor, and put it together with an average FBG temperature sensor and an electronic crystal temperature sensor into a stainless steel container filled by water to observe the room temperature change. By comparing their results in two weeks, we have found out that: although the high sensitive FBG temperature sensor is in much better agreement with the electronic crystal sensor than the average FBG sensor is, it has occurred some small drifts. Because the drifts appeared in the process of further pulling the FBG, it might be a result of the slip of the FBG fixing points. This contributes some good experiences to the application of FBG in high accuracy temperature measurement.

The application of fiber Bragg grating sensor to high precision temperature measurement [基于光纤光栅的高精度测温传感器研究]

Basing on the character that Fiber Bragg Grating (FBG) is sensitive to both temperature and strain, by using Al and Fe-Ni alloy’s bimetal structure, we successfully design and manufacture a high accuracy FBG temperature sensor for earthquake premonition. Furthermore, we analyze the accuracy of the FBG sensors with enhanced sensitivity for the first time, and get its accuracy is up to ±0.05℃ with highest resolution ever in all FBG temperature sensors (0.0014℃/pm). This work experimentally proves the feasibility of using FBG in the earthquake premonition monitoring, and builds the foundation for the application of optic technology in earthquake premonition monitoring.

Characterization of tight gas reservoir pore structure using USANS/SANS and gas adsorption analysis

Small-angle and ultra-small-angle neutron scattering (SANS and USANS) measurements were performed on samples from the Triassic Montney tight gas reservoir in Western Canada in order to determine the applicability of these techniques for characterizing the full pore size spectrum and to gain insight into the nature of the pore structure and its control on permeability. The subject tight gas reservoir consists of a finely laminated siltstone sequence; extensive cementation and moderate clay content are the primary causes of low permeability. SANS/USANS experiments run at ambient pressure and temperature conditions on lithologically-diverse sub-samples of three core plugs demonstrated that a broad pore size distribution could be interpreted from the data. Two interpretation methods were used to evaluate total porosity, pore size distribution and surface area and the results were compared to independent estimates derived from helium porosimetry (connected porosity) and low-pressure N2 and CO2 adsorption (accessible surface area and pore size distribution). The pore structure of the three samples as interpreted from SANS/USANS is fairly uniform, with small differences in the small-pore range (<2000 Å), possibly related to differences in degree of cementation, and mineralogy, in particular clay content. Total porosity interpreted from USANS/SANS is similar to (but systematically higher than) helium porosities measured on the whole core plug. Both methods were used to estimate the percentage of open porosity expressed here as a ratio of connected porosity, as established from helium adsorption, to the total porosity, as estimated from SANS/USANS techniques. Open porosity appears to control permeability (determined using pressure and pulse-decay techniques), with the highest permeability sample also having the highest percentage of open porosity. Surface area, as calculated from low-pressure N2 and CO2 adsorption, is significantly less than surface area estimates from SANS/USANS, which is due in part to limited accessibility of the gases to all pores. The similarity between N2 and CO2-accessible surface area suggests an absence of microporosity in these samples, which is in agreement with SANS analysis. A core gamma ray profile run on the same core from which the core plug samples were taken correlates to profile permeability measurements run on the slabbed core. This correlation is related to clay content, which possibly controls the percentage of open porosity. Continued study of these effects will prove useful in log-core calibration efforts for tight gas.

A USANS/SANS study of the accessibility of pores in the barnett shale to methane and water

Shale is an increasingly important source of natural gas in the United States. The gas is held in fine pores that need to be accessed by horizontal drilling and hydrofracturing techniques. Understanding the nature of the pores may provide clues to making gas extraction more efficient. We have investigated two Mississippian Barnett Shale samples, combining small-angle neutron scattering (SANS) and ultrasmall-angle neutron scattering (USANS) to determine the pore size distribution of the shale over the size range 10 nm to 10 μm. By adding deuterated methane (CD4) and, separately, deuterated water (D2O) to the shale, we have identified the fraction of pores that are accessible to these compounds over this size range. The total pore size distribution is essentially identical for the two samples. At pore sizes >250 nm, >85% of the pores in both samples are accessible to both CD4 and D2O. However, differences in accessibility to CD4 are observed in the smaller pore sizes (∼25 nm). In one sample, CD4 penetrated the smallest pores as effectively as it did the larger ones. In the other sample, less than 70% of the smallest pores (<25 nm) were accessible to CD4, but they were still largely penetrable by water, suggesting that small-scale heterogeneities in methane accessibility occur in the shale samples even though the total porosity does not differ. An additional study investigating the dependence of scattered intensity with pressure of CD4 allows for an accurate estimation of the pressure at which the scattered intensity is at a minimum. This study provides information about the composition of the material immediately surrounding the pores. Most of the accessible (open) pores in the 25 nm size range can be associated with either mineral matter or high reflectance organic material. However, a complementary scanning electron microscopy investigation shows that most of the pores in these shale samples are contained in the organic components. The neutron scattering results indicate that the pores are not equally proportioned in the different constituents within the shale. There is some indication from the SANS results that the composition of the pore-containing material varies with pore size; the pore size distribution associated with mineral matter is different from that associated with organic phases.

Permeability estimation conditioned to geophysical downhole log data in sandstones of the northern Galilee Basin, Queensland : methods and application

This study uses borehole geophysical log data of sonic velocity and electrical resistivity to estimate permeability in sandstones in the northern Galilee Basin, Queensland. The prior estimates of permeability are calculated according to the deterministic log–log linear empirical correlations between electrical resistivity and measured permeability. Both negative and positive relationships are influenced by the clay content. The prior estimates of permeability are updated in a Bayesian framework for three boreholes using both the cokriging (CK) method and a normal linear regression (NLR) approach to infer the likelihood function. The results show that the mean permeability estimated from the CK-based Bayesian method is in better agreement with the measured permeability when a fairly apparent linear relationship exists between the logarithm of permeability and sonic velocity. In contrast, the NLR-based Bayesian approach gives better estimates of permeability for boreholes where no linear relationship exists between logarithm permeability and sonic velocity.

Thermal-elastic stresses and the criticality of the continental crust

Heating or cooling can lead to high stresses in rocks due to the different thermal-elastic properties of minerals. In the upper 4 km of the crust, such internal stresses might cause fracturing. Yet it is unclear if thermal elasticity contributes significantly to critical stresses and failure deeper in Earth's continental crust, where ductile creep causes stress relaxation. We combined a heating experiment conducted in a Synchrotron microtomograph (Advanced Photon Source, USA) with numerical simulations to calculate the grain-scale stress field in granite generated by slow burial. We find that deviatoric stresses >100 MPa can be stored during burial, with relaxation times from 100's to 1000's ka, even in the ductile crust. Hence, grain-scale thermal-elastic stresses may serve as nuclei for instabilities, thus rendering the continental crust close to criticality.

Improved estimates of percolation and anisotropic permeability from 3-D X-ray microtomography using stochastic analyses and visualization

X-ray microtomography (micro-CT) with micron resolution enables new ways of characterizing microstructures and opens pathways for forward calculations of multiscale rock properties. A quantitative characterization of the microstructure is the first step in this challenge. We developed a new approach to extract scale-dependent characteristics of porosity, percolation, and anisotropic permeability from 3-D microstructural models of rocks. The Hoshen-Kopelman algorithm of percolation theory is employed for a standard percolation analysis. The anisotropy of permeability is calculated by means of the star volume distribution approach. The local porosity distribution and local percolation probability are obtained by using the local porosity theory. Additionally, the local anisotropy distribution is defined and analyzed through two empirical probability density functions, the isotropy index and the elongation index. For such a high-resolution data set, the typical data sizes of the CT images are on the order of gigabytes to tens of gigabytes; thus an extremely large number of calculations are required. To resolve this large memory problem parallelization in OpenMP was used to optimally harness the shared memory infrastructure on cache coherent Non-Uniform Memory Access architecture machines such as the iVEC SGI Altix 3700Bx2 Supercomputer. We see adequate visualization of the results as an important element in this first pioneering study.