10.08.2015 Views

Wind Erosion in Western Queensland Australia

Modelling Land Susceptibility to Wind Erosion in Western ... - Ninti One

Modelling Land Susceptibility to Wind Erosion in Western ... - Ninti One

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Chapter 2 – Land Erodibility Controls2 = u *tRfor a surface with roughness elements (2.31) = 2su *tSfor a bare surfacewhere u *tR and u *tS are the threshold friction velocities for rough and smooth surfacesrespectively. The threshold friction velocity ratio (R t ) can then be computed as:u* tS ' s 1Rt = = =(2.32)u * tR( 1)( 1+ )where the factor (1 + βλ) -1/2 accounts for the reduction <strong>in</strong> surface stress by shelter fromroughness elements, and (1 - σλ) -1/2 accounts for the “amplification’ of the stress (τ’s) on theerodible surface over the difference between total stress and that act<strong>in</strong>g on the roughnesselements (τ – τ R ) caused by the fraction of soil surface (σλ) occupied by the roughnesselements. This formulation was then modified to account for the fact that u *t is controlled bythe maximum stress (τ’’s) act<strong>in</strong>g at any po<strong>in</strong>t on the surface, rather than spatially averagedstress (τ’s). The modification was thus designed to account for non-uniformity around nonerodibleroughness elements, and is implemented by the <strong>in</strong>clusion of the m parameter (wherem ≤ 1). Written for the computation of the threshold friction velocity over vegetated surfaces(u *tR ) as an adjustment to the threshold friction velocity of a smooth surface, we have:( 1m)( m )u* tR= u*tS1+(2.33)Key simplifications <strong>in</strong> the model are that: 1) surfaces are only described <strong>in</strong> terms of theroughness density, and the model does not account for differences <strong>in</strong> element shape; and 2)the model does not account for the presence of turbulent air flow close to the surface.Limitations identified <strong>in</strong> the drag partition<strong>in</strong>g scheme revolve around def<strong>in</strong><strong>in</strong>g appropriatevalues for the m parameter, the applicability of the m parameter itself, and the roughnessdensity range for which the model is valid λ ≤ 0.1 (Ok<strong>in</strong> and Gillette, 2004). The dragpartition<strong>in</strong>g scheme presented by Raupach (1992) and Raupach et al. (1993) was revised byYang and Shao (2005) to account for multiple drag partition<strong>in</strong>g (e.g. between roughnesslayers at different heights), and for roughness densities where λ ≥ 0.1.60

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!