The influence of thermal and cyclic stressing on the strength of rocks from Mount St. Helens, Washington

Stratovolcanoes and lava domes are particularly susceptible to sector collapse resulting from wholesale rock failure as a consequence of decreasing rock strength. Here, we provide insights into the influence of thermal and cyclic stressing on the strength and mechanical properties of volcanic rocks. Specifically, this laboratory study examines the properties of samples from Mount St. Helens; chosen because its strength and stability have played a key role in its history, influencing the character of the infamous 1980 eruption. We find that thermal stressing exerts different effects on the strengths of different volcanic units; increasing the heterogeneity of rocks in situ. Increasing the uniaxial compressive stress generates cracking, the timing and magnitude of which was monitored via acoustic emission (AE) output during our experiments. AEs accelerated in the approach to failure, sometimes following the pattern predicted by the failure forecast method (Kilburn 2003). Crack damage during the experiments was tracked using the evolving static Young’s modulus and Poisson’s ratio, which represent the quasi-static deformation in volcanic edifices more accurately than dynamic elastic moduli which are usually implemented in volcanic models. Cyclic loading of these rocks resulted in a lower failure strength, confirming that volcanic rocks may be weakened by repeated inflation and deflation of the volcanic edifice. Additionally, volcanic rocks in this study undergo significant elastic hysteresis; in some instances, a material may fail at a stress lower than the peak stress which has previously been endured. Thus, a volcanic dome repeatedly inflated and deflated may progressively weaken, possibly inducing failure without necessarily exceeding earlier conditions.     

Controls on the geotechnical response of sedimentary rocks to weathering

Weathering reduces the strength of rocks and so is a key control on the stability of rock slopes. Recent research suggests that the geotechnical response of rocks to weathering varies with ambient stress conditions resulting from overburden loading and/or stress concentrations driven by near-surface topography. In addition, the stress history experienced by the rock can influence the degree to which current weathering processes cause rock breakdown. To address the combined effect of these potential controls, we conducted a set of weathering experiments on two sedimentary lithologies in laboratory and field conditions. We firstly defined the baseline geotechnical behaviour of each lithology, characterising surface hardness and stress–strain behaviour in unconfined compression. Weathering significantly reduced intact rock strength, but this was not evident in measurements of surface hardness. The ambient compressive stress applied to samples throughout the experiments did not cause any observable differences in the geotechnical behaviour of the samples. We created a stress history effect in sub-sets of samples by generating a population of microcracks that could be exploited by weathering processes. We also geometrically modified groups of samples to cause near-surface stress concentrations that may allow greater weathering efficacy. However, even these pronounced sample modifications resulted in insignificant changes in geotechnical behaviour when compared to unmodified samples. The observed reduction in rock strength changed the nature of failure of the samples, which developed post-peak strength and underwent multiple stages of brittle failure. Although weakened, these samples could sustain greater stress and strain following exceedance of peak strength. On this basis, the multi-stage failure style exhibited by weaker weathered rock may permit smaller-magnitude, higher-frequency events to trigger fracture through intact rock bridges as well as influencing the characteristics of pre-failure deformation. These findings are consistent with patterns of behaviour observed in field monitoring results. © 2019 John Wiley & Sons, Ltd.     

岩石破裂强度的温度和应变率效应及其对岩石圈流变结构的影响

研究了温度和应变率对岩石破裂强度的影响，得到了岩石圈中一些典型岩石破裂强度的新的经验规律.新的经验规律除考虑围压和标本尺度的影响外，还考虑了温度和应变率的影响，并增加了新岩石的结果，所以更能反映岩石圈内岩石破裂的真实状态.通过对鄂尔多斯平均流变结构的计算和对比研究表明：传统的忽略脆性破裂的流变模型过高地估计了流变强度，流变机制的分布也不尽合理.而考虑了脆性破裂机制的流变模型的结果表明脆性区分为两部分，浅部以摩擦滑动机制控制，深部以脆性破裂机制控制.由于新的经验规律考虑的代表性岩石更全面，并考虑了应变率的影响，得到的脆性区的范围进一步增大,流变强度进一步降低.     

