ZNEČIŠTĚNÍ OVZDUŠÍ NA ÚZEMÍ ČESKÉ REPUBLIKY V ROCE 2009
Český hydrometeorologický ústav - Úsek ochrany čistoty ovzduší



II.2 Mapování územního rozložení charakteristik kvality ovzduší

Česká legislativa přejímá ze směrnic EU obecné přístupy pro stanovení úrovně kvality ovzduší a případného překročení stanovených imisních limitů v zónách pro účely řízení kvality ovzduší s cílem dosáhnout ve stanovených termínech kvality ovzduší vyhovující imisním limitům a cílovým imisním limitům. Stanovuje, že posuzování úrovně znečištění ovzduší se provádí měřením v aglomeracích a v zónách, kde úroveň znečištění ovzduší znečišťující látkou dosahuje nebo přesahuje horní mez pro posuzování, a měřením v oblastech, kde úroveň znečištění ozonem překračuje dlouhodobé imisní cíle (během posledních 5 let); dále modelováním nebo odborným odhadem v oblastech, kde úroveň znečištění ovzduší znečišťující látkou nepřesahuje dolní mez pro posuzování; a dále kombinací měření a modelování v oblastech, kde úroveň znečištění ovzduší dosahuje nebo přesahuje dolní mez pro posuzování a současně je nižší než horní mez pro posuzování.

Stanovení úrovní znečištění má pokrývat celou hodnocenou oblast, nejen pouhé okolí monitorovací stanice. Hodnocení kvality ovzduší v zónách a aglomeracích – zejména identifikace a vymezení oblastí případného překračování imisních limitů na základě měření – je tak problém územního odhadu rozložení sledované míry znečištění ovzduší a spočívá v tom, jak zobecnit „bodová“ měření při dané hustotě a rozložení monitorovacích stanic a akceptovatelné chybě odhadu na celé hodnocené území. Zvýšení prostorového pokrytí měření lze dosáhnout provedením ověřovacích měření. Směrnice pro kvalitu venkovního ovzduší a souhlasně i národní legislativa nevyžaduje měření jako jediný nástroj určování úrovní v zóně a předpokládá, v závislosti na úrovni znečištění, využití modelování, případně odborných odhadů a jejich kombinaci. Předností modelů je, že ve srovnání s bodovými měřeními lépe vystihují pokrytí posuzovaného území, nicméně jsou všeobecně považovány za méně přesné než měření. Modelováním se především myslí kauzální modely rozptylu a transportu včetně chemických transformací znečišťujících příměsí. Nezastupitelnou roli však mají i empirické, matematicko-statistické modely odhadu časového či prostorového rozložení imisních charakteristik.

Mapy imisních charakteristik a atmosférické depozice jsou vytvářeny na základě propojení a integrace systému GIS, relační databáze naměřených imisí a chemického složení atmosférických srážek ISKO a modelových výpočtů zejména na základě emisí, což je umožněno využitím výkonné výpočetní techniky a moderního softwarového vybavení. Důležitou roli hraje i doplňování a korekce objektivních výpočtů na základě odborného odhadu pracovníků pověřené instituce. Těmito postupy je možné velmi kvalitně hodnotit znečištění ovzduší a vytvořit adekvátní vizualizace a prezentace vhodné pro uživatele, ať už jsou jimi pracovníci státní správy či odborná i laická veřejnost.

Kromě využití výsledků z přímého měření koncentrací znečišťujících látek jsou využity i výsledky modelování. Pro území ČR se používá gaussovského disperzního modelu SYMOS 97, který počítá koncentrace na základě podrobných emisních inventur a meteorologických podmínek relevantních pro období hodnoceného kalendářního roku. Pro účely výpočtu modelu je celé území státu rozděleno zejména na základě geomorfologického členění do 47 oblastí, které mají rozdílné meteorologické podmínky. Každou z oblastí charakterizuje příslušná větrná růžice, která je jedním ze vstupů do modelu. Do výpočtu jsou zahrnuty poslední dostupné informace o zdrojích znečišťování z emisní databáze ISKO a informace o emisích z liniových zdrojů. Kromě zdrojů v ČR jsou do výpočtu pravidelně zahrnovány i dostupné informace o emisích ze zahraničních zdrojů, které mají nezastupitelnou úlohu zejména při výpočtu koncentrací v pohraničních oblastech, mohou se však uplatnit i v regionech od hranic vzdálenějších.
Kromě rozptylového modelu je v některých případech (např. pro přízemní ozon) aplikován empirický model za využití veličin vykazujících regresní závislost s naměřenými koncentracemi (jako nadmořská výška).
Jedním z důležitých předpokladů při tvorbě polí rozložení koncentrací je pečlivý výběr měřicích stanic zahrnutých do hodnocení z hlediska jejich využívání, klasifikace a reprezentativnosti.

