Anisotropic P-wave attenuation measured from a multi-azimuth surface seismic reflection survey

A system of aligned vertical fractures produces azimuthal variations in stacking velocity and amplitude variation with offset, characteristics often reported in seismic reflection data for hydrocarbon exploration. Studies of associated attenuation anisotropy have been mostly theoretical, laboratory or vertical seismic profiling based. We used an 11 common‐midpoint‐long portion of each of four marine surface‐seismic reflection profiles, intersecting each other at 45° within circa 100 m of a common location, to measure the azimuthal variation of effective attenuation, Q⁻¹_eff and stacking velocity, in a shallow interval, about 100 m thick, in which consistently orientated vertical fracturing was expected due to an underlying salt diapirism. We found qualitative and quantitative consistency between the azimuthal variation in the attenuation and stacking velocity, and published amplitude variation with offset results. The 135° azimuth line showed the least apparent attenuation (1000 Q⁻¹_eff= 16 ± 7) and the fastest stacking velocity, hence we infer it to be closest to the fracture trend: the orthogonal 45° line showed the most apparent attenuation (1000Q⁻¹_eff= 52 ± 15) and slowest stacking velocity. The variation of Q⁻¹_eff with azimuth φ is well fitted by 1000Q⁻¹_eff = 34 − 18cos[2(φ+40°)] giving a fracture direction of 140 ± 23° (±1SD, derived from ‘bootstrapping’ fits to all 11⁴ combinations of individual common‐midpoint/azimuth measurements), compared to 134 ± 47° from published amplitude variation with offset data. The effects of short‐window spectral estimation and choices of spectral ratio bandwidth and offset ranges used in attenuation analysis, individually give uncertainties of up to ±13° in fracture direction. This magnitude of azimuthal variation can be produced by credible crack geometries (e.g., dry cracks, radius 6.5 m, aspect ratio 3 × 10⁻⁵, crack density 0.2) but we do not claim these to be the actual properties of the interval studied, because of the lack of well control (and its consequences for the choice of theoretical model and host rock physical properties) and the small number of azimuths available here.     

Coherence methods in mapping AVO anomalies and predicting P-wave and S-wave impedances

Filters for migrated offset substacks are designed by partial coherence analysis to predict ‘normal’ amplitude variation with offset (AVO) in an anomaly free area. The same prediction filters generate localized prediction errors when applied in an AVO‐anomalous interval. These prediction errors are quantitatively related to the AVO gradient anomalies in a background that is related to the minimum AVO anomaly detectable from the data. The prediction‐error section is thus used to define a reliability threshold for the identification of AVO anomalies. Coherence analysis also enables quality control of AVO analysis and inversion. For example, predictions that are non‐localized and/or do not show structural conformity may indicate spatial variations in amplitude–offset scaling, seismic wavelet or signal‐to‐noise (S/N) ratio content. Scaling and waveform variations can be identified from inspection of the prediction filters and their frequency responses. S/N ratios can be estimated via multiple coherence analysis. AVO inversion of seismic data is unstable if not constrained. However, the use of a constraint on the estimated parameters has the undesirable effect of introducing biases into the inverted results: an additional bias‐correction step is then needed to retrieve unbiased results. An alternative form of AVO inversion that avoids additional corrections is proposed. This inversion is also fast as it inverts only AVO anomalies. A spectral coherence matching technique is employed to transform a zero‐offset extrapolation or near‐offset substack into P‐wave impedance. The same technique is applied to the prediction‐error section obtained by means of partial coherence, in order to estimate S‐wave velocity to P‐wave velocity (V_S/V_P) ratios. Both techniques assume that accurate well ties, reliable density measurements and P‐wave and S‐wave velocity logs are available, and that impedance contrasts are not too strong. A full Zoeppritz inversion is required when impedance contrasts that are too high are encountered. An added assumption is made for the inversion to the V_S/V_P ratio, i.e. the Gassmann fluid‐substitution theory is valid within the reservoir area. One synthetic example and one real North Sea in‐line survey illustrate the application of the two coherence methods.     

