Underground mining of coal, in aseismic areas like the Upper Silesia Coal Basin (USCB) (Figure 1), leads to a disturbance of the rock mass equilibrium. Furthermore, it triggers or induces the accumulated elastic energy release as an earthquake that may result in perceptible shaking on the surface which can even damage the infrastructure. For hazard assessment purpose, it is necessary to predict scale of the ground motion based on the Ground-Motion Prediction Equation, GMPE [1,2,3,4]. The GMPE, in addition to the parameter representing the source (energy) and path (distance), takes into account also the influence of the surface layer at the measurement site (amplification) [1, 2]. The highest attenuation takes place in a medium with inelastic properties, i.e. loose rocks [5,6,7,8]. However, for specific geological conditions of near surface layers, the amplitude amplification may occur for the corresponding resonance frequency of loose layers which is related to the velocity of seismic waves and the thickness of these layers .
Therefore, the main aim of the research was to find the best model of GMPE relation in the mining area of “Ziemowit” Coal Mine, located in the Upper Silesian Coal Basin in the central area of the Main Syncline (Figure 1) representing one of the most active mining areas . A model of a forecasting equation was proposed which provides for local effects (amplification coefficient) at the location of seismometer stations. In addition, the parameter h, being the estimated depth of exploitation, was also taken into account in the calculation of the model [3,5].
2 Site characterization
The “Ziemowit” Coal Mine is located in Bieruń and Lędziny commune, and it belongs to the Polish Mining Group. The mine includes two mining areas: “Lędziny I” and “Imielin I” and their mining areas cover the entire city of Lędziny and partially: Imielin, Chełm Śąlski, Bieruń, Tychy, Katowice and Mysłowice. The “Lędziny I” area is located in the central part of the Main Syncline in USCB where the hanging wall block of “Książęcy” normal fault is present. The layers are inclined in the South-East direction at an angle of 2°-8°. Here, the coal seams are cut by a large number of faults that divide them into technical blocks . Moreover, the major tectonic disturbances in this area are represented by: the Lędziński fault zone, the Imieliński-Wanda fault zone, the fault zone at the contact of Lędziński and Smardzowicki faults. The minor tectonic zones are represented by Chełmski fault, Eastern fault, Western fault, Piastowski fault and Ławecki fault (Figure 2). The “Imielin I” mining area is located on the border of the central and eastern part of the Main Syncline (the hanging wall block). Layers extend in the directions from S-N to NWW-SEE. The dip varies from 3° -10° and its azimuth changes from the north to the south Main Syncline. In the northern part, it is rather S, E and SE, but in the southern part it is SW . The coal seams in this mining area are cut by many faults with throw reaching values of 160 m. The largest dislocations have the NW-SE course and NS directions creating the Imieliński fault zone .
The overburden consists of the formations of the Quaternary and Tertiary (Table 1) . Quaternary deposits are present in entire study area and cover Tertiary and Triassic layers which form the irregular arranged blocks. Below those rock, the Carboniferous creates a rigid basement. However, in some cases the Tertiary or Triassic formations are the rigid basement (Table 1).
Outline of the geological structure in the area of Ziemowit Coal Mine 
|Age||Type of rocks||Thickness||Range|
|Quaternary (fluvioglacial deposits of the Lower Pleistocene and Holocene:)||gray and yellow medium- and fine-grained sands, and gray moraine clays with gravel and silt inserts||the thickness varies from 0 (Bieruń area) up to 5 m in the eastern and north-eastern part of Ziemowit mining area; the largest thickness occurs in the southern part of the area and is 23.5m (near Chełm Wielki)||Quaternary covers the entire area.|
|Tertiary||these layers are represented by gray-green clays, sandy clays and marls||the thickness varies from 0 to the maximum value of 328.4m||appearing irregularly on uneven surfaces of Carboniferous and Triassic deposits|
|Triassic (Buntsandstein, Ret and Muschelkalk)||sandy and silty sediments with a thickness of 12-23m below which there are Muschelkalk deposits (marls) and Ret dolomites as well as gray (pellet) limestones, crinoid limestones with marl inserts.||the thickness ranges from 0m - 6m up to 54.4m and reaches a maximum of 122.1m||it lies irregularly on the Carboniferous, forming blocks, it is present mainly in the eastern and north-western part of the mining area|
|Carboniferous||Libiąż layers (Westphal D), Laziska layers (Westphal C), Orzesze layers (Westphal B)||-||A solid rock basement covering the entire area|
The seismic stations located (Figure 2) in the northern part of the study area (Dziećkowice, SUW, Imielin, W-II) are installed in similar geological conditions, represented by approximately 30 m thick formation of the Quaternary sandy deposit (Table 2). Below the Quaternary complex, there are Tertiary sediments composed mainly of clays and silts, their thickness reaches 90m and they cover Triassic marls formation. The central part of study area, where Rubinowa, Ziemowit, Goławice, Pompownia stations are located, is characterized by ca. 15 m sandy Quaternary deposit covering very thick formation of Tertiary formations (thickness reaches 170 m) (Table 2). In the vicinity of Fawent and Chełm stations, the Triassic marls are very shallow – ca. 25 m below the surface(Table 2), and they are covered by a 20-meters layer of sands (Quaternary) and a thin 5-meter layer of silts (Teriary) (Table 2). The southern part of the study area (Czerniny and Kopciowice stations) is located in the area where the Quaternary deposit is about 15m thick and covers a 35-meters formation of Tertiary sediments and Triassic marls (Table 2). Specific geology can be found below the Zapora station, there are 9-meter thick Quaternary sediments (sands and clays) covering a 9-meter thick Tertiary silt and this last layer is situated on Triassic dolomites (Table 2) .
The coordinates of seismic station in the Sucha Góra system together with the number of registrations taken for the purpose of the analysis and geological settings below stations 
|No.||Seismometer station||Coordinates in the Sucha Góra geodetic system||Number of data recorded at each station||Thickness [m] (main type of rock)|
|X||Y||Elevation [m a.s.l.]||Q||Tr||T|
|1||Fawent||–32 360||22 170||245||2943||22.0 (sands1)||4.0 (silts)||65.0 (limest.)|
|2||Goławiec||–33 878||19 987||240||1614||13.4 (sands2)||142.6 (silts)||12.7 (marls)|
|3||Chełm||–33 272||23 082||244||4497||22.0 (sands1)||4.0 (silts)||55.5 (limest.)|
|4||Imielin||–29 253||22 119||263||397||30.0 (till+)||210.0 (siltst.3)||20.0 (congl.)|
|5||Szyb W–II||–30 472||22 091||253||365||28.5 (sands+)||232.5 (siltst.3)||–|
|6||Rubinowa||–30 725||20 972||261||757||15.1 (till)||-||-|
|7||Pompownia||–31 485||24 626||235||1397||12.0 (sands+)||171.9 (marls4)||-|
|8||Zapora||–33 077||24 668||236||2219||9.0 (sands+)||9.0 (silts)||11.0 (limest.)|
|9||Ziemowit||–31 724||18 059||257||221||15.0 (till+)||-||-|
|10||MSK||–31 020||23 737||244,5||506||12.0 (sands+)||171.9 (marls4)||-|
|11||Kopciowice||–35 057||23 369||245||237||39.0 (sands)||19.0 (till)||53.5 (limest.)|
|12||Czerniny||–34 324||21 813||239||230||16.0 (sands)||34.5 (silts)||48.0 (limest.)|
|13||SUW||–28 652||24 415||260||114||31.0 (sands+)||89.7 (silts)||19.7 (dolom.)|
|14||Dziećkowice||–25 961||24 304||260||44||30.5 (sands+)||89.7 (silts)||19.8 (dolom.)|
3 Data and methods
3.1 Seismic data
In the period from January 2010 to December 2016, a total of 15 541 seismic events were recorded by 14 seismometer stations (triangles in Figure 2). The coordinates of the stations in the Sucha Góra geodetic system and the number of data taken for analysis are presented in Table 2, and their location is shown in Figure 2. The measurements were carried out using the AMAX-GSI equipment, that is a multichannel network of seismometers recording accelerations of seismic waves on the surface. Each station is composed of a three-component piezoelectric accelerometer integrated with a preamplifier, amplifier and a digital recorder connected to a PC with the SEJSGRAM and AMAX software installed and matched with a 16-bit AC converter card. The accelerometer is characterized by frequency range of 1-100Hz and dynamics of 70-80 dB with the maximum value of recorded acceleration equal to 3m/s2. The positioning of the sensors ensures that they have good contact with the ground. Their installation place depends on space availability and they were fixed either to the load-bearing wall on the ground floor of the building or to the foundation wall in in basement. This network allowed to collect the PGA data from the studied period 2010-2016.
The energies and locations of the recorded event were calculated using the Seismological Observation System (SOS) developed by the Central Mining Institute in Katowice. The SOS consists of the network of 37 underground geophones and 6 surface seismometers . Determination of energy and localization of mining tremors were carried out using the MultiLok program provided by the Central Mining Institute as well.
Both these systems (AMAX and SOS) allowed to complete the seismic catalog that was used in the GMPE model calculation. The catalog contained information on: time of event occurrence, location of the event (Sucha Góra coordinate system), event energy, and value of PGA recorded by the surface stations.
3.2 GMPE model and evaluation
The GMPE model, chosen for calculation, was proposed by Lasocki  and it represents a regression model considering energy, mechanisms of amplitude decrease over distance (the geometric spreading), and site effects (relative amplification coefficients). The regression model can be expressed as [1, 3]:
where α, β and γ are regression parameters corresponding to energy E and geometric spreading at the wave path (distance), respectively, the distance between source and site is represented by
The regression model was obtained using the Normal Equation System to find a solution of multiple regression. Statistical assessment of the model was also carried out and the fitting quality of the model was confirmed by: determination coefficient R2, dispersion plot, standard error of estimate SEE and analysis of residuals . Relative amplification coefficient is logarithmized, thus to obtain a full value, the ai is raised to power. Moreover, this value is relative and is calculated as a difference between the considered value and the reference value as follow :
In this case, the lowest ai is assumed as the reference value .
4 Results and discussion
The best model of GMPE was searched iteratively with respect to the h parameter representing depth. The h-value was changed from 200m to 1200m below surface with the step of 50m (Figure 3). The assessment of the best model was based on the values of R2 and SEE (Figure 3) and it indicated that the satisfactory results were obtained for h = 900m because SEE = 0.232, where it actually assumes the lowest value and it corresponds to the highest value of R2 = 0.74 (Table 3). Moreover, this depth below surface corresponds to the average depth of exploitations in the mine.
Results of multiple regression with statistical assessment and calculated AGMPE
|Site||Logarithmized amplification factor ai||Relative amplification AGMPE||α||β||γ||R2||SEE|
The calculated regression model parameters and descriptive statistics were compared in Table 3. A similar study was carried out for the “Ziemowit” mine in 2012 . Comparing the regression parameter values with those calculated in the paper by Golik and Mendecki  with those obtained here, only γ remains the same. It can be also treated as another qualitative confirmation of correctness of the model as it indicates that the coefficient of geometric scattering (associated with geology) has been preserved. The β value in this study is higher than the previous one . This parameter is related to the b-value of Gutenberg-Richter , thus the increase of the data set can change or refine the b-value. This can explain changes in β parameter. The data set presented here is larger than the one used in . The same reason might apply to the differences in the relative amplification coefficients. In this study, the number of seismic data and number of sites increase as well.
Next, the map of spatial distribution of the relative amplification coefficients AGMPE was prepared (Figure 4A) using the block kriging as the gridding method with a linear drift function [14,15]. The map was prepared in the Surfer 10 software. The semi-variogram model applied to the data was a power function. The grid was the size of 100 rows and 73 columns. The descriptive statistics for the analyzed data set of 14 relative amplification coefficients are characterized by the mean value of 2.1, standard deviation of 0.4, maximum value of 3.6 and minimum value of 1.0.
The calculated GMPE amplification values range from 1 to 3.7. The strongest site effects (AGMPE > 2) can be reported for the station areas in: Dzieckowice, Imielin, Chełm, MSK SUW and Zapora. At the last station (Zapora), AGMPE was the highest and reached 3.7 (Table 3). The GMPE model was also assessed by means of further statistical techniques which confirmed correctness and good fit of the model. The analysis of residuals indicated that residuals are focused around zero (Figure 4B) and the residual histogram has the Gaussian-shape distribution [2, 5]. The reliability of the model is also confirmed by the normality plot (Figure 4C) which is almost a straight line with minor discrepancies at the ends. In Figure 4D, the dispersion plot was placed. It represents the relation of the predicted PGA versus observed PGA. Its cigar-shape is also the confirmation that the model reproduces estimated (predicted) data well. According to Lasocki , the visual inspection of residuals characteristics (Figure 5) also allows to assume that the residuals do not depend on the station location (Figure 5A), because the mean residuals have very small values (near zero) and 95% confidence intervals associated with the station locations indicated the value zero falls well within the respective confidence intervals. Moreover, there is no distinct dependence of residuals on the epicentral distance (Figure 5B).
The statistical analysis of the proposed model (Figure 4) and inspection of the residuals (Figure 5) convince that the calculations made are appropriate and the GMPE relation can be used for seismic hazard study. Figure 6 shows the last comparisons between the observed amplitudes, predicted values that are the medians [1, 3] as well as 95% confidence intervals for prediction that were calculated according to the formula presented in [1, 3]. The values were sorted in a descending order of estimates and the amplitude values were unlogarithmized. Most of the observed amplitudes fall within the confidence intervals, and only 2.7% of cases exceed the upper limit (417 cases) as well as 2.2% of cases are below the lower limit (345 cases).
Thicknesses of the loose layers, in which the amplification may occurred, vary on the entire study area. Comparison of the calculated values of AGMPE with the thicknesses of Quaternary deposit or combination of Quaternary and Tertiary loose sediments (Table 4) indicated that there is no visible correlation. The correlation coefficient for Quaternary thickness and AGMPE reached value of 0.026. Here, the thickness varies in small range from 9 to 31 m and this layer is mainly composed of sand with tilt inserts, thus it can be assumed that there are similar geological conditions.
Comparison of geology settings with the AGMPE values
|No.||Seismometer station||Thickness [m] (main type of rock)||AGMPE|
|1||Fawent||22.0 (sands1)||4.0 (silts)||65.0 (limest.)||1.59|
|2||Goławiec||13.4 (sands2)||142.6 (silts)||12.7 (marls)||1.96|
|3||Chełm||22.0 (sands1 )||4.0 (silts)||55.5 (limest.)||2.35|
|4||Imielin||30.0 (till+)||210.0 (siltst.3)||20.0 (congl.)||2.31|
|5||Szyb W-II||28.5 (sands+)||232.5 (siltst.3)||-||1.71|
|7||Pompownia||12 .0 (sands+)||171.9 (marls4)||-||1.03|
|8||Zapora||9.0 (sands+)||9.0 (silts)||11.0 (limest.)||3.70|
|10||MSK||12.0 (sands+)||171.9 (marls4)||-||2.52|
|11||Kopciowice||39.0 (sands)||19.0 (till)||53.5 (limest.)||1.84|
|12||Czerniny||16.0 (sands)||34.5 (silts)||48.0 (limest.)||1.00|
|13||SUW||31.0 (sands+)||89.7 (silts)||19.7 (dolom.)||2.19|
|14||Dziećkowice||30.5 (sands+)||89.7 (silts)||19.8 (dolom.)||3.10|
Due to the lack of correlation for Quaternary deposits we assumed cumulative thickness of all loose sediments (Quaternary and Tertiary) and the correlation coefficient was still low (r = 0.153), what also indicated that there is a lack of visible relation between amplification coefficient and thickness. Therefore, it seems that this relation might be more complex and considers other petrophysical parameters such as wave velocities, densities or water content that should be taken into account. However, this parameters remain still unknown and require next tests. More common approach is to compare resonance frequency of loose layer with thickness and some authors [16,17,18] reported that there is a clear power-law relation between these two parameters, but in this case the GMPE relation did not provide such data. The best solution will be to carry out HVSR measurements introduced by Nakamura .
Additional problem is the boundary between amplifying layer and rigid substratum. The geology is complex thus there is difficulty to find the appropriate thickness of the loose layer that can be excited by seismic waves. This aspect could also have an influence on the lack of correlation.
The statistical evaluation of the proposed GMPE model allows to state that this model can be used in the future PGA prediction in the central part of the Main Syncline of the USCB. The reliability of the model was confirmed by:
- analysis of residuals that allows to check the quality of models;
- coefficient of determination that informs what percentage of the original variability of the dependent variable was explained by the model;
- standard error of estimation;
- visual inspection of residual characteristics that showed that the mean residuals possess very small values (near zero) and 95% confidence intervals associated with the station locations indicated the value of zero and the residuals did not depend on epicentral distances;
- and, finally, less than 5% of data exceeded the 95% confidence intervals.
Comparing the obtained relative amplification coefficients with the geological settings, it can be concluded that the highest amplification values and, thus, the seismic hazard occur in the outskirts of the Syncline (N, NE and SE), i.e. on the opposite sides of the basin. The lowest amplification values, therefor the lowest seismic hazard areas, are located in the central part of the studied area, in the zone of a block formed by faults: Lędziński, Imlieliski, Eastern and Chełmski which are in the main axis of the Syncline in the W-E direction.
The authors would like to thank the director and the employees of the Piast-Ziemowit Coal Mines for all support and consent to publish this work. The authors would like to highlight the significance of the reviewers’ comments that improved this paper and would like to thank them for their revisions.
Olszewska D., Attenuation relations of ground motion acceleration response spectra for the Polkowice Region. Publs. Inst. Geophys. Pol. Acad. Sci, 2006, 1000(395), 161-174.
Golik A., Mendecki M., Ground-motion prediction equations for induced seismicity in the main anticline and main syncline, Upper Silesian Coal Basin, Poland. Acta Geophysica, 2012, 60(2), 410-425.
Lasocki S., Site specific prediction equations for peak acceleration of ground motion due to earthquakes induced by underground mining in Legnica-Głogów Copper District in Poland. Acta Geophysica, 2013, 61(5), 1130-1155.
Chodacki J., New Ground Motion Prediction Equation for Peak Ground Velocity and Duration of Ground Motion for Mining Tremors in Upper Silesia. Acta Geophysica, 2014, 64(6), 2449-2470.
Lasocki S., Relacja tłumienia wartości szczytowej składowej poziomej przyspieszenia drgań gruntu w paśmie częstotliwości do 10Hz dla rejonu miasta Polkowice. Publs. Inst. Geophys. Pol. Acad Sci., 2002, M-27 (352), 79 – 90 (in Polish).
Campbell K.W., Strong-motion attenuation relations. In: Lee W.H.K., Kanamori H., Jennings P.C., Kisslinger C. (Eds.), International Handbook of Earthquake and Engineering Seismology, Part B, IASPEI, 2003, 1003-1012.
Frej A., Zuberek W., Dyskusja nad amplifikacją drgań sejsmicznych w rejonie niecki bytomskiej. WUG: Bezpieczeństwo Pracy i Ochrona Środowiska w Górnictwie, 2007.6, 35-37 (in Polish).
Frej A., Zuberek W. M., Local effects in peak accelerations caused by mining tremors in Bytom Syncline Region (Upper Silesia). Acta Geodyn. Geomater, 2008, 5(2), 150.
Bonnefoy-Claudet S., Cotton F., Bard P., The nature of noise wavefield and its applications for site effects studies: A literature review. Earth-Science Reviews, 2006, 79, 205–227.
Stec K., Characteristics of seismic activity of the Upper Silesian Coal Basin in Poland. Geophysical Journal International, 2007, 168(2), 757-768.
Kędzior S., Methane contents and coal-rank variability in the Upper Silesian Coal Basin, Poland. International Journal of Coal Geology, 2015, 139, 152-164.
Documentation of Ziemowit Coal Mine in Lędziny
Boore D.M., Atkinson G.M., Ground-motion prediction equations for the average horizontal component of PGA, PGV, and 5%-damped PSA at spectral periods between 0.01 s and 10.0 s. Earthquake Spectra, 2008, 24(1), 99-138.
Hudson G., Wackernagel H., Mapping temperature using kriging with external drift: theory and an example from Scotland. International journal of Climatology, 1994, 14(1), 77-91.
Ciotoli G., Stigliano F., Marconi F., Moscatelli M., Mancini M., Cavinato G.P., Mapping the anthropic backfill of the historical center of Rome (Italy) by using intrinsic random functions of order k (IRF-K). In: Murgante B., Gervasi O., Iglesias A., Taniar D., Apduhan B.O. (Eds.) International Conference on Computational Science and Its Applications, Springer, Berlin, Heidelberg, 2011, 92-102.
Delgado J., Lopez Casado C., Giner J., Estevez A., Cuenca A, Molina S., Microtremors as a Geophysical Exploration Tool: Applications and Limitations. Pure and Applied Geophysics, 2000, 157, 1445-1462.
Dinesh B.V., Nair G.J. Prasad A.G.V., Nakkeeran P.V., Radhakrishna M.C., Relation between sedimentary layer thickness and fundamental frequency of the H/V spectra for Bangalore city. Current Science, 2009, 97, 7, 1073-1077.
Mendecki M.J., Bieta B., Mycka M., Determination of the resonance frequency–thickness relation based on the ambient seismic noise records from Upper Silesia Coal Basin. Contemporary Trends in Geoscience, 2014, 3(1), 41-51.
Nakamura Y., A method for dynamic characteristics estimation of subsurface using microtremor on the ground surface. QR Railway Tech. Res. Inst., 1989, 30(1), 25-33.
clays inserts, limest