Magmatic-Hydrothermal Processes of Vein-Type Haman-Gunbuk-Daejang Copper Deposits in the Gyeongnam Metallogenic Belt in South Korea

Haman, Gunbuk, and Daejang deposits are neighboring vein-type hydrothermal Cu deposits located in the SE part of the Korean Peninsula. These three deposits are formed by magmatic-hydrothermal activity associated with a series of Cretaceous granodioritic intrusions of the Jindong Granitoids, which have created a series of veins and alterations in a hornfelsed shale formation. The copper deposits have common veining and alteration features: 1) a pervasive chlorite-epidote alteration, cut by 2) Cu-Pb-Zn-bearing quartz veins with a tourmaline-biotite alteration, and 3) the latest barren calcite veins. Chalcopyrite, pyrite, and pyrrhotite are common ore minerals in the three deposits. Whereas magnetite is a dominant mineral in the Haman and Gunbuk deposits, no magnetite is present, but sphalerite and galena are abundant in the Daejang deposit. Ore-bearing quartz veins have three types of fluid inclusions: 1) liquid-rich, 2) vapor-rich, and 3) brine inclusions. Hydrothermal temperatures obtained from the brine inclusion assemblages are about 340–600, 250–500, and 320–460°C in the Haman, Gunbuk, and Daejang deposits, respectively. The maximum temperatures (from 460 to 600°C) recorded in the fluid inclusions of the three deposits are higher than those of the Cu ore precipitating temperature of typical porphyry-like deposits (from 300 to 400°C). Raman spectroscopy of vapor inclusions showed the presence of CO2 and CH4 in the three deposits, which indicates relatively reduced hydrothermal conditions as compared with typical porphyry deposits. The Rb/Sr ratios and Cs concentrations of brine inclusions suggest that the Daejang deposit was formed by a later and more fractionated magma than the Haman and Gunbuk deposits, and the Daejang deposit has lower Fe/Mn ratios in brine inclusions than the Haman and Gunbuk deposits, which indicates contrasting redox conditions in hydrothermal fluids possibly caused by an interaction with a hosting shale formation. In brines, concentrations of base metals do not change significantly with temperature, which suggests that significant ore mineralization precipitation is unlikely below current exposure levels, especially at the Haman deposit. Ore and alteration mineral petrography and fluid inclusions suggest that the Haman deposit was formed near the top of the deep intrusion center, whereas the Gunbuk deposit was formed at a shallower intrusion periphery. The Daejang deposit was formed later at a shallow depth by relatively fractionated magma.


INTRODUCTION
The Haman, Gunbuk, and Daejang deposits occur within the Gyeongnam belt in the southeastern part of the Korean Peninsula, which is a Cu-bearing province in South Korea, and the area has been well investigated for Cu ores (Sillitoe, 1980;Jin et al., 1982;Lee et al., 2003). The Gyeongnam belt hosts numerous hydrothermal vein-and breccia-type copper deposits (Koh et al., 2003;Choi et al., 2006) associated with a porphyritic granodiorite-granite-quartz monzonite resulting from late Cretaceous igneous activities (Bulguksa orogeny: 115-70 Ma) (Min et al., 1982). The Haman, Gunbuk, and Daejang deposits are vein-type Cu deposits in the Gyeongnam metallogenic belt (Choi, 1985;Shin et al., 1987) and are geographically neighboring and hosted in a Cretaceous shale formation (Haman-Jindong shale formations) and Cretaceous granitoids (Jindong Granitoids) (Park et al., 1985;Shin et al., 1987). These three deposits have been studied intensively because the area is regarded to have potential porphyry-style Cu mineralization (Choi, 1985;Park et al., 1985). The Haman and Gunbuk deposits consist of eight orebodies with an ore grade of 1-6 Cu wt% and a tonnage of 1,450 kt (Kim, 1973). The Daejang deposit is located NW of the Haman and Gunbuk deposits and consists of two orebodies and prospects. The ore grades and tonnages of the Daejang deposit have not been reported.
Numerous geologic and geochemical studies have been conducted on the ore deposits in the Haman, Gunbuk, and Daejang deposit areas. Park et al. (1985) reported high metal contents (e.g. Cu, Zn, and Pb) in the Jindong Granitoids as compared with other nearby Cretaceous granitoids. The wide range of salinities (6-68 wt% NaCl equiv.) determined for oremineralized fluids in the Haman and Gunbuk deposits (Park et al., 1985) was suggested to have been caused by hydrothermal system changes from an early lithostatic-magmatic fluid into a later meteoric-dominated hydrostatic fluid system (Choi, 1985). Such a hydrothermal evolution in the area was supported by oxygen and hydrogen isotope analysis . In the Daejang deposit, Shin et al. (1987) presented a range of pressures up to 4.0 kb, which was estimated by FeS analysis on sphalerite, and temperatures of 250-430°C based on AsS analysis on arsenopyrite. Park and Kim (1988) subdivided mineralizations in the Daejang deposit into stages; 1) the earliest base-metal sulfides, 2) Sb-Ag minerals, and 3) the latest barren calcite, and also estimated temperatures (140-280°C) and pressures (fs 2 10 -13 -10 -20 bar) of Sb-Ag mineral by fluid inclusion microthermometry. Based on the decreasing trends of the hydrogen-oxygen isotope values observed, Choi (1998) suggested the involvements of meteoric fluids in the later hydrothermal stage in the Daejang deposit.
Whereas individual "deposit-scale" geological features in each of the deposits have been studied, no integrated "district-scale" magmatic-hydrothermal model has been established for the area. Based on systematic field observations and sample petrography studies, we confirmed the alteration patterns and ore mineralogies of the three deposits. After petrography, we performed a fluid inclusion study, which included systematic petrographic observations, microthermometry, laser Raman spectroscopy, and LA-ICP-MS microanalysis. By studying the fluid inclusions in mineralized quartz veins, we obtained information on pressure, temperature, and salinity and on trace element contents of the ore-bearing fluids in the three deposits. Based on the information obtained, we constructed a geologic-geochemical process for magmatic-hydrothermal fluids and ore mineralization in the deposits. Implications regarding mineralization potentials are discussed.

