Throughfall partitioning by trees

Although we know that rainfall interception (the rain caught, stored, and evaporated from aboveground vegetative surfaces and ground litter) is affected by rain and throughfall drop size, what was unknown until now is the relative proportion of each throughfall type (free throughfall, splash throughfall, canopy drip) beneath coniferous and broadleaved trees. Based on a multinational data set of >120 million throughfall drops, we found that the type, number, and volume of throughfall drops are different between coniferous and broadleaved tree species, leaf states, and timing within rain events. Compared with leafed broadleaved trees, conifers had a lower percentage of canopy drip (51% vs. 69% with respect to total throughfall volume) and slightly smaller diameter splash throughfall and canopy drip. Canopy drip from leafless broadleaved trees consisted of fewer and smaller diameter drops (D_{50_DR}, 50th cumulative drop volume percentile for canopy drip, of 2.24 mm) than leafed broadleaved trees (D_{50_DR} of 4.32 mm). Canopy drip was much larger in diameter under woody drip points (D_{50_DR} of 5.92 mm) than leafed broadleaved trees. Based on throughfall volume, the percentage of canopy drip was significantly different between conifers, leafed broadleaved trees, leafless broadleaved trees, and woody surface drip points (p ranged from <0.001 to 0.005). These findings are partly attributable to differences in canopy structure and plant surface characteristics between plant functional types and canopy state (leaf, leafless), among other factors. Hence, our results demonstrating the importance of drop‐size‐dependent partitioning between coniferous and broadleaved tree species could be useful to those requiring more detailed information on throughfall fluxes to the forest floor.     

Application of an optical disdrometer to characterize simulated rainfall and measure drop-size distribution

ABSTRACT

Knowledge of rainfall characteristics such as drop-size distribution is essential for the development of erosion-mitigation strategies and models. This research used an optical disdrometer to elucidate the relationships between raindrop-size distribution, median volume drop diameter (D₅₀), kinetic energy and radar reflectivity (dBz) of simulated rainfall of different intensities. The D₅₀ values were higher for the simulated rain than for natural rain at almost all rainfall intensities, perhaps due to variations in rainfall types and the turbulence in natural rain that breaks up large drops. The kinetic energy ranged from 26.67 to 5955.51 J m^?2 h ^?1, while the median volume drop diameter (D₅₀) was in the range 1.94–7.25 mm, for intensities between 1.5 and 202.6 mm h^?1. The relationship between radar reflectivity (Z) and the intensity (R) of the simulated rain was best described by a power law function (Z = aR^b), with a and b coefficients in the ranges 162–706 and 0.94–2.46, respectively, throughout the range of rainfall intensities (1.5–202.6 mm h^?1).     

Multifractal comparison of the outputs of two optical disdrometers

ABSTRACT

In this paper a universal multifractals comparison of the outputs of two types of collocated optical disdrometers installed on the roof of the Ecole des Ponts ParisTech is performed. A Campbell Scientific PWS100 which analyses the light scattered by the hydrometeors and an OTT Parsivel² which analyses the portion of occluded light are deployed. Both devices provide a binned distribution of drops according to their size and velocity. Various fields are studied across a range of scales: rain rate (R), liquid water content (ρ), polarimetric weather radar quantities such the horizontal reflectivity (Z_h) and the specific differential phase (K_dp), and drop size distribution (DSD) parameters such as the total drop concentration (N_t) and the mass-weighted diameter (D_m). For both devices, good scaling is retrieved on the whole range of available scales (2?h–30?s), except for the DSD parameters for which the scaling only holds down to few minutes. For R, the universal multifractal parameters are found to equal 1.5 and 0.2 for α and C₁, respectively. Results are interpreted with the help of the classical Z_h–R and R–K_dp radar relations.

Editor D. Koutsoyiannis; Associate editor E. Volpi     

The coalescence of water drops I. A theoretical model of approaching drops

