共查询到20条相似文献，搜索用时 31 毫秒
1.
DEM smoothing is a common pre-processing technique used to remove undesirable roughness from a DEM. However, it is hypothesized that smoothing straightens and reduces the length of overland flow paths, which is an important factor controlling modelled time-to-peak flow. Currently, there is a lack of research comparing how different smoothing techniques alter the distribution of overland flow path length. Four low-pass filtering techniques were applied to three fine-resolution LiDAR DEMs of varying relief: the mean filter, the median filter, the Gaussian filter, and the feature-preserving DEM smoothing (FPDEMS) filter, each with different degrees of smoothing. Downslope-distance-to-stream distributions were then derived using D8 and D∞ flow directions and statistically compared to distributions derived from the unsmoothed DEM for each study site. The results indicate that the alteration of flow path length distributions as a result of smoothing is complex. Mean flow path lengths may decrease or increase in response to smoothing, depending on landscape relief and the derivation of flow directions, and generalized flow paths may become longer. The largest increase in mean flow path lengths was 19.2 m using the 21 × 21 median filter and D8 flow directions in the high-relief study site, relative to an unsmoothed mean length of 138.6 m in this site. The largest decrease in mean flow path length was 48.9 m using the 21 × 21 mean filter and D∞ flow directions in the low-relief study site, relative to an unsmoothed mean length of 290.9 m in this site. Furthermore, minimal flow path length alterations were achieved with the Gaussian filter when gentle smoothing is required, and with the FPDEMS filter when moderate to aggressive smoothing is required. These results suggest that an appropriate smoothing method should be chosen based on the relief of the landscape and the degree of smoothing required. 相似文献
2.
Two‐dimensional nonlinear diffusive numerical simulation of geomorphic modifications to cinder cones
The temporal evolution of simple landforms such as cinder cones by nonlinear diffusive processes is studied through the use of a new 2D numerical model using well‐established and accurate numerical mathematics and high‐resolution digital elevation models (DEMs). Extending 1D (profile) nonlinear diffusion analyses used in cinder cone, hillslope and fault scarp evolution studies, we have implemented a 2D numerical model with a spatially and temporally varying sediment transport rate coefficient scaled nonlinearly by the ratio of local slope to critical slope. The high accuracy and efficient numerical implementation are documented in the paper and the MATLAB toolkit developed is used to solve for the developmentof an initial 2D cone form. First, we examine the nonlinear transport rule and suggest a refinement that accounts explicitly for flux at threshold slopes. We find that the maximum diffusion (necessarily introduced in the numerical model to avoid infinite rates) at the critical slope controls the final morphology, especially approaching steady state. Secondly, solving the landscape evolution problem in 2D enables a natural accounting for sediment flux convergence or divergence in the profile. Thirdly, the boundary behavior of a given landscape element controls much of what happens in that domain and so we allow for arbitrary flux magnitude or elevation boundary conditions. Fourthly, landscapes are heterogeneous in their surface cover and so we allow for spatially and temporally varying transport rate k and we permit an arbitrary vertical displacement field within the model domain. To test the new formulation for the nonlinear term, the effect of variable diffusivity k and the numerical schemes implemented, we apply the model to cinder cones built on the flanks of Mount Etna in 2001 and 2002–2003. We explore the effects of DEM resolution with data from the 2001 cone and the utility of spatially variable diffusivity to explain the variation in erosion measured by differencing repeat light detection and ranging (LIDAR) surveys gathered in 2004 and 2007 over the 2002–2003 cone complex. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
3.
Understanding how nitrogen transport across the landscape varies with landscape characteristics is important for developing sound nitrogen management policies. We used a spatially referenced regression analysis (SPARROW) to examine landscape characteristics influencing delivery of nitrogen from sources in a watershed to stream channels. Modelled landscape delivery ratio varies widely (by a factor of 4) among watersheds in the southeastern United States—higher in the western part (Tennessee, Alabama, and Mississippi) than in the eastern part, and the average value for the region is lower compared to other parts of the nation. When we model landscape delivery ratio as a continuous function of local‐scale landscape characteristics, we estimate a spatial pattern that varies as a function of soil and climate characteristics but exhibits spatial structure in residuals (observed load minus predicted load). The spatial pattern of modelled landscape delivery ratio and the spatial pattern of residuals coincide spatially with Level III ecoregions and also with hydrologic landscape regions. Subsequent incorporation into the model of these frameworks as regional scale variables improves estimation of landscape delivery ratio, evidenced by reduced spatial bias in residuals, and suggests that cross‐scale processes affect nitrogen attenuation on the landscape. The model‐fitted coefficient values are logically consistent with the hypothesis that broad‐scale classifications of hydrologic response help to explain differential rates of nitrogen attenuation, controlling for local‐scale landscape characteristics. Negative model coefficients for hydrologic landscape regions where the primary flow path is shallow ground water suggest that a lower fraction of nitrogen mass will be delivered to streams; this relation is reversed for regions where the primary flow path is overland flow. Published in 2009 by John Wiley & Sons, Ltd. 相似文献
4.
This paper presents improvements to the global D8 (GD8) method for calculating single flow directions in a grid digital elevation model. Flow directions computed from grid digital elevation models serve as the foundation for much of the analysis and modeling of hydrological processes that are driven by topographic gradients. The literature includes both single flow direction methods, where flow goes to only one downslope cell, and multiple flow direction methods that apportion flow among multiple downslope cells. Among single flow direction methods, the standard D8 method, in which the flow direction is set based on the steepest local slope, results in bias on surfaces that do not align with the grid directions. Efforts to address this problem have led to the development of extended methods that account for elevation values further upslope in determining flow directions. We have identified discrepancies in one such method, GD8, and have examined ways to resolve these discrepancies. An improvement to GD8, named iGD8, is presented that allows replacing a reference cell from which path deviations are accumulated and that considers horizontal path deviation rather than global slope as a flow direction criterion. The improved method is found to be effective in resolving the problems encountered with GD8 and to be more efficient than a previously proposed alternative method (least transversal deviation (LTD) based D8, namely D8‐LTD) that uses recursive searching for the largest upstream area when multiple flow paths converge. The proposed improved GD8 method offers the opportunity for improved analysis and modeling of topographically driven hydrological processes by providing better foundational flow directions for these analyses. 相似文献
5.
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. 相似文献
6.
Linking landscape patterns to specific ecological processes has been and will continue to be the key topic in landscape ecology.However,this goal is difficult to achieve by using the traditional landscape metric based on the Patch-Mosaic Model(PMM),as they don't integrate ecological processes with landscape patterns.In this paper,we proposed a conceptual model,i.e.,the Source-Pathway-Sink Model(SPSM),which designates the role of a landscape unit into "source", "sink",or "pathway" based on specific ecological processes during the landscape pattern analysis.While the traditional landscape metrics derived from the PMM model is visual-or geometrical-oriented and lack of linkage to ecological significance,the SPSM model is process-oriented,dynamic,and scale dependent.A comparison between the PMM and the SPSM models shows that the SPSM model is complementary to the PMM model,and can provide a simple and dynamic perspective on landscape pattern analysis.The SPSM model may represent a conceptual innovation in landscape ecology. 相似文献
7.
This work is the third in a series of papers on the thermodynamically constrained averaging theory (TCAT) approach to modeling flow and transport phenomena in multiscale porous medium systems. Building upon the general TCAT framework and the mathematical foundation presented in previous works in this series, we demonstrate the TCAT approach for the case of single-fluid-phase flow. The formulated model is based upon conservation equations for mass, momentum, and energy and a general entropy inequality constraint, which is developed to guide model closure. A specific example of a closed model is derived under limiting assumptions using a linearization approach and these results are compared and contrasted with the traditional single-phase-flow model. Potential extensions to this work are discussed. Specific advancements in this work beyond previous averaging theory approaches to single-phase flow include use of macroscale thermodynamics that is averaged from the microscale, the use of derived equilibrium conditions to guide a flux–force pair approach to simplification, use of a general Lagrange multiplier approach to connect conservation equation constraints to the entropy inequality, and a focus on producing complete, closed models that are solvable. 相似文献
8.
Watershed delineation is a required step when conducting any spatially distributed hydrological modelling. Automated approaches are often proposed to delineate a watershed based on a river network extracted from the digital elevation model (DEM) using the deterministic eight‐neighbour (D8) method. However, a realistic river network cannot be derived from conventional DEM processing methods for a large flat area with a complex network of rivers, lakes, reservoirs, and polders, referred to as a plain river network region (PRNR). In this study, a new approach, which uses both hydrographic features and DEM, has been developed to address the problems of watershed delineation in PRNR. It extracts the river nodes and determines the flow directions of the river network based on a vector‐based hydrographic feature data model. The river network, lakes, reservoirs, and polders are then used to modify the flow directions of grid cells determined by D8 approach. The watershed is eventually delineated into four types of catchments including lakes, reservoirs, polders, and overland catchments based on the flow direction matrix and the location of river nodes. Multiple flow directions of grid cells are represented using a multi‐direction encoding method, and multiple outflows of catchments are also reflected in the topology of catchments. The proposed approach is applied to the western Taihu watershed in China. Comparisons between the results obtained from the D8 approach, the ‘stream burning’ approach, and those from the proposed approach clearly demonstrate an improvement of the new approach over the conventional approaches. This approach will benefit the development of distributed hydrological models in PRNR for the consideration of different types and multiple inlets and outlets of catchments. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
9.
Raphaël Bunel Nicolas Lecoq Yoann Copard Eric Guérin Marco Van de Wiel Nicolas Massei 《地球表面变化过程与地形》2021,46(3):593-610
Understanding, analysing, and predicting the erosion mechanisms and sedimentary flows produced by catchments plays a key role in environmental conservation and restoration management and policies. Numerical case-testing studies are generally undertaken to analyse the sensitivity of flood and soil erosion processes to the physical characteristics of catchments. Most analyses are conducted on simple virtual catchments with physical characteristics that, unlike real catchments, are perfectly controlled. Virtual catchments generally correspond to V-shaped valley catchments. However, although these catchments are suitable for methodical analysis of the results, they do not provide a realistic representation of the spatial structures of the landscape and field conditions. They can, therefore, lead to potential modelling errors and can make it difficult to extend or generalize their results. Our proposed method bridges the gap between real and traditional virtual catchments by creating realistic virtual catchments with perfectly controllable physical characteristics. Our approach represents a real alternative to traditional test case procedures and provides a new framework for geomorphological and hydrological communities. It combines a field procedural generation approach, geographic information system processing procedures, and the CAESAR-Lisflood landscape evolution model. We illustrate how each of these components acts in the process of generating virtual catchments. Five physical parameters were adjusted and tested for each virtual catchment: drainage density, hypsometric integral, mean slope of the main channel, granulometry, and land use. One of our virtual catchments is compared with a real catchment and a virtual catchment produced by a standard method. This comparison indicates that our approach can produce more realistic virtual catchments than those produced by more traditional methods, while a high degree of controllability is maintained. This new method of generating virtual catchments therefore offers significant research potential to identify the impacts of the physical characteristics of catchments on hydro-sedimentary dynamics and responses. 相似文献
10.
Generation of permeability field in a reservoir model that matchs historical dynamic production data requires an inverse calculation. A gradient method is typically used to solve the inverse minimization problem and requires sensitivity coefficients of reservoir responses, e.g. fractional flow rate or pressure, with respect to the change in the permeability. This paper presents a novel semi-analytical streamline-based method for computing such sensitivity coefficients under the framework of two-phase (oil-water) flow conditions. This method is shown to be significantly faster and generate permeability fields with lower objective function than the traditional perturbation method. The method decomposes the multiple-dimensional full flow problem into multiple 1D problems along streamlines. The sensitivity of fractional flow rate at the production well is directly related to the sensitivity of time-of-flight (TOF) along each individual streamline and the sensitivity of pressure at grid cells along the streamline. The sensitivity of TOF of a streamline can be obtained analytically. The sensitivity of pressure is obtained as part of a fast single phase flow simulation. The proposed method is implemented in a geostatistically based inverse technique, called the sequential self-calibration (SSC) method. Results for fractional flow rate sensitivities are presented and compared with the traditional perturbation method. This new method can be easily extended to compute sensitivity coefficients of saturation (concentration) data. 相似文献
11.
A parameterized stochastic model of near‐fault ground motion in two orthogonal horizontal directions is developed. The major characteristics of recorded near‐fault ground motions are represented. These include near‐fault effects of directivity and fling step; temporal and spectral non‐stationarity; intensity, duration, and frequency content characteristics; directionality of components; and the natural variability of ground motions. Not all near‐fault ground motions contain a forward directivity pulse, even when the conditions for such a pulse are favorable. The proposed model accounts for both pulse‐like and non‐pulse‐like cases. The model is fitted to recorded near‐fault ground motions by matching important characteristics, thus generating an ‘observed’ set of model parameters for different earthquake source and site characteristics. A method to generate and post‐process synthetic motions for specified model parameters is also presented. Synthetic ground motion time series are generated using fitted parameter values. They are compared with corresponding recorded motions to validate the proposed model and simulation procedure. The use of synthetic motions in addition to or in place of recorded motions is desirable in performance‐based earthquake engineering applications, particularly when recorded motions are scarce or when they are unavailable for a specified design scenario. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
12.
Application of distributed hydrological models for predictions in ungauged basins: a method to quantify predictive uncertainty 下载免费PDF全文
下载免费PDF全文 Stream flow predictions in ungauged basins are one of the most challenging tasks in surface water hydrology because of nonavailability of data and system heterogeneity. This study proposes a method to quantify stream flow predictive uncertainty of distributed hydrologic models for ungauged basins. The method is based on the concepts of deriving probability distribution of model's sensitive parameters by using measured data from a gauged basin and transferring the distribution to hydrologically similar ungauged basins for stream flow predictions. A Monte Carlo simulation of the hydrologic model using sampled parameter sets with assumed probability distribution is conducted. The posterior probability distributions of the sensitive parameters are then computed using a Bayesian approach. In addition, preselected threshold values of likelihood measure of simulations are employed for sizing the parameter range, which helps reduce the predictive uncertainty. The proposed method is illustrated through two case studies using two hydrologically independent sub‐basins in the Cedar Creek watershed located in Texas, USA, using the Soil and Water Assessment Tool (SWAT) model. The probability distribution of the SWAT parameters is derived from the data from one of the sub‐basins and is applied for simulation in the other sub‐basin considered as pseudo‐ungauged. In order to assess the robustness of the method, the numerical exercise is repeated by reversing the gauged and pseudo‐ungauged basins. The results are subsequently compared with the measured stream flow from the sub‐basins. It is observed that the measured stream flow in the pseudo‐ungauged basin lies well within the estimated confidence band of predicted stream flow. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
13.
Integrated river basin models should provide a spatially distributed representation of basin hydrology and transport processes to allow for spatially implementing specific management and conservation measures. To accomplish this, the Soil and Water Assessment Tool (SWAT) was modified by integrating a landscape routing model to simulate water flow across discretized routing units. This paper presents a grid‐based version of the SWAT landscape model that has been developed to enhance the spatial representation of hydrology and transport processes. The modified model uses a new flow separation index that considers topographic features and soil properties to capture channel and landscape flow processes related to specific landscape positions. The resulting model is spatially fully distributed and includes surface, lateral and groundwater fluxes in each grid cell of the watershed. Furthermore, it more closely represents the spatially heterogeneous distributed flow and transport processes in a watershed. The model was calibrated and validated for the Little River Watershed (LRW) near Tifton, Georgia (USA). Water balance simulations as well as the spatial distribution of surface runoff, subsurface flow and evapotranspiration are examined. Model results indicate that groundwater flow is the dominant landscape process in the LRW. Results are promising, and satisfactory output was obtained with the presented grid‐based SWAT landscape model. Nash–Sutcliffe model efficiencies for daily stream flow were 0.59 and 0.63 for calibration and validation periods, and the model reasonably simulates the impact of the landscape position on surface runoff, subsurface flow and evapotranspiration. Additional revision of the model will likely be necessary to adequately represent temporal variations of transport and flow processes in a watershed. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
14.
Peter M. Allen Jeffrey G. Arnold Lance Auguste Joseph White John Dunbar 《地球表面变化过程与地形》2018,43(1):202-217
Gully erosion begins in streambanks and uplands as a consequence of adjustments in driving forces on the landscape imposed by changes in land use or climate. The deleterious effects of gullies worldwide have led to many site‐specific studies of gully form and function. In the continental United States, gully erosion in agricultural land has destroyed valuable farmland yet, prediction of gully processes remains problematic on a national scale. This research has proposed a simple method to predict gully headcut advance. When combined with SWAT hydrologic flow routines, the model predicted gully headcut advance with reasonable accuracy on a daily time step for time periods exceeding two decades. The model was tested in two distinct land resource areas of the United States with differing climate, soils, cover and drainage. The inputs for the headcut model have been kept simple as the model will be applied over large areas. Model inputs consist of headcut height, headcut resistance (based on soil erodibility and a root‐cover factor), and daily flow. The model is compared with an annual time step model used in assessment of headcut advance and appears to offer a better way to assess gully headcut advance. Copyright © 2017 John Wiley & Sons, Ltd. 相似文献
15.
Suspended sediment yield modeling in Mahanadi River,India by multi-objective optimization hybridizing artificial intelligence algorithms 总被引：2，自引：1，他引：1
River sediment produced through weathering is one of the principal landscape modification processes on earth.Rivers are an integral part of the hydrologic cycle and are the major geologic agents that erode the continents and transport water and sediments to the oceans.Estimation of suspended sediment yield is always a key parameter for planning and management of any river system.It is always challenging to model sediment yield using traditional mathematical models because they are incapable of handling the complex non-linearity and non-stationarity.The suspended sediment modeling of the river depends on the number of factors such as rock type,relief,rainfall,temperature,water discharge and catchment area.In this study,we proposed a hybrid genetic algorithm-based multi-objective optimization with artificial neural network(GA-MOO-ANN)with automated parameter tuning model using these factors to estimate the suspended sediment yield in the entire Mahanadi River basin.The model was validated by comparing statistically with other models,and it appeared that the GA-MOO-ANN model has the lowest root mean squared error(0.009)and highest coefficient of correlation(0.885)values among all comparative models(traditional neural network,multiple linear regression,and sediment rating curve)for all stations.It was also observed that the proposed model is the least biased(0.001)model.Thus,the proposed GA-MOOANN is the most capable model,compared to other studied models,for estimating the suspended sediment yield in the entire Mahanadi river basin,India.The results also suggested that the proposed GA-MOO-ANN model is unable to estimate suspended sediment yield satisfactorily at gauge stations having very small catchment areas whereas performing satisfactorily on locations having moderate to the large catchment area.The models provide the best result at Tikarapara,the gauge station location in the extreme downstream,having the largest catchment area. 相似文献
16.
A laboratory study of the rheology of mudflows in Hangzhou Bay, China, is reported in this paper. Both the steady and oscillatory （dynamic） rheological properties are studied using RMS-605 rheometer. A Dual-Bingham model is proposed for analyzing flow curves and compared with Worrall-Tuliani model. It is found that Dual-Bingham plastic rheological model is easier to implement than Worrall-Tuliani model and can provide satisfactory representations of the steady mudflows in Hangzhou Bay and other published data. The dependence of the yield stress and viscosity on sediment concentration is discussed based on the data from Hangzhou Bay mud and other published data. For the dynamic rheologieal properties of Hangzhou Bay mud, empirical expressions for elastic modulus and dynamic viscosity are provided in the form of exponential functions of sediment volume concentration, and comparisons with other published data also discussed. 相似文献
17.
Modeling the spread of subsurface contaminants requires coupling a groundwater flow model with a contaminant transport model. Such coupling may provide accurate estimates of future subsurface hydrologic states if essential flow and contaminant data are assimilated in the model. Assuming perfect flow, an ensemble Kalman filter (EnKF) can be used for direct data assimilation into the transport model. This is, however, a crude assumption as flow models can be subject to many sources of uncertainty. If the flow is not accurately simulated, contaminant predictions will likely be inaccurate even after successive Kalman updates of the contaminant model with the data. The problem is better handled when both flow and contaminant states are concurrently estimated using the traditional joint state augmentation approach. In this paper, we introduce a dual estimation strategy for data assimilation into a one-way coupled system by treating the flow and the contaminant models separately while intertwining a pair of distinct EnKFs, one for each model. The presented strategy only deals with the estimation of state variables but it can also be used for state and parameter estimation problems. This EnKF-based dual state-state estimation procedure presents a number of novel features: (i) it allows for simultaneous estimation of both flow and contaminant states in parallel; (ii) it provides a time consistent sequential updating scheme between the two models (first flow, then transport); (iii) it simplifies the implementation of the filtering system; and (iv) it yields more stable and accurate solutions than does the standard joint approach. We conducted synthetic numerical experiments based on various time stepping and observation strategies to evaluate the dual EnKF approach and compare its performance with the joint state augmentation approach. Experimental results show that on average, the dual strategy could reduce the estimation error of the coupled states by 15% compared with the joint approach. Furthermore, the dual estimation is proven to be very effective computationally, recovering accurate estimates at a reasonable cost. 相似文献
18.
Scott W. Anderson 《地球表面变化过程与地形》2019,44(5):1015-1033
Topographic surveys inevitably contain error, introducing uncertainty into estimates of volumetric or mean change based on the differencing of repeated surveys. In the geomorphic community, uncertainty has often been framed as a problem of separating out real change from apparent change due purely to error, and addressed by removing measured change considered indistinguishable from random noise from analyses (thresholding). Thresholding is important when quantifying gross changes (i.e. total erosion or total deposition), which are systematically biased by random errors in stable parts of a landscape. However, net change estimates are not substantially influenced by those same random errors, and the use of thresholds results in inherently biased, and potentially misleading, estimates of net change and uncertainty. More generally, thresholding is unrelated to the important process of propagating uncertainty in order to place uncertainty bounds around final estimates. Error propagation methods for uncorrelated, correlated, and systematic errors are presented. Those equations demonstrate that uncertainties in modern net change analyses, as well as in gross change analyses using reasonable thresholds, are likely to be dominated by low-magnitude but highly correlated or systematic errors, even after careful attempts to reduce those errors. In contrast, random errors with little to no correlation largely cancel to negligible levels when averaged or summed. Propagated uncertainty is then typically insensitive to the precision of individual measurements, and is instead defined by the relative mean error (accuracy) over the area of interest. Given that real-world mean elevation changes in many landscape settings are often similar in magnitude to potential mean errors in repeat topographic analyses, reducing highly correlated or systematic errors will be central to obtaining accurate change estimates, while placing uncertainty bounds around those results provides essential context for their interpretation. Published 2018. This article is a U.S. Government work and is in the public domain in the USA. 相似文献
19.
G. R. Hancock 《地球表面变化过程与地形》2008,33(3):459-474
The current generation of landscape evolution models use a digital elevation model for landscape representation. These programs also contain a hydrological model that defines overland flow with the drainage network routed to an outlet. One of the issues with landscape evolution modelling is the hydrological correctness of the digital elevation model used for the simulations. Despite the wide use and increased quality of digital elevation models, data pits and depressions in the elevation data are a common feature and their removal will remain a necessary step for many data sets. This study examines whether a digital elevation model can be hydrologically correct (i.e. all depressions removed so that all water can run downslope) before use in a landscape evolution model and what effect depression removal has on long‐term geomorphology and hydrology. The impact on sediment transport rates is also examined. The study was conducted using a field catchment and a proposed landform for a post‐mining landscape. The results show that there is little difference in catchment geomorphology and hydrology for the non‐depression removed and depression removed data sets. The non‐depression removed and depression removed digital elevation models were also evaluated as input to a landscape evolution model for a 50 000 year simulation period. The results show that after 1000 years there is little difference between the data sets, although sediment transport rates did vary considerably early on in the simulation. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
20.
This work introduces water–air two‐phase flow into integrated surface–subsurface flow by simulating rainfall infiltration and run‐off production on a soil slope with the finite element method. The numerical model is formulated by partial differential equations for hydrostatic shallow flow and water–air two‐phase flow in the shallow subsurface. Finite element computing formats and solution strategies are presented to obtain a numerical solution for the coupled model. An unsaturated seepage flow process is first simulated by water–air two‐phase flow under the atmospheric pressure boundary condition to obtain the rainfall infiltration rate. Then, the rainfall infiltration rate is used as an input parameter to solve the surface run‐off equations and determine the value of the surface run‐off depth. In the next iteration, the pressure boundary condition of unsaturated seepage flow is adjusted by the surface run‐off depth. The coupling process is achieved by updating the rainfall infiltration rate and surface run‐off depth sequentially until the convergence criteria are reached in a time step. A well‐conducted surface run‐off experiment and traditional surface–subsurface model are used to validate the new model. Comparisons with the traditional surface–subsurface model show that the initiation time of surface run‐off calculated by the proposed model is earlier and that the water depth is larger, thus providing values that are closer to the experimental results. 相似文献