桂林-阳朔地区DEM地形特征与岩性相关性分析及分类研究

地形地貌是岩性解译的重要信息,地形因子作为描述DEM数字曲面几何特征的定量指标参数,可用来定量化表达不同岩性所在地区地形地貌特征。本文以桂林-阳朔地区为研究区,研究地形因子数学、地质意义,建立岩性与地形因子组合间的定量关联,进而实现岩石类型划分。本文基于ASTERGDEM提取坡度、起伏度等12个地形因子,在分析各个地形因子地质意义基础上,通过聚类分析及方差分析的多元统计分析方法,研究各岩性地形因子特性及其关联性,建立研究区岩性之间的定量差异;此外,利用因子分析方法研究岩性分类过程中的主导因素,确定适宜岩性分类方法以实现定量化岩性分类。实验结果表明：不同岩性、不同地形地貌的地形因子（组合）之间具有显著差异,基于因子分析得到的宏观地形复杂度指数（MTI）以及微观曲率指数（MCI）对岩石类型的分类精度达77.36%。研究表明,地形复杂度等地形因子可用于岩性分类,采用因子分析方法可获取反映地形地貌宏观、微观特征的定量指标,且岩性分类效果良好。     

Revisiting the existence of an effective stress for wet granular soils with micromechanics

A possible effective stress variable for wet granular materials is numerically investigated based on an adapted discrete element method (DEM) model for an ideal three‐phase system. The DEM simulations consider granular materials made of nearly monodisperse spherical particles, in the pendular regime with the pore fluid mixture consisting of distinct water menisci bridging particle pairs. The contact force‐related stress contribution to the total stresses is isolated and tested as the effective stress candidate for dense or loose systems. It is first recalled that this contact stress tensor is indeed an adequate effective stress that describes stress limit states of wet samples with the same Mohr‐Coulomb criterion associated with their dry counterparts. As for constitutive relationships, it is demonstrated that the contact stress tensor used in conjunction with dry constitutive relations does describe the strains of wet samples during an initial strain regime but not beyond. Outside this so‐called quasi‐static strain regime, whose extent is much greater for dense than loose materials, dramatic changes in the contact network prevent macroscale contact stress‐strain relationships to apply in the same manner to dry and unsaturated conditions. The presented numerical results also reveal unexpected constitutive bifurcations for the loose material, related to stick‐slip macrobehavior.     

Determination of subcatchment and watershed boundaries in a complex and highly urbanized landscape

Urban development significantly alters the landscape by introducing widespread impervious surfaces, which quickly convey surface run‐off to streams via stormwater sewer networks, resulting in “flashy” hydrological responses. Here, we present the inadequacies of using raster‐based digital elevation models and flow‐direction algorithms to delineate large and highly urbanized watersheds and propose an alternative approach that accounts for the influence of anthropogenically modified land cover. We use a semi‐automated approach that incorporates conventional drainage networks into overland flow paths and define the maximal run‐off contributing area. In this approach, stormwater pipes are clustered according to their slope attributes, which define flow direction. Land areas drained by each cluster and contributing (or exporting) flow to a topographically delineated catchment were determined. These land masses were subsequently added or removed from the catchment, modifying both the shape and the size. Our results in a highly urbanized Toronto, Canada, area watershed indicate a moderate net increase in the directly connected watershed area by 3% relative to a topographically forced method; however, differences across three smaller scale subcatchments are greater. Compared to topographic delineation, the directly connected watershed areas of both the upper and middle subcatchments decrease by 5% and 8%, respectively, whereas the lower subcatchment area increases by 15%. This is directly related to subsurface storm sewer pipes that cross topographic boundaries. When directly connected subcatchment area is plotted against total streamflow and flashiness indices using this method, the coefficients of variation are greater (0.93 to 0.97) compared to the use of digital elevation model‐derived subcatchment areas (0.78 to 0.85). The accurate identification of watershed and subcatchment boundaries should incorporate ancillary data such as stormwater sewer networks and retention basin drainage areas to reduce water budget errors in urban systems.     

Quantifying terrain controls on runoff retention and routing in the Northern Prairies

