Secular Variation of Early Cretaceous Granitoids in Kyushu, SW Japan: The Role of Mélange Rocks as a Possible Magma Source

The Early Cretaceous volcanic-arc granitic rocks from Kyushu, SW Japan are contemporaneous with the granitic rocks of the Yanshan Orogeny (SE China) along the eastern Eurasian continental margin. The secular geochemical variations of the whole-rock major elemental and Sr–Nd isotope data of the Early Cretaceous granitic rocks from Kyushu, SW Japan, as well as the zircon and apatite saturation temperatures, shows distinct changes during the Albian (∼115 to ∼100 Ma) as: (1) the mASI value of the rocks (i.e., Shiraishino granodiorites) decreases below 1, (2) the Sr–Nd isotopic data are relatively constant [87Sr/86Sri = 0.70471 to 0.70573; εNd(t) = + 0.2 to + 1.9] within different rock types including granites, granodiorites, tonalites, and adakitic rocks (i.e., the Shiraishino granodiorites), following the increase of 87Sr/86Sri and decrease of εNd(t) from Berriasian, and (3) higher maximum temperatures at ∼105 Ma. The secular changes indicate that important geodynamic changes occurred in the arc system of SW Japan as it changed from subduction-accretion during the Jurassic to continental arc during the Early Cretaceous. Thermodynamic partial melting modeling demonstrates that the Albian granitic rocks can be derived from mélange rocks, such as chlorite-actinolite schists, at moderate depth and variable redox conditions. It is concluded that the genesis of the Early Cretaceous granitic rocks from Kyushu, SW Japan, may be related to upwelling of the asthenosphere and hot corner flow into the mantle wedge caused by slab rollback, which followed a shallowing of the subduction angle and subsequent flat-slab subduction during the Late Jurassic. The resultant heat induced the partial melting of the mélange rocks that formed on and were transported from the subducted plate interface.


INTRODUCTION
Convergent margin magmatism occurred along the eastern Eurasian continental margin during the Late Jurassic to Late Cretaceous as the Paleo-Pacific plate subducted. Magmatism extended from the Korean Peninsula in the north to Vietnam in the south and represents a major period of juvenile crust formation and recycling in East Asia (Xu et al., 2002;Nguyen et al., 2004;Zhou et al., 2006;Kiminami and Imaoka, 2013;Shellnutt et al., 2013;Kim et al., 2016). Alternating periods of crustal shortening and extension accompanied by magmatism during the Late Mesozoic are collectively referred to as the Yanshan Orogeny (Wu, 2005;Zhou et al., 2006;Li and Li, 2007;Dong et al., 2018). The Yanshan Orogeny consists of an early period (∼170 to ∼135 Ma) of crustal shortening and a late (∼135 to ∼90 Ma) period of extension followed by a brief period (∼80 Ma) of weak crustal shortening (Dong et al., 2018). The majority of Early Cretaceous magmatic rocks from SE China, NE China, Korea, and SW Japan are compositionally similar to volcanic-arc granites and it appears that their emplacement migrated eastward toward the Eurasia margin over time (Chen and Jahn, 1998;Zhou and Li, 2000;Xu et al., 2002;Zhou et al., 2006;Kiminami and Imaoka, 2013;Kim et al., 2016).
The formation of Late Yanshanian granitic rocks is debated as there are different models proposed for different regions. On the one hand, plutons from SE China are thought to be related to extension-induced melting of mafic igneous protoliths at midlower crustal depths, due to the heat advection resulting from the underplating of basaltic magmas along the active continental margin (Zhou and Li, 2000;Zhou et al., 2006). On the other hand, the Early Cretaceous granitic rocks from the Korean Peninsula and SW Japan are explained by partial melting of lower continental crust, due to slab rollback (Kiminami and Imaoka, 2013;Kim et al., 2016), whereas, for NE China, by interaction of delaminated lower mafic continental crust with an upwelling asthenosphere (Xu et al., 2002;Hou et al., 2007;Lan et al., 2011). The difference in the melting regimes and tectonomagmatic models may be, in part, related to the paleoposition of Korea and Japan relative to China. Although there is along strike subduction, Early Cretaceous magmatism in NE Asia may be related to subduction of the Izanagi plate, whereas in eastern China, it is related to subduction of the Paleo-Pacific plate (Maruyama et al., 1997). Moreover, it is still debated whether the paleogeographic position of Korea and SW Japan belongs to the North China craton or the South China craton (Osanai et al., , 2006Ishiwatari and Tsujimori, 2003;Ernst et al., 2007;Oh and Kusky, 2007;Omori and Isozaki, 2011;Suga et al., 2017).
The upwelling of hot asthenosphere into the mantle wedge caused by slab rollback may have played a vital role in the formation of Early Cretaceous granitic rocks in NE Asia as this is thought to be the primary process that produces the adakitic rocks of SW Japan (Kamei, 2004;Imaoka et al., 2014;Kim et al., 2016). Alternatively, it is possible that the combination of asthenospheric upwelling and hot corner flow caused by slab rollback may have induced partial melting of mélange rocks, such as chlorite-actinolite schists, producing intermediate to silicic melts (Marschall and Schumacher, 2012;Hao et al., 2016). The Early Cretaceous granitic rocks, including adakitic rocks, in Kyushu, SW Japan, therefore, potentially provide an opportunity to investigate the genesis of granitic rocks related to the combination of asthenospheric upwelling and hot corner flow caused by slab rollback.
In this study, whole rock geochemical data from ∼150 granitic samples across Kyushu, SW Japan, are compiled in order to determine if there is a secular variation in the nature of the plutons and to assess the role subduction-related mélange rocks may play in their origin. Thermodynamic modeling using the program Rhyolite-MELTS is used to evaluate whether partial melting of mélange rocks is either partially or wholly responsible for the genesis of some Early Cretaceous felsic igneous rocks, including the adakitic plutons.