A laboratory study of seismic velocity and attenuation anisotropy in near-surface sedimentary rocks

The laboratory ultrasonic pulse‐echo method was used to collect accurate P‐ and S‐wave velocity (±0.3%) and attenuation (±10%) data at differential pressures of 5–50 MPa on water‐saturated core samples of sandstone, limestone and siltstone that were cut parallel and perpendicular to the vertical borehole axis. The results, when expressed in terms of the P‐ and S‐wave velocity and attenuation anisotropy parameters for weakly transversely isotropic media (ɛ, γ, ɛ_Q, γ_Q) show complex variations with pressure and lithology. In general, attenuation anisotropy is stronger and more sensitive to pressure changes than velocity anisotropy, regardless of lithology. Anisotropy is greatest (over 20% for velocity, over 70% for attenuation) in rocks with visible clay/organic matter laminations in hand specimens. Pressure sensitivities are attributed to the opening of microcracks with decreasing pressure. Changes in magnitude of velocity and attenuation anisotropy with effective pressure show similar trends, although they can show different signs (positive or negative values of ɛ, ɛ_Q, γ, γ_Q). We conclude that attenuation anisotropy in particular could prove useful to seismic monitoring of reservoir pressure changes if frequency‐dependent effects can be quantified and modelled.     

Source/receiver array directivity effects on marine seismic attenuation measurements

Attenuation of seismic waves, quantified by the seismic quality factor Q, holds important information for seismic interpretation, due to its sensitivity to rock and fluid properties. A recently published study of Q, based on surface seismic reflection data, used a modified spectral ratio approach (QVO), but both source and receiver responses were treated as isotropic, based on simple raypath arguments. Here, this assumption has been tested by computing apparent attenuation generated by frequency-dependent directivity of typical marine source and receiver arrays and acquisition geometries. Synthetic wavelet spectra were computed for reflected rays, summed over the first Fresnel zone, from the base of a single interval, 50–3000 m thick and velocity 2000 m/s, overlying a 2200 m/s half-space, and for offsets of 71–2071 m. The source and receiver geometry were those of an actual survey. The modelled spectra are clearly affected by directivity, most strongly because of surface ghosts. In general, the strong high-frequency component, produced by the array design, leads to apparently negative attenuation in individual reflection events, though this is dependent on offset and target depth. For shallow targets (less than 400–500 ms two-way traveltime (TWT) depth), apparent Q-values as extreme as ?50 to ?100 were obtained. For deeper target depths, the directivity effect is far smaller. The implications of the model study were tested on real data. QVO was applied to 20 true-spectrum-processed CMPs, in a shallow (405–730 ms TWT) and a deeper (1000–1300 ms TWT) interval, firstly using a measured far-field source signature (effectively isotropic), and secondly using computed directivity effects instead. Mean interval Q^?1-values for the deeper interval, 0.029 ± 0.011 and 0.027 ± 0.018 for conventional and directional processing, respectively, suggested no directivity influence on attenuation estimation. For the shallow interval (despite poor spectral signal-to-noise ratios and hence scattered attenuation estimates), directional processing removed directivity-generated irregularities from the spectral ratios, resulting in an improvement from Q^?1_int = ?0.036 ± 0.130 to a realistic Q^?1_int = 0.012 ± 0.030: different at 94% confidence level. Equivalent Q-values are: for the deeper interval, 35 and 37 for conventional and directional processing, respectively, and ?28 and 86 for the shallow interval. These results support the conclusions of the model studies, i.e. that source/receiver directivity has a negligible effect except for shallow targets (e.g. TWT depth ≤ 500 ms) imaged with conventional acquisition geometry. In such cases directivity corrections to spectra are strongly recommended.     

Experimental study on transient and steady‐state dynamics of bedforms in supply limited configuration

