Browsing by Author "Orlandini, S"

Now showing 1 - 9 of 9

Open Access
Analytical basis for determining slope lines in grid digital elevation models
(Water Resources Research, 2014-01-01) Orlandini, S; Moretti, G; Gavioli, A
An analytical basis for the determination of slope lines in grid digital elevation models is provided by using the D8-LTD method (eight slope directions, least transverse deviation). The D8-LTD method's capability to predict consistently exact slope lines as the grid cell size goes to zero is shown analytically by applying mathematical analysis methods. The use of cumulative, least transverse deviations is found to be the key factor allowing for globally unbiased approximations of slope lines. The D8-LTD method's properties are also demonstrated numerically by using digital elevation models of a synthetic sloping surface obtained from the Himmelblau function. It is shown that slope lines obtained from the D8-LTD method can approximate the exact slope lines as close as desired by selecting a grid cell size that is small enough. In contrast, the standard D8 method is found to produce significantly biased results even when small grid cells are used. The D8-LTD method outperforms the D8 method over a wide range of grid cell sizes (up to 20 m in this application), beyond which grid cell size becomes too large to validly represent the underlying sloping surface. It is therefore concluded that the D8-LTD method should be used in preference to the standard D8 method in order to obtain slope lines that are only limited in reliability by the detail of topographic data, and not by the accuracy of the slope direction method applied. Key Points The standard D8 method is inadequate to predict slope lines The D8-LTD method predicts consistently exact slope lines Dispersive methods provide patterns developing around slope lines ©2013. American Geophysical Union. All Rights Reserved.
Open Access
Automatic delineation of drainage basins from contour elevation data using skeleton construction techniques
(Water Resources Research, 2008-05-01) Moretti, G; Orlandini, S
New methods for automatic delineation of drainage basins from contour elevation data are presented. As a fundamental preprocessing step, the points defining a set of contour lines are used to compute the Delaunay triangulation, the Voronoi diagram, and other structures known in computational geometry as the crust and the skeleton (or medial axis transform). By exploiting the skeleton extracted from contour lines, a recursive algorithm is then developed to solve critical topographic structures such as ridges, saddles, and peaks in a fully automated and accurate manner. Finally, the algorithm is further extended to deal with the construction of flow nets. Numerical experiments based on high-accuracy contour elevation data of real terrains show that the proposed methods are able to process automatically complex topographic structures and to produce results comparable to those that can be interpreted visually from contour lines. The gain in accuracy over current state-of-the-art solutions is generally found to be significant and to increase as the contour interval increases. Copyright 2008 by the American Geophysical Union.
Open Access
Control of coupling mass balance error in a process-based numerical model of surface-subsurface flow interaction
(Water Resources Research, 2015-07-01) Fiorentini, M; Orlandini, S; Paniconi, C
© 2015. American Geophysical Union. All Rights Reserved.A process-based numerical model of integrated surface-subsurface flow is analyzed in order to identify, track, and reduce the mass balance errors affiliated with the model's coupling scheme. The sources of coupling error include a surface-subsurface grid interface that requires node-to-cell and cell-to-node interpolation of exchange fluxes and ponding heads, and a sequential iterative time matching procedure that includes a time lag in these same exchange terms. Based on numerical experiments carried out for two synthetic test cases and for a complex drainage basin in northern Italy, it is shown that the coupling mass balance error increases during the flood recession limb when the rate of change in the fluxes exchanged between the surface and subsurface is highest. A dimensionless index that quantifies the degree of coupling and a saturated area index are introduced to monitor the sensitivity of the model to coupling error. Error reduction is achieved through improvements to the heuristic procedure used to control and adapt the time step interval and to the interpolation algorithm used to pass exchange variables from nodes to cells. The analysis presented illustrates the trade-offs between a flexible description of surface and subsurface flow processes and the numerical errors inherent in sequential iterative coupling with staggered nodal points at the land surface interface, and it reveals mitigation strategies that are applicable to all integrated models sharing this coupling and discretization approach.
Open Access
Determination of surface flow paths from gridded elevation data
(Water Resources Research, 2009-03-01) Orlandini, S; Moretti, G
[1] Surface flow paths are obtained from gridded elevation data by connecting grid cell centers along predetermined flow directions. These flow directions are commonly determined using single and multiple flow direction algorithms. It remains, however, unclear whether multiple flow direction algorithms, which introduce artificial dispersion, can be used to describe surface flow paths and gravity-driven processes across a terrain without causing unrealistic flow dispersion. To explore this issue, a unified algorithm for the determination of flow directions has been developed, and new methods for the validation of the resulting surface flow paths are introduced. The unified algorithm makes it possible, by setting appropriate parameters, to perform local or path-based analyses and to experiment with different combinations of single and multiple flow directions in a morphologically significant manner. The new validation methods use drainage systems delineated from contour elevation data as a reference and take into consideration the overlap between these systems and those obtained from gridded elevation data. The unified algorithm is presented, and the results are evaluated for selected case studies in order to provide guidance on the use of surface flow path algorithms based on gridded elevation data. © 2009 by the American Geophysical Union.
Open Access
Evaluation of flow direction methods against field observations of overland flow dispersion
(Water Resources Research, 2012-10-29) Orlandini, S; Moretti, G; Corticelli, MA; Santangelo, PE; Capra, A; Rivola, R; Albertson, JD
The D8, D8-LTD, D-LTD, D, MD, and MD8 flow direction methods are evaluated against field observations of overland flow dispersion obtained from novel experimental methods. Thin flows of cold water were released at selected points on a warmer slope and individual overland flow patterns originating from each of these points were observed using a terrestrial laser scanner and a thermal imaging camera. Land microtopography was determined by using laser returns from the dry land surface, whereas overland flow patterns were determined by using either laser returns or infrared emissions from the wetted portions of the land surface. Planar overland flow dispersion is found to play an important role in the region lying immediately downslope of the point source, but attenuates rapidly as flow propagates downslope. In contrast, existing dispersive flow direction methods are found to provide a continued dispersion with distance downslope. Predicted propagation patterns, for all methods considered here, depend critically on the size h of grid cells involved. All methods are found to be poorly sensitive in extremely fine grids (h 2 cm), and to be poorly specific in coarse grids (h = 2 m). Satisfactory results are, however, obtained in grids having resolutions h that approach the average flow width (50 cm), with the best performances displayed by the MD8 method in the finest grids (5 h 20 cm), and by the MD, D, and D-LTD methods in the coarsest grids (20 cm < h 1 m). © 2012. American Geophysical Union. All Rights Reserved.
Open Access
Evidence of an emerging levee failure mechanism causing disastrous floods in Italy
(Water Resources Research, 2015-10-01) Orlandini, S; Moretti, G; Albertson, JD
© 2015. The Authors.A levee failure occurred along the Secchia River, Northern Italy, on 19 January 2014, resulting in flood damage in excess of $500 million. In response to this failure, immediate surveillance of other levees in the region led to the identification of a second breach developing on the neighboring Panaro River, where rapid mitigation efforts were successful in averting a full levee failure. The paired breach events that occurred along the Secchia and Panaro Rivers provided an excellent window on an emerging levee failure mechanism. In the Secchia River, by combining the information content of photographs taken from helicopters in the early stage of breach development and 10 cm resolution aerial photographs taken in 2010 and 2012, animal burrows were found to exist in the precise levee location where the breach originated. In the Panaro River, internal erosion was observed to occur at a location where a crested porcupine den was known to exist and this erosion led to the collapse of the levee top. This paper uses detailed numerical modeling of rainfall, river flow, and variably saturated flow in the levee to explore the hydraulic and geotechnical mechanisms that were triggered along the Secchia and Panaro Rivers by activities of burrowing animals leading to levee failures. As habitats become more fragmented and constrained along river corridors, it is possible that this failure mechanism could become more prevalent and, therefore, will demand greater attention in both the design and maintenance of earthen hydraulic structures as well as in wildlife management.
Open Access
On the prediction of channel heads in a complex alpine terrain using gridded elevation data
(Water Resources Research, 2011-03-09) Orlandini, S; Tarolli, P; Moretti, G; Dalla Fontana, G
Threshold conditions for channel initiation are evaluated by using gridded elevation data derived from a lidar survey, a reliable algorithm for the determination of surface flow paths, and field observations of channel heads for a study area located in the eastern Italian Alps. These threshold conditions are determined by considering the channel heads observed across a portion of the study area and computing the related values of (1) drainage area A, (2) area-slope function AS2, with S being the local slope, and (3) Strahler order ω* of surface flow paths extracted from gridded elevation data. Attention is focused on the dependence of the obtained threshold values on the size of grid cells involved and on the ability of the identified threshold conditions to provide reliable predictions of channel heads across the entire study area. The results indicate that the threshold values of A, AS2, and ω* are all significantly dependent on grid cell size, and the uncertainty in the determination of threshold values of ω* is significantly smaller than that affecting the determination of threshold values of A and AS2. The comparison between predicted and observed channel heads indicates that the considered methods display variable reliability and sensitivity over different drainage basins and grid cell sizes, with a general tendency to predict more channel heads than can be observed in the field. Acceptable predictions are normally obtained where channel heads are formed essentially by surface runoff. More comprehensive methods seem, however, to be needed to predict channel heads affected by groundwater seeping upward. Copyright 2011 by the American Geophysical Union.
Open Access
Robust numerical solution of the reservoir routing equation
(Advances in Water Resources, 2013-09-01) Fiorentini, M; Orlandini, S
The robustness of numerical methods for the solution of the reservoir routing equation is evaluated. The methods considered in this study are: (1) the Laurenson-Pilgrim method, (2) the fourth-order Runge-Kutta method, and (3) the fixed order Cash-Karp method. Method (1) is unable to handle nonmonotonic outflow rating curves. Method (2) is found to fail under critical conditions occurring, especially at the end of inflow recession limbs, when large time steps (greater than 12. min in this application) are used. Method (3) is computationally intensive and it does not solve the limitations of method (2). The limitations of method (2) can be efficiently overcome by reducing the time step in the critical phases of the simulation so as to ensure that water level remains inside the domains of the storage function and the outflow rating curve. The incorporation of a simple backstepping procedure implementing this control into the method (2) yields a robust and accurate reservoir routing method that can be safely used in distributed time-continuous catchment models. © 2013 Elsevier Ltd.
Open Access
Surface-subsurface flow modeling with path-based runoff routing, boundary condition-based coupling, and assimilation of multisource observation data
(Water Resources Research, 2010-02-01) Camporese, M; Paniconi, C; Putti, M; Orlandini, S
A distributed physically based model incorporating novel approaches for the representation of surface-subsurface processes and interactions is presented. A path-based description of surface flow across the drainage basin is used, with several options for identifying flow directions, for separating channel cells from hillslope cells, and for representing stream channel hydraulic geometry. Lakes and other topographic depressions are identified and specially treated as part of the preprocessing procedures applied to the digital elevation data for the catchment. Threshold-based boundary condition switching is used to partition potential (atmospheric) fluxes into actual fluxes across the land surface and changes in surface storage, thus resolving the exchange fluxes, or coupling, between the surface and subsurface modules. Nested time stepping allows smaller steps to be taken for typically faster and explicitly solved surface runoff routing, while a mesh coarsening option allows larger grid elements to be used for typically slower and more compute-intensive subsurface flow. Sequential data assimilation schemes allow the model predictions to be updated with spatiotemporal observation data of surface and subsurface variables. These approaches are discussed in detail, and the physical and numerical behavior of the model is illustrated over catchment scales ranging from 0.0027 to 356 km2, addressing different hydrological processes and highlighting the importance of describing coupled surface-subsurface flow. Copyright 2010 by the American Geophysical Union.