---
title: LIGO/Virgo Sky Map Probabilities
url: https://www.emergentmind.com/topics/ligo-virgo-probability-sky-map
type: topic
---

# LIGO/Virgo Sky Map Probabilities

A LIGO/Virgo probability sky map is a discretized representation of the posterior probability distribution over the celestial sphere for the location of a gravitational wave (GW) source detected by the Advanced LIGO and Virgo observatories. These sky maps are the central data product for coupling gravitational-wave detection with electromagnetic follow-up campaigns. They encode, on a pixelized grid, the Bayesian probability that a GW source (typically a compact binary merger such as binary neutron star, BNS, or black hole mergers) lies in a given direction, and are derived from differences in arrival time, amplitude, and phase across the detector network, marginalized over all other source parameters. The construction, interpretation, and use of these sky maps are governed by both the statistical methodology underlying GW parameter estimation and the practical constraints of low-latency multi-messenger astronomy.

## 1. Bayesian Framework and Statistical Formalism

LIGO/Virgo sky maps are built via Bayesian inference, wherein the joint posterior distribution for all GW source parameters—including sky location $\Omega = (\alpha, \delta)$, intrinsic (masses, spins), extrinsic (distance, inclination), and waveform parameters—given the detector data $d$, is constructed as
$$
p(\theta|d) \propto L(d|\theta)\,\pi(\theta)
$$
where $L(d|\theta)$ is the likelihood under typically Gaussian noise assumptions and $\pi(\theta)$ encodes astrophysical or uninformative priors [1304.0670, 1312.6013, 1309.3273, 2510.21930]. For sky mapping, one focuses on the marginal posterior
$$
p(\Omega|d) = \int p(\theta|d)\,d\lambda
$$
where $\lambda$ are nuisance parameters marginalized out. The likelihood involves matched-filtering the strain data in each detector against parameterized templates; the detector response is a linear combination of $h_+$ and $h_\times$ polarisations weighted by antenna pattern functions $F_+$, $F_\times$ and shifted by sky-position-dependent time delays.

Posterior evaluation is achieved via either:
- **Fully coherent Bayesian stochastic sampling** (e.g., Markov Chain Monte Carlo or nested sampling, as in LALInference), marginalizing over all parameters, yielding highest fidelity but $\gtrsim 10$-hour latency [1309.3273, 1312.6013, 2110.01874].
- **Low-latency analytic or semianalytic approximation** (e.g., BAYESTAR), where sky location and distance are explored after fixing intrinsic parameters at their detection-pipeline maximum likelihood values; timing, amplitude, and phase matches are modeled by Gaussian-residual likelihoods [1404.5623, 2510.21930, 1312.6013, 2110.01874].

Bayesian priors are generally isotropic on the celestial sphere, volumetric in distance ($\propto D^2$), and astrophysically motivated for mass, spin, and orientation [1309.3273, 2105.04263].

## 2. Sky Map Pixelization, Posterior Assignment, and Credible Regions

Posterior distributions are rendered on a HEALPix grid, dividing the sphere into $N_{\rm pix}$ equal-area pixels of solid angle $\Omega_{\rm pix} = 4\pi/N_{\rm pix}$. The posterior probability in pixel $p$ is
$$
P_p \approx p(\Omega_p|d)\, \Omega_{\rm pix}
$$
with normalization $\sum_p P_p = 1$ [1404.5623, 1304.0670]. High throughput is crucial for rapid follow-up; BAYESTAR produces these maps in $O(10\,{\rm s})$–$O(1\,{\rm min})$ using summary statistics (arrival times $t_i$, SNRs $\rho_i$, arrival phases $\phi_i$) [1404.5623, 2510.21930].

To construct a credible region at level $\alpha$ (e.g., 90%), pixels are ranked by $P_p$ and summed until cumulative probability $\sum_{i=1}^{k} P_{(i)} \geq \alpha$ is reached. The area of this region is $A_\alpha = k\cdot\Omega_{\rm pix}$ [1404.5623, 1312.6013]. This "smallest-area" construction ensures, over repeated detections, frequentist coverage and is the LIGO/Virgo standard for reporting localization precision.

## 3. Sky Map Algorithms: Fully Coherent vs. Low-Latency Approaches

Two principal algorithmic families exist for producing sky maps:

**A. Fully coherent Bayesian parameter estimation** (high-latency): Employs stochastic sampling (MCMC or nested sampling) to sample the full joint posterior, including all correlations between sky position and other parameters. Post-sampling, marginalization yields a cloud of $(\alpha_j, \delta_j)$, which are binned into HEALPix sky maps [1309.3273, 1312.6013, 1105.3184]. Median 90% area for BNS in a three-detector network is $\sim 10$–$13.5\,{\rm deg}^2$ for optimal events; with four sites (e.g., adding India or Australia), $A_{90}$ drops to $\sim 2$–$6\,{\rm deg}^2$ [1309.3273, 1105.3184].

**B. Triangulation-based, rapid approximations** (low-latency): BAYESTAR evaluates analytic likelihoods based on timing, amplitude, and phase via a Fisher-matrix-based Gaussian for each detector. No waveform regeneration occurs, enabling $<1$ minute latency at the cost of larger, often arc-shaped regions (median $A_{90} \sim 200$–$600\,{\rm deg}^2$ in early HLV networks) [1404.5623, 2510.21930, 1312.6013]. Semianalytical approaches, such as that in [2110.01874], reduce the dimensionality of necessary integrations, achieving similar performance with further speed improvements.

