Simulating winter wheat production potential under near-future climate change in arid regions of northeast Iran

Theoretical and Applied Climatology - Water scarcity is the key challenge in arid regions, which exacerbates under climate change (CC) and must be considered to assess the impacts of CC on cropping...     

Water leakage paths in the Doosti Dam,Turkmenistan and Iran

The Doosti Dam, with a reservoir capacity of 1,250 million cubic meters, was constructed on the Harirood River at the border of Turkmenistan and Iran. The reservoir is in direct contact with permeable formations on the right abutment of the dam including the Neyzar Sandstone, the Kalat Limestone and the Pesteleigh alternative layers of marlstone and sandstone. After the reservoir impoundment, several new springs and seepages emerged from these formations and the alluvium. The amount of leakage increased with the rise in reservoir water level. Fifteen kilograms of sodium fluorescein were injected into a 113-m deep borehole intersecting three permeable sandstone layers of the Pesteleigh formation. Dye was detected downstream of the grout curtain in boreholes and springs that were in direct contact with parts of the Pesteleigh formation having the same sandstone layers as the injection borehole. The dye velocity was in the range of diffuse flow, confirming the good performance of the grout curtain in the Pesteleigh formation. No dye was detected in the other formations because the injection borehole was not in direct contact with these formations. The hydraulic relation of the other formations with the reservoir was determined by considering direct contact of the formations with the reservoir, emergence of new springs and seepages after reservoir impoundment, correlation of time variations of the springs discharge and the borehole’s water level with the reservoir water level, and in some cases the hydrochemistry of the water. The results show that the Kalat and Neyzar formations are hydraulically connected to the reservoir, but the small amounts of leakage from these formations at a hydraulic gradient of 24% indicates good performance of the built grout curtain. The total reservoir leakage at maximum reservoir water level was 100 l/s which is insignificant compared with the 15 m³/s average annual release of the reservoir.     

A modi?ed failure criterion for transversely isotropic rocks

A modified failure criterion is proposed to determine the strength of transversely isotropic rocks. Me-chanical properties of some metamorphic and sedimentary rocks including gneiss, slate, marble, schist, shale, sandstone and limestone, which show transversely isotropic behavior, were taken into consider-ation. Afterward, introduced triaxial rock strength criterion was modified for transversely isotropic rocks. Through modification process an index was obtained that can be considered as a strength reduction parameter due to rock strength anisotropy. Comparison of the parameter with previous anisotropy in-dexes in literature showed reasonable results for the studied rock samples. The modified criterion was compared to modified Hoek-Brown and Ramamurthy criteria for different transversely isotropic rocks. It can be concluded that the modified failure criterion proposed in this study can be used for predicting the strength of transversely isotropic rocks.     

The effects of hydrometer reading times on the spatial variability of soil textures in southeast Iran

Soil texture is an important physical soil property that may contribute to variations in many soil functions as well as nutrient storage and availability, water retention, and soil erosion. Although several methods for determining the texture classes of soil particles have been proposed, differences among hydrometer reading times have presented challenges in determining the precise soil texture classes. Therefore, this study was conducted to evaluate the effects of hydrometer reading time on the spatial variability of soil textures in the Rafsanjan area, southeast Iran. To accomplish this, 77 soil samples were collected on a 500-m square sampling grid from depths of 0–40, 40–80, and 80–120 cm, and their particle sizes were determined through analysis for 40 s, 2 h, 6.5 h, and 8 h using the Bouyoucos hydrometer method. The results showed a strong spatial correlation in the soil particles among sampling soil layers and across the study area. Moreover, the differences among hydrometer reading times did not have a significant impact on determination of coarse soil texture classes, although they did influence determination of the finer classes. Although the 8 h reading time provided the most accurate response with respect to mechanical analysis of a soil, after 6.5 h the hydrometer could also largely (more than 80.0 %, on average) achieve this goal. Additionally, the 2 h hydrometer reading time could also be useful for the initial assessment or general overview of the soil texture in a certain region; however, it is not recommended for precision agriculture or site-specific management.     

Seismic piezocone interpretation for shear wave velocity (<Emphasis Type="Italic">V</Emphasis><Subscript>S</Subscript>) determination in the Persian Gulf

