The Role of Fault-Zone Architectural Elements on Pore Pressure Propagation and Induced Seismicity

We used hydrogeologic models to assess how fault-zone properties promote or inhibit the downward propagation of fluid overpressures from a basal reservoir injection well (150 m from fault zone, Q = 5000 m³/day) into the underlying crystalline basement rocks. We varied the permeability of the fault-zone architectural components and a crystalline basement weathered layer as part of a numerical sensitivity study. Realistic conduit-barrier style fault zones effectively transmit elevated pore pressures associated with 4 years of continuous injection to depths of approximately 2.5 km within the crystalline basement while compartmentalizing fluid flow within the injection reservoir. The presence of a laterally continuous, relatively low-permeability altered/weathered basement horizon (k_{altered layer} = 0.1 × k_basement) can limit the penetration depth of the pressure front to approximately 500 m. On the other hand, the presence of a discontinuous altered/weathered horizon that partially confines the injection reservoir without blocking the fault fluid conduit promotes downward propagation of pressures. Permeability enhancement via hydromechanical failure was found to increase the depth of early-time pressure front migration by a factor of 1.3 to 1.85. Dynamic permeability models may help explain seismicity at depths of greater than 10 km such as is observed within the Permian Basin, NM.     

Characterisation and multifaceted anisotropy assessment of Corvio sandstone for geological CO2 storage studies

We present a comprehensive characterisation of the physical, mineralogical, geomechanical, geophysical, and hydrodynamic properties of Corvio sandstone. This information, together with a detailed assessment of anisotropy, is needed to establish Corvio sandstone as a useful laboratory rock‐testing standard for well‐constrained studies of thermo–hydro–mechanical–chemical coupled phenomena associated with CO₂ storage practices and for geological reservoir studies in general. More than 200 core plugs of Corvio sandstone (38.1 and 50 mm diameters, 2:1 length‐to‐diameter ratio) were used in this characterisation study, with a rock porosity of 21.7 ± 1.2%, dry density 2036 ± 32 kg m^?3, and unconfined compressive and tensile strengths of 41 ± 3.28 and 2.3 ± 0.14 MPa, respectively. Geomechanical tests show that the rock behaves elastically between ～10 and ～18 MPa under unconfined conditions with associated Young's modulus and Poisson's ratio of 11.8 ± 2.8 GPa and 0.34 ± 0.01 GPa, respectively. Permeability abruptly decreases with confining pressure up to ～10 MPa and then stabilises at ～1 mD. Ultrasonic P‐ and S‐wave velocities vary from about 2.8–3.8 km s^?1 and 1.5–2.4 km s^?1, respectively, over confining and differential pressures between 0.1 and 35 MPa, allowing derivation of associated dynamic elastic moduli. Anisotropy was investigated using oriented core plugs for electrical resistivity, elastic wave velocity and attenuation, permeability, and tracer injection tests. Corvio sandstone shows weak transverse isotropy (symmetry axis normal to bedding) of <10% for velocity and <20% for attenuation.     

Searching for Regional Crustal Velocity-Density Relations with the Use of 2-D Gravity Modelling – Central Europe Case

In this paper we search for a reference relation between seismic P-wave velocity V and density ρ _ref for the continental crust. Based on the results of modern seismic experiments, we compiled 2-D seismic models into a network of four, each about 1100–1400 km long, continental-scale seismic transects cutting all main tectonic units in Central Europe. The Moho depth (about 52 km beneath the TESZ in SE Poland, to about 25 km beneath the Pannonian Basin) and the crustal structure of this area are characterised by a large variation. This structural variation provides an interesting basis for gravity studies and especially for analysing the difference of the density structure between two major tectonic provinces of distinctive age difference: Precambrian and Phanerozoic. The 2-D gravity modelling applied for crustal cross-sections representing the regional structure, based on a unified gravity anomaly map of the area, allows for a stable determination of some general features of the regional reference velocity-density relation for the continental crust. In general three major seismo-petrological types of rocks can be distinguished: sediments, crystalline crust and mantle. In compacted sediments the reference velocity-density relation is well described by the Gardner or Nafe-Drake model. Calculated gravity anomalies, using unified velocity-density relation for the whole crystalline crust, well describe observed anomalies, with an average difference of 14 mGal. However, calculated gravity anomalies, using separated velocity-density relations for the crystalline crust of Precambrian and Phanerozoic Europe, describe observed anomalies better than for the entire crust, with an average difference 12 mGal. The most important feature of these relations is the large differentiation of the derivative dρ _ref /dV in the crystalline crust, being about 0.3 g s/m⁴ for Precambrian, and about 0.1 g s/m⁴ for the Phanerozoic crystalline crust. The modelling suggests a very small density value in the uppermost mantle ρ = 3.11 g/cm³ below the younger area, while for the older area it is ρ = 3.3 g/cm³.     

