Relationship between seismic P-wave velocity and the composition of anhydrous igneous and meta-igneous rocks

Author Posting. © American Geophysical Union 2003. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 4 (2003): 1041, doi:10.1029/2002GC000393.This study presents a new approach to quantitatively assess the relationship between the composition and seismic P-wave velocity of anhydrous igneous and meta-igneous rocks. We perform thermodynamic calculations of the equilibrating phase assemblages predicted for all igneous composition space at various pressure and temperature conditions. Seismic velocities for each assemblage are then estimated from mixing theory using laboratory measurements of the elastic parameters for pure mineral phases. The resultant velocities are used to derive a direct relationship between Vp and major element composition valid to ±0.13 km/s for pressure and temperature conditions along a normal crustal geotherm in the depth range of 5–50 km and equilibration pressures ≤12 kbar. Finally, we use the calculated velocities to invert for major element chemistry as a function of P-wave velocity assuming only the in situ temperature and pressure conditions are known. Compiling typical velocity-depth profiles for the middle and lower continental and oceanic crust, we calculate compositional bounds for each of these geologic environments. We find that the acceptable compositional range for the middle (15–30 km) and lower continental (≥35 km) crust is broad, ranging from basaltic to dacitic compositions, and conclude that P-wave velocity measurements alone are insufficient to provide fundamental constraints on the composition of the middle and lower continental crust. However, because major oxides are correlated in igneous rocks, joint constraints on Vp and individual oxides can narrow the range of acceptable crustal compositions. In the case of the lower oceanic crust (≥2 km), observed velocities are 0.2–0.3 km/s lower than velocities calculated based on the average bulk composition of gabbros in drill cores and exposed ophiolite sequences. We attribute this discrepancy to a combination of residual porosity at crustal depths less than ∼10 km and hydrous alteration phases in the lower crust, and suggest caution when inferring mantle melting parameters from observed velocities in the lower oceanic crust.This research was supported by National Science Foundation Grants OCE- 9819666, EAR-9910899, and EAR-0087706 (P.B. Kelemen)

Evaluation of moisture sources for the Central European summer flood of May/June 2013 based on regional climate model simulations

Heavy precipitation affected Central Europe in May/June 2013, triggering damaging floods both on the Danube and the Elbe rivers. Based on a modelling approach with COSMO-CLM, moisture fluxes, backward trajectories, cyclone tracks and precipitation fields are evaluated for the relevant time period 30 May–2 June 2013. We identify potential moisture sources and quantify their contribution to the flood event focusing on the Danube basin through sensitivity experiments: Control simulations are performed with undisturbed ERA-Interim boundary conditions, while multiple sensitivity experiments are driven with modified evaporation characteristics over selected marine and land areas. Two relevant cyclones are identified both in reanalysis and in our simulations, which moved counter-clockwise in a retrograde path from Southeastern Europe over Eastern Europe towards the northern slopes of the Alps. The control simulations represent the synoptic evolution of the event reasonably well. The evolution of the precipitation event in the control simulations shows some differences in terms of its spatial and temporal characteristics compared to observations. The main precipitation event can be separated into two phases concerning the moisture sources. Our modelling results provide evidence that the two main sources contributing to the event were the continental evapotranspiration (moisture recycling; both phases) and the North Atlantic Ocean (first phase only). The Mediterranean Sea played only a minor role as a moisture source. This study confirms the importance of continental moisture recycling for heavy precipitation events over Central Europe during the summer half year

Diapirs as the source of the sediment signature in arc lavas

Author Posting. © The Author(s), 2011. This is the author's version of the work. It is posted here by permission of Nature Publishing Group for personal use, not for redistribution. The definitive version was published in Nature Geoscience 4 (2011): 641-646, doi:10.1038/ngeo1214Many arc lavas show evidence for the involvement of subducted sediment in the melting process. There is debate whether this “sediment melt” signature forms at relatively low temperature near the fluid-saturated solidus or at higher temperature beyond the breakdown of trace-element-rich accessory minerals. We present new geochemical data from high- to ultrahigh-pressure rocks that underwent subduction and show no significant depletion of key trace elements in the sediment melt component until peak metamorphic temperatures exceeded ~1050ºC from 2.7 to 5 GPa. These temperatures are higher than for the top of the subducting plate at similar pressures based on thermal models. To address this discrepancy, we use instability calculations for a non-Newtonian buoyant layer in a viscous half-space to show that, in typical subduction zones, solid-state sediment diapirs initiate at temperatures between 500–850ºC. Based on these calculations, we propose that the sediment melt component in arc magmas is produced by high degrees of dehydration melting in buoyant diapirs of metasediment that detach from the slab and rise into the hot mantle wedge. Efficient recycling of sediments into the wedge by this mechanism will alter volatile fluxes into the deep mantle compared to estimates based solely on devolatilization of the slab.Funding for this work was provided by NSF and WHOI’s Deep Ocean Exploration Institute

Reconstruction of the Talkeetna intraoceanic arc of Alaska through thermobarometry

