台湾集集地震近场地震动的上盘效应

1999年9月21日（当地时间）台湾集集7.6级地震是一个逆断层型地震．用回归分析法对台湾集集地震的加速度峰值数据进行分析，得出了这次地震的水平与垂直向的加速度峰值衰减关系．从残差分布上看，位于断层上盘和下盘上的加速度峰值与从衰减关系所得到的结果相比存在不同的系统偏差，断层上盘地表的加速度峰值较高，而下盘地表的加速度峰值较低．从这次地震的加速度峰值分布等值线图上也可以看出，加速度峰值的分布相对于断层呈现明显的不对称性，上盘衰减较慢而下盘衰减较快．在近断层强地面运动研究、地震危险性分析、设定地震研究与震害预测等工作中，应考虑可能地震的震源机制特点，以便使所用的衰减模型更能反映不同地震环境地区的地震动分布特征．      

龙门山断裂带金河磷矿剖面断层泥的低速至高速摩擦性质研究

本文对龙门山断裂带金河磷矿浅钻岩芯中的三种断层泥开展了低速到高速摩擦滑动的实验研究,并对实验变形样品开展了BET比表面积研究.摩擦实验在干燥和孔隙水压条件下开展,速率范围涵盖20 μm·s^-1~1.4 m·s^-1.实验结果显示,三种断层泥在干燥条件下的摩擦性质差别不大,但在孔隙水压条件下,三者的中低速摩擦强度与层状硅酸盐矿物的种类而非总含量紧密相关,蒙脱石和伊利石相比绿泥石更能有效地弱化断层.三种断层泥在孔隙水压条件下存在中低速率域的速度强化,暗示着对断层的加速滑动存在一定的阻碍作用.孔隙水压下,黄绿色和灰绿色断层泥的初始动态弱化非常迅速并伴随断层泥层的瞬时扩容,凹凸体急剧加热导致的局部热压作用可能是造成这种力学行为的物理机制.在经历高速滑动之后,三种断层泥在干、湿条件下的BET比表面积都显著降低,暗示着可能发生了颗粒烧结.中低速域内,孔隙水的存在使得断层泥呈现分散式的剪切变形,BET比表面积的增加因此比干燥条件下更加明显.对表面能的估算表明,颗粒磨碎所消耗的能量至多不超过摩擦力做功的8%,暗示着断层作用中颗粒磨碎所占的能量比例较低.     

Review and new interpretation for the propagation characteristics associated with the 1999 Chi‐Chi earthquake faulting event

The Chi‐Chi earthquake (M_W = 7.6) took place in central western Taiwan in 1999. The earthquake caused reactivation of the Chelungpu Fault and resulted in 100‐km‐long surface ruptures. The fault strikes mostly north–south to NNE–SSW; however, the northern tip of the southern segment of the surface ruptures rotates clockwise to define an east–west trend, then jumps to a shorter NNW‐trending rupture. The largest vertical displacement is recorded in the Shihkang area of the Shihkang–Shangchi Fault Zone, where vertical slips are up to 8–10 m. The Shihkang–Shangchi Fault Zone displays a complex fault pattern as a linkage damage zone between two fault segments with the greatest concentration of faults and fractures. Our new interpretation, based on recent published geometric, kinematic, and geophysical studies on the Chi‐Chi earthquake fault, suggests that the Shihkang–Shangchi Fault Zone is not a simple termination zone, but may be an ‘overstep zone’ or a ‘transfer zone’. Slip analysis along the surface ruptures indicates that they are composed of three fault segments and the amount of slip partly depends on the intersection angle between slip direction and fault strike. Our numerical modeling for the area indicates that Coulomb stress changes are mainly concentrated on tips and bends of the surface ruptures. Slip patterns indicate that the fault propagates toward the northeast. Therefore, this study suggests high potential for future earthquake activity along the unruptured Shangchi segment. Hence, future geohazard studies should focus on the Shangchi segment to evaluate potential earthquakes, determine recurrence intervals, and reduce future earthquake hazards.     

Pore-fluid pressures and frictional heating on a fault surface

