A theoretical approach to the propagation of interacting cracks

We propose a scheme to compute interaction effects between two randomly oriented cracks under compressive stresses and we discuss the role crack interactions play in the crack coalescence process. Stress intensity factors are computed by using an iterative technique based on the method of successive approximations. Once crack propagation occurs, curved wing cracks grow from the initial crack tips. The stress intensity factors at the wing crack tips are calculated as the sum of two terms: a component for a single wing crack subjected to both the applied stresses and the interaction effect, and a component due to the sliding of the initial crack. We have applied our procedure to various crack geometries. Our results show that interaction effects act on the crack propagation path. For cracks under tension, our approach correctly predicts the curving, hook-shaped paths of interacting cracks that have been observed in various materials. For en echelon compressive cracks, interaction effects depend on the geometry of stepping. For right-stepping cracks, no mode I crack coalescence occurs. A mixedmode propagation criterion may be introduced to check whether coalescing secondary shear fractures initiate. For left-stepping cracks, depending on whether or not there is overlapping, crack coalescence is achieved by tension wing cracks at the inner crack tips. Without overlapping, the growing wing cracks delimit a region where a tensile secondary fracture may develop and lead to coalescence. These results are consistent with previous work and show that our procedure may be now extended to a population of cracks.     

Non-hypersingular boundary integral equations for two-dimensional non-planar crack analysis

We derive a set of non-hypersingular boundary integral equations, both elastodynamic and elastostatic, for the analysis of arbitrarily shaped 2-D anti-plane and in-plane cracks located in an infinite homogeneous isotropic medium, rendered in a unified nomenclature for all cases. The hypersingularities that appear in the usual formulations for the dynamic cases, existent both at the source point and at the wavefront, are removed by way of a regularization technique based on integration by parts. The equations for the in-plane cases are presented in terms of a local Cartesian coordinate system, one of the axes of which is always held locally tangential to the crack trace. The expressions for the elastic field at any point on the model plane are also given.
Our formulations are shown to yield accurate numerical results, as long as appropriate stabilization measures are taken in the numerical scheme. The numerical applicability of our method to non-planar crack problems is illustrated by simulations of dynamic growth of a hackly crack which has small off-plane side-branches. The results imply that the branching of a crack brings about a significant decrease in the crack-tip stress concentration level and consequently may play an essential role in the arrest of earthquake rupturing.     

Crack propagation by differential insolation on desert surface clasts

In the southwest U.S., cracks in alluvial fan surface clasts have a preferred orientation independent of rock fabric and shape. In this paper, we show that differential insolation of incipient cracks of random orientations predicts a distribution of crack orientations consistent with field observations. In this model, crack growth by hydration and/or thermal weathering is primarily a function of local water content at the crack tip. Crack tips that experience minimal solar insolation maintain a greater average moisture and, hence, weather more rapidly than cracks that experience greater solar insolation. To show this, we used a numerical radiative transfer code to quantify the solar insolation of rectangular cracks at 35° N. latitude with a range of depths and orientations. The amount of solar energy reaching the bottom of each crack was calculated at 5-min intervals over the day for several days of the year to determine hourly, daily, seasonal, and annual energy deposition as a function of crack depth and orientation. By assuming that only crack orientations that effectively shield their interiors and minimize their water loss are able to grow, the pattern of cracks produced by the model is consistent with field observations. The annual average insolation, which controls water retention, is associated with the two primary modes of crack orientation. The effect of daily recharge by summer rains of the North American monsoon system is consistent with the observed deviations from these primary modes. Model results suggest that both the annual average insolation and the daily pattern of rainfall is recorded in the preferred crack orientations of surface clasts in the southwest U.S.     

Evaluation of stress and displacement fields due to an elliptical plane shear crack

Summary. Following the classic work of Eshelby, the slip and stress distributions due to an elliptical plane shear crack are evaluated. The relation between average (or maximum) slip on the crack and the (constant) static stress drop, for faults of different aspect ratios, is found. The slip vector is not parallel to the applied stress but makes a small angle to it, except when the stress is applied along the major or minor axis of the ellipse. The stress -distribution around the crack shows that in addition to the expected stress concentration along the crack edge, there are broad regions of stress increase off the crack plane for circular and elliptical cracks, similar to those known to exist for in-plane but not for antiplane shear cracks. Whether the stress- intensity factor at the end of one axis is greater or less than that at the end of the other axis ( k_a ≶ k_b ), depends on the condition: √ b/a ≶ (1 − v ) where a and b are the semi-axes of the ellipse, k_a and k_b are the stress-intensity factors at the end of the a- and b -axes and v is Poisson's ratio. The total stress-intensity factor varies smoothly along the edge of the ellipse from one axis to the other and it is found that this variation is rather small.     

