Adskillelse af effekter af herbicider og kvælstof på vegetation og leddyr i hegn og græslandsvegetation

3 Eksperimentel adskillelse af herbicid- og N-påvirkning

3.1 Materiale og metoder
   3.1.1 Forsøgsetablering
   3.1.2 Artssammensætning af eksperimentel vegetation
   3.1.3 Udsåning og etablering af vegetation
   3.1.4 Dataindsamling
   3.1.5 Databehandling
3.2 Resultater
   3.2.1 Jordbund og blok-effekter
   3.2.2 Etablering af plantearter
   3.2.3 Klassisk analyse af vegetationsdata
   3.2.4 Multivariat behandling af vegetationsdata
   3.2.5 Leddyr
3.3 Diskussion af felteksperimentets resultater
   3.3.1 Jordbund og forsøgsetablering
   3.3.2 Effekter af herbicidbehandling og N-tilførsel på vegetation og leddyr
   3.3.3 Hypotesene

Den eksperimentelle adskillelse af herbicidbehandling og N-tilførsel afspejler til en vis grad et successionsstudie under forskellige påvirkninger af N og herbicid. I relation til hegnsundersøgelserne kunne forsøgets tre første år afspejle billedet af naturlig indvandring af plantearter i nyplantede hegn på en dybdepløjet jord med forskellige grad af landbrugspåvirkning. I eksperimentet øgede vi dog 'successionshastigheden' ved at udså en lang række arter for hurtigt at få en ensartet vegetation til statistiske sammenligninger. Formålet med undersøgelsen var imidlertid ikke at belyse ændringerne over tid, så der vil for vegetationens vedkommende blive fokuseret på resultaterne efter tredje år. For leddyrenes vedkommende præsenteres resultater for alle tre år (2001-2003), men dog med størst fokus på 2003.

3.1 Materiale og metoder

3.1.1 Forsøgsetablering

En eksperimental vegetation blev etableret på nærings- og kalkfattig tør sandjord ved Ebdrup på Djursland. Marken havde tidligere været anvendt til korndyrkning, men havde ligget brak et år. Jorden blev dybdepløjet (60 cm dybt), harvet og tromlet for at begrave eksisterende frøbank fra pløjelaget. 120 permanente prøvefelter á 7x7m med 1,4m mellem prøvefelterne blev udlagt med 9m bufferzone rundt om eksperimentet, der er omgivet af trævækst på alle fire sider. 10 replikater af 12 unikke behandlinger blev etableret i et randomiseret blok design (Mead 1990). Behandlingene var 0, 1, 5 og 25 % af markdosis af glyphosat og 0, 25 og 100 kg N/ha/år.

For at undgå begrænsning af andre næringsstoffer fik alle felter hvert år i maj tilført 53 kg Fosfor, 141 kg Kalium, 50 kg Svovl og 0,7 kg Kobber /ha/år.

Glyphosat (RoundupBio, glyphosat-isopropylaminsalt, 360 g a.i./l) blev første gang udbragt efter etablering af tæt vegetation 24. august 2001 med vindhastighed på ca. 2 m/s, 50 cm over jordoverflade. Sprøjtebommen var 3 m bred med 50 cm mellem dyserne (Lurmark Lo-drift LD 015 Green) med et tryk på 2.0 bar. Anbefalet markdosis er ifølge producentens vejledning 4-6 l/ha. Vi valgte at bruge 4 l/ha (1440 g a.i./ha) som udgangspunkt for 0, 1, 5 og 25 % markdosis som resulterede i hhv. 0, 4, 20 og 100 ml glyphosat pr. 15 l vand. Anden og tredje glyphosatdosering fandt sted i april 2002 og 2003.

3.1.2 Artssammensætning af eksperimentel vegetation

I alt blev 30 arter med forskellige funktionelle egenskaber udsået. Udvælgelsen af arter er foretaget så objektivt som muligt med fokus på, at det skulle blive en græsdomineret vegetation, hvor de vigtigste livsstrategier (funktionelle grupper af planter) var repræsenteret. Denne sammensætning blev suppleret med enkelte kendte herbicid- og/eller nitrogenfølsomme arter, arter i tilbagegang i agerlandet og en sjælden art. Endvidere har aktuel forekomst i læhegn/græsland og tilgængelighed af kommercielle frø spillet en rolle

Udfra en bruttoliste på 255 arter blev der indhentet information om arternes livsstrategier. Karplanternes livsstrategier kan forklares med tre funktionelle syndromer, som betegnes CSR-strategier (Konkurrence-, Nøjsomheds- og Pionersyndromet): C for konkurrence (competition) tolerance og S for stresstolerante og R for ruderat/pionertilpassede eller forstyrrelsestolerante arter (Grime 1988, 1997; Ejrnæs 2000). Der findes ”rene” C-strateger som fx alm. kvik, men de fleste arter er en kombination af to eller tre syndromer. CSR-strategierne gives en værdi fra 0-12 for alle tre syndromer, hvor den maksimale sumværdi er12. For at få repræsenteret forskellige strategier i den udsåede vegetation blev arterne udtrukket stratificeret tilfældigt i de forskellige grupper, hvor der først udvalgtes en græsart for hver strategi for at sikre græsdominans. Derefter udvalgtes øvrige urter i de forskellige strategier.

3.1.3 Udsåning og etablering af vegetation

Spiringsevnen hos de udvalgte arter blev testet eksperimentelt i petriskåle før udsåning. På trods af at marken stort set var plan og tilsyneladende meget homogen, opstod der en vis vand- og vinderosion efter dybdepløjning og såning.

3.1.4 Dataindsamling

3.1.4.1 Jordbundsdata
Der blev indsamlet jordprøver i oktober 2001 og i august 2003 fra hver af de 120 felter. Kornstørrelsesfordeling blev kun foretaget i 2001-prøverne, mens kemiske analyser blev foretaget på 01 og 03-prøver. Seks delprøver jord taget 0-10 cm under A0 horisonten blev puljet til én prøve. Jorden blev sigtet i laboratoriet gennem 2 mm sigte for at fjerne rødder, sten, mv. Prøverne for 2001 blev inddelt i flg. fraktioner: grovsand: 200-2000µm, fint sand 63-200µm, groft silt 20-63 µm, silt 2-20 µm og ler <2 µm (Sørensen & Bülow-Olsen 1984).

Prøverne for 01 og 03 blev analyseret for pHH2O, ledningsevne, glødetab samt total N og P. Frisk jord extraheredes for uorganisk N og P i 1 time i 50 ml 0,5 M K2SO4. Ekstraktet blev filtreret gennem Whatman GF-D-filter og nedfrosset til analyse for NH4+-N med indophenol metoden, NO3--N med cadmium-reduktionsmetoden og PO4-P med molybdæn-blå-metoden (Allen 1989). PO4-P ekstrationen var meget lav og er derfor udeladt. En delmænge af jorden tørredes ved 70 grader C og blev knust. 4 gram heraf blev brændt ved 550 grader C. En anden fraktion (200 mg tørvægt) blev opløst i svovlsyre i en time for analyse af total N og P (Kedrowski 1983) med Hitachi U-2000 spektophotometer. C-total måltes med LECO CNS-1000 med metoden efter Sørensen & Bülow-Olsen (1984).

Kornstørrelsesfordeling blev foretaget af Forskningscenter Foulum og de kemiske analyser af Økologisk Kemisk Laboratorium, Botanisk Institut, Københavns Universitet.

3.1.4.2 Vegetation
Vegetationen blev i 2001 registreret før og efter udbringning af glyphosat (hhv. juli og oktober). I 2002 og 2003 blev vegetationen registreret i august måned. Efter tre vækstsæsoner blev vegetationen endvidere høstet i 9 systematisk placerede Raunkjær cirkler (0,1m²) i alle 120 prøvefelter. Alle levende karplanter (såvel såede som indvandrede) blev sorteret på artsniveau, tørret ved 70 i min. 48 timer og vejet. Dødt plantemateriale (førne) i cirklerne blev indsamlet og behandlet tilsvarende. Mosser blev indsamlet og artsbestemt ved mikroskopi, mens mos-frekvens blev bestemt ved tilstedeværelse i de 9 cirkler efter at karplanterne var høstet.

3.1.4.3 Leddyr
Leddyr blev indsamlet en gang hvert år, hhv. 19. september 2001, 9.-10. juli 2002 og 8.-10. juli 2003. Til prøveindsamlingen blev brugt en insektstøvsuger (d-vac, Dietrick 1961). Hver prøve bestod af ti sug à 10 sekunder (prøveareal ialt 0,9 m²) tilfældigt fordelt i hvert 7 m 7 m prøvetagningsfelt. Efter indsamling blev dyrene frosset ned og senere bestemt til arts- eller gruppeniveau, jvf. Bilag C1 (2001) og C2 (2002 og 2003).

3.1.5 Databehandling

