• Title/Summary/Keyword: Discharge water

Search Result 3,352, Processing Time 0.039 seconds

Migration Patterns of Brachymystax lenok tsinlingensis Using Radio Tags in the Upper Part of the Nakdong River (Radio tag을 이용한 낙동강 상류에 서식하는 열목어의 이동양상)

  • Yoon, Ju-Duk;Jang, Min-Ho
    • Korean Journal of Ecology and Environment
    • /
    • v.42 no.1
    • /
    • pp.58-66
    • /
    • 2009
  • The telemetry of eight adult manchurian trouts (Brachymystax lenok tsinlingensis) in the upper part of the Nakdong River, which is the southern limit of distribution of manchurian trout on the Korean peninsula, was used to examine migration patterns and evaluate characteristics of over-wintering and the spawning season between December, 2007 and May, 2008. Based on the tracking data, the tagged fish showed a limited migration between adjacent pools, moving only up to $8.6m\;day^{-1}$ during the winter season (December to February). Hydraulic conditions of over-wintering pool areas were, ca. 1m depth, slow moving surface water with areas of sand and gravel. The migration of tagged individuals was successful, moving up to $96.2m\;day^{-1}$ during the spawning season. Two tagged individuals (BL4, BL6) exhibited upstream migration, whereas others showed downstream movements. The timing of upstream migration of the two individuals was consistent with an increasing phases of water level and discharge. The fishes migrating toward the down stream moved to the wide pool areas downstream, where they spent the summer season for the growth.

Distributional characteristics of phytoplankton and nutrient limitation during spring season in Jinhae Bay (춘계 진해만에서 식물플랑크톤 증식과 제한영양염 분포특성)

  • Son, Moonho;Kim, Dongseon;Baek, Seung Ho
    • Journal of the Korea Academia-Industrial cooperation Society
    • /
    • v.15 no.5
    • /
    • pp.3345-3350
    • /
    • 2014
  • We investigated to assess the relationships between the major nutrients and phytoplankton dynamics during the spring season in 2010 and 2011 at 23 stations in Jinhae Bay, Korea. The bay is divided into four different zones based on pollutant sources and geographical characteristics. Nutrient limitation (>80%) was significant in Zone II, which is located in central bay and is influenced by the water well mixed from outer bay. The limited nutrient was followed in Zone III and IV that was occupying between 17% and 83%. However, the low levels are being kept below 35% in Zone I, which is characterized by the semi-enclosed eutrophic area of Masan and Haegam bays. Based on the PCA (principle component analysis) analysis, the nitrogen (N) sources in 2010 were particularly dominant and it may be due to the water mixing and wastewater formed from bottom layers and sewage. In 2011, major nutrients including nitrogen, silicon and phosphorus were dominant in the bay and are supplied by the river discharge after rainfalls with low salinity conditions. In particular, the N nutrients being supplied in 2010 are correlated with pennate diatoms Pseudo-nitzchia spp. and is not related to the phytoplankton population densities in 2011. The present study suggests that N sources play an important role in the proliferation of diatom, and the rapid nutrient uptakes by them are potential nutrient limitation factors in the bay.

