Observation of bioturbation and hyporheic flux in streambeds

In the Elkhorn River, burrows, tubes, and sediment mounds created by invertebrate bioturbation were observed in the exposed streambed and commonly concentrated on the fine-sediment patches, which consist of silt, clay, and organic matter. These invertebrate activities could loosen the thin layer of clogging sediments and result in an increase of pore size in the sediments, leading to greater vertical hydraulic conductivity of the streambed (K _v ). The measurements of the vertical hydraulic gradient across the submerged streambed show that vertical flux in the hyporheic zone can alter directions (upward versus downward) for two locations only a few meters apart. In situ permeameter tests show that streambed K _v in the upper sediment layer is much higher than that in the lower sediment layer, and the calculated K _v in the submerged streambed is consistently greater than that in the clogged sediments around the shorelines of the sand bars. Moreover, a phenomenon of gas bubble release at the water-sediment interface from the subsurface sediments was observed in the groundwater seepage zone where flow velocity is extremely small. The bursting of gas bubbles can potentially break the thin clogging layer of sediments and enhance the vertical hydraulic conductivity of the streambed.     

Study of flow formed by three coplanar impinging pipe jets at inclination angles of 30° and 45°

This paper presents an experimental study of the interactions of three fully-submerged, coplanar impinging jets issued from long pipes. The jets were neutrally buoyant and were arranged symmetrically about the axis of a central jet, with two side jets set to intersect with the central jet at two inclination angles (30° and 45°) and three Reynolds numbers (4240, 6400 and 8000). Measurements of the flow fields were performed using particle image velocimetry to examine the flow structures in various planes, i.e., the jet axis plane (X–Y), the jet normal plane (X–Z) and the cross-sectional plane (Y–Z). This flow configuration results in pronounced interactions among the three jets, and hence better mixing than that of a canonical single pipe jet as illustrated by augmented centreline velocity decay, spreading rate and turbulence level. The jets at the inclination angle of 45° impinge and mix more rapidly than those at 30°. For each case, the flow is highly 3-dimensional, and jet development displayed several distinct regions (converging, merging and combining) along the streamwise direction. The expansion of flow in the X–Y plane is similar to the shape of a hyperbola with necking formed immediately downstream of the impinging point, whereas that in the X–Z plane assumes the shape of a parabola with an open rim exhibiting a pronounced velocity deficit in the central part of the combining flow. Self-similarity of streamwise mean velocity is explored in the combining region of the flow on the two planes of symmetry (X–Y and X–Z). Flow development in the combining region is dominated by large-scale vortical structures, including von Kárman-like spanwise vortices in the X–Y plane and secondary circulation in the form of two pairs of counter-rotating streamwise vortices in the Y–Z plane.     

Direct measurements of the drag force over aligned arrays of cubes exposed to boundary-layer flows

Wind tunnel measurements of the total drag force for aligned arrays of cubes exposed to two different boundary-layer flows at three flow velocities are discussed. The drag force for eight different building packing densities λ _p (from 0.028 to 1) is measured with a standard load cell generating a novel dataset. Different λ _p are reproduced by increasing the number of buildings on the same lot area; this represents a real situation that an urban planner is faced with when a lot area of a given (fixed) size is allocated to the development of new built areas. It is assumed that the surrounding terrain is uniform and there is a transition from a given roughness (smooth) to a new roughness (rough). The approaching flow will adjust itself over the new surface within a distance that in general may be larger than the horizontal length covered by the array. We investigate the region where the flow adjustment occurs. The wide range of packing densities allowed us to analyse in detail the evolution of the drag force. The drag force increases with increasing packing densities until it reaches a maximum at an intermediate packing density (λ _p  = 0.25 in our case) followed by a slight decrease at larger packing densities. The value of the drag force depends on the flow adjustment along the array which is evaluated by introducing the parameter “drag area” to retrieve information about the drag distribution at different λ _p . Results clearly suggest a change of the distribution of the drag force, which is found to be relatively uniform at low packing densities, while most of the force acts on first rows of the arrays at large packing densities. The drag area constitutes the basis for the formulation of a new adjustment length scale defined as the ratio between the volume of the air within the array and the drag area. The proposed adjustment length scale automatically takes into account the change in drag distribution along the array for a better parameterization of urban effects in dispersion models.     

