Increased Terrigenous Supply to the Pelagic Panthalassa Superocean Across the Carnian Pluvial Episode: A Possible Link With Extensive Aridification in the Pangean Interior

Abstract

In the Late Triassic, a global environmental change called the Carnian Pluvial Episode (CPE) emerged, causing major biological turnover. The CPE has been recognized by siliciclastic input to sedimentary basins, multiple carbon isotope perturbations, and climate proxies for humidification. The CPE is considered to have been associated with increased atmospheric pCO₂ from eruptions of large igneous provinces. However, the nature of this global environmental perturbation on the continents is still not well understood. Here we present a geochemical analysis of a pelagic deep-sea bedded chert sequence across the CPE in the Jurassic accretionary complex of Mino terrane, central Japan. Fluctuations in terrigenous material supply were reconstructed using Principal Component Analysis of major element compositions. The first principal component positively correlates with elements enriched in clay minerals such as Al₂O₃, whereas it negatively correlates with CaO, P₂O₅, and MnO, derived from apatite and manganese. A sudden increase in terrigenous supply was detected around the Julian/Tuvalian boundary, suggesting that CPE-related siliciclastic input also occurred in the abyssal plain environment. The terrigenous supply returned to the pre-CPE state in the Tuvalian. Since the terrigenous material supplied to the abyssal plain is thought to be derived from eolian dust blown from continental arid regions, the increasing terrigenous supply detected in the pelagic deep-sea chert succession may indicate extensive aridification. This result seems to conflict with the common view of the CPE as a humidification event. This contradiction possibly suggests that the extensive aridification occurred within the interior of the supercontinent Pangea, while hydrological circulation enhanced on the coastal region during the CPE.

1 Introduction

1.1 Carnian Pluvial Episode

The Carnian Pluvial Episode (CPE) was a time of global environmental perturbations with significant biological turnover both in the ocean and on land (e.g., Simms and Ruffell, 1989; Dal Corso et al., 2020). The CPE was initially recognized in many European sedimentary basins as drastic lithologic changes, such as the demise of carbonate platforms and increased siliciclastic input (e.g., Gattolin et al., 2015). Multiple negative shifts of δ¹³C are recognized across the CPE, suggesting a link with volcanic activity of Wrangellia LIPs (e.g., Dal Corso et al., 2018). Recent studies revealed that the intensification of the hydrological cycle also occurred not only in the western Tethys realm but in worldwide sedimentary basins (e.g., North China Plate; Lu et al., 2021; Boreal realm; Mueller et al., 2016; Gondwana region; Bernardi et al., 2018).

In contrast, the terrestrial records of the Pangean inland region are still controversial. Paleoclimate proxies, such as clay mineral profiles and paleosols, suggest only minor climate change with a generally arid to semi-arid background climate across the CPE from the United Kingdom (Baranyi et al., 2019a) and from Spain (Barrenechea et al., 2018), while data from the Germanic Basin and the Danish Basin did not show any humidity signal across the CPE (Visscher et al., 1994; Lindström et al., 2017). Enhanced continental weathering related to a more humid climate is only suggested for the early stages of the CPE in the Julian 2 (Baranyi et al., 2019b). Such observed discrepancies in the environmental signal across the CPE seem to be related to the contrasts in the climate of different regions and the limited age constraint of the terrestrial strata.

Abyssal deep-sea sediment provides a relatively continuous and consistent archive (e.g., Ikeda and Tada., 2014). An increase in the relative abundance of smectite among Fe-bearing minerals was detected in Japanese pelagic deep-sea chert across Julian/Tuvalian boundary (Nakada et al., 2014). Studies of Cenozoic sediments show that the terrigenous materials supplied to abyssal plains are mostly eolian dust generated from continental arid regions (e.g., Asahara, 1999). It is also noteworthy that a relatively thick, muddy silicious claystone interval (CS-1; Sugiyama, 1997) is observed within the smectite interval (Nakada et al., 2014). Along with the presence of CS-1, the thickness of intercalated mudstones drastically increased above the CS-1, and this can be correlated at least regionally (Sugiyama, 1997; Sugiyama et al., 2001). Because chert mainly consists of radiolarian tests and eolian mudstone (e.g., Hori et al., 1993), such an increase in mudstone may be related to increased eolian supply and/or decreased productivity (Sugiyama, 1997). The paleoenvironmental significance of these lithologic changes in the Carnian bedded cherts remains unclear.

In this study, the Carnian bedded cherts from the type section of Japanese radiolarian zones were investigated in order to elucidate the influence of the CPE on the Pangean interior. We seek to extract terrestrial environmental conditions by investigating the geochemistry of the intercalated mudstone of the bedded cherts, which concentrate terrigenous materials much higher than chert layers. The terrigenous content fluctuation was reconstructed throughout the Carnian by multivariate statistical analyses of major element composition. The increase of terrigenous materials flux implying reinforcement of eolian supply occurred. Our results thus firstly provide evidence that extensive aridification occurred in the inland of the Pangea across the CPE.

2 Geological Setting and Lithology Description

2.1 Geological Background

Early Triassic to Early Jurassic bedded chert is exposed in the Inuyama area, central Japan, as a part of the Mino-Tamba belt, Jurassic accretionary complex (Figure 1). The bedded cherts and the overlying Middle Jurassic clastic rocks are repeatedly exposed as a chert-clastic sequence truncated by bedding-parallel thrusts (Figure 1B; Yao et al., 1980; Matsuda and Isozaki, 1991). Since this bedded chert had an extremely slow sedimentation rate of 1–3 mm/kyr, and lacks coarse-grained terrigenous material for over 60 Myr, it is considered to have been deposited in a distal setting (Matsuda and Isozaki, 1991; Hori et al., 1993). The latitude of the depositional site was estimated as 12.3° ± 5.6°N or S and 19.0° ± 10.5°N or S for 237 Ma and 223.5 Ma, respectively based on the inclination of magnetic component of hematite (Figure 1A; Oda and Suzuki, 2000).

