风成沙波纹的数值模拟及稳定性分析

本文基于颗粒的质量守恒原理,考虑了沙粒的跃移、溅移和蠕移运动,建立了风成沙波纹的非线性积分微分演化方程,并进而通过尺度变换得到了近似的微分方程。数值求解该微分方程,得到了与野外观测相一致的沙波纹斑图,同时再现了沙波纹生成、合并的非线性行为,为进一步风成地貌形成演化的理论分析提供了一种可行的数学物理方法。     

Concentration profiles for fine and coarse sediments suspended by waves over ripples: An analytical study with the 1-DV gradient diffusion model

Field and laboratory measurements of suspended sediments over wave ripples show, for time-averaged concentration profiles in semi-log plots, a contrast between upward convex profiles for fine sand and upward concave profiles for coarse sand. Careful examination of experimental data for coarse sand shows a near-bed upward convex profile beneath the main upward concave profile. Available models fail to predict these two profiles for coarse sediments. The 1-DV gradient diffusion model predicts the main upward concave profile for coarse sediments thanks to a suitable β(y)$\beta (y)$-function (where β$\beta$ is the inverse of the turbulent Schmidt number and y   is the distance from the bed). In order to predict the near-bed upward convex profile, an additional parameter α$\alpha$ is needed. This parameter could be related to settling velocity (α$\alpha$ equal to inverse of dimensionless settling velocity) or to convective sediment entrainment process. The profiles are interpreted by a relation between second derivative of the logarithm of concentration and derivative of the product between sediment diffusivity and α$\alpha$.     

An empirical model of subcritical bedform migration

The ability to predict bedform migration in rivers is critical for estimating bed material load, yet there is no relation for predicting bedform migration (downstream translation) that covers the full range of conditions under which subcritical bedforms develop. Here, the relation between bedform migration rates and transport stage is explored using a field and several flume data sets. Transport stage is defined as the non‐dimensional Shields stress divided by its value at the threshold for sediment entrainment. Statistically significant positive correlations between both ripple and dune migration rates and transport stage are found. Stratification of the data by the flow depth to grain‐size ratio improved the amount of variability in migration rates that was explained by transport stage to ca 70%. As transport stage increases for a given depth to grain‐size ratio, migration rates increase. For a given transport stage, the migration rate increases as the flow depth to grain‐size ratio gets smaller. In coarser sediment, bedforms move faster than in finer sediment at the same transport stage. Normalization of dune migration rates by the settling velocity of bed sediment partially collapses the data. Given the large amount of variability that arises from combining data sets from different sources, using different equipment, the partial collapse is remarkable and warrants further testing in the laboratory and field.     

Surface and subsurface sedimentary structures produced by salt crusts

The growth and subsequent dissolution of salts on or within sediment may alter sedimentary structures and textures to such an extent that it is difficult to identify the depositional origin of that sediment and, as a result, the sediment may be misinterpreted. To help to overcome such problems with investigating ancient successions, results are presented from a comprehensive study of the morphology and fabrics of three large areas of modern salt flats in SE Arabia: the Sabkhat Matti inland region and the At Taf coastal region, both in the Emirate of Abu Dhabi, and the Umm as Samim region in Oman. These salt flats are affected by tidal‐marine, alluvial and aeolian depositional processes and include both clastic‐ and carbonate‐dominated surficial sediments. The efflorescent and precipitated salt crusts in these areas can be grouped into two main types: thick crusts, with high relief (>10 cm) and a polygonal or blocky morphology; or thin crusts, with low relief (<10 cm) and a polygonal or blister‐like appearance. The thin crusts may assume the surface morphology of underlying features, such as ripples or biogenic mats. A variety of small‐scale textures were observed: pustular growths, hair‐like spikes and irregular wrinkles. Evolution of these crusts over time results in a variety of distinctive sedimentary fabrics produced by salt‐growth sediment deformation, salt‐solution sediment collapse, sediment aggradation and compound mixtures of these processes. Salt‐crust processes produce features that may be confused with aeolian adhesion structures. An example from the Lower Triassic Ormskirk Sandstone Formation of the Irish Sea Basin demonstrates how this knowledge of modern environments improves the interpretation of the rock record. A distinctive wavy‐laminated facies in this formation had previously been interpreted as the product of fluvial sheetfloods modified by soft‐sediment deformation and bioturbation. Close inspection of laminations seen in core reveals many of the same sedimentary fabrics seen in SE Arabia associated with salt crusts. This facies is the product of salt growth on aeolian sediment and is not of fluvial origin.     

