Size, not morphology, determines hydrodynamic performance of a kelp during peak flow

The morphology and shape of algae can affect their survival in wave-swept environments because of the hydrodynamic drag created by water flow. Studies of morphology and drag are typically conducted at relatively low water velocities, and the influence of algal morphology on drag, over the range of water velocities algae must cope with in their natural environment, remains unclear. Here, we tested the link between morphological variation and hydrodynamic drag for a dominant kelp with complex morphology (Ecklonia radiata), over a range of water velocities representative of conditions on wave-swept reefs. Our results indicated that kelps on subtidal reefs must withstand maximal orbital water velocities in excess of 2–3 m s^?1. Our measurements of drag, resulting from flows ranging from 1 to 3 m s^?1, revealed that shape- and width-related thallus and lamina characters were important to drag at low speed, but that total thallus area (or biomass) was the main determinant of drag at high flow. Drag coefficients converged at increasing speed suggesting that, at high flow, significant thallus reconfiguration (more streamlined shape) decoupled drag from morphology. This implies that, at peak velocities, only size (total area), not morphology, is important to drag and the probability of dislodgment.     

The nitrate time bomb: a numerical way to investigate nitrate storage and lag time in the unsaturated zone

Nitrate pollution in groundwater, which is mainly from agricultural activities, remains an international problem. It threatens the environment, economics and human health. There is a rising trend in nitrate concentrations in many UK groundwater bodies. Research has shown it can take decades for leached nitrate from the soil to discharge into groundwater and surface water due to the ‘store’ of nitrate and its potentially long travel time in the unsaturated and saturated zones. However, this time lag is rarely considered in current water nitrate management and policy development. The aim of this study was to develop a catchment-scale integrated numerical method to investigate the nitrate lag time in the groundwater system, and the Eden Valley, UK, was selected as a case study area. The method involves three models, namely the nitrate time bomb—a process-based model to simulate the nitrate transport in the unsaturated zone (USZ), GISGroundwater—a GISGroundwater flow model, and N-FM—a model to simulate the nitrate transport in the saturated zone. This study answers the scientific questions of when the nitrate currently in the groundwater was loaded into the unsaturated zones and eventually reached the water table; is the rising groundwater nitrate concentration in the study area caused by historic nitrate load; what caused the uneven distribution of groundwater nitrate concentration in the study area; and whether the historic peak nitrate loading has reached the water table in the area. The groundwater nitrate in the area was mainly from the 1980s to 2000s, whilst the groundwater nitrate in most of the source protection zones leached into the system during 1940s–1970s; the large and spatially variable thickness of the USZ is one of the major reasons for unevenly distributed groundwater nitrate concentrations in the study area; the peak nitrate loading around 1983 has affected most of the study area. For areas around the Bowscar, Beacon Edge, Low Plains, Nord Vue, Dale Springs, Gamblesby, Bankwood Springs, and Cliburn, the peak nitrate loading will arrive at the water table in the next 34 years; statistical analysis shows that 8.7 % of the Penrith Sandstone and 7.3 % of the St Bees Sandstone have not been affected by peak nitrate. This research can improve the scientific understanding of nitrate processes in the groundwater system and support the effective management of groundwater nitrate pollution for the study area. With a limited number of parameters, the method and models developed in this study are readily transferable to other areas.     

Alternating skimming flow over a stepped spillway

The study of stepped spillways in laboratory scales has been essentially focused on two separated sub-regimes within skimming flow. In this paper we investigate the appearance of an unclassified alternating skimming flow regime in a 0.5 m wide stepped spillway which does not fit on these earlier definitions, and which does not occur in a 0.3 m wide spillway. Our aim is to explain the genesis of this unclassified flow which is visualised in the physical stepped spillway, by using 3D numerical modelling. Flow depths and velocities are measured using an ultrasonic sensor and Bubble Image Velocimetry in the wider flume (0.5 m). The numerical model is validated with the experimental data from the 0.5 m wide spillway. After validation, the channel width of the same numerical model is reduced to 0.3 m wide spillway in order to characterise (compare) the case without (with) alternating skimming flow. Both cases are solved using Reynolds-Averaged Navier–Stokes equations together with the Volume-of-Fluid technique and SST k-\(\omega\) turbulence model. The experimental results reveal that the alternating skimming flow regime is characterised by an evident seesaw pattern of flow properties over consecutive steps. In turn, the numerical modelling clarified that this seesaw pattern is due to the presence of a complex system of cross waves along the spillway. These cross waves are also responsible for a mass and momentum exchange in the transversal direction and for the formation of the alternating skimming flow in the spillway.     

