Geology topics

D.J. Wald

Publications and source records attributed to D.J. Wald.

At least 19 recordsLinked to original sources

USGS2014

Estimating structural collapse fragility of generic building typologies using expert judgment

The structured expert elicitation process proposed by Cooke (1991), hereafter referred to as Cooke’s approach, is applied for the first time in the realm of structural collapse-fragility assessment for selected generic construction types. Cooke’s approach works on the principle of objective calibration scoring of judgments coupled with hypothesis testing used in classical statistics. The performance-based scoring system reflects the combined measure of an expert’s informativeness about variables in the problem area under consideration, and their ability to enumerate, in a statistically accurate way through expressing their true beliefs, the quantitative uncertainties associated with their assessments. We summarize the findings of an expert elicitation workshop in which a dozen earthquake-engineering professionals from around the world were engaged to estimate seismic collapse fragility for generic construction types. Development of seismic collapse fragility functions was accomplished by combining their judgments using weights derived from Cooke’s method. Although substantial effort was needed to elicit the inputs of these experts successfully, we anticipate that the elicitation strategy described here will gain momentum in a wide variety of earthquake seismology and engineering hazard and risk analyses where physical model and data limitations are inherent and objective professional judgment can fill gaps.

Conference Paper↗

USGS2013

Strategies for rapid global earthquake impact estimation: the Prompt Assessment of Global Earthquakes for Response (PAGER) system

This chapter summarizes the state-of-the-art for rapid earthquake impact estimation. It details the needs and challenges associated with quick estimation of earthquake losses following global earthquakes, and provides a brief literature review of various approaches that have been used in the past. With this background, the chapter introduces the operational earthquake loss estimation system developed by the U.S. Geological Survey (USGS) known as PAGER (for Prompt Assessment of Global Earthquakes for Response). It also details some of the ongoing developments of PAGER’s loss estimation models to better supplement the operational empirical models, and to produce value-added web content for a variety of PAGER users.

Book chapter↗

USGS2012

Using shakecast and shakemap for lifeline post-earthquake response and earthquake scenario planning

The U. S. Geological Survey (USGS) ShakeCast system is a freely-available, post-earthquake situational awareness application that automatically retrieves earthquake shaking data from ShakeMap, compares intensity measures against users' facilities, sends notifications of potential damage to responsible parties, and generates facility damage assessment maps and other Web-based products for emergency managers and responders. The ShakeMap/ShakeCast combined system allows a consistent approach for the evaluation of facility performance using the ShakeMap methodology by combining observations with ground motion predictions and covers major geographic regions and earthquakes, real-time, historical, and scenario, worldwide. We present common approaches to assessing facility impact estimates by combining measured or estimated ground-motion parameters with vulnerability functions. Vulnerability functions can be pre-computed, use one of these ground-motion parameters as input, and produce a multi-state discrete output of damage likelihood (or, inspection priorities). Two custom modules for evaluating Arias Intensity and Newmark landslide displacements are presented as examples of using the ShakeCast application with more advanced, customized damage functions. ?? 2009 ASCE.

Conference Paper↗

USGS2012

USGS "Did You Feel It?" internet-based macroseismic intensity maps

The U.S. Geological Survey (USGS) "Did You Feel It?" (DYFI) system is an automated approach for rapidly collecting macroseismic intensity data from Internet users' shaking and damage reports and generating intensity maps immediately following earthquakes; it has been operating for over a decade (1999-2011). DYFI-based intensity maps made rapidly available through the DYFI system fundamentally depart from more traditional maps made available in the past. The maps are made more quickly, provide more complete coverage and higher resolution, provide for citizen input and interaction, and allow data collection at rates and quantities never before considered. These aspects of Internet data collection, in turn, allow for data analyses, graphics, and ways to communicate with the public, opportunities not possible with traditional data-collection approaches. Yet web-based contributions also pose considerable challenges, as discussed herein. After a decade of operational experience with the DYFI system and users, we document refinements to the processing and algorithmic procedures since DYFI was first conceived. We also describe a number of automatic post-processing tools, operations, applications, and research directions, all of which utilize the extensive DYFI intensity datasets now gathered in near-real time. DYFI can be found online at the website http://earthquake.usgs.gov/dyfi/. ?? 2011 by the Istituto Nazionale di Geofisica e Vulcanologia.

Annals of Geophysics↗

USGS2012

Impact-based earthquake alerts with the U.S. Geological Survey's PAGER system: what's next?