Evaluation of grain shear stresses required to initiate movement of particles in natural rivers

Shear stresses were evaluated at different sites on two rivers. The first (the Rulles) is characterized by a pebbly bedload and a meandering bed with riffles and pools. The second (the Rouge Eau) has mainly a sandy rippled bed where meandering is well developed but also flat gravelly sectors without meandering system. Shear stresses calculated from friction velocities (τ*) using a redefined y₁ roughness height parameter were compared with total shear stresses calculated from the energy grade line and the hydraulic radius (τ), Divergence between these shear stresses seems to increase in the presence of bedforms and large-scale irregularities of the channel. The τ*/τ ratio is close to 0·5 in the gravelly sector of the Rouge Eau and reaches 0·65 in the riffles of the Rulles (generally located at the inflexion point of the meanders), while it is less than 0·3 in the pools of the same river (located in the loops) and only 0·2 in the sandy rippled sector of the Rouge Eau. Grain and bedform shear stresses were evaluated at these same sites by different methods. The grain shear stress (τ') represents on average 30 per cent of the total shear stress in the riffles of the Rulles and the gravelly sector of the Rouge Eau, but less than 15 per cent in the pools in the Rulles and the sandy sectors of the Rouge Eau. However, it emerges from experiments conducted with marked pebbles and in situ observations of erosion and transport of sandy and gravelly particles, that the grain shear stresses are underestimated and cannot explain the movements and modifications actually observed. Conversely, shear stresses calculated from friction velocities at the sites where erosion actually occurred (or failed to occur despite very high velocities) provide a better explanation of the observed movements.     

Morphodynamics and cobbles behavior in and near the surf zone

The evolution of an initially flat sandy slope and the dynamics of large objects (cobbles/mines) emplaced on it are studied in a laboratory wave tank under simulated surf conditions. Upon initiation of wave forcing, the initially flat beach undergoes bedform changes before reaching a quasi-steady morphology characterized by a system of sand ripples along the slope and a large bar near the break point. Although the incoming wave characteristics are held fixed, the bottom morphology never reaches a strict steady state, but rather slowly changes due to the migration of ripples and bar transformation. When the wave characteristics are changed, the bedform adjusts to a new quasi-steady state after a suitable adjustment time. Studies conducted by placing model cobbles/mines on the evolving sandy bottom subjected to wave forcing show four distinct scenarios: (i) periodic cobble oscillations with zero mean displacement and small scour around the cobbles, (ii) mean onshore motion of relatively light cobbles, (iii) periodic burial of relatively heavy cobbles when their sizes are comparable to those of sand ripples, and (iv) the burial of relatively large cobbles under the bar, when the bar migrates due to changes of incoming waves. Quantitative data on the characteristics and dynamics of the bedform, including ripple-formation front propagating down the slope, ripple growth and drift, and flow around ripples, are presented. Physical explanations are provided for the observations.     

An investigation of the velocity field under regular and irregular waves over a sand beach

Near-bed horizontal (cross-shore) and vertical velocity measurements were acquired in a laboratory wave flume over a 1:8 sloping sand beach of finite depth. Data were acquired using a three-component acoustic Doppler velocimeter to measure the velocity field close to, but at a fixed distance from the bed. The near-bed velocity field is examined as close as 1.5 cm above a trough and crest of a ripple under three different types of wave forcing (Stokes waves, Stokes groups, and irregular waves). Although both horizontal and vertical velocity measurements were made, attention is focused primarily on the vertical velocity. The results clearly indicate that the measured near-bed vertical velocity (which was outside the wave-bottom boundary layer) is distinctly nonzero and not well predicted by linear theory. Spectral and bispectral analysis techniques indicate that the vertical velocity responds differently depending on the location over a ripple, and that ripple-induced effects on the velocity field are present as high as 4–8 cm above the bed (for vortex ripples with wavelengths on the order of 8 cm and amplitudes on the order of 2 cm). At greater heights above the bed, the observed wave-induced motion is adequately predicted by the linear theory.     

Barchan‐shaped ripple marks in a wave flume

Barchans, isolated crescent‐shaped bedforms, are believed to be formed under almost unidirectional wind or water ?ows and limited sand supply. The formation of barchan morphologies under the action of purely oscillatory wave motion has not yet been fully investigated. The present study attempted to form barchan topography in a wave ?ume and to compare this with barchans in the ?eld. Barchan morphologies of ripple size, called the barchan ripples, were generated from a ?at bed by the action of waves. The horn width, the distance between horn tips, of the barchan ripples increased linearly with an increase in the total length, the overall length projected on the centre line of the barchan, with a coef?cient common to barchan dunes in deserts. The ratio of horn length to horn width of the barchan ripples was smaller than that of barchan dunes, but similar to that of subaqueous barchans in the ?eld. The longer the wave period was, the larger the ratio of the body length to horn width became. Most subaqueous barchans formed under waves (in the laboratory) and unidirectional ?ows (in the ?eld) had blunter horns than subaerial barchans. The shape of the barchan ripples changed with wave period. The outer rim became rounder with increasing wave period. The relationship between the base area and the height of barchan morphologies seems to be linear, with a constant coef?cient for the scale from ripples to dunes. The barchan ripples showed a linear relationship between the height and the horn width, similar to that for barchan dunes. The migration speed of the barchan ripples was proportional to the cube of the ?ow velocity and was inversely proportional to height. The same relation with a different value of the coef?cient was obtained for barchan dunes. Copyright © 2004 John Wiley & Sons, Ltd.     

Controls on the aerodynamic roughness length and the grain‐size dependence of aeolian sediment transport

Estimates of the wind shear stress exerted on Earth's surface using the fully rough form of the law‐of‐the‐wall are a function of the aerodynamic roughness length, z₀. Accurate prediction of aeolian sediment transport rates, therefore, often requires accurate estimates of z₀. The value of z₀ is determined by the surface roughness and the saltation intensity, both of which can be highly dynamic. Here we report field measurements of z₀ values derived from velocity profiles measured over an evolving topography (i.e. sand ripples). The topography was measured by terrestrial laser scanning and the saltation intensity was measured using a disdrometer. By measuring the topographic evolution and saltation intensity simultaneously and using available formulae to estimate the topographic contribution to z₀, we isolated the contribution of saltation intensity to z₀ and document that this component dominates over the topographic component for all but the lowest shear velocities. Our measurements indicate that the increase in z₀ during periods of saltation is approximately one to two orders of magnitude greater than the increase attributed to microtopography (i.e. evolving sand ripples). Our results also reveal differences in transport as a function of grain size. Each grain‐size fraction exhibited a different dependence on shear velocity, with the saltation intensity of fine particles (diameters ranging from 0.125 to 0.25 mm) saturating and eventually decreasing at high shear velocities, which we interpret to be the result of a limitation in the supply of fine particles from the bed at high shear velocities due to bed armoring. Our findings improve knowledge of the controls on the aerodynamic roughness length and the grain‐size dependence of aeolian sediment transport. The results should contribute to the development of improved sediment transport and dust emission models. © 2018 John Wiley & Sons, Ltd.     

沙波纹生成的耦合映射格子模型

应用耦合映射格子模型模拟了三维、不对称、带有Y形连结的沙波纹的形成与演化过程。模拟显示了在风力作用下,平坦床面上的沙粒通过自组织形成沙波纹,其波长通过合并逐渐增长,而成熟后的沙波纹具有长而缓的迎风坡和短而陡的背风坡,并以确定的速度向下风方向迁移,这些特征很接近自然界中的沙波纹。此外,还模拟了双风向下沙波纹形态的变化,结果与Rubin等人关于横向、斜向床面形态形成条件的论断相符合。