High-resolution shallow crustal shear wave velocity structure of Anyuan mining area and its adjacent region in Jiangxi Province, China

High-resolution seismic image is critically important for mining minerals. In this study, we collected seismic data from a local dense seismic array consisting of 154 stations around the Anyuan mining area and its adjacent region of Pingxiang City, Jiangxi Province in South China, and applied the ambient noise tomography (ANT) method to image the shear wave velocity structure in the study area. Shallow crustal velocities at depths less than 3.3 km were determined by direct inversion of Rayleigh wave group velocity dispersion curves at the period range of 0.5–5.0 s. Overall, the S-wave velocity structure has a tight correlation with surface geological and tectonic features in the study area. The shear wave velocity structure in the shallow crust of the Anyuan Mine and its adjacent areas displayed distinct low-velocity anomalies, which can be attributed to the depression of sedimentary structures and coal mining activities in the Pingxiang-Leping region. The zones surrounding the Anyuan fault (AYF) and Wangkeng fault (WKF) zones exhibited low-velocity anomalies from the ground surface to ~ 3.3 km underground. And the low-velocity anomalies at depths less than 1.2 km could be related to the sedimentary environment of coal mine and the coal mining activities, while the low-velocity anomalies at depths below 1.2 km are caused by the presence of fracture medium, oil and gas in the fault zone. The shear wave velocity changes sharply across the AYF, and the characteristics of the velocity change interface indicate that the AYF is inclined toward the northwest, with its extension reaching depths of approximately 3 km underground.


Introduction
In recent decades, the ambient noise tomography (ANT) method develops quickly to image the crustal velocity structure.The basic idea of this method is to extract Empirical Green Functions (EGFs) from cross-correlations of ambient noise and coda waves between station pairs with sufficiently long records (Campillo and Paul 2003;Shapiro and Campillo 2004;Sabra et al. 2005).Based on the time reverse theory, each station can be set as a virtual source or a signal-receiving station (Cassereau et al. 1992).Therefore, this method is no longer limited by the temporal and spatial distribution of seismic events, and it is very suitable for deciphering precise crust and mantle velocity structures in low seismicity regions (Li et al. 2016a;Zhang et al. 2016;Liu et al. 2018).Benefiting from this advantage, the ANT has been widely used to image crustal and upper-mantle velocity structures at various scales (Sabra et al. 2005;Yao et al. 2006Yao et al. , 2008;;Yang et al. 2007;Lin et al. 2007Lin et al. , 2008;;Bensen et al. 2007Bensen et al. , 2008;;Zheng et al. 2011;Li et al. 2009Li et al. , 2010Li et al. , 2012)).
In recent years, with the development of dense seismic array technology, the high-resolution structure of the shallow crust inversion at a small regional scale using high-frequency surface wave in ambient noise has developed significantly (Huang et al. 2010;Lin et al. 2013;Fang et al. 2015;Li et al. 2016b).Among these inversion methods, Fang et al. (2015) developed a direct surface wave imaging technology based on the typical ANT method.This technology can directly image the 3-D shear wave velocity from the surface wave dispersion between station pairs, omitting the intermediate steps of surface wave group velocity and phase velocity imaging.Moreover, the direct surface-wave tomography method takes influences of surface wave ray bending in complex medium, thus can provide a more accurate velocity structure of the shallow crust at small regional scale (Fang et al. 2015;Li et al. 2016b), and the advantages of this method have been validated by numerous tomographic studies (Wang et al. 2018;Gu et al. 2019;Xie et al. 2021;Huang et al. 2021).For instance, Wang et al. (2018) used this method imaged the shallow 3-D S-wave velocity structure of the Hutubi gas storage area in the southern margin of Junggar Basin and revealed the characteristics of near-surface structure and the gas storage.Gu et al. (2019)  The shallow crustal structure provides fundamental clues for analyzing metallogenic mechanisms and estimating the ore reserve.The Anyuan mining area located in the western section of the Pingxiang-Leping depression in Jiangxi Province, South China (Fig. 1), which is a key coal production base of Jiangxi Province with a mining history of over a century.The Pingxiang-Leping depression is a large depression oriented in the NEE direction, bounded by deep faults on the north and south sides, and the late Paleozoic and Mesozoic strata are widely exposed.The depression exhibits a variety of facies, including marine, continental, and transitional, from the Devonian to the Paleogene.Several sets of rock source strata, primarily composed of black shale, have formed in the depression, such as Middle Permian (Xiaojiangbian Formations-Maokou Formation Nangang section), Upper Permian (Leping Formation), and Upper Triassic (Anyuan Formation) (Li et al. 2003;Deng et al. 2005).Our study region belongs to Anyuan Formation, including three mining areas, named Baiyuan, Anyuan and Gaokeng minings (Fig. 2b).Among these, only the Anyuan mining area is currently operational (Xiao 2020).There are three NE-trending reverse faults distributed in the study region, named Luojiachong fault (LJCF), Anyuan fault (AYF), and Wangkeng fault (WKF) (Tan et al. 2001).
Currently, there are no dangerous rock, collapse, landslide, and other adverse geological phenomena in the mining area (Liu et al. 2021).However, due to the problems of resource exhaustion, gas leakage, ground fissures, ground collapse, and mine earthquake in century-old mine, it is particularly urgent to understand the 3-D structure of the shallow crust of the Anyuan Mining area and its adjacent region.Although Li et al. (2015) have studied the 3-D crustal velocity structure of Jiangxi Province, the detailed velocity structure of the shallow crust around the Anyuan mining area needs further investigation.
Therefore, it is essential to develop a high-resolution image of the shallow structure in the Anyuan mine region is of critical importance.In this study, we apply the direct surface-wave tomography method (Fang et al. 2015) to image the 3-D shear wave velocity structure of the shallow crust in Anyuan mining area and its adjacent region based on a local dense array, and our results could provide a significant scientific basis for coal mine resources exploring and seismic structure modeling.