The purpose of the present study is to investigate experimentally the development of bedforms in a configuration where the sediment supply is limited. The experimental setup is a rectangular closed duct combining an innovative system to control the rate of sediment supply Q_in , and a digitizing system to measure in real time the 3D bedform topography. We carried out different sets of experiments with two sediment sizes (100 µm and 500 µm) varying both the sediment supply and the water flow rate to obtain a total of 46 different configurations. After a transient phase, steady sub‐centimeter bedforms of various shapes have been observed: barchans dunes, straight transverse dune, linguoid transverse dunes and bedload sheets. Height, spacing, migration speed, and mean bed elevation of the equilibrium bedforms were measured. For a given flow rate, two regimes were identified with fine sediment: (i) a monotonic increasing regime where the equilibrium bedform height and velocity increase with the sediment supply rate Q_in and (ii) an invariant regime for which both parameters are almost independent of Q_in. For coarse sediment, only the first regime is observed. We interpret the saturation of height and velocity for fine sediment bedforms as the transition from a supply‐limited regime to a transport‐limited regime in which the bedload flux has reached its maximum value under the prevailing flow conditions. We also demonstrate that all experiments can be rescaled if the migration speed and height of the bedforms are, respectively, divided and multiplied by the cube of the shear velocity. This normalization is independent of grain size and of bedform morphology. These experimental results provide a new quantification of the factors controlling equilibrium height and migration speed of bedforms in supply‐limited conditions against which theoretical and numerical models can be tested.     

A new laboratory technique for determining the compressional wave properties of marine sediments at sonic frequencies and in situ pressures

We describe a new laboratory technique for measuring the compressional wave velocity and attenuation of jacketed samples of unconsolidated marine sediments within the acoustic (sonic) frequency range 1–10 kHz and at elevated differential (confining – pore) pressures up to 2.413 MPa (350 psi). The method is particularly well suited to attenuation studies because the large sample length (up to 0.6 m long, diameter 0.069 m) is equivalent to about one wavelength, thus giving representative bulk values for heterogeneous samples. Placing a sediment sample in a water‐filled, thick‐walled, stainless steel Pulse Tube causes the spectrum of a broadband acoustic pulse to be modified into a decaying series of maxima and minima, from which the Stoneley and compressional wave, velocity and attenuation of the sample can be determined. Experiments show that PVC and copper jackets have a negligible effect on the measured values of sediment velocity and attenuation, which are accurate to better than ± 1.5% for velocity and up to ± 5% for attenuation. Pulse Tube velocity and attenuation values for sand and silty‐clay samples agree well with published data for similar sediments, adjusted for pressure, temperature, salinity and frequency using standard equations. Attenuation in sand decreases with pressure to small values below Q^?1 = 0.01 (Q greater than 100) for differential pressures over 1.5 MPa, equivalent to sub‐seafloor depths of about 150 m. By contrast, attenuation in silty clay shows little pressure dependence and intermediate Q^?1 values between 0.0206–0.0235 (Q = 49–43). The attenuation results fill a notable gap in the grain size range of published data sets. Overall, we show that the Pulse Tube method gives reliable acoustic velocity and attenuation results for typical marine sediments.     

Sediment Thicknesses and Q<Subscript>s</Subscript> vs. Q<Subscript>p</Subscript> Relations in the Kachchh Rift Basin,Gujarat, India Using Sp Converted Phases

