Frictional Properties and Microstructure of Calcite-Rich Fault Gouges Sheared at Sub-Seismic Sliding Velocities

We report an experimental and microstructural study of the frictional properties of simulated fault gouges prepared from natural limestone (96 % CaCO₃) and pure calcite. Our experiments consisted of direct shear tests performed, under dry and wet conditions, at an effective normal stress of 50 MPa, at 18–150 °C and sliding velocities of 0.1–10 μm/s. Wet experiments used a pore water pressure of 10 MPa. Wet gouges typically showed a lower steady-state frictional strength (μ = 0.6) than dry gouges (μ = 0.7–0.8), particularly in the case of the pure calcite samples. All runs showed a transition from stable velocity strengthening to (potentially) unstable velocity weakening slip above 80–100 °C. All recovered samples showed patchy, mirror-like surfaces marking boundary shear planes. Optical study of sections cut normal to the shear plane and parallel to the shear direction showed both boundary and inclined shear bands, characterized by extreme grain comminution and a crystallographic preferred orientation. Cross-sections of boundary shears, cut normal to the shear direction using focused ion beam—SEM, from pure calcite gouges sheared at 18 and 150 °C, revealed dense arrays of rounded, ~0.3 μm-sized particles in the shear band core. Transmission electron microscopy showed that these particles consist of 5–20 nm sized calcite nanocrystals. All samples showed evidence for cataclasis and crystal plasticity. Comparing our results with previous models for gouge friction, we suggest that frictional behaviour was controlled by competition between crystal plastic and granular flow processes active in the shear bands, with water facilitating pressure solution, subcritical cracking and intergranular lubrication. Our data have important implications for the depth of the seismogenic zone in tectonically active limestone terrains. Contrary to recent claims, our data also demonstrate that nanocrystalline mirror-like slip surfaces in calcite(-rich) faults are not necessarily indicative of seismic slip rates.     

Deformation mechanisms in experimentally and naturally deformed amphiboles

Experimental deformation of an igneous hornblendite at 600–750°C, 10 kbars confining pressure and a strain rate of 10 ^?5 sec^?1 results in kink bands whose normals cluster about [001] and axes of external rotation subparallel to [010]. This is consistent with glide on the system T = (100), t = [001]. At temperatures of 800°C to the breakdown temperature of the hornblende, external rotation axes spread along the (100) plane although kink band boundary poles remain subparallel to [001]. This indicates glide on the (100) plane in variable directions. Analysis of bend zones in a naturally deformed actinolite indicates glide on the (100), [001] system. Mechanical (101) twinning was not observed in any of the samples, and we suggest that in these amphiboles the critical resolved shear stress for glide is lower than that for twinning.     

基于团簇的岩石动态CT图像可视化与结构变形分析

岩石的变形机制是研究岩石圈变形和局部构造活动的重要基础,岩石的微观结构是研究其变形机制的关键。微观层析成像（微观CT）技术可以无损地获取样品内部高精度三维数字化图像,为观测岩石内部结构提供了技术保证。新近发展起来的动态CT技术,使得观测岩石变形过程成为可能。本文利用一组莱塔石灰岩受差应力变形的三维动态CT图像数据,通过图像切割、图像分割等处理步骤,以团簇为基本表征单位,表征不同形状、大小孔隙的动态变化过程,分析对比全部孔隙、大团簇、小团簇和裂隙型小团簇的动态图像。结果表明:随着岩石的压缩变形,其大孔隙在不断地缩小;在变形的初始阶段,靠近施压活塞的样品上部出现了密集的新生小孔隙;当岩石所受差应力到达73 MPa时,其下半部出现了明显的剪切带,具体表现为大量的新生裂隙生成并且新生裂隙集中带与轴向成约45°交角。此前的研究观测到该样品中下部若干压缩带形成的过程;本文通过基于团簇的可视化技术,观测到了高孔隙度岩石变形过程中压缩带和剪切带共同形成的过程。同时,新生的沿剪切带分布的微小裂隙并没有显示与剪切带相同的方向,而是具有各方向均匀分布的特征。该特征可能与岩石内部复杂的局部孔隙结构有关,是值得进一步探讨的问题。      

