Magma accumulation beneath Santorini volcano, Greece, from P-wave tomography

Geology ◽  
2019 ◽  
Vol 48 (3) ◽  
pp. 231-235 ◽  
Author(s):  
B.G. McVey ◽  
E.E.E. Hooft ◽  
B.A. Heath ◽  
D.R. Toomey ◽  
M. Paulatto ◽  
...  
Keyword(s):  
Three Dimensional ◽  
Seismic Attenuation ◽  
P Wave ◽  
Structural Constraints ◽  
Low Velocity ◽  
Magma Body ◽  
Velocity Anomaly ◽  
Velocity Models ◽  
Santorini Volcano ◽  
Ray Bending

Abstract Despite multidisciplinary evidence for crustal magma accumulation below Santorini volcano, Greece, the structure and melt content of the shallow magmatic system remain poorly constrained. We use three-dimensional (3-D) velocity models from tomographic inversions of active-source seismic P-wave travel times to identify a pronounced low-velocity anomaly (–21%) from 2.8 km to 5 km depth localized below the northern caldera basin. This anomaly is consistent with depth estimates of pre-eruptive storage and a recent inflation episode, supporting the interpretation of a shallow magma body that causes seismic attenuation and ray bending. A suite of synthetic tests shows that the geometry is well recovered while a range of melt contents (4%–13% to fully molten) are allowable. A thin mush region (2%–7% to 3%–10% melt) extends from the main magma body toward the northeast, observed as low velocities confined by tectono-magmatic lineaments. This anomaly terminates northwest of Kolumbo; little to no melt underlies the seamount from 3 to 5 km depth. These structural constraints suggest that crustal extension and edifice loads control the geometry of magma accumulation and emphasize that the shallow crust remains conducive to melt storage shortly after a caldera-forming eruption.

scholarly journals Seismic Structure of Local Crustal Earthquakes beneath the Zipingpu Reservoir of Longmenshan Fault Zone

2011 ◽  
Vol 2011 ◽  
pp. 1-6 ◽  
Author(s):  
Haiou Li ◽  
Xiwei Xu ◽  
Wentao Ma ◽  
Ronghua Xie ◽  
Jingli Yuan ◽  
...  
Keyword(s):  
Fault Zone ◽  
Three Dimensional ◽  
P Wave ◽  
The Tibetan Plateau ◽  
Seismic Structure ◽  
Infiltration Depth ◽  
Low Velocity ◽  
Longmenshan Fault Zone ◽  
Velocity Models ◽  
Longmenshan Fault

Three-dimensional P wave velocity models under the Zipingpu reservoir in Longmenshan fault zone are obtained with a resolution of 2 km in the horizontal direction and 1 km in depth. We used a total of 8589 P wave arrival times from 1014 local earthquakes recorded by both the Zipingpu reservoir network and temporary stations deployed in the area. The 3-D velocity images at shallow depth show the low-velocity regions have strong correlation with the surface trace of the Zipingpu reservoir. According to the extension of those low-velocity regions, the infiltration depth directly from the Zipingpu reservoir itself is limited to 3.5 km depth, while the infiltration depth downwards along the Beichuan-Yingxiu fault in the study area is about 5.5 km depth. Results show the low-velocity region in the east part of the study area is related to the Proterozoic sedimentary rocks. The Guanxian-Anxian fault is well delineated by obvious velocity contrast and may mark the border between the Tibetan Plateau in the west and the Sichuan basin in the east.

Three-dimensional seismic structure of the lithosphere under Montana Lasa

1976 ◽  
Vol 66 (2) ◽  
pp. 501-524
Author(s):  
Keiiti Aki ◽  
Anders Christoffersson ◽  
Eystein S. Husebye
Keyword(s):  
Generalized Inverse ◽  
Random Medium ◽  
Three Dimensional ◽  
P Wave ◽  
Seismic Structure ◽  
Low Velocity ◽  
Medium Theory ◽  
Velocity Anomaly ◽  
Inversion Techniques ◽  
Teleseismic Events

