Flood frequency analysis based on simulated peak discharges

Flood frequency approaches vary from statistical methods, directly applied on the observed annual maximum flood series, to adopting rainfall–runoff simulation models that transform design rainfalls to flood discharges. Reliance on statistical flood frequency analysis depends on several factors such as the selected probability distribution function, estimation of the function parameters, possible outliers, and length of the observed flood series. Through adopting the simulation approach in this paper, watershed-average rainfalls of various occurrence probabilities were transformed into the corresponding peak discharges using a calibrated hydrological model. A Monte Carlo scheme was employed to consider the uncertainties involved in rainfall spatial patterns and antecedent soil moisture condition (AMC). For any given rainfall depth, realizations of rainfall spatial distribution and AMC conditions were entered as inputs to the model. Then, floods of different return periods were simulated by transforming rainfall to runoff. The approach was applied to Tangrah watershed in northeastern Iran. It was deduced that the spatial rainfall distribution and the AMCs exerted a varying influence on the peak discharge of different return periods. Comparing the results of the simulation approach with those of the statistical frequency analysis revealed that, for a given return period, flood quantiles based on the observed series were greater than the corresponding simulated discharges. It is also worthy to note that existence of outliers and the selection of the statistical distribution function has a major effect in increasing the differences between the results of the two approaches.     

Regional hydrologic mapping of flows in stream networks

River low flow (LF) is an important hydrological characteristic used in management of the quantity and quality of water resources. A common way to transpose at-station low flow quantiles to ungaged locations involves development of regional multivariate regression models. In site selection studies of water control facilities, need arises to automate the mapping of LF characteristics onto the stream network in a continuous manner. This would allow estimation of LF at any desired ungaged location. A mapping algorithm is introduced that determines the value of regression model's independent variables (input components) over the drainage area of each stream pixel and calculates the flow (output). The mapping algorithm relies mainly on the digital elevation model (DEM) and its derivatives, such as flow direction and flow accumulation. Furthermore, the contribution of each independent variable of the regional model to the total flow may be plotted to represent the flow profiles along the streams. A case study involving LF mapping in rivers of Gilan province, Iran, where LF quantiles are required for issuing water withdrawal permits as well as maintaining water quality standards, is also described in the paper. Overall, the maps and profiles of flow statistics in the region of interest provide convenient visualization and assessment tools for water resource and environmental engineers. The transposition algorithm may also be applied in regional mapping of other flow characteristics such as flood or average flows.     

Analysis of threshold and incipient conditions for sediment movement

Prediction of threshold conditions and incipient motion is the essential issue for the study of sediment transport. This work compares existing empirical threshold curves proposed for Shields diagram, a method based on the concept of probability of sediment movement, and an empirical method based on movability number. These methods are used to predict the incipient motion conditions for experimental runs taken from various studies. Most of the experimental data, used in this work, have not been used before in derivation of alternative formulations for Shields diagram and other methods. The empirical threshold curves based on the Shields entrainment function was the least successful at predicting the measured incipient motion conditions, while the use of the movability number gives good predictions of critical shear velocity compared with experimental data.     

Issues in Eulerian–Lagrangian modeling of sediment transport under saltation regime

The saltation regime is very important for understanding the sediment transport mechanism. However,there is no consensus on a model for the saltation regime. This study answers several questions raised with respect to the Eulerian-Lagrangian modeling of sediment transport. The first question is why the previous saltation models that use different combinations of hydrodynamic forces yielded acceptable results? The second question is which shear lift model(i.e. a shear lift expression and its coefficient) is more appropriate? Another important question is which hydrodynamic forces have greater contributions to the saltation characteristics of a sediment particle? The last question is what are the contributions of the turbulence fluctuations as well as effects of using two-and three-dimensional(2 D and 3 D) models on the simulation results? In order to fairly answer these questions, a systematic study was done by considering different scenarios. The current study is the first attempt to clearly discuss these issues. A comprehensive 3 D saltation model for non-cohesive sediment was developed that includes all the hydrodynamic forces acting on the particle. The random nature of sediment transport was included using turbulent flow and bed-particle collision models. The eddy interaction model was applied to generate a3 D turbulent flow field. Bed-particle collisions were considered using the concept of a contact zone and a corresponding contact point. The validation of the model was done using the available experimental data for a wide range of sediment size(0.03 to 4.8 cm). For the first question, the results indicated that some of the hydrodynamic effects show opposing trends and some have negligible effects. With these opposing effects it is possible to adjust the coefficients of different models to achieve acceptable agreement with the same experimental data while omitting some aspects of the physics of the process. A suitable model for the shear lift force was developed by linking the lift coefficient to the drag coefficient and the contributions of the hydrodynamic forces and turbulence fluctuations as well as the consequences of using of 2 D and 3 D models were studied. The results indicate that the shear lift force and turbulent flow fluctuations are important factors for the saltation of both sand and gravel, and they cannot be ignored.     

