Search USGS⌕ Search

Geology topics

Diego Melgar

Publications and source records attributed to Diego Melgar.

At least 19 recordsLinked to original sources

Detecting earthquakes in noisy real-time GNSS data with deep learning for improved PGD magnitude estimation

To disseminate accurate and useful warnings, earthquake early warning (EEW) systems must quickly determine the size and location of an earthquake to estimate expected shaking. Traditional seismic‐based algorithms tend to underestimate the true magnitudes of large earthquakes, a phenomenon known as magnitude saturation. This limitation motivated the recent inclusion of Global Navigation Satellite Systems (GNSS) data into the U.S. Geological Survey’s ShakeAlert EEW system with the Geodetic First Approximation of Size and Time (GFAST) algorithm because GNSS data do not saturate with large ground motions. However, the noise levels of GNSS data are very high compared with traditional seismic data, which obscures P ‐wave arrivals and can result in less accurate magnitude estimations if displacement amplitudes are low, such as for lower magnitude earthquakes or large source–station distances. In this study, we develop a deep‐learning model that detects earthquakes in GNSS data and use the Ridgecrest, California, earthquake sequence as a case study to demonstrate how the model could act as a filter to reduce the amount of low‐quality data that enters an algorithm like GFAST. To preserve our limited real earthquake data for model inference, we generated a training dataset composed of >700,000 synthetic displacement waveforms. We combined the synthetic waveforms with real‐time GNSS noise to produce realistically noisy training waveforms and then tested our model on additional synthetic data and performed inference using the real data that were held back. We discuss the performance of our trained model on both the unseen synthetic data and real inference data. Our model can be used to selectively filter only high‐quality data where an earthquake signal is observed for input into an algorithm like GFAST (outperforming a simple signal‐to‐noise ratio–based filter) to reduce the error in GFAST’s real‐time earthquake magnitude estimations.

California↗

The Cascadia Region Earthquake Science Center (CRESCENT): Advancing understanding of Cascadia's earthquake hazards

The Cascadia Subduction Zone (CSZ) is capable of producing large megathrust earthquakes and tsunamis, posing significant risks to the Pacific Northwest. This paper examines the seismic characteristics of the CSZ, including its offshore coupling, sparse seismicity, and a complex paleoseismic record. While recent advances in geophysical methods and data collection have improved our understanding of subduction zone processes, important uncertainties remain regarding rupture dynamics, ground shaking, and tsunami impacts. The Cascadia Region Earthquake Science Center (CRESCENT) has been established to address these challenges through coordinated research, education, and community engagement. By integrating scientific studies, fostering collaboration, and leveraging new technologies, CRESCENT seeks to improve hazard assessment and resilience in the region. This paper outlines the primary challenges and opportunities in advancing subduction zone science and enhancing societal preparedness for future CSZ events.

California, Oregon, Washington↗

Cascadia Subduction Zone science: Call for the next generation community seismic velocity model

The Cascadia subduction zone (CSZ) hosts major seismic and tsunami hazards, yet key questions persist about the relationship between margin structure, fluid distribution, episodic tremor and slip, shallow megathrust behavior, shaking and tsunamigenesis, and the resulting hazard estimates. Addressing these problems requires an empirically grounded, three‐dimensional seismic velocity model to illuminate subsurface structure and properties and to provide a basis for geophysical studies such as earthquake simulations and ground‐motion estimation. In May 2024, the National Science Foundation‐funded Cascadia Region Earthquake Science Center (CRESCENT) community velocity model (CVM) working group, with U.S. Geological Survey and regional partners, convened a workshop to identify priorities for such a model. Participants emphasized the features necessary for addressing key science questions, including implementing findability, accessibility, interoperability, and reusability (FAIR) access, capturing along‐strike and along‐dip structural heterogeneity, resolving shallow offshore–onshore structure, constraining elastic properties and quantifying their uncertainties for numerical wave propagation simulations, their validation benchmarks, and supporting associated accurate earthquake ground‐motion simulations and hazard assessments. This article describes the priorities defined in the workshop, and a description of how, guided by these needs, CRESCENT plans to develop multiple generations of a CVM to advance CSZ science and improve seismic and tsunami hazard modeling across the Pacific Northwest. The CVM will span the CSZ from the surface to ∼100 km depth, offshore and east of the Cascades into Idaho (∼132°–110° W) and the southern and northern tectonic regime transitions (∼36°–52° N) to capture the entire tectonic system as well as its surroundings.

