A new hybrid vertical coordinate ocean model and its application in the simulation of the Changjiang diluted water

Based on the analysis of the advantages and disadvantages of some vertical coordinates applied in the calculation of the Changjiang diluted water(CDW),a new hybrid vertical coordinate is designed,which uses σ coordinate for current and σ-z coordinate for salinity.To combine the current and salinity,the Eulerian-Lagrangian method is used for the salinity calculation,and the baroclinic pressure gradient(BPG) is calculated on the salinity sited layers.The new hybrid vertical coordinate is introduced to the widely used model of POM(Princeton Ocean Model) to make a new model of POM-σ-z.The BPG calculations of an ideal case show that POM-σ-z model brings smaller error than POM model does.The simulations of CDW also show that POM-σ-z model is better than POM model on simulating the salinity and its front.     

Introduction of accumulation-dispersal coefficient of marine organisms

Avoiding the subject for fish accumulation, the traditional view in fish population dynamics has ascribed immigration and emigration of fish to dispersal of fish. The main purpose of this paper is to find a quantity that represents the time rate of accumulation-dispersal of marine organisms, and also has some relation to the horizontal convergence of current velocity of the surrounding water. For this, the accumulation-dispersal coefficient is introduced not in the form of diffusion, but in the same form as the convergence. Since the accumulation-dispersal of organisms is a factor that changes its distribution density, all factors causing the change are first classified to locate the position occupied by the accumulation dispersal. Each factor corresponds to each coefficient appearing in a linearized equation describing the rate of change in the density, averaged over a region or a group. The immigration-emigration coefficient is divided into three coefficients of passage, accumulation dispersal and diffusion velocity. For the organisms ranging in a nearly horizontal layer, the accumulation-dispersal coefficient is shown to equal the area-averaged horizontal convergence of organismal velocity relative to land, which is the sum of the area-averaged horizontal convergences of swimming velocity relative to water and of current velocity. However, the area-averaged convergence of current velocity associated with the accumulation-dispersal coefficient for a region is shown to be somewhat different from the usual one.     

Modal evolution and inversion for seabed geoacoustic properties inweakly range-dependent shallow-water waveguides

In a shallow-water ocean environment, the range dependent variation of the geoacoustic properties of the seabed is one of the crucial factors affecting sound propagation. Since the local modes of propagation depend on the spatial changes in the bottom sediments, the local eigenvalues of these modes are useful as tools for examining the range dependence of the sediment properties. In order to extract the local eigenvalues from measurements of the pressure field in a laterally inhomogeneous waveguide, the zeroth-order asymptotic Hankel transform with a short sliding window is utilized. The local peak positions in the output spectra differ from the local eigenvalues due to both the range variation of the local modes and the interference of adjacent modes. The departure due to the former factor is evaluated analytically by using the stationary phase method. In order to reduce the error induced by the latter factor, mode filtering is utilized by incorporating data from a fixed vertical array of receivers. The methods developed are applied to simulated pressure field data as well as experimental field data, and it is shown that the range evolution of the local modes can be successfully estimated. In addition, field measurements are used to demonstrate that the modal trajectories in range can be used to infer the range-dependent geoacoustic properties of the seabed     

Numerical study on residual current and its impact on mass transport in the Hangzhou Bay and the Changjiang Estuary Ⅰ.A3-D joint model of the Hangzhou Bay and the Changjiang Estuary

A 3-D numerical model is set up in a large domain covering the Hang zhou Bay and the Changjiang Estuary based on the ECOM model in orthogonal curvilinear coordinates.The numerical schemes for baroclinic pressure gradient (BPG) terms and convective terms are improved in the paper according to the characteristics of velocity field and mass transport in the area.The model is validated by the simulations of residual current and salinity transport in the Hangzhou Bay and the Changjiang Estuary.     

Numerical study on residual current and its impact on mass transport in the Hangzhou Bay and the Changjiang Estuary I. A3-D joint model of the Hangzhou Bay and the Changjiang Estuary

