Anelastic degradation of acoustic pulses in rock

Measurements on acoustic pulses propagating in massive rock lead to a simple empirical relationship between the pulse rise time, τ and the time of propagation of a pulse, t:

$\text{τ=τ}{}_0\text{+C}\text{∫}\text{)}\text{T}\text{Q}{}^{\mathrm{?1}}\text{d}\text{t}$

where τ₀ is the initial rise time (at t = 0), Q is the anelastic parameter which may be expressed in terms of the fractional loss of energy per cycle of a sinusoidal wave, Q = 2π(ΔE/E)^?1, and is assumed to be essentially independent of frequency, and C is a constant whose value we estimate experimentally to be 0.53 ± 0.04. Of the linear theories of seismic pulse attenuation, model 2 of Azimi et al. (1968) is favoured. Pulse shapes computed from equations of Futterman (1962) also give C = 0.5, but the pulse arrives earlier than in a non-attenuating medium with the same elasticity and density. Pulse shapes calculated using Strick's (1967, 1970, 1971) theory give values of C incompatible with our results. The observations suggest that a method of estimating the Q-structure of the earth from seismic pulse rise times may have a particular advantage over the spectral ratio method.     

A note on pulse broadening and anelastic attenuation in near-surface rocks

A pulse rise-time method has been used to study pulse broadening on seismograms generated by a weight drop source at distances up to 600 m. Both source and receiver were placed on glacial overburden overlying a gneiss-monzonite rock body. One of two data sets showed a significant increase in pulse rise time, τ, as a function of travel time, T. This increase, if due to anelastic attenuation in the uppermost part of the rock body, implies a Q value of 243 ± 53, assuming a linear relationship between τ and T. The data were not capable of discriminating between the models of pulse broadening of Gladwin and Stacey (τ ∞ T) and Ricker ($\text{τ ∞ T}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$).     

The estimation of spectra from time‐frequency transforms for use in attenuation studies

We introduce the signal dependent time–frequency distribution, which is a time–frequency distribution that allows the user to optimize the tradeoff between joint time–frequency resolution and suppression of transform artefacts. The signal‐dependent time–frequency distribution, as well as the short‐time Fourier transform, Stockwell transform, and the Fourier transform are analysed for their ability to estimate the spectrum of a known wavelet used in a tuning wedge model. Next, the signal‐dependent time–frequency distribution, and fixed‐ and variable‐window transforms are used to estimate spectra from a zero‐offset synthetic seismogram. Attenuation is estimated from the associated spectral ratio curves, and the accuracy of the results is compared. The synthetic consisted of six pairs of strong reflections, based on real well‐log data, with a modeled intrinsic attenuation value of 1000/Q = 20. The signal‐dependent time–frequency distribution was the only time–frequency transform found to produce spectra that estimated consistent attenuation values, with an average of 1000/Q = 26±2; results from the fixed‐ and variable‐window transforms were 24±17 and 39±10, respectively. Finally, all three time–frequency transforms were used in a pre‐stack attenuation estimation method (the pre‐stack Q inversion algorithm) applied to a gather from a North Sea seismic dataset, to estimate attenuation between nine different strong reflections. In this case, the signal‐dependent time‐frequency distribution produced spectra more consistent with the constant‐Q model of attenuation assumed in the pre‐stack attenuation estimation algorithm: the average L1 residuals of the spectral ratio surfaces from the theoretical constant‐Q expectation for the signal‐dependent time‐frequency distribution, short‐time Fourier transform, and Stockwell transform were 0.12, 0.21, and 0.33, respectively. Based on the results shown, the signal‐dependent time‐frequency distribution is a time–frequency distribution that can provide more accurate and precise estimations of the amplitude spectrum of a reflection, due to a higher attainable time–frequency resolution.     

τ-p SEISMIC DATA FOR VISCOELASTIC MEDIA – PART 1: MODELLING1

