GBO algorithm for seismic source parameters inversion

The use of geodetic observation data for seismic fault parameters inversion is the research hotspot of geodetic inversion, and it is also the focus of studying the mechanism of earthquake occurrence. Seismic fault parameters inversion has nonlinear characteristics, and the gradient-based optimizer(GBO) has the characteristics of fast convergence speed and falling into local optimum hardly. This paper applies GBO algorithm to simulated earthquakes and real LuShan earthquakes in the nonlinear inversion of the Okada model to obtain the source parameters. The simulated earthquake experiment results show that the algorithm is stable, and the seismic source parameters obtained by GBO are slightly closer to the true value than the multi peak particle swarm optimization(MPSO). In the 2013 LuShan earthquake experiment, the root mean square error between the deformation after forwarding of fault parameters obtained by the introduced GBO algorithm and the surface observation deformation was 3.703 mm, slightly better than 3.708 mm calculated by the MPSO. Moreover, the inversion result of GBO algorithm is better than MPSO algorithm in stability. The above results show that the introduced GBO algorithm has a certain practical application value in seismic fault source parameters inversion.

Effect of source parameters on forward-directivity velocity pulse for vertical strike slip fault in half space

It has been found that the large velocity pulse is one of the most important characteristics of near-fault strong ground motions. Some statistical relationships between pulse period and the moment magnitude for near-fault strong ground motions have been established by Somerville （1998）; Alavi and Krawinkler （2000）; and Mavroeidis and Papageorgiou （2003）, where no variety of rupture velocity, fault depth, and fault distance, etc. were considered. Since near-fault ground motions are significantly influenced by the rupture process and source parameters, the effects of some source parameters on the amplitude and the period ofa forward-directivity velocity pulse in a half space are analyzed by the finite difference method combined with the kinematic source model in this paper. The study shows that the rupture velocity, fault depth, position of the initial rupture point and distribution of asperities are the most important parameters to the velocity pulse. Generally, the pulse period decreases and the pulse amplitude increases as the rupture velocity increases for shallow crustal earthquakes. In a definite region besides the fault trace, the pulse period increases as the fault depth increases. For a uniform strike slip fault, rupture initiating from one end of a fault and propagating to the other always generates a higher pulse amplitude and longer pulse period than in other cases.

Changes in source parameters of foreshocks and aftershocks of the 2001 M_S=6.0 Yajiang, Sichuan, earthquake

In this paper changes in focal mechanisms, parameters of wave spectra, and stress drops for the MS=5.0 foreshock and MS=6.0 mainshock in February 2001 in Yajiang County, Sichuan, and seismicity in epicentral region are studied. Comparison of focal mechanisms for the Yajiang earthquakes with distribution patterns of aftershocks, the nodal plane I, striking in the direction of NEN, of the Yajiang M=5.0 event is chosen as the faulting plane; the nodal plane II, striking in the direction of WNW, of the M=6.0 event as the faulting plane. The strikes of the two faulting planes are nearly perpendicular to each other. The level of stress drops in the epicentral region before the occurrence of the M=6.0 earthquake increases, which is consistent with increase of seismicity in the epicentral region. The rate decay of the Yajiang earthquake sequence, changes in wave spectra for foreshocks and aftershocks, and focal mechanisms are complex.

Inversion of source parameters for moderate and small earthquakes in Beijing region

According to the geological structural features, Beijing and the adjacent areas can be divided into two regions of plain in the east and mountain in the west. Among the stations covered by the telemetered digital seismic station network of Earthquake Administration of Beijing Municipality, the stations in the plain area are all borehole ones and the stations in the western mountainous region are all located on the surface bedrock. In the paper, 511 wave- form data recorded by the network from Oct. 2001 to Oct. 2004 are used in the researches for the entire Beijing region, the western mountainous region and the eastern plain area, respectively. The Q values are calculated for each area by Atkinson′s method and compared with the existed data. The reliability of the Q values and the reasons for the difference in the Q values are also discussed. Then, the source parameters and site response are inverted by the Moya′s method, in which two models are used. The first model uses the Q values, earthquakes and stations in the sub-areas and the second model uses the Q values, earthquakes and stations in the entire Beijing region. The results indicate that the source parameters and site responses obtained by two models are basically consistent with each other. It also indicates that the source parameters obtained by these methods are not affected by the size of station network.

