Mineralogy, Geochemistry, and Age Constraints on the Axinite-Bearing Gukjeon Pb–Zn Skarn Deposit in the Miryang Area, South Korea

: The axinite-bearing Gukjeon Pb–Zn deposit is hosted by the limestone, a member of the Jeonggaksan Formation, which, in turn, forms the part of the Jusasan subgroup of the Yucheon Group in the Gyeongsang Basin in the southeastern part of the Korean Peninsula. In this study, we attempted to interpret the spatial and temporal relationships among geologic events, including the mineralization of this deposit. We constructed a new 3D orebody model and suggested a relation-ship between skarn alteration and related mineralization. Mineralization timing was constrained using SHRIMP zircon age dating results combined with boron geochemistry on coeval intrusive rocks. Skarn alterations are restrictively found in several horizons of the limestone formation. The major skarn minerals are garnet (grossular), pyroxene (hedenbergite), amphibole (actinolite and ferro-actinolite), axinite (tizenite and ferro-axinite), and epidote (clinozoisite and epidote). The three stages of pre-skarn, syn-skarn, and post-skarn alteration are recognized within the deposit. The syn-skarn alteration is characterized by prograde metasomatic pyroxene and garnet, and the retrograde metasomatic amphibole, axinite, and epidote. Major skarn sulfide minerals are sphalerite, chalcopyrite, galena, and pyrite, which were predominantly precipitated during the retrograde stage and formed amphibole and axinite skarns. The skarn orebodies seem to be disc- or flat-shaped with a convex form at the central part of the orebodies. The vertical ascending and horizontal infiltration of boron-rich hydrothermal fluid probably controlled the geometry of the orebodies. Considering the whole-rock major, trace, and boron geochemical and geochronological results, the timing of Pb– Zn mineralization can be tightly constrained between the emplacement of boron-poor intrusion (fine-grained granodiorite, 82.8 Ma) and boron-rich intrusion (porphyritic andesite in Beomdori andesitic rocks, 83.8 Ma) in a back-arc basin setting. The boron for mineralization was sourced from late Cretaceous (Campanian), subduction-related magmatic rocks along the margin of the Pacific plate.


Introduction
The Gukjeon Pb-Zn deposit is located within the Gyeongsang Basin, in the southeastern Korean Peninsula. Unlike the Great Limestone Series in the Taebaeksan Basin in the northeastern part of the Korean Peninsula where limestone-hosted skarn deposits are common, this sedimentary basin is much less widely mineralized. It contains a thick succession of non-marine sedimentary rocks such as conglomerate, sandstone, mudstone, and minor carbonate and volcanoclastic rocks [1]. Carbonate-hosted skarn ore deposits have rarely been found in the Gyeongsang Basin. The Gukjeon deposit was reported as a unique manto-type skarn deposit which was formed by skarnization and skarn mineralization of the limestone beds interbedded within the Late Cretaceous Jeonggaksan Formation in the Gyeongsang Basin [1].
In this study, we attempted to identify an exact shape of orebodies and to construct a 3D model using geological information relating to the underground level, in order to more clearly represent the aspect of the mineralization caused in the limited horizons of the limestone formation. Zircon U-Pb geochronological and geochemical studies were carried out to interpret the spatial and temporal relationships among the geologic events, including the timing of mineralization.
This study is expected to provide the anatomical results for understanding the geological environment forming the Gukjeon Pb-Zn skarn deposit and for correlating the geological events of the Yucheon sub-basin (as well as the Gyeongsang Basin) in the southeastern part of the Korean Peninsula.

Regional Geology
The Gyeongsang Basin is tectonically subdivided into the Silla sub-basin, Yucheon sub-basin, and Yeongyang sub-basin from south to north ( Figure 1). The Nagdong trough is bounded to the west from these three sub-basins. All of these tectonic provinces were formed in a large pull-apart basin within the back-arc environment related to the subduction of the Pacific plate ( Figure 1) [7]. The study area, i.e., the Gukjeon Pb-Zn deposit, is located in the Miryang Block in the central part of the Yucheon sub-basin of the southern part of the Gyeongsang Basin ( Figure 1).
The Gyeongsang Basin comprises the mainly sedimentary Shindong Group, volcanogenic Hayang Group, and volcanic Yucheon Group ( Figure 1). About six hundred hydrothermal deposits, such as Au -Ag, Cu, Pb-Zn, and clay deposits, occur in this basin. The mineralization, showing progressive episodes from west to east within the Gyeongsang Basin, has been genetically associated with successive volcanism-plutonism -hydrothermal activities due to the subduction of the Pacific and Kula plates during the Late Cretaceous to the Early Tertiary (i.e., 100 to 60 Ma) [8].  [7][8][9]).

