Numerical Study of a Kyucho and a Bottom Intrusion in the Bungo Channel,Japan: Disturbances Generated by the Kuroshio Small Meanders

The generation and propagation mechanisms of a Kyucho and a bottom intrusion in the Bungo Channel, Japan, have been studied numerically using the hydrostatic primitive equations by assuming density stratification during summer. The experiments are designed to generate a Kuroshio small meander in Hyuga-Nada, which acts as a trigger for these disturbances. After the current speed of the Kuroshio is changed, a small meander is generated. At the head of the small meander, warm Kuroshio water is engulfed, and encounters the southwest coast of Shikoku. However, convergence of heat flux on the bump off Cape Ashizuri suppresses the generation of a warm disturbance, if the current speed is large. As the cold eddy associated with the small meander approaches Cape Ashizuri, the heat flux diverges on the bump. This heat source forces a warm disturbance, which intrudes along the east coast of the Bungo Channel as a baroclinic Kelvin wave (a Kyucho). After the cold eddy passes off Cape Ashizuri, the Kuroshio approaches the bump again. Strong convergence of heat flux then occurs on the bump, which forces a cold disturbance. This disturbance propagates as a topographic Rossby wave along the shelf break at the mouth of the channel. After the topographic wave reaches the west end of the shelf break, it intrudes along the bottom layer of the channel as a density current (a bottom intrusion). These results suggest that a Kyucho and a bottom intrusion are successive events associated with the propagation of the small meander.     

Upper Ocean Sensitivity to Wind Forcing in the South China Sea

The Naval Research Laboratory (NRL) Layered Ocean Model (NLOM) has been used to investigate the sensitivity of the upper South China Sea (SCS) circulation to various atmospheric wind forcing products. A 1/16° 6-layer, thermodynamic Pacific Ocean north of 20°S version of NLOM has been integrated using observed climatological monthly mean winds (Hellerman and Rosenstein, 1983) and climatologies based on two atmospheric prediction models: the European Centre for Medium-Range Weather Forecasts (ECMWF) and the National Centers for Environmental Prediction (NCEP). ECMWF products include the 10 meter winds (at both 1.125° and 2.5° resolution) and surface stresses (1.125°). The NCEP forcing (1.875°) is a surface stress product. Significant differences exist in the wind stress curl patterns and this is reflected in the upper ocean model response, which is compared to observational data. The model experiments suggest the generation of the West Luzon Eddy is controlled by positive wind stress curl. The degree of Kuroshio intrusion into the SCS, however, is not affected by wind stress curl but is governed by the coastline geometry of the island chain within Luzon Strait. The summertime offshore flow from the Vietnamese coast is present in all simulations but the dipole structure on either side of the jet is variable, even among experiments with similar wind stress curl patterns. The ECMWF surface stresses exhibit spurious coastal wind stress curl patterns, especially in locations with significant orographic features. This manifests itself in unrealistic small scale coastal gyres in NLOM. High resolution basin-scale and coastal models might be adversely affected by these stresses. This revised version was published online in July 2006 with corrections to the Cover Date.     

Correspondence between Lagrangian and Eulerian Velocity Statistics at the ASUKA Line

The 20–30% differences between drifter- and altimetry-derived mean velocities reported at the ASUKA line (Affiliated Surveys of the Kuroshio off Cape Ashizuri) are explained in terms of the non-Eulerian character of statistics of drifter trajectories' crossovers with the TOPEX/Poseidon track. Larger mean drifter velocities are shown to result from the fact that more water crosses the line when velocity is large, so that the “line” average velocity corresponds to the Eulerian quantity <V⋅|V_n|>/<|V_n|>, where V_n is a velocity component perpendicular to the line. Practical methods are suggested to correct the bias. General properties of the “line” statistic and improvement of dynamical balance by correct accounting for the acceleration term are discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Predictability of Interannual Variability in the Kuroshio Transport South of Japan Based on Wind Stress Data over the North Pacific

