Methoden zur Berechnung der Deposition von Schwermetallen im Umfeld der Affi

Kapitelende

Modellierung

Die Abteilung Luftuntersuchungen in der Umweltbehörde misst die Staubdeposition, die sich in einem Glasbecher auf einer Stange in 1.5 m Höhe innerhalb eines Monats sammelt. Welchen Umkreis repräsentiert eine an einem Punkt gesammelte Probe? In der alten "Technischen Anleitung Luft" (TA Luft, Bestandteil des Bundesimmissionsschutzgesetzes) wurde der Mittelwert der vier Eckpunkte eines Planquadrats für die Belastung der Fläche gesetzt. In diesem Sinne würde man die Fracht auf den Untersuchungsraum aus dem Mittel der 25 Messpunkte berechnen. Dies ist in unserm Fall nicht realistisch, weil die Staubkonzentration um eine niedrige Quelle (Erzumschlagplatz, Dachreiter) herum rasch abfällt. Der dreckigste Platz liegt unter der Quelle. Die Emissionen aus hohen Quellen, den Schornsteinen, schlagen sich erst in größerer Entfernung nieder - der sauberste Platz liegt am Fuss des Schornsteins - und werden vom Programm der BUG kaum erfasst. Zwischen den Messwerten an den Punkten 452a und 452b zwischen den Werksteilen und denen an den Rändern des Untersuchungsraumes liegen zwei Größenordnungen. Interpolierte man linear zwischen dem höchsten Messwert von bespielsweise 500 ug/m2*d und einem 2km entfernten Punkt mit 10 ug/m2*d, könnte man in der Mitte den Wert 255 ug/m2*d erwarten, findet jedoch typischerweise nur ca. 50 ug/m2*d. Über dem Untersuchungsgebiet liegt eine Staubglocke mit dem Maximum auf dem Werksgelände, die zu den Rändern flach ausläuft. Um diese "Glocke" realistisch zu berechnen, wurden die Punktwerte in Rasterschritten von 200m interpoliert mit einem Computerprogramm "Idrisi", einem geografischen Informationssystem (Idrisi for Windows, version 2), wobei das "Gewicht" eines Punktes mit der Entfernung mit einem Exponent X abnahm. Punkte ohne Messwert (weil im Messjahr nicht beprobt) wurden nicht einbezogen.

Am Beispiel der Arsen-Deposition 1987 wurden mehrere Varianten durchgespielt mit den Exponenten 1, 2, 3 und 4.

Das Programm Idrisi füllt von jedem Punkt die Pixel mit Zwischenwerten zu den nächstgelegenen 6 Punkten und vereint dann diese Resultate zu einer Fläche. Dadurch werden die ursprünglichen Messwerte auf ihrem zugehörigen Pixel durch einen für die Interpolation am besten passenden Wert ersetzt. Legt man Schnitte durch die "Staubglocke" und trägt die Pixelwerte entlang des Schitts auf, erhält man für den Schnitt von Nordwest nach Südost folgende Variationen:

Ohne flächendeckendesInterpolationsverfahren erhielte man die rosafarbene Verbindungslinie zwischen den Messwerten auf Pixel 1 (Messpunkt Nr. 397), Pixel 6 (Messpunkt Nr. 423), Pixel 11 (Messpunkt Nr. 452b) und Pixel 16 (Messpunkt Nr. 483). Die resultierende Fracht, das Volumen unter der "Staubpyramide", ist bei einem solchen Ansatz am höchsten. Die Interpolation mit Exponent 1, blaue Linie, drückt das Maximum, "verschmiert" die Belastung aber weiträumig. Weiterhin wurde mit den Exponenten 2 (grün), 3 (rot) und 4 (violett) gerechnet. Weil mit dem Exponent 4 die ursprünglichen Messwerte am besten erhalten bleiben und die Form der "Staubglocke" plausibel simuliert wird, wurde er für die weiteren Rechnungen gewählt. Die Fracht ist bei dieser Variante am geringsten.

Interpolation	ug/m2*d	kg/a
arithm. Mittelwert	65,8	600,5
Exponent 1	57,8	527,8
Exponent 2	52,3	476,8
Exponent 3	48,8	445,6
Exponent 4	47,0	428,9
Exponent 4 + Filter	43,8	400,1

In einem einfachen Median-Filter Durchgang wurde das Interpolationsbild geglättet, um die bei einer Modellierung über 25 Messpunkte auftretenden Artefakte zu eliminieren. Das ist zulässig, weil jedes Modell nur eine mehr oder weniger verzerrte Karikatur der Wirklichkeit ist. Erstes Ziel der Berechnungen war, vergleichbare Zeitreihen zu erzeugen, zweitrangig war die Bestimmung der absoluten Frachten.

Aussreisser

