			Home About us Contact

Lateral Variations (lateral + variation)

Distribution by Scientific Domains

Earth and Environmental Science	100%

Selected Abstracts

A transportable mb(Lg) scale for central Europe and implications for low-magnitude Ms,mb discrimination

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2005
Howard J. Patton
SUMMARY We have extended the mb(Lg) method of Nuttli using root-mean-square (rms) amplitudes corrected for noise and a ,,1 dependence for geometrical spreading. Lg waves recorded on the German Regional Seismic Network (GRSN) for earthquakes in south-central Europe were used to develop an mb(Lg) formula requiring a new calibration constant Crms to keep rms mb(Lg) on the same baseline as Nuttli's traditional formula based on 3rd-peak amplitudes. GRSN stations had to be calibrated for site terms and for Lg attenuation. Lateral variations in LgQ appear to be significant across the study area, and a regional Q model consisting of constant- Q partitions north, south and in the central Alps was developed using measurements based on interstation and two-event, single-station methods. When plotted against surface wave estimates of Mw, rms mb(Lg) measurements in central Europe are found to be consistent with Mw,mb(Lg) relationships for north America and southern Asia, thus supporting the transportability of our mb(Lg) formula. Frequency,wavenumber processing of Gräfenberg Array data enabled us to extract Rayleigh waves for small events, and regional Ms were measured using the Marshall and Basham formula. Our Ms,mb(Lg) relationship extends to Ms 2.5 and agrees well with observations in other regions including the western United States. The discrimination potential of Ms,mb(Lg) observations was examined under realistic monitoring conditions, where path corrections were inferred from earthquake data and applied uniformly to natural sources and explosions. Under these conditions, mb(Pn, P) are greater than mb(Lg) for large NTS explosions; however, Ms,mb scaling slopes are steeper for P waves than they are for Lg, and Ms,mb observations for NTS explosions converge near mb 4. Thus, allowing for measurement errors and additional uncertainty in mb(Pn) due to regional bias, there is little difference in the discrimination potential for Pn and Lg waves at small magnitudes. As such, a regional Ms,mb discriminant based on Lg might be preferred owing to the better detectability of Lg waves for small earthquakes. These results need to be confirmed for explosions at other test sites. Compared to teleseismic experience, regional Ms,mb observations extend the discrimination capability to lower magnitudes by at least one Ms unit. [source]

Evolution of an accretionary complex along the north arm of the Island of Sulawesi, Indonesia

ISLAND ARC, Issue 1 2004
Yusuf Surachman Djajadihardja
Abstract Seismic reflections across the accretionary prism of the North Sulawesi provide excellent images of the various structural domains landward of the frontal thrust. The structural domain in the accretionary prism area of the North Sulawesi Trench can be divided into four zones: (i) trench area; (ii) Zone A; (iii) Zone B; and (iv) Zone C. Zone A is an active imbrication zone where a decollement is well imaged. Zone B is dominated by out-of-sequence thrusts and small slope basins. Zone C is structurally high in the forearc basin, overlain by a thick sedimentary sequence. The subducted and accreted sedimentary packages are separated by the decollement. Topography of the oceanic basement is rough, both in the basin and beneath the wedge. The accretionary prism along the North Sulawesi Trench grew because of the collision between eastern Sulawesi and the Bangai,Sula microcontinent along the Sorong Fault in the middle Miocene. This collision produced a large rotation of the north arm of Sulawesi Island. Rotation and northward movement of the north arm of Sulawesi may have resulted in southward subduction and development of the accretionary wedge along North Sulawesi. Lateral variations are wider in the western areas relative to the eastern areas. This is due to greater convergence rates in the western area: 5 km/My for the west and 1.5 km/My for the east. An accretionary prism model indicates that the initiation of growth of the accretionary prism in the North Sulawesi Trench occurred approximately 5 Ma. A comparison between the North Sulawesi accretionary prism and the Nankai accretionary prism of Japan reveals similar internal structures, suggesting similar mechanical processes and structural evolution. [source]

Frontal accretion and thrust wedge evolution under very oblique plate convergence: Fiordland Basin, New Zealand

