Dynamic Model of Mesoscale Eddies     

Mikhail S. Dubovikov 《地球物理与天体物理流体动力学》2013,107(4):311-358

Oceanic mesoscale eddies which are analogs of well known synoptic eddies (cyclones and anticyclones), are studied on the basis of the turbulence model originated by Dubovikov (Dubovikov, M.S., "Dynamical model of turbulent eddies", Int. J. Mod. Phys. B7, 4631-4645 (1993).) and further developed by Canuto and Dubovikov (Canuto, V.M. and Dubovikov, M.S., "A dynamical model for turbulence: I. General formalism", Phys. Fluids 8, 571-586 (1996a) (CD96a); Canuto, V.M. and Dubovikov, M.S., "A dynamical model for turbulence: II. Sheardriven flows", Phys. Fluids 8, 587-598 (1996b) (CD96b); Canuto, V.M., Dubovikov, M.S., Cheng, Y. and Dienstfrey, A., "A dynamical model for turbulence: III. Numerical results", Phys. Fluids 8, 599-613 (1996c)(CD96c); Canuto, V.M., Dubovikov, M.S. and Dienstfrey, A., "A dynamical model for turbulence: IV. Buoyancy-driven flows", Phys. Fluids 9, 2118-2131 (1997a) (CD97a); Canuto, V.M. and Dubovikov, M.S., "A dynamical model for turbulence: V. The effect of rotation", Phys. Fluids 9, 2132-2140 (1997b) (CD97b); Canuto, V.M., Dubovikov, M.S. and Wielaard, D.J., "A dynamical model for turbulence: VI. Two dimensional turbulence", Phys. Fluids 9, 2141-2147 (1997c) (CD97c); Canuto, V.M. and Dubovikov, M.S., "Physical regimes and dimensional structure of rotating turbulence", Phys. Rev. Lett. 78, 666-669 (1997d) (CD97d); Canuto, V.M., Dubovikov, M.S. and Dienstfrey, A., "Turbulent convection in a spectral model", Phys. Rev. Lett. 78, 662-665 (1997e) (CD97e); Canuto, V.M. and Dubovikov, M.S., "A new approach to turbulence", Int. J. Mod. Phys. 12, 3121-3152 (1997f) (CD97f); Canuto, V.M. and Dubovikov, M.S., "Two scaling regimes for rotating Raleigh-Benard convection", Phys. Rev. Letters 78, 281-284, (1998) (CD98); Canuto, V.M. and Dubovikov, M.S., "A dynamical model for turbulence: VII. The five invariants for shear driven flows", Phys. Fluids 11, 659-664 (1999a) (CD99a); Canuto, V.M., Dubovikov, M.S. and Yu, G., "A dynamical model for turbulence: VIII. IR and UV Reynolds stress spectra for shear driven flows", Phys. Fluids 11, 656-677 (1999b) (CD99b); Canuto, V.M., Dubovikov, M.S. and Yu, G., "A dynamical model for turbulence: IX. The Reynolds stress for shear driven flows", Phys. Fluids 11, 678-694 (1999c) (CD99c).). The CD model derives from general principles and does not resort to any free parameters. Yet, it successfully describes a wide variety of quite different turbulent flows. In the present work we apply CD model to the compressible ocean. The model yields mesoscale eddies generated by the baroclinic instability. The latter, in turn, arises from the nonhorizontal orientation of the surfaces of the constant potential density (isopycnals). The obtained dynamic equations for eddy fields reduce to a vertical eigen value problem, an eigen value real part yielding an eddy radius, while an imaginary part - an eddy drift velocity. The size of the eddy is about 3r_d (where r_d is the Rossby deformation radius). The eddy dynamics has the following distinctive features: (1) the large scale potential energy feeds the eddy potential energy (EPE) at scales ～ r_d , (2) from r_d EPE cascades to the smaller scales down to ～ l ₁ determined from the condition that the spectral Rossby number Ro(q) ≡ qU'(q)f^?1 becomes ～ 1 (q is two-dimensional wave number within an isopycnal surface), (3) at scales ～ l ₁ EPE transforms into eddy kinetic energy (EKE) which cascades backwards to the larger scales up to ～ r_d , where it transforms back into EPE, thereby closing the energy flux circulation in a wavenumber space, (4) dissipation of the eddy energy (EE) occurs at scales ～ l ₁ since at those scales the fluctuating component of the vertical shear is maximal and equals to the Brunt-Vaisala frequency. The latter equality is the well known condition for generating the vertical turbulence which dissipates EE. The model enables to determine all turbulence characteristics, including the horizontal (isopycnal) diffusivity κ _h in terms of the large scale mean fields. From the typical values of the latter follow estimates for the parameters of an eddy which agree well with the observational and simulational data: k_h ～ 10³m²s^?1, EKE K ～ 10³m²s^?1, r_d ～ 3 × 10⁴m, l_I ～ 10. In what concerns the bolus velocity, it contains additional terms (as compared to the model of Gent and McWilliams (Gent, P.R. and McWilliams, J.C., "Isopycnal mixing in ocean circulation models", J. Phys. Oceanogr. 20, 150-155 (1990)) which result from the eddy fields advection by a mean velocity ū. Since the latter varies with depth, it is inevitable to differ from the eddy drift velocity that produces a shearing force eroding the eddy coherent structures and, therefore, contributing negatively to EE production. This is in contrast with the positive contribution from the GM term (which is due to the baroclinic instability). In those regions where the disruptive action is stronger, there is no eddy generation.  相似文献   

