Methodological approaches for delimitating mixing zones in rivers: establishing admissibility criteria and flow regime representation

The present paper addresses defining the extent of the mixing zones of effluents discharged into rivers, which is a problem that should be analyzed based on probabilistic terms, as it is governed by several random processes. A river’s flow regime is one of the main variables, and it has a high dependence on hydrological processes. Additionally, after calculating the extent of the mixing zone, it is necessary to determine if the resulting dimensions are admissible or not. Common practice includes the admissibility criteria associated with the geometry of the river. However, this practice does not consider the environmental characteristics of the river that make it capable of absorbing the impact of the effluent (the biological and hydromorphological status of the river and the presence of structures that can change the river flow conditions, ecologically sensitive area or threatened species). This paper presents work on two important topics: on the one hand, the development of a methodology to establish the admissible extent of the mixing zone as a function of the environmental vulnerability of the river to the discharged effluent and, on the other hand, the proposal of a procedure to perform the calculations of the effluent mixing considering the variability of the river’s flow regime. The proposed methodological approaches are illustrated with an application to a real case, including a numerical simulation of the hydrodynamic and effluent evolution of the river during a year, to test the proposed methodology’s suitability and demonstrate the important savings in computational effort that can be achieved.     

Chemical composition of rivers and streams in the forest and savanna zones of Nigeria

The aim of this study is to investigate the differences in the chemical conditions of lotic waterbodies in the two major ecosystems in Nigeria, the forest and savanna zones. The forest waters were slightly acidic (mean±SD pH = 6.72±0.58) while the savanna waters were slightly alkaline (pH = 7.11±0.33). The cationic order of dominance in the forest waters was Na⁺ > Ca²⁺ > Mg²⁺ > K⁺ in contrast to Ca²⁺ > Mg²⁺ > Na⁺ > K⁺ in savanna waters. The forest waters were chlorided (typical of coastal and/or marine waters) whereas the savanna waters were carbonated in nature, typical of the worldwide freshwater. Organic carbon was significantly higher in forest waters than in the savanna waters (p < 0.05) while nutrient compounds were significantly higher in savanna waters than in forest waters. The seasonal variation of the chemical parameters was generally more evident in savanna than in forest waters. The differences in water quality between the two major vegetation zones reflect the differences in the biogeochemical processes and nutrient cycling that characterise forest and savanna ecosystems.     

Cambial zones in gills of Bivalvia

Localization and intensity of cell division in the gills of 9 species of Bivalvia were studied. In addition, deoxyribonucleic acid (DNA) synthesis in the gills of two species of mussels (Crenomytilus grayanus and Mytilus edulis) and of the scallop Mizuhopecten yessoensis was investigated by means of radioautography using thymidine-H³. In all bivalves examined, the distribution of branchial cell division was very uneven. Distinctly visible zones exhibited considerably higher mitotic and DNA-synthesizing activities than the remaining branchial epithelium. Topography of these zones appeared similar in all Filibranchia and Eulamellibranchia examined. Areas with enhanced cell division in the gills of Bivalvia are considered as cambial zones or points of growth, apparently serving as a source of growth and/or renewal of the entire organ.     

Models predicting the fate of pollutants in rivers

The spill of 2,4,‐D in the Rhine river was used to show the evolution from simple to sophisticated models. The first simulation was done with an analytical solution of the dispersion‐advection equation without elimination. Elimination was introduced in a second simulation. The third simulation was carried out with a numerical model. This included elimination and variable dispersion. The lack of data limited the use of very detailed models.     

Quantifying vertical mixing in estuaries