FIGURE 1

Bedded cherts are composed of rhythmic alteration of chert layers (cm scale) and intercalated mudstones. Chert layers are mainly composed of microcrystalline quartz and biogenic siliceous tests, whereas intercalated mudstones are mainly composed of clay minerals derived from eolian dust brown from the continental region (Hori et al., 1993). The relationship of chemical composition of chert layers and intercalated mudstones has been recognized as silica dilution by radiolarian test through the observation of negative correlation of SiO2 with other elements (e.g., Hori et al., 2000).

Bedded cherts are considered to be formed by a diagenetic segregation process as silica dissolution and transportation (Tada, 1991). The lithology of bedded cherts in the Inuyama area is mainly divided into two types: F type and B type which denote “fine grained quartz” and “biogenic clasts”, respectively (Imoto 1983; Sugiyama, 1997). F-type cherts show a few siliceous biogenic tests scattered in the fine-grained quartz matrix under microscopic observation with very thin (ca. 0.5–2 mm) intercalated mudstone layers and laminations. B-type cherts show concentrated biogenic tests under microscopic observation and are usually associated with thick (ca. 2–20 mm) intercalated mudstones. The difference in their formation mechanism may be related to changes in silica and terrigenous contents of the precursor sediment of bedded chert, formed by cyclic changes of sedimentary environments and/or diagenetic conditions (Sugiyama, 1997).

In this study, Section R and Section Q of Sugiyama (1997) in the Inuyama area are selected for detailed geological mapping and geochemical study in order to reconstruct paleoclimatic conditions of the entire Carnian stage (Figures 1, 2). Section R is the type section of Japanese radiolarian zone TR5A, which contains the most complete and unaltered sequence of the TR5A zone (Sugiyama, 1997). Section Q is located to the west of Section R and is one of the most complete and unaltered sequences of TR5B zone in Inuyama area (Sugiyama, 1997).

FIGURE 2

The depositional age of the Carnian bedded chert is constrained by radiolarian and conodont biostratigraphy, cyclostratigraphy, and carbon isotope stratigraphy (Sugiyama, 1997; Ikeda and Tada., 2014; Nakada et al., 2014; Yamashita et al., 2018; Tomimatsu et al., 2021). Sugiyama (1997) correlated interval (ca. 20 m thick) from the upper part of the radiolarian zone TR4B (Spongoserrula delhi) to the lower part of TR6A (Capnodoce–Trialatus) based on the first occurrence (FO) of the Carnian radiolaria Praeheliostaurus laevis, the FO of Carnian conodont Neogondolella foliata, and the FO of Norian radiolaria Capnodoce gilis in Sections N, R and Q in Inuyama area. Nakada et al. (2014) assigned the Julian/Tuvalian boundary within the Elbistanium gracile radiolarian zone between their samples R3-1 and R5-3. Tomimatsu et al. (2021) detected two negative excursions in carbon isotopes, NCIE-α and NCIE-β, in the TR5A zone within the Nicoraella? budaensis conodont zone in Section N, and correlated them to the global carbon isotopic excursion across the CPE.

2.2 Field Description

We conducted a detailed three-dimensional outcrop observation and constructed a bed-by-bed scale continuous lithostratigraphy to avoid missing or duplicating the strata (Figure 3). The lithological column was made from the lower part of Section R to the upper part of Section Q, in ascending order. The F-type red chert is dominant with a white chert bed WC-2 of Sugiyama (1997) at the base. About 900 cm above WC-2, the chert changes to ca. 300 cm purple laminated interval. Above the purple chert interval, ca. 200 cm B-type chert is intercalated with the ca. 20 cm siliceous claystone bed CS-1 of Sugiyama (1997) at the base. The CS-1 between F-type chert to B-type chert is widely recognized, not only in Inuyama area, but also in Hisuikyo area, which is 25 km north to Inuyama area (Sugiyama, 1997; Sugiyama et al., 2001).

FIGURE 3

3 Materials and Methods

3.1 Geochemical Analysis

Eighty-one intercalated mudstone samples were collected from Sections R and Q to cover the Carnian interval (Figure 2). Fresh samples were handpicked for geochemical analyses. After cleaning in an ultrasonic bath, the samples were crushed using a tungsten carbide mortar and pestle then pulverized in an agate ball mill. The samples were heated at 550°C for 4 h and 950°C for 2 h to decompose organic matter and carbonate material respectively prior to analysis of major elements and loss on ignition was measured.

Major elements (SiO₂, TiO₂, Al₂O₃, Fe₂O₃, MnO, MgO, CaO, Na₂O, K₂O, and P₂O₅) were analyzed by X-ray fluorescence (ZSX Primus II, Rigaku Co., Japan) at Waseda University. To minimize matrix effects and mineral effects, the fused disk method was used with a sample to lithium borate flux ratio of 1:10. Calibration lines were created using reference samples GSJ (JDo-1, JSd-1, JSd-2, JSd-3, JSI-2, JLk-1, JCh-1, JA-2, JA-3, JB-1a, JB-2, JB-3, JG-1a, JG-2, JG-3, JGb-1, JP-1, JR-1, JR-2). The accuracy of calibration lines for all elements was better than 0.6%. The sample R39, which is located on ca. 60 cm below the bottom of CS-1, was hand-picked and measured twice because it showed very high CaO and P₂O₅ concentrations.

Some mudstone samples were observed using scanning electron microscopy (SEM-EDS, Hitachi S-3400N, Japan) at the Waseda University with an accelerating voltage of 15 kV. Soft X-ray photography was employed using a Naomi NX digital X-ray Transmission System (RF Co., Ltd, Japan) with tube voltage of 50 kV and tube current of 40 mA at the Waseda University.

