			Home About us Contact

Seismic Tomography (seismic + tomography)

Distribution by Scientific Domains

Earth and Environmental Science	90%

Selected Abstracts

Yardsticks for industrial tomography

GEOPHYSICAL PROSPECTING, Issue 4 2008
A. Vesnaver
ABSTRACT Seismic tomography has been developed and applied for decades in seismological applications and for basic research purposes. During the last decade, large-scale applications in the oil and gas industry became standard as tomostatics and velocity modelling for pre-stack depth migration. In this paper, I take a snapshot of some current industrial applications, quantifying practical aspects by yardsticks such as data and model size and I try to draw a road map for the current decade. [source]

2-D/3-D multiply transmitted, converted and reflected arrivals in complex layered media with the modified shortest path method

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2009
Chao-Ying Bai
SUMMARY Grid-cell based schemes for tracing seismic arrivals, such as the finite difference eikonal equation solver or the shortest path method (SPM), are conventionally confined to locating first arrivals only. However, later arrivals are numerous and sometimes of greater amplitude than the first arrivals, making them valuable information, with the potential to be used for precise earthquake location, high-resolution seismic tomography, real-time automatic onset picking and identification of multiple events on seismic exploration data. The purpose of this study is to introduce a modified SPM (MSPM) for tracking multiple arrivals comprising any kind of combination of transmissions, conversions and reflections in complex 2-D/3-D layered media. A practical approach known as the multistage scheme is incorporated into the MSPM to propagate seismic wave fronts from one interface (or subsurface structure for 3-D application) to the next. By treating each layer that the wave front enters as an independent computational domain, one obtains a transmitted and/or converted branch of later arrivals by reinitializing it in the adjacent layer, and a reflected and/or converted branch of later arrivals by reinitializing it in the incident layer. A simple local grid refinement scheme at the layer interface is used to maintain the same accuracy as in the one-stage MSPM application in tracing first arrivals. Benchmark tests against the multistage fast marching method are undertaken to assess the solution accuracy and the computational efficiency. Several examples are presented that demonstrate the viability of the multistage MSPM in highly complex layered media. Even in the presence of velocity variations, such as the Marmousi model, or interfaces exhibiting a relatively high curvature, later arrivals composed of any combination of the transmitted, converted and reflected events are tracked accurately. This is because the multistage MSPM retains the desirable properties of a single-stage MSPM: high computational efficiency and a high accuracy compared with the multistage FMM scheme. [source]

Effects of mode coupling and location of rotational axis on glacial induced rotational deformation in a laterally heterogeneous viscoelastic earth

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 2 2006
Patrick Wu
SUMMARY Lateral viscosity variations introduce two problems in the study of glacial induced rotational deformation. One of them is mode coupling. The other arises because centrifugal deformations are dependent on the location of the rotational axis relative to the lateral heterogeneities in the earth. The effects of mode coupling and displacement of the rotational axis on centrifugal induced potential perturbations are studied on a layered earth with lateral viscosity variations inferred from seismic tomography. The modified coupled Laplace-finite-element method is used to calculate the centrifugal deformations. It is found that mode coupling is weak because the lateral viscosity variation in the deep mantle is relatively small. The effect of displacing the rotational axis is to shift the spectral amplitude from degree 2 and order 0 harmonic to other orders within degree 2 as is required for changing the orientation of the applied centrifugal potential, but again the leakage of spectral amplitude to other harmonics are small. This implies that rotational motion in an earth with lateral viscosity variations can be approximated by the calculation of the corresponding laterally homogeneous earth. [source]

Near-lithostatic pore pressure at seismogenic depths: a thermoporoelastic model

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 3 2006
Francesca Zencher
SUMMARY A model is presented for pore pressure migration through a transition layer separating a meteoric aquifer at hydrostatic pressure from a deeper reservoir at lithostatic pressure. This configuration is thought to be pertinent to the South Iceland seismic zone (SISZ) and to other tectonically active regions of recent volcanism, where volatiles are continuously released by ascending magma below the brittle,ductile transition. Poroelastic parameters are computed for basaltic rock. The model is 1-D, the fluid viscosity is temperature dependent and rock permeability is assumed to be pressure dependent according to a dislocation model of a fractured medium. Environment conditions are considered, pertinent to basalt saturated with water at shallow depth (case I) and at mid-crustal depth (case II). If the intrinsic permeability of the rock is high, no significant effects are observed in the pressure field but advective heat transfer shifts the brittle,ductile transition to shallower depths. If the intrinsic permeability is low, the pressure-dependent permeability can propagate near-lithostatic pore pressures throughout most of the transition layer, while the temperature is practically unaffected by advective contributions so that the rock in the transition layer remains in brittle condition. Geometrical parameters characterizing the fracture distribution are important in determining the effective permeability: in particular, if an interconnected system of fractures develops within the transition layer, the effective permeability may increase by several orders of magnitude and near-lithostatic pore pressure propagates upwards. These modelling results have important bearings on our understanding of seismogenic processes in geothermal areas and are consistent with several geophysical observations in the SISZ, in connection with the two 2000 June M= 6.5 earthquakes, including: (i) fluid pressure pulses in deep wells, (ii) low resistivity at the base of the seismogenic layer, (iii) low VP/VS ratio and time-dependent seismic tomography, (iv) heterogeneity of focal mechanisms, (v) shear wave splitting, (vi) high b -value of deep foreshocks, (vii) triggered seismicity and (viii) Radon anomalies. [source]