Development of faults as zones of deformation bands and as slip surfaces in sandstone

Three forms of fault are recognized in Entrada and Navajo Sandstones in the San Rafael Desert, southeastern Utah; deformation bands, zones of deformation bands, and slip surfaces. Small faults occur asdeformation bands, about one millimeter thick, in which pores collapse and sand grains fracture, and along which there are shear displacements on the order of a few millimeters or centimeters. Two or more deformation bands adjacent to each other, which share the same average strike and dip, form azone of deformation bands. A zone becomes thicker by addition of new bands, side by side. Displacement across a zone is the sum of displacements on each individual band. The thickest zones are about 0.5 m and total displacement across a thick zone rarely exceeds 30 cm. Finally,slip surfaces, which are through-going surfaces of discontinuity in displacement, form at either edge of zones of highly concentrated deformation bands. In contrast with individual deformation bands and zones of deformation bands, slip surfaces accommodate large displacements, on the order of several meters in the San Rafael Desert.The sequence of development is from individual deformation bands, to zones, to slip surfaces, and each type of faulting apparently is controlled by somewhat different processes. The formation of zones apparently involves strain hardening, whereas the formation of slip surfaces probably involves strain softening of crushed sandstone.     

Compaction and Failure in High Porosity Carbonates: Mechanical Data and Microstructural Observations

We investigated systematically the micromechanics of compaction in two carbonates of porosity above 30%, Majella grainstone and Saint Maximin limestone. The composition, grain size and pore surface area of these rocks were determined. Hydrostatic compression experiments were performed under dry and wet conditions beyond the onset of grain crushing. A significant water weakening effect was observed in both rocks. A set of conventional triaxial experiments was also performed on both rocks under dry conditions at confining pressures ranging from 3 to 31 MPa. Microstructural observations were carried out on the deformed samples. The mechanical behavior of these high porosity carbonates is dominated by shear-enhanced compaction associated in most cases with strain hardening. Stress-induced cracking and grain crushing are the dominant micromechanisms of deformation in both rocks. In Majella grainstone, compactive shear bands appeared at low confinement, in qualitative agreement with the deformation bands observed in the field. At higher confining pressures, compaction localization was inhibited and homogeneous cataclastic flow developed. In Saint Maximin limestone, compaction localization was observed at all confining pressures. An increasing number of compactive shear bands at various orientations appeared with increasing strain. These new data suggest that compaction localization is important in the mechanical compaction of high porosity carbonates.     

Initiation and Propagation of Compaction Bands in Dry and Wet Bentheim Sandstone

We investigated initiation and propagation of compaction bands (CB) in six wet and four dry Bentheim sandstone samples deformed in axial compression tests with strain rates ranging from 3.2 × 10^?8 s^?1 to 3.2 × 10^?4 s^?1. Circumferential notches with 0.8-mm width and 5-mm depth served to initiate CB at mid-sample length. Wet samples were saturated with distilled water and deformed at 195 MPa confining pressure and 10 MPa pore pressure. Dry samples were deformed at 185 MPa confining pressure. Twelve P-wave sensors, eight S-wave sensors and two pairs of orthogonally oriented strain-gages were glued to the sample surface to monitor acoustic emission (AE), velocities and local strain during the loading process. Nucleation of compaction bands is indicated by AE clusters close to the notch tips. With progressive loading, AE activity increased and AE hypocenters indicated propagation of a single CB normal to the sample axis. CB propagation from the sample periphery towards the centre was monitored. Microstructural analysis of deformed samples shows excellent agreement between location of AE clusters and CBs. In both dry and wet samples the lateral propagation of CBs was about 100 times faster than axial shortening rates. At the slowest displacement rate, AE activity during band propagation was reduced and CB nucleation in wet samples occurred at 20% lower stresses. This may indicate an increasing contribution of stress corrosion processes to the formation of the compaction bands. In dry and wet samples inelastic compaction energy per area ranged between 16 and 80 kJ m^?2. This is in good agreement with previous estimates from laboratory and field studies.     

Interseismic Strain Localization in the San Jacinto Fault Zone

