Visualizing the Number of Existing and Future Gravitational-wave Detections from Merging Double Compact Objects

How many gravitational-wave observations from double compact object mergers have we seen to date? This seemingly simple question surprisingly yields a somewhat ambiguous answer that depends on the chosen data-analysis pipeline, detection threshold, and other underlying assumptions. To illustrate this we provide visualizations of the number of existing detections from double compact object mergers by the end of the third observing run (O3) based on recent results from the literature. Additionally, we visualize the expected number of observations from future-generation detectors, highlighting the possibility of up to millions of detections per year by the time next-generation ground-based detectors like Cosmic Explorer and Einstein Telescope come online. We present a publicly available code that highlights the exponential growth in gravitational-wave observations in the coming decades and the exciting prospects of gravitational-wave (astro)physics.


INTRODUCTION
Since the first detection of gravitational waves (GWs) from a binary black hole merger in 2015 (Abbott et al. 2016a), GW observatories Advanced LIGO (Aasi et al. 2015) and Advanced Virgo (Acernese et al. 2015) have rapidly increased the number of observed double-compact object mergers.The exact number of confident detections from double-compact object mergers to date, however, depends on the data-analysis pipelines used as well as underlying assumptions such as the detection or signal-to-noise threshold.For future observation runs and next-generation detectors, different expected detection numbers can arise from different assumptions on the detector capabilities or the underlying population properties of double compact object mergers in our Universe.
As we will discuss, there are a number of feasible choices for claiming the bona fide detection of a GW signal.We therefore spend the remainder of this manuscript shedding light on the number of detections in existing and future GW observation runs in two different ways.First, we outline a pedagogical overview of important concepts regarding the number of detected GW events in catalogs before exploring different choices present in the literature and their consequences.We additionally discuss the future growth of GW catalogs as we enter later generations of GW detector networks, and present a publicly available code to visualize the number of GW sources for outreach and research fsb2127@columbia.edusharan.banagiri@u.northwestern.eduepayne@caltech.eduarXiv:2303.17628v2[astro-ph.HE] 23 Apr 2024 purposes.We conclude with an outlook beyond the intricacies of event counting in the coming decades of future detector networks.

INTRODUCTION TO NUMBER OF DETECTIONS AND CATALOG CONSTRUCTION
When determining detection significance, a number of possible measures and criteria exist -a number of which have been adopted over the last decade.The most straight-forward calculation is to utilize the matched-filter signal-to-noise ratio (SNR) (Cutler & Flanagan 1994), to quantify a detection.Here, h(f ) is the Fourier transform of the measured strain, µ(f ) is the Fourier transform of the template signal model and S n (f ) is the power-spectral density of the interferometric noise.In the limit of Gaussian noise and "perfect" signal templates spanning the entire parameter space, the matched-filter SNR is the optimal detection statistic.Additionally, the noise distribution of SNRs follows a reduced χ 2 distribution (Allen et al. 2005) -providing an analytical expression of the distribution for the significance of an excursion from Gaussianity.However, GW signals are not the only signal present to cause non-Gaussianity (e.g., Akutsu et al. 2021;Davis et al. 2021).Glitches due to instrumental noise, seismic motion, and light scattering (to name a few), reduce the viability of a simple SNR threshold when analyzing real GW data (though such thresholds are often used to estimate detection probabilities when running a full search pipeline is infeasible, e.g., Abadie et al. 2010).This is the motivating factor behind many of the flagship searches, where significant effort is invested in determining the behavior of the background of noise transients.We point to Abbott et al. (2020a) for a guide to LVK data treatment and analysis methods.
Real search algorithms (e.g., Klimenko et al. 2016;Messick et al. 2017;Nitz et al. 2018;Sachdev et al. 2019;Aubin et al. 2021) 1 therefore employ different data-reduction techniques while analyzing the GW and output triggers which could potentially be statistically significant.Associated with each trigger is some ranking statistic x.The searches then compute the likelihood ratio between the probability of a trigger being attributed to a GW signal and the probability of being a noise transient, which can then be translated into a false-alarm rate (FAR).The FAR is a frequentist statistic -a close analogue of the p-value -which measures the frequency at which something as signal-like as the trigger can be generated by instrumental noise alone.A lower FAR hence implies a trigger that is statistically more significant.
FAR is a statistic that is particularly well-suited for the first detections, when the population of sources has not entirely come into focus.When we have several confident events that have mapped out regions of the parameter space, it is desirable that we also fold in information on the true alarm rate of the signals as well.The statistic developed to fold in this information is called p astro , defined as (Farr et al. 2015), where x is the trigger and S and N stand for the signal and noise hypotheses, respectively.Therefore, p astro is simply the probability that the trigger is astrophysical.We note that S assumes not just a signal model (e.g.compact object mergers) but also a specific distribution of sources as that impacts the likelihood of triggers under the signal hypothesis.
Using p astro allows us to distinguish between "vanilla" events which are expected and those that stand out as outliers when building population statistics.It also allows for a more apples-to-apples comparison between events that occur at different observing runs and hence at differing sensitivities (see Fig. 6 and the discussion in Sec.IV C in Abbott et al. 2023b).On the other hand, since p astro depends on some signal model S, in regions of the parameter space with a sparse set of detections, it could simply return back our prior assumptions in the population model.As a conservative choice in such cases, one can fall back on using a FAR threshold that is low enough to ensure minimal contamination.Once a few detections are obtained in such a region, the signal distribution can be inferred from them.

