Characterization of transient noise in Advanced LIGO relevant to gravitational wave signal GW150914

A gravitational wave signal, denoted GW150914, has been detected by the Advanced LIGO detectors [1]. The recovered waveform indicated the source was a binary black hole system with component masses ${36}_{-4}^{+5}$ ${M}_{\odot }$ and ${29}_{-4}^{+4}$ ${M}_{\odot }$, which coalesced at a distance of ${410}_{-180}^{+160}$ Mpc away from Earth. The significance of the GW150914 event was measured to be greater than $5.1$ σ, corresponding to a false-alarm rate of less than 1 event per $203\;000$ years [1]. The event, lasting $0.2$ s in Advanced LIGO's sensitive frequency range, was detected in independent searches for modeled compact binary coalescences (CBCs) and for unmodeled gravitational wave bursts [2, 3].

The US-based detectors, in Hanford, Washington (H1) and in Livingston, Louisiana (L1) jointly comprise the Laser Interferometer Gravitational-wave Observatory (LIGO). The detectors are designed to measure spacetime strain induced by passing gravitational waves using a modified Michelson interferometer with 4 km length arms, as described in [4–6]. The detectors were operating in their nominal configuration at the time of GW150914. The corresponding detector sensitivity is shown in figure 1; both detectors achieved a best sensitivity of ∼10⁻²³ Hz${}^{-1/2}$ between roughly 50 and 300 Hz. Peaks in the strain-equivalent noise amplitude spectral density are due largely to mechanical resonances, mains power harmonics, and injected signals used for calibration. Non-stationarity in the detector noise manifests as variations in the level and shape of these sensitivity curves over time.

Zoom In Zoom Out Reset image size

Even in their nominal state, the detectors' data contain non-Gaussian noise transients introduced by behavior of the instruments or complex interactions between the instruments and their environment. For LIGO, the fundamental signature of a transient gravitational wave signal is a near-simultaneous signal with consistent waveforms in the two detectors. The rate of coincident noise transients between the independent detector data sets is estimated by the astrophysical searches using time-shift techniques [2, 3]. A common time-shift method is to shift the data of one detector relative to the other detector's data by a time interval significantly greater than 10 ms, the maximum difference in signal arrival time between detectors. Coincident triggers in time-shifted data yield a distribution of background triggers produced solely by the chance coincidence of transient noise. This time-shifting of the data is performed many times to obtain a representative estimate of the expected rate of background triggers, as detailed in [2, 7]. The significance of a gravitational wave event is a measure of the probability that it is a false detection due to coincident noise. We study the characteristics of background triggers as well as correlations between the gravitational wave strain data and instrument or environment signals to guide further detector improvements and increase the sensitivity of the searches.

GW150914 occurred on 14 September 2015 09:50:45 UTC, 28 days into the eighth engineering run (ER8)¹³⁸ , 3 days into stable data collection with an accurate calibration, and 4 days preceding the scheduled start of the first observing run (O1).

After the event was identified as a highly significant candidate, the software and hardware configuration of each LIGO detector was held fixed until enough coincident data had been collected to set a sufficiently accurate upper bound on the false-alarm rate using the time-shift technique described above. It took roughly six weeks to collect the required ∼16 d of coincident data because low noise operation of the detectors is disrupted by noisy environmental conditions (such as storms, earthquakes, high ground motion, or anthropogenic noise sources). During this six week period we only performed non-invasive maintenance that was required for instrument stability.

The significance of GW150914 was calculated using data taken from 12 September 2015 00:00 through 20 October 2015 13:30 UTC. This data set was analyzed after removing time segments during which an identified instrumental or environmental noise source coupled to the gravitational wave strain signal. At these times, any triggered output of the astrophysical searches would likely be due to noise. These data quality vetoes were built on detector characterization efforts in earlier stages of testing and commissioning of the Advanced LIGO detectors, as reported in [8].

This paper summarizes detector characterization techniques for identification of transient noise (section 2). We then present examples of transient noise couplings that can impact the detectors (section 3) and discuss techniques used to mitigate the impact of known noise sources (section 4). We show that the selected analysis period provides an accurate estimate of the significance of GW150914 reported in [1] by discussing the stability of the search backgrounds, and presenting the impact of applied data quality vetoes relevant to GW150914 (section 5). We also detail the specific checks performed to rule out an instrumental or environmental noise-transient origin for GW150914, including potentially correlated noise sources such as global magnetic noise that would not be captured by time-shift background estimation techniques (section 6). Similar studies were also performed for the second most significant event in the CBC search over the analysis period, designated LVT151012¹³⁹ , observed with a false alarm probability of ∼2% [1, 2, 9].

