Shock metamorphism of plagioclase and amphibole (Experimental data)

The shock metamorphism of plagioclase and amphibole of various chemical compositions from amphibolite and granulite facies schists was studied in experiments with shock wave loading of samples in steel recovery ampoules of plane geometry. A maximum shock pressure was reached after a few circulations of waves in the sample (stepwise shock wave compression) and varied within 26–52 GPa. The recovered samples were examined by the methods of scanning electron microscopy and microprobe and X-ray phase analysis. It is established that an increase in the F, Ti, and K concentrations in amphibole and a decrease in the Ca concentration in plagioclase make these minerals more stable with respect to shock waves. It is shown that the migration of some chemical elements, starting already at the solid phase stage of transformation in plagioclase and amphibole, is intensified at the stage of melting. It is established that isotropization of plagioclase occurs through two different mechanisms. At relatively low pressures, it is caused by the fragmentation of substance at the microlevel and is accompanied by the formation of maskelynite, a typical mineral of meteorites and astrobleme rocks. At higher pressures, isotropization is associated with melting-induced amorphization.     

Shock reequilibration of fluid inclusions in Coconino sandstone from Meteor Crater,Arizona

This study examines the effects of natural shock metamorphism on fluid inclusions trapped in porous sedimentary target rocks and compares these results to previous experimental work on single crystal quartz. Samples of shock metamorphosed Coconino sandstone were collected from Barringer Meteorite Crater (Meteor Crater, Arizona) and classified based on their shock features into the six shock stages described by Kieffer [S.W. Kieffer, 1971. Shock metamorphism of the Coconino sandstone at Meteor Crater, Arizona, Journal of Geophysical Research 76, 5449-5473.]. The frequency of two-phase fluid inclusions decreases dramatically from unshocked samples of Coconino sandstone through shock stages 1a, 1b, and 2. No two-phase fluid inclusions were observed in shock stage 3 or 4 samples. However, the total number of grains containing fluid inclusions remains approximately the same for shock stages 1a–2, suggesting that two-phase fluid inclusions reequilibrated during impact to form single-phase inclusions. In shock stages 3 and 4, the total number of inclusions also decreases, indicating that at these higher shock pressures fluid inclusions are destroyed by plastic deformation and phase changes within the host mineral. Entrained quartz grains within a shock stage 5 sample contain two-phase inclusions, emphasizing the short duration of melting associated with the impact and the heterogeneous nature of impact processes. These results are similar to those observed in single-crystal experiments, although inclusions survive to slightly higher shock pressures in samples of naturally shocked Coconino sandstone. Results of this study suggest that the rarity of fluid inclusions in meteorites does not preclude the presence of fluids on meteorite parent bodies. Instead, fluid inclusions trapped during alteration events may have been destroyed due to shock processing. In addition, loss of fluids from inclusion vesicles along fractures and microcracks may lead to shock devolatilization, even in unsaturated target rocks.     

Surfatron particle acceleration in the near-Earth shock wave

A theoretical model, within the scope of which the process of surfatron particle acceleration (surfing) at the shock front was studied, is proposed to describe the near-Earth shock wave structure. The most complex shock front geometry, when it is assumed that the magnetic field vectors and the velocity of the plasma flow incident on the front are oriented arbitrarily with respect to the front plane, has been considered in detail. The theoretical ultimate energies of protons, accelerated due to surfing at the shock front, agree with the observational data.     

The mechanism and dynamics of rock fracture upon mechanical impact and electric discharge

The mechanism and dynamics of the deformation and fracture of quartz, granite, and marble samples under the striker blow on their surface and electric discharge inside them are studied by the fractoluminescence (FL), electromagnetic (EME), and acoustic emission (AE) methods with 10-ns resolution. The impact excites a forced deformation wave with a velocity within 0.8 to 2 km/s depending on the mineral. The atomic bonds rupture and microcracks are formed at the nodes of the wave, which leads to the emergence of the FL flashes and disruption of the time dependences of EME. Based on the intensity of the flashes, the dimensions of microcracks are estimated to vary from 2 to 70 µm depending on the mineral. In turn, the emergence of microcracks initiates additional deformation waves.The discharge inside the studied samples excites a pressure shock wave which transforms into the tension wave after reflection from the surface. According to the analysis of FL spectra, this leads to the breakdown of the rocks into positively charged ions and electrons. The shock wave velocity in granites is measured at 4.8 km/s, which is close to the velocity of the longitudinal acoustic vibrations ~5 km/s. The microcracks in the rock have not enough time to form with this loading velocity. It is supposed that the shock wave stretches the deformed interatomic bonds at the dislocation nuclei in the crystal lattices of the minerals up to their breakdown into positively charged ions.     

