Petrophysical constraints on the seismic properties of the Kaapvaal craton mantle root

Abstract. We calculated the seismic properties of 47 mantle xenoliths from 9 kimberlitic pipes in the Kaapvaal craton based on their modal composition, the crystal-preferred orientations (CPO) of olivine, ortho- and clinopyroxene, and garnet, the Fe content of olivine, and the pressures and temperatures at which the rocks were equilibrated. These data allow constraining the variation of seismic anisotropy and velocities within the cratonic mantle. The fastest P and S2 wave propagation directions and the polarization of fast split shear waves (S1) are always subparallel to olivine [100] axes of maximum concentration, which marks the lineation (fossil flow direction). Seismic anisotropy is higher for high olivine contents and stronger CPO. Maximum P wave azimuthal anisotropy (AVp) ranges between 2.5 and 10.2% and the maximum S wave polarization anisotropy (AVs), between 2.7 and 8%. Changes in olivine CPO symmetry result in minor variations in the seismic anisotropy patterns, mainly in the apparent isotropy directions for shear wave splitting. Seismic properties averaged over 20 km-thick depth sections are, therefore, very homogeneous. Based on these data, we predict the anisotropy that would be measured by SKS, Rayleigh (SV) and Love (SH) waves for five endmember orientations of the foliation and lineation. Comparison to seismic anisotropy data from the Kaapvaal shows that the coherent fast directions, but low delay times imaged by SKS studies, and the low azimuthal anisotropy with with the horizontally polarized S waves (SH) faster than the vertically polarized S wave (SV) measured using surface waves are best explained by homogeneously dipping (45°) foliations and lineations in the cratonic mantle lithosphere. Laterally or vertically varying foliation and lineation orientations with a dominantly NW–SE trend might also explain the low measured anisotropies, but this model should also result in backazimuthal variability of the SKS splitting data, not reported in the seismological data. The strong compositional heterogeneity of the Kaapvaal peridotite xenoliths results in up to 3% variation in density and in up to 2.3% variation of Vp, Vs, and Vp / Vs ratio. Fe depletion by melt extraction increases Vp and Vs, but decreases the Vp / Vs ratio and density. Orthopyroxene enrichment due to metasomatism decreases the density and Vp, strongly reducing the Vp / Vs ratio. Garnet enrichment, which was also attributed to metasomatism, increases the density, and in a lesser extent Vp and the Vp / Vs ratio. Comparison of density and seismic velocity profiles calculated using the xenoliths' compositions and equilibration conditions to seismological data in the Kaapvaal highlights that (i) the thickness of the craton is underestimated in some seismic studies and reaches at least 180 km, (ii) the deep sheared peridotites represent very local modifications caused and oversampled by kimberlites, and (iii) seismological models probably underestimate the compositional heterogeneity in the Kaapvaal mantle root, which occurs at a scale much smaller than the one that may be sampled seismologically.

Download Full-text

Petrophysical constraints on the seismic properties of the Kaapvaal craton mantle root

Solid Earth Discussions ◽

10.5194/sed-5-963-2013 ◽

2013 ◽

Vol 5 (2) ◽

pp. 963-1005 ◽

Cited By ~ 2

Author(s):

V. Baptiste ◽

A. Tommasi

Keyword(s):

Seismic Anisotropy ◽

Seismic Velocity ◽

Body Wave ◽

Azimuthal Anisotropy ◽

Mantle Xenoliths ◽

P Wave ◽

Kaapvaal Craton ◽

Compositional Heterogeneity ◽

P Waves ◽

Seismic Properties

