On the advection of sharp material interfaces in geodynamic problems: entrainment of the D″ layer

The D″ layer is a dense and chemically distinct layer at the base of the convecting mantle. Numerical modeling of the entrainment of this layer by mantle convection requires the solution of the advective transport equation without introducing numerical diffusion across sharp material boundaries. We use our improved second moment numerical method to solve the equation. The method conserves the amount of material and the first and second moments of material distribution in each control volume. We first consider two examples of isothermal Rayleigh–Taylor instability to illustrate the performance of our method by comparing our results with those of a number of field, tracer and marker chain methods. We show that the performance of our method in minimizing the numerical diffusion is better than the field methods and comparable to the tracer and marker chain methods. We then study the instability of the dense D″ layer and its interaction with the overlying mantle. A range of density contrast between the D″ layer and the mantle, layer thickness, and the Rayleigh number, Ra, is examined. We show that for higher values of these parameters, the amount of entrainment decreases and the layer remains stable over longer periods of time. For very thick D″ layers and high Ra values, internal convection can take place within the layer.     

New evidence for Jurassic continental rifting in the northern Sanandaj Sirjan Zone,western Iran: the Ghalaylan seamount,southwest Ghorveh

ABSTRACT

We address the growing controversy about the tectonic setting in which Jurassic magmatism of Iran occurred: arc or continental rift. In the Ghorveh area of the northern Sanandaj Sirjan zone (SaSZ), the Ghalayan metabasites are interlayered with marble and schist and locally cut by acidic dikes. Zircon U-Pb dating of the metabasitic rocks shows that these crystallized at ca. 145–144 Ma ago in the Late Jurassic (Tithonian). This complex was metamorphosed in the lower greenschist facies, however, some protolithic structures such as pillow lava and primary minerals are preserved. The metabasites are tholeiites with low SiO₂ (45.6–50.5 wt.%), moderate Al₂O₃ (11.3–17.0 wt.%), and high TiO₂ (0.7–2.9 wt.%) and Fe₂O₃ (9.4–14.1 wt.%). The Ghalayan metabasites are enriched in Light rare earth elements (LREEs) without significant Nb, Ta, Pb, Sr and Ba anomalies, similar to modern continental intra-plate tholeiitic basalts such as Afar and East African rifts. The Ghalaylan metabasites show wide ranges for ⁸⁷Sr/⁸⁶Sr_(i) (0.7039–0.7077) and positive ε_Nd(t) values (+0.1 to +4.6). These isotopic compositions are similar to those expected for slightly depleted subcontinental lithospheric mantle sources. Independently built discrimination diagrams indicate an intra-continental rifting regime for the source of Jurassic metabasites in the northern SaSZ. Geochemical and tectonic evidence suggests that rifting or a mantle plume was responsible for volcanic activity in the Upper Jurassic SaSZ. Considering the variation of ages of basaltic volcanism along the SaSZ, we suggest that Ghalayan basaltic magmatism reflected a submarine volcano that formed as part of the late stage continental rift, similar to Afar in the East African Rift system. Our results indicate that an extensional tectonic regime dominated SaSZ tectonics in the Middle to Late Jurassic.     

Envrionmental impact assessment of urban development plan by vulnerability model application

