High-resolution and high-precision correlation of dark and light layers in the Quaternary hemipelagic sediments of the Japan Sea recovered during IODP Expedition 346

The Quaternary hemipelagic sediments of the Japan Sea are characterized by centimeter- to decimeter-scale alternation of dark and light clay to silty clay, which are bio-siliceous and/or bio-calcareous to a various degree. Each of the dark and light layers are considered as deposited synchronously throughout the deeper (> 500 m) part of the sea. However, attempts for correlation and age estimation of individual layers are limited to the upper few tens of meters. In addition, the exact timing of the depositional onset of these dark and light layers and its synchronicity throughout the deeper part of the sea have not been explored previously, although the onset timing was roughly estimated as ~ 1.5 Ma based on the result of Ocean Drilling Program legs 127/128. Consequently, it is not certain exactly when their deposition started, whether deposition of dark and light layers was synchronous and whether they are correlatable also in the earlier part of their depositional history. The Quaternary hemipelagic sediments of the Japan Sea were drilled at seven sites during Integrated Ocean Drilling Program Expedition 346 in 2013. Alternation of dark and light layers was recovered at six sites whose water depths are > ~ 900 m, and continuous composite columns were constructed at each site. Here, we report our effort to correlate individual dark layers and estimate their ages based on a newly constructed age model at Site U1424 using the best available paleomagnetic datum and marker tephras. The age model is further tuned to LR04 δ18O curve using gamma ray attenuation density (GRA) since it reflects diatom contents that are higher during interglacial high-stands. The constructed age model for Site U1424 is projected to other sites using correlation of dark layers to form a high-resolution and high-precision paleo-observatory network that allows to reconstruct changes in material fluxes with high spatio-temporal resolutions.


Background
The Japan Sea is one of the marginal seas on the western rim of northern Pacific that were formed more or less simultaneously during the Cenozoic. It is the most isolated marginal sea that is connected to neighboring oceans only with shallow (< 130 m in sill depth) and narrow (< 90 km in width) straits (Fig. 1). Partly because of its semi-isolation, the Japan Sea has its own deep water (called Japan Sea Proper Water), which ventilates the deeper part of the sea with the residence time of1 00 years (Gamo et al. 2014), and thus, the sea is occasionally called a miniature ocean (Gamo et al. 2014). The deep water below the water depth of~2000 m is highly homogenous with low and nearly constant temperature of 0.04~0.07°C, slightly low salinity of 34.07‰, low phosphorous of 1.9~2.2 μmol/kg, and high dissolved oxygen of 200~230 μmol/kg (Talley et al. 2006;Gamo et al. 2014). Consequently, the modern Japan Sea is one of the most oxic ocean basins in the world and relatively poor in nutrients. However, the situation was not necessarily the same in the geological past when sea level and the properties of the inflowing waters were different (e.g., Gallagher et al. 2015). We understand that the marginal sea discussed here is known by different names, including the "Japan Sea" and "East Sea". In this paper, we use "Japan Sea" given the order of authorship.

