# Auriferous pyrite formed by episodic fluid inputs in the Akeshi and Kasuga high-sulfidation deposits, Southern Kyushu, Japan

Mizuki Ishida, Rurik Romero, Mathieu Leisen, Kazutaka Yasukawa, Kentaro

Nakamura, Fernando Barra, Martin Reich, Yasuhiro Kato

► To cite this version:

Mizuki Ishida, Rurik Romero, Mathieu Leisen, Kazutaka Yasukawa, Kentaro Nakamura, et al. . Auriferous pyrite formed by episodic fluid inputs in the Akeshi and Kasuga high-sulfidation deposits, Southern Kyushu, Japan. Mineralium Deposita, 2021, 57, pp.129-145. 10.1007/s00126-021-01053-4 .
hal-04902908

HAL Id: hal-04902908

https://ut3-toulouseinp.hal.science/hal-04902908v1

Submitted on 21 Jan 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Auriferous pyrite formed by episodic fluid inputs in the Akeshi and Kasuga high-sulfidation deposits, Southern Kyushu, Japan

Mizuki Ishida<fnref n="1" /> ID • Rurik Romero<fnref n="2" />,<fnref n="3" /> ID • Mathieu Leisen<fnref n="2" />,<fnref n="3" />,<fnref n="4" /> ID • Kazutaka Yasukawa<fnref n="1" />,<fnref n="5" />,<fnref n="6" /> ID • Kentaro Nakamura<fnref n="1" /> ID • Fernando Barra<fnref n="2" />,<fnref n="3" /> ID • Martin Reich<fnref n="2" />,<fnref n="3" /> ID • Yasuhiro Kato<fnref n="1" />,<fnref n="5" />,<fnref n="6" />,<fnref n="7" /> ID

Received: 30 October 2020 / Accepted: 15 March 2021 / Published online: 31 March 2021

© The Author(s) 2021

## Abstract

Pyrite geochemistry has proven useful for tracking changes in the composition and physico-chemical conditions of hydrothermal fluids in ore-forming environments. Here, we investigated the microtextural features and chemical composition of pyrite, a main Au-bearing phase in the Akeshi and Kasuga deposits (Southern Kyushu, Japan), to better constrain the ore-forming processes in these high-sulfidation epithermal Au deposits. Despite the widespread distribution of Au-bearing pyrite in both deposits, no visible Au minerals coexist with pyrite. However, in situ laser ablation inductively coupled plasma mass spectrometry results show that Au concentrations in pyrite vary from below the detection limit to 41 ppm and are positively correlated with Cu ($r = 0.4$; up to 7400 ppm) and Bi concentrations ($r = 0.44$; up to 640 ppm). In both deposits, high Cu and Au concentrations occur in small ($< 25 \ \mu\text{m}$) anhedral grains of pyrite, which are interpreted to have rapidly crystallized from the ore-forming hydrothermal fluid. In addition, dissolution-reprecipitation textures and thin, concentric, Cu-rich overgrowths were identified in a number of larger ($> 25 \ \mu\text{m}$) pyrite grains and aggregates. These abrupt changes in the trace element compositions of pyrite grains likely record episodic metal-rich fluid inputs. We also propose that gold adsorption onto growing pyrite surfaces played a key role in the mineralization of these deposits.

**Keywords** Epithermal Au deposits • Pyrite geochemistry • LA-ICP-MS • EPMA • Nansatsu district

## Introduction

Southern Kyushu, Japan, is known for its intense volcanic-hydrothermal activities related to the subduction of the Philippine Sea plate over the past 5 million years (Watanabe 2005; Mahony et al. 2011). Abundant epithermal Au

mineralization is associated with this volcanism, including all four currently mined epithermal Au deposits in Japan (Watanabe 2005; Garwin et al. 2005). Active hydrothermal systems, including Satsuma-Iwojima (Shinohara et al. 1993; Hedenquist et al. 1994a) and Kirishima (Hedenquist and Aoki 1991), have been studied recently, plus the Quaternary

Editorial handling: B. Lehmann

<figure><img src="image_2.png" /></figure>

Yasuhiro Kato
ykato@sys.t.u-tokyo.ac.jp

<footnote n="1">¹ Department of Systems Innovation, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan</footnote>

<footnote n="2">² Department of Geology and Andean Geothermal Center of Excellence (CEGA), FCFM, Universidad de Chile, Plaza Ercilla, 803 Santiago, Chile</footnote>

<footnote n="3">³ Millennium Nucleus for Metal Tracing Along Subduction, FCFM, Universidad de Chile, Plaza Ercilla, 803 Santiago, Chile</footnote>

<footnote n="4">⁴ Present address: Géosciences Environnement Toulouse(GET), Université de Toulouse, UPS, CNRS, IRD, CNES, 14 avenue E. Belin, 31400 Toulouse, France</footnote>

<footnote n="5">⁵ Frontier Research Center for Energy and Resources (FRCER), The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan</footnote>

<footnote n="6">⁶ Ocean Resources Research Center for Next Generation (ORCeNG), Chiba Institute of Technology, 2-17-1 Tsudanuma, Narashino, Chiba 275-0016, Japan</footnote>

<footnote n="7">⁷ Submarine Resources Research Center, Research Institute for Marine Resources Utilization, Japan Agency for Marine-Earth Science and Technology (JAMSTEC), 2-15 Natsushima-cho, Yokosuka, Kanagawa 237-0061, Japan</footnote>
Hishikari low-sulfidation Au deposit (e.g., Izawa et al. 1990; Faure et al. 2002; Morishita et al. 2018; Fig. 1). Other studies have focused on the Pliocene high-sulfidation Au deposits in the Nansatsu district, southern-most Kyushu. Earlier studies have reported ubiquitous pyrite in these deposits (Tokunaga 1954; Nakamura et al. 1994; Hedenquist et al. 1994b) as an important host of primary Au (Imai et al. 1996; Morishita et al. 2019). However, there has been few studies on the microtextural characteristics and the trace element composition of pyrite in these deposits.

Detailed microtextural observations of pyrite combined with in situ geochemical analyses have proven useful in unraveling Au mineralization processes in ore-forming environments and, particularly, for tracking changes in the composition and physico-chemical conditions of hydrothermal fluids (e.g., Deditius et al. 2008; Reich et al. 2013; Tardani et al. 2017; Román et al. 2019). This is because pyrite (i) is ubiquitous and stable over a wide range of physico-chemical conditions in ore deposits (Deditius et al. 2008, 2014; Reich et al. 2013) and (ii) can incorporate a variety of metals and

semi-metals (Cu, Co, Pb, Sb, As, Ag, Au, Ni, Zn, Se, Te, and Hg, among others) either in solid solution or as micro- to nano-sized inclusions depending on its formation conditions (Deditius et al. 2011, 2014; Román et al. 2019). Moreover, pyrite has been recognized as an important host of Au in ore deposits because pyrite can absorb and incorporate Au from fluids undersaturated with respect to native Au, increasing its scavenging potential in ore systems (Reich et al. 2005; Cook et al. 2009; Gopon et al. 2019; Pokrovski et al. 2019; Kusebauch et al. 2019).

Here, we present the microtextural occurrence and in situ chemical analyses of pyrite in the Akeshi and Kasuga high-sulfidation epithermal Au deposits. Microtextural observations using scanning electron microscopy (SEM) were integrated with major, minor, and trace element analyses of pyrite performed by electron probe microanalysis (EPMA) and laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS). By combining these data, we contribute additional information on the hydrothermal activities that led to Au mineralization in the Nansatsu district.

<figure><img src="image_3.png" /><figcaption>**Fig. 1** Location of the Akeshi and Kasuga deposits, Kagoshima, Japan</figcaption></figure>

<figure><img src="image_4.png" /></figure>
Deposit geology and samples

There are several Au-bearing silicified rocks and high-
sulfidation epithermal Au deposits in the southern part of the
Satsuma peninsula in southern Kyushu, Japan (Nansatsu
district; Fig. 2). These deposits, including Akeshi and
Kasuga, are referred to as the Nansatsu-type deposits. The
district is also considered to have potential for porphyry cop-
per mineralization, although no deposit has been discovered to
date (Izawa and Hayashi 2018; Watanabe et al. 2018). Alunite
K-Ar ages and paleomagnetic data from previous studies in-
dicate an approximate age of 4 ± 1 Ma for the formation of the
Nansatsu-type deposits (Izawa et al. 1984; Takeda et al. 2001;
Watanabe 2005).

Gold mineralization in the Nansatsu district was closely
associated with the eastward-shifting calc-alkaline volcanism
in this region (Izawa 1992; Watanabe 2005). The Kasuga
deposit is hosted by the Upper Formation of Nansatsu
Group, which comprises pyroxene andesite lavas and pyro-
clastic rocks, whereas the Akeshi deposit is hosted by the

Nansatsu Middle Volcanics, composed of hornblende andes-
ite lavas and pyroclastic rocks (Nakamura et al. 1994; Metal
Mining Agency of Japan 1993).

