Image transforms for determining fit-for-purpose complexity of geostatistical models in flow modeling

Increased diversity of water or energy resources has led to an increased complexity in models aimed at representing accurately dynamic behavior and geological variability in such systems. In terms of variability of properties at least, simple layered models have mostly been replaced with more complex geostatistical models. The newest trend is to replace covariance-based models with geologically more realistic models such as Boolean, multiple-point, surface-, or process-based models. In this paper, we address the following question: given some design purpose or a set of flow-based decision variables, does adding more complexity increase predictability and ultimately improve decisions based on such models? In this paper, we do not attempt to make any generalizing statements or answer this question with yes/no, but provide some initial ideas on practical workflows to discover the needed complexity. We do treat complexity only in the context of decision making under uncertainty. Two workflows are proposed: complexifying versus simplifying. In these workflows, we attempt to extract, using image transforms, relevant features of the variability between geostatistical realizations that are related to uncertainty in flow dynamics. A simple distance-based calibration between the static variability and dynamic variability provides a means to decide on what the relevant complexity of geostatistical models should be for the given purpose.     

Electrostatic solitons in rotating dusty plasmas with anisotropic ion pressure

The Zakharov-Kuznetsov (ZK) equation is derived for electrostatic wave in a rotating magnetoplasma with anisotropic ion pressure and in the presence of stationary charged dust particles. The anisotropic ion pressure is defined using double adiabatic Chew-Golberger-Low (CGL) theory. The reductive perturbation method is employed to study the dynamics of obliquely propagating low frequency ion acoustic wave with adiabatic ions. It is found that the ion pressure anisotropy, polarity, density of the dust particles and rotational frequency have significant effects on the formation nonlinear structures in rotating magnetized dusty plasmas. The numerical results are also presented for illustration.     

GIS-based DRASTIC model for groundwater vulnerability and pollution risk assessment in the Peshawar District,Pakistan

Groundwater is the most economic natural source of drinking in urban and rural areas which are degraded due to high population growth and increased industrial development. We applied a GIS-based DRASTIC model in a populated urban area of Pakistan (Peshawar) to assess groundwater vulnerability to pollution. Six input parameters—depth to phreatic/groundwater level, groundwater recharge, aquifer material, soil type, slope, and hydraulic conductivity—were used in the model to generate the groundwater vulnerable zones. Each parameter was divided into different ranges or media types, and ratings R?=?1?–?10 were assigned to each factor where 1 represented the very low impact on pollution potential and 10 represented very high impact. Weight multipliers W?=?1?–?5 were also used to balance and enhance the importance of each factor. The DRASTIC model scores obtained varied from 47 to 147, which were divided into three different zones: low, moderate, and high vulnerability to pollution. The final results indicate that about 31.22, 39.50, and 29.27% of the total area are under low, moderate, and high vulnerable zones, respectively. Our method presents a very simple and robust way to assess groundwater vulnerability to pollution and helps the decision-makers to select appropriate landfill sites for waste disposals, and manage groundwater pollution problems efficiently.     

Effect of increasing temperatures on cooling systems. A case of study: European greenhouse sector

The work attempts to assess the effects of global warming on the efficacy of current greenhouse cooling methods following a methodology previously proved for other agricultural buildings. The cooling potential of four greenhouse cooling techniques (natural ventilation, forced ventilation, fogging and shading) were simulated by computer modelling for five European locations, calculating the greenhouse internal air temperature from measured external climate data. Four 2080s scenarios were analysed in these five locations. They were constructed as a combination of General Circulation Models (Had CM3 and ECHAM4) downscaled for Europe with the HIRHAM and RCA3 regional models and driven by the A2 and B2 socio-economic scenarios. The crop considered as reference was tomato. The results showed that, in locations in southern Europe, adding evaporative cooling methods to ventilation and/or shading will be indispensable. In some areas of northern Europe, natural ventilation will no longer be sufficient, and shading or fogging will also be necessary. The economic consequences will be important, over all in the southern locations where water consumption, investment and working costs will be higher and necessary to ensure the crop production.     

Spatial Statistical Properties of Pressure Solution Seams in Clastic Rocks in Southwest Ireland

