USGS Science⌕ Search

Geology topics

R. E. Kayen

Publications and source records attributed to R. E. Kayen.

At least 19 recordsLinked to original sources

Spatial correlation of shear-wave velocity within San Francisco Bay Sediments

Sediment properties are spatially variable at all scales, and this variability at smaller scales influences high frequency ground motions. We show that surface shear-wave velocity is highly correlated within San Francisco Bay Area sediments using shear-wave velocity measurements from 210 seismic cone penetration tests. We use this correlation to estimate the surface sediment velocity structure using geostatistics. We find that the variance of the estimated shear-wave velocity is reduced using ordinary kriging, and that including this velocity structure in 2D ground motion simulations of a moderate sized earthquake improves the accuracy of the synthetics. Copyright ASCE 2006.

Conference Paper↗

A case study of alternative site response explanatory variables in Parkfield, California

The combination of densely-spaced strong-motion stations in Parkfield, California, and spectral analysis of surface waves (SASW) profiles provides an ideal dataset for assessing the accuracy of different site response explanatory variables. We judge accuracy in terms of spatial coverage and correlation with observations. The performance of the alternative models is period-dependent, but generally we observe that: (1) where a profile is available, the square-root-of-impedance method outperforms V S30 (average S-wave velocity to 30 m depth), and (2) where a profile is unavailable, the topographic-slope method outperforms surficial geology. The fundamental site frequency is a valuable site response explanatory variable, though less valuable than V S30 . However, given the expense and difficulty of obtaining reliable estimates of V S30 and the relative ease with which the fundamental site frequency can be computed, the fundamental site frequency may prove to be a valuable site response explanatory variable for many applications.

California↗

Multiscale site-response mapping: A case study of Parkfield, California

The scale of previously proposed methods for mapping site-response ranges from global coverage down to individual urban regions. Typically, spatial coverage and accuracy are inversely related.We use the densely spaced strong-motion stations in Parkfield, California, to estimate the accuracy of different site-response mapping methods and demonstrate a method for integrating multiple site-response estimates from the site to the global scale. This method is simply a weighted mean of a suite of different estimates, where the weights are the inverse of the variance of the individual estimates. Thus, the dominant site-response model varies in space as a function of the accuracy of the different models. For mapping applications, site-response models should be judged in terms of both spatial coverage and the degree of correlation with observed amplifications. Performance varies with period, but in general the Parkfield data show that: (1) where a velocity profile is available, the square-rootof- impedance (SRI) method outperforms the measured VS30 (30 m divided by the S-wave travel time to 30 m depth) and (2) where velocity profiles are unavailable, the topographic slope method outperforms surficial geology for short periods, but geology outperforms slope at longer periods. We develop new equations to estimate site response from topographic slope, derived from the Next Generation Attenuation (NGA) database.

Bulletin of the Seismological Society of America↗

Retesting of liquefaction and nonliquefaction case histories from the 1976 Tangshan earthquake

A field investigation was performed to retest liquefaction and nonliquefaction sites from the 1976 Tangshan earthquake in China. These sites were carefully investigated in 1978 and 1979 by using standard penetration test (SPT) and cone penetration test (CPT) equipment; however, the CPT measurements are obsolete because of the now nonstandard cone that was used at the time. In 2007, a modern cone was mobilized to retest 18 selected sites that are particularly important because of the intense ground shaking they sustained despite their high fines content and/or because the site did not liquefy. Of the sites reinvestigated and carefully reprocessed, 13 were considered accurate representative case histories. Two of the sites that were originally investigated for liquefaction have been reinvestigated for cyclic failure of fine-grained soil and removed from consideration for liquefaction triggering. The most important outcome of these field investigations was the collection of more accurate data for three nonliquefaction sites that experienced intense ground shaking. Data for these three case histories is now included in an area of the liquefaction triggering database that was poorly populated and will help constrain the upper bound of future liquefaction triggering curves.

Tangshan↗

A geostatistical approach to mapping site response spectral amplifications

