Three-dimensional VP and VP /VS models of the upper crust in the Friuli area (northeastern Italy)

3-D images of P velocity and P - to S -velocity ratio have been produced for the upper crust of the Friuli area (northeastern Italy) using local earthquake tomography. The data consist of 2565 P and 930 S arrival times of high quality. The best-fitting V _P and V _P / V _S 1-D models were computed before the 3-D inversion. V _P was measured on two rock samples representative of the investigated upper layers of the Friuli crust. The tomographic V _P model was used for modelling the gravity anomalies, by converting the velocity values into densities along three vertical cross-sections. The computed gravity anomalies were optimized with respect to the observed gravity anomalies. The crust investigated is characterized by sharp lateral and deep V _P and V _P / V _S anomalies that are associated with the complex geological structure. High V _P / V _S values are associated with highly fractured zones related to the main faulting pattern. The relocated seismicity is generally associated with sharp variations in the V _P / V _S anomalies. The V _P images show a high-velocity body below 6 km depth in the central part of the Friuli area, marked also by strong V _P / V _S heterogeneities, and this is interpreted as a tectonic wedge. Comparison with the distribution of earthquakes supports the hypothesis that the tectonic wedge controls most of the seismicity and can be considered to be the main seismogenic zone in the Friuli area.     

On the linear relation between mb and Ms for discrimination between explosions and earthquakes

Summary. The statistical capability of the m _b: M _s discriminant for the discrimination of earthquake and explosion populations is examined by application of discriminant functions to a group of 83 explosions and 72 earthquakes in Eurasia. Equations are derived for the probability that an event is an earthquake or an explosion. The positive sign of DIS in the decision index equation, DIS _i = 34.3383 – 11.9569 mb _t + 7.1161 M _si , indicates that the event i is an earthquake. Its negative sign indicates that event i is an explosion. The probability of correct classification for an event, P _i , is related to its DIS _i value, by P _i = [1-exp (DIS _i )]⁻¹, where a large, positive DIS indicates a high probability that an event is an earthquake and a large, negative DIS indicates a high probability that an event is an explosion. The discrimination line M _s = 1.680 m _b– 4.825, or m _b= 0.595 M _s+ 2.872 very successfully separates the explosion population from the earthquake population. The points on this line have an equal chance of being an earthquake or an explosion; moreover, for any event, the distance parallel to the M _s-axis from the point representing that event in the m _b: M _s plane to this line is a measure of the probability for the correct classification of that event.     

Stacking for the frequencies and Qs of 0S0 and 1S0

Summary. Using nine IDA records for the Indonesian earthquake of 1977 August 19, we have formed an optimal linear combination of the records and have measured the frequency and Q of ₀ S ₀ and ₁ S ₀. The frequency was measured using the moment ratio method. The attenuation was measured by the minimum width method and by the time-lapse method. The frequency and attenuation were measured simultaneously by varying them to obtain a best fit to the data. A 2000-hr stack, the sum of nine individual records, for ₀ S ₀ gave a frequency of 0.814664 mHz±4 ppm. The values for the Q of ₀ S ₀ for the three different methods of measurement were 5600,5833 and 5700, respectively. The error in the estimates of Q ^-1 is about 5 per cent for the minimum power method. For ₁ S ₀ a 300-hr stack yielded a frequency of 1.63151 mHz±30 ppm. The values of Q for this mode were 1960, 1800 and 1850, respectively, with an error in Q ^-1 of about 12 per cent for the minimum power method.     

A new derivation of the ratio Q(wobble)/Qμ (mantle)

Summary An extension of the Love-Larmor theory to a low-loss unelastic earth model, leads to the surprisingly simple approximation
   
where τ_s= 447.4 sidereal day is the static wobble period, τ_R= 306 sidereal day is the rigid-earth wobble period and τ_w= 433 sidereal day is the observed Chandler period. Q _W, Q _μ are the respective average Q values of the wobble and the Earth's mantle at τ_W. The known numerical factor F is only slightly dependent on the Earth structure.     

Run-up from impact tsunami

