花岗岩非稳态流变实验

脆塑性转化带对于研究岩石圈变形、断层强度和变形机制以及强震的孕育和发生具有重要意义.文中采用汶川地震震源区彭灌杂岩中具有代表性的细粒花岗岩样品,在固体压力介质三轴实验系统上开展了高温高压非稳态流变实验研究.实验设计模拟了汶川地震区地壳10～30km深度的实际温度和压力,温度为190～490℃,压力为250～750MPa...     

The distribution of earthquakes with depth and stress in subducting slabs

We explain the global variation of Benioff zone seismicity with depth and the orientation of stress axes of deep and intermediate earthquakes using numerical models of subducting slabs. Models that match the seismicity and stress require a barrier to flow at the 670 km seismic discontinuity. The barrier may be a viscosity increase of at least an order of magnitude or a chemical discontinuity. Instantaneous flow is subparallel to the slabs for models with a viscosity increase but contorted for models with a chemical barrier. Log N (number of earthquakes) decreases linearly to 250–300 km depth and increases thereafter. Stress magnitude in our models shows the same pattern, in accord with experiments showing N proportional to e^kσ, with k a constant and σ stress magnitude. The models predict downdip compression in the slabs at depths below 300–400 km, as observed for earthquake stress axes.     

Fluid Pressure Variation in a Sedimentary Geothermal Reservoir in the North German Basin: Case Study Groß Schönebeck

The Rotliegend of the North German basin is the target reservoir of an interdisciplinary investigation program to develop a technology for the generation of geothermal electricity from low-enthalpy reservoirs. An in situ downhole laboratory was established in the 4.3 km deep well Groβ Schönebeck with the purpose of developing appropriate stimulation methods to increase permeability of deep aquifers by enhancing or creating secondary porosity and flow paths. The goal is to learn how to enhance the inflow performance of a well from a variety of rock types in low permeable geothermal reservoirs. A change in effective stress due to fluid pressure was observed to be one of the key parameters influencing flow properties both downhole and in laboratory experiments on reservoir rocks. Fluid pressure variation was induced using proppant-gel-frac techniques as well as waterfrac techniques in several different new experiments in the borehole. A pressure step test indicates generation and extension of multiple fractures with closure pressures between 6 and 8.4 MPa above formation pressure. In a 24-hour production test 859 m³ water was produced from depth indicating an increase of productivity in comparison with former tests. Different depth sections and transmissibility values were observed in the borehole depending on fluid pressure. In addition, laboratory experiments were performed on core samples from the sandstone reservoir under uniaxial strain conditions, i.e., no lateral strain, constant axial load. The experiments on the borehole and the laboratory scale were realized on the same rock types under comparable stress conditions with similar pore pressure variations. Nevertheless, stress dependences of permeability are not easy to compare from scale to scale. Laboratory investigations reflect permeability variations due to microstructural heterogeneities and the behavior in the borehole is dominated by the generation of connections to large-scale structural patterns.     

Mantle rheology: radial and lateral viscosity variations inferred from microphysical creep laws

Inferences on the rheology of the mantle based on theoretical and experimental rate equations for steady state creep are discussed and compared with results from geophysical models. The radial increase of viscosity by one to three orders of magnitude across the mantle, required by inversion of postglacial rebound and geodynamic data, is confirmed by microphysical models based on the estimation of continuous and discontinuous changes of creep parameters with depth. The upper mantle (viscosity 10²⁰–10²¹ Pa s) is likely to show non-Newtonian rheology (power-law creep) for average grain sizes larger than 0.1 mm as an order of magnitude. Given the variability of both grain size and stress conditions, local regions of linear rheology can be present. The rheology of transition zone and lower mantle (viscosity 10²²–10²⁴ Pa s) cannot be definitely resolved at present. Estimation of creep parameters leads to possible nonlinear or mixed rheology, if grain sizes are not lower than 0.1 mm and flow conditions can be approximated by a constant strain rate of about 10⁻¹⁵ s⁻¹. This conclusion can be modified by different flow conditions (e.g. a decrease in strain rate or constant viscous dissipation). Furthermore, experiments on fine-grained garnetites and perovskite analogues have shown that diffusion creep is predominant at laboratory conditions. However, the pressure dependence of creep in these phases is unknown, and therefore direct extrapolation to lower mantle conditions is necessarily speculative. Lateral variations of viscosity, largest in the upper and lowermost mantle (up to 2–4 orders of magnitude) are predicted by models based on lateral temperature anomalies derived from seismic tomographic models.     