Delineation of the top sedimentary structure and its Q_s vs. Q_p relationship using the travel-time difference of direct S and converted Sp phase is key to understanding the seismic hazard of any sedimentary basin area. We constructed filtered displacement waveforms from local ETNA Episensor acceleration recordings as well as local velocity recordings of aftershocks of the 2001 Bhuj earthquake recorded by the Kachchh seismological network of the National Geophysical Research Institute (NGRI), Hyderabad, India during 2001–2004. Stations are within 15–70km of epicenters, and the resulting displacement waveforms are generally simple, displaying prominent P, Sp, and S wave pulses. Particle motion of P and S waves suggest near-vertical raypaths consistent with preliminary depth estimates. The direct S wave on the horizontal component is characterized by lower frequency content than the converted Sp phase on the vertical component. This difference in frequency content between S and Sp phases can be explained in terms of different attenuation effects for P and S waves in the unconsolidated sediments. The Sp phase is generated by S-to-P phase conversion at the base of Mesozoic sediments of the Kachchh basin. Travel-time inversion (VELEST) of 2565 P and 2380 S arrivals from 658 well located aftershocks recorded at 8–14 three-component local seismic stations led to 1 D velocity models indicated very slow sediments in the upper 0–2 km depth range (Vp: 2.92 km/s and Vs: 0.90 km/s) and an increasing trend of velocities with depth at 2–40 km depth. The estimated sediment thicknesses beneath 12 accelerograph and 6 seismograph sites from the estimated velocity model and the travel-time difference between S and converted Sp phases reaches a maximum of (1.534 ± 0.117) km beneath Bandri (near the location of 2001 Bhuj mainshock) and attains a minimum sediment thickness of (0.858 ± 0.104) km beneath Ramvav and Burudia. The spectral ratios between Sp and S from 159 three-component accelerograms have been used to study seismic wave attenuation in the Kachchh rift basin. The estimated Q_s vs. Q_p relations for 12 accelerograph sites vary from Q_s = 0.184 Q_p (at Chobari) to Q_s = 0.505 Q_p (at Dudhai). For stations Chobari, Chopdwa, Jahawarnagar, Vondh and Tapar, the spectral ratio slopes and hence the calculated Q_s vs. Q_p relations are effectively the same, and the correlation coefficients are quite high (0.91–0.93). Stations Adhoi, Manfara, New Dudhai, Dudhai and Sikara have similar Q_s vs. Q_p relationships to each other and also have high correlation coefficients (0.78–0.87). The spectral ratios for stations Anjar and Ramvav are small and poorly constrained, resulting in less reliable Q_s vs. Q_p relations. This could be due to noisy data, fewer available waveforms, or scattering due to velocity heterogeneities and/or interface irregularities.     

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.     

Attenuation of P- and S-waves in the Chamoli Region, Himalaya, India

The attenuation properties of the crust in the Chamoli region of Himalaya have been examined by estimating the frequency-dependent relationships of quality factors for P waves (Q_α) and for S waves (Q_β) in the frequency range 1.5–24 Hz. The extended coda normalization method has been applied on the waveforms of 25 aftershocks of the 1999 Chamoli earthquake (M 6.4) recorded at five stations. The average value of Q_α is found to be varied from 68 at 1.5 Hz to 588 at 24 Hz while it varies from 126 at 1.5 Hz to 868 at 24 Hz for Q_β. The estimated frequency-dependent relations for quality factors are Q_α = (44 ± 1)f^(0.82±.04) and Q_β = (87 ± 3)f^(0.71±.03). The rate of increase of Q(f) for P and S waves in the Chamoli region is comparable with the other regions of the world. The ratio Q_β/Q_α is greater than one in the region which along with the frequency dependence of quality factors indicates that scattering is an important factor contributing to the attenuation of body waves in the region. A comparison of attenuation relation for S wave estimated here (Q_β = 87f^0.71) with that of coda waves (Q_c = 30f^1.21) obtained by Mandal et al. (2001) for the same region shows that Q_c > Q_β for higher frequencies (>8 Hz) in the region. This indicates a possible high frequency coda enrichment which suggests that the scattering attenuation significantly influences the attenuation of S waves at frequencies >8 Hz. This observation may be further investigated using multiple scattering models. The attenuation relations for quality factors obtained here may be used for the estimation of source parameters and near-source simulation of earthquake ground motion of the earthquakes, which in turn are required for the assessment of seismic hazard in the region.     

Near surface velocity and QS structure of the Quaternary sediment in Bohai basin,China

