格网DEM上径流长度计算误差的定量分析

基于格网DEM的径流长度的计算误差主要来源于DEM数据和流径提取算法。通过数据独立的误差分析方法,定量分析流径提取算法模型产生的误差。首先选取5种数学曲面作为模拟DEM数据源,计算出径流长度,以由曲面公式推导的径流长度公式计算的径流长度作为理论真值,通过两者的对比,分析误差及其空间分布。结果表明,汇水区的径流长度计算精度明显高于分水区;分水区的径流长度含有一定的系统误差。     

粗糙元几何参数的交互作用对床面空气动力学粗糙度的影响

 为了全面地揭示粗糙元的所有几何参数的交互作用对空气动力学粗糙度的影响,利用风洞实验研究了粗糙元高度、密度、高度与间距比、孔隙度和方向比率等几何参数交互作用对空气动力学粗糙度的影响。结果表明,密实和孔隙粗糙元的无量纲空气动力学粗糙度（空气动力学粗糙度度/粗糙元高度）均可表示为粗糙元密度/等效密度的正比例函数,而比例系数反映了粗糙元几何参数交互作用。据此,该研究发展了一个全面反映粗糙元高度h、密度/等效密度λ、高度与行间距比Sp和方向比率AR等几何参数交互作用的空气动力学粗糙度模式:Z0h=-0.0028+0.5403S0.32p·AR-0.07·λ。该模式改进了模拟的精度,扩大了适用范围。     

The Effect of Stratification on the Aerodynamic Roughness Length and Displacement Height

The roughness length, z _0u , and displacement height, d _0u , characterise the resistance exerted by the roughness elements on turbulent flows and provide a conventional boundary condition for a wide range of turbulent-flow problems. Classical laboratory experiments and theories treat z _0u and d _0u as geometric parameters independent of the characteristics of the flow. In this paper, we demonstrate essential stability dependences—stronger for the roughness length (especially in stable stratification) and weaker but still pronounced for the displacement height. We develop a scaling-analysis model for these dependences and verify it against experimental data.     

Temporal and Spatial Variations of the Aerodynamic Roughness Length in the Ablation Zone of the Greenland Ice Sheet

To understand the response of the Greenland ice sheet to climate change the so-called ablation zone is of particular importance, since it accommodates the yearly net surface ice loss. In numerical models and for data analysis, the bulk aerodynamic method is often used to calculate the turbulent surface fluxes, for which the aerodynamic roughness length (z ₀) is a key parameter. We present, for the first time, spatial and temporal variations of z ₀ in the ablation area of the Greenland ice sheet using year-round data from three automatic weather stations and one eddy-correlation mast. The temporal variation of z ₀ is found to be very high in the lower ablation area (factor 500) with, at the end of the summer melt, a maximum in spatial variation for the whole ablation area of a factor 1000. The variation in time matches the onset of the accumulation and ablation season as recovered by sonic height rangers. During winter, snow accumulation and redistribution by snow drift lead to a uniform value of z ₀≈ 10⁻⁴ m throughout the ablation area. At the beginning of summer, snow melt uncovers ice hummocks and z ₀ quickly increases well above 10⁻² m in the lower ablation area. At the end of summer melt, hummocky ice dominates the surface with z ₀ > 5  ×  10⁻³ m up to 60 km from the ice edge. At the same time, the area close to the equilibrium line (about 90 km from the ice edge) remains very smooth with z ₀ = 10⁻⁵ m. At the beginning of winter, we observed that single snow events have the potential to lower z ₀ for a very rough ice surface by a factor of 20 to 50. The total surface drag of the abundant small-scale ice hummocks apparently dominates over the less frequent large domes and deep gullies. The latter results are verified by studying the individual drag contributions of hummocks and domes with a drag partition model.     

The Parameterisation of Scalar Transfer over Rough Ice

We test a surface renewal model that is widely used over snow and ice surfaces to calculate the scalar roughness length (z _s ), one of the key parameters in the bulk aerodynamic method. For the first time, the model is tested against observations that cover a wide range of aerodynamic roughness lengths (z ₀). During the experiments, performed in the ablation areas of the Greenland ice sheet and the Vatnajökull ice cap in Iceland, the surface varied from smooth snow to very rough hummocky ice. Over relatively smooth snow and ice with z ₀ below a threshold value of approximately 10^?3 m, the model performs well and in accord with earlier studies. However, with growing hummock size, z ₀ increases well above the threshold and the bulk aerodynamic flux becomes significantly smaller than the eddy-correlation flux (e.g. for z ₀ = 0.01 m, the bulk aerodynamic flux is about 50% smaller). Apparently, the model severely underpredicts z _s over hummocky ice. We argue that the surface renewal model does not account for the deep inhomogeneous roughness sublayer (RSL) that is generated by the hummocks. As a consequence, the homogeneous substrate ice grain cover becomes more efficiently ‘ventilated’. Calculations with an alternative model that includes the RSL and was adapted for use over hummocky ice, qualitatively confirms our observations. We suggest that, whenever exceedance of the threshold occurs (z ₀  >  10^?3 m, i.e., an ice surface covered with at least 0.3-m high hummocks), the following relation should be used to calculate scalar roughness lengths, ln (z _s /z ₀)  =  1.5  ? 0.2 ln (Re _*)  ? 0.11(ln (Re _*))².     