This study considers the effects of heat transfer and fluid flow on the thernal, hydrologic, and mechanical response of a fault surface during seismic failure. Numerical modeling techniques are used to account for the coupling of the thermal, fluid-pressure, and stress fields. Results indicate that during an earthquake the failure surface is heated to a tempeature required for the thermal expansion of pore fluids to balance the rate of fluid loss due to flow and the fluid-volume changes due to pore dilatation. Once this condition is established, the pore fluids pressurize and the shear strength decreases rapidly to a value sufficient to maintain the thermal pressurization of pore fluids at near-lithostatic values. If the initial fluid pressure is hydrostatic, the final temperature attained on the failure surface will increase with depth, because a greater pressure increase can occur before a near-lithostatic pressure is reached. The rate at which thermal pressurization proceeds depends primarily on the hydraulic characteristics of the surrounding porous medium, the coefficient of friction on the fault surface, and the slip velocity. If either the permeability exceeds 10^–15 m² or the porous medium compressibility exceeds 10^–8 Pa^–1, then frictional melting may occur on the fault surface before thermal pressurization becomes significant. If the coefficient of friction is less than 10^–1 and if the slip velocity is less than 10^–2 msec^–1, then it is doubtful that either thermal pressurization or frictional melting on the fault surface could cause a reduction in the dynamic shear strength of a fault during an earthquake event.     

地幔物质上涌引起断层强度的变化和分布

本文根据多孔介质弹性变形与流体扩散在力学上的相互作用,选取了耦合本构关系式作为数学模型并对其求解,研究了地壳中三种类型断层在地幔物质上涌作用下,所产生的附加断层强度的变化和分布。主要结果是:(1)在空间上,“上涌”面所造成的断层强度软化区位于隆起的边翼地区;(2)在时间上,“上涌”对构造地震的触发滞后于隆起过程。本文通过对唐山地震的分析认为,地幔物质上涌在地震孕育和发生过程中起着重要的作用,并估算了震源面上的强度软化值。最后,求出了在“上涌”作用下,地壳中所产生的最大剪切应力分布。结果表明:“上涌”造成的破裂区域最易发生在地壳中层的隆起边翼地带,而且此区域形成的断层很可能是高倾角的断层。     

A Simple Relation to Constrain Groundwater Models Using Surface Deformation

A simple relation between pore pressure change and one-dimensional surface deformation is presented. The relation is for pore pressure change in a confined aquifer that causes surface deformation. It can be applied to groundwater models of any discretization and is computationally efficient. The estimated surface deformation from model results can be compared to observed surface deformation through geodetic techniques such as Differential Interferometric Synthetic Aperture Radar. Model parameters then are constrained using the observed surface deformation. The validity of this relation is shown through constraint of model parameters for surface uplift due to pore pressure increase caused by wastewater disposal injection.     

Inversion of near-field waveform data for earthquake source rupture process ( II ): Inversion of Chi Chi earthquake source, 20 September, 1999, Taiwan

The new inversion algorithm developed based on the recent progress in the nonlinear programming study by us is used to invert the earthquake source process of Chi Chi earthquake M _w7.6, 20 Semptember, 1999, Taiwan. A curve fault model is constructed in our inversion to make the fault model close to the real rupturing fault to reduce the influence from the discrepancy between the constructed fault model and the real rupturing fault. The results show that (1) the rupture process of the Chi Chi earthquake source lasted about 32 seconds and the main faulting occurred between 6th to 21st second after the start of the ruptures and the high slip area were mainly located at the northern segment of the fault. (2) The slip was dominated by thrust faulting. The average rake angle was 64.5°, which was very consistent with those inverted by USGS, Harvard and CWB (Central Weather Bureau of Taiwan). The amount of the moment inverted in this paper was 7.76×10²⁰ NM, which was a slightly bigger than those inverted by USGS and Harvard. (3) A clear nucleation step existed in the source faulting process and it lasted about 6 seconds. The moment release rate accelerated obviously at the end of the nucleation step. The faulting started from the southern segment and mainly occurred at the northern segment after 10 seconds. At the end of this paper, we analyzed the reliability of the inversion result via comparing with the GPS observations and discussed its scientific signification.     

