Área de estudio

En la Figura 1 se muestra el área de estudio donde se instrumentaron cuatro sitios de estudio (98°56’-98°54’ O y 19°27’-19°28’ N), ubicados con dispositivo GPS Garmin ^® (Etrex 22x). El área se localiza en una transición entre un clima templado semiseco y subhúmedo, con una precipitación media anual de 700 mm (^{García, 1968}). Los suelos son profundos (>2.0 m), con pendientes de 1 a 2 %, de textura media y moderada capacidad de retención de agua (^{Jasso y Pimentel, 1985}). Para evaluar la partición de lluvia se trabajó con dos ejemplares de Prunus capuli (capulín) y de Casuarina equisetifolia (casuarina) mediante transectos, con acceso seguro a la red de energía eléctrica de corriente alterna a 120 V para alimentar los dataloggers.

Figura 1 Ubicación del área y los lotes de estudio. 

Para el registro de variables se diseñaron cuatro dataloggers (uno por sitio), a los que se les incorporó un sensor de humedad de suelo y tres balancines con interruptores magnéticos reed switch para estimar el flujo caulinar y de dosel, así como la precipitación directa (^{dos Santos et al., 2020}). Los componentes de cada datalogger incluyen una placa Arduino ^® Mega 2560, un chip de reloj en tiempo real (Módulo RTC DS3231), un módulo lector de memoria externa micro SD (para Arduino) y una pantalla OLED (SSD1306 128×64). Estos elementos se montaron en una placa de circuito impreso (PCB personalizada) diseñada a través de PCB Wizard y ensamblada. El Arduino se programó en el IDE de software Arduino de código abierto (C/C++) (^{Arduino, 2021}) a través de un cable USB para consultar las lecturas de los sensores cada 10 minutos y registrar los datos en una tarjeta SD. El suministro de energía fue a través de la red eléctrica.

Calibración de balancines en el laboratorio

El pluviómetro de balancín (diseño propio) se empleó para medir la precipitación en tiempo real, con tecnología Arduino, por medio de un interruptor magnético que se activa al inclinarse el balancín de PVC (Figura 2) cuando acumula 18 cm³ de agua (^{Strangeways, 2007}; ^{Santana et al., 2015}).

Figura 2 Diagrama del balancín construido en PVC. 

En el Laboratorio de Simulador de Lluvia del Colegio de Postgraduados (Campus Montecillo) (Figura 3) se calibró el volumen de vuelco registrado por el sensor magnético. En condiciones de campo se comparó la precisión de los balancines con mediciones de pluviómetros totalizadores que sirvieron como testigos para validar la precisión de los registros en condiciones reales (^{Krishnamurthi et al., 2015}).

Figura 3 Calibración de pluviómetros con el simulador de lluvia. 

La precipitación instantánea se midió al contabilizar las oscilaciones de los baldes del balancín a través de un imán y el interruptor magnético reed switch. El sensor, alojado en un tubo de vidrio, registró la tensión generada al cerrar un circuito eléctrico, información que fue transmitida al datalogger. En total, para los cuatro sitios, se utilizaron 12 balancines para medir el flujo caulinar, el flujo de dosel y la precipitación directa.

La calibración de los balancines se efectuó con base en los principios de calibración dinámica (^{Edwards et al., 1974}; ^{Lee et al., 1996}); es decir, se compararon los valores del simulador de lluvia con las mediciones del balancín para cuatro intensidades de lluvia: 0.33, 0.83, 1.08 y 1.62 mm min^-1, a partir de la metodología de ^{Morin et al. (1967)} y 10 repeticiones de cada una. Se consideró un volumen de volteo de 20 cm³ por balde y un área de captación de 660 cm².

Con los resultados obtenidos se elaboraron curvas de calibración y se comparó la tasa de intensidad de lluvia del simulador con la tasa registrada en el datalogger, a fin de conocer la precisión y confiabilidad de los datos.