Toppling failures from alpine cliffs on ben Lomond,Tasmania

Large cliff failures involving forward toppling over a basal hinge have occurred on more than half of the plateau edge of Ben Lomond, northeastern Tasmania. This mode of failure, which is readily identified from the columnar structure of the dolerite involved, has affected up to 10⁷ m³ of rock at a time and a total of more than 50 × 10⁶ m³ in all the cases which can still be identified. It represents perhaps the most important form of cliff retreat, amounting to a rate of 0.2 mm yr^?1 over the last 100,000 years. Topographic evidence and joint surveys suggest that two different mechanisms have produced the topples on Ben Lomond. One has involved failure in the sediments underlying the dolerite with consequent foundering and cambering of large sections of the plateau edge. This mechanism accounts for relatively few of the Ben Lomond topples, though it includes the largest individual cases. The second mechanism, dominant in most of the topples, involved slab failure in the cliffs. Both modes of failure have been facilitated by vertical weaknesses within the bedrock and both require an initially steep cliff profile. Because of the latter requirement, which is not met on the other mountains of northeastern Tasmania, large-scale topples are found only on Ben Lomond, and only there where glacial steepening of the cliff has been possible. Following the initial failure, topples of both types have migrated downslope by block sliding for distances up to 2 km.     

The effect of screening of the magnetotelluric field by high resistivity layers: Three-dimensional numerical modeling constraints

A 3-D numerical modeling approach is applied to the effect of screening of the magnetotelluric (MT) field by high resistance (≥10⁶ Ω m²) layers, which up to now has mainly been studied in terms of 2-D heterogeneous models of geological sections [Berdichevsky and Yakovlev, 1990]. Three-layer models above a poorly conducting basement (ρ=10⁴ Ω m) are used. The resistance of the screening layer was generally taken equal to 2×10⁶ Ω m² and its thickness was varied from 200 to 2000 m. The resistivity of the host medium was set at 10 Ω m. Heterogeneities of a low and a high resistivity (ρ=10^?3 and 10⁴ Ω m, respectively) ranging in horizontal size from 4 to 40 km and having a height of 1–2.8 km (a protrusion on the basement) were examined. Based on calculations of these models, 2-D and 3-D screening effects were compared. The 3-D modeling determined values of model parameters at which the screening properties are preserved both in the case of a screen of limited horizontal dimensions and in the presence of a rupture in the screen comparable in horizontal size with the 3-D heterogeneity. As follows from the modeling results, the screening effect of the high resistivity model layers seriously complicates the use of MT soundings for the identification of a local heterogeneity of both a low and a high resistivity (ρ=10^?3 and 10⁴ Ω m, respectively) if its lateral size is smaller than the 4-km thickness of the studied three-layer section in question and its height is 1–1.5 km. The regular patterns of the screening effect revealed in this work are of interest in electromagnetic sounding applications.     

THE INDUCTIVE RESPONSE OF A HORIZONTAL CONDUCTING CYLINDER BURIED IN A UNIFORM EARTH FOR A UNIFORM INDUCING FIELD*

