Characterizing groundwater recharge using oxygen and hydrogen isotopes: a case study in a temperate forested region,South Korea

In the cool temperate region of South Korea, oxygen and hydrogen isotopes of groundwater, lake water, and precipitation were studied to determine the season of groundwater recharge. All the groundwater samples, irrespective of season, on δ¹⁸O–δ²H scale plotted along the summer precipitation, suggesting summer precipitation largely modulates recharge. The deuterium excess values of groundwater (d-excess) show clear seasonal difference, higher in winter (> 18‰) and lower in summer (< 10‰). And its resemblance to the summer precipitation d-excess value further suggests dominant role of summer precipitation in groundwater recharge. Based on the mass balance equation, with end-member d-excess values of seasonal precipitation and groundwater as input variables, groundwater is composed of 66% summer and 34% winter precipitation. Despite the study area being heavily forested, summer rainfall contribution higher than winter suggests that evapotranspiration effect is minimal in the region; may be due to thin sand–gravel-based porous soil overlying highly weathered granitic rock system.     

Evaluating the impact of artificial groundwater recharge structures using geo-spatial techniques in the hard-rock terrain of Rajasthan,India

Groundwater levels in hard-rock areas in India have shown very large declines in the recent past. The situation is becoming more critical due to a paucity of rainfall, limited surface water resources and an increasing pattern of groundwater extraction in these areas. Consequently, the Ground Water Department with the aid of World Bank has implemented the water structuring programme to mitigate groundwater scarcity and to develop a viable solution for sustainable development in the region. The present study has been undertaken to assess the impact of artificial groundwater recharge structures in the hard-rock area of Rajasthan, India. In this study groundwater level data (pre-monsoon and post-monsoon) of 85 dug-wells are used, spread over an area of 413.59 km². The weathered and fractured gneissic basement rocks act as major aquifer in the area. Spatial maps for pre- and post-monsoon groundwater levels were prepared using the kriging interpolation technique with best fitted semi-variogram models (Spherical, Exponential and Gaussian). The groundwater recharge is calculated spatially using the water level fluctuation method. The entire study period (2004–2011) is divided into pre- (2004–2008) and post-intervention (2009–2011) periods. Based on the identical nature of total monsoon rainfall, two combinations of average (2007 and 2009) and more than average (2006 and 2010) rainfall years are selected from the pre- and post-intervention periods for further comparisons. All of the water harvesting structures are grouped into the following categories: as anicuts (masonry overflow structure); percolation tanks; subsurface barriers; and renovation of earthen ponds/nadis. A buffer of 100 m around the intervention site is taken for assessing the influence of these structures on groundwater recharge. The relationship between the monsoon rainfall and groundwater recharge is fitted by power and exponential functions for the periods of 2004–2008 and 2008–2011 with R ² values of 0.95 and 0.98, respectively. The average groundwater recharge is found to be 18% of total monsoon rainfall prior to intervention and it became 28% during the post-intervention period. About 70.9% (293.43 km²) of the area during average rainfall and more than 95% (396.26 km²) of the area during above-average rainfalls show an increase in groundwater recharge after construction of water harvesting structures. The groundwater recharge pattern indicates a positive impact within the vicinity of intervention sites during both average and above-average rainfall. The anicuts are found to be the most effective recharge structures during periods of above-average rainfall, while subsurface barriers are responded well during average rainfall periods. In the hard-rock terrain, water harvesting structures produce significant increases in groundwater recharge. The geo-spatial techniques that are used are effective for evaluating the response of different artificial groundwater recharge techniques.     

The stress regime in a Rotliegend reservoir of the Northeast German Basin

In-situ stresses have significant impact, either positive or negative, on the short and long term behaviour of fractured reservoirs. The knowledge of the stress conditions are therefore important for planning and utilization of man-made geothermal reservoirs. The geothermal field Groß Schönebeck (40 km north of Berlin/Germany) belongs to the key sites in the northeastern German Basin. We present a stress state determination for this Lower Permian (Rotliegend) reservoir by an integrated approach of 3D structural modelling, 3D fault mapping, stress ratio definition based on frictional constraints, and slip-tendency analysis. The results indicate stress ratios of the minimum horizontal stress S _hmin being equal or increasing 0.55 times the amount of the vertical stress S _V (S _hmin ≥ 0.55S _V ) and of the maximum horizontal stress S _Hmax ≤ 0.78–1.00S _V in stress regimes from normal to strike slip faulting. Thus, acting stresses in the 4,100-m deep reservoir are S _V  = 100 MPa, S _hmin = 55 MPa and S _Hmax = 78?100 MPa. Values from hydraulic fracturing support these results. Various fault sets of the reservoir are characterized in terms of their potential to conduct geothermal fluids based on their slip and dilatation tendency. This combined approach can be adopted to any other geothermal site investigation.     

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.     

