NUTRIENT LOAD ESTIMATION METHODS FOR RIVERS

Pollutant load estimation is essential for watershed management and water pollution control. For most watersheds, only sparse water quality measurements (e.g. monthly data) are available. The influence of input data on the accuracy of non-point source pollution load estimation is studied using the water quality and stream flow data from a small watershed in Hong Kong. Comparison and analysis of the results using 8 different methods show that the accuracy of stream-flow runoff is the single most important factor for the calculation of pollutant load. Based on the results, the stream flow correction coefficient is advanced to provide a more reliable load estimation. The improved method of pollutant load estimation can be easily applied in practice since the stream-flow runoff can be measured by hydrological station or estimated with various hydrological methods.     

Earthquake magnitude estimation using the τc and Pd method for earthquake early warning systems

Earthquake early warning （EEW） systems are one of the most effective ways to reduce earthquake disaster. Earthquake magnitude estimation is one of the most important and also the most difficult parts of the entire EEW system. In this paper, based on 142 earthquake events and 253 seismic records that were recorded by the KiK-net in Japan, and aftershocks of the large Wenchuan earthquake in Sichuan, we obtained earthquake magnitude estimation relationships using the τe and Pa methods. The standard variances of magnitude calculation of these two formulas are ±0.65 and ±0.56, respectively. The Pd value can also be used to estimate the peak ground motion of velocity, then warning information can be released to the public rapidly, according to the estimation results. In order to insure the stability and reliability of magnitude estimation results, we propose a compatibility test according to the natures of these two parameters. The reliability of the early warning information is significantly improved though this test.     

Study on the method of optimal location on explosion events at near-source site

