海岸风沙流中不同粒径组沙粒的垂向分布模式

通过对河北昌黎黄金海岸沙丘风沙流的野外实地观测与室内风洞模拟实验数据的数值模拟,探讨了我国典型海岸沙丘风沙流中不同粒径组沙粒输沙量的垂向分布模式。结果表明,河北昌黎黄金海岸沙丘表面风沙流中不同粒径组沙粒输沙量的垂向分布特征并不一致,其中细沙和中沙符合典型的指数递减分布规律,但粗沙则为负幂函数分布。究其原因,主要与不同粒径组沙粒输沙量的分布高度及运动方式差异有关。在实际非均匀沙床面上,粗沙主要集中分布于沙丘表面4cm高度内湍流发育的近地表层,运动方式以蠕移为主,沙丘表面湍流的主导作用使其输沙量随高度的变化满足负幂律关系,但中沙和细沙则以跃移运动为主,跃移沙粒输沙量的垂向分布呈现指数递减特征。     

VERTICAL DISTRIBUTION OF GRAIN-SIZE PARAMETERS OF DRIFTING PARTICLES DURING SAND STORMS IN THE TAKLAMAKAN DESERT,CENTRAL ASIA

The vertical distributions of grain-size parameters of drifting sand flux during sand storms in the Taklamakan Desert, Central Asia, were investigated through field observation and laboratory wind-tunnel tests. The results show that grain-size parameters of the drifting sands near the ground surface are similar to those of the source deposits. Fine and very fine sands in the flux layer near the ground surface (from zero to 20 cm high) comprise more than 85% of the sample's weight, with a mean diameter of 0.09 mm. The distributions of transported particles are mainly negatively skewed, with good sorting. The creep population consists of poorly sorted coarser particles. In the saltation layer, particle size becomes finer as a power function of height from 4 to 14 cm. Sorting in the saltation layer is good and improves as a logarithmic function of height. The higher the transport layer, the more negatively skewed the particle distribution becomes as a logarithmic function of height. It becomes more leptokurtic as a power function of height in the transport layer from zero to 12 cm high. Grain-size parameters change irregularly at heights above 14 cm. Regional differences in the vertical distribution of grain-size compositions are affected by source deposits, wind speeds, and transport rates, but source has more influence on the sorting of the drifting material than wind velocity.     

风沙流中不同粒径组沙粒的 输沙量垂向分布实验研究

在非均匀沙床面上, 风沙流中不同粒径组沙粒的输沙量垂向分布, 是非均匀风沙运动研究的重点。研究首先通过风洞实验, 收集了风洞中垂线垂向输沙量分布沙样, 然后对集沙沙样进行了沙粒粒度分析实验, 实验分析结果得出了不同粒径组沙粒的输沙量垂向分布规律, 基于稳定平衡风沙跃移运动模型和本文实验结果, 最后数值模拟研究了不同粒径组沙粒输沙量垂向分布, 与沙粒起跳速度和角度之间的关系。本文研究结果得出, 在非均匀风沙流中, 粗粒径组沙粒垂向输沙量上部符合指数递减分布但近床面区偏离指数分布, 呈现为偏大型分布, 粗粒径组对应的沙粒起跳速度和角度分布均为指数函数; 细粒径组沙粒垂向输沙量在整 个高度上均符合指数递减规律, 细粒径组沙粒对应的起跳速度分布为指数函数, 起跳角度分布为高斯函数。沙粒的平均起跳速度, 在0.4u*~2.2u* 之间变化, 随着气流风速(u*) 和沙粒粒径的增加而减小。     