Study on inelastic attenuation coefficient, site response and source parameters in Shanxi region

Based on 310 horizontal-component digital seismograms recorded at 14 seismic stations in Shanxi Digital Seis-mograph Network, the inelastic attenuation coefficient in Shanxi region is studied. By the methods of Atkinson and Moya, the site response of each station and several source parameters are obtained and the inversion results from both methods are compared and analyzed. The frequency-dependent inelastic attenuation coefficient Q is estimated as Q( f )=323.2 f 0.506. The site responses of 14 seismic stations do not show significant amplification, which is consistent with their basement on rock. We also found the dependence of corner frequency on seismic moment, seismic moment on stress drop, source radius on stress drop.

Earthquake source parameters of the 2009 M_W7.8 Fiordland (New Zealand) earthquake from L-band InSAR observations

The 2009 M W 7.8 Fiordland （New Zealand） earthquake is the largest to have occurred in New Zealand since the 1931 M W 7.8 Hawke’s Bay earthquake, 1 000 km to the northwest. In this paper two tracks of ALOS PALSAR interferograms （one ascending and one descending） are used to determine fault geometry and slip distribution of this large earthquake. Modeling the event as dislocation in an elastic half-space suggests that the earthquake resulted from slip on a SSW-NNE orientated thrust fault that is associated with the subduction between the Pacific and Australian Plates, with oblique displacement of up to 6.3 m. This finding is consistent with the preliminary studies undertaken by the USGS using seismic data.

A new approach for separating mixed model parameters:application to simultaneous inversion of earthquake source parameters

A method for simultaneous determination of mixed model parameters,which have different physical dimensions or different responses to data,is presented.Mixed parameter estimation from observed data within a single model space shows instabilities and trade-offs of the solutions. We separate the model space into N-subspaces based on their physical properties or computational convenience and solve the N-subspaces systems by damped least-squares and singular-value decomposition. Since the condition number of each subsystem is smaller than that of the single global system,the approach can greatly increase the stability of the inversion. We also introduce different damping factors into the subsystems to reduce the tradeoffs between the different parameters. The damping factors depend on the conditioning of the subsystems and may be adequately chosen in a range from 0.1 % to 10 % of the largest singular value. We illustrate the method with an example of simultaneous determination of source history,source geometry,and hypocentral location from regional seismograms,although it is applicable to any geophysical inversion.

Applying Source Parameter Imaging (SPI) to Aeromagnetic Data to Estimate Depth to Magnetic Sources in the Mamfe Sedimentary Basin

Aeromagnetic data over the Mamfe Basin have been processed. A regional magnetic gridded dataset was obtained from the Total Magnetic Intensity (TMI) data grid using a 3 × 3 convolution (Hanning) filter to remove regional trends. Major similarities in magnetic field orientation and intensities were observed at identical locations on both the regional and TMI data grids. From the regional and TMI gridded datasets, the residual dataset was generated which represents the very shallow geological features of the basin. Processing this residual data grid using the Source Parameter Imaging (SPI) for magnetic depth suggests that the estimated depths to magnetic sources in the basin range from about 271 m to 3552 m. The highest depths are located in two main locations somewhere around the central portion of the study area which correspond to the area with positive magnetic susceptibilities, as well as the areas extending outwards across the eastern boundary of the study area. Shallow magnetic depths are prominent towards the NW portion of the basin and also correspond to areas of negative magnetic susceptibilities. The basin generally exhibits a variation in depth of magnetic sources with high, average and shallow depths. The presence of intrusive igneous rocks was also observed in this basin. This characteristic is a pointer to the existence of geologic resources of interest for exploration in the basin.

Towards combining multiple geophysical datasets to determine earthquake source parameters in China

1.Difficulties of conventional seismic studies on earthquake source parameters Earthquake source parameters,including magnitude,location,focal mechanism,rupture process are key factors for understanding seismogenic environment,mitigating seismic hazards,estimating earthquake triggering,and tectonic analysis.Traditionally,source parameters are determined by seismological methods.For example,Fang L H et al.(2014)relocated the 2012 Ms6.6 Xinjiang Xinyuan earthquake sequence using local seismograms based on the double difference method,and obtained the distribution of

Source parameters of the Gonghe,Qinghai Province,China,earthquake from inversion of digital broadband waveform data

