南海东部次海盆地震背景噪声分析

背景噪声的强弱是影响地震台站观测的一个重要因素.获取背景噪声的分布特征对评估海底地震仪记录数据质量及对数据的降噪处理均具有重要的指示意义.利用概率密度函数方法获取台站数据的功率谱密度的概率分布特征并与全球背景噪声高值模型和低值模型进行对比是研究台站周围环境背景噪声水平的有利手段.本研究基于南海大规模的被动源海底地震仪台...     

南海北部陆缘OBS2018-H2测线地壳结构初步结果

南海北部陆缘洋陆转换带实施的OBS2018-H2测线的地壳速度结构, 将为探讨南海张裂-破裂机制提供重要证据。文章介绍了OBS2018-H2测线前期数据处理流程, 包括多道反射地震数据处理、海底地震仪OBS (Ocean Bottom Seismometer)数据格式转换、炮点和OBS位置校正, 以及OBS震相的初步识别, 并对地壳结构进行了初步分析。结果表明: 炮点和OBS位置校正效果良好; 多道反射地震数据为建立初始速度模型提供了良好约束; OBS综合地震剖面识别了多组清晰的P波震相, 包括Pw、Pg、PmP和Pn震相。根据测线西侧OBS36、OBS37两台站的震相分布特征初步估算台站下方地壳厚度约为6~7km, 与根据多道地震剖面LW3的双程走时估算的厚度6~9km大致相符。     

Reproducibility of spatial and temporal distribution of aseismic slips in Hyuga-nada of southwest Japan

The Hyuga-nada region of southwest Japan, which is located off the east coast of Kyushu Island, may have the potential to generate great interplate earthquakes along the Nankai trough in the future. In this area, thrust earthquakes of M = 6.7–7.2 have occurred with recurrence intervals of approximately 30 years. In association with these earthquakes, possible local heterogeneities of plate coupling may be expected within 100 km from the coast in the Hyuga-nada region. We investigate numerical experiments to determine the spatial and temporal resolution of slip on the plate interface beneath the Hyuga-nada offshore region. For this purpose, we calculated synthetic displacement data from the result of numerical simulation conducted for the afterslip following an M_w 6.8 earthquake, for existing global positioning system stations on land and planned ocean floor seismic network stations. The spatial and temporal distribution of fault slip is then estimated using a Kalman filter-based inversion. The slip distribution estimated by using ocean floor stations demonstrates that the heterogeneity of plate coupling is resolved approximately within 50 km from the coastal area. This heterogeneity corresponds to the coseismic area of an M_w 6.8 earthquake with a radius of 10 km. Our study quantitatively evaluates the spatial resolution of aseismic slip in the Hyuga-nada region. Analysis based on continuous ocean floor data is useful for resolving the spatial variations of heterogeneities in plate couplings.     

Seismic activity recorded by a single OBS/H near the active Longqi hydrothermal vent at the ultraslow spreading Southwest Indian Ridge (49°39′ E)

ABSTRACT

In this paper, we describe the results of a first pilot experiment of passive seismic monitoring near the active Longqi hydrothermal vent at the Southwest Indian Ridge during Chinese cruise DY115-30. During 102-day experiment, we recorded over 2,000 seismic events on a single ocean bottom seismometer with hydrophone close to the hydrothermal field. We classified these events into the following three classes based on their time–frequency characteristics: (1) volcano tectonic microearthquakes (VTMs); (2) regional earthquakes; (3) short-duration events (SDEs). We preliminarily located the 1,277 VTMs using a single-station location method. Our results revealed hypocenters primarily beneath the axial valley ridge which is about 10?km away from hydrothermal vent and extending to the upper mantle at a depth of 15?km. Two discrete swarms of VTMs occurred during our survey period. The SDEs show a complex frequency characteristics ascribed to clusters of monochromatic waves around different frequencies. This suggests that SDEs result from combination of harmonic waves resonating in fluid-filled pipes and cracks.     