有机污染物在土壤和地下水中迁移建模

地下水有机污染是由人类活动引起的各种物理、化学和生物等干扰过程造成有机污染物自土壤表面迁移至地下含水层的结果。为了解有机污染物如何自土表经过不饱和层（包气带）进入含水层，需要对不饱和层中所进行的各种物理、化学和生物等过程进行较准确的定量描述。本文就不饱和层和含水层中水分和有机污染物迁移建模的一些基本概念和方法进行综述，并列举了一些经典的和新颖的建模方法。     

A simulation model for projecting changes in salinity concentrations and species dominance in the coastal margin habitats of the Everglades

Sharp boundaries typically separate the salinity tolerant mangroves from the salinity intolerant hardwood hammock species, which occupy the similar geographical areas of southern Florida. Evidence of strong feedback between tree community-type and the salinity of the unsaturated (vadose) zone of the soil suggests that a severe disturbance that significantly tilts the salinity in the vadose zone might cause a shift from one vegetation type to the other. In this study, a model based upon the feedback dynamics between vegetation and salinity of the vadose zone of the soil was used to take account of storm surge events to investigate the mechanisms that by which this large-scale disturbance could affect the spatial pattern of hardwood hammocks and mangroves. Model simulation results indicated that a heavy storm surge that completely saturated the vadose zone at 30 ppt for 1 day could lead to a regime shift in which there is domination by mangroves of areas previously dominated by hardwood hammocks. Lighter storm surges that saturated the vadose zone at less than 7 ppt did not cause vegetation shifts. Investigations of model sensitivity analysis indicated that the thickness of the vadose zone, coupled with precipitation, influenced the residence time of high salinity in the vadose zone and therefore determined the rate of mangrove domination. The model was developed for a southern Florida coastal ecosystem, but its applicability may be much broader.     

Quantification of advective solute travel times and mass transport through hydrological catchments

This study has investigated and outlined the possible quantification and mapping of the distributions of advective solute travel times through hydrological catchments. These distributions are essential for understanding how local water flow and solute transport and attenuation processes affect the catchment-scale transport of solute, for instance with regard to biogeochemical cycling, contamination persistence and water quality. The spatial and statistical distributions of advective travel times have been quantified based on reported hydrological flow and mass-transport modeling results for two coastal Swedish catchments. The results show that the combined travel time distributions for the groundwater-stream network continuum in these catchments depend largely on the groundwater system and model representation, in particular regarding the spatial variability of groundwater hydraulic parameters (conductivity, porosity and gradient), and the possible contributions of slower/deeper groundwater flow components. Model assumptions about the spatial variability of groundwater hydraulic properties can thus greatly affect model results of catchment-scale solute spreading. The importance of advective travel time variability for the total mass delivery of naturally attenuated solute (tracer, nutrient, pollutant) from a catchment to its downstream water recipient depends on the product of catchment-average physical travel time and attenuation rate.     

Estimating river flow effects on water ages by hydrodynamic modeling in Little Manatee River estuary, Florida, USA

The water age in a tidal river in Florida, Little Manatee River, has been investigated in this study by the application of a three-dimensional hydrodynamic model. In response to a pulse dye release in the upper end of the river boundary, the hydrodynamic model determines the water age for a given location by recording the time for the dye to reach the given river location. The hydrodynamic model uses horizontal curvilinear orthogonal grids to represent the complex river system that includes several bayous and tributaries. The model has was calibrated and verified in previous study by using two continuous data sets for a 6 month period. Satisfactory model verifications indicate that the hydrodynamic model is capable of quantifying the mixing and transport process for calculating the water age in the tidal river. For 17 freshwater inflow scenarios in the Little Manatee River, the hydrodynamic model was applied to simulate water ages along the main channel of the river at 2-km interval. Flow rates in the 17 scenarios varying from 0.26 to 76.56 m³/s cover the range of the observed flows in the Little Manatee River. Water ages from model predictions range from the minimum 1.2 days under the maximum 76.56 m³/s inflow condition to the 50 days under the minimum 0.26 m³/s inflow condition. Empirical regression equations at three selected stations, with the correlation coefficient R² above 0.96, were derived from numerical model simulations to correlate water ages to freshwater inflows. The empirical water-age equation derived from hydrodynamic model simulations can be used to provide quick and low-cost estimations of water ages in response to various inflow scenarios for studying physical–chemical and biological processes in the river.     

