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.     

Turbulence and aeration in hydraulic jumps: free-surface fluctuation and integral turbulent scale measurements

In an open channel, a change from a supercritical to subcritical flow is a strong dissipative process called a hydraulic jump. Herein some new measurements of free-surface fluctuations of the impingement perimeter and integral turbulent time and length scales in the roller are presented with a focus on turbulence in hydraulic jumps with a marked roller. The observations highlighted the fluctuating nature of the impingement perimeter in terms of both longitudinal and transverse locations. The results showed further the close link between the production and detachment of large eddies in jump shear layer, and the longitudinal fluctuations of the jump toe. They highlighted the importance of the impingement perimeter as the origin of the developing shear layer and a source of vorticity. The air–water flow measurements emphasised the intense flow aeration. The turbulent velocity distributions presented a shape similar to a wall jet solution with a marked shear layer downstream of the impingement point. The integral turbulent length scale distributions exhibited a monotonic increase with increasing vertical elevation within 0.2 < L_z/d₁ < 0.8 in the shear layer, where L_z is the integral turbulent length scale and d₁ the inflow depth, while the integral turbulent time scales were about two orders of magnitude smaller than the period of impingement position longitudinal oscillations.     

A Quantitative Analysis of Energy Dissipation among Three Typical Air Entrainment Phenomena

This study investigates energy dissipation due to air bubble entrainment for three typical phenomena; a hydraulic jump, a 2-D vertical plunging jet and a vertical circular plunging jet into water. A simple model is presented here which enables to estimate the energy transformation and dissipation achieved by air bubbles quantitatively for three above phenomena. The average rate of energy dissipation by air bubbles obtained from the experimental data are 25%, 1.4%, and 2.15% with respect to total energy loss for the hydraulic jump, 2-D vertical plunging jet and vertical circular plunging jet, respectively.     

Stability and Mixing of a Vertical Axisymmetric Buoyant Jet in Shallow Water

The stability, mixing and effect of downstream control on axisymmetric turbulent buoyant jets discharging vertically into shallow stagnant water is studied using 3D Reynolds-averaged Navier–Stokes equations (RANS) combined with a buoyancy-extended k –ε model. The steady axisymmetric turbulent flow, temperature (or tracer concentration) and turbulence fields are computed using the finite volume method on a high resolution grid. The numerical predictions demonstrate two generic flow patterns for different turbulent heated jet discharges and environmental parameters (i) a stable buoyant discharge with the mixed fluid leaving the vertical jet region in a surface warm water layer; and (ii) an unstable buoyant discharge with flow recirculation and re-entrainment of heated water. A stratified counterflow region always appears in the far-field for both stable and unstable buoyant discharges. Provided that the domain radius L exceeds about 6H, the near field interaction and hence discharge stability is governed chiefly by the jet momentum length scale to depth ratio l_M/H, regardless of downstream control. The near field jet stability criterion is determined to be l_M/H = 3.5. A radial internal hydraulic jump always exists beyond the surface impingement region, with a 3- to 6-fold increase in dilution across the jump compared with vertical buoyant jet mixing. The predicted stability category, velocity and temperature/concentration fields are well-supported by experiments of all previous investigators.     

Turbulent air–water flows in hydraulic structures: dynamic similarity and scale effects

In hydraulic structures, free-surface aeration is commonly observed: i.e., the white waters. The air bubble entrainment may be localised (hydraulic jumps, plunging jets) or continuous along an interface (water jets, chutes). Despite recent advances, there are some basic concerns about the extrapolation of laboratory results to large size prototype structures. Herein the basic air bubble entrainment processes are reviewed and the relevant dynamic similarities are discussed. Traditionally, physical studies are conducted using a Froude similitude which implies drastically smaller laboratory Reynolds numbers than in the corresponding prototype flows. Basic dimensional analyses are developed for both singular and interfacial aeration processes. The results are discussed in the light of systematic investigations and they show that the notion of scale effects is closely linked with the selection of relevant characteristic air–water flow properties. Recent studies of local air–water flow properties highlight that turbulence levels, entrained bubble sizes and interfacial areas are improperly scaled based upon a Froude similitude even in large-size models operating with the so defined Reynolds numbers ρ _w × q _w/μ _w up to 5 E+5. In laboratory models, the dimensionless turbulence levels, air–water interfacial areas and mass transfer rates are drastically underestimated.     