地壳岩石变形行为的转变及其温压条件

岩石脆延性转化 (brittle ductiletransition)和脆塑性转化 (brittle plastictransition)是不同的概念。脆延性转化指从岩石的局部变形破坏到宏观均匀流动变形的转化 ,它与宏观结构和力学行为的变化相关。脆塑性转化指脆性向晶体塑性变形的转化 ,它与力学行为和微观机制的变化相关。通过地壳中最主要的石英、长石的实验室和野外变形温压条件对比发现 ,达到相同的变形特征 ,在实验室和野外所需温压条件不同。建立变形机制图使解决这一矛盾成为可能。但受实验资料的限制 ,目前几种主要岩石的变形机制图还无法建立。因此 ,通过对实验与自然环境下变形特征及微观机制对比 ,找出两者温压条件的差别 ,就成为将实验研究结果外推解决实际地质问题的有效途径     

On the state of stress in the near-surface of the earth's crust

Five models for near-surface crustal stresses induced by gravity and horizontal deformation and the influence of rock property contrasts, rock strength, and stress relaxation on these stresses are presented. Three of the models—the lateral constraint model, the model for crustal stresses caused by horizontal deformation, and the model for the effects of anisotropy—are linearly elastic. The other two models assume that crustal rocks are brittle or viscoelastic in order to account for the effects of rock strength and time on near-surface stresses. It is shown that the lateral constraint model is simply a special case of the combined gravity-and deformation-induced stress field when horizontal strains vanish and that the inclusion of the effect of rock anisotropy in the solution for crustal stresses caused by gravity and horizontal deformation broadens the range for predicted stresses. It is also shown that when stress levels in the crust reach the limits of brittle rock strength, these stresses become independent of strain rates and that stress relaxation in ductile crustal rocks subject to constant horizontal strain rates causes horizontal stresses to become independent of time in the long term.     

True Triaxial Stresses and the Brittle Fracture of Rock

This paper reviews the efforts made in the last 100 years to characterize the effect of the intermediate principal stress σ ₂ on brittle fracture of rocks, and on their strength criteria. The most common theories of failure in geomechanics, such as those of Coulomb, and Mohr, disregard σ ₂ and are typically based on triaxial testing of cylindrical rock samples subjected to equal minimum and intermediate principal stresses (σ ₃=σ ₂). However, as early as 1915 Böker conducted conventional triaxial extension tests (σ ₁=σ ₂) on the same Carrara marble tested earlier in conventional triaxial compression by von Kármán that showed a different strength behavior. Efforts to incorporate the effect of σ ₂ on rock strength continued in the second half of the last century through the work of Nadai, Drucker and Prager, Murrell, Handin, Wiebols and Cook, and others. In 1971 Mogi designed a high-capacity true triaxial testing machine, and was the first to obtain complete true triaxial strength criteria for several rocks based on experimental data. Following his pioneering work, several other laboratories developed equipment and conducted true triaxial tests revealing the extent of σ ₂ effect on rock strength (e.g., Takahashi and Koide, Michelis, Smart, Wawersik). Testing equipment emulating Mogi's but considerably more compact was developed at the University of Wisconsin and used for true triaxial testing of some very strong crystalline rocks. Test results revealed three distinct compressive failure mechanisms, depending on loading mode and rock type: shear faulting resulting from extensile microcrack localization, multiple splitting along the σ ₁ axis, and nondilatant shear failure. The true triaxial strength criterion for the KTB amphibolite derived from such tests was used in conjunction with logged breakout dimensions to estimate the maximum horizontal in situ stress in the KTB ultra deep scientific hole.     

Volcano collapse promoted by progressive strength reduction: new data from Mount St. Helens

