Reconstruction of the deformation history of an active fault: implications from the Atera Fault, Central Japan

Atera Fault clay gouges were collected for age dating near Kawaue, Nakatsugawa City, Central Japan, and the results integrated within its complex geological history. The results form an internally consistent data set constrained by extensive geochronological data (AFTA, ZFTA, CHIME) and support the application of gouge dating in constraining timing of brittle deformation in Central Japan. The Atera illite age data complete recently obtained limited illite fault gouge age data from underground exposure in the Toki Granite; the new illite age data are identical within error. The age of the heterogenous welded tuff breccia zone (Atera 1) ranges from 40.6 ± 1.0 Ma to 60.0 ± 1.4 Ma, whereas ages of the fault core gouge sample (Atera 2) range from 41.8 ± 1.0 Ma to 52.7 ± 1.2 Ma. The finest < 0.1 µm fraction of the fault breccia and fault core gouge yield ages around 41 Ma, identical within error. The new illite age data indicate brittle faulting and a following geothermal event occurring in the Paleogene–Eocene, similar to the nearby Toki Granite area and confirm they were both synchronous with a post-intrusive pluton exhumation. The Atera Fault illite age data provide additional insights into an integrated, regional-scale record of the tectonic displacement of Central Japan and might be influenced by large-scale tectonic processes such as the Emperor sea mount kink around 55 to 46 Ma with fault initiation around 50 Ma and brittle fault cessation or reactivation around 40 Ma in the Eocene.


Introduction
Clay fault gouges are present in active brittle faults formed under near surface conditions (e.g., Scholz 2002).Active fault zones often document numerous geological processes ranging from the disintegration of host rocks with subsequent grain-size reduction and authigenesis of clay minerals (chlorite, smectite, illite) and muscovite.
In addition, secondary heating and melting of host rocks by frictional deformation (pseudotachylytes) and fault-related fluid flow might result in mineral vein formation (e.g., Holdsworth et al. 1997;Tagami 2012).Authigenic illite in fault gouges can reliably constrain the timing of fault activity by K-Ar, 40 Ar/ 39 Ar and Rb-Sr dating (see summaries in Vrolijk et al. 2018;Tsukamoto et al. 2020).Despite more than 50 years of research, the timing of slip on brittle faults in the Earth's upper crust are difficult to constrain and radiometric dating of fault zones remains challenging.Over the last decade, K-Ar data of authigenic illite formed in brittle fault zones are integrated with apatite fission-track (AFTA) and zircon fission-track (ZFTA) data as the frictional heating caused by fault motion can partially or totally reset the "clock" of the dating methods (Zwingmann et al. 2004(Zwingmann et al. , 2010;;Zwingmann and Mancktelow 2004;Tagami 2005Tagami , 2012;;Yamasaki et al. 2013).Timing of fault formation and its corresponding geofluid sources can be investigated through the robust isotopic memory of neo-crystallized authigenic clay minerals as its crystal structure includes H (as OH) and K, enabling constraint of fluid sources through H isotope analysis and timing through Argon dating (e.g., Hetzel et al. 2013;Mancktelow et al. 2015Mancktelow et al. , 2016;;Lynch and van der Pluijm 2016;Heinecke et al. 2019).Limited studies apply a whole range of dating methods such as K-Ar, 40 Ar/ 39 Ar, and Rb-Sr dating of syntectonic fault-related illite (e.g., Middleton et. 2014;Uysal et al. 2020).U-Pb dating of carbonite geochronology can provide direct timing constraints for geological processes involving diagenesis, fluid flow, brittle faulting via carbonate slickenfibers and vein-related tectonic processes (e.g., Mottram et al. 2020), but was applied prior to in situ U-Pb dating of low-U fault-related carbonate (e.g., Nuriel et al. 2017;Roberts et al. 2020).However, some fault gouges lack calcite and quartz veins developed by fluid interaction and therefore cannot be investigated by in situ U-Pb geochronology (i.e., Kemp et al. 2019).Most studies focus on the thermal effect on clay gouge during brittle deformation, as fault zones might reactivate numerous times over their complex geological history, from formation to last slip event and subsequent grain-size reductions during multiple faulting events might affect its isotopic signature (e.g., Zwingmann et al. 2019).Clay fault gouge associated with the Atera Fault were collected within this study for K-Ar dating from a well-exposed river section near Kawaue, Nakatsugawa city, Central Japan (Fig. 1).

