Source Parameters of Weak Crustal Earthquakes of the Vrancea Region from Short-period Waveform Inversion

High-frequency records of nine low magnitude shallow earthquakes of the Vrancea (Romania) seismic region are inverted for the seismic moment tensor (MT). An approach is suggested regarding how to obtain at least a rough estimate of the MT when the information on the structure of the crust is poor. Here simple 1-D layered models are used in the Greens function synthesis despite the fact that the structure of the region is undoubtedly very complex. Different 1-D models were used for different source-station paths to approximate lateral variations. Record of a station located on a ray path which crosses a structure differing substantially from a 1-D model may however bias the retrieved MT essentially. Therefore, we did not collectively process all records, but subsets of stations separately. We check the consistency of the MTs resulting from these individual bootstrap solutions, and reject those which differ substantially, assuming that this is due to the oversimplification of the forward modeling. Thereafter, the averaged moment tensor yielded by the consistent subsets is accepted. Moreover, the distribution of the T, N and P axes from the moment tensors used for averaging provides a rough estimate of their reliability.Following this simplistic procedure, we found seven acceptably constrained solutions among nine events processed. Their P axes are compared with the general trend of the stress in the area: three comply with it, the others should be attributed to the complex stress field active in the region.Acknowledgment This research has been made possible by MURST (40% and 60%), by UNESCO-IGCP project 414 Realistic Modelling of Seismic Input for Megacities and Large Urban Areas and NATO SfP 972266. Partial support issued from the Grant Agency of Czech Rep. (Grant 205/02/0383) and from the National Agency for Science, Technology and Innovation of Romania (Grant 6185 GR/2000).     

Calibration of the Tibetan Plateau Using Regional Seismic Waveforms

We use the recordings from 51 earthquakes produced by a PASSCAL deployment in Tibet to develop a two-layer crustal model for the region. Starting with their ISC locations, we iteratively fit the P-arrival times to relocate the earthquakes and estimate mantle and crustal seismic parameters. An average crustal P velocity of 6.2–6.3 km/s is obtained for a crustal thickness of 65 km while the P velocity of the uppermost mantle is 8.1 km/s. The upper layer of the model is further fine-tuned by obtaining the best synthetic SH waveform match to an observed waveform for a well-located event. Green's functions from this model are then used to estimate the source parameters for those events using a grid search procedure. Average event relocation relative to the ISC locations, excluding two poorly located earthquakes, is 16 km. All but one earthquake are determined by the waveform inversion to be at depths between 5 and 15 km. This is 15 km shallower, on average, than depths reported by the ISC. The shallow seismicity cut-off depth and low crustal velocities suggest high temperatures in the lower crust. Thrust faulting source mechanisms dominate at the margins of the plateau. Within the plateau, at locations with surface elevations less than 5 km, source mechanisms are a mixture of strike-slip and thrust. Most events occurring in the high plateau where elevations are above 5 km show normal faulting. This indicates that a large portion of the plateau is under EW extension.     

Crustal Structure of the Korean Peninsula Using Surface Wave Dispersion and Numerical Modeling

Surface wave dispersion is studied to obtain the 1-D average velocity structure of the crust in the Korean Peninsula by inverting group- and phase-velocities jointly. Group velocities of short-period Rayleigh and Love waves are obtained from cross-correlations of seismic noise. Multiple-filter analysis is used to extract the group velocities at periods between 0.5 and 20 s. Phase velocities of Rayleigh waves in 10- and 50-s periods are obtained by applying the two-station method to teleseismic data. Dispersion curves of all group and phase velocities are jointly inverted for the 1-D average model of the Korean Peninsula. The resultant model from surface wave analysis can be used as an initial model for numerical modeling of observations of North Korean events for a velocity model appropriated to the Korean Peninsula. The iterative process is focused especially on the surface sedimentary layer in the numerical modeling. The final model, modified by numerical modeling from the initial model, indicates that the crust shear wave velocity increases with depth from 2.16 km/s for a 2-km-thick surface sedimentary layer to 3.79 km/s at a Moho depth of 33 km, and the upper mantle has a velocity of 4.70 km/s.     

How to Invert Multi-Band,Regional Phase Amplitudes for 2-D Attenuation and Source Parameters: Tests Using the USArray

