大凌河口水动力数值模拟分析

采用大凌河实测长系列河道径流及水下地形等数据, 通过建立水动力数值模型, 综合考虑潮流和河道径流的影响, 对大凌河口的潮流进行模拟, 探讨大凌河口及辽东湾北部海域水动力过程的主要特征. 结果表明: 大凌河口附近海域的海流以潮流为主, 具有明显的往复性质. 潮流总的运动趋势是, 涨潮主流向为北东, 落潮主流向为南西. 大潮流速大于小潮流速, 涨潮流历时与落潮流历时几乎相当. 最大涨潮流速约为0.52 m/s, 最大落潮流速约为0.4 m/s, 潮流涨落平均潮流强度的分布大体和该海域等深线相适应.     

珠江网河水沙分配变化及其对伶仃洋水沙场的影响

采用二维ECOMSED模型对径流来水来沙变化前后的1978年和1999年伶仃洋洪水期水动力场和悬沙场进行模拟,结果表明:伶仃洋水动力场整体增强,但不同地貌单元水动力变化具有明显差异.分流比变化后洪水期伶仃洋西槽、东槽和中滩水动力增强,涨潮和落潮流速普遍增大0.1m/s;西部各口门区和西滩涨潮流势减弱,径流优势更加明显,但西滩三个口门水道及滩面潮沟落潮流速增大约0.1 m/s;东滩水动力特征变化不大.悬沙浓度亦整体升高,平均增加了0.02~0.05 kg/m3,悬沙高值区范围扩大,0.10 kg/m3等值线由西槽西侧明显向东推进到沿东槽南北向分布,径流来沙的影响范围东扩.     

磨刀门水道咸潮上溯动力特性分析

为研究磨刀门水道咸潮上溯的动力特性,基于非结构网格海洋模型（Finite Volume Coastal Ocean Model,FVCOM）,构建了覆盖珠江河口及其上游网河区的高分辨率三维斜压数值模型,采用实测资料对其进行率定和验证,并开展了咸潮上溯的数值模拟计算。根据计算结果和实测资料,对磨刀门水道大、中、小潮期间的盐淡水分层与混合特征、盐分物质的分层输移机制进行分析,探讨其咸潮上溯强度时空分布差异的原因。结果表明:小潮期,底层累积盐通量明显大于表层,净输移方向为陆向;大潮期,表层累积盐通量明显大于底层,净输移方向为海向;而平衡点一般出现于中潮期,这就是磨刀门水道咸潮上溯最强和最弱时刻分别出现于小潮和大潮后的中潮期的原因所在。     

长江河口涨落潮槽沉积物特征及其动力响应

涨落潮槽是河口区的重要地貌单元,涨落潮槽的水动力有着明显的差异。通过对涨潮槽新桥水道和南小泓以及落潮槽南支主槽和南港表层沉积物的粒度、粘土矿物、重矿物以及磁性特征分析,发现落潮槽表层沉积物的粒径较粗,为粉砂质砂,涨潮槽沉积物主要是砂质粉砂。在双向水流的作用下,粘土矿物重新发生分配,涨落潮槽粘土矿物的组分变化不大。涨潮槽的重矿物颗粒百分含量中,稳定的不透明矿物比落潮槽有所减少,而比重小的片状矿物有所增加,碳酸盐含量较高。磁性矿物的含量在不同的地方相差很大,落潮槽中的亚铁磁性矿物含量高于涨潮槽。这些沉积特性的不同是对涨落潮槽内水动力差异响应的结果。     

磨刀门河口环流与咸淡水混合层化机制

为研究磨刀门盐水混合层化特征,基于SCHISM模型,建立了三维盐度数值模型,根据实测资料对其进行验证。结合水体势能异常理论,对枯季磨刀门河口混合层化的时空变化特征及深槽与浅滩的层化机制差异进行分析。结果表明:磨刀门河口小潮时水体层化最强,中潮时水体层化最弱,且拦门沙至挂定角段水体层化始终较强。磨刀门深槽水体层化主要受纵向平流、纵向水深平均应变和垂向混合影响,而浅滩水体层化则受横向平流、横向水深平均应变和垂向混合影响;磨刀门河口表、底层水体湍动能耗散率较高,而中间水层存在低耗散区,且涨潮时湍动能耗散率比落潮时大。     