3.2 Multivariate Statistics

3.2.1 Hierarchical Cluster Analysis

Q-mode cluster analysis was conducted to classify the compositional behavior of major elements. Variation matrix (Aitchison, 1986), which consists of variances of all possible logratios among components, was used as a distance matrix. Ward algorithm, which hierarchically groups the cluster by minimizing the total within-cluster variance, was used for the clustering criteria.

3.2.2 Principal Component Analysis

To investigate the latent variables that underlie the multivariate data, principal component analysis (PCA) was performed on the major element chemistry. PCA extracts the direction of maximum variance of the input dataset. Prior to PCA, the compositional data of 10 major elements are transformed to centered logratio (clr) space to avoid constant-sum constraint during the calculation without causing an increase in the dimension of the dataset (Aitchison, 1982; Ohta and Arai, 2007). The equation for clr transformation of major element compositional data x = (SiO₂, TiO₂, … , P₂O₅) is as follows:

To obtain robust results in the presence of outliers, robust PCA (‘Hubert method’, Hubert et al., 2005) were used for the current study. In the current dataset, some samples show a high abundance of specific elements (e.g., CaO, P₂O₅), indicating the anomalously high concentration of some mineral(s). Since PCA captures the direction of the highest variance of the given dataset, avoiding outliers before analysis is indispensable for investigating the latent variables, if any, of the structure of the dataset. The minimum covariance determinant method was adopted for outlier detection on the current dataset and the correlation matrix was used for the PCA calculation. Shapiro-Wilk test was used to evaluate normality of the clr variables after removing detected outliers.

The R package “compositions” (van den Boogaart and Tolosana-Delgado, 2008) and “rrcov” (Filzmoser and Todorov, 2013) were used for the cluster analysis and principal component analysis calculations.

4 Result

4.1 Major Element Geochemistry

The result of major element composition of the intercalated mudstone samples and data normalized by PAAS (Post Archean Australian Shale; Taylor and McLennan, 1985) are shown in Table 1 and Figure 4, and Supplementary Figure S1. To avoid the false correlation effect of silica dilution, TiO₂ was used to normalize the major element concentrations. The elemental abundance of the intercalated mudstones relative to TiO₂ is generally similar to PAAS, except for CaO, Na₂O, MnO (Figure 4). The value of MnO/TiO₂ of the studied samples is higher than that of PAAS, whereas CaO/TiO₂ and Na₂O/TiO₂ are distinctively depleted.

FIGURE 4

The presence of apatite has a significant effect on the CaO and P₂O₅ concentration. Figure 5 shows the Back Scatter Electron (BSE) image, Energy dispersive spectrometry, and scatter plot of CaO and P₂O₅ of the intercalated mudstones. Apatite fragments were observed in SEM. The concentration of CaO and P₂O₅ shows a clear positive correlation along with the P₂O₅/CaO of 0.75, the ratio of apatite (e.g. Takiguchi et al., 2006).

FIGURE 5

To evaluate the possible provenance of the studied materials, provenance discrimination diagram of Roser and Korsch (1988) were adopted (Figure 6). Almost all samples are plotted in the “Quartzose sedimentary provenance”.

FIGURE 6

To identify groups of elements with the same behavior, Q-mode cluster analysis were performed on the bulk chemical composition of the studied mudstones (Figure 7). Results of a cluster ordination of the different variables point to the different grouping of elements associated with siliciclastics (SiO₂, K₂O, MgO, Al₂O₃, Fe₂O₃, TiO₂, and Na₂O) and apatite elements (CaO and P₂O₅) along with MnO, which solely consists of a cluster.

FIGURE 7

4.2 Principal Component Analysis

To investigate latent information of the bulk chemical composition of the studied samples, PCA was performed (Table 2; Figure 8). The p-values of the Shapiro-Wilk test are 0.35 for clr (SiO₂), 0.44 for clr (TiO₂), 0.14 for clr (Al₂O₃), 0.45 for clr (Fe₂O₃), 0.86 for clr (MnO), 0.70 for clr (MgO), 0.09 for clr (CaO), 0.59 for clr (Na₂O), 0.04 for clr (K₂O), and 0.84 for clr (P₂O₅).

FIGURE 8

PC1 accounts for 38.9% of the variance of the clr-transformed dataset. PC1 has positive loading from Al₂O₃, K₂O, Na₂O, MgO, SiO₂, TiO₂, and Fe₂O₃, whereas it has negative loading from MnO, P₂O₅, and CaO. Elements positively loaded on PC1 belong to the siliciclastic cluster of above mentioned Q-mode cluster analysis. On the contrary, the elements belonging to the other two clusters (the cluster of manganese and apatite minerals) show negative loadings on PC1.

Stratigraphic variations of PC1 and Al₂O₃/TiO₂ ratio are shown in Figure 9. The PC1 and Al₂O₃/TiO₂ show a weak positive correlation (R = 0.3). Al₂O₃/TiO₂ ratio increased across the CS-1, and then decreased across the B type/F type boundary but this trend is within the fluctuation throughout the studied section. The stratigraphic variation of PC1 along with the other elemental percent data is shown in Supplementary Table S1.

FIGURE 9

PC2, which is statistically orthogonal with PC1, accounts for 23.1% of the variance of the clr-transformed dataset. PC2 is characterized by high positive loading from P₂O₅, CaO, and Na₂O, and negative loading from MnO and Fe₂O₃.

5 Discussion

5.1 Characterization of Intercalated Mudstone in the Bedded Cherts

