Wave velocity and stress/strain in rock brittle failure

This paper reports the relationship between wave velocity and stress/strain during rock brittle failure under compression. It is assumed that the rock is a strain-softening medium whose strength can be described by Weibull’s distribution. Based on renormalization group theory, it is found that the stress ratio λ _c (the ratio of the stress at the critical point to the peak stress) depends mainly on the homogeneity index or shape parameter m in the Weibull’s distribution for the rock. Measured stress–strain curves and corresponding wave velocity–strain curves in four case studies suggest that the changes caused by internal cracking are correlated to variations of rock wave velocity, and the stress at the critical point on the stress–strain curve is considered to correspond to the rapid decreasing point on the wave velocity–strain curve. For uniaxial compression, the critical stress ratio λ _c is derived to be about 75.4 %, which is very close to the measured value.     

An Experimental Study of the Fracture Coalescence Behaviour of Brittle Sandstone Specimens Containing Three Fissures

To analyse the fracture coalescence behaviour of rock, rectangular prismatic sandstone specimens (80?×?160?×?30?mm in size) containing three fissures were tested under uniaxial compression. The strength and deformation behaviours of the specimens are first analysed by investigating the effects of the ligament angle β₂ on the peak strength, peak strain and crack initiation stress of the specimens. To confirm the sequence of crack coalescence, a photographic monitoring technique is used throughout the entire period of deformation. Based on the results, the relationship between the real-time crack coalescence process and the axial stress–strain curve of brittle sandstone specimens is also developed, and this relationship can be used to evaluate the macroscopic deformation characteristics of pre-cracked rock. The equivalent strain evolution fields of the specimen, with α?=?β₁?=?45° and β₂?=?90°, are obtained using the digital image correlation technique and show good agreement with the experimental results of pre-cracked brittle sandstone. These experimental results are expected to improve the understanding of fracture mechanisms and be used in rock engineering with intermittent structures, such as deep underground excavated tunnels.     

Is there Link between the Type of the Volumetric Strain Curve and Elastic Constants, Porosity, Stress and Strain Characteristics ?

The stress [crack damage stress (σ _cd) and uniaxial compressive strength (σ _c)] and strain characteristics [maximum total volumetric strain (ε _cd), axial failure strain (ε _af)], porosity (n) and elastic constants [elastic modulus (E) and Poisson’s ratio (ν)] and their ratios were coordinated with the existence of two different types (type 1 and type 2) of volumetric strain curve. Type 1 volumetric strain curve has a reversal point and, therefore, σ _cd is less than the uniaxial compressive strength (σ _c). Type 2 has no reversal point, and the bulk volume of rock decreases until its failure occurs (i.e., σ _cd = σ _c). It is confirmed that the ratio between the elastic modulus (E) and the parameter λ = n/ε _cd strongly affects the crack damage stress (σ _cd) for both type 1 and type 2 volumetric strain curves. It is revealed that heterogeneous carbonate rock samples exhibit different types of the volumetric strain curve even within the same rock formation, and the range of σ _cd/σ _c = 0.54–1 for carbonate rocks is wider than the range (0.71 < σ _cd/σ _c < 0.84) obtained by other researchers for granites, sandstones and quartzite. It is established that there is no connection between the type of the volumetric strain curve and values of n, E, σ _cd, ν, E/(1 ? 2ν), M _R = E/σ _c and E/λ. On the other hand, the type of volumetric strain curve is connected with the values of λ and the ratio between the axial failure strain (ε _af) and the maximum total volumetric strain (ε _cd). It is argued that in case of small ε _af/ε _cd–small λ, volumetric strain curve follows the type 2.     

Experimental study on mechanical characteristics of sandstone containing arc fissures