The inductive response of a conducting horizontal cylinder embedded in a uniform earth is studied using numerical results obtained for an analytical solution for the problem of a conducting cylinder buried in a homogeneous earth for the case of a uniform inducing field. A check of the validity of the numerical results is made by a comparison with analogue model measurements for a number of cases. Numerical results for a range of cylinder radii (a = 1–10 km), depths of burial (d= 0–4 km), conductivity contrasts (σ₂= 10^?2-10 Sm^?1), and source frequencies (f= 10^?1-10^?4 Hz) of interest in the interpretation of magnetotelluric field measurements are presented. The results indicate that for a uniform inducing field the conductivity and depth of burial of a horizontal cylindrical inhomogeneity are best determined through a measurement of the amplitudes H_y, H_z and E_x and the phases φ_y and Ψ_x.     

Greenhouse gas emissions (CO2, CH4, and N2O) from several perialpine and alpine hydropower reservoirs by diffusion and loss in turbines

We investigated greenhouse gas emissions (CO₂, CH₄, and N₂O) from reservoirs located across an altitude gradient in Switzerland. These are the first results of greenhouse gas emissions from reservoirs at high elevations in the Alps. Depth profiles were taken in 11 reservoirs located at different altitudes between the years 2003 and 2006. Diffusive trace gas emissions were calculated using surface gas concentrations, wind speeds and transfer velocities. Additionally, methane entering with the inflowing water and methane loss at the turbine was assessed for a subset of the reservoirs. All reservoirs were emitters of carbon dioxide and methane with an average of 970?±?340?mg?m^?2?day^?1 (results only from four lowland and one subalpine reservoir) and 0.20?±?0.15?mg?m^?2?day^?1, respectively. One reservoir (Lake Wohlen) emitted methane at a much higher rate (1.8?±?0.9?mg?m^?2?day^?1) than the other investigated reservoirs. There was no significant difference in methane emissions across the altitude gradient, but average dissolved methane concentrations decreased with increasing elevation. Only lowland reservoirs were sources for N₂O (72?±?22???g?m^?2?day^?1), while the subalpine and alpine reservoirs were in equilibrium with atmospheric concentrations. These results indicate reservoirs from subalpine/alpine regions to be only minor contributors of greenhouse gases to the atmosphere compared to other reservoirs.     

Soil degradation in central Spain due to sheet water erosion by low‐intensity rainfall events

Runoff and sediment lost due to water erosion were recorded for 36 (1 m²) plots with varying types of vegetative cover located on sloping gypsiferous fields in the South of Madrid. 75% of the events had maximum 30‐minute intensity (I₃₀) less than 10 mm h^?1 in the period studied (1994–2005). As for the vegetative cover, maximum correlation between runoff and soil loss was found in the least protected plots (0–40% cover) during the most intense rainfall events; however, a significant positive correlation was also observed in plots with greater coverage (40–60%). If coverage exceeded 60%, rainfall erosivity declined. The average amount of sediment produced in high‐intensity events was significantly greater (approximately 7 g m^?2 per I₃₀ event >10 mm h^?1) than that produced in the rest of the moderate‐intensity events (approximately 3 g m^?2 per I₃₀ event <10 mm h^?1), but due to the high rate of occurrence of the latter throughout the year sediment loss during the period studied totaled 128 g m^?2. By comparison, only 40 g m^?2 was produced by the I₃₀ events greater than 10 mm h^?1. Even though the amount of soil lost is relatively insignificant from a quantitative standpoint, the organic matter content lost in the sediment (six times more than in the soil) is a permanent loss that threatens the development of the surface of the soil in this area when the vegetative cover is less than 40%. The soil here experiences a chronic loss of 0·02 mm annually as a consequence of frequent, moderate events, in addition to any loss produced by extraordinary events, which, though less frequent, are much more erosive. If moderate events are ignored, an important part of soil loss will be lost in the long run. Copyright © 2007 John Wiley & Sons, Ltd.     

Energy balance and evaporation loss of an agricultural reservoir in a semi‐arid climate (south‐eastern Spain)