Changes in flow and chemistry of groundwater heavily affected by human impacts in the Baiyangdian catchment of the North China Plain

To identify impacts of air pollution, sewage drainage, agricultural production, over-pumping and reservoir storage on groundwater, a field survey was conducted in the Baiyangdian catchment of the North China Plain. Major ions and water isotopes were measured. Results show that hydrological processes and hydrogeochemical evolution of shallow groundwater were greatly disturbed by human activities. Excessive pumping resulted in significant declines of groundwater levels over the study area. This also induced infiltration of surface water into groundwater. A groundwater depression cone was the conflux center of groundwater surrounded by recharge zones including alluvial fans and surface water in alluvial plain. Pumping almost was the only way to discharge groundwater. Emission of SO _x and NO _x contributed at least 11% of rock weathering by dissolving into infiltrating precipitation. Surface waters containing sewage replenished ambient groundwater with an average mixing ratio of 74 ± 17% due to groundwater level drawdown. As a result, groundwater had elevated concentrations of Na⁺ and SO₄ ^2? with Na⁺ exchanged into aquifer sediments. About 29 ± 16% of Na⁺ was exchanged from groundwater into soil matrix. Agriculture nitrate was high only in the recharge zones. The most important result is that the transformation of the study area from a place rich in water resource into an area lack of water just took several decades with the joint action of the heavily human activities. Our study also indicates that shallow groundwater could sensitively respond to and record environmental changes.     

Groundwater flow dynamics of weathered hard-rock aquifers under climate-change conditions: an illustrative example of numerical modeling through the equivalent porous media approach in the north-western Pyrenees (France)

A numerical groundwater model of the weathered crystalline aquifer of Ursuya (a major water source for the north-western Pyrenees region, south-western France) has been computed based on monitoring of hydrological, hydrodynamic and meteorological parameters over 3 years. The equivalent porous media model was used to simulate groundwater flow in the different layers of the weathered profile: from surface to depth, the weathered layer (5?·?10^?8?≤?K?≤?5?·? 10^?7 m s^?1), the transition layer (7?·?10^?8?≤?K?≤?1?·? 10^?5 m s^?1, the highest values being along major discontinuities), two fissured layers (3.5?·?10^?8?≤?K?≤?5?·?? 10^?4 m s^?1, depending on weathering profile conditions and on the existence of active fractures), and the hard-rock basement simulated with a negligible hydraulic conductivity (K = 1 10 ^?9 ). Hydrodynamic properties of these five calculation layers demonstrate both the impact of the weathering degree and of the discontinuities on the groundwater flow. The great agreement between simulated and observed hydraulic conditions allowed for validation of the methodology and its proposed use for application on analogous aquifers. With the aim of long-term management of this strategic aquifer, the model was then used to evaluate the impact of climate change on the groundwater resource. The simulations performed according to the most pessimistic climatic scenario until 2050 show a low sensitivity of the aquifer. The decreasing trend of the natural discharge is estimated at about ?360 m³ y^?1 for recharge decreasing at about ?5.6 mm y^?1 (0.8 % of annual recharge).     

Testing Study of Subcritical Crack Growth Mechanism During Water Rock Interaction

Subcritical crack growth plays an important role in evaluating the long-term stability of structures in rocks. By applying the constant-displacement-relax method, two groups of test specimens that one immersed in groundwater and the other in air were tested to get the parameters of subcritical crack growth in double torsion test. The relations of the stress intensity factor K _I versus the subcritical crack growth velocity V were obtained under the two different environments, and the behavior of subcritical crack growth was also analyzed. The results showed: the relations of lg K _I  ? lg V accorded with linear rules, which is in good agreement with Charles theory; Compared with specimens in nature state, the lg K _I  ? lg V curves of saturated water specimens locate top left comer of those of air specimens. The slope of curve is smaller, and the intercept is bigger, which shows that the water–rock interaction speeds up the subcritical crack growth. And A increases 2.9 × 10¹⁸ folds but n decreases from 85.12 to 40.83 because of the water–rock interaction. Meanwhile, the fracture toughness K _IC also decreases obviously from 2.55 in air to 2.26 in water due to water rock interaction. The testing results provide a basis for time-dependence of rock engineering stability.     