Wave speeds and attenuation of elastic waves in material containing cracks

Summary. Expressions now exist from which may be calculated the propagation constants of elastic waves travelling through material containing a distribution of cracks. The cracks are randomly distributed in position and may be randomly orientated. The wavelengths involved are assumed to be large compared with the size of the cracks and with their separation distances so that the formulae, based on the mean taken over a statistical ensemble, may reasonably be used to predict the properties of a single sample. The results are valid only for small concentrations of cracks.
Explicit expressions, correct to lowest order in the ratio of the crack size to a wavelength, are derived here for the overall elastic parameters and the overall wave speeds and attenuation of elastic waves in cracked materials where the mean crack is circular, and the cracks are either aligned or randomly orientated. The cracks may be empty or filled with solid or fluid material. These results are achieved on the basis of simply the static solution for an ellipsoidal inclusion under stress.
The extension to different distributions of orientation or to mixtures of different types of crack is quite straightforward.     

Simulation of the spontaneous growth of a dynamic crack without constraints on the crack tip path

The spontaneous growth of a dynamic in-plane shear crack is simulated using a newly developed method of analysis in which no a priori constraint is required for the crack tip path, unlike in other classical studies. We formulate the problem in terms of boundary integral equations; the hypersingularities of the integration kernels are removed by taking the finite parts. Our analysis shows that dynamic crack growth is spontaneously arrested soon after the bending of the crack tips, even in a uniformly stressed medium with homogeneously distributed fracture strengths. This shows that the dynamics of crack growth has a significant effect on forming the non-planar crack shape, and consequently plays an essential role in the arrest of earthquake rupturing.     

Acoustoelasticity in cracked solids

A new model that accounts for the stress dependence of the phase velocity of elastodynamic waves propagating in a cracked solid under compression is presented. The phase velocities of longitudinal and shear waves are derived from the effective elastic properties of a cracked solid, which are evaluated within the framework of Kachanov's approach. Following Kachanov, the extra-compliance tensor of the cracked solid is related to the crack compliances, which display a marked non-linear behaviour when subjected to a compressive load. Such non-linear behaviour is shown to be derived from the elastic interaction between the contacting crack faces under compression. This work does not address the effect of mutual interaction among cracks and the generation of higher harmonics due to the medium non-linearity. Numerical examples are presented that illustrate the phase velocity changes occurring in a solid with a random distribution of parallel cracks as a function of an external compressive load. A distinctive feature of the acoustoelastic effect in solids with large parallel fractures and in solids with distributions of aligned microcracks is also illustrated.     

Shallow earthquakes in a viscoelastic shear zone with depth-dependent friction and rheology

Summary. Most crustal earthquakes of the world are observed to occur within a seismogenic layer which extends from the Earth's surface to a depth of a few tens of kilometres at most. A model is proposed in which the shear zone along a transcurrent plate margin is represented as a viscoelastic medium with depth-dependent power-law rheology. A frictional resistance linearly increasing with depth is assumed on a vertical transcurrent fault within the shear zone. Such a model is able to reproduce a continuous transition from the brittle behaviour of the upper crust to the ductile behaviour at depth. Assuming that the shear zone is subjected to a constant strain rate from the opposite motions of the two adjacent plates, it is found that there exists a maximum depth H below which tectonic stress can never reach the frictional threshold: this may be identified as the maximum depth of earthquake nucleation. The value of H is consistent with observations for plausible values of the model parameters. The stress evolution in the shear zone is calculated in the linear approximation of the constitutive equation. A change in rigidity with depth, which is also introduced in the model, may reproduce the high vertical gradient of shear stress, which has been measured across the San Andreas fault, and the fact that most earthquakes are nucleated at some depth in the seismogenic layer. A crack which drops the ambient stress to the dynamic frictional level is then introduced in the model. To this aim, a crack solution is employed without a stress singularity at its edges, which is compatible with a frictional stress threshold criterion for fracture. A constraint on the vertical friction gradient is obtained if such cracks are assumed to be entirely confined within the seismogenic layer.     

Shear-wave polarizations near the North Anatolian Fault – III. Observations of temporal changes.

