基于新型高频声学多普勒流速剖面仪的河口近底部边界层观测初探

河口物质输运、能量交换与底边界层内的水动力过程密切相关,底边界层参数(如切应力、拖曳系数)的确定至关重要。挪威Nortek公司生产的新型声学多普勒流速剖面仪AD2CP相比传统ADCP具有高频、低噪的优点,可用于高频(16Hz)流速剖面观测,而被广泛应用于底边界层观测的ADV只能测量单点高频流速。本文采用AD2CP在长江口南槽最大浑浊带区域进行座底式观测,并与同步近底部三脚架上ADV的观测结果进行对比。结果表明,使用AD2CP测得的近底部平均流速与ADV的测量结果吻合良好;使用惯性耗散法计算了底切应力,基于ADV的单点高频流速数据计算结果为2.16×10~(-2)～5.69×10~(-1)N/m~2,基于AD2CP的结果为2.09×10~(-2)～4.26×10~(-1)N/m~2,二者范围大致相当。在此基础上,基于AD2CP数据计算出摩阻流速为4.55×10~(-3)～2.06×10~(-2)m/s、底拖曳系数范围为1.84×10~(-4)～2.49×10~(-3),与ADV的计算结果基本一致。此外,由于AD2CP可以获得高频的流速剖面数据,优于单点ADV,具备观测近底部边界层参数和边界层内湍流剖面的潜力。     

高、低潮数据的调和分析

本文给出了由高、低潮数据计算潮汐调和常数和由最大流速及转流数据计算潮流调和常数的方法。文章还讨论了不同频率分潮间的混淆效应,指出本方法用于混合型和全日型潮汐可获得比半日型潮汐好的效果。对两个港口实测高、低潮数据分析表明,所得调和常数与由逐时潮高分析所得的数值一般较接近,但对长周期分潮较差。本方法所得调和常数用于预报高、低潮或最大流速和转流时可得到很好的结果。     

非线性弱色散波内部流场的重构

基于势流理论和级数直接求逆方法,本文建立了基于Ｂｏｕｓｉｎｅｓｑ方程或Ｇｒｅｅｎ－Ｎａｇｈｄｉ方程给出的水深平均流速或某特征流速及波面信息重构非线性弱色散波内部流场的算法。以Ｂｏｕｓｉｎｅｓｑ方程的孤立波解为例,用本反演方法计算了孤立波的表面水平流速及底部水平流速。结果表明本算法是有效的。本反演算法可用于获取非线性弱色散波的内部流场的详细信息。     

斜压海洋动力学的一种三维数值模式 Ⅰ.动力学方程数值格式

依据自由海面海洋动力学原始方程建立了一种三维有限差分数值模式,可用于潮波、风暴潮和海流的数值模拟和预报。运动方程和连续方程的数值格式采用内、外模态分离的技术。外模态采用交替方向隐格式,用于计算海面高度和垂直平均流速,时间步长不受Ｃｏｕｒａｎｔ－Ｆｒｉｄｅｒｉｃｈｓ－Ｌｅｗｙ条件限制;内模态采用半隐格式,用于计算海流的垂直２颁布,其时间步长可大于外模态时间步长。模式的计算程度比一般显式模式可快１０倍     

河口地区航槽开挖后槽内流速变化

根据河口挖槽的特点进行水槽与挖槽设计，并在宽水槽试验的基础上，通过因次和线型分析，得出顺直挖槽和斜挖槽后统一的流速计算公式，公式计算值与实测值为接近，本文公式已用于实际航道回淤预报计算中。     

导管架浮筒浮运强度分析

本文介绍了钢质固定平台导管架浮运时总强度分析的准静力法,并以渤海埕北油田A平台导管架的计算为例说明了该方法的应用。准静力法虽然并没有考虑导管架在实际海况中所受到的载荷的随机特性,但其计算简单,在工程实际中仍被广泛的采用。文中介绍的方法及其相应的计算程序可用于中、小型导管架浮运的强度分析。     

潮流作用下洋山港水域悬沙和底沙的交换

通过对泥沙起动流速、淤积流速和不淤流速的计算,划定了悬沙和底沙交换的各个时段,提出了冲刷悬沙浓度、淤积悬沙浓度和不淤悬沙浓度的概念及其计算方法;继而导出了潮流作用下悬沙和底沙在交换层面上单位面积的日交换量计算式,并用于洋山港水域的冲淤计算。结果表明:自然状况下洋山港水域北部年淤积量为37.5 kg/m2,淤高了1.77 cm/a,属微淤;中部年冲刷量为10.0 kg/m2左右,刷深不超过0.52 cm/a,可认为冲淤平衡;南部年冲刷量为55.5 kg/m2,刷深了2.60 cm/a,属微冲。这一结果与该水域海床的自然冲淤状况相吻合。     