In addition to the gravitational wave strain data, h(t), each of the LIGO detectors also records over 200 000 auxiliary channels that monitor instrument behavior and environmental conditions. These channels witness a broad spectrum of potential coupling mechanisms, useful for diagnosing instrument faults and identifying noise correlations. Examples of instrument witness channels include measured angular drift of optics, light transmitted through a mirror as detected by a set of photodiodes, and actuation signals used to control optic position in order to maintain optical cavity resonance. In addition to candidate gravitational wave events, we study background triggers for correlation with trends or coincident transient noise in auxiliary channels on the broad scale of hours to days. We also identify correlations on the order of the duration of transient astrophysical signals; a fraction of a millisecond to a few seconds. Systematic correlations are used to generate data quality vetoes used by the astrophysical searches to reduce the background, as described in appendix A.

An important set of auxiliary channels are the physical environment monitor (PEM) sensors, which monitor the local surroundings for potential disturbances that may affect the gravitational wave strain data, such as motion of the ground or optics tables, magnetic field variations, acoustic disturbances, or potentially, cosmic ray showers [10]. A PEM sensor array is distributed throughout each detector site such that external environmental disturbances that could influence the detectors are witnessed with a significantly higher signal-to-noise ratio (SNR) in the PEM sensors than in h(t). The PEM sensors are detailed in appendix B.

The relationship between environmental noise as witnessed by the PEM sensor array and the gravitational wave strain signal h(t) is investigated using injection studies, where an intentional stimulus is introduced and the responses of both PEM sensors and the instrument are analyzed. These injections ensure that the environmental sensors are more sensitive to environmental disturbances than the detector is, and also quantify the coupling between the environment and h(t). Figure 2 illustrates a magnetic field injection test at the LIGO-Hanford detector that measured magnetic field coupling to h(t) as well as the response of the local magnetometer to the injected field. The frequency-dependent coupling between the local magnetic field and h(t) can be calculated from these measurements and used to accurately predict the response of h(t) to the presence of a magnetic field, as witnessed by the local magnetometers. Figure 2 shows an injection performed at one of the strongest coupling locations, in the building containing the beam splitter and most interferometer optics. Other magnetic field injection measurements identical to this test were also conducted for other locations throughout the detector site. Similar injection studies were also conducted for radio, acoustic, and mechanical vibration sources.

Zoom In Zoom Out Reset image size

Transient noise in h(t) must occur within the frequency range targeted by the transient astrophysical searches to affect the background. This range is dictated by the equivalent strain noise of the detectors, as shown in figure 1 for the Hanford and Livingston detectors during the analysis period.

Motivated by this sensitivity curve, the transient astrophysical searches generally limit the search frequency range to above 30 Hz and below 2–3 kHz, or roughly the human-audible range. For example, a binary black hole signal like GW150914 is expected to have power measurable by the Advanced LIGO detectors between roughly 35 and 250 Hz and sources of short-duration noise with similar frequency content could impact the background estimation of such events.

3.1. Uncorrelated noise

The following are examples of uncorrelated local noise features anticipated to be of particular interest or known to have a significant impact on the gravitational wave search backgrounds. The contribution of any uncorrelated noise sources is well estimated using time shifts.

