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.

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.

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 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.

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.

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.

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.

Software to Estimate Earthquake Spectral and Source Parameters

A software (EQK_SRC_PARA) has been developed to estimate spectral parameters of earthquake source spectrum, namely: low frequency displacement spectral level (Ω0), corner frequency above which spectrum decays with a rate of 2 (fc), the cut-off frequency above which the spectrum again decays (fmax) and the rate of decay above fmax (N). A Brune’s source model [1,2] that yield a fall-off of 2 beyond corner frequency is considered with high cut-off frequency factor presented by Boore [3] that fits well for frequencies greater than fmax. The software EQK_SRC_PARA is written in MATLAB and uses input data in Sesame ASCII Format (SAF) format. The obtained spectral parameters have been used to estimate source parameters (e.g., seismic moment, source dimension and stress drop etc.) and to develop scaling laws for the study region. The cut-off frequency “fmax” can also be studied and interpreted to confirm about its origin.

Variation characteristics of aquifer parameters induced by groundwater source heat pump operation under variable flow

The variation characteristics of aquifer parameters,induced by groundwater source heat pump(GWSHP) operation under variable flow,were theoretically analyzed through a case study,in which the characteristics of building air conditioning load were considered.The results,compared with the constant flow operation,indicate that the influence on the variations of porosity,hydraulic conductivity and confined water head is decreased by 48%,51% and 71%,respectively,under variable flow operation.The security of variable flow operation is superior to that of constant flow.It is also concluded that the climate region and function of the buildings are primary factors which affect the suitability of variable flow operation in GWSHP.

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.

Source Parameters Evaluation of the 28^(th) December 1908 and 16^(th) January 1975 Earthquakes in the Messina Straits(Southern Italy)

The Calabria-Peloritani Arc is one of the most seismically active regions of Southern Italy.The most damaging event(Ⅰ0=ⅪMSK)occurred on December 28th,1908 in the Messina Straits area.This event,which caused at least 85000 casualties and widespread damage,has been identified as the greatest European seismic disaster of the 20th century.

ThesourceparametersandsourcemechanismfortheJanuary5,1998JingyangearthquakewithMS=4.8

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

Interrelation between seismicity parameters and delimiting potential seismic sources in a seismic statistical region and its influence on seismic risk estimation

In the paper, we have discovered the abnormal area distribution features of maximum variation values of ground motion parameter uncertainty with different probabilities of exceedance in 50 years within the range of 100°-120°E, 29°-42°N for the purpose to solve the problem that abnormal areas of maximum variation values of ground motion parameter uncertainties emerge in a certain cities and towns caused by seismicity parameter uncertainty in a seismic statistical region in an inhomogeneous distribution model that considers tempo-spatial nonuniformity of seismic activity. And we have also approached the interrelation between the risk estimation uncertainty of a site caused by seismicity parameter uncertainty in a seismic statistical region and the delimitation of potential sources, as well as the reasons for forming abnormal areas. The results from the research indicate that the seismicity parameter uncertainty has unequal influence on the uncertainty of risk estimation at each site in a statistical region in the inhomogeneous distribution model, which relates to the scheme for delimiting potential sources. Abnormal areas of maximum variation values of ground motion parameter uncertainty often emerge in the potential sources of Mu greater than or equal 8 (Mu is upper limit of a potential source) and their vicinity. However, this kind of influence is equal in the homogeneous distribution model. The uncertainty of risk estimation of each site depends on its seat. Generally speaking, the sites located in the middle part of a statistical region are only related to the seismicity parameter uncertainty of the region, while the sites situated in or near the juncture of two or three statistical regions might be subject to the synthetic influences of seismicity parameter uncertainties of several statistical regions.

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.

Using Grey Relational Analysis to Evaluate Different Influencing Factors on Aquifer Parameters Variations with Groundwater Source Heat Pump Operation

The influence of particle migration induced by groundwater source heat pump( GWSHP) operation on aquifer parameters was theoretically analyzed. Then the sensitivity analysis of different influencing factors,which influenced aquifer parameters with GWSHP operation,was evaluated by using grey relational analysis method through a case study. The present study indicated that the erosion parameter, hydraulic conductivity and initial concentration were critical influencing factors which affected the aquifer parameters and the suitability of GWSHP. It also indicated that the discharge of single well should be controlled and the diameter of single well should be selected appropriately in order to decrease particle migration.

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 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 Log 10 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 dyne5cm （equal to an ML4.0 event）, stress drop is weakly dependent on M0; whereas, when M0>4×10 21 dyne5cm, stress drop is independent of M0.

Characterization of Source Parameters and Some Empirical Relations between Them for Kachchh Region, Gujarat, India: Implication of January 26, 2001 Bhuj Earthquake and Its Aftershock Sequence

Study of source parameters of small to moderate and large earthquakes is important to understand the differences and similarities between dynamic ruptures of different earthquakes and clarifying the scaling relations. In the present study, we have characterized source parameters and presented new and revised empirical relationships between various source parameters for Kachchh region of Gujarat, India to facilitate to draw first-order conclusions regarding the trends in the region. We have studied total 202 aftershocks of shallow-focus (hypo central depth less than 40 km) and moderate magnitude recorded over the Kachchh region during January 2001 to December 2012 by different seismological observatories of India Meteorological Department. We have adopted the spectral technique for source parameter estimation, where S-wave displacement spectra are considered and applied Fast Fourier Transform (FFT) to compute displacement spectra. We have followed the Brune’s source model for our estimation and the estimated values of source parameters show close approximation to the global values. While derived empirical relations between different source parameters, they demonstrate direct or inverse proportion to linear or power scale. Interrelation between seismic moment, rupture parameters, corner frequency and radiated seismic energy can be summarized as , , and ER μ M0 and ER μ Mw for our analysis. Stress drop distribution over the Kachchh region is very scattered and due to its peculiar behavior, it is difficult to derive its empirical relation with other source parameters. Sufficient accuracy on measuring source parameters like corner frequency stress drop, rupture dimensions, radiated seismic energy etc. helps to understand earthquake processes in the region. This is the first ever attempt to establish empirical relation between different source parameters for Kachchh region for longer aftershock sequence and they are useful to assess future earthquake potential over the region.

First Determination of Source Parameters of Moderate Earthquakes (4.1 ≤ M ≤ 5.1) in Morocco from Spectral Analysis

Recent installation of an array of broad band seismological stations in Morocco allowed us to study the records of five recent (2005-2008) moderate earthquakes (4.1 ≤ M ≤ 5.1) in order to determine their source parameters (seismic moment, fault slip, rupture area and stress drop) from P-wave spectra. We also studied the older Rissani events of 1992 using teleseismic data. Values of Mo, r, Δu and Δσ are, respectively, 1.1 × 1013 - 6 × 1016 Nm;0.50 - 3.9 km;0.8 - 5.8 cm and 0.3 - 1.49 MPa. The results are in accordance with the seismotectonic and geodynamic setting of Morocco as, for instance, the amount of slip along the faults with respect to the relative displacement of Nubia to Iberia (~4 mm·yr-1) determined from GPS data, taking into account the period of stress accumulation. However, some events show very variable corner frequency and low-frequency amplitude values which lead to considerably higher stress drop and fault slip values, especially at the nearest stations, which may reflect some site effects or uncertainties on depth and take-off angles.