Relationship between bed shear stress and suspended sediment concentration:annular flume experiments

In this study,annular flume experiments were carried out,using the sediment samples collected from the lower part of the inter-tidal zone at Xiaoyangkou,Jiangsu coast,China.The Ariathurai-Partheniades equation was used to determine the bed shear stress,by evaluating variations in the suspended sediment concentration within the water column.The derived relation between the bed shear stress and suspended sediment concentration shows that,at various stages of seabed erosion, suspended sediment concentration increases rapidly when the flow velocity is increased,but the pattern of change in the bed shear stress does not follow suit.At low concentrations,bed shear stress initially increases markedly with increasing flow velocity.However,when the concentration reaches an apparently critical level around 0.55 kg m"3,the rate of change in the bed shear stress abruptly slows down,or becomes almost constant,in response to further increases in the flow velocity.Results of experiments indicate that,from a critical level onward,suspended sediment concentration has a strong influence on the bed shear stress.     

COMBINED EXPERIMENTAL AND NUMERICAL INVESTIGATION OF FORMULAE FOR EROSION AND DEPOSITION OF COHESIVE SEDIMENTS

1 mTsonvCTIoxTIansport of fine cOhesive sedinds inevitably leads to the develoPment of sediment deposits inreservoirs or zones with small flow velocihes in lowland rivers. Such sediment dePosits often conshtUte aconsiderable ecological danger, as cOhesive sediInnts have the ability to bind POllutants like heavymetals, polycyclic aromatic hydrocrtons or chlorinated hydrocaIbons (F6f8iner, l989). In the case of anextreme hydraUlic event like a flood or a reservoir dePlehon, erosion of cone…     

Estuarine mud flocculation properties determined using an annular mini-flume and the LabSFLOC system

Most entrained estuarine sediment mass occurs as flocs. Parameterising flocculation has proven difficult as it is a dynamically active process dependent on a set of complex interactions between the sediment, fluid and the flow. However the natural variability in an estuary makes it difficult to study the factors that influence the behaviour of flocculation in a systematic manner. This paper presents preliminary results from a laboratory study that examined how floc properties of a natural estuarine mud from the Medway (UK), evolved in response to varying levels of suspended sediment concentration and induced turbulent shearing. The experiments utilised the LabSFLOC floc video camera system, in combination with an annular mini-flume to shear the suspended sediment slurries. The flows created in the mini-flume produced average shear stresses, at the floc sampling height, ranging from 0.01 N m⁻² to a peak of 1.03 N m⁻². Nominal suspended particulate matter concentrations of 100, 600 and 2000 mg l⁻¹ were introduced into the flume. The experimental runs produced individual flocs ranging in size from microflocs of 22.2 μm to macroflocs 583.7 μm in diameter. Average settling velocities ranged from 0.01 to 26.1 mm s⁻¹, whilst floc effective densities varied from 3.5 up to 2000 kg m⁻³. Low concentration and low shear stress were seen to produce an even distribution of floc mass between the macrofloc (>160 μm) and microfloc (<160 μm) fractions. As both concentration and stress rose, the proportion of macrofloc mass increased, until they represented over 80% of the suspended matter. A maximum average macrofloc settling velocity of 3.3 mm s⁻¹ was attained at a shear stress of 0.45 N m⁻². Peak turbulence conditions resulted in deflocculation, limiting the macrofloc fall velocity to only 1.1 mm s⁻¹ and placing over 60% of the mass in the microfloc size range. A statistical analysis of the data suggests that the combined influence of both suspended concentration and turbulent shear controls the settling velocity of the fragile, low density macroflocs.     

Modeling flocculation processes of fine-grained particles using a size-resolved method: Comparison with published laboratory experiments

