Experimental investigation of bubbly flow and turbulence in hydraulic jumps

Many environmental problems are linked to multiphase flows encompassing ecological issues, chemical processes and mixing or diffusion, with applications in different engineering fields. The transition from a supercritical flow to a subcritical motion constitutes a hydraulic jump. This flow regime is characterised by strong interactions between turbulence, free surface and air–water mixing. Although a hydraulic jump contributes to some dissipation of the flow kinetic energy, it is also associated with increases of turbulent shear stresses and the development of turbulent eddies with implications in terms of scour, erosion and sediment transport. Despite a number of experimental, theoretical and numerical studies, there is a lack of knowledge concerning the physical mechanisms involved in the diffusion and air–water mixing processes within hydraulic jumps, as well as on the interaction between the free-surface and turbulence. New experimental investigations were undertaken in hydraulic jumps with Froude numbers up to Fr = 8.3. Two-phase flow measurements were performed with phase-detection conductivity probes. Basic results related to the distributions of void fraction, bubble frequency and mean bubble chord length are presented. New developments are discussed for the interfacial bubble velocities and their fluctuations, characterizing the turbulence level and integral time scales of turbulence representing a “lifetime” of the longitudinal bubbly flow structures. The analyses show good agreement with previous studies in terms of the vertical profiles of void fraction, bubble frequency and mean bubble chord length. The dimensionless distributions of interfacial velocities compared favourably with wall-jet equations. Measurements showed high turbulence levels. Turbulence time scales were found to be dependent on the distance downstream of the toe as well as on the distance to the bottom showing the importance of the lower (channel bed) and upper (free surface) boundary conditions on the turbulence structure.     

Spatial–temporal structure of mixing interface turbulence at two large river confluences

Converging flows at stream confluences often produce highly turbulent conditions. The shear layer/mixing interface that develops within the confluence hydrodynamic zone (CHZ) is characterized by complex patterns of three-dimensional flow that vary both spatially and temporally. Previous research has examined in detail characteristics of mean flow and turbulence along mixing interfaces at small stream confluences and laboratory junctions; however few, if any, studies have examined these characteristics within mixing interfaces at large river confluences. This study investigates the structure of mean velocity profiles as well as spatial and temporal variations in velocity, backscatter intensity, and temperature within the mixing interfaces of two large river confluences. Velocity, temperature, and backscatter intensity data were obtained at stationary locations within the mixing interfaces and at several cross sections within the CHZ using acoustic Doppler current profilers. Results show that mean flow within the mixing interfaces accelerates over distance from the junction apex. Turbulent kinetic energy initially increases rapidly over distance, but the rate of increase diminishes downstream. Hilbert–Huang transform analysis of time series data at the stationary locations shows that multiple dominant modes of fluctuations exist within the original signals of velocity, backscatter intensity, and temperature. Frequencies of the largest dominant modes correspond well with predicted frequencies for shallow wake flows, suggesting that mixing-interface dynamics include wake vortex shedding—a finding consistent with spatial patterns of depth-averaged velocities at measured cross sections. Spatial patterns of temperature and backscatter intensity show that the converging flows at both confluences do not mix substantially, indicating that turbulent structures within the mixing interfaces are relatively ineffective at producing mixing of the flows in the CHZ.     

Flows through forest canopies in complex terrain

Recent progress on boundary layer flow within and above tall forest canopies in complex terrain is reviewed from the perspective of developing methods to interpret carbon dioxide fluxes from tower measurements in real terrain. Two examples of complex terrain are considered in detail: a forest edge, which exemplifies nonuniform forests, and hilly terrain, which can lead to drainage currents at night. Dynamical arguments show that, when boundary layer winds approach a forest edge, the mean wind adjusts on a length scale of approximately 3L(c), where L(c) is the canopy drag length scale, which depends inversely on the leaf area density of the forest. Over a further distance that also scales on L(c), turbulence in the flow adjusts, and the mixing and transport in the canopy approaches the homogeneous limit. Even low hills change the neutral flow within and above the forest canopy substantially. When the canopy is tall, pressure gradients drive flow up both the upwind and downwind slopes of the hill, leading to an ejection of air out of the top of the canopy just downwind of the crest. This flow at the crest can then advect scalar out of the top of the forest, leading to large variations in the flux of scalar across the hill. At night, when the air near the ground cools and becomes stably stratified, turbulence within the canopy can collapse, even when the flow above the canopy remains turbulent. This leads to a decoupling of the air motions within the canopy from those above. The air above the canopy can then continue to pass up and over the hill, as it does in the neutral case, but at the same time, air within the canopy drains down the hill slopes as drainage currents. These analyses will help us understand when flux towers are reliably measuring the net ecosystem exchange and suggest ways of correcting the flux tower data in more complex situations.     

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.

     

