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.     

Analysis of magnetic activity of the rapidly rotating stars He 373 and AP 225

Spectroscopic and photometric data for the two rapidly rotating members of the α Persei cluster He 373 and AP 225 are analyzed. Improved estimates have been obtained for the projected equatorial rotation velocities: v sin i = 164 km/s for He 323 and v sin i = 129 km/s for AP 225. Multi-band photometric mapping is used to map the spot distributions on the surfaces of the two stars. The fractional spotted areas S and mean temperature difference ΔT between the unspotted photosphere and the spots are estimated (S = 7% and ΔT = 1000 K for He 373; S = 9% and ΔT = 800 K for AP 225). The H _α line profiles of both stars have variable emission components whose widths are used to deduce the presence of extended regions of emission reaching the corotation radius.     

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.     

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.     

Spots and the activity of M dwarfs from observations with the Kepler Space Telescope

A method for estimating the spottedness parameter S (the spotted area as a fraction of the surface of an active star) proposed earlier is applied to an analysis of activity in 1570 M dwarf stars. The analysis is based on observational material obtained with the Kepler Space Telescope, as well as data on the fluxes of the studied objects in the near and far ultraviolet (NUV and FUV) based on data from the GALEX space telescope. The variations of S with the ages of the stars are studied (four groups with different ages are distinguished), as well as variations of S with their rotational periods. A diagram characterizing the relationship between S and the Rossby number Ro resembles the classical dependence of the X-ray luminosities of active stars on Ro, and a saturation regime is attained at the same value, Ro = 0.13. Moreover, objects with ages of more than 100 million years do not form a single sequence (and stars older than 900 million years possess surface spottednesses of order 1%). The S?Ro dependence obtained could expand possibilities for analyzing the dependence of the X-ray luminosities of active stars on their Rossby numbers, and could also be applied to refine parameters characterizing the action of dynamo mechanisms, such as the dynamo number N _D . A comparison of the GALEX NUV and FUV brightness estimates with the activity parameters of the stars suggests that younger, more rapidly rotating active stars are brighter in the NUV, and that the FUV flux grows and the difference of the FUV and NUV brightnesses decreases with increasing spottedness S.     

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.     

A comparative study of hydraulic conductivity estimations using geostatistics

Three approaches for estimating the hydraulic conductivity (K) of the Trifa aquifer, Morocco were investigated: (1) kriging of the K values obtained from pumping tests, (2) cokriging of the pumping test data with electrical resistivity data as a secondary variable, and (3) cokriging of the pumping test data with the slope of the water table. Gauss-transformed values of the variables are used because they provide more robust variograms and transformed values of the primary and secondary variables show correlations higher than the raw values, which is beneficial in cokriging. In cokriging with electrical resistivity, two zones are considered since the geological deposits are different from the north to the south of the aquifer, which is reflected in different correlations between the variables. Comparison of the three approaches is based mainly on the estimation errors, and to a lesser degree on the cross-validations of the corresponding variogram models and general considerations, like the measurements’ reliability and aquifer make-up. The best-estimated K is given by cokriging with the slope of the water table and is therefore preferred for further use in groundwater flow modeling. Thus, electrical resistivity or the slope of the water table can both be used as secondary variables to estimate K, especially in heterogeneous aquifers with lateral variations in lithology, as is the case of the Trifa aquifer.     

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.     

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.     

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.     

Spots and the Activity of Stars in the Hyades Cluster from Observations with the Kepler Space Telescope (K2)

Observations of the K2 mission (continuing the program of the Kepler Space Telescope) are used to estimate the spot coverage S (the fractional area of spots on the surface of an active star) for stars of the Hyades cluster. The analysis is based on data on the photometric variations of 47 confirmed single cluster members, together with their atmospheric parameters, masses, and rotation periods. The resulting values of S for these Hyades objects are lower than those stars of the Pleiades cluster (on average, by ΔS ~ 0.05?0.06). A comparison of the results of studies of cool, low-mass dwarfs in the Hyades and Pleiades clusters, as well as the results of a study of 1570 M stars from the main field observed in the Kepler SpaceMission, indicates that the Hyades stars are more evolved than the Pleiades stars, and demonstrate lower activity. The activity of seven solar-type Hyades stars (S = 0.013 ± 0.006) almost approaches the activity level of the present-day Sun, and is lower than the activity of solar-mass stars in the Pleiades (S = 0.031 ± 0.003). Solar-type stars in the Hyades rotate faster than the Sun (〈P〉 = 8.6 ^d ), but slower than similar Pleiades stars.     