The transport of fine-grained particles in estuarine and coastal waters is influenced by flocculation processes (aggregation and floc breakup). As a consequence, the particle size varies with time in the water column, and can be orders of magnitude larger than those of primary particles. In this study the variations in floc size is simulated using a size-resolved method, which approximates the real size distribution of particles by a range of size bins and solves a mass balance equation for each bin. To predict the size distribution both aggregation and breakup processes are included. The conventional rectilinear aggregation kernel is used which considers both turbulent shear and differential settling. The breakup kernel accounts for the fractal dimension of the flocs. A flocculation simulation is compared to the settling column lab experiments of Winterwerp [1998. A simple model for turbulence induced flocculation of cohesive sediment, Journal of Hydraulic Research, 36, 309–326], and a one-dimensional sediment transport model is verified with the observed variations in floc size and concentration over tidal cycles in a laboratory flume experiment of Bale et al. [2002. Direct observation of the formation and break-up of aggregates in an annular flume using laser reflectance particle sizing. In: Winterwerp, J.C., Kranenburg, C. (Eds.), Fine Sediment Dynamics in the Marine Environment. Elsevier, pp. 189–201]. The numerical simulations compare qualitatively and quantitatively well with the laboratory measurements, and the analysis of the two simulation results indicates that the median floc size can be correlated to the sediment concentration and Kolmogorov microscale. Sensitivity studies are conducted to explore the role of settling velocity and erosion rate. The results are not sensitive towards the formulation of settling velocity, but the parameterization of erosion flux is important. The studies show that for predicting the sediment deposition flux it is crucial to include flocculation processes.     

The influence of Pacific salmon decay products on near‐field streambed sediment and organic matter dynamics: a flume simulation

Pacific salmon are biogeomorphic agents shown to induce positive feedbacks on their natal watersheds. However, the literature documenting their ecological effects on in‐stream natal environments is more divisive. The disturbance salmon create during redd construction has the potential to reduce stream productivity. The pulse of salmon organic matter (SOM) and marine derived nutrients (MDNs) released during carcass decay has been reported as either stimulating in‐stream productivity or having no local effect. To evaluate the ecological costs and benefits of salmon spawning events, MDN delivery and storage processes need to be identified and quantified. A simulation was conducted in three flow‐through flumes (2 m × 2 m × 30 m) over a 33‐day period (consisting of 15 baseline, four MDN exposure, and 14 post‐exposure days) to assess near‐field sediment and organic matter dynamics during active and post‐spawn simulations. The objective of the study was to measure changes in the amounts and particle sizes of suspended and gravel‐stored fine sediment, in order to elucidate the process and significance of SOM recruitment to the gravel bed via sedimentation. Gravel beds in all flumes were enriched with SOM following treatments but the response was highest in the active spawn simulation. The more effective delivery in the active spawn simulation was attributed to its higher inorganic sediment concentration, which is known to enhance floc formation. Although the active spawn simulation delivered more SOM to the gravel bed, the post‐spawn phase may be equally important to natural streams because its decay phase is longer than the active spawn and consequently can provide SOM to the streambed as long as carcasses remain in‐stream. The delivery, and potential retention, of SOM to spawning streambeds and the intergravel environment may be particularly important for interior streams, which experience low flow conditions during the spawning phase and accordingly have the potential for hyporheic nutrient recruitment and storage. Copyright © 2014 John Wiley & Sons, Ltd.     

Sediment transport and resuspension due to combined motion of wave and current in the northern Adriatic Sea during a Bora event in January 2001: A numerical modelling study