Previous studies
In the summer of 1989, JOIDES Resolution entered the Japan Sea and drilled at six sites during Ocean Drilling Program (ODP) legs 127/128 (Tamaki et al. 1990) to explore its tectonic and paleoceanographic histories. Based on the results of these legs, Tada et al. (1992) reported that Quaternary sediments in the deeper part of the sea are characterized by centimeter-to meter-scale alternations of dark, organic-rich layers and light, organic-poor layers. They further demonstrated that their deposition was synchronous basin-wide and the dark (and light) layers can be correlated in-between sites. Tada et al. (1999) demonstrated that alternation of dark and light layers reflects millennial-scale changes in the contribution of nutrient-rich East China Sea coastal water to the influx through the Tsushima Strait, which probably reflects changes in Yangtze River discharge caused by the changes in East Asian Summer Monsoon intensity in association with the Dansgaard-Oeschger Cycle (DOC) with dark layers deposited during interstadials (see also Tada 2005Tada , 2012. Consequently, the onset timing of the dark and light layer deposition and temporal changes in the amplitude and frequency of their alternation Fig. 1 A map showing locations of studied sites (red circles) together with unstudied sites (white circles) drilled during IODP Expedition 346 in the Japan Sea and northern East China Sea. Surface currents within the sea, water depths of each sites, and sill depths of each straits are also indicated Tada et al. Progress in Earth and Planetary Science (2018) 5:19 collected wide attention. Tada (2005) estimated the onset timing of the deposition of dark and light layers at approximately 1.5 Ma. However, since drilling was conducted in a single hole, the correlation gaps existed between cores and correlation among the sites was not perfect. Consequently, correlation went back only to 0.7 Ma (Tada et al. 1992). In addition, the age model was of lower resolution and precision since drilling was in a single hole and only limited paleomagnetic and biostratigraphic events were available at that time (e.g., Tada 1994).

Objective of this study
In order to specify the onset timing of deposition of dark and light layers, to examine the origin and the nature of their alternations, and to explore their temporal changes in frequency and amplitude, Integrated Ocean Drilling Program (IODP) Expedition 346 was conducted in the summer of 2013 and drilled at seven sites in the Japan Sea and two sites in the northern tip of the East China Sea (Fig. 1). To accomplish the objectives described above, high-resolution and high-precision correlation of sedimentary sequences among the sites within the sea is critical, especially for the Quaternary interval. Here, we present a result of group effort to correlate the dark and light layers in the Quaternary interval among six sites located in the deeper (> 900 m) parts of the sea where dark and light color contrast is distinct.

Studied sites
Seven sites were drilled during the IODP Expedition 346. Sites were selected so that they cover wide ranges of latitude, longitude, and water depths to form a "paleo-observatory network". With respect to the present water depth, Site U1427 is the shallowest site at~320 m water depth, Sites U1426 and U1430 are shallower depth sites of around 1000 m, Sites U1425 and U1423 are intermediate depth sites of slightly less than 2000 m, and Sites U1422 and U1424 are the deepest sites in excess of 2800 m water depth (Fig. 1). Sites U1426 and U1430 are on the bank, Site U1425 is on a terrace within the failed rift in the middle of the Yamato Rise, Sites U1423 and U1427 are on gentle slopes, and Sites U1422 and U1424 are on the basin floor. Sites U1427, U1424, and U1423 are under the influence of the first branch of the Tsushima Warm Current (TWC), Site U1426 is under the influence of the second branch, and Sites U1430 and U1425 are under the influence of the third branch of TWC, respectively. Site U1422 is under the influence of the Liman Cold Current. Sites U1422, U1423, and U1424 were influenced by sea ice at least during the last glacial period (Ikehara 2003). Sites U1427, U1426, and U1425, as well as Sites U1424, U1423, and U1422 can form latitudinal transects; Sites U1427, U1426, U1425, and U1424, as well as Sites U1423 and U1422 can form depth transects; and U1430, U1425, and U1424 can form a longitudinal transect, respectively.
Three to five holes are drilled at each site at least to the depth of~100 m so that we can construct a continuous stratigraphic column at each site by core splicing, at least for the Quaternary interval (Tada et al. 2015). Since Sites U1424 and U1426 are drilled approximately at the same locations as ODP Sites 794 and 798, respectively, where previous drilling (legs 127/128) drilled deeper into Miocene, it is possible to use the cores from these sites to complement stratigraphic information on the older part.

Inter-hole correlations and construction of composite columns
It is important to construct perfectly continuous composite columns at all the sites to conduct high-resolution and high-precision correlations of dark and light layers inbetween the sites. We conducted inter-hole correlation of all the recovered cores covering the lithological Unit I (Tada et al. 2015) at each site for this purpose. In addition, interhole correlation of cores enables us to precisely correlate off-splice parts of cores to the composite columns (spliced sections) so that we can use samples taken from off-splice parts of cores in the same way as those taken from the spliced parts and thus can drastically increase the amount of samples utilized for high-resolution studies.
Preliminary inter-hole correlation of cores was conducted onboard by stratigraphic correlators based on RGB (redgreen-blue), GRA (gamma ray attenuation), and NGR (natural gamma ray) profiles taken immediately after core retrieval (Tada et al. 2015). However, there was no time onboard to examine core photographs in detail to crosscheck the validity of inter-hole correlations. For this reason, there were possibilities that stratigraphic gaps and/or repetitions of the sequence caused by small faults mostly formed during drilling were overlooked. Also, drilling disturbance may affect physical properties such as GRA and NGR and even color in some cases, but they are easily overlooked just by examining profiles of physical properties and/or color especially if disturbances are minor. If the core sediments are homogeneous or heavily bioturbated, such small faults or disturbances are difficult to identify. However, in the case of the Quaternary part of the Japan Sea sediments, centimeter-to meter-scale dark and light layers serve as good markers to detect such minor faults and disturbances. We went back to high-resolution core photos of Unit I to identify minor faults, core disturbances, and artificial modifications such as plugging sponges in interstitial water (IW) sampling intervals and modified tie points to avoid such intervals. Details of procedures of inter-hole correlation, revision of tie points, and construction of composite sections are presented in Irino et al. (in press), and we adopt their results (revised composite sections) in this study.

Inter-site correlations
Inter-site correlations of dark and light layers are conducted based on visual inspection of composite columns of the six sites. First, we make a continuous composite column at each site using high-resolution core photographs (Irino et al., in press) and pick up relatively thick and distinct dark layers as marker layers. We set the bottom of marker dark layers as correlation tie points since many of the dark layers are laminated especially at their bottoms; thus, the effect of bioturbation is minimal (Watanabe et al. 2007). We use a composite column of Site U1424 as a reference column since its preliminary age model (Tada et al. 2015) suggests that the linear sedimentation rate (LSR) is more or less constant throughout the Quaternary and no significant sedimentation break exists. In addition, magnetostratigraphic events are best identified (Xuan et al., 2018, personal communication) and the largest number of tephra layers is found at this site. Next, we visually correlate marker dark layers at Site U1424 with those at other sites from the top of the column downward. We use marker tephra layers (Sagawa et al., in press) and magnetostratigraphic events (Xuan et al., 2018, personal communication) as guides to correlate marker dark layers among the sites. After inter-site correlation is completed, photographs of composite columnar sections of each sites are adjusted in their lengths at the core top and the Brunhes/Matsuyama boundary so that we can visually check the appropriateness of inter-site correlation of marker dark layers (Additional file 1: Figure S1).

Construction of the age model
The master age model of Unit I is constructed for Site U1424 where sedimentation seems continuous with relatively constant LSR, magnetostratigraphic events are best defined, and the largest number of tephra are identified among the six deeper sites drilled in the Japan Sea during Exp. 346. We use the revised magnetostratigraphy of Xuan et al. 2018 (personal communication) and tephrochronology of Sagawa et al. (in press) and Satoguchi and Nagahashi (2012) to construct an age model, which we call a tephra-p-mag-based age model (Fig. 2). The datum used to construct the age model is listed in Table 1. A tephra-p-mag-tuned age model is further tuned to LR04 age model (Lisiecki and Raymo 2005) by visually correlating GRA and NGR composite profiles of U1424 to LR04 δ 18 O stack of benthic foraminifers. It is well known that GRA decreases with the increase in diatom contents because diatom frustules have high internal porosity (Tada and Iijima 1983). NGR also decreases with the increase in diatom contents because diatoms dilute detrital material. In case of the Quaternary sediments of the Japan Sea, diatom abundance strongly links to glacio-eustatic sea level because surface productivity of the sea is controlled by nutrient influx through the Tsushima Strait, which in turn has been modulated by glacio-eustatic sea level changes (Tada et al. 1992(Tada et al. , 1999Tada 2012). Because the Tsushima Strait is only 130 m deep and its width and depth increase with the rise in sea level, the influx of the TWC is expected to increase drastically with the rise in sea level. Thus, diatom abundance can be used as a sea level proxy.
Tuning of U1424 age to LR04 age is conducted for the last~3 m.y. based on visual correlation of the GRA and NGR composite profiles with LR04 δ 18 O stack profile. We put priority on GRA over NGR for correlation because GRA data has higher resolution than NGR. We then cross-check the correlation by comparing LR04 stack profile with NGR profile. For most of glacialinterglacial cycles, GRA and NGR profiles show patterns similar to each other without phase shift and thus confirm the correlation based on GRA. For several intervals where the GRA profile does not show a clear pattern similar to LR04 δ 18 O stack profile, NGR profile is used for correlation. The correlation of GRA and NGR with LR04 δ 18 O stack is shown in Additional file 2: Figure S2 and CCSF-D_patched depth and age for each correlation point are listed in Table 2 where CCSF-D_patched depth is the depth scale for the revised composite section. A resulted LR04-tuned age model is shown in Fig. 3.

Age estimation of marker dark layers
In this study, we first constructed perfectly continuous composite columns for the entire Quaternary sequences at the six deeper sites studied, and using these composite columns, we correlated distinct dark layers among the sites Page 4 of 10 for the interval covering the last 1.5 My corresponding to the lithological Subunit IA (Additional file 1: Figure S1). Then, we constructed a high-resolution and high-precision age model at Site U1424 based on correlation of GRA and NGR with LR04 δ 18 O stack (Additional file 2: Figure S2). By using this age model, it is possible to estimate the ages of individual marker dark layers and thus project their ages to other sites.

Coding of dark layers
As is described above, dark layers are useful markers to make a high-resolution correlation between sites and assign precise ages to the sedimentary sequences in the deeper part of the Japan Sea. Consequently, it is useful to give them code numbers so that we can easily identify dark layers of our interest and estimate their ages. As to the last glacial period, thinly laminated (TL) layer numbering was conducted in the way so that TL numbers agree with interstadial numbers of DOC (Tada et al. 1999;Dansgaard et al. 1993). No numbering has been made on dark layers deposited before Marine isotope stage (MIS) 6 except by Tada et al. (1992) who made tentative naming of distinct dark layers at ODP Sites 794, 795, and 797. However, their naming was not systematic and went back to only ca. 0.7 Ma. In this study, we conduct numbering of the marker dark layers in Subunit IA covering the last 1.45 Ma. We first identify relatively thick and distinct dark layers as marker layers for the purpose of inter-site correlation as is described above. Then, using these marker dark layers and tephra layers as guides, we further identify thinner but distinct dark layers that can be correlated among the deeper six sites and can be used as second order marker layers. The results are listed in Additional file 3: Table S1.  The relatively thick and distinct dark layers occur every several meters and they tend to be associated with quasicyclic changes in amplitude and frequency of sediment lightness (L*) in Subunit IA, which seems to be associated with glacial-interglacial sea level changes as will be discussed later. In brief, thicker dark layers tend to occur during glacial maxima especially when sea level was lowered by more than 100 m and euxinic deep water emerged due to development of low-salinity surface water (Tada et al. 1999;Kido et al. 2007;also Seki 2017). Thick dark layers accumulated under euxinic deep water were found during MIS 2, 6, 10, 12, 16, 20, and 22 (Seki 2017). Such dark layers are overlain by thick light layers which were deposited during interglacial maxima. These characteristics are less clear before MIS 22 (~0.9 Ma) when thick diatomaceous dark layers frequently deposited during interglacial maxima.
We define a cycle of such several meter-scale change in amplitude and frequency of lightness as starting from the top of a relatively thick and distinct dark layer downward to the base of the relatively thick light layer so that one cycle basically corresponds to one glacial-interglacial cycle. From the core top downward, cycle number starts from 0 and increases downward. Because we tune changes in sea level proxies such as GRA and NGR to LR04 δ 18 O stack so as to construct high-resolution age model for Unit I as described above, and because one cycle basically represents one glacial-interglacial cycle, there is a relation between a cycle number and MIS number. Namely,  Lisiecki and Raymo (2005). All tuned points fall in an age versus depth relationship based on marker tephras and magnetostratigraphic events within uncertainties where X is a cycle number and Y is an MIS number when the thick dark layer at the top of cycle X deposited. Exceptions are cycle 0 that corresponds to the Holocene when no thick dark layer exists at the top, cycle 1 whose top thick dark layer corresponds to MIS 2 and thus X = Y/2, cycle 20 whose thick and dark layer deposited during MIS 45 (interglacial period) and thus X = Y/2 − 2.5, and cycle 21 whose thick dark layer deposited at MIS 46 and thus X = Y/2 − 2. A relatively thick and distinct dark layer at a top of cycle X is named Dark Layer (DL)-X-1 (Additional file 1: Figure S1; Additional file 3: Table S1). Every cycle contains several thinner but still distinct dark layers that are possible to correlate in-between sites. We put numbers 2, 3, … on them in descending order and described as suffixes to the cycle number. Namely, Zth dark layer in cycle X is DL-X-Z. In Additional file 3: Table S1, we list ages at the base of each dark layer based on tephra-p-mag-tuned and LR04-tuned age models for the last~3 My. We also list depths and ages at the top of marker dark layer that define the top of each cycle, although it should be noted that their depths and ages may be less well defined and could be slightly diachronous due to the result of bioturbation.
Temporal and spatial changes in linear sedimentation rates (LSRs) in the Japan Sea during the last 1.45 My As is described in the previous section, we estimated the age at the base of each marker dark layer at Site 1424 and projected its age to other sites based on inter-site correlation described above (Table 2). By using these data, it is possible to calculate linear sedimentation rates (LSRs) at all the six sites with high time resolution (Additional file 4: Table S2). Additional file 5: Figure S3 shows temporal changes in LSRs at all sites. LSRs are more or less constant at4 and~3 cm/ky at Sites U1422 and U1424 (deeper sites); gradually increasing from~3 to~10 cm/ky, from~2 to~7 cm/ky, and from~2 to~8 cm/ky at Sites U1423, U1425, and U1430 (intermediate depths sites); and relatively high and more or less constant at~10 cm/ky at Site U1426 (shallow site), respectively. Exceptions are between~850 and~1200 ka when LSRs are slightly higher and variable and after~250 ka when LSRs increase upward. At Site U1430, extremely low LSR is observed at~800 ka.
In general, LSR tends to be lower at deeper sites, which probably reflects dissolution of carbonate (and biogenic silica to some extent). The upward increase in LSRs after~250 ka is probably because near-surface sediments are less compacted. The gradually increasing trend of LSR at intermediate depth sites can be partly due to lesser degree of compaction, but the increasing trends started from deeper parts. The upward increasing trend of LSRs could be explained by the upward increase in biogenic carbonate content due to deepening of calcium carbonate compensation depth (CCD) that crossed the intermediate depths of 1000 to 2000 m at~1.0 Ma. Variable LSRs between~850 and~1200 ka could be related to large amplitude changes of CCD during transition from shallower CCD before 1200 ka to deeper CCD after 850 ka. Mass accumulation rate (MAR) calculation of carbonate at these sites based on high-resolution XRF core scanner and gamma ray attenuation (GRA) porosity data combined with LSR data will be necessary to test this possibility.

Establishment of paleo-observatory network
Reconstruction of material cycling such as carbon, sulfur, and phosphorous within the earth surface is critical to understand dynamics and controlling mechanisms of climatic changes in the past and preparing for the future. Especially, changes in centennial to millennial time scales are of prime importance since they are beyond the coverage of observational records, but still, climatic changes of such time scales are relevant to the society. Spatial coverage is also important since spatial changes in mass flux often give clues to understand the processes within the system of interest and their controlling factors. However, it is generally difficult to correlate geological records from different localities with high time precision and resolution. In this respect, Quaternary hemipelagic sediments of the Japan Sea provide a rare opportunity to study spatiotemporal changes of mass flux with high time precision and resolution. Namely, marker dark layers and tephra layers can be correlated with the time precision of less than 100 years, and approximately 250 marker layers are defined during the last 1.45 Ma, providing~6 ky resolution in average. These time precision and resolution are an order of magnitude higher than those for other correlation methods such as biostratigraphy, paleomagnetism, and radiometric dating. In addition, dry bulk density of the sediments can be estimated from GRA data measured onboard. Thus, it is possible to reconstruct MARs with an average time resolution of~6 ky at the six deeper IODP sites that cover the area more than a half of the Japan Sea. Moreover, we can calculate MAR of individual chemical element and/or mineral by combining high-resolution MAR data with high-resolution chemical and/or mineralogical analyses data.
In this study, we propose to utilize the six IODP sites in the Japan Sea as paleo-observatory network since highprecision age data and high-resolution LSRs and GRA data are available; the sites cover wide latitude, longitude, and water depth ranges; proxy data at each site can be