Seismic sources and attenuation properties at the Campi Flegrei volcanic area

Microearthquake digital data collected at Campi Flegrei during the recent (1982–1985) ground uplift episode have been analyzed in order to infer source and medium seismic properties. The main results obtained from these analyses are:

1. Hypocenter distribution and the size of the seismic zone do not change with time and do not depend on the ground uplift rate. Events occurred clustered in time with no simple causal relations between the cluster occurrences and their energy.
2. Anelastic attenuation does not depend strongly on frequency, showing a constant pattern at high frequencies. The observed values of low and high frequency attenuation, due to the short source receiver distances, do not seriously affect the spectral content of signals radiated by the sources.
3. A constant Brune stress drop pattern (～4–5 bars) as a function of seismic moment is observed. This indicates that the manner of fracturing is almost independent on magnitude of earthquakes (hypothesis of self-similarity (Aki, 1967)). Seismic processes in a prefractured medium can explain the observed small stress drop values.
4. Focal mechanisms from moment tensor estimates show that radiation patterns are mostly well interpreted in terms of double couple source models.
5. The scaling of peak ground motion parameters (A _max andV _max vs seismic moment) can be explained by an ω² source model (constant stress drop) multiplied by an exponential function with a small decay parameter, which takes into account the measured attenuation.

These results support the hypothesis of earthquakes generated by simple shear fractures along prefractured structures as a response to changes in the stress field due to the ground deformation.     

Multiple scattering of seismic waves in fractured media: Velocity and effective attenuation of the coherent components ofP waves andS waves

We show that the multiple scattering by small fractures of seismic waves with wavelengths long compared to the fracture size and fracture spacing is indistinguishable from multiple-scattering effects produced by regular porosity, except for an orientation factor due to fracture alignment. The fractures reduce theP-wave andS-wave velocities and produce an effective attenuation of the coherent component of the seismic waves. The attenuation corresponds to 1000/Q of about unity for a Gaussian spectrum of fractures, and it varies with frequencyf asf ³. For a Kolmogorov spectrum of fractures of spectral index the attenuation is an order of magnitude or so larger and varies with frequency asf ^3-v The precise degree of attenuation depends upon the matrix properties, the fracture porosity, the degree of fracture anisotropy, the type of fluid filling the fractures, and the incidence angle of the wave.For fracture porosities less than about 15% theP-wave andS-wave velocities are decreased by the order of 5–10% with a lesser dependence on the type of fluid filling the fractures (gas, oil, or brine) and with a dependence on both the degree of anisotropy and the incident angle made by the wave. The tendency of fractures to occur perpendicularly to bedding suggests that the best way to measure seismically fractured rock behavior in situ is by using the travel-time delay and reflection amplitude. As both the offset and the azimuth of receivers vary from a shot, the travel-time delay and reflection amplitude should both show an elliptical pattern of behavior—the travel-time delay in response to the varying seismic speed, and the reflection amplitude in response to angular variations in the multiple scattering. Observations of attenuation at several frequencies should permit (a) determination of the spectrum of fractures (Gaussian versus Kolmogorovian) and (b) determination of the contribution of viscous damping to the effective attenuation.     

Attenuation coefficients of Rayleigh and <Emphasis Type="Italic">Lg</Emphasis> waves

Analysis of the frequency dependence of the attenuation coefficient leads to significant changes in interpretation of seismic attenuation data. Here, several published surface-wave attenuation studies are revisited from a uniform viewpoint of the temporal attenuation coefficient, denoted by χ. Theoretically, χ( f) is expected to be linear in frequency, with a generally non-zero intercept γ?=?χ(0) related to the variations of geometrical spreading, and slope dχ/df = π/Q _e caused by the effective attenuation of the medium. This phenomenological model allows a simple classification of χ( f) dependences as combinations of linear segments within several frequency bands. Such linear patterns are indeed observed for Rayleigh waves at 500–100-s and 100–10-s periods, and also for Lg from ~2 s to ~1.5 Hz. The Lg χ( f) branch overlaps with similar linear branches of body, Pn, and coda waves, which were described earlier and extend to ~100 Hz. For surface waves shorter than ~100 s, γ values recorded in areas of stable and active tectonics are separated by the levels of \(\gamma _{D} \approx 0.2 \times 10^{-3}\) s^???1 (for Rayleigh waves) and 8 ×10^???3 s^???1 (for Lg). The recently recognized discrepancy between the values of Q measured from long-period surface waves and normal-mode oscillations could also be explained by a slight positive bias in the geometrical spreading of surface waves. Similarly to the apparent χ, the corresponding linear variation with frequency is inferred for the intrinsic attenuation coefficient, χ _i , which combines the effects of geometrical spreading and dissipation within the medium. Frequency-dependent rheological or scattering Q is not required for explaining any of the attenuation observations considered in this study. The often-interpreted increase of Q with frequency may be apparent and caused by using the Q-based model of attenuation and following preferred Q( f) dependences while ignoring the true χ( f) trends within the individual frequency bands.     