Label
Total physics and astronomy to use a 5σ threshold on the p value, corresponding to a probability of 1 in 3 × 10 −7 .With the GW community, one common standard is to use a p astro = 0.5 threshold as the dividing line to pick candidates for detailed analysis (Abbott et al. 2016b(Abbott et al. , 2019a(Abbott et al. , 2023b(Abbott et al. , 2024)).Whatever the choice of the threshold, there would be the possibility of a false alarm in the number of events that pass it.Another way to count the number of detections is to simply sum the p astro values, whose mean should converge to the total number of detections, as was done, for instance, in Fig. 2 of the GWTC-2.1 paper (Abbott et al. 2024).
In addition to detection, classification of candidates into source types also requires adopting some (populationinclusive) signal model.Catalogs usually calculate probabilities of different sources-types, which can be ambiguous for some candidates.In such cases, we adopt the source type with the maximum probability for the visualizations used in this paper.

CURRENT DETECTIONS FROM O1, O2, AND O3
At the time of writing, the LIGO-Virgo-KAGRA (LVK) Collaboration has completed and shared the data and detections from observing runs O1, O2, O3a, and O3b in the Gravitational-wave Transient Catalog (GWTC) releases GWTC-1 (Abbott et al. 2019a), GWTC-2 (Abbott et al. 2021a), GWTC-2.1 (Abbott et al. 2024), and GWTC-3 (Abbott et al. 2023b).O1 ran from September 12th, 2015 to January 19th, 2016, O2 ran from November 30th, 2016 to August 25th, 2017, O3a ran from 1 April 2019 until 1 October 2019, and O3b ran between 1 November 2019 and 27 March 2020.Based on the data in these finished runs, several catalogs have been presented in the literature with different numbers of events, based on different underlying assumptions for the detection pipelines and thresholds.We discuss the main catalog options considered in the literature below and summarize them in Table 1 and Figure 1.

LVK GWTC
The LVK Collaboration published several cumulative event catalogs based on observing runs O1, O2, and O3 using a number of search pipelines.As the field of GW astronomy evolved and the number of detections increased, these catalogs were presented using different statistics and thresholds.In GWTC-1, which contains candidates from both O1 and O2 runs, a preliminary FAR threshold of 1 per 30 days in at least one search pipeline and a p astro > 0.5 threshold were applied.GWTC-1 presented 3 binary black hole (BHBH) mergers in O1, one detection of a binary neutron star (NSNS) in O2, and 7 more BHBH mergers in O2 which passed this threshold (Abbott et al. 2019a).
In GWTC-2, the LVK reported first a total of 39 compact object coalescence detection candidates from the O3a run, including the second NSNS merger, using a FAR threshold of 2 year −1 in at least one pipeline (Abbott et al. 2021a).While no separate p astro threshold was used, the LVK calculated the astrophysical probability of all BHBH candidates  1.The exact number of detections depends on the underlying assumptions such as the chosen threshold of assigning gravitationalwave events to be of (likely) astrophysical origin.The animation shows the number of detections as a function of time and is available at Broekgaarden (2024).
and found that they were greater than 50%.The LVK collaboration reported an updated, deeper list of candidates from O3a in GWTC-2.1 using a FAR threshold of 2 day −1 in at least one pipeline (Abbott et al. 2024).While 1201 candidates passed this FAR threshold, only 44 of them had a p astro > 0.5, including one NSNS.In GWTC-2.1,LVK found that three GW events from GWTC-2 fall now below a p astro > 0.5 threshold but eight additional candidates were added, due to improved calibration, change of threshold, and subtraction of excess noise.
In GWTC-3, the LVK Collaboration presented an update of the GWTC with detections from O3b.They found 1083 GW candidates which passed a preliminary FAR threshold of 2 day −1 of which 35 have a p astro > 0.5, including three black hole-neutron stars (BHNSs).All in all, this gives a total of 90 events, with 85 BHBHs, 3 BHNSs, and 2 NSNSs that have an astrophysical probability > 0.5.The combined LVK GWTC is presented and kept up-to-date as the GWTC by the Gravitational Wave Open Science Center (GWOSC) (Abbott et al. 2023c;LVK Scientific Collaboration et al. 2023) 2 .Noteworthy, while the BHNS merger GW200105 162426 is included in the publication presenting the first two detections of a BHNS by the LVK Collaboration (Abbott et al. 2021b) as well as the GWTC-3 population inference catalogs (which use a FAR instead of a p astro threshold; see below), it was found to have an astrophysical probability of < 0.5.Future BHNS detections will help characterize the properties of the BHNS population and might address whether GW200105 162426 is likely of astrophysical origin.We also note that in the LVK GWTC (and all other catalogs mentioned below), GW190814 and GW200210 092254 with secondary masses of m 2 ≈ 2.5 − 3 are most likely classified as BHBHs as the secondary masses are above the hypothesized maximum mass of a neutron star (Essick & Landry 2020;Huang et al. 2020;Most et al. 2020;Shao et al. 2020;Tsokaros et al. 2020;Godzieba et al. 2021;Lim et al. 2021;Tews et al. 2021).For this reason we also classify both as BHBHs for the visualizations, but the possibility exists that they are instead BHNS sources (Abbott et al. 2020b(Abbott et al. , 2023a,b),b).
It is common that GW-related studies use subsets of the LVK GWTC, and in the context of the number of GW detections two subsets presented in the LVK population analysis Abbott et al. (2023a) are often used in the literature.First, in Abbott et al. (2023a), the LVK Collaboration presents the population properties of compact object mergers based on the data from O1, O2, and O3.The authors use an subset of events from the LVK GWTC based on the chosen false-alarm threshold of < 0.25 yr −1 to assure high purity, particularly among mergers with a neutron star 3 .This omits several candidates of moderate significance including GW190403 051519 and GW200220 061928.At this threshold LVK expects approximately one event to not be of astrophysical origin.They present the population catalog in their Table 1, consisting of 63 BHBH, 2 BHNS, and 2 NSNS mergers.Second, in the same paper the authors instead also discusses the inclusion of all events with FAR < 1 yr −1 , which is a slightly less strict threshold.At this threshold LVK expects approximately 4.6 events to not be of astrophysical origin.Using this threshold, LVK presents this population in Table 1 with 76 events consisting of 70 BHBH, 4 BHNS, and 2 NSNS mergers (Abbott et al. 2023a).This population catalog is used by the authors to analyze the redshift evolution of BHBH mergers.
3.2.4-OGC: independent catalog events with p astro > 0.5 In Nitz et al. (2023), the authors present an independent data analysis of the O1, O2, and O3 data with slightly different data analysis pipelines.Specifically, both the LVK analyses and Nitz et al. (2023) use the PyCBC opensource package, albeit with different settings, such as template bank construction.Compared to the earlier OGC catalogs (1-OGC, 2-OGC, and 3-OGC), the authors find seven more BHBH events compared to GWTC-3, and also find that three earlier BHBH events now pass the threshold for detection, based on updated assumptions.We take all the detections with p astro > 0.5 (which all have an inverse false alarm rate > 100 years -an additional threshold mentioned by the authors) from Table 2 in Nitz et al. (2023), which includes a total of 94 observed events with 90 BHBH, 2 BHNS 4 , and 2 NSNS mergers.
3.3.Optimistic catalog that combines events from the LVK GWTC, OGC and IAS analyses with p astro > 0.5) In addition to the catalogs presented by the LVK Collaboration and by the OGC Collaboration, the IAS Collaboration has also developed independent analysis of the public LIGO-Virgo data and presented additional events.The IAS pipeline focuses on specializing their pipeline to find new events instead of publishing full catalogs.For the purpose of visualization, we decide to make an 'optimistic' catalog that consists of all events from the LVK GWTC (GWOSC), 4-OGC, and IAS studies, that have an astrophysical origin of > 0.5.
This includes in O1 one additional BHBH event 5 with astrophysical probability p astro ∼ 0.71 (Venumadhav et al. 2019;Zackay et al. 2019).In O2 three additional BHBH events with likely astrophysical origin (p astro ≥ 0.98) have been reported by Venumadhav et al. (2020).These events were found independently with p astro ≥ 0.7 in Nitz et al. (2020).Three other BHBH events with astrophysical probabilities between 0.5-0.8 were also reported in Venumadhav et al. (2020), but with astrophysical probabilities < 0.25 by 2-OGC in (Nitz et al. 2020), but then updated to astrophysical probabilities of 0.37, 0.86, and 0.72 in Nitz et al. (2023).In Nitz et al. (2023), the authors presented with 4-OGC seven additional events with p astro ≥ 0.5 which we also add to the catalog.We also include GW170817A, which has been reported to have an astrophysical probability of 0.86 (Zackay et al. 2021).The same study found the potential gravitational-wave source GWC170402, but we do not include it, as the event cannot fully be described by a binary black hole waveform (see App. F of Abbott et al. (2023b) for an extended discussion on finding events not described by a compact binary coalescence waveform).We also do not include the intermediate mass binary black hole candidate GW170502, which is considered a marginal trigger in all literature as far as we are aware. 6 In Olsen et al. (2022) and Mehta et al. (2023) the authors re-analyze the Hanford and Livingston coincident data and find 10, and 5 new BHBH events with p astro > 0.5, in the O3a and O3b data, respectively.Among these, their events GW190704 104834, GW190910 012619, and GW200316 235947 could also be BHNS events, but we include 3 Where the low-significance event GW190531 was not included, lacking parameter inference results for use in hierarchical analyses (Abbott et al. 2023a): GW190531 was found in GWTC-2.1 with a FAR of 0.41 yr −1 in the GstLAL pipeline, passing the threshold of a FAR of less than 2 per year, but had pastro < 0.5, so no parameter inference was performed (Abbott et al. 2024), leading to its exclusion in the LVK population analysis even though this uses a FAR threshold of 1 yr −1 (Abbott et al. 2023a). 4Namely, GW200105 162426 and GW200115 042309. 5Where we note that this event, GW151205, with pastro = 0.53 was reported by (Nitz et al. 2020)  them here as BHBH sources.In Wadekar et al. (2023), the authors analyze LIGO-Virgo O3 data using an improved IAS pipeline that includes higher harmonics in the GW templates and down-weighting noise transients to improve the search sensitivity to high-mass and high-redshift merger events.They find 14 new BHBH events with p astro ≥ 0.5 that we include.We refer readers to the tables in, e.g., Olsen et al. (2022); Mehta et al. (2023); Wadekar et al. (2023) for overviews of the different p astro values between the OGC, GWTC and IAS analysis.We combine all additional events mentioned in this paragraph with the LVK GWTC (GWOSC) and 4-OGC catalogs to create a new 'Optimistic' catalog containing a total of 125 events with 119 BHBH, 4 BHNS, and 2 NSNS.

FUTURE DETECTIONS
We show a visualization of the possible expected observations in the coming decades in Figure 2. We emphasize that these are very rough estimates for the detection number and uncertainty.We provide publicly available code that allows the user to change these assumptions.For outreach purposes, it can be useful to show a video with a linear y-axis because it can be difficult for people to correctly perceive exponential growth (Wagenaar & Sagaria 1975;Wagenaar & Timmers 1979) leading people to vastly underestimate exponential growth (Levy & Tasoff 2017).We provide a video with a linear y-scale online and show a snapshot of the video in Figure 3.

Expected Detections for O4 and O5
The next upcoming LVK observing runs are O4 in 2023 and O5 around 2027.O4 started on May 24, 2023, and it is planned to have a total of 18-calendar month observing run, and LIGO, Virgo, and KAGRA detector sensitivities for binary neutron stars of about 160 Mpc, 40-80 Mpc, and 1-10 Mpc respectively7 .O5 will start around 2027 and is planned to run for about two years with a targeted LIGO, Virgo, and KAGRA sensitivities of 240-325 Mpc, 150-260 Mpc, and 25-128 Mpc, respectively.
We discuss three recent publications reporting the number of detections expected for O4 and O5.These expected number of detections are based on an SNR threshold of 10, rather than the FAR or p astro for simplicity and ease of projecting onto future detector network sensitivities.This SNR threshold is expected to be equivalent to a FAR ≲ 1 per year (see the discussion about comparing O1, O2, and O3 injections in  3, 4, and 5 for BHBH, BHNS, and NSNS detections.We add these also to Table 2. Third, in Kiendrebeogo et al. (2023) the authors perform a detailed analysis for O4 and O5 with based on two different assumptions for the underlying BHBH, BHNS and NSNS population (called LRR and GWTC-3).We add their rates for O4 and O5 from their Table 3.
Table 2. Summary of expected annual detection numbers for O4 and O5.Acronyms are: FAR (False Alarm Rate), GWTC (Gravitational Wave Transient Catalog), GWOSC (Gravitational Wave Open Science Center).For simplicity we do not show the provided uncertainties that often exceed unity.References: (a) taken from Abbott et al. (2018) using their updated table from https://emfollow.docs.ligo.org/and scaling down by a factor 3.375 to convert to confident detections (representing SNR≳ 12; corresponding to a single-detector SNR threshold of 8).(b) Numbers are taken from Tables 3, 4, and 5 from Iacovelli et al. (2022).(c) Numbers are taken from Table 3 from Kiendrebeogo et al. (2023).

Label
Total BHBH BHNS NSNS Ref It is important to realize that the exact number of BHBH, BHNS, and NSNS mergers that will be detected in O4 and O5 is highly uncertain.This number will depend on many factors including the actual sensitivity and observing duty cycle of the instruments and the Poisson statistics of the events.Another important uncertainty is the underlying population properties.One example is the mass distributions of the double compact object merger populations.Since the detector sensitivity scales strongly with the source masses, the number of detected events is currently heavily biased towards more massive events, making the frequency and properties of lower-mass black holes more uncertain (Abbott et al. 2023a).This makes it challenging to estimate the number of detections for future observing runs that probe larger distances and therefore see more of the underlying astrophysical population.Another important uncertainty is the redshift evolution of the rate and masses of the compact object merger events.From theory and observations, it is thought that the rate and properties of the double-compact objects change at higher redshifts (e.g.van Son et al. 2022;Abbott et al. 2023a;Callister & Farr 2023).Moreover, different formation channels with different properties might be at play (Ng et al. 2021;Zevin et al. 2021).This can all impact the number of detections in future detectors.

Detections for O4
Given the many uncertain factors that can impact the number of detected events in future detectors, we decide to take a very simplistic approach for the purpose of the visualizations presented here.We caveat to the reader that these visualizations are to sketch a potential landscape of the number of expected gravitational-wave detections and that the actual expected number is still highly uncertain and requires detailed analysis.We instead simply use the LVK GWTC (GWOSC) catalog as a proxy for the ratio between BHBH, BHNS, and NSNS detections as 85 : 3 : 2, and assume a total number of 75 BHBH mergers to be somewhat representative for a 11-month observing run at O3 sensitivity 10 .To obtain the expected number of detections for O4 and O5, we then scale the O3 detection number to the larger detection volume by using the LIGO O4 and O5 NSNS ranges and multiplying them by the planned observing run periods of 18 months and 2 years for O4 and O5, respectively.We also assume a uniform distribution of sources in an Euclidean universe.This gives an extremely rough estimate for O4 and O5 for BHBH sources of N in O4 = 75 × (18/11) × (160/120) 3 = 291 total BHBH detections in 18 months, N in O5 = 75 × (24/11) × (282.5/120) 3 = 2135 total BHBH detections in 24 months, with 120 Mpc as proxy for the LIGO O3 mean NSNS distance, 160 Mpc as a proxy for the average LIGO NSNS distance during O4, and 282.5 Mpc being the mean of the estimated LIGO NSNS distance for O5 11 Doing the same for the BHNS and NSNS rates using 3 BHNS and 2 NSNS detections in O3 gives 10 (75) BHNS and 7 (50) BNS expected in O4 (O5).For the uncertainty, we simply multiply and divide these numbers by a factor of three.We note that the 10 Where 75 was chosen to exclude the 10 events from O1 and O2 11 See https://observing.docs.ligo.org/plan/,accessed on 17/01/2024; the NSNS ranges could be updated over time.
Table 3. Summary of annual detection numbers expected with next-generation ground-based detectors Einstein Telescope and Cosmic Explorer (and variations).Acronyms are: FAR (false alarm rate), GWTC (Gravitational Wave Transient Catalog), GWOSC (Gravitational Wave Open Science Center, ET (Einstein Telescope), CE (Cosmic Explorer) N/A (not applicable; authors did not mention rates for this flavor).References: (a) Taken from the minimum and maximum numbers for ET-B and CE from Table 2 in Baibhav et al. (2019).(b) Taken from Table 4 from Borhanian & Sathyaprakash (2022) based events with SNR ≳ 10.We take the lowest and highest (ECS) rates for the pessimistic and optimistic rates, relatively.(c) Taken from Tables 3, 4, and 5 from Iacovelli et al. (2022), where we use the 'Analysed' events that have SNR≳ 12.(d) Taken from the lowest (pessimistic) and highest (optimistic) bounds in Section 5.4 from Iacovelli et al. (2022).(e) Taken from Table 6 in Ronchini et al. (2022).(f) Taken from column 4 (events with SNR≳ 10) in Table IV in Gupta et al. (2023).We include error bars and take the lowest and highest values for the pessimistic and optimistic detection numbers, respectively.1.8 × 10 6 1.5 × 10 5 4.4 × 10 5 1.2 × 10 6 (f) minimum 3.4 × 10 4 3.1 × 10 4 1.9 × 10 3 9.1 × 10 2 maximum 1.8 × 10 6 5.4 × 10 5 4.4 × 10 5 1.2 × 10 6 numbers quoted in Table 2 fall within our visualized uncertainty ranges, especially when including the uncertainties discussed in the studies (see Kiendrebeogo et al. 2023 for a detailed example of the effect of underlying assumptions).

Expected Detections for A ♯ , Cosmic Explorer, and Einstein Telescope
The exact sensitivity, and hence the expected number of detections that will be made with future upgrades to the LIGO and Virgo detectors beyond O5 is uncertain12 .In Baibhav et al. (2019); Adhikari et al. (2020); Gupta et al. (2023), the authors mention that a future best sensitivity of the LIGO/Virgo detectors, sometimes called A ♯ (or Voyager) might detect about a hundred times more sources compared to Advanced LIGO/Advanced Virgo/KAGRA (sometimes referred to as A + /AdV + ).For simplicity, we decide to use a factor of about 100 times more sources in A ♯ compared to our O3 estimate in Figure 2. We also include a factor 3 for the uncertainty estimate.
The number of detections expected with Cosmic Explorer (CE) and Einstein Telescope (ET) is also still highly uncertain.A source of uncertainty includes the (local) merger rate, which is used to calculate the underlying merger rate population.We include for the visualizations the following recent work calculating the expected number of detections.Gupta et al. (2023) present expected detection numbers for GW networks from A ♯ to networks with three next-generation detectors.We take their highest and lowest detection numbers for networks with at least one next-generation detector from their Table IV and add it to table 3 as a proxy for a CE/ET era detector.In Iacovelli et al. (2022) present expectations for the number of BHBH, BHNS, and NSNS detections with an ET or ET and two CE-like interferometer networks (see their Tables 3, 4, and 5).We add these rates to Table 3.In Section 5. 4 Iacovelli et al. (2022) discuss the impact on the detection number by changing the rate of the underlying population.We take the most optimistic and pessimistic numbers from these ranges and add this also to Table 3. Expected numbers with next-generation detectors have also been calculated by Borhanian & Sathyaprakash (2022); Ronchini et al. (2022), whose thorough analysis of the prospects for CE and ET are aligned with the bounds by Iacovelli et al. (2022) and Gupta et al. (2023), and we add their rates to Table 3. Baibhav et al. ( 2019) calculated expected detection rates for BHBH, BHNS, and NSNS for future detector networks based on population synthesis models.We add their minimum and maximum expected rates for an ET or CE detector to Table 3.The authors in these papers use slightly different models for the local merger rates, detection threshold (SNR), detector network, and underlying population, leading to the different answers, see appendix B in Iacovelli et al. (2022) for a detailed discussion.For our visualizations in Figure 2 we use the minimum and maximum detection numbers from the values above, as shown in the last two rows in Table 3.These numbers are aligned with expectations from earlier work such as Baibhav et al. (2019).

Caveats
For visualization purposes, we focused in this paper on simplified approximations for the (expected) detected number of GW sources.For example, the studies used for the visualization (Tables 1, 2, 3) often use different SNR or FAR thresholds and different underlying population properties -sometimes even within the same catalog.As a result, we note that the visualization in Figure 2 uses different thresholds for the existing and future detections.This has minimal impact on the visualization given the many other uncertainties.We provide an extra version of Figure 2 in our online repository that for the detections in O1, O2, and O3 uses GWTC-3 (GWOSC) events with 'matched-filter' SNR ≥ 10 (which contains 55 BBH, 1 BHNS, and 2 BNS mergers) instead of the p astro > 0.5 threshold in Figure 2 to align more with the threshold typically used in the studies for future detections (see also Section 2).Moreover, the detection analyses don't account for uncertainties in the population distribution of mergers.Additionally, with multiple pipelines analyzing the same underlying GW data, there are non-trivial trial factors that most analyses do not fully account for.There is significant ongoing work in the GW community on making these statistical foundations more robust; by better measuring our sensitivity to signals (Essick 2023;Essick & Fishbach 2023;Talbot & Golomb 2023), by combining and comparing different catalogs in a more self-consistent way (Sutton 2009;Biswas et al. 2012;Banagiri et al. 2023;Mould et al. 2023), and by development of techniques for simultaneous inference of detection statistics and population properties (Roulet et al. 2020;Galaudage et al. 2020).In addition, future work should investigate in more detail the effect on the number of detections from uncertain source classification (i.e.whether a detection classified as a BHBH could be a BHNS or NSNS).A larger number of detected events will help increase our understanding of the underlying merger population and guide individual source classification (Abbott et al. 2023a).

Future Outlook: beyond detector numbers
Future detectors such as A # , Einstein Telescope, and Cosmic Explorer will transform astronomy much beyond just increasing the number of detections, which has been the focus of this paper.While our understanding of the population improves as the square-root of the number of observations, these next-generation detectors provide data with everincreasing precision, at wider accessible gravitational-wave frequency bands, and out to larger distances, allowing detections of gravitational-wave sources possibly beyond redshifts of ten.This will allow us to address the redshift evolution of these sources, and probe compact objects and their progenitor (stars) in environments vastly different from the Milky Way.Moreover, an important development is that future detectors beyond increasing detection numbers will also detect many more events with higher precision (SNR), which will allow scientists to infer the source properties with more precision (see e.  2023)), which is important to address many open questions.For example.future detectors could determine the neutron star equation of state (e.g., Coupechoux et al. 2023), can provide early warning detections of binary neutron star mergers (e.g., Magee & Borhanian 2022), provide improved sky localization allowing for more multi-messenger observations (e.g., Mills et al. 2018;Petrov et al. 2022), might be able to observe the existence of black holes before star formation took place (e.g., Chen & Huang 2020;Ng et al. 2022), and can provide improved tests of General Relativity and estimates for the graviton mass (e.g.Perkins et al. 2021).We refer the reader to Reitze et al. (2019); Maggiore et al. (2020); Evans et al. (2021, 2023), and Gupta et al. (2023) and the references therein for more details and a detailed description of the science case for next-generation detectors.
All in all, we present in this paper publicly available visualizations to highlight the many detections to date in gravitational-wave astronomy as well as the expected detections in the exciting decades to come! 6. SOFTWARE AND THIRD PARTY DATA REPOSITORY CITATIONS All data and code to reproduce the results, visualizations, and figures in this paper are publicly available at Broekgaarden (2024).This code makes use of the publicly available GWOSC catalog from Abbott et al. (2023c);LVK Scientific Collaboration et al. (2023), the 4-OGC catalog from Nitz et al. (2023), and the GWTC-3 population catalogs based on Abbott et al. (2023a) retrieved from LVK Scientific Collaboration et al. (2021).The IAS data is based on Venumadhav et al. (2019Venumadhav et al. ( , 2020)); Zackay et al. (2019); Olsen et al. (2022); Mehta et al. (2023), andWadekar et al. (2023); see references in the papers for the repositories of the underlying IAS data.The code made use of Python from the Python Software Foundation available at http://www.python.org(van Rossum 1995).In addition, the following Python packages were used: Matplotlib (Hunter 2007), NumPy (Harris et al. 2020), SciPy (Virtanen et al. 2020), IPython/Jupyter (Perez & Granger 2007;Kluyver et al. 2016) FSB is grateful for the many suggestions and feedback received for this work from a list of people including Will Farr, Edo Berger, Carl-Johan Haster, Karan Jani, Alex Nitz, Soheb Mandhai, Eric Thrane, Ilya Mandel, Yvette Cendes, Maryam Hussaini, and Lieke van Son.All authors especially like to thank Christopher Berry for many discussions and feedback.We thank the anonymous referee for the insightful and valuable comments on the manuscript.FSB was supported by NASA FINESST scholarship 80NSSC22K1601, awarded to the President and Fellows of Harvard College.FSB is also grateful for the support on visualizing research received through the Harvard Horizons Scholarship program.FSB is also supported by Simons Foundation award 1141468 hosted at Columbia University.LIGO was constructed by the California Institute of Technology and Massachusetts Institute of Technology with funding from the National Science Foundation and operates under cooperative agreement PHY-1764464 .This paper carries LIGO Document Number LIGO-P2300435.This research has made use of data or software obtained from the Gravitational Wave Open Science Center (gwosc.org), a service of the LIGO Scientific Collaboration, the Virgo Collaboration, and KAGRA.This material is based upon work supported by NSF's LIGO Laboratory which is a major facility fully funded by the National Science Foundation, 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.Virgo is funded, through the European Gravitational Observatory (EGO), by the French Centre National de Recherche Scientifique (CNRS), the Italian Istituto Nazionale di Fisica Nucleare (INFN) and the Dutch Nikhef, with contributions by institutions from Belgium, Germany, Greece, Hungary, Ireland, Japan, Monaco, Poland, Portugal, Spain.KAGRA is supported by Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan Society for the Promotion of Science (JSPS) in Japan; National Research Foundation (NRF) and Ministry of Science and ICT (MSIT) in Korea; Academia Sinica (AS) and National Science and Technology Council (NSTC) in Taiwan.

Figure 1 .
Figure 1.Cumulative histogram of the number of detected mergers from binary black holes (BH-BH; black), black hole-neutron stars (BHNS; teal), and binary neutron stars (NS-NS; orange) for different data analysis pipelines in the literature from Table1.The exact number of detections depends on the underlying assumptions such as the chosen threshold of assigning gravitationalwave events to be of (likely) astrophysical origin.The animation shows the number of detections as a function of time and is available atBroekgaarden (2024).

Figure 2 .
Figure2.Landscape of the number of detected mergers from binary black holes (BHBH; black), black hole-neutron stars (BHNS; teal), and binary neutron stars (NSNS; orange).See text for more details.The animation shows the catalog as a function of time and is available atBroekgaarden (2024).

Figure 3 .
Figure3.Still of our animation that visualizes the growth of the gravitational-wave catalog in the coming decades on a linear scale.Available atBroekgaarden (2024).The linear y-scale emphasizes how the detections expected in the next decades will completely dwarf the number of detections observed today.