abstract Using P-wave residuals for teleseismic events observed at the Montana Large Aperture Seismic Array (LASA), we have determined the three-dimensional seismic structure of the lithosphere under the array to a depth of 140 km. The root-mean-square velocity fluctuation was found to be at least 3.2 per cent which may be compared to estimate of ca. 2 per cent based on the Chernov random medium theory. The solutions are given by both the generalized inverse and stochastic inverse methods in order to demonstrate the relative merit of different inversion techniques. The most conspicuous feature of the lithosphere under LASA is a low-velocity anomaly in the central and northeast part of the array siting area with the N60°E trend and persisting from the upper crust to depths greater than 100 km. We interpret this low-velocity anomaly as a zone of weakness caused by faulting and shearing associated with the building of the Rocky Mountains.

Seismic velocity models for heat zones in Athabasca tar sands

Geophysics ◽  
1990 ◽  
Vol 55 (8) ◽  
pp. 1108-1112 ◽  
Author(s):  
Larry R. Lines ◽  
Ronald Jackson ◽  
James D. Covey
Keyword(s):  
Seismic Velocity ◽  
Three Dimensional ◽  
Velocity Model ◽  
Steam Injection ◽  
P Wave ◽  
Injection Velocity ◽  
Low Velocity ◽  
Borehole Data ◽  
Tar Sands ◽  
Velocity Models

Recent laboratory and field studies indicate that the P-wave velocity in Athabasca tar sands decreases when temperature increases during steam injection. In this paper we derive time variant velocity models from seismic traveltime inversions of both reflection and borehole data. Prior to steam injection, three‐dimensional (3-D) reflector velocity‐depth models are established using image‐ray conversions of traveltimes to depth. The changes in velocity due to steam injection are modeled by inverting traveltime data from seismic monitor surveys after steam injection and comparing these results to velocities computed prior to steam injection. Velocity models are essentially determined by traveltimes from the 3-D seismic reflection survey. The surface‐to‐wellbore data traveltimes show the expected delay caused by steam injection but do not significantly alter the velocity model produced by reflection traveltimes. For seismic monitor surveys, low‐velocity zones show a very good correlation with zones of temperature increase at injector well positions. The results indicate that velocity models obtained from seismic traveltimes may prove useful in detecting steam fronts in tar sands.

scholarly journals P-wave crustal tomography of Greece with use of an accurate two-point ray tracer

1997 ◽  
Vol 40 (1) ◽  
Author(s):  
G. Drakatos ◽  
G. Karantonis ◽  
G. N. Stavrakakis
Keyword(s):  
Large Scale ◽  
Velocity Structure ◽  
Three Dimensional ◽  
Seismic Energy ◽  
Aegean Sea ◽  
P Wave ◽  
Low Velocity ◽  
Aegean Area ◽  
Arrival Times ◽  
Horizontal Variations

The three-dimensional velocity structure of the crust in the Aegean sea and the surrounding regions (34.0º-42.OºN, 19.0ºE-29.0ºE) is investigated by inversion of about 10000 residuals of arrival times of P-wave from local events. The resulting velocity structure shows strong horizontal variations due to the complicated crustal structure and the variations of crustal thickness. The northern part of the region generally shows high velocities. In the inner part of the volcanic arc (Southern Aegean area), relatively low velocities are observed, suggesting a large-scale absorption of seismic energy as confirmed by the low seismicity of the region. A low velocity zone was observed along the subduction zone of the region, up to a depth of 4 km. The existence of such a zone could be due to granitic or other intrusions in the crust during the uplift of the region during Alpidic orogenesis.

scholarly journals Preliminary Results of Receiver Function Forward Velocity Modelling at Merapi Volcano

2021 ◽  
Vol 873 (1) ◽  
pp. 012056
Author(s):  
M F R Auly ◽  
A K Ilahi ◽  
I Madrinovella ◽  
S Widyanti ◽  
S K Suhardja ◽  
...  
Keyword(s):  
Receiver Function ◽  
Tectonic Setting ◽  
Vertical Component ◽  
P Wave ◽  
Function Technique ◽  
Boundary Current ◽  
Low Velocity ◽  
Mantle Boundary ◽  
Velocity Models ◽  
Synthetic Model

