html5-img
1 / 42

The Rupture Process of the August 23, 2011 Virginia Earthquake

The Rupture Process of the August 23, 2011 Virginia Earthquake . Martin Chapman Virginia Tech. Circles show mainshock and early (August 24-26) aftershock epicenters from a 7 station temporary deployment of instruments by VPI &SU. The beachball diagram indicates the

barb
Download Presentation

The Rupture Process of the August 23, 2011 Virginia Earthquake

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. The Rupture Process of the August 23, 2011 Virginia Earthquake Martin Chapman Virginia Tech

  2. Circles show mainshock and early (August 24-26) aftershock epicenters from a 7 station temporary deployment of instruments by VPI &SU. The beachball diagram indicates the mainshock focal mechanism (USGS/SLU). Colors show major geologic units.

  3. North Anna Power Station 21 km A Aftershocks between August 24 and September 2 2011

  4. Aftershock hypocenters, Aug. 24 - Sept. 2, 2011 Looking broadside at the fault plane (strike N29E, dip 51 deg) Looking approximately N29E

  5. Strong motion recordings on the foundation base mat of the unit 1 containment structure at North Anna. The records are remarkable for the very short duration, and pulse-like character of the largest motion--- the S wave(s). Note the modulated character of the Fourier amplitude spectra.

  6. Recordings from Unit 1, Dominion, Inc., North Anna power station 23 km from epicenter A weak initiation phase and two strong pulses are apparent on the transverse component

  7. Observed Transverse Full wave-field synthetic Chapman and Godbee, BSSA 2012.

  8. 17 stations at distances beyond 20 degrees showed clear sub-event arrivals

  9. North Anna Transverse Record And Source function used to Model the teleseismic data Broadband velocity recordings at Selected Teleseismic Stations Observed data (transverse) 1.22 x 1024 dyne-cm 300 bars 0.25 x 1024 dyne-cm 50 bars 2.25 x 1024 dyne-cm 300 bars

  10. We wish to determine the actual locations and origin times of the 3 principle subevents comprising the major moment release. The process used is as follows: 1) determine arrival times of the two later subevents relative to the initial (small) subevent at local and teleseismic receiver locations. 2) Assume that all subevents occur on the same fault plane, strike N29oE, dip 51o, and solve for the locations of the two later events relative to the initial subevent. 3) Locate the epicenter of the mainshock initial subevent relative to the well-determined epicenters of the best recorded aftershocks, using the handful of local stations that recorded the mainshock as well as the aftershocks. 4) Determine the focal depth of the mainshock by comparing actual teleseismic waveforms with synthetics. (all of this requires an accurate velocity model of the shallow crust)

  11. Determination of P and S wave Velocities in the Aftershock Zone Steps 1) Assume a half-space structure with an initial estimate of P and S wave velocities. 2) Get initial estimates of 36 aftershock hypocenter locations using the initial velocity model and P and S arrival time data. 3) Determine joint estimates of aftershock origin times and P/S velocity ratio, following the Wadati approach. 4) Use the initial epicenter locations and the joint estimates of the origin times from step (3) to obtain joint estimates of the P wave velocity and the focal depths of the aftershocks. 5) Update the velocity model using the estimate of P-wave velocity from 4), and velocity ratio from step 3), fix the aftershock origin times using the results of step 3) and repeat the process from step 2). After a few iterations, the process converged to an estimate of P-wave velocity (the values of the parameter estimates of subsequent iterations stopped changing).

  12. Joint Wadati-type inversion for origin times and P/S velocity ratio solve for slope and intercept terms for i = 1 to 36 aftershocks using least squares. Joint T2 versus X2 inversion for P-wave velocity and focal depth solve for the slope and intercept terms for i =1 to 36 aftershocks

  13. Joint estimate of Vp/Vs= 1.69 Plot of S-P arrival time intervals versus P arrival time for 36 aftershocks P-wave velocity = 5.96 km/s Plot of residuals from the Joint estimate of P-wave velocity and focal depth.

  14. s = strike d = dip a = epicenter - station azimuth a = P-wave velocity, P = ray parameter Vop = position vector of subevent w/r to origin at location of initial subevent. Vr = unit vector in direction of ray. tp - to = difference in origin times of later subevent and initial subevent. tp-o = observed arrival time of later subevent - observed arrival time of initial subevent. The unknowns are tp - to , I, and J, defining the origin time difference and the location of the later subevent relative to the initial subevent in the fault plane. With observations of tp-o at 4 or more stations, the unknowns can be determined by the method of least squares.

  15. Results of 500 solutions using 1/2 the total data set, randomly sampled with substitution (bootstrap) Locations constrained to lie on a fault plane striking N29oE, dip 51o location of initial subevent, at origin up-dip direction strike direction location of 2nd subevent X = 0.98 +/0.12 km, Y = 0.85 +/- 0.17 km origin time difference: 0.75 +/- 0.01 sec. location of 3rd subevent X = 1.48 +/- 0.27 km, Y = 1.33 +/- 0.39 km origin time difference: 1.57 +/- 0.03 sec.

  16. Permanent Local Stations that Recorded the Mainshock and Larger Aftershocks

  17. Calculation of Synthetics • Assume Vp = 5.96 km/sec Vs = 3.52 km/sec, a focal depth and focal mechanism • strike, dip and rake. (N29E, 51 deg to the southeast, 113 deg. • Determine Ray Parameter P for surface focus, at given epicentral distance • For given focal depth, determine the time intervals between the direct P and the • surface reflections TsP - TP and TpP - TP, assuming constant ray parameter. • Generate a source pulse for a given subevent in the frequency domain , S(w). • Use a causal attenuation operator in the frequency domain... A(w) Signal (w) = S(w) A(w) [Rp + RsPSP exp(iw(Tsp-TP) + RpPPPexp(iw(TpP-TP)] (Geometrical Spreading) Where PP and SP are the free surface reflection coefficients, and Rp, RsP and RpP are the radiation pattern values for direct P, sP and pP. The process is repeated for each subevent, with time shifts accounting for the relative locations and origin times of the subevents, and the results are summed. Transform to time domain and compare with data

  18. Velocity and Acceleration waveforms at Eskdalemuir, Scotland, compared with synthetics for 3 different depths.

  19. Eskdalemuir, Scotland, NE az Alert, Canadian Arctic, due North az Yellowknife, Northwest Territory, NW az

  20. Early aftershock locations using data from 7 VPI&SU stations: (August 26 - September 2, 2011) Looking N29oE, in strike direction subevent 3 subevent 2 mainshock initiation 1 Looking N61oW, normal to strike 3 2 1

  21. Conclusions • Most of the moment release occurred in 3 distinct slip events. • The mainshock occurred on the southwestern end of the early aftershock zone. • The aftershocks define a tabular zone with orientation in good agreement with the • mainshock focal mechanism. • Most of the mainshock moment release is from two sub-events in the depth range • 7.0 to 7.5 km, deeper than most of the early aftershocks which lie in the range from • near-surface to 7 km, averaging about 4 km. These aftershocks appear to occur in a • part of the planar zone that lies in the direction of mainshock rupture propagation. • The rupture sequence initiated at approximately 8 km, and proceeded up-dip and • and to the northeast along strike with a velocity of less than 2.0 km/sec, based • on the estimates of origin time and location of the subevents. • The implied fault rupture area is small (approximately 2.0 km between subevent 1 • and subevent 3), but the total rise time (based on the North Anna recording) is • about 3.0 seconds.

  22. Preliminary Results Concerning High-Frequency Ground Motion From the Aftershock Sequence

  23. Vertical component (acceleration) recordings of the 16:54, Sept. 5, 2011 UTC aftershock recorded by the AIDA profile stations

  24. 42.6 km 3.8 km (Left) Acceleration recorded at AIDA profile station 3120 at 3.8 km from the epicenter of the 16:54, Sept. 5, 2011 UTC aftershock. (Right) Acceleration recorded at station 3510, 42.6 km from the epicenter.

  25. (Left) Fourier amplitude spectra of the S wave and pre-P wave noise at AIDA profile station 3120. (Right) station 3510.

  26. AIDA Profile: Blue triangles YC network (IRIS Ramp Stations): Squares XY network (VT stations): Triangles

  27. The profile data are important because they provide a substantial amount of control in regard to source radiation pattern effects. The aftershocks are less than magnitude 3.0, and the spectra are very rich in high-frequency information. • Geometrical spreading near the source • Kappa • Q

  28. Band-Pass Filtered Max. Amplitudes of the 16:54, Sept. 5, 2011 UTC aftershock PGA, geometric mean of the horizontal comp. 8 -16 Hz pass-band, Geometric mean horizonal AIDA Profile: Black YC network: Green XY network: Red

  29. Band-Pass Filtered Max. Amplitudes recordings of the 16:54, Sept. 5, 2011 UTC aftershock 4 - 8 Hz pass-band, Geometric mean horizontal 2 - 4 Hz pass-band, Geometric mean horizontal bogus bogus AIDA Profile: Black YC network: Green XY network: Red

  30. Fourier Amplitude of the E-W component S wave along 50 km profile from near the 2011 epicenter to Corbin, VA 2 Hz Fourier Amp. 4 Hz Fourier Amp. r-1.1 r-1 r-1 r-1.1 Red symbols are Q corrected Using Q=893f0.32 (A&B 2006)

  31. 8 Hz Fourier Amp. 16 Hz Fourier Amp. r-1.0 r-1.0 r-1.6 r-1.8 Red symbols are Q corrected Using Q=893f0.32 (A&B 2006)

  32. Theoretical Geometrical Spreading, from Chapman and Godbee (2012) BSSA, in press

  33. Theoretical Geometrical Spreading, from Chapman and Godbee (2012) BSSA, in press The modeling results indicate that geometrical spreading around 1 Hz is approximately r-1.0, whereas it is approximately r-1.6 at higher frequencies for “rock” velocity models. This appears to be consistent with the VA aftershock data.

  34. 20110903_2110 E-W comp. Pre- P noise window S wave window as recorded S-wave window corrected for r-1.6 geometrical spreading

  35. 20110903_2110 N-S comp. Pre- P noise window S wave window as recorded S-wave window corrected for r-1.6 geometrical spreading

  36. 20110903_2110 Vertical comp. Pre- P noise window S wave window as recorded S-wave window corrected for r-1.6 geometrical spreading

  37. Preliminary analysis suggests: • Geometrical spreading from shallow (e.g., 5-6 km) focal depth • aftershocks appears to be frequency dependent: • r -1.0 for < 4 Hz, • r -1.6 for > 4 Hz • 2) The high-frequency acceleration spectra appear flat to 40Hz. • Apparent Kappa is very small, Q is apparently very high* • *These results are preliminary and conditional on spectral ratio • analysis for confirmation.

  38. Thank You!

More Related