Waveform inversion and analysis of an unusual earthquake swarm in the Boshan mining area, Shandong Province, China

Moment tensor solutions were retrieved for the earthquake swarm that occurred during November and December 2010 in the Boshan mining area, Shandong Province, China. The results showed that the double-couple components in the source mechanisms were higher at the beginning of the swarm and consisted mainly of shear faulting controlled by tectonic stress. The subsequent events had significant non-double-couple components, indicating tensile faulting. The double-couple components predominately presented as normal faulting and the P axes were orientated almost vertically. The slip vectors of the swarm events were relatively stable. With reference to the tectonic features near the epicenter, we concluded that the swarm was a result of subordinate fault motion related to the Wangmu Mountain fault and that high-pressure pore fluids played a crucial role in the activity of the earthquake swarm.


Introduction
In 2010, a series of earthquake occurred near Boshan, a mining area city in the past in Shandong Province, China. The earthquake sequence continued over 1 month, from November 24, 2010 to December 27, 2010 and were measured by the Shandong seismic network. The tremors from these earthquakes propagated widely through these mining areas and affected almost all the villages in Boshan, with consequent social impacts. There were six earthquakes of M L C 2.0 and the maximum event of M L 3.4 occurred on December 20, 2010. This event had the highest magnitude of all the earthquakes that have occurred in the mining areas of Shandong Province in recent years.
The term ''mine earthquake'' is used here to refer to all earthquakes occurring in mining areas, whether they were caused by rock bursts, mine collapses, or tectonics. Mine earthquakes can be divided into three types: (1) earthquakes in which a rupture occurs in the surrounding rock mass of the mining areas-these are caused by the adjustment of stress in the Earth's crust and result in slippage along a few minor faults; (2) shocks caused by mine collapse, roof falls, or rock bursts; and (3) shocks caused by the detonation of gas or powdered coal in mine tunnels (Zhang et al. 1993). The source of the earthquake is obviously different in each of these three categories. The source type can be indicated by the moment tensor. For earthquakes caused by tectonic activity resulting from stress adjustment in mining areas, the source mechanism follows the double-couple (DC) shear dislocation model and is close to that of natural earthquakes. The source of earthquakes caused by mine collapse, roof falls, or rock bursts is primarily uniaxial compression, and there is a large compensated linear vector dipole (CLVD) component in the moment tensor solutions and a characteristic implosion mechanism resulting from the sudden closure of goafs. For mine earthquakes resulting from the detonation of gas or powdered coal, the explosive source (i.e., the isotropic volume component) plays a large part in the source mechanism. Some studies have shown that many mine earthquakes have non-DC mechanisms and may contain a large volume or linear component under some circumstances (Gibowicz and Kijko 1998;McGarr 1992a, b;Feignier and Young 1992;Phillips et al. 1999;Fletcher and McGarr 2005;Stec 2007;Julià et al. 2009). There have been a number of studies in China in recent years of the source mechanisms of mining-induced earthquakes He et al. 2007;Chen et al. 2009), although there has been little specific research into the moment tensor of mine earthquakes.
Non-tectonic earthquakes and natural earthquakes are technically distinguishable as they differ in many aspects, such as the first motion of the seismic waves, the ground motion frequency, and the rate of attenuation. Many researchers have conducted in-depth research into the differences between non-tectonic and natural earthquakes using seismic waveforms. For example, Lin et al. (1990) compared the seismic phase characteristics of collapse earthquakes and tectonic earthquakes and suggested that collapse earthquakes have a different wave propagation velocity, attenuation, and amplitude. Gibowicz (1995) studied the relationship between the seismic moment of mining-induced events and the P and S wave corner frequencies in various regions using statistical comparisons. Oye et al. (2005) studied seismic source parameters such as the coda attenuation and apparent stress of a mine earthquake swarm in Finland using the spectral integral and multiple empirical Green's functions. He et al. (2006) and Liu et al. (2005) solved the recognition factors for mine and tectonic earthquakes based on wavelet packet decomposition of the waveform. Shen et al. (2006) studied the rupture and nucleation processes of mine earthquakes based on micro-tremor records and micro-seismic waveforms in mining areas. Zhang et al. (2009) compared and analyzed the time-frequency characteristics of the spectra of mine and tectonic earthquakes in the Three Gorges area of China.
During emergency work after this swarm of earthquakes in 2010, we noticed that these particular mine earthquake events differed from ordinary mine earthquakes and that there was a need for further analysis. Moment tensor analysis is a powerful tool for source mechanism studies. Therefore, we analyzed and inverted all the digital waveform recordings of these mine earthquake swarm events to determine the source moment tensor solutions.