Total gold production in the Akeshi deposit has reached 12.7 metric tonnes (t) since 1912 (average Au grade of 5.8 g/t; H. Kon per. commun. 2020) and 10.7 t since 1910 in the Kasuga deposit (average Au grade of 2.7 g/t; K. Shinagawa per. commun. 2020). Hydrothermal activity contemporaneous with volcanism resulted in strong leaching of the volcanic host rocks by acidic fluids (Izawa 1992; Hedenquist et al. 1994b). The preferential dissolution of phenocrysts and pumice fragments produced a vuggy texture in residual quartz (Hedenquist et al. 1994b). Secondary deposition of quartz is common in vugs and veinlets, locally resulting in massive siliceous rock (Hedenquist et al. 1994b). The ore zone, with gold grades of 1 g/t or more, is restricted to mushroom-shaped residual quartz bodies (Online Resource 1, Figs. S1 and S2) on the order of 100 m wide, not including the advanced argillic alteration zone around the orebody (Tokunaga 1954). Published works provided detailed descriptions of the hydro-

<figure><img src="image_5.png" /><figcaption>**Fig. 2** Geological map of the Nansatsu district. Modified from the Ministry of International Trade and Industry (1985)</figcaption></figure>
thermal alteration and gangue minerals associated with the orebody (Tokunaga 1954; Nakamura et al. 1994; Hedenquist et al. 1994b).

Limited fluid inclusion studies have been conducted in the Akeshi and Kasuga deposits due to the small size of quartz grains. All these analyses were performed on secondary quartz, formed after the initial leaching event (Hedenquist et al. 1994b; Nakamura et al. 1994). Homogenization temperatures range from 160 to 350 °C, and the occurrence of vapor-rich inclusions suggests that fluids were locally trapped from boiling solutions (Takenouchi 1983; Nakamura et al. 1994; Hedenquist et al. 1994b). Generally, the salinities are low (< 1 wt% NaCl<sub>eq</sub>) in both deposits, but higher salinities—based on the presence of daughter minerals—were reported in one sample from the Kasuga deposit (Hedenquist et al. 1994b).

The Akeshi and Kasuga ore zones typically comprise around 90 wt% SiO₂ as fine-grained quartz (Nakamura et al. 1994; Fig. 3). The pyrite content of the ore is reported to increase with depth, partly reflecting a decrease in post-depositional oxidation (Hedenquist et al. 1994b). Today, pyrite is ubiquitous throughout the orebody, possibly because the upper part of the orebody has already been completely mined. Alunite and native sulfur are also common in the orebody (Hedenquist et al. 1994b; Nakamura et al. 1994). According to Imai et al. (1996), the Akeshi and Kasuga ores can be classified into three types: (i) unoxidized pyrite–gold ore, (ii) unoxidized gold–copper ore, and (iii) oxidized high-grade gold ore. The unoxidized pyrite–gold ore type is the most abundant, and is composed of pyrite and rutile dissemi-nated in fine-grained quartz. The unoxidized gold–copper ore, which comprises pyrite, enargite (in Kasuga) or luzonite (in Akeshi), covellite, scorodite, and rare native gold within a fine-grained quartz matrix, is found mainly in the upper sec-tion of the Akeshi deposit (Tokunaga 1954; Nakamura et al. 1994). The relation between these first two ore types is not well constrained. The oxidized high-grade Au ore is related to white quartz veinlets with native gold and rust-colored silici-fied rock with iron oxides (goethite, scorodite, and tripuhyte), barite, clay, and native gold (Nakamura et al. 1994; Imai et al. 1996). This ore type has the highest Au grades locally reaching 1000 g/t or more (Imai et al. 1996) and is likely related to hypogene and/or supergene oxidation and weathering of unoxidized ore types.

Samples analyzed in this study were collected from the residual quartz orebodies in both deposits. Nine samples were obtained from three different locations in the open pit in the Akeshi deposit, and seven samples were selected from two different areas in the open pit in the Kasuga deposit (Table 1). Sampling locations, which were the active mining sites at the time of sample collection, are shown in Online Resource 1 (Figs. S1 and S2). Pyrite-rich

<figure><img src="image_6.png" /><figcaption>**Fig. 3** Photographs of representative ore samples (residual quartz rock). **a** Gray-colored area is pyrite-rich, whereas the reddish zone corresponds to iron oxides. White area consists of fine-grained quartz. Akeshi deposit. **b** Detailed photo of **a** showing native sulfur in the pyrite-rich zone, fine-grained quartz, and iron oxides in the upper right corner. **c** Pyrite-rich and iron oxides zone in hand sample from the Kasuga deposit</figcaption></figure>

parts of each sample (Fig. 3) was selected to avoid the effect of supergene oxidation and prepared as a 2-cm-diameter polished sections. Analytical methods are reported in the **Appendix**.

<figure><img src="image_7.png" /></figure>
# Results

## Scanning electron microscope observations

Pyrite is the main sulfide mineral in all samples analyzed (Fig. 4, Table 1). In some samples, pyrite is spatially associated with enargite (Fig. 4a, Table 1) or rutile (Fig. 4b; Table 1), indicating that the studied samples are dominated by unoxidized pyrite–gold ore (Imai et al. 1996). Pyrite has various morphologies, from euhedral crystals (Fig. 4b, d, f) to anhedral aggregates (Fig. 4c, e). Their grain sizes also varied significantly, from a few microns to up to hundreds of microns (Fig. 4). Since pyrite grains are finely dissemi- nated within the vuggy quartz matrix, it was not possible to discriminate pyrite associated with the initial leaching stage from pyrite formed during secondary quartz deposition. Textural features observed in eight out of 16 samples (Table 1, Fig. 5a–c, Online Resource 1 Fig. S3) indicate coupled dissolution–reprecipitation (CDR) processes (i.e., fluid-mediated mineral replacement reactions; Putnis 2009). Two common characteristics of CDR were recog- nized in the studied pyrite grains: (i) a porous phase that preserves the external morphology of the original grain and (ii) sharp reaction fronts between the primary grain and product phases (Fig. 5c). Some grains have multiple reaction fronts that are parallel to each other (Fig. 5a, b).

## Electron probe microanalysis and X-ray maps

EPMA measurements were carried out on primary, inclusion- and pore-free pyrite crystals. All EPMA analyses are reported in Online Resource 2. Cobalt and Bi data are not discussed here, because the measured concentrations may have been affected by spectral interferences from Fe and S, respectively. Overall, trace element concentrations in pyrite as detected by EPMA range from a few wt% to detection limits (~100–900 ppm). Copper, Se, and Sb were detected in more than a quarter of the total analyzed points. For other elements, concentra- tions were below detection in more than three quarters of the analyzed points. Copper showed the greatest variation (medi- an: 0.06 wt%), ranging from a maximum concentration of 3.9 wt% to below the detection limit (251 ppm). Selenium also showed variation (median: 0.06 wt%), ranging from 0.6 wt% to below the detection limit (216 ppm); higher Se concentrations were mainly attributed to two samples from the Kasuga deposits (samples KA15 and KA19). The maximum concentration of As, a key element commonly contained in high abundance in pyrite from a large variety of Au deposits, was 0.4 wt% (median: 0.05 wt%).

WDS X-ray maps for selected elements in pyrite grains are shown in Online Resource 3. Due to the relatively high detec- tion limits of the EPMA technique, no clear zonation was observed in the WDS maps for most elements. However, Cu

<table><caption>**Table 1** List of the polished section used in this study, the summary of pyrite morphology, and mineral assemblage</caption><thead><tr><th rowspan="2">Sample name</th><th rowspan="2">Deposit</th><th rowspan="2">Orebody</th><th rowspan="2">Depth (mL)</th><th colspan="3">Pyrite</th><th rowspan="2">Rutile</th><th rowspan="2">Enargite</th><th rowspan="2">CDR</th></tr><tr><th>Euhedral</th><th>Subhedral</th><th>Anhedral</th></tr></thead><tbody><tr><td>AK01 01-2</td><td>Akeshi</td><td>#1</td><td>15</td><td>+</td><td></td><td></td><td></td><td></td><td>+</td></tr><tr><td>AK01-9</td><td>Akeshi</td><td>#1</td><td>15</td><td>+</td><td>+</td><td></td><td></td><td></td><td>+</td></tr><tr><td>AK02-03</td><td>Akeshi</td><td>#2</td><td>5</td><td>+</td><td>+</td><td>+++</td><td></td><td>+</td><td>+</td></tr><tr><td>AK02-06</td><td>Akeshi</td><td>#2</td><td>5</td><td>+</td><td>+</td><td></td><td></td><td></td><td>+</td></tr><tr><td>AK04-0</td><td>Akeshi</td><td>#4</td><td>20</td><td>+</td><td></td><td></td><td></td><td></td><td>+</td></tr><tr><td>AK04-04</td><td>Akeshi</td><td>#4</td><td>15</td><td>+++</td><td>+</td><td>+</td><td>+</td><td></td><td>+</td></tr><tr><td>AK15112705</td><td>Akeshi</td><td>#4</td><td>10</td><td>+</td><td></td><td></td><td>+</td><td></td><td>+</td></tr><tr><td>AK15112706</td><td>Akeshi</td><td>#4</td><td>10</td><td>+++</td><td></td><td>+</td><td></td><td></td><td>+</td></tr><tr><td>AK15112710</td><td>Akeshi</td><td>#4</td><td>20</td><td>+</td><td></td><td></td><td></td><td></td><td>+</td></tr><tr><td>KA05</td><td>Kasuga</td><td>Honko</td><td>75</td><td></td><td>+</td><td>+</td><td></td><td></td><td>+</td></tr><tr><td>KA08</td><td>Kasuga</td><td>Honko</td><td>75</td><td>+</td><td>+</td><td>+</td><td></td><td></td><td>+</td></tr><tr><td>KA09</td><td>Kasuga</td><td>Honko</td><td>75</td><td></td><td>+</td><td></td><td></td><td></td><td>+</td></tr><tr><td>KA10</td><td>Kasuga</td><td>Honko</td><td>75</td><td></td><td>+</td><td>+++</td><td></td><td></td><td>+</td></tr><tr><td>KA15</td><td>Kasuga</td><td>Honko</td><td>75</td><td>+</td><td>+</td><td>+</td><td></td><td></td><td>+</td></tr><tr><td>KA19</td><td>Kasuga</td><td>Honko</td><td>75</td><td>+</td><td>+</td><td>+</td><td></td><td></td><td>+</td></tr><tr><td>KA27</td><td>Kasuga</td><td>Honko</td><td>130</td><td></td><td>+</td><td>+++</td><td></td><td></td><td>+</td></tr></tbody></table>

