Complementary use of tracer and pumping tests to characterize a heterogeneous channelized aquifer system in New Zealand

The combined use of pumping and tracer test data enabled the derivation of equivalent average hydraulic conductivities (K_avg) for each test in a heterogeneous channelized alluvial aquifer, whereas K values of the preferential flow paths were two orders of magnitude higher. Greater and earlier drawdown was generally observed along preferential flow lines in a pumping test, within an array of 21 wells. The study aim was to characterize hydraulic properties of a channelized aquifer system in New Zealand by combining tracer and pumping test data. Estimates were able to be made of the percentage of highly permeable channels within the profile (～1.2%), effective porosity that reflected the maximum fraction of highly permeable channels within the aquifer (?_eff–pc ～0.0038), and flows through highly permeable channels (～98%) and the sandy gravel matrix material (～2%). Using ?_eff–pc, a tracer test K_avg value (～93 m/day) was estimated that was equivalent to pumping test values (～100 m/day), but two orders of magnitude smaller than K calculated solely from transport through permeable channels (K_pc ～8,400 m/day). Derived K values of permeable and matrix material were similar to values derived from grain size distribution using the Kozeny-Carman equation.     

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.     

Spatial information of soil hydraulic conductivity and performance of cokriging over kriging in a semi-arid basin scale

Unlike the studies in small parcels by systematic measurements, the spatial variability of soil properties is expected to increase in those over relatively large areas or scales. Spatial variability of soil hydraulic conductivity (K _h) is of significance for the environmental processes, such as soil erosion, plant growth, transport of the plant nutrients in a soil profile and ground water levels. However, its variability is not much and sufficiently known at basin scale. A study of testing the performance of cokriging of K _h compared with that of kriging was conducted in the catchment area of Sarayköy II Irrigation Dam in Cank?r?, Turkey. A total of 300 soil surface samples (0–10 cm) were collected from the catchment with irregular intervals. Of the selected soil properties, because the water-stable aggregates (WSA) indicated the highest relationship with the hydraulic conductivity by the Pearson correlation analysis, it is used as an auxiliary variable to predict K _h by the cokriging procedure. In addition, the sampling density was reduced randomly to n = 175, n = 150, n = 75 and n = 50 for K _h to determine if the superiority of cokriging over kriging would exist. Statistically, the results showed that all reduced K _h was as good as the complete K _h when its auxiliary relations with WSA were used in cokriging. Particularly, the results of the “Relative Reduction in MSE” (RMSE) revealed that the reduced data set of n = 75 produced the most accurate map than the others. In this basin-scaled study, there was a clear superiority of the cokriging procedure by the reduction in data although a very undulating topography and topographically different aspects, two different land uses with non-uniform vegetation density, different parent materials and soil textures were present in the area. Hence, using the statistically significant auxiliary relationship between K _h and WSA might bring about a very useful data set for watershed hydrological researches.     

Hydrogeologic and hydraulic characterization of aquifer and nonaquifer layers in a lateritic terrain (West Bengal,India)

The hydrogeologic and hydraulic characteristics of a lateritic terrain in West Bengal, India, were investigated. Test drilling was conducted at ten sites and grain-size distribution curves (GSDCs) were prepared for 275 geologic samples. Performance evaluation of eight grain-size-analysis (GSA) methods was carried out to estimate the hydraulic conductivity (K) of subsurface formations. Finally, the GSA results were validated against pumping-test data. The GSDCs indicated that shallow aquifer layers are coarser than the deeper aquifer layers (uniformity coefficient 0.19–11.4). Stratigraphy analysis revealed that both shallow and deep aquifers of varying thickness exist at depths 9–40 and 40–79 m, respectively. The mean K estimates by the GSA methods are 3.62–292.86 m/day for shallow aquifer layers and 0.97–209.93 m/day for the deeper aquifer layers, suggesting significant aquifer heterogeneity. Pumping-test data indicated that the deeper aquifers are leaky confined with transmissivity 122.69–693.79 m²/day, storage coefficient 1.01?×?10^?7–2.13?×?10^?4 and leakance 2.01?×?10^?7–34.56?×?10^?2 day^?1. Although the K values yielded by the GSA methods are generally larger than those obtained from the pumping tests, the Slichter, Harleman and US Bureau Reclamation (USBR) GSA methods yielded reasonable values at most of the sites (1–3 times higher than K estimates by the pumping-test method). In conclusion, more reliable aquifers exist at deeper depths that can be tapped for dependable water supply. GSA methods such as Slichter, Harleman and USBR can be used for the preliminary assessment of K in lateritic terrains in the absence of reliable field methods.     