Wave mixing rise inferred from Lyapunov exponents

We study the horizontal surface mixing and the transport induced by waves in a coastal environment. A comparative study is addressed by computing the Lagrangian coherent structures, via Finite Size Lyapunov Exponents, that arise in two different numerical settings: with and without wave coupled to currents. In general, we observe that mixing is increased in the area due to waves. Besides, the methodology presented here is tested by deploying a set of eight Lagrangian drifters at different locations. This dynamical approach is shown as a valuable tool to extract information about transport, mixing and residence embedded in the Eulerian time dependent velocity fields obtained from numerical models.     

Adsorption distribution impact on preferential transport within horizontal flow constructed wetland (HFCW)

Predicting preferential water flow and solute transport in soils is in the interest of scientists and engineers in the fields of agricultural soil, forest hydrology, soil physics and wastewater treatment by constructed wetland. In artificial wetlands, any preferential pathway induces an insufficient residence time of pollutants in the soil, making an incomplete and unfinished biodegradation process, a wrong evaluation of the hydraulic residence time of the system which would hinder its management in a complete system with several entities of treatment and a non-homogeneous growth of the biofilm in the solid filter mass. This paper is a contribution in tracer experiment data analysis within a horizontal flow constructed wetland built in a storm water basin. A two-dimensional numerical model was used to simulate flow and reactive solute transport. The flow model was successfully calibrated in very dry soil conditions. The adsorption profiles used in the reactive transport modeling are those of five molecules: metolachlor, atrazine, deethylatrazine (DEA), deisopropylatrazine (DIA), and hydroxyatrazine (HA). We show that the adsorption distribution is an internal factor of soil which is responsible to the “preferential” pathway transport in a homogeneous gravel texture. The mean residence time of pollutants within the filter is strongly correlated with the average value of the adsorption coefficient. Moreover, we note a lack of significant impact of the heterogeneity of the medium on the statistical moments of breakthrough curve. It appears that a uni-modal breakthrough curve is not significant to the absence of preferential flow in the medium and at least a two-dimensional display can provide sufficient evidence on the presence or absence of preferential pathways. Finally, using a PLS regression, it is possible to perfectly discriminate the number of peaks of concentration and the asymmetry of the breakthrough curve.     

The wave-induced solute flux from submerged sediment

The issue of the transport of dissolved nutrients and contaminants between the sediment in the bottom of a lake or reservoir and the body of water above it is an important one for many reasons. In particular the biological and chemical condition of the body of water is intricately linked to these mass transport processes. As the review by Boudreau (Rev Geophys 38(3):389–416, 2000) clearly demonstrates those transport processes are very complex involving mechanisms as diverse as the wave-induced flux between the sediment and the overlying water and the effect of burrowing animals on the transport within the sediment as well as basic diffusion mechanisms. The present paper focuses on one facet of these transport processes; we re-examine the balance of diffusion and wave-induced advection and demonstrate that the wave-induced flux of a solute from submerged sediment is not necessarily purely diffusive as suggested by Harrison et al. (J Geophys Res 88:7617–7622, 1983) but can be dominated by a mean or time-averaged flux induced by the advective fluid motion into and out of the sediment caused by the fluctuating pressure waves associated with wave motion. Indeed along the subtidal shoreline where the fluctuating bottom pressures are greatest, wave-induced advection will dominate the mean, time-averaged transport of solute into or out of the sediment as suggested in the work of Riedl et al. (Mar Biol 13:210–221, 1972). However, the present calculations also indicate that this advective flux decreases rapidly with increasing depth so that further away from the shoreline the advective flux becomes negligible relative to the diffusive flux and therefore the latter dominates in deeper water.     

Supply-side ecology of the brown mussel,<Emphasis Type="Italic"> Perna perna</Emphasis>: an investigation of spatial and temporal variation in,and coupling between,gamete release and larval supply