Regional Geology
The Gyeongsang basin in the SE of the Korean Peninsula is a Cretaceous sedimentary basin formed in the Precambrian metamorphic basement (Figure 1). This basin is a tectonic pull-apart basin formed by northward oblique subduction of an ancient Izanagi oceanic plate into the Eurasian Plate during the Early Cretaceous. Subsequently, in the Late Cretaceous, the subduction angle became northwestward, and orthogonal subduction created an arc volcanism, called a Bulguksa orogeny (Chun and Chough, 1992;Chough et al., 2000;Lee et al., 2003;Ryu et al., 2006;Chough and Sohn, 2010). Bulguksa igneous rocks, generally granodioritic to granitic in composition, intrude into the Gyeongsang sedimentary basin. These are fine to coarse-grained porphyritic rocks mainly comprised of quartz, plagioclase, K-feldspar, biotite, and hornblende, and are classified as I-type magnetite series granitoid (Uchida et al., 2012). Ryu et al. (2006) divided the Gyeongnam metallogenic belt into five stratigraphic units based on integrated stratigraphy, paleontology, paleomagnetism, geophysics, and age dating of igneous activities. They suggested a process for the formation of the Gyeongsang basin involving igneous activity caused by the direction and angle of subduction of the Izanagi plate. Koh et al. (2003) divided the Gyeongnam metallogenic belt into 3-period groups according to mineralization ages and directions of ore veins and suggested the direction of mineralization in the Gyeongsang basin proceeded from west to east connecting with the three groups with a belt of igneous activities in the SE part of China, the SE part of the Korean Peninsula, and the SW part of Japan.
The Haman, Gunbuk, and Daejang deposits in the Gyeongnam metallogenic belt compose an Early to Late Cretaceous Sindong-Hayang sedimentary group, which is intruded by a series of Late Cretaceous Bulguksa magmas called Jindong granitoids ( Figure 1). The Haman and Jindong shale formations belong to the Sindong-Hayang sedimentary group and are hosted in the Haman, Gunbuk, and Daejang deposits. The Haman formation consists of conglomeratesandstone, black sandstone, and reddish mudstone (Chang, 1977;Chang et al., 1998). The Jindong formation, which conformably overlays the Haman formation, consists of black shale, black mudstone, and siltstone (Kim, 1973). The Haman-Jindong shale formations are tilted at a strike of N45 ∼ 70 E and a dip of 10 ∼ 15 SE (Moon et al., 1970). The Jindong granitoids intrude into both Haman and Jindong formations and are composed of porphyritic-plutonic textured rocks such as granite, granodiorite, and quartz diorite. Granodiorite is the most widely observed igneous rock among the Jindong granitoids. The whole-rock K-Ar (Park et al., 1985;Shin et al., 1987) and zircon U-Pb (Kim et al., 2016a) ages of the Jindong granitoids (granodiorite) have been reported to be 88-97 and 90-91 Ma, respectively. The Jindong granitoid intrusion in the Gyeongsang basin is a magnetite-bearing I-type intrusion, which has been suggested to have caused copper mineralization in the Gyengnam metallogenic belt, including the Haman, Gunbuk, and Daejang deposits (Park et al., 1983;Choi, 1985;Park et al., 1985;Shin et al., 1987;Lee, 1989;Koh et al., 2003). The Jingdong granitoids have lower Rb and Ba contents and higher Sr contents than the neighboring Pb-Zn related granitoids (Lee and Moon, 1989), which suggests fertile granitoids with higher Cu-Pb-Zn and Cl contents of biotite than barren granitoids . The shale in the Haman and Jindong formations in studied areas is widely metamorphosed into hornfels by intrusions.

