Mineral chemical composition and geodynamic significance of peridotites from Nain ophiolite,central Iran

The peridotites from north of the town of Nain in central Iran consist of clinopyroxene-bearing harzburgite and lherzolite with small lenses of dunite and chromitite pods. The lherzolite contains aluminous spinel with a Cr number (Cr^# = Cr/[Cr + Al]) of 0.17. The Cr number of spinels in harzburgite and chromitite is 0.38–0.42 and 0.62, respectively. This shows that the lherzolite and harzburgite resulted from <18% of partial melting of the source materials. The estimated temperature is 1100 ± 200 °C for peridotites, the estimated pressure is <15 ± 2.3 kbar for harzburgites and >16 ± 2.3 kbar for lherzolites and estimated fo₂ is 10^?1±0.5 for peridotites. Discriminant geochemical diagrams based on mineral chemistry of harzburgites indicate a supra-subduction zone (SSZ) to mid-oceanic ridge (MOR) setting for these rocks. On the basis of their Cr^#, the harzburgite and lherzolite spinels are analogous to those from abyssal peridotites and oceanic ophiolites, whereas the chromites in the chromitite (on the basis of Cr^# and boninitic nature of parental melts) resemble those from SSZ ophiolitic sequences. Therefore, the Nain ophiolite complex most likely originated in an oceanic crust related to supra-subduction zone, i.e. back arc basin. Field observations and mineral chemistry of the Nain peridotites, indicating the suture between the central Iran micro-continent (CIM) block and the Sanandaj–Sirjan zone, show that these peridotites mark the site of the Nain–Baft seaway, which opened with a slow rate of ocean-floor spreading behind the Mesozoic arc of the Sanandaj–Sirjan zone as a result of change of Neo Tethyan subduction régime during middle Cretaceous.     

Electron Temperatures and Flow Speeds of the Low Solar Corona: MACS Results from the Total Solar Eclipse of 29 March 2006 in Libya

An experiment was conducted in conjunction with the total solar eclipse on 29 March 2006 in Libya to measure both the electron temperature and its flow speed simultaneously at multiple locations in the low solar corona by measuring the visible K-coronal spectrum. Coronal model spectra incorporating the effects of electron temperature and its flow speed were matched with the measured K-coronal spectra to interpret the observations. Results show electron temperatures of (1.10±0.05) MK, (0.70±0.08) MK, and (0.98±0.12) MK, at 1.1 R _⊙ from Sun center in the solar north, east and west, respectively, and (0.93±0.12) MK, at 1.2 R _⊙ from Sun center in the solar west. The corresponding outflow speeds obtained from the spectral fit are (103±92) km s⁻¹, (0+10) km s⁻¹, (0+10) km s⁻¹, and (0+10) km s⁻¹. Since the observations were taken only at 1.1 R _⊙ and 1.2 R _⊙ from Sun center, these speeds, consistent with zero outflow, are in agreement with expectations and provide additional confirmation that the spectral fitting method is working. The electron temperature at 1.1 R _⊙ from Sun center is larger at the north (polar region) than the east and west (equatorial region).     

The impact of land use and land cover changes on soil erosion in western Iran

Natural Hazards - Analysis of long-term land use and land cover (LULC) changes requires up-to-date remotely sensed data to assess their effects on erosion. This is a particularly important...     

Impact of climate change on runoff in Lake Urmia basin,Iran

Theoretical and Applied Climatology - Investigation of the impact of climate change on water resources is very necessary in dry and arid regions. In the first part of this paper, the climate model...     

Socio-hydrological framework for investigating farmers’ activities affecting the shrinkage of Urmia Lake; hybrid data mining and agent-based modelling

ABSTRACT

Developing a general framework to capture the complexities associated with the non-linear and adaptive nature of farmers facing water resources scarcity is a challenging problem. This paper integrates agent-based modelling (ABM) and a data mining method to develop a hybrid socio-hydrological framework to provide future insights for policy-makers. The data associated with the farmers’ main characteristics were collected through field surveys and interviews. Afterwards, the association rule was employed to discover the main patterns representing the farmers’ agricultural decisions. The discovered patterns were then used as the behavioural rules in ABM to simulate the agricultural activities. The proposed framework has been was applied to explore the interactions between agricultural activities and the main river feeding the Urmia-Lake, Iran. The outcomes indicate that farmers’ acquisitive traits and belongings have significant impacts on their socio-hydrological interactions. The reported values of the efficiency criteria may support the satisfactory performance of the proposed framework.     

Selection of Less Biased Threshold Angles for SAM Classification Using the Real Value–Area Fractal Technique

The accuracy of classification of the Spectral Angle Mapping (SAM) is warranted by choosing the appropriate threshold angles, which are normally defined by the user. Trial‐and‐error and statistical methods are commonly applied to determine threshold angles. In this paper, we discuss a real value–area (RV–A) technique based on the established concentration–area (C–A) fractal model to determine less biased threshold angles for SAM classification of multispectral images. Short wave infrared (SWIR) bands of the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) images were used over and around the Sar Cheshmeh porphyry Cu deposit and Seridune porphyry Cu prospect. Reference spectra from the known hydrothermal alteration zones in each study area were chosen for producing respective rule images. Segmentation of each rule image resulted in a RV–A curve. Hydrothermal alteration mapping based on threshold values of each RV–A curve showed that the first break in each curve is practical for selection of optimum threshold angles. The hydrothermal alteration maps of the study areas were evaluated by field and laboratory studies including X–ray diffraction analysis, spectral analysis, and thin section study of rock samples. The accuracy of the SAM classification was evaluated by using an error matrix. Overall accuracies of 80.62% and 75.45% were acquired in the Sar Cheshmeh and Seridune areas, respectively. We also used different threshold angles obtained by some statistical techniques to evaluate the efficiency of the proposed RV–A technique. Threshold angles provided by statistical techniques could not enhance the hydrothermal alteration zones around the known deposits, as good as threshold angles obtained by the RV–A technique. Since no arbitrary parameter is defined by the user in the application of the RV‐A technique, its application prevents introduction of human bias to the selection of optimum threshold angle for SAM classification.     