If quantitative estimates of the seismic properties do not exist at a location of interest then the site response spectral amplifications must be estimated from data collected at other locations. Currently, the most common approach employs correlations of site class with maps of surficial geology. Analogously, correlations of site class with topographic slope can be employed where the surficial geology is unknown. Our goal is to identify and validate a method to estimate site response with greater spatial resolution and accuracy for regions where additional effort is warranted. This method consists of three components: region-specific data collection, a spatial model for interpolating seismic properties, and a theoretical method for computing spectral amplifications from the interpolated seismic properties. We consider three spatial interpolation schemes: correlations with surficial geology, termed the geologic trend (GT), ordinary kriging (OK), and kriging with a trend (KT). We estimate the spectral amplifications from seismic properties using the square root of impedance method, thereby linking the frequency-dependent spectral amplifications to the depth-dependent seismic properties. Thus, the range of periods for which this method is applicable is limited by the depth of exploration. A dense survey of near-surface S-wave slowness (Ss) throughout Kobe, Japan shows that the geostatistical methods give more accurate estimates of Ss than the topographic slope and GT methods, and the OK and KT methods perform equally well. We prefer the KT model because it can be seamlessly integrated with geologic maps that cover larger regions. Empirical spectral amplifications show that the region-specific data achieve more accurate estimates of observed median short-period amplifications than the topographic slope method. ?? 2010 Elsevier B.V.

Engineering Geology↗

Impediments to predicting site response: Seismic property estimation and modeling simplifications

We compare estimates of the empirical transfer function (ETF) to the plane SH-wave theoretical transfer function (TTF) within a laterally constant medium for invasive and noninvasive estimates of the seismic shear-wave slownesses at 13 Kiban-Kyoshin network stations throughout Japan. The difference between the ETF and either of the TTFs is substantially larger than the difference between the two TTFs computed from different estimates of the seismic properties. We show that the plane SH-wave TTF through a laterally homogeneous medium at vertical incidence inadequately models observed amplifications at most sites for both slowness estimates, obtained via downhole measurements and the spectral analysis of surface waves. Strategies to improve the predictions can be separated into two broad categories: improving the measurement of soil properties and improving the theory that maps the 1D soil profile onto spectral amplification. Using an example site where the 1D plane SH-wave formulation poorly predicts the ETF, we find a more satisfactory fit to the ETF by modeling the full wavefield and incorporating spatially correlated variability of the seismic properties. We conclude that our ability to model the observed site response transfer function is limited largely by the assumptions of the theoretical formulation rather than the uncertainty of the soil property estimates.

Bulletin of the Seismological Society of America↗

Database for earthquake strong motion studies in Italy

We describe an Italian database of strong ground motion recordings and databanks delineating conditions at the instrument sites and characteristics of the seismic sources. The strong motion database consists of 247 corrected recordings from 89 earthquakes and 101 recording stations. Uncorrected recordings were drawn from public web sites and processed on a record-by-record basis using a procedure utilized in the Next-Generation Attenuation (NGA) project to remove instrument resonances, minimize noise effects through low- and high-pass filtering, and baseline correction. The number of available uncorrected recordings was reduced by 52% (mostly because of s-triggers) to arrive at the 247 recordings in the database. The site databank includes for every recording site the surface geology, a measurement or estimate of average shear wave velocity in the upper 30 m (V s30 ), and information on instrument housing. Of the 89 sites, 39 have on-site velocity measurements (17 of which were performed as part of this study using SASW techniques). For remaining sites, we estimate V s30 based on measurements on similar geologic conditions where available. Where no local velocity measurements are available, correlations with surface geology are used. Source parameters are drawn from databanks maintained (and recently updated) by Istituto Nazionale di Geofisica e Vulcanologia and include hypocenter location and magnitude for small events (M< ??? 5.5) and finite source parameters for larger events. ?? 2009 A.S. Elnashai & N.N. Ambraseys.

Journal of Earthquake Engineering↗

Use of airborne and terrestrial lidar to detect ground displacement hazards to water systems

