Hysteresis‐based analysis of overland metal transport

Introducing a concept of equivalent mass depth of flow, this study describes the phenomenon of non‐point source pollutant (metal) transport for pavement (or overland) flow in analogy with wave propagation in wide open channels. Hysteretic and normal mass rating curves are developed for runoff rate and mass of 12 dissolved and particulate‐bound metal elements (pollutants) using the rainfall‐runoff and water quality data of the 15 × 20 m² instrumented pavement in Cincinnati, USA. Normal mass rating curves developed for easy computation of pollutant load are found to be of a form similar to Manning's, which is valid for open channel flows. Based on the hysteresis analysis, wave types for dissolution and mixing of particulate‐bound metals are identified. The analysis finds that the second‐order partial‐differential equation normally used for metal transport does not have the efficacy to describe fully the strong non‐linear phenomena such as is described for various metal elements by dynamic waves. In addition, the proportionality concept of the popular SCS‐CN concept is extended for determining the potential maximum metal mass M_p of all the 12 elements transported by a rain storm and related to the antecedent dry period (ADP). For the primary metal zinc element, M_p is found to increase with the ADP and vice versa. Copyright © 2003 John Wiley & Sons, Ltd.     

Modelling the effect of low soil temperatures on transpiration by Scots pine

For ecosystem modelling of the Boreal forest it is important to include processes associated with low soil temperature during spring‐early summer, as these affect the tree water uptake. The COUP model, a physically based SVAT model, was tested with 2 years of soil and snow physical measurements and sap flow measurements in a 70‐year‐old Scots pine stand in the boreal zone of northern Sweden. During the first year the extent and duration of soil frost was manipulated in the field. The model was successful in reproducing the timing of the soil warming after the snowmelt and frost thaw. A delayed soil warming, into the growing season, severely reduced the transpiration. We demonstrated the potential for considerable overestimation of transpiration by the model if the reduction of the trees' capacity to transpire due to low soil temperatures is not taken into account. We also demonstrated that the accumulated effect of aboveground conditions could be included when simulating the relationship between soil temperature and tree water uptake. This improved the estimated transpiration for the control plot and when soil warming was delayed into the growing season. The study illustrates the need of including antecedent conditions on root growth in the model in order to catch these effects on transpiration. The COUP model is a promising tool for predicting transpiration in high‐latitude stands. Copyright © 2006 John Wiley & Sons, Ltd.     

The influence of macrophyte growth,typical of eutrophic conditions,on river flow velocities and turbulence production

The influence of emergent and submerged macrophytes on flow velocity and turbulence production is demonstrated in a 140 m reach of the River Blackwater in Farnborough, Hampshire, UK. Macrophyte growth occurs in patches and is dominated by Sparganium erectum and Sparganium emersum. In May 2001, patches of S. erectum were already established and occupied 18% of the channel area. The flow adjusted to these (predominantly lateral) patches by being channelled through a narrower cross‐section. The measured velocity profiles showed a logarithmic form, with deviations attributable to topographic control. The channel bed was the main source of turbulence. In September 2001, in‐stream macrophytes occupied 27% of the channel, and overhanging bank vegetation affected 32% of the area. Overall flow resistance, described by Manning's n, showed a threefold increase that could be attributed to the growth of S. emersum in the middle of the channel. Velocity profiles showed different characteristic forms depending on their position relative to plant stems and leaves. The overall velocity field had a three‐dimensional structure. Turbulence intensities were generally higher and turbulence profiles tended to mirror the velocity profiles. Evidence for the generation of coherent eddies was provided by ratios of the root mean square velocities. Spectral analysis identified deviations from the Kolmogorov ?5/3 power law and provided statistical evidence for a spectral short‐cut, indicative of additional turbulence production. This was most marked for the submerged vegetation and, in some instances, the overhanging bank vegetation. The long strap‐like leaves of S. emersum being aligned approximately parallel to the flow and the highly variable velocity field created by the patch arrangement of macrophytes suggest that the dominant mechanism for turbulence production is vortex shedding along shear zones. Wake production around individual stems of S. emersum close to the bed may also be important locally. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical study of enhanced energy dissipation near a seamount

Using a two-dimensional primitive equation model, we examine nonlinear responses of a semidiurnal tidal flow impinging on a seamount with a background Garrett-Munk-like (GM-like) internal wavefield. It is found that horizontally elongated pancake-like structures of high vertical wavenumber near-inertial current shear are created both in the near-field (the region over the slope of the seamount) and far-field (the region over the flat bottom of the ocean). An important distinction is that the high vertical wavenumber near-inertial current shear is amplified only at mid-latitudes in the far-field (owing to a parametric subharmonic instability (PSI)), whereas it is amplified both at mid-and high-latitudes (above the latitude where PSI can occur) in the near-field. In order to clarify the generating mechanism for the strong shear in the near-field, additional numerical experiments are carried out with the GM-like background internal waves removed. The experiments show that the strong shear is also created, indicating that it is not caused by the interaction between the background GM-like internal waves and the semidiurnal internal tides. One possible explanation is proposed for the amplification of high vertical wavenumber near-inertial current shear in the near-field where tide residual flow resulting from tide-topography interaction plays an important role in transferring energy from high-mode internal tides to near-inertial internal waves.     

The dimensions of sand ripples in full-scale oscillatory flows

