The dimensions of sand ripples in full-scale oscillatory flows

New large-scale experiments have been carried out in two oscillatory flow tunnels to study ripple regime sand suspension and net sand transport processes in full-scale oscillatory flows. The paper focuses on ripple dimensions and the new data are combined with existing data to make a large dataset of ripple heights and lengths for flows with field-scale amplitudes and periods. A feature of the new experiments is a focus on the effect of flow irregularity. The combined dataset is analysed to examine the range of hydraulic conditions under which oscillatory flow ripples occur, to examine the effects of flow irregularity and ripple three-dimensionality on ripple dimensions and to test and improve existing methods for predicting ripple dimensions.The following are the main conclusions. (1) The highest velocities in a flow time-series play an important role in determining the type of bedform occurring in oscillatory flow. Bedform regime is well characterised by mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth peak velocities in the case of irregular flow. (2) For field-scale flows, sand size is the primary factor determining whether equilibrium ripples will be 2D or 3D. 2D ripples occur when the sand D₅₀ ≥ 0.30 mm and 3D ripples occur when D₅₀ ≤ 0.22 mm (except when the flow orbital diameter is low). (3) Ripple type (2D or 3D) is the same for regular and irregular flows and ripple dimensions produced by equivalent regular and irregular flows follow a similar functional dependence on mobility number, with mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth velocities in the case of irregular flow. For much of the ripple regime, ripple dimensions have weak dependency on mobility number and ripple dimensions are similar for regular and irregular flows with the same flow orbital amplitude. However, differences in ripples produced by equivalent regular and irregular flows become significant at the high mobility end of the ripple regime. (4) Ripple dimensions predicted using the Wiberg and Harris formulae are in poor agreement with measured ripple dimensions from the large-scale experiments. Predictions based on the Mogridge et al. and the Nielsen formulae show better overall agreement with the data but also show systematic differences in cases of 3D ripples and ripples generated by irregular flows. (5) Based on the combined large-scale data, modifications to the Nielsen ripple dimension equations are proposed for the heights and lengths of 2D ripples. The same equations apply to regular and irregular flows, but with mobility number appropriately defined. 3D ripples are generally smaller than 2D ripples and estimates of 3D ripple height and length may be obtained by applying multipliers of 0.55 and 0.73 respectively to the 2D formulae. The proposed modified Nielsen formulae provide an improved fit to the large-scale data, accounting for flow irregularity and ripple three-dimensionality.     

Second-order moment advection scheme applied to Arctic Ocean simulation

We apply the second-order moment (SOM) advection scheme of (Prather, M.J. 1986. Numerical advection by conservation of second-order moments. J. Geophys. Res. 91, 6671–6681.) to the simulation of the large-scale circulation of the Arctic Ocean with a coupled ocean–sea-ice model. Compared to three other advection schemes commonly employed in ocean simulations (centred differences, flux corrected transport, and multidimensional positive definite advection transport), the SOM method helps preserve the vertical structure of Arctic water masses. The depth, thickness and hydrographic properties of the Arctic Surface Water and the Arctic Atlantic Layer are better represented with SOM than with any of the other three advection algorithms. We also present a convenient method for calculating the implicit numerical diffusivity of upstream based schemes, such as the SOM method, and discuss three approaches for improving the monotonicity properties of the SOM algorithm.     

Electronic atlas of the Russian Arctic coastal zone

A set of digital maps including geology, Quaternary sediments, landscapes, engineering-geological, vegetation, geocryological and the series of regional sources have been selected to characterize the Russian Arctic coast. Based on this data, new maps of engineering geocryological zoning and zoning of the coast with respect to the intensity of exogenous geological processes and risk of technogenic impacts have been generated at the scales of 1:4,000,000–1:8,000,000. These maps are a tool to assess the impact of industry on the Arctic coast of the country.     

Modified mild slope equation and open boundary conditions

A numerical model to compute wave field is developed. It is based on the Berkhoff diffraction-refraction equation, in which an energy dissipation term is added, to take into account the breaking and the bottom friction phenomena. The energy dissipation function, by breaking and by bottom friction, is introduced in the Berkhoff equation to obtain a new equation of propagation.The resolution is done with the hybrid finite element method, where lagrangians elements are used.     

Preferential concentration of marine particles in isotropic turbulence

The effect of small-scale turbulence on marine and aquatic particle transport has traditionally been to act as a means for creating homogeneous distributions. However, previous numerical simulations of heavy particle transport in turbulent flows have shown that particles are preferentially concentrated by turbulence and that effects of preferential concentration are most pronounced for particle parameters comparable to the Kolmogorov scales. Therefore, the focus of the present work is examination of the preferential concentration of marine particles. Application of Kolmogorov scaling indicates that effects of preferential concentration may be important for marine particles with diameters of order 1 mm in the upper mixed layer. Numerical simulations of unstratified isotropic turbulence are then used to support the notion that preferential concentration of particles possessing material characteristics representative of those encountered in marine environments can occur. In the simulations, particles of order 1 mm diameter are idealized as rigid spheres with a density ratio of 1.005. Simulation results demonstrate preferential concentration with peak particle number densities ranging from 10 to 60 times the global mean value. Implications of preferential concentration are also discussed, together with the limitations of the approach employed in the present study.     

Simulation of the Atlantic Ocean seasonal variability, using a quasi-isopycnic model

The paper discusses the data derived from a numerical experiment on the ocean’s response (between the equator and 64°N) to the seasonal variability of the atmospheric forcing (wind and heat flux through the ocean surface). A multilayer (7 layers) non-linear model is used incorporating the upper mixed layer interacting with the internal layers in the regimes of entraining and subduction. The restructuring of the layer composition, the currents and temperature variability, as well as the alternation of the entrainment and subduction regimes are analysed. Translated by Vladimir A. Puchkin.     

Effects of benthic biology on bottom boundary layer processes, Dry Tortugas Bank, Florida Keys

 An instrumented tetrapod was deployed for three weeks on the Dry Tortugas Bank at a depth of 26 m in February 1995. Bottom roughness was dominated by shrimp burrows and worm mounds with rms roughness amplitudes ranging from 0.47 to 1.75 cm. Logarithmic velocity profiles show apparent total roughness heights ranging from 0.30 to 1.45 cm, values consistent with observed biological roughness. The bed sediments were weakly bound by an algal crust at the sediment–water interface. When this bound layer was scraped away by a mooring that was accidentally dragged, sharp-crested wave-induced ripples appeared within the resulting swath. We conclude that physically induced roughness is biologically suppressed, but if dominant, would be significantly higher than the prevailing biological roughness.