Characterizing land displacement in complex hydrogeological and geological settings: a case study in the Beijing Plain,China

Characterization of land displacement induced by long-term overexploitation of groundwater is necessary to ensure sustainable water supply in Beijing, China. The northern part of the Beijing Plain is an important water source area and is also designed for groundwater recharge from South-to-North Water Diversion Project. We aim to depict the process of characterizing land displacement under complex hydrogeological and geological context in the region using remote sensing and geographic information system. Interferometric synthetic aperture radar time-series analysis was used to detect land displacement from 2003 to 2010. Statistic linear regression equations between groundwater level and land displacement were built based on linear consolidation principle. The spatial difference of Pearson correlation coefficient (R) and slope (k) were discriminated to quantify the response of land displacement to groundwater level change. The results show that there are two major displacement cones with annual rates up to ?40 and ?24 mm year^?1. R and k had a negative and positive correlation with increasing land displacement, respectively. A larger R reflects that the groundwater level has a closer relation with the occurrence of land displacement. The weak correlation is due to the delay in the propagation of the pressure drawdown in the fine-sediment layers or lens from the pumped aquifers where the pressure is measured. Thick compressible layer has more potential for land displacement. Results of this study are necessary to clarify the land displacement characteristics, to make full use of abundant spatial–temporal dataset, and ultimately to support hazard prevention and mitigation decisions.     

Sequence signal processing of 2012 Ahar–Varzaghan earthquake (<Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript> 6.4) of NW Iran

Before starting seismic cycle of Ahar–Varzaghan 2012 event, a partial gap in the form of a pre-seismic calm sequence (seismicity rate, r = 0.46 event/year, b = 1.4) with duration of 303 days spatially has dominated over the entire seismogenic area. From April 17, 2012, to May 31, 2012, r significantly increased to 2.16, indicating strong foreshock sequence, and b value changed to 1.9, remarkably. In the last two months before the mainshock, foreshocks have partially migrated toward the earthquake fault (with a decrease in size, b = 2.0). Significantly, high rate of seismicity and low V _P /V _S (1.64) in the foreshocks sequence and also very high seismicity rate (17.3) and high V _P /V _S (1.76) in the aftershocks sequence make substantial differences between the seismic cycle and the background seismicity. Moreover, a significant E–W migration of the microseismicity was confirmed in the study area.     

The transfer of oxygen isotopic signals from precipitation to drip water and modern calcite on the seasonal time scale in Yongxing Cave,central China

Stable isotope data of precipitation (δ¹⁸O_p and deuterium excess), drip water (δ¹⁸O_d), and modern calcite precipitates (δ¹⁸O_c and δ¹³C_c) from Yongxing Cave, central China, are presented, with monthly sampling intervals from June 2013 to September 2016. Moderate correlations between the monthly variation of δ¹⁸O_p values (from ??11.5 to ??0.7‰) and precipitation amount (r = ??0.59, n?=?34, p?<?0.01) and deuterium excess (r?=?0.39, n?=?31, p?<?0.01) imply a combined effect of changes in precipitation amount and atmospheric circulation. At five drip sites, the δ¹⁸O_d values have a much smaller variability (from ??9.1 to ??7.5‰), without seasonal signals, probably a consequence of the mixing in the karst reservoir with a deep aquifer. The mean δ¹⁸O_d value (??8.4‰) for all drip waters is significantly more negative than the mean δ¹⁸O_p value (??6.9‰) weighted by precipitation amount, but close to the wet season (May to September) mean value (??8.3‰), suggesting that a threshold of precipitation amount must be exceeded to provide recharge. Calculation based on the equilibrium fractionation factor indicates that the δ¹⁸O_c values are not in isotopic equilibrium with their corresponding drip waters, with a range of disequilibrium effects from 0.4 to 1.4‰. The δ¹⁸O_c and δ¹³C_c values generally increase progressively away from the locus of precipitation on glass plates. The disequilibrium effects in the cave are likely caused by progressive calcite precipitation and CO₂ degassing related to a high gradient of CO₂ concentration between drip waters and cave air. Our study provides an important reference to interpret δ¹⁸O_c records from the monsoon region of China.     

