Souhrn
Cílem této studie je regionalizace České republiky z hlediska rizika výskytu sucha v jednotlivých povodích. Regionalizace České republiky byla provedena pro sadu 133 (mezi-)povodí pro období 1901–2015. Primárním indexem popisujícím sucho v rámci této studie jsou nedostatkové objemy vzhledem k prahu odpovídajícímu 20% kvantilu měsíčních průtoků. Na základě srážek, výparu, celkového a základního odtoku a hydrogeologických rajonů byla pomocí shlukové analýzy provedena regionalizace České republiky z hlediska chování v době sucha. Tato regionalizace byla následně expertně revidována. Nedostatkové objemy v simulaci modelu Bilan byly vyčísleny a byl vytvořen statistický model pro odhad N-letých nedostatkových objemů. Charakteristiky sucha v simulaci modelu Bilan pro povodí s dostupnými pozorovanými daty i výsledky statistického modelu byly úspěšně validovány.
Úvod
Odhad charakteristik sucha, včetně regionalizace jeho výskytu, je zatížen značnou nejistotou. Tato nejistota je, oproti analýze extrémních srážek, navíc znásobena skutečností, že sucho se nevyskytuje každý rok a efektivní délka standardně používaných 30–50letých časových řad je ve skutečnosti přibližně poloviční, významnější sucha se vyskytují pouze v malém zlomku let. Ve studii popsané v tomto článku jsou proto využity dvě strategie pro snížení nejistot odhadů charakteristik sucha.
První spočívá ve využití dlouhých časových řad srážek a teploty (1901–2015), které jsou dostupné v měsíčním kroku a prostorovém rozlišení 0,5° pro celý svět. Tyto řady, opravené o systematické chyby, slouží jako vstupy do modelu Bilan a výsledná simulace odtoku je analyzována z hlediska charakteristik sucha. Druhá strategie je založena na využití regionální frekvenční analýzy. Tato metoda předpokládá, že rozdělení normovaných extrémů (např. maximálních ročních deficitních objemů) je v homogenních oblastech stejné pro všechna povodí. Normovaným rozdělením extrémů se rozumí rozdělení hodnot vydělených normovacím faktorem (např. průměr nebo jiný parametr polohy daného rozdělení). Kombinací dat z různých povodí tak lze částečně nejistoty odhadu redukovat.
Následující kapitola stručně popisuje použitá data, studované charakteristiky nedostatkových objemů a jejich validaci. Další kapitola popisuje regionalizaci ČR z hlediska rizika sucha a poslední kapitola prezentuje odhady N-letých nedostatkových objemů.
Obr. 1. Porovnání charakteristik sucha pozorovaného a modelovaného průtoku; na ose x je vynesena pravděpodobnost z distribuční funkce jednotlivých charakteristik, na ose y jsou uvedeny hodnoty jednotlivých charakteristik sucha, barevné plochy jsou rozsahy 25% a 75% kvantilů hodnot charakteristik ze souboru povodí
Fig. 1. Comparison of drought characteristics for observed and simulated runoff; the horizontal axis shows the probability from the distribution of the individual characteristics, the vertical axis indicate the values of the drought characteristics, the colored polygons correspond to the range between the 25th and 75th percentile of the value of drought characteristics in the set of catchments
Charakteristiky sucha
Pro účely studie je sucho definováno na základě nedostatkových objemů odvozených ze simulované hydrologické bilance pro sadu mezipovodí pro období 1901–2015. Odvozené časové řady měsíčních srážek a teplot na plochy mezipovodí jsou založeny na kombinaci HadCRU-TS3.21 [1], staničních dat a gridovaného datasetu srážek a teploty [2]. Gridovaná data byla převedena na plochu povodí pomocí váženého průměru (podle velikosti průniku plochy povodí a příslušných grid boxů), následně byla u teploty provedena korekce na nadmořskou výšku na základě porovnání nadmořské výšky povodí a průměrné nadmořské výšky grid boxů a korekce srážek na základě vrstvy dlouhodobých průměrných srážek v období 1981–2010. Pomocí těchto časových řad a parametrů modelu Bilan získaných kalibrací na pozorovaných datech byla simulována hydrologická bilance pro celé období 1901–2015.
Nedostatkové objemy jsou definovány jako suma doplňků měsíčních odtoků souvisle pod 20% kvantilem rozdělení odtoku. Mezi základní charakteristiky nedostatkových objemů lze zařadit:
- velikost nedostatkového objemu, D [mm nebo m3],
- doba trvání deficitní události, L [měsíce],
- intenzita sucha, I = D/L [mm/měsíc nebo m3/měsíc],
- poměr deficitního objemu k průměrnému měsíčnímu odtoku, rD [-],
- relativní intenzita, rI = rD/L [t-1].
Podle průměrných hodnot jednotlivých charakteristik sucha v tabulce 1 je zřejmé, že jsou na sobě jednotlivé charakteristiky lineárně závislé. Také lze pozorovat postupný růst jednotlivých charakteristik v průběhu první (1901–1930) až třetí periody (1961–1990), v poslední periodě (1991–2015) lze pozorovat jejich pokles.
Tabulka 1. Porovnání průměrných hodnot charakteristik sucha podle jednotlivých period na základě modelovaného celkového odtoku (D [mm] velikost nedostatkových objemů, I [mm/měsíc] intenzita sucha, L [měsíc] délka deficitního období, rD [-] poměr deficitního objemu k průměrnému měsíčnímu odtoku, rI [t-1] relativní intenzita)
Table 1. Average drought characteristics for different time periods based on simulated runoff (D [mm] deficit volume, I [mm/month] drought intensity, L [month] length of the deficit period, rD [-] ratio of deficit volume to average monthly discharge, rI [t-1] relative intensity)
Validace deficitních objemů
Validace deficitních objemů probíhala pomocí selekce uzavřených povodí se známým pozorovaným odtokem tak, aby bylo možné porovnat jednotlivé charakteristiky sucha vypočítané na základě modelovaného a zároveň pozorovaného průtoku. Validaci bylo možno provést na 65 povodích (ze 133).
Jak je možné vidět v tabulce 2, všechny charakteristiky sucha vypočítané na základě simulovaného odtoku jsou pouze mírně vyšší oproti charakteristikám získaným pomocí pozorovaného odtoku.
Obrázek 1 prezentuje porovnání charakteristiky jednotlivých deficitních událostí. Je zřejmé, že model Bilan simuluje důležité charakteristiky sucha dostatečně přesně.
Tabulka 2. Porovnání celkových průměrů charakteristik sucha pro vybraných 65 povodí na základě pozorovaného a simulovaného odtoku (D [mm] velikost nedostatkových objemů, I [mm/měsíc] intenzita sucha, L [měsíc] délka deficitního období, rD [-] poměr deficitního objemu k průměrnému měsíčnímu odtoku, rI [t-1] relativní intenzita)
Table 2. Comparison of average drought characteristics in simulated and observed data for 65 catchment with available data (D [mm] deficit volume, I [mm/month] drought intensity, L [month] length of the deficit period, rD [-] ratio of deficit volume to average monthly discharge, rI [t-1] relative intensity)
Regionalizace ČR z hlediska výskytu sucha
Ze všech fází regionalizace je identifikace homogenních oblastí zpravidla nejobtížnější a vyžaduje velké množství subjektivních rozhodnutí. Pro identifikaci homogenních oblastí byl využit algoritmus K-means [3], jehož cílem je nalezení shluků v datasetu. Vstupem algoritmu je množina bodů, které jsou definované souřadnicemi v n-rozměrném prostoru, a číslo k, určující požadovaný počet shluků. Všechny shluky jsou reprezentované svými středy a každý bod náleží do shluku, jehož střed je mu nejblíže. Souřadnice středů se určují iterativním způsobem, tento iterativní algoritmus neustále zmenšuje chybu, definovanou jako součet vzdáleností všech bodů od středů svých shluků, a spěje tak k lokálně-optimálnímu řešení. Shluková analýza pro povrchový odtok proběhla se škálovanými daty průměrných odtoků a rozdílů odtoků a potenciální evapotranspirace (tj. ve dvourozměrné metrice [odtok; odtok – PET]), algoritmus proběhl celkově desetkrát, pokaždé se startem v jiné, náhodné pozici. Do padesáti iterací každý běh konvergoval k lokálně-optimálnímu řešení, výsledné shluky a jejich středy jsou zobrazeny na obr. 2 a výsledná regionalizace pak na obr. 3.
Obr. 2. Výsledné shluky (barevně) a jejich středy (lokální optima)
Fig. 2. Resulting clusters (colors) and their centers (bold dots)
Regionální extremální model pro nedostatkové objemy
Odhad N-letých nedostatkových objemů je zatížen značnou nejistotou, vyplývající mj. z délky dostupných řad. Podobně jako pro extrémní průtoky a srážky lze i pro nedostatkové objemy využít tzv. index-flood metodu. Předpokladem je, že nedostatkové objemy v každém (mezi-)povodí v rámci předem definované oblasti mohou být normovány tak, že rozdělení těchto normovaných nedostatkových objemů je v dané oblasti stejné. Normovací faktor, který je určen pro jednotlivá (mezi-)povodí, je zpravidla označován index-flood, stejně jako tato metoda [4]. Primárně uvažujeme, že nedostatkové objemy v jednotlivých povodích lze popsat pomocí zobecněného rozdělení extrémních hodnot (Generalized Extreme Value Distribution – GEV). Rozdělení je korigováno vzhledem k výskytu let bez sucha na základě vztahu, který uvádí Engeland [5].
Výsledkem optimalizace modelu jsou normovací faktory pro jednotlivá povodí a regionální kvantilová funkce. Detaily modelu popisuje Hanel [6]. Regionální kvantilové funkce spolu s normovanými nedostatkovými objemy pro jednotlivá povodí uvádí obr. 4. Pravděpodobnosti p jsou transformované pomocí vztahu –log(-log (p)). Tato transformace zajišťuje lepší zobrazení extrémů s vysokou dobou opakování.
Na obr. 5 jsou uvedeny odhady 2, 5, 10, 20 a 50letého nedostatkového objemu [mm] pro celkový odtok. Je zjevné, že nedostatkové objemy jsou do jisté míry korelované s celkovým odtokem. V relativním vyjádření vzhledem k průměrnému celkovému odtoku (obr. 6) nedostatkové objemy rostou s dobou opakování (platí i v absolutní hodnotě) a zřetelně je vidět větší dopad na povodí v oblasti 3 (tj. v oblasti s nejvyšším rizikem výskytu sucha).
Obr. 3. Výsledky regionalizace
Fig. 3. Resulting regionalization
Obr. 4. Regionální kvantilové funkce pro jednotlivé oblasti (tlustá čára) spolu s normovanými kvantilovými funkcemi pro pozorovaná data z jednotlivých (mezi-)povodí
Fig. 4. Regional quantile functions for individual regions (bold line) together with standardized quantile functions for the observed data from individual catchments
 |
 |
 |
 |
 |
