We use interferometric synthetic aperture radar (InSAR) and broadband seismic waveform data to estimate a source model of the 11th July, 2004 M W 6.2 Zhongba earthquake, Tibet of China. This event occurred within th...We use interferometric synthetic aperture radar (InSAR) and broadband seismic waveform data to estimate a source model of the 11th July, 2004 M W 6.2 Zhongba earthquake, Tibet of China. This event occurred within the seismically active zone of southwestern Tibetan Plateau where the east-west extension of the upper crust is observed. Because of limitations in one pair of InSAR data available, there are trade-offs among centroid depth, rupture area and amount of slip. Available seismic data tightly constrain the focal mechanism and centroid depth of the earthquake but not the horizontal location. Together, two complementary data sets can be used to identify the actual fault plane, better constrain the slip model and event location. We first use regional seismic waveform to estimate point source mechanism, then InSAR data is used to obtain better location. Finally, a joint inversion of teleseismic P-waves and InSAR data is performed to obtain a distributed model. Our preferred point source mechanism indicates a seismic moment of ~2.2×10 18 N·m (~M W 6.2), a fault plane solution of 171° (342 ° )/42 ° (48 ° )/-83 ° (-97 ° ), corresponding to strike/dip/rake, and a depth of 11 km. The fault plane with strike of 171 ? and dip of 42° is identified as the ruptured fault with the aid of InSAR data. The preferred source model features compact area of slips between depth of 5–11 km and 10 km along strike with maximum slip amplitude of about 1.5 m.展开更多
Recent events beneath Central America have produced excellent sets of inner core reflection (PKiKP phase) at high frequency recorded by USArray ranging from 18° to 30°. However, the amplitude of this phase...Recent events beneath Central America have produced excellent sets of inner core reflection (PKiKP phase) at high frequency recorded by USArray ranging from 18° to 30°. However, the amplitude of this phase displays considerable scatter with a factor of six or more. Such scatter has been attributed to upper-mantle scattering and the Inner Core Boundary (ICB) in combination. Here, we show that neighboring events share upper-mantle scatterers beneath the receivers, and their ratio allows a clearer image of deep earth structure. Alter confirming some of the measured variation is indeed due to deep structure, we stacked nearby traces to reduce fine scale variations which are mostly due to shallow structure. Then, the remaining relatively large scale variation pattern of PKiKP phase is caused by the inner core boundary, as demonstrated by numerical experiments. After migration of data to the 1CB, we observe a consistent image. We find such a pattern can be explained by a patch of mushy material of a few kilometers high where the material changes gradually from that of the outer core to that of the inner core.展开更多
Earthquake magnitude and rupture duration are among the most important parameters characterizing an earthquake for the purpose of early tsunami warning. While they can be routinely determined from broadband P waveform...Earthquake magnitude and rupture duration are among the most important parameters characterizing an earthquake for the purpose of early tsunami warning. While they can be routinely determined from broadband P waveforms with iterative inver- sion procedures, the inversion procedures may fail when the rupture either lasts longer than the interval between P and later arrivals or requires too much time or human intervention. Little contaminated by later arrivals, high frequency P waves are useful for modeling earthquake source processes, though the envelope waveform is affected by strong scattering in lithosphere. With high frequency envelopes from aftershocks as Empirical Green's Function (EGF), the coda effects can be removed and more accurate relative source time function (RSTF) of the main shock can be obtained. Assuming that RSTFs cannot be negative, we use the projected Landweber deconvolution method (PLD) to obtain high frequency RSTFs because PLD method has the advantage of non-negativity, causality, and compactness (finite duration). We are able to determine rupture durations of four large earthquakes: the 2004 Sumatra-Andaman earthquake, the 2005 Nias event, the 2006 Java event, and the 2011 Tokuko earthquake. The rupture durations of the Sumatra-Andaman, Nias, and Tohuko events are found to be around 550, 110, and 120 s respectively, consistent with previous studies. The rupture duration of the Java event is about 130 s, supporting that the Java event is a tsunami earthquake. The magnitudes of these earthquakes are found to depend on both the amplitude and the duration of the deconvolved waveforms, and can be approximated by integrating these waveforms.展开更多
基金study was supported jointly by National Natural Science Foundation of China (Nos.40821160549 and 41074032)CAS Knowledge Innovation Program (No. KZCX2-YW-116-1)Joint Seismological Science Fundation of China(Nos.200808078 and 200708035)
文摘We use interferometric synthetic aperture radar (InSAR) and broadband seismic waveform data to estimate a source model of the 11th July, 2004 M W 6.2 Zhongba earthquake, Tibet of China. This event occurred within the seismically active zone of southwestern Tibetan Plateau where the east-west extension of the upper crust is observed. Because of limitations in one pair of InSAR data available, there are trade-offs among centroid depth, rupture area and amount of slip. Available seismic data tightly constrain the focal mechanism and centroid depth of the earthquake but not the horizontal location. Together, two complementary data sets can be used to identify the actual fault plane, better constrain the slip model and event location. We first use regional seismic waveform to estimate point source mechanism, then InSAR data is used to obtain better location. Finally, a joint inversion of teleseismic P-waves and InSAR data is performed to obtain a distributed model. Our preferred point source mechanism indicates a seismic moment of ~2.2×10 18 N·m (~M W 6.2), a fault plane solution of 171° (342 ° )/42 ° (48 ° )/-83 ° (-97 ° ), corresponding to strike/dip/rake, and a depth of 11 km. The fault plane with strike of 171 ? and dip of 42° is identified as the ruptured fault with the aid of InSAR data. The preferred source model features compact area of slips between depth of 5–11 km and 10 km along strike with maximum slip amplitude of about 1.5 m.
基金supported by NSF EAR-1053064 and CSEDI EAR-1161046 at CalTech with partial support of D. Sun at USC under EAR-0809023
文摘Recent events beneath Central America have produced excellent sets of inner core reflection (PKiKP phase) at high frequency recorded by USArray ranging from 18° to 30°. However, the amplitude of this phase displays considerable scatter with a factor of six or more. Such scatter has been attributed to upper-mantle scattering and the Inner Core Boundary (ICB) in combination. Here, we show that neighboring events share upper-mantle scatterers beneath the receivers, and their ratio allows a clearer image of deep earth structure. Alter confirming some of the measured variation is indeed due to deep structure, we stacked nearby traces to reduce fine scale variations which are mostly due to shallow structure. Then, the remaining relatively large scale variation pattern of PKiKP phase is caused by the inner core boundary, as demonstrated by numerical experiments. After migration of data to the 1CB, we observe a consistent image. We find such a pattern can be explained by a patch of mushy material of a few kilometers high where the material changes gradually from that of the outer core to that of the inner core.
基金supported by Chinese Academy of Sciences (Grant No. KZCX2-YE-142)National Natural Science Foundation of China(Grant Nos. 40974034, 41174086, 41021003)Chinese Academy of Sciences (Grant No. KZCX2-YW-116)
文摘Earthquake magnitude and rupture duration are among the most important parameters characterizing an earthquake for the purpose of early tsunami warning. While they can be routinely determined from broadband P waveforms with iterative inver- sion procedures, the inversion procedures may fail when the rupture either lasts longer than the interval between P and later arrivals or requires too much time or human intervention. Little contaminated by later arrivals, high frequency P waves are useful for modeling earthquake source processes, though the envelope waveform is affected by strong scattering in lithosphere. With high frequency envelopes from aftershocks as Empirical Green's Function (EGF), the coda effects can be removed and more accurate relative source time function (RSTF) of the main shock can be obtained. Assuming that RSTFs cannot be negative, we use the projected Landweber deconvolution method (PLD) to obtain high frequency RSTFs because PLD method has the advantage of non-negativity, causality, and compactness (finite duration). We are able to determine rupture durations of four large earthquakes: the 2004 Sumatra-Andaman earthquake, the 2005 Nias event, the 2006 Java event, and the 2011 Tokuko earthquake. The rupture durations of the Sumatra-Andaman, Nias, and Tohuko events are found to be around 550, 110, and 120 s respectively, consistent with previous studies. The rupture duration of the Java event is about 130 s, supporting that the Java event is a tsunami earthquake. The magnitudes of these earthquakes are found to depend on both the amplitude and the duration of the deconvolved waveforms, and can be approximated by integrating these waveforms.