Application of X-band polarimetric radar observation for flood forecasting in Japan

  • Kim, Sun-Min;Yorozu, Kazuaki;Tachikawa, Yasuto;Shiiba, Michiharu
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2011.05a
    • /
    • pp.15-15
    • /
    • 2011
  • The radar observation system in Japan is operated by two governmental groups: Japan Meteorological Agency (JMA) and the Ministry of Land, Infrastructure, Transport and Tourism (MLIT) of Japan. The JMA radar observation network is comprised of 20 C-band radars (with a wavelength of 5.6 cm), which cover most of the Japan Islands and observe rainfall intensity and distribution. And the MLIT's radar observation system is composed of 26 C-band radars throughout Japan. The observed radar echo from each radar unit is first modified, and then sent to the National Bureau of Synthesis Process within the MLIT. Through several steps for homogenizing observation accuracy, including distance and elevation correction, synthesized rainfall intensity maps for the entire nation of Japan are generated every 5 minutes. The MLIT has recently launched a new radar observation network system designed for flash flood observation and forecasting in small river basins within urban areas. It is called the X-band multi parameter radar network, and is distinguished by its dual polarimetric wave pulses of short length (3cm). Attenuation problems resulting from the short wave length of radar echo are strengthened by polarimetric wavelengths and very dense radar networks. Currently, the network is established within four areas. Each area is observed using 3-4 X-band radars with very fine resolution in spatial (250 m) and temporal (1 minute intervals). This study provides a series of utilization procedures for the new input data into a real-time forecasting system. First of all, the accuracy of the X-band radar observation was determined by comparing its results with the rainfall intensities as observed by ground gauge stations. It was also compared with conventional C-band radar observation. The rainfall information from the new radar network was then provided to a distributed hydrologic model to simulate river discharges. The simulated river discharges were evaluated again using the observed river discharge to estimate the applicability of the new observation network in the context of operations regarding flood forecasting. It was able to determine that the newly equipped X-band polarimetric radar network shows somewhat improved observation accuracy compared to conventional C-band radar observation. However, it has a tendency to underestimate the rainfall, and the accuracy is not always superior to that of the C-band radar. The accuracy evaluation of the X-band radar observation in this study was conducted using only limited rainfall events, and more cases should be examined for developing a broader understanding of the general behavior of the X-band radar and for improving observation accuracy.

  • PDF

Synthesis and Characterization of IPA-co-HDO-co-(TPA/MA) Anion-Exchange Membrane for All-Vanadium Redox Flow Battery (전바나듐계 레독스-흐름 전지용 IPA-co-HDO-co-(TPA/MA) 음이온교환막의 합성 및 특성)

  • Jung, Jae-Chul;Kwak, Noh-Seok;Hwang, Taek-Sung
    • Polymer(Korea)
    • /
    • v.35 no.6
    • /
    • pp.593-598
    • /
    • 2011
  • The IPA-co-HDO-co-(TPA/MA) copolymers for all-vanadium redox flow battery were synthesized by melt condensation polymerization using isophthalic acid(IPA), 1,6-hexandiol (HDO), terephthalic acid(TPA) and maleic anhydride(MA). The amination of chloromethylated IPA-co- HDO-co-(TPA/MA)(CIHTM) copolymer was carried out using trimethylamine, and the anion exchange membrane was also prepared by UV crosslinking reaction. The structure and thermal stability of IHTM copolymers were confirmed by FTIR, $^1H$ NMR, and TGA analysis. The anion membrane properties such as water uptake, ion exchange capacity, electric resistance and electrical conductivity, were measured by gravimetry, titration and LCR meter. The efficiency of the all-vanadium redox flow battery was analyzed. The ion exchange capacity, electric resistance and electrical conductivity were 1.10 meq/g, $1.98{\Omega}{\cdot}cm^2$, and 0.009 S/cm, respectively. The efficiency of charge-discharge, voltage, and energy for the allvanadium redox flow battery were 96.5, 74.6, 70.0%, respectively.

