Changes in surface water masses in the northern East China Sea since the Last Glacial Maximum based on diatom assemblages

We investigated diatom assemblages for the last 19 ka in a piston core KY07-04 PC01 collected from the northern Okinawa Trough, East China Sea. Most of the diatom taxa in the Core KY07-04 PC01 were subtropical to temperate marine species throughout the core. Pronounced turnover of major diatom groups from periphytic to planktonic taxa occurred during deglaciation between 17 and 14 ka ago. During the Last Glacial Maximum (LGM) to early deglaciation, the representative taxa were Paralia spp., a periphytic diatom abundant in continental shelves of the modern East China Sea covering the Chinese Coastal Water. Relative abundances of Paralia spp. rapidly decreased during 15–14 ka ago, corresponding to sea-level rise, displaced by planktonic taxa such as Thalassionema. Thalassionema nitzschioides varieties living in subtropical water increased since 12 ka and have reached the present level at 8 ka ago, suggesting strengthening the influences of Kuroshio water in the northern East China Sea during the early to middle Holocene. Warm planktonic diatom taxa such as Nitzschia bicapitata also gradually increased from 12 to 8 ka ago. Near the core top, a sudden increase in N. bicapitata might relate to the intensifying Kuroshio during the latest Holocene.

(Changjiang) discharge into the East China Sea, providing much fresh water and suspended particles. The Yellow River is the second largest river globally in suspended sediment load (Milliman and Meade 1983). The Yangtze River supplies about 80% of the total discharge of freshwater from rivers flowing into the East China Sea, playing a critical role in sea-surface salinity distribution of the continental shelves of the East China Sea (Ichikawa and Beardsley 2002).
Bathymetry constrains the distribution of surface water masses and oceanic fronts in the East China Sea (Ichikawa and Beardsley 2002). The Kuroshio is a western boundary current of the North Pacific, playing an essential role in poleward heat transport (Qiu 2019). The Kuroshio flows northward and enters the East China Sea through the channel between Taiwan and Yonaguni-Jima. In the East China Sea, the Kuroshio path follows closely

Progress in Earth and
Planetary Science