CDR: coupled dissolution and reprecipitation

+: the mineral phase or the feature is observed in the sample

+++: the mineral phase or the feature was very common in the sample
**Fig. 4** Representative back-scattered electron images of pyrite (Py) disseminated in fine-grained quartz (Qz). **a** Anhedral pyrite grains with minor enargite (En). **b** Euhedral pyrite grains in a quartz matrix. Also observed is anhedral rutile (Rt). **c** Small aggregate of anhedral pyrite grains. **d** Fine and coarse euhedral pyrite grains. **e** Aggregate of small anhedral grains. **f** Aggregate of euhedral pyrite grains. **b, d, f** Euhedral, relatively coarse-grained pyrites may indicate precipitation under relatively steady physico-chemical conditions. **c, e** Anhedral, fine-grained pyrites may represent rapid crystallization. See text for the discussion. **a–d, f** Akeshi deposit. **e** Kasuga deposit

<figure><img src="image_8.png" /></figure>

is an exception and shows sharp compositional contrasts (Figs. 6 and 7). For example, zones enriched in Cu were observed in 14 out of 21 EPMA pyrite maps. These include thin (< 10 µm) Cu- and Co-rich rims in pyrite grains (e.g., Fig. 6b, d), as well as Cu-rich zones with porous, CDR-like textures (Fig. 7a, b).

## LA-ICP-MS trace element analysis

We conducted LA-ICP-MS spot analyses on primary, inclusion- and pore-free pyrite crystals, as well as in pyrite aggregates and pyrite grains showing porous CDR textures. Thus, our LA-ICP-MS data represent an “efficient microsampling of small volumes” (Cook et al. 2009) of pyrite in the Akeshi and Kasuga deposits, and hence, we interpret the laser ablation signal profiles considering contributions from

elements incorporated in solid solution and as (micro-/nano-)inclusions within pyrite.

A statistical summary of LA-ICP-MS analyses is shown in Fig. 8 (results are listed in Online Resource 4). Copper, Ni, Pb, and Co were above detection in more than three quarters of the analyzed spots. Among these four elements, Cu showed the highest concentration (above detection median of 138.5 ppm; maximum: 7400 ppm). Arsenic, Au, Bi, Tl, Te, and Sb were detected in more than half of the spots. The maximum concentration of As (6500 ppm) was comparable to that of Cu, but the overall concentrations were lower than Cu (above detection median of As concentration was 61.5 ppm).

Scatter diagrams of selected elements are shown in Fig. 9. Gold shows positive correlations with Bi (*r* = 0.44) and Cu (*r* = 0.4; Fig. 9a). Bismuth, Te, and Sb, on the other hand, display positive correlations between each other (Fig. 9b, d). Cobalt also shows a positive correlation with Ni (*r* = 0.57).

<figure><img src="image_9.png" /></figure>
**Fig. 5** Representative back-scattered electron images of pyrites (Py) with textures indicating coupled dissolution-reprecipitation (CDR). **a** Pyrite grains from the Kasuga deposit with multiple layers of CDR textures. The white square is the area of elemental mapping in Fig. 7. **b** Pyrite from the Akeshi deposit showing a light gray rim. **c** A close-up image of the contact between pristine primary and altered pyrite. Akeshi deposit. For **a** and **b**, inferred reaction fronts are shown in Online Resource 1 Fig. S3

<figure><img src="image_10.png" /></figure>

Arsenic and Tl did not show clear correlations with any of these elements. Correlation coefficients of Cu, Ni, Pb, Co, As, Au, Bi, Tl, Te, and Sb are shown in Online Resource 1 Fig. S4.

## Discussion

### Polymetallic Au-bearing nanoparticles in pyrite

Our LA-ICP-MS depth profiles suggest that Au in the analyzed pyrites is present both in solid solution and as sub-micron polymetallic Au-bearing inclusions (Figs. 10 and 11). Although LA-ICP-MS depth profiles alone cannot confirm the presence of metal nanoparticles, the composition of the apparent inclusions (discussed below) seems to be similar to those reported from other high-sulfidation epithermal Au deposits (Deditius et al. 2011).

Gold and Ag signals appear at about half of the ablation time, coinciding with increasing Cu and Bi signals (gray band in Fig. 10a, at ~20–35 s; see Online Resource

1 Fig. S5 for BSE image of this analysis spot). This profile indicates that Au is distributed heterogeneously in pyrite, probably forming clusters of Au–Ag–Bi–Cu-bearing nanoparticles.

Profiles of Au and Te (Fig. 10b) display matching sharp peaks (gray bands). For example, at ~35 s, signal peaks of Au and Te correspond to peaks in the Cu and Sb profile. These observations suggest the presence of Au–Te and Au–Te–Sb–Cu nanoparticles at the analyzed spot. These nanoparticles might be related to CDR processes (Cook et al. 2009; Deditius et al. 2011; Tooth et al. 2011), since several CDR-like textures are observed near this spot (See Online Resource 1 Fig. S6 for BSE image of this analysis spot).

The presence of Pb and Sb–Au nanoparticles is inferred based on similar patterns (gray bands; Fig. 11b). These nanoparticles may also have been formed via CDR processes, although there is no conclusive evidence—for instance, preservation of the external morphology of the parent phase—to support our suggestion that porous zones in the corresponding spot (Fig. 11a) were actually formed via CDR processes.
<figure><img src="image_11.png" /><figcaption>**Fig. 6** a, c Back-scattered electron images of representative pyrite grains from the Akeshi deposit. b, d Corresponding wavelength-dispersive spectrometry X-ray elemental mapping of Cu in pyrite (Py) showing the Cu-rich rim</figcaption></figure>

## Au–As relation in pyrite

Gold and As in pyrite commonly show a positive correlation in many hydrothermal Au deposits (Reich et al. 2005; Deditius et al. 2008, 2014; Morishita et al. 2018), which has led previous studies to propose a key role of As to facilitate Au incorporation into the pyrite structure (e.g., Kusebauch et al. 2019; Pokrovski et al. 2019). The Au–As correlation in pyrite was also examined by Morishita et al. (2019) for the Akeshi and Kasuga deposits. Based on 21 secondary ion mass spectrometry (SIMS) analyses, they concluded that the Au–As relationship is not clear.

Our significantly larger dataset (617 LA-ICP-MS analyses) complements a previous study of the Akeshi and Kasuga deposits (Morishita et al. 2019), allowing a broader evaluation of the Au–As relation. Our LA-ICP-MS data for Au and As in pyrite (Fig. 12) are consistent with concentrations reported in previous SIMS analyses (gray area in Fig. 12, from Morishita et al. 2019). Gold and As do not show a clear correlation,

which may be due to the variable forms of incorporation of these two elements in the studied pyrite grains. Gold incorporation has been thought to be enhanced by As¹⁻, which substitutes for S in pyrite (e.g., Fleet and Mumin 1997). Recently, more oxidized forms of As in pyrite, such as As³⁺, which substitutes for Fe (Deditius et al. 2008), and As⁰, which exists as “liquid-like” inclusions (Deditius et al. 2009b), have also been reported in pyrite from other high-sulfidation deposits (e.g., Pueblo Viejo, Dominican Republic). As³⁺ is reported to correlate with Au in the Yanacocha high-sulfidation epithermal gold deposit in Peru (Deditius et al. 2008), but it remains unclear whether this correlation truly reflects the ability of As³⁺ in pyrite to enhance Au incorporation (i.e., higher pyrite-fluid partition coefficient for Au). Instead, it may simply reflect the composition of the mineralizing fluid and/or the depositional environment (e.g., pressure, temperature) as suggested by several other studies (e.g., Reich and Becker 2006; Deditius et al. 2009a; Reich et al. 2013; Tardani et al. 2017; Román et al. 2019; Herazo et al. 2021). Since we do not have

