Tschermak fractionation in calc-alkaline magmas: the Eocene Sabzevar volcanism (NE Iran)

Calc-alkaline arc magmatism at convergent plate margins is volumetrically dominated by metaluminous andesites. Many studies highlighted the importance of differentiation via fractionation processes of arc magmas, but only in the last decades, it has been demonstrated that not all rock-forming minerals may affect the evolution of calc-alkaline suites. In particular, a major role exerted by Al-rich hornblende amphibole as fractionating mineral phase has been documented in many volcanic arc settings. The aim of this work is to understand the role of the Tschermak molecule (CaAlAlSiO₆) hosted in the hornblende and plagioclase fractionation assemblage in driving magma differentiation in calc-alkaline magmatic suites. We explore this issue by applying replenishment–fractional crystallization (RFC) and rare earth element–Rayleigh fractional crystallization (REE-FC) modeling to the Sabzevar Eocene (ca. 45–47 Ma) calc-alkaline volcanism of NE Central Iran, where hornblende-controlled fractionation has been demonstrated. Major element mass balance modeling indicates RFC dominated by a fractionating assemblage made of Hbl_52.0–52.5 + Pl_44.1–44.2 + Ttn_3.3–3.9 (phases are expressed on total crystallized assemblage). REE-FC modeling shows, instead, a lower degree of fractionation with respect to RFC models that is interpreted as due to hornblende and plagioclase resorption by the residual melt. Calculations demonstrate that fractionation of the Tschermak molecule can readily produce dacite and rhyolite magmas starting from a calc-alkaline andesite source (FC = ca. 30 %). In particular, the Tschermak molecule controls both the heavy rare earth elements (HREE) and light rare earth element (LREE) budgets in calc-alkaline differentiation trends.     

Normalized residual displacement spectra for post-mainshock assessment of structures subjected to aftershocks

Residual displacement, as a significant measure of structural inelasticity, is effectively used in post-earthquake seismic assessment of structures. This demand can be considered for seismic evaluation of structures under multiple earthquakes. This study introduces a simple and novel index to predict the residual displacement of mainshock-damaged structures against subsequent aftershock. The proposed index is defined as a ratio between residual displacement of damaged structures against aftershock and peak inelastic displacement of intact structures under mainshock. In this study, constantstrength spectra based on the index are developed considering the effects of important structural characteristics and also significant seismic parameters. Moreover, analytical equations are presented to predict the proposed index for bi-linear single-degree-of-freedom(SDOF) systems considering both the effects of positive and negative polarities of aftershock. Furthermore, an equation is suggested to estimate the peak inelastic displacement of intact systems under mainshock, which is required to compute the index.     

An alternative method for density variation modeling of the crust based on 3-D gravity inversion

In this paper we are proposing an alternative method for determination of density variations of the crust from constrained inversion of the terrestrial gravity data. The main features of the method can be summarized as follows: (i) Constructing a band-pass filter to remove the long and short wavelength signals from the terrestrial gravity data. (ii) Using an iterative method for stabilization and solution of the inverse problem. The mentioned regularization method is first validated by simulated gravity data and next the methodology is used for development of a new regional density variation model of the crust in three layers based on real gravity data in geographical area of Iran. Application of the band-pass filter to the latter data resulted the residual gravitation variations in the range of − 300 to 50 (mGal) which next based on the iterative method resulted following ranges for residual densities: −120 to 40 (kg/m³) in first layer, −40 to 40 (kg/m³) in second layer, and − 40 to 40 (kg/m³) in third layer.     

Flood hazard mapping using synthesis hydraulic and geomorphic properties at watershed scale

The occurrence of natural phenomena such as floods has caused serious consequences for human societies. The simulation of flood hazard maps and its depth in a river is one of the most complex processes in hydrology. In fact both geomorphological and hydraulic procedures for deriving the flood hazard maps and depth are imperfect at watershed scale. In this study, a combination of both procedures, using a probabilistic approach is used. Flood inundation maps for 2-, 10-, 25-,50- and 100-return period floods using flood routine within HEC-RAS in combination of Arc-GIS and topographic wetness index (TWI) map were produced. TWI threshold was identified using a maximum likelihood method in order to produce flood prone areas and calibrated over the reach of Zirab City. The correlation between TWI threshold and the flood depth was carried out and simple linear regression developed for various return periods. The resulting regression model is used in order to create flood hazard maps with various return periods at watershed scale.     

Inversion of surface gravity data for 3-D density modeling of geologic structures using total variation regularization

