Observations of turbulence under weakly and highly stratified conditions in the Ariake Sea

Observations of turbulence, stratification, and mean current were made using a microstructure profiler and an acoustic Doppler current profiler (ADCP) during four cruises at a central location in the Ariake Sea, under weakly and strongly stratified conditions. Continuous measurements of the dissipation rate of turbulent kinetic energy (TKE), ε, were made. These revealed that frictional bed turbulence with quarterdiurnal variation in the bottom boundary layer (BBL) was one of the most energetic sources of vertical mixing in the sea. Thickness of the BBL was strongly confined by the stable stratification. We investigate a relationship between the BBL height h and the Ozmidov scale. We present a systematic argument that describes the vertical structure and characteristic scales of velocity and turbulence inside the frictional BBL, where the stratification persisted. Considerable deviation of observed vertical shear from the law of the wall indicated a modification of turbulent scales by the stratification. Shear stress calculated from the velocity data using vertical integration of the equation of motion was found to decrease approximately linearly with height. The TKE production rate P, estimated using the shear stress, was highly correlated with the dissipation rate. The buoyancy contribution to TKE balance in the BBL was quantified in terms of the flux Richardson number R _f as R _f?=?0.12.     

夏季长江口及其邻近海域湍流特征分析

利用2019年7月在长江口科学考察实验研究夏季航段(NORC2019-03-02)中获得的MSS90L湍流剖面仪的直接观测数据,本文计算并分析了该断面的湍动能耗散率ε和垂向湍扩散系数K_Z的分布情况。湍动能耗散率的大小为1.72×10?10～2.95×10?5 W/kg;垂向湍扩散系数的大小为3.24×10?7～4.55×10?2 m2/s。湍动能耗散率和垂向湍扩散系数的分布相似,均为上层最强,底层次之,中层最弱。上层由于风应力的作用,使得湍动能耗散率和垂向湍扩散系数较大;温跃层处层化较强,抑制了湍动能的耗散和垂向上的湍混合。盐度锋面的次级环流会促使低盐水团脱离,锋面引起的垂向环流会加强海洋的湍混合。低盐水团与外界的能量交换较少,湍动能耗散率较弱。长江口海区存在明显的上升流和下降流,它们是由锋面的次级环流产生的;上升流和下降流的存在促进湍动能的耗散与湍混合。     

On the turbulent prandtl number in a stably stratified atmospheric boundary layer

The results obtained from both atmospheric and laboratory measurements and from LES data show that, in the stably stratified flows of the atmospheric boundary layer, turbulent mixing occurs at gradient Richardson numbers Ri _g that significantly exceed one: the inverse turbulent Prandtl number Pr _t ⁻¹ decreases with an increase in the thermal flow stability. The decreasing trend of the inverse turbulent Ptandtl number is reproduced in a stably stratified atmospheric boundary layer in agreement with measurement data with the aid of an improved three-parameter turbulence model. In this model, a modified model that takes into account the effect of stratification in the expression for the time scale of the scalar field is used for the pressure-scalar correlation.     

Dynamics of surf-zone turbulence in a strong plunging breaker

The characteristics of turbulence created by a plunging breaker on a 1 on 35 plane slope have been studied experimentally in a two-dimensional wave tank. The experiments involved detailed measurements of fluid velocities below trough level and water surface elevations in the surf zone using a fibre-optic laser-Doppler anemometer and a capacitance wave gage. The dynamical role of turbulence is examined making use of the transport equation for turbulent kinetic energy (the k-equation). The results show that turbulence under a plunging breaker is dominated by large-scale motions and has certain unique features that are associated with its wave condition. It was found that the nature of turbulence transport in the inner surf zone depends on a particular wave condition and it is not similar for different types of breakers. Turbulent kinetic energy is transported landward under a plunging breaker and dissipated within one wave cycle. This is different from spilling breakers where turbulent kinetic energy is transported seaward and the dissipation rate is much slower. The analysis of the k-equation shows that advective and diffusive transport of turbulence play a major role in the distribution of turbulence under a plunging breaker, while production and dissipation are not in local equilibrium but are of the same order of magnitude. Based on certain approximate analytical approaches and experimental measurements it is shown that turbulence production and viscous dissipation below trough level amount to only a small portion of the wave energy loss caused by wave breaking. It is suggested that the onshore sediment transport produced by swell waves may be tied in a direct way to the unique characteristics of turbulent flows in these waves.     