Figure 2.
Geologic map and sample locality of the study area (modified from [10,11]).
The major Pb-Zn orebodies in the study area are mainly hosted by the limestone member of the Jeonggaksan Formation, and minor orebodies are distributed in the dacite of a member of the Unmunsa subgroup to the east of the Gukjeon Pb-Zn deposit ( Figure  2). The deposit area is composed of hornfelsic siltstone and sandstone, andesitic tuff, andesite, skarns of the Jeonggaksan Formation, and the andesite and andesitic tuff of the Beomdori andesite. The mineralized interbedded limestones intermittently occur as several beds of 2-5 m in width and circa 50 m in length in the tuffaceous rocks of the Jeonggaksan Formation.

Miryang Andesite
Miryang andesite occurs throughout the southern part of the study area ( Figure 2) and is mainly composed of andesitic lava and porphyritic andesite. Aphanitic andesite lava showing a typical trachytic flow structure consists of small plagioclase laths and a small amount of epidote ( Figure S1A,B). These show gray, dark gray, and grayish-green colors and are carbonatized and epidotized, where epidotization is prominent near the pyrite veinlets ( Figure S1A). Porphyritic andesite is characterized by plagioclase and pyroxene phenocrysts, and the plagioclase shows a compositional zonal variation ( Figure  S1C,D).
Andesitic rocks are the main host rocks of the Miryang pyrophyllite deposit far from the study area (about 14 km to the west). This deposit has been reported to be a hydrothermal pyrophyllite deposit, which was formed by Late Cretaceous post-volcanic activities in the Gyeongsang Basin [12]. The K-Ar age of sericite occurring in the deposit has been reported as 69.7 ± 2.1 Ma [13].

Jeonggaksan Formation
The 500 m thick Jeonggaksan Formation is the most widely exposed formation in the eastern and northern parts of the study area ( Figure 2). This formation unconformably overlies the Miryang andesite ( Figure S2A) and is underlain by the Beomdori andesite. It is mainly composed of green, gray, and purple tuffaceous shale, sandstone, conglomerate, and interbedded limestone, as well as intercalated andesitic crystal-lithic tuff ( Figure S2C-F). Well-bedded tuffaceous sandstone and siltstone are the dominant components of this formation ( Figure S2C). Andesitic tuff is also a main component of this formation, which is interbedded with tuffaceous sedimentary rock. This rock is a crystal lithic andesitic tuff with abundant feldspar crystals and essential andesitic rock fragments ( Figure S2E,F). Some feldspars have been altered to epidote ( Figure S2F). The limestone of this formation, with a thickness of several meters, is the host rock of the Pb-Zn skarn deposit ( Figure S2B). Most limestones are skarnitized and the original composition is strongly changed, but the bedding is weakly preserved near the deposit.

Beomdori Andesite
The Beomdori andesite is distributed in the central part of the study area ( Figure 2). This rock, which overlies the Jeonggaksan Formation and has an aphyric texture and dark greenish color, is mainly composed of lath-like plagioclase, chlorite, calcite, epidote, and subordinate with accessory minerals [10]. This andesite is interpreted as the latest volcanic product of the Jusasan subgroup in the study area. The rock near the deposit consists of aphanitic andesite lava and porphyritic andesite. The porphyritic andesite is composed of dominant phenocrysts of plagioclase and hedenbergite, and rarely orthoclase, (ferro-) actinolite, and cryptocrystalline groundmass ( Figure S3A,B). Hedenbergite and (ferro-) actinolite phenocrysts have mostly been altered to chlorite.

Dacite
The dacite distributed in the northern part of the study area ( Figure 2) has a prominent porphyritic texture ( Figure S3C) with a cryptocrystalline groundmass. This rock is composed of dominant plagioclase phenocrysts, commonly of biotite, and rarely of orthoclase, and a microcrystalline groundmass of small plagioclase lath and quartz ( Figure  S3D).