Large aperture seismic imaging at a convergent margin: Techniques and results from the Costa Rica seismogenic zone

In March and April 1995 a cooperative German, Costa Rican, and United States research team recorded onshore-offshore seismic data sets along the Pacific margin of Costa Rica using the R/V Ewing. Off the Nicoya Peninsula we used a linear array of ocean bottom seismometers and hydrophones (OBS/H) with onshore seismometers extending across much of the isthmus. In the central area we deployed an OBS/H areal array consisting of 30 instruments over a 9 km by 35-km area and had land stations on the Nicoya Peninsula adjacent to this marine array and also extending northeast on the main Costa Rican landmass. Our goal in these experiments was to determine the crustal velocity structure along different portions of this convergent margin and to use the dense instrument deployments to create migrated reflection images of the plate boundary zone and the subducting Cocos Plate. Our specific goal in the central area was to determine whether a subducted seamount is present at the location of the 1990, M 7 earthquake off the Nicoya Peninsula and can thus be linked to its nucleation. Subsequently we have processed the data to improve reflection signals, used the data to calculate crustal velocity models, and developed several wide-aperture migration techniques, based on a Kirchhoff algorithm, to produce reflection images. Along the northern transect we used the ocean bottom data to construct a detailed crustal velocity model, but reflections from the plate boundary and top and bottom of the subducting Cocos plate are difficult to identify and have so far produced poor images. In contrast, the land stations along this same transect recorded clear reflections from the top of the subducting plate or plate boundary, within the seismogenic zone, and we have constructed a clear image from this reflector beneath the Nicoya shelf. Data from the 3-D seismic experiment suffer from high-amplitude, coherent noise (arrivals other than reflections), and we have tried many techniques to enhance the signal to noise ratio of reflected arrivals. Due to the noise, an apparent lack of strong reflections from the plate boundary zone, and probable structural complexity, the resulting 3-D images only poorly resolve the top of the subducting Cocos Plate. The images are not able to provide compelling evidence of whether there is a subducting seamount at the 1990 earthquake hypocenter. Our results do show that OBS surveys are capable of creating images of the plate boundary zone and the subducting plate well into the seismogenic zone if coherent reflections are recorded at 1.8 km instrument spacing (2-D) and 5 km inline by 1 km crossline spacing for 3-D acquisition. However, due to typical high amplitude coherent noise, imaging results may be poorer than expected, especially in unfavorable geologic settings such as our 3-D survey area. More effective noise reduction in acquisition, possibly with the use of vertical hydrophone arrays, and in processing, with advanced multiple removal and possibly depth filtering, is required to achieve the desired detailed images of the seismogenic plate boundary zone.     

基于背景噪声互相关技术的海底地震仪的时钟校正方法在西南印度洋脊的应用

Seismic monitoring using ocean bottom seismometers(OBS) is an efficient method for investigating earthquakes in mid-ocean ridge far away from land. Clock synchronization among the OBSs is difficult without direct communication because electromagnetic signals cannot propagate efficiently in water. Time correction can be estimated through global positioning system(GPS) synchronization if clock drift is linear before and after the deployment. However, some OBSs in the experiments at the southwest Indian ridge(SWIR) on the Chinese DY125-34 cruise had not been re-synchronized from GPS after recovery. So we attempted to estimate clock drift between each station pairs using time symmetry analysis(TSA) based on ambient noise cross-correlation. We tested the feasibility of the TSA method by analyzing daily noise cross-correlation functions(NCFs) that extract from the data of another OBS experiment on the Chinese DY125-40 cruise with known clock drift and the same deployment site. The results suggest that the NCFs' travel time of surface wave between any two stations are symmetrical and have an opposite growing direction with the date. The influence of different band-pass filters,different components and different normalized methods was discussed. The TSA method appeared to be optimal for the hydrophone data within the period band of 2–5 s in dozens of km-scale interstation distances. A significant clock drift of ~2 s was estimated between OBSs sets through linear regression during a 108-d deployment on the Chinese cruise DY125-34. Time correction of the OBS by the ambient noise cross-correlation was demonstrated as a practical approach with the appropriate parameters in case of no GPS re-synchronization.     