Theory and application of equivalent transformation relationships between plane wave and spherical wave

Based on the governing equations and the equivalent models,we propose an equivalent transformation relationships between a plane wave in a one-dimensional medium and a spherical wave in globular geometry with radially inhomogeneous properties.These equivalent relationships can help us to obtain the analytical solutions of the elastodynamic issues in an inhomogeneous medium.The physical essence of the presented equivalent transformations is the equivalent relationships between the geometry and the material properties.It indicates that the spherical wave problem in globular geometry can be transformed into the plane wave problem in the bar with variable property fields,and its inverse transformation is valid as well.Four different examples of wave motion problems in the inhomogeneous media are solved based on the presented equivalent relationships.We obtain two basic analytical solution forms in Examples Ⅰ and Ⅱ,investigate the reflection behavior of inhomogeneous half-space in Example Ⅲ,and exhibit a special inhomogeneity in Example Ⅳ,which can keep the traveling spherical wave in constant amplitude.This study implies that our idea makes solving the associated problem easier.     

Hugoniot equation of state of olivine and its geodynamic implications

Large olivine samples were hot-pressed synthesized for shock wave experiments. The shock wave experiments were carried out at pressure range between 11 and 42 GPa. Shock data on olivine sample yielded a linear relationship between shock wave velocity D and particle velocity u described by D=3.56(?0.13)+2.57(?0.12)u. The shock temperature is determined by an energy relationship which is approximately 790°C at pressure 28 GPa. Due to low temperature and short experimental duration, we suggest that no phase change occurred in our sample below 30 GPa and olivine persisted well beyond its equilibrium boundary in metastable phase. The densities of metastable olivine are in agreement with the results of static compression. At the depth shallower than 410 km, the densities of metastable olivine are higher than those of the PREM model, facilitating cold slab to sink into the mantle transition zone. However, in entire mantle transition zone, the shock densities are lower than those of the PREM model, hampering cold slab to flow across the "660 km" phase boundary.     

Compressional wave events in the dawn plasma sheet observed by Interball-1

Compressional waves with periods greater than 2 min (about 10–30 min) at low geomagnetic latitudes, namely compressional Pc5 waves, are studied. The data set obtained with magnetometer MIF-M and plasma analyzer instrument CORALL on board the Interball-1 are analyzed. Measurements performed in October 1995 and October 1996 in the dawn plasma sheet at −30 R_E ≤ X_GSM and |Z_GSM| ≤ 10 R_E are considered. Anti-phase variations of magnetic field and ion plasma pressures are analyzed by searching for morphological similarities in the two time series. It is found that longitudinal and transverse magnetic field variations with respect to the background magnetic field are of the same order of magnitude. Plasma velocities are processed for each time period of the local dissimilarity in the pressure time series. VeloCity disturbances occur mainly transversely to the local field line. The data reveal the rotation of the veloCity vector. Because of the field line curvature, there is no fixed position of the rotational plane in the space. These vortices are localized in the regions of anti-phase variations of the magnetic field and plasma pressures, and the vortical flows are associated with the compressional Pc5 wave process. A theoretical model is proposed to explain the main features of the nonlinear wave processes. Our main goal is to study coupling of drift Alfven wave and magnetosonic wave in a warm inhomogeneous plasma. A vortex is the partial solution of the set of the equations when the compression is neglected. A compression effect gives rise to a nonlinear soliton-like solution.     

Shock and post-shock temperatures in an ice–quartz mixture: implications for melting during planetary impact events

Melting of H₂O ice during planetary impact events is a widespread phenomenon. On planetary surfaces, ice is often mixed with other materials; yet, at present, the partitioning of energy between the components of a shocked mixture is still an open question in the shock physics community. Knowledge of how much energy is partitioned into the ice component is necessary to predict and interpret a wide range of processes, including shock-induced melting and chemistry. In this work, we construct a conceptual framework for the thermodynamic pathways of the components in a shocked hydrodynamic mixture by defining three broad regimes based on the characteristic length scale of the mixture compared to the thickness of the shock front: (1) small length scale mixtures where pressure and temperature equilibrate immediately behind the shock front; (2) intermediate length scales where pressure but not thermal equilibration is achieved behind the shock front; and (3) long length scales where pressure equilibration requires multiple shock wave reflections. We conduct shock wave experiments, reaching pressures from 8 to 23 GPa, in an H₂O ice–SiO₂ quartz mixture in the intermediate length scale regime. In each experiment, all the parameters required to address the question of energy partitioning were determined: the shock velocity in the mixture, the shock front thickness, and the shock and post-shock temperatures of the H₂O component. The measured pressure is in agreement with the bulk compressibility of the mixture. The shock and post-shock temperatures of the H₂O component indicate that the ice was shocked close to the principal Hugoniot. Therefore, in the intermediate length scale regime, the partitioning of shock energy is defined initially by the Hugoniots of the components at the equilibrated pressure. We discuss energy partitioning in mixtures over the wide range of length and time scales encountered during planetary impact events and identify the current challenges in calculating the volume of melted ice. In some cases, the criteria for shock-induced melting of ice in a mixture are the same as for pure ice.     