The depositional site of the studied bedded chert is located on the deep-sea plane, which was isolated from continental fluvial transportation by trench (Figure 1; Matsuda and Isozaki, 1991). Thus, terrigenous materials contained by the bedded cherts are considered to supplied from aeolian transportation (Asahara, 1999). The provenance discrimination diagram of Roser and Korsch (1988) suggests that the provenance of the studied samples is from “Quartzose sedimentary provenance” (Figure 6). This result is concordant with the supposed transportation mode of the bedded cherts, that is, the terrigenous materials were blown in from continental interiors (Hori et al., 1993).

The major elemental compositions of the studied intercalated mudstone samples reflect the varying proportions of lithogenic, biogenic, and hydrogenic components in the sediments. The results of Q-mode cluster analysis point to the different grouping of “manganese” (Mn), “apatite” (Ca, P), and “terrestrial” elements (SiO₂, K₂O, MgO, Al₂O₃, Fe₂O₃, TiO₂, and Na₂O) (Figure 7).

Mn in our samples is enriched relative to PAAS (Figure 4). MnO along with Fe₂O₃ shows large negative loadings in PC2. Mn hydroxides or oxides (e.g., MnO₂) are deposited rapidly in particulate form under oxic conditions, whereas under suboxic and anoxic conditions Mn is reduced to divalent (Mn²⁺) and forms soluble cations (Sholkovitz et al., 1992; Calvert and Pedersen, 1993). It is thus suggested that the enrichment of Mn may reflect an oxic depositional environment of this section. Alternatively, additional Mn may form under surface oxidized Mn-rich zones during diagenesis (e.g., Pedersen et al., 1986).

The variation of CaO and P₂O₅ in our dataset is controlled by apatite content because of the strong positive correlation of CaO and P₂O₅ with the P₂O₅/CaO ratio of 0.27–0.69 similar to apatite (Figure 5). Takiguchi et al. (2006) interpreted the main carrier of apatite in this succession was authigenic in origin and formed during early diagenesis based on their P₂O₅/CaO ratio lower than that of conodont apatite (∼0.74; e.g., Ruttenberg and Berner, 1993; Van Cappellen and Berner, 1991). The variation of CaO, P2O5, and MnO characterize PC2, which contains 23.1% of the total variance of the studied dataset. CaO and P₂O₅ show large positive loading whereas MnO shows negative loading, indicating the opposite influence of phosphate and manganese. This relation agrees with the chemistry variation under manganese redirection (Dickens and Owen, 1994; Betzler et al., 2016; Bialik et al., 2020).

The total concentrations of “manganese” and “apatite” elements are generally less than 5%, except for the extremely low PC1 samples at ca. 20 cm below the thick claystone CS-1 (Sample Ra39-1, Ra39-2, Ra40). These levels correspond to the unusually high values of CaO and P₂O₅ contents with “apatite” elements content of ca. 7–18%, which can be recognized as an event bed (Figure 5 and Figure 9). The P₂O₅/CaO ratio of these three samples are within 0.6–0.7, lower than that of conodont apatite (Takiguchi et al., 2006). This unusually high apatite concentration may be caused by a pelagic marine environmental shift which led to precipitate abiotic apatite in the sediments, but this topic remains outside the scope of this paper.

For “terrestrial” elements, although biogenic Si is the main component of radiolarian chert, Si in intercalated mudstone is considered to be mainly derived from siliciclastic materials due to migration of biogenic silica to adjacent chert during diagenesis segregation (e.g., Tada, 1991). During this diagenetic segregation, opal-A in siliceous fossils are extensively dissolved and reprecipitated as tiny opal-CT particles within the pore space mainly in higher dissolved Si area, and thus Si migrated from proto-mudstone intervals to adjacent proto-chert intervals (e.g., Tada, 1991). We cannot rule out contamination of biogenic Si in the intercalated mudstones, but the contribution should be negligible because only well-sorted mudstones were selected for the geochemical analysis.

5.2 Reconstruction of Terrigenous Fluctuations to the Panthalassa Superocean Through the Carnian

PC1 gradually increases below CS-1 from -1.5 to 1, except for the extremely negative values ca. 20 cm below the CS-1 (Sample Ra39-1, Ra39-2, Ra40). The PC1 score suddenly increases to 3.4 at the CS-1 horizon and remains constantly high values around 1.5 at the following ca. 1 m thick B-type bedded chert interval. The PC1 values decrease across the B type and F type boundary, and then seem to maintain a steady-state with higher variations than Julian to the end of Tuvalian.

Since PC1 has positive loading on all the “terrestrial” elements and negative loading on “manganese” and “apatite” elements, PC1 can be regarded as a latent variable that denotes the degree of relative contribution of terrigenous materials against marine-derived components. Besides, PC1 scores do not show a negative correlation with Al₂O₃/TiO₂ ratio (R = 0.3). Al₂O₃/TiO₂ ratio is known to reflect biological productivity on the abyssal fine-grained marine sediments since excess Al is scavenged by the biological pump (Murray and Leinen, 1996; Dymond et al., 1997). The studied bedded cherts are estimated to be deposited at the low latitude of abyssal deep-sea Panthalassa (Ando et al., 2001; Uno et al., 2015) thus the fluctuation of Al₂O₃/TiO₂ ratio of the current dataset is expected to reflect productivity. Therefore, the non-negative correlation observed between PC1 and Al₂O₃/TiO₂ suggests that the changes in PC1 are not mainly regulated by biological productivity fluctuations.

In addition, lithology of the bedded chert suddenly changes from F-type to B-type just above CS-1, with corresponding increased mudstone thickness (Sugiyama, 1997) and an increase in PC1 (Figure 9). Two possible factors that regulate the type of bedded chert (F type or B type) that is, the changes in burial flux of biogenic silica and/or siliciclastics (Sugiyama, 1997). Because of the increased Al₂O₃/TiO₂ ratio across the CS-1, a coincident increase in PC1 score may reflect increased terrigenous supply rather than decreased productivity (Figure 9).

