A robust mesh optimisation method for multiphase porous media flows

Flows of multiple fluid phases are common in many subsurface reservoirs. Numerical simulation of these flows can be challenging and computationally expensive. Dynamic adaptive mesh optimisation and related approaches, such as adaptive grid refinement can increase solution accuracy at reduced computational cost. However, in models or parts of the model domain, where the local Courant number is large, the solution may propagate beyond the region in which the mesh is refined, resulting in reduced solution accuracy, which can never be recovered. A methodology is presented here to modify the mesh within the non-linear solver. The method allows efficient application of dynamic mesh adaptivity techniques even with high Courant numbers. These high Courant numbers may not be desired but a consequence of the heterogeneity of the domain. Therefore, the method presented can be considered as a more robust and accurate version of the standard dynamic mesh adaptivity techniques.     

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.     

Short-term greenhouse emission lowering effect of biochars from solid organic municipal wastes

Qatar economy has been growing rapidly during the last two decades during which waste generation and greenhouse gas emissions increased exponentially making them among the main environmental challenges facing the country. Production of biochar from municipal solid organic wastes (SOWs) for soil application may offer a sustainable waste management strategy while improving crop productivity and sequestering carbon. This study was conducted to (1) investigate the physicochemical parameters of biochars for SOW, (2) select the best-performing biochars for soil fertility, and (3) evaluate the potential benefits of these biochars in lowering greenhouse gases (GHGs) during soil incubation. Biochars were produced from SOW at pyrolysis temperatures of 300–750 °C and residence times of 2–6 h. Biochars were characterized before use in soil incubation to select the best-performing treatment and evaluation of potential GHG-lowering effect using CO₂ emission as proxy. Here, soil–biochar mixtures (0–2%w/w) were incubated in greenhouse settings for 120 days at 10% soil moisture. Soil properties, such as pH, EC, TC, and WHC, were significantly improved after soil amendment with biochar. Two biochars produced from mixed materials at 300–500 °C for 2 h and used at 0.5–1% application rate performed the best in enhancing soil fertility parameters. A significant decrease in CO₂ emission was observed in vials with soil–biochar mixtures, especially for biochars produced at 500 °C compared the corresponding raw materials which exhibited an exponential increase in the CO₂ emission. Hence, application of biochar to agricultural soils could be beneficial for simultaneously improving soil fertility/crop productivity while sequestering carbon, thereby reducing anthropogenic emissions of GHGs.     

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.     

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.     

2D seismic reflection contribution to structural and geometric study of Cenozoic aquifer systems in the Sisseb El-Alem basin,central-eastern Tunisia

The Sisseb-ElAlem basin in central-eastern Tunisia contains several aquifers intensely used in agriculture. Highly exploited shallow aquifers are vulnerable to pollution and drought and reveal high salinity. Consequently, the new hydrogeological surveys were oriented to deep aquifers. The present work corresponds to study of 2D seismic sections, calibrated by water and petroleum wells, and geological outcrops. It highlights geometry and faulting of buried and preserved structures in depth and leads to identify the spatial distribution of Cenozoic potential reservoir layers in the different paleogeographic domains. The Sisseb-ElAlem basin is affected mainly by near E- and N-trending faults. Iso-depth and iso-thickness mapping were used to follow lateral and vertical variations of sedimentary series related to tectonics. The study area consists of Sisseb, Sbikha and El Ktifa uplifts, and El Alem, and Bled Saadia subsiding synclines. The Oligocene-Miocene reservoir series are exploited locally at the high borders. They reveal thickening and deepening southward in the El Alem-Bled Saadia domain, where they can be exploited at about 400 to 2000 m depths.This study should contribute to the orientation of water research to deeper aquifers, and therefore, ensure sufficient drinking and irrigation water supply mainly in areas called for sustainable development.     

Investigation of monotonic and cyclic behavior of sand using a bounding surface plasticity model

Developing the pore water pressures in loose to medium sands below the water table may lead to liquefaction during earthquakes. The simulation of liquefaction (cyclic mobility and flow liquefaction) in sandy soils is one of the major challenges in constitutive modeling of soils. This paper presents the simulation of sand behavior using a critical state bounding surface plasticity model (Dafalias and Manzari’s model, 2004) during monotonic and cyclic loading. The drained, undrained, and cyclic triaxial tests were simulated using Dafalias and Manzari’s model. The simulation results showed that the model predicts behavior of sand, reasonably well. Also, for CSR?<?0.2, number of cycles for liquefaction is significantly increased. The residual strength of Babolsar sand is produced when it is deformed to an axial strain of 20 to 25%.     

Characterization of aerosol type based on aerosol optical properties over Baghdad,Iraq

Aerosol optical depth (AOD), Angstrom exponent (AE), and ozone monitoring instrument aerosols index (OMI-AI) data, derived from MODerate Resolution Imaging Spectroradiometer (MODIS) and OMI sensor on board NASA’s Aqua satellite and NASA-Aura satellite platforms, have been analyzed and classified over Baghdad, Iraq, for an 8-year period (2008–2015). In order to give an obvious understanding of temporal inconsistency in the characteristics and classification of aerosols during each season separately, PREDE POM-02 sky radiometer measurements of AOD, carried out during a 2-year period (2014–2015), were compared with MODIS–Aqua AODs. On seasonal bases, MODIS–Aqua AODs corroborate well with ground-based measurements, with correlation coefficients ranging between 0.74 and 0.8 and RMSE ranging from 0.097 to 0.062 during spring and autumn seasons respectively. The overall satellite- and ground-based measurement comparisons showed a good agreement with correlation coefficients of 0.78 and RMSE of 0.066. These results suggest that MODIS–Aqua gives a good estimate of AOD. Analysis of MODIS–Aqua data for the 8-year period showed that the overall mean AOD, AE, and OMI-AI over Baghdad were 0.44?±?0.16, 0.77?±?0.29, and 1.34?±?0.33 respectively. AOD records presented a unique peak which was extended from mid-spring (April) to mid-summer (July) while the AE annual variability indicated a more complicated behavior with minimum values during the period from late spring (May) to early autumn (September). The maximum AOD and OMI-AI values occurred during summer while their minimum values occurred during winter. The AE showed an opposite behavior to that of AOD such that the highest AE values occurred during autumn and winter and the lowest values happened during spring and summer. This behavior may be attributed to the domination of coarse aerosol particles during autumn and winter seasons and fine aerosol particles during spring and summer seasons. A Hybrid Single-Particle Lagrangian Integrated Trajectory model was utilized to determine the source of air mass transport and to recognize the variability of aerosol origin regions. Finally, AOD, AE, and OMI-AI values have been employed to identify several aerosol types and to present seasonal heterogeneity in their contribution based on their origins.