Subduction zone earthquakes and stress in slabs

Summary The pattern of seismicity as a function of depth in the world, and the orientation of stress axes of deep and intermediate earthquakes, are explained using viscous fluid models of subducting slabs, with a barrier in the mantle at 670 km. 670 km is the depth of a seismic discontinuity, and also the depth below which earthquakes do not occur. The barrier in the models can be a viscosity increase of an order of magnitude or more, or a chemical discontinuity where vertical velocity is zero. LongN versus depth, whereN is the number of earthquakes, shows (1) a linear decrease to about 250–300 km depth, (2) a minimum near that depth, and (3) an increase thereafter. Stress magnitude in a subducting slab versus depth, for a wide variety of models, shows the same pattern. Since there is some experimental evidence thatN is proportional toe , where is a constant and is the stress magnitude, the agreement is encouraging. In addition, the models predict down-dip compression in the slab at depths below 400 km. This has been observed in earlier studies of earthquake stress axes, and we have confirmed it via a survey of events occurring since 1977 which have been analysed by moment tensor inversion. At intermediate depths, the models predict an approximate but not precise state of down-dip tension when the slab is dipping. The observations do not show an unambiguous state of down-dip tension at intermediate depths, but in the majority of regions the state of stress is decidedly closer to down-dip tension than it is to down-dip compression. Chemical discontinuities above 670 km, or phase transitions with an elevation of the boundary in the slab, predict, when incorporated into the models, stress peaks which are not mirrored in the profile of seismicity versus depth. Models with an asthenosphere and mesosphere of appropriate viscosity can not only explain the state of stress observed in double Benioff zones, but also yield stress magnitude profiles consistent with observed seismicity. Models where a nonlinear rheology is used are qualitatively consistent with the linear models.     

The storage and release of magma on Mount Etna

The eruptive history of Etna during the past 450 years provides data on effusion rates, volumes of magma involved, and the nature of the eruptive conduits. These data are interpreted in terms of a two-part intravolcanic magma reservoir which feeds the flank eruptions through dike-like conduits. The structural framework of the volcano which controls the spatial distribution of eruptive sites is partly inherited from the basement and partly controlled by the central magma column and the surrounding caldera boundary faults. Hydraulic fracturing theory predicts that the central magma column will fail at depths below 1 km if the tensile strength of the conduit rocks is about 100 bars and that a peak fracturing capability will be reached between 1 and 2 km depth. This inference agrees well with the peak of flank eruptive activity at 1.4 km below the summit observed in the data on the loci of eruptions. The average flank-eruption feeding dike is defined and shown to be capable of the observed maximum effusion rates (20–100 m³ s⁻¹) from magmatic pressure differences of 30–150 bars     

亚稳态橄榄石对俯冲带负浮力的影响及其动力学意义

利用有限元方法,计算了不同俯冲速度及热传导系数俯冲岩石层的负浮力.在h-lt;400km和h-gt;740km的深度范围,低温高密俯冲带的负浮力随深度单调增加.因为尖晶石相到后尖晶石相的相变有负的克兰帕龙斜率,俯冲带冷的物质在660km间断面以下不到100km的深度范围内仍以低密度的低压物质相存在.所以在该深度范围,俯冲带受到了周围地幔阻止其插入下地幔的浮力作用.在400km-lt;h-lt;660km深度范围,由于受橄榄石相变的影响,不同计算模型负浮力随深度的变化有明显的不同.对于可能的相变动力学模型,亚稳态橄榄石的存在使负浮力随深度的增加值减小.其作用是不利于俯冲带直接穿透660km间断面,并引起该深度范围俯冲带沿俯冲方向压应力分布的变化.     