Impact of consistent boundary layer mixing approaches between NAM and CMAQ

Discrepancies in grid structure, dynamics and physics packages in the offline coupled NWS/NCEP NAM meteorological model with the U.S. Environmental Protection Agency Community Multiscale Air Quality (CMAQ) model can give rise to inconsistencies. This study investigates the use of three vertical mixing schemes to drive chemistry tracers in the National Air Quality Forecast Capability (NAQFC). The three schemes evaluated in this study represent various degrees of coupling to improve the commonality in turbulence parameterization between the meteorological and chemistry models. The methods tested include: (1) using NAM predicted TKE-based planetary boundary height, h, as the prime parameter to derive CMAQ vertical diffusivity; (2) using the NAM mixed layer depth to determine h and then proceeding as in (1); and (3) using NAM predicted vertical diffusivity directly to parameterize turbulence mixing within CMAQ. A two week period with elevated surface O₃ concentrations during the summer 2006 has been selected to test these schemes in a sensitivity study. The study results are verified and evaluated using the EPA AIRNow monitoring network and other ozonesonde data. The third method is preferred a priori as it represents the tightest coupling option studied in this work for turbulent mixing processes between the meteorological and air quality models. It was found to accurately reproduce the upper bounds of turbulent mixing and provide the best agreement between predicted h and ozonesonde observed relative humidity profile inferred h for sites investigated in this study. However, this did not translate into the best agreement in surface O₃ concentrations. Overall verification results during the test period of two weeks in August 2006, did not show superiority of this method over the other 2 methods in all regions of the continental U.S. Further efforts in model improvement for the parameterizations of turbulent mixing and other surface O₃ forecast related processes are warranted.     

Analysis of the velocity field in a large rectangular channel with lateral shockwave

In this work the authors describe the main characteristics of the velocity field of hydraulic jumps in a very large channel where lateral shockwaves occur. Experiments were carried out at the Coastal Engineering Laboratory of the Water Engineering and Chemistry Department of the Technical University of Bari (Italy). Extensive flow velocity measurements were investigated in order to have a clearer understanding of both hydraulic jump development and lateral shockwave formation in a very large channel. Eight experiments were performed in a 4m wide rectangular channel; the experiments differed in the inlet Froude number F ₀ and the jump type. Seven tests were carried out with undular jumps and one with a roller jump. The flow velocity and the flow free surface measurements were taken using a two-dimensional Acoustic Doppler Velocimeter (ADV) and an ultrasonic profiler, respectively. The experimental results can be summarized as follow: (i) the formation of well developed lateral shockwaves similar to those of oblique jumps were observed; (ii) the comparison of the experimental and theoretical data shows that the classic shockwave theory is sufficiently confirmed in the analyzed range of Reynolds number, taking into account the experimental errors and the difference between the theoretical and experimental assumptions; (iii) the transversal flow velocity profiles in the recirculating zone show a good agreement with the numerical simulations presented in literature in the case of a separated turbulent boundary layer over a flat plate. This conclusion enables us to confirm the hypothesis that the lateral shockwaves in the channel are the result of a boundary layer which, as observed, forms on the channel sidewalls.     

Modelling hydraulic jump using the bubbly two-phase flow method

