Geology Reports⌕ Search

Geology topics

Yuehua Zeng

Publications and source records attributed to Yuehua Zeng.

36 records · Page 2Linked to original sources

A fault-based model for crustal deformation, fault slip-rates and off-fault strain rate in California

We invert Global Positioning System (GPS) velocity data to estimate fault slip rates in California using a fault‐based crustal deformation model with geologic constraints. The model assumes buried elastic dislocations across the region using Uniform California Earthquake Rupture Forecast Version 3 (UCERF3) fault geometries. New GPS velocity and geologic slip‐rate data were compiled by the UCERF3 deformation working group. The result of least‐squares inversion shows that the San Andreas fault slips at 19–22  mm/yr along Santa Cruz to the North Coast, 25–28  mm/yr along the central California creeping segment to the Carrizo Plain, 20–22  mm/yr along the Mojave, and 20–24  mm/yr along the Coachella to the Imperial Valley. Modeled slip rates are 7–16  mm/yr lower than the preferred geologic rates from the central California creeping section to the San Bernardino North section. For the Bartlett Springs section, fault slip rates of 7–9  mm/yr fall within the geologic bounds but are twice the preferred geologic rates. For the central and eastern Garlock, inverted slip rates of 7.5 and 4.9  mm/yr, respectively, match closely with the geologic rates. For the western Garlock, however, our result suggests a low slip rate of 1.7  mm/yr. Along the eastern California shear zone and southern Walker Lane, our model shows a cumulative slip rate of 6.2–6.9  mm/yr across its east–west transects, which is ∼1  mm/yr increase of the geologic estimates. For the off‐coast faults of central California, from Hosgri to San Gregorio, fault slips are modeled at 1–5  mm/yr, similar to the lower geologic bounds. For the off‐fault deformation, the total moment rate amounts to 0.88×10 19   N·m/yr, with fast straining regions found around the Mendocino triple junction, Transverse Ranges and Garlock fault zones, Landers and Brawley seismic zones, and farther south. The overall California moment rate is 2.76×10 19   N·m/yr, which is a 16% increase compared with the UCERF2 model.

Bulletin of the Seismological Society of America↗

Seismic hazard in the Intermountain West

The 2014 national seismic-hazard model for the conterminous United States incorporates new scientific results and important model adjustments. The current model includes updates to the historical catalog, which is spatially smoothed using both fixed-length and adaptive-length smoothing kernels. Fault-source characterization improved by adding faults, revising rates of activity, and incorporating new results from combined inversions of geologic and geodetic data. The update also includes a new suite of published ground motion models. Changes in probabilistic ground motion are generally less than 10% in most of the Intermountain West compared to the prior assessment, and ground-motion hazard in four Intermountain West cities illustrates the range and magnitude of change in the region. Seismic hazard at reference sites in Boise and Reno increased as much as 10%, whereas hazard in Salt Lake City decreased 5–6%. The largest change was in Las Vegas, where hazard increased 32–35%.

Earthquake Spectra↗

Seismic source characterization for the 2014 update of the U.S. National Seismic Hazard Model

We present the updated seismic source characterization (SSC) for the 2014 update of the National Seismic Hazard Model (NSHM) for the conterminous United States. Construction of the seismic source models employs the methodology that was developed for the 1996 NSHM but includes new and updated data, data types, source models, and source parameters that reflect the current state of knowledge of earthquake occurrence and state of practice for seismic hazard analyses. We review the SSC parameterization and describe the methods used to estimate earthquake rates, magnitudes, locations, and geometries for all seismic source models, with an emphasis on new source model components. We highlight the effects that two new model components—incorporation of slip rates from combined geodetic-geologic inversions and the incorporation of adaptively smoothed seismicity models—have on probabilistic ground motions, because these sources span multiple regions of the conterminous United States and provide important additional epistemic uncertainty for the 2014 NSHM.

Earthquake Spectra↗

