Data set containing results from constant mean stress - constant Lode angle true triaxial compression tests performed on Castlegate Sandstone. From the test preformed, the bedding plane and the strain type inside the band of sedimentary rocks can be related to stress histories. The goal of these tests are to understand the conditions that lead to localized deformation in porous sandstone which has geotechnical applications such as oil and natural gas production, carbon dioxide sequestration, and hazardous waste storage.

Includes test data from Confined Brazilian Tests. The tests were aim at investigating the rock failure envelope of four different lithologies including Indiana limestone, Scioto sandstone, Grey Berea sandstone, and Tennessee sandstone. Tables include the strength data obtained using load and damage controlled Brazilian tests on the six lithologies. Confine Brazilian tests are preformed by wrapping a sample in a copper jacket so that the confining fluid does not directly interact with the sample. Then a constant confining pressure is applied to the fluid to create a state of triaxial stress in the disc such that the disc center is under three nonzero and unequal principal stress components. By increasing the confining pressure, the least principal stress changes from tensile to compressive so that rock failure can be investigated over a wide range of stress conditions.

Whole and sidewall core inventories, analyses, photos, and thin section photos for Lawnichak 9-33 in Dover 33 reef, Chester 8-16 and Chester 6-16 in Chester 16 reef.

Directional Cooling-Induced Fracturing (DCIF) experiments were conducted on three rectangular Westerly granite blocks (width=depth=4.0", height=2.0") which were preheated to 200, 400, and 600 degree C to induce damage (microcracks) with varying degrees. Liquid nitrogen was poured in a small, 1"-diameter copper cup attached to the top of the sample, and the resulting acoustic emissions (AEs) and temperature changes on the surface of the sample were monitored. The experiments were conducted under one selected biaxial stress (5.8MPa). The obtained AEs were used to determine the microcracking source locations and amplitude, and the associated moment tensors. The onset time of the AEs was correlated with the cooling temperature, which was used to show that the temperature at the onset of microcracking is not affected significantly by the preexisting damage, compared to the impact of the stress in the sample. Included in this submission are the animations of the AE locations and graphics displaying the measured temperature-AE activity changes for samples with different degrees of microcrack damage.

Directional Cooling-Induced Fracturing (DCIF) experiments were conducted on rectangular Westerly granite blocks (width=depth=4.0", height=2.0"). Liquid nitrogen was poured in a small, 1"-diameter copper cup attached to the top of the sample, and the resulting acoustic emissions (AEs) and temperature changes on the surface of the sample were monitored. Several confining stresses were applied bi-axially to the sides of the samples so that the onset of AE activity and the stress applied to the sample were correlated. The obtained AEs were used to determine the microcracking source locations and amplitude, and the associated moment tensors. Included in this submission are the animations of the AE locations and graphics displaying the measured temperature-AE activity changes for different stresses.

Directional Cooling-Induced Fracturing (DCIF) experiments were conducted on a short, cylindrical sample of Westerly granite (diameter = 4 inches, height ~ 2 inches). Liquid nitrogen was poured in a copper cup attached to the top of the sample, and the resulting acoustic emissions (AEs) and temperature changes on the surface of the sample were monitored. The obtained AEs were used to determine the microcracking source locations and amplitude, and the associated moment tensors. Included in this submission is an animation of the AEs, a graphic displaying the temperature changes, and the measured data.

Data includes Directional Cooling-Induced Fracturing (DCIF) testing data using westerly granite blocks. This submission includes data from two samples of westerly granite, lab sample 7 and 8. Files contain stress, temperature and acoustic emission data acquired during polyaxial, laboratory testing of westerly granite blocks for each sample. FILES: .tradb -- files containing acoustic emission waveforms; sqlite3 database .pridb -- files containing basic acoustic emission information (no waveforms); sqlite3 database .geom -- geometry of the AE sensor network (ASCII)

Characterizing the stimulation mode of a fracture is critical to assess the hydraulic efficiency and the seismic risk related to deep fluid manipulations. We have monitored the three-dimensional displacements of a fluid-driven fracture during water injections in a borehole at ~1.5 km depth in the crystalline rock of the Sanford Underground Research Facility (USA). The fracture initiates at 61% of the minimum horizontal stress by micro-shearing of the borehole on a foliation plane. As the fluid pressure increases further, borehole axial and radial displacements increase with injection time highlighting the opening and sliding of a new hydrofracture growing ~10 m away from the borehole, in accordance with the ambient normal stress regime and in alignment with the microseismicity. Our study reveals how fluid-driven fracture stimulation can be facilitated by a mixed-mode process controlled by the complex hydromechanical evolution of the growing fracture. The data presented in this submission refer to the SIMFIP measurements and analyses of the stimulation tests conducted on the 164 ft (50 m) notch of the Sanford Underground Research Facility (SURF), during the EGS-Collab test 1. In addition to the datafiles, there is the draft of a manuscript submitted to Geophysical Research Letters (GRL).

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.

