Merged bibliography of the quadrennial Seismology Report (1983-1986) to the International Union of Geodesy and Geophysics
Explore the source record for details and available documents.
Open-File Report↗
SEARCH · Search USGS
Results for “International Geophysics”
Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
A trans-dimensional Bayesian Markov chain Monte Carlo algorithm for model assessment using frequency-domain electromagnetic data
A meaningful interpretation of geophysical measurements requires an assessment of the space of models that are consistent with the data, rather than just a single, ‘best’ model which does not convey information about parameter uncertainty. For this purpose, a trans-dimensional Bayesian Markov chain Monte Carlo (MCMC) algorithm is developed for assessing frequencydomain electromagnetic (FDEM) data acquired from airborne or ground-based systems. By sampling the distribution of models that are consistent with measured data and any prior knowledge, valuable inferences can be made about parameter values such as the likely depth to an interface, the distribution of possible resistivity values as a function of depth and non-unique relationships between parameters. The trans-dimensional aspect of the algorithm allows the number of layers to be a free parameter that is controlled by the data, where models with fewer layers are inherently favoured, which provides a natural measure of parsimony and a significant degree of flexibility in parametrization. The MCMC algorithm is used with synthetic examples to illustrate how the distribution of acceptable models is affected by the choice of prior information, the system geometry and configuration and the uncertainty in the measured system elevation. An airborne FDEM data set that was acquired for the purpose of hydrogeological characterization is also studied. The results compare favorably with traditional least-squares analysis, borehole resistivity and lithology logs from the site, and also provide new information about parameter uncertainty necessary for model assessment.
Geophysical Journal International↗
High-frequency Born synthetic seismograms based on coupled normal modes
High-frequency and full waveform synthetic seismograms on a 3-D laterally heterogeneous earth model are simulated using the theory of coupled normal modes. The set of coupled integral equations that describe the 3-D response are simplified into a set of uncoupled integral equations by using the Born approximation to calculate scattered wavefields and the pure-path approximation to modulate the phase of incident and scattered wavefields. This depends upon a decomposition of the aspherical structure into smooth and rough components. The uncoupled integral equations are discretized and solved in the frequency domain, and time domain results are obtained by inverse Fourier transform. Examples show the utility of the normal mode approach to synthesize the seismic wavefields resulting from interaction with a combination of rough and smooth structural heterogeneities. This approach is applied to an ∼4 Hz shallow crustal wave propagation around the site of the San Andreas Fault Observatory at Depth (SAFOD).
Geophysical Journal International↗
The value of geophysical methods in ground‐water studies
Two meanings are unfortunately given to the term geophysics. In the broad sense, based on the etymology of the term, geophysics means the physics of the Earth. This is its significance in the names “Geophysical Laboratory of the Carnegie Institution of Washington,” “International Union of Geodesy and Geophysics,” and “American Geophysical Union.” In this sense the sciences of geology and hydrology are largely geophysics, although in part they are geochemistry or biology. Geologists and hydraulic engineers who are elected to membership in the American Geophysical Union are, however, often surprised to find that they are regarded as “geophysicists.” I hope that when the Union establisnes a journal, its name will include the expression “Earth physics” or “Physics of the Earth,” in order that there will be no misunderstanding as to its scope.
Eos, Transactions, American Geophysical Union↗
Locations and magnitudes of historical earthquakes in the Sierra of Ecuador (1587–1996)
The whole territory of Ecuador is exposed to seismic hazard. Great earthquakes can occur in the subduction zone (e.g. Esmeraldas, 1906, M w 8.8), whereas lower magnitude but shallower and potentially more destructive earthquakes can occur in the highlands. This study focuses on the historical crustal earthquakes of the Andean Cordillera. Several large cities are located in the Interandean Valley, among them Quito, the capital (∼2.5 millions inhabitants). A total population of ∼6 millions inhabitants currently live in the highlands, raising the seismic risk. At present, precise instrumental data for the Ecuadorian territory is not available for periods earlier than 1990 (beginning date of the revised instrumental Ecuadorian seismic catalogue); therefore historical data are of utmost importance for assessing seismic hazard. In this study, the Bakun & Wentworth method is applied in order to determine magnitudes, locations, and associated uncertainties for historical earthquakes of the Sierra over the period 1587–1976. An intensity-magnitude equation is derived from the four most reliable instrumental earthquakes ( M w between 5.3 and 7.1). Intensity data available per historical earthquake vary between 10 (Quito, 1587, Intensity ≥VI) and 117 (Riobamba, 1797, Intensity ≥III). The bootstrap resampling technique is coupled to the B&W method for deriving geographical confidence contours for the intensity centre depending on the data set of each earthquake, as well as confidence intervals for the magnitude. The extension of the area delineating the intensity centre location at the 67 per cent confidence level (±1σ) depends on the amount of intensity data, on their internal coherence, on the number of intensity degrees available, and on their spatial distribution. Special attention is dedicated to the few earthquakes described by intensities reaching IX, X and XI degrees. Twenty-five events are studied, and nineteen new epicentral locations are obtained, yielding equivalent moment magnitudes between 5.0 and 7.6. Large earthquakes seem to be related to strike slip faults between the North Andean Block and stable South America to the east, while moderate earthquakes ( M w ≤ 6) seem to be associated with to thrust faults located on the western internal slopes of the Interandean Valley.
Geophysical Journal International↗
Contemporary deformation in the Yakima fold and thrust belt estimated with GPS
Geodetic, geologic and palaeomagnetic data reveal that Oregon (western USA) rotates clockwise at 0.3 to 1.0° Ma −1 (relative to North America) about an axis near the Idaho–Oregon–Washington border, while northeast Washington is relatively fixed. This rotation has been going on for at least 15 Ma. The Yakima fold and thrust belt (YFTB) forms the boundary between northern Oregon and central Washington where convergence of the clockwise-rotating Oregon block is apparently accommodated. North–south shortening across the YFTB has been thought to occur in a fan-like manner, increasing in rate to the west. We obtained high-accuracy, high-density geodetic GPS measurements in 2012–2014 that are used with earlier GPS measurements from the 1990s to characterize YFTB kinematics. The new results show that the deformation associated with the YFTB starts at the Blue Mountains Anticline in northern Oregon and extends north beyond the Frenchman Hills in Washington, past the epicentre of the 1872 M w 7.0 Entiat earthquake to 49°N. The north–south strain rate across the region is 2 to 3 × 10 −9 yr −1 between the volcanic arc and the eastern edge of the YFTB (241.0°E); east of there it drops to about 10 −9 yr −1 . At the eastern boundary of the YFTB, faults and earthquake activity are truncated by a north-trending, narrow zone of deformation that runs along the Pasco Basin and Moses Lake regions near 240.9°E. This zone, abutting the Department of Energy Hanford Nuclear Reservation, accommodates about 0.5 mm yr −1 of east to northeast shortening. A similar zone of N-trending transpression is seen along 239.9°E where there is a change in the strike of the Yakima folds. The modern deformation of the YFTB is about 600 km wide from south to north and internally may be controlled by pre-existing crustal structure.
Geophysical Journal International↗
National- to continental-scale governmental geophysical efforts for critical mineral mapping, USA
The U.S. Geological Survey (USGS) has established robust collaborations with domestic state and international geological surveys to provide geophysical and other types of earth science data that act to underpin critical mineral research efforts across the United States, Canada, and Australia. The Earth Mapping Resource Initiative (EMRI) is a national-scale collaborative effort with state geological surveys to improve geophysical and geological data to advance our understanding of the United States’ critical mineral endowment. The Critical Mineral Mapping Initiative (CMMI) is a tri-national collaboration with the federal geological surveys of Canada and Australia to conduct research that will aid in identifying new areas with potential for critical mineral deposits across all three countries. This study describes the important interplay between the EMRI and CMMI and how each act in a complementary fashion to advance critical mineral research. We present examples that illustrate how magnetic anomaly data are used to define critical mineral prospectivity for Mississippi Valley-type (MVT) Zn-Pb mineral systems and illustrate how CMMI magnetic derivative maps were considered into USGS’ EMRI efforts to acquire modern high-resolution airborne geophysical data over a large area within the US Midcontinent.
Conference Paper↗
Body-wave traveltime and amplitude shifts from asymptotic travelling wave coupling
We explore the sensitivity of finite-frequency body-wave traveltimes and amplitudes to perturbations in 3-D seismic velocity structure relative to a spherically symmetric model. Using the approach of coupled travelling wave theory, we consider the effect of a structural perturbation on an isolated portion of the seismogram. By convolving the spectrum of the differential seismogram with the spectrum of a narrow window taper, and using a Taylor's series expansion for wavenumber as a function of frequency on a mode dispersion branch, we derive semi-analytic expressions for the sensitivity kernels. Far-field effects of wave interactions with the free surface or internal discontinuities are implicitly included, as are wave conversions upon scattering. The kernels may be computed rapidly for the purpose of structural inversions. We give examples of traveltime sensitivity kernels for regional wave propagation at 1 Hz. For the direct SV wave in a simple crustal velocity model, they are generally complicated because of interfering waves generated by interactions with the free surface and the Mohorovic??ic?? discontinuity. A large part of the interference effects may be eliminated by restricting the travelling wave basis set to those waves within a certain range of horizontal phase velocity. ?? Journal compilation ?? 2006 RAS.
Geophysical Journal International↗
Seismotectonic framework of the 2010 February 27 M w 8.8 Maule, Chile earthquake sequence
After the 2010 M w 8.8 Maule earthquake, an international collaboration involving teams and instruments from Chile, the US, the UK, France and Germany established the International Maule Aftershock Deployment temporary network over the source region of the event to facilitate detailed, open-access studies of the aftershock sequence. Using data from the first 9-months of this deployment, we have analyzed the detailed spatial distribution of over 2500 well-recorded aftershocks. All earthquakes have been relocated using a hypocentral decomposition algorithm to study the details of and uncertainties in both their relative and absolute locations. We have computed regional moment tensor solutions for the largest of these events to produce a catalogue of 465 mechanisms, and have used all of these data to study the spatial distribution of the aftershock sequence with respect to the Chilean megathrust. We refine models of co-seismic slip distribution of the Maule earthquake, and show how small changes in fault geometries assumed in teleseismic finite fault modelling significantly improve fits to regional GPS data, implying that the accuracy of rapid teleseismic fault models can be substantially improved by consideration of existing fault geometry model databases. We interpret all of these data in an integrated seismotectonic framework for the Maule earthquake rupture and its aftershock sequence, and discuss the relationships between co-seismic rupture and aftershock distributions. While the majority of aftershocks are interplate thrust events located away from regions of maximum co-seismic slip, interesting clusters of aftershocks are identified in the lower plate at both ends of the main shock rupture, implying internal deformation of the slab in response to large slip on the plate boundary interface. We also perform Coulomb stress transfer calculations to compare aftershock locations and mechanisms to static stress changes following the Maule rupture. Without the incorporation of uncertainties in earthquake locations, just 55 per cent of aftershock nodal planes align with faults promoted towards failure by co-seismic slip. When epicentral uncertainties are considered (on the order of just ±2–3 km), 90 per cent of aftershocks are consistent with occurring along faults demonstrating positive stress transfer. These results imply large sensitivities of Coulomb stress transfer calculations to uncertainties in both earthquake locations and models of slip distributions, particularly when applied to aftershocks close to a heterogeneous fault rupture; such uncertainties should therefore be considered in similar studies used to argue for or against models of static stress triggering.
Maule↗
International heat flow commission celebrates 40 years
The outflow of heat from the Earth's interior is, in terms of energy, the most impressive terrestrial phenomenon. Its present rate of about 1021 joules per year is order‐of‐magnitudes greater than the energy dissipation of earthquakes or heat loss from volcanic eruptions. The study of the Earth's internal heat plays an important role in understanding the Earth's origin, internal constitution, and plate tectonics. The thermal structure of the Earth, generally referred to as the geothermics, or also as the (terrestrial) heat flow, has been studied for a long time. However, modern geothermics, one of the fundamental geophysical disciplines, is relatively young. The International Heat Flow Commission (IHFC), under which the academic geothermal research on the international scale is organized, was created only in 1963. The IHFC, operating under the International Association of Seismology and Physics of the Earth's Interior (IASPEI), covers a broad scope of geophysical studies, and links the activities of other associations of the International Union of Geodesy and Geophysics (IUGG).
Eos, Earth and Space Science News↗
RMT focal plane sensitivity to seismic network geometry and faulting style
Modern tectonic studies often use regional moment tensors (RMTs) to interpret the seismotectonic framework of an earthquake or earthquake sequence; however, despite extensive use, little existing work addresses RMT parameter uncertainty. Here, we quantify how network geometry and faulting style affect RMT sensitivity. We examine how data-model fits change with fault plane geometry (strike and dip) for varying station configurations. We calculate the relative data fit for incrementally varying geometries about a best-fitting solution, applying our workflow to real and synthetic seismograms for both real and hypothetical station distributions and earthquakes. Initially, we conduct purely observational tests, computing RMTs from synthetic seismograms for hypothetical earthquakes and a series of well-behaved network geometries. We then incorporate real data and station distributions from the International Maule Aftershock Deployment (IMAD), which recorded aftershocks of the 2010 M W 8.8 Maule earthquake, and a set of regional stations capturing the ongoing earthquake sequence in Oklahoma and southern Kansas. We consider RMTs computed under three scenarios: (1) real seismic records selected for high data quality; (2) synthetic seismic records with noise computed for the observed source-station pairings and (3) synthetic seismic records with noise computed for all possible station-source pairings. To assess RMT sensitivity for each test, we observe the ‘fit falloff’, which portrays how relative fit changes when strike or dip varies incrementally; we then derive the ranges of acceptable strikes and dips by identifying the span of solutions with relative fits larger than 90 per cent of the best fit. For the azimuthally incomplete IMAD network, Scenario 3 best constrains fault geometry, with average ranges of 45° and 31° for strike and dip, respectively. In Oklahoma, Scenario 3 best constrains fault dip with an average range of 46°; however, strike is best constrained by Scenario 1, with a range of 26°. We draw two main conclusions from this study. (1) Station distribution impacts our ability to constrain RMTs using waveform time-series; however, in some tectonic settings, faulting style also plays a significant role and (2) increasing station density and data quantity (both the number of stations and the number of individual channels) does not necessarily improve RMT constraint. These results may be useful when organizing future seismic deployments (e.g. by concentrating stations in alignment with anticipated nodal planes), and in computing RMTs, either by guiding a more rigorous data selection process for input data or informing variable weighting among the selected data (e.g. by eliminating the transverse component when strike-slip mechanisms are expected).
Geophysical Journal International↗
Study of a prehistoric landslide using seismic reflection methods integrated with geological data in the Wasatch Mountains, Utah, USA
An integration of geological and geophysical techniques characterizes the internal and basal structure of a landslide along the western margin of the Wasatch Mountains in northern Utah, USA. The study area is within a region of planned and continuing residential development. The Little Valley Landslide is a prehistoric landslide as old as 13??ka B.P. Drilling and trenching at the site indicate that the landslide consists of chaotic and disturbed weathered volcanic material derived from Tertiary age volcanic rocks that comprise a great portion of the Wasatch Range. Five short high-resolution common mid-point seismic reflection profiles over selected portions of the site examine the feasibility of using seismic reflection to study prehistoric landslides in the Wasatch Mountain region. Due to the expected complexity of the near-surface geology, we have pursued an experimental approach in the data processing, examining the effects of muting first arrivals, frequency filtering, model-based static corrections, and seismic migration. The results provide a framework for understanding the overall configuration of the landslide, its basal (failure) surface, and the structure immediately underlying this surface. A glide surface or de??collement is interpreted to underlie the landslide suggesting a large mass movement. The interpretation of a glide surface is based on the onset of coherent reflectivity, calibrated by information from a borehole located along one of the seismic profiles. The glide surface is deepest in the center portion of the landslide and shallows up slope, suggesting a trough-like feature. This study shows that seismic reflection techniques can be successfully used in complex alpine landslide regions to (1) provide a framework in which to link geological data and (2) reduce the need for an extensive trenching and drilling program. ?? 2007 Elsevier B.V. All rights reserved.
Engineering Geology↗
Three-dimenstional crustal velocity structure beneath the Strait of Georgia, British Columbia
The Strait of Georgia is a topographic depression straddling the boundary between the Insular and Coast belts in southwestern British Columbia. Two shallow earthquakes located within the strait ( M = 4.6 in 1997 and M = 5.0 in 1975) and felt throughout the Vancouver area illustrate the seismic potential of this region. As part of the 1998 Seismic Hazards Investigation of Puget Sound (SHIPS) experiment, seismic instruments were placed in and around the Strait of Georgia to record shots from a marine source within the strait. We apply a tomographic inversion procedure to first-arrival traveltime data to derive a minimum-structure 3-D P -wave velocity model for the upper crust to about 13 km depth. We also present a 2-D velocity model for a profile orientated across the Strait of Georgia derived using a minimum-parameter traveltime inversion approach. This paper represents the first detailed look at crustal velocity variations within the major Cretaceous to Cenozoic Georgia Basin, which underlies the Strait of Georgia. The 3-D velocity model clearly delineates the structure of the Georgia Basin. Taking the 6 km s −1 isovelocity contour to represent the top of the underlying basement, the basin thickens from between 2 and 4 km in the northwestern half of the strait to between 8 and 9 km at the southeastern end of the study region. Basin velocities in the northeastern half are 4.5–6 km s −1 and primarily represent the Upper Cretaceous Nanaimo Group. Velocities to the south are lower (3–6 km s −1 ) because of the additional presence of the overlying Tertiary Huntingdon Formation and more recent sediments, including glacial and modern Fraser River deposits. In contrast to the relatively smoothly varying velocity structure of the basin, velocities of the basement rocks, which comprise primarily Palaeozoic to Jurassic rocks of the Wrangellia Terrane and possibly Jurassic to mid-Cretaceous granitic rocks of the Coast Belt, show significantly more structure, probably an indication of the varying basement rock lithologies. The 2-D velocity model more clearly reveals the velocity layering associated with the recent sediments, Huntingdon Formation and Nanaimo Group of the southern Georgia Basin, as well as the underlying basement. We interpret lateral variation in sub-basin velocities of the 2-D model as a transition from Wrangellian to Coast Belt basement rocks. The effect of the narrow, onshore–offshore recording geometry of the seismic experiment on model resolution was tested to allow a critical assessment of the validity of the 3-D velocity model. Lateral resolution throughout the model to a depth of 3–5 km below the top of the basement is generally 10–20 km.
British Columbia↗
Palaeomagnetic dating of widespread remagnetization on the southeastern border of the French Massif Central and implications for fluid and Mississippi Valley-type mineralization
Palaeomagnetic dating techniques have been applied to determine the age of fluid migration that produced the Mississippi Valley-type (MVT) Pb-Zn-Ba-F deposits in the Ce??vennes region of southern France. 15 sampling sites in two gently deformed areas around the Largentie??re and Croix-de-Pallie??res mines on the Ce??vennes border were selected for palaeomagnetic study. They yielded a very well-defined direction of remagnetization corresponding to an Early-Middle Eocene age. This remagnetization cannot be related to the formation of magnetic as a result of the transformation of smectite to illite because the latter has been well dated as a Mesozoic event. The magnetic overprint in this area is related to a chemical phenomenon during fluid migration. The age of remagnetization corresponds to a major uplift in the Pyre??ne??es mountains, located to the south of the Ce??vennes. This implies that fluid migration occurred from the south to the north as a result of hydraulic head established in the Pyre??ne??es orogenic belt during orogenesis and suggests that the MVT deposits in the Ce??vennes region formed from a gravity-driven fluid system as described by Garven & Freeze (1984a,b).
Geophysical Journal International↗
Seismic tomography shows that upwelling beneath Iceland is confined to the upper mantle
We report the results of the highest-resolution teleseismic tomography study yet performed of the upper mantle beneath Iceland. The experiment used data gathered by the Iceland Hotspot Project, which operated a 35-station network of continuously recording, digital, broad-band seismometers over all of Iceland 1996-1998. The structure of the upper mantle was determined using the ACH damped least-squares method and involved 42 stations, 3159 P-wave, and 1338 S-wave arrival times, including the phases P, pP, sP, PP, SP, PcP, PKIKP, pPKIKP, S, sS, SS, SKS and Sdiff. Artefacts, both perceptual and parametric, were minimized by well-tested smoothing techniques involving layer thinning and offset-and-averaging. Resolution is good beneath most of Iceland from ??? 60 km depth to a maximum of ??? 450 km depth and beneath the Tjornes Fracture Zone and near-shore parts of the Reykjanes ridge. The results reveal a coherent, negative wave-speed anomaly with a diameter of 200-250 km and anomalies in P-wave speed, Vp, as strong as -2.7 per cent and in S-wave speed, Vs, as strong as -4.9 per cent. The anomaly extends from the surface to the limit of good resolution at ??? 450 km depth. In the upper ??? 250 km it is centred beneath the eastern part of the Middle Volcanic Zone, coincident with the centre of the ??? 100 mGal Bouguer gravity low over Iceland, and a lower crustal low-velocity zone identified by receiver functions. This is probably the true centre of the Iceland hotspot. In the upper ??? 200 km, the low-wave-speed body extends along the Reykjanes ridge but is sharply truncated beneath the Tjornes Fracture Zone. This suggests that material may flow unimpeded along the Reykjanes ridge from beneath Iceland but is blocked beneath the Tjornes Fracture Zone. The magnitudes of the Vp, Vs and Vp/Vs anomalies cannot be explained by elevated temperature alone, but favour a model of maximum temperature anomalies <200 K, along with up to ??? 2 per cent of partial melt in the depth range ??? 100-300 km beneath east-central Iceland. The anomalous body is approximately cylindrical in the top 250 km but tabular in shape at greater depth, elongated north-south and generally underlying the spreading plate boundary. Such a morphological change and its relationship to surface rift zones are predicted to occur in convective upwellings driven by basal heating, passive upwelling in response to plate separation and lateral temperature gradients. Although we cannot resolve structure deeper than ??? 450 km, and do not detect a bottom to the anomaly, these models suggest that it extends no deeper than the mantle transition zone. Such models thus suggest a shallow origin for the Iceland hotspot rather than a deep mantle plume, and imply that the hotspot has been located on the spreading ridge in the centre of the north Atlantic for its entire history, and is not fixed relative to other Atlantic hotspots. The results are consistent with recent, regional full-thickness mantle tomography and whole-mantle tomography images that show a strong, low-wave-speed anomaly beneath the Iceland region that is confined to the upper mantle and thus do not require a plume in the lower mantle. Seismic and geochemical observations that are interpreted as indicating a lower mantle, or core-mantle boundary origin for the North Atlantic Igneous Province and the Iceland hotspot should be re-examined to consider whether they are consistent with upper mantle processes.
Geophysical Journal International↗
Three-dimensional seismic structure and moment tensors of non-double-couple earthquakes at the Hengill-Grensdalur volcanic complex, Iceland
The volcanic and geothermal areas of Iceland are rich sources of non-double-couple (non-DC) earthquakes. A state-of-the-art digital seismometer network deployed at the Hengill–Grensdalur volcanic complex in 1991 recorded 4000 small earthquakes. We used the best recorded of these to determine 3-D V P and V P /V S structure tomographically and accurate earthquake moment tensors. The V P field is dominated by high seismic wave speed bodies interpreted as solidified intrusions. A widespread negative (−4 per cent) V P /V S anomaly in the upper 4km correlates with the geothermal field, but is too strong to be caused solely by the effect of temperature upon liquid water or the presence of vapour, and requires in addition mineralogical or lithological differences between the geothermal reservoir and its surroundings. These may be caused by geothermal alteration. Well-constrained moment tensors were obtained for 70 of the best-recorded events by applying linear programming methods to P - and S -wave polarities and amplitude ratios. About 25 per cent of the mechanisms are, within observational error, consistent with DC mechanisms consistent with shear faulting. The other 75 per cent have significantly non-DC mechanisms. Many have substantial explosive components, one has a substantial implosive component, and the deviatoric component of many is strongly non-DC. Many of the non-DC mechanisms are consistent, within observational error, with simultaneous tensile and shear faulting. However, the mechanisms occupy a continuum in source-type parameter space and probably at least one additional source process is occurring. This may be fluid flow into newly formed cracks, causing partial compensation of the volumetric component. Studying non-shear earthquakes such as these has great potential for improving our understanding of geothermal processes and earthquake source processes in general.
Geophysical Journal International↗
Three-dimensional interface modelling with two-dimensional seismic data: the Alpine crust-mantle boundary
We present a new approach to determine the 3-D topography and lateral continuity of seismic interfaces using 2-D-derived controlled-source seismic reflector data. The aim of the approach is to give the simplest possible structure consistent with all reflector data and error estimates. We define simplicity of seismic intrafaces by the degree of interface continuity (ie shortest length of offsets) and by the degree of interface roughness (least surface roughness). The method is applied to structural information of the crust-mantle boundary (Moho) obtained from over 250 controlled-source seismic reflection and refraction profiles in the greater Alpine region. The reflected and refracted phases from the Moho interface and their interpretation regarding crustal thickness are reviewed and their reliability weighted. Weights assigned to each reflector element are transformed to depth errors considering Fresnel volumes. The 2-D-derived reflector elements are relocated in space (3-D migration) and interpolation is performed between the observed reflector elements to obtain continuity of model parameters. Interface offsets are intoduced only where required according to the prinipal of simplicity. The resulting 3-D model of the ALpine crust-mantle boundary shows two offsets that eivide the interface into a European, an Adriatic and a Ligurian Moho, with the European Moho subducting below the Adriatic Moho, and with the Adriatic Moho underthrusting the Ligurian Moho. Each sub-interface depicts the smoothest possible (ie simplest) surface, fitting the reflector data within their assigned errors. The results are consistent with previous studies for those regions with dense and reliable controlled-source seismic data. The newly derived Alpine Moho interface, however, surpasses earlier studies by its lateral extent over an area of about 600km by 600km, by quantifying reliability estimates along the interface, and by obeying the priciple of being consistently as simple as possible.
Geophysical Journal International↗
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 3 He (Sorey et al. 1993) and the emission of about 500 tonnes day −1 of CO 2 (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 2 , 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.
Geophysical Journal International↗