Spots and activity of Pleiades stars from observations with the Kepler Space Telescope (K2)

Observations of the K2 continuation of Kepler Space Telescope program are used to estimate the spot coverage S (the fractional spotted area on the surface of an active star) for stars of the Pleiades cluster. The analysis is based on data on photometric variations of 759 confirmed clustermembers, together with their atmospheric parameters, masses, and rotation periods. The relationship between the activity (S) of these Pleiades stars and their effective temperatures shows considerable change in S for stars with temperatures T _eff less than 6100 K (this can be considered the limiting value for which spot formation activity begins) and a monotonic increase in S for cooler objects (a change in the slope for stars with Teff ~ 3700 K). The scatter in this parameter ΔS about its mean dependence on the (V ?Ks)0 color index remains approximately the same over the entire (V?K _s )0 range, including cool, fully convective dwarfs. The computated S values do not indicate differences between slowly rotating and rapidly rotating stars with color indices 1.1 < (V?K _s )0 < 3.7. The main results of this study include measurements of the activity of a large number of stars having the same age (759 members of the Pleiades cluster), resulting in the first determination of the relationship between the spot-forming activity and masses of stars. For 27 stars with masses differing from the solarmass by nomore than 0.1M⊙, themean spot coverage is S = 0.031±0.003, suggesting that the activity of candidate young Suns is more pronounced than that of the present-day Sun. These stars rotate considerably faster than the Sun, with an average rotation period of 4.3^d. The results of this study of cool, low-mass dwarfs of the Pleiades cluster are compared to results from an earlier study of 1570 M stars.     

The impact of reflection effects on the parameters of the old pre-cataclysmic variables MS Peg and LM Com

We have determined the main parameters of the old precataclysmic variable stars MS Peg and LM Com. The radial velocities of the components, reflection effects in the spectra, and light curves of the systems are studied based on model stellar atmospheres subject to external irradiation. Forty-seven moderate-resolution spectra for MS Peg and 57 for LM Com obtained with the 6-m telescope of the Special Astrophysical Observatory are used to derive the refined orbital periods of 0.1736660 days and 0.2586873 days, respectively; the orbital eccentricities do not exceed e=0.04. The mass (M _w =0.49e_⊙) and radius (e _w =0.015R_⊙) of the MS Peg primary calculated using the gravitational redshift correspond to those for a cooling carbon white dwarf with a thin hydrogen envelope. The parameters of the red dwarf (M _r =0.19M_⊙, T_eff=3560 K, R _r =0.18R_⊙) are close to those derived from evolutionary tracks for main-sequence M stars with solar chemical composition. The radius (R _r =0.22R_⊙) and temperature (T_eff=3650 K) of the LM Com secondary exceed theoretical estimates for main-sequence stars with masses of M _r =0.17M_⊙. The luminosity excess of the red dwarf in LM Com can be explained by a prolonged (T>5×10⁶ yrs) relaxation of the M star to its normal state after the binary leaves the common-envelope stage. For both systems, theoretical U, B, V, and R light curves and spectra calculated using the adopted sets of parameters are generally consistent with the observations. This confirms the radiative origin of the hot spots, the unimportance of horizontal radiative transport, and the absence of large-scale velocity fields with high values (V_trans>50 km/s) at the surfaces of the secondaries. Most of the emission lines in the spectra of these objects are formed under conditions close to thermalization, enabling modeling of their pro files in an LTE approximation. A strong λ3905 Å emission line has been identified as the 3s²3p4s 1P⁰-3s²3p² 1S SiI λ3905.52 Å line formed in the atmosphere of the hot spot. The observed intensity can be explained by non-LTE “superionization” of SiI atoms by soft UV radiation from the white dwarf. We suggest a technique for identifying binaries whose cool components are subject to UV irradiation based on observations of λ3905 Å emission in their spectra.     

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.     