BASIN RESEARCH, Issue 4 2002
P. M. Barnes
ABSTRACT A thrust wedge with unusual geometry has developed under very oblique (50,60°) convergence between the Pacific and Australian Plates, along the 240-km length of the Fiordland margin, New Zealand. The narrow (25 km-wide) wedge comprises three overlapping components, lying west of the offshore section of the Alpine Fault, and straddles a change of > 30° in the regional strike of the plate boundary. Swath bathymetry, marine seismic reflection profiles, and dated samples together reveal the stratigraphy, structure, and evolution of the wedge and the underthrusting, continental, Caswell High (Australian Plate). Lateral variations in the composition and structure of the accretionary wedge, and the depth of the décollement thrust, result partly from variations in crustal structure and basement relief of the underthrust plate, and from associated variations in the thickness of turbidites available for frontal accretion. In the southern Fiordland Basin the underthrust plate is undergoing flexural uplift and extension, and a thick turbidite section is available for accretion. Along-strike, a structurally elevated portion of the underthrust plate is very obliquely colliding with the central part of the accretionary wedge, the turbidite section available for accretion is condensed, and structural inversion occurs in the underthrust plate. ,Growth of the thrust wedge is inferred to have commenced in the Pliocene prior to 3 ± 1 Ma, but much of the wedge developed in the Quaternary. The spatial distribution of thrusting has varied through time, with most late Quaternary shortening occurring on structures within 10 km of the right-stepping deformation front. Estimates of the magnitude and rates of deformation indicate that the wedge accommodates a significant component of the oblique convergence between the Pacific and Australian Plates. Shortening of up to 7.3 ± 1.4 km and 9.1 ± 1.8 km within the southern and central parts of the wedge, respectively, represent about 5,15% of the total 70,140 km of shortening predicted across the plate boundary since 6.4 Ma, and about 10,30% since 3 Ma. Late Quaternary shortening rates of the order of 1,5 mm yr,1, estimated across both the northern and southern parts of the wedge, represent about 10,50 and 5,21% of the total NUVEL-1 A shortening across the plate boundary at these respective latitudes, implying that most shortening is occurring onshore. Furthermore, possible oblique-slip thrusting within the wedge may be accommodating boundary-parallel displacement of 0,6 mm yr,1, representing 0,17% of the total predicted within the plate boundary. [source]

Streamer tomography velocity models for the Gulf of Corinth and Gulf of Itea, Greece

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2004
Barry C. Zelt
SUMMARY The Gulf of Corinth (GOC), Greece is a rapidly extending, active continental rift with a record of large, damaging earthquakes. An extensive multichannel seismic (MCS) survey of the GOC conducted in 2001 provided, in addition to the processed MCS images, the opportunity to constrain velocity structure using refracted arrivals recorded along a 6-km-long streamer. We use first-arrival traveltimes to derive tomographic P -wave velocity models for several profiles collected in the central portion of the GOC. Eight of the profiles are closely spaced, north,south lines crossing the GOC and extending into the Gulf of Itea (GOI); a ninth profile is an east,west-oriented tie line. The N,S profiles image the relatively simple velocity structure of the deep Corinth rift basin and more complicated structure of the northern margin of the currently active rift. Integration of the velocity models with migrated MCS sections shows that south of the GOI the basement, which comprises Mesozoic nappes, occurs at a velocity of 4.5 km s,1 in the velocity models, although the actual velocity at, or just below, the top of basement is probably closer to 5,5.5 km s,1. The maximum sediment thickness in the Corinth basin is 2.2 km. The basement shallows to the north into a fault-bounded terrace in the central region between the two gulfs. Sediment cover in this central region decreases in thickness from west to east. Beneath the GOI, low average velocities beneath the rift-onset reflector indicate the presence of pre-rift sediments. The pre-rift velocity structure in the GOI is complex, with significant lateral variation from west to east. The E,W line shows that high-velocity basement is shallow (,1 km depth) and flat to the west of the GOI but dips ,20° east down to ,1.5 km beneath the pre-rift sediments of the GOI. [source]