Is the small-scale turbulence an exclusive breaking product of oceanic internal waves

On the basis of the theoretical research results by the author and the literature published up to date, the analysis and the justification presented in this paper show that the breaking products of oceanic internal waves are not only turbulence, but also the fine-scale near-inertial internal waves (the oceanic reversible finestructure) for inertial waves and the internal solitary waves for internal tides respectively. It was found that the oceanic reversible finestructure may be induced by the effect of the horizontal component f (f = 2Ωcosφ) of the rotation vector on inertial waves. And a new instability of the theoretical shear and strain spectra due to the effect of f occurs at critical vertical wavenumber β c ≈ 0.1 cpm. It happens when the levels of shear and strain of the reversible finestructure are higher than those of inertial waves, which is induced by the effect of f along an "iso-potential-pycnal" of internal wave. If all breaking products of internal waves are taken into account, the average kinetic energy dissipation rate is an order of magnitude larger than the values of turbulence observed by microstructure measurements. The author’s theoretical research results are basically in agreement with those observed in IWEX, DRIFTER and PATCHEX experiments. An important impersonal fact is that on the mean temporal scale of thermohaline circulation these breaking products of internal waves exist simultaneously with turbulence. Because inertial waves are generated by winds at the surface, and internal tides are generated by strong tide-topography interactions, the analysis and justification in this paper support in principle the abyssal recipes Ⅱ:energetics of tidal and wind mixing by Munk & Wunsch in 1998, in despite of the results of microstructure measurements for the turbulent kinetic energy dissipation rate and the diapycnal turbulent eddy diffusivity.     

吕宋海峡不同湍流估算方法的应用对比研究

由于湍流混合直接观测技术难度大、成本高,很大程度限制湍流混合的研究,所以基于温、盐、流资料估算海洋湍流混合的方法应运而生。本文应用在吕宋海峡观测到的23个自由下降微结构湍流剖面仪观测数据和水文观测数据,首次对目前常用的Gregg-Henyey-Polzin(GHP)细尺度参数化、Mackinnon and Gregg(MG)参数化和Thorpe尺度方法进行比较研究,评估它们的适用性。发现GHP参数化方法能够很好地估算吕宋海峡的湍流混合。虽然GHP参数化方法估算的耗散率总体上要偏弱于观测的结果,但估算和观测的差异在2倍以内的结果占71%,与微结构湍流剖面仪观测到的耗散率在水平分布和垂向分布上呈现出相同的分布特征。基于MG参数化方法发现估算的吕宋海峡西侧1200m以浅的耗散率比观测值大,但总体上呈现出相同的分布特征。另外,MG参数化估算与观测差异在2倍以内的结果占58%。表明相比于GHP参数化方法,MG参数化方法的估算值更偏离观测值。Thorpe尺度方法在估算吕宋海峡的耗散率时,估算和观测的差异在2倍以内的结果仅有30%,70%的估算结果与观测结果相差高出1个量级,空间分布上与观测结果差异较大。对比吕宋海峡湍流混合参数化方法的结果表明GHP参数化方法最优,MG参数化方法其次,Thorpe尺度方法相对GHP和MG参数化方法较差。     

Microstructure measurements in natural waters: Methodology and applications

Modern approaches to microstructure data processing, including wavelet denoising, are discussed. The wavelet procedure is applied to small-scale shear signals before estimating the dissipation rate ε and to the temperature/density profiles used to calculate Thorpe scales. Microstructure data obtained on the Mediterranean shelf of Catalonia are used to illustrate various approaches to the Thorpe displacement calculations. It is suggested that the Weibull probability function is an appropriate model for the Thorpe scale distribution. Microstructure measurements from the upper layer of the Boadella reservoir (Catalonia, Spain) support this finding.A new analytical approximation for the 1D Panchev–Kesich spectrum is deduced and the results of ε computation are compared with spectral fitting by the widely used Nasmyth spectrum. Applying the Kraichnan spectral model to compute ε from temperature spectra in the convective-viscous sub-range is examined as an alternative to the Batchelor spectrum. Microstructure measurements taken in Lake Banyoles (Catalonia, Spain) and in the North Atlantic were used for spectral calculations.Statistical analysis of eddy K_b and thermal K_θ diffusivities measured on a shallow shelf of the Black Sea shows the importance of process-orientated domain averaging of the diffusivities in obtaining good correspondence between K_b and K_θ in active turbulent regions. In weakly turbulent, stratified interior layers, the averaged K_b and K_θ differ significantly, which may point to the inapplicability of isotropic formulae used for ε and temperature dissipation χ_θ estimates, as well as to a dependence of the mixing efficiency γ on the Richardson number or in some cases on regions of fossil turbulence.     

