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.

Hybrid Slip Model for Near-Field Ground Motion Estimation Based on Uncertainty of Source Parameters

The hybrid slip model used to generate a finite fault model for near-field ground motion estimation and seismic hazard assessment was improved to express the uncertainty of the source form of a future earthquake.In this process, source parameters were treated as normal random variables, and the Fortran code of hybrid slip model was modified by adding a random number generator so that the code could generate many finite fault models with different dimensions and slip distributions for a given magnitude.Furth...

Study on Inelastic Attenuation and Source Parameters of the Shidian Swarm in 2001

Based on horizontal-component digital seismograms recorded on 6 stations of the Yunnan Regional Digital Network, we inversed the inelastic attenuation in the source region of the Shidian swarm using the Atkinson method and the site responses of the 6 stations 200kin around the Shidian epicenters using the Moya method. The observational seismic waveform data were corrected by removing the propagation, instrument and site effects before the source parameters of the Shidian swarm in 2001 were determined using genetic algorithms. The results are as follows： （ 1 ） There is a linear relation between seismic moment and local magnitude. The seismic moment is between 1×10^12 and 10^14 N·m. The rupture radius of the seismic focus varies from 157m to 973m. The seismic moment and the rupture radius maintain a linear correlation. （2） The corner frequency increases as the seismic moment decreases. Based on the expression between corner frequency and seismic moment using least squares fitting, we can obtain the estimated value of the corner frequency. The time-varying value of the calculated corner frequency minus the estimated corner frequency shows that there were continuous high and low anomalies before the strong aftershocks. （3） The seismic stress drop is in the range of 0.07-1.55MPa. The stress drop seems independent of the local magnitude. The variation of stress drops is high before the occurrence of the strong aftershocks. （4） The depth of aftershocks is mostly in a range from 5kin to 10km, which means that energy release of aftershocks is mainly concentrated in this range of depth.

Research on the Q-value, Site Response and Seismic Source Parameters in Jiangsu and Its Adjacent Areas

Based on 49 digital seismograms recorded by 73 seismic stations in the Jiangsu Telemetered Seismic Network,the paper uses Atkinson's method to calculate the inelastic attenuation coefficient of the Jiangsu area. We find that the frequency-dependent Q in the Jiangsu region is Q( f) = 272. 1·f^(0. 5575). We also use Moya's method to invert the 63 stations' site responses. The results show that the site responses of the 25 stations in Jiangsu are approximately 1 at a range between 1Hz and 20 Hz, which is consistent with their basements on rocks. The response curves of the site responses of the 14 underground stations are similar to each other. Their site responses show an amplification at low frequencies and minimization at high frequencies. The calculation of the Brune model on the waveform data of M_L≥2. 5 earthquakes from Jiangsu Digital Seismic Network between October 2010 and May 2015 in terms of seismic source parameters of 58 seismic waves shows that there are good correlations between seismic magnitude and other source parameters such as seismic moment, source radius and corner frequency, while the correlations between seismic magnitude and stress drop,and stress drop and source radius are not so good.

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.

Study on the Source Parameters of Earthquakes in the Middle Eastern Area of North Tianshan in Xinjiang

By applying the empirical Green＇s function method, we first proved that the source spectra in the middle eastern area of North Tianshan follow the ω^-2 Brune model, and then, on this basis, obtained the source spectra for 105 earthquakes of ML 2.5 - 5.7 in the studied area upon removing the noise, instrument, propagation and site effects in the observational seismic S waveforms. Finally, we determined the source parameters such as the seismic moment, stress drop, source dimension, etc., based on the Brune model. The results show that there is a good linear correlation between scalar seismic moment, M0, and ML, which is Log10 M0 = 1.10 ML ＋ 17.20; The source radii range between 100- 1500m; the stress drop, 1-16MPa with a predominant range of 1 - 10MPa. The relationship of the seismic moment versus corner frequencies indicates that there may exist two source scaling, that is, when M0 〈 4×10^21 dyne. cm （equal to an ML4.0 event）, stress drop is weakly dependent on M0： whereas, when M0 〉 4×10^21 dyne·cm, stress drop is independent of M0.

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.

Study on the Inelastic Attenuation Coefficient,Site Response and Source Parameters in the Southeast Area of Gansu Province

