Development of a multi-scale seismic imaging method and its application to Newberry volcano
Abstract
While some shallow magma reservoirs beneath land volcanoes have been seismically imaged, our understanding of their mid-to lower-crustal magma plumbing systems remains limited. Active source techniques allow for optimizing source geometry to image shallow crustal magmatic systems, but typically do not resolve depths greater than ~6-10 km, since the velocity gradient in the mid-to lower crust limits the energy that turns at these depths. For closely spaced seismic stations, teleseismic data provide constraints on deeper crustal heterogeneity, particularly its lateral variation. On the other hand, teleseismic studies do not resolve vertical variations in crustal structure well because their ray paths are almost vertical, which results in limited crossing rays within the crust. Furthermore, due to the generally lower frequency content of teleseismic arrivals, these data primarily image longer wavelength heterogeneity. Here we develop a multi-scale seismic imaging method that combines high-frequency active source data with lower frequency teleseismic data, and test its usefulness for imaging the crustal structure beneath Newberry Volcano in central Oregon. To constrain the P wave structure throughout the crust, we develop an iterative method that includes 3-D sensitivity kernels and 3-D raytracing. The use of sensitivity kernels provides a physically motivated method that accounts for the different resolving capability of teleseismic and active source data. Various approximations to the sensitivity kernel are explored, with the goal of maximizing both computational efficiency and inversion accuracy. 3-D raytracing determines an accurate ray path, which is important in the highly heterogeneous crust. We verify our approach with synthetic data calculated using finite difference waveform modeling. We discuss our inversion results and our ability to resolve the velocity structure throughout the crust. We will investigate the presence of high and low velocity regions that may be associated with cooled intrusives and the magma plumbing system, respectively. We will test whether our approach can resolve the deeper magma plumbing system and its relation to a previously imaged magma reservoir at 3-6 km depth (Beachly et al., 2012).
- Publication:
-
AGU Fall Meeting Abstracts
- Pub Date:
- December 2013
- Bibcode:
- 2013AGUFM.S23A2477H
- Keywords:
-
- 7270 SEISMOLOGY Tomography