Hydrodynamics and turbulence in emergent and sparsely vegetated open channel flow

This present study reports the results of an experimental study characterizing thorough variation of turbulent hydrodynamics and flow distribution in emergent and sparsely vegetated open channel flow. An emergent and rigid sparse vegetation patch with regular spacing between stems along the flow and transverse directions was fixed in the central region of the cross-section of open channel. Experiments were conducted in subcritical flow conditions and velocity measurements were obtained with an acoustic Doppler Velocimetry system. Large variations of the turbulence intensities, Reynolds shear stress, turbulent kinetic energy and vortical motions are found in and around the vegetation patch. At any cross-section through the interior of the vegetation patch, streamwise velocity decreases with increase in streamwise length and the velocity profiles converge from the log-law to a linear profile with increasing slope. Time-averaged lateral and vertical velocities inside the vegetation patch increase with increasing streamwise distance and converge from negative values to positive values. Turbulence intensities interior of the sparse vegetation patch are more than those of without the vegetation patch. Similar to the trend of streamwise velocity profiles inside the vegetation, turbulence intensities and longitudinal-normal Reynolds shear stress profile decreases with streamwise direction. In the interior of the vegetation patch and downstream of the trailing edge, turbulent kinetic energy profiles are exhibiting irregular fluctuations and the maximum values are occurring in the outer layer. Analysis of flow distribution confirms sparse vegetation patch is inducing a serpentine flow pattern in its vicinity. At the leading edge, flow is rushing towards the right hand sidewall, and at the trailing edge, flow is turning to the left hand sidewall. In between the leading and trailing edges, the streamlines are following a zig-zag fashion at varied degree along the streamwise and lateral directions. Immediate upstream of the leading edge and in the interior of the vegetation patch, vortex motion is clearly visible and the vortices are stretched along the width of the channel with streamwise direction.     

Numerical simulation and optimization study of the wind flow through a porous fence

Three turbulence closure models (RNG k-ε, SST k-ω and RSM) were used to investigate the flow characteristics around a two-dimensional isolated porous fence. The comparison between the numerical results and the experimental measurements indicated that RSM model shows a better performance than the other two models. The aim of this paper is to accurately and efficiently determine the optimum porosity that attain the best shelter effect of the wind fence in the near wake region (0–4h_b) and in the far wake region (4h_b–10h_b) respectively, where h_b is the height of the fence. The gradient algorithm was adopted as the optimization algorithm and the RSM model was used to model turbulent features of the flow. The shelter effect was parameterized by the peak velocity ratio involving velocity and turbulence. The objective was to reduce the peak velocity ratio in the near or far wake region by changing the design variable porosity (?) of the fence, which ranged between 2 and 60%. The results revealed that a porosity of 10.2% was found as the optimum value giving rise to the best shelter effect in the near wake region, and ? = 22.1% was determined in the case of the far wake region. In addition, based on the proposed optimization method, it is found that the recirculating bubble behind the fence can only be detected when ? < 29.9%.     

Empirical equation for transverse dispersion coefficient based on theoretical background in river bends

There are different approaches to estimating the transverse dispersion coefficient in river mixing. Theoretical approaches have derived the dispersion coefficient from the concept of shear flow, which has dominant effects on the transverse mixing. Empirical approaches have developed an equation using the hydraulic and geometric data of rivers through dimensional analysis and regression techniques. These two equations interact closely with each other. For example, the complicated theoretical equation can be simplified by empirical approaches, and the functional relationships of the empirical equation can be derived from theoretical bases. In this study, a new empirical equation for the transverse dispersion coefficient has been developed based on the theoretical background in river bends. As a regression method, the least-square iterative method was used because the equation was a nonlinear model. The estimated dispersion coefficients derived by the new equation were compared with observed transverse dispersion coefficients acquired from natural rivers and coefficients calculated by the other existing empirical equations. From a comparison of the existing transverse dispersion equations and the proposed equation, it appears that the behavior of the existing formula in a relative sense is very much dependent on the flow condition and the river geometry. Moreover, the proposed equation does not vary widely according to variation of flow conditions. Also, it was revealed that the equation proposed in this study becomes an asymptotic curve as the curvature effect increases.     