In September 2010, the USGS began publicly releasing earthquake alerts for significant earthquakes around the globe based on estimates of potential casualties and economic losses with its Prompt Assessment of Global Earthquakes for Response (PAGER) system. These estimates significantly enhanced the utility of the USGS PAGER system which had been, since 2006, providing estimated population exposures to specific shaking intensities. Quantifying earthquake impacts and communicating estimated losses (and their uncertainties) to the public, the media, humanitarian, and response communities required a new protocol—necessitating the development of an Earthquake Impact Scale—described herein and now deployed with the PAGER system. After two years of PAGER-based impact alerting, we now review operations, hazard calculations, loss models, alerting protocols, and our success rate for recent (2010-2011) events. This review prompts analyses of the strengths, limitations, opportunities, and pressures, allowing clearer definition of future research and development priorities for the PAGER system.

Conference Paper↗

USGS2012

Improving PAGER's real-time earthquake casualty and loss estimation toolkit: a challenge

We describe the on-going developments of PAGER’s loss estimation models, and discuss value-added web content that can be generated related to exposure, damage and loss outputs for a variety of PAGER users. These developments include identifying vulnerable building types in any given area, estimating earthquake-induced damage and loss statistics by building type, and developing visualization aids that help locate areas of concern for improving post-earthquake response efforts. While detailed exposure and damage information is highly useful and desirable, significant improvements are still necessary in order to improve underlying building stock and vulnerability data at a global scale. Existing efforts with the GEM’s GED4GEM and GVC consortia will help achieve some of these objectives. This will benefit PAGER especially in regions where PAGER’s empirical model is less-well constrained; there, the semi-empirical and analytical models will provide robust estimates of damage and losses. Finally, we outline some of the challenges associated with rapid casualty and loss estimation that we experienced while responding to recent large earthquakes worldwide.

Conference Paper↗

USGS2012

Developing ShakeCast statistical fragility analysis framework for rapid post-earthquake assessment

When an earthquake occurs, the U. S. Geological Survey (USGS) ShakeMap estimates the extent of potentially damaging shaking and provides overall information regarding the affected areas. The USGS ShakeCast system is a freely-available, post-earthquake situational awareness application that automatically retrieves earthquake shaking data from ShakeMap, compares intensity measures against users’ facilities, sends notifications of potential damage to responsible parties, and generates facility damage assessment maps and other web-based products for emergency managers and responders. We describe notable improvements of the ShakeMap and the ShakeCast applications. We present a design for comprehensive fragility implementation, integrating spatially-varying ground-motion uncertainties into fragility curves for ShakeCast operations. For each facility, an overall inspection priority (or damage assessment) is assigned on the basis of combined component-based fragility curves using pre-defined logic. While regular ShakeCast users receive overall inspection priority designations for each facility, engineers can access the full fragility analyses for further evaluation.

Conference Paper↗

USGS2012

Developing Vs30 site-condition maps by combining observations with geologic and topographic constraints

Despite obvious limitations as a proxy for site amplification, the use of time-averaged shear-wave velocity over the top 30 m (V S30 ) remains widely practiced, most notably through its use as an explanatory variable in ground motion prediction equations (and thus hazard maps and ShakeMaps, among other applications). As such, we are developing an improved strategy for producing V S30 maps given the common observational constraints. Using the abundant V S30 measurements in Taiwan, we compare alternative mapping methods that combine topographic slope, surface geology, and spatial correlation structure. The different V S30 mapping algorithms are distinguished by the way that slope and geology are combined to define a spatial model of V S30 . We consider the globally applicable slope-only model as a baseline to which we compare two methods of combining both slope and geology. For both hybrid approaches, we model spatial correlation structure of the residuals using the kriging-with-a-trend technique, which brings the map into closer agreement with the observations. Cross validation indicates that we can reduce the uncertainty of the V S30 map by up to 16% relative to the slope-only approach.

Conference Paper↗

USGS2012

ShakeMap Atlas 2.0: an improved suite of recent historical earthquake ShakeMaps for global hazard analyses and loss model calibration

We introduce the second version of the U.S. Geological Survey ShakeMap Atlas, which is an openly-available compilation of nearly 8,000 ShakeMaps of the most significant global earthquakes between 1973 and 2011. This revision of the Atlas includes: (1) a new version of the ShakeMap software that improves data usage and uncertainty estimations; (2) an updated earthquake source catalogue that includes regional locations and finite fault models; (3) a refined strategy to select prediction and conversion equations based on a new seismotectonic regionalization scheme; and (4) vastly more macroseismic intensity and ground-motion data from regional agencies All these changes make the new Atlas a self-consistent, calibrated ShakeMap catalogue that constitutes an invaluable resource for investigating near-source strong ground-motion, as well as for seismic hazard, scenario, risk, and loss-model development. To this end, the Atlas will provide a hazard base layer for PAGER loss calibration and for the Earthquake Consequences Database within the Global Earthquake Model initiative.