Rock shear strength plays a fundamental role in volcano flank collapse, yet pertinent data from modern collapse surfaces are rare. Using samples collected from the inferred failure surface of the massive 1980 collapse of Mount St. Helens (MSH), we determined rock shear strength via laboratory tests designed to mimic conditions in the pre-collapse edifice. We observed that the 1980 failure shear surfaces formed primarily in pervasively shattered older dome rocks; failure was not localized in sloping volcanic strata or in weak, hydrothermally altered rocks. Our test results show that rock shear strength under large confining stresses is reduced ∼20% as a result of large quasi-static shear strain, as preceded the 1980 collapse of MSH. Using quasi-3D slope-stability modeling, we demonstrate that this mechanical weakening could have provoked edifice collapse, even in the absence of transiently elevated pore-fluid pressures or earthquake ground shaking. Progressive strength reduction could promote collapses at other volcanic edifices.     

应力途径对岩石脆性-延性变化的影响

用两种加载方式对岩石的破裂进行了实验研究。一种方式是在一定的围压下增加轴压使岩石破裂（A型）;另一种方式是在一定围压下增加轴压直到破裂前某一应力状态,然后停止加轴压转而减小围压使岩石破裂（B型）。所用的岩石样品为济南辉长岩和山东掖县白大理岩。着重研究应力途径对岩石脆性-延性变化的影响。辉长岩在1.5千巴以下两种应力途径下的破裂都表现为脆性,但是对于同样的应力状态,B型实验比A型实验显得更脆。随着围压增加到200-250巴之间,大理岩由脆性转变为延性。围压250巴以上,大理岩的A型实验发生延性破裂。样品承受载荷的能力是逐渐丧失的。破裂过程中声发射率极低,听不到破裂声响。最后在样品中形成了剪切断面,但破裂很慢。然而,在大理岩的B型实验中,围压在250巴以上发生了脆性破裂,其表现为轴向应力突然下降,伴随着脆性破裂的声响并有声发射率剧增的前兆。看来,B型应力途径对岩石起了一个脆化的作用。     

Permeability Evolution and Rock Brittle Failure

This paper reports an experimental study of the evolution of permeability during rock brittle failure and a theoretical analysis of rock critical stress level. It is assumed that the rock is a strain-softening medium whose strength can be described by Weibull’s distribution. Based on the two-dimensional renormalization group theory, it is found that the stress level λ_c (the ratio of the stress at the critical point to the peak stress) depends mainly on the homogeneity index or shape parameter m in the Weibull’s distribution for the rock. Experimental results show that the evolution of permeability is closely related to rock deformation stages: the permeability has a rapid increase with the growth of cracks and their surface areas (i.e., onset of fracture coalescence point), and reaches the maximum at rock failure. Both the experimental and analytical results show that this point of rapid increase in permeability on the permeability-pressure curve corresponds to the critical point on the stress-strain curve; for rock compression, the stress at this point is approximately 80% of the peak strength. Thus, monitoring the evolution of permeability may provide a new means of identifying the critical point of rock brittle fracture.     

Pore fluid pressure effects in anisotropic rocks: Mechanisms of induced seismicity and weak faults

Many observations and studies indicate that pore fluid pressure in the crustal rocks plays an important role in deformation, faulting, and earthquake processes. Conventional models of pore pressure effects often assume isotropic porous rocks and yield the nondeviatoric pressure effects which seem insufficient to explain diverse phenomena related to pore pressure variation, such as fluid-extraction induced seismicity and crustal weak faults. We derive the anisotropic effective stress law especially for transversely-isotropic and orthotropic rocks, and propose that the deviatoric effects of pore fluid pressure in anisotropic rocks not only affect rock effective strength but also cause variation of shear stresses. Such shear stress variations induced by either pore pressure buildup or pore pressure decline may lead to faulting instability and trigger earthquakes, and provide mechanisms for the failure of crustal weak faults with low level of shear stresses. We believe that the deviatoric effects of pore fluid pressure in anisotropic rocks are of wide application in studies of earthquake precursors and aftershocks, oil and gas reservoir characterization, enhanced oil recovery, and hydraulic fracturing.     

Estimates of possible diffusional effects in trace element separations