The character of scaling earthquake source spectra for Kamchatka in the 3.5–6.5 magnitude range

The properties of the source spectra of local shallow-focus earthquakes on Kamchatka in the range of magnitudes M _w = 3.5–6.5 are studied using 460 records of S-waves obtained at the PET station. The family of average source spectra is constructed; the spectra are used to study the relationship between M _w and the key quasi-dimensionless source parameters: stress drop Δσ and apparent stress σ_a. It is found that the parameter Δσ is almost stable, while σ_a grows steadily as the magnitude M _w increases, indicating that the similarity is violated. It is known that at sufficiently large M _w the similarity hypothesis is approximately valid: both parameters Δσ and σ_a do not show any noticeable magnitude dependence. It has been established that M _w ≈ 5.7 is the threshold value of the magnitude when the change in regimes described occurs for the conditions on Kamchatka.     

Spots and activity of solar-type stars from Kepler observations

The spot coverages S for 2846 solar-type stars with effective temperatures from 5700 K to 5800 K and gravities from 4.4 to 4.5 have been measured. An analysis based on the MAST catalog, which presents photometric measurements obtained with the Kepler Space Telescope during Q9 is presented. The existence of two groups of solar-type stars, with S values between 0.001 and 0.007 and with S > 0.007, is inferred. The second group (active stars) contains 279 stars (about 10% of the total number of stars analyzed). The mean S parameter for the entire sample is 0.004, comparable to the mean spot coverage of the Sun. In general, the dependence of S on the rotation period for solar-type stars has characteristics similar to those found earlier for stars with exoplanets. For the vast majority of the stars in the sample, the activity is constant, and independent of age. The activity of the small number of active stars with S > 0.007 decreases with age. The age variations of the chromospheric activity index R′_HK are compared to variations of the spot coverage S. The relations analyzed have common characteristic features. It is likely that both the spot activity level and the chromospheric activity level abruptly decrease for stars older than 4 billion yrs.     

Long-term Spot-Coverage Variations of 13 BY Dra G–K Dwarfs

The results of spot-coverage modeling for 13 active G–K dwarf stars based on many-year photometric observations are presented. The results of UBVRI observations of eight stars performed at the Crimean Astrophysical Observatory were used together with data from the literature in this analysis. The spot-coverage parameters for 13 selected BY Dra active red dwarfs have been redetermined to improve the zonal spot-coverage model for the stellar photospheres, which currently allows for the presence of two active longitudes. Time variations of the spot-activity characteristics of these systems were analyzed with the aim of searching for possible cyclic variations. All the stars, with the exception of OU Gem and BE Cet, show fairly strong correlations between variations in the spot latitudes and spot areas, with absolute values of the correlation coefficients, R(〈?〉, S), ranging from 0.38 to 0.92. For five stars, an anti-correlation between the mean latitude and area of the spots was found (R(〈?〉, S) from–0.24 to–0.73). This behavior may reflect the drift of spots toward the equator in the course of their development. Eight stars feature positive correlations, i.e. the spots drift towards the pole as their areas increase. Nine stars demonstrate activity cycles, which are reflected in photometric variations as well as variations of the spot areas and mean latitudes. The periods of the latitude drift of the spots are found for five stars; the magnitudes of the spot-latitude drift rates are lower than the corresponding value for sunspots by a factor of 1.5–3.     

Quality control in geochemistry from a comparison of four central tendency and five dispersion estimators and example of a geochemical reference material