Summary. Almost all shear-waves from local earthquakes recorded on closely-spaced three-component seismometer networks deployed near the North Anatolian Fault, Turkey, in two experiments in 1979 and 1980, display shear-wave splitting. The observations are consistent with the presence of EDA (extensive-dilatancy anisotropy), distributions of fluid-filled cracks and microcracks aligned by the regional stress field. Temporal changes in the stress-field, which may occur before an earthquake, may modify the geometry and possibly the orientation of the EDA-microcracks, and lead to corresponding changes in the behaviour of the split shear-waves. A third experiment was undertaken in 1984 to investigate EDA further and to search for possible temporal variations of the polarization of the leading split shear-wave and the time delay between split shear-waves. Observations indicate that the polarization alignments, which are parallel to the strike of the parallel vertical EDA-cracks, are unaltered between 1979 and 1984, implying that the direction of the regional stress field has not changed significantly. Temporal changes in the stress field are more likely to cause changes in the crack density and/or aspect ratio, which would result in a corresponding change in time delay between the split shear-waves. We examine observations of time delay in relation to their propagation path with respect to the crack geometry since it is then possible to separate the effects of changes in crack density and changes in aspect ratio. With this procedure, a small temporal variation of time delays is found between 1979 and 1984, consistent with a decrease in crack density, and consequently a relaxation of stress, in this time period. No evidence was found for any observable variation of time delay over a six month observation span in 1984. We suggest that analysis of repeated shear-wave VSPs offers a technique for monitoring stress changes before earthquakes.     

Scattering of waves in elastic media containing passive and active cracks

Scattering of wavefields in a 3-D medium that includes passive and/or active structures, is numerically solved by using the boundary integral equation method (BIEM). The passive structures are velocity anomalies that generate scattered waves upon incidence, and the active structures contain endogenous fracture sources, which are dynamically triggered by the dynamic load due to the incident waves. Simple models are adopted to represent these structures: passive cracks act as scatterers and active cracks as fracture sources. We form cracks using circular boundaries, which consist of many boundary elements. Scattering of elastic waves by the boundaries of passive cracks is treated as an exterior problem in BIEM. In the case of active cracks, both the exterior and interior problems need to be solved, because elastic waves are generated by fracturing with stress drop, and the growing crack boundaries scatter the incident waves from the outside of the cracks. The passive cracks and/or active cracks are randomly distributed in an infinite homogeneous elastic medium. Calculations of the complete waveform considering a single scatter show that the active crack has weak influence on the attenuation of first arrivals but strong influence on the amplitudes of coda waves, as compared with those due to the passive crack. In the active structures, multiple scattering between cracks and the waves triggered by fracturing strongly affect the amplitudes of first arrivals and coda waves. Compared to the case of the passive structures, the attenuation of initial phase is weak and the coda amplitudes decrease slowly.     

Shear-wave polarization anisotropy in the Pacific Basin

Summary. Inversion of 295. Love- and Rayleigh-wave phase travel times across the Pacific Basin has yielded a structure which has a channel that is anisotropic with respect to the polarization of shear waves. The velocity of SH waves is approximately 4.24 km/s, and the velocity of SV waves is approximately 4.10 km/s in the low-velocity channel. The lid to the channel is isotropic with respect to the polarization of S waves and the velocity is approximately 4.60 km/s. The lid to the low-velocity channel increases in thickness with lithospheric age at the expense of the channel, and its thickness is apparently still increasing at a sea-floor age of 150 Myr.
These results can be explained in terms of a model with both randomly-and preferentially-oriented, liquid-filled cracks in the channel. In the model, it is assumed that the liquid-filled cracks are due to partial melting in the channel, and that any preferred orientation is caused by a shear-flow gradient resulting from differential motion between the lid and the deeper parts of the mantle.     

Hydraulic crack propagation in a porous medium

We develop a model for the propagation of a fluid-filled crack in a porous medium. The problem is motivated by the mechanism whereby drainage networks may form in partially molten rock below the Earth's lithosphere. Other applications include the propagation of hydraulic fractures in jointed rocks and in oil drilling operations, and the formation of dessication cracks in soils. Motivated by the lithospheric problem, we study a situation in which gravity acts in the direction of crack propagation. The model couples the elastoliydrodynamic equations of crack propagation with a pore pressure field in the porous rock, which drives the fluid flow which supplies the crack. The effect of the pore flow is to include a diffusional term in the evolution equation for the crack width, thus allowing a crack initiated at the base of the lithosphere to propagate down into the asthenosphere. Asymptotic and numerical solutions are presented for this crack evolution. However, the predicted drainage of melt into this crack is tiny compared with the upward percolative melt migration, and the predicted width of cracks (millimetres) is much too small to allow propagation of melt into the lithosphere without freezing. As a mechanism to explain magma fracturing in the lithosphere, the process described here therefore requires further refinement.     