Arc fissure is one of the basic forms of defective rock, where the expansion and evolution mechanism plays an important role in the stability of engineering rock mass under the external load action. Uniaxial compression experiments of sandstone samples that contained various angles of arc fissures (sandstone sample was 80 mm?×?160 mm?×?30 mm) were performed in order to investigate the effect that arc angle α had on the mechanical properties, the failure mode, and the fracture evolution process of sandstone. The results showed that when arc angle α was increased, the peak strength and the strain of the sandstone samples initially decreased before increasing and the minimum peak strength and strain were reached when α?=?15°. The deterioration of the bearing capacity and the number of cracks that appeared during the sandstone loading process decreased as the arc angle of the fissure increased. The arc fissure destruction was primarily initiated from the fragile area of the arch tip. The tensile cracks appeared on the fissure tip and non-tip as the axial force increased. The various arc angle α played an important role in the initiation stress and the rupture evolution of the specimen.     

Rainfall-induced landslide in the active frontal fold–thrust belt of Northwestern Himalaya,Jammu: dynamics inferred by geological evidences and Ground Penetrating Radar

Landslides are the most established geological hazards in the frontal fold–thrust belt of Northwestern Himalaya comprising of Siwaliks and Murree strata. The continuous rainfall from 2 to 6 September, 2014 caused a massive landslide at village Sadal in Udhampur district of Jammu and Kashmir state. The landslide occurred in the early morning of September 6, 2014, destroying entire Sadal habitation comprising 45 houses, and killing 41 people and more than 500 domestic animals. Google earth images of pre and post-landslide events along with the field measurements show the kinematics of upper and lower parts of the slide. Horizontal and vertical components of displacement and mode of failure suggest the landslide as of complex nature. The shallow subsurface geophysical imaging through Ground Penetrating Radar (GPR) survey shows the failure plane composed of friable mudstone bed underlain by massive mudstone and overlain by cross-bedded sandstone. The depth of debris material above the failure plane ranges from 6 m at Site S1a-b to 10 m at Site-S2b and 20 m at Site S3a. The velocity analysis of Site-3 shows four thick layers represented from bottom to surface by L1—sandstone (V?=?0.16 m/ns, travel time?=?356.36 ns), L2—mudstone (V?=?0.17 m/ns, travel time?=?288.48 ns), L3—massive mudstone (V?=?0.19 m/ns, travel time 220.68 ns), and L4—cross-laminated sandstone (V?=?0.20 m/ns, travel time?=?77.58 ns) overlaying the failure plane. The study shows the landslide occur along the western limb of a fold identified during the present work. We mapped an old landslide on the same limb which shows 5–6 m-thick subsurface debris material with thick rock fragments involved in the landslide process. The detailed geological and geophysical investigations suggest that both the landslides were triggered by extreme rain fall events.     

Characteristics of landslides in western Colorado,USA

Mass movement can be activated by earthquakes, rapid snowmelt, or intense rainstorms in conjunction with gravity. Whereas mass movement plays a major role in the evolution of a hillslope by modifying slope morphology and transporting material from the slope to the valley, it is also a potential natural hazard. Determining the relationships of frequency and magnitude of landslides are fundamental to understanding the role of landslides in the study of landscape evolution, hazard assessment, and determination of the rate of hillslope denudation. We mapped 735 shallow and active landslides in the Paonia to McClure Pass area of western Colorado from aerial photographs and field surveys. The study area covers ～815 km². The frequency–magnitude relationships of the landslides illustrate the flux of debris by mass movement in the area. The comparison of the probability density of the landslides with the double Pareto curve, defined by power scaling for negative slope (α), power scaling for positive slope (β), and location of rollover (t), shows that α?=?1.1, β?=?1.9, and t?=?1,600 m² for areas of landslides and α?=?1.15, β?=?1.8, and t?=?1,900 m³ for volumes of landslides. The total area of landslides is 4.8?×?10⁶ m² and the total volume of the landslides is 1.4?×?10⁷ m³. The areas (A) and the volumes (V) of landslides are related by V?=?0.0254?×?A ^1.45. The frequency–magnitude analysis shows that landslides with areas ranging in size from 1,600 to 20,000 m² are the most hazardous landslides in the study area. These landslides are the most frequent and also do a significant amount of geomorphic work. We also developed a conceptual model of hillslope development to upland plateau driven by river incision, shallow landsliding, and deep-seated large landsliding. The gentle slope to flat upland plateau that dominated the Quaternary landscape of the study area was modified to the present steep and rugged topography by the combined action of fluvial incision and glacial processes in response to rock uplift, very-frequent shallow landsliding, and less-frequent deep-seated landsliding.     