We investigate the use of multiepoch airborne and terrestrial lidar to detect and measure ground displacements of sufficient magnitude to damage buried pipelines and other water system facilities that might result, for example, from earthquake or rainfall-induced landslides. Lidar scans are performed at three sites with coincident measurements by total station surveying. Relative horizontal accuracy is evaluated by measurements of lateral dimensions of well defined objects such as buildings and tanks; we find misfits ranging from approximately 5 to 12 cm, which is consistent with previous work. The bias and dispersion of lidar elevation measurements, relative to total station surveying, is assessed at two sites: (1) a power plant site (PP2) with vegetated steeply sloping terrain; and (2) a relatively flat and unvegetated site before and after trenching operations were performed. At PP2, airborne lidar showed minimal elevation bias and a standard deviation of approximately 70 cm, whereas terrestrial lidar did not produce useful results due to beam divergence issues and inadequate sampling of the study region. At the trench site, airborne lidar showed minimal elevation bias and reduced standard deviation relative to PP2 (6-20 cm), whereas terrestrial lidar was nearly unbiased with very low dispersion (4-6 cm). Pre- and posttrench bias-adjusted normalized residuals showed minimal to negligible correlation, but elevation change was affected by relative bias between epochs. The mean of elevation change bias essentially matches the difference in means of pre- and posttrench elevation bias, whereas elevation change standard deviation is sensitive to the dispersion of individual epoch elevations and their correlation coefficient. The observed lidar bias and standard deviations enable reliable detection of damaging ground displacements for some pipelines types (e.g., welded steel) but not all (e.g., concrete with unwelded, mortared joints). ?? ASCE 2009.

Journal of Surveying Engineering↗

New Orleans and Hurricane Katrina. IV: Orleans East Bank (Metro) protected basin

This paper addresses damage caused by Hurricane Katrina to the main Orleans East Bank protected basin. This basin represented the heart of New Orleans, and contained the main downtown area, the historic French Quarter, the Garden District, and the sprawling Lakefront and Canal Districts. Nearly half of the loss of life during this hurricane, and a similar fraction of the overall damages, occurred in this heavily populated basin. There are a number of important geotechnical lessons, as well as geo-forensic lessons, associated with the flooding of this basin. These include the difficulties associated with the creation and operation of regional-scale flood protection systems requiring federal and local cooperation and funding over prolonged periods of time. There are also a number of engineering and policy lessons regarding (1) the accuracy and reliability of current analytical methods; (2) the shortcomings and potential dangers involved in decisions that reduced short-term capital outlays in exchange for increased risk of potential system failures; (3) the difficulties associated with integrating local issues with a flood risk reduction project; and (4) the need to design and maintain levees as systems; with each of the many individual project elements being required to mesh seamlessly. These lessons are of interest and importance for similar flood protection systems throughout numerous other regions of the United States and the world. ?? 2008 ACSE.

Journal of Geotechnical and Geoenvironmental Engin↗

Submarine slope failures near Seward, Alaska, during the M9.2 1964 earthquake

Following the 1964 M9.2 megathrust earthquake in southern Alaska, Seward was the only town hit by tsunamis generated from both submarine landslides and tectonic sources. Within 45 seconds of the start of the earthquake, a 1.2-km-long section of waterfront began sliding seaward, and soon after, ~6-8-m high waves inundated the town. Studies soon after the earthquake concluded that submarine landslides along the Seward waterfront generated the tsunamis that occurred immediately after the earthquake. We analyze pre- and post-earthquake bathymetry data to assess the location and extent of submarine mass failures and sediment transport. New NOAA multibeam bathymetry shows the morphology of the entire fjord at 15 m resolution. We also assembled all older soundings from smooth sheets for comparison to the multibeam dataset. We gridded the sounding data, applied corrections for coseismic subsidence, post-seismic rebound, unrecovered co-seismic subsidence, sea-level rise (vertical datum shift), and measurement errors. The difference grids show changes resulting from the 1964 earthquake. We estimate the total volume of slide material to be about 211 million m3. Most of this material was transported to a deep, flat area, which we refer to as &ldquo;the bathtub&rdquo;, about 6 to 13 km south of Seward. Sub-bottom profiling of the bathtub shows an acoustically transparent unit, which we interpret as a sediment flow deposit resulting from the submarine landslides. The scale of the submarine landslides and the distance over which sediment was transported is much larger than previously appreciated.

Alaska↗

Reassessment of seismically induced, tsunamigenic submarine slope failures in Port Valdez, Alaska, USA

The M9.2 Alaska earthquake of 1964 caused major damage to the port facilities and town of Valdez, most of it the result of submarine landslide and the consequent tsunamis. Recent bathymetric multibeam surveys, high-resolution subbottom profiles, and dated sediment cores in Port Valdez supply new information about the morphology and character of the landslide deposits. A comparison of pre- and post-earthquake bathymetry provides an estimate of the net volume of landslide debris deposited in the basin and the volume of sediment removed from the source region. Landslide features include (1) large blocks (up to 40-m high) near the location of the greatest tsunamiwave runup (~50 m), (2) two debris lobes associated with the blocks, (3) a series of gullies, channels and talus, near the fjord-head delta and badly damaged old town of Valdez, and (4) the front of a debris lobe that flowed half-way down the fjord from the east end.

