Constraints on hydraulic parameters and implications for groundwater flux across the upland–estuary interface

Analyses of independent laboratory- and field-scale measurements from two sites on Sapelo Island, Georgia reveal heterogeneity in hydraulic parameters across the upland–estuary interface. Regardless of the method used (short-duration pumping tests, amplitude attenuation of tidal pumping data, sediment grain size distributions, and falling head permeameter tests), we obtain hydraulic conductivity of 10⁻⁴ m s⁻¹ for the fine-grained, well-sorted, clean sands that make up the upland areas. Proximal to the upland–estuary boundary, the tidal pumping analyses and permeameter tests suggest that hydraulic conductivities decrease by more than two orders of magnitude, a result consistent with the presence of a clogging layer. Such a clogging layer may arise due to a variety of physical, chemical, or biological processes. The extent and orientation of the layers of reduced hydraulic conductivity near the upland–estuary boundary influence the nature of the aquifer's response to tidal forcing. Where the lower conductivity layer forms a relatively flat creek bank, tidal pumping produces a primarily mechanical response in the adjacent aquifer. Where the creek bank is nearly vertical, there is a more direct hydraulic connection between the tidal creek and the adjacent aquifer. The clogging layer likely contributes to the development of complicated flow pathways across the upland–estuary boundary. Effective flow paths calculated from tidal pumping data terminate within the marsh, beyond the boundary of the upland aquifer, suggesting a diffuse regime of groundwater discharge in the marsh. We postulate that, in many settings, submarsh flow may be as important as seepage faces for groundwater discharge into the marsh–estuary complex.     

Clogging process by suspended solids during groundwater artificial recharge: Evidence from lab simulations and numerical modelling

Permeability reduction of infiltration media due to suspended solid (SS) clogging is the bane of groundwater artificial recharge. To overcome the clogging problem and advance the understanding of the process‐based spatial‐temporal evolution of SS clogging, a 1D laboratory column simulation was carried out, followed by numerical modelling of the experimental data in this study. It was found that clogging caused a reduction in the hydraulic conductivity (K) in the upper layer at the beginning and extended deeper to approximately 50 cm, and no reduction in K was detected below 52 cm throughout the experimental period of 129 hr. The most clogged layer spanned from the surface to a depth of 11 cm, and the middle 11–52 cm was characterized by a slight decrease in K. The clogging rates of the different layers decreased with the depth, which was based on data analysis, with the largest value of 0.038 hr^?1 in the upper 1 cm. The overall K began to decrease from the surface layer and was increasingly affected by clogging with time. A mathematical model was established to simulate the SS clogging process evolution based on considerations of the attachments and detachments of particles. Then the model was applied to perform several scenario analyses after calibration and validation using the data obtained in the experiment. The simulation results indicated that the SS concentration was much more sensitive than the groundwater depth below the land surface, and 10 days of constant recharge is recommended as the disposal cycle of the clogged layer under the given conditions.     

Comparison between gradient based UCODE_2005 and the ensemble Kalman Filter for transient groundwater flow inverse modeling

Gradient based UCODE_2005 and data assimilation based on the Ensemble Kalman Filter(EnKF) are two different inverse methods. A synthetic two-dimensional flow case with four no-flow boundaries is used to compare the UCODE_2005 with the Ensemble Kalman Filter(EnKF) for their efficiency to inversely calculate and calibrate a hydraulic conductivity field based on hydraulic head data. A zonal, random heterogeneous conductivity field is calibrated by assimilating the time series of heads observed in monitoring wells. The study results indicate that the two inverse methods, UCODE_2005 and EnKF, could be used to calibrate the hydraulic conductivity field to a certain degree. More available observations and information about the conductivity field, more accurate inverse results will be obtained for the UCODE_2005. On the other hand, for a realistic zonal heterogeneous hydraulic conductivity field, EnKF can only efficiently determine the hydraulic conductivity field at the first several assimilated time steps. The results obtained by the UCODE_2005 look better than those by the EnKF. This is possibly due to the fact that the UCODE_2005 uses observed head data at every time step, while EnKF can only use observed heads at first several steps due to the filter divergence problem.     