GEOLOGICAL BACKGROUND
The Japanese Islands comprise a number of subduction-related orogenic belts that were produced beneath the Eurasian continent over more than c. 500 Myr (Maruyama et al., 1997;Isozaki et al., 2010). The geotectonic framework in the inner zone of southwest Japan is made up of, from north to south, the Hida-Oki, Renge, Akiyoshi, Suo, Maizuru, Ultra-Tamba, Mino-Tanba, and Ryoke belts, with a tectonically downward-younging polarity that resulted from stepwise accretions during the Paleozoic to Mesozoic (Nishimura, 1998; Figure 1). In Kyushu, SW Japan, the Renge, Akiyoshi, Suo, and Ryoke belts occur on the northern side of Usuki-Yatsushiro tectonic line, where Cretaceous igneous rocks are distributed (Figure 1).
The Nagasaki high-P metamorphic rocks are distributed in western Kyushu (Figure 1). The rocks are a western extension of the high-P Sanbagawa metamorphic rocks of southwest Japan (Faure et al., 1988;Nishimura et al., 2004), which are the higher P/T members of a Cretaceous paired metamorphic belt that joins the high-T metamorphic rocks of the Ryoke plutono-metamorphic belt along the Median Tectonic Line. Toward the west, the Cretaceous plutonic rocks and related high-T metamorphic rocks, including the Higo high-T metamorphic rocks, are distributed throughout central to northern Kyushu.
By comparison, the diorites (n = 24), for the exception of 3 samples, are metaluminous, magnesian (Fe * = 0.36-0.50), and calc-alkalic to calcic in composition (Figures 2-4). They have higher Mg# ranging from 64 to 76 and lower SiO 2 contents of 52.2-59.4 wt.%, and have similar K 2 O/Na 2 O ratios of 0.3-0.9, compared to the granites, granodiorites, and tonalites. There is a weak positive trend of Na 2 O + K 2 O against SiO 2 . The trend line is extrapolated into the corresponding lines for the granites, granodiolites, and tonalites. LOI contents are 0.5-2.1 wt.%, which are similar to those of the granites, granodiorites, and tonalites.
The primitive mantle normalization patterns show enrichment of the less incompatible elements over the incompatible elements but they also have prominent depletions of Ba, Nb, and Ti ( Figure 5). The granodiorites have a subgroup of Hf depleted rocks but it is not certain if it is unique as Hf data was not reported for the diorites. The chondrite-normalized rare-earth element patterns of the granodiorites and tonalites show similar light rareearth element (LREE) enrichment and tend to have high La/Yb N ratios (5.6 to 35.4), whereas the diorites have moderately to slight LREE (La/Yb N < 5) enrichment (Figure 6). Rare earth elemental data are not available for the granites. The rocks generally do not display prominent depletion of Eu as the tonalites, diorites and most of the granodiorites (all but 4) have Eu/Eu * values ≥ 0.8. The relatively high Eu/Eu * values indicate that the rocks probably did not experience significant plagioclase fractionation.
Strontium and Nd isotopic data are only available for the granites, granodiorites, and tonalites (Supplementary Table S2).

