风成沙波纹脊线提取与应用计算

风成沙波纹是沙质地表在风力作用下形成的最小地貌单元,对研究风沙的起动过程和运动过程极其重要,同时,沙波纹形态特征的研究为了解大尺度范围风沙地貌形态演变提供理论基础。然而由于风成沙波纹尺度较小,形成速度较快,导致对其形态特征的观测比较困难。近些年,随着计算机图形学的迅猛发展,数字图像处理方法得到了较大发展,使得测量和计算更加便捷。本文基于高清相机拍摄的风成沙波纹图像,借助于Matlab平台,采用数字图像处理技术,对沙波纹脊线进行提取,并应用于沙波纹形态参数计算。本文选取腾格里沙漠东南缘的沙波纹进行验证分析,得到风成沙波纹波长的正态分布规律,沙波纹的波长随时间逐渐增长,在40 min的时间范围内,波长由不足1 cm逐步发展到接近10 cm。最后,采用数字图像处理方法,计算了沙波纹脊线长度和波长。采用数字图像处理方法,波长等数据更易于获取和统计,数据采集效率大大提高,为风成沙波纹的研究提供了新的方法。     

农田近地表风沙流风程效应变化特征研究

风沙流的风程效应研究是定量获取风沙流沿程变化的核心和难点,风程效应是指输沙率随沙床表面或地块长度的增加而不断增大,而后趋于稳定的变化特征,饱和输沙率（f_max）和饱和路径长度（L_sat）是风程效应的重要参数。采用自动连续称重式集沙仪,以河北坝上地区康保县境内典型旱作农田为研究对象,观测了2017、2018年和2021年内4次典型风蚀事件,分析近地表5 cm高度风沙流的风程效应在5 min时间尺度下的变化特征。结果表明：（1） 近地表输沙通量随风程距离的增大而增大。（2） 4次风蚀事件中L_sat的变化范围在11~280 m之间,并存在明显差异,其变化与风速无关。（3） 近地表风沙流的f_max与风速（U）呈幂函数关系。（4） 风程效应的变化特征与地表可蚀性因子、地表微地貌变化有着紧密联系,未来应对不同的土壤类型和质地农田的风程效应进行深入研究。     

风成沙波纹数学模型综述

风成沙波纹是风力作用下地表形成的最基本的地貌类型。近年来,随着计算机技术的迅猛发展,风成沙波纹数值模拟研究越来越引起重视。本文简要介绍了风成沙波纹数学建模的理论基础,并将已有的风成沙波纹数学模型划分为3类,着重介绍了几个经典的风成沙波纹数学模型,详细列举了这些模型的假设条件和数学表达,分析了他们的优缺点。目前,风成沙波纹数学建模的两个基本理论依据为:一个是地表坡度变化导致的蠕移沙通量的变化;另一个是跃移长度随沙床面高度变化的规律,但目前关于此理论的研究并不多。最后,总结和展望了这一领域的研究方向和前景。     

稳态风沙流中瞬态输沙特征

风沙流中沙粒运动在来流风速不变时也会表现出非稳态特征。在风洞内利用粒子图像测速系统（PIV）测量了风沙运动的时间序列，并基于PIV测量技术提出风沙流中沙粒平均直径、数密度、平均水平速度和输沙通量等参数在某一时刻的计算方法，其中输沙通量的计算考虑沙粒大小垂向分布的影响。结果表明：来流风速不变时，沙粒平均直径、数密度、平均水平速度和输沙通量随时间具有明显的波动性；沙粒平均直径和平均水平速度的标准偏差一般随高度增加而增加，沙粒数密度和输沙通量标准偏差随高度增大而减小；这些参数的相对标准偏差均随高度增加而增大。     

应用近景摄影法研究沙纹的移动

近景摄影能够在无干扰的条件下,精确地测定正在发育阶段的风成沙纹,以波动方式整体移动的平均移动速度V_R(cm·min^-1),和相应的沙纹波长L(cm)波高Δh(mm)及其平均值。同时,求得沙纹移动速度V_R与局地风速V_L之间的实验关系式为V_R=158×(V_L-5.5)^0.67,沙纹移动是缓慢的,其数量级变化于10^-1~10¹cm·min^-1之间。沙纹波长与波高与直接测定值相当吻合。     

湍流度随高度的变化及其对城市宏观地形的依赖