<figure><img src="image_12.png" /></figure>
<figure><img src="image_13.png" /><figcaption>**Fig. 7** a Detailed back-scattered electron image of a coupled dissolution-reprecipitation rim. The area corresponds to the white rectangle in Fig. 5a.
b Corresponding wavelength-dispersive spectrometry X-ray elemental mapping of Cu</figcaption></figure>

clear evidence for As speciation in pyrite in the studied de-
posits, we consider it difficult to evaluate the role of As in Au
uptake in the Akeshi and the Kasuga deposits. Future spectro-
scopic studies (e.g., synchrotron XANES and EXFAS) should
provide further evidence on As oxidation state(s) and incor-
poration into pyrite (e.g., Simon et al. 1999a; Savage et al.
2000; Manceau et al. 2020) and its influence on Au
enrichment.

Cu enrichment events recorded in pyrite

Cu as an indicator of changes in ore-forming environments

In the Akeshi and Kasuga pyrites, the Cu content shows a larger variation than other elements. In addition, Au and Bi contents show positive correlations with Cu content (Fig. 9). Here, we discuss the Cu concentration and distribution in pyrite and evaluate its use to track changes in elemental supply and/or physico-chemical conditions during the genesis of these deposits.

Copper negatively correlates with the Fe content in
pyrite (Fig. 13; *r* = −0.52). This negative correlation
between Fe and Cu is a common feature between the
Akeshi (orange triangles) and Kasuga (blue circles) de-
posits. Since inclusion-free pyrite grains were selected
for microprobe analyses, the effect of nanoparticles, espe-
cially those related to CDR and porous textures, is small.
Therefore, simple substitution of Fe²⁺ for Cu²⁺ in pyrite is
likely to be the major mechanism for Cu incorporation
into the Akeshi and Kasuga pyrite, as has been reported
in several hydrothermal deposits (e.g., Pacevski et al.
2008; Reich et al. 2013). However, the negative correla-
tion between Fe and Cu is not necessarily conclusive ev-
idence for the elemental substitution process, as it is com-
mon for Cu to occur as both solid solution and micro- to
nano-sized inclusions (Pacevski et al. 2008; Reich et al.
2013), and hence, their influence cannot be disregarded.

Previous studies have suggested that the trace element
content in the pyrite structure reflects the trace element
composition of the pyrite-forming fluid (Deditius et al.
2009a; Reich et al. 2013; Tardani et al. 2017). Based on
these findings, the Cu content in the Akeshi and Kasuga
pyrites may reflect changes in the Cu content of the ore-
forming fluid. High-temperature fluids may be responsible
for the Cu enrichment, considering that the solubility of
Cu in Cl-bearing hydrothermal fluids is enhanced at
higher fluid temperature (Liu and McPhail 2005;
Kouzmanov and Pokrovski 2012).

Rapid crystallization of Cu-rich pyrite

The trace element composition of pyrite can be influenced by
its growth rate, in addition to the composition of the pyrite-
forming fluid (e.g., Chouinard et al. 2005a). The fact that
relatively large euhedral grains (Fig. 4b, d, f) and aggregates
of small (< 25 µm; smaller than our LA-ICP-MS spot)
anhedral pyrite grains (Fig. 4c, e) both exist in the Akeshi
and Kasuga deposits suggests that the pyrite growth rate fluc-
tuated during ore formation. Among these morphologies, the
small grain aggregates are suggestive of rapid pyrite crystalli-
zation from the hydrothermal fluid (Huston et al. 1995; Simon
et al. 1999b; Román et al. 2019). All the aggregates analyzed
**Fig. 8** Graphical summary of LA-ICP-MS spot analyses (n = 617) of pyrite from the Akeshi and Kasuga deposits. Each boxplot illustrates the minimum, median, and maximum concentrations. The number of analyses above detection limit for each element is labeled above each box

<figure><img src="image_14.png" /></figure>

in this study (n = 55) show higher Cu concentrations (median, 1900 ppm) than the other pyrite types (Fig. 14). Moreover, the range of Cu contents in the aggregates is two orders of magnitude less than in all other pyrite types analyzed herein,

which is not observed for other elements (Fig. 14). Therefore, the processes that formed the pyrite aggregates (i.e., rapid crystallization) in the Akeshi and Kasuga deposits likely triggered efficient Cu precipitation.

**Fig. 9** Scatter diagrams of a Cu vs. Au, b Sb vs. Bi, c Cu vs. Bi, and d Sb vs. Te for pyrite from the Akeshi and Kasuga deposits. See text for the discussion

<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /></figure>
<figure><img src="image_17.png" /><figcaption>**Fig. 10** Representative depth-concentration profiles (time versus intensity) of selected isotopes obtained by LA-ICP-MS analysis of pyrite from the Kasuga deposit. The profiles correspond to the spots shown in the light blue circles in **a** Online Resource 1 Fig. S5 and **b** Online Resource 1 Fig. S6a. **a** Note that the signal intensity of Au and Ag increases at ~20 s coupled to an increase in the Cu and Bi signals, indicating a heterogeneous distribution of these elements. **b** The spiky signal of Au corresponds with Te (in the left gray zone) and Sb, Cu, and Te (in the right gray zone), indicating the presence of Au–Te and Au–Te–Sb–Cu nanoparticles</figcaption></figure>

In addition to Cu, the pyrite aggregates contain higher concentrations of Pb, Sb, Te, Bi, and Au than other pyrite types (Fig. 14). The 1st quartile elemental contents of pyrite aggregates are higher than the 3rd quartile elemental contents of other pyrite types. This may reflect abundant micro- to nano-sized sulfosalt inclusions and polymetallic nanoparticles in pyrite (e.g., Deditius et al. 2011; Ciobanu et al. 2012), probably due to supersaturation of the fluid with respect to these elements. In contrast, the concentration of Co, Ni, As, and Tl in all pyrites are similar (Fig. 14).

<figure><img src="image_18.png" /><figcaption>**Fig. 11** **a** Back-scattered electron image of porous pyrites (Py) from the Akeshi deposit presumably created via coupled dissolution–reprecipitation. The light blue circle denotes the location of the LA-ICP-MS analysis. **b** LA-ICP-MS depth–concentration profiles (time versus intensity) of selected isotopes in the area in **a**. The spiky peaks of Pb, Sb, and Au indicate the presence of micro- and/or nano-inclusions or aggregates of particles</figcaption></figure>

## Cu-rich overgrowth in pyrite

Sector and concentric zoning in pyrite can be related to two different types of control on trace element incorporation into pyrite (Chouinard et al. 2005a), either stereochemical influence of the crystal surface structure on the incorporation of trace elements, or concentric zoning that is believed to reflect temporal variation of the pyrite-forming fluid condition (e.g., temperature and/or fluid composition) (Chouinard et al. 2005a; Deditius et al. 2009a; Reich et al. 2013; Tardani et al. 2017).

The WDS X-ray elemental maps (Fig. 6) illustrate pyrite with Cu-rich overgrowths. Concentric micrometer-wide zones enriched in Cu were identified in some pyrite grains (Fig. 6b, d). Such sharp transition from the trace element-poor core
<figure><img src="image_19.png" /><figcaption>**Fig. 12** Scatter diagram of As versus Au content in pyrites from the Akeshi and Kasuga deposits. Black curve represents the solubility limit of Au as a function of As (Reich et al. 2005). Gray-shaded area is the range of Au and As concentration in pyrite from the Akeshi and Kasuga deposits reported by Morishita et al. (2019).</figcaption></figure>

to micron-scale, Cu-rich rim in pyrite from Au deposits in
epithermal environments has been suggested to form as a con-
sequence of rapid changes in the composition and/or physico-
chemical condition of the pyrite-forming fluid (Deditius et al.
2009a; Peterson and Mavrogenes 2014; Sanchez-Alfaro et al.

<figure><img src="image_20.png" /><figcaption>**Fig. 13** Scatter diagram of Fe versus Cu content (wt%) in pyrites of the Akeshi and Kasuga deposits. Data from electron probe microanalysis (EPMA)</figcaption></figure>

2016; Tardani et al. 2017; Wu et al. 2019b). Hence, it seems
reasonable to assume that the Cu content of the ore-forming
fluid in the Akeshi and Kasuga deposits experienced a rela-
tively sudden increase, which we refer to as a Cu enrichment
event. Since Au in pyrite positively correlates with Cu (Fig.
9a), the Cu enrichment event is likely related to at least a
portion of the Au mineralization in the Akeshi and Kasuga
deposits. In addition, the rapid crystallization of Cu-rich pyrite
suggested by the anhedral aggregates also seems to be closely
associated with the Cu enrichment event.

Replacement of primary pyrite by Cu-rich pyrite

The porous zones in some pyrite grains probably formed via CDR (Fig. 7a), as they usually match the Cu-rich zones identified in the EPMA elemental maps (Fig. 7b). This suggests that primary (early) pyrite was replaced by Cu-rich pyrite via CDR. However, a caveat is that our WDS X-ray maps yield relative concentrations, not absolute values. In addition, we were not able to measure the actual composition of the porous zones because most were not smooth enough for microprobe analysis; also, other mineral inclusions or nanoparticles may also be present (e.g., the extremely Cu-rich area in Fig. 6b).

