日本东北M_W9.0地震海啸在港池及邻近区域诱发的涡流危险性计算与评估分析

海啸造成的灾害与损失并非都与淹没有关,特别是港口中海啸诱导的强流会对船只及海事设施产生重要的影响及损害.由于海啸流观测数据稀缺及海啸诱导涡流机制的不确定性,过去60年海啸科学主要集中于对海啸波特征及淹没过程的研究与分析,海啸流模拟及验证工作开展较少,导致对海啸流基本特征及其造成灾害现象的曲解.开展海啸诱导的涡流研究及预警服务显得尤为重要及紧迫.考虑快速海啸预警需要,综合对比海啸诱导涡流的物理框架及模型方法,探索兼顾效率与计算精度的海啸流模拟方法是本文的核心工作及出发点.通过分析浅层湍流相干结构(TCS)产生的主要物理耗散机制,确定了考虑2D水平耗散机制的非线性浅水方程可用于海啸涡流的模拟分析.基于高精、高分辨率有限体积模型Geoclaw建立了三个精细化的港口海啸流模型,模型分辨率为5m.利用基于海啸浮标反演的海啸源模型作为初始条件,模拟分析了日本东北地震海啸在远场的海啸波流特征.海啸波流特征模拟结果与观测吻合较好,结果可信.对比发现:波驱动的自由表面流,小的位相或波幅误差就会导致大的流速误差,流的模拟和预报相对波幅来说更具挑战性.研究了海啸波流能量在港池中的分布特征,得到:港池入口及防波堤两端常被强流控制,具有极高的危险性;相对于波幅的空间变化,海啸流具有更强的空间敏感性;所建立的高分辨率海啸模型模拟再现了日本海啸在近场的涡旋结构,给出了与观测基本一致的涡流特征.最后,引入海啸流危险等级标准,分析了港口海啸流危险性等级分布、船只疏散的安全深度及回港的时间周期.针对港口、海湾同时考虑海啸波流特征的海啸预警与评估对于港口应急管理者科学决策具有重要意义.     

Turbulent dissipation and mixing in Prydz Bay

In this paper,we present measurements of velocity,temperature,salinity,and turbulence collected in Prydz Bay,Antarctica,during February,2005.The dissipation rates of turbulent kinetic energy(ε) and diapycnal diffusivities(K z) were estimated along a section in front of the Amery Ice Shelf.The dissipation rates and diapycnal diffusivities were spatially non-uniform,with higher values found in the western half of the section where ε reached 10-7 W/kg and K z reached 10-2 m 2 /s,about two and three orders of magnitude higher than those in the open ocean,respectively.In the western half of the section both the dissipation rates and diffusivities showed a high-low-high vertical structure.This vertical structure may have been determined by internal waves in the upper layer,where the ice shelf draft acts as a possible energy source,and by bottom-generated internal waves in the lower layer,where both tides and geostrophic currents are possible energy sources.The intense diapycnal mixing revealed in our observations could contribute to the production of Antarctic Bottom Water in Prydz Bay.     

海山当地物源和南极底层水对富钴结壳成矿作用的影响——来自海山周围水柱化学分析的证据

为了了解海山富钴结壳中金属元素的来源及其形成环境,对中太平洋某海山周围两个站位(一个位于海山斜坡上部,另一个位于海山外围较远处)进行了海水水柱综合水化学分析,发现两个站位溶解氧、pH及营养盐的分布是一致的,反映了两个站位海水化学环境的相似性,而剖面上金属元素含量却存在很大差异,位于海山斜坡上部站位的金属元素含量明显高于离海山较远的站位,说明海山结壳的主要金属来源于海山当地的水岩反应甚至低温热液活动。尽管存在差异,但在两个不同站位均存在近底层金属元素含量高,这印证了通过势温分布推断的本区存在南极底层水(AABW)影响的判断,说明AABW在本研究区确实存在,并且是维持海山斜坡强氧化环境的主要因素,AABW可能同时起着输运部分成矿元素的作用,但两个站位水柱金属元素含量的差异主要揭示了海山区存在水岩反应或者低温热液等当地物质的输入。     

The rate of sedimentation from turbulent suspension: An experimental model with application to pyroclastic density currents and discussion on the grain‐size dependence of flow runout