Deposit Overview and Sampling Strategy
The Haman and Gunbuk deposits are neighboring deposits spatially located in a small district (about 25 km 2 ) in the SW part of the Gyeongnam metallogenic belt (35°13′ 07.9309″N, 128°21′ 41.7172″E ∼ 35°11′21.7″N, 128°23′15.1″E) ( Figure 1). They are hydrothermal vein-type Cu deposits with by-products of Au and Ag (Park et al., 1985;Heo et al., 2003). The estimated reserve of Cu ore in the Haman deposit is about 510 kt with an ore grade of 2-3 wt% of Cu (Kim, 1973). There are six orebodies in the Haman deposit, which are named Manse, Dundeok, Namgok-Bukgok, Deokgok, Gilgok, and Ogok (Figures 2A,B). The strike/ dip of ore veins in the Manse orebody are N15 ∼ 70°W to N80 ∼ 85°E and in the Ogok orebody are N20°∼ 30°E to 80°∼ 85°SE (Kim, 1973). The Gunbuk deposit is located 2 km to the NW of the Haman deposit. The estimated reserve of Cu ore in the Gunbuk deposit is about 960 kt (Kim, 1973). Joint and fault systems with a strike of N70°S to N15 ∼ 30°W are related to ore distribution in the Gunbuk deposit (Kim, 1973). There are two major orebodies in this deposit, named the Gunbuk and Baeknyun orebodies ( Figure 2C). The strike/slip of ore veins of the Gunbuk orebody is N5 ∼ 20 to N80°∼ 90°W, and ore grades are 5-10 wt% for Cu, 15-20 g/t for Au, and 70-300 g/t for Ag (Kim, 1973). The strike/dip of the veins in the Baeknyun orebody are NS, N15°to S40°∼ 80°W, and the ore grades are 1-2 wt% for Cu, 7 g/t of Au, and 100-260 g/t for Ag (Kim, 1973). The Daejang deposit is a hydrothermal vein-type Cu deposit with minor Au and Ag components (Shin et al., 1987;Choi, 1998) (35°14′ 05.9125″N, 128°18′ 05.4025″E) located 5 km NW of the Gunbuk deposit. There is no record of ore reserves for the Daejang deposit; joint and fault systems with a strike of NS to N20°W are related to ore distribution in this deposit (Kim et al., 2016b). The Daejang deposit has two orebodies, that is, North Cu-Pb-Zn-Au-Ag and South Cu-Au-Ag orebodies ( Figures  2D,E). The strike/slip of ore veins of the North orebody is N15°W to 65°∼ 85°NE, with ore grades of ∼0.11 wt% for Cu, 3-7 g/t for Pb, 2-5 g/t for Zn, ∼ 0.4 g/t for Au, and 380-770 g/t for Ag (Kim et al., 2016b). The strike/slip of ore veins of the South orebody is N20°W to 80°SW, and ore grades are ∼3.5 wt% for Cu, 1.4 g/t for Au, and 10-55 g/t for Ag (Kim et al., 2016b).
In the three deposits, we collected quartz vein samples showing a well-preserved texture of ore-alteration mineralogy and a crosscutting relationship. In the Haman deposit, we collected 31 quartz vein samples from the Manse and Ogok orebodies. In the Gunbuk deposit, we collected nine quartz vein samples from the Gunbuk and Baeknyun orebodies and some granodiorite outcrops. In the Daejang deposit, 10 quartz veins from the North and South orebodies, and some granodiorite outcrops were collected.

Ore and Alteration Minerals
Hornfels is widely formed (ca. 30°km 2 ) in the Haman and Jindong formations around Jindong granitoid intrusions ( Figure 1) and contain disseminated pyrite and magnetite. Magnetic susceptibilities (×10 −3 SI units) of the hornfels and granodiorite of Jindong granitoids in the region are about 50-60 and 17, respectively. The vein-bearing hornfels we studied commonly had pervasive chlorite and epidote alteration (CE alteration) before multiple veining events. In this section, we describe temporal relationships between hydrothermal alterations and paragenesis of ore-gangue mineral assemblages (Figures 3-6).
The major ore minerals of the Haman deposit include magnetite, pyrite, pyrrhotite, arsenopyrite, pyrite, and chalcopyrite, and minor ore minerals include scheelite, molybdenite, sphalerite, and galena. Quartz, tourmaline, actinolite, chlorite, and epidotes were observed as gangue minerals. We classified four veining stages in the Haman deposit: the first three quartz veining stages and the latest calcite veining stage. Hornfels in the Haman deposit has an early pervasive CE alteration. The hornfels with CE alteration was crosscut subsequently by a series of quartz veins hosting alterations of albite-actinolite (AA alteration) and tourmaline-biotite (TB alteration) ( Figure 3E). Quartz veins accompanying an AA alteration are barren. Quartz veins are curvy and have thicknesses of 1-3 cm ( Figure 3E). Quartz veins with TB alterations contain ore minerals and crosscut the earlier AA alteration ( Figures 3C,D). These veins have thicknesses of 2-3 cm, but some veinlets are less than 1 cm. We classified two separate stages among the quartz veins with TB alterations. Earlier ones have massive oxide and sulfide ore minerals of magnetite-tourmaline-chalcopyrite-pyrite with disseminated scheelite-molybdenite in the vein. Later ones accompany minor K-feldspar with TB alteration, which has similar mineralogy but is richer in sulfide minerals. The veins show tourmaline-pyrrhotite-chalcopyrite-pyrite-arsenopyrite with In the Gunbuk deposit, we observed magnetite, arsenopyrite, pyrrhotite, pyrite, and chalcopyrite as major ore minerals with sphalerite and galena as minor minerals. Quartz, tourmaline, and calcite were observed as gangue minerals. There are two veining stages in the Gunbuk deposit; a mineralized quartz vein and a postdating barren calcite vein. A pervasive CE alteration also appeared in hornfels before the quartz vein with TB alteration ( Figure 4B). Unlike the Haman deposit, AA alteration was absent in Gunbuk veins. The quartz veins are 1-2 cm thick and have ore minerals including massive magnetite, pyrrhotite, chalcopyrite, pyrite, arsenopyrite, sphalerite, and galena. Unlike the Haman deposit, the oxide-sulfide ores in the quartz veins occur together in the Gunbuk deposit. Sphalerite and galena were more commonly observed in the Gunbuk veins than in the Haman veins. Second stage calcite veins were barren and crosscut previous quartz veins ( Supplementary Table S1 and Figures 4, 6B).
The Daejang deposit, which is geographically removed from the Haman and Gunbuk deposits, contains a distinct ore mineral assemblage. They are Pb-Zn rich and contain no magnetite. Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 752908 5 Pyrrhotite, arsenopyrite, pyrite, chalcopyrite, sphalerite, and galena were observed as ore minerals, and quartz, and calcite as gangue minerals. Two stages of veins were identified in the Daejang deposit: earlier quartz veins with TB alterations and a later calcite vein stage. A pervasive CE alteration was also observed in hornfels, crosscut by a series of mineralized quartz veins with TB alteration. The quartz veins were 1-3 cm thick and hosted ore mineralization. In the Daejang deposits, TB alterations that developed along quartz veins were generally thinner (<1 cm) than in the Haman and Gunbuk deposits ( Figure 5B). The quartz veins have massive chalcopyrite-pyrite-sphalerite-galena with minor pyrrhotite-arsenopyrite ( Figure 5C). Unlike the Haman and Gunbuk deposits, the Daejang deposit hosts massive sphalerite and galena in quartz veins. It has been reported that some Sb-bearing minerals such as gudmundite and boulangerite are also present in the Daejan deposit (Shin et al., 1987;Park and Kim, 1988;Choi, 1998). The latest thin (<1 cm) and barren calcite veins crosscut former quartz veins (Supplementary Table S1 and Figure 6C).
A summary of the mineral and alteration paragenesis in the three deposits is provided in Figure 6D. The three deposits show common veining and alteration features: earlier CE alterations in hornfels, crosscutting quartz veins with TB alteration associated with oxidesulfide mineralization, and the latest calcite veins.
We observed varying distributions and proportions of fluid inclusion types in each deposit and orebody (Supplementary Table S2) Brine inclusions predominated in the Haman deposit (ca. > 50%), whereas liquid and vapor inclusions accounted for less than about 20 and 30%, respectively. In the Gunbuk deposit, brine inclusions accounted for ∼30%, liquid inclusions ∼40% and for ∼30% of vapor inclusions. In The chlorite-epidote altered host rock was overprinted by albite-actinolite alterations and cut by quartz veins with a tourmaline-magnetite-biotite alteration. The green dotted lines represent a chlorite-epidote "CE" alteration, the yellow dotted lines represent an albite-actinolite "AA" alteration, and the red dotted lines represent a tourmaline-(magnetite)-biotite "TB" alteration.

