Geology Reports⌕ Search

Geology topics

Oliver S. Boyd

Publications and source records attributed to Oliver S. Boyd.

71 records · Page 4Linked to original sources

Comparison of smoothing methods for the development of a smoothed seismicity model for Alaska and the implications for seismic hazard

In anticipation of the update of the Alaska seismic hazard maps (ASHMs) by the U. S. Geological Survey, we report progress on the comparison of smoothed seismicity models developed using fixed and adaptive smoothing algorithms, and investigate the sensitivity of seismic hazard to the models. While fault-based sources, such as those for great earthquakes in the Alaska-Aleutian subduction zone and for the ~10 shallow crustal faults within Alaska, dominate the seismic hazard estimates for locations near to the sources, smoothed seismicity rates make important contributions to seismic hazard away from fault-based sources and where knowledge of recurrence and magnitude is not sufficient for use in hazard studies. Recent developments in adaptive smoothing methods and statistical tests for evaluating and comparing rate models prompt us to investigate the appropriateness of adaptive smoothing for the ASHMs. We develop smoothed seismicity models for Alaska using fixed and adaptive smoothing methods and compare the resulting models by calculating and evaluating the joint likelihood test. We use the earthquake catalog, and associated completeness levels, developed for the 2007 ASHM to produce fixed-bandwidth-smoothed models with smoothing distances varying from 10 to 100 km and adaptively smoothed models. Adaptive smoothing follows the method of Helmstetter et al. and defines a unique smoothing distance for each earthquake epicenter from the distance to the nth nearest neighbor. The consequence of the adaptive smoothing methods is to reduce smoothing distances, causing locally increased seismicity rates, where seismicity rates are high and to increase smoothing distances where seismicity is sparse. We follow guidance from previous studies to optimize the neighbor number (n-value) by comparing model likelihood values, which estimate the likelihood that the observed earthquake epicenters from the recent catalog are derived from the smoothed rate models. We compare likelihood values from all rate models to rank the smoothing methods. We find that adaptively smoothed seismicity models yield better likelihood values than the fixed smoothing models. Holding all other (source and ground motion) models constant, we calculate seismic hazard curves for all points across Alaska on a 0.1 degree grid, using the adaptively smoothed and fixed smoothed seismicity models separately. Because adaptively smoothed models concentrate seismicity near the earthquake epicenters where seismicity rates are high, the corresponding hazard values are higher, locally, but reduced with distance from observed seismicity, relative to the hazard from fixed-bandwidth models. We suggest that adaptively smoothed seismicity models be considered for implementation in the update to the ASHMs because of their improved likelihood estimates relative to fixed smoothing methods; however, concomitant increases in seismic hazard will cause significant changes in regions of high seismicity, such as near the subduction zone, northeast of Kotzebue, and along the NNE trending zone of seismicity in the Alaskan interior.

Alaska↗

2014 update of the U.S. national seismic hazard maps

We held 8 regional and topical workshops across the U.S. to gather information for these maps. The maps were available to the public for comment during a 60-day period. A Steering Committee (9 experts) was assembled to review the inputs and results and provide additional insights. The maps have been presented at several professional meetings. In this talk we discuss: (1) CEUS, (2) WUS (outside California), (3) California, (4) Ground motion models We discuss the hazard maps and hazard curves used in the building code. The maps in this presentation are for 2% probability of exceedance in 50 years on a uniform firm rock site condition (V S30 =760 m/s). Maps are for 0.2 s (5 Hz) spectral acceleration with 5% damping. We consider difference maps (2014-2008) and ratio maps (2014/2008).

Book↗

Aftershock risks such as those demonstrated by the recent events in New Zealand and Japan