Large‐scale experiments generating ground‐hugging multiphase flows were carried out with the aim of modelling the rate of sedimentation, of pyroclastic density currents. The current was initiated by the impact on the ground of a dense gas‐particle fountain issuing from a vertical conduit. On impact, a thick massive deposit was formed. The grain size of the massive deposit was almost identical to that of the mixture feeding the fountain, suggesting that similar layers formed at the impact of a natural volcanic fountain should be representative of the parent grain‐size distribution of the eruption. The flow evolved laterally into a turbulent suspension current that sedimented a thin, tractive layer. A good correlation was found between the ratio of transported/sedimented load and the normalized Rouse number of the turbulent current. A model of the sedimentation rate was developed, which shows a relationship between grain size and flow runout. A current fed with coarser particles has a higher sedimentation rate, a larger grain‐size selectivity and runs shorter than a current fed with finer particles. Application of the model to pyroclastic deposits of Vesuvius and Campi Flegrei of Southern Italy resulted in sedimentation rates falling inside the range of experiments and allowed definition of the duration of pyroclastic density currents which add important information on the hazard of such dangerous flows. The model could possibly be extended, in the future, to other geological density currents as, for example, turbidity currents.     

Contour current imprints and contourite drifts in the Bahamian archipelago

New data collected along the slopes of Little and Great Bahama Bank and the abyssal plain of the Bahama Escarpment provides new insights about contour current‐related erosive structures and associated deposits. The Bahamian slope shows abundant evidence of bottom current activity such as furrows, comet‐like structures, sediment waves and drifts. At a seismic scale, large erosion surfaces and main periods of drift growth resulted from current acceleration related to plate tectonic processes and progressive opening and closure of gateways and long‐term palaeoclimate evolution. At present‐day, erosion features and contourite drifts are either related to relatively shallow currents (<1000 m water depth) or to deep currents (>2500 m water depth). It appears that the carbonate nature of the drifts does not impact the drift morphology at the resolution addressed in the present study. Classical drift morphologies defined in siliciclastic environments are found, such as mounded, plastered and separated drifts. In core, contourite sequences show a bi‐gradational trend that resembles classical contourite sequences in siliciclastic deposits showing a direct relationship with a change in current velocity at the sea floor. However, in a carbonate system the peak in grain size is associated with increased winnowing rather than increased sediment supply as in siliciclastic environments. In addition, the carbonate contourite sequence is usually thinner than in siliciclastics because of lower sediment supply rates. Little Bahama Bank and Great Bahama Bank contourites contain open‐ocean input and slope‐derived debris from glacial episodes. Inner platform, platform edge and open ocean pelagic input characterize the classical periplatform ooze during interglacials. In all studied examples, the drift composition depends on the sea floor topography surrounding the drift location and the type of sediment supply. Carbonate particles are derived from either the slope or the platform in slope and toe of slope drifts, very deep contourites have distant siliciclastic sources of sediment supply. The recent discovery of the importance of a large downslope gravitary system along Bahamian slopes suggests frequent interactions between downslope and along‐slope (contour currents) processes. The interlayering of mass flow deposits and contourites at a seismic scale or the presence of surface structures associated with both contour currents and mass flow processes shows that both processes act at the same location. Finally, contour currents have an important impact on the repartition of deep‐water coral mounds. Currents can actively interact with mounds as a nutrient and oxygen supplier or have a passive interaction, with mounds solely being obstacles orienting erosion and deposition.     

Pyroclastic dune bedforms: macroscale structures and lateral variations. Examples from the 2006 pyroclastic currents at Tungurahua (Ecuador)

Pyroclastic currents are catastrophic flows of gas and particles triggered by explosive volcanic eruptions. For much of their dynamics, they behave as particulate density currents and share similarities with turbidity currents. Pyroclastic currents occasionally deposit dune bedforms with peculiar lamination patterns, from what is thought to represent the dilute low concentration and fluid‐turbulence supported end member of the pyroclastic currents. This article presents a high resolution dataset of sediment plates (lacquer peels) with several closely spaced lateral profiles representing sections through single pyroclastic bedforms from the August 2006 eruption of Tungurahua (Ecuador). Most of the sedimentary features contain backset bedding and preferential stoss‐face deposition. From the ripple scale (a few centimetres) to the largest dune bedform scale (several metres in length), similar patterns of erosive‐based backset beds are evidenced. Recurrent trains of sub‐vertical truncations on the stoss side of structures reshape and steepen the bedforms. In contrast, sporadic coarse‐grained lenses and lensoidal layers flatten bedforms by filling troughs. The coarsest (clasts up to 10 cm), least sorted and massive structures still exhibit lineation patterns that follow the general backset bedding trend. The stratal architecture exhibits strong lateral variations within tens of centimetres, with very local truncations both in flow‐perpendicular and flow‐parallel directions. This study infers that the sedimentary patterns of bedforms result from four formation mechanisms: (i) differential draping; (ii) slope‐influenced saltation; (iii) truncative bursts; and (iv) granular‐based events. Whereas most of the literature makes a straightforward link between backset bedding and Froude‐supercritical flows, this interpretation is reconsidered here. Indeed, features that would be diagnostic of subcritical dunes, antidunes and ‘chute and pools’ can be found on the same horizon and in a single bedform, only laterally separated by short distances (tens of centimetres). These data stress the influence of the pulsating and highly turbulent nature of the currents and the possible role of coherent flow structures such as Görtler vortices. Backset bedding is interpreted here as a consequence of a very high sedimentation environment of weak and waning currents that interact with the pre‐existing morphology. Quantification of near‐bed flow velocities is made via comparison with wind tunnel experiments. It is estimated that shear velocities of ca 0·30 m.s^?1 (equivalent to pure wind velocity of 6 to 8 m.s^?1 at 10 cm above the bed) could emplace the constructive bedsets, whereas the truncative phases would result from bursts with impacting wind velocities of at least 30 to 40 m.s^?1.     

