A method of alternating characteristics with application to advection-dominated environmental systems

We present a numerical method for solving a class of systems of partial differential equations (PDEs) that arises in modeling environmental processes undergoing advection and biogeochemical reactions. The salient feature of these PDEs is that all partial derivatives appear in linear expressions. As a result, the system can be viewed as a set of ordinary differential equations (ODEs), albeit each one along a different characteristic. The method then consists of alternating between equations and integrating each one step-wise along its own characteristic, thus creating a customized grid on which solutions are computed. Since the solutions of such PDEs are generally smoother along their characteristics, the method offers the potential of using larger time steps while maintaining accuracy and reducing numerical dispersion. The advantages in efficiency and accuracy of the proposed method are demonstrated in two illustrative examples that simulate depth-resolved reactive transport and soil carbon cycling.     

Seismic stratigraphy of the Mianwali and Bannu depressions,north-western Indus foreland basin

International Journal of Earth Sciences - Regional seismic reflection profiles, deep exploratory wells, and outcrop data have been used to study the structure and stratigraphic architecture of the...     

Numerical solution and convergence analysis of steam injection in heavy oil reservoirs

In this paper, the numerical methods for solving the problem of steam injection in the heavy oil reservoirs are presented. We consider a 3-dimensional model of 3-phase flow, oil, water, and steam, with the effect of 3-phase relative permeability. Interphase mass transfer of water and steam is considered; oil is assumed nonvolatile. We apply the simultaneous solution approach to solve the corresponding nonlinear discretized partial differential equation in the fully implicit form. The convergence of finite difference scheme is proved by the Rosinger theorem. The heuristic Jacobian-Free-Newton-Krylov (HJFNK) method is proposed for solving the system of algebraic equations. The result of this proposed numerical method is well compared with some experimental results. Our numerical results show that the first iteration of the full approximation scheme (FAS) provides a good initial guess for the Newton method. Therefore, we propose a new hybrid-FAS-HJFNK method while there is no steam in the reservoir. The numerical results show that the hybrid-FAS-HJFNK method converges faster than the HJFNK method.     

Subduction-related mafic to felsic magmatism in the Malayer–Boroujerd plutonic complex,western Iran

The Malayer–Boroujerd plutonic complex (MBPC) in western Iran, consists of a portion of a magmatic arc built by the northeast verging subduction of the Neo-Tethys plate beneath the Central Iranian Microcontinent (CIMC). Middle Jurassic-aged felsic magmatic activity in MBPC is manifested by I-type and S-type granites. The mafic rocks include gabbroic intrusions and dykes and intermediate rocks are dioritic dykes and minor intrusions, as well as mafic microgranular enclaves (MMEs). MBPC Jurassic-aged rocks exhibit arc-like geochemical signatures, as they are LILE- and LREE-enriched and HFSE- and HREE-depleted and display negative Nb–Ta anomalies. The gabbro dykes and intrusions originated from metasomatically enriched garnet-spinel lherzolite [Degree of melting (f_mel) ~ 15%] and exhibit negative Nd and positive to slightly negative ε_Hf(T) (+ 3.0 to ? 1.6). The data reveal that evolution of Middle Jurassic magmatism occurred in two stages: (1) deep mantle-crust interplay zone and (2) the shallow level upper crustal magma chamber. The geochemical and isotopic data, as well as trace element modeling, indicate the parent magma for the MBPC S-type granites are products of upper crustal greywacke (f_mel: 0.2), while I-type granites formed by partial melting of amphibolitic lower crust (f_mel: 0.25) and mixing with upper crustal greywacke melt in a shallow level magma chamber [Degree of mixing (f_mix): 0.3]. Mixing between andesitic melt leaving behind a refractory dense cumulates during partial crystallization of mantle-derived magma and lower crustal partial melt most likely produced MMEs (f_mix: 0.2). However, enriched and moderately variable ε_Nd(T) (? 3.21 to ? 4.33) and high (⁸⁷Sr/⁸⁶Sr)_i (0.7085–0.7092) in dioritic intrusions indicate that these magmas are likely experienced assimilation of upper crustal materials. The interpretations of magmatic activity in the MBPC is consistent with the role considered for mantle-derived magma as heat and mass supplier for initiation and evolution of magmatism in continental arc setting, elsewhere.     

Health risk of heavy metals from vegetables irrigated with sewage water in peri-urban of Dera Ismail Khan,Pakistan

