Magma Flow Instabilities in a Volcanic Conduit: Implications for Long-Period Seismicity

Silicic volcanic eruptions are typically accompanied by repetitive Long-Period (LP) seismicity that originates from a small region of the upper conduit. These signals have the capability to advance eruption prediction, since they commonly precede a change in the eruption vigour. Shear bands forming along the conduit wall, where the shear stresses are highest, have been linked to providing the seismic trigger. However, existing computational models are unable to generate shear bands at the depths where the LP signals originate using simple magma strength models. Presented here is a model in which the magma strength is determined from a constitutive relationship dependent upon crystallinity and pressure. This results in a depth-dependent magma strength, analogous to planetary lithospheres. Hence, in shallow highly-crystalline regions a macroscopically discontinuous brittle type of deformation will prevail, whilst in deeper crystal-poor regions there will be a macroscopically continuous plastic deformation mechanism. This will result in a depth where the brittle-ductile transition occurs, and here shear bands disconnected from the free-surface may develop. We utilize the Finite Element Method and use axi-symmetric coordinates to model magma flow as a viscoplastic material, simulating quasi-static shear bands along the walls of a volcanic conduit. Model results constrained to the Soufrière Hills Volcano, Montserrat, show the generation of two types of shear bands: upper-conduit shear bands that form between the free-surface to a few 100 metres below it and discrete shear bands that form at the depths where LP seismicity is measured to occur corresponding to the brittle-ductile transition and the plastic shear region. It is beyond the limitation of the model to simulate a seismic event, although the modelled viscosity within the discrete shear bands suggests a failure and healing cycle time that supports the observed LP seismicity repeat times. However, due to the paucity of data and large parameter space available these results can only be considered to be qualitative rather than quantitative at this stage.

Coaxial flattening at deep levels of orogenic belts: evidence from blueschists and eclogites on Syros and Sifnos (Cyclades, Greece)

This work presents new Structural data from a high-pressure/low-temperature (HP/LT) metamorphic terrane exposed on the islands of Syros and Sifnos (Cyclades, Greece). The structure and the metamorphism of a relatively coherent HP/LT rock section were studied in order to elucidate how strain was accommodated at deep crustal levels during the formation and exhumation of HP/LT rocks. At least three deformation phases associated with eclogite- and blueschist-facies conditions (P = 8-15 kbar; T = 400-550 degreesC) were recognised. The earliest deformation fabric (S1), preserved as inclusion trails within garnet porphyroblasts, is aligned to define a sub-vertical schistosity (at present orientation), which is frequently orthogonal to the flat matrix schistosity (S2), and may indicate that deep crustal thickening involved upright folding. The currently dominant fabric in the HP rock section, S2, is Usually moderately dipping and locally contains NW-trending glaucophane lineations, symmetric pressure-shadows and eclogitic boudins. The symmetric structures associated with this fabric seem to indicate coaxial vertical thinning, although the existence of non-coaxial structures out of the study area cannot be excluded. Glaucophane-bearing shear bands (S3), with top-to-NW sense of shearing, locally crosscut the earlier structures. The latest recognised fabric (D4) is scarce and often absent within the HP rocks. It is associated with top-to-NE kinematic criteria that formed at greenschist-facies conditions (P = 4-7 kbar; T = 400-450 degreesC). Based on these observations, it is suggested that partitioning of strain occurred at different crustal levels and at different times. Deep crustal deformation was governed by thickening via upright folding followed by coaxial vertical thinning, whereas non-coaxial shearing occurred when the rocks were already exhumed to relatively shallow crustal levels. The earliest fabrics (D1 to D3) pertain to Alpine orogenesis and possibly to syn-orogenic extension, whereas the latest correspond to whole-crust back-are extension. (C) 2002 Elsevier Science Ltd. All rights reserved.

Mineral chemistry, whole-rock compositions, and petrogenesis of leg 176 gabbros: Data and discussion