The Adriatic Sea general circulation model coupled to a third generation wave model SWAN and a sediment transport model was implemented in the Adriatic Sea to study the dynamics of the sediment transport and resuspension in the northern Adriatic Sea (NAS) during the Bora event in January 2001. The bottom boundary layer (BBL) was resolved by the coupled model with high vertical resolution, and the mechanism of the wave–current interaction in the BBL was also represented in the model. The study found that, during the Bora event of 13–17 January 2001, large waves with significant wave height 2 m and period of 5 s were generated by strong winds in the northwestern shelf of the Adriatic where the direction of wave propagation was orthogonal to the current. The combined motion of the wave and current in the BBL increased the bottom stress over the western Adriatic shelf, resulting in stronger sediment resuspension there. Combining stronger bottom resuspension and strong upward vertical flux of resuspended sediments due to turbulent mixing, the model predicted that sediment concentration near the Po River was much higher than that predicted by the model run without wave forcing. The study also shows that wave–current interaction in the BBL reduced the western Adriatic Coastal Currents (WACCs) in the shallower north. It is concluded that wave forcing significantly changed the sediment distributions and increased the total horizontal fluxes over the western shelf. These results signified wave effect on sediment flux and distribution in the NAS, and suggested that waves cannot be neglected in the study of dynamics of sediment transport and resuspension in the shallow coastal seas. By including the tidal forcing in the coupled model, we also examined the effect of tides on the sediment transport dynamics in the NAS.     

CO2-induced dissolution of low permeability carbonates. Part II: Numerical modeling of experiments

We used the 3D continuum-scale reactive transport models to simulate eight core flood experiments for two different carbonate rocks. In these experiments the core samples were reacted with brines equilibrated with pCO₂ = 3, 2, 1, 0.5 MPa (Smith et al., 2013 [27]). The carbonate rocks were from specific Marly dolostone and Vuggy limestone flow units at the IEAGHG Weyburn-Midale CO₂ Monitoring and Storage Project in south-eastern Saskatchewan, Canada. Initial model porosity, permeability, mineral, and surface area distributions were constructed from micro tomography and microscopy characterization data. We constrained model reaction kinetics and porosity–permeability equations with the experimental data. The experimental data included time-dependent solution chemistry and differential pressure measured across the core, and the initial and final pore space and mineral distribution. Calibration of the model with the experimental data allowed investigation of effects of carbonate reactivity, flow velocity, effective permeability, and time on the development and consequences of stable and unstable dissolution fronts.The continuum scale model captured the evolution of distinct dissolution fronts that developed as a consequence of carbonate mineral dissolution and pore scale transport properties. The results show that initial heterogeneity and porosity contrast control the development of the dissolution fronts in these highly reactive systems. This finding is consistent with linear stability analysis and the known positive feedback between mineral dissolution and fluid flow in carbonate formations. Differences in the carbonate kinetic drivers resulting from the range of pCO₂ used in the experiments and the different proportions of more reactive calcite and less reactive dolomite contributed to the development of new pore space, but not to the type of dissolution fronts observed for the two different rock types. The development of the dissolution front was much more dependent on the physical heterogeneity of the carbonate rock. The observed stable dissolution fronts with small but visible dissolution fingers were a consequence of the clustering of a small percentage of larger pores in an otherwise homogeneous Marly dolostone. The observed wormholes in the heterogeneous Vuggy limestone initiated and developed in areas of greater porosity and permeability contrast, following pre-existing preferential flow paths.Model calibration of core flood experiments is one way to specifically constrain parameter input used for specific sites for larger scale simulations. Calibration of the governing rate equations and constants for Vuggy limestones showed that dissolution rate constants reasonably agree with published values. However the calcite dissolution rate constants fitted to the Marly dolostone experiments are much lower than those suggested by literature. The differences in fitted calcite rate constants between the two rock types reflect uncertainty associated with measured reactive surface area and appropriately scaling heterogeneous distribution of less abundant reactive minerals. Calibration of the power-law based porosity–permeability equations was sensitive to the overall heterogeneity of the cores. Stable dissolution fronts of the more homogeneous Marly dolostone could be fit with the exponent n = 3 consistent with the traditional Kozeny–Carman equation developed for porous sandstones. More impermeable and heterogeneous cores required larger n values (n = 6–8).     

Entropy model to assess sediment resuspension probability and trap efficiency of small dams