Geological framework
Active inland faults are widespread in Japan and investigations are important for earthquake forecasting and seismic hazard assessments (Toda 2016).Comprehensive geological summaries of fault zones in Japan and the Geology of Japan, including the study area, can be found elsewhere (Moreno et al. 2016).The study area is located in Central Japan approximately 70 km NE of Nagoya (Fig. 1).The Atera Fault is one of numerous large strikeslip left-lateral active faults in Central Japan and extends for ca.55 km from Nakatsugawa City to Gero City in Gifu Prefecture in a general NW to SE trend (Takagi et al. 2005).As discussed by Niwa et al. (2009), the outcrop comprises fault breccia, clay-rich fault gouge, and cataclasite with significant asymmetry due to different protolith lithology.The Atera Fault zone outcrop record several stages of deformation which are based on a fieldbased and microstructural study by Niwa et al. (2009).The last recent displacement of the Atera Fault occurred during a large earthquake in 1586 (Toda et al. 1996;Niwa et al. 2009).
Central Japan has a long and complex deformation history documented by neotectonic processes and interactions between topography and erosion.Surface exposures enable the study of fault gouges which are often poorly preserved due to extended weathering.The ages of the basement and host rocks in the study area are well constrained and supplemented by a large set of regional thermochronological data.However, only a few studies on localized retrograde authigenic clay mineral formation in discrete fault zones in Japan have been reported in the literature (Yamasaki et al. 2013), although such brittle faults are widespread and similar to studies reported within the European Alps (Zwingmann and Mancktelow 2004;Zwingmann et al. 2010a).
The investigated fault outcrops are located in Kawaue, Nakatsugawa city (Fig. 1c; Niwa et al. 2009;2015).K-Ar dating of Atera Fault gouge is discussed in a report by Yamada et al. (1992) for limited < 2 μm clay fractions from different locations.The mineral compositions of the investigated gouges by Yamada et al. (1992) were determined by qualitative XRD analysis indicating mica clay mineral in addition to montmorillonite, quartz, kaolinite, chlorite, and feldspar.The reported ages range from 22.9 ± 0.8 Ma to 54.6 ± 1.8 Ma covering a large age range of 33.7 Ma for 7 samples.Based on a biotite K-Ar age from host rock granite at the same location, the 22.9 ± 0.8 Ma fault gouge sample suggests alteration.According to Niwa et al. (2009), the Kawaue basement comprises the Inagawa Granite, the Gero Ash-Flow sheet of the Nohi Rhyolite, the Naegi-Agematsu Granite, the Toki Sand and Gravel Bed, and the Ueno Basaltic Rocks (Fig. 2).The Inagawa Granite relates to the extensive Ryoke Granites, which intruded during the Cretaceous and its age was constrained by CHIME analyses in two localities with ages of 82.6 ± 1.8 Ma and 81.9 ± 1.4 Ma (2σ) which are identical within error (Suzuki and Adachi 1998).The Gero Ash-Flow Sheet represents a welded tuff with a younger 75.2 ± 6.2 Ma (1σ) Rb-Sr whole-rock isochron age (Shirahase 2005;Yamada and Koido 2005) and is part of the Late Cretaceous Nohi Rhyolite composite volcanism in Central Japan (Yamada et al. 2005).The Naegi-Agematsu Granite was dated at 67.2 ± 3.2 Ma (CHIME, 2σ) (Suzuki et al. 1994) and 71.3 ± 1.6 Ma (SHRIMP, 2σ) (Nakajima et al. 1993) and intrudes the Inagawa Granite and the Gero Ash-Flow Sheet.The terrestrial Toki Sand and Gravel Beds are of Upper Miocene to Lower Pleistocene age (Editorial Committee of CHUBU II 1988).
The Ueno Basaltic Rocks are extrusive volcano rocks above the Toki Sand and Gravel Bed.K-Ar ages of the Ueno Basaltic Rocks at Sakashita range from 1.41 ± 0.12 Ma (Uto and Yamada 1985), 1.68 ± 0.11 Ma (Shimizu et al. 1988), 1.64 ± 0.08 Ma to 1.54 ± 0.08 Ma (Ujike et al. 1992) with an overall average of 1.57 ± 0.05 Ma (2σ) based on error propagation indicating a Late Pliocene to Early Pleistocene origin.
The Atera Fault in the Kawaue area is separated by the welded tuff of the Gero Ash-Flow Sheet from the Naegi-Agematsu Granite (Fig. 2a, b).The petrography of the welded tuff is described in detail by Niwa et al. (2009) highlighting rhyolitic and phenocryst-rich material, with phenocrysts of quartz, plagioclase, K-feldspar, and biotite in a glassy groundmass supplemented by fine-grained chlorite and quartz.The feldspars remnants are often sericitized.Groundmass alteration results in clay minerals and calcite formation as secondary minerals.The Naegi-Agematsu Granite petrography highlights quartz, plagioclase, K-feldspar mineral phases, and biotite which is partially chloritized and fine-grained (Niwa et al. 2009;2015).
The Atera Fault outcrops consist of an ~ 42-m-wide fault zone, divided into numerous damage zones derived from welded tuff and granite, and a narrow 1.2-m-fault core (Fig. 2).
In the western part of the fault zone the dominant welded tuff damage zone is composed of a 10-m-wide zone of welded tuff and a 14-m-wide fault breccia.In the eastern part of the Atera Fault the granite contains a ~ 1.2-m-wide damage zone and a 5-m-wide granite cataclasite zone.The 1.2-m-wide fault core is clay-rich and contains fragments of welded tuff, granite, and basalt.According to Niwa et al. (2009) remaining basalt fragments correlate with the neighboring volcanic rocks of the Ueno Basaltic Rocks allowing to constraint on the maximum age of fault core formation.The fault core trends NW-SE to NNW-SSE with a subvertical dip.A sinistral sense of shear can be constrained as several rounded fragments have been rotated displaying asymmetric tails.These structural features of the fault core are consistent with the general structure of the Atera Fault (Niwa et al. 2009;2015).In addition to the 1.2-m-wide fault core, a 30-cm-wide fault overprints the granite cataclasite (Fig. 2c).The fault gouge is clay-rich and includes granite fragments.The trend of the gouge is subparallel to the 1.2-m-wide fault core with NW-SE to NNW-SSE strike, and a subvertical dip and slickenlines on the gouge plunge gently south (Niwa et al. 2015).