Viscoelastic modelling reveals that the interaction of compressional-wave velocity C_p, compressional-wave quality factor Q_p, shear-wave velocity C_s, shear-wave quality factor Q_s and Poisson's ratio as a function of time intercept τ and ray parameter p, is complicated; however, distinct, potentially diagnostic behaviours are seen for different combinations of viscoelastic parameters. Synthetic seismograms for three viscoelastic reservoir models show that variations in the Poisson's ratio produce visible differences when compared to the corresponding elastic synthetic seismograms; these differences are attributable to interaction of the elastic parameters with Q_p and Q_s. When the P-wave acoustic impedance contrast is small, viscoelastic effects become more apparent and more useful for interpretation purposes. The corresponding amplitude and net phase spectra reveal significant differences between the elastic and the viscoelastic responses. When P-wave reflectivities are large, they tend to dominate the total response and to mask the Q reflectivity effects. The attenuation effects are manifested as an amplitude decay that increases with both time and ray parameter. The sensitivity of the computed seismic responses for various combinations of viscoelastic parameters suggests the opportunity for diagnostic interpretation of τ-p seismic data. The interpretation of the viscoelastic parameters can permit a better understanding of the rock types and pore fluid distribution existing in the subsurface.     

Intrinsic Qp seismic attenuation from the rise time of microearthquakes: a local scale application at Rio‐Antirrio,Western Greece

We determine the attenuation structure of a three‐dimensional medium based on first pulse‐width measurements from microearthquake data. Ninety‐five microearthquakes from a seventy stations local network were considered in this study. Measurements of the first half cycle of the wave, the so‐called rise time τ were carried out on high quality velocity seismograms and inverted to estimate the P‐waves intrinsic quality factor Q_p. The results of this investigation indicate that first pulse width data from a local microearthquake network permit retrieval with sufficient accuracy of the heterogeneous Q_p structure. The inferred attenuation variability corresponds to the known geological formations in the region.     

An attenuation study using earthquakes from the 1997 Umbria-Marche sequence

We studied spatial and temporal characteristics of seismic attenuation inCentral Italy using S- and coda- waves recorded by the MarchesanSeismograph Network from earthquakes located in the epicentral area ofthe 1997 Umbria-Marche sequence. The amplitude decay of the S waveswith distance was defined calculating empirical attenuation functions at 15frequencies between 1 and 25 Hz. We analyzed separately foreshocks andaftershocks and we found the same attenuation functions, suggesting thatthe possible temporal variations could be confined in a small area. Thefrequency dependence of Q _S was approximated by the equation Q _S=18 · f ^2.0between 1 and 10 Hz. At higher frequencies (10–25 Hz), the frequencydependence of Q _s weakens, having an average value of Q _S=990. We also estimated Q from coda waves (Q _C) using the single-scattering models of Aki andChouet (1975) and Sato (1977). We found that Q _C=77 · f ^0.6, (between 2 and 20Hz) at the western side of the mountain chain, using either foreshocks oraftershocks. This relation is consistent with previous estimates of Q _Creported for the Central Apennines. For a volume sampling the Colfioritobasin, the Apennines and the Marche region we found that Q _C=55 · f ^0.8,indicating highattenuation below the mountain belt. To detect small temporal changes ofQ, we calculated spectral ratios of 5 temporal doublets located in theepicentral area and recorded at the closest station. We found temporalchanges of Q that vary from 27% to 56%, depending on the locationof the doublets. This variability suggests that the temporal change ofattenuation may depend on the spatial variation of Q and perhaps on thespatial distribution of tectonic stress in the epicentral area.     

Attenuation of P and pP waves in Vrancea area – Romania

The seismic attenuation in the Vrancea region (Romania) is investigated from teleseismic recordings of P and pP waves during the four major, intermediate-depth Romanian events that occurred since the onset of digital instrumentation. Most stations are located in Canada and in the United States, being equipped with a variety of sensors, especially short-period ones. The amplitude spectral ratio method is used, assuming no frequency dependence of the Q _P factor in the range 0.2–2 Hz. No apparent correlation between the derived attenuation value and the type of recording sensor is observed. Lateral variations of the attenuation are obtained, with a very low Q _P area (values down to 33) located in the northwestern part of the Vrancea seismogenic volume. For the stations with different azimuth angles in relation to the epicentral area, Q _P values routinely exceed 200. Most likely, the low attenuation values are related to an upwelling mantle material located immediately beneath the crust, but limited in depth to at least 100 km.     

