浪致混合对亚热带冬季海洋混合强度的影响

上层海洋在全球气候系统中起着至关重要的作用。对上层海洋层结及混合的模拟偏差一直是海洋和气候数值模式发展中悬而未决的问题。本文首先评估了CMIP5中45个模式对上层海洋层结模拟的偏差,确认了冬季亚热带地区海洋模式垂向混合偏强。随后,基于自然资源部第一海洋研究所地球系统模式(FIO-ESM v1.0),分别开展了1986−2005年期间包含和关闭海浪垂向混合情况下的数值实验,分析浪致混合对亚热带冬季海洋混合强度模拟的影响及机制。发现浪致混合使得气候模式中亚热带海域冬季的海洋层结增强,增强的层结使上层海洋更加稳定。首次揭示了增加浪致混合反而降低了海洋总体的垂向混合率：浪致混合使北半球冬季亚热带海域混合率从无浪实验的227 cm²/s降低到有浪实验的178 cm²/s,降低了21.6%;南半球冬季亚热带海域混合率从无浪实验的189 cm²/s降低到有浪实验的165 cm²/s,降低了12.7%。进一步分析发现,浪致混合主要是通过增加冬季亚热带海域上层海洋的热含量从而强化了海洋的层结,最终改善了气候模式对上层海洋混合的模拟。     

Sensitive study of the long and short surface wave-induced vertical mixing in a wave-circulation coupled model

The previous studies by the MASNUM research team have shown the effectiveness of the wave-induced mixing (Bv) in improving the simulation of upper-ocean thermal structure. The mechanisms of Bv are further investigated by incorporating different Bv products into the MASNUM wave-circulation coupled model. First, experiments were designed to explore the effects of Bv, which contain the contributions at different wave lengths (l). The results of three experiments, the non-Bv case, the short-wave case (l <300 m), and the long-wave case (l >300 m) are compared, and it is found that the long waves are the most important component for Bv to generate mixing in the upper ocean. As the swell plays dominant role in mixing, the parameterization of Bv into wind may be not a proper way. Second, Bv effects at different time-scales, including daily and monthly, were examined. The results show that the monthly averaged Bv has larger impact than the daily averaged Bv, especially in summer.     

海浪搅拌混合对北太平洋海表面温度模拟的影响

利用NCEP再分析风场驱动WAVEWATCH III海浪模式对北太平洋海域的海浪过程进行模拟,利用浮标观测资料对模拟出的海浪要素有效波高进行验证,发现他们之间具有很好的一致性。基于模式输出的有效波高等波浪要素,利用特征波参数化理论,在海洋环流模式中引入海浪搅拌混合作用,分析其对北太平洋海表面温度模拟的影响,初步数值模拟结果表明,sbPOM模式在考虑海浪搅拌混合作用以后,模拟精度进一步提升,这对提供一个准确的大气模式下边界条件具有重要作用。     

Effects of the surface wave-induced mixing on circulation in an isopycnal-coordinate oceanic circulation model

The influence of the nonbreaking surface wave-induced mixing under the mixed layer on the oceanic circulation was investigated using an isopycnal-coordinate oceanic circulation model. The effect of the wave-induced mixing within the mixed layer was eliminated via a bulk mixed layer model. The results show that the wave-induced mixing can penetrate through the mixed layer and into the oceanic interior. The wave-induced mixing under the mixed layer has an important effect on the distribution of temperature of the upper ocean at middle and high latitudes in summer, especially the structure of the seasonal thermocline. Moreover, the wave-induced mixing can affect the oceanic circulation, such as western boundary currents and the North Equatorial Currents through changes of sea surface height associated with the variation of the thermal structure of the upper ocean.     

浪致混合对2016年北太平洋海表温度季节性预测的影响

