Caustics in a periodically layered transversely isotropic medium with vertical symmetry axis

Wave propagation in a finely layered medium is a very important topic in seismic modelling and inversion. Here we analyse non‐vertical wave propagation in a periodically layered transversely isotropic (VTI) medium and show that the evanescent (attenuation) zones in the frequency‐horizontal slowness domain result in caustics in the group velocity domain. These caustics, which may appear for both the quasi‐compressional (qP) and quasi‐shear (qSV) wave surfaces are frequency dependent but display weak dependence at low frequencies. The caustics computed for a specific frequency differ from those observed at the low‐ and high‐frequency limits. We illustrate these caustics with a few numerical examples and snapshots computed for both qP‐ and qSV‐wave types.     

Approximate Relation Between the Ray Vector and the Wave Normal in Weakly Anisotropic Media

Determination of the ray vector (the unit vector specifying the direction of the group velocity vector) corresponding to a given wave normal (the unit vector parallel to the phase velocity vector or slowness vector) in an arbitrary anisotropic medium can be performed using the exact formula following from the ray tracing equations. The determination of the wave normal from the ray vector is, generally, a more complicated task, which is usually solved iteratively. We present a first-order perturbation formula for the approximate determination of the ray vector from a given wave normal and vice versa. The formula is applicable to qP as well as qS waves in directions, in which the waves can be dealt with separately (i.e. outside singular directions of qS waves). Performance of the approximate formulae is illustrated on models of transversely isotropic and orthorhombic symmetry. We show that the formula for the determination of the ray vector from the wave normal yields rather accurate results even for strong anisotropy. The formula for the determination of the wave normal from the ray vector works reasonably well in directions, in which the considered waves have convex slowness surfaces. Otherwise, it can yield, especially for stronger anisotropy, rather distorted results.     

SLANT STACKING AND ITS SIGNIFICANCE FOR ANISOTROPY*

Slant stacking transforms seismic data, recorded as a function of source-receiver offset and traveltime, into the domain of intercept time τ and ray parameter p. The shape of the τ-p-curves thus obtained is closely related to the slowness surfaces of the layers. A layer-stripping operation in the τ-p-domain removes all effects of the layers above the target layer. The resulting curve is equal to the slowness surface of the layer except for a scaling factor containing the thickness and dip of the layer. The slowness surface is a characteristic surface for anisotropic media. This makes the τ-p-domain very suitable for detecting and describing anisotropic layers. The relationship between the shape of τ-p-curves, the slowness surfaces, and the geometry of the layers is derived. Synthetic τ-p-curves calculated with the reflectivity method show some difficulties that can arise in determining the shape of the curves and in applying the stripping operation. It is shown that the effects of vertical inhomogeneity on the interpretation of τ-p-curves in terms of anisotropy are small.     

Inversion of Travel Times in Weakly Anisotropic Rock Samples

Based on the perturbation theory, inversion formulae for travel time of qP and qS waves in arbitrary weak anisotropy media are presented. The inversion formulae are linear expressions of elastic parameters expressed in terms of weak anisotropy (WA) parameters. The formulae of qS₁ and qS₂ waves have the same form and they can be used without identifying which wave is considered. A synthetic experiment similar to the measurement of rock sample in the laboratory is carried out to illustrate the efficiency of the presented inversion formulae. Two data sets for qP wave travel time from rock samples in the laboratory are inverted and 15 WA parameters are obtained.     

On occurrence of point singularities in orthorhombic media

Degeneracies of the slowness surfaces of shear (and compressional) waves in low-symmetry anisotropic media (such as orthorhombic), known as point singularities, pose difficulties during modelling and inversion, but can be potentially used in the latter as model parameter constraints. I analyse the quantity and spatial arrangement of point singularities in orthorhombic media, as well as their relation to the overall strength of velocity anisotropy. A classification scheme based on the number and spatial distribution of singularity directions is proposed. In normal orthorhombic models (where the principal shear moduli are smaller than the principal compressional moduli), point singularities can only be arranged in three distinct patterns, and media with the theoretical minimum (0) and maximum (16) number of singularities are not possible. In orthorhombic models resulting from embedding vertical fractures in transversely isotropic background, only two singularity distributions are possible, in contrast to what was previously thought. Although the total number of singularities is independent of the overall anisotropy strength, for general (non-normal) orthorhombic models, different spatial distributions of singularities become more probable with increasing magnitude of anisotropy.     