Did satellite imagery supersede aerial imagery? A perspective from 3D geopositioning accuracy

In this study, the geometric accuracy comparison of aerial photos and WorldView-2 satellite stereo image data is evaluated with the different number and the distribution of the ground control points (GCPs) on the basis of large scale map production. Also, the current situation of rivalry between airborne and satelliteborne imagery was mentioned. The geometric accuracy of Microsoft UltraCam X 45 cm ground sampling distance (GSD) aerial imagery and WorldView-2 data both with and without GCPs are also separately analyzed. The aerial photos without any GCP by only using global navigation satellite system (GNSS) and inertial measurement unit (IMU) data with tie points give an accuracy of ±1.17 m in planimetry and ±0.71 m in vertical that means nearly two times better accuracy than the rational polynomial coefficient (RPC) of stereo WorldView-2. Using one GCP affects the accuracies of aerial photos and WorldView-2 in different ways. While this situation distorts the aerial photo block, it corrects the shift effect of RPC in WorldView-2 and increases the accuracy. By using four or more GCPs, ½?pixel (～0.23 m) accuracy in aerial photos and 1 pixel (～0.50 m) accuracy in WorldView-2 can be achieved in horizontal. In vertical, aerial photos have 1 pixel (～0.55 m) and WorldView-2 has 1.5 pixels (～0.85 m) accuracy. These results show that Worldview-2 imagery can be used in the production of class I 1:5000 scale maps according to the ASPRS Accuracy Standards for Digital Geospatial Data in terms of geometric accuracy. It is concluded that the rivalry between aerial and satellite imagery will continue for some time in the future.     

High resolution facies architecture and digital outcrop modeling of the Sandakan formation sandstone reservoir,Borneo:Implications for reservoir characterization and flow simulation

Advances in photogrammetry have eased the acquisition of high-resolution digital information from outcrops, enabling faster, non-destructive data capturing and improved reservoir modeling. Geocellular models for flow dynamics with in the virtual outcrop in siliciclastic deposits at different sets of sandstone facies architecture remain, however, a challenge. Digital maps of bedding, lithological contrast, spatial-temporal variations of bedding and permeability characteristics make it more easy to understand flow tortuosity in a particular architecture. An ability to precisely model these properties can improve reservoir characterization and flow modeling at different scales. Here we demonstrate the construction of realistic 2 D sandstone facies based models for a pragmatic simulation of flow dynamics using a combination of digital point clouds dataset acquired from LiDAR and field investigation of the Sandakan Formation, Sabah, Borneo.Additionally, we present methods for enhancing the accuracy of outcrop digital datasets for producing high resolution flow simulation. A well-exposed outcrop from the Sandakan Formation, Sabah, northwest Borneo having a lateral extent of 750 m was chosen in order to implement our research approach. Sandstone facies and its connectivity are well constrained by outcrop observations, data from air-permeability measurements, bilinear interpolation of permeability, grid construction and water vector analysis for flow dynamics.These proportions were then enumerated in terms of static digital outcrop model(DOM) and facies model based on sandstone facies bedding characteristics. Flow simulation of water vector analysis through each of the four sandstone facies types show persistent spatial correlation of permeability that align with either cross-bedded orientation or straight with more dispersion high quality sandstone(porosity 21.25%-41.2%and permeability 1265.20-5986.25 mD) and moderate quality sandstone(porosity 10.44%-28.75% and permeability 21.44-1023.33 mD). Whereas, in more heterolithic sandstone(wavy-to flaser-bedded and bioturbated sandstone), lateral variations in permeability show spatially non-correlated patterns over centimeters to tens of meters with mostly of low quality sandstone(porosity 3.4%-12.31% and permeability < 1 mD to 3.21 mD). These variations reflect the lateral juxtaposition in flow dynamics. It has also been resulted that the vertical connectivity and heterogeneities in terms of flow are mostly pragmatic due to the interconnected sandstone rather than the quality of sandstone.     

Computing derivative information of sequentially coupled subsurface models

A generic framework for the computation of derivative information required for gradient-based optimization using sequentially coupled subsurface simulation models is presented. The proposed approach allows for the computation of any derivative information with no modification of the mathematical framework. It only requires the forward model Jacobians and the objective function to be appropriately defined. The flexibility of the framework is demonstrated by its application in different reservoir management studies. The performance of the gradient computation strategy is demonstrated in a synthetic water-flooding model, where the forward model is constructed based on a sequentially coupled flow-transport system. The methodology is illustrated for a synthetic model, with different types of applications of data assimilation and life-cycle optimization. Results are compared with the classical fully coupled (FIM) forward simulation. Based on the presented numerical examples, it is demonstrated how, without any modifications of the basic framework, the solution of gradient-based optimization models can be obtained for any given set of coupled equations. The sequential derivative computation methods deliver similar results compared to FIM methods, while being computationally more efficient.