Recent earthquakes in New Zealand and Japan show that it is important to consider the spatial and temporal distribution of aftershocks following large magnitude events since the probability of high intensity ground motions from aftershocks, which are capable of causing significant societal impact, can be considerable. This is due to the fact that a mainshock will have many aftershocks, some of which may occur closer to populated areas and may be large enough to cause damage. When a large magnitude event strikes a region, the chance that aftershocks will cause damage can be significant as was observed after the 2011 Tohoku and 2010 Canterbury earthquakes (e.g., damage caused by M w 6.6 April 11, 2011 Fukushima-Hamadori earthquake following Tohoku earthquake or by M w 6.3 February 22, 2011 Christchurch earthquake following Canterbury earthquake). Aftershock events may further damage already damaged buildings, thereby further complicating assessments of risk to the built environment. In this paper, the issue of aftershock risk is addressed by summarizing current research regarding: (1) aftershock hazard, (2) structural fragility/vulnerability before and after the mainshock, and (3) change in risk due to aftershocks.

Conference Paper↗

Toward a time-dependent probabilistic seismic hazard analysis for Alaska

We report on a time-dependent seismic hazard analysis for Alaska and the Aleutians to complement our recently completed time-independent map. Whereas the time-independent map treats all sources as statistically independent, the time-dependent analysis is based on calculations of the conditional probability of occurrence for the next 50 years by using a Brownian Passage Time model for the seismic sources judged to be characteristic. We then consider how those probabilities are modified by coseismic and postseismic stress changes resulting from large regional earthquakes occurring from 1938 to 2002. Recombining the time-dependent probabilities with time-independent truncated Gutenberg–Richter and smoothed seismicity sources leads to our time-dependent probabilistic seismic hazard results. We find that when accounting for time dependence without stress changes, earthquake probabilities can be significantly altered, reducing probabilities to near zero or increasing them to several times the time-independent values. In addition, accounting for coseismic stress changes tends to have a local influence on earthquake probabilities, whereas postseismic effects can be far-reaching in both time and space. In sum, however, since we combine time-dependent and time-independent sources, the modification to seismic hazard is relatively minor, increasing or decreasing hazard adjacent to characteristic faults by about 10%. Most cities, located far from characteristic faults, are little affected.

Alaska↗

Seismic velocity model of the central United States (Version 1): Description and simulation of the 18 April 2008 Mt. Carmel, Illinois, Earthquake

We have developed a new three‐dimensional seismic velocity model of the central United States (CUSVM) that includes the New Madrid Seismic Zone (NMSZ) and covers parts of Arkansas, Mississippi, Alabama, Illinois, Missouri, Kentucky, and Tennessee. The model represents a compilation of decades of crustal research consisting of seismic, aeromagnetic, and gravity profiles; geologic mapping; geophysical and geological borehole logs; and inversions of the regional seismic properties. The density, P ‐ and S ‐wave velocities are synthesized in a stand‐alone spatial database that can be queried to generate the required input for numerical seismic‐wave propagation simulations. We test and calibrate the CUSVM by simulating ground motions of the 18 April 2008 M w 5.4 Mt. Carmel, Illinois, earthquake and comparing the results with observed records within the model area. The selected stations in the comparisons reflect different geological site conditions and cover distances ranging from 10 to 430 km from the epicenter. The results, based on a qualitative and quantitative goodness‐of‐fit (GOF) characterization, indicate that both within and outside the Mississippi Embayment the CUSVM reasonably reproduces: (1) the body and surface‐wave arrival times and (2) the observed regional variations in ground‐motion amplitude, cumulative energy, duration, and frequency content up to a frequency of 1.0 Hz. In addition, we discuss the probable structural causes for the ground‐motion patterns in the central United States that we observed in the recorded motions of the 18 April Mt. Carmel earthquake.

Illinois↗

Including foreshocks and aftershocks in time-independent probabilistic seismic hazard analyses