We develop an inversion procedure using the total variation (TV) regularization method as a stabilizing function to invert surface gravity data to retrieve 3-D density models of geologic structures with sharp boundaries. The developed inversion procedure combines several effective algorithms to solve the TV regularized problem. First, a matrix form of the gradient vector is designed using the Kronecker product to numerically approximate the 3-D TV function. The piecewise polynomial truncated singular value decomposition (PP-TSVD) algorithm is then used to solve the TV regularized inverse problem. To obtain a density model with depth resolution, we use a sensitivity-based depth weighting function. Finally, we apply the Genetic Algorithm (GA) to select the best combination of the PP-TSVD algorithm and the depth weighting function parameters. 3-D simulations conducted with synthetic data show that this approach produces sub-surface images in which the structures are well separated in terms of sharp boundaries, without the need of a priori detailed density model. The method applied to a real dataset from a micro-gravimetry survey of Gotvand Dam, southwestern Iran, clearly delineates subsurface cavities starting from a depth of 40 m within the area of the dam reservoir.     

Experimental study on local scour around single oblong pier fitted with a collar in a 180 degree flume bend

In this work, investigation on the development of local scour around an oblong pier in a 180 degree flume bend is presented. Scour hole can cause failure of the bridge especially during the river floods. In this study, the use of oblong collars for reducing the effects of local scour at a bridge pier is presented together with the time aspect of the scour development. Tests were conducted using one oblong pier in positions of 60degree under one flow conditions. The study was conducted using a physical hydraulic model operated under clear-water conditions in cohesionless bed material. In this study, the time development of the local scour around the oblong pier fitted with and without collar plates was studied. Investigated was the effect of size and elevation collar on the time development of scour and its efficacy at preventing scour at a bridge pier. The time development of the scour hole around the model pier with and without a collar installed was compared with similar studies on bridge piers. The results of the model study indicated that the maximum depth of scour is highly dependent on the experimental duration. It was observed that, as the minimum depth of scour occurs for the square collar at width of 3B placed at elevation of 0.1B below the bed and the size of a collar plate increases, the scour decreases. Measuring depth of scouring based on experimental observation, an empirical relation is developed with regression coefficient 95%.     

Metal concentrations in selected tissues and main prey species of the annulated sea snake (Hydrophis cyanocinctus) in the Hara Protected Area, northeastern coast of the Persian Gulf, Iran

This study is the first detailed ecotoxicological study of the annulated sea snake, Hydrophis cyanocinctus. Concentrations of lead, cadmium, nickel and vanadium were evaluated in muscle, liver, kidney, skin and blood of the annulated sea snake (H. cyanocinctus) and in the whole bodies of its main prey species (Periophthalmus waltoni and Boleophthalmus dussumieri) in the Hara Protected Area, the Persian Gulf. The mean concentrations of lead and vanadium were highest in the kidney, which identified the kidney as a target organ for metals in sea snakes as it is in other reptilian groups. Mean concentrations of cadmium and nickel were highest in the liver and skin, respectively. Mean cadmium concentrations were significantly higher in the liver compared to prey species, which indicated that prey items may be a source of cadmium for the annulated sea snake in the study area. Data presented here may be considered as a baseline for further ecotoxicological studies in sea snakes.     

Phenolic endocrine disrupting chemicals (EDCs) in Anzali Wetland, Iran: elevated concentrations of 4-nonylphenol, octhylphenol and bisphenol A

We have studied the distribution and value of phenolic endocrine disrupting chemicals (EDCs) in surface sediment samples taken from Anzali Wetland, Iran. These samples were collected from 22 stations during the time span of June-May 2010. In each of the sampling stations, we detected 4-nonylphenol (4-NP), octylphenol (OP), and bisphenol A (BPA) with maximal concentrations of 29, 4.3, and 7 μg g(-1) dry weight (dw), respectively. High levels of alkylphenols (APs) and BPA were also found near urban areas. Furthermore there were no significant differences between those stations in terms of the detected levels. One of the important factors in controlling the fate of these compounds in the aquatic environment appeared to be Total Organic Carbon (TOC). Hierarchical cluster analysis showed differences in the biomarker characteristics of EDCs and TOC between the stations. Our findings indicate that EDCs are ubiquitous in sediments from northeast Wetlands of Iran, contaminating the aquatic habitats in this area.     

Biomonitoring of Ni and V Contamination Using Oysters (Saccostrea cucullata) at Lengeh Port,Persian Gulf,Iran

Vanadium and nickel are two important indicators of oil pollution. Lengthy exposure to these elements causes serious harmful effects in human health, different harsh allergies being examples. The accumulation of two heavy metals (Ni and V) in sediment and the soft and hard tissues of Saccostrea cucullata were analyzed at three sampling sites along the coast of Lengeh Port, Persian Gulf. Results indicated at all the sampling sites; the Ni levels in soft tissues (STs) were higher than in the shells (SHs) and sediments, whereas the V levels were higher in the sediments. In addition, meaningful relationship (r = 0.65; p < 0.05) was observed across Ni levels in ST of S. cucullata and sediment, while for V concentrations a strong relationship (r = 0.83; p < 0.01) was found in SH of S. cucullata and sediment. This indicates that ST and SH of oyster can be considered as a biomonitoring agent for Ni and V levels, in coastal waters, respectively. The exposure of the consumer is compared directly to minimal risk level and provisional maximum tolerable daily intake. Result indicated that levels of Ni and V were within the safety limits for human consumption.