Assessment of floodplain landuse and channel morphology within meandering reach of the Talar River in Iran using GIS and aerial photographs

Abstract

The main objective of this study is to investigate and monitor the landuse and morphological changes in the floodplain part of the Talar River, northern parts of Iran. In the present study, the aerial photographs have been used to produce landuse maps of the floodplain for three periods including 1968, 1994, and 2013. The quantitative analysis of the produced landuse maps showed that the floodplain has undergone substantial landuse changes. Moreover, the sediment bar and the beach area have been decreased about 97 and 90%, respectively, during the 45-year period. Substantial increases of 192 and 622% have been observed for orchards and residential areas, respectively. On the other hand, not only the forest and riparian vegetation were decreased but also the average width of river was decreased about 25.5 m. In addition, flow length of the study reach increased about 8 m. The RNCI was about ?0.7 m per year indicating sedimentation process. During the period of 1968–2013, Caspian Sea has retreated about 150 m and the delta of Talar River was changed. This study showed that morphological actions during first 26 years (1968–1994) were the stable and last 19 years had the change period, especially sedimentation (bar).     

Soil erosion and sediment yield modeling using RS and GIS techniques: a case study, Iran

Water erosion is a serious and continuous environmental problem in many parts of the world. The need to quantify the amount of erosion, sediment delivery, and sediment yield in a spatially distributed form has become essential at the watershed scale and in the implementation of conservation efforts. In this study, an effort to predict potential annual soil loss and sediment yield is conducted by using the Revised Universal Soil Loss Equation (RUSLE) model with adaptation in a geographic information system (GIS). The rainfall erosivity, soil erosivity, slope length, steepness, plant cover, and management practice and conservation support practice factors are among the basic factors that are obtained from monthly and annual rainfall data, soil map of the region, 50-m digital elevation model, remote sensing (RS) techniques (with use of Normalized Difference Vegetation Index), and GIS, respectively. The Ilam dam watershed which is located southeast part of Ilam province in western Iran is considered as study area. The study indicates that the slope length and steepness of the RUSLE model are the most effective factors controlling soil erosion in the region. The mean annual soil loss and sediment yield are also predicted. Moreover, the results indicated that 45.25%, 12.18%, 12.44%, 10.79%, and 19.34% of the study area are under minimal, low, moderate, high, and extreme actual erosion risks, respectively. Since 30.13% of the region is under high and extreme erosion risk, adoption of suitable conservation measures seems to be inevitable. So, the RUSLE model integrated with RS and GIS techniques has a great potential for producing accurate and inexpensive erosion and sediment yield risk maps in Iran.     

Land subsidence pattern controlled by old alpine basement faults in the Kashmar Valley, northeast Iran: results from InSAR and levelling

Fluid storage systems, such as oil, gas, magma or water reservoirs, are often controlled by the host rock structure and faulted terrain. In sedimentary basins, where no direct information about underlying structure is available, the pattern of ground deformation may allow us to assess the buried fault arrangement. We provide an example in the semi-arid area of Iran, in the Kashmar Valley, a region subject to land subsidence due to water overexploitation. Geodetically determined subsidence rates in the Kashmar Valley exceed 15–30 cm yr⁻¹. The pattern of surface deformation is strongly non-uniform and displays NE–SW elongated bowls of subsidence. The trend resembles old Cretaceous-to-Tertiary faults that evolved during early alpine tectonic deformation. Although these early alpine structures are considered tectonically inactive in the present day, the observed land subsidence pattern indicates significant structural control on the geometry of the aquifer basin and its deformation during reservoir drainage.     

Three-phase measurements of oil and gas trapping in sand packs

We measure the trapped saturations of oil and gas as a function of initial saturation in water-wet sand packs. We start with a water-saturated column and inject octane (oil), while water and oil are produced from the bottom. Once water production has ceased, air (gas) then enters from the top, allowing oil and gas to drain under gravity for different times. Finally water is then injected from the bottom to trap both oil and gas. The columns are sliced and the fluids analyzed using gas chromatography. We find that for high initial gas saturations more gas can be trapped in the presence of oil than in a two-phase (gas/water) system. The residual gas saturation can be over 20% compared to 14% in two-phase flow [Al Mansoori SK, Iglauer S, Pentland CH, Bijeljic B, Blunt MJ. Measurements of non-wetting phase trapping applied to carbon dioxide storage. Energy Procedia 2009;1(1):3173–80]. This is unlike previous measurements on consolidated media, where the trapped gas saturation is either similar or lower to that reached in an equivalent two-phase experiment. For lower initial gas saturation, the amount of trapping follows the initial-residual trend seen in two-phase experiments. The amount of oil trapped is insensitive to initial gas saturation or the amount of gas that is trapped, again in contrast to measurements on consolidated media. More oil is trapped than would be predicted from an equivalent two-phase (oil/water) system, although the trapped saturation is never larger than the maximum reached in two-phase flow (around 11%) [Pentland CH, Al Mansoori SK, Iglauer S, Bijeljic B, Blunt MJ. Measurement of non-wetting phase trapping in sand packs. In: SPE 115697, proceedings of the SPE annual technical conference and exhibition, Denver, Colorado, USA; 21–24 September 2008]. These initially surprising results are explained in the context of oil layer stability and the competition between snap-off and piston-like advance. In two-phase systems, displacement is principally by cooperative piston-like advance with relatively little trapping, whereas in consolidated media snap-off is generally more significant. However, oil layer collapse events during three-phase waterflooding rapidly trap the oil which acts as a barrier to direct water/gas displacement, except by snap-off, leading to enhanced gas trapping.     