• Some anthropogenic noise sources are likely to produce short duration transients in h(t), such as human activity within one of the rooms that houses the vacuum chambers or infrequent strong ground motion or noise from other nearby locations. To reduce such vibrational or acoustic noise, detector staff do not enter the rooms containing the optical components of the detectors when the detectors are taking data. Any anthropogenic noise that could influence the detector is monitored by an array of accelerometers, seismometers, and microphones.
• Earthquakes can produce ground motion at the detectors with frequencies from approximately 0.03 to 0.1 Hz or higher if the epicenter is nearby [10]. R-waves, the highest amplitude component of seismic waves from an earthquake [11], are the most likely to adversely impact data quality by rendering the detectors inoperable or inducing low frequency optic motion that up-converts to higher frequencies in h(t) via mechanisms such as bilinear coupling of angular motion or light scattering [12]. A network of seismometers installed at the LIGO detectors can easily identify earthquake disturbances.
• Radio frequency (RF) modulation sidebands are used to sense and control a variety of optical cavities within the detector. Two modulations are applied to the input laser field at 9 and 45 MHz [6]. Since the beginning of the analysis period, sporadic periods of a high rate of loud noise transients have been observed at LIGO-Hanford due to a fault in the 45 MHz electro-optic modulator driver system, which then couples to the gravitational wave channel between 10 and 2000 Hz, covering the entire frequency range analyzed by the CBC searches. Data associated with this electronic fault were vetoed and not analyzed. The engineering of this veto, as applied to the GW150914 analysis period, is detailed in appendix A.
• Blip transients are short noise transients that appear in the gravitational wave strain channel h(t) as a symmetric 'teardrop' shape in time–frequency space, typically between 30 and 250 Hz, with the majority of the power appearing at the lowest frequencies, as seen in figure 3. They appear in both detectors independently with modest amplitude. The single detector burst identification algorithm Omicron, which identifies excess power transients using a generic sine-Gaussian time–frequency projection [13, 14], will resolve such noise transients with an SNR ratio of 10–100. No clear correlation to any auxiliary channel has yet been identified. As a result, there is currently no veto available to remove these noise transients from the astrophysical searches. Blip transients contribute to some of the most significant background triggers in both the unmodeled burst and modeled CBC searches. The noise transient shown in figure 3 is one example¹⁴⁰ .

Zoom In Zoom Out Reset image size

The impact of noise sources on the astrophysical searches is discussed in section 5.2.

3.2. Correlated noise

Ideally, when a noise source is identified, the instrument hardware or software is modified to reduce the coupling of the noise to h(t) such that it no longer impacts astrophysical searches. If mitigating the noise source is not viable, as in the case of data collected prior to an instrumental improvement, periods of time in which there are significant problems with the quality of the data are omitted, or vetoed, from transient gravitational wave searches through a procedure similar to those utilized in previous LIGO analyses [21].

There are two different types of data quality products that can be applied as vetoes. Data quality flags typically exclude periods of data on the order of seconds to hours when some reproducible criterion associated with known noise couplings is met [21–24]. For example, a data quality flag might be defined for periods when any of the photodiodes used to sense the laser field in the detector were overflowing their analog-to-digital converters. Data quality triggers are short duration vetoes generated by algorithms that identify significant statistical correlations between a transient in h(t) and transient noise in auxiliary channels [25–28].

Data quality products are applied as vetoes in different categories that depend on the severity of the problem or the impact of individual data quality products on a search's background. Data quality flags used in category 1 collectively indicate times when data should not be analyzed due to a critical issue with a key detector component not operating in its nominal configuration. Since category-1-flagged times indicate major known problems with an instrument they are identically defined across all transient searches. Data quality flags used in category 2 collectively indicate times when a noise source with known physical coupling to h(t) is active. Category 2 vetoes are typically applied after the initial processing of data for a specific search. This approach renders more data useable by the searches because they require unbroken strides of continuous data of up to 620 s for the coherent burst search and up to 2064 s for the CBC searches. There are three considerations for applying a data quality product as a category 2 veto to an astrophysical search: the physical noise coupling mechanism must be understood, the associated veto must have a demonstrated advantageous effect on the background of that search, and the veto must be safe.

The safety of a veto is a measure of the likelihood that the veto criteria would accidentally remove a true gravitational wave signal. Veto safety is measured using hardware injection tests, where a signal is injected into h(t) by inducing motion of the optics [25, 26, 29]. If any auxiliary channels witness a corresponding response to a number of injected signals greater than expected by chance, these channels are considered unsafe and are not used in the definition of any applied veto.

The effectiveness of each data quality product in reducing the background is measured by the ratio of its efficiency, or the fraction of background triggers it removes from a search, to its introduced deadtime, or the fraction of time a particular flag will remove from the total duration of the set of analyzable data. Data quality flags used as category 2 vetoes have an efficiency-to-deadtime ratio for high SNR triggers significantly greater than 1, or the value expected for random behavior. An example is described in appendix A.