不同温压条件下蛇纹岩和角闪岩中波速与衰减的各向异性

在实验室中研究了蛇纹岩和角闪岩样品在不同温压条件下的纵、横波速度和Q值．这两种岩样对应的主要组成矿物叶蛇纹石和普通角闪石都具有很强的晶格优选方位（LPO）．随着围压的增加,波速和Q值均增大,但是在相互正交的三个方向上（垂直或平行于层理面及线理方向）增大的速度并不相同,这与微裂隙的逐渐闭合密切相关．在600MPa的围压下升高温度直到600℃以上,由于微裂隙的热扩张受到约束,波速和Q值下降幅度很小．观测到的波速和Q值的各向异性具有不同的机理,波速各向异性主要与定向分布的微裂隙和主要矿物的LPO等构造因素有关;高围压下纵波Q值各向异性与速度各向异性正好相反,可能是由于形成层理面的定向排列的平板状矿物晶体沿不同方向边界之间接触程度不同造成的．     

A note on three-dimensional scattering and diffraction by a hemispherical canyon–I: Vertically incident plane P-wave

The three-dimensional scattering by a hemi-spherical canyon in an elastic half-space subjected to seismic plane and spherical waves has long been a challenging boundary-value problem. It has been studied by earthquake engineers and strong-motion seismologists to understand the amplification effects caused by surface topography. The scattered and diffracted waves will, in all cases, consist of both longitudinal (P-) and shear (S-) shear waves. Together, at the half-space surface, these waves are not orthogonal over the infinite plane boundary of the half-space. Thus, to simultaneously satisfy both zero normal and shear stresses on the plane boundary numerical approximation of the geometry and/or wave functions were required, or in some cases, relaxed (disregarded). This paper re-examines this boundary-value problem from the applied mathematics point of view, and aims to redefine the proper form of the orthogonal spherical-wave functions for both the longitudinal and shear waves, so that they can together simultaneously satisfy the zero-stress boundary conditions at the half-space surface. With the zero-stress boundary conditions satisfied at the half-space surface, the most difficult part of the problem will be solved, and the remaining boundary conditions at the finite canyon surface will be easy to satisfy.     

A METHOD OF TACKLING SOME SEISMIC PROBLEMS AND ITS APPLICATIONS *

The theory is explained and practical applications are shown for a numerical procedure in seismology. Particularly the problems concerning the generation of waves under the action of external pressures, and their propagation, in non-homogeneous, both elastic and absorbing media, have been carried out. These problems have been assumed mono-dimensional and refer to plane and spherical waves. The procedure is based on the solution, by means of series, of the wave differential equation, non-homogeneous, and with non-constant coefficients. It is a direct numerical method whose advantage is, mainly, the possibility of tackling, without great difficulties, problems regarding non-homogeneous elastic and absorbing media. On the contrary the methods which require the theoretical expression of the solutions by means of formulae, generally, present conceptual and numerical difficulties. As examples of application of this procedure, the following cases have been carried out by means of numerical calculations. a) Propagation of a wave, initially of symmetrical shape, in a viscoelastic medium; from the results it appears that the wave propagates without losing its symmetry, i.e. without sensible dispersion. A theoretical analysis has been carried out to justify this result, showing that the dispersion in viscoelastic media is noticeable only for relatively high frequencies. It seems that the practical absence of dispersion in field experiments do not exclude necessarily the viscoelastic character of absorption. b) Generation of plane waves under the effect of a uniform pressure distributed on the plane surface of a medium. The way the length and the shape of the generated wave depends not only on the type of pressure acting on the surface but also on the near surface impedance variations has been studied. c) Generation of a spherical wave under the action of a pressure in a spherical hole. The examples treated show how the length and shape of the wave depends on the radius of the hole. Particularly the frequencies of the wave spectrum are proportional to this radius, for a given type of pressure acting in the hole. The characteristics of this procedure would also permit the study of media for which the stress-strain relations are not univocal and linear (non linear absorption). This study, interesting for the wave propagation in rocks, is worth while to be carried out in a special paper.     