The role of hummocky terrain in governing runoff routing and focussing groundwater recharge in the Northern Prairies of North America is widely recognised. However, most hydrological studies in the region have not effectively utilised information on the surficial geology and associated landforms in large-scale hydrological characterization. The present study uses an automated digital elevation model (DEM) analysis of a 6500-km² area in the Northern Prairies to quantify hydrologically relevant terrain parameters for the common types of terrains in the prairies with different surficial deposits widespread in the prairies, namely, moraines and glaciolacustrine deposits. Runoff retention (and storage) capacity within depressions varies greatly between different surficial deposits and is comparable in magnitude with a typical amount of seasonal snowmelt runoff generation. The terrain constraint on potential runoff retention varies from a few millimetres in areas classified as moraine to tens of millimetres in areas classified as stagnant ice moraine deposits. Fluted moraine and glaciolacustrine deposits have intermediate storage capacity values. The study also identified the probability density function describing a number of immediate upstream neighbours for each depression in a fill-and-spill network. A relationship between depression parameters and surficial deposits, as well as identified depression network structure, allows parametrisation of hydrologic models outside of the high-resolution DEM coverage, which can still account for terrain variation in the Prairies.     

山地叶面积指数反演理论、方法与研究进展

叶面积指数LAI（Leaf Area Index）是表征叶片疏密程度和冠层结构特征的重要植被参数，在气候变化、作物生长模型以及碳、水循环研究中发挥着重要作用。遥感是获取区域及全球尺度LAI的一个重要手段，当前LAI产品主要基于遥感数据反演得到，但是多数LAI产品算法并未考虑地形特征的影响，导致山地LAI遥感反演精度不确定性大。提高山地LAI遥感反演精度亟需考虑地形因子对冠层反射率的影响，其中山地冠层反射率模型和遥感数据地形校正是提升山地LAI遥感反演精度的关键。本文围绕山地LAI遥感反演理论与方法，综合分析了国内外山地冠层反射率模型和地形校正模型的研究进展，总结了目前山地LAI遥感反演存在的问题，并讨论了未来研究的发展趋势。     

基于2 m格网点云数据的高精度DEM制作第三次国土调查坡度图

地形地貌随着时间的推移时刻发生变化，以高精度数字高程模型（EDE）数据制作坡度图，进而计算田坎系数，更精准地进行耕地面积计算及统计。以2 m格网点云数据生成的DEM数据为基础，按照第三次全国国土调查的相关技术要求，根据耕地坡度分级要求进行分级，生成坡度分级栅格数据图。对坡度分级栅格数据图进行矢量化，生成坡度分级矢量化数据。对矢量化数据进行图斑综合、界线平滑、拓扑重建、数据裁切等处理，制作完成调查区域坡度图及相关属性数据制作。并对生产的坡度图成果进行分析，找出不同尺度格网生产的坡度图的技术差异，并对高精度坡度的应用进行了展望。     

基于“DEM-NDVI-土地覆盖分类”的天山博格达自然遗产地山地垂直带提取与变化分析

本文基于Landsat影像数据获取天山博格达自然遗产地土地覆盖分类,结合归一化植被指数（NDVI）和数字高程模型（DEM）构建“DEM-NDVI-土地覆盖分类”散点图分析研究区植被受海拔和坡向的水热空间变化影响的分布特征,通过概率统计分析提取博格达遗产地山地垂直带,并结合研究区的气温、降水数据和NDVI变化特征分析垂直带变化的原因。研究结果表明：① 本文利用“DEM-NDVI-土地覆盖分类”散点图,揭示了研究区1989年和2016年的NDVI值和分类类别随着海拔上升的变化特征,其中NDVI值随着海拔上升呈现“倒U形”变化,而不同分类类别在一定的海拔区间内呈现出聚集效应,且不同分类类别有明显的高程界限。② 1989年和2016年博格达遗产地山地垂直带分带上限分别为：1278 m和1185 m（温带荒漠草原带）、1784 m和1759 m（山地草原带）、2706 m和2730 m（山地针叶林带）、3272 m和3293 m（高山草甸带）、3636 m和3690 m（高山垫状植被带）。③ 博格达遗产地1989年和2016年山地垂直带受区域气温升高和降雨增加的影响有较为明显的改变,其中温带荒漠草原带最为敏感,其上限变化最大,向下收缩93 m;山地针叶林带的分布范围则向两侧扩张49 m;山地草甸带带宽基本保持不变,但整体上移了约20 m;冰雪带则受到全球气候变暖的影响向上退缩54 m。     

