Abstract
Accurately capturing the Kuroshio’s main axis and its multiscale frontal variations remains challenging due to the constraints of traditional fixed-section extraction methods. Here, we develop an adaptive iterative tracking algorithm utilizing high-resolution reanalysis data (2002–2024) that dynamically adjusts search directions and cross-sections via local velocity vectors, integrated with a dynamic step size and two-dimensional validation. Applying a multiscale variability decomposition framework across four key regions reveals distinct spatiotemporal dynamics. The North Equatorial Current bifurcation zone exhibits a significant strengthening trend driven by seasonal zonal and decadal meridional flows. Conversely, the Kuroshio east of Taiwan is dominated by high-frequency mesoscale processes (~70%) with a semi-annual cycle and no long-term trend. The East China Sea front maintains a highly stable seasonal meridional signal (25%). Crucially, the Luzon Strait intrusion shows a significant long-term weakening trend (~0.0029 m·s−1·a−1, p < 0.01), characterized by eastward strengthening and northward weakening, with ENSO significantly modulating its seasonal cycle. This approach substantially reduces systematic extraction errors compared to traditional fixed-section methods, as independently verified using satellite SST frontal gradients (median deviation < 0.2°), providing critical observational evidence for understanding western boundary current–marginal sea interactions and their dynamical responses under global warming.