In the Akeshi and Kasuga deposits, multiple CDR textures
can be recognized in a single pyrite grain of several hundred
microns. This observation suggests repetitive Cu enrichment
events in a relatively short period. Although no detailed data
regarding the timescale of pyrite CDR are available, previous
experimental studies showed that CDR could be a relatively
fast process, on the order of seconds to minutes (Renard et al.
2019) or days to weeks (Putnis 2009).

Implications for the mineralization processes

Possible cause of Cu enrichment events

According to Deditius et al. (2009a), abrupt and intermittent
changes in the composition of the pyrite-forming fluid, the Cu
enrichment event, are unlikely to occur by wall-rock reaction
or mineral deposition; rather, such changes are more likely to
occur when a (magmatic) fluid separates into two different
fluids or by mixing with another fluid. Within this framework,
one possible scenario is intermittent introduction of magmatic
fluids into the brittle, lower-pressured domain, where meteoric
water circulates at hydrostatic pressure (Fournier 1999). The
fluid injection itself can also promote hydraulic fracturing and
brecciation in the brittle domain and reduce fluid pressure and
temperature, which further enhances phase separation and/or
fluid mixing between magmatic and meteoric fluid (Sillitoe
1985; Fournier 1999). Hence, steep pressure and temperature
gradients are required to explain the Cu enrichment events.
Although direct evidence of magmatic fluid is lacking, most
previous studies agree that magmatic fluids supply a major

<figure><img src="image_21.png" /></figure>
<figure><img src="image_22.png" /><figcaption>**Fig. 14** Boxplots showing the compositional difference of pyrite aggregates (indicating rapid crystallization) and other pyrites in this study. For Cu, Pb, Sb, Te, Bi, and Au, concentrations are much higher in “Aggregates.” Horizontal lines in the boxes denote the median, and the bottom and top of each box denote the 1st and 3rd quartiles, respectively. The number of analyses above detection limit for each element is labeled below each box</figcaption></figure>

portion of the Au and Cu in high-sulfidation deposits (Vennemann et al. 1993; Arribas 1995; Heinrich 2005). This interpretation may also be consistent with findings in other hydrothermal Au deposits that pressure-driven, episodic hydrothermal fluid movements in brittle rheological zones are recorded in pyrite grains, sometimes bearing CDR textures (Peterson and Mavrogenes 2014; Rottier et al. 2016; Sanchez-Alfaro et al. 2016; Wu et al. 2019a).

## Role of pyrite in mineralization

Deposition of Au in epithermal deposits is commonly regarded as a consequence of H₂S removal, destabilization of Au hydrosulfide complexes (the dominant Au species in epithermal environments), and supersaturation of the fluid with respect to Au⁰ or other Au minerals (e.g., Drummond and Ohmoto 1985; Krupp and Seward 1987; Stoffregen 1987; Spycher and Reed 1989). Major processes that remove H₂S in epithermal environments include boiling, fluid oxidation, and sulfide precipitation (Pokrovski et al. 2014 and references therein).

Indeed, the Cu enrichment events in the Akeshi and Kasuga deposits are likely to have caused the rapid crystallization of pyrite, consuming some H₂S from the fluid. However, these events do not appear to have always resulted in the immediate supersaturation of Au in the fluid; this is evidenced by the fact that no pyrite grain in this study, including the rapidly crystallized aggregates, coexists with discrete or visible Au minerals. In addition, more than half of our LA-ICP-MS analyses detected Au in pyrite. Considering that pyrite can precipitate over a wide range of temperatures and physico-chemical conditions (e.g., Reich et al. 2013), it is plausible to consider a continuous formation of Au-bearing pyrite after the Cu enrichment events in the Akeshi and

Kasuga deposits. Our findings indicate, at least qualitatively, that the hydrothermal fluid was largely undersaturated (e.g., Simmons et al. 2016 and references therein) with respect to native Au for a considerable period during the mineralization process.

When the fluid is undersaturated with respect to Au, processes operating at the mineral surfaces can play an important role in Au mineralization (e.g., Kusebauch et al. 2019); the Au hydrosulfide complexes in the fluid can be adsorbed onto the growing pyrite surface, and consequently, Au can be incorporated in pyrite either as solid solution or nanoparticles (Widler and Seward 2002; Reich et al. 2005; Pokrovski et al. 2019). Adsorption on growing pyrite surfaces, a potentially effective mechanism for scavenging Au from a relatively Au-poor fluid, is also proposed to have played a key role in mineralization in other high-sulfidation deposits where pyrite is the main host of Au (Chouinard et al. 2005b; King et al. 2014). Therefore, we suggest that the same mechanism may also have contributed to Au mineralization in the Akeshi and Kasuga deposits. This is consistent with the Au–As plot (Fig. 12) where 80% of the data points plot below the empirical Au solubility line defined by Reich et al. (2005). This indicates that most of the Au in pyrite occurs in solid solution, probably as a consequence of adsorption and incorporation from solutions undersaturated with respect to native Au (or other Au-bearing mineral species).

## Conclusions

This study combined SEM, EPMA, and LA-ICP-MS analyses to unravel the microtextural characteristics and in situ chemical composition of pyrite, a main Au-bearing phase in the Akeshi and Kasuga high-sulfidation deposits in Kagoshima,
Japan. Overall, pyrite from the Akeshi and Kasuga deposits shows various morphologies, from euhedral crystals to anhedral aggregates of small pyrite grains; the latter indicating rapid crystallization from a hydrothermal fluid. In addition, multiple coupled dissolution-reprecipitation (CDR) textures are commonly recorded within single pyrite grains a few hundred micrometers in diameter. The CDR textures indicate the replacement of primary pyrite by Cu-rich pyrite.

Several trace elements including As, Au, Bi, Co, Cu, Ni, Pb, Sb, Te, and Tl were detected in more than half of the LA-ICP-MS analyses of all pyrite types. Pyrite aggregates formed by rapid crystallization have higher Cu, Au, Pb, Sb, Te, and Bi concentrations than large, euhedral pyrite grains, whereas Cu-rich overgrowths likely reflect rapid changes in the pyrite-forming fluid composition. In general, the Cu content in pyrite has a large variation compared to other elements and positively correlates with the Au content in pyrite. Most Au in pyrite is likely to exist as solid solution, whereas Au-bearing, polymetallic micron- to nanometer-sized inclusions may be present in some portions of pyrite grains. Gold and As in pyrite of the Akeshi and Kasuga deposits do not show a clear correlation, presumably reflecting multiple valence states of these elements.

Finally, we suggest that episodic input of metal-rich fluid, likely magmatic in origin, may be recorded in pyrite grains as sharp changes in trace element composition. Gold adsorption onto growing pyrite surfaces may also have played a key role in mineralization by scavenging Au from a relatively Au-poor fluid.

## Appendix: Analytical methods

Microtextural observations of the pyrite-bearing samples were carried out by using a JEOL JXA-8500F scanning electron microscope (SEM) at the Japan Agency for Marine-Earth Science and Technology (JAMSTEC) and a FEI Quanta 250 SEM at the Andean Geothermal Center of Excellence (CEGA) in the Department of Geology, University of Chile, Santiago.

Electron probe microanalysis (EPMA) of pyrite grains (362 spot analyses in total) was performed at the JAMSTEC using a field-emission JEOL JXA-8500F instrument. The beam energy was 20 keV, the beam current was 12 nA, and the beam size was 0.1 to 3 µm. Elements were detected using the following analyzing crystals: LiF for Mn Kα, Fe Kα, Cu Kα, Zn Kα; LiFH for Co Kα, Ni Kα; PETJ for S Kα, Pb Mα, Sb Lα; PETH for Bi Mα, Ag Lα; and TAP for As Lα, Se Lα. The standards used (Astimex MINM-25-53 and METM-25-44) included natural and synthetic sulfides, arsenides, tellurides, selenides, and native metals. Wavelength-dispersive spectrometry (WDS) X-ray maps were also acquired, with beam

energy 20 keV, beam current 50 nA, and beam size 0.1 to 3 µm.

LA-ICP-MS spot analyses of selected pyrite grains (617 points in total) were conducted using a 193 nm ArF excimer laser (Photon Machines Analyte 193) coupled to a quadrupole ICP-MS (Thermo Fisher Scientific iCAP Q) at the Mass Spectrometry Laboratory of the Andean Geothermal Center of Excellence (CEGA), Department of Geology, University of Chile. The analytical protocol is based on Román et al. (2019), with some modifications. Prior to each analytical session, the ICP-MS was tuned by ablating a NIST SRM 610 glass reference material (i.e., ²³⁸U⁺/²³²Th⁺ between 0.95 and 1.05), oxide production (ThO⁺/Th⁺ < 0.5%), and double-charged production (²²M⁺/⁴⁴Ca⁺⁺ < 0.01%). Ablation was carried out using a laser pulse frequency of 4 Hz, an energy density of ~1.5 J/cm², and a spot size of 25, 30, or 40 µm. Pure He was used as a carrier gas. Each spot was ablated for 30 s following 30 s of gas background collection. The following isotopes were monitored: ³⁴S, ⁵¹V, ⁵²Cr, ⁵³Cr, ⁵⁵Mn, ⁵⁷Fe, ⁵⁹Co, ⁶⁰Ni, ⁶³Cu, ⁶⁵Cu, ⁶⁶Zn, ⁶⁹Ga, ⁷²Ge, ⁷³Ge, ⁷⁵As, ⁷⁷Se, ⁸²Se, ⁹⁵Mo, ⁹⁷Mo, ¹⁰⁷Ag, ¹⁰⁹Ag, ¹¹¹Cd, ¹¹⁵In, ¹¹⁸Sn, ¹²⁰Sn, ¹²¹Sb, ¹²³Sb, ¹²⁵Te, ¹²⁸W, ¹⁹⁷Au, ²⁰²Hg, ²⁰⁵Tl, ²⁰⁶Pb, ²⁰⁷Pb, ²⁰⁸Pb, and ²⁰⁹Bi, using a total quadrupole sweep time of 0.55 s.