Impact of pore water conductivity and porosity on the electrical conductivity of kaolinite

The purpose of this study is to quantify the magnitudes of surface conduction and pore water conduction from the measured electrical conductivity of kaolinite, with the ultimate goal of estimating the electrical conductivity of kaolinite with a wide range of pore water conductivities (σ _w = 0.013–3.356 S/m) and porosities (n = 0.368–1.0). Therefore, the theoretical background of the electrical conductivity in soils was reviewed, and electrical conductivity measurements on kaolinite were performed using both slurry and consolidation tests in this study. The results of this study demonstrate that the variations of measured electrical conductivity (σ _mix) with n are debatable according to the values of σ _w, because a decrease in n results in both an increase in surface conduction (K _s) and a decrease in pore water conduction (K _w); this causes the relative magnitude of K _s compared to that of K _w to vary with σ _w and n. Consequently, this study develops the relation between the porosity-normalized K _s/K _w and 1/σ _w. Additionally, the surface conductivity of the tested kaolinite is back-calculated and compared with the previous relationship between K _s and zeta potential of kaolinite. The measured and estimated σ _mix values are compared with the varying pore water conductivity and porosity values.     

Upscaling of transmissivity,derived from specific capacity: a hydrogeomorphological approach applied to the Doon Valley aquifer system in India

Transmissivity (T) is one of the most important parameters in groundwater studies, and is generally estimated from pumping tests. T can also be deduced from abundantly available specific-capacity (Q/s) data by using analytical and/or empirical approaches, further upscaled by geostatistical methods. A different, remote sensing based, hydrogeomorphological approach is proposed, to upscale T from point- or well-scale to aquifer-scale, and it is applied to the piedmont alluvial aquifer system of Doon Valley in India. In the first step, Q/s and T data-pairs available from aquifer tests were used to establish an empirical, logarithmic relation. Subsequently, satellite imagery along with available data from published and unpublished maps, literature sources and ground surveys were used to divide the aquifer system into major hydrogeomorphological domains that control the groundwater occurrence and flow. Then, the T data derived from Q/s (using the empirical relation) and that available from pumping tests at well-scale were upscaled to aquifer-scale by averaging the T values within each hydrogeomorphological domain. Such a stratification approach is especially useful in areas where availability of only a few data-pairs of known Q/s and T limit the use of geostatistical techniques. A comparative study of the published empirical relations between Q/s and T in various hydrogeologic settings revealed that the relation obtained for Doon Valley aquifer system is close to that found for a similar alluvial aquifer system in Morocco.     

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.     

Toward reliable calibration of aquifer hydrodynamic parameters: characterizing and optimization of arid groundwater system using swarm intelligence optimization algorithm

Despite advanced development in computational techniques, the issue of how to adequately calibrate and minimize misfit between system properties and corresponding measurements remains a challenging task in groundwater modeling. Two important features of the groundwater regime, hydraulic conductivity (k) and specific yield (S _y), that control aquifer dynamic vary spatially within an aquifer system due to geologic heterogeneity. This paper provides the first attempt in using an advanced swarm-intelligence-based optimization algorithm (cuckoo optimization algorithm, COA) coupled with a distributed hydrogeology model (i.e., MODFLOW) to calibrate aquifer hydrodynamic parameters (S _y and k) over an arid groundwater system in east Iran. Our optimization approach was posed in a single-objective manner by the trade-off between sum of absolute error and the adherent swarm optimization approach. The COA optimization algorithm further yielded both hydraulic conductivity and specific yield parameters with high performance and the least error. Estimation of depth to water table revealed skillful prediction for a set of cells located at the middle of the aquifer system whereas showed unskillful prediction at the headwater due to frequent water storage changes at the inflow boundary. Groundwater depth reduced from east toward west and southwest parts of the aquifer because of extensive pumping activities that caused a smoothening influence on the shape of the simulated head curve. The results demonstrated a clear need to optimize arid aquifer parameters and to compute groundwater response across an arid region.     

