Summery Intra-Tidal Variations of Suspended Sediment Transportation–Topographical Response and Dynamical Mechanism in the Aoshan Bay and Surrounding Area,Shandong Peninsula

As the second largest bay in Qingdao,the Aoshan Bay and its adjacent sea area play an important role in aquaculture development and urban land and sea coordination for the eastern gulf type of city in the Qingdao Blue Silicon Valley Core Area(QBSVCA).Based on in-situ sedimentary dynamical observation and previous research results,the thermohaline structure,the transportation of suspended sediment and its mechanism,and the coastal geomorphic response were elaborated and analyzed in detail in this paper.The result indicated that the thermohaline and density distribution have obvious intra-tidal characteristics in the QBSVCA and the adjacent waters of the islands,during summer neap tide stage.The development of the bottom high suspended sediment concentration(SSC)layer was slightly enhanced in flood slack at each of the four stations.Suspended sediment transportation near the QBSVCA is related closely with the vertical mixing-stratification mechanism.Combined with previous research results,this study once again showed that submarine topography and the grain size of sea bed sediments would respond to hydrodynamic forces.The medians of the bottom E and D50 in the Aoshan Bay were the highest,followed by those in the Daguan Island and Xiaoguan Island,and the data in the Laoshan Bay were the lowest.This showed that the capacity of suspended sediment transportation in the bottom water layer of the Aoshan Bay was stronger than that of the adjacent sea area.The re-suspension and migration of fine sediments lead to the strong coarsening of sediments in this area.

Tracing suspended and bed sediments during high and low water periods using geochemical characteristics-Case study:Vazrood watershed,northern Iran

Complete and comprehensive information about sediment dynamic and identification of hotspots of sediment production and transport are necessary for understanding the erosion processes and increasing the efficiency of soil and water conservation practices.Numerous studies used the sediment fingerprint techniques to investigate the contribution of different sources in suspended and bed sediment yield of the watersheds.However,the contribution of various land use/land covers in suspended and bed sediment yield for the great Caspian Sea basin is in an aura of ambiguity and the present study was conducted to gather information about an important part of this area in northern Iran,where rangelands are located upstream of Hyrcanian forests and dense agricultural lands are located downstream.The surface soil of different land use/land covers including forest,rangeland,agriculture and streambank lands were sampled in 30 points.Suspended and bed sediments were sampled in the watershed outlet in two high and low water periods.Geochemical characteristics of soil and sediment samples containing 59 elements were measured using ICP-OES GBC Integra.The reliable and suitable tracers from 59 elements were then selected using Range test,Kruskal-Wallis and Discriminant Function Analysis,respectively,in FingerPro package of R software.The results showed that for suspended sediment,streambank and rangeland had the highest contributions of 86.2%and 47.5%,respectively,in two high and low water periods.For bed sediment,in two high and low water periods,rangeland and streambank had the highest contributions of 73.8%and 84.4%,respectively.Land use change and especially human activities such as agriculture,road construction and development of residential areas along the main river riparian zone has led to a significant increase in suspended and bed sediments.

A new formula of recovery factor for non-equilibrium transport of graded suspended sediment in the Middle Yangtze River

Suspended sediment concentrations in the Middle Yangtze River(MYR)reduced greatly after the Three Gorges Project operation,causing the composition of bed material to coarsen continuously.However,little is known about the non-equilibrium transport of graded suspended sediment owing to different bed material compositions(BMCs)along the MYR,and it is necessary to determine the magnitude of recovery factor.Using the Markov stochastic process in conjunction with the hiding-exposure effect of non-uniform bed-material,a new formula is proposed for calculating the recovery factor including the effect of different BMCs,and it is incorporated into the non-equilibrium transport equation to simulate the recovery processes of suspended load in both sand-gravel bed and sand bed reaches of the MYR.The results show that:(i)the recovery rate of graded sediment concentrations at Zhicheng was slower than that at Shashi during the period 2003-2007;(ii)the mean recovery factors of the coarse,medium,and fine sediment fractions in the ZhichengShashi reach were 0.152,0.0012,and 0.0005,respectively,and the coarse sediment recovered up to the maximum sediment concentration of 0.138 kg/m3over a distance of 15 km;and(iii)the results of the new formula that can consider the effect of bed material composition are in general agreement with the field observations,and the spatial and temporal delay effects are inversely related to particle size and BMC.Consequently,the BMC effect on the nonequilibrium sediment transport in different reaches of the MYR needs to be considered for higher simulation accuracy.