Adequacy Evaluation and Analysis of Discharge Measurement in Tributary Gaging Station by Analyzing Loop Rating Pattern Considering Backwater - For the Nam River and the Geumho River - (배수영향을 고려한 고리형 수위-유량관계 분석을 통한 지류 수위관측소의 유량관측 적절성 평가 및 분석 - 남강 및 금호강을 대상으로 -)

  • Kim, Kyung Dong;Kim, Dongsu
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2019.05a
    • /
    • pp.96-97
    • /
    • 2019
  • 최근 홍수파 또는 댐 및 보 등의 하천구조물과, 본류와 지류의 합류로 인해 발생하는 배수영향 등으로 하천에 고리형 수위-유량관계가 나타나 수위-유량관계식의 신뢰도가 저하된다고 판단되는 대하천 본류와 합류부 인근 지류에서 초음파유속계 측정결과를 지표로 활용하는 자동유량장치가 현재 58개소에 도입되어 운용되고 있다. 그러나, 최근 4대강 사업으로 다기능보 설치 등 하천의 수리적 특성이 변동하였고, 지류에서 주로 운용되는 수위-유량관계 기반 유량관측소에 본류-지류합류로 인한 배수영향으로 고리형 수위-유량관계가 형성될 경우, 지류 수위관측소를 자동유량관측소로 대체할 필요성이 제기되고 있다. 또한, 기존 자동유량관측소의 경우 홍수량이 기준이하로 나타날 경우 수위-유량 관계식으로 대체하는 경우도 있는데 적용기준이 명확하지 않다. 본 연구에서는 신규 자동유량장치 설치, 기존 수위관측소의 자동유량관측소대체, 자동유량관측소에서 수위-유량관계곡선의 활용 기준 마련 등을 명확하게 판단하기 위해 고리형 수위-유량관계 패턴을 분석하고 배수영향구간을 산정하여 제시하고 지류 수위관측소의 존속여부, 유량자료의 적절성을 판단할 수 있는 방법론을 제시 하고자한다. 이를 위해, 다기능보 설치 및 지류-본류로 인한 지류하천의 배수영향과 홍수파에 의한 고리형 수위-유량관계의 양상을 다양한 수문사상에 대해 분석하였다. 분석방법으로는 대상유역을 선정하고, 대상유역에 대한 부정류 모형을 구축하기 위해 HEC-RAS 모형을 사용하고, 대상유역에 대한 HEC-RAS 지형자료와 유량입력자료를 구축한다. HEC-RAS 부정류 모형 구축 후 다양한 본류와 지류의 유량을 가정하여 고리형 수위-유량 발생 패턴을 분석한다. 지류-본류 합류 영향의 경우, 본류영향의 유무를 고려하여 상승된 수위의 비율로부터 일정 유의수준을 두어 배수영향구간을 산정한다. 다기능보의 영향은 다기능보 설치 전후를 가정하여 영향을 분석한다. 본 연구에서는 시범구역으로 하류에 창녕함안보가 위치한 낙동강과 합류하는 남강과 하류에 달성보가 위치한 낙동강과 합류하는 금호강을 대상하천으로 하여 위 분석방법을 적용하였다. 분석결과 남강의 경우 지류 유량이 작은 경우에도 본류 유량이 상대적으로 큰 경우 배수영향에 의한 고리형이 나타났고, 다기능보에 의한 배수영향은 매우 적은 것으로 나타났다. 또한, 지류의 유량이 큰 경우 홍수파에 의한 고리형 관계도 발생함을 알 수 있었다. 금호강의 경우 또한 마찬가지로 본류의 유량에 의해 발생하는 고려형 관계가 지배적임을 알 수 있었다. 따라서 이러한 결과를 바탕으로 95%의 유의 수준을 기반으로 고리형 발생 여부를 확일할 수 있는 모노그래프를 남강, 금호강에 대해 작성하여 수위-유량관계 기반 유량 산정의 적절성을 판단하는 기준으로 제안하였다.

  • PDF

Analysis of Intercepted Discharge Considering the Change of Street Inlet Size (빗물받이 유입구 크기 변화에 따른 차집유량 분석)

  • Kim, Jung Soo;Ryu, Teak Hee;Rim, Chang Soo;Yoon, Sei Eui
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2015.05a
    • /
    • pp.375-375
    • /
    • 2015
  • 최근 도시화에 따른 도로구역의 증가는 도심지의 불투수면적의 확대 및 도로의 수로화를 야기하여 집중호우 발생 시 노면 및 저지대의 침수를 가중시키고 있으며, 이로 인한 인명피해 및 재산 피해는 점차 증가하는 추세이다. 이에 기후변화에 따른 강우량의 변화 및 국지성 집중호우 등을 고려하고 도심지 침수 방지 및 방재성능 강화를 위하여 과거에 간선 10년, 지선 5년에서 간선 30년, 지선 10년으로 관거의 설계빈도가 상향조정 되었다. 도심지의 도로 및 유역 유출량의 배수에 상당한 영향을 미치는 빗물받이 및 빗물받이 유입구는 도심지의 방재성능강화 목적에 맞는 설계가 이루어져야 하나 현재 국내의 빗물받이 설치기준은 하수도시설기준(2009, 환경부)에서 '빗물받이 크기별, 도로 차선별 적정 빗물받이 설치간격' 으로 설계빈도 5년을 반영하여 제시된 설치기준이므로 이에 대한 수정 및 개정이 필요한 실정이다. 이에 본 연구에서는 관거시설의 설계빈도 상향을 고려한 도로 빗물받이의 유입량 산정을 위해 빗물받이 유입구 크기별 유입량을 도로의 차선, 종경사, 측구 횡경사 및 설계빈도의 변화를 고려하여 도로에서의 유출량을 산정하였으며, Froude 상사법칙을 이용하여 수리실험모형을 제작하여 실험을 실시하였다. 빗물받이 제원은 현재 대부분의 국도에 설치되는 크기인 $40{\times}50cm$, $40{\times}100cm$$40{\times}150cm$를 축소 제작하였다. 측구의 유량은 도로의 차선(2~4차선), 경사(도로 종경사 2~10%, 측구 횡경사 2~10%) 및 설계빈도(최대 30년)을 고려하였다. 실험 결과 측구의 횡경사가 커질수록 빗물받이로 유입되는 유량은 증가하였으며, 빗물받이 유입부의 길이가 증가함에 따라 유입부 측면부를 통한 횡유입량을 증가시켜 빗물받이 유입부의 차집율을 증가시켰다. 또한, 도로의 조건 변화 및 유입량의 변화에 따라서 흐름폭이 감소할수록 차집율이 증가하였다. 수리실험결과를 토대로 도로의 조건변화(차선, 도로 종경사, 측구 횡경사) 및 설계빈도 변화를 고려한 빗물받이의 차집율을 제시하였고 이를 회귀분석 하여 빗물받이 유입구 크기별 유입량 산정식을 도출하였다. 이는 설계빈도 상향에 따른 국내 빗물받이 유입부의 설계 기준 제시에 기초자료로 활용이 가능할 것으로 판단된다.

  • PDF