Time‐independent probabilistic seismic‐hazard analysis treats each source as being temporally and spatially independent; hence foreshocks and aftershocks, which are both spatially and temporally dependent on the mainshock, are removed from earthquake catalogs. Yet, intuitively, these earthquakes should be considered part of the seismic hazard, capable of producing damaging ground motions. In this study, I consider the mainshock and its dependents as a time‐independent cluster, each cluster being temporally and spatially independent from any other. The cluster has a recurrence time of the mainshock; and, by considering the earthquakes in the cluster as a union of events, dependent events have an opportunity to contribute to seismic ground motions and hazard. Based on the methods of the U.S. Geological Survey for a high‐hazard site, the inclusion of dependent events causes ground motions that are exceeded at probability levels of engineering interest to increase by about 10% but could be as high as 20% if variations in aftershock productivity can be accounted for reliably.

Bulletin of the Seismological Society of America↗

Understanding earthquake hazards in urban areas - Evansville Area Earthquake Hazards Mapping Project

The region surrounding Evansville, Indiana, has experienced minor damage from earthquakes several times in the past 200 years. Because of this history and the proximity of Evansville to the Wabash Valley and New Madrid seismic zones, there is concern among nearby communities about hazards from earthquakes. Earthquakes currently cannot be predicted, but scientists can estimate how strongly the ground is likely to shake as a result of an earthquake and are able to design structures to withstand this estimated ground shaking. Earthquake-hazard maps provide one way of conveying such information and can help the region of Evansville prepare for future earthquakes and reduce earthquake-caused loss of life and financial and structural loss. The Evansville Area Earthquake Hazards Mapping Project (EAEHMP) has produced three types of hazard maps for the Evansville area: (1) probabilistic seismic-hazard maps show the ground motion that is expected to be exceeded with a given probability within a given period of time; (2) scenario ground-shaking maps show the expected shaking from two specific scenario earthquakes; (3) liquefaction-potential maps show how likely the strong ground shaking from the scenario earthquakes is to produce liquefaction. These maps complement the U.S. Geological Survey's National Seismic Hazard Maps but are more detailed regionally and take into account surficial geology, soil thickness, and soil stiffness; these elements greatly affect ground shaking.

Indiana↗

Liquefaction hazard for the region of Evansville, Indiana

We calculated liquefaction potential index for a grid of sites in the Evansville, Indiana area for two scenario earthquakes-a magnitude 7.7 in the New Madrid seismic zone and a M6.8 in the Wabash Valley seismic zone. For the latter event, peak ground accelerations range from 0.13 gravity to 0.81 gravity, sufficiently high to be of concern for liquefaction. Recently acquired cone-penetrometer test data at 58 sites were used to estimate the factor of safety against liquefaction and liquefaction potential index at each site. To extend the estimation of liquefaction hazard to a grid of sites in the area, the soil columns at these grid sites were divided into three categories, and for each category a sufficient number of cone-penetrometer test sites were available to characterize statistically each group's cone-penetrometer test tip resistance and sleeve friction. At each grid site, Monte Carlo sampling was used to generate values for these two parameters at 2-meter intervals for depths down to 20 meters or bedrock. The groundwater table at each grid site was likewise sampled from a mean value and group-dependent standard deviation. For each grid site, 25,000 realizations of the soil profile were generated and a probability distribution of liquefaction potential index values was obtained. Maps of liquefaction hazard for each scenario earthquake present (1) Mean liquefaction potential index at each site, and (2) Probabilities that liquefaction potential index values exceed 5 (threshold for expression of surface liquefaction) and 12 (threshold for lateral spreading). Values for the liquefaction potential index are high in the River alluvium group, where the soil profiles are predominantly sand, while values in the Lacustrine terrace group are lower, owing to the predominance of clay. Liquefaction potential index values in the Outwash terrace group are less consistent because the soil profiles contain highly variable sequences of silty sand, clayey sand, and sandy clay, justifying the use of the Monte Carlo procedure to capture the consequences of this complexity.

Evansville, Indiana↗

Probabilistic seismic hazard assessment including site effects for Evansville, Indiana, and the surrounding region