The approach of two water drops in the absence of air flow around them is theoretically investigated. By assuming deformation criteria it is possible to solve the equation of motion of the drops under the influence of a variety of forces. These forces include the viscous force exerted by the air between the two deformed surfaces, the London-Van Der Waals forces and the force of gravity. It is found that the viscous forces dominate over the whole distance of the interaction. The equations have analytical solutions when a head-on approach is considered and when the deformation of the drops is assumed constant during the interaction. The equations were solved numerically for other deformation criteria and for non head-on approaches.The results of the present model are used in the following paper to compute the coalescence efficiencies of water drops. The model is primarily applicable to situations in which the large drop is stationary and the small one approaches it from below. However, it could also be used for interaction between freely falling drops as long as their relative velocities exceed about 13 cm/sec.Appendix: List of symbols C constant of the motion - D distance between the deformed surfaces of the drops - D _o initial value ofD - D _m the value at which the viscous force is maximum - D _N normalized distance - D _s the distance at which the velocity of approach vanishes - F _c centrifugal force - F _g force due to gravity - F _N normalized viscous force - F _LV force due to London-Van der Waals effect - F _R radial component of the force - F _V viscous force - F _t tangential component of the force - g acceleration due to gravity - M _L mass of large drop - m _s mass of small drop - p ratio of radii of interacting drops - R radius of an arbitrary drop - r distance between the centers of mass of the two drops - R _D radius of deformation - R _L radius of larger drop - R _s radius of smaller drop - t time - u defined in equation 20 — has the meaning of kinetic energy - v relative velocity of the deformed surfaces - v ₀ initial value ofv - V ₀ initial relative velocity of the centers of the drops - V _c critical impact velocity - V _i impact velocity - V _N ,v _n normalized velocity - V _t tangential component of the velocity - W _i velocity of the small drop at infinity for it to reach the pointD ₀ at velocityV ₀ - x instantaneous impact distance -  average critical impact distance for coalescence - x ₀ initial value of the impact distance - x _c critical impact distance for coalescence - coefficient of deformation - _i impact angle according toWhelpdale andList (1971) - coefficient of deformation - viscosity - surface tension - F _s sum of forces acting on the small drop - F _L sum of forces acting on the large drop - time constant - _R Rayleigh's oscillation period On sabbatical leave (1976–77) from the Department of Geophysics and Planetary Sciences, Tel Aviv University, Ramat Aviv, Israel.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Characterization of differential throughfall drop size distributions beneath European beech and Norway spruce

Forest canopies present irregular surfaces that alter both the quantity and spatiotemporal variability of precipitation inputs. The drop size distribution (DSD) of rainfall varies with rainfall event characteristics and is altered substantially by the forest stand properties. Yet, the influence of two major European tree species, European beech (Fagus sylvatica L.) and Norway spruce (Picea abies (L.) H. Karst), on throughfall DSD is largely unknown. In order to assess the impact of these two species with differing canopy structures on throughfall DSD, two optical disdrometers, one above and one below the canopy of each European beech and Norway spruce, measured DSD of both incident rainfall and throughfall over 2 months at a 10‐s resolution. Fractions of different throughfall categories were analysed for single‐precipitation events of different intensities. While penetrating the canopies, clear shifts in drop size and temporal distributions of incoming rainfall were observed. Beech and spruce, however, had different DSD, behaved differently in their effect on diameter volume percentiles as well as width of drop spectrum. The maximum drop sizes under beech were higher than under spruce. The mean ± standard deviation of the median volume drops size (D50) over all rain events was 2.7 ± 0.28 mm for beech and 0.80 ± 0.04 mm for spruce, respectively. In general, there was a high‐DSD variability within events indicating varying amounts of the different throughfall fractions. These findings help to better understand the effects of different tree species on rainfall partitioning processes and small‐scale variations in subcanopy rainfall inputs, thereby demonstrating the need for further research in high‐resolution spatial and temporal properties of rainfall and throughfall.     

The coalescence of water drops II. The coalescence problem and coalescence efficiency