In the Brazilian drylands, there are tens of thousands of small dams. Despite their paramount importance to the rural population, they are rarely monitored. Water demand increases with time while, simultaneously, siltation reduces reservoir water capacity and availability. Reservoir siltation models are, therefore, mandatory to manage the numerous ungauged small dams in these drylands. The objective of the current study is to improve sediment trap efficiency (TE) modeling by including resuspension as a key probabilistic process. The Shannon entropy was used to estimate the sediment resuspension probability, which was merged with the Camp model and generated the RETSED entropy model. To validate the RETSED model, an experimental check dam (ECD, 300 m³), located in the Gilbués desertification site, Brazil, was monitored hourly during one hydrological year (July 2018–July 2019). Measurements show that the annual volumetric decline of the check dam was 12%; and that the average trap efficiency equaled 86%. Only 9.5% of the hourly sediment concentration outflow exceeded the average plus one standard deviation, showing that the reservoir is well mixed; a fact which highlights the relevance of the resuspension process. Three empirical models failed to mimic the experimental results: Churchill (TE = 99%), Brune (TE = 75%), and Maryland (TE = 94%). According to the RETSED entropy model, the resuspension probability during the experiment was 10% and TE = 81%, a value only 6% below the measured one. The Camp model simulated TE = 89%, only 3% higher than the measured value, but confirmed the Camp model's tendency to overestimate TE due to a disregard of resuspension. The deterministic model showed low sensitivity concerning the hydrodynamic effects, whereas the entropy formulation proved to be more consistent with physical behavior: the resuspension probability clearly increased and trap efficiency decreased with rising reservoir discharge.     

Modelling scour and deposition in harbours due to complex tsunami-induced currents

In this paper, a set of models responsible for hydrodynamics, sediment transport, and morphological evolution are introduced with their theoretical backgrounds, and it is explained how they are fully connected through a two-way coupling to yield an integrated sediment transport model applicable to tsunami cases. In particular, a fully nonlinear Boussinesq model with bottom shear-induced rotational terms is chosen for the hydrodynamic model in order to provide a better physical approximation of tsunami-related, near-bed hydrodynamics in the nearshore. A finite-volume scheme, stable and suitable for phase-resolving model runs longer than 10 simulated hours, is adopted in the numerical discretization. The accuracy and applicability of the developed model are investigated through numerical tests on various sediment problems in the shallow region. Calculated results agree well with existing experimental records. Finally, an ocean-wide, field-scale simulation of the 2011 Tohoku-oki tsunami is attempted, with a focus on the localized effects of tsunami-induced morphological changes at Crescent City Harbor and Santa Cruz Harbor (USA). Consistent with the reported observations, strong and vortical velocity fields are generated through the model and result in significant changes in morphological configurations. Depth variations and areas of scouring and deposition are compared between modelled and observed records, and the results are discussed. © 2019 John Wiley & Sons, Ltd.     

Two-dimensional modeling of fine sediment transport with mixed sediment and consolidation: Application to the Gironde Estuary,France

In order to optimize ship navigation in the macrotidal Gironde Estuary, a recent project funded by the port of Bordeaux aims at better understand and forecast hydrodynamic and fine sediment transport within the estuary. In the framework of this project, a two-dimensional hydro-sedimentary model is built. The model includes hydrodynamic forcings, mixed-sediment transport, and consolidation processes. The harmonic analysis of the astronomical tides reveals a strong distortion of the tidal wave inducing the growth of overtide constituents and the non-significant effect of tide-surge interactions in annual-scale prediction. Depending on hydrological conditions, river discharge can considerably alter the model accuracy due to the migration of the turbidity maximum zone modifying the bottom roughness. Comparison with measurements shows the ability of the model to reproduce suspended-sediment concentrations in the central Estuary. Sensitivity of the model to sediment features has also been discussed in regard of suspended-sediment concentrations and fluid mud deposits. The model will be further coupled with ship squat predictions and a morphodynamic model.     

Numerical modeling of bedload and suspended load contributions to morphological evolution of the Seine Estuary (France)