Fluid Inclusion Microthermometry
Fluid inclusion microthermometry was performed on doubly polished quartz chips from the 13 quartz vein samples obtained from each of the three deposits. A total of 133 FIAs (about 300 single fluid inclusions) were studied, that is, 50 from the Haman deposit, 41 from the Gunbuk deposit, and 42 from the Daejang deposit. Fluid inclusion microthermometry was carried out using a Linkam FTIR 600 heating-freezing stage at Inha University, which was calibrated using synthetic fluid inclusions of H 2 O and CO 2 -H 2 O. The freezing and heating procedures were conducted at accuracies of ±0.1 and ±1°C, respectively. Final ice melting and final halite melting temperatures were measured to calculate apparent salinities (NaCl equivalent wt%) (Bodnar and Vityk, 1994). Calculated salinities and homogenization temperatures (T h ) were combined to calculate a minimum trapping pressure and an isochore using a model NaCl-H 2 O system (Driesner and Heinrich, 2007). An average ±1 sigma was calculated for single fluid inclusions in the FIA (Supplementary Table S3).

Raman Spectroscopy
Laser Raman spectroscopy was performed on the bubble parts in vapor inclusions to identify gases in the fluids. A LabRAm HR Evolution confocal Raman spectrometer (HORIBA) in Inha University was used, which was equipped with an 1800 g mm −1 grating, a 100x confocal lens, and an edge filter. Raman spectra were obtained at a laser excitation wavelength of 532 nm (green laser), a laser power of 2 mW, filtering between 25 and 50%, a beam diameter of 100 μm, an acquisition time of 5 s, accumulation of 3-20 times, and a spectral range of 600-4,200 cm −1 . Twelve vapor inclusions from the Haman deposit, 21 from the Gunbuk deposit, and 13 from the Daejang deposit were analyzed.  Tourmaline and biotite were confirmed by microscopy. The green dotted lines represent a chlorite-epidote "CE" alteration, and the red dotted lines represent a tourmaline-(magnetite)-biotite "TB" alteration.
Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 752908 brine inclusions. The laser energy density used was 7-8 J/cm 2 and 3-10 Hz repetition rates were applied. To analyze inclusions and hosting quartz, we used 10 ms dwell times for the following isotopes; Sr, 133 Cs, and 208 Pb. A standard reference material (SRM) of NIST 610 glass was used as an external standard to calculate element ratios. Empirical estimation of Na concentrations from fluid inclusion salinities (Heinrich et al., 2003) was applied as an internal standard. Data reductions and quantifications from convoluted transient signals of fluid inclusions (Figure 8) were carried out using SILLS software (Guillong et al., 2008). During the ablations of hosting quartz, trace elements, including Ti and Al, in the quartz were simultaneously analyzed.