Study of the Reservoir Damage during Cyclic Injection and Extraction Process in the Hutubi Gas Storage

During the injection and extraction process in gas storage,cyclic loading stress is applied to the reservoir pore structure,affecting the capacity and service life of gas storage. To study the influence of cyclic injection and extraction on reservoir physical property,a finite-element model of reservoir elastoplastic damage under the effect of pore pressure is established based on elastoplastic mechanics and fatigue-damage mechanics. Considering the mechanical property of the Hutubi reservoir and the injection-extraction mode of gas storage,the finite-element model is used to study the pore pressure and porosity change,fault safety,and ground surface deformation. The results show that the change in pore pressure and porosity,the deformation of the Hutubi fault,as well as the uplift and sink of the grounds surface are controlled by the fatigue damage of the reservoir pore structure during the injection-extraction process. After a 6-year service life, the cumulative deformation of the reservoir pore reaches its limit,causing pore structure collapse and porosity decrease. As a result,storage capacity declines and the ground surface starts to sink. Through this study,the safety of the Hutubi gas storage is evaluated,and some suggestions for field operations are offered.     

Fractal analysis of normal faults in northwestern Aegean area, Greece

Normal faults within the Ptolemais coal field and large seismogenic faults in the northwestern Aegean remain fractal for displacement values larger than about 1m. The kinematic parameters on reverse drag profiles such as length of rollover, footwall uplift and wavelength of footwall uplift show that all three parameters have a power law relationship, expressed by a c exponent of about 1, with the maximum displacement which take place across the fault. Footwall uplift/hanging wall subsidence ratio is about 1/2.The displacement analysis help us to propose a growth model for larger seismogenic faults in the NW Aegean, as is the ‘Hepiros fault set’ and the ‘Aliakmon fault zone’. Faults within the ‘Aliakmon fault zone’ were independently developed, at the first stages of deformation, by tip line deformation and out-of plane bifurcation, whereas later, deformation continued by segment linkage. One of these faults the ‘Sarakina fault’ was reactivated during the 1995 earthquake to produce a 25 km long surface rupture. A long term slip rate of about 0.3 mm a⁻¹ has been estimated by taking into consideration that over the past 6 Ma a maximum displacement of 1700 m across this fault has taken place.     

Characteristics of earthquake ground motions and the H/V of microtremors in the southwestern part of Taiwan

The relationships between the spectral characteristics of earthquake ground motions and those of micro‐tremors are investigated using the observed data from a dense strong‐motion network consisting of 108 stations in the Yun‐Li, Chia‐Yi and Tai‐Nan areas in southwestern Taiwan. Many high‐quality recordings, including those of the 921 Chi‐Chi earthquake (M_w=7.6), the 1022 Chia‐Yi mainshock (M_L=6.4), the 1022 major aftershock (M_L=6.0), as well as some weak motion events are selected to evaluate site responses. Microtremor measurements are also performed at most ground motion stations. With many stations in the area located on an alluvium structure, however, it is difficult to find good reference stations on rock sites, which therefore necessitates the calculation of single‐station H/V ratios. The predominant frequencies obtained from H/V ratios are consistent with those from spectral ratios. The site characteristics between the strong and weak events are different, however. This implies that a nonlinear effect probably occurred with the strong‐motion events. The main peak in the H/V spectra of the microtremors is in good agreement with the first peak obtained from the spectra of earthquake ground motions. It is reasonable to claim that the main peak reflects the deep underground structure. On the basis of the H/V ratios of the microtremors, it is concluded that the lower predominant frequencies appear in the plain area, while the higher values are near the mountainous region. Copyright © 2002 John Wiley & Sons, Ltd.     

Propagation of slip along frictional surfaces

The propagation of slip along a pre-existing frictional plane is formulated for the faults containing interstitial fluid. When normal and shear stresses satisfy the effective frictional law, a frictional sliding occurs stably or unstably, depending on the inhomogeneities of the surface. For a slip to rapidly sweep the surface, the pore pressure of fluid must exceed the critical value that is related to the physical or geometrical irregularities of the surface. If this condition fails, stable sliding is expected, analogous to a seismic fault creep. This prediction makes the role of water in seismic faulting more clear.     