On the eddy mixing and energetics of turbulence in a stable atmospheric boundary layer

Some changes in the eddy mixing in the atmospheric boundary layer (ABL) are investigated with the use of the mesoscale RANS turbulence model. It is found that the behavior of parameters of the eddy turbulence mixing is in compliance with the recently obtained data of laboratory and atmospheric measurements. In particular, the flow Richardson number (Ri _f ) during the transient flow to a strongly stable state can behave nonmonotonically, growing with the increasing gradient Richardson number (Ri _g ) to the state of saturation at a certain gradient Richardson number (Ri _g ? 1), which separates two different turbulent regimes: the regimes of strong mixing and weak mixing. An analysis of the energetics based on the balance equations of kinetic and potential turbulence energies shows, in particular, that the weak mixing (Ri _g > 1) is quite capable of transferring momentum. This phenomenon can be explained not only by the fact that the flow is sustained by propagating internal waves, which effectively transfer momentum under strong stratification conditions, but also by the fact that turbulence permanently arises in the free atmosphere and in the deep ocean at Ri _g ? 1.     

Validating a turbulence closure against estuarine microstructure measurements

A simple k–ε turbulence closure is introduced which has no stability functions but instead a Richardson number-dependent turbulent Prandtl number. Its free parameters are determined in a comparison with microstructure observations from a stratified and sheared tidal estuary and laboratory measurements. The closure is able to simulate observed turbulent dissipation rates (ε) and turbulent length scales (l_th) in regions of strong mean shear and small gradient Richardson number (R_g) to within factors of 2–3. It fails in regions of small shear and large R_g, presumably because of the dominance of internal wave-driven mixing. Additional simulations with a k–ε closure with stability functions taken from Canuto et al. [Canuto, V.M., Howard, A., Cheng, Y., Dubovikov, M.S., 2001. Ocean turbulence I: one-point closure model. Momentum and heat vertical diffusivities. J. Phys. Oceanogr. 31, 1413–1426] and with the closure of Baumert and Peters [Baumert, H., Peters, H., 2004. Turbulence closure, steady state, and collapse into waves. J. Phys. Oceanogr. 34, 505–512] show poor performance. Establishing a valid 1:1 comparison of simulated and observed ε and l_th requires nudging the model velocity and density toward observed values because free model integrations quickly diverge from the observations. Steady state gradient Richardson numbers are constrained to a range of 0.18–0.25, while flux Richardson numbers are constrained to the range of 0.1–0.22. The closure output is rather insensitive to such parameter variations. The simulations are sensitive, however, to the treatment of the observed velocity and density used to nudge the model. Good closure performance requires averaging the measured tidal flow over about an hour, a time scale for which conventional numerical models of estuarine circulations should be able to match observed shears. In the closure simulations the TKE balance stays close to a production–dissipation balance. The time rate of change and vertical diffusion of TKE are small, of the same order of magnitude, and vary in magnitude relative to each other systematically across the water column.     

Tidally driven mixing in a numerical model of the ocean general circulation

Astronomical data reveals that approximately 3.5 terawatts (TW) of tidal energy is dissipated in the ocean. Tidal models and satellite altimetry suggest that 1 TW of this energy is converted from the barotropic to internal tides in the deep ocean, predominantly around regions of rough topography such as mid-ocean ridges. A global tidal model is used to compute turbulent energy levels associated with the dissipation of internal tides, and the diapycnal mixing supported by this energy flux is computed using a simple parameterization.The mixing parameterization has been incorporated into a coarse resolution numerical model of the global ocean. This parameterization offers an energetically consistent and practical means of improving the representation of ocean mixing processes in climate models. Novel features of this implementation are that the model explicitly accounts for the tidal energy source for mixing, and that the mixing evolves both spatially and temporally with the model state. At equilibrium, the globally averaged diffusivity profile ranges from 0.3 cm² s⁻¹ at thermocline depths to 7.7 cm² s⁻¹ in the abyss with a depth average of 0.9 cm² s⁻¹, in close agreement with inferences from global balances. Water properties are strongly influenced by the combination of weak mixing in the main thermocline and enhanced mixing in the deep ocean. Climatological comparisons show that the parameterized mixing scheme results in a substantial reduction of temperature/salinity bias relative to model solutions with either a uniform vertical diffusivity of 0.9 cm² s⁻¹ or a horizontally uniform bottom-intensified arctangent mixing profile. This suggests that spatially varying bottom intensified mixing is an essential component of the balances required for the maintenance of the ocean’s abyssal stratification.     