Extracting drainage networks and their connectivity using LiDAR data

Policies, measures, and models geared towards flood prevention and managing surface waters benefit from high quality data on the presence and characteristics of drainage ditches. As a cost and labour effective alternative for acquiring such data through field surveys, we propose a method (a) to extract vector data representing ditch drainage networks based on local morphologic features derived from high resolution digital elevation models (DEM) and (b) to identify possible connections in the ditch network by calculating a probability of the connectivity using a logistic regression where the predictor variables are characteristics of the ditch centre lines or derived from the DEM. Using Light Detection and Ranging (LiDAR) derived DEMs with a 1 m resolution, the method was developed and tested for a mixed agricultural residential area in north‐eastern Belgium. The derived ditch segments had an error of omission of 8% and an error of commission of 5%. The original positional accuracy of the centre lines of the extracted ditches was 0.6 m and could be improved to 0.4 m by shifting each vertex to the position of the lowest LiDAR point located within a radius equal to the spatial resolution of the used DEM. About 69% of the false disconnections in the network were identified and corrected leading to a reduction of the unconnected parts of the ditch network by 71%. The extracted and connected network approximated the reference ditch network fairly well.     

THE APPLICATION OF MINIATURE UNMANNED AERIAL VEHICLE IN 25 NOVEMBER 2016 ARKETAO MW6.6 EARTHQUAKE

In order to complete the field investigation to the 25 November 2016 Arketao M_W6.6 earthquake, ultra-low altitude remote-sensing data were obtained from miniature unmanned aerial vehicle. The surface rupture surveying has important significance for earthquake research. This paper selects the macro-epicenter of Arketao as the study area. The pictures were obtained with DJI Phantom 3 professional input into the software, the Digital Elevation Model (DEM), Digital Orthophoto Map (DOM) were acquired based on photogrammetry method using the overlapped optical remote-sensing images of UAV. Using these data, we can identify surface ruptures that have vertical dislocation. We selected six feature points and drew the elevation profile. In the elevation profile map, we chose smooth part of the surface rupture sides and obtained the trend line. A stable point in the surface rupture was selected and the abscissa of the point was taken into the equation of two straight lines. Then subtracting the results of the two equations, we can get the vertical dislocation of the surface rupture. On this basis, we chose six feature points and determined their vertical dislocation, which are between 4.4cm and 10.4cm. What's more, taking Bulungkou Xiang in Xinjiang Uygur Autonomous Region for example, we speculated some surface ruptures that have vertical dislocation. It can provide a new method for identifying surface rupture in the field. In addition, we get DEM data of the Bulunkou area where ambient conditions are very poor, by using miniature unmanned aerial vehicle and taking 255 photos. Putting those photos into the EasyUAV software, we got the area digital elevation of 2cm resolution. Comparing these data with RTK data, we summarized some practical problems and solutions in the practical operation and evaluated the accuracy of miniature unmanned aerial vehicle data. The Pearson Correlation Coefficient is 0.996 6. In terms of absolute elevation, the average result of UAV and RTK differs by 156.96m. In terms of relative elevation, the average result of UAV and RTK differs by 9.74m. Compared with the previous test of Pishan County, there is a notable divergence in the results. It shows that the data accuracy will be affected to some extent in the cold weather in high elevations. The specific impact needs further exploration.     

基于LiDAR-DP的机载激光雷达数据处理

基于国产机载LiDAR数据处理软件LiDAR-DP,总结了机载LiDAR数据处理的技术流程,并根据该流程对LiDARDP的相关功能进行了研究。利用LiDAR-DP对实际数据进行了DEM生产试验,并与同类商业软件进行了功能对比。结果表明,LiDAR-DP是一款优秀的国产机载LiDAR数据处理软件,完全能满足DEM的生产需求。