Investigating spatial variability of vertical water fluxes through the streambed in distinctive stream morphologies using temperature and head data

Investigating the interaction of groundwater and surface water is key to understanding the hyporheic processes. The vertical water fluxes through a streambed were determined using Darcian flux calculations and vertical sediment temperature profiles to assess the pattern and magnitude of groundwater/surface-water interaction in Beiluo River, China. Field measurements were taken in January 2015 at three different stream morphologies including a meander bend, an anabranching channel and a straight stream channel. Despite the differences of flux direction and magnitude, flux directions based on vertical temperature profiles are in good agreement with results from Darcian flux calculations at the anabranching channel, and the Kruskal-Wallis tests show no significant differences between the estimated upward fluxes based on the two methods at each site. Also, the upward fluxes based on the two methods show similar spatial distributions on the streambed, indicating (1) that higher water fluxes at the meander bend occur from the center of the channel towards the erosional bank, (2) that water fluxes at the anabranching channel are higher near the erosional bank and in the center of the channel, and (3) that in the straight channel, higher water fluxes appear from the center of the channel towards the depositional bank. It is noted that higher fluxes generally occur at certain locations with higher streambed vertical hydraulic conductivity (K _v) or where a higher vertical hydraulic gradient is observed. Moreover, differences of grain size, induced by stream morphology and contrasting erosional and depositional conditions, have significant effects on streambed K _v and water fluxes.     

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).     

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.     

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.     

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.     

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.     

Masses of the visual components and black holes in X-ray novae: Effects of proximity of the components

It is shown that the approximation of the complex, tidally distorted shape of a star as a circular disc with local line profiles and a linear limb-darkening law, which is usually applied when deriving equatorial stellar rotation velocities from line profiles, leads to overestimation of the equatorial velocity V _rot sin i and underestimation of the component mass ratio q = M _x /M _v . A formula enabling correction of the effect of these simplifying assumptions on the shape of a star is used to re-determine the mass ratios q and the masses of the black holes M _x and visual components M _v in low-mass X-ray binary systems containing black holes. Taking into account the tidal–rotational distortion of the stellar shape can significantly increase the mass ratios q = M _x /M _v , reducing M _v , while M _x changes only slightly. The resulting distribution of M _v attains its maximum near M _v ? 0.35M _⊙, in disagreement with the results of population synthesis computations realizing standard models for Galactic X-ray novae with black holes. Possible ways to overcome this inconsistency are discussed. The derived distribution of M _x also differs strongly from the mass distribution for massive stars in the Galaxy.     

Measurement of streambed hydraulic conductivity and its anisotropy

 A method is described for the measurement of streambed hydraulic conductivity. Unlike permeameter methods, this method applies straight and l-shaped standpipes directly to streambeds for measurements of in-situ hydraulic conductivity in the vertical (K _v ) and horizontal (K _h ) directions, as well as in other oblique directions (K _s ). This method has advantages in determination of K _v values over grain-size analysis, permeameter tests, or slug test methods. Also unique to this method is that it provides K _s values of a streambed. The measured results can be used to construct a hydraulic conductivity ellipse and to evaluate the anisotropy of streambed sediments. Field examples from the Republican River, Nebraska, demonstrated the usefulness of this method in the determination of streambed hydraulic conductivity and anisotropy along or across a river channel. Results indicate that the K _h is about three to four times larger than K _v , whereas K _s values are larger than K _v but smaller than K _h . Received: 6 March 2000 · Accepted: 18 April 2000     

Distribution rule of the in situ stress state and its influence on the permeability of a coal reservoir in the southern Qinshui Basin,China

The recent development of the coalbed methane (CBM) industry has a significant role in advancing hydraulic fracturing theory and technology. However, further development requires a better understanding of how fractures influence reservoir permeability. In situ stress data from 54 CBM wells in the southern Qinshui Basin, China, were obtained by the injection/falloff test method to analyse the effect of in situ stress on the permeability of the CBM reservoir. The types of in situ stress states were classified, and the coal reservoir permeability under different in situ stress states was analysed. The results indicate that the maximum horizontal principal stress (σ_H), minimum horizontal principal stress (σ_h) and vertical principal stress (σ_v) all have positive linear relationships with the coal seam burial depth. Three in situ stress states were observed from the shallow to deep regions of the CBM reservoir in the study area: σ_H?>?σ_h?>?σ_v, σ_H?>?σ_v?>?σ_h and σ_v?>?σ_H?>?σ_h, which account for 9, 76 and 15% of the test wells, respectively. Coal reservoir permeability decreases with increasing horizontal principal stress, whereas it first decreases with increasing σ_v, then increases and finally decreases. The variation in permeability with σ_v is due to the conversion of the in situ stress states. Coal reservoir permeability has obvious differences under different in situ stress states. The permeability is the largest when σ_v?>?σ_H?>?σ_h, followed by σ_H?>?σ_h?>?σ_v and smallest when σ_H?>?σ_v?>?σ_h. The permeability differences are caused by the fracture propagation shape of the rock strata under different in situ stress states.     