We report mineral chemistry, whole-rock major element compositions, and trace element analyses on Hole 735B samples drilled and selected during Leg 176. We discuss these data, together with Leg 176 shipboard data and Leg 118 sample data from the literature, in terms of primary igneous petrogenesis. Despite mineral compositional variation in a given sample, major constituent minerals in Hole 735B gabbroic rocks display good chemical equilibrium as shown by significant correlations among Mg# (= Mg/[Mg+Fe2+]) of olivine, clinopyroxene, and orthopyroxene and An (=Ca/[Ca+Na]) of plagioclase. This indicates that the mineral assemblages olivine + plagioclase in troctolite, plagioclase + clinopyroxene in gabbro, plagioclases + clinopyroxene + olivine in olivine gabbro, and plagioclase + clinopyroxene + olivine + orthopyroxene in gabbronorite, and so on, have all coprecipitated from their respective parental melts. Fe-Ti oxides (ilmenite and titanomagnetite), which are ubiquitous in most of these rocks, are not in chemical equilibrium with olivine, clinopyroxene, and plagioclase, but precipitated later at lower temperatures. Disseminated oxides in some samples may have precipitated from trapped Fe-Ti–rich melts. Oxides that concentrate along shear bands/zones may mark zones of melt coalescence/transport expelled from the cumulate sequence as a result of compaction or filter pressing. Bulk Hole 735B is of cumulate composition. The most primitive olivine, with Fo = 0.842, in Hole 735B suggests that the most primitive melt parental to Hole 735B lithologies must have Mg# ≤ 0.637, which is significantly less than Mg# = 0.714 of bulk Hole 735B.

In situ TiB-reinforced Cu-based bulk metallic glass composites

Cu-based bulk metallic glass (BMG) composites containing in situ TiB particles were successfully fabricated. The reinforcing TiB particles with a size of 5-10 mu m are uniformly distributed in the amorphous matrix. The particles have a good bonding to the matrix with a reaction layer. The BMG composites exhibit an obvious ductility with a plastic strain of 2% for the 17.5 vol.% TiB sample due to the suppression of shear band propagation and the generation of multiple shear bands during compressive testing. The hardness of the materials is increased from Hv543 for monolithic BMG to Hv650 for 23.6 vol.% TiB-containing BMG composite. (c) 2006 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Bed shear stress measurements in dam break and swash flows

A novel shear plate was used to make direct bed shear stress measurements in laboratory dam break and swash flows on smooth, fixed, impermeable beds. The pressure gradient due to the slope of the fluid free-surface across the plate was measured using pressure transducers. Surface elevation was measured at five locations using acoustic displacement sensors. Flow velocity was measured using an Acoustic-Doppler Velocimeter and calculated using the ANUGA inundation model. The measured bed shear stress at the dam break fluid tip for an initially dry, horizontal bed was close to twice that estimated using steady flow theory. The temporal variation of swash bed shear stress showed a large peak in landward directed stress at the uprush tip, followed by a rapid decay throughout the uprush flow interior. The peak seaward directed stress during the backwash phase was less than half that measured in the uprush. Close to the still water line, in the region of bore collapse and at the time of initial uprush, favourable pressure gradients were measured. In the lower swash region predominately weak adverse pressure gradients were measured.

Determination of coalescence kernels for high-shear granulation using DEM simulations

Discrete element method (DEM) modeling is used in parallel with a model for coalescence of deformable surface wet granules. This produces a method capable of predicting both collision rates and coalescence efficiencies for use in derivation of an overall coalescence kernel. These coalescence kernels can then be used in computationally efficient meso-scale models such as population balance equation (PBE) models. A soft-sphere DEM model using periodic boundary conditions and a unique boxing scheme was utilized to simulate particle flow inside a high-shear mixer. Analysis of the simulation results provided collision frequency, aggregation frequency, kinetic energy, coalescence efficiency and compaction rates for the granulation process. This information can be used to bridge the gap in multi-scale modeling of granulation processes between the micro-scale DEM/coalescence modeling approach and a meso-scale PBE modeling approach.

Isotropic singularities in shear-free perfect fluid cosmologies

We investigate barotropic perfect fluid cosmologies which admit an isotropic singularity. From the General Vorticity Result of Scott, it is known that these cosmologies must be irrotational. In this paper we prove, using two different methods, that if we make the additional assumption that the perfect fluid is shear-free, then the fluid flow must be geodesic. This then implies that the only shear-free, barotropic, perfect fluid cosmologies which admit an isotropic singularity are the FRW models.

A director theory for visco-elastic folding instabilities in multilayered rock