Sampling
Three clay-rich gouge samples were obtained from the main outcrop shown in Fig. 2c and selected based on the host rock mineralogy.An undeformed host rock sample of the welded tuff from the Gero Ash-Flow Sheet of the Nohi Rhyolite was collected about 3 km to the west (Tsukechi River 1) (Fig. 1c).Sample Atera 1 (3369) originates from the welded tuff fault breccia, sample Atera 2 (3370) from the center of the 1.2-m-wide fault core and sample Kawaue 1 (3367) from the 30-cm-wide fault gouge within the granite cataclasite.The location of the outcrop is shown in Fig. 1c, and GPS location data summarized in Table 1.The procedures applied in this study are based on Zwingmann et al. (2010aZwingmann et al. ( , 2011)).
About 50 to 100 g of fresh gouge and cataclasite material was collected from each sample site.Sample disintegration involved a gentle repetitive freezing and heating technique which avoids artificial reduction of K-bearing minerals rock components (micas or K-feldspars, see Liewig et al. 1987;Zwingmann et al. 2017).Separated grain-size fractions comprise < 2, 2-6 and 6-10 µm obtained according to Stokes' law.Fine grained < 0.1 µm fractions were separated using a high-speed centrifuge, but only limited amount could be obtained.For wholerock host rock analysis a small amount of the sample rock was crushed gently by hand in an agate mortar.Mineralogical characterization of the sample fractions involved

XRD
The mineralogy of the size fractions was determined by XRD on randomly oriented powder samples.Results for the < 0.1 µm fractions could not be obtained due to very limited sample amount and these fractions were mainly investigated by TEM (see below).Diffractograms of the samples were obtained and analyzed using a Rigaku Ultima IV with a two-axis goniometer, CuKα radiation, a Kβ filter, and a one-dimensional X-ray detector (D/teX Ultra), 40 kV, 30 mA, 0.5° and 10 mm divergence slits, and 8 mm anti-scattering slit.The samples were scanned over the range 2θ = 3-70° with a scan rate of 1° min −1 and a step size of 0.01°.The XRD results are listed in Table 2.