A third veto category (category 3), applied in the same way as category 2, is generally reserved for data quality triggers, which are statistically generated, and data quality flags where the coupling mechanism is not understood.

During the GW150914 analysis period, data quality triggers were applied as category 3 by burst searches. Times during hardware injection tests were also flagged and removed from the transient searches.

Modeled CBC searches, which use matched filtering techniques [2], apply additional mitigation methods to target loud noise transients with a duration on the order of a second or less that are particularly damaging. An accurate power spectral density (PSD) estimate is required to calculate the amount of signal power that matches a template waveform. Consequently, noise transients with a large amount of broadband power can corrupt the analyzed data up to the duration of the strain-equivalent noise PSD estimate, ±8 s from the time of the noise transient. Additionally, a loud, short-duration noise transient can act as a delta function, which may imprint the impulse response of the matched filter on the output data, generating triggers. As a result, before analyzing the data the CBC searches apply a technique called gating that smoothly rolls the input data stream off to zero for short-duration excursions identified as too loud to be consistent with an astrophysical signal [2].

The data set used to calculate the significance of GW150914 is appropriate in both the stability of the search backgrounds over the analysis period and the judicious application of data quality vetoes.

5.1. Stability of the period analyzed for GW150914

To illustrate the level of variability of detector performance over the several weeks of data collected for the analyzed time, figure 4 shows the maximum sensitive distance of each of the detectors for the coalescence of a binary black hole system with the same spin and mass parameters as GW150914 in the detector frame ($70$ ${M}_{\odot }$, $0.7$). This is calculated as the distance from Earth at which the coalescence of a binary object pair produces an SNR of 8 in a single detector using matched filtering, assuming optimal sky location and source orientation. LIGO-Hanford had a mean maximum sensitive distance to GW150914-like signals of 1906 Mpc during the analysis period, and LIGO-Livingston had a mean of 1697 Mpc.

Zoom In Zoom Out Reset image size

LIGO-Hanford's maximum sensitive distance exhibited a 90% range of ∼1800–2000 Mpc, and LIGO-Livingston's a 90% range of ∼1500–1900, which was sufficiently stable to provide a reliable estimate of the CBC search background throughout the analysis period. These small variations are due to a variety of fluctuations in the detectors and their environment, such as optic alignment variations or changing low frequency ground motion. Figure 5 shows the single-interferometer background trigger rate over time for the PyCBC search [7] with two different thresholds on the detection statistic, ${\chi }^{2}$-weighted SNR¹⁴¹ [2, 30, 31]. Triggers with a ${\chi }^{2}$-weighted SNR $\geqslant$ 6.5 (shown in green) comprise the bulk of the distribution and indicate the overall trigger rate from the search: ∼ 1–10 Hz. Triggers with ${\chi }^{2}$-weighted SNR $\geqslant$ 8 (shown in blue) are fairly rare, typically showing up at a rate < 0.01 Hz during the analysis period.

Zoom In Zoom Out Reset image size

The burst search background was also stable throughout the analysis containing GW150914. Figure 6 shows the behavior of background triggers from the coherent all-sky burst search cWB (coherent WaveBurst) [32, 33] during the analysis period. In contrast to the single-interferometer CBC triggers shown in figure 5, the coherent burst search requires coherent signal between multiple detectors to produce triggers, so the cWB background distribution is generated using time-shifted data. The main features of the background remain constant throughout the analyzed six weeks, particularly the domination of lower frequency triggers. Week 6 shows a small excess of triggers, ∼ 3% of total triggers, at lower than 60 Hz, which is below the majority of the power in event GW150914.

Zoom In Zoom Out Reset image size

Variations in the environmental conditions and instrumental state throughout the analysis time, as captured in the range variation seen in figure 4, did not have a significant impact on the PyCBC or cWB background distributions.

5.2. The impact of data quality flags on the transient searches

Data quality flags were generated independently for each detector in response to instrumental problems that demonstrated a well-defined, repeatable correlation with transient noise in h(t). Figure 7 shows the CBC background trigger distributions from each detector with and without data quality products applied. The LIGO-Hanford background distribution was dramatically improved by the application of data quality vetoes, dominated by the effect of a single data quality flag. This flag was designed to indicate a fault in the phase modulation system used to create optical cavity control feedback signals, as discussed in appendix A. LIGO-Livingston exhibits a longer tail of unvetoed background events which is largely composed of the blip noise transients discussed in section 3. The total time removed from the CBC search by vetoes is summarized for each detector by veto category in table 1.