Equally, considering the negligible contribution of “manganese” and “apatite” elements in samples other than Ra39 and Ra40, the increased mudstone bed thickness across CS-1 likely reflects increased flux of terrigenous material, probably eolian dust in the equatorial pelagic Panthalassa (Hori et al., 1993).

5.3 Continental Environmental Condition Inferred by Terrigenous Supply and Clay Mineral Content Fluctuations of Deep-Sea Bedded Chert

The Carnian Pluvial Episode has been characterized by abrupt siliciclastic input in Tethyan regions (Schlager and Schollnberger, 1974; Simms and Ruffell, 1989; Simms and Ruffell, 1990). This study suggests an increased terrigenous input also occurred in the abyssal deep-sea Panthalassa across the CPE. The increases in siliciclastic input at continental and coastal regions are considered to have been a consequence of enhanced hydrological cycling and consequent increased runoff. However, as discussed above, the terrigenous input in the abyssal plane reflects increasing eolian supply (Hori et al., 1993). Since this eolian dust is likely blown from arid continental areas, the increasing terrigenous supply detected in this study may imply the extensive aridification of the Pangean supercontinent.

This idea seems to contradict the existence of smectite across the CPE which was interpreted as evidence of continental humidification (Nakada et al., 2014). However, a megamonsoonal climate was known to prevail during the Triassic owing to the configuration of supercontinent-superocean (e.g., Kutzbach and Gallimore, 1989). Observational and climate modeling studies at the present day have indicated the humidification of tropical regions and aridification of subtropical regions under global warming (i.e., “Dry gets drier, wet gets wetter”; Held and Soden, 2006). The expansion of arid regions and associated aridification of previously humid regions could account for production of eolian dust containing smectite. The arid and seasonal influences around Carnian are detected in some Gondwana successions (Ruffell et al., 2016).

Expansion of arid regions could have also increased the emission of eolian dust with smectite across CS-1. Eolian deposits temporarily appeared in the Fundy (Canada), Iberian (Spain) and Argana (Morocco) basins during the Carnian, possibly near the end of the CPE, based on the limited biostratigraphy of terrestrial vertebrates and plant fossils (Arche and López-Gómez, 2014; Sues and Olsen, 2015). Interestingly, the onset of the smectite-rich interval slightly preceded the increased mudstone bed thickness across the CS-1 (Figure 10; Sugiyama, 1997; Sugiyama et al., 2001; Nakada et al., 2014). This temporal relation is consistent with the formation of smectite during the CPE and subsequent increase in eolian dust flux owing to the expansion of arid regions in the previously humid regions. Besides, weathering of mafic volcanic sources could also account for the smectite.

FIGURE 10

Increasing wind gustiness by global cooling and increased latitudinal temperature difference is also known as a main cause of increased eolian flux increase on glacial-interglacial timescales (e.g., McGee et al., 2010). Across the CPE, high-resolution sea surface temperature profile derived from oxygen-isotope of conodont apatite show a possible two-pulse warming event in the Nanpanjiang Basin in South China Block (Sun et al., 2016), while other low-resolution records from western Tethys show only one pulse (e.g., Hornung et al., 2007). In the South China Block, the onset of the CPE, namely NCIE-1, coincided with the first pulse of warming of ∼4°C, and followed by a short cooling phase across the Julian/Tuvalian boundary, and then by a prolonged and more intense second phase of warming (Figure 10; Sun et al., 2016). According to the biostratigraphy and carbon isotope correlation (Figure 10; Tomimatsu et al., 2021), the CS-1 level can be correlated with this short cooling phase across the Julian/Tuvalian boundary. Thus, cooling and associated wind gustiness, in addition to the expansion of arid regions, could have also enhanced the eolian dust emission recorded as the lithologic change in this study section (Sugiyama, 1997).

6 Conclusion

The cluster analysis and principal component analysis for the major element composition of intercalated mudstone within bedded chert were performed to investigate the continental environment fluctuation in the Carnian. The PC1 represents the terrigenous content fluctuation whereas PC2 reflects the relative contribution of manganese over apatite. The stratigraphic variation of PC1, coupled with an increase in mudstone thickness, suggests that terrigenous content in abyssal deep-sea sediments from Japan increased at around the Julian/Tuvalian boundary (i.e., across the CPE), and then return to pre-CPE level in Tuvalian. A delayed increase in dust flux relative to an increase in smectite content in deep-sea chert implied expansion of arid regions and associated aridification of previously humid regions. The appearance of eolian successions in low-latitude Pangea are consistent with this idea, but these age models have large uncertainty. This study thus provides deep-sea evidence for widespread aridification that occurred in the Pangean interior across the Julian/Tuvalian boundary, potentially related to CPE.

In addition, an event bed with very high apatite content was found immediately before the high terrigenous supply interval, implying environmental perturbation occurred in pelagic Panthalassa, potentially related to the CPE.

References

• 1

  AitchisonJ. (1982). The Statistical Analysis of Compositional Data. J. R. Stat. Soc. Ser. B Methodol.44, 139–160. 10.1111/j.2517-6161.1982.tb01195.xCrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 2

  AitchisonJ. (1986). The Statistical Analysis of Compositional Data/J. Aitchison. London: Chapman & Hall. Google Scholar

  • Google Scholar