This numerical modeling study (i) assesses the influence of the sediment erosion process on the sediment dynamics and subsequent morphological changes of a mixed-sediment environment, the macrotidal Seine estuary, when non-cohesive particles are dominant within bed mixtures (non-cohesive regime), and (ii) investigates respective contributions of bedload and suspended load in these dynamics. A three dimensional (3D) process-based morphodynamic model was set up and run under realistic forcings (including tide, waves, wind, and river discharge) during a 1-year period. Applying erosion homogeneously to bed sediment in the non-cohesive regime, i.e., average erosion parameters in the erosion law (especially the erodibility parameter, E₀), leads to higher resuspension of fine sediment due to the presence of coarser fractions within mixtures, compared to the case of an independent treatment of erosion for each sediment class. This results in more pronounced horizontal sediment flux (two-fold increase for sand, +30% for mud) and erosion/deposition patterns (up to a two-fold increase in erosion over shoals, generally associated with some coarsening of bed sediment). Compared to observed bathymetric changes, more relevant erosion/deposition patterns are derived from the model when independent resuspension fluxes are considered in the non-cohesive regime. These results suggest that this kind of approach may be more relevant when local grain-size distributions become heterogeneous and multimodal for non-cohesive particles. Bedload transport appears to be a non-dominant but significant contributor to the sediment dynamics of the Seine Estuary mouth. The residual bedload flux represents, on average, between 17 and 38% of the suspended sand flux, its contribution generally increasing when bed sediment becomes coarser (can become dominant at specific locations). The average orientation of residual fluxes and erosion/deposition patterns caused by bedload generally follow those resulting from suspended sediment dynamics. Sediment mass budgets cumulated over the simulated year reveal a relative contribution of bedload to total mass budgets around 25% over large erosion areas of shoals, which can even become higher in sedimentation zones. However, bedload-induced dynamics can locally differ from the dynamics related to suspended load, resulting in specific residual transport, erosion/deposition patterns, and changes in seabed nature.     

Three-dimensional numerical modeling of cohesive sediment transport and wind wave impact in a shallow oxbow lake

It was observed that in some closed inland lakes sediment transport was dominated by wind-induced currents, and the sediment resuspension was primarily driven by wind-induced waves. This paper presents the development and application of a three-dimensional numerical model for simulating cohesive sediment transport in water bodies where wind-induced currents and waves are important. In the model, the bottom shear stresses induced by currents and waves were calculated, and the processes of resuspension (erosion), deposition, settling, etc. were considered. This model was first verified by a simple test case consisting of the movement of a non-conservative tracer in a prismatic channel with uniform flow, and the model output agreed well with the analytical solution. Then it was applied to Deep Hollow Lake, a small oxbow lake in Mississippi. Simulated sediment concentrations were compared with available field observations, with generally good agreement. The transport and resuspension processes of cohesive sediment due to wind-induced current and wave in Deep Hollow Lake were also discussed.     

Modelling of hydrodynamics and cohesive sediment transport in Tanshui River estuarine system,Taiwan

A laterally averaged two-dimensional numerical model is used to simulate hydrodynamics and cohesive sediment transport in the Tanshui River estuarine system. The model handles tributaries as well as the main stem of the estuarine system. Observed time series of salinity data and tidally averaged salinity distributions have been compared with model results to calibrate the turbulent diffusion coefficients. The overall model verification is achieved with comparisons of residual currents and salinity distribution. The model reproduces the prototype water surface elevation, currents and salinity distributions. Comparisons of the suspended cohesive sediment concentrations calculated by the numerical model and the field data at various stations show good agreement. The validated model is applied to investigate the tidally averaged salinity distributions, residual circulation and suspended sediment concentration under low flow conditions in the Tanshui River estuarine system. The model results show that the limit of salt intrusion in the mainstem estuary is located at Hsin-Hai bridge in Tahan Stream, 26 km from the River mouth under Q75 flow. The null point is located at the head of salt intrusion, using 1 ppt isohaline as an indicator. The tidally averaged sediment concentration distribution exhibits a local maximum around the null point.     