Assessment of textural variational pattern and electrical conduction of economic and accessible Quaternary hydrolithofacies via geoelectric and laboratory methods in SE Nigeria: A case study of select locations in Akwa Ibom State

Textural variational pattern of economic and accessible Quaternary aquifer repositories and its conductivity in the south-eastern Nigeria have been assessed through the integration of vertical electrical sounding and laboratory measurements. The results have shown the lithological attributes, pore-water and amount of residual clay minerals in the assumed clean sand; mechanism of charge fixation at the fluid - surface interface; intricate geometry of pores and pore channels; formation’s ability to transmit pore-water and cation exchange capacity.The connections of electrical and hydraulic properties and their distributions have been established. The average interface conductivity contributed by residual clay minerals in assumed clean sands of the aquifer repositories in the study area have been estimated as 30µS/m. Intrinsic average porosity and formation factor have been respectively deduced as 12% and 14.75. Comparing the simulated aquifer formation factor obtained from the observed porosity data with the observed aquifer formation factor, indicates the that study area has 0.5 ≤ a ≤ 0.8 pore geometry factor and 1.5 ≤ m ≤ 2.0 cementation factor as the best fitting values. The interrelations between aquifer parameters have been established through different plots and the aquifer have been empirically proved to be associated with residual clay minerals as the interface conductivity C_q is not equal to zero. The wide ranges of parameters estimated are an indication of variations in grain size. The estimated intrinsic average porosity, formation factor and the average BQ_v are viable in characterizing the aquifer flow dynamics and contaminant modelling in the associated aquifer sands For low pore geometry factors a (0.2) and low cementation factor m (0.5) the formation factor remains fairly constant. However, marked variability is noticed at higher a (1.0) and m (2.5). Despite the observed variability in formation factors at the indicated porosities, the spatial or geometrical spread of the formation factor remains unchanged in the aquifer units. The Tables for geoelectric and petrophysical parameters and the associated mathematical models generated in this study can be used for groundwater contaminant modelling and simulation of pore space parameters with reasonable accuracy.     

Impact of model complexity and multi-scale data integration on the estimation of hydrogeological parameters in a dual-porosity aquifer

This study investigates the impact of model complexity and multi-scale prior hydrogeological data on the interpretation of pumping test data in a dual-porosity aquifer (the Chalk aquifer in England, UK). In order to characterize the hydrogeological properties, different approaches ranging from a traditional analytical solution (Theis approach) to more sophisticated numerical models with automatically calibrated input parameters are applied. Comparisons of results from the different approaches show that neither traditional analytical solutions nor a numerical model assuming a homogenous and isotropic aquifer can adequately explain the observed drawdowns. A better reproduction of the observed drawdowns in all seven monitoring locations is instead achieved when medium and local-scale prior information about the vertical hydraulic conductivity (K) distribution is used to constrain the model calibration process. In particular, the integration of medium-scale vertical K variations based on flowmeter measurements lead to an improvement in the goodness-of-fit of the simulated drawdowns of about 30%. Further improvements (up to 70%) were observed when a simple upscaling approach was used to integrate small-scale K data to constrain the automatic calibration process of the numerical model. Although the analysis focuses on a specific case study, these results provide insights about the representativeness of the estimates of hydrogeological properties based on different interpretations of pumping test data, and promote the integration of multi-scale data for the characterization of heterogeneous aquifers in complex hydrogeological settings.     

Time constant of hydraulic-head response in aquifers subjected to sudden recharge change: application to large basins