It is expected that a roughly two-year forecast of the Kuroshio transport variation can be made from a past record of wind stress data over the ocean, since it takes nearly ten years for the first-mode baroclinic Rossby wave to traverse the entire basin in the midlatitude North Pacific (∼30°N). We therefore investigated the predictability using an ocean general circulation model driven by the wind stress data from the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis. Referring to a hindcast experiment as the control run, we carried out fifteen forecast experiments, the initial conditions of which are taken from the hindcast experiment at intervals of two years during the period from the end of 1969 to the end of 1997. Each of the forecast experiments is driven only by wind stress in the year preceding each experiment. The forecasted Kuroshio transport anomaly south of Japan agrees better with the hindcasted one during the first two years of the forecast in most cases. In some cases, however, significant disagreements occur, most of which are likely due to larger unpredictable variations caused by wind stress anomalies near Japan. At the end of forecast year 2, the anomaly correlation coefficient is about 0.7, and rms of the forecast error is smaller than rms of the hindcasted anomaly. These results indicate that the prediction of the interannual variability in the Kuroshio transport could be made two years in advance at a statistically significant level. This revised version was published online in July 2006 with corrections to the Cover Date.     

2002年4～5月琉球群岛两侧海流的研究

基于日本气象厅“长风丸”调查船在2002年4～5月航次期间的CTD资料,结合卫星风场资料,采用改进逆方法计算了琉球群岛两侧海域各断面的流速和流量分布,并分析卫星跟踪浮标资料和同期的卫星高度计资料,得出下面一些主要结论:(1)黑潮流速在PN断面上只有一个流核.通过断面PN的净东北向流量约为34.7×106m3/s,此流量包括台湾暖流、东海黑潮和黑潮以东的反气旋涡的流量.(2)黑潮流速在断面TK上有两个流核,通过断面TK净东向的流量为25.6×106m3/s,黑潮通过海峡后流向断面ASUKA.(3)冲绳岛东南海区琉球海流的流量约为8.8×106m3/s,并流向断面AM.(4)奄美大岛以东的北向海流的流量为12.7×106m3/s,并流向断面ASUKA.在断面ASUKA东南部出现一个中尺度反气旋涡,直径约240 km,其流量约为28.5×106m3/s.(5)四国以南黑潮第一层水体基本来源于通过吐噶喇海峡的黑潮,第二、三层水体来自吐噶喇海峡和奄美大岛以东海域的流量大致相当,而第四层的流量则主要来自于奄美大岛以东海域.(6)浮标资料显示,奄美大岛以东的海流部分来自于断面AM以东海区,并通过断面ASUKA.     

POM and Two-Way Nesting POM Study of Kuroshio Damping Phenomenon Caused by a Strong Wind

A Kuroshio damping phenomenon of a few days scale caused by a strong wind was studied using the Princeton Ocean Model (POM) and a two-way nesting model for the POM. We simulated an idealized eastward Kuroshio in a zonal channel in terms of the inflow/outflow condition at the western/eastern boundaries and examined responses of the idealized Kuroshio to a strong easterly wind. This study was motivated by observations of JMA research vessel Shumpu Maru which reported deformation of isopycnals and damping of the Kuroshio before and after the attack of Typhoon 6804. Deformation of isopycnals and damping of the Kuroshio are found to depend on magnitude of wind stress, intensity of stratification, distance of the Kuroshio from coast, and grid resolution. The deformation of isopycnals is related to the vertical circulation pattern caused by the divergence/convergence of the onshore Ekman transport. A simple stratification model composed of Niiler (1969)'s modification of the Ekman transport and the coastal boundary is proposed to explain the damping phenomenon. An idealized cape was added in other experiments in order to study whether the Kuroshio damping mechanism discussed here works in three dimensions. The newly developed two-way nesting model for the POM was applied for this experiment and made it possible to clarify more detailed features of response in the nested area than did the coarse grid model.     

Spatio-temporal variability in sea surface wind stress near and off the east coast of Korea

Sea surface wind stress variabilities near and off the east coast of Korea, are examined using 7 kinds of wind datasets from measurements at 2 coastal (land) stations and 2 ocean buoys,satellite scatterometer (QuikSCAT), and global reanalyzed products (ECMWF,NOGAPS,and NCEP/NCAR). Temporal variabilities are analyzed at 3 frequency bands; synoptic (2-20 d), intra-seasonal (20-90 d),and seasonal (>90 d).Synoptic and intra-seasonal     

Inter-Comparison of Numerical Model Generated Surface Winds with QuikSCAT Winds over the Indian Ocean