声相关海流剖面(ACCP)测量的流速估值方法

声相关海流剖面（ACCP）测量技术对水体回波信号进行空间相关处理来估计流速值，与声学多普勒海流剖面仪（ADCP）相比ACCP的工作频工，更适合深海测流应用，ACCP测量系统对回波信号进行时空相关处理，根据相关函数的位置来估计被测水团的流速值，ACCP从信号中提取流速值要运用参数估计技术，本文介绍了采用极大似然法和最小二乘法进行流速参数的处理方法，给出一些典型流速情形的数值计算结果。     

频率低于惯性频率的潮流之垂直结构

Sverdrup在研究西伯利亚大陆架潮汐时,得到了有关半日潮流垂直结构的重要结论:最大流速随着深度的增加而减小;最大流速方向随深度增加而右偏;最大流速发生时刻则随深度增加而提前。长时间以来,这些结论被人们用以描述海洋中潮流垂直分布的特征。后来,有不少研究者进一步做了研究,得出了一些不同的结果。例如,陈宗镛的计算表明:在渤海中20多米水深     

山东荣成湾月湖口门落潮干道的推移质输运

改进了Gao等 ( 1 994)方法 ,并用以计算潮汐汊道口门落潮干道的垂线平均流速。利用月湖潮汐汊道系统冬、夏季各一个月的潮位资料 ,计算落潮干道的潮流流速和推移质输运率。结果表明 ,流速的计算值和实测值呈良好的相关性 ,落潮干道内落潮流的输沙能力高于涨潮流输沙能力     

Rhines效应对高度计观测的海洋中尺度涡空间分布特征的影响

Rhines效应是指Rossby波和大湍流（中尺度涡）相互作用,将涡动能量以波的形式传播出去,从而使中尺度涡发生形变,最终消亡的一种动力学机制。本文通过比较海洋里涡特征速度和Rossby长波波速的方法,研究了一种广义的Rhines效应对高度计观测的海洋中尺度涡空间分布特征的影响。结果显示,广义Rhines效应比只考虑行星涡度梯度的传统形式对中尺度涡的分布具有更显著的影响。大部分中尺度涡分布在涡特征速度（U_e）大于由广义Beta值计算的Rossby长波波速（U_cg）的区域。这些涡可以由动能反向串级过程获取能量,成长为振幅和空间尺度较大的涡。热带海域以外的“涡旋沙漠”区域,中尺度涡的数量稀少,强度很弱,大都分布于U_ecg的海域。广义Rhines效应可能是这些海域中尺度涡难以成长的动力学机制。     

石岛地震台远震记录反演研究

利用石岛地震台的远震体波记录,采用旋转相关函数法和接收函数法分别反演了台站下方介质的各向异性特征和速度结构.(1)对震中距25°~35°且记录良好的5次地震的ScS震相,采用旋转相关函数法反演了岩石圈的剪切波分裂参数.对深源地震的反演结果表明,石岛地震台快波偏振方向为N94°E,这意味着西沙附近处于近东西向微偏南的拉张或地壳下方的地幔流方向为近东西微偏南,西沙地区地壳是过渡性的,其底部的驱动力主要来自与欧亚板块运动一致的物质流.快慢波时间延迟为1.3 s,估算各向异性层厚度为100 km左右.(2)对震中距20°~60°的9次远震P波波形三分向记录,采用接收函数法反演了地壳和上地幔的S波速度结构.反演结果表明,石岛地震台下方地壳分为3层:约5 km以上有一速度梯度带,S波速度从1.5 km/s逐渐增加到3.5 km/s,其间有若干小的分层;在5~16 km的平均速度为3.8 km/s左右,其间有若干小的分层;在16.0~26.5 km的速度为3.6 km/s左右,这是一个明显的低速层;莫霍面埋深为26.5 km,莫霍面以下平均速度为4.7 km/s,也有若干小的分层,尤其是在莫霍面之下有一个明显的低速层.根据转换波到时分析和速度剖面左右摆动现象,认为反演结果中的小分层可能是不真实的,但在16.0~26.5 km的低速层的真实程度还是较高的,表明下地壳具有一定的塑性.     

一种改进的基于逆Omega-K算法的海面场景SAR原始数据仿真方法