We report on calculations of the on-shore run-up of waves that might be generated by the impact of subkilometre asteroids into the deep ocean. The calculations were done with the COULWAVE code, which models the propagation and shore-interaction of non-linear moderate- to long-wavelength waves  ( kh < π)  using the extended Boussinesq approximation. We carried out run-up calculations for several different situations: (1) laboratory-scale monochromatic wave trains onto simple slopes; (2) 10–100 m monochromatic wave trains onto simple slopes; (3) 10–100 m monochromatic wave trains onto a compound slope representing a typical bathymetric profile of the Pacific coast of North America; (4) time-variable scaled trains generated by the collapse of an impact cavity in deep water onto simple slopes and (5) full-amplitude trains onto the Pacific coast profile. For the last case, we also investigated the effects of bottom friction on the run-up. For all cases, we compare our results with the so-called 'Irribaren scaling': The relative run-up   R / H ₀=ξ= s ( H ₀/ L ₀)^−1/2  , where the run-up is   R , H ₀  is the deep-water waveheight, L ₀ is the deep-water wavelength, s is the slope and ξ is a dimensionless quantity known as the Irribaren number. Our results suggest that Irribaren scaling breaks down for shallow slopes   s ≤ 0.01  when  ξ < 0.1 − 0.2  , below which   R / H ₀  is approximately constant. This regime corresponds to steep waves and very shallow slopes, which are the most relevant for impact tsunami, but also the most difficult to access experimentally.     

The attenuation mechanism of seismic waves in northwestern Himalayas

We analysed local earthquake waveforms recorded on a broad-band seismic network in northwestern Himalayas to compute the intrinsic and scattered attenuation parameters from coda waves. Similar to other tectonically active and heterogeneous regions, attenuation-frequency relation for western Himalaya is   Q ⁻¹ _c = (113 ± 7)  f ^(1.01±0.05)  where   Q_c   is the coda Q parameter. Intrinsic  ( Q ⁻¹ _i )  and scattering  ( Q ⁻¹ _s )  attenuations was separated using   Q_c   and direct S -wave Q data  ( Q_d )  . It is observed that estimated   Q ⁻¹ _c   is close to   Q ⁻¹ _i   and both of them are much larger than   Q ⁻¹ _s   suggesting that coda decay is predominantly caused by intrinsic attenuation. At higher frequencies, both the attenuation parameters   Q_c   and,   Q_d   are similar indicating that coda is predominantly composed of back-scattered S waves at these frequencies.     

Core–mantle boundary deformations and J2 variations resulting from the 2004 Sumatra earthquake

The deformation at the core–mantle boundary produced by the 2004 Sumatra earthquake is investigated by means of a semi-analytic theoretical model of global coseismic and postseismic deformation, predicting a millimetric coseismic perturbation over a large portion of the core–mantle boundary. Spectral features of such deformations are analysed and discussed. The time-dependent postseismic evolution of the elliptical part of the gravity field ( J ₂) is also computed for different asthenosphere viscosity models. Our results show that, for asthenospheric viscosities smaller than 10¹⁸ Pa s, the postseismic J ₂ variation in the next years is expected to leave a detectable signal in geodetic observations.     

Seismic motion in urban sites consisting of blocks in welded contact with a soft layer overlying a hard half-space

We address the problem of the response to a seismic wave of an urban site consisting of   N _b   blocks overlying a soft layer underlain by a hard substratum. The results of a theoretical analysis, appealing to a space–frequency mode-matching (MM) technique, are compared to those obtained by a space–time finite-element (FE) technique. The two methods are shown to give rise to the same prediction of the seismic response for   N _b = 1  , 2 and 40 blocks. The mechanism of the interaction between blocks and the ground, as well as that of the mutual interaction between blocks, are studied. It is shown, in the first part of this paper, that the presence of a small number of blocks modifies the seismic disturbance in a manner which evokes qualitatively, but not quantitatively, what was observed during the 1985 Michoacan earthquake in Mexico City. Anomalous earthquake response at a much greater level, in terms of duration, peak and cumulative amplitude of motion, is shown, by a theoretical and numerical analysis in the second part of this paper, to be induced by the presence of a large (≥10) number of identical equi-spaced blocks that are present in certain districts of many cities.     

The solution of the inverse gravity potential problem for cylindrical bodies

Summary. The inverse gravity potential problem consists in the determination of the form and the density of the body by its exterior gravity potential. We describe two similar classes of bodies for which this problem has a unique constructive solution.
(1) The first class contains the cylindrical bodies with finite length, arbitrary form of section and ρ( R , ø, z) =ρ₁( z )ρ₂( R , ø) density distribution, where z is the cylindrical coordinate; R , ø are the polar coordinates in a section plane. This class is important for prospecting geophysics in that it allows us to determine in a unique and constructive way, the function ρ₁( R , ø), the length, form and orientation of the cylinder if we know the function ρ₁( z ) and the exterior potential. The classical moment problem of functions is the basis for the solution of this problem.
(2) The analogous problem for the class of the spherical cylinders, or bodies bounded by arbitrary similar sections of two different concentric spheres and the radial lateral surface, appears when bodies of planetary size are studied. (An example of these bodies would be the Moon mascons.) The density distribution of these cylinders is ρ(τ, θ, ø) =ρ₁(τ)ρ₂(θ, ø) where τ, θ, ø are the spherical coordinates. The function ρ₁(θ, ø), length and form of spherical sections can be uniquely determined by exterior potential if we know the function ρ₁(τ). We propose a new constructive method for harmonic continuation of the gravity potential into the region containing the perturbing masses for the solution of the problem.     

