Temporal variations of net Kuroshio transport based on a repeated hydrographic section along 137°E

Temporal variations of net Kuroshio transport are examined for 1972–2018 based on a repeated hydrographic section along 137°E, which is maintained by the Japan Meteorological Agency. The net Kuroshio transport obtained by integration of geostrophic current velocity relative to 1000 dbar depth fluctuates on inter-annual and decadal timescales. The predominant timescale of the net Kuroshio transport changes with time; the inter-annual variation is pronounced in 1972–1990 and 2000–2018, and the decadal variation is detected only before 2000. We find that a winter wind stress curl variation in the central North Pacific which reflects meridional movements of the Aleutian Low and intensity fluctuations of the North Pacific subtropical high on an inter-annual timescale and intensity fluctuations of the Aleutian Low on a decadal timescale, causes the net Kuroshio transport variation. In addition to the inter-annual and decadal variations, we further pointed out a bi-decadal-scale variation of the net Kuroshio transport and its possible link to the Aleutian Low intensity fluctuation. Moreover, our results indicate that during large net Kuroshio transport, sea surface temperature around the Kuroshio and Kuroshio Extension region tends to increase, resulting in vigorous upward sensible and latent heat release.


Introduction
The Kuroshio is the western boundary current of the North Pacific subtropical gyre. It flows northeastward along the south coast of Japan and transports large amounts of heat from the western tropical Pacific to the mid-latitude North Pacific. The Kuroshio transport can substantially influence regional and basin-scale climate variability by modulating sea surface temperature (SST) in the western boundary region (e.g., Qiu 2000;Wu et al. 2012). Improved satellite observations and numerical simulations have revealed significant influences of the Kuroshio and its extension, the Kuroshio Extension (KE) on the overlying atmosphere (e.g., Kuwano-Yoshida and Minobe 2017;Nakamura et al. 2012;Nonaka and Xie 2003;Sugimoto 2014;Sugimoto et al. 2017;Tanimoto et al. 2011).
Since the Kuroshio is a part of the wind-driven circulation in the North Pacific subtropical gyre, many authors have investigated causes of Kuroshio transport variations, focusing on basin-scale wind stress curl (WSC) fields (e.g., Deser et al. 1999;Yasuda and Kitamura 2003). The wind stress (WS) field over the mid-to-high latitude in the North Pacific largely reflects the Aleutian Low (AL) activity in the northern North Pacific (e.g., Trenberth and Hurrell 1994) (Fig. 1). Past works pointed out the AL shows two types of variations; one is intensity variations and the other is meridional position changes (e.g., Sugimoto and Hanawa 2009). The AL variations have significant impacts on the Kuroshio transport as a result of oceanic baroclinic responses (e.g., Sugimoto et al. 2010;Yasuda and Hanawa 1997). The North Pacific subtropical high (NPSH) is distributed in the eastern North Pacific (Fig. 1). Recent study 1 3 of Qiu et al. (2020) pointed out that the WSC changes around the NPSH region have significant influences on the KE path state. However, roles of the NPSH on the Kuroshio transport have not been revealed yet.
The Kuroshio transport varies across a wide range of timescales from inter-annual to bi-decadal (e.g., Qiu and Joyce 1992;Sugimoto et al. 2010). However, the length of the Kuroshio transport time series used in the previous studies is about 35 years at most, which is not sufficient to investigate low-frequency variations like decadal and bi-decadal variations. The Japan Meteorological Agency (JMA) started hydrographic observations along 137°E across the Kuroshio in 1967 (Fig. 2). This repeated hydrography is powerful tool to understand a temporal variation of the Kuroshio transport. In this study, using the 137°E sections, we produce a time series of Kuroshio transport which is 10 years longer than that in most existing studies based on integration of geostrophic current velocity relative to 1000 dbar, and then examine causes of the Kuroshio transport variation in terms of the AL and NPSH. We also consider the impact of the Kuroshio transport on SST and air-sea heat exchanges.