Experiment of Discharge Supply Verification in River Experiment Center (하천실증연구센터 공급 유량 검증 실험)

  • Kim, Jong Min;Lee, Du Han;Kim, Sung Joong;Kim, Hyung Suk
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2018.05a
    • /
    • pp.77-77
    • /
    • 2018
  • 한국 건설기술연구원의 하천실증연구센터는 기존에 불가능했던 준 실규모 하천실험, 생태실험, 치수분야 대형 모형실험 및 하천관련 기준수립을 위한 치수분야 수요와 하천 구조물의 품질성능평가 등을 목적으로 설립된 대규모 하천실증연구센터이다. 안동 하천실증연구센터는 총면적 $193,051m^2$의 부지에 하천의 일반적인 형상을 갖고 있는 완경사, 급경사, 만곡수로가 있으며, 수리량 측정 실험을 위한 초음파 유속계, 프로펠러 유속계, 전자식 유속계 등의 유속 측정 장비를 보유하고 있고, Total Station, Lidar, RTK-GPS 등 지형 측량에 사용할 수 있는 장비들도 보유하고 있어, 하천관련 실험을 수행하기에 충분한 여건을 보유하고 있다. 하천실증연구센터는 2013년부터 내 외부활용에서 활용할 수 있도록 지원을 수행하였으며, 총 118건의 하천 수리 및 생태분야, 지반분야 등 다양한 토목분야의 연구자들이 하천실증연구센터를 방문하여 실험을 수행하였고, 국내 기관 및 대학뿐만 아니라 미국의 Iowa 대학, Idaho 대학, USGS, 네덜란드 Deltares, 프랑스 Irstea 등 외국의 대학 및 기관에서도 방문하여 실험 및 공동연구를 수행하였다. 안동 하천실증연구센터의 활용 범위와 활용 건수가 지속적으로 증가함에 따라 자체적인 검토 및 외부의 의견에 따라 센터의 실험시설에 대한 개선이 필요한 사항이 발생하게 되었다. 하천실증연구센터의 수로는 기존의 지반위에 성토를 하여 인공적으로 건설하였기 때문에 일반적인 실내의 실험수로와는 다르게 하상에서의 침투가 발생하며, 이에 따라 하류방향으로 흐름이 진행될수록 유량의 변동이 발생하게 된다. 이에 따라 침투로 인해 발생하는 유량의 변동이 하도 유량 공급의 안정화에 미치는 영향에 대한 신뢰성 검증의 필요성이 필요하다고 판단되어 펌프에서의 유량 공급 및 하도내에서의 유량 안정화에 대한 검증 실험을 수행하고자 한다. 유량 공급에 대한 검증 실험은 크게 두 단계로 구성하였다. 첫 번째 단계는 펌프의 유량 공급에 대한 검증실험을 수행하고, 두 번째 단계에서는 동일한 조건내에서 하류의 지하수 수심, 공급수조의 수심, 하도내의 수심과 ADV 및 ADCP로 측정된 유량과의 관계를 통해 하도내의 유량이 얼마나 일정하게 공급되고, 하류방향으로 갈수록 유량 손실의 경향에 대한 실험을 계획하였으며, 전체적인 실험은 장기적으로 각 수로에 대하여 수행하는 것을 목표로 하였으며, 본 실험단계에서는 완경사 수로를 대상으로 하였다.

  • PDF