Analytical formulae are proposed to describe the first-order temporal evolution of the head in large groundwater systems (such as those found in North Africa or eastern Australia) that are subjected to drastic modifications of their recharge conditions (such as those in Pleistocene and Holocene times). The mathematical model is based on the hydrodynamics of a mixed-aquifer system composed of a confined aquifer connected to an unconfined one with a large storage capacity. The transient behaviour of the head following a sudden change of recharge conditions is computed with Laplace transforms for linear one-dimensional and cylindrical geometries. This transient evolution closely follows an exponential trend exp(?t/τ). The time constant τ is expressed analytically as a function of the various parameters characterizing the system. In many commonly occurring situations, τ depends on only four parameters: the width a _c of the main confined aquifer, its transmissivity T _c, the integrated storage situated upstream in the unconfined aquifer M?=?S _u a _u, and a curvature parameter accounting for convergence/divergence effects. This model is applied to the natural decay of large aquifer basins of the Sahara and Australia following the end of the mid-Holocene humid period. The observed persistence of the resource is discussed on the basis of the time constant estimated with the system parameters. This comparison confirms the role of the upstream water reserve, which is modelled as an unconfined aquifer, and highlights the significant increase of the time constant in case of converging flow.     

Estimating groundwater evapotranspiration rates using diurnal water-table fluctuations in a semi-arid riparian zone

In semi-arid climates, phreatophytes draw on shallow aquifers, and groundwater evapotranspiration (ET_G) is a principal component of groundwater budgets. Diurnal water table fluctuations, which often are a product of ET_G, were monitored in the riparian zone of Red Canyon Creek, Wyoming, USA. These fluctuations were higher in a riparian wetland (2–36 mm) than a grass-covered meadow (1–6 mm). The onset and cessation of water-table fluctuations correspond to daily temperatures relative to freezing. Spatial differences were due to vegetation type and specific yield, while temporal changes were due to vegetation dormancy. Ratios of ET_G to potential evapotranspiration (PET), K _c,GW, were similar to ratios of actual evapotranspiration (ET) to PET, K _c, in semi-arid rangelands. Before vegetation senescence, K _c,GW increased between precipitation events, suggesting phreatophytes pull more water from the saturated zone as soil moisture decreases. In contrast, K _c decreases with soil moisture following precipitation events as ET becomes increasingly water-limited. Error in ET_G is primarily from estimates of specific yield (S _y), which is difficult to quantify in heterogeneous sediments. ET_G values may be more reliable because the range of acceptable S _y is smaller than K _c and S _y does not change with vegetation type or soil moisture.     

Nonlinear consolidation in randomly heterogeneous highly compressible aquitards

Severe land subsidence due to groundwater extraction may occur in multiaquifer systems where highly compressible aquitards are present. The highly compressible nature of the aquitards leads to nonlinear consolidation where the groundwater flow parameters are stress-dependent. The case is further complicated by the heterogeneity of the hydrogeologic and geotechnical properties of the aquitards. The effect of realistic vertical heterogeneity of hydrogeologic and geotechnical parameters on the consolidation of highly compressible aquitards is investigated by means of one-dimensional Monte Carlo numerical simulations where the lower boundary represents the effect of an instant drop in hydraulic head due to groundwater pumping. Two thousand realizations are generated for each of the following parameters: hydraulic conductivity (K), compression index (C _c), void ratio (e) and m (an empirical parameter relating hydraulic conductivity and void ratio). The correlation structure, the mean and the variance for each parameter were obtained from a literature review about field studies in the lacustrine sediments of Mexico City. The results indicate that among the parameters considered, random K has the largest effect on the ensemble average behavior of the system when compared to a nonlinear consolidation model with deterministic initial parameters. The deterministic solution underestimates the ensemble average of total settlement when initial K is random. In addition, random K leads to the largest variance (and therefore largest uncertainty) of total settlement, groundwater flux and time to reach steady-state conditions.     

Automatic estimation of aquifer parameters using long-term water supply pumping and injection records

