Channelized melt flow in downwelling mantle: Implications for 226Ra-210Pb disequilibria in arc magmas

We present the results of an analytical model of porous flow of viscous melt into a steadily dilating ‘‘channel’’ (defined as a cluster of smaller veins) in downwelling subarc mantle. The model predicts the pressure drop in the mantle wedge matrix surrounding the channel needed to drive melt flow as a function of position and time. Melt is sucked toward the dilatant region at a near-constant velocity (105 s1) until veins comprising the channel stop opening (t = t). Fluid elements that complete their journey within the time span t < t arrive at a channel. Our results make it possible to calculate the region of influence sampled by melt that surrounds the channel. This region is large compared to the model size of the channelized region driving flow. For a baseline dilation time of 1 year and channel half width of 2 m, melt can be sampled over an 80-m radius and has the opportunity to sample matrix material with potentially contrasting chemistry on geologically short timescales. Our mechanical results are consistent with a downgoing arc mantle wedge source region where melting and melt extraction by porous flow to a channel network are sufficiently rapid to preserve source-derived 238U-230Th-226Ra, and potentially also 226 Ra-210Pb, disequilibria, prior to magma ascent to the surface. Since this is the rate-determining step in the overall process, it allows the possibility that such short-lived disequilibria measured in arc rocks at the surface are derived from deep in the mantle wedge. Stresses due to partial melting do not appear capable of producing the desired sucking effect, while the order of magnitude rate of shear required to drive dilation of 107 s1 is much larger than values resulting from steady state subduction. We conclude that local deformation rates in excess of background plate tectonic rates are needed to ‘‘switch on’’ the dilatant channel network and to initiate the sucking effect

Subduction factory: 4. Depth‐dependent flux of H 2 O from subducting slabs worldwide

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/95328/1/jgrb16727.pd

Grain-size distribution in the mantle wedge of subduction zones

Author Posting. © American Geophysical Union, 2011. 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 116 (2011): B10203, doi:10.1029/2011JB008294.Mineral grain size plays an important role in controlling many processes in the mantle wedge of subduction zones, including mantle flow and fluid migration. To investigate the grain-size distribution in the mantle wedge, we coupled a two-dimensional (2-D) steady state finite element thermal and mantle-flow model with a laboratory-derived grain-size evolution model. In our coupled model, the mantle wedge has a composite olivine rheology that incorporates grain-size-dependent diffusion creep and grain-size-independent dislocation creep. Our results show that all subduction settings lead to a characteristic grain-size distribution, in which grain size increases from 10 to 100 μm at the most trenchward part of the creeping region to a few centimeters in the subarc mantle. Despite the large variation in grain size, its effect on the mantle rheology and flow is very small, as >90% of the deformation in the flowing part of the creeping region is accommodated by grain-size-independent dislocation creep. The predicted grain-size distribution leads to a downdip increase in permeability by ∼5 orders of magnitude. This increase is likely to promote greater upward migration of aqueous fluids and melts where the slab reaches ∼100 km depth compared with shallower depths, potentially providing an explanation for the relatively uniform subarc slab depth. Seismic attenuation derived from the predicted grain-size distribution and thermal field is consistent with the observed seismic structure in the mantle wedge at many subduction zones, without requiring a significant contribution by the presence of melt.Funding for this research was provided by the National Science Foundation through a MARGINS Postdoctoral Fellowship (NSF OCE‐0840800) and NSF grant EAR‐0854673

Upper mantle electrical resistivity structure beneath the central Mariana subduction system

Author Posting. © American Geophysical Union, 2010. 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 11 (2010): Q09003, doi:10.1029/2010GC003101.This paper reports on a magnetotelluric (MT) survey across the central Mariana subduction system, providing a comprehensive electrical resistivity image of the upper mantle to address issues of mantle dynamics in the mantle wedge and beneath the slow back-arc spreading ridge. After calculation of MT response functions and their correction for topographic distortion, two-dimensional electrical resistivity structures were generated using an inversion algorithm with a smoothness constraint and with additional restrictions imposed by the subducting slab. The resultant isotropic electrical resistivity structure contains several key features. There is an uppermost resistive layer with a thickness of up to 150 km beneath the Pacific Ocean Basin, 80–100 km beneath the Mariana Trough, and 60 km beneath the Parece Vela Basin along with a conductive mantle beneath the resistive layer. A resistive region down to 60 km depth and a conductive region at greater depth are inferred beneath the volcanic arc in the mantle wedge. There is no evidence for a conductive feature beneath the back-arc spreading center. Sensitivity tests were applied to these features through inversion of synthetic data. The uppermost resistive layer is the cool, dry residual from the plate accretion process. Its thickness beneath the Pacific Ocean Basin is controlled mainly by temperature, whereas the roughly constant thickness beneath the Mariana Trough and beneath the Parece Vela Basin regardless of seafloor age is controlled by composition. The conductive mantle beneath the uppermost resistive layer requires hydration of olivine and/or melting of the mantle. The resistive region beneath the volcanic arc down to 60 km suggests that fluids such as melt or free water are not well connected or are highly three-dimensional and of limited size. In contrast, the conductive region beneath the volcanic arc below 60 km depth reflects melting and hydration driven by water release from the subducting slab. The resistive region beneath the back-arc spreading center can be explained by dry mantle with typical temperatures, suggesting that any melt present is either poorly connected or distributed discontinuously along the strike of the ridge. Evidence for electrical anisotropy in the central Mariana upper mantle is weak.Japanese participation in the Marianas experiment was supported by Japan Society for the Promotion of Science for Grant-In-Aid for Scientific Research (15340149 and 12440116), Japan-U.S. Integrated Action Program and the 21st Century COE Program of Origin and Evolution of Planetary Systems, and by the Ministry of Education, Culture, Sports, Science, and Technology for the Stagnant Slab Project, Grant-in Aid for Scientific Research on Priority Areas (17037003 and 16075204). U.S. participation was supported by NSF grant OCE0405641. Australian support came from Flinders University. T. M. is supported by the Postdoctoral Scholar Program at the Woods Hole Oceanographic Institution, with funding provided by the Deep Ocean Exploration Institute