Abstract. We calculated the seismic properties of 47 mantle xenoliths from 9 kimberlitic pipes in the Kaapvaal craton based on their modal composition, the crystal preferred orientations (CPO) of olivine, ortho- and clinopyroxene, and garnet, the Fe content of olivine, and the pressures and temperatures at which the rocks were equilibrated. These data allow constraining the variation of seismic anisotropy and velocities with depth. The fastest P wave and fast split shear wave (S1) polarization direction is always close to olivine [100] maximum. Changes in olivine CPO symmetry result in minor variations in the seismic anisotropy patterns. Seismic anisotropy is higher for high olivine contents and stronger CPO. Maximum P waves azimuthal anisotropy (AVp) ranges between 2.5 and 10.2% and S waves polarization anisotropy (AVs) between 2.7 and 8%. Seismic properties averaged in 20 km thick intervals depth are, however, very homogeneous. Based on these data, we predict the anisotropy that would be measured by SKS, Rayleigh (SV) and Love (SH) waves for 5 end-member orientations of the foliation and lineation. Comparison to seismic anisotropy data in the Kaapvaal shows that the coherent fast directions, but low delay times imaged by SKS studies and the low azimuthal anisotropy and SH faster than SV measured using surface waves may only be consistently explained by dipping foliations and lineations. The strong compositional heterogeneity of the Kaapvaal peridotite xenoliths results in up to 3% variation in density and in up to 2.3% of variation Vp, Vs and the Vp/Vs ratio. Fe depletion by melt extraction increases Vp and Vs, but decreases the Vp/Vs ratio and density. Orthopyroxene enrichment decreases the density and Vp, but increases Vs, strongly reducing the Vp/Vs ratio. Garnet enrichment increases the density, and in a lesser manner Vp and the Vp/Vs ratio, but it has little to no effect on Vs. These compositionally-induced variations are slightly higher than the velocity perturbations imaged by body-wave tomography, but cannot explain the strong velocity anomalies reported by surface wave studies. Comparison of density and seismic velocity profiles calculated using the xenoliths' compositions and equilibrium conditions to seismological data in the Kaapvaal highlights that: (i) the thickness of the craton is underestimated in some seismic studies and reaches at least 180 km, (ii) the deep sheared peridotites represent very local modifications caused and oversampled by kimberlites, and (iii) seismological models probably underestimate the compositional heterogeneity in the Kaapvaal mantle root, which occurs at a scale much smaller than the one that may be sampled seismologically.

Download Full-text

A database of plagioclase crystal preferred orientations (CPO) and microstructures – implications for CPO origin, strength, symmetry and seismic anisotropy

Solid Earth Discussions ◽

10.5194/sed-5-1191-2013 ◽

2013 ◽

Vol 5 (2) ◽

pp. 1191-1257 ◽

Cited By ~ 1

Author(s):

T. Satsukawa ◽

B. Ildefonse ◽

D. Mainprice ◽

L. F. G. Morales ◽

K. Michibayashi ◽

...

Keyword(s):

Plastic Deformation ◽

Seismic Anisotropy ◽

P Wave ◽

S Wave ◽

Seismic Properties ◽

Strong Point ◽

Pole Figures ◽

Significant Difference ◽

Sampling Statistics ◽

P Type

Abstract. This study presents a unique database of 170 plagioclase Crystallographic Preferred Orientations (CPO) of variously deformed gabbroic rocks. The CPO characteristics as a function of the deformation regime (magmatic or crystal-plastic) are outlined and discussed. The studied samples are dominantly from slow- and fast-spread present-day ocean crust, as well as from the Oman ophiolite. Plagioclase is the dominant mineral phase in the studied samples. Plagioclase CPOs are grouped in three main categories: Axial-B, a strong point alignment of (010) with a girdle distribution of [100]; Axial-A, a strong point maximum concentration of [100] with parallel girdle distributions of (010) and (001); and P-type, point maxima of [100], (010), and (001). A majority of CPO patterns are Axial-B and P-type, in samples showing either magmatic or crystal-plastic deformation textures. Axial-A CPOs are less common; they represent 21% of the samples deformed by crystal-plastic flow. Although fabric strength (ODF J-index) does not show any consistent variation as a function of the CPO patterns, there is a significant difference in the relationship between the ODF and pole figures J-indices; the magmatic type microstructures have high (010) pole figures J-indices, which increase linearly with ODF J-index, whereas the high [100] pole figures J-indices of plastically deformed samples vary in a more scattered manner with ODF J-index. The multistage nature of plastic deformation superposed on a magmatic structure compared with magmatic flow, and the large number of possible slip-systems in plagioclase probably account for these differences. Calculated seismic properties (P wave and S wave velocities and anisotropies) of plagioclase aggregates show that anisotropy (up to 12% for P wave and 14% for S wave) tends to increase as a function of ODF J-index. In comparison with the olivine 1998 CPO database, the magnitude of P wave anisotropy for a given J-index is much less than olivine, whereas it is similar for S wave anisotropy. Despite a large variation of fabric patterns and geodynamic setting, seismic properties of plagioclase-rich rocks have similar magnitudes of anisotropy. There is a small difference in the aggregate elastic symmetry, with magmatic microstructures having higher orthorhombic and hexagonal components, whereas plastic deformation microstructures have a slightly higher monoclinic component, possibly correlated with predominant monoclinic simple shear flow in plastically-deformed samples. Overall, plots for CPO strength (ODF J-index), pole figure strength, CPO symmetry and seismic anisotropy show significant scattering. This could be related to sampling statistics, although our database is a factor of ten higher than the olivine database of 1998, or it could be related to the low symmetry (triclinic) structure of plagioclase resulting in the addition of degrees of freedom in the processes creating the CPOs.