Lateral variations of coda Q and attenuation of seismic waves in Southwest Anatolia

The seismic quality factor (Q _c) and the attenuation coefficient (δ) in the earth’s crust in southwest (SW) Anatolia are estimated by using the coda wave method based on the decrease of coda wave amplitude by time on the seismogram. The quality factor Q _o, the value of Q _c at 1 Hz, and its frequency dependency η are determined from this method depending on the attenuation properties of scattered coda waves. δ is determined from the observations of amplitude variations of seismic waves. In applying the coda wave method, firstly, a type curve representing the average pattern of the individual coda decay curves for 0.75, 1.5, 3.0, 6.0, 12.0, and 24.0 Hz values was estimated. Secondly, lateral variation of coda Q and the attenuation coefficients for three main tectonic patterns are estimated. The shape of the type curve is controlled by the scattering and attenuation in the crustal volume sampled by the coda waves. The Q _o and η values vary from 30 to 180 and from 0.55 to 1.25, respectively for SW Anatolia. In SW Anatolia, coda Q–f relation is described by and δ = 0.008 km⁻¹. These results are expected to help in understanding the degree of tectonic complexity of the crust in SW Anatolia.     

A model for attenuation and scattering in the Earth's crust

The mechanisms contributing to the attenuation of earthquake ground motion in the distance range of 10 to 200 km are studied with the aid of laboratory data, coda wavesRg attenuation, strong motion attenuation measurements in the northeast United States and Canada, and theoretical models. The frequency range 1–10 Hz has been studied. The relative contributions to attenuation of anelasticity of crustal rocks (constantQ), fluid flow and scattering are evaluated. Scattering is found to be strong with an albedoB ₀=0.8–0.9 and a scattering extinction length of 17–32 km. The albedo is defined as the ratio of the total extinction length to the scattering extinction length. TheRg results indicate thatQ increases with depth in the upper kilometer or two of the crust, at least in New England. CodaQ appears to be equivalent to intrinsic (anelastic)Q and indicates that thisQ increases with frequency asQ=Q _o f ⁿ , wheren is in the range of 0.2–0.9. The intrinsic attenuation in the crust can be explained by a high constantQ (500Q _o2000) and a frequency dependent mechanism most likely due to fluid effects in rocks and cracks. A fluid-flow attenuation model gives a frequency dependence (QQ _o f ^0.5) similar to those determined from the analysis of coda waves of regional seismograms.Q is low near the surface and high in the body of the crust.     

Borehole study of compressional and shear attenuation of basalt flows penetrated by the Brugdan and William wells on the Faroes shelf

We investigated the seismic attenuation of compressional (P‐) and converted shear (S‐) waves through stacked basalt flows using short‐offset vertical seismic profile (VSP) recordings from the Brugdan (6104/21–1) and William (6005/13–1A) wells in the Faroe‐Shetland Trough. The seismic quality factors (Q) were evaluated with the classical spectral ratio method and a root‐mean‐square time‐domain amplitude technique. We found the latter method showed more robust results when analysing signals within the basalt sequence. For the Brugdan well we calculated effective Q estimates of 22–26 and 13–17 for P‐ and S‐waves, respectively, and 25–33 for P‐waves in the William well. An effective Q_S/Q_P ratio of 0.50–0.77 was found from a depth interval in the basalt flow sequence where we expect fully saturated rocks. P‐wave quality factor estimates are consistent with results from other VSP experiments in the North Atlantic Margin, while the S‐wave quality factor is one of the first estimates from a stacked basalt formation using VSP data. Synthetic modelling demonstrates that seismic attenuation for P‐ and S‐waves in the stacked basalt flow sequence is mainly caused by one‐dimensional scattering, while intrinsic absorption is small.     