Conference Paper↗

USGS2010

A revised ground-motion and intensity interpolation scheme for shakemap

We describe a weighted-average approach for incorporating various types of data (observed peak ground motions and intensities and estimates from groundmotion prediction equations) into the ShakeMap ground motion and intensity mapping framework. This approach represents a fundamental revision of our existing ShakeMap methodology. In addition, the increased availability of near-real-time macroseismic intensity data, the development of newrelationships between intensity and peak ground motions, and new relationships to directly predict intensity from earthquake source information have facilitated the inclusion of intensity measurements directly into ShakeMap computations. Our approach allows for the combination of (1) direct observations (ground-motion measurements or reported intensities), (2) observations converted from intensity to ground motion (or vice versa), and (3) estimated ground motions and intensities from prediction equations or numerical models. Critically, each of the aforementioned data types must include an estimate of its uncertainties, including those caused by scaling the influence of observations to surrounding grid points and those associated with estimates given an unknown fault geometry. The ShakeMap ground-motion and intensity estimates are an uncertainty-weighted combination of these various data and estimates. A natural by-product of this interpolation process is an estimate of total uncertainty at each point on the map, which can be vital for comprehensive inventory loss calculations. We perform a number of tests to validate this new methodology and find that it produces a substantial improvement in the accuracy of ground-motion predictions over empirical prediction equations alone.

Bulletin of the Seismological Society of America↗

USGS2010

PAGER--Rapid assessment of an earthquakes impact

PAGER (Prompt Assessment of Global Earthquakes for Response) is an automated system that produces content concerning the impact of significant earthquakes around the world, informing emergency responders, government and aid agencies, and the media of the scope of the potential disaster. PAGER rapidly assesses earthquake impacts by comparing the population exposed to each level of shaking intensity with models of economic and fatality losses based on past earthquakes in each country or region of the world. Earthquake alerts--which were formerly sent based only on event magnitude and location, or population exposure to shaking--now will also be generated based on the estimated range of fatalities and economic losses.

Fact Sheet↗

USGS2009

Advancing techniques to constrain the geometry of the seismic rupture plane on subduction interfaces a priori: Higher-order functional fits

Ongoing developments in earthquake source inversions incorporate nonplanar fault geometries as inputs to the inversion process, improving previous approaches that relied solely on planar fault surfaces. This evolution motivates advancing the existing framework for constraining fault geometry, particularly in subduction zones where plate boundary surfaces that host highly hazardous earthquakes are clearly nonplanar. Here, we improve upon the existing framework for the constraint of the seismic rupture plane of subduction interfaces by incorporating active seismic and seafloor sediment thickness data with existing independent data sets and inverting for the most probable nonplanar subduction geometry. Constraining the rupture interface a priori with independent geological and seismological information reduces the uncertainty in the derived earthquake source inversion parameters over models that rely on simpler assumptions, such as the moment tensor inferred fault plane. Examples are shown for a number of wellconstrained global locations. We expand the coverage of previous analyses to a more uniform global data set and show that even in areas of sparse data this approach is able to accurately constrain the approximate subduction geometry, particularly when aided with the addition of data from local active seismic surveys. In addition, we show an example of the integration of many two-dimensional profiles into a threedimensional surface for the Sunda subduction zone and introduce the development of a new global threedimensional subduction interface model: Slab1.0. ?? 2009 by the American Geophysical Union.

Geochemistry, Geophysics, Geosystems↗

USGS2009

An Atlas of ShakeMaps and population exposure catalog for earthquake loss modeling