The similarity solution for turbulent mixing of two-layer stratified fluid

Experiments are reviewed in which a two-layer salt-stratified tank of water was mixed by turbulence. The density profile began as a single step and evolved to a smooth mixed profile. The turbulence was generated by many excursions of a horizontally moving vertical rod with Richardson number Ri >  0.9 and Reynolds Number Re  >  600. There was almost perfect collapse of all the profiles to one universal profile as a function of a similarity variable. We develop a theoretical model for a simple mixing law with a buoyancy flux that is a function of internal Richardson number Rii. A similarity equation is found. A flux law that increases with small Rii and decreases with large Rii is considered next. Since no analytical solution is known, the similarity concept is tested by numerically integrating the equations in space and time. With buoyancy flux monotonically increasing with internal Richardson number, the similarity approach is valid for a profile starting from a slightly smoothed step. However, a shock forms for a mixing law with higher initial Rii (so that buoyancy flux decreases with Richardson number) and the similarity approach is invalid for those initial conditions.     

Experimental assessment of characteristic turbulent scales in two-phase flow of hydraulic jump: from bottom to free surface

A hydraulic jump is a turbulent shear flow with a free-surface roller. The turbulent flow pattern is characterised by the development of instantaneous three-dimensional turbulent structures throughout the air–water column up to the free surface. The length and time scales of the turbulent structures are key information to describe the turbulent processes, which is of significant importance for the improvement of numerical models and physical measurement techniques. However, few physical data are available so far due to the complexity of the measurement. This paper presents an investigation of a series of characteristic turbulent scales for hydraulic jumps, covering the length and time scales of turbulent flow structures in bubbly flow, on free surface and at the impingement point. The bubbly-flow turbulent scales are obtained for Fr = 7.5 with 3.4 × 10⁴ < Re < 1.4 × 10⁵ in both longitudinal and transverse directions, and are compared with the free-surface scales. The results highlight three-dimensional flow patterns with anisotropic turbulence field. The turbulent structures are observed with different length and time scales respectively in the shear flow region and free-surface recirculation region. The bubbly structures next to the roller surface and the free-surface fluctuation structures show comparable length and time scales, both larger than the scales of vortical structures in the shear flow and smaller than the scales of impingement perimeter at the jump toe. A decomposition of physical signals indicates that the large turbulent scales are related to the unsteady motion of the flow in the upper part of the roller, while the high-frequency velocity turbulence dominates in the lower part of the roller. Scale effects cannot be ignored for Reynolds number smaller than 4 × 10⁴, mainly linked to the formation of large eddies in the shear layer. The present study provides a comprehensive assessment of turbulent scales in hydraulic jump, including the analyses of first data set of longitudinal bubbly-flow integral scales and transverse jump toe perimeter integral scales.     

Modeling the buoyant flow of heated water discharged from surface and submerged side outfalls in shallow and deep water with a cross flow

In this study, a three-dimensional model was used to numerically study the buoyant flow, along with its mixing characteristics, of heated water discharged from the surface and submerged side outfalls in shallow and deep water with a cross flow. Hydraulic experimental data were used to evaluate the applicability of the model. The simulation results agree well with the experimental results, particularly, the jet trajectories, the dimensions of the recirculating zone, and the distribution of the dimensionless excess temperature. The level of accuracy of the simulation results of the present study is nearly identical to that of the results conducted by McGuirk and Rodi (1978). If the heated water is discharged into shallow water where the momentum flux ratio and the discharge densimetric Froude number are high, the submerged discharge method is better than the surface discharge method in terms of the scale of the recirculating zone and the minimum dilution. In deep water, where the momentum flux ratio and discharge densimetric Froude number are low, however, the submerged discharge method had few advantages. In shallow water, the discharge jet is deflected by the ambient cross flow, while forcing the ambient flow to bend towards the far bank for the full depth. For a submerged discharge in shallow water, the recirculating zone is the largest in the lowest layer but becomes smaller in the upper layer. As the water depth increases, the ambient flow goes over the jet and diminishes the blocking effect, thereby decreasing the bending of the jet.     

