Robert L. Wesson

USGS National Seismic Hazard Maps

The U.S. Geological Survey (USGS) recently completed new probabilistic seismic hazard maps for the United States, including Alaska and Hawaii. These hazard maps form the basis of the probabilistic component of the design maps used in the 1997 edition of the NEHRP Recommended Provisions for Seismic Regulations for New Buildings and Other Structures, prepared by the Building Seismic Safety Council and published by FEMA. The hazard maps depict peak horizontal ground acceleration and spectral response at 0.2, 0.3, and 1.0 sec periods, with 10%, 5%, and 2% probabilities of exceedance in 50 years, corresponding to return times of about 500, 1000, and 2500 years, respectively. In this paper we outline the methodology used to construct the hazard maps. There are three basic components to the maps. First, we use spatially smoothed historic seismicity as one portion of the hazard calculation. In this model, we apply the general observation that moderate and large earthquakes tend to occur near areas of previous small or moderate events, with some notable exceptions. Second, we consider large background source zones based on broad geologic criteria to quantify hazard in areas with little or no historic seismicity, but with the potential for generating large events. Third, we include the hazard from specific fault sources. We use about 450 faults in the western United States (WUS) and derive recurrence times from either geologic slip rates or the dating of pre-historic earthquakes from trenching of faults or other paleoseismic methods. Recurrence estimates for large earthquakes in New Madrid and Charleston, South Carolina, were taken from recent paleoliquefaction studies. We used logic trees to incorporate different seismicity models, fault recurrence models, Cascadia great earthquake scenarios, and ground-motion attenuation relations. We present disaggregation plots showing the contribution to hazard at four cities from potential earthquakes with various magnitudes and distances.

Triggered earthquakes and deep well activities

Earthquakes can be triggered by any significant perturbation of the hydrologic regime. In areas where potentially active faults are already close to failure, the increased pore pressure resulting from fluid injection, or, alternatively, the massive extraction of fluid or gas, can induce sufficient stress and/or strain changes that, with time, can lead to sudden catastrophic failure in a major earthquake. Injection-induced earthquakes typically result from the reduction in frictional strength along preexisting, nearby faults caused by the increased formation fluid pressure. Earthquakes associated with production appear to respond to more complex mechanisms of subsidence, crustal unloading, and poroelastic changes in response to applied strains induced by the massive withdrawal of subsurface material. As each of these different types of triggered events can occur up to several years after well activities have begun (or even several years after all well activities have stopped), this suggests that the actual triggering process may be a very complex combination of effects, particularly if both fluid extraction and injection have taken place locally. To date, more than thirty cases of earthquakes triggered by well activities can be documented throughout the United States and Canada. Based on these case histories, it is evident that, owing to preexisting stress conditions in the upper crust, certain areas tend to have higher probabilities of exhibiting such induced seismicity.

Overlapping fault planes of the 1971 San Fernando and 1994 Northridge, California earthquakes

Aftershocks of the 1971 San Fernando and 1994 Northridge earthquakes were relocated using a three-dimensional velocity model that was derived from inverting P-wave travel time data. The hypocenters show clear orientations of the dipping fault planes. The San Fernando aftershocks form a plane extending from a depth of 15 km to the surface, dipping toward the northeast at about 40°. The Northridge aftershocks delineate a fault extending from a depth of 18 km up to about 5 km, dipping toward the southwest at about 40°. In the region the aftershocks overlap in map view, the San Fernando plane cuts off the Northridge plane at a depth of 5 to 8 km, preventing it from reaching the surface. The similar but oppositely dipping fault planes suggest a pair of conjugate planes reflecting a horizontal northeast-southwest compression.

Variation of P-Wave Velocity before the Bear Valley, California, Earthquake of 24 February 1972

Residuals for P-wave traveltimes at a seismnograph station near Bear Valley, California, for small, precisely located local earthquakes at distances of 20 to 70 kilometers show a sharp increase of nearly 0.3 second about 2 months before a magnitude 5.0 earthquake that occurred within a few kilometers of the station. This indicates that velocity changes observed elsewhere premonitory to earthquakes, possibly related to dilatancy, occur along the central section of the San Andreas fault system.