We inverted for laterally varying attenuation, absolute site terms, moments and apparent stress using over 460,000 Lg amplitudes recorded by the USArray for frequencies between 0.5 and 16 Hz. Corner frequencies of Wells, Nevada, aftershocks, obtained by independent analysis of coda spectral ratios, controlled the tradeoff between attenuation and stress, while independently determined moments from St. Louis University and the University of California constrained absolute levels. The quality factor, Q, was low for coastal regions and interior volcanic and tectonic areas, and high for stable regions such as the Great Plains, and Colorado and Columbia Plateaus. Q increased with frequency, and the rate of increase correlated inversely with 1-Hz Q, with highest rates in low-Q tectonic regions, and lowest rates in high-Q stable areas. Moments matched independently determined moments with a scatter of 0.2 NM. Apparent stress ranged from below 0.01 to above 1 MPa, with means of 0.1 MPa for smaller events, and 0.3 MPa for larger events. Stress was observed to be spatially coherent in some areas; for example, stress was lower along the San Andreas fault through central and northern California, and higher in the Walker Lane, and for isolated sequences such as Wells. Variance reduction relative to 1-D models ranged from 50 to 90 % depending on band and inversion method. Parameterizing frequency dependent Q as a power law produced little misfit relative to a collection of independent, multi-band Q models, and performed better than the omega-square source parameterization in that sense. Amplitude residuals showed modest, but regionally coherent patterns that varied from event to event, even between those with similar source mechanisms, indicating a combination of focal mechanism, and near source propagation effects played a role. An exception was the Wells mainshock, which produced dramatic amplitude patterns due to its directivity, and was thus excluded from the inversions. The 2-D Q plus absolute site models can be used for high accuracy, broad area source spectra, magnitude and yield estimation, and, in combination with models for all regional phases, can be used to improve discrimination, in particular for intermediate bands that allow coverage to be extended beyond that available for high frequency P-to-S discriminants.     

Tsunami Source of the 2010 Mentawai, Indonesia Earthquake Inferred from Tsunami Field Survey and Waveform Modeling

The 2010 Mentawai earthquake (magnitude 7.7) generated a destructive tsunami that caused more than 500 casualties in the Mentawai Islands, west of Sumatra, Indonesia. Seismological analyses indicate that this earthquake was an unusual “tsunami earthquake,” which produces much larger tsunamis than expected from the seismic magnitude. We carried out a field survey to measure tsunami heights and inundation distances, an inversion of tsunami waveforms to estimate the slip distribution on the fault, and inundation modeling to compare the measured and simulated tsunami heights. The measured tsunami heights at eight locations on the west coasts of North and South Pagai Island ranged from 2.5 to 9.3 m, but were mostly in the 4–7 m range. At three villages, the tsunami inundation extended more than 300 m. Interviews of local residents indicated that the earthquake ground shaking was less intense than during previous large earthquakes and did not cause any damage. Inversion of tsunami waveforms recorded at nine coastal tide gauges, a nearby GPS buoy, and a DART station indicated a large slip (maximum 6.1 m) on a shallower part of the fault near the trench axis, a distribution similar to other tsunami earthquakes. The total seismic moment estimated from tsunami waveform inversion was 1.0 × 10²¹ Nm, which corresponded to M_w 7.9. Computed coastal tsunami heights from this tsunami source model using linear equations are similar to the measured tsunami heights. The inundation heights computed by using detailed bathymetry and topography data and nonlinear equations including inundation were smaller than the measured ones. This may have been partly due to the limited resolution and accuracy of publically available bathymetry and topography data. One-dimensional run-up computations using our surveyed topography profiles showed that the computed heights were roughly similar to the measured ones.     

多震相走时和波形三维地震成象方法及在祁连山地壳结构研究中的应用

本文的理论方法是以几何射线理论为基础发展起来的、天然地震走时反演技术及天然地震层析成像技术。它存在着震源函数与介质参数的解耦问题。本研究采取了五种方法来改善反演结果。包括,利用,Pg,Pn等震相增大约束条件;用已有精度较高的人工地震测深结果作速度约柬：用波形反演来修改模型,把诸多物理量开发出来互为约束,以修改后的模型再作反演,使解的稳定性大大提高：采用最优化过程,选择遗传算法。可以进行震源定位,走时反演,波形反演：得到任意深度的速度分布及从地表到Moho面的速度剖面。用于在祁连山地区的结果表明,这些层析剖面对认识大地构造、重大深部事件动力学是很有益的。     

