<blockquote>
  <p><em>Updated to use reported timing confidence intervals instead of trying
  to infer those from reported sky location resolution uncertainties.
  The latter approach was misleading because the sky location has been further
  resolved using more information than just the arrival time difference.</em></p>
</blockquote>

<p>Beware that I am neither an astronomer nor involved with aLIGO.
I just use classic euclidean geometry here,
combined with occasional transformations between cartesian
and spherical coordinates.</p>

<p>We have two observatories, $L$ (Livingston) and $H$ (Hanford)
and a far-away point $E$ emitting a (gravitational wave) signal.
Let us assume that the signal propagates at constant speed $c$
through low-density matter, including planet Earth.
The straight-line distance between $L$ and $H$ (passing through Earth)
is $|HL| \approx c\cdot 10.0\cdot10^{-3}\,\mathrm{s}$.</p>

<p>$L$ receives the signal from $E$ about
$t_{ELH}\approx 7\cdot10^{-3}\,\mathrm{s}$ before $H$.
This implies that $E$ is a distance
$d = c\,t_{ELH}$ closer to $L$ than to $H$.</p>

<p>Now assuming large-scale euclidean geometry and looking at the
plane containing $E,L,H$, the set of all candidate points $E$
where $|EH|-|EL| = d$ is one branch of a hyperbola, namely
the branch closer to $L$:</p>

<p><a href="https://i.sstatic.net/GDk39.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/GDk39.png" alt="Hyperbola"></a></p>

<p>In fact, setting
$$\begin{align}
    r &amp;= \frac{d}{2}    &amp; r_1 &amp;= \frac{|HL|}{2} &amp; R &amp;= |LE|
\\  L &amp;= (r_1,0)        &amp; H &amp;= (-r_1,0)         &amp; E &amp;= (x,y)
\end{align}$$
we get the cartesian hyperbola equation
$$\frac{x^2}{r^2} - \frac{y^2}{r_1^2-r^2} = 1$$
and in polar coordinates, viewed from $L$:
$$R = \frac{r_1^2 - r^2}{r - r_1\cos\alpha}$$
For very large $R/r_1$, $|\alpha|$ must be very close to its asympotic limit:
$$\begin{align}
    \cos\alpha &amp;\approx \frac{r}{r_1} = \frac{d}{|HL|} = \frac{t_{ELH}}{t_{HL}}
\\  \text{where}\quad
    t_{HL} &amp;= \frac{|HL|}{c}
    \approx 10.0\cdot 10^{-3}\,\mathrm{s}
    \qquad\text{(time directly from $H$ to $L$)}
\end{align}$$
In 3D this means that $E$ lies on one of the two sheets of a hyperboloid
of revolution whose axis of revolution is the line $HL$.
In particular, in $L$'s view of the sky, partially hidden under the horizon,
there should be a fairly thin circular band that covers all projected
candidate points for $E$.</p>

<p>The solid angle of view enclosed by the circular band is
$$  A_{\text{ste}} = 2\pi(1 - \cos\alpha)
\approx 2\pi\left(1 - \frac{t_{ELH}}{t_{HL}}\right)$$
Uncertainties in $t_{ELH}$ yield uncertainties in $\cos\alpha$
which widen the circular band.
Thus the circle becomes an annulus.
The solid angle of view covered by that annulus measures
$$\begin{align}
    A_{\text{ste.max}} - A_{\text{ste.min}}
    = 2\pi(\cos\alpha_{\text{min}} &amp;{}- \cos\alpha_{\text{max}})
    \approx 2\pi\frac{t_{ELH\text{.max}} - t_{ELH\text{.min}}}{t_{HL}}
\\  \text{briefly:}\quad
    \Delta A_{\text{ste}} &amp;\approx 2\pi\frac{\Delta t_{ELH}}{t_{HL}}
\end{align}$$
Substituting reported data for $t_{ELH}$, with $\text{min}$ and $\text{max}$
referring to 90% probability intervals, we have:
$$\begin{array}{rcrrr}
    \text{entity} &amp; \text{unit} &amp; \text{nom} &amp; \text{min} &amp; \text{max} &amp; \Delta
\\\hline
    t_{ELH} &amp; 10^{-3}\,\mathrm{s}   &amp; 6.9   &amp; 6.5   &amp; 7.4   &amp; 0.9
\\  \cos\alpha &amp; 1                  &amp; 0.69  &amp; 0.65  &amp; 0.74  &amp; 0.09
\\  90^\circ - \alpha &amp; {}^\circ    &amp; 44    &amp; 41    &amp; 48    &amp; 7
\end{array}$$
And thus the annulus, when projected to a unit (celestial) sphere, covers
$$\Delta A_{\text{ste}} \approx 0.18\pi \approx (1.9\cdot10^3)^{\circ\circ}$$
(The ${}^{\circ\circ}$ means
<a href="https://en.wikipedia.org/wiki/Square_degree" rel="nofollow noreferrer">square degrees</a>.)
In contrast, the reported value of the <em>credible region</em> varies between
$140^{\circ\circ}$ (50% probability) and $590^{\circ\circ}$ (90% probablity).
This is because large parts of the annulus could be ruled out based on
additional information such as signal strength ratios.
But those refinements are not the topic of this thread.</p>