A typical agricultural water reservoir (AWR) of 2400 m² area and 5 m depth, located in a semi‐arid area (southern Spain), was surveyed on a daily basis for 1 year. The annual evaporation flux was 102·7 W m^?2, equivalent to an evaporated water depth of 1310 mm year^?1. The heat storage rate G exhibited a clear annual cycle with a peak gain in April (G ～ 45 W m^?2) and a peak loss in November (G ～ 40 W m^?2), leading to a marked annual hysteretic trend when evaporation (λE) was related to net radiation (R_n). λE was strongly correlated with the available energy A, representing 91% of the annual AWR energy loss. The sensible heat flux H accounted for the remaining 9%, leading to an annual Bowen ratio in the order of 0·10. The equilibrium and advective evaporation terms of the Penman formula represented 76 and 24%, respectively, of the total evaporation, corresponding to a annual value of the Priestley–Taylor (P–T) coefficient (α) of 1·32. The P–T coefficient presented a clear seasonal pattern, with a minimum of 1·23 (July) and a maximum of 1·65 (December), indicating that, during periods of limited available energy, AWR evaporation increased above the potential evaporation as a result of the advection process. Overall, the results stressed that accurate prediction of monthly evaporation by means of the P–T formula requires accounting for both the annual cycle of storage and the advective component. Some alternative approaches to estimating R_n, G and α are proposed and discussed. Copyright © 2009 John Wiley & Sons, Ltd.     

Computation of Solar Radiative Fluxes by 1D and 3D Methods Using Cloudy Atmospheres Inferred from A-train Satellite Data

This study used realistic representations of cloudy atmospheres to assess errors in solar flux estimates associated with 1D radiative transfer models. A scene construction algorithm, developed for the EarthCARE mission, was applied to CloudSat, CALIPSO and MODIS satellite data thus producing 3D cloudy atmospheres measuring 61 km wide by 14,000 km long at 1 km grid-spacing. Broadband solar fluxes and radiances were then computed by a Monte Carlo photon transfer model run in both full 3D and 1D independent column approximation modes. Results were averaged into 1,303 (50 km)² domains. For domains with total cloud fractions A _c  < 0.7 top-of-atmosphere (TOA) albedos tend to be largest for 3D transfer with differences increasing with solar zenith angle. Differences are largest for A _c  > 0.7 and characterized by small bias yet large random errors. Regardless of A _c , differences between 3D and 1D transfer rarely exceed ±30 W m^?2 for net TOA and surface fluxes and ±10 W m^?2 for atmospheric absorption. Horizontal fluxes through domain sides depend on A _c with ～20% of cases exceeding ±30 W m^?2; the largest values occur for A _c  > 0.7. Conversely, heating rate differences rarely exceed ±20%. As a cursory test of TOA radiative closure, fluxes produced by the 3D model were averaged up to (20 km)² and compared to values measured by CERES. While relatively little attention was paid to optical properties of ice crystals and surfaces, and aerosols were neglected entirely, ～30% of the differences between 3D model estimates and measurements fall within ±10 W m^?2; this is the target agreement set for EarthCARE. This, coupled with the aforementioned comparison between 3D and 1D transfer, leads to the recommendation that EarthCARE employ a 3D transport model when attempting TOA radiative closure.     

Submarine pyroclastic deposits formed during the 20th May 2006 dome collapse of the Soufrière Hills Volcano, Montserrat

The 20th May 2006 lava dome collapse of the Soufrière Hills Volcano, Montserrat, had a total non-dense rock equivalent (non-DRE) collapse volume of approximately 115?×?10⁶?m³. The majority of this volume was deposited into the ocean. The collapse was rapid, 85% of the mobilized volume being removed in just 35?min, giving peak pyroclastic flow flux of 66?×?10³?m³?s^?1. Channel and levee facies on the submarine flanks of the volcano and formation of a thick, steep-sided ridge, suggest that the largest and most dense blocks were transported proximally as a high concentration granular flow. Of the submerged volume, 30% was deposited from the base of this granular flow, forming a linear, high-relief ridge that extends 7?km from shore. The remaining 70% of the submerged volume comprises the finer grain sizes, which were transported at least 40?km by turbidity currents on gradients of <2°. At several localities, the May 2006 distal turbidity currents ran up 200?m of topography and eroded up to 20?cm of underlying substrate. Multiple turbidites are preserved, representing current reflection from the graben margins and deflection around topography. The high energy of the May 2006 collapse resulted in longer submarine run out than the larger (210?×?10⁶?m³) Soufrière Hills dome collapse in July 2003.     