Analyses are presented of long-term hydrographs perturbed by variable pumping/injection events in a confined aquifer at a municipal water-supply well field in the Region of Waterloo, Ontario (Canada). Such records are typically not considered for aquifer test analysis. Here, the water-level variations are fingerprinted to pumping/injection rate changes using the Theis model implemented in the WELLS code coupled with PEST. Analyses of these records yield a set of transmissivity (T) and storativity (S) estimates between each monitoring and production borehole. These individual estimates are found to poorly predict water-level variations at nearby monitoring boreholes not used in the calibration effort. On the other hand, the geometric means of the individual T and S estimates are similar to those obtained from previous pumping tests conducted at the same site and adequately predict water-level variations in other boreholes. The analyses reveal that long-term municipal water-level records are amenable to analyses using a simple analytical solution to estimate aquifer parameters. However, uniform parameters estimated with analytical solutions should be considered as first rough estimates. More accurate hydraulic parameters should be obtained by calibrating a three-dimensional numerical model that rigorously captures the complexities of the site with these data.     

Analytical expressions for three-phase generalized relative permeabilities in water- and oil-wet capillary tubes

We analyze three-phase flow of immiscible fluids taking place within an elementary capillary tube with circular cross-section under water- and oil-wet conditions. We account explicitly for momentum transfer between the moving phases, which leads to the phenomenon of viscous coupling, by imposing continuity of velocity and shear stress at fluid-fluid interfaces. The macroscopic flow model which describes the system at the Darcy scale includes three-phase effective relative permeabilities, K _{i j,r} , accounting for the flux of the ith phase due to the presence of the jth phase. These effective parameters strongly depend on phase saturations, fluid viscosities, and wettability of the solid matrix. In the considered flow setting, K _{i j,r} reduce to a set of nine scalar quantities, K _{i j,r} . Our results show that K _{i j,r} of the wetting phase is a function only of the fluid phase own saturation. Otherwise, K _{i j,r} of the non-wetting phase depends on the saturation of all fluids in the system and on oil and water viscosities. Viscous coupling effects (encapsulated in K _{i j,r} with i ≠ j) can be significantly relevant in both water- and oil-wet systems. Wettability conditions influence oil flow at a rate that increases linearly with viscosity ratio between oil and water phases.     

Non-Darcian flow to a partially penetrating well in a confined aquifer with a finite-thickness skin

Non-Darcian flow to a partially penetrating well in a confined aquifer with a finite-thickness skin was investigated. The Izbash equation is used to describe the non-Darcian flow in the horizontal direction, and the vertical flow is described as Darcian. The solution for the newly developed non-Darcian flow model can be obtained by applying the linearization procedure in conjunction with the Laplace transform and the finite Fourier cosine transform. The flow model combines the effects of the non-Darcian flow, partial penetration of the well, and the finite thickness of the well skin. The results show that the depression cone spread is larger for the Darcian flow than for the non-Darcian flow. The drawdowns within the skin zone for a fully penetrating well are smaller than those for the partially penetrating well. The skin type and skin thickness have great impact on the drawdown in the skin zone, while they have little influence on drawdown in the formation zone. The sensitivity analysis indicates that the drawdown in the formation zone is sensitive to the power index (n), the length of well screen (w), the apparent radial hydraulic conductivity of the formation zone (K _r2), and the specific storage of the formation zone (S _s2) at early times, and it is very sensitive to the parameters n, w and K _r2 at late times, especially to n, while it is not sensitive to the skin thickness (r _s).     

High-pressure thermo-elastic properties of beryl (Al<Subscript>4</Subscript>Be<Subscript>6</Subscript>Si<Subscript>12</Subscript>O<Subscript>36</Subscript>) from ab initio calculations,and observations about the source of thermal expansion

Ab initio calculations of thermo-elastic properties of beryl (Al₄Be₆Si₁₂O₃₆) have been carried out at the hybrid HF/DFT level by using the B3LYP and WC1LYP Hamiltonians. Static geometries and vibrational frequencies were calculated at different values of the unit cell volume to get static pressure and mode-γ Grüneisen’s parameters. Zero point and thermal pressures were calculated by following a standard statistical-thermodynamics approach, within the limit of the quasi-harmonic approximation, and added to the static pressure at each volume, to get the total pressure (P) as a function of both temperature (T) and cell volume (V). The resulting P(V, T) curves were fitted by appropriate EoS’, to get bulk modulus (K ₀) and its derivative (K′), at different temperatures. The calculation successfully reproduced the available experimental data concerning compressibility at room temperature (the WC1LYP Hamiltonian provided K ₀ and K′ values of 180.2 Gpa and 4.0, respectively) and the low values observed for the thermal expansion coefficient. A zone-centre soft mode \( P6/mcc \to P\bar{1} \) phase transition was predicted to occur at a pressure of about 14 GPa; the reduction of the frequency of the soft vibrational mode, as the pressure is increased, and the similar behaviour of the majority of the low-frequency modes, provided an explanation of the thermal behaviour of the crystal, which is consistent with the RUM model (Rigid Unit Model; Dove et al. in Miner Mag 59:629–639, 1995), where the negative contribution to thermal expansion is ascribed to a geometric effect connected to the tilting of rigid polyhedra in framework silicates.     