Zoom In Zoom Out Reset image size

Table 1.  The deadtime introduced by each data quality (DQ) veto category, as discussed in section 4, for the CBC search during the analyzed period for LIGO-Hanford (left) and LIGO-Livingston (right).

	Hanford
DQ veto	Total	% of total
category	deadtime (s)	coincident time
1	73446	4.62%
2	5522	0.35%

	Livingston
DQ veto	Total	% of total
category	deadtime (s)	coincident time
1	1066	0.07%
2	87	0.01%

For GW150914, the reported false-alarm probability was not significantly affected by these data quality vetoes. GW150914 was the loudest recovered event during the analysis period—significantly louder than every background event even without data quality products applied.

For less significant triggers, the application of data quality vetoes is more important [34]. As an example, the false-alarm probability of the second most significant trigger (LVT151012) was 2%. Without the inclusion of data quality vetoes, the false-alarm probability would have been 14%, increased by roughly a factor of 7.

Figure 8 shows the impact of data-quality vetoes on the coherent burst search background, as well as the signal-consistency cut that requires resolved signals to have a time–frequency morphology consistent with expected astrophysical sources [3]. The data quality flag with the highest efficiency-to-deadtime ratio for the coherent burst search background indicated large excursions in h(t). This effective veto was defined using digital-to-analog overflows of the optic motion actuation signal used to stabilize the differential arm motion of the interferometer. This veto removed three of the loudest cWB background triggers during the analysis period. The remaining outliers with vetoes applied are blip-like noise transients of unknown instrumental origin.

Zoom In Zoom Out Reset image size

The total coincident time removed by each veto category from the burst search is summarized for each detector in table 2. Category 1 was defined identically between the burst and CBC searches, but there were some differences in the definition of category 2 largely due to differences in the observed impact of individual data quality products on the searches. For example, the CBC search used a data quality flag indicating periods of excess 10–30 Hz ground motion at LIGO-Hanford at category 2, but it was not applied to the burst search because it did not have a significant impact. The coherent burst search also applied a set of data quality triggers [25] at category 3, whereas the CBC search did not find this data quality product effective in reducing the background. A complete description of all data quality vetoes applied to the transient searches during the analysis period is reported in [35].

Table 2.  The deadtime introduced by each data quality (DQ) veto category for the coherent burst search during the analyzed period for LIGO-Hanford (left) and LIGO-Livingston (right).

	Hanford
DQ veto	Total	% of total
category	deadtime (s)	coincident time
1	73446	4.62%
2	1900	0.12%
3	12815	0.81%

	Livingston
DQ veto	Total	% of total
category	deadtime (s)	coincident time
1	1066	0.07%
2	736	0.05%
3	1319	0.08%

Figure 9 shows the effect of data quality vetoes on Omicron triggers from each detector. Since flags are tuned for specific problems at each detector, the impact on single-detector Omicron triggers is much more apparent than on the coherent burst search background in figure 8, where the search requirement of a high degree of signal correlation between multiple detectors is effective in reducing the background.

Zoom In Zoom Out Reset image size

Figure 9(a) shows that the same category 1 data quality veto that dominated the reduction in the LIGO-Hanford CBC background distribution only impacted noise transients up to an SNR of roughly 100. The higher SNR Omicron triggers vetoed at category 2 from both detectors are mostly large excursions in h(t) that are witnessed by overflows in the digital-to-analog conversion of the actuation signal controlling major optics, as mentioned for a data quality flag used effectively at category 2 for the coherent burst search. Blip noise transients are the main contributor to the unvetoed high SNR tail at both detectors along with 60–200 Hz non-stationarity that was persistent throughout the analysis period at LIGO-Livingston with an undetermined instrumental coupling.

The GW150914 event produced a strong gravitational wave signal in the Advanced LIGO detectors that shows the expected form of a binary black hole coalescence, as shown in figure 10 [1, 36]. Immediately around the event the data are clean and stationary.

Zoom In Zoom Out Reset image size

