A model for moisture and heat flow in fractured unsaturated porous media

The present study investigates propagation of a cohesive crack in non‐isothermal unsaturated porous medium under mode I conditions. Basic points of skeleton deformation, moisture, and heat transfer for unsaturated porous medium are presented. Boundary conditions on the crack surface that consist of mechanical interaction of the crack and the porous medium, water, and heat flows through the crack are taken into consideration. For spatial discretization, the extended finite element method is used. This method uses enriched shape functions in addition to ordinary shape functions for approximation of displacement, pressure, and temperature fields. The Heaviside step function and the distance function are exploited as enrichment functions for representing the crack surfaces displacement and the discontinuous vertical gradients of the pressure and temperature fields along the crack, respectively. For temporal discretization, backward finite difference scheme is applied. Problems solved from the literature show the validity of the model as well as the dependency of structural response on the material properties and loading. Copyright © 2016 John Wiley & Sons, Ltd.     

Lower Oligocene Calc-Alkaline Spessartitic Lamprophyres from Central Iran (East of Anarak Area); an Evidence from the Eastern Branch of Neotethys Subduction-Related Mantle Enrichment

Geotectonics - The Lower Oligocene Kal-e-kafi (East of Anarak, Central Iran) lamprophyres occur as stocks and dikes, which cross-cut the Eocene volcanic and Cretaceous sedimentary rocks. The...     

Evolution of lithospheric mantle in the north of Nain‐Baft oceanic crust (Neo‐Tethyan ophiolite of Ashin,Central Iran)

