副热带东北太平洋混合层深度及其对潜沉的影响

The present climate simulations of the mixed layer depth(MLD) and the subduction rate in the subtropical Northeast Pacific are investigated based on nine of the CMIP5 models. Compared with the observation data,spatial patterns of the MLD and the subduction rate are well simulated in these models. The spatial pattern of the MLD is nonuniform, with a local maximum MLD(140 m) region centered at(28°N, 135°W) in late winter. The nonuniform MLD pattern causes a strong MLD front on the south of the MLD maximum region, controls the lateral induction rate pattern, and then decides the nonuniform distribution of the subduction rate. Due to the inter-regional difference of the MLD, we divide this area into two regions. The relatively uniform Ekman pumping has little effect on the nonuniform subduction spatial pattern, though it is nearly equal to the lateral induction in values. In the south region, the northward warm Ekman advection(–1.75×10~(–7) K/s) controls the ocean horizontal temperature advection(–0.85×10~(–7) K/s), and prevents the deepening of the MLD. In the ensemble mean, the contribution of the ocean advection to the MLD is about –29.0 m/month, offsetting the sea surface net heat flux contribution(33.9 m/month). While in the north region, the southward cold advection deepens the MLD(21.4 m/month) as similar as the heat flux(30.4 m/month). In conclusion, the nonuniform MLD pattern is dominated by the nonuniform ocean horizontal temperature advection. This new finding indicates that the upper ocean current play an important role in the variability of the winter MLD and the subduction rate.     

Response of the North Pacific subtropical countercurrent and its variability to global warming

Response of the North Pacific subtropical countercurrent (STCC) and its variability to global warming is examined in a state-of-the-art coupled model that is forced by increasing greenhouse gas concentrations. Compared with the present climate, the upper ocean is more stratified, and the mixed layer depth (MLD) shoals in warmer climate. The maximum change of winter MLD appears in the Kuroshio–Oyashio extension (KOE) region, where the mean MLD is the deepest in the North Pacific. This weakens the MLD front and reduces lateral induction. As a result of the reduced subduction rate and a decrease in sea surface density in KOE, mode waters form on lighter isopycnals with reduced thickness. Advected southward, the weakened mode waters decelerate the STCC. On decadal timescales, the dominant mode of sea surface height in the central subtropical gyre represents STCC variability. This STCC mode decays as CO₂ concentrations double in the twenty-first century, owing both to weakened mode waters in the mean state and to reduced variability in mode waters. The reduced mode-water variability can be traced upstream to reduced variations in winter MLD front and subduction in the KOE region where mode water forms.     

Mixed layer depth front and subduction of low potential vorticity water under seasonal forcings in an idealized OGCM

The mixed layer depth (MLD) front and subduction under seasonal variability are investigated using an idealized ocean general circulation model (OGCM) with simple seasonal forcings. A sharp MLD front develops and subduction occurs at the front from late winter to early spring. The position of the MLD front agrees with the curve where \({\rm D}T_{\rm s}/{\rm D}t = \partial T_{\rm s} /\partial t + {\user2{u}}_{\rm g} \cdot \nabla T_{\rm s} = 0\) is satisfied (t is time, \({\user2{u}}_{\rm g}\) is the upper-ocean geostrophic velocity, \(T_{\rm s}\) is the sea surface temperature (SST), and \(\nabla\) is the horizontal gradient operator), indicating that thick mixed-layer water is subducted there parallel to the SST contour. This is a generalization of the past result that the MLD front coincides with the curve \({\user2{u}}_{\rm g} \cdot \nabla T_{\rm s} = 0\) when the forcing is steady. Irreversible subduction at the MLD front is limited to about 1 month, where the beginning of the irreversible subduction period agrees with the first coincidence of the MLD front and \({\rm D}T_{\rm s}/{\rm D}t =0\) in late winter, and the end of the period roughly corresponds to the disappearance of the MLD front in early spring. Subduction volume at the MLD front during this period is similar to that during 1 year in the steady-forcing model. Since the cooling of the deep mixed-layer water occurs only in winter and SST can not fully catch up with the seasonally varying reference temperature of restoring, the cooling rate of SST is reduced and the zonal gradient of the SST in the northwestern subtropical gyre is a little altered in the seasonal-forcing case. These effects result in slightly lower densities of subducted water and the eastward shift of the MLD front.     