We provide a probabilistic seismic hazard assessment for the Evansville, Indiana region incorporating information from new surficial geologic mapping efforts on the part of the U.S. Geological Survey (USGS) and the Kentucky and Indiana State Geological Surveys, as well as information on the thickness and properties of near surface soils and their associated uncertainties. The subsurface information has been compiled to determine bedrock elevation and reference depth-dependent shear-wave velocity models for the different soil types. The probabilistic seismic hazard calculation applied here follows the method used for the 2008 U.S. Geological Survey National Seismic Hazard Maps, with modifications to incorporate estimates of local site conditions and their uncertainties, in a completely probabilistic manner. The resulting analysis shows strong local variations of acceleration with 2 percent probability of exceedance in 50 years, particularly for 0.2-second (s) period spectral acceleration (SA), that are clearly correlated with variations in the thickness of unconsolidated soils above bedrock. These values are much greater than the USGS national seismic hazard map values, which assume B/C site conditions. When compared to the national maps with an assumed uniform site D class amplification factor applied, the high-resolution seismic hazard maps have higher amplitudes for peak ground acceleration and 0.2-s SA for most of the map region. However, deamplification relative to the D class national seismic hazard maps appears to play an important role within the limits of the ancient bedrock valley underlying Evansville where soils are thickest. For 1.0-s SA, the new high-resolution seismic hazard maps show levels consistent with D class site response within the limits of this ancient bedrock valley, but levels consistent with B/C site conditions outside.

Open-File Report↗

Earthquake scenario ground motions for the urban area of Evansville, Indiana

The Wabash Valley seismic zone and the New Madrid seismic zone are the closest large earthquake source zones to Evansville, Indiana. The New Madrid earthquakes of 1811-1812, over 180 kilometers (km) from Evansville, produced ground motions with a Modified Mercalli Intensity of VII near Evansville, the highest intensity observed in Indiana. Liquefaction evidence has been documented less than 40 km away from Evansville resulting from two large earthquakes in the past 12,000 years in the Wabash Valley. Two earthquake scenarios are described in this paper that demonstrate the expected ground motions for a 33×42-km region around Evansville based on a repeat earthquake from each of these source regions. We perform a one-dimensional analysis for a grid of sites that takes into account the amplification or deamplification of ground motion in the unconsolidated soil layer using a new three-dimensional model of seismic velocity and bedrock depth. There are significant differences in the calculated amplification from that expected for National Earthquake Hazard Reduction Program site class D conditions, with deamplification at many locations within the ancient bedrock valley underlying Evansville. Ground motions relative to the acceleration of gravity (g) in the Evansville area from a simulation of a magnitude (M) 7.7 New Madrid earthquake range from 0.15 to 0.25 g for peak ground acceleration, 0.14 to 0.7 g for 0.2-second (s) spectral acceleration, and 0.05 to 0.25 g for 1.0-s spectral acceleration. Ground motions from a M6.8 Wabash Valley earthquake centered 40 km northwest of the city produce ground motions that decrease with distance from 1.5 to 0.3 g for 0.2-s spectral acceleration when they reach the main part of Evansville, but then increase in amplitude from 0.3 to 0.6 g south of the city and the Ohio River. The densest urbanization in Evansville and Henderson, Ky., is within the area of preferential amplification at 1.0-s period for both scenarios, but the area experiences relatively less amplification than surrounding areas at 0.2 s, consistent with expected resonance periods based on the soil profiles.

Indiana↗

Seismic-Hazard Maps for the Conterminous United States, 2008

Probabilistic seismic-hazard maps were prepared for the conterminous United States portraying peak horizontal acceleration and horizontal spectral response acceleration for 0.2- and 1.0-second periods with probabilities of exceedance of 10 percent in 50 years and 2 percent in 50 years. All of the maps were prepared by combining the hazard derived from spatially smoothed historic seismicity with the hazard from fault-specific sources. The acceleration values contoured are the random horizontal component. The reference site condition is firm rock, defined as having an average shear-wave velocity of 760 m/s in the top 30 meters corresponding to the boundary between NEHRP (National Earthquake Hazards Reduction program) site classes B and C. This data set represents the results of calculations of hazard curves for a grid of points with a spacing of 0.05 degrees in latitude and longitude. The grid of points were contoured to produce the final representation of the seismic-hazard. These maps are intended to summarize the available quantitative information about seismic ground motion hazard for the conterminous United States from geologic and geophysical source.