风沙流输沙通量垂向分布研究 ——以塔克拉玛干沙漠南缘流沙地表风沙流观测为例

 于2010年4—6月在塔克拉玛干沙漠南缘策勒绿洲外缘流沙地进行风速、风向和输沙量同步观测,利用观测到的28组数据对输沙通量的垂向分布进行研究。结果表明,在4.5～8 m·s-1的风速区间内,总输沙率随风速增长速率较慢,大于8  m·s-1时增长迅速。在较小风速时,4个高度层(1 cm、3 cm、9 cm和19 cm)的输沙比率变化不明显,当风速大于6.6 m·s-1时,1 cm和19 cm高度输沙比率的调整较为明显,前者减小而后者增加,中间两层(3 cm、9 cm)则维持一定波动。不同的蚀积条件对于输沙通量的垂向分布具有一定的影响,平均跃移高度则在一定程度上反映了地表沙粒的吹蚀和堆积过程。在4种常见的描述输沙通量垂向分布的模型中,Zingg的修正幂函数具有最好的拟合程度,但反映的风沙流信息最少,简单幂函数的参数包涵了最为丰富的风沙流信息但拟合程度最差,指数函数同时具有较好的拟合程度和参数意义,Fryrear的幂函数则与指数函数恰恰相反。在沙物质以极细沙为主的塔克拉玛干沙漠南缘,指数函数较之简单幂函数能更好地描述该区域输沙通量的垂向分布,验证并补充了风洞实验的结果。综合比较认为,指数函数能更好地描述塔克拉玛干沙漠南缘输沙通量的垂向分布。     

兰新高铁烟墩风区戈壁近地表风沙流跃移质垂直分布特性

跃移质作为风沙流的主体,其近地表垂直分布规律是风沙物理学的重要研究内容,对防沙工程具有重要的指导意义。受研究条件与观测仪器限制,戈壁特别是极端大风区近地表风沙流结构特性研究较为薄弱。利用多梯度风蚀传感器与阶梯式集沙仪对兰新高铁烟墩风区戈壁近地表风沙流跃移质的垂直分布特性进行了观测研究。结果表明：兰新高铁烟墩风区戈壁沙粒发生跃移运动的2 m高临界风速达12 m·s^-1;戈壁近地表风沙流具有明显的阵性特征,沙粒跃移发生的时间比例在50%以下,与平均风速成正相关关系,与风速脉动强度无显著相关关系;2 m高阵风7级风速下,戈壁跃移沙粒主要集中于地表50 cm范围内,近地表风沙流结构呈"象鼻效应",跃移质最大质量通量出现在地表2.5~5 cm高度处,沙粒最大跃移高度可达2 m,且沙粒跃移高度随2 m高风速的增加呈指数规律递增。因此,兰新高铁烟墩风区2 m高阻沙栅栏不足以完全阻截戈壁风沙流,是造成烟墩风区兰新高铁轨道积沙的重要原因之一。     

乌兰布和沙漠东北缘地表风沙流结构特征

在国家林业局磴口荒漠生态站长期监测的基础上,利用多种积沙仪,对乌兰布和沙漠东北缘流动沙丘、油蒿半固定沙丘、白刺半固定沙丘、油蒿固定沙丘、白刺固定沙丘5种典型下垫面近地面（0～100 cm）的风沙流输沙量进行了实地观测和对比分析。结果表明：（1） 输沙率（q）随高度（h）增加呈幂函数（q=ah^-b,R²≥0.8409）规律衰减,随风速（v）增大呈幂函数（q=av^b,R²≥0.9256）规律增加,42.8%～70.7%的输沙量分布在10 cm高度内,67.6%～90.0%的输沙量分布于30 cm高度内。当地表植被盖度达到40%以上时,输沙率下降至无植被覆盖地表输沙率的6.6%以下,可有效阻止地表风蚀。（2） 沙物质主要由粒径为50～250 μm的细沙和极细沙构成,各高度层风蚀物粒度组成服从单峰态分布,峰值在100～250 μm。随高度增加,风蚀物粒径范围趋于变窄,粒径趋于更细。（3） 起沙风多出现在WSW和NW方向,占全年起沙风的53.19%。风沙流中跃移输沙、蠕移输沙的空间分布在理论上应与风向频率分布基本一致,差异性主要由各方位风的强度和持续时间等因素导致。研究结果可为该区域防沙工程设计提供理论参考。     

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.     

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.     

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.     

Sonic anemometers in aeolian sediment transport research