Health risks of heavy metals in vegetables irrigated with sewage water were investigated in the present study. The findings indicated a massive accumulation of heavy metals in soil and vegetables collected from Dera Ismail khan, Pakistan. The concentration levels of heavy metal in vegetables grown on soil irrigated with untreated sewage water were significantly higher at (P ≤ 0.001) than in vegetables grown on fresh-water-irrigated soil and proceeded the recommended limits of World health organization. Moreover, the findings also indicated that the adults and children consuming such vegetables ingested a large proportion of the selected toxic metals. Health risk index was greater than one for Pb and Cd in all the selected vegetables and was greater than one for Ni in three vegetables like Spinacia oleracea, Benincasa fistulosa and Lactuca sativa. Health risk assessment would be a useful tool for information regarding any threats of heavy metals contamination in vegetables.     

Impact of drilling mud on chemistry and microbiology of an Upper Triassic groundwater after drilling and testing an exploration well for aquifer thermal energy storage in Berlin (Germany)

After completion of an exploration well, sandstones of the Exter Formation were hydraulically tested to determine the hydraulic properties and to evaluate chemical and microbial processes caused by drilling and water production. The aim was to determine the suitability of the formation as a reservoir for aquifer thermal energy storage. The tests revealed a hydraulic conductivity of 1–2 E-5 m/s of the reservoir, resulting in a productivity index of 0.6–1 m³/h/bar. A hydraulic connection of the Exter Formation to the overlaying, artesian “Rupelbasissand” cannot be excluded. Water samples were collected for chemical and microbiological analyses. The water was similarly composed as sea water with a maximum salinity of 24.9 g/L, dominated by NaCl (15.6 g/L Cl and 7.8 g/L Na). Until the end of the tests, the water was affected by drilling mud as indicated by the high pH (8.9) and high bicarbonate concentration (359 mg/L) that both resulted from the impact of sodium carbonate (Na₂CO₃) additives. The high amount of dissolved organic matter (>?58 mg/L) and its molecular-weight distribution pattern indicated that residues of cellulose, an ingredient of the drilling mud, were still present at the end of the tests. Clear evidence of this contamination gave the measured uranine that was added as a tracer into the drilling mud. During fluid production, the microbial community structure and abundance changed and correlated with the content of drilling mud. Eight taxa of sulfate-reducing bacteria, key organisms in processes like bio-corrosion and bio-clogging, were identified. It can be assumed that their activity will be affected during usage of the reservoir.     

Geospatial assessment of soil erosion intensity and sediment yield: a case study of Potohar Region,Pakistan

Estimation of spatial extent of soil erosion, one of the most serious forms of land degradation, is critical because soil erosion has serious implications on soil fertility, water ecosystem, crop productivity and landscape beauty. The primary objective of the current study was to assess and map the soil erosion intensity and sedimentation yield of Potohar region of Pakistan. Potohar is the rainfed region with truncated and complex topography lying at the top of the Indus Basin, the world’s largest irrigation networks of canals and barrages. Spatially explicit Revised Universal Soil Loss Equation (RUSLE) Model integrated with Remote Sensing-GIS techniques was used for detecting/mapping of erosion prone areas and quantification of soil losses. The results show that the Potohar region is highly susceptible to soil erosion with an average annual soil loss of 19 tons ha^?1 year^?1 of which the maximum erosion (70–208 tons ha^?1 year^?1) was near the river channels and hilly areas. The sediment yield due to the erosion is as high as 148 tons ha^?1 year^?1 with an average of 4.3 tons ha^?1 year^?1. It was found that 2.06% of the total area falls under severe soil erosion, 13.34% under high erosion, 15.35% under moderate soil erosion while 69.25% of the area lies in the low (tolerable) soil erosion. Chakwal and Jhelum districts of the region are seriously affected by erosion owing to their topography and soil properties. The information generated in this study is a step forward towards proper planning and implementation of strategies to control the erosion and for protection of natural resources. It is, hence, necessary that suitable water harvesting structures be made to control water to prevent soil erosion and provision of water in the lean season in this region. Tree plantation and other erosion control practices such as strip cropping can also minimize soil erosion in this region.     

Development of an efficient surrogate model based on aquifer dimensions to prevent seawater intrusion in anisotropic coastal aquifers,case study: the Qom aquifer in Iran

