Three-dimensional numerical simulation of obfique wave action on vertical walls

A three-dimensional numerical model based on the potential theory was developed to study the oblique wave action on vertical walls. A source term inside the domain was used to generate incident waves and outgoing waves were dissipated by sponge layers and transmitted by radiation boundaries. The finite difference method was used to solve the governing equations and boundary conditions in the regular transformed domain in σ-coordinate. Satisfactory agreements between the numerical predictions and experimental results of wave force were obtained. It is concluded that the maximum wave force acting on the vertical walls is induced by the obliquely incident waves rather than the normally incident waves.     

利用归一化层速度预测和研究储层物性的方法

归一化层速度是对任意深度处的层速度进行归一化处理的一种速度形式，归一化处理的目的是把任意深度处的层速度转换为与深度无关的格式，归一化后层速度的横向变化主要取决于地层岩性，孔隙度等物性的变化，利用这一特点可以方便地利用已知钻井资料对其进行标定，比较简单地进行地层岩性，孔隙度等特性的预测。实际应用表明，其预测结果与钻井揭示吻合较好，预测效果令人满意。     

具有内波特性的渤海南部次重力波

次重力波(infragravity waves)是频率范围介于局地平均Brunt-V¾is¾l¾频率和0.05Hz(涌浪频率)之间的一种海洋波动,其能量在深水大洋中一般甚低,但在浅海中却很强.该波动的研究对于浅海海洋内部混合、近岸海洋工程、海洋生态学、海洋沉积学及军事海洋学均有重要的意义.Munk[1]很早就开始研究这类波动,然而主要由于观测技术的限制,对于该类波动的动力学机制的研究进展迟缓.     

辽东湾地区下第三系地震速度—岩性预测模型研究

地震速度－岩性预测模型由四个子模型组成：即（１）砂泥岩压实模型；（２）地震层速度转换模型；（３）速度校正模型；（４）砂岩指数转换模型。辽东湾地区的实际资料证明，上述四个子模型不但决定了岩性预测的具体方法，而且严格控制岩性预测的精度。     

Wave Defending Effects of V-Type Bottom-Mounted Breakwaters

An analytical method is developed for the study of the wave defending effects of the V-type bottom-mounted breakwater. The breakwater is assumed to be rigid, thin, impermeable and vertically located in water of constant depth. The fluid domain is divided into three sub-regions by an imaginary interface. The velocity potential in each region is expanded by eigenfunctions. By satisfying the corresponding boundary conditions and matching conditions in and between sub-regions, a set of hnear algebraic equations can be obtained to determine the unknown coetfficients for the eigenfunction expansions for each sub-region. The accuracy of the present model is verified by a comparison with existing results for the case of an isolated breakwater. Numerical results, in the form of contour maps of the relative wave amplitude around the breakwater, are presented for a range of wave and breakwater parameters. The results show that the V-type bottommounted breakwater is generally effective in defending against waves. In general, the wave height in the protected area is about 20-50 percent of the incident wave height.     

Elucidating the dynamics and mixing agents of a shallow fjord through age tracer modelling

The mixing agents and their role in the dynamics of a shallow fjord are elucidated through an Eulerian implementation of artificial tracers in a three-dimensional hydrodynamic model. The time scales of vertical mixing in this shallow estuary are short, and the artificial tracers are utilized in order to reveal information not detectable in the temperature or salinity fields. The fjord's response to external forcing is investigated through a series of model experiments in which we quantify vertical mixing, transport time scales of fresh water runoff and estuarine circulation in relation to external forcing.Using age tracers released at surface and bottom, we quantify the time scales of downward mixing of surface water and upward mixing of bottom water. Wind is shown to be the major agent for vertical mixing at nearly all depth levels in the fjord, whereas the tide or external sea level forcing is a minor agent and only occasionally more important just close to the bottom. The time scale of vertical mixing of surface water to the bottom or ventilation time scale of bottom water is estimated to be in the range 0.7 h to 9.0 days, with an average age of 2.7 days for the year 2004.The fjord receives fresh water from two streams entering the innermost part of the fjord, and the distribution and age of this water are studied using both ageing and conservative tracers. The salinity variations outside this fjord are large, and in contrast to the salinity, the artificial tracers provide a straight forward analysis of river water content. The ageing tracer is used to estimate transport time scales of river water (i.e. the time elapsed since the water left the river mouth). In May 2004, the typical age of river water leaving the fjord mouth is 5 days. As the major vertical mixing agent is wind, it controls the estuarine circulation and export of river water. When the wind stress is set to zero, the vertical mixing is reduced and the vertical salinity stratification is increased, and the river water can be effectively exported out of the fjord.We also analyse the river tracer fields and salinity field in relation to along estuary winds in order to detect signs of wind-induced straining of the along estuary density gradient. We find that events of down estuary winds are primarily associated with a reduced along estuary salinity gradient due to increased surface salinity in the innermost part of the fjord, and with an overall decrease in vertical stratification and river water content at the surface. Thus, our results show no apparent signs of wind-induced straining in this shallow fjord but instead they indicate increased levels of vertical mixing or upwelling during down estuary wind events.     