3.1.5.1 Vegetation
Statistiske overvejelser er beskrevet overordnet i afsnit 1.3.1. Der blev foretaget en variansanalyse af datasættet (antal plantearter, mosfrekvens, total biomasse og førne biomasse) med den såkaldte Lineære Mixed-Effect (LME, Pinheiro & Bates 2000) i S-PLUS, som er parallel til Proc-Mixed i SAS. Blok-replikater var tilfældig variabel mens glyphosat- og N-tilførsel inkluderedes som systematiske variable. Inden analyse blev data log-transformeret før videre behandling for at opfylde forudsætningerne med varianshomogenitet og normalitet, jf. Sokal & Rohlf (1995). Biomassen er opdelt i arter og deres CSR strategi-værdi er udregnet som et vægtet gennemsnit (af biomasse) af de tilstedeværende planter i de enkelte felter. Der er testet med et signifikansniveau på 5 %.

DCA analysen af den eksperimentelle vegetation blev foretaget med følgende nødvendige forholdsregler, som anbefalet af Økland (1990). I alt 20 sjældne arter blev udeladt (forekom i mindre end 5 prøvefelter), ét prøvefelt blev udeladt (< fem arter til stede). Efter første DCA kørsel skilte fire prøveflader sig klart ud ift andre felter med > 1 s.d. enhed. De fire outliers blev fjernet fra analysen som anbefalet af Økland (1990). Der indgår derfor i alt indgår 115 prøvefelter (af 120 i alt) og i alt 29 arter. Ordinationen blev foretaget i PC-ORD (McCune & Mefford 1999) med standardopsætning. Kendall Tau-korrelationer blev beregnet for diskrete variable og produkt-moment korrelationer for kontinuerte variable. DCA-plottene er afbildet med middelværdien for en behandlingsgruppe som centrum (radius = 1 S.E.) i en cirkel.

CCA-analyse blev foretaget på samme datamatrice med CANOCO version 4.02 (ter Braak & Smilauer 1997-99) efterfulgt af fortløbende udvælgelse og Monte Carlo permutationstest for at vise den relative betydning af de enkelte variable og hvorvidt de forklarende variable afviger fra permuterede tilfældige variable (Sokal & Rohlf 1995). Der blev anvendt 999 permutationer uden restriktioner. De forskellige variable blev standardiseret for at sikre ens vægtning af alle variable.

For at analysere responsen af enkeltarter blev der lavet en indikator-artsanalyse (INSPAN) (Dufrêne and Legendre 1997) i PC-ORD. INSPAN blev udført på tre undergrupper af data, hvor enten herbicidbehandling eller N-tilførsel blev holdt konstant pga signifikante interaktionseffekter af de to faktorer i variansanalysen.

3.1.5.2 Leddyr
Variansanalyser er kun udført på overordnede grupper af leddyr (summer), dvs. ikke på artsniveau. Inden de statistiske analyser blev udført, blev leddyrdata transformeret (kvadratrod x + 1) for at kunne opfylde forudsætninger om normalitet og varianshomogeneitet. Alle data fra de eksperimentelle prøvefelter blev analyseret med Proc MIXED (SAS Institute 1989). Det samlede datasæt for alle tre år blev analyseret med glyphosat, kvælstof og glyphosatkvælstof som systematiske variable og blok og år som tilfældige variable (se 1.3.1). Data fra 2003 blev desuden analyseret alene, med glyphosat og kvælstof som systematiske variable og blok som tilfældig variabel. For de fleste dyregrupper var interaktionen ikke signifikant og blev derfor i de tilfælde udeladt i den endelige test.

Det fulde datasæt, hvor tæger, jordlopper og snudebiller er bestemt til art, blev undersøgt vha. den multivariate programpakke PRIMER version 5.2.9 (se 1.3.1). Data for jordbunds- og plantevariable blev inddraget som forklaringsvariable. Data blev om nødvendigt transformeret inden analyserne (primært for at bringe alle variable til nogenlunde samme størrelsesorden). Den indbyrdes korrelation mellem de tilgængelige data for plantebiomasse og –diversitet samt jordbundsforhold blev undersøgt, og for stærkt korrelerede variables vedkommende blev en enkelt af disse valgt ud. PI en 120120 matrice blev korrelationen mellem artssammensætningen og de udvalgte miljøvariable (relateret til jordbund og planter) undersøgt vha. BIOENV. Ved hjælp af MDS blev fremstillet plots af artssammensætningen i relation til behandlingsvariablene (glyphosat og kvælstof). To-vejs ANOSIM blev anvendt til at teste, om artssammensætningen af leddyrene var signifikant forskellig mellem prøvefelter behandlet med forskellige doser glyphosat og kvælstof. SIMPER blev brugt til at undersøge de enkelte arters bidrag til forskelligheden i artssammensætningen mellem prøvetagningsfelter med forskellig behandling undersøgt. Alle tests er vurderet på et 5 % signifikansniveau.

3.2 Resultater

3.2.1 Jordbund og blok-effekter

En variansanalyse på jordbundskarakteristika viste signifikante effekter af kvælstoftilførslen på en række målte jordbundskarakteristika (ledningsevne, pH, NO3 og NH4 og uorganisk N). Det kunne dermed konstateres, at N-tilførslen har givet forskelle i tilgængeligt næringstofmængder. Endvidere viste der sig visse forskelle i kornstørrelsesfordelingen i forsøgsområdet (Figur 3.1).

Figur 3.1. Fordeling af fint sand (venstre) og ler (højre) i forsøgsarealet. Langs 'x- aksen' er de 12 behandlinger fordelt på de 10 tilfældigt placerede replica på 'y-aksen'

Figur 3.1. Fordeling af fint sand (venstre) og ler (højre) i forsøgsarealet. Langs 'x- aksen' er de 12 behandlinger fordelt på de 10 tilfældigt placerede replica på 'y-aksen'

3.2.2 Etablering af plantearter

I Bilag B findes fotos af den eksperimentelle vegetation. Tre udsåede arter (langbladet vortemælk Euphorbia esula, almindelig agermåne Agrimonia eupatoria og almindelig torskemund, Linaria vulgaris) spirede aldrig i petri-skålene, men blev etableret i felten. Fire måneder efter såning i 2001 var 23 ud af de 30 (77 %) udsåede arter etableret. Efter 3 vækstsæsoner var over 90 % af de udsåede arter etableret og alle fire grupper (C, S, R og blandet) af plantestrategier var vel repræsenteret (Tabel 3.1).

Tabel 3.1. De udsåede arters højde, etablering og total biomasse efter 3 vækstsæsoner. (C = Competitiv strategi, S = Stress tolerant strategi, R = Ruderat strategi) ifølge Grime (1988) og Ejrnæs and Bruun (2000). Højde iflg. Lid (1987). Frøspiring er procent spiring af 20 frø i petriskål efter 59 dage (n=2). Etablering er udtrykt som procent af prøvefelter med succes fire mdr. og 2,5 år efter udsåning. De udsåede arters relative dominans er udtrykt som % af total høstet biomasse.

Klik her for at se Tabel 3.1

Mindst 3 arter fra hver gruppe blev etableret i min. 40 % af felterne efter 4 mdr. En græsart fra alle fire funktionelle grupper etablerede sig, men kvik (Elytriga repens) blev dog dårligt etableret. Dette blev delvist kompenseret ved en kraftig etablering af stortoppet hvene (Agrostis gigantea med en CR strategi), der var iblandet frø af alm. hvene (A. capillaris). De tre dominerende græsser udgjorde 86 % af den totale høstede biomasse. Af de enårige blev almindelig dværgløvefod (Aphanes arvensis) aldrig etableret, mens enårig rapgræs (Poa annua) og kornblomst (Centaurea cyanus) blev etableret, men forsvandt stort set igen efter 2,5 år. Derimod har to andre enårige, salomons lysestage (Lepidium campestre) og mark-forglemmigej (Myosotis arvensis), formået at forblive omend kun som en sporadisk del af biomassen.

31 arter indvandrede til prøvefelterne, heraf 11 mosser, men ingen af de indvandrede arter formåede at dominere vegetationen (99,5 % af biomassen stammer fra de udsåede arter). De fem mest almindelige mosser var Ceratodon purpureus (104 prøvefelter), Bryum sp. (24 prøvefelter), Brachythecium rutabulum (17 prøvefelter), Ditrichum heteromallum (9 prøvefelter) og Polytrichyum piliferum (8 prøvefelter).

3.2.3 Klassisk analyse af vegetationsdata

Resultater af variansanalysen på vegetationsdata viste en signifikant interaktion af herbicidbehandling og N-tilførsel (Tabel 3.2). Interaktionen bestod i, at forskellene mellem N-niverauerne mindskedes med øget herbiciddosering, eller sagt på en anden måde, at herbicideffekten (skønt overordnet set ikke signifikant) steg med N-tilførslen (Figur 3.2). For førnemængden var der ligeledes signifikant interaktion mellem herbicid- og kvælstofeffekterne (Tabel 3.2). Der var en markant stigning i førneproduktionen med øget N- tilførsel ved lave herbicidtilførsler, og ligesom for den totale biomasse øgedes herbicideffekten ved stigende N-tilførsel (Figur 3.2).