A model for finely layered visco-elastic rock proposed by us in previous papers is revisited and generalized to include couple stresses. We begin with an outline of the governing equations for the standard continuum case and apply a computational simulation scheme suitable for problems involving very large deformations. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered beam under compression. We analyse folding up to 40% shortening. The standard continuum solution becomes unstable for extreme values of the shear/normal viscosity ratio. The instability is a consequence of the neglect of the bending stiffness/viscosity in the standard continuum model. We suggest considering these effects within the framework of a couple stress theory. Couple stress theories involve second order spatial derivatives of the velocities/displacements in the virtual work principle. To avoid C-1 continuity in the finite element formulation we introduce the spin of the cross sections of the individual layers as an independent variable and enforce equality to the spin of the unit normal vector to the layers (-the director of the layer system-) by means of a penalty method. We illustrate the convergence of the penalty method by means of numerical solutions of simple shears of an infinite layer for increasing values of the penalty parameter. For the shear problem we present solutions assuming that the internal layering is oriented orthogonal to the surfaces of the shear layer initially. For high values of the ratio of the normal-to the shear viscosity the deformation concentrates in thin bands around to the layer surfaces. The effect of couple stresses on the evolution of folds in layered structures is also investigated. (C) 2002 Elsevier Science Ltd. All rights reserved.

Effect of shear stress on expression of a recombinant protein by Chinese hamster ovary cells

A flow chamber was used to impart a steady laminar shear stress on a recombinant Chinese hamster ovary (CHO) cell line expressing human growth hormone (hGH). The cells were subjected to shear stress ranging from 0.005 to 0.80 N m(-2). The effect of shear stress on the cell specific glucose uptake, cell specific hGH, and lactate productivity rates were calculated. No morphological changes to the cells were observed over the range of shear stresses examined. When the cells were subjected to 0.10 N m(-2) shear in protein-free media without Pluronic F-68, recombinant protein production ceased with no change in cell morphology, whereas control cultures were expressing hGH at 0.35 mug/10(6) cells/h. Upon addition of the shear protectants, Pluronic F-68 (0.2% [w/v]) or fetal bovine serum (1.0% [v/v] FBS), the productivity of the cells was restored. The effect of increasing shear stress on the cells in protein-free medium containing Pluronic F-68 was also investigated. Cell specific metabolic rates were calculated for cells under shear stress and for no-shear control cultures performed in parallel, with shear stress rates expressed as a percentage of those obtained for control cultures. Upon increasing shear from 0.005 to 0.80 N m(-2), the cell specific hGH productivity decreased from 100% at 0.005 N m(-2) to 49% at 0.80 N m(-2) relative to the no-shear control. A concurrent increase in the glucose uptake rate from 115% at 0.01 N m(-2) to 142% at 0.80 N m(-2), and decreased lactate productivity from 92% to 50%, revealed a change in the yield of products from glucose compared with the static control. It was shown that shear stress, at sublytic levels in medium containing Pluronic F-68, could decrease hGH specific productivity. (C) 2002 Wiley Periodicals, Inc.

Texture alignment in simple shear

We illustrate the flow behaviour of fluids with isotropic and anisotropic microstructure (internal length, layering with bending stiffness) by means of numerical simulations of silo discharge and flow alignment in simple shear. The Cosserat theory is used to provide an internal length in the constitutive model through bending stiffness to describe isotropic microstructure and this theory is coupled to a director theory to add specific orientation of grains to describe anisotropic microstructure. The numerical solution is based on an implicit form of the Material Point Method developed by Moresi et al. [1].

Apparent strain localization and shear wave dispersion in elastic fault gouge with microrotations

Shear deformation of fault gouge or other particulate materials often results in observed strain localization, or more precisely, the localization of measured deformation gradients. In conventional elastic materials the strain localization cannot take place therefore this phenomenon is attributed to special types of non-elastic constitutive behaviour. For particulate materials however the Cosserat continuum which takes care of microrotations independent of displacements is a more appropriate model. In elastic Cosserat continuum the localization in displacement gradients is possible under some combinations of the generalized Cosserat elastic moduli. The same combinations of parameters also correspond to a considerable dispersion in shear wave propagation which can be used for independent experimental verification of the proposed mechanism of apparent strain localization in fault gouge.

Shear deformation at 29% solid during solidification of magnesium alloy AZ91 and aluminium alloy A356

Partially solid commercial Al-Si and Mg-Al alloys have been deformed in shear during solidification using vane rheometry. The dendritic mush was deformed for a short period at 29% solid and allowed to cool naturally after deformation. Both alloys exhibited yield point behaviour and deformation was highly localised at the surface of maximum shear stress. The short period of deformation was found to have a distinct impact on the as-cast microstructure leading to fragmented dendrites in the deformation region of both alloys. In the case of the Mg-Al alloy, a concentrated region of interdendritic porosity was also observed in the deformation region. Concentrated porosity was not observed in the Al-Si alloy. (c) 2005 Elsevier B.V. All rights reserved.