Subtropical Mode Water in a recent persisting Kuroshio large-meander period: part II—formation and temporal evolution in the Kuroshio recirculation gyre off Shikoku

Since August 2017, the large meander (LM) path of the Kuroshio south of Japan has continued to date for more than five years, with a local Kuroshio recirculation gyre off Shikoku (RGOS). The formation and temporal evolution of Subtropical Mode Water (STMW) in the isolated RGOS were examined by using Argo profiling float and shipboard observation data. Since the late winter of 2018, which was the first winter of the current LM event, multiple STMW layers with different densities overlapping with each other have been observed for more than three years inside the RGOS. The deeper STMW layer of 18.0 °C, which was formed in the late winter of 2017 before the current LM event began, has survived at least until September 2021 while its thickness and horizontal extent has decreased gradually. We may have succeeded in discovering the longest survived STMW ever observed due to a suitable environment for STMW conservation inside the RGOS. On the other hand, the shallower STMW layer(s) of 19.0–19.8 °C, which was originally formed over the deeper STMW layer in the late winter of 2018, has been repeatedly renewed or dissipated in the following years. The oxygen utilization rate at the core of the deeper STMW at a depth of 500 dbar was estimated to be 4.6 μmol kg−1 yr−1, which is believed to be explained solely by oxygen consumption due to remineralization.


Introduction
The North Pacific Subtropical Mode Water (STMW; Masuzawa 1969), which is characterized by a thick thermostad with a potential temperature (θ) between 16 and 19.5 °C (Oka 2009), is formed as deep winter mixed layers (ML) south of the Kuroshio and Kuroshio Extension and north of 28ºN (Suga and Hanawa 1990;Oka and Suga 2003).The colder variety of STMW with θ = 16 − 18 °C, which accounts for more than 80% of the total STMW volume, is formed east of 140ºE, i.e., south of the KE, while the warmer variety of STMW with θ = 18 − 19.5 °C is formed mainly south of the Kuroshio to the west of 140°E (Oka et al. 2021).After formation, both varieties are advected southwestward by the mean flow in the southern part of the Kuroshio recirculation gyre (Suga and Hanawa 1995a) and by mesoscale eddy activity (Uehara et al. 2003;Nishikawa et al. 2010;Xu et al. 2016) to spread as far as 20°N to the south and just east of Taiwan to the west within 2-3 years (Oka 2009).STMW transmits signals resulting from air-sea interaction such as temperature anomalies to the ocean subsurface through its subduction process (Yasuda and Hanawa 1997;Oka and Qiu 2012;Newman et al. 2016).These anomalies then reappear at the sea surface in the formation region (Hanawa and Sugimoto 2004) or in a remote area (Sugimoto and Hanawa 2005), possibly affecting atmosphere.
The formation and advection of STMW are significantly affected by the Kuroshio path variations south of Japan among three typical paths: the large meander (LM) path and the offshore and nearshore non-large-meaner (NLM) paths (Fig. 1; Kawabe 1995).When the Kuroshio takes an LM path, the westward advection of the colder variety of STMW formed east of 140°E is inhibited, drastically reducing its volume on the western side of the LM (Suga and Hanawa 1995a, b;Oka et al. 2021).In addition, when the Kuroshio takes a meandering path (i.e., LM and offshore NLM paths), the warm STMW exceeding 19 °C is formed in the local Kuroshio's recirculation gyre off Shikoku (RGOS), which is separated from the recirculation gyre south of the KE and isolated to the southwest of the meandering path of the Kuroshio (Nishiyama et al. 1980(Nishiyama et al. , 1981;;Oka 2009;Sugimoto and Hanawa 2014;Fig. 1).By using Argo profiling float data from 2005 to 2011, Sugimoto and Hanawa (2014) showed a double-layer structure of STMW in the isolated RGOS, which composed of the warmer STMW exceeding 19 °C developed in the RGOS and the colder one formed east of 140ºE in previous winters and advected from there.However, each meandering path event during 2005 − 2011 (i.e., the LM event in 2005 and offshore NLM events in 2007 and 2009) lasted for less than one year, and the behaviors of STMWs in the isolated RGOS over multiple years have remained unclear.
In August 2017, the Kuroshio LM occurred 12 years after the previous event in July 2004 − August 2005, and has continued to date (as of February 2023), which is the longest record in the observation history since 1950 (Qiu and Chen 2021).This is also the first multi-year LM event in the era of satellite altimeter and Argo, which enabled us to investigate water mass processes in relation to currents and mesoscale eddies in detail.To clarify the formation and temporal evolution of STMW in the isolated, multi-year RGOS, we analyze Argo profiling floats, satellite-derived altimeter data, and shipboard observations in this paper.Section 2 explains the data and the analysis method.Section 3 describes the STMW structure and its time evolution in the RGOS observed by Argo float and shipboard observations.Section 4 discusses the estimated oxygen utilization rate in the STMW.Finally, Sect. 5 gives our summary.

