基于ZY3 DEM和SRTM DEM的北京平原区地貌特征探析

研究区主要包括北京东部的通州、平谷、顺义、燕郊、大厂和三河,位于华北沉降带北部、阴山东西向构造带南缘。该区地势总体上西北略高,东南略低,为大型缓倾斜平原,地形不仅受本区内北东向断裂的影响,还受北西向断裂的影响。资源三号卫星立体影像(ZY3)DEM和NASA航天飞机雷达地形任务(SRTM)DEM在北京地区地貌特征研究中显示出了许多优于其他DEM的优势。文章利用包括水系偏转、河长坡降指标(SL)、面积-高程积分值(HI)、等值线分级色彩显示和地形粗糙度等几种地形参数,分析了地形特征与断层构造活动性的关系。水系的右旋偏转证实了该区北东向断裂的右旋走滑活动,如夏垫断裂和其他相邻断裂。SL指数表明,河床纵剖面在横跨断层处均是陡坡;HI的高值主要集中在南口-孙河断裂与该区其他北东向断裂的交汇处,夏垫断裂的中段,以及香河县西北部,反映了该区强烈的构造活动,同时也反映出可能存在隐伏断层。另外,研究表明,面积、高程积分值的高值区与低值区相邻衔接处,往往是地震危险性较大的地方。此外,平原区相同值域的面积-高程积分所反映的构造活动强度应该与山区不同,因此地貌演化阶段划分标准也应与山区有所不同。等值线分级色彩显示对地形变化非常敏感,当与已识别的断层重叠时,大部分地形突变的地方都与断层活动有关。地形粗糙度分析显示不同断裂带的地形粗糙度与断裂活动强度呈正相关关系。     

Optimization approach to retrieve soil surface parameters from single-acquisition single-configuration SAR data

This study suggests a novel approach to the retrieval of soil surface parameters using a single-acquisition single-configuration synthetic-aperture radar (SAR) system. Soil surface parameters such as soil moisture and surface roughness are key elements for many environmental studies, including Earth surface water cycles, energy exchange, agriculture, and geology. Remote sensing techniques, especially SAR data, are commonly used to retrieve such soil surface parameters over large areas. Several backscattering models have been proposed for soil surface parameters retrieval from SAR data. However, commonly, these backscattering models require multi configuration SAR data, including multi-polarization, multi-frequency, and multi-incidence angle. Here we propose a methodology that employs single-acquisition single-configuration SAR data for the retrieval of soil surface parameters. The originality is to use single-acquisition single-configuration SAR data to retrieve the soil surface parameters using an optimization approach by the genetic algorithm (GA); we have used the modified Dubois model (MDM) in HH polarization as the backscattering model. Three HH polarization and C band data sets from Quebec (Radarsat-1), Ontario (SIR-C), and Oklahoma (AIRSAR) were analyzed. The retrieved values of soil moisture and soil surface roughness were then compared to ground truth measurements with corresponding parameters. We employed diverse criteria, including the mean absolute error (MAE), the root mean square error (RMSE), the coefficient of performance (CP), and the correlation coefficient to investigate the performance of the proposed methodology. This analysis suggests the capability of the GA for the retrieval of soil surface parameters. Based on our findings, this method presents a viable alternative approach to the retrieval of soil surface parameters when only single-acquisition single-configuration SAR data is available.     

流动型态对曼宁糙率系数的影响研究

根据明渠水流流动型态的概念,通过对实验资料的分析,得到均匀流时糙率系数随水深增加而减小,形成M1型水面线的非均匀流时糙率系数随水深、水力坡度的增大而增大的规律。提出的糙率系数的二步计算法,可以解决因非均匀流水深沿程变化,难以建立糙率系数与水深关系的困难,为在水力计算中修正糙率系数提供了依据。     

Sources of Uncertainty in Shear Stress and Roughness Length Estimates Derived from Velocity Profiles*

Shear stress is a fundamental parameter in many sediment transport expressions. It is commonly estimated from information contained in the velocity profile, specifically, shear velocity, u_*, and roughness length, z₀ Under ideal conditions, the functional relationship between elevation above the bed and flow velocity is log-linear, as expressed by the “law of the wall.” Unless the field data conform exactly to this ideal relationship, there is uncertainty in estimates of u_* and z₀. derived from velocity profiles. This uncertainty depends on the degree of scatter or deviation from the assumed log-linear model. Expressions are presented to quantify the error and to correct for it. Guidelines are suggested to minimize potential uncertainty, especially with regard to instrument deployment and velocity profile analysis. Calculation of confidence intervals around estimates of u_* and z₀, as proposed by Wilkinson (1984), is necessary and recommended.     

Hydraulic conditions for rill incision under simulated rainfall: A laboratory experiment

A series of controlled laboratory experiments were conducted in order to obtain precise data on the hydraulic and sediment transport conditions during rill formation. Tests were carried out using a crusting-prone binary mixed soil in a 15 m long flume at an average slope of 0·087 under simulated rainfall. Rainfall intensities varied from 30–35 mm h^?1 and developed about 70 per cent of the kinetic energy of natural rainfall of similar intensity. Runoff and sediment discharge measured at the downstream weir were strongly influenced by rill forming processes. Essentially, rill incision reduced runoff discharge as a result of increased percolation in rill channels but increased sediment discharge. Secondary entrainment processes, such as bank collapse, also increased sediment discharge at the weir. Knickpoint bifurcation and colluvial deposition, however, decreased sediment discharge. Rills always developed through the formation of a knickpoint. The critical condition for knickpoint initiation was the development of supercritical flow and waves which mould and incise the bed. Prior smoothing of the soil surface by entrainment and redistribution of sediment facilitated supercritical flow. Statistical analysis showed that hydraulic and sediment transport conditions differed significantly in rilled and unrilled flows. The relationship between sediment discharge, rill erosion, and flow hydraulics was found to be nonlinear, conforming to a standard power function in the form y = ax^b. Rills were also associated with significantly increased sediment transport capacities. However, rill initiation was not clearly defined by any specific hydraulic threshold. Instead, rilled and unrilled flows were separated by zones of transition within which both types of flow occur.     