上层海洋通过海气交换影响大气-海洋耦合系统,海浪引起的垂向混合影响上层海洋结构,从而在气候预测过程中发挥着重要的作用。本文基于国家海洋局第一海洋研究所地球系统模式（FIO-ESM）,以2016年为例,分别开展了耦合和关闭海浪模式情况下的短期气候预测实验,分析浪致混合对北太平洋海表温度（SST）季节性预测的影响。通过对模式预测的SST异常（SSTA）进行定量评估发现,浪致混合能够显著降低北太平洋高纬度海区预测误差,在（45°N,150°E）附近海区SSTA改善可达1℃,气候模式能够更好地预测SSTA的经向分布特征,特别是能够准确地反映25°~45°N海区SSTA分布特征。通过分析有浪和无浪两个实验的热收支贡献发现,垂向混合是导致上层海洋温度差异的主导影响因子。海浪通过改变垂向混合,使2016年北太平洋SST在高纬度海区大幅降低,在低纬度海区略有升高,最终提升了模式对北太平洋SST的季节性预测能力。     

Comparison of turbulence schemes for prediction of wave-induced near-bed sediment suspension above a plane bed

Based on a wave bottom boundary layer model and a sediment advection-diffusion model, seven turbulence schemes are compared regarding their performances in prediction of near-bed sediment suspension beneath waves above a plane bed. These turbulence algorithms include six empirical eddy viscosity schemes and one standard two-equation k-ε model. In particular, different combinations of typical empirical formulas for the eddy viscosity profile and for the wave friction factor are examined. Numerical results are compared with four laboratory data sets, consisting of one wave boundary layer hydrodynamics experiment and three sediment suspension experiments under linear waves and the Stokes second-order waves. It is shown that predictions of near-bed sediment suspension are very sensitive to the choices of the empirical formulas in turbulence schemes. Simple empirical turbulence schemes are possible to perform equally well as the two-equation k-ε model. Among the empirical schemes, the turbulence scheme, combining the exponential formula for eddy viscosity and Swart formula for wave friction factor, is the most accurate. It maintains the simplicity and yields identically good predictions as the k-ε model does in terms of the wave-averaged sediment concentration.     

Numerical study on water waves and wave-induced longshore currents in Obaky coastal water

In this paper, the water waves and wave-induced longshore currents in Obaky coastal water which is located at the Mediterranean coast of Turkey were numerically studied. The numerical model is based on the parabolic mild-slope equation for coastal water waves and the nonlinear shallow water equation for the wave-induced currents. The wave transformation under the effects of shoaling, refraction, diffraction and breaking is considered, and the wave provides radiation stresses for driving currents in the model. The numerical results for the water wave-induced longshore currents were validated by the measured data to demonstrate the efficiency of the numerical model. Then the water waves and longshore currents induced by the waves from main directions were numerically simulated and analyzed based on the numerical results. The numerical results show that the movement of the longshore currents was different while the wave propagated to a coastal zone from different directions.     

两种湍流参数化方案针对水体温度演变过程的模拟与比较

依据山东跋山水库的观测数据,采用MY2.5参数化方案和κ-ε参数化方案在不同的表面边界条件下模拟观测的温度剖面演化过程。模拟结果显示,在风应力强迫的表面边界条件下,采用这两种参数化方案的模拟结果差别不大;当考虑Mellor和Blumberg的波浪破碎表面边界条件时,采用MY2.5参数化方案的模拟结果没有明显的改善,但采用κ-ε参数化方案的模拟结果有明显的改善,得到与观测相符的模拟结果。     

两种湍流参数化方案针对水体温度演变过程的模拟与比较

依据山东跋山水库的观测数据,采用MY2.5参数化方案和κ-ε参数化方案在不同的表面边界条件下模拟观测的温度剖面演化过程。模拟结果显示,在风应力强迫的表面边界条件下,采用这两种参数化方案的模拟结果差别不大;当考虑Mellor和Blumberg的波浪破碎表面边界条件时,采用MY2.5参数化方案的模拟结果没有明显的改善,但采用κ-ε参数化方案的模拟结果有明显的改善,得到与观测相符的模拟结果。     