Data collection and processing procedure
The vertical component continuous data from 154 seismic stations around the Anyuan mining area (Fig. 2a) deployed by East China University of Technology and the Jiangxi earthquake Agency between 27 July 2020 and 31 August 2020 are used in this study.This temporary dense seismic array consists of 13 broadband seismographs of  (Zhang et al. 2013).The red rectangle in the bottom map corresponding the study area as shown in Fig. 2a and the black lines denote the faults GL-PCS60 type (the blue triangles in Fig. 2a), 98 shortperiod seismographs of EPS type (the red triangles in Fig. 2a), and 43 short-period seismographs of QS-05A type (the yellow triangles in Fig. 2a).The seismic array spans an area of about 18 km × 18 km with an average station spacing of about 1 km in the center area and the minimum station spacing is ~ 0.5 km.The collected data sets are recorded with a simple rate of 100 Hz.
The data processing procedure follows the methodology proposed by Bensen et al. (2007).Seismograms are first cut into daily segments.After the removal of the mean, trend, and instrument response, a 0.5-10 s band-pass filter is applied.To reduce the effect of earthquakes and instrumental irregularities on crosscorrelations, the seismograms are normalized using a time-frequency normalization method (Ekström et al. 2009;Shen et al. 2012).Then daily cross-correlations for each station pair are calculated and these are stacked using the normalized linear stacking method to obtain cross-correlation functions (CCFs).Finally, the empirical green functions (EGFs) are derived from the time derivative of the CCFs (Sabra et al. 2005;Yao et al. 2006).
Theoretically, surface wave signals should appear symmetrically on both the causal and acausal time lags of the EGFs.However, the uneven and varied spatial distribution of noise sources can affect the EGFs and make the signal emerge on the positive and negative time lags asymmetric (Stehly et al. 2007;Zhan and Clayton 2010).Therefore, to enhance the signal-to-noise ratio (SNR), the EGFs are reversed and stacked to create a symmetric component before the dispersion measurement extraction (Yang et al. 2007;Yao et al. 2011).Figure 3a illustrates the path distribution of the station pairs between the L124 station and other stations.Figure 3b displays the CCFs of station pairs with 0.5-5 s band-pass filtered and aligned by the station-pair distance.The fundamental order signals of the Rayleigh wave were clearly observed, and the surface wave group velocity ranges from approximately 1.7 to 3.6 km/s (Fig. 3b).

