Application of a 3D photorealistic model for the geological analysis of the Permian carbonates (Khuff Formation) in Saudi Arabia

Three dimensional (3D) photorealistic models of geological outcrops have the potential to enhance the teaching of earth sciences by providing scale models in a virtual reality environment. These models can be run on low-cost desktop computers. Photorealistic models for geological outcrops are a digital illustration of outcrop photographs with either a point cloud representation or Triangular Irregular Network (TIN) mesh of the outcrop surface. The level of detail for these models is dependent on the target resolutions (physical and optical) that were used during data acquisition. In addition, the technique in which the data is rendered as a digital model affects the level of detail that can be observed by the geologists. A colored point cloud representation is suitable for large-scale features, but fine details are lost when the geologist zooms in to view the model close up. In contrast, a photorealistic model that is constructed from photographs draped onto a triangle mesh surface derived from Light Detection and Ranging (LiDAR) point clouds provides a level of detail that is restricted only by the resolution of the photographs.     

3D geomechanical modeling and estimating the compaction and subsidence of Fahlian reservoir formation (X-field in SW of Iran)

A geomechanical model can reveal the mechanical behavior of rocks and be used to manage the reservoir programs in a better mode. Fluid pressure will be reduced during hydrocarbon production from a reservoir. This reduction of pressure will increase the effective stress due to overburden sediments and will cause porous media compaction and surface subsidence. In some oil fields, the compacting reservoir can support oil and gas production. However, the phenomena can also cause the loss of wells and reduced production and also cause irreparable damage to the surface structures and affect the surrounding environment. For a detailed study of the geomechanical behavior of a hydrocarbon field, a 3D numerical model to describe the reservoir geomechanical characteristics is essential. During this study, using available data and information, a coupled fluid flow-geomechanic model of Fahlian reservoir formation in X-field in SW of Iran was constructed to estimate the amount of land subsidence. According to the prepared model, in this field, the maximum amount of the vertical stress is 110 MPa and the maximum amount of the horizontal stress is 94 MPa. At last, this model is used for the prediction of reservoir compaction and subsidence of the surface. The maximum value of estimated ground subsidence in the study equals to 29 mm. It is considered that according to the obtained values of horizontal and vertical movement in the wall of different wells, those movements are not problematic for casing and well production and also the surrounding environment.     

Estimating geometrical parameters of cylindrical targets detected by ground-penetrating radar using template matching algorithm

In the current research, the ground-penetrating radar (GPR) method has been employed to identify physical and geometrical parameters of buried cylindrical structures using the pattern recognition approach. To achieve this goal, the well-established mathematical relationships between geometrical parameters of cylindrical target (radius, burial depth, and horizontal location) and the associated GPR hyperbolic response characteristics are employed using the template matching method. In order to validate the applicability of the template matching method in providing estimates of such parameters, the method is first examined on GPR responses of synthetic models with known geometrical parameters followed by applying on real data using two different similarity criteria including 2-D spatial convolution and normalized cross correlation in the wave number domain. In the first step, the GPR responses of 71 synthetic models encompassing one, two, and three horizontal cylinders were produced using the improved 2-D finite difference in frequency domain. Then, appropriate preprocessing sequences to reduce random noise caused by forward modeling were applied on synthetic data. The proposed algorithm applied on several synthetic model responses could estimate the known geometrical parameters of the buried cylinders with acceptable accuracy (maximum error of 15%). The template matching algorithm was also used to extract geometrical parameters of water and wastewater pipes buried in Imam Hossein Square, Isfahan city, as real GPR data. Depending on environmental conditions and subsurface host formation, the real GPR data normally contain a variety of noises; therefore, a series of appropriate objective preprocessing and processing stages were designed in order to apply on real GPR images before deploying template matching algorithm. The applicability of the template matching algorithm on real data and validity of the estimated parameters were proved based on assessing the accuracy of the estimated geometrical parameters of respective pipes through GPR response versus the measured parameters. The proposed algorithm was designed in such a way that all steps of estimating geometrical parameters of buried cylindrical targets are automatically carried out.     

Impact of rock-water interactions and recharge on water resources quality of the Agadir-Essaouira basin,southwestern Morocco

The Agadir-Essaouira area in the occidental High Atlas Mountains of Morocco is characterized by a semi-arid climate. The scarcity and quality of water resources, exacerbated by long drought periods, constitute a major problem for a sustainable development of this region. Groundwater resources of carbonate units within Jurassic and Cretaceous aquifers are requested for drinking and irrigation purposes. In this study, we collected 84 samples from wells, boreholes, springs, and rivers. Hydrochemical and isotopic data were used to examine the mineralization and origin of water, which control groundwater quality. The chemical composition of water seems to be controlled by water-rock interactions, such as dissolution of carbonates (calcite and dolomite), weathering of gypsum, as well as ion exchange processes, which explain the observed variability. Stable isotopes results show that groundwater from the mainly marly Cretaceous aquifer are submitted to an evaporation effect, while samples from the chiefly calcareous Jurassic aquifer indicate a meteoric origin, due to a rapid infiltration of recharge runoff through the karstic outcrops. The low values of δ¹⁸O and δ²H suggest a local recharge from areas with elevations ranging from 400 to 1200 m for the Cretaceous aquifer and from 800 to 1500 m for the Jurassic units.     

