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.

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.

Timing, origin and emplacement dynamics of mass flows offshore of SE Montserrat in the last 110 ka : implications for landslide and tsunami hazards, eruption history, and volcanic island evolution

Mass flows on volcanic islands generated by volcanic lava dome collapse and by larger-volume flank collapse can be highly dangerous locally and may generate tsunamis that threaten a wider area. It is therefore important to understand their frequency, emplacement dynamics, and relationship to volcanic eruption cycles. The best record of mass flow on volcanic islands may be found offshore, where most material is deposited and where intervening hemipelagic sediment aids dating. Here we analyze what is arguably the most comprehensive sediment core data set collected offshore from a volcanic island. The cores are located southeast of Montserrat, on which the Soufriere Hills volcano has been erupting since 1995. The cores provide a record of mass flow events during the last 110 thousand years. Older mass flow deposits differ significantly from those generated by the repeated lava dome collapses observed since 1995. The oldest mass flow deposit originated through collapse of the basaltic South Soufriere Hills at 103-110 ka, some 20-30 ka after eruptions formed this volcanic center. A ∼1.8 km3 blocky debris avalanche deposit that extends from a chute in the island shelf records a particularly deep-seated failure. It likely formed from a collapse of almost equal amounts of volcanic edifice and coeval carbonate shelf, emplacing a mixed bioclastic-andesitic turbidite in a complex series of stages. This study illustrates how volcanic island growth and collapse involved extensive, large-volume submarine mass flows with highly variable composition. Runout turbidites indicate that mass flows are emplaced either in multiple stages or as single events.

Late Pleistocene stratigraphy of IODP Site U1396 and compiled chronology offshore of south and south west Montserrat, Lesser Antilles

Marine sediments around volcanic islands contain an archive of volcaniclastic deposits, which can be used to reconstruct the volcanic history of an area. Such records hold many advantages over often incomplete terrestrial datasets. This includes the potential for precise and continuous dating of intervening sediment packages, which allow a correlatable and temporally-constrained stratigraphic framework to be constructed across multiple marine sediment cores. Here, we discuss a marine record of eruptive and mass-wasting events spanning ~250 ka offshore of Montserrat, using new data from IODP Expedition 340, as well as previously collected cores. By using a combination of high-resolution oxygen isotope stratigraphy, AMS radiocarbon dating, biostratigraphy of foraminifera and calcareous nannofossils and clast componentry, we identify five major events at Soufriere Hills volcano since 250 ka. Lateral correlation of these events across sediment cores collected offshore of the south and south west of Montserrat, have improved our understanding of the timing, extent and associations between events in this area. Correlations reveal that powerful and potentially erosive density-currents travelled at least 33 km offshore, and demonstrate that marine deposits, produced by eruption-fed and mass-wasting events on volcanic islands, are heterogeneous in their spatial distribution. Thus, multiple drilling/coring sites are needed to reconstruct the full chronostratigraphy of volcanic islands. This multidisciplinary study will be vital to interpreting the chaotic records of submarine landslides at other sites drilled during Expedition 340 and provides a framework that can be applied to the stratigraphic analysis of sediments surrounding other volcanic islands.

Submarine record of volcanic island construction and collapse in the Lesser Antilles arc: First scientific drilling of submarine volcanic island landslides by IODP Expedition 340