Compound channel flow with a longitudinal transition in hydraulic roughness over the floodplains

Flows in a compound open-channel (two-stage geometry with a main channel and adjacent floodplains) with a longitudinal transition in roughness over the floodplains are experimentally investigated in an 18 m long and 3 m wide flume. Transitions from submerged dense vegetation (meadow) to emergent rigid vegetation (wood) and vice versa are modelled using plastic grass and vertical wooden cylinders. For a given roughness transition, the upstream discharge distribution between main channel and floodplain (called subsections) is also varied, keeping the total flow rate constant. The flows with a roughness transition are compared to flows with a uniformly distributed roughness over the whole length of the flume. Besides the influence of the downstream boundary condition, the longitudinal profiles of water depth are controlled by the upstream discharge distribution. The latter also strongly influences the magnitude of the lateral net mass exchanges between subsections, especially upstream from the roughness transition. Irrespective of flow conditions, the inflection point in the mean velocity profile across the mixing layer is always observed at the interface between subsections. The longitudinal velocity at the main channel/floodplain interface, denoted \(U_{int}\), appeared to be a key parameter for characterising the flows. First, the mean velocity profiles across the mixing layer, normalised using \(U_{int}\), are superimposed irrespective of downstream position, flow depth, floodplain roughness type and lateral mass transfers. However, the profiles of turbulence quantities do not coincide, indicating that the flows are not fully self-similar and that the eddy viscosity assumption is not valid in this case. Second, the depth-averaged turbulent intensities and Reynolds stresses, when scaled by the depth-averaged velocity \(U_{d,int}\) exhibit two plateau values, each related to a roughness type, meadow or wood. Lastly, the same results hold when scaling by \(U_{d,int}\) the depth-averaged lateral flux of momentum due to secondary currents. Turbulence production and magnitude of secondary currents are increased by the presence of emergent rigid elements over the floodplains. The autocorrelation functions show that the length of the coherent structures scales with the mixing layer width for all flow cases. It is suggested that coherent structures tend to a state where the magnitude of velocity fluctuations (of both horizontal vortices and secondary currents) and the spatial extension of the structures are in equilibrium.     

Semi-parametric Bayesian density estimation using ranked set sample in the presence of ranking error

In this paper, we propose a Bayesian method to estimate the underlying density function of a study variable Y using a ranked set sample in which an auxiliary variable X is used to rank the sampling units. The amount of association between X and Y is not known, resulting in an unknown degree of ranking error. We assume that (X, Y) follows a Morgenstern family of distributions. The study variable Y is assumed to have a parametric distribution, with the distribution of the parameters having a Dirichlet process prior. A Markov chain Monte Carlo procedure is developed to obtain a Bayesian estimator of the desired density function as well as of the ranking error. A simulation study is used to evaluate the performance of the proposed method. An example from forestry is used to illustrate a real-life application of the proposed methodology.     

Diffusional mass transfer coefficient at the water–sediment interface for wind-induced flow in very shallow lagoons