Sampling of recruitment-associated variables of Perna perna was done approximately monthly for 14 months at intertidal locations 500 m apart, nested within sites 25 km apart. Paired with intertidal locations were nearshore locations, 600 m to sea. Sampling assessed spawning, densities of larvae in the water column and densities of late plantigrades and juveniles on the shore. Major events in each variable were synchronous over larger scales (10s of kilometres) while subsidiary events were synchronised at smaller scales, varying within sites (100s of metres) or even within locations (metres). This suggests that the processes driving major events operated over large scales while processes operating at much more local scales drove less intense, more localised events. A major spawning event occurred at all locations in May–June 1998. Weaker spawning events occurred at different times in different locations. Larvae were found on 80% of sampling occasions, densities peaking in January–March 1998 and 1999 at all locations. Plantigrades and juveniles showed less clear patterns, with considerable residual variation. There was no sign of strong coupling among variables with few significant direct or cross correlations. The major sources of variability shifted from time to space as one progressed from spawning, to plantigrade density to juvenile density. For spawning, time was the most important source (58%) of heterogeneity and space accounted for little (8%) of the total variance. For larvae and late plantigrades, time was still the most important source of variability (41% and 33%, respectively), but space was a much more substantial component. For juveniles, small-scale (residual) spatial variability dominated total variability (75%). This strongly suggests the importance of hydrography and its effects on variation in delivery of larvae to the intertidal from offshore. These findings also indicate greater spatial heterogeneity as recruits age, reflecting small-scale variations in larval delivery and the increasing importance of post-settlement mortality.Communicated by G.F. Humphrey, Sydney     

The dynamics of phytoplankton blooms in puget sound a fjord in the Northwestern United States

This paper describes a quantitative investigation of relationships between the growth of phytoplankton, and climatic and hydrodynamci conditions in temperate fjords with marked tides, as exemplified by Puget Sound, Washington (USA). Algal growth in the open waters of the central basin of the Sound is dominated by a number of intense blooms beginning in late April or May and recurring throughout the summer. Rarely, and only briefly, does nitrate become exhausted. The phytoplankton production rate in the central basin of Puget Sound is about 465 g C m^-2 year^-1. During the springs of 1966 and 1967, oceanographic measurements were carried out at a mid-channel station with sufficient frequency to allow investigation of physical and biological processes with time scales of the order of a day. The principal investigative tool is a numerical model in which the hydrodynamical conditions are represented by an approximate analysis of the gravitational convection mode of circulation. Algal concentration is represented as a continous function of space and time in the model which ascribes changes in phytoplankton density to variations in photosynthetic and respiratory activity, algal sinking, grazing by herbivores, and to mixing and advection. Computations adequately reproduce the principal features of phytoplankton concentrations observed during 75 days and 35 days in the springs of 1966 and 1967, respectively. Numerical experiments assess the relative importance of various processes which govern the level of primary production in Puget Sound. It is concluded that phytoplankton growth is limited by a combination of factors, including vertical advection and turbulence, modulation of underwater light intensity by self-shading and inorganic particulates, sinking of algal cells, and occasional rapid horizontal advection of the population from the area by sustained winds. The high primary productivity of the Sound is due to intensive upward transport of nitrate by the estuarine mechanism. These results should be generally applicable to other temperate fjords because of the largely conventional choice of the biological functions.     

Hydrology and pollution of lake vistonis in Northern Greece

The basic objective of this study is the understanding of the hydrology and hydraulic characteristics of Lake Vistonis, a very shallow lake in Northern Greece, connected to the Aegean Sea through an irregular shallow channel. In this paper we estimate the pollutant residence time and the discharge of lake to the Aegean Sea as a function of lake elevation. Also we determine the seasonal variations of salinity of the lake.     

Spatial and temporal heterogeneities in the contact behaviour of rabbits

