大同地震台远震波形走时残差的计算和应用

选取大同地震台记录的2010—2019年531个M_S≥5.0地震事件,根据发震区域的不同对其波形走时残差进行分析,得到不同发震区域地震波形的平均走时残差。结果表明,日本、喜马拉雅地区地震波形的走时残差总体为负值,平均残差分别为-1.51 s、-0.70 s;斐济、中国台湾地区、南美洲西海岸地区大多为正值,平均残差为0.79 s、1.83 s、1.86 s。同时发现,在有干扰的情况下参考不同区域平均走时残差计算结果,通过判断实际初至震相与理论初至震相的位置关系,可以更好地开展初至波震相识别工作。     

Travel time residuals in southeastern Europe

The Pn travel time relative residuals, in respect to a crustal model of the Aegean area, have been determined for 103 permanent seismological stations in southeastern Europe, western Turkey and the Middle East. The values of these residuals are considered to depend mainly on the crustal thickness beneath the seismological stations. Based on these values seven regions with different crustal thickness, varying between 31 Km and 42 Km, have been defined. The crust in these regions is continental. A region with very high negative residuals has been defined in the Middle East (Egypt, Israel, Lebanon). These negative residuals are attributed to different crustal structure of the eastern Mediterranean (oceanic crust with an extra thick sedimentary layer) and not the crustal thickness at the station sites.Independently from the interpretation, these Pn residuals can be used successfully to considerably improve (up to 2 Km) the determination of the earthquake foci locations.     

A regional study of mantle velocity variations beneath eastern Australia and the southwestern Pacific using short-period recordings of P,S, PcP,ScP and ScS waves produced by Tongan deep earthquakes

Deep earthquakes located in the Tonga-Kermadec region produce exceptionally clear and sharp short-period P, S, PcP, ScP, and ScS phases which are recorded at many stations at distances of less than 60°. The data used in this study are produced by short-period stations located in oceanic-type regions (Fiji and New Caledonia), a mobile continental region (eastern Australia) and a shield region (central Australia). Differential travel-time residuals of the above phases at these stations are investigated to determine the contribution to the differential residuals from: (1) the upper part of the mantle (S-P residuals); (2) the core-to-station portion of the mantle (ScS-ScP residuals); and (3) the hypocenter-to core portion of the mantle (ScP-PcP residuals). The use of differential travel-time residuals considerably reduces near-station effects and effects due to inaccurate determination of the source parameters, and hence the results can be interpreted as due to variations along the propagation paths. The results show that (S-P) residuals from phases traveling along event-to-station paths are about 7 s smaller at the shield station than at the oceanic stations. This correlation with surface tectonic environments is equally strong for the (ScS-ScP) residuals, with the shield/oceanic station difference being about 4 s. Moreover, the data suggest that this correlation between differential residuals and surface tectonic environments is caused by variations in shear velocity within the upper part of the mantle. However, the data cannot uniquely resolve the required depth of these variations within the mantle. For example, if the shear velocity variations extend to a depth of 400 km beneath the recording stations, then the average shear velocity difference between shield- and oceanic-type environments is about 4%. However, if the variations extend only to a depth of 200 km, this difference is more than 8%.(ScP-PcP) and (ScS-PcS) residuals vary from about +1 to about +4 s at the different stations, apparently because of compressional velocity variations in the mantle along the Pc path. If the variation in compressional velocity within the mantle below a depth of about 600 km is about 10% and occurs near the source region, these results suggest that, in the vicinity of deep earthquake zones, variations in compressional velocity extend to a depth of about 1000 km. However, these results can equally be explained by a 1% variation in compressional velocity, evenly distributed along the entire Pc path. An estimate of Q determined from the observed predominant frequency of ScS waves, as recorded at the shield station, suggests that the average 〈Q_s〉 of the mantle beneath about 600 km is about 1050 at frequencies of about 1 Hz.     

Mapping of upper mantle structure from differential (PP — P) travel-time residuals