Temporal variability of deep-sea zooplankton in the Arabian Sea

Mesozooplankton samples from two stations in the Arabian Sea (WAST, 4,050 m, 16°15′N, 60°20′E; CAST, 3,950 m, 14°30′N, 64°30′E) were collected from the surface down to 20 m above bottom during three monsoon periods: the autumn intermonsoon in October 1995, the spring intermonsoon in April 1997, and the NE monsoon in February 1998. The main goal of this study is to enhance our knowledge on the effect of spatial and temporal differences in primary production and particle flux rates on the abundance and distribution of mesozooplankton, with special attention to the deep sea. Literature data indicate episodically high rates of primary production and particle flux in the region during the SW monsoon and the autumn intermonsoon. Set in this context, the zooplankton showed an in-phase coupling in biomass and abundance with the primary production in the surface 150 m. In the mesopelagic realm (150–1,050 m), the seasonal coupling was less clear. In the bathypelagic zone, below 1,050 m, there was no evidence of in-phase coupling, though temporal differences in the distribution of zooplankton abundance and biomass with depth between seasons could be shown by an analysis of covariance and an a posteriori test. The results suggest that the bathypelagic community responds to increased particle flux rates, but with longer time gaps than in the epipelagic zone. This is probably due to longer development and response times of zooplankton in the cold, deep-water environment, independent of possible lateral advection processes.     

Potential environmental impact of dissolution of raw phosphate in sea water of the Gulf of Aqaba

Dissolution of raw phosphate (apatite) in sea water of the Gulf of Aqaba was investigated through lab incubation experiments. Three types from three different sources (Al-Hasa, Al-Abyad and Esh-Shydiya) have been used for these experiments. Impact of quantity, grain size, and source (type) of raw phosphate on dissolution rate were studied. Statistical analysis shows significant differences between the results obtained from comparing each two weights; as weight of apatite increased, dissolved inorganic phosphate-phosphorus (DIP) and fluoride in sea water solution increased. The differences between the dissolution rates of raw phosphate from the three sites were not significant while the differences between the different grain size fractions were significant. Dissolution rates were inversely related to particle size. Using a worst-case scenario, a conservative estimate of the maximum increase in DIP in seawater of the Gulf of Aqaba due to the apatite particles lost to the sea during ship loading resulted in DIP concentrations of 0.03 μM per year. As the residence time of the water in the Gulf of Aqaba is about one year, the DIP concentration will not increase by more than 0.03 μM under the estimated annual quantity of exported phosphate. Fluoride will not increase by more than 0.03 mg/l under the same conditions.     

Long-term behavior at foraging sites of adult female loggerhead sea turtles (Caretta caretta) from three Florida rookeries

We used satellite telemetry to study behavior at foraging sites of 40 adult female loggerhead sea turtles (Caretta caretta) from three Florida (USA) rookeries. Foraging sites were located in four countries (USA, Mexico, the Bahamas, and Cuba). We were able to determine home range for 32 of the loggerheads. One turtle moved through several temporary residence areas, but the rest had a primary residence area in which they spent all or most of their time (usually >11 months per year). Twenty-four had a primary residence area that was <500 km² (mean = 191). Seven had a primary residence area that was ≥500 km² (range = 573–1,907). Primary residence areas were mostly restricted to depths <100 m. Loggerheads appeared to favor areas with larger-grained sediment (gravel and rock) over areas with smaller-grained sediment (mud). Short-term departures from primary residence areas were either looping excursions, typically involving 1–2 weeks of continuous travel, or movement to a secondary residence area where turtles spent 25–45 days before returning to their primary residence area. Ten turtles had a secondary residence area, and six used it as an overwintering site. For those six turtles, the primary residence area was in shallow water (<17 m) in the northern half of the Gulf of Mexico (GOM), and overwintering sites were farther offshore or farther south. We documented long winter dive times (>4 h) for the first time in the GOM. Characterizing behaviors at foraging sites helps inform and assess loggerhead recovery efforts.     

Seasonal and decadal changes in distribution patterns of Halobates (Hemiptera: Gerridae) populations in the eastern tropical Pacific