Cascadia Subduction Zone↗

Wavelet Inversion for SliP (WISP): Open-source earthquake slip modeling software

Models of the spatiotemporal evolution of earthquake slip, termed finite-fault models, are a critical component of rapid earthquake and tsunami response, earthquake forecasting, seismic ground-motion estimates, and studies of earthquake kinematics. Here, we detail a newly released finite-fault modeling software, Wavelet Inversion for SliP (WISP), in use at the U.S. Geological Survey’s National Earthquake Information Center (NEIC) and available to the public. WISP version 1.1.0 allows inversion of teleseismic body and surface waves, as well as local strong-motion, static and dynamic Global Navigation Satellite System, and satellite imagery (e.g., Interferometric Synthetic Aperture Radar) observations on single or multiple planar fault segments. The software is used in NEIC rapid response of earthquakes M w ≥ 7, generally resulting in a published model within the first few hours after the event origin time. The rupture location and dimensions are then used as inputs to downstream products to estimate earthquake shaking, predict loss, and model the likelihood of secondary hazards, namely landslides and liquefaction. WISP is also used in research studies to evaluate the characteristics of complex ruptures including multifault ruptures and earthquake doublets, among others. The WISP version 1.1.0 software release is composed of Python-wrapped FORTRAN code to accomplish the inversion procedure. A simple command line interface facilitates ease of use even for those with only a cursory knowledge of Python scripting. WISP version 1.1.0 includes a Jupyter Notebook tutorial demonstrating use of the software for modeling the 2015 M w 8.3 Illapel, Chile, earthquake. In parallel with the tutorial, we demonstrate the typical usage of the WISP software using the M w 8.3 Illapel earthquake example here.

Seismological Research Letters↗

Slow rupture, long rise times, and multi-fault geometry: The 2020 M6.4 southwestern Puerto Rico mainshock

The M 6.4 mainshock of the southwestern Puerto Rico seismic sequence on 7 January 2020, was one of the most impactful modern earthquakes in the northeastern Caribbean. Due to its offshore location and complex aftershock distribution, its source kinematics remain poorly constrained. This active sequence illuminated a complex set of previously unrecognized structures that indicate multiple causative faults may have slipped during its rupture. Here, we utilize seismic and geodetic observations to enhance model resolution, estimate the finite slip of the mainshock, and test a multi-segment, geologically realistic fault geometry. Our refined model finds a lower rupture velocity and longer rise times than typical for an event of this magnitude. This indicates a slow-evolving rupture process that resembles characteristics of a tsunami earthquake. Although this normal/strike-slip faulting event was not tsunamigenic, these qualities, if pervasive for this region, have important implications for future seismic monitoring and hazards in southwestern Puerto Rico.

Puerto Rico↗

Increased flood exposure in the Pacific Northwest following earthquake-driven subsidence and sea-level rise

Climate-driven sea-level rise is increasing the frequency of coastal flooding worldwide, exacerbated locally by factors like land subsidence from groundwater and resource extraction. However, a process rarely considered in future sea-level rise scenarios is sudden (over minutes) land subsidence associated with great (>M8) earthquakes, which can exceed 1 m. Along the Washington, Oregon, and northern California coasts, the next great Cascadia subduction zone earthquake could cause up to 2 m of sudden coastal subsidence, dramatically raising sea level, expanding floodplains, and increasing the flood risk to local communities. Here, we quantify the potential expansion of the 1% floodplain (i.e., the area with an annual flood risk of 1%) under low (~0.5 m), medium (~1 m), and high (~2 m) earthquake-driven subsidence scenarios at 24 Cascadia estuaries. If a great earthquake occurred today, floodplains could expand by 90 km 2 (low), 160 km 2 (medium), or 300 km 2 (high subsidence), more than doubling the flooding exposure of residents, structures, and roads under the high subsidence scenario. By 2100, when climate-driven sea-level rise will compound the hazard, a great earthquake could expand floodplains by 170 km 2 (low), 240 km 2 (medium), or 370 km 2 (high subsidence), more than tripling the flooding exposure of residents, structures, and roads under the high subsidence scenario compared to the 2023 floodplain. Our findings can support decision-makers and coastal communities along the Cascadia subduction zone as they prepare for compound hazards from the earthquake cycle and climate-driven sea-level rise and provide critical insights for tectonically active coastlines globally.