Estuarine turbulence is notable in that both the dissipation rate and the buoyancy frequency extend to much higher values than in other natural environments. The high dissipation rates lead to a distinct inertial subrange in the velocity and scalar spectra, which can be exploited for quantifying the turbulence quantities. However, high buoyancy frequencies lead to small Ozmidov scales, which require high sampling rates and small spatial aperture to resolve the turbulent fluxes. A set of observations in a highly stratified estuary demonstrate the effectiveness of a vessel-mounted turbulence array for resolving turbulent processes, and for relating the turbulence to the forcing by the Reynolds-averaged flow. The observations focus on the ebb, when most of the buoyancy flux occurs. Three stages of mixing are observed: (1) intermittent and localized but intense shear instability during the early ebb; (2) continuous and relatively homogeneous shear-induced mixing during the mid-ebb, and weakly stratified, boundary-layer mixing during the late ebb. The mixing efficiency as quantified by the flux Richardson number Rf was frequently observed to be higher than the canonical value of 0.15 from Osborn (J Phys Oceanogr 10:83–89, 1980). The high efficiency may be linked to the temporal–spatial evolution of shear instabilities.     

Optimum swimming pathways of fish spawning migrations in rivers

Fishes that swim upstream in rivers to spawn must navigate complex fluvial velocity fields to arrive at their ultimate locations. One hypothesis with substantial implications is that fish traverse pathways that minimize their energy expenditure during migration. Here we present the methodological and theoretical developments necessary to test this and similar hypotheses. First, a cost function is derived for upstream migration that relates work done by a fish to swimming drag. The energetic cost scales with the cube of a fish's relative velocity integrated along its path. By normalizing to the energy requirements of holding a position in the slowest waters at the path's origin, a cost function is derived that depends only on the physical environment and not on specifics of individual fish. Then, as an example, we demonstrate the analysis of a migration pathway of a telemetrically tracked pallid sturgeon (Scaphirhynchus albus) in the Missouri River (USA). The actual pathway cost is lower than 10(5) random paths through the surveyed reach and is consistent with the optimization hypothesis. The implication--subject to more extensive validation--is that reproductive success in managed rivers could be increased through manipulation of reservoir releases or channel morphology to increase abundance of lower-cost migration pathways.     

Development of drone-type float for surface-velocity measurement in rivers

Environmental Fluid Mechanics - The present study develops a safety survey system for measuring natural river discharge. Monitoring of rivers is very important for river environment conservation...     

Industrial cooling water disposal in rivers and canals

The geometry of the warm layers of stratified streams is studied on the basis of the one‐dimensional model. Particularly, the derived equations of the dimensionless length and shape of the warm layers are expressed as function of easily measured flow and density parameters, as well as, function of the interfacial friction coefficient, f_i . These analytical models are also expressed in graphic forms, suitable for practical applications. Moreover, when f_i cannot be estimated, approximated forms of the analytical models are also presented graphically.     

Racial mixing in South African honeybees: the effects of genotype mixing on reproductive traits of workers

To test the hypothesis that the honeybee hybrid zone in South Africa is a tension zone due to increased reproductive conflict in colonies that contain both Apis mellifera capensis and Apis mellifera scutellata worker genotypes, we constructed mixed subspecies and hybrid colonies via a combination of artificial and natural matings. We measured emergence weight, ovary activation, and the presence/absence of a spermatheca on workers of different genotypes. We show that the measured characteristics were all affected by genotype with some traits also affected by the social environment in which the worker was reared. Workers with both an A. m. capensis mother and father had the highest emergence weight. When workers had an A. m. capensis mother, paternity affected emergence weight with A. m. capensis fathers producing heavier workers. When the queen was A. m. scutellata, paternity had less effect on weight. Presence of spermatheca was highest in mixed colonies irrespective of maternity and colonies containing pure A. m. capensis workers only. Paternity had a significant effect on the presence of a spermatheca within mixed colonies, with workers that had an A. m. capensis father being more likely to possess a spermatheca. Rates of ovary activation were highest in colonies with an A. m. scutellata queen mated to drones of both genotypes, suggesting that mixed subspecies colonies likely suffer increased reproductive strife among workers. Our results provide support for the hypothesis that the South African honeybee hybrid zone is a tension zone arising from reduced fitness of genetically mixed colonies.     

Scalar mixing in an urban canyon