Tabel 3.2. Resultater af variansanalyse (Generaliseret lineær model LME -Linear Mixed Effect model) for vegetation på de eksperimentelle forsøgsfelter som funktion af herbicid behandling x produktivitet (N-tilførsel). Fed angiver statistisk signifikans. Ved signifikant interaktion mellem herbcid- og kvælstofeffekt giver det ikke mening at se på de separate effekter.

  Varians-
komponent %
p-værdier
Forklarende variabel Enhed Transfor-
mationer
Blok Residual Herbicid-
tilførsel
N-tilførsel Herbicid ×
N-tilførsel
Total plantebiomasse g/m² log (X+100) <1 100 0,39 0,0025 0,0039
Førnebiomasse g/m² log (X+100) 10 90 0,0002 <0,0001 0,0001
Karplantediversitet no. log (X+100) 1 99 0,039 <0,0001 0,019
Frekvens af mosser - log (X+100) <1 100 0,074 <0,0001 0,0001
C- værdier (vægtet) - - 3 97 <0,0001 0,006 0,06
S- værdier (vægtet) - - 3 97 <0,0001 0,02 0,32
R- værdier (vægtet) - - 3 97 <0,0001 0,21 0,79

For karplantediversiteten bestod interaktionen mellem de to variable i, at artsantallet faldt med øget N-tilførsel, når der ingen herbicidtilførsel var, mens N-effekten faldt ved øget herbiciddosering (Figur 3.2). Sagt på en anden måde var der en negativ effekt af øget herbiciddosering i de ugødede felter, mens der ved 25 kg N/ha stort set ingen herbicideffekt var, og ved 100 kg N/ha var der en positiv effekt af øget herbiciddosering på antallet af karplantearter. Mosfrekvensen var signifikant påvirket af inteaktionen mellem N-tilførsel og herbiciddosering (Tabel 3.2, Figur 3.2), idet der i de ugødede felter var en negativ effekt af øget herbicidpåvirkning (især 25 % markdosis), mens der ved 25 kg N/ha ingen herbicideffekt var, og ved 100 kg N/ha var der en positiv effekt af herbiciddoseringen på mosfrekvensen.

Figur 3.2. Effekter af glyphosatbehandling (markdosis 1440 g a.i./ha) på totalbiomasse, førneproduktion, antal karplanter og mosfrekvens samt effekter på konkurrenceplanter (C-strrateger), nøjsomhedsplanter (S-strateger) og ruderatplanter (R-strateger) beregnet som vægtet biomassse ved tre niveauer af N-tilførsel (0, 25 og 100 kg N/ha/år).

Figur 3.2. Effekter af glyphosatbehandling (markdosis 1440 g a.i./ha) på totalbiomasse, førneproduktion, antal karplanter og mosfrekvens samt effekter på konkurrenceplanter (C-strrateger), nøjsomhedsplanter (S-strateger) og ruderatplanter (R-strateger) beregnet som vægtet biomassse ved tre niveauer af N-tilførsel (0, 25 og 100 kg N/ha/år).

For de tre hovedgrupper af plantestrategier blev der fundet en signifikant effekt af både herbicidbehandling og N-tilførsel på den fremherskende plantestrategi (Tabel 3.2). Herbicideffekten var for alle grupper uafhængig af N-tilførslen. Den biomassevægtede C- og R-værdi faldt med øget herbiciddosering, mens S-værdien steg (Figur 3.2). Det betyder, at øget pesticidanvendelse giver lavere C- og R-værdier; men højere S-værdi i plantesamfundet.

Glyphosat kunne alene reducerede artsantallet med maks. 23 %, og kvælstofgødskning reducerede artsantallet med maks. 29 % (Figur 3.2). De to faktorer indvirker således klart og nogenlunde lige meget rent kvantitativt på diversiteten. N alene virkede stærkt på plantebiomassen (fra ca. 500 til 830 g/m²), mens glyphosat ved højt N kunne reducere biomassen tilsvarende.

Af Figur 3.2 kan det endvidere ses, at ved 0, 1 og 25 % herbicid var der ingen klar effekt af kvælstofgødskning på de tre strategier, mens der ved 5 % glyphosat var en tendens til at ved øget N-tilførsel steg andelene af C og R, mens S faldt.Ved 0 Herbicid: R var næsten konstant ved øget N-tilførsel, mens C udkonkurrerede S ved højere N.Ved 25 % herbicid: R næsten konstant ved øget N-tilførsel, S udkonkurrerede C ved højere N.

S-strateger udkonkurrerede altså C- og til dels R-strateger ved høj herbicidbelastning ved såvel oligo- som eutrofe forhold. Her skal man dog være opmærksom på, at de enkelte strategiers biomasseproduktion primært udgjordes af fåresvingel, stortoppet hvene og rejnfan. En tilsvarede situation med få arter, der dominerer, er dog ikke ualmindeligt i agerlandet.

En analyse af gennemsnitsværdier af total biomasse og diversitet for de 12 forskellige kombinationsbehandlinger viste, at uden N-tilførsel var der maksimal biomasse og diversitet ved ingen eller 1 % herbicid (Figur 3.3). Endvidere var der ved 25 kg N-tillførsel maximal biomasse ved 5 % herbicid, mens artsantallet var relativt ens for de fire herbicid-niveauer. Ved 100 kg N-tilførsel var der højst artsantal ved 25 % herbicid, men biomasse maximum ved 0 herbicid.

Figur 3.3: Gennemsnitsværdier af total biomasse for de 12 forskellige kombinationsbehandlinger (markdosis 1440 g a.i./ha) med indikation af artsrigdom. Symbol størrelsen er proportional med artsantallet, cirkler = 0 N, trekanter = 25 Kg N ha-1,firkant 100 kg N ha-1

Figur 3.3: Gennemsnitsværdier af total biomasse for de 12 forskellige kombinationsbehandlinger (markdosis 1440 g a.i./ha) med indikation af artsrigdom. Symbol størrelsen er proportional med artsantallet, cirkler = 0 N, trekanter = 25 Kg N ha-1,firkant 100 kg N ha-1

Analysen af plantebiomasse, artsrigdom og mosfrekvens mod de forklarende variable viste, at biomassen korrelerede kraftigt og positivt med førne, C- og R-strategier samt ler, men korrelerer negativt med S- strategien (Tabel 3.3). Artsrigdommen og mosfrekvensen korrelerede ikke stærkt med nogle af de viste variable, men mosser havde dog en vis positiv korrelation med fint sand og negativ korrelation med førne.

Tabel 3.3. Korrelationskoefficienter mellem biomasse, artsdiversitet, mosfrekvens og forklarende variable. Grime indikator værdier er baseret på vægtet gennemsnit af prøvefelter udregnet på biomassen.

Forklarende variable Biomasse Artsrigdom total Frekvens af mosser
Herbicid tilførsel -0,30 -0,08 0,06
Nitrogen tilførsel 0,33 -0,23 -0,34
Biomasse (g m-2) 1 0,07 -0,08
Førne (g m-2) 0,45 -0,12 -0,39
Fint sand -0,08 -0,20 0,38
Ler 0,46 0,23 -0,10
Ledningsevne 0,15 -0,21 -0,20
pH 0,21 0,16 0,13
Total P (2001) 0,01 0,18 0,33
Total N (2001) 0,31 0,26 0,01
C-værdi (Grime) 0,60 0,18 -0,19
S-værdi (Grime) -0,59 -0,18 0,16
R-værdi (Grime) 0,52 0,17 -0,09

3.2.4 Multivariat behandling af vegetationsdata

3.2.4.1 DCA-analyse
En DCA analyse af vegetationsdata viste, at længden af akserne var hhv. 2,5; 1,4 og 1,2 for de tre vigtigste gradienter. Førsteaksen viste meget stærke negative korrelationer for C- og R-strateger, for biomasse og førne og i faldende grad en række jordbundskarakteristika. Førsteaksen var også kraftig positivt korreleret med pesticidbehandlingen (Tabel 3.4).

Tabel 3.4. Korrellationskoefficienter mellem de to første ordinationsakser og forklarende variable. Alle de viste variable er signifikant korrelleret til en af de to akser med p-værdier < 0.01. De forklarende variable er sorterede ift værdien i DCA 1 søjlen. Fed indikerer de variable med størst korrellation med DCA 1 og 2.

  DCA1 DCA2
Total biomasse -0,70 -0,03
Førne -0,62 -0,06
Konduktivitet_03 -0,57 0,15
Glødetab_03 -0,41 0,04
Silt -0,35 -0,07
Total N_01 -0,32 -0,14
Ler -0,30 0,05
Total N_03 -0,26 0,08
Nitrogen tilførsel -0,13 0,25
PH_03 0,01 -0,30
Fint sand 0,09 -0,26
NO3_01 0,16 0,34
Herbicid behandling 0,49 0,13

Dette afspejledes tydeligt i DCA plottene, når de fire glyphosatbehandlinger afbildedes (Figur 3.4a). Her fremgår det, at DCA 1 primært er en glyphosatrelateret gradient, men også at glyphosat i et vist omfang påvirker anden-aksen, hvor 25 % herbicid placerer sig relativt isoleret fra de øvrige herbicidbehandlinger.

