CHDS: conflict handling in direct sampling for stochastic simulation of spatial variables

Stochastic Environmental Research and Risk Assessment - In recent years, multiple-point geostatistical (MPS) approaches have gained significant popularity for modeling subsurface heterogeneity in...     

A GIS-based comparative study of hybrid fuzzy-gene expression programming and hybrid fuzzy-artificial neural network for land subsidence susceptibility modeling

Stochastic Environmental Research and Risk Assessment - Land subsidence is a complicated hazard that artificial intelligence models can model it without approximation and simplification. In this...     

Optimal Location of Wells for Storage and Recovery of Surplus Desalinated Water in Coastal Aquifers

Storage of water in aquifers using injection wells is an efficient way for utilizing excess desalinated water in arid regions. In this investigation we estimate the benefits of optimally recharging seasonal surplus desalinated water into a strategic coastal aquifer already benefitting from natural recharge of flash-floods water by a recharge dam. Since, usually the buyers of desalinated water commit to purchase surplus desalinated water under take-or-pay contracts, any attempt in utilizing the paid water is beneficial. Coastal cities are observing an increased urbanization leaving limited space for aquifer recharge infrastructure. In order to determine the optimal location of wells and maximize the use of surplus desalinated water available in winter period, a decision tool combining a numerical groundwater flow simulation model (MODFLOW) with an optimization model is developed. The results of this study show that increasing the number of wells from the existing 45 wells to 173 would allow storing 31.4 million cubic meter per year of excess desalinated water into the aquifer that can be used during later during summer months. The net benefit would reach US$55 million/year while the cost of drilling the new wells is US$5.11 million.     

Socio-hydrological framework for investigating farmers’ activities affecting the shrinkage of Urmia Lake; hybrid data mining and agent-based modelling

ABSTRACT

Developing a general framework to capture the complexities associated with the non-linear and adaptive nature of farmers facing water resources scarcity is a challenging problem. This paper integrates agent-based modelling (ABM) and a data mining method to develop a hybrid socio-hydrological framework to provide future insights for policy-makers. The data associated with the farmers’ main characteristics were collected through field surveys and interviews. Afterwards, the association rule was employed to discover the main patterns representing the farmers’ agricultural decisions. The discovered patterns were then used as the behavioural rules in ABM to simulate the agricultural activities. The proposed framework has been was applied to explore the interactions between agricultural activities and the main river feeding the Urmia-Lake, Iran. The outcomes indicate that farmers’ acquisitive traits and belongings have significant impacts on their socio-hydrological interactions. The reported values of the efficiency criteria may support the satisfactory performance of the proposed framework.     

Seismic performance and damage incurred by monolithic concrete self-centering rocking walls under the effect of axial stress ratio

Self-centering rocking walls offer the possibility of minimizing repair costs and downtimes, and also nullify the residual drift after seismic events, thanks to their self-centering properties. In this paper, the effect of axial stress ratio on the behavior of monolithic self-centering rocking walls is investigated by utilizing a developed finite element model. To verify the validity of the finite element model, results and observed damage in the model are compared with those of a full-scale wall test. The axial stress ratio is varied from 0.024 to 0.30 while keeping the other structural parameters constant. For qualitative damage evaluation, the observed damage in the model compared with expected damage states of desired performance levels. In order to evaluate the incurred damage quantitatively, the amount of crushing and damage in the wall is calculated by utilizing several ratios (crushing ratio and damage ratio). Furthermore, seismic response factors (i.e., μ, R and C_d) are calculated for different axial stress ratio values. The obtained results showed that, in order to satisfy the requirements of desired performance levels, the maximum axial stress ratio should be approximately within the range of 0.10–0.15. In addition, the maximum overall damage ratio and crushing ratio are suggested to be less than 5%. For axial stress ratio higher than 0.15, the flag-shaped pattern of hysteresis curves completely disappeared and the variation of displacement ductility is less sensitive to axial stress ratio. Considering the maximum axial stress ratio limited to 0.150, values of 4 and 3.5 are conservatively proposed as a period-independent response modification factor and displacement modification factor of the investigated structural wall, respectively.     

Three-parameter Radon transform in layered transversely isotropic media