Estimation of Flood Discharge using Satellite-derived Rainfall in Abroad Watershed - A Case Study of Pasig-Marakina, Phillippines - (위성강우를 이용한 해외 유역 홍수량 추정 - 필리핀 파시그-마라키나강 유역을 대상으로 -)

  • Kim, Joo Hun;Choi, Yun Seok;Kim, Kyeong Tak
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2018.05a
    • /
    • pp.398-398
    • /
    • 2018
  • OECD 발표에 의하면 물산업 관련 인프라 투자 전망은 전세계 GDP 대비 2010~2020년 약 1.01%에서 2020~2030년 약 1.03%로 확대될 전망으로 다른 통신, 전력, 철도 인프라 투자수요보다 많을 것으로 전망하고 있다(파이넨셜 뉴스, 2013.3.21.). 우리나라는 2005년 베트남 홍강종합개발사업을 시작으로 2015년 기준으로 세계 35개국에 진출하고 있다. 그러나 대부분의 물 산업 진출 대상 국가는 미계측 유역이 많고 지상에서 계측된 수문 자료가 부족한 실정이다. Namgung and Lee(2014)에 의하면 네팔의 수력발전소 건설에 관측된 강우량 자료가 없어 발전소 하류 10km 지점의 유하량 자료를 이용하여 자료의 정확도 검증을 대신하여 적용한 바 있다. 이와 같이 계측자료가 없거나 부족한 지역에 대하여 기상 위성을 이용하여 추정된 강수량 자료가 해당 지역의 강수 특성을 파악하는데 중요한 자료로 이용될 수 있다. 글로벌 위성 기반의 강수량 관측에 대한 역사는 1979년에 IR방법에 의해 위성으로부터 강우자료를 유도하는 개념이 도입된 이후 1987년 다중 채널의 마이크로파(MW) 복사계를 이용한 방법, 이후 두 IR과 MW를 혼합한 방법에서, 1997년 TRMM위성의 PR(Precpipitation Radar)의 레이더를 이용하는 방법, 그리고 2014년 GPM 핵심 위성(GPM Core Observatory)에 탑재된 Dual PR에 의한 방법으로 위성강수의 정확도를 매우 높여가고 있다. 본 연구는 KOICA 사업으로 진행중인 필리핀 메트로 마닐라 홍수조기경보 및 모니터링 체계 구축사업 중 파시그-마라키나강(Pasig-Marakina) 유역의 2012년 8월의 홍수사상에 대한 위성강우 및 글로벌 지형자료를 이용하여 홍수 유출량을 추정하는 것으로 목적으로 하고 있다. 유역내 6개 관측소의 일일 강우량 자료와 GPM IMERG 일강우량 자료 상관분석 결과 약 0.623, Bias는 -0.147, RMSE는 15.7정도로 분석되었다. 홍수량 분석은 2012년 8월 홍수가 발생한 시기인 2012년 8월 1일 00(UTC)부터 2012년 8월 16일 00(UTC)까지의 1시간 간격의 위성강우자료와 글로벌 지형자료를 이용하였고, 한국건설기술연구원의 MapWindow 기반 GRM 모형(mwGRM)을 이용하였다. 분석 결과 첨부홍수가 발생한 시기는 8월 7일 18:00(UTC)였고, 첨두 홍수량은 $4,073.9m^3/sec$로 분석되었다. 향후 수위-유량 관계식에 의해 정확도평가를 수행할 계획이다.

  • PDF

