Revised age of proximal deposits in the Zagros foreland basin and implications for Cenozoic evolution of the High Zagros

The regionally extensive, coarse-grained Bakhtiyari Formation represents the youngest synorogenic fill in the Zagros foreland basin of Iran. The Bakhtiyari is present throughout the Zagros fold-thrust belt and consists of conglomerate with subordinate sandstone and marl. The formation is up to 3000 m thick and was deposited in foredeep and wedge-top depocenters flanked by fold-thrust structures. Although the Bakhtiyari concordantly overlies Miocene deposits in foreland regions, an angular unconformity above tilted Paleozoic to Miocene rocks is expressed in the hinterland (High Zagros).

The Bakhtiyari Formation has been widely considered to be a regional sheet of Pliocene–Pleistocene conglomerate deposited during and after major late Miocene–Pliocene shortening. It is further believed that rapid fold growth and Bakhtiyari deposition commenced simultaneously across the fold-thrust belt, with limited migration from hinterland (NE) to foreland (SW). Thus, the Bakhtiyari is generally interpreted as an unmistakable time indicator for shortening and surface uplift across the Zagros. However, new structural and stratigraphic data show that the most-proximal Bakhtiyari exposures, in the High Zagros south of Shahr-kord, were deposited during the early Miocene and probably Oligocene. In this locality, a coarse-grained Bakhtiyari succession several hundred meters thick contains gray marl, limestone, and sandstone with diagnostic marine pelecypod, gastropod, coral, and coralline algae fossils. Foraminiferal and palynological species indicate deposition during early Miocene time. However, the lower Miocene marine interval lies in angular unconformity above ~ 150 m of Bakhtiyari conglomerate that, in turn, unconformably caps an Oligocene marine sequence. These relationships attest to syndepositional deformation and suggest that the oldest Bakhtiyari conglomerate could be Oligocene in age.

The new age information constrains the timing of initial foreland-basin development and proximal Bakhtiyari deposition in the Zagros hinterland. These findings reveal that structural evolution of the High Zagros was underway by early Miocene and probably Oligocene time, earlier than commonly envisioned. The age of the Bakhtiyari Formation in the High Zagros contrasts significantly with the Pliocene–Quaternary Bakhtiyari deposits near the modern deformation front, suggesting a long-term (> 20 Myr) advance of deformation toward the foreland.     

Detrital zircon provenance of Neoproterozoic to Cenozoic deposits in Iran: Implications for chronostratigraphy and collisional tectonics

Ion-microprobe U–Pb analyses of 589 detrital zircon grains from 14 sandstones of the Alborz mountains, Zagros mountains, and central Iranian plateau provide an initial framework for understanding the Neoproterozoic to Cenozoic provenance history of Iran. The results place improved chronological constraints on the age of earliest sediment accumulation during Neoproterozoic–Cambrian time, the timing of the Mesozoic Iran–Eurasia collision and Cenozoic Arabia–Eurasia collision, and the contribution of various sediment sources of Gondwanan and Eurasian affinity during opening and closure of the Paleotethys and Neotethys oceans. The zircon age populations suggest that deposition of the extensive ~ 1 km-thick clastic sequence at the base of the cover succession commenced in latest Neoproterozoic and terminated by Middle Cambrian time. Comparison of the geochronological data with detrital zircon ages for northern Gondwana reveals that sediment principally derived from the East African orogen covered a vast region encompassing northern Africa and the Middle East. Although most previous studies propose a simple passive-margin setting for Paleozoic Iran, detrital zircon age spectra indicate Late Devonian–Early Permian and Cambrian–Ordovician magmatism. These data suggest that Iran was affiliated with Eurasian magmatic arcs or that rift-related magmatic activity during opening of Paleotethys and Neotethys was more pronounced than thought along the northern Gondwanan passive-margin. For a Triassic–Jurassic clastic overlap assemblage (Shemshak Formation) in the Alborz mountains, U–Pb zircon ages provide chronostratigraphic age control requiring collision of Iran with Eurasia by late Carnian–early Norian time (220–210 Ma). Finally, Cenozoic strata yield abundant zircons of Eocene age, consistent with derivation from arc magmatic rocks related to late-stage subduction and/or breakoff of the Neotethys slab. Together with the timing of foreland basin sedimentation in the Zagros, these detrital zircon ages help bracket the onset of the Arabia–Eurasia collision in Iran between middle Eocene and late Oligocene time.     

Calcium isotope (δCa) fractionation along hydrothermal pathways, Logatchev field (Mid-Atlantic Ridge, 14°45′N)