Study on wall-type gravel drains as liquefaction countermeasure for underground structures

One of the most dramatic causes of damage to engineering structures during earthquakes has been the development of soil liquefaction beneath and around the structures. In order to dissipate the excess pore water pressures near structures, gravel drains are usually employed. In this study, the use of recycled concrete crushed stones as gravel drain materials is addressed. In order to investigate the performance of wall-type gravel drains, two series of shaking table tests were performed. The test results showed that gravel drains, when appropriate grain size distribution is considered, effectively dissipate the excess pore water pressure underneath the structure, and consequently reduce the magnitude of uplift. To supplement the laboratory tests, finite element analyses were also performed. For specified structure, ground and earthquake conditions, there is a critical width of gravel drain at which no uplift of structure will occur. The results of the model tests and the finite element analyses were then employed in developing design charts for determining the critical width of gravel drain to prevent buoyant rise of structure when the surrounding soil mass liquefies.     

Uncertainty in applying the linear poroelasticity model to field situations as a result of periodic loading in heterogeneous aquifers

The solid Earth's surface frequently experience changes in total stresses as a result of periodic loading. When the fluid‐saturated porous media deform in response to changes in stress, the induced variations in pore volume affect the pore water pressure. The fluid flow therefore occurs in response to the gradient in the induced excess pore water pressure. This work aims at quantifying the spatial variability in excess pressure head produced by the periodic loading accounting for the variation of log hydraulic conductivity (lnK). It is important for the rational management of groundwater resources. A closed‐form expression is developed by the nonstationary spectral approach to analyse the influence of the statistical properties of lnK process, the hydraulic parameters, and the spatial position. The general stochastic framework outlined in this work provides a basis for assessing the impact of statistical properties of input aquifer parameters on the output variability (or uncertainty). Copyright © 2014 John Wiley & Sons, Ltd.     

太湖草源性"湖泛"水域沉积物营养盐释放估算

于太湖草源性"湖泛"暴发期,采集柱状沉积物并应用peeper被动采样装置获得"湖泛"区原位沉积物间隙水.泥水样品分析表明:"湖泛"发生水域表层(0~7 cm)沉积物的含水率、孔隙度和有机质含量均明显高于对照区,其中有机质含量更是对照区样品的4倍左右,沉水植物残体促使表层沉积物物化性质改变的作用明显;"湖泛"发生水域表层沉积物间隙水中铵态氮(NH+4-N)、溶解性反应磷(SRP)及Fe2+含量远高于未发生区,植物残体降解对沉积物厌氧环境的营造显著.运用分子扩散模型对沉积物释放通量估算:"湖泛"发生区沉积物NH+4-N、SRP和Fe2+的释放速率分别是对照区的49.8、15.3和123.1倍.研究认为,草源性"湖泛"水体氮、磷等营养物含量升高的主要原因是沉积物的释放,而"湖泛"所营造的厌氧环境是氮、磷释放急剧增加的主要驱动因素.     

Poroelasticity: Finite element modelling of anomalous tilt and pore pressure caused by pumping in a sedimentary half space with fault

Extraction of groundwater or hydrocarbons causes pore pressure gradients and soil deformation due to poroelastic coupling. Recent studies show that high-resolution engineering tiltmeters installed at shallow depth between 2 and 10 m resolve this deformation. Models using poroelasticity can describe the relationship between fluid extraction, pore pressure gradients and induced tilt for homogeneous and layered sedimentary half spaces. Faults intersecting a stack of sedimentary layers, for example in the Lower-Rhine-Embayment, are of fundamental impact to the groundwater flow system of an area. However, the fault’s hydromechanical effect on pump induced tilt and the pore pressure regime is still poorly investigated. We chose a comparatively simple approach to quantify anomalous pump induced tilt and pore pressure observed near a fault and close to the surface in a sedimentary subsoil. A PC-based Finite Element software is used to model poroelastic deformation, i.e. modelling vertical tilt and excess pore pressure in response to fluid extraction through a singular well. We compare numerical solutions for models with and without faults and show that a fault can modify symmetry and amplitude of the deformation field by more than a magnitude. We conclude that tilt and pore pressure measurements also at shallow depth can thus be biased by large subsurface structures like faults. Vice versa, these measurements may provide means to quantify hydromechanical effects caused by subsurface structures. However, depending on the geological setting, i.e. if pathways are established by a fault, the anomaly caused by the fault can also be small and hard to detect. Therefore, faults and geological structures like material boundaries have to be considered in poroelastic models carefully. For tilt surveys with a limited number of instruments in geologically well constrained areas these models allow the preselection of potential positions for tiltmeters where prominent field anomalies are expected.     