利用海底地震仪数据分析台风对海底环境噪音的影响

在海底布设的海底地震仪（OBS）能比较清晰地记录到海底的环境噪音,而台风可以直接或间接的产生在海底传播的弹性波,从而影响海底的环境噪音,并在较大程度上影响OBS的数据记录。本文通过分析台风对工作区的整个影响过程中OBS记录数据的振幅变化,再选择合适的滤波方式,首次发现台风产生的风浪及涌浪在短周期海底地震仪的记录数据上有良好的表现特征,指出了台风对海底环境噪音的另一种可能的影响方式,并由此得出：1）台风产生的风浪和涌浪对海底环境噪音的影响模式不同;2）风浪和涌浪所加强的海底环境噪音的范围和程度不同;3）短周期OBS可以比较清晰的记录涌浪信息,其周期主要是6—8 s,且能量稳定（简称“8秒现象”）。这三点结论为后期的海洋地震研究和海洋学其他研究提供经验与借鉴。     

Relocation of two earthquakes in the Southwest Indian Ridge area combining land seismic stations’ with OBSs’ data

Two earthquakes were recorded by 20 ocean bottom seismometers (OBS) deployed in the Southwest Indian Ridge (SWIR) area during a three-dimensional seismic survey in 2010. Their magnitudes (both M _b = 4.4) and hypocenters have been determined by National Earthquake Information Center (NEIC) only using land seismic stations onset times. After the frequency analysis and the band-pass filtering of the OBSs’ data, 7 and 13 P-phase onset times from OBSs were successfully picked for these two events, respectively. Then these two events were relocated by HYPOSAT program with onset times together from OBSs and land seismic stations using different velocity models. These relocation experiments confirm both the importance of adding OBSs’ onset data and the need to apply a local oceanic velocity model for the location of these two events happened on the SWIR. This research has accumulated a wealth of experience for earthquakes observation and research using OBSs in the ocean.     

Evidence that biological activity affects Ocean Bottom Seismograph recordings

Brief and impulsive signals of uncertain origin appear regularly on records from Ocean Bottom Seismographs (OBS) of several institutions. These signals have been recorded on nearly all deployments of the Texas OBS, including sites at depths greater than 7000 m. At some sites, they account for over 90% of the events recorded. They are of short duration (usually 0.5–4.0 s) and have a characteristic frequency (usually in the range of 4–18 Hz) that differs from site to site. When networks of OBS instruments are deployed, the signals are not recorded simultaneously by different instruments. Neither the frequency content nor the distribution of durations of these signals is similar to what is observed for known earthquake events.We present evidence suggesting that the signals are of biological origin, perhaps caused by animals touching the OBS units. (1) The distribution of these signals on instruments deployed at depths shallower than 1000 m shows a 24 h periodicity, while there is a 24 h periodic pattern on instruments deployed at sites deeper than 1000 m (where there is no visible light). (2) The frequency of occurrence of signals is similar to the vertical distribution of biomass in the oceans, i.e., they appear most frequently on OBS instruments deployed at very shallow depths. (3) Biological material has been found attached to several OBS units upon recovery.University of Texas Institute for Geophysics contribution number 468.     

国产OBS的时间矫正方法——以2019年台湾海峡地壳结构海陆探测实验为例