Andenaksen korrelerede svagt positivt med nitrogentilførsel og nitrat i jorden og svagt negativt med pH i 2003 og andelen af fint sand. Dette afspejledes mere tydeligt i DCA-plottene (Figur 3.4b) hvor de tre N-tilførselsniveauer blev afbildet. Der var en tydelig N-gradient i materialet på begge akser, hvor især 0 kg N isolerede sig fra 25 og 100 kg N (interaktion mellem herbiciddosering og kvælstoftilførsel).

Figur 3.4. DCA analyser af vegetationsdata for eksperimentet. Hver ring har middelværdien af de 10 replika som centrum og 1 standardfejl som radius. A viser herbicidbehandlingen (markdosis 1440 g a.i./ha), B de forskellige N-tilførselsniveauer.

Figur 3.4. DCA analyser af vegetationsdata for eksperimentet. Hver ring har middelværdien af de 10 replika som centrum og 1 standardfejl som radius. A viser herbicidbehandlingen (markdosis 1440 g a.i./ha), B de forskellige N-tilførselsniveauer.

3.2.4.2 CCA-analyse
Fem ud af de 27 målte variable forklarede signifikante dele af variationen i artsssammensætningen (Tabel 3.5). Den stærkeste forklarende variabel var klart herbicidtilførslen, der alene forklarede 19,6 % af variationen. Til sammenligning forklarede N-tilførslen kun 2,5 % af variationen i artssammensætningen. Ledningsevnen forklarede ganske meget af variationen (7,3 %), men også fint sand og ler havde signifikant betydning (hhv. 3,4 og 2,1 %).

Tabel 3.5. Forklarende variable og ranket varians forklaret af de enkelte variable beregnet med fortløbende selektion af variable enkeltvis efterfulgt at Monte Carlo Permutations test. Variable der forklarer < 1 % er opført alfabetisk. N.S.= Ikke signifikant, * betyder p<0,05, *** betyder p<0,001.

Klik her for at se Tabel 3.5

3.2.4.3 Indikatorartsanalyse
Når vi under kvælstoffattige (0 kg N) og kvælstofrige (100 kg N) forhold sammenlignede effekten af herbicider på de enkelte plantearter i eksperimentet viste det sig, at visse arter reagerer nogenlunde ens på herbicidtilførsler og andre forskelligt afhængigt af N-niveauet (Tabel 3.6 og 3.7). Eksempelvis reagerede rejnfan (Tanacetum vulgare) signifikant negativt på øget herbicidbehandling både ved lav og højt N-niveau, mens fx fåresvingel (Festuca ovina) reagerede signifikant positivt i begge tilfælde. Rejnfan var altså særdeles gly-

phosatfølsom, mens fåresvingel var meget robust over for glyphosat.

Tabel 3.6. Indikatorartsanalyse i de forskellige herbicid-doseringsgrupper (markdosis 1440 g a.i./ha) under oligotrofe (0 kg N) forhold. Kun arter med Indikatorværdi (IV) >20 er medtaget. i indikerer arter der ikke er udsået. Fed angiver signifikans (p < 0,05).

  Relativ frekvens i herbicid-
dosis gruppe (%)
  C S R IV-værdi Maks.
behand-
lings-
gruppe
p-værdi 0% 1% 5% 25 %
Linaria vulgaris 6 0 6 21 0 0,672 43 24 32 1
Verbascum thapsus 2 2 8 23 0 0,914 38 17 2 44
Ceratodon purpureusi       28 0 0,535 28 26 27 19
Hieracium pilosella 2 8 2 31 0 0,809 35 10 31 25
Leucanthemum vulgare 8 2 2 39 0 0,552 39 10 21 30
Tanacetum vulgare 9 0 3 70 0 0,012 77 13 8 1
Agrostis capillaris 4 4 4 35 1 0,081 28 35 26 11
Agrostis giganatea 6 0 6 35 1 0,118 24 35 21 20
Hypochoeris radicata 4 4 4 48 1 0,088 26 48 16 10
Betula pubescensi       24 25 0,094 2 6 12 80
Hypericum perforatum 5 2 5 28 25 0,311 11 35 11 43
Lotus corniculatus 2 8 2 26 25 0,453 45 2 2 52
Festuca ovina 0 12 0 57 25 0,001 12 8 23 57

Tabel 3.7. Indikatorartsanalyse i de forskellige herbicid-doseringsgrupper (markdosis 1440 g a.i./ha) under eutrofe (100 kg N) forhold. Kun arter med Indikatorværdi (IV) >20 er medtaget. i indikerer arter der ikke er udsået. Fed angiver signifikans (p < 0,05).

  Relativ frekvens i herbicid-
dosis gruppe (%)
  C S R IV-værdi Maks.
behand-
lings-
gruppe
p-værdi 0 % 1 % 5 % 25 %
Lepidium campestre 2 2 8 28 0 0,309 46 23 25 6
Agrostis giganatea 6 0 6 32 0 0,19 32 32 26 11
Elytrigia repens 9 0 3 37 0 0,09 74 13 9 4
Agrostis capillaris 4 4 4 39 0 0,101 39 23 31 7
Tanacetum vulgare 9 0 3 58 0 0,013 58 19 23 1
Artemisia vulgaris 9 0 3 28 1 0,286 1 94 0 5
Myosotis arvensis 0 3 9 21 5 0,22 0 4 70 26
Linaria vulgaris 6 0 6 30 5 0,428 16 37 43 4
Leucanthemum vulgare 8 2 2 48 5 0,296 13 13 60 14
Chenopodium albumi 3 0 9 30 25 0,048 0 0 0 100
Rumex acetosellai 2 5 5 30 25 0,05 0 0 0 100
Oenothera biennis 2 2 8 39 25 0,082 2 24 25 49
Ceratodon purpureusi - - - 42 25 0,007 15 15 24 47
Verbascum thapsus 2 2 8 45 25 0,114 41 0 2 56
Festuca ovina 0 12 0 57 25 0,002 6 15 23 57
Hypochoeris radicata 4 4 4 62 25 0,006 7 4 20 69

Andre arter reagerede forskelligt under de to N-niveauer – omend billedet ikke var så statistisk klart. Alm. kongepen (Hypochoeris radiata), hvidmelet gåsefod (Chenopodium album) og mosset Ceratodon purpureus klarede sig signifikant bedre i konkurrencen med højt glyphosatniveau og høj N-tilførsel, mens dette ikke var tilfældet ved lavt N-niveau, hvor der var en ikke-signifikant tendens til det modsatte.

Uden herbicidbehandling kunne vi finde indikatorarter med klar respons på forskellig N- tilførsel (Tabel 3.8). Alm. kongepen (Hypochoeris radiata), alm. kællingetand (Lotus coniculatus) og mosset Ceratodon purpureus klarede sig signifikant bedre i konkurrencen uden N. Kvik (Elytrigia repens) og rejnfan (Tanacetum vulgare) reagerede derimod positivt på N- tilførsel. Prikbladet perikon (Hypericon perforatum) reagerede intermediært ved at have optimum ved 25 kg N.

Tabel 3.8. Indikatorartsanalyse i de forskellige N-tilførselsniveauer uden herbicidbehandling. Kun arter med Indikatorværdi (IV) >20 er medtaget. i indikerer arter der ikke er udsået. Fed angiver signifikans.

  Relativ frekvens i N-dosis
gruppe (%)
Species name C S R IV-
værdi
Maks.
behand-
lings-
gruppe
p-værdi 0
kgNha-1
25
kgN ha-1
100
kgN ha-1
Ceratodon purpureusi - - - 49 0 0,006 49 38 13
Hieracium pilosella 2 8 2 46 0 0,187 51 32 17
Hypochoeris radicata 4 4 4 70 0 0,01 70 26 4
Leucanthemum vulgare 8 2 2 52 0 0,286 52 42 6
Lotus corniculatus 2 8 2 40 0 0,022 100 0 0
Rumex acetosellai 2 5 5 30 0 0,087 100 0 0
Artemisia vulgaris 9 0 3 29 25 0,359 0 96 4
Brachythecium rutabulumi - - - 26 25 0,155 22 64 14
Festuca ovina 0 12 0 44 25 0,399 40 44 16
Hypericum perforatum 5 2 5 38 25 0,028 25 75 0
Agrostis capillaris 4 4 4 37 100 0,849 32 31 37
Agrostis giganatea 6 0 6 41 100 0,282 19 40 41
Elytrigia repens 9 0 3 49 100 0,009 1 2 97
Lepidium campestre 2 2 8 43 100 0,155 9 20 71
Linaria vulgaris 6 0 6 29 100 0,765 28 24 48
Oenothera biennis 2 2 8 21 100 0,696 26 22 52
Tanacetum vulgare 9 0 3 59 100 0,079 19 22 59

3.2.5 Leddyr