Alaska↗

Spatial correlation of shear-wave velocity in the San Francisco Bay Area sediments

Ground motions recorded within sedimentary basins are variable over short distances. One important cause of the variability is that local soil properties are variable at all scales. Regional hazard maps developed for predicting site effects are generally derived from maps of surficial geology; however, recent studies have shown that mapped geologic units do not correlate well with the average shear-wave velocity of the upper 30 m, V s (30). We model the horizontal variability of near-surface soil shear-wave velocity in the San Francisco Bay Area to estimate values in unsampled locations in order to account for site effects in a continuous manner. Previous geostatistical studies of soil properties have shown horizontal correlations at the scale of meters to tens of meters while the vertical correlations are on the order of centimeters. In this paper we analyze shear-wave velocity data over regional distances and find that surface shear-wave velocity is correlated at horizontal distances up to 4 km based on data from seismic cone penetration tests and the spectral analysis of surface waves. We propose a method to map site effects by using geostatistical methods based on the shear-wave velocity correlation structure within a sedimentary basin. If used in conjunction with densely spaced shear-wave velocity profiles in regions of high seismic risk, geostatistical methods can produce reliable continuous maps of site effects.

California↗

CPT-based probabilistic and deterministic assessment of in situ seismic soil liquefaction potential

This paper presents a complete methodology for both probabilistic and deterministic assessment of seismic soil liquefaction triggering potential based on the cone penetration test (CPT). A comprehensive worldwide set of CPT-based liquefaction field case histories were compiled and back analyzed, and the data then used to develop probabilistic triggering correlations. Issues investigated in this study include improved normalization of CPT resistance measurements for the influence of effective overburden stress, and adjustment to CPT tip resistance for the potential influence of "thin" liquefiable layers. The effects of soil type and soil character (i.e., "fines" adjustment) for the new correlations are based on a combination of CPT tip and sleeve resistance. To quantify probability for performancebased engineering applications, Bayesian "regression" methods were used, and the uncertainties of all variables comprising both the seismic demand and the liquefaction resistance were estimated and included in the analysis. The resulting correlations were developed using a Bayesian framework and are presented in both probabilistic and deterministic formats. The results are compared to previous probabilistic and deterministic correlations. ?? 2006 ASCE.

Journal of Geotechnical and Geoenvironmental Engin↗

Varieties of submarine failure morphologies of seismically-induced landslides in Alaskan fjords

The Great Alaska Earthquake of 1964 caused major damage and 43 deaths in the coastal communities of Seward and Valdez. Most of these losses were caused by tsunamis that occurred immediately after the earthquake and were most likely induced by local submarine landslides. Recent NOAA multibcam bathymetric surveys near Seward and Valdez provide detailed information about the morphology of landslide deposits in both areas. High-resolution (chirp) surveys were conducted by the USGS, and sediment samples were taken over apparent landslide debris. Landslide deposits near Seward typically take the form of a series of large and small blocks lying directly off the front of the town, although there are indications of sandy and muddy debris flows occurring off river deltas. Near Valdez, landslide morphologies include at least three forms: a field of large blocks (up to 40 m high), an intricate series of gullies, channels, and talus near the fjord-head delta and a broad debris lobe that apparently flowed half-way down the fjord and stopped. The 1964 landslide tsunamis may have been composites resulting from a number of landslide events.

Conference Paper↗

Probabilistic liquefaction triggering based on the cone penetration test

Performance-based earthquake engineering requires a probabilistic treatment of potential failure modes in order to accurately quantify the overall stability of the system. This paper is a summary of the application portions of the probabilistic liquefaction triggering correlations proposed recently proposed by Moss and co-workers. To enable probabilistic treatment of liquefaction triggering, the variables comprising the seismic load and the liquefaction resistance were treated as inherently uncertain. Supporting data from an extensive Cone Penetration Test (CPT)-based liquefaction case history database were used to develop a probabilistic correlation. The methods used to measure the uncertainty of the load and resistance variables, how the interactions of these variables were treated using Bayesian updating, and how reliability analysis was applied to produce curves of equal probability of liquefaction are presented. The normalization for effective overburden stress, the magnitude correlated duration weighting factor, and the non-linear shear mass participation factor used are also discussed.

