Characterization of water masses around the southern Ryukyu Islands based on isotopic compositions

We investigated the water-mass structure on the Okinawa Trough and Pacific sides of the southern Ryukyu Island Arc (Yonaguni, Iriomote, and Ishigaki subareas) using the Nd isotope composition (143Nd/144Nd ratios; expressed as εNd values) of benthic foraminiferal tests in surface sediments, which reflect bottom-water composition, along with hydrogen and oxygen isotope compositions (δD and δ18O values, respectively) and physical properties (temperature and salinity) of seawater. The Okinawa Trough side has lower εNd values than the Pacific side due to continental/island material inputs characterized by relatively low εNd values. Moreover, within the Okinawa Trough, other processes control the Nd behavior of seawater and primarily affect the Yonaguni and Iriomote subareas, as follows. (1) Surface and subsurface waters are influenced by Taiwanese river discharge combined with temporospatial variations in oceanographic conditions including Kuroshio Current meandering. (2) Intermediate water is characterized by low εNd values (down to − 8.2), possibly attributable to sediment plumes and turbiditic fluxes. (3) The εNd values of bottom water indicate upwelling and vertical mixing, with composition therefore being similar to those of intermediate water. The εNd profiles are better defined on the Pacific side. High εNd values occur in surface and subsurface (< 300 m depth, potential density < 25.0 kg m−3) waters, and low values (down to − 7.0) occur in subsurface–core-intermediate water (400–600 m depth, 26–27 kg m−3). εNd values increase slightly to − 4.0 below 750 m depth and remain constant down to about 2000 m depth, below which deep water shows a slight decrease in εNd values. Intermediate and bottom/deep waters are distinguished from upper layers by their lower δD and δ18O values.


