黄土高原半干旱区非均一下垫面粗糙度分析

利用2007年4月17日-2008年4月16日兰州大学半干旱气候与环境观测站边界层气象塔的风速、 风向、 温度、 气压、 湿度等观测资料, 采用经典的廓线法和风速、 风向标准差法, 分别计算了中性大气层结下观测站下垫面粗糙度长度, 并得到了具有黄土高原地理特征的地表粗糙度及其时空变化特征。计算结果表明, 季节变化对粗糙度的影响幅度可达0.159 m, 空间非均一性对粗糙度的影响幅度可达0.155 m。测站附近粗糙度春季为0.017 m, 夏季为0.062 m, 秋季为0.065 m, 冬季为0.018 m。测站西北方向上游粗糙度春季为0.17 m, 夏季为0.22 m, 秋季为0.34 m, 冬季为0.05 m。测站东南方向上游粗糙度春季为0.11 m, 夏季为0.17 m, 秋季为0.19 m, 冬季为0.05 m。该站下垫面粗糙度计算宜选用风速为6±1.5 m·s-1, 风向变化30°范围内的数据。     

Effect of bed surface roughness on longitudinal dispersion in artificial open channels

In several empirical and modelling studies on river hydraulics, dispersion was negatively correlated to surface roughness. In this study, it was aimed to investigate the influence of surface roughness on longitudinal dispersion under controlled conditions. In artificial flow channels with a length of 104 m, tracer experiments with variations in channel bed material were performed. By use of measured tracer breakthrough curves, average flow velocity, mean longitudinal dispersion, and mean longitudinal dispersivity were calculated. Longitudinal dispersion coefficients ranged from 0·018 m² s^?1 in channels with smooth bed surface up to 0·209 m² s^?1 in channels with coarse gravel as bed material. Longitudinal dispersion was linearly related to mean flow velocity. Accordingly, longitudinal dispersivities ranged between 0·152 ± 0·017 m in channels with smooth bed surface and 0·584 ± 0·015 m in identical channels with a coarse gravel substrate. Grain size and surface roughness of the channel bed were found to correlate positively to longitudinal dispersion. This finding contradicts several existing relations between surface roughness and dispersion. Future studies should include further variation in surface roughness to derive a better‐founded empirical equation forecasting longitudinal dispersion from surface roughness. Copyright © 2011 John Wiley & Sons, Ltd.     

A real-time, event-triggered storm surge forecasting system for the state of North Carolina

A new real-time, event-triggered storm surge prediction system has been developed for the State of North Carolina to assist emergency managers, policy-makers and other government officials with evacuation planning, decision-making and resource deployment during tropical storm landfall and flood inundation events. The North Carolina Forecast System (NCFS) was designed and built to provide a rapid response assessment of hurricane threat, accomplished by driving a high-resolution, two-dimensional, depth-integrated version of the ADCIRC (Advanced Circulation) coastal ocean model with winds from a synthetic asymmetric gradient wind vortex. These parametric winds, calculated at exact finite-element mesh node locations and directly coupled to the ocean model at every time step, are generated from National Hurricane Center (NHC) forecast advisories the moment they are inserted into the real-time weather data stream, maximizing the number of hours of forecast utility. Tidal harmonic constituents are prescribed at the open water boundaries and applied as tidal potentials in the interior of the ocean model domain. A directional surface roughness parameterization that modulates the wind speed at a given location based on the types of land cover encountered upwind, a forest canopy sheltering effect, and a spatially varying distribution of Manning’s–n friction coefficient used for computing the bottom/channel bed friction are also included in the storm surge model. Comparisons of the simulated wind speeds and phases against their real meteorological counterparts, of model elevations against actual sea surface elevations measured by NOAA tide gauges along the NC coast, and of simulated depth-averaged current velocities against Acoustic Doppler Current Profiler (ADCP) data, indicate that this new system produces remarkably realistic predictions of winds and storm surge.     