Application of electrical resistivity method in estimating geohydraulic properties of a sandy hydrolithofacies: a case study of Ajali Sandstone in Ninth Mile,Enugu State,Nigeria

Geoelectric investigation using vertical electrical sounding (VES) (Schlumberger electrode configuration) was carried out in 14 locations at Ninth Mile area, southeastern Nigeria to determine the variations and interrelationship of some geoelectric and geohydraulic parameters of a sandstone hydrolithofacies. The measured resistivity data were interpreted using manual and computer software packages, which gave the resistivity, depth, and thickness for each layer within the maximum current electrodes separation. The aquifer resistivity values range from 86.56 to 4753.0 Ωm with 1669.40 Ωm average value. The values of water resistivity from borehole locations close to the sounding points range from 79.49 to 454 .55 Ωm and averaging about 264.7 Ωm. Porosity values of the sandy aquifer range from 30.19 to 34.20%. Fractional porosity values range from 0.3019 to 0.3292, while the tortuosity values vary between 2.91 and 22.85. The geohydraulic parameters estimated vary across the study area. Formation factor ranges from 0.28 to 15.29, hydraulic conductivity ranges from 1.21 to 66.54 m/day which, however, influences the natural flow of water in the aquifer while tortuosity values range from 2.91 to 23.27. The contour maps clearly show the variation of these parameters in the subsurface and the plots show their relationship and high correlation coefficients with one another. The results of this study have revealed the geological characteristics of the subsurface aquifer, established the influence on the amount of groundwater, and proposed a strategy for the management and exploitation of groundwater resources in the area and other aquiferous formations.     

Application of WetSpass model to estimate groundwater recharge variability in the Nile Delta aquifer

Long-term groundwater recharge from rainfall in the Nile Delta is needed as an input for integrated groundwater modelling in the Nile Delta aquifer for more accurate simulation. The main objective is to estimate the spatial and temporal variation of groundwater recharge from rainfall in the Nile Delta aquifer. Water and Energy Transfer between Soil, Plants and Atmosphere under quasi-Steady State (WetSpass) model parameters were identified for the Nile Delta based on the available meteorological data for the area collected in 1991 and 2000. The collected data were rainfall, temperature, wind speed and evapotranspiration. Geomorphological characteristics, such as soil type, topography, groundwater depth and slope, were also collected as input data for the WetSpass model. ENVI software was used to come up with land use classification based on available land cover images of the Nile Delta for 1972, 1984, 1990, 2000 and 2009. The WetSpass model was calibrated by comparing the simulated groundwater recharge with the calculated one by using the water balance equation model. The results indicated close agreement in groundwater recharge between the two model outputs with R ² of 0.99 and 0.94, while the root-mean-square errors (RMSEs) were 4.86 and 9.39 mm for 1991 and 2000, respectively. The WetSpass model was then applied in respect of 1970, 1980, 1990 and 2010 for the purpose of validation. The overall RMSE and R ² for the 6 years were 8.83 mm and 0.88, respectively. The results of the WetSpass calibrated model provide information to support integrated groundwater modelling. The results reveal that WetSpass works well in simulating the components of the hydrological balance in the Nile Delta.     

Coral reefs as a tool for monitoring environmental and radiation hazards

Analysis of coral reef species enables characterization of the coral reefs and the degree to which human activities influence the reef composition. Geochemical characteristics, X-ray diffraction, and natural radionuclides analyses of four common coral reef species representing the branching and massive forms were conducted at the three areas under threat (Quseir Harbor, Safaga Harbor, and El-Esh area) along the Red Sea coast. Branching recorded higher carbonates and lower organic matter than the massive forms. Stylophora pistillata at Safaga Harbor and El-Esh area and Acropora humilis at Qusier Harbor recorded the highest carbonate percentages (96.88, 96.24, and 96.89%, respectively) meanwhile Platygyra daedalea at Safaga Harbor recorded the highest organic matter contents (5.07 and 4.91%). The highest Ca contents were observed in S. pistillata and Porites solida (65.87 and 64.87%), the highest Mg was in Acropora humilis and P. daedalea (1.06 and 0.88%) at Qusier Harbor, whilst the highest Sr was in S. pistillata and P. daedalea at Safaga Harbor. Then highest value of ²²⁶Ra recorded in A. humilis and P. solida (71 ± 3.5 Bq/kg and 63 ± 3.1 Bq/k) and ²³²Th in S. pistillata at El-Esh area may attribute to petroleum exploration activities. A. humilis and P. daedalea at Safaga Harbor recorded significant high ⁴⁰K values (505 ± 30 and 472 ± 13 Bq/kg, respectively) relative to the other localities. The low Ca and High Sr as well as the highest averages of ²³²Th and ⁴⁰K at Safaga Harbor indicated that the coral reef communities were highly affected by the terrestrial inputs from phosphate shipments.     