California, Oregon, Washington↗

An unexplained tsunami: Was there megathrust slip during the 2020 Mw7.6 Sand Point, Alaska, earthquake?

On October 19, 2020, the M w 7.6 Sand Point earthquake struck south of the Shumagin Islands in Alaska. Moment tensors indicate the earthquake was primarily strike-slip, yet the event produced an enigmatic tsunami that was larger and more widespread than expected for an earthquake of that magnitude and mechanism. Using a suite of hydrodynamic, seismic, and geodetic modeling techniques, we explore plausible causes of the tsunami. We find that strike-slip models consistent with the moment tensor orientation cannot produce the observed tsunami. Hydrodynamic inversion of sea surface deformation from deep ocean and tide gauge data suggest seafloor deformation more closely matches a megathrust, rather than a strike-slip, source. Static slip inversions, using sea level and Global Navigation Satellite System data, allow for a portion of co-seismic megathrust slip that can explain tsunamigenesis. Combining all available geophysical datasets to model the kinematic rupture, we show that considerable, relatively slow, megathrust slip is allowable in the Shumagin segment, concurrent with strike-slip faulting. We hypothesize that the slow megathrust rupture does not contribute much seismic radiation allowing it to previously go unnoticed with traditional seismic monitoring.

Alaska↗

Combining multisite tsunami and deformation modeling to constrain slip distributions for the 1700 C.E. Cascadia earthquake

A major earthquake ruptured the Cascadia subduction zone (CSZ) on 26 January 1700. Key paleoseismic evidence associated with this event include tsunami deposits, stratigraphic evidence of coastal coseismic subsidence, written Japanese records of a tsunami unaccompanied by earthquake shaking, and margin‐wide turbidites found offshore and in lacustrine environments. Despite this wealth of independent clues, important details about this event remain unresolved. Dating uncertainties do not conclusively establish whether the proxies are from one earthquake or a sequence of them, and we have limited knowledge of the likely slip distributions of the event or events. Here, we use a catalog of 37,500 candidate synthetic ruptures between M w 7.8 and 9.2 and simulate their resulting coseismic deformation and tsunami inundation. Each model is then compared against estimated Japan tsunami arrivals, regional coastal subsidence records, and local paleotsunami deposits mapped at six different coastal marshes and one coastal lake along the CSZ. We find that seven full‐margin ruptures with a median magnitude of M w 9.1 satisfy all three constraints. We favor one M w 9.11 model that best matches all site paleoseismic observations and suggests that the Cascadia megathrust slipped up to ∼30 m and must have shallow geodetic coupling. We also find that some sequences composed of three or four ruptures can still satisfy the observations, yet no sequences of two ruptures can. Sequences are differentiated into three groups based on whether they contain a mainshock rupture located in the south (>44° N) or further north. All sequences contain unruptured portions of the megathrust and most contain mainshocks with peak slip above 40 m. The fit of the geologic evidence from sequences is poor in comparison to single‐event models. Therefore, sequences are generally less favored compared to full‐margin events.

British Columbia, California, Oregon, Washington↗

Evaluation of an open earthquake early warning system in Mexico, and laboratory tests of their sensors

Microelectromechanical system (MEMS) accelerometers are useful for seismological and engineering applications because of their ability to record unsaturated large seismic signals. Recent advances in MEMS technologies enable the design of instruments with improved capabilities that also allow the recording of small signals. As a result, MEMS can be useful across a broad dynamic range and for both major earthquakes and smaller magnitude events. Here, we analyze the performance of a MEMS‐based accelerometer with the capability of real‐time data transmission. We discuss the self‐noise level, dynamic range, and sensitivity, along with the timing precision and data transmission latency, and compare these parameters to other MEMS‐based instruments. We also summarize observations from a regional network deployed in southern Mexico over a period of 3+ yr for the purpose of earthquake early warning. In addition to the sensor evaluation, we present a large, openly available data set of strong motion data from the Mexican network that comprises continuous ground‐motion records from 24 accelerometers in the period between 2017 and 2022.

Seismological Research Letters↗

Rapid estimation of single-station earthquake magnitudes with machine learning on a global scale