Study on the Sediment Transport Flux and Mechanism in the Bohai Strait at the Tidal and Monthly Scales in Summer

Based on the data of tidal currents and suspended sediment concentrations observed synchronously at 11 stations in the Bohai Strait lasting for 25 hours,the temporal and spatial variations of currents and suspended sediment concentrations in the Bohai Strait in summer were analyzed.The Study preliminarily discussed the transport mechanism,transport trend and transport flux of suspended sediments in summer,using flux-mechanism decomposition method and numerical simulation.The suspended sediment transport was mainly controlled by advection and next influenced by vertical net circulation,while resuspension is relatively weak in Bohai Strait.The single-width fluxes of investigation stations varied from 3.8 to 89.1 gm^(−1) s^(−1),with the maximum value in Miaodao Strait.The suspended sediment transport trends in Laotieshan channel along the vertical section are obviously distinct.The waters mainly flow out of the Bohai Sea in surface layer,while into the Bohai Sea in bottom layer.However,the transport trends of other channels in the centre and south are consistent vertically.The sediments in the Bohai Strait follows the transport pattern of moving outward from the south and inward from the north in summer,i.e.,the sediments are carried out of the Bohai Sea through the Laotieshan channel,while into the Bohai Sea through other channels.And the outflow flux exceeds the inflow flux in August with the net water flux of 1.4×10^(10)m^(3),basically same as the deliveries of the rivers into the Bohai Sea.Moreover,the suspended sediment flux is 0.33 Mt under the action of tidal residual currents in the Yellow Sea in August.

Influence of asymmetric tidal mixing on sediment dynamics in a partially mixed estuary

To investigate the influence of asymmetric tidal mixing(ATM) on sediment dynamics in tidal estuaries, we developed a vertically one-dimensional idealized analytical model, in which the M_2 tidal flow, residual flow and suspended sediment concentration(SSC) are described. Model solutions are obtained in terms of tidallyaveraged, and tidally-varying components(M_2 and M_4) of both hydrodynamics and sediment dynamics. The effect of ATM was considered with a time-varying eddy viscosity and time-varying eddy diffusivity of SSC. For the first time, an analytical solution for SSC variation driven by varying diffusivity could be derived. The model was applied to York River Estuary, where higher(or lower) eddy diffusivity was observed during flood(or ebb) in a previous study. The model results agreed well with the observation in both hydrodynamics and sediment dynamics. The vertical sediment distribution under the influence of ATM was analyzed in terms of the phase lag of the M_2 component of SSC relative to tidal flow. The phase lag increases significantly in estuaries with typical ATM(higher diffusivity during flood and lower diffusivity during ebb) for the case of seaward-directed net bottom shear stress(e.g., strong river discharge). In contrary, the phase lag is reduced by ATM, if the tidally-averaged bottom shear stress is landward(e.g., strong horizontal density gradient). The dynamics of sediment transport was analyzed as a function of ATM phase lag to identify the time of highest sediment diffusivity, as well as a function of the residual flow, to evaluate the relative importance of seaward and landward residual flows. In estuaries with relative strong fresh water discharge or weak tidal forcing(in case of flood season or neap tide), the near bottom SSC could be higher during ebb than during flood, since the bottom shear stress is higher during ebb due to seaward residual flow. However, landward net sediment transport can be expected in these estuaries in case of a typical ATM, because higher diffusivity causes higher SSC and landward transport during the flood period, while both SSC and seaward transport could be lower during ebb. On the contrary, seaward sediment transport can be expected in estuaries with landward tidally mean bottom shear stress in case of a reverse ATM,where sediment diffusivity is higher during the ebb.

Numerical Simulation of Bed Load and Suspended Load Sediment Transport Using Well-Balanced Numerical Schemes

Sediment transport can be modelled using hydrodynamic models based on shallow water equations coupled with the sediment concentration conservation equation and the bed con-servation equation.The complete system of equations is made up of the energy balance law and the Exner equations.The numerical solution for this complete system is done in a seg-regated manner.First,the hyperbolic part of the system of balance laws is solved using a finite volume scheme.Three ways to compute the numerical flux have been considered,the Q-scheme of van Leer,the HLLCS approximate Riemann solver,and the last one takes into account the presence of non-conservative products in the model.The discretisation of the source terms is carried out according to the numerical flux chosen.In the second stage,the bed conservation equation is solved by using the approximation computed for the system of balance laws.The numerical schemes have been validated making comparisons between the obtained numerical results and the experimental data for some physical experiments.The numerical results show a good agreement with the experimental data.