Influences of loading direction and intermediate principal stress ratio on the initiation of strain localization in cross-anisotropic sand

Since cross-anisotropic sand behaves differently when the loading direction or the stress state changes, the influences of the loading direction and the intermediate principal stress ratio (b = (σ ₂ ? σ ₃)/(σ ₁ ? σ ₃)) on the initiation of strain localization need study. According to the loading angle (angle between the major principal stress direction and the normal of bedding plane), a 3D non-coaxial non-associated elasto-plasticity hardening model was proposed by modifying Lode angle formulation of the Mohr–Coulomb yield function and the stress–dilatancy function. By using bifurcation analysis, the model was used to predict the initiation of strain localization under plane strain and true triaxial conditions. The predictions of the plane strain tests show that the major principal strain at the bifurcation points increases with the loading angle, while the stress ratio decreases with the loading angle. According to the loading angle and the intermediate principal stress ratio, the true triaxial tests were analyzed in three sectors. The stress–strain behavior and the volumetric strain in each sector can be well captured by the proposed model. Strain localization occurs in most b value conditions in all three sectors except for those which are close to triaxial compression condition (b = 0). The difference between the peak shear strength corresponding to the strain localization and the ultimate shear strength corresponding to plastic limit becomes obvious when the b value is near 0.4. The influence of bifurcation on the shear strength becomes weak when the loading direction changes from perpendicular to the bedding plane to parallel. The bifurcation analysis based on the proposed model gives out major principal strain and peak shear strength at the initiation of strain localization; the given results are consistent with experiments.     

Testing Study of Subcritical Crack Growth Mechanism During Water Rock Interaction

Subcritical crack growth plays an important role in evaluating the long-term stability of structures in rocks. By applying the constant-displacement-relax method, two groups of test specimens that one immersed in groundwater and the other in air were tested to get the parameters of subcritical crack growth in double torsion test. The relations of the stress intensity factor K _I versus the subcritical crack growth velocity V were obtained under the two different environments, and the behavior of subcritical crack growth was also analyzed. The results showed: the relations of lg K _I  ? lg V accorded with linear rules, which is in good agreement with Charles theory; Compared with specimens in nature state, the lg K _I  ? lg V curves of saturated water specimens locate top left comer of those of air specimens. The slope of curve is smaller, and the intercept is bigger, which shows that the water–rock interaction speeds up the subcritical crack growth. And A increases 2.9 × 10¹⁸ folds but n decreases from 85.12 to 40.83 because of the water–rock interaction. Meanwhile, the fracture toughness K _IC also decreases obviously from 2.55 in air to 2.26 in water due to water rock interaction. The testing results provide a basis for time-dependence of rock engineering stability.     

Estimation of source parameters and attenuation using digital waveforms of Al-Ays 2009 earthquake,Saudi Arabia