整治工程影响下的河口拦门沙航道回淤特征及成因——以射阳港航道为例

通过不同时期实测地形资料,对比分析射阳港拦门沙航道一期整治工程建设后航道的回淤分布特征;利用实测水文泥沙资料和潮流数学模型计算分析射阳河口拦门沙航道开挖后的水流、含沙量特征;探讨了射阳港拦门沙航道淤积的泥沙来源和淤积成因。研究表明:航道开挖后沿程普遍淤积,尤其是口门及堤根区段航道淤积严重。导致航道淤积的主要原因为:有掩护航道内涨潮流速大、落潮流速小,口门外高含沙量水体进入航道后泥沙很难被带出航道,导致航道内普遍淤积;口门内2 km范围内回淤较大是由于越堤流输沙和口门内侧涨潮期间存在明显回流;导堤根部航道段回淤较大是由于射阳河与黄沙港河汇流后河道放宽,落潮期间出河口的水体挟沙力突然降低。     

长江河口崇明东滩周期性淹水区域水流的基本特征

潮间带周期性淹水区域水深、流速的变化过程是潮滩水动力过程的基本组成部分,也是潮流与泥沙相互作用的基础。通过2002年4月至2003年1月4个季节的野外实测,获得了平静天气条件下,崇明东滩滩面潮流水深、流速与流向的变化过程数据。结果表明,崇明东滩盐沼和邻近光滩处涨潮历时均小于落潮历时,水深过程变化呈现出“陡涨缓落”的特点。光滩与盐沼交界处光滩一侧流速过程呈“双峰型”特征,涨落潮均出现流速峰值;盐沼（植物生长期）流速过程具有“单峰型”特点,仅在涨潮初出现峰值。研究区潮流不对称性明显,主要表现为涨潮优势,且由光滩向盐沼上部不断增强,潮沼植物和地形变化是加强盐沼区涨潮优势的主要原因。流速变化过程的差异和潮流不对称性使盐沼区域发生稳定的泥沙淤积,盐沼前缘光滩则会出现较频繁的冲淤变化,平静天气条件下,它们是控制崇明东滩泥沙输移和潮滩动力地貌过程的动力基础。     

磨刀门水道咸潮入侵及其变异分析

随着区域经济及人口的进一步增长,珠江三角洲磨刀门水道已成为江门、中山、珠海、澳门等城市重要的水源地,相应咸潮入侵对供水安全造成的危害也越来越严重。本文依据实测资料和现场调研,对磨刀门水道咸潮入侵的潮相变化、径流—咸潮的响应模式、风对咸潮入侵的影响及咸潮入侵的历史演变趋势及其变异、最小压咸流量及最佳压咸时机等进行了研究。研究结果表明,影响咸潮入侵的主要因素是径流,有利于咸潮入侵的风向为北风～东北风。以思贤滘流量(即马口水文站与三水水文站流量之和)为参照,为保障珠海、澳门供水安全,磨刀门水道最小压咸流量范围为2200-2700m3/s,平均流量为2450m3/s,最佳压咸补淡的时机为大潮转小潮期,也就是朔望(半月)周期的落潮期。     

长江口动量系数分布特征研究

基于长江口1996年3月和9月共23个测点实测水流资料,计算了动量系数的时空变化,并应用线性回归法,计算了特征时间段内的平均动量系数,获得了长江口动量系数时空分布的基本特征:时间上,涨潮大于落潮,洪季大于枯季,而大、中、小潮则依次递增;空间上,自口外向口内,北支呈逐渐增大,而南支呈先增大后减小、口门附近最大的基本特征.由动量系数分析了潮流的摩阻特性,为数值模拟合理选取曼宁系数提供参考.     