This study is focused on a plagioclase‐bearing spinel lherzolite from Chah Loqeh area in the Neo‐Tethyan Ashin ophiolite. It is exposed along the west of left‐lateral strike‐slip Dorouneh Fault in the northwest of Central‐East Iranian Microcontinent. Mineral chemistry (Mg#^olivine < ~ 90, Cr#^{clinopyroxene} < ~ 0.2, Cr#^spinel < ~ 0.5, Al₂O₃^{orthopyroxene} > ~ 2.5 wt%, Al₂O₃^{clinopyroxene} > ~ 4.5 wt%, Al₂O₃^spinel > ~ 41.5 wt%, Na₂O^{clinopyroxene} > ~ 0.11 wt%, and TiO₂^{clinopyroxene} > ~ 0.04 wt%) confirms Ashin lherzolite was originally a mid‐oceanic ridge peridotite with low degrees of partial melting at spinel‐peridotite facies in a lithospheric mantle level. However, some Ashin lherzolites record mantle upwelling and tectonic exhumation at plagioclase‐peridotite facies during oceanic extension and diapiric motion of mantle along Nain‐Baft suture zone. This mantle upwelling is evidenced by some modifications in the modal composition (i.e. subsolidus recrystallization of plagioclase and olivine between pyroxene and spinel) and mineral chemistry (e.g. increase in TiO₂ and Na₂O of clinopyroxene, and TiO₂ and Cr# of spinel and decrease in Mg# of olivine), as a consequence of decompression during a progressive upwelling of mantle. Previous geochronological and geochemical data and increasing the depth of subsolidus plagioclase formation at plagioclase‐peridotite facies from Nain ophiolite (~ 16 km) to Ashin ophiolite (~ 35 km) suggest a south to north closure for the Nain‐Baft oceanic crust in the northwest of Central‐East Iranian Microcontinent.     

Coping with water scarcity in Kashafroud G-WADI Basin,Iran: climate change or growing demands?

This paper assesses the various factors contributing to climate change in the region of the Kashafroud G-WADI Basin in Iran; quantifies the local impacts of climate change, especially local water scarcity; and simulates and discusses several proposed methods to combat these impacts. Hydrologic and climatic data are statistically analyzed and VENSIM modeling is used for various simulations of water resources in the basin. Results show that the natural climate changes affecting Kashafroud Basin include increased temperature, less rainfall, more frequent droughts, and changes in rainfall patterns, all of which are local symptoms of climate change in recent years. However, the most important challenge in the basin is the overexploitation of surface and groundwater resources to meet the growing water demands, especially domestic needs. Changes in land use, reallocation of water uses, groundwater depletion, and degradation of the quality of surface waters have all contributed to significant changes in the environmental features of this basin, and are the main reason why water demands now exceed the renewal capacity of the basin. Proposed response measures include reallocation of resources among different uses, inter-basin water transfers, drawing water from six small dams on the Kashafroud River, reducing groundwater extraction, and replacing groundwater extraction for agriculture by reuse of urban wastewater. This study concludes that although changes in global climatic parameters have altered environmental features in the basin, local factors, such as water utilization beyond the renewable capacity of the basin, are more significant in worsening the impacts of climate change.     

Assessment of soil liquefaction potential based on numerical method- a case study of an urban railway project

Paying special attention to geotechnical hazards such as liquefaction in huge civil projects like urban railways especially in susceptible regions to liquefaction is of great importance. A number of approaches to evaluate the potential for initiation of liquefaction, such as Seed and Idriss simplified method have been developed over the years. Although simplified methods are available in calculating the liquefaction potential of a soil deposit and shear stresses induced at any point in the ground due to earthquake loading, these methods cannot be applied to all earthquakes with the same accuracy, also they lack the potential to predict the pore pressure developed in the soil. Therefore, it is necessary to carry out a ground response analysis to obtain pore pressures and shear stresses in the soil due to earthquake loading. Using soil historical, geological and compositional criteria, a zone of the corridor of Tabriz urban railway line 2 susceptible to liquefaction was recognized. Then, using numerical analysis and cyclic stress method using QUAKE/W finite element code, soil liquefaction potential in susceptible zone was evaluated based on design earthquake.     

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.     

Cretaceous tectonic evolution of the Neo-Tethys in Central Iran:Evidence from petrology and age of the Nain-Ashin ophiolitic basalts

The Nain and Ashin ophiolites consist of Mesozoic melange units that were emplaced in the Late Cretaceous onto the continental basement of the Central-East Iran microcontinent(CEIM).They largely consist of serpentinized peridotites slices;nonetheless,minor tectonic slices of sheeted dykes and pillow lavas-locally stratigraphically associated with radiolarian cherts-can be found in these ophiolitic melanges.Based on their whole rock geochemistry and mineral chemistry,these rocks can be divided into two geochemical groups.The sheeted dykes and most of the pillow lavas show island arc tholeiitic(IAT)affinity,whereas a few pillow lavas from the Nain ophiolites show calc-alkaline(CA)affinity.Petrogenetic modeling based on trace elements composition indicates that both IAT and CA rocks derived from partial melting of depleted mantle sources that underwent enrichment in subduction-derived components prior to melting.Petrogenetic modeling shows that these components were represented by pure aqueous fluids,or sediment melts,or a combination of both,suggesting that the studied rocks were formed in an arc-forearc tectonic setting.Our new biostratigraphic data indicate this arc-forearc setting was active in the Early Cretaceous.Previous tectonic interpretations suggested that the Nain ophiolites formed,in a Late Cretaceous backarc basin located in the south of the CEIM(the so-called Nain-Baft basin).However,recent studies showed that the CEIM underwent a counter-clockwise rotation in the Cenozoic,which displaced the Nain and Ashin ophiolites in their present day position from an original northeastward location.This evidence combined with our new data and a comparison of the chemical features of volcanic rocks from different ophiolites around the CEIM allow us to suggest that the Nain-Ashin volcanic rocks and dykes were formed in a volcanic arc that developed on the northern margin of the CEIM during the Early Cretaceous in association with the subduction,below the CEIM,of a Neo-Tethys oceanic branch that was existing between the CEIM and the southern margin of Eurasia.As a major conclusion of this paper,a new geodynamic model for the Cretaceous evolution of the CEIM and surrounding Neo-Tethyan oceanic basins is proposed.     

Alkaline basalt from the Central Iran, a mark of previously subducted Paleo-Tethys oceanic crust

In western part of the CEIM (Central-East Iranian Microcontinent) (Bayazeh area, Isfahan province, Iran), a series of Paleozoic basaltic rocks, occur. Major minerals of these basalts are olivine, clinopyroxene (diopside, augite), plagioclase (albite), sanidine, amphibole (kaersutite), phlogopite, ilmenite and magnetite. Secondary minerals include epidote, pumpellyite, albite, calcite and chlorite. Olivine and clinopyroxene are as phenocryst, while feldspars are restricted to groundmass. Chemical composition of clinopyroxenes indicates crystallization during ascending of magma. Geochemical analysis of whole rock samples shows that these rocks are characterized by low SiO₂ (43.21–48.45 wt %), high TiO₂ (1.81–3.00 wt %) and P₂O₅ (0.18–0.34 wt %). Petrography, chemistry of clinopyroxenes and whole rock analyses reveal an alkaline nature of these basalts. They are enriched in alkalis (Na₂O + K₂O = 4.1–7.7 wt %), LILE, HFSE and LREE. The Bayazeh alkali-basalts present strong enrichment in LREE relative to HREE (La/Lu ratio = 77.6–119.6) and were dominantly derived from partial melting of a metasomatized asthenospheric garnet-amphibole lherzolite. Field relationships reveal that junction of faults in west of the Bayazeh prepared a suitable path for ascending of magma from deep regions to surface and intra-plate continental magmatism. The Paleo-Tethys subduction from lower to upper Paleozoic is too enough for mantle enrichment in volatiles and basaltic alkaline magmatisrn in upper Paleozoic of Bayazeh area.     

Chalk–fluid interactions with glycol and brines

Mechanical properties of high porosity chalks are strongly dependent on the type of fluid in the pores. Dry chalk is considerably stronger than water-saturated chalk, and this phenomenon is often referred to as the water-weakening effect. To address the problem of chalk–fluid interactions, several series of tests have been performed with glycol and high concentration brines as saturating fluids. Glycol is a fluid that in many aspects resembles oil, except that glycol is miscible with water. Glycol-saturated chalk turns out to have properties very similar to oil-saturated chalk. Compared to dry chalk, both oil and glycol make the chalk somewhat weaker, but this weakening effect is much less than with water. Several series of tests with brines with high concentrations of calcium chloride or sodium chloride show that the water-weakening effect is considerably reduced in high ionic strength solutions. Most tests were performed as quasi-hydrostatic tests, with a constant stress ratio of 0.9. In such tests, the yield point marks the onset of accelerated pore collapse, and the yield value is close to the hydrostatic yield stress. In addition to these compressive tests, a series of Brazilian tests were performed, revealing the same trend. The variations in mechanical strength have been correlated with the activity of water in the brines. Within the experimental accuracy of the compressive tests, there is a linear trend between reduction of water activity and the corresponding increase in strength. This leads to the hypothesis that the water activity may be a key parameter in the water-weakening mechanism. But this conclusion also indicates that water weakening may be a special case of general chalk–fluid interactions where the degree of weakening depends on the strength of adsorption of the fluid molecules to the calcite surfaces.