Stephen Hartzell

ORCID: 0000-0003-0858-9043
Publications
Citations
Views
---
Saved
---
About
Contact & Profiles
Research Areas
  • Seismic Waves and Analysis
  • earthquake and tectonic studies
  • Seismic Performance and Analysis
  • High-pressure geophysics and materials
  • Seismology and Earthquake Studies
  • Structural Health Monitoring Techniques
  • Seismic Imaging and Inversion Techniques
  • Landslides and related hazards
  • Earthquake Detection and Analysis
  • Geophysics and Sensor Technology
  • Geological and Geochemical Analysis
  • Geotechnical Engineering and Underground Structures
  • Underwater Acoustics Research
  • Geotechnical Engineering and Soil Mechanics
  • Geotechnical and Geomechanical Engineering
  • Robotics and Sensor-Based Localization
  • Vibration Control and Rheological Fluids
  • Meteorological Phenomena and Simulations
  • Machine Fault Diagnosis Techniques
  • Structural Response to Dynamic Loads
  • Tree Root and Stability Studies
  • Polydiacetylene-based materials and applications
  • Wind and Air Flow Studies
  • Geophysical Methods and Applications
  • Sparse and Compressive Sensing Techniques

United States Geological Survey
2013-2024

Denver Federal Center
2011-2024

Geologic Hazards Science Center
2013

U.S. Air Force Institute of Technology
2012

California Institute of Technology
1979-1989

Scripps Institution of Oceanography
1977-1979

University of California, San Diego
1977-1979

Abstract A least-squares point-by-point inversion of strong ground motion and teleseismic body waves is used to infer the fault rupture history 1979 Imperial Valley, California, earthquake. The represented by a plane embedded in half-space where elastic properties vary with depth. yields both spatial temporal variations dislocation on for right-lateral strike-slip normal dip-slip components motion. Inversions are run different dips constant variable velocity models. Effects data sets also...

10.1785/bssa07306a1553 article EN Bulletin of the Seismological Society of America 1983-12-01

A method is presented for modeling earthquake strong ground motion, which uses the aftershocks associated with a large as Green's functions. major earthquake, rupture surface, modeled by collection of point sources distributed over fault plane. The response each source approximated motion closest aftershock. By using aftershock responses, effects true earth structure are included in process. This used to model El Centro displacement record 1940 Imperial Valley earthquake.

10.1029/gl005i001p00001 article EN Geophysical Research Letters 1978-01-01

We have compared aftershock patterns following several moderate to large earthquakes with the corresponding distributions of coseismic slip obtained from previous analyses recorded strong ground motion and teleseismic waveforms. Well-located hypocenters are projected onto main shock fault plane, their positions examined relative zones displacement indicated by estimated slip. also examine focal mechanisms, when these data available, in an attempt identify possible secondary faulting within...

10.1785/bssa0780041438 article EN Bulletin of the Seismological Society of America 1988-08-01

Abstract This article compares techniques for calculating broadband time histories of ground motion in the near field a finite fault by comparing synthetics with strong-motion data set 1994 Northridge earthquake. Based on this comparison, preferred methodology is presented. Ground-motion-simulation are divided into two general methods: kinematic- and composite-fault models. Green's functions three types evaluated: stochastic, empirical, theoretical. A hybrid scheme found to give best fit...

10.1785/bssa0890061484 article EN Bulletin of the Seismological Society of America 1999-12-01

Abstract Near-source strong motion velocity records and teleseismic short-period P waveforms are modeled to obtain the spatial temporal distribution of slip for 1984 Morgan Hill earthquake. Both forward modeling constrained, least-squares inversion techniques used interpret in frequency range approximately 0.2 2.0 Hz. These data support a nearly unilateral rupture southeast with propagation nine-tenths local S-wave velocity. The majority occurs over fault length 25 km first approximation can...

10.1785/bssa0760030649 article EN Bulletin of the Seismological Society of America 1986-06-01

We present a new method for calculating broadband time histories of ground motion based on hybrid low-frequency/high-frequency approach with correlated source parameters. Using finite-difference we calculate low- frequency synthetics (< ∼1 Hz) in 3D velocity structure. also compute 1D model using frequency-wavenumber method. The low frequencies from the calculation are combined high by matched filtering at crossover 1 Hz. description, common to both and synthetics, is random distributions...

10.1785/0120060036 article EN Bulletin of the Seismological Society of America 2006-12-01

This study examines the rationale of Hilbert-Huang transform (HHT) for analyzing dynamic and earthquake motion recordings in studies seismology engineering. In particular, this paper first provides fundamentals HHT method, which consist empirical mode decomposition (EMD) Hilbert spectral analysis. It then uses to analyze hypothetical real wave motion, results are compared with obtained by Fourier data processing technique. The analysis two indicates that method is able extract some...