Při tvorbě mapových podkladů imisního a depozičního zatížení území ze znečišťování ovzduší jsou při odhadech polí imisních a depozičních charakteristik na podkladě staničních měření využívány geostatistické postupy a nástroje mapové algebry geografického informačního systému (GIS).

Pro vytvoření výsledných map je použita asimilace naměřených a modelových (resp. dalších doplňkových) dat pomocí lineární regresní závislosti příslušných veličin (měření a modelu, resp. nadmořské výšky) s následnou interpolací reziduí této regrese. Při interpolaci je aplikována metoda IDW (interpolace pomocí váženého průměru hodnot naměřených v okolí interpolovaného bodu), či interpolační metoda kriging (interpolace pomocí váženého průměru hodnot naměřených v okolí interpolovaného bodu, kde váhy jsou funkcí statistické struktury pole imisních nebo depozičních charakteristik). Obě tyto interpolační metody umožňují odhadnout hodnotu sledované charakteristiky v libovolném místě pole. Za předpokladu, že vyšetřované pole je statisticky homogenní [1], je odhad získaný metodou kriging optimální v tom smyslu, že je nestranný a jeho střední kvadratická chyba je minimální. Programové vybavení geografického informačního systému umožňuje při aplikaci interpolační metody kriging vypočítat chyby odhadu interpolovaných hodnot. Hodnoty těchto chyb ukazují mimo jiné na účelnost zahuštění staniční sítě a naopak.
Při konstrukci map u většiny znečišťujících látek je výše popsaná metodika aplikována zvlášť pro městská a zvlášť pro venkovská území (s využitím klasifikace stanic), výsledná mapa je vytvořena sloučením pomocí gridu populační hustoty.

Při konstrukci prostorového rozložení koncentrací PM10 byl použit empirický model, který kombinuje rozptylové modely SYMOS, CAMx, evropský model EMEP a nadmořskou výšku s naměřenými koncentracemi na pozaďových stanicích za pomoci metodiky vyvíjené v rámci ETC/ACC [28]. Aplikace samotného modelu SYMOS by byla v případě znečištění PM10 nedostatečná, jelikož v modelu jsou započítány pouze emise z primárních zdrojů. Významný podíl ve znečištění ovzduší PM10 mají ovšem sekundární částice1 a resuspendované částice2 , které v emisích z primárních zdrojů zahrnuty nejsou, zohledňují je však právě modely EMEP a CAMx.

Základním určením míry reprezentativnosti je klasifikace stanic. Pozaďové stanice (typ „venkovská“ [rural], případně „městská pozaďová“ [urban background], či „předměstská pozaďová“ [suburban background]) s velkou reprezentativností (desítky kilometrů) jsou stanice ovlivňované pouze vzdálenými zdroji; pro vystižení lokálních poměrů jsou zohledňovány dopravní a průmyslové stanice (typ „dopravní“ [traffic] a „průmyslová“ [industrial]) s nejmenším poloměrem reprezentativnosti přímo ovlivňované místními zdroji.

Od roku 1994 se začaly vytvářet základní geografické a tematické vrstvy ve standardizované projekci (konformní Gauss-Krügerovo zobrazení). Z podkladů digitálních vrstev DMÚ 200, DMR-2 a nově DMÚ25 byly vytvořeny základní vrstvy pro geografický informační systém: orografie, nejvýznamnější vodní toky a vodní plochy, sídla, hranice okresů – správní zřízení, silniční síť a rostlinný kryt. Nejnověji byly vytvořeny základní vrstvy administrativního členění z geografických podkladů poskytnutých Českým statistickým úřadem.


1Sekundární částice jsou částice vzniklé reakcemi mezi oxidy síry a dusíku a mezi amoniakem a organickými sloučeninami v ovzduší. Definice dle: http://glossary.eea.eu.int/EEAGlossary/S/secondary_particles, viz též [22].
2Resuspendované částice jsou částice původně usazené na zemský povrch zvířené větrem nebo pohybem vozidel.