The accurate surface wind in the equatorial Indian Ocean is crucial for modeling ocean circulation over this region. In this study, the surface wind analysis generated at the European Center for Medium Range Weather Forecasts (ECMWF) and the National Centers for Environmental Prediction (NCEP) are compared with NASA QuikSCAT satellite derived Level2B (swath level) and Level3 (gridded) surface winds for the year 2005. It is observed that the ECMWF winds exhibit speed bias of 1.5 m/s with respect to QuikSCAT Level3 in the southern equatorial Indian Ocean. The NCEP winds are found to exhibit speed bias (1.0–1.5 m/s) in the southern equatorial Indian Ocean specifically during January–February 2005. The biases are also observed in the analysis when compared with Level2B product as well; however, it is less in comparison to Level3 products. The amplitude of daily variations of both ECMWF and NCEP wind speed in Bay of Bengal and parts of the Arabian Sea is about 80% of that in QuikSCAT, while in the equatorial Indian Ocean it is about 60% of that of QuikSCAT.     

Evaluation of relative performance of QuikSCAT and NCEP re-analysis winds through simulations by an OGCM

The response of an ocean general circulation model (OGCM) to two different wind products, viz., NCEP/NCAR reanalysis and QuikSCAT scatterometer, was examined. OGCM-simulated thermodynamic variables from the two simulations, hereafter referred to as NCEP-R (NCEP/NCAR wind forced) and QS-R (QuikSCAT wind forced) were intercompared and also were compared against observations for a period of 3 years (2000–2002). In the tropical Indian Ocean (IO), the sea-level anomaly (SLA) simulated by QS-R has less root mean square error (RMSE) and higher correlation with respect to TOPEX/Poseidon SLA observations than SLA simulated by NCEP-R. Intraseasonal variability of currents observed by TRITON buoy in the IO was closely captured by QS-R, although the magnitudes are somewhat underestimated. Surface currents simulated by QS-R have less RMSE than those simulated by NCEP-R in the Pacific. However, the sub-surface currents are much weaker in magnitude in both the solutions, possibly because of deficiencies in the diffusion and viscosity parameterization. Sea-surface temperature (SST) simulated by QS-R has a cooler bias. The RMSE of SST simulated by NCEP-R is less than the RMSE of SST simulated by QS-R, with the latter capturing the variabilities more realistically. The large differences between SST simulated by QS-R and observations could be partly due to physical inconsistency between the momentum and heat fluxes. Scatterometer-forced model simulations of 20^oC thermocline depths (D20) are in better agreement with in situ-derived D20 than the D20 simulated by NCEP-R. Variations in the mixed layer depth at the TRITON buoy are better captured by QS-R than by NCEP-R. Speed of Kelvin and Rossby waves and the strength of upwelling/downwelling features in the IO are closer to observations in QS-R than in NCEP-R simulations.     

Comparison of surface wind stress characteristics over the Tropical Atlantic (10°N–40°S) in fields derived from the UWM/COADS, NCEP/NCAR and QuikSCAT datasets

A comparison of monthly wind stress derived from winds of NCEP/NCAR (National Centers for Environmental Prediction/National Center for Atmospheric Research) reanalysis and UWM/COADS (The University of Wisconsin-Milwaukee/Comprehensive Ocean-Atmosphere Data Set) dataset (1950–1993), and of NCEP/NCAR reanalysis and satellite-based QuikSCAT dataset (2000–2006), is made over the South Atlantic (10°N–40°S). On a mean seasonal scale, the comparison shows that these three wind stress datasets have qualitatively similar patterns. Quantitatively, in general, from about the equator to 20°S in the mid-Atlantic the wind stress values are stronger in NCEP/NCAR data than those in UWM/COADS data. On the other hand, in the Intertropical Convergence Zone (ITCZ) area the wind stress values in NCEP/NCAR data are slightly weaker than those in UWM/COADS data. In the South Atlantic, between 20° S–40°S, the QuikSCAT dataset presents complex circulation structures which are not present in NCEP/NCAR and UWM/COADS data. The wind stress is used in a numerical ocean model to simulate ocean currents, which are compared to a drifting-buoy observed climatology. The modeled South Equatorial Current agrees better with observations between March–May and June–August. Between December–February, the South Equatorial Current from UWM/COADS and QuikSCAT experiments is stronger and more developed than that from NCEP/NCAR experiment. The Brazil Current, in turn, is better represented in the QuikSCAT experiment. Comparison of the annual migration of ITCZ at 20° and 30°W in UWM/COADS and NCEP/NCAR data sources show that the southernmost position of ITCZ at 30°W in February, March and April coincides with the rainy season in NE Brazil, while the northernmost position of ITCZ at 20°W in August coincides with the maximum rainfall of Northwest Africa.     