Microthermometry and Calculation of Pressure
Ranges of salinities and T h values obtained from liquid inclusion assemblages were found to be highly overlapping ( Figures 9A,  10). T h values from liquid inclusion assemblages were 96-360°C in the Haman deposit, 200-310°C in the Gunbuk deposit, and 150-290°C in the Daejang deposit. Salinities in liquid inclusion assemblages were 2-19 wt% in the Haman deposit, 3-18 wt% in the Gunbuk deposit, and 6-17 wt% in the Daejang deposit (Supplementary Table S3).
T h values of brine inclusion assemblages were 340-600°C in the Haman deposit, 250-470°C in the Gunbuk deposit, and 320-460°C in the Daejang deposit. Salinities in the brine inclusion assemblages are 33-62 wt% in the Haman deposit, 34-59 wt% in the Gunbuk deposit, and 39-55 wt% in the Daejang deposit. Fluid pressures of each deposit were calculated from salinities and T h values obtained from brine inclusions. Maximum pressures of the Haman, Gunbuk, and Daejang deposits were 500, 330, and 270 bar, respectively. T h values and calculated pressures were highest for the Haman deposit and similar for the Gunbuk and Daejang deposits ( Figure 9B).

Origin of Alterations and Ore Mineral Assemblages
Based on petrographic features such as vein-cross cutting, alteration, and ore mineral assemblages in the Haman, Gunbuk, and Daejang deposits, we classify their common paragenetic sequences as follows: 1) an early barren pervasive CE alteration overprinting hornfels, 2) main-stage quartz veins of the TB alteration mineralized with chalcopyrite-pyrite-arsenopyrite (±magnetite-pyrrhotite-sphaleritegalena), and 3) the latest barren calcite vein. The main-stage TB alteration was often associated with the AA alteration in the Haman deposit and with K-feldspar in the Haman, Gunbuk, and Daejang deposits.
The Haman-Jindong hornfels was formed by magmatism of the Jindong granitoids (Park et al., 1985;Shin et al., 1987). Magnetic susceptibilities (SI unit) in the hornfels were significantly higher (50-60 × 10 -3 ) than in granitoids (granodiorite; 15-17 × 10 -3 and pegmatite; 2-4 × 10 -3 ). Although we cannot observe under a microscope, we estimated that the hornfels contained disseminated microparticles of magnetite. Observations of disseminated pyrite in Haman-Jindong hornfels ( Figure 5A) and unmetamorphosed Haman-Jindong shale (Kim and Paik, 2001;Kim et al., 2018) indicated that the Fe from the Haman-Jindong parent shale was    (Klemm et al., 2007;Landtwing et al., 2010). The gray circles represent porphyry deposits, the white diamonds represent the Haman deposit, the gray triangles represent the Gunbuk deposit, and the black triangles represent the Daejang deposit. Ranges of the Gunbuk and Daejang deposits partially overlap with porphyry deposits, maximum temperatures were exceptionally higher in the Haman deposit than in porphyry deposits. (B) Plots of salinity versus calculated pressure for brine inclusion assemblages in the three deposits. The circles represent liquid inclusion assemblages, and the rhombus and triangle symbols represent brine inclusion assemblages. White circles and rhombus symbols represent T h and salinity of the Haman deposit, gray circles and triangles represent T h and salinity of the Gunbuk deposit, and black circles and triangles represent T h and salinity of the Daejang deposit.
Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 752908 probably the source of the formation of disseminated magnetite in hornfels. Whereas massive precipitation of magnetite was present in the quartz veins of the Haman and Gunbuk deposits ( Figure 3C and Figure 4A), and magnetite did not precipitate in the Daejang deposit, which may have been due to different temperature and sulfidation conditions of the magmatic-hydrothermal fluids (Einaudi et al., 2003) in the Haman, Gunbuk, and Daejang deposits. The quartz veins with the TB alteration in the Daejang deposit contained massive sulfide minerals such as pyrrhotite, pyrite, galena, and sphalerite, instead of magnetite ( Figure 5C), whereas Fe concentrations in the brines of all three deposits were comparable (Supplementary Table S5). These contrasting Fe mineral precipitations might be due to relatively lower hydrothermal temperatures, which agrees with our fluid inclusion results (Figures 9, 10), and suggests the possibility of higher sulfidation states of magmatic-hydrothermal fluids in the Daejang deposit as compared with the Haman and Gunduk deposits ( Figure 14) (Einaudi et al., 2003;Reed and Palandri, 2006). The CE alteration which occurs exclusively along the margin of the Jindong granitoids (Choi, 1998;Heo et al., 2003) indicates that the alteration is strongly associated with the intrusions. A CE alteration can be formed either by a convection of external fluids driven by magmatism (Seedorff et al., 2005;Ayuso et al., 2010;Sillitoe, 2010), or by a cooled magmatic fluids (Orovan et al., 2018;Pacey et al., 2020). The geochemical studies on the alteration minerals (e.g. chlorite, epidote) are required to trace the origin of the alteration and use it as an exploration target in the area .
An AA alteration ( Figure 3A) occurs exclusively in the Haman deposit. The AA alteration, which is associated with Na-Ca rich alteration minerals, is often observed in the barren and lower part of Cu ore bodies of porphyry deposits (Carten, 1986;Dilles and Einaudi, 1992;Dilles et al., 2000;Sillitoe, 2010) or in some Iron-Oxide-Copper-Gold (IOCG) deposits ( Barton and Johnson, 2004;Williams et al., 2005), in which replacing K-feldspar into Na plagioclase or epidote, and biotite or hornblende into Ca amphibole (actinolite) (Dilles and Einaudi, 1992). The existence of the AA alteration suggests that the Haman deposit may have been formed deeper than the Gunbuk and Daejang deposits, which is in line with our fluid inclusion results. The origins of AA alteration, which can be referred to as a Na-Ca alteration, occurring in some of the magmatic-hydrothermal ore deposits have been suggested to be either externally-derived salty fluids (Dilles and Einaudi, 1992;Barton and Johnson, 1996;Xavier et al., 2008) or magmatically-derived fluids (Oliver, 1996;Pollard et al., 1997a;Mark and Foster, 2000;Perring et al., 2000). In the Haman deposit, pervasive CE alterations occurred after the metamorphism of hornfels due to granitoid intrusions, and quartz veins with AA alterations cut through the CE alternations ( Figure 3A). The alteration and crosscutting relationship feature in the Haman deposit suggests that the AA alteration might have formed due to a magma-derived fluid from Jindong granitoids, but we cannot rule out the Occurrence of scheelite-molybdenite ( Figure 3D) in the Ogok orebody of the Haman deposit indicates an elevated W-Mo concentration in magmatic-hydrothermal fluids, which can be related to a significantly fractionated magma (Kooiman et al., 1986;Seedorff et al., 2005;Sinclair, 2007). This differentiated, possibly granitic, magma might be occur at depth, while outcrops of the Jindong granitoids show a rather intermediate magmatic composition ranging from gabbro to granodiorite (Wee et al., 2007).
Like biotite, tourmaline was a commonly observed alteration mineral in the mineralized veins of all deposits (TB alteration: Figures 3C, 4C, 5C), and the quartz veins associated with tourmaline contained boiling assemblages indicating fluid phase separation (Supplementary Table S3). Tourmaline is commonly associated with many alteration types such as Na-Ca, K, and phyllic alterations in porphyry deposits (Sillitoe and Sawkins, 1971;Gustafson and Hunt, 1975;Dilles and Einaudi, 1992). Lynch and Ortega (1997) suggested that the precipitation of tourmaline can be promoted by phase separation of magmatichydrothermal fluids rich in boric acid derived from a fractionated B-rich granite. B is abundant in oceanic sediments and in altered oceanic crust and has a relatively high solubility in fluids, which results in B-rich arc magma (Moran et al., 1992;Leeman et al., 1994;Dreyer et al., 2010). Abundant tourmaline in the deposits studied suggests a B-rich causative magma and B-rich magmatic fluids in the area. The Gukjeon Pb-Zn deposit, which is located 55 km NE from the Haman, Gunbuk, and Daejang deposits, contains abundant axinite [(Ca, Fe, Mn) 3 Al 2 BO 3 Si 4 O 12 (OH)] as a major skarn mineral, and its B has been suggested to originate from Late Cretaceous arc magma (Kim et al., 2021).