Coda wave analysis for distinguishing attenuation due to isotropic scattering from attenuation due to absorption

When the quality factorQ is taken into account in attenuation studies, it is necessary to know the relative losses of wave energy due to scattering and to anelastic absorption. The coda is the most important phenomenon now known which is related to elastic scattering of seismic waves. Utilizing coda, this study presents relationships which give theQ factors of the medium around the recording station and discriminate between attenuations arising from elastic scattering (under the assumption of isotropic scattering) and those arising from anelastic absorption. This work proposes a technique for separately determining the attenuation due to isotropic scattering and that due to absorption from the observed envelope of coda waves.     

Body-wave Attenuation in the Region of Garda, Italy

We analyzed the spectral amplitude decay with hypocentral distance of P and S waves generated by 76 small magnitude earthquakes (M_L 0.9–3.8) located in the Garda region, Central-Eastern Alps, Italy. These events were recorded by 18 stations with velocity sensors, in a distance range between 8 and 120 km. We calculated nonparametric attenuation functions (NAF) and estimated the quality factor Q of both body waves at 17 different frequencies between 2 and 25 Hz. Assuming a homogeneous model we found that the Q frequency dependence of P and S can be approximated with the functions Q _P = 65 f ^0.9 and Q _S = 160 f ^0.6 , respectively. At 2 Hz the Q _S /Q _P ratio reaches the highest value of 2.8. At higher frequencies Q _S /Q _P varies between 0.7 and 1.7, suggesting that for this frequency band scattering may be an important attenuation mechanism in the region of Garda. To explore the variation of Q in depth, we estimated Q at short (r ≤ 30 km) and intermediate (35–90 km) distance paths. We found that in the shallow crust P waves attenuate more than S (1.3 < Q _S /Q _P < 2.5). Moreover, P waves traveling along paths in the lower crust (depths approximately greater than 30 km) attenuate more than S waves. To quantify the observed variability of Q in depth we considered a three-layer model and inverted the NAF to estimate Q in each layer. We found that in the crust Q increases with depth. However, in the upper mantle (~40–50 km depth) Q decreases and in particular the high frequency Q _S (f > 9 Hz) has values similar to those estimated for the shallow layer of the crust.     

Distribution of<Emphasis Type="Italic">L</Emphasis><Subscript>g</Subscript> coda<Emphasis Type="Italic">Q</Emphasis> in the Chinese continent and its adjacent regions

785 traces of vertical components from shallow earthquakes recorded by 10 CDSN (Chinese Digital Seismographic Network) stations and 5 GSN (Global Seismographic Network) stations were collected to study the attenuation characteristics ofL _g coda in the Chinese continent and its adjacent regions. The records were processed first using the stack spectral ratio method to obtain the average values ofQ ₀ (Q at 1Hz) and η, the frequency dependence, ofL _g coda in the ellipses corresponding to the paths. The back-projection technique was then employed to obtain the tomographic maps ofQ ₀ and η values, and the distribution of their errors. Results indicate that in the studied areaQ ₀ varies between 200 and 500. The lowest value ofQ ₀ exists in the Yun-nan-Tibetan region, while the highest value ofQ ₀ occurs in the southern edge of Siberian platform. η varies between 0.3 and 0.8. For most part of the studied area η varies inversely withQ ₀.     

Crustal attenuation for Jamaica, West Indies