Optimization of the Sefid-Roud Dam desiltation process using a sophisticated one-dimensional numerical model

Although water and soil conservation activities reduce reservoir sedimentation, it is inevitable that reservoirs fed by rivers transporting high amounts of sediment will experience sedimentation. The Ghezel-Ozan and Shah-Roud rivers, which flow to the Sefld-Roud reservoir dam, are both highly sediment-laden and transport significant amounts of sediment in both bed load and suspended load forms to the reservoir. Hence, it seems that the only practical way to remove the sediment from the reservoir is to flush it out using the Chasse method. In the present paper, field measurements of Chasse operation characteristics taken in previous years are presented, and a numerical model that simulates this process is introduced. After calibrating the model using field measured data, the calculated results （for reservoir pressure flushing and released sediment volume） of the numerical model were compared with other measured data for the same Chasse operation and the results agree well. Finally, using the numerical simulation results, the best approaches to ensure highly effective flushing while conserving reservoir water are presented （at least for the Sefid-Roud dam）. The operation of the bottom outlet gates, the shape of the output hydrograph, and the reservoir water level variation during flushing were optimized. In addition, the numerical model and related parameters, which need to be calibrated, are discussed.     

Numerical simulation of sedimentation processes in reservoirs as a function of outlet location

This paper builds on a recently published one-dimensional moving-boundary model of the coevolution of topset, foreset and bottomset in a reservoir that captures the dynamics of the internal muddy pond typical to reservoirs. This model was modified to account for different outlet locations at the reservoir＇s downstream end. This model considers a river carrying two sustained phases of sediments： coarse （sand） and fine （mud）. The coarse phase deposits in the topset and delta foreset, while the fine phase forms a dilute suspension of wash load in the river. As the river enters the reservoir, the muddy water plunges on the foreset to form a Froude-supercritical （purely depositional） turbidity current. This turbidity current emplaces the bottomset. The modified numerical model was tested against five laboratory experiments previously reported by the author. The model successfully locates the muddy-water/clear-water interface. In addition, modeled and measured bed deposits are in good agreement. Results clearly indicate that the location of the internal hydraulic jump plays a key role in the final bed deposit.     

Analytical study of the transverse distribution of along-channel and transverse residual flows in tidal estuaries

We present an analytical model to decompose complex along-channel and transverse residual flows into components induced by individual mechanisms. The model describes the transverse distribution of residual flows in tidally dominated estuaries. Scaling and perturbation techniques are used to obtain analytical solutions for residual flows over arbitrary across-channel bed profiles. The flows are induced by horizontal density gradients, tidal rectification processes, river discharge, wind, channel curvature and the earth's rotation. These rectification processes induce residual flows that are up-estuary to the right and down-estuary to the left of an estuarine channel (looking up-estuary in the northern hemisphere). The tidal rectification processes fundamentally change the transverse structure of along-channel residual flows in many tidal estuaries, as these processes cause the flows to be internally asymmetric about the mid-axis of the channel for relatively large tidal velocities, steep channels or narrow estuaries. In addition, velocity scales are derived from the analytical solutions to estimate the relative importance of the various residual flow mechanisms from estuarine parameters. A case study of a transect across the Upper Chesapeake Bay showed that important features of the residual flow observed in that transect are reproduced and explained by the analytical model. The velocity scales were able to identify the relevant residual flow mechanisms as well. The tidal rectification processes considered here result from advection of along-channel tidal momentum by Coriolis-induced transverse tidal currents.