New Insight into the Genetic Mechanism of Shear Zone Type Gold Deposits from Muping-Rushan Metallogenic Belt (Jiaodong Peninsula of Eastern China)

: Most gold deposits are genetically controlled by shear zones, which are called shear zone type gold deposits (SZTGD). A better understanding of kinematics of shear zones and its constraint on the ore-forming process is critical to reveal the genetic mechanism of the SZTGD and favorable to mineral exploration. By conducting detailed structural analysis including ﬁeld and microscopic observations and electron backscatter di ﬀ raction (EBSD) and fractal dimension analysis in the Muping-Rushan shear zone (MR) as well as several gold deposits, the kinematic characteristics of the MR are well recognized and the metallogenic process of the SZTGD are discussed. The main conclusions are as follows: (1) petrology, geometry, kinematics, macro- and micro-structures imply that the MR has experienced a progressive shearing history exhumed via middle crust to subsurface level under the NW-SE extensional regime from late Jurassic to early Cretaceous; (2) in the MR, gold may precipitate both in the brittle fractures at middle crust level and brittle deformation part at shallow crust level during the stress-chemical process and (3) comparison of gold deposits between the MR and other areas show that the SZTGD has a uniform metallogenic mechanism, which is from (multi-stage) pluton emplacement, hydrothermal ﬂuid action, shearing action, brittle fracturing, sudden reduction of ﬂuid pressure, ﬂash vaporization to (gold) mineralization.


Introduction
Shear zones are one type of strain localization structures [1], which are widely developed in various structural environments. As weak zones in regions, shear zones are always concentrated for intense magma activities, fluid actions and metamorphism, hence provide favorable transport channels and precipitation spaces for ore-forming elements, e.g., Au, Ag, Cu, Pb, Zn, W, Sn and Mo [2][3][4][5][6][7]. Many gold deposits are host in, or spatially associated with first-and second-order shear zones [8][9][10][11][12][13][14][15], which we call the shear zone type gold deposits (SZTGD). In the SZTGD, the formation, transport and precipitation of gold are strictly controlled by shearing activities [16,17]. This kind of gold deposit has a wide distribution in the world, e.g., Abitibi greenstone belt in Canada, the Norseman-Wiluna Muping-Rushan metallogenic belt marks the eastern part of Jiaodong Peninsula. Exposed rocks of this region comprise the Archean Jiaodong Group, the Paleoproterozoic Jingshan Group and widespread Mesozoic granitoid plutons ( Figure 1). Evidence from petrological, geochemical and isotopic data suggests two major periods of magmatism during the Mesozoic. The late Jurassic Kunyushan granitoids, emplaced in 160-155 Ma, are interpreted as originating from partial melting of thickening crust [51]. The early Cretaceous Sanfoshan granitoids intruded the Kunyushan granitoids in 115-111 Ma and were formed by crust-mantle interaction during the lithosphere thinning [52]. At late Cretaceous (110-100 Ma), these two kinds of granitoids were intruded by numerous mafic dykes and pyroxene diorite stocks that originated from a rich lithospheric mantle contemporaneously with the granites [53].
Gold deposits in this belt develop along NE-NNE-striking shear zones ( Figure 2). From north to south, many large gold deposits developed there including Denggezhuang, Jinniushan, Rushan, Yinggezhuang and Sanjia. The main mineralization type is quartz-vein type gold deposits, and the individual ore vein is usually 40-350 m in length, 1-10 m in thickness and extends 40-70 m even kilometers in length [26]. Zircon U-Pb and K-feldspar 40 Ar/ 39 Ar constrain the age of mineralization for Denggezhuang deposit at 120 ± 3. 6 Ma and 123 ± 0. 53 Ma, respectively [51,54]. Zircon U-Pb and Sericite 40 Ar/ 39 Ar dating show 117 ± 3 Ma and 129 Ma of Rushan gold deposit, respectively [54,55]. Besides, K-feldspar 40 Ar/ 39 Ar dating shows 111.24 ± 0. 39 Ma of Yinggezhuang gold deposit [56]. Muping-Rushan metallogenic belt marks the eastern part of Jiaodong Peninsula. Exposed rocks of this region comprise the Archean Jiaodong Group, the Paleoproterozoic Jingshan Group and widespread Mesozoic granitoid plutons ( Figure 1). Evidence from petrological, geochemical and isotopic data suggests two major periods of magmatism during the Mesozoic. The late Jurassic Kunyushan granitoids, emplaced in 160-155 Ma, are interpreted as originating from partial melting of thickening crust [51]. The early Cretaceous Sanfoshan granitoids intruded the Kunyushan granitoids in 115-111 Ma and were formed by crust-mantle interaction during the lithosphere thinning [52]. At late Cretaceous (110-100 Ma), these two kinds of granitoids were intruded by numerous mafic dykes and pyroxene diorite stocks that originated from a rich lithospheric mantle contemporaneously with the granites [53].
Gold deposits in this belt develop along NE-NNE-striking shear zones ( Figure 2). From north to south, many large gold deposits developed there including Denggezhuang, Jinniushan, Rushan, Yinggezhuang and Sanjia. The main mineralization type is quartz-vein type gold deposits, and the individual ore vein is usually 40-350 m in length, 1-10 m in thickness and extends 40-70 m even kilometers in length [26]. Zircon U-Pb and K-feldspar 40 Ar/ 39 Ar constrain the age of mineralization for Denggezhuang deposit at 120 ± 3. 6 Ma and 123 ± 0.53 Ma, respectively [51,54]. Zircon U-Pb and Sericite 40 Ar/ 39 Ar dating show 117 ± 3 Ma and 129 Ma of Rushan gold deposit, respectively [54,55]. Besides, K-feldspar 40 Ar/ 39 Ar dating shows 111.24 ± 0.39 Ma of Yinggezhuang gold deposit [56]. These results were consistent with the regional gold metallogenic events (120 ± 10 Ma), which suggest gold mineralization was formed in Mesozoic under the condition of tectonic regime transition, strong thinning of lithosphere, activation of a series of detachment shear zones and metamorphic core complexes (MCCs) and large scale metallogeny [15,26,54]. These results were consistent with the regional gold metallogenic events (120 ± 10 Ma), which suggest gold mineralization was formed in Mesozoic under the condition of tectonic regime transition, strong thinning of lithosphere, activation of a series of detachment shear zones and metamorphic core complexes (MCCs) and large scale metallogeny [15,26,54]. The Muping-Rushan shear zone (MR) is considered as the second-order fracture of strike-slipping Tan-Lu fault, which is a crust scale fault along the East China continental margin (ECCM) and experienced a transition from sinistral at Late Jurassic to normal faulting at Early Cretaceous [58]. The MR is mainly located along the boundary between the Paleoproterozoic metamorphic rocks and the Kunyushan pluton and comprises four NNE-striking shear zones including Qinghushan-Tangjiagou (QT), Shigou-Wushan (SW), Jiangjunshi-Quhezhuang (JQ) and Gekou (GK) shear zones from west to east [52]. This study focused on the first three shear zones due to their close relationship with gold mineralization. Field works show that these shear zones mostly strike NE about 0-15° and dip SE with an angle around 65-85°, locally dip NW [57]. Most of them are considered as normal faults with the brittle deformation regime [15,52,57]. The Muping-Rushan shear zone (MR) is considered as the second-order fracture of strike-slipping Tan-Lu fault, which is a crust scale fault along the East China continental margin (ECCM) and experienced a transition from sinistral at Late Jurassic to normal faulting at Early Cretaceous [58]. The MR is mainly located along the boundary between the Paleoproterozoic metamorphic rocks and the Kunyushan pluton and comprises four NNE-striking shear zones including Qinghushan-Tangjiagou (QT), Shigou-Wushan (SW), Jiangjunshi-Quhezhuang (JQ) and Gekou (GK) shear zones from west to east [52]. This study focused on the first three shear zones due to their close relationship with gold mineralization. Field works show that these shear zones mostly strike NE about 0-15 • and dip SE with an angle around 65-85 • , locally dip NW [57]. Most of them are considered as normal faults with the brittle deformation regime [15,52,57].