Conference Paper↗

Standard penetration test-based probabilistic and deterministic assessment of seismic soil liquefaction potential

This paper presents'new correlations for assessment of the likelihood of initiation (or triggering) of soil liquefaction. These new correlations eliminate several sources of bias intrinsic to previous, similar correlations, and provide greatly reduced overall uncertainty and variance. Key elements in the development of these new correlations are (1) accumulation of a significantly expanded database of field performance case histories; (2) use of improved knowledge and understanding of factors affecting interpretation of standard penetration test data; (3) incorporation of improved understanding of factors affecting site-specific earthquake ground motions (including directivity effects, site-specific response, etc.); (4) use of improved methods for assessment of in situ cyclic shear stress ratio; (5) screening of field data case histories on a quality/uncertainty basis; and (6) use of high-order probabilistic tools (Bayesian updating). The resulting relationships not only provide greatly reduced uncertainty, they also help to resolve a number of corollary issues that have long been difficult and controversial including: (1) magnitude-correlated duration weighting factors, (2) adjustments for fines content, and (3) corrections for overburden stress. ?? ASCE.

Journal of Geotechnical and Geoenvironmental Engin↗

Influence of near-surface stratigraphy on coastal landslides at Sleeping Bear Dunes National Lakeshore, Lake Michigan, USA

Lake-level change and landslides are primary controls on the development of coastal environments along the coast of northeastern Lake Michigan. The late Quaternary geology of Sleeping Bear Dunes National Lakeshore was examined with high-resolution seismic reflection profiles, ground-penetrating radar (GPR), and boreholes. Based on sequence-stratigraphic principles, this study recognizes ten stratigraphic units and three major unconformities that were formed by late Pleistocene glaciation and postglacial lake-level changes. Locally high sediment supply, and reworking by two regressions and a transgression have produced a complex stratigraphy that is prone to episodic failure. In 1995, a large landslide deposited approximately 1 million m3 of sediment on the lake floor. The highly deformed landslide deposits, up to 18 m thick, extend 3-4 km offshore and unconformably overlie well-stratified glacial and lacustrine sediment. The landslide-prone bluff is underlain by channel-fill deposits that are oriented nearly perpendicular to the shoreline. The paleochannels are at least 10 m deep and 400 m wide and probably represent stream incision during a lake-level lowstand about 10.3 ka B.P. The channels filled with sediment during the subsequent transgression and lake-level highstand, which climaxed about 4.5 ka B.P. As lake level fell from the highstand, the formation of beach ridges and sand dunes sealed off the channel and isolated a small inland lake (Glen Lake), which lies 5 m above the level of Lake Michigan and may be a source of piped groundwater. Our hypothesis is that the paleochannels act as conduits for pore water flow, and thereby locally reduce soil strength and promote slope failure.

Michigan↗

Landslides and liquefaction triggered by the M 7.9 denali fault earthquake of 3 November 2002

The moment magnitude (M) 7.9 Denali Fault earthquake in Alaska of 3 November 2002 triggered an unusual pattern of landslides and liquefaction effects. The landslides were primarily rock falls and rock slides that ranged in volume from a few cubic meters to the 40 million-cubic-meter rock avalanche that covered much of the McGinnis Glacier. Landslides were concentrated in a narrow zone ???30 km wide that straddled the fault rupture zone over its entire 300 km length. Large rock avalanches all clustered at the western end of the rupture zone where acceleration levels are reported to have been the highest. Liquefaction effects, consisting of sand blows, lateral spreads, and settlement, were widespread within susceptible alluvial deposits extending from Fairbanks eastward several hundred kilometers. The liquefaction effects displayed a pattern of increasing concentration and severity from west to east and extended well beyond the zone of landslides, which is unusual. The contrasting patterns formed by the distributions of landslides and liquefaction effects initially seemed to be inconsistent; however, preliminary analyses of strong-motion records from the earthquake offer a possible explanation for the unusual ground-failure patterns that are related to three subevents that have been discerned from the earthquake records.

GSA Today↗