Using inverse modeling and hierarchical cluster analysis for hydrochemical characterization of springs and Talkhab River in Tang-Bijar oilfield,Iran

Hierarchical cluster analysis (HCA) and inverse modeling (PH REdox EQuilibrium (in C language) (PHREEQC)) were simultaneously useful approaches in interpreting surface water hydrochemistry within Talkhab River in the Tang-Bijar oilfield, Iran, where large uncertainties exist in the understanding of the water quality system. Q-mode HCA applied to the data revealed three major surface water associations distinguished on the basis of the major causes of variation in the hydrochemistry. The three water groups were classified as upstream waters (group 1: Ca–SO₄), intermediate waters (group 2: Ca–SO₄–Cl), and downstream waters (group 3: Na–Cl). Geochemical reaction models were constructed using PHREEQC to establish the reactions associated with the different mineral phases through inverse modeling. The hydrochemical compositions of the water groups and the mass balance calculations indicate that the dominant processes and reactions responsible for the hydrochemical evolution in the system are (1) dissolution of evaporites, (2) precipitation of carbonate minerals, (3) silicate weathering reactions, (4) limited mixing with saline water, and (5) ion exchange.     

Hydrogeological assessment of the Tabarteh fault zone by physicochemical analysis and multivariate statistical methods

Hydrogeologically, faults may impede, conduit, exert no influence, or may play a combination of these roles on groundwater flow. The object of this paper is to study the hydrogeological role of the Tabarteh fault, which is located on the border of Zagros and Central Iran tectonic zones in an alluvial aquifer. The recorded data of water table levels, chemical parameters, and discharge rate of wells, in addition to geological maps and geophysical results, were collected and evaluated. The outcrop of travertine in limited areas and the emergence of a few small springs within the alluvium show a barrier role of the fault in the groundwater flow. The spatial analysis of chemical components, head time series, and groundwater flow direction assessment demonstrated that the fault acts as both a barrier and a non-barrier in different sections. The multivariate statistical methods of cluster and discriminant analyses also confirm the dual role of the fault.     

Three-dimensional supervirtual seismic refraction interferometry

Improved supervirtual interferometry (SVI) method that enhance the signal-to-noise ratio (SNR) of noisy seismic refraction data is presented. This 3D supervirtual algorithm is helpful in removal of side lobes caused by convolution in the conventional 3D SVI method. In this study, we extend the supervirtual seismic refraction interferometry method to 3D geometries commonly used in active seismic exploration. To achieve this objective, synthetic 3D seismic refraction data were created using single patch orthogonal geometry. Simple two-layers velocity model was used. A zero phase wavelet was convolved with the refraction travel time to create the seismic wiggle trace. The refraction arrival is computed using the equation for horizontal refractor case and is not adapted for dipping refractor. Noise of Gaussian distribution with zero mean and 0.25 standard deviation was added to simulate a case of moderate ambient noise. The supervirtual algorithm, consisting of crosscorrelation, alignment, summation, and first arrival calculation is performed. By aligning and summation of all the correlogram, the stationary position of source-receiver pairs as required in the 2D supervirtual method are eliminated in this case. Synthetic data presented in this study shows accurate first arrivals after the application of the 3D SVI and traces with much better SNR than the actual traces.     

Anomalous optical variations in the grossular garnet from the Eden Mills,Belvidere Mountain (Vermont,USA)

The optical anomalies, and surface and lamellar textures of a birefringent grossular garnet crystal from the Eden Mills, Belvidere Mountain, Vermont, USA, have been investigated by optical polarizing microscope, electron-probe micro-analyzer, back-scattered electron imaging, infrared spectroscopy, and single-crystal X-ray diffractometer from the standpoint of crystal growth. This grossular shows one-to-one correlation between natural surface features and its internal textures under crossed polarizers. Electron-probe micro-analyzer (EPMA) gave average chemical composition in (110) thin section, of bright lamella {Ca_2.97Mn_0.06}_∑3.03 [Al_1.59Fe_0.37Ti_0.01]_∑1.97(Si_3.00)_∑3.00 (Gros_79.5And_18.9Sps_1.6) and of dark host {Ca_2.99Mn _0.06}_∑3.05 [Al_1.73Fe_0.26 Ti_0.01]_∑2.00(Si_2.97OH_0.03)_∑3.00 (Gros_85.4And₁₃Sps_1.6). The correspondence of surface features and the internal textures with spiral or pyramidal growth mechanism suggest that the internal textures of the Eden Mills grossular are formed during growth process. The optical vibrational orientations and the growth steps inclination along [001] and \( \left[\overset{-}{1}10\right] \) directions predict monoclinic symmetry. With X-ray diffractometer (XRD) method, pseudocubic parameters are a = 11.839(2) Å, b = 11.855(1) Å, and c = 11.868(2) Å with interaxial angles α = 90.00(1)°, β = 89.99(1)°, and γ = 90.02(2)° that show orthorhombic symmetry of this crystal. Lamellar texture of Al³⁺-rich host with Fe³⁺-rich lamella infers cation ordering at octahedral site of the garnet structure. IR data favors the non-cubic orientation of [(OH) ₄] at tetrahedral position in this grossular structure.     