Investigation on regional attenuation of Vrancea (Romania) intermediate-depth earthquakes

This paper aims at investigating possible regional attenuation patterns in the case of Vrancea(Romania) intermediate-depth earthquakes.Almost 500 pairs of horizontal components recorded during 13 intermediate-depth Vrancea earthquakes are employed in order to evaluate the regional attenuation patterns.The recordings are grouped according to the azimuth with regard to the Vrancea seismic source and subsequently,Q models are computed for each azimuthal zone assuming similar geometrical spreading.Moreover,the local soil amplification which was disregarded in a previous analysis performed for Vrancea intermediate-depth earthquakes is now clearly evaluated.The results show minor differences between the four regions situated in front of the Carpathian Mountains and considerable differences in attenuation of seismic waves between the forearc and backarc regions(with regard to the Carpathian Mountains).Consequently,an average Q model of the type Q(f) = 115×f~(1.25) is obtained for the four forearc regions,while a separate Q model of the type Q(f) = 70×f~(0.90) is computed for the backarc region.These results highlight the need to evaluate the seismic hazard of Romania by using ground motion models which take into account the different attenuation between the forearc/backarc regions.     

The 1986 Kermadec earthquake and its relation to plate segmentation

To evaluate the tectonic significance of the October 20, 1986 Kermadec earthquake (M _w =7.7), we performed a comprehensive analysis of source parameters using surface waves, body waves, and relocated aftershocks. Amplitude and phase spectra from up to 93 Rayleigh waves were inverted for centroid time, depth, and moment tensor in a two-step algorithm. In some of the inversions, the time function was parameterized to include information from the body-wave time function. The resulting source parameters were stable with respect to variations in the velocity and attenuation models assumed, the parameterization of the time function, and the set of Rayleigh waves included. The surface wave focal mechanism derived (=275°, =61°, =156°) is an oblique-compressional mechanism that is not easy to interpret in terms of subduction tectonics. A seismic moment of 4.5×10²⁰ N-m, a centroid depth of 45±5 km, and a centroid time of 13±3 s were obtained. Directivity was not resolvable from the surface waves. The short source duration is in significant contrast to many large earthquakes.We performed a simultaneous inversion ofP andSH body waves for focal mechanism and time function. The focal mechanism agreed roughly with the surface wave mechanism. Multiple focal mechanisms remain a possibility, but could not be resolved. The body waves indicate a short duration of slip (15 to 20 s), with secondary moment release 60s later. Seismically radiated energy was computed from the body-wave source spectrum. The stress drop computed from the seismic energy is about 30 bars. Sixty aftershocks that occurred within three months of the mainshock were relocated using the method of Joint Hypocentral Determination (JHD). Most of the aftershocks have underthrusting focal mechanisms and appear to represent triggered slip on the main thrust interface. The depth, relatively high stress drop, short duration of slip, and paucity of true aftershocks are consistent with intraplate faulting within the downgoing plate. Although it is not clear on which nodal plane slip occurred, several factors favor the roughly E-W trending plane. The event occurred near a major segmentation in the downgoing plate at depth, near a bend in the trench, and near a right-lateral offset of the volcanic are by 80 km along an E-W direction. Also, all events in the region from 1977 to 1991 with CMT focal mechanisms similar to that of the Mainshock occurred near the mainshock epicenter, rather than forming an elongate zone parallel to the trench as did the aftershock activity. We interpret this event as part of the process of segmentation or tearing of the subducting slab. This segmentation appears to be related to the subduction of the Louisville Ridge, which may act as an obstacle to subduction through its buoyancy.     

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.     