IODP Expedition 340 successfully drilled a series of sites offshore Montserrat, Martinique and Dominica in the Lesser Antilles from March to April 2012. These are among the few drill sites gathered around volcanic islands, and the first scientific drilling of large and likely tsunamigenic volcanic island-arc landslide deposits. These cores provide evidence and tests of previous hypotheses for the composition and origin of those deposits. Sites U1394, U1399, and U1400 that penetrated landslide deposits recovered exclusively seafloor sediment, comprising mainly turbidites and hemipelagic deposits, and lacked debris avalanche deposits. This supports the concepts that i/ volcanic debris avalanches tend to stop at the slope break, and ii/ widespread and voluminous failures of preexisting low-gradient seafloor sediment can be triggered by initial emplacement of material from the volcano. Offshore Martinique (U1399 and 1400), the landslide deposits comprised blocks of parallel strata that were tilted or microfaulted, sometimes separated by intervals of homogenized sediment (intense shearing), while Site U1394 offshore Montserrat penetrated a flat-lying block of intact strata. The most likely mechanism for generating these large-scale seafloor sediment failures appears to be propagation of a decollement from proximal areas loaded and incised by a volcanic debris avalanche. These results have implications for the magnitude of tsunami generation. Under some conditions, volcanic island landslide deposits composed of mainly seafloor sediment will tend to form smaller magnitude tsunamis than equivalent volumes of subaerial block-rich mass flows rapidly entering water. Expedition 340 also successfully drilled sites to access the undisturbed record of eruption fallout layers intercalated with marine sediment which provide an outstanding high-resolution data set to analyze eruption and landslides cycles, improve understanding of magmatic evolution as well as offshore sedimentation processes.

Numerical modeling of the contaminant Transport through stratified soil

The transport processes of the dissolved chemicals in stratified or layered soils have been studied for several decades. In case of the solute transport through stratified layers, interface condition plays an important role in determining appropriate transport parameters. First‐ type and third‐ type interface conditions are generally used in the literature. A first‐type interface condition will result in a continuous concentration profile across the interface at the expense of solute mass balance. On the other hand, a discontinuity in concentration develops when a third‐ type interface condition is used. To overcome this problem, a combined first‐ and third‐ type condition at the interface has been widely employed which yields second‐ type condition. This results in a similar break‐through curve irrespective of the layering order, which is non‐physical. In this work, an interface condition is proposed which satisfies the mass balance implicitly and brings the distinction between the breakthrough curves for different layering sequence corroborating with the experimental observations. This is in disagreement with the earlier work by H. M. Selim and co‐workers but, well agreement with the hypothetical result by Bosma and van der Zee; and Van der Zee.

A comparative study of ordinary kriging and support vector machine models for the spatial variability of rock depth in Bangalore

In this paper, reduced level of rock at Bangalore, India is arrived from the 652 boreholes data in the area covering 220 sq.km. In the context of prediction of reduced level of rock in the subsurface of Bangalore and to study the spatial variability of the rock depth, ordinary kriging and Support Vector Machine (SVM) models have been developed. In ordinary kriging, the knowledge of the semivariogram of the reduced level of rock from 652 points in Bangalore is used to predict the reduced level of rock at any point in the subsurface of Bangalore, where field measurements are not available. A cross validation (Q1 and Q2) analysis is also done for the developed ordinary kriging model. The SVM is a novel type of learning machine based on statistical learning theory, uses regression technique by introducing e-insensitive loss function has been used to predict the reduced level of rock from a large set of data. A comparison between ordinary kriging and SVM model demonstrates that the SVM is superior to ordinary kriging in predicting rock depth.

Interlaboratory comparison of magnesium isotopic compositions of 12 felsic to ultramafic igneous rock standards analyzed by MC-ICPMS

To evaluate the interlaboratory mass bias for high-precision stable Mg isotopic analysis of natural materials, a suite of silicate standards ranging in composition from felsic to ultramafic were analyzed in five laboratories by using three types of multicollector inductively coupled plasma mass spectrometer (MC-ICPMS). Magnesium isotopic compositions from all labs are in agreement for most rocks within quoted uncertainties but are significantly (up to 0.3 parts per thousand in Mg-26/Mg-24, > 4 times of uncertainties) different for some mafic samples. The interlaboratory mass bias does not correlate with matrix element/Mg ratios, and the mechanism for producing it is uncertain but very likely arises from column chemistry. Our results suggest that standards with different matrices are needed to calibrate the efficiency of column chemistry and caution should be taken when dealing with samples with complicated matrices. Well-calibrated standards with matrix elements matching samples should be used to reduce the interlaboratory mass bias.

Condicionantes Geomorfológicos e Análise Geoecológica em Unidades de Conservação: Estudo de caso no Parque Natural Municipal Montanhas de Teresópolis