Radon transform is a powerful tool with many applications in different stages of seismic data processing, because of its capability to focus seismic events in the transform domain. Three-parameter Radon transform can optimally focus and separate different seismic events, if its basis functions accurately match the events. In anisotropic media, the conventional hyperbolic or shifted hyperbolic basis functions lose their accuracy and cannot preserve data fidelity, especially at large offsets. To address this issue, we propose an accurate traveltime approximation for transversely isotropic media with vertical symmetry axis, and derive two versions of Radon basis functions, time-variant and time-invariant. A time-variant basis function can be used in time domain Radon transform algorithms while a time-invariant version can be used in, generally more efficient, frequency domain algorithms. Comparing the time-variant and time-invariant Radon transform by the proposed basis functions, the time-invariant version can better focus different seismic events; it is also more accurate, especially in presence of vertical heterogeneity. However, the proposed time-invariant basis functions are suitable for a specific type of layered anisotropic media, known as factorized media. We test the proposed methods and illustrate successful applications of them for trace interpolation and coherent noise attenuation.     

The Effect of Hydraulic Gradient and Pattern of Conduit Systems on Tracing Tests: Bench-Scale Modeling

Tracer breakthrough curves provide valuable information about the traced media, especially in inherently heterogeneous karst aquifers. In order to study the effect of variations in hydraulic gradient and conduit systems on breakthrough curves, a bench scale karst model was constructed. The bench scale karst model contains both matrix and a conduit. Eight tracing tests were conducted under a wide range of hydraulic gradients from 1 to greater than 5 for branchwork and network-conduit systems. Sampling points at varying distances from the injection point were utilized. Results demonstrate that mean tracer velocities, tracer mass recovery and linear rising slope of the breakthrough curves were directly controlled by hydraulic gradient. As hydraulic gradient increased, both one half the time for peak concentration and one fifth the time for peak concentration decreased. The results demonstrate the variations in one half the time for peak concentration and one fifth the time for peak concentration of the descending limb for different sampling points under differing hydraulic gradients are mainly controlled by the interactions of advection with dispersion. The results are discussed from three perspectives: different conduit systems, different hydraulic-gradient conditions, and different sampling points. The research confirmed the undeniable role of hydrogeological setting (i.e., hydraulic gradient and conduit system) on the shape of the breakthrough curve. The extracted parameters (mobile-fluid velocity, tracer-mass recovery, linear rising limb, one half the time for peak concentration, and one fifth the time for peak concentration) allow for differentiating hydrogeological settings and enhance interpretations the tracing tests in karst aquifers.     

Polynomial chaos expansion and response surface method for nonlinear modelling of reference evapotranspiration

The feasibility of polynomial chaos expansion (PCE) and response surface method (RSM) models is investigated for modelling reference evapotranspiration (ET₀). The modelling results of the proposed models are validated against the M5 model tree and multi-layer perceptron neural network (MLPNN) methods. Two meteorological stations, Isparta and Antalya, in the Mediterranean region of Turkey, are inspected. Various input combinations of daily air temperature, solar radiation, wind speed and relative humidity are constructed as input attributes for the ET₀. Generally, the modelling accuracy is increased by increasing the number of inputs. Including wind speed in the model inputs considerably increases their accuracy in modelling ET₀. Mean absolute error (MAE), root mean square error (RMSE), agreement index (d) and Nash-Sutcliffe efficiency (NSE) are used as comparison criteria. The PCE is the most accurate model in estimating daily ET₀, giving the lowest MAE (0.036 and 0.037 mm) and RMSE (0.047 and 0.050 mm) and the highest d (0.9998 and 0.9999) and NSE (0.9992 and 0.9996) with the four-input PCE models for Isparta and Antalya, respectively.     

Effects of rainfall patterns on runoff and rainfall-induced erosion

Rainfall-induced erosion involves the detachment of soil particles by raindrop impact and their transport by the combined action of the shallow surface runoff and raindrop impact.Although temporal variation in rainfall intensity(pattern)during natural rainstorms is a common phenomenon,the available information is inadequate to understand its effects on runoff and rainfall-induced erosion processes.To address this issue,four simulated rainfall patterns(constant,increasing,decreasing,and increasing-decreasing)with the same total kinetic energy were designed.Two soil types(sandy and sandy loam)were subjected to simulated rainfall using 15 cm×30 cm long detachment trays under infiltration conditions.For each simulation,runoff and sediment concentration were sampled at regular intervals.No obvious difference was observed in runoff across the two soil types,but there were significant differences in soil losses among the different rainfall patterns and stages.For varying-intensity rainfall patterns,the dominant sediment transport mechanism was not only influenced by raindrop detachment but also was affected by raindrop-induced shallow flow transport.Moreover,the efficiency of equations that predict the interrill erosion rate increased when the integrated raindrop impact and surface runoff rate were applied.Although the processes of interrill erosion are complex,the findings in this study may provide useful insight for developing models that predict the effects of rainfall pattern on runoff and erosion.     

Wave breaking and bubble formation associate energy dissipation and wave setup

Ocean Dynamics - In this work, theoretical development of void fraction effect over wave energy dissipation and wave setup in a surf zone is made. The formulation incorporates simple energy...