Abstract The tectonic setting of Java island, located at southwestern edge of the Eurasia continent, is dominated by the subduction of Indo-Australia plate. One of the characteristics of active subduction is active seismicity, the generation of arc magmatism and volcanic activity. Mt. Merapi is one example of active volcano related with the subduction process. It is one of the most active volcanoes with location close to high population area. To better understand this area, we employed the Receiver Function technique, a method to image sub surface structure by removing the vertical component from horizontal component. First, we collected high magnitude events and processed RF with water level deconvolution method. Then, we constructed synthetic model with initial velocity input from previous tomography model. Note that we used reflectivity method in generating synthetic model with input parameters matched with parameters from real data processing. Next, we adjusted velocity inputs mainly on tops sediments (1-3 km) to include sediment layers and volcanic rocks, mid-depth low velocity zone that may be related with magma chamber and depth of crust-mantle boundary. Current forward velocity models show a relatively good agreement from 3 stations (ME25, ME32 and ME36). We estimate a thin layer of sediments followed a zone of velocity layer at a depth of 10-15 km and crust-mantle boundary ranging from 26-29 km. In this study, simulated that the signal of sediments layer and low velocity layers interfere main crust mantle boundary that supposed to be highest signal after the P wave in the typical receiver function study.

scholarly journals Crustal velocity structure of the Apennines (Italy) from P-wave travel time tomography

1996 ◽  
Vol 39 (6) ◽  
Author(s):  
C. Chiarabba ◽  
A. Amato
Keyword(s):  
Velocity Structure ◽  
P Wave ◽  
Depth Interval ◽  
Low Velocity ◽  
Data Set ◽  
Arrival Times ◽  
Velocity Anomaly ◽  
Minimum Number ◽  
Velocity Images ◽  
Lower Crustal

In this paper we provide P-wave velocity images of the crust underneath the Apennines (Italy), focusing on the lower crustal structure and the Moho topography. We inverted P-wave arrival times of earthquakes which occurred from 1986 to 1993 within the Apenninic area. To overcome inversion instabilities due to noisy data (we used bulletin data) we decided to resolve a minimum number of velocity parameters, inverting for only two layers in the crust and one in the uppermost mantle underneath the Moho. A partial inversion of only 55% of the overall dataset yields velocity images similar to those obtained with the whole data set, indicating that the depicted tomograms are stable and fairly insensitive to the number of data used. We find a low-velocity anomaly in the lower crust extending underneath the whole Apenninic belt. This feature is segmented by a relative high-velocity zone in correspondence with the Ortona-Roccamonfina line, that separates the northern from the southern Apenninic arcs. The Moho has a variable depth in the study area, and is deeper (more than 37 km) in the Adriatic side of the Northern Apennines with respect to the Tyrrhenian side, where it is found in the depth interval 22-34 km.

scholarly journals Velocity-attenuation model from check-shot drift trends in North Sea well data

Geophysics ◽  
2020 ◽  
Vol 85 (2) ◽  
pp. D65-D74 ◽  
Author(s):  
Andrew J. Carter ◽  
Veronica A. Torres Caceres ◽  
Kenneth Duffaut ◽  
Alexey Stovas
Keyword(s):  
North Sea ◽  
Velocity Dispersion ◽  
Seismic Attenuation ◽  
P Wave ◽  
Linear Wave ◽  
Velocity Models ◽  
Velocity Attenuation ◽  
Attenuation Models ◽  
Well Data ◽  
Attenuation Values

Seismic attenuation distorts phase and narrows bandwidth in seismic surveys. It is also an exploration attribute, as, for example, gas or overpressure, may create attenuation anomalies. Compensating attenuation in imaging requires accurate models. Detailed attenuation models may be obtained using full-waveform inversion (FWI) or attenuation tomography, but their accuracy benefits from reliable starting models and/or constraints. Seismic attenuation and velocity dispersion are necessarily linked for causal linear wave propagation such that higher frequencies travel faster than lower frequencies in an attenuative medium. In publicly released well data from the Norwegian North Sea, we have observed systematic positive linear trends in check-shot drift when comparing (lower frequency) time-depth curves with (higher frequency) integrated sonic transit times. We observe velocity dispersion consistent with layers having constant seismic attenuation. Adapting a previously published method, and assuming an attenuation-dispersion relationship, we use drift gradients, measured over thick stratigraphic units, to estimate interval P-wave attenuation and tentatively interpret its variation in terms of porosity and fluid mobility. Reflectivity modeling predicts a very low attenuation contribution from peg-leg multiples. We use the attenuation values to develop a simple regional relationship between P-wave velocity and attenuation. Observed low drift gradients in some shallower units lead to an arch-shaped model that predicts low attenuation at both low and high velocities. The attenuation estimates were broadly comparable with published effective attenuation values obtained independently nearby. This general methodology for quickly deriving a regional velocity-attenuation relationship could be used anywhere that coincident velocity models are available at seismic and sonic frequencies. Such relationships can be used for fast derivation (from velocities) of starting attenuation models for FWI or tomography, constraining or linking velocity and attenuation in inversion, deriving models for attenuation compensation in time processing, or deriving background trends in screening for attenuation anomalies in exploration.