Relations between Sea Surface Temperature and Air-Sea Heat Flux at Periods from 1 Day to 1 Year Observed at Ocean Buoy Stations around Japan

Relations between sea-surface temperature (T _s) and heat flux at the sea surface (F) have been investigated using data from ocean observation buoys located off Shikoku in the Sea of Japan and in the East China Sea. Wavelet transformation decomposed F and T _s to wavelet coefficients (WLC) in the period-time domain. Assuming one-dimensional heat transfer by eddy diffusion in the upper ocean, the phase difference (δθ) defined as the difference between the phase of the temporal change rate of T _s, and the phase of F ranges statistically from 0 to +π/4 when F changes T _s, and is around −π/2 when heat convergence in the sea (Av) forces T _s. The δθ values are distributed from 0 to +π/4 at one-day and one-year periods at all buoys. WLC amplitude (WLCA) of F at periods from 16 to 32 day periods, which may be caused by the atmospheric ridge-trough systems, maintains energy longer than WLCA at periods from 2 to 16 days, which may be caused by monsoonal surges. At periods from 2 to 64 days, δθ values distribute from 0 to +π/4 or around −π/2 at each event, reflecting the surroundings of each ocean, i.e., Kuroshio recirculation in the off-Shikoku area, water-temperature front in the Sea of Japan, and water exchange in the continental shelf edge in the East China Sea. We demonstrate that the wavelet analysis can characterize the correspondence between irregular signals of F and T _s in various time scales and locations. This revised version was published online in August 2006 with corrections to the Cover Date.     

Water mass exchange and diapycnal mixing at Bussol’ Strait revealed by water mass properties

Intensive CTD observations that resolve the mean and tidal components were done with a total of 129 casts in summer of 2001 at Bussol’ Strait. Based on these data and all the available historical data, we have revealed the outflow from Bussol’ Strait to the Pacific and the significant diapycnal mixing in the strait. In the range 27.0−27.3σ _θ , the water property in Bussol’ Strait is almost identical to that of the Kuril Basin Water (KBW). The KBW out of Bussol’ Strait forms a water mass front with the East Kamchatka Current Water (EKCW). This front also corresponds to the front of the Oyashio Current. In the lower part of the intermediate layer (27.3−27.6σ _θ ), part of the water in the strait is characterized by lower temperature, lower salinity, and higher dissolved oxygen than that of KBW and EKCW, which can be explained only by the diapycnal mixing. The strong diapycnal mixing in the strait can also be shown by the density inversion, occurrence frequency of which corresponds well to the amplitude distribution of the diurnal current. In the density range 26.7−26.8σ _θ , the water in Bussol’ Strait has the lowest potential vorticity, suggesting that it is a source region of the low potential vorticity water. Seasonal change of the water can reach up to a density of 26.8σ _θ around Bussol’ Strait. This leads us to propose that the combination of winter convection and local tidal mixing leads to effective ventilation of the intermediate layer.     

Pumping rates of the mud shrimpCallianassa japonica

Pumping rate of a mud shrimp,Callianassa japonica, in its burrow was measured by continuous monitoring of dye concentration in the burrow water. Measurement of dilution in two directions from stained overlying seawater to normal burrow water andvice versa, gave no significant difference in results. The rate of exchange (v) of burrow water was estimated from,v=(u _t−u ₀)V/(m−u _t−1)t, whereV is volume of burrow water,u ₀,u _t−1 andu _tis dye concentration of burrow water at time 0,t−1 andt, respectively, andm is dye concentration of overlying water. The pumping rate ranged from 0.63 to 5.46 ml min⁻¹, which corresponded to a turnover time for the burrow water of 7–51 min. Short term changes in the pumping rate were correlated to intermittent behaviour of the shrimp in the burrow.     

副热带东北太平洋混合层深度及潜沉率季节循环对全球变暖的响应