Very shallow lagoons that are a few centimeters deep are common in the arid Andes of Northern Chile, Argentina, Bolivia and Perú. The dynamics of these lagoons are dominated by the water–sediment interface (WSI) and strong afternoon winds. Although many studies have examined the diffusional mass transfer coefficients (k _t ) of open channel flows, estimates for wind-induced flows are still unknown. The aim of this article is to propose and validate an analytical expression for computing k _t at the WSI for wind-induced flow. The laboratory measurements were conducted in a wind tunnel with a water tank of variable depth located at its downwind end. Natural muddy sediments were placed in the middle of the tank so that the dissolved oxygen (DO) was consumed in the sediments. The diffusional mass transfer coefficient that characterizes the DO uptake in the sediment was obtained from DO micro-profiles measured with an OX-25 Unisense microelectrode. Water velocity profiles were measured with a 2D side-view Sontek acoustic doppler velocimetry (ADV), and the wind shear velocity was computed based on wind velocity profiles that were measured with an Extech hot-wire anemometer. A total of 16 experiments were conducted with different water depths and wind shear stresses. The constants required by the model were determined from these experiments, and the analytical expression was successfully validated by the laboratory observations. The analytical expression obtained for computing k _t was also validated with field observations that were conducted in October, 2012, in Salar del Huasco, Northern Chile (20.274° S, 68.883° W, 3800 m above sea level). The comparison between the observed and predicted values of k _t provides a determination coefficient of r ² = 0.48 and a p value < 0.01. The results show that the value of k _t for wind-induced flow is proportional to the wind shear velocity and the inverse of the Reynolds number of the wind-induced current.     

Flow and drag force around a free surface piercing cylinder for environmental applications

This paper investigates flows around a free surface piercing cylinder with Froude number F > 0.5 and Reynolds number around Re = 50,000. The aim of this work is to gain a better understanding of the flow behaviour in environmental systems such as fishways. The advances are based upon experimental and numerical results. Several flow discharges and slopes are tested to obtain both subcritical and supercritical flows. The drag force exerted on the cylinder is measured with the help of a torque gauge while the velocity field is obtained using particle velocimetry. For the numerical part, two URANS turbulence models are tested, the k-\(\omega\) SST and the RNG k-\(\varepsilon\) models using the OpenFOAM software suite for subcritical cases, and then compared with the corresponding experimental results. With fishways applications in mind, the changes in drag coefficient \(C_d\) versus Froude number and water depth are studied and experimental correlations proposed. We conclude that the most suitable URANS turbulence model for reproducing this kind of flow is the k-\(\omega\) SST model.     

Experimental and numerical study on the shear velocity distribution along one or two dunes in tandem

This paper investigates, experimentally and numerically, the shear velocity distribution along a single transverse dune and along two closely spaced dunes, analyzing the flow effects of one dune upon the other. The paper focuses on two-dimensional models simulating transverse sand dunes. The shape of the two pile geometries studied is described by sinusoidal curves, one having a maximum slope of $32^{\circ }$ and the other $27.6^{\circ }$ , with leeward flow separation. The tests were carried out for two undisturbed wind speeds and the experimental data obtained through wind-tunnel modeling encompass flow visualization and shear-velocity results. A generally good agreement is observed between the experimental measurements and computational results. From the inquiry between shear velocity distributions and published eroded contours for the same geometries, it appears the Bagnold’s approach is insufficient in the prediction of threshold conditions in wake flows formed in the dune’s leeward side.     

Effects of a dynamic membrane formed with polyethylene glycol on the ultrafiltration of natural organic matter

The formation of a dynamic membrane (DM) was investigated using polyethylene glycol (PEG) (molecular weight of 35000 g/mol, concentration of 1 g/L). Two natural organic matters (NOM), Dongbok Lake NOM (DLNOM) and Suwannee River NOM (SRNOM) were used in the ultrafiltration experiments along with PEG. To evaluate the effects of the DM with PEG on ultrafiltration, various transport experiments were conducted, and the analyses of the NOM in the membrane feed and permeate were performed using high performance size exclusion chromatography, and the effective pore size distribution (effective PSD) and effective molecular weight cut off (effective MWCO) were determined. The advantages of DM formed with PEG can be summarized as follows: (1) PEG interferes with NOM transmission through the ultrafiltration membrane pores by increasing the retention coefficient of NOM in UF membranes, and (2) low removal of NOM by the DM is affected by external factors, such as pressure increases during UF membrane filtration, which decreases the effective PSD and effective MWCO of UF membranes. However, a disadvantage of the DM with PEG was severe flux decline; thus, one must be mindful of both the positive and negative influences of the DM when optimizing the UF performance of the membrane.     

