Separation of magnetic anomalies into induced and remanent magnetization contributions

Inversion of magnetic data is complicated by the presence of remanent magnetization, and it provides limited information about the magnetic source because of the insufficiency of data and constraint information. We propose a Fourier domain transformation allowing the separation of magnetic anomalies into the components caused by induced and remanent magnetizations. The approach is based on the hypothesis that each isolated source is homogeneous with a uniform and specific Koenigsberger ratio. The distributions of susceptibility and remanent magnetization are subsequently recovered from the separated anomalies. Anomaly components, susceptibility distribution and distribution of the remanent and total magnetization vectors (direction and intensity) can be achieved through the processing of the anomaly components. The proposed method therefore provides a procedure to test the hypotheses about target source and magnetic field, by verifying these models based on available information or a priori information from geology. We test our methods using synthetic and real data acquired over the Zhangfushan iron-ore deposit and the Yeshan polymetallic deposit in eastern China. All the tests yield favourable results and the obtained models are helpful for the geological interpretation.     

Lithology-dependent seismic anisotropic amplitude variation with offset correction in transversely isotropic media

Analysis of amplitude variation with offset is an essential step for reservoir characterization. For an accurate reservoir characterization, the amplitude obtained with an isotropic assumption of the reservoir must be corrected for the anisotropic effects. The objective is seismic anisotropic amplitude correction in an effective medium, and, to this end, values and signs of anisotropic parameter differences (Δδ and Δε) across the reflection interfaces are needed. These parameters can be identified by seismic and well log data. A new technique for anisotropic amplitude correction was developed to modify amplitude changes in seismic data in transversely isotropic media with a vertical axis of symmetry. The results show that characteristics of pre-stack seismic data, that is, amplitude variation with offset gradient, can be potentially related to the sign of anisotropic parameter differences (Δδ and Δε) between two layers of the reflection boundary. The proposed methodology is designed to attain a proper fit between modelled and observed amplitude variation with offset responses, after anisotropic correction, for all possible lithofacies at the reservoir boundary. We first estimate anisotropic parameters, that is, δ and ε, away from the wells through Backus averaging of elastic properties resulted from the first pass of isotropic pre-stack seismic inversion, on input data with no amplitude correction. Next, we estimate the anisotropic parameter differences at reflection interfaces (values and signs of Δδ and Δε). We then generate seismic angle gather data after anisotropic amplitude correction using Rüger's equation for the P-P reflection coefficient. The second pass of isotropic pre-stack seismic inversion is then performed on the amplitude-corrected data, and elastic properties are estimated. Final outcome demonstrates how introduced methodology helps to reduce the uncertainty of elastic property prediction. Pre-stack seismic inversion on amplitude-corrected seismic data results in more accurate elastic property prediction than what can be obtained from non-corrected data. Moreover, a new anisotropy attribute (ν) is presented for improvement of lithology identification.     

Latitudinal and Longitudinal Variations of Earth’s Magnetic Force on Electrified Hydrometeors

Geomagnetism and Aeronomy - In this study, a hypothesis is proposed about the possible effect of Geomagnetic field (GMF) on the charge structure of a thundercloud based on Lorentz force equation...     

An integrated object-based image analysis and CA-Markov model approach for modeling land use/land cover trends in the Sarab plain

The present paper is an attempt to integrate a semi-automated object-based image analysis (OBIA) classification framework and a cellular automata-Markov model to study land use/land cover (LULC) changes. Land use maps for the Sarab plain in Iran for the years 2000, 2006, and 2014 were created from Landsat satellite data, by applying an OBIA classification using the normalized difference vegetation index, salinity index, moisture stress index, soil-adjusted vegetation index, and elevation and slope indicators. The classifications yielded overall accuracies of 91, 93, and 94% for 2000, 2006, and 2014, respectively. Finally, using the transition matrix, the spatial distribution of land use was simulated for 2020. The results of the study revealed that the number of orchards with irrigated agriculture and dry-farm agriculture in the Sarab plain is increasing, while the amount of bare land is decreasing. The results of this research are of great importance for regional authorities and decision makers in strategic land use planning.     

Extraction of forward-directivity velocity pulses using S-Transform-based signal decomposition technique

Robust methods for time-frequency analysis of time series, which provide local information of signals, allow earthquake engineers to study both the input and output of dynamic time history analysis with more reliability. Moreover, time-frequency representations (TFRs) have a major role in the analysis of non-stationary seismic signals exhibiting significant time variation of frequency content. S-Transform (ST) is a modern TFR, which can measure local characteristics of a signal such as amplitude, frequency, and phase at any time instant. This paper presents a new method for decomposition of ground motion signals. A modified version of ST-based technique, originally employed to decompose signals of gearbox vibration, is introduced and applied to the extraction and characterization of pulse-like part of near-fault velocity records, which is contributed to the directivity effects. In addition, a new definition based on ST analysis is used to identify pulse period. The results of implementation of proposed procedure on a database of pulse-like ground motion recordings belonging to the different ranges of magnitude demonstrate the efficiency of proposed method compared with other available approaches. The results, also, indicate that simple approximation of distinct pulses using single-period waveforms, unlike the extracted pulses, cannot represent the impulsive nature of real records adequately.     

