Surface wave ray tracing and azimuthal anisotropy: a generalized spherical harmonic approach     

Lapo Boschi  John H. Woodhouse 《Geophysical Journal International》2006,164(3):569-578

  相似文献   

Global multiresolution models of surface wave propagation: comparing equivalently regularized Born and ray theoretical solutions     

Lapo Boschi 《Geophysical Journal International》2006,167(1):238-252

I invert a large set of teleseismic phase-anomaly observations, to derive tomographic maps of fundamental-mode surface wave phase velocity, first via ray theory, then accounting for finite-frequency effects through scattering theory, in the far-field approximation and neglecting mode coupling. I make use of a multiple-resolution pixel parametrization which, in the assumption of sufficient data coverage, should be adequate to represent strongly oscillatory Fréchet kernels. The parametrization is finer over North America, a region particularly well covered by the data. For each surface-wave mode where phase-anomaly observations are available, I derive a wide spectrum of plausible, differently damped solutions; I then conduct a trade-off analysis, and select as optimal solution model the one associated with the point of maximum curvature on the trade-off curve. I repeat this exercise in both theoretical frameworks, to find that selected scattering and ray theoretical phase-velocity maps are coincident in pattern, and differ only slightly in amplitude.  相似文献   

Surface wave tomography: global membrane waves and adjoint methods     

D. Peter  C. Tape  L. Boschi  J. H. Woodhouse 《Geophysical Journal International》2007,171(3):1098-1117

We implement the wave equation on a spherical membrane, with a finite-difference algorithm that accounts for finite-frequency effects in the smooth-Earth approximation, and use the resulting 'membrane waves' as an analogue for surface wave propagation in the Earth. In this formulation, we derive fully numerical 2-D sensitivity kernels for phase anomaly measurements, and employ them in a preliminary tomographic application. To speed up the computation of kernels, so that it is practical to formulate the inverse problem also with respect to a laterally heterogeneous starting model, we calculate them via the adjoint method, based on backpropagation, and parallelize our software on a Linux cluster. Our method is a step forward from ray theory, as it surpasses the inherent infinite-frequency approximation. It differs from analytical Born theory in that it does not involve a far-field approximation, and accounts, in principle, for non-linear effects like multiple scattering and wave front healing. It is much cheaper than the more accurate, fully 3-D numerical solution of the Earth's equations of motion, which has not yet been applied to large-scale tomography. Our tomographic results and trade-off analysis are compatible with those found in the ray- and analytical-Born-theory approaches.  相似文献   

Strain measurements in the Ancona area     

P. Baldi  D. Postpischl  M. Unguendoli  E. Boschi 《Geophysical Journal International》1977,51(1):253-258

The purpose of this short note is to describe a plano-altimetric network recently installed in the Ancona (Italy) seismic area. The results of the first two series of surveys performed (1975 September, 1976 September) are also shown and discussed.  相似文献   

On the interpretation of slow ground deformation precursory to the 1976 Friuli earthquake     

Michele Dragoni  Maurizio Bonafede  Enzo Boschi 《Pure and Applied Geophysics》1984,122(6):781-792

During three years preceding the 1976 Friuli earthquake, a continuous southward ground tilt was recorded by a tiltmeter placed near Tolmezzo, 15 km north-west of the epicentre of the impending earthquake. The cumulative ground tilt amounted to as much as 3 minutes of arc. Since the tiltmeter was placed in the proximity of an active fault, such a tilt can be explained if the fault slipped aseismically on its shallower section during the same three year period. Aseismic slip on the fault might have been caused by the same mechanism which concentrated stress in the region and eventually produced the 1976 earthquake.  相似文献   

Downslope flow models of a Bingham liquid: Implications for lava flows     

Michele Dragoni  Maurizio Bonafede  Enzo Boschi 《Journal of Volcanology and Geothermal Research》1986,30(3-4)

