基于高分辨率风场的海洋近惯性能通量计算——时空特征及其影响因素

基于高分辨率CFSR(climate forecast system reanalysis)风场资料、气候态海洋混合层厚度资料和卫星高度计海面高度异常资料,本文估计了大气风场向全球海洋混合层的近惯性能通量和近惯性能量输入功率,并探究了混合层厚度、风场时间分辨率、经验衰减系数和中尺度涡旋涡度对近惯性能通量和能量输入功率的影响。浮标实测风场和流速表明,本文所用的风场和阻尼平板模型可用于估计风场向全球海洋的近惯性能通量。本文计算得到的大气向全球海洋输入近惯性能量的功率为0.56TW(1TW=10~(12)W),其中北半球贡献0.22TW,南半球贡献0.34TW。在时间上,风场的近惯性能通量呈现各个半球冬季最强、夏季最弱的特征,这和西风带风场的季节变化有关。在空间上,近惯性能通量的高值海域为南、北半球西风带海洋,尤其是南大洋。混合层厚度和风场空间不均匀性使得西风带近惯性能通量呈现纬向变化,即海盆西部强于海盆东部。风场时间分辨率对近惯性能通量的估计至关重要,低时间分辨率风场对近惯性能通量的低估达到13%—30%。阻尼平板模型中的经验衰减系数对近惯性能通量估计的影响不超过5%。中尺度涡旋涡度仅改变近惯性能通量的空间分布,而对全球近惯性能量输入功率的影响可以忽略。     

An approximate solution of wave-modified Ekman current for gradually varying eddy viscosity

An approximate steady solution of the wave-modified Ekman current is presented for gradually varying eddy viscosity by using the WKB method with the variation of parameters technique. The parameters involved in the solution can be determined by the two-dimensional wavenumber spectrum of ocean waves, wind speed, the Coriolis parameter and the densities of air and water. The solution reduces to the exact solution when the eddy viscosity is taken as a constant. As illustrative examples, for a fully developed wind-generated sea with different wind speeds and a few proposed gradually varying eddy viscosities, the current profiles calculated from the approximate solutions are compared with those of the exact solutions or numerical ones by using the Donelan and Pierson wavenumber spectrum, the WAM wave model formulation for wind input energy to waves, and wave energy dissipation converted to currents. It is shown that the approximate solution presented has an elegant form and yet would be valid for any given gradually varying eddy viscosity. The applicability of the solution method to the real ocean is discussed following the comparisons with published observational data and with the results from a large eddy simulation of the Ekman layer.     

Coupling winds to ocean surface currents over the global ocean

A Wind stress–Current Coupled System (WCCS) consisting of the HYbrid Coordinate Ocean Model (HYCOM) and an improved wind stress algorithm based on Donelan et al. [Donelan, W.M., Drennan, Katsaros, K.B., 1997. The air–sea momentum flux in mixed wind sea and swell conditions. J. Phys. Oceanogr. 27, 2087–2099] is developed by using the Earth System Modeling Framework (ESMF). The WCCS is applied to the global ocean to study the interactions between the wind stress and the ocean surface currents. In this study, the ocean surface current velocity is taken into consideration in the wind stress calculation and air–sea heat flux calculation. The wind stress that contains the effect of ocean surface current velocity will be used to force the HYCOM. The results indicate that the ocean surface velocity exerts an important influence on the wind stress, which, in turn, significantly affects the global ocean surface currents, air–sea heat fluxes, and the thickness of ocean surface boundary layer. Comparison with the TOGA TAO buoy data, the sea surface temperature from the wind–current coupled simulation showed noticeable improvement over the stand-alone HYCOM simulation.     

北太平洋北部的风能输入门户

风不仅驱动了上层海洋的环流,也是深层海洋运动的主要能量来源.本文主要研究了北太平洋北部的风能输入的季节性分布特征和年际变化趋势,包括风向表面波、表层地转流和表层非地转流的能量输入.基于SODA3数据的结果表明,风能输入门户随季节变化显著,其中黑潮延伸区是冬季门户,副极地流涡是春、秋季门户,大洋东边界则是夏季门户,能量输...     

A theoretical reason to expect inviscid western boundary currents in realistic oceans