We present an Atlas of ShakeMaps and a catalog of human population exposures to moderate-to-strong ground shaking (EXPO-CAT) for recent historical earthquakes (1973-2007). The common purpose of the Atlas and exposure catalog is to calibrate earthquake loss models to be used in the US Geological Survey's Prompt Assessment of Global Earthquakes for Response (PAGER). The full ShakeMap Atlas currently comprises over 5,600 earthquakes from January 1973 through December 2007, with almost 500 of these maps constrained-to varying degrees-by instrumental ground motions, macroseismic intensity data, community internet intensity observations, and published earthquake rupture models. The catalog of human exposures is derived using current PAGER methodologies. Exposure to discrete levels of shaking intensity is obtained by correlating Atlas ShakeMaps with a global population database. Combining this population exposure dataset with historical earthquake loss data, such as PAGER-CAT, provides a useful resource for calibrating loss methodologies against a systematically-derived set of ShakeMap hazard outputs. We illustrate two example uses for EXPO-CAT; (1) simple objective ranking of country vulnerability to earthquakes, and; (2) the influence of time-of-day on earthquake mortality. In general, we observe that countries in similar geographic regions with similar construction practices tend to cluster spatially in terms of relative vulnerability. We also find little quantitative evidence to suggest that time-of-day is a significant factor in earthquake mortality. Moreover, earthquake mortality appears to be more systematically linked to the population exposed to severe ground shaking (Modified Mercalli Intensity VIII+). Finally, equipped with the full Atlas of ShakeMaps, we merge each of these maps and find the maximum estimated peak ground acceleration at any grid point in the world for the past 35 years. We subsequently compare this "composite ShakeMap" with existing global hazard models, calculating the spatial area of the existing hazard maps exceeded by the combined ShakeMap ground motions. In general, these analyses suggest that existing global, and regional, hazard maps tend to overestimate hazard. Both the Atlas of ShakeMaps and EXPO-CAT have many potential uses for examining earthquake risk and epidemiology. All of the datasets discussed herein are available for download on the PAGER Web page ( http://earthquake.usgs.gov/ eqcenter/pager/prodandref/ ). ?? 2009 Springer Science+Business Media B.V.

Bulletin of Earthquake Engineering↗

USGS2009

Developing framework to constrain the geometry of the seismic rupture plane on subduction interfaces a priori - A probabilistic approach

A key step in many earthquake source inversions requires knowledge of the geometry of the fault surface on which the earthquake occurred. Our knowledge of this surface is often uncertain, however, and as a result fault geometry misinterpretation can map into significant error in the final temporal and spatial slip patterns of these inversions. Relying solely on an initial hypocentre and CMT mechanism can be problematic when establishing rupture characteristics needed for rapid tsunami and ground shaking estimates. Here, we attempt to improve the quality of fast finite-fault inversion results by combining several independent and complementary data sets to more accurately constrain the geometry of the seismic rupture plane of subducting slabs. Unlike previous analyses aimed at defining the general form of the plate interface, we require mechanisms and locations of the seismicity considered in our inversions to be consistent with their occurrence on the plate interface, by limiting events to those with well-constrained depths and with CMT solutions indicative of shallow-dip thrust faulting. We construct probability density functions about each location based on formal assumptions of their depth uncertainty and use these constraints to solve for the ‘most-likely’ fault plane. Examples are shown for the trench in the source region of the M w 8.6 Southern Sumatra earthquake of March 2005, and for the Northern Chile Trench in the source region of the November 2007 Antofagasta earthquake. We also show examples using only the historic catalogues in regions without recent great earthquakes, such as the Japan and Kamchatka Trenches. In most cases, this method produces a fault plane that is more consistent with all of the data available than is the plane implied by the initial hypocentre and CMT mechanism. Using the aggregated data sets, we have developed an algorithm to rapidly determine more accurate initial fault plane geometries for source inversions of future earthquakes.

Geophysical Journal International↗

USGS2009

On the use of high-resolution topographic data as a proxy for seismic site conditions (V S30 )

An alternative method has recently been proposed for evaluating global seismic site conditions, or the average shear velocity to 30 m depth (V S30 ), from the Shuttle Radar Topography Mission (SRTM) 30 arcsec digital elevation models (DEMs). The basic premise of the method is that the topographic slope can be used as a reliable proxy for V S30 in the absence of geologically and geotechnically based site-condition maps through correlations between VS30 measurements and topographic gradient. Here we evaluate the use of higher-resolution (3 and 9 arcsec) DEMs to examine whether we are able to resolve V S30 in more detail than can be achieved using the lower-resolution SRTM data. High-quality DEMs at resolutions greater than 30 arcsec are not uniformly available at the global scale. However, in many regions where such data exist, they may be employed to resolve finer-scale variations in topographic gradient, and consequently, V S30 . We use the U.S. Geological Survey Earth Resources Observation and Science (EROS) Data Center's National Elevation Dataset (NED) to investigate the use of high-resolution DEMs for estimating V S30 in several regions across the United States, including the San Francisco Bay area in California, Los Angeles, California, and St. Louis, Missouri. We compare these results with an example from Taipei, Taiwan, that uses 9 arcsec SRTM data, which are globally available. The use of higher-resolution NED data recovers finer-scale variations in topographic gradient, which better correlate to geological and geomorphic features, in particular, at the transition between hills and basins, warranting their use over 30 arcsec SRTM data where available. However, statistical analyses indicate little to no improvement over lower-resolution topography when compared to V S30 measurements, suggesting that some topographic smoothing may provide more stable V S30 estimates. Furthermore, we find that elevation variability in canopy-based SRTM measurements at resolutions greater than 30 arcsec are too large to resolve reliable slopes, particularly in low-gradient sedimentary basins.