Eonyang Granite (Bulguksa Intrusives)
The Eonyang granite distributed in the southern part of the study area ( Figure 2) belongs to the Bulguksa intrusives. This rock has a typical granophyric texture with an intergrowth of alkali feldspar and vermicular quartz ( Figure S3E,F). Small amounts of biotite and pyrite are associated with it. Most of the rocks in the study area belong to granophyre, which indicates the marginal facies of the stock of the Eonyang granite.
The underground map was taken from a previous study [3]. We could only access level 0 in the Gukjeon deposit, as the other levels (1-9) were flooded and inaccessible. For level 0, we obtained 3D cloud point data using a 3D laser scanning system (Riegl VZ-400i, RIEGL Laser Measurement Systems Gmbh, Riedenburgstraße, Austria) and reproduced the 3D underground mining model. The numerical data for the rest of the levels were obtained by digitizing from the 2D raster underground map. The final 3D geology and orebody geometry were made by merging the 3D scanned data for level 0 and digitized numerical data for levels 1-9.
The whole rock major and trace elements were analyzed on the molten bead by induced coupled plasma (ICP)-optical emission spectrometry and ICP-mass spectrometry at Activation Laboratory Ltd (Actlabs, Ancaster, ON, Canada) in Canada. Boron elements were analyzed by prompt-gamma neutron activation analysis using pulverized powder at Actlabs and the detection limit was 10 ppm.
Zircon separation was carried out according to the following procedure. Crushed rock powders smaller than 200 μm caused overflow of light mineral grains out of the beaker and magnetic minerals (e.g., magnetite, etc.) were removed by conventional magnetic separation, after which zircon grains were handpicked under a binocular microscope. The analytical points were decided after backscattered electron and cathodoluminescence (CL) imaging using a JEOL 6610LV (JEOL Ltd, Tokyo, Japan) scanning electron microscope at the Korea Basic Science Institute (KBSI). The SL13 (Sri Lankan Gem zircon) and FC1 (Duluth gabbroic anorthosite, 1099 Ma) [14] standard zircons were employed for U concentration (238 ppm) and U-Pb aging. The U-Th-Pb isotopic compositions were measured using the SHRIMP IIe/MC at KBSI. The analytical procedure of [15] was followed. Terra-Wasserburg diagrams and weighted mean age calculations were performed using Isoplot software (version 3.7-Berkeley Geochornology Center, CA, USA [16]). In this study, we used 208 Pb corrected U and Pb isotopic data for the Terra-Wasserburg diagrams and 207 Pb corrected 206 Pb/ 238 U ages for calculating the weighted mean ages.

Pre-Skarn Alteration
Contact metamorphic alteration before the skarnization strongly affected the tuffaceous and sedimentary rocks of the Jeonggaksan Formation. In particular, argillaceous and arenaceous sedimentary rocks were completely altered to hornfels. Most the hornfelsic rocks are distributed near the skarn orebodies of the Gukjeon Pb-Zn deposit. They are commonly found within the skarn zone of the underground and outcrops. The dense and hard hornfelsic rocks are mainly composed of small crystalline quartz, plagioclase, muscovite, chloritized biotite, and calcite ( Figure 3B). In addition, the matrix consists of very fine quartz and sericite. Despite the strong alteration, the original bedding and/or lamination of the host sedimentary rocks are well preserved in most areas ( Figure 3A).

Skarn Alteration
Skarn alteration of the Gukjeon Pb-Zn deposit is restricted to limestone interbedded with tuffaceous sedimentary rocks in the Jeonggaksan Formation. The representative skarn minerals are garnet, pyroxene, amphibole, axinite, and epidote. On the basis of optical microscope textural analyses, the skarn alteration can be divided into prograde and retrograde stages. The prograde stage is the metasomatic anhydrous stage. Garnet and pyroxene are formed in this stage and these commonly replaces the primary calcite and are replaced by the second stage retrograde phases, including amphibole, axinite and epidote, etc. during the following retrograde skarnization. The second stage is the retrograde hydrous stage. The most abundant skarn minerals of this stage are amphibole, axinite, and epidote. The detailed occurrences and geochemical characteristics of skarn minerals are as follows.

Pyroxene
Most pyroxenes occur as fragmented aggregates and discrete grains of small crystal relics in the skarn zone. They are replaced by anisotropic garnet ( Figure 3C) and large bladed axinite or replace large crystals of calcite ( Figure 3D). Judging from these replacement processes and occurrences, pyroxene seems to be the first prograde skarn mineral in the deposit. EPMA analytical results for pyroxenes at 19 points are shown in Table S1. Plotting to ternary diagrams (Mg-Ca-Fe 2+ ) for the analytical results of pyroxenes, they were mostly classified as hedenbergite ( Figure 4).

