DIGITALNA ARHIVA ŠUMARSKOG LISTA
prilagođeno pretraživanje po punom tekstu
ŠUMARSKI LIST 13/2005 str. 221 <-- 221 --> PDF |
IZLAGANJE NA ZNANSTVENOM SKUPU PRESENTATION AT THE INTERNATIONAL SYMPOSIUM Šumarski list SUPLEMENT (2005), 219-228 UDK 630* 181.3 PROCJENA RIZIKA KOD EKSTREMNIH HIDROLOSKIH STANJA SLUČAJ SLIVOVA VSEMINKE I DREVNICE, REPUBLIKA ČEŠKA RISK ASSESSMENT OF EXTREME HYDROLOGICAL SITUATIONS CASE STUDY OF THE VSEMINKAAND DREVNICE CATCHMENTS, CZECH REPUBLIC Pavel KOVÄR* SAŽETAK: Primjena modela vodne ravnoteže daje mogućnost identifikacije ekstrema u procesima padavina/otjecanja. Naime, sadržaj vlage u tlu u gornjoj zoni malih vodnih slivova može hiti pokazatelj potencijalnih poplava ili suša. Ove se zone gornjih slivova (0.20 m—0.60m) obično nazivaju "aktivnim zonama " i jako utječu na oblikovanje izravnog otjecanja, posebice kada je njihov sadržaj vlage blizu zasićenja ili vodi prema suši, kada je iscrpljen sadržaj vode. I šumsko gospodarenje, te dobra šumarska praksa, igraju značajnu ulogu u procesima padavina/otjecanja tj. u procesima održavanja vode u gornjim slojevima slivova. Ova je uloga još važnija u malim slivovima s jako fluktuirajućim izlijevanjima vode. Članak prikazuje kako se može izvesti simulacija ovih hidroloških procesa primjenom modeliranja vodne ravnoteže. U tu svrhu primijenjen je model WBCM-5 (Water Balance Conceptual Model version 5), u jednodnevnim fazama na nekoliko pokusnih slivova u Republici Češkoj. Prikazana su dva predpoplavna razdoblja tik pred katastrofalnim poplavama u srpnju 1997. Isto tako, iz povijesnih su zapisa identificirana razdoblja suše. Sušni se trendovi mogu naći primjenom modela, pokazujući kako su osjetljivi slivovi u bezpadavinskim razdobljima. Konceptualni sadržaj vlage u tlu odnosi se na API indeks (Antecedent Precipitation Index). Kvantitativna studija također je izrađena o utjecaju uporabe zemljišta i gospodarenja na izravnu formaciju otjecanja i na nadoknadu podzemne rezerve. U radu se prikazuju rezultati modela za poplavno razdoblje iz srpnja 1997. u dva pokusna moravska sliva, Vseminka (48 %šume) i Drevnice (81 % šume) te ih uspoređuje s onima iz sliva Rokvtenka (24 % šume) u istočnoj Češkoj. Vidljivo je da dubina aktivne zone u slivu ima dominantan utjecaj na stupanj rizika od poplave i/ili suše. Ovi se hidrološki ekstremi mogu reducirati djelomično i samo u maloj mjeri. Međutim, racionalno gospodarenje slivom i dobro orgnaizirana uporaba zemljišta temeljena više na zaštitnu funkciju nego na proizvodnju, preduvjet je za zaštitu od opasnih hidroloških pojava s katastrofalnim posljedicama. KIjučne riječi: model vodne ravnoteže, prethodni uvjeti vlage, hidrološki ekstremi, deficit vlage u tlu, scenarij simulacije * Pavel Koväf, Češko poljoprivredno sveučilište Prag, Fakultet šumarstva i okoliša, Kamycka 129, 165 21 Prague 6-Suchdol, Czech Republic, E-mail: kovar@fle.czu.cz |
ŠUMARSKI LIST 13/2005 str. 222 <-- 222 --> PDF |
P. Kovan PROCJENA RIZIKA KOD EKSTREMNIH MIDROLOŠKIH STANJA SLUČAJ SLIVOVA ... Šumarski list SUPLEMENT (2005), 219-228 UVOD - Introduction Vodna se ravnoteža nekog sliva može dobro kvantificirati analizom komponenata vodne ravnoteže s obzirom na neposredno otjecanje i dopounjavanje podzemne vode. Posebice je gornji sloj sliva od ogromne važnosti, jer njegov sadržaj vlage predodređuje moguće poplave ili suše - prema deficitnoj vrijednosti njegove vlage u tlu(Maidment 1993). Stoga je taj sloj (obično 0.20 m-0.50 m u dubinu) poznat kao aktivna zona. Utjecaj uporabe zemljišta na vodnu ravnotežu u aktiv noj zoni je također stvar od velike važnosti u gospodarenju vodnim resursima i kontroli poplave (Falken mark, 1999). Studija temeljena na izračunu vodne ravnoteže izrađena je pomoću podataka s dva mala slivna područja, Vseminka i Drevnice u Moravskoj. U širem smislu, cilj je ove studije istražiti utjecaje tla, fiziografske čimbenike i uporabe zemljišta na kapacitete vodne retencije u aktivnim zonama slivnog područja. EKSPERIMENTALNI SLIVOVI - Experimental catchments Oba su pokusna slivna područja sa sličnom površi2001.) Slika 1 daje opći prikaz situacije dvaju slivova. nom nedavno pretrpjela barem dvije poplave (1997. i Slika 1. Situacija eksperimentalnih slivova Figure I Situation of experimental catchments Tablica la. Karakteristike tla u slivovima Table I a. Catchment soil characteristics Prosječni kapacitet Ime sliva Površina (km") terena (KT) Catchment Area (km´) Average field name capacity FC Vseminka 21,51 0,40-0,42 Drevnice 22,58 0,40-0,43 Tablica lb Hidrološke grupe tala u slivovima Table 1 b. Catchment hydrological soil groups Grupa tla - Soil group Vseminka A 0,0 % B 10,8 % C 46,1 % D 43,1 % Drevnice 5,1% 30,5 % 39,0 % 25,4 % Oblik je Vseminkc Prosječna poroznost POR Average porosity POR 0,45-0,50 0,44-0,48 Dubina aktivne zone DROT (mm) Depth of active zone DROT 300,0 300,0 longitudinalan, gornji dijelovi N Sadržaj saturacije QSAT (mm) Saturation content QSÄT 123,0 125,0 središnje doline pošumljeni su sve do vodnog sliva, a uporaba zemljišta na ostalim dijelovima je različita. Sliv Drevnice je u obliku lepeze, vrlo je pošumljen, a samo njegov donji dio je pod boljom strukturom uporabe zemljišta. Karakteristike tla slivova prikazane su u Tablici la, b. Grupe tala podijeljene su prema najšire primijenjenim metodama (U.S. SCS, 1986). Uporabe zemljišta i fiziografski čimbenici predstavljeni su u Tablici 2. |
ŠUMARSKI LIST 13/2005 str. 223 <-- 223 --> PDF |
P. Kovar: PROCJENA RIZIKA KOD EKSTREMNIH HIDROLOSKIH STANJA SLUČAJ SLIVOVA ... Šumarski list - SUPLEMENT (2005). 2 19-228 MODEL - Model Model WBCM primijenjen s ciljem kvantificiranja dinamike vlage u tlu pokusnih slivova je model s mogućnošću paramctarske distriubcije na toj površini. On se temelji na pristupu integrirane zalihe. Svaki element zalihe predstavlja intercepciju prirodne zalihe, površinu tla i aktivnu zonu. U toj verziji nisu primijenjene sve nezasićene zone i zone podzemne vode. Međutim, uzimala se u obzir samo nadopuna vode u njih kao output iz aktivne zone u obliku duboke infiltracije. Model s dnevnim stopama inputa i outputa simulira sljedeće procese (Kovar 1999. i Kovar 2000.): potencijalna evapotranspiracija, intercepcija i neposredna padavina na tlo izravno otjecanje dinamika aktivne zone vlage tla Opaska: Sadržaj vlage tla čitave nezasićene zone, dinamike podzemnih voda, bazno otjecanje i ukupno otjecanje obračunati su samo za proces parametarske kalibracije. Modificirana metoda Monteith-Penman primijenjena je za izračun potencijalne evapotranspiracije. Linearna distribucija lokalnih kapaciteta intercepcije je u načelu uzeta za aktualnu intercepciju i izravnu padavinu. Za kvantificiranje izravnog otjecanja primijenjena je metoda US Soil Conservation Service (SCS) temeljena na određivanju broja krivulje otjecanja (CN). Nadopunjavanje aktivne zone i njeno pražnjenje najviše ovise o parametrima tla (kapacitet terena, poroznost, hidrauličke provodljivosti, KS), sadržaju vlage u tlu i potencijalnim uvjetima evapotranspiracije. Za ovaj izračun primijenjena je jednodimenzionalna Richardova jednadžba u obliku konačne diferencije (Kovar 2001). Tablica 2. Uporaba slivnog zemljišta i fiziografski faktori Table 2 Catchment land use and physiographic factors Uporaba zemljišta -Land use Vseminka Drevnice Obradivo zemljište -Arable land Livade, pašnjaci -Meadows, pastures Divlje zelenilo- Wild greenery Suma -Forest Ostalo -Others Fiziografski faktori -Physiographic factors Dužina rijeke - Length of river Prosječni nagib rijeke - Average river slope Prosječni nagib sliva -A verage catchment slope 9,3 % 2,7 % 24,2 % 9,0 % 14,0% 0,8 % 48,2 % 81,0% 4,3 % 6,5 % Vseminka Drevnice 9,2 km 9,1 km 3,6 % 3,5 % 19,4% 23,0 % Prosječna altituda sliva - Average catchment altitude 400 m n.v. 495 m n.v. Tablica 3. Glavni podaci o poplavama u srpnju 1997. Table 3 Major data onfloodsin 07/1997 Sliv Datum Catchment Date 4/7-10/7 Vseminka 18/7-10/7 4/7-10/7 Drevnice 18/7-27/7 Vrhunac poplave Qmax (mV) Flood peak Qmax (m3.s´) 11,4 7,3 11,2 5,3 Kauzalna oborina (mm) Causal rainfall (mm) 123,3 70,4 160,6 90,1 REZULTATI I DISKUSIJA Ovaj je model primijenjivan od 1992. za izračun komponenata vodne ravnoteže. Na taj su način parametri modela već bili kalibrirani. Cilj ovog članka je proučiti dva predpoplavna razdoblja u srpnju 1997. s obzirom na sadržaj vlage u tlu u aktivnim zonama obaju pokusnih slivova. Dva poplavna vala preko teritorija istočne Češke i Moravske dogodila su se od 4. do 10. Početni SMD Initial SMD 37,8 10,8 36,1 12,4 Indeks prethodnih oborina API30 (mm) A n tec eden t precip ita Hon index APL3n (mm) 30,0 87,8 35,9 127,2 - Results and discussion Završni SMD (mm) Final SMD (mm) 1,2 0,0 0,8 0,0 srpnja, te od 18. do 27. srpnja 1997. Glavni podaci o ovim poplavama prikazani su u Tablici 3. Zatim je primijenjen model WBCM-5 za izračun fluktuacije sadržaja vlage u tlu u obliku vrijednosti dnevnog deficita vlage (SMD), kao i dnevni indeks prethodnog poplavljivanja (API). Razdoblje od 1. lipnja do 3. rujna 1997. pokazuje kako se opasnom čini poplavna situacija kada |
ŠUMARSKI LIST 13/2005 str. 224 <-- 224 --> PDF |
P. Kovan PROCJENA RIZIKA KOD EKSTREMNIH H1DROLOŠKIH STANJA SLUČAJ SLIVOVA Šumarski list - SUPLLMENT (2005). 219-228 vrijednosti SM D padnu ispod 10 mm i još niže, a slivno područje više nema kapacitet retencije za zalihu sljedeće kiše. Posebice je bilo opasno razdoblje od 5. srpnja do 7. kolovoza 1997. na Vseminki i Drevnicama. Došlo je do najveće poplave u posljednjih 100 godina, poplave koja je odnijela mnoge ljudske živote i uzrokovala veliku materijalnu štetu. Drugi dio ovog istraživanja bio je odabir dvaju ekstremna razdoblja: suhog i mokrog, otprilike istog trajanja. Iz zapisa podataka je za suhu godinu s najmanje oborina i najvišom temperaturom zraka u vegetacijskom razdoblju (travanj-listopad) u posljednjih 20 godina izabrana godina 1992. Slično tome je za vlažnu godinu izabrana godina 1997. s najvišim oborinama određenim iz zapisa podataka istih stanica lociranih odgovarajuće za oba eksperimentalna sliva. Međutim, sljedeći cilj ovog članka je izabrati dva razdoblja, jedno tipično suho, te drugo tipično vlažno, oba u trajanju od oko 6 do 7 tjedana, kada nije bilo nikakve kiše u suhoj sezoni, niti bitnih olujnih kiša kombiniranih s regionalnim kišama s njihovom visokom periodičnosti (p<0,01, N>100 godina) u vlažnoj sezoni. Tako su ovi zapisi o dnevnim kišama i odgovarajućim dnevnim otjecanjima analizirani za obje karakteristične godine, 1992. i 1997., te su za slivove Vseminka i Drevnice odgovarajuće odre đena sljedeća razdoblja: Suho razdoblje 1992.: 16. srpanj do 31. kolovoz 1992. (47 dana) Vlažno razdoblje 1997: 26. lipanj do 14. kolovoz 1997. (50 dana) Ova su razdoblja dalje analizirana, te su rekonstruirane komponente njihove dnevne vodne ravnoteže. Za rekonstrukciju izabranih suhih i vlažnih razdoblja primijenjeni su ranije kalibrirani modeli parametra WBCM za čitava vegetacijska razdoblja za godine 1992. odnosno 1997 (Kovar et al. 2001). Osim prethodne analize situacije u godini 1997., situacija 1992. je bila potpuno drukčija, a aktivna se zona zatvorila do točke sušenja. Dnevne vrijendosti SMD (mm) su recipročne dnevnim vrijednostima prethodnih oborina API (mm). Kada uspoređujemo oba eksperimentalna sliva, evidentno je daje sliv Vseminka bolji "držač vode" od sliva Drevnice, unatoč većoj šumskoj površini Drevnica. Konačni cilj ovoga rada je simulirati hipotetične situacije gdje oba sliva odgovaraju na signifikantnu kišu u suhim i vlažnim uvjetima. Ove su kiše bile projektnog karaktera s vremenom ponavljanja N (godine): 2, 5, 10, 20, 50, te 100 godina. Te su kiše izvedene za stanicu Velikovä iz publikacije Samaj , et al. 1983.: Povratno razdoblje N - godina -Return period N -year N= 2 5 St. Velikovä (mm) 38,´ 51,7 Zbog kratkoće izvješća, ovdje dajemo samo projektirane 100-godišnje rezultate kiše. Ovaj je kišni scenarij trebao pasti krajem suhog razdoblja, 20. kolovoza 1992. te na kraju vlažnog razdoblja 27. srpnja 1997, jednakomjerno raspoređen na oba slivna područja Vseminke i Drevnice. 10 20 59,9 68,4 50 100 78,9 87,1 U suhom razdoblju 1992. u slivu Vseminka, utjecaj 100-godišnje kiše povećao je izravno otjecanja s 0,5 mm na 7,3 mm, a u slivu Drevnice u istom razdoblju te vrijednosti variraju od 4,0 mm do 32,0 mm. Tablica 4 Komparativna vodna ravnoteža u slivu Vseminka u suhom razdoblju od 16. srpnja do 31. kolovoza 1992. te vlažnom razdoblju od 26. lipnja do 14. kolovoza 1997. Table 4 Comparative water balance in the Vseminka catchment dry period 16/7 to 31/8/1992, wet period 26/6 to 14/8/1997 Komponenta vodne ravnoteže (mm) Water balance component (mm) Kiša, RAIN Rainfall, RAIN Ukupno otjecanje, STF Total runoff, STF Direktno otjecanje, SOF Direct runoff, SOF Stvarna evapotranspir., AE Actual evapotranspir, AE Promjene u zalihi AW Change in storage AW Suha situacija Vlažna situacija Scenarij sa 100-ljetnom kišom 1992. (mm) 1997. (mm) 1992 (mm) Dry situation Wet situation Scenario1992 (mm) 1997 (mm) 1992 (mm) 18,2 389,9 105,3 1,0 151,9 7,8 0,5 103,1 7,3 90,6 120,0 101,9 -71,8 118,0 -4,4 1997 (mm) with 100-year rain 1997 (mm) 477,0 208,2 148,2 121,4 147,4 |
ŠUMARSKI LIST 13/2005 str. 225 <-- 225 --> PDF |
P. Kovär: PROCJENA RIZIKA KOD EKSTREMNIH HIDROEOŠKIII STANJA SLUČAJ SLIVOVA . Šumarski list - SUPLEMENT (2005), 219-228 Tablica 5 Komparativna vodna ravnoteža u slivu Drevnice u suhom razdoblju od 16. srpnja do 31. kolovoza 1992. te vlažnom razdoblju od 26. lipnja do 14. kolovoza 1997. Table 5 Comparative water balance on the Drevnice catchment dry period 16/7 to 31/8/1992, wet period 26/6 to 14/8/1997 Suha situacija Komponenta vodne ravnoteže (mm) 1992. (mm) Water balance component (mm) Dry situation 1992 (mm) Kiša, RAIN Rainfall, RAIN 18,2 Ukupno otjecanje, STF Total runoff, STF 4,0 Direktno otjecanje, SOF Direct runoff, SOF 0,6 Stvarna cvapotranspir, AE Actual evapotranspir., AE 90,8 Promjene u zalihi AW Change in storage AW -76,6 Ni u jednom slivu ovo povećanje izravnog otjecanja nema katastrofalni karakter. U slivu Vseminka to je bilo +6,8 mm, a u Drevnicama samo +28,0mm. Međutim, za vrijeme vlažnog razdoblja ovo je povećanje signifikantnije. U Vseminki to je bilo +45,1 mm, a u Drevnici čak +81,3 mm. To znači da u drugom slučaju u slivu plitkog tla s potpuno natopljenom aktivnom zonom gotovo 93 % kiše može dovesti sliv u stanje izravnog otjecanja. To bi nesumnjivo dovelo do katastrofe. Povećanje potpore podzemne zalihe također je značajno u rezultatima svih četiri scenarija. Tu su prikazane 100 Vlažna situacija Scenarij sa 100-ljetnom kišom 1997. (mm) 1992 (mm) 1997 (mm) Wet situation Scenario with 100-year rain 1997 (mm) 1992 (mm) 1997 (mm) 389,9 105,3 477,0 232,0 32,0 308,2 207,2 26,1 288,5 117,5 102,2 119,8 39,4 -28,9 147,4 postojeća situacija i ona kada padne projektirana kiša od 87,1 mm. Niske vrijednosti deficita vlage tla, SMD (te odgovarajuće visoke vrijednosti indeksa prethodne oborine, API) predstavljaju situaciju gdje više nema retencije da nadopune podzemne zalihe u slivu. Slična je situacija u slivu Drevnice, čak još više zbog brze fluktuacije uslijed pliće aktivne zone. Tablice 4 i 5, te Slike 2 i 3 prikazuju rezultate vrijednosti sastojaka komparativne vodne ravnoteže. Slike što predstavljaju sliv Drevnice slične su, ali nisu priložene ovom članku zbog limitiranog prostora. 90 80 70 E E 60 c ´5 I 50 40 30 20 10 date i ram I rain N=100y Qobs. . Q comp.sim. N=100y Slika 2. Hidrografi u slivu Vseminka u suhom razdoblju 1992. Figure 2 Hydrograph* in the Vseminka catchment in the dry period 1992 |
ŠUMARSKI LIST 13/2005 str. 226 <-- 226 --> PDF |
P. Kovar: PROCJENA RIZIKA KOD EKSTREMNIH HIDROLOSKIH STANJA SLUČAJ SLIVOVA Šumarski list - SUPLEMENT (2005), 219-228 I date I rain N=100y Qobs. Q comp.sim. N=100y Slika 3. Hidrografi u slivu Vseminka u vlažnom razdoblju 1997. Figure 3 Hydrographs in the Vseminka catchment in the wet period 1997 ZAKLJUČCI - Podaci o hidrometeorološkim opservacijama, tlu, uporabi zemljišta i gospodarenju slivnim područjem provjereni su s obzirom na pouzdanost, te su zatim prostorno obrađeni s ARC/INFO u slivovima Vseminka i Drevnice. Zatim, nakon primjene modela WBCM-5 na te podatke, dobiveni se outputi analiziraju te se mogu donijeti sljedeći rezultati: Komparacijom padavina/otjecanja srpnja 1997 u pokusnim slivovima, unatoč relativno većem postotku pošumljenosti u slivu Drevnice (81,0 % prema 48,2 % u Vseminki), količina oborina (Drevnice oko 30 % više nego u Vseminki) i njen intenzitet su važniji čimbenici od uporabe zemljišta s obzirom na poplavu. Modeliranje vodne ravnoteže kao što je primjena samo "gornjih slojeva" (tj. aktivnih zona) u slivovima, može jednostavno kvantificirati vrijednosti dnevnog deficita vlage u tlu (SMD). U suhom razdoblju godine 1992. izravno otjecanje iz scenarijem projektiranih oborina od 87,1 mm bilo je Conclusions u pokusnom slivu Vseminka 13 %, a u Drevnicama 32 %. U vlažnom razdoblju godine 1997. izravno otjecanje iz tih oborina bilo je 52 % u Vseminki, a u Drevnicama čak 93 %. Stupanj situacije izražene bilo s API ili SMD vrlo je znatan u procesu izravnog oblikovanja otjecanja. Promjene u podzemoj vodi su logične, velika poboljšanja u nadoknadi vode u suhim razdobljima ili obim slivovima su znatna, prekomjernost u vlažnim radobljima nije opasna, jer je tok podzemnih voda iz njihovih podzemnih zaliha procijenjen kao ne suviše brz. Male promjene aktualne evapotranspiracije ovise o nekoj aktivnoj zoni i čitavom nenatopljenom vodnom sadržaju tla. Ove promjene nisu bitne. Modelom WBCM može se koristiti samo za kvantifikaciju individualnih komponenata vodne ravnoteže, ali isto tako za identificaciju potencijalne štete od poplave i/ili suše. ZAHVALA - Acknowledgement Izražavamo zahvalnost za podršku Ureda za znanvećanja ekološke stabilnosti, recencije i akumulacije stvenu podršku Ministarstva obrazovanja, Republika vodnih kapaciteta u krajoliku ". Češka, Projekt br. MSM 414100008 "Mogućnosti po |
ŠUMARSKI LIST 13/2005 str. 227 <-- 227 --> PDF |
P. Kovär: PROCJENA RIZIKA KOD EKSTREMNIH HIDROLOSKIH STANJA SLUČAJ SLIVOVA ... Šumarski list - SUPI.F.MENT (2005), 219-228 LITERATURA - References Cudlin, P., F. Zemek, M. Hečman, P. Kovär, (2000): Evaluation of the extreme hydrological phenomena conditions in landscape using methods of mathematical modelling and GIS (in Czech), Workshop Proceedings 2000, pp. 261-272, Publishing House: CVUT Prague, ISBN 80-01-02318-4, Grant GACR 103/99/1470 Extrcmni hydrologicke jevy v povodich,. Falken m ark, M., (1999): Water, a reflection of land use, 128 pp., Published by Swedish Natural Science Research Council and UNESCO, ISBN 91546 03587,. US SCS, (1986): Urban hydrology for small watersheds. Technical Release 55: 13, 48 pgs. GACR Report, (1996): Implementation of Hydrological Models on the Rokytenka Catchment. In: Transport Processes in Head-Water Areas, GACR 103-94-0653, 66pgs. Ko vär, P., (1999): Modelling of Significant Rainfall- Runoff Events. In: Hydrological Processes Mo delling, pp. 181-207, Kluwer Academic Publisher. Amsterdam/Dordrecht, NL. ISBN-0-62855881- 8. Kovär, P., P. Cudlin, M. Korytär, F. Zemek, M. Hečman, (2001): Comparative study of water balance on the experimental catchments Všeminka and Drevnice, pp. 260-266, Rostlinnä vyroba, 47, 2001 (6). Maidment, D.R., (1993): Handbook of Hydrology. Published by McGraw-Hill, Inc. n ISBN 0-07039732- 5. Šamaj,F., R. Bräzdi 1, J. Val ović, (1983): Dennć ührny zräzok s mimoriadnou vydatnoscou v ČSSR v obdobi 1901-80. (in Slovak). (Daily rainfalls of maximum intensity in ČSSR in the period 1901-80, Proceedings of Slovak Hydrometeorological Institute). Sbornik präc. SHMU, Nakladatelstvi Alfa, Bratislava. |
ŠUMARSKI LIST 13/2005 str. 228 <-- 228 --> PDF |
PRESENTATION AT THE INTERNATIONAL SYMPOSIUM Sumarsk i list - SUPLEMENT (2005). 219-228 RISK ASSESSMENT OF EXTREME HYDROLOGICAL SITUATIONS CASE STUDY OF THE VSEMINKA AND DREVNICE CATCHMENTS, CZECH REPUBLIC Pavel KOVÄR* SUMMARY: Use of water balance models gives an opportunity to identify extremes in the rainfall-runoff processes, namely soil moisture content in the upper soil zone of small head-water catchments can indicate potential floods or drought. These upper catchment zones (0.20-0.60 m) are usually called "active zones" and strongly influence the formation of direct runoff, particularly when their soil moisture content is close to saturation, or lead to drought when the water content is exhausted. Also forest management and good forestry practices play a significant role in rainfall-runoff processes, namely in the water-holding processes of upper catchment layers. This role is even more important on small catchments with highly fluctuating water discharges. The paper shows how a simulation of these hydrological processes can be performed using water balance modelling. The WBCM-5 model (Water Balance Conceptual Model, version 5) running in one-day steps on several experimental catchments in the Czech Republic was implemented for this purpose. The two preflooding periods just before the catastrophic floods in July 1997 are clearly shown. Similarly, periods of drought can also be well identified from the historical records. The trends of "drying" can be found by implementing the model, and they show how sensitive a catchment is in rainless periods. Conceptual soil moisture content is related to the Antecedent Precipitation Index (API). A quantitative study is also made of the impact of land use and management on direct runoff formation and on subsurface storage replenishment. The paper gives the model results for the flood period of July 1997 on two Moravian experimental catchments, Vseminka (48 %forested) and Drevnice (81 % forested), and compares them with those from the Rokytenka catchment (24 % forested) in Eastern Bohemia. It is evident that the depth of the active zone on a catchment has a predominant influence on the degree of flood and/or drought risk. These hydrological extremes can be reduced partially and only to a small extent. However, rational catchment management and well organised land use and forest management based more on the protection function than on production is prerequisite for dangerous hydrological events with catastrophical consequences. Key words: water balance model, antecedent moisture conditions, hydrological extremes, soil moisture deficit, scenario simulation. INTRODUCTION Water regimes of a catchment can be well quantified particular, the upper soil layer in a catchment is extreby an analysis of water balance components with refemely important, as its moisture content predetermines rence to direct runoff and subsurface water recharge. In possible floods or droughts according to its soil moistu re deficit value (Maidment , 1993). Therefore, this * Pavel Koväf, Czech University of Agriculture Prague, layer (usually 0.20 m to 0.50 m in depth) is known as Faculty of Forestry and Environment, Kamycka 129, the active zone. The impact of land use on the water 165 21 Prague 6-Suchdol, Czech Republic E-mail: kovar@nc.czu.cz balance in an active zone is also a matter of great im |
ŠUMARSKI LIST 13/2005 str. 229 <-- 229 --> PDF |
P. Kovär: RISK ASSESSMENT OF EXTREME HYDROLOGICAL SITUATIONS CASE STUDY ... Šumarski list SUPLEMENT (2005). 2I9-22S portancc in water resources management and flood control (Falkenmark, 1999). A case study based on water balance computation has been carried out using data from two small catchments Vseminka and Drevnice in Moravia. The aim of EXPERIMENTALBoth Moravian experimental catchments covering similar areas have suffered at least twice from floods during recent years (1997, 2001). Fig. 1 gives a general view of the situation of the catchments. The shape of the Vseminka catchment is longitudinal, the upper sides of the central valley are forested down from the water divide, and land use on the other parts is diversified. The Drevnice catchment is fan-shaped, highly forested and only its lower part has a more structured land use. The soil characteristics of the catchments are given in Table 1 a, b. The soil groups have been selected in accordance with methods widely used (U.S. SCS, 1986). The land use and physiographic factors are presented in Table 2. The WBCM model that was implemented with the aim of quantifying the soil moisture dynamics of the experimental catchments is a lumped model with probability parameter distribution over the area. It is based on the integrated storage approach. Each storage element represents the natural storage interception, soil surface and active zone. In this version not all the unsaturated and ground water zones were used, however only a recharge to them was considered as an output from the active zone in the form of deep infiltration. The mo- RESULTS ANDThe model has been implemented annually since 1992 to compute the water balance components. Thus the model parameters had already been calibrated. The aim of this paper is to study two pre-flooding periods in July 1997 with reference to soil moisture content in the active zones of both experimental catchments. There were two flood waves over the territory of Eastern Bohemia and Moravia (4/7 to 10/7 and 18/7 to 27/7/ 1997). The major data on these floods is shown in Table 3. Then, the WBCM-5 model was implemented to compute soil moisture content fluctuation in the form of daily moisture deficit values (SMD), as well as the daily antecedent precipitation index (API). The period of June 1 to September 3, 1997, shows how dangerous flooding situation appears when the SMD values drop below 10 mm and even lower and the catchments have no more retention capacity to store next rainfall. In particular, the period 5/7 to 7/8/ 1997 on the Vseminka and Drevnice was the most hazardous. The result was a flood close to 100-year return period, which caused major loss in human life and in property. The second this study, in a broader sense, is to investigate the effect of soil, physiographical factors and land use on water retention capacities in the active zones of the catchments. CATCHMENTS del with daily step input/output rates simulates the following processes (Kovar 1999, Kovar 2000): potential evapotranspiration, interception and throughfall direct runoff active soil moisture zone dynamics Note: Soil moisture content of the whole unsaturated zone, ground water dynamics, base flow and total flow were computed only for the process of parameter calibration. The modified Monteith-Penman method was used for computing potential evapotranspiration. A linear distribution of local interception capacities was in principle used for actual interception and throughfall. For quantifying direct runoff, the US Soil Conservation Service, SCS method based on runoff curve number (CN) assessment was used. The recharge of the active zone and its depletion depend mostly on soil parameters (field capacity FC, porosity POR, hydraulic conductivities, KS), on soil moisture content and on potential evapotranspiration conditions. For this computation, the one-dimensional Richards equation is used in the finite difference form (Kovar, 2001). DISCUSSION part of this research was to select two extreme periods: dry and wet ones of about the same duration. As a dry year with the least precipitation and the highest air temperature in vegetation period (April - October) in the last 20 years the year 1992 was selected in the data record. Similarly, as a wet year was selected the year 1997 with the highest precipitation determined from the data record from the same stations located appropriately for both experimental catchments. However, the next aim of this paper was to select two periods, one typically dry and the other typically wet, both of the duration about 6 to 7 weeks when there was either no rain in dry season or substantial torrential storms combined with regional rainfalls with their high periodicity (p< 0.01, N > 100 years) in the wet season. Thus these daily rain and corresponding daily runoff records were analysed in both characteristic years 1992 and 1997 and for the Vseminka and Drevnice catchments the following periods were adequately selected: Dry period 1992: July 16 to August 31,1992 (47 days) Wet period 1997: June 26 to August 14,1997 (50 days) |
ŠUMARSKI LIST 13/2005 str. 230 <-- 230 --> PDF |
P- Kovär: RISK ASSESSMENT OF EXTREME HYDROLOG1CAL SITUATIONS CASE STUDY . Šumarski list SUPLEMENT (2005), 219-228 These periods were further analysed and their daily water balance components reconstructed. For the reconstruction of the selected dry and wet periods were used the WBCM model parameters calibrated earlier for the whole vegetation periods of the years 1992 and 1997 resp. (Kovar etal, 2001). Next to the previous analysis of the situation in 1997 the situation in 1992 was completely different and the active zone was closed to the wilting point. Daily SMD (mm) values are reciprocally to daily antecedent precipitation index values, API (mm). When we compare both experimental catchments, it is evident that the Vseminka catchment is more "water keeping" than the Drevnice catchment in spite of larger forested area on the latter case. The final aim of this work was to simulate the hypothetical situations how both experimental catchments respond to significant rainfall in dry and wet conditions. These rainfalls were those of the design character with the recuiTcnce time N (years): 2, 5, 10, 20, 50, and 100 years.These rainfalls were derived for the Velikovä station from the publication (Samaj, et al,1983): For the sake of brevity, only the 100-year design rain results are given here, in the paper. This scenario- design rain was supposed to fall at the end of dry period, on August 20, 1992 as well as at the end of wet period on July 27, 1997 uniformly distributed over both catchment areas of Vseminka and Drevnice. In dry pe riod 1992 on the Vseminka catchment, the impact of 100-year rain has increased direct runoff from 0.5 mm to 7.3 mm, and on the Drevnice catchment in the same period it varies from 4.0 mm to 32.0 mm. On both catchments this direct runoff increase has no catastrophical character, on the Vseminka catchment it was +6.8 mm, on the Drevnice catchment + 28.0 mm only. However, during the wet period this increase was more significant. On Vseminka was +45. lmm and on Drevnice even +81.3 mm. It means that on the latter case on a shallow-soil catchment with fully saturated active zone, almost 93 % of rainfall can leave the catchment in the form of direct runoff. This could undoubtedly bring a catastrophical event. The increase of subsidy to subsurface storage is also remarkable in all four scenario results. There the existing situation as well as the situation when the design rainfall of 87.1 mm has fallen are presented. Low values soil moisture deficit, SMD (and correspondingly high values of antecedent precipitation index, API) represent the situation when there is no more retention to replenish subsurface storages on the catchment. The situation on the Drevnice catchment is similar, even more sensitive for faster fluctuation because of the shallower active zone. Tables 4 and 5 and also Figures 2 and 3 give results of the comparative water balance component values. Figures representing the Drevnice catchment are similar but not attached to this paper because of its limited extent. CONCLUSIONS The data on hydrometeorological observations, soil, land use and on watershed management were checked for reliability and then spatially processed by ARC/INFO on the catchments of Vseminka and Drevnice. Then, after implementing the WBCM-5 model on that data, the resulting outputs were analysed and the following conclusions can be made: Comparing rainfall-runoff events in July 1997 on the experimental catchments, in spite of the relatively greater percentage of forestation in the Drevnice catchment (81.0 % against 48.2 % on Vseminka), the depth of rainfall (Drevnice about 30 % higher than those on Vseminka) and its intensity are more important factors than land use in a flood consideration. Water balance modelling such as implementing only "upper layers" (i.e. active zones) on catchments can easily quantify daily soil moisture deficit (SMD) values which signalize a degree of water saturation In dry period of the year 1992 was direct runoff from the scenario-design rainfall of 87.1 mm on the experimental catchment Vseminka 13 % and on Drevnice 32 %. In wet period of the year 1997 was direct runoff from this rainfall on Vseminka 52 % and on Drevnice as much as 93 %. A degree of situation expressed either by API or SMD is, besides the design rainfall characteristics its depth, duration and intensity, very substantial in a direct runoff formation process. Changes in subsurface water subsidy are logical, the great recharge improvements in dry periods or both catchments are remarkable, their excess in wet periods is not dangerous as groundwater flow from their subsurface storages is supposed as not too rapid. Small changes of actual evapotranspiration depend on an active zone and whole unsaturated soil water content. These changes are not essential. The WBCM model can be used only for individual water balance components quantification but also for an identification of the potential harm of floods and/or droughts. ACKNOWLEDGEMENT Support from the Scientific Grant Agency of the Mi- lity, retention and accumulation capacities of water in nistry of Education, Czech Republic, Project No. MSM the landscape " is gratefully acknowledged. 414100008 "Possibilities to increase ecological stabi |