Download Full-text

An analytic model for the stress-induced anisotropy of dry rocks

Geophysics ◽

10.1190/1.3567950 ◽

2011 ◽

Vol 76 (3) ◽

pp. WA125-WA133 ◽

Cited By ~ 41

Author(s):

Boris Gurevich ◽

Marina Pervukhina ◽

Dina Makarynska

Keyword(s):

Seismic Anisotropy ◽

Uniaxial Stress ◽

Azimuthal Anisotropy ◽

P Wave ◽

Analytic Model ◽

Induced Anisotropy ◽

S Wave ◽

Isotropic Rock ◽

Model Predictions ◽

Stress Induced Anisotropy

One of the main causes of azimuthal anisotropy in sedimentary rocks is anisotropy of tectonic stresses in the earth’s crust. We have developed an analytic model for seismic anisotropy caused by the application of a small anisotropic stress. We first considered an isotropic linearly elastic medium (porous or nonporous) permeated by a distribution of discontinuities with random (isotropic) orientation (such as randomly oriented compliant grain contacts or cracks). The geometry of individual discontinuities is not specified. Instead, their behavior is defined by a ratio B of the normal to tangential excess compliances. When this isotropic rock is subjected to a small compressive stress (isotropic or anisotropic), the number of cracks along a particular plane is reduced in proportion to the normal stress traction acting on that plane. This effect is modeled using the Sayers-Kachanov noninteractive approximation. The model predicts that such anisotropic crack closure yields elliptical anisotropy, regardless of the value of the compliance ratio B. It also predicts the ratio of Thomsen’s anisotropy parameters [Formula: see text] as a function of the compliance ratio B and Poisson’s ratio of the unstressed rock. A comparison of the model predictions with the results of laboratory measurements indicates a reasonable agreement for moderate magnitudes of uniaxial stress (as high as 30 MPa). These results can be used for differentiating stress-induced anisotropy from that caused by aligned fractures. Conversely, if the cause of anisotropy is known, then the anisotropy pattern allows one to estimate P-wave anisotropy from S-wave anisotropy.

Download Full-text

A database of plagioclase crystal preferred orientations (CPO) and microstructures – implications for CPO origin, strength, symmetry and seismic anisotropy in gabbroic rocks

Solid Earth ◽

10.5194/se-4-511-2013 ◽

2013 ◽

Vol 4 (2) ◽

pp. 511-542 ◽

Cited By ~ 26

Author(s):

T. Satsukawa ◽

B. Ildefonse ◽

D. Mainprice ◽

L. F. G. Morales ◽

K. Michibayashi ◽

...

Keyword(s):

Plastic Deformation ◽

Seismic Anisotropy ◽

P Wave ◽

S Wave ◽

Seismic Properties ◽

Strong Point ◽

Pole Figures ◽

Significant Difference ◽

Sampling Statistics ◽

P Type

