As part of the geophysical characterization suite for the first EGS Collab tesbed, here are the baseline cross-well seismic data and resultant models. The campaign seismic data have been organized, concatenated with geometry and compressional (P-) & and shear (S-) wave picks, and submitted as SGY files. P-wave data were collected and analyzed in both 2D and 3D, while S-wave data were collected and analyzed in 2D only. Inversion models are provided as point volumes; the volumes have been culled to include only the points within source/receiver array coverage. The full models space volumes are also included, if relevant. An AGU 2018 poster by Linneman et al. is included that provides visualizations/descriptions of the cross-well seismic characterization method, elastic moduli calculations, and images of model inversion results.

This data submission includes the raw time-lapse ERT (electrical resistivity tomography) monitoring data, flow system data, operator logs, E4D (https://e4d.pnnl.gov) inversion files, and metadata necessary to reproduce the 4D ERT inversion for the Oct. 24 through Nov. 7 2018 post-stimulation flow test in test bed 1. The tests were done at the Sanford Underground Research Facility at Homestake Mine in South Dakota.

The link points to a website at NCEDC to download the full moment tensors inversion software The moment tensor analysis conducted in the current project is based on the full moment tensor model described in Minson and Dreger (2008). The software including source, examples and tutorial can be obtained from ftp://ncedc.org/outgoing/dreger (download file pasi-nov282012.tar.gz). Performance criteria, mathematics and test results are provided by Minson and Dreger (2008), Ford et al. (2008, 2009, 2010, 2012) and Saikia (1994). References: Ford, S., D. Dreger and W. Walter (2008). Source Characterization of the August 6, 2007 Crandall Canyon Mine Seismic Event in Central Utah, Seism. Res. Lett., 79, 637-644. Ford, S. R., D. S. Dreger and W. R. Walter (2009). Identifying isotropic events using a regional moment tensor inversion, J. Geophys. Res., 114, B01306, doi:10.1029/2008JB005743. Ford, S. R., D. S. Dreger and W. R. Walter (2010). Network sensitivity solutions for regional moment tensor inversions, Bull. Seism. Soc. Am., 100, p. 1962-1970. Ford, S. R., W. R. Walter, and D. S. Dreger (2012). Event discrimination using regional moment 665 tensors with teleseismic-P constraints, Bull. Seism. Soc. Am. 102, 867-872. Minson, S. and D. Dreger (2008), Stable Inversions for Complete Moment Tensors, Geophys. J. Int., 174, 585-592. Saikia, C.K. (1994), Modified Frequency-Wavenumber Algorithm for Regional Seismograms using Filons Quadrature: Modeling of Lg Waves in Eastern North America. Geophys. J. Int., 118, 142-158.

Graph theory is useful for estimating time-dependent model parameters via weighted least-squares using interferometric synthetic aperture radar (InSAR) data. Plotting acquisition dates (epochs) as vertices and pair-wise interferometric combinations as edges defines an incidence graph. The edge-vertex incidence matrix and the normalized edge Laplacian matrix are factors in the covariance matrix for the pair-wise data. Using empirical measures of residual scatter in the pair-wise observations, we estimate the variance at each epoch by inverting the covariance of the pair-wise data. We evaluate the rank deficiency of the corresponding least-squares problem via the edge-vertex incidence matrix. We implement our method in a MATLAB software package called GraphTreeTA available on GitHub (https://github.com/feigl/gipht). We apply temporal adjustment to the data set described in Lu et al. (2005) at Okmok volcano, Alaska, which erupted most recently in 1997 and 2008. The data set contains 44 differential volumetric changes and uncertainties estimated from interferograms between 1997 and 2004. Estimates show that approximately half of the magma volume lost during the 1997 eruption was recovered by the summer of 2003. Between June 2002 and September 2003, the estimated rate of volumetric increase is (6.2 +/- 0.6) x 10^6 m^3/yr. Our preferred model provides a reasonable fit that is compatible with viscoelastic relaxation in the five years following the 1997 eruption. Although we demonstrate the approach using volumetric rates of change, our formulation in terms of incidence graphs applies to any quantity derived from pair-wise differences, such as wrapped phase or wrapped residuals. Date of final oral examination: 05/19/2016 This thesis is approved by the following members of the Final Oral Committee: Kurt L. Feigl, Professor, Geoscience Michael Cardiff, Assistant Professor, Geoscience Clifford H. Thurber, Vilas Distinguished Professor, Geoscience

Using LiDAR to identify structural and volcanic evolution of a Miocene-Pleistocene age bimodal volcanic complex and implications for geothermal potential. The file includes an updated geologic map, methods, and preliminary results.

The PoroTomo team has completed inverse modeling of the three data sets (seismology, geodesy, and hydrology) individually, as described previously. The estimated values of the material properties are registered on a three-dimensional grid with a spacing of 25 meters between nodes. The material properties are listed an Excel file. Figures show planar slices in three sets: horizontal slices in a planes normal to the vertical Z axis (Z normal), vertical slices in planes perpendicular to the dominant strike of the fault system (X normal), and vertical slices in planes parallel to the dominant strike of the fault system (Y normal). The results agree on the following points. The material is unconsolidated and/or fractured, especially in the shallow layers. The structural trends follow the fault system in strike and dip. The geodetic measurements favor the hypothesis of thermal contraction. Temporal changes in pressure, subsidence rate, and seismic amplitude are associated with changes in pumping rates during the four stages of the deployment in 2016. The modeled hydraulic conductivity is high in fault damage zones. All the observations are consistent with the conceptual model: highly permeable conduits along faults channel fluids from shallow aquifers to the deep geothermal reservoir tapped by the production wells.

Multiple data sets converge on a geologic structural model for Glass Buttes, Oregon geothermal prospect, Patrick Walsh, Brigette Martini, Chet Lide, Darrick Boschmann, John DIlles, Andrew Meigs, 2010 Ormat Nevada, Zonge Geophysical, Oregon State University American Geophysical Union, Poster Session

These data are Pacific Northwest National Lab inversions of an amalgamation of two surface gravity datasets: Davenport-Newberry gravity collected prior to 2012 stimulations and Zonge International gravity collected for the project "Novel use of 4D Monitoring Techniques to Improve Reservoir Longevity and Productivity in Enhanced Geothermal Systems" in 2012. Inversions of surface gravity recover a 3D distribution of density contrast from which intrusive igneous bodies are identified. The data indicate a body name, body type, point type, UTM X and Y coordinates, Z data is specified as meters below sea level (negative values then indicate elevations above sea level), thickness of the body in meters, suscept, density anomaly in g/cc, background density in g/cc, and density in g/cc. The model was created using a commercial gravity inversion software called ModelVision 12.0 (http://www.tensor-research.com.au/Geophysical-Products/ModelVision). The initial model is based on the seismic tomography interpretation (Beachly et al., 2012). All the gravity data used to constrain this model are on the GDR: https://gdr.openei.org/submissions/760.

This submission contains 167 deviatoric moment tensor (MT) solutions for the seismicity observed two years prior and three years post start of injection activities at The Geysers Prati 32 EGS Demonstration. Also included is a statistical representation of the properties of 751 fractures derived from the analysis of seismicity observed two years prior and three years post start of injection activities at The Geysers Prati 32 EGS Demonstration Project. The locations of the fractures are taken from microseismic hypocenters, the fracture areas are derived from moment magnitudes via scaling relationships, and the azimuths (sigma 1) and dips (sigma 3) are derived from the results of stress analyses.

This submission contains 167 full moment tensor (MT) solutions for the seismicity observed two years prior and three years post start of injection activities. Also included are the azimuth and plunge angles for the three main stress directions sigma1, sigma 2 and sigma 3 at the Prati32 EGS demonstration site in the northwest Geysers geothermal reservoir. The data are divided into 15 time periods spanning a range of five years, including two years prior to start of injection until three years post start of injection activities.

The site characterization data used to develop the conceptual geologic model for the Snake River Plain site in Idaho, as part of phase 1 of the Frontier Observatory for Research in Geothermal Energy (FORGE) initiative. This collection includes data on seismic events, groundwater, geomechanical models, gravity surveys, magnetics, resistivity, magnetotellurics (MT), rock physics, stress, the geologic setting, and supporting documentation, including several papers. Also included are 3D models (Petrel and Jewelsuite) of the proposed site. Data for wells INEL-1, WO-2, and USGS-142 have been included as links to separate data collections. These data have been assembled by the Snake River Geothermal Consortium (SRGC), a team of collaborators that includes members from national laboratories, universities, industry, and federal agencies, lead by the Idaho National Laboratory (INL). Other contributors include the National Renewable Energy Laboratory (NREL), Lawrence Livermore National Laboratory (LLNL), the Center for Advanced Energy Studies (CEAS), the University of Idaho, Idaho State University, Boise State University, University of Wyoming, University of Oklahoma, Energy and Geoscience Institute-University of Utah, US Geothermal, Baker Hughes Campbell Scientific Inc., Chena Power, US Geological Survey (USGS), Idaho Department of Water Resources, Idaho Geological Survey, and Mink GeoHydro.

These resources describe the 3D geophysical inversion modeling of gravity data at the FORGE site near Milford, Utah. FORGE is the Frontier Observatory for Research in Geothermal Energy and the site in Utah has been selected by the U.S. Dept. of Energy for a 5-year R&D program to test technologies for the development of Engineered Geothermal Systems (EGS). 3D modelling of gravity data at the FORGE site is to help characterize the subsurface geologic framework. Specifically, modelling of gravity data in 3D, used in conjunction with rock density measurements and other subsurface geologic information can provide an independent test of an existing 3D geologic model (e.g. Witter et al., 2018). Such an exercise can be useful for reducing uncertainty in 3D geologic models (Witter et al, 2019).

This is the 3D version of the MT data for the West Flank Coso FORGE area. The Coso geothermal field has had three Magnetotelluric (MT) datasets collected including surveys in 2003, 2006, and 2011. The final collection, in 2011, expanded the survey to the west and covers the West Flank of FORGE area.This most recent data set was collected by Schlumberger/WesternGeco and inverted by the WesternGeco GeoSolutions Integrated EM Center of Excellence in Milan, Italy; the 2003 and 2006 data were integrated for these inversions in the present study.