Flow controls on lowland river macrophytes: A review

We review the current status of knowledge regarding the role that flow parameters play in controlling the macrophyte communities of temperate lowland rivers. We consider both direct and indirect effects and the interaction with other factors known to control macrophyte communities. Knowledge gaps are identified and implications for the management of river systems considered. The main factors and processes controlling the status of macrophytes in lowland rivers are velocity (hence also discharge), light, substrate, competition, nutrient status and river management practices. We suggest that whilst the characteristics of any particular macrophyte community reflect the integral effects of a combination of the factors, fundamental importance can be attributed to the role of discharge and velocity in controlling instream macrophyte colonisation, establishment and persistence. Velocity and discharge also appear to control the relative influence of some of the other controlling factors. Despite the apparent importance of velocity in determining the status of macrophyte communities in lowland rivers, relatively little is understood about the nature of the processes controlling this relationship. Quantitative knowledge is particularly lacking. Consequently, the ability to predict macrophyte abundance and distribution in rivers is still limited. This is further complicated by the likely existence of feedback effects between the growth of macrophytes and velocity. Demand for water resources increases the pressure on lowland aquatic ecosystems. Despite growing recognition of the need to allocate water for the needs of instream biota, the inability to assess the flow requirements of macrophyte communities limits the scope to achieve this. This increases the likelihood of overexploitation of the water resource as other users, whose demands are quantifiable, are prioritised. (C) 2008 Elsevier B.V. All rights reserved.

Riparian shading controls instream spring phytoplankton and benthic diatom growth

Dissolved oxygen (DO) concentrations showed a striking pattern in a multi-year study of the River Enborne, a small river in SE England. In each of three years (2010-2012), maximum DO concentrations were attained in mid-April, preceded by a period of steadily increasing diurnal amplitudes, followed by a steady reduction in both amplitude and concentration. Flow events during the reduction period reduce DO to low concentrations until the following spring. Evidence is presented that this pattern is mainly due to benthic algal growth which is eventually supressed by the growth of the riparian tree canopy. Nitrate and silicate concentrations are too high to inhibit the growth of either benthic algae or phytoplankton, but phosphate concentrations might have started to reduce growth if the tree canopy development had been delayed. This interpretation is supported by evidence from weekly flow cytometry measurements and analysis of the diurnal, seasonal and annual patterns of nutrient concentrations. As the tree canopy develops, the river switches from an autotrophic to a heterotrophic state. The results support the use of riparian shading to help control algal growth, and highlight the risks of reducing riparian shade.

The influence of viscous and latent heating on crystal-rich magma flow in a conduit

The flow dynamics of crystal-rich high-viscosity magma is likely to be strongly influenced by viscous and latent heat release. Viscous heating is observed to play an important role in the dynamics of fluids with temperature-dependent viscosities. The growth of microlite crystals and the accompanying release of latent heat should play a similar role in raising fluid temperatures. Earlier models of viscous heating in magmas have shown the potential for unstable (thermal runaway) flow as described by a Gruntfest number, using an Arrhenius temperature dependence for the viscosity, but have not considered crystal growth or latent heating. We present a theoretical model for magma flow in an axisymmetric conduit and consider both heating effects using Finite Element Method techniques. We consider a constant mass flux in a 1-D infinitesimal conduit segment with isothermal and adiabatic boundary conditions and Newtonian and non-Newtonian magma flow properties. We find that the growth of crystals acts to stabilize the flow field and make the magma less likely to experience a thermal runaway. The additional heating influences crystal growth and can counteract supercooling from degassing-induced crystallization and drive the residual melt composition back towards the liquidus temperature. We illustrate the models with results generated using parameters appropriate for the andesite lava dome-forming eruption at Soufriere Hills Volcano, Montserrat. These results emphasize the radial variability of the magma. Both viscous and latent heating effects are shown to be capable of playing a significant role in the eruption dynamics of Soufriere Hills Volcano. Latent heating is a factor in the top two kilometres of the conduit and may be responsible for relatively short-term (days) transients. Viscous heating is less restricted spatially, but because thermal runaway requires periods of hundreds of days to be achieved, the process is likely to be interrupted. Our models show that thermal evolution of the conduit walls could lead to an increase in the effective diameter of flow and an increase in flux at constant magma pressure.