Indication of active Faults in the New Madrid Seismic Zone from Precise Location of Hypocenters

1988 ◽  
Vol 59 (4) ◽  
pp. 123-131 ◽  
Author(s):  
L. Himes ◽  
W. Stauder ◽  
R. B. Herrmann
Keyword(s):  
Active Faults ◽  
P Wave ◽  
New Madrid Seismic Zone ◽  
Supporting Evidence ◽  
Seismic Zone ◽  
Low Velocity ◽  
Focal Mechanism Solutions ◽  
Velocity Models ◽  
Planar Surfaces ◽  
New Madrid

Abstract The hypocenter locations of the larger and better recorded earthquakes of the New Madrid seismic zone are examined in order to determine how closely the hypocenters lie along planar surfaces, thus relating the foci to active fault surfaces. For this purpose more than 500 earthquakes of the region have been selected for study, based on the number (7 or more) of observing stations used in the initial hypocenter location and on the quality of the P-wave onset. These events are relocated using a joint hypocenter-velocity-depth (JHVD) algorithm. The relocated earthquakes are separated geographically into three trends: ARK, the southwest trending zone from Caruthersville, Missouri, to Marked Tree, Arkansas; DWM, the northeast trending zone from New Madrid to Charleston, Missouri; and CEN, the central, left-stepping offset zone from Ridgely, Tennessee, to New Madrid, Missouri. Vertical profiles taken along and across the ARK and DWM trends verify the strike and dip of dominantly strike slip motion on near vertical active faults along these trends. These results agree with previously determined composite focal mechanism solutions for these trends. No coherent picture has been obtained for the CEN trend. As a by-product of the study, velocity models from the JHVD inversion are found to be reasonably uniform throughout the New Madrid seismic zone, and to offer supporting evidence for the presence of a shallow low velocity zone in the central portion of the Mississippi embayment.

scholarly journals The Local Earthquake Tomography of Erzurum (Turkey) Geothermal Area

2019 ◽  
Vol 23 (3) ◽  
pp. 209-223 ◽  
Author(s):  
Caglar Ozer ◽  
Mehmet Ozyazicioglu
Keyword(s):  
Wave Velocity ◽  
Seismic Velocity ◽  
Three Dimensional ◽  
P Wave ◽  
Geothermal Area ◽  
S Wave ◽  
Velocity Models ◽  
Local Earthquake Tomography ◽  
P Wave Velocity ◽  
Local Earthquake

Erzurum and its surroundings are one of the seismically active and hydrothermal areas in the Eastern part of Turkey. This study is the first approach to characterize the crust by seismic features by using the local earthquake tomography method. The earthquake source location and the three dimensional seismic velocity structures are solved simultaneously by an iterative tomographic algorithm, LOTOS-12. Data from a combined permanent network comprising comprises of 59 seismometers which was installed by Ataturk University-Earthquake Research Center and Earthquake Department of the Disaster and Emergency Management Authority  to monitor the seismic activity in the Eastern Anatolia, In this paper, three-dimensional Vp and Vp/Vs characteristics of Erzurum geothermal area were investigated down to 30 km by using 1685 well-located earthquakes with 29.894 arrival times, consisting of 17.298 P- wave and 12.596 S- wave arrivals. We develop new high-resolution depth-cross sections through Erzurum and its surroundings to provide the subsurface geological structure of seismogenic layers and geothermal areas. We applied various size horizontal and vertical checkerboard resolution tests to determine the quality of our inversion process. The basin models are traceable down to 3 km depth, in terms of P-wave velocity models. The higher P-wave velocity areas in surface layers are related to the metamorphic and magmatic compact materials. We report that the low Vp and high Vp/Vs values are observed in Yedisu, Kaynarpinar, Askale, Cimenozu, Kaplica, Ovacik, Yigitler, E part of Icmeler, Koprukoy, Uzunahmet, Budakli, Soylemez, Koprukoy, Gunduzu, Karayazi, Icmesu, E part of Horasan and Kaynak regions indicated geothermal reservoir.

Sign in / Sign up

Export Citation Format

Share Document