By the use of the model of approaching drops (Arbel and Levin, 1977) the coalescence efficiencies of drops are computed. It is found that for interactions of drops at their terminal velocities the coalescence depends both on the size of the large drop and on the size ratio of the interacting drops in agreement with the experimental results of Whelpdale and List (1971) and Levin and Machnes (1977).The results were found to be sensitive to the assumption of the drops deformation and to the critical separation distance. This distance is defined as the distance at which the drops begin to merge. The variations of the coalescence efficiency with these parameters is discussed.Appendix: List of symbols D distance between the deformed surfaces of the drops - D _o initial value ofD - D _s stop distance, the distance at which the impact velocity vanishes - D _c critical coalescence distance - E collection efficiency - E ₁ collision efficiency - E ₂ coalescence efficiency - E _2R coalescence efficiency for collisions with stationary targets - F _c centrifugal force - p ratio of the radii of the interacting drops - r _o initial distance between drops' centers - R _L radius of larger drop - R _s radius of smaller drop - R _D radius of deformation - v approach velocity of two deformed surfaces - v _o initial value ofv - V _i impact velocity (given negative sign when drops approach each other) - V _c critical impact velocity - W _i velocity of the smaller drop at infinity for it to reachD _o with velocityv _o - x _i impact distance, the distance between the trajectories of the two drops - x _c critical impact distance for coalescence -  average critical impact distance for coalescence - X _c critical impact distance for collisions - coefficient of deformation given in equation 1 - _i impact angle defined byWhelpdale andList (1971) given also inArbel andLevin (1977) - coefficient of deformation given in equation 2 - viscosity of air - _i impact angle used inArbel andLevin (1977) and here - _c critical angle for coalescence - average critical angle for coalescence On sabbatical leave (1976–77) from the Department of Geophysics and Planetary Sciences, Tel Aviv University, Ramat Aviv, Israel.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Simulation of the size distribution and erosivity of raindrops and throughfall drops

This paper presents a model that simulates the size distribution and erosivity of raindrops and throughfall drops. It utilizes existing models of rainfall drop size distribution and fall velocity and combines them with newly collated evidence of throughfall drop size distributions. A sensitivity analysis reveals that the model is sensitive to parameters that are easily measured or estimated: rainfall intensity, the mean volume drop diameter of the intercepted throughfall, canopy cover, and canopy height. The results of the model may be used at two levels. Firstly, to calculate specifically the size and fall velocity of individual drops, parameters that are needed in studies examining the response of soil surfaces to forces applied by rainfall. Secondly, to produce erosivity indices, based on rainfall intensity but which take account of the effects of a vegetation canopy. The paper shows that while the kinetic energy of rainfall (E(0), J mm^?1 m^?2) may be calculated from an equation of the familiar form: the kinetic energy of throughfall under any canopy may be calculated by combining this equation with another that relates the energy of drops under a 100 per cent canopy cover (E(100)) and the canopy height: .     

Rain properties controlling soil splash detachment

Laboratory experiments have been conducted to study the effects of various rain properties on sand detachment resulting from raindrop impact. Splash cups were exposed to simulated rainfall intensities ranging between 10 and 140 mm h⁻¹. The detached sand was collected and weighed whereas rain intensity, equivalent drop diameter and fall velocity of raindrops were measured with an optical spectro‐pluviometer (OSP). The properties of the simulated rain (i.e. median volume diameter and kinetic energy) were compared with those observed in natural conditions. Statistical analyses have been undertaken in order to evaluate which rain variable best predicts the mass of sand detached. Linear and non‐linear correlations between the mass of detached sediment and the product of drop size (d) by drop velocity (v), i.e. D^αV^β, with values of α varying between 1 and 6 and β between 0 and 3, have been computed. The results indicate that the coefficient of determination (R²) for α ranging between 3 and 5 and β lower or equal to 2 are satisfying. Although kinetic energy (D³V²) described splash detachment relatively well, the product of momentum by drop diameter (D⁴V) was slightly superior in describing splash detachment. Therefore, the momentum multiplied by the drop diameter is recommended as the best rain variable to describe splash detachment. Copyright © 2000 John Wiley & Sons, Ltd.     

Effects of collision-induced breakup on drop size distributions in steady state rainshafts

New equations and techniques for dealing with drop breakups are developed and applied to the modelling of the evolution of raindrop spectra in rainshafts. Breakup experiments byMcTaggart-Cowan andList (1975) served as data base.No matter what the original size distribution, the spectrum evolution will always lead to a Marshall-Palmer type equilibrium di tributionN=N ₀e^–D, with =constant andN ₀ proportional to the rainfall rateR. (D stands for raindrop diameter.) ForR29 mm h^–1 and an original Marshall-Palmer distribution, the required fall height to reach equilibrium is 2 km.The equilibrium distributions are characterized by linear relationships betweenR, the radar reflectivity factorZ, the liquid water content LWC and theN ₀ of the Marshall-Palmer distribution. Possible explanations for the discrepancy with observations are given.The fact that the all-water processes cannot produce drops withD2.5 mm (as confirmed by observations) leads to the conclusion that observed large raindrops withD5 mm represent melted hailstones and have not yet reached an equilibrium distribution. These latter conclusions were reached within the original assumption of videspread, steady state precipitation.     