The S and coda wave spectra of small earthquakes on the island of Jamaica were used to determine the near surface and coda Q attenuation, and Q _c,respectively. Q _c determined by the single-station method was found in the range of 1 to 10 Hertz to be given by the relationship, Q _c= 60 ± 5f ^0.87±0.05. This suggests that the Jamaican crust is highly attenuating which is further supported by the observation of rapid intensity fall-off with distance for earthquakes that have affected the island in the past. , determined from S-wave spectra with short travel times was found to be 0.058 ± 0.012 on the central crustal block, which makes up nearly two-thirds of the island, and 0.080 ± 0.014 in surrounding belt sub-regions. The pattern of values seems to fit with the surface geology in that the central block has areas of exposed outcrops of older and harder rock than the belts, which are characterized by thicker sedimentary sequences as well as intense fracturing and faulting.Atkinson and Boore (1998) and Atkinson(2001) presented an alternative method to stochastic modelling for ground motion in Eastern North America, whereby California attenuation relationships were modified to account for crustal differences invelocity-depth profile, Q and between both regions. Following their example, the California spectral attenuation relation of Boore, Joyner and Fumal (1997) was modified to account for differences between the California and Jamaica crust, resulting in an attenuation relation that is deemed to be more appropriate for Jamaica. Spectral accelerations for Jamaica when compared to California, are especially reduced beyond 20 km from the source and at high frequencies, f 1 hertz.The study concludes that the Jamaican crust, although having an oceanic composition is highly attenuating, which may be a result of intensive tectonic processes, whereas is consistent with near-rock conditions on the central block and soft rock conditions elsewhere on the island.     

Intrinsic and Scattering Seismic Attenuation in Southwestern Anatolia

The intrinsic dissipation and scattering attenuation in southwestern (SW) Anatolia, which is a tectonically active region, is studied using the coda waves. First the coda quality factor (Q_c) assuming single scattering is estimated from the slope of the coda-wave amplitude decay. Then the Multiple Lapse Time Window (MLTW) analysis is performed with a uniform earth model. Three non-overlapping temporal data windows are used to calculate the scattered seismic energy densities against the source-receiver distances, which, in turn, are used to calculate separate estimates of the intrinsic and scattering factors. In order to explore the frequency dependency, the observed seismograms are band pass-filtered at the center frequencies of 0.75, 1.5, 3.0, 6.0 and 12.0. The scattering attenuation (Q_s⁻¹) is found lower than the intrinsic attenuation (Q_i⁻¹) at all frequencies except at 0.75 Hz where the opposite is observed. Overall the intrinsic attenuation dominates over the scattering attenuation in the SW Anatolia region. The integrated energy curves obtained for the first energy window (i.e., 0–15 s) are somewhat irregular with distance while the second (i.e., 15–30 s) and third (i.e., 30–45 s) data windows exhibit more regular change with distance at most frequencies. The seismic albedo B₀ is determined as 0.61 at 0.75 Hz and 0.34 at 12.0 Hz while the total attenuation factor denoted by L_e⁻¹ changes in the range 0.034–0.017. For the source-station range 20–180 km considered the scattering attenuation is found strongly frequency dependent given by the power law Q_s⁻¹ = 0.010*f^−1.508. The same relations for Q_i⁻¹, Q_t⁻¹ (total), Q_c⁻¹ and (expected) hold as Q_i⁻¹ = 0.0090*f^−1.17, Q_t⁻¹ = 0.019*f^−1.31, Q_c⁻¹ = 0.008*f^−0.84 and respectively. Compared to the other attenuation factors Q_c⁻¹ and are less dependent on the frequency.     

FREQUENCY WAVENUMBER APPROACH OF THE τ-p TRANSFORM: SOME APPLICATIONS IN SEISMIC DATA PROCESSING*

The τ-p transform is an invertible transformation of seismic shot records expressed as a function of time and offset into the τ (intercept time) and p (ray parameter) domain. The τ-p transform is derived from the solution of the wave equation for a point source in a three-dimensional, vertically non-homogeneous medium and therefore is a true amplitude process for the assumed model. The main advantage of this transformation is to present a point source shot record as a series of plane wave experiments. The asymptotic expansion of this transformation is found to be useful in reflection seismic data processing. The τ-p and frequency-wavenumber (or f-k) processes are closely related. Indeed, the τ-p process embodies the frequency-wavenumber transformation, so the use of this technique suffers the same limitations as the f-k technique. In particular, the wavefield must be sampled with sufficient spatial density to avoid wavenumber aliasing. The computation of this transform and its inverse transform consists of a two-dimensional Fast Fourier Transform followed by an interpolation, then by an inverse-time Fast Fourier Transform. This technique is extended from a vertically inhomogeneous three-dimensional medium to a vertically and laterally inhomogeneous three-dimensional medium. The τ-p transform may create artifacts (truncation and aliasing effects) which can be reduced by a finer spatial density of geophone groups by a balancing of the seismic data and by a tapering of the extremities of the seismic data. The τ-p domain is used as a temporary domain where the attack of coherent noise is well addressed; this technique can be viewed as ‘time-variant f-k filtering’. In addition, the process of deconvolution and multiple suppression in the τ-p domain is at least as well addressed as in the time-offset domain.     

