Gravimetric determination of an intrusive complex under the Island of Faial (Azores): some methodological improvements

We present some improvements of a gravity inversion method to determine the geometry of the anomalous bodies for priori density contrasts. The 3-D method is based on an exploratory process applied, not for the global model, but for the steps of a growth approach. The (positive and/or negative) anomalous structure is described by successive aggregation of cells, while its corresponding gravity field remains nearly proportional to the observed one. Moreover, a simple (e.g. linear) regional trend can be simultaneously adjusted. The corresponding program is applied to new gravity data on the volcanic island of Faial (Azores archipelago). The inversion approach shows a subsurface anomalous structure for the island, the main feature being an elongated high-density body. The body is interpreted as a compact sheeted dyke swarm, emplaced along Faial-Pico Fracture Zone, a leaky transform structure that forms the current boundary between Eurasian and African plates in the Azores area. The new results in this paper are (1) a Bouguer gravity anomaly map, (2) several improvements in the inversion process (robust process, optimal balance fitness/model magnitude), (3) a new gravimetric method for estimating the mean terrain density, (4) a 3-D model for subsurface mass anomalies in Faial and (5) some interpretative conclusions about a main intrusive complex detected under the island as a wall-like structure extending from a depth of 0.5 to 6 km b.s.l., with a N100°E trend and corresponding to an early fissural volcanic episode controlled by the regional tectonics.     

Sequential integrated inversion of refraction and wide-angle reflection traveltimes and gravity data for two-dimensional velocity structures

A new algorithm is presented for the integrated 2-D inversion of seismic traveltime and gravity data. The algorithm adopts the 'maximum likelihood' regularization scheme. We construct a 'probability density function' which includes three kinds of information: information derived from gravity measurements; information derived from the seismic traveltime inversion procedure applied to the model; and information on the physical correlation among the density and the velocity parameters. We assume a linear relation between density and velocity, which can be node-dependent; that is, we can choose different relationships for different parts of the velocity–density grid. In addition, our procedure allows us to consider a covariance matrix related to the error propagation in linking density to velocity. We use seismic data to estimate starting velocity values and the position of boundary nodes. Subsequently, the sequential integrated inversion (SII) optimizes the layer velocities and densities for our models. The procedure is applicable, as an additional step, to any type of seismic tomographic inversion.
We illustrate the method by comparing the velocity models recovered from a standard seismic traveltime inversion with those retrieved using our algorithm. The inversion of synthetic data calculated for a 2-D isotropic, laterally inhomogeneous model shows the stability and accuracy of this procedure, demonstrates the improvements to the recovery of true velocity anomalies, and proves that this technique can efficiently overcome some of the limitations of both gravity and seismic traveltime inversions, when they are used independently.
An interpretation of field data from the 1994 Vesuvius test experiment is also presented. At depths down to 4.5 km, the model retrieved after a SII shows a more detailed structure than the model obtained from an interpretation of seismic traveltime only, and yields additional information for a further study of the area.     

A probabilistic approach to the inversion of data from a seismic array and its application to volcanic signals

Array techniques are particularly well‐suited for detecting and quantifying the complex seismic wavefields associated with volcanic activity such as volcanic tremor and long‐period events. The methods based on the analysis of the signal in the frequency domain, or spectral methods, have the main advantages of both resolving closely spaced sources and reducing the necessary computer time, but may severely fail in the analysis of monochromatic, non‐stationary signals. Conversely, the time‐domain methods, based on the maximization of a multichannel coherence estimate, can be applied even for short‐duration pulses. However, for both the time and the frequency domain approaches, an exhaustive definition of the errors associated with the slowness vector estimate is not yet available. Such a definition become crucial once the slowness vector estimates are used to infer source location and extent. In this work we develop a method based on a probabilistic formalism, which allows for a complete definition of the uncertainties associated with the estimate of frequency–slowness power spectra from measurement of the zero‐lag cross‐correlation. The method is based on the estimate of the theoretical frequency–slowness power spectrum, which is expressed as the convolution of the true signal slowness with the array response pattern. Using a Bayesian formalism, the a posteriori probability density function for signal slowness is expressed as the difference, in the least‐squares sense, between the model spectrum and that derived from application of the zero‐lag cross‐correlation technique. The method is tested using synthetic waveforms resembling the quasi‐monochromatic signals often associated with the volcanic activity. Examples of application to data from Stromboli volcano, Italy, allow for the estimate of source location and extent of the explosive activity.     

Towards understanding the lithospheric structure of the southern Chilean subduction zone (36°S–42°S) and its role in the gravity field

The Southern Andes differ significantly from the Central Andes with respect to topography and crustal structures and are, from a geophysical point of view, less well known. In order to provide insight into the along-strike segmentation of the Andean mountain belt, an integrated 3-D density model was developed for the area between latitudes 36°S and 42°S. The model is based on geophysical and geological data acquired in the region over the past years and was constructed using forward density modelling. In general, the gravity field of the South American margin is characterized by a relatively continuous positive anomaly along the coastline and the forearc region, and by negative anomalies along the trench and the volcanic arc. However, in the forearc region of the central part of the study area, located just to the south of the epicentre of the largest ever recorded earthquake (Valdivia, 1960), the trench-parallel positive anomaly is disrupted. The forearc gravity anomaly differences thus allow the study area to be divided into three segments, the northern Arauco-Lonquimay, the middle Valdivia-Liquiñe, and the southern Bahía-Mansa-Osorno segment, which are also evident in geology. In the proposed model, the observed negative gravity anomaly in the middle segment is reproduced by an approximately 5 km greater depth to the top of the slab beneath the forearc region. The depth to the slab is, however, dependent upon the density of the upper plate structures. Therefore, both the upper and lower plates and their interaction have a significant impact on the subduction-zone gravity field.     

