Regional distribution pattern of groundwater heavy metals resulting from agricultural activities

Contaminations of groundwater by heavy metals due to agricultural activities are growing recently. The objective of this study was to evaluate and map regional patterns of heavy metals (Cd, Zn and Cu) in groundwater on a plain with high agricultural activities. The study was conducted to investigate the concentration of heavy metals and distribution in groundwater in regions of Shush Danial and Andimeshk aquifers in the southern part of Iran. Presently, groundwater is the only appropriate and widely used source of drinking water for rural and urban communities in this region. The region covers an area of 1,100 km² between the Dez and Karkhe rivers, which lead to the Persian Gulf. For this study, the region was divided into four sub-regions A, B, C and D. Additionally, 168 groundwater samples were collected from 42 water wells during the earlier months of 2004. The flame atomic absorption spectrometry (AAS-Flame) was used to measure the concentration of heavy metals in water samples and the Surfer software was used for determination of the contour map of metal distribution. The results demonstrated that in all of the samples, Cd and Zn concentrations were below the EPA MCLG and EPA secondary standard, respectively. However, the Cu contents of 4.8 % of all samples were higher than EPA MCL. It is also indicated that the concentrations of metals were more pronounced at the southern part of the studied region than at the others. The analysis of fertilizers applied for agricultural activities at this region also indicated that a great majority of the above-mentioned heavy metals were discharged into the environment. Absence of confining layers, proximity to land surface, excess agricultural activities in the southern part and groundwater flow direction that is generally from the north to the southern parts in this area make the southern region of the Shush plain especially vulnerable to pollution by heavy metals than by other contaminants.     

Shape optimization of two-dimensional cavitators in supercavitating flows, using NSGA II algorithm

The reduction of energy consumption of high speed submersible bodies is an important challenge in hydrodynamic researches. In this paper, shape optimization of two-dimensional cavitators in supercavitating flows is studied. A two dimensional supercavitation potential flow passes a symmetric two dimensional cavitator, which is placed perpendicular to the flow in a channel of infinite length and immediately a cavity is formed behind the cavitator. This is because of the generation of a gas or vapor cavity between the body and the surrounding liquid due to the change in a high speed flow direction passing the cavitator. Drag force acting on this supercavitating body dictates the thrust requirements for the propulsion system, to maintain a required cavity at the operating speed. Therefore, any reduction in the drag force, by modifying the shape of the cavitator, will lead to decrease this force. This study concentrates on the optimization of two dimensional cavitators in order to decrease drag coefficient for a specified after body length and velocity in a potential flow. To achieve this goal a multi-objective optimization problem is defined to optimize cavitator shapes in supercavitating flow. The so-called NSGA II (Non-dominated Sorting Genetic Algorithm) algorithm is used as an optimization method. Design parameters and constraints are obtained according to supercavitating flow characteristics and cavitator modeling and objective functions are generated using Linear Regression Method. The obtained results are compared with other classic optimization methods, like the weighted sum method, for validation.     

Eocene continental dyke swarm from Central Iran (Khur area)

The Eocene dyke swarm with east-west general trend intrudes the Cretaceous sedimentary rocks in ～25 km north of the Khur city (Central Iran). Some of the studied dykes can be followed for over 7 km, but the majority of exposures in the area are less than 5 km long. The dykes commonly exhibit a chilled contact with the wall rocks. These dykes are trachybasalt and basalt in composition. The trachybasalt dykes are much more abundant. The basaltic dykes cross cut the trachybasalt dykes in some locations, indicating that trachybasalt dykes are older than the basaltic ones. Primary igneous minerals of the basaltic dykes are olivine (chrysolite), clinopyroxene (diopside, augite), plagioclase (labradorite), sanidine, magnetite, orthopyroxene (enstatite), spinel and phlogopite, and secondary minerals are zeolite (natrolite and mesolite), chlorite (diabantite), calcite and serpentine. The trachybasalt dykes are composed of clinopyroxene (diopside), plagioclase (labradorite), sanidine, mica (biotite and phlogopite), amphibole (magnesio-hastingsite) and magnetite as primary minerals, and chlorite and calcite as secondary ones. Whole rocks geochemical data of the studied dykes indicate their basic and calc-alkaline nature and suggest that these two set of dykes were derived from the same parental magma. The chondrite-normalized REE patterns and the primitive mantle-normalized multi-elemental diagram of the Khur dykes show enrichment of light rare earth elements (LREE) relative to heavy rare earth elements (HREE), and negative anomalies of high field strength elements (HFSE) (e.g. Ti, Nb and Ta). These rocks show enrichment of the large ion lithophile elements (LILE) (e.g. Cs, Ba, Th and U) and depletion of the HREE and Y relative to MREE, Zr and Hf. In the chondrite-normalized REE diagram, the basalts show elevated REE abundances relative to the trachybasalt samples. Geochemical analyses of the studied samples suggest a spinel lherzolite from the mantle as the source rock and confirm the role of subduction in their generation. The chemical characteristics of the Khur dykes resemble those of continental arc rocks, and they were possibly formed by subduction of the Central-East Iranian microcontinent (CEIM) confining oceanic crust and decompression melting of a lithospheric subcontinental mantle spinel lherzolite enriched by subduction.     