Source Scaling of Inland Crustal Earthquake Sequences in Japan Using the S-Wave Coda Spectral Ratio Method

We estimate corner frequencies and stress drops for 298 events ranging from M _w 3.2–7.0 in 17 inland crustal earthquake sequences in Japan to investigate the source scaling and variation in stress drops. We obtain the source spectral ratio from observed records by the S-wave coda spectral ratio method. The advantage of using the S-wave coda is in obtaining much more stable source spectral ratios than using direct S-waves. We carefully examine the common shape of the decay of coda envelopes between event pair records. The corner frequency and stress drop are estimated by modeling the observed source spectral ratio with the omega-square source spectral model. We investigate the dependences of stress drops on some tectonic effects such as regionality, focal mechanism, and source depth. The principal findings are as follows: (1) a break in self-similar source scaling is found in our dataset. Events larger than M _w 4.5 show larger stress drops than those of smaller events. (2) Stress drops of aftershocks are mostly smaller than those of mainshocks in each sequence. (3) There are no systematic differences between stress drops of events occurring inside and outside the Niigata-Kobe Tectonic Zone in Japan. (4) Clear dependence of the faulting type on stress drops cannot be seen. (5) Stress drops of aftershocks depend on their source depth. (6) The crack size obtained from the corner frequency corresponds to the total rupture area of heterogeneous slip models for large events.     

The vP Velocity Structures of Sub-blocks in the North China Craton from Active Source Waveform Modeling and Tectonic Implications

We present the 1-D crustal velocity structure of the major tectonic blocks of the North China Craton(NCC) along 36°N based on synthetic seismogram modeling of long-range wide-angle reflection/refraction data. This profile extends from southwest Yan'an of central Shaanxi Province of China(109.47°E), across the southern Trans-North China Orogen(TNCO), the southwestern part of the North China Plain(NCP), the Luxi Uplift(LU) and the Sulu Orogen(SLO), ending at Qingdao City of Shandong Province, the eastern margin of China(120.12°E) along 36°N. We utilized reflectivity synthetic seismogram modeling of the active source data to develop 1-D velocity structures of the sub-blocks of the NCC. Our final model shows that the NCC crust varies remarkably among the tectonic units with different velocity structure features. Higher lower crustal velocity and Moho depth ～42 km is a major feature of the crust beneath southern Ordos Blockt. The TNCO which is composed of Lyuliangshan Mountains(LM), Shanxi Graben(SXG) and Taihangshan Mountains(TM) shows dominant trans-orogenic features. The NCP shows a dominant thickening of sediments, sharp crust thinning with Moho depth ～32 km and significant lower average velocity. The SLO and the LU shows a stratified crust, higher average velocity and crust thinning with Moho depth of ～35 km. Our model shows the coincidence between the deep structure and the surface geology among all the tectonic sub-blocks of the NCC.     

多震相走时和波形三维地震成像方法及在朝鲜半岛和祁连山地壳结构研究的应用

本项成果包括：提出天然地震走时反演层析成像技术,采用下列方法使得处理结果得以改善：1)利用Pg,Sg,Pm,Sm,Pn,Sn等震相增大约束条件;2)用已有精度较高的人工地震测深结果作速度约束;3)用波形反演来修改模型,把波源,介质吸收,散射等全部物理特征集中反映在记录中,把诸多物理量开发出来互为约束,以修改后的模型再作反演,使解的稳定性大大提高;4)采用最优化过程,选择遗传算法。可以进行震源定位,走时反演,波形反演;5)得到任意深度的速度分布及从地表到Moho面的速度剖面。主要应用结果：对于朝鲜半岛南部,划分为5(沿纬度)*６(沿经度)*8(沿深度)块,对于中部分288块.得到从地表到M面的8个水平切面;中部地区沿纬度13个二维剖面及其Moho面深度分布。上述方法也用于祁连山中东段地壳三维结构成像加上地震台网数字记录,反演。该区属塔里木-阿拉善地块走廊过渡带与北祁连褶皱带;从剖面可看出该地块上地壳低速层厚,下地壳有-低速层。北祁连褶皱带盆地与隆起构造之间的起伏差异,显示古浪断裂与金强河断裂之间的深部差异与界线。两地区的结果表明,这些剖面对认识大地构造、地质结构的稳定性,深部事件的性质是很有益的,对地球动力学研究也有重要意义。     

