Geology topics

Jonathan P. Stewart

Publications and source records attributed to Jonathan P. Stewart.

29 records · Page 2Linked to original sources

USGS2018

A proposed seismic velocity profile database model

We describe the data model that we intend to use in a publicly available site profile database under development for the United States. The initial implementation of the database contains data from California. Currently, our prototype data model consists of JavaScript Object Notation (JSON) format files for storing metadata and data. For a site to be included in the database, the minimum metadata requirements are geodetic coordinates and elevation values, and the minimum data requirement is a shear-wave velocity profile. The JSON files are structured in a hierarchal manner to store metadata and data using a nested structure consisting of location, velocity profiles, dispersion curve data (for surface-wave methods), geotechnical data, and horizontal-to-vertical spectral ratios. The database schema at the current stage of the project, and as we continue to develop the data model we will consider including other relevant data, as well as evaluate other file formats to increase the efficiency of data storage and querying. In the current data model, location information includes site geodetic values (latitude, longitude, and elevation) and various site descriptors related to surface geology, geomorphic terrain category, slope gradient at various resolutions, and a geotechnical site category. Velocity data include the geophysical method(s) used to obtain the shear-wave velocity profile, type of data recorded, modeled primary- and shear-wave velocity as a function of depth, modeled profile maximum depth, and the calculated VS30 value. In the case of surface-wave based data, dispersion curve data can be recorded in data structure as phase velocity versus either wavelength or frequency. Geotechnical data includes boring logs penetration resistance, cone penetration test sounding logs, and laboratory index test results. Horizontal-to-vertical spectral ratio plots are given as a function of frequency.

Conference Paper↗

USGS2017

Development of the next generation of seismic design value maps for the 2020 NEHRP Provisions

During the period January 2015 through August 2018, a joint committee of U.S. Geological Survey (USGS) representatives and National Institute of Building Sciences Building Seismic Safety Council (BSSC) volunteers and staff formed a committee to conduct Project 17. The purpose of Project 17 was to formulate recommendations for the rules by which next-generation seismic design value maps, derived from USGS national seismic hazard models, will be developed for adoption by the 2020 National Earthquake Hazard Reduction Program Recommended Provisions for Seismic Regulations of New Buildings and Other Structures (NEHRP Provisions), ASCE 7-22 and the 2024 International Building Code. Two similar projects, Project 97 and Project 07, had been conducted in the past. Each of these projects established, for a period of approximately ten years, the rules by which design ground motion values referenced by the building codes would be developed both by USGS and by private consultants engaged in site-specific studies. Project 17 was originally commissioned in response to issues identified in adopting the 2014 edition of the USGS national seismic hazard model and the design procedures that reference them for use, including the NEHRP Provisions, building codes and referenced standards. Specific issues included: the engineering profession’s discontent with the fluctuating design values portrayed by successive map editions; discovery that the standard spectral shape referenced by the design provisions did not adequately represent ground motion amplitude and spectral character on some sites; and a change in seismologic characterization of the possible size of earthquakes originating on various faults and source zones. Project 17 was funded by the Federal Emergency Management Agency (FEMA), and supported by the USGS with some collaborating experts. An initial planning committee met throughout calendar year 2015 to identify key issues to be considered and to develop a work plan for addressing these as part of the 2020 NEHRP Provisions update cycle. The planning committee recommended an effort of approximately 30-months duration during which the USGS would develop draft design maps based on the rules proposed, to allow evaluation and refinement of the recommendations. A Project 17 Committee (P17C) was empaneled and four task subcommittees were formed, each charged with evaluating one of the key issues identified in the planning effort: Stabilizing mapped values; Definition of Acceptable Risk; Development of multi-period spectral parameter data; and, Definition of procedures for computing deterministic caps, should it be necessary to continue use of such caps in development of the maps. A fifth task subcommittee was formed in 2017 to look at ways to stabilize the seismic design category as an extended effort to stabilize mapped values. The P17C met three times per year throughout 2016, 2017 and 2018 to resolve these issues and develop recommendations for an updated technical basis and procedures to be followed in preparing next-generation seismic design value maps for inclusion in the NEHRP Provisions. The P17C documented these in the form of draft proposals for revision of the NEHRP Provisions. In August 2018, the P17C passed these recommendations to the Provisions Update Committee (PUC) for completion, development of consensus and adoption as appropriate.

Report↗

USGS2016

NGA-West2 equations for predicting vertical-component PGA, PGV, and 5%-damped PSA from shallow crustal earthquakes