Land assessment for flood spreading site selection using geospatial information system

Irregular patterns of precipitations from temporal as well as spatial perspectives not only cause destructions but also waste surface water resources. Hence, controlling surface water and leading the flood to underground stores improve the efficiency of water usage. Selecting appropriate sites for optimal use of water floods is one of the most important factors in recharging underground water tables in dry lands where the agricultural and rangelands are vulnerable. Traditional methods of site selections are, however, time consuming and error prone. This paper attempts to; analyze existing schemes of site selection; introduces an appropriate method of locating flood-spreading sites using Geospatial Information System; implements the strategy in a case study; and scientifically assesses its results. The study area of this research is Samal sub-basin covering 31571.7 ha of Ahrom basin in Boushehr province. In the present research, factors pertain to earth sciences (quaternary units, slope and landuse) and hydrology (runoff infiltration rate and aquifers’ depth) are considered. Information layers are weighted, classified and integrated through several models such as boolean logics, index overlay and fuzzy logics. The results are then checked against the existing sites to estimate their accuracy. The results of this research demonstrated that fuzzy logic operators including gamma=0.1, gamma=0.2 and products of fuzzy logics yield the best when compared to control fields and therefore, the models are introduced as the most suitable site selection strategies for flood spreading.     

Wildfire impacts on education and healthcare: Paradise,California, after the Camp Fire

Natural Hazards - The 2018 Camp Fire caused significant damages to the education and healthcare systems in the town of Paradise, CA. This paper presents the findings of a qualitative case study...     

Simulation of drilling‐induced compaction bands using discrete element method

Rock failure is observed around boreholes often with certain types of failure zones, which are called breakouts. Laboratory‐scale drilling tests in some high‐porosity quartz‐rich sandstone have shown breakouts in the form of narrow localized compacted zones in the minimum horizontal stress direction. They are called fracture‐like breakouts. Such compaction bands may affect hydrocarbon extraction by forming barriers that inhibit fluid flow and may also be a source of sand production. This paper presents the results of numerical simulations of borehole breakouts using 3D discrete element method to investigate the mechanism of the fracture‐like breakouts and to identify the role of far‐field stresses on the breakout dimensions. The numerical tool was first verified against analytical solutions. It was then utilized to investigate the failure mechanism and breakout geometry for drilled cubic rock samples of Castlegate sandstone subjected to different pre‐existing far‐field stresses. Results show that failure occurs in the zones of the highest concentration of tangential stress around the borehole. It is concluded that fracture‐like breakout develops as a result of a nondilatant failure mechanism consisting of localized grain debonding and repacking and grain crushing that lead to the formation of a compaction band in the minimum horizontal stress direction. In addition, it is found that the length of fracture‐like breakouts depends on both the mean stress and stress anisotropy. However, the width of the breakout is not significantly changed by the far‐field stresses. Copyright © 2013 John Wiley & Sons, Ltd.     

Adsorption of diethyl and dibutyl phthalates onto activated carbon produced from <Emphasis Type="Italic">Albizia julibrissin pods:</Emphasis> kinetics and isotherms