Seismic image of a CO2 reservoir beneath a seismically active volcano

Mammoth Mountain is a seismically active volcano 200 000 to 50 000 years old, situated on the southwestern rim of Long Valley caldera, California. Since 1989 it has shown evidence of unrest in the form of earthquake swarms (Hill et al. 1990), volcanic 'long-period' earthquakes (Pitt & Hill 1994), increased output of magmatic ³He (Sorey et al. 1993) and the emission of about 500 tonnes day ⁻¹ of CO₂ (Farrar et al. 1995; Hill 1996; M. Sorey, personal communication, 1997), which has killed trees and poses a threat to human safety. Local-earthquake tomography shows that in mid-1989 areas of subsequent tree-kill were underlain by extensive regions where the ratio of the compressional and shear elastic-wave speeds V_P/V_S was about 9 per cent lower than in the surrounding rocks. Theory (Mavko & Mukerji 1995), experiment (Ito, DeVilbiss & Nur 1979), and experience at other geothermal/volcanic areas (Julian et al. 1996) and at petroleum reservoirs (Harris et al. 1996) indicate that V_P/V_S is sensitive to pore-fluid compressibility, through its effect on V_P . The observed V_P/V_S anomaly is probably caused directly by CO₂, and seismic V_P/V_S tomography is thus a promising tool for monitoring gas concentration and movement in volcanoes, which may in turn be related to volcanic activity.     

Coupled normal-mode sensitivity to inner-core shear velocity and attenuation

We have studied the response of normal modes to perturbations in inner-core shear velocity and attenuation, using fully coupled mode synthetics. Our results indicate that (i) mode pairs   _n S _l – _{n ±1} S _l   are strongly coupled by anelasticity, (ii) this coupling causes shear velocity perturbations to strongly affect the Q values of modes through exchange of inner-core characteristics, (iii) there is no evidence for a weakly attenuating inner core in shear, and (iv) the discrepancy between attenuation models returned from normal modes and body waves is small. These results suggest that inversions for inner-core attenuation and shear velocity should be performed simultaneously and should take account of the strong cross-coupling due to attenuation.     

Focal mechanisms of recent earthquakes in the Southern Korean Peninsula

We evaluate the stress field in and around the southern Korean Peninsula with focal mechanism solutions, using the data collected from 71 earthquakes ( M_L = 1.9–5.2) between 1999 and 2004. For this, the hypocentres were relocated and well-constrained fault plane solutions were obtained from the data set of 1270 clear P -wave polarities and 46 SH / P amplitude ratios. The focal mechanism solutions indicate that the prevailing faulting types in South Korea are strike-slip-dominant-oblique-slip faultings with minor reverse-slip component. The maximum principal stresses (σ₁) estimated from fault-slip inversion analysis of the focal mechanism solutions show a similar orientation with E–W trend (269°–275°) and low-angle plunge (10°–25°) for all tectonic provinces in South Korea, consistent with the E–W trending maximum horizontal stress (σ_Hmax) of the Amurian microplate reported from in situ stress measurements and earthquake focal mechanisms. The directions of the intermediate (σ₂) and minimum (σ₃) principal stresses of the Gyeongsang Basin are, however, about 90 deg off from those of the other tectonic provinces on a common σ₂–σ₃ plane, suggesting a permutation of σ₂ and σ₃. Our results incorporated with those from the kinematic studies of the Quaternary faults imply that NNW- to NE-striking faults (dextral strike-slip or oblique-slip with a reverse-slip component) are highly likely to generate earthquakes in South Korea.     

Palaeomagnetism of the continental sector of the Cameroon Volcanic Line, West Africa