Tectonomagmatic Classification
According to the classification scheme of Frost et al. (2001), the Early Cretaceous granitic rocks from Kyushu, are primarily magnesian, peraluminous to metaluminous, calc-alkalic to calcic in composition (Figures 2-4). Rocks of this composition are considered to be most commonly associated with subductionzone settings. There are some rocks in the data set that are ferroan but they are also highly silicic (SiO 2 > 72 wt.%). Magnesian and ferroan rocks tend to overlap at highly silicic compositions within cogenetic rock suites due to chemical differentiation as intermediate magnesian magmas eventually become ferroan after extensive fractional crystallization (Frost et al., 2001). Given that the ferroan rocks in this study are highly silicic and that they fall along the same negative or positive trends of TiO 2 , Al 2 O 3 , CaO, Mg#, and Na 2 O + K 2 O against SiO 2 as the magnesian granitic rocks, they are likely the most evolved compositions of an original magnesian parental magma (Figure 4). The trace element compositions of all samples fall unambiguously within the volcanic-arc/I-type granite field of the Y + Nb-Rb discrimination diagram (Pearce et al., 1984;Christiansen and Keith, 1996) (Figure 7). Moreover, the results indicate that the chemical differentiation trends identified in some of the plutons did not have a significant effect on their Y, Nb, and Rb concentrations and thus did not cause 'elemental drift' across the field boundaries (Forster et al., 1997). In other words, the major and trace elemental compositions are mutually supportive of each other, and consistent with the fact that the Japanese Islands comprise a number of subduction-related orogenic belts over more than c. 500 Myr (Maruyama et al., 1997;Isozaki et al., 2010). Thus, there is little doubt that the Early Cretaceous granitic rocks of Kyushu were exclusively emplaced at a volcanic-arc setting.