Stress wave patterns and the size of the non-elastic zones of an explosive source

Summary The travel-time curves of characteristic moments of stress wave patterns were investigated. The arrival time of the onset of stress wavet ₀, of the maximum amplitudet ₁ and of the moment terminating the pressure part (the first half-wave) of this wavet ₂ were taken to be representative. These travel-time curves were used to determine the velocities of propagation of these moments (V ₀,V ₁ andV ₂) as functions of the distance from the source. According to their variations it is possible to appoint the size of the cavity created by the explosion and to determine the distance at which the elastic source surface of stress waves is to be found. The radius of the cavity is given by 1.5 times the distance at whichV ₁ separates fromV ₀. The elastic source surface of stress waves is defined by the distance from whichV ₀ is constant and at whichV ₁ andV ₂ have minimum values. These two distances determine the points at which the character of the stress wave pattern changes: from shock wave in gases, through shock wave in a solid medium to a seismic wave.Paper presented at the General Assembly of the European Seismological Commission, Luxembourg, Sept. 21–29, 1970.     

Deep drainage sensitivity to climate,edaphic factors,and woody encroachment,Oklahoma, USA

Because groundwater is Earth's largest pool of freshwater, understanding the sensitivity of deep drainage to climate, soils, and land cover is critical in managing water resources. To better understand controls on this critical flux in the context of woody encroachment, we determined the sensitivity of deep drainage to climate, soil texture, soil compaction, rooting depth, growing season duration, and plant–water stress response using Hydrus‐1D to simulate deep drainage. To evaluate the simulation results, we compared these results with ground measurements at two anchor sites. At both anchor sites, Hydrus‐1D predictions of deep drainage matched measured values within the errors inherent in ground measurements. Sensitivity analysis suggested greatest sensitivity of deep drainage to climate (24 mm yr^?1) and rooting depth (12 mm yr^?1), moderate sensitivity to growing season duration (5 mm yr^?1) and soil texture (4 mm yr^?1), and lowest sensitivity to topsoil compaction and plant–water stress response (3 mm yr^?1). The sensitivity analysis indicated the relative importance of the plant‐related factors considered, which, in decreasing order, were rooting depth, growing season duration, and plant–water stress response – factors that change concomitantly as a result of forestation or woody encroachment. Further ground‐truth measurements of woody encroachment effects on deep drainage are needed to confirm or refine the results of this simulation modelling study. Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

Quantitative estimates of interplate coupling inferred from outer rise earthquakes

Interplate coupling plays an important role in the seismogenesis of great interplate earthquakes at subduction zones. The spatial and temporal variations of such coupling control the patterns of subduction zone seismicity. We calculate stresses in the outer rise based on a model of oceanic plate bending and coupling at the interplate contact, to quantitatively estimate the degree of interplate coupling for the Tonga, New Hebrides, Kurile, Kamchatka, and Marianas subduction zones. Depths and focal mechanisms of outer rise earthquakes are used to constrain the stress models. We perform waveform modeling of body waves from the GDSN network to obtain reliable focal depth estimates for 24 outer rise earthquakes. A propagator matrix technique is used to calculate outer rise stresses in a bending 2-D elastic plate floating on a weak mantle. The modeling of normal and tangential loads simulates the total vertical and shear forces acting on the subducting plate. We estimate the interplate coupling by searching for an optimal tangential load at the plate interface that causes the corresponding stress regime within the plate to best fit the earthquake mechanisms in depth and location.We find the estimated mean tangential load over 125–200 km width ranging between 166 and 671 bars for Tonga, the New Hebrides, the Kuriles, and Kamchatka. This magnitude of the coupling stress is generally compatible with the predicted shear stress at the plate contact from thermal-mechanical plate models byMolnar andEngland (1990), andVan den Buekel andWortel (1988). The estimated tectonic coupling,F _tc , is on the order of 10¹²–10¹³ N/m for all the subduction zones.F _tc for Tonga and New Hebrides is about twice as high as in the Kurile and Kamchatka arcs. The corresponding earthquake coupling forceF _ec appears to be 1–10% of the tectonic coupling from our estimates. There seems to be no definitive correlation of the degree of seismic coupling with the estimated tectonic coupling. We find that outer rise earthquakes in the Marianas can be modeled using zero tangential load.     