IntroductionEarthquake location is one of the inverse problems in seismology. The accurate earthquake source location is the basis of studies of earthquake images, the relation analysis between seismicity and tectonics, and the important parameters of earthquakes. It is also the basic data used to study the space distribution of aftershocks after a strong earthquake. The earthquake location accuracy mainly depends on reliability of the velocity models for the earth(s crust, the reasonable desi…     

汶川地震余震活动固体潮调制的统计分析

A statistical analysis on the Wenchuan aftershock activity triggered by tidal forces is systematically studied based on Schusters test, including earthquakes triggered by tidal force, tidal stress and tidal coulomb failure stress. The results show that a group of strong aftershocks which occurred at the end of July to early August in 2008 at the north of Wenchuan were obviously triggered by earth tide, the same conclusion is drawn by Schusters smooth test of the tidal force, tidal stress and tidal coulomb failure stress. In addition, the Wenchuan aftershock activity is obviously triggered by fortnight tide. In the north, the aftershocks happened more frequently in the first and last quarters of the moon, and in the south, the aftershocks happened more frequently in the first and last quarters of the moon and during the full moon.     

Effect of Aftershocks on Assessment of Seismic Hazard for Hebei Seismic Zone,China

In this paper,the earthquake hazard parameters(λm,b and Mmax) from the maximum likelihood method for the raw catalogue and declustered catalogues have been used to discuss the effect of the aftershocks on the earthquake hazard estimation.The declusterd catalogues have been compoiled from the raw catalogues by deleting the aftershocks in different time interval after main shock according to the criteria for the aftershock activity period.As an example,Hebei seismic zone is taken to show the effect of the aftershock on the earthquake hazard assessment because three strong earthquakes with the aftershocks occurred from 1966,to 1976,The results have been shown that the effect of the aftershocks on the parameters λm,b is significant.The difference between the clustered and declustered catalogues has reached in seismic activity rate and recurrence period,at most,45% and 90%,respectively,But the diffeence in Mmax is smaller,Based on this,the suggestion that the aftershocks should be omitted in the estimation of the earthquake hazard could be made,but how long the aftershocks activity duration is still left to the future research.     

Relocation of aftershocks of the Wenchuan M_S8.0 earthquake and its implication to seismotectonics

From 14:28 (GMT+8) on May 12th,2008,the origin time of MS8.0 Wenchuan earthquake,to December 31th,2008,more than 10 000 aftershocks (M≥2.0) had been recorded by the seismic networks in Sichuan and surrounding areas.Using double difference algorithm,the main shock and more than 7 000 aftershocks were relocated.The aftershocks distribute about 350 km long.The depths of aftershocks are mainly between 10 km and 20 km.The average depth of aftershocks is about 13 km after relocation.In the southwest,the distribut...     

Research on the Result of Calibration for the Location of the Seismic Array in Shanghai

The paper proposes the calibration theory and methods for the Shanghai seismic array and analyzes the calibration results. As a result, the calibration results for the seismic array based on 2 typical earthquakes have been drawn; the difference of calibration results between Hokkaido and Honshu region in Japan is investigated. And calibration results of different directions, different epicenter distances and different magnitudes are probed into. The result shows that the location of earthquakes on the Shanghai seismic array is greatly improved.     

Accurate Location of Synthetic Acoustic Emissions and Location Sensitivity to Relocation Methods, Velocity Perturbations, and Seismic Anisotropy

Acoustic emission (AE) monitoring is a non-invasive method of monitoring fracturing both in situ, and in experimental rock deformation studies. Until recently, the major impediment for imaging brittle failure within a rock mass is the accuracy at which the hypocenters may be located. However, recent advances in the location of regional scale earthquakes have successfully reduced hypocentral uncertainties by an order of magnitude. The least-squares Geiger, master event relocation, and double difference methods have been considered in a series of synthetic experiments which investigate their ability to resolve AE hypocentral locations. The effect of AE hypocenter location accuracy due to seismic velocity perturbations, uncertainty in the first arrival pick, array geometry and the inversion of a seismically anisotropic structure with an isotropic velocity model were tested. Hypocenters determined using the Geiger procedure for a homogeneous, isotropic sample with a known velocity model gave a RMS error for the hypocenter locations of 2.6 mm; in contrast the double difference method is capable of reducing the location error of these hypocenters by an order of magnitude. We test uncertainties in velocity model of up to ±10% and show that the double difference method can attain the same RMS error as using the standard Geiger procedure with a known velocity model. The double difference method is also capable of precise locations even in a 40% anisotropic velocity structure using an isotropic model for location and attains a RMS mislocation error of 2.6 mm that is comparable to a RMS mislocation error produced with an isotropic known velocity model using the Geiger approach. We test the effect of sensor geometry on location accuracy and find that, even when sensors are missing, the double difference method is capable of a 1.43 mm total RMS mislocation compared to 4.58 mm for the Geiger method. The accuracy of automatic picking algorithms used for AE studies is ±0.5 μs (1 time sample when the sampling rate is 0.2 μs). We investigate how AE locations are effected by the accuracy of first arrival picking by randomly delaying the actual first arrival by up to 5 time samples. We find that even when noise levels are set to 5 time samples the double difference method successfully relocates the synthetic AE.     

Relocation of the M8.0 Wenchuan earthquake and its aftershock sequence

We relocated M8.0 Wenchuan earthquake and 2706 aftershocks with M≥2.0 using double-difference algorithm and obtained relocations of 2553 events. To reduce the influence of lateral variation in crustal and upper mantle velocity structure, we used different velocity models for the east and west side of Longmenshan fault zone. In the relocation process, we added seismic data from portable seismic sta-tions close to the shocks to constrain focal depths. The precisions in E-W, N-S, and U-D directions after relocation are 0.6, 0.7, and 2.5 km respectively. The relocation results show that the aftershock epi-centers of Wenchuan earthquake were distributed in NE-SW direction, with a total length of about 330 km. The aftershocks were concentrated on the west side of the central fault of Longmenshan fault zone, excluding those on the north of Qingchuan, which obviously deviated from the surface fault and passed through Pingwu-Qingchuan fault in the north. The dominant focal depths of the aftershocks are between 5 and 20 km, the average depth is 13.3 km, and the depth of the relocated main shock is 16.0 km. The depth profile reveals that focal depth distribution in some of the areas is characterized by high-angle westward dipping. The rupture mode of the main shock features reverse faulting in the south, with a large strike-slip component in the north.     

Strong motion observations and recordings from the great Wenchuan Earthquake

The National Strong Motion Observation Network System （NSMONS） of China is briefly introduced in this paper. The NSMONS consists of permanent free-field stations, special observation arrays, mobile observatories and a network management system. During the Wenchuan Earthquake, over 1,400 components of acceleration records were obtained from 460 permanent free-field stations and three arrays for topographical effect and structural response observation in the network system from the main shock, and over 20,000 components of acceleration records from strong aftershocks occurred before August 1, 2008 were also obtained by permanent free-field stations of the NSMONS and 59 mobile instruments quickly deployed after the main shock. The strong motion recordings from the main shock and strong aftershocks are summarized in this paper. In the ground motion recordings, there are over 560 components with peak ground acceleration （PGA） over 10 Gal, the largest being 957.7 Gal. The largest PGA recorded during the aftershock exceeds 300 Gal.     

上海地震台阵的标定方法

地震台阵标定可以校正视慢度和后方位角，从而进一步校正震源位置，有助于提高地震定位精度及地震类型鉴别研究等。此外，台阵标定也可校正由于台阵接收系统而导致视慢度和后方位角的偏差。     

Locating Small Aftershocks Using a Small-Aperture Temporary Seismic Array

In this paper, we developed a specialized method to locate small aftershocks using a small-aperture temporary seismic array. The array location technique uses the first P arrival times to determine the horizontal slowness vector of the incoming P wave, then combines it with S–P times to determine the event location. In order to reduce the influence of lateral velocity variation on the location determinations, we generated slowness corrections using events well-located by the permanent broadband network as calibration events, then we applied the corrections to the estimated slownesses. Applications of slowness corrections significantly improved event locations. This method can be a useful tool to locate events recorded by temporary fault-zone arrays in the near field but unlocated by the regional permanent seismic network. As a test, we first applied this method to 64 well-located aftershocks of the 1992 Landers, California, earthquake, recorded by both the Caltech/USGS Southern California Seismic Network and a small-aperture, temporary seismic array. The average horizontal and vertical separations between our locations and the well-determined catalogue locations are 1.35 and 1.75 km, respectively. We then applied this method to 132 unlocated aftershocks recorded only by the temporary seismic array. The locations show a clear tendency to follow the surface traces of the mainshock rupture.     

Slowness–azimuth corrections of teleseismic events for IMS primary arrays in China

The seismic arrays at Hailar (HILR) and at Lanzhou (LZDM) in China are both primary stations of the International Monitoring System for verifying compliance with the Comprehensive Nuclear Test Ban Treaty. These two stations became operational in 2002 and have since then provided continuous data. In this study, the so-called slowness–azimuth station corrections (SASC) were derived and used to improve the location accuracy of the two arrays. The SASC are found by comparing the back-azimuth and slownesses obtained from array processing to the theoretical values calculated from the reported event locations and the corresponding seismic velocity model. Events reported by the National Earthquake International Center in the time period 2002 to 2006 were used as reference events, and the IASP91 was used as the theoretical velocity model. Small correction vectors with random orientation were found for HILR. Larger correction vectors with systematic vector biases were found for LZDM. The LZDM correction vectors seem to point to the same direction in a large part of the slowness space and may be attributed to local structure. After introducing the SASC for HILR, the standard deviations of back-azimuth and slowness residuals drop from 7.1° to 4.6° and from 1.0 to 0.6 s/°, respectively. For LZDM, these values drop from 22.3° to 10.2° and from 2.9 to 1.1 s/°, respectively. The variations of back-azimuth and slowness residuals were reduced by 32% and 30.2%, respectively, for HILR after SASC and the reductions were 21% and 40.2% for LZDM. The improvements were 77% in back-azimuth and 67% in slowness location for HILR and were 79% and 81% for LZDM after SASC. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

格尔木地震台阵勘址数据分析与台阵布局设计

为了增强西部地区的地震监测能力, 中国拟在格尔木地区建立一个小孔径地震台阵。 本文对台阵勘址数据进行了噪声与信号的相关性分析, 得到地震监测的最佳台间距。 结果表明, 对于近震和区域震的监测, 该台阵子台间距最好小于500 m, 对于远震的监测, 台站组合间距应为1500～2000 m。 最后将勘址布设的台阵作为初选台阵进行了台阵响应计算, 计算显示, 台阵响应的主瓣在NW-ES方向较窄, 表明对来自该方向事件的慢度分辨率较高; 由于呈“L”型分布, 该初选台阵确定某些方向地震的方位角较好, 但检测其他方向事件的方位角精度不高, 这可以通过台阵校正进行改善; 台阵响应中出现的多处侧瓣是由于子台间距较大造成的。     

上海地震台阵的地震定位方法

上海地震台阵数据处理软件系统的地震定位方法,采用台阵的聚束方法得到地震方位角和视慢度,根据统计得到的视慢度——震中距表推算震中距。并结合了地震台网的定位方法,由单台记录的各类主要震相从JB走时表得到震中距,然后进行地震定位。该定位方法可对近震、远震进行定位处理,并由深震相得到震源深度。     

Slowness Corrections — One Way to Improve IDC Products

—?The first step to identify and locate a seismic event is the association of observed onsets with common seismic sources. This is especially important in the context of monitoring the Comprehensive Nuclear-Test-Ban Treaty (CTBT) at the International Data Center (IDC) being developed in Vienna, Austria. Well-defined slowness measurements are very useful for associating seismic phases to presumed seismic events.¶Shortly after installation of the first seismic arrays, systematic discrepancies between measured and theoretically predicted slowness values were observed, and therefore slowness measurements of seismic stations should be calibrated. The observed slownesses measured with small aperture arrays, some of which will be included in the International Monitoring System (IMS) now being implemented for verifying compliance with the CTBT, show large scatter and deviations from theoretically expected values. However, in this study a method is presented, by which mean slowness corrections can be derived, which show relatively stable patterns specific to each array.¶The correction of measured slowness values of these arrays clearly improved the single array location capabilities. Applying slowness corrections with seismic phases observed by ARCES, FINES, GERES, and NORES, and associated to seismic events in the bulletins of the prototype International Data Center (pIDC) in Arlington, VA, also clearly demonstrates the advantages of these corrections. For arrays with large slowness deviations that are due to the influence of a dipping layer, the corrections were modeled with a sine function depending on the measured azimuth. In addition, the measured values can be weighted with the corresponding uncertainties known from the process of deriving the mean corrections.     

Fault zone structures of northern and southern portions of the main central fault generated by the 2008 Wenchuan earthquake using fault zone trapped waves

The rupture process of the May 12,2008 MS8.0 Wenchuan earthquake was very complex. To study the rupture zones generated by this earthquake,four dense temporary seismic arrays across the two surface breaking traces of the main-shock were deployed in July and recorded a great amount of aftershocks. This paper focuses on the data interpretation of two arrays across the central main fault,the northern array line 1 and southern array line 3. The fault zone trapped waves recorded by the two arrays were used to st...     

地震台阵上信号方位角和慢度的时、频域估计方法比较

地震台阵常用频率-波数分析法来估计地震信号的后方位角和慢度。尽管有关的算法是众所周知的,但不同的实现在某些情况下可能导致不同的结果。如应用在宽频带的澳大利亚WRA台阵的记录时,标准的f-k分析方法往往给出不正确的结果。本文发现错误是由FFT的频谱泄漏效应引起的,如在进行FFT之前对原始数据进行高通滤波则可以有效地避免这样的错误。进一步对在时域中直接估算信号方位角和慢度的一种算法进行了分析,比较发现时域方法的计算速度和精度都不低于频域方法,且在某些低信噪比的情况下,前者可以给出更可靠的结果。     

Slowness and azimuth determination for Bucovina array (BURAR) applying multiple signal techniques

The BURAR seismic array, located in Northern Romania (Bucovina region), is designed to monitor events located in an area poorly covered by other existing seismic stations. In order to use the BURAR array for single-station locations, it is crucial to calibrate the azimuth and slowness parameters, which are currently used in array techniques to locate earthquakes, blasts or nuclear explosions. The goal of this study is to apply “f–k” and plane wave fit techniques in order to constrain the slowness and azimuth parameters at BURAR for teleseismic, regional and local events. The analysis was carried out using P and S waves recorded for events occurred between 2004 and 2008 within a radius of 50° around BURAR. The azimuth values obtained applying both methods strongly deviated from the theoretical values for regions like Central Turkey, Bulgaria, Dodecanese Islands and other parts of Greece, while the ray parameter deviations with respect to a 1-D IASP91 reference model are less significant. For the local events, the anomalies are smaller, except the particular case of Vrancea intermediate-depth earthquakes for which strong azimuth deviations (33.5°), both positive and negative, are observed. We investigate how these systematic deviations in azimuth are explained by the structure lateral heterogeneities which characterize the study region.     

上海地震台阵的F-K分析

频率－波数（F－K）谱分析是地震台阵数据处理基本方法之一。采用该分析方法，可以从地脉动背景中提取有用的地震信号，提高地震的定位精度。本文列举了F－K功率谱的不同计算方法，并将该方法应用于上海地震台阵所记录的地震事件。对地震和地脉动的F－K分析结果得出，地震的方位角和慢度值分别对记录长度的分布呈线性特征，而地脉动的方位角和慢度值则呈不规则分布。