Air-Sea Gas Transfer Velocity in Stormy Winter Estimated from Radon Deficiency

Thirteen vertical profiles of ²²⁶Ra and ²²²Rn in the near-surface water were obtained in the western North Pacific in winter, and the gas transfer velocities across the air-sea interface were estimated. The transfer velocities found by applying a steady state model varied widely from 2.1 to 30.2 m day⁻¹ with a mean of 9.4 m day⁻¹. The mean value is almost 5 times higher than that in summer in other oceans, and the maximum value is a record high for world oceans. This is partly due to the inadequacy of the steady state model, which overestimates when stronger winds blow in more recent days than the ²²²Rn half-life of about 4 days. In fact, a strong low pressure zone passed through the station about 2 days earlier, which was one of the low pressure zones that with a period of develop once a week or so in the northwestern North Pacific in winter. Instead of steady-state removal, if half of the radon removal occurred sporadically every 7 days, and the last removal took place two days before the observation, the transfer velocity would be 26 m day⁻¹. Our mean transfer velocity, which is less than 20% different from the steady state value including both overestimated and underestimated values, 9.4 ± 4.8 m day⁻¹, seems to represent the mean state of this region in winter. This suggests that the gas exchange fluxes under extremely rough conditions in the open ocean are larger than those estimated by using a transfer velocity equation with a linear or quadratic relationship with wind speed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Correction method for full-depth current velocity with lowered acoustic Doppler current profiler (LADCP)

A new method is presented to process and correct full-depth current velocity data obtained from a lowered acoustic Doppler current profiler (LADCP). The analysis shows that, except near the surface, the echo intensity of a reflected sound pulse is closely correlated with the magnitude of the difference in vertical shear of velocity between downcast and upcast, indicating an error in velocity shear. The present method features the use of echo intensity for the correction of velocity shear. The correction values are determined as to fit LADCP velocity to shipboard ADCP (SADCP) and LADCP bottom-tracked velocities. The method is as follows. Initially, a profile of velocity relative to the sea surface is obtained by integrating vertical shears of velocity after low-quality data are rejected. Second, the relative velocity is fitted to the velocity at 100–800 dbar measured by SADCP to obtain an “absolute” velocity profile. Third, the velocity shear is corrected using the relationship between the errors in velocity shears and echo intensity, in order to adjust the velocity at sea bottom to the bottom-tracked velocity measured by LADCP. Finally, the velocity profile is obtained from the SADCP-fitted velocity at depths less than 800 dbar and the corrected velocity shear at depths greater than 800 dbar. This method is valid for a full-depth LADCP cast throughout which the echo intensity is relatively high (greater than 75 dB in the present analysis). Although the processed velocity may include errors of 1–2 cm s⁻¹, this method produced qualitatively good current structures in the Northeast Pacific Basin that were consistent with the deep current structures inferred from silicate distribution, and the averaged velocities were significantly different from those calculated by the Visbeck (2002) method.     

Temporal change of clustered distribution of planktonic ciliates in Toyama Bay in summers of 1989 and 1990

Temporal change of clustered distribution in vertical profiles of three nutritional groups of planktonic ciliates, e.g. heterotrophic naked ciliates, mixotrophic naked ciliates and heterotrophic loricated ciliates, was investigated by following a drifting buoy in Toyama Bay on the Japan Sea coast of central Japan in summers of 1989 and 1990. Clustered distribution, represented as the mode of population density in the vertical plane, occurred mainly in the oligotrophic upper layer (0–50 m depth) above the subsurface chlorophyll-a maximum layer. Its clustered degree was stronger when the mode of population density in the vertical plane was formed at shallower depth, while its longevity was shorter as mentioned above. Vertical distribution of ciliates during summer in Toyama Bay is characterized by ephemeral clustered distribution, or in other wards, by rapid alternations of appearance and disappearance of the clustered distribution.     

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

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