Forest cover change assessment in conflict-affected areas of northwest Pakistan: the case of swat and shangla districts

This study examines the spatial and temporal forest cover changes in Swat and Shangla districts to understand the deforestation pattern in context of the recent security conflict in these districts. We used multi-resolution satellite images to assess the long term deforestation from 2001 to 2009 and also to identify episodic forest cutting areas appeared during the conflict period of Oct. 2007 - Oct. 2008. There are only 58 ha of deforestation identified during the conflict period while 1268 ha of gross annual deforestation were assessed during last eight years. Most of the deforestation patches persist around the administrative boundaries at sub-district levels (tehsils) which can be attributed to ambiguity in unclear jurisdiction between the forest official. The results highlight that the forest cutting appeared in Swat and Shangla during the conflict period is not as significant when compared with the long term deforestation pattern in the area. On the one side the results of the study are supportive to the picture that emerges from international studies which report high rate of deforestation in the country and on the other side it negates any relation between the security situation and the increasing deforestation in the north western Pakistan. The study concludes that deforestation assessments require verification by independent sources of data, such as satellite imagery to improve our understanding of deforestation processes.     

Higher-order compositional modeling with Fickian diffusion in unstructured and anisotropic media

We present advances in compositional modeling of two-phase multi-component flow through highly complex porous media. Higher-order methods are used to approximate both mass transport and the velocity and pressure fields. We employ the Mixed Hybrid Finite Element (MHFE) method to simultaneously solve, to the same order, the pressure equation and Darcy's law for the velocity. The species balance equation is approximated by the discontinuous Galerkin (DG) approach, combined with a slope limiter. In this work we present an improved DG scheme where phase splitting is analyzed at all element vertices in the two-phase regions, rather than only as element averages. This approximation is higher-order than the commonly employed finite volume method and earlier DG approximations. The method reduces numerical dispersion, allowing for an accurate capture of shock fronts and lower dependence on mesh quality and orientation. Further new features are the extension to unstructured grids and support for arbitrary permeability tensors (allowing for both scalar heterogeneity, and shear anisotropy). The most important advancement in this work is the self-consistent modeling of two-phase multi-component Fickian diffusion. We present several numerical examples to illustrate the powerful features of our combined MHFE–dg method with respect to lower-order calculations, ranging from simple two component fluids to more challenging real problems regarding CO₂ injection into a vertical domain saturated with a multi-component petroleum fluid.     

Application of wavelet-artificial intelligence hybrid models for water quality prediction: a case study in Aji-Chay River,Iran

The accuracy of Artificial Neural Network (ANN), Adaptive Neuro-Fuzzy Inference System (ANFIS), wavelet-ANN and wavelet-ANFIS in predicting monthly water salinity levels of northwest Iran’s Aji-Chay River was assessed. The models were calibrated, validated and tested using different subsets of monthly records (October 1983 to September 2011) of individual solute (Ca²⁺, Mg²⁺, Na⁺, SO₄ ^2? and Cl^?) concentrations (input parameters, meq L^?1), and electrical conductivity-based salinity levels (output parameter, µS cm^?1), collected by the East Azarbaijan regional water authority. Based on the statistical criteria of coefficient of determination (R²), normalized root mean square error (NRMSE), Nash–Sutcliffe efficiency coefficient (NSC) and threshold statistics (TS) the ANFIS model was found to outperform the ANN model. To develop coupled wavelet-AI models, the original observed data series was decomposed into sub-time series using Daubechies, Symlet or Haar mother wavelets of different lengths (order), each implemented at three levels. To predict salinity input parameter series were used as input variables in different wavelet order/level-AI model combinations. Hybrid wavelet-ANFIS (R² = 0.9967, NRMSE = 2.9 × 10^?5 and NSC = 0.9951) and wavelet-ANN (R² = 0.996, NRMSE = 3.77 × 10^?5 and NSC = 0.9946) models implementing the db4 mother wavelet decomposition outperformed the ANFIS (R² = 0.9954, NRMSE = 3.77 × 10^?5 and NSC = 0.9914) and ANN (R² = 0.9936, NRMSE = 3.99 × 10^?5 and NSC = 0.9903) models.     

Hydrogeologic framework of the Maku area basalts, northwestern Iran

The Maku area in northwestern Iran is characterized by young lava flows which erupted from Mount Ararat in Turkey. These fractured volcanic rocks overlie alluvium associated with pre-existing rivers and form a good basalt-alluvium aquifer over an area of 650 km². Groundwater discharge occurs from 12 large springs, ranging from 20 to 4,000 L s^?1, and from some extraction wells. Permian and Oligo-Miocene age limestones along the northern boundary of the Bazargan and Poldasht Plains basalts are intensively karstified and groundwater from these high lands easily enters the basalt-alluvium aquifers. The transmissivity of the basalt-alluvium aquifer ranges from 24 to 870 m² d^?1, indicating heterogeneity. Groundwater of the aquifer is a sodium-bicarbonate and mixed cation-bicarbonate type and the concentration of fluoride is higher than the universal maximum admissible concentrations for drinking. In order to determine the chemical composition and identify the source of the high fluoride concentrations in the groundwater of the basaltic area, water samples from the springs, wells and rivers were analyzed. The results indicate that the high fluoride water enters the study area from the Sari Su River.