Lithosphere structure of Europe and Northern Atlantic from regional three-dimensional gravity modelling

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2002
T. P. Yegorova
Summary Large-scale 3D gravity modelling using data averaged on a 1° grid has been performed for the whole European continent and part of the Northern Atlantic. The model consists of two regional layers of variable thickness,the sediments and the crystalline crust, bounded by reliable seismic horizons,the ,seismic' basement and the Moho surface. Inner heterogeneity of the model layers was taken into account in the form of lateral variation of average density depending on the type of geotectonic unit. Density parametrization of the layers was made using correlation functions between velocity and density. For sediments, sediment consolidation with depth was taken into account. Offshore a sea water layer was included in the model. As a result of the modelling, gravity effects of the whole model and its layers were calculated. Along with the gravity modelling an estimation of isostatic equilibrium state has been carried out for the whole model as well as for its separate units. Residual gravity anomalies, obtained by subtracting the gravity effect of the crust from the observed field, reach some hundred mGal (10,5 m s,2) in amplitude; they are mainly caused by density heterogeneities in the upper mantle. A mantle origin of the residual anomalies is substantiated by their correlation with the upper-mantle heterogeneities revealed by both seismological and geothermal studies. Regarding the character of the mantle gravity anomalies, type of isostatic compensation, crustal structure, age and supposed type of endogenic regime, a classification of main geotectonic units of the continent was made. As a result of the modelling a clear division of the continent into two large blocks,Precambrian East-European platform (EEP) and Variscan Western Europe,has been confirmed by their specific mantle gravity anomalies (0 ÷ 50 × 10,5 m s,2 and ,100 ÷,150 × 10,5 m s,2 correspondingly). This division coincides with the Tornquist,Teisseyre Zone (TTZ), marked by a gradient zone of mantle anomalies. In the central part of the EEP (over the Russian plate) an extensive positive mantle anomaly, probably indicating a core of ancient consolidation of the EEP, has been distinguished. To the west and to the east of this anomaly positive mantle anomalies occur, which coincide with a deep suture zone (TTZ) and an orogenic belt (the Urals). Positive mantle anomalies of the Alps, the Adriatic plate and the Calabrian Arc, correlating well with both high-velocity domains in the upper mantle and reduced temperatures at the subcrustal layer, are caused by thickened lithosphere below these structures. Negative mantle anomalies, revealed in the Western Mediterranean Basin and in the Pannonian Basin, are the result of thermal expansion of the asthenosphere shallowing to near-Moho depths below these basins. [source]

Upward continuation with topographic datuming operator: the integrated wave equation datuming scheme revised

GEOPHYSICAL PROSPECTING, Issue 6 2009
Kai Yang
ABSTRACT The integrated wave equation datuming scheme is an integrated datuming method to implement the wave equation velocity replacement with a one-way extrapolator for overthrust data. The integrated wave equation datuming scheme is accomplished by downward continuing the data from the topography to a non-planar base of the weathering layer, followed by an upward continuation from the base of the weathering layer to a final planar datum. Both the downward and the upward continuation are accomplished by a one-way extrapolator by a classical ,two-step' implementation. The topographic datuming operator is a distinctive technique by which the shots and receivers are downward continued simultaneously. Although its inadaptability to non-planar datum and the strong lateral variation of the near-surface structure constrains its application to the downward continuation of the overthrust data, we found that the ,one-step' feature of the topographic datuming operator is extremely suitable for the upward continuation. In this paper, the topographic datuming operator is incorporated into the integrated wave equation datuming scheme to perform upward continuation instead of a one-way extrapolator , thus a revised integrated wave equation datuming scheme is constructed. Compared with the previous scheme, the computatiol cost of the revised integrated wave equation datuming scheme is greatly reduced, thus the applicability of the integrated wave equation datuming scheme is further improved. The synthetic and real data examples demonstrate its effectiveness and efficiency. [source]

P- and S-wave velocities of consolidated sediments from a seafloor seismic survey in the North Celtic Sea Basin, offshore Ireland