Quasiperiodic small meanders of the Kuroshio off Cape Ashizuri and their inter-annual modulation caused by quasiperiodic arrivals of mesoscale eddies

A quasiperiodic variation of 100–110 days in the Kuroshio path off Cape Ashizuri, resulting from the passage of small meanders, was detected by observation with moored current meters during 1993–1995. TOPEX/POSEIDON altimeter data covering 9 years showed that the quasiperiodic variation period was not persistent and modulated twice, with a ∼110-day period from mid-1993 to late 1996, a ∼150-day period from late 1996 to mid-1999, and a ∼110-day period from mid-1999 to late 2001. The quasiperiodic variations of the Kuroshio path migration were contemporaneous with the quasiperiodic arrivals of mesoscale eddies from the east along 27–32°N over the same ∼110- and ∼150-day period quasiperiodic variations. The periodic arrivals of the eddies configure the periodic variations of the Kuroshio path and its inter-annual modulation.     

基于浮标、NCEP2及OAFlux产品的新海表感热和潜热通量对比分析

New satellite-derived latent and sensible heat fluxes are performed by using Wind Sat wind speed, Wind Sat sea surface temperature, the European Centre for Medium-range Weather Forecasting(ECMWF) air humidity, and ECMWF air temperature from 2004 to 2014. The 55 moored buoys are used to validate them by using the 30 min and 25 km collocation window. Furthermore, the objectively analyzed air-sea heat fluxes(OAFlux) products and the National Centers for Environmental Prediction-National Center for Atmospheric Research reanalysis 2(NCEP2) products are also used for global comparisons. The mean biases of sensible and latent heat fluxes between Wind Sat flux results and buoy flux data are –0.39 and –8.09 W/m~2, respectively. In addition, the rootmean-square(RMS) errors of the sensible and latent heat fluxes between them are 5.53 and 24.69 W/m~2,respectively. The RMS errors of sensible and latent heat fluxes are observed to gradually increase with an increasing buoy wind speed. The difference shows different characteristics with an increasing sea surface temperature, air humidity, and air temperature. The zonal average latent fluxes have some high regions which are mainly located in the trade wind zones where strong winds carry dry air in January, and the maximum value centers are found in the eastern waters of Japan and on the US east coast. Overall, the seasonal variability is pronounced in the Indian Ocean, the Pacific Ocean, and the Atlantic Ocean. The three sensible and latent heat fluxes have similar latitudinal dependencies; however, some differences are found in some local regions.     

Variation of the Kuroshio in the Tokara Strait Induced by Meso-Scale Eddies

Temporal variations of the Kuroshio volume transport in the Tokara Strait and at the ASUKA line are decomposed by phase-propagating Complex EOF modes of high-resolution sea surface dynamic topography (SSDT) field during the first tandem period of TOPEX/POSEIDON and ERS-1 (from October 1992 to December 1993). Both variations are dominated by a mode with nearly semi-annual cycle, which indicates a series of interactions between the Kuroshio and meso-scale eddies. Namely, northern part of a westward-propagating meso-scale eddy at 23°N is captured into the southern side of the Kuroshio at the south of Okinawa, then it moves downstream along the Kuroshio path passing the Tokara Strait, and reaches to the ASUKA line where it merges with another eddy propagating from the east at 30°N. The variation at the ASUKA line is, however, less dominated by this mode; instead, it includes the SSDT variations in the south of Shikoku and the east of Kyushu which would be directly affected by eddies from the east without passing the Tokara Strait. On the other hand, the same analysis for movements of the Kuroshio axis in the Strait indicates that they are governed by short-term variations locally confined to the Kuroshio in the East China Sea without being induced by meso-scale eddies. This results, however, seem to depend strongly on a time scale of interest. It is suggested that the long-term movements of the Kuroshio axis in the Strait would demonstrate coincidence with SSDT variation in the south of Japan.     

Intraannual variability of sea level around Tosa Bay