A straightforward approach for using single time domain simulations to assess characteristic extreme responses

Short-term wave design approach of marine structures, using nonlinear time domain simulations, is a design procedure that is recognized by various modern standard codes. One of the most challenging points of this approach is the evaluation of the characteristic extreme values for response parameters used in the design check equations. The most straightforward and recommended way to evaluate a response characteristic value is by fitting an extreme value probability distribution to the N-sample of extreme values extracted from N independent time domain simulations with duration equal to the short-term period indicated by the code, which is usually taken as 3 h. However, this procedure would not be practical for some types of marine structures, such as risers and mooring lines, under numerous design load cases and demanding huge finite element models. A more feasible approach would be to assess the response extreme value distribution using only a single short-term time domain simulation with duration shorter than 3 h. But reduced time simulations always introduce some additional statistical uncertainty into the extreme values estimates. This paper discusses a workable way of properly taking into account the statistical uncertainty associated with the simulation length in the assessment of a characteristic short-term extreme response value based on a single time series.     

独立坐标系模型对长度变形影响的研究

针对海拔较高的城市建立或改造独立坐标系,存在着如何选择模型及各模型对长度变形的影响等问题,通过算例对常见椭球变换和比例缩放模型进行较为全面系统的分析,归纳出常见模型在不同投影带区域、不同投影面高度对长度变形的影响差异,并得出有益的结论。     

杭州湾淡化对虾养殖池水中氮磷营养盐的存在特征

利用2002年5～7月对杭州湾淡化对虾养殖虾池水环境中氮磷营养盐的调查资料,分析了该养殖水环境中氮磷营养盐的存在形态与行为,结果表明:(1)在该养殖水环境中,IP含量随时间变化不大,相对较稳定;(2)IN含量的变化主要受生物活动控制,在养殖过程中变化较大,NO3-N为IN的主要存在形态,但在养殖后期,NO2-N含量有所增加;(3)该养殖水环境中氮磷营养盐含量均超过富营养化阈值,因而已无法判断浮游植物生长的限制因子;(4)叶绿素a与IN含量成显著负相关,而与IP含量的相关性不显著。     

城市下垫面空气动力学参数的确定

分析近10年北京325m气象塔常年观测资料,研究城市化发展对城市大气边界层动力学结构和特征的影响.结果表明,空气动力学粗糙度和零平面位移总体是逐年递增的.气象塔西南方向的高层建筑群是20世纪90年代初期逐渐建成,在SW方向,空气动力学粗糙度和零平面位移在90年代增加明显,增加的幅度1997年和1999年之间较大,这一点与气象塔西南方向城市化的加快相符合.而在NE、SE和NW方向,从1987年到1994年空气动力学粗糙度和零平面位移增加明显,1994年以后变化较小.     

Influence of length effect on embankment slope reliability in 3D

Embankment slopes composed of spatially variable soils have a variety of different failure modes that are affected by the correlation distances of the material properties and the geometry and total length of the slope. This paper examines the reliability of soil slopes for embankments of different length and uses parallel computing to analyse very long embankments (up to 100 times the embankment height) for a clay soil characterised by a spatially varying undrained shear strength. Based on a series of analyses using the 3D random finite element method (RFEM), it is first shown that the reliability of slopes of various length can be efficiently computed by combining simple probability theory with a detailed 3D RFEM analysis of a representative shorter slope of length 10 times the slope height. RFEM predictions of reliability indices for longer slopes are then compared with results obtained using Vanmarcke's (1977a) simplified 3D method and Calle's (1985) extended 2D approach. It is shown that these methods can give significantly different results, depending on the horizontal scale of fluctuation relative to the slope length, with RFEM predicting a lower slope reliability than the Vanmarcke and Calle solutions in all cases. The differences in the solutions are evaluated and attributed to differences in the assumed and computed failure surface geometries.