ＩＮＴＲＯＤＵＣＴＩＯＮＲｅｓｉｄｕａｌｃｕｒｒｅｎｔａｎｄｉｔｓｉｍｐａｃｔｏｎｍａｓｓｔｒａｎｓｐｏｒｔａｒｅｉｍｐｏｒｔａｎｔｔｏｔｈｅｓｔｕｄｙｏｆｃｏａｓｔａｌｅｎ ｖｉｒｏｎｍｅｎｔ.Ａｌｔｈｏｕｇｈｌｏｔｓｏｆｒｅｓｅａｒｃｈｅｓｈａｖｅｂｅｅｎｄｏｎｅｏｎｔｈｅｒｅｓｉｄｕａｌｃｕｒｒｅｎｔａｎｄｍａｓｓｔｒａｎｓ ｐｏｒｔｉｎｔｈｅＨａｎｇｚｈｏｕＢａｙａｎｄｔｈｅＣｈａｎｇｊｉａｎｇＥｓｔｕａｒｙ (Ｃａｏ ,1 989;ＣａｏａｎｄＦａｎｇ ,1 986 ;Ｃｈｅｎｅｔａｌ.,1 992 ;ＨｕａｎｄＨ…     

Hybrid N-order Lagrangian Interpolation Eulerian?Lagrangian Method for Salinity Calculation

The Eulerian?Lagrangian method (ELM) has been used by many ocean models as the solution of the advection equation, but the numerical error caused by interpolation imposes restriction on its accuracy. In the present study, hybrid N-order Lagrangian interpolation ELM (LiELM) is put forward in which the N-order Lagrangian interpolation is used at first, then the lower order Lagrangian interpolation is applied in the points where the interpolation results are abnormally higher or lower. The calculation results of a step-shaped salinity advection model are analyzed, which show that higher order (N=3?8) LiELM can reduce the mean numerical error of salinity calculation, but the numerical oscillation error is still significant. Even number order LiELM makes larger numerical oscillation error than its adjacent odd number order LiELM. Hybrid N-order LiELM can remove numerical oscillation, and it significantly reduces the mean numerical error when N is even and the current is in fixed direction, while it makes less effect on mean numerical error when N is odd or the current direction changes periodically. Hybrid odd number order LiELM makes less mean numerical error than its adjacent even number order LiELM when the current is in the fixed direction, while the mean numerical error decreases as N increases when the current direction changes periodically, so odd number of N may be better for application. Among various types of Hybrid N-order LiELM, the scheme reducing N-order directly to 1st-order may be the optimal for synthetic selection of accuracy and computational efficiency.     

A numerical model on sea surface wind of typhoon and its hindcasting calibration

Ａｎｕｍｅｒｉｃａｌｍｏｄｅｌｏｎｓｅａｓｕｒｆａｃｅｗｉｎｄｏｆｔｙｐｈｏｏｎａｎｄｉｔｓｈｉｎｄｃａｓｔｉｎｇｃａｌｉｂｒａｔｉｏｎ￥ＳｈｅＪｕｎ；ＹｕａｎＹｅｌｉａｎｄＰａｎＺｅｎｇｄｉ（ＲｅｃｅｉｖｅｄＡｐｒｉｌ１４，１９９３；ａｃｃｅｐｔｅ...     

Lateral stratification and dynamic balance at the Matane transect in the lower Saint Lawrence estuary

The lateral stratification and momentum balance in the lower Saint Lawrence estuary (LSLE) were investigated over a period of six days from current, temperature, and salinity measurements taken in 1979 at the Matane transect. Asymptotic singular decomposition analysis was used to show that the spatial variability of hydrodynamic parameters at that cross-section was more pronounced than their temporal or tidal variability. Classification of the section in the lateral circulation/stratification diagram revealed that the LSLE is one of the most laterally-stratified estuaries examined in the world. Estimates of various acceleration terms in the lateral dynamic balance equation showed that, when averaged over many tidal cycles, the motion can be considered as being quasigeostrophic to a first approximation, with contributions of the same order of magnitude from the baroclinic and the barotropic pressure gradients. Convection accelerations were found to be important only near the north shore, where strong cross-channel currents were recorded, whilst local and centrifugal accelerations were found to be of secondary importance.     

A note on average cross-isopycnal mixing in the North Atlantic ocean

Area averaged mixing is inferred from the difference between cross-isopycnal advection and air-sea forcing, using a density budget in isopycnal layers. The results suggest that mixing is weak (indistinguishable from zero) in the thermocline in qualitative agreement with local and regional measurements. Mixing increases in warmer density classes, but becomes indeterminate at low densities, where climatological flux errors are largest. At high density, mixing reappears owing to the heat flux required to balance the formation of dense water and reflects mixing with intermediate and bottom water.     

Rapid changes of the oceanic circulation in a hierarchy of ocean models

The response time of the large-scale oceanic circulation due to freshwater perturbations is investigated with models of different complexity. A three-dimensional ocean general circulation model (OGCM) and a zonally averaged ocean model are employed. In order to distinguish advection and diffusion from fast baroclinic processes (e.g. waves in the OGCM) a color tracer is injected at the same time and location as the freshwater discharges. In spite of the inability of the zonally averaged model to represent wave processes in a realistic way similarities with the OGCM are found for the propagation patterns of density anomalies and of color tracer. In the OGCM as well as in the zonally averaged model, density anomalies propagate faster than anomalies of a passive color tracer in the case of vertical density stratification. The progression of density anomalies leads to changes of the oceanic circulation, and both oceanic models exhibit circulation changes in areas distant from the discharge places long before the passively entrained color tracer has reached these regions. The fact that a zonally averaged model simulates baroclinic processes faster than advection even if internal gravity waves are not represented due to neglected acceleration terms, is clarified with a conceptual box model.     

多目标渔业资源科学调查的分层方案设计优化

渔业资源科学调查常用于收集高质量的生物学和生态学数据以进行渔业资源评估与管理。渔业资源调查设计优化十分必要,其有助于在经济有效的采样努力量下提高调查估计量的精度。对于估计单鱼种资源量指数和物种多样性指数为主的多目标渔业资源调查设计,本研究应用模拟方法评价和优化了分层随机采样设计的分层方案。在不同月份,对于不同调查目标,分别比较了不同分层方案设计的表现。对于大多数指标,与简单随机采样设计相比,分层方案设计可以提高调查估计量的精度。目前采用的具有5层的分层随机采样设计表现最好。通过分层方案设计可以补偿由于采样努力量降低造成的估计量精度的下降,采样努力量的减少有助于降低调查成本、减轻调查拖网对于种群数量较低鱼种的不利影响。本研究表明对于不同的调查目标,最优化的分层方案设计不同。调查后分析有助于改善渔业资源调查的分层方案设计。     

Hybrid <Emphasis Type="Italic">N</Emphasis>-order Lagrangian interpolation Eulerian-Lagrangian method for salinity calculation

The Eulerian-Lagrangian method (ELM) has been used by many ocean models as the solution of the advection equation, but the numerical error caused by interpolation imposes restriction on its accuracy. In the present study, hybrid N-order Lagrangian interpolation ELM (LiELM) is put forward in which the N-order Lagrangian interpolation is used at first, then the lower order Lagrangian interpolation is applied in the points where the interpolation results are abnormally higher or lower. The calculation results of a step-shaped salinity advection model are analyzed, which show that higher order (N=3-8) LiELM can reduce the mean numerical error of salinity calculation, but the numerical oscillation error is still significant. Even number order LiELM makes larger numerical oscillation error than its adjacent odd number order LiELM. Hybrid N-order LiELM can remove numerical oscillation, and it significantly reduces the mean numerical error when N is even and the current is in fixed direction, while it makes less effect on mean numerical error when N is odd or the current direction changes periodically. Hybrid odd number order LiELM makes less mean numerical error than its adjacent even number order LiELM when the current is in the fixed direction, while the mean numerical error decreases as N increases when the current direction changes periodically, so odd number of N may be better for application. Among various types of Hybrid N-order LiELM, the scheme reducing N-order directly to 1st-order may be the optimal for synthetic selection of accuracy and computational efficiency.     

基于海空磁力测量数据的向下延拓方法比较

地磁场的向下延拓是实现磁场数据空间转换的主要手段,是多源海洋磁力测量数据融合和构建三维海洋磁空间背景场模型的关键技术。在分析频域位场向下延拓方法原理基础上,分别采用了4种向下延拓方法将300 m高度磁测数据延拓至200 m高度,并与同一海区200 m高度实测数据对比,以检验各延拓方法的计算效果,为构建三维海洋磁空间背景场模型提供参考。结果表明:Tikhonov正则化法延拓结果与实测数据最接近,积分迭代法容易放大实测数据中高频噪声,延拓误差最大,迭代Tikhonov正则化法与Landweber迭代法延拓结果等值线最光滑,各方法计算误差相差不大,均在4 nT左右,但各方法计算值缺乏局部细节信息,与实测数据仍存在一定差距。     

钱塘江河口尖山河段整治规划线的方案试验研究

根据实测资料分析和二维数值计算 ,借助实体模型 ,通过定床与动床试验 ,比较分析了两种整治规划线方案。结果表明 ,推荐方案对削弱尖山以上江道的潮势 ,改善萧绍平原的排涝条件 ,效果显著     

Comparison of Sea Surface Heights Derived from Satellite Altimetry and from Ocean Bottom Pressure Gauges: The SW Pacific MOTEVAS Project

A bottom pressure gauge (BPG) was installed in proximity (3.7 km at closest approach) of Jason-1 and formerly TOPEX/Poseidon (T/P) ground track No. 238 at the Wusi site, located ∼ 10 km offshore off the west coast of Santo Island, Vanuatu, Southwest (SW) Pacific. Sea level variations are inferred from the bottom pressure, seawater temperature, and salinity, corrected for the measured surface atmospheric pressure. The expansion of the water column (steric increase in sea surface height, SSH) due to temperature and salinity changes is approximated by the equation of state. We compare time series of SSH derived from T/P Side B altimeter Geophysical Data Records (GDR) and Jason-1 Interim Geophysical Data Records (IGDR), with the gauge-inferred sea level variations. Since altimeter SSH is a geocentric measurement, whereas the gauge-inferred observation is a relative sea level measurement, SSH comparison is conducted with the means of both series removed in this study. In addition, high-rate (1-Hz) bottom pressure implied wave heights (H_1/3) are compared with the significant wave height (SWH) measured by Jason-1. Noticeable discrepancy is found in this comparison for high waves, however the differences do not contribute significantly to the difference in sea level variations observed between the altimeter and the pressure gauge. In situ atmospheric pressure measurements are also used to verify the inverse barometer (IB) and the dry troposphere corrections (DTC) used in the Jason IGDR. We observe a bias between the IGDR corrections and those derived from the local sensors. Standard deviations of the sea level differences between T/P and BPG is 52 mm and is 48 mm between Jason and BPG, indicating that both altimeters have similar performance at the Wusi site and that it is feasible to conduct long-term monitoring of altimetry at such a site.     

The annual cycle of vertical mixing and restratification in the Northern Gulf of Eilat/Aqaba (Red Sea) based on high temporal and vertical resolution observations

The stratification in the Northern Gulf of Eilat/Aqaba follows a well-known annual cycle of well-mixed conditions in winter, surface warming in spring and summer, maximum vertical temperature gradient in late summer, and erosion of stratification in fall. The strength and structure of the stratification influences the diverse coral reef ecosystem and also affects the strength of the semi-diurnal tidal currents. Long-term (13 months) moored thermistor data, combined with high temporal and vertical resolution density profiles in deep water, show that transitions from summer to fall and winter to spring/summer occur in unpredictable, pulses and are not slow and gradual, as previously deduced from monthly hydrographic measurements and numerical simulations forced by monthly climatologies. The cooling and deepening of the surface layer in fall is marked by a transition to large amplitude, semi-diurnal isotherm displacements in the stratified intermediate layer. Stratification is rebuilt in spring and summer by intermittent pulses of warm, buoyant water that can increase the upper 100–150 m by 2 °C that force surface waters down 100–150 m over a matter of days. The stratification also varies in response to short-lived eddies and diurnal motions during winter. Thus, the variability in the stratification exhibits strong depth and seasonal dependence and occurs over range of timescales: from tidal to seasonal. We show that monthly or weekly single-cast hydrographic data under-samples the variability of the stratification in the Gulf and we estimate the error associated with single-cast assessments of the stratification.     

Evaluation of Seasonal Sea Level Variation at Syowa Station, Antarctica, Using GPS Observations

Global Positioning System (GPS) measurements of the vertical displacement of fast ice near Syowa Station, Antarctica, were conducted between April and December 1998 to evaluate measurements of sea level variation derived with a conventional bottom pressure gauge (BPG). The BPG-derived sea level revealed a seasonal variation of about 0.13 m, with a high in April–June and a low in November–December. The GPS-derived sea level, combined with observed sea ice thickness, supported the BPG result, with an RMS error of 0.007 m. Our result also demonstrates that GPS is a powerful technique for monitoring sea level variations. This revised version was published online in July 2006 with corrections to the Cover Date.     

TY01和O02两种海面空气动力粗糙度方案的比较

本文采用美国国家浮标站(44008)2003年1～3月的资料，通过COARE算法(版本2．6b)，比较了O02和TY01这两种海面空气动力粗糙度长度的参数化方案。通过对摩擦速度、拖曳系数、海面粗糙度及风应力等物理因子的计算得出：在粗糙的海面上，TY01和O02两种参数化方案的计算结果是比较一致的。在考虑浪的信息方面，TY01和O02都是很好的参数化方案．它们都可以适用于不同的风速条件，适用于各种尺度的海洋及湖泊。但是这两种方案在处理幼波时存在不连续的缺点。并且，对于风速较小的光滑海面，尽管它们计算的结果较一致，但是仍然存在偏差。据此，本文的结果对于理论分析和数值计算如何正确使用上述两种海面空气动力粗糙度参数化方案，可提供必要的参考价值。     

基于机器学习的内孤立波波要素关系研究

内孤立波在海洋中的传播会携带能量和动量,不同振幅的内孤立波对海洋中的能量交换及海上工程等影响也不同,因此,研究内孤立波振幅与半波宽度、水深、分层条件、密度等水文特征参量之间的关系显得尤为重要。以往在研究中建立内孤立波振幅与它们之间的关系时,会受到不同理论有效适用范围的限制。本文借助实验室的水槽方法,设计了不同的水深、分层及密度条件下的内孤立波系列综合实验,发现内孤立波的振幅与半波宽度、水深、分层条件以及水体密度等参量之间并非简单线性关系。因此,利用机器学习的方法建立内孤立波振幅与上述参量之间的非线性关系,建立了支持向量机（SVM）和随机森林（RF）两种机器学习模型。将1 266组实验数据建立样本库,其中包含训练集970组,测试集296组,对模型进行参数调优,最终通过测试集验证,SVM模型的平均相对误差为17.3%,RF模型的平均相对误差为15.5%。该方法适用于多种不同的水文条件,有效解决先前理论存在的适用性问题。     

Comparison of Sea Surface Heights Derived from Satellite Altimetry and from Ocean Bottom Pressure Gauges: The SW Pacific MOTEVAS Project

A bottom pressure gauge (BPG) was installed in proximity (3.7 km at closest approach) of Jason-1 and formerly TOPEX/Poseidon (T/P) ground track No. 238 at the Wusi site, located ～ 10 km offshore off the west coast of Santo Island, Vanuatu, Southwest (SW) Pacific. Sea level variations are inferred from the bottom pressure, seawater temperature, and salinity, corrected for the measured surface atmospheric pressure. The expansion of the water column (steric increase in sea surface height, SSH) due to temperature and salinity changes is approximated by the equation of state. We compare time series of SSH derived from T/P Side B altimeter Geophysical Data Records (GDR) and Jason-1 Interim Geophysical Data Records (IGDR), with the gauge-inferred sea level variations. Since altimeter SSH is a geocentric measurement, whereas the gauge-inferred observation is a relative sea level measurement, SSH comparison is conducted with the means of both series removed in this study. In addition, high-rate (1-Hz) bottom pressure implied wave heights (H _1/3 ) are compared with the significant wave height (SWH) measured by Jason-1. Noticeable discrepancy is found in this comparison for high waves, however the differences do not contribute significantly to the difference in sea level variations observed between the altimeter and the pressure gauge. In situ atmospheric pressure measurements are also used to verify the inverse barometer (IB) and the dry troposphere corrections (DTC) used in the Jason IGDR. We observe a bias between the IGDR corrections and those derived from the local sensors. Standard deviations of the sea level differences between T/P and BPG is 52 mm and is 48 mm between Jason and BPG, indicating that both altimeters have similar performance at the Wusi site and that it is feasible to conduct long-term monitoring of altimetry at such a site.