Spatial and temporal variations of pCO2, dissolved inorganic carbon and stable isotopes along a temperate karstic watercourse

This study investigated CO₂ degassing and related carbon isotope fractionation effects in the Wiesent River that drains a catchment in the karst terrain of the Franconian Alb, Southern Germany. The river was investigated by physico‐chemical and stable isotope analyses of water and dissolved inorganic carbon during all seasons along 65‐km long downstream transects between source and mouth. Calculated pCO₂ values at the source were 21 400 ± 2400 µatm. The pCO₂ rapidly decreased in the river water and dropped to an average of 1240 ± 330 µatm near the mouth. About 90% of this decrease occurred within the first 6 km of the river. The river was supersaturated with respect to CO₂ over its entire course and must have acted as a continuous year‐round CO₂ source to the atmosphere. The average CO₂ flux from the karst river was estimated with 450 mmol m^?2 day^?1 with higher fluxes up to 5680 mmol m^?2 day^?1 at the source. At the source, δ¹³C_DIC values showed no seasonal variations with an average of ?14.2 ± 0.2‰. This indicated that groundwater retained high pCO₂ mainly from soil CO₂. The contribution of soil CO₂ to dissolved inorganic carbon was estimated at 65% to 72%. The downstream CO₂ loss caused a positive shift in δ¹³C_DIC values of 2‰ between source and mouth because of the preferential loss of the ¹²C isotope during degassing. Considering the findings of this study and the fact that carbonate lithology covers a significant part of the earth's surface, CO₂ evasion from karst regions might contribute notably to the annual carbon dioxide release from global freshwater systems. Copyright © 2015 John Wiley & Sons, Ltd.     

Zur ökonomischen Bedeutung der biologischen Selbstreinigung

In twelve flowing water (Q 1.04… 30.4 m³/s, 0.5… 1.2 m/s) the self-purification efficiency in 110 river sections (L 1.3… 26.4 km) is determined as the load difference in kg · d^?1 COD-Mn as well as in g · m^?3 referred to the daily passage. In the economical comparison, the self-purification efficiency is valued as substitution for wastewater treatment plants of the same capacity. At a mean specific capacity of 8.6 g · m^?3, the self-purification efficiency of the investigated waters is equivalent to economical values of 124 · 10³… 534 · 10³ M km^?1 investment costs, 11 · 10³… 80 · 10³ M · a^?1km^?1 operating costs and 5… 81 MWh · a^?1 km^?1 expenditure of energy. The specific capacity in g · m^?3 COD-Mn shows an exponential regression to the degree of saprobity (L = 0.015 · exp (1.7358 · S)). From this empirical model the limits of the self-purification capacity of aerobic waters by oxygen input can be detected: in respective examples more than 50% of the required oxygen input are due to weirs.     

On the use of the groundwater replenishment periods approach to determining aquifer parameters

Chalk crops out from a wide belt around the Paris basin, France, covering an area of about 70 000 km². In this region, the chalk presents the most important unconfined aquifer because of its extent and the size of its resources ((11–12) × 10⁹ m³ year^?1). The assessment of underground outflow depends on the vertical feeding, the infiltration and the hydrometry. This paper analyses the regional structural map, interprets the groundwater reaction under rainfall, explains the water circulation in such media where the reservoir geometry plays an essential role on aquifer response, and determines different aquifer physical parameters. Attempts are made to identify the rapid transfer of groundwater at the level of faults and their important fissures. Based on a quantitative study of the seasonal and interannual piezometry fluctuations, it is noted that the mode and the piezometric chronology events are controlled by major parameters of geological and hydrogeologic contexts, aquifer hydraulic characteristics, the position of upstream and downstream basin limits, groundwater depth and replenishment time. This paper ends with determination of groundwater physical parameters of the diffusivity (T/S, T and S) values by examination of the groundwater replenishments periods. Copyright © 2005 John Wiley & Sons, Ltd.     

Zur Rolle der Stickstoffixation im Stoffhaushalt eines geschichteten Sees