P and PP travel-time residuals have been collected from 23 quarters of the regional catalog of earthquakes of the International Seismological Centre (ISC). The difference in travel-time residuals between PP and P are used as an indicator of the P-delay time near the PP reflection point. Global maps of these differential travel-time residuals are used to study lateral heterogeneity in the crust and upper mantle. Two of the regions studied have mean differential travel-time residuals significantly different (at greater than the 95% confidence level) from the global mean of the differential travel-time residuals. In the central Pacific, southeast of Hawaii, differential PP — P travel-time residuals of ? 3.8 ± 0.36 s are observed. These observations imply an average P-delay time of ? 1.1 ± 0.38 s. Average P-delay times of 0.67 ± 0.20 s are observed beneath the Tibetan Plateau and the western portion of the People's Republic of China. These travel-time variations basically agree in location and sign with the surface wave results of Tanimoto and Anderson.     

利用中国地震台网(CSN)记录到的PKP波研究中太平洋下地幔底部的小尺度P波速度变化

PKP震相包含了下地幔底部P波速度结构的重要信息.中国地震台网（CSN）台站记录到的南美洲地震的PKP波的射线,对中太平洋下的D″层有很好的采样.本研究采用这些PKP波的AB和DF两个分支的走时差,研究了中太平洋下地幔底部P波速度的小尺度变化.AB DF的走时差减小了上地幔横向不均匀性的影响,而对下地幔底部P波速度的横向变化十分敏感.与此同时,AB DF的走时差也减小了地震定位误差的影响,消除地震发震时间测定误差的影响.本研究的结果表明,在中太平洋的地幔底部存在着大范围的AB DF走时正残差,也即低速异常区,这可能是太平洋下超大地幔热柱的源处.观测到的P波速度异常的空间分布总体上与Grand通过层析成像得到的CMB的S波的速度异常相一致,并在变化的幅度上很好地相关,P波速度的扰动值（在D″层大约为2%）是Grand 速度模型中的S波异常的36%.这一结果有助于太平洋下超大地幔热柱的结构和性质的进一步研究.     

Weak velocity anomaly in the Earth’s outer core from seismic data

Experimental data on the differential travel time t _BC–t _DF of seismic waves PKP_DF and PKP_BC in the Earth’s core under Africa and Australia are analyzed. The differential travel-time residuals beneath Africa in a narrow range of angles from 21° to 25° between the direction of the seismic ray in the core and the Earth’s rotation axis exhibit a scoop-shaped peculiarity not accounted for by cylindrical anisotropy in the inner core. A model with a 0.2–0.8% P-wave velocity anomaly with a radius of 1375 km in the cylindrical region in the outer core is proposed, which closely fits the experimental data. We suggest that the velocity anomaly is generated by the dynamical processes occurring in the outer core, namely, the growth of the inner core and the convection in the outer core, both leading to the formation of a low-density anomaly in the outer core.     

The effect of the descending lithosphere beneath the Tonga island arc on P-wave travel-time residuals at the Warramunga Seismic Array

P-wave travel-time residuals at the Warramunga Seismic Array (WRA) in the Northern Territory, Australia, have been studied from 49 earthquakes with epicenters south of 19°S in the Fiji-Tonga region. Focal depths are between 42 and 679 km as determined from pP-P. Using the Jeffreys-Bullen and the Herrin travel-time tables the epicentral parameters have been redetermined by considering only “normal” seismic stations in the location procedure. These are those stations where P-wave travel times are probably not affected by lateral heterogeneities caused by the lithosphere descending beneath the Tonga trench. Epicenters of deep earthquakes below 300 km have been relocated by using stations at Δ > 25° only. Epicenters from shallower-depth earthquakes have been recalculated without using stations between 35 < Δ < 75° epicentral distance. In both cases focal depths were determined from pP-P times. The resulting pattern of P-residuals at WRA does not show any significant change with depth below 350 km. The residuals become more negative for shallower earthquakes above about 250 km. P-waves to WRA are advanced by approximately 2 s compared with those from deep earthquakes. The results do not essentially differ for the two different travel-time tables used. The observations can be interpreted by P-wave velocities that are higher in the sinking slab down to 350–400 km by 5±2% than in both the Jeffreys-Bullen and Herrin models. Without considering possible elevations of phase boundaries this estimate yields a temperature contrast of 1000±450°C between slab and normal mantle material in this depth range.     