Even though the routine data quality checks did not indicate any problems with the data, in-depth checks of potential noise sources were performed around the time of GW150914. Potential noise couplings were considered from sources internal to the detector and local to each site, as well as common, coincident sources external to the detectors. All checks returned negative results for any pollution or interference large enough to have caused GW150914. Activities of personnel at the detectors, both locally and via remote internet connections, were confirmed to have no potential to induce transient noise in h(t). Because GW150914 occurred during the early morning hours at both detectors, the only people on-site were the control room operators. Signs of any anomalous activity nearby and the state of signal hardware injections were also investigated. These checks came back conclusively negative [37]. No data quality vetoes were active within an hour of the event. Rigorous checks of the data calibration were also performed [38].

The results of a key subset of checks intended to demonstrate nominal detector performance, quiet environment behavior, and clean data quality around the event are reported here.

For example, the U S Geological Survey [39] reported two magnitude 2.1 earthquakes within 20 min of GW150914; one with an epicenter off the coast of Alaska and another 70 miles south–west of Seattle. The earthquakes produced minimal vertical ground motion at 0.03–0.1 Hz at the time of arrival; roughly 10 nm s⁻¹ as measured by local seismometers at both detectors, which is an order of magnitude too small to produce an impact on the detector data.

6.1. Checks for potentially coincident noise sources

6.2. Checks of auxiliary channels for noise coincident with GW150914

6.3. Vetting of channels with identified excess power near the event time

A by-eye examination of spectrograms of every auxiliary channel identified a small subset of auxiliary channels that exhibited excess power within one second of GW150914, however, we found no evidence of noise that could generate GW150914 at either detector. In addition to the magnetometer events discussed above in relation to potentially coincident sources, there were 4 excess power events identified in magnetometers that monitor electromagnetically noisy electronics rooms. The observed magnetic fields would have had to have been at least 20 times stronger to account for the amplitude of GW150914 through coupling to the electronics. Channels from a seismometer and an accelerometer at LIGO-Hanford and two accelerometers at LIGO-Livingston also exhibited excess power. These vibrational disturbances were at least 17 times too small to account for the amplitude of GW150914. None of the environmental events matched GW150914 in time and frequency behavior.

The excess power triggers in the seismometer channels at LIGO-Hanford were likely due to a nearby air compressor with degraded vibration isolation that was running about 100 m away from optical components during the detection of GW150914. This excess ground motion, shown in figure 11, lasted for approximately three minutes at multiples of about 14 Hz (28, 42, 56 Hz). During the second containing GW150914, the largest disturbance detected by the seismometer (at ∼56 Hz) was at least 30 times too small to account for the amplitude of GW150914.

Zoom In Zoom Out Reset image size

There was also excess noise in the Livingston input mode cleaner [6] that was ruled out as a potential indication of noise that might mimic GW150914. This noise had time–frequency morphology that was inconsistent with any potential coupling mechanism. In particular, all power was below 8 Hz and the noise duration was nearly one second. Such a long transient would be unlikely to couple from the input mode cleaner to h(t) with duration comparable to GW150914 (∼200 ms).

6.4. Investigation of noise transients with similar morphology to CBC waveforms

Both detectors occasionally record short noise transients of unknown origin consisting of a few cycles around 100 Hz, including blip noise transients, discussed in section 3. None have ever been observed to occur in coincidence between detectors and follow-up examination of many of these transients confirmed an instrumental origin. While these transients are in the same frequency band as the candidate event, they have a characteristic time-symmetric waveform with significantly less frequency evolution, and are thus clearly distinct from the candidate event.

To illustrate this, figure 12 shows a blip transient that produced one of the most significant CBC background triggers associated with blip transients (${\chi }^{2}$-weighted SNR ≳ 9; compare to figure 7) during the analysis period and the neutron-star-black-hole (NSBH) binary template waveform it most closely matched. Although these noise transients do have significant overlap with regions of the CBC parameter space that produce very short waveforms, such as very high total mass binaries with extreme anti-aligned spins, they do not have a time domain morphology that matches CBC templates with similar character to GW150914.

Zoom In Zoom Out Reset image size

The potential impact of any accidental coincidence between such noise transients on the sensitivity of the searches is accounted for in the reported background distribution. No noise transients identified to have similar morphology elements to CBC signals [53], including blip transients, produced nearly as high a ${\chi }^{2}$-weighted SNR as GW150914.

6.5. LVT151012