An earthquake of M S=6.9 occurred in Gonghe County, Qinghai Province, China on April 26, 1990.This earthquake was followed by three larger aftershocks of M S=5.5 on May 7, 1990, M S=6.0 on Jan.3, 1994, and M S=5.7 on Feb.16, 1994, consecutively. The moment tensors of these earthquakes as function of time were obtained by the technique of moment tensor inversion in frequency domain . The results inverted indicate that these earthquakes had a very similar focal mechanism of predominantly reverse faulting on a plane striking NWW, dipping to SSW.The scalar seismic moments of these earthquakes are M 0=9.4×10 18 Nm for the M S=6.9 event, 8.0×10 16 Nm for the M S=5.5 event, 4.9×10 17 Nm for the M S =6.0 event and 2.9×10 17 Nm for the M S=5.7 event, respectively. The results inverted also show that the source processes of these events were significantly different. The main shock had a very complex process, consisting of two distinct sub events with comparable sizes. The first sub event occurred in the first 12s, having a seismic moment of 4.7×10 18 Nm, and the second one continued from 31s to 41s, having a seismic moment of 2.5×10 18 Nm. In addition, a much smaller sub event, having a seismic moment of about 2.1×10 18 Nm, may exist in the interval of 12 s and 31 s, In contrast, the source processes of the three aftershocks are quite simple. The source time function of each of aftershocks is a single impulse, suggestting that each of aftershocks consists of a mainly uninterrupted rupture. The rise times and total rupture durations are 4 s and 11 s for the M S=5.5 event, 6 s and 16 s for the M S= 6.0 event and 6 s and 13 s for the M S=5.7 event, respectively.

ThesourceparametersandsourcemechanismfortheJanuary5,1998JingyangearthquakewithMS=4.8

ＡｎｅａｒｔｈｑｕａｋｅｗｉｔｈＭＳ＝４．８ｏｃｕｒｅｄｎｅａｒＹｏｎｇｌｅＴｏｗｎｉｎＪｉｎｇｙａｎｇＣｏｕｎｔｙ，ＳｈａａｎｘｉＰｒｏｖｉｎｃｅｏｎＪａｎｕａｒｙ５，１９９８．Ｉｔｓｅｐｉｃｅｎｔｅｒｉｓ３４．５°Ｎ，１０８．９°Ｅ，ｏｒｉｇｉｎ...

Regional Finite-Fault Source Model for Development of Ground Motion Attenuation Relationship in Sichuan, China

The attenuation relationship of ground motion based on seismology has always been a front subject of engineering earthquake.Among them,the regional finite-fault source model is very important.In view of this point,the general characteristics of regional seism-tectonics,including the dip and depth of the fault plane,are emphasized.According to the statistics of regional seism-tectonics and focal mechanisms in Sichuan,China,and the sensitivity of estimated peak ground acceleration(PGA)attenuation is analyzed,and the dip angle is taken as an average of 70°.Based the statistics of the upper crustal structure and the focal depth of regional earthquakes,the bottom boundary of the sedimentary cover can be used as the upper limit for estimating the depth of upper-edge.The analysis shows that this value is sensitive to PGA.Based on the analysis of geometric relations,the corresponding calculation formula is used,and a set of concepts and steps for building the regional finite-fault source model is proposed.The estimation of source parameters takes into account the uncertainty,the geometric relationship among parameters and the total energy conservation.Meanwhile,a set of reasonable models is developed,which lay a foundation for the further study of regional ground motion attenuation based on seismology.

The K Method for Estimating Earthquake Activity Parameters and Effect of the Boundary Uncertainty of the Source Region:Discussion on the Seismic Zoning Method

Two aspects of a new method,which can be used for seismic zoning,are introduced in this paper.On the one hand,the approach to estimate b value and annual activity rate proposed by Kijko and Sellevoll needs to use the earthquake catalogue.The existing earthquake catalogue contains both historical and recent instrumental data sets and it is inadequate to use only one part.Combining the large number of historical events with recent complete records and taking the magnitude uncertainty into account,Kijko’s method gives the maximum likelihood estimation of b value and annual activity rate,which might be more realistic.On the other hand,this method considers the source zone boundary uncertainty in seismic hazard analysis,which means the earthquake activity rate across a boundary of a source zone changes smoothly instead of abruptly and avoids too large a gradient in the calculated results.