Hydraulic jumps have complex flow structures, characterised by strong turbulence and large air contents. It is difficult to numerically predict the flows. It is necessary to bolster the existing computer models to emphasise the gas phase in hydraulic jumps, and avoid the pitfall of treating the phenomenon as a single-phase water flow. This paper aims to improve predictions of hydraulic jumps as bubbly two-phase flow. We allow for airflow above the free surface and air mass entrained across it. We use the Reynolds-averaged Navier–Stokes equations to describe fluid motion, the volume of fluid method to track the interface, and the k–ε model for turbulence closure. A shear layer is shown to form between the bottom jet flow and the upper recirculation flow. The key to success in predicting the jet flow lies in formulating appropriate bottom boundary conditions. The majority of entrained air bubbles are advected downstream through the shear layer. Predictions of the recirculation region’s length and air volume fraction within the layer are validated by available measurements. The predictions show a linear growth of the shear layer. There is strong turbulence at the impingement, and the bulk of the turbulence kinetic energy is advected to the recirculation region via the shear layer. The predicted bottom-shear-stress distribution, with a peak value upstream of the toe of the jump and a decaying trend downstream, is realistic. This paper reveals a significant transient bottom shear stress associated with temporal fluctuations of mainly flow velocity in the jump. The prediction method discussed is useful for modelling hydraulic jumps and advancing the understanding of the complex flow phenomenon.     

Sensitivity of air quality model prediction to parameterization of vertical eddy diffusivity

This paper investigates the effects of vertical eddy diffusivities derived from the 3 different planetary boundary layer (PBL) schemes on predictions of chemical components in the troposphere of East Asia. Three PBL schemes were incorporated into a regional air quality model (RAQM) to represent vertical mixing process and sensitivity simulations were conducted with the three schemes while other options are identical. At altitudes <2km, all schemes exhibit similar skill for predicting SO₂ and O₃, but more difference in the predicted NO_x concentration. The Gayno–Seaman scheme produces the smallest vertical eddy diffusivity (Kz) among all schemes, leading to higher SO₂ and NO_x concentrations near the surface than that from the other 2 schemes. However, the effect of vertical mixing on O₃ concentration is complex and varies spatially due to chemistry. The Gayno–Seaman scheme predicts lower O₃ concentrations than the other two schemes in the parts of northern China (generally VOC-limited) and higher ones in most parts of southern China (NO_x-limited). The Byun and Dennis scheme produces the largest mixing depth in the daytime, which bring more NO_x into upper levels, and the mixing depth predicted by the Gayno–Seaman scheme is the smallest, leading to higher NO_x and lower O₃ concentrations near the surface over intensive emission regions.     

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.     

Self-aeration in the rapidly- and gradually-varying flow regions of steep smooth and stepped spillways

In high-velocity chute flows, free-surface aeration is often observed. The phenomenon is called self-aeration or white waters. When the turbulent shear stresses next to the free-surface are large enough, air bubbles are entrained throughout the entire air–water column. A rapidly-varied flow region is observed immediately downstream of the inception point of free-surface aeration. An analytical solution of the air diffusion equation is proposed and the results compare well with new experimental data. Both experiments and theory indicate that the flow bulking spans over approximately 3–4 step cavities downstream of the inception point of free-surface aeration on a stepped chute. Further downstream the void fraction distributions follow closely earlier solutions of the air diffusion equation. The application of the diffusion equation solution to prototype and laboratory data shows air bubble diffusivities typically larger than the momentum transfer coefficient. The result highlights however a marked decrease in the ratio of air bubble diffusivity to eddy viscosity with increasing Reynolds number. The finding might indicate some limitation of laboratory investigations of air bubble diffusion process in self-aerated flows and of their extrapolation to full-scale prototype applications.     

Experimental Study of Gas-Liquid Interfacial Properties in a Stepped Cascade Flow