远震反演中新的地壳异常体改正方法

近年来，远震走时层析成像方法有了长足的发展，在地下结构反演的研究中获得了众多成果。针对射线在台站下方覆盖率较差而导致远震反演方法对地壳速度约束不足的问题，本文提出了一种新的地壳异常体改正方法——对同一台站的相对走时残差进行求和平均去均值，消除了地壳中复杂的速度异常体对上地幔速度结构反演的干扰，反演结果具有更优的相对走时残差分布，同时反演得到的速度模型具有更小的数据方差。     

Travel time tomography inversion and the deep structure of the southern Tibet

ＴｒａｖｅｌｔｉｍｅｔｏｍｏｇｒａｐｈｙｉｎｖｅｒｓｉｏｎａｎｄｔｈｅｄｅｅｐｓｔｒｕｃｔｕｒｅｏｆｔｈｅｓｏｕｔｈｅｒｎＴｉｂｅｔＱＩＮＧ－ＴＩＡＮＬＵ１）（吕庆田）,ＭＥＩ－ＪＩＡＮＧ１）（姜枚）,ＫＡＩ－ＹＩＭＡ１）（马开义）,Ａ．Ｈｉｒｎ２...     

Regional travel-time and slowness measurements for P-waves in the distance range 40–57°

A method for analysing travel times measured at a large array or a network of seismographs from many earthquakes within a specific region has been developed. Approximate relative station corrections are calculated from the residuals on a least-squares line or least-squares quadratic form fitted through the times for each earthquake, and may be improved by iteration after a preliminary travel-time curve has been derived. Accurate relative baseline corrections for each earthquake are also calculated iteratively, and an optimum slowness-distance curve is determined from the combined corrected travel times from all earthquakes using a trade-off procedure. Calculations using synthetic travel-time data suggest that abrupt changes in slowness of ～ 0.4 s deg^?1 due to the presence of triplications are generally resolvable, provided that the effects of lateral variations are small, even with random epicentre mislocations in the range ± 0.5°. Slowness measurements at a network of temporary stations deployed across Australia do not show any discontinuities in slowness greater than 0.2 s deg^?1 in the distance range 45–54°. Similar measurements at the Warramunga array from the same source regions, however, suggest the presence of complexity in the slowness curve at distances close to 50°. Relative arrival times at the temporary network generally have standard deviations less than 0.25 s, thus suggesting that details of structure finer than those derived from conventional travel-time studies can be resolved.     

北京地震台地震走时表的研制

长期地震观测发现,北京地震台实际地震记录的震相走时与全球地震走时表IASP91之间存在走时偏差。结合北京地震台及周边地下40km精细速度结构的研究成果和AK135模型,构建其地震走时表的计算模型。通过对该地震台2003--2008年NCDSN地震观测实际走时比对,对模型结构进行修正,最终建立北京地震台地震走时表。     

Seismic detectability of anomalous structure at the top of the Earth’s outer core with broadband array analysis of SmKS phases

I have examined precisely the differential travel times and waveforms of SmKS seismic phases propagating under the southern Indian Ocean obtained from African broadband seismic arrays. The SmKS phases analyzed in this study travel in the mantle with weak heterogeneity confirmed by a global tomographic study for the distance range of 115-135°. The SmKS differential times were obtained from a vespagram (a stack intensity on a time-slowness diagram), and comparison with the vespagram created from synthetic waveforms with PREM gives the travel-time residual for each event-array pair. Although the residuals of S3KS-S2KS times exhibit apparently a systematic dependence on epicentral distance, this is likely due to small-scale heterogeneity beneath the Oceania where is covered by the SmKS ray entering points at the CMB. Waveform modeling was applied to a record section with a small travel-time residual that suggests a small effect from the mantle heterogeneity on the data set, I found that a low-velocity zone in the outermost 50 km in the core rather than PREM can explain an additional arrival detected just after the S3KS phase. This result is still inconclusive because of the small number of data and non-uniqueness of the model and ambiguity due to mantle structure. However, accumulation of the precise measurement described in this study may help the reduction of uncertainty and trade-offs.     