GW150914 was by far the most significant event in all transient search results over the sixteen days of analyzed data. The CBC search also identified the second most interesting event on 12 October 2015. This trigger most closely matched the waveform of a binary black hole system with masses ${23}_{-6}^{+18}$ ${M}_{\odot }$ and ${13}_{-5}^{+4}$ ${M}_{\odot }$, producing a trigger with a false-alarm rate of 1 event per 2.3 years; far too high to be a strong detection candidate [1, 2, 54].

We performed similar in-depth checks of potential noise sources for this trigger. For LIGO-Livingston data, LVT151012 is in coincidence with significant excess power at 10 Hz lasting roughly three seconds, a portion of which can be seen in figure 13. There is no obvious indication of upconversion to the frequency range analyzed by the transient searches, so the low frequency noise is not thought to have caused the signal associated with LVT151012 in the Livingston detector.

Zoom In Zoom Out Reset image size

The data around this event were found to be significantly more non-stationary than those around GW150914. The noise transient rate in the hours around LVT151012 was significantly higher than usual at both LIGO detectors, seen in the Omicron trigger rate even on a broad time scale for LIGO-Livingston in particular, as illustrated in figure 14. This was likely due to increased low frequency ground motion associated with ocean waves [55]. The elevated noise transient rate at both sites induced a higher rate of background triggers around the time of LVT151012.

Zoom In Zoom Out Reset image size

No detector characterization studies to date indicate that LVT151012 was caused by a noise artifact.

6.6. Noise transient rate

At the time of GW150914, the LIGO detectors were operating in a low-noise state with nominal environmental and instrumental noise. Following the event, the detectors were maintained in the same configuration to ensure that detector changes would not cause unanticipated consequences which might bias the background estimation for the event. The backgrounds measured by the transient searches were stable throughout this analyzed period. Data quality vetoes were produced for each detector in response to instrumental or environmental noise sources. We conclude that the selected analysis period provides an accurate estimation of the significance of GW150914.

Additionally, thorough investigations found no evidence that environmental influences or non-Gaussian detector noise at either LIGO site might have caused the observed gravitational wave signal GW150914. A detailed study of environmental influences conclusively ruled out all postulated potential sources of correlated detector output at the time of the event, except for a binary black hole gravitational wave signal.

Characterization of the LIGO detectors via investigations of noise types that most impact the astrophysical searches and mitigation of noise couplings will continue to play a critical role in gravitational wave astronomy. Reducing the rate of high-significance background events and increasing search sensitivity is particularly important for near-threshold events such as LVT151012. Detector characterization will effectively expand the range of astrophysical sources that the gravitational wave detectors are sensitive to, providing a significantly greater number, and perhaps also variety, of events from which we can draw confident physical inferences.

The authors gratefully acknowledge the support of the United States National Science Foundation (NSF) for the construction and operation of the LIGO Laboratory and Advanced LIGO as well as the Science and Technology Facilities Council (STFC) of the United Kingdom, the Max-Planck-Society (MPS), and the State of Niedersachsen/Germany for support of the construction of Advanced LIGO and construction and operation of the GEO600 detector. Additional support for Advanced LIGO was provided by the Australian Research Council. The authors gratefully acknowledge the Italian Istituto Nazionale di Fisica Nucleare (INFN), the French Centre National de la Recherche Scientifique (CNRS) and the Foundation for Fundamental Research on Matter supported by the Netherlands Organisation for Scientific Research, for the construction and operation of the Virgo detector and the creation and support of the EGO consortium. The authors also gratefully acknowledge research support from these agencies as well as by the Council of Scientific and Industrial Research of India, Department of Science and Technology, India, Science & Engineering Research Board (SERB), India, Ministry of Human Resource Development, India, the Spanish Ministerio de Economía y Competitividad, the Conselleria d'Economia i Competitivitat and Conselleria d'Educació Cultura i Universitats of the Govern de les Illes Balears, the National Science Centre of Poland, the European Commission, the Royal Society, the Scottish Funding Council, the Scottish Universities Physics Alliance, the Hungarian Scientific Research Fund (OTKA), the Lyon Institute of Origins (LIO), the National Research Foundation of Korea, Industry Canada and the Province of Ontario through the Ministry of Economic Development and Innovation, the Natural Science and Engineering Research Council Canada, Canadian Institute for Advanced Research, the Brazilian Ministry of Science, Technology, and Innovation, Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), Russian Foundation for Basic Research, the Leverhulme Trust, the Research Corporation, Ministry of Science and Technology (MOST), Taiwan and the Kavli Foundation. The authors gratefully acknowledge the support of the NSF, STFC, MPS, INFN, CNRS and the State of Niedersachsen/Germany for provision of computational resources.