Observation of wave-enhanced turbulence in the near-surface layer of the ocean during TOGA COARE

Dissipation rate statistics in the near-surface layer of the ocean were obtained during the month-long COARE Enhanced Monitoring cruise with a microstructure sensor system mounted on the bow of the research vessel. The vibration contamination was cancelled with the Wiener filter. The experimental technique provides an effective separation between surface waves and turbulence, using the difference in spatial scales of the energy-containing surface waves and small-scale turbulence. The data are interpreted in the coordinate system fixed to the ocean surface. Under moderate and high wind-speed conditions, we observed the average dissipation rate of the turbulent kinetic energy in the upper few meters of the ocean to be 3–20 times larger than the logarithmic layer prediction. The Craig and Banner (J. Phys. Oceanogr. 24 (1994) 2546) model of wave-enhanced turbulence with the surface roughness length from the water side z₀ parameterized according to the Terray et al. (J. Phys. Oceanogr. 26 (1996) 792) formula z₀=cH_s provides a reasonable fit to the experimental dissipation profile, where z is the depth (defined here as the distance to the ocean surface), c≈0.6, and H_s is the significant wave height. In the wave-stirred layer, however, the average dissipation profile deviates from the model (supposedly because of extensive removing of the bubble-disturbed areas close to the ocean surface). Though the scatter of individual experimental dissipation rates (10-min averages) is significant, their statistics are consistent with the Kolmogorov's concept of intermittent turbulence and with previous studies of turbulence in the upper ocean mixed layer.     

夏季黄东海硝酸盐垂向扩散通量的分析

湍流扩散过程导致的硝酸盐垂向输运对海水表层的浮游植物生长和初级生产力的大小有着重要影响。本文基于2018年夏季黄、东海水文环境、硝酸盐浓度和湍动能耗散率的同步、原位数据,分析了海域温度、盐度和硝酸盐的空间分布特征,结果表明营养盐含量丰富的黄海冷水团、长江冲淡水、东海北部底层混合水与黑潮次表层水是影响研究海域硝酸盐分布的主要水团。利用垂向湍扩散硝酸盐通量公式,计算了三个选定断面上的硝酸盐垂向扩散通量,其高值区与湍流扩散系数的高值区的位置基本一致。针对存在明显硝酸盐跃层的站位,计算得到跨硝酸盐跃层的垂向通量F_ND的范围在-9.78—36.60mmol/（m²·d）之间,在黄海冷水团区,夏季温跃层限制了该区营养盐向近表层的湍流垂向扩散;东海北部底层混合水区,湍流垂向扩散向上层补充了大量硝酸盐,促进了跃层之上浮游植物的生长;黑潮次表层水影响海区,夏季中层水体混合较弱,跨跃层的垂向通量也普遍偏低。开展硝酸盐垂向扩散通量的计算与分析,对进一步明确营养盐的输运机制有着重要研究意义。     

Distributions of copepod nauplii and turbulence on the southern flank of Georges Bank: implications for feeding by larval cod (Gadus morhua)