本文利用第五次国际耦合模式比较计划（CMIP5）中的地球系统模式（ESM2M）,结合Argo观测数据和由Ishii等整理的再分析数据集,分析现在气候背景和辐射强迫极端增强下副热带东北太平洋海域（10°～40°N,110°～160°W）混合层深度（MLD）和潜沉率的季节变化特征,研究其对全球变暖的响应。在现在气候背景下,二者最大值均出现在冬季。潜沉率的主要贡献项存在显著的季节变化差异,1?5月主要受侧向潜沉率的变化控制,6?12月则由风应力旋度导致的埃克曼抽吸速度变化主控。全球变暖后,季节循环信号的主控要素不变。但受风应力旋度等要素变化的影响,各季节的MLD减小,大值区范围收缩。由于冬季减小幅度远大于夏季,MLD季节波动幅度（振幅）显著变小。长期看,MLD呈现持续变浅的趋势,其空间不均匀性减弱引起的MLD锋面减弱是控制侧向潜沉率减弱,最终导致总潜沉率减弱的关键。由于埃克曼抽吸速度的季节变化信号对全球变暖的响应较小,因此总潜沉率在冬季受全球变暖的影响最为强烈。上述结果表明,构成潜沉率的两个关键要素对总潜沉率的贡献比例是随着季节变化而改变的:冬季MLD锋面强盛时期,侧向潜沉率的影响将显著增强。全球变暖前后二者截然不同的变化会显著改变潜沉率的季节循环振幅,可能对该区域模态水的形成和输运产生深远的影响。     

北太平洋副热带东部模态水现在和未来的模拟分析

The present climate simulation and future projection of the Eastern Subtropical Mode Water(ESTMW) in the North Pacific are investigated based on the Geophysical Fluid Dynamics Laboratory Earth System Model(GFDL-ESM2M). Spatial patterns of the mixed layer depth(MLD) in the eastern subtropical North Pacific and the ESTMW are well simulated using this model. Compared with historical simulation, the ESTMW is produced at lighter isopycnal surfaces and its total volume is decreased in the RCP8.5 runs, because the subduction rate of the ESTMW decreases by 0.82×10-6 m/s during February–March. In addition, it is found that the lateral induction decreasing is approximately four times more than the Ekman pumping, and thus it plays a dominant role in the decreased subduction rate associated with global warming. Moreover, the MLD during February–March is banded shoaling in response to global warming, extending northeastward from the east of the Hawaii Islands(20°N, 155°W) to the west coast of North America(30°N, 125°W), with a maximum shoaling of 50 m, and then leads to the lateral induction reduction. Meanwhile, the increased northeastward surface warm current to the east of Hawaii helps strengthen of the local upper ocean stratification and induces the banded shoaling MLD under warmer climate. This new finding indicates that the ocean surface currents play an important role in the response of the MLD and the ESTMW to global warming.     

Response of the mixed layer depth and subduction rate in the subtropical Northeast Pacific to global warming

The response of the mixed layer depth(MLD) and subduction rate in the subtropical Northeast Pacific to global warming is investigated based on 9 CMIP5 models. Compared with the present climate in the 9 models, the response of the MLD in the subtropical Northeast Pacific to the increased radiation forcing is spatially nonuniform, with the maximum shoaling about 50 m in the ensemble mean result. The inter-model differences of MLD change are non-negligible, which depend on the various dominated mechanisms. On the north of the MLD front, MLD shallows largely and is influenced by Ekman pumping, heat flux, and upper-ocean cold advection changes. On the south of the MLD front, MLD changes a little in the warmer climate, which is mainly due to the upper-ocean warm advection change. As a result, the MLD front intensity weakens obviously from 0.24 m/km to0.15 m/km(about 33.9%) in the ensemble mean, not only due to the maximum of MLD shoaling but also dependent on the MLD non-uniform spatial variability. The spatially non-uniform decrease of the subduction rate is primarily dominated by the lateral induction reduction(about 85% in ensemble mean) due to the significant weakening of the MLD front. This research indicates that the ocean advection change impacts the MLD spatially non-uniform change greatly, and then plays an important role in the response of the MLD front and the subduction process to global warming.     

