DOI: 10.5880/gfz.4.2.2018.002
title: Supplementary material to "Sensitivity and stability analysis of coda quality factors at The Geysers geothermal field, California"

description: This data set is supplementary to the BSSA research article of Blanke et al. (2019), in which the local S-wave coda quality factor at The Geysers geothermal field, California, is investigated. Over 700 induced microseismic events recorded between June 2009 and March 2015 at 31 short-period stations of the Berkeley-Geysers Seismic Network were used to estimate the frequency-dependent coda quality factor (Q_C) using the method of Phillips (1985). A sensitivity analysis was performed to different input parameters (magnitude range, lapse time, moving window width, total coda length and seismic sensor component) to gain a better overview on how these parameters influence Q_C estimates. Tested parameters mainly show a low impact on the outcome whereas applied quality criteria like signal-to-noise ratio and allowed uncertainties of Q_C estimates were found to be the most sensitive factors.
      Frequency-dependent mean-Q_C curves were calculated from seismograms of induced earthquakes for each station located at The Geysers using the tested favored input parameters. The final results were tested in the context of spatio-temporal behavior of Q_C in the reservoir considering distance-, azimuth and geothermal production rate variations. A distance and azimuthal dependence was found which is related to the reservoir anisotropy, lithological-, and structural features. By contrast, variations in geothermal production rates do not influence the estimates. In addition, the final results were compared with previous estimated frequency-independent intrinsic direct S-wave quality factors (Q_D) of Kwiatek et al. (2015). A match of Q_D was observed with Q_C estimates obtained at 7 Hz center-frequency, suggesting that Q_D might not be of an intrinsic but of scattering origin at The Geysers. Additionally, Q_C estimates feature lower spreading of values and thus a higher stability.
      The Geysers geothermal field is located approximately 110 km northwest of San Francisco, California in the Mayacamas Mountains. It is the largest steam-dominated geothermal reservoir operating since the 1960s. The local seismicity is clearly related to the water injections and steam production with magnitudes up to ~5 occurring down to 5 km depth, reaching the high temperature zone (up to 360°C). The whole study area is underlain by a felsite (granitic intrusion) that shows an elevation towards the southeast and subsides towards northwest. A fracture network induces anisotropy into the otherwise isotropic rocks featuring different orientations. Moreover, shear-wave splitting and high attenuating seismic signals are observed and motivate to analyze the frequency-dependent coda quality factor.
      Two data sets were analyzed: one distinct cluster located in the northwest (NW) close to injection wells Prati-9 and Prati-29, and the other one southeast (SE) of The Geysers, California, USA, close to station TCH (38° 50′ 08.2″ N, 122° 49′ 33.7″ W and 38° 46′ 59.5″ N, 122° 44′ 13.2″ W, respectively).
      The frequency-dependent coda quality factor is estimated from the seismic S-wave coda by applying the moving window method and regression analysis of Phillips (1985). Different input parameters including moving widow width, lapse time and total coda length are used to obtain Q_C estimates and associated uncertainties. Within a sensitivity analysis we investigated the influence of these parameters and also of magnitude ranges and seismic sensor components on Q_C estimates. The coda analysis was performed for each event at each sensor component of each station. The seismograms were filtered in predefined octave-width frequency bands with center-frequencies ranging from 1-69 Hz. The moving window method was applied starting in the early coda (after the S-onset) for each frequency band measuring the decay of Power Spectral Density spectra. The decay of coda amplitudes was fitted with a regression line and Q_C estimates were calculated from its decay slope for each frequency band. In a final step a mean-Q_C curve was calculated for each available station within the study area resulting in different curves dependent on event location sites in the northwest and southeast.
      Data Description
      The data contain final mean-Q_C estimates of the NW and SE Geysers, coda Q estimates at 7 Hz center-frequency calculated by using the NW cluster, and initial direct Q estimates of Kwiatek et al. (2015) using the same data of the NW cluster. Table S1 shows final mean coda quality factor estimates obtained from the NW cluster at injection wells Prati-9 and Prati-29. The column headers show stations (station), center-frequencies of octave-width frequency bands in Hertz (f[Hz]), mean coda Q estimates (meanQc) and related standard deviations (std), all obtained by coda analysis. Table S2 shows the final mean coda quality factor estimates obtained from additional selected 100 events in the SE Geysers. Column headers correspond to those in Table S1. Table S3 shows coda Q estimates related to 7 Hz center-frequency. The column headers show stations (station), center-frequency of octave-width frequency bands in Hertz (f[Hz]), coda Q estimates at 7 Hz center-frequency (Q_C) and related standard deviations (std2sigma; 95% confidence level), all obtained by coda analysis. Table S4 shows selected direct S-wave quality factors of Kwiatek et al. (2015) obtained by spectral fitting. The column headers show stations (station) and direct S-wave Q estimates (Q_D). The four tables are provided in tab separated txt format.
      Tables S3 and S4 are used for a comparative study and displayed in Figure 12 of the BSSA article mentioned above.