The 2014 United States National Seismic Hazard Model

New seismic hazard maps have been developed for the conterminous United States using the latest data, models, and methods available for assessing earthquake hazard. The hazard models incorporate new information on earthquake rupture behavior observed in recent earthquakes; fault studies that use both geologic and geodetic strain rate data; earthquake catalogs through 2012 that include new assessments of locations and magnitudes; earthquake adaptive smoothing models that more fully account for the spatial clustering of earthquakes; and 22 ground motion models, some of which consider more than double the shaking data applied previously. Alternative input models account for larger earthquakes, more complicated ruptures, and more varied ground shaking estimates than assumed in earlier models. The ground motions, for levels applied in building codes, differ from the previous version by less than ±10% over 60% of the country, but can differ by ±50% in localized areas. The models are incorporated in insurance rates, risk assessments, and as input into the U.S. building code provisions for earthquake ground shaking.

Earthquake Spectra↗

Crustal deformation in the New Madrid seismic zone and the role of postseismic processes

Global Navigation Satellite System data across the New Madrid seismic zone (NMSZ) in the central United States over the period from 2000 through 2014 are analyzed and modeled with several deformation mechanisms including the following: (1) creep on subsurface dislocations, (2) postseismic frictional afterslip and viscoelastic relaxation from the 1811–1812 and 1450 earthquakes in the NMSZ, and (3) regional strain. In agreement with previous studies, a dislocation creeping at about 4 mm/yr between 12 and 20 km depth along the downdip extension of the Reelfoot fault reproduces the observations well. We find that a dynamic model of postseismic frictional afterslip from the 1450 and February 1812 Reelfoot fault events can explain this creep. Kinematic and dynamic models involving the Cottonwood Grove fault provide minimal predictive power. This is likely due to the smaller size of the December 1811 event on the Cottonwood Grove fault and a distribution of stations better suited to constrain localized strain across the Reelfoot fault. Regional compressive strain across the NMSZ is found to be less than 3 × 10 −9 /yr. If much of the present-day surface deformation results from afterslip, it is likely that many of the earthquakes we see today in the NMSZ are aftershocks from the 1811–1812 New Madrid earthquakes. Despite this conclusion, our results are consistent with observations and models of intraplate earthquake clustering. Given this and the recent paleoseismic history of the region, we suggest that seismic hazard is likely to remain significant.

Journal of Geophysical Research B: Solid Earth↗

Long‐term time‐dependent probabilities for the third Uniform California Earthquake Rupture Forecast (UCERF3)

The 2014 Working Group on California Earthquake Probabilities (WGCEP 2014) presents time-dependent earthquake probabilities for the third Uniform California Earthquake Rupture Forecast (UCERF3). Building on the UCERF3 time-independent model, published previously, renewal models are utilized to represent elastic-rebound-implied probabilities. A new methodology has been developed that solves applicability issues in the previous approach for un-segmented models. The new methodology also supports magnitude-dependent aperiodicity and accounts for the historic open interval on faults that lack a date-of-last-event constraint. Epistemic uncertainties are represented with a logic tree, producing 5,760 different forecasts. Results for a variety of evaluation metrics are presented, including logic-tree sensitivity analyses and comparisons to the previous model (UCERF2). For 30-year M≥6.7 probabilities, the most significant changes from UCERF2 are a threefold increase on the Calaveras fault and a threefold decrease on the San Jacinto fault. Such changes are due mostly to differences in the time-independent models (e.g., fault slip rates), with relaxation of segmentation and inclusion of multi-fault ruptures being particularly influential. In fact, some UCERF2 faults were simply too long to produce M 6.7 sized events given the segmentation assumptions in that study. Probability model differences are also influential, with the implied gains (relative to a Poisson model) being generally higher in UCERF3. Accounting for the historic open interval is one reason. Another is an effective 27% increase in the total elastic-rebound-model weight. The exact factors influencing differences between UCERF2 and UCERF3, as well as the relative importance of logic-tree branches, vary throughout the region, and depend on the evaluation metric of interest. For example, M≥6.7 probabilities may not be a good proxy for other hazard or loss measures. This sensitivity, coupled with the approximate nature of the model and known limitations, means the applicability of UCERF3 should be evaluated on a case-by-case basis.

