SANICLAY: simple anisotropic clay plasticity model

SANICLAY is a new simple anisotropic clay plasticity model that builds on a modification of an earlier model with an associated flow rule, in order to include simulations of softening response under undrained compression following K_o consolidation. Non‐associativity is introduced by adopting a yield surface different than the plastic potential surface. Besides, the isotropic hardening of the yield surface both surfaces evolve according to a combined distortional and rotational hardening rule, simulating the evolving anisotropy. Although built on the general premises of critical state soil mechanics, the model induces a critical state line in the void ratio–mean effective stress space, which is a function of anisotropy. To ease interpretation, the model formulation is presented firstly in the triaxial stress space and subsequently, its multiaxial generalization is developed systematically, in a form appropriate for implementation in numerical codes. The SANICLAY is shown to provide successful simulation of both undrained and drained rate‐independent behaviour of normally consolidated sensitive clays, and to a satisfactory degree of accuracy of overconsolidated clays. The new model requires merely three constants more than those of the modified Cam clay model, all of which are easily calibrated from well‐established laboratory tests following a meticulously presented procedure. Copyright © 2006 John Wiley & Sons, Ltd.     

Experimental constraints on the effects of pressure and H2O on the fractional crystallization of high-Mg island arc basalt

H₂O-undersaturated melting experiments of synthesized basalt (SiO₂ = 50.7 wt.%, MgO = 8.3 wt.%, Mg# = 60) were conducted at fO₂ corresponding to NNO+1 and NNO−1 to clarify the effects of pressure (2–7 kbar) and H₂O on fractional crystallization in island arcs. H₂O content was ranged from nominally anhydrous to 4.4 wt.%. Differentiation trends, namely the liquid lines of descent, change sensitively according to pressure-H₂O relations. Tholeiitic differentiation trends are reproduced with H₂O ≤ ∼2 wt.% in primary magma. With such quantities of H₂O, fractional crystallization is controlled by olivine + plagioclase at 2 kbar. Increasing the pressure from 2 to ≥4 kbar induces early crystallization of orthopyroxene instead of olivine and therefore SiO₂ enrichment in the residual melts is suppressed. Increasing H₂O (≥ ∼2 wt.% in primary magma) stabilizes clinopyroxene relative to orthopyroxene and/or magnetite. Although the phase relations and proportions strongly depend on fO₂ and H₂O content, differentiation trends are always calc-alkaline.     

Pressures of Crystallization of Icelandic Magmas

Iceland lies astride the Mid-Atlantic Ridge and was createdby seafloor spreading that began about 55 Ma. The crust is anomalouslythick (20–40 km), indicating higher melt productivityin the underlying mantle compared with normal ridge segmentsas a result of the presence of a mantle plume or upwelling centeredbeneath the northwestern edge of the Vatnajökull ice sheet.Seismic and volcanic activity is concentrated in 50 km wideneovolcanic or rift zones, which mark the subaerial Mid-AtlanticRidge, and in three flank zones. Geodetic and geophysical studiesprovide evidence for magma chambers located over a range ofdepths (1·5–21 km) in the crust, with shallow magmachambers beneath some volcanic centers (Katla, Grimsvötn,Eyjafjallajökull), and both shallow and deep chambers beneathothers (e.g. Krafla and Askja). We have compiled analyses ofbasalt glass with geochemical characteristics indicating crystallizationof ol–plag–cpx from 28 volcanic centers in the Western,Northern and Eastern rift zones as well as from the SouthernFlank Zone. Pressures of crystallization were calculated forthese glasses, and confirm that Icelandic magmas crystallizeover a wide range of pressures (0·001 to 1 GPa), equivalentto depths of 0–35 km. This range partly reflects crystallizationof melts en route to the surface, probably in dikes and conduits,after they leave intracrustal chambers. We find no evidencefor a shallow chamber beneath Katla, which probably indicatesthat the shallow chamber identified in other studies containssilica-rich magma rather than basalt. There is reasonably goodcorrelation between the depths of deep chambers (> 17 km)and geophysical estimates of Moho depth, indicating that magmaponds at the crust–mantle boundary. Shallow chambers (<7·1 km) are located in the upper crust, and probablyform at a level of neutral buoyancy. There are also discretechambers at intermediate depths (11 km beneath the rift zones),and there is strong evidence for cooling and crystallizing magmabodies or pockets throughout the middle and lower crust thatmight resemble a crystal mush. The results suggest that themiddle and lower crust is relatively hot and porous. It is suggestedthat crustal accretion occurs over a range of depths similarto those in recent models for accretionary processes at mid-oceanridges. The presence of multiple stacked chambers and hot, porouscrust suggests that magma evolution is complex and involvespolybaric crystallization, magma mixing, and assimilation. KEY WORDS: Iceland rift zones; cotectic crystallization; pressure; depth; magma chamber; volcanic glass     

Geochemistry and petrogenesis of the post-orogenic bimodal dyke swarms in NW Sinai, Egypt: constraints on the magmatic–tectonic processes during the late Precambrian

The study area in the northwest Sinai represents one of the most significant regions in the Egyptian basement intensely invaded by post-orogenic calc-alkaline dyke swarms. Two post-orogenic dyke swarms have been recognized in NW Sinai namely: (1) mafic dykes of basalt, basaltic andesite and andesite composition and (2) felsic dykes of dacite, rhyodacite and rhyolite composition. These basaltic to rhyolitic dykes intruded contemporaneously and shortly after the intrusion of the post-orogenic leucogranite. The mafic and felsic dykes are enriched in incompatible elements, especially in the large ion lithophile elements (e.g. K, Rb, Ba) and depleted in high field strength elements with negative P, Ti and Nb anomalies. Major and trace element geochemistry indicates that investigated mafic and felsic magma types are not related via fractional crystallization. The protoliths of the mafic and felsic dykes appear to have evolved from different parental magmas. The incompatible trace element patterns favour a derivation of the mafic dykes from melting of a garnet peridotite source followed by fractional crystallization of olivine, clinopyroxene amphibole and zircon. The felsic dykes, on the other hand, could be generated by melting of garnet-free source modified subsequently by fractional crystallization of plagioclase, apatite and titanomagnetite. This implies variable source characteristics at the end of the Pan-African in the NW Sinai.The mafic and felsic dykes can be related to an intracontinental setting and that this was accompanied by a chemical evolution of the subcontinental lithosphere. Magma generation and ascent in the area was favoured by extensional movements, which is already known from other areas in NE Africa.     

Computation of Large Anisotropic Seismic Heterogeneities (CLASH)

A general tomographic technique is designed in order (i) to operate in anisotropic media; (ii) to account for the uneven seismic sampling and (iii) to handle massive data sets in a reasonable computing time. One modus operandi to compute a 3-D body wave velocity model relies on surface wave phase velocity measurements. An intermediate step, shared by other approaches, consists in translating, for each period of a given mode branch, the phase velocities integrated along ray paths into local velocity perturbations. To this end, we develop a method, which accounts for the azimuthal anisotropy in its comprehensive form. The weakly non-linear forward problem allows to use a conjugate gradient optimization. The Earth's surface is regularly discretized and the partial derivatives are assigned to the individual grid points. Possible lack of lateral resolution, due to the inescapable uneven ray path coverage, is taken into account through the a priori covariances on parameters with laterally variable correlation lengths. This method allows to efficiently separate the 2ψ and the 4ψ anisotropic effects from the isotropic perturbations. Fundamental mode and overtone phase velocity maps, derived with real Rayleigh wave data sets, are presented and compared with previous maps. The isotropic models concur well with the results of Trampert & Woodhouse. Large 4ψ heterogeneities are located in the tectonically active regions and over the continental lithospheres such as North America, Antarctica or Australia. At various periods, a significant 4ψ signature is correlated with the Hawaii hotspot track. Finally, concurring with the conclusions of Trampert & Woodhouse, our phase velocity maps show that Rayleigh wave data sets do need both 2ψ and 4ψ anisotropic terms.     

裂隙型单斜介质中弹性系数的计算及波的传播特性研究

根据Hudson等关于裂隙介质弹性系数计算的扰动理论及Bond变换矩阵原理，给出了各向同性介质中含多组垂直裂隙时等效弹性系数的计算方法。计算了含两组斜交的垂直裂隙形成的单斜各向异性介质中的等效弹性系数，并根据Christoffel方程推导、得出水平面内平面波传播的相速度和群速度随方位变化的特性。