Formation Process of North Pacific Intermediate Water Revealed by Profiling Floats Set to Drift on 26.7 σθ Isopycnal Surface

Six newly developed floats, which were set to drift on the 26.7 σ_θ isopycnal surface and to profile temperature, salinity and pressure above 1000 dbar once a week, were deployed in the Oyashio and Kuroshio Extension (KE) in order to examine the circulation, formation site and time scale of newly formed North Pacific Intermediate Water (NPIW). The floats were deployed in February or May 2001, and the data from their deployments to December 2002 are analyzed here. Four of the six floats were deployed near the KE axis at around the first meander crest, and they moved eastward to 157°E–176°W at latitudes of 30°N–45°N. The other two floats deployed in the Oyashio water with low-potential vorticity near the south coast of Hokkaido moved southward to reach the KE front and then moved eastward to the same region as the first four floats. The temperature and salinity at 26.7 σ_θ measured by the profiling floats indicate that the source waters of NPIW, Oyashio and Kuroshio waters are drastically mixed and modified in the mixed water region west of 160°E. The floats were separated into the three paths east of 160°E between the Kuroshio Extension front and the north of Water-Mass front (nearly subarctic front). New NPIW is judged to be formed along these three paths since the vertical profiles of temperature and salinity are quite smooth, having a salinity minimum at about 26.7σ_θ along each path. Kuroshio-Oyashio isopycnal mixing ratios of the new NPIW are 7:3, 6:4 and 5:5 at 26.7σ_θ along the southern, middle and northern paths, respectively. Potential vorticity converges to about 14–15 × 10⁻¹¹ m⁻¹s⁻¹ along these paths. The time scale of new NPIW formation is estimated to be 1–1.5 years from the merger of Oyashio and Kuroshio waters to the formation of the new NPIW. This revised version was published online in July 2006 with corrections to the Cover Date.     

New type of pycnostad in the western subtropical-subarctic transition region of the North Pacific: Transition Region Mode Water

A new type of pycnostad has been identified in the western subtropical-subarctic transition region of the North Pacific, based on the intensive hydrographic survey carried out in July, 2002. The potential density, temperature and salinity of the pycnostad were found to be 26.5–26.7 σ _θ , 5°–7°C and 33.5–33.9 psu respectively. The pycnostad is denser, colder and fresher than those of the North Pacific Central Mode Water and different from those of other known mode waters in the North Pacific. The thickness of the pycnostad is comparable to that of other mode waters, spreading over an area of at least 650 × 500 km around 43°N and 160°E in the western transition region. Hence, we refer to the pycnostad as Transition Region Mode Water (TRMW). Oxygen data, geostrophic current speed and climatology of mixed layer depth in the winter suggest that the TRMW is formed regularly in the deep winter mixed layer near the region where it was observed. Analysis of surface heat flux also supports the idea and suggests that there is significant interannual variability in the property of the TRMW. The TRMW is consistently distributed between the Subarctic Boundary and the Subarctic Front. It is also characterized by a wide T-S range with similar density, which is the characteristic of such a transition region between subtropical and subarctic water masses, which forms a density-compensating temperature and salinity front. The frontal nature also tends to cause isopycnal intrusions within the pycnostad of the TRMW.     

Temporal variability of winter mixed layer in the mid-to high-latitude North Pacific

Temperature and salinity data from 2001 through 2005 from Argo profiling floats have been analyzed to examine the time evolution of the mixed layer depth (MLD) and density in the late fall to early spring in mid to high latitudes of the North Pacific. To examine MLD variations on various time scales from several days to seasonal, relatively small criteria (0.03 kg m⁻³ in density and 0.2°C in temperature) are used to determine MLD. Our analysis emphasizes that maximum MLD in some regions occurs much earlier than expected. We also observe systematic differences in timing between maximum mixed layer depth and density. Specifically, in the formation regions of the Subtropical and Central Mode Waters and in the Bering Sea, where the winter mixed layer is deep, MLD reaches its maximum in late winter (February and March), as expected. In the eastern subarctic North Pacific, however, the shallow, strong, permanent halocline prevents the mixed layer from deepening after early January, resulting in a range of timings of maximum MLD between January and April. In the southern subtropics from 20° to 30°N, where the winter mixed layer is relatively shallow, MLD reaches a maximum even earlier in December–January. In each region, MLD fluctuates on short time scales as it increases from late fall through early winter. Corresponding to this short-term variation, maximum MLD almost always occurs 0 to 100 days earlier than maximum mixed layer density in all regions.     