Abstract. This study presents a unique database of 172 plagioclase Crystallographic Preferred Orientations (CPO) of variously deformed gabbroic rocks. The CPO characteristics as a function of the deformation regime (magmatic or crystal-plastic) are outlined and discussed. The studied samples are dominantly from slow- and fast-spread present-day ocean crust, as well as from the Oman ophiolite. Plagioclase is the dominant mineral phase in the studied samples. Plagioclase CPOs are grouped into three main categories: Axial-B, a strong point alignment of (010) with a girdle distribution of [100]; Axial-A, a strong point maximum concentration of [100] with parallel girdle distributions of (010) and (001); and P-type, point maxima of [100], (010), and (001). A majority of CPO patterns are Axial-B and P-type, in samples showing either magmatic or crystal-plastic deformation textures. Axial-A CPOs are less common; they represent 21% of the samples deformed by crystal-plastic flow. Although fabric strength (ODF J index) does not show any consistent variation as a function of the CPO patterns, there is a significant difference in the relationship between the ODF and pole figures J indices; the magmatic type microstructures have high (010) pole figures J indices, which increase linearly with ODF J index, whereas the high [100] pole figures J indices of plastically deformed samples vary in a more scattered manner with ODF J index. The multistage nature of plastic deformation superposed on a magmatic structure compared with magmatic flow, and the large number of possible slip-systems in plagioclase probably account for these differences. Calculated seismic properties (P wave and S wave velocities and anisotropies) of plagioclase aggregates show that anisotropy (up to 12% for P wave and 14% for S wave) tends to increase as a function of ODF J index. In comparison with the olivine 1998 CPO database, the magnitude of P wave anisotropy for a given J index is much less than olivine, whereas it is similar for S wave anisotropy. Despite a large variation of fabric patterns and geodynamic setting, seismic properties of plagioclase-rich rocks have similar magnitudes of anisotropy. There is a small difference in the aggregate elastic symmetry, with magmatic microstructures having higher orthorhombic and hexagonal components, whereas plastic deformation microstructures have a slightly higher monoclinic component, possibly correlated with predominant monoclinic simple shear flow in plastically deformed samples. Overall, plots for CPO strength (ODF J index), pole figure strength, CPO symmetry and seismic anisotropy show significant scattering. This could be related to sampling statistics, although our database is a factor of ten higher than the olivine database of 1998, or it could be related to the low symmetry (triclinic) structure of plagioclase resulting in the addition of degrees of freedom in the processes creating the CPOs.

Download Full-text

Fault parameters determined by near- and far-field data: The Wakasa Bay earthquake of March 26, 1963

Bulletin of the Seismological Society of America ◽

10.1785/bssa0640051369 ◽

1974 ◽

Vol 64 (5) ◽

pp. 1369-1382 ◽

Cited By ~ 1

Author(s):

Katsuyuki Abe

Keyword(s):

Effective Stress ◽

Stress Drop ◽

P Wave ◽

Slip Angle ◽

Polarization Angle ◽

S Wave ◽

Seismological Data ◽

Fault Parameters ◽

Wakasa Bay ◽

Dip Angle

Abstract The source process of the Wakasa Bay earthquake (M = 6.9, 35.80°N, 135.76°E, depth 4 km) which occurred near the west coast of Honshu Island, Japan, on March 26, 1963, is studied on the basis of the seismological data. Dynamic and static parameters of the faulting are determined by directly comparing synthetic seismograms with observed seismograms recorded at seismic near and far distances. The De Hoop-Haskell method is used for the synthesis. The average dislocation is determined to be 60 cm. The overall dislocation velocity is estimated to be 30 cm/sec, the rise time of the slip dislocation being determined as 2 sec. The other fault parameters determined, with supplementary data on the P-wave first motion, the S-wave polarization angle, and the aftershocks, are: source geometry, dip direction N 144°E, dip angle 68°, slip angle 22° (right-lateral strike-slip motion with some dip-slip component); fault dimension, 20 km length by 8 km width; rupture velocity, 2.3 km/sec (bilateral); seismic moment, 3.3 × 1025 dyne-cm; stress drop, 32 bars. The effective stress available to accelerate the fault motion is estimated to be about 40 bars. The approximate agreement between the effective stress and the stress drop suggests that most of the effective stress was released at the time of the earthquake.

Download Full-text

The Australian Plate and underlying mantle from waveform tomography with massive datasets

10.5194/egusphere-egu21-5272 ◽

2021 ◽

Author(s):

Janneke de Laat ◽

Sergei Lebedev ◽

Bruna Chagas de Melo ◽

Nicolas Celli ◽

Raffaele Bonadio

Keyword(s):

Structural Information ◽

Azimuthal Anisotropy ◽

P Wave ◽

Depth Range ◽

S Wave ◽

Crust And Upper Mantle ◽

Southern Boundary ◽

Velocity Anomaly ◽

Australian Plate ◽