Surface wave dispersion curve extraction
In general, when using the ANT method to image the surface wave velocity structure over a long time period, the inter-station distance of at least three wavelengths is required to obtain reliable frequency dispersion curves (Yao et al. 2006;Zheng et al. 2008).Yao et al. (2011) and Luo et al. (2015) demonstrated that dense seismic arrays can yield accurate dispersion curves in shortperiod surface wave inversion, even when the interstation distance is only one times the wavelengths.So, we use multiple filtering techniques (Dziewonski et al. 1969;Yao et al. 2011) to extract the high-quality dispersion data from the EGFs with inter-station distances greater than one wavelength and SNR > 5.The SNR was defined as the ratio of the maximum amplitude of the signal window of EGFs (− 15 to ~ 15 s) and the average absolute amplitude of the noise window (− 30 to ~ − 15 s and 15 to ~ 30 s).
Furthermore, we removed the measurements that exceed the average value of group velocity at each period plus or minus 0.5 km/s and eliminated the data points with too high or negative slopes.Ultimately, we retained 4538 high-quality group velocity dispersion curves out of the 11,781 theoretically available dispersion curves from 154 seismic stations (Fig. 4). Figure 4 illustrates the variation of surface wave group velocity Fig. 3 Cross-correlation functions between station L124 and other stations.a Ray-path coverage map between station L124 and other stations, respectively.b Shows cross-correlation functions with 0.5-5 s band filter between station L124 and other stations, the dash inclined lines denote the ranges of velocity (1.7 to 3.6 km/s) Fig. 4 Group velocity dispersion curves extract from the EGFs.The black thin solid lines denote the selected group velocity dispersion curves and the gray thin solid lines denote the removed dispersion curves.The black thick solid line shows the number of the dispersion at each period and the thick gray dashed line presents the average group velocity dispersion curve ranging from approximately 1.8 to 3 km/s, with the majority of periods distributed between 0.6 and 3.5 s.