Space-time variations in the parameters of split S waves excited by local earthquakes in southern Kamchatka

An analysis of the distribution of split S-wave parameters due to small local earthquakes occurring down to a depth of 200 km beneath the PET IRIS station (in the town of Petropavlovsk-Kamchatskii) for the period 1993–2002 revealed an inhomogeneous distribution of anisotropic properties and time-dependent changes in the state of stress and strain in the sinking slab. Three time intervals have been identified from the behavior of the split-wave parameters: 1993–1995, 1996–1998, and 1999–2002. The most orderly orientation of fast-wave polarization azimuths at all depths and lower stresses in the earth was observed in 1996–1998. Regions of relatively stable and unstable parameter behavior consistent with increased and decreased P and S velocities have been identified. The central block down to a depth of 120 km and the lower part of the Benioff zone are classified as stable regions dominated by E-SE fast-wave azimuths. Unstable behavior of fast-wave azimuths was noticed for depths of 60 to 90 km and around the top of the Benioff zone. The highest values of S-wave delay time and regions where the delay times were migrating mostly occur at contacts between more rigid and weakened blocks.     

Transient creep and semibrittle behavior of crystalline rocks

We review transient creep and semibrittle behavior of crystalline solids. The results are expected to be pertinent to crystalline rocks undergoing deformation in the depth range 5 to 20 km, corresponding to depths of focus of many major earthquakes. Transient creep data for crystalline rocks at elevated temperatures are analyzed but are poorly understood because of lack of information on the deformation processes which, at low to moderate pressure, are likely to be semibrittle in nature. Activation energies for transient creep at high effective confining pressure are much higher than those found for atmospheric pressure tests in which thermally-activated microfracturing probably dominates the creep rate. Empirical transient creep equations are extrapolated at 200° to 600°C, stresses from 0.1 to 1.0 kbar, to times ranging from 3.17×10² to 3.17×10⁸ years. At the higher temperatures, appreciable transient creep strains may take place but the physical significance of the results is in question because the flow mechanisms have not been determined. The purpose of this paper is to stimulate careful research on this important topic.     

Infiltration,runoff and sediment production in blanket peat catchments: implications of field rainfall simulation experiments

Blanket peat covers the headwaters of many major European rivers. Runoff production in upland blanket peat catchments is flashy with large flood peaks and short lag times; there is minimal baseflow. Little is known about the exact processes of infiltration and runoff generation within these upland headwaters. This paper presents results from a set of rainfall simulation experiments performed on the blanket peat moorland of the North Pennines, UK. Rainfall was simulated at low intensities (3–12 mm h^?1), typical of natural rainfall, on bare and vegetated peat surfaces. Runoff response shows that infiltration rate increases with rainfall intensity; the use of low‐intensity rainfall therefore allows a more realistic evaluation of infiltration rates and flow processes than previous studies. Overland flow is shown to be common on both vegetated and bare peat surfaces although surface cover does exert some control. Most runoff is produced within the top few centimetres of the peat and runoff response decreases rapidly with depth. Little vertical percolation takes place to depths greater than 10 cm owing to the saturation of the peat mass. This study provides evidence that the quickflow response of upland blanket peat catchments is a result of saturation‐excess overland flow generation. Rainfall–runoff response from small plots varies with season. Following warm, dry weather, rainfall tends to infiltrate more readily into blanket peat, not just initially but to the extent that steady‐state surface runoff rates are reduced and more flow takes place within the peat, albeit at shallow depth. Sediment erosion from bare peat plots tends to be supply limited. Seasonal weather conditions may affect this in that after a warm, dry spell, surface desiccation allows sediment erosion to become transport limited. Copyright © 2002 John Wiley & Sons, Ltd.     