南海浅层经向翻转环流及相关的内部水体流动

The structure of the annual-mean shallow meridional overturning circulation(SMOC) in the South China Sea(SCS) and the related water movement are investigated,using simple ocean data assimilation(SODA) outputs.The distinct clockwise SMOC is present above 400 m in the SCS on the climatologically annual-mean scale,which consists of downwelling in the northern SCS,a southward subsurface branch supplying upwelling at around 10°N and a northward surface flow,with a strength of about 1×10~6 m~3/s.The formation mechanisms of its branches are studied separately.The zonal component of the annual-mean wind stress is predominantly westward and causes northward Ekman transport above 50 m.The annual-mean Ekman transport across 18°N is about 1.2×10~6 m~3/s.An annual-mean subduction rate is calculated by estimating the net volume flux entering the thermocline from the mixed layer in a Lagrangian framework.An annual subduction rate of about 0.66×10~6m~3/s is obtained between 17° and 20°N,of which 87% is due to vertical pumping and 13% is due to lateral induction.The subduction rate implies that the subdution contributes significantly to the downwelling branch.The pathways of traced parcels released at the base of the February mixed layer show that after subduction water moves southward to as far as 11°N within the western boundary current before returning northward.The velocity field at the base of mixed layer and a meridional velocity section in winter also confirm that the southward flow in the subsurface layer is mainly by strong western boundary currents.Significant upwelling mainly occurs off the Vietnam coast in the southern SCS.An upper bound for the annual-mean net upwelling rate between 10° and 15°N is 0.7×10~6m~3/s,of which a large portion is contributed by summer upwelling,with both the alongshore component of the southwest wind and its offshore increase causing great upwelling.     

Simulation and future projection of the mixed layer depth and subduction process in the subtropical Southeast Pacific

The present climate simulation and future projection of the mixed layer depth (MLD) and subduction process in the subtropical Southeast Pacific are investigated based on the geophysical fluid dynamics laboratory earth system model (GFDL-ESM2M). The MLD deepens from May and reaches its maximum (>160 m) near (24°S, 104°W) in September in the historical simulation. The MLD spatial pattern in September is non-uniform in the present climate, which shows three characteristics: (1) the deep MLD extends from the Southeast Pacific to the West Pacific and leads to a "deep tongue" until 135°W; (2) the northern boundary of the MLD maximum is smoothly near 18°S, and MLD shallows sharply to the northeast; (3) there is a relatively shallow MLD zone inserted into the MLD maximum eastern boundary near (26°S, 80°W) as a weak "shallow tongue". The MLD non-uniform spatial pattern generates three strong MLD fronts respectively in the three key regions, promoting the subduction rate. After global warming, the variability of MLD spatial patterns is remarkably diverse, rather than deepening consistently. In all the key regions, the MLD deepens in the south but shoals in the north, strengthing the MLD front. As a result, the subduction rate enhances in these areas. This MLD antisymmetric variability is mainly influenced by various factors, especially the potential-density horizontal advection non-uniform changes. Notice that the freshwater flux change helps to deepen the MLD uniformly in the whole basin, so it hardly works on the regional MLD variability. The study highlights that there are regional differences in the mechanisms of the MLD change, and the MLD front change caused by MLD non-uniform variability is the crucial factor in the subduction response to global warming.     

南大洋混合层的时空变化特征