Through analysis of monthly in situ hydrographic, tide gauge, altimetry and Kuroshio axis data for the years 1993–2001, the intraannual variability of sea level around Tosa Bay, Japan, with periods of 2–12 months is examined together with the intraannual variability of the Kuroshio south of the bay. It is shown that the intraannual variation of steric height on the slope in Tosa Bay can account for that of sea level at the coast around the bay as well as on this slope. It is found that the steric height (or sea level) variation on the slope in this bay is mainly controlled by the subsurface thermal variation correlated with the Kuroshio variation off Cape Ashizuri, the western edge of Tosa Bay. That is, when the nearshore Kuroshio velocity south of the cape is intensified [weakened] concurrently with the northward [southward] displacement of the current axis, temperature in an entire water column decreases [increases] simultaneously, mainly due to the upward [downward] displacement of isotherms, coincident with that of the main thermocline. It follows that the steric height (or sea level) decreases [increases].     

Sea level response to wind field fluctuation around the tip of the Izu Peninsula

The mechanism of a characteristic sea level response (barotropic coastal ocean response) to wind field fluctuation around the tip of the Izu Peninsula observed during the middle of December 2000 to the middle of January 2001 was investigated based on three types of numerical experiments using the Princeton Ocean Model with various parameters. The response was characterized by the relaxation of sea level falling (rising) during eastward upwelling (westward downwelling) favorable wind regime. Analyses of quasi-realistic numerical model results in terms of the vertically integrated momentum balances and vorticity balance for the barotropic mode revealed that: 1) development/abatement of two anomalous circulations generated around the tip of the Izu Peninsula controls the sea level response through the acceleration/deceleration of a quasi-geostrophic barotropic coastal current between the circulations; 2) nonlinear vorticity advection by the Kuroshio Current and by the coastal current, coupled with vorticity diffusion, decelerates the quasi-geostrophic coastal current in the latter half of the wind regimes, which induces the relaxation of sea level rise/fall. The results of the quasi-realistic numerical experiment suggest that an analysis of the vorticity balance for the barotropic mode contributes to a better understanding of sea level responses to wind in coastal regions with strong currents and complex topography. In addition, a numerical experiment with idealized spatially uniform density stratification and a quasi-realistic wind field shows that if the Kuroshio Current had been shifted far offshore from the Izu Peninsula during the observation period, westward propagating continental shelf waves would have controlled the coastal sea level response.     

Estimation of mean and extreme waves in the East China Seas

Accurately estimating the mean and extreme wave statistics and better understanding their directional and seasonal variations are of great importance in the planning and designing of ocean and coastal engineering works. Due to the lack of long-term wave measurement data, the analysis of extreme waves is often based on the numerical wave hind-casting results. In this study, the wave climate in the East China Seas (including the Bohai Sea, the Yellow Sea and the East China Sea) for the past 35 years (1979–2013) is hind-casted using a third generation wave model – WAMC4 (Cycle 4 version of WAM model). Two sets of reanalysis wind data from NCEP (National Centers for Environmental Prediction, USA) and ECMWF (European Centre for Medium-range Weather Forecasts) are used to drive the wave model to generate the long-term wave climate. The hind-casted waves are then analysed to study the mean and extreme wave statistics in the study area. The results show that the mean wave heights decrease from south to north and from sea to land in general. The extreme wave heights with return periods of 50 and 100 years in the summer and autumn seasons are significantly higher than those in the other two seasons, mainly due to the effect of typhoon events. The mean wave heights in the winter season have the highest values, mainly due to the effect of winter monsoon winds. The comparison of extreme wave statistics from both wind fields with the field measurements at several nearshore wave observation stations shows that the extreme waves generated by the ECMWF winds are better than those generated by the NCEP winds. The comparison also shows the extreme waves in deep waters are better reproduced than those in shallow waters, which is partly attributed to the limitations of the wave model used. The results presented in this paper provide useful insight into the wave climate in the area of the East China Seas, as well as the effect of wind data resolution on the simulation of long-term waves.     

太平洋海浪场时空特征分析