Shear wave velocity structure inversion
We use the direct surface wave tomography method (Fang et al. 2015) to directly invert the mixed path Rayleigh wave group velocity dispersion data into a 3-D velocity structure.This method first assumes an initial 3-D velocity model and then performs surface wave ray tracing at each period using a fast marching method (Rawlinson and Sambridge 2004) and iteratively update the sensitivity kernels of period-dependent dispersion measurements to updates of the 3-D velocity model.To mitigate the impact of the uneven path coverage, a wavelet-based sparsity-constrained seismic tomography method (Simons et al. 2011;Fang and Zhang 2014) is employed in the 3-D wave velocity inversion.
As is well-known, the essence of inversion is to find a model m that minimizes the differences δt i (ω) between the observed arrival times t obs i (ω) and theoretical arrival time t i (ω) for all frequencies ω .The travel-time pertur- bation for the path i at frequency ω is expressed by Fang et al. (2015) as Eq. ( 1): where t i (ω) denote the calculated travel-time using a ref- erence model that can be updated in the inversion, and (1) , C k (ω ) and δC k (ω) present the phase velocity and its per- turbation of the kth grid node at the angular frequency ω , respectively.v ik is the bilinear interpolation coefficients along the ray-path associated with the ith travel-time data.
Generally speaking, surface wave dispersion is most sensitive to shear wave velocity ( β ), but short-period Ray- leigh wave dispersion is also very sensitive to compressional wave velocity ( α ) in shallow crust.Therefore, when we use the empirical formula (Brocher 2005) to relate the compressional wave velocity and mass density ( ρ ) to the shear wave velocity, Eq. ( 1) can be rewritten to Eq. (2) as Fang et al. (2015) described: where k present the 1-D initial model at the kth surface grid point on the surface and α k (z j ) , β k (z j ) , and ρ k (z j ) are the compression wave velocity, shear wave velocity, and mass density at the jth depth grid node, respectively.R α and R β present the scaling factor, and J is the number of grid points in the depth direction.The total grid points (2) where the first term on the right-hand side gives the l 2 -norm data misfit and the second term denotes the l 2 - norm model regularization term.L is the model smoothing factor and λ is the weight factor of data fitting and model regularization (Simons et al. 2002;Aster et al. 2013).Thus, the solution m can be deduced from following equation (Paige andSaunders 1982a, 1982b): (3 (5 To obtain a robust initial model, two methods are used to determine the 1-D mean velocity model.The first one is based on empirical formulas (Shearer et al. 2009;Fang et al. 2015), which involved using 1.1 times the average measured group velocities at a depth of 1/3 times the wavelength as the velocity model (Fig. 4a, hereafter, called AVG model).Subsequently, the AVG model is used as the input model and the iteratively damped least square method developed by Herrmann and Ammon (2002) in the CPS330 package (Herrmann 2013) is applied to obtain the other 1-D velocity model (Fig. 4b, hereafter, called CPS model).After that, the AVG and CPS models are used as the initial models for the direct 3-D surface wave inversion program to calculate the 1-D mean shear wave velocity model, respectively (Fig. 5).
Despite the significant differences between the two input models (AVG and GPS models), the inversion results were found to be essentially the same (Fig. 5a, b).This suggests that the two velocity models have minimal impact on the inversion results in this study.Previous studies by Yu et al. (2020) and Xie et al. (2021) have also proved that short-period surface wave inversion has relatively low dependence on the initial model.Thus, in view of the smaller difference between the CPS model and the inverse model, the CPS velocity model is selected as the initial model for the direct 3-D surface wave velocity structure inversion.Figure 5c shows the fitting of the dispersion curve (solid black line) and the original dispersion curve (black circle).
In the direct inversion method, the weighting factor λ in Eq. ( 4) is particularly important in the residual term of the data measurement and the regularization term of the model (Paige andSaunders 1982a, 1982b).If the λ is too small, the result may be over-fitted and the error will increase.Conversely, if the λ is too large, the inversion model will be excessively smooth, making the inversion result highly reliant on the initial model (Luo et al. 2019).Hence, the L-curve analysis method (Hansen 2007) is used to determine the optimal λ, and the result reveals the optimal weighting factor λ is 3.
During the inversion process, the depth sensitivity kernel is used to assess the sensitivity of surface wave group velocity on shear wave velocity in the depth direction (Yao et al. 2008;Zheng et al. 2008).Figure 6 displays the depth sensitivity kernel functions at periods of 0.6, 1.0, 2.0, 3.0, 4.0, and 4.5 s for the fundamental Rayleigh wave group velocity calculated from the CPS model as shown in Fig. 5b.It can be observed that the surface wave dispersion data between 0.6 and 4.5 s have the capability to resolve depths down to ~ 5 km underground.

Results and discussion
To enhance the robustness of the 3-D shear wave inversion results, the checkerboard resolution tests are performed to assess the lateral and vertical resolutions.Theoretically, the recovery ability of the checkerboards depends primarily on the coverage of surface-wave raypath.Thus, we first check the ray-path coverage at different periods and then performed the checkerboard resolution test.Figure 7 illustrates the ray-path coverage of group velocity measurements at six selected periods, obtained from the final velocity model.The results indicate that the study region has a relatively high density of surface-wave ray-path coverage, suggesting that the structure of the region can be well imaged by the data set.We also observe that the ray-path number decreases with increasing period, such as the number of the ray-path at 4 s period is only about 27% of that at 1.5 s period.
To assess both the lateral and vertical resolution, we use a 3-D checkerboard model.The 3-D reference velocity model is constructed based on the CPS model as shown in Fig. 5b.After several comparative tests, the study region is meshed with 25 by 25 grid points with an interval of 0.01° in latitude and longitude.And 18 grid nodes are set along the depth direction from the ground surface to 13 km underground, according to the interface depth of the velocity structure.For the lateral resolution test, the anomaly size is approximately 0.025° × 0.025°, while for the vertical resolution tests, the grid size is 0.025° along the horizontal direction, gradually increasing along the vertical direction with increasing crustal depth, ranging from approximately 0.8 km to 1.5 km. Figure 8a-d shows the recovery results of the lateral checkerboard at depths of 0.4, 1.2, 2.2 and 3.3 km.The surface locations and the inverted results of vertical profiles MM' and NN' are shown in Fig. 8d and Fig. 8f, h.The recovery results indicate that the recovery of the detection board is influenced by the density of ray-path coverage.The checkerboard test is well recovered in the area with dense ray-path coverage, while the recovery is poor in the area with poor ray-path coverage.Overall, the results suggest that the structures of the Anyuan mine and its adjacent During the inversion process, we use the same velocity model and the parameters as those in the checkerboard resolution test.The inversion process is relatively stable and the standard deviation of residual error (STD) of surface-wave travel-time decreasing as the number of iterations times increased (Fig. 9) The STD decreases rapidly during the first three iterations, followed by a slower decrease until convergence.After 10 times iterations, the average travel-time residuals decreased from 0.454 s to 0.354 s, indicating the final velocity model fits the observed data better.
The coal-bearing strata in the study region belong to the Anyuan Formation of the Upper Triassic, which can be divided into Zijiachong member, Sanjiachong member, and Sanqiutian member from bottom to top (Fig. 10).The lithofacies assemblages of the Anyuan Formation consist mainly gray to light gray, light purple sandstone, siltstone, sandy shale, and carbonaceous shale (Li et al. 2016c).And the sedimentary environment is primarily marine-continental transitional facies and the coal seams exhibit multiple layers, complex structures, and uneven distribution (Liu and Chen 2021).
Figure 11 shows the shear wave velocity structure at different depths obtained from group velocity dispersion data through the direct surface inversion.The patterns reveal significant variations in shear wave velocity structure at the same depth, indicating strong lateral heterogeneity in the shallow crustal structure of the study region.At the depth of 0.4 km underground, the basin and coal mining areas show low-velocity anomalies, while the uplifted mountain areas without coal mine distribution exhibit high-velocity anomalies.The lowvelocity anomalies in the basin areas can be attributed to the sedimentary layer on the ground surface, while the high-velocity anomalies in the uplifted mountain areas indicate missing sedimentary layers or exposed bedrock.However, the Anyuan mine area located in the uplifted mountainous, exhibits low-velocity anomalies in its shear wave velocity structure.This abnormal phenomenon could be related to the sedimentary environment of coal seam and coal mining activities.The Anyuan mining area and its adjacent regions have a sedimentary environment characterized by the interaction of marine and continental facies with sedimentary basins, containing significant amounts of sand and silt (Xi et al. 1999;Li et al. 2003).Additionally, the Anyuan coal mine has a mining history of over one hundred years, resulting in the production of a large amount of coal mine gob area and coalbed methane during the mining process (Li et al. 2003;Qiu and Li 2019).All these factors contribute to the presence of lowvelocity anomalies in the shear wave velocity structure of the Anyuan mine area.
The distribution of shear wave velocity structure in the depth ranges of 1.2-2.2km exhibits similar characteristics.The Anyuan and Gaokeng mining areas still show low-velocity anomalies, but the ranges of the anomalies decrease with the increase of depths, gradually approaching the WKF and AYF zones.It is also observed that the shear wave velocity around the AYF zone and the piedmont basin forms a distinct "Y-shaped" zone of low-velocity anomalies.Overall, in the depth ranges of 0.4-1.2km underground, the shear wave velocity in the entire Anyuan mining area exhibits significant low-velocity anomalies.However, at a depth of 2.2 km underground, the shear wave velocity structure of the Anyuan mining area is bounded by the WKF zone.The shear wave velocity structure in the northern part of the mining area still presents low-velocity anomalies, while the southern part gradually transitions to high-velocity anomalies, consistent with the distribution characteristics of shear wave velocity structure in high mountain areas.
At the depth of 3.3 km underground, the shear wave velocity structure is well consistent with the surface and topographic features.Mountainous areas with high terrain, such as the areas between the LJCF and AYF belts and the southern part of the WKF belt exhibit highvelocity anomalies, while basins and piedmont basins exhibit low-velocity anomalies.It is worth noting that despite the high elevation of the AYF and WKF zones, there are still low-velocity anomalies in the shear wave velocity structure around the fault zones.Previous studies on the Pingxiang-Leping depression in Jiangxi Province have shown that there are abundant Coalbed methane resources beneath the Anyuan mining area and its adjacent region (Li et al. 2003;Qiu and Li 2019).Therefore, the low-velocity anomalies observed at depths below 2 km around the AYF and WKF regions may indicate the presence of oil and gas between the fault zones.
Figure 12 shows the positions of the four sections and their corresponding shear wave velocity profiles.The AA' profile shows that the basin areas to the west of the AYF exhibit low-velocity anomalies, and the depths of these anomalies are negatively correlated with the surface elevation.Along the profiles AA' and CC' , there is a sharp change in the shear wave velocity structure across the AYF zone, with the velocity transition interface inclined towards the northwest, consistent with the trend of the AYF (Tan et al. 2001;Liu et al. 2004).The velocity transition interface reaches depths of approximately 3 km underground, which may indicate the extension depth of the AYF.The BB' and CC' profiles show a noticeable low-velocity anomaly on the northwest side of the AYF, possibly caused by the presence of sediment and debris in the piedmont floodplain (Tan et al. 2001;Li et al. 2016c), which significantly reduces the propagation speed of shear waves.
Along the AA' , BB' , CC' and DD' profiles, it can be seen that the shear wave velocity structure in the shallow surface layer (< 1.2 km) presents low-velocity anomalies, which could be related to the sedimentary environment of the coal seam.Xi et al. (1999) and Li et al. (2003) have suggested that the Anyuan mining area and its adjacent region are primarily composed of sea-land interaction deposition, with a maximum deposition thickness of up to 1.3 km.We observed that the low-velocity anomalies in the shallow surface layer (< 1.2 km) in the Baiyuan, Anyuan, and Gaokeng mine areas are particularly significant, which may be caused by coal mining activities.Because the presence of coalmine gob area and coalbed methane production during coal mining can slow down the shear wave velocity.Additionally, the low-velocity anomalies observed around the AYF and WKF zones extend to depths of approximately 3 km (Figs.11 and 12), which could be associated with fracture media and the presence of oil and gas in the fault zone.The top coal mining area is often composed of swamp facies or tidal flat facies mudstone and silty mudstone, which has a strong sealing ability for coal reservoirs (Qiu and Li 2019).Consequently, the oil and gas produced during coal mining could migrate into the faults and gaps between the fault zones, resulting in low-velocity anomalies in the shear wave.

Conclusion
We collected about 1 month of seismic data recorded by a local dense array and applied direct surface-wave tomography to image the high-resolution shallow crustal shear wave structure around the Anyuan mining area and its adjacent region of Pingxiang City, Jiangxi Province in South China.The results show that the basins and piedmont basins areas presented low-velocity anomalies, while the uplifted mountainous areas without coal mining exhibited high-velocity anomalies.Influenced by the sedimentary structural depressions and coal mining activities, the mining areas exhibited significant low-velocity anomalies in the shallow surface (< 1.2 km).Specifically, the AYF and WKF zones exhibited low-velocity anomalies from the ground surface to a depth of ~ 3 km underground.The low-velocity anomalies at depths less than 1.2 km could be related to the sedimentary environment of coal mines and the coal mining activities, while the low-velocity anomalies at depths below 1.2 km were caused by the presence of fracture medium, oil and gas in the faults zone.Furthermore, the result also revealed the AYF inclined towards the northwest and extended to a depth of ~ 3 km underground.In summary, our results agreed well with the topography and regional geological structure of the study areas, which could reveal information on the coal seam distribution, coal mining activity, and fault geometry to a certain extent.
imaged the detailed shallow subsurface structure around the Tanlu fault zone by deploying 53 stations with an average spacing close to 5 km.Xie et al. (2021) imaged a 3-D crustal velocity structure within 5.0 km underground the mining area and provided an effective basis for exploring deep-buried gold deposits in the Qingchengzi eastern of Liaoning Province.Based on this method, Huang et al. (2021) imaged a 3-D S-wave velocity structure within 1 km underground the well field area and predicted the distribution characteristics of gas content by combining petrologic experiment in the Sijiazhuang Coal mine in Shanxi Provinces of China.

Fig. 1
Fig. 1 Tectonic map of South China block and the Pingxiang-Leping depression areas, the bottom map is enlarged of the black rectangle area in the top map.The thick black dash line presents the boundary of the South China block and the thin black solid line is the boundary between the Yangtze Block and the Cathaysia Block in the top map (Zhang et al. 2013).The red rectangle in the bottom map corresponding the study area as shown in Fig. 2a and the black lines denote the faults

Fig. 2
Fig.2The distribution of seismic stations and the geological map of the study area.a The three different types of seismic instrument (GL-PCS60, EPS, and QS-05A) locations are shown as the blue, red, and yellow triangles, respectively.White closed lines represent the distributions of Anyuan, Baiyuan, and Gaokeng mining exploitation areas (advised from Xiao 2020).The main region of the inversion result is given by the blue dotted square.The black lines denote the three main faults in the study region (according to the 1:200,000 Pingxiang fault of Jiangxi Province).b The geological background of the Anyuan mining area and its adjacent region.The brown rectangle denotes major coal deposits (modified according to the 1:1,000,000 East China geological map and 1:500,000 Jiangxi Province geological map from http:// www.ngac.cn/).(For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

Fig. 5
Fig. 5 Comparison of inversion results of different initial models.a The green line represents the AVG model deduced from the average dispersion curves.The gray lines and the red line represent the intermediate iteration and the final model of the inversion.b The blue dashed line represents the CPS model computed from the CPS330 inversion.The gray and red lines represent the intermediate iteration and the final model of the inversion, respectively.c The fitting situation of the observed dispersion curve (green line) and the fitting dispersion curve (blue dashed line)

Fig. 6
Fig.6Depth sensitivity kernel functions at periods of 0.6, 1.0, 2.0, 3.0, 4.0, and 4.5 s for the fundamental Rayleigh wave group velocity calculated from the CPS model in Fig.4b

Fig. 7
Fig. 7 Ray-path coverage maps of Rayleigh wave group velocity measurement from the final 3-D model at different periods

Fig. 9
Fig. 9 Variation of the standard deviation of residual errors of surface-wave travel-time with 10 iterations of the CPS model.The insert image on the upper right corner is the distribution of the residuals

Fig. 11
Fig. 11 Shear wave velocity slices at different depths: a 0.4 km, b 1.2 km, d 2.2 km, f 3.3 km.Black lines are the three main faults in the study area.White closed lines represent the distribution of mining exploitation areas

Fig. 12
Fig. 12 Selected vertical profiles locations and their shear wave velocity distributions: a AA′, b B-B′, c C-C′ and d D-D′