Improved barometric and loading efficiency estimates using packers in monitoring wells

Measurement of barometric efficiency (BE) from open monitoring wells or loading efficiency (LE) from formation pore pressures provides valuable information about the hydraulic properties and confinement of a formation. Drained compressibility (α) can be calculated from LE (or BE) in confined and semi-confined formations and used to calculate specific storage (S _s). S _s and α are important for predicting the effects of groundwater extraction and therefore for sustainable extraction management. However, in low hydraulic conductivity (K) formations or large diameter monitoring wells, time lags caused by well storage may be so long that BE cannot be properly assessed in open monitoring wells in confined or unconfined settings. This study demonstrates the use of packers to reduce monitoring-well time lags and enable reliable assessments of LE. In one example from a confined, high-K formation, estimates of BE in the open monitoring well were in good agreement with shut-in LE estimates. In a second example, from a low-K confining clay layer, BE could not be adequately assessed in the open monitoring well due to time lag. Sealing the monitoring well with a packer reduced the time lag sufficiently that a reliable assessment of LE could be made from a 24-day monitoring period. The shut-in response confirmed confined conditions at the well screen and provided confidence in the assessment of hydraulic parameters. A short (time-lag-dependent) period of high-frequency shut-in monitoring can therefore enhance understanding of hydrogeological systems and potentially provide hydraulic parameters to improve conceptual/numerical groundwater models.     

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.     

Predicting saturated hydraulic conductivity using particle swarm optimization and genetic algorithm

Soil saturated hydraulic conductivity (K_s) is considered as soil basic hydraulic property, and its precision estimation is a key element in modeling water flow and solute transport processes both in the saturated and vadose zones. Although some predictive methods (e.g., pedotransfer functions, PTFs) have been proposed to indirectly predict K_s, the accuracy of these methods still needs to be improved. In this study, some easily available soil properties (e.g., particle size distribution, organic carbon, calcium carbonate content, electrical conductivity, and soil bulk density) are employed as input variables to predict K_s using a fuzzy inference system (FIS) trained by two different optimization techniques: particle swarm optimization (PSO) and genetic algorithm (GA). To verify the derived FIS, 113 soil samples were taken, and their required physical properties were measured (113 sample points?×?7 factors?=?791 input data). The initial FIS is compared with two methods: FIS trained by PSO (PSO-FIS) and FIS trained by GA (GA-FIS). Based on experimental results, all three methods are compared according to some evaluation criteria including correlation coefficient (r), modeling efficiency (EF), coefficient of determination (CD), root mean square error (RMSE), and maximum error (ME) statistics. The results showed that the PSO-FIS model achieved a higher level of modeling efficiency and coefficient of determination (R²) in comparison with the initial FIS and the GA-FIS model. EF and R² values obtained by the developed PSO-FIS model were 0.69 and 0.72, whereas they were 0.63 and 0.54 for the GA-FIS model. Moreover, the results of ME and RMSE indices showed that the PSO-FIS model can estimate soil saturated hydraulic conductivity more accurate than the GA-FIS model with ME?=?10.4 versus 11.5 and RMSE?=?5.2 versus 5.5 for PSO-FIS and GA-FIS, respectively.     

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.     

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.     

Development of a technique for IR spectroscopic determination of nitrogen content and aggregation degree in diamond crystals

A technique for IR spectroscopic determination of the total nitrogen content N _S in the form of A-and B ₁-defects is suggested. It provides for the computer processing and decomposition of IR spectra into constituent bands, calculation of the total absorption band area S _N and individual areas S _A and S _B1 and their normalization with respect to the total area of the diamond intrinsic absorption S ₀, with the normalization coefficients K _S , K _A , and K _B1 being calculated. Based on the analysis of the IR spectra of 60 octahedral diamond crystals from the Mir and Yubileinaya pipes (Sakha-Yakutiya), the empirical functions N _S = 911.85 K _S ^0.9919 ppm (R ² = 0.9859), N _A = 1185.6 K _A ^1.1511 ppm (R ² = 0.8703), and N _B1 = 911.85 K _S ^0.9919 ? 1185.6 K _A ^1.1511 ppm have been defined.     