Isotope evidence for quantifying river evaporation and recharge processes in the lower reaches of the Yellow River

Evaporation and recharge are important hydrological processes in the water cycle. However, accurately quantifying these two processes of river remains to be difficult due to their spatial heterogeneity and the limitations of traditional methods. In this study, a more reliable method of stable isotopes of δ ¹⁸O and δ ²H based on the Rayleigh distillation equation and mass conservation was used to estimate the evaporation and recharge of the rivers in the lower reaches of the Yellow River, North China Plain. Comprehensive sampling campaigns including 30 surface water samples from 10 rivers, 33 groundwater samples from domestic and observation wells, and two Yellow River water samples were conducted. The results showed that the evaporation proportion of the rivers based on δ ¹⁸O and δ ²H both averaged 14.4%. The evaporation proportions in each river did not completely follow a linear increasing trend along the flow path. This phenomenon could be mainly explained by the different proportions of recharge from groundwater and Yellow River water. With closer to the Yellow river, evaporation of the rivers decreased while the recharge by the Yellow River increased. Regression equations based on δ ¹⁸O, δ ²H, and their average revealed that the evaporation proportion respectively increased by 1.02, 0.79, and 0.90% with the increase in the distance to the Yellow River per 10 km. On the contrary, the recharge proportion decreased by 7.68, 5.51, and 6.59%, respectively. In addition, using δ ¹⁸O rather than δ ²H was more reliable in studying the spatial influence of the Yellow River on evaporation and recharge. Sensitivity analysis showed that the evaporation model was most sensitive to isotopic composition, rather than to air temperature or relative humidity. The results of this study provide insights into the determination of river hydrological processes and the management of water resources.     

Comparing data-driven landslide susceptibility models based on participatory landslide inventory mapping in Purwosari area,Yogyakarta, Java

There are different approaches and techniques for landslide susceptibility mapping. However, no agreement has been reached in both the procedure and the use of specific controlling factors employed in the landslide susceptibility mapping. Each model has its own assumption, and the result may differ from place to place. Different landslide controlling factors and the completeness of landslide inventory may also affect the different result. Incomplete landslide inventory may produce significance error in the interpretation of the relationship between landslide and controlling factor. Comparing landslide susceptibility models using complete inventory is essential in order to identify the most realistic landslide susceptibility approach applied typically in the tropical region Indonesia. Purwosari area, Java, which has total 182 landslides occurred from 1979 to 2011, was selected as study area to evaluate three data-driven landslide susceptibility models, i.e., weight of evidence, logistic regression, and artificial neural network. Landslide in the study area is usually affected by rainfall and anthropogenic activities. The landslide typology consists of shallow translational and rotational slide. The elevation, slope, aspect, plan curvature, profile curvature, stream power index, topographic wetness index, distance to river, land use, and distance to road were selected as landslide controlling factors for the analysis. Considering the accuracy and the precision evaluations, the weight of evidence represents considerably the most realistic prediction capacities (79%) when comparing with the logistic regression (72%) and artificial neural network (71%). The linear model shows more powerful result than the nonlinear models because it fits to the area where complete landslide inventory is available, the landscape is not varied, and the occurence of landslide is evenly distributed to the class of controlling factor.     

The diagenetic geochemistry and contamination assessment of iron,cadmium, and lead in the sediments from the Shuangtaizi estuary,China

Sediment and pore water samples have been collected from the coastal tidal flat in the Shuangtaizi estuary, China, in order to investigate the geochemical behavior of iron, cadmium, and lead during diagenesis and to assess the degree of contamination. The calculated enrichment factors and geoaccumulation indices for separate elements show that anthropogenic activities have had no significant influence on the distribution of Fe and Pb in the study area, whereas the distribution of Cd has been closely influenced in this way. The high percentage of exchangeable Cd (average of 56.34%) suggests that Cd represents a potential hazard to benthic organisms in the estuary. The calculated diffusive fluxes of metals show that the most mobilized metal is Fe (9.22 mg m^?2 a^?1), followed by Cd (0.54 mg m^?2 a^?1) and Pb (0.42 mg m^?2 a^?1). Low Fe²⁺ contents in surface pore water, alongside high chromium-reducible sulfur contents, and low acid-volatile sulfur, and elemental sulfur contents at 0–25 cm depth in sediments show that Fe²⁺ is formed by the reduction of Fe oxides and is transformed first to a solid phase of iron monosulfides (FeS) and eventually to pyrite (FeS₂). The release of adsorbed Pb due to reductive dissolution of Fe/Mn oxides during early diagenesis could be a source of Pb²⁺ in pore water. From the relatively low total organic carbon contents measured in sediments (0.46–1.28%, with an average of 0.94%) and the vertical variation of Cd²⁺ in pore water, sulfide or Fe/Mn oxides (instead of organic matter) are presumed to exert a significant influence on carrying or releasing Cd by the sediments.     