A data quality veto is generally constructed using an auxiliary channel which is strongly correlated with an instrumental problem. A notable example from the analyzed period was observed at LIGO-Hanford; intermittent periods with a significantly elevated transient noise rate in h(t). This behavior began suddenly five days before GW150914, independent of any activities taking place on site. The behavior was traced back to the 45 MHz electro-optic modulator driver system used to generate optical cavity control feedback signals [6]. To find the auxiliary channel which best correlates with non-stationary data in h(t), auxiliary channels recording interferometric cavity readouts and control signals associated with this driver were examined for excursions coincident with h(t) noise transients. A channel monitoring amplitude fluctuations in the signal used to generate the 45 MHz optical sidebands was found to be the best indicator of this non-stationary behavior.

Spikes in this auxiliary channel correlate well with a high rate of noise transients seen in h(t). However, the mean value of this channel varies significantly over time, meaning a simple threshold on the timeseries was not suitable for defining a data quality veto. Instead, band-limited root mean square values of this witness channel over minute strides were used. The effectiveness of different thresholds was tested using an 11 d subset of the analysis period. An example of the behavior of this veto over a 6 h time scale can be seen in figure A1 . With the selected threshold, this data quality veto removed 56% of noise transients with a SNR > 20, while only introducing 3% of deadtime over the 11 d of data. Figure A2 shows the distribution of Omicron triggers identified and removed, over the 11 days, by this veto.

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

This data quality flag was applied as a category 1 veto to the transient gravitational wave searches, responsible for removing $2.62 \%$ of the total coincident time from the analysis period.

The environment can influence the detector by mechanical force, electromagnetic waves, static electric and magnetic fields, and possibly high-energy radiation from cosmic rays. Mechanical forces, due to ground motion, temperature fluctuations, or air pressure fluctuations, are transmitted through structures that house and support interferometer optics and other key instrumentation.

Certain global-scale environmental effects could influence both detectors within 10 ms, which is the light travel time between the LIGO detectors and the maximum time delay for a gravitational wave signal of astrophysical origin. A network of sensors is employed such that global-scale environmental disturbances that could influence the detectors, such as electromagnetic disturbances in the atmosphere or transient fluctuations in the power grid, are redundantly monitored using PEM sensors that are significantly more sensitive to these disturbances than the detectors themselves.

By monitoring the immediate environment for disturbances that can be transmitted to the detector strain signal, we cover a large variety of environmental effects that can influence the detector data. For example, wind can couple through vibrations in the ground and air, and its behavior is witnessed by seismometers, accelerometers, and microphones (audio and infrasound frequencies). Lightning could couple by magnetic fields and electromagnetic waves at frequencies that we demodulate into the detection band for optic cavity control [6] and is monitored by magnetometers and RF receivers.

Figure B1 shows how these sensors are distributed at key locations throughout the LIGO-Livingston detector site (the LIGO-Hanford layout is very similar) [10]. Each building is equipped with seismometers and ground tilt sensors to monitor the motion of the concrete slab on which vacuum chambers and optical tables are mounted. Each of these buildings also contains an infrasound microphone and a set of audio-frequency microphones, including a microphone near the electronics that control the detector feedback loops and acquire auxiliary channel data. Power voltage monitors are installed in the electronics room of each building. Fluxgate magnetometers sense disturbances in the local magnetic field in all electronics rooms as well as a nearby subset of vacuum chambers. Accelerometers are mounted on vacuum chamber walls as well as on in-air optics tables and the concrete slab of each building. External to the detector buildings are RF receivers as well as wind speed sensors and outdoor weather stations. The PEM system at the Hanford detector includes a cosmic ray detector located underneath one of the test masses.

Zoom In Zoom Out Reset image size

There are a total of 173 PEM channels at LIGO-Hanford and 130 at LIGO-Livingston, where a greater number of channels at Hanford is due to additional redundancy in sensors as well as the cosmic ray detector.