Coastal aquifers are usually exposed to saltwater intrusion. Therefore, groundwater extracted from these aquifers should be regulated considering their dimensions and effective parameters. In this paper, optimum discharge from a large number of exploitation wells is evaluated according to variations of width, length, and anisotropy coefficient in the Qom aquifer near the salt lake in central Iran. First, the wells are divided into clusters to decrease the number of decision variables. Then, the location and discharge from each cluster is obtained using SEAWAT and charged system search (CSS) simulation–optimization model with the assumption of three-dimensional variable density flow. The maximum discharge considering various anisotropy rates is computed based on different values of lengths and widths of the aquifer. Finally, an M5-tree model is trained using the obtained samples to derive a linear relationship between input and output data. Based on the results, for various ranges of width and length of an aquifer with impermeable boundaries, different linear equations for optimum discharge are obtained. Also, it was found that for an aquifer with a small width, the critical discharge is a function of the length while the effect of the boundaries is negligible. Sensitivity analysis of the anisotropy coefficient reveals that with increasing the anisotropy rate, thickness and slope of the transition zone decrease and as the maximum discharge increases consequently. However, the sensitivity of the discharge to anisotropy rate is not remarkable. A comparison between the results of this study with those of the analytical method based on sharp interface assumption is carried out. For the critical condition, the best agreement between analytical equation (\(\overline {L} =0.87\overline {W} +0.62\)) and proposed method (\(\overline {L} =0.83\overline {W} - 1.41\)) is achieved for the anisotropic aquifer when the 50% isochlor is assumed as the measure of salt water intrusion.     

Ghaleh-khargushi rhyodacite and Gorid andesite from Iran: characterization,uses, and durability

Durability of building stones is an important issue in sustainable development. Crystallization of soluble salts is recognized as one of the most destructive weathering agents of building stones. For this reason, durability of Ghaleh-khargushi rhyodacite and Gorid andesite from Iran was investigated against sodium sulfate crystallization aging test. Petrographic and physico-mechanical properties and pore size distribution of these stones were examined before and after the aging test. The characteristics of the microcracks were quantified with fluorescence-impregnated thin sections. Durability and physico-mechanical characteristics of Ghaleh-khargushi rhyodacite are mainly influenced by preferentially oriented preexisting microcracks. Stress induced by salt crystallization led to the widening of preexisting microcracks in Ghaleh-khargushi rhyodacite, as confirmed by the pore size distributions before and after the aging test. The preexisting microcracks of Gorid andesite were attributed to the mechanical stress induced by contraction of lava during cooling. The number of transcrystalline microcracks was significantly increased after the aging test. The degree of plagioclase microcracking was proportional to its size. Durability of the studied stones depends on initial physico-mechanical properties, pore size distribution, and orientation of microcracks. Initial effective porosity is found to be a good indicator of the stones’ durability. Salt crystallization resulted in an increase in the effective porosity with a parallel decrease in the wave velocities. Surface microroughness parameters increased with the development of salt crystallization-induced microcracking. Gorid andesite showed higher quality and durability than Ghaleh-khargushi rhyodacite.     

Hydrothermal alteration and evolution of Zr-Th-U-REE mineralization in the microgranite of Wadi Ras Abda,North Eastern Desert,Egypt

Ras Abda plutonic suite, North Eastern Desert of Egypt, consists predominantly of Neoproterozoic calc-alkaline older granites. Minor exposures of pink microgranite are occurring along Wadi Ras Abda within the older granites. Previous studies on this area demonstrated that the microgranite is altered in some parts and contains anomalous concentrations of rare metal elements (Zr, Th, and U). These altered and mineralized zones are re-assessed using field observations, chemical analysis, and by the application of various transmitted light and electron microscopic techniques. The rare metals exist as mineral segregation grew freely into open cavities of the microgranite and concordant with the NNE strike-slip fault movement. The mineralized zones contain an assemblage of secondary magnetite, zircon, uranothorite, columbite-(Mn), fergusonite-(Y), and allanite-(Ce). The extreme abundance of zircon in the mineralized zone, along with other evidence, indicates a hydrothermal origin of this zircon together with associated rare metals. The geochemical investigation and mass balance calculations revealed extreme enrichment of Zr, Th, U, Y, Nb, Ta, and REE. Post-magmatic hydrothermal alterations resulted in such pronounced chemical and mineralogical heterogeneity. The hydrothermal fluids are thought to be oxidizing, alkaline and of medium temperature (>?250 °C). The average contents of the elements Zr (1606 ppm), Th (1639 ppm), U (306 ppm), Nb (955 ppm), and REE (1710 ppm) in the mineralized microgranite reach sub-economic levels and could be a potential source of these elements.