Geological and mineralogical study of eclogite and glaucophane schists in the Naga Hills Ophiolite,Northeast India

A variety of low‐ to high‐pressure metamorphic assemblages occur in the metabasic rocks and metachert in the Upper Cretaceous–Eocene ophiolite belt of the central part of the Naga Hills, an area in the northern sector of the Indo–Myanmar Ranges in the Indo–Eurasian collision zone. The ophiolite suite includes peridotite tectonite containing garnet lherzolite xenoliths, layered ultramafic–mafic cumulates, metabasic rocks, basaltic lava, volcaniclastics, plagiogranite, and pelagic sediments emplaced as dismembered and imbricated bodies at thrust contacts between moderately metamorphosed accretionary rocks/basement (Nimi Formation/Naga Metamorphics) and marine sediments (Disang Flysch). It is overlain by coarse clastic Paleogene sediments of ophiolite‐derived rocks (Jopi/Phokphur Formation). The metabasic rocks, including high‐grade barroisite/glaucophane‐bearing epidote eclogite and glaucophane schist, and low‐grade greenschist and prehnite–clinochlore schist, are associated with lava flows and ultramafic cumulates at the western thrust contact. Chemically, the metabasites show a low‐K tholeiitic affinity that favors derivation from a depleted mantle source as in the case of mid‐ocean ridge basalt. Thermobarometry indicates peak P–T conditions of about 20 kb and 525°C. Retrogression related to uplift is marked by replacement of barroisite and omphacite by glaucophane followed by secondary actinolite, albite, and chlorite formation. A metabasic lens with an eclogite core surrounded by successive layers of glaucophane schist and greenschist provides field evidence of retrogression and uplift. Presence of S‐C mylonite in garnet lherzolite and ‘mica fish’ in glaucophane schist indicates ductile deformation in the shear zone along which the ophiolite was emplaced.     

Formation,uplift and dissection of planation surfaces at passive continental margins – a new approach

The usefulness of large‐scale, low‐relief, high‐level landscapes as markers of uplift events has become a subject of disagreement among geomorphologists. We argue that the formation of low‐relief surfaces over areas of large extent and cutting across bedrock of different age and resistance must have been guided by distinct base levels. In the absence of other options the most likely base level is sea level. We have analysed West Greenland landscapes in a recent study by combining the cooling history from apatite fission‐track analysis (AFTA) data with the denudation history from landscape analysis and the stratigraphic record. An important difference between our approach and that of classical geomorphology is that we now have the ability to document when thick sections of rocks have been deposited and then removed. The present‐day high‐level plateau in West Greenland is the remnant of a planation surface that was formed by denudation that lasted c. 20 million years during which up to 1 km of cover was removed after maximum burial at the Eocene–Oligocene transition. Here we present additional AFTA data to show that the planation surface is the end‐product of Cenozoic denudation even in basement areas and argue that Phanerozoic sediments – most likely of Cretaceous–Palaeogene age – must have been present prior to denudation. The planation surface was offset by reactivated faults and uplifted to present‐day altitudes of up to 2 km. The uplift occurred in two late Neogene phases that caused incision of valleys below the planation surface and their subsequent uplift. We therefore find that the elevated and deeply dissected plateau is evidence of episodic post‐rift uplift that took place millions of years after cessation of sea‐floor spreading west of Greenland. We suggest that other margins with similar morphology may also be characterized by episodic post‐rift uplift unrelated to the processes of rifting and continental separation, rather than being permanently uplifted since the time of rifting, as is commonly assumed. Copyright © 2009 John Wiley & Sons, Ltd.     