This is a final technical report for the project: Foam Fracturing Study for Stimulation Development of Enhanced Geothermal Systems (EGS). The goal is to demonstrate the feasibility of foam fracturing in EGS applications. The project, led by Oak Ridge National Laboratory (ORNL), was conducted in collaboration with Temple University. The report describes the research activities with Task 1 at ORNL: foam fracturing testing system development and experimental study on foam fracturing, and Task 2 at Temple University: foam testing and foam characterization. Main findings are: 1. A foam fracturing test system has been developed at ORNL, which can be used to perform foam fracturing under pressure up to 6,000 psi. The system monitors foam density during fracturing online and is capable of testing materials in both monotonic and cyclic (up to 50 Hz) injections. 2. Foam fracturing tests were carried out on Charcoal black granite specimens with a blind borehole to the middle length. Two diameters of blind borehole were tested; G2 series: 9.53 mm and G3 series: 4.76 mm. N2-in-water foam was used with AOS as a surfactant. 3. There was a hole-size effect on fracture initiation pressure. The effect is smaller in the case of foam, which was influenced by the high penetrability of gas in foam. Breakdown pressure showed a behavior just as that of fracture pressure; namely an increased value for small hole samples, while the effect in water fracture was more impressive than in foam fracture. 4. Water mass was reduced in foam fracturing within similar range of breakdown pressures. In G2 series, it was decreased from 10.44 g for water fracturing to 5.17 g, representing more than 50% water reduction. Therefore, there is the potential to reduce water use in EGS stimulation through foam fracturing. 5. Use of cyclic injection has the potential to reduce the breakdown pressure and seismicity in EGS application. Experiments using 4-s cycle period found that specimens can be fractured with a low number of cycles. The fatigue pressure was approximately 64 - 77% of monotonic breakdown pressure for water fracturing and 58 - 94% of the breakdown pressure for foam fracturing. 6. A foam stability testing system has been developed that can test foam at 220 Deg C to 2,000 psi. Tested components of candidate foams included two gases: N2 and CO2; 4 surfactants: AOS, SDS, NP-40 and CTAC; 5 stabilizing agents: guar, bentonite clay, borate salt, silica NPs, and GO. 7. N2 and AOS provided the most stable performance over the tested ranges. Furthermore, the AOS foam with stabilizing agents of guar and borate salt (crosslinker) offered the highest half-life of 20 minutes at 200 Deg C and 1,000 psi. 8. Arrhenius equation and modified power law have been demonstrated to fit well the half-time vs. temperature and pressure data, respectively. These relations can be useful to provide the suggestion for future foam stability study. This submission contains the supporting data developed during the project: 1) A final technical report 2) Granite fracturing data in monotonic and cyclic injections with water and N2 foam Foam performance data in various temperatures and pressures, including half-time, is submitted separately.

Processed SEGY files from the FutureGen2 project. Data includes 5 categories of OVSP data with 15 lines in each category (all in TWT). VSP data was processed from Sterling Seismic.

To prepare for its third phase, the Hawaii Play Fairway project conducted groundwater sampling and analyses in ten locations in the Hawaiian islands, magnetotelluric (MT) and gravity surveys, as well as calculations of 3D subsurface stress due to the weight of the rock underlying the topography of the volcano. The subsurface stresses were used to evaluate the potential for fracture-induced permeability. Inversions of the MT and gravity data produce 3D models of resistivity and density, respectively, on Lanai, across Haleakala's SW rift (Maui), and surrounding Mauna Kea (Hawaii Island). The project developed and applied a new method for incorporating depth information about resistivity, density, and potential for fracture-induced permeability into the statistical method for computing resource probability in these three focus areas. The project then incorporated the new groundwater results with the new geophysical results and the calculations of potential for fracture-induced permeability to produce updated maps of resource probability and confidence. These results were used to identify target sites for exploratory drilling. Spreadsheet information: Each sheet contains data for a particular depth in kilometers. Positive depths are above sea level, and negative below. For more information, go to the Hawaii Groundwater and Geothermal Resources Center website linked in the resources.

Hydraulic Fracturing and Associated Stress Modeling for the Eastern Gas Shale Project - Final Report; November 1980

The Global Groundwater Information System (GGIS) is an interactive, web-based portal to groundwater-related information and knowledge. The GGIS consists of several modules structured around various themes. Each module has its own map-based viewer with underlying database to allow storing and visualizing geospatial data in a systematic way. Data sets include global data on transboundary aquifers, global groundwater data by aquifer, and country disaggregation, global groundwater stress (based on GRACE data), global groundwater quality data. There is also specific regional/national data focusing on the following aquifers: Dinaric Karst (Balkans), Ramotswa and Stampriet aquifers (Southern Africa), Esquipulas-Ocotepeque-Citala (Central Amerca), Pretashkent Aquifer (Central Asia). It also provides access to SADC Groundwater Information Portal, and groundwater on Small Island States

In this paper, we present an analysis using unsupervised machine learning (ML) to identify the key geologic factors that contribute to the geothermal production in Brady geothermal field. Brady is a hydrothermal system in northwestern Nevada that supports both electricity production and direct use of hydrothermal fluids. Transmissive fuid-fow pathways are relatively rare in the subsurface, but are critical components of hydrothermal systems like Brady and many other types of fuid-fow systems in fractured rock. Here, we analyze geologic data with ML methods to unravel the local geologic controls on these pathways. The ML method, non-negative matrix factorization with k-means clustering (NMFk), is applied to a library of 14 3D geologic characteristics hypothesized to control hydrothermal circulation in the Brady geothermal field. Our results indicate that macro-scale faults and a local step-over in the fault system preferentially occur along production wells when compared to injection wells and non-productive wells. We infer that these are the key geologic characteristics that control the through-going hydrothermal transmission pathways at Brady. Our results demonstrate: (1) the specific geologic controls on the Brady hydrothermal system and (2) the efficacy of pairing ML techniques with 3D geologic characterization to enhance the understanding of subsurface processes. This submission includes the published journal article detailing this work, the published 3D geologic map of the Brady Geothermal Area used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity, 3D well data, along which geologic data were sampled for PCA analyses, and associated metadata file. This work was done using the GeoThermalCloud framework, which is part of SmartTensors (both are linked below).