Igneous material dredged from the Rio Grande rise, South Atlantic Ocean, includes basaltic rocks, some having mafic nodules and megacrysts, and volcanic breccias composed largely of basaltic fragments. These samples represent the only volcanic rocks recovered from this aseismic rise. Bulk compositions show alkalic basalt, trachybasalt, and trachyandesite; the rock types are similar to those of nearby Tristan da Cunha, Gough, and the Walvis ridge. Microprobe analyses show basaltic groundmass to have olivine, Fo₈₅, pyroxene, Fs₁₃Wo₄₆, feldspar, An₇₁, plus interstitial alkali feldspar. Mafic nodules and megacrysts have olivine, Fo_86–90 and pyroxene Fs_6–7.5Wo_45–46; Al₂O₃ 2.5–4 wt.%.The Rio Grande rise rocks have compositional characteristics of an alkalic basaltic suite, and not of mid-ocean ridge tholeiite. Based on mineral compositions, nodules and megacrysts in basalt are interpreted as cognate inclusions. Because oceanic alkalic basaltic rocks are almost invariably associated with islands and seamounts, the Rio Grande rise probably represents a series of alkalic-basalt islands that formed and eventually subsided during rifting of the South Atlantic; the dredged volcanic breccias are probably slump deposits from those volcanoes. This interpretation lends support to the Rio Grande rise having formed at a hot spot, but the possibility of alkalic rocks having formed along fracture zones should not be discounted.     

Damage evolution and fluid flow in poroelastic rock

We present a formulation for mechanical modeling of the interaction between fracture and fluid flow. Our model combines the classic Biot poroelastic theory and a damage rheology model. The model provides an internally consistent framework for simulating coupled evolution of fractures and fluid flow together with gradual transition from brittle fracture to cataclastic flow in high-porosity rocks. The theoretical analysis, based on thermodynamic principles, leads to a system of coupled kinetic equations for the evolution of damage and porosity. A significant advantage of the model is the ability to reproduce the entire yield curve, including positive and negative slopes, in high-porosity rocks by a unified formulation. A transition from positive to negative values in the yield curve, referred to as a yield cap, is determined by the competition between the two thermodynamic forces associated with damage and porosity evolution. Numerical simulations of triaxial compression tests reproduce the gradual transition from localized brittle failure to distributed cataclastic flow with increasing pressure in high-porosity rocks and fit well experimentally measured yield stress for Berea sandstone samples. We modified a widely used permeability porosity relation by accounting for the effect of damage intensity on the connectivity. The new damage-permeability relation, together with the coupled kinetics of damage and porosity evolution, reproduces a wide range of realistic features of rock behavior. We constrain the model variables by comparisons of the theoretical predictions with laboratory results reporting porosity and permeability variation in rock samples during isotropic and anisotropic loading. The new damage-porosity-permeability relation enables simulation of coupled evolution of fractures and fluid flow and provides a possible explanation for permeability measurements in high-porosity rocks, referred to as the “apparent permeability paradox.” The text was submitted by the authors in English.     

The dykes and structural setting of the volcanic front in the Lesser Antilles island arc

The orientations of dykes from many of the islands of the Lesser Antilles island arc have been mapped. Most of these dykes can be interpreted in terms of local or regional swarms derived from specific volcanoes of known age, with distinct preferred orientations. Dykes are known from all Cenozoic epochs except the Palaeocene, but are most common in Pliocene, Miocene and Oligocene rocks. A majority of the sampled dykes are basaltic, intrude volcaniclastic host rocks and show a preference for widths of 1–1.25 m. Locally, dyke swarms dilate their hosts by up to 9% over hundreds of metres and up to 2% over distances of kilometres. The azimuths of dykes of all ages show a general NE-SW preferred orientation with a second NW-SE mode particularly in the Miocene rocks of Martinique. The regional setting for these minor intrusions is a volcanic front above a subduction zone composed of three segments: Saba-Montserrat, Guadeloupe-Martinique, St. Lucia-Grenada. The spacing of volcanic centres along this front is interpreted in terms of rising plumes of basaltic magma spaced about 30 km apart. This magma is normally intercepted at crustal depths by dioritic plutons and andesitic/dacitic magma generated there. Plumes which intersect transverse fracture systems or which migrate along the front can avoid these crustal traps. Throughout its history the volcanic front as a whole has migrated, episodically, towards the backarc at an average velocity of about 1 km/Ma. The local direction of plate convergence is negatively correlated with the local preferred orientation of dykes. The dominant NE-SW azimuth mode corresponds closely to the direction of faulting in the sedimentary cover of the backarc and the inferred tectonic fabric of the oceanic crust on which the arc is founded. A generalised model of the regional stress field that controls dyke intrusion outside of the immediate vicinity of central volcanic vents is proposed, in which the maximum horizontal stress parallels the volcanic front except in the northern segment where subduction of the Barracuda Rise perturbs the stress field. There is also evidence of specific temporal changes in the stress field that are probably due to large scale plate kinematics.     