福建省东北部沿海罗源湾互花米草盐沼环境下粘性沉积物的侵蚀-沉降过程

根据在福建罗源湾互花米草盐沼内的沉积动力要素(水深、流速、悬沙浓度和沉降速率)的观测及表层沉积物的粒度分析结果,罗源湾互花米草盐沼表层沉积物组份以粉砂和粘土为主,为粘性沉积物;互花米草盐沼内的潮流流速很小,最大值不超过7cm/s,相应地,底部切应力也相对较小,悬沙在潮周期内多处于沉降状态,表层沉积物在涨潮初期和落潮后期多发生侵蚀.在台风显著影响下,潮流流速平均值变化很小,但水流紊动能量和底部切应力显著增大,导致表层沉积物多处于侵蚀状态,在台风登陆期间潮周期内净侵蚀通量可达到40.099kg/m2;台风影响过后,底部切应力逐渐减小,沉积物逐渐不发生侵蚀,悬沙多发生沉降,整个潮周期内以悬沙沉降为主.     

A study of the dynamics of a marine sandwave

The movement of the crest of a sandwave was studied using cross-sectional profiles obtained from lines of sea-bed reference stakes. Measurements were made, over a six month period, before and after flood and ebb tides in relation to both spring and neap tides and surface wave conditions. Additional observations were obtained on a daily basis, over an equinoctial neap to spring to neap tidal period, in conjunction with boundary layer flow measurements. Tracer experiments were conducted to study the dispersion of sediment from the sandwave crest. The results showed that the sandwave was relatively stable at neap tides, whilst at higher tidal ranges, the crest position oscillated with successive flood and ebb tides. Net flank erosion occurred on the less steep, upstream slope during the dominant ebb tide. This, together with increased deposition on the lee slope, caused the crest to advance. It was not possible to extrapolate sandwave migration over long periods as the tidal dynamic trends were interrupted by wind stress and surface wave activity. High particle orbital velocities, generated at the sea-bed by storm waves, caused major reductions in crestal heights. Calculated volumes of sediment eroded and accreted were used, with boundary layer flow measurements, to calculate threshold velocities for the movement of the sediment and sediment transport rates.     

Circulation and mixing due to tidal forcing in the Bertioga Channel, São Paulo, Brazil

Bertioga Channel is a partially mixed (type 2) tidal estuary on the coastal plain of São Paulo, Brazil. Hourly current and salinity measurements during neap and spring tides in July 1991 yielded information about the physical structure of the system. Peak along-channel velocities varied from 40 cm s^?1 to 60 cm s^?1 during flood tides and from 70 cm s^?1 to 100 cm s^?1 during ebb tides. Net vertical velocity profiles indicate that the net current reverses directions at a depth of 2.5–3.0 m in the halocline. Due to appreciable fortnightly tidal modulation, the estuary alternates from being highly stratified (type 2b) during neap tides, with advection and diffusion contributing equally to the net upstream salt flux, to being moderately stratified (type 2a) during spring tides, when 90% of the net upstream salt transport is the result of effective tidal diffusion. Decomposition of the salt flux indicates that the relative contribution to the upstream salt transport by gravitational circulation shear is greater than the oscillatory tidal flux by a factor of 2.6 during neap tides. The oscillatory tidal flux is generated by the correlation of the tidal components of the u-velocity and salinity and is responsible for approximately the same amount of upstream salt transport, during neap and spring tides. However, during spring tides, this oscillatory term is greater than the other salt flux terms by a factor of 1.4. The total salt transport, through a unit width of the section perpendicular to the flow, was within 2% of the sum of the seven major decomposed, advective and dispersive terms. On the assumption that the Bertioga Channel is laterally homogeneous, the results also indicate that the estuary is not in steady state with respect to salt flux.     

Tidally-induced shear stress variability above intertidal mudflats in the macrotidal seine estuary

