Microseismic event catalog and 3-D P-wave velocity (Vp) model at San Emidio, Nevada, developed using seismic data collected from San Emidio from April to May 2022 as part of the WHOLESCALE project. SanEmidio_2022MicroseismicEventCatalog.csv:This is seismic event catalog, including following fields: Origin Time: In UTC    Magnitude: coda duration magnitude. "NaN" signifies not a number.    Latitude: Relocated latitude in decimal degrees on WGS84Longitude: Relocated longitude in decimal degrees on WGS84E(m): UTM zone 11 easting in meters N(m): UTM zone 11 northing in metersH(m): Elevation above sea level (geoid) in metersUncertainty_E(m): uncertainty in meters of estimated location in easting Uncertainty_N(m): uncertainty in meters of estimated location in northing Uncertainty_H(m): uncertainty in meters of estimated location in elevation SanEmidio_2022Vpmodel.csv:This is 3-D tomographic Vp model, including location (latitude, longitude, E, N, H) and velocity (Vp) values at each tomographic node. The last column represents the model resolvability, estimated from the checkerboard resolution test. Regions with resolvability values above 0.7 are considered well-resolved. Details on the development of the event catalog and Vp model are provided in a manuscript by Hao Guo et al., titled “Microseismicity Modulation Due to Changes in Geothermal Production at San Emidio, Nevada, USA,” which is currently under revision at Geophysical Research Letters.

Microseismic event catalog and 3-D P-wave velocity (Vp) model at San Emidio, Nevada, developed using seismic data collected from San Emidio from April to May 2022 as part of the WHOLESCALE project. SanEmidio_2022MicroseismicEventCatalog.csv:This is seismic event catalog, including following fields: Origin Time: In UTC    Magnitude: coda duration magnitude. "NaN" signifies not a number.    Latitude: Relocated latitude in decimal degrees on WGS84Longitude: Relocated longitude in decimal degrees on WGS84E(m): UTM zone 11 easting in meters N(m): UTM zone 11 northing in metersH(m): Elevation above sea level (geoid) in metersUncertainty_E(m): uncertainty in meters of estimated location in easting Uncertainty_N(m): uncertainty in meters of estimated location in northing Uncertainty_H(m): uncertainty in meters of estimated location in elevation SanEmidio_2022Vpmodel.csv:This is 3-D tomographic Vp model, including location (latitude, longitude, E, N, H) and velocity (Vp) values at each tomographic node. The last column represents the model resolvability, estimated from the checkerboard resolution test. Regions with resolvability values above 0.7 are considered well-resolved. Details on the development of the event catalog and Vp model are provided in a manuscript by Hao Guo et al., titled “Microseismicity Modulation Due to Changes in Geothermal Production at San Emidio, Nevada, USA,” which is currently under revision at Geophysical Research Letters.

This submission includes a high-precision seismic event catalog estimated from seismic data collected at San Emidio, Nevada from April to May 2022. The catalog lists the precise time, location, and magnitude of microseismic events recorded during this period. Both the seismic data and microseismic event catalog were produced as part of the Water & Hole Observations Leverage Effective Stress Calculations and Lessen Expenses (WHOLESCALE) project. Attached here are the microseismic event catalog, a link to the raw seismic data, and a detailed description of methods used to create the catalog.

This is a collection of metadata from the 2021 seismic survey at the San Emidio geothermal field in Nevada. In April and May 2021, 37 tri-axial short period seismographs were deployed in a 1.8km diameter cluster centered on 40.367278 deg N, 119.409019 deg W. The first data record started at 2021-04-06T07:09:10Z UTC and the last record ended 2021-05-11T02:58:52Z UTC. The 37 stations included 29 SmartSolo IGU-16HR 3C all-in-one 5 Hz seismographs and 8 DataCube seismographs with 4.5 Hz HGS HG-6(B coil) tri-axial geophones. The raw format (level 0) data includes 353 GB of 500 sps SmartSolo data in native DLD format, 113 GB of 400 sps DataCube data in native DataCube format, and 3.4 GB of GPS data collected during the RTK GPS survey. The SAC data (level 1) files include 564 GB of hourly SAC files. The experiment was designed to monitor seismic activity before, during, and after the planned three-day plant maintenance shutdown April 19-21, 2021. The pumping stopped at 2021/04/19 12:51:45 UTC and resumed about 2021/04/21 21:00:00 UTC. The raw and processed data are in an associated GDR submission, linked below. The metadata here includes files containing experiment details, station locations, seismic data logger specifications, instrumentation specifications, and descriptions of data. Also included are data and metadata from a 2016 seismic survey at the same site.

This dataset includes raw and processed seismic data from the 2021 seismic survey at the San Emidio geothermal field in Nevada. In April and May 2021, 37 tri-axial short period seismographs were deployed in a 1.8km diameter cluster centered on 40.367278 N, 119.409019 W. The first data record started at 2021-04-06T07:09:10Z UTC and the last record ended at 2021-05-11T02:58:52Z UTC. The 37 stations included 29 SmartSolo IGU-16HR 3C all-in-one 5 Hz seismographs and 8 DataCube seismographs with 4.5 Hz HGS HG-6(B coil) tri-axial geophones. The raw format (level 0) data includes 353 GB of 500 sps SmartSolo data in native DLD format, 113 GB of 400 sps DataCube data in native DataCube format, and 3.4 GB of GPS data collected during the RTK GPS survey. The SAC data (level 1) files include 564 GB of hourly SAC files. The experiment was designed to monitor seismic activity before, during, and after the planned three-day plant maintenance shutdown April 19-21, 2021. The pumping stopped at 2021/04/19 12:51:45 UTC and resumed about 2021/04/21 21:00:00 UTC. The dataset is hosted in an AWS data lake, with an associated GDR data set providing the metadata. Links to both of these resources are included below. Additionally, this collection features data and metadata from a 2016 seismic survey at the same site for comparative analysis.