Garnet
Garnets are one of the most commonly occurring skarn minerals in the deposit. They are strongly fragmented and scattered as very small grains ( Figure 5B). Generally, garnet replaces calcite and is replaced by epidote ( Figure 5A,C,D) and amphibole. The garnet is associated with fragmented pyroxene aggregates and is abundant in the marginal zone of the Pb-Zn main orebody. Judging from these replacement processes and occurrences, garnet is the second prograde skarn mineral (after pyroxene) in the deposit. EPMA analytical results for garnet at 13 points are shown in Table S2. Except for one sample, all of the analyzed garnets have a high grossular proportion (XGrs) ranging from 78.45 to 83.90% and an andradite proportion (XAnd) ranging from 6.28% to 18.80%. The garnet grains from sample 19KJU-4-1 have higher andradite proportions (81.86% and 90.67%) than the other samples and lower grossular proportions (6.17% and 12.75%). As shown in Figure 6, whole analyzed garnets are plotted in the field of skarns/low-grade metabasic rock/calc-silicate granulite ( Figure 6).

Amphibole
Amphibole is the most abundant retrograde skarn mineral in the deposit. It generally occurs as long prismatic, small, flat, and/or fibrous aggregates ( Figure 7A,B). It overprints most of the primary minerals, such as calcite, and the prograde skarn minerals, such as pyroxene and garnet ( Figure 7B,C). They are sometimes injected by axinite veinlets. Judging from these replacement processes and occurrences, amphibole is the first retrograde skarn mineral in the deposit. EPMA analytical results for amphibole at 11 points are shown in Table S3. Plotting to the discrimination diagram of amphibole, they are classified as actinolite and ferro-actinolite ( Figure 8).

Axinite
Axinite is the most characteristic skarn mineral and occurs abundantly throughout the deposit. Axinites occur mostly as large crystals with a wedge or axe-head shape, or sometimes a fragmented irregular shape. They occur as cavity filling ( Figure 9A) and veinlet ( Figure 9B). Axinite replaces primary calcite and the retrograde skarn mineral, amphibole ( Figure 9D), and is replaced by small grains of epidote ( Figure 9B). Sometimes, axinite veinlets are injected into the amphibole skarn ( Figure 9C,D). Judging from these replacement processes and occurrences, axinite is the second retrograde skarn mineral (after amphibole) in the deposit. EPMA analytical results for axinite at 10 points are shown in Table  S4. Plotting to a discrimination diagram of axinite, they are classified as Mn-and Fe-rich tizenite and ferroaxinite ( Figure 10).

Epidote
Epidote is one of the most abundant skarn minerals in the deposit. It occurs as small to large lath forms and fills the matrix as small aggregates. It generally occurs as veinlets and has a layered form. Epidote commonly replaces the primary calcite ( Figure 11A) and cross-cuts pre-existing skarn minerals such as pyroxene, garnet, amphibole, and axinite ( Figure 11B-D). These replacement textures and occurrences indicate that the epidote is the final retrograde skarn mineral (after amphibole and axinite) in the deposit. EPMA analytical results for epidote at 32 points are shown in Table S5. Plotting to a discrimination diagram of epidote, they are classified as mostly epidote-clinozoisite ( Figure 12).

Post-Skarn Alteration
The post-skarn alterations of the deposit mostly consist of chloritization and silicification. Chloritization took place in favor of the fracture zone ( Figure 13A) and in the selective layering ( Figure 13B).

Orebody Geometry
The geology of the underground mining level is composed of hornfelsic rock, andesitic tuff, skarn rocks, and Pb-Zn orebody belonging to the Jeonggaksan Formation and Beomdori andesite ( Figure 14). The Beomdori andesite shows a fault contact with epidote skarn of the Jeonggaksan Formation. Most hornfelsic rocks preserve the original structures, such as lamination and bedding. The orebody is distributed concordantly with the bedding (NW strike and NE dip) of the Jeonggaksan Formation and is surrounded by hornfelsic rock and epidote-skarn ( Figure 14). The mineralization would have been caused within the limestone beds, selectively and limitedly, along the bedding and fracture zones. The shape of the orebody identified in the underground mining level of the deposit seems to be of a disk-or plate-type, with a convex form at the center. It was characterized as a geometric shape with a wide width in the south and a narrow edge in the north and it shows a NW strike and a gentle NE dip ( Figure 14). This disk-type orebody would be formed, limitedly and selectively, due to structural control throughout the vertical rise and lateral migration of the hydrothermal fluids derived from the concealed related intrusive.