Measurement of samples from 154 sites in the continental sector of the Cameroon Volcanic Line yielded six palaeomagnetic poles, at 243.6°E, 84.6°N, α ₉₅ = 6.8°; 224.3°E, 81.2°N, α ₉₅ = 8.4°; 176.1°E, 82.0°N, α ₉₅ = 8.5°; 164.3°E, 86.4°N, α ₉₅ = 3.4°; 169.4°E, 82.6°N, α ₉₅ = 4.6° and 174.7°E, 72.8°N, α ₉₅ = 9.5°, belonging to rocks which have been dated by the K–Ar method at 0.4–0.9  Ma, 2.6  Ma, 6.5–11  Ma, 12–17  Ma, 20–24  Ma and 28–31  Ma, respectively. The results are in general agreement with other palaeomagnetic poles from Oligocene to Recent formations in Africa.
  The first three poles for rocks formed between 0.4 and 11  Ma are not significantly different from the present geographical pole. Together with other African poles for the same period, this suggests that the African continent has moved very little relative to the pole since 11  Ma. The other three poles for rocks dated between 12 and 31  Ma are significantly different from the present geographical pole, showing a 5° polar deviation from the present pole in the Miocene and 13° in the Middle Oligocene.     

Observation of Coriolis coupled modes below 1 mHz

We present observations of spectral energy at toroidal mode frequencies in vertical seismic recordings of the 1998 Balleny Islands earthquake. Since toroidal modes on a spherically symmetric, nonrotating Earth have horizontally polarized particle motion these observations call for an explanation. We first rule out local and instrumental effects as being responsible for the verticalcomponent signal of the toroidal modes ₀ T ₃ (0.59 mHz) and ₀ T ₄ (0.77 mHz). The global effects that we consider are general heterogeneous mantle structure, ellipticity of figure and rotation. We find that rotation through Coriolis coupling of loworder spheroidal and toroidal oscillations is the dominant mechanism.     

Tomographic inversion using ℓ1-norm regularization of wavelet coefficients

We propose the use of ℓ₁ regularization in a wavelet basis for the solution of linearized seismic tomography problems   A m = d   , allowing for the possibility of sharp discontinuities superimposed on a smoothly varying background. An iterative method is used to find a sparse solution m that contains no more fine-scale structure than is necessary to fit the data d to within its assigned errors.     

Thermochronometry and palaeomagnetism of the Archaean Nelshoogte Pluton, South Africa

²⁰⁷Pb/²⁰⁶Pb single-grain zircon, ⁴⁰Ar/³⁹Ar single-grain hornblende and biotite, and ⁴⁰Ar/³⁹Ar bulk-sample muscovite and biotite ages from the Nelshoogte trondhjemite pluton located in eastern Transvaal, South Africa, show that this granitoid had a protracted thermal history spanning 3213±4  Ma to about 3000  Ma. Whole-rock ⁴⁰Ar/³⁹Ar ages from cross-cutting dolerite dykes indicate that these were intruded at about 1900  Ma. There is no evidence of this or other, later events significantly affecting the argon systematics of the minerals from the pluton dated by the ⁴⁰Ar/³⁹Ar method.
  The pluton has a well-defined palaeomagnetic pole which is dated at 3179±18 (2 σ ) Ma by ⁴⁰Ar/³⁹Ar dating of hornblende. This pole (18°N, 310°E, A ₉₅=9°) yields a palaeolatitude of 0°, significantly different from other Archaean poles from the Kaapvaal Craton. The palaeolatitude difference implies that there was significant apparent polar wander during the Archaean. A second, overprinting magnetization seen in the pluton is also seen in the lower-Proterozoic dolerite dykes, and is consistent with other lower-Proterozoic (2150–1950  Ma) poles for southern Africa.     

Depositional and tectonic evolution of a supradetachment basin: 40Ar/39Ar geochronology of the Nova Formation, Panamint Range, California

The Nova Basin contains an upper Miocene to Pliocene supradetachment sedimentary succession that records the unroofing of the Panamint metamorphic core complex, west of Death Valley, California. Basin stratigraphy reflects the evolution of sedimentation processes from landslide emplacement during basin initiation to the development of alluvial fans composed of reworked, uplifted sections of the basin fill. ⁴⁰Ar/³⁹Ar geochronology of volcanic units in middle and lower parts of the sequence provide age control on the tectonic and depositional evolution of the basin and, more generally, insights regarding the rate of change of depositional environments in supradetachment basins. Our work, along with earlier research, indicate basin deposition from 11.38 Ma to 3.35 Ma. The data imply sedimentation rates, uncorrected for compaction, of ~100 m Myr⁻¹ in the lower, high-energy part to ~1000 m Myr⁻¹ in the middle part characterized by debris-flow fan deposition. The observed variation in sediment flux rate during basin evolution suggests that supradetachment basins have complex depositional histories involving rapid transitions in both the style and rate of sedimentation.     