合成孔径雷达(SAR)海面场景原始数据仿真是研究海洋动力参数(表面波浪、风矢量和洋流)的有效工具。目前海面场景原始数据仿真方法已经基于逆Omega-K算法实现了海洋运动参数的空间变化。但是目前仅仅讨论了正侧视情况下的海面场景仿真,应用范围有限,同时没有考虑Stokes漂流以及Bragg相速度的影响,而这两者都是存在于真实海面的。通常情况下为了反演得到海面流场的二维速度矢量,雷达需要从两个不同的方位方向观察海面的同一个区域,因此这就需要考虑大斜视的雷达波束,同时Stokes漂流和Bragg相速度是SAR海表面流场观测不容忽视的两种运动。本文在不改变原有正侧视逆Omega-K算法的情况下,通过增加重新计算零方位时刻的斜视波束中心位置坐标,并据此确定SAR原始数据在多普勒域的位置来将其扩展到大斜侧视逆Omega-K算法,并通过时域Stokes漂流公式到频域内离散化Stokes漂流公式的推导来加入Stokes漂流,以及根据Bragg散射机制加入了Bragg相速度。仿真结果表明,经过聚焦成像后的SAR图像很好的体现了真实海面波浪场的形状,同时能够很好地反演出设定的雷达径向流场速度,且流速精度误差控制在6%以内。最后也证明了Bragg相速度以及Stokes漂流对于海面流场的影响不可忽视。     

多普勒雷达资料在近海强台风模拟中的同化试验

利用中尺度数值模式(WRF),并同化了多普勒雷达反射率和径向速度资料以及非常规的观测资料,对近几年登陆于浙闽沿海的4例强台风进行了数值模拟。通过高时空分辨率的模拟结果对比分析表明:雷达资料的同化,对近海登陆台风路径和降水模拟以及中尺度降水特征都有进一步改进的效果;模拟较好的揭示了台风近中心螺旋云带中的强中尺度对流系统。通过模拟分析表明,在台风近中心的螺旋云带中,低层有一条强辐合线存在,它与实况多普勒雷达给出的低层平显(PPI)强度回波带有较好的对应关系,也与沿海地区中尺度暴雨系统紧密联系,并由此看到近海海域降水带和强对流区的存在。     

P-Vector inverse method evaluated using the modular ocean model (MOM)

Several major inverse methods (Stommel-Schott method, Wunsch method, and Bernoulli method) have been successfully developed to quantitatively estimate the geostrophic velocity at the reference level from hydrographic data. No matter the different appeance, they are based on the same dynamical sophistication: geostrophy, hydrostatic, and potential density (ϱ) conservation (Davis, 1978). The current inverse methods are all based on two conservation principles: potential density and potential vorticity (q=f∂ϱ/∂z) and require β-turning. Thus, two necessary conditions can be incorporated into any inverse methods: (1) non-coincidence of potential density and potential vorticity surfaces and (2) existence of vertical turning of the velocity (β-turning). This can be done using the P-Vector, a unit vector in the direction of ▽ϱ×▽q (Chu, 1994, 1995). The first necessary condition becomes the existence of the P-vector, and the second necessary condition leads to the existence of the P-vector turning in the water column. Along this line, we developed the P-vector inverse emthod with a pre-requirement check-up. The method was verified in this study using the Modular Ocean Model (MOM) from Pacanowskiet al. (1991) version of Bryan-Cox-Semtner ocean general circulation model (OGCM), which is based on the work of Bryan (1969). The statistically steady solutions of temperature and salinity from MOM are used as a “no-error data” set for computing absolute geostrophic velocities by the P-vector inverse method. Circulations are similar between the MOM statistically steady solutions and the P-vector solutions. Furthermore, the quantitative analysis shows that this inverse method has capability of picking up the major signal of the velocity field.     

0505号台风"海棠"特殊移动路径的分析

“海棠”台风在弱环境场引导气流中，受台湾东南面生成的次生低压的吸引作用，产生了打转效应．台风形成的不对称结构影响其移向和降水分布，路径突转时极易产生短时强降水．食指状和“人”字形的外螺旋雨带、螺旋雨带最大回波强度区的移动与台风前进方向一致．在台风停滞打转期间台风上下层的正负速度中心极值有一个调整过程，螺旋雨带都很接近同心圆，是台风风力最大的时候．多普勒雷达提供的风廓线产品资料能够较好地反映出垂直风场的连续变化，同时也给台风定位工作提供了有价值的参考资料．     

SAR-imaged spiral eddies in Mutsu Bay and their dynamic and kinematic models

