Abstract
Crustal response induced by the subduction of Paleo-Pacific Plate, is crucial to study the metallogenic law and prediction of Jiaodong type gold deposit. Gold deposits are preferentially hosted in the lithospheric-scale second-order faults, which provide migration channels, mineralization corridors, and the physicochemical conditions. Super-large gold deposits in northwestern Jiaobei Terrane are generally controlled by regional linear faults, such as Sansandao, Jiaojia, Zhaoping and Qixia faults. However, Zhaoping fault is exposed sporadically (in Xiadian and Jiangjiayao gold deposits) beneath a thick cover of Quaternary basalts, which indicating that limited information can be obtained by traditional geological survey. In this study, high-precision geophysical surveys such as gravity and controlled source audio frequency magnetotellurics (CSAMT) are conducted in the Shanhou gold deposit, which hosted in southern part of Zhaoping fault. Three integrated geophysical profiles results further identified the spatial distribution of second-order oblique-slip faults of Zhaoping fault and established the relationships between brittle deformation and mineralization. The spatial variability of stratum (high density 2.82 g/cm3 and intermediate to high specific resistivity 4000 Ω m), brittle-slip faults (intermediate density 2.73 g/cm3 and low specific resistivity 200 Ω m) and magmatite (low density 2.57 g/cm3 high specific resistivity 8000 Ω m) are identified by two inversion techniques, including the man-machine interactive inversion technology and nonlinear conjugate gradient technology. In addition, the geophysical inversion pseudo sections delicately portray a smooth-out waveform low specific resistivity anomaly, which strongly couples with the spatial ore-hosted locations and enrichment places. The gold mineralized alteration zones generally dips to southeast and the angle changing from 45° to 30°, which are strong coupling with the large-scale acid magmatic intrusion and extensional tectonics occurred in the late Mesozoic. Finally, according to the deep prospecting prediction, a gold ore body with 10 m in thickness was revealed by drilling below the depth of 1,015 m, which is hosting in the brittle deformation of Zhaoping fracture zone.
Introduction
Jiaodong type gold deposit, which is the primary type of gold endowment in the North China Craton (NCC), generally formed after prolonged stabilization of the craton (, ). This type deposit occurs under transpressinal condition from a transpressional to transtensional tectonic setting, which induced by Late Mesozoic breakoff of the subduction slab and rollback of the Paleo-Pacific Plate (; ). The metallogenic law has an evidently relationship with the NNE-NE-trending oblique-slip faults (; ; ). Generally, the formation of cratons remained stable in the early Precambrian and hardly triggered large-scale tectono-magmatic activities and gold mineralization (). However, multi-disciplinary observations indicate that the lithosphere of NCC has experienced unsteady flow. Two important gold metallogenic belts (the East gold belt and the West gold belt, Figure 1) formed during the metasomatic transformation of mantle-derived fluids and magmatic activities (; ) (Figure 1).
FIGURE 1
The Jiaodong gold province, with more than 5000 t of proven gold resources, is considered as the only global district that preserved giant gold resources in the East gold belt (
Here we present high-precision gravity and CSAMT methods to provide insights into spatial distribution of stratum, structure and magmatite for the Shanhou gold deposit hosting in southern part of Zhaoping fault. Moreover, relatively advanced data processing and inversion techniques help us to elucidate the coupling relationship between gold mineralized alteration zones and smooth-out waveform faults.
Geological setting
Regional geology
The Jiaodong gold province extends between the eastern NCC and northern SuLu metamorphic belt (Figure 1A) (
Shanhou gold deposit is located in the southern part of the Zhaoping fault and is dominated by middle Early Cretaceous Guojialing granite, Archean Jiaodong Group gneiss, and Paleoproterozoic metamorphic rocks of the Jinshan Group (Figure 2). Three gold deposits (Xiadian, Beipo and Shanhou) are hosted in the granitic cataclastic rock and sericitized granitic cataclastic rock (Figure 2). Gold ore-bodies are typically controlled by brittle deformation structures, and show ductile-brittle deformation with sheared ores, altered rocks or pyrite aggregates, such as the elongated pyrite aggregates (Figure 2) (
FIGURE 2

Simplified regional geological map of Shanhou gold deposit. 1- Quaternary; 2- Paleoproterozoic Jinshan Group; 3- Late Jurassic Linglong granite; 4- Paleoproterozoic Shuangdingshan gneiss; 5- Paleoproterozoic Penglai metabasite; 6- Archean Qixia gneiss; 7-Archean Luanjiazhai metabasite; 8- Deep-seated fault; 9-Reservoir; 10- Secondary fracture; 11- Geophysical profile; 12- Gold deposit.
Geophysical analytical techniques
Petrophysical parameters
Significant petrophysical parameters differences are the prerequisite of surveying geophysical research. Table 1 lists the petrophysical parameters in density and specific resistivity. The collected samples include the ore-bearing rocks and wall rocks in the Shanhou gold deposit. Samples from the Paleoproterozoic Jinshan Group yielded the relatively high density (2.82 g/cm3) and high specific resistivity (4000 Ω m). The Linglong granite was characterized by low density (2.57 g/cm3) and high specific resistivity (8000 Ω m), and the ore-bearing rocks by intermediate density (2.73 g/cm3) and lowest specific resistivity (200 Ω m) (Table 1). This results are generally consistent with region petrophysical parameters statistics in Zhaoping fault (
TABLE 1
| Lithology | Number | Density | Specific resistivity |
|---|---|---|---|
| Common value | Common value | ||
| g/cm3 | (Ω•m) | ||
| Paleoproterozoic Jinshan Group | 23 | 2.82 | 4000 |
| Linglong granite | 30 | 2.57 | 8000 |
| Archean Qixia gneiss | 21 | 2.77 | 3000 |
| Paleoproterozoic Lianzhou metabasite | 6 | 2.93 | 7000 |
| Ore-bearing rock | 5 | 2.73 | 200 |
Density and specific resistivity data for the samples from the ore-bearing rocks and wall rocks in the Shanhou gold deposit.
High-precision gravity
Three high-precision gravimetric observation curves (Lines 16, 20 and 24) were collected along a section of 4.80 km long. These gravity profiles passes through the northeast–southwest–trending Zhaoping fault with the trend of NW. Each section included 33 continental gravity survey points at intervals of 50 m. It was collected by CG-5 type Gravimeter (Scintrex Corporation, Canada), which has the resolution less than 1 μGal, repeatability less than 5 μGal, and the static drift less than 0.02 mGal per day.
The two point five dimension man machine interactive inversion technology was applied to probe the deep-seated gold deposit (Figure 3). Firstly, the initial models are build according to the exiting gold metallogenic model. Then, the gravity forward modeling curves are subsequently calculated according to the initial models (
FIGURE 3

Schematic diagram of the two point five dimension level prism, modified after
Controlled source audio frequency magnetotellurics
CSAMT profile data were collected along a northwest–southeast-trending section of 3.10 km long (Figure 4). This section included 63 CSAMT survey points at intervals of 50 m. It was surveyed with broadband, multi-channel digital electromagnetic exploration system named GDP-32II (Zonge Corporation, the United States of America). The GDP-32II receiver is an 16-bit±1/2 low power consumption, high-accuracy data collected system, with one 250 mv/nT pass-band sensitivity ANT-6 senor, and 16 solid non-polarized poles for collecting electromagnetic signal instruments. The horizontal electric dipole transmitter was setting at E 120°18′12.02″, N 37°05′36.32″, and the emission current was about 30 A.
FIGURE 4

The CSAMT field observation system.
The results of smooth Cagniard resistivity and impedance phase indicate that the field survey data is true and effective (Figure 5). The commercial software gather the observed data and build a new model section, including line annotation, survey configuration, and inversion control. After multiple iterations, the initial inversion resistivity is obtain. Then the software automatically updates the model section and obtains the further inversion resistivity. Finally, the inversion errors (about 3.47%) satisfy the requirements (less than 7%).
FIGURE 5

The pseudosection of the CSAMT profile; (A) Cagniard resistivity; (B) Impedance phase. The X-axis shows the survey point and the Y-axis shows the apparent depth.
Results
High-precision gravity
The amplitudes of high-precision gravity data are determined by two factors: lithology and fault. The overall trend of three gravity profiles are rising and degrading in some regional survey points. The black, red and green lines represent the field surveyed gravity curves. The yellow lines are the theoretical gravity curves which are calculated by two point five dimensional man-machine interactive inversion technology (Figures 6A–C). The assumed geological bodies were adjusted until the theoretical gravity curve is consistent with the measured gravity data (
FIGURE 6

Man machine interactive inversion of three gravity profiles. (A–C) Observation curves and fitting curves; (A’–C’) Geologic interpretation. 1- Observation curve 16 line; 2- Observation curve 20 line; Observation curve 24 line; 4- Fitting curve; 5- Paleoproterozoic Jinshan Group; 6- Archean Qixia gneiss; 7-Late Jurassic Linglong granite; 8- Paleoproterozoic Lianzhou metabasite; 9- gold ore bodies.
Controlled source audio frequency magnetotellurics profile
The CSAMT profile vertically passes through the trend of different geological bodies, producing a high resolution geoelectric structure (Figure 7). Laterally, two dimensional NLCG inversion results can be divided into three segments. The first one (sites 10–13) with intermediate to high apparent resistivity (∼103.6 Ω m) should be corresponding to Late Jurassic Linglong granite. The second one (sites 15–28) with intermediate apparent resistivity (∼103.2 Ω m) is the Archean gneisses. At survey point 14, an abrupt change in the apparent resistivity (104 to 102.2 Ω m) is interpreted as a regional brittle fracture (named Zhaoping fault), which separates Linglong Granite from the Archean Qixia gneiss (Figure 7A). The third segment (sites 29–42) with highest apparent resistivity (∼104.2 Ω m) is the Paleoproterozoic Jinshan Group. At survey point 28, Heihu fault with the low apparent resistivity (∼102.2 Ω m) separates Archean Jiaodong Group from Paleoproterozoic Jinshan Group (Figure 7A).
FIGURE 7

Inversion resistivity model of CSAMT by NLCG. (A) Line 20; (B) Line 18; (C) Line 22; 1- Ore-controlling brittle fracture; 2- The proven superficial gold bodies; 3- Drilled hole; 4- Recommended drilling.
Vertically, the inversion results can be segmented into three layers. The first layer (depths of 500 m) represents intermediate to high apparent resistivity (103.5–104.2 Ω m) geological bodies, including Late Jurassic Linglong granite, Archean Qixia gneiss and Paleoproterozoic Jinshan Group. The second layer (depths of 500–1000 m) with the lowest apparent resistivity (∼102.0 Ω m) should be corresponding to deep-seated Zhaoping fault. This undulating fault no only controls the superficial Shanhou gold deposit, but also controls the deep-seated blind gold ore-bodies. According to the location of surface exposure and the apparent resistivity high value, the deeper layer should be consisted by Late Jurassic Linglong granite (Figure 7–9).
FIGURE 8

Geological map of tunnel samplings. (A) At the depth of −50 m; (A) At the depth of −100 m; 1- Archean plagioclase gneiss; 2- Sericitization granite; 3- Fault gouge and mylonite; 4- gold ore bodies; 5- tunnel; 6- geological profile.
FIGURE 9

Geological profile of exploration in Shanhou gold deposit. (A) Line 20; (B) Line 18; (C) Line 22; 1- Quaternary; 2- Archean plagioclase gneiss; 3- Archean biotite granulite; 4- Paleoproterozoic dolomitic marble; 5- Sericitized cataclastic granite; 6- Pyrite sericite quartz; 7- Alteration zone; 8- gold bodies; 9- Occurrence ; 10- Completed drilling; 11- Recommended drilling.
Discussion
Spatial correlation between gold mineralization and fault
The majority of Jiaodong type gold deposit, are characterized by massive gold ore bodies in deep-seated faults, such as Jiaodong Sanshandao gold deposit occurring in Sansandao fault, Liaodong Wulong gold deposit occurring in Jixingou fault, and Jilin Jiapigou gold deposit occurring in Jiapigou fault (
The regional NNE-NE-trending faults in the Jiaodong gold province are regarded as the subsidiary fault of Tan-Lu fault. They are initially determined as the main gold ore-guided structures instead of the ore-controlling structures, such as Sanshandao, Jiaojia, Zhaoping, Qixia from west to east, respectively (
Geodynamics of Jiaodong type gold deposit in Jiaodong gold province
Tectonic altered rock type or dissemination type gold deposits are commonly coupling with NNE-NE trending brittle faults, which are sub-parallel to the Tan-Lu fault to the west (
At plate tectonic scale, the Paleo-Pacific Plate had a slow northwesterly subduction earlier than ∼135 Ma and then gradually deflected in a clockwise direction during 135–125 Ma. At the period of 125–100 Ma, the abrupt breakoff and roll-back of the subduction slab triggered the destruction of NCC, resulted in extension and thinning of lithosphere (
FIGURE 10

Cartoon showing the geodynamics of Jiaodong type gold deposit in Jiaodong gold province. (A) Metallogenic model of Jiaodong gold province; (B) Lithospheric and crustal architecture for the effective deposition of auriferous magma, modified from
The far-field stresses of subduction of Paleo-Pacific Plate triggered the temporal evolution of igneous activity and the structural geometries of gold deposits (
Conclusion
1) Based on the results of the high-precision gravity and CSAMT profiles, smooth-out waveform geophysical anomaly regarded as ore-bearing zone was identified.
2) The smooth-out waveform fault is essential to the effective deposition of auriferous magma.
3) The mineralization of Jiaodong type gold deposits was commonly related to subduction of Paleo-Pacific Plate.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author contributions
Z-HX designed the surveys and carried them out. G-WG performed the geophysical data processing. Z-JD and G-WG performed the geophysical data processing. J-YJ, LW, and X-GN prepared the manuscript with contributions from all co-authors. Z-JD prepared the manuscript with contributions from all co-authors.
Funding
This research benefited from the support of Fundamental Research Funds for the Central Universities (ZY2299211143), Shandong Provincial Engineering Laboratory of Application and Development of Big Data for Deep Gold Exploration Funds (SDK202221), and Science and Technology Research Project in Higher Learning Institutions of Hebei Province (ZC2022106), and Engineering Research Center of Geothermal Resources Development Technology and Engineering (22006).
Acknowledgments
We thank Frontiers in Earth Science Topic Editors for their critical reviews of an earlier version of this paper, which helped us to improve the quality. We thank Jilin University Fengyue Sun, LW and Jilin Exploration Geophysics Institute Chief Engineer Fuwen Li for their helpful suggestions during manuscript preparation, This research has been funded by the Engineering Research Center of Geothermal Resources DevelopmentTechnology and Engineering, Ministry of Education, Jilin University.
Conflict of interest
Author X-GN was employed by Inner Mongolia Non-ferrous Geological Mining Industry.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
CaiY. C.FanH. R.SantoshM.HuF. F.YangK. F.LiX. H. (2018). Decratonic gold mineralization: Evidence from the Shangzhuang gold deposit, eastern North China Craton. Gondwana Res.54, 1–22. 10.1016/j.gr.2017.09.009
2
CoxS. F. (2010). The application of failure mode diagrams for exploring the roles of fluid pressure and stress states in controlling styles of fracture-controlled permeability enhancement in faults and shear zones. Geofluids10, 217–233. 10.1111/j.1468-8123.2010.00281.x
3
DengJ.WangC. M.BagasL.CarranzaE. J. M.LuY. (2015). Cretaceous–cenozoic tectonic history of the Jiaojia fault and gold mineralization in the Jiaodong Peninsula, China: Constraints from zircon U–Pb, illite K–Ar, and apatite fission track thermochronometry. Min. Depos.50, 987–1006. 10.1007/s00126-015-0584-1
4
DengJ.YangL. Q.GrovesD. I.ZhangL.QiuK. F.WangQ. F. (2020). An integrated mineral system model for the gold deposits of the giant Jiaodong province, eastern China. Earth-Science Rev.208 (2), 103274. 10.1016/j.earscirev.2020.103274
5
DengJ.YangL. Q.LiR. H.GrovesD. I.SantoshM.WangZ. L.et al (2019). Regional structural control on the distribution of world class gold deposits: An overview from the giant Jiaodong gold province, China. Geol. J.54, 378–391. 10.1002/gj.3186
6
FanZ. G.HuangX. Z.TanL.YangX.ZhangH.ZhouD.et al (2014). A study of iron deposits in the Anshan area, China based on interactive inversion technique of gravity and magnetic anomalies. Ore Geol. Rev.57, 618–627. 10.1016/j.oregeorev.2013.09.017
7
GoldfarbR.GrovesD. I.GardollS. (2001). Orogenic gold and geologic time: A global synthesis. Ore Geol. Rev.18, 1–75. 10.1016/s0169-1368(01)00016-6
8
GrovesD. I.BierleinF. P. (2007). Geodynamic settings of mineral deposit systems. J. Geol. Soc. Lond.164, 19–30. 10.1144/0016-76492006-065
9
GrovesD. I.SantoshM.ZhangL. (2020). A scale-integrated exploration model for orogenic gold deposits based on a mineral system approach. Geosci. Front.11 (3), 719–738. 10.1016/j.gsf.2019.12.007
10
GrovesD. I.VielreicherR. M.GoldfarbR. J.CondieK. C. (2005). Controls on the heterogeneous distribution of mineral deposits through time. Geol. Soc. Lond. Spec. Publ.248 (1), 71–101. 10.1144/gsl.sp.2005.248.01.04
11
GuoL. N.GoldfarbR. J.WangZ. L.LiR. H.ChenB. H.LiJ. L. (2017). A comparison of jiaojia-and linglong-type gold deposit ore-forming fluids: Do they differ?Ore Geol. Rev.88, 511–533. 10.1016/j.oregeorev.2016.12.003
12
HanJ. L.DengJ.ZhangY.SunJ.WangQ.ZhangY.et al (2021). Au mineralization-related magmatism in the giant Jiapigou mining district of Northeast China. Ore Geol. Rev.141, 104638. 10.1016/j.oregeorev.2021.104638
13
LiC.LiL.LiC. R.SantoshM.ShenJ. F. (2022). Geochemistry of hydrothermal zircon as a proxy to fingerprint ore fluids in late Mesozoic decratonic gold deposits. Ore Geol. Rev.143, 104703. 10.1016/j.oregeorev.2022.104703
14
LiJ., W.BiS. J.SelbyD.ChenL.VasconcelosP.ThiedeD.et al (2012). Giant Mesozoic gold provinces related to the destruction of the North China craton. Earth Planet. Sci. Lett.349-350, 26–37. 10.1016/j.epsl.2012.06.058
15
LiX. C.FanH. R.SantoshM.HuF. F.YangK. F.LanT. G. (2013). Hydrothermal alteration associated with Mesozoic granite-hosted gold mineralization at the Sanshandao deposit, Jiaodong Gold Province, China. Ore Geol. Rev.53, 403–421. 10.1016/j.oregeorev.2013.01.020
16
LiangY. Y.DengJ.LiuX.WangQ.MaY.GaoT.et al (2019). Water contents of early cretaceous mafic dikes in the Jiaodong Peninsula, eastern north China craton: Insights into an enriched lithospheric mantle source metasomatized by paleo–Pacific plate subduction–related fluids. J. Geol.127 (3), 343–362. 10.1086/702648
17
SuiS. W.YuC. C.YaoC. L. (2004). The semi-intelligent processing and interpretation software for gravity and magnetic anomalies along the profile of rolling topography and its application. Geophys. Geochem. Explor.28 (1), 65–68. 10.3969/j.issn.1000-8918.2004.01.017
18
WangZ. L.YangL. Q.DengJ.SantoshM.ZhangH. F.LiuY.et al (2014). Gold-hosting high Ba-Sr granitoids in the Xincheng gold deposit, Jiaodong Peninsula, East China: Petrogenesis and tectonic setting. J. Asian Earth Sci.95, 274–299. 10.1016/j.jseaes.2014.03.001
19
WeiW. S.MuT.ZhangH. H. (2001). Geological map of Jiaodong gold Province showing the metallogenic belts, ore deposit types and distribution of the gold deposits. Lang Fang, Hebei: Gold Geology Research Institute of China Armed Police Force.
20
WuJ. J.ZengQ. D.SantoshM.FanH.BaiR.LiX.et al (2022). Deep ore-forming fluid characteristics of the Jiaodong gold province: Evidence from the Qianchen gold deposit in the Jiaojia gold belt. Ore Geol. Rev.145, 104911. 10.1016/j.oregeorev.2022.104911
21
XuZ. H.SunF.XinW.SunN.LiF.NiuJ.et al (2019). formation and evolution of paleoproterozoic orogenic belt in southern Jilin, jiao–liao–ji belt, north China craton: Constraints from Geophysics. Precambrian Res.333, 105433. 10.1016/j.precamres.2019.105433
22
YangJ. H.WuF. Y.WildeS. A. (2003). A review of the geodynamic setting of large-scale Late Mesozoic gold mineralization in the North China Craton: An association with lithospheric thinning. Ore Geol. Rev.23, 125–152. 10.1016/s0169-1368(03)00033-7
23
YangL. Q.DengJ.WangZ. L.GuoL. N.LiR. H.GrovesD. I.et al (2016). Relationships between gold and pyrite at the xincheng gold deposit, Jiaodong Peninsula, China: Implications for gold source and deposition in a brittle epizonal environment. Econ. Geol.111, 105–126. 10.2113/econgeo.111.1.105
24
YangL. Q.DilekY.WanZ. L.WeinbergR. F.LiuY. (2018). Late jurassic, high Ba–Sr Linglong granites in the Jiaodong Peninsula, east China: Lower crustal melting products in the eastern north China craton. Geol. Mag.155, 1040–1062. 10.1017/s0016756816001230
25
YaoC. L.LiY. S.GuanZ. N. (1998). Improvement of real-time and visualized forward and inversion of gravity andmagnetic anomalies. Geoscience12 (1), 115–122.
26
YinJ. N.LindsayM.TengS. R. (2018). Mineral prospectivity analysis for BIF iron deposits: A case study in the Anshan-Benxi area, Liaoning province, North-East China. Ore Geol. Rev.120, 102746. 10.1016/j.oregeorev.2018.11.019
27
ZhangL.GrovesD. I.YangL. Q.WangG. W.LiuX. D.LiD. P.et al (2020). Relative roles of formation and preservation on gold endowment along the Sanshandao gold belt in the Jiaodong gold province, China: Importance for province- to district-scale gold exploration. Min. Depos.55 (2), 325–344. 10.1007/s00126-019-00908-1
28
ZhangZ.WangG.CarranzaE.ZhangJ.TaoG.ZengQ.et al (2019). Metallogenic model of the Wulong gold district, China, and associated assessment of exploration criteria based on multi-scale geoscience datasets. Ore Geol. Rev.114, 103138. 10.1016/j.oregeorev.2019.103138
29
ZhuR.SunW. (2021). The big mantle wedge and decratonic gold deposits. Sci. China Earth Sci.64 (9), 1451–1462. 10.1007/s11430-020-9733-1
30
ZhuR. X.ChenL.WuF. Y.LiuJ. (2011). Timing, scale and mechanism of the destruction of the North China Craton. Sci. China Earth Sci.547, 789–797. 10.1007/s11430-011-4203-4
31
ZhuR. X.FanH. R.LiJ. W.MengQ.LiS.ZengQ. (2015). Decratonic gold deposits. Sci. China Earth Sci.58, 1523–1537. 10.1007/s11430-015-5139-x
Summary
Keywords
deep exploration, decratonic gold deposit, subduction of Paleo-Pacific Plate, Jiaodong gold province, geophysical method, cratonic destruction
Citation
Xu Z-H, Ding Z-J, Gu G-W, Jiang J-Y, Wang L and Niu X-G (2022) Deep exploration of Jiaodong type gold deposit, taking Shanhou gold deposit, southern part of Zhaoping fault as an example. Front. Earth Sci. 10:939375. doi: 10.3389/feart.2022.939375
Received
09 May 2022
Accepted
25 July 2022
Published
19 September 2022
Volume
10 - 2022
Edited by
Liang Zhang, China University of Geosciences, China
Reviewed by
Ahmed M. Eldosouky, Suez University, Egypt
Sirui Wang, China University of Geosciences, China
Updates

Check for updates
Copyright
© 2022 Xu, Ding, Gu, Jiang, Wang and Niu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zheng-Jiang Ding, 110357960@qq.com; Guan-Wen Gu, sun_ggw@163.com
This article was submitted to Economic Geology, a section of the journal Frontiers in Earth Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.