Surface motion of a fluid planet induced by impacts

In order to approximate the free-surface motion of an Earth-sized planet subjected to a giant impact, we have described the excitation of body and surface waves in a spherical compressible fluid planet without gravity or intrinsic material attenuation for a buried explosion source. Using the mode summation method, we obtained an analytical solution for the surface motion of the fluid planet in terms of an infinite series involving the products of spherical Bessel functions and Legendre polynomials. We established a closed form expression for the mode summation excitation coefficient for a spherical buried explosion source, and then calculated the surface motion for different spherical explosion source radii a (for cases of   a / R = 0.001  to 0.035, R is the radius of the Earth) We also studied the effect of placing the explosion source at different radii r ₀ (for cases of   r ₀/ R = 0.90  to 0.96) from the centre of the planet. The amplitude of the quasi-surface waves depends substantially on a / R , and slightly on   r ₀/ R   . For example, in our base-line case,   a / R = 0.03, r ₀/ R = 0.96  , the free-surface velocity above the source is 0.26 c , whereas antipodal to the source, the peak free surface velocity is 0.19 c . Here c is the acoustic velocity of the fluid planet. These results can then be applied to studies of atmosphere erosion via blow-off caused by asteroid impacts.     

ScSp observed on North Island, New Zealand: implications for subducting plate structure

Seismic phase conversions provide important constraints on the layered nature of subduction zone structures. Recordings from digital stations in North Island, New Zealand, have been examined for converted ScS ‐to‐ p ( ScSp ) arrivals from deep (>150 km) Tonga–Kermadec earthquakes to image layering in the underlying Hikurangi subduction zone. Consistent P ‐wave energy prior to ScS has been identified from stations in eastern and southern North Island, where the subducted plate interface is at a depth of between 15 and 30 km. Two ScS precursors are observed. Ray tracing indicates that the initial precursor ( ScSp ₁) corresponds to conversion from the base of an 11–14 km thick subducting Pacific crust. The second precursor is interpreted as a conversion from the top of the subducting plate. The amplitude ratio, ScSp ₁: ScS , increases from 0.10 to 0.19 from northern to southern North Island. This is within the range expected from a simple first‐order velocity discontinuity at an oceanic Moho. A 1–2 km thick layer of low‐velocity sediment at the top of the subducting plate is required to explain the remaining ScSp waveform. Our results imply that the abnormally thick Hikurangi–Chatham Plateau has been subducting beneath New Zealand for at least 2.9 Myr, thus explaining the high uplift rates observed across eastern North Island.     

U–Pb zircon and 40Ar–39Ar K-feldspar dating of syn-sedimentary volcanism of the Neoproterozoic Maricá Formation: constraining the age of foreland basin inception and inversion in the Camaquã Basin of southern Brazil

New U–Pb zircon and ⁴⁰Ar–³⁹Ar K-feldspar data are presented for syn-sedimentary volcanogenic rocks from the Neoproterozoic Maricá Formation, located in the southern Brazilian shield. Seven (of nine) U–Pb sensitive high-resolution ion microprobe analyses of zircons from pyroclastic cobbles yield an age of 630.2±3.4 Ma (2σ), interpreted as the age of syn-sedimentary volcanism, and thus of the deposition itself. This result indicates that the Maricá Formation was deposited during the main collisional phase (640–620 Ma) of the Brasiliano II orogenic system, probably as a forebulge or back-bulge, craton-derived foreland succession. Thus, this unit is possibly correlative of younger portions of the Porongos, Brusque, Passo Feio, Abapã (Itaiacoca) and Lavalleja (Fuente del Puma) metamorphic complexes. Well-defined, step-heating ⁴⁰Ar–³⁹Ar K-feldspar plateau ages obtained from volcanogenic beds and pyroclastic cobbles of the lower and upper successions of the Maricá Formation yielded 507.3±1.8 Ma and 506.7±1.4 Ma (2σ), respectively. These data are interpreted to reflect total isotopic resetting during deep burial and thermal effects related to magmatic events. Late Middle Cambrian cooling below ca . 200 °C, probably related to uplift, is tentatively associated with intraplate effects of the Rio Doce and/or Pampean orogenies (Brasiliano III system). In the southern Brazilian shield, these intraplate stresses are possibly related to the dominantly extensional opening of a rift or a pull-apart basin, where sedimentary rocks of the Camaquã Group (Santa Bárbara and Guaritas Formations) accumulated.