The Quaternary Kurobegawa Granite: an example of a deeply dissected resurgent pluton

1 Vol.:(0123456789) Scientific Reports | (2021) 11:22059 | https://doi.org/10.1038/s41598-021-01562-2 www.nature.com/scientificreports The Quaternary Kurobegawa Granite: an example of a deeply dissected resurgent pluton Hisatoshi Ito1*, Yoshiko Adachi1, Aitor Cambeses2, Fernando Bea2, Mayuko Fukuyama3, Koji Fukuma4, Ryuji Yamada5, Takashi Kubo6, Mami Takehara7 & Kenji Horie7,8 The Quaternary Kurobegawa Granite, central Japan, is not only the youngest known granitic pluton exposed on the Earth’s surface, it is one of few localities where both Quaternary volcanics and related plutons are well exposed. Here, we present new zircon U–Pb ages together with whole rock and mineral geochemical data, revealing that the Kurobegawa Granite is a resurgent pluton that was emplaced following the caldera-forming eruption of the Jiigatake Volcanics at 1.55 ± 0.09 Ma. Following the eruption, the remnant magma chamber progressively cooled forming the voluminous Kurobegawa pluton in the upper crust (~ 6 km depth) until ~ 0.7 Ma when resurgence caused rapid uplift and erosion in the region. This is the first study to document the detailed spatiotemporal evolution of resurgent pluton for a Quaternary caldera system. Our new findings may contribute significantly to understanding the fate of active caldera systems that can produce supereruptions. Plutons are an important component of Earth’s lithosphere and it is vital to better understand how they are assembled in the crust. Two different processes have been proposed to explain how plutons are emplaced in the upper crust: (1) plutons represent the remains of large-volume magma chambers that are genetically linked to caldera-forming eruptions1,2; or (2) plutons are emplaced incrementally through amalgamation of intrusions over millions of years and consequently not formed by the crystallization of large-volume magma chambers3–5. Seismic data may be consistent with protracted timescales of plutonic suite construction because geophysical studies have failed to locate large volumes of melt beneath active volcanic regions6–10, although this does not necessarily preclude the existence of voluminous melt11. Here, we focus on Earth’s youngest exposed pluton, the Kurobegawa Granite in the Hida Mountain Range (HMR), central Japan (Fig.1), to better understand pluton assembly. The Kurobegawa Granite intrudes the Jiigatake Volcanics, which together are inferred to be petrogenetically linked and constitute a volcano–plutonic complex13 (Fig.1b). Harayama etal.13 noted that the Kurobegawa Granite may be associated with a resurgent dome14 emplaced following the voluminous eruption of the Jiigatake Volcanics, although geochronological evidence to support their idea was lacking. Harayama etal.13 and Harayama15 also suggested that the rapid uplift and erosion of the Kurobegawa pluton was due to the prevailing E–W compressional stress regime, resulting in the pluton’s near-vertical eastward tilt. Ito etal.16 provided zircon U–Pb age data from 35 localities associated with the Kurobegawa Granite (sensu lato) and revealed that it is the Earth’s youngest exposed pluton amalgamated starting at 10Ma with multiple episodes of magmatic intrusion; the youngest at ca. 0.8Ma. In this study, we refer to the Kurobegawa Granite as a Quaternary pluton that intruded after the Jiigatake eruption (sensu stricto). Ito etal.17 analyzed modern fluvial (= detrital) zircons collected in the vicinity of the northwestern margin of the Kurobegawa pluton, revealing two magmatic pulses at ~ 2.3Ma and ~ 0.9Ma. Spencer etal.18 obtained high and low temperature thermochronogical data for the HMR and calculated the highest rates and magnitude of plutonic rock exhumation in Japan—one of the highest worldwide, which they associated with collision of the Izu–Bonin oceanic arc on the Philippine Sea Plate with the Japan arc (Fig.1a). Here, we provide new zircon U–Pb ages and geochemical data to better understand the spatiotemporal petrogenetic intrusion history of the Kurobegawa pluton and its predecessor, the Jiigatake Volcanics. We argue that OPEN 1Central Research Institute of Electric Power Industry, Chiba 270-1194, Japan. 2Department of Mineralogy and Petrology, University of Granada, 18002 Granada, Spain. 3Graduate School of Engineering Science, Akita University, Akita 010-8502, Japan. 4Department of Environmental System Science, Doshisha University, Kyotanabe 610-0394, Japan. 5National Research Institute for Earth Science and Disaster Resilience, Tsukuba, Ibaraki 305-0006, Japan. 6Asahi Town Board of Education, Toyama 939-0743, Japan. 7National Institute of Polar Research, Tokyo 190-8518, Japan. 8Department of Polar Sciences, The Graduate University for Advanced Studies, SOKENDAI, Tokyo 190-8518, Japan. *email: [email protected]en.or.jp 2 Vol:.(1234567890) Scientific Reports | (2021) 11:22059 | https://doi.org/10.1038/s41598-021-01562-2 www.nature.com/scientificreports/ the Kurobegawa Granite is a rare example of a deeply dissected resurgent pluton associated with a Quaternary caldera system and that the resurgence itself played a major role in uplift and erosion. Brief description of the Kurobegawa pluton and its associates The Kurobegawa Granite in the HMR is located to the west of the Itoigawa–Shizuoka Tectonic Line (ISTL), an active plate-bounding fault zone, and ~ 30km north of another Quaternary sizable pluton, the Takidani Granodiorite (Fig.1b). Both plutons are also bounded by N–S tectonic lineaments to the west19. Moreover, both plutons intrude volcanic rocks produced by caldera-forming eruptions and are interpreted to be the intrusive remnants of the Jiigatake–Kurobegawa and Hotaka–Takidani volcano–plutonic complexes13,19–21, although no topographic features of caldera margins remain due to rapid uplift and erosion. Two stages of mountain building have been suggested for the Quaternary tectonics in the HMR: a period of magmatic intrusion from 2.7 to 1.5Ma, followed by a period of E–W compressional stress from 1.4 to 0.5 Ma13,15. The Takasegawa fault (Fig.1b), a tectonic lineament > 30km in length and with a < 200m wide fracture zone, lies parallel to the ISTL in the middle of the HMR with mylonitized granitic rocks exposed along the lineament15,22. It is thought that E–W compressive stress since 1.4Ma focused along this lineament caused the eastern part of the HMR to be uplifted and tilted eastward, exposing deeper sections of the granite13,15; the credibility of this hypothesis will be discussed in this paper. Figure1. (a) Map of Japan with tectonic boundaries (white lines) and the location of the study site shown by a white rectangular. The base map is from Google Earth. (b) Geological map of the Hida Mountain Range showing the location of two Quaternary granitic plutons: the Kurobegawa Granite and the Takidani Granodiorite. The plutons lie between N–S tectonic lineaments (red dashed lines). ISTL Itoigawa–Shizuoka Tectonic Line, EP Eurasian Plate, NAP North American Plate, PSP Philippine Sea Plate, PP Pacific Plate. Yellow arrows in (a) indicate plate motion vectors (numbers indicating mm/year) relative to the North American Plate12. A black triangle in (b) is an active volcano (Tateyama Volcano). Geological map is modified from Seamless digital geological map of Japan 1: 200,000 (Geological Survey of Japan, AIST) using QGIS software (version 3.18.3) (https:// qgis. org/ ja/ site/) and Inkscape (version 1.0.2) (https:// inksc ape. org/). The base map in (b) is the Chiriin Tile (Japan and Its Surroundings) (http:// maps. gsi. go. jp) from the Geospatial Information Authority of Japan. 3 Vol.:(0123456789) Scientific Reports | (2021) 11:22059 | https://doi.org/10.1038/s41598-021-01562-2 www.nature.com/scientificreports/ The Kurobegawa Granite is a plutonic complex exposed over an area of ~ 100 km2 with vertical exposure from 700 to 2900m in elevation in the northern HMR. It is divided into three units (upper, middle, and lower) based on its geochemistry and texture23 (Fig.2) and locally intruded the Jiigatake Volcanics at its southeastern border. Along this contact, the granite becomes leucocratic with contact-metamorphism observable in the Jiigatake Volcanics15,23,24. The Jiigatake Volcanics are exposed as intra-caldera volcanic deposits of Quaternary age that extend 16km to the north–south and are locally > 3500m in thickness13. We estimate the volume of Jiigatake Volcanics to be ~ 300 km3 of dense rock equivalent (See “Methods” section for the volume estimate). The Jiigatake eruption occurred at least twice15,25, forming a composite caldera (the Jiigatake–Shirasawatengu Caldera Complex) coevally. The Kurobegawa Granite is interpreted to be the remnant of a felsic magma chamber that was zoned as a result of fractional crystallization23. It contains abundant microgranular mafic enclaves (MMEs) suggesting that mafic magma was repeatedly intruded into the base of a crystal mush magma body, eventually forming MMEs23. Sampling was performed to cover all the three units and also to supplement areas that are sparsely dated in the literature16–18 (Fig.2). Zircon U–Pb dating was performed using laser ablation–inductively coupled plasma–mass spectrometry (LA–ICP–MS). Geochemical analyses, including whole rock major and minor element analyses and trace element analyses in zircon, were performed to better understand magmatic history of the Jiigatake–Kurobegawa volcano–plutonic complex. Detailed investigation of a sample that includes both felsic and mafic components and an MME sample (samples name: KRW, Fig.3) was also performed using sensitive high-resolution ion microprobe (SHRIMP) for U–Pb dating and various geochemical techniques. Figure2. Geological map of the area surrounding the Kurobegawa Granite, showing the sample locations and zircon U–Pb ages. Dotted area denotes the Quaternary Kurobegawa Granite including ~ 2Ma portion. Large circles and small squares are ages in this study and the literature16–18, respectively. Dashed and dash-dot areas are youngest (< 1.0Ma) and oldest (> 2.0Ma) portions, respectively. Two blue lines indicate the granite becomes younger to the west. The three lithological units (upper, middle, and lower shown as Ku, Km, and Kl, respectively) of Wada etal.23 is delineated by yellow lines. Hida marginal belt comprises ultramafic rocks, sediments and metamorphic rocks. Geological map is modified from Seamless digital geological map of Japan 1: 200,000 (Geological Survey of Japan, AIST) using QGIS software (version 3.18.3) (https:// qgis. org/ ja/ site/) and Inkscape (version 1.0.2)(https:// inksc ape. org/). 4 Vol:.(1234567890) Scientific Reports | (2021) 11:22059 | https://doi.org/10.1038/s41598-021-01562-2 www.nature.com/scientificreports/ Results Zircon U–Pb analyses. LA-ICP-MS zircon U–Pb ages from 35 sites were newly obtained in this study (Fig.2; Supplementary TableS2). Representative cathodoluminescence (CL) images with U–Pb ages are shown in Supplementary Fig.S2. The CL images generally show oscillatory and patchy zoning. Four samples from the Jiigatake Volcanics were dated, yielding 1.62 ± 0.13Ma (KRB29; see Supplementary Fig.S1 for location), 1.51 ± 0.23Ma (KRB31), 1.25 ± 0.70Ma (KRB37), and 1.47 ± 0.18Ma (KRG40). These four ages overlap within uncertainty and yield a weighted mean age of 1.55 ± 0.09Ma (error: 95% confidence; MSWD = 0.93). Notably, the Jiigatake Volcanics contain some elongated zircons rich in melt inclusions (e.g., KRG40-9). The granitic rocks we analyzed ranged from ~ 250 to 0.7Ma. Assuming that the age of the Kurobegawa Granite is younger than the age of the Jiigatake Volcanics (1.55 ± 0.09Ma), we identify 20 of the 30 granite samples to be associated with the Kurobegawa Granite. Older samples include the following: A granite (KRB20) at the northeastern margin of the Kurobegawa Granite yielded 2.10 ± 0.16Ma. From the western margin of the Kurobegawa Granite, two granite samples yielded 3.50 ± 0.14Ma (KRG12) and 3.66 ± 0.18Ma (YH126), four granite samples (KRB08; KRB23; KRG14; KRG24) yielded 10–8Ma, and a metamorphosed granite (KRG16) yielded 198.0 ± 3.5Ma. Along the northeastern margin of the Kurobegawa Granite, a granite (KRB19) and a rhyolitic Figure3. (a) Photographs of a sample KRW that includes both felsic and mafic components and an MME (2: inset in a) that were analyzed in detail at Granada University. The scale in (a) is cm. The sampling site is shown as KRW3 in Supplementary Fig. S1. Photographs (b–d) are magnified images for subsamples KWR-1, KWR-2, and KWR-3, respectively. Subsamples KWR-4, -5, -5b, and -7 are also shown as stars in (a). Photos were taken by Aitor Cambeses. 5 Vol.:(0123456789) Scientific Reports | (2021) 11:22059 | https://doi.org/10.1038/s41598-021-01562-2 www.nature.com/scientificreports/ welled tuff (KRB25) yielded 63.6 ± 1.3Ma and 62.4 ± 2.20Ma, respectively. Finally, a granite from the northern margin of the Kurobegawa Granite (KRB11) yielded 253.1 ± 9.2Ma. For the Kurobegawa Granite, a sample collected in the northern part of the pluton (KRB18) containing both the felsic and mafic components were separately dated, which yielded 0.91 ± 0.20Ma and 1.06 ± 0.12Ma, respectively. Because the ages overlap within error and both dates were obtained from the same sample, they are shown as one point and a weighted mean age of 1.02Ma on the age distribution map (Fig.2). As a test of the accuracy of our LA–ICP–MS ages, samples KRB04 and KRB08 were collected from the sites of JP08 and JP16 of Spencer etal.18, respectively. New LA–ICP–MS zircon U–Pb ages for KRB04 and KRB08 are 1.05 ± 0.07Ma and 9.00 ± 0.12Ma, respectively. Previously reported zircon U–Pb ages for JP08 and JP16 are 1.14 ± 0.14Ma and 9.11 ± 0.29Ma, respectively, measured using a different LA–ICP–MS system at Curtin University, Australia18. The consistency of these ages supports the robustness of the U–Pb ages in this study. A granite portion (KRW-3) of the KRW sample (Fig.3) from the northern part of the pluton was dated by both LA–ICP–MS and SHRIMP methods, yielding 0.73 ± 0.03Ma (Supplementary TableS2) and 0.77 ± 0.04Ma (Supplementary TableS3), respectively. Representative CL images and SHRIMP U–Pb ages are shown in Supplementary Fig.S3. Thermobarometric results. Zircon saturation temperatures for the Kurobegawa Granite were calculated following Watson and Harrison26 (Supplementary Fig.S4a), yielding a range from 700 to 850°C. Zircon saturation temperatures for the Jiigatake Volcanics, a 2Ma granite, and a 10Ma granite were also plotted together with the Kurobegawa Granite (Supplementary Fig.S4b). The calculated temperatures of the Jiigatake Volcanics are generally warmer than the Kurobegawa Granite. Zircon and apatite saturation temperatures were also calculated using the KRW sample and Wada etal.23 data (Fig.4a). In general, the apatite saturation temperatures are 150–200°C higher than the zircon saturation temperatures, which is in agreement with the observation that apatites are an early crystallizing phase and occur as inclusions in zircon. MMEs in the KRW sample were further analyzed by amphibole-plagioclase thermobarometry to evaluate the depth of emplacement of the KRW. The magmatic temperature and pressure of the MMEs were estimated at 748 ± 22°C and 1.7 ± 0.2kbar, respectively (Fig.4b; see “Supplementary Information” for further explanation). The relatively low emplacement pressure (or shallow depth) of 1.7 kbar is in accordance with the fact that the Kurobegawa Granite is mostly porphyritic (Fig.5a–c) with miarolitic cavities extensively developed throughout the pluton23. The depth estimate is within the optimal emplacement window for shallow magma chambers, estimated using thermo-mechanical models33 (Fig.4b). Petrography and geochemical results. We studied in detail a Kurowegawa Granite sample (KRW) that includes both felsic and MME components (KRW-1, -2; -3, -4, 5 and 5b) and an MME sample (KRW-7) (Fig.3). We processed the sample to obtain representative portions of MME (KRW-1, -5 and -7), host granite (KRW-3 and -4) and intermingled zone (KRW-2 and -5b) in order to better quantify the textures and geochemical compositions. Whole rock major and trace element analytical results for these samples are shown in Supplementary TableS4. The MMEs are tonalites characterized by a medium-grained, equigranular hypidiomorphic texture. They contain plagioclase, amphibole, biotite and interstitial alkali feldspar and quartz. Accessory phases include apatite, sphene, epidote, zircon, magnetite and pyrite. The texture of the host granite is coarseto medium-grained and inequigranular. It is composed of large crystals of plagioclase, alkali feldspar, biotite and quartz in a fine-grained matrix with the same mineralogy and minor amphibole. Accessory phases in the granites are apatite, zircon, monazite, epidote, ilmenite and magnetite. In the intermingled zone, there are recognizable phases from both MMEs and host granites. However, in these areas, main phases show disequilibrium textures, including: normal and reversed oscillatory-zoned plagioclase, spongy-cellular plagioclase included in large plagioclase, patchy-zoned amphibole, and ocelar quartz megacrysts showing a corona of mafic minerals overgrowth with small grain quartz. Textures and mineral chemistry allow identification of different types of phases that formed during different crystallization stages. Chondrite normalized rare earth element (REE) patterns for whole rock and zircon from the Kurobegawa Granite, the Jiigatake Volcanics, and older granites are shown in Supplementary Fig.S5 and Supplementary Fig.S6, respectively. Supplementary Fig.S5a shows that the Jiigatake Volcanics (dotted lines) are enriched in REE relative to the Kurobegawa Granite, with the exception of KRB30. Sample KRB30, collected near the Jiigatake Volcanics, is the most REE enriched, indicating that the magma from which the KRB30 zircon crystallized at ~ 1.3Ma was most evolved. This sample also contains a high proportion of xenocrystic (> 50Ma) and antecrystic (> 2Ma) zircons (Supplementary Fig.S2; Supplementary TableS2), which may indicate that this magma assimilated older intrusive material in the shallow upper crust. See Supplementary Fig.S7 for a summary of zircon Hf isotope compositions for the Kurobegawa Granite compared to the measured U–Pb ages. The initial 176Hf/177Hf ratios, expressed as εHf(t), predominantly range from − 4 to − 1, which are comparable with the results from Spencer etal.18. Spencer etal.18 also measured older (~ 65Ma and ~ 10Ma) granites from surrounding areas. The overall agreement with the Kurobegawa Granite compositions and the older granites indicates that the Kurobegawa magma assimilated similar upper crustal material. The εHf(t) value for the mafic part (KRB18M) ranged from − 3.4 to − 0.9, which is more tightly clustered than the felsic part (KRB18F) ranging from − 2.9 to + 3.2, suggesting that the MME derives from a more homogenies (and deeper) source than the felsic part of the granite. Larger variation of εHf(t) in zircon from the felsic part shows the variable degree of assimilation during the formation of the Kurobegawa pluton. 6 Vol:.(1234567890) Scientific Reports | (2021) 11:22059 | https://doi.org/10.1038/s41598-021-01562-2 www.nature.com/scientificreports/ Discussion Although the newly obtained U–Pb ages generally agree with those in the literature16–18, we can now more accurately constrain the history of the Jiigatake–Kurobegawa volcano–plutonic complex. Based on our new results, the U–Pb age of the Jiigatake Volcanics is more tightly constrained at 1.55 ± 0.09Ma. Previously, the published ages for this formation were 1.55 ± 0.28Ma (K–Ar groundmass)34 and 1.90 ± 0.50Ma (U–Pb zircon)16. The accordance of our new U–Pb zircon with previously published K–Ar groundmass ages confirms that the Jiigatake magma accumulated and erupted at 1.55 ± 0.09Ma. Figure4. Temperature and pressure estimates for the Kurobegawa Granite. (a) Calculated zircon saturation temperature (TZrn) vs apatite saturation temperature (TAp) using the KRW sample (Fig.3) and data from Wada etal.23 (b) Pressure vs temperature for various amphibole (Amp)-plagioclase (Pl) pairs in MMEs. We grouped four Amp-Pl pairs based on total Al content (atoms per formula unit: apfu) in Amp and anorthite (An) content in Pl. Temperatures were calculated using the Amp-Pl thermometer of Molina etal.27 for groups 1 and 2 and Holland and Blundy28 for groups 3 and 4. Pressures were estimated using the Al in amphibole barometer29. The pressure range by the Al in sphene (Sph) barometry30 and the temperature range by apatite (Ap) saturation thermometry are also plotted. The red dotted line shows the experimentally determined water-saturated solidus for tonalite31 and the red solid line is the water-saturated solidus calculated using rhyolite-MELTS for diorite32. The transparent-grey area represents the proposed optimal emplacement window of shallow magma chambers (P = 2 ± 0.5 kbar), estimated using thermo-mechanical models33. 7 Vol.:(0123456789) Scientific Reports | (2021) 11:22059 | https://doi.org/10.1038/s41598-021-01562-2 www.nature.com/scientificreports/ The geochronology for the remnant magma chamber that, upon cooling, eventually became Kurobegawa Granite was not so simple. Based on our new and previous ages, magmatic activity along the north-eastern margin of the pluton (dash-dot lines in Fig.2) started between 2.2 and 2.1Ma; ~ 0.6 million years earlier than the Jiigatake volcanic activity. Following the eruption of the Jiigatake Volcanics at 1.55Ma, Kurobegawa magmatism migrated westward, defined by ages for the Kurobegawa Granite that become younger toward the west (blue arrows in Fig.2). Kurobegawa magmatism continued until 0.8–0.7Ma along the north-western margin (dashed lines in Fig.2). The youngest ages correspond to the areas of modern high temperature (> 60°C) hot springs unrelated to volcanic activity34, suggesting remnant heat from the cooling Kurobegawa Granite is likely driving hydrothermal activity. In contrast, the southern part of the Kurobegawa Granite yields ages mostly 1.4–1.0Ma, indicating that magma chamber was replenished soon after the voluminous 1.55Ma Jiigatake eruption and until ~ 1.0Ma. The weighted mean of all the ages from the Kurobegawa pluton (including ~ 2Ma portion) yields 0.99 ± 0.04Ma (95% confidence level; MSWD = 16; n = 187) (Fig.6). Therefore, zircon forming magmatic activity may have culminated 0.5–0.6 million years after the Jiigatake eruption in a broad sense. However, the relatively high MSWD of 16 for this age suggests it is not a single population. The probability density plot of zircon crystallization ages for the Kurobegawa pluton ages (Fig.7a) has multiple age peaks, compared to the well-defined single age peak for the Jiigatake Volcanics (Fig.7b). We have deconvoluted the peak spectra for the Kurobegawa pluton ages using the Unmix Ages program35 assuming three components. If we take the grouping of zircon crystallization ages to represent timing of magma intrusion and/or reheating of pre-existing mush, the result shows that the Kurobegawa pluton experienced three magmatic pulses at 1.25 ± 0.02Ma (fraction: 36%), 0.99 ± 0.03Ma (35%), and 0.70 ± 0.02Ma (29%) (Fig.7a). A similar result is also obtained using IsoplotR’s36 radial plot (Supplementary Fig.S8) confirming these populations of ages are robust. Although it is possible that the Figure5. Photomicrographs of the Kurobegawa Granite (a–c) and a ~ 10Ma granite (d) taken in crosspolarized light. (a) Typical porphyritic textures with micro crystals of quartz, plagioclase and biotite at the northwestern margin of the granite at Loc. KRB16. (b) Groundmass quartz shows a strong dynamic recrystallization and fracturing of a large plagioclase phenocryst near the Jiigatake Volcanics at Loc. KRB30. Many xenocrystic zircons were identified in this sample. (c) Groundmass quartz shows a weak dynamic recrystallization, and a quartz phenocryst shows dissolution at the south of the granite at Loc. KRB27. (d) Equigranular granite showing indications of neither tectonic shear nor fluid migration located to the west of the northern Kurobegawa Granite at Loc. KRB23. Qz quartz, Pl plagioclase, Bt biotite, Kfs potassium feldspar. 8 Vol:.(1234567890) Scientific Reports | (2021) 11:22059 | https://doi.org/10.1038/s41598-021-01562-2 www.nature.com/scientificreports/ Figure6. 238U–206Pb age distributions for the Kurobegawa Granite (including ~ 2Ma portion) zircons younger than 10Ma (7 grains > 10Ma are omitted). Individual grain ages with 2σ uncertainty are arranged in rank order. Analyses in grey represent statistical outliers and were excluded. Black and red horizontal bars represent mean U–Pb age of the Kurobegawa Granite and the ~ 1.55Ma Jiigatake eruption, respectively. MSWD mean square weighted deviation. Figure7. 238U–206Pb age distributions (histogram, probability density plot) for zircon younger than 10Ma for (a) the Kurobegawa Granite (including ~ 2Ma portion) and (b) the Jiigatake Volcanics. The inset in (a) shows the result of age peak deconvolution using Isoplot35. The deconvolution was performed for ages < 1.55Ma to analyze grouping of zircon ages after the ~ 1.55Ma Jiigatake eruption. In (b), individual grain ages with 2σ uncertainty are also plotted. n = number of grains plotted. 9 Vol.:(0123456789) Scientific Reports | (2021) 11:22059 | https://doi.org/10.1038/s41598-021-01562-2 www.nature.com/scientificreports/ Kurobegawa pluton experienced continuous, prolonged intrusion of magma from ca. 1.4 to 0.7Ma, the large number of zircon analyses presented here allow us to resolve three statistically distinct pulses that we interpret as magma intrusion events. Harayama etal.13 assumed nearly-vertical eastward tilting of the entire Jiigatake–Kurobegawa volcano–plutonic complex due to E–W compressional stress regime. One main reason is that mylonitic textures are observed in the western margin of the Kurobegawa Granite. However, this is unlikely because at the northwestern margin of the Kurobegawa Granite, no mylonitic texture has been identified in the Kurobegawa Granite (KRB16), nor in the 10Ma granite to the west (KRB23) (Fig.5a,d; Supplementary Fig.S1). In this area, mylonite exists to the west of the 10Ma granite. Therefore, we assume that the N–S mylonite zone formed before 10Ma and that both the 10Ma granite and the Kurobegawa Granite intruded along the mylonite zone. Harayama etal.13 identified mylonitic texture in the western margin of the Kurobegawa Granite at different localities from ours, but based on their descriptions we assume they may have misidentified the Kurobegawa Granite. Resurgence of the Kurobegawa granite. As noted above, it is unlikely that the Kurobegawa Granite tilted almost vertically to the east by the regional E–W compression13,15. Similarly, the suggestion that upper units lie progressively to the east23 (Fig.2) is doubtful because this assertion is based on the hypothesis of eastward tilting. It is also unlikely that the entire Kurobegawa Granite experienced a melting phase as a whole as envisaged by Wada etal.23 and a large chamber model may not be applicable here based on our new geochronology interpretation. Instead, the Kurobegawa Granite magma likely accumulated and progressively filled the void space caused by the 1.55Ma Jiigatake eruption from the south and east to the north and west, in a broad sense (Fig.2). Intrusion of this magma resulted in resurgence, which caused rapid uplift and erosion. Hence, we propose the Kurobegawa Granite is the result of resurgent magmatism because it is spatially restricted to the presumed caldera area (Fig.1b), the location of which suggested by previous authors13,15,19. Additionally, magmatism occurred at shallow crustal levels immediately after the voluminous (~ 300 km3) Jiigatake eruption. Assuming that the outcropping Kurobegawa Granite magma solidified at ~ 6km depth (equivalent to 1.7kbar assuming a mean crustal density of 2.75g cm−3) at ~ 1Ma, the uplift/erosion rate is ~ 6mm/year. Here, we interpret magmatic resurgence to be the dominant driving force for the observed rapid uplift37. This is similar, for example, to Toba Caldera, Indonesia, where resurgence caused parts of the Samosir Island to be uplifted 700m over the past 33,000year, equivalent to 21mm/year38. Recently, the intrusion depth and uplift/erosion rate of a ~ 5Ma granite to the west of the Kurobegawa Granite were estimated ~ 6–11km and 0.93–2.5mm/year, respectively25. Therefore, the uplift/erosion rate of the Kurobegawa pluton was higher than the surrounding areas, which also supports the resurgent uplift of the pluton. Prior to this study, identification of magmas driving resurgence in Quaternary systems could only be imaged by geophysical tools such as seismicity, gravity, and resistivity39,40. The Kurobegawa Granite is a rare case globally where the deep interior of a Quaternary caldera system can be investigated in the field. On the contrary, the rapid uplift is not restricted to the Kurobegawa Granite exposure, but it extends to the whole HMR18. Therefore, we suspect a deeper and more voluminous magmatic reservoir exists beneath the HMR (i.e., larger than the areal extent of the Kurobegawa Granite exposure) as in the case for other large caldera systems such as Toba41. Matsubara etal.6 found two low seismic velocity zones at depths of 2–4km and 12–20km beneath Tateyama Volcano, an active volcano to the west of the Kurobegawa Granite (Fig.1b) and assumed they correspond to partially molten rocks or magmatic reservoirs. We assume deep magmatic reservoirs such as the 12–20km magmatic reservoir imaged by Matsubara etal.6 likely contribute to the regional uplift of the HMR as a whole. This deep reservoir likely feeds hotter, mafic magma into the base of the shallower Kurobegawa reservoir, similar to models envisaged by other researchers42,43. Tectonic and magmatic evolution in the Hida Mountain Range (HMR). We present a model for the tectonic and magmatic evolution in the HMR as schematic snapshots at ~ 10Ma, ~ 2Ma, and ~ 0.7Ma (Fig.8). In the HMR, especially around the Kurobegawa Granite area, granite-producing magmatism prevailed at ~ 65 Ma16 (Fig.8a). Ito etal.16 assumed that the Kurobegawa Granite intruded ~ 65Ma granite because of the existence of ~ 65Ma xenocrystic zircons in the Kurobegawa pluton. This is further corroborated by the existence of ~ 60Ma xenocrysts at the southeastern part of the pluton (KRB30; Supplementary Fig.S2). N–S trending fracture zones, such as ISTL and the Takasegawa Fault (Fig.1), developed as a result of (1) back-arc spreading (Japan Sea opening) which culminated at ~ 15 Ma44 and/or (2) collision of the Izu-Bonin oceanic arc with the Japan arc since ~ 13Ma at the latest45. At ~ 10Ma (Fig.8a), granites intruded along N–S trending faults especially along the Takasegawa Fault and its branch faults (Fig.2). Granite intrusion occurred episodically to ~ 2 Ma16. This magmatism associated with the Takasegawa Fault occurred mainly to the west of the presently exposed Kurobegawa Granite area, whereas less voluminous magmatism also occurred to the east of the Kurobegawa Granite (Fig.8b) based on the identification of ~ 5Ma and ~ 2Ma granites near the northeastern margin of the Kurobegawa Granite (Fig.2). We have also analyzed 9–4Ma xenocrystic zircons in samples from the eastern margin of the pluton and adjacent Jiigatake Volcanics (e.g. samples KRB29, KRB30 in Supplementary TableS2). Therefore, we assume the existence of another fracture zone (Fracture zone2 in Fig.8a) to the east of the Kurobegawa Granite. At 1.55 ± 0.09Ma, a caldera-forming Jiigatake eruption occurred where the Jiigatake–Kurobegawa complex is presently exposed. Subsequently, resurgence occurred until ~ 0.7Ma accompanied by rapid uplift and erosion together with more regional magmatism beneath the resurgent pluton. This history has resulted in the presentday exposure of the Kurobegawa Granite (Fig.8c). Regionally, the HMR experienced three voluminous eruptions in the Quaternary: two eruptions at ca. 1.76Ma and 1.75Ma related to the Hotaka Andesite46 (Fig.1b) in addition to the ~ 1.55Ma eruption producing