Badania wykonane w latach 2005–2008 obejmowały studium materiałów publikowanych i archiwalnych, wyznaczenie miejsc pobierania próbek gleb na mapach topograficznych w skali 1:10 000, pobranie próbek i pomiary współrzędnych w miejscach ich lokalizacji, analizy chemiczne próbek, utworzenie baz danych terenowych i laboratoryjnych, opracowanie wektorowego podkładu topograficznego, obliczenia statystyczne, opracowanie map geochemicznych i mapy geologicznej oraz interpretację wyników. Kolejność prac ilustruje załączony schemat (fig. 1).
PRACE TERENOWE
Próbki gleb pobierano w regularnej siatce 250×250 m (16 próbek/km²). Łącznie pobrano próbki gleb z 1289 miejsc. W każdym miejscu pobierano gleby z głębokości: 0,0–0,3 m i 0,8–1,0 m (lub z mniejszej głębokości w przypadku płytszego zalegania skał macierzystych). Próbki gleb (o masie ok. 500 g), pobierane za pomocą ręcznej sondy o średnicy 60 mm, umieszczano w płóciennych woreczkach opatrzonych odpowiednimi numerami i wstępnie suszono na drewnianych paletach w magazynie terenowym.
Próbki osadów i wód powierzchniowych pobierano z rzek, strumieni, rowów melioracyjnych, kanałów, sadzawek i stawów. Odległość między miejscami pobierania próbek na ciekach wynosiła ok. 250 m. Próbki osadów o masie ok. 500 g (i możliwie najdrobniejszej frakcji) pobierano z brzegów zbiorników za pomocą czerpaka i umieszczano w pojemnikach plastikowych o pojemności 500 ml, opisanych odpowiednimi numerami.
Próbki wód powierzchniowych pobierano z tych samych miejsc, z których pobrano próbki osadów wodnych. Przewodność elektryczną właściwą wód (EC) i ich odczyn (pH) mierzono w terenie. Do pomiarów EC stosowano konduktometr z automatyczną kompensacją temperaturową, przyjmując temperaturę referencyjną 25°C. Próbki wód były filtrowane w terenie przez filtry Millipore 0,45 µm, a po napełnieniu butelek o objętości 30 ml, zakwaszane kwasem azotowym. Butelki opisywano odpowiednimi numerami.
Miejsca pobrania wszystkich próbek zaznaczano na mapach topograficznych w skali 1:10 000 i opisywano odpowiednimi numerami.
Położenie miejsc pobierania próbek określono techniką GPS (Global Positioning System), stosując urządzenie wyposażone w zewnętrzną antenę oraz komputer umożliwiający obok pomiarów współrzędnych rejestrację dodatkowych informacji (wartości pH i EC wód, danych o zabudowie i użytkowaniu terenu oraz o charakterze litologicznym próbek). Pomiar współrzędnych rejestrowano z dokładnością ±2–10 m. Przed wyjazdem w teren do pamięci urządzenia GPS wprowadzano sieć współrzędnych miejsc pobierania próbek. Kolejne miejsca pobierania próbek były wyszukiwane w terenie metodą nawigacji satelitarnej. Dla większego bezpieczeństwa wszystkie dane terenowe notowane były również na specjalnie przygotowanych kartach (fig. 2).
PRACE LABORATORYJNE
Przygotowanie próbek. Po przewiezieniu do laboratorium próbki gleb były suszone w temperaturze pokojowej i przesiewane przez sita nylonowe o oczkach 2 mm (ISO 11464). Każdą próbkę gleby z głębokości 0,0–0,3 m po przesianiu i kwartowaniu dzielono na trzy podpróbki: jedną przeznaczoną do analizy chemicznej, drugą do analizy granulometrycznej i trzecią archiwalną. Każdą próbkę gleby z głębokości 0,8–1,0 m po przesianiu i kwartowaniu, dzielono na dwie podpróbki: jedną przeznaczoną do analizy chemicznej, drugą archiwalną (fig. 1). Próbki gleb przeznaczone do analiz chemicznych ucierano do frakcji <0,06 mm w agatowych młynach kulowych.
Próbki osadów wodnych były suszone w temperaturze pokojowej, a następnie przesiewane przez sita nylonowe o oczkach 0,2 mm. Frakcję <0,2 mm po kwartowaniu dzielono na dwie podróbki: jedną przeznaczoną do analizy chemicznej, drugą archiwalną (fig. 1).
Wszystkie próbki archiwalne zmagazynowano w Państwowym Instytucie Geologicznym-Państwowym Instytucie Badawczym w Warszawie.
Analizy chemiczne wykonano w laboratorium chemicznym Państwowego Instytutu Geologicznego-Państwowego Instytutu Badawczego w Warszawie.
Roztwarzanie próbek gleb i osadów wodnych przeprowadzono w wodzie królewskiej (1 g próbki do końcowej objętości 50 ml) przez 1 godz. w temp. 95°C w termostatowanym bloku aluminiowym.
Oznaczenia zawartości Ag, Al, As, Ba, Ca, Cd, Co, Cr, Cu, Fe, Mg, Mn, Ni, P, Pb, S, Sr, Ti, V i Zn w glebach oraz osadach wodnych wykonano za pomocą spektrometrii emisyjnej ze wzbudzeniem plazmowym (ICP-AES). Analizy Hg w próbkach gleb i osadów wodnych wykonano metodą spektrometrii absorpcji atomowej z techniką zimnych par (CV-AAS) z systemem przepływowym. Odczyn gleb oznaczono pehametrycznie w wyciągach wodnych, a zawartość w nich węgla organicznego metodą kulometryczną. Oznaczenia zawartości B, Ba, Ca, Cr, Fe, K, Mg, Mn, Na, P, SiO₂, SO₄, Sr, Ti i Zn w wodach powierzchniowych przeprowadzono metodą ICP-AES, a zawartości Ag, Al, As, Cd, Cl, Co, Cu, Li, Mo, Ni, Pb, Rb, Sb, Tl i U metodą ICP-MS. Zestawienie metod analitycznych i granic oznaczalności pierwiastków przedstawiono w tabeli 1.
Poprawność wykonywanych oznaczeń chemicznych sprawdzano poprzez analizę próbek podwójnych (około 5% ogólnej liczby próbek), analizę materiałów odniesienia z atestowaną zawartością badanych pierwiastków (2% ogólnej liczby próbek) oraz analizę wewnętrznych próbek kontrolnych potwierdzających prawidłowe wykonywanie pomiarów instrumentalnych (5% ogólnej liczby próbek). Stosowano „ślepe próbki odczynnikowe” jak również „ślepe próbki proceduralne”. Czystość odczynników i naczyń kontrolowano za pomocą „ślepych próbek odczynnikowych”. „Ślepe próbki proceduralne” (sea sand extra pure Merck) stosowano do sprawdzania zanieczyszczeń możliwych do wprowadzenia podczas kolejnych etapów przygotowania próbki.
Dla próbek stałych precyzja oznaczeń wynosi ±10–15% (na podstawie analiz próbek podwójnych). Dla próbek wód precyzja wynosi ±10–20% (zależnie od zawartości pierwiastka).
Analizy granulometryczne gleb z głębokości 0,0–0,3 m wykonano w Zakładzie Hydrogeologii i Geologii Inżynierskiej Państwowego Instytutu Geologicznego-Państwowego Instytutu Badawczego w Warszawie, łącząc analizę sitową z metodą laserowego pomiaru wielkości cząstek.
Badania składu ziarnowego przeprowadzono metodą odbiegającą od analiz tradycyjnych (wykonywanych metodą sitowo-sedymentacyjną zgodnie z klasyfikacją przyjętą w podziałach międzynarodowych FAO i USDA). Porównania wyników składu ziarnowego uzyskane metodą sitowo-sedymentacyjną i metodą laserową, wskazują na znaczne różnice w zawartości poszczególnych frakcji (Issmer, 2000; Kasza, 1992). Bezpośrednie wykorzystanie wyników uzyskanych za pomocą tej metody nie umożliwia zatem podziału gleb według kryteriów gleboznawczych. Dane są natomiast bardzo pomocne przy interpretacji wyników badań geochemicznych.
Próbki przesiewano przez zestaw sit o oczkach 2 mm, 1 mm i 0,5 mm. Próbki niektórych gleb gliniastych rozdrabniano w moździerzu porcelanowym przed przesiewaniem. Ważono uzyskane frakcje: 2–1 mm, 1,0–0,5 mm i <0,5 mm. Pomiary wielkości ziaren we frakcji <0,5 mm przeprowadzono za pomocą laserowego miernika wielkości cząstek Analysette-22 firmy Fritsch. Wyniki analiz granulometrycznych (po przeliczeniu na udziały procentowe) przedstawiono na mapach z uwzględnieniem wybranych klas ziarnowych: 1,0–0,1 mm – frakcja piaszczysta, 0,1–0,02 mm – frakcja pylasta, <0,02 mm – frakcja ilasta (tabl. 4–6).
BAZY DANYCH I KONSTRUKCJA MAP GEOCHEMICZNYCH
Podkład topograficzny. Jako podkładu map geochemicznych w skali 1:25 000 użyto najbardziej aktualnej mapy topograficznej w skali 1:50 000 w układzie współrzędnych 1992, arkusz Jaworzno M-34-63-B (zapis wektorowy VMap L2). Mapa topograficzna zawiera następujące wektorowe warstwy informacyjne:
– rzeźba terenu,
– hydrografia (z podziałem na rzeki, strumienie, rowy i zbiorniki wód stojących),
– sieć komunikacji drogowej (z podziałem według klas),
– sieć kolejowa,
– zabudowa terenu (z podziałem na zabudowę wiejską, miejską oraz przemysłową),
– lasy,
– tereny przemysłowe (obiekty przemysłowe, wyrobiska kopalniane, hałdy i osadniki).
Mapa geologiczna. Dla ilustracji budowy geologicznej badanego obszaru wykorzystano Szczegółową mapę geologiczną Polski 1:50 000,arkusz Jaworzno (Kurek i in., 1999). Korzystając z plików ArcInfo, tworzono obrazy wektorowe poszczególnych elementów mapy geologicznej zakrytej, które połączono z podkładem topograficznym w formie mapy geologicznej w skali 1:25 000 (tabl. 1).
Bazy danych. Utworzono oddzielne bazy danych dla:
– gleb z głębokości 0,0–0,3 m,
– gleb z głębokości 0,8–1,0 m,
– osadów wodnych,
– wód powierzchniowych.
Bazy danych gleb zawierają: numery próbek, wyniki pomiarów współrzędnych w miejscach pobierania próbek, zapis obserwacji terenowych (rodzaj zabudowy, użytkowanie terenu, gatunek gleby, lokalizacja miejsca pobierania próbek – powiat, gmina, miejscowość), datę i nazwisko osoby pobierającej próbki oraz wyniki analiz chemicznych.
Bazy danych osadów wodnych i wód powierzchniowych zawierają: numery próbek, wyniki pomiarów współrzędnych w miejscach pobierania próbek, zapis obserwacji terenowych (rodzaj zabudowy, użytkowanie terenu, rodzaj zbiornika wodnego, gatunek osadu, lokalizacja miejsca pobierania próbek – powiat, gmina, miejscowość), datę i nazwisko osoby pobierającej próbki) oraz wyniki analiz chemicznych.
Obliczenia statystyczne. Wyniki zgromadzone w bazach danych służyły do wydzielania podzbiorów do obliczeń statystycznych według różnych kryteriów środowiska, na przykład zawartości pierwiastków w glebach terenów przemysłowych, glebach leśnych, glebach miejskich oraz w osadach wodnych i w wodach poszczególnych zbiorników, jak również do tworzenia map geochemicznych. Obliczenia parametrów statystycznych wykonano zarówno dla całych zbiorów, jak i podzbiorów dla gleb, osadów wodnych i wód powierzchniowych. Do obliczeń statystycznych przyjmowano zawartość równą połowie limitu detekcji danej metody analitycznej, w przypadku zawartości pierwiastków poniżej granicy oznaczalności tej metody. Wyliczano średnią arytmetyczną, średnią geometryczną, medianę oraz wartość minimalną i maksymalną. Dane te dla poszczególnych pierwiastków i wskaźników zestawiono w tabelach 2–5 oraz zamieszczono na mapach geochemicznych. Obliczone wskaźniki służyły do wyznaczania klas zawartości przy tworzeniu map geochemicznych.
Opracowanie map. Dla arkusza Jaworzno opracowano następujące mapy (tabl. 2–63):
– zabudowa terenu,
– użytkowanie terenu,
– zawartość węgla organicznego oraz frakcji piaszczystej, pylastej i ilastej w glebach z głębokości 0,0–0,3 m,
– pH w glebach z głębokości 0,0–0,3 m i 0,8–1,0 m,
– zawartość Ag, Al, As, Ba, Ca, Cd, Co, Cr, Cu, Fe, Hg, Mg, Mn, Ni, P, Pb, S, Sr, Ti, V i Zn w glebach z głębokości 0,0–0,3 m i 0,8–1,0 m oraz w osadach wodnych,
– pH i EC oraz zawartość Ag, Al, As, B, Ba, Ca, Cd, Cl, Co, Cu, Fe, K, Li, Mg, Mn, Mo, Na, Ni, P, Pb, Rb, Sb, SiO₂, SO₄, Sr, Ti, Tl, U i Zn w wodach powierzchniowych,
– klasyfikacja gleb z głębokości 0,0–0,3 m wskazująca właściwy sposób ich użytkowania (z wydzieleniem grup użytkowania gleb na podstawie Rozporządzenia..., 2002).
Zabudowę i użytkowanie terenu oraz klasyfikację gleb z zakresu głębokości 0,0-0,3 m w odniesieniu do sposobu ich użytkowania, przedstawiono w postaci map punktowych (tabl. 2, 3 i 63).
Do przedstawienia rozkładu klas ziarnowych (tabl. 4–6) oraz zawartości pierwiastków w glebach wybrano izoliniową (obszarową) metodę opracowania map ze względu na jej przejrzystość i czytelność. Izoliniowe mapy geochemiczne utworzono z użyciem programu Surfer, stosując metodę Inverse Distance to a Power. Klasy zawartości pierwiastków dobierano stosując postęp geometryczny.
Odczyn gleb (tabl. 7 i 8) przedstawiono według skali przyjętej w gleboznawstwie (z podziałem na gleby kwaśne, obojętne i zasadowe).
Mapy geochemiczne gleb opracowano dla zbioru wyników analiz chemicznych arkusza Jaworzno i arkuszy sąsiednich w skali 1:25 000, aby uniknąć niezgodności na granicy arkuszy. Z utworzonych map monopierwiastkowych wycinano arkusz Jaworzno i łączono z podkładem topograficznym.
Mapy geochemiczne osadów wodnych i wód powierzchniowych opracowano oddzielnie dla arkusza Jaworzno. Utworzono je w formie kartodiagramów kołowych, przypisując ich odpowiednie średnice do poszczególnych klas zawartości, ułożonych najczęściej w postępie geometrycznym.
Sporządzając mapę klasyfikacji gleb powierzchniowych (tabl. 63), wskazującą właściwy sposób ich użytkowania, wyniki badań geochemicznych odniesiono do wartości stężeń dopuszczalnych metali określonych w Rozporządzeniu Ministra Środowiska (2002), zgodnie z zaleceniem „glebę lub ziemię uznaje się za zanieczyszczoną, gdy stężenie co najmniej jednej substancji przekracza wartość dopuszczalną”.
Na podstawie zawartości każdego z badanych metali (wymienionych w Rozporządzeniu..., 2002) każda z badanych próbek gleb została zaklasyfikowana do grupy A, B lub C. W przypadku jednakowych granic wartości dopuszczalnych dla grupy A i B (wyznaczonych w Rozporządzeniu..., 2002 dla arsenu, baru i kobaltu) przyjęto zasadę klasyfikacji gleb do grupy A, korzystniejszej dla użytkownika i pozwalającej na wielofunkcyjne wykorzystywanie terenu.
Dla celów publikacji mapy geochemiczne opracowano poprzez połączenie parami na jednej tablicy mapy geochemicznej gleb z głębokości 0,0–0,3 m z mapą geochemiczną osadów wodnych oraz mapy geochemicznej gleb z głębokości 0,8–1,0 m z mapą geochemiczną wód powierzchniowych. Ten sposób prezentacji pozwala na bezpośrednie porównanie obrazów geochemicznych różnych środowisk. Kierując się wygodą użytkowania, mapy (opatrzone skalą liniową) wydrukowano w formacie nieco pomniejszonym (A3). Zabieg ten nie spowodował pominięcia żadnego szczegółu treści map. Dla zainteresowanych mapami w skali 1:25 000 istnieje możliwość udostępnienia całości opracowania lub poszczególnych map drukowanych na ploterze.