Gas-liquid interface measurements were conducted in a strongly turbulent free-surface flow (i.e., stepped cascade). Local void fractions, bubble count rates, bubble size distributions and gas-liquid interface areas were measured simultaneously in the air-water flow region using resistivity probes. The results highlight the air-water mass transfer potential of a stepped cascade with measured specific interface area over 650 m^–1 and depth-average specific area up to 310 m^–1. A comparison between single-tip and double-tip resistivity probes suggests that simple robust single-tip probes may provide accurate, although conservative, gas-liquid interfacial properties. The latter device may be used in the field and in prototype plants. Notation a = specific interface area (m^–1); a _mean = depth-average specific interface area (m^–1): a _mean=frac1Y ₉₀limits sup> Y ₉₀ sup ₀(1–C)dy; C = local void fraction; C _gas = dissolved gas concentration (kg m^–3); C _mean = depth-average mean air concentration defined as: C _mean=1–d/Y ₉₀; C _s = saturation concentration (kg m^–3); D = dimensionless air bubble diffusivity (defined by [1]); d = equivalent clear-water flow depth (m): d=limits sup> Y ₉₀ sup ₀(1–C) dy; d_ab = air bubble diameter (m); d_c = critical flow depth (m); for a rectangular channel: d _c=sqrt[3]q _w ²/g; F = air bubble count rate (Hz); F _max = maximum bubble count rate (Hz), often observed for C=50%; g = gravity acceleration (m s^–2); h = step height (m); K _L = liquid film coefficient (m s^–1); K = integration constant defined as: K=tanh ^–1 sqrt0.1)+(2D)^–1 [1]; L = chute length (m); N = velocity distribution exponent; ———– ^*Corresponding author, E-mail: h.chanson@mailbox.uq.edu.au Q _w = water discharge (m³ s^–1); q _w = water discharge per unit width m² s^–1); t = time (s); V = local velocity (m s^–1); V _c = critical flow velocity (m s^–1); for a rectangular channel: V _c=sqrt[3]q _w g V _max = maximum air-water velocity (m s^–1); V ₉₀ = characteristic air-water velocity (m s^–1) where C = 90%; W = channel width (m); x = longitudinal distance (m) measured along the flow direction (i.e., parallel to the pseudo-bottom formed by the step edges); y = distance (m) normal to the pseudo-bottom formed by the step edges; Y₉₀ = characteristic distance (m) where C=0.90; Y ₉₈ = characteristic distance (m) where C=0.98; = slope of pseudo-bottom by the step edges; = diameter (m).     

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.     

Calculation of Entrainment in Density Jumps

The density jump in a two-layer channel flow of miscible fluids, in which one of the layers is infinitely deep and at rest, is analyzed using the momentum flux and mass flux conservation equations. The analysis yields simple equations relating the heights upstream and downstream of the jump with the upstream Froude number and the rate of entrainment into the moving layer, as well as a relation between the maximum possible entrainment and the upstream Froude number. The analysis also shows that when the flow down-stream of the jump is controlled by an obstruction or a contraction, the entrainment rate into the jump depends solely on the dimensionless obstruction height contraction ratio and the upstream Froude number.     

Hydrodynamics,diffusion-boundary layers and photosynthesis of the seagrasses Thalassia testudinum and Cymodocea nodosa

Photosynthetic rates of aquatic plants frequently increase with increasing current velocities. This is presumably due to a reduction in the thickness of the diffusion boundary-layer which allows for a higher carbon availability on the plant surface. Blades of the seagrasses Thalassia testudinum and Cymodocea nodosa exposed to different current velocities under controlled laboratory conditions, showed increased photosynthetic rates with increasing flow only at low current velocities (expressed as blade friction velocities, u _*). Carbon saturation of photosynthetic processes occurred at a relatively low u _* level (0.25 cm s^-1) for T. testudinum collected from a calm environment compared to C. nodosa (0.64 cm s^-1) collected from a surf zone. No further enhancement of photosynthetic rates was observed at higher u _* levels, suggesting limitations in carbon diffusion through the boundary layer below critical u ^* levels and possible limitations in carbon fixation by the enzymatic system at higher u _* levels. These results, as well as those of previous theoretical studies, assumed the flow on the immediate seagrass-blade surface to be hydrodynamically smooth. The presence of epiphytes and attached debris causes the surface of in situ seagrass blades to be exposed to flows ranging from smooth to rough-turbulent. As a consequence, the boundary-layer thickness on moderately epiphytized blades under medium to high flow-conditions is not continuous, but fluctuates in time and space, enhancing carbon transport. In situ u _* levels measured directly on blades of seagrasses indicate that T. testudinum and C. nodosa can be exposed to conditions under which the boundary layer limits photosynthesis during short periods of time (milliseconds) during low-energy events. As waves cause the thickness of the diffusion boundary-layer to fluctuate constantly, carbon-limiting conditions do not persist for prolonged periods.     

