Implications of grain size evolution on the seismic structure of the oceanic upper mantle

Author Posting. © Elsevier B.V., 2009. 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 282 (2009): 178-189, doi:10.1016/j.epsl.2009.03.014.We construct a 1-D steady-state channel flow model for grain size evolution in the oceanic upper mantle using a composite diffusion-dislocation creep rheology. Grain size evolution is calculated assuming that grain size is controlled by a competition between dynamic recrystallization and grain growth. Applying this grain size evolution model to the oceanic upper mantle we calculate grain size as a function of depth, seafloor age, and mantle water content. The resulting grain size structure is used to predict shear wave velocity (VS) and seismic quality factor (Q). For a plate age of 60 Myr and an olivine water content of 1000 H/106Si, we find that grain size reaches a minimum of ~15 mm at ~150 km depth and then increases to ~20–30 mm at a depth of 400 km. This grain size structure produces a good fit to the low seismic shear wave velocity zone (LVZ) in oceanic upper mantle observed by surface wave studies assuming that the influence of hydrogen on anelastic behavior is similar to that observed for steady state creep. Further it predicts a viscosity of ~1019 Pa s at 150 km depth and dislocation creep to be the dominant deformation mechanism throughout the oceanic upper mantle, consistent with geophysical observations. We predict larger grain sizes than proposed in recent studies, in which the LVZ was explained by a dry mantle and a minimum grain size of 1 mm. However, we show that for a 1 mm grain size, diffusion creep is the dominant deformation mechanism above 100– 200 km depth, inconsistent with abundant observations of seismic anisotropy from surface wave studies. We therefore conclude that a combination of grain size evolution and a hydrated upper mantle is the most likely explanation for both the isotropic and anisotropic seismic structure of the oceanic upper mantle. Our results also suggest that melt extraction from the mantle will be significantly more efficient than predicted in previous modeling studies that assumed grain sizes of ~1 mm.Funding for this research was provided by NSF Grants EAR-06-52707 and EAR-07-38880

Phase diagram of the random field Ising model on the Bethe lattice

The phase diagram of the random field Ising model on the Bethe lattice with a symmetric dichotomous random field is closely investigated with respect to the transition between the ferromagnetic and paramagnetic regime. Refining arguments of Bleher, Ruiz and Zagrebnov [J. Stat. Phys. 93, 33 (1998)] an exact upper bound for the existence of a unique paramagnetic phase is found which considerably improves the earlier results. Several numerical estimates of transition lines between a ferromagnetic and a paramagnetic regime are presented. The obtained results do not coincide with a lower bound for the onset of ferromagnetism proposed by Bruinsma [Phys. Rev. B 30, 289 (1984)]. If the latter one proves correct this would hint to a region of coexistence of stable ferromagnetic phases and a stable paramagnetic phase.Comment: Article has been condensed and reorganized; Figs 3,5,6 merged; Fig 4 omitted; Some discussion added at end of Sec. III; 9 pages, 5 figs, RevTeX4, AMSTe

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

La culture du vanillier

We use observations of ice sheet surface motion from a Global Positioning System network operating from 2006 to 2014 around North Lake in west Greenland to investigate the dynamical response of the Greenland Ice Sheet's ablation area to interannual variability in surface melting. We find no statistically significant relationship between runoff season characteristics and ice flow velocities within a given year or season. Over the 7 year time series, annual velocities at North Lake decrease at an average rate of −0.9 ± 1.1 m yr−2, consistent with the negative trend in annual velocities observed in neighboring regions over recent decades. We find that net runoff integrated over several preceding years has a negative correlation with annual velocities, similar to findings from the two other available decadal records of ice velocity in western Greenland. However, we argue that this correlation is not necessarily evidence for a direct hydrologic mechanism acting on the timescale of multiple years but could be a statistical construct. Finally, we stress that neither the decadal slowdown trend nor the negative correlation between velocity and integrated runoff is predicted by current ice-sheet models, underscoring that these models do not yet capture all the relevant feedbacks between runoff and ice dynamics needed to predict long-term trends in ice sheet flow