Analytical wavefront curvature correction to plane‐wave reflection coefficients for a weak‐contrast interface

Most amplitude versus offset (AVO) analysis and inversion techniques are based on the Zoeppritz equations for plane‐wave reflection coefficients or their approximations. Real seismic surveys use localized sources that produce spherical waves, rather than plane waves. In the far‐field, the AVO response for a spherical wave reflected from a plane interface can be well approximated by a plane‐wave response. However this approximation breaks down in the vicinity of the critical angle. Conventional AVO analysis ignores this problem and always utilizes the plane‐wave response. This approach is sufficiently accurate as long as the angles of incidence are much smaller than the critical angle. Such moderate angles are more than sufficient for the standard estimation of the AVO intercept and gradient. However, when independent estimation of the formation density is required, it may be important to use large incidence angles close to the critical angle, where spherical wave effects become important. For the amplitude of a spherical wave reflected from a plane fluid‐fluid interface, an analytical approximation is known, which provides a correction to the plane‐wave reflection coefficients for all angles. For the amplitude of a spherical wave reflected from a solid/solid interface, we propose a formula that combines this analytical approximation with the linearized plane‐wave AVO equation. The proposed approximation shows reasonable agreement with numerical simulations for a range of frequencies. Using this solution, we constructed a two‐layer three‐parameter least‐squares inversion algorithm. Application of this algorithm to synthetic data for a single plane interface shows an improvement compared to the use of plane‐wave reflection coefficients.     

大麻坪橄榄石冲击波实验研究与高压态物态方程参数的确定

利用冲击波进行的动高压实验对研究超高压下的物质性质也是十分有效的.本文报告了对大麻坪采集的橄榄石进行压力在10~45GPa范围的冲击波动高压实验结果.结合前人等温冷压实验结果,确定了实验过程的温度,对于冲击波实验压力从10GPa变化到30GPa时,温度在摄氏几十度到摄氏800多度之间.测量了岩石超高压下密度变化,在3.627~4.009g·cm~(-3)之间.通过回收实验和确定的温度,表明小于30GPa压力实验条件下,没有发生相变过程.同时也确定了状态方程的参数.最后,指出了实验结果在上地幔地球内部物质运动过程的含义,即冷板块中的亚临界橄榄岩可以存在地幔转换带中.     

An Integrated Study on Physical Properties of a KTB Gneiss Sample and Marble from Portugal: Pressure Dependence of the Permeability and Frequency Dependence of the Complex Electrical Impedance

— Pressure-induced variations in pore geometry were studied on dry- and fluid- saturated samples by means of electrical impedance spectroscopy and permeability measurements. Hydrostatic pressures (up to 120 MPa) and uniaxial pressures (up to failure) were applied. Hydrostatic pressures reduce the aspect ratio of cracks and thus cause a decrease of permeability and electrical bulk conductivity. The opposite was observed in uniaxial pressure experiments where new cracks were formed and consequently permeability and electrical conductivity were increased. More specific informations of these generated observations were derived from the interpretation of the frequency dispersion of the complex electrical conductivity. This least-squares-refinement considers electrochemical interactions between the fluid pore electrolyte and the inner surface of the sample, thus providing informations on the pore geometry and pressure-induced variations. Consequently changes in aspect ratio, size and geometry of the pore system can be detected by means of impedance spectroscopy.     

岩石超高压状态参数的实验研究

岩石的超高压状态参数对研究地球动力学问题和岩石爆破工程问题具有十分重要的意义.本实验应用平面波炸药透镜驱动飞板撞击靶板产生超高压冲击波,并用探针测量冲击波在岩石试件中的传播速度和粒子速度,再由 Rankine-Hugoniot 方程得到岩石的密度与压力之间的关系.所研究的岩石为石灰岩和花岗岩.实验的最高压力超过70GPa.得到两种岩石的超高压状态方程是:石灰岩————=3.22+0.016p;花岗岩————=3.20+0.024p.其中的单位为 g/cm3,p 的单位为 GPa.      

Modeling the Magnetic Anomaly of the Bosumtwi (Ghana) Complex Meteorite Crater by Taking Into Account the Impact Demagnetization and Morphological Features

Izvestiya, Physics of the Solid Earth - Abstract—The formation of impact craters on the Earth’s surface is accompanied by the effect of shock waves on rocks. The shock wave compression...     