Observed Wind Profiles and Turbulence Fluxes over an ice Surface with Changing Surface Roughness

Wind profile and eddy-correlation data obtained at two sites on a melting glacier surface in Iceland during the summer of 1996 are presented. Throughout the experiment the surface roughness increased rapidly from smooth to very rough, with the largest roughness element height obtained being about 1.7 m. In a layer close to the rough surface we find that the wind speed profiles were disturbed showing horizontal inhomogeneities as in a roughness sublayer. Its height was approximately two times the height of the main roughness elements (h) at both sites throughout the experiment. From the wind profiles and eddy-correlation data we calculated corrections for the displaced zero plane as a function of time and compared these with results obtained from a drag partitioning model. In general, the agreement was reasonable considering the ranges of uncertainty but the results indicate that the increasing horizontal anisotropy of the surface probably limits the use of the model. The values obtained for the roughness lengths are in good agreement with those calculated from a simple linear model, i.e., z₀/h = 0.5 with the frontal area index. Above the roughness sublayer the wind profiles, normalised standard deviations of wind speed, and the balance of the turbulence kinetic energy budget behaved as over an ideal homogeneous surface thereby confirming similarity of the flow.     

Surface roughness signatures of summer arctic snow-covered sea ice in X-band dual-polarimetric SAR

ABSTRACT

Surface roughness of sea ice is primary information for understanding sea ice dynamics and air–ice–ocean interactions. Synthetic aperture radar (SAR) is a powerful tool for investigating sea ice surface roughness owing to the high sensitivity of its signal to surface structures. In this study, we explored the surface roughness signatures of the summer Arctic snow-covered first-year sea ice in X-band dual-polarimetric SAR in terms of the root mean square (RMS) height. Two ice campaigns were conducted for the first-year sea ice with dry snow cover in the marginal ice zone of the Chukchi Sea in August 2017 and August 2018, from which high-resolution (4 cm) digital surface models (DSMs) of the sea ice were derived with the help of a terrestrial laser scanner to obtain the in situ RMS height. X-band dual-polarimetric (HH and VV) SAR data (3 m spatial resolution) were obtained for the 2017 campaign, at a high incidence angle (49.5°) of TerraSAR-X, and for the 2018 campaign, at a mid-incidence angle (36.1°) of TanDEM-X 1–2 days after the acquisition of the DSMs. The sea ice drifted during the time between the SAR and DSM acquisitions. As it is difficult to directly co-register the DSM to SAR owing to the difference in spatial resolution, the two datasets were geometrically matched using unmanned aerial vehicle (4 cm resolution) and helicopter-borne (30 cm resolution) photographs acquired as part of the ice campaigns. A total of five dual-polarimetric SAR features―backscattering coefficients at HH and VV polarizations, co-polarization ratio, co-polarization phase difference, and co-polarization correlation coefficient ―were computed from the dual-polarimetric SAR data and compared to the RMS height of the sea ice, which showed macroscale surface roughness. All the SAR features obtained at the high incidence angle were statistically weakly correlated with the RMS height of the sea ice, possibly influenced by the low backscattering close to the noise level that is attributed to the high incidence angle. The SAR features at the mid-incidence angle showed a statistically significant correlation with the RMS height of the sea ice, with Spearman’s correlation coefficient being higher than 0.7, except for the co-polarization ratio. Among the intensity-based and polarimetry-based SAR features, HH-polarized backscattering and co-polarization phase difference were analyzed to be the most sensitive to the macroscale RMS height of the sea ice. Our results show that the X-band dual-polarimetric SAR at mid-incidence angle exhibits potential for estimation of the macroscale surface roughness of the first-year sea ice with dry snow cover in summer.     

Horizontal Scale of Roughness Variations for Realistic Landscapes

A method is described to infer the horizontal scale of roughness variations from land-use maps, for use in flux aggregation schemes designed to account for sub-grid scale terrain heterogeneity in numerical models of the atmosphere. Based on simple statistical-geometrical considerations, it is shown that this horizontal scale is inversely proportional to edge density, which can be readily derived from high-resolution land-use maps. The method is demonstrated for a highly heterogeneous landscape in northern Belgium, the resulting length scale values approximately ranging between 200 and 1500 m. Finally, the extension towards satellite imagery and coarse land-use data is discussed.     

Effect of depression storage capacity on overland‐flow generation for rough horizontal surfaces: water transfer distance and scaling

Overland‐flow triggering on rough surfaces was investigated using an understanding‐oriented model. The model was based on conditioned‐walker technique and developed to simulate and analyse the evolution of puddle connection on numerically generated rough surfaces. The percolation theory gave a theoretical framework to formalize model outputs and to study overland‐flow scaling. Overland‐flow triggering appeared consistent with a percolation process. A scale‐change exponent was suggested. New insights based on the concept of transfer distance of water were emphasized. Transfer distance enabled us to analyse the water redistribution inside a field and helped to define rainfall efficiency when infiltration occurred. Copyright © 2002 John Wiley & Sons, Ltd.