Assessing the recharge process and importance of montane water to adjacent tectonic valley-plain groundwater using a ternary end-member mixing analysis based on isotopic and chemical tracers

A study in eastern Taiwan evaluated the importance of montane water contribution (MC) to adjacent valley-plain groundwater (VPG) in a tectonic suture zone. The evaluation used a ternary natural-tracer-based end-member mixing analysis (EMMA). With this purpose, VPG and three end-member water samples of plain precipitation (PP), mountain-front recharge (MFR), and mountain-block recharge (MBR) were collected and analyzed for stable isotopic compositions (δ²H and δ¹⁸O) and chemical concentrations (electrical conductivity (EC) and Cl^?). After evaluation, Cl^? is deemed unsuitable for EMMA in this study, and the contribution fractions of respective end members derived by the δ¹⁸O–EC pair are similar to those derived by the δ²H–EC pair. EMMA results indicate that the MC, including MFR and MBR, contributes at least 70% (679?×?10⁶ m³ water volume) of the VPG, significantly greater than the approximately 30% of PP contribution, and greater than the 20–50% in equivalent humid regions worldwide. The large MC is attributable to highly fractured strata and the steep topography of studied catchments caused by active tectonism. Furthermore, the contribution fractions derived by EMMA reflect the unique hydrogeological conditions in the respective study sub-regions. A region with a large MBR fraction is indicative of active lateral groundwater flow as a result of highly fractured strata in montane catchments. On the other hand, a region characterized by a large MFR fraction may possess high-permeability stream beds or high stream gradients. Those hydrogeological implications are helpful for water resource management and protection authorities of the studied regions.     

High-pressure displacive phase transition of a natural Mg-rich pigeonite

A high-pressure single-crystal X-ray diffraction study has been carried out on a P2₁/c natural Mg-rich pigeonite sample with composition ca. Wo₆En₇₆Fs₁₈ using a diamond anvil-cell. The unit-cell parameters were determined at 14 different pressures to 7.14 GPa. The sudden disappearance of the b-type reflections (h + k = odd) and a strong discontinuity (about 2.8%) in the unit-cell volume indicated a first-order P2₁/c–C2/c phase transition between 4.66 and 4.88 GPa. The P(V) data of the P2₁/c phase were fitted to 4.66 GPa by a third-order Birch–Murnaghan equation of state (BM3 EoS), whereas the limited number of experimental data collected within the C2/c phase between 4.88 and 7.14 GPa were fitted using the same equation of state but with K′ constrained to the value obtained for the P2₁/c fitting. The equation of state coefficients are V ₀ = 424.66(6) ų, K _T0 = 104(2) GPa and K′ = 8(1) for the P2₁/c phase, and V ₀ = 423.6(1) ų, K _T0 = 112.4(8) GPa, and K′ fixed to 8(1) for the C2/c phase. The axial moduli for a, b, and c for the P2₁/c phase were obtained using also a BM3-EoS, while for the C2/c phase only a linear calculation could be performed, and therefore the same approach was applied for comparison also to the P2₁/c phase. In general the C2/c phase exhibits axial compressibilities (β _c > β _a >> β _b) lower than those of the P2₁/c phase (β _b > β _c ≈ β _a; similar to those found in previous studies in clinopyroxenes and orthopyroxenes). The lower compressibility of the C2/c phase compared with that of the P2₁/c could be ascribed to the greater stiffness along the b direction. A previously published relationship between P _c and M2 average cation radius (i.r.) has been updated using all the literature data on P2₁/c clinopyroxene containing large cations at M2 site and our new data. The following weighted regression was obtained: P _c (GPa) = 26(4) ? 28(5) ×  i.r (Å), R ² = 0.97. This improved equation can be used to predict the critical pressure of natural P2₁/c clinopyroxene samples just knowing the composition at M2 site.     

Permeability of granular soil employing flexible wall permeameter