Submarine basaltic volcanism: Morphometric parameters for discriminating hyaloclastites from hyalotuffs

Hyaloclastites are volcaniclastic rocks generated by non-explosive granulation of volcanic glass which takes place when basaltic magmas are quenched by contact with water; hyaloclastites are common products of deep submarine basaltic central volcanoes (seamounts and guyots). We suggested (Honnorlz, 1966) calling hyalotuffs the actual pyroclastic rocks which are generated by phreatomagmatic and phreatic explosions taking place when basaltic volcanoes crupt in shallow waters; hyalotuffs are restricted to shallow subaqueous conditions since no volcanic explosion can occur in deep seas. The distinction between hyaloclastites and hyalotuffs is therefore a useful tool when reconstructing the paleogeography of ancient submarine volcanic edifices and the mechanism by which their lavas were emplaced. We propose using two sets of morphometric parameters to discriminate the hyaloclastites from the hyalotuffs. The granules making up these two volcaniclastic rock types plot in different areas of either 1) a ternary diagram the apexes of which represent the grain planarity (P), convexity (V) and concavity (C) % («roundness» according to Szadeczky-Kardoss); or 2) a binary diagram relating the number of grain corners (N) to their planarity (P) %.     

The effect of water on stress relaxation of faulted and unfaulted sandstone

A series of stress relaxation experiments have been carried out on faulted and intact Tennessee sandstone to explore the influence of pore water on strength at different strain rates. Temperatures employed were 20, 300 and 400°C, effective confining pressure was 1.5 kb and strain rates as low as 10^–10 sec^–1 were achieved. Most samples were prefaulted at 2.5 kb confining pressure and room temperature. This is thought to have secured a reproducible initial microstructure.The strength of the dry rock was almost totally insensitive to strain rate in the range 10^–4 to 10^–10 sec^–1. In contrast, the strength of the wet rock decreased rapidly with strain rate at rates less than 10^–6 sec^–1. Brittle fracture of the quartz grains which constitute this rock is the most characteristic mode of failure under the test conditions used.The experimental data are discussed in terms of the possible deformation rate controlling processes, and it is suggested that in the wet experiments at intermediate to high strain rates (10^–7 to 10^–4 sec^–1) the observed deformation rate is controlled by the kinetics of water assisted stress corrosion, whilst deformation at low strain rates (ca. 10^–9 sec^–1) is controlled by a pressure solution process.The results have implications for the rheology of fault rocks at depths of perhaps 10 to 15 km in sialic crust.     

岩石加速破裂行为的物理自相似律

掌握岩石变形破坏过程中体积膨胀点至峰值强度点之间加速破裂行为的演化规律,是实现地质灾害物理预测的关键.本文考虑裂纹张开和闭合两种情况,基于断裂力学建立了三轴应力作用下裂纹扩展临界尺度与等效应力的关系.对微元体破坏概率,分别采用以等效应力表达的Weibull分布函数和裂纹尺度分形函数,通过对比导出了形状参数m与裂纹分布分维D_f关系的表达式.一个有趣的发现是,岩石峰值强度点与体积膨胀点应变比仅与m或D_f有关.对岩石蠕变或准蠕变破坏,合理的m值范围为[1.0,4.0],在此范围内应变比近似为常数1.48,该常数是描述不同尺度岩石加速破裂规律的物理自相似常数.实例分析表明,基于岩石加速破裂规律构建的多锁固段脆性破裂理论,其适用性广,尤其在崩滑和大地震预测领域,具有良好应用前景.此外,本文给出了b值与m值定量关系,以解释b值的物理意义,并探讨将其用于地震预测的可行性.     

