Modeling Soil Erosion by Water on Agricultural Land in Cenei, Timiş County, Romania
Total Page:16
File Type:pdf, Size:1020Kb
Advances in Environmental Development, Geomatics Engineering and Tourism Modeling soil erosion by water on agricultural land in Cenei, Timiş County, Romania GEORGE NARCIS PELEA, IOANA ALINA COSTESCU, ERIKA BEILICCI, TEODOR EUGEN MAN, ROBERT BEILICCI Faculty of Civil Engineering, Hydrotehnical Engineering Department University Politehnica of Timişoara G. Enescu Street No. 1A 300022 Timişoara, Timiş County ROMANIA [email protected] http://www.upt.ro Abstract: This paper aims to present soil erosion by water on agricultural land with low slope up to 5%. In the study were analyzedthe agricultural land of Cenei, in Timiş County, Romania, on part of drainage systems Răuţi – Sânmihaiu German and Checea – Jimbolia, and on the course section of river BegaVeche. The model analyzes soil erosion by water and sediment transport on the main drainage channels and the course section of river BegaVeche. Key-Words:soil erosion by water, sediment transport, drainage system, earth channels,hydrographic area, low slope erosion. 1 Introduction In order to evaluate the soil erosion by water on the Through the water cycle in nature appears the studied surface due to the small slope of phenomenon of soil erosionby water on agricultural surface,which is up to 5%., for the analysis was used land. This phenomenon affects the condition of the a modelling program -WEPP (Water Erosion drainage systems by clogging the channels section. Prediction Project) [5] Due to this problem it is brought to discussion the WEPP is software for the prediction of erosion possibility of analysis of soil erosionby water on processes in watersheds, developed by USDA agricultural land with small slopes. Forest Service, Agricultural Research Service, The model analyzes agricultural land of Ceneiwhich Natural Resources Conservation Service, is located in Timiş County, on western part of Department of Interior's Bureau of Land Romania[1]. Management and Geological Survey from USA, Agricultural land in the territorial area ofCenei now days, the most used model in numerous overlaps with parts of drainage systems Răuţi – countries. WEPP has been under a permanent Sânmihaiu German and Checea – Jimbolia development process since 1985. administered by the National Agency for Land The WEPP model may be used in both hillslope and Reclamation, Timis Branch. Both drainage systems watershed applications. The model is a distributed were designed and built in the '70s, and consists in parameter, continuous simulation and erosion main drainage channels, secondary channels and prediction, implemented as a set of computer tertiary channels. Themain channels are named CPE programs for personal computers. and they transport water to the drainage pumping The hillslope component of the WEPP erosion stations, SP Uivar and SP Cenei, located on the left model requires a minimum of four input data files to and right bank of the river BegaVechewhich flows run: climate file, slope file, soil file and beside Cenei.[2][3] plant/management file. The watershed component requires a minimum of seven input data files: each hillslope information file, structure file, slope file, soil file, management 2 Problem Formulation file, climate file and channel file. The sediment formation consist of processes in WEPP considering the hillslope consist in numerous which soil particles are detached and transported on parallel rills; the surface erosion occurs on interrill soil surface by external dynamic agents surfaces and the dislocated soil particles are (precipitation, wind) and deposited at different distances from the detachment place. [4] ISBN: 978-960-474-385-8 20 Advances in Environmental Development, Geomatics Engineering and Tourism transported downhill by rill flow (rill erosion is also Ie – effective rainfall intensity (m/s) considered in calculus). σir – interrill runoff rate (m/s) WEPP produces many different kinds of output, in SDRRR – sediment delivery ratio (function various quantities, depending upon the user’s needs. of random roughness, side slope, particle size The basic output contains the runoff and erosion distribution) summary information, which may be produced on a Fnozzle – irrigation adjustment factor (if is storm by storm, monthly, annual or average annual case) basis. The time – integrated estimates of runoff, Rs – rill spacing (m) erosion, sediment delivery and sediment enrichment w – rill width. are contained in this output, as well as the spatial distribution of erosion on the hillslope (fig.1). Rill erosion is divided into two parts, detachment and deposition. Rill detachment is predicted when the flow shear stress exceeds the critical shear stress of the soil and when the sediment load is below the calculated sediment transport capacity. For the calculus of Drthe following relation is used: G Dr Dc1−= (3) Tc where: Dc – detachment capacity by rill flow (kg/s.m2) Tc – sediment transport capacity (kg/s.m). Dc Kr ( τ−τ= cf ) (4) where: Kr– rill erodibility parameter τ f – flow shear stress acting on soil (Pa) τ c – critical shear stress of soil (Pa). Deposition in a rill will be predicted when the Fig.1 WEPP Model sediment load is greater than the sediment transport capacity. For its calculus, we use the following The model has two basic components: relation: βV Erosion component f Dr = ()Tc − G (5) It is based on the continuity equation: q dG where: Vf – effective fall velocity for the sediment Dr += D (1) dx i (m/s) q – flow discharge per unit width (m2/s) where: G – sediment load (kg/s.m) β – rainfall turbulence coefficient. x – distance downslope (m) Dr – rill erosion rate (+ for detachment, - for The shear stress factor τ for the rill flow is deposition) calculated with the help of the following equation: 2 Di – interrill sediment delivery (kg/s.m ). f s αγ=τ For the calculus of D following relation is used: R sin (6) i ft RS 2 2 Di KiadjI σ= ire SDR RR Fnozzle (2) where: γ – specific weight of water (kg/m .s ) w R – hydraulic radius (m) α – average slope angle for uniform where: Kiadj – interrillerodibility segment ISBN: 978-960-474-385-8 21 Advances in Environmental Development, Geomatics Engineering and Tourism fs/ft – ratio of friction factor for the soil to S – average capillary suction at the wetting total rill friction factor (accounts for shear stress front acting on surface cover like residue etc.) IMD – initial moisture deficit Ks – saturated hydraulic conductivity of For the calculus of hydraulic radius, a rectangular soil. rill is assumed and the flow depth in the rill is calculated using the Darcy – Weisbach equation: Depresional storage is related by following relation: 8gRS 2 f = (7) Sd =0.112 ⋅ rr +3.1 ⋅ rr −1.2 ⋅ rr ⋅S0 (11) V2 where: rr – random roughness (m) where: g – gravity (m/s2) S0 – slope of flow surface. R – hydraulic radius (m) S – average slope (%) Peak discharge (Pr) of surface flow is calculated V – flow velocity (m/s). using a kinematic wave model: The transport capacity is calculated by applying the ∂h ∂q + =i − f (12) full Yalin (1963) equation to the end of an overland ∂t ∂x flow element (OFE) or hillslope. where: i – rainfall intensity (m/s) T =k τ 2/3 f – infiltration intensity (m/s) c t f (8) h – depth of flow (m) 2 0,5 2 0,5 q – discharge per unit with of plane (m /s where: k – transport coefficient (m s /kg ), t t – time (s) depends on hill slope x – distance (m). τ – hydraulic shear stress acting on the soil f (Pa). q= C ⋅ S 5.0 h 5.1 0 (13) Hydrology component where: C – Chezy coefficient The hydrology component of the WEPP model S0 – slope value. produces four inputs for the erosion component: peak runoff rate, effective runoff duration, effective Effective duration of rainfall (te) (rainfall intensity > rainfall intensity, effective rainfall duration. infiltration intensity) is calculated using the Infiltration is calculated using the Green - Ampt and following relation: Mein – Larson (GAML) equations: For F < Fs: 5.1/1 L S⋅ IMD te = (14) F = for i> K 5.0 5.0 s s CS⋅()i − f f= i and i / Ks −1 (9) 0 no calculation of Fs for i≤ Ks where: L – length of the overland flow surface (m). For F ≥ Fs: The effective runoff duration (tr) is calculated as: S⋅ IMD Vt f= fp and fp =Ks1 + (10) tr = (15) F Pr 3 where: f – infiltration rate where: Vt – total runoff volume (m ) 3 fp – infiltration capacity Pr – peak runoff rate (m /s). i – rainfall intensity F – cumulative infiltration volume The effective rainfall intensity (Ie) (rainfall intensity Fs – cumulative infiltration volume required > infiltration intensity) is calculated as: to cause surface saturation ISBN: 978-960-474-385-8 22 Advances in Environmental Development, Geomatics Engineering and Tourism ∫ Idt To create de model were introduced the land surface Ie = (16) with slops to the drainage channels (fig.4). For te analyzing the surface was essential the correlation between field data consisting in land reclamation where: I – break point rainfall intensity (m/s). works and land surveying. 3 Problem Solution Cenei locality is located in the western part of Romania in Timiş County, near the border with Serbia (fig.2). The location is in a plain area with relatively small slopes up to 5%.[6] Fig. 4 WEPP model with surfaces and slopes Fig. 2 Cenei location and the analyzed area For obtaining relevant data into the program were The analyzed surface overlaps with parts of the entered two types of crops: cereals and corn (fig.5). drainange systems Răuţi – Sânmihaiu German and Essential at this stage was entering the stratification Checea – Jimbolia (fig.3). of terrain and soil type for each surface. Fig. 3 WEPP model containing the analyzed surface Fig. 5 WEPP model with crop management ISBN: 978-960-474-385-8 23 Advances in Environmental Development, Geomatics Engineering and Tourism The climate data file was created usingdata from t/ha/yrwith a discharge volume of 11509206 Jimbolia meteorological station or where it was m3/yr.The precipitation volume in contributing area necessary fromTimisoara meteorological station is about 60331420 m3/yr.