Petrography
The Atera 1 and 2 fault gouge and the host rock samples were investigated by thin section petrography.Petrographic information for sample Kawaue 1 (3367) from the 30-cm-wide fault gouge within the granite cataclasite can be found in Niwa et al. (2009).
A low-viscosity epoxy (Struers Epofix) was applied to the fault rock samples before cutting.As the fault rocks might contain clay minerals with different water content, synthetic hydrocarbon oil (ENEOS METAL WORK S) was used during thin section preparation, cutting and grinding.To mount fault sample chips to glass slides Petropoxy 154 was used and mounted sections reduced to ~ 35 µm thickness and polished.Freshly broken surfaces of sample chips were carbon-coated and examined by SEM in secondary electron mode using a HITACHI TD-1000 SEM equipped with an energy dispersive system X-ray analyzer (EDS; Horiba EMAX) at the JAEA Tono Geoscience Center, Japan.A JEOL F200 (CF-HR) TEM operating at 200 kV equipped with a Gatan OneView camera was used for comprehensive, grain-by-grain morphological characterization of selected < 0.1 µm clay fractions, and to investigate the grain-size distribution at the CMCA, UWA, Perth, Australia.For TEM investigation, samples were prepared by placing one drop of clay solution on a Cu TEM grid with a thin, amorphous carbon support film and drying in air.The composition of individual particles was investigated with a JEOL SD100WL 100 mm2 SDD (Silicon Drift Detector) EDS detector run through Oxford Instruments' AZtec software.

K-Ar dating
K-Ar dating follows standard methods given in detail by Dalrymple and Lanphere (1969) and carried out at CSIRO Energy, Perth, Australia, using a procedure similar to that described by Bonhomme et al. (1975).Potassium content was measured by ICP-OES (Agilent Technologies 725 ICP-OES) on ~ 20 mg of sample which was dissolved using HF and HNO 3 and diluted to 10-50 ppm K.The pooled error of K determination of all samples and standards is better than 2%.
For argon analyses the mineral fraction were fused within a vacuum line connected to a 38 Ar spike pipette.An online VG3600 mass spectrometer was used to measure the isotopic composition of the spiked Ar.The 38 Ar was calibrated against standard biotite GA1550 ( McDougall and Roksandic 1974).Special care was taken in the preparation of both K and Ar sample splits as clays are hygroscopic.For Argon analysis the sample splits were loaded into clean Mo-foil (Goodfellow molybdenum foil, thickness 0.025 mm, purity 99.9%), weighed and preheated to 80°C for 8 h to remove moisture, and reweighted using a Mettler AT20 balance and the measured dry weight used in the K-Ar age calculation.A desiccator was used for sample storage prior to loading to the argon sample holder.The mass discrimination factor was

XRD
The fault clay gouge mineralogy was obtained by XRD analysis.The representative XRD spectra of the < 2 µm fractions and the host rock powder are summarized in Fig. 3.The fault gouge samples comprise: smectite, illite, quartz, K-feldspar and plagioclase (Table 2).Amounts of smectite and illite increase with decreasing particle size and amounts of quartz, K-feldspar and plagioclase decrease with decreasing particle size.For sample Atera 1, the < 2 µm fraction contains the highest amount of illite with 27%, decreasing to 20% in the 2-6 µm and 12% in the coarsest 6-10 µm fraction.This trend is similar for the Atera 2 sample with illite enriched in the < 2 µm fraction with 30% decreasing to 17% in the 2-6 µm fraction.
The smectite content of both samples is decreasing in a similar manner for both fault gouge samples from 44% in the < 2 µm fraction, 22% in the 2-6 µm fraction to 13% in the coarse 6-10 µm fraction for Atera 1 and, 57% for the < 2 µm fraction to 43% in the 2-6 µm fraction for the Atera 2 sample, respectively.In both samples the K-feldspar content increases with grain size from 5 to 8% for Atera 1 fractions whereas, for the Atera 2 fractions, K-feldspar content is homogenous ~ 3 to 4%.The Kawaue 1 sample contains a high amount of illite at 33%, a low amount of smectite at 9% and a significant K-feldspar content of 11%.The Nohi Rhyolite host rock is rich in quartz, plagioclase, K-feldspar, and chlorite with no illite present.Illite polytype analysis results are inconclusive due to low amounts of sample and mixture of illite and feldspar in the < 2 µm and coarser fractions.