The vertical distributions of copecod nauplii and water properties were sampled at well-mixed and stratified sites on Georges Bank using a pumping system, CTD and in vivo fluorometer during a four day period in late May 1992. At each stratified station at least one sample was taken within the thermocline and the fluorescence maximum, which usually co-occurred. Well-mixed sites had low average concentrations of nauplii, ca 41⁻¹, and showed little variation of abundance with depth. Stratified sites had from 4 to 16 times the integrated (0–50 m) abundance of nauplii compared to well-mixed sites and showed strong vertical patterns of distribution. Maximum concentrations of nauplii, up to 1601⁻¹, were associated with the thermocline at 7 of the 9 stratified stations. At the two remaining stratified sites the naupliar maximum was in the upper mixed layer, sampled at 5 m depth. The encounter rate between early feeding cod (Gadus morhua) larvae and their naupliar prey was calculated with and without turbulence. Turbulence was estimated from two sources: wind stress in the upper layer (calculated from wind observations during our cruise) and tidal shear in the lower layer (estimated initially from a tidal mixing equation). We ultimately replaced the lower layer estimates with turbulence data from a series of measurements made in 1995. The latter are more robust and had the advantage of providing dissipation rates for the pycnocline as well as the lower layer. Theory predicts an increase in encounters between a predator and its prey with the addition of turbulence parameters into standard models of encounter. We combined turbulence profiles with the vertical distribution of nauplii to examine the potential contribution of turbulent kinetic energy to predator-prey encounter rates at various depths in stratified and mixed water columns. Our calculations suggest the following increases due to turbulence at stratified sites on Georges Bank during our cruise: from 34 to 219% in the upper mixed layer, depending on wind speed and depth; approximately 8% in the pycnocline; and approximately 110% below the pycnocline. Mixed sites experience increases of at least 110% (tide only), but greater increases (118–192% in this study) occur when the wind blows because of the combined (spatially overlapped) effects of wind and tidal mixing at these sites. The absolute values for encounter rates and their modification by turbulence are sensitive to a number of assumptions in the models. We used a series of stated assumptions to generate estimates that range from 0.6 to 26.5 prey h⁻¹, depending on geographical location, physical forcing and depth.     

Turbulent mixing in a stratified estuarine tidal channel: Hikapu Reach,Pelorus Sound,New Zealand

Channel constrictions within an estuary can influence overall estuary-sea exchange of salt or suspended/dissolved material. The exchange is modulated by turbulent mixing through its effect on density stratification. Here we quantify turbulent mixing in Hikapu Reach, an estuarine channel in the Marlborough Sounds, New Zealand. The focus is on a period of relatively low freshwater input but where density stratification still persists throughout the tidal cycle, although the strength of stratification and its vertical structure vary substantially. The density stratification increases through the ebb tide, and decreases through the flood tide. During the spring tides observed here, ebb tidal flow speeds reached 0.7?m?s^?1 and the buoyancy frequency squared was in the range 10^?5 to 10^?3?s^?2. Turbulence parameters were estimated using both shear microstructure and velocimeter-derived inertial dissipation which compared favourably. The rate of dissipation of turbulent kinetic energy reached 1?×?10^?6?m²?s^?3 late in the ebb tide, and estimates of the gradient Richardson number (the ratio of stability to shear) fell as low as 0.1 (i.e. unstable) although the results show that bottom-boundary driven turbulence can dominate for periods. The implication, based on scaling, is that the mixing within the channel does not homogenise the water column within a tidal cycle. Scaling, developed to characterise the tidal advection relative to the channel length, shows how riverine-driven buoyancy fluxes can pass through the tidal channel section and the stratification can remain partially intact.     

Diapycnal mixing by meso-scale eddies

The mean available potential energy released by baroclinic instability into the meso-scale eddy field has to be dissipated in some way and Tandon and Garrett [Tandon, A., Garrett, C., 1996. On a recent parameterization of mesoscale eddies. J. Phys. Oceanogr. 26 (3), 406–416] suggested that this dissipation could ultimately involve irreversible mixing of buoyancy by molecular processes at the small-scale end of the turbulence cascade. We revisit this idea and argue that the presence of dissipation within the thermocline automatically requires that a component of the eddy flux associated with meso-scale eddies must be associated with irreversible mixing of buoyancy within the thermocline. We offer a parameterisation of the implied diapycnal diffusivity based on (i) the dissipation rate for eddy kinetic energy given by the meso-scale eddy closure of Eden and Greatbatch [Eden, C., Greatbatch, R.J., 2008. Towards a meso-scale eddy closure. Ocean Modell. 20, 223–239.] and (ii) a fixed mixing efficiency. The implied eddy-induced diapycnal diffusivity (κ) is implemented in a coarse resolution model of the North Atlantic. In contrast to the vertical diffusivity given by a standard vertical mixing scheme, large lateral inhomogeneities can be found for κ in the interior of the ocean. In general, κ is large, i.e. up to o(10) cm²/s, near the western boundaries and almost vanishing in the interior of the ocean.     

Numerical study of turbulence and wave damping induced by vegetation canopies