Secular Compositional and Temperature Variations
The whole-rock major elemental and Sr-Nd isotope data, as well as zircon (T Zrc ) and apatite (T Apt ) saturation temperatures, of the granitic rocks are plotted against time (Ma), in order to evaluate the secular compositional and temperature variations during the Cretaceous (Figures 8, 9 and Supplementary Tables S1, S2). The ages of the rocks are reported from high closure temperature minerals (U-Pb zircon, Rb-Sr whole-rock isochoron, K-Ar hornblende) and are interpreted to be crystallization ages. The zircon (T Zrc ) and apatite (T Apt ) saturation temperatures were calculated using the whole rock major and trace element data available for the granitic rocks (Harrison and Watson, 1984;Boehnke et al., 2013).
There is clear secular variation in the major element compositions (Figure 8). The most distinct major element change is the decrease in the alumina saturation index (mASI) value at ∼105 Ma (i.e., Shiraishino granodiorites). All other intrusions have near constant mASI values of ∼1.1 to ∼1.2, whereas the Shiraishino granodiorites have a wide range but drop below 1 and become metaluminous ( Figure 8A). In fact the Shiraishino granodiorite is the most distinct intrusion as it has high TiO 2 (∼0.66 wt.%) and low SiO 2 (∼63.8 wt.%) contents and low K 2 O/Na 2 O ratios (∼0.3) compared to the other granitic rocks.
However, in spite of the Shiraishino granodiorite, the overall trends show an increase in the Mg# and TiO 2 and Al 2 O 3 contents and a decrease in the K 2 O/Na 2 O ratio and SiO 2 over time (Figures 8B-F).
The Sr-Nd isotopes show two secular trends (Figures 9A,B). There are overall trends for both Sr and Nd, but they both show changes during the Albian (∼115 to ∼105 Ma). The 87 Sr/ 86 Sr i ratio increases from the Early Cretaceous to the Late Cretaceous, which is mirrored by an overall negative trend in the ε Nd (t) values. However, between ∼115 and ∼105 Ma it is clear that the 87 Sr/ 86 Sr i ratio flattens and becomes relatively constant (∼0.7047 to ∼0.7057) until it increases (∼0.7065) after 105 Ma and then drops (∼0.7051) at ∼90 Ma. Although there is less Nd isotopic data, a similar pattern as the Sr isotopes is observed at ∼115 to ∼105 Ma. The ε Nd (t) values of the tonalites, granodiorites and some granites are constant (∼ + 0.2 to ∼ + 1.9) before dropping to the lowest values (−4.7) at ∼95 Ma and increasing again (+ 0.4) at ∼85 Ma. Unlike the major elements where most of the compositional variation at ∼105 Ma is attributed to the Shiraishino granodiorites, the isotopic consistency between 115 and 105 Ma occurs within different rock types, implying a relatively uniform source for the Albian rocks. After ∼105 Ma, the increase of the 87 Sr/ 86 Sr i ratio and decrease of the ε Nd (t) values suggests a greater contribution of less radiogenic material was involved in the origin of the granites. It could be that these primitive magmas experienced crustal contamination during emplacement; or that the granites were derived by partial melting of the crust; or the source became more enriched due to an increase in subducted sediment into the mantle wedge. It would seem the first two options are more likely as the Sr and Nd isotopic values after 95 Ma return to the values observed during the Albian and some of the granites have higher Th/Nb PM values (>20) (Figure 10).
The calculated T Zrc and T Apt for the granites, granodiorites, and tonalites are similar: T Zrc = 653-740 • C and T Apt = 806-927 • C for the granites, T Zrc = 666-756 • C and T Apt = 803-1022 • C for the granodiorites, T Zrc = 674-739 • C and T Apt = 779-943 • C for the tonalites. The diorites did not yield meaningful T Zrc values, because their M values were outside (>1.9) of the experimental range (Miller et al., 2003), but T Apt (584-896 • C) estimates are the lowest of all rock types. The estimates of T Zrc for each rock type are lower than those of T Apt , implying the magma sources were undersaturated in Zr during partial melting and crystallization. In other words, the estimates of T Apt represent the maximum temperatures experienced by the granitic rocks, whereas, for T Zrc , are probably the minimum temperatures.
Over this time interval, the T Zrc and T Apt estimates of the granites, granodiorites, and tonalites increased from T Zrc = ∼660 • C and T Apt = ∼810 • C at 135 Ma and reached a peak of T Zrc = ∼750 • C and T Apt = ∼960 • C at ∼105 Ma. This is followed by lower T Apt temperatures to 840 • C at 81 Ma. Both the higher peaks of T Zrc and T Apt estimates are recorded in the granodiorites (i.e., Shiraishino granodiorite).
The secular variations observed within the whole-rock major elemental and Sr-Nd isotope data, and zircon (T Zrc ) and apatite (T Apt ) saturation temperatures reveal important geodynamic changes in the development of the Cretaceous arc system of  Christiansen and Keith (1996). Symbols as in Figure 2. VAG, volcanic arc granites; ORG, ocean ridge granites; syn-COLG, syn-collisional granites; WPG, within-plate granites; OR, ocean ridge. SW Japan. The rocks are principally magnesian and calc-alkalic, and have relatively high Eu/Eu * values (≥0.8). The magnesian nature of the rocks implies that their parental magmas were likely oxidized as this would lead to early crystallization of magnetite or Ti-rich magnetite and a relatively high Mg# (Osborn, 1959;Buddington and Lindsley, 1964). Moreover, the rocks are calcalkalic to calcic and have relatively high Eu/Eu * values which is likely a consequence of water-rich magma systems as plagioclase crystallization is suppressed under hydrous conditions (Arculus, 2003). In this regard, the relative uniformity of the rocks is supportive of a long term volcanic-arc setting as oxidizing and water-rich conditions are expected at subduction zones (Arculus, 2003;Plank et al., 2013).
It is clear that the Albian rocks, and perhaps the older (Valanginian) and youngest (Coniacian) rocks, are derived from a relatively uniform source as they have similar 87 Sr/ 86 Sr i ratios and ε Nd (t) values, as described above. The two exceptions are the oldest rocks (Berriasian) that are the most radiogenic and the ∼95 Ma rocks (Cenomanian) that are the most unradiogenic rocks. As mentioned previously, the Sr-Nd isotopic data of the ∼95 Ma rocks may reflect crustal contamination during emplacement, or partial melting of the crust. The fact that the Albian rocks are all less radiogenic than the Berriasian-Valanginian rocks in spite of their different rock types (granodiorite, tonalite, granite) indicates they are probably derived from a source that became slightly less radiogenic over time as a result of the material that was brought into the subduction system/mantle wedge. The T Zrc and T Apt estimates may also be supportive of the change in the subduction and melting geodynamics as the Albian rocks (i.e., Shiraishino granodiorite) have the higher mineral saturation temperatures, whereas older and younger rocks appear to have lower mineral saturation temperatures (Figure 9D). The higher temperatures for the granodiorites at Albian may imply the heat input into the magma source was higher during this time.
The magma source of the Shiraishino granodiorites was proposed to be subducted oceanic crust that melted due to asthenospheric upwelling and hot corner flow caused by slab rollback (Kamei, 2004;Kiminami and Imaoka, 2013;Imaoka et al., 2014;Kim et al., 2016). The proposed petrogenetic model is broadly consistent with the classification of the granodiorites as being related to slab failure using the Nb + Y-Nb/Y discrimination diagram (Figure 11), which implies that they were derived from partial melting of the upper layers of oceanic crust (Whalen and Hildebrand, 2019). In contrast, the other granitic rocks are thought to be derived by partial melting of mafic lower continental crust (Kamei, 2004), tonalitic middle crust (Kamei, 2002), and young and hot subducted oceanic crust . However, the other granitic rocks are also classified as slab failure-related (Figure 11), suggesting they were derived from the same source as the Shiraishino granodiorites. The model proposed for the genesis of the Shiraishino granodiorites, as well as the other granitic rocks, can be tested as there are constraints on the type of source material that is representative the upper part of subducted oceanic crust in the region. For example, the proximal Nishisonogi chlorite-actinolite schist is interpreted to be a mélange sequence that was derived from the upper part of subducted oceanic crust and melting of these rock types are known to produce silicic melts (Marschall and Schumacher, 2012;Hao et al., 2016). In the following section, we present thermodynamic modeling in order to determine if granitic melts similar to the Albian rocks can be derived by partial melting of the Nishisonogi chlorite-actinolite schist.
Equilibrium partial melting models were calculated using the Nishisonogi chlorite-actinolite schist as a proxy for mélange rocks in the region. Our goal is to determine if it is possible to derive some or all of the Albian granitic rocks from the Nishisonogo schist including: the Tsutsugatake granite, the Shiraishino granodiorites, the Manzaka tonalite, the Kunisaki and Kita-Taku high-Mg diorites, and the Kanmon andesites. We selected a starting composition that is typical of the Nishisonogi chlorite-actinolite schist (NK25; Mori et al., 2014; Table 1).
The thermodynamic modeling of granitic rocks can be modeled using the program Rhyolite-MELTS (Gualda et al., 2012). Rhyolite-MELTS as it is calibrated to the SiO 2 -TiO 2 -Al 2 O 3 -Fe 2 O 3 -Cr 2 O 3 -FeO-MnO-MgO-CaO-Na 2 O-K 2 O-P 2 O 5 -H 2 O bulk composition system. The software is optimized for silicic systems and enables the user to adjust modeling parameters, such FIGURE 10 | Sr/Y vs. primitive mantle-normalized Th/Nb (Th/Nb PM ) plot of the Cretaceous granitic (granites, granodiorites, tonalites, diorites) rocks from Kyushu. Primitive mantle values are from Sun and McDonough (1989). Symbols as in Figure 2.
as pressure (bars), relative oxidation state (f O 2 ), and water (wt.%) content of the system being modeled, in order to constrain the possible conditions and processes (i.e., fractional crystallization or partial melting). The crust of Kyushu is comprised of the Paleozoic-Mesozoic high-P metamorphic rocks and Jurassic accretionary rocks, and its thickness is estimated to be ∼26 km (Arai et al., 2009). However, it is likely that the crust was thinner during the Cretaceous as it would not include the additional material added since the Cenozoic. Consequently, the initial pressure used for the modeling is 0.5 GPa (i.e., 5 kbar, ∼20 km thickness) which is based on the likely depth of melting for the rocks. The relative oxidation state of the parental rock assemblage is not constrained and therefore we ran models using oxidation conditions that range from relatively reducing to relatively oxidizing (i.e., FMQ − 1, FMQ + 1, and FMQ + 3) with water contents of 0.5 wt%.
Although the models show that there is good agreement between the liquid evolution curves and the bulk composition of the granitic rocks, the K 2 O content in each model is consistently too low (Figure 12). The explanations for the discrepancy between the model curves and the K 2 O concentration are either related to the highly mobile nature of K + during metamorphism to hydrothermal alteration or to the fact that the Nishisonogi chlorite-actinolite schist is not a perfect analog of the mélange rocks that melted during the Albian (Wood et al., 1976;Kogiso et al., 1997;Zheng et al., 2011).