Bank filtration in the sandy littoral zone of Lake Tegel (Berlin): Structure and dynamics of the biological active filter zone and clogging processes

This study focused on clogging processes and on the benthic microalgal and meiofaunal assemblage in the sandy littoral zone of Lake Tegel, which are significantly involved in bank filtration, in a long-term. Our approach combined field studies and “in situ” experiments to highlight the structure of the biological active filter zone as well as the mechanisms and effects of clogging in the interstices that influence the infiltration process.Campaigns to measure “in situ” infiltration rates and hydraulic potential were conducted monthly from March 2004 to April 2005. Meiofaunal abundances and fine particulate organic matter (FPOM) were determined every 6 weeks in freeze cores down to depths of 50 cm. In parallel, concentrations of carbon, nitrogen and chlorophyll a were measured in samples of unfrozen sediment cores, that were divided in 1-cm steps down to depths of ≥10 cm. Similar sediment profiles were generated for analysis of colloidal carbohydrates, extracellular polymeric substances (EPS) and proteins between December 2005 and June 2006. Electron microscopy was used to visualize biofilm structure. Long-term experiments with natural FPOM and melamine resin particles as fluorescent tracers were performed to study “in situ” particle retention and transport, respectively. Additionally seston input was quantified during a 1-week period in April 2005.Infiltration rates showed a high temporal and spatial variability, but were not correlated with hydraulic conductivities as hydraulic gradients changed a lot. Likewise a correlation between infiltration rates and hydraulic potentials was not observed, indicating clogging processes. These are triggered to a high extend by biological compounds. In addition, seston input and intermittent gas intrusion are considered to reduce the hydraulic conductivity considerably. No significant “in situ” transport of inert natural fluorescent tracers was observed. However, a complete and permanent clogging of the sandy sediment does not occur, and daily infiltration rates of 0.7-27 L m⁻² h⁻¹ (mean 9 L m⁻² h⁻¹) guarantee a sufficient water supply by bank filtration for decades.     

Mechanism of Suspended Kaolinite Particle Clogging in Porous Media During Managed Aquifer Recharge

Managed aquifer recharge is an effective strategy for urban stormwater management. Chemical ions are normally retained in stormwater and groundwater and may accelerate clogging during the recharge process. However, the effect of water chemistry on physical clogging has not previously been investigated. In this study, we investigated the hydrogeochemical mechanism of saturated porous media clogging in a series of column experiments. The column was packed with river sand and added suspensions of kaolinite particles. Calcium chloride and sodium chloride are used as representative ions to study chemical effects. We found that an increase in ionic strength resulted in retention of kaolinite solids in the column, with a breakthrough peak of C/C₀ value of 1 to 0.2. The corresponding hydraulic conductivity decreased with increased solids clogging. Divalent cations were also found to have a greater influence on kaolinite particle clogging than monovalent cations. The enhanced hydrochemical-related clogging was caused by kaolinite solids flocculating and increasing the deposition rate coefficient by 1 to 2 times in high ionic strength conditions. Three clogging mechanisms of kaolinite solids are proposed: surface filtration, inner blocking, and attachment. This study further deepens the understanding of the mechanisms of solids clogging during aquifer recharge and demonstrates the significance of ionic strength on recharge clogging risk assessments.     

Stochastic analysis of groundwater flow in aquifers with leakage

 Stochastic analysis of one- and two-dimensional flow through a shallow semi-confined aquifer with spatially variable hydraulic conductivity K represented by a stationary (statistically homogeneous) random process is carried out by using the spectral technique. The hydraulic head covariance functions for flows in a semi-confined aquifer bounded by a leaky layer above and an impervious stratum below are derived by assuming that the randomness forcing the head variation to originate from the hydraulic conductivity field of the aquifer. The head covariance functions are studied using two convenient forms of the logarithmic hydraulic conductivity process. The results demonstrate the significant reduction in the head variances and covariances due to the presence of a leaky layer. The hydraulic head correlation distance is also reduced greatly due to the presence of the leaky layer.     