Effects of salt solutions on the hydro-mechanical behavior of compacted GMZ01 Bentonite

In this study, the effects of salinity of infiltrating solutions on the swelling strain, compressibility, and hydraulic conductivity of compacted GMZ01 Bentonite were investigated. After swelling under vertical load using either distilled water or NaCl solutions with concentrations of 0.1, 0.5 M, and 1 M, laboratory oedometer tests were conducted on the compacted GMZ01 Bentonite. Based on the oedometer test results, hydraulic conductivity was determined using the Casagrande’s method. Results show that the swelling strain of highly compacted GMZ01 Bentonite decreases as the concentration of NaCl solution increases. The compression index C _c ^* increases and then turns to decrease with an increase in the vertical stress or a decrease in the void ratio for different solutions, and the C _c ^* decreases as the concentration of NaCl solution increases. The secondary consolidation coefficient C _α increases linearly with the increase of the compression index C _c ^* . Furthermore, a bi-linear relationship between the swelling index C _s ^* and the secondary consolidation coefficient C _α can be characterized clearly. The hydraulic conductivity increases as the concentration of NaCl solution increases, however, this increase can be prevented if a high confining stress is applied.     

The mass of the compact object in the X-ray binary her X-1/HZ her

We have obtained the first estimates of the masses of the components of the Her X-1/HZ Her X-ray binary system taking into account non-LTE effects in the formation of the H _γ absorption line: m _x = 1.8 M _⊙ and m _v = 2.5 M _⊙. These mass estimates were made in a Roche model based on the observed radial-velocity curve of the optical star, HZ Her. The masses for the X-ray pulsar and optical star obtained for an LTE model lie are m _x = 0.85 ± 0.15 M _⊙ and m _v = 1.87 ± 0.13 M _⊙. These mass estimates for the components of Her X-1/HZ Her derived from the radial-velocity curve should be considered tentative. Further mass estimates from high-precision observations of the orbital variability of the absorption profiles in a non-LTE model for the atmosphere of the optical component should be made.     

Statistical and geostatistical features of streambed hydraulic conductivities in the Platte River, Nebraska

This paper presents streambed hydraulic conductivities of the Platte River from south-central to eastern Nebraska. The hydraulic conductivities were determined from river channels using permeameter tests. The vertical hydraulic conductivities (K _v ) from seven test sites along this river in south-central Nebraska belong to one statistical population. Its mean value is 40.2 m/d. However, the vertical hydraulic conductivities along four transects of the Ashland test site in eastern Nebraska have lower mean values, are statistically different from the K _v values in south-central Nebraska, and belong to two different populations with mean values of 20.7 and 9.1 m/d, respectively. Finer sediments carried from the Loup River and Elkhorn River watersheds to the eastern reach of the Platte River lowers the vertical hydraulic conductivity of the streambed. Correlation coefficients between water depth and K _v values along a test transect indicates a positive correlation – a larger K _v usually occurs in the part of channel with deeper water. Experimental variograms derived from the vertical hydraulic conductivities for several transects across the channels of the Platte River show periodicity of spatial correlation, which likely result from periodic variation of water depth across the channels. The sandy to gravelly streambed contains very local silt and clay layers; spatially continuous low-permeability streambed was not observed in the river channels. The horizontal hydraulic conductivities were larger than the vertical hydraulic conductivities for the same test locations.     

Determination of the anisotropy of an upper streambed layer in east-central Nebraska, USA

Information on the anisotropy of streambed hydraulic conductivity (K) is a necessity for analyses of water exchange and solute transport in the hyporheic zone. An approach is proposed for the determination of K, developed from existing in-situ permeameter test methods. The approach is based on determination of vertical and horizontal hydraulic conductivity of streambed sediments on-site and eliminates the effects of vertical flow in the hyporheic zone and stream-stage fluctuation, which normally influence in situ permeameter tests. The approach was applied to seven study sites on four tributaries of the Platte River in east-central Nebraska, USA. On-site permeameter tests conducted on about 172 streambed cores for the determination of vertical hydraulic conductivity (K _v) and horizontal hydraulic conductivity (K _h) at the study sites indicate that the study sites have relatively small anisotropic ratios, ranging from 0.74 to 2.40. The ratios of K _h to K _v from individual locations within a study site show greater variation than the anisotropic ratios from the mean or median K at each of the study sites.     

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.