Data quality control in geochemistry constitutes a fundamental problem that is still to be solved from the application of statistics and computation. We used refined Monte Carlo simulations of 10,000 replications and 190 independent experiments for sample sizes of 5 to 100. Statistical contaminations of 1 to 4 observations were used to compare 9 statistical parameters (4 central tendency—mean, median, trimean, and Gastwirth mean, and 5 dispersion estimates—standard deviation, median absolute deviation, S _n , Q _n , and \( {\widehat{\sigma}}_n \)). The presence of discordant observations in the data arrays rendered the outlier-based and robust parameters to disagree with each other. However, when the mean and standard deviation (outlier-based parameters) were estimated from censored data arrays obtained after the identification and separation of outlying observations, they generally provided a better estimate of the population than the robust estimates obtained from the original data arrays. This inference is contrary to the general belief, and therefore, reasons for the better performance of the outlier-based methods as compared to the robust methods are suggested. However, when all parameters were estimated from censored arrays and appropriate precise and accurate correction factors put forth in this work were applied, all of them became fully consistent, i.e., the mean agreed with the median, trimean and Gastwirth mean, and the standard deviation with the median absolute deviation, S _n , Q _n , and \( {\widehat{\sigma}}_n \). An example of inter-laboratory chemical data for a Hawaiian reference material BHVO-1 included sample sizes from 5 to 100, which showed that small samples of up to 20 provide inconsistent estimates, whereas larger samples of 20–100, especially >40, were more appropriate for estimating statistical parameters through robust or outlier-based methods. Although all statistical estimators provided consistent results, our simulation study shows that it is better to use the censored sample mean and population standard deviation as the best estimates.     

Physical simulation research on evolution laws of clay aquifuge stability during slice mining

The aquifuge stability is the key to study the impacts of coal mining on the aquifer. Based on the geological conditions of a mine in Yili of Xinjiang, China, this paper has studied the stability evolution laws of clay aquifuge during extremely thick coal seam mining by similar material simulation experiment in the laboratory. For the water-swelling and expansion property of clay aquifuge, the reasonable proportion of the similar material is firstly determined by taking the uniaxial compressive strength and the permeability coefficient as core indexes. Then, the overlying strata movement coupled solid–liquid physical model is established. In addition, the aquifuge deformation, the water level changes of the aquifer, and the height of fracture zone in overburden are analyzed during the slice mining. The research results indicate that the clay aquifuge will gradually occur instability failure during the mining of the working face, and the aquifuge stability has the threshold effects. When the ratio of the vertical displacement of the aquifuge to the thickness is D_v/T?≤?58.0%, the ratio of the horizontal displacement to the thickness is D_h/T?≤?17.0%, and the height of fractured zone in overburden is below the aquifuge, the mining-induced fractures may be closed and the aquifuge stability could be maintained. If D_v/T?≥?75.0%, D_h/T?≥?23.9%, and the height of fractured zone in overburden is within the aquifuge, the fractures will develop and connect the aquifuge and the stability failure of the aquifuge will occur, which has a direct correlation with the mining height.     

Estimating hydraulic parameters of a heterogeneous aquitard using long-term multi-extensometer and groundwater level data

The classical aquitard-drainage model COMPAC has been modified to simulate the compaction process of a heterogeneous aquitard consisting of multiple sub-units (Multi-COMPAC). By coupling Multi-COMPAC with the parameter estimation code PEST++, the vertical hydraulic conductivity (K _v) and elastic (S _ske) and inelastic (S _skp) skeletal specific-storage values of each sub-unit can be estimated using observed long-term multi-extensometer and groundwater level data. The approach was first tested through a synthetic case with known parameters. Results of the synthetic case revealed that it was possible to accurately estimate the three parameters for each sub-unit. Next, the methodology was applied to a field site located in Changzhou city, China. Based on the detailed stratigraphic information and extensometer data, the aquitard of interest was subdivided into three sub-units. Parameters K _v, S _ske and S _skp of each sub-unit were estimated simultaneously and then were compared with laboratory results and with bulk values and geologic data from previous studies, demonstrating the reliability of parameter estimates. Estimated S _skp values ranged within the magnitude of 10^?4 m^?1, while K _v ranged over 10^?10–10^?8 m/s, suggesting moderately high heterogeneity of the aquitard. However, the elastic deformation of the third sub-unit, consisting of soft plastic silty clay, is masked by delayed drainage, and the inverse procedure leads to large uncertainty in the S _ske estimate for this sub-unit.