Experimenteel ontwerp en plotvoorbereiding
Het experiment werd uitgevoerd in juli en begin augustus 2023 en 2024. Gewoonlijk vindt het sojaboonenseizoen in Ohio ongeveer plaats van april tot november. De locaties die tijdens het groeiseizoen van 2023 werden gebruikt, bevonden zich in Clinton County (39.488, -83.769) en Union County (40.296, -83.283), Ohio. De locaties die tijdens het groeiseizoen van 2024 werden gebruikt, bevonden zich in Clinton County (39.489, -83.762), Preble County (39.893, -84.611) en Licking County (39.948, -82.501), Ohio.
Het experiment maakte gebruik van een willekeurig volledig blokontwerp om ruimtelijke effecten binnen het sojaboonveld te controleren. Binnen elk blok werden 5 percelen aangewezen voor een specifieke behandeling. De behandelingen bestonden uit twee groeistadia (R3 en R5) en twee aanbrengingstijden (midden in de ochtend en midden op de dag), en een onbehandelde controle. Dit resulteerde in de volgende 4 pesticidentoepassingen op percelen binnen elk blok: R3 midden in de ochtend, R3 midden op de dag, R5 op het midden van de ochtend, R5 op het middaguur. Toepassingen midden in de ochtend of midden op de dag vonden plaats binnen een venster van ~2 uur, met aanvragen rond 10:00 uur en middagapplicaties rond 12:30 uur. Naast de 5 behandelingspercelen binnen elk blok werden 3 percelen aangewezen als reservepercelen die in het geval van slechte vestiging of groei van sojabonen konden worden overgeplaatst naar de behandelingspercelen. De stadia van sojaboonplanten werden bepaald met behulp van stadiegidsen zoals die van de University of Minnesota Extension39.
Blokken sojabonenpercelen werden ongeveer 1,016 m (40 inch) uit elkaar geplaatst zodat onderzoekers gemakkelijk tussen hen konden lopen. De sojapercelen waren 1,9 m bij 8,5 m (0,0016 ha; 6,25 ft by 28 ft), en de sojaboonvariëteit die voor alle locaties werd gebruikt, was Pioneer P35T155E (Corteva, Indianapolis, IN). Op elk perceel werden 6 rijen geplant en stonden ze 38 cm uit elkaar. Er werd een zaaicapaciteit van 370.658 zaden per ha (150.000 zaden per acre) gebruikt. De bijen die bij het experiment betrokken waren, kwamen allemaal uit lokale populaties. Onderzoekers voegden geen aanvullende bijen toe. De auteurs hadden geen kennis van zowel stationaire als migrerende bijenboeren nabij de experimentlocaties. Alle materialen die voor het experiment zijn gebruikt, staan vermeld in de Materiaaltabel.
Pesticidebereiding
De gebruikte pesticiden waren Fitness-fungicide (41,8% propiconazol-werkzaam) en Fastac EC insecticide (10,9% alfa-cypermethrine werkzame stof). De toedieningssnelheden van pesticiden waren 6 fl oz/ac (438,47 mL/ha) voor het SBI-fungicide en 3,8 fl oz/ac (277,70 mL/ha) voor het pyrethroïde insecticide. Om de pesticiden voor toepassing voor te bereiden, werden 9,375 mL SBI-fungicide en 5,937 mL pyrethroïde insecticide gemengd met water om een eindvolume van 3 L als tankmengsel te bereiken. Chemicaliënresistente handschoenen werden gedragen bij het bereiden van het mengsel. De uiteindelijke concentratie SBI-fungicide in de tank was 3,125 mL per L. De uiteindelijke concentratie pyrethroïde insecticide in het aquarium was 1,979 mL per L.
Toepassing van pesticiden
Dit mengsel werd aangebracht op individuele sojaboonpercelen binnen elk blok met behulp van een op maat gemaakte handzame rugzakspuit. Het pesticidemengsel werd aangebracht bovenop het sojaboondak, met een snelheid van ~5 km per uur. Na het starten van het pesticide-aanbrengproces voor een bepaalde behandeling duurde het ~10 minuten voordat alle percelen de toepassing hadden ontvangen. Voor elke 3 L tank kan een oppervlakte van ~0,027 ha worden bedekt met het pesticide-mengsel. Het verbruikte spuitvolume per hectare is ongeveer 113 L. Dit betekent dat ~0,18 L van het tankmengsel per grafiek is aangebracht. Om te voorkomen dat pesticiden door wind naar andere percelen verdwijnen, werden toepassingen op winderige dagen vermeden. Volledige persoonlijke beschermingsmiddelen (shirt met lange mouwen en lange broek, chemicaliënresistente handschoenen, schoenen plus sokken en beschermende bril) werd gedragen bij het aanbrengen van het pesticide. Als er na afloop van een pesticide-toepassing een teveel pesticidemengsel was, werd dit volgens aanbevelingen van het staatsmilieubeheersingsagentschap verwijderd.
Opstelling/gebruik van de blokfluiter
De algemene opstelling en het gebruik van audiorecorders lijkt op een protocol beschreven door andere auteurs40. Handzame audiorecorders (Sony ICD-PX370) werden bevestigd aan plastic elektrische hekpennen met hakenband, en open-cel schuim werd als voorruit over de microfoon geplaatst. Een kegelvormige 3D-geprinte regenhoes werd ~3 cm boven de recorder aan de paal bevestigd om deze tegen het weer te beschermen.
Audiorecorders werden vervolgens ~24 uur voor de eerste pesticide toegediend in het midden van elk perceel geplaatst en minstens 24 uur na elke toepassing laten draaien. Twaalf recorders werden ingezet voor elke toepassingsdag, vier voor elke tijdstipgroep (mid-op-dag, midden-dag, controle) in zowel R3 als R5.
Audiogegevensverwerking
Alle audio-opnamen werden verwerkt via de machine learning-tool "buzzdetect"41 , die te vinden is op: https://github.com/OSU-Bee-Lab/buzzdetect. Dit hulpmiddel is getraind om het zoemend geluid van insecten te onderscheiden van andere geluiden die vaak voorkomen in de experimentele omgeving, zoals auto's, combines of vliegtuigen. Het programma is niet in staat het gezoem van Apis mellifera te onderscheiden van andere bijensoorten. BuzzDetect kan echter zoemgeluiden onderscheiden van andere niet-zoemende insectengeluiden, hoewel specifieke geluiden (zoals krekelgeluiden, meestal 's nachts) nog steeds moeilijk te onderscheiden zijn. Meer specifieke instructies voor het uitvoeren van buzzDetect zijn te vinden op: https://buzzdetect.readthedocs.io/en/latest/gui/. Na audio-analyse werd R42 gebruikt om het detectiepercentage binnen intervallen van 1 uur en 8 uur te berekenen.
Ontbladeringsgegevens
Bladerverwijderingsbeoordelingen werden uitgevoerd door bladeren van het bovenste sojaboonbladerdak visueel te beoordelen op 10 planten in elk perceel dat bij het experiment betrokken was, inclusief alle behandelingen en onbehandelde controlepercelen zoals aanbevolen door de University of Nebraska Extension43. Deze beoordelingen werden 24 uur na het gebruik van pesticiden uitgevoerd tijdens de R5-groeifase.
Opbrengstgegevens
Sojabonen werden geoogst en de opbrengst werd gemeten in bushels per acre met behulp van een maaidorser en een graanmeterweegsysteem. De opbrengstmetingen werden aangepast naar 13% vocht. De percelen hadden elk 6 rijen; alleen de 4 middelste rijen werden geoogst voor opbrengstmetingen.
Statistische analyse
Alle statistieken zijn uitgevoerd in R42. De gegevens werden geanalyseerd met ofwel lineaire modellen met gemengde effecten via het lme4-pakket44 of bèta-regressie via het glmmTMB-pakket45. De ruwe gegevens van BuzzDetect zijn beschikbaar in het aanvullende materiaal. Voor deze studie werden drie hoofdtypen gegevens geanalyseerd. Gemiddelde bijenactiviteitsgegevens over 1 uur en 8 uur na het toedienen van pesticiden, de opbrengst van elk sojaperceel en de bladloof van sojabonen waargenomen in de R5-groeifase.
De gemiddelde bijactiviteit werd verkregen door het aandeel van de positieve tijd voor bijenzoem in een venster van 8 uur of 1 uur te berekenen. Het venster van één uur is bedoeld om de acute effecten van pesticidetoepassing op bijenactiviteit te schatten, terwijl het venster van acht uur bedoeld is om de acute en residuele effecten te schatten. Resterende toxiciteitsgegevens zijn niet beschikbaar voor alfa-cypermethrine, maar de resterende toxiciteit voor andere pyrethroiden is gerapporteerd op ~8 uur na de eerste toepassing46. Omdat de berekening van bijenactiviteit een verhouding is, werden verschillen in activiteit geanalyseerd met bètaregressie uitgevoerd met het R-pakket "glmmTMB". De modellen onderzochten de gemiddelde detectiepercentages voor 1 uur of 8 uur na het toedienen van een pesticide en hoe deze varieerden bij het vergelijken van behandelde percelen met controlepercelen (de "gespoten" conditie), sojaboongroeistadium, het tijdstip van de dag waarop het pesticide werd aangebracht, of een groeistadium en tijdstip van de dag interactie. De vergelijking tussen behandelde en controlegrafieken wordt weergegeven met de variabele "gespoten" in zowel de 1-uurs- als de 8-uursmodelformules.
In het 1-uur venster model waren de vaste effecten de "gespoten" conditie, groeistadium, tijdstip van de dag en een groeifase en tijdstip van de interactie. Het model had zowel een jaar als een jaar per locatie interactieterm die als willekeurige effecten werden ingebracht. De volledige modelformule was:
(1)
De variabele detectiesnelheid weerspiegelt de gemiddelde bijendetecties die werden waargenomen ten opzichte van het totale waarnemingsvenster—in dit geval 1 uur. De variabele die wordt gespoten komt overeen met de vraag of een perceel een pesticide heeft ontvangen of dat het een onbehandeld controleperceel was. De variabele time_of_day komt overeen met de vraag of de detecties zijn vastgelegd na de pesticidentoediening midden in de ochtend of het middag gedurende de gewenste tijdsbestek. De variabele fase komt overeen met de groeifase waarin de detecties zijn vastgelegd.
In het 8-uursvenster model waren de vaste effecten de groeifase en de "gespoten" conditie. Het 8-uursmodel bevat nachturen voor de middagsproeitijd, wat zorgt voor een oneerlijke vergelijking tussen de mid-ochtend- en middagsproeitijden. Dit komt doordat de aanbrengingstijden rond het middaguur niet de volledige 8 uur van residuele toxiciteit testen doordat bijen 's avonds inactief zijn, een groot deel van deze 8-uur vensters. Daarom werden middagapplicaties uitgesloten van het model, evenals de termen "tijd van de dag" en "tijd van de dag:podium". De volledige formule voor dit model was:
(2)
De variabelen detectiesnelheid, gespoten en stadium hebben in deze formule dezelfde betekenis als in de formule die wordt gebruikt voor de 1-uur vensters. Het enige verschil is dat de variabele detectiesnelheid wordt berekend voor een venster van 8 uur. De opbrengstgegevens werden geanalyseerd met een lineair model met gemengde effecten. Het model keek naar bushels per acre per perceel en hoe de opbrengst per behandeling varieerde. Het vaste effect in het model was de behandeling, en locatie en jaar werden als willekeurige effecten gepast. De volledige modelformule wordt geschreven als:
(3)
Bladverlies van sojabonen werd gemeten als het gemiddelde percentage ontbrekende bladoppervlak per perceel. Een lineair model met gemengde effecten werd gebruikt om het effect van het tijdstip van de pesticide op het ontbladeringspercentage te vergelijken. Het vaste effect in het model was het tijdstip waarop de percelen werden bespoten. Locatie en jaar werden als willekeurige effecten ingebracht. De volledige modelformule is:
(4)
Alle protocollen die in deze studie zijn gebruikt, zijn samengevoegd in Aanvullend Bestand 1.