Numerical modeling of surface runoff and erosion due to moving rainstorms at the drainage basin scale

A physically-based distributed erosion model (MEFIDIS) was applied to evaluate the consequences of storm movement on runoff and erosion from the Alenquer basin in Portugal. Controlled soil flume laboratory experiments were also used to test the model. Nine synthetic circular storms were used, combining three storm diameters (0.5, 1 and 2 times the Alenquer basin’s axial length) with three speeds of storm movement (0.5, 1 and 2 m/s); storm intensities were synthesized in order to maintain a constant rainfall depth of 50 mm. The model was applied to storms moving downstream as well as upstream along the basin’s axis. In all tests, downstream-moving storms caused significantly higher peak runoff (56.5%) and net erosion (9.1%) than did upstream-moving storms. The consequences for peak runoff were amplified as the storm intensity increased. The hydrograph shapes were also different: for downstream-moving storms, runoff started later and the rising limb was steeper, whereas for upstream moving storms, runoff started early and the rising limb was less steep. Both laboratory and model simulations on the Alenquer basin showed that the direction of storm movement, especially in case of extreme rainfall events, significantly affected runoff and soil loss.     

Modelling and experimental study of coupled exchange of p,p′‐DDE and colloids in a flume simulated stream–streambed system

Stream–subsurface exchange plays a significant role in the fate and transport of contaminants in streams. It has been modelled explicitly by considering fundamental processes such as hydraulic exchange, colloid filtration, and contaminant interactions with streambed sediments and colloids. The models have been successfully applied to simulate the transport of inorganic metals and nutrients. In this study, laboratory experiments were conducted in a recirculating flume to investigate the exchange of a hydrophobic organic contaminant, p,p′‐dichloro‐diphenyl‐dichloroethane (DDE), between a stream and a quartz sand bed. A previously developed process‐based multiphase exchange model was modified by accounting for the p,p′‐DDE kinetic adsorption to and desorption from the bed sediments/colloids and was applied to interpret the experimental results. Model input parameters were obtained by conducting independent small‐scale batch experiments. Results indicate that the immobilization of p,p′‐DDE in the quartz sand bed can occur under representative natural stream conditions. The observed p,p′‐DDE exchange was successfully simulated by the process‐based model. The model sensitivity analysis results show that the exchange of p,p′‐DDE can be sensitive to either the sediment sorption/desorption parameters or colloidal parameters depending on the experimental conditions tested. For the experimental conditions employed here, the effect of colloids on contaminant transport is expected to be minimal, and the stream–subsurface exchange of p,p′‐DDE is dominated by the interaction of p,p′‐DDE with bed sediment. The work presented here contributes to a better mechanistic understanding of the complex transport process that hydrophobic organic contaminants undergo in natural streams and to the development of reliable, predictive models for the assessment of impacted streams. Copyright © 2015 John Wiley & Sons, Ltd.     

Experimental investigation of the morphodynamic response of riffles and pools to unsteady flow and increased sediment supply