Turbulent transport of magnetic fields. V. Distribution of magnetic energy in a simple α2-dynamo     

J. H. G. M. Van Geffen  P. Hoyng 《地球物理与天体物理流体动力学》2013,107(1-2):109-123

Abstract

In this paper we analyse the stationary mean energy density tensor T_ij = B_iB_j for the x ²-sphere. This model is one of the simplest possible turbulent dynamos, originally due to Krause and Steenbeck (1967): a conducting sphere of radius R with homogeneous, isotropic and stationary turbulent convection, no differential rotation and negligible resistivity. The stationary solution of the (linear) equation for T_ij is found analytically. Only T_rr , T _θθ and T _φφ are unequal to zero, and we present their dependence on the radial distance r.

The stationary solution depends on two coefficients describing the turbulent state: the diffusion coefficient β≈?u²?τ_c/3 and the vorticity coefficient γ ≈ ?|?×u|²?τ_c/3 where u(r, t) is the turbulent velocity and _c its correlation time. But the solution is independent of the dynamo coefficient α≈??u·?×u?τ_c/3 although α does occur in the equation for T_ij . This result confirms earlier conclusions that helicity is not required for magnetic field generation. In the stationary state, magnetic energy is generated by the vorticity and transported to the boundary, where it escapes at the same rate. The solution presented contains one free parameter that is connected with the distribution of B over spatial scales at the boundary, about which T_ij gives no information. We regard this investigation as a first step towards the analysis of more complicated, solar-type dynamos.  相似文献   

Forces on spheres moving horizontally in a rotating stratified fluid     

P. J. Mason 《地球物理与天体物理流体动力学》2013,107(1):137-154

Abstract

Measurements have been made of the net horizontal force F acting on a sphere moving with horizontal velocity U (Reynolds numbers in the range 10²-10⁴) through a stratified fluid rotating about a vertical axis with uniform angular velocity Ω. In both homogeneous and stratified rotating fluids with small Rossby number R(R = U/Ωa ? 1 where a is the radius of the sphere) the force F is of magnitude 2ΩρUV (where ρ is the density of the fluid and V is the volume of the sphere). In a homogeneous fluid the relative directions of F and U were found to depend on the quantity F = 8Ωa ²/UD (where D is the depth of the fluid in which the object is placed (Mason, 1975)). In a rotating stratified fluid the relative directions of F and U are found to depend on the inverse Froude number k(k = Na/U where N ² = (g/δ)?ρ/?z) provided D > 4aΩ/N. In a homogeneous fluid with F ? 1 the force F is mainly in the U direction (a drag force due to inertial wave radiation) and is ～ ?0.4 |MX 2ΩρUV For F ? 1 a “Taylor column” occurs and the force, in correspondence with theoretical expectations, is ～ - 2Ω |MX UρV In a rotating stratified fluid with N ～2Ω and k ? 1 the force F is mainly in the U direction but is roughly one half of that occurring in the homogeneous situation with F ? 1 (tentatively explained as due to the evanescence of inertia-gravity disturbances). In a rotating stratified fluid with k ? 1 the flow should have no vertical motion (as with F ? 1) and again in correspondence with theoretical expectations the drag is ～ ?2 Ω |MX UρV. In a non-rotating stratified fluid the drag coefficient C _D(C _D = F U/½?ρU ²) was measured in the range k = 0.1 to 10 and had a maximum value ～ 1.2 for k ～ 3.  相似文献   