分别构建广州主建成区垂直比例尺为1﹕2 000、1﹕1 000和1﹕500的3个建筑物模型，利用大型边界层风洞，在西北和东南两风向下，基于中性流模拟分析了复杂城市地形下湍流度随高度的变化及其对宏观地形的依赖。结果表明：风廓线指数α与不同高度的湍流度之间的关系密切，利用现有模型，根据4类粗糙度边界层和不同垂直比例尺，可确定相应的湍流度随高度变化模型的主要系数，预测精度高。城市地形下最大湍流度面发育在0~0.2 h之间狭窄的范围内。用湍流度形态指数β来表征湍流度随高度的变化，无论城市屋脊还是平坦地形，随着风程区的延伸，廓线的指数α升高，湍流度形态指数β降低。表明同一高度湍流度值具有由迎风区、丘顶区向背风区增高，沿风程逐渐增大的规律，对地形部位和风程的依赖性强，与来流翻越简单地形时的特征一致。     

人工卵石床面风沙流粒度特征

沉积物粒径分布对风沙沉积物的起动、输送和沉降过程十分重要。不同粒径的风沙沉积物空气动力学特征不同,从而导致其起动机理、输送过程和沉降模式不同。因此,粒度可以作为风沙沉积过程的指示器。戈壁是西北地区的一种重要地貌景观,戈壁风沙流是风沙物理研究的内容之一,但目前研究较少,特别是对戈壁风沙流粒度分布的研究,几乎没有报道。本文利用人工卵石床面模拟戈壁地表,对0.25、0.5、1 m和2 m高度的风沙流粒度特征进行观测,结果表明：沙源、地表状况、风程效应等影响沙粒粒径随高度的变化。沙源近,平均粒径随高度先减小后增加;沙源远,平均粒径随高度增加减小。戈壁表面的不同区域,沙粒粒径级配不同：随高度增加,粗砂,中砂含量降低,细砂、极细砂和粉砂含量增加;随距离增加,中砂含量降低,细砂、极细砂和粉砂含量增加。平均粒径越大,分选越好,偏度越趋于正偏,峰度越趋于变宽。偏度、峰度随分选系数增加而增加。峰度随偏度增加而增加。     

科尔沁沙地大型沙波纹的初步研究

通过实地考察、观测和采样,利用Google Earth遥感影像和激光粒度分析仪,研究了科尔沁沙地大型沙波纹（large-scale ripples, LSR）的空间分布、单体和群体的基本形态和粒度特征,探讨了LSR粗颗粒的来源、它同普通沙波纹和普通风成沙在颗粒组成、形态特征和内部沉积结构方面的区别。结果表明：（1）科尔沁沙地的LSR主要分布于翁牛特旗中部和北部地区。（2）LSR空间单元的空间形态具有片状、斑块状和条带状3种类型,分别发育于宽阔的丘间地、沙丘中上部、槽形低地3类地形部位。（3）LSR平均长度为6.32 m,总体走向为东北—西南,平均波长为1.68 m,空间分异较为明显;单体LSR的前后坡不对称。（4）LSR的颗粒为中砂-粗砂粒级,其中粗颗粒主要来自旧河道的河流冲积层、下伏Q3河湖相地层、剥蚀残山的风化壳和山麓洪积物。（5）LSR同普通沙波纹在外观、物质组成、几何形态和内部沉积结构方面有明显差异。该项研究将有助于促进风沙地貌学的理论发展和实践。     

海岸新月形沙丘移动与形态变化的典型研究

依据2006~2008年连续3 a共9次采用RTK GPS技术与测量方法对一个典型海岸新月形沙丘形态的高精度测量数据,分析了海岸新月形沙丘的移动方向、方式、速度以及形态变化特点。结果表明,海岸新月形沙丘具有缓慢、向陆往复式前进的移动特点,形态变化则具有随季节增减变化中高度、宽度、长度、断面面积与体积增加的加积特征,究其原因是区域风况、海岸地表覆被、沙丘形态及人类活动等共同作用的结果。     

毛乌素沙地不同下垫面的风沙运动特征

 基于野外实地观测,研究了毛乌素沙地不同下垫面的风沙运动特征。平坦流沙地上的风速随高度呈对数分布; 沙丘上的风速梯度在高度上的变化随植被盖度增大而增强; 相近风速下,0~60 cm高度范围内的输沙强度从大到小依次为流动沙丘丘顶＞平坦流沙地＞半流动沙丘丘顶＞裸露黄土梁地＞黄土梁翻耕地＞玉米留茬地＞半固定沙丘丘顶。同一植被盖度下,沙丘顶部输沙率随风速增加而增大,粗糙度随风速增加而减小。当风速小于起沙风速时,粗糙度随植被盖度的增加而增大; 当风速大于起沙风速时,平坦流沙地和沙丘顶部的粗糙度随风速增大有增加趋势,产生这种现象的可能原因是跃移层沙物质对气流的阻力随风速增大而增大。根据各类沙源地的起沙强度和总面积,风沙灾害防治的重点是植被盖度较小的沙地和黄土梁地。     