Data and information
2.1 Distribution of the earthquake swarm and seismic stations A total of six events with M L C 2.0 were recorded in the Boshan earthquake swarm. In addition, one event with M L 2.7 occurred near the epicenter of the Boshan swarm on September 12, 2010 and, while the swarm was active, one event with M L 2.7 occurred 16 km to the south of the epicenter on November 25, 2010. As these two events were closely related to the Boshan swarm both tectonically and in time, they were also included in our study. All the broadband digital seismic stations within 150 km of the epicenters of the largest magnitude events in the swarm recorded relatively clear waveform signals. Figure 1 shows After several mine earthquakes in late November 2010, the Earthquake Networks Center of Shandong Province and the Earthquake Administration of Zibo City set up four mobile seismic stations in the vicinity of the epicenter. These four stations plus the BSH station, which is \10 km from the epicenter, formed a temporary seismic network of five stations and recorded several mine earthquake events. Figure 2 shows the relocated epicenters of the Boshan swarm, and the temporary stations at the Boshan area.
The data from the temporary stations show that the S-P value is [0.5 s at the closest station L02; the relocated results show that the focal depths of all these events ranged from 4 to 5 km. The depths of mine earthquake events such as rock bursts and collapses caused by industrial mining activities are usually of the order of hundreds of meters, although they may exceed 2 km in a few deep mine areas [e.g., deep gold mines in South Africa (McGarr et al. 2009)]. However, according to data from the Zibo Mining Group Co. Ltd (Jiao et al. 1990), the Zibo coalfield contains Carboniferous-Permian coal measures with [20 layers of coal, 12 of which are minable in part or in whole; the main coal-bearing strata are 250 m thick and the mining depths are normally in the range of 200-300 m, with a maximum of \400 m. Therefore, the cause of this set of events may not be due purely to mining operations, such as rock bursts. As the epicenter of this earthquake swarm lies in the coal mining area, some researchers have called it a ''mine earthquake'' swarm; however, we refer to it as an ''earthquake swarm'' in this study. Figure 3 shows the velocity model used in our study, which is derived from geophysical exploration in Shandong area (Jiang et al. 2000). The epicenter locations obtained through relocation was used in the inversion of the source moment tensor. The search was conducted in the depth direction with a starting depth of 1 km and a step length of 1 km. In total, there were 23 search depths to a maximum of 23 km.
2.2 Waveform time-frequency analysis of Boshan swarm Zhang et al. (2009) have shown that there are a few clear differences between the waveform time-frequency analysis results for mine and tectonic earthquakes. First, the main frequencies of the waves from tectonic earthquakes are higher than those of mine earthquakes. Second, the frequency components of waves from tectonic earthquakes are richer than those of mine earthquakes and their energy density spectrum is spread over a wider range of frequencies-the bandwidth of tectonic earthquakes is four times as wide as that of mine earthquakes. Third, the peak energy intensity of mine earthquakes over the entire wave train appears much earlier than that of tectonic earthquakes, although the attenuation in energy is more rapid. The timefrequency spectrogram of the Boshan earthquake swarm events shows that the events have distinctive features. The events share some characteristics with general mine earthquake events, such as a low main frequency and rapid attenuation. However, the energy density spectrum of the events is widely spread along the frequency axis and the dominant frequency of the S waves extends from about 7 to 17 Hz (Fig. 4), indicating the rich frequency components of the event.
3 Theory and methods

Moment tensor decomposition of mine earthquakes
In mathematics, the moment tensor describes the effectiveness of any type of source and can be simply divided into three parts corresponding to different types of sources: an isotropic (ISO) component, a CLVD subcomponent, and a DC component. The DC tensor describes the shear slip in isotropic media along fault planes and the ISO component corresponds to explosive or implosion sources. The ISO and CLVD components are referred to as the non-DC components and are attributable to the opening or closing of faults. Many studies on the source moment tensor of mine earthquakes have shown that, as a result of the volume change caused by the closure of goafs, the source mechanisms of mine earthquakes (also referred to as ''shock bump'' or ''rock burst'') can be decomposed into a combination of two modes: shear slip accompanied by a certain amount of a co-seismic implosion component. McGarr (1994) studied the mechanisms of mine earthquakes and found that the combination of a shear faulting model plus an implosion component could explain the waveforms of similar events observed in deep gold mines in South Africa (McGarr 1992a, b). Increasing numbers of studies are showing that many shallow, naturally occurring mine collapse events have an implosion mechanism as the seismic source (Gibowicz and Kijko 1998;Miller et al. 1998). Yang et al. (1998) studied explosion-induced mine collapse. The results of source moment tensor inversion showed that the source mechanism of collapse can be expressed by horizontal opening and an implosion. Analysis showed that the unique source characteristics inducing collapse result from the fact that, unlike synchronous collapse, ordinary collapse is initially caused by extension fractures.
Unlike tectonic earthquakes caused by ordinary shear faulting, the fault slip vectors of mine earthquake events are not limited to the fault planes and may deviate from them. For normal seismic sources, the moment tensor can be expressed as a function of a slip vector. Aki and Richards (1980) devised a formula for the point source moment tensor of faults in isotropic media: where k and l are the Lamé constants of the faults, d kl is the Kronecker symbol, ½u is the slip vector, and n is the normal fault plane. If then the moment tensor may be in the following form: where a½Àp=2; p=2 is the included angle between the slip vector and the fault plane. When a [ 0, the source is an extension fracture; when a\0, it is a compression fracture; and when a ¼ 0, it is pure DC shear faulting. When using the principal stress axis coordinate system, Eq. (3) can be written in the following diagonal form: Then, the matrix trace is as follows: Under the constraints that have physical significance, (l [ 0; k=l [ À 2=3), the trace of matrix M is only related to the included angle a between the slip vector and the fault plane. Therefore, for an extension-type source, the trace of matrix M is positive; otherwise, it is a compression-type source. We designate this by k=l ¼ j as the moment tensor can be decomposed into We can obtain the various components in the moment tensor in the following proportions: where M jmaxj is the maximum absolute value of the characteristic value of the moment tensor M and e is the DC component corresponding to the CLVD component, defined as where M Ã jminj and M Ã jmaxj are, respectively, the minimum and maximum absolute values of the characteristic value of the deviatoric moment tensor M Ã .
Based on Eqs. (5), (7), and (8), Vavryčuk (2001) devised the following equation: We can infer from the formula that, based on the proportion of different moment tensor components, the specific value of the Lamé coefficient of the fracture fault and the dip of the slip vector of the fault in relation to fault plane can be calculated.

Method for moment tensor inversion
We used ISOLA package to do the waveform moment tensor inversion in time domain Zahradník 2013, 2008). For small earthquakes, the synthetic seismogram can be approximated with a combination of basic seismograms corresponding to six basic types of source mechanism: where E i ðtÞ; i ¼ 1; 2; . . .; 6, is six basic seismograms; a i is the coefficient to be solved, which is related to the moment tensor, M ij , in the following geographical coordinates (north, east, and up are positive): The use of the least-squares method to solve the inversion problem minimizes the degree of misfit between the observed waveform and the synthetic waveform. The misfit function is defined as follows: where i = 1, 2,…, n; n = number of stations 9 number of vector components 9 number of frequency of use. Green's functions were calculated using Bouchon's discrete wavenumber method (Bouchon 1981) and the original waveform recordings were calibrated experimentally. The appropriate frequency band was selected for band-pass filtering, e.g., 0.10-0.25 Hz for M L 3 event. 5 % of cosine taper is used for FFT. A grid search method for a preset trial hypocentral location and time displacement was used to determine the optimum result.

Results
We obtained the source moment tensor solutions of eight events in this earthquake swarm (Tables 1, 2, 3). The j value in Table 2 shows that all the inversion results are in a physically significant range. Figure 5 shows a beach ball representation of the moment tensor solutions.
4.1 December 20, 2010 M L 3.4 event As this event was the largest, all the broadband digital seismic stations within 150 km of the epicenter recorded clear waveform signals. We chose the data from six stations (LQU, WLS, YSH, NLA, XIT, and TIA) with epicentral distances within 100 km for the moment tensor studies.
The calculated scalar seismic moment for this event was 1.588 9 10 14 N m, the focal depth was 5 km, and the optimized DC source mechanism solutions were nodal plane A strike 302°, dip 38°, rake -49°; and nodal plane B strike 75°, dip 61°, rake -116°; see Fig. 5 for moment tensor solutions.
The correlation coefficient of the optimum synthetic waveform solution was nearly 0.8, showing a good inversion result. The moment tensor solution results showed that the DC component in the seismic source was relatively small, whereas the non-DC components were about 70 % and the ISO and CLVD subcomponents were dramatic, showing that the swarm events were not shear dislocation induced purely by tectonic causes. The source mechanism solutions in Fig. 6 show that the seismic source contains a large uniaxial tensile component. Figure 7 shows the match between the synthetic and filtered observed waveform at five stations (WLS, YSH, LQU, XIT, TIA). The three-direction average correlation coefficients of these stations are fairly good, supporting the reliability of the moment tensor solutions obtained in this work.
To explain the authenticity of the non-DC components in the moment tensor solutions, inversion was carried out again with the source mechanisms constrained to pure DC solutions; the solutions at different depths are shown in Fig. 8. This figure shows that the degree of waveform correlation of the optimum pure DC solutions at various depths is significantly lower than the moment tensor solutions, explaining the authenticity and reliability of the non-DC components in the moment tensor solutions. Table 2 shows that, except for the three initial events, the proportion of the DC components in the events was very small (\60 %). Trifu and Shumila (2002) found that large non-DC components appeared in some mine   were 15 %-20 %; their analysis showed that the moment tensor results were reliable. Dahm et al. (2000) studied the accuracy of the seismic source mechanism by comparing the absolute and relative moment tensor solutions; they estimated the error caused by the medium model and inaccurate amplitude data. Dahm et al. (2000) suggested that the significantly high level of similar non-DC components was not caused by inaccurate moment tensor inversion and that a large proportion of these non-DC components reflected the special properties of the seismic sources or media. Miller et al. (1998) used statistics to show that many mine earthquakes contain non-DC mechanisms and that the first motion is inconsistent with a normal quadrantal distribution. Stickney and Sprenke (1993) studied 21 mine earthquakes (rock burst) in the Coeur d'Alene mining area, Idaho, USA, and found that 90 % of the P-wave first motions observed were downward and that the first motions of 11 of the 21 events were negative. The first motion polarities from 25 stations for the Dec 20, 2010 M L = 3.4 event were obtained from the waveform data recorded by the Shandong Seismic Networks. In contrast to ordinary mine earthquakes, most of the first motions were upward and even the downward first motions observed by three stations were not obvious; the possibility of misreading cannot be ruled out. Figure 9 shows the correspondence between the source mechanism solutions obtained by MT inversion and the first motion polarities; the results include nine stations with an epicentral distance within 200 km, but which could not determine a clear first motion. Figure 9 shows that because the volume components in this seismic source were large, in the first motion expansion zone in which the first motion should originally be downward, the first motions observed by many stations were blurred. Many positive first motion polarities were observed and this further supports the reliability of solutions from moment tensor inversion.

Discussion
For non-DC components in seismic moment tensor solutions, some seismic moment tensor analysis has shown Trifu and Shumila 2002;Horálek et al. 2010) that, under certain circumstances (e.g., mine earthquakes due to reservoir filling, water injection in oil production, volcanic activity and geothermal eruption, landslides, collapses, and subsidence), the non-DC components (the CLVD and ISO components) may form a large proportion of the total. Julian et al. (1998)   force system and systematically analyzed the possible physical processes of non-DC earthquakes. They gave the theoretically reasonable explanations that the non-DC rupture process may have its roots in the shear faulting and extension fracture of complex faults, or in shear faulting or polymorphic phase transition in anisotropic media.
The results for the solution of typical mine earthquake focal mechanisms induced by mining operations often show a large implosion component and, as the goafs collapsed, non-DC components in moment tensor solutions. The ISO and CLVD components are negative, indicating the existence of a large amount of uniaxial compression accompanying volume reduction (McGarr 1992a, b;Julià et al. 2009). The computation of this Boshan earthquake swarm events showed that, at the beginning of swarm activity, the DC components formed a large proportion of the source mechanisms ([75 %), whereas the non-DC components may be attributed, to some extent, to inversion error due to constraints in the observed data. However, the moment tensor in the later events showed very large non-DC components and the ISO and CLVD components were both positive, indicating that the seismic sources were tensile cracks.

Stress state and orientation of P, T axes
Previous studies of the source mechanisms of typical mine earthquake events have shown (McGarr 1992a, b;Julià et al. 2009) that the maximum stress in mine earthquake sources is compressive stress in a nearly vertical direction, accompanied by a reduction in volume and normal faulting; this is consistent with the closure of goafs as a result of gravity. Figure 10 shows the orientation of the P, T axes in eight of the mine earthquake events studied in this work. Except for one tectonic earthquake event (201011241356), the P axis in the other seven events is orientated in a nearly vertical direction and the DC components show normal faulting. However, in contrast with the findings of Julià et al. (2009), the maximum stress in these events was tensile stress. The T axis was orientated to the southwest in the projection direction of the focal sphere, indicating that the source rupture in these mine earthquake events was largely caused by the action of tensile stress in the southeast in a nearly horizontal direction. The projection of the slip vectors on the focal sphere shows that, apart from the event 16 km south of the earthquake swarm on November 25, 2010, the projection direction of the seismic source slip vectors of the earthquake swarm was relatively stable. Feignier and Young (1992) created many non-DC mechanism earthquakes by explosions in a 3.5-m-wide mine tunnel cut into complete granite in the Manitoba URL Laboratory, Canada. They calculated the moment tensors of 33 recorded microquakes with moment magnitudes ranging from -2 to -4 and found that tensile events mainly occurred around the caved areas at the top of the mine tunnels and that the implosion events usually occurred in front of active surfaces. Their computation showed that the maximum principal stress in the area was characterized by a positive tensile stress and that the tectonic stress background could explain the existence of a tensile source.
In Boshan District, Zibo City is considered to be a seismic window for Shandong because there are frequent small swarms of earthquakes and microquakes. This area often had an obvious concentrated activity of small earthquakes before the occurrence of moderately strong earthquakes in the zone (within about 500 km). (Su 1997) Boshan District lies around the top of the arc of a horseshoe-shaped structure in west Shandong where stress tends P T α to concentrate. Faults are developing in this structure, mostly in the form of high-level normal faults; reverse faults are rare. The data show that (Jiao et al. 1990) the epicenters are mainly in the Zibo syncline; this strikes c. 40°-50°NE and trends to NW at an angle of 15°-18°. The direction of tensile stress caused by pure rock gravity therefore does not fit with the computation. The epicenter of this earthquake swarm was mainly the Wangmu Mountain fault in a northeasterly direction. According to Shi et al. (2009), the Wangmu Mountain fault trends west with a dip of about 75°and the activity is mainly in the form of normal faults. Although these fracture properties are close to our moment tensor for tensile cracking, the strike of this fault is inconsistent with the inversion result. The Yuwang Mountain fault lies 6 km to the west of this earthquake swarm and strikes 2°-10°NE, trending approximately east with a dip of 60°-80°, and is mainly characterized by a tensile structure. About 5 km to the south of the epicenter, the Shentou-Xihe fault strikes 40°-60°NE, trending to the SE with a dip of 60°-70°. Therefore, within the small structural unit enclosed by the above three areas of extension tectonics, there is also tension in the main stress direction. The small earthquake moment tensor we obtained through inversion is consistent with this analytical result. According to data from the Mining Bureau of Zibo (Jiao et al. 1990), there are two groups of secondary faults in research area besides large fault shown in Fig. 2. The 1st group spread approximately in EW direction; the 2nd group arranged in NW direction with normal faulting. As no more accurate data about these secondary faults are available, it is impossible to judge which group of faults these swarm events are related to. However, we are certain that the swarm is caused by fracture tectonic activity related to the Wangmu Mountain fault under the conditions of stress adjustment in this small area of Boshan District.

Analysis of causes for Boshan earthquake swarm activities
There are many possible sources for non-DC source mechanisms. If multiple shear faulting occurs in conterminous, but different, fault planes, or the shear faulting of non-planar faults is mistakenly considered as an event, then the real DC tensor and moment tensor of the combined events are usually non-DC tensors (Frohlich 1994). In this case, shear faulting may generate a non-DC mechanism; however, there is no isotropic volume component. As the trace of the combined moment tensor is equal to the sum of the traces of its components, the combination of DC tensors cannot generate non-zero ISO components (Horálek et al. 2010).
The heterogeneity or anisotropism of the underground media may also result in the appearance of non-DC components in the moment tensor solutions of shear faulting. The experiments of Kawasaki and Tanimoto (1981) and the studies of Vavryčuk (2002) have shown that the anisotropism and heterogeneity of the underground media are not the main causes of higher non-DC mechanisms. Considering that the data from field mobile stations were used in our studies of the earthquakes that occurred in the period December 2-5, 2010 within a small range with a radius of \10 km, the properties of the underground media can be assumed to be uniform and the influence of lateral heterogeneity is negligible.
Theoretical and experimental studies on the source causes of tensile earthquakes have been reported previously (Walter and Brune 1993;Scholz 1990). Studies have shown that tensile earthquakes usually occur in areas with large amounts of underground fluids, such as in a geothermal field or near volcanoes Baisch et al. 2009). Julian et al. (1998) proposed that highpressure underground fluids are the primary cause of tensile cracking because the existence of fluid pressure may offset rock stresses in the Earth's interior. Fluids under high pressure may rapidly fill the rupture zones when the rocks burst, giving rise to large tensile components (Fig. 11). Walter et al. (2009) studied the source mechanisms in icequake events caused by the opening of crevasses on the surface layer of the Gornergletscher glacier, Switzerland, using the full waveform inversion method. The moment tensor results showed extensional fractures, indicating the reliability of moment tensor inversion. Walter et al. (2010) showed that the involvement of high-pressure liquids in the rupture process may lead to tensile cracking by studying the moment tensor solutions in icequake events that occurred at the base of glaciers.
In Boshan District, Zibo City has been long known as the Second Spring City of Shandong because there are four groups of springs (Shentou, Qiugu, Zhulong, and Liangzhuang) and more than 40 recorded springs. In addition, a number of geothermal areas and abnormalities in P is the fault plane, v is the normal to the fault, ½u * is the slip vector, and a is the inclination of the slip vector from the fault plane hydrochemistry exist in this area (Yang et al. 2005), which matches the observations in earlier studies that tensile earthquakes with non-DC mechanisms mainly occur in areas where there are large underground geothermal fluids or volcanoes. Therefore, the large non-DC components in the moment tensor solutions of the Boshan swarm may be caused by the involvement of underground fluids in the rupture process in the source area.
Moment tensor analysis showed that (except for the three events that occurred on September 12, November 24, and November 25, 2010) the mechanism in the early stages of the swarm was primarily tensile cracking. As massive underground fluids are present in this area, the main cause of the earthquakes was shear faulting at the beginning of swarm activity. After fracturing, underground fluids with a high pore pressure filled the burst area, reducing the friction stress and making a second rupture easier, which resulted in a continuously active earthquake swarm. The high pore pressure also gives rise to the source mechanism of tensile cracking. The results of rock crushing tests support tensile cracking (Byerlee 1990). Š pičák and coworkers (2000) and Horálek et al. (2002) studied the West Bohemian region of the Czech Republic and they proposed that high fluid pressures may be a key factor in the frequent swarm activities in this region. This feature is similar to that of the Boshan seismic window. Therefore, we conclude that underground fluids with high pore pressures were involved in the rupture of the hypocentral region and played an important role in the activity of the Boshan earthquake swarm.
It is usually thought that full moment tensor solution retrieved by inverting seismograms may include errors due to incomplete knowledge of the structure along the propagation path. Studies in gold mine in South Africa shows that the full MT solutions inverted from close stations within a few kilometers to the source are believable, even in homogeneous model (Š ílený and Milev 2006). In our study, the farthest temporary station to the swarm is no more than 10 km. Despite a few local stations, the final result should still reliable.

Conclusions
Using data from the Shandong seismic station networks and waveforms recorded by the temporary station network, we determined the source mechanisms of eight events in a mine earthquake swarm in November and December, 2010 in Boshan District using directional instrumentation and the time domain moment tensor inversion method. We also determined the source mechanism of an earthquake of M L 2.7 that occurred on September 12, 2010. The following conclusions were drawn.
(1) At the beginning of the swarm activities, the DC components formed a large proportion of the source mechanisms, which may be primarily attributable to shear faulting under tectonic stresses. However, the moment tensors in the following events showed very large non-DC components and the ISO and CLVD components were both positive, indicating that the seismic source was tensile cracking.
(2) The source mechanisms of the swarm events showed that-except for the earthquake event 16 km to the south on November 24, 2010, which showed reverse faulting-the DC components of the events showed normal faulting. The P axis was almost vertical and the maximum stress was tensile; the projection direction of the source slip vector of the swarm events was relatively reliable. (3) We analyzed the epicentral tectonic features and suggest that this swarm event was caused by fracture tectonic activity related to the Wangmu Mountain fault undergoing stress adjustment in the small area of Boshan District.
Based on analyses of the reliability of the moment tensor solutions, we proposed that underground fluids with high pore pressures were involved in the rupture process in the source area and played an important role in these Boshan swarm activities.
The results obtained by this study, constrained by the original data and velocity model, showed that the average waveform correlation coefficient of an individual event (e.g., 11291539) was not high; however, the results are still acceptable for microquakes of M L 2 (Fig. 1).
Although these findings showed that the source type at the beginning of the Boshan earthquake swarm activities was mainly fault rupture under tectonic stress, these activities occurred in a mining area where coal goafs have altered the stress distribution inside the underground formations. However, we can make no clear judgment about whether mining activities have caused these swarms. The moment tensor results for this swarm show that the source activities can be divided into two stages. Rupture in the earlier stage may have an influence on the activity in the next stage, although the physical mechanism of their interactions still unclear and requires further study.
creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.