We investigate interseismic deformation across the San Jacinto fault at Anza, California where previous geodetic observations have indicated an anomalously high shear strain rate. We present an updated set of secular velocities from GPS and InSAR observations that reveal a 2–3 km wide shear zone deforming at a rate that exceeds the background strain rate by more than a factor of two. GPS occupations of an alignment array installed in 1990 across the fault trace at Anza allow us to rule out shallow creep as a possible contributor to the observed strain rate. Using a dislocation model in a heterogeneous elastic half space, we show that a reduction in shear modulus within the fault zone by a factor of 1.2–1.6 as imaged tomographically by Allam and Ben-Zion (Geophys J Int 190:1181–1196, 2012) can explain about 50 % of the observed anomalous strain rate. However, the best-fitting locking depth in this case (10.4 ± 1.3 km) is significantly less than the local depth extent of seismicity (14–18 km). We show that a deep fault zone with a shear modulus reduction of at least a factor of 2.4 would be required to explain fully the geodetic strain rate, assuming the locking depth is 15 km. Two alternative possibilities include fault creep at a substantial fraction of the long-term slip rate within the region of deep microseismicity, or a reduced yield strength within the upper fault zone leading to distributed plastic failure during the interseismic period.     

Formation of systems of incompact bands parallel to the compression axis in the unconsolidated sedimentary rocks: A model

Uniaxial compression of poorly lithified rocks leads to the formation of thin incompact layers (or bands, in the two-dimensional case) parallel to the compression axis, which are characterized by increased porosity. The standard model of the formation of such bands, as well as deformation bands of other types, associates them with the narrow zones of localization of plastic deformations. In the case of decompaction, it is assumed that transverse tensile deformations are localized within the band, which cause the band to dilate. Here, the formation of a band of localized deformations is treated as a loss-of-stability phenomenon. Based on observations, we propose a fundamentally different model of incompact bands formation, according to which the microdefects in sediment packing (pores) rather than the deformations are localized in the narrow zones. The localization of pores, which are initially randomly distributed in the medium, occurs as a result of their migration through the geomaterial. The migration and subsequent localization of pores are driven by a common mechanism, namely, a trend of a system to lower its total energy (small variations in total energy are equal to the increment of free energy minus the work of external forces). Migration of a single pore in a granular sedimentary rock is caused by the force f driving the defect. This force was introduced by J. Eshelby (1951; 1970). An important feature of our model is that the formation of an incompact band here does not have a sense of a loss of stability. Quite the contrary, the formation of incompact bands is treated as a gradual process spread over time. In this context, the origination of incompact band systems directly follows from our model itself, without any a priori assumptions postulating the existence of such systems and without any special tuning of the model parameters. Moreover, based on the proposed model, we can predict the incompact bands to always occur in the form of systems rather than as individual structures. A single incompact band may only be formed when the force resisting the pore motion, f _c , is absent.     

Characterization of Damage in Sandstones along the Mojave Section of the San Andreas Fault: Implications for the Shallow Extent of Damage Generation

Following theoretical calculations that suggest shallow generation of rock damage during an earthquake rupture, we measure the degree of fracture damage in young sedimentary rocks from the Juniper Hills Formation (JHF) that were displaced 21 km along the Mojave section of San Andreas Fault (SAF) and were not exhumed significantly during their displacement. In exposures adjacent to the fault, the JHF typically displays original sedimentary fabrics and little evidence of bulk shear strain at the mesoscopic scale. The formation is, however, pervasively fractured at the microscopic scale over a zone that is about a 100 m wide on the southwest side of the SAF near Little Rock. The abundance of open fractures, the poor consolidation, and the shallow inferred burial depth imply that the damage was generated close to the surface of the Earth. The spatial correlation of this damage with a seismically active trace of the SAF suggests that it was generated by SAF slip events that by assumption were of a seismic nature throughout the displacement history of the JHF. Thus the JHF provides a very shallow upper bound for the generation of brittle damage in a seismic fault zone. The fracture fabric is characterized by preferred orientations of fractures that split grains between contact points and is consistent with overall deformation under directed compression. However, the available results cannot be used to distinguish between proposed off-fault damage mechanisms. Fracture orientations are compatible with a maximum compressive stress oriented at a high angle to the fault at about 10 m, and at a lower, more variable angle farther away from the fault. The fracture distribution and fabric are consistent with observations made of the microscale damage characteristics of the Hungry Valley Formation in the northwestern section of the SAF in the Mojave, and with previous observations of exhumed, ancestral strands of the SAF.     

