The 1998 Papua New Guinea Earthquake and its Fault Plane Estimated from Relocated Aftershocks

— The 1998 Papua New Guinea earthquake of M _w 7.0 occurred near the Wewak trench where the North Bismarck plate is subducting beneath the Australian plate. Its mechanism is thrust-type, and one of the nodal planes is almost parallel to the plate interface. To determine which of the two nodal planes of the main shock is the fault plane, we relocated the main shock and aftershocks using a method of modified joint hypocenter determination. We combined and employed two types of data in this study. Firstly, we used data reported by the National Earthquake Information Center (NEIC) of the U.S. Geological Survey (USGS), which includes three stations at the northeastern edge of Irian Jaya and one station in northern Papua New Guinea, from which the epicentral distances are less than 2 degrees. Secondly, in addition to the above permanent-station data, we used data from temporary aftershock observations near the epicentral area around the Sissano Lagoon carried out by Tsuji et al. (1998). Using three-component seismometers, they carried out observations from August 2 to October 2, 1998 at three sites. Although the network did not record the main shock and immediate aftershocks, the data obtained by temporary observation sites can clearly assist in identifying their absolute locations, since it is possible to apply the joint hypocenter determination (JHD) method. Hypocenters were relocated between the coastline and the Wewak trench, distributed along a nodal plane dipping shallowly to the southwest. Therefore, we can conclude that this nodal plane is the main shock fault and that the 1998 Papua New Guinea earthquake was an interplate earthquake between the North Bismarck and Australian plates.     

类新西兰地震波形分析方法

对泰安地震台记录的2010年9月4日新西兰地震波形,进行震相、震源深度等地震参数辨析和分析,以具有代表性的2009年7月15日新西兰南岛西海岸远海地震波形和2010年4月5日墨西哥地震波形分析为辅,提出震中距大于83°、小于105°的类新西兰地震分析方法.根据两个辅震的波形特性,结合中国地震台网中心定位结果,对震级误差...     

内蒙古区域地震台网监测能力研究

国际上新近发展的“基于概率的完整性震级”(PMC)方法,具有可考察地震定位中由于台站人为选择等造成的台网监测能力下降,以及避免传统基于G-R关系的统计算法因地震数目过少而无法评估等优点.本研究利用PMC方法,计算得到内蒙古区域地震台网39个台站对周边地震事件的检测概率及台网检测概率.单台检测概率结果显示:PMC方法能够客观地反映39个台站对地震事件的检测能力;因台网布局等影响,内蒙古区域地震台网中西部和中东部地区的台站检测能力较强,而靠近蒙古、俄罗斯边境的台站, 阿拉善右旗附近地区的台站,以及邻近吉林、黑龙江等地区的台站检测能力较低.合成检测概率结果显示,由于邻省台站的引入,全区80%的地区基于概率的最小完整性震级M_P达到2.2左右,其余地区M_P达到3.3左右.为提高地震台网监测能力,建议在监测能力较弱的中蒙交界地区、东北部地区,以及阿拉善左旗以西地区适度加密台站,进一步优化台网布局.      

Subducted slabs beneath the eastern Indonesia-Tonga region: insights from tomography

Tomographic images of mantle structure beneath the region north and northeast of Australia show a number of anomalously fast regions. These are interpreted using a recent plate tectonic reconstruction in terms of current and former subduction systems. Several strong anomalies are related to current subduction. The inferred slab lengths and positions are consistent with Neogene subduction beneath the New Britain and Halmahera arcs, and at the Tonga and the New Hebrides trenches where there has been rapid rollback of subduction hinges since about 10 Ma. There are several deeper flat-lying anomalies which are not related to present subduction and we interpret them as former subduction zones overridden by Australia since 25 Ma. Beneath the Bird’s Head and Arafura Sea is an anomaly interpreted to be due to north-dipping subduction beneath the Philippines-Halmahera arc between 45 and 25 Ma. A very large anomaly extending from the Papuan peninsula to the New Hebrides, and from the Solomon Islands to the east Australian margin, is interpreted to be the remnant of south-dipping subduction beneath the Melanesian arc between 45 and 25 Ma. This interpretation implies that a flat-lying slab can survive for many tens of millions of years at the bottom of the upper mantle. In the lower mantle there is a huge anomaly beneath the Gulf of Carpentaria and east Papua New Guinea. This is located above the position where the tectonic model interprets a change in polarity of subduction from north-dipping to south-dipping between 45 and 25 Ma. We suggest this deep anomaly may be a slab subducted beneath eastern Australian during the Cretaceous, or subducted north of Australia during the Cenozoic before 45 Ma. The tomography also supports the tectonic interpretation which suggests little Neogene subduction beneath western New Guinea since no slab is imaged south of the New Guinea trench. However, one subduction zone in the tectonic model and many others, that associated with the Trobriand trough east of Papua New Guinea and the Miocene Maramuni arc, is not seen in the tomographic images and may require reconsideration of currently accepted tectonic interpretations.     