Analytical and 3-D numerical modelling of Mt. Etna (Italy) volcano inflation

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 2 2005
A. Bonaccorso
SUMMARY Since 1993, geodetic data obtained by different techniques (GPS, EDM, SAR, levelling) have detected a consistent inflation of the Mt. Etna volcano. The inflation, culminating with the 1998,2001 strong explosive activity from summit craters and recent 2001 and 2002 flank eruptions, is interpreted in terms of magma ascent and refilling of the volcanic plumbing system and reservoirs. We have modelled the 1993,1997 EDM and GPS data by 3-D pressurized sources to infer the position and dimension of the magma reservoir. We have performed analytical inversions of the observed deformation using both spheroidal and ellipsoidal sources embedded in a homogeneous elastic half-space and by applying different inversion methods. Solutions for these types of sources show evidence of a vertically elongated magma reservoir located 6 km beneath the summit craters. The maximum elevation of topography is comparable to such depth and strong heterogeneities are inferred from seismic tomography; in order to assess their importance, further 3-D numerical models, employing source parameters extracted from analytical models, have been developed using the finite-element technique. The deformation predicted by all the models considered shows a general agreement with the 1993,1997 data, suggesting the primary role of a pressure source, while the complexities of the medium play a minor role under elastic conditions. However, major discrepancies between data and models are located in the SE sector, suggesting that sliding along potential detachment surfaces may contribute to amplify deformation during the inflation. For the first time realistic features of Mt. Etna are studied by a 3-D numerical model characterized by the topography and lateral variations of elastic structure, providing a framework for a deeper insight into the relationships between internal sources and tectonic structures. [source]

Upwelling plumes, superswells and true polar wander

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 3 2004
Marianne Greff-Lefftz
SUMMARY The geological evolution of the rotational axis of the Earth is most likely controlled by internal mass redistribution within the mantle. Palaeomagnetic observations suggest that it is episodic in nature, with periods of quasi-standstill alternating with periods of faster wander. Here, we investigate two models for the influence of mantle plumes that vary at different spatial wavelengths on the time variations of the rotational axis (true polar wander; TPW). In the first model, we represent an upwelling plume as a sphere whose radius varies as a function of the flux of material in the conduit and that traverses the mantle at the Stokes velocity. Such a plume produces very little wander of the rotational axis. We then study the effects of two superswells that mimic the ones observed with seismic tomography and conclude that a doming regime within the mantle involves significant polar wander. Some of the features of this TPW that are directly linked to the periodicity of doming are reminiscent of observed phases of slow and fast TPW, with similar peak velocities. [source]

Simulated geomagnetic reversals and preferred virtual geomagnetic pole paths

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 3 2004
C. Kutzner
SUMMARY The question of whether virtual geomagnetic poles (VGPs) recorded during reversals and excursions show a longitudinal preference is a controversial one amongst palaeomagnetists. One possible mechanism for such VGP clustering is the heterogeneity of heat flux at the core,mantle boundary (CMB). We use 3-D convection-driven numerical dynamo models with imposed non-uniform CMB heat flow that show stochastic reversals of the dipole field. We calculate transitional VGPs for a large number of token sites at the Earth's surface. In a model with a simple heat flux variation given by a Y22 harmonic, the VGP density maps for individual reversals differ substantially from each other, but the VGPs have a tendency to fall around a longitude of high heat flow. The mean VGP density for many reversals and excursions shows a statistically significant correlation with the heat flow. In a model with an imposed heat flux pattern derived from seismic tomography we find maxima of the mean VGP density at American and East Asian longitudes, roughly consistent with the VGP paths seen in several palaeomagnetic studies. We find that low-latitude regions of high heat flow are centres of magnetic activity where intense magnetic flux bundles are generated. They contribute to the equatorial dipole component and bias its orientation in longitude. During reversals the equatorial dipole part is not necessarily dominant at the Earth's surface, but is strong enough to explain the longitudinal preference of VGPs as seen from different sites. [source]