Gravel-bed rivers characteristically exhibit shallow riffles in wide sections and deeper pools where the channel becomes constricted and narrow. While rivers can adjust to changing flow and sediment supply through some combination of adjustments of channel slope, bed-surface sorting, and channel shape, the degree to which riffle-pools may adopt these changes in response to changing flows and sediment supplies remains unclear. This article presents results from a flume experiment investigating how constant- and variable-width channels adjust their morphology in response to changing flow and increased sediment supply. Two flume geometries were used: (1) constant-width and (2) variable-width, characterized by a sinusoidal pattern with a mean width equal to that of the first channel. The variable-width channel developed bed undulations in phase with the width, representing riffle-pools. The experiment consisted of three phases for each flume geometry: (1) steady flow, constant sediment supply; (2) unsteady flow, constant sediment supply; and (3) unsteady flow, doubled sediment supply. Unsteady flow was implemented in the form of repeated symmetrical stepped hydrographs, with a mean discharge equal to that in the steady flow phase. In all phases the bed and sediment supply were composed of a sand/gravel mixture ranging from 1 to 8 mm. In both the straight and variable-width channels, transitioning from steady flow to repeated hydrographs did not result in significant changes in bed morphology. The two channel geometries had different responses to increased sediment supply: the slope of the straight channel increased nearly 40%, while the variable-width channel reduced the relief between bars and pools and decreased the variability in cross-sectional elevation with a slight slope increase. Bar-pool relief varied with repeat discharge hydrographs. Pool elevation changed twice the distance of bar elevations, emphasizing the relevance of pool scour for riffle-pool self-maintenance in channels with width variations.     

鲫(Carassius auratus)对浅水湖泊水体浊度、营养水平和浮游生物生物量的影响——基于中宇宙模拟实验

鲫(Carassius auratus)是我国各类淡水水体的优势鱼类之一.作为底栖杂食性鱼类,一方面,鲫可以通过排泄和扰动沉积物影响湖泊营养和光照水平,通过"上行效应"促进浮游植物生长;另一方面,鲫也可以捕食浮游动物,通过"下行控制"影响藻类生长以及营养盐循环.对于浅水湖泊,两种途径对于生态系统影响的相对重要性仍有待研究.本研究设计了一个两因素户外中宇宙实验,通过在沉积物表面添加隔网的方式,比较两种情况下(能、否接触沉积物),鲫对水体浊度、营养盐和浮游生物生物量的影响.实验在16个大型钢化玻璃桶(400 L)中进行,持续36 d(2019年8月6日—9月11日).研究结果表明:1)在能接触沉积物的条件下,鲫显著促进了沉积物再悬浮,表现为水体的总悬浮物(TSS)和无机悬浮物(ISS)浓度大幅升高;水体的光衰减系数(Kd)增加,总氮(TN)和总磷(TP)浓度明显升高; 2)在不能接触沉积物的条件下,鲫对水体悬浮物(TSS和ISS)浓度和Kd的影响不明显,但是显著降低了水体TN和TP浓度; 3)在两种情况下,鲫对浮游植物叶绿素a浓度以及浮游动物生物量的影响均不显著.本研究表明鲫只有在能够接触沉积物的条件下,才会显著提高水体浊度和营养水平.因此,在缺乏沉水植被的浅水湖泊中,鲫扰动沉积物产生的"上行效应"可能是其对生态系统产生负面影响的主要途径.     

Infiltration through series of soil aggregates: Neutron radiography and modeling

Soils are often structured as fine-porous aggregates separated by large inter-aggregate pores. Under unsaturated conditions, water is mostly stored in the aggregates and water flow depends on the properties of the aggregates as well as on those of the contacts between aggregates. The goal of this study is to model and evaluate the hydraulic properties of the contacts. We used neutron radiography to monitor the infiltration of water through series of aggregates. The flow process was numerically simulated by considering the hydraulically conducting contact area between aggregates as a variable that depends on the capillary pressure. This contact area was evaluated by matching the observed and simulated water flow across aggregates. We determined the conductivity of the contacts assuming that it scales with the contact area. We also measured the equivalent conductivity of series of aggregates. We found that during drainage the hydraulically conducting contact area drastically decreases and the conductivity of the contacts becomes much smaller than that of the aggregates. We also found that the equivalent conductivity of the aggregate series decreases as the conductivity of the contacts. We concluded that the contacts control the flow: they are highly conductive when wet, but act as bottle-necks under drained conditions. The abrupt transition between these two limiting cases indicates that the contact region is more rapidly drained than the aggregate interior. Our findings might be extended to describe water transfer to an evaporating boundary, infiltration and storage of water in unsaturated aggregated soils.