Crustal deformation before the 2008 Wenchuan MS8.0 earthquake studied using GPS data

We used GPS velocities from approximately 700 stations in western China to study the crustal deformation before the Wenchuan M_S8.0 earthquake. The processing methods included analyses of the strain rate field, inversion of fault locking and the GPS velocity profiles. The GPS strain rate in the E-W direction in the Qinghai-Tibet block shows that extensional deformation was dominant in the western region of the block (west of 92.5° E), while compressive deformation predominated in the eastern region of the block (from 92.5° E to 100° E). On a regional scale, the hypocentral region of the Wenchuan earthquake was located at the edge of an intense compression deformation zone of about 1.9 × 10⁻⁸/a in an east-west direction. The characteristic deformation in the seismogenic fault was compressive with a dextral component. The compression deformation rate was greater in the fault's western region than in its eastern region, and the strain accumulation was very slow on the fault scale. The results of a fault locking inversion show that the locking fraction and slip deficit was greater in the middle-northern section of the seismogenic fault than in the southern section. The GPS velocity profile before the Wenchuan earthquake shows that the compression deformation was smaller than the dextral deformation, which is asymmetrical with respect to the distribution of co-seismic displacement. These deformation characteristics should provide some clues to the Wenchuan earthquake which occurred in the later period of the earthquake cycle.     

蒙古-贝加尔裂谷区地壳应变场及其地球动力学涵义

蒙古-贝加尔地区是现今构造最活跃的大陆地区之一,其地壳构造运动及变形对我们理解大陆动力学问题具有重要的科学意义.基于融合的这一区域的GPS速度场,本文计算了该区应变率场和应变能变化率场.结果显示,蒙古褶皱带以南区域表现为NNE-SSW方向的压缩状态,主压应变率约为-2.0×10^-9/a,剪应变及面膨胀均较弱,表明蒙古褶皱带比较稳定.贝加尔裂谷整体处于拉张状态且伴有较强的剪应变和面膨胀,暗示可能有多种机制控制裂谷的张裂过程.蒙古高原西部有两条高应变率的构造带,结合深部存在地幔热柱等证据,我们认为这两条构造带及所围限的区域共同构成Amurian板块的西部边界—一条弥散变形的边界带.蒙古-贝加尔地区剪应变分布与0~40 km的地震活动性基本一致,表明该地区形变在地壳尺度耦合程度较高.地幔对流拖曳力场与主应变轴方向及应变率场的一致性表明,地幔对流可能是蒙古—贝加尔地区区域构造动力学过程主要控制因素之一.     

Shear localisation in upper mantle peridotites

Upper mantle peridotite bodies at the earth's surface contain relict structures and microstructures which provide direct information on the role and the mechanisms of shear localisation in the upper mantle. Deformation which occurred at high temperatures (T>950±50°C) is relatively homogeneous within domains ranging in scale from a few kilometres to a few tens of kilometres. Below 950±50°C strain is localised into centimetre to several hundred metre wide shear zones which commonly contain hydrated mylonitic peridotites. The microstructures developed in the peridotites suggest there is a correlation between the occurrence of shear localisation and the occurrence of strain softening and brittle deformation processes. The most important strain softening processes are inferred to be structural and reaction induced softening. Structural softening processes include dynamic recrystallisation and strain-induced transitions from dislocation creep to some form of grain-size-sensitive (GSS) creep. Reaction induced softening is related to the formation of fine grained polyphase reaction products which deform by GSS creep and the formation of weak sheet silicates such as phlogopite, chlorite, talc and antigorite. From experimental studies these softening processes and brittle deformation processes are inferred to occur mainly at temperatures less than about 910±160°C. This temperature range is inferred to be a significant rheological transition in the upper mantle. Below 910±160°C deformation during orogenesis may be accommodated by an anastomosing network of hydrated mylonitic shear zones with a distinct, perhaps weak, rheology. At higher temperatures strain is accommodated in much wider deformation zones.On the scale of the lithosphere the degree of localisation may be different to that determined at the scale of the periodotite massif. An anastomosing network of hundred metre wide mylonitic shear zones forming 0.05–0.3 by volume fraction of the mantle lithosphere atT<950°C could accommodate inhomogeneous or homogeneous bulk deformation depending on the spatial distribution and ordering of the mylonite zones. The higher temperature deformation at deeper levels in the mantle could be markedly inhomogeneous being concentrated in shear zones with widths in the range of 2–20 km, alternatively these zones may widen significantly during deformation, resulting in a decrease in the degree of localisation with increasing bulk strain.     