In the dimict lake Arend (5.1 km², 146 hm³, 49.5 m z_max), nitrogen is production-limiting with concentrations below the detection limit during the production period. Phytoplankton achieves biomasses of up to 18 mg/l fresh matter, essential contributions being made by Aphanizomenon with 2 mg/l and Anabaena with up to 10 mg/l. Nitrogen fixation was measured by the ethine reduction technique (acetylene reduction) during periods of the occurrence of heterocystforming Cyanophyceae and achieved peak values up to 6.59 μg N² · h^?1 · l^?1 or 14.87 m^?2 · h^?1 g N₂ · m^?2 · h^?1. The rates of fixation show a safe correlation with the biomass of heterocyst-containing Cyanophyceae (r = 0.88), their development beginning at values below the N : P-ratio of 2.66.     

Geochemical characteristics of pyrolysis gas from epimetamorphic rocks in the northern basement of Songliao Basin,Northeast China

The Xushen gas field, located in the north of Songliao Basin, is a potential giant gas area for China in the future. Its proved reserves have exceeded 1000×10⁸ m³ by the end of 2005. But, the origin of natural gases from the deep strata is still in debating. Epimetamorphic rocks as a potential gas source are widely spreading in the northern basement of Songliao Basin. According to pyrolysis experiments for these rocks in the semi-confined system, gas production and geochemistry of alkane gases are discussed in this paper. The Carboniferous-Permian epimetamorphic rocks were heated from 300°C to 550°C, with temperature interval of 50°C. The gas production was quantified and measured for chemical and carbon isotopic compositions. Results show that δ ¹³C₁ is less than −20‰, carbon isotope trend of alkane gas is δ ¹³C₁<δ ¹³C₂<δ ¹³C₃ or δ ¹³C₁<δ ¹³C₂>δ ¹³C₃, these features suggest that the gas would be coal-type gas at high-over maturity, not be inorganic gas with reversal trend of gaseous alkanes (δ ¹³C₁>δ ¹³C₂>δ ¹³C₃). These characteristics of carbon isotopes are similar with the natural gas from the basin basement, but disagree with gas from the Xingcheng reservoir. Thus, the mixing gases from the pyrolysis gas with coal-typed gases at high-over maturity or oil-typed gases do not cause the reversal trend of carbon isotopes. The gas generation intensity for epimetamorphic rocks is 3.0×10⁸–23.8×10⁸ m³/km², corresponding to R _o from 2.0% to 3.5% for organic matter.

     

Velocity Structure of the West-Bengal Sedimentary Basin, India along the Palashi-Kandi Profile Using a Travel-time Inversion of Wide-angle Seismic data and Gravity Modeling – An Update

The 2-D shallow velocity structure along the north-south Palashi-Kandi profile in the West Bengal sedimentary basin has been updated by travel-time inversion of seismic refraction, wide-angle reflection and gravity data. A six-layer shallow model up to a depth of about 7 km has been derived. The first layer, which has an average velocity of 2.0 kms^?1, represents the alluvium deposit, which rests over the shale formation with average velocity of 3.0 kms^?1. The thin (200 m) Sylhet limestone, observed at a nearby Palashi well, remains hidden in the present data set. Hence a 200-m thin layer with a velocity of 3.7 kms^?1, corresponding to the Sylhet limestone, has been assumed to be present throughout the profile. The fourth layer with a velocity of 4.5–4.7 kms^?1 at a depth of 1.7–2.4 km represents the Rajmahal traps. The ‘skip’ phenomenon and rapid amplitude decay of first arrivals indicate a low-velocity zone (LVZ) in the study area. Using the ‘skip’ phenomena and wide-angle reflection data, identified on seismograms, the LVZ with a velocity of 4.0 kms^?1, indicating the Gondwana sediments, has been delineated below the Rajmahal traps. The next layer with a velocity 5.4–5.6 kms^?1 overlying the crystalline basement (5.8–6.25 kms^?1) may be associated with the Singhbhum group of meta volcanic rock that has been exposed in the western part of the basin. The basement lies at a variable depth of 4.9 to 6.8 km. The overall uncertainties of various velocity and boundary nodes are ± 0.15 kms^?1 and ± 0.5 km, respectively. The elevated basement feature in the north might have acted as a structural barrier for the deposition of Sylhet limestone during the Eocene epoch. The seismically derived shallow structure correctly explains the observed Bouguer gravity anomaly along the profile. The addition of reflections in the present analysis provides a stronger control on the depths and velocities of basement and overlying sedimentary formations, compared to the earlier model derived mainly by the first arrival seismic data.     