This data set includes the magnetotelluric (MT) data collected from October 21 to November 9, 2016 over the San Emidio geothermal field in Nevada by Quantec Geoscience USA Inc. on behalf of US Geothermal Inc. as part of a project entitled "A Novel Approach to Map Permeability Using Passive Seismic Emission Tomography". This data set includes descriptions of the instrumentation, data acquisition and processing procedures, as well as the final processed data and digital archive formats. A total of 81 MT locations were surveyed (52 profile sites, and 29 MT sites). Data were processed and inspected for quality assurance on site, and reviewed daily by the geophysicist in charge of the project.

The primary objective of this project is to develop a three-blade MHK rotor with low manufacturing and maintenance costs. The proposed program will design, fabricate and test a novel half-scale low cost, net shape fabricated single piece three-blade MHK rotor with integrated health management technology to demonstrate significant Capital Expenditures (CAPEX) and Operational Expenditures (OPEX) cost reductions due to the novel design and manufacturing process. The proposed project is divided into three major tasks: Task 1: Single Piece Three-blade Kinetic Hydropower System (KHPS) Rotor Full-Scale and Half-Scale Design; Task 2: Composite Manufacturing Trials and Half-Scale Prototype Rotor Fabrication; and Task 3: Material Characterization and Half-Scale Prototype Test and Evaluation. These three tasks include design and analysis of full-scale and half-scale three-blade rotor prototypes using computational fluid dynamics (CFD) and finite-element analysis (FEA), demonstration of a novel half-scale net shape fabrication process, determination of a fatigue threshold composite strain allowable, three-blade rotor mold design, manufacture of half-scale rotor clam shell mold, three-blade rotor test fixture design and fabrication, development of final manufacturing and test plans, manufacture of the half-scale net shape composite single blade and three-blade prototypes, and test and evaluation of the half-scale rotor.

The utility of passive seismic emission tomography for mapping geothermal permeability has been tested at San Emidio in Nevada. The San Emidio study area overlaps a geothermal field in production since 1987 and another resource to the south of the production field. Passive seismic data collections were completed at San Emidio in late 2016 by Microseismic Inc as part of a DOE project. The PSET results are being analyzed as part of the WHOLESCALE project. This submission includes P-wave velocity model data, and the passive seismic data with more information on each bellow.

Processed SEGY files from the Perch 3D seismic area in the Michigan Basin. Data was used to make estimates of the stress field in the study area. Data was processed from Texseis, Inc.

Combined geochemical and geophysical data, weighted and ranked for geothermal prospect favorability.

As part of the planning for stimulation of the Newberry Volcano Enhanced Geothermal Systems (EGS) Demonstration project in Oregon, a high-resolution borehole televiewer (BHTV) log was acquired using the ALT ABI85 BHTV tool in the slightly deviated NWG 55-29 well. The image log reveals an extensive network of fractures in a conjugate set striking approximately N-S and dipping 50 deg that are well oriented for normal slip and are consistent with surface-breaking regional normal faults in the vicinity. Similarly, breakouts indicate a consistent minimum horizontal stress, Shmin, azimuth of 092.3 +/- 17.3 deg. In conjunction with a suite of geophysical logs, a model of the stress magnitudes constrained by the width of breakouts at depth and a model of rock strength independently indicates a predominantly normal faulting stress regime.

This entry contains raw data files from experiments performed on the Vulcan beamline at the Spallation Neutron Source at Oak Ridge National Laboratory using a pressure cell. Cylindrical granite and marble samples were subjected to confining pressures of either 0 psi or approximately 2500 psi and internal pressures of either 0 psi, 1500 psi or 2500 psi through a blind axial hole at the center of one end of the sample. The sample diameters were 1.5" and the sample lengths were 6". The blind hole was 0.25" in diameter and 3" deep. One set of experiments measured strains at points located circumferentially around the center of the sample with identical radii to determine if there was strain variability (this would not be expected for a homogeneous material based on the symmetry of loading). Another set of experiments measured load variation across the radius of the sample at a fixed axial and circumferential location. Raw neutron diffraction intensity files and experimental parameter descriptions are included.

EGS field projects have not sustained production at rates greater than 1/2 of what is needed for economic viability. The primary limitation that makes commercial EGS infeasible is our current inability to cost-effectively create high-permeability reservoirs from impermeable, igneous rock within the 3,000-10,000 ft depth range. Our goal is to develop a novel fracturing fluid technology that maximizes reservoir permeability while reducing stimulation cost and environmental impact. Laboratory equipment development to advance laboratory characterization/monitoring is also a priority of this project to study and optimize the physicochemical properties of these fracturing fluids in a range of reservoir conditions. Barrier G is the primarily intended GTO barrier to be addressed as well as support addressing barriers D, E and I.

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.

In December 2016, 1301 vertical-component seismic instruments were deployed at the San Emidio Geothermal field in Nevada. The first record starts at 2016-12-05T02:00:00.000000Z (UTC) and the last record ends at 2016-12-11T14:00:59.998000Z (UTC). Data are stored in individual files in one-minute increments in SEGD and MSEED formats. See the metadata in GDR submission (linked below as "Seismic Survey 2016 Metadata at San Emidio Nevada") for details about the seismic station locations, seismic data logger specifications, instrumentation specifications, descriptions of data, a fracture finding summary, and the final report for the 2016 seismic survey done in San Emidio, Nevada.

1301 Vertical Component seismic instruments were deployed at San Emidio Geothermal field in Nevada in December 2016. The first record starts at 2016-12-05T02:00:00.000000Z (UTC) and the last record ends at 2016-12-11T14:00:59.998000Z (UTC). Data are stored in individual files in one-minute increments. Data includes seismic station locations, seismic data logger specifications, instrumentation specifications, descriptions of data, a fracture finding summary and the final report for the 2016 WHOLESCALE seismic survey done in San Emidio, Nevada.