Heavily populated by Beijing and Tianjin cities, Bohai basin is a seismically active Cenozoic basin suffering from huge lost by devastating earthquakes, such as Tangshan earthquake. The attenuation (Q_P and Q_S) of the surficial Quaternary sediment has not been studied at natural seismic frequency (1?10 Hz), which is crucial to earthquake hazards study. Borehole seismic records of micro earthquake provide us a good way to study the velocity and attenuation of the surficial structure (0?500 m). We found that there are two pulses well separated with simple waveforms on borehole seismic records from the 2006 M_W4.9 Wen'an earthquake sequence. Then we performed waveform modeling with generalized ray theory (GRT) to confirm that the two pulses are direct wave and surface reflected wave, and found that the average v_P and v_S of the top 300 m in this region are about 1.8 km/s and 0.42 km/s, leading to high v_P/v_S ratio of 4.3. We also modeled surface reflected wave with propagating matrix method to constrain QS and the near surface velocity structure. Our modeling indicates that Q_S is at least 30, or probably up to 100, much larger than the typically assumed extremely low Q (~10), but consistent with Q_S modeling in Mississippi embayment. Also, the velocity gradient just beneath the free surface (0?50 m) is very large and velocity increases gradually at larger depth. Our modeling demonstrates the importance of borehole seismic records in resolving shallow velocity and attenuation structure, and hence may help in earthquake hazard simulation.     

Determining Q of near-surface materials from Rayleigh waves

High-frequency (≥2 Hz) Rayleigh wave phase velocities can be inverted to shear (S)-wave velocities for a layered earth model up to 30 m below the ground surface in many settings. Given S-wave velocity (V_S), compressional (P)-wave velocity (V_P), and Rayleigh wave phase velocities, it is feasible to solve for P-wave quality factor Q_P and S-wave quality factor Q_S in a layered earth model by inverting Rayleigh wave attenuation coefficients. Model results demonstrate the plausibility of inverting Q_S from Rayleigh wave attenuation coefficients. Contributions to the Rayleigh wave attenuation coefficients from Q_P cannot be ignored when Vs/V_P reaches 0.45, which is not uncommon in near-surface settings. It is possible to invert Q_P from Rayleigh wave attenuation coefficients in some geological setting, a concept that differs from the common perception that Rayleigh wave attenuation coefficients are always far less sensitive to Q_P than to Q_S. Sixty-channel surface wave data were acquired in an Arizona desert. For a 10-layer model with a thickness of over 20 m, the data were first inverted to obtain S-wave velocities by the multichannel analysis of surface waves (MASW) method and then quality factors were determined by inverting attenuation coefficients.     

Coda wave attenuation in the Parecis Basin,Amazon Craton,Brazil: sensitivity to basement depth