Australian Continent

We present a new S-wave velocity tomographic model of the Australian Plate, Aus21.&#160; It is constrained by waveforms of 0.9 million seismograms with both the corresponding sources and stations located within the half-hemisphere centred at the Australian continent. Waveform inversion extracts structural information from surface, S- and multiple S-waves on the seismograms in the form of a set of linear equations. These equations are then combined in a large linear system and inverted jointly to obtain a tomographic model of S- and P-wave speeds and S-wave azimuthal anisotropy of the crust and upper mantle. The model has been validated by resolution tests and, for particular locations in Australia with notable differences with previous models, by independent inter-station measurements of surface-wave phase velocities, which we performed using available array data.&#160;&#160;Aus21 offers new insights into the structure and evolution of the Australian Plate and its boundaries. The Australian cratonic lithosphere occupies nearly all of the western and central Australia but shows substantial lateral heterogeneity. It extends up to the northern edge of the plate, where it is colliding with island arcs, without subducting. The rugged eastern boundary of the cratonic lithosphere provides a lithospheric definition of the Tasman Line. The thin, warm lithosphere below the eastern part of the continent, east of the Tasman Line, underlies the Cenozoic volcanism locations in the area. The lithosphere is also thin and warm below much of the Tasman Sea, underlying the Lord Howe hotspot and the submerged part of western Zealandia. A low velocity anomaly that may indicate the single source of the Lord Howe and Tasmanid hotspots is observed in the transition zone offshore the Australian continent, possibly also sourcing the East Australia hotspot. Another potential hotspot source is identified below the Kermadec Trench, causing an apparent slab gap in the overlying slab and possibly related to the Samoa Hotspot to the north. Below a portion of the South East Indian Ridge (the southern boundary of the Australian Plate) a pronounced high velocity anomaly is present in the 200-400 km depth range just east of the Australian-Antarctic Discordance (AAD), probably linked to the evolution of this chaotic ridge system.

Download Full-text

The 1998 Valhall microseismic data set: An integrated study of relocated sources, seismic multiplets, and S-wave splitting

Geophysics ◽

10.1190/1.3205028 ◽

2009 ◽

Vol 74 (5) ◽

pp. B183-B195 ◽

Cited By ~ 40

Author(s):

K. De Meersman ◽

J.-M. Kendall ◽

M. van der Baan

Keyword(s):

Seismic Anisotropy ◽

Amplitude Ratio ◽

Temporal Variations ◽

Oil Field ◽

P Wave ◽

Wave Form ◽

S Wave ◽

Data Set ◽

Wave Splitting ◽

Source Mechanisms

We relocate 303 microseismic events recorded in 1998 by sensors in a single borehole in the North Sea Valhall oil field. A semiautomated array analysis method repicks the P- and S-wave arrival times and P-wave polarizations, which are needed to locate these events. The relocated sources are confined predominantly to a [Formula: see text]-thick zone just above the reservoir, and location uncertainties are half those of previous efforts. Multiplet analysis identifies 40 multiplet groups, which include 208 of the 303 events. The largest group contains 24 events, and five groups contain 10 or more events. Within each multiplet group, we further improve arrival-time picking through crosscorrelation, which enhances the relative accuracy of the relocated events and reveals that more than 99% of the seismic activity lies spatially in three distinct clusters. The spatial distribution of events and wave-form similarities reveal two faultlike structures that match well with north-northwest–south-southeast-trending fault planes interpreted from 3D surface seismic data. Most waveform differences between multiplet groups located on these faults can be attributed to S-wave phase content and polarity or P-to-S amplitude ratio. The range in P-to-S amplitude ratios observed on the faults is explained best in terms of varying source mechanisms. We also find a correlation between multiplet groups and temporal variations in seismic anisotropy, as revealed by S-wave splitting analysis. We explain these findings in the context of a cyclic recharge and dissipation of cap-rock stresses in response to production-driven compaction of the underlying oil reservoir. The cyclic nature of this mechanism drives the short-term variations in seismic anisotropy and the reactivation of microseismic source mechanisms over time.

Download Full-text

Petrophysical properties variation of bitumen-bearing carbonates at increasing temperatures from laboratory to model

Geophysics ◽

10.1190/geo2019-0790.1 ◽

2020 ◽