基于能量耗散原理的非均质岩石损伤破坏元胞自动机模型研究

基于能量耗散理论建立非均质岩石的动态损伤破坏元胞自动机,分析单轴压缩试验中岩石破坏截面的损伤状态,得到该情况下的岩石裂隙微观动态发展过程、损伤演化关系以及全程应力应变曲线。研究发现：①加载过程中,均质度较低的岩样裂纹的萌生和扩展较为分散,损伤速率低;均质度较高的岩样裂纹的萌生和扩展非常集中,表现为脆性破坏。②损伤演化曲线呈3阶段S形发展。③随着岩石均质度参数的增加,岩石的峰值强度和峰值应变都有所提高,峰后曲线越来越陡。④给出的衡量岩石脆性破坏强弱程度的指标参数,能够较好地描述岩石破坏形式随着均值度参数m变化的规律。     

DEFORMATION OF THE BRITTLE-PLASTIC TRANSITION ZONE AT THE POST-SEISMIC RELAXATION PERIOD: A CASE STUDY OF THE RED RIVER FAULT

The transition from microscopic brittle deformation to microscopic plastic deformation is called brittle-plastic transition, which is considered as a key layer for determining the limit of lower continental crust seismicity. The depth and deformation mechanism of the brittle-plastic transition zone is controlled mainly by temperature. Besides, the strain rate and fluid pore pressure also affect the transition during the different deformation stages at the seismic cycle. In this paper, microstructure observation of catalcastic samples collected from the Red River Fault was carried out using optical polarized microscopy and scanning electron microscopy. The morphology, microstructures of deformation characteristics, mineral composition, water-rock reaction, pressure solution, exsolution, crack healing in the samples were systematically observed. The mineral components quantitative analyses were examined using the EDS. Water-rock reaction and pressure solution were systematically observed under SEM. The fabric of the main minerals in the samples was measured using electron backscattered diffraction(EBSD). Based on these analyses, the deformation mode was setup for the brittle-plastic transition zone of the fault during the post-seismic relaxation period. Both brittle deformation and plastic deformation were developed in the cataclastic samples. EBSD data shows that the c axial fabrics of quartz present low-temperature plastic deformation characteristics. The feldspar deformed as cataclastic rock, and the micro-fracture in feldspar was healed by static recrystallized quartz and calcite veins. The calcite vein underwent plastic deformation, which represents the post-seismic relaxation deformation. Based on the analysis of deformation mechanism of cataclastic samples in brittle-plastic transition zone of the Red River Fault, and combined with previous studies, we concluded that the brittle fracture and fracture healing is the main deformation mode at brittle-plastic transition zone in the post-seismic relaxation. High stress and high strain rate at post-seismic relaxation lead to brittle fracture of high-strength minerals such as feldspar in rocks. Plastic deformation occurs in low-strength minerals such as quartz and mica. Under the fluid condition, micro-fractures were healed by quartz and calcite. The minerals such as quartz and calcite in the fracture transformed from static recrystallization to dynamic recrystallization with stress gradually accumulating. With fracture healing and stress accumulation, the fault strength gradually increases which could accumulate energy for the next earthquake.     

A new model for the fracture energy budget of phreatomagmatic explosions

The fragmentation of magma and of the hosting country rocks is a major process in explosive eruptions. It is important to quantify the mechanical energy needed for fragmentation in order to assess the physical processes of this volcanic phenomenon. This paper presents a method to calculate the fragmentation energy of country rock using granulometry data of a typical phreatomagmatic Eifel maar volcano explosion. The total fracture area of country rock fragments in one tephra layer was quantified and related to the critical fragmentation energy of these country rocks. The rock parameters critical shear stress and critical fragmentation energy were determined experimentally, whereas the pre-volcanic crack inventory was measured in the field. The paper concludes with the calculation of the energy balance (i.e. partitioning of thermal energy into kinetical energy and mechanical energy of the fragmentation) of one Eifel maar volcanic explosion.     

图们江流域新生代火山岩Sr、Nd同位素初步研究

中国地震局地质研究所; 中国地震局地质研究所 北京