A emergência da temática ambiental no cenário nacional, com repercussões locais, levou a criação de inúmeras Unidades de Conservação nas últimas décadas, com o intuito de preservar seu espaço territorial e seus recursos ambientais. Entretanto, ao se delimitar uma área de conservação, é necessária a compreensão acerca das relações homemmeio, visto que ambos estão integrados dentro do mesmo sistema ambiental, e por isso, suscetíveis a mesma dinâmica de funcionamento. Neste contexto, este estudo tem como objetivo realizar uma análise geoecológica, com ênfase na geomorfologia, da área do Parque Natural Municipal Montanhas de Teresópolis. Para isto, serão analisados alguns parâmetros morfométricos das bacias hidrográficas da região na qual o parque está inserido, tal como o seu entorno. Localizado na região serrana do estado do Rio de Janeiro, durante a delimitação da área deste parque, separouse o homem da natureza, desconsiderando a integração de ambos em um só sistema. Nesse sentido, a partir de um histórico sobre a visão do homem em relação a natureza, pautamos nossa analise nas perspectivas dos Geossistemas e da Geoecologia da Paisagem. Considerando o sistema como o conjunto de elementos que se encontram em relação entre si, e que formam uma determinada integridade, a forma como este processo ocorreu evidencia as consequências de distúrbios no equilíbrio dos seus componentes, quando estes sofrem interferências humanas de maneira equivocada. Esse desequilíbrio no sistema ambiental no qual o PNMMT está inserido ficou evidente durante os eventos catastróficos que ocorreram na região em janeiro de 2011, no qual mais de 390 pessoas morreram apenas em Teresópolis. Nesse evento, o sistema ambiental da área do parque foi atingido, tanto áreas protegidas quanto áreas de ocupação urbana, negligenciadas por politicas habitacionais ao longo da historia do município, e também, durante a delimitação da área do parque, Evidencia-se assim, o quão prejudicial podem ser as intervenções humanas, quando realizadas de maneira equivocada. Através da presente análise, pretendemos diagnosticar possíveis falhas durante a delimitação desta área. Por meio das análises morfométricas, realizadas com o auxílio de ferramentas do geoprocessamento, observamos que a maior parte da área do PNMMT está inserida no domínio montanhoso e, a medida que há uma atuação deste em conjunto com outros fatores hidrológicos, a área como um todo, assim como sua zona de amortecimento, torna-se mais susceptível à processos erosivos. Deste modo, serão sugeridas medidas que busquem tanto uma conservação mais efetiva dos remanescentes da Mata Atlântica da região, quanto locais menos susceptíveis a eventos tais como os de 2011, para a implantação de politicas habitacionais. Uma vez que, diante dos resultados encontrados, fica evidente que existem uma área maior a ser conservada e, que esta não é compatível com moradias.

Mechanical and hydraulic properties of Nankai accretionary prism sediments: Effect of stress path

We have conducted triaxial deformation experiments along different loading paths on prism sediments from the Nankai Trough. Different load paths of isotropic loading, uniaxial strain loading, triaxial compression (at constant confining pressure, Pc), undrained Pc reduction, drained Pc reduction, and triaxial unloading at constant Pc, were used to understand the evolution of mechanical and hydraulic properties under complicated stress states and loading histories in accretionary subduction zones. Five deformation experiments were conducted on three sediment core samples for the Nankai prism, specifically from older accreted sediments at the forearc basin, underthrust slope sediments beneath the megasplay fault, and overthrust Upper Shikoku Basin sediments along the frontal thrust. Yield envelopes for each sample were constructed based on the stress paths of Pc-reduction using the modified Cam-clay model, and in situ stress states of the prism were constrained using the results from the other load paths and accounting for horizontal stress. Results suggest that the sediments in the vicinity of the megasplay fault and frontal thrust are highly overconsolidated, and thus likely to deform brittle rather than ductile. The porosity of sediments decreases as the yield envelope expands, while the reduction in permeability mainly depends on the effective mean stress before yield, and the differential stress after yield. An improved understanding of sediment yield strength and hydromechanical properties along different load paths is necessary to treat accurately the coupling of deformation and fluid flow in accretionary subduction zones. © 2012 American Geophysical Union All Rights Reserved.