Økotoksikologisk vurdering af begroningshindrende biocider og biocidfrie bundmalinger Bilag 1. Model til beregning af eksponeringskoncentrationer (PEC)1. IndledningI dette bilag beskrives den eksponeringsmodel, der blev anvendt til beregning af eksponeringskoncentrationer PEC (Predicted Environmental Concentration) af DCOI og zinkpyrithion (ZPT) samt deres vigtigste nedbrydningsprodukter. Den opstillede model benytter modelleringsprincipper, der generelt anvendes ved bestemmelse af PEC. Talrige modeller til eksponeringsberegninger er beskrevet i litteraturen. Her kan bl.a. nævnes:
2. Opstilling af beregningsmodelGenerelt består eksponeringsvurderinger af følgende elementer:
Disse elementer blev også anvendt i nærværende eksponeringsberegninger. 2.1 ScenarierDer er opstillet 2 standardscenarier til beregning af eksponeringskoncentrationer (PEC):
Begge standardscenarier er således placeret i Roskilde Fjord. Hvert scenarie blev karakteriseret m.h.t.:
Netto vandudvekslingen mellem Roskilde Fjord og udløbet ved Isefjord blev antaget at svare til netto ferskvandstilførslen til fjorden, som er angivet til at være ca. 1,25 × 10-4 m3 × s-1 pr. m af fjordens længde (Harremoës & Malmgren-Hansen 1989). Da Roskilde Fjord er ca. 38 km lang (Harremoës & Malmgren-Hansen 1989), tilføres i alt ca. 410.400 m3 vand/dag. Med et overfladeareal på ca. 125 km2 svarer dette til en netto vandudveksling på 0,003 m3/m2/dag. Netto vandudvekslingen er således meget lav. DHI (1994) angiver således også, at vandstanden i Roskilde Fjord primært er bestemt af vindforhold og tidevandsvariationer. Vandstandsvariationerne er bestemmende for strømforholdene og dermed for fjordens vandskifte. De største vandstandsvariationer induceres af vinden, men tidevandet er bestemmende for de regelmæssige minimale variationer og dermed for hvilket mindste vandskifte, der forekommer på kort sigt. Ved Hundested i nord er tidevandets normale størrelse ca. 20 cm, hvilket er tilsvarende den daglige variation ved roligt vejr, som er observeret af brovagten ved Kronprins Frederiks Bro ved Frederikssund. I bunden af fjorden ved Roskilde er tidevandsvariationen nede på 6 7 cm (DHI 1994). På grund af vinden er almindeligt forekommende variationer over længere perioder dog langt større.
For begge scenarier er det antaget, at stofkoncentrationen i det vand, der transporteres ind i det betragtede vandområde, er ubetydelig.
Antal både i det betragtede vandområde.
Gennemsnitlig opholdstid for bådens tyngdepunkt i vandområdet:
I tabel B1.1 er der givet en oversigt over de parametre, hvormed de to scenarier er karakteriseret. Disse parametre er anvendt i basisberegningerne. Tabel B1.1 Se her Opstilling af standardscenarie. 2.2 Emission af antibegroningsmidlerHastigheden, hvormed antibegroningsmidlet frigives til vandmiljøet, udtrykkes som : U = [frigivet masse pr. areal pr. tidsenhed]. Måling af realistiske frigivelseshastigheder af antibegroningsmidler er behæftet med store problemer, da frigivelseshastigheden er afhængig af en række faktorer, som
Der foreligger et udkast til standard (ISO 1999), hvorefter frigivelseshastigheden af antibegroningsmidler kan beregnes. Frigivelseshastigheden bestemmes ud fra et estimat over malingens levetid, hvor det simpelt antages, at al begroningsmiddel frigives i løbet af levetiden. Der regnes med en højere frigivelsesrate i de første 14 dage, efter at malingen er påført båden. Efter 14 dage antages frigivelsesraten at være konstant. Standarden tager ikke højde for, at frigivelsen under sejlads sandsynligvis er højere, end når båden ligger i havn, og vil derfor sandsynligvis overestimere frigivelseshastigheden, når båden ligger i havn og under-estimere frigivelseshastigheden ved sejlads. Udkastet til standarden foreslår endvidere typiske filmtykkelser og levetider for forskellige malinger (ISO 1999). Tykkelsen af det afslidte filmlag på 6 måneder (svarende til en sejlsæson) kan ud fra standardforslagets foreslåede filmtykkelser og levetider, beregnes til 42 mm (opløselig matrix), 38 mm (uopløselig matrix), 45 mm (tinbaseret selvpolerende maling) og 50 mm (tinfri selvpolerende maling). Disse lagtykkelser er i god overensstemmelse med, at Hempel har estimeret det gennemsnitlige afslidte filmlag til 42 mm pr. sejlsæson for lystbåde i Danmark (Hempel 1999c). Hempel har baseret beregningerne på hvor meget bundmaling, der sælges til det danske marked pr. år samt antallet af lystbåde/motorbåde over 6 m (svarende til dem, der males) samt deres gennemsnitlige bundareal. Til simulering af den forøgede frigivelseshastighed ved sejlads antages det, at 60 mm af filmlaget slides af ved 6 måneders konstant sejlads. Til simulering af den lavere frigivelseshastighed, når båden ligger i havn, antages det, at 30 mm af filmlaget slides af, hvis båden konstant ligger i havn i 6 måneder. Disse antagelser svarer til, at bådene er ude at sejle i ca. 2 måneder af en sejlsæson og ligger i havn i de resterende 4 måneder, samt at frigivelseshastigheden er dobbelt så stor ved sejlads sammenlignet med, når båden ligger stille. Ud fra fortrolige oplysninger fra Hempel (Hempel 1999b) omkring indhold af antibegroningsmidler, tørstofindhold, massefylde kan den gennemsnitlige frigivelseshastighed for de to typer antibegroningsmidler beregnes. Resultatet af disse beregninger er angivet i tabel B1.2. Tabel B1.2 Beregnede gennemsnitlige frigivelseshastigheder.
Den samlede frigivelse af aktivstoffet til vandmasserne pr. tidsenhed udtrykkes i modellen ved:
hvor
2.3 PEC-ModelModellen deles op i følgende dele:
Omsætning i vandsøjlen Følgende forhold i vandsøjlen blev taget i betragtning:
Følgende forhold i sedimentet blev taget i betragtning:
En massebalance for vandsøjlen og sedimentet for både udgangsstoffet og de vigtigste nedbrydningsprodukter blev opstillet. Følgende tre PEC (Predicted Environmental Concentration) blev beregnet:
De tre koncentrationer blev sat lig med "steady-state" koncentrationen, dvs. den koncentration, som de beregnede koncentrationer nærmer sig med tiden, når der simuleres en kontinuer frigivelse af stoffet til vandmiljøet. For alle stoffer blev baggrundsniveauet antaget at være 0. 2.4 Data for aktivstoffer2.4.1 DCOI Aerob nedbrydning Figur B1.1 viser det forsimplede nedbrydningsmønster, som blev antaget for DCOI. DCOI blev i første omgang antaget at blive nedbrudt til N-(n-octyl) malonaminsyre, som under frigivelse af CO2 omdannes til N-(n-octyl) acetamid samt til N-(n-octyl) b hydroxypropionamid. Disse to nedbrydningsprodukter antages at blive omdannet til en lang række forskellige organiske forbindelser, som er samlet under "Andre nedbrydningsprodukter". Halveringstiden for denne pseudo-reaktion er antaget at være den samme for N-(n-octyl) acetamid og N-(n-octyl) b hydroxypropionamid. Disse nedbrydningsprodukter vil i et vist omfang blive mineraliseret under dannelse af CO2 (halveringstid for denne omdannelse er antaget ens for alle "andre nedbrydningsprodukter"). Halveringstiden for omdannelsen af DCOI til N-(n-octyl) malonaminsyre blev bestemt på baggrund af et forsøg, hvor fjernelsen af DCOI blev målt i havvand fra Jyllinge Lystbådehavn over en periode på 72 timer ved 12° C (Jacobson & Kramer 1999). Ved at minimere summen af arealerne af de relative rest af DCOI (RRSQ), der angives af Jacobson & Kramer (1999), ved anvendelse af følgende ligning: hvor i er en numerisk reference til observationen yi(observeret) er den målte nedbrydning (%) yi(estimeret) er den estimerede nedbrydning (%) ved antagelse af første ordens kinetik og ved antagelse af, at yi(estimeret) til tiden 0 = yi(observeret) til tiden 0 kan halveringstiden beregnes til 12,8 timer (for replikat 1) og 15,3 timer (for replikat 2) med en gennemsnitlig halveringstid på 14,1 timer (ved 12° C). De øvrige halveringstider blev estimeret ud fra de kvantiteter, som på baggrund af eksperimentelle undersøgelser (Mazza 1993) blev vurderet til at være til stede efter 30 dages aerob nedbrydning (se figur B1.1). De estimerede halveringstider for den aerobe omdannelse af DCOI ved 25° C er givet i tabel B1.3. Tabel B1.3 Estimerede halveringstider (dage) ved 25°C for aerob nedbrydning af DCOI.
Anaerob nedbrydning Den anaerobe nedbrydning af DCOI blev antaget at følge det samme reaktionsmønster som den aerobe nedbrydning. Nedbrydningshastighederne blev dog antaget at være langsommere for den anaerobe nedbrydning. Omdannelsen af DCOI ved aerobe og anaerobe forsøgsbetingelser er vist som funktion af tiden i figur B1.2. Data fra hovedrapportens tabel 3.2 (aerobe forhold) og 3.4 (anaerobe forhold; koncentrationer normeret til 100%) er anvendt. Det skal bemærkes, at tidsaksen for afbildningen af de anaerobe forsøg er 4,5 gange længere end tidsaksen for de aerobe forsøg. Det fremgår af figur B1.2, at der herved er en rimelig overensstemmelse mellem de målte koncentrationer for de aerobe hhv. de anaerobe forsøg. Det blev derfor antaget i beregningerne, at halveringstiderne for reaktionerne under anaerobe forhold er 4,5 gange længere end under aerobe forhold. Figur B1.1 Se her Nedbrydning af DCOI. Figur B1.2 Se her Nedbrydningsforløb af DCOI ved aerobe og anaerobe forhold. Stofegenskaber Udvalgte egenskaber for DCOI og dets nedbrydningsprodukter er givet i tabel B1.4. DCOI og de tre nedbrydningsprodukter blev antaget ikke at forefindes på ioniseret form ved pH = 7. Tabel B1.4
* Beregnet vha. KowWin (Syracuse Research Corporation 1996).* Calculated by means of KowWin (Syracuse Research Corporation 1996). ** Målt (data angivet i hovedrapporten), measured (data stated in the main report). 2.4.2 Zinkpyrithion (ZPT) Bionedbrydning Figur B1.3 angiver de forenklede biologiske nedbrydningsveje for ZPT, som blev simuleret i eksponeringsberegningerne. Det er en stærk forenkling i forhold til de meget komplicerede nedbrydningsveje for ZPT. Følgende forkortelser er benyttet:
Andre énringede, heterocykliske nedbrydningsprodukter karakteriseres med betegnelserne NP1-NP5 (jf. hovedrapporten, kapitel 4). Identiteten af NP1-NP5 er kendt af VKI. To hovednedbrydningsveje er antaget: 1) Primært under aerobe forhold:
2) Primært under anaerobe forhold:
Det blev antaget, at OMDS, NP3, PSoA og OMSo omdannes videre til andre forbindelser, der mineraliseres i mindre omfang. Halveringstiden for den primære reaktion (ZPT ® PT® OMDS + NP3) blev sat til 0,5 dage. De øvrige halveringstider blev estimeret ud fra de kvantiteter, som blev fundet i de aerobe og anaerobe nedbrydningsforsøg, hvor stofkoncentrationerne blev målt som funktion af tiden (disse forsøg er omtalt i hovedrapporten). De estimerede halveringstider for den aerobe og anaerobe nedbrydning er givet i tabel B1.5 og B1.6. Målte og beregnede koncentrationer er gengivet i figur B1.4. Tabel B1.5 Modelsimulering af aerob biologisk nedbrydning af zinkpyrithion. Estimerede
halveringstider (dage) ved 25°C.
Tabel B1.6
Fotolyse Som nævnt tidligere inkluderes fotolytisk nedbrydning af ZPT i modelsimuleringerne. ZPT antages at blive omdannet til NP3 ved den fotolytiske nedbrydning. Den fotolytiske halveringstid for ZPT i dagslys ved middagstid er blevet bestemt til:
Forsøgene blev foretaget i september ved en nordlig breddegrad på Den målte fotolytiske hastighedskonstant ved skyfri himmel er antaget til at kunne beskrives ved (se Schwarzenbach et al. 1993): hvor F er det såkaldte quantum yield, som her er antaget uafhængig af bølgelængden [-]. F angiver, hvor stor en del af de molekyler, som er exciteret af lyset, der omdannes til et andet stof.
l er bølgelængde [nm]. W(l) er lysintensiten ved de forskellige bølgelængder [milliEinstein/cm2/tid/nm]. Lysintensiteten ved breddegraden 40o og efterår er hentet fra Zepp & Cline (1977). D(l) er forholdet mellem den gennemsnitlige banelængde for lyset og den dybde af et vandelement, som kan antages fuldstændigt opblandet [-]. D(l ) er her antaget lig med 1 for alle bølgelængder. e (l) er den såkaldte molære ekstinktionskoefficient [(mol/l- 1/cm]. Værdier for disse er for ZPT taget fra (Fenn 1999). Der er taget hensyn til kuvettens lysafbøjning. F blev herved bestemt til 0,07. Det amerikanske program GCSOLAR (U.S. EPA 1999) blev herefter anvendt til beregning af den fotolytiske halveringstid. Programmet anvender de såkaldte attenuationkoefficienter, a1 , der benyttes til at beregne, hvor meget vandet absorberer lyset, som funktion af dybden. Attenuationkoefficienten for vandet ved Kronprins Frederiks bro, er bestemt ud fra målte sigtedybder ved 2 stationer i nærheden af broen (Roskilde Amt og Frederiksborg Amt, 1997). Den mindste sigtedybde i sommerperioden er på ca. 2,5 m. Ved anvendelse af data fra Calkins (1977) er følgende sammenhæng mellem sigtedybden og attenuation-koefficienten a ved sigtedybder under 4 m fundet: a [m-1] = 3,05 - 0,57× Sigtedybden [m]. Attenuationkoefficienten for Jyllinge Lystbådehavn er fundet ud fra litteraturdata i kystnære område (Zafirioiu 1977). Værdierne angives her som funktion af bølgelængde og med en minimum og maximum værdi. De højeste værdier er anvendt i nærværende beregninger. GCSOLAR tager ikke højde for skydækkets indflydelse på fotolysehas-tigheden. Det amerikanske program EXAMS (Burns et al. 1981), som også kan anvendes til simulering af fotolysen af stoffer, tager højde for skydækkets indflydelse på fotolysehastigheden. Ved brug af EXAMS er det beregnet, at halveringstiden i forhold til halveringstiden ved klar himmel bliver ca. 50% højere ved et skydække på 60%, som er det gennemsnitlige skydække i Danmark i månederne april-september (Danmarks Statistik 1996). GCSOLAR kan beregne den gennemsnitlige fotolytiske halveringstid for hver af årstiderne: Forår, sommer, efterår og vinter og ved forskellige breddegrader (dog kun breddegrader delelige med 10). Den gennemsnitlige fotolytiske halveringstid for ZPT i årstiderne forår, sommer, efterår samt breddegraderne 50o og 60o blev fundet til 9,8 timer for Jyllinge Lystbådehavn (6,5 timer, hvis der intet skydække er) og 6,6 timer for Kronprins Frederiks Bro (4,4 timer, hvis der intet skydække er). Den fotolytiske nedbrydning er svagt afhængig af temperaturen. Der er dog ikke inkluderet en temperaturafhængighed af den fotolytiske nedbrydning i nærværende beregninger. Den fotolytiske halveringstid er bestemt for et åbent vandområde, hvor der er taget hensyn til skydække og den faldende lysintensitet ned gennem vandsøjlen, men hvor der er ikke taget hensyn til, at de både, der er i havnen, samt molen vil afskære en vis del af sollyset. Der er derfor udført beregninger, dels hvor der er taget højde for den fotolytiske nedbrydning og dels hvor den fotolytiske nedbrydning er negligeret. De aktuelle forhold vil sandsynligvis ligge et sted mellem de to betragtninger, men det er ikke umiddelbart muligt at kvantificere, hvor stor betydning skyggepåvirkningen fra bådene og molen har på den mængde lys, der rent faktisk rammer vandoverfladen. For den befærdede sejlrute under Kronprins Frederiks bro vil der lige under broen være begrænset sol-lysadgang, medens der i den øvrige del af sejlruten ikke vil være større skyggevirkninger. Stofegenskaber Stofegenskaberne for zinkpyrithion og dets nedbrydningsprodukter er angivet i tabel B1.7. De forskellige énringede, heterocykliske nedbrydningsprodukter har en lav beregnet log KOW, hvilket indikerer en høj vandopløselighed. Sulfonsyrerne forventes også at være meget stærke syrer, hvorfor de sandsynligvis er fuldt dissocieret ved den herskende pH i de to vandområder. Bindingen til sediment og suspenderet materiale for disse to forbindelser forventes derfor at være lav. De beregnede KOC-værdier er dog anvendt til beregning af bindingen til suspenderet materiale og til sedimentet. Tabel B1.7 Egenskaber for zinkpyrithion og dets nedbrydningsprodukter.
* Beregnet vha. KowWin (Syracuse Research Corporation 1996).* Calculated by means of KowWin (Syracuse Research Corporation 1996). ** Målt (data angivet i hovedrapporten), measured (data stated in the main report). *** Fiktivt sat, fictive Figur B1.3 Se her Nedbrydningningsveje for zinkpyritihion (både aerobt og anaerobt). Figur B1.4a Se her Beregnede og målte stofkoncentrationer for zinkpyrithion og dets
nedbrydningsprodukter. Aerobe forsøgsbetingelser. Figur B1.4b Se her Beregnede og målte stofkoncentrationer for zinkpyrithion og dets nedbrydningsprodukter. Anaerobe forsøgsbetingelser. Calculated and measured concentrations of zinc pyrithione and its metabolites. Anaerobic experimental conditions. 3. BeregningsresultaterI tabel B1.8 er de beregnede PEC-værdier for de to scenarier og de forskellige stoffer angivet. Koncentrationerne er som tidligere nævnt "steady-state" koncentrationer. Det fremgår af tabel B1.8, at de højeste beregnede koncentrationer findes for lystbådehavnen. De beregnede koncentrationer er her ca. 100 gange højere end de koncentrationer, der beregnes for den befærdede sejlrute. For udgangsstofferne er følgende "steady-state" koncentrationer for vandfasen - PEC(vand) - beregnet:
Tabel B1.8a
Tabel B1.8b
4. FølsomhedsanalyseBeregningsresultaterne er bl.a. betinget af de værdier, som de forskellige parametre er tildelt. En følsomhedsanalyse af følgende parametres betydning for de beregnede koncentrationer af udgangsstofferne (DCOI og ZPT) er udført:
4.1 HavnescenarierDer er udført supplerende beregninger for 5 andre lystbådehavne. Karakteriseringen af disse er givet i tabel B1.9. Disse data er indhentet af Hempel og videregivet til VKI. Vandudvekslingen i havnene er sat til 0,6 m3/m2/dag for alle havne undtagen for Svendborg og for Horsens havn. For Horsens Lystbådehavn er vandudvekslingen sat til 0,8 m3/m2/dag. Svendborg Lystbådehavn står på pæle i Svendborg Sund. Strømforholdene er derfor antaget til at være tilsvarende strømforholdene i Svendborg Sund. En gennemsnitlig strømhastighed på 0,5 m/s er antaget for Svendborg Lystbådehavn, hvilket svarer til amplituden i den kraftige periodiske hastighedsvariation, som tidevandet bevirker i Svendborg Sund (Harremoës & Malmgren-Hansen 1989). Tabel B1.9 Se her De beregnede "steady-state" koncentrationer i disse havnescenarier er givet i tabel B10. Tabel B1.10
Det fremgår af tabel B1.10, at Jyllinge Lystbådehavn resulterer i de højeste beregnede koncentrationer. Den primære grund til dette er, at der er forholdsmæssigt flere både i Jyllinge Lystbådehavn i forhold til det vandvolumen, der skal fortynde det frigivne kemikalie (se tabel B1.10). Det anvendte scenarie er således et konservativt scenarie som forudsat ved udpegningen af Jyllinge Lystbådehavn. 4.2 Følsomhedsanalyse af øvrige parametreTabellerne B1.11a og B1.11b viser forholdet mellem den beregnede koncentration af udgangsstoffet i vandfasen efter ændring af parameter og den beregnede koncentration i standardscenariet. Det fremgår af tabellerne B1.11a og B1.11b, at inden for den variation, der er tillagt de enkelte parametre, er det den samlede udludningshastighed, som forårsager de største variationer i de beregnede koncentrationer. Sedimentationshastigheden har kun en meget lille betydning for de beregnede koncentrationer af udgangsstofferne. Tabel B1.11a Følsomhedsanalyse af PEC (vand) for DCOI. Tallene angiver forholdet mellem den beregnede koncentration af udgangsstoffet i vandfasen efter ændring af parameter og den beregnede koncentration i standardscenariet. Sensitivity analysis of PEC (water) of DCOI. The figures indicate the relation between the calculated concentration of the basic substance in the water phase, after changing the parameter, and the calculated concentration in the standard scenario.
Tabel B1.11b Følsomhedsanalyse af PEC (vand) for ZPT. Tallene angiver forholdet mellem den beregnede koncentration af udgangsstoffet i vandfasen efter ændring af parameter og den beregnede koncentration i standardscenariet. Sensitivity analysis of PEC (water) of ZPT. The figures indicate the relation between the calculated concentration of the basic substance in the water phase, after changing the parameter, and the calculated concentration in the standard scenario.
Udtrykkes lysintensiteten ved overfladen som I0 og lysintensiteten i dybden z som I, beskrives sammenhængen mellem I og I0 : log10(I/I0) = -a × z |
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||