The foundation of earthquake monitoring is the ability to rapidly detect, locate, and estimate the size of seismic sources. Earthquake magnitudes are particularly difficult to rapidly characterize because magnitude types are only applicable to specific magnitude ranges, and location errors propagate to substantial magnitude errors. We developed a method for rapid estimation of single‐station earthquake magnitudes using raw three‐component P waveforms observed at local to teleseismic distances, independent of prior size or location information. We used the MagNet regression model architecture ( Mousavi and Beroza, 2020b ), which combines convolutional and recurrent neural networks. We trained our model using ∼2.4 million P ‐phase arrivals labeled by the authoritative magnitude assigned by the U.S. Geological Survey. We tested input data parameters (e.g., window length) that could affect the performance of our model in near‐real‐time monitoring applications. At the longest waveform window length of 114 s, our model (Artificial Intelligence Magnitude [AIMag]) is accurate (median estimated magnitude within ±0.5 magnitude units from catalog magnitude) between M 2.3 and 7.6. However, magnitudes above M ∼7 are more underestimated as true magnitude increases. As the windows are shortened down to 1 s, the point at which higher magnitudes begin to be underestimated moves toward lower magnitudes, and the degree of underestimation increases. The over and underestimation of magnitudes for the smallest and largest earthquakes, respectively, are potentially related to the limited number of events in these ranges within the training data, as well as magnitude saturation effects related to not capturing the full source time function of large earthquakes. Importantly, AIMag can determine earthquake magnitudes with individual stations’ waveforms without instrument response correction or knowledge of an earthquake’s source‐station distance. This work may enable monitoring agencies to more rapidly recognize large, potentially tsunamigenic global earthquakes from few stations, allowing for faster event processing and reporting. This is critical for timely warnings for seismic‐related hazards.

Bulletin of the Seismological Society of America↗

Beyond the teleseism: Introducing regional seismic and geodetic data into routine USGS finite‐fault modeling

The U.S. Geological Survey (USGS) National Earthquake Information Center (NEIC) routinely produces finite‐fault models following significant earthquakes. These models are spatiotemporal estimates of coseismic slip critical to constraining downstream response products such as ShakeMap ground motion estimates, Prompt Assessment of Global Earthquake for Response loss estimates, and ground failure assessments. Because large earthquakes can involve slip over tens to hundreds of kilometers, point‐source approximations are insufficient, and it is vital to rapidly assess the amount, timing, and location of slip along the fault. Initially, the USGS finite‐fault products were computed in the first several hours after a significant earthquake, using teleseismic body wave and surface wave observations. With only teleseismic waveforms, it is generally possible to obtain a reliable model for earthquakes of magnitude 7 and larger. Here, we detail newly implemented updates to NEIC’s modeling capabilities, specifically to allow joint modeling of local‐to‐regional strong‐motion accelerometer, Global Navigation Satellite System (GNSS), and Interferometric Synthetic Aperture Radar (InSAR) observations in addition to teleseismic waveforms. We present joint inversion results for the 2015 M w ">

Seismological Research Letters↗

Zoning Verification in Mexico City Using Strong Motions of the M7.1 Puebla-Morelos Earthquake of September 19, 2017