Observation of abnormal thermal and infrasound signals prior to the earthquakes: a study on Bonin Island earthquake M7.8 (May 30, 2015)

Scientists from all over the world try to incorporate multi-disciplinary precursors to forecast the earthquake on a short-term basis. The authors here have analyzed outgoing longwave radiations acquired from polar-orbiting National Oceanic and Atmospheric Adminstration (NOAA) satellites and long-period infrasound waves recorded by the ground observatories in China prior to the recent Bonin Islands, Japan region earthquake which occurred on May 30, 2015 with the magnitude of 7.8. The anomalous outgoing longwave radiation (OLR) was observed on May 15, 2015, and was recorded by “NOAA 18” satellite during its “night pass.” Similarly, an abnormal infrasound spike was recorded at the Beijing station on May 17, 2015. The delay in observing anomalous infrasound waves compared to the OLR anomaly is due to these low-frequency waves traveling at low speed with the velocity range of 10–15 m/s. From the analysis of the results, it can be inferred that there is a substantial relation between parameters like OLR and infrasound waves; hence, the authors conclude that it is possible to forecast the earthquake on a short-term basis with reasonable accuracy.     

Recovery of groundwater in the Sanriku region contaminated by the tsunami inundation from the 2011 Tohoku earthquake

Many wells in the Sanriku region used as sources for water supply systems were heavily contaminated by the tsunami of the 2011 great Tohoku earthquake on March 11 in 2011. To better understand the nature of the groundwater contamination by the tsunami inundation and to clarify the recovery process of contaminated groundwater at the study wells, groundwater monitoring has been conducted once or twice yearly since early summer in 2011. High and abnormal values of electric conductivity (EC), chloride ion concentration (CIC), Na⁺, Ca⁺, heavy metal ions, and heavier isotopes of the contaminated groundwater were also obtained in April and June 2011. The chemical elements have rapidly and exponentially decreased as a result of effective pumping of the contaminated groundwater from the study wells and because of abundant rainfall in 2011. In April 2015 (about 4 years after the tsunami inundation), the CIC and EC of the contaminated groundwater of two study wells in Minamisanriku town had reached pre-inundation values. The estimated residence times of groundwater of the two study wells were 105–118 days in the full-day pumping stage and 910–1000 days in the daytime-only pumping stage.     

The human impact on the transformation of juniper forest landscape in the western part of the Pamir-Alay range (Tajikistan)

Detailed analyses were conducted of human impact on juniper forest landscapes occurring within the Zarafshan Range (Pamir-Alay). Juniperus seravschanica and J. semiglobosa belong to forest-forming species in Central Asia. At present, juniper forests all over Tajikistan are seriously threatened as a result of excessive logging and cattle grazing. The aim of this paper is to present juniper forest transformation as a result of human activities as well as the diversity of soil properties in the organic and humus horizons in the altitudinal system of soil zonation. Three groups of phytocoenoses were distinguished: those with a dominant share of Juniperus seravschanica; those with a dominant share of J. semiglobosa; and mixed. Associations with Juniperus seravschanica and J. semiglobosa feature several variants of phytocoenoses with dominant species: Artemisia lehmanniana, A. dracunculus, Eremurus olgae, Festuca sulcata, Ligularia thomsonii, Stipa turkestanica, Thymus seravschanicus, and Ziziphora pamiroalaica. The collected soil samples differ in their granulometric composition. Gravelly cobble fractions >2 mm are dominant; the share of sandy particles <2 mm is much lower (about 10–20%). Fraction 0.5–0.05 attains 35% on average. The Corg content of the soil varied from 0.26 to 11.40% in the humus horizon (A) and from 4.3 to 25% in the organic (O). Similar relationships were reported in the case of Ntot concentration. A clear relationship can be observed between concentrations of Corg and Ntot. Soil pH varied, ranging from very low acidic (pH 5.5) to neutral (pH 8.5). The content of available P varied; high concentrations were noted in organic (O) (40.46–211 mg kg^?1) and mixed horizons (OA) (2.61–119 mg kg^?1). Maximum accumulations of Pavail (1739.6 mg kg^?1) and Ptot (9696 mg kg^?1) were observed at a site heavily affected by intense grazing. Concentrations of Mgavail varied from 116 to 964 mg kg^?1. Most of the analysed soil profiles lacked an organic horizon; only thin humus occurred.