Numerical Modeling of Unsaturated Flow in Wastewater Soil Absorption Systems

It is common practice in the United States to use wastewater soil absorption systems (WSAS) to treat domestic wastewater. WSAS are expected to provide efficient, long-term removal of wastewater contaminants prior to ground water recharge. Soil clogging at the infiltrative surface of WSAS occurs due to the accumulation of suspended solids, organic matter, and chemical precipitates during continued wastewater infiltration. This clogging zone (CZ) creates an impedance to flow, restricting the hydraulic conductivity and rate of infiltration. A certain degree of clogging may improve the treatment of wastewater by enhancing purification processes, in part because unsaturated flow is induced and residence times are significantly increased. However, if clogging becomes excessive, the wastewater pond height at the infiltrative surface can rise to a level where system failure occurs. The numerical model HYDRUS-2D is used to simulate unsaturated flow within WSAS to better understand the effect of CZs on unsaturated flow behavior and hydraulic retention times in sandy and silty soil. The simulations indicate that sand-based WSAS with mature CZs are characterized by a more widely distributed flow regime and longer hydraulic retention times. The impact of clogging on water flow within the silt is not as substantial. For sand, increasing the hydraulic resistance of the CZ by a factor of three to four requires an increase in the pond height by as much as a factor of five to achieve the same wastewater loading. Because the degree of CZ resistance directly influences the pond height within a system, understanding the influence of the CZ on flow regimes in WSAS is critical in optimizing system design to achieve the desired pollutant-treatment efficiency and to prolong system life.     

In situ Mixing of Organic Matter Decreases Hydraulic Conductivity of Denitrification Walls in Sand Aquifers

In a previous study, a denitrification wall was constructed in a sand aquifer using sawdust as the carbon substrate. Ground water bypassed around this sawdust wall due to reduced hydraulic conductivity. We investigated potential reasons for this by testing two new walls and conducting laboratory studies. The first wall was constructed by mixing aquifer material in situ without substrate addition to investigate the effects of the construction technique (mixed wall). A second, biochip wall, was constructed using coarse wood chips to determine the effect of size of the particles in the amendment on hydraulic conductivity. The aquifer hydraulic conductivity was 35.4 m/d, while in the mixed wall it was 2.8 m/d and in the biochip wall 3.4 m/d. This indicated that the mixing of the aquifer sands below the water table allowed the particles to re-sort themselves into a matrix with a significantly lower hydraulic conductivity than the process that originally formed the aquifer. The addition of a coarser substrate in the biochip wall significantly increased total porosity and decreased bulk density, but hydraulic conductivity remained low compared to the aquifer. Laboratory cores of aquifer sand mixed under dry and wet conditions mimicked the reduction in hydraulic conductivity observed in the field within the mixed wall. The addition of sawdust to the laboratory cores resulted in a significantly higher hydraulic conductivity when mixed dry compared to cores mixed wet. This reduction in the hydraulic conductivity of the sand/sawdust cores mixed under saturated conditions repeated what occurred in the field in the original sawdust wall. This indicated that laboratory investigations can be a useful tool to highlight potential reductions in field hydraulic conductivities that may occur when differing materials are mixed under field conditions.     

Using a pore‐scale model to quantify the effect of particle re‐arrangement on pore structure and hydraulic properties