Impact of Tibetan Orography and Heating on the Summer Flow over Asia

The influence on the summer flow over Asia of both the orographic and thermal forcing of the Tibetan Plateau is investigated using a sequence of idealised experiments with a global primitive equation model. The zonally averaged flow is prescribed and both realistic and idealised orography and heating are used. There is some similarity between the responses to the two forcings when applied separately. The upper tropospheric Tibetan anticyclone is predominantly forced by the heating but also weakly by the orography. Below this, both forcings tend to give air descending in an equatorward anticyclonic circulation down the isentropes to the west and rising in a similar poleward circulation to the east. However the heating-only response has a strong ascending southwesterly flow that is guided around the south and south-east of the orography when it is included. On the northern side, the westerly flow over the orography gives ascent on the upslope and descent on the downslope. It is found that heating over the Plateau leads to a potential vorticity (PV) minimum and that if it is sufficiently strong the flow is unstable, producing a quasi-biweekly oscillation. During this oscillation the Tibetan anticyclone changes between a single centre over the southwestern side of the Plateau and a split/double structure with centres over China and the Middle East. These characteristics are similar to observed variability in the region. Associated with this quasi-biweekly oscillation are significant variations in the strength of the ascent over the Plateau and the Rossby wave pattern over the North Pacific. The origin of the variability is instability associated with the zonally extended potential vorticity PV minimum on a θ-surface, as proposed by Hsu and Plumb (2000). This minimum is due to the tendency to reduce the PV above the heating over the Plateau and to advection by the consequent anticyclone of high PV around from the east and low PV to the west. The deep convection to the south and southeast of the Plateau tends to suppress the quasi-biweekly oscillation because the low PV produced above it acts to reduce the meridional PV gradient reversal. The occurrence of the oscillation depends on the relative magnitude of the heating in the two regions.

Spatial variability of flow statistics within regular building arrays

Turbulence statistics obtained by direct numerical simulations are analysed to investigate spatial heterogeneity within regular arrays of building-like cubical obstacles. Two different array layouts are studied, staggered and square, both at a packing density of λp=0.25 . The flow statistics analysed are mean streamwise velocity ( u− ), shear stress ( u′w′−−−− ), turbulent kinetic energy (k) and dispersive stress fraction ( u˜w˜ ). The spatial flow patterns and spatial distribution of these statistics in the two arrays are found to be very different. Local regions of high spatial variability are identified. The overall spatial variances of the statistics are shown to be generally very significant in comparison with their spatial averages within the arrays. Above the arrays the spatial variances as well as dispersive stresses decay rapidly to zero. The heterogeneity is explored further by separately considering six different flow regimes identified within the arrays, described here as: channelling region, constricted region, intersection region, building wake region, canyon region and front-recirculation region. It is found that the flow in the first three regions is relatively homogeneous, but that spatial variances in the latter three regions are large, especially in the building wake and canyon regions. The implication is that, in general, the flow immediately behind (and, to a lesser extent, in front of) a building is much more heterogeneous than elsewhere, even in the relatively dense arrays considered here. Most of the dispersive stress is concentrated in these regions. Considering the experimental difficulties of obtaining enough point measurements to form a representative spatial average, the error incurred by degrading the sampling resolution is investigated. It is found that a good estimate for both area and line averages can be obtained using a relatively small number of strategically located sampling points.

Structure of turbulent flow over arrays of cubical roughness