Fast-response wind and turbulence instruments, including sonic anemometers, are used more and more in aeolian sediment transport research. These instruments give information on mean wind, but also on fluctuations and turbulent statistics, such as the uw covariance, which is a direct measure of Reynolds' stress (RS) and friction velocity. This paper discusses the interpretation of sonic anemometer data, the transformations needed to get proper results and turbulence spectra, and how they are influenced by instrument size, sampling frequency, and measurement height.Turbulence spectra characterize how much the different frequencies in the turbulent signals contribute to the variance of wind speed, or to the covariance of horizontal and vertical wind speed. They are important in determining the measurement strategy when working with fast-response instruments, such as sonic anemometers, and are useful for interpreting the measurement results. Choices on the type of sonic anemometer, observation height, sampling period, sampling frequency, and filtering can be made on the basis of expected high and low-frequency losses in turbulent signals, which are affected by those variables, as well as wind speed and atmospheric stability.Friction velocity and RS, important variables in aeolian sediment transport research, are very sensitive to tilt or slope errors. During a field experiment, the slope sensitivity of the RS was established as 9% per degree of slope, which is 1.5 times the value reported in literature on the basis of theoretical considerations. An important reason for the difference probably is the large influence of streamline curvature on turbulence statistics and thereby on the slope sensitivity of the RS. An error of 9% per degree of slope in the RS will translate into an error of approximately 4% per degree of slope in the calculated friction velocity.Space–time correlation of the horizontal wind speed is much larger than that of the vertical wind speed and the instantaneous RS. This largely explains why, in previous studies, a poor correlation was found between instantaneous RS measured at 3 m height and saltation flux near the surface, whereas the correlation between wind speed at some height and saltation flux was much better. Therefore, the poor correlation between RS away from the surface and saltation flux does not contradict that saltation flux is caused by RS.     

Dune sand transport as influenced by wind directions, speed and frequencies in the Ordos Plateau, China

The Ordos Plateau in China is a region with extensive wind erosion, severe desertification and various aeolian sand hazards. In order to determine aeolian sand transport in this region, the relationship between the sand transport rate and wind speed at 10 min frequencies was established by field observation in both the Qubqi Sand Desert and the Mu Us Sandy Land. Threshold wind speeds (2 m above the ground) for mobile, semi-fixed and fixed dune surfaces were estimated by field observations. The sand transport rate increased with the increase of the bare land ratio and near-bed wind speed. High-resolution meteorological 10 min average wind velocity data at 10 m above the ground were converted into velocity values at a height of 2 m to calculate sand transport potential based on three specific parameters decisive for sand transport: wind speed, duration and direction. The quantity of aeolian sand transported was calculated for various wind speed levels and directions, and the overall characteristics of sand transport on different dune surface types were determined by vector operation techniques. Sand transporting winds took place mainly in springtime. The prevailing wind directions were W, WNW and NW, with a frequency of more than 60% in total, and sand transport in these directions made up more than 70% of the total transport, corresponding to a general southeastward encroachment of aeolian sand in the study area. The relationship between wind frequency and speed can be expressed by a power function. High magnitude strong winds had a low frequency, but they played a dominant role in aeolian sand transport.     

新月形沙丘表面100cm高度内风沙流输沙量垂直分布函数分段拟合

为研究新月形沙丘表面不同层位风沙流输沙量的垂直分布函数,实测了塔克拉玛干沙漠腹地典型新月形沙丘表面100 cm高度内(以1 cm分隔)的输沙量。分段拟合分析表明:新月形沙丘迎风坡脚输沙量垂直分布规律不完全服从指数函数,出现与戈壁风沙流结构特征相似的"象鼻效应",在0~3 cm区间内输沙量逐渐增大,3cm以上输沙量随高度呈指数函数衰减;沙丘顶部0~10 cm区间输沙量随高度呈指数函数衰减,10 cm以上呈二次函数衰减;沙丘左翼端输沙量随高度呈幂函数分布,沙丘右翼端0~20 cm内以指数函数衰减,20 cm以上呈三次函数衰减;沙丘背风坡脚风沙流输沙量在0~60 cm和60 cm以上分别呈不同形式的三次函数分布。     

The flux profile of a blowing sand cloud: a wind tunnel investigation