Both external and internal standard calibrations were applied (Longerich et al. 1996). The MASS-1 pressed synthetic sulfide reference material (Wilson et al. 2002) was used as a primary standard, and total Fe concentrations obtained by EPMA were used as the internal standard. Additionally, the GSE-1G glass reference material was employed as a secondary standard using the concentration values reported by Jochum et al. (2005). External standard measurements were performed at the beginning and end of each analysis round of 20 spot analyses. Data integration and reduction was carried out using LAtools, an open-source Python package (Branson et al. 2019), with some modifications. Detection limits were calculated based on the equations of Longerich et al. (1996). For some elements with very low background, values from Heinrich et al. (2003) were used as the detection limit. The MASS-1 reference material may contain some heterogeneities, especially for Au (Wilson et al. 2002). Therefore, LA-ICP-MS depth profiles of the MASS-1 reference material were thoroughly inspected before calibration to avoid introducing Au heterogeneities, and MASS-1 readings were made in duplicate each analysis round. Except for the MASS-1 reference material, we did not exclude analyses affected by possible inclusions.

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1007/s00126-021-01053-4.

**Acknowledgements** We thank H. Kon and R. Maekawa at Mitsui Kushikino Mining Co., Ltd., and K. Shinagawa, M. Inoue, and K. Yamamiya at Kasuga Mines Co., Ltd., for the sampling support in the

<figure><img src="image_23.png" /></figure>
Akeshi and Kasuga deposits. We gratefully acknowledge J. Torimoto, K. Uematsu, M. Yano, K. Fujinaga, and T. Nozaki for the help with sample preparation and EPMA analyses. We also thank faculty and student members of the Millennium Nucleus for Metal Tracing, University of Chile, for their dedicated support in LA-ICP-MS analysis and technical discussion. We appreciate J.W. Hedenquist and A. Arribas for their thorough and constructive reviews that significantly improved this manuscript.

**Funding** This research was financially supported by Japan Society for the Promotion of Science (JSPS) KAKENHI Grants No. JP15H05771 and JP20H05658 to Y.K., No. JP17H01361 to N.K., No. JP18K14168 to K.Y., and No. JP20J22038 to M.I. LA-ICP-MS analysis was supported by ANID-FONDEQUIP instrumentation grant EQM120098. Additional funding was provided by ANID through Millennium Science Initiative Program (NCN13_065). M.I. is grateful to the University of Tokyo School of Engineering Study Abroad Musha Shugyo Program for their financial support for travel to and stay in University of Chile.

**Data and materials availability** Not applicable.

**Code availability** Not applicable.

**Declarations**

**Conflict of interest** The authors declare no competing interests.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

**References**

Arribas A (1995) Characteristics of high-sulfidation epithermal deposits, and their relation to magmatic fluid. Mineral Assoc Can Short Course 23:419–454

Branson O, Fehrenbacher J, Vetter L, Sadekov AY, Eggins SM, Spero HJ (2019) LAtools: a data analysis package for the reproducible reduction of LA-ICPMS data. Chem Geol 504:83–95. https://doi.org/10.1016/j.chemgeo.2018.10.029

Chouinard A, Paquette J, Williams-Jones AE (2005a) Crystallographic controls on trace-element incorporation in auriferous pyrite from the Pascua epithermal high-sulfidation deposit, Chile-Argentina. Can Mineral 43:951–963. https://doi.org/10.2113/gscanmin.43.3.951

Chouinard A, Williams-Jones AE, Leonardson RW, Hodgson CJ, Silva P, Téllez C, Vega J, Rojas F (2005b) Geology and genesis of the multistage high-sulfidation epithermal Pascua Au-Ag-Cu deposit, Chile and Argentina. Econ Geol 100:463–490. https://doi.org/10.2113/gsecongo.100.3.463

Ciobanu CL, Cook NJ, Utsunomiya S, Kogagwa M, Green L, Gilbert S, Wade B (2012) Gold-telluride nanoparticles revealed in arsenic-free pyrite. Am Mineral 97:1515–1518. https://doi.org/10.2138/am.2012.4207

Cook NJ, Ciobanu CL, Mao J (2009) Textural control on gold distribution in As-free pyrite from the Dongping, Huangtuliang and Hougou gold deposits, North China Craton (Hebei Province, China). Chem Geol 264:101–121. https://doi.org/10.1016/j.chemgeo.2009.02.020

Deditius AP, Utsunomiya S, Renock D, Ewing RC, Ramana CV, Becker U, Kesler SE (2008) A proposed new type of arsenian pyrite: composition, nanostructure and geological significance. Geochim Cosmochim Acta 72:2919–2933. https://doi.org/10.1016/j.GCA.2008.03.014

Deditius AP, Utsunomiya S, Ewing RC, Chryssoulis SL, Venter D, Kesler SE (2009a) Decoupled geochemical behavior of As and Cu in hydrothermal systems. Geology 37:707–710. https://doi.org/10.1130/G25781A.1

Deditius AP, Utsunomiya S, Ewing RC, Kesler SE (2009b) Nanoscale “liquid” inclusions of As-Fe-S in arsenian pyrite. Am Mineral 94: 391–394. https://doi.org/10.2138/am.2009.3116

Deditius AP, Utsunomiya S, Reich M, Kesler SE, Ewing RC, Hough R, Walshe J (2011) Trace metal nanoparticles in pyrite. Ore Geol Rev 42:32–46. https://doi.org/10.1016/j.oregeorev.2011.03.003

Deditius AP, Reich M, Kesler SE, Utsunomiya S, Chryssoulis SL, Walshe J, Ewing RC (2014) The coupled geochemistry of Au and As in pyrite from hydrothermal ore deposits. Geochim Cosmochim Acta 140:644–670. https://doi.org/10.1016/j.gca.2014.05.045

Drummond SE, Ohmoto H (1985) Chemical evolution and mineral deposition in boiling hydrothermal systems. Econ Geol 80:126–147. https://doi.org/10.2113/gsecongo.80.1.126

Faure K, Matsuhisa Y, Metsugi H, Mizota C, Hayashi S (2002) The Hishikari Au-Ag epithermal deposit, Japan: oxygen and hydrogen isotope evidence in determining the source of paleohydrothermal fluids. Econ Geol 97:481–498. https://doi.org/10.2113/gsecongo.97.3.481

Fleet ME, Mumin AH (1997) Gold-bearing arsenian pyrite and marcasite and arsenopyrite from Carlin Trend gold deposits and laboratory synthesis. Am Mineral 82:182–193. https://doi.org/10.2138/am-1997-1-220

Foumier RO (1999) Hydrothermal processes related to movement of fluid from plastic into brittle rock in the magmatic-epithermal environment. Econ Geol 94:1193–1211. https://doi.org/10.2113/gsecongo.94.8.1193

Garwin S, Hall R, Watanabe Y (2005) Tectonic setting, geology, and gold and copper mineralization in Cenozoic magmatic arcs of Southeast Asia and the West Pacific. In: Hedenquist JW, Thompson JFH, Goldfarb RJ, Richards JP (eds) Economic Geology One Hundredth Anniversary Volume. Society of Economic Geologists, pp 891–930. https://doi.org/10.5382/AV100.27

Gopon P, Douglas JO, Auger MA, Hansen L, Wade J, Cline JS, Robb LJ, Moody MP (2019) A nanoscale investigation of Carlin-type gold deposits: an atom-scale elemental and isotopic perspective. Econ Geol 114:1123–1133. https://doi.org/10.5382/econgoe.4676

Hedenquist JW, Aoki M (1991) Meteoric interaction with magmatic discharges in Japan and the significance for mineralization. Geology 19:1041–1044. https://doi.org/10.1130/0091-7613(1991)019<1041:MIWMDI>2.3.CO;2

Hedenquist JW, Aoki M, Shinohara H (1994a) Flux of volatiles and ore-forming metals from the magmatic-hydrothermal system of Satsuma Iwojima volcano. Geology 22:585–588. https://doi.org/10.1130/0091-7613(1994)022<0585:FOVAOF>2.3.CO;2

Hedenquist JW, Matsuhisa Y, Izawa E, White NC, Giggenbach WF, Aoki M (1994b) Geology, geochemistry, and origin of high sulfidation Cu-Au mineralization in the Nansatsu district, Japan. Econ Geol 89:1–30. https://doi.org/10.2113/gsecongoe.89.1.1