Tidal currents and the spatial variability of tidally-induced shear stress were studied during a tidal cycle on four intertidal mudflats from the fluvial to the marine part of the Seine estuary. Measurements were carried out during low water discharge (<400 m³ s⁻¹) in neap and spring tide conditions. Turbulent kinetic energy, covariance, and logarithmic profile methods were used and compared for the determination of shear stress. The c_TKE coefficient value of 0.19 cited in the literature was confirmed. Shear stress values were shown to decrease above mudflats from the mouth to the fluvial part of the estuary due to dissipation of the tidal energy, from 1 to 0.2 N m⁻² for spring tides and 0.8 to 0.05 N m⁻² for neap tides. Flood currents dominate tidally-induced shear stress in the marine and lower fluvial estuary during neap and spring tides and in the upper fluvial part during spring tides. Ebb currents control tidally-induced shear stress in the upper fluvial part of the estuary during neap tides. These results revealed a linear relationship between friction velocities and current velocities. Bed roughness length values were calculated from the empirical relationship given by Mitchener and Torfs (1996) for each site; these values are in agreement with the modes of the sediment particle-size distribution. The influence of tidal currents on the mudflat dynamics of the Seine estuary was examined by comparing the tidally-induced bed shear stress and the critical erosion shear stress estimated from bed sediment properties. Bed sediment resuspension induced by tidal currents was shown to occur only in the lower part of the estuary.     

Effect of a broad shallow sill on tidal circulation and salt transport in the entrance to a coastal plain estuary (Mira—Vila Nova de Milfontes, Portugal)

We describe the tidal circulation and salinity regime of a coastal plain estuary that connects to the ocean through a flood tide delta. The delta acts as a sill, and we examine the mechanisms through which the sill affects exchange of estuarine water with the ocean. Given enough buoyancy, the dynamics of tidal intrusion fronts across the sill and selective withdrawal (aspiration) in the deeper channel landward appear to control the exchange of seawater with estuarine water. Comparison of currents on the sill and stratification in the channel reveals aspiration depths smaller than channel depth during neap tide. During neap tide and strong vertical stratification, seawater plunges beneath the less dense estuarine water somewhere on the sill. Turbulence in the intruding bottom layer on the sill promotes entrainment of fluid from the surface layer, and the seawater along the sill bottom is diluted with estuarine water. During ebb flow, salt is effectively trapped landward of the sill in a stagnant zone between the aspiration depth and the bottom where it can be advected farther upstream by flood currents. During spring tide, the plunge point moves landward and off the sill, stratification is weakened in the deep channel, and aspiration during ebb extends to the bottom. This prevents the formation of stagnant water near the bottom, and the estuary is flooded with high salinity water far inland. The neapspring cycle of tidal intrusion fronts on flood coupled with aspiration during ebb interacts with the sill to play an important role in the transport and retention of salt within the estuary.     

On the linkages among density, flow, and bathymetry gradients at the entrance to the Chesapeake Bay

Linkages among density, flow, and bathymetry gradients were explored at the entrance to the Chesapeake Bay with underway measurements of density and flow profiles. Four tidal cycles were sampled along a transect that crossed the bay entrance during cruises in April–May of 1997 and in July of 1997. The April–May cruise coincided with neap tides, while the July cruise occurred during spring tides. The bathymetry of the bay entrance transect featured a broad Chesapeake Channel, 8 km wide and 17 m deep, and a narrow North Channel, 2 km wide and 14 m deep. The two channels were separated by an area with typical depths of 7 m. Linkages among flows, bathymetry, and water density were best established over the North Channel during both cruises. Over this channel, greatest convergence rates alternated from the left (looking into the estuary) slope of the channel during ebb to the right slope during flood as a result of the coupling between bathymetry and tidal flow through bottom friction. These convergences were linked to the strongest transverse shears in the along-estuary tidal flow and to the appearance of salinity fronts, most markedly during ebb periods. In the wide channel, the Chesapeake Channel, frontogenesis mechanisms over the northern slope of the channel were similar to those in the North Channel only in July, when buoyancy was relatively weak and tidal forcing was relatively strong. In April–May, when buoyancy was relatively large and tidal forcing was relatively weak, the recurrence of fronts over the same northern slope of the Chesapeake Channel was independent of the tidal phase. The distinct frontogenesis in the Chesapeake Channel during the increased buoyancy period was attributed to a strong pycnocline that insulated the surface tidal flow from the effects of bottom friction, which tends to decrease the strength of the tidal flow over relatively shallow areas.     