Bulletin of the Seismological Society of America↗

USGS2009

PAGER-CAT: A composite earthquake catalog for calibrating global fatality models

We have described the compilation and contents of PAGER-CAT, an earthquake catalog developed principally for calibrating earthquake fatality models. It brings together information from a range of sources in a comprehensive, easy to use digital format. Earthquake source information (e.g., origin time, hypocenter, and magnitude) contained in PAGER-CAT has been used to develop an Atlas of Shake Maps of historical earthquakes (Allen et al. 2008) that can subsequently be used to estimate the population exposed to various levels of ground shaking (Wald et al. 2008). These measures will ultimately yield improved earthquake loss models employing the uniform hazard mapping methods of ShakeMap. Currently PAGER-CAT does not consistently contain indicators of landslide and liquefaction occurrence prior to 1973. In future PAGER-CAT releases we plan to better document the incidence of these secondary hazards. This information is contained in some existing global catalogs but is far from complete and often difficult to parse. Landslide and liquefaction hazards can be important factors contributing to earthquake losses (e.g., Marano et al. unpublished). Consequently, the absence of secondary hazard indicators in PAGER-CAT, particularly for events prior to 1973, could be misleading to sorne users concerned with ground-shaking-related losses. We have applied our best judgment in the selection of PAGER-CAT's preferred source parameters and earthquake effects. We acknowledge the creation of a composite catalog always requires subjective decisions, but we believe PAGER-CAT represents a significant step forward in bringing together the best available estimates of earthquake source parameters and reports of earthquake effects. All information considered in PAGER-CAT is stored as provided in its native catalog so that other users can modify PAGER preferred parameters based on their specific needs or opinions. As with all catalogs, the values of some parameters listed in PAGER-CAT are highly uncertain, particularly the casualty numbers, which must be regarded as estimates rather than firm numbers for many earthquakes. Consequently, we encourage contributions from the seismology and earthquake engineering communities to further improve this resource via the Wikipedia page and personal communications, for the benefit of the whole community.

Seismological Research Letters↗

USGS2008

Rapid Assessment of earthquake-induced landsliding

The Pacific Northwest in the United States including Seattle, Washington, experienced unusually heavy rainfall in the winters of 1995/1996 and 1996/1997, which caused numerous landslides. Following these two winters, the City of Seattle resolved to reduce future landslide losses within its jurisdiction. By coincidence, in 1997 the U.S. Geological Survey (USGS) began a five-year project designed to increase scientific understanding of all the natural hazards facing the Seattle area, including landslides. This article documents the convergence of the efforts of the City of Seattle, the State of Washington, and the USGS to understand the landslide hazards facing the Seattle area and to implement hazard reduction and mitigation of damage from landslides. Specifically, this article identifies how and where USGS information on landslide hazards fit in with efforts already undertaken by the City. This combined effort enabled the City of Seattle to understand its landslide hazards, to formulate public policy, to educate its citizens, and to coordinate the activities of the numerous departments under its jurisdiction with the ultimate goal of lessening the impact of landslides. The City used USGS maps and reports to guide decisions that prioritized funds for siting, maintaining and constructing public facilities in addition to meeting Washington State laws that mandated growth management and environmental review.

Conference Paper↗

USGS2008

WHE-PAGER Project: A new initiative in estimating global building inventory and its seismic vulnerability

The U.S. Geological Survey’s Prompt Assessment of Global Earthquake’s Response (PAGER) Project and the Earthquake Engineering Research Institute’s World Housing Encyclopedia (WHE) are creating a global database of building stocks and their earthquake vulnerability. The WHE already represents a growing, community-developed public database of global housing and its detailed structural characteristics. It currently contains more than 135 reports on particular housing types in 40 countries. The WHE-PAGER effort extends the WHE in several ways: (1) by addressing non-residential construction; (2) by quantifying the prevalence of each building type in both rural and urban areas; (3) by addressing day and night occupancy patterns, (4) by adding quantitative vulnerability estimates from judgment or statistical observation; and (5) by analytically deriving alternative vulnerability estimates using in part laboratory testing.

Conference Paper↗