Blandt de grupper af leddyr, der er sorteret (dvs. herbivore gruppe, løbebiller og rovbiller), var tæger de talrigeste med i gennemsnit 67 dyr per prøve i 2003, hvoraf de fleste var blomstertæger (48 per prøve) efterfulgt af frøtæger (12 per prøve). Cikaderne var også talrige (gennemsnitligt 56 per prøve), rovbiller forekom med et gennemsnit på 7 per prøve, Clavicornia-gruppen (skimmelbiller, mariehøns og glimmerbøsser) havde et gennemsnit på 4 per prøve, hvoraf 3 var skimmelbiller, og de øvrige grupper (bier/hvepse, svirrefluer, sommerfugle, bladlopper, løbebiller, smeldere, bladbiller, jordlopper og snudebiller) forekom alle i et gennemsnitligt antal under 1 per prøve.

3.2.5.1 Udvikling over år
Antallet af leddyr steg fra 2001 til 2002 og fra 2002 til 2003 (hhv. 94, 129 og 136 i gennemsnit per prøve). Antallet af arter/grupper steg fra 2001 til 2002, men ikke fra 2002 til 2003 (hhv. 13, 21 og 20 arter/grupper i gennemsnit per prøve à 0,9 m²). Data for de forskellige gødnings- og pesticidbehandlinger for hvert af de tre år er vist i Figur 3.5. De dominerende leddyrgrupper varierede ikke meget mellem årene; dog var antallet af tæger meget lavere (16 per prøve) i 2001 end de to følgende år, mens der samme år var mange sommerfuglelarver, og i 2002 var der mange skimmelbiller (11 per prøve) og glimmerbøsser (4 per pr. prøve) sammenlignet med de to øvrige år.

Set over alle tre år var der signifikante pesticideffekter på bier og hvepse, svirrefluer, cikader, de fleste tægegrupper, jordlopper, små rovbiller, rovbillelarver, alle Clavicorniagrupper (skimmelbiller, glimmerbøsser og mariehøns), det samlede antal leddyr og antallet af arter (Tabel 3.9). Der var i alle tilfælde tale om et fald i antallet ved højeste pesticiddosering sammenlignet med de usprøjtede felter (se Figur 3.5). For nogle grupper var der tale om et begyndende fald i antalllerede ved 1 % af markdosis, mens der for andre først var en negativ effekt ved 25 % af markdosis. Effekten af kvælstof var signifikant for næsten de samme grupper som pesticideffekten, med undtagelse af bier/hvepse og svirrefluer, som ikke var påvirket af kvælstofniveauet, og sommerfugle og løbebillelarver, som var signifikant påvirket af kvælstof, men ikke af herbicidet (Tabel 3.9). For de fleste grupper steg det gennemsnitlige antal med stigende kvælstoftilførsel, men for Scutelleridae-/Pentatomidae-tægenymfer var effekten omvendt, og løbebillelarvernes antal steg ved 25 kg N/ha sammenlignet med 0, men faldt så ved 100 kg N/ha.

Tabel 3.9. P-værdier for effekt af pesticid- og N-behandling på leddyrfaunaen. Testen er udført på data fra årene 2001-2003. I de tilfælde, hvor interaktionleddet har en signifikant effekt, kan effekten af de enkelte variable ikke vurderes ud fra tabellen.

  Glyphosat Kvælstof Glyphosat×
kvælstof
Bier og hvepse 0.0011 0.2915 0.170
Svirrefluer 0.0192 0.9456 0.401
Sommerfugle 0.7693 0.0472 0.369
Sommerfuglelarver 0.1740 0.0935 0.569
Cikader 0.0001 <0.0001 0.564
Bladlopper 0.0926 0.8719 0.101
Scutelleridae-/Pentatomidaenymfer 0.2428 0.0219 0.216
Frøtæger <0.0001 0.9109 0.257
Nymfetæger <0.0001 0.0067 0.154
Blomstertæger* <0.0001 <0.0001 0.030
Tæger ialt <0.0001 0.0007 0.060
Bladbiller 0.6885 0.7062 0.563
Jordlopper 0.0263 <0.0001 0.478
Snudebiller 0.1570 0.5030 0.702
Løbebiller < 5 mm 0.7547 0.1886 0.636
Løbebiller > 5 mm 0.0590 0.7944 0.716
Løbebillelarver 0.0635 0.0410 0.833
Rovbiller < 5 mm* <0.0001 <0.0001 0.004
Rovbiller > 5 mm* 0.2255 0.1627 0.031
Rovbillelarver* <0.0001 <0.0001 0.048
Skimmelbiller 0.0001 0.0031 0.563
Glimmerbøsser <0.0001 0.4540 0.858
Clavicornia ialt <0.0001 0.0046 0.715
Smeldere 0.3681 0.0830 0.190
Dyr minus bladlus og trips <0.0001 <0.0001 0.353
Dyr ialt <0.0001 <0.0001 0.353
Arter/grupper <0.0001 <0.0001 0.323

* Testet med interaktionsleddet

Af Tabel 3.10 fremgår, at blokeffekten maksimalt udgjorde 10 % af den samlede tilfældige varians, mens forskelle i forekomsten af leddyr mellem årene gennemgående udgjorde en større del (op til 56 %). Blokeffekten varierede mellem årene (interaktion), idet der i 2001 og 2003 var en kraftigere blokeffekt end i 2002. Variationen mellem felter inden for blokken (felt(blok) i Tabel 3.10) var af cirka samme størrelse som variationen mellem blokke.

Tabel 3.10. Procentvis bidrag til den samlede ”tilfældige” variation for det samlede leddyrdatasæt (2001-2003). Data er testet uden interaktionsleddet glyphosat*NH4 som systematisk variabel.

  Blok År År*blok Felt(blok) Residual
Bier og hvepse 0,0 17,2 0,2 4,4 78,2
Svirrefluer 0,1 3,9 5,4 1,4 89,2
Sommerfugle 0,0 6,8 5,2 0,0 87,9
Sommerfuglelarver 0,0 34,4 11,0 0,3 54,3
Cikader 10,2 8,6 16,4 10,9 53,8
Bladlopper 0,1 35,4 5,0 0,0 59,5
Scutelleridae-/Pentatomidaenymfer 1,1 45,1 3,3 0,0 50,4
Frøtæger 0,4 49,1 3,5 0,0 47,0
Nymfetæger 1,0 40,5 8,6 1,7 48,3
Blomstertæger 1,0 48,2 2,2 2,5 46,1
Tæger ialt 1,3 55,7 2,5 3,5 36,9
Bladbiller 0,0 0,5 0,0 0,0 99,5
Jordlopper 6,4 12,8 6,3 3,5 71,0
Snudebiller 2,6 11,9 3,5 6,9 75,0
Løbebiller < 5 mm 0,0 11,6 0,0 0,0 88,4
Løbebiller > 5 mm 0,0 9,5 0,0 0,3 90,2
Løbebillelarver 0,8 6,8 12,7 0,1 79,6
Rovbiller < 5 mm 0,0 25,4 11,7 0,0 62,8
Rovbiller > 5 mm 0,0 0,0 5,8 3,9 90,3
Rovbillelarver 0,4 40,0 0,0 0,0 59,6
Skimmelbiller 1,7 30,9 4,7 0,0 62,8
Glimmerbøsser 1,1 30,7 5,2 0,0 63,1
Clavicornia ialt 1,4 37,4 4,3 0,0 56,8
Smeldere 0,0 3,7 0,0 0,0 96,3
Dyr minus bladlus og trips* 5,7 9,7 11,6 7,7 65,3
Dyr ialt** 5,7 9,7 11,6 7,7 65,3
Arter/grupper 4,8 35,5 4,8 4,5 50,5

* Data for alle år
** Data for 2001 kun for enkelte prøvefelter

3.2.5.2 Leddyr år 2003
Ved en særskilt analyse af data fra 2003 viste det sig (Tabel 3.11), at der for de fleste grupper ikke var en interaktiv effekt af pesticid- og gødningsbehandlingen. For frøtæger, blomstertæger, tæger ialt, store løbebiller, små rovbiller, skimmelbiller, antal dyr og antal arter var en signifikant effekt af glyphosatbehandlingen. Kvælstoftilføslen havde signifikante effekter på cikader, Scutelleridae-/Pentatomidaetægenymfer, nymfetæger, blomstertæger, små rovbiller, rovbillelarver, skimmelbiller, Clavicornia ialt, antal dyr og antal arter.

Tabel 3.11. P-værdier (Proc MIXED) for effekt af glyphosatbehandling og kvælstoftilførsel på leddyr indsamlet i 2003. Effekten af hver af de to variable kan kun vurderes i de tilfælde, hvor interaktionen ikke var signifikant. Tal i parentes indikerer problemer med testen (for få dyr, manglende konvergens el. lign.).

  Glyphosat Kvælstof Glyphosat×