Author Posting. © American Geophysical Union, 2008. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 113 (2008): B03204, doi:10.1029/2007JB005208.The Talkeetna arc is one of two intraoceanic arcs where much of the section from the upper mantle through the volcanic carapace is well exposed. We reconstruct the vertical section of the Talkeetna arc by determining the (re)crystallization pressures at various structural levels. The thermobarometry shows that the tonalites and quartz diorites intruded at ∼5–9 km into a volcanic section estimated from stratigraphy to be 7 km thick. The shallowest, Tazlina and Barnette, gabbros crystallized at ∼17–24 km; the Klanelneechena Klippe crystallized at ∼24–26 km; and the base of the arc crystallized at ∼35 km depth. The arc had a volcanic:plutonic ratio of ∼1:3–1:4. However, many or most of the felsic plutonic rocks may represent crystallized liquids rather than cumulates so that the liquid:cumulate ratio might be 1:2 or larger. The current 5- to 7-km structural thickness of the plutonic section of the arc is ∼15–30% of the original 23- to 28-km thickness. The bulk composition of the original Talkeetna arc section was ∼51–58 wt % SiO2.Funded by NSF EAR-9910899

Role of melt supply in oceanic detachment faulting and formation of megamullions

Author Posting. © Geological Society of America, 2008. This article is posted here by permission of Geological Society of America for personal use, not for redistribution. The definitive version was published in Geology 36 (2008): 455-458, doi:10.1130/G24639A.1.Normal faults are ubiquitous on mid-ocean ridges and are expected to develop increasing offset with reduced spreading rate as the proportion of tectonic extension increases. Numerous long-lived detachment faults that form megamullions with large-scale corrugations have been identified on magma-poor mid-ocean ridges, but recent studies suggest, counterintuitively, that they may be associated with elevated magmatism. We present numerical models and geological data to show that these detachments occur when ~30%–50% of total extension is accommodated by magmatic accretion and that there is significant magmatic accretion in the fault footwalls. Under these low-melt conditions, magmatism may focus unevenly along the spreading axis to create an irregular brittle-plastic transition where detachments root, thus explaining the origin of the enigmatic corrugations. Morphological and compositional characteristics of the oceanic lithosphere suggested by this study provide important new constraints to assess the distribution of magmatic versus tectonic extension along mid-ocean ridges.This research was supported by the National Science Foundation and by the Henry Bryant Bigelow Chair in Oceanography to Tucholke at Woods Hole Oceanographic Institution

Timescales for the growth of sediment diapirs in subduction zones

Author Posting. © The Author(s), 2012. This article is posted here by permission of John Wiley & Sons for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 190 (2012): 1361–1377, doi:10.1111/j.1365-246X.2012.05565.x.In this study, we calculate timescales for the growth of gravitational instabilities forming in the sediment layer on the downgoing slab at subduction zones. Subducted metasediments are buoyant with respect to the overlying mantle and may form diapirs that detach from the slab and rise upwards into the mantle wedge. We use a particle-in-cell, finite-difference method to calculate growth rates for instabilities forming within a buoyant, wet-quartz metasediment layer underlying a dense mantle half-space composed of wet olivine. These growth rates are used to determine where sediment diapirs initiate and detach from the slab over a range of subduction zone thermal structures. We find that, given a sufficient layer thickness (200–800 m, depending on slab-surface and mantle-wedge temperatures), sediment diapirs begin to grow rapidly at depths of ∼80 km and detach from the slab within 1–3 Myr at temperatures ≤900 °C and at depths roughly corresponding to the location of the slab beneath the arc. Diapir growth is most sensitive to absolute slab temperature, however it is also affected by the viscosity ratio between the sediment layer and the mantle wedge and the length-scale over which viscosity decays above the slab. These secondary affects are most pronounced in colder subduction systems with old slabs and faster subduction rates. For a broad range of subduction zone thermal conditions, we find that diapirs can efficiently transport sediments into the mantle wedge, where they would melt and be incorporated into arc magmas. Thus, we conclude that sediment diapirism is a common feature of many subduction zones, providing a potential explanation for the ‘sediment signature’ in the chemistry of arc magmas.This work was supported by NSF Grant EAR-0652707 and a WHOI Deep Ocean Exploration Institute Fellowship to MB

Analysis of reported error in Monte Carlo rendered images

Evaluating image quality in Monte Carlo rendered images is an important aspect of the rendering process as we often need to determine the relative quality between images computed using different algorithms and with varying amounts of computation. The use of a gold-standard, reference image, or ground truth (GT) is a common method to provide a baseline with which to compare experimental results. We show that if not chosen carefully the reference image can skew results leading to significant misreporting of error. We present an analysis of error in Monte Carlo rendered images and discuss practices to avoid or be aware of when designing an experiment

Mid-ocean ridge jumps associated with hotspot magmatism