A particle conveyor belt process in the Columbia River estuary: Evidence from chlorophyll<Emphasis Type="Italic">a</Emphasis> and particulate organic carbon

Using both the photosynthetically active chlorophylla (chla) content of the organic carbon fraction of suspended particulate matter (chla/POC) and the percentage of photosynthetically, active chla in fluorometrically measured chla plus pheophytina (% chla), we determined that under specified hydrodynamic conditions, neap-spring tidal differentiation in particle dynamics could be observed in the Columbia River estuary. During summer time neap tides, when river discharge was moderate, bottom chla/POC remained relatively unchanged from riverine chla/POC over the full 0–30 psu salinity range, suggesting a benign trapping environment. During summertime spring tides, bottom chla/POC decreased at mid range salinities indicating resuspension of chla-poor POC during flood-ebb transitions. Bottom % chla during neap tides tended to average higher than that during spring tides, suggesting that neap particles were more recently hydrodynamically trapped than those on the spring tides. Such differentiation supported the possibility of operation of a particle conveyor belt process, a process in which low-amplitude neap tides favor selective particle trapping in estuarine turbidity maxima (ETM)., while high-amplitude spring tides favor particle resuspension from the ETM. Untrapped river-derived particles at the surface would continue through the estuary to the coastal ocean on the neap tide; during spring tide some particles eroded from the ETM would combine with unsettled riverine particles in transit toward the ocean. Because in tensified biogeochemical activity is associated with ETM, these neap-spring differences may be critical to maintenance and renewal of populations and processes in the estuary. Very high river discharge (15, 000 m³ s⁻¹) tended to overwhelm neap-spring differences, and significant oceanic input during very low river discharge (5,000 m³ s⁻¹) tended to do the same in the estuarine channel most exposed to ocean input. During heavy springtime phytoplankton blooms, development of a thick bottom fluff layer rich in chla also appeared to negate neapspring differentiation because spring tides apparently acted to resuspend the same rich bottom material that was laid down during neap tides. When photosynthetic assimilation numbers [μgC (μgchl,a)⁻¹h⁻¹] were measured across, the full salinity range, no neap-spring differences and no river discharge effects occurred, indicating that within our suite of measurements the compositional distinction of suspended particulate material was mainly a function of chla/POC, and to a lesser extent % chla. Even though these measurements suggest the existence of a conveyor belt process, proof of actual operation of this phenomenon requires scalar flux measurements of chla properties in and out of the ETM on both neap and spring tides.     

磨刀门河口洪季波生流及其泄洪影响

磨刀门已由"径流型"向"径流-波浪型"河口转变,波浪已是该河口主要动力之一,但波浪对河口洪季水流及泄洪的影响缺少研究。在2-D潮流数学模型中添加随潮位实时变化的波浪辐射应力,建立波浪潮流耦合数学模型;波浪求解采用缓坡方程,背景水深由潮流模型实时提供,可通过比较考虑和未考虑波浪影响的河口流场来分析波浪对泄洪的影响。在年均常浪作用下,磨刀门河口洪季涨落潮阶段均有明显的波生环流结构。由于波浪作用方向向陆,波生流减弱了浅滩区的向海余流,增大了浅滩向陆余流;受浅滩向海余流减弱影响,河口动力自调整后形成归槽水流,促使深槽内向海余流增大。波浪有顶托河口泄洪之势,可改变滩槽泄洪分配比例;年均常浪的波高较小,其对潮流及泄洪的影响区域限制在浅水区,故对泄洪的负面影响有限。