Drag coefficient for rigid vegetation in subcritical open-channel flow

Drag coefficient has been commonly used as a quantifying parameter to represent the vegetative drag, i.e., resistance to the flow by vegetation. In this study, the measured data on the drag coefficient for rigid vegetation in subcritical open-channel flow reported in previous studies are collected and preprocessed for multi-parameter analysis. The effect of Froude number (Fr) on the drag coefficient for rigid vegetation in subcritical flow cannot be ignored, especially when \(Fr < 0.12\). The drag coefficient is observed to exponentially decrease with the stem Reynolds number (R _d ) and logarithmically decreased with the vegetation density (λ) when \(0.012 < \lambda < 0.12\). The relative submergence (h ^* ) has a significant effect on the drag coefficient, and a positive logarithmic relationship is summarized. A simplified three-stage empirical formula is obtained based on the divisions of Fr. Laboratory tests (with \(Fr < 0.02\)) prove that the present empirical model has higher precision compared with existing models.     

Multiple lineages and absence of panmixia in the “circumpolar” crinoid <Emphasis Type="Italic">Promachocrinus kerguelensis</Emphasis> from the Atlantic sector of Antarctica

Despite considerable interest in physiology, evolution and life history of Antarctic marine invertebrates, only a limited number of studies have examined the genetic variability and diversity patterns of these organisms. Moreover, understanding and characterizing patterns of Antarctic biodiversity has taken on a degree of urgency because of potential impacts of global warming. To expand an understanding of the evolutionary history of Antarctic marine invertebrates, the genetic diversity of the crinoid Promachocrinuskerguelensis Carpenter, 1888 was investigated, which is documented to have a circumpolar distribution extending to subantarctic islands. Specimens of P. kerguelensis were collected from the western side of the Antarctic Peninsula, and the subantarctic islands South Georgia, South Sandwich and Bouvetøya Island from 2001 to 2004. P. kerguelensis was previously subject to morphological review that confirmed the taxonomic recognition of only one species. The wide distribution and reported high dispersal capability for P. kerguelensis predicts one large panmictic population. In contrast, nucleotide sequence data from mitochondrial cytochrome oxidase subunit I and cytochrome b genes, collected herein, reveal distinct genetic structure and cryptic speciation within P.kerguelensis. In the Antarctic Atlantic sector alone, there were at least five “species-level” clades. Some of these clades are geographically limited, and most exist in sympatry. The largest and most widespread of these clades was examined to help elucidate connectivity along the subantarctic islands of the Scotia Arc and the Antarctic Peninsula. Within this clade, most genetic diversity was contained within populations, but significant differences were present between regions (Antarctic Peninsula, South Sandwich Is., South Georgia, Bouvetøya Is.), suggesting a corresponding lack of gene flow. Given that P. “kerguelensis” is a well-studied taxon, the finding of considerable genetic diversity within the Atlantic sector alone suggests that the recognized diversity of Antarctica’s benthic marine life may be underestimated, and will rise dramatically with phylogeographic analyses of putative widespread species.     

A two-dimensional inviscid model of the gravity-current head by Lagrangian block simulation

A two-dimensional inviscid model of the gravity-current head produced by the release of a relatively small volume of dense fluid from behind a tall lock gate is constructed by Lagrangian block simulation. Three numerical experiments are conducted for the lock’s height-to-length aspect ratios H/L _o  = 8, 4 and 2. The front speeds obtained by the simulations agree with the laboratory observation for a similar range of aspect ratios. The floor velocity in the wake behind these heads is found to be greater than their front speed. The high floor velocity is caused by the impingement of the coherent wake vortex on the floor. It is a condition that permits these gravity-current heads to maintain their structural integrity so that the fine sediments can travel with the head over long distances on the ocean floor. The structural coherence of the current head depends on the lock aspect ratio. The gravity-current head produced by the release from the lock with the highest aspect ratio of H/L _o  = 8 is most coherent and relatively has the greatest floor velocity and the least trailing current behind the head.     