Bilateral propagation of a spontaneous two-dimensional anti-plane shear crack under the influence of cohesion

Summary. We consider the problem of the bilateral extension of a two-dimensional anti-plane shear crack that initiates spontaneously at a point and extends under the influence of cohesive forces at the edges. An approximation to the stresses in the regions beyond the edges of the crack has been found that simplifies the calculation. The exact stresses in these regions are also found iteratively. In the cases of uniformly propagating cracks, the estimates of cohesive forces obtained from the approximation are close to the exact values for high crack speeds but are significantly different for low crack speeds. It is also found that if healing is initiated due to the encounter of one end of a uniformly propagating crack with an unbreakable barrier, the static stress drop in the torn region is constant but may either overshoot or undershoot the dynamical stress drop. In these cases, the final static slip distribution is obtained by freezing the dynamic solution along a characteristic line through the location of the barrier. We find that the crack length cannot be unambiguously derived from the far field spectral properties.     

Numerical simulation of the propagation of P waves in fractured media

We study the propagation of P waves through media containing open fractures by performing numerical simulations. The important parameter in such problems is the ratio between crack length and incident wavelength. When the wavelength of the incident wavefield is close to or shorter than the crack length, the scattered waves are efficiently excited and the attenuation of the primary waves can be observed on synthetic seismograms. On the other hand, when the incident wavelength is greater than the crack length, we can simulate the anisotropic behaviour of fractured media resulting from the scattering of seismic waves by the cracks through the time delay of the arrival of the transmitted wave. The method of calculation used is a boundary element method in which the Green's functions are computed by the discrete wavenumber method. For simplicity, the 2-D elastodynamic diffraction problem is considered. The rock matrix is supposed to be elastic, isotropic and homogeneous, while the cracks are all empty and have the same length and strike direction. An iterative method of calculation of the diffracted wavefield is developed in the case where a large number of cracks are present in order to reduce the computation time. The attenuation factor Q ⁻¹ of the direct waves passing through a fractured zone is measured in several frequency bands. We observe that the attenuation factor Q ⁻¹ of the direct P wave peaks around kd = 2, where k is the incident wavenumber and d the crack length, and decreases proportionally to ( kd ) ⁻¹ in the high-wavenumber range. In the long-wavelength domain, the velocity of the direct P wave measured for two different crack realizations is very close to the value predicted by Hudson's theory on the overall elastic properties of fractured materials.     

The Oblique Seismic Experiment on DSDP Leg 52

Summary. The Oblique Seismic Experiment (OSE) has been proposed to increase the usefulness of the IPOD crustal borehole as a means of investigating layer 2 of oceanic crust. Specific objectives are: to determine the lateral extent of the structure intersected by the borehole, to analyse the role of cracks in the velocity structure of layer 2, to look for anisotropy which may be caused by large cracks with a preferred orientation and, finally, to measure attenuation in oceanic crust.
The first successful Oblique Seismic Experiment in oceanic crust was carried out in 1977 March in a hole 400 miles north of Puerto Rico. An adequate study of lateral velocity variations was impossible because the hole was not deep enough, the hole was inadequately logged, and the small scale basement topography was not known. In general both P - and S -wave velocity profiles suggest that the crack density decreases with depth in layer 2. Velocities at the bottom of layer 2 are the same as matrix velocities for basalt, implying that crack density may be negligible at this depth. No convincing evidence for anisotropy in either layer 2 or 3 is found from travel time analysis. The hole was not deep enough to measure attenuation from normal incidence shots and amplitudes were not consistent enough to obtain a measure of attenuation from long range shots.     

Shear-wave polarizations near the North Anatolian Fault – II. Interpretation in terms of crack-induced anisotropy

Summary. The polarizations of shear waves recorded by networks of digital three-component seismometers immediately above small earthquakes near the North Anatolian Fault in Turkey display shear-wave splitting on almost all shear-wave seismograms recorded within the shear-wave window. This splitting is incompatible with source radiation-patterns propagating through simple isotropic structures but is compatible with effective anisotropy of the internal structure of the rock along the ray paths. This paper interprets the phenomena in terms of widespread crack-induced anisotropy. Distributions of stress-induced cracks model many features of the observations, and synthetic polarization diagrams calculated for propagation through simulated cracked rock are similar to the observed patterns. This evidence for widespread crack-induced anisotropy lends strong support to the hypothesis of extensive-dilatancy anisotropy (EDA) suggested by laboratory experiments in subcritical crack-growth. The crucial evidence confirming some form of EDA would be observations of temporal changes in shear-wave splitting as the stress field alters the crack density and crack geometry. There is some weak evidence for such temporal changes at one site, but further analysis of suitable digital three-component seismometer networks in seismic areas is required to confirm EDA.     