<p>For easier visualization, let us move to a point $B$ on Earth's surface
where the center of the annulus is in the zenith,
and the annulus itself is at constant altitude $90^\circ - \alpha$,
as given in the above table.
To that end, let $O$ be the location of Earth's center and require
that the 3D ray $OB$ points in the same direction as the ray $HL$.
This determines the latitude and longitude of $B$.
Note that the parallel shift from $HL$ to $OB$ incurs a parallax,
but that is negligible because the distance to $E$ is so large.</p>

<p>Idealizing our planet to a sphere and doing the calculations
(transforming the locations of $H$ and $L$ from spherical coordinates to
cartesian coordinates, computing the difference vector, normalizing it,
and transforming back to spherical coordinates), I find $B$ to be at
$$B\colon\ 27^\circ 21'54.22''\mathrm{S}\ 38^\circ 36'33.45''\mathrm{W}
= 27.365061^\circ\mathrm{S}\ 38.609293^\circ\mathrm{W}$$
which is <a href="https://www.google.de/maps/place/27%C2%B021%2754.2%22S+38%C2%B036%2733.5%22W/@-27.0915311,-38.9169065,4z/data=!4m2!3m1!1s0x0:0x0" rel="nofollow noreferrer">some 700km southeast of Rio de Janeiro in the Atlantic</a>.
If you had been there at <code>2015-09-14 09:50:45 UTC</code>
and had spun around yourself, drawing a band at constant altitude
$90^\circ - \alpha \approx 44^\circ$ and width $7^\circ$ in the sky,
that would have marked all the directions
where the signal source $E$ could have been.</p>

<p>To get the sky view simulated, I have written a tiny script <code>gw150914.sts</code> for
<a href="http://nightshadesoftware.org/projects/nightshade-legacy/wiki" rel="nofollow noreferrer">Nightshade</a>:</p>

<pre><code># Display the sky at the time of the aLIGO GW150914 observation
# from an Earth-based position chosen such that the set of possible sky
# positions of the source form a circular annulus with azimuth-independent
# altitude = asin(delta_t/10ms), about 44 deg.
clear
timerate rate 0
date utc 2015-09-14T09:50:45
flag show_tui_datetime on
set home_planet "Earth"
moveto lat -27.365061 lon -38.609293 alt 0 heading 0 duration 2
wait duration 2
flag azimuthal_grid on
zoom fov 120 duration 2
</code></pre>

<p>This gives the view below. I have added dashed green circles to delimit
the annulus.</p>

<p><a href="https://i.sstatic.net/CF9KI.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/CF9KI.png" alt="Nightshade view with added indication of the annulus"></a></p>

<p>But do not expect to find likely candidate spots within that annulus.
The annulus hides the radial depth of the projected space volume.
A range for that radial depth is given by the estimated bounds
for the distance $R=|LE|$. The corresponding flat space volume is
$$\Delta V = \frac{1}{3}\Delta A_{\text{ste}}
\left(R_{\text{max}}^3 - R_{\text{min}}^3\right)$$
Plugging in the reported values, again from 90% probability intervals,
$$\begin{align}
\Delta A_{\text{ste}} &amp;\approx 590^{\circ\circ} \approx 0.0572\pi
&amp;   R_{\text{min}} &amp;\approx 0.23\cdot10^{9}\,\mathrm{pc}
&amp;   R_{\text{max}} &amp;\approx 0.57\cdot10^{9}\,\mathrm{pc}
\end{align}$$
yields
$$\Delta V \approx 10^{25}\,\mathrm{pc}^3$$
To get a very rough estimate of the number of galaxies in that volume,
let us multiply $\Delta V$ with a rough galaxy density estimate of $10^{-20}\,\mathrm{pc}^{-3}$,
which results in about $10^5$ galaxies in that volume.
My conclusion from this back-of-the-envelope calculation is
that we cannot hope to localize $E$ further, unless given more information.</p>

<p>Things would have been a bit different if we had had a third aLIGO observatory
distant from the line $HL$ operating and capable of detecting the signal
from $E$.
Then we could draw three annuli, take intersections and thus
determine the direction of the signal within one reasonably small spot
on the celestial sphere.</p>

<h2>Further reading</h2>

<ul>
<li>As <a href="https://www.ligo.caltech.edu/image/ligo20160211b" rel="nofollow noreferrer">this LIGO post</a> tells
us, parts of the annulus could be ruled out by considering additional features
of the data such as varying signal strength.</li>
<li>A more detailed description of the probability distributions for the results
is given in <a href="https://dcc.ligo.org/public/0122/P1500218/012/GW150914_parameter_estimation_v13.pdf" rel="nofollow noreferrer">this LIGO paper</a>.</li>
<li>More information about the methods used to further narrow the credible
sky regions of the source location is given
in <a href="http://dx.doi.org/10.1007/lrr-2016-1" rel="nofollow noreferrer">this paper</a>.</li>
</ul>