Small local earthquakes from two aftershock sequences in Porto dos Gaúchos, Amazon craton—Brazil, were used to estimate the coda wave attenuation in the frequency band of 1 to 24 Hz. The time-domain coda-decay method of a single backscattering model is employed to estimate frequency dependence of the quality factor (Q _c) of coda waves modeled using Q_c = Q₀ f^hQ_{\rm c} =Q_{\rm 0} f^\eta , where Q ₀ is the coda quality factor at frequency of 1 Hz and η is the frequency parameter. We also used the independent frequency model approach (Morozov, Geophys J Int, 175:239–252, 2008), based in the temporal attenuation coefficient, χ(f) instead of Q(f), given by the equation c(f)=g+\fracpfQ_e \chi (f)\!=\!\gamma \!+\!\frac{\pi f}{Q_{\rm e} }, for the calculation of the geometrical attenuation (γ) and effective attenuation (Q_e^-1 )(Q_{\rm e}^{-1} ). Q _c values have been computed at central frequencies (and band) of 1.5 (1–2), 3.0 (2–4), 6.0 (4–8), 9.0 (6–12), 12 (8–16), and 18 (12–24) Hz for five different datasets selected according to the geotectonic environment as well as the ability to sample shallow or deeper structures, particularly the sediments of the Parecis basin and the crystalline basement of the Amazon craton. For the Parecis basin Q_c = (98±12)f^(1.14±0.08)Q_{\rm c} =(98\pm 12)f^{(1.14\pm 0.08)}, for the surrounding shield Q_c = (167±46)f^(1.03±0.04)Q_{\rm c} =(167\pm 46)f^{(1.03\pm 0.04)}, and for the whole region of Porto dos Gaúchos Q_c = (99±19)f^(1.17±0.02)Q_{\rm c} =(99\pm 19)f^{(1.17\pm 0.02)}. Using the independent frequency model, we found: for the cratonic zone, γ = 0.014 s^− 1, Q_e^-1 = 0.0001Q_{\rm e}^{-1} =0.0001, ν ≈ 1.12; for the basin zone with sediments of ~500 m, γ = 0.031 s^− 1, Q_e^-1 = 0.0003Q_{\rm e}^{-1} =0.0003, ν ≈ 1.27; and for the Parecis basin with sediments of ~1,000 m, γ = 0.047 s^− 1, Q_e^-1 = 0.0005Q_{\rm e}^{-1} =0.0005, ν ≈ 1.42. Analysis of the attenuation factor (Q _c) for different values of the geometrical spreading parameter (ν) indicated that an increase of ν generally causes an increase in Q _c, both in the basin as well as in the craton. But the differences in the attenuation between different geological environments are maintained for different models of geometrical spreading. It was shown that the energy of coda waves is attenuated more strongly in the sediments, Q_c = (78±23)f^(1.17±0.14)Q_{\rm c} =(78\pm 23)f^{(1.17\pm 0.14)} (in the deepest part of the basin), than in the basement, Q_c = (167±46)f^(1.03±0.04)Q_{\rm c} =(167\pm 46)f^{(1.03\pm 0.04)} (in the craton). Thus, the coda wave analysis can contribute to studies of geological structures in the upper crust, as the average coda quality factor is dependent on the thickness of sedimentary layer.     

Trend analyses with river sediment rating curves

Sediment rating curves, which are fitted relationships between river discharge (Q) and suspended‐sediment concentration (C), are commonly used to assess patterns and trends in river water quality. In many of these studies, it is assumed that rating curves have a power‐law form (i.e. C = aQ^b, where a and b are fitted parameters). Two fundamental questions about the utility of these techniques are assessed in this paper: (i) how well to the parameters, a and b, characterize trends in the data, and (ii) are trends in rating curves diagnostic of changes to river water or sediment discharge? As noted in previous research, the offset parameter, a, is not an independent variable for most rivers but rather strongly dependent on b and Q. Here, it is shown that a is a poor metric for trends in the vertical offset of a rating curve, and a new parameter, â, as determined by the discharge‐normalized power function [C = â (Q/Q_GM)^b], where Q_GM is the geometric mean of the Q‐values sampled, provides a better characterization of trends. However, these techniques must be applied carefully, because curvature in the relationship between log(Q) and log(C), which exists for many rivers, can produce false trends in â and b. Also, it is shown that trends in â and b are not uniquely diagnostic of river water or sediment supply conditions. For example, an increase in â can be caused by an increase in sediment supply, a decrease in water supply or a combination of these conditions. Large changes in water and sediment supplies can occur without any change in the parameters, â and b. Thus, trend analyses using sediment rating curves must include additional assessments of the time‐dependent rates and trends of river water, sediment concentrations and sediment discharge. Published 2014. This article is a U.S. Government work and is in the public domain in the USA. Hydrological Processes published by John Wiley & Sons Ltd.     

Frequency-Dependent Attenuation of Coda Waves in the Crust in Southwest Anatolia (Turkey)