Locating Small Aftershocks Using a Small-Aperture Temporary Seismic Array

In this paper, we developed a specialized method to locate small aftershocks using a small-aperture temporary seismic array. The array location technique uses the first P arrival times to determine the horizontal slowness vector of the incoming P wave, then combines it with S–P times to determine the event location. In order to reduce the influence of lateral velocity variation on the location determinations, we generated slowness corrections using events well-located by the permanent broadband network as calibration events, then we applied the corrections to the estimated slownesses. Applications of slowness corrections significantly improved event locations. This method can be a useful tool to locate events recorded by temporary fault-zone arrays in the near field but unlocated by the regional permanent seismic network. As a test, we first applied this method to 64 well-located aftershocks of the 1992 Landers, California, earthquake, recorded by both the Caltech/USGS Southern California Seismic Network and a small-aperture, temporary seismic array. The average horizontal and vertical separations between our locations and the well-determined catalogue locations are 1.35 and 1.75 km, respectively. We then applied this method to 132 unlocated aftershocks recorded only by the temporary seismic array. The locations show a clear tendency to follow the surface traces of the mainshock rupture.     

Reflection and refraction of waves at the interface of an isotropic medium over a highly anisotropic medium

In this paper, we have considered the reflection and refraction of a plane wave at an interface between two half-spaces. The lower half-spaces is composed of highly anisotropic triclinic crystalline material and the upper half-space is homogeneous and isotropic. It has been assumed that due to incidence of a plane quasi-P (qP) wave, three types of waves, namely, quasi-P (qP), quasi-SV (qSV) and quasi-SH (qSH), will be generated in the lower half space whereas P and S waves will be generated in the upper half space. The phase velocities of all the quasi waves have been calculated. It has been assumed that the direction of particle motion is neither parallel nor perpendicular to the direction of propagation. Some specific relations have been established between directions of motion and propagation, respectively. The expressions for reflection coefficients of qP, qSV, qSH and refracted coefficients of P and SV waves are obtained. Results of reflection and refraction coefficients are presented.     

The effects of tillage induced surface roughness,slope and discharge rate on soil detachment by concentrated flow: An experimental study

Soil erosion in sloping cropland is a key water and soil conservation issue in the Loess Plateau region, China. How surface roughness influences soil detachment remains unclear due to the inconsistent results obtained from existing studies. The objectives of the present study were to evaluate the effects of tillage practices on soil detachment rate in sloping cropland and establish an accurate empirical model for the prediction of soil detachment rates. A series of movable bed experiments were conducted on sloping surfaces under three different tillage practices (manual dibbling, manual hoeing, and contour drilling), with a smooth surface (non-tillage) as a control. The research indicated that soil detachment rate significantly increased with roughness (p < 0.05) since the average soil detachment rate was the highest under the contour drilling treatment (6.762 g m⁻² s⁻¹), followed by manual hoeing (4.180 g m⁻² s⁻¹), and manual dibbling (3.334 g m⁻² s⁻¹); the lowest detachment rate was observed under the non-tillage treatment (3.214 g m⁻² s⁻¹). Slope gradient and unit discharge rate were positively correlated with soil detachment rate and proved to be more influential than soil surface roughness. Four composite hydraulic parameters were introduced to estimate soil detachment rate on tilled surfaces. Regression analyses revealed that stream power was the most effective predictor of soil detachment rate compared with unit length shear force, shear stress, and unit stream power. By integrating surface roughness as a variable, the detachment rate could be accurately described as a nonlinear function of stream power and surface roughness. The results of the present study indicate that tillage practice could influence soil loss on sloping cropland, considering the higher soil detachment rates under all tillage practices tested compared with non-tillage. The results are attributed mainly to concentrated flow caused by the high water storage levels on tilled surfaces, which could damage surface microtopography and, subsequently, the development of headcuts.     

Effect of varying normal stress on stability and dynamic motion of a spring-slider system with rate- and state-dependent friction