Changes of microbial composition during wastewater reclamation and distribution systems revealed by high-throughput sequencing analyses

This study employed 454-pyrosequencing to investigate microbial and pathogenic communities in two wastewater reclamation and distribution systems. A total of 11972 effective 16S rRNA sequences were acquired from these two reclamation systems, and then designated to relevant taxonomic ranks by using RDP classifier. The Chao index and Shannon diversity index showed that the diversities of microbial communities decreased along wastewater reclamation processes. Proteobacteria was the most dominant phylum in reclaimed water after disinfection, which accounted for 83% and 88% in two systems, respectively. Human opportunistic pathogens, including Clostridium, Escherichia, Shigella, Pseudomonas and Mycobacterium, were selected and enriched by disinfection processes. The total chlorine and nutrients (TOC, NH₃-N and NO₃-N) significantly affected the microbial and pathogenic communities during reclaimed water storage and distribution processes. Our results indicated that the disinfectant-resistant pathogens should be controlled in reclaimed water, since the increases in relative abundances of pathogenic bacteria after disinfection implicate the potential public health associated with reclaimed water.     

Steady Shallow-Water Current and Solute Transport around a Semi-Conical Headland

Laboratory experiments have been carried out to investigate the effects of a sloping wall headland on the flow characteristics and the associated concentration distributions from a point source around the headland. A semi-conical headland with a slope of 1:2 was set up in a flow basin, 4.8 m long and 3.8 m wide. In this paper, the experimental results of a steady shallow-water current are reported. Three dimensional flow velocities in the basin were measured using Sontek-ADV instrument. The dye concentration levels in the basin were measured by two fluorometers. The experimental results showed a large-scale re-circulation region behind the semi-conical headland. The peak turbulence energy, at about 53% of the local kinetic flow energy, coincides with the region of high velocity gradient. Significant vertical flows were observed around the area near the downhill slope of the headland, with a maximum ratio of vertical to horizontal velocities being about 22%. Such relatively significant vertical scouring velocities, coupled with strong turbulence energy and high horizontal velocity gradients in the same region, could cause severe bed erosion. The experimental results have also been compared with the predicted results of a depth-averaged numerical model. The predicted eddy structure and the concentration distribution in the re-circulation area were found to compare favourably with the experimental results. However, the discrepancies in the flow velocities and the concentration levels near the headland were apparent. It was observed that the dye concentration continued to spread in the cross-stream direction after passing the headland, whereas only a limited extent of the lateral spreading was predicted by the numerical model further downstream of the headland.     

Ancient settlement dynamics and predictive archaeological models for the Metapontum coastal area in Basilicata,southern Italy: from geomorphological survey to spatial analysis