Velocity profile of a sand cloud blowing over a gravel surface

Particle dynamic analyzer (PDA) measurement technology was used to study the turbulent characteristics and the variation with height of the mean horizontal (in the downwind direction) and vertical (in the upward direction) particle velocity of a sand cloud blowing over a gravel surface. The results show that the mean horizontal particle velocity of the cloud increases with height, while the mean vertical velocity decreases with height. The variation of the mean horizontal velocity with height is, to some extent, similar to the wind profile that increases logarithmically with height in the turbulent boundary layer. The variation of the mean vertical velocity with height is much more complex than that of the mean horizontal velocity. The increase of the resultant mean velocity with height can be expressed by a modified power function. Particle turbulence in the downwind direction decreases with height, while that in the vertical direction is complex. For fine sands (0.2–0.3 mm and 0.3–0.4 mm), there is a tendency for the particle turbulence to increase with height. In the very near-surface layer (<4 mm), the movement of blown sand particles is very complex due to the rebound of particles on the bed and the interparticle collisions in the air. Wind starts to accelerate particle movement about 4 mm from the surface. The initial rebound on the bed and the interparticle collisions in the air have a profound effect on particle movement below that height, where particle concentration is very high and wind velocity is very low.     

The blown sand flux over a sandy surface: a wind tunnel investigation on the fetch effect

Detailed wind tunnel tests were conducted to examine the fetch effect of a sandy surface on a sand cloud blowing over it. The results suggest that the fetch length of a sandy surface has a significant effect on both the vertical flux profile and total horizontal flux. The sand flux over a sandy surface increases with height in the very near surface layer, but then decays exponentially. In agreement with the widely accepted conclusion, the decay function can be expressed by q=aexp(−h/b), where q is the sand flux at height h. Coefficient a that tends to increase with wind speed implies the influence of wind, while coefficient b that defines the relative decay rate shows the influence of both the fetch and wind. The relative decay rate increases with fetch when the fetch length is short, then becomes constant when the fetch reaches a certain length. The threshold fetch length over which the relative decay rate keeps constant increases with wind speed. The average saltation height generally increases with fetch. Both the relative decay rate and average saltation height show that the fetch effect on the flux profile becomes more significant when the wind speed increases. The total sand transport equation for the total fetch can be expressed by Q=C(1−U_t/U)²U³(ρ/g), where Q is the total sand transport rate, U and U_t are the wind velocity and threshold wind velocity at the centerline height of the wind tunnel, respectively, g is gravitational acceleration, ρ is the density of air, and C is a proportionality coefficient that increases with the fetch length, implying that the total sand flux increases with the fetch length.     

Effects of fetch and surface texture on aeolian sand transport on two nourished beaches

Aeolian sand transport on two nourished beaches was related to the fetch of the wind over the beach sand and to surface characteristics. Meteorological and hydrological conditions were recorded for 2 months. The fetch of wind over beach sand was estimated from wind direction, water level, wave height and beach topography. Aeolian sand transport was determined with sand traps. Sediment flux was found to increase with fetch, although this relation was especially affected by the variability in surface characteristics. On one of the beaches sediment supply was limited as a result of shells, forming a lag deposit.     

Controls on the height and spacing of eolian ripples and transverse dunes: A numerical modeling investigation