On 19 May 2009, an earthquake sequence of M _w?=?4.8 occurred at 25.20°N 37.76°E about 60 km onshore of the Red Sea coastline, Saudi Arabia. In the present study, the digital waveform data from the largest four events were used to estimate the source parameters and attenuation characteristics along the source-to-station path in the Arabian Shield. A grid search technique, combined with an assumption of circular source model, was applied to find the best-fit spectral amplitude over the space parameters: long period spectral level (Ω ₀), corner frequency (f ₀) and asymptotic high-frequency fall-off (γ). Consequently, the spectral parameters were used to estimate source parameters: seismic moment, fault radius (assumed circular rupture model) and stress drop. Seismic moments are founded to be within the range of 2.34E+14 to 2.83E+16 Nm and their corresponding moment magnitudes range from 3.5 to 4.8; the fault radius ranges from 369 to 1,498 m, and stress drops are observed in the range of 8.7 to 32.0 b. The spectral slopes beyond the corner frequency displayed ω ^?2.4 to ω ^?2.6 behaviours in contrast with Brune's source model of ω ^?2. This finding requires more detailed investigations on large data sets to distinguish the behaviour mechanism of the spectral slopes at high frequencies. By taking the ratio between observed and calculated spectra, the attenuation curves for P and S waves were derived along the source-to-station paths. The preliminarily results exhibited high quality factors of Q _α?=?3,883 and Q _β?=?3,530 for P and S waves, respectively. To this end, the ratio Q _β/Q _α is founded to be slightly less than unity indicating that the body waves from source-to-station paths crossed a crustal volume that is partially saturated with fluids causing lower attenuation effect on P waves than on S waves in the Arabian Shield.     

Effect of particle friction and polydispersity on the macroscopic stress–strain relations of granular materials

The macroscopic mechanical behavior of granular materials inherently depends on the properties of particles that compose them. Using the discrete element method, the effect of particle contact friction and polydispersity on the macroscopic stress response of 3D sphere packings is studied. The analytical expressions for the pressure, coordination number and fraction of rattlers proposed for isotropically deformed frictionless systems also hold when the interparticle coefficient of friction is finite; however, the numerical values of the parameters such as the jamming volume fraction change with varying microscopic contact and particle properties. The macroscopic response under deviatoric loading is studied with triaxial test simulations. Concerning the shear strength, our results agree with previous studies showing that the deviatoric stress ratio increases with particle coefficient of friction μ starting from a nonzero value for μ = 0 and saturating for large μ. On the other hand, the volumetric strain does not have a monotonic dependence on the particle contact friction. Most notably, maximum compaction is reached at an intermediate value of the coefficient of friction μ ≈ 0.3. The effect of polydispersity on the macroscopic stress–strain relationship cannot be studied independent of initial packing conditions. The shear strength increases with polydispersity when the initial volume fraction is fixed, but the effect of polydispersity is much less pronounced when the initial pressure of the packings is fixed. Finally, a simple hypoplastic constitutive model is calibrated with numerical test results following an established procedure to ascertain the relation between particle properties and material coefficients of the macroscopic model. The calibrated model is in good qualitative agreement with simulation results.     

In-situ Rock Spalling Strength near Excavation Boundaries

It is widely accepted that the in-situ strength of massive rocks is approximately 0.4 ± 0.1 UCS, where UCS is the uniaxial compressive strength obtained from unconfined tests using diamond drilling core samples with a diameter around 50 mm. In addition, it has been suggested that the in-situ rock spalling strength, i.e., the strength of the wall of an excavation when spalling initiates, can be set to the crack initiation stress determined from laboratory tests or field microseismic monitoring. These findings were supported by back-analysis of case histories where failure had been carefully documented, using either Kirsch’s solution (with approximated circular tunnel geometry and hence σ _max = 3σ ₁ ?σ ₃) or simplified numerical stress modeling (with a smooth tunnel wall boundary) to approximate the maximum tangential stress σ _max at the excavation boundary. The ratio of σ _max /UCS is related to the observed depth of failure and failure initiation occurs when σ _max is roughly equal to 0.4 ± 0.1 UCS. In this article, it is suggested that these approaches ignore one of the most important factors, the irregularity of the excavation boundary, when interpreting the in-situ rock strength. It is demonstrated that the “actual” in-situ spalling strength of massive rocks is not equal to 0.4 ± 0.1 UCS, but can be as high as 0.8 ± 0.05 UCS when surface irregularities are considered. It is demonstrated using the Mine-by tunnel notch breakout example that when the realistic “as-built” excavation boundary condition is honored, the “actual” in-situ rock strength, given by 0.8 UCS, can be applied to simulate progressive brittle rock failure process satisfactorily. The interpreted, reduced in-situ rock strength of 0.4 ± 0.1 UCS without considering geometry irregularity is therefore only an “apparent” rock strength.     