Characteristics of cloud drop spectra in tropical warm clouds

Characteristics of cloud drop spectra were studied using 400 samples obtained from 120 warm cumulus clouds formed during the summer monsoon season.The total concentration of cloud drops (N _T) varied from 384 to 884 cm^–3 and the maximum concentration was observed in the layer below the cloud-top. The width of the drop spectrum was broader in the cloud-base region and in the region below the cloud-top. The spectrum was multimodal at all levels except in the cloud-top region where it was unimodal. The concentration of drops with diameter greater than 50 m (N _L) varied from 0.0 to 0.674 cm^–3.N _L was larger in the cloud-base region.N _L decreased with height up to the middle level and thereafter showed an increase. In the cloud-top region no large drops were present. The computed values of the liquid water varied between 0.132 and 0.536 g m^–3 and the mean volume diameter (MVD) varied between 8.1 and 12.0 m. The LWC and MVD showed a decrease with height except in the middle region of the cloud where the values were higher than the adjacent levels. The dispersion of the cloud drops was lower (0.65) in the cloud-top region and higher (1.01) in the cloud-base region.The observed cloud microphysical characteristics were attributed to vertical mixing in clouds induced by the cloud-top gravity oscillations (buoyancy oscillations) generated by the intensification of turbulent eddies due to the buoyant production of energy by the microscale-fractional-condensation (MFC) in turbulent eddies.     

Non-axisymmetric α2Ω-dynamo waves in thin stellar shells

Linear α²Ω-dynamo waves are investigated in a thin turbulent, differentially rotating convective stellar shell. A simplified one-dimensional model is considered and an asymptotic solution constructed based on the small aspect ratio of the shell. In a previous paper Griffiths et al. (Griffiths, G.L., Bassom, A.P., Soward, A.M. and Kuzanyan, K.M., Nonlinear α²Ω-dynamo waves in stellar shells, Geophys. Astrophys. Fluid Dynam., 2001, 94, 85–133) considered the modulation of dynamo waves, linked to a latitudinal-dependent local α-effect and radial gradient of the zonal shear flow. These effects are measured at latitude θ by the magnetic Reynolds numbers R _α f(θ) and R _Ω g(θ). The modulated Parker wave, which propagates towards the equator, is localised at some mid-latitude θ_p under a Gaussian envelope. In this article, we include the influence of a latitudinal-dependent zonal flow possessing angular velocity Ω_*(θ) and consider the possibility of non-axisymmetric dynamo waves with azimuthal wave number m. We find that the critical dynamo number D _c?=?R _α R _Ω is minimised by axisymmetric modes in the αΩ-limit (Rα→0). On the other hand, when Rα?≠?0 there may exist a band of wave numbers 0?m?m ^? for which the non-axisymmetric modes have a smaller D _c than in the axisymmetric case. Here m ^? is regarded as a continuous function of R _α with the property m?→0 as R _α→0 and the band is only non-empty when m??>1, which happens for sufficiently large R _α. The preference for non-axisymmetric modes is possible because the wind-up of the non-axisymmetric structures can be compensated by phase mixing inherent to the α²Ω-dynamo. For parameter values resembling solar conditions, the Parker wave of maximum dynamo activity at latitude θ_p not only propagates equatorwards but also westwards relative to the local angular velocity Ω_*(θ _p ). Since the critical dynamo number D _c?=?R _α R _Ω is O (1) for small R _α, the condition m ^??>?1 for non-axisymmetric mode preference imposes an upper limit on the size of |dΩ_*/dθ|.     

Drop size distribution and kinetic energy load of rainfall events in the highlands of the Central Rift Valley,Ethiopia

Abstract