Mexico City suffers extensive damage from large earthquakes that originate at far distances due mainly to densely built areas on a filled lakebed. Seismic design codes in Mexico recognize the site-specific issues in Mexico City by acknowledging zones that represent lakebed as being riskier compared to other Mexico City areas. At the time of the 19 September 1985 M8.1 Michoacán Earthquake, (a) the seismic zoning map comprised only three zones (Hills [now Zone I], Transition [now Zone II] and Lake Zone [now Zone III] and (b) a limited number of permanent strong motion stations (~6) existed. Since, 1985 the seismic zonation maps have evolved: (a) the Lake Zone (now Zone III) into 4 sub-zones (a, b, c and d) and (b) the number of strong motion stations has increased significantly (~25) and recorded the 19 September 2017 M7.1 Puebla-Morelos earthquake. To identify the site frequency [f] (period [T]) at the specific location of a strong motion station within a particular zone or sub-zone (Zone I, II and III a,b,c,d), two well-known methods (H/H and H/V) for determining spectral ratios are employed using these records. The predominant site periods identified by both methods are in good agreement with those interpolated from the zoning-site period maps of Mexico City in the 2004 Seismic Design Code. Such verification is a testament to the benefits of improved zoning and greater distribution of strong motion network stations, and sets a good example for applications in other seismic regions.

Mexico City↗

A ground motion model for GNSS peak ground displacement

We present an updated ground‐motion model (GMM) for M w "> M w Mw 6–9 earthquakes using Global Navigation Satellite Systems (GNSS) observations of the peak ground displacement (PGD). Earthquake GMMs inform a range of Earth science and engineering applications, including source characterization, seismic hazard evaluations, loss estimates, and seismic design standards. A typical GMM is characterized by simplified metrics describing the earthquake source (magnitude), observation distance, and site terms. Most often, GMMs are derived from broadband seismometer and accelerometer observations, yet during strong shaking, these traditional seismic instruments are affected by baseline offsets, leading to inaccurate recordings of low‐frequency ground motions such as displacement. The incorporation of geodetic data sources, particularly for characterizing the unsaturated ground displacement of large‐magnitude events, has proven valuable as a complement to traditional seismic approaches and led to the development of an initial point‐source GMM based on PGD estimated from high‐rate GNSS data. Here, we improve the existing GMM to more effectively account for fault finiteness, slip heterogeneity, and observation distance. We evaluate the limitations of the currently available GNSS earthquake data set to calibrate the GMM. In particular, the observed earthquake data set is lacking in observations within 100 km of large‐magnitude events ( ⁠ M w > 8 "> M w > 8 Mw>8 ⁠ ), inhibiting evaluation of fault dimensions for earthquakes too large to be represented as point sources in the near field. To that end, we separately consider previously validated synthetic GNSS waveforms within 10–1000 km of M w "> M w Mw 7.8–9.3 Cascadia subduction zone scenario ruptures. The synthetic data highlight the importance of fault distance rather than point‐source metrics and improve our preparedness for large‐magnitude earthquakes with spatiotemporal qualities unlike those in our existing data set.

Bulletin of the Seismological Society of America↗

Characterizing large earthquakes before rupture is complete

Whether large and very large earthquakes are distinguishable from each other early on in the rupture process has been a subject often debated over the past several decades. Studies have shown that the frequency content of radiated seismic energy in the first few seconds of an earthquake scales with the final magnitude of the event, implying determinism. Other studies have shown that the recordings of ground displacement from small-to-moderate sized earthquakes are indistinguishable, and thus earthquakes share a universal early rupture process. Regardless of how earthquakes start, however, at some point in the rupture process events of different sizes must be distinguishable from one another. If that difference occurs early - i.e., before the rupture duration of the smaller event - this implies that earthquakes demonstrate some level of determinism. Here we show through analysis of a large database of source time functions and near-source displacement records that after an initiation phase, ruptures of M7-9 earthquakes organize into a slip pulse, the kinematic properties of which scale with magnitude. As such, early in the rupture process - after about 10s - large and very large earthquakes demonstrate different properties and can thus be distinguished.

Science Advances↗

Ground motions from the 7 and 19 September, 2017 Tehuantepec and Puebla-Morelos, Mexico, earthquakes

The 2017 M "> M 8.2 Tehuantepec and M "> M 7.1 Puebla‐Morelos earthquakes were deep inslab normal‐faulting events that caused significant damage to several central‐to‐southern regions of Mexico. Inslab earthquakes are an important component of seismicity and seismic hazard in Mexico. Ground‐motion prediction equations (GMPEs) are an integral part of seismic hazard assessment as well as risk and rapid‐response products. This work examines the observed ground motions from these two events in comparison to the predicted median ground motions from four GMPEs. The residuals between the observed and modeled ground motions allow us to study regional differences in shaking, the effects of each earthquake, and basin effects in Mexico City, Puebla, and Oaxaca. We find that the ground motions from these two earthquakes are generally well modeled by the GMPEs. However, the Tehuantepec event shows larger than expected ground motions at greater distances and longer periods, which suggests a waveguide effect from the subduction zone geometry. Finally, Mexico City and the cities of Puebla and Oaxaca exhibit very large ground motions, indicative of well‐known site and basin effects that are much stronger than the basin terms included in some of the GMPEs. Simple and rapid ground‐motion parameter estimates that include site effects are key for hazard and real‐time risk assessments in regions such as Mexico, where the vast majority of the population lives in areas where the aforementioned effects are relevant. However, GMPEs based on site correction terms dependent on topographic slope proxies underestimate, at least in the three cities tackled in this work, the observed amplification. Therefore, there is a need to improve models of seismic amplification in basins that could be included in GMPEs.