GEOPHYSICAL PROSPECTING, Issue 2 2008
Donna J. Shillington
ABSTRACT A geophysical survey was conducted over a hydrocarbon prospect in the North Celtic Sea Basin using a small array of ocean-bottom seismographs (OBSs). The purpose of this study was to determine the ratio of compressional (P)- to shear (S)-wave velocity of consolidated sedimentary rocks in order to constrain possible subsurface variations in pore-fluid content. The ratio of VP and VS is known to be particularly sensitive to lithology, porosity and pore-fluid content, making it a useful parameter for evaluating hydrocarbon prospects. OBSs offer a relatively cheap and time-effective means of acquiring multi-component data compared with ocean-bottom cables. In this contribution, we demonstrate the ability of an OBS survey comprising three pairs of two OBSs spaced at 1.6 km to recover lateral variations in the VP/VS ratio. A key requirement of this type of study is that S waves will be generated by mode conversions in the subsurface, since they cannot be generated in nor travel through fluids. In this survey, the contrast in physical properties of the hard seabed of the North Celtic Sea Basin provided a means of generating converted S waves. Two-dimensional ray-tracing and forward modelling was used to create both VP and VS models along a profile crossing the Blackrock prospect in the North Celtic Sea Basin. These models comprise four layers and extend to a maximum depth of 1.1 km. The observed northward decrease in the VP/VS ratio at depths of 500,1000 m below the seafloor in the study area is interpreted to represent lateral variation in the amount of gas present in the pore space of Upper Cretaceous chalks and shales overlying the prospective reservoir. [source]

Steering of experimental channels by lateral basin tilting

BASIN RESEARCH, Issue 3 2010
Wonsuck Kim
ABSTRACT A major issue in tectonics and sedimentation is the role of cross-stream tectonic tilting in steering channels. The general idea is that channels will be attracted to lateral maxima in subsidence rate. A physical experiment performed in 1999 at the St. Anthony Falls Laboratory, however, was in conflict with the idea and showed that fluvial channels and resulting stratigraphy can be insensitive to even relatively strong lateral variation in subsidence. Here, we present results from an experiment which uses a simplified relay-ramp geometry with laterally variable uplift and subsidence to test a hypothesis developed from the earlier experiment: Tectonic tilting steers channels only when the ratio of the time scales describing lateral channel mobility to tectonic deformation is sufficiently large. Occupation time by experimental channels and sand fraction in the deposit (a proxy for channel deposition) both increase with subsidence rate indicating strong steering of channels by tectonic forcing. We also found that, due to local incision, uplift lengthened the time scale for lateral channel migration relative to subsidence. Comparing channel mobility at the beginning of the experiment, with no tectonic forcing, to later tectonic stages of the experiment indicates that active tectonics increased the channel time scale. The interplay of channel steering with uplift and subsidence led to cyclic appearance and disappearance of an autogenic lake in the hanging-wall basin. This lake was associated with alternation between channels going around vs. across the adjoining upstream uplifted footwall region. This creation and filling of the lake under constant tectonic forcing (constant fault slip rate) in the hanging wall created subaerial fan-delta parasequences separated by fluvial deposits. [source]

Effects of Near-Surface Absorption on Reflection Characteristics of Continental Interbedded Strata: the Dagang Oilfield as an Example

ACTA GEOLOGICA SINICA (ENGLISH EDITION), Issue 5 2010
LI Guofa
Abstract: Due to the effects of seismic wave field interference, the reflection events generated from interbedded and superposed sand and shale strata no longer have an explicit corresponding relationship with the geological interface. The absorption of the near-surface layer decreases the resolution of the seismic wavelet, intensifies the interference of seismic reflections from different sand bodies, and makes seismic data interpretation of thin interbedded strata more complex and difficult. In order to concretely investigate and analyze the effects of the near-surface absorption on seismic reflection characteristics of interbedded strata, and to make clear the ability of current technologies to compensate the near-surface absorption, a geological model of continental interbedded strata with near-surface absorption was designed, and the prestack seismic wave field was numerically simulated with wave equations. Then, the simulated wave field was processed by the prestack time migration, the effects of near-surface absorption on prestack and poststack reflection characteristics were analyzed, and the near-surface absorption was compensated for by inverse Q -filtering. The model test shows that: (1) the reliability of prediction and delineation of a continental reservoir with AVO inversion is degraded due to the lateral variation of the near-surface structure; (2) the corresponding relationships between seismic reflection events and geological interfaces are further weakened as a result of near-surface absorption; and (3) the current technology of absorption compensation probably results in false geological structure and anomaly. Based on the model experiment, the real seismic data of the Dagang Oil Field were analyzed and processed. The seismic reflection characteristics of continental interbedded strata were improved, and the reliability of geological interpretation from seismic data was enhanced. [source]