Incorporating rate and state friction laws, stability of linearly stable (i.e., with stiffness greater than the critical value) spring-slider systems subjected to triggering perturbations was analyzed under variable normal stress condition, and comparison was made between our results and that of fixed normal stress cases revealed in previous studies. For systems associated with the slip law, the critical magnitude of rate steps for triggering unstable slips are found to have a similar pattern to the fixed normal stress case, and the critical velocity steps scale with a/(b - a) when k = k_cr for both cases. The rate-step boundaries for the variable normal stress cases are revealed to be lower than the fixed normal stress case by 7 %–16 % for a relatively large α = 0.56 with (b - a)/a ranging from 0.25 to 1, indicating easier triggering under the variable normal stress condition with rate steps. The difference between fixed and variable normal stress cases decreases when the α value is smaller. In the same slip-law-type systems, critical displacements to trigger instability are revealed to be little affected by the variable normal stress condition. When k ≥ k_cr(V_*), a spring-slider system with the slowness law is much more stable than with the slip law, suggesting that the slowness law fits experimental data better when a single state variable is adopted. In stick-slip motions, the variable normal stress case has larger stress drops than the constant normal stress case. The variable normal stress has little effect on the range of slip velocity in systems associated with the slowness law, whereas systems associated with the slip law have a slowest slip velocity immensely smaller than the fixed normal stress case, by ~10 orders of magnitude.     

Inhomogeneous Harmonic Plane Waves in Viscoelastic Anisotropic Media

In viscoelastic media, the slowness vector p of plane waves is complex-valued, p = P + iA. The real-valued vectors P and A are usually called the propagation and the attenuation vector, repectively. For P and A nonparallel, the plane wave is called inhomogeneousThree basic approaches to the determination of the slowness vector of an inhomogeneous plane wave propagating in a homogeneous viscoelastic anisotropic medium are discussed. They differ in the specification of the mathematical form of the slowness vector p. We speak of directional specification, componental specification and mixed specification of the slowness vector. Individual specifications lead to the eigenvalue problems for 3 × 3 or 6 × 6 complex-valued matrices.In the directional specification of the slowness vector, the real-valued unit vectors N and M in the direction of P and A are assumed to be known. This has been the most common specification of the slowness vector used in the seismological literature. In the componental specification, the real-valued unit vectors N and M are not known in advance. Instead, the complex-valued vactorial component p of slowness vector p into an arbitrary plane with unit normal n is assumed to be known. Finally, the mixed specification is a special case of the componental specification with p purely imaginary. In the mixed specification, plane represents the plane of constant phase, so that N = ±n. Consequently, unit vector N is known, similarly as in the directional specification. Instead of unit vector M, however, the vectorial component d of the attenuation vector in the plane of constant phase is known.The simplest, most straightforward and transparent algorithms to determine the phase velocities and slowness vectors of inhomogeneous plane waves propagating in viscoelastic anisotropic media are obtained, if the mixed specification of the slowness vector is used. These algorithms are based on the solution of a conventional eigenvalue problem for 6 × 6 complex-valued matrices. The derived equations are quite general and universal. They can be used both for homogeneous and inhomogeneous plane waves, propagating in elastic or viscoelastic, isotropic or anisotropic media. Contrary to the mixed specififcation, the directional specification can hardly be used to determine the slowness vector of inhomogeneous plane waves propagating in viscoelastic anisotropic media. Although the procedure is based on 3 × 3 complex-valued matrices, it yields a cumbersome system of two coupled equations.     

Microtopographic evolution of mineral surfaces as a tool to identify and date young fault scarps in bedrock