Bulletin of the Seismological Society of America↗

Development of a geodetic component for the U.S. West Coast Earthquake Early Warning System

An earthquake early warning (EEW) system, ShakeAlert, is under development for the West Coast of the United States. This system currently uses the first few seconds of waveforms recorded by seismic instrumentation to rapidly characterize earthquake magnitude, location, and origin time; ShakeAlert recently added a seismic line source algorithm. For large to great earthquakes, magnitudes estimated from the earliest seismic data alone generally saturate. Real‐time Global Navigation Satellite System (GNSS) data can directly measure large displacements, enabling accurate magnitude estimates for M w 7 + "> M w 7 + events, possibly before rupture termination. GNSS‐measured displacements also track evolving slip and, alone or in combination with seismic data, constrain finite‐fault models. Particularly for large‐magnitude, long‐rupture events, GNSS‐based magnitude and rupture extent estimates can improve updates to predicted shaking and thus alert accuracy. GNSS data processing centers at ShakeAlert partner institutions provide real‐time streams to the EEW system, and three geodetic EEW algorithms have been developed through the ShakeAlert collaboration. These algorithms will undergo initial testing within ShakeAlert’s computational architecture using a suite of input data that includes simulated real‐time displacements from synthetic earthquakes and GNSS recordings from recent earthquakes worldwide. Performance will be evaluated using metrics and standards consistent with those adopted for ShakeAlert overall. This initial assessment will guide method refinement and synthesis of the most successful features into a candidate geodetic algorithm for the ShakeAlert production system. In parallel, improvements to geodetic networks and streamlining approaches to data processing and exchange will ensure robust geodetic data availability in the event of an earthquake.

West coast↗

The 19 September 2017 M 7.1 Puebla‐Morelos earthquake: Spectral ratios confirm Mexico City zoning

One important element of understanding basin response to strong shaking is the analysis of spectral ratios, which may provide information about the dominant frequency of ground motion at specific locations. Spectral ratios computed from accelerations recorded by strong‐motion stations in Mexico City during the mainshock of the 19 September 2017 M "> M M 7.1 Puebla‐Morelos earthquake reveal predominate periods consistent with those mapped in the 2004 Mexican seismic design code. Furthermore, the predominant periods thus computed validate those studies using mainshock and aftershock recordings of the handful strong‐motion stations that recorded the 19 September 1985 M "> M M 8.1 Michoacán earthquake. Even though the number of stations in each of the zones (zones I, II, IIIa, b, c, and d) is not the same, they still allow confirmation of site frequencies (periods) attributable to the specific zones (particularly those in zones IIIa, b, c, and d). Spectral ratios are computed with two different methods: (1) horizontal to horizontal (H/H) ratio of smoothed amplitude spectrum of a horizontal channel in direction X of a station with respect to the smoothed amplitude spectrum of the horizontal channel in the same X direction of a reference stiff soil (or rock) station and (2) horizontal to vertical (H/V) ratio (or also known as the Nakamura method) of both horizontal (H) and vertical (V) channels of the same station. We show a comparison of the identified frequencies (periods) derived by both methods and find they are very similar and in good agreement with those indicated in the zoning maps of Mexico City in the 2004 seismic design code.

Mexico City↗

Research to improve ShakeAlert earthquake early warning products and their utility