P-T Conditions and the Geochemistry of Hydrothermal Fluids
The presence of boiling assemblages suggests fluid phase separation due to a rapid depressurization, which commonly occurs in magmatic-hydrothermal deposits (Audétat et al., 1998;Heinrich et al., 1999). Many boiling assemblages have been reported in quartz veins in the three deposits, and these could provide the actual P-T conditions of hydrothermal fluids (Roedder and Bodnar, 1980). A density contrast between brine and vapor in boiling assemblages indicates degrees of depressurization during phase separation. Pressures calculated from brine analysis results of boiling assemblages can be converted into paleodepths based on lithostatic (density of 2.7 g/cm 3 ) or hydrostatic (density of 1.0 g/cm 3 ) pressure regimes. In the Haman deposit, the apparent paleodepths calculated from lithostatic and hydrostatic pressure regimes were 0.5-2.0 and 1.3-5.4 km, in the Gunbuk deposit were 0.1-1.2 and 0.2-3.2 km, and in the Daejang deposit were 0.3-1.0 and 0.8-2.7 km, respectively (Supplementary Table S3).
Since the Fe (Fe 2+ → Fe 3+ ) oxidation has a faster reaction rate than Mn (Mn 2+ → Mn 4+ ), Fe/Mn ratio can be used as a redox indicator of fluid passing through a fluid-rock buffered system (i.e. the Fe/Mn ratio is lower for an oxidizing fluid) (Boctor, 1985;Naeher et al., 2013). Bottrell and Yardley (1991) reported an increasing trend for Fe/Mn ratio when fluids at equilibrium with chlorite react with the reduced host rock, and Seo et al. (2012) showed a reducing fluid evolution from Cu-stage to Mo-stage in a porphyry deposit when Fe/Mn ratios in brines were increased. Fe/ Mn ratios in brine inclusions were found to be lower in the Daejang deposit than in the Haman and Gunbuk deposits ( Figure 13C), which suggests a relatively oxidizing fluid in the Daejang deposit. CO 2 was detected by Raman spectroscopy in some vapor inclusions (Figure 11), whereas clathrate was not observed by fluid inclusion microthermometry. Choi (1985) and Heo et al. (2003) reported that CO 2 was not detected in fluid inclusions of the Haman and Gunbuk deposits, whereas Choi (1998) reported a carbonic liquid phase in fluid inclusions of the Daejang deposit. CO 2 contents in fluids might have increased fluid carbonate activity in the fluids and formed the latest calcite vein (Choi, 1998). CH 4 was detected in fluid inclusions in the Gunbuk and the Daejang deposits by Raman spectroscopy (Figure 11). The ratios of mole fractions of CO 2 and CH 4 in fluid inclusions have been used to reconstruct fluid redox conditions (Giggenbach, 1987;Dubessy et al., 1989; Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 752908 14 Giggenbach, 1993). Relative redox conditions in deposits can be compared using ratios of their Raman peak heights. The CO 2 / CH 4 ratios of Raman peak heights from the Gunbuk and Daejang deposits were 4 and 9, respectively, which suggested more reducing fluid in the Gunbuk deposit. However, CO 2 /CH 4 ratios measurements are limited by vapor inclusions, due to a possibility of post-entrapment ion or H 2 diffusion (Li et al., 2009;Zajacz et al., 2009;Seo and Heinrich, 2013) of fluid inclusions in quartz, which could modify gas ratios. The fluids in the three deposits were reduced more than those in typical porphyry-type deposits, comparisons of Fe/Mn ratios in brine inclusions and CO 2 /CH 4 ratios in vapor inclusions in the three deposits suggested that fluids in the Daejang deposit were more oxidized than those in the Haman and Gunbuk deposits.
Since the three deposits are hosted within a shale formation (Figure 1), which is rich in organic carbon (Kim and Paik, 2001;Kim et al., 2018), differences between the redox states of fluids may have been caused by different degrees of interactions with the hosting shale. It was also noteworthy that, regardless of the redox conditions of the original magmatic ore fluids, fluid-rock interactions would effectively reduce fluids and precipitate ore minerals. Shale formation, therefore, in the study area might prove useful for exploring residual deposits. Rusk and Reed (2008) suggested hydrothermal fluid acidity controls Al concentrations in quartz, and thus, lower fluid pH values, and that this is highly pronounced in quartz formed at lower temperatures (<350°C) (Rusk and Reed, 2008). Similar Al concentrations in the quartz of Haman and Gunbuk deposits indicated similar pH conditions of ore-forming fluids ( Figure 12A). The Ti content in quartz has been applied as a Ti-in-quartz geothermometer (Wark and Watson, 2006;Huang and Audétat, 2012) in TiO 2 saturated or rutile-bearing system. Although this geothermometer cannot be applied directly because the activities of TiO 2 in hydrothermal systems may not have achieved equilibrium (Huang and Audétat, 2012), temperature conditions can be compared. The maximum Ti concentration in the Haman deposit was higher than in the Gunbuk and Daejang deposits, which suggests a higher temperature in the Haman deposit and concurs with our microthermometry results ( Figure 12B). If not in equilibrium, kinetically accumulated Al and Ti in hydrothermal quartz may represent the growth rate of the quartz in hydrothermal fluids (Mu€ller et al., 2002;Huang and Audétat, 2012), which suggests higher growth rates of quartz veins in the Haman deposit than in the Daejang deposit ( Figure 12).
The Rb/Sr ratios of magmatic-hydrothermal fluids are direct analogs of Rb/Sr ratios in causative magma. Since the Rb/Sr ratios in magma controlled by fractionation of K-feldspar and plagioclase, these ratios indicate a magma source or a degree of magma differentiation (Nockolds and Allen, 1953;Olade and Fletcher, 1975). Cs concentration in magmatic-hydrothermal fluids has been used to determine the fractionation degree of magma Klemm et al., 2008;Lüders et al., 2009). The ranges of Rb/Sr ratios and Cs concentrations observed in brine inclusion assemblages ( Figure 13B) suggest that magmatichydrothermal fluids that formed the Haman and Gunbuk deposits originated from magmas with similar degrees of differentiation, whereas those of the Daejang deposit indicate a more fractionated magma ( Figure 13B). Relatively higher base metal (Zn and Pb) concentrations in the fluids of the Daejang deposit ( Figures 13D,E) further support that this deposit may have been formed from a more differentiated magma than the Haman and Gunbuk deposits.