• 3

  AndoA.KodamaK.KojimaS. (2001). Low-latitude and Southern Hemisphere Origin of Anisian (Triassic) Bedded Chert in the Inuyama Area, Mino Terrane, Central Japan. J. Geophys. Res.106, 1973–1986. 10.1029/2000jb900305CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 4

  ArcheA.López-GómezJ. (2014). The Carnian Pluvial Event in Western Europe: New Data from Iberia and Correlation with the Western Neotethys and Eastern North America-NW Africa Regions. Earth-Science Rev.128, 196–231. 10.1016/j.earscirev.2013.10.012CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 5

  AsaharaY. (1999). 87Sr/86Sr Variation in North Pacific Sediments: a Record of the Milankovitch Cycle in the Past 3 Million Years. Earth Planet. Sci. Lett.171, 453–464. 10.1016/s0012-821x(99)00158-2CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 6

  BaranyiV.MillerC. S.RuffellA.HounslowM. W.KürschnerW. M. (2019a). A Continental Record of the Carnian Pluvial Episode (CPE) from the Mercia Mudstone Group (UK): Palynology and Climatic Implications. J. Geol. Soc.176, 149–166. 10.1144/jgs2017-150CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 7

  BaranyiV.RostásiÁ.RaucsikB.KürschnerW. M. (2019b). Palynology and Weathering Proxies Reveal Climatic Fluctuations during the Carnian Pluvial Episode (CPE) (Late Triassic) from Marine Successions in the Transdanubian Range (Western Hungary). Glob. Planet. Change177, 157–172. 10.1016/j.gloplacha.2019.01.018CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 8

  BarrenecheaJ. F.López-GómezJ.De La HorraR. (2018). Sedimentology, Clay Mineralogy and Palaeosols of the Mid-carnian Pluvial Episode in Eastern Spain: Insights into Humidity and Sea-Level Variations. J. Geol. Soc.175, 993–1003. 10.1144/jgs2018-024CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 9

  BernardiM.GianollaP.PettiF. M.MiettoP.BentonM. J. (2018). Dinosaur Diversification Linked with the Carnian Pluvial Episode. Nat. Commun.9, 1499. 10.1038/s41467-018-03996-1PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 10

  BetzlerC.EberliG. P.KroonD.WrightJ. D.SwartP. K.NathB. N.et al (2016). The Abrupt Onset of the Modern South Asian Monsoon Winds. Sci. Rep.6, 29838. 10.1038/srep29838PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 11

  BialikO. M.ReolidJ.BetzlerC.EberliG. P.WaldmannN. D. (2020). Source Shifts to Periplatform Deposits during the Early to Middle Miocene in Response to Climatic and Oceanographic Forcing, Maldives, Western Indian Ocean. Palaeogeogr. Palaeoclimatol. Palaeoecol.559, 109969. 10.1016/j.palaeo.2020.109969CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 12

  CalvertS. E.PedersenT. F. (1993). Geochemistry of Recent Oxic and Anoxic Marine Sediments: Implications for the Geological Record. Mar. Geol.113, 67–88. 10.1016/0025-3227(93)90150-tCrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 13

  CappellenP. V.BernerR. A. (1991). Fluorapatite Crystal Growth from Modified Seawater Solutions. Geochimica Cosmochimica Acta55, 1219–1234. 10.1016/0016-7037(91)90302-lCrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 14

  Dal CorsoJ.MiettoP.NewtonR. J.PancostR. D.PretoN.RoghiG.et al (2012). Discovery of a Major Negative 13C Spike in the Carnian (Late Triassic) Linked to the Eruption of Wrangellia Flood Basalts. Geology40, 79–82. 10.1130/g32473.1CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 15

  Dal CorsoJ.GianollaP.NewtonR. J.FranceschiM.RoghiG.CaggiatiM.et al (2015). Carbon Isotope Records Reveal Synchronicity between Carbon Cycle Perturbation and the “Carnian Pluvial Event” in the Tethys Realm (Late Triassic). Glob. Planet. Change127, 79–90. 10.1016/j.gloplacha.2015.01.013CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 16

  Dal CorsoJ.GianollaP.RigoM.FranceschiM.RoghiG.MiettoP.et al (2018). Multiple Negative Carbon-Isotope Excursions during the Carnian Pluvial Episode (Late Triassic). Earth-science Rev.185, 732–750. 10.1016/j.earscirev.2018.07.004CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 17

  Dal CorsoJ.BernardiM.SunY.SongH.SeyfullahL. J.PretoN.et al (2020). Extinction and Dawn of the Modern World in the Carnian (Late Triassic). Sci. Adv.6, eaba0099. 10.1126/sciadv.aba0099PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 18

  DickensG. R.OwenR. M. (1994). Late Miocene-Early Pliocene Manganese Redirection in the Central Indian Ocean: Expansion of the Intermediate Water Oxygen Minimum Zone. Paleoceanography9, 169–181. 10.1029/93pa02699CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 19

  DymondJ.CollierR.McManusJ.HonjoS.ManganiniS. (1997). Can the Aluminum and Titanium Contents of Ocean Sediments Be Used to Determine the Paleoproductivity of the Oceans?Paleoceanography12, 586–593. 10.1029/97pa01135CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 20

  FilzmoserP.TodorovV. (2013). Robust Tools for the Imperfect World. Inf. Sci.245, 4–20. 10.1016/j.ins.2012.10.017CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 21

  GattolinG.PretoN.BredaA.FranceschiM.IsottonM.GianollaP. (2015). Sequence Stratigraphy after the Demise of a High-Relief Carbonate Platform (Carnian of the Dolomites): Sea-Level and Climate Disentangled. Palaeogeogr. Palaeoclimatol. Palaeoecol.423, 1–17. 10.1016/j.palaeo.2015.01.017CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 22

  HeldI. M.SodenB. J. (2006). Robust Responses of the Hydrological Cycle to Global Warming. J. Clim.19, 5686–5699. 10.1175/jcli3990.1CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 23

  HoriR. S.ChoC.-F.UmedaH. (1993). Origin of Cyclicity in Triassic-Jurassic Radiolarian Bedded Cherts of the Mino Accretionary Complex from Japan. Isl. Arc2, 170–180. 10.1111/j.1440-1738.1993.tb00084.xCrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 24

  HoriS. R.HiguchiY.FujikiT. (2000). Chemical Composition and Their Environmental Records of Bedded Cherts from Accretionary Complexes in Japan. Mem. Geol. Soc. Jpn.55, 43–59. (in Japanese). Google Scholar

  • Google Scholar