Seismic velocity model AKEP2020 uses earthquake travel-time and ambient noise group velocity data to update Eberhart-Phillips et al. (2006: AKEP2006) The model is provided in a table: vlAKEP2020xyzltlnSFDRE.tbl.txt and in the simul output from velocity inversion: vlAKEP2020.out.txt Map plots of Vp and Vp/Vs are also provided, with lines denoting limits of adequate data. Velocity Inversion Procedure Notes for AKEP2020 model The 2006 AK model and the 2020 updated model both use Transverse Mercator coordinate transformation with central meridian= -150 and counterclockwise rotation of 162.1. Earth-flattening transformation is used for velocity during ray-tracing. The depths are relative to sea-level and station elevations are used. For the group-velocity data, the surface is taken as the 30-km median filtered topography. Velocity with the 3D gridded model is defined by linearly interpolating between nodes. A gradational inversion approach was used with earthquake and shot travel-time data, and group velocity observations for periods 6-15 s. The table provides, from the computed resolution matrix, the diagonal resolution element (DRE), and the spread function (SF), for each Vp and Vp/Vs node. This material is based upon work supported by the U.S. Geological Survey under Grant No. G19AP00019. Note that an earlier model, AKEP2018, from the first year of this funded project was reported in Eberhart-Phillips et al. (2019).

The views and conclusions contained in this document are those of the authors and should not be interpreted as representing the opinions or policies of the U.S. Geological Survey. Mention of trade names or commercial products does not constitute their endorsement by the U.S. Geological Survey

Seismic velocity model AKEP2020 uses earthquake travel-time and ambient noise group velocity data to update Eberhart-Phillips et al. (2006: AKEP2006) The model is provided in a table: vlAKEP2020xyzltlnSFDRE.tbl.txt and in the simul output from velocity inversion: vlAKEP2020.out.txt Map plots of Vp and Vp/Vs are also provided, with lines denoting limits of adequate data. Velocity Inversion Procedure Notes for AKEP2020 model The 2006 AK model and the 2020 updated model both use Transverse Mercator coordinate transformation with central meridian= -150 and counterclockwise rotation of 162.1. Earth-flattening transformation is used for velocity during ray-tracing. The depths are relative to sea-level and station elevations are used. For the group-velocity data, the surface is taken as the 30-km median filtered topography. Velocity with the 3D gridded model is defined by linearly interpolating between nodes. A gradational inversion approach was used with earthquake and shot travel-time data, and group velocity observations for periods 6-15 s. The table provides, from the computed resolution matrix, the diagonal resolution element (DRE), and the spread function (SF), for each Vp and Vp/Vs node. This material is based upon work supported by the U.S. Geological Survey under Grant No. G19AP00019. Note that an earlier model, AKEP2018, from the first year of this funded project was reported in Eberhart-Phillips et al. (2019).

The views and conclusions contained in this document are those of the authors and should not be interpreted as representing the opinions or policies of the U.S. Geological Survey. Mention of trade names or commercial products does not constitute their endorsement by the U.S. Geological Survey

This package contains a 3D Seismic velocity model and an updated microseismic catalog obtained for a double-difference seismic tomography study. The 3D_seismic_velocity_model text file contains x (m), y(m), z(m), P-wave velocity (km/s), P-wave velocity quality indicator (1 for well-constrained; 0 for poorly constrained), S-wave velocity (km/s), and S-wave velocity quality indicator (1 for well-constrained; 0 for poorly constrained). The Updated_MEQ_catalog text file contains event origin time, x(m), y(m), z(m), error in x (m), error in y (m), error in z (m), and RMS misfit (millisecond). The 3D_seismic_P-wave_velocity_model animation file shows slices of the 3D P-wave velocity model. The 3D_seismic_S-wave_velocity_model animation file shows slices of the 3D S-wave velocity model. The Interactive_MEQ_location_comparison API file is an interactive visualization of the updated microseismic event locations and the original seismic catalog. The visualization allows users to view and compare the event locations by dragging, rotating, and zooming in. An updated version of 3D_seismic_velocity_model and associated animations were included, which were calculated with a more strict assumption for quality indicators.

Relocated Earthquake Catalog for Eastern Indonesia (April 2009 to November 2018). Please refer to: Supendi, P., Nugraha, A.D., Widiyantoro, S., Pesicek, J.D., Thurber, C.H., Abdullah, C.I., Daryono, D., Wiyono, S.H., Shiddiqi, H.A., and Rosalia, S. (2020). Relocated aftershocks and background seismicity in Eastern Indonesia shed light on the 2018 Lombok and Palu earthquake sequences, Geophysical Journal International, ggaa118, https://doi.org/10.1093/gji/ggaa118

Relocated Earthquake Catalog for Eastern Indonesia (April 2009 to November 2018). Please refer to: Supendi, P., Nugraha, A.D., Widiyantoro, S., Pesicek, J.D., Thurber, C.H., Abdullah, C.I., Daryono, D., Wiyono, S.H., Shiddiqi, H.A., and Rosalia, S. (2020). Relocated aftershocks and background seismicity in Eastern Indonesia shed light on the 2018 Lombok and Palu earthquake sequences, Geophysical Journal International, ggaa118, https://doi.org/10.1093/gji/ggaa118