Geological Reconstruction of the Mineralization Processes
Concentrations of metals remain unchanged by changing the P-T conditions of fluids ( Figures 13D,E), and the temperature ranges of ore-forming fluids are relatively higher than that required for copper mineralization (<350°C) (Ayuso et al., 2010) in magmatic-hydrothermal deposits (Figure 9), indicating that a substantial ore at or below the current level of exposure may not be expected (Helgeson, 1969;Ulrich and Gǖnther, 2001;Landtwing et al., 2005). Ore metals in the fluids might have precipitated above the current exposure, but we would expect that most of the economic portions of such deposits might have already been eroded, especially at the Haman deposit where hightemperature alteration minerals such as actinolite and hornblende prevail. The presence of a fluid inclusion FIGURE 14 | The plot of hydrothermal temperatures versus sulfidation states of ore mineral assemblages, as described by Einaudi (2003) with some modification. Hydrothermal fluids of the Daejang deposit showed a higher sulfidation state (arsenopyrite-pyrrhotite-pyrite-chalcopyrite, no magnetite) and a lower temperature range than the Haman and Gunbuk deposits (arsenopyrite-pyrrhotite-pyrite-chalcopyrite, massive magnetite).
Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 752908 suggests that a deep (hydrostatically >5 km) hot (max. 600°C) magmatic-hydrothermal system exists in the Haman deposit that might be near the upper part of an intrusion center. However, the Gunbuk and Daejang deposits were formed shallower (hydrostatically <3 km) and under cooler conditions in Gunbuk and Daejang deposits (max. 500 and 460°C, respectively) than in the Haman deposit, possibly at the periphery or margin of intrusions. The fluids in the Haman and Gunbuk deposits were derived from magma of similar level of differentiation, while the fluids in the Daejang deposit were derived from a rather fractionated or physically separated magma. Since the three deposits are neighboring and possibly interacted with an organic-rich shale formation, a varying degree of fluid-rock interaction with the shale could accounts for the different redox states of fluids in the area (i.e. Fe/Mn ratios in Figure 13C).
Geological models for the area can therefore be established. 1) The hornfels in the Haman-Jindong shale were formed around the pluton of late Cretaceous Jindong granitoids ( Figure 15A). 2) A CE alteration in the hornfels was formed pervasively around the Jindong granitoids ( Figure 15B). 3) External fluids (or magmatic fluids) formed the quartz veins accompanying the AA alteration (only in the Haman deposit) ( Figure 15B). 4) Magmatic fluids formed the quartz veins accompanying TB alterations with Cu mineralization ( Figure 15C). The Haman deposit was formed near the upper part of a sub-volcanic magma stock, while the Gunbuk deposit was formed shallower at the periphery ( Figures 15D,E). In the case of the Daejang deposit, magmatic-hydrothermal fluid is more differentiated than in the Haman-Gunbuk deposits. In the Daejang deposit, the ore precipitated shallower and at the periphery of magma ( Figures 15D,E). The ore mineralogy in the Daejang deposit indicates a relatively higher sulfidation state than the magnetiterich Haman and Gunbuk deposits. 5) The deposits were exposed to the current surface level by denudation, and a portion of economic orebodies appeared to have been substantially eroded ( Figure 15F).