An RES-Based Model for Risk Assessment and Prediction of Backbreak in Bench Blasting

Most blasting operations are associated with various forms of energy loss, emerging as environmental side effects of rock blasting, such as flyrock, vibration, airblast, and backbreak. Backbreak is an adverse phenomenon in rock blasting operations, which imposes risk and increases operation expenses because of safety reduction due to the instability of walls, poor fragmentation, and uneven burden in subsequent blasts. In this paper, based on the basic concepts of a rock engineering systems (RES) approach, a new model for the prediction of backbreak and the risk associated with a blast is presented. The newly suggested model involves 16 effective parameters on backbreak due to blasting, while retaining simplicity as well. The data for 30 blasts, carried out at Sungun copper mine, western Iran, were used to predict backbreak and the level of risk corresponding to each blast by the RES-based model. The results obtained were compared with the backbreak measured for each blast, which showed that the level of risk achieved is in consistence with the backbreak measured. The maximum level of risk [vulnerability index (VI) = 60] was associated with blast No. 2, for which the corresponding average backbreak was the highest achieved (9.25 m). Also, for blasts with levels of risk under 40, the minimum average backbreaks (<4 m) were observed. Furthermore, to evaluate the model performance for backbreak prediction, the coefficient of correlation (R ²) and root mean square error (RMSE) of the model were calculated (R ² = 0.8; RMSE = 1.07), indicating the good performance of the model.     

Experimental Anelastic Strain Recovery Compliance of Three Typical Rocks

The experimental determination of anelastic strain recovery (ASR) compliances for three types of rocks (granite, marble, and sandstone) was performed in the laboratory. Preloading of specimens for uniaxial compression creep tests was at 50 % of the uniaxial compressive strength (UCS) for each rock type. We obtained the shear mode Jas(t) and volumetric mode Jav(t) ASR compliances and calculated the ratio of Jas(t) to Jav(t). The Kelvin model for rock rheology was then applied in numerical simulations and the results were in good agreement with the measured data for Jas(t) and Jav(t). These results showed that both the magnitude and rate of increase of the ASR compliances are strongly dependent on the rock type, and the values of the Jas(t)/Jav(t) ratio for a loading of 50 % of the UCS showed a trend leading to different constants for each of the three rock types. Further experimental and numerical analyses showed approximate power-law relationships between the ASR compliances at 50 % of UCS, and both the UCS and the tangential Young’s modulus at 50 % of UCS (E _t50). These relationships may be useful for the preliminary estimation of ASR compliances.     

非线性硬化与非线性软化的巷、隧道围岩塑性分析

在符合岩体实际的弹性、应变非线性硬化和具拐点的非线性软化的三段式光滑连接应力-应变关系条件下,求得静水压力下圆形巷、隧道围岩塑性硬化区和软化区的应力分布表达式。证明了所求得的巷道围岩径向和切向应力曲线在硬化区半径 和软化区半径 处均光滑连接,特别是所求得的切向应力 - 曲线没有尖峰向上的应力集中,此结果与日本学者久武胜保的试验结果一致。指出基于理想弹、塑性本构模型的Kastener解有明显缺陷,其 - 曲线尖峰向上的应力集中与实际不符。     

Strain Rate Dependency of Coarse Crystal Marble Under Uniaxial Compression: Strength,Deformation and Strain Energy