2015年3月新不列颠MS7.4地震震源及邻区构造应力场特征

2015年3月30日至5月15日,巴布亚新几内亚-新不列颠地区发生了一系列地震.为研究该地区的构造应力环境及孕震背景,本文基于Global CMT目录,对新不列颠区域浅部进行构造应力场反演,拟得到高精度的应力图像.反演结果显示：（1）沿着南、北俾斯麦块体边界的区域构造应力场呈走滑体系,最大主压应力轴方位呈SWW-NEE向.（2）所罗门海的NW和NE走向的海沟处于压缩状态,所罗门海块体向新不列颠和所罗门群岛俯冲的板块弯曲部分是局部拉张.（3）受俯冲带的北向推挤,南俾斯麦板块顺时针旋转的挤压,太平洋板块向西部运动汇聚作用,新不列颠岛东北部与新爱尔兰岛南部交汇区域呈现明显非均匀应力状态.（4）此次地震序列的大多数走滑型和逆冲型地震,可能是所罗门海块体俯冲运动,和南俾斯麦块体与太平洋板块的近EW向挤压作用共同引发.     

Seismicity of the Indian Peninsular Shield from Regional Earthquake Data

--In this comprehensive study of seismicity and seismotectonics of the peninsular Indian shield region, seismic data of regional earthquakes spanning two decades (1978-1997), obtained at Gauribidanur seismic array (India) and integrated where necessary with data from other seismological stations in the region, have been analyzed in detail. With a slow rate of stress accumulation, the shield is found to have low to moderate seismicity that takes into account a couple of earthquakes of magnitude slightly larger than 6. The frequency-magnitude analysis of the data set gives a b value of 1.18. The spatio-temporal pattern of occurrences of the earthquakes combined with their magnitude and seismic energy distribution is consistent with the view that the peninsular seismicity is low to moderate and episodic in nature. Regions of moderate seismicity and its low-grade counterpart constituted by microearthquakes (magnitude less than 3), appear correlated to the areas traversed by known geologic faults and subfaults, shear zones, and other such tectonic features. Microearthquakes represent about two-thirds of the total regional seismic events during the past two decades.     

区域数字地震台网近距离面波震级MS的测定

使用福建区域数字台网7个宽带台观测到的47个发生于台湾地区的浅源强震(M_S ≥ 5)的近距离面波记录,量算了水平向乐夫波的最大合成振幅及震中距离,测定了所有这些强震的近距离面波震级MS。并将结果与国家台网公布的远距离平均MS震级进行比较和统计分析。结果表明,近距离测定的面波震级比远距离测定的M_S 有较大偏差。震级越小、距离越近,差异越大。对于区域地震台网测定的M_S 震级,需要做必要的修正,才能使其测定结果尽可能接近于国家台网公布的平均M_S 震级。文中还讨论了近距离和远距离面波震级测定差异的原因,及其修正办法。可供区域数字宽带台网速报本区域发生的近距离强震面波震级时参考。     

中国地震台网监测能力评估和台站检测能力评分(2008—2015年)