Data
The JMA has maintained hydrographic observations along 137°E traversing the western North Pacific from 34°N to 3°N or further south for more than 50 years (Fig. 2). The observation has been repeated every winter (usually in January) since 1967, and every summer (usually in July-August) since 1972. Details of the 137°E observation are reviewed in Oka et al. (2018). We construct a vertical-meridional gridded point values of winter and summer sections along 137°E for 0-2000 dbar by applying an optimal interpolation (OI) method for Temperature/Salinity (T/S) profiles taken by a reversing thermometer, Nansen bottles, conductivity-temperature-depth (CTD) profiler with Niskin bottles, expendable bathythermographs (XBT), expendable CTD (XCTD), and digital bathythermograph (DBT) (Table 1), as follows: (1) We produce an average cross section for 0-2000 dbar in each season, using data for 1992-2017 when observations deeper than 2000 dbar have been conducted, except for summer of 1992, 1993, 2002, and 2009 and winter of 1994 and 1995 because of shallow observations and a lack of observations (Fig. 2a); (2) we vertically interpolate T/S profiles from a reversing thermometer and Nansen bottles into 1 dbar by a method of Akima (1970), and vertically smooth those from a CTD profiler with Niskin bottles, XBT, XCTD, and DBT by Gaussian filter with an e-folding scale of 10 dbar. Note that, in summer of 2009, CTD profiles observed within ± 2° (longitude) from Argo floats, research cruise by the Fisheries Research Agency (FRA), and another cruise by the JMA are also used (Table 1) since observations between 15°N and 34°N are missing (Fig. 2a) because of severe weather; (3) we apply the OI with an exponential function with a horizontal and vertical decay scale of 160 km and 400 dbar with the average section in (1) as a first guess; and (4) finally we obtain summer and winter cross section with a horizontal interval of 1/3° for 31-34°N, 1/2° for 30-31°N, and 1° for 3-30°N (Fig. 2b) and a vertical interval of 1 dbar. For cruises whose observations are mainly in upper 1250 dbar (mainly before 1992, Fig. 2a), grid point values are calculated in upper 1250 dbar. Based on this OI gridded T/S dataset, we calculate dynamic height and zonal geostrophic current velocity relative to the reference level of 1000 dbar as a depth of no-motion. We checked the following results are insensitive to the choice of the reference level; for instance, almost identical results were obtained by using 1250 dbar as the reference level.
We also use monthly WS and sea level pressure (SLP) data from Japanese 55-year Reanalysis (JRA-55; Kobayashi et al. 2015) with a horizontal resolution of 1.25 ° (longitude) × 1.25° (latitude). WSC is calculated over the North Pacific. Since the WSC is defined by spatial derivative of WS, it tends to include small-scale features. To examine large-scale variations, we spatially smooth the WSC field by applying a 9-point average (with a weight of 1.0 for the center grid, 0.5 for side grids, and 0.3 for corner grids). The latitudinal position of the AL is defined as the latitude of a grid point with the minimum winter SLP within the region of 30-65°N, 140°E-140°W and AL intensity is defined as the minimum winter SLP within this region. The intensity of the NPSH is defined as the maximum winter SLP within 20-40°N, 120°-150°W.

Calculation of net Kuroshio transport
The Kuroshio takes various paths (e.g., Kawabe 1995), and westward flows associated with a cold-core eddy and the Kuroshio Counter Current (KCC) are located in the north and south. Therefore, a calculation of zonal transport in a fixed meridional range (e.g., regarding 30-34°N along the 137°E section as the Kuroshio) is not appropriate to estimate Kuroshio transport. We calculate the transport based on a method of Sugimoto et al. (2010); (1) we determine the northern and southern boundaries of westward flow of the cold-core eddy, Kuroshio, and KCC in each section based on a distribution of dynamic height relative to 1000 dbar (Fig. 3); (2) and we obtain cold-core eddy transport, Kuroshio transport, and KCC transport by integrating zonal geostrophic current velocity from sea surface to 1000 dbar depth within each region (Fig. 4). In this study, we especially focus on net eastward transport in the northwestern subtropical gyre and refer to it as net Kuroshio transport. We define the net Kuroshio transport as a sum of cold-core eddy transport, Table 1 Observational data used to construct a grid dataset of the 137°E hydrographic section Before 1990, temperature was mainly measured by a reversing thermometer and water samples were obtained from Nansen bottles at standard depths. The vertical interval of the measurement becomes larger with a depth from about 10 m around the sea surface to 250 m at the maximum in deep layers. Since 1990, the observation has been conducted by using a CTD profiler with Niskin bottles. The vertical resolution of the T/S profiles is 1 dbar Kuroshio transport, and KCC transport, and explore its annual-mean variation since 1972 when the summer cruise started (Fig. 5a). The net Kuroshio transport is controlled by atmospheric forcing through Sverdrup balance. Therefore, we focus on a relationship between the net Kuroshio transport variation and atmospheric variations. Each component of the net Kuroshio transport would be influenced by not only atmospheric forcing but also various oceanic processes such as thermohaline variations in surface water masses and interaction with oceanic eddies. Their variations should be investigated by using eddy-resolving ocean models in future studies.