Optical microscope petrography
Optical microscopic observation confirms that all Atera Fault gouge samples are clay-rich with distinct foliation and a preferred alignment of minerals (Fig. 4).The Atera 1 sample shows significant development of microfractures in lithic fragments of welded tuff and mineral fragments of quartz and feldspar.Clay minerals are unevenly distributed and concentrated in larger fractures confirming previous studies by Niwa et al. (2009).
The Atera 2 sample originates from the main fault core and contains polymineralic fragments of welded tuff, granite, and andesite/basalt, and monomineralic fragments of quartz, plagioclase, and K-feldspar in a clay matrix.The fragments are subangular to subrounded.Fine fragments of well aligned phyllosilicates and quartz are observed in the clay matrix.Additional petrographic information of the Atera Fault zone can be found in Niwa et al. (2009).
Contrary to the Atera fault samples, the thin section petrography of the Nohi Rhyolite host rock indicates no

Table 3 K-Ar standards and airshot data
References HD-B1 (Hess and Lippolt 1994), LP6 (Odin et al., 1982) foliation and few microfractures.The host rock sample contains large ~ 1-3 mm quartz and plagioclase minerals as well as chloritized mafic minerals.The fine-grained groundmass is mainly composed of quartz and plagioclase.The quartz phenocrysts are often corroded and plagioclase mostly sericitized.

SEM results
SEM investigations reveal that smectite/illite minerals in the fault gouge zones document preferred orientation and exhibit well crystallized euhedral grain boundaries, related to an authigenic and synkinematic origin, compared to host rock samples as shown in Fig. 5 (Clauer and Chaudhuri 1995).However, no direct anastomosing foliation defined by illite/smectite (I/S) coatings around clasts of host rock can be observed.The observed gouge minerals are interpreted to be of synkinematic origin due to the preferred orientation.Irregular remaining host rock grains are surrounded by a finer grained matrix and might be a result of mechanical comminution in contrast to forms more typical of well-crystallized, authigenic clay minerals comprising fibrous, lath-shaped or platelet morphologies (Inoue et al. 1988;Pevear 1999).The petrographic observations are similar to a study focusing on retrograde faults from the Swiss Alps (Mancktelow et al. 2015).Features revealed by the SEM are different from the fault gouge and Nohi Rhyolite host rock.Euhedral mineral grains mostly composed of quartz and feldspar are filled tightly, without fragmentation.The SEM observed and investigated particles comprise: Al, Si, O, K and Fe as determined by EDS, and correspond to smectite/illite minerals as identified by XRD.

TEM results
TEM investigations of the selected < 0.1 µm clay fractions indicate distinct idiomorphic platy illite flakes, with well crystallized euhedral edges (Fig. 6), and minor electron dense dark particles with diffuse and irregular edges often interpreted as detrital in origin (Clauer and Chaudhuri 1995) (Fig. 6).The particles investigated by TEM EDS contain mainly: Al, Si, O, K and Fe and correspond to smectite/illite minerals as confirmed by SEM and XRD.Detrital clay particles are often volumetrically larger than neoformed illite particles and are characterized by different morphological shapes and often exhibit irregular edges, while neoformed illite particles

K-Ar data
Nine new K-Ar analyses covering one whole rock sample and eight clay fractions from < 0.1 to 6-10 µm, were obtained in this study (Table 4).The whole rock analysis As reported in numerous studies involving clay-size fractions, a consistent decrease in age with decreasing grain size can be observed (Fig. 7) (see summaries in Vrolijk et al. 2018;Tsukamoto et al. 2020).The ages of the < 2 µm fractions from the three fault gouge samples vary significantly: from 50.2 ± 1.2 Ma for Atera 1, to 47.0 ± 1.1 Ma for Kawaue 1 and 44.1 ± 1.0 Ma for sample Atera 2. Ages of the Atera 1 sample originating within the welded tuff breccia range from 40.6 ± 1.0 Ma for the finest < 0.1 µm fraction to 60.0 ± 1.4 Ma for the coarsest 6-10 µm fraction.Ages of the Atera 2 samples originating from the fault gouge yield an age of 41.8 ± 1.0 Ma for the finest < 0.1 µm fraction (identical within error of the Atera 1 < 0.1 µm sample) and a significantly younger age of 52.7 ± 1.2 Ma for the coarser 2-6 µm fraction.K concentration ranges from 2.74 to 4.98% and radiogenic 40 Ar content from 78.9 to 95.7% indicating no significant contamination with atmospheric 40 Ar.