Faulting that results in surface ruptures through bedrock can be particularly difficult to date. For example, stratigraphic control on the age of faulting, based on the age of the bedrock, often leaves unacceptably large uncertainty on the age of the faulting. From a paleoseismological perspective, there is a clear need to determine if a bedrock fault scarp is actually a young feature. For young fault ruptures that create fresh mineral surfaces, analysis of microtopography developed by weathering of the mineral surface may provide a quantifiable method for determining the fault age. The direct quantitative measurement of mineral surface microtopography using Atomic Force Microscopy affords a novel method to study the rupture ages of active faults. The method for using microtopographic evolution of mineral surfaces depends on three conditions. The first condition is that freshly exposed mineral cleavage surfaces, which can be described geometrically as planes, are formed during a rupture event. The formation of these fresh surfaces is analogous to the initiation of a weathering ‘clock’ that defines time t=0. Following cleavage formation dissolution of the planar mineral surface occurs. The rate of dissolution for a mineral species under given climatic conditions, governs the rate of mineral surface alteration. Thus as dissolution proceeds, the roughness of the mineral surface increases. We suggest that the progression of microtopographic roughness over time, which can be estimated by computing quantitative statistics derived from digital mineral surface topography, will systematically vary until a steady state surface topography is reached. The fractal dimension, D_f, is one such measure of surface roughness where, D_f at time t=0 is 2. The dissolution of the mineral surface increases the fractal dimension as the removal of material proceeds. We posit that somewhere between D_f=2 and D_f=3, the microtopography reaches a steady state. Therefore, in the pre-steady state stage of surface roughness, the quantitative measure of roughness of the mineral may serve as a measure of time elapsed since faulting. The period of time this initial stage of surface roughening represents is dependent on the mineral and as a consequence, its dissolution rate, in a specific set of environmental conditions. The time elapsed since fault rupture and grain cleavage can also be estimated from the measurement of the volume of material removed through dissolution. If part of the original cleavage surface remains and can be identified then AFM measurements of the surface microtopography can be used to calculate the dissolved volume per unit area.     

GUIDED LOW-FREQUENCY NOISE FROM AIRGUN SOURCES*

The presence of the water layer in marine seismic prospecting provides an effective waveguide for acoustic energy trapped between the sea-bed and the sea-surface. This energy persists to large ranges and can be the dominant early feature on far-offset traces. On airgun records, there is commonly a lower frequency set of arrivals following the water-trapped waves. These arrivals are not as obvious with higher frequency watergun sources. By using a combination of intercept-time/slowness (τ—p) mapping on observational data and theoretical modelling, we are able to identify the origin of the events. If a very rapid increase in a seismic wavespeed occurs beneath the sea-bed sediments, a new waveguide is formed bounded by the sea surface and this transition zone. The low frequency waves are principally guided within this thicker waveguide. Numerical filtering in the τ—p domain followed by trace reconstruction is very effective in removing the low frequency noise.     

THE INSENSITIVITY OF REFLECTED SH WAVES TO ANISOTROPY IN AN UNDERLYING LAYERED MEDIUM1

Propagation in the plane of mirror symmetry of a monoclinic medium, with displacement normal to the plane, is the most general circumstance in anisotropic media for which pure shear-wave propagation can occur at all angles. Because the pure shear mode is uncoupled from the other two modes, its slowness surface in the plane is an ellipse. When the mirror symmetry plane is vertical the pure shear waves in this plane are SH waves and the elliptical SH sheet of the slowness surface is, in general, tilted with respect to the vertical axis. Consider a half-space of such a monoclinic medium, called medium M, overlain by a half-space of isotropic medium I with plane SH waves incident on medium M propagating in the vertical symmetry plane of M. Contrary to the appearance of a lack of symmetry about the vertical axis due to the tilt of the SH-wave slowness ellipse, the reflection and transmission coefficients are symmetrical functions of the angle of incidence, and further, there exists an isotropic medium E with uniquely determined density and shear speed which gives exactly the same reflection and transmission coefficients underlying medium J as does monoclinic medium M. This means that the underlying monoclinic medium M can be replaced by isotropic medium E without changing the reflection and transmission coefficients for all values of the angle of incidence. Thus no set of SH seismic experiments performed in the isotropic medium in the symmetry plane of the underlying half-space can reveal anything about the monoclinic anisotropy of that underlying half-space. Moreover, even when the underlying monoclinic half-space is stratified, there exists a stratified isotropic half-space that gives the identical reflection coefficient as the stratified monoclinic half-space for all angles of incidence and all frequencies.     

The dynamic response of a fluid‐filled borehole under a normal point load on the free surface