A pore‐scale model based on measured particle size distributions has been used to quantify the changes in pore space geometry of packed soil columns resulting from a dilution in electrolyte concentration from 500 to 1 mmol l^?1 NaCl during leaching. This was applied to examine the effects of particle release and re‐deposition on pore structure and hydraulic properties. Two different soils, an agricultural soil and a mining residue, were investigated with respect to the change in hydraulic properties. The mining residue was much more affected by this process with the water saturated hydraulic conductivity decreasing to 0·4% of the initial value and the air‐entry value changing from 20 to 50 cm. For agricultural soil, there was little detectable shift in the water retention curve but the saturated hydraulic conductivity decreased to 8·5% of the initial value. This was attributed to localized pore clogging (similar to a surface seal) affecting hydraulic conductivity, but not the microscopically measured pore‐size distribution or water retention. We modelled the soil structure at the pore scale to explain the different responses of the two soils to the experimental conditions. The size of the pores was determined as a function of deposited clay particles. The modal pore size of the agricultural soil as indicated by the constant water retention curve was 45 µm and was not affected by the leaching process. In the case of the mining residue, the mode changed from 75 to 45 µm. This reduction of pore size corresponds to an increase of capillary forces that is related to the measured shift of the water retention curve. Copyright © 2007 John Wiley & Sons, Ltd.     

Evaluating the effect of using artificial pore water on the quality of laboratory hydraulic conductivity measurements of peat

Pore dilation, the compaction of humic acids on peat fibres due to the process of flocculation, causes the hydraulic conductivity of peat to increase with increasing pore water electrical conductivity. This is a reversible process and a reduction in the pore water conductivity produces a decrease in the hydraulic conductivity due to the constriction of pores. We verify how this dilation and constriction of pores, resulting from the application of artificial pore water (primarily deionized water), affects laboratory measurements of the hydraulic conductivity of peat. Repeat measurements of the hydraulic conductivity were performed on samples of Sphagnum peat. It is shown that the application of deionized water during constant head permeameter tests causes a significant decrease in the hydraulic conductivity. Between tests, the hydraulic conductivity of the peat continues to decline without an associate decrease in the pore water electrical conductivity because of a lagged pore constriction effect. We suggest that the use of artificially high or low pore water electrical conductivities, during laboratory hydraulic conductivity measurements, is likely to lead to significant errors. Experimental protocols must, therefore, be revised to take better account of the pore water chemistry. The ionic concentrations of the natural pore fluid should be replicated during hydraulic conductivity tests, either by using pore fluid extracted from the study site or by artificially replicating the major ionic composition of the natural pore fluid. In addition, prior to the hydraulic conductivity measurements, peat samples should be flushed with this solution until the hydraulic conductivity stabilizes and the samples subsequently allowed to equilibrate. Copyright © 2010 John Wiley & Sons, Ltd.     

Temporal and spatial variation of infiltration in urban green infrastructure

Infiltration is the primary mechanism in green stormwater infrastructure (GSI) systems to reduce the runoff volume from urbanized areas. Soil hydraulic conductivity is most important in influencing GSI infiltration rates. Saturated hydraulic conductivity (Ksat) is a critical parameter for GSI design and post-construction performance. However, Ksat measurement in the field is problematic due to temporal and spatial variability and measurement errors. This review paper focuses on a comparison of methods for in-situ Ksat measurement and the causes of temporal and spatial variations of Ksat within GSI systems. Automated infiltration testing methods, such as the Modified Philip–Dunne (MPD) and SATURO infiltrometers, show promise for efficient Ksat measurements. Soil Ksat values can change over time and substantially vary throughout a GSI, which can be attributed to multiple factors, including but not limited to temperature changes, soil composition and properties, soil compaction level, plant root morphology and distribution, biological and macrofauna activities in the soil, inflow sediment characteristics, quality of infiltrating water, and measurement errors. There is evidence that infiltration rates in vegetated urban GSI systems are sustained given an appropriate GSI design, reasonable concentration of suspended sediments in the inflow runoff, and routine maintenance procedures. These observations indicate that clogging can be counteracted by processes that tend to increase the soil hydraulic conductivity (e.g., plant root and biological activities). This self-sustainability underlines that infiltration-based GSI systems are a reliable long-term stormwater management solution. Recommendations on how to incorporate the temporal changes of Ksat in GSI design and on obtaining a spatially-representative Ksat for the GSI design are presented.     