Source modeling of the Hsingtai,China earthquakes of March 1966

The Hsingtai, China earthquakes of March 1966 were a series of destructive earthquakes associated with the Shu-lu graben. Five strong shocks of M_s ≥ 6 occurred within a period of less than a month, the largest of which was M_s 7.2. Body and surface waves over the period range from several to 100 s have been modeled for the four largest events using synthetic seismograms in the time domain and spectral analysis in the frequency domain. Data from ground deformation, local geology, regional seismic network, and teleseismic joint epicenter determination have also been used to constrain the source model and the rupture process.The fault mechanism of the Hsingtai sequence was mainly strike-slip with a small component of normal dip-slip. The strikes of the four largest shocks range from ～ N26° to 30°E, approximately along strike of the major faults of the Shu-lu graben and the aftershock distribution. The source mechanisms can be explained with a NNW-SSE extensional stress and a NEE-SWW compressional stress acting in the area. The major shocks all had focal depths ～ 10 km.The four largest shocks in the sequence were characterized by a relatively simple and smooth dislocation time history. The durations of the far-field source time functions ranged from 3.5 to 5 s, while the rise times were all ～ 1 s. The seismic moments of the four largest earthquakes ranged from 1.43 × 10²⁵ to 1.51 × 10²⁶ dyne cm^?1. The fault sizes of the four events were very close. Assuming circular faults, the diameters of the four events were determined to be between 10 and 14 km. Stress drops varied from ～ 52 to 194 bars. A trend of increasing stress drop with earthquake size was observed.A survey of stress drop determinations for 15 major intraplate earthquakes shows that on the average the magnitude of stress drop of oceanic intraplate earthquakes and passive continental margin events is higher (～ 200 to several hundred bars) than that of continental intraplate earthquakes (～ 100 bars or less).     

日本海及中国东北地震的深度分布及其应力状态

本文分析了日本海及中国东北的地震深度分布。证实了日本本州北部至中国东北的贝尼奥夫带(Benioff)基本是连续的,该带的倾向约为北85°西,倾角约为29°,深度在150公里以下贝尼奥夫带厚度约为20公里。研究了日本本州北部至中国东北的震级M_b≥5.0地震的震源机制解,发现中国东北地壳应力场与日本海地壳的应力场方向一致,来源于太平洋板块的挤压。在俯冲带上,深度在100公里到200公里之间的情况较为复杂,大多数地震显示的主压应力方向与贝尼奥夫带的倾向、倾角一致,有的T轴取向与贝尼奥夫带的倾向、倾角一致,有的特征方向与贝尼奥夫带倾向、倾角均不一致。深度在200公里至500公里之间,主压应力方向近于水平,并与贝尼奥夫带走向垂直,张应力轴相对集中。深度大于500公里时,主压应力方向与贝尼奥夫带的倾向、倾角一致,张应力轴相对集中     

帕米尔—兴都库什地区构造应力场反演及拆离板片应力形因子特征研究

从Global CMT目录搜集了1976年1月至2016年6月之间的震源深度大于70km的255个震源机制解,用阻尼应力反演方法,分70～160km和170～310km两个深度,计算了帕米尔—兴都库什地区的构造应力场;同时以10km为间隔计算了兴都库什地区深度介于70～310km之间的应力形因子.得到以下初步结论:兴都库什板片向下俯冲和帕米尔地区断裂带的横向拉张,可能是导致应力场不同的原因.兴都库什俯冲带与帕米尔俯冲带碰撞,导致交汇地区(37°N—37.5°N)的应力场参数突变.兴都库什俯冲板片受到深部温度、压力等因素,出现薄弱面进而形成拆离板片.其脱离了主俯冲板片的束缚后,重力的上下拉张作用导致空区附近张轴倾伏角接近90°,拆离板片俯冲至上地幔不连续面,导致板片部分熔融进而应力形因子随着深度变小.而拆离板片受到地幔挤压其内部发生破碎,其压应力轴由西部的NS到东部NW-SE方向偏转及纵向张应力轴倾伏角变小.     