用S波研究雅江地震序列震源波谱

利用成都区域数字化地震台网记录的雅江地震序列的S波段资料，使用波谱分析方法计算雅江地震序列的震源参数。在对不同时段的震源参数进行初步分析后，得到如下结果：①给出了适合于雅江地震序列的地震矩M0和矩震级Mw与里氏震级ML的关系；②拐角频率从5．0级地震前的1．4左右下降到5．0级地震后的1．3左右，6．0级地震之后又下降到1．2左右；③前震和余震的一些震源参数存在明显差异，可作为预测未来强震的重要指标之一。     

甘肃岷县漳县6.6级地震前后区域地壳形变分析

利用陆态网GNSS连续站及流动站观测数据处理结果,从坐标时间序列、基线时间序列、速度场、形变场等方面分析岷县漳县地震前后青藏高原东北缘地区GPS形变场变化,得出以下结论:岷县漳县地震前GNSS测站时间序列没有明显变化,震前异常不突出;地震引起不同测点不同程度的永久性同震位移,最大值达到16mm;地震前后基线变化不明显,未发现明显的前兆异常;阿拉善及附近的祁连构造带存在顺时针运动的迹象,即运动方向由东侧的南西向逐渐变到西部的北西向,运动大小的范围基本保持在7mm/a之内;区域主压应变优势方向为东北-西南向,最大主压应变区呈带状分布于研究区域南部。     

Waveform Retrieval and Phase Identification for Seismic Data from the CASS Experiment

The little destruction to the deployment site and high repeatability of the Controlled Accurate Seismic Source (CASS) shows its potential for investigating seismic wave velocities in the Earth's crust. However, the difficulty in retrieving impulsive seismic waveforms from the CASS data and identifying the seismic phases substantially prevents its wide applications. For example, identification of the seismic phases and accurate measurement of travel times are essential for resolving the spatial distribution of seismic velocities in the crust. Until now, it still remains a challenging task to estimate the accurate travel times of different seismic phases from the CASS data which features extended wave trains, unlike processing of the waveforms from impulsive events such as earthquakes or explosive sources. In this study, we introduce a time-frequency analysis method to process the CASS data, and try to retrieve the seismic waveforms and identify the major seismic phases traveling through the crust. We adopt the Wigner-Ville Distribution (WVD) approach which has been used in signal detection and parameter estimation for linear frequency modulation (LFM) signals, and proves to feature the best time-frequency convergence capability. The Wigner-Hough transform (WHT) is applied to retrieve the impulsive waveforms from multi-component LFM signals, which comprise seismic phases with different arrival times. We processed the seismic data of the 40-ton CASS in the field experiment around the Xinfengjiang reservoir with the WVD and WHT methods. The results demonstrate that these methods are effective in waveform retrieval and phase identification, especially for high frequency seismic phases such as PmP and SmS with strong amplitudes in large epicenter distance of 80–120 km. Further studies are still needed to improve the accuracy on travel time estimation, so as to further promote applicability of the CASS for and imaging the seismic velocity structure.     

Characterization and Source History Modeling Using Low‐k Zone Profiles at Two Source Areas

Like tree rings, high‐resolution soil sampling of low‐permeability (low‐k) zones can be used to evaluate the style of source history at contaminated sites (i.e., historical pattern of concentration and composition vs. time since releases occurred at the interface with the low‐k zone). This is valuable for the development of conceptual site model (CSM) and can serve as an important line of evidence supporting monitored natural attenuation (MNA) as a long‐term remedy. Source histories were successfully reconstructed at two sites at Naval Air Station Jacksonville using a simple one‐dimensional (1D) model. The plume arrival time and historical composition were reconstructed from the time initial releases that were suspected to occur decades earlier. At the first site (Building 106), the source reconstructions showed relatively constant source concentrations, but significant attenuation over time in the downgradient plume in the transmissive zone, suggesting MNA may not be an appropriate remedy if source control is a requirement, but attenuation processes are clearly helping to maintain plume stability and reduce risk. At the second site (Building 780), source concentrations in the transmissive zone showed an approximately a one order of magnitude over time, but apparently less attenuation in the downgradient plume. The source reconstruction method appeared to reflect site remediation efforts (excavation, soil vapor extraction) implemented in the 1990s. Finally, a detailed analysis using molecular biological tools, carbon isotopes, and by‐products suggests that most degradation activity is associated with high‐k zones but not with low‐k zones at these source areas. Overall, the source reconstruction methodology provided insight into historical concentration trends not obtainable otherwise given the limited long‐term monitoring data.     