Introduction
The Kuroshio Current (KC) enters the East China Sea (ECS) through the strait between Taiwan and Yonaguni Jima Island (TY), flows northeastward in the Okinawa Trough (OT) along the Ryukyu Island Arc (RIA), and bifurcates to the southwest of Yaku Shima Island (Fig. 1). The main current changes direction, veering towards the Pacific in the Tokara Strait (TS) area, where the Tokara Gap (TG) lies between Amami Oshima and Yaku Shima islands. The KC is a warm current that transports heat from the tropics to the Ryukyu Islands (the Ryukyus) and the southern half of mainland Japan.
It has enabled the development of reefs with highly diversified hermatypic coral fauna throughout the Ryukyus, despite their relatively high latitudes within the reef province (Veron 1992) since the opening of the OT and subsequent influx of the KC into the backarc basin (Iryu et al. 2006). Knowledge of the pathway, strength, and physicochemical properties of the KC is critical for understanding modern oceanographic and climatic conditions in East Asia, and a detailed study of water-mass characteristics around the RIA is expected to contribute to this knowledge.
Planktic foraminiferal tests record the εNd values of ambient seawater at various depths during their lifespan. After planktic foraminifers die and their tests are transported to the sediment-water interface, the original εNd values are modified to reach isotopic equilibrium with bottom waters over at least 300 yrs (Palmer 1985;Roberts et al. 2012). However, benthic foraminiferal tests comprise calcite precipitated in Nd isotope equilibrium with ambient seawater during both their life and death stages. Because benthic foraminifers also have wide temporospatial distributions, the εNd values of benthic foraminiferal tests have been used to detect variations in the εNd values of bottom water due to water-mass structures (Murray 2013). There is no significant difference in εNd values between a sample of a single benthic foraminiferal species and a sample of multiple foraminiferal species (Haley et al. 2005;Roberts et al. 2012;Wu et al. 2015a); therefore, when single-species samples are limited or unavailable, the εNd values of multi-species samples can be used as a reliable proxy for identifying a water mass and its temporospatial variation (e.g., Palmer 1985;Lynch-Stieglitz and Marchitto 2014).
The ECS is influenced by multiple sources of Nd. The primary source is the Philippine Sea (PS), where εNd values range from − 4.4 in deep water to − 2.0 in surface-subsurface water (Wu 2014;Behrens et al. 2018). Although in lower proportions, some influence of the South China Sea (SCS) is known to reach the ECS through the TY and the Kerama Gap (KG; Chen 2005;Nakamura et al. 2013). The εNd values in the northeastern SCS vary from − 4.5 in deep water to − 3.6 in subsurface water, with a surface-water εNd value of − 4.5. Many rivers in China flow into the ECS, with the most significant influence being from the Changjiang (εNd = − 11.6), Yellow (− 13.9 to − 9.6), and Yangtze (− 14.3 to − 8.7) rivers (Meng et al. 2008;Dou et al. 2012). The OT is also influenced by rivers in Taiwan (εNd = − 11.1; Lan et al. 2002).
There have been few Nd isotope studies focusing on the OT, and even fewer on the KC. Amakawa et al. (2000) reported εNd values of Kuroshio Surface Water (KSW) at < 5 m depth southeast of Taiwan (εNd = − 4.8 at station PA-S-18; 21.83°N, 121.11°E; Fig. 1) and in the western-central OT (− 4.9 at station PA-S-19; 26.30°N, 125.54°E; Fig. 1), with these values being lower than those measured to the east of Luzon Strait in the PS (− 4.2 at station PA-S-1; 19.54°N, 131.27°E; Fig. 1). A Nd isotope profile at station LM-6/11 (34.16°N, 142°E; Fig. 1) exhibits lower values than those measured to the east of Taiwan (Amakawa et al. 2004b), with εNd values of − 5.6 at 5 m depth, − 8.7 at 251 m, − 6.0 at 503 m, and − 5.9 at 1989 m. This SW-NE decrease in εNd values in the KC has been attributed to riverine sediment inputs from the ECS (Amakawa et al. 2004b;Behrens et al. 2018) and transportation of less radiogenic continental Nd (relatively low εNd values) into the central Pacific by the KC (Amakawa et al. 2004a). Dou et al. (2012) identified a post-glacial increase in εNd values during the last deglaciation (from − 11.4 at 14 ka to − 10.7 at 7 ka) in sediments collected at 766 m water depth on the slope of the ECS shelf, and they interpreted this as reflecting a reduced influence of riverine input from China due to strengthening of the KC. However, no studies have determined the Nd isotope composition in the KC around the RIA at high spatial resolution and over a wide range of depths, with most previous studies focusing on the KC before entry or after exit of the OT.
Stable hydrogen and oxygen isotope compositions, D/H and 18 O/ 16 O expressed as δD and δ 18 O values, respectively (Sect. 3.2), provide crucial oceanographic information regarding thermohaline circulation and the formation, mixing, and ventilation of water masses. δD and δ 18 O values vary spatially in the oceans depending on advection processes and the mixing of water masses with different isotopic compositions (Xu et al. 2012;Rohling 2013). Transition phases of the hydrological cycle, turbulence fluxes, and Rayleigh distillation affect the δD and δ 18 O values of surface water, which change progressively from low to high latitudes. For example, heavier isotopic compositions occur in areas of high evaporation near the Equator, whereas lighter compositions are observed near the poles (Oka and Kawabe 1998;Schmidt et al. 1999;Rohling 2013). The relationships of temperature and salinity to δD and δ 18 O values have been broadly used to identify water masses and ocean dynamics (e.g., Horibe and Ogura 1968;Schmidt et al. 1999;Bigg and Rohling, 2000;LeGrande and Schmidt 2006;Sengupta et al. 2013;Wu et al. 2021). Global datasets indicate lower δ 18 O values near coastal areas such as the ECS and the PS, which have δ 18 O values of -0.1‰ and + 0.1‰, respectively (Bigg and Rohling 2000;LeGrande and Schmidt 2006).
Evaporation and precipitation control salinity-δD and salinity-δ 18 O relationships of shallow water distal to continental inputs. When evaporation exceeds precipitation and freshwater discharge, δ 18 O and δD values and salinity increase. Conversely, when precipitation and freshwater discharge exceed evaporation, δ 18 O and δD values and salinity decrease (Rohling 2013;Sengupta et al. 2013). Furthermore, physical oceanographic processes such as sinking, advection, and upwelling may cause different water masses to have distinct slopes and intercepts China C o n t i n e n t a l i n p u t T a i w a n SCS R y u k y u  (Amakawa et al. 2004b); PA-S-1, PA-S-18, and PA-S-19 (Amakawa et al. 2000); 2, 3, 4, and 5 (Behrens et al. 2018); LS C3a, P2, SCS C5, and WPS N18-05 (Wu 2014 Horibe and Ogura (1968) (Additional file 5: Fig. S3). The horizontal yellow line in the Luzon Strait indicates a transect studied by Mensah et al. (2014) Page 4 of 20 Cruz Salmeron et al. Progress in Earth and Planetary Science (2022) 9:44 on salinity-δD and salinity-δ 18 O cross-plots (Sengupta et al. 2013). The KC is a source of warm oligotrophic water, so its influx into the OT has been important for the evolution of oceanographic conditions and climate states in Far East Asia. However, past flow paths of the KC and the timing of KC entrance to the OT are poorly constrained. Existing proxies, such as δ 13 C and δ 18 O values of microfossils, numerical models, and sediment provenance show some discrepancies (e.g., Ujiié et al. 2003;Kao et al. 2006;Dou et al. 2010Dou et al. , 2012Lee et al. 2013). This study aimed to establish a new proxy of Nd isotope composition that may be used to reconstruct the past water-mass structures of the OT and around the RIA. From a physical oceanography perspective, the OT is a challenging study area due to its complex bathymetry, seasonal-tointerannual variability in seawater physical properties, occurrence of several water masses, and riverine inputs. Therefore, Nd isotopes are expected to provide vital information on the current system and water-mass structure within the OT.
In this study, εNd values were determined by analysis of a small number of modern benthic foraminifers, and the values are directly compared with εNd data from previous studies. This is the first study to delineate variations in εNd values with wide spatial coverage in the KC and related water masses on the OT and Pacific sides of the southern RIA at high spatial resolution, based on εNd, δD, and δ 18 O values and physical properties of seawater (e.g., temperature and salinity). The study area is located near the KC entrance to the OT, so it is an ideal site for characterization of the Nd isotope signature of the KC immediately after it interacts with terrigenous materials from the Eurasian continent and Taiwan.

Geological settings
The Ryukyus include more than 200 islands and islets to the southwest of mainland Japan, extending from Tanegashima Island (30.73°N, 131°E) in the northeast to Yonaguni Jima Island (24.45°N, 123°E) in the southwest (Fig. 1). The islands are arranged in an arc known as the RIA, an active island arc related to subduction of the Philippine Sea Plate along the Ryukyu Trench beneath the Eurasian Plate. The RIA is bounded by the ECS to the northwest and the Pacific Ocean to the southeast and is geographically divided into three regions (North, Central, and South Ryukyus;Iryu et al. 2006;Arai et al. 2014) by the TG and KG (Fig. 1). The OT is a back-arc basin separating the RIA from the ECS shelf. The study area is located in the western part of the South Ryukyus. The study area on the OT side is divided into three subareas ( Fig. 2): Yonaguni (Yonaguni OT), Iriomote (Iriomote OT), and Ishigaki (Ishigaki OT). On the Pacific side, it is also divided into three subareas: Yonaguni (Yonaguni Pacific), Iriomote (Iriomote Pacific), and Ishigaki (Ishigaki Pacific).
The OT has a NE-SW length of ~ 1000 km and a NW-SE width of 200 km. The bathymetry of the trough ranges from shallower depths of ~ 600 m in the north to maximum depths of ~ 2100 m in the Yaeyama and Miyako rifts in the south. Three channels connect the OT with the PS; i.e., the TY, KG, and TG, with maximum water depths of 1700, 2100, and 1500 m and sill depths of 775, ~ 1100, and 700 m, respectively (Arai et al. 2013;Nakamura et al. 2013;Nishizawa et al. 2019).
The opening of the OT and subsequent influx of the KC into the back-arc basin was a key event in the Cenozoic environmental and oceanographic evolution of Far East Asia including the Ryukyus, with the KC being a source of warm oligotrophic water. The initial formation of coral reefs in the Ryukyus is considered to record the influx of the KC into the ECS (i.e., the back-arc side of the RIA). This changed the oceanographic conditions from a "mud sea, " represented by siliciclastic deposits of the upper Miocene-lower Pleistocene Shimajiri Group, to a "coral sea, " represented by carbonates of the Pleistocene Ryukyu Group that formed in reefs and associated fore-reefs and island shelves (Iryu et al. 2006;Gallagher et al. 2015;Imai et al. 2017).

Oceanographic settings
Here we introduce general aspects of the ocean circulation and water masses in the study area, referring to oceanographic data for the three subareas on the OT and Pacific sides of the RIA (Additional file 1: Tables S1, Additional file 2: Table S2; Fig. 2, Additional file 3: Fig.  S1) provided by the Japan Oceanographic Data Center (Japan Oceanographic Data Center, 2021). The JODC temperature and salinity data used in this study were collected from 1 January 1990 to 31 December 2019. Long-term (30 yr) averaging of the data has eliminated short-term variations in temperature and salinity, such as those caused by seasonality, the El Niño-Southern Oscillation, and variability in the KC flow path (Nitani 1975). Therefore, the JODC data are more suitable for characterization of the water masses investigated here than onboard data obtained during the GK19 cruise (Additional file 1: Tables S1, Additional file 2: Table S2). We identified water masses following water-mass classifications of previous studies (Nakamura 1996;Wang and Chen 1998;Qu et al. 1999;Ichikawa and Chaen 2000;Amakawa et al. 2004b;Chen 2005;Talley 2008;Kawabe and Fujio 2010;Talley et al. 2011;Amakawa et al. 2013;Qi et al. 2014;Wu 2014;Yang et al. 2015;Behrens et al. 2018;Zuo et al. 2019). Using the JODC salinity and Page 5 of 20 Cruz Salmeron et al. Progress in Earth and Planetary Science (2022) 9:44 temperature data, we calculated potential temperatures and densities and then correlated our results using previous data for the OT and PS (Tables 1 and 2).

Okinawa trough side of the Ryukyu Island Arc
The North Equatorial Current (NEC) bifurcates into the northward-flowing KC and southward-flowing Mindanao Current to the east of the Philippine coast. The KC enters the OT, flows northeastward along the RIA, and drifts out to the Pacific through the TS. The core of the KC has an average thickness of 500-700 m but reaches depths of up to 1350 m in some areas (Chen 2005;Kubota et al. 2017). Based on physical and chemical properties (e.g., temperature, salinity, nutrients, and pH), three water masses have been identified within the OT in the vertical profile of the KC path: KSW (< 150 m depth), Kuroshio Tropical Water (KTW; 150-300 m depth), and Kuroshio Intermediate Water (KIW; 500-700 m depth; Chen et al. 1995;Wang and Chen 1998;Chen 2005;Mensah et al. 2014). Distinctive KC temperature and salinity fronts exist between the continental shelf region and its offshore region in the OT, with salinity being lower in the former due to freshwater discharge from rivers in China (Oka and Kawabe 1998).
The JODC temperature and salinity data indicate the presence of KSW above ~ 100 m depth, with salinity of 34.47-34.78 and potential density (σ θ ) of < 23.0 kg m −3 . Page 6 of 20 Cruz Salmeron et al. Progress in Earth and Planetary Science (2022) Na et al. 2014). As a result, the salinity of KIW falls to levels between those of NPIW and SCSIW. KIW exists at 500-755 m depth in the OT side, characterized by its low salinity of 34.29-34.36 and σ θ of 26.5-27.0 kg m −3 . The intermediate water is more saline by ~ 0.11 on the OT side of the RIA than on the Pacific side, indicating mixing with more saline water within the OT, as for SCSIW (Fig. 3). In the southern OT, the mixing proportions of SCSIW and NPIW are 55% and 45%, respectively (Nakamura et al. 2013). A salinity front separating western high-salinity water from eastern low-salinity water has been detected to the east of Taiwan at 21.45°N near 122.5°E (Chen 2005), suggesting that SCSIW encounters NPIW in this area ( Fig. 1). Bottom waters in the OT have salinity of 34.63 and σ θ of > 27.5 kg m −3 (Table 1; Fig. 4).