kvælstof
Bier og hvepse (0,0133) 0,8202 0,1828
Svirrefluer 0,5929 0,8926 0,1046
Sommerfugle 0,9598 0,1811 0,9059
Sommerfuglelarver 0,5686 0,1889 0,2413
Cikader 0,3064 0,0002 0,0894
Bladlopper 0,5902 0,3378 0,2592
Scutelleridae-/Pentatomidaenymfer 0,3719 0,0006 0,4546
Frøtæger <0,0001 0,5368 0,5604
Nymfetæger 0,2138 0,0380 0,3001
Blomstertæger <0,0001 0,0314 0,5148
Tæger ialt <0,0001 0,1900 0,4584
Bladbiller 0,3960 0,3714 0,4298
Jordlopper 0,3365 0,0001 0,1780
Snudebiller 0,1382 0,8517 0,2875
Løbebiller < 5 mm 0,2215 0,1940 0,7298
Løbebiller > 5 mm (0,0272) (1,0000) (1,0000)
Løbebillelarver 0,3960 0,3714 0,4298
Rovbiller < 5 mm * 0,0067 <0,0001 0,0052
Rovbiller > 5 mm * 0,0891 0,1814 0,0022
Rovbillelarver 0,3050 <0,0001 0,4544
Skimmelbiller* 0,0079 0,0001 0,0374
Glimmerbøsser 0,7949 0,8850 0,4985
Clavicornia ialt 0,0503 0,0003 0,0333
Smeldere 0,3960 0,3714 0,4298
Dyr minus bladlus og trips <0,0001 0,0010 0,1490
Dyr ialt <0,0001 0,0010 0,1490
Arter/grupper * 0,0485 0,0252 0,0343

* Testet med interaktionsleddet glyphosat*kvælstof som systematisk variabel.

Af Figur 3.5 ses, at antal arter og antal dyr per prøve er reduceret ved den højeste glyphosatdosering. Felterne, der modtog 100 kg N/ha, havde flere dyr og arter i den usprøjtede situation end de andre kvælstofniveauer, men faldet i antal dyr og arter som funktion af stigende pesticiddosering var også størst i 100 kg N-felterne (Figur 3.5). Den maksimale herbicideffekt var 47 % reduktion i antallet af dyr og 26 % reduktion af artsantallet. Kvælstoftilførsel resulterede i en øgning af antal leddyr på maks. 55 % og maks. 21 % flere arter (cf. Figur 3.5). Den pesticideffekt, der ses for antal dyr ialt, går igen for de fleste øvrige grupper med signifikant pesticideffekt. For bier/hvepse var der dog ingen klar tendens, mens store løbebiller (for få dyr), og skimmelbiller forekom i mindre antal ved 1 og 5 % markdosis end i de usprøjtede plots, men i højest antal i plots behandlet med 25 % af markdosis (middelværdier hhv. 3,1; 2,3; 1,8 og 4,1). Antallet af de forskellige grupper steg generelt med stigende kvælstoftilførsel, med undtagelse af Scutelleridae-/Pentatomidaetægenymfer, som forekom i faldende antal med stigende kvælstoftilførsel.

Den tilfældige variation mellem prøverne udgjordes for op til 30 %s vedkommende af variation mellem blokke (Tabel 3.12), mens variationen mellem felter inden for blokke udgjorde op til 64 %. Begge variationsparametre varierede meget mellem de forskellige grupper af leddyr.

Figur 3.5. Antal leddyr (minus bladlus og trips) og antal leddyrarter per prøve i årene 2001-2003 som funktion af glyphosatdosis (markdosis 1440 g a.i./ha) ved de tre forskellige gødskningsniveauer. Punkterne repræsenterer rmiddelværdier af 10 prøver, de vertikale linjer er standardfejlen på middeltallet.

Figur 3.5. Antal leddyr (minus bladlus og trips) og antal leddyrarter per prøve i årene 2001-2003 som funktion af glyphosatdosis (markdosis 1440 g a.i./ha) ved de tre forskellige gødskningsniveauer. Punkterne repræsenterer rmiddelværdier af 10 prøver, de vertikale linjer er standardfejlen på middeltallet.

MDS-plottene (Figur 3.6) viser, at der er en sammenfaldende gradient i artssammensætning på den ene side og glyphosatbehandling eller plantebiomasse på den anden, mens det ikke er tilfældet for kvælstoftilførslen. Med andre ord er forskelle i artssammensætning bedre korreleret med pesticidbehandlingen og plantebiomasse end med kvælstoftilførslen.

Den mest forklarende enkeltvariabel blandt de undersøgte jordbunds- og plantevariable (miljøparametre) var total plantebiomasse (korrelation 0,247 i BIOENV), efterfulgt af pesticidbehandlingen (korrelation 0,202). Ingen andre miljøvariable havde en korrelation på over 0,2 med artssammensætningen af leddyr. Når det maksimale antal forklarende variable blev sat til tre, var artssammensætningen af leddyr bedst korreleret med pesticidbehandlingen, den totale plantebiomasse og førnebiomassen (korrelation 0,351). Da der var en vis sammenhæng mellem de to biomassemål, blev førnebiomassen efterfølgende udelukket fra analysen, hvorefter pesticidbehandlingen, den totale plantebiomasse og det i 2003 målte NO3-N forklarede artsfordelingen bedst (korrelation 0,336). Når maksimalt fem forklarende variable indgik i analysen, var korrelationen bedst med en kombination af pesticidbehandling, kvælstoftilførsel, antal plantearter, førnebiomasse og bloknummer (korrelation 0,396). En opdeling af plantebiomassen i græsser, øvrige urter (ikke-græsser), mosser og træer og efterfølgende analyse viste, at artssammensætningen af leddyr var bedst forklaret af græsbiomassem (korrelation 0,14 i BIOENV), efterfulgt af den samlede plantebiomasse (0,127) og biomassen af øvrige urter (0,077).

Figur 3.6. Multidimensionel scalings-plots med hhv. pesticidbehandling (A), gødskning (B) og total plantebiomasse (C) lagt over, dvs. afstanden mellem boblerne er proportional med forskellen i artssammensætning, og størrelsen af boblerne svarer til behandlingsintensiteten/plantebiomassen.

Figur 3.6. Multidimensionel scalings-plots med hhv. pesticidbehandling (A), gødskning (B) og total plantebiomasse (C) lagt over, dvs. afstanden mellem boblerne er proportional med forskellen i artssammensætning, og størrelsen af boblerne svarer til behandlingsintensiteten/plantebiomassen.

To-vejs ANOSIM (Tabel 3.13) gav en overordnet R på 0,188 for effekten af glyphosatbehandlingen (maks. 1,0) på artssammensætningen og et signifikansniveau på 0,001 (testet på tværs af kvælstoftilførselsniveauerne). Artssammensætningen i de usprøjtede felter samt felter behandlet med 1 og 5 % af markdosis var signifikant forskellig fra felterne behandlet med 25 % af markdosis. For kvælstoftilførslen var der en R på 0,177 og et overordnet signifikansniveau på 0,001 (testet på tværs af pesticidniveauerne). Der var signifikante forskelle i artssammensætningen mellem alle kvælstoftilførselsniveauer.

Ingen enkeltarter bidrog iøjnefaldende til forskellene i artsammensætning mellem felter behandlet med forskellige glyphosat- og kvælstofdoser (SIMPER). Det maksimale bidrag var 11 % af forskellen mellem felter sprøjtet med 25 og 5 % af markdosis forklaret af Stenodema-nymfer (blomstertæge, lever på græsser), baseret på en gennemsnitlig frekvens på hhv. 4,21 og 6,12 i de to pesticidniveauer. Denne type tæge var i alle, undtagen et, tilfælde den mest forklarende for artsammensætningsforskelle mellem såvel pesticid- som kvælstoftilførselsniveauer, efterfulgt af andre tæger, cikader og små rovbiller i varierende rækkefølge.

Tabel 3.12. Procentvis bidrag til den samlede ”tilfældige” variation for leddyr indsamlet i 2003.

  Blok Glyphosat×
kvælstof(blok)
Residual
Bier og hvepse 0,0 6,3 93,7
Svirrefluer 0,0 2,2 97,8
Sommerfugle 7,3 0,0 92,7
Sommerfuglelarver 0,0 3,6 96,4
Cikader 30,4 34,6 35,0
Bladlopper 8,8 0,0 91,2
Scutelleridae/Pentatomidaenymfer 12,2 0,0 87,8
Frøtæger 11,1 48,0 40,8
Nymfetæger 11,1 0,0 88,9
Blomstertæger 12,5 57,4 30,1
Tæger ialt 13,6 60,1 26,3
Bladbiller 0,0 0,8 99,2
Jordlopper 5,3 0,0 94,7
Snudebiller 6,3 0,0 93,7
Løbebiller < 5 mm 0,0 1,6 98,4
Løbebiller > 5 mm 0,0 0,4 99,6
Løbebillelarver 0,0 0,1 99,9
Rovbiller < 5 mm 16,6 22,2 61,2
Rovbiller > 5 mm 2,0 0,0 98,0
Rovbillelarver 3,3 22,3 74,5
Skimmelbiller 8,9 11,8 79,3
Glimmerbøsser 15,8 0,1 84,1
Clavicornia ialt 6,4 26,6 67,0
Smeldere 0,0 0,1 99,9
Dyr minus bladlus og trips 18,9 64,3 16,8
Dyr ialt 18,9 64,3 16,8
Arter/grupper 8,5 0,0 91,5