Strain rate during testing, uniaxial or triaxial, has important influence on the measured mechanical properties of rocks. Uniaxial compression tests were performed at nine pre-specified static-to-quasistatic strain rates (ranging from 1 × 10^?5 to 1 × 10^?1 s^?1) on coarse crystal marble. The aim is to gain deep insight into the influence of strain rate on characteristic stresses, deformation properties and conversion of strain energy of such rock. It is found that the strain rate of 5 × 10^?3 s^?1 is the threshold to delineate the failure modes the tested coarse marble behaves in. At a strain rate less than this threshold, single-plane shear and conjugate X-shaped shear are the main failure modes, while beyond this threshold, extensile and splitting failures are dominant. The stress for crack initiation, the critical stress for dilation, the peak stress, and Young’s modulus are all found to increase with strain rate, with an exception that the above stresses and modulus appear relatively low compared to the strain rate in the range of between 1 × 10^?4 and 5 × 10^?3 s^?1. The pre-peak absorbed strain energy, damage strain energy and elastic strain energy are found to increase with strain rate. In addition, the elastic strain energy stored before peak point favors brittle failure of the specimen, as the more stored elastic energy in the specimen, the stronger the fragmenting.     

Multiscale remote-sensing retrieval in the evapotranspiration of Haloxylon ammodendron in the Gurbantunggut desert, China

Transpiration, an essential component of surface evapotranspiration, is particularly important in the research of surface evapotranspiration in arid areas. The paper explores the spectral information of the arid vegetal evapotranspiration from a semi-empirical perspective by the measured data and the up-scaling method. The paper inverted the transpiration of Haloxylon ammodendron at the canopy, pixel and regional scales in the southern edge of the Gurbantunggut desert in Xinjiang, China. The results are as follows: at the canopy scale, the optimal exponential model of the sap flow rate based on the hyperspectrum is y = 0.0015e^3.8922x , R ² = 0.806. At the pixel scale, there was a good linear relationship between the sap flow and the SR index, with a relationship of y = ?1197.38x ³ + 1048.43x ² ? 305.47x + 455.15, R ² = 0.845. At the regional scale, based on the optimal exponential model and the EO-1 Hyperion remote-sensing data, the transpiration of the study area was inverted. Comparing the results of the SEBAL and SEBS models, the errors of the simulation results were 7.78 and 8.80 %. The paper made full use of the knowledge flow at different scales, bridging the scale difference in canopy and remote-sensing images to avoid the information bottleneck in the up-scaling. However, there are many constraints in the data acquirement, the efficiency of the models, the endmembers determination, the temporal–spatial up-scaling, and the accuracy assessment, which would be improved in future studies.     

单轴拉伸岩样破坏过程及尺寸效应数值模拟

由于从实验及理论角度研究岩样单轴拉伸条件下的破坏全过程及尺寸效应难度都很大。因此采用拉格朗日元法来研究这些问题。在峰值强度之前后,岩石材料的本构模型分别取为线弹性及拉破坏线性应变软化模型。为了使拉伸塑性区不出现在试样的端部,在试样的两侧面中部预制了2个凹槽。数值模拟结果表明,全程拉应力-拉应变曲线分为峰前和峰后阶段。在接近峰值的峰前阶段,由于两凹槽附近具有明显的拉应力集中现象,拉伸塑性区最先出现在两凹槽附近。随着轴向拉应变的增加,发生拉伸破坏的单元的数目增加,新发生拉伸破坏的单元越来越接近试样的中心,直到两块拉伸塑性区在应变软化阶段贯通。两凹槽连线上各单元拉应力的分布呈现3个阶段,“澡盆型”（“U型”）阶段,“双峰型”（“M型”）阶段及“单峰型”（“Π型”）阶段。“澡盆型”阶段对应于全程拉应力-拉应变曲线的弹性阶段。“双峰型”阶段及“单峰型”阶段对应于全程拉应力-拉应变曲线的非弹性阶段（包括峰值强度之前的一小段,即应变硬化阶段及峰后的应变软化阶段）。增加试样高度及降低试样宽度,拉应力-拉应变曲线的软化段变得越来越陡峭,因而试样越容易发生失稳破坏。由于试样宽度较大时,试样内部的单元并非处于单向拉应力状态,因此,增加试样宽度,全程拉应力-拉应变曲线的峰值强度增加。当试样宽度较小时,从出现塑性区,到塑性区贯通所需要的时间步较小,或应变范围较窄。这说明试样的脆性较强,前兆不明显。前兆不明显的脆性破坏对应常见的是洞室岩爆、冲击地压及地震等灾害。     