涪陵地区井下页岩岩芯裂缝体密度与声学性质关系实验研究

本研究对涪陵地区井下产气的龙马溪组页岩岩芯,采用三轴压机压裂制造裂缝,用工业CT扫描压裂前后岩芯,应用图像识别技术统计出裂缝体密度,又用超声脉冲透射法测定样品破裂前后裂缝方向上的纵波、单偏振方向与裂缝不同夹角的横波,来研究裂缝对页岩岩芯声学性质的影响.实验结果表明:压裂后样品的纵波速度略微降低,只有含较多内部裂隙的150#样品纵波速度减小幅度明显.压裂前后样品的纵波波形差别不大,纵波主频随裂缝体密度呈下降趋势,即压裂后纵波频谱主频向低频端移动.压裂前横波速度随自身与裂隙方位角变化而变化,与0°和180°相比,在45°和135°时略微减小,在90°时速度降低幅度最为明显并且发生相位反转.典型样品的横波主频随偏振方向与裂缝夹角的增大而逐渐向低频移动;压裂后,横波频谱杂乱,出现多处局部峰值,速度和主频较压裂前更低,平均横波波速随裂缝体密度呈明显减小趋势.平均纵横波速比随裂缝密度呈近线性增加,表明其与裂缝体密度有较强相关性.     

Velocity anisotropy and attenuation of shale in under- and overpressured conditions

The seismic velocity and attenuation of fully saturated shales were measured for the first time under overpressured conditions, using the ultrasonic reflection technique. Shale cores from naturally overpressured horizons in the North Sea were tested in the laboratory, at confining and pore pressures relevant to in situ conditions.
A single-frequency tone-burst pulse wave was used to determine the seismic wave velocities and quality factors of the shale samples, with errors less than 0.3% and 0.1 dB/cm, respectively, at a frequency of 0.75 MHz. Sample length changes with varying confining and pore pressure were measured and the pore pressure equilibration time was monitored for each sample.
The anisotropy of the seismic attributes ( V _p, V _s, Q _p and Q _s) was determined over a range of differential pressures from 5 to 60 MPa, with respect to the predominant foliation. The ultrasonic velocity data followed a transversely isotropic pattern depending on the direction of wave propagation with respect to the laminations. The Poisson's ratio was found to rise by 5% as the shale material progressed from a normally pressured to an overpressured state. The quality factor ( Q ) characteristics were interpreted in terms of pore geometry and connectivity as well as the directional permeability of the transversely isotropic shale material. The results were converted to bulk and shear loss modulus defects, and a positive bulk loss was observed for waves propagating perpendicular to the lamination plane even above differential pressures of 20 MPa. This indicates different levels of Biot-flow and squirt-flow attenuation mechanisms acting within the shale structure, depending on the wave propagation and vibration directions.     

山西裂谷带地壳岩石波速的研究

通过实验建立了岩石成分与波速的关系,波速和岩石的各向异性,矿物相变在波速曲线中的反映以及实验温度和压力对波速的影响,尤其得到了高温下的新资料.用实验数据结合地球物理、地质资料建立了山西裂谷带与相邻太行山隆起区地壳各层圈的物质组成,对比说明,两构造单元下地壳的岩石组成是不同的.最后分析了盆地下地壳中部的组成、水和环境条件,认为低速带是由岩石部分熔融造成的.     

孔隙压力对含随机缺陷岩石破坏过程及全部变形特征的影响

对于平面应变压缩条件下含有随机缺陷的岩样,利用FLAC研究了孔隙压力对岩样破坏过程、全部变形特征及前兆的影响。以前编写的若干FISH函数,被用于生成缺陷和计算轴向、侧向、体积应变及侧向应变与轴向应变比值的负值（计算得到的泊松比）。在峰值应力之后,密实岩石单元服从线性应变软化行为及随后的理想塑性行为,而材料缺陷呈现理想塑性行为。当孔隙压力较高时,应力-侧向应变曲线具有一平台;破坏的前兆更明显;变形后岩样的体积总是大于原始体积;在初始加栽阶段、均匀变形阶段及峰后变形阶段,由于明显的侧向膨胀,计算得到的泊松比远大于0．5。当孔隙压力较低时,在峰值应力之前,变形后岩样的体积小于原始体积,体积扩容出现于峰值应力之后,引起了负的体积应变。利用广义虎克定律,解释了平面应变弹性状态下数值结果的合理性。对岩样进行带状区域扫描后,确认随机缺陷的初始分布与岩样的最终破坏形态紧密相关。