Numerical method of moments for solute transport in a nonstationary flow field conditioned on hydraulic conductivity and head measurements

The unconditional stochastic studies on groundwater flow and solute transport in a nonstationary conductivity field show that the standard deviations of the hydraulic head and solute flux are very large in comparison with their mean values (Zhang et al. in Water Resour Res 36:2107–2120, 2000; Wu et al. in J Hydrol 275:208–228, 2003; Hu et al. in Adv Water Resour 26:513–531, 2003). In this study, we develop a numerical method of moments conditioning on measurements of hydraulic conductivity and head to reduce the variances of the head and the solute flux. A Lagrangian perturbation method is applied to develop the framework for solute transport in a nonstationary flow field. Since analytically derived moments equations are too complicated to solve analytically, a numerical finite difference method is implemented to obtain the solutions. Instead of using an unconditional conductivity field as an input to calculate groundwater velocity, we combine a geostatistical method and a method of moment for flow to conditionally simulate the distributions of head and velocity based on the measurements of hydraulic conductivity and head at some points. The developed theory is applied in several case studies to investigate the influences of the measurements of hydraulic conductivity and/or the hydraulic head on the variances of the predictive head and the solute flux in nonstationary flow fields. The study results show that the conditional calculation will significantly reduce the head variance. Since the hydraulic head measurement points are treated as the interior boundary (Dirichlet boundary) conditions, conditioning on both the hydraulic conductivity and the head measurements is much better than conditioning only on conductivity measurements for reduction of head variance. However, for solute flux, variance reduction by the conditional study is not so significant.     

Steady state flow with hydraulic conductivity change around large diameter wells

Darcian flow law in aquifers assumes that the aquifer hydraulic conductivity is constant and the groundwater movement is due only to the piezometric level changes through hydraulic gradient. In practice, after the well development the aquifer just around the well has comparatively larger hydraulic conductivity and gradient. Patchy aquifer solutions in the literature consider sudden hydraulic conductivity changes with distance for the steady state flow. The change of transmissivity is demonstrated by the application of slope‐matching procedure to actual field data. It is the main purpose of this paper to derive simple analytical expressions for aquifer parameter evaluations with steadily decreasing hydraulic conductivity around the well. Spatial nonlinear hydraulic conductivity changes around a large‐diameter well within the depression cone of a confined aquifer are considered as exponentially decreasing functions of the radial distance. Copyright © 2010 John Wiley & Sons, Ltd.     

On the relation between steep monoclinal flexure zones and steep hydraulic gradients

Steep hydraulic gradients are found in association with steep monoclinal flexures. However, the physics of the reduction of the hydraulic conductivity, which is responsible for the steep gradients, has seldom been studied. We present results of hydrological and mechanical modeling aiming to study the effect of such steep hydraulic gradients demonstrated in the Judea Group Aquifer system, Israel. The hydrological configuration of steep dips and anisotropy between flows parallel and perpendicular to the bedding planes was simulated using the FEFLOW code. It exhibited a situation whereby part of the flow is oblique to the bedding planes and therefore some steepening of the hydraulic gradients occurred due to actual conductivity reduction. However, this reduction is not enough to account for the steeper gradients observed. The effect of a deep-seated reverse fault under the monocline on the permeability distribution within the structure was examined by numerical mechanical simulations. It exhibited a compressional stress distribution in the steep part of the monocline, which, due to shortening and closure of joints and voids, is presumably responsible for a significant pressure-induced permeability reduction. This process by itself in a layered structure, including interlayering of thin marl layers, could be responsible for the steep hydraulic gradients in the steep part of the monocline.     

Variably saturated groundwater modelling for optimizing managed aquifer recharge using trench infiltration