针对中国地震台网"十五"项目建成后的地震监测能力科学评估的需求,为进一步优化台网布局、提升边疆海域等重点地区监测能力,本研究利用"基于概率的完整性震级"（PMC）方法,对中国地震台网1001个台站以及2008-10-01-2015-09-17期间实际产出的地震观测资料进行了研究,分析了指定震级档下的检测概率P_E和最小完整性震级M_P的分布.除台网整体监测能力分布外,为直观地用单分值表述逐个台站的地震检测能力,本文发展了基于等振幅曲线的"地震检测能力评分表",给出了国家台和区域台每个台站的地震检测能力评分统计特征和空间分布特征.此外,研究中还采用设定"最佳"地震监测能力目标函数的方式,模拟了通过改进观测条件可获得的地震台网监测能力提升的理论结果.研究结果表明,我国华北和东南沿海等东部地区地震监测能力较高,西部尤其是青藏高原南部地区M_p仅约为4.5,近海海域M_p仅约为3.5;从单个地震台站的运行效益角度,台网运行水平和地震观测资料的分析程度对台站的实际的地震检测能力影响显著,新疆等部分台站稀疏地区地震检测能力较高,而中等台站密度的贵州等部分区域相对较低;国家台的地震检测能力评分D_score系统优于区域台,新疆等西部边疆地区,以及福建等东南沿海地区的D_score明显高于台站密集的东部地区;模拟结果显示,在现有台站布局条件下,通过台站优化改造和提升运维管理水平,可显著提升对内蒙古西部、四川西部、甘肃-青海的北部交界地区、鄂尔多斯地块内部、贵州大部分地区,以及我国近海海域、朝鲜半岛北部和中南半岛北部地区的地震监测能力.     

New Zealand Earthquake Forecast Testing Centre

The New Zealand Earthquake Forecast Testing Centre is being established as one of several similar regional testing centres under the umbrella of the Collaboratory for the Study of Earthquake Predictability (CSEP). The Centre aims to encourage the development of testable models of time-varying earthquake occurrence in the New Zealand region, and to conduct verifiable prospective tests of their performance over a period of five or more years. The test region, data-collection region and requirements for testing are described herein. Models must specify in advance the expected number of earthquakes with epicentral depths h ≤ 40 km in bins of time, magnitude and location within the test region. Short-term models will be tested using 24-h time bins at magnitude M ≥ 4. Intermediate-term models and long-term models will be tested at M ≥ 5 using 3-month, 6-month and 5-year bins, respectively. The tests applied will be the same as at other CSEP testing centres: the so-called N test of the total number of earthquakes expected over the test period; the L test of the likelihood of the earthquake catalogue under the model; and the R test of the ratio of the likelihoods under alternative models. Four long-term, three intermediate-term and two short-term models have been installed to date in the testing centre, with tests of these models commencing on the New Zealand earthquake catalogue from the beginning of 2008. Submission of models is open to researchers worldwide. New models can be submitted at any time. The New Zealand testing centre makes extensive use of software produced by the CSEP testing centre in California. It is envisaged that, in time, the scope of the testing centre will be expanded to include new testing methods and differently-specified models, nonetheless that the New Zealand testing centre will develop in parallel with other regional testing centres through the CSEP international collaborative process.     

Large-scale volcanic cone collapse: The 1888 slope failure of Ritter volcano,and other examples from Papua New Guinea

A largely submarine avalanche amphitheatre that formed catastrophically in 1888 on Ritter volcano has been identified from a bathymetric survey. Collapse of the volcano in 1888 therefore is considered to have been caused by rapid, large-scale slope failure, rather than by cauldron subsidence, as previously supposed. Escarpments of pre-historic slope failures are common on other Papua New Guinea volcanoes. Directions of avalanching on some volcanoes in the Bismarck volcanic arc appear to be controlled by a regional stress pattern, and those for some volcanoes in the Fly-Highlands province on mainland Papua New Guinea point away from the regional centre of Pliocene uplift. Large amphitheatres such as at Doma Peaks in the Fly-High-lands province probably originated by multiple collapses.     

Observation and Modeling of the January 2009 West Papua, Indonesia Tsunami

We modeled a tsunami from the West Papua, Indonesia earthquakes on January 3, 2009 (M _w?=?7.7). After the first earthquake, tsunami alerts were issued in Indonesia and Japan. The tsunami was recorded at many stations located in and around the Pacific Ocean. In particular, at Kushimoto on Kii Peninsula, the maximum amplitude was 43?cm, larger than that at Manokwari on New Guinea Island, near the epicenter. The tsunami was recorded on near-shore wave gauges, offshore GPS sensors and deep-sea bottom pressure sensors. We have collected more than 150 records and used 72 stations?? data with clear tsunami signals for the tsunami source modeling. We assumed two fault models (single fault and five subfaults) which are located to cover the aftershock area. The estimated average slip on the single fault model (80?×?40?km) is 0.64?m, which yields a seismic moment of 1.02?×?10²⁰?Nm (M _w?=?7.3). The observed tsunami waveforms at most stations are well explained by this model.     

Eyewitness Accounts of the Impact of the 1998 Aitape Tsunami, and of Other Tsunamis in Living Memory, in the Region from Jayapura, Indonesia, to Vanimo, Papua New Guinea