Crustal Composition of China Continent Constrained from Heat Flow Data and Helium Isotope Ratio of Underground Fluid

ACTA GEOLOGICA SINICA (ENGLISH EDITION), Issue 1 2010
Yang WANG
Abstract: Based on conservation of energy principle and heat flow data in China continent, the upper limit of 1.3 ,W/m3 heat production is obtained for continental crust in China. Furthermore, using the data of heat flow and helium isotope ratio of underground fluid, the heat productions of different tectonic units in China continent are estimated in range of 0.58,1.12 ,W/m3 with a median of 0.85 ,W/m3. Accordingly, the contents of U, Th and K2O in China crust are in ranges of 0.83,1.76 ,g/g, 3.16,6.69 ,g/g, and 1.0%,2.12%, respectively. These results indicate that the abundance of radioactive elements in the crust of China continent is much higher than that of Archean crust; and this fact implies China's continental crust is much evolved in chemical composition. Meanwhile, significant lateral variation of crustal composition is also exhibited among different tectonic units in China continent. The crust of eastern China is much enriched in incompatible elements such as U, Th and K than that of western China; and the crust of orogenic belts is more enriched than that of platform regions. It can also be inferred that the crusts of eastern China and orogenic belts are much felsic than those of western China and platform regions, respectively, derived from the positive correlation between the heat production and SiO2 content of bulk crust. This deduction is consistent with the results derived from the crustal seismic velocity data in China. According to the facts of the lower seismic velocity of China than the average value of global crust, and the higher heat production of China continent compared with global crust composition models published by previous studies, it is deduced that the average composition models of global continent crust by Rudnick and Fountain (1995), Rudnick and Gao (2003), Weaver and Tarney (1984), Shaw et al. (1986), and Wedepohl (1995) overestimate the abundance of incompatible elements such as U, Th and K of continental crust. [source]

A glassy lowermost outer core

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2009
Vernon F. Cormier
SUMMARY New theories for the viscosity of metallic melts at core pressures and temperatures, together with observations of translational modes of oscillation of Earth's solid inner core, suggest a rapid increase in the dynamic viscosity near the bottom of the liquid outer core. If the viscosity of the lowermost outer core (F region) is sufficiently high, it may be in a glassy state, characterized by a frequency dependent shear modulus and increased viscoselastic attenuation. In testing this hypothesis, the amplitudes of high-frequency PKiKP waves are found to be consistent with an upper bound to shear velocity in the lowermost outer core of 0.5 km s,1 at 1 Hz. The fit of a Maxwell rheology to the frequency dependent shear modulus constrained by seismic observations at both low and high-frequency favours a model of the F region as a 400-km-thick chemical boundary layer. This layer has both a higher density and higher viscosity than the bulk of the outer core, with a peak viscosity on the order of 109 Pa s or higher near the inner core boundary. If lateral variations in the F region are confirmed to correlate with lateral variations observed in the structure of the uppermost inner core, they may be used to map differences in the solidification process of the inner core and flow in the lowermost outer core. [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]

The lowermost mantle beneath northern Asia,II.

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2002
Evidence for lower-mantle anisotropy
Summary We have analysed prediffracted S -waves with turning points beneath northern Siberia in a study of anisotropy in the lowermost mantle. Deep-focus earthquakes beneath the Marianas, Izu Bonin and the Sea of Japan recorded at stations in western Europe are used. A correction for upper-mantle anisotropy is applied to the data. Comparisons of the data with synthetic data for models with and without a high velocity D, layer suggest that there is a velocity discontinuity at the top of the D, region and that the style of anisotropy is transversely isotropic in this region. Time separations between S -waves on the radial and transverse component show a weak trend where the separation increases with epicentral distance. A normalization of this separation with the travel distance within D, (300 km thick in this region) suggests that the anisotropy is uniformly distributed within this layer and has an average value of 0.5 per cent. A combination of different studies which investigate the structure of the lowermost mantle beneath Europe and northern Siberia reveals a complicated picture. Tomographic models from this area and evidence of D, anisotropy, lower mantle scatterers, reflections from a D, discontinuity and ultra-low-velocity zones suggest two distinct regions. One exhibits high velocities, D, anisotropy, a D, discontinuity and no evidence of scatterers or ultra-low-velocity zones. These features are likely associated with the palaeosubduction of the Izanagi plate well into the lowermost mantle. The other region has a lower overall velocity and shows evidence of scatterers and ultra-low-velocity zones, perhaps suggesting the presence of partial melt. These results suggest dramatic lateral variations in the nature of the lowermost mantle beneath northern Asia over a length scale of roughly 30 degrees. [source]