Obr. 5. Mapa N-letých nedostatkových objemů pro (simulovaný) celkový odtok pro 2, 5, 10, 20, 50 let
Fig. 5. Map of N-year deficit volumes for (simulated) total runoff for 2, 5, 10, 20, 50 years
Závěr
Hydrologický model Bilan byl použit pro simulaci hydrologické bilance 133 mezipovodí v ČR pro období 1901–2015. Bylo ukázáno, že simulované charakteristiky sucha odpovídají dobře charakteristikám pozorovaným (pro povodí, pro která jsou dostupná data). Území ČR bylo rozděleno do tří homogenních oblastí označovaných „bez rizika“ (horské polohy), „mírně ohrožené“ (Vysočina a střední polohy) a „ohrožené“ (nížiny – zejména jižní Morava a střední Čechy). Byly odhadnuty parametry statistického modelu umožňujícího vyčíslení N-letých nedostatkových objemů. Z hlediska rozdělení nedostatkových objemů se jednotlivé oblasti liší zejména variabilitou sucha – oblasti označované „bez rizika“ mají poměrně vyrovnaný režim a nedostatkové objemy v jednotlivých suchých obdobích se příliš neliší. Naopak v oblastech označovaných jako „ohrožené“ je velký rozdíl mezi různě extrémními suchy – tj. rozdíl mezi 2letým a 10letým suchem je zde podstatně větší než v ostatních oblastech.
Dalšími kroky v regionalizaci hydrologického sucha v ČR je zjemnění prostorového měřítka použitého pro regionalizaci a odhad N-letých nedostatkových objemů až na povodí útvarů povrchových vod, vyhodnocení jednotlivých historických případů sucha, prodloužení rekonstrukce dále do minulosti, simulace podmínek ovlivněných klimatickou změnou a vyhodnocení nejistot.
Obr. 6. Poměr deficitního objemu k průměrnému měsíčnímu odtoku pro doby opakování 2, 5, 10, 20, 50 let
Fig. 6. Ratio of deficit volume to mean monthly runoff for 2, 5, 10, 20, 50 years return periods
Poděkování
Studie byla provedena v rámci úkolu Činnosti k podpoře výkonu státní správy v problematice sucho v roce 2016 a vnitřních grantů Fakulty životního prostředí ČZU č. 20164230 (FS) a 20174227 (VM).