Laboratory measurements of spatial fluctuation and attenuation of elastic waves by scattering due to random heterogeneities

To study the effects of strong scattering on elastic waves, spatial fluctuation and scattering attenuation ofP waves were examined by laboratory experiments for 2-D models of random media approximately characterized by a triangular correlation function in the range of 2<ka<33, wherek is the wave number anda is the correlation distance of the heterogeneities, i.e., the heterogeneity size. The results obtained are as follows: (1) Forka>10, both the intensity and the correlation distance of the amplitude fluctuation are approximate for any phase of theP-wave train. The correlation distance nearly agrees with the heterogeneity size. These fluctuation properties are quite consistent with the theoretical prediction by the forward-scattering approximation. (2) For 3<ka<6, the fluctuation intensity becomes stronger in later phases of theP-wave train. This shows that scattering is approximately isotropic, and therefore, the scattered energy increases with time within theP-wave train. The correlation distance of the amplitude fluctuation disagrees with the heterogeneity size, and it shows a frequency-dependent property decreasing from 7a to 4a with the increase ofka from 3 to 6. These properties for 3<ka<6 have not yet been predicted theoretically. (3) Forka<3, though the fluctuation is considerably smaller compared with that ofka>10 and 3<ka<6, the fluctuation property is considered similar to that of 3<ka<6. (4) The observed scattering attenuation,Q ^–1, increases withka forka<3, has a peak aroundka=35, and then decreases withka. (5) When _min = 15° and = 0.075, the theoreticalQ ^–1 curve, predicted by the approximate theory of Wu, roughly matches the observedQ ^–1 values, where _min is the minimum scattering angle measured from the propagation direction of theP waves and is the rms of fractional velocity fluctuation. This suggests that the energy scattered in the range of >15° is lost from theP waves, while the energy scattered in the range of <15° is retained; and that the approximate theory overestimates by about three times the value of the model media used owing to the neglect of multiple scattering. (6) When the size of velocity heterogeneities responsible for forward scattering at 3<ka<6 is estimated from the _min value of 15° on the basis of Wu's theory, it nearly agrees with the correlation distance for the initial phase of theP-wave train.     

Measuring changes in fracture properties from temporal variations in anisotropic attenuation of microseismic waveforms

We investigate fracture‐induced attenuation anisotropy in a cluster of events from a microseismic dataset acquired during hydraulic fracture stimulation. The dataset contains 888 events of magnitude ?3.0 to 0.0. We use a log‐spectral‐amplitude‐ratio method to estimate change in over a half‐hour time period where fluid is being injected and an increase in fracturing from S‐wave splitting analysis has been previously inferred. A Pearson's correlation analysis is used to assess whether or not changes in attenuation with time are statistically significant. P‐waves show no systematic change in during this time. In contrast, S‐waves polarised perpendicular to the fractures show a clear and statistically significant increase with time, whereas S‐waves polarised parallel to the fractures show a weak negative trend. We also compare between the two S‐waves, finding an increase in with time. A poroelastic rock physics model of fracture‐induced attenuation anisotropy is used to interpret the results. This model suggests that the observed changes in t* are related to an increase in fracture density of up to 0.04. This is much higher than previous estimates of 0.025 ± 0.002 based on S‐wave velocity anisotropy, but there is considerably more scatter in the attenuation measurements. This could be due to the added sensitivity of attenuation measurement to non‐aligned fractures, fracture shape, and fluid properties. Nevertheless, this pilot study shows that attenuation measurements are sensitive to fracture properties such as fracture density and aspect ratio.     

Travel times for a linear variation of the velocity of bodily waves

Summary The expression of the travel timet of bodily waves propagating in a layered earth, as a function of the emerging angle is obtained, to complement the one given in a previous paper [1]²) of the epicentral distance . The waves velocityV is again assumed to vary linearly with the depth. Having both expressions it becomes possible to construct a theoretical functiont=t() and to explore the true values ofV by adjusting it to the experimental travel time tables.     

Source parameters and scaling relationships in the Friuli-Venezia Giulia (Northeastern Italy) region