Spectral-element simulations of global seismic wave propagation,II.

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2002
Three-dimensional models, oceans, rotation, self-gravitation
Summary We simulate global seismic wave propagation based upon a spectral-element method. We include the full complexity of 3-D Earth models, i.e. lateral variations in compressional-wave velocity, shear-wave velocity and density, a 3-D crustal model, ellipticity, as well as topography and bathymetry. We also include the effects of the oceans, rotation and self-gravitation in the context of the Cowling approximation. For the oceans we introduce a formulation based upon an equivalent load in which the oceans do not need to be meshed explicitly. Some of these effects, which are often considered negligible in global seismology, can in fact play a significant role for certain source,receiver configurations. Anisotropy and attenuation, which were introduced and validated in a previous paper, are also incorporated in this study. The complex phenomena that are taken into account are introduced in such a way that we preserve the main advantages of the spectral-element method, which are an exactly diagonal mass matrix and very high computational efficiency on parallel computers. For self-gravitation and the oceans we benchmark spectral-element synthetic seismograms against normal-mode synthetics for the spherically symmetric reference model PREM. The two methods are in excellent agreement for all body- and surface-wave arrivals with periods greater than about 20 s in the case of self-gravitation and 25 s in the case of the oceans. At long periods the effect of gravity on multiorbit surface waves up to R4 is correctly reproduced. We subsequently present results of simulations for two real earthquakes in fully 3-D Earth models for which the fit to the data is significantly improved compared with classical normal-mode calculations based upon PREM. For example, we show that for trans-Pacific paths the Rayleigh wave can arrive more than a minute earlier than in PREM, and that the Love wave is much shorter in duration. [source]

Geodetic imaging: reservoir monitoring using satellite interferometry

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 3 2002
D. W. Vasco
Summary Fluid fluxes within subsurface reservoirs give rise to surface displacements, particularly over periods of a year or more. Observations of such deformation provide a powerful tool for mapping fluid migration within the Earth, providing new insights into reservoir dynamics. In this paper we use Interferometric Synthetic Aperture Radar (InSAR) range changes to infer subsurface fluid volume strain at the Coso geothermal field. Furthermore, we conduct a complete model assessment, using an iterative approach to compute model parameter resolution and covariance matrices. The method is a generalization of a Lanczos-based technique which allows us to include fairly general regularization, such as roughness penalties. We find that we can resolve quite detailed lateral variations in volume strain both within the reservoir depth range (0.4,2.5 km) and below the geothermal production zone (2.5,5.0 km). The fractional volume change in all three layers of the model exceeds the estimated model parameter uncertainty by a factor of two or more. In the reservoir depth interval (0.4,2.5 km), the predominant volume change is associated with northerly and westerly oriented faults and their intersections. However, below the geothermal production zone proper [the depth range 2.5,5.0 km], there is the suggestion that both north- and northeast-trending faults may act as conduits for fluid flow. [source]

