Geology topics
Publications and source records attributed to David J. Wald.
ShakeCast is a freely available, post-earthquake situational awareness application that automatically retrieves earthquake shaking data from ShakeMap, compares intensity measures against users? facilities, and generates potential damage assessment notifications, facility damage maps, and other Web-based products for emergency managers and responders.
Earthquakes have claimed approximately 8 million lives over the last 2,000 years (Dunbar, Lockridge and others, 1992) and fatality rates are likely to continue to rise with increased population and urbanizations of global settlements especially in developing countries. More than 75% of earthquake-related human casualties are caused by the collapse of buildings or structures (Coburn and Spence, 2002). It is disheartening to note that large fractions of the world's population still reside in informal, poorly-constructed & non-engineered dwellings which have high susceptibility to collapse during earthquakes. Moreover, with increasing urbanization half of world's population now lives in urban areas (United Nations, 2001), and half of these urban centers are located in earthquake-prone regions (Bilham, 2004). The poor performance of most building stocks during earthquakes remains a primary societal concern. However, despite this dark history and bleaker future trends, there are no comprehensive global building inventories of sufficient quality and coverage to adequately address and characterize future earthquake losses. Such an inventory is vital both for earthquake loss mitigation and for earthquake disaster response purposes. While the latter purpose is the motivation of this work, we hope that the global building inventory database described herein will find widespread use for other mitigation efforts as well. For a real-time earthquake impact alert system, such as U.S. Geological Survey's (USGS) Prompt Assessment of Global Earthquakes for Response (PAGER), (Wald, Earle and others, 2006), we seek to rapidly evaluate potential casualties associated with earthquake ground shaking for any region of the world. The casualty estimation is based primarily on (1) rapid estimation of the ground shaking hazard, (2) aggregating the population exposure within different building types, and (3) estimating the casualties from the collapse of vulnerable buildings. Thus, the contribution of building stock, its relative vulnerability, and distribution are vital components for determining the extent of casualties during an earthquake. It is evident from large deadly historical earthquakes that the distribution of vulnerable structures and their occupancy level during an earthquake control the severity of human losses. For example, though the number of strong earthquakes in California is comparable to that of Iran, the total earthquake-related casualties in California during the last 100 years are dramatically lower than the casualties from several individual Iranian earthquakes. The relatively low casualties count in California is attributed mainly to the fact that more than 90 percent of the building stock in California is made of wood and is designed to withstand moderate to large earthquakes (Kircher, Seligson and others, 2006). In contrast, the 80 percent adobe and or non-engineered masonry building stock with poor lateral load resisting systems in Iran succumbs even for moderate levels of ground shaking. Consequently, the heavy death toll for the 2003 Bam, Iran earthquake, which claimed 31,828 lives (Ghafory-Ashtiany and Mousavi, 2005), is directly attributable to such poorly resistant construction, and future events will produce comparable losses unless practices change. Similarly, multistory, precast-concrete framed buildings caused heavy casualties in the 1988 Spitak, Armenia earthquake (Bertero, 1989); weaker masonry and reinforced-concrete framed construction designed for gravity loads with soft first stories dominated losses in the Bhuj, India earthquake of 2001 (Madabhushi and Haigh, 2005); and adobe and weak masonry dwellings in Peru controlled the death toll in the Peru earthquake of 2007 (Taucer, J. and others, 2007). Spence (2007) after conducting a brief survey of most lethal earthquakes since 1960 found that building collapses remains a major cause of earthquake mortality and unreinforced masonry buildings are one of the mos
An atlas of maps of peak ground motions and intensity 'ShakeMaps' has been developed for almost 5,000 recent and historical global earthquakes. These maps are produced using established ShakeMap methodology (Wald and others, 1999c; Wald and others, 2005) and constraints from macroseismic intensity data, instrumental ground motions, regional topographically-based site amplifications, and published earthquake-rupture models. Applying the ShakeMap methodology allows a consistent approach to combine point observations with ground-motion predictions to produce descriptions of peak ground motions and intensity for each event. We also calculate an estimated ground-motion uncertainty grid for each earthquake. The Atlas of ShakeMaps provides a consistent and quantitative description of the distribution and intensity of shaking for recent global earthquakes (1973-2007) as well as selected historic events. As such, the Atlas was developed specifically for calibrating global earthquake loss estimation methodologies to be used in the U.S. Geological Survey Prompt Assessment of Global Earthquakes for Response (PAGER) Project. PAGER will employ these loss models to rapidly estimate the impact of global earthquakes as part of the USGS National Earthquake Information Center's earthquake-response protocol. The development of the Atlas of ShakeMaps has also led to several key improvements to the Global ShakeMap system. The key upgrades include: addition of uncertainties in the ground motion mapping, introduction of modern ground-motion prediction equations, improved estimates of global seismic-site conditions (VS30), and improved definition of stable continental region polygons. Finally, we have merged all of the ShakeMaps in the Atlas to provide a global perspective of earthquake ground shaking for the past 35 years, allowing comparison with probabilistic hazard maps. The online Atlas and supporting databases can be found at http://earthquake.usgs.gov/eqcenter/shakemap/atlas.php/.
We describe algorithms for quantifying and qualifying uncertainties associated with USGS ShakeMap ground motions. The uncertainty values computed consist of latitude/longitude grid-based multiplicative factors that scale the standard deviation associated with the ground motion prediction equation (GMPE) used within the ShakeMap algorithm for estimating ground motions. The resulting grid-based 'uncertainty map' is essential for evaluation of losses derived using ShakeMaps as the hazard input. For ShakeMap, ground motion uncertainty at any point is dominated by two main factors: (i) the influence of any proximal ground motion observations, and (ii) the uncertainty of estimating ground motions from the GMPE, most notably, elevated uncertainty due to initial, unconstrained source rupture geometry. The uncertainty is highest for larger magnitude earthquakes when source finiteness is not yet constrained and, hence, the distance to rupture is also uncertain. In addition to a spatially-dependant, quantitative assessment, many users may prefer a simple, qualitative grading for the entire ShakeMap. We developed a grading scale that allows one to quickly gauge the appropriate level of confidence when using rapidly produced ShakeMaps as part of the post-earthquake decision-making process or for qualitative assessments of archived or historical earthquake ShakeMaps. We describe an uncertainty letter grading ('A' through 'F', for high to poor quality, respectively) based on the uncertainty map. A middle-range ('C') grade corresponds to a ShakeMap for a moderate-magnitude earthquake suitably represented with a point-source location. Lower grades 'D' and 'F' are assigned for larger events (M>6) where finite-source dimensions are not yet constrained. The addition of ground motion observations (or observed macroseismic intensities) reduces uncertainties over data-constrained portions of the map. Higher grades ('A' and 'B') correspond to ShakeMaps with constrained fault dimensions and numerous stations, depending on the density of station/data coverage. Due to these dependencies, the letter grade can change with subsequent ShakeMap revisions if more data are added or when finite-faulting dimensions are added. We emphasize that the greatest uncertainties are associated with unconstrained source dimensions for large earthquakes where the distance term in the GMPE is most uncertain; this uncertainty thus scales with magnitude (and consequently rupture dimension). Since this distance uncertainty produces potentially large uncertainties in ShakeMap ground-motion estimates, this factor dominates over compensating constraints for all but the most dense station distributions.
At the beginning of 2006, the U.S. Geological Survey (USGS) Earthquake Hazards Program (EHP) introduced a new automated Earthquake Notification Service (ENS) to take the place of the National Earthquake Information Center (NEIC) "Bigquake" system and the various other individual EHP e-mail list-servers for separate regions in the United States. These included northern California, southern California, and the central and eastern United States. ENS is a "one-stop shopping" system that allows Internet users to subscribe to flexible and customizable notifications for earthquakes anywhere in the world. The customization capability allows users to define the what (magnitude threshold), the when (day and night thresholds), and the where (specific regions) for their notifications. Customization is achieved by employing a per-user based request profile, allowing the notifications to be tailored for each individual's requirements. Such earthquake-parameter-specific custom delivery was not possible with simple e-mail list-servers. Now that event and user profiles are in a structured query language (SQL) database, additional flexibility is possible. At the time of this writing, ENS had more than 114,000 subscribers, with more than 200,000 separate user profiles. On a typical day, more than 188,000 messages get sent to a variety of widely distributed users for a wide range of earthquake locations and magnitudes. The purpose of this article is to describe how ENS works, highlight the features it offers, and summarize plans for future developments.
PAGER (Prompt Assessment of Global Earthquakes for Response) is an automated system to rapidly assess the number of people and regions exposed to severe shaking by an earthquake, and inform emergency responders, government agencies, and the media to the scope of the potential disaster. PAGER monitors the U.S. Geological Survey?s near real-time U.S. and global earthquake detections and automatically identifies events that are of societal importance, well in advance of ground-truth or news accounts.
Executive Summary It is well-known that large global earthquakes can have a dramatic effect on local communities and the built environment. Moreover, ground motions amplified by surficial materials can exacerbate the situation, often making the difference between minor and major damage. For a real-time earthquake impact alert system, such as Prompt Assessment of Global Earthquakes for Response (PAGER) (Wald and others, 2006), we seek to rapidly evaluate potential ground shaking in the source region and subsequently provide an estimate of the population exposure to potentially fatal levels of ground shaking in any region of the world. The contribution of surficial geology (particularly soft sediments) to the amplification of ground shaking is an important component in predicting the levels of ground motion observed at any site. Unfortunately, the availability of information regarding seismic siteconditions is only available at a few sites around the globe. Herein, we describe a methodology for deriving maps of seismic site-conditions anywhere in the world using topographic slope as a proxy. Average shear-velocity down to 30 m (or VS30) measurements are correlated against topographic slope to develop two sets of coefficients for predicting VS30: one for active tectonic regions that possess dynamic topographic relief, and one for stable continental regions where changes in topography are more subdued. These coefficients have been applied to the continental United States, in addition to other regions around the world. They are subsequently compared to existing site-condition maps based on geology and observed VS30 measurements, where available. The application of the topographic slope method in regions with abundant VS30 measurements (for example California, Memphis, and Taiwan) indicates that this method provides site condition-maps of similar quality, or in some cases, maps superior to those developed from more traditional techniques. Having a first-order assessment of seismic site-conditions anywhere in the world provides a valuable tool to rapidly estimate ground motions following any global earthquake, the primary motivation for this research. These VS30 maps will enable us to better quantify possible ground shaking and rapidly deliver these predictions to emergency managers and responders. In addition, the VS30 maps for the globe will also have practical applications for numerous related probabilistic- and scenario-based studies. To date, several researchers have requested maps or have used the approach outlined herein for their own applications (for example Cagnan and Kariptas, written commun., 2007; Harmandar and others, 2007). Given that we anticipate a significant demand for these products, we have developed an internet delivery service so that users can download maps and grids of seismic site-conditions for specified regions. To some extent, these grids can also be customized by the user if they disagree with the predefined correlations derived using the methodologies described within this report. Finally, this report represents a more comprehensive account of this technique and provides a more fully illustrated global description of results than that given in Wald and Allen (2007), which has been accepted for publication in the Bulletin of the Seismological Society of America.
Since the early 1990s, the magnitude and location of an earthquake have been available within minutes on the Internet. Now, as a result of work by the U.S. Geological Survey and with the cooperation of various regional seismic networks, people who experience an earthquake can go online and share information about its effects to help create a map of shaking intensities and damage. Such “Community Internet Intensity Maps” (CIIMs) contribute greatly toward the quick assessment of the scope of an earthquake emergency and provide valuable data for earthquake research.
PAGER (Prompt Assessment of Global Earthquakes for Response) is an automated alarm system being developed to rapidly and accurately assess the severity of damage caused by an earthquake and to provide emergency relief organizations, government agencies, and the media with an estimate of the societal impact from the potential catastrophe.
ShakeMap (http://earthquake.usgs.gov/shakemap) --rapidly, automatically generated shaking and intensity maps--combines instrumental measurements of shaking with information about local geology and earthquake location and magnitude to estimate shaking variations throughout a geographic area. The results are rapidly available via the Web through a variety of map formats, including Geographic Information System (GIS) coverages. These maps have become a valuable tool for emergency response, public information, loss estimation, earthquake planning, and post-earthquake engineering and scientific analyses. With the adoption of ShakeMap as a standard tool for a wide array of users and uses came an impressive demand for up-to-date technical documentation and more general guidelines for users and software developers. This manual is meant to address this need. ShakeMap, and associated Web and data products, are rapidly evolving as new advances in communications, earthquake science, and user needs drive improvements. As such, this documentation is organic in nature. We will make every effort to keep it current, but undoubtedly necessary changes in operational systems take precedence over producing and making documentation publishable.
[1] We investigate the rupture process of the 1999 Chi‐Chi, Taiwan, earthquake using extensive near‐source observations, including three‐component velocity waveforms at 36 strong motion stations and 119 GPS measurements. A three‐plane fault geometry derived from our previous inversion using only static data [ Ji et al. , 2001 ] is applied. The slip amplitude, rake angle, rupture initiation time, and risetime function are inverted simultaneously with a recently developed finite fault inverse method that combines a wavelet transform approach with a simulated annealing algorithm [ Ji et al. , 2002b ]. The inversion results are validated by the forward prediction of an independent data set, the teleseismic P and SH ground velocities, with notable agreement. The results show that the total seismic moment release of this earthquake is 2.7 × 10 20 N m and that most of the slip occurred in a triangular‐shaped asperity involving two fault segments, which is consistent with our previous static inversion. The rupture front propagates with an average rupture velocity of ∼2.0 km s −1 , and the average slip duration (risetime) is 7.2 s. Several interesting observations related to the temporal evolution of the Chi‐Chi earthquake are also investigated, including (1) the strong effect of the sinuous fault plane of the Chelungpu fault on spatial and temporal variations in slip history, (2) the intersection of fault 1 and fault 2 not being a strong impediment to the rupture propagation, and (3) the observation that the peak slip velocity near the surface is, in general, higher than on the deeper portion of the fault plane, as predicted by dynamic modeling.
We present isoseismal maps, macroseismic intensities, and community summaries of damage for the MW=6.8 Nisqually, Washington, earthquake of 28 February, 2001. For many communities, two types of macroseismic intensity are assigned, the traditional U.S. Geological Survey Modified Mercalli Intensities (USGS MMI) and a type of intensity newly introduced with this paper, the USGS Reviewed Community Internet Intensity (RCII). For most communities, the RCII is a reviewed version of the Community Internet Intensity (CII) of Wald and others (1999). For some communities, RCII is assigned from such non-CII sources as press reports, engineering reports, and field reconnaissance observations. We summarize differences between procedures used to assign RCII and USGS MMI, and we show that the two types of intensity are nonetheless very similar for the Nisqually earthquake. We do not see evidence for systematic differences between RCII and USGS MMI that would approach one intensity unit, at any level of shaking, but we document a tendency for the RCII to be slightly lower than MMI in regions of low intensity and slightly higher than MMI in regions of high intensity. The highest RCII calculated for the Nisqually earthquake is 7.6, calculated for zip code 98134, which includes the ?south of downtown? (Sodo) area of Seattle and Harbor Island. By comparison, we assigned a traditional USGS MMI 8 to the Sodo area of Seattle. In all, RCII of 6.5 and higher were assigned to 58 zip-code regions. At the lowest intensities, the Nisqually earthquake was felt over an area of approximately 350,000 square km (approximately 135,000 square miles) in Washington, Oregon, Idaho, Montana, and southern British Columbia, Canada. On the basis of macroseismic effects, we infer that shaking in the southern Puget Sound region was somewhat less for the 2001 Nisqually earthquake than for the Puget Sound earthquake of April 13, 1949, which had nearly the same hypocenter and magnitude. Allowing for differences in hypocenter, shaking in the 2001 earthquake was very similar to that produced by the Puget Sound earthquake of April 25, 1965. First-person accounts of the effects of the 2001 earthquake on individual households are given for some communities.
The San Andreas fault is the longest fault in California and one of the longest strikeslip faults in the world, yet little is known about the aftershocks following the most recent great event on the San Andreas, the M 7.8 San Francisco earthquake, on 18 April 1906. This open-file report is a compilation of first-hand accounts (felt reports) describing aftershocks and triggered events of the 1906 earthquake, for the first twenty months of the aftershock sequence (through December 1907). The report includes a chronological catalog. For the larger events, Modified Mercalli intensities (MMIs) have been assigned based on the descriptions judged to be the most reliable.
Rapid (3-5 minutes) generation of maps of instrumental ground-motion and shaking intensity is accomplished through advances in real-time seismographic data acquisition combined with newly developed relationships between recorded ground-motion parameters and expected shaking intensity values. Estimation of shaking over the entire regional extent of southern California is obtained by the spatial interpolation of the measured ground motions with geologically based frequency and amplitude-dependent site corrections. Production of the maps is automatic, triggered by any significant earthquake in southern California. Maps are now made available within several minutes of the earthquake for public and scientific consumption via the World Wide Web; they will be made available with dedicated communications for emergency response agencies and critical users.
We have developed regression relationships between Modified Mercalli Intensity ( I mm ) and peak ground acceleration (PGA) and velocity (PGV) by comparing horizontal peak ground motions to observed intensities for eight significant California earthquakes. For the limited range of Modified Mercalli intensities ( I mm ), we find that for peak acceleration with V ≤ I mm ≤ VIII, I mm = 3.66 log ( PGA ) − 1.66, and for peak velocity with V ≤ I mm ≤ IX, I mm = 3.47 log ( PGV ) + 2.35. From comparison with observed intensity maps, we find that a combined regression based on peak velocity for intensity > VII and on peak acceleration for intensity < VII is most suitable for reproducing observed I mm patterns, consistent with high intensities being related to damage (proportional to ground velocity) and with lower intensities determined by felt accounts (most sensitive to higher-frequency ground acceleration). These new I mm relationships are significantly different from the Trifunac and Brady (1975) correlations, which have been used extensively in loss estimation.
No abstract available.
Explore the source record for details and available documents.
The October 17, 1989, Loma Prieta, Calif., earthquake (0004:15.2 G.m.t. October 18; lat 37.036º N., long 121.883º W.; 19-km depth) had a local magnitude (M L ) of about 6.7, a surface-wave magnitude (M S ) of 7.1, a seismic moment of 2.2x10 19 N-m to 3.5x10 19 N-m, a source duration of 6 to 15 s, and an average stress drop of at least 50 bars. Slip occurred on a dipping fault surface about 35 km long and was largely confined to a depth of about 7 to 20 km. The slip vector had a large vertical component, and slip was distributed in two main regions situated northwest and southeast of the hypocenter. This slip distribution caused about half of the earthquake's energy to be focused toward the urbanized San Francisco Bay region, while the other half was focused toward the southeast. Had the rupture initiated at the southeast end of the aftershock zone, shaking in the bay region would have been both longer and stronger. These source parameters suggest that the earthquake was not a typical shallow San Andreas-type event but a deeper event on a different fault with a recurrence interval of many hundreds of years. Therefore, the potential for a damaging shallow event on the San Andreas fault in the Santa Cruz Mountains may still exist.