Slip and Dilation Tendency in focus areas Critically stressed fault segments have a relatively high likelihood of acting as fluid flow conduits (Sibson, 1994). As such, the tendency of a fault segment to slip (slip tendency; Ts; Morris et al., 1996) or to dilate (dilation tendency; Td; Ferrill et al., 1999) provides an indication of which faults or fault segments within a geothermal system are critically stressed and therefore likely to transmit geothermal fluids. The slip tendency of a surface is defined by the ratio of shear stress to normal stress on that surface: Ts = T / on (Morris et al., 1996). Dilation tendency is defined by the stress acting normal to a given surface: Td = (o1-on) / (o1-o3) (Ferrill et al., 1999). Slip and dilation were calculated using 3DStress (Southwest Research Institute). Slip and dilation tendency are both unitless ratios of the resolved stresses applied to the fault plane by ambient stress conditions. Values range from a maximum of 1, a fault plane ideally oriented to slip or dilate under ambient stress conditions to zero, a fault plane with no potential to slip or dilate. Slip and dilation tendency values were calculated for each fault in the focus study areas at, McGinness Hills, Neal Hot Springs, Patua, Salt Wells, San Emidio, and Tuscarora on fault traces. As dip is not well constrained or unknown for many faults mapped in within these we made these calculations using the dip for each fault that would yield the maximum slip tendency or dilation tendency. As such, these results should be viewed as maximum tendency of each fault to slip or dilate. The resulting along-fault and fault-to-fault variation in slip or dilation potential is a proxy for along fault and fault-to-fault variation in fluid flow conduit potential. Stress Magnitudes and directions Stress field variation within each focus area was approximated based on regional published data and the world stress database (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2010; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012; Moeck et al., 2010; Moos and Ronne, 2010 and Reinecker et al., 2005) as well as local stress information if applicable. For faults within these focus systems we applied either a normal faulting stress regime where the vertical stress (sv) is larger than the maximum horizontal stress (shmax) which is larger than the minimum horizontal stress (sv>shmax>shmin) or strike-slip faulting stress regime where the maximum horizontal stress (shmax) is larger than the vertical stress (sv) which is larger than the minimum horizontal stress (shmax >sv>shmin) depending on the general tectonic province of the system. Based on visual inspection of the limited stress magnitude data in the Great Basin we used magnitudes such that shmin/shmax = .527 and shmin/sv= .46, which are consistent with complete and partial stress field determinations from Desert Peak, Coso, the Fallon area and Dixie valley (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2011; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012). Slip and dilation tendency for the McGinness Hills geothermal field was calculated based on the faults mapped McGinness Hills area (Siler 2012, unpublished). The McGinness Hills area lies in the Basin and Range Province, as such we applied a normal faulting stress regime to the McGinness area faults, with a minimum horizontal stress direction oriented 115, based on inspection of local and regional stress determinations, as explained above. Under these stress conditions north-northeast striking, steeply dipping fault segments have the highest dilation tendency, while north-northeast striking 60 degrees dipping fault segments have the highest tendency to slip. The McGinness Hills geothermal system is characterized by a left-step in a north-northeast striking west-dipping fault system within a north northeast striking accommodation zone. As such, the normal faults that define these two structures are well oriented for both slip and dilation, including the west dipping faults that are exploited for both production and injection. Interestingly, although there is pressure communication between production and injection wells at McGinness Hills (B. Delwiche, personal comm.) the northwest striking fault, which creates hard linkage between the production and injection locations, is poorly oriented for both slip and dilation and therefore unlikely to host permeability. NOTE: 'o' is used in this description to represent lowercase sigma.

Slip and Dilation Tendency in focus areas Critically stressed fault segments have a relatively high likelihood of acting as fluid flow conduits (Sibson, 1994). As such, the tendency of a fault segment to slip (slip tendency; Ts; Morris et al., 1996) or to dilate (dilation tendency; Td; Ferrill et al., 1999) provides an indication of which faults or fault segments within a geothermal system are critically stressed and therefore likely to transmit geothermal fluids. The slip tendency of a surface is defined by the ratio of shear stress to normal stress on that surface: Ts = T / on (Morris et al., 1996). Dilation tendency is defined by the stress acting normal to a given surface: Td = (o1-on) / (o1-o3) (Ferrill et al., 1999). Slip and dilation were calculated using 3DStress (Southwest Research Institute). Slip and dilation tendency are both unitless ratios of the resolved stresses applied to the fault plane by ambient stress conditions. Values range from a maximum of 1, a fault plane ideally oriented to slip or dilate under ambient stress conditions to zero, a fault plane with no potential to slip or dilate. Slip and dilation tendency values were calculated for each fault in the focus study areas at, McGinness Hills, Neal Hot Springs, Patua, Salt Wells, San Emidio, and Tuscarora on fault traces. As dip is not well constrained or unknown for many faults mapped in within these we made these calculations using the dip for each fault that would yield the maximum slip tendency or dilation tendency. As such, these results should be viewed as maximum tendency of each fault to slip or dilate. The resulting along-fault and fault-to-fault variation in slip or dilation potential is a proxy for along fault and fault-to-fault variation in fluid flow conduit potential. Stress Magnitudes and directions Stress field variation within each focus area was approximated based on regional published data and the world stress database (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2010; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012; Moeck et al., 2010; Moos and Ronne, 2010 and Reinecker et al., 2005) as well as local stress information if applicable. For faults within these focus systems we applied either a normal faulting stress regime where the vertical stress (sv) is larger than the maximum horizontal stress (shmax) which is larger than the minimum horizontal stress (sv>shmax>shmin) or strike-slip faulting stress regime where the maximum horizontal stress (shmax) is larger than the vertical stress (sv) which is larger than the minimum horizontal stress (shmax >sv>shmin) depending on the general tectonic province of the system. Based on visual inspection of the limited stress magnitude data in the Great Basin we used magnitudes such that shmin/shmax = .527 and shmin/sv= .46, which are consistent with complete and partial stress field determinations from Desert Peak, Coso, the Fallon area and Dixie valley (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2011; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012). Based on inversion of fault kinematic data, Edwards (2013) interpreted that two discrete stress orientations are preserved at Neal Hot Springs. An older episode of east-west directed extension and a younger episode of southwest-northeast directed sinistral, oblique -normal extension. This interpretation is consistent with the evolution of Cenozoic tectonics in the region (Edwards, 2013). As such we applied a southwest-northeast (060) directed normal faulting stress regime, consistent with the younger extensional episode, to the Neal Hot Springs faults. Under these stress conditions northeast striking steeply dipping fault segments have the highest tendency to dilate and northeast striking 60 degrees dipping fault segments have the highest tendency to slip. Under these stress conditions, both the Neal Fault and Sugarloaf Butte faults area well-oriented for both slip and dilation and thus for fracture permeability. In addition, several subsidiary faults on the eastern side and within the step-over between the Neal fault and Sugarloaf Butte fault are well oriented for slip and dilation as well. NOTE: 'o' is used in this description to represent lowercase sigma.