Regional teleseismic tomography of the western Lachlan Orogen and the Newer Volcanic Province, southeast Australia

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 2 2002
Frank M. Graeber
Summary From 1998 May to September a portable array of 40 short-period digital seismograph stations was operated in western Victoria, southeast Australia, across the western end of the mid-Paleozoic Lachlan Foldbelt and the Newer Volcanic Province. Consisting of four parallel, almost W,E-oriented receiver lines, the array covered an area of about 270 × 150 km2. The major aim of the LF98 (Lachlan Foldbelt survey 1998) project is to map lateral variations in P -wave speeds (Vp) in the crust and upper mantle using teleseismic arrival time tomography, primarily in order to investigate whether the major surface structural zones are associated with seismic velocity signatures at depth. Little a priori information from seismic profiling is available. We invert 4067 relative arrival time residuals for a minimum structure Vp model in the upper few hundred km using non-linear iteration and 3-D ray tracing. The most prominent negative anomaly (,3.8 per cent) in Vp is found at a depth of about 45 km underneath the eastern part of the Newer Volcanic Province. It correlates spatially with the highest density of Pliocene and Pleistocene eruption centres northwest of Melbourne, and is therefore interpreted as a hotspot-related high-temperature anomaly causing reduced mantle velocities. The related coherent volume of significantly lower than average velocities extends down to depths greater than 100 km in the east, and extends west underneath the Newer Volcanic Province. A strong velocity contrast, with average velocities ,2 per cent greater in the west, is found down to about 100 km across the Moyston Fault Zone, which forms the major structural boundary between the early-Paleozoic Delamerian Orogen in the west and the Lachlan Orogen in the east. This result suggests that the Moyston Fault Zone should be seen as a major lithospheric boundary. In the south this boundary is also expressed by a distinct discontinuity in Sr-isotopic ratios of xenoliths (the so-called Mortlake discontinuity) and a change in the geochemistry of plutons of similar age. However, if the east to west velocity contrast originally existed in this southern zone, it is now overprinted by the thermally reduced mantle velocities beneath the Newer Volcanic Province. [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]

Non-uniqueness with refraction inversion , the Mt Bulga shear zone

GEOPHYSICAL PROSPECTING, Issue 4 2010
Derecke Palmer
ABSTRACT The tau-p inversion algorithm is widely employed to generate starting models with many computer programs that implement refraction tomography. However, this algorithm can frequently fail to detect even major lateral variations in seismic velocities, such as a 50 m wide shear zone, which is the subject of this study. By contrast, the shear zone is successfully defined with the inversion algorithms of the generalized reciprocal method. The shear zone is confirmed with a 2D analysis of the head wave amplitudes, a spectral analysis of the refraction convolution section and with numerous closely spaced orthogonal seismic profiles recorded for a later 3D refraction investigation. Further improvements in resolution, which facilitate the recognition of additional zones with moderate reductions in seismic velocity, are achieved with a novel application of the Hilbert transform to the refractor velocity analysis algorithm. However, the improved resolution also requires the use of a lower average vertical seismic velocity, which accommodates a velocity reversal in the weathering. The lower seismic velocity is derived with the generalized reciprocal method, whereas most refraction tomography programs assume vertical velocity gradients as the default. Although all of the tomograms are consistent with the traveltime data, the resolution of each tomogram is comparable only with that of the starting model. Therefore, it is essential to employ inversion algorithms that can generate detailed starting models, where detailed lateral resolution is the objective. Non-uniqueness can often be readily resolved with head wave amplitudes, attribute processing of the refraction convolution section and additional seismic traverses, prior to the acquisition of any borehole data. It is concluded that, unless specific measures are taken to address non-uniqueness, the production of a single refraction tomogram that fits the traveltime data to sufficient accuracy does not necessarily demonstrate that the result is either correct, or even the most probable. [source]

Stacking velocities in the presence of overburden velocity anomalies

GEOPHYSICAL PROSPECTING, Issue 3 2009
Emil Blias
ABSTRACT Lateral velocity changes (velocity anomalies) in the overburden may cause significant oscillations in normal moveout velocities. Explicit analytical moveout formulas are presented and provide a direct explanation of these lateral fluctuations and other phenomena for a subsurface with gentle deep structures and shallow overburden anomalies. The analytical conditions for this have been derived for a depth-velocity model with gentle structures with dips not exceeding 12°. The influence of lateral interval velocity changes and curvilinear overburden velocity boundaries can be estimated and analysed using these formulas. An analytical approach to normal moveout velocity analysis in a laterally inhomogeneous medium provides an understanding of the connection between lateral interval velocity changes and normal moveout velocities. In the presence of uncorrected shallow velocity anomalies, the difference between root-mean-square and stacking velocity can be arbitrarily large to the extent of reversing the normal moveout function around normal incidence traveltimes. The main reason for anomalous stacking velocity behaviour is non-linear lateral variations in the shallow overburden interval velocities or the velocity boundaries. A special technique has been developed to determine and remove shallow velocity anomaly effects. This technique includes automatic continuous velocity picking, an inversion method for the determination of shallow velocity anomalies, improving the depth-velocity model by an optimization approach to traveltime inversion (layered reflection tomography) and shallow velocity anomaly replacement. Model and field data examples are used to illustrate this technique. [source]