Results
The net Kuroshio transport fluctuates on inter-annual and decadal timescales (Fig. 5). The inter-annual variation is obvious for 1972-1990 and 2000-2018. On the other hand, the decadal-scale variation is clearly detected before 2000, but absent after 2000. We investigate causes of the net Kuroshio transport variation with a focus on winter (December-February) when the atmospheric forcing is strong. A lag-correlation analysis shows that the net Kuroshio transport responds to changes of winter WSC in the central North Pacific (24-34°N, 150°E-175°W; hereafter, we refer to this region as the forcing region), characterized by the negative peak in winter WSC climatology after 2 years (Fig. 6): the WSC in the forcing region also has inter-annual and decadal timescales. The lag of 2 years is consistent with the time to reach the Kuroshio region south of Japan from the forcing region with a propagation speed of oceanic Rossby waves. We also detect positive correlation coefficients between the net Kuroshio transport and the WSC near the center of the AL (Fig. 6a). This result implies a close relationship between the net Kuroshio transport variation and AL 1 Sv = 10 6 m 3 s −1 ], b Kuroshio transport [Sv], and c KCC transport [Sv]. Positive values mean eastward transport activities as pointed out in the previous studies and discussed in the following paragraphs. We investigate the winter WSC variation in the forcing region from a view point of large-scale atmospheric circulation. We perform an Empirical Orthogonal Function (EOF) analysis for winter WSC in the subtropical North Pacific (15-40°N, 140°E-120°W) (Fig. 7a, b, d, e). The first mode is dominant along the westerlies and easterlies in the north and south of the forcing region (Fig. 7b). This mode well reflects meridional movements of the AL (Fig. 7a) (correlation coefficient r = 0.68, exceeding a 1% significance level), which show an inter-annual timescale before 1985 and after 1995, and decadal timescale before 2000 (Fig. 8b). The second mode has larger variance in the center of the North Pacific subtropical gyre (Fig. 7e), and well represents intensity fluctuations of the AL (Fig. 7d) (r = -0.82, exceeding a 1% significance level) with a decadal timescale for 1970-1990 (Fig. 8a). It is worth noting that the AL intensity is related to an occurrence of El Nino-Southern Oscillation (ENSO) events especially before 1980s and after 2000s (Fig. 7d), supporting Sugimoto and Hanawa (2009).
The net Kuroshio transport changes its predominant timescale with time (Fig. 5). It is therefore expected that the dominant forcing field on the Kuroshio would change with time. We investigate relationships between the winter WSC in the forcing region and the first two leading EOF modes by performing a running correlation analysis (Fig. 7c, f). The WSC is significantly correlated with the second mode throughout the whole period. Especially, during a period around  when the AL intensity fluctuates on a decadal timescale, the high correlation coefficients are obtained. Also, around 2000s when the AL intensity shows sharp positive and negative peaks  1975 1980 1985 1990 1995 2000 2005 2010 2015 Transport [Sv] Net Kuroshio transport  Fig. 8a, b indicates an inter-annual timescale variation of the AL is not significant after 2000 when significant inter-annual-scale signals are detected in the net Kuroshio transport. The NPSH is distributed in the region southeast of the AL, the latitudinal band of which is consistent with the Kuroshio region. The winter NPSH intensity clearly shows decadal timescale before 2000 and inter-annual timescale after 2000 (Fig. 8c). This NPSH intensity fluctuation is not linked to a meridional movement of the AL (r = 0.02). Furthermore, we found that the NPSH intensity is significantly correlated to the first EOF mode of winter WSC especially after 1970 (r = 0.57, exceeding a 1% significance level) (see Fig. 7a). Therefore, we strongly suggest that the NPSH intensity would play an important role on the net Kuroshio transport through modulation of the WSC in the forcing region by shifting westerlies meridionally.