Understanding the changes in permeability of soil, when soil is subjected to high confining pressure and flow pressure, which may alter the textural and geomechanical characteristics of soil, is of great importance to many geo-engineering activities such as, construction of high-rise buildings near the coast or the water bodies, earthen dams, pavement subgrades, reservoir, and shallow repositories. It is now possible to evaluate the changes in permeability of soil samples under varying conditions of confining pressure and flow pressure using flexible wall permeameter (FWP). In the present study, investigation was carried out on a cylindrical sample of granular soil employing FWP under varied conditions of confining pressure (σ₃)—50–300 kPa, which can simulate the stress conditions equivalent to depth of about 20 m under the earth’s crust, and a flow pressure (f_p)—20–120 kPa, which is mainly present near the small earthen embankment dams, landfill liners, and slurry walls near the soft granular soil with high groundwater table. The obtained results indicate a linear relationship between hydraulic conductivity (k) with effective confining pressure (σ_eff.), k, decreasing linearly with an incremental change in σ_eff.. Further, k increases significantly with an increase in f_p corresponding to each σ_eff., and q increases significantly with increase in the f_p corresponding to each (σ₃). It was also observed that corresponding to the low f_p of 20 kPa, the reduction in k is nonlinear with σ_3. The percentage reduction in k is observed to be 9, 13, and 27% corresponding to σ₃ of 50–100, 100–200, and 200-300 kPa, respectively.     

Accuracy assessment of near-shore bathymetry information retrieved from Landsat-8 imagery

The improvement in the capabilities of Landsat-8 imagery to retrieve bathymetric information in shallow coastal waters was examined. Landsat-8 images have an additional band named coastal/aerosol, Band 1: 435–451 nm in comparison with former generation of Landsat imagery. The selected Landsat-8 operational land image (OLI) was of Chabahar Bay, located in the southern part of Iran (acquired on February 22, 2014 in calm weather and relatively low turbidity). Accurate and high resolution bathymetric data from the study area, produced by field surveys using a single beam echo-sounder, were selected for calibrating the models and validating the results. Three methods, including traditional linear and ratio transform techniques, as well as a novel proposed integrated method, were used to determine depth values. All possible combinations of the three bands [coastal/aerosol (CB), blue (B), and green (G)] have been considered (11 options) using the traditional linear and ratio transform techniques, together with five model options for the integrated method. The accuracy of each model was assessed by comparing the determined bathymetric information with field measured values. The standard error of the estimates, correlation coefficients (R ² ) for both calibration and validation points, and root mean square errors (RMSE) were calculated for all cases. When compared with the ratio transform method, the method employing linear transformation with a combination of CB, B, and G bands yielded more accurate results (standard error = 1.712 m, R ² _calibration = 0.594, R ² _validation = 0.551, and RMSE =1.80 m). Adding the CB band to the ratio transform methodology also dramatically increased the accuracy of the estimated depths, whereas this increment was not statistically significant when using the linear transform methodology. The integrated transform method in form of Depth = b ₀  + b ₁ X _CB  + b ₂ X _B  + b ₅ ln(R _CB )/ln(R _G ) + b ₆ ln(R _B )/ln(R _G ) yielded the highest accuracy (standard error = 1.634 m, R ² _calibration = 0.634, R ² _validation = 0.595, and RMSE = 1.71 m), where R _i (i = CB, B, or G) refers to atmospherically corrected reflectance values in the i ^th band [X _i  = ln(R _i -R _{deep water})].     

Elastic behavior of vanadinite,Pb<Subscript>10</Subscript>(VO<Subscript>4</Subscript>)<Subscript>6</Subscript>Cl<Subscript>2</Subscript>, a microporous non-zeolitic mineral