We present ground motion prediction equations (GMPEs) for computing natural log means and standard deviations of vertical-component intensity measures (IMs) for shallow crustal earthquakes in active tectonic regions. The equations were derived from a global database with M 3.0–7.9 events. The functions are similar to those for our horizontal GMPEs. We derive equations for the primary M - and distance-dependence of peak acceleration, peak velocity, and 5%-damped pseudo-spectral accelerations at oscillator periods between 0.01–10 s. We observe pronounced M -dependent geometric spreading and region-dependent anelastic attenuation for high-frequency IMs. We do not observe significant region-dependence in site amplification. Aleatory uncertainty is found to decrease with increasing magnitude; within-event variability is independent of distance. Compared to our horizontal-component GMPEs, attenuation rates are broadly comparable (somewhat slower geometric spreading, faster apparent anelastic attenuation), V S 30 -scaling is reduced, nonlinear site response is much weaker, within-event variability is comparable, and between-event variability is greater.

Earthquake Spectra↗

USGS2014

Validation of the SCEC broadband platform V14.3 simulation methods using pseudo spectral acceleration data

This paper summarizes the evaluation of ground motion simulation methods implemented on the SCEC Broadband Platform (BBP), version 14.3 (as of March 2014). A seven-member panel, the authorship of this article, was formed to evaluate those methods for the prediction of pseudo-‐spectral accelerations (PSAs) of ground motion. The panel’s mandate was to evaluate the methods using tools developed through the validation exercise (Goulet et al. ,2014), and to define validation metrics for the assessment of the methods’ performance. This paper summarizes the evaluation process and conclusions from the panel. The five broadband, finite-source simulation methods on the BBP include two deterministic approaches herein referred to as CSM (Anderson, 2014) and UCSB (Crempien and Archuleta, 2014); a band-‐limited stochastic white noise method called EXSIM (Atkinson and Assatourians, 2014); and two hybrid approaches, referred to as G&P (Graves and Pitarka, 2014) and SDSU (Olsen and Takedatsu, 2014), which utilize a deterministic Green’s function approach for periods longer than 1 second and stochastic methods for periods shorter than 1 second. Two acceptance tests were defined to validate the broadband finite‐source ground methods (Goulet et al., 2014). Part A compared observed and simulated PSAs for periods from 0.01 to 10 seconds for 12 moderate to large earthquakes located in California, Japan, and the eastern US. Part B compared the median simulated PSAs to published NGA-‐West1 (Abrahamson and Silva, 2008; Boore and Atkinson, 2008; Campbell and Bozorgnia, 2008; and Chiou and Youngs, 2008) ground motion prediction equations (GMPEs) for specific magnitude and distance cases using a pass-‐fail criteria based on a defined acceptable range around the spectral shape of the GMPEs. For the initial Part A and Part B validation exercises during the summer of 2013, the software for the five methods was locked in at version 13.6 (see Maechling et al., 2014). In the spring of 2014, additional moderate events were considered for the Part A validation, and additional magnitude and distance cases were considered for the Part B validation, for the software locked in at version 14.3. Several of the simulation procedures, specifically UCSB and SDSU, changed significantly between versions 13.6 and 14.3. The CSM code was not submitted in time for the v14.3 evaluation and its detailed performance is not addressed in this paper. As described in Goulet et al. (2014) and Maechling et al. (2014), the BBP generates a variety of products, including three-‐component acceleration time series. A series of post-‐processing codes were developed to provide individual component PSAs and average median horizontal-‐component PSA (referred to as RotD50; Boore, 2010) for oscillator periods ranging from 0.01 to 10 seconds, as well as median PSA values computed using the NGA-‐West 1 GMPEs. The BBP was also configured to provide statistical analysis of simulation results relative to recordings (Part A) and GMPEs (Part B) as described further in sections below. As part of our evaluation, we reviewed documentation provided by each of the developers, which included the technical basis behind the methods and the developer’s self-‐assessments regarding the extrapolation capabilities (in terms of magnitude and distance ranges) of their methods. Two workshops were held in which methods and results were presented, and the panel was given the opportunity to question the developers and to have detailed technical discussions. A SCEC report (Dreger et al., 2013) describes the results of this review for BBP version 13.6. This paper summarizes that work and presents results for the more recent BBP 14.3 validation.

Seismological Research Letters↗

USGS2014

NGA-West2 Research Project