Improved estimators of correlation and R2 for skewed hydrologic data     

《水文科学杂志》2012,57(1):87-101

ABSTRACT

The coefficient of determination R² and Pearson correlation coefficient ρ = R are standard metrics in hydrology for the evaluation of the goodness of fit between model simulations and observations, and as measures of the degree of dependence of one variable upon another. We show that the standard product moment estimator of ρ, termed r, while well-behaved for bivariate normal data, is upward biased and highly variable for bivariate non-normal data. We introduce three alternative estimators of ρ which are nearly unbiased and exhibit much less variability than r for non-normal data. We also document remarkable upward bias and tremendous increases in variability associated with r using both synthetic data and daily streamflow simulations from 905 calibrated rainfall–runoff models. We show that estimators of ρ = R accounting for skewness are needed for daily streamflow series because they exhibit high variability and skewness compared to, for example, monthly/annual series, where r should perform well.  相似文献   

A laboratory study of baroclinic instability     

J. E. Hart 《地球物理与天体物理流体动力学》2013,107(1):181-209

Abstract

Experimental and theoretical results are presented for a simple system which exhibits baroclinic instability. We consider the motion of two immiscible fluids with densities ρ ₁ and ρ ₂ contained in a cylinder rotating with angular frequency ω. The motion is driven by a contact lid rotating with frequency ω + ω. In this paper ω, ω, 2(ρ ₂ – ρ ₁)/(ρ ₂ + ρ ₁), and the geometry are such that the interface does not intersect the “ground” (e.g. an almost horizontal boundary). The motions are described by two-layer quasi-geostrophic equations which are identical, except perhaps for the presence of interfacial friction and tension, with those used in meteorology and oceanography. For small enough internal Froude number F = 4ω² L ²/(g(Δρ/ρ)H) or small enough Rossby number ? = ω/2ω the flow is steady and axisymmetric, the velocity field in each layer being determined primarily by frictional effects in top, bottom, and interfacial Ekman layers. For certain (F, ?) the flow becomes non-axisymmetric. The transition points for the case where the basic potential vorticity gradient is due to interface slope alone have been carefully measured and are in very good agreement with a linear instability theory which neglects sidewall effects. Some preliminary observations of supercritical motion, which include repeatable amplitude and wavenumber vacillation, are reported.  相似文献   

On the growth of disturbances to forced and dissipated barotropic flows     

Jieping Zou  John Fyfe 《地球物理与天体物理流体动力学》2013,107(1-4):41-56

Abstract

We consider the growth of disturbances to large-scale zonally-asymmetric steady states in a truncated spectral model for forced and dissipated barotropic flow. A variant of the energy method is developed to optimize the instantaneous disturbance energy growth rate. The method involves solving a matrix eigenvalue problem amenable to standard numerical techniques. Two applications are discussed. (1) The global stability of a family of steady states is assessed in terms of the Ekman damping coefficient r. It is shown that monotonic global stability (i.e., every disturbances energy monotonically decays to zero) prevails when r≥r_c . (2) Initially fastest-growing disturbances are constructed in the r<r_c regime. Particular attention is paid to a subregion of the r<r_c regime where initially-growing disturbances exist despite stability with respect to normal modes. Nonlinear time-dependent simulations are performed in order to appraise the time evolution of various disturbances.  相似文献   

Two-and three-dimensional linear convection in a rotating annulus     

G. T. Greed  K. Zhang 《地球物理与天体物理流体动力学》2013,107(1-2):23-34

Abstract