Source Characteristics of the 2012 Ahar-Varzaghan Earthquake

The August 11, 2012 Mw 6.4 earthquakes in northwestern Iran occurred as a result of oblique strike-slip faulting in the shallow crust of the Eurasia plate, approximately 300 kilometer east of the plate boundary between the Eurasia and Arabia plates. The seismotectonics of this region is controlled by the collision of the Arabia and Eurasia plates;at the latitude of the earthquakes, the Arabia plate moves almost due north with respect to the Eurasia plate at a rate of approximately 26 millimeter per year. Over the past forty years, seven earthquakes of Mw 6 or greater have occurred within 300 kilometer of today’s events. The nearest was a Mw 6.1 earthquake in February of 1997, approximately 100 kilometer to the east, which caused 1100 fatalities. The studied area encourage the authors to determine the focal mechanism, source time function and sub events’ features which can guide us to reliable judges.

UBSS and blind parameters estimation algorithms for synchronous orthogonal FH signals

By using the sparsity of frequency hopping（FH） signals,an underdetermined blind source separation（UBSS） algorithm is presented. Firstly, the short time Fourier transform（STFT） is performed on the mixed signals. Then, the mixing matrix, hopping frequencies, hopping instants and the hooping rate can be estimated by the K-means clustering algorithm. With the estimated mixing matrix, the directions of arrival（DOA） of source signals can be obtained. Then, the FH signals are sorted and the FH pattern is obtained. Finally, the shortest path algorithm is adopted to recover the time domain signals. Simulation results show that the correlation coefficient between the estimated FH signal and the source signal is above 0.9 when the signal-to-noise ratio（SNR） is higher than 0 d B and hopping parameters of multiple FH signals in the synchronous orthogonal FH network can be accurately estimated and sorted under the underdetermined conditions.

A single-station coda solution for source,attenuation and site factors

ased on the empirical formulation formed from coda observations, the digital waveforms from 33 local events with magnitude ML ranging between 2.1 and 3.5, recorded at Changli station of Beijing Telemetered Seismographic Network from 1989 to 1991, are used to compute coda Q for the Changli region and the source factors of all earthquakes by single-station coda method. Furthermore, assuming a certain source model, we have obtained the station site frequency response and source spectra, as well as source parameters such as corner frequencies, seismic moments and stress drops and so on. Their variations with time are monitored before and after larger earthquakes. Because the coda method can more effectively reduce the influence of source radiation pattern and a particular propagation path than direct wave method, more data can be used and more accurate results can be obtained, which provided a possible approach to study the source properties and reveal the variation of source parameters before larger earthquakes.

Studies on seismic source

During the period of 1999~2002, the Chinese seismologists made a serious of developments in the study on seismic sources including observations, experiments and theory. In the field of observation, the methods of the accuracy location of earthquake sources, the inversion of seismic moment tensor and the mechanism of earthquake source are improved and developed. A lot of important earthquake events are studied by using these methods. The rupture processes of these events are inverted and investigated combined with the local stress fields and the tectonic moment by using the measurements of surface deformation. In the fields of experiments and theory, many developments are obtained in cause of seismic formation, condition of stress and tectonics, dynamics of earthquake rupture, rock fracture and nucleation of strong earthquakes.

Discharge Characteristics of Large-Area High-Power RF Ion Source for Positive and Negative Neutral Beam Injectors

A large-area high-power radio-frequency（RF） driven ion source was developed for positive and negative neutral beam injectors at the Korea Atomic Energy Research Institute（KAERI）. The RF ion source consists of a driver region, including a helical antenna and a discharge chamber, and an expansion region. RF power can be transferred at up to 10 kW with a fixed frequency of 2 MHz through an optimized RF matching system. An actively water-cooled Faraday shield is located inside the driver region of the ion source for the stable and steady-state operations of high-power RF discharge. Plasma ignition of the ion source is initiated by the injection of argongas without a starter-filament heating, and the argon-gas is then slowly exchanged by the injection of hydrogen-gas to produce pure hydrogen plasmas. The uniformities of the plasma parameter,such as a plasma density and an electron temperature, are measured at the lowest area of the driver region using two RF-compensated electrostatic probes along the direction of the shortand long-dimensions of the driver region. The plasma parameters will be compared with those obtained at the lowest area of the expansion bucket to analyze the plasma expansion properties from the driver region to the expansion region.