Mineralization
Skarn ore minerals are mainly composed of sphalerite, chalcopyrite, galena, arsenopyrite, and pyrite. Ore minerals occur commonly as fracture filling and rarely disseminated minerals, and also in veinlets. Fracture-filling is the most common characteristic. The dominant mineral assemblages are sphalerite-chalcopyrite-galena (-arsenopyrite), sphalerite-galena-chalcopyrite-pyrite, and sphalerite-galena-pyrite (-arsenopyrite). Sphalerite, the major skarn mineral, occurs as a large irregular crystal or an aggregate and fills the fractures in the cataclastic matrix. It commonly contains small grains of chalcopyrite and irregular galena ( Figure 15A,B). Chalcopyrite characteristically occurs as an inclusion within sphalerite ( Figure 15A). The intergrowth of chalcopyrite blebs within sphalerite-termed "chalcopyrite disease", first introduced by [22]-is a common and characteristic occurrence in the deposit. Pyrite is the second dominant sulfide mineral, occurring as large euhedral single grains and small anhedral/subhedral aggregates. It has commonly replaced the sphalerite and filled the fractures of fragmented sphalerite ( Figure  15C). Pyrite has also been replaced by sphalerite ( Figure 15D). Therefore, pyrite seems to have formed before and after sphalerite. This indicates that the pyrite was precipitated throughout all skarn stages. Galena occurs as an inclusion within the sphalerite ( Figure  15A,B) and commonly replaced sphalerite ( Figure 15E,F). It also replaces pyrite ( Figure  15G). Arsenopyrite showing a wedge shape occurs rarely and contains sphalerite inclusions ( Figure 15H) and, therefore, seems to have formed after sphalerite. The paragenetic sequence of the major skarn sulfide minerals in the deposit is thought to be fine anhedral/subhedral pyrite, sphalerite, chalcopyrite, galena, and large euhedral pyrite. Considering the above-mentioned occurrences of the major sulfide minerals (i.e., sphalerite, chalcopyrite, galena, and pyrite) and the skarn minerals (i.e., pyroxene, garnet, amphibole, axinite, and epidote), the sulfide minerals were mostly precipitated in the amphibole and axinite skarns of the second retrograde stage.
The volcanic rocks of the study area exhibit a wide compositional range, from basaltic andesite, trachy-andesite, and andesite to dacite; in particular, the andesitic rocks of the Miryang andesite and Beomdori andesite show the same compositional range ( Figure  16A). The Eonyang granite shows SiO2 over-saturated characteristics, which indicates the latest composition of marginal facies of the granite stock, which was well-recognized by the micrographic texture ( Figure S3F). The volcanic rocks of the study area belong to the high-K calc-alkaline to calc-alkaline series, while the granitic rocks belong to the high-K calc-alkaline series ( Figure 16B). The analyzed andesite, granite, and dacite rocks are enriched in large-ion-lithophileelements (LILEs) and depleted in high-field-strength elements (HFSEs), although the negative troughs for Sr, P, and Ti are more prominent in the Eonyang granite ( Figure 17A). In the C1 chondrite normalized REE patterns ( Figure 17B), the Eonyang granite shows remarkable negative Eu anomalies (0.17-0.27) compared to those of Miryang andesite (0.53-0.99), Beomdori andesite (0.72-1.01), and dacite (0.87). The analyzed rocks are all very similar in terms of LREE enrichments (La/Lu)N: 5.1-11.3 for Miryang andesite; 4.4-6.9 for Beomdori andesite; 3.8-6.0 for Eonyang granite; 7.5 for dacite. In the tectonic discrimination diagrams (Figure 18), most of the analyzed rocks are plotted in the volcanic arc granite field. The Nb, Ta, and Ti negative anomalies in the trace element patterns ( Figure 17A) and tectonic discrimination diagrams ( Figure 18A,B) show that the Miryang and Beomdori andesites, Eonyang granite, and dacite in the study area were extruded or intruded in a subduction-related arc tectonic setting.