The attenuation of coda waves in the earth’s crust in southwest (SW) Anatolia is estimated by using the coda wave method, which is based on the decrease of coda wave amplitude in time and distance. A total of 159 earthquakes were recorded between 1997 and 2010 by 11 stations belonging to the KOERI array. The coda quality factor Q _c is determined from the properties of scattered coda waves in a heterogeneous medium. Firstly, the quality factor Q ₀ (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 for frequencies of 1.5, 3.0, 6.0, 8.0, 12 and 20 Hz. Secondly, the attenuation coefficients (δ) are estimated. The shape of the curve is controlled by the scattering and attenuation in the crustal volume sampled by the coda waves. The average Q _c values vary from 110 ± 15 to 1,436 ± 202 for the frequencies above. The Q ₀ and η values vary from 63 ± 7 to 95 ± 10 and from 0.87 ± 0.03 to 1.04 ± 0.09, respectively, for SW Anatolia. In this region, the average coda Q–f relation is described by Q _c = (78 ± 9)f ^0.98±0.07 and δ = 0.012 km^?1. The low Q ₀ and high η are consistent with a region characterized by high tectonic activity. The Q _c values were correlated with the tectonic pattern in SW Anatolia.     

The estimation of shear-wave statics using in situ measurements in marine near-surface sediments

Shear‐wave statics in marine seismic exploration data are routinely too large to be estimated using conventional techniques. Near‐surface unconsolidated sediments are often characterized by low values of V_s and steep velocity gradients. Minor variations in sediment properties at these depths correspond to variations in the shear‐wave velocity and will produce significant static shifts. It is suggested that a significant proportion of the shear‐wave statics solution can be estimated by performing a separate high‐resolution survey to target near‐surface unconsolidated sediments. Love‐wave, shear‐wave refraction and geotechnical measurements were individually used to form high‐resolution near‐surface shear‐wave velocity models to estimate the shear‐wave statics for a designated survey line. Comparisons with predicted statics revealed that shear‐wave statics could not be estimated using a velocity model predicted by substituting geotechnical measurements into empirical relationships. Empirical relationships represent a vast simplification of the factors that control V_s and are therefore not sufficiently sensitive to estimate shear‐wave statics. Refraction measurements are potentially sensitive to short‐wavelength variations in sediment properties when combined with accurate navigational data. Statics estimated from Love‐wave data are less sensitive, and sometimes smoothed in appearance, since interpreted velocity values represent an average both laterally and vertically over the receiver array and the frequency–depth sensitivity range, respectively. For the survey site, statics estimated from near‐surface irregularities using shear‐wave refraction measurements represent almost half the total statics solution. More often, this proportion will be greater when bedrock relief is less.     

Systematic retrieval of ejecta velocities and gas fluxes at Etna volcano using L-Band Doppler radar

Strombolian-type volcanic activity is characterized by a series of gas bubbles bursting at the top of a magma column and leading to the ejection of lava clots and gas emission at the surface. The quantitative analysis of physical parameters (e.g., velocity, size, and mass fluxes) controlling the emission dynamics of these volcanic products is very important for the understanding of eruption source mechanisms but remains difficult to obtain in a systematic fashion. Ground-based Doppler radar is found to be a very effective tool for measuring ejecta velocities at a high acquisition rate and close to the emission source. We present here a series of measurements carried out at Mt. Etna’s Southeast crater, using an L-band volcanological Doppler radar, during the 4 July 2001 Strombolian eruptions. Doppler radar data are supplemented by the analysis of video snapshots recorded simultaneously. We provide here a set of physical parameters systematically retrieved from 247 Strombolian explosions spanning 15 min and occurring during the paroxysm of the eruption from 21:30 to 21:45 UT. The time-average values give a maximum particle velocity of V_max^p = 94.7±24 \textm/s V_{{\max }}^p = {94}.{7}\pm {24} {\text{m/s}} , a bulk lava jet velocity of V_{\textPW - rad} = 37.6±1.9 \textm/s {V_{{{\text{PW - rad}}}}} = {37}.{6}\pm {1}.{9} {\text{m/s}} , and an initial gas velocity at the source vent of V₀^g = 118.4±36 \textm/s V_0^g = {118}.{4}\pm {36} {\text{m/s}} . The time-averaged particle diameter is found to be about D_{\textPW - rad} = 4.2±2.1 \textcm {D_{{{\text{PW - rad}}}}} = {4}.{2}\pm {2}.{1} {\text{cm}} . The volume and mass gas fluxes are estimated from time-averaged source gas velocities over the sequence duration at Q_v^g = 3 - 11 ×10³\textm³\text/s Q_v^g = {3} - {11} \times {1}{0^{{3}}}{{\text{m}}^{{3}}}{\text{/s}} and Q_m^g = 0.5 - 2 ×10³\textkg/s Q_m^g = 0.{5} - {2} \times {1}{0^{{3}}}{\text{kg/s}} , respectively.     

Experimental verification of the fracture density and shear‐wave splitting relationship using synthetic silica cemented sandstones with a controlled fracture geometry

We present laboratory ultrasonic measurements of shear‐wave splitting from two synthetic silica cemented sandstones. The manufacturing process, which enabled silica cementation of quartz sand grains, was found to produce realistic sandstones of average porosity 29.7 ± 0.5% and average permeability 29.4 ± 11.3 mD. One sample was made with a regular distribution of aligned, penny‐shaped voids to simulate meso‐scale fractures in reservoir rocks, while the other was left blank. Ultrasonic shear waves were measured with a propagation direction of 90° to the coincident bedding plane and fracture normal. In the water saturated blank sample, shear‐wave splitting, the percentage velocity difference between the fast and slow shear waves, of <0.5% was measured due to the bedding planes (or layering) introduced during sample preparation. In the fractured sample, shear‐wave splitting (corrected for layering anisotropy) of 2.72 ± 0.58% for water, 2.80 ± 0.58% for air and 3.21 ± 0.58% for glycerin saturation at a net pressure of 40 MPa was measured. Analysis of X‐ray CT scan images was used to determine a fracture density of 0.0298 ± 0.077 in the fractured sample. This supports theoretical predictions that shear‐wave splitting (SWS) can be used as a good estimate for fracture density in porous rocks (i.e., SWS = 100ε_f, where ε_f is fracture density) regardless of pore fluid type, for wave propagation at 90° to the fracture normal.     

Validation of 3-D basin structure models for long-period ground motion simulation in the Osaka basin, western Japan

We studied the applicability of two types of existing three-dimensional (3-D) basin velocity structure models of the Osaka basin, western Japan for long-period ground motion simulations. We synthesized long-period (3–20 s) ground motions in the Osaka basin during a M6.5 earthquake that occurred near the hypothetical Tonankai earthquake source area, approximately 200 km from Osaka. The simulations were performed using a 3-D finite-difference method with nonuniform staggered grids using the two basin velocity structure models. To study the ground motion characteristics inside the basin, we evaluated the wave field inside the basin using the transfer functions derived from the synthetics at the basin and a reference rock site outside the basin. The synthetic waveforms at the basin site were obtained by a convolution of the calculated transfer function and the observed waveform at the reference rock site. First, we estimated the appropriate Q values for the sediment layers. Assuming that the Q value depends on the S wave velocity V _S and period T, it was set to Q = (1/3V _S)(T ₀/T) where V _S is in m/s and the reference period T ₀ is 3.0 s. Second, we compared the synthetics and the observations using waveforms and pseudovelocity response spectra, together with a comparison of the velocity structures of the two basin models. We also introduced a goodness-of-fit factor to the pseudovelocity response spectra as an objective index. The synthetics of both the models reproduced the observations reasonably well at most of the stations in the central part the basin. At some stations, however, especially where the bedrock depth varies sharply, there were noticeable discrepancies in the simulation results of the models, and the synthetics did not accurately reproduce the observation. Our results indicate that the superiority of one model over the other cannot be determined and that an improvement in the basin velocity structure models based on simulation studies is required, especially along the basin edges. We also conclude that our transfer function method can be used to examine the applicability of the basin velocity structure models for long-period ground motion simulations.