Structural Properties and Deformation Patterns of Evolving Strike-slip Faults: Numerical Simulations Incorporating Damage Rheology

We present results on evolving geometrical and material properties of large strike-slip fault zones and associated deformation fields, using 3-D numerical simulations in a rheologically-layered model with a seismogenic upper crust governed by a continuum brittle damage framework over a viscoelastic substrate. The damage healing parameters we employ are constrained using results of test models and geophysical observations of healing along active faults. The model simulations exhibit several results that are likely to have general applicability. The fault zones form initially as complex segmented structures and evolve overall with continuing deformation toward contiguous, simpler structures. Along relatively-straight mature segments, the models produce flower structures with depth consisting of a broad damage zone in the top few kilometers of the crust and highly localized damage at depth. The flower structures form during an early evolutionary stage of the fault system (before a total offset of about 0.05 to 0.1 km has accumulated), and persist as continued deformation localizes further along narrow slip zones. The tectonic strain at seismogenic depths is concentrated along the highly damaged cores of the main fault zones, although at shallow depths a small portion of the strain is accommodated over a broader region. This broader domain corresponds to shallow damage (or compliant) zones which have been identified in several seismic and geodetic studies of active faults. The models produce releasing stepovers between fault zone segments that are locations of ongoing interseismic deformation. Material within the fault stepovers remains damaged during the entire earthquake cycle (with significantly reduced rigidity and shear-wave velocity) to depths of 10 to 15 km. These persistent damage zones should be detectable by geophysical imaging studies and could have important implications for earthquake dynamics and seismic hazard.     

Geometry of the Nojima Fault at Nojima-Hirabayashi, Japan – I. A Simple Damage Structure Inferred from Borehole Core Permeability

The 1995 Kobe (Hyogo-ken Nanbu) earthquake, M = 7.2, ruptured the Nojima fault in southwest Japan. We have studied core samples taken from two scientific drillholes that crossed the fault zone SW of the epicentral region on Awaji Island. The shallower hole, drilled by the Geological Survey of Japan (GSJ), was started 75 m to the SE of the surface trace of the Nojima fault and crossed the fault at a depth of 624 m. A deeper hole, drilled by the National Research Institute for Earth Science and Disaster Prevention (NIED) was started 302 m to the SE of the fault and crossed fault strands below a depth of 1140 m. We have measured strength and matrix permeability of core samples taken from these two drillholes. We find a strong correlation between permeability and proximity to the fault zone shear axes. The half-width of the high permeability zone (approximately 15 to 25 m) is in good agreement with the fault zone width inferred from trapped seismic wave analysis and other evidence. The fault zone core or shear axis contains clays with permeabilities of approximately 0.1 to 1 microdarcy at 50 MPa effective confining pressure (10 to 30 microdarcy at in situ pressures). Within a few meters of the fault zone core, the rock is highly fractured but has sustained little net shear. Matrix permeability of this zone is approximately 30 to 60 microdarcy at 50 MPa effective confining pressure (300 to 1000 microdarcy at in situ pressures). Outside this damage zone, matrix permeability drops below 0.01 microdarcy. The clay-rich core material has the lowest strength with a coefficient of friction of approximately 0.55. Shear strength increases with distance from the shear axis. These permeability and strength observations reveal a simple fault zone structure with a relatively weak fine-grained core surrounded by a damage zone of fractured rock. In this case, the damage zone will act as a high-permeability conduit for vertical and horizontal flow in the plane of the fault. The fine-grained core region, however, will impede fluid flow across the fault.     