We investigate the Logatchev Hydrothermal Field at the Mid-Atlantic Ridge, 14°45′N to constrain the calcium isotope hydrothermal flux into the ocean. During the transformation of seawater to a hydrothermal solution, the Ca concentration of pristine seawater ([Ca]_SW) increases from about 10 mM to about 32 mM in the hydrothermal fluid endmember ([Ca]_HydEnd) and thereby adopts a δ^44/40Ca_HydEnd of −0.95 ± 0.07‰ relative to seawater (SW) and a ⁸⁷Sr/⁸⁶Sr isotope ratio of 0.7034(4). We demonstrate that δ^44/40Ca_HydEnd is higher than that of the bedrock at the Logatchev field. From mass balance calculations, we deduce a δ^44/40Ca of −1.17 ± 0.04‰ (SW) for the host-rocks in the reaction zone and −1.45 ± 0.05‰ (SW) for the isotopic composition of the entire hydrothermal cell of the Logatchev field. The values are isotopically lighter than the currently assumed δ^44/40Ca for Bulk Earth of −0.92 ± 0.18‰ (SW) [Skulan J., DePaolo D. J. and Owens T. L. (1997) Biological control of calcium isotopic abundances in the global calcium cycle. Geochim. Cosmochim. Acta61,(12) 2505-2510] and challenge previous assumptions of no Ca isotope fractionation between hydrothermal fluid and the oceanic crust [Zhu P. and Macdougall J. D. (1998) Calcium isotopes in the marine environment and the oceanic calcium cycle. Geochim. Cosmochim. Acta62,(10) 1691-1698; Schmitt A. -D., Chabeaux F. and Stille P. (2003) The calcium riverine and hydrothermal isotopic fluxes and the oceanic calcium mass balance. Earth Planet. Sci. Lett. 6731, 1-16]. Here we propose that Ca isotope fractionation along the fluid flow pathway of the Logatchev field occurs during the precipitation of anhydrite. Two anhydrite samples from the Logatchev Hydrothermal Field show an average fractionation of about Δ^44/40Ca = −0.5‰ relative to their assumed parental solutions. Ca isotope ratios in aragonites from carbonate veins from ODP drill cores indicate aragonite precipitation directly from seawater at low temperatures with an average δ^44/40Ca of −1.54 ± 0.08‰ (SW). The relatively large fractionation between the aragonite precipitates and seawater in combination with their frequent abundance in weathered mafic and ultramafic rocks suggest a reconsideration of the marine Ca isotope budget, in particular with regard to ocean crust alteration.     

Chemical characteristics of groundwater in parts of mountainous region, Alvand, Hamadan, Iran

Eighty-seven groundwater samples have been collected from a mountainous region (Alvand, Iran) for hydrochemical investigations to understand the sources of dissolved ions and assess the chemical quality of the groundwater. Most water quality parameters are within World Health Organization acceptable limits set for drinking water. The least mineralized water is found closest to the main recharge zones and the salinity of water increased towards the north of the basin. The most prevalent water type is Ca–HCO₃ followed by water types Ca–NO₃, Ca–Cl, Ca–SO₄ and Mg–HCO₃. The Ca–NO₃ water type is associated with high nitrate pollution. Agricultural and industrial activities were associated with elevated level of NO₃⁻. Mineral dissolution/weathering of evaporites dominates the major element hydrochemistry of the area. Chemical properties of groundwater in Alvand region are controlled both by natural geochemical processes and anthropogenic activities.     

Integrating cellular automata,artificial neural network,and fuzzy set theory to simulate threatened orchards: application to Maragheh,Iran

Urbanization processes challenge the growth of orchards in many cities in Iran. In Maragheh, orchards are crucial ecological, economical, and tourist sources. To explore orchards threatened by urban expansion, this study first aims to develop a new model by coupling cellular automata (CA) and artificial neural network with fuzzy set theory (CA–ANN–Fuzzy). While fuzzy set theory captures the uncertainty associated with transition rules, the ANN considers spatial and temporal nonlinearities of the driving forces underlying the urban growth processes. Second, the CA–ANN–Fuzzy model is compared with two existing approaches, namely a basic CA and a CA coupled with an ANN (CA–ANN). Third, we quantify the amount of orchard loss during the last three decades as well as for the upcoming years up to 2025. Results show that CA–ANN–Fuzzy with 83% kappa coefficient performs significantly better than conventional CA (with 51% kappa coefficient) and CA–ANN (with 79% kappa coefficient) models in simulating orchard loss. The historical data shows a considerable loss of 26% during the last three decades, while the CA–ANN–Fuzzy simulation reveals a considerable future loss of 7% of Maragheh’s orchards in 2025 due to urbanization. These areas require special attention and must be protected by the local government and decision-makers.     

Prediction of water inflow to mechanized tunnels during tunnel-boring-machine advance using numerical simulation

An accurate estimate of the groundwater inflow to a tunnel is one of the most challenging but essential tasks in tunnel design and construction. Most of the numerical or analytical solutions that have been developed ignore tunnel seepage conditions, material properties and hydraulic-head changes along the tunnel route during the excavation process, leading to inaccurate prediction of inflow rates. A method is introduced that uses MODFLOW code of GMS software to predict inflow rate as the tunnel boring machine (TBM) gradually advances. In this method, the tunnel boundary condition is conceptualized and defined using Drain package, which is simulated by dividing the drilling process into a series of successive intervals based on the tunnel excavation rates. In addition, the drain elevations are specified as the respective tunnel elevations, and the conductance parameters are assigned to intervals, depending on the TBM type and the tunnel seepage condition. The Qomroud water conveyance tunnel, located in Lorestan province of Iran, is 36 km in length. Since the Qomroud tunnel involved groundwater inrush during excavating, it is considered as a good case study to evaluate the presented method. The groundwater inflow to this tunnel during the TBM advance is simulated using the proposed method and the predicted rates are compared with observed rates. The results show that the presented method can satisfactorily predict the inflow rates as the TBM advances.