The scalar dynamics within a unit-aspect-ratio street canyon are studied using large-eddy simulation. The key processes of ventilation and mixing are analysed with the canyon-averaged concentration, mean tracer age and variance. The results are sensitive to the source location and can be classified according to the streamline geometry. The canyon-averaged concentrations for the corner vortices, vortex sea and central vortex do not converge to the same value at large times, though the mean decay rates do. The variance measured with respect to the canyon average shows two distinct decay regimes: the early regime reflects large-scale straining and enhanced diffusion across streamlines, while the late regime is associated with escape from the canyon, i.e., ventilation. Analytical predictions for the variance-decay or mixing time scales are verified for the early regime. It is argued that the presence of an open boundary at the roof level suppresses rapid mixing of the scalar field and is responsible for differences with respect to scalar dynamics within closed domains.

     

Extremely high secondary production of introduced snails in rivers.

The functional importance of invasive animals may be measured as the degree to which they dominate secondary production, relative to native animals. We used this approach to examine dominance of invertebrate secondary production by invasive New Zealand mudsnails (Potamopyrgus antipodarum) in rivers. We measured secondary production of mudsnails and native invertebrates in three rivers in the Greater Yellowstone Area (Wyoming, USA): Gibbon River, Firehole River, and Polecat Creek. Potamopyrgus production was estimated by measuring in situ growth rates and multiplying by monthly biomass; native invertebrate production was estimated using size frequency and instantaneous growth methods. Mudsnail growth rates were high (up to 0.06 d(-1)) for juvenile snails and much lower for adult females (0.003 d(-1)). Potamopyrgus production in Polecat Creek (194 g x m(-2) x yr(-1)) was one of the highest values ever reported for a stream invertebrate. Native invertebrate production ranged from 4.4 to 51 g x m(-2) x yr(-1). Potamopyrgus was the most productive taxon and constituted 65-92% of total invertebrate productivity. Native invertebrate production was low in all streams. Based on a survey of production measures from uninvaded rivers, the distribution of secondary production across taxa was much more highly skewed toward the invasive dominant Potamopyrgus in the three rivers. We suggest that this invasive herbivorous snail is sequestering a large fraction of the carbon available for invertebrate production and altering food web function.     

Preliminary survey of pharmaceutical residues in some important Romanian rivers

A number of important rivers of Romania were investigated for pharmaceutical and antifungal residues. In a post target approach analysis, a single-stage, high-resolution mass spectrometry was successfully applied for the screening of 43 pharmaceuticals and fungicides in water, using U-HPLC-Exactive Orbitrap MS at 50,000 full width at half maximum resolution. Detection was based on accurate masses and retention times. The study confirmed the presence of pharmaceuticals and antifungals such as diclofenac, sulfamethoxazole, carbamazepine, trimethoprim, thiabendazole, and clotrimazole, in water. Diclofenac and carbamazepine were the most frequently detected (eight samples); concentrations ranged from 5 to 50 ng L^?1. Low concentrations of griseofulvin were detected in the Prut River. Further monitoring studies should be carried out in order to fill the gaps of knowledge concerning the presence of pharmaceutical residues in Romanian water environment and to improve public understanding regarding the environmental risk of pharmaceutical contamination.     

Definition and calculation of uncertainty in ecological process models

We present a method of multi-criteria assessment for the analysis of process model uncertainty that combines analysis of model structure, parameters and data requirements. There are three components in calculation and definition of uncertainty.

(1)

Assessment criteria: Uncertainty in a process model is reduced as the model can simultaneously simulate an increased number of assessment criteria selected to test specific aspects of the theory being investigated, and within acceptable limits set for those criteria. This reduces incomplete specification of the model—the characteristic that a model may explain some, but not all, of the observed features of a phenomenon. The calculation required is computation of the Pareto set which provides the list of simultaneously achieved criteria within specified ranges.     

Ecosystem health assessment of urban rivers and lakes