Pacific side of the RIA
North Pacific Subtropical Gyre water flows northeastward across the PS on the Pacific side of the RIA, and the subsurface layer has been identified as the Ryukyu Current System (Nakamura et al. 2007; Thoppil et al. 2016).
The JODC data indicate that the salinity and σ θ of surface water (< 83 m depth) are 34.31-34.74 and < 23.00 kg m −3 on the Pacific side of the RIA, respectively. However, these values vary seasonally, with salinity being lower in summer and higher in winter (Rudnick et al. 2011). Moreover, the salinity of surface-subsurface water increases in an E-W direction (Chen 2005;Mensah et al. 2014). Salinity maxima correspond to NPTW at 99-216 m depth in the study area (34.67-34.90; σ θ = 23.5-24.7 kg m −3 ). Central Water, formed in the subtropics and Equatorial Pacific, overlies North Pacific Intermediate Water (NPIW) and is characterized by a wide range of potential temperature (θ, °C) and salinity. Central Water spreads over the thermocline through advection and diffusion processes over the isopyncals after subduction from the surface mixed layer. For North Pacific Central Water, upper and lower limits correspond to North Pacific Subtropical Mode Water (σ θ = 25.5 kg m −3 ) and North Pacific Central Mode Water (σ θ = 26.3 kg m −3 ), respectively (Nakamura, 1996;Talley et al., 2011). Here we refer to this water mass as being WNPCW, identified as having salinity of 34.37-34.82 and σ θ of 24.9-26.0 kg m −3 at 222-476 m depth (Fig. 4).
NPIW has been identified as the minimum-salinity layer in the western PS (Wang and Chen 1998;Qu et al. 1999;Chen 2005;Behrens et al. 2018; Table 2). In the study area, NPIW has salinity of 34.18-34.39 and σ θ of 26.3-27.2 kg m −3 .
In the Pacific Ocean, the southern-sourced Lower Circumpolar Deep Water and Upper Circumpolar Deep Water (UCDW) flow in a S-N direction. UCDW upwells at high latitudes in the North Pacific and flows southward as Pacific Deep Water (PDW; Talley et al. 2011). In the PS, PDW can be traced at > 2000 m depth with salinity of 34.6-34.7, θ of 1.6 °C-2.1 °C, and σ θ of > 27.6 kg m −3 (Amakawa et al. 2004b(Amakawa et al. , 2013Talley 2008;Wu et al. 2015a; Table 2). PDW flows above UCDW to the east of Luzon Strait in the western PS (Wu et al, 2015b;Wan et al. 2018), whereas it is at 1907-2620 m depth in the study area, where salinity is 34.59-34.64 and σ θ > 27.7 kg m −3 .

Sampling and sample processing
A total of 189 bottom-water and sediment samples were collected using a K-grab sampler from the area between Ishigaki Jima and Yonaguni Jima islands over a depth range of 68-2620 m during the GK19 cruise of R/V Kaiyo-maru No. 1 undertaken during June 23-July 22, 2019, by the Geological Survey of Japan of the National Institute of Advanced Industrial Science and Technology (AIST/GSJ; Additional file 1: Tables S1, Additional file 2: Table S2; Fig. 2). Bottom-water samples were collected from ~ 7 m above the seafloor using a Niskin sampler attached to a K-grab sampler. Some physical and chemical properties were measured using a conductivity-temperature-depth (CTD) profiler (CTD 90 M; Sea & Sun Technology GmbH, Germany) attached to the grab sampler. Temperature, salinity, turbidity, pH, and dissolved-oxygen data are presented in Additional file 1: Tables S1 and Additional file 2: Table S2. Surface sediment samples (5-10 cm 3 ; 0-2 cm depth) collected from 42 sites were washed using ultra-pure water and a 63 µm mesh sieve. Particles of > 63 µm size were air-dried at 40 °C for 12 h and sieved using 125 and 250 µm sieves. Samples were divided into three size fractions: 63-125, 125-250, and > 250 µm. Benthic foraminiferal tests (50-70) were handpicked from the > 250 µm fraction under a binocular microscope (Stemi 2000-C; Zeiss, Germany). Although the ages of the samples were not determined, the samples were obtained from 0-2 cm below the seafloor, and previous studies have dated such seafloor sediment at < 1 ka near the study site (Chang et al. 2008;Dou et al. 2012 (2000) Page 9 of 20 Cruz Salmeron et al. Progress in Earth and Planetary Science (2022)

Sediment sample preparation and Nd isotope analyses
Forty-two sediment samples were selected for Nd isotope analysis. The cleaning of bulk benthic foraminifers followed Lea et al. (2000) and Sagawa et al. (2005). Bulk benthic foraminiferal samples were collected in a 1.5 mL micro-centrifuge tube and decontaminated using a multistep cleaning procedure, including clay removal using ultrapure water and methanol, reductive cleaning to remove oxidized metal coatings using hydrazine/ammonium citrate, and oxidative cleaning to remove organic matter using H 2 O 2 buffered by NaOH. The cleaning was undertaken at the Institute of Geology and Paleontology, Graduate School of Science, Tohoku University, Japan. Nd isotope ( 143 Nd/ 144 Nd) analysis of sub-nanogramnanogram Nd samples followed Wakaki and Ishikawa (2018). The Nd in bulk benthic foraminiferal samples was chemically separated and purified using a three-step extraction and ion-exchange chromatographic procedure. First, Nd was separated from major elements and other REEs using TRU and Ln resins (Eichrom, USA) following a procedure modified from that of Pin and Santos (1997). Sample solutions were first loaded onto Teflon columns filled with 0.2 mL of TRU resin. Major elements were eluted with 4 mL 3 M HNO 3 , before light REEs (LREEs) were eluted with 2 mL H 2 O. Eluted LREEs were loaded directly onto a second Teflon column filled with 0.5 mL Ln resin. La, Ce, and Pr were eluted with 3.5 mL 0.2 M HCl, and the Nd fraction with 3 mL 0.2 M HCl. Finally, the separated Nd fraction was loaded onto a third Teflon column filled with 0.1 mL cation-exchange resin (AG50W-X12; Bio-Rad Laboratories, Inc., Japan) for Nd purification. Impurities, remaining cations, and organic materials that may affect Nd ionization during thermal ionization mass spectrometry (TIMS) were eluted with 1.5 mL 2.5 M HCl, and the purified Nd fraction was collected in 2 mL 6 M HCl and evaporated to dryness.
The purified Nd was dissolved in 2.5 µL 3 M HNO 3 and loaded onto degassed Re double filaments with 0.8 µL 0.3 M H 3 PO 4 . Nd isotope ratios were determined using TIMS (Triton, Thermo-Finnigan, USA) at the Kochi Core Center, Japan, using the total evaporation normalization (TEN)-TIMS technique, with Nd + ions used to analyze the Nd isotope composition of a small amount of sample (Wakaki and Ishikawa 2018). The 143 Nd/ 144 Nd ratio was corrected by internal normalization using the exponential law with an 146 Nd/ 144 Nd ratio of 0.7219 (O'Nions et al. 1977) to correct for mass discrimination; 143 Nd/ 144 Nd ratios were corrected for a small systematic bias based on the recommended 143 Nd/ 144 Nd value of standard JNdi-1 (0.512115; Tanaka et al. 2000). The final analytical error for the normalized 143 Nd/ 144 Nd ratio, expressed as 2σ (Additional file 1: Tables S1, Additional file 2:

Water sample preparation and hydrogen and oxygen isotope analyses
Bottom-water samples for H and O isotope analyses were collected at 189 sites at depths of 68-2620 m during the GK19 cruise. Samples were filtered with a 0.45 µm membrane filter and collected in 30 mL glass vials, which were sealed immediately after collection and stored in a cool, dark place. The hydrogen (D/H) and oxygen ( 18 O/ 16 O) isotope ratios of the water samples were determined using cavity ring-down spectroscopy (L2120-I Analyzer; Picarro, USA) at the Atmosphere and Ocean Research Institute, The University of Tokyo, Japan, without chemical preprocessing. Measured H and O isotope ratios were calibrated against an in-house standard and converted into the conventional Vienna standard mean ocean water (V-SMOW) scale (per mil; ‰) based on the following equations: External precision for H-O isotope analyses, based on replicate measurements (n = 38) of the in-house standard, was ± 0.61‰ for δD and ± 0.10‰ for δ 18 O.
Temperature and salinity were measured using the CTD profiler at the bottom-water sampling sites; the pH of the bottom-water samples was measured using a compact pH meter (LAQUAtwin-pH-33; Horiba Scientific, Japan) immediately after samples were collected. The external precision (1σ) of pH values, based on replicate measurements of a standard solution (pH = 7.41), was ± 0.02.

Neodymium isotopes
The εNd values ranged from − 8.2 to − 2.2 on the OT side of the RIA (Fig. 3). Comparison of εNd profiles in the three subareas shows (1) an eastward increase at < 300 m depth in the Iriomote and Ishigaki OT; (2) the occurrence of a water mass with a wide range of εNd values (− 5.9 to − 2.8) at 300-600 m depth; and (3)  Page 10 of 20 Cruz Salmeron et al. Progress in Earth and Planetary Science (2022) 9:44 εNd values (− 5.9 to − 4.8) at > 600 m depth, with two highly negative values of − 7.6 and − 8.2 being recorded in the Yonaguni OT at depths of 848 and 1026 m, respectively.
The εNd values ranged from − 7.0 to − 1.1 on the Pacific side of the RIA (Fig. 3). The Nd isotope profiles of the three subareas on the Pacific side had similar values and trends from surface to deeper waters, with εNd values of − 2.8 to − 1.1 at < 157 m depth (with an outlier of − 4.1 at 125 m in the Iriomote Pacific). At > 222 m depth, εNd values decreased to − 5.5 at 656 m depth, with an outlier of − 7.0 at 497 m at Yonaguni Pacific, and increased to − 4.0 at 1239 m depth (Fig. 3). A value of − 5.3 was recorded at 2264 m depth at the Iriomote Pacific site.

Hydrogen isotopes
On the OT side of the RIA, seawater δD values had similar depth trends for the three subareas (Fig. 3), except at 250-500 m depth in the Ishigaki OT where higher values were recorded. δD values were higher at 117-388 m depth (up to 4.9‰), decreasing steadily to a minimum value of − 1.0‰ at 1093 m depth. Although average δD values were almost constant at greater depths, an outlier of > 1.3‰ was recorded at 1800 m depth.
The three subareas on the Pacific side also displayed similar trends, with values of up to 5.3‰ above 380 m depth. δD values decreased to their lowest value of − 0.7‰ at 976 m depth, before gradually increasing to 1.8‰ at 2128 m depth. A value of − 0.8‰ was recorded at 2620 m depth (Fig. 3).

Oxygen isotopes
On the OT side of the RIA, seawater δ 18 O values had similar depth trends in the three subareas (Fig. 3), with higher values of up to 0.7‰ being recorded above 298 m depth. Below this depth, values decreased to a minimum of − 0.2‰ at 760 m depth, remaining roughly constant at ~ 0.0‰ at greater depths with an outlier of 0.3‰ at 1800 m. The δ 18 O profiles on the Pacific side also had similar trends in the three subareas, with δ 18 O values being higher above 285 m depth (up to 0.7‰) and decreasing to − 0.2‰ at 726-976 m depth before increasing to 2.0‰ at 1850 m depth (Fig. 3).
The εNd values of the KC decrease as it flows to higher latitudes in the OT and then to the southeast of Honshu, as shown by low εNd values of − 7.3 to − 3.9 (< 200 m depth, station 2; Behrens et al. 2018) and − 8.7 to − 5.6 (< 455 m depth, station LM-6/11; Amakawa et al. 2004b;Figs. 1 and 5).
When Northwest Pacific waters flow near Taiwan or into the ECS, they are exposed to the influence of sediment and water with low εNd values derived from the Eurasian continent and Taiwan, as indicated by εNd values as low as − 6.2 at 0-500 m depth (σ θ < 26.4 kg m −3 ) on the OT side of the RIA. Possible low-εNd sources include continental rocks of the Eurasian continent (China) and Taiwan, which are enriched in REEs and have strongly negative εNd values of − 15.2 to − 8.2 (Lan et al. 2002;Li et al. 2013;He et al. 2015). Sediments derived from these rocks are transported to the ocean by numerous rivers, including the Changjiang (εNd = − 11.6; Dou et al. 2012), Yellow (− 13.9 to − 9.6; Meng et al. 2008), and Yangtze (− 14.3 to − 8.7; Meng et al. 2008) rivers on the Eurasian continent and Taiwan rivers (− 11.1; Lan et al. 2002). In contrast, surface water (εNd = − 1.7 to − 1.2) and NPTW (− 3.4 to − 1.1) on the Pacific side of the study area retain the Nd isotope signatures of those from the PS, for which similar εNd values (Fig. 5, Table 2) have been reported; e.g., − 3.2 to − 2.0 (Behrens et al. 2018;Wu 2014).
The observation of some high εNd values on the OT side and low values on the Pacific side of the RIA may result from direct interaction with continental inputs in combination with physical processes, as follows.
(1) Ocean circulation is likely to change near continental zones and shelves owing to local topographical effects and physical processes. For example, areas with narrow bathymetry (i.e., the OT) may promote abrupt hydrodynamic behavior that increases horizontal fluxes to values of up to three times vertical fluxes (Huthnance 1995;Lacan and Jeandel 2005;Lambelet et al. 2016).
(2) Westward-moving eddies to the east of Taiwan trigger KC meandering. Differing εNd values between similar locations and depths have been reported in areas Page 11 of 20 Cruz Salmeron et al. Progress in Earth and Planetary Science (2022) 9:44 under the influence of mesoscale eddies, which play an essential role in vertical and horizontal mixing processes (Grasse et al. 2012). The KC also undergoes periods of high and low transport. At the southeastern tip of Taiwan, a considerable portion of the KC spreads northeastward during low-transport periods, then flows back to the OT between the South Ryukyus (Gawarkiewicz et al. 2011). Drifters released to the east of Taiwan indicate complex spread patterns within the OT and western PS, usually reaching the Pacific side of the study area and, to a lesser extent, as far as Okinawa Jima Island (Gawarkiewicz et al. 2011;Rudnick et al. 2015). Particulates discharged by eastern Taiwan rivers likely follow the same patterns, being deposited on the OT and Pacific sides of the Ryukyus.
(3) Chemical processes controlling the εNd distribution play a significant role in the variable composition of the OT and Pacific sides, as discussed in the following subsections.

Intermediate water
Intermediate water (mainly of 500-700 m depth; σ θ = 26.5-27.0 kg m −3 ) on both the OT and Pacific sides of the RIA has εNd values mainly in the range of − 5.5 to − 4.8, at least 1.7 epsilon units lower than those reported for NPIW (σ θ = 26.67-27.11 kg m −3 ) at similar depths (− 3.1 to − 2.7; Behrens et al. 2018;− 2.8, Wu 2014. Table 2). Although Wu (2014) reported that SCSIW is characterized at similar depths by εNd values ranging from − 4.1 (station LS C3a) to − 3.6 (station SCS C5; Fig. 1, Table 2), there remains a difference of at least 0.7 Fig. 5 εNd profiles for the study area, together with those reported for the northwestern Pacific (Amakawa et al. 2004b;Wu, 2014;Behrens et al. 2018) Page 12 of 20 Cruz Salmeron et al. Progress in Earth and Planetary Science (2022) 9:44 epsilon units that cannot be explained by the influence of SCSIW around the study area. The presence of a low-εNd source is indicated by lower εNd values at depths of > 1000 m in the three subareas on the OT side of the RIA (Figs. 3 and 5; e.g., − 8.2 at 1026 m depth in the Yonaguni OT, − 5.1 at 1545 m in the Iriomote OT, and − 5.9 at 1433 m in the Ishigaki OT). OT circulation is controlled mainly by Pacific Ocean water flowing throughout the western margin of the PS and being carried by the KC from subtropical areas, with εNd values of surface-intermediate water in the PS being greater than − 3.7 (Table 2; Behrens et al. 2018). To a lesser extent, SCS surface-subsurface water also flows into the OT (Chen 2005), with εNd values of − 4.5 to − 2.9 (Wu 2014; Fig. 6). Even if SCS water fully controlled circulation within the OT, it could be responsible only for εNd values above − 4.5, which is not the case, so an external source must provide the lower εNd values.
It is apparent that riverine sediments from Taiwan contribute to some extent to low εNd values down to the bottom of the OT. Studies have indicated that up to 60% of the sediment deposited at 1274 m depth at the southwestern corner of the OT is derived from Taiwan (Diekmann et al. 2008;Dou et al. 2016), although from an oceanographic perspective, sediment from Taiwan cannot cross the KC. As the Yonaguni OT is characterized by the occurrence of intermediate water with the lowest εNd values of all water columns investigated in this study, some processes (e.g., turbiditic fluxes) must contribute vertical transport from low-εNd sources to much greater depths. Sediment fluxes in the southern OT increase with depth, suggesting horizontal advection in mid-depth water. Moreover, aluminosilicate plumes identified at depths of 500-700 m may result from resuspended sediments of the Ilan Shelf and Ridge (Hsu et al. 2004).

Deep water
The  Page 13 of 20 Cruz Salmeron et al. Progress in Earth and Planetary Science (2022) 9:44 Previous studies have shown that εNd values decrease from the PS toward the SCS, as follows: western PS, εNd = − 2.8 at 1500 m (σ θ = 27.58 kg m −3 ) at station 5 of Behrens et al. (2018); Luzon Strait, − 3.5 at 1600 m depth at station P2; northeastern SCS, − 3.9 at 1600 m depth at station LS C3a; and − 4.4 at 1600 m depth at station SCS C5 of Wu (2014; Fig. 5). This confirms that SCSIW flows out through the Luzon Strait into the PS at 350-1350 m depth (Wang and Chen 1998). Subsequently, SCSIW encounters NPIW, resulting in KIW having salinity and εNd values between those of SCSIW and NPIW. Deep waters in the OT are ventilated by intermediate water from the western PS, which enters via the KG at ~ 1100 m depth (Nakamura et al. 2013). However, due to the limited Nd isotope data available for depths > 1500 m (σ θ > 27.35 kg m −3 ) in the OT, it is difficult to confirm whether its Nd isotope composition is similar to that of upper deep water (~ 2000 m, σ θ = 27.66 kg m −3 ) on the Pacific side of the RIA. Several factors may explain the low εNd values, as follows.
Our two εNd values for the deepest sites in the OT (− 5.9 at 1433 m, σ θ = 27.36 kg m −3 in the Ishigaki OT; − 5.1 at 1545 m, σ θ = 27.37 kg m −3 in the Iriomote OT) suggest that the upper layer (immediately above 1500 m, σ θ < 27.35 kg m −3 ) ventilates and imparts a low εNd signature to the lower layer (Fig. 3). It is likely that vertical mixing between intermediate and deeper waters homogenizes the Nd isotope composition below 500 m depth, as indicated by the uniform εNd values (− 5.9 to − 4.8) in the Iriomote and Ishigaki OT. Moreover, salinity and δD and δ 18 O values are similar at > 750 m depth (σ θ > 27.0 kg·m −3 ; Fig. 3). Horibe and Ogura (1968) reported that considerably high δD values characterize not only surface water but also bottom water around the KC path, compared with adjacent open-sea values (Additional file 5: Fig. S3). They suggested that as the KC flows northeastward, the influence of surface and subsurface waters with high δD values reaches deeper waters (e.g., intermediate water had δD values of 4.4‰ at 433 m depth and 3.5‰ at 647 m depth at station 58; Fig. 1). Surface and bottom-water mixing due to a KC-induced flow may reach > 500 m depth. Our results indicate considerably high δD values in deep waters of the study area, with average values of 0.4‰ (σ = 0.8‰) at 1500-1938 m depth on the OT side and 0.7‰ (σ = 0.9‰) at 1907-2620 m depth on the Pacific side ( Fig. 3 and Additional file 5: Fig. S3). Therefore, the presence of low εNd values even below the intermediate-water domain must be promoted to some extent by vertical mixing, which decreases with depth as shown by the differences in δD values below 750 m depth (Fig. 3, Table 3).
Another possible explanation of the negative εNd values in the deeper areas is that turbidites, reported at ~ 1400 m depth to the east of Taiwan in the southern OT (Huh et al. 2006), carry sediment with highly negative εNd values to deeper areas. In addition, near continents or marginal seas, the εNd distribution is more complex than that in the open ocean. Therefore, near the continents, the εNd value of a water mass cannot be explained only on the basis of conservative Page 14 of 20 Cruz Salmeron et al. Progress in Earth and Planetary Science (2022) 9:44 behavior and water mass mixing/advection, so local sources must be involved (Grasse et al. 2012). Boundary exchange and benthic Nd fluxes may affect the εNd values of seawater in complex ways (Lacan and Jeandel 2005;Abbott et al. 2015). However, as we do not have data to determine whether this process is critical, it is not discussed further here. Deep-water (2000-3000 m depth, σ θ = 27.65-27.73 kg m −3 ) εNd values decrease from east to west in the PS from station 5 (− 3.5 to − 2.8) and station 4 (− 3.8 to − 3.5) of Behrens et al. (2018) to station P2 (− 4.2 to − 3.9) and station WPS N18-05 (− 4.3 to − 3.6) of Wu (2014;Figs. 1, 5). Deep water with low εNd signatures of the western PS then flows northeastward to higher latitudes along the western PS margin, as shown by data for the Iriomote Pacific (− 5.3; this study) and station LM-6/11 (− 6.0; Amakawa et al. 2004b). Our results indicate that the western PS is more affected by UCDW with low εNd values (i.e., − 6.5 to − 4.6; Amakawa et al. 2013) than the eastern PS below the intermediate layer, as shown by lower εNd values in the Iriomote Pacific and Ishigaki Pacific. UCDW flows in a SE-NW direction, bordering Papua New Guinea. About 10 × 10 6 m 3 s −1 (= 10 Sv) of seawater enters the PS via the East and West Caroline basins and the East Mariana Basin. Of the total UCDW influx (10 × 10 6 m 3 s −1 ) into the PS, at least 7.5 × 10 6 m 3 s −1 of UCDW does not flow out of the PS, as indicated by the western boundary deep-water flowing to the north at 2650 m depth to the east of Okinawa Jima Island (Kawabe and Fujio 2010). Therefore, the remaining flux must be upwelled, imprinting low εNd signatures on the upper layers by vertical upwelling and mixing.

Hydrogen and oxygen isotopes
The Our δD data are snapshot values rather than long-term data, such as the JODC salinity and temperature data, so they may not necessarily represent average δD values of water masses. However, as the same is true for most of the δD data in previous studies (Benetti et al. 2017;Belem et al. 2019), it is appropriate to use our δD data here in delineating water masses around the RIA. Therefore, the argument based here on δD data is less rigorous than that based on Nd and δ 18 O data.
On the OT side of the study area, we measured the δD and δ 18 O values of one sample from KSW in Ishigaki OT. The values were 4.5‰ and 0.5‰ (Additional file 1: Table S1), with average values for KTW of 3.6‰ and 0.4‰, respectively. On the Pacific side of the RIA, surface water had average δD and δ 18 O values of 3.3‰ and 0.4‰, and NPTW had average values of 3.4‰ and 0.4‰, respectively (Fig. 3, Table 3). These values indicate that high δ 18 O values characterize the KC, as shown in previous studies (e.g., Horibe and Ogura 1968). Our δD values are similar to those reported previously for NWP waters at low latitudes. For example, Horibe and Ogura (1968) reported δD values of surface water (2.7‰-3.7‰ at < 130 m depth) and tropical water (2.8‰-3.4‰ at 130-174 m) from east of Luzon Island near the KC domain in the western PS (station 85, Fig. 1).
The intermediate-water core on the OT side had average δD and δ 18 O values of 1.0‰ and 0.2‰, which are lower by ~ 2.1‰ and ~ 0.2‰ than those of WNPCW, respectively (Table 3) On the OT side of the study area, the δD and δ 18 O values for KSW (4.5‰ and 0.5‰,) KTW (2.4‰-4.8‰ and 0.1‰-0.6‰), and WNPCW (1.0‰-5.2‰ and 0.1‰-0.7‰) overlap, so the three water masses cannot be distinguished using these isotopes. This is also the case for the Pacific side of the RIA, where the ranges in δD and δ 18 O values of surface water, NPTW, and WNPCW also overlap.
In contrast, average δD and δ 18 O values differ between subsurface water (δD ≥ 3.1‰ and δ 18 O ≥ 0.4‰) and intermediate water (δD ≤ 1.0‰ and δ 18 O ≤ 0.2‰) on the OT side of the RIA (Table 3), as is also the case for Page 15 of 20 Cruz Salmeron et al. Progress in Earth and Planetary Science (2022) (Table 3).

Relationships among temperature, salinity, and hydrogen and oxygen isotopes
The δD-δ 18 O cross-plots for surface-subsurface water (Fig. 7) indicate the influence of runoff near continental areas, with the intercepts representing freshwater inputs (Benetti et al. 2017). The δD and δ 18 O values of seawater decrease as the influence of freshwater and precipitation increases because freshwater and meteoric water are depleted in D and 18 O. Salinity-δD/δ 18 O and temperature-δD/δ 18 O relationships (Additional file 6: Fig. S4) have been used to characterize water masses, delineate water-mass mixing, and identify the influence of precipitation/evaporation on surface waters. The three subareas on the Pacific side of the RIA display similar trends in cross-plots for surface to deep waters (Fig. 7c), with none of the three subareas forming a discrete cluster of values. This reflects a more stable circulation, which could be attributed to the long-distance high-speed KC and the complex fluxes generated. Intermediate water in Ishigaki Pacific (orange rhombuses) subarea is clustered to the left of the Iriomote and Yonaguni subareas, as it is characterized by the lowest salinities (Additional file 6: Fig. S4e). In addition, in Fig. 4, the Ishigaki Pacific corresponds to the low salinity curve. This correspondence is likely to be attributed to the difference in salinities between NPIW and SCSIW. Therefore, the Ishigaki Pacific better retains the physical properties of NPIW, which is represented by a distinctive curve with the lowest salinities at about 600-700 m water depth, between σ θ = 26 and 27 kg m −3 . Intermediate water on the OT side of the RIA is more saline by ~ 0.11 than that on the Pacific side, for which there are two possible explanations. First, PDW is more saline than intermediate water, and as it approaches the RIA with decreasing depth, it may mix with the intermediate water before entering the OT. However, this is not plausible, because deep water in the OT is topographically isolated from western PDW (UCDW and North Pacific deep water; Nakamura et al. 2013;Behrens et al. 2018). Although the ventilation process in the OT is not well understood, σ θ at > 1500 m depth on both the OT and Pacific sides of the RIA (σ θ = 27.4 kg m −3 ) indicates that the deeper intermediate water flowing from the northern PS into the OT via the KG at ~ 1100 m depth  Trough (a, b) and Pacific (c) sides of the study area. a Cross-plots for < 433 m depth; b, c plots for seawater collected from all depths. Outliers for which δ 18 O and/or δD values are not within the mean ± 2σ range for individual water masses are not shown. For abbreviations, see Fig. 3 Page 16 of 20 Cruz Salmeron et al. Progress in Earth and Planetary Science (2022) 9:44 may ventilate as far as the southern OT (at 2000 m depth; Nakamura et al. 2013). The second and more plausible explanation for the higher salinity is that the SCSIW (Wang and Chen 1998;Chen 2005), which is modified KIW becoming more saline in the southern SCS (salinity 34.6; Qu et al. 1999), flows out through the Luzon Strait at 350-1350 m depth (Wang and Chen 1998) with salinity of 34.4-34.6 (Gong et al. 1992;Qu et al. 2000). The mixed water comprising SCSIW (55%; wine-colored symbols in Fig. 3) and NPIW (45%) enters the OT via the TY (Nakamura et al. 2013). A salinity front separating western high-salinity water from eastern low-salinity water has been detected to the east of Taiwan at 21.45°N (Chen 2005), suggesting that SCSIW interacts with NPIW in this area. The KC has a turbulent flux that causes surface and subsurface waters to reach down to 500 m depth (Horibe and Ogura 1968). Bottom-water upwelling has been reported from the southern OT (Nakamura et al. 2013). Therefore, the minimum salinity layer within the OT may increase through vertical mixing with higher salinity of upper and lower layers. From surface to bottom water, the Ishigaki OT has an intercept of 0.82, Iriomote OT 0.27, and Yonaguni OT 0.15. The very different pattern for Yonaguni OT is caused by the lack of data for surface water, resulting in a slope considerably lower than those of the Iriomote and Ishigaki OT (Fig. 7b). The δ 18 O-δD relationship at < 433 m depth on the OT side indicates that the Ishigaki OT has the highest intercept (δD = 2.56‰), followed by Iriomote (2.04‰) and Yonaguni OT (0.97‰). This indicates the more significant influence of riverine discharge in the Yonaguni and Iriomote subareas than in the Ishigaki OT (Fig. 7a). Although the proximity of Yonaguni OT to Taiwan inputs is clear, additional processes must control vertical circulation, as shown by the low εNd values and similar salinities and δD and δ 18 O values below 750 m depth. On the OT side, temperature and salinity are lower by ~ 1.0 °C-1.5 °C and ~ 0.03-05 in the Yonaguni OT than in the Iriomote OT and Ishigaki OT, respectively. In contrast, the Iriomote OT shows higher salinities at > 250 m depth (Fig. 3). However, Fig. 4 shows similar θ-salinity curves for the three subareas at σ θ = 23-25 kg m −3 . These results suggest that changes in θ and salinity within the three subareas may be related to vertical variations in physical properties (e.g., a vertical shift of the thermocline). Thermocline shifts are likely to be driven by mesoscale eddies (Chen and Chen 2020), which in turn affect KC behavior (e.g., meandering and upwelling) and sea surface height. Areas with sea surface height anomalies are related to eddies, which promote geostrophic circulation (Itoh and Yasuda 2009). The presence of eddies from the east coast to the northeast of Taiwan (Zhang et al. 2001;Gawarkiewicz et al. 2011) may cause vertical shifts in the physical properties of water around the Yonaguni OT.
Another possibility is as follows. The cross-plots indicate that the δD and δ 18 O values of intermediate water in the Yonaguni OT overlap with those of the upper and lower layers. Such a trend is not confirmed for the Iriomote and Ishigaki OT subareas (Fig. 7c-f ). The interaction of the KC with the topography of the southern OT may promote cross-shelf and downslope transport, enhancing vertical mixing in this area (Huh et al. 2006). Therefore, a different circulation pattern should govern the Yonaguni OT, causing intermediate water to ventilate down to bottom water.

Conclusions
The water-mass structure was delineated on the OT and Pacific sides of the southern RIA (off Yonaguni Jima, Iriomote Jima, and Ishigaki Jima islands) based on εNd values of benthic foraminifers, δD and δ 18 O values, and physical properties (temperature and salinity) of seawater. δD and δ 18 O values were determined for seawater collected at 69-1938 m depth on the OT side and 68-2620 m on the Pacific side of the study area. εNd values were determined for benthic foraminifers collected at 69-1545 m depth on the OT side and 74-2264 m on the Pacific side of the study area. The main results of this study are as follows.
1. Riverine sediment discharging into the ECS, including the southwestern OT, is derived mainly from old continental rocks of the Eurasian continent and Taiwan with low εNd values. The low Nd isotope signature is imprinted in water masses flowing within the OT, resulting in a wide range of εNd values of − 8.2 to − 2.2, with slightly lower values than those of the Pacific side for seawater at the same depth. 2. The OT-side circulation is disturbed when the highspeed KC interacts with the complex and narrow bathymetry. As a result, there are overlapping εNd values between different seawater layers. The εNd values on the OT side suggest that additional processes controlling the Nd behavior of seawater occur in individual seawater layers. Surface and subsurface waters are influenced by Taiwanese river discharge associated with temporospatial variations in oceanographic conditions such as KC meandering; intermediate water has the lowest εNd values (as low as − 8.2), possibly resulting from sediment plumes and turbiditic fluxes; bottom water has εNd values similar to those of intermediate water, suggesting active upwelling and vertical mixing. 3. Unlike the OT side, the Pacific side has betterdefined εNd profiles, with more pronounced changes