Heinrich CA (2005) The physical and chemical evolution of low-salinity magmatic fluids at the porphyry to epithermal transition: a thermodynamic study. Mineral Deposita 39:864–889. https://doi.org/10.1007/s00126-004-0461-9
Heinrich CA, Pettke T, Halter WE, Aigner-Torres M, Audétat A, Günther D, Hattendorf B, Bleiner D, Guillong M, Horn I (2003) Quantitative multi-element analysis of minerals, fluid and melt inclusions by laser-ablation inductively-coupled-plasma mass-spectrometry. *Geochim Cosmochim Acta* 67:3473–3497. https://doi.org/10.1016/S0016-7037(03)00084-X

Herazo A, Reich M, Barra F, Morata D, Id R (2021) Trace element geochemistry of pyrite from bitumen-bearing stratabound Cu-(Ag) deposits, Northern Chile. ACS Earth Space Chem 5:566–579. https://doi.org/10.1021/acsearthspacechem.0c00321

Huston DL, Sie SH, Suter GF, Cooke DR, Both RA (1995) Trace elements in sulfide minerals from eastern Australian volcanic-hosted massive sulfide deposits: part I. Proton microprobe analyses of pyrite, chalcopyrite, and sphalerite, and part II. Selenium levels in pyrite: comparison with δ³⁴S values and implications for the source of sulfur in volcanogenic hydrothermal systems. Econ Geol 90: 1167–1196. https://doi.org/10.2113/gsecongeol.90.5.1167

Imai H, Takahashi M, Yamaguchi M (1996) Relation between volcanic-hosted precious-and base-metal deposits and geothermal systems. *Resour Geol* 46:73–94. https://doi.org/10.11456/shigenchishitsu1992.46.73

Izawa E (1992) Evolution of volcanic hydrothermal systems in southern Kyushu. Geol Surv Jpn Rep 279:84–88. https://www.gsj.jp/data/rep-gsj/No279.pdf

Izawa E, Hayashi T (2018) Potential for porphyry copper mineralization below the Kasuga lithocap, Nansatsu district, Japan. Resour Geol 68:181–194. https://doi.org/10.1111/rge.12163

Izawa E, Urashima Y, Okubo Y (1984) Age of mineralization of the Nansatsu type gold deposits, Kagoshima, Japan-K-Ar dating of al-unite from Kasuga, Iwato and Akeshi. Min Geol 34:343–351. https://doi.org/10.11456/shigenchishitsu1951.34.343

Izawa E, Urashima Y, Ibaraki K, Suzuki R, Yokoyama T, Kawasaki K, Koga A, Taguchi S (1990) The Hishikari gold deposit: high-grade epithermal veins in Quaternary volcanics of southern Kyushu, Japan. J Geochem Explor 36:1–56. https://doi.org/10.1016/0375-6742(90)90050-K

Jochum KP, Willbold M, Raczek I, Stoll B, Herwig K (2005) Chemical characterisation of the USGS reference glasses GSA-1G, GSC-1G, GSD-1G, GSE-1G, BCR-2G, BHVO-2G and BIR-1G ysing EPMA, ID-TIMS, ID-ICP-MS and LA-ICP-MS. Geostand Geoanal Res 29:285–302. https://doi.org/10.1111/j.1751-908x.2005.tb00901.x

King J, Williams-Jones AE, van Hinsberg V, Williams-Jones G (2014) High-sulfidation epithermal pyrite-hosted Au (Ag-Cu) ore formation by condensed magmatic vapors on Sangihe Island, Indonesia. Econ Geol 109:1705–1733. https://doi.org/10.2113/econgeol.109.6.1705

Kouzmanov K, Pokrovski GS (2012) Hydrothermal Controls on Metal Distribution in Porphyry Cu (-Mo-Au) Systems In: Hedenquist JW, Harris M, Camus F (eds) Geology and Genesis of Major Copper Deposits and Districts of the World: A Tribute to Richard H Sillitoe. Society of Economic Geologists, pp 573-618. https://doi.org/10.5382/SP.16.22

Krupp RE, Seward TM (1987) The Rotokawa geothermal system, New Zealand: an active epithermal gold-depositing environment. Econ Geol 82:1109–1129. https://doi.org/10.2113/gsecongeol.82.5.1109

Kusebauch C, Gleeson SA, Oelze M (2019) Coupled partitioning of Au and As into pyrite controls formation of giant Au deposits. Sci Adv 5:eaav5891. https://doi.org/10.1126/sciadv.aav5891

Liu W, McPhail DC (2005) Thermodynamic properties of copper chloride complexes and copper transport in magmatic-hydrothermal solutions. Chem Geol 221:21–39. https://doi.org/10.1016/j.chemgeo.2005.04.009

Longerich HP, Jackson SE, Günther D (1996) Laser ablation inductively coupled plasma mass spectrometric transient signal data acquisition

and analyte concentration calculation. J Anal At Spectrom 11:899–904. https://doi.org/10.1039/JA9961100899

Mahony SH, Wallace LM, Miyoshi M, Villamor P, Sparks RJ, Hasenaka T (2011) Volcano-pectonic interactions during rapid plate-boundary evolution in the Kyushu region, SW Japan. Bull Geol Soc Am 123: 2201–2223. https://doi.org/10.1130/B30408.1

Manceau A, Merkulova M, Mathon O, Glatzel P, Murdzek M, Batanova V, Simionovici A, Steinmann SN, Paktunc D (2020) The mode of incorporation of As(–I) and Se(–I) in natural pyrite revisited. ACS Earth Space Chem 4:379–390. https://doi.org/10.1021/acsearthspacechem.9b00301

Metal Mining Agency of Japan (1993) Report on the detailed geological survey in the Nansatsu area in the fiscal year Heisei 4.

Morishita Y, Shimada N, Shimada K (2018) Invisible gold in arsenian pyrite from the high-grade Hishikari gold deposit, Japan: significance of variation and distribution of Au/As ratios in pyrite. Ore Geol Rev 95:79–93. https://doi.org/10.1016/J.OREGOREV.2018.02.029

Morishita Y, Hammond NQ, Momii K, Konagaya R, Sano Y, Takahata N, Ueno H (2019) Invisible gold in pyrite from epithermal, banded-iron-formation-hosted, and sedimentary gold deposits: evidence of hydrothermal influence. Minerals 9:447–447. https://doi.org/10.3390/min9070447

Nakamura K, Yamanaka K, Yamasaki T (1994) Gold mineralization of the Akeshi mine, Kagoshima prefecture. Resour Geol 44:155–171. https://doi.org/10.11456/shigenchishitsu1992.44.155

Pacevski A, Libowitzky E, Zivkovic P, Dimitrijevic R, Cvetkovic L (2008) Copper-bearing pyrite from the Čoka Marin polymetallic deposit, Serbia: mineral inclusions or true solid-solution? Can Mineral 46:249–261. https://doi.org/10.3749/canmin.46.1.249

Peterson EC, Mavrogenes JA (2014) Linking high-grade gold mineralization to earthquake-induced faultvalve processes in the Porgera gold deposit, Papua New Guinea. Geology 42:383–386. https://doi.org/10.1130/G35286.1

Pokrovski GS, Akinfiev NN, Borisova AY, Zotov AV, Kouzmanov K (2014) Gold speciation and transport in geological fluids: insights from experiments and physical-chemical modelling. Geol Soc Lond, Spec Publ 402:9–70. https://doi.org/10.1144/SP402.4

Pokrovski GS, Kokh MA, Proux O, Hazemann JL, Bazarkina EF, Testemate D, Escoda C, Boiron MC, Blanchard M, Aigouy T, Gouy S, de Parseval P, Thibaut M (2019) The nature and partitioning of invisible gold in the pyrite-fluid system. Ore Geol Rev 109:545–563. https://doi.org/10.1016/j.oregeorev.2019.04.024

Putnis A (2009) Mineral replacement reactions. Rev Mineral Geochem 70:87–124. https://doi.org/10.2138/rmg.2009.70.3

Reich M, Becker U (2006) First-principles calculations of the thermodynamic mixing properties of arsenic incorporation into pyrite and marcasite. Chem Geol 225:278–290. https://doi.org/10.1016/j.chemgeo.2005.08.021

Reich M, Kesler SE, Utsunomiya S, Palenik CS, Chryssoulis SL, Ewing RC (2005) Solubility of gold in arsenian pyrite. Geochim Cosmochim Acta 69:2781–2796. https://doi.org/10.1016/j.gca.2005.01.011

Reich M, Deditius A, Chryssoulis S, Li J-W, Ma C-Q, Parada MA, Barra F, Mittermayr F (2013) Pyrite as a record of hydrothermal fluid evolution in a porphyry copper system: a SIMS/EMPA trace element study. Geochim Cosmochim Acta 104:42–62. https://doi.org/10.1016/j.gca.2012.11.006

Renard F, Røyne A, Putnis CV (2019) Timescales of interface-coupled dissolution-precipitation reactions on carbonates. Geosci Front 10: 17–27. https://doi.org/10.1016/j.gsf.2018.02.013