为了深入了解太平洋海浪场特别是涌浪场的时空分布特征,利用44 a(1958~2001年)ERA-40海浪再分析资料对南、北太平洋风浪和涌浪的波高和波向进行了统计分析,结果表明:北太平洋海浪场比南太平洋具有更明显的季节变化特征,四季中南太平洋涌浪均有明显的越赤道北传过程.南、北半球西风带海浪波高随时间呈线性增长趋势,且涌浪分别存在2.4~3.7 a 和2.9 a 左右的显著周期,风浪和混合浪波高存在6.5 a 和5.2 a 的共同周期     

Climatic variability of transport in the upper layer of the Antarctic Circumpolar Current from hydrological and satellite data

Based on the satellite altimetry dataset of sea level anomalies, the climatic hydrological database World Ocean Atlas-2009, ocean reanalysis ECMWF ORA-S3, and wind velocity components from NCEP/NCAR reanalysis, the interannual variability of Antarctic Circumpolar Current (ACC) transport in the ocean upper layer is investigated for the period 1959–2008, and estimations of correlative connections between ACC transport and wind velocity components are performed. It has been revealed that the maximum (by absolute value) linear trends of ACC transport over the last 50 years are observed in the date-line region, in the Western and Eastern Atlantic and the western part of the Indian Ocean. The greatest increase in wind velocity for this period for the zonal component is observed in Drake Passage, at Greenwich meridian, in the Indian Ocean near 90° E, and in the date-line region; for the meridional component, it is in the Western and Eastern Pacific, in Drake Passage, and to the south of Africa. It has been shown that the basic energy-carrying frequencies of interannual variability of ACC transport and wind velocity components, as well as their correlative connections, correspond to the periods of basic large-scale modes of atmospheric circulation: multidecadal and interdecadal oscillations, Antarctic Circumpolar Wave, Southern Annual Mode, and Southern Oscillation. A significant influence of the wind field on the interannual variability of ACC transport is observed in the Western Pacific (140° E–160° W) and Eastern Pacific; Drake Passage and Western Atlantic (90°–30° W); in the Eastern Atlantic and Western Indian Ocean (10°–70° E). It has been shown in the Pacific Ocean that the ACC transport responds to changes of the meridional wind more promptly than to changes of the zonal wind.     

Long-term variability of oceanic frontal zones associated with large-scale atmospheric forcing

Global satellite sea surface temperature (SST) measurements and NCEP/NCAR reanalysis wind data for the period of 1982–2009 have been used to study the relationship between long-term variability of oceanic frontal zones (OFZ) and large-scale atmospheric forcing. Statistically significant positive correlations between the maximum magnitude of the meridional gradient of zonally averaged SST and meridional shear of zonal wind (which is an estimate of the Ekman convergence intensity) were found for all subpolar and subtropical OFZ of the World Ocean. Variability of the latitudinal position of OFZ cores may be associated with Ekman advection variability due to zonal wind variations. Strengthening of zonal wind results in a shift of subpolar OFZ cores to the south/north in the Northern/Southern hemispheres.     

Verification of the wind-driven transport in the North Pacific subtropical gyre using gridded wind-stress products

Wind-stress products supplied by satellite scatterometers carried the European Remote-sensing Satellite (ERS) and QuikSCAT (QSCAT), together with numerical weather predictions from the European Centre for Medium Range Weather Forecasting (ECMWF) and the National Centre for Environmental Prediction (NCEP) were used to estimate wind-driven transports of the North Pacific subtropical gyre. At 30°N, we compared the wind-driven transports with geostrophic transports calculated from World Ocean Database 2005. The wind-driven transports for QSCAT and NCEP are in good agreement with the geostrophic transport within reasonable error, except for a regional difference in the eastern part of the section. The difference in the eastern part suggests an anti-cyclonic deviation of the geostrophic transport, resulting from an anti-cyclonic anomalous flow in the surface layer. It is suggested that this anomalous flow is the Eastern Gyral, produced by the thermohaline process associated with the formation of the Eastern Subtropical Mode Water. To investigate the validity of QSCAT and NCEP data, we examined whether or not the Sverdrup transports for these products are consistent with the transport of the western boundary current estimated by past studies. The net southward transport, given by the sum of the Sverdrup transport for QSCAT and NCEP and the thermohaline transport, agrees well with the net northward transport of the western boundary current. From this result, together with the fact that the wind-driven transports for these products are in good agreement with the geostrophic transport, we conclude that the Sverdrup balance can hold in the North Pacific subtropical gyre.