A High-Order Conservative Semi-Lagrangian Solver for 3D Free Surface Flows with Sediment Transport on Voronoi Meshes

In this paper,we present a conservative semi-Lagrangian scheme designed for the numeri-cal solution of 3D hydrostatic free surface flows involving sediment transport on unstruc-tured Voronoi meshes.A high-order reconstruction procedure is employed for obtaining a piecewise polynomial representation of the velocity field and sediment concentration within each control volume.This is subsequently exploited for the numerical integration of the Lagrangian trajectories needed for the discretization of the nonlinear convective and viscous terms.The presented method is fully conservative by construction,since the transported quantity or the vector field is integrated for each cell over the deformed vol-ume obtained at the foot of the characteristics that arises from all the vertexes defining the computational element.The semi-Lagrangian approach allows the numerical scheme to be unconditionally stable for what concerns the advection part of the governing equations.Furthermore,a semi-implicit discretization permits to relax the time step restriction due to the acoustic impedance,hence yielding a stability condition which depends only on the explicit discretization of the viscous terms.A decoupled approach is then employed for the hydrostatic fluid solver and the transport of suspended sediment,which is assumed to be passive.The accuracy and the robustness of the resulting conservative semi-Lagrangian scheme are assessed through a suite of test cases and compared against the analytical solu-tion whenever is known.The new numerical scheme can reach up to fourth order of accu-racy on general orthogonal meshes composed by Voronoi polygons.

Effects of shore-normal coastal structure on medium-to long-term embayed shoreline evolution

Based on high-tide shoreline data extracted from 87 Landsat satellite images from 1986 to 2019 as well as using the linear regression rate and performing a Mann-Kendall(M–K)trend test,this study analyzes the linear characteristics and nonlinear behavior of the medium-to long-term shoreline evolution of Jinghai Bay,eastern Guangdong Province.In particular,shoreline rotation caused by a shore-normal coastal structure is emphasized.The results show that the overall shoreline evolution over the past 30 years is characterized by erosion on the southwest beach,with an average erosion rate of 3.1 m/a,and significant accretion on the northeast beach,with an average accretion rate of 5.6 m/a.Results of the M–K trend test indicate that significant shoreline changes occurred in early 2006,which can be attributed to shore-normal engineering.Prior to that engineering construction,the shorelines are slightly eroded,where the average erosion rate is 0.7 m/a.However,after shore-normal engineering is performed,the shoreline is characterized by significant erosion(3.2 m/a)on the southwest beach and significant accretion(8.5 m/a)on the northeast beach,thus indicating that the shore-normal engineering at the updrift headland contributes to clockwise shoreline rotation.Further analysis shows that the clockwise shoreline rotation is promoted not only by longshore sediment transport processes from southwest to northeast,but also by cross-shore sediment transport processes.These findings are crucial for beach erosion risk management,coastal disaster zoning,regional sediment budget assessments,and further observations and predictions of beach morphodynamics.

Sediment Transport in the Nearshore Area of Phoenix Island

Based on the measured data, suspended sediment concentration, surface sediment grain size, current and waves, the sediment transport mechanisms and pathways in the Phoenix Island area were analyzed using methods of flux decomposition and Grain Size Trend Analysis(GSTA). The results show that net suspended sediment is mainly transported by average current, Stokes drift, and gravitational circulation. The transport direction of suspended sediment is varying and basically following the direction of residual tidal currents. Surface sediment transport pathways are primarily parallel to the coastline along with two convergent centers. Waves and longshore currents have a significant influence on sediment transport, but the influence is limited due to a steep and deep underwater bank. Tidal current is the main controlling factor for sediment transport, especially in the deep water area. Neither suspended nor surface sediment is transported towards the southwest. The South Shandong Coastal Current(SSCC) has little effect on sediment transport processes in the nearshore area of Phoenix Island.

Scale effects of eroded sediment transport in Wujiang River Basin, Guizhou Province, China