过去对南大洋的研究受限于长期观测的缺乏,而现在地转海洋学实时观测阵(Arrayfor Real-timeGeostrophicOceanography,Argo)项目自开始以来持续提供了高质量的温度盐度观测,使系统地研究南大洋海洋上层结构成为可能。本研究使用2000—2018年的Argo浮标观测数据,分析了南大洋混合层深度(Mixed Layer Depth, MLD)的时空分布特征。结果表明:南大洋混合层存在明显的季节变化,冬春两季MLD在副南极锋面北侧达到最高值并呈带状分布,夏秋两季由于海表加热导致混合层变浅,季节变化幅度达到400m以上;在年际尺度上,MLD受南半球环状模(Southern HemisphereAnnularMode,SAM)调制,呈现纬向不对称空间分布特征,这与前人结果一致;本文指出在所研究时段,南大洋混合层在90°E以东,180°以西有加深趋势,而在60°W以西,180°以东有变浅趋势,显示出偶极子分布特征,并且这种趋势特征主要是风场的作用。     

Effects of mass source/sink at the western boundary on the wind-driven gyres: Implications for the ventilation of the north pacific intermediate layer through convection in the Okhotsk Sea and tidal mixing at the Kuril Straits

A steady quasi-geostrophic 2.5-layer model, forced by both Ekman pumping and a mass source/sink situated at the western boundary has been constructed to investigate the effect of diapycnal transport due to convection in the Okhotsk Sea and tidal mixing at the Kuril Straits on the intermediate layer in the North Pacific. The model illustrates a combined effect of the wind-driven and mass-driven circulations. First, net mass input induces a “barotropic” mode inter-gyre flow along the western boundary through the dynamical influence of Kelvin waves. This flow creates characteristic curves (geostrophic contours) that facilitate inter-gyre communication through the western boundary layer from the location of the mass source to the subtropical gyre. Due to the effect of wind-driven circulation, the offshore part turns eastward into the interior, encircles the outer rim of the region (which would otherwise be the pool region in the absence of mass input), and then encounters the western boundary. Eventually, the water fed into the lower layer flows mostly along this path and later flows away to the equatorial region. Conversely, in the upper layer, water is fed from the equator to the subtropics, and to the subpolar interior region through the western boundary current. The water then circulates along the outer rim and is absorbed into the mass sink. The model is controlled mainly by three nondimensional parameters: (1) the ratio of net mass input rate to the maximum Sverdrup transport (Q/T _Sv ^max ), which affects the inter-gyre communication by altering the paths of geostrophic contours, (2) the ratio of a mass input rate into the lower layer to that in total (Q ₂/Q), which controls the vertical structure of the inter-gyre flow, and (3) the measure of the wind forcing effect relative to the β effect, which determines the horizontal extent of the area influenced by the mass input. The other parameter regimes with respect to Q/T _Sv ^max and Q ₂/Q are also presented.     

Determination of the surface temperature and the temperature drop in the skin layer from a moving vessel using IR data

In this paper methods are given of the IR radiometric determination of the temperature of the radiating oceanic layer,T _s, and the determination on its basis of the temperature drop in the skin layer,dt, from on board a ship underway, which allows for the emissivity of the sea surface and the contribution of the atmosphere to the upward radiation of the sea surface. The seawater temperatureT _w was registered by a contact sensor towed along the ship at a depth of 0.1–0.4 m. The valuedt=T _s–T _w was found to be within the limits of 0.07 to –0.85 K, being on average 0.33 K. The near-surface wind speedv (v=1–13 m s^–1) had the best relationship withdt;dt=–0.66+0.06v. The empirical values ofdt were compared with the model ones obtained by the well-known formulae for the conditions of forced and free convection of the surface layer.Translated by Mikhail M. Trufanov.     

Circulation and mixing at the New England shelfbreak front: Results of purposeful tracer experiments