In the ancient territory of the Greek settlement of Metaponto and its Chora (Ionian sector of Basilicata, southern Italy) geomorphological mapping, GIS-supported statistics, and landform topographic features were investigated and combined to extract settlement rules and site dynamics. Analyses of environmental dynamics as well as of spatial and temporal settlement evolution were carried through an integrated approach that, starting from a detailed geomorphological analysis, tried to extract the spatial relationships between archaeological site locations and landform/landscape features. Spatial analysis was used to investigate the relationships between environmental parameters and archaeological sites. The advantages to develop a predictive model using Archeological Spatial Evaluation (ASE) analysis lay in the possibility to deeply study the relationships between man and the environment. In the present study, the analysis was carried out in a pioneering way in the ~400 km² Metaponto territory in order to robustly validate - in the aftermath data - human and territorial dynamics. A model describing interactions between sites and parameters such as elevation, slope, aspect, landforms, land use, and distances from rivers was constructed. Sensibility maps were produced, which will help archaeologists to know how many and which are the parameters increasing the probability to find new archaeological sites as well as how the main geomorphological features of the study area, contextualized within the most likely paleoclimate scenario, affect ancient site arrangement though time. The preferential occupation of mid-altitude marine terraces, and the consequence spreading of agriculture on these territories, is likely due to the existence of well-developed soil profiles on them. Clearly colonists recognized that on these landform units there were the better conditions for the development of massive agricultural practices: since the SP age (4,1 to 2,0 Kyr B.P.) the settlers expands inland looking for most fertile territories as well as for areas sufficiently protected (being far away from the coast). The increase in farmhouses on the top of marine terraces from the FC up to EC-EE periods and the modifications of settlement distribution (gradual abandonment of alluvial) in the Metaponto area is likely related to the acceleration of alluvial processes. The progressive decrease of human occupation during the Hellenistic up to the Roman Age is clearly consequence of the Roman conquest. However, a role played by the increase in flooding occurrence in the coastal plain/floodplains of the main rivers in triggering the abandonment of these territories and marking the beginning of the decadent phases should not be neglected. Productive areas (locally associated to farmhouses) are preferential set along the fluvial incisions (sometime associated to farmhouses), thus implying that productive area setting is strongly linked to the presence of rivers, and of their lower rank tributaries, both as an intrinsic need of manufacturing and to facilitate the spreading of products. From the Classical to the Hellenistic Age until the Roman Age, Necropolis are spread in the entire area, from marine terraces to alluvial plain. On the other hand, Sacred areas - which appearance is recorded starting from the Colony Foundation Age - exhibits a settlement continuity in the study area starting from this time up to the entire Roman period, and the sacred areas setting is strictly associated to the category of Settlement.     

Characteristics of residual organochlorine pesticides in soils under different land-use types on a coastal plain of the Yellow River Delta

The residual levels of organochlorine pesticides (OCPs) were examined in soils covering five types of land use along a salinity gradient on the Yellow River Delta. The most prominent OCPs were dichlorodiphenyltrichloroethane (∑DDT, arithmetic mean = 5.11 μg kg^?1), hexachlorocyclohexane (∑HCH, 1.69 μg kg^?1) and ∑endosulfan (10.4 μg kg^?1). The spatial variability of OCPs composition shifted from γ-HCH and o,p′-DDT dominated pesticides in coastal soils to p,p′-DDE dominated pesticides in inland soils. In different land-use types, the percentages of β-HCH and p,p′-DDE are characterized by more recalcitrant components in decreasing order of vegetable fields, cereal fields, cotton fields, wetlands and tidal flats with increasing soil salinity. However, the less recalcitrant components, γ-HCH and o,p′-DDT, showed an opposite trend. Endosulfan sulfate predominated in all land-use types. Residual levels of β-HCH were affected by soil organic matter. The correlations between γ-HCH and clay content and between p,p′-DDE, o,p′-DDT and salinity might associate with the influence of sediment cotransport by the Yellow River and the density of anthropogenic activities in coastal region. Depth distribution of the OCPs in typical soil profiles also implied that local historical usage and sediment transport by the Yellow River both affected the OCPs residual in this region.     

Study on the concentration distribution in a trapezoidal open-channel flow with a side discharge

The number and distribution of pollutant concentration in a trapezoidal open channel flow with a side discharge is calculated and effects of the bank gradient are investigated in this paper. A sigma-coordinate water quality numerical model is used to simulate the process of both water and pollutant transportation in the trapezoidal channel open flow. The diffusion coefficient used in the prediction is determined by two methods including constant coefficient and the depth-averaged k-epsilon turbulence closure model. The change of the concentration with the bank gradient is acquired based on the simulation of cases with different bank gradients. An analytical formula is derived by using the mirror image method and related diffusion theories, ignoring the discharge momentum and the influence of the opposite bank. The formula can predict the number and distribution of pollutant concentration with some acceptable errors. The results demonstrate that the bank gradient has great influence on the concentration distribution which will decrease with the increase of the bank gradient approximately following a hyperbolic law.