Ripples and transverse dunes in areas of abundant sand supply increase in height and spacing as a function of time, grain size, and excess shear velocity. How and why each of these factors influence ripple and transverse dune size, however, is not precisely known. In this paper, the controls on the height and spacing of ripples and transverse dunes in areas of abundant sand supply are investigated using a numerical model for the formation of eolian bedforms from an initially flat surface. This bedform evolution model combines the basic elements of Werner's [Werner, B.T., 1995. Eolian dunes: Computer simulations and attractor interpretation. Geology 23, 1107–1110.] cellular automaton model of dune formation with a model for boundary layer flow over complex topography. Particular attention is paid to the relationship between bed shear stress and slope on the windward (stoss) side of evolving bedforms. Nonlinear boundary layer model results indicate that bed shear stresses on stoss slopes increase with increasing slope angle up to approximately 20°, then decrease with increasing slope angle as backpressure effects become limiting. In the bedform evolution model, the linear boundary layer flow model of Jackson and Hunt [Jackson, P.S., Hunt, J.C.R., 1975. Turbulent wind flow over a low hill. Quarterly Journal of the Royal Meteorological Society 101, 929–955.], generalized to 3D, is modified to include the nonlinear relationship between bed shear stress and slope. Bed shear stresses predicted by the modified Jackson and Hunt flow model are then used to predict rates of erosion and deposition iteratively through time within a mass-conservative framework similar to Werner [Werner, B.T., 1995. Eolian dunes: Computer simulations and attractor interpretation. Geology 23, 1107–1110.]. Beginning with a flat bed, the model forms ripples that grow in height and spacing until a dynamic steady-state condition is achieved in which bedforms migrate downwind without further growth. The steady-state ripple spacing predicted by this model is approximately 3000 times greater than the aerodynamic roughness length of the initially flat surface, which is a function of grain size and excess shear velocity. Once steady-state ripples form, they become the dominant roughness element of the surface. The increase in roughness associated with ripple formation triggers the same bedform instability that created ripples, causing dunes to form at a larger scale. In this way, the numerical model of this paper suggests that ripples and dunes are genetically linked. Transverse dunes in this model have a steady-state height and spacing that is controlled by the effective roughness length of the rippled surface, which is shown to be on the order of 500 times greater than the original roughness length, but varies significantly with the details of ripple morphology. The model predictions for ripple and dune spacing and their controlling variables are consistent with field measurements from the published literature. The model of this paper provides a preliminary process-based understanding of the granulometric control of ripples and dunes in areas of abundant sand supply and unidirectional prevailing winds, and it argues for a genetic linkage between ripples and dunes via a scaling relationship between eolian bedform size and the aerodynamic roughness length.     

Experimental investigation of the concentration profile of a blowing sand cloud

Detailed wind tunnel tests were carried out to establish the mean downwind velocity and transport rate of different-sized loose dry sand at different free-stream wind velocities and heights, as well as to investigate the vertical variation in the concentration of blowing sand in a cloud. Particle dynamic analyzer (PDA) technology was used to measure the vertical variation in mean downwind velocity of a sand cloud in a wind tunnel. The results reveal that within the near-surface layer, the decay of blown sand flux with height can be expressed using an exponential function. In general, the mean downwind velocity increases with height and free-stream wind velocity, but decreases with grain size. The vertical variation in mean downwind velocity can be expressed by a power function. The concentration profile of sand within the saltation layer, calculated according to its flux profile and mean downwind profile, can be expressed using the exponential function: c_z=ae^−bz, where c_z is the blown sand concentration at height z, and a and bare parameters changing regularly with wind velocity and sand size. The concentration profiles are converted to rays of straight lines by plotting logarithmic concentration values against height. The slope of the straight lines, representing the relative decay rate of concentration with height, decreases with an increase in free-stream wind velocity and grain size, implying that more blown sand is transported to greater heights as grain size and wind speed increase.     

Grain size and transport characteristics of non-uniform sand in aeolian saltation

Size frequency distributions of sediment particles in a wind tunnel containing a bed of non-uniform sand are investigated by re-interpreting existing experimental data using particle-size analysis. Each particle sample is classified into one of eight groups according to its size grading. The analysis reveals that the modal shape of the particle-size frequency distributions of the saltating sand at different elevations or longitudinal distances is similar to that of the mixed sand in the bed once the boundary layer is fully developed. The standard deviation of the grain-size frequency distribution increases with increasing elevation above the bed then stays constant, whereas its skewness decreases. The mean grain size decays exponentially with elevation. The aeolian sand mass flux is determined for each size grading at different vertical and horizontal measurement locations. The vertical profile of aeolian horizontal mass flux depends on the size grading. The distribution of the sand transport rate according to the mean grain size in each grading fits the normal distribution. A parameter w_i is defined to reflect the likelihood of saltation for sand particles of the i-th size grading, and the mean sand size corresponding to the maximum value of w_i is found to be 0.2 mm. In addition, wind velocity strongly influences the magnitudes of the particle-size distribution and the sand mass flux distribution in both vertical and longitudinal directions.     

Aeolian sand transport along beaches

Storm‐generated ephemeral transverse sand ridges were observed developing along the beach fronting Sir Richard Peninsula, South Australia during 25 March 1984. The ridges displayed a mean height of 10 cm and a wavelength of 12 m; their breadth was approximately 2 m, and length varied with beach width but ranged up to 40 m over 10 km of coastline. The steeper sides of the ridges faced upwind due to erosion after initial ridge development. Damp sand in the swales inhibited wind scour and restricted sand supply, but provided a firm substrate over which the sand could saltate. Approximately 5000 m³ of sand were incorporated into the ridges by westerly winds blowing at velocities between 45 and 69 km/hr. This observation emphasises the role of alongshore winds in transporting beach sediments and developing essentially ephemeral forms, which might, nevertheless, be preserved in the geological record. The significance of these forms varies with coastal orientation and local wind regimes.