10.1061/(asce)0733-9399(2003)129:8(861) article EN Journal of Engineering Mechanics 2003-07-22

Abstract We have inverted the teleseismic P waveforms recorded by stations of Global Digital Seismograph Network for 8 July 1986 North Palm Springs, California, 28 October 1983 Borah Peak, Idaho, and 19 September 1985 Michoacan, Mexico, earthquakes to recover distribution slip on each faults using a point-by-point inversion method with smoothing positivity constraints. In procedure, fault plane fixed strike dip is placed in region earthquake hypocenter divided into large number subfaults....

10.1785/bssa0780031092 article EN Bulletin of the Seismological Society of America 1988-06-01

Large subduction earthquakes on the Cascadia zone pose a potential seismic hazard. Very young oceanic lithosphere (10 million years old) is being subducted beneath North America at rate of approximately 4 centimeters per year. The shares many characteristics with zones in southern Chile, southwestern Japan, and Colombia, where comparably also subducting. large earthquakes, ranging energy magnitude (M(w)) between 8 9.5, have occurred along these other zones. If storing elastic energy,...

10.1126/science.236.4798.162 article EN Science 1987-04-10

A hybrid global search algorithm is used to solve the nonlinear problem of calculating slip amplitude, rake, risetime, and rupture time on a finite fault. Thirty‐five strong motion velocity records are inverted by this method over frequency band from 0.1 1.0 Hz for Northridge earthquake. Four regions larger‐amplitude identified: one near hypocenter at depth 17 km, second west about same depth, third updip 10 fourth northwest. The results further show an initial fast with 2.8 3.0 km/s...

10.1029/96jb01883 article EN Journal of Geophysical Research Atmospheres 1996-09-10