Scientific Investigations Map↗

Evansville Area Earthquake Hazards Mapping Project (EAEHMP)— Progress report, 2008

Maps of surficial geology, deterministic and probabilistic seismic hazard, and liquefaction potential index have been prepared by various members of the Evansville Area Earthquake Hazard Mapping Project for seven quadrangles in the Evansville, Indiana, and Henderson, Kentucky, metropolitan areas. The surficial geologic maps feature 23 types of surficial geologic deposits, artificial fill, and undifferentiated bedrock outcrop and include alluvial and lake deposits of the Ohio River valley. Probabilistic and deterministic seismic hazard and liquefaction hazard mapping is made possible by drawing on a wealth of information including surficial geologic maps, water well logs, and in-situ testing profiles using the cone penetration test, standard penetration test, down-hole shear wave velocity tests, and seismic refraction tests. These data were compiled and collected with contributions from the Indiana Geological Survey, Kentucky Geological Survey, Illinois State Geological Survey, United States Geological Survey, and Purdue University. Hazard map products are in progress and are expected to be completed by the end of 2009, with a public roll out in early 2010. Preliminary results suggest that there is a 2 percent probability that peak ground accelerations of about 0.3 g will be exceeded in much of the study area within 50 years, which is similar to the 2002 USGS National Seismic Hazard Maps for a firm rock site value. Accelerations as high as 0.4-0.5 g may be exceeded along the edge of the Ohio River basin. Most of the region outside of the river basin has a low liquefaction potential index (LPI), where the probability that LPI is greater than 5 (that is, there is a high potential for liquefaction) for a M7.7 New Madrid type event is only 20-30 percent. Within the river basin, most of the region has high LPI, where the probability that LPI is greater than 5 for a New Madrid type event is 80-100 percent.

Indiana, Kentucky↗

Documentation for the 2008 update of the United States National Seismic Hazard Maps