A traversing system to measure bottom boundary layer hydraulic properties

This study describes a new convenient and robust system developed to measure benthic boundary layer properties, with emphasis placed on the determination of bed shear stress and roughness height distribution within estuarine systems by using velocity measurements. This system consisted of a remotely operated motorised traverser that allowed a single ADV to collect data between 0 and 1 m above the bed. As a case study, we applied the proposed traversing system to investigate bottom boundary layer (BBL) hydraulic properties within Coombabah Creek, Queensland, Australia. Four commonly-employed techniques: (1) Log-Profile (LP); (2) Reynolds stress (RS); (3) Turbulent Kinetic Energy (TKE); and (4) Inertial Dissipation (ID) used to estimate bed shear stresses from velocity measurements were compared. Bed shear stresses estimated with these four methods agreed reasonably well; of these, the LP method was found to be most useful and reliable. Additionally, the LP method permits the calculation of roughness height, which the other three methods do not. An average value of bed shear stress of 0.46 N/m², roughness height of 4.3 mm, and drag coefficient of 0.0054 were observed within Coombabah Creek. Results are consistent with that reported for several other silty bed estuaries.     

基于卫星遥感数据改进区域尺度的显热通量模型

以西北半干旱地区的内蒙古农牧交错地带锡林郭勒盟为研究区,以不同退化程度的草地为研究对象,通过叶面积指数计算各像元的动力传输粗糙度长度,以改进显热通量算法;并结合植被／土壤组分温度分别计算基于亚像元的植被／土壤的显热通量。由敏感性分析和模型验证表明,改进的显热通量算法能够提高区域尺度的能量平衡计算精度。     

Numerical modeling investigation on turbulent oscillatory flow over a plane rough bed composed by randomly arrayed particles

A three-dimensional numerical model is established to simulate the turbulent oscillatory boundary layer over a fixed and rough bed composed by randomly arrayed solid spheres based on the lattice Boltzmann method and the large eddy simulation model.The equivalent roughness height,the location of the theoretical bed and the time variation of the friction velocity are investigated using the log-fit method.The time series of turbulent intensity and Reynolds stress are also investigated.The equivalent roughness height of cases with Reynolds numbers of 1×10~4–6×10~4 is approximately 2.81 d(grain size).The time variation of the friction velocity in an oscillatory cycle exhibits sinusoidal-like behavior.The friction factor depends on the relative roughness in the rough turbulent regime,and the pattern of solid particles arrayed as the rough bed in the numerical simulations has no obvious effect on the friction factor.     

机载LiDAR地形特征信息快速提取

为更有效地获取地形特征信息,提出一种机载LiDAR地形特征信息快速提取算法。首先,通过构建二次曲面拟合模型,建立实测LiDAR地形数据与拟合曲面的几何规则;然后,采用LM算法迭代参数寻优,获得最优化结果下的地形拟合参数,计算拟合时间及拟合精度;最后,以地形拟合模型为基础,进行地形特征信息的快速提取。通过机载LiDAR实测数据验证,当最优搜索半径为2 m时,地形曲面的拟合时间仅为0.02 s,RMSE仅为5.09 cm。该算法保证了地形特征信息提取效率和精度,能够有效满足机载LiDAR科学研究和工程应用的技术需求。     

Smoothing and re-roughening processes: The geometric evolution of a single fault zone