Urban development plans have significant potential impacts on affected environment according to the various activities and due to their different construction and operation phases. Establishment and development of the new-established Municipality of Tehran District-22, with about 10,000 hectares area has been planned as the last limit of continuous extension of Great Tehran and as the largest urban city of the country. It is based on the general lay out of Tehran master development plan, ratified in 1970 and 1991. In terms of current facilities, it has been anticipated necessary spaces to establish various activities including residential, commercial and urban services, establishing green spaces, man-made lake and other necessary infrastructures of urban development. Griding method was used in order to determine the environment vulnerability of the development plan of the region. In this method, the studied district has been divided into 1 to 10 square kilometers grids. To improve the accuracy of the study, the grids dimension is smaller than the common size. After griding, the most appropriate layers out of the required layers have been used to determine the vulnerability. These layers are, land slope, geographical aspect, pedology, geology, plant density and ground water vulnerability. The limited codes of each layer have also been defined. Summary of the above-mentioned layers results have been analyzed in the map form 1:50,000 scale by Arc-View software in GIS system in order to determine the specification of impact units. After the determination of limitation, the most prevailing limited code of each grid was selected. The obtained results were transferred to excel program and the final amount of vulnerability were then calculated by using equation E =∑ (a-b)/4. The final vulnerability was substitute by each grid code as a new data and the last map of vulnerability was set. Regarding to this map at the studied region affected by current limitation of environmental condition, the north part of the district especially A₉, B₁₀, B₁₁ B₁₂, C₉, C₁₁, C₁₈, D₁₂ and E₁₄ showed the highest amount of vulnerability. According to the final vulnerability map, it is distinguished that the north parts of the district mainly in parts 1 or 2 vulnerability codes have had high vulnerability.     

Agricultural activities impact on groundwater nitrate pollution

Concern over agricultural diffuse pollution sources in integrated water quality management has been growing recently. High nitrogen fertilizers application rates may increase the potential groundwater pollution. These effects were investigated in Andimeshk and Susa plains that cover an area of 1100 km² between the Dez and Karkhe rivers in north of Khozestan-Iran. This region divided to 4 sub-regions A, B, C, and D. Additionally 168 groundwater samples were collected from 42 water wells during the months April, May, August, and September of 2004. The Hackspectrophotometer nitrate test was used to measure the NO₃ ^? concentration in water samples. Information about further nitrate data was obtained. A questionnaire procedure was used for collection N-fertilizers application rate data in studied area. The results demonstrated that all of the groundwater samples have NO₃ ^t- concentration below the EPA MCL (44.27 mg/l) and WHO guideline (50 mg/l). The mean nitrate concentrations are 16.1, 19.5, 13.3, and 7.9 mg/l in sub-regions A, B, C, and D respectively. There are different amount of N-fertilizers applied in sub-regions A, B, C, and D. Correlation between NO₃t- concentrations and N-fertilizers rate suggests a inverse correlation between Nfertilizers application rate and ground waters nitrate concentrations in studied area (r=?0.69).     

Classification of Satellite Images in Assessing Urban Land Use Change Using Scale Optimization in Object-Oriented Processes (A Case Study: Ardabil City,Iran)

By using satellite imagery, the recognition and evaluation of various phenomena and extraction of information necessary for the planning of land resources or other purposes are easily accomplished. The purpose of this study is to compare the efficiency of seven commonly used methods of monitored classification of satellite data to evaluate land use changes using TM and OLI Landsat, IRS, Spot5 and Quick Bird bands as well as different color combinations of these images to detect agricultural land, residential areas and aquatic areas using object-oriented processing. Digital processing of satellite images was carried out in 1998 and 2016 using advanced methods. Training samples were extracted in five user classes by eCognition software using segmentation scale optimization, different color combinations and coefficients of shape and compression. The appropriate segmentation scale for arable land, human complications and the blue areas were, respectively, 50, 8 and 10. Then each image was classified separately using seven methods and extracted samples, and efficiency of each classification method was obtained by calculating two general health and Kappa coefficients. The results show that the accuracy of each classification method and the neural network with a total accuracy of 94.475 and Kappa coefficient of 92.095 were selected as the most accurate classification method. These results show that the sampling of educational samples with proper precision of the classes in the images and dependency probability of each satellite images pixel can be useful in classifying group available in helpful area.     

Sensor modeling,self-calibration and accuracy testing of panoramic cameras and laser scanners