It is becoming increasingly evident that the social behaviour of many group-living species is more complex than previously assumed and that free mixing of individuals, even within social groups, is rare. This has important implications for ecological processes, such as disease transmission, which are dependent on interactions between individuals. European wild rabbits (Oryctolagus cuniculus) have been considered traditionally as highly sociable animals that mix freely within groups but interact less frequently between groups. We deployed proximity logging devices to quantify the intra- and inter-group contact behaviour of free-living wild rabbits in two populations in a temperate region of Australia. Altogether, 126 rabbits were fitted with proximity loggers at least once during the study. Radio-tracking was carried out alongside proximity data collection to determine the space use of rabbit social groups within the study sites. On average, a rabbit made only 1.54 ± 0.23 (SE) (median = 0.54) contacts per day with each other rabbit carrying a proximity logger in its social group, and the mean daily contact duration was 202 ± 38 s/day (SE) (median = 29 s). Despite the high degree of home range overlap between the neighbouring social groups, inter-group contacts were highly infrequent and brief. Our results demonstrated considerable spatial and temporal heterogeneities in the contact behaviour between individual rabbits, both between populations and between and within social groups in the same population. Such variations in the social organisation of rabbits are likely to create complex patterns of disease transmission through direct contact and may contribute towards observed heterogeneities in the effects of disease on wild rabbit populations.     

Variations in the elemental compositions of plants,and computer‐aided sampling in ecosystems∗

To obtain comparable results of multi‐element analysis of plant materials by different laboratories, a harmonized sampling procedure for terrestrial and marine ecosystems is essential. The heterogeneous distribution of chemical elements in living organisms is influenced by different biological parameters. These parameters are mainly characterized by genetic predetermination, seasonal changes, edaphic and climatic conditions, and delocalization processes of chemical substances by metabolic activities.

The biological variations of the element content in plants were divided into 5 systematic levels, which are: 1. the plant species; 2. the population; 3. the stand (within an ecosystem); 4. the individual; and 5. the plant compartment. Each of these systematic levels can be related to: 1. genetic variabilities; 2. different climatic, edaphic and anthropogenic influences; 3. microclimatic or microedaphic conditions; 4. age of plants (stage of development), exposure to environmental influences (light, wind, pollution etc.), seasonal changes; and 5. transport and deposition of substances within the different plant compartments (organs, tissues, cells, organelles).

An expert system for random and systematic sampling for multi‐element analysis of environmental materials, such as plants, soils and precipitation is presented. After statistical division of the research area, the program provides advice for contamination‐free collection of environmental samples.     

Linking chronic wasting disease to mule deer movement scales: a hierarchical Bayesian approach.

Observed spatial patterns in natural systems may result from processes acting across multiple spatial and temporal scales. Although spatially explicit data on processes that generate ecological patterns, such as the distribution of disease over a landscape, are frequently unavailable, information about the scales over which processes operate can be used to understand the link between pattern and process. Our goal was to identify scales of mule deer (Odocoileus hemionus) movement and mixing that exerted the greatest influence on the spatial pattern of chronic wasting disease (CWD) in northcentral Colorado, USA. We hypothesized that three scales of mixing (individual, winter subpopulation, or summer subpopulation) might control spatial variation in disease prevalence. We developed a fully Bayesian hierarchical model to compare the strength of evidence for each mixing scale. We found strong evidence that the finest mixing scale corresponded best to the spatial distribution of CWD infection. There was also evidence that land ownership and habitat use play a role in exacerbating the disease, along with the known effects of sex and age. Our analysis demonstrates how information on the scales of spatial processes that generate observed patterns can be used to gain insight when process data are sparse or unavailable.     

Population dynamics of the barnacle <Emphasis Type="Italic">Chthamalus montagui</Emphasis> at two spatial and temporal scales in northern Spain

This study describes spatial and temporal patterns of variability in population parameters in the barnacle Chthamalus montagui Southward in three localities of northern Spain and evaluates whether density-dependent settlement may regulate population dynamics. The sampling design considered two spatial scales, localities and sites within localities, and two temporal scales, years and six month intervals. Density, amount of free space, mortality, growth rate and magnitude of settlement (both absolute and per unit of free space) were obtained from photographs of permanent quadrats and from direct counts in the field. The number of settlers in scraped and untouched quadrats was used to estimate the importance of the presence of conspecifics in settlement. Significant variation at the two spatial and temporal scales was found for most parameters. Large spatial and temporal variations in adult mortality rate, density, and settlement were observed. Patterns of mortality were not consistent with differences in density among localities. Differences in settlement among localities were maintained through time. We suggest that magnitude of settlement is regulated by persistent features such as topography or local water circulation. We assume that early post-settlement mortality does not differ among localities. In the absence of differential mortality, settlement determines average population density within localities. Within localities, settlement was independent of density and free space. No consistent evidence was found on preferential settlement at the vicinity of conspecifics. The main conclusion is that density-dependent settlement is not relevant for the regulation of the populations of C. montagui in the northern Spain. Regulation might occur by density-dependent processes within the adult fraction of the population and/or the larval phase before settlement.     

