The morphodynamic modelling of tidal sand waves on the shoreface

An artificial sand wave on the Dutch shoreface of the North Sea has been studied in conditions with relatively strong tidal currents in the range of 0.5 to 1 m/s and sediments in the medium sand size range of 0.2 to 0.5 mm. The sand wave is perpendicular to the tidal current and has a maximum height and length of the order of 5 m and 1 km, respectively. The sand wave is dynamically active and shows migration rates of the order of a few metres per year. A numerical morphodynamic model (DELFT3D model) has been used to simulate the morphological behaviour of the sand wave in the North Sea. This model approach is based on the numerical solution of the three-dimensional shallow water equations in combination with a surface wave propagation model (wind waves) and the advection–diffusion equation for the sediment particles with online bed updating after each time step. The model results show that the sand wave grows in the case of dominant bed-load transport (weak tidal currents; relatively coarse sediment; small roughness height; low waves) and that the sand wave decays in the case of dominant suspended transport (strong currents, relatively fine sediment, large roughness height; storm waves).     

Data assimilation for a coastal area morphodynamic model: Morecambe Bay

This paper investigates the use of data assimilation in coastal area morphodynamic modelling using Morecambe Bay as a study site. A simple model of the bay has been enhanced with a data assimilation scheme to better predict large-scale changes in bathymetry observed in the bay over a 3-year period. The 2DH decoupled morphodynamic model developed for the work is described, as is the optimal interpolation scheme used to assimilate waterline observations into the model run. Each waterline was acquired from a SAR satellite image and is essentially a contour of the bathymetry at some level within the inter-tidal zone of the bay. For model parameters calibrated against validation observations, model performance is good, even without data assimilation. However the use of data assimilation successfully compensates for a particular failing of the model, and helps to keep the model bathymetry on track. It also improves the ability of the model to predict future bathymetry. Although the benefits of data assimilation are demonstrated using waterline observations, any observations of morphology could potentially be used. These results suggest that data assimilation should be considered for use in future coastal area morphodynamic models.     

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.     

Wind-induced drift of objects at sea: The leeway field method

A method for conducting leeway field experiments to establish the drift properties of small objects (0.1-25 m) is described. The objective is to define a standardized and unambiguous procedure for condensing the drift properties down to a set of coefficients that may be incorporated into existing stochastic trajectory forecast models for drifting objects of concern to search and rescue operations and other activities involving vessels lost at sea such as containers with hazardous material.An operational definition of the slip or wind and wave-induced motion of a drifting object relative to the ambient current is proposed. This definition taken together with a strict adherence to a 10 m wind speed allows us to refer unambiguously to the leeway of a drifting object. We recommend that all objects if possible be studied using what we term the direct method, where the object’s leeway is studied directly using an attached current meter.We establish a minimum set of parameters that should be estimated for a drifting object for it to be included in the operational forecast models used for prediction of search areas for drifting objects.We divide drifting objects into four categories, depending on their size. For the smaller objects (less than 0.5 m), an indirect method of measuring the object’s motion relative to the ambient current must be used. For larger objects, direct measurement of the motion through the near-surface water masses is strongly recommended. Larger objects are categorized according to the ability to attach current meters and wind monitoring systems to them.The leeway field method proposed here is illustrated with results from field work where three objects were studied in their distress configuration; a 1:3.3 sized model of a 40-foot Shipping container, a World War II mine and a 220 l (55-gallon) oil drum.     

The continent-ocean transition at the southeastern margin of the South China Sea

The South China Sea formed by magma-poor, or intermediate volcanic rifting in the Paleogene. We investigate the structure of the continent-ocean transition (COT) at its southern margin, off NW Palawan between the continental blocks of Reed Bank and the islands of Palawan and Calamian. Several surveys, recorded by the BGR from 1979 to 2008, established a comprehensive database of regional seismic lines, accompanied with magnetic and gravity profiles.We interpret two major rifted basins, extending in the NE direction across the shelf and slope, separated by a structural high of non volcanic origin.The continent-ocean transition is interpreted at the seaward limit of the continental crust, when magnetic spreading anomalies terminate some 80-100 km farther north. The area in between displays extensive volcanism - as manifest by extrusions that occasionally reach and cut the seafloor, by dykes, and by presumed basaltic lava flows - occurring after break-up.The COT is highly variable along the NW Palawan slope: One type shows a distinct outer ridge at the COT with a steep modern seafloor relief. The other type is characterised by rotated fault blocks, bounded by listric normal faults ramping down to a common detachment surface. Half-grabens developed above a strongly eroded pre-rift basement. The seafloor relief is smooth across this other type of COT.We suggest the pre-rift lithospheric configuration had major influence on the formation of the COT, besides transfer zones. Volcanic domains, confined to the north of competent crustal blocks correlate with the style of the COT.Gravity modelling revealed an extremely thinned crust across the shelf. We propose a depth-dependent extension model with crust being decoupled from mantle lithosphere, explaining the discrepancy of subsidence observed across the South China Sea region.     

贡湖及其相关水域风生流模拟研究

在太湖风生流的数值模拟研究的基础上，设计模拟贡湖及其相关水域风生流的套网格数值模式。应用该模式研究不同风情下，梅梁湖，太湖主体湖水对贡湖的影响。     

Three-dimensional Layer-integrated Modelling of Estuarine Flows with Flooding and Drying

Details are given of the refinement and application of a thee-dimensional (3-D) layer-integrated numerical model of tidal circulation, with the aim of simulating severe tidal conditions for practical engineering applications. The mode splitting strategy has been used in the model. A set of depth-integrated 2-D equations are first solved to give the pressure gradient, and the layer-integrated 3-D equations are then solved to obtain the vertical distributions of the flow velocities. Attention has been given to maintaining consistency of the physical quantities derived from the 2-D and 3-D equations. A TWO=layer mixing length turbulence model for the vertical shear stress distribution has been included in the model. Emphasis has been focused on applying the model to a real estuary, which is geometrically complicated and has large tidal ranges giving rise to extensive flooding and drying. The model has been applied to three examples, including: wind-driven flow in a rectangular lake, tidal circulation in a model rectangular harbour, and tidal circulation in a large estuary. Favourable results have been obtained for both the simple and complex flow beds.     

海洋调查数据的可视化建模及实现方法

科学计算可视化作为新一代的分析工具 ,正在进入物理海洋学领域。文中介绍对海洋调查数据进行可视化的关键问题——数据建模 ,并从分析海洋调查数据的特点出发提出针对海洋研究中可视化建模的原则和有效途径 ,给出相关的实现方法。其中 ,重点讨论局部体样条建模中必须面对的区域分割优化问题 ,引入遗传算法并给出较好的解决方案。