Bei der monatlichen Beprobung werden gelegentlich sehr hohe Werte gefunden, die den Jahresdurchschnitt deutlich nach oben treiben. Die Behörde hat hierfür keine plausible Erklärung. Bei unseren Modellen der Staubglocke wurden die Ausreisser nicht berücksichtigt. Nichtsdestoweniger ist es zu fordern, dass Behörde und Affi sich um eine Klärung des Sachverhalts bemühen, denn die Extremwerte sind nicht statistisch auf alle Messpunkte und Parameter gestreut. In der Tabelle sind alle Ausreisser dargestellt. Der Parameter ist mit Farbe unterlegt (As magenta, Cu orange, Cd gelb, Pb grau). Der Zahlenwert gibt das Jahresmittel am Messpunkt unter Berücksichtigung des Ausreissers an, in Klammern ohne diesen. Beim Parameter Arsen findet man im Vergleich zu den anderen häufiger Extremwerte. Das Jahr 1996 fällt durch Ausreisser an mehreren Punkten und mit allen vier Schadstoffen am Punkt 483 (südöstlich vom Werk Ost) auf.

Frachten

Aus dem Inhalt der „Staubglocke“ lässt sich die Fracht, die auf die 25 km2 Untersuchungsfläche nieder rieselt, einfach berechnen. Zusätzlich wurde der Werksbereich und die ummittelbare Umgebung ausgestanzt. Der Verlauf der werksnahen und werksfernen Deposition gibt Hinweise auf die Quelle der Emission und ihre Höhe.
as
cd
cu
pb

Werksbereich und Wasserflächen im Untersuchunggebiet

Soweit vorhanden wurden Emissionsmessungen der Behörde und der Affi den Tabellen angefügt. Es wurden offenbar nicht für jedes Jahr solche Untersuchungen durchgeführt bzw. veröffentlicht, und bei den Angaben der Affi aus jüngerer Zeit fehlt die Aufschlüsselung nach Quellhöhen. Generell sind die Emissionsfrachten, die an den Quellen gemessen wurden, höher als die durch die Deposition in der Umgebung bestimmten Frachten. Das hat mehrere Ursachen:

• die Messmethoden sind sehr verschieden

• das Abluftvolumen an den beprobten Quellen, wie Dachreitern, ist schwer zu bestimmen und wird wahrscheinlich überschätzt

• an den Quellen werden auch grobe Staubkörner erfasst, die noch auf dem Werksgelände niedergehen, und nicht mehr zu den (Fein)Staubsammlern der BUG gelangen. Stünden mitten in den beiden Werksteilen noch jeweils Staubsammler, würde auch das Depositionsmessprogramm der BUG höhere Frachten ergeben.

• die Frachten aus hohen Quellen (100 – 159 m) schlagen ausserhalb des Untersuchungsraumes der BUG sich nieder.

Beide Messverfahren sollten abgestimmt aufeinander fortgesetzt werden, um die Wirksamkeit von Emissionsminderungsmaßnahmen der Affi besser beurteilen zu können.

Tabelle 1: Übersetzer

Tabelle 2: Arsen

Tabelle 3: Cadmium

Tabelle 4: Kupfer

Tabelle 5: Blei

Als „Nebeneffekt“ dieser Arbeit wurde berechnet, wieviel Schadstoffe aus der Affi direkt in Gewässerflächen eingetragen werden, um der BUG bei der Wasserrahmenrichtlinie auf die Sprünge zu helfen. Die BUG nämlich hielt es nicht der Mühe wert, sich Gedanken über diffuse Einträge in Gewässer zu machen. Allein der direkte Aufschlag auf Wasserflächen übertrifft die wie auch immer „gemessenen“ Einleitungsfrachten der BUG, die sie als „signifikant“ deklariert. Man muss jedoch noch hinzu rechnen, welche Mengen auf befestigte Flächen sich ablagern und mit dem nächsten Regen in Gewässer gespült werden (wobei auf dem Affi-Werksgelände das Regenwasser erst durch eine Reinigungsanlage läuft).

Schwermetallbelastung 2001 durch direkte Einleitung in (Angaben BUG, Zwischenberichte WRRL) und Deposition auf Gewässer (Berechnung RdE)

stoff	Einleitungsfracht lt. BUG kg/a	Eintrag durch Deposition auf Wasserflächen kg/a
as	28	27
cd	19	5
cu	531	1359
pb	47	371

Lesestoff:

• Freie und Hansestadt Hamburg, Umweltbehörde: Ergebnisse der Staubniederschlags- und Schwebstaubmessungen mit anschließender Bestimmung von Staubinhaltsstoffen im Zeitraum 1993 bis 1997 in Hamburg; Hamburger Umweltberichte 54/97

• Freie und Hansestadt Hamburg, Umweltbehörde: Bericht über die Emissionen der Norddeutschen Affinerie und die Immissionen in ihrem Umfeld; Hamburg, Sept. 1999

• Norddeutsche Affinerie AG: Menschen und Umwelt; Hamburg, Jan. 2004

• EPER (European Pollutant Emission Register)

• Idrisi

Kapitel oben
erstellt März 2004