Current horizontal strain field in Chinese mainland derived from GPS data

Introduction In the years when the reliable data could not be obtained and in the analysis of strain property and magnitude in history, the intensity, property and activity pattern of strain field were mainly inferred on the bases of geometric characters of surface traces and behaviors (especially the faults) as well as the characteristics of petrology (XIE, et al, 1993; Molnar, Tapponnier, 1975, 1977; Tapponnier, Molnar, 1977; FU, et al, 2000). However, they are the averaged results accumu…     

Spatial variation of crustal strain in the Kachchh region,India: Implication on the Bhuj earthquake of 2001

The Kachchh province of Western India is a major seismic domain in an intraplate set-up. This seismic zone is located in a rift basin, which was developed during the early Jurassic break-up of the Gondwanaland. The crustal strain determined from the GPS velocity data of post-seismic time period following the 2001 Bhuj earthquake indicates a maximum strain rate of ∼266 × 10⁻⁹ per year along N013°. Focal mechanism solutions of the main event of 26 January 2001 and the aftershocks show that the maximum principal stress axis is close to this high strain direction. Maximum shear strain rate determined from the GPS data of the area has similar orientation. The unusually high strain rate is comparable in magnitude to the continental rift systems. The partitioning of the regional NE–SW horizontal stress (S_Hmax) by the pre-existing EW-striking boundary fault developed the strike–slip components parallel to the regional faults, the normal components perpendicular to the faults, NE-striking conjugate Riedel shear fractures and tension fractures. The partitioned normal component of the stress is considered to be the major cause for compression across the regional EW faults and development of the second-order conjugate shear fractures striking NE–SW and NW–SE. The NE-striking transverse faults parallel to the anti-Riedel shear planes have become critical under these conditions. These anti-Riedel planes are interpreted to be critical for the seismicity of the Kachchh region. The high strain rate in this area of low to moderate surface heat flow is responsible for deeper position of the brittle–ductile transition and development of deep seated seismic events in this intraplate region.     

Development of composite planar fabric in mylonites along the Median Tectonic Line, southwest Japan

Abstract Mylonites along the Median Tectonic Line, southwest Japan commonly contain shear bands comprising S(-C)-Ss fabrics. This paper stresses the lithologic control on the orientation, dimension and development of shear bands by comparing the microstructure of the shear bands in different rock types (P mylonites, F mylonites, micaceous phyllonite and quartzose phyllonite). There is no significant change of the α angles (average 21–24°) between Ss and S toward the centre of the shear zone (viz. increasing the intensity of mylonitization) and it is different from the S-C relationship in a narrow sense.
The generation of the composite planar fabric can be classified into four different strain partitioning models: S only type without any slip surface (model A); S-C type (model B); S-Ss type with Ss-slip precedence (model C), and S-Ss type with S-slip precedence (model D). Model C is proposed in this paper and is similar to the model for the generation of Riedel shears in brittle shear zones. An unstable slip between porphyroclasts and the matrix during ductile flow can easily initiate shear bands. Formation of a composite planar fabric is initiated according to model A, followed by model C in conditions of increasing strain, and then model D when the angle between S and the shear zone boundary becomes small enough (α/2 = 10°) to produce S-slip. Thus the generation of the shear bands probably begins in the early stages of shear deformation and continues until the latest stages.     

Failure Mode and Spatial Distribution of Damage in Rothbach Sandstone in the Brittle-ductile Transition

— To elucidate the spatial complexity of damage and evolution of localized failure in the transitional regime from brittle faulting to cataclastic ductile flow in a porous sandstone, we performed a series of triaxial compression experiments on Rothbach sandstone (20% porosity). Quantitative microstructural analysis and X-ray computed tomography (CT) imaging were conducted on deformed samples. Localized failure was observed in samples at effective pressures ranging from 5 MPa to 130 MPa. In the brittle faulting regime, dilating shear bands were observed. The CT images and stereological measurements reveal the geometric complexity and spatial heterogeneity of damage in the failed samples. In the transitional regime (at effective pressures between 45 MPa and 130 MPa), compacting shear bands at high angles and compaction bands perpendicular to the maximum compression direction were observed. The laboratory results suggest that these complex localized features can be pervasive in sandstone formations, not just limited to the very porous aeolian sandstone in which they were first documented. The microstructural observations are in qualitative agreement with theoretical predictions of bifurcation analyses, except for the occurrence of compaction bands in the sample deformed at effective pressure of 130 MPa. The bifurcation analysis with the constitutive model used in this paper is nonadequate to predict compaction band formation, may be due to the neglect of bedding anisotropy of the rock and multiple yield mechanisms in the constitutive model.     