Spreading‐basin methods have resulted in more than 130 million cubic metres of recharge to the unconfined Navajo Sandstone of southern Utah in the past decade, but infiltration rates have slowed in recent years because of reduced hydraulic gradients and (or) clogging. Trench infiltration is a promising alternative technique for increasing recharge and minimizing evaporation. This paper uses a variably saturated flow model to further investigate the relative importance of the following variables on rates of trench infiltration to unconfined aquifers: saturated hydraulic conductivity, trench spacing and dimensions, initial water‐table depth, alternate wet/dry periods, and number of parallel trenches. Modelling results showed (1) increased infiltration with higher hydraulic conductivity, deeper initial water tables, and larger spacing between parallel trenches, (2) deeper or wider trenches do not substantially increase infiltration, (3) alternating wet/dry periods result in less overall infiltration than keeping the trenches continuously full, and (4) larger numbers of parallel trenches within a fixed area increases infiltration but with a diminishing effect as trench spacing becomes tighter. An empirical equation for estimating expected trench infiltration rates as a function of hydraulic conductivity and initial water‐table depth was derived and can be used for evaluating feasibility of trench infiltration in other hydrogeologic settings. Copyright © 2014 John Wiley & Sons, Ltd.     

Water balance prediction in stormwater infiltration basins using 2-D modeling: An application to evaluate the clogging process

Sustainable urban drainage systems are built along roads and in urban areas to collect urban runoff and avoid flooding, and to filter water pollutants. Sediment collected by runoff is deposited in the stormwater basin and progressively reduces water infiltration efficiency, leading to the clogging of the basin. To help stormwater basin managers and stakeholders better understand and predict clogging rates in order to elaborate maintenance plans and schedules, water transport prediction models are necessary. However,because of the heterogeneous sediment hydrodynamic properties inside the stormwater basin, a twodimensional(2-D) water flow model is required to predict water levels and possible overflow as accurately as possible. Saturated hydraulic conductivity(Ks) and sediment water retention curves were measured in the overall sediment layer of the stormwater basin, in addition to sediment layer thickness and organic matter content(11 sampling points). Sediment depth was used to predict organic matter(OM) content, and the OM was used to predict Ks. Water height in the basin was modeled with the HYDRUS-2 D model by taking into account the sediment hydrodynamic properties distribution. The HYDRUS-2 D model gave a satisfactory representation of the measured data. Scenarios of the hydraulic properties of stormwater basin sediment were tested over time, and hydraulic resistance, R, was calculated to assess the stormwater basin performance. Presently, after 20 years of functioning, the stormwater basin still ensures efficient water infiltration, but the first outflow(Hydraulic resistance,R 24 h)) is expected to appear in the next 5 years, and clogging(R 47 h) in the next 13 years. This 2-D water balance model makes it possible to integrate the hydrodynamic heterogeneity of a stormwater basin. It gives interesting perspectives to better predict 2-D/3-D contaminant transport.     

Pressure variations in peat as a result of gas bubble dynamics

Transient high pore‐water pressures, up to 50 cm higher than ambient pressure, developed over the summer season at various depths in a shallow (1 m) fen peat. The excess pressures had a pattern of gradual increases and sharp drops, and their initiation and release typically corresponded to abrupt changes in atmospheric pressure. We conclude that these phenomena depend on gas bubbles (probably methane) generated by biological activity, both by clogging pores and by building up pressure as they grow. These transient and spatially discontinuous high‐pressure zones were found using pressure transducers in sealed (backfilled) pits, but not in piezometers open to the atmosphere. Piezometers may provide a conduit for the release of gas and pressure, thus rendering them unsuitable for measuring this phenomenon. Although the development of localized zones of high pressure causes erratic and unpredictable hydraulic gradients, we suggest that their effect on the flow of water or solutes is offset by the reduced permeability caused by the bubbles, which allows them to be sustained. These zones, however, probably deflect flows driven by the dominant hydraulic gradients. Furthermore, they may cause the peat volume to adjust (swell). The use and interpretation of traditional methods for estimating hydraulic head and conductivity in peat soils thus require great caution. Copyright © 2004 John Wiley & Sons, Ltd.     