Critically stressed fault segments have a relatively high likelihood of acting as fluid flow conduits (Sibson, 1994). As such, the tendency of a fault segment to slip or to dilate provides an indication of which faults or fault segments within a geothermal system are critically stressed and therefore likely to transmit geothermal fluids. Slip and dilation were calculated using 3DStress (Southwest Research Institute). Slip and dilation tendency are both unitless ratios of the resolved stresses applied to the fault plane by ambient stress conditions. Values range from a maximum of 1, a fault plane ideally oriented to slip or dilate under ambient stress conditions to zero, a fault plane with no potential to slip or dilate. Slip and dilation tendency values were calculated for each fault in the focus study areas at, McGinness Hills, Neal Hot Springs, Patua, Salt Wells, San Emidio, and Tuscarora on fault traces. Stress Magnitudes and directions Stress field variation within each focus area was approximated based on regional published data and the world stress database as well as local stress information. Slip and dilation tendency analysis for the Patua geothermal system was calculated based on faults mapped in the Hazen Quadrangle (Faulds et al., 2011). Patua lies near the margin between the Basin and Range province, which is characterized by west-northwest directed extension and the Walker Lane province, characterized by west-northwest directed dextral shear. As such, the Patua area likely has been affected by tectonic stress associated with either or both of stress regimes over geologic time. In order to characterize this stress variation we calculated slip tendency at Patua for both normal faulting and strike slip faulting stress regimes. Dilation tendency results for a strike-slip faulting stress regime and for a normal faulting stress regime are virtually identical, so we present one result for dilation tendency applicable to both strike-slip and normal faulting stress conditions along with slip tendency for both a normal faulting and a strike-slip faulting stress regime. Under these stress conditions, north-northeast striking steeply dipping fault segments have the highest dilation tendency. Under the strike-slip faulting stress regime, north-northwest and east-northeast striking, steeply dipping fault have the highest slip tendency, while under normal faulting conditions north northeast striking, 60 degrees dipping faults have the highest slip tendency.

Critically stressed fault segments have a relatively high likelihood of acting as fluid flow conduits (Sibson, 1994). As such, the tendency of a fault segment to slip (slip tendency; Ts; Morris et al., 1996) or to dilate (dilation tendency; Td; Ferrill et al., 1999) provides an indication of which faults or fault segments within a geothermal system are critically stressed and therefore likely to transmit geothermal fluids. The slip tendency of a surface is defined by the ratio of shear stress to normal stress on that surface: Ts = T / on (Morris et al., 1996). Dilation tendency is defined by the stress acting normal to a given surface: Td = (o1-on) / (o1-o3) (Ferrill et al., 1999). Slip and dilation were calculated using 3DStress (Southwest Research Institute). Slip and dilation tendency are both unitless ratios of the resolved stresses applied to the fault plane by ambient stress conditions. Values range from a maximum of 1, a fault plane ideally oriented to slip or dilate under ambient stress conditions to zero, a fault plane with no potential to slip or dilate. Slip and dilation tendency values were calculated for each fault in the focus study areas at, McGinness Hills, Neal Hot Springs, Patua, Salt Wells, San Emidio, and Tuscarora on fault traces. As dip is not well constrained or unknown for many faults mapped in within these we made these calculations using the dip for each fault that would yield the maximum slip tendency or dilation tendency. As such, these results should be viewed as maximum tendency of each fault to slip or dilate. The resulting along-fault and fault-to-fault variation in slip or dilation potential is a proxy for along fault and fault-to-fault variation in fluid flow conduit potential. Stress Magnitudes and directions Stress field variation within each focus area was approximated based on regional published data and the world stress database (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2010; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012; Moeck et al., 2010; Moos and Ronne, 2010 and Reinecker et al., 2005) as well as local stress information if applicable. For faults within these focus systems we applied either a normal faulting stress regime where the vertical stress (sv) is larger than the maximum horizontal stress (shmax) which is larger than the minimum horizontal stress (sv>shmax>shmin) or strike-slip faulting stress regime where the maximum horizontal stress (shmax) is larger than the vertical stress (sv) which is larger than the minimum horizontal stress (shmax >sv>shmin) depending on the general tectonic province of the system. Based on visual inspection of the limited stress magnitude data in the Great Basin we used magnitudes such that shmin/shmax = .527 and shmin/sv= .46, which are consistent with complete and partial stress field determinations from Desert Peak, Coso, the Fallon area and Dixie valley (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2011; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012). Slip and dilation tendency for the Salt Wells geothermal field was calculated based on the faults mapped in the Bunejug Mountains quadrangle (Hinz et al., 2011). The Salt Wells area lies in the Basin and Range Province (N. Hinz personal comm.) As such we applied a normal faulting stress regime to the Salt Wells area faults, with a minimum horizontal stress direction oriented 105, based on inspection of local and regional stress determinations. Under these stress conditions north-northeast striking, steeply dipping fault segments have the highest dilation tendency, while north-northeast striking 60 degrees dipping fault segments have the highest tendency to slip. Several such faults intersect in high density in the core of the accommodation zone in the Bunejug Mountains and local to the Salt Wells geothermal .