Five species of the marine insect Halobates share similar ecology but have distinct biogeographic ranges in the eastern tropical Pacific, a region from approximately 75°W–160°W and 10°S–35°N. Between 2001 and 2010, the Sea Education Association collected Halobates from 682 neuston tows (surface net 1 m × 0.5 m, 335-μm mesh) during fifteen cruises between San Diego, USA, Mexico and Tahiti. Total Halobates spp. densities varied substantially from year to year, but our data do not show a sustained change from a data set collected 40 years earlier from 1967 to 1968 (Cheng and Shulenberger in Fish Bull 78(3):579–591, 1980). Halobates are sensitive to sea surface temperature and we observed significant differences in species distributions over time, but these were not due to differences in water temperature or climate change. Our analyses show that the patterns observed are attributable to substantial but previously undescribed seasonal shifts that occur each year in the ranges for both Halobates sobrinus and Halobates micans. There is substantial overlap in ranges during seasonal shifts, but very little co-occurrence of H. sobrinus and H. micans in individual net tows, suggesting biological mechanisms rather than physical factors are restricting distribution and co-occurrence of these two species.     

Mass transfer limitation of photosynthesis of coral reef algal turfs

Algal turfs are the major primary producing component on many coral reefs and this production supports higher levels in the complex reef trophic web. Rates of metabolism of algal turfs are related positively to water motion, consistent with limitation by the diffusion of a substance through a boundary layer. Based on engineering mass transfer theory, we hypothesized that photosynthesis of algal turfs is controlled by rates of mass transfer and responses of photosynthesis to increasing flow speed should be predicted by engineering correlations. This hypothesis was tested in ten experiments where photosynthesis was estimated in a flume/respirometer from changes in dissolved oxygen at eight flow speeds between 0.08 and 0.52 m/s. Flow in the flume and over the reef at Kaneohe Bay, Oahu, Hawaii was estimated using hot-film thermistor and electromagnetic current meters. Rates of photosynthesis were related positively to flow in all experiments and plots of the log of the average Sherwood number (Sh _meas) versus log Reynolds number (Re _D) for each experiment are lower than predicted for mass transfer through a turbulent boundary layer. Algal turf-covered plates are characterized as hydrodynamically transitional to fully rough surfaces and the lower than predicted slopes suggest that roughness reduces rates of mass transfer. A negative correlation between algal turf biomass and slopes of the log Sh _meas−log Re _D plots suggests that mass transfer to algal turfs is affected significantly by the physical structure of the algal community. Patterns of photosynthesis based on changes in dissolved oxygen and dissolved inorganic carbon concentrations (DIC) indicate that the flow speed effect is not the result of increased flux of oxygen from the algal turfs, and combined with the short response time to flow speed, suggest that DIC may limit rates of photosynthesis. Although there are differences between flow in the flume and flow over algal turfs on the reef, these results suggest that photosynthesis is controlled, at least in part, by mass transfer. The chemical engineering approach provides a framework to pose further testable hypotheses about how algal canopy height, flow oscillation, turbulence, and substratum roughness may modulate rates of metabolism of coral reef algal turfs.     

Sediment processes and flow reversal in the undular tidal bore of the Garonne River (France)

A tidal bore is a series of waves propagating upstream as the tidal flow turns to rising, and the bore front corresponds to the leading edge of the tidal wave in a funnel shaped estuarine zone with macro-tidal conditions. Some field observations were conducted in the tidal bore of the Garonne River on 7 June 2012 in the Arcins channel, a few weeks after a major flood. The tidal bore was a flat undular bore with a Froude number close to unity: $\hbox {Fr}_{1} = 1.02$ and 1.19 (morning and afternoon respectively). A key feature of the study was the simultaneous recording of the water elevation, instantaneous velocity components and suspended sediment concentration (SSC) estimates, together with a detailed characterisation of the sediment bed materials. The sediment was some silty material ( $\hbox {d}_{50} \approx 13~\upmu \hbox {m}$ ) which exhibited some non-Newtonion thixotropic behaviour. The velocity and SSC estimate were recorded simultaneously at high frequency, enabling a quantitative estimate of the suspended sediment flux at the end of the ebb tide and during the early flood tide. The net sediment flux per unit area was directed upstream after the bore, and its magnitude was much larger than that at end of ebb tide. The field observations highlighted a number of unusual features on the morning of 7 June 2012. These included (a) a slight rise in water elevation starting about 70 s prior to the front, (b) a delayed flow reversal about 50 s after the bore front, (c) some large fluctuations in suspended sediment concentration (SSC) about 100 s after the bore front and (d) a transient water elevation lowering about 10 min after the bore front passage. The measurements of water temperature and salinity showed nearly identical results before and after the tidal bore, with no evidence of saline and thermal front during the study.     