An exceptional case to the model-independent theory of Knobloch (1995) is presented, by investigating a rotating cylindrical annulus of height H and side wall radii r _o and r _i, with non-slip, perfectly thermally conducting side walls and thermally insulating stress-free ends. Radial heating permits the possibility of either two- or three-dimensional convective solutions being the preferred mode. An analytical solution is obtained for the two-dimensional case and a numerical solution for the three-dimensional solution, which is also applied to the two-dimensional solution. It is shown that both two- and three-dimensional solutions can be realized depending on the aspect ratio, γ = H/d, where d = r _o-r _i is the thickness of the annulus, the radii ratio λ = r _i/r _o and the rotation rate of the model. For γ = O(1) and λ = 0.4, the preferred convective solution is three-dimensional when the Taylor number, T < 10² and two-dimensional for T > 10². For small aspect ratios, γ ? 1, the preferred mode is two-dimensional for all rotation rates.  相似文献   

Resistive instabilities in rapidly rotating fluids: Linear theory of the convective modes     

Weijia Kuang  Paul H. Roberts 《地球物理与天体物理流体动力学》2013,107(1-4):129-162

Abstract

This paper analyzes the linear stability of a rapidly-rotating, stratified sheet pinch in a gravitational field, g, perpendicular to the sheet. The sheet pinch is a layer (O ? z ? d) of inviscid, Boussinesq fluid of electrical conductivity σ, magnetic permeability μ, and almost uniform density ρ _o; z is height. The prevailing magnetic field. B _o(z), is horizontal at each z level, but varies in direction with z. The angular velocity, Ω, is vertical and large (Ω ? V_A/d, where V_A = B₀√(μρ₀) is the Alfvén velocity). The Elsasser number, Λ = σB² ₀/2Ωρ₀, measures σ. A (modified) Rayleigh number, R = gβd²/ρ₀V² _A, measures the buoyancy force, where β is the imposed density gradient, antiparallel to g. A Prandtl number, P_K = μσK, measures the diffusivity, k, of density differences.  相似文献   

Treatment of anchor pixels in the METRIC model for improved estimation of sensible and latent heat fluxes     

Ramesh K. Singh  Ayse Irmak 《水文科学杂志》2013,58(5):895-906

Abstract

Reliable estimation of sensible heat flux (H) is important in energy balance models for quantifying evapotranspiration (ET). This study was conducted to evaluate the value of adding the Priestley-Taylor (PT) equation to the METRIC (Mapping Evapotranspiration at high Resolution with Internalized Calibration) model. METRIC was used to estimate energy fluxes for 10 Landsat images from the 2005, 2006 and 2007 crop growing seasons in south-central Nebraska, USA, where each image owing to recent rainfall exhibited high residual moisture content even at the hot pixel. The METRIC model performed satisfactorily for net radiation (R_n ) and soil heat flux (G) estimation with a root mean square error (RMSE) of 52 and 24 W m^-2, respectively. A RMSE of 122 W m^-2 for H indicated the limitation of the METRIC model in estimating H for high residual moisture content of the hot pixel (Alfalfa reference ET fraction, ET _r F > 0.15). The modified METRIC model (wet METRIC or wMETRIC) incorporating the PT equation was applied to calculate H at the anchor pixels (hot and cold) for high residual moisture content of the hot pixel. The α coefficient of the PT equation was locally calibrated using hourly meteorological data from an automatic weather station and R_n and G data from a Bowen ratio flux tower. The mean α coefficient value was 1.14. The wMETRIC model reduced the RMSE of H from 122 to 64 W m^-2 and that of latent heat flux, LE, from 163 to 106 W m^-2. The RMSE of daily ET decreased from 1.7 to 1.1 mm d^-1 with wMETRIC. The results indicate that treatment of anchor pixels for high residual moisture content with the PT approach gives improved estimation of H, LE and daily ET. It is recommended that the wMETRIC model be used for estimating ET if the hot pixel has high residual moisture (i.e. reference ET fraction > 0.15).

Citation Singh, R. K. & Irmak, A. (2011) Treatment of anchor pixels in the METRIC model for improved estimation of sensible and latent heat fluxes. Hydrol. Sci. J. 56(5), 895–906.  相似文献   

Multi‐decadal water table manipulation alters peatland hydraulic structure and moisture retention          下载免费PDF全文

P.A. Moore  P.J. Morris  J.M. Waddington 《水文研究》2015,29(13):2970-2982