上混合中剪切湍流和朗缪尔环流动力特征差异

Large eddy simulation(LES) is used to investigate contrasting dynamic characteristics of shear turbulence(ST)and Langmuir circulation(LC) in the surface mixed layer(SML). ST is usually induced by wind forcing in SML. LC can be driven by wave-current interaction that includes the roles of wind, wave and vortex forcing. The LES results show that LC suppresses the horizontal velocity and greatly modifies the downwind velocity profile, but increases the vertical velocity. The strong downwelling jets of LC accelerate and increase the downward transport of energy as compared to ST. The vertical eddy viscosity Km of LC is much larger than that of ST. Strong mixing induced by LC has two locations. They are located in the 2ds–3ds(Stokes depth scale) and the lower layer of the SML,respectively. Its value and position change periodically with time. In contrast, maximum Km induced by ST is located in the middle depth of the SML. The turbulent kinetic energy(TKE) generated by LC is larger than that by ST. The differences in vertical distributions of TKE and Km are evident. Therefore, the parameterization of LC cannot be solely based on TKE. For deep SML, the convection of large-scale eddies in LC plays a main role in downward transport of energy and LC can induce stronger velocity shear(S2) near the SML base. In addition, the large-scale eddies and S2 induced by LC is changing all the time, which needs to be fully considered in the parameterization of LC.     

The role of surface waves in the ocean mixed layer

Previously, most ocean circulation models have overlooked the role of the surface waves. As a result, these models have produced insufficient vertical mixing, with an under - prediction of the ,nixing layer （ML） depth and an over - prediction of the sea surface temperature （SST）, particularly during the summer season. As the ocean surface layer determines the lower boundary conditions of the atmosphere, this deficiency has severely limited the performance of the coupled ocean - atmospheric models and hence the climate studies. To overcome this shortcoming, a new parameterization for the wave effects in the ML model that will correct this systematic error of insufficient mixing. The new scheme has enabled the mixing layer to deepen, the surface excessive heating to be corrected, and an excellent agreement with observed global climatologic data. The study indicates that the surface waves are essential for ML formation, and that they are the primer drivers of the upper ocean dynamics; therefore, they are critical for climate studies.     

Influence of surface forcing on near-surface and mixing layer turbulence in the tropical Indian Ocean

An autonomous upwardly-moving microstructure profiler was used to collect measurements of the rate of dissipation of turbulent kinetic energy (ε) in the tropical Indian Ocean during a single diurnal cycle, from about 50 m depth to the sea surface. This dataset is one of only a few to resolve upper ocean ε over a diurnal cycle from below the active mixing layer up to the air–sea interface. Wind speed was weak with an average value of ~5 m s⁻¹ and the wave field was swell-dominated. Within the wind and wave affected surface layer (WWSL), ε values were on the order of 10⁻⁷–10⁻⁶ W kg⁻¹ at a depth of 0.75 m and when averaged, were almost a factor of two above classical law of the wall theory, possibly indicative of an additional source of energy from the wave field. Below this depth, ε values were closer to wall layer scaling, suggesting that the work of the Reynolds stress on the wind-induced vertical shear was the major source of turbulence within this layer. No evidence of persistent elevated near-surface ε characteristic of wave-breaking conditions was found. Profiles collected during night-time displayed relatively constant ε values at depths between the WWSL and the base of the mixing layer, characteristic of mixing by convective overturning. Within the remnant layer, depth-averaged values of ε started decaying exponentially with an e-folding time of 47 min, about 30 min after the reversal of the total surface net heat flux from oceanic loss to gain.     

Energetics of wind-induced turbulent mixing in the ocean