Longshore sediment transport in the surf zone based on different formulae: a case study along the central west coast of India

Understanding longshore sediment transport (LST) is a prerequisite for designing an effective coastal zone management strategy. The present study estimates the LST along the central west coast of India based on four bulk LST formulae: (1) the Coastal Engineering Research Center (CERC) formula, (2) the Walton and Bruno formula, (3) the Kamphuis formula and (4) the Komar formula. The Delft3D–wave module is used for estimating nearshore waves with measured directional wave data at a water depth of 9 m as the input parameter. Wave data for the validation of the nearshore wave transformation model is measured using the InterOcean S4DW wave gauge. The model results show that waves approach from the south 90 % of the time in a year and that they generate predominantly northerly longshore currents. Upon comparison with the measured data, the findings show that the estimates based on the Kamphuis formula agree with the field data. A high ratio (~1) of the monthly net and gross transport rates indicates that the LST is dominating in one direction in all months except February and July. The study shows that a slight change in the angle of the wave approach during the Asian summer monsoon period (JJAS) can significantly alter the direction and magnitude of the LST. Inter-annual variations in the LST based on the data for 2009 and 2011 show that the variations in the annual net and gross LST rates in different years are less than 7 %.     

Complete desalination of seawater using a novel polyvinylidene fluoride/zeolite membrane

Water scarcity arises rapidly with increasing human population, urbanization and industrialization. Therefore, water production from natural sources using efficient and environmental friendly technologies is gaining attention. Pervaporative desalination is a promising technology based upon membrane filtration. Here, we fabricated and tested a novel hybrid-composite membrane: polyvinylidene fluoride (PVDF)/polyvinyl pyrrolidone/(PVP)/3A zeolite. Desalination performance was evaluated in terms of flux, salt rejection, zeolite incorporation within the membrane and temperature. Results show that all membranes exhibited salt rejection higher than 99.5%. The highest flux of 2.5 kg m^?2 h^?1 and salt rejection of 99.9% were obtained using 10 wt% 3A incorporated PVDF/PVP membrane at 40 °C. The PVDF/PVP/3A zeolite membrane preserved its stability during 55-h desalination operating without performance decline.     

Influence of varying nutrient and pesticide mixtures on abatement efficiency using a vegetated free water surface constructed wetland mesocosm

The nutrient and pesticide abatement efficiency of varying mixtures was examined in a vegetated free water surface constructed wetland. Three different agricultural chemical pollutant mixture conditions were assessed: nutrients only (N and P); pesticides only (atrazine, S-metolachlor and permethrin); and a mixture of nutrients and pesticides. With nutrients only, 672 h nutrient mitigation of 77–91% total phosphorous (TP) and 74–98% total nitrogen (TN) was associated with distance from the injection point and rainfall, whereas with nutrient and pesticide mixtures, 672 h nutrient mitigation of 11–71% TP and 84–98% TN were associated with distance and time. With pesticides only, 672 h pesticide mitigation of 50–99% was associated with distance and time, whereas with nutrients and pesticide mixtures, 672 h pesticide mitigation of 48–99% was associated primarily with distance. Dissipation half-lives were 2–10 times greater for P and 1.5–5 times greater for N when pesticides were present. Pesticide dissipation half-lives showed no clear differences with or without nutrients. While vegetated free water surface constructed wetlands can be effective best management practice tools to trap and abate agricultural run-off during rainfall events, efficiencies can be affected by different types of complex pollutant mixtures and wetland design and implementation should accommodate varying efficiencies.     

漓江上游纳污最适总量及水质变化规律研究

为确保 2 0 0 0 - 2 0 10年桂林漓江上游进入城区水质能保持国家II类水质标准进行纳污总量控制及其水质变化规律研究。采用系数法对 2 0 0 0 - 2 0 10年陆源排污总量进行科学预测 ,按河流功能区划要求 ,首次利用动态数学模型计算河流在不同流量、不同纳污总量条件下保持II类水质的最适纳污值 ,探讨了漓江水质变化规律。为环境管理打下基础