Length, height, thickness and spacing measurements of pressure solution seams at outcrop, hand sample and thin section scale were taken from clastic rocks located in the southwest of Ireland. The lengths and spacings of pressure solution seams have similarly shaped (approximately log-hyperbolic) distributions at the observed scales suggesting that length and spacing distributions are scale-independent over the scales studied with a fractal dimension in the range of 1.4 to 1.6. Pressure solution seam lengths and thicknesses are related by a power-law and their spacings have a linear relationship to bed thickness. Although pressure solution seams are often considered as anticracks (forming under the same remote stresses as joints, but with opposite sign) we describe how the mechanism of pressure solution differs substantially from that of jointing. We use an existing mechanical model to show that stresses around pressure solution seam tips are much lower than those for joints under equal but opposite loading conditions. Pressure solution seams also have a decreasing tendency to lengthen as they grow, which is reflected in their length distributions. We propose that pressure solution seams, unlike joints, do not reach fracture saturation spacing because of transverse coalescence.     

Reservoir operation using system dynamics under climate change impacts: a case study of Yamchi reservoir,Iran

This paper proposes a decision support system for Yamchi reservoir operation in semi-arid region of Iran. The paper consists of the following steps: Firstly, the potential impacts of climate change on the streamflow are predicted. The study then presents the projections of future changes in temperature and precipitation under A2 scenario using the LARS-WG downscaling model and under RCP2.6, RCP4.5, and RCP8.5 using the statistical downscaling model (SDSM) in the northwestern of Iran. To do so, a general circulation model of HadCM3 is downscaled by using the LARS-WG model. As a result, the average temperature, for the horizon 2030 (2011–2030), will increase by 0.77 °C and precipitation will decrease by 11 mm. Secondly, the downscaled variables are used as input to the artificial neural network to investigate the possible impact of climate change on the runoffs. Thirdly, the system dynamics model is employed to model different scenarios for reservoir operation using the Vensim software. System dynamics is an effective approach for understanding the behavior of complex systems. Simulation results demonstrate that the water shortage in different sectors (including agriculture, domestic, industry, and environmental users) will be enormously increased in the case of business-as-usual strategy. In this research, by providing innovative management strategies, including deficit irrigation, the vulnerability of reservoir operation is reduced. The methodology is evaluated by using different modeling tests which then motivates using the methodology for other arid/semi-arid regions.     

Petrographical and geochemical signatures of Jurassic rocks of Chari Formation,Western India: implications for provenance and tectonic setting

The sandstones of the Ridge and Athleta members of Chari Formation (Callovian-Oxfordian) exposed at Jara have been analyzed for their petrographical and geochemical studies. Texturally, these sandstones are medium to coarse grained, poorly to well sorted, sub-angular to sub-rounded, and show low to medium sphericity. These sandstones were derived from a mixed provenance including granites, granite-gneisses, low and high-grade metamorphic, and some basic rocks of Aravalli range and Nagarparkar massif. The petrofacies analysis reveals that these sandstones belong to the continental block and recycled orogen tectonic regime. The studied sandstones are modified by paleoclimate, distance of transport, and diagenesis. Mineralogically and geochemically, sandstones are classified as quartzarenite, subarkose, arkose, sublithic arenite, and wacke, respectively. The A-CN-K ternary plot and CIA, CIW, PIA, and ICV values suggest that the similar source rocks suffered moderate to high chemical weathering under a hot-humid climate in an acidic environment with higher \({\text{P}}_{{{\text{CO}}_{ 2} }}\). Generally good to strong correlations between Al₂O₃ and other oxides in these sediments indicate clay mineral control. The K₂O/Na₂O versus SiO₂ diagram indicates that the studied samples occupy passive margin fields but the SiO₂/Al₂O₃ versus K₂O/Na₂O plot suggests that the Athleta Sandstone and Ridge Sandstone fall within the passive margin field, while Ridge Shale falls within the active continental margin field.     

Mapping subvertical discontinuities in rock cuts using a 400-MHz ground penetrating radar antenna