La eficacia de la curva de calibración se evaluó a través de la Eficiencia de Nash-Sutcliffe (NSE, por sus siglas en inglés) y PBIAS, según se expone en las Ecuaciones 1 y 2, respectivamente (^{Somavilla et al., 2019}).

Donde:

n = Número de evaluaciones

Y ₁ ^ad = Caudal añadido

Y ₁ ^est = Caudal estimado por el instrumento

Y ₁ ^mean = Media del caudal añadido

Donde:

Y ₁ ^ad = Caudal añadido

Y ₁ ^est = Caudal estimado por el instrumento

Calibración de los sensores de humedad en el laboratorio

Para calibrar los sensores capacitivos de humedad del suelo se recolectaron muestras de suelo a 15 cm de profundidad con humedad conocida en los cuatro sitios de estudio. Las muestras se colocaron en frascos perforados con base en la metodología de ^{Gil-Marin et al. (2022)}. En el laboratorio se midió el peso inicial de las muestras (báscula Mettler ^® PC 2000) y se registraron lecturas con el sensor capacitivo (Módulo Higrómetro Capacitivo v1.2 para Arduino) (Figura 4), que funciona mediante la constante dieléctrica aparente del suelo.

Figura 4 Determinación de las lecturas con el sensor capacitivo. 

Los frascos se sellaron parcialmente para facilitar la entrada de agua por capilaridad y se sumergieron en bandejas con agua destilada a diferentes contenidos de humedad (10 a 80 mL). Después un día, al alcanzar la saturación, cada muestra se pesó y se midió la humedad gravimétrica con el sensor. El contenido gravimétrico del agua se determinó por la diferencia de peso posterior al secado, 10 veces con dos repeticiones. Se relacionaron las lecturas del sensor y las mediciones del contenido de humedad gravimétrica del suelo, con análisis de regresión lineal se determinó la precisión de los resultados a través del oeficiente de determinación (R ² ) (^{Gil-Marin et al., 2022}).

Humedad del suelo en condiciones de lluvia real

Una vez calibrados los sensores capacitivos de humedad de suelo, la evaluación de los niveles de humedad del suelo en los cuatro sitios de estudio se realizó con las salidas digitales registradas en el datalogger. El sensor se colocó a 10 cm de profundidad a fin de medir la humedad a nivel de las raíces de los árboles.

Flujo caulinar en condiciones de lluvia real

La medición del flujo caulinar (stemflow) se realizó mediante un sistema de recolección instalado en los troncos de los árboles, que incluyó un anillo colector a 1 m de altura, fijado al tronco para dirigir el agua hacia una manguera flexible, reforzada con soga de henequén y plastilina, que descargaba a un balancín previamente calibrado a unidades volumétricas (Figura 5). Para obtener la lámina de agua (mm) se utilizó la Ecuación 3, propuesta por ^{Barbosa et al. (2016)}.

Figura 5 Arreglo empleado para medir el flujo caulinar en Prunus capuli Cav. 

Donde:

F _c = Flujo caulinar (mm)

v _pb = Volumen promedio captado por los baldes del balancín por volteo (m³)

A _c = Área de la copa (m²)

n _b = Número de volteos de los baldes del balancín durante el evento de lluvia

Flujo de dosel en condiciones de lluvia real

Para la medición del agua que atraviesa el dosel (throughfall) se instalaron dos canaletas, de forma triangular (30 cm por lado y 130 cm de longitud), abajo del dosel de los árboles (Figura 6). El cálculo de la lámina de agua por evento de lluvia se realizó con la Ecuación 4 (^{Barbosa et al., 2016}).

Figura 6 Arreglo de las canaletas usadas para estimar el flujo de dosel. 

Donde:

F _d = Flujo del dosel (mm)

v _c = Volumen captado en la canaleta por evento de lluvia (m³)

A _cc = Área de captación de la canaleta (m²)

Para conocer la lámina instantánea de agua a lo largo del evento de lluvia, se instaló una segunda canaleta que conducía el flujo de agua hacia un balancín previamente calibrado.