The conventional picture of an ocean gyre, based on an ocean with vertical sidewalls, assumes a balance between an input of vorticity by wind stress curl, and a viscous flux of vorticity through the boundary at the same latitude, resulting from a viscous western boundary current which may be significantly modified by nonlinear terms. Potential interactions with topography are also commonly acknowledged as a possible complicating factor. In this idealized picture, the zonal momentum balance is taken to be geostrophic, as numerous model analyses confirm. A theoretical argument is given here which shows that, in an ocean with sloping sidewalls, this geostrophic balance results in bottom pressure torques which balance the wind stress curl at each latitude. This removes the requirement for a viscous western boundary current at each latitude suggesting that the dynamics within a western boundary current may be essentially inviscid. While inviscid western boundary currents have already been found in certain idealized systems, and in one set of diagnostics from an eddy permitting model, the generality of the argument presented here gives a strong reason to believe that these are not special cases. Inviscid western boundary currents are in fact the rule, and the vertical sidewall case is an unrealistic exception.     

Global estimates of wind energy input to subinertial motions in the Ekman-Stokes layer

By incorporating the wave-induced Coriolis-Stokes forcing into the classical Ekman model, the wind energy input to the Ekman-Stokes layer is investigated, with an emphasis on the surface wave effects when the direction of Stokes drift deviates from that of wind stress. Theoretical analysis of the kinetic energy balance of the Ekman-Stokes layer shows that the total wind energy input consists of the direct wind energy input and the wave-induced energy input. Details of the direct wind and wave-induced energy input are discussed. Based on the ECMWF ERA-40 Re-Analysis wind stress and surface wave data, the global total wind energy input to subinertial motions in the Ekman-Stokes layer is estimated at 2.19 TW, including 0.26 TW (12%) wave-induced energy input and 1.93 TW (88%) direct wind energy input. The effect of sea-ice coverage on the energy input to the Ekman-Stokes layer is also considered. It is shown that the global total energy input could be overestimated by 0.08 TW (about 4%) without taking the sea-ice coverage into account.     

黄海海洋对台风“灿鸿”外围过程响应的观测研究

2015年9月台风“灿鸿”经过黄海海域的QF111浮标,该浮标观测到黄海海洋的温度和流速响应。受台风影响,黄海海表温度有明显的下降,降温幅度在2~4°C,且在台风路径右侧形成了海表面温度低于20°C的低温斑块。流速的响应主要表现为台风过境期间流速增强,达1.2~1.5 m/s,呈现出全水深强化的特征。风一直对海洋有能量输入,且在台风过境期间风输入海洋的能量最大,旋转功率谱的分析结果表明台风输入海洋的能量向海洋内部传播。黄海海域由台风激发出的近惯性振荡衰减较快,近惯性振荡e折时间尺度约为2 d。     

基于SODA资料的南海表层风能输入的空间分布与长期趋势研究

海面风不仅是驱动上层海洋运动的主要动力, 其能量也是维持海洋表层流动的主要机械能来源。为了分析南海表层流风能输入的变化, 用SODA(Simple Ocean Data Assimilation)(1901—2010)资料估算了风向南海表层流(表层地转流+表层非地转流)的能量输入。结果表明, 风向南海表层流、表层地转流和表层非地转流输入的能量总体均呈减少趋势, 110年间分别减小了约56%、65%和49%。导致风能输入减小的最主要因素是风应力的减弱(减小了35%)。由于南海受季风系统的控制, 风向表层流及其各成分输入的能量呈现出显著的季节性变化。冬季风能输入最强, 高值区位于南海西部及北部区域, 呈一个显著的“回力镖”状结构。这些结果对深入认识南海环流具有理论意义。     

海气动量通量研究综述

海气界面动量通量也称为风应力,是海流和表面海浪的主要驱动力,是海洋从大气获得动量的重要途径。因此,合理可靠的海洋表面风应力的参数化对于海洋、大气和波浪以及气候模式的准确预报都具有非常重要的科学意义和实用价值。对风应力拖曳系数的参数化是风应力参数化的主要内容。近来的观测发现,风应力拖曳系数随着风速的增加出现了先增后减的趋势,同时还与海面的波浪状态以及海流有关。基于观测或理论分析,目前已经得到了一系列的风应力拖曳系数计算方法或公式,有的考虑了海浪的作用,有的没有,但这些方案大都是适合中低风速,在高风速下的适用性还有待检验。本文回顾了目前在海气动量通量观测和参数化方面的研究进展,并建议应增加高风速下风速、海流以及海浪等的同步观测,以进一步完善风应力参数化方案。     

2000—2015年夏秋季孟加拉湾湾口区涡流相互作用能量学特征