The NGA-West2 project is a large multidisciplinary, multi-year research program on the Next Generation Attenuation (NGA) models for shallow crustal earthquakes in active tectonic regions. The research project has been coordinated by the Pacific Earthquake Engineering Research Center (PEER), with extensive technical interactions among many individuals and organizations. NGA-West2 addresses several key issues in ground-motion seismic hazard, including updating the NGA database for a magnitude range of 3.0–7.9; updating NGA ground-motion prediction equations (GMPEs) for the “average” horizontal component; scaling response spectra for damping values other than 5%; quantifying the effects of directivity and directionality for horizontal ground motion; resolving discrepancies between the NGA and the National Earthquake Hazards Reduction Program (NEHRP) site amplification factors; analysis of epistemic uncertainty for NGA GMPEs; and developing GMPEs for vertical ground motion. This paper presents an overview of the NGA-West2 research program and its subprojects.

Earthquake Spectra↗

USGS2013

Calibration of semi-stochastic procedure for simulating high-frequency ground motions

Broadband ground motion simulation procedures typically utilize physics-based modeling at low frequencies, coupled with semi-stochastic procedures at high frequencies. The high-frequency procedure considered here combines deterministic Fourier amplitude spectra (dependent on source, path, and site models) with random phase. Previous work showed that high-frequency intensity measures from this simulation methodology attenuate faster with distance and have lower intra-event dispersion than in empirical equations. We address these issues by increasing crustal damping (Q) to reduce distance attenuation bias and by introducing random site-to-site variations to Fourier amplitudes using a lognormal standard deviation ranging from 0.45 for M w < 7 to zero for M w 8. Ground motions simulated with the updated parameterization exhibit significantly reduced distance attenuation bias and revised dispersion terms are more compatible with those from empirical models but remain lower at large distances (e.g., > 100 km).

Earthquake Spectra↗

USGS2013

Implications of the Mw9.0 Tohoku-Oki earthquake for ground motion scaling with source, path, and site parameters

The M w 9.0 Tohoku-oki Japan earthquake produced approximately 2,000 ground motion recordings. We consider 1,238 three-component accelerograms corrected with component-specific low-cut filters. The recordings have rupture distances between 44 km and 1,000 km, time-averaged shear wave velocities of V S 30 = 90 m/s to 1,900 m/s, and usable response spectral periods of 0.01 sec to >10 sec. The data support the notion that the increase of ground motions with magnitude saturates at large magnitudes. High-frequency ground motions demonstrate faster attenuation with distance in backarc than in forearc regions, which is only captured by one of the four considered ground motion prediction equations for subduction earthquakes. Recordings within 100 km of the fault are used to estimate event terms, which are generally positive (indicating model underprediction) at short periods and zero or negative (overprediction) at long periods. We find site amplification to scale minimally with V S 30 at high frequencies, in contrast with other active tectonic regions, but to scale strongly with V S 30 at low frequencies.

Earthquake Spectra↗

USGS2011

Representation of bidirectional ground motions for design spectra in building codes

The 2009 NEHRP Provisions modified the definition of horizontal ground motion from the geometric mean of spectral accelerations for two components to the peak response of a single lumped mass oscillator regardless of direction. These maximum-direction (MD) ground motions operate under the assumption that the dynamic properties of the structure (e.g., stiffness, strength) are identical in all directions. This assumption may be true for some in-plan symmetric structures, however, the response of most structures is dominated by modes of vibration along specific axes (e.g., longitudinal and transverse axes in a building), and often the dynamic properties (especially stiffness) along those axes are distinct. In order to achieve structural designs consistent with the collapse risk level given in the NEHRP documents, we argue that design spectra should be compatible with expected levels of ground motion along those principal response axes. The use of MD ground motions effectively assumes that the azimuth of maximum ground motion coincides with the directions of principal structural response. Because this is unlikely, design ground motions have lower probability of occurrence than intended, with significant societal costs. We recommend adjustments to make design ground motions compatible with target risk levels.

Earthquake Spectra↗

USGS2011

Liquefaction and other ground failures in Imperial County, California, from the April 4, 2010, El Mayor-Cucapah earthquake