En este caso, la estimación del flujo fue con el volumen captado por el balancín y el número de volteos mediante la Ecuación 5 (^{Barbosa et al., 2016}).

Donde:

F _db = Flujo del dosel con balancín (mm)

v _pb = Volumen promedio captado por los baldes del balancín por volteo (m³)

n _b = Número de volteos del balancín durante el evento de lluvia

A _cc = Área de captación de la canaleta (m²)

Precipitación directa en condiciones de lluvia real

Para la medición de la precipitación directa por sitio se utilizó un pluviómetro compuesto por dos baldes de plástico sobrepuestos de 29 cm de diámetro y 36 cm de altura; el diámetro inferior del embudo de lámina de acetato es de 3 mm. Al balde superior se le adaptó uno de los balancines calibrados y se perforó su fondo a fin de captar en el balde inferior el agua vertida por el balancín (Figura 7). Con la relación entre el volumen de agua colectado en el balde inferior y la superficie del pluviómetro se determinó la lámina de precipitación. La lámina de precipitación con balancín se obtuvo a partir del número de volteos registrados en el datalogger mediante la Ecuación 6 (^{Barbosa et al., 2016}).

Figura 7 Arreglo de los pluviómetros usados para medir la precipitación. 

Donde:

P _d2 = Precipitación directa con balancín (mm)

v _pb = Volumen captado por el balancín por volteo (m³)

n _b = Número de volteos del balancín durante el evento de lluvia

A _pl = Área de captación del pluviómetro (m²)

La cantidad de agua captada, en unidades volumétricas, se determinó con base en la Ecuación 7 (^{Barbosa et al., 2016}):

Donde:

P _d1 = Precipitación directa (mm)

v _pl = Volumen captado en el pluviómetro por evento de lluvia (m³)

A _pl = Área de captación del pluviómetro (m²)

Los registros de lluvia se efectuaron entre el 16 de julio de 2022 y el 18 de octubre de 2022, periodo que permitió recopilar datos de lluvia bajo una amplia gama de condiciones de viento e intensidad de preecipitación.

Calibración de balancines en el laboratorio (simulador de lluvia)

Los doce balancines empleados en el presente estudio se calibraron en un simulador de lluvia con 10 repeticiones por intensidad, para un total de 480 eventos. Se comparó la precipitación real con las lecturas de los balancines para calcular los coeficientes de determinación (^{Hudson, 1993}). Los resultados entre las mediciones de los balancines y los valores reales mostraron una alta correlación (R ² =0.9556 a R ² =0.9878) (Cuadro 1), y confirmó la eficacia de los balancines controlados por Arduino para monitorear la partición de la lluvia /en coberturas arbóreas.

Los resultados demuestran que los balancines controlados por Arduino ofrecen altas correlaciones en la medición de la lluvia a diferentes intensidades de precipitación (^{Fankhauser, 1997}; ^{Habib et al., 2008}; ^{Shedekar et al., 2016}; ^{Segovia-Cardozo et al., 2023}). Estas correlaciones consideran la influencia de factores como la estructura del pluviómetro y su configuración en las mediciones (^{Fankhauser, 1997}; ^{Shedekar et al., 2009}).

Se observó una variación significativa en los coeficientes de determinación entre los balancines, lo que indica que las condiciones específicas de construcción y montaje pueden afectar el rendimiento, lo cual se corrobora con los coeficientes de PBIAS (0.56 a 3.83). A pesar de estas variaciones, los balancines mostraron ser robustos y fiables (NSE=0.9763 a NSE=0.9907), aunque confirma su utilidad incluso frente a variaciones en condiciones de construcción y diferentes intensidades de precipitación (^{Chrit, 2022}).

Estos hallazgos subrayan el potencial de los balancines controlados por Arduino como herramientas valiosas para la investigación y monitoreo de la precipitación.

Calibración del sensor de humedad en el laboratorio

En la calibración de los sensores de humedad, para las muestras de suelo provenientes de los cuatro sitios de estudio, se relacionaron las diferentes lecturas del sensor con el porcentaje de la humedad gravimétrica de las muestras de suelo.