The dynamic response of a semi‐infinite fluid‐filled borehole embedded in an elastic half‐space under a concentrated normal surface load is analysed in the long‐wavelength limit. The solution of the problem is obtained with integral transforms in the form of a double integral with respect to the slowness and frequency. The partial P‐ and SV‐wave responses are further transformed to path integrals along Cagniard paths in the complex slowness plane. Unlike the traditional Cagniard‐de Hoop technique based on the Laplace transform of time dependence, this paper is based on the Fourier transform. The tube‐wave response is presented as a causal integral over a slowness range. The resultant representation in the time‐domain is suitable for the numerical evaluation of the complete response in the fluid‐filled borehole, especially at large distances. Asymptotic analysis of seismic phases arising in the borehole is performed on the basis of the obtained solution. The complete asymptotic wavefield consists in P‐ and SV‐waves, the Rayleigh wave and the low‐frequency Stoneley (tube) wave. Pressure synthetics obtained by the use of the asymptotic formulas are shown to be in good agreement with straightforward calculations.     

Local Determination of Weak Anisotropy Parameters from qP-wave Slowness and Particle Motion Measurements

— We propose an algorithm for local evaluation of weak anisotropy (WA) parameters from measurements of slowness vector components and/or of particle motions of q P waves at individual receivers in a borehole in a multi-azimuthal multiple-source offset VSP experiment. As a byproduct the algorithm yields approximate angular variation of q P-wave phase velocity. The formulae are derived under assumption of weak but arbitrary anisotropy and lateral inhomogeneity of the medium. The algorithm is thus independent of structural complexities between the source and the receiver. If complete slowness vector is determinable from observed data, then the information about polarization can be used as an independent additional constraint. If only the component of the slowness along the borehole can be determined from observations (which is mostly the case), the inversion without information about polarization is impossible. We present several systems of equations which can be used when different numbers of components of the slowness vector are available. The SVD algorithm is used to solve an overdetermined system of linear equations for WA parameters for two test examples of synthetic multi-azimuthal multiple-source offset VSP data. The system of equations results from approximate first-order perturbation equations for the slowness and polarization vectors of the q P wave. Analysis of singular values and of variances of WA parameters is used for the estimation of chances to recover the sought parameters. Effects of varying number of profiles with sources and of noise added to “observed” data are illustrated. An important observation is that although, due to insufficient data, we often cannot recover all individual WA parameters with sufficient accuracy, angular phase velocity variation can be recovered rather well.     

A numerical study of comparison of two one-state-variable,rate- and state-dependent friction evolution laws

The two one-state-variable, rate- and state-dependent friction laws, i.e., the slip and slowness laws, are compared on the basis of dynamical behavior of a one-degree-of-freedom spring-slider model through numerical simulations. Results show that two (normalized) model parameters, i.e., Δ (the normalized characteristic slip distance) and β?α (the difference in two normalized parameters of friction laws), control the solutions. From given values of Δ, β, and α, for the slowness laws, the solution exists and the unique non-zero fixed point is stable when Δ>(β?α), yet not when Δ < (β?α). For the slip law, the solution exists for large ranges of model parameters and the number and stability of the non-zero fixed points change from one case to another. Results suggest that the slip law is more appropriate for controlling earthquake dynamics than the slowness law.     

Reflection/transmission laws for slowness vectors in viscoelastic anisotropic media

The reflection/transmission laws (R/T laws) of plane waves at a plane interface between two homogeneous anisotropic viscoelastic (dissipative) halfspaces are discussed. Algorithms for determining the slowness vectors of reflected/transmitted plane waves from the known slowness vector of the incident wave are proposed. In viscoelastic media, the slowness vectors of plane waves are complex-valued, p = P + iA, where P is the propagation vector, and A the attenuation vector. The proposed algorithms may be applied to bulk plane waves (A = 0), homogeneous plane waves (A ≠ 0, P and A parallel), and inhomogeneous plane waves (A ≠ 0, P and A non-parallel). The manner, in which the slowness vector is specified, plays an important role in the algorithms. For unrestricted anisotropy and viscoelasticity, the algorithms require an algebraic equation of the sixth degree to be solved in each halfspace. The degree of the algebraic equation decreases to four or two for simpler cases (isotropic media, plane waves in symmetry planes of anisotropic media). The physical consequences of the proposed algorithms are discussed in detail. vcerveny@seis.karlov.mff.cuni.cz     