Characteristics of groundwater flow field after land creation engineering in the hilly and gully area of the Loess Plateau

In order to study the degree of influence and control mechanism to groundwater flow field caused by land creation engineering in the hilly and gully area of the Loess Plateau, based on the geological and engineering conditions of the first stage project of Yan’an new district in China, numerical simulation of groundwater flow is carried out by the Feflow and GIS technologies. From the simulation, punning measure relatively reduces infiltration recharge and artificial gravel drain increases groundwater seepage. The basic characteristics of groundwater flow field is controlled by the old and new topographies in the whole study area, and artificial gravel drain plays an auxiliary role in accelerating groundwater drawdown upstream and promotes groundwater rise downstream. According to differences of groundwater level and declining percentages of hydraulic gradient in the main and secondary gullies, dewatering of artificial gravel drain in the secondary gully is more effective than that in the main gully, which will yet play an important role in the future. The study results will make contributions to understand groundwater response to land creation engineering and will be beneficial to take necessary measures to prevent collapse of loess and failure of building foundation in the hilly and gully area of the Loess Plateau.     

Effect of different soil data on hydrological process modeling in Weihe River basin of Northwest China

Hydrological process modeling depends on the soil data spatial resolution of the watershed. Especially, in a large-scale watershed, could a higher resolution of soil data contribute to a more accurate result? In this study, two soil datasets with different classification systems FAO (World Reference Base) and GSCC (the Genetic Soil Classification of China) were used as inputs for the SWAT model to study the effects of soil datasets on hydrological process modeling in Weihe River basin, China. Results show that the discharge simulated using FAO soil data was better than one simulated using GSCC soil data before model calibration, which indicates that FAO soil data needed less effort to calibrate. After model calibration, discharges were simulated better by both of FAO and GSCC soil data but statistical parameters demonstrate that we can make a relatively more accurate estimation of discharge using the GSCC rather than FAO soil data. Soil water content (SW) simulated using GSCC soil data was statistically significantly higher than those simulated using FAO soil data. However, variations in other hydrological components (surface runoff (SURQ), actual evapotranspiration (ET), and water yield (WYLD) were not statistically significant. This might be because SW is more sensitive to soil properties. For studies aiming to simulate or compare SW, merely calibrating and validating models using river discharge observations is not enough. The hydrological modelers need to identify the key hydrological components intrinsic to their study and weigh the advantages and disadvantages before selecting suitable soil data.     

Stochastic groundwater modeling of a sedimentary aquifer: evaluation of the impacts of abstraction scenarios under conditions of reduced recharge

A transient finite difference groundwater flow model has been calibrated for the Nasia sub-catchment of the White Volta Basin. This model has been validated through a stochastic parameter randomization process and used to evaluate the impacts of groundwater abstraction scenarios on resource sustainability in the basin. A total of 1500 equally likely model realizations of the same terrain based on 1500 equally likely combinations of the data of the key aquifer input parameters were calibrated and used for the scenario analysis. This was done to evaluate model non-uniqueness arising from uncertainties in the key aquifer parameters especially hydraulic conductivity and recharge by comparing the realizations and statistically determining the degree to which they differ from each other. Parameter standard deviations, computed from the calibrated data of the key parameters of hydraulic conductivity and recharge, were used as a yardstick for evaluating model non-uniqueness. All model realizations suggest horizontal hydraulic conductivity estimates in the range of 0.03–78.4 m/day, although over 70 % of the area has values in the range of 0.03–14 m/day. Low standard deviations of the horizontal hydraulic conductivity estimates from the 1500 solutions suggest that this range adequately reflects the properties of the material in the terrain. Lateral groundwater inflows and outflows appear to constitute significant components of the groundwater budgets in the terrain, although estimated direct vertical recharge from precipitation amounts to about 7 % of annual precipitation. High potential for groundwater development has been suggested in the simulations, corroborating earlier estimates of groundwater recharge. Simulation of groundwater abstraction scenarios suggests that the domain can sustain abstraction rates of up to 200 % of the current estimated abstraction rates of 12,960 m³/day under the current recharge rates. Decreasing groundwater recharge by 10 % over a 20-year period will not significantly alter the results of this abstraction scenario. However, increasing abstraction rates by 300 % over the period with decreasing recharge by 10 % will lead to drastic drawdowns in the hydraulic head over the entire terrain by up to 6 m and could cause reversals of flow in most parts of the terrain.     

Textural characteristics of sediments along the southern Red Sea coastal areas,Saudi Arabia

A lack of understanding exists of the origin and textural characteristics of Saudi Arabian Red Sea coastal sediments. This paper concerns the southern coastline of Jizan on the Saudi Red Sea. It is some 160 km long characterised by either narrow rocky headlands with intermittent pocket beaches or wide low-lying beaches dissected by wadis. Granulometric testing of samples from 135 locations showed that beach sand size was mainly very fine to medium grained (M _z = 3.93 Ø), sorting ranged from 1.65 to 0.41 and skewness values from ?051 to 0.39, being mainly negative; dune sands were medium to fine grained (M _z = 1.13 Ø; average sorting 2.8), while skewness variations within dune samples indicated symmetrical to fine skewed values (б _Ι = 0.55 to 0.89). Most foreshore samples were derived from wadis. Wadi mud levels can be high, e.g. Baysh (84%), and wadi Samrah (90%) with mean grain size ranging from very fine to medium sand (M _z = 3.9 Ø), sorting being well to poor (0.45 to 1.52) due to sediment influxes. Sabkha had a wide range of sand/mud and significantly higher carbonate percentages than other environments. Sediment source differences and littoral reworking contributed to grain size variation. The carbonate content varied between 1.5 and 31.5% due to hinterland contributions, and spatial analysis showed increasing quantities of carbonate minerals towards the south. On the wider geographical front, findings from Jizan are similar to those of the Northern United Arab Emirates (UAE), including sabkhas, being composed of sand, skeletal carbonate, fine fluvial material and wind-blown silt and clay components of wadi origin. Further work on the northeastern Red Sea edge can hopefully confirm these findings.     

Using a Thermal Proxy to Examine Sediment–Water Exchange in Mid-Continental Shelf Sandy Sediments

Fluid exchange across the sediment–water interface in a sandy open continental shelf setting was studied using heat as a tracer. Summertime tidal oscillation of cross-shelf thermal fronts on the South Atlantic Bight provided a sufficient signal at the sediment–water interface to trace the advective and conductive transport of heat into and out of the seabed, indicating rapid flushing of ocean water through the upper 10–40 cm of the sandy seafloor. A newly developed transport model was applied to the in situ temperature data set to estimate the extent to which heat was transported by advection rather than conduction. Heat transported by shallow 3-D porewater flow processes was accounted for in the model by using a dispersion term, the depth and intensity of which reflected the depth and intensity of shallow flushing. Similar to the results of past studies in shallower and more energetic nearshore settings, transport of heat was greater when higher near-bed velocities and shear stresses occurred over a rippled bed. However, boundary layer processes by themselves were insufficient to promote non-conductive heat transport. Advective heat transport only occurred when both larger boundary layer stresses and thermal instabilities within the porespace were present. The latter process is dependent on shelf-scale heating and cooling of bottom water associated with upwelling events that are not coupled to local-scale boundary layer processes.     

Si and O self-diffusion in hydrous forsterite and iron-bearing olivine from the perspective of defect chemistry

We discuss the experimental results of silicon and oxygen self-diffusion coefficients in forsterite and iron-bearing olivine from the perspective of defect chemistry. Silicon diffusion is dominated by V_O ^··-associated V_Si″″, whereas oxygen diffusion is dominated by hopping of V_O ^·· under anhydrous conditions, and by (OH)_O ^· under hydrous conditions. By considering the charge neutrality condition of [(OH)_O ^·] = 2[V_Me″] in hydrous forsterite and iron-bearing olivine, we get D _Si ∝ (\(C_{{{\text{H}}_{2} {\text{O}}}}\))^1/3 and D _O ∝ (\(C_{{{\text{H}}_{2} {\text{O}}}}\))⁰, which explains the experimental results of water effects on oxygen and silicon self-diffusion rates (Fei et al. in Nature 498:213–215, 2013; J Geophys Res 119:7598–7606, 2014). The \(C_{{{\text{H}}_{2} {\text{O}}}}\) dependence of creep rate in the Earth’s mantle should be close to that given by Si and O self-diffusion coefficients obtained under water unsaturated conditions.