The high-pressure behavior of a vanadinite (Pb₁₀(VO₄)₆Cl₂, a = b = 10.3254(5), c = 7.3450(4) Å, space group P6₃/m), a natural microporous mineral, has been investigated using in-situ HP-synchrotron X-ray powder diffraction up to 7.67 GPa with a diamond anvil cell under hydrostatic conditions. No phase transition has been observed within the pressure range investigated. Axial and volume isothermal Equations of State (EoS) of vanadinite were determined. Fitting the P–V data with a third-order Birch-Murnaghan (BM) EoS, using the data weighted by the uncertainties in P and V, we obtained: V ₀ = 681(1) Å³, K ₀ = 41(5) GPa, and K′ = 12.5(2.5). The evolution of the lattice constants with P shows a strong anisotropic compression pattern. The axial bulk moduli were calculated with a third-order “linearized” BM-EoS. The EoS parameters are: a ₀ = 10.3302(2) Å, K ₀(a) = 35(2) GPa and K′(a) = 10(1) for the a-axis; c ₀ = 7.3520(3) Å, K ₀(c) = 98(4) GPa, and K′(c) = 9(2) for the c-axis (K ₀(a):K ₀(c) = 1:2.80). Axial and volume Eulerian-finite strain (fe) at different normalized stress (Fe) were calculated. The weighted linear regression through the data points yields the following intercept values: Fe _a (0) = 35(2) GPa for the a-axis, Fe _c (0) = 98(4) GPa for the c-axis and Fe _V (0) = 45(2) GPa for the unit-cell volume. The slope of the regression lines gives rise to K′ values of 10(1) for the a-axis, 9(2) for the c-axis and 11(1) for the unit cell-volume. A comparison between the HP-elastic response of vanadinite and the iso-structural apatite is carried out. The possible reasons of the elastic anisotropy are discussed.     

Characterization of flow and transport dynamics in karst aquifers by analyzing tracer test results in conduits and recharge areas (the Egino Massif,Basque Country,Spain): environmental and management implications

Karst aquifers contribute to supplying drinking water to almost a quarter of the world´s population. Their complex dynamics requires specific approaches aimed at recognizing their singularities, analyzing its vulnerability, and ensuring water resources quality. In this paper, the results of processing and modeling five breakthrough tracer curves obtained under different hydrodynamic conditions in the main conduit of Egino karst aquifer (Basque Country, Spain) are analyzed together with those involving pressure injections of the tracer in the saturated zone of the karst massif recharge area. In the conduit, transport is immediate and highly efficient (recovery rates above 84% and dispersion coefficients from 15.04 to 84.35 m²/min); tracer retentions increase as flow rates decrease and no significant contributions to its surroundings are observed. In contrast, tracer transport from the massif recharge area is more complex: after injection at a pressure of 1 MPa, most tracer remains in the surrounding of the injection borehole, retained in a saturated medium of low effective fracture porosity (? _f ?=?1.02?×?10^?4, assuming a radial divergent flow model); subsequently, the main tracer mobilization to the spring was registered with the first rains, with 0.088 m/min mean velocity and high concentrations per unit mass being injected (C _p /M₀?=?0.03 mg/L/kg), which is evidence that the tracer reaches soon the karst conduit network. In any case, a decreasing tracer presence is registered at the injection zone during a hydrological cycle. In both cases, the observed non-linearity of transport processes should be considered in the development of vulnerability approaches, modeling efforts, and mapping. Furthermore, in the case of karst massif recharge areas, as the presence of pollutants may have a significant impact on the springs and persist over time, their management and protection needs must be revised in each specific site. Simultaneously, quality-monitoring programs at the springs must be adapted to the aquifers recognized dynamics.     

Identifying groundwater recharge connections in the Moscow (USA) sub-basin using isotopic tracers and a soil moisture routing model

Globally, aquifers are suffering from large abstractions resulting in groundwater level declines. These declines can be caused by excessive abstraction for drinking water, irrigation purposes or industrial use. Basaltic aquifers also face these conflicts. A large flood basalt area (1.1?×?10⁵ km²) can be found in the Northwest of the USA. This Columbia River Basalt Group (CRBG) consists of a thick series of basalt flows of Miocene age. The two major hydrogeological units (Wanapum and Grand Ronde formations) are widely used for water abstraction. The mean decline over recent decades has been 0.6 m year^?1. At present day, abstraction wells are drying up, and base flow of rivers is reduced. At the eastern part of CRBG, the Moscow sub-basin on the Idaho/Washington State border can be found. Although a thick poorly permeable clay layer exists on top of the basalt aquifer, groundwater level dynamics suggest that groundwater recharge occurs at certain locations. A set of wells and springs has been monitored bi-weekly for 9 months for δ¹⁸O and δ²H. Large isotopic fluctuations and d-excess values close to the meteoric water line in some wells are indicating that recharge occurs at the granite/basalt interface through lateral flow paths in and below the clay. A soil moisture routing (SMR) model showed that most recharge occurs on the granitic mountains. The basaltic aquifer receives recharge from these sedimentary zones around the granite/basalt interface. The identification of these types of areas is of major importance for future managed-aquifer recharge solutions to solve problems of groundwater depletion.