• 25

  HornungT.BrandnerR.KrystynL.JoachimskiM. M.KeimL. (2007). Multistratigraphic Constraints on the NW Tethyan “Carnian Crisis”. Glob. Triassic41, 59–67. Google Scholar

  • Google Scholar

• 26

  HubertM.RousseeuwP. J.Vanden BrandenK. (2005). ROBPCA: a New Approach to Robust Principal Component Analysis. Technometrics47, 64–79. 10.1198/004017004000000563CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 27

  IkedaM.TadaR. (2014). A 70 Million Year Astronomical Time Scale for the Deep-Sea Bedded Chert Sequence (Inuyama, Japan): Implications for Triassic-Jurassic Geochronology. Earth Planet. Sci. Lett.399, 30–43. 10.1016/j.epsl.2014.04.031CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 28

  ImotoN. (1983). “Chapter 22 Sedimentary Structures of Permian-Triassic Cherts in the Tamba District, Southwest Japan,” in Anonymous Developments in Sedimentology (Elsevier), 377–393. 10.1016/s0070-4571(08)70101-8CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 29

  KimuraK.HoriR. (1993). Offscraping Accretion of Jurassic Chert-Clastic Complexes in the Mino-Tamba Belt, Central Japan. J. Struct. Geol.15, 145–161. 10.1016/0191-8141(93)90092-oCrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 30

  KutzbachJ. E.GallimoreR. G. (1989). Pangaean Climates: Megamonsoons of the Megacontinent. J. Geophys. Res.94, 3341–3357. 10.1029/jd094id03p03341CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 31

  LindströmS.ErlströmM.PiaseckiS.NielsenL. H.MathiesenA. (2017). Palynology and Terrestrial Ecosystem Change of the Middle Triassic to Lowermost Jurassic Succession of the Eastern Danish Basin. Rev. Palaeobot. Palynol.244, 65–95. 10.1016/j.revpalbo.2017.04.007CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 32

  LuJ.ZhangP.Dal CorsoJ.YangM.WignallP. B.GreeneS. E.et al (2021). Volcanically Driven Lacustrine Ecosystem Changes during the Carnian Pluvial Episode (Late Triassic). Proc. Natl. Acad. Sci. U. S. A.118, 5118. 10.1073/pnas.2109895118CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 33

  MacArthurR. H. (1957). On the Relative Abundance of Bird Species. Proc. Natl. Acad. Sci. U.S.A.43, 293–295. 10.1073/pnas.43.3.293PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 34

  MatsudaT.IsozakiY. (1991). Well-documented Travel History of Mesozoic Pelagic Chert in Japan: From Remote Ocean to Subduction Zone. Tectonics10, 475–499. 10.1029/90tc02134CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 35

  MatthewsK. J.MaloneyK. T.ZahirovicS.WilliamsS. E.SetonM.MüllerR. D. (2016). Global Plate Boundary Evolution and Kinematics since the Late Paleozoic. Glob. Planet. Change146, 226–250. 10.1016/j.gloplacha.2016.10.002CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 36

  McGeeD.BroeckerW. S.WincklerG. (2010). Gustiness: The Driver of Glacial Dustiness?Quat. Sci. Rev.29, 2340–2350. 10.1016/j.quascirev.2010.06.009CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 37

  MuellerS.HounslowM. W.KürschnerW. M. (2016). Integrated Stratigraphy and Palaeoclimate History of the Carnian Pluvial Event in the Boreal Realm; New Data from the Upper Triassic Kapp Toscana Group in Central Spitsbergen (Norway). J. Geol. Soc.173, 186–202. 10.1144/jgs2015-028CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 38

  MurrayR. W.LeinenM. (1996). Scavenged Excess Aluminum and its Relationship to Bulk Titanium in Biogenic Sediment from the Central Equatorial Pacific Ocean. Geochimica Cosmochimica Acta60, 3869–3878. 10.1016/0016-7037(96)00236-0CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 39

  NakadaR.OgawaK.SuzukiN.TakahashiS.TakahashiY. (2014). Late Triassic Compositional Changes of Aeolian Dusts in the Pelagic Panthalassa: Response to the Continental Climatic Change. Palaeogeogr. Palaeoclimatol. Palaeoecol.393, 61–75. 10.1016/j.palaeo.2013.10.014CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 40

  OdaH.SuzukiH. (2000). Paleomagnetism of Triassic and Jurassic Red Bedded Chert of the Inuyama Area, Central Japan. J. Geophys. Res.105, 25743–25767. 10.1029/2000jb900267CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 41

  OhtaT.AraiH. (2007). Statistical Empirical Index of Chemical Weathering in Igneous Rocks: A New Tool for Evaluating the Degree of Weathering. Chem. Geol.240, 280–297. 10.1016/j.chemgeo.2007.02.017CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 42

  PedersenT. F.VogelJ. S.SouthonJ. R. (1986). Copper and Manganese in Hemipelagic Sediments at 21 °N, East Pacific Rise: Diagenetic Contrasts. Geochimica Cosmochimica Acta50, 2019–2031. 10.1016/0016-7037(86)90256-5CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 43

  RoserB. P.KorschR. J. (1988). Provenance Signatures of Sandstone-Mudstone Suites Determined Using Discriminant Function Analysis of Major-Element Data. Chem. Geol.67, 119–139. 10.1016/0009-2541(88)90010-1CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 44

  RuffellA.SimmsM. J.WignallP. B. (2016). The Carnian Humid Episode of the Late Triassic: a Review. Geol. Mag.153, 271–284. 10.1017/s0016756815000424CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 45

  RuttenbergK. C.BernerR. A. (1993). Authigenic Apatite Formation and Burial in Sediments from Non-upwelling, Continental Margin Environments. Geochimica Cosmochimica Acta57, 991–1007. 10.1016/0016-7037(93)90035-uCrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 46

  SchlagerW.SchöllnbergerW. (1974). Das Prinzip stratigraphischer Wenden in der Schichtfolge der Nördlichen Kalkalpen. Mitt. Geol. Ges. Wien. 66, 165–193. Google Scholar

  • Google Scholar