A peatland complex disturbed by berm construction in the 1950s was used to examine the long‐term impact of water table (WT) manipulation on peatland hydraulic properties and moisture retention at three adjacent sites with increasing depth to WT (WET, INTermediate reference and DRY). Saturated hydraulic conductivity (K_s) was found to decrease with depth by several orders of magnitude over a depth of 1–1.5 m at all sites. The depth dependence of WT response to rainfall was similar across sites: WT response increased from 1 : 1 at the surface, to 5 : 1 at 50 cm depth. While surface specific yield (S_y) values were similar across all sites, it decreased with depth at a rate of 0.014 cm^?1 in hollows and 0.007 cm^?1 in hummocks. Bulk density (ρ_b) exhibited similar depth‐dependent trends as S_y and explains a high amount of variance (r² > 0.69) in moisture retention across a range of pore water pressures (?15 to ?500 cm H₂O). Because of higher ρ_b, hollow peat had greater moisture retention, where site effects were minimal. However, the estimated residual water content for surface Sphagnum samples, while on average lower in hummocks (0.082 m³ m^?3) versus hollows (0.087 m³ m^?3), increased from WET (0.058 m³ m^?3) to INT (0.088 m³ m^?3) to DRY (0.108 m³ m^?3) which has important implications for moisture stress under conditions of persistent WT drawdown. Given the potential importance of microtopographic succession for altering peatland hydraulic structure, our findings point to the need for a better understanding of what controls the relative height and proportional coverage of hummocks in relation to long‐term disturbance‐response dynamics. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

Blown‐sand transport rate     

Jinren Ni  Zhenshan Li  Cesar Mendoza 《地球表面变化过程与地形》2004,29(1):1-14

  相似文献   

Predicting new snow density in the Italian Alps: A variability analysis based on 10 years of measurements          下载免费PDF全文

Mauro Valt  Nicolas Guyennon  Franco Salerno  Anna B. Petrangeli  Rosamaria Salvatori  Paola Cianfarra  Emanuele Romano 《水文研究》2018,32(20):3174-3187

Despite its strong impact on the time evolution of the snowpack, current estimation of new snow density (ρ_hn) is usually accomplished either by using local empirical techniques or by assuming a constant snow density. Faced with the lack of an estimation model of ρ_hn valid for a wide spatial scale and supported by a suitable number of observations, this study aims to develop simple monthly linear regression models at scale of the entire Italian Alpine chain based on 12,112 snowfall observations at 122 stations, using only air temperature as predictor. Moreover, the remaining variance is investigated in both time and space, also considering some qualitative features of the snowfall events. The daily ρ_hn measurements present a mean value of 115 kg m^?3 (105 and 159 kg m^?3 for dry and wet conditions, respectively). The mean air temperature of the 24 hr preceding the snowfall event has been found to be the best predictor of the ρ_hn, within 31% of uncertainty. The analysis of associated residues allows supporting the idea that the adoption of a more local approach than the one analysed here is not able to substantially increase the predictive capabilities of the model. In fact, the main factor explaining the remaining variance over the air temperature is the wind, but in a complex orography, as mountain regions are, supplying realistic local wind fields is particularly challenging. Therefore, we conclude that using only the daily mean temperature as predictor is a good choice for estimating daily new snow density at scale of Italian Alpine chain, as well as at more regional scale.  相似文献   

Note on the scalar dynamo model     

R. Kaiser 《地球物理与天体物理流体动力学》2013,107(1-2):125-135

Abstract

Bayly (1993) introduced and investigated the equation (? _t + v·▽-η ▽²)S = RS as a scalar analogue of the magnetic induction equation. Here, S(r, t) is a scalar function and the flow field v(r, t) and “stretching” function R(r, t) are given independently. This equation is much easier to handle than the corresponding vector equation and, although not of much relevance to the (vector) kinematic dynamo problem, it helps to study some features of the fast dynamo problem. In this note the scalar equation is considered for linear flow and a harmonic potential as stretching function. The steady equation separates into one-dimensional equations, which can be completely solved and therefore allow one to monitor the behaviour of the spectrum in the limit of vanishing diffusivity. For more general homogeneous flows a scaling argument is given which ensures fast dynamo action for certain powers of the harmonic potential. Our results stress the singular behaviour of eigenfunctions in the limit of vanishing diffusivity and the importance of stagnation points in the flow for fast dynamo action.  相似文献   

Seismic vulnerability evaluation of axially loaded steel built-up laced members Ⅱ： evaluations   总被引：3，自引：2，他引：1  