We estimated the source parameters of 53 local earthquakes (2.0<M_L<5.7) of the Friuli-Venezia Giulia (Northeastern Italy) area, recorded by the short-period local seismic network of the Istituto Nazionale di Oceanografia e Geofisica Sperimentale (OGS), in the period 1995-2003. Data were selected on the basis of high quality locations and focal mechanisms. Standard H/V spectral ratios (HVRS) of the three-component stations of the network were performed in order to assess local amplifications, and only stations showing HVRS not exceeding two were considered for the source parameters estimation. Both velocity and acceleration data were used to compute the SH-wave spectra. Observed spectra were corrected for attenuation effects using an independent regional estimate of the quality factor Q and a station dependent estimate of the spectral decay parameter k. Only earthquakes with M_L>3.0 recorded with a sampling rate of 125 cps were used to compute k, thus allowing to visualize a linear trend of the high frequency acceleration spectrum up to 40-50 Hz. SH-wave spectra, corrected for attenuation, showed an ω⁻² shape allowing a good fit with the Brune model. Seismic moments and Brune radii ranged between 1.5×10¹² and 1.1×10¹⁷  N m and between 0.1 and 2.7 km respectively. We obtained M_o=1.1×10¹⁷  N m for the seismic moment of the Kobarid (SLO) main shock, in good agreement with the Harvard CMT solution (M_o=3.5×10¹⁷  N m). Brune stress drops were confined to the range from 0.07 to 5.31 MPa, with an average value of 0.73 MPa and seem to be approximately constant over five orders of magnitude of seismic moment. Radiated seismic energy computed from two nearby stations scales with seismic moment according to , and apparent stress values are between 0.02 and 4.26 MPa. The observed scatter of Brune stress drop data allowed to hypothesize a scaling relation between seismic moment and corner frequency in order to accommodate both Brune stress drop and apparent stress scalings. No systematic differences are evidenced between stress parameters of earthquakes with different focal mechanisms. As a consequence, a relation of the seismic stress release with the strength of rocks can be hypothesized. A high correlation (r>0.9) of Brune stress drop is found with both apparent stress and RMS stress drop, according to and respectively.     

Seismic Wave Attenuation in the Greater Cairo Region, Egypt

In the present study, a digital waveform dataset of 216 local earthquakes recorded by the Egyptian National Seismic Network (ENSN) was used to estimate the attenuation of seismic wave energy in the greater Cairo region. The quality factor and the frequency dependence for Coda waves and S-waves were estimated and clarified. The Coda waves (Q _c) and S-waves (Q _d) quality factor were estimated by applying the single scattering model and Coda Normalization method, respectively, to bandpass-filtered seismograms of frequency bands centering at 1.5, 3, 6, 12, 18 and 24?Hz. Lapse time dependence was also studied for the area, with the Coda waves analyzed through four lapse time windows (10, 20, 30 and 40?s). The average quality factor as function of frequency is found to be Q _c?=?35?±?9f ^0.9±0.02 and Q _d?=?10?±?2f ^0.9±0.02 for Coda and S-waves, respectively. This behavior is usually correlated with the degree of tectonic complexity and the presence of heterogeneities at several scales. The variation of Q _c with frequency and lapse time shows that the lithosphere becomes more homogeneous with depth. In fact, by using the Coda Normalization method we obtained low Q _d values as expected for a heterogeneous and active zone. The intrinsic quality factor (Q _i ^?1 ) was separated from the scattering quality factor (Q _s ^?1 ) by applying the Multiple Lapse Time Domain Window Analysis (MLTWA) method under the assumption of multiple isotropic scattering with uniform distribution of scatters. The obtained results suggest that the contribution of the intrinsic attenuation (Q _i ^?1 ) prevails on the scattering attenuation (Q _s ^?1 ) at frequencies higher than 3?Hz.     

Analysis of various large-scale disturbances and the associated zonal mean motions in the atmosphere

Summary In this article, we present a scale analysis of planetary waves, extended long waves, and long waves. (We mean the extended long waves to be the disturbances whose east-west length is of order 10⁶ m and north-south extension 10⁷ m). We find for the extended long waves the two terms, the interaction between kinetic and available potential energy of the disturbances, and the interaction between the zonal mean available potential energy, and the eddy available potential energy, are of two orders of magnitude larger than the kinetic energy interaction between the disturbances and the associated zonal mean flow. This theoretical result concerning the relative importance of the various interaction terms may be of use in explaining the observational findings thus far available.It is also shown theoretically that the kinetic energy interaction between the planetary waves, the horizontal size of which is 10⁷ m, and the long waves, whose horizontal size is 10⁶ m, is of the same order as the interaction of kinetic energy between the zonal mean motion and the disturbances. This agrees fairly well with the observational estimates thus far obtained.     