Data and method
We used temperature (T) and salinity (S) profiles from Argo floats inside the RGOS in 2017-2021, which were downloaded from the ftp site of the Argo Global Data Assembly Center (ftp:// usgod ae.org/ pub/ outgo ing/ argo, ftp:// ftp.ifrem er.fr/ ifrem er/ argos, https:// doi.org/ 10. 17882/ 42182) and edited as outlined in Oka et al. (2007).Each profile was vertically interpolated onto a 1-dbar grid using the Akima spline (Akima 1970).To obtain vertical cross sections in the RGOS, we averaged profiles in every 10-km bin from the RGOS center, which is detected as the highest sea surface height (SSH) from satellite altimeter (detailed below).We also used T, S, and dissolved oxygen profiles from two research cruises.The first one is observations by a conductivity-temperature-depth profiler (CTD) with a RINKO rapid-response dissolved oxygen sensor and expendable bathythermograph (XBT) across the RGOS, which were conducted in 16-17 June 2018 in the KS1805 cruise of R/V Keifu-maru by the Japan Meteorological Agency.Since XBT only measures T, we estimated S from XBT T profile based on the T − S relationship obtained from CTD (Stommel 1947;Goes et al. 2018).The second one is intensive shipboard observations during 28-29 May 2021 in the KS-21-9 cruise of the R/V Shinsei-maru that we carried out to explore the STMW in the RGOS.T and S profiles were obtained at intervals of 10 min in latitude and longitude crossing the RGOS meridionally and zonally by expendable CTD (XCTD).T, S, and dissolved oxygen profiles near the center of the RGOS were also obtained by a conductivitytemperature-depth-oxygen profiler (CTDO 2 ).
We calculated θ, potential density ( ), and potential vorticity ( Q ) from the Argo and shipboard T/S profiles.Here, Q is defined as Q = gf ∕ p , neglecting relative vorticity (Qiu et al. 2006), where g is the gravity acceleration, f is the Coriolis parameter, and p is pressure.STMW was detected as layers with Q < 1.5 ×10 −10 m −1 s −1 and θ = 16 − 20.5 °C.The core of STMW was defined as a local vertical minimum of Q in each profile.We defined mixed layer depth as the depth at which increases by 0.03 kg m −3 from its 10-dbar depth value (Boyer Montégut et al. 2004;Oka et al. 2007).After the saturated dissolved oxygen concentrations were estimated for each depth from the observed T and S (Weiss 1970), apparent oxygen utilization (AOU) was computed as To assess conditions of the Kuroshio and the RGOS and detect the RGOS center, the absolute SSH data from the satellite altimeter, Global Ocean Gridded Level 4 Sea Surface Heights (https:// doi.org/ 10. 48670/ moi-00148, https:// doi.org/ 10. 48670/ moi-00149), produced and distributed by the Copernicus Marine and Environment Monitoring Service (CMEMS; http:// marine.coper nicus.eu), were also used.The

Temporal evolution of STMW structure and long survival of deeper STMW in the RGOS
In late March 2017 before the current LM began, the Kuroshio took an offshore NLM path, and the RGOS was being established (SSH map in Fig. 2a).An Argo float located near the center of the RGOS in formation observed a deep winter mixed layer approaching 500 dbar with θ = 18.0˚C, S = 34.89,and = 25.18kg m −3 (Fig. 2a), which were close to the properties of the colder variety of STMW formed south of the Kuroshio Extension (Oka et al. 2021).In August 2017 when the current LM began, a vertically uniform STMW layer with a core θ of 17.9 °C, a core S of 34.89, and a core of 25.21 kg m −3 was found in the subsurface at depths of 150-600 dbar (Fig. 2b).From April to September of 2017, the STMW layer existed over 200 km from the RGOS center with a maximum thickness of 400 dbar (Fig. 3).Hereafter, the STMW formed in the late winter 2017 is referred to as the STMW-2017, and STMWs formed in other years are also denoted in the same manner.While Argo floats were not distributed near the center of the RGOS from fall of 2017 to spring of 2018 (Fig. 3), shipboard observations conducted in June 2018 captured a double layer structure of STMW across the RGOS (Figs. 4  and 5).The shallower STMW layer at depths of 100 − 300 dbar with a width of 200 km had a core θ of 19.0 °C, a core S of 34.93, and a core of 24.98 kg m −3 , while the deeper STMW layer at 350 − 620 dbar with a width of 100 km had a core θ of 18.0 °C, a core S of 34.89, and a core of 25.20 kg m −3 .The deeper STMW layer is believed to be formed in the previous year because the water properties are almost the same with those of the STMW-2017 (Fig. 2a, b).Therefore, the shallower STMW layer was probably formed in the late winter of 2018 (STMW-18).STMW-2018 was warmer than STMW-2017 by 1 ˚C, having properties close to the warmer variety of STMW formed south of Japan (Oka et al. 2021).The double STMW layers were also observed by an Argo float at depths of 100-300 and 350-600 dbar in August 2018 (Fig. 2c).The shallower and deeper STMW were distributed up to 150 km and 90 km from the RGOS center in the summer of 2018, respectively (Fig. 3).Such a double STMW layer structure was maintained until the end of 2018 while their horizontal extent gradually decreased.
Similar multiple structure of STMW has been frequently observed (e.g., Taneda et al. 2000;Oka et al. 2011;Liu et al. 2017).The earlier study inferred that the multiple structure south of the Kuroshio Extension was formed due to interleaving of STMWs with slightly different densities (Oka et al. 2011).A shipboard survey in the western North Atlantic observed pronounced double STMW layers with a core θ difference of 1.7 °C, which were formed due to the vertical alignment of two warm-core rings with low-Q water pinched off from the Gulf Stream meander (Belkin et al. 2020).On the other hand, profiling float observations in the western North Pacific revealed that in a westward propagating anticyclonic eddy, a warmer STMW was formed south of Japan, overlying a colder STMW formed in the region south of the KE and advected from there within the eddy (Xu et al. 2017;Liu et al. 2019).Such formation of multiple layers was similar to that in the RGOS.
In early March 2019 (Fig. 2d), the upper part of the STMW-2018 was entrained into the winter mixed layer that developed down to a depth of 240 dbar, and new STMW, i.e., STMW-2019, was being formed in the mixed layer, establishing a triple STMW layer structure.The θ, S, and of the winter mixed layer, i.e., shallow STMW layer in progress of formation, were 19.4 °C, 34.86, and 24.83 kg m −3 , respectively (Fig. 2d).The core θ, S, and of middle and deep STMW layers at depths of 309 and 493 dbar were 19.0 °C, 34.94, and 24.98 kg m −3 , and 18.0 °C, 34.91, and 25.21 kg m −3 , which were almost the same as those of the STMW-2018 and STMW-2017, respectively.In November 2019, the core of the middle STMW-2018 was mostly disappeared, and the double STMW layer structure was again established inside the RGOS (Fig. 2e).The shallow STMW-2019 had a core θ of 19.1 °C, a core S of 34.87, and a core of 24.90 kg m −3 .The deep STMW-2017 with a core θ of 17.9 °C, a core S of 34.90, and a core of 25.21 kg m −3 , the thickness and radius of which is about 250 dbar and 70 km, was well preserved within the RGOS, although it was difficult to see in the summer 2019 due to lack of profiles near the RGOS center (Fig. 3).
In May 2020, triple STMW layers consisting of the shallow STMW with a core θ, S, and of 19.5 °C, 34.81, and 24.76 kg m −3 , the middle STMW-2019 with 19.1 °C, 34.89, and 24.92 kg m −3 , and the deep STMW-2017 with 18.0 °C, 34.91, and 25.20 kg m −3 were observed near the RGOS center (Fig. 2f).The shallow STMW was believed to be formed in the late winter 2020 (STMW-2020) as a result of the renewal of the upper part of the STMW-2019 by the entrainment into the winter mixed layer.The middle STMW-2019 had a thickness of 150 dbar and spread within 120 km from the RGOS center (Fig. 3).The deep STMW-2017 still survived while its thickness and radius decreased to 100 dbar and 50 km, respectively (Figs. 2f and 3).
The Kuroshio LM became unstable after June 2020, detaching an intense cold eddy in October 2020 (Fig. 6).The detached cold eddy migrated westward for four months and merged with the Kuroshio south of Kyushu in February 2021 again.During the migration period, the RGOS was transformed significantly.How did such transformation influence the STMW structure in the RGOS, especially the deep STMW-2017?To observe the STMW structure after the RGOS deformation event, we carried out intensive CTDO 2 and XCTD observations during the KS-21-9 cruise of R/V Shinsei-maru in May 2021 (Figs.7a), and succeeded in observing triple STMW layers near the center of the RGOS (Figs. 7 and 8).The shallow STMW layer with a thickness of 250 dbar had a core θ of 19.8 °C, a core S of 34.82, and a core of 24.67 kg m −3 at a depth of 247 dbar (Fig. 8).Although two PV minima were found at depth of 82 and 247 dbar within the shallow STMW layer, they were considered within the same layer because the water properties at each depth of PV minima were quite similar.As these properties were different from those of the STMW-2020, the shallow STMW layer must have been newly formed in the late winter 2021 (STMW-2021).The middle STMW layer at depths of 350-400 dbar had a core θ of 19.1 °C, a core S of 34.87, and a core of 24.89 kg m −3 , which agree with those of  (18.0 °C, 34.87, and 25.18 kg m −3 ) were close to those of STMW-2017 (Fig. 2b).Thus, the STMW-2017 was found to survive for more than four years while its thickness and width had shrunk to 80 m and 40 km, respectively (Fig. 7).Furthermore, in September 2021, STMW with a core θ, S, and of 17.9 °C, 34.84, and 25.18 kg m −3 at a depth of 522 dbar was captured by an Argo float located 30 km from the RGOS center (Figs. 3 and 9).Since these properties were also similar to those of STMW-2017, it is expected that STMW-2017 has survived until at least September 2021.Most of STMW reaches the western boundary in two years after its formation, and then the properties dissipate (e.g., Oka 2009).A series of results indicate that the isolated RGOS is a suitable environment for STMW conservation, and we may have succeeded in observing the longest survived STMW.

Estimation of oxygen utilization rate
Time series observations of STMW-2017 in the RGOS for four years provide a unique opportunity to explore the oxygen utilization rate (OUR) in the ocean.AOU has been used to estimate the age of water masses (e.g., Suga et al. 1989;Oka and Suga 2005).It tends to increase in time due to the oxygen consumption through remineralization of organic matter after isolation of the water mass from the atmosphere (Jenkins 1980) and due to mixing with surrounding older water.OUR is an important indicator for biological processes in the interior ocean, although it is difficult to precisely estimate OUR since oxygen concentrations also fluctuate due to physical processes such as mixing.
Fortunately, the STMW-2017 had been trapped within the RGOS and isolated from its surroundings as described in Sect.3. The shipboard observation in June 2018 demonstrated that the core of the deeper STMW-2017 had an AOU of 18.4 μmol kg −1 and Q of 0.05 × 10 -10 m −1 s −1 at a depth of 545 dbar (orange curve in Fig. 5), while the shipboard observation in May 2021 showed that the core of the deep STMW-2017 had an AOU of 32.1 μmol kg −1 and Q of 0.348 × 10 -10 m −1 s −1 at a depth of 520 dbar (orange curve in Fig. 8).Since STMW-2017 observed in both 2018 and 2021 was characterized by a low core Q , it would be lit- tle affected by diffusion.Therefore, the AOU difference of 13.7 μmol kg −1 can be explained solely by oxygen consumption due to remineralization, and is divided by three years to yield an OUR of 4.6 μmol kg −1 yr −1 at depth of 500 dbar.This OUR is consistent with those from previous works: 4.5 μmol kg −1 yr −1 averaged at depths of 200-1000 m, which was derived from the observed AOU and the water age based on chlorofluorocarbons and 14 C in the western North Pacific by Feely et al. (2004), and 4.2 ± 0.27 μmol kg −1 yr −1 at depths of 100-500 m, which was estimated using water age based on tritium and helium isotope observation in the South China Sea by Xie et al. (2021).
While most of previous OUR estimates, including those mentioned above, have relied on chemical tracers, Billheimer et al. (2021) recently succeeded in estimate of OUR using high resolution oxygen data from biogeochemical profiling floats deployed under the CLIVER Mode Water Dynamics Experiment (CLIMODE; Marshall et al. 2009) in the Sargasso Sea.They also clarified the depth dependence and seasonal variation of OUR, which are valuable for diagnosis of the carbon cycle in the Sargasso Sea.Currently, biogeochemical floats with an oxygen sensor are being expanded in the North Pacific region.Their data are expected to reveal the detailed spatial-temporal distribution and variability of OUR and contribute to our improved understanding of the carbon cycle in the North Pacific.

Summary
Argo floats and shipboard observation data in 2017 − 2021 were analyzed to investigate the formation and temporal evolution of STMW in the isolated RGOS during the current Kuroshio LM period that has lasted since August 2017.In the late winter of 2017 before the current LM event, when the Kuroshio took an offshore NLM path, vertically uniform STMW-2017 of 18.0 °C with a thickness of 500 dbar was formed in the RGOS.In the late winter of 2018, the upper part of the STMW-2017 was entrained into the winter mixed layer, which newly became STMW-2018 of 19.0 °C.Since then, the STMW inside the RGOS has maintained a multi-layer structure for more than three years.The shallower STMW layer(s) originated in 2018 was partially renewed or dissipated in the following years.The deeper STMW-2017 layer has survived at least until September 2021 for 4.5 years, while its thickness and width had shrunk, particularly due to the transformation of the RGOS in October 2020-February 2021 in association with the migration of intense cold eddy detached from the LM.We may have succeeded in observing the longest survived STMW due to a suitable environment for STMW conservation in the RGOS.
The core of STMW-2017 was observed to have AOU values of 18.4 μmol kg −1 in June 2018 and 32.1 μmol kg −1 in May 2021, which yielded an OUR of 4.6 μmol kg −1 yr −1 at a depth of 500 dbar.This value agrees with previous estimates derived from chemical tracers, and is believed to be explained solely by oxygen consumption due to remineralization because the STMW had been trapped within the RGOS and little diffused.
The Kuroshio LM is still ongoing, which enables us to monitor the temporal evolution of the multi-layer STMW structure in the RGOS.As the deployment of Argo floats with an oxygen sensor is being expanded, we expect to reveal more detailed spatial-temporal distribution and variability of not only STMW but also OUR.Furthermore, it will be an interesting topic to investigate how the multi-layer STMW structure in the RGOS, which has been maintained during the current LM period, will change after the LM terminates.The southwestward advection of the warm STMW during the current LM period and afterward is also an interesting theme, as the warm STMW formed south of Japan during the previous LM period in 2005 has been reported to be advected to just east of Taiwan (Oka 2009).

Fig. 1
Fig. 1 Sea surface height (SSH) maps in March of 2020 (top), 2012 (middle), and 2016 (bottom) corresponding to the periods of the LM, offshore NLM, and nearshore NLM paths of the Kuroshio, respectively

Fig. 2
Fig. 2 Vertical profiles of Q (green), S (red), θ (blue), and (black) obtained by Argo floats (upper panels) on a 25 March 2017, b 13 August 2017, c 13 August 2018, d 1 March 2019, e 4 November 2019, and f 17 May 2020 at the position indicated by a white dot in the SSH maps (lower panels).In the upper panels, vertical green dashed line denotes Q = 1.5 ×10 −10 m −1 s −1 , while horizontal black dashed line represents the mixed layer depth

Fig. 4 a
Fig. 4 a Locations of CTD (white circle) and XBT (black cross) stations in the KS1805 cruise.Black contours with color denote SSH on 16 June 2018.b-d Cross sections of (b) Q , (c) θ, and (d) along the observation line indicated in (a).The "v" symbols in black and red indicate XBT and CTD observation points, respectively

Fig. 6
Fig. 6 Monthly SSH maps south of Japan during a LM deformation period from January 2020 to June 2021.The detached cold eddy is enclosed by black circle

Fig. 7 a
Fig. 7 a Locations of CTD (white circle) and XCTD (black cross) stations in the KS-21-9 cruise.Contours with color denote SSH on 29 May 2021.(b-d) Vertical cross section of b Q , c θ, and d along the observation lines indicated in (a).The "v" symbols in black and red indicate XCTD and CTD observation points

Fig. 9
Fig. 9 Vertical profiles (upper panel) and the SSH map (lower panel) on 27 September 2021, otherwise following Fig. 2