Assessment of geomorphological bank evolution of the alluvial threshold rivers based on entropy concept parameters

The complex stream bank profiles in alluvial channels and rivers that are formed after reaching equilibrium has been a popular topic of research for many geomorphologists and river engineers. The entropy theory has recently been successfully applied to this problem. However, the existing methods restrict the further application of the entropy parameter to determine the cross-section slope of the river banks. To solve this limitation, we introduce a novel approach in the extraction of the equation based on the calculation of the entropy parameter (λ) and the transverse slope of the bank profile at threshold channel conditions. The effects of different hydraulic and geometric parameters are evaluated on a variation of the entropy parameter. Sensitivity analysis on the parameters affecting the entropy parameter shows that the most effective parameter on the λ-slope multiplier is the maximum slope of the bank profile and the dimensionless lateral distance of the river banks.     

Evaluation of four supervised learning methods for groundwater spring potential mapping in Khalkhal region (Iran) using GIS-based features

One important tool for water resources management in arid and semi-arid areas is groundwater potential mapping. In this study, four data-mining models including K-nearest neighbor (KNN), linear discriminant analysis (LDA), multivariate adaptive regression splines (MARS), and quadric discriminant analysis (QDA) were used for groundwater potential mapping to get better and more accurate groundwater potential maps (GPMs). For this purpose, 14 groundwater influence factors were considered, such as altitude, slope angle, slope aspect, plan curvature, profile curvature, slope length, topographic wetness index (TWI), stream power index, distance from rivers, river density, distance from faults, fault density, land use, and lithology. From 842 springs in the study area, in the Khalkhal region of Iran, 70 % (589 springs) were considered for training and 30 % (253 springs) were used as a validation dataset. Then, KNN, LDA, MARS, and QDA models were applied in the R statistical software and the results were mapped as GPMs. Finally, the receiver operating characteristics (ROC) curve was implemented to evaluate the performance of the models. According to the results, the area under the curve of ROCs were calculated as 81.4, 80.5, 79.6, and 79.2 % for MARS, QDA, KNN, and LDA, respectively. So, it can be concluded that the performances of KNN and LDA were acceptable and the performances of MARS and QDA were excellent. Also, the results depicted high contribution of altitude, TWI, slope angle, and fault density, while plan curvature and land use were seen to be the least important factors.     

Effects of over-consolidation history on flow instability of clean and silty sands

Acta Geotechnica - This article presents an experimental program under the constant volume condition to investigate the influence of over-consolidation on flow instability of clean and silty sand...     

Prediction of Rock Abrasivity Index (RAI) and Uniaxial Compressive Strength (UCS) of Granite Building Stones Using Nondestructive Tests

Geotechnical and Geological Engineering - Rock abrasivity index (RAI) and uniaxial compressive strength (UCS) are two key parameters for assessing abrasivity and durability of building stones,...     

Mineral-Potential Mapping for MVT Deposits with Limited Data Sets Using Landsat Data and Geological Evidence in the Borden Basin, Northern Baffin Island, Nunavut, Canada

Posterior probabilities of occurrence for Zn-Pb Mississippi Valley Type (MVT) mineralization were calculated based on evidence maps derived from regional geology, Landsat-TM, RADARSAT-1, a digital elevation model and aeromagnetic data sets in the Borden Basin of northern Baffin Island, Canada. The vector representation of geological contacts and fault traces were refined according to their characteristics identified in Landsat-TM, RADARSAT-1, DEM, slope, aspect, and shaded relief data layers. Within the study area, there is an association between the occurrence of MVT mineralization and proximity to the contact of platformal carbonates and shale units of the adjacent geological formation. A spatial association also tends to exist between mineralization and proximity to E-W and NW-SE trending faults. The relationships of known MVT occurrences with the geological features were investigated by spatial statistical techniques to generate evidence maps. Supervised classification and filtering were applied to Landsat-TM data to divide the Society Cliffs Formation into major stratigraphic subunits. Because iron oxides have been observed at some of the MVT occurrences within the Borden Basin, Landsat-TM data band ratio (3/1) was calculated to highlight the potential presence of iron-oxides as another evidence map. Processed Landsat-TM data and other derived geological evidence maps provided useful indicators for identifying areas of potential MVT mineralization. Weights of evidence and logistic regression were used independently to integrate and generate posterior probability maps showing areas of potential mineralization based on all derived evidence maps. Results indicate that in spite of the lack of important data sets such as stream or lake sediment geochemistry, Landsat-TM data and regional geological data can be useful for MVT mineral-potential mapping.