<h2>Addendum</h2>

<p>An equivalent script <code>gw150914.ssc</code>
for <a href="http://stellarium.sourceforge.net" rel="nofollow noreferrer">Stellarium</a>:</p>

<pre><code>// Display the sky at the time of the aLIGO GW150914 observation
// from an Earth-based position chosen such that the set of possible sky
// positions of the source form a circular annulus with azimuth-independent
// altitude = asin(delta_t/10ms), about 44 deg.

var obs_lat = -27.365061;
var obs_long = -38.609293;
var obs_alt = 0;
var obs_locname = "Atlantic, 700km SE from Rio de Janeiro";
var obs_planet = "Earth";
var obs_date = "2015-09-14T09:50:45";
var duration = 2;

core.clear("deepspace");
core.setTimeRate(0);
core.setDate(obs_date, "utc");
core.setObserverLocation(obs_long, obs_lat, obs_alt, duration,
                         obs_locname, obs_planet);
core.wait(duration);
core.moveToAltAzi("90", "0", duration);
core.wait(duration);
GridLinesMgr.setFlagAzimuthalGrid(true);
StelMovementMgr.zoomTo(120, duration);
</code></pre>

<p>And yes, <code>core.setObserverLocation</code> expects longitude before latitude.</p>

CC BY-SA 3.0
2016-03-20T13:37:10.827Z
I think there is a typo on your assumed galaxy density. Impressive nonethless. Perhaps you could summarise in a sentence or two what the signal strength considerations are that diminish the search area by a factor three. Note also that the 600 sq deg is not an absolute error box, it is the high probability region.
ProfRob
CC BY-SA 3.0
2016-03-20T13:48:10.173Z
@RobJeffries: In the beginning I took the data from the [1-page GW150914 fact sheet](https://losc.ligo.org/s/events/GW150914/GW150914-FactSheet-BW.pdf). There was no $\Delta t_{ELH}$ and the $\Delta A_{\text{ste}}$ was only specified to be ~$600^{\circ\circ}$. Then I came across a [more detailed report](https://dcc.ligo.org/public/0122/P1500218/012/GW150914_parameter_estimation_v13.pdf) where the actual probability densities were described. I have now updated the post to match the latter data, mostly using the 90% confidence intervals, as indicated by now. Yes, it's no error box.
ccorn
CC BY-SA 3.0
2016-03-20T13:58:32.943Z
In the text I have given the reciprocal of the galaxy density. The [more detailed report](https://dcc.ligo.org/public/0122/P1500218/012/GW150914_parameter_estimation_v13.pdf) uses a comoving galaxy density estimate of ~$10^7$ per cubicgigaparsec, which converted and reciprocated means about $3.5\cdot10^{21}\,\mathrm{ly}^3$ per galaxy. Mine is $4.2$ instead of $3.5$, so in the same ball park. But I'll update to use their density.
ccorn
CC BY-SA 3.0
2016-03-12T23:29:02.033Z
2016-03-20T16:16:18.883Z
2016-03-20T16:16:18.883Z
ccorn
<p>The addition of more detectors will immensly help event localization.  To get an understanding of the calculation look at the old LORAN navigation system.  This system used synchronized transmitters, and the receiver measured the difference in time.</p>

<p>The receiving device measured the time between signals from multiple sources, that chirp at the same time. If you hear the chirps at the same time, from two of the synchronized transmitters then you are on the plane equidistant between the sources.  The plane is perpendicular to the line between the sources, and at the halfway point.  If you plot the intersection of the plane, and the surface of the earth it creates a great circle.  You are on the great circle.  </p>

<p>Signals that have a difference in time of reception form a curved surface.  As it cuts through the earth it defines a curved line on the surface (hyperbolic).</p>

<p>If the signals come at different times then you are closer to one than the other.  The charts used for determining position had lines of constant time difference plotted.  You listened to one pair of transmitters, got your time difference, found the line on the chart that represents that delay.  Do this for another pair of transmitters, get that difference, find the line for that pair with that delay, and where the lines cross is about where you are.</p>

<p>This is repeated for multiple pairs, and finding all the intersections (the best two lines are those crossing at right angles to reduce the error possible in your timers).</p>

<p>This is the same thing in reverse.  </p>

<p>Get more listener pairs and we can start to isolate the region of space of the signal source.</p>

<p>I hope that helps.</p>

CC BY-SA 3.0
2016-02-13T19:40:17.320Z
2016-02-13T19:40:17.320Z
Michael F.
<p>The signal does not travel from Livingston to Hanford. The signal comes in with an angle of about 45 degrees to the line joining the two, and has the same amplitude on planes perpendicular to the direction of travel (the peaks of water waves are not along the direction of travel of the water waves but perpendicular to the direction of travel). Thus if the signal came in at 90 degrees of the line joining Hanford and Livingston, the peaks of the signal would arrive at both at the same time. If it came in in the direction along the path, then it would take the signal 10 ms later at Hanford than at Livingston. At 45 degrees, the peak at Livingston would arrive <span class="math-container">$\cos(45°) \times 10$</span>ms <span class="math-container">$= 7$</span>ms later at Hanford than at Livingston. Thus the direction of the signal is on a cone with apex angle very close to 45 degrees. That is all one could tell from the timing.</p>
<p>Hanford and Livingston detectors were build to be as close to parallel to each other as possible. If we lived on a flat earth, then the signal at both would be identical. The earth is curved. Thus the detector at Hanford is tilted from the straight line joining them by about 27 degrees with respect to Livingston. This means that the sensitivity of the detector is not the same for incoming directions all along the cone. In some directions, Hanford is about 50% of the sensitivity of Livingston while for others they are nearly the same. By comparing the amplitudes at Hanford and Livingston, one can get a very rough idea of where along that 45 degree cone the signal came from, limiting the directions to something like half the cone.</p>
<p>When Virgo came online, even though its sensitivity is much less than either of the LIGO detectors, that extra timing information is enough to strongly limit the directions the signal could have come from.</p>

CC BY-SA 4.0
2020-07-22T04:35:23.550Z
2020-07-22T07:15:15.537Z
2020-07-22T07:15:15.537Z
W Unruh
<p>On Feb 12 James wrote:"A difference of 7ms would imply a minimum separation of the sites of about 2000 km, not a maximum."</p>

<p>The calculated minimum separation of the sites is around 3001.65 km compared with the estimated maximum of 3030.12 km indicated by gollum that same day. The angular separation between Livingston and Hanford detectors is around 27.25 degrees from the center of our planet using  an average Earth's radius of 6371.0 Km and gives this result. If the detected GW150914 gravitational waves propagated outwards at more speed that the speed of light a difference of 7ms is in accord with this MINIMUM possible detector distance. </p>

<p>The GW's velocity on 2015/09/14/09:50:45 UT based on this argument would be:</p>

<p>(3001.65 km/2099.0 km)*c ~ 1.43 times the speed of light </p>

CC BY-SA 3.0
2016-04-08T20:14:27.487Z
2016-04-08T20:14:27.487Z
Juan Zapata-Arauco
<p>Here's how my non-scientist mind envisions it. I draw a straight line between the two LIGO sites on a map. Then I take another straight line (like a straight edge/ruler) that represents the GW coming in.  If the GW line is coming from the south exactly parallel with the line I drew on the map, then both sites would detect the "chirp" at exactly the same time (0 milliseconds difference). If the GW line comes from the south exactly perpendicular to my drawn line, hitting the LA site first, then chirp would be detected at each site at approx 10 milliseconds (?) apart (assuming that's how long it takes for light to travel between the two sites in a straight line).</p>

<p>So I envision the GW "line" came in from the south somewhere between parallel and perpendicular (closer to perpendicular) in relation to my drawn line. Someone else smarter than I can figure out the maths and correct angles, numbers, etc.</p>

<p>Does that make any sense?</p>

CC BY-SA 3.0
2016-02-12T14:24:51.733Z
Yes, it makes sense, and this is how they know an approximate direction — from the time lag between the two detections. Like the time lag between the detections of a sound by our ears allow our brain to calculate the direction to the cause of the sound.
pela
CC BY-SA 3.0
2016-02-12T13:39:38.797Z
2016-02-12T13:39:38.797Z
iMerchant
<p>Quote from <a href="https://www.ligo.caltech.edu/news/ligo20160211" rel="nofollow noreferrer">LIGO's news release</a>:</p>
<blockquote>
<p>By looking at the time of arrival of the signals—the detector in
Livingston recorded the event <strong>7 milliseconds</strong> before the detector in
Hanford—scientists can say that the source was located in the Southern
Hemisphere.</p>
</blockquote>
<p>GPS coordinates of LIGO Livingston: <a href="https://www.google.at/maps/place/LIGO+Livingston/@30.5630018,-90.7763949,17z/data=!3m1!4b1!4m2!3m1!1s0x8626e835b7d97141:0x7f7b0db13ffa9864" rel="nofollow noreferrer">30.5630018,-90.7763949</a><br />
GPS coordinates of LIGO Hanford: <a href="https://www.google.at/maps/place/LIGO+Hanford+Observatory/@46.4554032,-119.4092701,831m/data=!3m1!1e3!4m2!3m1!1s0x54985dc1a03257c3:0x70ac2e7c3e7fc470!6m1!1e1" rel="nofollow noreferrer">46.4554032,-119.4092701</a></p>
<p>Distance: <a href="http://boulter.com/gps/distance/?from=46.4554032%20-119.4092701&amp;to=30.5630018%20-90.7763949&amp;units=k&amp;sort=distance" rel="nofollow noreferrer">3030.13 km</a></p>
<p>Now if the gravitational waves propagate outwards at the speed of light a difference of 7ms implies the following maximum possible detector distance:</p>
<pre><code>299792458 * 0.007 / 1000 = 2099 km
</code></pre>
<p>Where am I wrong ?</p>
<p><strong>Update:</strong></p>
<p>As pointed out by <a href="https://astronomy.stackexchange.com/users/2365/">james-k</a> in the comments below, the calculated 2099 km is the minimum distance and not the maximum ... these waves seem to have caused some mental derangement ;)</p>