Transport and Diffusion of Ozone in the Nocturnal and Morning Planetary Boundary Layer of the Phoenix Valley

The evolution of ozone (O ₃) in the nocturnal and morning-transitional planetary boundary layer (PBL) of the Phoenix valley was measured as a part of the `Phoenix Sunrise Experiment 2001' of the U.S. Department of Energy conducted in June 2001. The goal of the field program was to study the transport, distribution and storage of ozone and its precursors in the urban boundary layer over a diurnal cycle. The ground level O <sub>3</sub> as well as mean meteorological variables and turbulence were measured over the entire period, and vertical profiling (using a tethered balloon) was made during the morning transition period. Approximately half of the observational days showed the usual diurnal cycle of high O <sub>3</sub> during the day and low O <sub>3</sub> at night, with nitrogen oxides (NO <sub>x</sub> = NO <sub>2</sub> + NO) showing an out of phase relationship with O <sub>3</sub>. The rest of the days were signified by an anomalous increase of O <sub>3</sub> in the late evening ( 2200 LST), concomitant with a sudden drop of temperature, an enhancement of wind speed and Reynolds stresses, a positive heat flux and a change of wind direction. NO <sub>x</sub> measurements indicated the simultaneous arrival of an `aged' air mass, which was corroborated by the wind predictions of a mesoscale numerical model. In all, the results indicate that the recirculation of O ₃ rich air masses is responsible for the said high-O ₃ events. Such air masses are produced during the transport of O ₃ precursors by the upslope flow toward mountainous suburbs during the day, and they return back to the city at night via downslope winds (i.e. mountain breeze). The corresponding flow patterns, and hence the high-O ₃ events, are determined by background meteorological conditions. The vertical profiling of O ₃ and flow variables during the morning transition points to a myriad of transport, mixing and chemical processes that determine the fate of tropospheric O ₃. How well such processes are incorporated and resolved in predictive O ₃ models should determine the accuracy of their predictions.     

Mélange vertical et capacité photosynthétique du phytoplancton estuarien (estuaire du Saint-Laurent)

Time series of chlorophyll a, photosynthetic capacity and many physical parameters were sampled hourly for 167 h in August, 1975, at an anchor station located in the Middle Estuary of the St. Lawrence River, Canada. Sampling was carried out during the transition from neap tides to spring tides. The long-and short-term variations in chlorophyll a are coupled with the advection of water masses which depends on tidal currents. Vertical mixing also influences the chlorophyll a concentration of the cells, since it modifies the physiological state of the phytoplankton. Furthermore, circadian periodicities were observed in the photosynthetic capacity, suggesting that the phytoplankton of this area have a homogeneous light history due to strong vertical mixing. Under these conditions, the photosynthetic capacity is adapted to the mean light intensity in the mixed layer. The semimonthly (M _f) variations of the mean light intensity in the mixed layer depend on the M _f variations in the vertical mixing, whereas in the short-term, the variations in mean light intensity in the mixed layer are circadian.

---

Mélange vertical et capacité photosynthétique du phytoplancton estuarien (estuaire du Saint-Laurent)

---