The pattern and magnitude of the global ocean overturning circulation is believed to be strongly controlled by the distribution of diapycnal diffusivity below 1000 m depth. Although wind stress fluctuation is a candidate for the major energy sources of diapycnal mixing processes, the global distribution of wind-induced diapycnal diffusivity is still uncertain. It has been believed that internal waves generated by wind stress fluctuations at middle and high latitudes propagate equatorward until their frequency is twice the local inertial frequency and break down via parametric subharmonic instabilities, causing diapycnal mixing. In order to check the proposed scenario, we use a vertically two-dimensional primitive equation model to examine the spatial distribution of “mixing hotspots” caused by wind stress fluctuations. It is shown that most of the wind-induced energy fed into the ocean interior is dissipated within the top 1000 m depth in the wind-forced area and the energy dissipation rate at low latitudes is very small. Consequently, the energy supplied to diapycnal mixing processes below 1000 m depth falls short of the level required to sustain the global ocean overturning circulation.     

珊瑚礁礁坪宽度对波浪传播变形及增水影响的实验研究

通过在波浪水槽中进行一系列物理模型实验, 研究珊瑚礁礁坪宽度变化对珊瑚礁海岸附近波浪传播变形及礁坪上波浪增水的影响。物理实验采用理想化的珊瑚礁模型, 测试了3种礁坪宽度下的一系列不规则波工况。实验结果分析表明: 波浪沿礁传播过程中, 短波持续衰减; 低频长波波高沿礁逐渐增大, 直到海岸线附近达到最大; 随着礁坪宽度的增加, 海岸线附近的短波波高呈下降趋势, 低频长波波高的变化规律不显著; 礁坪上的波浪增水受礁坪宽度变化的影响不明显; 通过对海岸线附近的波浪进行频谱分析发现, 礁坪上低频长波的运动存在着一阶共振模式, 且共振放大效应强度受礁坪水深、入射波峰周期和礁坪宽度共同影响。     

Influence of the earth rotation on the interfacial wave solutions and wave-induced tangential stress in a twolayer fluid system

Interfacial waves and wave-induced tangential stress are studied for geostrophic small amplitude waves of two-layer fluid with a top free surface and a flat bottom. The solutions were deduced from the general form of linear fluid dynamic equations of two-layer fluid under the f-plane approximation, and wave-induced tangential stress were estimated based on the solutions obtained. As expected, the solutions derived from the present work include as special cases those obtained by Sun et al. (2004. Science in China, Ser. D, 47(12):1147-1154) for geostrophic small amplitude surface wave solutions and wave-induced tangential stress if the density of the upper layer is much smaller than that of the lower layer. The results show that the interface and the surface will oscillate synchronously, and the influence of the earth''s rotation both on the surface wave solutions and the interfacial wave solutions should be considered.     

基于CESM1.2.2模拟的浪致混合对末次冰盛期和工业革命前气候的影响

海洋上层垂向混合在模式中发挥重要的作用,以往的研究表明垂向混合的不足使得模拟的海洋温度和混合层深度与观测存在显著偏差。前人提出一种修正方案,考虑波浪产生的垂向混合,将由表面风作用下产生的波浪这样一个实际物理过程的湍混合进行参数化,其结果被证实能够显著提高模式模拟和预报的准确性。本文首次将浪致混合引入海气耦合的古气候模式,基于末次冰盛期和工业革命前2种不同的气候条件,探究浪致混合在海气耦合模式中的作用。在不同气候背景下,由于风场强度的不同,导致末次冰盛期浪致混合的强度小于工业革命前,但2个气候时期都体现出中纬度混合强度最大的特点。将浪致混合加入到气候模式中,模拟结果表明：中纬度海域2个时期都出现海表面降温而次表层升温的现象,但末次冰盛期的表面降温强度弱于工业革命前状态;不同月份下的模拟结果显示,在南北半球的夏季,海洋表层温度的降温最为显著。中纬度海域海洋上混合层深度在年平均条件下2个气候背景时期都出现加深现象,但末次冰盛期的加深程度弱于工业革命前;不同月份下的模拟结果显示,在南北半球的冬季,混合层加深的变化达到极值。另一方面,在高纬度海域,末次冰盛期的海表面温度出现了显著升高,这是由于浪...     

