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.     

Fully nonlinear numerical wave tank (NWT) simulations and wave run-up prediction around 3-D structures

A finite-difference scheme and a modified marker-and-cell (MAC) algorithm have been developed to investigate the interactions of fully nonlinear waves with two- or three-dimensional structures of arbitrary shape. The Navier–Stokes (NS) and continuity equations are solved in the computational domain and the boundary values are updated at each time step by the finite-difference time-marching scheme in the framework of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique developed for two fluid layers.To demonstrate the capability and accuracy of the present method, the numerical simulation of backstep flows with free-surface, and the numerical tests of the MDF technique with limit functions are conducted. The 3D program was then applied to nonlinear wave interactions with conical gravity platforms of circular and octagonal cross-sections. The numerical prediction of maximum wave run-up on arctic structures is compared with the prediction of the Shore Protection Manual (SPM) method and those of linear and second-order diffraction analyses based on potential theory and boundary element method (BEM). Through this comparison, the effects of non-linearity and viscosity on wave loading and run-up are discussed.     

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.     

Statistical evaluation of methods for the calculation of static formation temperatures in geothermal and oil wells using an extension of the error propagation theory

The feasibility of using the more sophisticated weighted least-squares (WLS) model, as opposed to the traditional ordinary least-squares (OLS), in linear regressions of BHT data to estimate the static formation temperatures (SFT) was investigated. The most commonly used analytical methods (line-source; spherical and radial heat flow; and cylindrical heat source) were applied. Error propagation equations were derived to calculate errors in the time function of each method. These errors were combined with the BHT measurement errors to compute weighting factors for applying the WLS. Intercept uncertainties were estimated for all regressions using sets of synthetic and actual borehole logs taken from geothermal and oil applications. SFT computed with the spherical and radial heat flow method were generally greater than those from the other two methods.     

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.