Description
This dataset contains a three-dimensional seismic slowness model and MATLAB processing code derived from P-wave seismic data collected at Newberry Volcano in central Oregon. The associated study combines active-source first-arrival travel times with teleseismic P-wave delay times to jointly image upper-crustal velocity structure beneath the volcano (targeting the shallow crust to ~10 km depth).
The included workflow uses a one-dimensional velocity model and a three-dimensional slowness model to compute a fractional velocity perturbation (dv/v) volume and to plot horizontal slices and vertical sections. The script reads the 1D model, interpolates it onto the depth grid of the 3D model, converts the 1D velocity profile to slowness, and then forms dv/v by comparing the 1D slowness volume to the 3D slowness model before plotting results. A python version of the workflow is also included, which was used to create the provided 3D model grid CSV files.
These data have been used by and uploaded to the Geothermal Data Repository (GDR) as part of the DE-risking Exploration of geothermal Plays in magmatic ENvironments (DEEPEN) project. The data were originally produced through the 2008 Newberry Seismic Experiment, funded by a USGS Venture Capital grant and NSF grants (EAR-207670 and EAR-207671). Seismic data from the 2008 experiment and from earlier 1983 and 1984 USGS experiments are archived at the IRIS DMC. Please refer to the associated manuscript (linked below) for additional details on the original data collection and research acknowledgements.
Resources
| Name | Format | Description | Link |
|---|---|---|---|
| Associated Publication | HTML | Research article "Imaging the magmatic system of Newberry Volcano using joint active source and teleseismic tomography" from Heath et al. (2015). The article describes a joint inversion of active-source travel times and teleseismic delay times at Newberry Volcano, including acquisition context, inversion approach, and interpretation of recovered velocity anomalies. | Open resource |
| 3D Grid Clipped.csv | CSV | This CSV contains a clipped subset of the 3D seismic model. The dataset includes P-wave slowness values and derived fractional P-wave velocity perturbations (dv/v). Please refer to the included README for information on how this file was created. | Open resource |
| 1D Reference Velocity Model.csv | CSV | One-dimensional velocity model used as the reference profile in the processing/anomaly workflow. | Open resource |
| CSV_README.txt | Text | README describing the derivation of, and each column in the "3D Grid Clipped" and "3D Grid Full" CSV files. | Open resource |
| MATLAB Processing Script.m | Unknown | MATLAB script that loads the 3D slowness model and the 1D velocity model, constructs a 1D slowness volume on the 3D grid, computes dv/v, and generates map-view slices and vertical cross-sections. | Open resource |
| Python Processing Script.ipynb | Python | Translated workflow to Python from "MATLAB Processing Script.m", with added CSV exporting and XY clipping logic. Used to create both 3D grid CSV files in this dataset. | Open resource |
| 3D Grid Full.csv | CSV | This CSV contains the full 3D seismic model. The dataset includes P-wave slowness values and derived fractional P-wave velocity perturbations (dv/v). Please refer to the included README for information on how this file was created. | Open resource |
| 3D Slowness Model.mat | MAT | MATLAB data file containing the 3D slowness model and its associated grids/dimensions used by the processing workflow. | Open resource |