GEOPHYSICS
Crowdsourced earthquake early warning
Sarah E. Minson,
1,2
Benjamin A. Brooks,
1
* Craig L. Glennie,
3
Jessica R. Murray,
1
John O. Langbein,
1
Susan E. Owen,
4
Thomas H. Heaton,
2
Robert A. Iannucci,
5
Darren L. Hauser
3
Earthquake early warning (EEW) can reduce harm to people and infrastructure from earthquakes and tsunamis, but it
has not been implemented in most high earthquake-risk
regions because of prohibitive cost. Common consumer
devices such as smartphones contain low-cost versions of the sensors used in EEW. Although less accurate than
scientific-grade instruments, these s
ensors are globally ubiquitous. Through controlled tests of consumer devices,
simulation of an
M
w
(moment magnitude) 7 earthquake on California
’
s Hayward fault, and real data from the
M
w
9 Tohoku-oki earthquake, we demonstrate that EEW could be achieved via crowdsourcing.
INTRODUCTION
Earthquake early warning (EEW) strives to detect an earthquake
’
sinitia-
tion, estimate its location and magnitude, and alert people and automated
systems to imminent shaking (
1
,
2
). EEW has had encouraging initial
results (
3
), although two issues impede performance and wider implemen-
tation. First, magnitude estimates are unreliable for larger earthquakes [mo-
ment magnitude (
M
w
)
>
7] when based solely on brief observations of the
earliest seismic waves (
4
). This can be overcome by using observations from
Global Navigation Satellite Systems (GNSS) such as the Global Positioning
System (GPS) (
5
). Second, it is expensive to install and operate the required
dense seismic and GNSS networks. Nevertheless, even well-monitored re-
gions such as California, Oregon, and Washington require extensive expan-
sion and upgrade of existing instrumentation, including installing
hundreds of new instruments, to implement EEW (
6
). Consequently, seis-
mic EEW is operational in a handful of regions, and only a few of those
(Japan, Mexico, and the United States) are incorporating GNSS data into
their systems (Fig. 1) (
7
). Much of the global population exposed to high
seismic risk, especially in poorer countries, does not benefit from EEW.
Commercial demand for personal m
obile navigation has led to a pro-
liferation of devices that use the same, albeit lower-quality, GNSS and In-
ertial Navigation Systems (INS) sensors used for EEW (
8
). Smartphones
alone currently number 1 billion worldwide and will increase to ~5.9 billion
by 2019 (
9
). The ubiquity of consumer devices raises the possibility that
operational EEW could be achieved via crowdsourcing (
10
–
12
). For a
global population exposed to ever
-increasing earthquake risk (
13
), the sig-
nificantly reduced costs associated with crowdsourcing could facilitate
widespread EEW implementation and substantially reduce the impact of
future earthquakes.
Although there has already been exploration of the potential of scien-
tific seismic and GPS data (
14
), as well as consumer-quality acceler-
ometers (
15
), for EEW and tsunami early warning, the potential of
consumer-quality GNSS receivers remains untapped (
16
–
18
). As we will
show, these data are rapidly improving in quality and are much less noisy
than consumer-quality accelerometers, at only a fraction of the cost of scientific-
quality instruments. Furthermore, GNSS data, which are displacement
observations, are particularly well suited to monitoring large earth-
quakes and to produce magnitudes that do not saturate. Finally, we
present a conceptual strategy for a
crowdsourced EEW system that in-
cludes device use, data processing a
nd quality control, earthquake de-
tection, false alarm suppression, earthquake location, and magnitude
determination that could be used to extend the benefits of EEW
worldwide.
We assess the potential performa
nce of smartphones and a crowd-
sourced EEW system in three ways. Fi
rst, we perform cont
rolled tests on
a variety of consumer devices to determine their noise character and
displacement detection capability. Second, for a scenario
M
w
7 rupture
on northern California
’
s Hayward fault, we generate synthetic smart-
phone accelerometer and GNSS time series for different data types
we might obtain from crowdsourcing at random locations based on
census population data. Finally, we
consider GPS position time series
of an actual earthquake, the
M
w
9.0 Tohoku-Oki event, obtained using
positioning data of the type found on consumer devices.
RESULTS
To obtain surface displacement observations for EEW, it is the
instrument
’
s change in position that must be measured accurately; ab-
solute position is used only for spati
al reference. Consumer devices typ-
ically use single-frequency, C/A (coarse acquisition) code methods for
GNSS positioning rather than the mo
re precise and accurate locations
derived from dual-frequency, carrier phase
–
based algorithms used in
scientific applications (
19
). These C/A code positions can be substantial-
ly improved by using differential corrections via satellite-based augmen-
tation systems (SBAS) (
19
), tracking the more precise GNSS carrier
phase and using it to filter the C/A code data (
“
phase smoothing
”
)
(
20
), or by combination with independent INS data in a Kalman filter
(
21
). Today
’
s smartphones have some or all of these capabilities.
Two data types represent the range of data most likely for crowd-
sourced EEW. First, we investigated the least-precise data recorded by
consumer devices: raw C/A code posi
tions and low-quality accelerom-
eter time series. These data can be u
sed alone or in combination using a
Kalman filtering approach (
21
). The Kalman filter takes advantage of
the fact that the accelerometer is recording the second derivative of
the same displacement time series that the GNSS receiver is recording,
and from these two data sources produces a unified estimate of
displacement with much less noise and bias than either of the original
time series. Second, we considered the most sophisticated GNSS receiver
commonly found in consumer devices: one capable of recording raw C/A
1
U.S. Geological Survey, Menlo Park, CA 94025, USA.
2
California Institute of Technology,
Pasadena, CA 91106, USA.
3
National Center for Airborne Laser Mapping, University of
Houston, Houston, TX 77204, USA.
4
Jet Propulsion Laboratory, La Cañada Flintridge, Pasadena,
CA 91109, USA.
5
Carnegie Mellon University
–
Silicon Valley, Moffett Field, CA 94035, USA.
*Corresponding author. E-mail: bbrooks@usgs.gov
2015 © The Authors, some rights reserved;
exclusive licensee American Association for
the Advancement of Science. Distributed
under a Creative Commons Attribution
NonCommercial License 4.0 (CC BY-NC).
10.1126/sciadv.1500036
RESEARCH ARTICLE
Minson
et al
. Sci. Adv. 2015;1:e1500036 10 April 2015
1of7
data as well as real-time SBAS differential corrections while applying
phase smoothing. To represent these scenarios (a poor GNSS receiver
supplemented by a poor accelerometer, and a good GNSS receiver),
we studied a latest-generation smart
phone (Google Nexus 5) that contains
a C/A code receiver and accelerometer, and a u-blox consumer GNSS re-
ceiver that is capable of recording SBAS and performing phase smoothing.
We then subjected these two devices, along with a scientific-grade INS, to a
series of displacements ranging from ~0.1 to 2.0 m (Fig. 2). The displace-
ment time series from the scientific-grade INS system may be considered
to be the true motion. Figure 2A shows the response of these two data
types, as well as the twice-integrated acceleration and C/A code time series
input to the Kalman filter, to a simple time series of motion. Both data types
reproduce the time history of displ
acement with very high fidelity.
To demonstrate the capability of co
nsumer-quality sensors to record
earthquake ground motions, we compare the noise from various
instruments to observed displacement time series for recent earth-
quakes (Fig. 2B) (
22
–
24
). Raw C/A code positions alone recorded by
any GNSS-equipped consumer device (light pink curve) are capable
of recovering displacements from great earthquakes. The light blue
curve represents observations from typical smartphones; generally,
smartphones do not report positions from raw C/A code data, but in-
stead use a Kalman filter to combine raw C/A code positions with the
device
’
s accelerometer. The filter operates on the GNSS chipset and is
optimized for consumer navigational needs, such as vehicular and ur-
ban canyon positioning, rather than recording higher-frequency
earthquake motions. If the same raw C/A code positions and accelera-
tion data recorded by a smartphone were combined in a Kalman filter
tuned to optimize earthquake surface displacements, the noise level
would be reduced to that shown by the cyan curve. Some consumer
devices are now beginning to use better GNSS positioning hardware
(SBAS-capable receivers, red curve) and algorithms (such as phase
smoothing, magenta curve). Either of these improvements reduces
noise so significantly that Kalman filtering is no longer required and
the detection threshold approaches smaller (M6-7) earthquakes.
We compute the minimum magnitude
earthquake observable with a
signal-to-noise ratio of at least 10 f
or all data types at different source-
receiver distances (Fig. 2C). These results show that consumer C/A code
GNSS data augmented by phase smoothing, SBAS, or Kalman filtering
with twice-integrated acceleration data could be used for measuring
ground displacement from ~
M
w
≥
6 earthquakes within ~100 km of
the source (
7
). Data from consumer accelerometers could be used to detect
smaller earthquakes but not to recover their associated ground displace-
ment (Fig. 2, B and C). This is because twice-integrated acceleration records
from both consumer and scientific devices are subject to drifts (usually
from tilting) that could render them useless for larger earthquakes after a
small number of seconds. However, Kalman filtering accelerometer data
with consumer C/A code-data yield d
isplacement time series with much
less drift (Fig. 2B).
To assess the performance of a crowdsourced EEW system composed
of consumer devices, we consider two events: an
M
w
7 scenario rupture on
the Hayward fault, and the 2011
M
w
9 Tohoku-Oki earthquake. For the
Hayward fault scenario, we assume that a device is triggered if it and its
four nearest neighbors record displacements greater than 5 cm. If at least
100 devices are triggered, we declare that an event has been detected. In
practice, the availability and noise cha
racteristics of crowdsourced observa-
tions will vary greatly depending on the device and how it is used. For ex-
ample, at any time, a fraction of all
smartphones are turned off, out of
communication, or subject to anthropogenic motion. These devices, how-
ever, produce enough information about data quality, connectivity, and
background motion to permit discrimination of sensors suitable for
EEW. For example, using this system information, it may be desirable to
only include observations from devices that are not otherwise being used by
their owners (and thus not subject to anthropogenic noise), are not
operating on battery power, have a lar
ge-bandwidth telemetry connection,
and whose GNSS position is derived from a sufficient number of satellites.
Therefore, we assume that only a very s
mall subset of potential consumer
devices will be useful during the event. With 0.2% (
n
= 4696) of the pop-
ulation reporting, the earthquake co
uld be detected in 5 s, a sufficient
amount of time to issue a warning to m
ajor population centers (San Fran-
cisco and San Jose) before damaging S waves arrive (Fig. 3). These data are
sufficient to estimate the epicentral location to within 5 km (using a power
law fit to the observed displacements) as soon as the system is triggered, and
to estimate the real-time magnitude evolution with high accuracy and very
little latency (using an analytical finite fault slip inversion) (Fig. 3, D and
̊
0
8
1
E
̊
0
9
̊
0
W
̊
0
9
̊
0
8
1
60 ̊S
0 ̊
60 ̊N
̊
0
8
1
E
̊
0
9
̊
0
W
̊
0
9
̊
0
8
1
60 ̊S
0 ̊
60 ̊N
̊
0
8
1
E
̊
0
9
̊
0
W
̊
0
9
̊
0
8
1
60 ̊S
0 ̊
60 ̊N
̊
0
8
1
E
̊
0
9
̊
0
W
̊
0
9
̊
0
8
1
60 ̊S
0 ̊
60 ̊N
̊
0
8
1
E
̊
0
9
̊
0
W
̊
0
9
̊
0
8
1
60 ̊S
0 ̊
60 ̊N
̊
0
8
1
E
̊
0
9
̊
0
W
̊
0
9
̊
0
8
1
60 ̊S
0 ̊
60 ̊N
̊
0
8
1
E
̊
0
9
̊
0
W
̊
0
9
̊
0
8
1
60 ̊S
0 ̊
60 ̊N
̊
0
8
1
E
̊
0
9
̊
0
W
̊
0
9
̊
0
8
1
60 ̊S
0 ̊
60 ̊N
̊
0
8
1
E
̊
0
9
̊
0
W
̊
0
9
̊
0
8
1
60 ̊S
0 ̊
60 ̊N
012345
PGA (m/s
s
)
̊
0
8
1
E
̊
0
9
̊
0
W
̊
0
9
̊
0
8
1
60 ̊S
0 ̊
60 ̊N
Seismic + GPS EEW
̊
0
8
1
E
̊
0
9
̊
0
W
̊
0
9
̊
0
8
1
60 ̊S
0 ̊
60 ̊N
̊
0
8
1
E
̊
0
9
̊
0
W
̊
0
9
̊
0
8
1
60 ̊S
0 ̊
60 ̊N
Seismic EEW
̊
0
8
1
E
̊
0
9
̊
0
W
̊
0
9
̊
0
8
1
60 ̊S
0 ̊
60 ̊N
̊
0
8
1
E
̊
0
9
̊
0
W
̊
0
9
̊
0
8
1
60 ̊S
0 ̊
60 ̊N
Fig. 1. Global seismic hazard and extent of EEW.
Symbols show the few regions of the world where public citizens and organizations currently
receive earthquake warnings and the types of data used to generate those warnings (
7
). Background color is peak ground acceleration with 10%
probability of exceedance in 50 years from the Global Seismic Hazard Assessment Program.
RESEARCH ARTICLE
Minson
et al
. Sci. Adv. 2015;1:e1500036 10 April 2015
2of7
E) (
7
,
25
). Note that although the earthquake epicenter is a property of
the rupture initiation, moment release
evolves over the duration of the rup-
ture. Thus, the performance of an EEW approach should be evaluated on
how quickly after origin and how accurately the epicenter is obtained, and
how accurately and with how little latency the methodology infers the mo-
ment release of the earthquake as a function of time (
7
,
25
). By both of these
metrics, the crowdsource-derived ep
icenter and moment estimates are ex-
tremely good. Decreasing the partici
pation rate to as little as 0.0125% (
n
=
294) of the population does not degrad
e the quality of the source parameter
estimates, although sparser observa
tion sets require more time to accumu-
late 100 triggers (Fig. 3C). Additionally, it is important to note that these
results hold whether the simulated obs
ervations are from Kalman-filtered
GNSS and INS data (Fig. 3) or from phase-smoothed GNSS data with
SBAS corrections (fig. S4) (
7
).
A
−1 m
1 m
0 s
10 s
20 s
−1 m
1 m
0 s
10 s
20 s
−1 m
1 m
0 s
10 s
20 s
−1 m
1 m
0 s
10 s
20 s
−1 m
1 m
0 s
10 s
20 s
−1 m
1 m
0 s
10 s
20 s
−1 m
1 m
0 s
10 s
20 s
Reference
Accelerometer
GNSS (C/A code)
GNSS (C/A code + p-s + SBAS)
Kalman filter
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
Mw 8.3 Tokachi-oki
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
Mw 9.0 Tohoku-oki
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
Mw 8.8 Maule
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
Mw 6.0 Parkfield
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
Mw 7.2 El Mayor
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
Strong motion seismometer
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
Consumer acc.
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
Preferred KF
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
C/A code
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
C/A code + p-s
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
C/A code + p-s + SBAS
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
Real-time scientific GNSS
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
Smartphone onboard KF
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
B
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
M5
M6
M7
M8
M9
0 km
100 km
200 km
Consumer acc. (acc.)
M5
M6
M7
M8
M9
0 km
100 km
200 km
M5
M6
M7
M8
M9
0 km
100 km
200 km
Strong motion seismometer
M5
M6
M7
M8
M9
0 km
100 km
200 km
M5
M6
M7
M8
M9
0 km
100 km
200 km
Consumer acc. (disp.)
M5
M6
M7
M8
M9
0 km
100 km
200 km
M5
M6
M7
M8
M9
0 km
100 km
200 km
Preferred KF
M5
M6
M7
M8
M9
0 km
100 km
200 km
M5
M6
M7
M8
M9
0 km
100 km
200 km
C/A code
M5
M6
M7
M8
M9
0 km
100 km
200 km
M5
M6
M7
M8
M9
0 km
100 km
200 km
C/A code + p-s
M5
M6
M7
M8
M9
0 km
100 km
200 km
M5
M6
M7
M8
M9
0 km
100 km
200 km
C/A code + p-s + SBAS
M5
M6
M7
M8
M9
0 km
100 km
200 km
M5
M6
M7
M8
M9
0 km
100 km
200 km
Real-time scientific GNSS
M5
M6
M7
M8
M9
0 km
100 km
200 km
M5
M6
M7
M8
M9
0 km
100 km
200 km
Smartphone KF
M5
M6
M7
M8
M9
0 km
100 km
200 km
C
M5
M6
M7
M8
M9
0 km
100 km
200 km
0.001 m
0.01 m
0.1 m
1 m
10 m
s
0
0
1
s
0
5
s
0
Fig. 2. Device tests.
(
A
) Comparison of displacements obtained from
consumer GNSS receivers with and without phase smoothing (p-s) and SBAS,
by twice integrating smartphone acce
leration and by Kalman filtering accel-
eration and GNSS data. Almost any sma
rtphone or similar consumer device
would generate the displacement an
d acceleration dat
a shown with gray
lines. Although these time series individually do a poor job of reproducing
the true time history of motion (shown in black), they can be combined using
a Kalman filter. This process produces
one unified estimate of displacement
(cyan) that is much less noisy than the original acceleration and displacement
data used as inputs to the filter. However, the best GNSS hardware found in
consumer devices (shown in red) is such high quality that there is no need to
supplement the data with acceleration ob
servations, and in fact, the displace-
ment time series could be degraded by doing so. (
B
) Drift of position obtained
from various devices (GNSS, double-int
egrated accelerometers, and Kalman
filtering thereof) compared to obs
erved earthquake displacements (
7
). Neither
GNSS displacement observations nor a
cceleration data double integrated to
displacement are stable over long perio
ds. For GNSS data, this is because the
inherent noise in the observations is not w
hite noise. For acceleration data, this
is because small tilts or steps in the observed acceleration cause large drifts
when integrated. Thus, over time, the apparent position of sensors drifts,
obscuring the true displacement of the instrument. The color curves show
the apparent drift expected for each sensor and data type based on controlled
tests. The black lines are observed displ
acement time series for earthquakes of
different magnitudes. Thus, anywhere t
hat a colored line is below a black line,
the signal-to-noise ratio for that data
type is greater than 1. In a crowdsourced
setting, we expect to obtain data ranging in quality from a Kalman filter of C/A
code data with acceleration data (cya
n line created by combining data from
light blue line with light green line), t
o C/A code data that have been phase-
smoothed (
“
C/A code + p-s,
”
magenta line), to C/A code data that have been
phase-smoothed and supplemented with SBAS (
“
C/A code + p-s + SBAS,
”
red
line). Although all of these data types are significantly noisier than scientific-
quality GNSS data (blue line), they are sensitive enough to record M6-7 earth-
quakes. (
C
) Using the drift curves shown in (B) and the peak ground displacement
expected as function of magnitude
and distance from the source (
27
), we
can calculate the minimum magnitude earthquake observable with a signal-
to-noise ratio of 10. Dotted line shows sensitivity of acceleration recorded
on a smartphone. Dashed lines show sensitivity of displacement data ob-
tained by twice integrating consumer and scientific acceleration data. At
very close distances, the highest-quality consumer devices can observe
earthquakes as small as M6 with a signal-to-noise ratio of at least 10.
RESEARCH ARTICLE
Minson
et al
. Sci. Adv. 2015;1:e1500036 10 April 2015
3of7
We assess the feasibility of crowds
ourced EEW for an actual earthquake
by analyzing C/A code positions for the 2011
M
w
9 Tohoku-Oki event
obtained from data recorded by the GPS Earth Observation Network
(GEONET) (
24
)(Fig.4A).Thedatafromt
hese 462 GEONET stations
represent the scenario of collecting crowdsourced observations from
~0.0004% of the population of Japan. Although the C/A code positioning
technique is identical to that which is used in the consumer devices, it
should be noted that these data were re
corded at scientific GPS stations
equipped with more sophisticated antennas than consumer electronics.
The better hardware decreases multipath noise effects, and so, the C/A
code positions obtained from the GEONET stations should be less noisy
than C/A code time series from consumer electronics. Figure S2, how-
ever, shows that the observed noise level on the GEONET C/A code
time series is equivalent to that obtained from a consumer device when the
C/A code position is supplemented by phase smoothing, and the GEONET
C/A positions are significantly noisier than consumer time series with both
phase smoothing and SBAS. Thus, the Tohoku
time series are, in fact, representative of the data
quality that we might obtain from consumer
devices. In addition to the C/A code position-
ing, we also compute position time series using
scientific-grade processing. The estimated static
displacements from t
hese two processing
methods compare well, following a 1:1
relationship over a range of horizontal displa-
cements from ~0.5 to 4 m, further demonstrat-
ing that consumer-quality GNSS data are
sufficient for EEW (
7
)(fig.S6).Weusesimilar
detection and location criteria and the same
magnitude estimation approach as for the
Hayward fault scenario earthquake. However,
in this case, we examine pre-event C/A code
time series to determine how often a device
might trigger due to noise. We then express
the number of triggers as SDs from the
background triggering rate and use the very
conservative criterion that we will not issue a
warning until the number of triggers exceeds
5
s
of the background triggering rate, which
would, theoretically, restrict the chances of is-
suing a false alarm to about one in 2 million.
With the use of this conservative detection
criterion, the earthquake is detected at 77 s
after the origin time (Fig. 4). At ~100 s, the
solution yields a location of comparable accu-
racy to one obtained by scientific-grade
instruments. The
M
w
estimate follows the
best estimate of the actual moment release
for the earthquake obtained from an inde-
pendent kinematic rupture model (
24
), quick-
ly growing to a final maximum value of ~
M
w
8.8 instead of saturating near
M
w
8likethe
seismic-only EEW estimate. Although not
fast enough to provide a warning for cities
closest to the offshore rupture, this infor-
mation would have permitted a warning
to be issued before the earthquake
’
sdamaging
S waves reached metropolitan Tokyo, and a
tsunami warning could have been issued
minutesbeforethetsunamimadelandfall.
DISCUSSION
Our results demonstrate that the GNSS
and INS navigational sensors built into
A
0 s
30 s
60 s
90 s
0 m
1 m
2 m
E
D
C
123 ̊W
122.5 ̊W
122 ̊W
37.5 ̊N
38 ̊N
38.5 ̊N
123 ̊W
122.5 ̊W
122 ̊W
37.5 ̊N
38 ̊N
38.5 ̊N
123 ̊W
122.5 ̊W
122 ̊W
37.5 ̊N
38 ̊N
38.5 ̊N
123 ̊W
122.5 ̊W
122 ̊W
37.5 ̊N
38 ̊N
38.5 ̊N
123 ̊W
122.5 ̊W
122 ̊W
37.5 ̊N
38 ̊N
38.5 ̊N
Rodgers
C
reek fault
123 ̊W
122.5 ̊W
122 ̊W
37.5 ̊N
38 ̊N
38.5 ̊N
Hayw
ard
fault
123 ̊W
122.5 ̊W
122 ̊W
37.5 ̊N
38 ̊N
38.5 ̊N
San
A
ndreas fault
123 ̊W
122.5 ̊W
122 ̊W
37.5 ̊N
38 ̊N
38.5 ̊N
Pacific Ocean
123 ̊W
122.5 ̊W
122 ̊W
37.5 ̊N
38 ̊N
38.5 ̊N
123 ̊W
122.5 ̊W
122 ̊W
37.5 ̊N
38 ̊N
38.5 ̊N
123 ̊W
122.5 ̊W
122 ̊W
37.5 ̊N
38 ̊N
38.5 ̊N
123 ̊W
122.5 ̊W
122 ̊W
37.5 ̊N
38 ̊N
38.5 ̊N
123 ̊W
122.5 ̊W
122 ̊W
37.5 ̊N
38 ̊N
38.5 ̊N
5
s
123 ̊W
122.5 ̊W
122 ̊W
37.5 ̊N
38 ̊N
38.5 ̊N
SF
SJ
OAK
123 ̊W
122.5 ̊W
122 ̊W
37.5 ̊N
38 ̊N
38.5 ̊N
B
123 ̊W
122.5 ̊W
122 ̊W
37.5 ̊N
38 ̊N
38.5 ̊N
Observation location
Estimated epicenter
Actual hypocenter
0
30
60
90
Time (s)
0
2000
4000
People/km
2
25 km
100
300
500
0 s
2 s
4 s
6 s
8 s
100
300
500
0.0125%
0.025%
0.05%
0.1%
0.2%
Population
Triggers
2 km
4 km
6 km
8 km
0 s
30 s
60 s
90 s
0 s
30 s
60 s
90 s
0.0125%
0.025%
0.05%
0.1%
0.2%
Population
Epicenter mislocation
M6
M7
M6
M7
0.0125%
0.025%
0.05%
0.1%
0.2%
Population
Magnitude
123 ̊W
122.5 ̊W
122 ̊W
37.5 ̊N
38 ̊N
38.5 ̊N
Fig. 3. Hayward fault earthquake scenario.
(
A
) Representative displacement time series from Hayward
fault rupture scenario. Black line: true displacement. Blue line: simulated smartphone C/A code GNSS. Green
line: simulated smartphone accelerometer, twice int
egrated. Red line: Kalman filter combining GNSS and
accelerometer. The red line is representative of the data we expect to observe with the least sophisticated
consumer devices, yet it still does a good job of recovering the true ground motion shown in black. (
B
)
Diamonds showing estimated epicentral location colored by time after origin. As soon as the earthquake
is detected (at 5 s after origin), its epicenter can be e
stimated with an error of less than 5 km using consumer-
quality data. Contour: S wave position when detection criterion is satisfied. Yellow text denotes major cities:
SF, San Francisco; SJ, San Jose; OK, Oakland. Blue dots
denote observer locations assuming 0.2% of the pop-
ulation within the blue box contribute data. (
C
) Number of observers who have detected a potential
earthquake trigger as a function of time. The higher t
he density of observations, the sooner the detection
criteria of a hundred triggers is reached. With just 0.2% of the population contributing data, the earthquake
can be detected in 5 s. (
D
) Epicenter location error as a function of time. The error on the epicenter location is
always <5 km even with very small percentages of
the population contributing observations. (
E
)Estimated
moment magnitude as a function of time for different pa
rticipation levels. Black l
ine: true magnitude. The
estimated magnitude release almost perfectly reprod
uced the actual time history of moment release with
very little latency, even for very low participation rates. The accuracy and low latency of the detection, lo-
cation, and magnitude estimate of the earthquake based on very small numbers of consumer-quality ob-
servations suggest that a crowdsourced EEW system is feasible.
RESEARCH ARTICLE
Minson
et al
. Sci. Adv. 2015;1:e1500036 10 April 2015
4of7
consumer devices such as smartphones
are capable of detecting surface displace-
ments from moderate and larger earth-
quakes. This economical approach
warrants further development, although
we do not suggest that it is a substitute
where monitoring of smaller, but still poten-
tially destructive, earthquakes is required. In
regions where resources cannot be allocated
for scientific-grade EEW due to limited
financial resources or less frequent occur-
rence of destructive earthquakes, crowd-
sourced EEW may be the best option.
For example, large regions of central and
south America, the Caribbean, the Pacific
rim, and south Asia have high seismic
hazard but no early warning capabilities
(Fig. 1). Whether the data from consumer
devices are retrievable from participants
in a crowdsourced monitoring system,
however, depends on each device
’
s
operating system and the levels of access
to raw data permitted by device vendors.
This is a challenge, in fact, germane to
many crowdsourcing endeavors, where
the crowdsource observational objective
may not align with the original commer-
cial intent of the device.
How close are we, then, to operational
crowdsourced EEW? As we have shown,
current smartphones could be used imme-
diately to provide warnings for the largest
earthquakes, such as those associated with
subduction zones, worldwide. Additionally,
because of the enormous number of poten-
tially available devices, a crowdsourced ap-
proach could be very conservative in terms
of data quality control (for example, exclud-
ing devices with poor sky view or subject to
unwanted accelerations such as from auto-
mobile usage) without sacrificing perfor-
mance. Ultimately, to detect a wider
range of earthquakes, access to unfiltered
raw C/A code data through these devices
’
application programming interfaces will
be necessary. This is a trivial technical task
requiring only a software change and no
hardware changes. However, this change
would need to be made in cooperation
with device manufacturers because it
could have commercia
l ramifications.
Because we have shown that inex-
pensive additions, such as adding SBAS
or phase smoothing capability, permit
smartphones to detect moderate to large
earthquakes, an interim solution would
be to deploy these sensors in extremely
low-cost monitoring networks. Although
A
0 s
100 s
200 s
300 s
1 m
3 m
5 m
E
D
C
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
7
7
s
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
1
0
0
s
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
Tokyo
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
B
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
Estimated epicenter
Actual hypocenter
Tsunami arrival
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
1 m offset
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
0
50
100
150
200
250
300
Time (s)
5
15
25
35
45
55
Time (min)
1000
3000
5000
People/km
2
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
25 km
1
σ
3
σ
5
σ
7
σ
0 s
50 s
100 s
1
σ
3
σ
5
σ
7
σ
Triggers
0 km
100 km
200 km
300 km
400 km
0 s
100 s
200 s
300 s
0 s
100 s
200 s
300 s
0 km
100 km
200 km
300 km
400 km
0 km
100 km
200 km
300 km
400 km
Epicenter mislocation
M6
M7
M8
M9
M6
M7
M8
M9
M6
M7
M8
M9
M6
M7
M8
M9
Magnitude
139 ̊E
140 ̊E
141 ̊E
142 ̊E
143 ̊E
144 ̊E
34 ̊N
35 ̊N
36 ̊N
37 ̊N
38 ̊N
39 ̊N
Fig. 4. Tohoku-okiearthquakeexample.
(
A
) Representative displacement time series observed for Tohoku-oki
earthquake. Black line: scientific-grade GPS. Red line: consumer-grade (C/A code) GPS. C/A code GPS
positions are the worst type of data we expect to obtain from consumer devices. However, even these
data do a good job of recovering the actual displacement time series as shown by the scientific-grade GPS
data. (
B
) Diamonds showing estimated epicentral location colored by time after origin. Waves indicate
tsunami arrival times (
28
). Blue contour: S wave position when detection criterion is satisfied. Cyan con-
tour: S wave position when S wave reaches Tokyo. Although there is higher latency in this example than
the Hayward fault example due to the offshore location of the earthquake and the noisier data used, the
proposed crowdsourcing approach could detect and locate the Tohoku earthquake before strong shaking
reaches Tokyo and before the tsunami makes landfall. (
C
) Number of potential earthquake triggers versus
time. We looked at the time series of C/A code positions before the earthquake to determine the frequen-
cy with which a trigger might be observed due to noise. We then expressed the number of triggers as SDs
from that background triggering rate and then, to be very conservative, do not issue a warning until the
number of observed triggers exceeds 5
s
of the background triggering rate. (
D
) Red: location error of our
estimated epicenter relative to the epicenter of (
29
). Purple: error associated with locations reported by
Japan Meteorological Agency (JMA) EEW system. Brown: first location available from global monitoring
(
30
). Although significantly slower than scientific-quality EEW (which includes offshore near-source obser-
vations from ocean-bottom seismometers), the consumer-quality data are capable of determining the
earthquake
’
s location just as accurately as the scientific-quality JMA EEW system and do so significantly
faster than an epicenter could be obtained from global scientific seismic data. (
E
) Red: estimated magni-
tude release as a function of time. Purple:
M
j
values reported by JMA
’
s EEW system. Brown: first
M
w
estimate available from global monitoring (
30
). Black: true magnitude from independent kinematic rup-
ture model (
24
). Again, although there is more latency in the magnitude estimated using only onshore
consumer-quality data than offshore scientific-quality data, the proposed crowdsourced EEW system is
significantly faster than the global response to the earthquake. Also, note that the consumer-quality mag-
nitude, which is based on GNSS data, does not saturate like the seismic magnitudes estimated from
scientific-quality seismic data.
RESEARCH ARTICLE
Minson
et al
. Sci. Adv. 2015;1:e1500036 10 April 2015
5of7
this would require some capital investment for installation and
maintenance, these costs are markedly less than would be required
to build a similar network of scientific instruments. A crowdsourced
EEW system would, of course, have almost zero hardware and main-
tenance cost because the participants would purchase and care for
the sensors (
26
). However, a formal monitoring system built on
consumer sensors would not only cost much less than an equivalent
scientific network but also cost less than the public would spend on a
crowdsourced system because the monitoring network could be built
from the hardware components found within consumer devices. These
components constitute only a small part of the cost of a smartphone or
similar device. Further, this low-cost
network strategy has several attract-
ive aspects in common with scientific-grade networks, including the abil-
ity to concentrate sensors near hazardous faults.
Crowdsourcing is an important phenomenon that has only begun to
be used across the sciences and must be considered seriously. Given the
long repeat times between earthquakes and tsunamis and limited bud-
gets with which to take preventive measures, crowdsourcing may be an
important part of building, maintaining, and operating warning
systems. Crowdsourcing drastically reduces the marginal costs asso-
ciated with EEW because sensor and communication costs would be
assumed by the system
’
s beneficiaries. Further, the commercial push
for ever-greater positioning performance would ensure that a crowd-
sourced EEW network would always incorporate the latest technology
withoutneedforlargeperiodiccapit
al outlays for equipment upgrades.
Finally, by encouraging inclusion of consumer devices into EEW, the
devices can be used not only to gather the observations used to issue
warningsbutalsotodeliverthesewarningstothepublic.Thiswillper-
mit alerts to be customized according to a user
’
s location and should
enhance system efficacy via a feedb
ack process: The more that users en-
gage with the system, the more effective it becomes at reducing the fu-
ture impact from earthquakes and tsunamis.
SUPPLEMENTARY MATERIALS
Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/
full/1/3/e1500036/DC1
Text
Fig. S1. Background position noise for various GNSS receivers found on consumer devices.
Fig. S2. Spectra of drift of position time series for various GNSS receivers found on consumer
devices.
Fig. S3. Observed time series from consumer accelerometers and GNSS receivers.
Fig. S4. Hayward fault earthquake scenario.
Fig. S5. Epicenter location uncertainty for Hayward fault scenario rupture.
Fig. S6. Tohoku-oki earthquake example.
Fig. S7. Epicenter location uncertainty for Tohoku-oki earthquake.
Table S1. Description of observed GPS earthquake displacement time series shown in Fig. 2B.
Table S2. Number of data used and detection response times for Hayward fault simulation.
Table S3. Locations of GEONET GPS stations used in analysis of Tohoku-oki earthquake.
References (
31
–
34
)
REFERENCES AND NOTES
1. R. M. Allen, Seconds before the big one.
Sci. Am.
,75
–
79 (2011).
2. T. Heaton, A model for a seismic computerized alert network.
Science
228
, 987
–
990
(1985).
3. M. Yamada, Warnings work, but must be better.
Nature
473
, 148 (2011).
4. R. M. Allen, P. Gasparini, O. Kamigaichi, M. Böse, The status of earthquake early warning
around the world: An introductory overview.
Seismol. Res. Lett.
80
, 682
–
693 (2009).
5. S. Colombelli, R. M. Allen, A. Zollo, Application of real-time GPS to earthquake early
warning in subduction and strike-slip environments.
J. Geophys. Res.
118
, 3448
–
3461
(2013).
6. D. D. Given, E. S. Cochran, T. Heaton, E. Hauksson, R. Allen, P. Hellweg, J. Vidale, P. Bodin,
Technical Implementation Plan for the ShakeAlert Production System
—
An Earthquake Early
Warning System for the West Coast of the United States. Open-File Report 2014
–
1097
(U.S.
Geological Survey, Reston, VA, 2014), p. 25.
7. See the supplementary materials in
Science Advances
Online for details.
8. A. Smith,
Smartphone Ownership
—
2013 Update
(PewResearchCenter, Washington, DC,
2013).
9. K. Fitchard, Ericsson:
Global Smartphone Penetration Will Reach 60% in 2019
(2013);
http://gigaom.com/2013/11/11/ericsson-global-smartphone-penetration-will-reach-60-in-
2019/.
10. E. Cochran, J. Lawrence, C. Christensen, R. Jakka, The Quake-Catcher Network: Citizen science
expanding seismic horizons.
Seismol. Res. Lett.
80
, 26 (2009).
11. R. W. Clayton, T. Heaton, M. Chandy, A. Krause, M. Kohler, J. Bunn, R. Guy, M. Olson, M. Faulkner,
M. H. Cheng, L. Strand, R. Chandy, D. Obenshain, A. Liu, M. Aivazis, Community Seismic Network.
Ann. Geophys.
54
, 6 (2011).
12. P. S. Earle, D. C. Bowden, M. Guy, Twitter earthquake detection: Earthquake monitoring in a
social world.
Ann. Geophys.
54
, 708
–
715 (2011).
13. T. L. Holzer, J. C. Savage, Global earthquake fatalities and population.
Earthq. Spectra
29
,
155
–
175 (2013).
14. Y. Bock, D. Melgar, B. W. Crowell, Real-time strong-motion broadband displacements from
collocated GPS and accelerometers.
Bull. Seismol. Soc. Am.
101
, 2904
–
2925 (2011).
15. J.F.Lawrence,E.S.Cochran,A.Chung,A.Kaiser,C.M.Christensen,R.Allen,J.W.Baker,B.Fry,
T. Heaton, D. Kilb, M. D. Kohler, M. Taufer, Rapid earthquake characterization using MEMS
accelerometers and volunteer hosts following the
M
7.2 Darfield, New Zealand, earthquake.
Bull.Seismol.Soc.Am.
104
, 184
–
192 (2014).
16. G. Blewitt, C. Kreemer, W. C. Hammond, H.-P. Plag, S. Stein, E. Okal, Rapid determination of
earthquake magnitude using GPS for tsunami warning systems.
Geophys. Res. Lett.
33
(2006).
17. G. Colosimo, M. Crespi, and A. Mazzoni, Real-time GPS seismology with a stand-alone re-
ceiver: A preliminary feasibility demonstration.
J. Geophys. Res.
116
(2011).
18. R. Tu, R. Wang, M. Ge, T. R. Walter, M. Ramatschi, C. Milkereit, D. Bindi, T. Dahm, Cost-
effective monitoring of ground motion related to earthquakes, landslides, or volcanic ac-
tivity by joint use of a single-frequency GPS and a MEMS accelerometer.
Geophys. Res. Lett.
40
, 3825
–
3829 (2013).
19. F. S. T. Van Diggelen,
A-GPS: Assisted GPS, GNSS, and SBAS
(Artech House, Norwood, MA,
2009), pp. 350.
20. R. R. Hatch,
“
The synergism of GPS code and carrier measurements,
”
Proceedings of the
Third International Geodetic Symposium on satellite Doppler Positioning
, Las Cruces, NM, 8 to
12 February 1982, pp. 1213
–
1232.
21. A. Gelb,
Applied Optimal Estimation
(MIT Press, Cambridge, MA, 1974), pp. 142.
22. B. W. Crowell, Y. Bock, M. B. Squibb, Demonstration of earthquake early warning using
total displacement waveforms from real-time GPS networks.
Seismol. Res. Lett.
80
,
772
–
782 (2009).
23. J. Langbein, J. R. Murray, H. A. Snyder, Coseismic and initial postseismic deformation from the
2004 Parkfield, California, earthquake, observed by Global Positioning System, electronic dis-
tance meter, creepmeters, and borehole strainmeters.
Bull. Seismol. Soc. Am.
96
, S304
–
S320
(2006).
24. S. E. Minson, M. Simons, J. L. Beck, F. Ortega, J. Jiang, S. E. Owen, A. W. Moore, A. Inbal, A. Sladen,
Bayesian inversion for finite fault earthquake source models
—
II. The 2011 great Tohoku-oki, Japan
earthquake.
Geophys. J. Int.
198
, 922
–
940 (2014).
25. S. E. Minson, J. R. Murray, J. O. Langbein, J. S. Gomberg, Real-time inversions for finite fault
slip models and rupture geometry based on high-rate GPS data.
J. Geophys. Res.
119
,
3201
–
3231 (2014).
26. J. Rifkin,
The Zero Marginal Cost Society: The Internet of Things, the Collaborative Commons,
and the Eclipse of Capitalism
(Palgrave Macmillan, New York, 2014).
27. K. W. Campbell, Y. Bozorgnia, NGA ground motion model for the geometric mean hori-
zontal component of PGA, PGV, PGD and 5% damped linear elastic response spectra for
periods ranging from 0.01 to 10 s.
Earthq. Spectra
24
, 139
–
171 (2008).
28. K. Satake, Y. Fujii, T. Harada, Y. Namegaya, Time and space distribution of coseismic slip of
the 2011 Tohoku earthquake as inferred from tsunami waveform data.
Bull. Seismol. Soc.
Am
103
, 1473
–
1492 (2013).
29. R. Chu, S. Wei, D. V. Helmberger, Z. Zhan, L. Zhu, H. Kanamori, Initiation of the great M
w
9.0
Tohoku
–
Oki earthquake.
Earth Planet. Sci. Lett.
308
, 277
–
283 (2011).
30. G. P. Hayes, P. S. Earle, H. M. Benz, D. J. Wald, R. W. Briggs, 88 hours: The U.S. Geological
Survey National Earthquake Information Center Response to the 11 March 2011 M
w
9.0
Tohoku earthquake.
Seismol. Res. Lett.
82
, 481
–
493 (2011).
31. J. Langbein, Noise in two-color electronic distance meter measurements revisited.
J. Geophys.
Res.
109
, B04406 (2004).
RESEARCH ARTICLE
Minson
et al
. Sci. Adv. 2015;1:e1500036 10 April 2015
6of7
32. B. T. Aagaard, R. W. Graves, D. P. Schwartz, D. A. Ponce, R. W. Graymer, Ground-motion
modeling of Hayward fault scenario earthquakes, part I: Construction of the suite of scenarios.
Bull. Seism. Soc. Am.
100
, 2927
–
2944 (2010).
33. https://gipsy-oasis.jpl.nasa.gov.
34. M. Simons, S. E. Minson, A. Sladen, F. Ortega, J. Jiang, S. E. Owen, L. Meng, J. P. Ampuero,
S. Wei, R. Chu, D. V. Helmberger, H. Kanamori, E. Hetland, A. W. Moore, F. H. Webb, The
2011 magnitude 9.0 Tohoku-Oki earthquake: Mosaicking the megathrust from seconds
to centuries.
Science
332
, 1421
–
1425 (2011).
Acknowledgments:
We thankB.Atwater,E.Cochran,B.Ellsworth,J.Foster,G.Fryer,R.Stein,C.Wolfe,
and two anonymous reviewers for discussion and reviews.
Funding:
This work was funded in part by
the U.S. Geological Survey Innovation Center for Earth Sciences, U.S. Department of Transportation
Office of the Assistant Secretary for Research and Technology grant RITARS-14-H-HOU awarded
to the University of Houston, and the Gordon and Betty Moore Foundation.
Author contributions:
S.E.M. wrote the manuscript, created the figures, and designed and performed the analyses; B.A.B. led the
project, cowrote the manuscript, performed device tests, and helped design the analyses; C.L.G. per-
formed device tests, helped design the analyses, and wrote the Kalman filter code; J.R.M. performed
device tests and helped design the analyses; J.O.L. performed device tests and the noise analysis; S.E.O.
processed the Tohoku C/A code data; T.H.H. helped design the analyses; R.A.I. helped design the
analyses; and D.L.H. performed device tests.
Competinginterests:
Any use of trade, firm, or product
names is for descriptive purposes only and does not imply endorsement by the U.S. government.
Submitted 10 January 2015
Accepted 6 March 2015
Published 10 April 2015
10.1126/sciadv.1500036
Citation:
S.E.Minson,B.A.Brooks,C.L.Glennie,J.R.Murray,J.O.Langbein,S.E.Owen,
T. H. Heaton, R. A. Iannucci, D. L. Hauser, Crowdsourced earthquake early warning.
Sci. Adv.
1
,
e1500036 (2015).
RESEARCH ARTICLE
Minson
et al
. Sci. Adv. 2015;1:e1500036 10 April 2015
7of7