The assessment of the ecosystem health of urban rivers and lakes is the scientific basis for their management and ecological restoration. This study developed a three-level indicator system for its assessment. The results indicated that: Zhonghai and Nanhai are in the state of transition from unhealthy to critical state and all the other lakes are in unhealthy states. Water environmental quality, structure and function of the aquatic ecosystem, and the structure of waterfront areas were the constraints. Nanhai was ranked as poor and the others were all ranked as very poor. However, the ecological environment of Zhonghai and Nanhai were better than the others, the sums of the degree of membership to the healthy state and critical state were all close to 0.6. and the restorations of these lakes were moderate. The sums of the degree of membership to the healthy state and critical state of the other lakes were under 0.3, as it was difficult to restore these lakes. Some suggestions on scientific management and ecological restoration of the six lakes were proposed: $`To control non-point pollution and to improve the water quality of six lakes and the water entering into these lakes; $aTo improve the hydrological conditions of six lakes; $bTo rehabilitate the aquatic ecosystem and waterfront areas.     

The degradation of plastic litter in rivers: Implications for beaches

Polythene sheeting is a major litter component on estuarine beaches and river banks. Sanitary towel backing strips, which are one of the commonest items of sewage related debris found on beaches, enter the riverine system via combined sewer overflows. Investigations on these items, positioned at natural riverine stranding levels, showed that after an initial rapid breakdown little further loss of tensile strength occurred. Experiments carried out on backing strips, buried in the bank suspended from a tree and tethered to the bank, showed significant change. Buried samples retained the greatest tensile strength retention, dropping no lower than 90 %, the other samples showed similar retention rates at 80%. The difference is probably due to photodegradation as biodegradation effects were minimal. Probably, the longevity of such plastics is a major reason for their abundance and widespread distribution both on river banks and beaches.     

Deterministic particle tracking simulation of pollutant discharges in rivers and estuaries

A two-dimensional deterministic particle tracking model, in which the anisotropic-dispersive process is described by a particle strength exchange scheme, is established for the simulation of pollutant transport in vertically well-mixed rivers and estuaries. By simulating two benchmark problems with analytic solutions, the PSE scheme is shown to be accurate even if the anisotropic ratio of dispersion coefficients is very high. Further simulations of two specific problems concerning the optimal effluent discharge location and procedure are presented. The major conclusion is that in a tidal estuary with a relatively large fresh-water flow, setting the discharge position at the mixing center and making the discharge rate proportional to flow speed may minimize the peaks of concentration.     

渔业养殖区水环境中汞的形态分布

为研究渔业养殖对水生生态环境中汞形态变化的影响,采用气相色谱-冷原子荧光法对浙江省舟山市六横岛海水样品中不同形态汞含量进行测定分析.结果表明,对于总汞(THg)、溶解态汞(DHg)和颗粒态汞(PHg),养殖区上浮水体与对照区上浮水体之间无显著差异(P0.05).纵向分布看,底层水体PHg含量要显著高于中间层和表层(P0.05),主要是由于沉积物向上浮水体释放大量颗粒物,而颗粒物本身对汞的吸附能力很强.对于上浮水体甲基汞含量分布,总甲基汞(TMe Hg)主要以(PMe Hg)的形态存在,且TMe Hg和PMe Hg在颗粒物含量较高的养殖区的含量要高于在对照区.但是,底层水体中可被微生物利用的溶解态甲基汞(DMe Hg)含量在对照区(0.066±0.009 pmol·L-1)要显著高于养殖区(0.031±0.007 pmol·L-1)(P0.05),这可能是由于常年养殖活动致使沉积物中大量有机物积累,有机物与无机汞离子络合,减少了汞的生物有效性,从而抑制了甲基汞的生成.     

城市污染河道沉积物有机质的氧化稳定性

为了研究城市河道沉积物有机质的氧化稳定性,以氧化稳定系数(Kos)和腐殖化程度(HA/FA)作为评价指标,在用H2O2处理沉积物的基础上,探讨城市污染河道有机质含量、有机质组成等对氧化稳定性的影响.结果显示,污染河道沉积物有机质含量明显高于自然水生态系统沉积物有机质含量;富里酸(FA)、胡敏酸(HA)、胡敏素(HM)是污染河道沉积物有机质主要赋存形式,其中胡敏素是主体,占全部有机质的65%以上;易氧化有机质(OMr)和难氧化有机质(OMd)与总有机质含量均显著正相关(P0.05);OMd的含量可以更直观地反映河道污染程度;重污染河道沉积物有机质去除量较高;随着有机质的去除,沉积物Kos在总体上有上升的趋势;不同形态有机组分其Kos差异性明显(P0.05),以HM为主的紧结态腐殖质的Kos值最高;用HA/FA和OMr可从大体上表征沉积物的氧化稳定性;HA、FA是影响城市河道沉积物氧化稳定性的主要因素.     

Integrating wetlands and riparian zones in river basin modelling

Wetlands, and in particular riparian wetlands, represent an interface between the catchment area and the aquatic environment. They control the exchange of water and related chemical fluxes from the upper catchment area to surface waters like streams and lakes. Their influence on water and nutrient balances has been investigated mainly at the patch scale. In this study an attempt was made (a) to integrate riparian zones and wetlands into eco-hydrological river basin modelling, and (b) to quantify the impacts of riparian wetland processes on water and nutrient fluxes in a meso-scale catchment located in the northeastern German lowland. The investigation was performed by analysing hydro-chemical field data and applying the eco-hydrological model SWIM (Soil and Water Integrated Model), which was extended to reproduce the relevant water and nutrient flows and retention processes at the catchment scale in general, and in riparian zones and wetlands in particular. The main extensions introduced in the model were: (1) implementation of daily groundwater table dynamics at the hydrotope level, (2) implementation of water and nutrient uptake by plants from groundwater in riparian zones and wetlands, and (3) assessment of nutrient retention in groundwater and interflow. The simulation results indicate that wetlands, though they represent relatively small parts of the total catchment area, may have a significant impact on the overall water and nutrient balances of the catchment. The uncertainty of the simulation results is considerably high, with the main sources of uncertainty being the model parameters representing the geo-hydrology and the input data for land use management.     

Initial mixing of inclined dense jet in perpendicular crossflow

A comprehensive experimental investigation for an inclined ( $60^{\circ }$ to vertical) dense jet in perpendicular crossflow—with a three-dimensional trajectory—is reported. The detailed tracer concentration field in the vertical cross-section of the bent-over jet is measured by the laser-induced fluorescence technique for a wide range of jet densimetric Froude number $Fr$ and ambient to jet velocity ratios $U_r$ . The jet trajectory and dilution determined from a large number of cross-sectional scalar fields are interpreted by the Lagrangian model over the entire range of jet-dominated to crossflow-dominated regimes. The mixing during the ascent phase of the dense jet resembles that of an advected jet or line puff and changes to a negatively buoyant thermal on descent. It is found that the mixing behavior is governed by a crossflow Froude number $\mathbf{F} = U_r Fr$ . For $\mathbf{F} < 0.8$ , the mixing is jet-dominated and governed by shear entrainment; significant detrainment occurs and the maximum height of rise $Z_{max}$ is under-predicted as in the case of a dense jet in stagnant fluid. While the jet trajectory in the horizontal momentum plane is well-predicted, the measurements indicate a greater rise and slower descent. For $\mathbf{F} \ge 0.8$ the dense jet becomes significantly bent-over during its ascent phase; the jet mixing is dominated by vortex entrainment. For $\mathbf{F} \ge 2$ , the detrainment ceases to have any effect on the jet behavior. The jet trajectory in both the horizontal momentum and buoyancy planes are well predicted by the model. Despite the under-prediction of terminal rise, the jet dilution at a large number of cross-sections covering the ascent and descent of the dense jet are well-predicted. Both the terminal rise and the initial dilution for the inclined jet in perpendicular crossflow are smaller than those of a corresponding vertical jet. Both the maximum terminal rise $Z_{max}$ and horizontal lateral penetration $Y_{max}$ follow a $\mathbf{F}^{-1/2}$ dependence in the crossflow-dominated regime. The initial dilution at terminal rise follows a $S \sim \mathbf{F}^{1/3}$ dependence.