Variation of seismic source parameters and stress drops within a descending slab and its implications in plate mechanics

A least-squares searching technique has been developed to estimate the source dimensions of intermediate and deep focus earthquakes using azimuthal variations of body wave pulse-widths. With this method and also amplitude data, modes of rupture propagation, seismic moments, and stress drops of 17 intermediate and deep focus earthquakes in the Tonga-Kermadec region have been determined in order to investigate variations in source properties and the state of stress within the descending slab there. Three different modes of rupture; unilateral, bilateral, and circular faults, are compared and tested against observations. Results indicate that the unilateral fault is the best model for most of the earthquakes studied. Stress drops of the 17 events vary within a very large range, from 20 bar to about 4.6 kbar, and change significantly with depth. The magnitude of stress drops for depths between 220 and 430 km is very much higher than at shallower depths. This change in stress drop magnitude at about 220 km-depth seems to reflect a change in material properties both in the mantle and within the slab. Two regions of high stress drop are observed at depths of about 360 and 640 km. A relative minimum of stress drop is found at about 450–560 km where the earthquake frequency is particularly high. Earthquakes at the northern end to the Tonga arc, where the Benioff zone is laterally bent, show systematically higher stress drops than other events at comparable depths, but away from the bend. Also, events in regions of low seismicity appear to have higher stress drops than those in regions of high seismicity. The upper bound of seismic efficiency is found to decrease with depth, implying an increase of frictional force with depth at the earthquake source.     

Rheology of the lower crust beneath the northern part of North China: Inferences from lower crustal xenoliths from Hannuoba basalts, Hebei Province, China

Lower crustal xenoliths brought up rapidly by basaltic magma onto the earth surface may provide di-rect information on the lower crust. The main purpose of this research is to gain an insight into the rheology of the lower crust through the detailed study of lower crustal xenoliths collected from the Hannuoba basalt, North China. The lower crustal xenoliths in this area consist mainly of two pyroxene granulite, garnet granulite, and light-colored granulite, with a few exception of felsic granulite. The equilibration temperature and pressure of these xenoliths are estimated by using geothermometers and geobarometers suitable for lower crustal xenoliths. The obtained results show that the equilibration temperature of these xenoliths is within the range of 785―900℃, and the equilibrium pressure is within the range of 0.8―1.2 GPa, corresponding to a depth range of 28―42 km. These results have been used to modify the previously constructed lower crust-upper mantle geotherm for the studied area. The dif-ferential stress during the deformation process of the lower crustal xenoliths is estimated by using recrystallized grain-size paleo-piezometer to be in the range of 14―20 MPa. Comparing the available steady state flow laws for lower crustal rocks, it is confirmed that the flow law proposed by Wilks et al. in 1990 is applicable to the lower crustal xenoliths studied in this paper. The strain rate of the lower crust estimated by using this flow law is within the range of 10-13―10-11 s-1, higher than the strain rate of the upper mantle estimated previously for the studied area (10-17―10-13 s-1); the equivalent viscosity is estimated to be within the range of 1017―1019Pa·s, lower than that of the upper mantle (1019―1021 Pa·s). The constructed rheological profiles of the lower crust indicate that the differential stress shows no significant linear relation with depth, while the strain rate increases with depth and equivalent vis-cosity decrease with depth. The results support the viewpoint of weak lower continental crust.     

滇西北地区的现代构造应力场

通过区内活断层性质、盆地生成和水系分布极向等的研究,进一步从地质上论证了该区北北西—近南北向的现代构造应力场。认为这一局部应力场是川滇菱块作南向滑移引起的地壳上部的构造变形事件。晚新生代以来的断陷盆地是这一滑移过程造成的地壳表层的拉分盆地。在地壳的下部有可能受到大区域板块北(偏东)向推挤作用的影响。中小地震和强震在深度分布、破裂特征及应力方向上的区别有可能是这种不同深度层次上应力差异作用的某些信息     

