Diseño experimental y preparación de la trama
El experimento se llevó a cabo en julio y principios de agosto de 2023 y 2024. Normalmente, la temporada de cultivo de la soja en Ohio tiene lugar aproximadamente de abril a noviembre. Los sitios utilizados durante la temporada de crecimiento de 2023 se situaron en el condado de Clinton (39.488, -83.769) y el condado de Union (40.296, -83.283), Ohio. Los sitios utilizados durante la temporada de cultivo de 2024 se encontraron en el condado de Clinton (39.489, -83.762), el condado de Preble (39.893, -84.611) y el condado de Licking (39.948, -82.501), Ohio.
El experimento empleó un diseño aleatorizado de bloques completos para controlar los efectos espaciales dentro del campo de soja. Dentro de cada bloque, se designaron 5 parcelas para un tratamiento específico. Los tratamientos incluyeron dos etapas de crecimiento (R3 y R5) y dos tiempos de aplicación (media mañana y mediodía), además de un control sin tratar. Esto resultó en las siguientes 4 aplicaciones de pesticidas en parcelas dentro de cada bloque: R3 a media mañana, R3 a mediodía, R5 a media mañana, R5 a mediodía. Las solicitudes durante la madrugada o mediodía se realizaban dentro de una ventana de ~2 horas, con las solicitudes a media mañana aproximadamente centradas alrededor de las 10:00 de la mañana y las de mediodía alrededor de las 12:30 PM. Además de las 5 parcelas de tratamiento dentro de cada bloque, se designaron 3 parcelas como repuestos que podían ser trasladadas a las parcelas de tratamiento en caso de un mal establecimiento o crecimiento de la soja. Las etapas de las plantas de soja se determinaron con la ayuda de guías de clasificación como la publicada por la Extensión39 de la Universidad de Minnesota.
Los bloques de parcelas de soja estaban separados aproximadamente 1,016 m (40 pulgadas) para permitir que los investigadores pudieran caminar fácilmente entre ellas. Las parcelas de soja medían 1,9 m por 8,5 m (0,0016 ha; 6,25 pies por 28 pies), y la variedad de soja utilizada en todos los sitios era Pioneer P35T155E (Corteva, Indianápolis, IN). Se plantaron 6 hileras en cada parcela y se separaron 38 cm (15 pulgadas). Se utilizó una tasa de siembra de 370.658 semillas por ha (150.000 semillas por acre). Las abejas involucradas en el experimento eran todas de poblaciones locales. No se añadieron abejas suplementarias por parte de los investigadores. Los autores no tenían conocimiento ni de colmenarios estacionarios ni migratorios cerca de los sitios experimentales. Todos los materiales utilizados para el experimento están listados en la Tabla de Materiales.
Preparación de pesticidas
Los pesticidas utilizados fueron el fungicida Fitness (principio activo de propiconazol al 41,8%) y el insecticida Fastac EC (principio activo alfa-cipermetrina al 10,9%). Las tasas de aplicación del pesticida fueron de 6 fl oz/ac (438,47 mL/ha) para el fungicida SBI y 3,8 fl oz/ac (277,70 mL/ha) para el insecticida piretroide. Para preparar los pesticidas para su aplicación, se mezclaron 9,375 mL de fungicida SBI y 5,937 mL de insecticida piretroide con agua para obtener un volumen final de 3 L como mezcla del tanque. Se usaban guantes resistentes a productos químicos durante la preparación de la mezcla. La concentración final de fungicida SBI en el acuario fue de 3,125 mL por L. La concentración final de insecticida piretroide en el acuario fue de 1,979 mL por L.
Aplicación de pesticidas
Esta mezcla se aplicó a parcelas individuales de soja dentro de cada bloque utilizando un pulverizador portátil personalizado para mochila. La mezcla de pesticidas se aplicó en la parte superior del dosel de la soja, a una velocidad de ~5 km por hora. Tras iniciar el proceso de aplicación del pesticida para un tratamiento determinado, se tardaron ~10 minutos en que todas las parcelas recibieran la aplicación. Por cada acuario de 3 L, se puede cubrir un área de ~0,027 ha con la mezcla de pesticidas. El volumen de pulverización consumido por hectárea es aproximadamente de 113 L. Esto significa que se aplicó ~0,18 L de la mezcla del tanque por parcela. Para evitar la deriva de pesticidas a otras parcelas debido al viento, se evitaron las aplicaciones en días ventosos. Se usaba equipo de protección individual completo (camisa de manga larga y pantalones largos, guantes resistentes a productos químicos, zapatos más calcetines y gafas protectoras) durante la aplicación del pesticida. Si había exceso de mezcla de pesticidas tras la conclusión de la aplicación, se eliminaba según las recomendaciones de la agencia estatal de control ambiental.
Configuración/uso de la grabadora
La configuración general y el uso de las grabadoras de audio es similar a un protocolo descrito por otrosautores 40. Las grabadoras de audio portátiles (Sony ICD-PX370) se fijaban a estacas de plástico de valla eléctrica con cinta de gancho y lazo, y se colocaba espuma de celda abierta sobre el micrófono como parabrisas. Se fijó una cubierta impermeable impresa en forma de cono en forma de 3D al estaca ~3 cm por encima del grabador para protegerla del clima.
Los grabadores de audio se colocaban en el centro de cada parcela ~24 horas antes de la primera aplicación de pesticidas y se dejaban funcionando al menos 24 horas después de cada aplicación. Se desplegaron doce grabadores para cada día de aplicación, cuatro para cada grupo de hora del día (media mañana, mediodía, control) tanto en R3 como en R5.
Procesamiento de datos de audio
Todas las grabaciones de audio se procesaban a través de la herramienta de aprendizaje automático "buzzdetect"41 , que se puede encontrar en: https://github.com/OSU-Bee-Lab/buzzdetect. Esta herramienta ha sido entrenada para distinguir los zumbidos de insectos de otros sonidos que ocurren con frecuencia en entornos experimentales, como coches, cosechadoras o aviones. El programa no es capaz de distinguir el zumbido de Apis mellifera de otras especies de abejas. Sin embargo, buzzDetect es capaz de distinguir ruidos de zumbido de otros insectos que no zumban, aunque ruidos específicos (como los cantos de grillos, principalmente durante la noche) siguen siendo difíciles de distinguir. Instrucciones más específicas sobre cómo ejecutar buzzDetect se pueden encontrar en: https://buzzdetect.readthedocs.io/en/latest/gui/. Tras el análisis de audio, se utilizó R42 para calcular la tasa de detección en intervalos de 1 hora y 8 horas.
Datos de defoliación
Las evaluaciones de defoliación se realizaron evaluando visualmente las hojas del dosel superior de la soja en 10 plantas de cada parcela involucrada en el experimento, incluyendo todos los tratamientos y parcelas de control no tratadas, según lo recomendado por la Extensión43 de la Universidad de Nebraska. Estas evaluaciones se realizaron 24 horas después de la aplicación del pesticida durante la fase de crecimiento R5.
Datos de rendimiento
Se cosechaban sojas y el rendimiento se medía en fanegas por acre usando una cosechadora de parcelas y un sistema de pesaje de grano. Las mediciones de rendimiento se ajustaron al 13% de humedad. Las parcelas tenían 6 filas cada una; solo se cosecharon las 4 filas centrales para medir el rendimiento.
Análisis estadístico
Todas las estadísticas se realizaron en R42. Los datos se analizaron con modelos lineales de efectos mixtos mediante el paquetelme4 44 o regresión beta mediante el paqueteglmmTMB 45. Los datos en bruto de BuzzDetect están disponibles en los materiales complementarios. Se analizaron tres tipos principales de datos para este estudio. Datos medios de actividad de abejas durante 1 y 8 horas tras la aplicación de un pesticida, el rendimiento de cada parcela de soja y la defoliación observada en la etapa de crecimiento R5.
La actividad media de las abejas se obtuvo calculando la proporción de tiempo positivo para el zumbido de abejas en una ventana de 8 o 1 hora. La ventana de 1 hora está destinada a estimar los efectos agudos de la aplicación de pesticidas sobre la actividad de las abejas, mientras que la ventana de 8 horas está destinada a estimar los efectos agudos y residuales. No se dispone de datos de toxicidad residual para alfa-cipermetrina, pero se ha reportado que la toxicidad residual para otros piretroides es de ~8 horas después de la aplicacióninicial 46. Dado que el cálculo de la actividad de las abejas es una proporción, las diferencias de actividad se analizaron mediante regresión beta realizada con el paquete R "glmmTMB". Los modelos examinaron tasas medias de detección para 1 u 8 horas tras la aplicación de un pesticida y cómo variaban al comparar parcelas tratadas con parcelas control (la condición "pulverizada"), la etapa de crecimiento de la soja, la hora del día en que se aplicó el pesticida o una etapa de crecimiento y una hora del día de interacción. La comparación entre los gráficos tratados y de control se referencia mediante la variable "pulverizada" tanto en las fórmulas de modelo de 1 hora como de 8 horas.
En el modelo de ventana de 1 hora, los efectos fijos fueron la condición de "pulverización", la etapa de crecimiento, la hora del día y un término de interacción en etapa de crecimiento y hora del día. El modelo tenía tanto un término de interacción por año como uno por sitio ajustados como efectos aleatorios. La fórmula completa del modelo era:
(1)
La tasa variable de detección refleja la media de detecciones de abejas observadas en relación con la ventana total de observación—en este caso, 1 h. La variable pulverizada corresponde a si una parcela recibió la aplicación de pesticidas o si era una parcela de control sin tratar. La variable time_of_day corresponde a si las detecciones se capturaron después de la aplicación del pesticida a media mañana o a mediodía durante el periodo deseado. La etapa variable corresponde a la etapa de crecimiento en la que se capturaron las detecciones.
En el modelo de ventana de 8 horas, los efectos fijos eran la etapa de crecimiento y la condición de "pulverización". El modelo de 8 horas incluye horas nocturnas para la pulverización al mediodía, lo que supone una comparación injusta entre los horarios de pulverización a media mañana y mediodía. Esto se debe a que los tiempos de aplicación a mediodía no analizan las 8 horas completas de toxicidad residual debido a que las abejas están inactivas por las tardes, una gran parte de estas ventanas de 8 horas. Por lo tanto, las aplicaciones de mediodía se excluyeron del modelo, al igual que los términos de "hora del día" y "hora del día:etapa". La fórmula completa para este modelo era:
(2)
Las variables tasa de detección, pulverización y etapa tienen los mismos significados en esta fórmula que en la fórmula utilizada para las ventanas de 1 hora. La única diferencia es que la tasa variable de detección se calcula para una ventana de 8 horas. Los datos de rendimiento se analizaron con un modelo lineal de efectos mixtos. El modelo analizó los bushels por acre para cada parcela y cómo variaba el rendimiento según el tratamiento. El efecto fijo en el modelo fue el tratamiento, y el sitio y el año se ajustaron como efectos aleatorios. La fórmula completa del modelo se escribe como:
(3)
La defoliación de hojas de soja se midió como el porcentaje medio de superficie foliar faltante por parcela. Se utilizó un modelo lineal de efectos mixtos para comparar el efecto del momento del día de aplicación del pesticida sobre el porcentaje de defoliación. El efecto fijo en el modelo era la hora del día en la que se pulverizaban los terrenos. El sitio y el año se ajustaron como efectos aleatorios. La fórmula completa del modelo es:
(4)
Todos los protocolos utilizados en este estudio se recopilan en el Archivo Suplementario 1.