A phase-averaged Boussinesq model with effective description of carrier wave group and associated long wave evolution

From the phase-resolving improved Boussinesq equations (Beji and Nadaoka, Ocean Engineering 23 (1996) 691), a phase-averaged Boussinesq model for water waves is derived by more effectively describing carrier wave groups and accompanying long wave evolution with less CPU time. Linear shoaling characteristics of carrier wave equations are investigated and found to agree exactly with the analytical expression obtained from the constancy of energy flux for the improved Boussinesq equations themselves, showing that the present model equations are the results of a consistent derivation procedure regarding energy considerations. Numerical simulations of the derived equations for the single wave group and narrow-banded random waves show the validity of the present model and its high performance, especially on the CPU time.     

长江口南支潮流的一个数值模拟

为了能够精确地拟合天然水域不规则的水陆边界,并能保持有限差分数值计算在矩形网格中进行,本文利用Thompson的数值网格生成技术,在长江口南支七丫口至横沙河段设计了一个椭圆型边界拟合坐标系,并在这个坐标系下建立了平面二维潮流数学模型。通过数值求解此数学模型,实现了计算区域内现状潮流的一个数值模拟,模拟结果表明,模式的设计是成功的。     

星载SAR浅海水下地形和水深测量模拟仿真──水下地形高度、坡度和方向与可测水深分析

根据星载合成孔径雷达(SAR)浅海水下地形和水深成像机理,建立了浅海水下地形和水深雷达后向散射截面仿真模型.利用该模型模拟并分析了不同地形条件下,浅海水下地形的雷达后向散射截面.分析结果表明,水下地形高度越高,SAR可测量的水深越深;水下地形坡度越大,越易被SAR所观测.水下地形的星载SAR测量还与水下地形的方向有关,与卫星飞行方向平行的水下地形最易被SAR观测,与卫星飞行方向垂直的水下地形最不易被SAR观测.     

渤海、黄海强温跃层数值预报模式

文中对前人关于水温垂直结构计算、模拟和预报研究成果进行了评价，并在此基础上建立了考虑海面吸收辐射和透射辐射、地形、风、海流、界面摩擦及其引起的混合作用的强温跃层三维数值预报模式。该模式运行的驱动量仅为风和气温场资料，因而具有较好的实用性；此外，从试报的结果看，效果是令人满意的。     

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.     

Risk Assessment for Tuzla Naval Base Breakwater

1 .IntroductionTurkeywasstruckbytwomajorearthquakeeventsonAugust 1 7thandNovember 1 2th ,1 999,namedIzmit (Mw=7.4 )andD櫣ｚｃｅ (Mw=7.2 )earthquakes,respectively .ThestationsoperatedbytheGeneralDirectorateofDisasterAffairs,theKandilliObservatoryandEarthquakeResearchInstituteofIstanbulTechnicalUniversitymeasuredatleast 2 7stronggroundmotionsfortheIzmitearthquakewithin 2 0 0kmofthefault.AsignificantsegmentofthefaultrupturedintheareabetweenthewestofGolcukandtheeastofLakeSapanca .Inthesou…     

Model simulations of carbon sequestration in the northwest Pacific by patch fertilization

Iron fertilization of nutrient-rich surface waters of the ocean is one possible way to help slow the rising levels of atmospheric CO₂ by sequestering it in the oceans via biological carbon export. Here, I use an ocean general circulation model to simulate a patch of nutrient depletion in the subpolar northwest Pacific under various scenarios. Model results confirm that surface fertilization is an inefficient way to sequester carbon from the atmosphere (Gnanadesikan et al., 2003), since only about 20% of the exported carbon comes initially from the atmosphere. Fertilization reduces future production and thus CO₂ uptake by utilizing nutrients that would otherwise be available later. Effectively, this can be considered as leakage when compared to a control run. This “effective” leakage and the actual leakage of sequestered CO₂ cause a significant, rapid decrease in carbon retention (only 30–45% retained after 10 years and less than 20% after 50 years). This contrasts markedly with the almost 100% retention efficiency for the same duration using the same model, when carbon is disposed directly into the northwest Pacific (Matsumoto and Mignone, 2005). As a consequence, the economic effectiveness of patch fertilization is poor in two limiting cases of the future price path of carbon. Sequestered carbon in patch fertilization is lost to the atmosphere at increasingly remote places as time passes, which would make monitoring exceedingly difficult. If all organic carbon from one-time fertilization reached the ocean bottom and remineralized there, acidification would be about −0.05 pH unit with O₂ depletion about −20 μmol kg⁻¹. These anomalies are probably too small to seriously threaten deep sea biota, but they are underestimated in the model because of its large grid size. The results from this study offer little to advocate purposeful surface fertilization as a serious means to address the anthropogenic carbon problem.     

Modified PIC Method for Sea Ice Dynamics

The sea ice cover displays various dynamical characteristics such as breakup, rafting, and ridging under external forces. To model the ice dynamic process accurately, the effective numerical modeling method should be established. In this paper, a modified particle-in-cell （PIC） method for sea ice dynamics is developed coupling the finite difference （FD） method and smoothed particle hydrodynamics （SPH）. In this method, the ice cover is first discretized into a series of lagrangian ice particles which have their own sizes, thicknesses, concentrations and velocities. The ice thickness and concentration at Eulerian grid positions are obtained by interpolation with the Gaussian function from their surrounding ice particles. The momentum of ice cover is solved with FD approach to obtain the Eulerian cell velocity, which is used to estimate the ice particle velocity with the Gaussian function also. The thickness and concentration of ice particles are adjnsted with particle mass density and smooth length, which are adjusted with the redistribution of ice particles. With the above modified PIC method, numerical simulations for ice motion in an idealized rectangular basin and the ice dynamics in the Bohai Sea are carried out. These simulations show that this modified PIC method is applicable to sea ice dynamics simulation.     

波峰特征量的统计分布

本文在作者另文数值模拟得到的三维海浪基础上 ,进一步分析给出了三维波峰长度、高度及方向角度等特征量的统计分布。发现考虑海浪的方向性质后波峰分布不再是均匀连续的 ;波峰的长度分布受方向函数影响 ,方向分布越宽 ,波峰的平均长度越短 ,波峰的方向角分布越宽 ;波峰高度和长度在波峰高度较小时有很大的相关性 ,而在波峰高度很大时无关     

海南东寨港红树林港湾潮汐动力研究

在海南东寨港林市村附近红树林港湾潮汐动力测量的基础上，建立二维可移动边界有限元数值模型，研究红树林港湾浅水湖滩潮流场特征。实测与计算结果表明，潮沟区流速最大，具驻波特征，呈明显时间－速度不对称，东潮沟涨潮优势，百潮沟落潮优势，整个港湾水体呈顺时针循环。白滩与林地周期性淹没，白滩流速较湖沟略小。林地曼宁系数为白滩的１０倍，流速小于１０ｃｍ·ｓ（－１），约为相应潮沟流速的１／１０。林地如去掉林木，流速可增加３倍以上，流向也有变化。