Splay fault slip during the Mw 8.8 2010 Maule Chile earthquake

Splay faults are thrusts that emerge from the plate boundaries of subduction zones. Such structures have been mapped at several convergent margins and their activity commonly ascribed to large megathrust earthquakes. However, the behavior of splay faults during the earthquake cycle is poorly constrained because typically these structures are located offshore and are difficult to access. Here we use geologic mapping combined with space and land geodesy, as well as offshore sonar data, to document surface-fault ruptures and coastal uplift at Isla Santa Maria in south-central Chile (37°S) caused by the 27 February 2010 Maule earthquake ( M w 8.8). During the earthquake, the island was tilted parallel to the margin, and normal faults ruptured the surface and adjacent ocean bottom. We associate tilt and crestal normal faulting with growth of an anticline above a blind reverse fault rooted in the Nazca–South America plate boundary, which slipped during the Maule earthquake. The splay fault system has formed in an area of reduced coseismic plate-boundary slip, suggesting that anelastic deformation in the upper plate may have restrained the 2010 megathrust rupture. Surface fault breaks were accompanied by prominent discharge of fluids. Our field observations support the notion that splay faulting may frequently complement and influence the rupture of subduction-zone earthquakes.

Spatial Correlation of Probabilistic Earthquake Ground Motion and Loss

Spatial correlation of annual earthquake ground motions and losses can be used to estimate the variance of annual losses to a portfolio of properties exposed to earthquakes. A direct method is described for the calculation of the spatial correlation of earthquake ground motions and losses. Calculations for the direct method can be carried out using either numerical quadrature or a discrete, matrix-based approach. Numerical results for this method are compared with those calculated from a simple Monte Carlo simulation. Spatial correlation of ground motion and loss is induced by the systematic attenuation of ground motion with distance from the source, by common site conditions, and by the finite length of fault ruptures. Spatial correlation is also strongly dependent on the partitioning of the variability, given an event, into interevent and intraevent components. Intraevent variability reduces the spatial correlation of losses. Interevent variability increases spatial correlation of losses. The higher the spatial correlation, the larger the variance in losses to a portfolio, and the more likely extreme values become. This result underscores the importance of accurately determining the relative magnitudes of intraevent and interevent variability in ground-motion studies, because of the strong impact in estimating earthquake losses to a portfolio. The direct method offers an alternative to simulation for calculating the variance of losses to a portfolio, which may reduce the amount of calculation required.

Seismic Hazard in Hawaii: High Rate of Large Earthquakes and Probabilistic Ground-Motion Maps

The seismic hazard and earthquake occurrence rates in Hawaii are locally as high as that near the most hazardous faults elsewhere in the United States. We have generated maps of peak ground acceleration (PGA) and spectral acceleration (SA) (at 0.2, 0.3 and 1.0 sec, 5% critical damping) at 2% and 10% exceedance probabilities in 50 years. The highest hazard is on the south side of Hawaii Island, as indicated by the M I 7.0, M S 7.2, and M I 7.9 earthquakes, which occurred there since 1868. Probabilistic values of horizontal PGA (2% in 50 years) on Hawaiis south coast exceed 1.75 g . Because some large earthquake aftershock zones and the geometry of flank blocks slipping on subhorizontal decollement faults are known, we use a combination of spatially uniform sources in active flank blocks and smoothed seismicity in other areas to model seismicity. Rates of earthquakes are derived from magnitude distributions of the modern (1959–1997) catalog of the Hawaiian Volcano Observatorys seismic network supplemented by the historic (1868–1959) catalog. Modern magnitudes are M L measured on a Wood-Anderson seismograph or M S. Historic magnitudes may add M L measured on a Milne-Shaw or Bosch-Omori seismograph or M I derived from calibrated areas of MM intensities. Active flank areas, which by far account for the highest hazard, are characterized by distributions with b slopes of about 1.0 below M 5.0 and about 0.6 above M 5.0. The kinked distribution means that large earthquake rates would be grossly underestimated by extrapolating small earthquake rates, and that longer catalogs are essential for estimating or verifying the rates of large earthquakes. Flank earthquakes thus follow a semicharacteristic model, which is a combination of background seismicity and an excess number of large earthquakes. Flank earthquakes are geometrically confined to rupture zones on the volcano flanks by barriers such as rift zones and the seaward edge of the volcano, which may be expressed by a magnitude distribution similar to that including characteristic earthquakes. The island chain northwest of Hawaii Island is seismically and volcanically much less active. We model its seismic hazard with a combination of a linearly decaying ramp fit to the cataloged seismicity and spatially smoothed seismicity with a smoothing half-width of 10 km. We use a combination of up to four attenuation relations for each map because for either PGA or SA, there is no single relation that represents ground motion for all distance and magnitude ranges. Great slumps and landslides visible on the ocean floor correspond to catastrophes with effective energy magnitudes M E above 8.0. A crude estimate of their frequency suggests that the probabilistic earthquake hazard is at least an order of magnitude higher for flank earthquakes than that from submarine slumps.