Las curvas de calibración se ajustaron a un modelo logarítmico con altos coeficientes de determinación: 0.9198, 0.9192, 0.9849 y 0.9264 de los sitios 1, 2, 3 y 4, respectivamente (Cuadro 2). Los valores obtenidos en la presente investigación son similares a los referidos por ^{Placidi et al. (2020)}, quienes realizaron un ajuste de su curva de calibración a una función de segundo orden y con un coeficiente de determinación de 0.945.

El análisis reveló que los sensores de humedad con Arduino son una herramienta confiable para medir con precisión la humedad del suelo, a pesar de las variaciones de los sitios de análisis. Los altos valores de correlación indican que las mediciones hechas con los sensores de humedad alcanzan una alta precisión y consistencia (^{Lloyd y Marques, 1988}; ^{Levia y Frost, 2006}).

Evaluación de los balancines en condiciones de campo

En la evaluación de los balancines, previamente calibrados en laboratorio, se continuó registrando una elevada correlación lineal positiva, bajo condiciones de campo, entre el flujo de dosel con balancín y el volumen de agua testigo. Los resultados reunidos en el Cuadro 3 confirman la precisión y eficacia de los balancines en la medición del flujo del dosel y lluvia directa de los cuatros sitios bajo estudio.

Las fuertes correlaciones en todos los sitios (Cuadro 3) sugieren que los balancines son una herramienta eficaz para medir el flujo de dosel en diferentes tipos de cobertura arbórea, lo cual es crucial para entender la partición de la precipitación en dichas especies (^{Staelens et al., 2006}). Asimismo, resultan ser instrumentos eficaces para medir precipitaciones y una alternativa en estudios meteorológicos, climáticos y aquellos que evalúen el impacto hidrológico de las características dasonómicas de los árboles (densidad, disposición, textura y forma de las hojas) en el escurrimiento e infiltración de la lluvia (^{Levia y Frost, 2006}).

Partición de lluvia directa (flujo caulinar y dosel) y humedad del suelo

El análisis de los datos recolectados de los balancines mostró una alta correlación positiva entre el flujo caulinar y la precipitación directa (Figura 8, ejemplo correspondiente a Casuarina equisetifolia del Sitio 1), en el cual las correlaciones fueron 0.8345, 0.9163, 0.8471 y 0.9087 para los sitios 1, 2, 3 y 4, respectivamente (Cuadro 4). Esto significa que a medida que aumenta la lluvia directa, también lo hace el flujo caulinar.

Figura 8 Relación entre precipitación directa y flujo caulinar del Sitio 1. 

En la Figura 8 se aprecia que el flujo caulinar representa alrededor de 3 % de las precipitaciones, patrón que se repite en los demás sitios, similar a lo registrado por ^{Gómez-Tagle et al. (2015)}, quienes determinaron una proporción de escorrentía cortical de 2.4 % para un bosque secundario de pino-encino en el noreste del estado de Michoacán.

A pesar de mostrar mayor variabilidad, la correlación entre la precipitación directa y el flujo del dosel fueron altos (Cuadro 4): Sitio 1 (0.7723), Sitio 2 (0.9155), Sitio 3 (0.9177) y Sitio 4 (0.8332), e indicaron el buen desempeño de los balancines en campo.

Las relaciones entre la lluvia directa y la humedad del suelo tuvieron coeficientes de determinación (Cuadro 4) similares (0.5561, 0.5416, 0.5377 y 0.5479 para los sitios 1, 2, 3 y 4, respectivamente), lo que implica un incremento de humedad relativamente bajo entre eventos de lluvia, a medida que son muy intensas. Los coeficientes de determinación fueron similares a los calculados por ^{Cheng et al. (2020)}, quienes obtuvieron un R ² =0.3518 para cultivos abandonados y un R ² =0.5015 para matorrales.

La comparación de costos de los dataloggers frente a los precios de instrumentos de medición análogos es de USD $1 396 (Cuadro 5).