Discussion
K-Ar dating of authigenic clays and illite in general, including fault gouge, requires careful investigations as discussed by Clauer et al. ( 2012) and Meunier and Velde (2004).The potential and limitations of K-Ar illite fault gouge dating is discussed elsewhere (see Vrolijk et al. 2018;Tsukamoto et al. 2020).The grain-size fractions of the separated authigenic illite are mixtures of illite particles formed at different times during growth and the formation history is normally investigated by dating a range of different grain-size fractions.Fault gouge illite forms in low-temperature retrograde environments governed by kinetic processes compared to equilibrium thermodynamics of mineral formation processes (e.g., Vrolijk and van der Pluijm 1999;Zwingmann and Mancktelow 2004).
We interpret the obtained ages as maximum ages as the authigenic illite within the finest separated particle size represents the most recently grown illite.
Increasing K-Ar ages with increasing grain size are illustrated in Figs.7 and 8 and related to increased contamination in the coarse fractions compared to the finer fractions as confirmed by XRD data (Table 2).The Atera 1 gouge sample originates from a heterogenous welded tuff breccia zone and the finest < 0.1 µm fraction yields an age of 40.6 ± 1.0 Ma and increases significantly by ~ 10 Ma to 50.2 ± 1.2 Ma for the coarser < 2 µm fraction.For the coarser, 2-6 µm and 6-10 µm fraction, the ages increase again by a steep 10 Ma to 59.1 ± 1.4 Ma and 60.0 ± 1.4 Ma, respectively.A similar but less severe age increase is observed for the Atera 2 gouge sample originating from the fault core with the finest < 0.1 µm fraction yielding an age of 41.8 ± 1.0 Ma increasing by about 2 Ma to 44.1 ± 1.0 Ma for the coarser < 2 µm fraction.The age increases again by a significant 10 Ma to 52.7 ± 1.2 Ma for the coarser 2-6 µm fraction compared to the < 0.1 µm fraction.For the Kawaue 1 gouge sample (3367) obtained from within the granite cataclasite zone, only a < 2 µm fraction was analyzed which yielded an age of 47.0 ± 1.1 Ma.This age is intermediate between the Atera 1 and 2, < 2 µm ages.The obtained whole rock host rock age is similar to the coarser 2-6 and 6-10 µm fractions from the fault zone samples (Fig. 7) and suggesting a thermal overprint or hydrothermal overprint might be associated with the fault initiation around 41.2 Ma and, might have partially influenced and reset the isotopic signatures of the Nohi Rhyolite host rock.The K-Ar ages of < 2 μm fractions from the Atera fault gouges from our study fall within the age range reported by Yamada et al. (1992), ranging from 38 to 55 Ma.Yamada et al. (1992) only reported a limited < 2 μm fraction age data without detailed description of the investigated fault rocks, whereas this study includes detailed gouge K-Ar age data from different size fractions including finer, < 0.1 μm fractions in a well-described fault zone outcrop of the Atera Fault.
The increase in ages of both fault gouge samples can be explained by increasing K-feldspar contamination (Table 2) and highlights that the most reliable isotopic ages for authigenic illite can be obtained from the finest size fractions (< 0.1 µm or finer) of gouge rock samples due to minimal contamination as discussed by Zwingmann et al. (2010a, 2011), Torgersen et al. (2015) and Viola et al. (2016).The average age of the two < 0.1 µm Atera gouge fractions yields a mean of 41.2 ± 0.7 Ma, which is significantly younger than the < 2 µm fraction ages of 47.2 ± 0.8 Ma, and correspond, within error, to the age obtained from the gouge sample within the granite cataclasite zone (Kawaue 1, 3367).The ages of all size fractions are relatively heterogenous suggesting that: (1) different contamination sources might be present in the fault gouge samples, or (2) the fault core sample experienced a more severe last faulting event suggesting an increasing effect of mechanical comminution toward the fault core (Fukuchi et al. 2014).