Using digital seismograms recorded by the Gansu digital seismic network, the inelastic attenuation coefficient is calculated based on a genetic algorithm and the method proposed by Atkinson. Then, the site response and source parameters are investigated by the Moya method. The inversion results indicate the frequency-dependent inelastic attenuation, Q value, in the southeastern Gansu is estimated as Q （f） ： 404.2f^0.264 . Except for the Tianshui station, the site responses of the other stations do not show significant amplifications, which is consistent with their basement on rocks. The stress drops of all 39 earthquakes range between 1 × 10^5 and 7 × 10^6 Pa. We also found the dependence of corner frequency on seismic moment and seismic magnitude.

Temporal-spatial Variations of Dynamic Source Parameters in the Capital Circle Region and Its Surrounding Areas before the 2006 M_S5.1 Wenan Earthquake in Hebei Province

Source spectra,corner frequency and zero frequency amplitudes in near-source conditions were measured using waveform data from 989 earthquakes with magnitudes larger than ML2.0 observed by the Beijing Digital Telemetry Seismic Network in the Capital Circle Region of China and its surrounding areas from January 2002 to June 2006 by the Brune model.Relevant formulas that were used for the calculation of dynamic source parameters include rupture radius,seismic moment,seismic energy,stress drop,and apparent stress.Scaling relations and characteristics of temporal-spatial variations of these dynamic parameters before the MS5.1 Wenan earthquake in Hebei Province that occurred on July 20,2006 were analyzed.Results show that apparent stress,stress drop,and the ratio of seismic energy to the rupture radius had relatively high values in some areas before the Wenan earthquake.These high-value concentration areas were mainly distributed in the North China Plain seismic zone.As is seen from the time curves,parameters,such as apparent stress,stress drop,and ratio of seismic energy to rupture radius underwent significant ascending processes before the Wenan earthquake,but the variation in the corner frequency showed a descending trend.This result might be related to the enhancement of stress in the North China Plain seismic zone before the earthquake.

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.

The plate contact geometry investigation based on earthquake source parameters at the Burma arc subduction zone

Accurately characterizing the threedimensional geometric contacts between the crust of the Chinese mainland and adjacent regions is important for understanding the dynamics of this part of Asia from the viewpoint of global plate systems. In this pa per, a method is introduced to investigate the geometric contacts between the Eurasian and Indian plates at the Burma arc sub duction zone using earthquake source parameters based on the Slabl.0 model of Hayes et al. （2009, 2010）. The distribution of earthquake focus depths positioned in 166 sections along the Burma Arc subduction zone boundary has been investigated. Linear plane fitting and curved surface fitting has been performed on each section. Threedimensional geometric contacts and the extent of subduction are defined quantitatively. Finally, the focal depth distribution is outlined for six typical sections along the Burma arc subduction zone, combining focal mechanisms with background knowledge of geologic structure. Possible dy namic interaction patterns are presented and discussed. This paper provides an elementary method for studying the geometric contact of the Chinese mainland crust with adjacent plates and serves as a global reference for dynamic interactions between plates and related geodynamic investigations.

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

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

Effects of atomization parameters of dust removal nozzles on the de-dusting results for different dust sources

In order to obtain appropriate spray pressure and enhance the spraying and dust removal efficiency, various factors including the dust characteristics, nozzle spraying angle, effective spraying range, water consumption and droplet size are taken into account. The dust characteristics from different mines and atomization parameters of different pressure nozzles were measured. It was found that the internal pressure of coal cutters and roadheaders should be kept at 2 MPa, which could ensure large droplet size, large spraying angle and low water consumption and hence realizing a large-area covering and capture for large particle dusts. However, the external spray pressure should be kept at 4 MPa for smaller droplet size and longer effective spraying range, leading to effective dust removal in the operator zone. The spray pressure of support moving, drawing opening, and stage loader on a fully mechanized caving face and stage loader on a fully mechanized driving face should be kept at 8 MPa, under which the nozzles have long effective spraying range, high water flow and small droplet size for the rapid capture of instantaneous, high-concentration and small size dust groups. From the applications on the caving and driving faces in the coal mines, it is indicated that the optimization of spray pressure in different spraying positions could effectively enhance dust removal efficiency. Selecting appropriate nozzles according to the dust characteristics at different positions is also favorable for dust removal efficiency. With the selected nozzles under optimal pressures, the removal rates of both total dust and respirable dust could reach over70%, showing a significant de-dusting effect.

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.