Boron, Lead, and Zinc Geochemistry
As mentioned above, axinite is the most important skarn mineral in the Gukjoen Pb-Zn deposit. Axinite in this deposit is closely related to the Pb-Zn ore minerals. Axinite can have a much higher boron (B2O3) content than other skarn minerals (pyroxene, garnet, amphibole, and epidote). Boron is a lithophile element, commonly occurring as a minor component in granite melt, while a total B2O3 concentration of about 1 wt.% can be reached in highly evolved peraluminous tourmaline-bearing granite [27]. It also can be easily transported during the fluid-rock interaction and has significance for tracing the fluid source in hydrothermal ore deposits [28].
To trace the source of the boron, we checked the boron contents of the fresh volcanic and igneous rocks neighboring the Gukjeon Pb-Zn deposit. We further analyzed the boron contents of Beomdori andesite-hosted argillic altered rocks in the southern and eastern sides of the Gukjeon deposit to compare the fluid sources of the skarn Pb-Zn mineralization in the Gukjeon deposit and the argillic alteration near the deposit.
The altered andesitic rocks from the Miryang and Jeonggaksan andesites have anomalous high boron contents (11.4 ppm and 82.6 ppm for Jeonggaksan andesite and 28.3 to 370 ppm for Miryang andesite). Argillic altered rocks including pyrophyllite and kaolin have moderately anomalous boron contents, ranging from 0.9 to 44.3 ppm (Table 1). Both fine-grained granodiorite and Eonyang granite have relatively low boron contents of 5.3 and 6.9 ppm for the fine-grained granodiorite and 0.6 to 5.4 ppm for the Eonyang granite (Table 1). Considering the boron contents, the altered andesitic rocks of the Jeonggaksan Formation and Miryang andesite, as well as the axinite skarn in the deposit, would be affected by boron-bearing fluid. Therefore, a boron-rich fluid probably affected the argillic alteration of the Beomdori andesite as well as the skarn alteration in the deposit. The granitoids such as fine-grained granodiorite and Eonyang granite can be excluded from the boron source due to their low boron concentrations.