时间服务系统对利用走时层析成像方法进行地下介质速度结构反演至关重要。海底地震仪(ocean bottom seismometer, OBS)工作期间由于没有GPS时间接入, 其时间误差(包括守时误差和授时误差)主要来源于内部石英晶振的准确程度, 受到外部环境变化以及开关机等因素影响。长期实践发现, 部分国产OBS在记录气枪信号以及天然地震信号时存在较大的时间偏差。本文对2019年福建及台湾海峡地壳结构海陆探测实验所获得的53台次国产OBS记录进行了时间服务系统矫正。其中, 针对OBS授时误差, 利用出海前不中断采集的一致性试验和运输船运输过程中产生的晃动互相关进行时间矫正; 针对守时误差, 采用计算实际采样频率与理论采样频率偏差进行矫正; 通过对比矫正前后OBS记录到的天然地震信号, 进行秒级别的检测。结果表明, 经过以上步骤矫正的OBS数据, 其时间记录的准确性得到了显著提高, 从而降低了震相识别、走时拾取的时间误差, 为标准化国产OBS数据采集作业流程提供了重要参考。     

Ocean bottom seismometer detection of earthquakes on the Gorda Rise near 41.5° north

The spatial and temporal distribution of microearthquake activity associated with the southern section of the Gorda Rise was monitored during a three week period using an array of four ocean bottom seismometers (OBS's) located at 41.5° N, 127.5° W on the rise axis. Sixty-two events were detected, of which roughly half occurred during a seismic swarm following an event of magnitude 4.9 that was well-recorded on land. A total of eighteen earthquakes could be located, seventeen of which originated from the valley floor, walls, and crestal mountains while the remaining event was located to the east of the rise in the Gorda plate. Well-constrained focal depths for four earthquakes ranged from roughly 4 to 12 km. Depths in excess of 10 km for two of the events provide strong evidence for the absence of a shallow steady-state magma chamber beneath the slow-spreading portion of the Gorda rise. Most of the events are believed to be the result of uplift of the crustal blocks comprising the valley walls. Four events located at an offset of the rise near 41.5° N may be related to strike-slip movement on a developing transform fault. Two of the events that occurred during the survey were large enough to be detected and located by land-based instruments. Comparison of these locations with the OBS locations indicates that well-recorded events having magnitudes of at least 4.5 are routinely mislocated 40–50 km to the east of their actual locations in this area while lower magnitude, poorly recorded earthquakes may have significantly larger location errors.     

Anomalously large seismic amplifications in the seafloor area off the Kii peninsula

Seismic wave amplifications were investigated using strong-motion data obtained from the ground’s surface (K-net) on the Kii peninsula (southwestern Japan) and from the network of twenty seismic stations on the seafloor (DONET) located off the peninsula near the Nankai trough. Observed seismograms show that seismic signals at DONET stations are significantly larger than those at K-net stations, independent of epicentral distances. In order to investigate the cause of such amplifications, seismic wavefields for local events were simulated using the finite-difference method, in which a realistic 3D velocity structure in and around the peninsula was incorporated. Our simulation results demonstrate that seismic waves are significantly amplified at DONET stations in relation to the presence of underlying low-velocity sediment layers with a total thickness of up to 10 km. Our simulations also show considerable variations in the degree of amplification among DONET stations, which is attributed to differences in the thickness of the sediment layers. The degree of amplification is relatively low at stations above thin sediment layers near the trough axis, but seismic signals are much more amplified at stations closer to the Kii peninsula, where sediment layers are thicker than those at the trough axis. Simulation results are consistent with observations. This study, based on seafloor observations and simulations, indicates that because seismic signals are amplified due to the ocean-specific structures, the magnitude of earthquakes would be overestimated if procedures applied to data observed at land stations are used without corrections.     

台湾海峡强震活动及其动力环境

本文根据１９９４年９月１６日台湾海峡７．３级强震序列的空间活动图象变化特征，结合历史地震资料，比较分析了台湾海峡区域史今地震活动的总体演变过程，表明该区强震发生与周边地区的台湾，日本及菲律宾大地地震，在时间进程，强度变化和地域分布关系上存在关地球物理场变化的影响。     

Relationship of instruments and policy in the Hawaii warning system

The tsunami which follows a strong local earthquake occurs within a few minutes of the origin time. This requires that any warning system for such local tsunamis be decentralized, as in the system used by the Japan Meterological Agency. Thus, decisions are being made by local officials rather than in a central office. This usually means that the level of training, the education, and the preparation is not as intense in such a local office as in a central office. Therefore, the decision making should be assisted by instrumentation that discriminates tsunamigenic earthquakes from non‐tsunamigenic earthquakes. This is not yet possible in real time; however, since only very large earthquakes (more than 6.5) generate significant tsunamis, an instrument to inform the local official of the tsunami prospects can be implemented.

An instrument for assisting the local decision maker has been developed. This instrument consists of an analog computer (an inverted pendulum having a period of 0.75 sec and damping about 0.3) and a digital computer (a hardwired signal‐recognition circuit), providing output to a display of status or alarm. The level of displacement (or velocity or acceleration) , the number of times that the threshold must be exceeded, and the time window within which the excedance must occur are all adjustable. Initial settings require about 0.06 g three times, not more than ten seconds apart. Battery back‐up and test circuitry are provided.

Since the operating instructions require that motion be felt before the instrument alarm be considered valid, great weight is given to assuring that the instrument will operate when required. The instru‐ment is always “ON,”; with the earthquake turning it “OFF”; thus, it is continually self‐testing.

Twelve units have been installed in police stations, fire stations, or similar locations around the State of Hawaii, which funded the system developed at Indiana University. This tsunami seismic trigger should be considered whenever the primary objective is to trigger an alarm rather than to record data; the emphasis in design and development has been on reminding the local official when there may be a tsunami hazard and not on recording research data.     

Correlation between the intensity of tsunami in the Black Sea and the magnitude of underwater earthquakes

We perform the analysis of tsunami waves in the shelf zone of the Crimean peninsula generated by underwater earthquakes whose epicentres are located near the lower boundary of the continental slope. For this purpose, we use a one-dimensional nonlinear dissipative numerical model of long waves. The investigated distributions of the depth of the basin correspond to four points of the south coast of the Crimean peninsula. We use the empirical dependences of parameters of the sources of tsunami waves on the magnitude of the earthquake obtained earlier for the Pacific Ocean. It is shown that the height, vertical climb, and duration of tsunami waves increase with the magnitude of the earthquake. For equal magnitudes of the earthquake, the highest tsunamis on the south coast of the Crimea are observed between Alushta and Yalta. We also deduced a generalized regression dependence of the height of tsunami waves near the coast on the magnitude of the earthquake. Translated by Peter V. Malyshev and Dmitry V. Malyshev     

The estimation of ocean bottom seismographs' coupling characteristics by means of microshock recordings

A transient technique was used for estimating the bottom-instrument response function in experiments with bottom seismographs (OBS) in deep ocean. The sharp mechanical impacts on bottom seismograph usually induced by bottom displacements under the instrument (microshocks) are suggested as rough analogues of the impulses for the bottom-instrument system transient calibration. It has been found that microshocks has usually sufficiently small duration to be used for coupling characteristics estimates. Test measurements have shown that in many cases this method makes it possible to distinguish spectral features characterising earthquakes and seismic noise wavetrains from those which are caused by coupling resonances of the OBS-sediment inter-face.     

Variations in tidal response of the medium in the spatiotemporal neighborhood of catastrophic earthquake in Japan

The predictability of catastrophic earthquakes according to data on time variations in tidal response is discussed. We present the results of a harmonic analysis of (1) a one-year record of tide-gauge observations with the use of two long-term vertical seismometers at the MAJO station on Honshu Island near the epicenter of the M _w = 9 catastrophic earthquake of March 11, 2011, and (2) a 0.5-year record of horizontal pendulum observations at the ERM station near Hokkaido Island from March 13, 2010, to March 12, 2011, and from April 26, 2010, to August 3, 2010, respectively. The relative changes in the tidal amplitudes are found to be about 3 and 4%. The identification of distinct prognostic attributes awaits much more representative spatiotemporal statistics. These results can be adequately and reliably interpreted by analyzing “background” variations in the tidal amplitudes at these stations and, in particular, their possible seasonal variations. This requires at least several years of observations.     

Analysis of thermohaline and current distribution characteristics of Zhejiang and Fujian waters in summer

This paper analyzed the distribution of thermohaline and circulation characteristics of Zhejiang and Fujian waters, based on the cross-sectional thermohaline data and on current data (up to 30 d duration) at fixed-point moorings, collected in the summer of 2006. We also performed low-pass filtering and spectrum analysis on the mooring submersible buoy data. Based on that analysis, we discussed the characteristics of low frequency currents and time-variations in these waters. The main conclusions are as follows. (1) There is a low salinity pinnate area near the Hangzhou Bay in summer, and outside the low salinity area, an obvious salinity front is present from surface to bottom near 123°E. There is also a temperature front below the surface at a corresponding position. (2) Bottom water of the Taiwan Warm Current comes from the subsurface of Kuroshio. (3) The direction of low frequency current at fixed anchor stations is N-NE or S, which mainly depends on the interaction of control currents in this waters. (4) Significant spectral peaks at all mooring stations are typically semidiurnal and diurnal tides. Semidiurnal tidal waves are the main ones in these waters, and have more energy closer to the shore. (5) Significant energy spectral peaks of middle period (3 to 8 d) of currents are responses to weather frequency. (6) Significant energy spectral peaks of long periods at the surface or bottom are probably responses to seasonal wind or bottom friction, while, the long period peaks of other depths can reflect cyclical changes of interactions between currents. We conclude that the pulsation period of the Taiwan Warm Current in these waters is 10-17 d.     

Possible triggering relationship of six M_w>6 earthquakes in 2018–2019 at Philippine archipelago

Philippine archipelago(PA) has strong background seismicity, but there is no systematic study of earthquake triggering in this region. There are six earthquakes(M_w 6) occurred between 2018/12/29 and 2019/09/29 in PA,which provides an excellent opportunity to investigate the triggering relationship among these events. We calculate the static Coulomb stress changes of the first five events, and find that the local seismicity after the2018/12/29 M_w 7.0 earthquake is mostly associated with positive Coulomb stress changes, including the2019/05/31 M_w 6.1 event, suggesting a possible triggering relationship. However, we cannot rule out the dynamic triggering mechanism, due to increased microseismicity in both positive and negative stress change regions, and an incomplete local catalog, especially right after the first M_w 7.0 mainshock. The dynamic stresses from these M_w 6 events are large enough(from 5 kPa to 3 532 kPa) to trigger subsequent events, but a lack of seismicity and waveform evidence does not support delayed dynamic triggering among these events, even the shortest time interval is less than 24 hours. In the past 45 years, the released seismic energy shows certain peaks every 5–10 years. However, earthquakes with M_w 6.0 were relatively infrequent between 2004 and 2018 at PA. Hence, it is possible that several regions are relatively late in their earthquake cycles, which would enhance their susceptibility of being triggered by earthquakes at nearby and regional distances.     

The detectability of shallow slow earthquakes by the Dense Oceanfloor Network system for Earthquakes and Tsunamis (DONET) in Tonankai district,Japan

In order to understand the characteristics of shallow very low-frequency (VLF) events as revealed by recent ocean-floor observation studies, we perform a trial simulation of earthquake cycles in the Tonankai district by taking the characteristics of the 1944 Tonankai earthquake and assuming that slow earthquakes occur on numerous small asperities. Our simulation results show that the increase of moment release rate of shallower VLF events in the pre-seismic stage of a megathrust earthquake is higher than that of deeper VLF events. This increase may make leveling change due to VLF swarms detectable at Dense Oceanfloor Network system for Earthquakes and Tsunamis (DONET). We also introduce the time series of hydraulic pressure data at DONET, comparing with the leveling change expected from our numerical simulation. Since leveling change due to shallower VLF swarms is so local as to be incoherent, removal of the moving-averaged data from the data stacked by four nearby observation points in the same node may be useful to detect the short-term local leveling change.