CC BY-SA 3.0
2016-02-12T14:08:00.527Z
A difference of 7ms would imply a minimum separation of the sites of about 2000km, not a maximum.
James K
CC BY-SA 3.0
2016-02-12T16:11:55.367Z
You can use the straight line, not great circle, distance between the two sites, assume gravity moves at c, and calculate a circle along which the source must lie. Too bad we didn't have 3 detectors up, or it'd be a line.
Wayfaring Stranger
CC BY-SA 3.0
2016-04-07T15:41:21.397Z
@AlbertBriffoz Nothing can travel faster than the speed of light, its the corner stone of special relativity and is yet to be proven wrong (even with the best efforts of the neutrino community! (https://en.wikipedia.org/wiki/Faster-than-light_neutrino_anomaly)
Dean
CC BY-SA 3.0
2016-04-07T16:40:02.800Z
@AlbertBriffoz Since you're new here, you should know that we work differently than forums. Please take the [Tour](http://astronomy.stackexchange.com/tour) and see the [Help Center](http://astronomy.stackexchange.com/help). If you're reading a post and have a follow up question, you should post it as a new question on the site using the [Ask Question](http://astronomy.stackexchange.com/questions/ask) link in the top right corner of the browser view. The answer section is not for open-ended discussion but specifically for direct answers to the question above.
called2voyage
<gravitational-waves><ligo>
Gravitational wave detection time difference between LIGO Livingston and LIGO Hanford
gollum
# Gravitational wave detection time difference between LIGO Livingston and LIGO Hanford

(Asked by: username_0 on 2016-02-12)

Quote from LIGO's news release:

By looking at the time of arrival of the signals—the detector in
Livingston recorded the event 7 milliseconds before the detector in
Hanford—scientists can say that the source was located in the Southern
Hemisphere.

GPS coordinates of LIGO Livingston: 30.5630018,-90.7763949
GPS coordinates of LIGO Hanford: 46.4554032,-119.4092701
Distance: 3030.13 km
Now if the gravitational waves propagate outwards at the speed of light a difference of 7ms implies the following maximum possible detector distance:
299792458 * 0.007 / 1000 = 2099 km

Where am I wrong ?
Update:
As pointed out by james-k in the comments below, the calculated 2099 km is the minimum distance and not the maximum ... these waves seem to have caused some mental derangement ;)

**Tags:** gravitational-waves, ligo

> **Comments:**
> **username_1**: A difference of 7ms would imply a minimum separation of the sites of about 2000km, not a maximum.  
> **username_2**: You can use the straight line, not great circle, distance between the two sites, assume gravity moves at c, and calculate a circle along which the source must lie. Too bad we didn't have 3 detectors up, or it'd be a line.  
> **username_3**: @AlbertBriffoz Nothing can travel faster than the speed of light, its the corner stone of special relativity and is yet to be proven wrong (even with the best efforts of the neutrino community! (https://en.wikipedia.org/wiki/Faster-than-light_neutrino_anomaly)  
> **username_4**: @AlbertBriffoz Since you're new here, you should know that we work differently than forums. Please take the [Tour](http://astronomy.stackexchange.com/tour) and see the [Help Center](http://astronomy.stackexchange.com/help). If you're reading a post and have a follow up question, you should post it as a new question on the site using the [Ask Question](http://astronomy.stackexchange.com/questions/ask) link in the top right corner of the browser view. The answer section is not for open-ended discussion but specifically for direct answers to the question above.  

---

## 5 answers

### Answer by username_10 on 2016-02-12. Score: 6 | Accepted answer

Here's how my non-scientist mind envisions it. I draw a straight line between the two LIGO sites on a map. Then I take another straight line (like a straight edge/ruler) that represents the GW coming in. If the GW line is coming from the south exactly parallel with the line I drew on the map, then both sites would detect the "chirp" at exactly the same time (0 milliseconds difference). If the GW line comes from the south exactly perpendicular to my drawn line, hitting the LA site first, then chirp would be detected at each site at approx 10 milliseconds (?) apart (assuming that's how long it takes for light to travel between the two sites in a straight line).

So I envision the GW "line" came in from the south somewhere between parallel and perpendicular (closer to perpendicular) in relation to my drawn line. Someone else smarter than I can figure out the maths and correct angles, numbers, etc.

Does that make any sense?

> **Comments:**
> **username_11**: Yes, it makes sense, and this is how they know an approximate direction — from the time lag between the two detections. Like the time lag between the detections of a sound by our ears allow our brain to calculate the direction to the cause of the sound.  

---

### Answer by username_5 on 2016-03-12. Score: 7

Updated to use reported timing confidence intervals instead of trying
  to infer those from reported sky location resolution uncertainties.
  The latter approach was misleading because the sky location has been further
  resolved using more information than just the arrival time difference.

Beware that I am neither an astronomer nor involved with aLIGO.
I just use classic euclidean geometry here,
combined with occasional transformations between cartesian
and spherical coordinates.
We have two observatories, $L$ (Livingston) and $H$ (Hanford)
and a far-away point $E$ emitting a (gravitational wave) signal.
Let us assume that the signal propagates at constant speed $c$
through low-density matter, including planet Earth.
The straight-line distance between $L$ and $H$ (passing through Earth)
is $|HL| \approx c\cdot 10.0\cdot10^{-3}\,\mathrm{s}$.
$L$ receives the signal from $E$ about
$t_{ELH}\approx 7\cdot10^{-3}\,\mathrm{s}$ before $H$.
This implies that $E$ is a distance
$d = c\,t_{ELH}$ closer to $L$ than to $H$.
Now assuming large-scale euclidean geometry and looking at the
plane containing $E,L,H$, the set of all candidate points $E$
where $|EH|-|EL| = d$ is one branch of a hyperbola, namely
the branch closer to $L$:

In fact, setting
$$\begin{align}
    r &= \frac{d}{2}    & r_1 &= \frac{|HL|}{2} & R &= |LE|
\\  L &= (r_1,0)        & H &= (-r_1,0)         & E &= (x,y)
\end{align}$$
we get the cartesian hyperbola equation
$$\frac{x^2}{r^2} - \frac{y^2}{r_1^2-r^2} = 1$$
and in polar coordinates, viewed from $L$:
$$R = \frac{r_1^2 - r^2}{r - r_1\cos\alpha}$$
For very large $R/r_1$, $|\alpha|$ must be very close to its asympotic limit:
$$\begin{align}
    \cos\alpha &\approx \frac{r}{r_1} = \frac{d}{|HL|} = \frac{t_{ELH}}{t_{HL}}
\\  \text{where}\quad
    t_{HL} &= \frac{|HL|}{c}
    \approx 10.0\cdot 10^{-3}\,\mathrm{s}
    \qquad\text{(time directly from $H$ to $L$)}
\end{align}$$
In 3D this means that $E$ lies on one of the two sheets of a hyperboloid
of revolution whose axis of revolution is the line $HL$.
In particular, in $L$'s view of the sky, partially hidden under the horizon,
there should be a fairly thin circular band that covers all projected
candidate points for $E$.
The solid angle of view enclosed by the circular band is
$$  A_{\text{ste}} = 2\pi(1 - \cos\alpha)
\approx 2\pi\left(1 - \frac{t_{ELH}}{t_{HL}}\right)$$
Uncertainties in $t_{ELH}$ yield uncertainties in $\cos\alpha$
which widen the circular band.
Thus the circle becomes an annulus.
The solid angle of view covered by that annulus measures
$$\begin{align}
    A_{\text{ste.max}} - A_{\text{ste.min}}
    = 2\pi(\cos\alpha_{\text{min}} &{}- \cos\alpha_{\text{max}})
    \approx 2\pi\frac{t_{ELH\text{.max}} - t_{ELH\text{.min}}}{t_{HL}}
\\  \text{briefly:}\quad
    \Delta A_{\text{ste}} &\approx 2\pi\frac{\Delta t_{ELH}}{t_{HL}}
\end{align}$$
Substituting reported data for $t_{ELH}$, with $\text{min}$ and $\text{max}$
referring to 90% probability intervals, we have:
$$\begin{array}{rcrrr}
    \text{entity} & \text{unit} & \text{nom} & \text{min} & \text{max} & \Delta
\\\hline
    t_{ELH} & 10^{-3}\,\mathrm{s}   & 6.9   & 6.5   & 7.4   & 0.9
\\  \cos\alpha & 1                  & 0.69  & 0.65  & 0.74  & 0.09
\\  90^\circ - \alpha & {}^\circ    & 44    & 41    & 48    & 7
\end{array}$$
And thus the annulus, when projected to a unit (celestial) sphere, covers
$$\Delta A_{\text{ste}} \approx 0.18\pi \approx (1.9\cdot10^3)^{\circ\circ}$$
(The ${}^{\circ\circ}$ means
square degrees.)
In contrast, the reported value of the credible region varies between
$140^{\circ\circ}$ (50% probability) and $590^{\circ\circ}$ (90% probablity).
This is because large parts of the annulus could be ruled out based on
additional information such as signal strength ratios.
But those refinements are not the topic of this thread.
For easier visualization, let us move to a point $B$ on Earth's surface
where the center of the annulus is in the zenith,
and the annulus itself is at constant altitude $90^\circ - \alpha$,
as given in the above table.
To that end, let $O$ be the location of Earth's center and require
that the 3D ray $OB$ points in the same direction as the ray $HL$.
This determines the latitude and longitude of $B$.
Note that the parallel shift from $HL$ to $OB$ incurs a parallax,
but that is negligible because the distance to $E$ is so large.
Idealizing our planet to a sphere and doing the calculations
(transforming the locations of $H$ and $L$ from spherical coordinates to
cartesian coordinates, computing the difference vector, normalizing it,
and transforming back to spherical coordinates), I find $B$ to be at
$$B\colon\ 27^\circ 21'54.22''\mathrm{S}\ 38^\circ 36'33.45''\mathrm{W}
= 27.365061^\circ\mathrm{S}\ 38.609293^\circ\mathrm{W}$$
which is some 700km southeast of Rio de Janeiro in the Atlantic.
If you had been there at 2015-09-14 09:50:45 UTC
and had spun around yourself, drawing a band at constant altitude
$90^\circ - \alpha \approx 44^\circ$ and width $7^\circ$ in the sky,
that would have marked all the directions
where the signal source $E$ could have been.
To get the sky view simulated, I have written a tiny script gw150914.sts for
Nightshade:
# Display the sky at the time of the aLIGO GW150914 observation
# from an Earth-based position chosen such that the set of possible sky
# positions of the source form a circular annulus with azimuth-independent
# altitude = asin(delta_t/10ms), about 44 deg.
clear
timerate rate 0
date utc 2015-09-14T09:50:45
flag show_tui_datetime on
set home_planet "Earth"
moveto lat -27.365061 lon -38.609293 alt 0 heading 0 duration 2
wait duration 2
flag azimuthal_grid on
zoom fov 120 duration 2

This gives the view below. I have added dashed green circles to delimit
the annulus.

But do not expect to find likely candidate spots within that annulus.
The annulus hides the radial depth of the projected space volume.
A range for that radial depth is given by the estimated bounds
for the distance $R=|LE|$. The corresponding flat space volume is
$$\Delta V = \frac{1}{3}\Delta A_{\text{ste}}
\left(R_{\text{max}}^3 - R_{\text{min}}^3\right)$$
Plugging in the reported values, again from 90% probability intervals,
$$\begin{align}
\Delta A_{\text{ste}} &\approx 590^{\circ\circ} \approx 0.0572\pi
&   R_{\text{min}} &\approx 0.23\cdot10^{9}\,\mathrm{pc}
&   R_{\text{max}} &\approx 0.57\cdot10^{9}\,\mathrm{pc}
\end{align}$$
yields
$$\Delta V \approx 10^{25}\,\mathrm{pc}^3$$
To get a very rough estimate of the number of galaxies in that volume,
let us multiply $\Delta V$ with a rough galaxy density estimate of $10^{-20}\,\mathrm{pc}^{-3}$,
which results in about $10^5$ galaxies in that volume.
My conclusion from this back-of-the-envelope calculation is
that we cannot hope to localize $E$ further, unless given more information.
Things would have been a bit different if we had had a third aLIGO observatory
distant from the line $HL$ operating and capable of detecting the signal
from $E$.
Then we could draw three annuli, take intersections and thus
determine the direction of the signal within one reasonably small spot
on the celestial sphere.
Further reading

As this LIGO post tells
us, parts of the annulus could be ruled out by considering additional features
of the data such as varying signal strength.
A more detailed description of the probability distributions for the results
is given in this LIGO paper.
More information about the methods used to further narrow the credible
sky regions of the source location is given
in this paper.

Addendum
An equivalent script gw150914.ssc
for Stellarium:
// Display the sky at the time of the aLIGO GW150914 observation
// from an Earth-based position chosen such that the set of possible sky
// positions of the source form a circular annulus with azimuth-independent
// altitude = asin(delta_t/10ms), about 44 deg.

var obs_lat = -27.365061;
var obs_long = -38.609293;
var obs_alt = 0;
var obs_locname = "Atlantic, 700km SE from Rio de Janeiro";
var obs_planet = "Earth";
var obs_date = "2015-09-14T09:50:45";
var duration = 2;

core.clear("deepspace");
core.setTimeRate(0);
core.setDate(obs_date, "utc");
core.setObserverLocation(obs_long, obs_lat, obs_alt, duration,
                         obs_locname, obs_planet);
core.wait(duration);
core.moveToAltAzi("90", "0", duration);
core.wait(duration);
GridLinesMgr.setFlagAzimuthalGrid(true);
StelMovementMgr.zoomTo(120, duration);

And yes, core.setObserverLocation expects longitude before latitude.

> **Comments:**
> **username_6**: I think there is a typo on your assumed galaxy density. Impressive nonethless. Perhaps you could summarise in a sentence or two what the signal strength considerations are that diminish the search area by a factor three. Note also that the 600 sq deg is not an absolute error box, it is the high probability region.  
> **username_5**: @RobJeffries: In the beginning I took the data from the [1-page GW150914 fact sheet](https://losc.ligo.org/s/events/GW150914/GW150914-FactSheet-BW.pdf). There was no $\Delta t_{ELH}$ and the $\Delta A_{\text{ste}}$ was only specified to be ~$600^{\circ\circ}$. Then I came across a [more detailed report](https://dcc.ligo.org/public/0122/P1500218/012/GW150914_parameter_estimation_v13.pdf) where the actual probability densities were described. I have now updated the post to match the latter data, mostly using the 90% confidence intervals, as indicated by now. Yes, it's no error box.  
> **username_5**: In the text I have given the reciprocal of the galaxy density. The [more detailed report](https://dcc.ligo.org/public/0122/P1500218/012/GW150914_parameter_estimation_v13.pdf) uses a comoving galaxy density estimate of ~$10^7$ per cubicgigaparsec, which converted and reciprocated means about $3.5\cdot10^{21}\,\mathrm{ly}^3$ per galaxy. Mine is $4.2$ instead of $3.5$, so in the same ball park. But I'll update to use their density.  

---

### Answer by username_7 on 2016-02-13. Score: 3

The addition of more detectors will immensly help event localization. To get an understanding of the calculation look at the old LORAN navigation system. This system used synchronized transmitters, and the receiver measured the difference in time.

The receiving device measured the time between signals from multiple sources, that chirp at the same time. If you hear the chirps at the same time, from two of the synchronized transmitters then you are on the plane equidistant between the sources. The plane is perpendicular to the line between the sources, and at the halfway point. If you plot the intersection of the plane, and the surface of the earth it creates a great circle. You are on the great circle.

Signals that have a difference in time of reception form a curved surface. As it cuts through the earth it defines a curved line on the surface (hyperbolic).

If the signals come at different times then you are closer to one than the other. The charts used for determining position had lines of constant time difference plotted. You listened to one pair of transmitters, got your time difference, found the line on the chart that represents that delay. Do this for another pair of transmitters, get that difference, find the line for that pair with that delay, and where the lines cross is about where you are.

This is repeated for multiple pairs, and finding all the intersections (the best two lines are those crossing at right angles to reduce the error possible in your timers).

This is the same thing in reverse.

Get more listener pairs and we can start to isolate the region of space of the signal source.

I hope that helps.


---

### Answer by username_8 on 2020-07-22. Score: 3

The signal does not travel from Livingston to Hanford. The signal comes in with an angle of about 45 degrees to the line joining the two, and has the same amplitude on planes perpendicular to the direction of travel (the peaks of water waves are not along the direction of travel of the water waves but perpendicular to the direction of travel). Thus if the signal came in at 90 degrees of the line joining Hanford and Livingston, the peaks of the signal would arrive at both at the same time. If it came in in the direction along the path, then it would take the signal 10 ms later at Hanford than at Livingston. At 45 degrees, the peak at Livingston would arrive $\cos(45°) \times 10$ms $= 7$ms later at Hanford than at Livingston. Thus the direction of the signal is on a cone with apex angle very close to 45 degrees. That is all one could tell from the timing.

Hanford and Livingston detectors were build to be as close to parallel to each other as possible. If we lived on a flat earth, then the signal at both would be identical. The earth is curved. Thus the detector at Hanford is tilted from the straight line joining them by about 27 degrees with respect to Livingston. This means that the sensitivity of the detector is not the same for incoming directions all along the cone. In some directions, Hanford is about 50% of the sensitivity of Livingston while for others they are nearly the same. By comparing the amplitudes at Hanford and Livingston, one can get a very rough idea of where along that 45 degree cone the signal came from, limiting the directions to something like half the cone.

When Virgo came online, even though its sensitivity is much less than either of the LIGO detectors, that extra timing information is enough to strongly limit the directions the signal could have come from.


---

### Answer by username_9 on 2016-04-08. Score: 0

On Feb 12 James wrote:"A difference of 7ms would imply a minimum separation of the sites of about 2000 km, not a maximum."

The calculated minimum separation of the sites is around 3001.65 km compared with the estimated maximum of 3030.12 km indicated by username_0 that same day. The angular separation between Livingston and Hanford detectors is around 27.25 degrees from the center of our planet using an average Earth's radius of 6371.0 Km and gives this result. If the detected GW150914 gravitational waves propagated outwards at more speed that the speed of light a difference of 7ms is in accord with this MINIMUM possible detector distance.

The GW's velocity on 2015/09/14/09:50:45 UT based on this argument would be:

(3001.65 km/2099.0 km)\*c ~ 1.43 times the speed of light


---