Tabel 3.13. ANOSIM-test af effekten af glyphosat (markdosis 1440 g a.i./ha) og kvælstof på artssammensætningen af leddyr.

Variabel Overordnet test Parvis sammenligning p-værdi
Glyphosat
% markdosis
R = 0,188
p = 0,001
0 - 1
0 - 5
0 - 25
1 - 5
1 - 25
5 - 25
0,36
0,23
0,001
0,047
0,001
0,001
Kvælstof
kg/ha/år
R = 0,177
p = 0,001
0 - 25
0 - 100
25 - 100
0,041
0,001
0,001

3.3 Diskussion af felteksperimentets resultater

3.3.1 Jordbund og forsøgsetablering

Eksperimentet har generelt været meget vellykket i forhold til etablering og udvikling af vegetationen. Enkelte arter havde svært ved at etablere sig fra frø (fx kvik, stor nælde, gråbynke og agertidsel som er udbredte og konkurrencestærke hegnsplanter). Såbedet var stærkt soleksponeret, og få dages udtørring af frø kan have haft afgørende betydning for nogle arters etablering. Flere af C-strategerne formerer sig vegetativt, så de er vanskelige at etablere fra frø. De udsåede arter dominerer vegetationen og græsser opnåede den ønskede dominans – mens indvandringen af arter, specielt mosser, gav yderligere dimensioner til eksperimentet mht. koloniseringsmuligheder for planter.

Selvom arealet var udvalgt som værende homogen mht. jordbund, viste der sig dog visse forskelle i jordbunden. Den tilførte kvælstof har derudover haft tilsigtede signifikante effekter efter 3 år på kvælstofniveauet og dertil knyttede sideeffekter på pH og ledningsevne.

3.3.2 Effekter af herbicidbehandling og N-tilførsel på vegetation og leddyr

Hovedformålet med eksperimentet var at belyse effekten af herbicidafdrift og kvælstoftilførsel. Eksperimentet har påvist klare interaktioner mellem N- og glyphosatpåvirkningen for plantebiomasse, førne, karplantediversitet, mosfrekvenser og dominansen af konkurrencestærke plantearter. For plantebiomasse og førne var der en meget markant reduktion ved 25 og 100 kg N i kombination med 25 % glyphosat, men denne effekt udeblev ved 0 kg N. For mosserne var frekvensen ved høje kvælstofniveauer klart højere ved 25 % glyphosat end ved lavere doser, hvorimod det omvendte var tilfældet ved lavere kvælstofdoser (Figur 3.2).

Glyphosat ændrer totalt konkurrenceforholdene og dermed artssammensætningen af planter. Det skyldes sandsynligvis, at fåresvingel er både nøjsom og tilsyneladende relativt glyphosat-tolerant, således at den på trods af 25 % markdosis opnår en rimelig biomasse og dominerer felterne (Tabel 3.6), hvorimod de herbicidfølsomme arter forsvinder. Glyphosat forhindrer en række arter i at etablere sig og mindsker artsdiversiteten ved 0 kg N-tilførsel. Ved 25 kg N er der stort set ingen effekt af glyphosatbehandlingen, mens der ved 100 kg N er en positiv effekt af stigende glyphosatdosering, således at artsantallet ved 25 % af markdosis er sammenligneligt med antallet ved 0 og 25 kg N (Figur 3.2). Dette skyldes formentlig, at glyphosaten virker så negativt på nogle arter (konkurrencestrategerne), at deres tilbagegang giver nicher til relativt mere glyphosattolerante arter (bl.a. R- og S-strateger). Vi skal dog bemærke, at det især er meget få arter fra hver strategi, der dominerer vegetationen, så billedet kunne være anderledes med andre arter, selvom få arter, der dominerer i hver strategi, ikke vil være usædvanligt i hegnsvegetation. Glyphosat kan alene reducere artsantallet med ca ¼, og en tilsvarende reduktion i artsantal ses af kvælstofgødskning (Figur 3.2). Begge faktorer indvirker således klart og nogenlunde lige meget rent kvantitativt på diversiteten. Glyphosatens isolerede betydning på plantebiomassen er meget klar ved højt N-nivau, men næsten fraværende ved 0 kg N.

En mulig forklaring på, at C-strategerne tilsyneladende påvirkes mest af glyphosatbehandlingen kunne være, at C-strategerne oftest er højere og derfor måske kan 'skygge' for glyphosatafdrift på R- og S-strategerne. Vi har ingen målinger af afsat mængde glyphosat på de forskellige planter, men det er ikke vort indtryk, at der har været væsentlig forskel i afsætningen af glyphosat på de enkelte arter. Udbringningen er sket om foråret (maj), hvor de fleste planter endnu er små. Den potentielle forskel i eksponering bliver først mere udpræget senere i vækstsæsonen. Måske er konkurrenceplanterne generelt mere følsomme end de øvrige, fx pga. deres højere vækstrate. En sådan forklaring stemmer overens med, at ved stigende glyphosatdosering forsvinder den positive effekt af N-tilførsel på den totale plantebiomasse (Figur 3.2). En anden mulig forklaring på den interaktive effekt af N-tilførsel og glyphosatdosering er, at selv beskedne glyphosatdoser hindrer en konkurrencerelateret uddøen (competitive exclusion), som ellers typisk er resultatet af tre år med øget N-tilførsel, fordi de dominerende nitrofile arter bortskygger små og lave arter. Glyphosatbehandlingen vil på den måde kunne åbne vegetationen for mindre konkurrencestærke arter selv ved 100 kg N/ha/år.

Med de multivariate analyseteknikker viste det sig, at glyphosat er den væsentligste variabel til ændringer i vegetationssammensætningen (Tabel 3.6). Jordbundskarakteristika og N-tilførsel viste sig således at være af sekundær betydning i forhold til de anvendte glyphosatkoncentrationer. Det skyldes, at glyphosat ved 0 kg kvælstof øger uddøen af herbicidfølsomme arter og ved 100 kg N virker glyphosat især ved at reducere biomasse og førne. Dette åbner mulighed for øget indvandring/etablering af nye (nitrofile) arter, som ikke findes ved 0 kg N. Glyphosat kan således virke både reducerende og fremmende afhængigt af N-niveauet.

Indikatorartsanalysen viste, at nogle arter reagerer meget forskelligt på herbicidbehandling afhængig af N-tilførslen. Dette kan have været medvirkende årsag til, at de forskellige funktionelle plantegrupper har reageret så markant, idet deres biomasse ofte er domineret af en eller meget få arter. Hvis vi eksempelvis havde udsået en anden, mere glyphosatfølsom nøjsomhedgræs end fåresvingel, havde billedet måske set anderledes ud. Hvis stor nælde, agertidsel, gråbynke eller kvik var slået an, kunne billedet også have været anderledes. Derfor skal dette resultat tolkes med varsomhed. Det ændrer dog ikke ved, at arter har forskellige reaktionsmønstre på glyphosat og næringsberigelse, og at en vegetation, der udsættes for herbicidafdrift, meget vel kan være domineret af en eller få arter.

For leddyrene var der signifikante effekter af glyphosatbehandlingen på en lang række grupper, på antal dyr og på diversiteten, uafhængigt af N-tilførslen (Tabel 3.9, Figur 3.5). Specielt 25 % glyphosat slog kraftigt igennem på leddyrene. Eftersom glyphosat formentlig ikke er giftigt for leddyrene, må der være tale om en indirekte effekt via effekten på planterne, førnen og (indirekte) mikroklimaet. Dette underbygges af, at artsammensætningen af leddyrene i høj grad er korreleret med plantebiomasse og førnebiomasse. Den overvejende negative effekt af herbicidbehandlingen på leddyrene stemmer overens med fx resultaterne i Reddersen m.fl. (1998).

Der findes kun ganske få studier, der tidligere har beskæftiget sig med effekter af lav- herbiciddosis på vegetation sammensat af arter, der findes i (halv-) naturen og småbiotoper i agerlandet (Bicheludvalget 1999). Langt den overvejende del af pesticidforskningen har været rettet mod arter på markfladen, såvel ukrudtsarter som skadedyr og svampe. Der er lavet talrige dosis-respons studier af de forskellige midlers effekter på veldefinerede målorganismer som en del af godkendelsesprocedurer for pesticider. Derimod er denne markante effekt på dyre- og plantesamfund ikke påvist før. Endvidere har der været lavet en række studier af indirekte effekter af herbicider på leddyr og fugle i agerlandet (se review i Bicheludvalget 1999), men kun et enkelt studie (Kleijn & Snoeijing 1997) har tidligere forsøgt at kvantificere effekten af et herbicid (Starane 200; fluoroxypyr) på agerlandets kantvegetation. Kleijn & Snoeijing (1997) fandt et fald i artsrigdom med øget fluoroxypyrdosering, mens biomasse produktionen steg for græssers vedkommende og faldt for andre urter under højproduktive forhold. De fandt endvidere, at 5 og 10 % fluoroxypyr (af markdosis) reducerede biomasseproduktionen af indvandrende urter og øgede ekstinktionsraten af de tilstedeværende arter. Forsøget anvendte årlig slåning og fjernelse af biomassen, hvilket også påvirkede konkurrencen signifikant, så derfor kan resultaterne ikke direkte sammenlignes med nærværende forsøg.