Kelvin–Helmholtz Instability as a Boundary-Value Problem

The Kelvin–Helmholtz (KH) instability is traditionally viewed as an initial-value problem, wherein wave perturbations of a two-layer shear flow grow over time into billows and eventually generate vertical mixing. Yet, the instability can also be viewed as a boundary-value problem. In such a framework, there exists an upstream condition where a lighter fluid flows over a denser fluid, wave perturbations grow downstream to eventually overturn some distance away from the point of origin. As the reverse of the traditional problem, this flow is periodic in time and exhibits instability in space. A natural application is the mixing of a warmer river emptying into a colder lake or reservoir, or the salt-wedge estuary. This study of the KH instability from the perspective of a boundary-value problem is divided into two parts. Firstly, the instability theory is conducted with a real frequency and complex horizontal wavenumber, and the main result is that the critical wavelength at the instability threshold is longer in the boundary-value than in the initial-value situation. Secondly, mass, momentum and energy budgets are performed between the upstream, unmixed state on one side, and the downstream, mixed state on the other, to determine under which condition mixing is energetically possible. Cases with a rigid lid and free surface are treated separately. And, although the algebra is somewhat complicated, both end results are identical to the criterion for complete mixing in the initial-value problem.     

Detection of scale-specific community dynamics using wavelets

The response of ecological communities to anthropogenic disturbance is of both scientific and practical interest. Communities where all species respond to disturbance in a similar fashion (synchrony) will exhibit large fluctuations in total biomass and dramatic changes in ecosystem function. Communities where some species increase in abundance while others decrease after disturbance (compensation) can maintain total biomass and ecosystem function in the face of anthropogenic change. We examined dynamics of the Little Rock Lake (Wisconsin, USA) zooplankton community in the context of an experimental pH manipulation conducted in one basin of the lake. A novel application of wavelets was used to partition patterns of synchrony and compensation by time scale. We find interestingly that some time series show both patterns of synchrony and compensation depending on the scale of analysis. Within the unmanipulated basin, we found subtle patterns of synchrony and compensation within the community, largely at a one-year time scale corresponding to seasonal variation. Within the acidified lake basin, dynamics shifted to longer time scales corresponding to the pattern of pH manipulation. Comparisons between pairs of species in different functional groups showed both strong compensatory and synchronous responses to disturbance. The strongest compensatory signal was observed for two species of Daphnia whose life history traits lead to synchrony at annual time scales, but whose differential sensitivity to acidification led to compensation at multiannual time scales. The separation of time scales inherent in the wavelet method greatly facilitated interpretation as patterns resulting from seasonal drivers could be separated from patterns driven by pH manipulation.     

Mercury cycling in litter and soil in different forest types in the Adirondack region, New York, USA.

The fate of mercury in decomposing leaf litter and soil is key to understanding the biogeochemistry of mercury in forested ecosystems. We quantified mercury dynamics in decomposing leaf litter and measured fluxes and pools of mercury in litterfall, throughfall, and soil in two forest types of the Adirondack region, New York, USA. The mean content of total mercury in leaf litter increased to 134% of its original mass during two years of decomposition. The accumulation pattern was seasonal, with significant increases in mercury mass during the growing season (+4.9% per month). Litterfall dominated mercury fluxes into the soil in the deciduous forest, whereas throughfall dominated fluxes into the coniferous forest. The increase in mercury mass in decomposing deciduous litter during the growing season was greater than could be accounted for by throughfall inputs during the growing season (P < 0.05), suggesting translocation of mercury from the soil to the decomposing deciduous litter. This internal recycling mechanism concentrates mercury in the organic horizons and retards transport through the soil, thereby increasing the residence time of mercury in the forest floor. A mass balance assessment suggests that the ultimate fate of mercury in the landscape depends upon forest type and associated differences in the delivery and incorporation of mercury into the soil. Our results show that incorporation of mercury into decaying leaf litter increases its residence time in the landscape and may further delay the recovery of surface waters, fish, and associated biota following control of mercury emissions to the atmosphere.     