Cretaceous,Tertiary geodynamics: a North Atlantic exercise

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 3 2001
Trond H. Torsvik
Summary New reconstructions are presented for the Cretaceous,Early Tertiary North Atlantic using a combination of palaeomagnetic, hotspot and magnetic anomaly data. We utilize these reconstructions in an analysis of previously described misfits between the North Atlantic Plate elements at successive intervals during this time period. We are able to achieve reasonable overlap between the hotspot and palaeomagnetic reconstructions between 40 and 95 Ma and thus are able to support the idea that the Indo,Atlantic hotspots are relatively stationary. Small, but systematic discrepancies for this time interval can readily be modelled with a long-term, octopole non-dipole field contribution (G3 = g30/g10 = 0.08). However, hotspot and palaeomagnetic reconstructions for the Early Cretaceous North Atlantic show substantial differences that cannot be explained by constant, non-dipole fields and we favour an explanation for these discrepancies in terms of true polar wander (TPW) triggered by mantle instabilities between 125 and 95 Ma; this constitutes the only identifiable event of significant TPW since the Early Cretaceous. Taken in the context of available geochronological and geological data and seismic tomography from the region, the 95,40 Ma reconstructions and their time-consequent geological products are interpreted in terms of specific conditions of mantle-crust coupling and global plate motions/tectonic activity. Highlights from these reconstructions show uniform NE movement of the coupled North American, Greenland and Eurasian plates from 95 to 80 Ma; a marked cusp in the paths for all three elements at 80 Ma where the three plates simultaneously change direction and follow a uniform NW-directed motion until c. 20 Ma when Eurasia diverges NE, away from the still-NW-moving Greenland and North American elements. Positioning of the Iceland plume beneath the spreading-ridge at 20 Ma may have increased upwelling below the ridge, increased the ridge-push, and caused a NE shift in the absolute direction of Eurasia. [source]

Convection in the Earth's core driven by lateral variations in the core,mantle boundary heat flux

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 2 2000
Steven John Gibbons
Summary Moving core fluid maintains an isothermal core,mantle boundary (CMB), so lateral variations in the CMB heat flow result from mantle convection. Such variations will drive thermal winds, even if the top of the core is stably stratified. These flows may contribute to the magnetic secular variation and are investigated here using a simple, non-magnetic numerical model of the core. The results depend on the equatorial symmetry of the boundary heat flux variation. Large-scale equatorially symmetric (ES) heat flux variations at the outer surface of a rapidly rotating spherical shell drive deeply penetrating flows that are strongly suppressed in stratified fluid. Smaller-scale ES heat flux variations drive flows less dominated by rotation and so less inhibited by stratification. Equatorially anti-symmetric flux variations drive flows an order of magnitude less energetic than those driven by ES patterns but, due to the nature of the Coriolis force, are less suppressed by stratification. The response of the rotating core fluid to a general CMB heat flow pattern will then depend strongly on the subadiabatic temperature profile. Imposing a lateral heat flux variation linearly related to a model of seismic tomography in the lowermost mantle drives flow in a density stratified fluid that reproduces some features found in flows inverted from geomagnetic data. [source]

Geophysical surveys designed to delineate the altitudinal limit of mountain permafrost: an example from Jotunheimen, Norway

PERMAFROST AND PERIGLACIAL PROCESSES, Issue 3 2004
Christian Hauck
Abstract Three geophysical methods have been applied to delineate the altitudinal limit of permafrost at Juvvasshøe/Jotunheimen in southern Norway. By using each method in a complementary way according to its applicability, the permafrost distribution could be analysed on both large and small scales. In addition, temperature-based methods such as the BTS method (bottom temperature of snow cover) were used to validate the results. On a large scale, electromagnetic induction profiling with the EM-31 was used to detect the location of the transition area between frozen and unfrozen ground. Within this area direct-current (DC) resistivity and refraction seismic tomography were also applied to further characterize and visualize the permafrost transition. The large scale surveys revealed a sharp increase in conductivity, indicating an increase in unfrozen water content, near 1400,m a.s.l., with a possible transition zone between 1500,m a.s.l. and 1380,m a.s.l. The small scale investigation delineated the altitudinal limit of permafrost between 1470,m a.s.l. and 1410,m a.s.l., with sporadic ground ice occurrences and temperatures near the freezing point. Copyright © 2004 John Wiley & Sons, Ltd. [source]