Earthquake early warning (EEW) is the rapid detection of an earthquake and issuance of an alert or notification to people and vulnerable systems likely to experience potentially damaging ground shaking. The level of ground shaking that is considered damaging is defined by the specific application; for example, manufacturing equipment may experience damage at a lower intensity ground shaking than would cause damage to a building. Along the West Coast of the United States, the warning times for ground shaking could range as high as tens of seconds for moderate levels of ground shaking, or potentially longer, if a lower ground-shaking threshold is used to issue alerts. However, it is not always possible to provide advance warning of ground shaking, particularly for locations close to an earthquake that are most likely to experience very strong ground shaking. EEW alerts may be useful to individuals who can use a few seconds to move to a safe zone and to electromechanical systems that can take automatic actions to reduce damage and injuries. An EEW system, ShakeAlert, has been under development in the United States since 2006. Federal and State governments, as well as the private sector, are now investing in the ShakeAlert prototype system that will, when completed, become an operational public system for the West Coast of the United States. While the current prototype is delivering alerts to test users, improvements to the accuracy, timeliness, and utility of the alerts are needed. For this reason, it is essential that the ShakeAlert system be continuously improved through targeted research, involving not only the current ShakeAlert partner organizations, but also the broader scientific, engineering, and emergencyresponse communities. To this end, this report describes the opportunities for improvement that can be addressed through research and development over the next 5 years. Our recommendations are organized into four areas: (1) understand EEW capabilities and user needs, (2) make alerts as fast and accurate as possible, (3) ensure reliability when it counts, and (4) explore the use of new instrumentation. The first challenge is to understand EEW capabilities and user needs. EEW must deliver actionable information to people and to automated systems to mitigate short- and longterm impacts of damaging ground shaking, so development of EEW must be motivated by the needs of users. Within this challenge, we must study the technical capabilities and limitations of EEW in general, and the ShakeAlert system specifically. This includes development of performance metrics that assess the timeliness and accuracy of alerts to understand the value and utility of the ShakeAlert EEW product(s) for various user groups, including different industry sectors, emergency-management agencies, and the public. Research is needed to define the alerting choices that maximize the utility of the system for users and to determine what the available communication pathways are for providing timely alert information. Additionally, we engage users to assess how alerts will be used by different sectors to mitigate losses and to inform EEW product design. Further, social-science research is needed to develop alert messaging, including what relevant prior and follow-up information are required, to ensure effective use of alerts. The second challenge is to make alerts as fast and as accurate as possible. The timeliness and accuracy of an EEW alert is important because it will set in motion a series of actions and downstream products. An EEW alert will trigger notification across emergency-alert systems and across multiple communication channels to populations in impacted regions. The EEW alert region may grow as the earthquake fault-rupture length increases, and the EEW system’s characterization of it, evolves. We must continue research into new or improved seismic and geodetic waveform-processing methods necessary to rapidly characterize the expected ground shaking and associated uncertainties. It is important to thoroughly evaluate whether new methods improve alerts through more accurate ground-motion estimates and (or) reduced latencies (that is, longer warning times). New methods could include tracking the extent of a large rupture in real time (known as finite-fault algorithms) and ground-motionbased EEW algorithms. Additionally, ground motion predictions could be optimized for each earthquake as the earthquake fault rupture progresses by using, for example, event terms to shift ground-motion curves for more (or less) energetic ruptures. The third challenge is to ensure reliability when it counts. This challenge requires us to explore approaches that assess the expected performance of ShakeAlert across the range of earthquake magnitudes, locations, and depths that may occur within the alerting region. Large, damaging earthquakes and their associated aftershock sequences matter most for hazard and for EEW, but these large-earthquake sequences occur infrequently. We expect ShakeAlert to respond robustly to these large-earthquake sequences despite potentially long periods of relative seismic quiescence in the intervening years, and in spite of inevitable communication challenges that arise during and after a large earthquake. We must develop methods to utilize the broadest available datasets to test EEW performance, including ground-motion data recorded in other parts of the world. The observational period for large, damaging earthquakes in any particular region has been short in comparison to estimated large-earthquake recurrence times. Ground-motion records for very large, damaging western United States events and major aftershock sequences do not yet exist, nor do data exist for all potential sources of noise and spurious signals that ShakeAlert must be “tuned” to reject. In addition, robust synthetic data could provide the flexibility to test a wider range of earthquake magnitude, tectonic-setting, and noise scenarios than are covered by existing observational data. Synthetic ground-motion data must be thoroughly vetted against records of smaller magnitude earthquakes to ensure that they accurately capture both the onset and the amplitude of the ground shaking. The final challenge is to explore the use of new instrumentation. The development of EEW around the world to date has focused on the use of high-quality, scientific-grade seismic and geodetic instrumentation. The use of additional types of instrumentation or information may also improve EEW products by filling gaps in sensor coverage in countries that already have dense seismic networks or enable EEW in countries without such networks. We must keep up with these developments and continuously assess their value in supplementing existing EEW systems, such as ShakeAlert, or enabling EEW where such systems do not exist. Such developments include low-cost instrumentation with microelectromechanical system (MEMS) sensors and global positioning system (GPS)/global navigation satellite system (GNSS) antennas embedded in low-cost consumer electronics, sea-floor seismometers, geodetic instrumentation deployed along the Cascadia and Alaska megathrust margins of western North America, and borehole strainmeters that are already deployed across the region.

Open-File Report↗