Estimation of flood peak discharge using flood marks (홍수흔적을 이용한 첨두홍수량 추정)

  • Lee, Tae Hee;Lee, Jung Hoon;Kang, Jong Wan;Roh, Youngsin
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2019.05a
    • /
    • pp.75-75
    • /
    • 2019
  • 첨두홍수량 자료는 홍수예경보 및 치수계획수립 등 하천관리에 있어서 매우 중요한 요소이다. 그러나 대규모 홍수가 발생 시 악천후가 동반된 기상상황이나, 현장 접근이 어려운 환경적 조건과 예산 및 인력 부족 등에 의한 불가피한 문제로 첨두홍수량을 측정하는데 어려움 있다. 따라서 일반적으로 수위-유량관계곡선식을 이용하여 첨두홍수량을 산정하지만 단순 고수위 외삽 추정을 통해 개발된 곡선식을 이용한 첨두홍수량 산정에 있어서는 주의가 필요하다. 이러한 경우 홍수가 지나간 후 현장조사를 통해 획득한 위치, 표고, 횡단면적 등 홍수흔적(flood marks)을 가지고 경사면적법(slope-area method)과 같은 간접적인 방법으로 첨두홍수량을 추정할 수 있다. 본 연구에서는 2018년 큰 호우사상이 발생한 내성천의 지류인 서천의 영주시(월호교) 지점과 남강의 산청군(하촌리) 지점에서 홍수흔적 조사를 통해 지점별 두 개의 단면을 선정하였다. 영주시(월호교) 지점의 두 단면 간 거리는 약 90m, 높이차는 약 0.21m로 조사되었고, 산청군(하촌리) 지점의 두 단면 간 거리는 약 330m, 높이차는 약 0.47m로 조사되었다. 경사면적법을 이용한 첨두 홍수량 추정에 적용된 조도계수는 '서천 하천기본계획(2014)', '남강 하천기본계획(2013)'에서 계획 홍수량 산정에 적용된 조도계수 0.029와 0.025를 적용하였다. 영주시(월호교) 지점은 2018년 9월 4일 발생한 호우사상의 첨두수위 5.59m에서 수위-유량관계곡선식을 이용하여 산정된 유량은 $1,127.8m^3/s$이고 경사면적법을 이용하여 추정된 유량은 $1,105.9m^3/s$로 약 -1.98%의 편차율이 발생하였다. 산청군(하촌리) 지점은 2018년 8월 26일 발생한 호우사상의 첨두수위 6.75m에서 수위-유량관계곡선식을 이용하여 산정된 유량은 $3,435.0m^3/s$이고 경사면적법을 이용하여 추정된 유량은 $3,233.3m^3/s$로 약 -6.24%의 편차율이 발생하였다. 경사면적법을 이용하여 추정된 첨두홍수량은 수위-유량관계곡선식을 이용하여 산정된 유량과 편차율이 지점별 ${\pm}10%$ 이내의 근사한 범위로 산정되었다. 따라서 경사면적법을 이용한 첨두홍수량 추정 방법의 적용에 있어서 적절한 것으로 판단된다.

  • PDF

Estimation of Nonpoint Source Pollutant Loads for Rural Watershed by AvSWAT (AvSWAT를 이용한 농촌유역 비점원 오염물질 부하량 예측)

  • Kim, Jin-Ho;Lee, Jong-Sik;Kim, Won-Il;Jung, Goo-Bok;Han, Kuk-Heon;Ruy, Jong-Su;Kim, Suk-Cheol;Yun, Sun-Gang;Lee, Jeong-Taek;Kwun, Soon-Kuk
    • Korean Journal of Soil Science and Fertilizer
    • /
    • v.40 no.1
    • /
    • pp.12-17
    • /
    • 2007
  • This study was conducted to evaluate the characteristics of nonpoint source pollutants discharge from a small rural watershed. A typical rural area in Gongju City, Korea, was selected as the research site. Water quality and quantity in streams and rainfall samples were analyzed periodically from May to October 2005. Pollutant loads were estimated from a nonpoint source pollution model (AvSWAT, Arcview Soil and Water Assessment Tool). During the rainy season, from June 26 to 30 September 2005 and the dry season, before 26 June and after 30 September 2005, biological oxygen demands and chemical oxygen demands accounted for 91.3% and 93.7% of annual load, respectively, while total-N and total-P were 97.1% and 91.1% of annual load, respectively. The observed stream flow was $66.5m^3sec^{-1}$, while simulation stream flow was $66.2m^3sec^{-1}$. That can be assumed that simulation can be used to estimate the stream flow without practical measurement. However, the runoff trend following the occurrence of a storm event was not recorded properly.