It is widely recognized that lavas behave as Bingham liquids, which are characterized by a yield stress σ and a plastic viscosity η. We consider two models describing downslope flows of a Bingham liquid with different aspect ratios A (= flow height/flow width): model 1 with A 1 and model 2 with A ≈ 1. Sufficiently uphill with respect to the front, such flows can be considered as laminar and locally isothermal. For both models, we obtain analytically the steady-state solution of the Navier-Stokes equations and the constitutive equation for a Bingham liquid. We study the flow height and velocity as functions of flow rate, rheological parameters and ground slope. It is found that such flows remain in the Newtonian regime at low yield stresses (σ 10³dyne/cm²), but the transition to the Bingham regime also depends on flow rate and occurs at higher values of σ for higher flow rates: for instance, a high aspect ratio flow (model 2) is still very close to the Newtonian regime at σ = 10⁴ dyne/cm², if the flow rate is greater than 10⁵ g/s. In the Bingham regime, flow heights are generally greater and flow velocities are smaller than in the Newtonian regime; moreover, flow heights are independent of flow rate, so that a change in flow rate results exclusively in a velocity change. After assuming a specific temperature dependence of σ and η between the solidus and the liquidus temperatures of an ideal Bingham liquid (1000°C and 1200 °C respectively), flow heights and velocities are examined as functions of temperature along the flow. Several effects observed in lava flows are predicted by these models and allow a more quantitative insight into the behaviour of lava flows.  相似文献   

Applying the relative hypocentre location approach: where was the 1980 November 23 Irpinia earthquake?     

Domenico Giardini  Alberto Basili  Enzo Boschi 《Geophysical Journal International》1996,127(3):605-615

The purpose of this work is to evaluate under what conditions it is feasible and with what accuracy it is possible to locate the nucleation point of a large earthquake, given the availability of aftershocks located with high precision by the deployment of a local network. We experiment with several approaches and apply them lo the location of the epicentre of the 1980 November 23 Irpinia earthquake ( M ^w= 6.9).
First we use local P ^g phases selected to optimize the azimuthal coverage, obtaining a well-constrained location with a small statistical error, which typically underestimates the true hypocentre uncertainty.
We then exploit the relative location technique, obtaining stable, almost coincident solutions under three conditions: (1) using multiple independent master events to derive an average epicentre; (2) fitting simultaneously the larger data set for all available master events, using a forward approach; (3) conducting an a priori evaluation of the statistics of station and master events to separate model uncertainties and improve the statistical accuracy of the relative locations. Moreover, only by introducing station statistics can we achieve the desired accuracy of ≅ 1 km in constraining the rupture nucleation point of this large earthquake, and we show that the application of the relative location technique to uncleaned, unweighted data for a single master event provides only a crude epicentre with a confidence ellipse deceivingly smaller than the true hypocentre uncertainty.
The revised epicentre for the 1980 November 23 Irpinia earthquake (48.803 °N-15.302°E) validates the class of multidisciplinary reconstructions of the source process such as the model of Valensise et al. (1989), based on the hypocentre of Westaway & Jackson (1987), and is shifted by almost 13 km to the NW of the epicentre recently proposed by Westaway (1992).  相似文献   

Magma convection and mixing dynamics as a source of Ultra-Long-Period oscillations     

Antonella Longo  Paolo Papale  Melissa Vassalli  Gilberto Saccorotti  Chiara P. Montagna  Andrea Cassioli  Salvatore Giudice  Enzo Boschi 《Bulletin of Volcanology》2012,74(4):873-880

Many volcanic eruptions are shortly preceded by injection of new magma into a pre-existing, shallow (<10 km) magma chamber, causing convection and mixing between the incoming and resident magmas. These processes may trigger dyke propagation and further magma rise, inducing long-term (days to months) volcano deformation, seismic swarms, gravity anomalies, and changes in the composition of volcanic plumes and fumaroles, eventually culminating in an eruption. Although new magma injection into shallow magma chambers can lead to hazardous event, such injection is still not systematically detected and recognized. Here, we present the results of numerical simulations of magma convection and mixing in geometrically complex magmatic systems, and describe the multiparametric dynamics associated with buoyant magma injection. Our results reveal unexpected pressure trends and pressure oscillations in the Ultra-Long-Period (ULP) range of minutes, related to the generation of discrete plumes of rising magma. Very long pressure oscillation wavelengths translate into comparably ULP ground displacements with amplitudes of order 10⁻⁴–10⁻² m. Thus, new magma injection into magma chambers beneath volcanoes can be revealed by ULP ground displacement measured at the surface.  相似文献   