In recent years, research on spatial scale and scale transformation of eroded sediment transport has become a forefront field in current soil erosion research, but there are very few studies on the scale effect problem in Karst regions of China. Here we quantitatively extracted five main factors influencing soil erosion, namely rainfall erosivity, soil erodibility, vegetative cover and management, soil and water conservation, and slope length and steepness. Regression relations were built between these factors and also the sediment transport modulus and drainage area, so as to initially analyze and discuss scale effects on sediment transport in the Wujiang River Basin(WRB). The size and extent of soil erosion influencing factors in the WRB were gauged from: Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model(ASTER GDEM), precipitation data, land use, soil type and Normalized Difference Vegetation Index(NDVI) data from Global Inventory Modeling and Mapping Studies(GIMMS) or Advanced Very High Resolution Radiometer(AVHRR), and observed data from hydrometric stations. We find that scaling effects exist between the sediment transport modulus and the drainage area. Scaling effects are expressed after logarithmic transformation by a quadratic function regression relationship where the sediment transport modulus increases before decreasing, alongside changes in the drainage area. Among the five factors influencing soil erosion, slope length and steepness increases first and then decreases, alongside changes in the drainage area, and are the main factors determining the relationship between sediment transport modulus and drainage area. To eliminate the influence of scale effects on our results, we mapped the sediment yield modulus of the entire WRB, adopting a 1 000 km^2 standard area with a smaller fitting error for all sub-basins, and using the common Kriging interpolation method.

The controlling factors of high suspended sediment concentration in the intertidal flat off the Huanghe River Estuary

The Huanghe River(Yellow River)is known by its high suspended sediment concentration(SSC)in its river mouth tidal flat.However,the factors controlling the high SSC over there are not well understood.Therefore,we conducted 7-d hydrodynamic observations(water depth,wave height,and current velocity)and SSC measurements on the tidal flat off the Huanghe River Mouth.The data shows that in most of time,under the calm sea condition,the SSC ranges 0.1–3.5 g/L,and sediment discharge from the river is the main source.However,when hydrodynamics are enhanced in a tidal cycle and large-scale erosion occurs on the seafloor,resuspended sediment becomes the main source,and the SSC in the water column reaches 17.3 g/L.We find the suspended sediment flux is mainly controlled by the tidal current and Stokes drift,while the wave-induced shear stress could also affect the variation of suspended sediment flux.During the observation period,when sea under calm-rippled conditions,the current-induced resuspended sediment concentration(RSC)was greater than the wave-induced RSC.In contrast,in smooth-wavelet sea conditions,the wave-induced RSC was greater than the current-induced RSC,for instance,a single wave event was found to cause 11.8 cm seabed erosion within 6 h.This study reveals different controlling factors for the high SSC near a river-influenced tidal flat,and helps us get a better understanding of a delta's depositional and erosional mechanisms.

Sediment Characteristics, Sources, and Transport Patterns in Kompong Som Bay, Cambodia: Indications from Grain Size and Heavy Minerals

In this paper, we analyzed the grain size and heavy mineral compositions of 52 surface sediment samples collected from the Kompong Som Bay of Cambodia and the adjacent rivers to depict the marine sedimentary environments and transport processes. Heavy minerals in sediments are dominated by authigenic pyrite, siderite, and tourmaline, with average percentages of 36.52%, 29.02%, and 10.94%, respectively. Two provinces can be divided according to the spatial similarity of minerals. The sediments from Province I, covered by silt grains in the northern bay, are characterized by autogenic pyrite, indicating a weakly reducing environment;whereas in Province II, covered by sand grains in the southern bay, the siderite-tourmaline-authigenic pyrite-zircon-hornblende assemblage occurs, indicating a mild reducing environment and locally oxidizing environment. Most of the sediments in the Kompong Som Bay are introduced from the Preak Piphot River and Srae Ambel River, except that some of them in the south areas come from coastal erosion. Generally, the sediments are difficult to be transported because of the low sediment loads entering the sea and weak hydrodynamic conditions. However, they are transported from the north to the south during the tide ebbing when the hydrodynamic force is much stronger. The sediment distribution and transport patterns are controlled by many factors, including submarine topography, hydrodynamic conditions, the southwest monsoon, land contours, and sediment supply.

Sediment generation and soil mound denudation in areas of high-density tree throw along a river valley in the Jura Mountains, Switzerland