Critically stressed fault segments have a relatively high likelihood of acting as fluid flow conduits (Sibson, 1994). As such, the tendency of a fault segment to slip (slip tendency; Ts; Morris et al., 1996) or to dilate (dilation tendency; Td; Ferrill et al., 1999) provides an indication of which faults or fault segments within a geothermal system are critically stressed and therefore likely to transmit geothermal fluids. The slip tendency of a surface is defined by the ratio of shear stress to normal stress on that surface: Ts = T / on (Morris et al., 1996). Dilation tendency is defined by the stress acting normal to a given surface: Td = (o1-on) / (o1-o3) (Ferrill et al., 1999). Slip and dilation were calculated using 3DStress (Southwest Research Institute). Slip and dilation tendency are both unitless ratios of the resolved stresses applied to the fault plane by ambient stress conditions. Values range from a maximum of 1, a fault plane ideally oriented to slip or dilate under ambient stress conditions to zero, a fault plane with no potential to slip or dilate. Slip and dilation tendency values were calculated for each fault in the focus study areas at, McGinness Hills, Neal Hot Springs, Patua, Salt Wells, San Emidio, and Tuscarora on fault traces. As dip is not well constrained or unknown for many faults mapped in within these we made these calculations using the dip for each fault that would yield the maximum slip tendency or dilation tendency. As such, these results should be viewed as maximum tendency of each fault to slip or dilate. The resulting along-fault and fault-to-fault variation in slip or dilation potential is a proxy for along fault and fault-to-fault variation in fluid flow conduit potential. Stress Magnitudes and directions Stress field variation within each focus area was approximated based on regional published data and the world stress database (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2010; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012; Moeck et al., 2010; Moos and Ronne, 2010 and Reinecker et al., 2005) as well as local stress information if applicable. For faults within these focus systems we applied either a normal faulting stress regime where the vertical stress (sv) is larger than the maximum horizontal stress (shmax) which is larger than the minimum horizontal stress (sv>shmax>shmin) or strike-slip faulting stress regime where the maximum horizontal stress (shmax) is larger than the vertical stress (sv) which is larger than the minimum horizontal stress (shmax >sv>shmin) depending on the general tectonic province of the system. Based on visual inspection of the limited stress magnitude data in the Great Basin we used magnitudes such that shmin/shmax = .527 and shmin/sv= .46, which are consistent with complete and partial stress field determinations from Desert Peak, Coso, the Fallon area and Dixie valley (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2011; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012). Slip and dilation tendency for the San Emidio geothermal field was calculated based on the faults mapped Tuscarora area (Rhodes, 2011). The San Emidio area lies in the Basin and Range Province, as such we applied a normal faulting stress regime to the San Emidio area faults, with a minimum horizontal stress direction oriented 115, based on inspection of local and regional stress determinations, as explained above. This is consistent with the shmin determined through inversion of fault data by Rhodes (2011). Under these stress conditions north-northeast striking, steeply dipping fault segments have the highest dilation tendency, while north-northeast striking 60 degrees dipping fault segments have the highest tendency to slip. Interesting, the San Emidio geothermal field lies in an area of primarily north striking faults, which have moderate dilation tendency and moderate to low slip tendency. NOTE: 'o' is used in this description to represent lowercase sigma.