Examination of a Rock Failure Criterion Based on Circumferential Tensile Strain

—Uniaxial compression, triaxial compression and Brazialian tests were conducted on several kinds of rock, with particular attention directed to the principal tensile strain. In this paper we aim to clarify the effects of the experimental environment—such as confining pressure, loading rate, water content and anisotropy—on the critical tensile strain, i.e., the measured principal tensile strain at peak load.¶It was determined that the chain-type extensometer is a most suitable method for measuring the critical tensile strain in uniaxial compression tests. It is also shown that the paper-based strain gage, whose effective length is less than or equal to a tenth of the specimen’s diameter and glued on with a rubber-type adhesive, can be effectively used in the Brazilian tests.¶The effect of confining pressure P _C on the critical tensile strain ? _TC in the brittle failure region was between ?0.02 × 10^?10 Pa^?1 and 0.77 × 10^?10 Pa^?1. This pressure sensitivity is small compared to the critical tensile strain values of around ?0.5 × 10^?2. The strain rate sensitivities ?? _TC /?{log(d|?|/dt)} were observed in the same way as the strength constants in other failure criteria. They were found to be from ?0.10 × 10^?3 to ?0.52 × 10^?3 per order of magnitude in strain rate in the triaxial tests. The average magnitude of the critical tensile strain ? _TC increased due to the presence of water by 4% to 20% for some rocks, and decreased by 22% for sandstone. It can at least be said that the critical tensile strain is less sensitive to water content than the uniaxial compressive strength under the experimental conditions reported here. An obvious anisotropy was observed in the P-wave velocity and in the uniaxial compressive strength of Pombetsu sandstone. It was not observed, however, in the critical tensile strain, although the data do show some variation.¶A "tensile strain criterion" was proposed, based on the above experimental results. This criterion signifies that stress begins to drop when the principal tensile strain reaches the critical tensile strain. The criterion is limited to use within the brittle failure region. The critical tensile strain contains an inelastic strain component as well as an elastic one. It is affected by the strain rate, however, it is relatively insensitive to the confining pressure, the presence of water and anisotropy.     

The reservoir-associated earthquakes of April 1983 in western Thailand: Source modeling and implications for induced seismicity

On 22 April 1983, a very large area of Thailand and part of Burma were strongly shaken by a rare earthquake (m _b=5.8,M _s=5.9). The epicenter was located at the Srinagarind reservoir about 190 km northwest of Bangkok, a relatively stable continental region that experienced little previous seismicity. The main shock was preceded by some foreshocks and followed by numerous aftershocks. The largest foreshock ofm _b=5.2 occurred 1 week before the main shock, and the largest aftershock ofm _b=5.3 took place about 3 hours after the main shock. Focal mechanisms of the three largest events in this earthquake sequence have been studied by other seismologists using firts-motion data. However, the solutions for the main shock and the largest aftershock showed significant inconsistency with known surface geology and regional tectonics. We reexamined the mechanisms of these three events by using teleseismicP-andS-waveforms and through careful readings ofP-wave first motions. The directions of theP axes in our study range from NNW-SSE to NNE-SSW, and nodal planes strike in the NW-SE to about E-W in agreement with regional tectonics and surface geology. The main shock mechanism strikes 255°, dips 48°, and slips 63.5°. The fault motions during the main shock and the foreshock are mainly thrust, while the largest aftershock has a large strike-slip component. The seismic moment and the stress drop of the mainshock are determined to be 3.86×10²⁴ dyne-cm and 180 bars, respectively. The occurrence of these thrust events appears to correlate with the unloading of the Srinagarind reservoir. The focal depths of the largest foreshock, the main shock, and the largest aftershock are determined to be 5.4 km, 8 km, and 22.7 km, respectively, from waveform modeling and relative location showing a downward migration of hypocenters of the three largest events during the earthquake sequence. Other characteristics of this reservoir-induced earthquake sequence are also discussed.