Summary and discussion
We produced a gridded dataset of repeated hydrographic sections along 137°E based on an OI method, and then investigated net Kuroshio transport for 1972-2018. The net Kuroshio transport indicated two marked timescales with a period dependency; an inter-annual variation before 1990 and after 2000 and a decadal variation before 2000. Such dominant timescale changes in the net Kuroshio transport are never previously reported. We further showed that the net Kuroshio transport variation resulted from a combination of the meridional movements of the AL, intensity fluctuations of the AL, and intensity fluctuations of the NPSH. The AL meridional movements and NPSH intensity fluctuations (AL intensity fluctuations) are responsible for an inter-annual (decadal) net Kuroshio transport variation: such atmospheric variations are reflected to the winter WSC in the central North Pacific and result in the net Kuroshio transport variation with a lag of 2 years through a westward propagation of oceanic Rossby waves from the central North Pacific to the Kuroshio region south of Japan. Numerous studies have demonstrated impacts of the AL activities on the Kuroshio and wind-driven North Pacific subtropical gyre (e.g., Sugimoto et al. 2010), but those of the NPSH variations have been less understood yet. One important contribution of this paper is indicating a possible link between the net Kuroshio transport variation and NPSH intensity fluctuation. A relationship between NPSH variations and the Kuroshio and large-scale ocean circulation should be further investigated in future work. Since the NPSH dominates the North Pacific in a warm season, to evaluate the NPSH influence on the Kuroshio, summer atmospheric forcing may also need to be considered. As previously pointed out (e.g., Sugimoto and Hanawa 2009), the AL intensity exhibits a bi-decadal-scale variation with positive peaks around 1970,1990, and 2010 and negative peaks around 1980 and 2000 (Fig. 8a) while the bidecadal-scale signals are mainly within a cone of influence in a wavelet analysis. Interestingly, the bi-decadal feature is also detectable in the net Kuroshio transport throughout the analysis period (Fig. 5) with negative peaks around 1975, 1995, and 2015 and positive peaks around 1985 and 2005, although this is also within a cone of influence. These imply a possible link between the net Kuroshio transport and AL intensity on a bi-decadal timescale. Further accumulation of observational data will enable us to investigate the bidecadal variation of net Kuroshio transport in the future. It should be noted that changes in a Kuroshio path can affect the Kuroshio transport estimated by our method (i.e., an integration of zonal velocity): the transport can be underestimated when the Kuroshio takes a meander path and crosses the 137°E section with a slanted angle. To examine the influence of a Kuroshio path on the transport, we investigated absolute and zonal velocity for 1993-2018 based on satellite measurements (Copernicus Marine Environment Monitoring Service; marine.copernicus.eu). We found that the zonal velocity is mostly consistent with the absolute velocity ( Fig. 9) whereas their difference becomes large reaching 20% in 2017-2018 of the latest Kuroshio large meander event, implying an underestimate of the transport. However, a horizontal resolution of satellite-derived surface current data is 1/4° (longitude) × 1/4° (latitude), which is not sufficient to resolve the Kuroshio. Eddy-resolving ocean general circulation models would be useful for more accurate estimation of the Kuroshio transport and investigation of its variability.
The Kuroshio plays an important role for heat transport from low-to mid-latitude North Pacific, and would modulate SST in the western boundary region. A composite map represents that an increase of late-winter (February-March) SST in the Kuroshio/KE region associated with large net Kuroshio transport, and then, the vigorous heat is released from the ocean to atmosphere (Fig. 10a-c). These features are clearly observed in recently developed high-resolution dataset of ERA5 (Copernicus Climate Change Service 2017) ( Fig. 10d-f). These results indicate active influence of the Kuroshio on the overlying atmosphere. We suggest that the Kuroshio transport is a key factor for the establishment of air-sea interaction framework. Previously, influences of current path variations such as the Kuroshio large meander and KE dynamic state changes on SST and overlying atmosphere are well discussed (e.g., Qiu et al. 2020;Sugimoto et al. 2020Sugimoto et al. , 2021. In future work, in addition to such current path changes, impacts of the net Kuroshio transport variation should be focused.

Conflict of interest We have no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.