Seasonal variation in the balance between physiological mechanisms of feeding and digestion in Mytilus edulis (Bivalvia: Mollusca)

Seasonal variations in the balance between physiological mechanisms of feeding and digestion are considered among 45 to 57 mm shell length Mytilus edulis L. from southwest England. Mussels had been acclimated to standardised conditions of food availability in March, June and October, 1981. Results indicate that despite showing a negative relationship (P<0.001) with the efficiency of ¹⁵N absorption, coincident alterations of both ingestion rate (0.17 to 0.54 mg total Phaeodactylum tricornutum h^-1) and organic gut content (1.00 to 1.80 mg) represent significant mechanisms effecting the seasonal regulation of nutrient acquisition. Concurrent changes in absorption efficiency alone clearly exerted a relatively minor influence upon seasonal rates of absorption. Associated gut residence times of ¹⁵N-labelled material varied between 4 and 10 h. Furthermore, despite the ensuing deleterious effects, together with increased gut contents, absorption efficiencies among mussels that had been starved throughout acclimation were consistently maintained equivalent to those among individuals feeding normally. Such constancy was at least partially achieved by an increased residence time of material within diverticulae, relative to the digestive system as a whole. Finally, comparison of isotopic absorption efficiencies with those recorded for total organics and organic nitrogen has demonstrated differences of as much as 30% between gross and net uptake due, we suggest, to excretion of metabolically derived material within the alimentary canal.     

Acanthochromis polyacanthus,a fish lacking larval dispersal,has genetically differentiated populations at local and regional scales on the Great Barrier Reef

Acanthochromis Gill is a monotypic genus within the damselfish family Pomacentridae, erected for an unusual species [A. polyacanthus (Bleeker)] that uniquely lacks larval dispersal. Instead, offspring are reared in the parental territory, in the manner of cichlids, and fledged into the surrounding habitat. Phenotypic and genotypic variation was surveyed on the basis of body colouration and 7 polymorphic loci in 19 populations from 5 regions of the central and southern Great Barrier Reef (GBR). Variation in both characters was found at regional and local scales. Two colour morphs were recognised: a bicoloured morph from the three northern regions and a uniform dark morph from the two southern regions. Isozyme analysis showed a similar pattern with greatest variation between the different morphs, but also with significant variation at both regional and local scales within morphotypes. Heterozygosity was maximal in the central populations, which, together with other measures of variability, suggests a mixing of separate gene pools in this region and denies species status to the two morphotypes despite numerous fixed differences in allele frequencies between the most distant populations. The presence of fixed differences in multiple alleles between populations separated by 1000 km indicates negligible gene flow over such distances and long isolation of these gene pools. These patterns may reflect recolonisation of the GBR after the last sea-level rise by fish from two stocks. Founder effects and random drift in small populations after colonisation are probably the major sources of the local and regional variations observed at smaller spatial scales. This diversity has been maintained among populations at all scales by the very low levels of gene flow possible without an effective strategy for larval dispersal between coral reefs.     

Influence of human activity and monsoon dynamics on spatial and temporal hydrochemistry in tropical coastal waters (Sanya Bay,South China Sea)

Human activities and oceanic influences like mixing and/or upwelling define the hydrochemical and biological characteristics of coastal regions. In Sanya Bay, northern South China Sea, the balance between these two influences on spatial and temporal scales is poorly understood. The influence from human activities was due to discharge from Sanya River, whereas the most important marine influence was related to seasonal changes in hydrodynamics. Spatial differences (p<0.05) in total nitrogen (TN) and chlorophyll-a (Chl-a) were observed. Seasonal differences were observed for temperature and nutrients. Human activities are the dominant factor influemcing hydrochemistry in the inner bay. This region exhibited the maximum influence of discharge from Saya River as estimated by higher nutrient levels. Oceanic influences like upwelling and mixing caused by monsoons are the dominant factors influencing hydrochemistry in the central and outer bay. Both human activities and oceanic influences play important roles in coastal ecosystems in Sanya Bay. This study has defined these characteristics so that better management policies can be enacted.