Curie point depth from spectral analysis of aeromagnetic data from Cerro Prieto geothermal area,Baja California,México

Using aeromagnetic data acquired in the area from the Cerro Prieto geothermal field, we estimated the depth to the Curie point isotherm, interpreted as the base of the magnetic sources, following statistical spectral-based techniques. According to our results the Curie point isotherm is located at a depths ranging from 14 to 17 km. Our result is somewhat deeper than that obtained previously based only in 2-D and 3-D forward modeling of previous low-quality data. However, our results are supported by independent information comprising geothermal gradients, seismicity distribution in the crust, and gravity determined crustal thickness. Our results imply a high thermal gradient (ranging between 33 and 38 °C/km) and high heat flow (of about 100 mW/m²) for the study area. The thermal regime for the area is inferred to be similar to that from the Salton trough.     

In-situ fluid-pressure measurements for earthquake prediction: An example from a deep well at Hi Vista,California

Short-term earthquake prediction requires sensitive instruments for measuring the small anomalous changes in stress and strain that precede earthquakes. Instruments installed at or near the surface have proven too noisy for measuring anomalies of the size expected to occur, and it is now recognized that even to have the possibility of a reliable earthquake-prediction system will require instruments installed in drill holes at depths sufficient to reduce the background noise to a level below that of the expected premonitory signals. We are conducting experiments to determine the maximum signal-to-noise improvement that can be obtained in drill holes. In a 592 m well in the Mojave Desert near Hi Vista, California, we measured water-level changes with amplitudes greater than 10 cm, induced by earth tides. By removing the effects of barometric pressure and the stress related to earth tides, we have achieved a sensitivity to volumetric strain rates of 10^–9 to 10^–10 per day. Further improvement may be possible, and it appears that a successful earthquake-prediction capability may be achieved with an array of instruments installed in drill holes at depths of about 1 km, assuming that the premonitory strain signals are, in fact, present.     

日本海-鄂霍次克海下俯冲带的应力状态及其深部形变

通过地震分布及地震机制解所反映的日本海-鄂霍次克海俯冲带的形态及应力状态,研究了俯冲带深部形变及650km间断面的穿透问题.日本海Benioff带较直,连续性较好;鄂霍次克海Benioff带弯度稍大,220-320km深度之间地震很少.两俯冲带在浅部及深部地震密集,100-200km深度之间有双地震层.应力状态随深度变化,200km深度以下P,T轴方向相对集中,P轴接近俯冲方向,在约100-200km深度附近,P,T轴均接近俯冲方向.观测和理论地震图拟合分析表明,地震断层面走向接近俯冲带走向,断裂的结果使俯冲带在深部倾角变小.     

Deformation mechanism maps for polycrystalline galena

A deformation mechanism map, depicting the fields of stress and temperature in which modes of plastic flow are dominant (i.e. provide the fastest strain rate) has been calculated for polycrystalline, stoichiometric galena for each of two grain sizes; namely, 10 and 10³ μm. The deformation mechanisms considered were dislocation glide, dislocation creep (i.e. creep involving dislocation climb), Nabarro-Herring creep and Coble creep.During folding and related tectonic deformation in the earth's crust steady-state flow of galena may occur by either dislocation or diffusion creep at very low differential stresses (typically ranging from 100 bars down to 10^?2 bar, or less). The dislocation creep field will be enlarged at the expense of that for diffusion creep, however, if (a) the stress dependence of strain rate for dislocation creep decreases at low stresses and if the grain size is greater than 10 μm, or (b) diffusion creep rates decline at high strains due to the presence of second-phase particles in the grain boundaries. It is probable that dislocation glide will be the dominant deformation mechanism in galena only at low temperatures and fast strain rates and it is unlikely to be the dominant mechanism during folding. Kinking may have some potential as an indicator of these temperatures and strain rates.