Field investigations in 1999 confirmed that the tsunami that struck the Aitape coast of Papua New Guinea on 17 July, 1998 caused damage at points as far as 230 km to the west-northwest, particularly at locations where the coast is indented. Eyewitnesses saw the sea withdraw (in most cases), then surge to levels around 2 m higher than normal in a series of three waves. In some cases the time of arrival of the waves is known approximately by reference to the onset of darkness and to felt earthquakes. Seiche waves followed in some bays, notably in Yos Sudarso Bay, Indonesia, where waves persisted for 3–5 days. Damage was caused by the backwash from the waves. Bodies presumed to be those of Aitape victims were seen floating at sea off Jayapura five days after the tsunami. We record the recollections of people in the Yos Sudarso Bay area who experienced a number of tsunamis in the past 60 years; people that we interviewed on the Papua New Guinea side of the border recollected few or none.     

极值理论在福建沿海地区地震预测中的研究与应用

利用陈培善等人对极值理论修改后的极值分布函数模型对1971年1月-2012年9月福建沿海地区的地震资料进行统计分析,单位时间内最小与最大地震的确定分别考虑相关区域内台网的监控能力以及删除余震后的实际情况.根据修正后的极值理论统计计算出该地区地震的复发周期及截至2015年12月可能发生的相应地震的次数与发震概率,并对有关结果进行映震能力分析,为今后福建沿海地区地震趋势分析提供可遵循的理论依据.     

美国新马德里地震带的震间应变积累

使用全球定位系统（ＧＰＳ）卫星接收仪在美国新马德里地震带重新测量了历史上的三角网，并据此估算了震间应变速率．使用的观测资料包括５０年代初期在这个地区２０多个测站上获取的三角测量数据，以及在这些测站上于１９９１年和１９９３年所取得的ＧＰＳ测量数据．在整个网内的平均应变速率为每年０．１４１±０．０２５微应变．最大主压缩轴方向为ＮＥ８９°±６°．这一实测平均应变速率的数量级与考古地震研究和地震活动速率关于大地震重复周期在新马德里地区大约为５００－１０００ａ的估计相吻合．在美国中部进行的原地应力测量指出，该地区的最大主任应力轴为ＮＥＥ向，本文关于最大主压缩轴方向的结果与此相－致．     

震级国家新标准在内蒙古测震台网速报中的应用

选取内蒙古测震台网2016年1月—2018年9月速报的地震,按照新的国家标准《地震震级的规定》(GB17740—2017)和现行地震速报规定分别重新测定震级,并测定矩震级。结果表明:新国标测定的地震数可较好地满足G—R线性关系,且所测震级更接近矩震级。震级新的国家标准应用后,速报的震级下限降低,地震数明显增加,可有效满足社会发展对有感地震监测服务的需求,并更加科学地为地震预报研究提供基础数据资料。     

A preliminary investigation of strong-motion data from the French Antilles

Strong-motion networks have been operating in the Caribbean region since the 1970s, however, until the mid-1990s only a few analogue stations were operational and the quantity of data recorded was very low. Since the mid-1990s, digital accelerometric networks have been established on islands within the region. At present there are thought to be about 160 stations operating in this region with a handful on Cuba, 65 on the French Antilles (mainly Guadeloupe and Martinique), eight on Jamaica, 78 on Puerto Rico (plus others on adjacent islands) and four on Trinidad.After briefly summarising the available data from the Caribbean islands, this article is mainly concerned with analysing the data that has been recorded by the networks operating on the French Antilles in terms of their distribution with respect to magnitude, source-to-site distance, focal depth and event type; site effects at certain stations; and also with respect to their predictability by ground motion estimation equations developed using data from different regions of the world. More than 300 good quality triaxial acceleration time-histories have been recorded on Guadeloupe and Martinique at a large number of stations from earthquakes with magnitudes larger than 4.8, however, most of the records are from considerable source-to-site distances. From the data available it is found that many of the commonly-used ground motion estimation equations for shallow crustal earthquakes poorly estimate the observed ground motions on the two islands; ground motions on Guadeloupe and Martinique have smaller amplitudes and are more variable than expected. This difference could be due to regional dependence of ground motions because of, for example, differing tectonics or crustal structures or because the ground motions so far recorded are, in general, from smaller earthquakes and greater distances than the range of applicability of the investigated equations.     