The geometry of a fault zone exerts a major control on earthquake rupture processes and source parameters. Observations previously compiled from multiple faults suggest that fault surface shape evolves with displacement, but the specific processes driving the evolution of fault geometry within a single fault zone are not well understood. Here, we characterize the deformation history and geometry of an extraordinarily well-exposed fault using maps of cross-sectional exposures constructed with the Structure from Motion photogrammetric method. The La Quinta Fault, located in southern California, experienced at least three phases of deformation. Multiple layers of ultracataclasite formed during the most recent phase. Crosscutting relations between the layers define the evolution of the structures and demonstrate that new layers formed successively during the deformation history. Wear processes such as grain plucking from one layer into a younger layer and truncation of asperities at layer edges indicate that the layers were slip zones and the contacts between them slip surfaces. Slip surfaces that were not reactivated or modified after they were abandoned exhibit self-affine geometry, preserving the fault roughness from different stages of faulting. Roughness varies little between surfaces, except the last slip zone to form in the fault, which is the smoothest. This layer contains a distinct mineral assemblage, indicating that the composition of the fault rock exerts a control on roughness. In contrast, the similar roughness of the older slip zones, which have comparable mineralogy but clearly crosscut one another, suggests that as the fault matured the roughness of the active slip surface stayed approximately constant. Wear processes affected these layers, so for roughness to stay constant the roughening and smoothing effects of fault slip must have been approximately balanced. These observations suggest fault surface evolution occurs by nucleation of new surfaces and wear by competing smoothing and re-roughening processes.     

冰盖下矩形河道的综合糙率

综合糙率是采用曼宁公式确定河道水位和流量关系的关键参数。在河道冰封期,冰盖的出现增加了流动的阻力,明流条件下确定的综合糙率不再适用,需要重新估算。基于Einstein阻力划分过流断面的原理,冰盖下矩形河道的过水断面可划分为冰盖区、河床区和边壁区。根据总流的连续性方程,在确定各分区糙率系数、水力半径和断面面积的基础上,提出了冰盖下矩形河道综合糙率的计算公式。采用已有的试验水槽测量数据和天然河道实测资料对公式进行了验证,结果表明:公式计算的综合糙率与实测值吻合较好,与Einstein公式和Sabaneev公式相比,计算精度更高;对于冰封水流,宽浅河道采用分区水深代替水力半径进行简化计算的条件有别于明渠水流,在宽深比大于20时,计算结果才满足精度要求。     

汶川地震地表破裂面形貌特征

准确描述破裂面形貌对于我们理解地震断层作用是非常重要的,破裂面的形貌特征包含许多关于地震和断层机制的有用信息。在Mw7.9 2008汶川地震中断层活动产生了两个新鲜的破裂面,八角庙破裂面和沙坝破裂面。我们使用3D便携式激光扫描仪(Tri mble GX)对两个破裂表面进行测量,在野外微观尺度上研究了破裂面形貌特征。通过能谱密度和均方值两个方法分析破裂面形貌,新鲜的破裂面表现为自相仿性,能谱密度和均方值均与剖面长度存在幂律关系。在能谱密度与空间频率的对数图中,能谱密度曲线存在明显的拐点,该拐点所对应的波长称为"特征波长",表明单一分形不能准确描述破裂面形貌。八角庙破裂面在平行滑动方向上的特征波长为7 mm,在垂直方向上特征波长略大一些(区域Ⅰ为10 mm,区域Ⅱ为9 mm);沙坝破裂面在平行滑动方向上的特征波长为8 mm,但垂直方向上特征波长略小(6 mm)。均方值曲线的最小二乘拟合直线的斜率为Hurst指数,该指数依赖于剖面线方向并描述破裂面形貌的各向异性,H指数的最小值和最大值分别与平行擦痕和垂直擦痕方向对应,这与野外断层面擦痕测量结果一致。沙坝破裂面的H指数极坐标图中存在次级H指数峰值(对应剖面线方向为85°和160°),这揭示破裂面上存在一组隐匿擦痕。该组隐匿擦痕为汶川地震之前断层活动中形成的,但这还不足以推测上一次断层活动的时间和规模。另外,通过比较新鲜节理面和破裂面表明H指数是否大于0.8反映了断层类型。在整个空间频率域上,能谱密度曲线斜率(-α)和均方值曲线斜率(H)的线性拟合关系为α=1.22+1.72H,并不严格满足两者之间理论关系式,α=1+2H。这个差异是由于测量信号噪音、破裂面的多分形性和分析方法的差异造成的。