Comparison With Porphyry-Style Deposits and Exploration Potential in the Area
Rock-vein textures, alteration styles, and fluid inclusion types in the Haman, Gunbuk, and Daejang deposits were compared with the key characteristics of porphyry-style deposits. Geologic models of globally representative porphyry-style deposits, including an El Teniente Cu-Mo deposit (Klemm et al., 2007), Bingham Canyon Cu-Mo-Au deposit (Landtwing et al., 2005;Landtwing et al., 2010;Seo et al., 2012), and Ann Mason porphyry Cu deposit (Dilles and Einaudi, 1992), were applied to determine the possibility of a hidden economic ore in the studied area. The CE alteration in the study area corresponds to propylitic alteration of porphyry deposits. The CE alteration was pervasively distributed but only occurs in the hornfels of the Haman-Jindong shale. The CE alteration in hornfels can thus be used to vector the magmatic-hydrothermal or mineralization Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 752908 center. Economic portions of porphyry deposits are typically associated with K-feldspar-biotite (potassic) and quartz-sericite (phyllic) alterations in deposits (Sillitoe, 2010). Main stage orebearing veins in the Haman, Gunbuk, and Daejang deposits are associated with the TB alteration (similar with the potassic alteration). Phyllic alteration was absent in studied deposits possibly because of large formation depths. Our fluid inclusion results show that the TB alteration in the studied deposits formed at a higher temperature than the temperature range of typical potassic alterations in porphyry deposits. T h and pressure in brine inclusion assemblages in the El Teniente deposit are 320-410°C and max. 300 bar (Klemm et al., 2007), and about 350-410°C and max. 300 bar in the Bingham Canyon deposit (Landtwing et al., 2010). These P-T conditions are somewhat lower than those of the Haman deposit (340-600°C and max. 500 bar) but overlap with those of the Gunbuk (250-500°C and max. 330 bar) and Daejang deposits (320-460°C and max. 270 bar) ( Figure 9A). The AA alteration that appears in the Haman deposit can be compared with the sodic-calcic alteration (Na-Ca alteration) of some porphyry deposits (Dilles and Einaudi, 1992;Sillitoe, 2010). P-T conditions and alteration styles suggest that the Haman deposit might represent a deep and barren part of a porphyry-style deposit. Although the maximum temperatures are generally higher than the range of porphyry deposits, the Gunbuk and Daejang deposits are closer to the "economic" temperature range.
The alteration patterns and range of fluid temperatures observed suggest little possibility of an economic ore in the Haman deposit, but a possibility in the area of the Gunbuk and Daejang deposits. More studies of petrology and geochemistry such as on zircon fertility and stable isotopes in the area are required to fully understand its oreforming processes and explore hidden economic orebodies.

CONCLUSION
Based on the petrography of the vein-bearing samples and the fluid inclusion studies conducted on Haman, Gunbuk, and Daejang deposits, we derive the following conclusions.
1) In the three deposits, pervasive CE alteration overprinting hornfels, TB alteration with ore mineralized quartz veins, and the latest barren calcite veins commonly occur. 2) The hydrothermal fluids of the Haman and Gunbuk deposits have relatively lower sulfidation and oxidation states than the Daejang deposit.
3) The Daejang deposit appears to have been formed from a relatively later differentiation stage magma than the Haman and Gunbuk deposits, which could result in distinctive ore mineral differences between the Daejang and Haman-Gunbuk deposits. 4) The Haman deposit was formed at the central-upper part of crystalizing magmatic stock under P-T conditions higher than those typical of economic porphyry Cu deposits. The Gunbuk and Daejang deposits were formed at a relatively shallow peripheral location and partially overlapped with the P-T conditions of porphyry deposits. 5) Because of the high fluid pressure and temperature, no largescale Cu mineralization could be expected underneath the current exposure, especially in the Haman deposit. If there is an orebody remaining underneath, it might be placed in the vicinity of the Gunbuk and Daejang deposits.

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 author.

AUTHOR CONTRIBUTIONS
TL and JS wrote the manuscript. TL, JS, BY, BL, and JL did fieldwork in the area. TL performed petrography of samples and fluid inclusion microthermometry. TL, JS, YY, and SH performed an LA-ICP-MS microanalysis. TL performed laser Raman spectroscopy. All authors discuss the results and correct the manuscript.