文章主要使用全球简单海洋资料同化分析系统(Simple Ocean Data Assimilation, SODA)产出的海洋再分析数据产品和美国国家环境预报中心(National Centers for Environmental Prediction, NCEP)发布的风场资料, 通过能量学方法分析了2000—2015年夏季至秋季(6—11月)孟加拉湾涡-流相互作用特征在不同印度洋偶极子(Indian Ocean Dipole, IOD)事件发生年的表现。结果表明, 在IOD负位相年更强的西南季风背景下, 涡动能和涡势能的量值均较大, 海洋不稳定过程更多地将平均流场的能量输向涡旋场, IOD正位相年反之。另外, 研究发现孟加拉湾湾口区的涡动能在个别年份会发展出一种与气候态存在显著异常的空间分布, 即在个别年份湾口中央海域异常出现涡动能极大值。通过对出现该异常现象最显著的2010年进行个例分析, 发现当年的孟加拉湾海表风场发展出一个气旋式环流异常, 显著地改变了海洋上层环流形态, 极大地影响了平均流场与涡旋场之间的相互作用。进一步对维持涡动能平衡的各做功项进行诊断后发现, 湾口异常海域涡动能年际变化的主要影响因素为海洋内部的压强做功, 其次是正压不稳定过程和平流的做功, 海表风应力做功项贡献较小。     

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.     

Development of a regional ocean model for New Zealand

A regional ocean model with a horizontal resolution of 1/6° encompassing the New Zealand Exclusive Economic Zone is described. The regional model successfully downscaled solutions from a high resolution, global, coupled model HadCEM. Transport estimates from the global and regional models were compared with observations, and both models supported largely consistent, climatological mean solutions. The regional model used monthly mean forcing at the surface. Nevertheless, the regional model eddy kinetic energy (EKE) spatial patterns compared favourably with long‐term mean satellite altimetric estimates, although the modelled background EKE amplitudes were much lower than observed. A series of permanent eddies associated with the western boundary current system around the top of the North Island of New Zealand were reproduced, and an eddy adjacent to Norfolk Ridge was identified in both the global and regional models. The western boundary current system around the North Island of New Zealand and the associated eddies were the most sensitive components of the model solutions, being influenced by initial conditions, wind forcing, and the model domain size.     

Seasonal and interannual variations of the upper ocean energetics between Tasmania and Antarctica

Nine years of Topex/Poseidon and ERS satellite altimetry and XBT data from the SURVOSTRAL program were used to analyze the seasonal and interannual variations of the eddy energetics in terms of its spatial distribution and relation with the upper ocean heat content. Eddy kinetic energy is calculated in two frequency bands one associated with transient and the other with low-frequency variability. The two eddy components have distinct geographical distribution. At the SURVOSTRAL line, the transient eddy energy is twice the low-frequency energy, with maximum transient energy occurring during the austral summer period and maximum low-frequency energy in winter. The site is one of growing eddy energy. Eddy momentum flux is northward over the SURVOSTRAL line, and the summertime eddy heat flux is poleward across the Subantarctic and Subtropical Fronts, and equatorward either side of the fronts. Eddy fluxes are strongly influenced by their position relative to the bathymetry and the mean current.     

风生近惯性能通量和地形粗糙度对海洋内部混合的影响

使用细尺度参数化方法和2015—2019年全球的Argo温盐剖面资料,分析了风生近惯性能通量和地形粗糙度对全球海洋300—600m深度的涡流扩散系数的影响。结果表明,在30°—45°N区域,月均涡流扩散系数与月平均风生近惯性能通量随时间的变化趋势较为一致,相关系数可达0.43,前者滞后1个月,与后者的相关系数可达0.65,但在其他区域二者的变化趋势相差较大;相较于中纬度和北半球,低纬度和南半球的地形粗糙度与扩散系数的相关关系更好。基于这些分析结果,拟合并建立了30°—45°N区域300—600m深度的涡流扩散系数与风生近惯性能通量和地形粗糙度的关系式。而且,用此关系式和细尺度参数化方法计算出来的扩散系数平均量级差异为0.47,且91%的值偏差小于一个量级。据此,我们给出了1—12月30°—45°N太平洋区域的涡流扩散系数的网格化结果。     

应用“跨零-能量”法估算海洋波浪再生能资源