The inversion of 3-D crustal structure and hypocenter location in the Beijing-Tianjin-Tangshan-Zhangjiakou area by genetic algorithm

Ｔｈｅｉｎｖｅｒｓｉｏｎｏｆ３┐ＤｃｒｕｓｔａｌｓｔｒｕｃｔｕｒｅａｎｄｈｙｐｏｃｅｎｔｅｒｌｏｃａｔｉｏｎｉｎｔｈｅＢｅｉｊｉｎｇ┐Ｔｉａｎｊｉｎ┐Ｔａｎｇｓｈａｎ┐Ｚｈａｎｇｊｉａｋｏｕａｒｅａｂｙｇｅｎｅｔｉｃａｌｇｏ┐ｒｉｔｈｍＹＯＮＧ－Ｇ...     

Elevation of the olivine-spinel transition in subducted lithosphere: Seismic evidence

The top of the olivine-spinel phase change in subducted oceanic lithosphere can be located by the travel times of seismic waves which have propagated through the slab. P-wave travel-time residuals from deep earthquakes in the Tonga island are observed at Australian seismic stations are grouped according to the depth of the earthquake. The change in mean residual with a change in earthquake depth is related to the velocity contrast between slab and normal mantle at that depth. The curve mean residual versus earthquake depth displays a region of markedly increased slope between earthquake depths of about 250 and 350 km. The most probable explanation of this observation is an elevation by 100 km of the olivine-spinel phase change within the relatively cooler slab. No evidence was found for vertical displacements within the slab of any deeper phase changes.A temperature contrast between slab and normal mantle of about 1,000°C at 250 km depth is implied. This finding confirms current thermal models for subducted lithosphere but is inconsistent with the global intraplate stress field unless only a few percent of the negative buoyancy force at subduction zones is transmitted to the surface plates.     

Short scale heterogeneity in the lowermost mantle: insights from PcP-P and ScS-S data

We compare lateral variations at the base of the mantle as inferred from a global dataset of PcP-P travel time residuals, measured on broadband records, and existing P and S tomographic velocity models, as well as ScS-S travel time data in some selected regions. In many regions, the PcP-P dataset implies short scale lateral variations that are not resolved by global tomographic models, except under eastern Eurasia, where data and models describe a broad region of fast velocity anomalies across which variations appear to be of thermal origin. In other regions, such as central America and southeastern Africa, correlated short scale lateral variations (several hundred kilometers) are observed in PcP and ScS, implying large but not excessive values for the ratio R=∂ ln V_s/∂ ln V_p (∼2.5). On the other hand, in at least two instances, in the heart of the African Plume and on the edge of the Pacific Plume, variations in P and S velocities appear to be incompatible, implying strong lateral gradients across compositionally different domains, possibly also involving topography on the core-mantle boundary. One should be cautious in estimating R at the base of the mantle from global datasets, as different smoothing and sampling of P and S datasets may result in strong biases and meaningless results.     

福建及台湾海峡地区新一维速度模型的适用性分析

通过福建及台湾海峡地区的新一维速度模型与现有华南速度模型的对比,讨论了新一维速度模型在福建地震观测台网的适用性。理论走时分析结果表明,尽管两个速度模型差异明显,但震中距在0—100 km范围内的震相理论走时相差较小,一定程度上说明两速度模型所给出的本区域地壳平均速度差异较小。对利用18次人工定点爆破记录的地震定位结果的分析表明:当震源深度不受约束时,应用华南速度模型的定位结果精度稍优于新一维速度模型;将震源深度固定为0 km后,应用新一维速度模型的定位结果精度则明显优于华南模型。对19个仙游震群序列事件进行定位的结果显示,由于华南地区速度结构的横向变化较小,应用两模型的地震定位精度结果基本相当,但新一维速度模型定位的发震时刻较华南速度模型普遍早0.61 s左右,因此使得事件定位残差显著增大。     