Bulletin of the Seismological Society of America↗

Seismic-hazard maps for the conterminous United States, 2014

The maps presented here provide an update to the 2008 data contained in U.S Geological Survey Scientific Investigations Map 3195 ( http://pubs.usgs.gov/sim/3195/ ).Probabilistic seismic-hazard maps were prepared for the conterminous United States for 2014 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 historical 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 was 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 sources.

Scientific Investigations Map↗

Geodesy- and geology-based slip-rate models for the Western United States (excluding California) national seismic hazard maps

The 2014 National Seismic Hazard Maps for the conterminous United States incorporate additional uncertainty in fault slip-rate parameter that controls the earthquake-activity rates than was applied in previous versions of the hazard maps. This additional uncertainty is accounted for by new geodesy- and geology-based slip-rate models for the Western United States. Models that were considered include an updated geologic model based on expert opinion and four combined inversion models informed by both geologic and geodetic input. The two block models considered indicate significantly higher slip rates than the expert opinion and the two fault-based combined inversion models. For the hazard maps, we apply 20 percent weight with equal weighting for the two fault-based models. Off-fault geodetic-based models were not considered in this version of the maps. Resulting changes to the hazard maps are generally less than 0.05 g (acceleration of gravity). Future research will improve the maps and interpret differences between the new models.

Arizona;Colorado;Idaho;Montana;New Mexico;Nevada;O↗

Documentation for the 2014 update of the United States national seismic hazard maps

The national seismic hazard maps for the conterminous United States have been updated to account for new methods, models, and data that have been obtained since the 2008 maps were released (Petersen and others, 2008). The input models are improved from those implemented in 2008 by using new ground motion models that have incorporated about twice as many earthquake strong ground shaking data and by incorporating many additional scientific studies that indicate broader ranges of earthquake source and ground motion models. These time-independent maps are shown for 2-percent and 10-percent probability of exceedance in 50 years for peak horizontal ground acceleration as well as 5-hertz and 1-hertz spectral accelerations with 5-percent damping on a uniform firm rock site condition (760 meters per second shear wave velocity in the upper 30 m, V S30 ). In this report, the 2014 updated maps are compared with the 2008 version of the maps and indicate changes of plus or minus 20 percent over wide areas, with larger changes locally, caused by the modifications to the seismic source and ground motion inputs.

Open-File Report↗

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↗

Rupture model of the 2011 Mineral, Virginia, earthquake from teleseismic and regional waveforms

We independently invert teleseismic P waveforms and regional crustal phases to examine the finite fault slip model for the 2011 M w 5.8 Mineral, Virginia, earthquake. Theoretical and empirical Green's functions are used for the teleseismic and regional models, respectively. Both solutions show two distinct sources each about 2 km across and separated by 2.5 km. The source at the hypocenter is more localized in the regional model leading to a higher peak slip of 130 cm and higher average stress drop of 250 bars compared with 86 cm and 150 bars for the same source in the teleseismic model. Both sources are centered at approximately 8 km depth in the regional model, largely below the aftershock distribution. In the teleseismic model, the sources extend updip to approximately 6 km depth, into the depth range of the aftershocks. The rupture velocity is not well resolved but appears to be near 2.7 km/s.

Geophysical Research Letters↗

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↗

Uniform California earthquake rupture forecast, version 3 (UCERF3): the time-independent model

In this report we present the time-independent component of the Uniform California Earthquake Rupture Forecast, Version 3 (UCERF3), which provides authoritative estimates of the magnitude, location, and time-averaged frequency of potentially damaging earthquakes in California. The primary achievements have been to relax fault segmentation assumptions and to include multifault ruptures, both limitations of the previous model (UCERF2). The rates of all earthquakes are solved for simultaneously, and from a broader range of data, using a system-level "grand inversion" that is both conceptually simple and extensible. The inverse problem is large and underdetermined, so a range of models is sampled using an efficient simulated annealing algorithm. The approach is more derivative than prescriptive (for example, magnitude-frequency distributions are no longer assumed), so new analysis tools were developed for exploring solutions. Epistemic uncertainties were also accounted for using 1,440 alternative logic tree branches, necessitating access to supercomputers. The most influential uncertainties include alternative deformation models (fault slip rates), a new smoothed seismicity algorithm, alternative values for the total rate of M≥5 events, and different scaling relationships, virtually all of which are new. As a notable first, three deformation models are based on kinematically consistent inversions of geodetic and geologic data, also providing slip-rate constraints on faults previously excluded because of lack of geologic data. The grand inversion constitutes a system-level framework for testing hypotheses and balancing the influence of different experts. For example, we demonstrate serious challenges with the Gutenberg-Richter hypothesis for individual faults. UCERF3 is still an approximation of the system, however, and the range of models is limited (for example, constrained to stay close to UCERF2). Nevertheless, UCERF3 removes the apparent UCERF2 overprediction of M6.5–7 earthquake rates and also includes types of multifault ruptures seen in nature. Although UCERF3 fits the data better than UCERF2 overall, there may be areas that warrant further site-specific investigation. Supporting products may be of general interest, and we list key assumptions and avenues for future model improvements.

California↗

Model and parametric uncertainty in source-based kinematic models of earthquake ground motion

Four independent ground-motion simulation codes are used to model the strong ground motion for three earthquakes: 1994 M w 6.7 Northridge, 1989 M w 6.9 Loma Prieta, and 1999 M w 7.5 Izmit. These 12 sets of synthetics are used to make estimates of the variability in ground-motion predictions. In addition, ground-motion predictions over a grid of sites are used to estimate parametric uncertainty for changes in rupture velocity. We find that the combined model uncertainty and random variability of the simulations is in the same range as the variability of regional empirical ground-motion data sets. The majority of the standard deviations lie between 0.5 and 0.7 natural-log units for response spectra and 0.5 and 0.8 for Fourier spectra. The estimate of model epistemic uncertainty, based on the different model predictions, lies between 0.2 and 0.4, which is about one-half of the estimates for the standard deviation of the combined model uncertainty and random variability. Parametric uncertainty, based on variation of just the average rupture velocity, is shown to be consistent in amplitude with previous estimates, showing percentage changes in ground motion from 50% to 300% when rupture velocity changes from 2.5 to 2.9 km/s. In addition, there is some evidence that mean biases can be reduced by averaging ground-motion estimates from different methods.

Bulletin of the Seismological Society of America↗

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↗

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↗

Spatial seismicity rates and maximum magnitudes for background earthquakes

The background seismicity model is included to account for M 5.0 - 6.5 earthquakes on faults and for random M 5.0 ? 7.0 earthquakes that do not occur on faults included in the model (as in earlier models of Frankel et al., 1996, 2002 and Petersen et al., 1996). We include four different classes of earthquake sources in the California background seismicity model: (1) gridded (smoothed) seismicity, (2) regional background zones, (3) special fault zone models, and (4) shear zones (also referred to as C zones). The gridded (smoothed) seismicity model, the regional background zone model, and the special fault zones use a declustered earthquake catalog for calculation of earthquake rates. Earthquake rates in shear zones are estimated from the geodetically determined rate of deformation across an area of high strain rate. We use a truncated exponential (Gutenberg-Richter, 1944) magnitude-frequency distribution to account for earthquakes in the background models.

Open-File Report↗