We present the results of six dye tracer experiments that measured the mixing and circulation at the shelfbreak front on the New England Shelf. The last three were conducted during the New England Shelfbreak Productivity Experiment (NESPEX) with concurrent isopycnal float deployments. The results are consistent with the Chapman and Lentz [Chapman, D.C., and Lentz, S.J. (1994). Trapping of a coastal density front by the bottom boundary layer. Journal of Physical Oceanography, 24, 1465–1479.] model prediction of the separation and upwelling along the shelfbreak front of bottom boundary layer (BBL) water forced by an Ekman buoyancy flux, but show considerable variability. Cross-shelf velocities at the detachment point are 2–3 × 10⁻² m/s. But seaward, over the slope region, dye tagged water was sheared from the main patch into small filaments that upwelled along the front with cross-shelf speeds up to 0.1 m/s. Cross-shelf diffusion was of order 10 m²/s in the mixed bottom layer and 1 m²/s in the interior along the front. Within the stratified front, the mean vertical diffusivity was K_z  4 × 10⁻⁶ m²/s. The dispersion of shelfwater in the slope region is effected by turbulent flow with advective speeds exceeding the small scale diffusive mixing. The mean flux of the detached BBL water is sufficient to account for the net loss of shelf water during its transit from Cape Cod to Cape Hatteras.     

Convergence,divergence and vertical velocity at a tidal front in Hiuchi-Nada,Japan

The detailed flow structure around a tical front in Hiuchi-Nada, Japan was observed with the use of ADCP (Acoustic Doppler Current Profiler). The surface convergence region is observed at the transition zone between vertically well mixed area and the stratified area. The surface divergence regions exist next to the surface convergence region. The strong downward current is estimated in the middle layer just below the surface convergence region. The maximum surface convergence and the maximum downward velocity in the middle layer are 1.0×10^–4 s^–1 and 0.12 cm s^–1, respectively.     

Ventilation of the North Pacific subtropical pycnocline and mode water formation

The annual subduction rate of the North Pacific was calculated based on isopycnally averaged hydrographic climatology (HydroBase), high-resolution winter mixed-layer climatology (NWMLC), and various wind stress climatologies from ship reports, numerical weather prediction products, and satellite products. The calculation was performed using Lagrangian coordinates in the same manner as in previous works, except a less smoothed oceanic climatology (HydroBase and NWMLC) was used instead of a World Ocean Atlas. Differences in the wind stress climatologies have very little effect on subduction rate estimates. The subduction rate census for density classes showed peaks corresponding to subtropical mode water (STMW), central mode water (CMW), and eastern subtropical mode water (ESTMW). The deeper mixed layer and the associated sharper mixed-layer fronts in the present climatology resulted in a larger lateral induction, which boosted the subduction rate, especially for the potential density anomaly (σ_θ) range of the lighter STMW (25.0 < σ_θ < 25.2 kg m⁻³) and lighter CMW (26.0 < σ_θ < 26.2 kg m⁻³), compared to previous estimates. The renewal time of permanent pycnocline water was estimated as the volume of water divided by the subduction rate for each σ_θ class: 2–4 years for ESTMW (24.5 < σ_θ < 25.2 kg m⁻³), 2 years for the lighter STMW (25.0 < σ_θ < 25.3 kg m⁻³), 5–9 years for the denser STMW (25.3 < σ_θ < 25.6 kg m⁻³), 10–20 years for the lighter CMW (26.0 < σ_θ < 26.2 kg m⁻³), 20–30 years for the middle CMW (26.2 < σ_θ < 26.3 kg m⁻³), and 60 years or longer for the denser CMW (26.3 < σ_θ < 26.6 kg m⁻³). A comparison of the water volume and subduction rate in potential temperature–salinity (θ–S) space indicated that the upper permanent pycnocline water (25.0 < σ_θ < 26.2 kg m⁻³) was directly maintained by nondiffusive subduction of winter surface water, including STMW and lighter CMW. The lower permanent pycnocline water (26.2 < σ_θ < 26.6 kg m⁻³) may be maintained through the subduction of fresher and colder water from the subarctic–subtropical transition region and subsequent mixing with saltier and warmer water. Diagnosis of the potential vorticity (PV) of the subducted water demonstrated that the low PV of STMW was mainly due to the large subduction rate, whereas that of both ESTMW and CMW was due mainly to the small density advection rate (cross-isopycnal flow). Additionally, a relatively large subduction rate probably contributes to the low PV of part of the lighter CMW (ESTMW) formed in the region around 38°N and 170°W (28°N and 145°W), which is characterized by a relatively thick winter mixed layer and an associated mixed-layer front, causing a large lateral induction rate.