Investigation into the change of the porosity parameter during transition from laboratory to reservoir conditions

The purpose of the investigation is to reveal the dependences of P _p = f(C _p) on reservoir conditions and the lithological composition of rocks. The samples were studied using a set of lithological-petrographic investigations. To obtain the dependences on thermobaric conditions, 90 samples of different porosities (C _p) of 15, 20, and 25% were collected. As a result, the general pattern of the change in the rock resistivity during the transition from atmospheric conditions of measuring to reservoir conditions was established. Dependences of porosity parameter P _p on porosity coefficient C _p were obtained for three values of formation water salinity and three reservoir conditions. The measurement errors of the porosity parameter P _p were calculated using dependences obtained under atmospheric conditions.     

Modelling the effects of porous media deformation on the propagation of water-table waves in a sandy unconfined aquifer

This paper examines the influence of porous media deformation on water-table wave dispersion in an unconfined aquifer using a numerical model which couples Richards’ equation to the poro-elastic model. The study was motivated by the findings of Shoushtari et al. (J Hydrol 533:412–440, 2016) who were unable to reproduce the observed wave dispersion in their sand flume data with either numerical Richards’ equation models (assuming rigid porous media) or existing analytic solutions. The water-table wave dispersion is quantified via the complex wave number extracted from the predicted amplitude and phase profiles. A sensitivity analysis was performed to establish the influence of the main parameters in the poro-elastic model, namely Young’s modulus (E) and Poisson’s ratio (ν). For a short oscillation period (T?=?16.4 s), the phase lag increase rate (k _i) is sensitive to the chosen values of E and ν, demonstrating an inverse relationship with both parameters. Changes in the amplitude decay rate (k _r), however, were negligible. For a longer oscillation period (T?=?908.6 s), variations in the values of E and ν resulted in only small changes in both k _r and k _i. In both the short and long period cases, the poro-elastic model is unable to reproduce the observed wave dispersion in the existing laboratory data. Hence porous media deformation cannot explain the additional energy dissipation in the laboratory data. Shoushtari SMH, Cartwright N, Perrochet P, Nielsen P (2016) The effects of oscillation period on groundwater wave dispersion in a sandy unconfined aquifer: sand flume experiments and modelling. J Hydrol 533:412–440.     

Compressibility of strontium orthophosphate Sr<Subscript>3</Subscript>(PO<Subscript>4</Subscript>)<Subscript>2</Subscript> at high pressure

High pressure in situ synchrotron X-ray diffraction experiment of strontium orthophosphate Sr₃(PO₄)₂ has been carried out to 20.0 GPa at room temperature using multianvil apparatus. Fitting a third-order Birch–Murnaghan equation of state to the P–V data yields a volume of V ₀ = 498.0 ± 0.1 Å³, an isothermal bulk modulus of K _T  = 89.5 ± 1.7 GPa, and first pressure derivative of K _T ′ = 6.57 ± 0.34. If K _T ′ is fixed at 4, K _T is obtained as 104.4 ± 1.2 GPa. Analysis of axial compressible modulus shows that the a-axis (K _a  = 79.6 ± 3.2 GPa) is more compressible than the c-axis (K _c  = 116.4 ± 4.3 GPa). Based on the high pressure Raman spectroscopic results, the mode Grüneisen parameters are determined and the average mode Grüneisen parameter of PO₄ vibrations of Sr₃(PO₄)₂ is calculated to be 0.30(2).