GPS-based crustal deformations in Azerbaijan and their influence on seismicity and mud volcanism

Using Shen’s method (Shen et al., 1996), deformations of the Earth’s crust in Azerbaijan were studied based on GPS measurements. For estimating the rate of deformation, we used the field of velocity vectors for Azerbaijan, Iran, Georgia, and Armenia that were derived from GPS measurements during 1998–2012. It is established that compression is observable along the Greater Caucasus, in Gobustan, the Kura depression, Nakhchyvan Autonomous Republic, and adjacent areas of Iran. The axes of compression/contraction of the crust in the Greater Caucasus region are oriented in the S-NE direction. The maximum strain rate (approximately 200 × 10^?9 per annum) is documented in the zone of mud volcanism at the SHIK site (Shykhlar), which is marked by a sharp change in the direction of the compression axes (SW-NE). It is revealed that the deformation field also includes the zones where strain rates are very low approximating 5 × 10^?9 per annum. These zones include the Caspian-Guba and northern Gobustan areas, characterized by extensive development of mud volcanism. The extension zones are confined to the Lesser Caucasus and are revealed in the Gedabek (GEDA) and Shusha (SHOU) areas, as well as in the zone located between the DAMO and PIRM sites (Iran), where the deformation rate amounts to 100 × 10^?9 per annum. It is concluded that the predominant factor responsible for the eruption of mud volcanoes is the intensity of gas-generation processes in the earth’s interior, while deformation processes play the role of a trigger. The zone of the epicenters of strong earthquakes is correlated to the gradient zone in the crustal strain rates.     

Rheology of arc dacite lavas: experimental determination at low strain rates

Andesitic–dacitic volcanoes exhibit a large variety of eruption styles, including explosive eruptions, endogenous and exogenous dome growth, and kilometer-long lava flows. The rheology of these lavas can be investigated through field observations of flow and dome morphology, but this approach integrates the properties of lava over a wide range of temperatures. Another approach is through laboratory experiments; however, previous studies have used higher shear stresses and strain rates than are appropriate to lava flows. We measured the apparent viscosity of several lavas from Santiaguito and Bezymianny volcanoes by uniaxial compression, between 1,109 and 1,315?K, at low shear stress (0.085 to 0.42?MPa), low strain rate (between 1.1?×?10^?8 and 1.9?×?10^?5?s^?1), and up to 43.7 % total deformation. The results show a strong variability of the apparent viscosity between different samples, which can be ascribed to differences in initial porosity and crystallinity. Deformation occurs primarily by compaction, with some cracking and/or vesicle coalescence. Our experiments yield apparent viscosities more than 1 order of magnitude lower than predicted by models based on experiments at higher strain rates. At lava flow conditions, no evidence of a yield strength is observed, and the apparent viscosity is best approached by a strain rate- and temperature-dependent power law equation. The best fit for Santiaguito lava, for temperatures between 1,164 and 1,226?K and strain rates lower than 1.8?×?10^?4?s^?1, is $ \log {\eta_{\text{app}}} = - 0.738 + 9.24 \times {10^3}{/}T(K) - 0.654 \cdot \log \dot{\varepsilon } $ where η _app is apparent viscosity and $ \dot{\varepsilon } $ is strain rate. This equation also reproduced 45 data for a sample from Bezymianny with a root mean square deviation of 0.19 log unit Pa?s. Applying the rheological model to lava flow conditions at Santiaguito yields calculated apparent viscosities that are in reasonable agreement with field observations and suggests that internal shear heating may be significant ongoing heat source within these flows, enabling highly viscous lava to travel long distances.