Flexible aggregation in multiple attribute decision making: Application to the Kuranda Range Road Upgrade

The conventional method of aggregating the satisfaction of transport projects with respect to multiple attributes is commonly some variant of Simple Additive Weighting (SAW), which involves the sum of products of standardized outcomes of projects with respect to attributes and attribute importance weights. It is suggested that alternative forms of aggregation might be more useful, in particular, the Ordered Weighted Averaging (OWA) operator introduced by Yager (1988). Attribute importance weights and satisfaction of attributes by projects may be aggregated prior to aggregation via an OWA operator. In this case OWA operator weights may be based on the "attitudinal character of the decision maker expressed in terms of the degree of "orness and "andness of the aggregation. A well-known approach is maximum entropy aggregation, in which weights are derived to be as "even (or as minimally dispersed) as a possible subject to satisfying a given "orness or "andness constraint. Recently, aggregation processes have been proposed by Larsen (199920022003) which have several desirable properties and also may be considered as alternative forms of aggregation. An example is given relating to the Kuranda Range Road upgrade (Queensland, Australia) which is limited by grade, poor overtaking opportunities, poor horizontal alignment, and other constraints, and the road is expected to become increasingly congested over the next few years. A more flexible Multiple Attribute Decision Making is used to identify a "best project from a set of four alternative projects

SOLPS5 modelling of ELMing H-mode on TCV

Although ohmic H-modes have long been produced on the TCV tokamak and the effects of ELMs at the divertor targets studied in some detail [1], no attempt has yet been made to model the scrape-off layer (SOL) in these plasmas. This contribution describes details of the first such efforts to do so. Simulations with the coupled fluid-Monte Carlo SOLPS5 code are constrained by careful upstream Thomson scattering and fast reciprocating Langmuir probe profiles and the results compared w ith measurements at the divertor targets. Recent experiments with high power ECRH at the third harmonic have produced large, possibly Type I ELMs on TCV for the first time, but in standard ohmic H-modes operating close to the L-H transition threshold power, only Type III ELMs are obtained. Typical single null lower H-mode discharges have I p = 400 kA, n/n GW ~ 0.3 and steady ELMing phases with f ELM ~ 100 Hz, where each ELM exhausts only a few 100 J of plasma stored energy. As such, these ELMs cannot be compared with larger events typical elsewhere regarding the magnitude of target power fluxes etc., but their behaviou r with respect to transport in the SOL and interaction with the targets is no different. In fact, their benign nature makes the fluid plasma simulation in some ways more appropriate since the ELMs are insu fficiently large to require a true kinetic simulation and are likely to be less perturbing in the sense of parallel heat flux limits and variations in sheath heat transmission coefficients. The modelling attempts described here broadly follow the approach in [2], seeking first the closest match to upstream experimental profiles during inter-ELM phases using a step-like ansatz for the perpendicular particle and heat diffusivities (D and χ i,e ) in the edge and SOL regions, but also introducing a poloidal variation of the transport coefficients both in the main chamber and dive rtor. This is extremely important in TCV, where the unconventional divertor geometry means that care must be taken in the presence of steep H-mode edge barriers to tailor differently the transport in this region compared with the core. Similar reasoning applies even more to the ELM itself, which is known to burst into the SOL in the outboard, unfavourable curvature region and is thus extremely poloidally localised. This has also been accounted for in the simulations which, as in the earlier SOLPS5 attempts to simulate Type ELMs in JET [3], m odel the ELM as an instantaneous local increase in the transport coefficients and simulate the subsequent SOL transport in a time dependent way. Good agreement is found between code and experime nt in the inter-ELM phase. Experimentally, these Type III ELMs are too rapid for comparison on an individual ELM basis. Instead, many similar events are coherently averaged and simulation comp ared with power and par ticle flux measurements at the divertor targets and line integrated observations of recycling in the divertor volume obtained using a new fast AXUV diode camera system filtered for Lyman-alpha emission