Knowledge of rainfall characteristics is important for estimating soil erosion in arid areas. We determined basic rainfall characteristics (raindrop size distribution, intensity and kinetic energy), evaluated the erosivity of rainfall events, and established a relationship between rainfall intensity I and volume-specific kinetic energy KE_vol for the Central Rift Valley area of the Ethiopian highlands. We collected raindrops on dyed filter paper and calculated KE_vol and erosivity values for each rainfall event. For most rainfall intensities the median volume drop diameter (D₅₀) was higher than expected, or reported in most studies. Rainfall intensity in the region was not high, with 8% of rain events exceeding 30 mm h^-1. We calculated soil erosion from storm energy and maximum 30-min intensity for soils of different erodibility under conditions of fallow (unprotected soil), steep slope (about 9%) and no cover and management practice on the surface, and determined that 3 MJ mm ha^-1 h^-1 is the threshold erosivity, while erosivity of >7 MJ mm ha^-1 h^-1 could cause substantial erosion in all soil types in the area.

Editor Z.W. Kundzewicz; Associate Editor Q. Zhang     

Heat Waves in the South Moravian Region During the Period 1961-1995

Heat waves (periods of extremely hot summer weather) in the region of south Moravia are in the focus of this study. The introduced definition consists of three requirements imposed on the period that is considered a heat wave: at least three days with T _MAX 30.0°C must be observed; the mean T _MAX over the whole period is at least 30.0°C; and T _MAX must not drop below 25.0°C. To compare the severity of the individual heat waves, various characteristics (duration, number of tropical days, peak temperature, cumulative temperature excess, precipitation amount) are examined. The heat wave index HWI is defined to express the severity of heat waves in the most comprehensive way. An extraordinary heat wave occurred in July and August 1994; it lasted more than a month at several stations, while the duration of a typical heat wave is only 4 - 7 days. The extremely long unbroken period of tropical days, and even of days with T _MAX 32.0°C, represents the most distinct feature of the severe 1994 heat wave. With regard to heat wave characteristics, the summer temperature exceptionality of the early 1990s is indubitable.     

Pan evaporation modelling and changing attribution analysis on the Tibetan Plateau (1970–2012)

Pan evaporation (E_p) is an important indicator of water and energy and the decline of E_p has been reported in many regions over the last decades. The climate and E_p are dependent on each other. In this study, the temporal trends of E_p and main E_p drivers, namely mean air temperature (T_a), wind speed (u), global solar radiation (R_s), net long‐wave radiation(R_nl) and vapour pressure deficit (D) from 1970 to 2012, were calculated on the basis of 26 meteorological stations on the Tibetan Plateau. The arithmetic average of E_p from 26 stations decreased with the rate of ?11.91 mm a^?2; the trends of R_s, R_nl, T_a, u and D were ?1.434 w m^?2 decade^?1, 0.2511 w m^?2 decade^?1, 0.3590°C decade^?1, ?0.2376 m s^?1 decade^?1 and 9.523 Pa decade^?1, respectively. The diffuse irradiance is an essential parameter to model E_p and quantify the contribution of climatic factors to changing E_p. 60 724 observations of R_s and diffuse solar irradiance (R_d) from seven of the 26 stations were used to develop the correlation between the diffuse fraction (R_d/R_s), and the clearness index (R_s/R_o). On the basis of the estimation of the diffuse component of R_s and climatic data, we modified the PenPan model to estimate Chinese micro‐pan evaporation (E_p) and assess the attribution of E_p dynamics using partial derivatives. The results showed that there was a good agreement between the observed and calculated daily E_p values. The observed decrease in E_p was mostly due to declining wind speed (?13.7 mm a^?2) with some contributions from decreasing solar irradiance (?3.1 mm a^?2); and the increase of temperature had a large positive effect (4.55 mm a^?2) in total whilst the increase of R_nl had insignificant effect (0.35 mm a^?2) on E_p rates. The change of E_p is the net result of all the climatic variables. Copyright © 2014 John Wiley & Sons, Ltd.     

Variation of lunar and solar geomagnetic tides with sunspot and geomagnetic activity

Summary Solar and lunar geomagnetic tides inH at Alibag have been determined by spectral analysis of discrete Fourier transforms following the method of Black and the well-known Chapman-Miller method. The seasonal variation inL ₂(H) is opposite to that inL ₂(D) with maximum in thed season and minimum in thej season. In bothH andD the enhancement due to sunspot activity is larger in lunar tide than in solar tide. Surprisingly, the enhancement due to magnetic activity is greater inL ₂(H) than inS ₁(H), while the contrary is true for declination. It is inferred that there is a local time component of the storm time variation contrary to the view expressed by Green and Malin. The enhancements in amplitudesL ₂ andS ₁ inH andD, due to sunspot activity and due to magnetic activity, have been separated. The results show that the amplitude at zero sunspot number increases with magnetic activity in all the four parameters, while the enhancement due to sunspot activity at different levels of magnetic activity decreases with increase ofK _p. But if bothK _p andR are increasing, whenK _p increases enhancement due toR decreases and whenR increases enhancement due toK _p decreases.     