Constraints on the state of in situ effective stress and the mechanical behavior of ODP Leg 186 claystones in the Japan Trench forearc

Abstract The Japan Trench forearc offshore Honshu Island in northeast Japan, where the 130‐m.y.‐old Pacific oceanic plate is presently subducted, was drilled during the Ocean Drilling Program Leg 186. Results from mechanical and sedimentological studies of claystones recovered from Sites 1150 and 1151 in the overlying erosional forearc wedge are reported in the present study. Although many physical properties are similar in the seismic (Site 1150) and aseismic portion (Site 1151) of the shallow forearc, Site 1150 displayed a higher abundance of open fractures, two prominent fault zones and enigmatic pore fluid signatures in the claystones. The abundance of weak mineral phases, together with high smectite contents (from X‐ray diffraction), control the low friction coefficients of 0.33–0.39 of the claystones in ring‐shear experiments. Results from triaxial testing proposed overall low magnitudes of in situ effective vertical stress, with somewhat lower values at Site 1150 than at Site 1151. Similarly, samples from Site 1150 displayed slightly higher pore fluid pressures than those at Site 1151. The high sediment porosities, which are in part also a result of intact diatom tests (from scanning electron microscope), together with the anomalous fluid signatures and elevated pore fluid pressures, could very likely result from upward migration and influx of deep‐seated waters. Dewatering reactions at depth result in enhanced pore fluid pressure transients along out‐of‐sequence thrusts and consequently lower effective stress. At depths greater than that of Leg 186 drilling, elevated pressure–temperature conditions trigger mineral transformation and cementation, which result in increasing friction, unstable sliding and seismic rupture. Such earthquakes could have repeatedly disaggregated the consolidated claystone fabrics at the seismic site, and could be responsible for differences in yield strength and cementation when compared to the aseismic Site 1151.     

Influence of the hyporheic zone on the phosphorus dynamics of a large gravel‐bed river,Garonne River,France

Phosphorus (P) concentrations in sediments and in surface and interstitial water from three gravel bars in a large river (Garonne River, southern France) were measured daily, downstream of a wastewater treatment plant for a city of 740 000 inhabitants (Toulouse). Measurements were made of vertical hydraulic gradient (VHG), total dissolved phosphorus (TDP), soluble reactive phosphorus (SRP) and total phosphorus (TP) in water and of three extractable forms of phosphorus (water extractable, NaOH extractable and H₂SO₄ extractable) in hyporheic sediments from the gravel bars. Dissolved phosphorus was the major contributor to TP (74–79%) in both interstitial and surface waters on all sampling dates, and in most cases surface water P concentrations were significantly higher than interstitial concentrations. Hyporheic sediment TP concentrations ranged between 269 and 465 µg g^?1 and were highest in fine sediment fractions. Acid‐extractable P, a non‐bioavailable form, represented at least 95% of sediment TP. A positive relationship was observed between VHG and TP in two of the gravel bars, with wells that were strongly downwelling having lower TP concentrations. These results suggest that in downwelling zones, hyporheic sediments can trap surface‐derived dissolved P, and that much of this P becomes stored in refractory particulate forms. Bioavailable P is mainly present in dissolved form and only occupies a small fraction of total P, with particulate P comprising the majority of total P. Copyright © 2009 John Wiley & Sons, Ltd.     

Test of Hack's slope to bed material relationship in the Southern Eifel Uplands,Germany

The local reach gradient of small gravel bed rivers (drainage area 0-8-110 km²) in the Eifel, West Germany, is adjusted to transport the river bed sediments. Transport of gravel becomes possible under high flow conditions (Shields entrainment factor ≈-03). Mean bed material size for riffle sections increases with distance downstream. For small drainage areas channel slope is a negative exponential function of drainage area, while for the larger region the additional influence of bedload size has to be considered. Good agreement with Hack's data (1957) for Virginia and Maryland, U.S.A., is achieved (S = 0.0066 (D₅₀/A)^{- 40.}, r = 0.67).