Introduction
The East China Sea is a marginal sea of the western North Pacific bordered by the Kyushu, the Ryukyu Islands, and Taiwan. The East China Sea connects with the Japan Sea via the Tsushima Strait (~ 130 m sill depth) and the South China Sea via the Taiwan Strait (less than 100 m in depths). The East China Sea is composed of two areas: continental shelves (< 200 m water depth) in the northwest, including the Yellow Sea, and the Okinawa Trough (1000-2200 m water depth) in the southeast (Fig. 1). The Yellow River (Huanghe) and the Yangtze River along the steep continental slope and changes the direction to the east around ~ 30°N, flowing out to the Philippine Sea through the Tokara Strait. As a result, warm, saline, and nutrient-depleted Kuroshio water occupies the southern and central Okinawa Trough. The Tsushima Warm Current flows northward and enters the Japan Sea through the Tsushima Strait. The Taiwan Warm Current flowing northward through the Taiwan Strait and onshore Kuroshio intrusion across the East China Sea shelf break is the source of the Tsushima Warm Current (Isobe 2008). This warm water prevails in the northern Okinawa Trough. The Chinese Coastal Current flows southward along the Chinese coast on the continental shelves of the western East China Sea, mixing the Yangtze River discharge (Ichikawa and Beardsley 2002). Thus, cold, less saline, and nutrient-rich Chinese Coastal water distributes in the continental shelf regions of the East China Sea. The East China Sea has experienced marked geographic changes with the sea-level rise since the LGM. During the LGM between 21 and 17 ka ago, sea-level was estimated to be − 125 to − 130 m (Yokoyama et al. 2018), leading to aerial exposure of most East China Sea continental shelves, closure of the Taiwan Strait, and constriction of the Tsushima Strait. Such major geographic changes must have led to changes in the distribution and circulation of surface water in the East China Sea. There is a line of evidence for low sea-surface temperature (SST) widely observed over the glacial East China Sea (Ijiri et al. 2005;Sun et al. 2005;Chang et al. 2009a;Yu et al. 2009;Chen et al. 2010;Kubota et al. 2010Kubota et al. , 2015bXu et al. 2013;Ujiié et al. 2016;Matsuzaki et al. 2019), expansion of the Chinese Coastal Water (Xu and Oda 1999), advancement of river mouths of the Yellow River (Xu et al. 2014;Zheng et al. 2016;Zhao et al. 2017Zhao et al. , 2018Beny et al. 2018) and the Yangtze River (Dou et al. 2010a(Dou et al. , b, 2012. Most previous studies indicated weakened the Kuroshio Current during the LGM. At the same time, Kawahata and Ohshima (2004) suggested an intrusion of the Kuroshio Current into the northern East China Sea during the LGM based on the occurrence of pollen Phyllocladus, celery pines originated in New Guinea and Southeast Asia. The glacial path is controversial as reviewed by Gallagher et al. (2015) whether the southward migration occurred (Ujiié et al. 2003;Su and Wei 2005;Kao et al. 2006;Dou et al. 2010aDou et al. , b, 2012Dou et al. , 2016Lim et al. 2017;Li et al. 2018;Ma et al. 2021) or the path remained similar to the modern position (Lee et al. 2013;Zheng et al. 2016;Vogt-Vincent and Mitarai 2020).
Diatoms are unicellular photosynthetic eukaryotes producing biogenic opal frustules found in nearly all marine and freshwater habitats. They are one of the most major primary producers playing essential roles in the carbon and silicon cycles of the ocean (Smol and Stoermer 2010). A great diversity of species is living by sensitively responding to the aquatic environmental changes. In the modern East China Sea, characteristic diatom assemblages are dwelling in response to the water masses changes and the with pronounced seasonality (Asaoka 1975(Asaoka , 1980Oh and Koh 1995;Furuya et al. 1996;Chiang et al. 1999;Ishikawa and Furuya 2004;Guo et al. 2014;Li et al. 2014;Chen et al. 2019b). Although the East China Sea is befitting for reconstructing changes in surface water masses, a few paleoceanographic studies employing diatom assemblages are available. Chang et al. (2009b) presented a diatom assemblage record in Core MD01-2404 from the central part of the Okinawa Trough for the past 100 ka at millennial timescales. The diatom record in Core MD01-2404 seemed to be synchronous with the millennial-scale climate changes in the high-latitude Northern Hemisphere during the glacial-deglacial periods. Tanimura et al. (2002) employed a periphytic diatom, Paralia sulcata, in Core KH82-4 St-14 from the northern Okinawa Trough to discuss the expansion and retreat of the Chinese Coastal Water since the LGM with centennial timescales. They found abundant Paralia sulcata during the LGM to early glacial termination, which suggested significant expansion of the Chinese Coastal Water during then. Li et al. (2012Li et al. ( , 2015 investigated diatom assemblages in Core MD05-2908 from the southern Okinawa Trough at decadal resolutions for the past six ka. These studies have revealed that diatom assemblages are helpful to reconstruct past surface environmental changes of the East China Sea. However, a high-resolution whole assemblage record since the LGM has not yet been available in the East China Sea. Unfortunately, these previous studies (both modern and paleo) presented minimal taxonomic information and photos using a light microscope (LM). Only Li et al. (2014) presented scanning electron microscope (SEM) images for several Thalassiosira taxa. Page 3 of 17 Shirota et al. Progress in Earth and Planetary Science (2021)  Page 4 of 17 Shirota et al. Progress in Earth and Planetary Science (2021) 8:66 We present the whole diatom assemblages record and their LM and SEM images in Core KY07-04 PC01 from the northern Okinawa Trough since the LGM at centennial timescales. Available high-resolution geochemical proxy records such as oxygen isotopes and Mg/Ca SST in Core KY07-04 PC01 allow us to directly compare the diatom assemblage and the geochemical dataset (Kubota et al. 2010(Kubota et al. , 2015b. The purpose of this study is to clarify changes in surface water masses such as the Kuroshio Current and the Chinese Coastal Water from the semiclosed LGM to the present via deglacial sea-level rise.

Core KY07-04 PC01 and its age model
A piston core sample KY07-04 PC01 was obtained from a water depth of 758 m on the northern slope of the Okinawa Trough, East China Sea (31° 38.35ʹ N, 128° 56.64ʹ E) during the R/V Kaiyo KY07-04 cruise in February 2007 (Fig. 1). Core KY07-04 PC01 consists of homogeneous olive to gray silty clay with two ash layers (Kubota et al. 2010). The widespread tephra identified the upper ash layer between 408.5 and 488.5 cm in core depth as the Kikai-Akahoya (K-Ah; erupted at 7.3 ka ago).
The age model of KY07-04 PC01 was established based on 18 radiocarbon dating of planktic foraminiferal shells and tephrochronology by Kubota et al. (2010). This study revised the Core KY07-04 PC01 age model with 26 radiocarbon dating of planktic foraminiferal shells (Table 1). New radiocarbon dates for eight samples during the middle to late Holocene interval were measured by the National Electrostatics Corporation 250 kV single-stage accelerator mass spectrometry at Atmosphere and Ocean Research Institute, The University of Tokyo (Yokoyama et al. 2019). We conducted radiocarbon calibration to convert radiocarbon age to calendar age using CALIB 8.2 (Stuiver et al. 2021) with the Marine20 calibration curve (Heaton et al. 2020). Local marine reservoir age (ΔR) was estimated as − 200 ± 100 yr from 33 data from the East China Sea coast (Konishi et al. 1982;Kong and Lee 2005;Hideshima et al. 2001;Southon et al. 2002;Yoneda et al.  2007; Hirabayashi et al. 2017)) in 14CHRONO Marine20 reservoir correction database. The K-Ah tephra age was 7252 ± 56 yr, estimated from Lake Suigetsu varve sediment (McLean et al. 2018). Based on the radiocarbon and K-Ah tephra dates, the revised age model of Core KY07-04 PC01 was established (Table 1; Fig. 2). The revised age model is consistent with the previous one by Kubota et al. (2010). Its continuous sedimentation with a high sedimentation rate of ~ 80 cm kyr −1 ensures Core KY07-04 PC01 for centennial to millennial scales paleoceanographic study.

Microscopic observation for diatoms
Core KY07-04 PC01 samples were sliced at 2.5 cm intervals and stored in a refrigerator. A total of 72 subsamples taken at ~ 20 cm intervals were taken for diatom assemblage work and freeze-dried with an EYELA FDU-1200 freeze-dryer. Approximately 1 mg of the freeze-dried sample was treated with 25 mL of 30% hydrogen peroxide at hydrogen peroxide on a hot plate with top plate temperature set to 100 °C for about 1 h to digest organic matter. After the reaction was complete, Millipore-Elix water was added to a total volume of 200 mL and allowed to stand for 8 h. The removal of supernatant was removed, and Millipore-Elix water was added again. The following is a procedure of slides for LM observation.
The pretreated suspended solution of each sediment sample was drawn by 2 mL plastic pipette, dropped on cover glass (24 × 40 mm), and dried on a hot plate with top plate temperature set to 40-50 °C. Pleurax (Mount Media, Wako Pure Chemical Industries, Tokyo, Japan) was dropped on the cover glass and mounted on the glass slide at 120 °C. LM observations were performed at magnifications of × 1000 using Olympus BX53 and BX50 bright-field microscopes to identify and count diatoms in Core KY07-04 PC01 samples. Counting of diatom valves was performed as follows: (1) for centric diatoms, a valve with a central part was counted as one valve; and (2) for pennate diatoms, a valve with an apical part was counted as 0.5 valves. More than 200 diatom valves were counted for diatom assemblage analyses of Core KY07-04 PC01 samples (Additional file 1: Tables S1 and S2). Sample preparation for SEM observation is described as follows. The pretreated suspended solution of each sediment sample was filtered through an Advantec membrane filter (0.45 µm pore size, 47 mm diameter), rinsed with Millipore-Elix water, and dried. The filter samples were cut into 6 mm × 6 mm squares, mounted onto brass stubs with carbon tape, and coated with OsO 4 using an osmium coater (Neoc-STB; Meiwafosis Co. Ltd., Tokyo, Japan). A field-emission scanning electron microscope (FE-SEM, JMS-7001F; JEOL Ltd., Akishima, Japan) was used for detailed taxonomic work of the diatoms. The acceleration voltage was set at 5 kV.

Results
One hundred twenty-five diatom taxa were observed in Core KY07-04 PC01 samples (Table 2). Paralia, Thalassionema, and Nitzschia were the most abundant diatom genus in Core KY07-04 PC01 samples (Fig. 3)   (periphytic) was the dominant genus before 14 ka ago, whereas Thalassionema and Nitzschia (planktonic) were abundant after 14 ka ago. At the species level, downcore profiles of diatom taxa with greater than 4% at least in a sample were shown (Fig. 3). Paralia spp., periphytic coastal diatom taxa, was representative in Core KY07-04 PC01 samples before 14 ka ago, exceeding 30% during the LGM and early deglaciation from 19 to 15 ka ago (Fig. 3). Another periphytic coastal diatom species, Delphineis surirella, showed a similar pattern with Paralia spp., with high relative abundances before 16 ka ago (Fig. 3). Relative abundances of Nitzschia bicapitata, a major planktonic warm species, were high from the deglacial period to middle Holocene and the core top of Core KY07-04 PC01 (Fig. 3). Thalassionema species and varieties presented various downcore profile patterns in Core KY07-04 PC01 (Fig. 3). Morphological varieties of Thalassionema nitzschioides were persistently high percentages throughout the core. Among the T. nitzschioides varieties, three varieties (var. incurvatum, var. inflatum, and var. parvum) increased during the middle-late Holocene. On the other hand, Thalassionema frauenfeldii sensu stricto (s.s.) and Thalassionema synedriforme were relatively high during the deglacial period. Figure 4 shows diatom groups by habitats (freshwater and brackish) and geographic distributions (coastal, cold, warm) in Core KY07-04 PC01. During the LGM, periphytic diatoms, mostly Paralia spp., were predominant, whereas planktonic diatoms were the minority. During the deglacial pronounced turnover between 18 and 14 ka ago, Planktonic diatoms became the majority in Core KY07-04 PC01. With the turnover, coastal diatoms decreased significantly. Note that contributions of freshwater and brackish water diatoms have been small in total assemblage even during the LGM and early deglaciation (Fig. 4). Marine diatoms were the majority in the diatom assemblages in Core KY07-04 PC01 since the LGM. Planktonic warm water diatom group gradually increased from 18 to 12 ka ago, occupying ~ 40% for the last 12 ka.

Diatom assemblages as paleoenvironmental indicators
Throughout the core, marine diatoms were dominant, suggesting that the northern Okinawa Trough has been a marine environment and limited fresh/brackish water influences even during the LGM. This is consistent with the previous study reporting the continuous occurrences of planktic foraminifera (Xu and Oda 1999;Ijiri et al. 2005) and radiolarians (Matsuzaki et al. 2019) (Asaoka 1980;Oh and Koh 1995;Furuya et al. 1996;Tanimura et al. 2002;Chang et al. 2009b;Ishikawa and Furuya 2004;Guo et al. 2014;Yuan et al. 2014;Chen et al. 2019b). Paralia spp. distributes on muddy sediments of the Japanese coast with salinity ranging from 11 to 34 (Sawai et al. 2005). In the modern East China Sea and the Yellow Sea, Paralia spp. dwells on the continental shelves as a representative coastal water indicator throughout the year (Asaoka 1980;Furuya et al. 1996;Ishikawa and Furuya 2004;Guo et al. 2014;Yuan et al. 2014;Chen et al. 2019b), including the eutrophic coastal area of the Yellow Sea Liu et al. 2013). Asaoka (1980) revealed that Paralia spp. distributions in summer were almost synchronous with suspended mud deposited on the central part of East China Sea continental shelves. Furuya et al. (1996) and Ishikawa and Furuya (2004) reported dominant Paralia spp. in the inner shelf area (< 70 m water depth) by active vertical mixing in winter and spring, respectively. Guo et al. (2014) presented the seasonality of diatom assemblages and the abundances in the continental shelf area   Table 2). Yellow triangles and a red diamond on the right denote radiocarbon dates and K-Ah tephra, respectively Page 10 of 17 Shirota et al. Progress in Earth and Planetary Science (2021) 8:66 of the East China Sea. Based on their study, Paralia spp. was one of the most abundant diatom taxa in spring, fall, and winter. On the other hand, the occurrence of Paralia spp. was low in summer. These modern observations have clarified that abundant Paralia spp. occurrences in the nutrient-rich shelf water of the East China Sea with temperatures ranging from 11 to 20 °C and salinity from 33.5 to 34.5. This water-mass property corresponds to the water mass properties of Chinese coastal water (Ichikawa and Beardsley 2002). Paralia spp. is transported with suspended mud from the coastal area/continental shelves to the northeastern East China Sea. Tanimura et al. (2002) investigated Paralia spp. abundances in KH82-4 St-14 piston core (31° 44.4ʹ N, 129° 02.1ʹ E, 740 m water depth) adjacent to KY07-04 PC01. They found abundant Paralia spp. occurrences from 19 to 13 cal ka BP (16-10 ka BP in 14 C age shown in Tanimura et al. (2002), suggesting a significant expansion of Chinese coastal water sea-level drop during the LGM to early deglaciation. A similar but much smaller (up to 20%) abundance peak during the LGM to early deglaciation was observed at Core MD01-2404 in the southwestern Okinawa Trough (Chang et al. 2009b). Delphineis surirella, a periphytic coastal diatom species that dwell in the subtidal zone (Zong and Horton 1999;Sawai et al. 2009), showed relatively high abundances during 19-14 ka ago, presenting a similar pattern with Paralia spp. (Fig. 3). These consistently support the expansion of Chinese coastal water during the LGM to early deglaciation.

. Gray hatches indicate MWP-A and YD periods
Page 11 of 17 Shirota et al. Progress in Earth and Planetary Science (2021) 8:66 the Kuroshio (the shortest distance from the trap station to the Kuroshio path). When Kuroshio migrated to the north, the sediment trap station was covered with subtropical water where the three T. nitzschioides varieties were abundant. In Core KY-07-04 PC01, T-ratio has increased since 12 ka ago and has reached its present level since 8 ka ago (Fig. 5), which is consistent with the timings of northward migration of Kuroshio in the western North Pacific (Takemoto and Oda 1997;Tanimura 1999;Yamamoto et al. 2005) and intensification of Tsushima Warm Current flowing into the Japan Sea (Gallagher et al. 2015;Yokoyama et al. 2006), respectively. A typical dark color layer in many deep-sea sediment cores from the Japan Sea with high organic carbon contents, called TL1, is seen when the Tsushima Warm Current intrusion was re-established since the LGM (Oba et al. 1991;Tada et al. 1999;Yokoyama et al. 2006). Thus, the T-ratio increase suggests strengthening the influence of Kuroshio Current water in the northern East China Sea. Tanimura et al. (2007) recognized nine species, varieties, and morphological types of Thalassionema nitzschioides. Akiba et al. (2014) described the Nd morphological type of Thalassionema nitzschioides sensu lato (s.l.) as a new species Thalassionema umitakae. The distribution of T. umitakae in the North Pacific is mainly in the subarctic region but a few occurrences in the subarctic boundary of the western North Pacific (Tanimura et al. 2007). Abundant occurrences of T. umitakae were reported in the glacial Japan Sea (Akiba et al. 2014). In Core KY07-04 PC01, T. umitakae occurred throughout the core since the LGM with no pronounced trend (Fig. 3), implying no significantly cold surface water covered in the northern Okinawa Trough even during the LGM.  (Sun et al. 2005) and Core DGKS9604 (Yu et al. 2009); and the northern Okinawa Trough: Core KY07-04 PC01 (Kubota et al. 2010(Kubota et al. , 2015b and Core MD98-2195 (Ijiri et al. 2005) as well as residual δ 18 O as a salinity proxy in Core KY07-04 PC01 (Kubota et al. 2010(Kubota et al. , 2015a. SST and residual δ 18 O at Core KY07-04 PC01 and the errors (± 1σ) were derived by Paleo-Seawater Uncertainty Solver (PSU Solver) (Thirumalai et al. 2016) with the methodology of Kubota et al. (2019). Downcore profiles of relative abundances (%) of warm diatom indicators and Thalassionema variety index in Core KY07-04 PC01 were also shown. Gray hatches indicate MWP-A and YD periods. Horizontal broken lines denote pronounced decreases in warm diatom indicators during Holocene Page 12 of 17 Shirota et al. Progress in Earth and Planetary Science (2021) 8:66 Thalassionema frauenfeldii was relatively abundant in the western margin of the North Pacific (Tanimura et al. 2007). In this study, we counted T. frauenfeldii as separated the two morphological types T. frauenfeldii s.l. and T. frauenfeldii s.s. Simonsen (1974) described Thalassiothrix frauenfeldii (= Thalassionema frauenfeldii s.l.), T. frauenfeldii s.l. is reported from the open ocean of the western North Pacific (Hasle 2001;Tanimura et al. 2007). On the other hand, Thalassionema frauenfeldii s.s. is likely to be a coastal species because Hallegraeff (1986) and Sar et al. (2007) collected this from the Argentine and Korean coasts. In Core KY07-04 PC01, T. frauenfeldii s.l. gradually increased since the LGM, whereas T. frauenfeldii s.s. was relatively high abundances during the LGM to the deglacial period with a gentle peak centered at ~ 12 ka ago (Fig. 3). The overall trend suggests a transition from semi-closed to open ocean environments during the deglacial period.
Warm water taxa dwelling in the open ocean gradually increased since the LGM and reached ~ 30% at 8 ka ago (Fig. 3). Among them, Nitzschia bicapitata is the major open-ocean species abundant in tropical to temperate oceans (Fryxell 2000), including in the western Pacific warm pool in the western equatorial Pacific with > 28 °C sea-surface temperature (Kobayashi and Takahashi 2002) and Kuroshio region (Tanimura 1992). In Core KY07-04 PC01, relative abundances of N. bicapitata were high between 12 and 5 ka ago, with a rapid increase near the core top (Fig. 3). Li et al. (2012Li et al. ( , 2015 also reported a high percentage of N. bicapitata for the last 400 ka in Core MD05-2908 from the southern Okinawa Trough off northeast of Taiwan. Because of the restricted N. bicapitata distribution in the open ocean, abundant N. bicapitata suggests strengthening Kuroshio inflow into the East China Sea. During the early to middle Holocene, Kuroshio intensified due to strengthening the Subtropical Gyre by precessional forcing (Isono et al. 2009;Yamamoto 2009).

Changes in surface water mass in the East China Sea since the LGM
Diatom assemblages dominated by Paralia spp. in Core KY07-04 PC01 indicate intensified Chinese coastal water whereas weak Kuroshio water in the northern East China Sea during the LGM. Planktic foraminiferal assemblage records also suggested enhanced coastal water during the LGM (Xu and Oda 1999;Ijiri et al. 2005). Reconstructed SSTs during the LGM were ranging from 19 to 22 °C in the northern East China Sea based on planktic foraminiferal Mg/Ca (Kubota et al. 2010;Ujiié et al. 2016), alkenone (Ijiri et al. 2005), and radiolarian transfer function (Matsuzaki et al. 2019) (Fig. 5). In the central and southern Okinawa Trough, reconstructed SSTs during the LGM were between 21 and 24 °C (Sun et al. 2005;Chang et al. 2009a;Yu et al. 2009;Chen et al. 2010;Xu et al. 2013) (Fig. 5). These SSTs were 2-5 °C lower than the present ones but still within temperate to subtropical ranges, consistent with continuous occurrences of warm open-water diatom taxa in Core KY07-04 PC01 since the LGM. In addition, the regional high-resolution climate model study by (Vogt-Vincent and Mitarai 2020) supports scenario-based proxy records. Kuroshio kept flowing into the Okinawa Trough during the LGM but significantly weakened. However, Kuroshio's path and intensity during the LGM varied depending on the models. During the LGM, the Kuroshio path migrated to the south of the Ryukyu Islands suggested by a model (Kao et al. 2006). On the contrary, a simulation indicates stronger Kuroshio flowing into the Okinawa Trough (Lee et al. 2013).
Expansion of the nutrient-rich Chinese coastal water led to enhanced productivity in the glacial East China Sea (Ijiri et al. 2005;Chang et al. 2009a;Kubota et al. 2017;Matsuzaki et al. 2019). Dominant Paralia spp. in Core KY07-04 PC01 as an indicator of nutrient-rich coastal water (McQuoid and Nordberg 2003;Di et al. 2013;Liu et al. 2013) demonstrated high productivity in the glacial northern East China Sea. The Chinese coastal water prevailing in the northern East China Sea was an advancement of the river mouth of the Yellow River with sea-level drop during the LGM. Clay mineral compositions and Sr-Nd isotopes of the detrital matter revealed an ample detrital supply from the old Yellow River to the northern East China Sea during the LGM (Xu et al. 2014;Zhao et al. 2017Zhao et al. , 2018Beny et al. 2018). Tsushima Warm Current (TSWC) gradually intensified since 16 ka ago, elevated at 11 ka, and fully evolved at ~ 7 ka after the K-Ah eruption in the northern Okinawa Trough (Xu et al. 2014;Zhao et al. 2017). These timings were consistent with the increasing patterns of %warm diatoms and T-ratio (Fig. 5). In the central Okinawa Trough, the Yangtze River played a critical role in detrital supply (Diekmann et al. 2008;Dou et al. 2010aDou et al. , b, 2012. Relative abundances of Paralia spp., periphytic diatom taxa, in Core KY07-04 PC01 gradually decreased since the LGM, subsequently rapidly decreased from 15 to 14 ka ago (Fig. 3). One of the most pronounced deglacial events is meltwater pulse 1A (MWP-1A) between 14.6 and 14.3 ka ago when relative sea-level rose at Tahiti, Barbados, and Sunda Shelf about 10-20 m from − 110 to − 90 m (Deschamps et al. 2012). Due to the East China Sea location as categorized as a far-field (i.e., sites far away from past and present ice-covered regions (Yokoyama and Purcell 2021), a similar magnitude of relative sea-level rise is expected. With the rapid sea-level rise during MWP-1A, river mouths of the Yellow River and the Yangtze River must have retreated.
Page 13 of 17 Shirota et al. Progress in Earth and Planetary Science (2021) 8:66 The supporting evidence is by records of clay mineral compositions and Sr-Nd isotopes of detrital matter, suggesting switch of the source region from the river to the continental shelves (Dou et al. 2010a(Dou et al. , b, 2012Xu et al. 2014;Zhao et al. 2017;Hu et al. 2020). Gradual increases in warm diatoms, T-ratio, and Mg/Ca SST in Core KY07-04 PC01 indicate a gradual intensification of Kuroshio during the deglacial period. During the Younger Dryas between 12.9 and 11.6 ka ago (Cheng et al. 2020), warm diatoms kept increasing, whereas Mg/Ca SST was relatively low (Fig. 5).
In the northern Okinawa Trough, Kuroshio intensified between 12 and 5 ka ago, reaching its maximum at ~ 8 ka ago based on diatom taxa living in warm, open water (Fig. 5). Multi-proxy records also indicate the intensified Kuroshio during late deglacial to middle Holocene, i.e., SST rise (Ijiri et al. 2005;Kubota et al. 2010Kubota et al. , 2015bRuan et al. 2015;Zhao et al. 2015), deepening of the thermocline , planktic foraminiferal assemblage (Xu and Oda 1999;Ijiri et al. 2005), floral changes of calcareous nannoplankton assemblage (Su and Wei 2005), sortable silt (Diekmann et al. 2008), terrigenous input from Taiwan (Dou et al. 2010a(Dou et al. , b, 2012Xu et al. 2014Xu et al. , 2019Zheng et al. 2016;Zhao et al. 2017;Hu et al. 2020), organic carbon and reactive phosphorus (Shao et al. 2016;Li et al. 2018), mercury concentration (Lim et al. 2017), and redox sensitive metals (Zou et al. 2020). Warm diatoms showed several rapid decreases based on the direct comparison between diatom assemblage and Mg/Ca SST in Core KY07-04 PC01. Some of these events accompany Mg/Ca SST drops during the Holocene (Fig. 5). An event at 7 ka ago was the most prominent. Subsequently, Kuroshio intensity might have weakened during the late Holocene due to its southward migration by precessional forcing controlling the strength of the North Pacific subtropical gyre (Isono et al. 2009;Yamamoto 2009). However, a sudden increase in N. bicapitata, representative diatom species dwelling in the open ocean with high SST, was observed near the core top of KY07-04 PC01 (Fig. 3). High-resolution diatom assemblage record during the late Holocene in the southern Okinawa Trough showed high abundances of N. bicapitata for the last ~ 400 yrs . Such high abundances of N. bicapitata over the Okinawa Trough is consistent with the finding of intensification of Kuroshio during the latest Holocene (Zhang et al. 2019), implying persistent poleward heat transport regardless of the weakening of the Atlantic Meridional Overturning Circulation (Summerhayes and Zalasiewicz 2018). Alternatively, the high N. bicapitata abundances might be related to the response to the global warming (Chen et al. 2019a) and southward course shift of the Yellow River discharged into the Yellow Sea between 1546 and 1855 A.D. (Saito and Yang 1995).

Conclusions
We conducted microscopic observations for diatom assemblages since the LGM in Core KY07-04 PC01 from the northern Okinawa Trough, East China Sea. The produced diatom assemblage data are the first centennialscale resolution since the LGM in the East China Sea with SEM and LM photographs of the representative taxa. Most of the diatoms found in Core KY07-04 PC01 were marine taxa throughout the core, with a few freshwater/ brackish water taxa. There was a substantial turnover of major diatom groups from periphytic to planktonic taxa during deglaciation between 17 and 14 ka ago. During the LGM and early deglaciation, the representative periphytic taxa were Paralia spp., abundant in the modern continental shelf region of the East China Sea covered by the Chinese coastal water. The advancement of the Yellow River mouth with a sea-level drop during the LGM was the plausible reason for the prevailing Chinese coastal water in the northern East China Sea. Rapid decrease in the relative abundances of Paralia spp. during 15-14 ka ago corresponded to MWP-1A, a major deglacial event with the rapid sea-level rise during 14.6-14.3 ka ago. Planktonic cold diatom taxa were less than 10% throughout the core, suggesting the oceanic environment in the northern Okinawa Trough has been temperate to subtropical even during the LGM. Planktonic diatom Thalassionema was the dominant diatom genus since 14 ka ago. A ratio of three T. nitzschioides varieties (var. incurvatum, var. inflatum, and var. parvum) in Thalassionema nitzschioides complex (T-ratio) increased since 12 ka ago and has reached its present level since 8 ka ago, responding to strengthening and northward migration of Kuroshio during early to middle Holocene by precessional forcing. Warm planktonic diatom taxa also showed a gradual increase from 12 to 8 ka ago. A sudden increase in N. bicapitata, warm open ocean species, was found near the core top, which might be related to the intensifying Kuroshio during the latest Holocene.