On the frontal condition for finite amplitude α2Ω-dynamo wave trains in stellar shells

Abstract

In a previous paper, Bassom et al. (Proc. R. Soc. Lond. A, 455, 1443–1481, 1999) (BKS) investigated finite amplitude αΩ-dynamo wave trains in a thin turbulent, differentially rotating convective stellar shell; nonlinearity arose from α-quenching. There asymptotic solutions were developed based upon the small aspect ratio ε of the shell. Specifically, as a consequence of a prescribed latitudinally dependent α-effect and zonal shear flow, the wave trains have smooth amplitude modulation but are terminated abruptly across a front at some high latitude θ_F. Generally, the linear WKB-solution ahead of the front is characterised by the vanishing of the complex group velocity at a nearby point θ_f; this is essentially the Dee–Langer criterion, which determines both the wave frequency and front location.

Recently, Griffiths et al. (Geophys. Astrophys. Fluid Dynam. 94, 85–133, 2001) (GBSK) obtained solutions to the α²Ω-extension of the model by application of the Dee—Langer criterion. Its justification depends on the linear solution in a narrow layer ahead of the front on the short O(θ_f—θ_F) length scale; here conventional WKB-theory, used to describe the solution elsewhere, is inadequate because of mode coalescence. This becomes a highly sensitive issue, when considering the transition from the linear solution, which occurs when the dynamo number D takes its critical value D _c corresponding to the onset of kinematic dynamo action, to the fully nonlinear solutions, for which the Dee—Langer criterion pertains.

In this paper we investigate the nature of the narrow layer for α²Ω-dynamos in the limit of relatively small but finite α-effect Reynolds numbers R _α, explicitly ε^½ ? R ² _α ? 1. Though there is a multiplicity of solutions, our results show that the space occupied by the corresponding wave train is generally maximised by a solution with θ_f—θ_F small; such solutions are preferred as evinced by numerical simulations. This feature justifies the application by GBSK of the Dee—Langer criterion for all D down to the minimum D _min that the condition admits. Significantly, the frontal solutions are subcritical in the sense that |D _min| ≤ |D _c|; equality occurs as the α-effect Reynolds number tends to zero. We demonstrate that the critical linear solution is not connected by any parameter track to the preferred nonlinear solution associated with D _min. By implication, a complicated bifurcation sequence is required to make the connection between the linear and nonlinear states. This feature is in stark contrast to the corresponding results for αΩ-dynamos obtained by BKS valid in the limit R _{2 α} ? ε^½, which, though exhibiting a weak subcriticality, showed that the connection follows a clearly identifiable nonbifurcating track.     

Transport and dispersal of organic debris (peat blocks) in upland fluvial systems

This paper assesses the mechanisms and pathways by which peat blocks are eroded and transported in upland fluvial systems. Observations and experiments from the north Pennines (UK) have been carried out on two contrasting river systems. Mapping of peat block distributions and appraisal of reach‐based sediment budgets clearly demonstrates that macro‐size peat is an important stream load component. In small streams block sizes can approximate the channel width and much of the peat is transported overbank. Local ‘peat jams’ and associated mineral deposition may provide an important component of channel storage. In larger systems peat blocks rapidly move down‐channel and undergo frequent exchanges between bed and bank storage. Results of peat block tracing using painted blocks indicate that once submerged, blocks of all sizes are easily transported and blocks break down rapidly by abrasion. Vegetation and bars play an important role in trapping mobile peat. In smaller streams large block transport is limited by channel jams. Smaller blocks are transported overbank but exhibit little evidence of downstream fining. In larger rivers peat blocks are more actively sorted and show downstream reduction in size from source. A simple model relating peat block diameter (D_p) to average flow depth (d) suggests three limiting transport conditions: flotation (D_p < d), rolling (d < D_p > d/2) and deposition (D_p > d/2). Experiments demonstrate that peat block transport occurs largely by rolling and floating and the transport mechanism is probably controlled by relative flow depth (d/D_p ratio). Transport velocity varies with transport mechanism (rolling is the slowest mode) and transport lengths increase as flow depth increases. Abrasion rates vary with the transport mechanism. Rolling produces greater abrasion rates and more rounded blocks. Abrasion rates vary from 0 to 10 g m^?1 for blocks ranging in mass from 10 to 6000 g. Copyright © 2001 John Wiley & Sons, Ltd.     