Diagnostic value of a new sensitive membrane based technique for instantaneous D-dimer evaluation in patients with clinically suspected deep venous thrombosis.

International audienceBACKGROUND: Plasma D-Dimer analysis, using ELISA assays, has demonstrated in previous studies a high sensitivity, suggesting its utility in excluding deep venous thrombosis (DVT). Aim: To assess the performance of a new rapid plasma D-Dimer ELISA measurement in suspected DVT patients with recent clinical signs, not exceeding one week. METHODS: A prospective study of patients admitted for a suspected recent DVT. Contrast venography or compression ultrasonography were performed within 24 h of admission. A new membrane based ELISA technique, which uses an immunofiltration and two complementary monoclonal antibodies was tested. Results were expressed as positive or negative. A standard plasma D-Dimer ELISA measurement was also performed. D-Dimer performances were assessed at the end of the study. RESULTS: 265/448 patients had a proven DVT (72 distal, 193 proximal). The sensitivity of the instantaneous method in the diagnosis of overall DVT is 92 +/- 3.4% (95% CI), and specificity is 36.6 +/- 6.9%. Positive predictive value is 67.7 +/- 4.8% and negative predictive value is 76.1 +/- 8.9%. Sensitivity and negative predictive values reach 97.9 and 94.3% in the diagnosis of proximal DVT, but only 76.3 and 79.7% in the diagnosis of distal DVT. Similar results are observed with the standard ELISA method. CONCLUSION: This new rapid plasma D-Dimer measurement appears highly sensitive, and could substitute the older ELISA methods. Both methods provide lower sensitivity in the case of a distal DVT location

A community benchmark for subduction zone modeling

We have developed a suite of benchmarks to facilitate the comparison of numerical models for the dynamics and thermal structure of subduction zones. The benchmark cases are based on a thermomechanical approach in which the slab is prescribed kinematically and the wedge flow is computed dynamically. We propose various cases to investigate the influence of boundary conditions and rheology on wedge flow and resulting thermal structure. A comparison between the codes suggest that accurate modeling of the thermal field requires a good implementation of the velocity discontinuity along the seismogenic zone and high resolution in the thermal boundary layers. A minor modification to the boundary conditions of the wedge flow is also necessary to avoid a pressure singularity that exists in analytical solutions of the cornerflow model. © 2008 Elsevier B.V. All rights reserved

Fluid-assisted deformation of the subduction interface: Coupled and decoupled regimes from 2-D hydromechanical modeling

Shear deformation, accompanied with fluid activity inside the subduction interface, is related to many tectonic energy-releasing events, including regular and slow earthquakes. We have numerically examined the fluid-rock interactions inside a deforming subduction interface using state-of-the-art 2-D hydromechanical numerical models, which incorporate the rock fracturing behavior as a plastic rheology which is dependent on the pore fluid pressure. Our modeling results suggest that two typical dynamical regimes of the deforming subduction interface exist, namely, a “coupled” and a “decoupled” regime. In the coupled regime the subduction interface is subdivided into multiple rigid blocks, each separated by a narrow shear zone inclined at an angle of 15–20° with respect to the slab surface. In contrast, in the decoupled regime the subduction interface is divided into two distinct layers moving relative to each other along a pervasive slab surface-parallel shear zone. Through a systematic parameter study, we observe that the tensile strength (cohesion) of the material within the subduction interface dictates the resulting style of deformation within the interface: high cohesion (~60 MPa) results in the coupled regime, while low cohesion (~10 MPa) leads to the decoupled regime. We also demonstrate that the lithostatic pressure and inflow/outflow fluid fluxes (i.e., fluid-fluxed boundary condition) influence the location and orientation of faults. Predictions from our numerical models are supported by experimental laboratory studies, geological data, and geophysical observations from modern subduction settings