Attenuation of High Frequency P and S Waves in the Gujarat Region,India

The local earthquake waveforms recorded on broadband seismograph network of Institute of Seismological Research in Gujarat, India have been analyzed to understand the attenuation of high frequency (2–25 Hz) P and S waves in the region. The frequency dependent relationships for quality factors for P (Q _P) and S (Q _S) waves have been obtained using the spectral ratio method for three regions namely, Kachchh, Saurashtra and Mainland Gujarat. The earthquakes recorded at nine stations of Kachchh, five stations of Saurashtra and one station in mainland Gujarat have been used for this analysis. The estimated relations for average Q _P and Q _S are: Q _P = (105 ± 2) f ^{0.82 ± 0.01}, Q _S = (74 ± 2) f ^{1.06 ± 0.01} for Kachchh region; Q _P = (148 ± 2) f ^{0.92 ± 0.01}, Q _S = (149 ± 14) f ^{1.43 ± 0.05} for Saurashtra region and Q _P = (163 ± 7) f ^{0.77 ± 0.03}, Q _S = (118 ± 34) f ^{0.65 ± 0.14} for mainland Gujarat region. The low Q (<200) and high exponent of f (>0.5) as obtained from present analysis indicate the predominant seismic activities in the region. The lowest Q values obtained for the Kachchh region implies that the area is relatively more attenuative and heterogeneous than other two regions. A comparison between Q _S estimated in this study and coda Q (Qc) previously reported by others for Kachchh region shows that Q _C > Q _S for the frequency range of interest showing the enrichment of coda waves and the importance of scattering attenuation to the attenuation of S waves in the Kachchh region infested with faults and fractures. The Q _S/Q _P ratio is found to be less than 1 for Kachchh and Mainland Gujarat regions and close to unity for Saurashtra region. This reflects the difference in the geological composition of rocks in the regions. The frequency dependent relations developed in this study could be used for the estimation of earthquake source parameters as well as for simulating the strong earthquake ground motions in the region.     

Changes in source and site effects compared to codaQ −1 temporal variations using microearthquakes doublets in California

We analyse spectral ratio of the coda of doublets of microearthquakes. Our purpose is to find evidence for temporal changes of the attenuation in the crust before a large magnitude earthquake. A Moving Window Cross Spectral analysis of the coda of doublets gives a plot of the spectral ratio as a function of lapse time along the seismogram, for several frequency bands (SR(T, f) plot). From a certain pattern in theSR(T, f) plot, we should infer a temporal change in coda attenuation. Several doublets recorded in Central California by the USGS network are analysed.Using events very close in time from one another, we show that the radiation pattern can be different enough to induce important variations in the spectral ratio of the first arrivals and of the coda.Another doublet exhibits a strong variation of the low frequencies for stations in the region of Hollister (California), wherePhillips andAki (1986) have noted a strong amplification of low frequencies, that they attribute to site effects on unconsolidated sediments in the fault zone. These variations could be related to slight changes in local conditions (creep, or water table).On the other hand, some doublets, in the vicinity and close in time to the August 1979 Coyote Lake Earthquake (M=5.9), show no variation inSR(T, f) related to this earthquake: this proves that there was no major change in attenuation in the crust preceding this large shock. If a change occurred, it should have been confined to a very limited region, which was not sampled by the many paths we studied.The employed method probably provides today the most accurate estimation of spectral ratio in the coda of microearthquakes. It shows that there are numerous kinds of variations and that it is not straightforward to relate them to coda attenuation changes only. Also, we must be circumspect when dealing with the coda of microearthquakes in relation to forecasting earthquakes.     