Neo-Deterministic Seismic Hazard and Pattern Recognition Techniques: Time-Dependent Scenarios for North-Eastern Italy

An integrated neo-deterministic approach to seismic hazard assessment has been developed that combines different pattern recognition techniques, designed for the space?Ctime identification of impending strong earthquakes, with algorithms for the realistic modeling of seismic ground motion. The integrated approach allows for a time-dependent definition of the seismic input, through the routine updating of earthquake predictions. The scenarios of expected ground motion, associated with the alarmed areas, are defined by means of full waveform modeling. A set of neo-deterministic scenarios of ground motion is defined at regional and local scales, thus providing a prioritization tool for timely preparedness and mitigation actions. Constraints about the space and time of occurrence of the impending strong earthquakes are provided by three formally defined and globally tested algorithms, which have been developed according to a pattern recognition scheme. Two algorithms, namely CN and M8, are routinely used for intermediate-term middle-range earthquake predictions, while a third algorithm allows for the identification of the areas prone to large events. These independent procedures have been combined to better constrain the alarmed area. The pattern recognition of earthquake-prone areas does not belong to the family of earthquake prediction algorithms since it does not provide any information about the time of occurrence of the expected earthquakes. Nevertheless, it can be considered as the termless zero-approximation, which restrains the alerted areas (e.g. defined by CN or M8) to the more precise potential location of large events. Italy is the only region of moderate seismic activity where the two different prediction algorithms, CN and M8S (i.e. a spatially stabilized variant of M8), are applied simultaneously and a real-time test of predictions, for earthquakes with magnitude larger than a given threshold (namely 5.4 and 5.6 for CN algorithm, and 5.5 for M8S algorithm) has been ongoing since 2003. The application of the CN to the Adriatic region, which is relevant for seismic hazard assessment in the northeastern part of the Italian territory, is also discussed. Examples of neo-deterministic scenarios are provided, at regional and local scale and for the cities of Trieste and Nimis (Friuli Venezia Giulia region), where the knowledge of the local geological conditions permitted a detailed evaluation of the expected ground motion.     

河北易县M_L3.8、3.3地震同震形变响应及其震源机制

2010年7月30日河北省易县相继发生ML3.8、3.3地震,震中距离易县地震台仅7 km,震中附近主要以形变观测台站为主.虽然地震震级相对较小,但是12个形变测项均有不同程度的同震响应,这在中小地震的同震响应分析中是较为突出的震例.根据河北台网的地震波形资料对地震进行了重新定位.根据速报震中和重新定位震中分别计算两次地震的震源机制,震源断层初步断定为NEE向断层,断层表现为正断兼右旋走滑.     

TOPOGRAPHIC AND ISOSTATIC CORRECTIONS TO GRAVITY SURVEYS IN MOUNTAINOUS AREAS*

In mountainous areas it is an inadequate procedure to reduce gravity observations by merely subtracting the effect of an infinite flat slab of material between the station and sea-level, and adding a terrain correction. A programme is described which directly computes the effect of masses above sea-level, and mass-deficiencies below it; as well as the effect of compensating masses under the Pratt and Airy isostatic systems. As an example, the method has been applied to a regional gravity survey of Papua and New Guinea where it is seen to remove the usually high correlation of the Bouguer anomaly with local topography.     

广西测震台网台站地方性震级偏差分析

选用广西数字化测震台网2008—2020年记录的广西及邻区1 251个M_L ≥ 2.0地震事件震相资料,采用震级残差统计方法,获取广西51个地震台单台震级与台网平均震级的平均偏差和标准差,分析观测仪器类型、台基场地响应和震中距等因素对震级偏差的影响,进而提出台基校正值,修正区域量规函数。结果发现,南宁、天等、北海等15个地震台站的震级偏大,桂林、贺州、天峨等16个地震台站的震级偏小,崇左、灵山、平果等20个地震台站震级适中;在单台震级共14 784个样本中,震级偏差集中在-0.4—0.4,呈正态分布;震级偏差受台基场地响应和震中距因素影响明显,当场地响应值大于1.0时震级偏大,小于1.0时震级偏小,等于1.0时震级较为适合,当震中距在70 km范围内时单台震级偏小,当震中距在70—300 km时单台震级较为适合,当震中距在300—550 km时,单台震级稍有偏大。修正后的区域量规函数和台基校正值,可有效提高广西测震台网测定震级的准确性。