Author Posting. © Elsevier B.V., 2007. This is the author's version of the work. It is posted here by permission of Elsevier B.V. for personal use, not for redistribution. The definitive version was published in Earth and Planetary Science Letters 266 (2008): 256-270, doi:10.1016/j.epsl.2007.10.055.Hotspot-ridge interaction produces a wide range of phenomena including excess crustal thickness, geochemical anomalies, off-axis volcanic ridges and ridge relocations or jumps. Ridges are recorded to have jumped toward many hotspots including, Iceland, Discovery, Galapagos, Kerguelen and Tristan de Cuhna. The causes of ridge jumps likely involve a number of interacting processes related to hotspots. One such process is reheating of the lithosphere as magma penetrates it to feed near-axis volcanism. We study this effect by using the hybrid, finite-element code, FLAC, to simulate two-dimensional (2-D, cross-section) viscous mantle flow, elasto-plastic deformation of the lithosphere and heat transport in a ridge setting near an off-axis hotspot. Heating due to magma transport through the lithosphere is implemented within a hotspot region of fixed width. To determine the conditions necessary to initiate a ridge jump, we vary four parameters: hotspot magmatic heating rate, spreading rate, seafloor age at the location of the hotspot and ridge migration rate. Our results indicate that the hotspot magmatic heating rate required to initiate a ridge jump increases non-linearly with increasing spreading rate and seafloor age. Models predict that magmatic heating, itself, is most likely to cause jumps at slow spreading rates such as at the Mid-Atlantic Ridge on Iceland. In contrast, despite the higher magma flux at the Galapagos hotspot, magmatic heating alone is probably insufficient to induce a ridge jump at the present-day due to the intermediate ridge spreading rate of the Galapagos Spreading Center. The time required to achieve a ridge jump, for fixed or migrating ridges, is found to be on the order of 105-106 years. Simulations that incorporate ridge migration predict that after a ridge jump occurs the hotspot and ridge migrate together for time periods that increase with magma flux. Model results also suggest a mechanism for ridge reorganizations not related to hotspots such as ridge jumps in back-arc settings and ridge segment propagation along the Mid-Atlantic Ridge.Mittelstaedt, Ito and Behn were funded by NSF grant OCE03-51234 and OCE05-48672

Long-range chemical sensitivity in the sulfur K-edge X-ray absorption spectra of substituted thiophenes

© 2014 American Chemical Society. Thiophenes are the simplest aromatic sulfur-containing compounds and are stable and widespread in fossil fuels. Regulation of sulfur levels in fuels and emissions has become and continues to be ever more stringent as part of governments' efforts to address negative environmental impacts of sulfur dioxide. In turn, more effective removal methods are continually being sought. In a chemical sense, thiophenes are somewhat obdurate and hence their removal from fossil fuels poses problems for the industrial chemist. Sulfur K-edge X-ray absorption spectroscopy provides key information on thiophenic components in fuels. Here we present a systematic study of the spectroscopic sensitivity to chemical modifications of the thiophene system. We conclude that while the utility of sulfur K-edge X-ray absorption spectra in understanding the chemical composition of sulfur-containing fossil fuels has already been demonstrated, care must be exercised in interpreting these spectra because the assumption of an invariant spectrum for thiophenic forms may not always be valid

Multisite Investigation of Outcomes With Implementation of CYP2C19 Genotype-Guided Antiplatelet Therapy After Percutaneous Coronary Intervention

OBJECTIVES: This multicenter pragmatic investigation assessed outcomes following clinical implementation of CYP2C19 genotype-guided antiplatelet therapy after percutaneous coronary intervention (PCI). BACKGROUND: CYP2C19 loss-of-function alleles impair clopidogrel effectiveness after PCI. METHODS: After clinical genotyping, each institution recommended alternative antiplatelet therapy (prasugrel, ticagrelor) in PCI patients with a loss-of-function allele. Major adverse cardiovascular events (defined as myocardial infarction, stroke, or death) within 12 months of PCI were compared between patients with a loss-of-function allele prescribed clopidogrel versus alternative therapy. Risk was also compared between patients without a loss-of-function allele and loss-of-function allele carriers prescribed alternative therapy. Cox regression was performed, adjusting for group differences with inverse probability of treatment weights. RESULTS: Among 1,815 patients, 572 (31.5%) had a loss-of-function allele. The risk for major adverse cardiovascular events was significantly higher in patients with a loss-of-function allele prescribed clopidogrel versus alternative therapy (23.4 vs. 8.7 per 100 patient-years; adjusted hazard ratio: 2.26; 95% confidence interval: 1.18 to 4.32; p = 0.013). Similar results were observed among 1,210 patients with acute coronary syndromes at the time of PCI (adjusted hazard ratio: 2.87; 95% confidence interval: 1.35 to 6.09; p = 0.013). There was no difference in major adverse cardiovascular events between patients without a loss-of-function allele and loss-of-function allele carriers prescribed alternative therapy (adjusted hazard ratio: 1.14; 95% confidence interval: 0.69 to 1.88; p = 0.60). CONCLUSIONS: These data from real-world observations demonstrate a higher risk for cardiovascular events in patients with a CYP2C19 loss-of-function allele if clopidogrel versus alternative therapy is prescribed. A future randomized study of genotype-guided antiplatelet therapy may be of value