A pattern of slick streaks winding into a spiral, known as a spiral eddy, was identified in 5 images taken by the ERS-1/2 synthetic aperture radar (SAR) in Mutsu Bay (Japan); dynamic and kinematic models of these spiral eddies have been proposed. Common characteristics of the five spiral eddies are: 1) an eddy diameter of about 15 km; 2) their location in the western part of the bay; and 3) their cyclonic direction of rotation. Moreover, the wind conditions over the bay were common: prior to acquiring the images, a strong easterly wind continued blowing for more than one day. The wind field on the bay is known to be orographically steered and has strong windstress vorticity, which generates cyclonic circulation. The diameter and location of the circulation simulated with a numerical ocean model corresponded well to those of the identified spiral eddies. Based on these facts, we propose a dynamic model for the movement of a slick streak, and a kinematic model for the formation of a spiral eddy. We have assumed calm air, a microlayer and seawater with a cyclonic circulation in the dynamic model. The balance of forces is established in the microlayer among the frictional force from the seawater, the frictional force from the calm air, the gravitational force, and the Coriolis force. As a result, the velocity vector of the microlayer deflects slightly towards the center of the cyclonic circulation. We have assumed a point source of the microlayer in the kinematic model. The shapes of a slick streak simulated with the models agree well with the identified patterns in the SAR images.     

Estimating hypocentral uncertainties for marine microearthquake surveys: A comparison of the generalized inverse and grid search methods

For microearthquake surveys conducted with small networks in regions where the seismic velocity structure has large vertical gradients, the formal errors accompanying hypocentral solutions obtained by a generalized inverse method may be misleading since they do not incorporate the effects of nonlinearity in travel times. An alternative method for estimating uncertainties involves calculating travel time residuals over a regular grid and using the F statistic to contour confidence volumes. We present a statistical expression for the latter confidence limits that is applicable when an independent estimate of arrival time errors is available from observations accumulated for a number of earthquakes. Synthetic experiments comparing the results of the grid search and generalized inverse methods show that in cases where solutions are obtained either without S wave information or for epicenters which lie well outside the network, the effects of nonlinearity on the shape of the confidence regions may be significant. However, for the well-observed events both methods yield comparable confidence volumes in good agreement with the distribution of hypocenters obtained from repeated locations incorporating random errors. The generalized inverse method has the advantage that it requires fewer calculations, so the examination of systematic errors in hypocentral parameters produced by uncertainties in the seismic velocity structure can be studied in a more computationally efficient manner. Except in the cases of poorly resolved earthquakes, the effects of nonlinearity on uncertainties in hypocentral parameters can be observed by the application of the F statistic to the variation of the generalized inverse travel time residuals with focal depth.     

Bottom water circulation in Cascadia Basin

A combination of beta spiral and minimum length inverse methods, along with a compilation of historical and recent high-resolution CTD data, are used to produce a quantitative estimate of the subthermocline circulation in Cascadia Basin. Flow in the North Pacific Deep Water, from 900-1900 m, is characterized by a basin-scale anticyclonic gyre. Below 2000 m, two water masses are present within the basin interior, distinguished by different potential temperature-salinity lines. These water masses, referred to as Cascadia Basin Bottom Water (CBBW) and Cascadia Basin Deep Water (CBDW), are separated by a transition zone at about 2400 m depth. Below the depth where it freely communicates with the broader North Pacific, Cascadia Basin is renewed by northward flow through deep gaps in the Blanco Fracture Zone that feeds the lower limb of a vertical circulation cell within the CBBW. Lower CBBW gradually warms and returns to the south at lighter density. Isopycnal layer renewal times, based on combined lateral and diapycnal advective fluxes, increase upwards from the bottom. The densest layer, existing in the southeast quadrant of the basin below 2850 m, has an advective flushing time of 0.6 years. The total volume flushing time for the entire CBBW is 2.4 years, corresponding to an average water parcel residence time of 4.7 years. Geothermal heating at the Cascadia Basin seafloor produces a characteristic bottom-intensified temperature anomaly and plays an important role in the conversion of cold bottom water to lighter density within the CBBW. Although covering only about 0.05% of the global seafloor, the combined effects of bottom heat flux and diapycnal mixing within Cascadia Basin provide about 2-3% of the total required global input to the upward branch of the global thermohaline circulation.     

Dynamic parameters in models of atmospheric vortex structures

Vortex simulation and the computation of fields of dynamic parameters of vortex structures (velocity, rotor velocity, and helicity) are carried out with the use of exact hydrodynamic equations in a cylindrical coordinate system. Components of centripetal and Coriolis accelerations are taken into account in the initial equations. Internal and external solutions are defined. Internal solutions ignore the disturbances of the pressure field, but they are considered in external solutions. The simulation is carried out so that the effect of accounting for spatial coordinates on the structure of the above fields is pronounced. It is shown that the initial kinetic energy of rotating motion transforms into the kinetic energy of radial and vertical velocity components in models with centripetal acceleration. In models with Coriolis acceleration, the Rossby effect is clearly pronounced. The method of an “inverse problem” is used for finding external solutions, i.e., reconstruction of the pressure field at specified velocity components. Computations have shown that tangential components mainly contribute to the velocity and helicity vortex moduli at the initial stage.