设定场源参数的区域重力微变化特征与模型实验

流动重力测量作为地震前兆监测手段之一, 所得到的区域重力场微变化信号反映地壳内部密度变化、 地壳变形、 流体运移等信息。 采用定量化和模型化方法研究与地震孕育相关的重力场变化信号时空特征, 对于构建地震预报定量化指标体系具有重要意义。 通过设定场源参数构建了数值模型, 正演得到了多组不同场源条件下的重力场变化模型; 通过模型对比, 分析了重力变化与重力梯度变化对不同类型场源参数的关系。 研究结果表明: 重力垂直梯度变化对异常体边界比较敏感, 特对浅部异常更为敏感, 适用于圈定异常体范围。 以此为基础提出了观测系统的优化测量方案, 以及对观测精度的需求。 研究结果对于构建地震预报指标参数, 优化流动重力观测系统具有一定参考意义。     

山西运城Ms4.3、3.4两次地震波形特征和震源新参数对比分析

对山西运城2016年3月12日Ms4.3和27日Ms3.4两次地震的波形特征和震源参数进行对比分析。结果显示,波形一致性很好,PmP波可能由于震源深度的差异,导致在同一台站的特征有所不同;震源机制解非常接近,说明两次地震的应力场作用方式一致。     

Inversion of Source Parameters for Moderate Earthquakes Using Short-Period Teleseismic P Waves

In this paper, we introduce a new method for estimating the source parameters of moderate earthquakes (M _w ~5.0) by modeling short-period teleseismic waveforms. This method uses a grid-search algorithm to minimize misfits between observed data and synthetic seismograms in depth, magnitude, and mechanism domain in a relative high-frequency range of 0.8–2.0 Hz, similar to the traditional cut-and-paste method used in regional modeling (Zhu and Helmberger, Bull Sesimol Soc Am 86:1634–1641, 1996). In this frequency range, a significant challenge is determining the initial P-wave polarity because of a low signal-to-noise ratio (SNR). Therefore we first determine source properties for a master earthquake with a relative strong SNR. Both the travel time and amplitude corrections are developed relative to the reference 1D model along each path used in inverting the master event. We then applied these corrections to other earthquakes clustered in the same area to constrain the initial P polarities. Thus the focal mechanisms can be determined reasonably well. We inverted focal mechanisms for a small set of events beneath Qeshm Island in southern Iran and demonstrate the importance of radiation pattern at short periods.     

推断2000年三宅岛火山喷发前后的地形变与形变源

2000年6月26日,三宅岛火山开始活动,随即通过GPS连续观测网对其进行了准实时的地形变监测。这是日本国内继北海道有珠山火山喷发后,在火山喷发前进行GPS连续观测的第二例,取得了不少有关岩浆流动和喷发机制的新认识。例如,通过对伴随     

云南省区域数字地震遥测台网远震参数测定

利用云南省区域数字地震遥测台网记录到的资料对远震的定位方法进行了研究。采用“分离变量法”,将震源参数分为发震时间、震中位置和震源深度三类变量,相对独立进行计算,这就减少了Geiger法中,各类变量因不同权重而产生的互相干扰以及因初始解不好引起的迭代不收敛的问题,改善了定位精度。由于大地震的地震波初动清楚,可精确测定初动到时,仅用经过台站高程走时改正和台站综合走时改正后的初动到时即可精确测定零深度时的震中位置。再根据测定的零深度震中位置,查出台网最小到时子台的S波到时S0它与该台的实际S波到时Sh的差可有效判定地震是否深源地震和计算震源深度。为了保持地震震级定义的一致性,采用数字地震仪系统的传递函数(包括幅频特性和相频特性)对实测数字地震记录进行仿真,得到真正的地动速度,再由计算机自动检测定地动速度最大面波振幅并计算震级Ms。