Anelasticity of the crust and upper mantle beneath north and central India from the inversion of observed Love and Rayleigh wave attenuation data

The fundamental mode Love and Rayleigh waves generated by earthquakes occurring in Kashmir, Nepal Himalaya, northeast India and Burma and recorded at Hyderabad, New Delhi and Kodaikanal seismic stations are analysed. Love and Rayleigh wave attenuation coefficients are obtained at time periods of 15–100 seconds, using the spectral amplitude of these waves for 23 different paths along northern (across Burma to New Delhi) and central (across Kashmir, Nepal Himalaya and northeast India to Hyderabad and Kodaikanal) India. Love wave attenuation coefficients are found to vary from 0.0003 to 0.0022 km^–1 for northern India and 0.00003 km^–1 to 0.00016 km^–1 for central India. Similarly, Rayleigh wave attenuation coefficients vary from 0.0002 km^–1 to 0.0016 km^–1 for northern India and 0.00001 km^–1 to 0.0009 km^–1 for central India. Backus and Gilbert inversion theory is applied to these surface wave attenuation data to obtainQ ^–1 models for the crust and uppermost mantle beneath northern and central India. Inversion of Love and Rayleigh wave attenuation data shows a highly attenuating zone centred at a depth of 20–80 km with lowQ for northern India. Similarly, inversion of Love and Rayleigh wave attenuation data shows a high attenuation zone below a depth of 100 km. The inferred lowQ value at mid-crustal depth (high attenuating zone) in the model for northern India can be by underthrusting of the Indian plate beneath the Eurasian plate which has caused a low velocity zone at this shallow depth. The gradual increase ofQ ^–1 from shallow to deeper depth shows that the lithosphere-asthenosphere boundary is not sharply defined beneath central India, but rather it represents a gradual transformation, which starts beneath the uppermost mantle. The lithospheric thickness is 100 km beneath central India and below that the asthenosphere shows higher attenuation, a factor of about two greater than that in the lithosphere. The very lowQ can be explained by changes in the chemical constitution taking place in the uppermost mantle.     

Temporal correlations and clustering of landslides

This paper examines temporal correlations and temporal clustering of a proxy historical landslide time series, 2255 reported landslides 1951–2002, for an area in the Emilia‐Romagna Region, Italy. Landslide intensity is measured by the number of reported landslides in a day (D_L) and in an ‘event’ (S_event) of consecutive days with landsliding. The non‐zero values in both time series D_L and S_event are unequally spaced in time, and have heavy‐tailed frequency‐size distributions. To examine temporal correlations, we use power‐spectral analysis (Lomb periodogram) and surrogate data analysis, confronting our original D_L and S_event time series with 1000 shuffled (uncorrelated) versions. We conclude that the landslide intensity series D_L has strong temporal correlations and S_event has likely temporal correlations. To examine temporal clustering in D_L and S_event, we consider extremes over different landslide intensity thresholds. We first examine the statistical distribution of interextreme occurrence times, τ, and find Weibull distributions with parameter γ << 1·0 [D_L] and γ < 1·0 [S_event]; thus D_L and S_event each have temporal correlations, but S_event to a lesser degree. We next examine correlations between successive interextreme occurrence times, τ. Using autocorrelation analysis applied to τ, combined with surrogate data analysis, we find for D_L linear correlations in τ, but for S_event inconclusive results. However, using Kendall's rank correlation analysis we find for both D_L and S_event the series of τ are strongly correlated. Finally, we apply Fano Factor analysis, finding for both D_L and S_event the timings of extremes over a given threshold exhibit a fractal structure and are clustered in time. In this paper, we provide a framework for examining time series where the non‐zero values are strongly unequally spaced and heavy‐tailed, particularly important in the Earth Sciences due to their common occurrence, and find that landslide intensity time series exhibit temporal correlations and clustering. Many landslide models currently are designed under the assumption that landslides are uncorrelated in time, which we show is false. Copyright © 2010 John Wiley & Sons, Ltd.