Re-evaluating the US Geological Survey’s pumping tests (1967) in the Punjab region of Pakistan for use in groundwater studies

In 1967, the US Geological Survey (USGS) published the results of 141 pumping tests carried out throughout the Pakistani Punjab to establish representative hydraulic parameters of its large aquifer. Many authors have since concluded that the USGS had over-estimated the horizontal hydraulic conductivity (k _r) by 25–100 %, leaving vertical anisotropy and aquifer depth unresolved. No test wells have ever been drilled below 450 m to reach the base of the aquifer, although petroleum explorations mention depths between 1,500 and 4,500 m. After comparison and re-evaluation of all related papers, this study concludes that the USGS interpretation was correct, that its hydraulic values still stand without change, and that the USGS’s applied distance drawdown interpretation is valid to prevent influence of partial penetration on the results. This study also uniquely resolved vertical anisotropy and aquifer thickness by using early- and late-time drawdowns separately and proper scaling of the coordinates, which has often been omitted. With appropriate scaling, all interpretations match the data. The representative hydraulic aquifer values are: k _r?=?65 m/d, vertical anisotropy k _r/k _z?=?25 and aquifer depth 500–1,500 m. The conclusion is that these values can be used, at least as first estimates, for groundwater studies in the Pakistani Punjab.     

Investigation of rock quality of Shirinrud dam site by engineering seismology

At this paper, we studied about the rock quality of Shirinrud dam site by engineering seismology. Shirinrud dam site is located 80 km far from Kerman and 18 km far from Hojadk village. The dam and its constructions are established in the Bidu Formation which consists of seven rock units, and the refraction profiles were surveyed on Jb_3/2, Jb₄, and Jb₅ rock units. To evaluate the rock mass quality and basement topography at this site, nine refraction seismic profiles by primary waves and two refraction seismic profiles by secondary waves were surveyed. We used some methods such as Palmer method, the reciprocal method, plus–minus method, etc. to process and interpret data. Based on investigations, primary wave velocity in unit Jb_3/2 varies between 2,100 and 2,200 m/s, in unit Jb₄ is between 2,100 and 4,200 m/s, and in unit Jb₅ is between 2,500 and 3,000 m/s. The Q values on these three units are 0.05, 1.2, and 1.9, and the rock mass rating (RMR) values are 27.1, 40.5, and 33.5, respectively. With respect to wave velocity, Q, and RMR values, the units Jb_3/2, Jb₄, and Jb₅ are evaluated as very weak, intermediate, and weak, respectively.     

Transport properties of olivine grain boundaries from electrical conductivity experiments

Grain boundary processes contribute significantly to electronic and ionic transports in materials within Earth’s interior. We report a novel experimental study of grain boundary conductivity in highly strained olivine aggregates that demonstrates the importance of misorientation angle between adjacent grains on aggregate transport properties. We performed electrical conductivity measurements of melt-free polycrystalline olivine (Fo₉₀) samples that had been previously deformed at 1200 °C and 0.3 GPa to shear strains up to γ?=?7.3. The electrical conductivity and anisotropy were measured at 2.8 GPa over the temperature range 700–1400 °C. We observed that (1) the electrical conductivity of samples with a small grain size (3–6 µm) and strong crystallographic preferred orientation produced by dynamic recrystallization during large-strain shear deformation is a factor of 10 or more larger than that measured on coarse-grained samples, (2) the sample deformed to the highest strain is the most conductive even though it does not have the smallest grain size, and (3) conductivity is up to a factor of ~?4 larger in the direction of shear than normal to the shear plane. Based on these results combined with electrical conductivity data for coarse-grained, polycrystalline olivine and for single crystals, we propose that the electrical conductivity of our fine-grained samples is dominated by grain boundary paths. In addition, the electrical anisotropy results from preferential alignment of higher-conductivity grain boundaries associated with the development of a strong crystallographic preferred orientation of the grains.