The flux profile of a blowing sand cloud, or the variation of blown sand flux with height, is the reflection of blown sand particles that move in different trajectories, and also the basis for checking drifting sand. Here we report the wind tunnel results of systematic tests of the flux profiles of different sized sands at different free-stream wind velocities. The results reveal that within the 60-cm near-surface layer, the decay of blown sand flux with height can be expressed by an exponential function: q_h=aexp(−h/b), where, q_h is the blown sand transport rate at height h, a and b are parameters that vary with wind velocity and sand size. The significance of coefficient a and b in the function is defined: a represents the transport rate in true creep and b implies the relative decay rate with height of the blown sand transport rate. The true creep fraction, the ratio of the sand transported on the surface (h=0) to the total transport varies widely, decreasing with both sand size and wind speed. The flux profiles are converted to straight lines by plotting sand transport rate, q_h, on a log-scale. The slope of the straight lines that represents the relative decay rate with height of sand transport rate decreases with an increase in free-stream wind velocity and sand grain size, implying that relatively more of the blown sand is transported to greater heights as grain size and wind speed increase. The average saltating height represented by the height where 50% of the cumulative flux percentage occurs increases with both wind speed and grain size, implying that saltation becomes more intense as grain size and/or wind velocity increase.     

塔克拉玛干沙漠南缘风沙跃移运动研究——以策勒为例

 利用Sensit传感器、BSNE集沙仪、全方位跃移集沙仪结合2 m气象塔,于2010年7月5日至8月24日在策勒绿洲-荒漠过渡带开展了风沙跃移运动强化观测试验,分析了该区域风沙跃移运动的部分特征。结果表明:①试验点2 m高度临界跃移起动1 h平均风速为5.0 m·s-1;跃移颗粒数与风速呈现较好的幂函数关系;②跃移输沙量的方向分布以偏西方向为主,主要集中在SW、WSW、W、WNW、NW5个方位,共计占到80.9%;③跃移运动在一日之内的任何时间段都可能发生,发生频率较高的时段集中在白天的11:00—19:00;④观测期间0～40 cm高度范围的跃移输沙通量为199.4 kg·m-1。     

野外近地表风沙流脉动特征分析

基于风沙流实时测量系统对流沙地表的输沙强度以及不同高度风速进行实时测量所获得的数据,分析野外实际风场中的风沙流脉动特征。结果表明:沙粒的存在削弱平均流场,但对风速的高阶统计矩和分布形式影响不大,不同高度、不同来流平均风速条件下,风沙流中风速仍近似符合高斯分布;瞬时输沙强度服从指数分布,滞后风速约1 s,脉动强烈;30 min输沙强度平均值和脉动标准差均随平均风速的增大而呈幂函数增加,但随风速脉动强度的变化呈先增大后减小的趋势。这意味着脉动风速对输沙率具有重要影响,需要在输沙率预测公式中予以考虑。     

海岸湿沙表面风沙传输特征的风洞实验研究

通过湿沙表面风沙传输的风洞实验,研究了湿润海岸风沙流的垂直结构、输沙率随风速与表层湿度的变化规律。输沙测量使用60 cm高直立式积沙仪,湿度（M）为沙面表层1 mm厚的重量湿度值。结果表明,湿沙表面的输沙量和高度呈指数关系。一般,湿度增大,整体输沙率降低,高湿度床面的沙粒有相对更大的比例被传输到更高的位置。比起低湿度（M＜0.587%）沙粒,高湿度（0.587%＜M＜1.448%）沙粒的垂直运动对湿度变化的影响更加敏感,尤其是在跃移层的底部,当M＞1.448%时,输沙率已经很低,小于0.99 g·cm-1·s-1。伴随湿度0.587%和1.448%的过渡,风沙流垂向分布被分为3个不同坡度的区域,曲线坡度反映了沙粒间不同水分存在形式的影响差异。对于跃移沙粒,高湿度表面（M＞1.448%）仅起到了一个传输平台的作用;当表面变干到某种程度（M=0.587%）之前,表面湿度是跃移运动的主要控制因子,然后风速才重新开始影响输沙。     

The distribution of velocity and energy of saltating sand grains in a wind tunnel

Sand transport by wind is a special case of two-phase flow of gas and solids, with saltating grains accounting for about 75% of the transport rate. This form of flow is not only the main external agent moulding aeolian landforms but also the motive force responsible for transport, sorting and deposition of aeolian sediments. High-speed multiflash photography is an effective method of studying the distribution of velocity and energy of saltating grains within the boundary layer of wind tunnel. The experimental wind shear velocities were set at 0.63, 0.64, 0.74 and 0.81 ms⁻¹. The statistical study of the results showed that there is a power function relation between mean velocity and height of saltating grains. As the height is divided into 0.5-cm intervals, the sand grain velocities at various levels are consistent with the Pearson VII distribution pattern. The variations in kinetic energy and total energy of sand grains with height accord with the pulse peak modified with power term (Pulsepow) law; the maximum values occur at heights of 6 cm or so and tend to shift upward with increasing wind velocity.     