Vegetation canopies control mean and turbulent flow structure as well as surface wave processes in coastal regions. A non-hydrostatic RANS model based on NHWAVE (Ma et al., 2012) is developed to study turbulent mixing, surface wave attenuation and nearshore circulation induced by vegetation. A nonlinear k − ϵ model accounting for vegetation-induced turbulence production is implemented to study turbulent flow within the vegetation field. The model is calibrated and validated using experimental data from vegetated open channel flow, as well as nonbreaking and breaking random wave propagation in vegetation fields. It is found that the drag-related coefficients in the k − ϵ model C_fk and C_fϵ can greatly affect turbulent flow structure, but seldom change the wave attenuation rate. The bulk drag coefficient C_D is the major parameter controlling surface wave damping by vegetation canopies. Using the empirical formula of Mendez and Losada (2004), the present model provides accurate predictions of vegetation-induced wave energy dissipation. Wave propagation through a finite patch of vegetation in the surf zone is investigated as well. It is found that the presence of a finite patch of vegetation may generate strong pressure-driven nearshore currents, with an onshore mean flow in the unvegetated zone and an offshore return flow in the vegetated zone.     

The effect of standing biomass on flow velocity and turbulence in Spartina alterniflora canopies

Plant-flow interactions on the surface of tidal wetlands result in flow characteristics that are profoundly different from non-vegetated flows. Reductions in mean flow velocity and turbulence, especially the vertical components, limit vertical mixing and may impact a wide range of processes including geochemical exchanges at the sediment water interface, larval recruitment and dispersion, and sediment deposition and retention. The goal of this paper is to quantify horizontal and vertical components of velocity, turbulence intensity and total turbulent kinetic energy in Spartina alterniflora canopies in southeastern North Carolina and to relate flow characteristics to particulate transport on the marsh surface. Another aim of this paper is to assess the extent to which the distribution of standing biomass affects mean flow and turbulence by comparing S. alterniflora data to other canopy types and through a series of canopy manipulations which altered canopy height and stem densities.The results of this study indicate that flow velocity, turbulence intensity, and total turbulent kinetic energy (TKE) are significantly reduced within the vegetated canopy and that this reduction is inversely related to the amount of biomass present in the water column. Within the canopy, approximately 50% of the initial mean velocity and TKE is reduced within 5 m of the canopy edge. Within the canopy, mean velocity and TKE_horiz usually exceeded vertical velocity or TKE_vert and the vertical components of flow were attenuated more strongly than the horizontal. These results suggest that within the vegetation, turbulence contributes more to lateral advection than to vertical mixing. As a result, total suspended solid concentrations were shown to decrease logarithmically with distance from the canopy edge and to decrease at a faster rate in more densely vegetated regions of the canopy (i.e. lower TKE_vert) as compared to areas of sparser vegetation (i.e. higher TKE_vert).     

Assessment of turbulence closure models for resonant inertial response in the oceanic mixed layer using a large eddy simulation model

Large eddy simulation (LES) of the resonant inertial response of the upper ocean to strong wind forcing is carried out; the results are used to evaluate the performance of each of the two second-order turbulence closure models presented by Mellor and Yamada (Rev Geophys Space Phys 20:851–875, 1982) (MY) and by Nakanishi and Niino (J Meteorol Soc Jpn 87:895–912, 2009) (NN). The major difference between MY and NN is in the formulation of the stability functions and the turbulent length scale, both strongly linked with turbulent fluxes; in particular, the turbulent length scale in NN, unlike that in MY, is allowed to decrease with increasing density stratification. We find that MY underestimates and NN overestimates the development of mixed layer features, for example, the strong entrainment at the base of the oceanic mixed layer and the accompanying decrease of sea surface temperature. Considering that the stability functions in NN perform better than those in MY in reproducing the vertical structure of turbulent heat flux, we slightly modify NN to find that the discrepancy between LES and NN can be reduced by more strongly restricting the turbulent length scale with increasing density stratification.     

南海南部上层的湍流混合特征

A turbulent microstructure experiment was undertaken at a low latitude of 10°N in the South China Sea in late August 2012. The characteristics of the eddy diffusivity above 650 m were analyzed, which is one order of magnitude larger than that in the open ocean at that low latitude. Enhanced eddy diffusivities by strong shears and sharp changes in topography were observed. The strongest eddy diffusivity occurred in the mixed layer, and it reached O(10–2 m2/s). Strong stratification in the thermocline inhibited the penetration of surface eddy diffusivities through the thermocline, where the mixing was weakest. Below the thermocline, where the background eddy diffusivity was approximately O(10–6 m2/s), the eddy diffusivity increased with depth, and its largest value was O(10–3 m2/s).