Zircon SHRIMP U-Pb Geochronology
The representative CL images and SHRIMP U-Pb results are listed in Figure 19 and Table S7, respectively. Two andesitic lavas  and one andesite porphyry (19KJ-18-1) from the Miryang Andesite were selected for the SHRIMP zircon U-Pb analyses. Zircon grains from the andesitic lava  are equant to prismatic and euhedral crystals with oscillatory or banded zoning in CL images ( Figure 19). The U and Th contents range from 62 to 308 ppm and 56 to 280 ppm, respectively (Table S7) Two samples were selected from the Jeonggaksan Formation (20KJ2-1 and JK-1). Zircon grains from the tuffaceous sandstone (20KJ2-1) are generally prismatic crystals with oscillatory zoning (Figure 19). Their U and Th contents range from 55 to 5958 ppm and 20 to 2860 ppm, respectively. Their Th/U ratios vary from 0.25 to 0.92. One of these grains has an unzoned dark-CL emission of the inherited core, yielding an apparent age of 111. 4 Ma. Except for one datum, the analyzed data form one cluster and yield a weighted mean age of 90.7 ± 1.1 Ma (n = 23, MSWD = 2.6; Figure 20), which is interpreted as the maximum depositional age for the tuffaceous sandstone. Zircons from the andesite tuff (JK-1) are euhedral and stubby prisms with typical oscillatory zoning ( Figure 19). The weighted mean age from all the analyzed data yields 89.8 ± 2.0 Ma (n = 15, MSWD = 1.5; Figure 20), which is the eruption age of the andesite tuff.
We selected three porphyritic andesites (19KJ-22, JK-4, and JK-6) and one andesitic tuff (JK-3) from the Beomdori andesite for SHRIMP zircon U-Pb age dating. Zircon grains from the porphyritic andesite  are generally euhedral grains with banded or oscillatory zoning ( Figure 19). One data point (point 10.1) yields 155.1±3 Ma, interpreted as the inherited core of this porphyritic andesite. The rest of the analyzed data points yield a weighted mean age of 84.88 ± 0.94 Ma (n = 26, MSWD = 2.8; Figure 21). The U and Th concentrations range from 103 to 520 ppm and 61 to 633 ppm, respectively. Their Th/U ratios vary from 0.59 to 1.28. Zircon crystals from the porphyritic andesite (JK-6) are generally euhedral, prismatic grains showing a banded, combining sector, and oscillatory zoning ( Figure 19). The analyzed data points yield a weighted mean age of 83.8 ± 1.6 Ma (n = 14, MSWD = 1.3; Figure 21). Their U and Th contents and Th/U ratios range from 57 to 1348 ppm, 38 to 1421 ppm, and 0.51 to 1.31 ppm, respectively. Zircons from the porphyritic andesite (JK-4) are very similar to those of the previous two porphyritic andesites . All the analyzed data points make one cluster on the Terra-Wasserburg diagram and yield a weighted mean age of 83.8 ± 1.3 Ma (n = 16, MSWD = 1.03; Figure  21). The U content (81-436 ppm), Th content (49-546 ppm), and Th/U ratios (0.60-1.36) are also in a very similar range to those of samples 19KJ-22 and JK-6. Considering all analyzed data from the three porphyritic andesites of the Beomdori andesite (19JK-22, JK-6, and JK-4), the porphyritic andesites appear to have erupted between 88.88 and 83.8 Ma. Zircon crystals from the andesitic tuff (JK-3) are euhedral, prismatic grains generally showing combined oscillatory and sector zoning ( Figure 19). Their U and Th contents range from 85 to 654 ppm and 58 to 1075 ppm, respectively. Their Th/U ratios vary from 0.54 to 1.64. All analyzed data form one cluster in the Terra-Wasserburg diagram and yield a weighted mean age of 81.1 ± 1.1 Ma (n = 18, MSWD = 1.02; Figure 21), which was interpreted as the eruption age of the andesitic tuff.
Zircon grains from two samples (20KJ-16 & 20KJ-17) of fine-grained granodiorite show a similar CL texture. Those zircons are prismatic to hexagonal grains with an oscillatory or banded texture. The U and Th contents and Th/U ratios from sample 20KJ-16 range from 130 to 686 ppm, 96 to 1094 ppm, and 0.74 to 1.62, respectively, while those for sample 20KJ-17 vary between 85 and 831 ppm, 62 and 1564 ppm, and 0.58 and 1.88, respectively. The weighted mean ages for the two samples are 82.79 ± 1.80 Ma (n = 35, MSWD = 2.5) and 81.64 ± 0.41 Ma (n = 38, MSWD = 1.6). Together with the CL texture and Th/U ratios, fine-grained granodiorite intruded around 82 Ma ( Figure 21).
Zircons from the granitic rock of the Eonyang granite (19KJ-1012) are euhedral, prismatic grains with oscillatory zoning (Figure 19). All of the analyzed data form one cluster on the Terra-Wasserburg diagram and yield 72.0 ± 1.0 Ma (n = 22, MSWD = 1.5: Figure  21). Their U and Th contents and Th/U ratios range from 202 to 1158 ppm, 104 to 866 ppm, and 0.49 to 1.39, respectively. Considering their CL texture and Th/U ratios, Eonyang granite intruded at 72.0 Ma.
Zircon grains from dacite (19KJ-1001) are equant to prismatic crystals showing oscillatory, patch, and sector zoning ( Figure 19). Their weighted mean age is calculated as 68.62 ± 0.93 Ma (n = 15, MSWD = 0.97; Figure 21). Their U and Th contents and Th/U ratios vary from 37 to 940 ppm, 45 to 1822 ppm, and 0.49 to 2.29, respectively. Considering zircon texture and Th/U ratios, dacite erupted at 68.62 Ma.   We integrated the age dating results with the previous geochronological study [29]. The volcanism and sedimentation occurred as the successive stages or pulses, forming the Milyang andesite from 94.3 to 93.1 Ma, the Jeonggaksan Formation (andesite, andesitic tuff and tuffaceous sedimentary rocks) from 91.0 to 88.7 Ma, the Beomdori andesite from 84.9 to 78.4 Ma, and dacite and granite (Eonyang granite) from 72.0 to 68.6 Ma in the study area (Table 2 and Figure 22). The Jeonggaksan Formation, which was formed under submarine and subaerial environments, including carbonate beds, is the host formation of the Gukjeon Pb-Zn deposit ( Figure 22). Based on the age dating results and field occurrences, we suggest a simplified model showing the geological events including the mineralization ( Figure 22).