Ductile Deformation of the MR
Owing to the strong weathering and cultivation on the surface, there are lack of high-quality outcrops and well-documented deformation. However, we tracked along the MR and found some good outcrops. Evidence for ductile deformation in the MR is preserved by kinematic indicators involving S-C fabrics, sigmoidal porphyroid structures, domino structures and asymmetric folds ( Figure 3). gneissic granite (Figure 3a). The thin-banded mylonitic foliation formed by feldspar and quartz aggregates dips to the NW about 295° with a dip angle around 30°. The lineation of the mylonite and the hanging-wall are uneasy to observe in the field due to strong erosion. Another well-documented outcrop was found in the middle part of the MR with the ductile formation (near location 7 in Figure 2). The sub-horizontal (about 15°) thin-banded mylonite was intruded by felsic veins exhibiting domino structure (Figure 3b). The lineation formed by elongated quartz bands can be observed on the foliation surface and plunges towards 300° with plunge angles between 12 and 15°. The combination of the lineation and the domino structure of felsic veins also suggest the top-to-the-NW sense of shear. In another outcrop (location 1 in Figure 2), the gneissic foliation formed by feldspar, quartz and biotite dips to the NW about 287° with a dip angle around 37° and the lineation plunges towards 305° on the foliation surface ( Figure 3c). Aggregates including quartz and feldspar exhibit sigmoidal shapes, indicating a top-to-the-NW sense of shear (Figure 3d). Some asymmetric folds in calcite veins with fibers perpendicular to oblique to the vein walls also show the same shear sense, as illustrated in Figure 3e.  There is a narrow shear zone about 30 cm wide in location 10 ( Figure 2) above the weakly gneissic granite (Figure 3a). The thin-banded mylonitic foliation formed by feldspar and quartz aggregates dips to the NW about 295 • with a dip angle around 30 • . The lineation of the mylonite and the hanging-wall are uneasy to observe in the field due to strong erosion. Another well-documented outcrop was found in the middle part of the MR with the ductile formation (near location 7 in Figure 2). The sub-horizontal (about 15 • ) thin-banded mylonite was intruded by felsic veins exhibiting domino structure (Figure 3b).
The lineation formed by elongated quartz bands can be observed on the foliation surface and plunges towards 300 • with plunge angles between 12 and 15 • . The combination of the lineation and the domino structure of felsic veins also suggest the top-to-the-NW sense of shear. In another outcrop (location 1 in Figure 2), the gneissic foliation formed by feldspar, quartz and biotite dips to the NW about 287 • with a dip angle around 37 • and the lineation plunges towards 305 • on the foliation surface ( Figure 3c). Aggregates including quartz and feldspar exhibit sigmoidal shapes, indicating a top-to-the-NW sense of shear (Figure 3d). Some asymmetric folds in calcite veins with fibers perpendicular to oblique to the vein walls also show the same shear sense, as illustrated in Figure 3e The principle extension direction was determined for the development of the MR by using fault-slip data. Following the calculation method proposed by [59,60], the principal stress orientations (σ 1 , σ 2 and σ 3 ) were inverted from the measured fault-slip datasets, using software SG2PS developed by [61]. From this technique we can constrain the principle extension direction (σ 3 ) recorded on the study faults, although it is still in debate [49]. All the average angular difference (ANG) between predicted and actual striations are less than 25 • , which suggests that the calculated stress orientations are reliable [62]. Calculating results ( Figure 2) show that the principle extension direction is consistent with the development of the brittle faults in the MR, which is under NW-SE extension. It should be noted that the extension direction of normal faults is also consistent with that indicated by mineral elongation lineation in the MR, suggesting that both ductile shear zones and brittle faults in the MR developed under a consistent stress field.