Electrokinetic Migration of Nitrate Through Heterogeneous Granular Porous Media

This study investigates and quantifies the influence of physical heterogeneity in granular porous media, represented by materials with different hydraulic conductivity, on the migration of nitrate, used as an amendment to enhance bioremediation, under an electric field. Laboratory experiments were conducted in a bench‐scale test cell under a low applied direct current using glass bead and clay mixes and synthetic groundwater to represent ideal conditions. The experiments included bromide tracer tests in homogeneous settings to deduce controls on electrokinetic transport of inorganic solutes in the different materials, and comparison of nitrate migration under homogeneous and heterogeneous scenarios. The results indicate that physical heterogeneity of subsurface materials, represented by a contrast between a higher‐hydraulic conductivity and lower‐hydraulic conductivity material normal to the direction of the applied electric field exerts the following controls on nitrate migration: (1) a spatial change in nitrate migration rate due to changes in effective ionic mobility and subsequent accumulation of nitrate at the interface between these materials; and (2) a spatial change in the voltage gradient distribution across the hydraulic conductivity contrast, due to the inverse relationship with effective ionic mobility. These factors will contribute to higher mass transport of nitrate through low hydraulic conductivity zones in heterogeneous porous media, relative to homogeneous host materials. Overall electrokinetic migration of amendments such as nitrate can be increased in heterogeneous granular porous media to enhance the in situ bioremediation of organic contaminants present in low hydraulic conductivity zones.     

Modifications to the Conduit Flow Process Mode 2 for MODFLOW-2005

As a result of rock dissolution processes, karst aquifers exhibit highly conductive features such as caves and conduits. Within these structures, groundwater flow can become turbulent and therefore be described by nonlinear gradient functions. Some numerical groundwater flow models explicitly account for pipe hydraulics by coupling the continuum model with a pipe network that represents the conduit system. In contrast, the Conduit Flow Process Mode 2 (CFPM2) for MODFLOW-2005 approximates turbulent flow by reducing the hydraulic conductivity within the existing linear head gradient of the MODFLOW continuum model. This approach reduces the practical as well as numerical efforts for simulating turbulence. The original formulation was for large pore aquifers where the onset of turbulence is at low Reynolds numbers (1 to 100) and not for conduits or pipes. In addition, the existing code requires multiple time steps for convergence due to iterative adjustment of the hydraulic conductivity. Modifications to the existing CFPM2 were made by implementing a generalized power function with a user-defined exponent. This allows for matching turbulence in porous media or pipes and eliminates the time steps required for iterative adjustment of hydraulic conductivity. The modified CFPM2 successfully replicated simple benchmark test problems.     

岱海水体及沉积物细菌多样性及群落组成特征

采用高通量测序技术,研究了内蒙古岱海流域入湖河流、湖水及沉积物细菌多样性及群落组成.结果显示,细菌多样性从高到低依次为：沉积物>河流>湖泊.聚类分析表明入湖河流、湖水和沉积物细菌群落可分为明显不同的3支,说明这3种生境中细菌群落结构有较大差异.物种注释结果表明,河流中优势细菌菌群为髌骨细菌(Patescibacteria)、变形菌门(Proteobacteria)、拟杆菌门(Bacteroidetes)和放线菌门(Actinobacteria);湖水中优势细菌菌群为放线菌门(Actinobacteria);而沉积物中优势细菌菌群为变形菌门和绿弯菌门(Chloroflexi).典范对应分析及Monte Carlo检验表明,电导率和悬浮物含量对水体中(河流与湖泊)细菌群落影响显著,二者共解释了细菌群落变化的86.5%;而冗余分析及Monte Carlo检验表明,泥深、磁化率和总有机碳对沉积物中细菌群落影响显著,三者共解释了细菌群落变化的47.9%.近30年来,岱海地区气候变化和人类活动导致湖水咸化,沉积物碳氮指标显著增长.岱海水体及沉积物细菌多样性及群落组成的差异及其主要驱动因...