Three-dimensional semi-analytical solution to groundwater flow in confined and unconfined wedge-shaped aquifers

The Laplace domain solutions have been obtained for three-dimensional groundwater flow to a well in confined and unconfined wedge-shaped aquifers. The solutions take into account partial penetration effects, instantaneous drainage or delayed yield, vertical anisotropy and the water table boundary condition. As a basis, the Laplace domain solutions for drawdown created by a point source in uniform, anisotropic confined and unconfined wedge-shaped aquifers are first derived. Then, by the principle of superposition the point source solutions are extended to the cases of partially and fully penetrating wells. Unlike the previous solution for the confined aquifer that contains improper integrals arising from the Hankel transform [Yeh HD, Chang YC. New analytical solutions for groundwater flow in wedge-shaped aquifers with various topographic boundary conditions. Adv Water Resour 2006;26:471–80], numerical evaluation of our solution is relatively easy using well known numerical Laplace inversion methods. The effects of wedge angle, pumping well location and observation point location on drawdown and the effects of partial penetration, screen location and delay index on the wedge boundary hydraulic gradient in unconfined aquifers have also been investigated. The results are presented in the form of dimensionless drawdown-time and boundary gradient-time type curves. The curves are useful for parameter identification, calculation of stream depletion rates and the assessment of water budgets in river basins.     

Nonparametric estimation of wave dispersion in high‐rise buildings by seismic interferometry

Interferometric identification and health monitoring of high‐rise buildings has been gaining increasing interest in recent years. The wave dispersion in the structure has been largely ignored in these efforts but needs to be considered to further develop these methods. In this paper, (i) the goodness of estimation of vertical wave velocity in buildings, as function of frequency, by two nonparametric interferometric techniques is examined, using realistic fixed‐base Timoshenko beam benchmark models. Such models are convenient because the variation of phase and group velocities with frequency can be derived theoretically. The models are those of the NS and EW responses of Millikan Library. One of the techniques, deconvolution interferometry, estimates the phase velocity on a frequency band from phase difference between motions at two locations in the structure, while the other one estimates it approximately at the resonant frequencies based on standing wave patterns. The paper also (ii) examines the modeling error in wave velocity profiles identified by fitting layered shear beam in broader band impulse response functions of buildings with significant bending flexibility. This error may affect inferences on the spatial distribution of damage from detected changes in such velocity profiles. Copyright © 2014 John Wiley & Sons, Ltd.     

Ground motion record simulation for structural analysis by consideration of spectral acceleration autocorrelation pattern

A novel approach is introduced to generate simulated ground motion records by considering spectral acceleration correlations at multiple periods. Most of the current reliable Ground Motion Record(GMR) simulation procedures use a seismological model including source, path and site characteristics. However, the response spectrum of simulated GMR is somewhat different when compared with the response spectrum based on recorded GMRs. More specifi cally, the correlation between the spectral values at multiple periods is a characteristic of a record which is usually different between simulated and recorded GMRs. As this correlation has a signifi cant infl uence on the structural response, it is needed to investigate the consistency of the simulated ground motions with actual records. This issue has been investigated in this study by incorporating an optimization algorithm within the Boore simulation technique. Eight seismological key parameters were optimized in order to achieve approximately the same correlation coeffi cients and spectral acceleration between two sets of real and simulated records. The results show that the acceleration response spectra of the synthetic ground motions also have good agreement with the real recorded response spectra by implementation of the proposed optimized values.     

Semi-two-dimensional numerical model for river morphological change prediction: theory and concepts

This paper presents a new numerical model for river morphological predictions. This tool predicts vertical and lateral cross-section variations for alluvial rivers, which is an important task in predicting the associated hazard zone after a flood event. The Model for the HYdraulics of SEdiments in Rivers, version 1.0 (MHYSER 1.0) is a semi-two-dimensional model using the stream tubes concept to achieve lateral variations of velocity, flow stresses, and sediment transport rates. Each stream tube has the same conveyance as the other ones. In MHYSER 1.0, the uncoupled approach is used to solve the set of conservation equations. After the backwater calculation, the river is divided into a finite number of stream tubes of equal conveyances. The sediment routing and bed adjustments calculations are accomplished separately along each stream tube taking into account lateral mass exchanges. The determination of depth and width adjustments is based on the minimum stream power theory. Moreover, MHYSER 1.0 offers two options to treat riverbank stability. The first one is based on the angle of repose. The bank slope should not be allowed to increase beyond a certain critical value supplied to MHYSER 1.0. The second one is based on the modified Bishop’s method to determine a safety factor evaluating the potential risk of a landslide along the river bank.     

Deep learning neural networks for spatially explicit prediction of flash flood probability

Flood probability maps are essential for a range of applications, including land use planning and developing mitigation strategies and early warning systems. Th...