Shear wave velocity is one of the important factors representing the dynamic characteristics of soil layers. Hence, many researchers have focused their studies on determining shear wave velocity by direct field measurements or expressions developed by other soil parameters. The shear module and damping ratio of the soil layers also play a similar role in the majority of dynamic soil response analyses. Nevertheless, since they have to be measured in the laboratory by resonant column or cyclic triaxial tests on undisturbed samples, the possibility of preparing such samples and the reliability of the obtained results are of great concerns. In the present study, great effort has been made to determine the above dynamic factors by means of field data obtained from a versatile instrument, namely the seismic piezocone (SPCT_U), and to derive expressions correlating them with some parameters obtainable by much simpler instruments. The reliability of laboratory measurements on undisturbed samples is also evaluated. The seismic piezocone test apparatus has been employed to evaluate the soil properties at 1-m depth intervals by means of measuring tip resistance, sleeve resistance, pore pressure and shear wave velocity. The shear module and the damping ratio are calculated using field data. Meanwhile, in order to assess the laboratory measurements of these parameters, some resonant column tests and cyclic triaxial tests on undisturbed samples of the same soil layers have been carried out. In order to compare the field results of shear modulus and damping ratios with those obtained from laboratory tests, the influences of the soil nature and sample disturbance on the conventional laboratory methods are evaluated and discussed. The shear wave velocity is correlated to overburden pressure and the corrected tip resistance for two groups of fine soils, namely silty clays and carbonate clayey silts, which mainly cover the areas under study in this project, are located in southern parts of Iran near the Persian Gulf. According to the results of the present study, there are narrow limits of shear modulus regarding soils for which the laboratory tests and the field measurements yield approximately the same shear modulus. This limit of shear modulus is about 30–50(MPa) for clay deposits and 70–100 (MPa) for sandy deposits. Also the shear wave velocity can be calculated by a simple expression from total overburden pressure and the tip resistance of simple cone penetration test results conventionally available in many soil explorations prior to engineering practices. However, if the pore pressure inside the saturated soil deposits can be measured by a piezocone apparatus, the shear wave velocity may be calculated using another suggested equation in terms of effective overburden pressure in the present study. Regarding the shear module and the damping ratio, due to the disturbances of the stiff deposits in the sampling process and great deviations of laboratory results from field results, the laboratory measurements of these parameters out of the above limits are not recommended.     

Shoreline spatial and temporal response to natural and human effects in Boujagh National Park,Iran

Shoreline variation and river deltas are among the most dynamic systems in marine environments. The related different variations in spatial and temporal scales play significant roles in land planning and different management applications. Modeling the dynamics of seashore of Boujagh National Park (BNP) which is located on the southern coast of the Caspian Sea in the Sefidrud Delta (SD), considering natural and anthropogenic factors, was the main objective of the current study. To achieve this goal, a combination of remote sensing data, historical data, and numerical simulations was utilized. The BNP covers an area of 3,270 ha and includes two international wetlands, Boujagh and Kiashahr. In earlier periods, this area faced severe morphological changes whereas recently its shoreline has experienced gradual variations. Accordingly, at the first stage, the shoreline variation from 2006 to 2017 was extracted by processing and classifying Operational Land Imager (OLI) and Thematic Mapper (TM) images from Landsat satellites using the Maximum Likelihood approach. In the second stage, the two dimensional MIKE21 model was utilized to identify wave and coastal current patterns and parameters for the year 2015. Morphologically, the results showed that, the shoreline of the BNP is affected by several natural and anthropogenic factors. Seaward advancement of the shoreline occurred in zones A (east zone) and C (west zone) due to Caspian Sea Level drop and sedimentation while retreating occurred at Zone B (north zone) influenced by wave and current patterns and reduction of the Sefidrud River flows. Also, the results imply that maintaining the existing conditions results in the disappearance of a considerable part of the ecological area in the BNP. Hence, to manage and preserve the coastline of the BNP complying with the current anthropogenic and natural factors, it is vital to take necessary management measures.     

A framework for predicting abrasion rupture of crusts in wind erosion

Crusts play a crucial role in the reduction or control of wind erosion. In this regard, the resilience and durability of crusts are of prime importance. Crusts have high resilience and durability against wind flow shear stresses; however, they are prone to abrasion induced by saltating particles. Therefore, estimating crust durability in abrasion rupture has practical importance. In this study, a cyanocrust and a biocemented sand crust were subjected to a controlled flux of saltating particles for different sandblasting periods to provide a framework for predicting crust rupture. The velocity and pre- and post-collision energy of the saltating particles were measured using high-speed photography. The changes in the strength of the crusts after different periods of sandblasting were determined using a scratch test. The results suggested that the average strength of the cyanocrust and biocemented sand crust became 0.25 and 0.7 of their corresponding initial values after 30 min of sandblasting. Also, the average stiffness of the cyanocrust and biocemented sand crust decreased to 0.5 and 0.9 of their initial values, respectively. Furthermore, the amount of impact energy absorbed by the crusts increased by the deterioration of the crusts. Compiling the results of the wind tunnel experiment and scratch tests yielded an exponential equation which can be used to estimate crust durability in a given condition of saltation. Based on this equation, the cyanocrust and biocemented sand crust will break down entirely after 23 and 449 min, respectively, at a wind velocity of 6.8 m/s and a saltation flux of 1 g/s/m.     