Surface coseismic gravity changes caused by dislocations in a 3-D heterogeneous earth

This study describes an examination of surface gravity changes caused by dislocations within a 3-D heterogeneous earth. This new theory is described using six independent dislocations: a vertical strike-slip, two vertical dip-slips perpendicular to each other, and three tensile openings on three perpendicular planes. A combination of the six independent dislocations is useful to compute coseismic gravity changes resulting from an arbitrary seismic source at an arbitrary position. Based on the 3-D lateral inhomogeneous P -wave velocity model, we deduce the 3-D density and S -wave velocity models using the relation of Karato. Finally, numerical computations are performed for a location south of Japan (30°N, 135°E). We calculate the coseismic gravity changes resulting from the six independent dislocations for source depths of 100, 300 and 637 km, respectively. Numerical results show that the maximum 3-D effect varies concomitantly with the dislocation type and the source depth. For seismic problems, the effect of elastic parameter  μ  is dominant.     

3-D magnetotelluric inversion and model validation with gravity data for the investigation of flood basalts and associated volcanic rifted margins

20 magnetotelluric (MT) soundings were collected on the Isle of Skye, Scotland to provide a high-resolution three-dimensional (3-D) electrical resistivity model of a volcanic province within the framework of a project jointly interpreting gravity, seismic, geological and MT data. The full 3-D inversion of the MT data jointly interpreted with gravity data reveals upper crustal structure. The main features of the model are interpreted in conjunction with previous geological mapping and borehole data. Our model extends to 13 km depth, several kilometres below the top of the Lewisian basement. The top of the Lewisian basement is at approximately 7–8 km depth and the topography of its surface was controlled by Precambrian rifting, during which a 4.5 km thick sequence of Torridonian sediments was deposited. The Mesozoic sediments above, which can reach up to 2.2 km thick, have small-scale depocentres and are covered by up to 600 m of Tertiary lava flows. The interpretation of the resistivity model shows that 3-D MT inversion is an appropriate tool to image sedimentary structures beneath extrusive basalt units, where conventional seismic reflection methods may fail.     

Estimating gambling venue catchments for impact assessment using a calibrated gravity model

Gambling using electronic gaming machines (EGMs) has emerged as a significant public health issue. While social impact assessments are required prior to the granting of new gaming machine licenses in Australia, there are a few established techniques for estimating the spatial distribution of a venue’s clientele. To this end, we calibrated a Huff model of gambling venue catchments based on a geocoded postal survey (n = 7040). We investigated the impact of different venue attractiveness measures, distance measures, distance decay functions, levels of spatial aggregation and venue types on model fit and results. We then compared model estimates for different behavioural subgroups. Our calibrated spatial model is a significant improvement on previously published models, increasing R² from 0.23 to 0.64. Venue catchments differ radically in size and intensity. As different population subgroups are attracted to different venues, there is no single best index of venue attractiveness applicable to all subpopulations. The calibrated Huff model represents a useful regulatory tool for predicting the extent and composition of gambling venue catchments. It may assist in decision-making with regard to new license applications and evaluating the impact of health interventions such as mandated reductions in EGM numbers. Our calibrated parameters may be used to improve model accuracy in other jurisdictions.     

P-SH conversions in a flat-layered medium with anisotropy of arbitrary orientation

P-SH conversion is commonly observed in teleseismic P waves, and is often attributed to dipping interfaces beneath the receiver. Our modelling suggests an alternative explanation in terms of flat-layered anisotropy. We use reflectivity techniques to compute three-component synthetic seismograms in a 1-D anisotropic layered medium. For each layer of the medium, we prescribe values of seismic velocities and hexagonally symmetric anisotropy about a common symmetry axis of arbitrary orientation. A compressional wave in an anisotropic velocity structure suffers conversion to both SV -and SH -polarized shear waves, unless the axis of symmetry is everywhere vertical or the wave travels parallel to all symmetry axes. The P-SV conversion forms the basis of the widely used 'receiver function' technique. The P-SH conversion occurs at interfaces where one or both layers are anisotropic. A tilted axis of symmetry and a dipping interface in isotropic media produce similar amplitudes of both direct ( P ) and converted ( Ps ) phases, leaving the backazimuth variation of the P-Ps delay as the main discriminant. Seismic anisotropy with a tilted symmetry axis leads to complex synthetic seismograms in velocity models composed of just a few flat homogeneous layers. It is possible therefore to model observations of P coda with prominent transverse components with relatively simple 1-D velocity structures. Successful retrieval of salient model characteristics appears possible using multiple realizations of a genetic-algorithm (GA) inversion of P coda from several backazimuths. Using GA inversion, we determine that six P coda recorded at station ARU in central Russia are consistent with models that possess strong (> 10 per cent) anisotropy in the top 5 km and between 30 and 43 km depth. The symmetry axes are tilted, and appear aligned with the seismic anisotropy orientation in the mantle under ARU suggested by SKS splitting.