Vol 85 (5) ◽

pp. MR297-MR308

Author(s):

Roberta Ruggieri ◽

Fabio Trippetta

Keyword(s):

Central Italy ◽

Carbonate Reservoir ◽

P Wave ◽

Petrophysical Properties ◽

S Wave ◽

Reservoir Properties ◽

Velocity Change ◽

Seismic Properties ◽

Wave Velocities ◽

Velocity Changes

Variations in reservoir seismic properties can be correlated to changes in saturated-fluid properties. Thus, the determination of variation in petrophysical properties of carbonate-bearing rocks is of interest to the oil exploration industry because unconventional oils, such as bitumen (HHC), are emerging as an alternative hydrocarbon reserve. We have investigated the temperature effects on laboratory seismic wave velocities of HHC-bearing carbonate rocks belonging to the Bolognano Formation (Majella Mountain, central Italy), which can be defined as a natural laboratory to study carbonate reservoir properties. We conduct an initial characterization in terms of porosity and density for the carbonate-bearing samples and then density and viscosity measurements for the residual HHC, extracted by HCl dissolution of the hosting rock. Acoustic wave velocities are recorded from ambient temperature to 90°C. Our acoustic velocity data point out an inverse relationship with temperature, and compressional (P) and shear (S) wave velocities show a distinct trend with increasing temperature depending on the amount of HHC content. Indeed, samples with the highest HHC content show a larger gradient of velocity changes in the temperature range of approximately 50°C–60°C, suggesting that the bitumen can be in a fluid state. Conversely, below approximately 50°C, the velocity gradient is lower because, at this temperature, bitumen can change its phase in a solid state. We also propose a theoretical model to predict the P-wave velocity change at different initial porosities for HHC-saturated samples suggesting that the velocity change mainly is related to the absolute volume of HHC.

Download Full-text

P-wave and S-wave azimuthal anisotropy at a naturally fractured gas reservoir, Bluebell‐Altamont Field, Utah

Geophysics ◽

10.1190/1.1444636 ◽

1999 ◽

Vol 64 (4) ◽

pp. 1312-1328 ◽

Cited By ~ 37

Author(s):

Heloise B. Lynn ◽

Wallace E. Beckham ◽

K. Michele Simon ◽

C. Richard Bates ◽

M. Layman ◽

...

Keyword(s):

Azimuthal Anisotropy ◽

P Wave ◽

Fracture Density ◽

Wave Reflections ◽

S Wave ◽

Gas Reservoir ◽

Green River ◽

Wave Data ◽

Reflection Data ◽

Naturally Fractured

Reflection P- and S-wave data were used in an investigation to determine the relative merits and strengths of these two data sets to characterize a naturally fractured gas reservoir in the Tertiary Upper Green River formation. The objective is to evaluate the viability of P-wave seismic to detect the presence of gas‐filled fractures, estimate fracture density and orientation, and compare the results with estimates obtained from the S-wave data. The P-wave response to vertical fractures must be evaluated at different source‐receiver azimuths (travelpaths) relative to fracture strike. Two perpendicular lines of multicomponent reflection data were acquired approximately parallel and normal to the dominant strike of Upper Green River fractures as obtained from outcrop, core analysis, and borehole image logs. The P-wave amplitude response is extracted from prestack amplitude variation with offset (AVO) analysis, which is compared to isotropic‐model AVO responses of gas sand versus brine sand in the Upper Green River. A nine‐component vertical seismic profile (VSP) was also obtained for calibration of S-wave reflections with P-wave reflections, and support of reflection S-wave results. The direction of the fast (S1) shear‐wave component from the reflection data and the VSP coincides with the northwest orientation of Upper Green River fractures, and the direction of maximum horizontal in‐situ stress as determined from borehole ellipticity logs. Significant differences were observed in the P-wave AVO gradient measured parallel and perpendicular to the orientation of Upper Green River fractures. Positive AVO gradients were associated with gas‐producing fractured intervals for propagation normal to fractures. AVO gradients measured normal to fractures at known waterwet zones were near zero or negative. A proportional relationship was observed between the azimuthal variation of the P-wave AVO gradient as measured at the tops of fractured intervals, and the fractional difference between the vertical traveltimes of split S-waves (the “S-wave anisotropy”) of the intervals.

Download Full-text