共查询到20条相似文献，搜索用时 15 毫秒
1.
Coarse clastic sediments (boulders) on coastlines have seen a groundswell in geomorphic research interest over recent years, associated in part with the potential of boulder evidence for interpreting characteristics of high‐energy wave processes. Yet, the fundamental property of boulder volume is normally difficult to measure accurately owing to complex clast morphology and irregular surface texture. To tackle this problem, this paper concentrates on creating precise, measurable and textured three‐dimensional (3D) models of coastal boulders without physical contact with the object, based on multi‐view image measurement techniques. This method has several advantages over traditional measurements that are inaccurate or alternative solutions using costly techniques such as terrestrial laser scanning. Our methods propose the use of low‐cost equipment (digital cameras) that can be used in various coastal environments to easily acquire numerous images of the object of interest. Initial results can be rapidly assessed in the field for immediate quality control. Resulting 3D models, built from overlapping multi‐view digital photographs, allow the reconstruction of realistic‐looking and textured boulder surfaces. A particular interest in this task is the family of algorithms known as structure from motion (SFM). The work presents analysis of SFM techniques by examining 3D models of boulders observed at a coastal field site on Lu Dao Island in south‐eastern Taiwan. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
2.
Testing the use of an image‐based technique to measure gully erosion at Sparacia experimental area 
下载免费PDF全文
下载免费PDF全文 Costanza Di Stefano Vito Ferro Vincenzo Palmeri Vincenzo Pampalone Fabrizio Agnello 《水文研究》2017,31(3):573-585
The first part of this investigation was aimed at testing the use of a three‐dimensional (3D) digital terrain model and a quasi‐tridimensional (2.5D) digital elevation model obtained by a large series of oblique images of eroded channels taken from consumer un‐calibrated and non‐metric cameras. For two closed earth channels having a different sinuosity, the ground measurement of some cross sections by a profilometer (P) was carried out and their real volume was also measured. The comparison among the three methods (3D, 2.5D, and P) pointed out that a limited underestimation of the total volume always occurs and that the 3D method is characterized by the minimum difference between measured and real volume. For this reason, 3D model can be used as benchmark. In the subsequent part of the investigation, the three ground measurement methods were applied for surveying of an ephemeral gully (EG) channel at the Sparacia area. The morphological and hydraulic variable values of the 24 surveyed cross sections determined by both 2.5D model and profilometer were compared. This comparison showed that the estimate error is generally less than ±10%. The EG measurements carried out by the three methods supported the applicability both of the empirical relationship between EG length and its eroded volume and the theoretical dimensionless relationship among the morphological variables describing the channelized erosion process. Finally, it was demonstrated that the effect of the distance interval on the EG volume measurement by 3D and 2.5D models is negligible for the investigated EG. 相似文献
3.
Saffet Erdogan 《地球表面变化过程与地形》2009,34(3):366-376
Digital elevation models have been used in many applications since they came into use in the late 1950s. It is an essential tool for applications that are concerned with the Earth's surface such as hydrology, geology, cartography, geomorphology, engineering applications, landscape architecture and so on. However, there are some differences in assessing the accuracy of digital elevation models for specific applications. Different applications require different levels of accuracy from digital elevation models. In this study, the magnitudes and spatial patterning of elevation errors were therefore examined, using different interpolation methods. Measurements were performed with theodolite and levelling. Previous research has demonstrated the effects of interpolation methods and the nature of errors in digital elevation models obtained with indirect survey methods for small‐scale areas. The purpose of this study was therefore to investigate the size and spatial patterning of errors in digital elevation models obtained with direct survey methods for large‐scale areas, comparing Inverse Distance Weighting, Radial Basis Functions and Kriging interpolation methods to generate digital elevation models. The study is important because it shows how the accuracy of the digital elevation model is related to data density and the interpolation algorithm used. Cross validation, split‐sample and jack‐knifing validation methods were used to evaluate the errors. Global and local spatial auto‐correlation indices were then used to examine the error clustering. Finally, slope and curvature parameters of the area were modelled depending on the error residuals using ordinary least regression analyses. In this case, the best results were obtained using the thin plate spline algorithm. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
4.
Manuel Nitsche Jens M. Turowski Alexandre Badoux Dieter Rickenmann Tobias K. Kohoutek Michael Pauli James W. Kirchner 《地球表面变化过程与地形》2013,38(8):810-825
Topographic measurements are essential for the study of earth surface processes. Three‐dimensional data have been conventionally obtained through terrestrial laser scanning or photogrammetric methods. However, particularly in steep and rough terrain, high‐resolution field measurements remain challenging and often require new creative approaches. In this paper, range imaging is evaluated as an alternative method for obtaining surface data in such complex environments. Range imaging is an emerging time‐of‐flight technology, using phase shift measurements on a multi‐pixel sensor to generate a distance image of a surface. Its suitability for field measurements has yet not been tested. We found ambient light and surface reflectivity to be the main factors affecting error in distance measurements. Low‐reflectivity surfaces and strong illumination contrasts under direct exposure to sunlight lead to noisy distance measurements. However, regardless of lighting conditions, the accuracy of range imaging was markedly improved by averaging multiple images of the same scene. For medium ambient lighting (shade) and a light‐coloured surface the measurement uncertainty was approximately 9 mm. To further test the suitability of range imaging for field applications we measured a reach of a steep mountain stream with a horizontal resolution of approximately 1 cm (in the focal plane of the camera), allowing for the interpolation of a digital elevation model on a 2 cm grid. Comparison with an elevation model obtained from terrestrial laser scanning for the same site revealed that both models show similar degrees of topographic detail. Despite limitations in measurement range and accuracy, particularly at bright ambient lighting, range imaging offers three‐dimensional data in real time and video mode without the need of post‐processing. Therefore, range imaging is a useful complement or alternative to existing methods for high‐resolution measurements in small‐ to medium‐scale field sites. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
5.
German Hoecht Patrice Ricarte Steffen Bergler Evgeny Landa 《Geophysical Prospecting》2009,57(6):957-979
In common‐reflection‐surface imaging the reflection arrival time field is parameterized by operators that are of higher dimension or order than in conventional methods. Using the common‐reflection‐surface approach locally in the unmigrated prestack data domain opens a potential for trace regularization and interpolation. In most data interpolation methods based on local coherency estimation, a single operator is designed for a target sample and the output amplitude is defined as a weighted average along the operator. This approach may fail in presence of interfering events or strong amplitude and phase variations. In this paper we introduce an alternative scheme in which there is no need for an operator to be defined at the target sample itself. Instead, the amplitude at a target sample is constructed from multiple operators estimated at different positions. In this case one operator may contribute to the construction of several target samples. Vice versa, a target sample might receive contributions from different operators. Operators are determined on a grid which can be sparser than the output grid. This allows to dramatically decrease the computational costs. In addition, the use of multiple operators for a single target sample stabilizes the interpolation results and implicitly allows several contributions in case of interfering events. Due to the considerable computational expense, common‐reflection‐surface interpolation is limited to work in subsets of the prestack data. We present the general workflow of a common‐reflection‐surface‐based regularization/interpolation for 3D data volumes. This workflow has been applied to an OBC common‐receiver volume and binned common‐offset subsets of a 3D marine data set. The impact of a common‐reflection‐surface regularization is demonstrated by means of a subsequent time migration. In comparison to the time migrations of the original and DMO‐interpolated data, the results show particular improvements in view of the continuity of reflections events. This gain is confirmed by an automatic picking of a horizon in the stacked time migrations. 相似文献
6.
7.
Distributed hydrologic models based on triangulated irregular networks (TIN) provide a means for computational efficiency in small to large‐scale watershed modelling through an adaptive, multiple resolution representation of complex basin topography. Despite previous research with TIN‐based hydrology models, the effect of triangulated terrain resolution on basin hydrologic response has received surprisingly little attention. Evaluating the impact of adaptive gridding on hydrologic response is important for determining the level of detail required in a terrain model. In this study, we address the spatial sensitivity of the TIN‐based Real‐time Integrated Basin Simulator (tRIBS) in order to assess the variability in the basin‐averaged and distributed hydrologic response (water balance, runoff mechanisms, surface saturation, groundwater dynamics) with respect to changes in topographic resolution. Prior to hydrologic simulations, we describe the generation of TIN models that effectively capture topographic and hydrographic variability from grid digital elevation models. In addition, we discuss the sampling methods and performance metrics utilized in the spatial aggregation of triangulated terrain models. For a 64 km2 catchment in northeastern Oklahoma, we conduct a multiple resolution validation experiment by utilizing the tRIBS model over a wide range of spatial aggregation levels. Hydrologic performance is assessed as a function of the terrain resolution, with the variability in basin response attributed to variations in the coupled surface–subsurface dynamics. In particular, resolving the near‐stream, variable source area is found to be a key determinant of model behaviour as it controls the dynamic saturation pattern and its effect on rainfall partitioning. A relationship between the hydrologic sensitivity to resolution and the spatial aggregation of terrain attributes is presented as an effective means for selecting the model resolution. Finally, the study highlights the important effects of terrain resolution on distributed hydrologic model response and provides insight into the multiple resolution calibration and validation of TIN‐based hydrology models. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
8.
David L. Hochstetler Warren Barrash Carsten Leven Michael Cardiff Francesco Chidichimo Peter K. Kitanidis 《Ground water》2016,54(2):171-185
Hydraulic tomography is an emerging field and modeling method that provides a continuous hydraulic conductivity (K) distribution for an investigated region. Characterization approaches that rely on interpolation between one‐dimensional (1D) profiles have limited ability to accurately identify high‐K channels, juxtapositions of lenses with high K contrast, and breaches in layers or channels between such profiles. However, locating these features is especially important for groundwater flow and transport modeling, and for design and operation of in situ remediation in complex hydrogeologic environments. We use transient hydraulic tomography to estimate 3D K in a volume of 15‐m diameter by 20‐m saturated thickness in a highly heterogeneous unconfined alluvial (clay to sand‐and‐gravel) aquifer with a K range of approximately seven orders of magnitude at an active industrial site in Assemini, Sardinia, Italy. A modified Levenberg‐Marquardt algorithm was used for geostatistical inversion to deal with the nonlinear nature of the highly heterogeneous system. The imaging results are validated with pumping tests not used in the tomographic inversion. These tests were conducted from three of five clusters of continuous multichannel tubing (CMTs) installed for observation in the tomographic testing. Locations of high‐K continuity and discontinuity, juxtaposition of very high‐K and very low‐K lenses, and low‐K “plugs” are evident in regions of the investigated volume where they likely would not have been identified with interpolation from 1D profiles at the positions of the pumping well and five CMT clusters. Quality assessment methods identified a suspect high‐K feature between the tested volume and a lateral boundary of the model. 相似文献
9.
《地震工程与结构动力学》2018,47(5):1212-1228
This paper investigates the 3D response of a slab supported by wobbling columns. The columns are not allowed neither to slide nor to roll out of their initial position. An analytical model is proposed, the equations of motion are derived, and they are solved numerically. The paper concludes that the addition of the slab makes the columns more stable. In fact, the system is almost equivalent to the response of a solitary column with the same aspect ratio yet larger size. However, it is also shown that the system is less stable than its planar counterpart and that planar analysis can only qualitatively describe the behavior of 3D structures. A case study shows that the concept could be used as a seismic isolation technique for bridges. However, more research need to be performed on defining proper intensity measures for uplifting structures, as it is shown that there is large record‐to‐record variability, even when intensity measures developed for rocking structures are used. 相似文献
10.
Near‐surface problem is a common challenge faced by land seismic data processing, where often, due to near‐surface anomalies, events of interest are obscured. One method to handle this challenge is near‐surface layer replacement, which is a wavefield reconstruction process based on downward wavefield extrapolation with the near‐surface velocity model and upward wavefield extrapolation with a replacement velocity model. This requires, in theory, that the original wavefield should be densely sampled. In reality, data acquisition is always sparse due to economic reasons, and as a result in the near‐surface layer replacement data interpolation should be resorted to. For datasets with near‐surface challenges, because of the complex event behaviour, a suitable interpolation scheme by itself is a challenging problem, and this, in turn, makes it difficult to carry out the near‐surface layer replacement. In this research note, we first point out that the final objective of the near‐surface layer replacement is not to obtain a newly reconstructed wavefield but to obtain a better final image. Next, based upon this finding, we propose a new thinking, interpolation‐free near‐surface layer replacement, which can handle complex datasets without any interpolation. Data volume expansion is the key idea, and with its help, the interpolation‐free near‐surface layer replacement is capable of preserving the valuable information of areas of interest in the original dataset. Two datasets, i.e., a two‐dimensional synthetic dataset and a three‐dimensional field dataset, are used to demonstrate this idea. One conclusion that can be drawn is that an attempt to interpolate data before layer replacement may deteriorate the final image after layer replacement, whereas interpolation‐free near‐surface layer replacement preserves all image details in the subsurface. 相似文献
11.
Antoine Guitton 《Geophysical Prospecting》2006,54(2):135-152
Surface‐related multiples are attenuated for one sail line and one streamer of a 3D data set (courtesy of Compagnie Générale de Géophysique). The survey was carried out in the Gulf of Mexico in the Green Canyon area where salt intrusions close to the water‐bottom are present. Because of the complexity of the subsurface, a wavefield method incorporating the full 3D volume of the data for multiple removal is necessary. This method comprises modelling of the multiples, where the data are used as a prediction operator, and a subtraction step, where the model of the multiples is adaptively removed from the data with matching filters. The accuracy of the multiple model depends on the source/receiver coverage at the surface. When this coverage is not dense enough, the multiple model contains errors that make successful subtraction more difficult. In these circumstances, one can either (1) improve the modelling step by interpolating the missing traces, (2) improve the subtraction step by designing methods that are less sensitive to modelling errors, or (3) both. For this data set, the second option is investigated by predicting the multiples in a 2D sense (as opposed to 3D) and performing the subtraction with a pattern‐based approach. Because some traces and shots are missing for the 2D prediction, the data are interpolated in the in‐line direction using a hyperbolic Radon transform with and without sparseness constraints. The interpolation with a sparseness constraint yields the best multiple model. For the subtraction, the pattern‐based technique is compared with a more standard, adaptive‐subtraction scheme. The pattern‐based approach is based on the estimation of 3D prediction‐error filters for the primaries and the multiples, followed by a least‐squares estimation of the primaries. Both methods are compared before and after prestack depth migration. These results suggest that, when the multiple model is not accurate, the pattern‐based method is more effective than adaptive subtraction at removing surface‐related multiples while preserving the primaries. 相似文献
12.
Prior to hydrologic modelling, topographic features of a surface are derived, and the surface is divided into sub‐basins. Surface delineation can be described as a procedure, which leads to the quantitative rendition of surface topography. Different approaches have been developed for surface delineation, but most of them may not be applicable to depression‐dominated surfaces. The main objective of this study is to introduce a new depression‐dominated delineation (D‐cubed) method and highlight its unique features by applying it to different topographic surfaces. The D‐cubed method accounts for the hierarchical relationships of depressions and channels by introducing the concept of channel‐based unit (CBU) and its connection with the concept of puddle‐based unit (PBU). This new delineation method implements a set of new algorithms to determine flow directions and accumulations for puddle‐related flats. The D‐cubed method creates a unique cascaded channel‐puddle drainage system based on the channel segmentation algorithm. To demonstrate the capabilities of the D‐cubed method, a small laboratory‐scale surface and 2 natural surfaces in North Dakota were delineated. The results indicated that the new method delineated different surfaces with and without the presence of depressional areas. Stepwise changes in depression storage and ponding area were observed for the 3 selected surfaces. These stepwise changes highlighted the dynamic filling, spilling, and merging processes of depressions, which need to be considered in hydrologic modelling for depression‐dominated areas. Comparisons between the D‐cubed method and other methods emphasized the potential consequences of use of artificial channels through the flats created by the depression‐filling process in the traditional approaches. In contrast, in the D‐cubed method, sub‐basins were further divided into a number of smaller CBUs and PBUs, creating a channel‐puddle drainage network. The testing of the D‐cubed method also demonstrated its applicability to a wide range of digital elevation model resolutions. Consideration of CBUs, PBUs, and their connection provides the opportunity to incorporate the D‐cubed method into different hydrologic models and improve their simulation of topography‐controlled runoff processes, especially for depression‐dominated areas. 相似文献
13.
Many joint inversion schemes use 1D forward modelling in the integrated interpretation of various geophysical data. In extending the joint inversion approach to the investigation of 2D structures, the discretization of the model parameters and the appropriate choice of the forward‐modelling procedure play a very important role. In this paper, a hybrid seismic–geoelectric joint inversion method is proposed for the investigation of 2D near‐surface geological structures. The electric and seismic models are coupled together through the use of common boundaries between the adjacent layers. Assuming a 2D model composed of homogeneous layers with curved boundaries, a fast ray‐tracing algorithm is used for the calculation of refraction seismic traveltime data. In the geoelectric forward modelling, a locally 1D approximation is used. The boundary surfaces are written in the form of series expansion; the inversion algorithms are formulated for the expansion coefficients and the petrophysical parameters as unknowns. Two versions of the inversion method are proposed: in versions A and B, interval‐wise constant functions and Chebyshev polynomials are, respectively, used as basis functions of the series expansion. The versions are tested by means of synthetic and in situ measured data. The tests show that both methods are stable and accurate. 相似文献
14.
The magnitude of kriging errors varies in accordance with the surface properties. The purpose of this paper is to determine the association of ordinary kriging (OK) estimated errors with the local variability of surface roughness, and to analyse the suitability of probabilistic models for predicting the magnitude of OK errors from surface parameters. This task includes determining the terrain parameters in order to explain the variation in the magnitude of OK errors. The results of this research indicate that the higher order regression models, with complex interaction terms, were able to explain 95 per cent of the variation in the OK error magnitude using the least number of predictors. In addition, the results underscore the importance of the role of the local diversity of relief properties in increasing or decreasing the magnitude of interpolation errors. The newly developed dissectivity parameters provide useful information for terrain analysis. Our study also provides constructive guides to understanding the local variation of interpolation errors and their dependence on surface dissectivity. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
15.
Xiaoxian Zeng George A. McMechan Janok P. Bhattacharya Carlos L.V. Aiken Xueming Xu William S. Hammon III Rucsandra M. Corbeanu 《Geophysical Prospecting》2004,52(3):151-163
Most existing reservoir models are based on 2D outcrop studies; 3D aspects are inferred from correlation between wells, and so are inadequately constrained for reservoir simulations. To overcome these deficiencies, we have initiated a multidimensional characterization of reservoir analogues in the Cretaceous Ferron Sandstone in Utah. Detailed sedimentary facies maps of cliff faces define the geometry and distribution of reservoir flow units, barriers and baffles at the outcrop. High‐resolution 2D and 3D ground‐penetrating radar (GPR) images extend these reservoir characteristics into 3D to allow the development of realistic 3D reservoir models. Models use geometric information from mapping and the GPR data, combined with petrophysical data from surface and cliff‐face outcrops, and laboratory analyses of outcrop and core samples. The site of the field work is Corbula Gulch, on the western flank of the San Rafael Swell, in east‐central Utah. The outcrop consists of an 8–17 m thick sandstone body which contains various sedimentary structures, such as cross‐bedding, inclined stratification and erosional surfaces, which range in scale from less than a metre to hundreds of metres. 3D depth migration of the common‐offset GPR data produces data volumes within which the inclined surfaces and erosional surfaces are visible. Correlation between fluid permeability, clay content, instantaneous frequency and instantaneous amplitude of the GPR data provides estimates of the 3D distribution of fluid permeability and clay content. 相似文献
16.
How do sediment yields from post‐wildfire debris‐laden flows depend on terrain slope,soil burn severity class,and drainage basin area? Insights from airborne‐LiDAR change detection 下载免费PDF全文
下载免费PDF全文 We derived a high‐resolution, spatially continuous map of erosion and deposition associated with the debris‐laden flows triggered by the 2011 Las Conchas wildfire and subsequent rainstorms over a 197 km2 area in New Mexico, USA. This map was produced using airborne‐LiDAR‐derived bare‐earth digital elevation models (DEMs) acquired approximately one year before and one year after the wildfire. Differencing of the pre‐wildfire and post‐wildfire‐and‐rainstorm bare‐earth DEMs yielded a DEM‐of‐difference (DoD) map that quantifies the magnitude of ground‐surface elevation changes due to erosion/deposition within each 1 m2 pixel. We applied a 0.3 m threshold filter to our DoD to remove changes that could have been due to artifacts and/or imperfect georeferencing. The 0.3 m value for the threshold filter was chosen based on the stated accuracy of the LiDAR as well as a comparison of areas of significant topographic change mapped in aerial photographs with those predicted using a range of candidate threshold values for the DoD filter. We developed an automated procedure that accepts the DoD map as input and computes, for every pixel in the DEM, the net sediment volume exported through each pixel by colluvial and/or fluvial processes using a digital routing algorithm. An analysis of the resulting sediment volume map for the Las Conchas fire demonstrates that sediment volume is proportional to upstream contributing area. After normalized by contributing area, the average sediment yield (defined as the sediment volume divided by the contributing area) increases as a power‐law functions of the average terrain slope and soil burn severity class (SBSC) with exponents equal to approximately 1.5. Our analysis quantifies the relationships among sediment yield, average terrain slope, and average soil burn severity class at the watershed scale and should prove useful for predicting the geomorphic response of wildfire‐affected drainage basins. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
17.
本文基于显式有限元法研究了地震波垂直入射时三维复杂山谷地形对地震地面运动的影响,在数值分析中应用了三维化二维的解法和黏弹性人工边界的处理方法,实现了地震波垂直输入下三维复杂场地地震动数值模拟,并验证了该方法的合理性.以四川桃坪地区一山谷地形作为研究对象,基于地表高程数据分别建立了二维和三维场地模型,对比研究表明:在复杂地形情况下考虑二、三维模型时具有明显差异,三维模型能更真实地反映地形变化对地震动的影响,复杂地形条件下有必要考虑三维实际场地模型.本文对边界自由场的处理方法也可用于处理三维复杂场地地震动斜入射问题,为三维复杂地形场地地震效应研究提供参考. 相似文献
18.
A combined algorithm for automated drainage network extraction from digital elevation models 下载免费PDF全文
下载免费PDF全文 Drainage networks are the basis for segmentation of watersheds, an essential component in hydrological modelling, biogeochemical applications, and resource management plans. With the rapidly increasing availability of topographic information as digital elevation models (DEMs), there have been many studies on DEM‐based drainage network extraction algorithms. Most of traditional drainage network extraction methods require preprocessing of the DEM in order to remove “spurious” sink, which can cause unrealistic results due to removal of real sinks as well. The least cost path (LCP) algorithm can deal with flow routing over sinks without altering data. However, the existing LCP implementations can only simulate either single flow direction or multiple flow direction over terrain surfaces. Nevertheless, terrain surfaces in the real world are usually very complicated including both convergent and divergent flow patterns. The triangular form‐based multiple flow (TFM) algorithm, one of the traditional drainage network extraction methods, can estimate both single flow and multiple flow patterns. Thus, in this paper, it is proposed to combine the advantages of the LCP algorithm and the TFM algorithm in order to improve the accuracy of drainage network extraction from the DEM. The proposed algorithm is evaluated by implementing a data‐independent assessment method based on four mathematical surfaces and validated against “true” stream networks from aerial photograph, respectively. The results show that when compared with other commonly used algorithms, the new algorithm provides better flow estimation and is able to estimate both convergent and divergent flow patterns well regarding the mathematical surfaces and the real‐world DEM. 相似文献
19.
Landform‐oriented flow‐routing algorithm for the dual‐structure loess terrain based on digital elevation models 下载免费PDF全文
下载免费PDF全文 The loess landform in the Loess Plateau of China is with typical dual structure, namely, the upper smooth positive terrain and the lower cliffy negative terrain (P–N terrain for short). Obvious differences in their morphological feature, geomorphological mechanism, and hydrological process could be found in the both areas. Based on the differences, a flow‐routing algorithm that separately addresses the dual‐structure terrain would be necessary to encompass this spatial variation in their hydrological behaviour. This paper proposes a mixed flow‐routing algorithm to address aforementioned problems. First, the loess landform surface is divided into P–N terrains based on digital elevation models. Then, specific catchment area is calculated with the new algorithm to simulate the water flows in both positive and negative terrain areas. The mixed algorithm consists of the multiple flow‐routing algorithm (multiple‐flow direction) for positive areas and the D8 algorithm for negative areas, respectively. The approach is validated in two typical geomorphologic areas with low hills and dense gullies in the northern Shaanxi Loess Plateau. Four indices are used to examine the results, which show that the new algorithm is more suitable for loess terrain in simulating the spatial distribution of water accumulation, as well as in modeling the flow characteristics of the true surface by considering the morphological structures of the terrain. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
20.
V. PEREYRA 《Geophysical Prospecting》1992,40(3):267-287