Theoretical study on the influence of CMB topography on the core reflection ScS

Computing synthetic seismograms for media with localized heterogeneous regions can be performed using hybrid methods. Here, a combination of a finite-difference (FD) technique and a frequency-wavenumber (ω − k) filtering is applied to model wave reflection at different kinds of core-mantle boundary (CMB) topography. The FD method is only applied in the neighbourhood of the CMB, while the ω − k filter is used to continue the reflected wavefield to the Earth's surface. Synthetic SH-seismograms for ScS with a dominant frequency of 0.5 Hz are computed at epicentral distances from 44° to 69°. The topography varies in amplitude (maximum amplitude of 1.0–2.7 km) and in its wavenumber spectrum; it is either monochromatic (wavelengths from 55 to 270 km) or statistical (coloured noise). The seismograms for a CMB with topography are compared with those for a plane CMB. We observe that monochromatic topography with short wavelengths (less than 100 km) results in amplitude reduction and shorter travel times than in the case of a plane CMB, but no variations with epicentral distance appear, whereas greater wavelengths exhibit amplitude variations with distance as well as travel time residuals, which both correlate with the CMB topography. Statistical models show amplitude variations with epicentral distance, while the travel time residuals are very small (less than 0.1 s). All synthetics illustrate that wavefront healing occurs along the ray path from the CMB to the Earth's surface. While the seismograms at the CMB exhibit strong fluctuations, the fluctuations at the surface are smoothed and reduced. This demonstrates that it is necessary to use wave theoretical methods for computing synthetic seismograms for complicated structures at greater depth. It also follows that travel times are less sensitive to the structure than the amplitudes.     

Crustal transit time monitoring using PKP

The arrival time difference for the AB branch of PKP from deep Tonga earthquakes is constant over years with a standard deviation of ±0.05 seconds at seismographs located 10 to 50 km from each other. If published travel time curves are used to calculate the relative residuals of PKP the standard deviation from the constant mean is improved by approximately 0.01 seconds for AB branch data. For the BC branch, standard deviations of relative travel times of ±0.06 seconds are reduced to less than ±0.05 seconds by calculating relative residuals. We conclude that changes of crustal transit time forP-waves could be resolved, based on careful PKP arrival time measurement at two or more neighboring stations if the changes exceed 0.05 sec and last for more than one year. The conditions for achieving this result are that PKP from Tonga earthquakes is clearly recorded, and that time-keeping is accurate. The data on which these conclusions are based were obtained from the Graefenberg seismograph array, which is located in West Gemany and consists of 13 stations separated by distances of 10 km to 100 km. We propose that relative arrival times of PKP from Tonga could be used in the Mediterranean - Middle East area to search for precursory travel time changes before large earthquakes.     

Mantle lithosphere,asthenosphere and transition zone beneath Eastern Anatolia

Journal of Seismology - By using P and S wave receiver functions and P and S wave travel time residuals, we have found velocity models for 16 seismograph stations in Eastern Anatolia. Our study is...     

On the sensitivity of tidal residuals off the west coast of Britain to mesh resolution

An irregular mesh model of the west coast of Britain is used to examine the sensitivity of tidal residuals to mesh resolution in the region. Computed residuals are compared with earlier published results determined with a high resolution (1 km grid) finite difference model of the eastern Irish Sea. Initial calculations show that tidal residuals are largest in nearshore regions particularly in the vicinity of headlands. Local refinement of the mesh in these regions leads to a more detailed picture of the flow field, particularly adjacent to the coast. Although large scale offshore features of the flow can be resolved using the high resolution finite difference model, such an approach leads to a “stair case” representation of the coastal boundary with an adjacent near coastal region of spurious tidal residuals. By using an irregular mesh that follows the coast, this effect is removed. In the Mersey river region the tidal residual is resolved with a mesh resolution of 120 m, although calculations show that its distribution is particularly sensitive to small scale features of the topography. A variable mesh that can accurately represent the lateral variations in river width and details of topography in both the nearshore and estuarine environment appears essential in modelling the coastal spread of freshwater plumes from rivers and pollutants discharged into the near coastal environment.