Kangmin Lee  ;Michel Bruneau 《地震工程与工程振动(英文版)》2008,7(2):125-136

The test results described in Part 1 of this paper （Lee and Bruneau, 2008） on twelve steel built-up laced members （BLMs） subjected to quasi-static loading are analyzed to provide better knowledge on their seismic behavior. Strength capacity of the BLM specimens is correlated with the strength predicted by the AISC LRFD Specifications. Assessments of hysteretic properties such as ductility capacity, energy dissipation capacity, and strength degradation after buckling of the specimen are performed. The compressive strength of BLMs is found to be relatively well predicted by the AISC LRFD Specifications. BLMs with smaller kl/r were ductile but failed to reach the target ductility of 3.0 before starting to fracture, while those with larger kl/r could meet the ductility demand in most cases. The normalized energy dissipation ratio, EC/ET and the normalized compressive strength degradation, Cr″/Cr of BLMs typically decrease as normalized displacements δ/δb,exp increase, and the ratios for specimens with larger kl/r dropped more rapidly than for specimens with smaller kl/r; similar trends were observed for the monolithic braces. The BLMs with a smaller slenderness ratio, kl/r, and width-to-thickness ratio, b/t, experienced a larger number of inelastic cycles than those with larger ratios.  相似文献   

Thermal and magnetic instabilities in a rapidly rotating fluid sphere     

D. R. Fearn 《地球物理与天体物理流体动力学》2013,107(1):103-126

Abstract

A linear analysis is used to study the stability of a rapidly rotating, electrically-conducting, self-gravitating fluid sphere of radius r ₀, containing a uniform distribution of heat sources and under the influence of an azimuthal magnetic field whose strength is proportional to the distance from the rotation axis. The Lorentz force is of a magnitude comparable with that of the Coriolis force and so convective motions are fully three-dimensional, filling the entire sphere. We are primarily interested in the limit where the ratio q of the thermal diffusivity κ to the magnetic diffusivity η is much smaller than unity since this is possibly of the greatest geophysical relevance.

Thermal convection sets in when the temperature gradient exceeds some critical value as measured by the modified Rayleigh number R_c. The critical temperature gradient is smallest (R_c reaches a minimum) when the magnetic field strength parameter Λ ? 1. [R_c and Λ are defined in (2.3).] The instability takes the form of a very slow wave with frequency of order κ/r ² ₀ and its direction of propagation changes from eastward to westward as Λ increases through Λ _c ? 4.

When the fluid is sufficiently stably stratified and when Λ > Λ_m ? 22 a new mode of instability sets in. It is magnetically driven but requires some stratification before the energy stored in the magnetic field can be released. The instability takes the form of an eastward propagating wave with azimuthal wavenumber m = 1.  相似文献   

MAGNETOTELLURIC RESPONSE ON VERTICALLY INHOMOGENEOUS EARTH HAVING CONDUCTIVITY VARYING LINEARLY WITH DEPTH*     

D. KAO 《Geophysical Prospecting》1982,30(6):866-878

Magnetotelluric response is studied for an inhomogeneous medium having conductivity varying linearly with depth as σ(z) =σ₁+αz. For a medium having conductivity increasing linearly with depth, the phase of the impedance approaches 60° at long periods and the apparent resistivity becomes log (ρ_a) = 2 log (1.36/α^1/3) — 1/3 log (T'). The asymptote of log (ρ_a, T'→∞) when plotted against log (T') has a constant gradient —1/3 and has an intercept on the log (T') axis, which equals 6 log (1.36/α^1/3). When a homogeneous layer with a moderate thickness overlies an inhomogeneous half-space, this layer does not affect the asymptote, but it affects the cut-off period and pushes this toward the long period direction. For a medium having conductivity decreasing linearly with depth, the impedance is equivalent to that of a Cagniard two-layer model; the intercept period related to the thickness is T'₀=σ₁(h₂/2)². Homogeneous multilayer approximations to an inhomogeneous layer are also investigated, and it is shown that the fit to the model variation depends on the number of layers and the layer parameters chosen.  相似文献   

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

Derege Tsegaye Meshesha  Atsushi Tsunekawa  Mitsuru Tsubo  Nigussie Haregeweyn  Enyew Adgo 《水文科学杂志》2013,58(12):2203-2215

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.