Discussion and Conclusions
Skarn formation was restricted to the limestone member of the Jeonggaksan Formation, where the major skarn minerals are garnet (mainly grossular with minor andradite), hedenbergite, amphibole (actinolite and ferro-actinolite), axinite (tizenite and ferroaxinite), and epidote (clinozoisite and epidote). The most characteristic skarn mineral was abundantly occurring axinite, which can be formed under a plentiful supply of boron derived from magma or a hydrothermal solution. Boron content is a good indicator for understanding the magmatic-hydrothermal evolutionary processes forming metallic deposits [30]. The occurrence of abundant axinite in the deposit indicates the continuous boron supply from the hydrothermal system forming the deposit. Skarn alteration in the deposit is characterized by prograde metasomatic pyroxene and garnet and retrograde metasomatic amphibole, axinite, and epidote. In the final stage, weak chloritization and silicification after the skarnization occurred in the deposit. The sequence of the skarn and alteration minerals in the deposit was interpreted as being in the order of pyroxene→ garnet→amphibole→axinite→epidote→chlorite throughout the metasomatic and hydrothermal processes in the system. The major skarn sulfide minerals of sphalerite, chalcopyrite, galena, and pyrite were mostly deposited during the retrograde stage which formed amphibole and axinite skarns. The major skarn sulfide minerals would have been formed in the order of small anhedral/subhedral pyrite, sphalerite, chalcopyrite, galena, and large euhedral pyrite. The chalcopyrite disease texture, with the inclusion of chalcopyrite blebs within sphalerite, can be explained by the selective dissolution of Fe-rich sphalerite and the addition of Cu (or both Cu and Fe) from fluid sources [31] or sphalerite-chalcopyrite co-precipitation by exsolution [32,33].
The Jeonggaksan Formation was formed from the Turonian to the Coniacian. The andesite (sample 090615-8, 78.4 ± 1.5 Ma, MSWD = 1.0, n = 3) from Beomdori andesite [29] is distinct from our results (porphyritic andesite, 84.9 ± 0.9, 83.8 ± 1.6, and 83.8 ± 1.3 Ma; andesitic tuff, 81.1 ± 1.1 Ma). In addition, the authors of [29] did not report detailed information for the analyzed Beomdori andesite such as petrography and geochemistry. Furthermore, the authors of [34] pointed out that the age dating results in [29] for Beomdori andesite were not matched to the lithological sequence of the neighboring area. We therefore suggest that the exact extrusion or intrusion age for the Beomdori andesitic rocks should be circa 84.9 to 83.8 Ma. The Eonyang granite intruded the Jusasan and Unmunsa subgroups and the emplacement age is 72.0 ± 1.1 Ma. This emplacement age is coincident with the crystallization ages ( [29]; 71.1 ± 0.8, 72.4 ± 0.6, 72.4 ± 0.7, and 71.8 ± 0.6 Ma) within the uncertainties. The final volcanic eruption recorded in the dacite (68.6 ± 0.9 Ma) in the study area occurs in the Maastrichtian.
As mentioned above, the occurrence of axinite, a major boron-bearing mineral, gives an important clue to understanding the skarn mineralization. The andesitic rocks from the Miryang andesite and Jeonggaksan Formation show geochemical affinities for B-Pb-Zn by the participation of boron-rich ore-bearing fluid in the system ( Figure 23). This probably indicates that the Pb-Zn mineralization in the deposit was coeval with formation of the axinite skarn. Boron is a major component of the ore-bearing fluid in the study area. Despite the absence of age dating results for the ore minerals, the timing of mineralization could be constrained by boron-rich hidden magmatic rocks during late Cretaceous (Campanian) (Figure 23). These results are different from the conclusions for the skarn mineralization associated with the Bulguksa granitoids (e.g., Eonyang granite) suggested by previous studies [6,7] and implicate hidden sources under the current surface level. The skarn alteration zone occurs in the limestone horizon in the Jeonggaksan Formation and the skarnization processes can be divided into prograde, retrograde, and final chloritization and silicification stages. The occurrence of axinite in the axinite skarn gives an important clue to the participation of boron-rich fluid in the skarn alteration and related mineralization. Concluding from the boron contents and geochronological results for the neighboring magmatic rocks, the ore minerals of the Gukjeon Pb-Zn deposit were probably precipitated by the infiltration and migration of boron-rich fluids derived from hidden pluton in the Campanian of the Late Cretaceous.
Supplementary Materials: The following are available online at www.mdpi.com/2075-163X/11/6/619/s1. Figure S1: Rock slabs (A and B) and microphotographs (B and D) of andesitic rocks from the Miryang andesite. A and B are andesite lava with small plagioclase lath and small epidotized pyroxene (Ep) phenocrysts, which show a trachytic texture. C and D are porphyritic andesites consisting of phenocrysts of zonal plagioclase (Pl), pyroxene (Px), and glassy groundmass.  Table S1: Chemical compositions of major and trace elements of igneous rocks near Gukjeon Pb-Zn deposit. Table S2: SHRIMP zircon U-Th-Pb isotope data of the study area.