• 47

  ScoteseC. R. (2014). Atlas of Middle & Late Permian and Triassic Paleogeographic Maps. Paleomap Proj.3, 49. 10.13140/2.1.2609.9209CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 48

  SholkovitzE. R.ShawT. J.SchneiderD. L. (1992). The Geochemistry of Rare Earth Elements in the Seasonally Anoxic Water Column and Porewaters of Chesapeake Bay. Geochimica Cosmochimica Acta56, 3389–3402. 10.1016/0016-7037(92)90386-wCrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 49

  SimmsM. J.RuffellA. H. (1989). Synchroneity of Climatic Change and Extinctions in the Late Triassic. Geol17, 265–268. 10.1130/0091-7613(1989)017<0265:soccae>2.3.co;2CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 50

  SimmsM. J.RuffellA. H. (1990). Climatic and Biotic Change in the Late Triassic. J. Geol. Soc.147, 321–327. 10.1144/gsjgs.147.2.0321CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 51

  SuesH.OlsenP. (2015). Stratigraphic and Temporal Context and Faunal Diversity of Permian-Jurassic Continental Tetrapod Assemblages from the Fundy Rift Basin, Eastern Canada. Atl. Geol.51, 139–205. CrossRef Full Text | Google Scholar

  • Google Scholar

• 52

  SugiyamaK.KawakamiS.TakanoM. (2001). Aerial Photograph and Stratigraphic Correlation of the Triassic to Lower Jurassic Siliceous Claystone and Bedded Chert Units of the Inuyama and Hisuikyo Areas, Central Japan. Osaka, Japan: News of Osaka Micropaleontologists, 145–157. Google Scholar

  • Google Scholar

• 53

  SugiyamaK. (1997). Triassic and Lower Jurassic Radiolarian Biostratigraphy in the Siliceous Claystone and Bedded Chert Units of the Southeastern Mino Terrane, Central Japan. Bull. Mizunami Fossil Mus.24, 79–193. Google Scholar

  • Google Scholar

• 54

  SunY. D.WignallP. B.JoachimskiM. M.BondD. P. G.GrasbyS. E.LaiX. L.et al (2016). Climate Warming, Euxinia and Carbon Isotope Perturbations during the Carnian (Triassic) Crisis in South China. Earth Planet. Sci. Lett.444, 88–100. 10.1016/j.epsl.2016.03.037CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 55

  TadaR. (1991). in Compaction and Cementation in Siliceous Rocks and Their Possible Effect on Bedding Enhancement. Cycles and Events in Stratigraphy. Editor EinseleG.et al (New York: Springer-Verlag), 480–491. Google Scholar

  • Google Scholar

• 56

  TakiguchiT.SugitaniK.YamamotoK.SuzukiK. (2006). Biogeochemical Signatures Preserved in Ancient Siliceous Sediments; New Perspectives to Triassic Radiolarian Bedded Chert Compositions. Geochem. J.40, 33–45. 10.2343/geochemj.40.33CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 57

  TaylorS. R.McLennanS. M. (1985). The Continental Crust: Its Composition and Evolution. Oxford: Blackwellk, 312. Google Scholar

  • Google Scholar

• 58

  TomimatsuY.NozakiT.SatoH.TakayaY.KimuraJ.-I.ChangQ.et al (2021). Marine Osmium Isotope Record during the Carnian “Pluvial Episode” (Late Triassic) in the Pelagic Panthalassa Ocean. Glob. Planet. Change197, 103387. 10.1016/j.gloplacha.2020.103387CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 59

  UnoK.YamashitaD.OnoueT.UeharaD. (2015). Paleomagnetism of Triassic Bedded Chert from Japan for Determining the Age of an Impact Ejecta Layer Deposited on Peri-Equatorial Latitudes of the Paleo-Pacific Ocean: A Preliminary Analysis. Phys. Earth Planet. Interiors249, 59–67. 10.1016/j.pepi.2015.10.004CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 60

  Van den BoogaartK. G.Tolosana-DelgadoR. (2008). "compositions": A Unified R Package to Analyze Compositional Data. Comput. Geosciences34, 320–338. 10.1016/j.cageo.2006.11.017CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 61

  VisscherH.Van HouteM.BrugmanW. A.PoortR. J. (1994). Rejection of a Carnian (Late Triassic) “Pluvial Event” in Europe. Rev. Palaeobot. Palynology83, 217–226. 10.1016/0034-6667(94)90070-1CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 62

  YamashitaD.KatoH.OnoueT.SuzukiN. (2018). Integrated Upper Triassic Conodont and Radiolarian Biostratigraphies of the Panthalassa Ocean. Paleontological Res.22, 167–197. 10.2517/2017pr020CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 63

  YaoA.MtsudaT.IsozakiY. (1980). Triassic and Jurassic Radiolarians from the Inuyama Area, Central Japan. Osaka, Japan: Journal of Geosciences Osaka City University, 135–154. Google Scholar

  • Google Scholar