Modelling aftershock migration and afterslip of the San Juan Bautista, California, earthquake of October 3, 1972

Abstract The San Juan Bautista earthquake of October 3, 1972 ( M L = 4.8), located along the San Andreas fault in central California, initiated an aftershock sequence characterized by a subtle, but perceptible, tendency for aftershocks to spread to the northwest and southeast along the fault zone. The apparent dimension of the aftershock zone along strike increased from about 7–10 km within a few days of the earthquake, to about 20 km eight months later. In addition, the mainshock initiated a period of accelerated fault creep, which was observed at 2 creep meters situated astride the trace of the San Andreas fault within about 15 km of the epicenter of the mainshock. The creep rate gradually returned to the preearthquake rate after about 3 yrs. Both the spreading of the aftershocks and the rapid surface creep are interpreted as reflecting a period of rapid creep in the fault zone representing the readjustment of stress and displacement following the failure of a “stuck” patch or asperity during the San Juan Bautista earthquake. Numerical calculations suggest that the behavior of the fault zone is consistent with that of a material characterized by a viscosity of about 3.6×10 14 P, although the real rheology is likely to be more complicated. In this model, the mainshock represents the failure of an asperity that slips only during earthquakes. Aftershocks represent the failure of second-order asperities which are dragged along by the creeping fault zone.

Water-level fluctuations and earthquakes on the San Andreas fault zone

Observations of the water-level changes in a well drilled into the San Andreas fault zone have been under way since May 1971, with the objective of studying in situ pore-pressure changes in a zone of active tectonic creep and seismicity. Small water-level changes, characterized by a decrease and subsequent rise, have been followed by earthquakes of moderate size on the San Andreas fault zone. Compatibility of these observations with either a dilatancy-type behavior or a dislocation-type behavior for the pre-earthquake process can be demonstrated. Additional water-level observations at other sites in the fault zone are needed to examine the spatial and temporal extent of the actual preseismic process that is responsible for observed creep events, water-level changes, and tilt changes.

Five centuries of tsunamis and land-level changes in the overlapping rupture area of the 1960 and 2010 Chilean earthquakes

A combination of geological and historical records from south-central Chile provides a means to address general questions about the stability of megathrust rupture patches and the range of variation expected among earthquakes and tsunamis along a particular stretch of a subduction zone. The Tirua River estuary (38.3°S) records four large tsunamis and coseismic land-level changes over the past 450 years within the overlapping rupture zones of the great subduction-zone earthquakes of A.D. 1960 (Mw 9.5) and 2010 (Mw 8.8). Sand layers 2 km up the Tirua River represent the 2010 and 1960 tsunamis and two historical tsunamis, most likely in A.D. 1751 and 1575. Differing land-level changes during these earthquakes likely denote differences in the spatial distribution of slip on the megathrust in both the strike and dip directions within the overlapping rupture zone, with the uplift at Tirua in 1751 and 2010 probably caused by slip extending farther landward and to greater depth than in 1575 and 1960, which showed subsidence or little change.