Critically stressed fault segments have a relatively high likelihood of acting as fluid flow conduits (Sibson, 1994). As such, the tendency of a fault segment to slip (slip tendency; Ts; Morris et al., 1996) or to dilate (dilation tendency; Td; Ferrill et al., 1999) provides an indication of which faults or fault segments within a geothermal system are critically stressed and therefore likely to transmit geothermal fluids. The slip tendency of a surface is defined by the ratio of shear stress to normal stress on that surface: Ts = T / on (Morris et al., 1996). Dilation tendency is defined by the stress acting normal to a given surface: Td = (o1-on) / (o1-o3) (Ferrill et al., 1999). Slip and dilation were calculated using 3DStress (Southwest Research Institute). Slip and dilation tendency are both unitless ratios of the resolved stresses applied to the fault plane by ambient stress conditions. Values range from a maximum of 1, a fault plane ideally oriented to slip or dilate under ambient stress conditions to zero, a fault plane with no potential to slip or dilate. Slip and dilation tendency values were calculated for each fault in the focus study areas at, McGinness Hills, Neal Hot Springs, Patua, Salt Wells, San Emidio, and Tuscarora on fault traces. As dip is not well constrained or unknown for many faults mapped in within these we made these calculations using the dip for each fault that would yield the maximum slip tendency or dilation tendency. As such, these results should be viewed as maximum tendency of each fault to slip or dilate. The resulting along-fault and fault-to-fault variation in slip or dilation potential is a proxy for along fault and fault-to-fault variation in fluid flow conduit potential. Stress Magnitudes and directions Stress field variation within each focus area was approximated based on regional published data and the world stress database (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2010; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012; Moeck et al., 2010; Moos and Ronne, 2010 and Reinecker et al., 2005) as well as local stress information if applicable. For faults within these focus systems we applied either a normal faulting stress regime where the vertical stress (sv) is larger than the maximum horizontal stress (shmax) which is larger than the minimum horizontal stress (sv>shmax>shmin) or strike-slip faulting stress regime where the maximum horizontal stress (shmax) is larger than the vertical stress (sv) which is larger than the minimum horizontal stress (shmax >sv>shmin) depending on the general tectonic province of the system. Based on visual inspection of the limited stress magnitude data in the Great Basin we used magnitudes such that shmin/shmax = .527 and shmin/sv= .46, which are consistent with complete and partial stress field determinations from Desert Peak, Coso, the Fallon area and Dixie valley (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2011; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012). Slip and dilation tendency for the Tuscarora geothermal field was calculated based on the faults mapped Tuscarora area (Dering, 2013). The Tuscarora area lies in the Basin and Range Province, as such we applied a normal faulting stress regime to the Tuscarora area faults, with a minimum horizontal stress direction oriented 115, based on inspection of local and regional stress determinations, as explained above. Under these stress conditions north-northeast striking, steeply dipping fault segments have the highest dilation tendency, while north-northeast striking 60 degrees dipping fault segments have the highest tendency to slip. Tuscarora is defined by a left-step in a major north- to-north northeast striking, west-dipping range-bounding normal fault system. Faults within the broad step define an anticlinal accommodation zone wherein east-dipping faults mainly occupy western half of the accommodation zone and west-dipping faults lie in the eastern half of the accommodation zone. The geothermal system resides in the axial part of the accommodation, straddling the two fault dip domains. Within the axial part of the accommodation zone several west-dipping, north northeast-striking faults are well oriented for both slip and dilation, including fault strands that are exploited for both production and injection for the Tuscarora geothermal power plant. NOTE: 'o' is used in this description to represent lowercase sigma.

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.

The state of stress influences several potential risks associated with geologic carbon storage (GCS). Among these are the risk of induced seismicity and leakage due to unintentional hydraulic fracturing. In situ stress can be measured or inferred using several techniques, but these measurements are only made at discrete locations, and they generally provide an incomplete description of the full stress tensor. Because of this there is generally considerable uncertainty in the state of stress within a GCS reservoir and surrounding formations. The state of stress analysis tool (SOSAT) provides a method to use a Bayesian approach to calculate a probability distribution for in situ stress at a particular point using a variety of common data types. Using this stress state probability distribution, the SOSAT then calculates the probability of activating a critically-oriented fault at a specified range of pore pressures. The result is an easy to use physics-based tool for estimating the risk of induced seismicity and unintentional hydraulic fracturing.

This report describes research accomplishments achieved with funding provided through Department of Energy (DOE) Contract DE-FE0031686. The purpose of this DOE funding is to develop methods that can provide key information about stress fields that act on deep rocks without invading the earth to acquire that information. The research objective was to demonstrate methods that extract the azimuth directions of SHmax and SHmin stress in deep rocks from traditional seismic reflection data like the data that are used to explore for deep oil and gas reservoirs. Battelle performed two seismic investigations and achieved estimates of SHmax and SHmin azimuths in both efforts that agreed with local, non-seismic, ground-truth measurements of SHmax orientation. The first procedure utilized vertical seismic profiling (VSP) data; the second effort utilized three-dimensional (3D) seismic data.

This report presents a summary of Task 3 (Laboratory Experimental Characterization of Stress Dependent Wavespeed in Rocks and Application for In-Situ Stress Estimation) of the project titled “A Non-Invasive Approach for Elucidating the Spatial Distribution of In-Situ Stress in Deep Subsurface Geologic Formations Considered for CO2 Storage”. The main objective of this component of the project was to employ non-destructive techniques to determine how compressional (P)-wave and shear (S)-wave propagation velocities are related to triaxial-stress conditions in laboratory test samples of targeted reservoir rocks. The overarching principle is that such laboratory measurements can provide a basis whereby field measurements of wave-propagation velocities in rocks can be used to ascertain the stress state that would be expected to generate the observed wave-propagation behavior.

"This document presents a summary of TASK 4 (Field Testing) of the project entitled “A Non-Invasive Approach for Elucidating the Spatial Distribution of in-situ Stress in Deep Subsurface Geologic Formations Considered for CO2 Storage (FE0031686)”. This three-year project is part of the U.S. Department of Energy (DOE) National Energy Technology Laboratory (NETL) Carbon Storage program to research and address gaps that affect the economics of commercial CCS projects. One of these key gaps is the lack of certainty in predicting the geomechanical impacts of pressure migration due to CO2 injection into a storage complex. This report documents field testing that was conducted under TASK 4 to obtain input data needed to develop the site-specific geomechanical model in TASK 5 and to verify the model results (i.e. calculated stresses) for this site. The field work entailed collecting geophysical logs and core samples and conducting geomechanical stress tests in the Core Energy LLC State Otsego Lake (SOL 8-15A) well, located in Otsego County, Michigan."