First Results of DSD Measurement by Videodistrometer in the Czech Republic in 1998-1999

The first results of the almost one year drop size distribution (DSD) measurement in the Czech Republic are summarised in this study. The ESA-ESTEC 2D videodistrometer was used to measure the rain drop parameters. The average DSD is shown to be of the gamma type. One minute DSDs were evaluated to test the accuracy of analytical DSD models. Parameters of gamma distribution and exponential distribution functions were evaluated for the whole data set as well as for the various rain rate intervals. Regression technique and the method of moments were applied to estimate the parameters of DSD. It is shown that the parameter value strongly depends on the method of computation as well as on the rain type. Its average value is about 0.59 for the average (smooth) one minute DSD while an average value of un-smoothed DSD is 11.0 (moment method) or 5.4 (regression technique). The Joss's shape parameter and the Tokay-Short's parameter CS estimating roughly the rain type are also discussed (if CS>1, the event should be convective). The tendency of increasing numerical value of the CS parameter with the increasing rain rate was observed (the DSDs were distributed into classes respecting the rain rate value) and thus the idea of the convectivity occurrence bounded with the higher CS parameter value was supported. The study also compares the parameters of the average DSD with the averages of parameter values of all 4 183 one minute DSDs.     

Seasonal and annual throughfall and stemflow in Andean temperate rainforests

The partitioning of gross rainfall into throughfall, stemflow, and interception loss and their relationships with forest structure was studied for a period of four years (October 2002–September 2006) and two years (October 2005–September 2007) in seven experimental catchments of temperate rainforest ecosystems located in the Andes of south‐central Chile (39°37′S, 600–925 m a.s.l.). The amount of throughfall, stemflow, and interception loss was correlated with forest structure characteristics such as basal area, canopy cover, mean quadratic diameter (MQD), and tree species characteristics in evergreen and deciduous forests. Annual rainfall ranged from 4061 to 5308 mm at 815 m a.s.l. and from 3453 to 4660 mm at 714 m a.s.l. Throughfall ranged from 64 to 89% of gross rainfall. Stemflow contributed 0·3–3·4% of net precipitation. Interception losses ranged from 11 to 36% of gross rainfall and depended on the amount of rainfall and characteristics as well as on forest structure, particularly the MQD. For evergreen forests, strong correlations were found between stemflow per tree and tree characteristics such as diameter at breast height (R² = 0·92, P < 0·01) and crown projection area (R² = 0·65, P < 0·01). Stemflow per tree was also significantly correlated with epiphyte cover of trunks in the old‐growth evergreen forests (R² = 0·29, P < 0·05). The difference in the proportion of throughfall and interception loss among stands was significant only during winter. The reported relationships between rainfall partitioning and forest structure and composition provide valuable information for management practices, which aimed at producing other ecosystem services in addition to timber in native rainforests of southern Chile. Copyright © 2010 John Wiley & Sons, Ltd.     

Mixing efficiency in stably-stratified decaying turbulence

Abstract

The behavior of the flux Richardson number R _f, as a function of the overall Richardson number Ri ₀, was investigated for a stably stratified, grid-generated, turbulent flow evolving in a closed-loop water channel. The turbulent dissipation rate ε, the buoyancy or vertical mass flux p wbar; and the rms density fluctuation ρ′ were obtained from simultaneous single-point measurements of the horizontal and vertical velocity components and density fluctuations. From these, R _f and Ri ₀ were calculated at each point in the spatially evolving flow. The resulting curves of R _f vs. Ri ₀ exhibit the full range of behavior found in the very different case studied by Linden (1980). The length scale arguments of Gibson (1980) and Stillinger et al. (1983b) provide an underlying mechanism which successfully accounts for the shape of the R _f vs. Ri ₀ curve.