Terrestrial Linear Array CCD-based panoramic cameras have been used for purely imaging purposes, but they also have a high potential for use in high accuracy measurement applications. The imaging geometry and the high information content of those images make them suitable candidates for quantitative image analysis. For that a particular sensor model has to be established and the inherent accuracy potential has to be investigated. We developed a sensor model for terrestrial Linear Array-based panoramic cameras by means of a modified bundle adjustment with additional parameters, which models substantial deviations of a real camera from the ideal one. We used 3D straight-line information in addition to tie points to conduct a full calibration and orientation without control point information. Due to the similarity of the operation of laser scanners to panoramic cameras the sensor model of the panoramic cameras was extended for the self-calibration of laser scanners. We present the joint sensor model for panoramic cameras and laser scanners and the results of self-calibration, which indicate a subpixel accuracy level for such highly dynamic systems. Finally we demonstrate the systems’ accuracy of two typical panoramic cameras in 3D point positioning, using both a minimal number of control points and a free network adjustment. With these new panoramic imaging devices we have additional powerful sensors for image recording and efficient 3D object modeling.     

Energy partitioning during seismic slip in pseudotachylyte-bearing faults (Gole Larghe Fault, Adamello, Italy)

The determination of the earthquake energy budget remains a challenging issue for Earth scientists, as understanding the partitioning of energy is a key towards the understanding the physics of earthquakes. Here we estimate the partition of the mechanical work density into heat and surface energy (energy required to create new fracture surface) during seismic slip on a location along a fault. Earthquake energy partitioning is determined from field and microstructural analyses of a fault segment decorated by pseudotachylyte (solidified friction-induced melt produced during seismic slip) exhumed from a depth of ~ 10 km—typical for earthquake hypocenters in the continental crust. Frictional heat per unit fault area estimated from the thickness of pseudotachylytes is ~ 27 MJ m^− 2. Surface energy, estimated from microcrack density inside clast (i.e., cracked grains) entrapped in the pseudotachylyte and in the fault wall rock, ranges between 0.10 and 0.85 MJ m^− 2. Our estimates for the studied fault segment suggest that ~ 97–99% of the energy was dissipated as heat during seismic slip. We conclude that at 10 km depth, less than 3% of the total mechanical work density is adsorbed as surface energy on the fault plane during earthquake rupture.     

Origin of 1.8 Ga zircons in Post Eocene mafic dikes in the Roshtkhar area,NE Iran

The Roshtkhar area is located in the Khaf-Kashmar-Bardaskan volcano-plutonic belt to the northeastern Iran along the regional E–W trending Dorouneh Fault, northeastern of the Lut Block. There are several outcrops of subvolcanic rocks occurring mainly as dikes in the area, which intruded into Cenozoic intrusive rocks. We present U–Pb dating of zircons from a diabase dike and syenite rock using LA-ICP-MS that yielded an age of 1778 ± 10 Ma for the dike, indicating this Cenozoic dike has zircon xenocrysts inherited from deeper sources; and 38.0 ± 0.5 Ma, indicating an Late Eocene crystallization age for the syenite. Geochemically, the dikes typical of high-K calc-alkaline to shoshonitic magmas. Petrographic observations and major and trace element variations suggest that diabase melts underwent variable fractionation of clinopyroxene, olivine, and Fe-Ti oxides and minor crustal contamination during the differentiation process. Primitive mantle-normalized multi-element diagrams display enrichment in LILE, such as Rb, Ba, Th, U, and Sr compared to HFSE, as well as negative anomalies of Nb, Ta, P, and Ti, suggesting derivation from subduction-modified mantle. Chondrite-normalized REE plots show moderately LREE enriched patterns (<3.83 La_N/Yb_N <8.27), and no significant Eu anomalies. Geochemical modelling using Sm/Yb versus La/Yb and La/Sm ratios suggests a low-degree of batch melting (~1–3%) of a phlogopite-spinel peridotite source to generate the mafic dikes. The geochemical signatures suggest that the Roshtkhar mafic dikes cannot be related directly to subduction and likely resulted from melting of upper mantle in an extensional setting where the heat flow was provided from deeper levels. These dikes presumably derived the zircon xenocrysts from the assimilation of upper crust of Gondwanian basement. Processes responsible for partial melting of metasomatized lithospheric mantle and post-collision magmatism in NE Iran was triggered by heating due to asthenospheric upwelling in an extensional setting.