The Colorado River Delta region of southern Imperial Valley, California, and Mexicali Valley, Baja California, is a tectonically dynamic area characterized by numerous active faults and frequent large seismic events. Significant earthquakes that have been accompanied by surface fault rupture and/or soil liquefaction occurred in this region in 1892 (M7.1), 1915 (M6.3; M7.1), 1930 (M5.7), 1940 (M6.9), 1950 (M5.4), 1957 (M5.2), 1968 (6.5), 1979 (6.4), 1980 (M6.1), 1981 (M5.8), and 1987 (M6.2; M6.8). Following this trend, the M7.2 El Mayor-Cucapah earthquake of April 4, 2010, ruptured approximately 120 kilometers along several known faults in Baja California. Liquefaction caused by the M7.2 El Mayor-Cucapah earthquake was widespread throughout the southern Imperial Valley but concentrated in the southwest corner of the valley, southwest of the city centers of Calexico and El Centro where ground motions were highest. Although there are few strong motion recordings in the very western part of the area, the recordings that do exist indicate that ground motions were on the order of 0.3 to 0.6g where the majority of liquefaction occurrences were found. More distant liquefaction occurrences, at Fites Road southwest of Brawley and along Rosita Canal northwest of Holtville were triggered where ground motions were about 0.2 g. Damage to roads was associated mainly with liquefaction of sandy river deposits beneath bridge approach fills, and in some cases liquefaction within the fills. Liquefaction damage to canal and drain levees was not always accompanied by vented sand, but the nature of the damage leads the authors to infer that liquefaction was involved in the majority of observed cases. Liquefaction-related damage to several public facilities - Calexico Waste Water Treatment Plant, Fig Lagoon levee system, and Sunbeam Lake Dam in particular - appears to be extensive. The cost to repair these facilities to prevent future liquefaction damage will likely be prohibitive. As such, it is likely that liquefaction will recur at these facilities during the next large earthquake in this area.

California↗

USGS2008

Site characterization of Italian strong motion recording stations

A dataset of site conditions at 101 Italian ground motion stations with recorded motions has been compiled that includes geologic characteristics and seismic velocities. Geologic characterization is derived principally from local geologic investigations by ENEL that include detailed mapping and cross sections. For sites lacking such detailed geologic characterization, the geology maps of the by Servizio Geologico d'Italia are used. Seismic velocities are extracted from the literature and the files of consulting engineers, geologists and public agencies for 33 sites. Data sources utilized include post earthquake site investigations (Friuli and Irpinia events), microzonation studies, and miscellaneous investigations performed by researchers or consulting engineers/geologists. Additional seismic velocities are measured by the authors using the controlled source spectral analysis of surface waves (SASW) method for 18 sites that recorded the 1997–1998 Umbria Marche earthquake sequence. The compiled velocity measurements provide data for 51 of the 101 sites. For the remaining sites, the average seismic velocity in the upper 30 m (V s30 ) "> (V s30 ) is estimated using a hybrid approach. For young Quaternary alluvium, V s30 "> 𝑉 s30 an existing empirical relationship for California sites by Wills and Clahan (2006) is used, which we justify by validating this relationship against Italian data. For Tertiary Limestone and Italian Mesozoic rocks, empirical estimates of V s30 "> 𝑉 s30 are developed using the available data.

Conference Paper↗

USGS2008

Shear wave structure of Umbria and Marche, Italy, strong motion seismometer sites Affected by the 1997-98 Umbria-Marche, Italy, earthquake sequence

A long sequence of earthquakes, eight with magnitudes between 5 and 6, struck the Umbria and Marche regions of central Italy between September 26, 1997 and July 1998. The earthquake swarm caused severe structural damage, particularly to masonry buildings, and resulted in the loss of twelve lives and about 150 injuries. The source of the events was a single seismogenic structure that consists of several faults with a prevailing northwest-southeast strike and crosses the Umbria-Marche border. The focal mechanism of the largest shocks indicates that the events were the product of shallow extensional normal faulting along a NE-SW extension perpendicular to the trend of the Apennines. The network of analog seismometer stations in the Umbria and Marche regions recorded motions of the main September and October 1997 events and a dense array of mobile digital stations, installed since September 29, recorded most of the swarm. The permanent national network Rete Accelerometrica Nazionale (RAN) is administered and maintained by Dipartimento delle Protezione Civile (DPC: Civil Protection Department); the temporary array was managed by Servizio Sismico Nazionale (SSN) in cooperation with small agencies and Universities. ENEA, the operator of many seismometer stations in Umbria, is the public Italian National Agency for New Technologies, Energy and the Environment. Many of the temporary and permanent stations in the Italian seismic network have little or no characterization of seismic velocities. In this study, we investigated 17 Italian sites using an active-source approach that employs low frequency harmonic waves to measure the dispersive nature of surface waves in the ground. We used the Spectral Analysis of Surface Wave (SASW) approach, coupled with an array of harmonic-wave electro-mechanical sources that are driven in-phase to excite the ground. An inversion algorithm using a non-linear least-squares best-fit method is used to compute shear wave velocities for up to 100 meters of the soil column. A draft report was published in the summer of 2008, followed by a comment period, lengthy discussions with Italian colleagues, and improved knowledge of the subsurface at the sites from soil logs. Four of the sites were reprocessed in order to correct issues with phase unwrapping of the field dispersion curves that complicated the velocity profile calculations at the lowest velocity sites. This report presents the final results from the reprocessing effort.

Open-File Report↗