Abstract We examined short-period P waves to investigate if waveform data could be used determine which of two nodal planes was the actual fault plane for a small (ML 4.6) earthquake near Upland, California. removed path and site complications by choosing aftershock 2.7) as an empirical Green function. The main shock were deconvolved using function produce simple far-field displacement pulses. least-squares method invert these pulses slip distribution on finite fault. Both (strike 125°, dip...

10.1785/bssa0800030507 article EN Bulletin of the Seismological Society of America 1990-06-01

Abstract Following the 1989 Loma Prieta earthquake, a dense array of seven digitally recorded, three-component seismograph stations was deployed on Robinwood Ridge 7.3 km northwest epicenter. The purpose this to investigate cause high levels structural damage and ground cracking observed ridge crest. Aftershocks recorded by allow comparison motion up slope from base data present an extremely complicated pattern that demonstrates importance three-dimensionality problem. Slowness analysis P...

10.1785/bssa0840051336 article EN Bulletin of the Seismological Society of America 1994-10-01

The July 8, 1986, North Palm Springs earthquake is used as a basis for comparison of several different approaches to the solution rupture history finite fault. inversion waveform data considered; both teleseismic P waveforms and local strong ground motion records. Linear parametrizations slip amplitude are compared with nonlinear time. Inversions using synthetic empirical Green's functions considered. In general, accurate more readily calculable problem where simple ray theory flat‐layered...

10.1029/jb094ib06p07515 article EN Journal of Geophysical Research Atmospheres 1989-06-10

abstract Near-source ground motion at four azimuths but constant epicentral range is computed from a buried circular strike-slip fault in half-space. Particle acceleration, velocity, and displacement each station on the free surface frequency band 0.0 to 5.0 Hz. The assumed dislocation derived Kostrov (1964) function for continuously propagating stress relaxation. azimuthal variations amplitudes waveforms directly result spatially varying slip fault, radiation pattern over magnitude...

10.1785/bssa0710040939 article EN Bulletin of the Seismological Society of America 1981-08-01

Twelve three-component strong-motion displacement records are modeled for the 1979 Imperial Valley earthquake to recover distribution of slip on fault plane. The final model, which point source responses calculated by a discrete wavenumber/finite element technique, uses structure with gradients in material properties rather than layers. effects velocity gradient investigated comparing synthetics layer-over-a-half-space model using generalized rays. It is shown that uniform rupture...

10.1785/bssa0720020571 article EN Bulletin of the Seismological Society of America 1982-04-01

Abstract We simultaneously invert the strong-motion velocity records and long- intermediate-period teleseismic P waveforms of 19 September 1985 Michoacan, Mexico, earthquake to recover distribution slip on fault using a point-by-point constrained stabilized, least-squares inversion method. A plane with strike fixed at 300° dip 14° is placed in region hypocenter divided into 120 subfaults. Rupture assumed propagate 2.6 km/sec away from hypocenter. Synthetic near-source ground motions for pure...

10.1785/bssa0790030655 article EN Bulletin of the Seismological Society of America 1989-06-01

Abstract We estimate stress drops for earthquakes in and near the continental United States using method of spectral ratios. The ratio acceleration spectra between collocated recorded at a given station removes effects path recording site yields source parameters including corner frequency for, seismic moment between, two earthquakes. determine drop from these 1121 greater than M ∼3 60 earthquake clusters. find that average Brune few eastern (EUS) tectonic mainshocks studied (2.6–36 MPa) is...

10.1785/0120160219 article EN Bulletin of the Seismological Society of America 2017-02-21

Abstract Ground-motion records from aftershocks of the 1994 Northridge earthquake are used to estimate site response in urban Los Angeles area. Over 1300 shear-wave 61 sources and 90 sites a linear inversion for source site-response spectra. The methodology makes no assumptions about shape spectrum. To obtain stable unique inverse, Q model geometrical spreading factor assumed. In addition, at hard-rock is constrained be approximately 1.0 with kappa 0.02. spectra compare favorably results...

10.1785/bssa08601bs168 article EN Bulletin of the Seismological Society of America 1996-02-01

Abstract Aftershocks of the 1989 Loma Prieta, California, earthquake are used to estimate site response along San Francisco Peninsula. A total 215 shear-wave records from 24 sources and 21 sites in a linear inversion for source spectra. The methodology makes no assumptions about shape spectrum. However, obtain stable, unique inverse Q model geometrical spreading factor assumed, as well constraint on that sets averaged over two specific stations 1.0. Site responses calculated by this...

10.1785/bssa0820062308 article EN Bulletin of the Seismological Society of America 1992-12-01

Mathematical models of soil nonlinearity in common use and recently developed nonlinear codes are compared to investigate the range their predictions. We consider equivalent linear formulations with without frequency-dependent moduli damping ratios for total effective stress. Average velocity profiles 150 m depth midrange National Earthquake Hazards Reduction Program site classifications (B, BC, C, D, E) top 30 used compare response a wide conditions from rock soft soil. Nonlinear using...

10.1785/012003256 article EN Bulletin of the Seismological Society of America 2004-10-01

Abstract This article addresses the variability of site response in Los Angeles area and possible structural causes for observations. Aftershock records from 231 sites San Fernando basins surrounding mountains are used this study. Spectral ratios, taken with respect to a low-amplitude reference site, document variation amplification frequency range 2 6 Hz, both spatially backazimuth source. At higher frequencies (6 10 Hz), spectral ratios shown have greater spatial variability. Interstation...

10.1785/bssa0870061377 article EN Bulletin of the Seismological Society of America 1997-12-01

Abstract Strong motion records of the 1989 Loma Prieta earthquake are inverted to determine a model rupture history. Uncorrected horizontal and vertical accelerograms integrated particle velocity time histories for 38 stations within an epicentral range 75 km. The bandpassed filtered with corners at 0.05 1.0 Hz. These using nonlinear method solve distribution slip amplitudes times specified locations on fault plane. plane is priori: km long 17 wide, extending from 3 19 depth constant dip...

10.1785/bssa0810051573 article EN Bulletin of the Seismological Society of America 1991-10-01

Abstract Broadband, digital, teleseismic body waves ( P and SH) are inverted to obtain the history of slip for 17 October 1989, Loma Prieta earthquake. Comparison finite-fault, waveform, inversion results using an L 1 norm 2 reveal some basic uncertainties in rupture that not usually appreciated earthquake source studies. However, more robust features can be identified less dependent on choice minimization norm. From a consideration first-motion data forward modeling waves, strike dip our...

10.1785/bssa0810051518 article EN Bulletin of the Seismological Society of America 1991-10-01

Abstract A pair of significant earthquakes occurred on conjugate faults in the western Imperial Valley involving through-going Superstition Hills fault and Elmore Ranch cross fault. The first event was located fault, Ms = 6.2, larger 6.6. latter is seen as a doublet teleseismically with amplitudes ratio 1:2 delayed by about 8 sec. This 8-sec delay also dozen strong-motion records. These records are used constrained least-squares inversion scheme to determine distribution slip 2-D Upon closer...

10.1785/bssa0800051079 article EN Bulletin of the Seismological Society of America 1990-10-01

Broadband vertical P-wave records are analyzed from 63 of the largest shallow subduction zone earthquakes which have occurred in circum-Pacific last 45 yr. Most studied come a common instrument, Pasadena, California, Benioff 1-90 seismometer. Propagation and instrument effects deconvolved using damped least-squares inversion to obtain teleseismic source time function. The has additional constraint that function be positive everywhere. period band over functions considered reliable is 2.5 50...

10.1785/bssa0750040965 article EN Bulletin of the Seismological Society of America 1985-08-01
Coming Soon ...