P- and S-wave velocities of consolidated sediments from a seafloor seismic survey in the North Celtic Sea Basin, offshore Ireland

Seismic attenuation in Faroe Islands basalts

GEOPHYSICAL PROSPECTING, Issue 1 2008
F. Shaw
ABSTRACT We analysed vertical seismic profiling (VSP) data from two boreholes at Glyvursnes and Vestmanna on the island of Streymoy, Faroe Islands, to determine the magnitude and causes of seismic attenuation in sequences of basalt flows. The work is part of SeiFaBa, a major project integrating data from vertical and offset VSP, surface seismic surveys, core samples and wireline log data from the two boreholes. Values of effective seismic quality factor (Q) obtained at Glyvursnes and Vestmanna are sufficiently low to significantly degrade the quality of a surface reflection seismic image. This observation is consistent with results from other VSP experiments in the North Atlantic region. We demonstrate that the most likely cause of the low values of effective Q at Glyvursnes and Vestmanna is a combination of 1D scattering and intrinsic attenuation due to seismic wave-induced fluid flow within pores and micro-cracks. Tests involving 3D elastic wave numerical modelling with a hypothetical basalt model based on field observations, indicate that little scattering attenuation is caused by lateral variations in basalt structure. [source]

Detection and analysis of LNAPL using the instantaneous amplitude and frequency of ground-penetrating radar data

GEOPHYSICAL PROSPECTING, Issue 1 2002
Luciana Orlando
This paper reports the results of using the ground-penetrating radar (GPR) method to detect light non-aqueous phase liquids (LNAPL) floating on the water table in an area where the thickness of LNAPL present ranges from a few centimetres to several decimetres. To understand the GPR response in this context, GPR theoretical models are calculated using information from the literature and hydrogeological field data. The study revealed that in the case of LNAPL floating on the water table in a static condition, there is an increase in the reflection amplitude from the water table due to the decrease in the capillary fringe. Nevertheless the amplitude of reflection from the water table can discriminate the contaminated from the non-contaminated zone. Apart from an analysis of the real traces, the analysis of some attributes of the complex trace, instantaneous amplitude, phase and frequency, are also good tools to detect hydrocarbons floating on the water table. Such attributes, depending on both the signal frequency and the hydrocarbon thickness, can also give information about the thickness of the hydrocarbon layer. It is concluded that analysing the lateral variations in signal amplitude of the real trace and in the amplitude, phase and instantaneous frequency of the complex signal permits the delimiting of the area polluted by the hydrocarbon. [source]

Magnetic signal prospecting using multiparameter measurements: the case study of the Gallic Site of Levroux

ARCHAEOLOGICAL PROSPECTION, Issue 3 2010
M. Pétronille
Abstract The ,magnetic signal' that combines both the induced (Ji) and the remanent (Jr) magnetization is widely used in archaeological and pedological prospecting. Magnetic prospecting recording the lateral variations of the total magnetization is the most frequently used measurement before in-phase magnetic susceptibility (Kph) and magnetic viscosity (Kqu) mapping. The work presented here brings together three types of prospecting technique: magnetic field survey and electromagnetic measurements with both frequency and time domain devices that measure magnetic susceptibility and viscosity respectively. The site studied, the Gallic town of Levroux (Indre, France), is particularly interesting because it includes features such as pits and ditches dug into the calcareous substratum partly filled with topsoil and with residues of different metallurgical and fire activities. The field results indicated anomalies with different types of characterization: (i) many compact features filled with magnetic, electrically conductive and minimally viscous materials; and (ii) elongated anomalies characterized by lower magnetic properties and electric conductivity but relatively higher magnetic viscosity than those of the compact features. In addition to the location of the features, the combination of the information brought by the different types of measurements allows us to evaluate the possible erosion of their upper parts by ploughing, to assess their depth (never deeper than 1.30,m) and to precise the nature of the feature's fill. Copyright © 2010 John Wiley & Sons, Ltd. [source]