Occurrence of water ponding on soil surfaces depending on infiltration rates on Mongolian rangeland

The occurrence of water ponding on soil surfaces during and after heavy rainfall produces surface run‐off or surface water accumulation in low‐lying areas, which might reduce the water supply to soils and result in a reduction of the soil water that plants can use, especially in arid climates. On Mongolian rangeland, we observed ponded water on the surface of a specific soil condition subjected to a heavy rainfall of 30 mm/hr. By contrast, ponded water was not observed for the same type of soil where livestock grazing had been removed for 6–8 years via a fence or for nearby soil containing less clay. We measured the infiltration rate (the saturated hydraulic conductivity of the surface soil, K_s) of the three sites by applying ponded water on the soil surface (an intake rate test). The results showed that K_s in the rangeland was lower than the rainfall intensity in the site where water ponded on the soil surface; however, K_s of the soil inside of the fence has recovered to 3 times that of the soil outside of the fence to exceed the rainfall intensity. Heavy rainfall that exceeds the infiltration rate occurs several times a year at the livestock grazing site where we observed ponded water. Slight water repellency of the soil reduces rain infiltration to increase the possibility of surface ponding for the soil.     

Breakdown of the slow manifold in the Shallow-Water equations

Abstract

Numerical solutions are obtained by implicit multigrid solvers for initial-value problems in the rotating Shallow-Water Equations (SWE) with spatially complex initial conditions. Companion solutions are also obtained with the Shallow-Water Balance Equations (SWBE), both to determine the initial conditions for the SWE and to provide a comparison solution that lies entirely on the slow, advective manifold. We make use of a control parameter (here the Rossby number, R) to regulate the degree of slowness and balance. While there are measurable discrepancies between the evolving SWE and SWBE solutions for all R, there is a distinct, spatially local breakdown both of the slow manifold in the SWE solution and in the closeness of correspondence between the SWE and SWBE solutions. This critical value for breakdown is only slightly smaller than the R values at which, first, the SWE evolution becomes singular (i.e., the fluid depth vanishes), or second, a consistent initial condition for the SWBE cannot be defined. This breakdown is most clearly evident in a sudden increase in vertical velocity near the center of a strong, cyclonic vortex; its behavior is primarily associated with an enhanced dissipation rather than an initiation of gravity-wave propagation. The numerical performance of the multigrid solvers is satisfactory even in the difficult circumstances near solution breakdown or singularity.     

Concentrated flow erosion rates reduced through biological geotextiles

Soil erosion by concentrated flow can cause serious environmental damage. Erosion‐control geotextiles have considerable potential for reducing concentrated flow erosion. However, limited data are available on the erosion‐reducing potential of geotextiles. In this study, the effectiveness of three biological geotextiles in reducing soil losses during concentrated flow is investigated. Hereto, runoff was simulated in a concentrated flow flume, filled with an erodible sandy loam on three slope gradients (13·5, 27·0 and 41·5%). Treatments included three biological geotextiles (borassus, buriti and bamboo) and one bare soil surface. Darcy–Weisbach friction coefficients ranged from 0·01 to 2·84. The highest values are observed for borassus covered soil surfaces, followed by buriti, bamboo and bare soil, respectively. The friction coefficients are linearly correlated with geotextile thickness. For the specific experimental conditions of this study, borassus geotextiles reduced soil detachment rate on average to 56%, buriti geotextiles to 59% and bamboo geotextiles to 66% of the soil detachment rate for bare soil surfaces. Total flow shear stress was the hydraulic parameter best predicting soil detachment rate for bare and geotextile covered surfaces (R² = 0·75–0·84, p <0·001, n = 12–15). The highest resistance against soil detachment was observed for the borassus covered soil surfaces, followed by buriti, bamboo and bare soil surfaces, respectively. Overall, biological geotextiles are less effective in controlling concentrated flow erosion compared with interrill erosion. Copyright © 2009 John Wiley & Sons, Ltd.