海洋波浪能是一种可持续和零污染的再生能源,随着自然能源日渐匮乏和全球气候变化,人类对波浪能的开发利用显得尤其重要。本文应用"跨零—能量"波浪分析的新方法,首次推导了不规则波浪的波能流垂向分布及其理论计算公式。本文研究发现,浅水波能流具有均匀的垂向分布,深水波能流集中于海表层的水体中,过渡区波能流的垂向分布介于浅水和深水波能的分布之间。研究还发现,目前海洋波能流的估算方法和现今波浪发电装置的波能采集深度缺少相关性,过高估算了海洋波浪的可发电资源。     

Aspects of Deep Ocean Mixing

The turbulent motions responsible for ocean mixing occur on scales much smaller than those resolved in numerical simulations of oceanic flows. Great progress has been made in understanding the sources of energy for mixing, the mechanisms, and the rates. On the other hand, we still do not have adequate answers to first order questions such as the extent to which the thermohaline circulation of the ocean, and hence the earth's climate, is sensitive to the present mixing rates in the ocean interior. Internal waves, generated by either wind or flow over topography, appear to be the principle cause of mixing. Mean and eddy flows over topography generate internal lee waves, while tidal flows over topography generate internal tides. The relative importance of these different internal wave sources is unknown. There are also great uncertainties about the spatial and temporal variation of mixing. Calculations of internal tide generation are becoming increasingly robust, but we do not know enough about the subsequent behavior of internal tides and their eventual breakdown into turbulence. It does seem, however, that most internal tide energy flux is radiated away from generation sites as low modes that propagate over basin scales. The mechanisms of wave-wave interaction and topographic scattering both act to transfer wave energy from low modes to smaller dissipative scales. This revised version was published online in August 2006 with corrections to the Cover Date.     

基于潜标观测的吕宋海峡以东深海海流的低频变异

本文利用在菲律宾海布放的一套锚系潜标获取的长时间海流和水温观测数据,分析了吕宋海峡以东的深海海洋环境特征,着重阐释了该海域海流的全水深垂向结构及其低频变化特征。研究表明,表层(100～160 m)平均流向为西偏北,流速约为12.5 cm/s;中层(810 m)的平均流为西向,流速为2.6 cm/s;深层(1 550 m和2 560 m)的平均流速在1 cm/s以内,近底(4 040 m)的流向为较稳定的西南向,流速为2.3 cm/s。上层海流的动能比中层和深层大1～2个量级,总动能、平均动能、涡动动能均在表层最大,中层次之、深层最小,各层次涡动能均大于平均动能。中上层海流的低频变化具有极高的相似性,全年为81～85 d的周期振荡;近底层海流则不同,变化周期约为51 d。     

基于海面多普勒谱模型的海浪反演精度分析

Microwave remote sensing is one of the most useful methods for observing the ocean parameters. The Doppler frequency or interferometric phase of the radar echoes can be used for an ocean surface current speed retrieval,which is widely used in spaceborne and airborne radars. While the effect of the ocean currents and waves is interactional. It is impossible to retrieve the ocean surface current speed from Doppler frequency shift directly. In order to study the relationship between the ocean surface current speed and the Doppler frequency shift, a numerical ocean surface Doppler spectrum model is established and validated with a reference. The input parameters of ocean Doppler spectrum include an ocean wave elevation model, a directional distribution function, and wind speed and direction. The suitable ocean wave elevation spectrum and the directional distribution function are selected by comparing the ocean Doppler spectrum in C band with an empirical geophysical model function(CDOP). What is more, the error sensitivities of ocean surface current speed to the wind speed and direction are analyzed. All these simulations are in Ku band. The simulation results show that the ocean surface current speed error is sensitive to the wind speed and direction errors. With VV polarization, the ocean surface current speed error is about 0.15 m/s when the wind speed error is 2 m/s, and the ocean surface current speed error is smaller than 0.3 m/s when the wind direction error is within 20° in the cross wind direction.     

海洋对热带气旋响应的一种改进模式

建立一个改进的二层非线性原始方程海洋模式，研究海洋对热带气旋的响应。采用湍流动能收支参数化风应力产生的垂直混合（夹卷），其中考虑了盐度对层结强度的影响。通过海洋对７００２号台风响应的数值模拟，结果表明，在引起海表温度下降的各热通量分量中，夹卷约占了８３％，余下的海表面热通量占了１７％。在台风路径转向的右侧，海洋出现强烈的降温，表现出明显的右偏性。降温的幅度、范围和形状均与观测结果较为一致。