Gradual Fault Weakening with Seismic Slip: Inferences from the Seismic Sequences of L’Aquila, 2009, and Northridge, 1994   总被引：1，自引：0，他引：1  

Luca Malagnini  Irene Munafo’  Massimo Cocco  Stefan Nielsen  Kevin Mayeda  Enzo Boschi 《Pure and Applied Geophysics》2014,171(10):2709-2730

We estimate seismological fracture energies from two subsets of events selected from the seismic sequences of L’Aquila (2009), and Northridge (1994): 57 and 16 selected events, respectively, including the main shocks. Following Abercrombie and Rice (Geophys J Int 162: 406–424, 2005), we postulate that fracture energy (G) represents the post-failure integral of the dynamic weakening curve, which is described by the evolution of shear traction as a function of slip. Following a direct-wave approach, we compute mainshock-/aftershock-source spectral ratios, and analyze them using the approach proposed by Malagnini et al. (Pure Appl. Geophys., this issue, 2014) to infer corner frequencies and seismic moment. Our estimates of source parameters (including fracture energies) are based on best-fit grid-searches performed over empirical source spectral ratios. We quantify the source scaling of spectra from small and large earthquakes by using the MDAC formulation of Walter and Taylor (A revised Magnitude and Distance Amplitude Correction (MDAC2) procedure for regional seismic discriminants, 2001). The source parameters presented in this paper must be considered as point-source estimates representing averages calculated over specific ruptured portions of the fault area. In order to constrain the scaling of fracture energy with coseismic slip, we investigate two different slip-weakening functions to model the shear traction as a function of slip: (i) a power law, as suggested by Abercrombie and Rice (Geophys J Int 162: 406–424, 2005), and (ii) an exponential decay. Our results show that the exponential decay of stress on the fault allows a good fit between measured and predicted fracture energies, both for the main events and for their aftershocks, regardless of the significant differences in the energy budgets between the large (main) and small earthquakes (aftershocks). Using the power-law slip-weakening function would lead us to a very different situation: in our two investigated sequences, if the aftershock scaling is extrapolated to events with large slips, a power law (a la Abercrombie and Rice) would predict unrealistically large stress drops for large, main earthquakes. We conclude that the exponential stress evolution law has the advantage of avoiding unrealistic stress drops and unbounded fracture energies at large slip values, while still describing the abrupt shear-stress degradation observed in high-velocity laboratory experiments (e.g., Di Toro et al., Fault lubrication during earthquakes, Nature 2011).  相似文献   

Lava effusion — A slow fuse for paroxysms at Stromboli volcano?     

S. Calvari  L. Spampinato  A. Bonaccorso  C. Oppenheimer  E. Rivalta  E. Boschi 《Earth and Planetary Science Letters》2011,301(1-2):317-323

The 2007 effusive eruption of Stromboli followed a similar pattern to the previous 2002–2003 episode. In both cases, magma ascent led to breaching of the uppermost part of the conduit forming an eruptive fissure that discharged lava down the Sciara del Fuoco depression. Both eruptions also displayed a ‘paroxysmal’ explosive event during lava flow output. From daily effusion rate measurements retrieved from helicopter- and satellite-based infrared imaging, we deduce that the cumulative volume of lava erupted before each of the two paroxysms was similar. Based on this finding, we propose a conceptual model to explain why both paroxysms occurred after this ‘threshold’ cumulative volume of magma was erupted. The gradual decompression of the deep plumbing system induced by magma withdrawal and eruption, drew deeper volatile-rich magma into the conduit, leading to the paroxysms. The proposed model might provide a basis for forecasting paroxysmal explosions during future effusive eruptions of Stromboli.  相似文献