The structure of turbulent flow over large roughness consisting of regular arrays of cubical obstacles is investigated numerically under constant pressure gradient conditions. Results are analysed in terms of first- and second-order statistics, by visualization of instantaneous flow fields and by conditional averaging. The accuracy of the simulations is established by detailed comparisons of first- and second-order statistics with wind-tunnel measurements. Coherent structures in the log region are investigated. Structure angles are computed from two-point correlations, and quadrant analysis is performed to determine the relative importance of Q2 and Q4 events (ejections and sweeps) as a function of height above the roughness. Flow visualization shows the existence of low-momentum regions (LMRs) as well as vortical structures throughout the log layer. Filtering techniques are used to reveal instantaneous examples of the association of the vortices with the LMRs, and linear stochastic estimation and conditional averaging are employed to deduce their statistical properties. The conditional averaging results reveal the presence of LMRs and regions of Q2 and Q4 events that appear to be associated with hairpin-like vortices, but a quantitative correspondence between the sizes of the vortices and those of the LMRs is difficult to establish; a simple estimate of the ratio of the vortex width to the LMR width gives a value that is several times larger than the corresponding ratio over smooth walls. The shape and inclination of the vortices and their spatial organization are compared to recent findings over smooth walls. Characteristic length scales are shown to scale linearly with height in the log region. Whilst there are striking qualitative similarities with smooth walls, there are also important differences in detail regarding: (i) structure angles and sizes and their dependence on distance from the rough surface; (ii) the flow structure close to the roughness; (iii) the roles of inflows into and outflows from cavities within the roughness; (iv) larger vortices on the rough wall compared to the smooth wall; (v) the effect of the different generation mechanism at the wall in setting the scales of structures.

The effect of flow field & geometry on the dynamic contact angle

A number of recent experiments suggest that, at a given wetting speed, the dynamic contact angle formed by an advancing liquid-gas interface with a solid substrate depends on the flow field and geometry near the moving contact line. In the present work, this effect is investigated in the framework of an earlier developed theory that was based on the fact that dynamic wetting is, by its very name, a process of formation of a new liquid-solid interface (newly “wetted” solid surface) and hence should be considered not as a singular problem but as a particular case from a general class of flows with forming or/and disappearing interfaces. The results demonstrate that, in the flow configuration of curtain coating, where a liquid sheet (“curtain”) impinges onto a moving solid substrate, the actual dynamic contact angle indeed depends not only on the wetting speed and material constants of the contacting media, as in the so-called slip models, but also on the inlet velocity of the curtain, its height, and the angle between the falling curtain and the solid surface. In other words, for the same wetting speed the dynamic contact angle can be varied by manipulating the flow field and geometry near the moving contact line. The obtained results have important experimental implications: given that the dynamic contact angle is determined by the values of the surface tensions at the contact line and hence depends on the distributions of the surface parameters along the interfaces, which can be influenced by the flow field, one can use the overall flow conditions and the contact angle as a macroscopic multiparametric signal-response pair that probes the dynamics of the liquid-solid interface. This approach would allow one to investigate experimentally such properties of the interface as, for example, its equation of state and the rheological properties involved in the interface’s response to an external torque, and would help to measure its parameters, such as the coefficient of sliding friction, the surface-tension relaxation time, and so on.

Inertia-Gravity Waves Emitted from Balanced Flow: Observations, Properties, and Consequences

This paper describes laboratory observations of inertia–gravity waves emitted from balanced fluid flow. In a rotating two-layer annulus experiment, the wavelength of the inertia–gravity waves is very close to the deformation radius. Their amplitude varies linearly with Rossby number in the range 0.05–0.14, at constant Burger number (or rotational Froude number). This linear scaling challenges the notion, suggested by several dynamical theories, that inertia–gravity waves generated by balanced motion will be exponentially small. It is estimated that the balanced flow leaks roughly 1% of its energy each rotation period into the inertia–gravity waves at the peak of their generation. The findings of this study imply an inevitable emission of inertia–gravity waves at Rossby numbers similar to those of the large-scale atmospheric and oceanic flow. Extrapolation of the results suggests that inertia–gravity waves might make a significant contribution to the energy budgets of the atmosphere and ocean. In particular, emission of inertia–gravity waves from mesoscale eddies may be an important source of energy for deep interior mixing in the ocean.