Eksperimentet i dette projekt har ved at tildele varierende doser af N og ensartede doser af P og K kunnet påvise en række forventede effekter af kvælstof på levende biomasse og på førne. Nitrogentilførsel alene (dvs. uden glyphosatbehandling) øgede produktiviteten af en sådan vegetation med ca. 65 %, når vegetationen som i dette tilfælde var N-begrænset. Konkurrencestrategerne (især repræsenteret ved stortoppet hvene) viste en stigende dominans med øget N-tilførsel, mens Nøjsomhedsplanterne (S-strateger især repræsenteret ved fåresvingel) viste en vigende tilstedeværelse. Ruderaterne viste ingen signifikante ændringer. Artsrigdommen for karplanter og frekvensen for mosser er højst ved 0 kg N/ha og stærkt faldende ved stigende N-tilførsel. Dette bekræfter således teorierne om tålegrænser for luftbårent kvælstof til N-begrænsede økosystemer (Bobbink m.fl. 1998).

Der er næppe tvivl om at effekterne af øget N-tilførsel på de undersøgte leddyrgrupper, som jo primært består af herbivore arter, er afledt af den øgede plante- og førnebiomasse og faldende plantediversitet ved øget N (Figur 3.2 og 3.5). De fleste arter har gavn af den øgede mængde tilgængelig føde, bedre mikroklima og mere stratificerede vegetationsstruktur, der opstår ved øget vækst, mens nogle få arter kan tænkes at være negativt påvirkede i kraft af, at plantediversiteten falder. I nærværende studie sås kun et konsekvent fald for en type tæger (Scutelleridae/Pentatomidae), som lever på græs, hvilket ikke umiddelbart kan forklares af ovenstående argumentation. Den overvejende positive effekt i dette eksperiment kunne meget vel hænge sammen med, at vegetationen endnu ikke var fuldt udviklet, og biomassetilvæksten derfor har været den mest betydende faktor for leddyrene. En sådan effekt stemmer overens med resultaterne i Esbjerg & Petersen (2002), som fandt en positiv korrelation mellem plantebiomasse og leddyrbiomasse for visse herbivore grupper. Man kan sagtens forestille sig, at den forarmning af floraen, der sker ved længere tids N-påvirkning, på lang sigt vil resultere i en tilsvarende forarmning af leddyrfaunaen. Andre studier (Koricheva m.fl. 2000, Siemann 1998 og Haddad m.fl. 2000) tyder imidlertid på, at kun meget stationære leddyr vil påvirkes af en lokal forarmning af floraen som følge af fx N-tilførsel. I et studie af kvik-domineret flora sammenlignet med en artsrig ”naturlig” flora fandt Lagerlöf & Wallin (1993), at der godt nok var flere leddyrarter i plots med ”naturlig” flora, men kvik viste sig også at have en stor betydning som skjul for mange leddyrgrupper, primært rovlevende insekter.

Størrelsen af afdriften af ammoniak fra markflader og fejlplaceret gødskning er meget dårligt dokumenteret, hvorimod effekten af gødskning er undersøgt meget grundigt i en lang række undersøgelser (se reviews af Waide m.fl. 1999, Bobbink m.fl. 1998). De fleste er enige, om at produktivitet påvirker diversitet, men de underliggende mekanismer og mønstre er man ikke enige om (Waide m.fl. 1999). En stor del af undersøgelserne opererer med NPK-gødskning, mens forskning i N-tilførsel alene hovedsageligt er lavet på N-følsom vegetation i kvælstofbegrænsede økosystemer (højmoser, heder, skove etc.) i relation til at fastsætte tålegrænser for N-tilførsel fra luften, (Bak m.fl. 1999). Kleijn & Snoeijing (1997) har i lighed med dette projekt undersøgt kombinationer af herbicid (fluoroxypyr, 0-50 % markdosis) med produktivitet (0, 25 og 50 % af markdosis NPK-tilførsel), hvor resultaterne som forventet gav meget signifikante udslag med øget biomasseproduktion og faldende artsdiversitet (specielt lave arter blev udkonkurreret) med øget produktivitet. Effekten her kan dog ikke umiddelbart sammenlignes med nærværende forsøg, da der er tale om et NPK-forsøg.

3.3.2.1 Funktionel fortolkning af resultaterne
Resultaterne af forsøget efter 2 ½ år kan opgøres på en række måder. Man kan opfatte glyphosatbehandlingen som en stressfaktor, der påvirker væksten af arterne forskelligt. Derved kommer resultaterne til at ligne eksperimenter med fx kobberstress, hvor en enkelt art (fx krybhvene) kan komme til at dominere en kobberforurenet grund (Strandberg m.fl. in prep., Strandsby & Strandberg 1998). Resultaterne af eksperimentet kan derved tolkes som klassisk økotoksikologisk stress-påvirkning (selvom vi mangler de højere doser for at få en egentlig dosis-responskurve).

Resultaterne af glyphosatbehandlingen kan derudover tolkes som en forstyrrelse af økosystemet, der også har forskellig virkning på de forskellige plantestrategier. Forstyrrelser (fx. græsning, brand, slåning, herbicider) har oftest været undersøgt som enten-eller begivenheder, og en del undersøgelser påpeger, at mange halv-naturlige økosystemer har behov for forstyrrelser for at kunne opretholde artsdiversiteten (Mackey & Currie 2000). En eng har behov for græsning og/eller høslet for at opretholde sin struktur og funktion og de tilknyttede og tilpassede arter. Mangel på forstyrrelser går især ud over rosetplanter (typisk nøjsomhedsplanter) og ruderater samt mosser i græsland (Jacquemyn m.fl. 2003). Nogle typer af forstyrrelser forbedrer med andre ord vilkårene for visse plantearter og andre forringer mulighederne (Jutila & Grace 2002). Hvis man opfatter glyphosatbehandlingen som en forstyrrelse, kan resultatet af 25 % markdosis forklares med forbedrede konkurrencevilkår for en relativ glyphosattolerant nøjsomhedsplante (fåresvingel,) idet konkurrenceplanterne ikke kan klare forstyrrelsen.

Resultaterne af den eksperimentelle undersøgelse viser, at indvandring (med mosser som eksempel) er størst ved lave N-niveauer pga. lav biomasse (og dækning) af vegetationen. Derimod har mange planter som udgangspunkt svært ved at etablere sig ved høj N-tilførsel pga. konkurrencen fra C-strategerne, som imidlertid viser sig mere følsomme over for glyphosat i eksperimentet. Derved får mosserne muligheden for at etablere sig ved moderate herbiciddoser. Ved 5 % glyphosat-dosis var effekten af kvælstofgødskning markant, således at andelen af C- og R-strateger stiger ved øget N-tilførsel , mens S-strategernes andel falder (Figur 3.2). Med andre ord: S-strateger udkonkurrerede altså C- og til dels R-strateger ved høj herbicidbelastning ved såvel lav som høj N-tilførsel.

Vi har i denne rapport valgt at fokusere på år 3 (det sidste projektår), hvilket under alle omstændigheder kun er et billede i en succession, som kræver yderligere opfølgning, fx efter 5 og 10 år, specielt hvis resultaterne af eksperimentet skal sammenholdes med hegnsundersøgelsen (se kap. 4). Endvidere kunne det være særdeles interessant at indlægge en yderligere variabel, fx slåning i eksperimentet idet én af mulighederne for at forvalte eutrofieret vegetation er slåning og fjernelse af biomasse. Dette kan gælde for vegetationen i enge såvel som for vejkants- og hegnsfodposevegetation.

3.3.3 Hypotesene

På baggrund af resultaterne af eksperimentet kan hypotesen om, at herbicidafdriften vil påvirke diversitet og dominans af en græslandsvegetation og de tilknyttede leddyr, bekræftes med den modifikation, at især 25 % af markdosis glyphosat gav de forventede effekter. I de fleste tilfælde var der ingen signifikant effekt af 1 % dosis, hvorimod 5 % fx gav signifikante ændringer i artsdominansen.

Hypotese 2 vedrørende reduktion i artsantallet som følge af ammoniakafdrift simuleret ved N-tilførsel blev delvist bekræftet af experimentet. N-tilførslen reducerede artsantallet af en græslandsvegetation ved at øge dominansen af flerårige, konkurrencestærke plantearter på bekostning af mangfoldighed af planter og leddyr. Som forventet var denne reaktion stærkest hos planterne, da leddyrene formentlig reagerede på planternes vækst og tilstedeværelse.

Eksperimentets hypotese H3, der forudsiger interaktioner mellem herbicid- og kvælstofpåvirkning, blev fuldt ud bekræftet af eksperimentet.

 



Version 1.0 Maj 2004, © Miljøstyrelsen.