Microscopic Observations in the MR
The rock deformation mechanisms (e.g., ductile and brittle deformation) and conditions (e.g., temperature) can be identified by the microstructural features of deformed minerals (especially for quartz and feldspar) [63,64]. Deformed samples from the MR were used for the microscopic observation.
In order to get better information about the deformation style and the kinematics of shear zones, several samples were chosen to approximate the principal section XZ of the finite strain ellipsoid so that the thin sections were cut perpendicular to the foliation and parallel to the lineation.
In the thin sections, both high-and low-temperature deformation microstructures of minerals were observed. In the granitic mylonite, myrmekitic textures were present along the boundary between K-feldspar and quartz ( Figure 4a). In the gneissic mylonite, quartz shows dynamic recrystallization with GBM ( Figure 4b). Additionally, S-C fabrics formed by elongated quartz band and oriented biotite indicate a top-to-the-NW shear sense. Chessboard subgrains were also present in some quartz grains ( Figure 4c). These characters suggest high temperature deformation (500-600 • C) [63]. Microstructures of the medium temperature (400-500 • C) deformation were abundant in collected samples. For example, deformation twins were common in plagioclase (

EBSD Analysis
The electron backscatter diffraction (EBSD) analysis is widely used in the study of quartz crystallographic preferred orientations (CPOs), deformation mechanism, dislocation and slip

EBSD Analysis
The electron backscatter diffraction (EBSD) analysis is widely used in the study of quartz crystallographic preferred orientations (CPOs), deformation mechanism, dislocation and slip system, pressure and temperature conditions of crystals [64][65][66][67]. The CPOs of quartz c-axis are consistent with the dominant slip system under different conditions, which can provide further information about the deformation of shear zones [63].
The EBSD analysis was performed on a FEI Quanta 450 field emission gun scanning electron microscope coupled with an Oxford Instruments HKL Nordlys II EBSD detector at the State Key Laboratory of Geological Processes and Mineral Resources, China University of Geosciences, Wuhan. XZ sections for EBSD mapping were chemical-mechanical polished with a 0.05 µm colloidal silica solution for about two hours. Analytical conditions were operated at 15 kV accelerating voltage and working distance of 25 mm. Samples were tilted at 70 • and a step size between 2 and 5 µm was used for EBSD mapping. All collected data ( Figure 5) was processed using the HKL Channel 5 software package and all pore figures were presented with an equal area, lower hemisphere projection. Grain boundary maps of quartz were used for the fractal analysis with a subgrain definition of 10 degrees misorientation [68,69].
information about the deformation of shear zones [63].
The EBSD analysis was performed on a FEI Quanta 450 field emission gun scanning electron microscope coupled with an Oxford Instruments HKL Nordlys II EBSD detector at the State Key Laboratory of Geological Processes and Mineral Resources, China University of Geosciences, Wuhan. XZ sections for EBSD mapping were chemical-mechanical polished with a 0.05 μm colloidal silica solution for about two hours. Analytical conditions were operated at 15 kV accelerating voltage and working distance of 25 mm. Samples were tilted at 70° and a step size between 2 and 5 μm was used for EBSD mapping. All collected data ( Figure 5) was processed using the HKL Channel 5 software package and all pore figures were presented with an equal area, lower hemisphere projection. Grain boundary maps of quartz were used for the fractal analysis with a subgrain definition of 10 degrees misorientation [68,69].
Results of the quartz c-axis fabrics from the MR mainly displayed point maximas ( Figure 5). Fabrics in samples 1030 and 3021 consisted of a maxima close to the Y-axis, as a result of the prism <a> slip under high temperature (500-650 °C). Maximas in sample 2917 were distributed between the Y-axis and base circle, indicative of the rhomb <a> slip under medium temperature (400-500 °C). The quartz c-axis fabrics in sample 2810 and 3019 show Z maxima, which is compatible with low-temperature basal <a> slip (300-400 °C). What is more, the CPOs of quartz c-axis in sample 2914 and 2916 show complex transitions, e.g., from Y maxima to maximas between the Y-axis and base circle in sample 2914, and from maximas between the Y-axis and base circle to Z maxima, respectively. These CPO patterns were attributed to deformation superposition during the progressive shearing of the MR, which was consistent with the microscopic observations. Remarkably, these asymmetric CPOs of quartz c-axis are common during non-coaxial progressive deformation, and angles between the X-axis and the tilt direction of the CPOs indicate the top-to-the-NW sense of shear [70][71][72].  Results of the quartz c-axis fabrics from the MR mainly displayed point maximas ( Figure 5). Fabrics in samples 1030 and 3021 consisted of a maxima close to the Y-axis, as a result of the prism <a> slip under high temperature (500-650 • C). Maximas in sample 2917 were distributed between the Y-axis and base circle, indicative of the rhomb <a> slip under medium temperature (400-500 • C). The quartz c-axis fabrics in sample 2810 and 3019 show Z maxima, which is compatible with low-temperature basal <a> slip (300-400 • C). What is more, the CPOs of quartz c-axis in sample 2914 and 2916 show complex transitions, e.g., from Y maxima to maximas between the Y-axis and base circle in sample 2914, and from maximas between the Y-axis and base circle to Z maxima, respectively. These CPO patterns were attributed to deformation superposition during the progressive shearing of the MR, which was consistent with the microscopic observations. Remarkably, these asymmetric CPOs of quartz c-axis are common during non-coaxial progressive deformation, and angles between the X-axis and the tilt direction of the CPOs indicate the top-to-the-NW sense of shear [70][71][72].

Fractal Analysis of Dynamically Recrystallized Quartz Grains
A better understanding of the deformation temperature and strain ratio is critical to the evolution history and kinematics of shear zones. Despite EBSD analysis, fractal analysis of dynamically recrystallized quartz grains is proved as an effective method to get the physical conditions of shear zones [73][74][75]. In a simple sense, the grain boundary (perimeter) and grain size (calculated as diameter of equivalent circle) of recrystallized quartz have a statistically linear correlation in two dimensions, and the slope represents fractal dimensions D. The D-value decreases with increasing temperature during formation of shear zones. So, on a statistical basis, D can be applied as a geothermometer for shear zones [76][77][78].
Four mylonite samples from the main shear zones of the study area are employed to do the analysis. All quartz grains have recrystallized boundaries with sub-grains rotation mechanism. Photos from microscopic analysis and grain boundary maps from EBSD analysis were used for the dense area of recrystallization, respectively. Next use drawing software to do the vectorization and count the actual perimeter P and area A for each deformed grain (more than 50 grains). The area A is used to calculate the diameter d of circle having the same area. Then by plotting lgP on the vertical axis and lgd on the horizontal axis, the fractal dimension D is represented as the slope of a least-squares fitted line. The range of D is 1 ≤ D ≤ 2 for two-dimensional analysis.
The fractal analysis results of the MR are shown in Table 1 and Figure 6. The correlation coefficient revealed a value more than 0.9, and D was in the range of 1.15-1.22, which indicate the boundary shape for the measured grain size was self-similar statistically. Comparing the Kural thermometer [73] and the relationship diagram between fractal dimension and deformation temperature (Figure 7), it was obvious that the deformation temperature was in the range of 400-600 • C, corresponding to the metamorphic condition from high greenschist facies to low amphibolite facies. Among them, the deformation temperature of data from the EBSD analysis were much higher than microscopic data, which might suggest that the EBSD data could give a much more accurate classification (Figures 6  and 7).      Figure 7. Relationship between fractal dimension and deformation temperature (after [73]).

Deposit Geology
Most gold deposits, e.g., Rushan, Denggezhuang and Yinggezhuang deposits, are distributed along the NNE-striking MR (Figure 1). Several large-scale gold veins were identified in the Rushan deposit located in the main shear zone (Figure 8a). While some other veins in the Denggezhuang and Yinggezhuang deposits were located in the secondary shear zones (Figure 8b,c). Moreover, section maps from drill-holes show that auriferous parts were located in the occurrence-change portion in a whole quartz vein, especially in the dip-turn-to-gentle positions (Figure 8d,e). The gold-quartz veins were characterized by strong silicification and sulfidation, together with sericitization and K-feldspar wallrock alteration (Figure 9a,b). In the middle of veins, quartz was fractured into breccia surrounded by disseminated sulfides, e.g., pyrite (Figure 9a). There were many small-scale auriferous veins near the main shear zone (Figure 9b) and the occurrence analysis from the limited data suggests that they were subsidiary fractures, including R and T fractures (Figure 9g,h). Most gold deposits, e.g., Rushan, Denggezhuang and Yinggezhuang deposits, are distributed along the NNE-striking MR (Figure 1). Several large-scale gold veins were identified in the Rushan deposit located in the main shear zone (Figure 8a). While some other veins in the Denggezhuang and Yinggezhuang deposits were located in the secondary shear zones (Figure 8b,c). Moreover, section maps from drill-holes show that auriferous parts were located in the occurrence-change portion in a whole quartz vein, especially in the dip-turn-to-gentle positions (Figure 8d,e). The gold-quartz veins were characterized by strong silicification and sulfidation, together with sericitization and K-feldspar wallrock alteration (Figure 9a,b). In the middle of veins, quartz was fractured into breccia surrounded by disseminated sulfides, e.g., pyrite (Figure 9a). There were many small-scale auriferous veins near the main shear zone (Figure 9b) and the occurrence analysis from the limited data suggests that they were subsidiary fractures, including R and T fractures (Figure 9g,h).   Microscopic observations show that ore minerals for these deposits include pyrite, chalcopyrite and native gold and most of them experienced strong fracturing (Figure 9c,d). Native gold is mainly preserved at open and healed fractures in pyrite (Figure 9c,d). Ductile deformation microstructures at medium temperature are abundant in the wallrocks close to ore bodies (Figure 9e,f), e.g., ribbon-shaped subgrain rotation recrystallization in quartz and deformation twins in plagioclase. Other evidence from the fluid inclusion research of ore bodies indicate the depositing temperature of 170-324 °C and pressure of 92-269 MPa [38,53]. Microscopic observations show that ore minerals for these deposits include pyrite, chalcopyrite and native gold and most of them experienced strong fracturing (Figure 9c,d). Native gold is mainly preserved at open and healed fractures in pyrite (Figure 9c,d). Ductile deformation microstructures at medium temperature are abundant in the wallrocks close to ore bodies (Figure 9e,f), e.g., ribbon-shaped subgrain rotation recrystallization in quartz and deformation twins in plagioclase. Other evidence from the fluid inclusion research of ore bodies indicate the depositing temperature of 170-324 • C and pressure of 92-269 MPa [38,53].

Geological Characteristics of the MR
It has been widely proposed that NCC has undergone large-scale lithosphere thinning in the extensional regime since late Mesozoic [29,67,79]. A series of NE-NNE-striking detachment shear zones and MCCs have been identified in Jiaodong Peninsula during this period, e.g., the Linglong MCC [80], the Wulian MCC [81], the Queshan MCC [67] and the Kunyushan MCC [15]. In this study, the MR was located at the western part the Kunyushan MCC (Figure 1b). It is well-known that the exhumation of MCCs is usually accompanied by the intrusion of plutons and ductile deformation along detachment shear zones. However, previous studies largely focused on the high-angle normal faults in the MR and considered them as a part of the hanging-wall of the detachment shear zone [82]. Evidence of ductile deformation and evolution history in the MR has been lacking until now.
Our detailed analysis from field and microscopic observations reveal that the MR has experienced progressive shearing. Given the lack of direct age data for the detachment shearing in the MR, it is need to combine the age data of intrusion for the Kunyushan and Sanfoshan plutons and the contact relationship between the plutons and shear zones to accurately constrain the exhumation of the MCC. Based on the structural analysis and previous dating results, the evolution of the MR is constructed into two stages: early stage high-temperature ductile shearing and late stage low-temperature brittle deformation.
(1) In late Jurassic (160-155 Ma), magmatism in the Muping-Rushan area started its activity and the Kunyushan pluton emplaced. During this period, the ductile detachment shear zone (the MR) of the Kunyushan MCC was initiated between the pluton and Precambrian basement under the top-to-the-NW sense of extension (Figure 10a). S-C fabrics, asymmetric folds and sigmoidal porphyroid structures were well developed in this shear zone. Evidence from microscopic deformation fabrics shows high temperature condition (500-600 • C). Given 33 • C/km as the geothermal gradient [77], the corresponding depth of shearing was in the range of 15-18 km, which belongs to the middle crust [3].
(2) From late Jurassic to early Cretaceous, the MR was progressively exhumed via middle crust to subsurface level in a NW-SE extensional regime (Figure 10b). A series of brittle faults superimposed on the MR. Domino structures, conjugate fractures, slickensides and steps are abundant in the shear zone and minerals experienced middle-to-low temperature deformation and commonly show cataclastic texture. At the late period of early Cretaceous (115-111 Ma), the second stage emplacement began its activity, forming Sanfoshan pluton. The kinematic characters of the MR are consistent with other regional shear zones in Jiaodong Peninsula during early Cretaceous, e.g., the Zhaoyuan-Pingdu shear zone in the northwest part of Jiaodong Peninsula [30,83] and the Queshan shear zone in the earth part of Jiaodong Peninsula [67]. Macro and micro observations show that these NE-NNE-striking shear zones experienced NW-SE extension during gold mineralization. It has been widely proposed that the westward subduction, roll-back and trench retreat of the Paleo-Pacific plate (the Izanagi plate) gave rise to tectonics of Jiaodong Peninsula while these shear zones were active [15,29]. Dehydration of the subducted slab, partial melting of the mantle and the following intense magmatism play important roles in the activities of these shear zones and gold mineralization.

Relationship Between the MR and Gold Mineralization
Previous studies have revealed that the gold mineralization is coeval with structural activities in Jiaodong deposits indicating the coupled behavior between them [15,29]. However, the The kinematic characters of the MR are consistent with other regional shear zones in Jiaodong Peninsula during early Cretaceous, e.g., the Zhaoyuan-Pingdu shear zone in the northwest part of Jiaodong Peninsula [30,83] and the Queshan shear zone in the earth part of Jiaodong Peninsula [67]. Macro and micro observations show that these NE-NNE-striking shear zones experienced NW-SE extension during gold mineralization. It has been widely proposed that the westward subduction, roll-back and trench retreat of the Paleo-Pacific plate (the Izanagi plate) gave rise to tectonics of Jiaodong Peninsula while these shear zones were active [15,29]. Dehydration of the subducted slab, partial melting of the mantle and the following intense magmatism play important roles in the activities of these shear zones and gold mineralization.

Relationship Between the MR and Gold Mineralization
Previous studies have revealed that the gold mineralization is coeval with structural activities in Jiaodong deposits indicating the coupled behavior between them [15,29]. However, the relationship between ore bodies and shear zones in the MR still need further study. Combining with previous researches and our study, it is suggested that most gold deposits are located in the brittle deformation part of the MR. More detailed observations indicate that ore bodies are concentrated in occurrence change portions of shear zones, e.g., the striking-turning portions in the geological maps and the dip turn-to-gentle positions in the sectional maps (Figure 8). Among them, a part of gold ore bodies are located in the secondary fractures of the MR, which belong to R and T fractures (Figure 9g,h). Others are formed in the main shear zones with brittle fractures and cataclasite (Figure 9a). It should be noticed that wallrocks close to ore bodies generally experienced medium-temperature ductile deformation (Figure 9e,f), while the ore bodies show brittle deformation (Figure 9c,d). Given these geological settings of the MR, there are two possible reasons for gold mineralization.
(1) At middle crust level, the development of ductile deformation of the MR makes it more permeable (Figure 10a). Therefore, the shear zone is more likely acting as a deep channel promoting the migration of ore-forming fluids, which were derived from cooling magma or mantle degassing during the period of craton destruction [29]. The fluid pressure was gradually accumulated at the stress concentration parts due to ductile deformation [84].
(2) At upper crust level, brittle deformation promoted the quick migration of fluids (Figure 10b). When they arrived at the occurrence change portions, the higher normal stress prevented the slip of shear zones, resulting in the stress concentration and the formation of high-pressure fluids [19].
The occurrence of high-pressure fluids in these two processes above can not only decrease rock strength [85], but also accelerate the strain rate [5]. Under high strain rate, brittle fractures (e.g., R and T fractures and cataclasite in the MR) occurred and decreased the fluid pressure rapidly causing flash vaporization. This process plays an important role to the physical-chemical state of the fluid, e.g., the solubility of quartz would decrease suddenly resulting in co-precipitation of silica with gold [86]. With the development of the MR, this process worked as a circulation and gradually caused gold accumulation in the shear zone.

Comparison to other Gold Deposits
Comparing the features of gold deposits both in Jiaodong Peninsula and all over the world, Goldfarb and Groves proposed the famous model of orogenic gold deposits and considered Jiaodong Peninsula as a typical example of this type of deposits [10,48,87,88]. Based on the definition in the literatures above, the typical orogenic gold deposits are formed in the orogenic process, including subduction and collision. According to the previous study, gold deposits in Jiaodong Peninsula were formed in the background of arc and back-arc extension during the subduction of Pacific Plate, which caused the destruction of North China Craton. From this view, someone also called it as the decratonic gold deposits and all shear zones were developed in the arc setting [29,89]. Considering all the characteristics above, there has been continuing debate on the genesis of gold deposits in Jiaodong Peninsula in accordance with the tectonic background. However, all models have one common feature that shear zones play important roles in the mineralization process. So, this paper prefers emphasizing the controlling effect of shear zones on gold mineralization, rather than the tectonic background. It is suggested that the term of shear zone type gold deposits is more appropriate.
As an essential economic source, the genetic of SZTGD has long been a major concern. Previous studies show that SZTGD can be developed in different structural environments, e.g., normal, thrust and strike-slip shearing. For the normal and strike-slip shearing, e.g., deposits in Jiaodong Peninsula, gold mineralization are dominant by the synchronous structure-magma-fluid activities [26,90]. Combining with our previous study [19] and this paper, most gold bodies are host in the main shear zones or secondary fractures although they belong to different mineralization type. For altered-rock type gold deposits, e.g., Jiaojia deposit in Zhaoyuan-Laizhou belt, gold ores are distributed both in the main shear zones with strong alteration, and secondary quartz veins recognized as R, T and R' fractures. For vein type gold deposits, e.g., Linglong deposit in Zhaoyuan-Laizhou belt and Denggezhuang and Rushan deposits in the Muping-Rushan belt, ore bodies are hosted in the main shear zones with strong breccia and cataclastic deformation, as well as secondary fractures. Another significant character is that all the gold bodies are located in the occurrence change portions of shear zones. As for thrust shearing, e.g., Abitibi greenstone belt in Canada and Norseman-Wiluna greenstone belt in west Australia, these sub-horizontal and extensional veins are recognized as T and R fractures, respectively [19]. Meanwhile, gold veins are distributed at the steeper portions of thrust shear zones. Moreover, the significant factor influencing gold precipitation causes widespread consideration. A fluid inclusion study shows that deposits from different structural environments are commonly formed in middle-temperature conditions (200-350 • C) [11,29,31,38]. Experimental researches show that gold solubility does not change much in this range [91]. In contrast, fluid pressure shows fluctuation in general and plays an important role in gold mineralization [15,86].
All the evidence above indicate that SZTGD are formed by a united genetic mechanism, which is called as the stress-chemical process from stress concentration, brittle fracturing, sudden reduction of fluid pressure, flash vaporization to gold precipitation [19]. Additionally, both vein-type and altered-rock type gold deposits may have a unified formation process. As an important mineralization type, altered-rock type gold deposits are widely distributed in shear zones, especially in Jiaodong Peninsula. According to our investigation for the typical altered-rock type deposits in Sanshandao, rocks were broken into breccia at fault center ( Figure 11a) and visible gold generally distributed along the boundary of breccia (Figure 11b). In some cases clasts of reworked cataclasite were cemented by fluid actions. recognized as R, T and R' fractures. For vein type gold deposits, e.g., Linglong deposit in Zhaoyuan-Laizhou belt and Denggezhuang and Rushan deposits in the Muping-Rushan belt, ore bodies are hosted in the main shear zones with strong breccia and cataclastic deformation, as well as secondary fractures. Another significant character is that all the gold bodies are located in the occurrence change portions of shear zones. As for thrust shearing, e.g., Abitibi greenstone belt in Canada and Norseman-Wiluna greenstone belt in west Australia, these sub-horizontal and extensional veins are recognized as T and R fractures, respectively [19]. Meanwhile, gold veins are distributed at the steeper portions of thrust shear zones. Moreover, the significant factor influencing gold precipitation causes widespread consideration. A fluid inclusion study shows that deposits from different structural environments are commonly formed in middle-temperature conditions (200-350 °C) [11,29,31,38]. Experimental researches show that gold solubility does not change much in this range [91]. In contrast, fluid pressure shows fluctuation in general and plays an important role in gold mineralization [15,86]. All the evidence above indicate that SZTGD are formed by a united genetic mechanism, which is called as the stress-chemical process from stress concentration, brittle fracturing, sudden reduction of fluid pressure, flash vaporization to gold precipitation [19]. Additionally, both vein-type and altered-rock type gold deposits may have a unified formation process. As an important mineralization type, altered-rock type gold deposits are widely distributed in shear zones, especially in Jiaodong Peninsula. According to our investigation for the typical altered-rock type deposits in Sanshandao, rocks were broken into breccia at fault center ( Figure 11a) and visible gold generally distributed along the boundary of breccia (Figure 11b). In some cases clasts of reworked cataclasite were cemented by fluid actions. Based on these distribution features above and the stress-chemical process, gold mineralization of these two types is genetically resulted from the stress concentration in shear zones. Firstly, microfractures occurred with the formation of high-pressure fluids at stress concentration parts (Figure 12a). Under the continuous effect of stress-chemical process, a cloud of microfractures grew (Figure 12b). When the density of these cracks reached a critical value, rocks were broken into cataclasite, which allowed shear displacement to occur and fluids fully contacted with wallrocks ( Figure 12c). In this condition, altered-rock type gold deposits might occur in shear zones. Otherwise, if the density of these cracks did not exceed a critical value (Figure 12d) and fractures were filled by quartz or other minerals, vein-type ones perhaps occurred (Figure 12e). Based on these distribution features above and the stress-chemical process, gold mineralization of these two types is genetically resulted from the stress concentration in shear zones. Firstly, microfractures occurred with the formation of high-pressure fluids at stress concentration parts (Figure 12a). Under the continuous effect of stress-chemical process, a cloud of microfractures grew (Figure 12b). When the density of these cracks reached a critical value, rocks were broken into cataclasite, which allowed shear displacement to occur and fluids fully contacted with wallrocks ( Figure 12c). In this condition, altered-rock type gold deposits might occur in shear zones. Otherwise, if the density of these cracks did not exceed a critical value (Figure 12d) and fractures were filled by quartz or other minerals, vein-type ones perhaps occurred (Figure 12e).

Genetic Mechanism for SZTGD
During the stress-chemical process, the coupling relationship among magma, fluids and shear zones plays an essential role in gold mineralization and determines special mineralization portions in shear zones. Extensive magmatism and multi-stage pluton emplacement can provide not only metallogenic material but also heat, which are beneficial for fluid activities. In Jiaodong Peninsula, the long-lasting magmatism is widely distributed and provides favorable conditions for gold mineralization [15,26]. Fluids derived from magmatism or metamorphic dehydration can make gold activate from wallrocks and migrate in forms of complexes, e.g., Au(HS)2 − along shear zones, especially crust-scale ones. Meanwhile, fluid action has an essential influence on rock deformation mechanism. The occurrence of high-pressure fluids facilitates the development of brittle fractures at different structural levels, including ductile, ductile-brittle and brittle shear zones. In this study, gold deposits in the MR (China) are probably hosted both in the brittle fractures at middle crust level (ductile deformation) and brittle deformation part at shallow crust level, while some deposits are distributed in ductile part (e.g., the Sigma Mine, Val d'Or, Canada, see [11]). Finally, rapid pressure drop occur in these fractures and gold begin to precipitate during the stress-chemical process. Considering shear-zone activity has a long-lived characteristic with consistent kinematics, later gold re-concentrations processes can occur along the shear-zones, associated with later activations without magmatic supply, increasing the tonnage.
From these processes above, we suggest that shear zones are the bonds connecting all activities together. That is the reason why most gold deposits host in shear zones. According to this metallogenic model, the stress concentration portions in shear zones are favorable areas for gold mines, and we may explore new deposits at deep ductile deformation part in the MR.

Genetic Mechanism for SZTGD
During the stress-chemical process, the coupling relationship among magma, fluids and shear zones plays an essential role in gold mineralization and determines special mineralization portions in shear zones. Extensive magmatism and multi-stage pluton emplacement can provide not only metallogenic material but also heat, which are beneficial for fluid activities. In Jiaodong Peninsula, the long-lasting magmatism is widely distributed and provides favorable conditions for gold mineralization [15,26]. Fluids derived from magmatism or metamorphic dehydration can make gold activate from wallrocks and migrate in forms of complexes, e.g., Au(HS) 2 − along shear zones, especially crust-scale ones.
Meanwhile, fluid action has an essential influence on rock deformation mechanism. The occurrence of high-pressure fluids facilitates the development of brittle fractures at different structural levels, including ductile, ductile-brittle and brittle shear zones. In this study, gold deposits in the MR (China) are probably hosted both in the brittle fractures at middle crust level (ductile deformation) and brittle deformation part at shallow crust level, while some deposits are distributed in ductile part (e.g., the Sigma Mine, Val d'Or, Canada, see [11]). Finally, rapid pressure drop occur in these fractures and gold begin to precipitate during the stress-chemical process. Considering shear-zone activity has a long-lived characteristic with consistent kinematics, later gold re-concentrations processes can occur along the shear-zones, associated with later activations without magmatic supply, increasing the tonnage. From these processes above, we suggest that shear zones are the bonds connecting all activities together. That is the reason why most gold deposits host in shear zones. According to this metallogenic model, the stress concentration portions in shear zones are favorable areas for gold mines, and we may explore new deposits at deep ductile deformation part in the MR.

Conclusions
(1) Petrology, geometry, kinematics, macro-and micro-structures imply that the MR had experienced a top-to-the-NW sense of extensional activity from late Jurassic to early Cretaceous. It had a progressive shearing history exhumed via middle crust to subsurface level.
(2) In the MR, gold might precipitate both in the brittle fractures at the middle crust level and brittle deformation part at the shallow crust level during the stress-chemical process.
(3) Comparison of gold deposits between the MR and other areas show that shear zone type gold deposits had a uniform metallogenic mechanism, which is from (multi-stage) pluton emplacement, hydrothermal fluid action, shearing action, brittle fracturing, sudden reduction of fluid pressure, flash vaporization to (gold) mineralization.