This report describes research accomplishments achieved with funding provided through Department of Energy (DOE) Contract DE-FE0031686, for the project “A Non-Invasive Approach for Elucidating the Spatial Distribution of In Situ Stress in Deep Subsurface Geologic Formations Considered for CO2 Storage.” The purpose of this DOE funding is to develop non-invasive methods to obtain important information about subsurface stresses. The overall research goal of the project was to develop and improve methods for determining the spatial distribution of stresses, including magnitude and orientation of the three principal stress components in the subsurface, based on new methods that extract stress information from seismic data combined with well measurements, extended with well tests and logs, and unified in a numerical model that permits computation of the full stress tensor throughout the domains under investigation.

Laboratory Experiments Data of TASK3 of DOE project DE- FE0031686. Each signal has its own excel file. 150 signals are sent through the rock sample per each load configuration. There were approximately 50 to 150 load combinations tested during the experiments. Experiments were repeated after changing the polarity direction of the transducers. Signals of these repeats are in separate folders called "Reverse". There are three formations tested from the Michigan Core and four formations tested from the FutureGen site. The "Code" folder contains the Python code developed to analyze the signals. Large Dataset, please contact EDXSupport@netl.doe.gov

Project data and Matlab scripts for Task 5 of the project “A Non-Invasive Approach for Elucidating the Spatial Distribution of In Situ Stress in Deep Subsurface Geologic Formations Considered for CO2 Storage”. Project includes the static earth models, grids and other data used for simulation and calibration of the two case studies addressed in the project, referred to herein as the “FutureGen2” model (Illinois) and the “Perch” model (Otsego county, Michigan). Large Dataset, please contact EDXSupport@netl.doe.gov

Mesh, properties, initial conditions, injection/withdrawal rates for modelling thermal, hydrological, and mechanical effects of fluid injection to and withdrawal from ground for Stockton University reservoir cooling system (aquifer storage cooling system), Galloway, New Jersey, for unscheduled two hour injection at 133 % designed capacity, on fine scale grid, with some results. Second simulation of J.T. Smith, E. Sonnenthal, P. Dobson, P. Nico, and M. Worthington, 2021. Thermal-hydrological-mechanical modeling of Stockton University reservoir cooling system, Proceedings of the 46th Workshop on Geothermal Reservoir Engineering, Stanford University, SGP-TR-218, from which Figures 6-9, pertain.

This is the Phase 3 native state model update. The Phase 3 numerical model represents a significant subsurface volume below the FORGE site footprint. The model domain of 4.0 km x 4.0 km x 4.2 km is located approximately between depths of 4000 to 4200 meters below land surface. This data archive consists of 10 files, 4 of which are simulation input files and the remaining 6 are simulation output files. There is an included readme.txt file that contains details on each of the data files. The input files include meshes, FALCON code inputs, tabulated data of water properties, temperature values, and model boundaries. The output files include simulation outfiles and point data of modeled material properties.

Six samples were evaluated in unconfined and triaxial compression, their data are included in separate excel spreadsheets, and summarized in the word document. Three samples were plugged along the axis of the core (presumed to be nominally vertical) and three samples were plugged perpendicular to the axis of the core. A designation of "V"indicates vertical or the long axis of the plugged sample is aligned with the axis of the core. Similarly, "H" indicates a sample that is nominally horizontal and cut orthogonal to the axis of the core. Stress-strain curves were made before and after the testing, and are included in the word doc. The confining pressure for this test was 2800 psi. A series of tests are being carried out on to define a failure envelope, to provide representative hydraulic fracture design parameters and for future geomechanical assessments. The samples are from well 52-21, which reaches a maximum depth of 3581 ft +/- 2 ft into a gneiss complex.

This submission contains information on thirty-six rock samples collected from San Emidio, Nevada during January, 2021 for Subtask 2.3 of the WHOLESCALE project. The following resources include a .zip of rock sample photos taken in the field, a .zip of rock sample photos taken in the laboratory at UW-Madison, and an excel catalog of rock samples with information on sample name, rock type, coordinates of sample location, structural measurements, field notes, observations for plug preparation (e.g., weathering, ability to be cut and cored), and rock descriptions. It should be noted that not every sample was photographed in the field. Names and descriptions of rock formation units are taken from Rhodes et al. (2011). The README.txt file is a description of this submission.

Routine and geomechanical well logs for 3 wells— State Otsego Lake 8-15A, Lawnichak 9-33, Chester 8-16.

These core photos and descriptions were taken from the five boreholes that were drilled as part of the kISMET SubTER project conducted at the Sanford Underground Research Facility (SURF) in Lead, SD. The boreholes are subvertical in orientation, and were drilled on the 4850 level of SURF on the West Drift, about 450 feet from Governor's Corner. The well heads for the five wells are in a line, but the outer two wells (k001 and k005) were deviated to form a five-spot configuration at 50 m depth. Four of the five boreholes have a nominal depth of 50 m and have HQ core - the fifth, located in the center (k003) was drilled to a depth of 100m and has NQ core. The central borehole was used for stress and hydraulic fracturing - the other four boreholes were used for monitoring purposes. Core logging was conducted by Paul Cook (LBNL), Bill Roggenthen (SDSMT), and Drew Siler (LBNL). All core consists of rocks from the Poorman Formation. Some of the core photos are missing. These have been documented in the included spreadsheets labeled with the well name and the word missing. The locations of the boreholes are documented on the included map and spreadsheet.

The files in this submission describes the results of a series of stress measurement and hydraulic fracturing experiments conducted at the Sanford Underground Research Facility (SURF) in Lead, SD. This report describes the accomplishments of the kISMET (permeability (k) and Induced Seismicity Management for Energy Technologies) project. Five near-vertical boreholes were drilled and cored on the 4850 level of SURF in phyllite of the Precambrian Poorman Formation, and a series of hydraulic fracture experiments and stress measurements were conducted in the central borehole: the outer boreholes were used for monitoring purposes.