The low-latency products are broadcast for electromagnetic follow-up (e.g., via GCN), and subsequently replaced as high-fidelity, full-sampling maps become available [1404.5623, 2105.04263].

## 4. Network Dependence, Degeneracies, and Statistical Diagnostics

Localization precision is strongly dependent on network geometry and sensitivity. In a two-detector Hanford–Livingston (HL) configuration, sky maps exhibit a pronounced “180° degeneracy” (antipodal arcs corresponding to equal arrival-time difference between the sites), and many maps are bimodal [1404.5623]. Addition of a third site (e.g., Virgo in HLV) breaks degeneracies and enables smaller, unimodal regions [2510.21930, 1404.5623].

However, the benefit depends on the SNR in the additional detector: if Virgo’s SNR is $\lesssim 5$, its inclusion can actually *degrade* the searched area in $\sim$14–20% of cases due to noise-induced likelihood mislocalization [2510.21930]. Therefore, several diagnostic metrics are recommended:
- **Jaccard Index ($J_{90}$)** for overlap between different network configurations’ 90% contours.
- **Kullback-Leibler divergence ($D_{\rm KL}$)**, quantifying whether added detectors sharpen or simply shift the posterior.
- **Percentile–Percentile plots and KS metrics**, testing calibration by comparing nominal coverage to actual inclusion rates.
- **Growth in credible region area with network size** as an anomaly flag.

Automatic diagnostic flags—such as nonoverlapping contours or area growth—are essential for rapid vetting, particularly in O3+ runs [2510.21930].

## 5. Three-Dimensional Sky Maps, Galaxy Targeting, and Volume Reduction

Contemporary rapid pipelines reconstruct not only 2D positions but also conditional luminosity distance posteriors for each sky pixel. The resulting product is a 3D posterior $P(\Omega, r|{\rm data})$, encoded via a per-pixel ansatz $p(r|n̂)=\frac{N(n̂)}{\sqrt{2\pi}\sigma(n̂)}\exp[-\frac{(r-\mu(n̂))^2}{2\sigma(n̂)^2}]\,r^2$ [1603.07333]. HEALPix-FITS products then provide (for each pixel) the sky probability $\rho_i$ and distance moments $(\mu_i, \sigma_i, N_i)$. This enables rapid computation of credible volumes, posterior per unit distance, or probabilistic ranking of host galaxy candidates.

To optimize electromagnetic follow-up, these 3D maps are cross-multiplied with galaxy catalog density (e.g., 2MASS Photometric Redshift Survey) within the reconstructed distance slice. The simple Bayesian update
$$
P_{\rm combined}(\Omega) \propto P_{\rm GW}(\Omega) \cdot N_{\rm gal}(\Omega)
$$
yields a reweighted probability map where $N_{\rm gal}(\Omega)$ is the line-of-sight integral over cataloged galaxy densities [1602.07710]. This approach typically halves or better the required follow-up area, number of galaxies to tile, and total telescope exposure time versus 2D-only tiling [1603.07333, 1602.07710].

## 6. Reported and Expected Localization Performance Across Network Runs

The localization area $A_{90}$ evolves as sensitivity and network size improve. Reported values include:

| Epoch / Network           | Median $A_{90}$ (BNS, deg$^2$)         | Notes                                          |
|--------------------------|-----------------------------------------|------------------------------------------------|
| 2015 (HL)                | 545 (Bayestar); 500+ (full PE) [1404.5623] | Elongated arcs, $\sim$180$^\circ$ degeneracy   |
| 2016 (HLV)               | 621 (Bayestar); 235 (full PE) [1404.5623] | Virgo breaks degeneracy; $\sim$2$\times$ shrink |
| O3 (HLV, full-sens)      | 270 (median, expected) [1304.0670]         |                                                |
| O4 (HLVK)                | 33 (median, expected) [1304.0670]          | Addition of KAGRA                              |
| O5 (HLVKI)               | $\sim$3–5 (median, forecast) [1304.0670]   | Global, 4–5-site network                       |
| Ideal full-coherence     | $\sim$5–13.5 ($68$–$95\%$) [1309.3273]     | SNR=20, HLV, full MCMC                         |

Probability sky maps routinely include diagnostic flags for reliability and provide both 2D and 3D data products to optimize multi-messenger astronomy. Localization areas scale approximately as $(\mathrm{SNR}\,\sigma_f)^{-2}$ with detector network SNR and signal bandwidth [1304.0670].

## 7. Limitations, Caveats, and Future Prospects

Several limiting assumptions warrant attention. The standard prior assumes GW hosts trace the general galaxy distribution, neglecting further astrophysical priors (such as host stellar mass, star-formation rate, or metallicity) [1602.07710]. Photometric redshift uncertainties in galaxy catalogs can limit 3D cross-matching precision, motivating the use of volume-complete spectroscopic surveys for nearby events. In areas obscured by Galactic extinction, zone-of-avoidance inpainting is used, but small hidden groups cannot be individually reconstructed. Algorithmic advances (reduced-order quadrature, prior optimization, analytical marginalization) continue to reduce latency and computational cost, with future low-latency 3D sky maps expected to accompany every routine GW alert with full reliability diagnostics [2110.01874, 2105.04263, 2510.21930].

A plausible implication is that as the detector network expands, the expected rapid localization area for typical BNS events will further decrease to a few square degrees, supporting real-time targeted electromagnetic counterpart searches and full-sky population studies.

Source: https://www.emergentmind.com/topics/ligo-virgo-probability-sky-map