Influence of the zone of weakness on dip angle and shear heating of subducted slabs

We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated.     

Overburden deformation associated with halokinesis in the Southern North Sea: implications for the origin of the Silverpit Crater

Gravitational collapse in response to the creation of salt swells and diapirs is a common phenomenon within regions such as the Southern North Sea. Although collapse faulting and slumping of linear salt swell flanks has previously been documented, the existence of concentric collapse structures surrounding diapirs and isolated salt lows is a recently recognised feature. 3-D seismic data from the Southern North Sea demonstrates that concentric collapse faulting around both diapirs and subcircular salt lows is common. The recognition of concentric collapse structures formed in response to salt withdrawal, and adjacent inflation, provides an alternative mechanism for the origin of the proposed Silverpit impact crater.     

Tetrahedral Plots of the Phase Relations for Basalts

The phase relations of quaternary systems are generally represented by projections onto ternary compositional planes. Such projections often obscure relationships that would only be evident in a three-dimensional tetrahedral plot. The tetrahedral plot requires that compositions of the minerals and melts be transformed into Cartesian coordinates. It is shown here how this transformation is carried out. The application is demonstrated by tetrahedral plots of experimental melt compositions of partially molten lherzolite. Furthermore, the plot can be used to evaluate whether or not a particular basaltic composition represents a primary melt. The methods are applicable to any four-component system.     

Fossil wood

Plants have been making wood (secondary xylem) for more than 370 million years. This familiar material is one of the keys to their massive success. Wood allows plants to attain breathtaking heights and maximize the capture of sunlight for growth. By creating complex, multi‐layered forests, the evolution of wood, has done more to shape life on land than almost any other innovation. Wood is one of the most common terrestrial fossils encountered in the geological record. Whether preserved as huge petrified logs or as minute chunks of charcoal, fossil wood is abundant in rocks of Late Devonian age and younger. It is of enormous scientific importance, shedding light on the identity and stature of ancient trees. As a record of growth over a sustained period, it also tells us much about the climates and environments in which those trees lived. In this article, I explain some aspects of the origin, evolution, preservation, and identification of fossil wood, and emphasise its great significance for geology.     

Constitutive and numerical framework for modeling joints and faults in rock

A three‐dimensional constitutive model for joints is described that incorporates nonlinear elasticity based on volumetric elastic strain, and plasticity for both compaction and shear with emphasis on compaction. The formulation is general in the sense that alternative specific functional forms and evolution equations can be easily incorporated. A corresponding numerical structure based on finite elements is provided so that a joint width can vary from a fraction of an element size to a width that occupies several elements. The latter case is particularly appropriate for modeling a fault, which is considered simply to be a joint with large width. For small joint widths, the requisite equilibrium and kinematic requirements within an element are satisfied numerically. The result is that if the constitutive equation for either the joint or the rock is changed, the numerical framework remains unchanged. A unique aspect of the general formulation is the capability to handle either pre‐existing gaps or the formation of gaps. Representative stress–strain plots are given to illustrate both the features of the model and the effects of changes in values of material parameters. Copyright © 2015 John Wiley & Sons, Ltd.     

Glacial conditions that contribute to the regeneration of Fountain Glacier proglacial icing,Bylot Island,Canada

Proglacial icings are one of the most common forms of extrusive ice found in the Canadian Arctic. However, the icing adjacent to Fountain Glacier, Bylot Island, is unique due to its annual cycle of growth and decay, and perennial existence without involving freezing point depression of water due to chemical characteristics. Its regeneration depends on the availability of subglacial water and on the balance between ice accretion and hydro‐thermal erosion. The storage and conduction of the glacial meltwater involved in the accretion of the icing were analyzed by conducting topographic and ground penetrating radar surveys in addition to the modelling of the subglacial drainage network and the thermal characteristics of the glacier base. The reflection power analysis of the geophysical data shows that some areas of the lower ablation zone have a high accumulation of liquid water, particularly beneath the centre part of the glacier along the main supraglacial stream. A dielectric permittivity model of the glacier – sediment interface suggests that a considerable portion of the glacier is warm based; allowing water to flow through unfrozen subglacial sediments towards the proglacial outwash plain. All these glacier‐related characteristics contribute to the annual regeneration of the proglacial icing and allow for portions of the icing to be perennial. Copyright © 2013 John Wiley & Sons, Ltd.     

Spatially controlled Fe and Si isotope variations: an alternative view on the formation of the Torres del Paine pluton

We present new Fe and Si isotope ratio data for the Torres del Paine igneous complex in southern Chile. The multi-composition pluton consists of an approximately 1 km vertical exposure of homogenous granite overlying a contemporaneous 250-m-thick mafic gabbro suite. This first-of-its-kind spatially dependent Fe and Si isotope investigation of a convergent margin-related pluton aims to understand the nature of granite and silicic igneous rock formation. Results collected by MC-ICP-MS show a trend of increasing δ⁵⁶Fe and δ³⁰Si with increasing silica content as well as a systematic increase in δ⁵⁶Fe away from the mafic base of the pluton. The marginal Torres del Paine granites have heavier Fe isotope signatures (δ⁵⁶Fe = +0.25 ± 0.02 2se) compared to granites found in the interior pluton (δ⁵⁶Fe = +0.17 ± 0.02 2se). Cerro Toro country rock values are isotopically light in both Fe and Si isotopic systems (δ⁵⁶Fe = +0.05 ± 0.02 ‰; δ³⁰Si = ?0.38 ± 0.07 ‰). The variations in the Fe and Si isotopic data cannot be accounted for by local assimilation of the wall rocks, in situ fractional crystallization, late-stage fluid exsolution or some combination of these processes. Instead, we conclude that thermal diffusion or source magma variation is the most likely process producing Fe isotope ratio variations in the Torres del Paine pluton.