Contribution au programme du Groupe interuniversitaire de recherches océanographiques du Québec (GIROQ)     

Double-average mean flow and local turbulence intensity profiles from PIV measurements for an open channel flow with rigid vegetation

This contribution presents particle image velocimetry measurements for an open channel stationary uniform and fully developed flow of water over a horizontal flat bed of uniform glass beads in presence of a staggered array of vertical cylindrical stems. The main objective was to explore and quantify the influence of the stems-to-flow relative submergence, h _v /h, over the mean flow and local turbulence intensities. A comparison with measurements for the non-vegetated flow over the same granular bed is presented. Results indicate a remarkable influence of h _v /h over the whole flow field. The time-average mean flow presents a strong spatial variation in the layer of the flow occupied by the stems. The local velocity fluctuations are strongly affected by the presence of the stems, with regions in between the stems where they reach peaks that are several times larger than those encountered in the flow in absence of vegetation. The turbulence intensity profiles are noticeably different when compared to those measured in the non-vegetated flow conditions. From previous works it was possible to derive an equation for the mean velocity, U _v , of the flow through the vegetated layer of height h _v . The prediction of this equation is in good agreement with the uniform value for the double-average longitudinal velocity profile in this layer. A final brief discussion about the possible impact of these vegetated-flow features on the sediment transport is presented.     

Vertical distribution of mesoplankton in the open area of the Black Sea

In April–May 1984 mesoplankton vertical distribution in the Black Sea was studied by sampling with a 150-l waterbottle, vertical plankton nets with mesh-sizes of 180 and 500 mkm and by direct counting of the jelly-fish Aurelia aurita, the ctenophore Pleurobrachia pileus, Calanus helgolandicus and the chaetognath Sagitta setosa from the manned submersible Argus. During daytime throughout the whole deep-water body of the sea near the lower oxycline boundary, plankton forms a layer of high concentration (from 2.5 to 38 g m^-3); its thickness varies from 5 to 10–20 m and it has an unchangeable vertical structure; its upper portion is formed by the ctenophore P. pileus, its middle portion by V–VI copepodites of C. helgolandicus, and its lower portion by the population of S. setosa. The lower boundary of this layer coincides with 0.4 to 0.5 ml O₂ l^-1 isooxygen surface, and the depth of its location varies in different areas of the sea from 150 to 50 m, depending on the depth location of 0.5 ml O₂ l^-1 isooxygen surface. By night the animals, which form this layer, migrate towards the surface.     

Copper, Lead and Zinc Mobility and Bioavailability in a River Sediment Contaminated with Paint Stripping Residue

A river embankment sediment beneath a road bridge in South Australia had been surveyed on a previous occasion to a depth of 10 cm and found to contain high levels of copper, lead and zinc deposited during previous lead and zinc-based primer paint stripping operations prior to repainting the bridge. Since the river embankment is geologically different from the surrounding non-embankment soils, and a sufficient distance from the bridge along the river embankment could not be traversed, the geological background Cu, Pb and Zn fraction was established initially using a solid phase sequential extraction (SPSE) procedure. Assessment of the degree of contamination and extent of vertical mobility was then established. Copper was a minor contaminant and present predominantly as the geological background species. Lead was partitioned evenly between the oxalate soluble fraction and residual fraction to a depth of 30 cm with a decrease in total Pb concentration from 497 to 141 mg Pb kg^–1 while Zn was predominantly in the oxalate soluble fraction to a depth of 15cm with a decrease in concentration from 1013 to 150 mg Zn kg^–1. The reduced rate of nitrification and the lower total concentration of NO₃ ^–-N together with the higher respiratory quotient, combined with the lower microbial quotient, are indicative of microbial stress due to heavy metal toxicity. In the case of Pb and Zn paint stripping residues, these changes in indicators of microbial health are possibly the result of the greater abundance and hence possible bioavailabilty of organo- and chloro-Pb compared to organo- and chloro-Zn complexes.