河北坝上高原土壤风蚀物垂直分布的初步研究

根据河北坝上高原农田土壤风蚀物的观测与采样分析结果,风蚀物含量随高度的增加,在0~20 cm高度内以指数函数规律递减,反映了以跃移质为主的风沙流结构;在20~100 cm高度内以幂函数规律递减,反映的是悬移质为主的风沙流结构。在风蚀物粒度组成中,随高度增加,砂级颗粒含量减少而粉砂及粘土含量增多;随风速与输沙率的增加,砂级颗粒含量增多而粉砂及粘土含量减少。对风蚀物各粒级含量与高度的相关分析表明:易以迁移形式运动的粒径0.25~0.1 mm土粒百分量以幂函数形式向上递减;沉降速度较低的粒径<0.1 mm土粒百分含量成幂函数形式向上递增。此外,风蚀物平均粒径随高度变细且在50 cm处细于0.1 mm,因此,土壤颗粒成为悬移物质而上升到50 cm以上高度时可能搬运到较远的地区。     

海滩养护影响下的海岸风沙堆积——以福建平潭岛龙凤头海滩养护工程为例

海滩养护工程实施后,在保护海岸抵御侵蚀的同时往往还产生其他的环境问题。尤其在强风区海岸,海滩养护造成了滩面环境的改变,使滩面风沙搬运和沉积特征都发生了显著变化。对福建平潭岛龙凤头海滩养护工程实施前后的风沙沉积对比研究表明,养护后海滩的干滩宽度、滩面沉积物粒径、表层湿度、海滩高程等环境要素变化对滩面风沙搬运与沉积过程都产生了重要的影响。主要结论为:(1)养护海滩干滩滩肩高程的提高增加了滩面风速,增大了滩面风沙输沙率;(2)养护海滩干滩宽度的增加既扩大了风区长度也为风沙搬运提供了充足的物源,且不受潮汐过程影响,提供了持续的风沙物源供给,增加了风沙作用时间;(3)养护海滩剖面形态的变化使得海滩滩面沉积物含水率减小,降低风沙起动风速,增强了滩面风沙搬运。     

GRAIN-SIZE DISTRIBUTIONS OF WIND-ERODED MATERIAL ABOVE A FLAT BARE SOIL

Grain-size distributions and fluxes of each size fraction of aeolian materials in the surface layer are basic to understanding and predicting soil loss, air quality, and visibility, both in time and in space. Vertical distributions of eroded soil materials in a 612-cm-deep flow layer above an Amarillo fine sandy loam soil in Big Spring, Texas were measured during sand-dust storms in the 1995 wind-erosion season. All 942 samples in 33 groups of measurements show that the airborne particles in the flow layer below 40 cm are bimodal in distribution, which is similar to the size distribution of the original soil. In the flow layer above 80 cm, the airborne material is distributed unimodally and consists of silt and very fine sand. The distribution of the mean diameter is a discontinuous function of height. Mean diameters vary by an asymmetric double sigmoid function of height. Size distribution of airborne particles in the 0- to 35-cm flow layer is controlled by saltation mechanisms, becoming larger with an increase in height, and moderately sorted, positively skewed, and platykurtic. In the 35- to 50-cm flow layer, particle size decreases very rapidly with increased height, becoming poorly sorted, symmetrical, or platykurtic, and thus uniformly distributed. This represents a transition zone from saltation to suspension. Grain-size parameters in the flow layer above 50 to 80 cm demonstrate a wavelike pattern with a moderately well to well-sorted, negatively skewed, and leptokurtic distribution. The weight percent of the silt content and phi₉₅ diameter increase by a natural logarithmic polynomial function of height. Inclusive standard deviations vary with height by an asymmetric logistic function. Total distributions of the skewness and kurtosis are wavelike. The fluctuation of grain-size parameters of eroded soil materials in the surface layer is related to saltation mechanisms and the turbulence of air flow.