Román N, Reich M, Leisen M, Morata D, Barra F, Deditius AP (2019) Geochemical and micro-textural fingerprints of boiling in pyrite. Geochim Cosmochim Acta 246:60–85. https://doi.org/10.1016/J.GCA.2018.11.034
Rottier B, Kouzmanov K, Wälle M, Bendezú R, Fontboté L (2016) Sulfide replacement processes revealed by textural and LA-ICP-MS trace element analyses: example from the early mineralization stages at Cerro de Pasco, Peru. Econ Geol 111:1347–1367. https://doi.org/10.2113/econgeo.111.6.1347

Sanchez-Alfaro P, Reich M, Driesner T, Cembrano J, Arancibia G, Pérez-Flores P, Heinrich CA, Rowland J, Tardani D, Lange D, Campos E (2016) The optimal windows for seismically-enhanced gold precipitation in the epithermal environment. Ore Geol Rev 79:463–473. https://doi.org/10.1016/j.oregorev.2016.06.005

Savage KS, Tingle TN, O'Day PA, Waychunas GA, Bird DK (2000) Arsenic speciation in pyrite and secondary weathering phases, Mother Lode Gold District, Tuolumne County, California. Appl Geochem 15:1219–1244. https://doi.org/10.1016/S0883-2927(99)00115-8

Shinohara H, Giggenbach WF, Kazahaya K, Hedenquist JW (1993) Geochemistry of volcanic gases and hot springs of Satsuma-Iwojima, Japan: following Matsuo. Geochem J 27:271–285. https://doi.org/10.2343/geochemj.27.271

Sillitoe RH (1985) Ore-related breccias in volcanoplutonic arcs. Econ Geol 80:1467–1514. https://doi.org/10.2113/gsecongeo.80.6.1467

Simmons SF, Brown KL, Tutolo BM (2016) Hydrothermal transport of Ag, Au, Cu, Pb, Te, Zn, and other metals and metalloids in New Zealand geothermal systems: spatial patterns, fluid-mineral equilibria, and implications for epithermal mineralization. Econ Geol 111: 589–618. https://doi.org/10.2113/econgeo.111.3.589

Simon G, Huang H, Penner-Hahn JE, Kesler SE, Kao L-S (1999a) Oxidation state of gold and arsenic in gold-bearing arsenian pyrite. Am Mineral 84:1071–1079. https://doi.org/10.2138/am-1999-7-809

Simon G, Kesler SE, Chryssoulis S (1999b) Geochemistry and textures of gold-bearing arsenian pyrite, Twin Creeks, Nevada: implications for deposition of gold in Carlin-type deposits. Econ Geol 94:405–422. https://doi.org/10.2113/gsecongeo.94.3.405

Spycher NF, Reed MH (1989) Evolution of a Broadlands-type epithermal ore fluid along alternative P-T paths: implications for the transport and deposition of base, precious, and volatile metals. Econ Geol 84: 328–359. https://doi.org/10.2113/gsecongeo.84.2.328

Stoffregen R (1987) Genesis of acid-sulfate alteration and Au-Cu-Ag mineralization at Summitville, Colorado. Econ Geol 82:1575–1591. https://doi.org/10.2113/gsecongeo.82.6.1575

Takeda T, Shimada N, Ueno H (2001) Paleomagnetic constraints on mineralization age of the Nansatsu-type gold deposits in Southern Kyushu, Japan. Resour Geol 51:239–248. https://doi.org/10.1111/j.1751-3928.2001.tb00095.x

Takenouchi S (1983) Fluid inclusion study of the Nansatsu-type gold deposits, Southern Kyushu. Min Geol 33:237–245. https://doi.org/10.11456/shigenchishitsu1951.33.237

Tardani D, Reich M, Deditius AP, Chryssoulis S, Sánchez-Alfaro P, Wrage J, Roberts MP (2017) Copper-arsenic decoupling in an active geothermal system: a link between pyrite and fluid composition. Geochim Cosmochim Acta 204:179–204. https://doi.org/10.1016/J.GCA.2017.01.044

Tokunaga M (1954) Geology and ore deposits of the Kasuga mine and Akeshi mine, in the Makurazaki district, Kagoshima Prefecture. Min Geol 4:205–212. https://doi.org/10.11456/shigenchishitsu1951.4.14_205

Tooth B, Ciobanu CL, Green L, O'Neill B, Brugger J (2011) Bi-melt formation and gold scavenging from hydrothermal fluids: an experimental study. Geochim Cosmochim Acta 75:5423–5443. https://doi.org/10.1016/j.gca.2011.07.020

Vennemann TW, Muntean JL, Kesler SE, O'Neil JR, Valley JW, Russell N (1993) Stable isotope evidence for magmatic fluids in the Pueblo Viejo epithermal acid sulfate Au-Ag deposit, Dominican Republic. Econ Geol 88:55–71. https://doi.org/10.2113/gsecongeo.88.1.55

Watanabe Y (2005) Late Cenozoic evolution of epithermal gold metallogenic provinces in Kyushu, Japan. Mineral Deposita 40: 307–323. https://doi.org/10.1007/s00126-005-0025-7

Watanabe Y, Sato R, Sulaksono A (2018) Role of potassic alteration for porphyry Cu mineralization: implication for the absence of porphyry Cu deposits in Japan. Resour Geol 68:195–207. https://doi.org/10.1111/rge.12165

Widler AM, Seward TM (2002) The adsorption of gold(I) hydrosulphide complexes by iron sulphide surfaces. Geochim Cosmochim Acta 66:383–402. https://doi.org/10.1016/S0016-7037(01)00791-8

Wilson SA, Ridley WI, Koenig AE (2002) Development of sulfide calibration standards for the laser ablation inductively-coupled plasma mass spectrometry technique. J Anal At Spectrom 17:406–409. https://doi.org/10.1039/b108787h

Wu Y-F, Evans K, Li J-W, Fougerouse D, Large RR, Guagliardo P (2019a) Metal remobilization and ore-fluid perturbation during episodic replacement of auriferous pyrite from an epizonal orogenic gold deposit. Geochim Cosmochim Acta 245:98–117. https://doi.org/10.1016/J.GCA.2018.10.031

Wu Y-F, Fougerouse D, Evans K, Reddy SM, Saxey DW, Guagliardo P, Li J-W (2019b) Gold, arsenic, and copper zoning in pyrite: a record of fluid chemistry and growth kinetics. Geology 47:641–644. https://doi.org/10.1130/G46114.1

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Terms and Conditions

Springer Nature journal content, brought to you courtesy of Springer Nature Customer Service Center GmbH (“Springer Nature”).

Springer Nature supports a reasonable amount of sharing of research papers by authors, subscribers and authorised users (“Users”), for small-scale personal, non-commercial use provided that all copyright, trade and service marks and other proprietary notices are maintained. By accessing, sharing, receiving or otherwise using the Springer Nature journal content you agree to these terms of use (“Terms”). For these purposes, Springer Nature considers academic use (by researchers and students) to be non-commercial.

These Terms are supplementary and will apply in addition to any applicable website terms and conditions, a relevant site licence or a personal subscription. These Terms will prevail over any conflict or ambiguity with regards to the relevant terms, a site licence or a personal subscription (to the extent of the conflict or ambiguity only). For Creative Commons-licensed articles, the terms of the Creative Commons license used will apply.

We collect and use personal data to provide access to the Springer Nature journal content. We may also use these personal data internally within ResearchGate and Springer Nature and as agreed share it, in an anonymised way, for purposes of tracking, analysis and reporting. We will not otherwise disclose your personal data outside the ResearchGate or the Springer Nature group of companies unless we have your permission as detailed in the Privacy Policy.

While Users may use the Springer Nature journal content for small scale, personal non-commercial use, it is important to note that Users may not:

1. use such content for the purpose of providing other users with access on a regular or large scale basis or as a means to circumvent access control;

2. use such content where to do so would be considered a criminal or statutory offence in any jurisdiction, or gives rise to civil liability, or is otherwise unlawful;

3. falsely or misleadingly imply or suggest endorsement, approval, sponsorship, or association unless explicitly agreed to by Springer Nature in writing;

4. use bots or other automated methods to access the content or redirect messages

5. override any security feature or exclusionary protocol; or

6. share the content in order to create substitute for Springer Nature products or services or a systematic database of Springer Nature journal content.

In line with the restriction against commercial use, Springer Nature does not permit the creation of a product or service that creates revenue, royalties, rent or income from our content or its inclusion as part of a paid for service or for other commercial gain. Springer Nature journal content cannot be used for inter-library loans and librarians may not upload Springer Nature journal content on a large scale into their, or any other, institutional repository.

These terms of use are reviewed regularly and may be amended at any time. Springer Nature is not obligated to publish any information or content on this website and may remove it or features or functionality at our sole discretion, at any time with or without notice. Springer Nature may revoke this licence to you at any time and remove access to any copies of the Springer Nature journal content which have been saved.

To the fullest extent permitted by law, Springer Nature makes no warranties, representations or guarantees to Users, either express or implied with respect to the Springer Nature journal content and all parties disclaim and waive any implied warranties or warranties imposed by law, including merchantability or fitness for any particular purpose.

Please note that these rights do not automatically extend to content, data or other material published by Springer Nature that may be licensed from third parties.

If you would like to use or distribute our Springer Nature journal content to a wider audience or on a regular basis or in any other manner not expressly permitted by these Terms, please contact Springer Nature at

onlineservice@springernature.com