Development of a New Index to Assess the Rock Mass Drillability

Knowledge of drillability of rock masses in engineering projects is very important in determining drilling costs. In drilling operations, so many parameters such as the properties of rock and the drilling equipment affect the drilling performance. In this study, after discussing the rock mass drillability process and identifying all the effective parameters, interaction matrixes based on the rock engineering systems, that analyze the interrelationship between the parameters affecting rock engineering activities, is introduced to study the rock mass drillability tribosystem. Given that interaction matrix codes are not unique numbers, and then possible interactive intensities are calculated for each matrix and a group decision-making method, Fuzzy–Delphi–AHP technique has been used to obtain appropriate weights. As a result, rock mass drillability index (RMDI) is presented to classify the rock mass drillability. The results indicate the ability of this method to analyze rock mass drillability procedure. Drilling data along with laboratory rock properties from Sungun copper mine were collected and were ranked according to the new classification system. Fifteen zones at the mine site were ranked based upon the new index RMDI and a reasonable correlation was obtained between measured drilling rate at the zones and RMDI data.     

The magmatic record in the Arghash region (northeast Iran) and tectonic implications

The area of Arghash in northeast Iran, prominent for its gold mineralization, was newly mapped on a scale of 1:20,000 with particular attention to the occurring generations of igneous rocks. In addition, geochronological and geochemical investigations were carried out. The oldest geological unit is a late Precambrian, hornblende-bearing diorite pluton with low-K composition and primitive isotope signatures. This diorite (U–Pb zircon age 554 ± 6 Ma) is most likely a remnant from a Peri-Gondwana island-arc or back-arc basin. About one-third of the map area is interpreted as an Upper Cretaceous magmatic arc consisting of a volcanic and a plutonic part. The plutonic part is represented by a suite of hornblende-bearing medium-K, I-type granitoids (minor diorite, mainly quartz–monzodiorite and granodiorite) dated at 92.8 ± 1.3 Ma (U–Pb zircon age). The volcanic part comprises medium-K andesite, dacite and tuffitic rocks and must be at least slightly older, because it is locally affected by contact metamorphism through the hornblende–granitoids. The Upper Cretaceous arc magmatism in the Arghash Massif is probably related to the northward subduction of the Sabzevar oceanic basin, which holds a back-arc position behind the main Neotethys subduction front. Small occurrences of pillow basalts and sediments (sandstone, conglomerate, limestone) tectonically intercalated in the older volcanic series may be relics of earlier Cretaceous or even pre-Cretaceous rocks. In the early Cenozoic, the Cretaceous magmatic arc was intruded by bodies of felsic, weakly peraluminous granite (U–Pb zircon age 55.4 ± 2.3 Ma). Another strong pulse of magmatism followed slightly later in the Eocene, producing large masses of andesitic to dacitic volcanic rocks. The geochemistry of this prominent Eocene volcanism is very distinct, with a high-K signature and trace element contents similar to shoshonitic series (high P, Zr, Cr, Sr and Ba). High Sr/Y ratios feature affinities to adakite magmas. The Eocene magmatism in the Arghash Massif is interpreted as related to thermal anomalies in crust and mantle that developed when the Sabzevar subduction system collapsed. The youngest magmatic activities in the Arghash Massif are lamprophyres and small intrusions of quartz–monzodiorite porphyries, which cut through all other rocks including an Oligocene–Miocene conglomerate cover series.     

Capability of ETM+ in estimates standing stock in beech stands of mountain forest

Stand volume is an important criterion in forest sciences for monitoring status and function of forests, estimation of productivity, prediction and modeling of forest disturbance, economic and environmental issues, and forest planning. The aim of this research is to evaluate the ETM+ sensor of Landsat 7 satellite data ability for forest timber volume estimation. For this purpose, 40 selective sample plots with 60?×?60 m dimension were selected, and in each sample plot, standing volume was calculated. Correspondent digital data to plots were extracted from spectral and considered as independent variables. Original stand volume data, square root, and logarithm of them were considered as dependent potential for stand volume estimation. Variables using stepwise regression, the best model with coefficient of determination, and adjusted coefficient of determination with regard to determining the appropriate model were modified Log?V=?07.546?+?0.85Greenness???10.564NDWI???28 0.063Band1. In this model, the coefficient of determination and adjusted coefficient of determination were obtained 0.778 and 0.770, respectively. Result showed that spectral data of the mentioned sensor have a moderate.