Tectonomagmatic Evolution of the Early Cretaceous Granitic Rocks
It is thought that slab rollback during the Early Cretaceous, following a shallowing of the subduction angle and subsequent flat-slab subduction during the Late Jurassic, changed the geotectonic setting of SW Japan from subduction-accretion to continental arc (Kiminami and Imaoka, 2013;Kim et al., 2016). As a result of the slab rollback, a slab gap could be created and filled with an upwelling asthenosphere, creating a 'hot corner' within the mantle wedge. The change in subduction geodynamics is expected to be reflected in the secular compositional and temperature variations of the granitic rocks.
The initiation of asthenospheric upwelling and hot corner flow caused by slab rollback likely began during the Berriasian as the tonalites have the most radiogenic Sr and Nd isotopic values [ 87 Sr/ 86 Sr i = 0.70430 to 0.70450; ε Nd (t) = + 0.5 to + 3.1; Figures 9A,B]. The magmatic source is likely juvenile or at least the rocks were derived from a radiogenic source under the lower continental crust or from subducted oceanic crust. Once slab rollback began, it is possible that the dip of the subducted slab increased and more sediments were transported along the subducted plate interface into the mantle wedge, likely causing a magma source to become more isotopically enriched. This is consistent with the increase of 87 Sr/ 86 Sr i ratios and decrease of ε Nd (t) values of granitic rocks during Beriasian to Albian (Figures 9A,B). The maximum temperatures experienced by the granitic rocks (T Apt ) also increased over time, indicating that there may be higher heat flow or melting conditions due to the input of upwelling asthenosphere into the mantle wedge and development of hot corner flow (Figures 9C,D). During the Albian and the Valanginian to Coniacian, the granitic rocks have similar 87 Sr/ 86 Sr i ratios and ε Nd (t) values independent of rock type (Figures 9A,B) as described above, although some granites at ∼95 Ma have higher 87 Sr/ 86 Sr i ratios and lower ε Nd (t) values, which might be explained by crustal contamination during emplacement or they were derived by partial melting of the crust as discussed above. The peak of maximum temperatures (T Apt ) occurred during the Albian, but the rate of temperature changes decreased in this period (Figures 9C,D). These lines of evidence suggest the inputs of heat and sediments into the mantle wedge may be close to 'steady state, ' probably due to the subducted slab reaching a stable angle. The stable subduction angle was likely the key development that enabled the generation of melts from a relatively uniform source for the granodiorites, tonalites, and possibly the least evolved granitic rocks.
The magmatic source during the Albian could be similar to mélange rocks, such as chlorite-actinolite schists, which are intensely mixed and metasomatised rocks on the subducting plate interface, and contain the chemical components characteristic for the generation of arc magma at subduction zone settings (Marschall and Schumacher, 2012). The mélange rocks may be transported as a diapir or plume into the hot corner of the mantle wedge where the dehydration and partial melting of the mélange rocks could generate a range of fluid and melt compositions (Marschall and Schumacher, 2012). Equilibrium partial melting modeling of chlorite-actinolite schists suggests that it is possible that they could produce magma similar to the Albian rocks at moderate depth and variable redox conditions, indicating that this may be a viable process that produced some of the Kyushu granitic rocks (Figure 12).

CONCLUSION
The Early Cretaceous granitic rocks from Kyushu, SW Japan, formed in a volcanic-arc setting. Important geodynamic changes in the heat supply and magmatic sources occurred during the Berriasian to Albian. Upwelling of the asthenosphere and hot corner flow into the mantle wedge caused by slab rollback, following a shallowing of the subduction angle and subsequent flat-slab subduction during the Late Jurassic, likely supplied the heat for partial melting of upper plate mélange rocks or the initiation of diapirs that formed on the subducted plate-mantle interface. This process was likely responsible for dehydration and partial melting of the subducted slab. The partial melting of the mélange rocks such as chlorite-actinolite schists at high to moderate pressure can yield melts that are similar to the bulk compositions observed in the Albian granitic rocks including the adakitic rocks of the Shiraishino granodiorites. Melting of the same source could also explain the relative Sr-Nd isotopic homogeneity of the Albian rocks.

DATA AVAILABILITY STATEMENT
All datasets generated for this study are included in the article/Supplementary Material.

AUTHOR CONTRIBUTIONS
KS contributed to conception and design of the study and wrote the manuscript. M-WY contributed to interpretation of data for the work. Both authors contributed to manuscript revision, read and approved the submitted version.