A preliminary field-based investigation was undertaken in a small(<10 km^(2))river valley located in the mountainous Jura region of northwest Switzerland.The aims of the work were to assess sediment generation and annual sediment transport rates by tree throw on forested hillslopes,and to document surface hydrology characteristics on four fresh soil mounds associated with recent tree throws over a 24-day monitoring period.For the soil mounds,average sediment recovery ranged from 7.7-28.2 g(dry weight),equivalent to a suspended sediment concentration of 145.2-327.8 g L^(-1),and runoff coefficients ranged from 1.0%-4.2%.Based on a soil bulk density value of 1,044 kg m^(-3),upslope runoff generation areas were denuded by an average 0.14 mm by the end of the 24-day monitoring period,representing an erosion rate equivalent to 2.1 mm yr^(-1).A ca.50 cm high soil mound could therefore feasibly persist for around 200-250 years.For tree throw work,the dimensions of 215 individual tree throws were measured and their locations mapped in 12 separate locations along the river valley representing a cumulative area equivalent to 5.3 ha(av.density,43 per ha).Tree throws generated a total of 20.1 m^(3) of fine-sediment(<2 mm diameter),or the equivalent of 3.8×10^(-4) m^(3) m^(-2).The process of tree throw was originally attributed to two extreme weather events that occurred in west and central Europe in late December 1999.Taking the 18-year period since both storms,this represents an annual sediment transport rate of 2.7×10^(-5) m^(3) m^(-1) yr^(-1).Exploring the relationship with wind on fall direction,65.5%of tree throws(143)generally fell in a downslope direction irrespective of hillslope aspect on which they were located.This infers that individual storms may not have been responsible for the majority of tree throws,but instead,could be associated with root failure.Given the high density of tree throws and their relative maturity(average age 41 years),we hypothesise that once trees attain a certain age in this river valley,their physiognomy(i.e.height,mass and centre of gravity)compromises their ability to remain securely anchored.We tentatively attribute this possibility to the presence of bedrock close to the surface,and to the shallow soil profile overlaying steep hillslopes.

Linking Fish Habitat Modelling and Sediment Transport in Running Waters

The assessment of ecological status for running waters is one of the major issues within an integrated river basin management and plays a key role with respect to the implementation of the European Water Frame- work Directive (WFD).One of the tools supporting the development of sustainable river management is physi- cal habitat modeling,e.g.,for fish,because fish population are one of the most important indicators for the e- colngical integrity of rivers.Within physical habitat models hydromorphological ...

Impact of utility-scale solar photovoltaic array on the aeolian sediment transport in Hobq Desert, China

Deserts are ideal places to develop ground-mounted large-scale solar photovoltaic (PV) powerstation. Unfortunately, solar energy production, operation, and maintenance are affected bygeomorphological changes caused by surface erosion that may occur after the construction of the solar PVpower station. In order to avoid damage to a solar PV power station in sandy areas, it is necessary toinvestigate the characteristics of wind-sand movement under the interference of solar PV array. The studywas undertaken by measuring sediment transport of different wind directions above shifting dunes andthree observation sites around the PV panels in the Hobq Desert, China. The results showed that the twoparameterexponential function provides better fit for the measured flux density profiles to the near-surfaceof solar PV array. However, the saltation height of sand particles changes with the intersection anglebetween the solar PV array and wind direction exceed 45°. The sediment transport rate above shifting duneswas always the greatest, while that around the test PV panels varied accordingly to the wind direction.Moreover, the aeolian sediment transport on the solar PV array was significantly affected by wind direction.The value of sand inhibition rate ranged from 35.46% to 88.51% at different wind directions. When theintersection angle exceeds 45°, the mean value of sediment transport rate above the solar PV array reducesto 82.58% compared with the shifting dunes. The results of our study expand our understanding of theformation and evolution of aeolian geomorphology at the solar PV footprint. This will facilitate the designand control engineering plans for solar PV array in sandy areas that operate according to the wind regime.

Bed Load Transport Rate in Scouring and Armoring Process of Non-uniform Sediment River Bed

Flume experiments were carried out to study bed load transport rate during rive bed scouring and ar- moring.A theoretical differential equation linking the transport rate to the probability of incipient motion of non-uniform sediment is solved.The transport rate is shown to decrease exponentially with time,according to the theory,which is in good agreement with the experiment data.

A three-dimensional,wave-current coupled,sediment transport model for POM

