This submission includes synthetic seismic modeling data for the Push-Pull project at Brady Hot Springs, NV. The synthetic seismic is all generated by finite-difference method regarding different fracture and rock properties.

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.

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.

This submission contains datasets and a final manuscript associated with a project simulating carbon dioxide push-pull into a conjugate fault system modeled after Dixie Valley- sensitivity analysis of significant parameters and uncertainty prediction by data-worth analysis. Datasets include: (1) Forward simulation runs of standard cases (push & pull phases), (2) Local sensitivity analyses (push & pull phases), and (3) Data-worth analysis (push & pull phases).

This submission contains datasets and a final manuscript associated with a project simulating carbon dioxide push-pull into a conjugate fault system modeled after Dixie Valley- sensitivity analysis of significant parameters and uncertainty prediction by data-worth analysis. Datasets include: (1) Forward simulation runs of standard cases (push & pull phases), (2) Local sensitivity analyses (push & pull phases), and (3) Data-worth analysis (push & pull phases).

ASCII text files containing grid-block name, X-Y-Z location, and multiple parameters from TOUGH2 simulation output of CO2 injection into an idealized single fault representing a dipping normal fault at the Desert Peak geothermal field (readable by GMS). The fault is composed of a damage zone, a fault gouge and a slip plane. The runs are described in detail in the following: Borgia A., Oldenburg C.M., Zhang R., Jung Y., Lee K.J., Doughty C., Daley T.M., Chugunov N., Altundas B, Ramakrishnan T.S., 2017. Carbon Dioxide Injection for Enhanced Characterization of Faults and Fractures in Geothermal Systems. Proceedings of the 42st Workshop on Geothermal Reservoir Engineering, Stanford University, Stanford, California, February 13-17.

ASCII text files containing grid-block name, X-Y-Z location, and multiple parameters from TOUGH2 simulation output of CO2 injection into an idealized single fault representing a dipping normal fault at the Desert Peak geothermal field (readable by GMS). The fault is composed of a damage zone, a fault gouge and a slip plane. The runs are described in detail in the following: Borgia A., Oldenburg C.M., Zhang R., Jung Y., Lee K.J., Doughty C., Daley T.M., Chugunov N., Altundas B, Ramakrishnan T.S., 2017. Carbon Dioxide Injection for Enhanced Characterization of Faults and Fractures in Geothermal Systems. Proceedings of the 42st Workshop on Geothermal Reservoir Engineering, Stanford University, Stanford, California, February 13-17.

We deployed a dense seismic array to image the shallow structure in the injection area of the Brady Hot Springs geothermal site in Nevada. The array was composed of 238 5 Hz, three-component nodal instruments and 8,700 m of distributed acoustic sensing (DAS) fiber-optic cable installed in surface trenches plus 400 m installed in a borehole. The geophone array had about 60 m instrument spacing in the target zone, whereas DAS channel separations were about 1 m. The acquisition systems provided 15 days of continuous records including active source and ambient noise signals. A large vibroseis truck (T-Rex) was operated at 196 locations exciting a swept-frequency signal from 5 to 80 Hz over 20 seconds using three vibration modes. Sweeps were repeated up to four times during different modes of geothermal plant operation: normal operation, shutdown, high and oscillatory injection and production, and normal operation again. The cross- correlation method was utilized to remove the sweep signal from the geophone records. The first P arrivals were automatically picked from the cross-correlation functions using a combination of methods, and the travel times were used to invert for the 3D P-wave velocity structure. Models with 50 m horizontal node spacing were obtained, with vertical node spacing of 10 to 50 m. The travel time data were fit to about 30 ms, close to our estimated picking uncertainty. Boundaries between high and low velocity zones agree with previous surveys of local faults and low velocity zones near the surface correspond to fumarole locations. A rapid increase in velocity at about 50 m depth fits with borehole data on the depth of the Quaternary sediments. There is some evidence for changes in the P-wave velocity during the experiment with slower travel times at the beginning of the experiment.

We deployed a dense seismic array to image the shallow structure in the injection area of the Brady Hot Springs geothermal site in Nevada. The array was composed of 238 5 Hz, three-component nodal instruments and 8,700 m of distributed acoustic sensing (DAS) fiber-optic cable installed in surface trenches plus 400 m installed in a borehole. The geophone array had about 60 m instrument spacing in the target zone, whereas DAS channel separations were about 1 m. The acquisition systems provided 15 days of continuous records including active source and ambient noise signals. A large vibroseis truck (T-Rex) was operated at 196 locations exciting a swept-frequency signal from 5 to 80 Hz over 20 seconds using three vibration modes. Sweeps were repeated up to four times during different modes of geothermal plant operation: normal operation, shutdown, high and oscillatory injection and production, and normal operation again. The cross- correlation method was utilized to remove the sweep signal from the geophone records. The first P arrivals were automatically picked from the cross-correlation functions using a combination of methods, and the travel times were used to invert for the 3D P-wave velocity structure. Models with 50 m horizontal node spacing were obtained, with vertical node spacing of 10 to 50 m. The travel time data were fit to about 30 ms, close to our estimated picking uncertainty. Boundaries between high and low velocity zones agree with previous surveys of local faults and low velocity zones near the surface correspond to fumarole locations. A rapid increase in velocity at about 50 m depth fits with borehole data on the depth of the Quaternary sediments. There is some evidence for changes in the P-wave velocity during the experiment with slower travel times at the beginning of the experiment.

Included are experimental data recorded from shear experiments that specifically explore the link between fluid-injection rate and seismic moment resulting from shear reactivation of laboratory faults. Raw mechanical data from three experiments are included alongside corresponding MATLAB scripts that import and plot the data, as well as use it to calculate shear and normal stress. Experiments are performed on 2.5-3 inch long granitoid cores from the Utah FORGE EGS demonstration site, containing a single inclined fracture with small-scale roughness added to the fracture surface. The raw data included here were recorded from an aluminum triaxial pressure vessel (TEMCO) configured with three independent servo-controlled pumps, with distilled water used as the working fluid. The pumps control confining pressure, upstream pore pressure, and axial pressure, with each pump connected to a LabView interface to record applied pressures, cumulative injected water volumes, and pump flow rates. The downstream outlet from the fracture is closed to allow pressurization, which is measured by an external pressure transducer. A linear variable differential transformer (LVDT) attached to the axial piston measures axial displacement, from which we calculate shear displacement along the fracture. Additionally, P-wave transducers are used to record acoustic signatures, where acoustic emission events and maximum amplitudes are compared against seismic moment and shear slip velocity. Fluid injection rates range between 0.05 mL/min, 0.25 mL/min, and 0.75 mL/min for each experiment. Along-fault pressure distributions are progressively less uniform as injection rates increase, representing a switch from steady-state to transient conditions. Triggered shear displacement is used as a proxy for seismic moment and is indexed against cumulative injection volume and rate. Each experiment is performed under constant shear stress conditions, and the sample is fully saturated with DI water. Axial and confining stresses are applied to 3 MPa through pressure-stepping in 500 kPa increments. The pore pressure is held constant at 200 kPa prior to initiating the experiment, and initial axial displacement is recorded. The axial stress is then increased to initiate shear mobilization during the loading phase (run-in) until a peak steady state is achieved. The initial shear stress is reduced to approximately 80% of the peak shear stress by decreasing the axial stress, then held constant for the duration of each experiment.