Numerical study of the meridional overturning circulation with “mixing hotspots” in the Pacific Ocean

Using an idealized ocean general circulation model, we examine the effect of “mixing hotspots” (localized regions of intense diapycnal mixing) predicted based on internal wave-wave interaction theory (Hibiya et al., 2006) on the meridional overturning circulation of the Pacific Ocean. Although the assumed diapycnal diffusivity in the mixing hotspots is a little larger than the predicted value, the upwelling in the mixing hotspots is not sufficient to balance the deep-water production; out of 17 Sv of the downwelled water along the southern boundary, only 9.2 Sv is found to upwell in the mixing hotspots. The imbalance as much as 7.8 Sv is compensated by entrainment into the surface mixed layer in the vicinity of the downwelling region. As a result, the northward transport of the deep water crossing the equator is limited to 5.5 Sv, much less than estimated from previous current meter moorings and hydrographic surveys. One plausible explanation for this is that the magnitude of the meridional overturning circulation of the Pacific Ocean has been overestimated by these observations. We raise doubts about the validity of the previous ocean general circulation models where diapycnal diffusivity is assigned ad hoc to attain the current magnitude suggested from current meter moorings and hydrographic surveys.     

西北太平洋暖池区台风对海表盐度的影响

分析了西北太平洋暖池区2002和2003年夏季ARGO浮标得到的次表层温度和盐度剖面,结果表明大多数台风经过暖池区时,会引起海面盐度下降,这与Kwon和Riser等在大西洋观测到的飓风过后海面盐度上升的结论不同,表明西北太平洋暖池区特有的上层结构以及台风在此海域的降雨与大西洋不同,结果对研究西北太平洋暖池区的混合层混合和热交换过程有重要意义.     

湍流数值模拟中封闭模式应用的局限性

通过对琼州海峡的潮流场特征进行数值模拟，指出了选择不同的特征混合长度表达式对数值模拟结果的影响，表明了基于特征混合长度理论的流封闭模式在近海湍流数值模拟中应用的局限性。     

近岸波浪和波生流环境中泥沙输运的数值模拟

Owing to lack of observational data and accurate definition,it is difficult to distinguish the Kuroshio intrusion water from the Pacific Ocean into the South China Sea(SCS).By using a passive tracer to identify the Kuroshio water based on an observation-validated three-dimensional numerical model MITgcm,the spatio-temporal variation of the Kuroshio intrusion water into the SCS has been investigated.Our result shows the Kuroshio intrusion is of distinct seasonal variation in both horizontal and vertical directions.In winter,the intruding Kuroshio water reaches the farthest,almost occupying the area from 18°N to 23°N and 114°E to 121°E,with a small branch flowing towards the Taiwan Strait.The intrusion region of the Kuroshio water decreases with depth gradually.However,in summer,the Kuroshio water is confined to the east of 118°E without any branch reaching the Taiwan Strait;meanwhile the intrusion region of the Kuroshio water increases from the surface to the depth about 205 m,then it decreases with depth.The estimated annual mean of Kuroshio Intrusion Transport(KIT) via the Luzon Strait is westward to the SCS in an amount of –3.86×106 m3/s,which is larger than the annual mean of Luzon Strait Transport(LST) of –3.15×106 m3/s.The KIT above 250 m accounts for 60%–80% of the LST throughout the entire water column.By analyzing interannual variation of the Kuroshio intrusion from the year 2003 to 2012,we find that the Kuroshio branch flowing into the Taiwan Strait is the weaker in winter of La Ni?a years than those in El Ni?o and normal years,which may be attributed to the wind stress curl off the southeast China then.Furthermore,the KIT correlates the Ni?o 3.4 index from 2003 to 2012 with a correlation coefficient of 0.41,which is lower than that of the LST with the Ni?o 3.4 index,i.e.,0.78.