The adsorption behavior study of diethyl and dibutyl phthalates was investigated onto a new activated carbon prepared from an abundant biomass “Albizia julibrissin pods,” treated chemically by H₃PO₄. A series of experiments were conducted in a batch system to estimate the effect of operating conditions such as the adsorbent nature, the dose of adsorbent, the contact time, the initial concentration and the temperature on the adsorption efficiency. The optimum operating conditions were found to be 0.1 and 0.05 g of adsorbent for diethyl and dibutyl phthalates, respectively, at 30 min equilibrium time, 150 mg g^?1 and 293 K. The adsorption isotherms for both phthalates were fit at different temperatures using the nonlinear regression of Langmuir, Freundlich, Dubinin–Radushkevich and Redlich–Peterson. The pseudo-first order, pseudo-second order by nonlinear regression and intraparticle diffusion models were used to describe the adsorption kinetic. The results show that the intraparticle diffusion model is not the limiting step governing the adsorption mechanism. The structural and textural characteristics of adsorbent surface were investigated. FTIR analysis of unloaded and phthalates-loaded adsorbent revealed that the aliphatic groups attached to phthalate esters are involved in adsorption mechanism.     

Constitutive model for cyclic behaviour of cohesionless sands

This paper presents a constitutive model for describing the stress-strain response of sands under cyclic loading. The model, formulated using the critical state theory within the bounding surface plasticity framework, is an upgraded version of an existing model developed for monotonic behaviour of cohesionless sands. With modification of the hardening law, plastic volumetric strain increment and unloading plastic modulus, the original model was modified to simulate cyclic loading. The proposed model was validated against triaxial cyclic loading tests for Fuji River sand, Toyoura sand and Nigata sand. Comparison between the measured and predicted results suggests that the proposed modified model can capture the main features of cohesionless sands under drained and undrained cyclic loading.     

Recent developments of the Middle East catalog

This article summarizes a recent study in the framework of the Global Earth model (GEM) and the Earthquake Model of the Middle East (EMME) project to establish the new catalog of seismicity for the Middle East, using all historical (pre-1900), early and modern instrumental events up to 2006. According to different seismicity, which depends on geophysical, geological, tectonic, and seismicity data, this region is subdivided to nine subregions, consisting of Alborz–Azerbaijan, Afghanistan–Pakistan, Saudi Arabia, Caucasus, Central Iran, Kopeh–Dagh, Makran, Zagros, and Turkey (Eastern Anatolia; after 30° E). After omitting the duplicate events, aftershocks, and foreshocks by using the Gruenthal method, and uniform all magnitude to Mw scale, 28,244 main events remain for the new catalog of Middle East from 1250 B.C. through 2006. The magnitude of completeness (Mc) was determined as 4.9 for five out of nine subregions, where the least values of Mc were found to be 4.2. The threshold of Mc is around 5.5, 5.0, 4.5, and 4.0, for the time after 1950, 1963, 1975, and 2000, respectively. The average of teleseismic depths in all regions is less than 15 km. Totally, majority of depth for Kopeh–Dagh and Central Iran, Zagros, and Alborz–Azerbaijan, approximately, is 15, 13, and 11 km and for Afghanistan–Pakistan, Caucasus, Makran, Turkey (after 30° E), and Saudi Arabia is about 9 km.     

Seismic evaluation of rocking structures through performance assessment and fragility analysis

Numerical studies have been conducted for low- and medium-rise rocking structures to investigate their efficiency as earthquake-resisting systems in comparison with conventional structures. Several non-linear time-history analyses have been performed to evaluate seismic performance of selected cases at desired ground shaking levels, based on key parameters such as total and flexural story drifts and residual deformations. The Far-field record set is selected as input ground motions and median peak values of key parameters are taken as best estimates of system response. In addition, in order to evaluate the probability of exceeding relevant damage states, analytical fragility curves have been developed based on the results of the incremental dynamic analysis procedure. Small exceedance probabilities and acceptable margins against collapse, together with minor associated damages in main structural members, can be considered as superior seismic performance for medium-rise rocking systems. Low-rise rocking systems could provide significant performance improvement over their conventional counterparts notwithstanding certain weaknesses in their seismic response.