High-frequency radiation from crack (stress drop) models of earthquake faulting

We present a theory for the radiation of high-frequency waves by earthquake faults. We model the fault as a planar region in which the stress drops to the kinematic friction during slip. This model is entirely equivalent to a shear crack. For two-dimensional fault models we show that the high frequencies originate from the stress and slip velocity concentrations in the vicinity of the fault's edges. These stress concentrations radiate when the crack expands with accelerated motion. The most efficient generation of high-frequency waves occurs when the rupture velocity changes abruptly. In this case, the displacement spectrum has an ω^-2 behaviour at high frequencies. The excitation is proportional to the intensity of the stress concentration near the crack tips and to the change in the focusing factor due to rupture velocity. We extend these two-dimensional results to more general three-dimensional fault models in the case when the rupture velocity changes simultaneously on the rupture front. Results are similar to those described for two-dimensional faults. We apply the theory to the case of a circular fault that grows at constant velocity and stops suddenly. The present theory is in excellent agreement with a numerical solution of the same problem.
Our results provide upper bounds to the high-frequency radiation from more realistic models in which rupture velocity does not change suddenly. The ω^-2 is the minimum possible decay at high frequencies for any crack model of the source.     

Seismic radiation from dynamic coalescence, and the reconstruction of dynamic source parameters on a planar fault

The dynamic coalescence of two mode II cracks on a planar fault is simulated here using the elastodynamic boundary integral equation method. We focus on the complexity of the resultant slip rate and seismic radiation in the crack coalescence model (CCM) and on the reconstruction of a single crack model (SCM) that can reproduce the CCM waveforms from heterogeneous source parameters rather than coalescence. Simulation results reveal that localized higher slip rates are generated by coalescence as a result of stress interaction between the approaching crack tips. The synthesized seismic radiation exhibits a distinct coalescence phase that has striking similarities to stopping phases in the radiation and propagation properties. The corresponding SCM yields a singular increase in the stress drop distribution, which is accompanied by a sudden decrease in it across the point of coalescence in the CCM. This implies that the generation of high-frequency radiation is more efficient from coalescence than from stopping, although both phenomena exhibit the same strong  ω⁻²  -type displacement spectra.     

A higher order approximation to the wave propagation constants for a cracked solid

Summary . Expressions are available in the literature for the propagation constants of a 'mean' wave travelling in a material with a random distribution of cracks. These are approximations for small crack density and crack dimensions small compared with a wavelength. The formulae provided to second order in crack density by the method of smoothing are here extended to the case where the cracks consist of two or more sets, aligned in different directions.     

Mid-Jurassic through mid-Cretaceous extension in the Central Graben of the North Sea - part 2: estimates from faulting observed on a seismic reflection line

Abstract We present an interpretation of the structure and faulting of an industry multichannel line across the Central North Sea Graben. We observe substantial faulting between the mid-Jurassic and mid-Cretaceous and on the base Zechstein (salt) reflector. To estimate the extension from these faults we consider movement along both planar and curved faults. We demonstrate that summing the heave (the horizontal displacement) overestimates the time measure of elongation for planar, ‘domino-type’, faulting. However, for high-angle normal faults and up to 70% extension (β= 1.7) the heave only overestimates the extension by 13%. In the absence of other information, summing the heave provides a useful estimate of extension in the case of domino-type faulting. For curved ‘listric’ faults the heave is only a true measure of the elongation if the antithetic faulting which removes the voids is vertical. Antithetic movement along inclined shear planes implies significantly more extension. We used the two models; of faulting to introduce progressively greater amounts of internal deformation in the crustal rocks and sediments to attempt to reconcile the estimate of extension necessary to give the observed subsidence and that given by analysing the faults visible on the seismic line. Estimates of extension obtained by allowing antithetic faulting along inclined shear planes are consistent with the range of estimates necessary to account for the post-mid-Jurassic subsidence. The estimates for the prior mid-Jurassic faulting are still substantially less than those necessary to explain the subsidence. However, this is not of major concern as there are many reasons as to why analysis of the faulting should underestimate the pre mid-Jurassic extension. Our interpretation of the seismic line implies curved faults bottoming in the lithologically weak Zechstein salt. These faults are decoupled from the region below and, hence, do not reflect the geometry of the faulting in the basement.