New large-scale experiments have been carried out in two oscillatory flow tunnels to study ripple regime sand suspension and net sand transport processes in full-scale oscillatory flows. The paper focuses on ripple dimensions and the new data are combined with existing data to make a large dataset of ripple heights and lengths for flows with field-scale amplitudes and periods. A feature of the new experiments is a focus on the effect of flow irregularity. The combined dataset is analysed to examine the range of hydraulic conditions under which oscillatory flow ripples occur, to examine the effects of flow irregularity and ripple three-dimensionality on ripple dimensions and to test and improve existing methods for predicting ripple dimensions.The following are the main conclusions. (1) The highest velocities in a flow time-series play an important role in determining the type of bedform occurring in oscillatory flow. Bedform regime is well characterised by mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth peak velocities in the case of irregular flow. (2) For field-scale flows, sand size is the primary factor determining whether equilibrium ripples will be 2D or 3D. 2D ripples occur when the sand D₅₀ ≥ 0.30 mm and 3D ripples occur when D₅₀ ≤ 0.22 mm (except when the flow orbital diameter is low). (3) Ripple type (2D or 3D) is the same for regular and irregular flows and ripple dimensions produced by equivalent regular and irregular flows follow a similar functional dependence on mobility number, with mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth velocities in the case of irregular flow. For much of the ripple regime, ripple dimensions have weak dependency on mobility number and ripple dimensions are similar for regular and irregular flows with the same flow orbital amplitude. However, differences in ripples produced by equivalent regular and irregular flows become significant at the high mobility end of the ripple regime. (4) Ripple dimensions predicted using the Wiberg and Harris formulae are in poor agreement with measured ripple dimensions from the large-scale experiments. Predictions based on the Mogridge et al. and the Nielsen formulae show better overall agreement with the data but also show systematic differences in cases of 3D ripples and ripples generated by irregular flows. (5) Based on the combined large-scale data, modifications to the Nielsen ripple dimension equations are proposed for the heights and lengths of 2D ripples. The same equations apply to regular and irregular flows, but with mobility number appropriately defined. 3D ripples are generally smaller than 2D ripples and estimates of 3D ripple height and length may be obtained by applying multipliers of 0.55 and 0.73 respectively to the 2D formulae. The proposed modified Nielsen formulae provide an improved fit to the large-scale data, accounting for flow irregularity and ripple three-dimensionality.     

Drag force on a circular cylinder in unseparated laminar high Reynolds number nonlinear random oscillatory flow

The stochastic properties of the drag force maxima on a circular cylinder subjected to nonlinear random waves are investigated. Unseparated laminar high Reynolds number flow is considered. A simplified approach based on second order Stokes waves is presented, including the sum-frequency effect only. It is demonstrated how a drag force formula valid for regular linear waves can be used to find the cumulative distribution function of individual drag force maxima for nonlinear irregular waves. Here the [Wang, 1968] drag force coefficient is used.     

Exchange flow of oil and sea-water in a ruptured submarine pipeline

The rupture of a submarine oil pipeline starts various mechanisms leading to an oil spill. Among these mechanisms the leakage of oil driven by the difference in specific gravities of oil and sea-water is difficult to estimate. A simple mathematical model has been developed and laboratory experiments have been carried out to obtain an insight into the density-driven exchange flow and to determine the leak rate. The mathematical model is predictive and takes account of the effects of friction, inclination of the pipeline, and inertia of the fluid. The experiments were done in a horizontal model pipeline. Theoretical and experimental results are in satisfactory agreement.     

Up-temperature flow of surface-derived fluids in the mid-crust: the role of pre-orogenic burial of hydrated fault rocks

The Walter‐Outalpa shear zone in the southern Curnamona Province of NE South Australia is an example of a shear zone that has undergone intensely focused fluid flow and alteration at mid‐crustal depths. Results from this study have demonstrated that the intense deformation and ductile shear zone reactivation, at amphibolite facies conditions of 534 ± 20 °C and 500 ± 82 MPa, that overprint the Proterozoic Willyama Supergroup occurred during the Delamerian Orogeny (c. 500 Ma) (EPMA monazite ages of 501 ± 16 and 491 ± 19 Ma). This is in contrast to the general belief that the majority of basement deformation and alteration in the southern Curnamona Province occurred during the waning stages of the Olarian Orogeny (c. 1610–1580 Ma). These shear zones contain hydrous mineral assemblages that cut wall rocks that have experienced amphibolite facies metamorphism during the Olarian Orogeny. The shear zone rock volumes have much lower δ¹⁸O values (as low as 1‰) than their unsheared counterparts (7–9‰), and calculated fluid δ¹⁸O values (5–8‰) consistent with a surface‐derived fluid source. Hydrous minerals show a decrease in δD_(H2O) from ?14 to ?22‰, for minerals outside the shear zones, to ?28 to ?40‰, for minerals within the shear zones consistent with a contribution from a meteoric source. It is unclear how near‐surface fluids initially under hydrostatic pressure penetrate into the middle crust where fluid pressures approach lithostatic, and where fluid flow is expected to be dominantly upward because of pressure gradients. We propose a mechanism whereby faulting during basin formation associated with the Adelaidean Rift Complex (c. 700 Ma) created broad hydrous zones containing mineral assemblages in equilibrium with surface waters. These panels of fault rock were subsequently buried to depths where the onset of metamorphism begins to dehydrate the fault rock volumes evolving a low δ¹⁸O fluid that is channelled through shear zones related to Delamerian Orogenic activity.