Implications for geology of Central Japan
As discussed by Zwingmann andMancktelow (2004, Zwingmann et al. (2010a;2010b) and Yamasaki et al. (2013) the K-Ar ages of fault gouge illite are bracketed by the zircon and apatite fission-track ages from host rocks (Fig. 8).Authigenic illite formation during brittle deformation is closely related to temperature ranges and occurs between the partial annealing zones of apatite, ranging from a low 60-120 °C (see references in Yamasaki et al. 2013), to a higher 170-390°C range for zircon  (Yamada et al. 1995).Similar to the limited brittle fault illite age data from the Toki Granite (Yamasaki et al. 2013), the age of main illite formation within the Atera Fault of ~ 47.1 Ma is younger than the reported fault gouge illite age data from the Nojima Fault ca.250 km to the SW (Zwingmann et al. 2010b).However, numerous Paleogene-Early Eocene illite fault gouge ages are reported in Japan along the Median Tectonic Line in the Shikoku and Kinki area about 200-350 km to the SW and the Mozumi-Sukenobu Fault about 100 km to the N (e.g., Yamasaki et al. 2013).The active faults are related to the Cretaceous granite margin and Jurassic accretionary complexes in the Inner Zone of SW Japan (Takagi et al. 2005).As discussed by Yamasaki et al. (2013) significant hydrothermal events accompany faults in SW Japan, ranging in age from ca. 60-50 Ma.Subsequent to clay formation hydrothermal fluid circulation is limited to the fault core and thus unable to form new clay minerals during fault reactivation.In our study, the age difference between the Nohi Rhyolite breccia zone sample (Atera 1), and the fault core sample (Atera 2), might indicate a later reactivation.This is similar to gouge samples in other case studies comprising the Deokpori Thrust, Korea, and the Nojima Fault, Japan, suggesting age variation with decreasing ages in grain size is caused by potential fault reactivation (Zwingmann et al. 2010b(Zwingmann et al. , 2011)).Illite age data reported by Yamasaki et al. (2013) for the nearby Toki Granite suggest that brittle faults were influenced by a later thermal event within the illite formation, the temperature ranging between 80 and 200 °C (see Zwingmann et al. 2019) yielding ages around 54 to 43 Ma and no illite ages related to an earlier 43 Ma (Paleogene-Middle Eocene) hydrothermal fluid event.
The thermochronological history and cooling curve of the study area for the nearby Toki Granite was summarized by Yamasaki et al. (2013) (Fig. 9A) and we adapted and approximated the thermochronological history and well-defined cooling curve for the Nohi Rhyolite and Naegi-Agematsu Granite (Fig. 9B).Ages are plotted from high closing to low closing temperature of minerals based on Yamada et al. (1992), Nakajima et al. (1993), Suzuki et al. (1994), Shirahase (2005), and Yamada et al. (2012).The age data of monazite CHIME, whole rock Rb-Sr, hornblende, biotite K-Ar, and ZFTA for host rocks from the Nohi Rhyolite and Naegi-Agematsu Granite are similar to the nearby Toki Granite.
The summarized complex thermochronological data support a model from Yamasaki et al. (2013), with a pluton intrusion at ~ 70 Ma and subsequent rapid cooling to ~ 300 °C.As outlined for the Toki Granite setting, the Naegi-Agematsu Granite experienced a similar rapid cooling after intrusion followed by slow cooling and associated with successive exhumation between the ZFTA and AFTA closure temperatures of zircon (~ 300°C) and apatite (~ 100°C), respectively, at a cooling rate of ~ 7-9°C/m.y. as suggested by Yamasaki and Umeda (2012).
The newly obtained K-Ar illite ages ranging from 40 to 60 Ma for the Atera 1 sample fraction and 41 to 53 Ma for Atera 2 fractions, plot on the 100-200°C temperature range of this cooling curve similar to the Toki Granite fault gouge formation (Fig. 9A).The < 0.1 µm fractions of both Atera fault gouge samples yield ages identical within error around 40 Ma, whereas the < 2 µm ages vary by 5 Ma from 50 to 44 Ma which might be caused by variable K-feldspar contamination.As the Kawaue 1 < 2 µm fraction contains the highest amount of K-feldspar contamination of ~ 11% but yields an intermediate age of 47 Ma, the overall impact of contamination might be negligible as the K-feldspar might be partially reset (Zwingmann et al. 2010a).The new Atera Fault gouge illite age data confirm that brittle deformation and thermal events are synchronous with the post-intrusive uplifting in this part of Central Japan (Yamasaki et al. 2013).The Atera Fault illite age data might provide new insights into an integrated regional-scale record of tectonic displacement in Central Japan and might be influenced by large-scale tectonic processes such as the Emperor sea mount kink around 50 to 46 Ma (Torsvik et al. 2017;Kubota et al. 2020) with fault initiation around 55 Ma In addition, similar illite K-Ar ages were reported for Nojima Fault gouges by Zwingmann et al. (2010b).

Conclusions
Complex fault cores and brittle fault rocks are rarely preserved in outcrops in Japan due to extensive weathering, erosion and/or vegetation cover.River outcrops and fresh exposures therefore provide important opportunities to investigate the fault zones and architecture of regionally important tectonic frameworks.The results presented from the Atera Fault in Central Japan form an internally consistent data set, supported by a large amount of geochronological data (AFTA, ZFTA, CHIME), demonstrating the application of gouge dating and provide new data constraining the timing of brittle deformation in Central Japan.The new illite age data from the well-exposed Atera Fault complete recently obtained, limited fault illite gouge age data from the Toki Granite (Yamasaki et al. 2013), and are, within error, identical.The age of the heterogenous welded tuff breccia zone (Atera 1) ranges from 40.6 ± 1.0 Ma to 60.0 ± 1.4 Ma (Eocene-Bartonian-Paleocene-Selandian).The age of the fault core gouge sample (Atera 2) ranges from 41.8 ± 1.0 Ma to 52.7 ± 1.2 Ma (Eocene-Lutetian-Ypresian).The finest < 0.1 µm fraction of the fault breccia and fault core gouge yield within error identical ages of around 41 Ma (Eocene-Bartonian).Similar to the nearby Toki Granite, the data constrain the timing of brittle deformation followed by a geothermal event in the Paleogene-Eocene.Brittle deformation and thermal event are related to the post-intrusive exhumation of Cretaceous granite.The new Atera Fault illite ages are comparable with results reported from other faults in the Inner Zone of SW Japan including the Median Tectonic Line, the Mozumi-Sukenobu Fault and the Nojima Fault (Shibata et al. 1989;Yamada et al. 1992;Zwingmann et al. 2010b) and might be related to large-scale tectonic processes, i.e., Emperor sea mount kink around 50 to 46 Ma (Torsvik et al. 2017;Kubota et al. 2020;Hu et al. 2022).

Fig. 1
Fig. 1 Study area.A Location map Japan with Median Tectonic Line (MTL) shown in red, N.A. North American Plate, B fault zones in Central Japan (modified from Takagi et al. 2005), C Atera Fault local geology

Fig. 2
Fig. 2 Outcrop details.A Atera Fault overview Kawaue area, B. studied outcrop sketch, C outcrop image and structural context of B (modified from Niwa et al. 2015).(1) to (4) in B and C show major lithological boundaries and correspond to each other

Fig. 8
Fig. 8 K-Ar age data in the geochronological context of the study area in Central Japan.Toki Granite host rock data based on Shibara and Ishihara (1979), Suzuki and Adachi (1998), Yamasaki and Umeda (2012), Yuguchi et al. (2019) and Yamasaki et al. (2013).B biotite, H hornblende, R Rb-Sr, C Chime (chemical U-Th total Pb isochron method), U U-Pb, K K-feldspar, S Shrimp U-Pb.Opening of Japan sea based on Hoshi et al. (2015)

Table 1
GPS outcrop location

Table 2
Cohen et al. (2013)1977)odically by "airshots" which are small amounts of air for40Ar/ 36 Ar ratio measurements.Blanks were systematically determined for the extraction line and mass spectrometer.Five international standards (3 HD-B1, 2 LP-6) and 4 airshots were analyzed during the course of the study (Table3).The 40 Ar/ 36 Ar values of the airshot yielded 295.63 ± 0.21 (1σ).The error for argon analyses is below 1.0% (2σ) based on the long-term precision of measurements of international Argon standards.The 40 K abundance and decay constants recommended bySteiger and Jäger (1977)were used for K-Ar age calculation.The age uncertainties consider the errors during sample weighing,38Ar/ 36 Ar and 40 Ar/ 38 Ar measurements, and K analysis.Ages are reported to the timescale ofCohen et al. (2013).
Semi-quantification based on integrated intensity ratio using reference of corundum, WR whole rock (Steiger and Jäger 1977)ger 1977)

Table 4 K
-Ar age data WR whole rock, timescale Cohen et al. (2013) K-Ar age versus grain-size plot.AFTA apatite fission-track ages, ZFTA zircon fission-track ages, WR whole rock