Hidden subvertical discontinuities oriented parallel to subparallel to the exposed faces of outcropping sandstone were effectively mapped at three different study sites in central Missouri using a ground-penetrating radar system (GPR) equipped with a 400-MHz monostatic antenna and a survey wheel. At each site, a suite of 2-D ground-penetrating radar profiles were acquired along multiple closely spaced traverses on relatively smooth exposed rock surfaces. Time-zero correction was applied to the raw GPR data which were then processed using band-pass filtering, range and display gain, color transformation, and deconvolution techniques. Pseudo 3D images of each identified discontinuity at each site were constructed based on the interpretation of the nonmigrated ground-penetrating radar profiles. These pseudo 3D images were hand-migrated and transformed into true 3D images which depict variable depths at “perpendicular horizontal distance” to each discontinuity relative to the exposed rock face. The results demonstrate that GPR can be used to detect and map hidden discontinuities. This information can then be used for rock slope stability analysis and rock engineering purposes.     

Geochemical and Sr–Nd–Pb isotopic compositions of the Eocene Dölek and Sariçiçek Plutons, Eastern Turkey: Implications for magma interaction in the genesis of high-K calc-alkaline granitoids in a post-collision extensional setting

The major and trace elements and Sr–Nd–Pb isotopes of the host rocks and the mafic microgranular enclaves (MME) gathered from the Dölek and Sariçiçek plutons, Eastern Turkey, were studied to understand the underlying petrogenesis and geodynamic setting. The plutons were emplaced at  43 Ma at shallow depths ( 5 to 9 km) as estimated from Al-in hornblende geobarometry. The host rocks consist of a variety of rock types ranging from diorite to granite (SiO₂ = 56.98–72.67 wt.%; Mg# = 36.8–50.0) populated by MMEs of gabbroic diorite to monzodiorite in composition (SiO₂ = 53.21–60.94 wt.%; Mg# = 44.4–53.5). All the rocks show a high-K calc-alkaline differentiation trend. Chondrite-normalized REE patterns are moderately fractionated and relatively flat [(La/Yb)_N = 5.11 to 8.51]. They display small negative Eu anomalies (Eu/Eu = 0.62 to 0.88), with enrichment of LILE and depletion of HFSE. Initial Nd–Sr isotopic compositions for the host rocks are ε_Nd(43 Ma) = − 0.6 to 0.8 and mostly I_Sr = 0.70482–0.70548. The Nd model ages (T_DM) vary from 0.84 to 0.99 Ga. The Pb isotopic ratios are (²⁰⁶Pb/²⁰⁴Pb) = 18.60–18.65, (²⁰⁷Pb/²⁰⁴Pb) = 15.61–15.66 and (²⁰⁸Pb/²⁰⁴Pb) = 38.69–38.85. Compared with the host rocks, the MMEs are relatively homogeneous in isotopic composition, with I_Sr ranging from 0.70485 to 0.70517, ε_Nd(43 Ma) − 0.1 to 0.8 and with Pb isotopic ratios of (²⁰⁶Pb/²⁰⁴Pb) = 18.58–18.64, (²⁰⁷Pb/²⁰⁴Pb) = 15.60–15.66 and (²⁰⁸Pb/²⁰⁴Pb) = 38.64–38.77. The MMEs have T_DM ranging from 0.86 to 1.36 Ga. The geochemical and isotopic similarities between the MMEs and their host rocks indicate that the enclaves are of mixed origin and are most probably formed by the interaction between the lower crust- and mantle-derived magmas. All the geochemical data, in conjunction with the geodynamic evidence, suggest that a basic magma derived from an enriched subcontinental lithospheric mantle, probably triggered by the upwelling of the asthenophere, and interacted with a crustal melt that originated from the dehydration melting of the mafic lower crust at deep crustal levels. Modeling based on the Sr–Nd isotope data indicates that  77–83% of the subcontinental lithospheric mantle involved in the genesis. Consequently, the interaction process played an important role in the genesis of the hybrid granitoid bodies, which subsequently underwent a fractional crystallization process along with minor amounts of crustal assimilation, en route to the upper crustal levels generating a wide variety of rock types ranging from diorite to granite in an extensional regime.