Threshold seismic energy and liquefaction distance limit during the 2008 Wenchuan earthquake

Estimating the possible region of liquefaction occurrence during a strong earthquake is highly valuable for economy loss estimation, reconnaissance efforts and site investigations after the event. This study identified and compiled a large amount of liquefaction case histories from the 2008 Wenchuan earthquake, China, to investigate the relationship between the attenuation of seismic wave energy and liquefaction distance limit during this earthquake. Firstly, we introduced the concept of energy absorption ratio, which is defined as the absorbed energy of soil divided by the imparted energy of seismic waves at a given site, and the relationship between the energy absorption ratio and the material damping ratio was established based on shear stress–strain loop of soil element and the seismic wave propagation process from the source to the site. Secondly, the threshold imparted seismic energy of liquefaction was obtained based on existing researches of absorbed energy required to trigger liquefaction of sandy soils and the ground motion attenuation characteristics of the 2008 Wenchuan earthquake, and the liquefaction distance limit of this earthquake was estimated according to the proposed magnitude–energy–distance relationship. Finally, the field liquefaction database of 209 sites of the 2008 Wenchuan earthquake was used to validate such an estimation, and the field observed threshold imparted seismic energy to cause liquefaction in recent major earthquakes worldwide was back-analyzed to check the predictability of the present method, and several possible mechanisms were discussed to explain the discrepancy between the field observations and the theoretical predictions. This study indicates that seismic energy attenuation and liquefaction distance limit are regional specific and earthquake dependent, and 382 J/m³ is the average level of threshold imparted seismic energy to cause liquefaction for loose saturated sandy soils, and the corresponding liquefaction distance limit is approximately 87.4 km in fault distance for a M_w?=?7.9 event in the Chengdu Plain. The proposed regional energy attenuation model and threshold imparted seismic energy may be considered as an approximate tool in evaluating the liquefaction hazard during potential earthquakes in this area.     

Effects of geomagnetic disturbances in power spectra of atmospheric waves in the dynamo region of the ionosphere

The dynamics of wave disturbances in the ionospheric E region in the band of periods of thermal tidal waves and waves of planetary scales (T = 48, 72, and 192 h) has been studied based on the variations in the horizontal component of the geomagnetic field, observed at Paratunka and Barrow observatories in September–October 1999. It has been found that, at midlatitudes during high geomagnetic activity, the intensity of oscillations in the power spectra with T = 24 and 12 h varies with a periodicity of 16 days different from the periodicity of changes in the ΣKp index. The maximal deviations of these periods from the values under quiet conditions coincide with the maximal changes in the ΣKp index. The variations in the 48–192 h band of periods (especially with T ～192 h) intensify simultaneously with increasing geomagnetic activity. The intensity of this harmonic is several times as high as that of the harmonic with T ～ 24 h. The periodicity of changes in the harmonics intensity within the 48–192 h band coincides with the periodicity of changes in the ΣKp index. In the polar ionosphere, the effect of high geomagnetic activity is observed as an increase in the variations with a quasi-period of T ～ 24 h and as an appearance of variations in the 48–192 h band with the periodicity coinciding with the maximums in the ΣKp index variations.     

Scaling of peak ground motions from digital recordings of small earthquakes at Campi Flegrei,southern Italy

The dependence of peak ground acceleration and velocity on seismic moment is studied for a set of small earthquakes (0.7<M _L<3.2) recorded digitally at distances of a few km in the Campi Flegrei volcanic area near Naples, Italy, during the ground uplift episode of 1982–1984. Numerical simulations, using the -square spectral model with constant stress drop and ane ^–kf high frequency decay, fit well both the velocity and acceleration data for an averagek=0.015. The observed ground motions in the 1–24 Hz frequency band appear to consist of radiation from simple sources modified only slightly by attenuation effects. Moreover, the scaling of peak values agrees closely with those determined in nonvolcanic areas, once the difference in stress drop is taken into account.