The 2008 U.S. Geological Survey (USGS) National Seismic Hazard Maps display earthquake ground motions for various probability levels across the United States and are applied in seismic provisions of building codes, insurance rate structures, risk assessments, and other public policy. This update of the maps incorporates new findings on earthquake ground shaking, faults, seismicity, and geodesy. The resulting maps are derived from seismic hazard curves calculated on a grid of sites across the United States that describe the frequency of exceeding a set of ground motions. The USGS National Seismic Hazard Mapping Project developed these maps by incorporating information on potential earthquakes and associated ground shaking obtained from interaction in science and engineering workshops involving hundreds of participants, review by several science organizations and State surveys, and advice from two expert panels. The new probabilistic hazard maps represent an update of the 2002 seismic hazard maps developed by Frankel and others (2002), which used the methodology developed for the 1996 version of the maps (Frankel and others, 1996). Algermissen and Perkins (1976) published the first probabilistic seismic hazard map of the United States which was updated in Algermissen and others (1990). The National Seismic Hazard Maps represent our assessment of the “best available science” in earthquake hazards estimation for the United States (maps of Alaska and Hawaii as well as further information on hazard across the United States are available on our Web site at http://earthquake.usgs.gov/research/hazmaps/ ).

Open-File Report↗

Revision of Time-Independent Probabilistic Seismic Hazard Maps for Alaska

We present here time-independent probabilistic seismic hazard maps of Alaska and the Aleutians for peak ground acceleration (PGA) and 0.1, 0.2, 0.3, 0.5, 1.0 and 2.0 second spectral acceleration at probability levels of 2 percent in 50 years (annual probability of 0.000404), 5 percent in 50 years (annual probability of 0.001026) and 10 percent in 50 years (annual probability of 0.0021). These maps represent a revision of existing maps based on newly obtained data and assumptions reflecting best current judgments about methodology and approach. These maps have been prepared following the procedures and assumptions made in the preparation of the 2002 National Seismic Hazard Maps for the lower 48 States. A significant improvement relative to the 2002 methodology is the ability to include variable slip rate along a fault where appropriate. These maps incorporate new data, the responses to comments received at workshops held in Fairbanks and Anchorage, Alaska, in May, 2005, and comments received after draft maps were posted on the National Seismic Hazard Mapping Web Site. These maps will be proposed for adoption in future revisions to the International Building Code. In this documentation we describe the maps and in particular explain and justify changes that have been made relative to the 1999 maps. We are also preparing a series of experimental maps of time-dependent hazard that will be described in future documents.

Open-File Report↗

Preliminary Earthquake Hazard Map of Afghanistan

Introduction Earthquakes represent a serious threat to the people and institutions of Afghanistan. As part of a United States Agency for International Development (USAID) effort to assess the resource potential and seismic hazards of Afghanistan, the Seismic Hazard Mapping group of the United States Geological Survey (USGS) has prepared a series of probabilistic seismic hazard maps that help quantify the expected frequency and strength of ground shaking nationwide. To construct the maps, we do a complete hazard analysis for each of ~35,000 sites in the study area. We use a probabilistic methodology that accounts for all potential seismic sources and their rates of earthquake activity, and we incorporate modeling uncertainty by using logic trees for source and ground-motion parameters. See the Appendix for an explanation of probabilistic seismic hazard analysis and discussion of seismic risk. Afghanistan occupies a southward-projecting, relatively stable promontory of the Eurasian tectonic plate (Ambraseys and Bilham, 2003; Wheeler and others, 2005). Active plate boundaries, however, surround Afghanistan on the west, south, and east. To the west, the Arabian plate moves northward relative to Eurasia at about 3 cm/yr. The active plate boundary trends northwestward through the Zagros region of southwestern Iran. Deformation is accommodated throughout the territory of Iran; major structures include several north-south-trending, right-lateral strike-slip fault systems in the east and, farther to the north, a series of east-west-trending reverse- and strike-slip faults. This deformation apparently does not cross the border into relatively stable western Afghanistan. In the east, the Indian plate moves northward relative to Eurasia at a rate of about 4 cm/yr. A broad, transpressional plate-boundary zone extends into eastern Afghanistan, trending southwestward from the Hindu Kush in northeast Afghanistan, through Kabul, and along the Afghanistan-Pakistan border. Deformation here is expressed as a belt of major, north-northeast-trending, left-lateral strike-slip faults and abundant seismicity. The seismicity intensifies farther to the northeast and includes a prominent zone of deep earthquakes associated with northward subduction of the Indian plate beneath Eurasia that extends beneath the Hindu Kush and Pamirs Mountains. Production of the seismic hazard maps is challenging because the geological and seismological data required to produce a seismic hazard model are limited. The data that are available for this project include historical seismicity and poorly constrained slip rates on only a few of the many active faults in the country. Much of the hazard is derived from a new catalog of historical earthquakes: from 1964 to the present, with magnitude equal to or greater than about 4.5, and with depth between 0 and 250 kilometers. We also include four specific faults in the model: the Chaman fault with an assigned slip rate of 10 mm/yr, the Central Badakhshan fault with an assigned slip rate of 12 mm/yr, the Darvaz fault with an assigned slip rate of 7 mm/yr, and the Hari Rud fault with an assigned slip rate of 2 mm/yr. For these faults and for shallow seismicity less than 50 km deep, we incorporate published ground-motion estimates from tectonically active regions of western North America, Europe, and the Middle East. Ground-motion estimates for deeper seismicity are derived from data in subduction environments. We apply estimates derived for tectonic regions where subduction is the main tectonic process for intermediate-depth seismicity between 50- and 250-km depth. Within the framework of these limitations, we have developed a preliminary probabilistic seismic-hazard assessment of Afghanistan, the type of analysis that underpins the seismic components of modern building codes in the United States. The assessment includes maps of estimated peak ground-acceleration (PGA), 0.2-second spectral acceleration (SA), and 1.0-secon

Open-File Report↗