In the high-energy environment of coastal seas and estuaries,strong sediment resuspension/ deposition events are driven by surface waves,tides,winds and buoyancy driven currents.In recent years,A POM based three-dimensional,wave-current coupled,sediment transport model has been developed by the University of New South Wales.This paper presents several examples of the model applications to study sediment dynamics in the environments where forcings such as waves,tides,and winds are equally important to affect sediment fluxes and distributions.Firstly,the sediment transport model coupled to the Yellow Sea general circulation model and a third generation wave model SWAN was implemented in the Yellow Sea to study the dynamics of the sediment transport and resuspension in the northern Jiangsu shoal-wate(rNJSW).The sediment distributions and fluxes and their inter-annual variability were studied by realistic numerical simulations.The study found that the surface waves played a dominant role over the tides to form the turbidity maxima along the muddy coast of NJSW. Secondly,the sediment transport model was used to explore the effect of suspended sediment-induced stratification in the bottom boundary laye(rBBL).The model uses a re-parameterized bottom drag coefficient Cd that incorporates a linear stability function of flux Richardson number Rf.The study has shown that the sediment induced stratification in the BBL reduces the vertical eddy viscosity and bottom shear stress in comparison with the model prediction in a neutrally stratified BBL.In response to these apparent reductions,the tidal current shear is increased and sediments are abnormally concentrated within a thin wall layer that is overlain by a thicker layer with much smaller concentration.The formation of this fluid-mud layer near the seabed has led to a significant reduction in the total sediment transport.This study contributes to the understanding of formations of tidal flats along the coasts of turbid seas and estuaries.

Sediment Transport in Rivers with Overbank Flow

Some concepts related to sediment transport in rivers with overbank flow are described.Following a description of the physical processes that are involved when a river inundates its floodplains,some simple com- putational methods are presented which permit the depth-averaged velocity and boundary shear stress to be pre- dicted within a cross section of variable,but prismatic shape.The methoda account for the strong transverse shear in velocity that occurs when the stage is just above bankfull,as well as ...

The Dynamics of Suspended Particulate Matter in the Oued Nachef(Tafna Basin,Algeria):Typology of Flood Events and Contribution to Sediment Transport

This article focuses on the study of sediment transport during flood events in the Oued Nachef watershed that feeds the Mefrouche dam.To understand the sediment dynamics in this watershed,ANRH data on instantaneous water discharges and the respective concentrations of suspended particulate matter were used.This enabled the selection of some of the largest flood events over a 24-year period in order to establish the log-log relationships between sediment load(concentration and flux)and water discharge.However,the discharge-concentration relationships revealed hysteresis phenomena that enabled a flood typology to be established and classified into seven categories,thus showing very different transfer dynamics in relation to flood events.The results showed that Category 6 floods presenting hysteresis in the form of a figure of eight exported almost 44%of the suspended particulate matter load while representing just 29%of the flow discharge.

Assessment of Runoff, Sediment Yields and Nutrient Loss Using the Swat Model in Upper Indus Basin of Pakistan

The main objective of this study is to understand the runoff, sediment yield and water quality of the Upper Indus River Basin of Pakistan. To achieve this goal, specific objectives have been met which include, setup of a hydrological model using Soil and Water Assessment Tool (SWAT) then calibration and validation of the hydrological model using river discharges and in the end investigating the performance of the hydrological model by SWAT. This research will have great impacts on socio-economic conditions of Pakistan because study of upper Indus River basin is imperative to provide data needed for its management, and to warrant that it is sustainable to support the increasing population and conservational flows. A set of programmable mapping components MapWindow Geographic Information System (GIS) was used which is an open source GIS based mapping application. It is SWAT used spatially distributed information on elevation, land use, slope and soil. The program Sequential Uncertainty Fitting ver.2 (SUFI-2) in a combination of uncertainty analysis and calibration of outputs was used in SWAT-CUP. SWAT model used input data, which have climate information to obtain results. The observed climate data of temperature gauges and rain gauge were used as input in the SWAT model;the calibration results for three discharge stations were produced. The initial P-factor value was satisfactory but more iteration to attempt narrow uncertainty band with improving goal function, resulted in small percentage of observed data within uncertainty band. A warm up period of three years (1979-1982) was used for simulation of SWAT model. The model was calibrated for selected three catchments for the period 1982-2000 and validated for period 2001-2010. Results are quite comparable with the observed flows.