text new page (beta)
Spanish (pdf)
Article in xml format
Article references
Automatic translation
Send this article by e-mail
Cited by SciELO 
Similars in
SciELO 
Permalink
Introducción
La radiación solar disponible en el sotobosque es uno de los principales factores que influyen en el establecimiento y desarrollo de vegetación bajo el dosel forestal (Lieffers, Messier, Stadt, Gendron y Comeau, 1999; Macdonald y Fenniak, 2007; Hu, Yan, Wu y Li, 2010). La radiación solar, tanto directa como difusa, que incide en un determinado lugar tiene influencia sobre los procesos fisiológicos, morfológicos y reproductivos de los organismos presentes bajo el dosel forestal, así como sobre los procesos químicos y físicos del ecosistema (Hardy et al., 2004; Promis, Schindler, Reif y Cruz, 2009; Sonohat, Balandier y Ruchaud, 2004). El balance entre la radiación solar incidente en el dosel y en el sotobosque, así como su distribución, dependerán de la transmisión e intercepción de la radiación solar por y a través de la cubierta vegetal (Webster, Rutter, Zahner y Jonas, 2016), la cual está determinada por la densidad, tamaño y distribución de las copas de los árboles (i.e. su estructura) (Canham y Burbank, 1994; Hardy et al., 2004; Jennings, Brown y Sheil, 1999). De esta forma, es posible establecer una relación respecto de la cantidad y heterogeneidad de la radiación solar como resultado de los patrones espaciales del dosel (Canham y Burbank, 1994; Denny y Nielsen, 2017). Con base en lo anterior, el modelo de Beer-Lambert-Bouguer ha sido modificado para ser usado en el dosel forestal, considerando que la cantidad de luz transmitida a través de este dependerá de la densidad de su estructura, y es descrito como un medio turbio con elementos distribuidos de manera aleatoria (Larsen y Kershaw, 1996; Monsi y Saeki, 2005).
En este estudio, se evaluó al dosel forestal por medio del uso de fotografías hemisféricas (FH), dado que la estructura del dosel afecta directamente la cantidad, calidad, distribución espacial y temporal de la luz dentro del bosque, además de que influye en la presencia o ausencia de repoblación, así como en el establecimiento y desarrollo de diversas especies de plantas a nivel del sotobosque (Macdonald y Fenniak, 2007). La estructura del dosel también modifica las condiciones ambientales y edáficas (e.g. microclima, disponibilidad de nutrientes, pH), así como las funciones ecofisiológicas de los organismos en el bosque (Jennings et al., 1999). El uso de FH se ha generalizado como una herramienta para caracterizar la estructura del dosel forestal, así como para estimar regímenes de radiación solar presente en bosques y selvas (Nobis y Hunziker, 2005; Seidel, Fleck, Leuschner y Hammett, 2011). El empleo de las FH se basa en el supuesto teórico de la interacción entre la radiación solar y el dosel (Canham y Burbank, 1994; Frazer, Canham y Lertzman, 1999; Schleppi y Paquette, 2017), de tal manera que es posible calcular el valor del coeficiente de extinción de luz (k), a partir de correlaciones entre el ángulo de inclinación de las hojas, la cobertura del dosel (Pekin y Macfarlane, 2009) y el índice de área foliar (Chianucci y Cutini, 2013; Schleppi y Paquette, 2017) a partir de FH. De igual forma, mediante FH es posible obtener un registro detallado de la distribución espacial, tamaño y forma de los claros del dosel en función del ángulo cenital (Gonsamo, D’odorico y Pellikka, 2013), permitiendo el cálculo del índice de transmisión de luz (GLI), el cual es relativamente similar al índice de luz en claros propuesto por Canham (1988) (Frazer et al., 1999). El procesamiento y análisis de FH digitales, a través de software especializado, comprende la transformación de la posición de cada pixel en coordenadas angulares, la discriminación entre pixeles sin obstrucción (cielo) y obstruidos (vegetación) y el cómputo de la distribución de brillo del cielo. Esta información, de manera combinada o individual, dará como resultado estimaciones de parámetros de la estructura del dosel, índices de transmisión de radiación, frecuencia y duración de haces de luz, así como proporción de luz directa y difusa (Frazer et al., 1999).
Desde la perspectiva de la silvicultura, con este estudio se busca estimar la transmisión de radiación solar en rodales de Abies religiosa (Kunth) Schltdl. et Cham., a partir de características dasométricas, de la estructura del dosel (apertura del dosel) y del rodal (área basal o densidad), Esto permitiría predecir la radiación incidente directa y difusa, a nivel del sotobosque, y proponer tratamientos silvícolas que ayuden a alcanzar diversos objetivos relacionados con los procesos de repoblación, supervivencia, desarrollo, sucesión forestal y diversidad (Hale, Edwards, Mason, Price y Peace, 2009); ya que la radiación solar transmitida a través del dosel, es un factor ambiental crucial en los procesos descritos.
Objetivos
Analizar los efectos del dosel forestal en la transmisión de radiación solar al sotobosque en rodales puros de Abies Religiosa y evaluar si la estructura del dosel, así como parámetros a nivel del rodal, permiten explicar la variación de la radiación solar bajo el dosel.
Materiales y métodos
Área de estudio
El presente estudio se realizó en el conjunto predial Reserva Forestal “El Manantial” S.C. ubicado en la región de la “Sierra Norte” del estado de Puebla, entre los 19°42'27.4" y 19°44'37.6" de latitud norte y los 98°00'43.9" y 97°58'56.5" de longitud oeste dentro de la subprovincia fisiográfica Lagos y Volcanes del Anáhuac, como parte de la provincia Eje Neovolcánico, con presencia de pendientes de entre 20% y 60% y altitudes desde los 2650 m hasta los 3178 m (Castaños y Castro, 2014). Predominan en el área los suelos de tipo Andosol (Organización de las Naciones Unidas para la Alimentación y la Agricultura - Organización de las Naciones Unidas para la Educación, la Ciencia y la Cultura [FAO-Unesco], 1998). Como resultado del relieve, altitud e influencia de los vientos alisios del Golfo de México, se identifican dos tipos de clima, C(w1) por debajo de los 3000 m de altitud y C(w2) por arriba de dicha cota altitudinal, difiriendo únicamente en la cantidad de precipitación anual (Castaños y Castro, 2014; García, 1964). La vegetación corresponde a masas mezcladas de pinos (Pinus patula, P. pseudostrobus, P. ayacahuite var. veitchii, P. rudis y P. teocote) y oyamel (Abies religiosa) en diferentes proporciones y masas puras de pino (P. patula) así como de oyamel (A. religiosa) (Castaños y Castro, 2014).
Sitios de muestreo
Se establecieron 16 unidades de muestreo circulares de 500 m2 (radio = 12.61 m) (t0.05 con 7 g.l, s2 = 43.02 y estimador de la media = 24.47) (Freese, 1976), que representan una intensidad de muestreo de 2.3% (Schreuder, Ernst y Ramírez-Maldonado, 2004) del área que cubren las masas puras de oyamel en el predio donde se obtuvieron los datos dasométricos del arbolado (Castaños y Castro, 2014). Las variables dasométricas medidas fueron diámetro normal (cm) y altura total (m) y se registró la especie a la que pertenecía cada individuo con diámetro normal ≥ 5.0 cm (Bonham, 2013). Se establecieron relaciones lineales entre parámetros como diámetro cuadrático (dc), área basal (G) y densidad (N) con el régimen de luz presente bajo el dosel (Hale et al., 2009; Comeau, 2001). Se utilizaron FH para estimar el porcentaje de claros, índice de área foliar y radiación trasmitida a través del dosel, la FH mide la fracción de huecos en múltiples ángulos a zenit (Seidel et al., 2011). Dentro de cada unidad de muestreo, con una distribución espacial al centro, norte, sur, este y oeste, se ubicaron cinco puntos (cuatro a una distancia del centro de 8.92 m) para realizar la captura de fotografías hemisféricas y estimar la radiación solar además de la fracción de cielo visible (Nackaerts, Coppin, Muys y Hermy, 2000), a fin de establecer su relación lineal con parámetros del rodal (Chazdon, Pearcy, Lee y Fetcher, 1996).
Obtención de imágenes hemisféricas
Cada una de las fotografías hemisféricas fue capturada con una cámara digital Canon® EOS SL1 APS-C sensor CMOS de 18 megapixeles, provista de una lente hemisférica “ojo de pescado” (Fisheye) Rokinon® 8 mm F3.5 HD y un nivel de burbuja Polaroid® de doble eje. El equipo se montó en un trípode para ubicar el lente a una altura de 1.30 m sobre el nivel del suelo, y facilitar su nivelación horizontal y orientación hacia el norte. Para capturar las imágenes se empleó el protocolo estándar para fotografías hemisféricas usadas en la estimación de la estructura del dosel (Zhang, Chen y Miller, 2005), el cual se basa en la determinación del grado de exposición a la luz, para lo cual se controla la apertura del lente y la velocidad de obturación (Leblanc, Chen, Fernandes, Deering y Conley, 2005). El valor de apertura del lente se determinó en f5.6, una velocidad de obturación de 1/125 segundos y a una sensibilidad ISO 400 (1/4 de segundo). La calidad de las imágenes capturadas fue de 5184 píxeles × 3456 píxeles con una resolución tanto vertical como horizontal de 138 puntos por centímetro (350 puntos por pulgada) en formato JPEG a color y RAW. Se obtuvo un total de 80 fotografías de todos los sitios de muestreo, esto es, cinco fotografías por sitio.
Procesamiento y análisis digital de imágenes
Métodos ópticos indirectos, como la FH, se sustentan en mediciones de la transmisión de radiación a través del dosel, para lo cual aplican la ley de Beer-Lambert-Bouguer (Ecuación 1). Esta ley considera la cantidad total de radiación interceptada por el dosel y depende de la radiación incidente, de la estructura del dosel y de las propiedades ópticas del mismo. Para estimar la transferencia de radiación solar a través del dosel forestal, en este estudio, se ha adoptado y modificado dicho modelo para describir la atenuación de la radiación solar a través de un dosel horizontal y homogéneo (Monsi y Saeki, 2005; Forrester et al., 2014).
Donde K es la radiación entrante o incidente arriba y bajo el dosel, f es el espesor del dosel, el cual dependerá del índice de área foliar y de variables como la fracción de apertura del dosel y el ángulo de incidencia de la radiación solar, mientras que e es el número de Euler (Hardy et al., 2004).
El procesamiento y análisis de las FH capturadas se realizó con el software Gap Light Analyzer© (GLA) v.2.0 conforme al procedimiento descrito en el manual del usuario (Frazer et al., 1999; Nobis y Hunziker, 2005). Para el análisis de cada fotografía, durante la etapa de registro de cada imagen, se suministró información particular del sitio donde se realizó la captura de la imagen como latitud, longitud, altitud, exposición y pendiente; además, se estableció el período de análisis a considerar (365 días). Durante la etapa de configuración, previo a la clasificación de cada imagen, se proporcionaron algunos parámetros de ajuste definidos previamente por el usuario (Tabla 1) (Frazer et al., 1999). El índice de nubosidad se fijó en 0.5, que es adecuado para describir la nubosidad anual en América del Norte (Frazer et al., 1999). El valor de la fracción espectral, es decir el porcentaje de radiación fotosintéticamente activa, corresponde aproximadamente a 45% de la radiación anual incidente en una superficie horizontal (Hardy et al., 2004; Lieffers et al., 1999; Papaioannou, Nikolidakis, Asimakopoulus y Retalis,1996). La fracción de radiación directa incidente, para la mayor parte de las regiones del norte del continente, es de aproximadamente 50% de la radiación global anual (Beaudet et al., 2011; Canham y Burbank, 1994; Frazer et al., 1999). El coeficiente de transmisión bajo cielo claro, para la mayor parte de las regiones libres de polvo se encuentra entre 0.6 µm y 0.7 µm (Frazer et al., 1999; Gates, 1980), por lo que en este trabajo se utilizó un valor medio de 0.65 µm (Beaudet et al., 2011; Scanga, 2014); si bien el valor de la constante solar se ha calculado en 1360.8 W·m-2 ± 0.5 W·m2 (Kopp y Lean, 2011), en este estudio se empleó el valor recomendado en el manual del usuario de GLA© (Frazer et al., 1999; Schleppi y Paquette, 2017).
Tabla 1 Variables suministradas por el usuario para el análisis de las fotografías hemisféricas empleando el software GLA©.
| Variable | Valor |
|---|---|
| Índice de nubosidad | 0.5 |
| Fracción del espectro electromagnético | 0.45 |
| Fracción de la radiación directa incidente | 0.5 |
| Coeficiente de transmisión bajo cielo claro (µm) | 0.65 |
| Valor de la constante solar (W.m2) | 1367 |
| Intensidad del sol en condiciones de cielo nublado | UOC |
UOC = Universal Overcast Sky, el cual considera que en condiciones de cielo nublado uniforme la radiación incidente es isotrópica (Schleppi y Paquette, 2017)
Por medio del software GLA©, se obtuvieron estimaciones de variables de estructura del dosel: fracción de apertura del dosel (FAD) e índice de área foliar (IAF) (Frazer, Fournier, Trofymow y Hall, 2001; Webster et al., 2016), además de estimadores de radiación solar en función de la radiación incidente por encima del dosel: radiación directa transmitida a través del dosel (TransDir) y luz difusa transmitida (TransDif). Adicionalmente, para describir la duración en minutos de la radiación solar que potencialmente podría alcanzar un punto en el sotobosque (Hu et al., 2010), se obtuvo: duración promedio acumulada de los haces de luz (DAS), duración media de un haz de luz (DPS) y frecuencia relativa de haces de luz (FRH) (Frazer et al., 1999). La distribución de frecuencia relativa de haces de luz, se realizó utilizando la clasificación propuesta por Smith y Berry (2013). Esta clasificación divide a los episodios de radiación solar en función de la duración e intensidad de las variaciones de radiación presentes tanto en el dosel forestal como en el sotobosque y se detalla a continuación: haz de luz o “sun-flecks” con una intensidad menor a 100% y duración menor a 8 minutos (SF); episodios con intensidades de 100% y duraciones mayores a 8 minutos como parche o mancha de luz “sun-patch” (SP), mayores a 60 minutos claro de luz “sun-gap” (SG) y aquellos mayores a 120 minutos como claros “clearing” (CL).
De igual forma, por medio del modelo básico de transmisión de luz a través de claros del software GLA©, se estimaron los porcentajes de radiación incidente directa, difusa y total transmitida a través de claros en el dosel hacia el sotobosque (Frazer et al., 1999), y se determinó el porcentaje de radiación transmitida, que en conjunto se denomina índice de transmisión de luz (GLI) (Frazer et al., 1999; Canham, 1988; Bartemucci, Messier y Canham, 2006). El GLI es relativamente similar al índice de luz en claros (Gap Light Index - GLI) (Ecuación 2) propuesto por Canham (1988), que especifica el porcentaje de radiación fotosintéticamente activa (RFA) incidente, transmitida a través de un claro del dosel hacia cualquier punto sobre el sotobosque durante un período determinado. Este índice permite relacionar, a pesar de la variación espacial y temporal de los niveles de luz incidente en el sotobosque, el tamaño y la forma de un claro con la transmisión RFA tanto directa como difusa a cualquier punto en particular dentro o fuera del claro.
Donde, P dif y P dir son las proporciones de RFA incidente en la parte superior del dosel como radiación difusa y radiación directa, respectivamente (P dif = 1 - P dir ); y T dif y T dir son las proporciones de radiación difusa y directa que son transmitidas a través de un claro en el dosel a un punto del sotobosque. Este índice toma valores desde 0, que es cuando no está definido un claro visible en el dosel, hasta un máximo de 100 que corresponde a un sitio en ausencia de dosel, o bien a cielo abierto (Canham, 1988; Bartemucci et al., 2006). Debido a que al incrementar la apertura del dosel (i.e. mayor presencia de claros), se incrementará el porcentaje de radiación tanto directa como difusa que será transmitido a través de este (Messier, Parent y Bergeron, 1998; Promis et al., 2009), se determinó la asociación entre los parámetros FAD y GLI mediante el coeficiente de correlación de Pearson.
Análisis estadístico
El análisis de las variables obtenidas se realizó con el paquete para computación y gráficos estadísticos R© (Venables y Smith, 2017), por medio del cual, para cada una de las variables de estructura del rodal y del dosel así como de radiación solar estimadas, se calcularon los valores medios, desviación estándar y coeficiente de variación (C.V.) por unidad muestral. Además, se realizaron análisis de correlación de Pearson entre las características estructurales y las variables de disponibilidad de luz bajo dosel, para después estimar asociaciones entre ellas por medio de análisis de regresión lineal simple (Ecuación 3).
Donde, β 1 es el intercepto teórico de la variable dependiente (Y) del modelo, β 2 es la pendiente teórica de la línea, X la variable independiente que afecta a Y, mientras que u es el error residual (Muenchen y Hilbe, 2010).
Resultados
El área basal (G) de las unidades de muestreo se distribuyó en un intervalo de 22.64 m2 a 53.55 m2 ha-1; mientras que el intervalo del diámetro cuadrático (dc) fue entre 14.0 cm y 35.0 cm; la densidad (N) presentó una distribución entre 460 y 1980 árboles por hectárea; la altura media entre los rodales fue entre 15.15 m y 23.97 m (Tabla 2). El análisis de correlación mostró la existencia de relaciones positivas significativas (p < 0.05) del diámetro cuadrático (dc) con los parámetros FAD (r = 0.53), DAS (r = 0.59), TranDir (r = 0.52) y GLI (r = 0.50) (Tabla 3). La relación de la densidad (N) fue también significativa (p < 0.05) pero negativa con DAS (R 2 = 29) y DPS (R2 = 0.27) (Fig. 1). El resto de los parámetros de estructura del rodal, régimen y distribución de luz, mostraron correlaciones nulas o débiles (p > 0.05), lo cual contrasta con otros trabajos que señalan relaciones significativas negativas no lineales entre la radiación bajo el dosel y el área basal del rodal (Buckley, 1999; Ligot et al., 2014); en particular, entre la fracción de radiación difusa transmitida (Comeau, 2001) y la capacidad de transmisión de luz a través del dosel con el área basal (Sonohat et al., 2004).
Tabla 2 Valores promedio por unidad muestral de los parámetros de estructura del rodal, estructura del dosel, transmisión de luz a través del dosel y régimen de luz en rodales de Abies religiosa.
| Parámetros/Unidad de muestreo | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Estructura del rodal | ||||||||||||||||
| Diámetro medio cuadrático (dg) cm | 21.8 | 22.0 | 30.7 | 26.6 | 19.6 | 34.9 | 29.4 | 28.4 | 21.6 | 34.2 | 30.2 | 34.2 | 24.5 | 14.0 | 15.2 | 18.0 |
| Altura media (Hm) m | 20.17 | 18.39 | 21.89 | 21.02 | 17.86 | 23.56 | 22.50 | 22.15 | 18.48 | 21.52 | 22.17 | 23.97 | 20.93 | 13.86 | 12.15 | 16.16 |
| Área basal (G) m2.ha-1 | 38.87 | 37.66 | 34.12 | 35.70 | 36.96 | 44.08 | 40.75 | 48.19 | 29.58 | 53.55 | 34.46 | 49.69 | 35.02 | 32.20 | 22.64 | 29.72 |
| Densidad (N) Ind.ha-1 | 1040 | 1000 | 460 | 640 | 1220 | 460 | 600 | 760 | 800 | 580 | 480 | 540 | 740 | 1980 | 1240 | 1160 |
| Estructura del dosel | ||||||||||||||||
| Fracción de apertura del dosel (FAD) | 0.23 | 0.22 | 0.24 | 0.36 | 0.24 | 0.34 | 0.31 | 0.37 | 0.44 | 0.38 | 0.41 | 0.35 | 0.31 | 0.29 | 0.23 | 0.19 |
| Índice de área foliar (IAF) m2.m-2 | 1.68 | 1.78 | 1.69 | 1.11 | 1.63 | 1.28 | 1.30 | 1.09 | 0.87 | 1.13 | 1.01 | 1.19 | 1.35 | 1.36 | 1.69 | 1.89 |
| Transmisión de luz | ||||||||||||||||
| Índice de transmisión de luz (GLI) % | 26.91 | 23.02 | 26.59 | 38.95 | 26.51 | 33.21 | 32.57 | 39.66 | 47.67 | 37.96 | 42.01 | 38.79 | 32.26 | 30.24 | 24.21 | 19.03 |
| Radiación directa transmitida a través del dosel (TransDir) | 5.04ac | 4.15ab | 4.78ab | 6.86ac | 4.78ab | 5.70ac | 5.80ac | 7.10bc | 8.58c | 6.77ac | 7.23bc | 7.38bc | 5.78ac | 5.15ac | 4.15ab | 3.18a |
| Radiación difusa transmitida a través del dosel (TransDif) | 4.47ab | 3.98ab | 4.62abc | 6.90bd | 4.59abc | 6.04ad | 5.71ad | 6.92bd | 8.27d | 6.65bd | 7.61cd | 6.33ad | 5.62ad | 5.54ad | 4.40ab | 3.55a |
| Dinámica de luz | ||||||||||||||||
| Duración media acumulada de luz (DAS) min.dia-1 | 178.92 | 170.01 | 183.70 | 264.10 | 174.16 | 233.50 | 231.55 | 282.36 | 326.54 | 271.08 | 279.75 | 295.34 | 233.81 | 194.25 | 157.35 | 114.07 |
| Número de "sunflecks" por día (NSD) | 33 | 36 | 32 | 26 | 28 | 39 | 33 | 33 | 32 | 43 | 23 | 30 | 38 | 40 | 37 | 19 |
| Duración promedio de haz de luz (DPS) min. | 5.23 | 4.71 | 5.82 | 10.33 | 6.10 | 5.86 | 7.07 | 8.55 | 10.32 | 6.35 | 12.85 | 10.09 | 6.10 | 4.79 | 4.22 | 6.04 |
| Proporción de haces de luz (SF) | 0.87 | 0.90 | 0.90 | 0.75 | 0.87 | 0.85 | 0.81 | 0.77 | 0.74 | 0.86 | 0.76 | 0.79 | 0.85 | 0.89 | 0.92 | 0.86 |
| Proporción de parches de luz (SP) | 0.12 | 0.10 | 0.10 | 0.22 | 0.12 | 0.14 | 0.18 | 0.22 | 0.23 | 0.14 | 0.20 | 0.18 | 0.14 | 0.11 | 0.08 | 0.14 |
| Proporción de claros de luz (SG) | 0.00 | 0.00 | 0.00 | 0.03 | 0.01 | 0.01 | 0.01 | 0.01 | 0.03 | 0.00 | 0.02 | 0.01 | 0.01 | 0.00 | 0.00 | 0.01 |
| Proporción de claros (CL) | 0.00 | 0.00 | 0.01 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.01 | 0.00 | 0.03 | 0.01 | 0.00 | 0.00 | 0.00 | 0.00 |
Unidades muestrales 1) Chichicaxtla Bis I, 2) Chichicaxtla Bis II, 3) Chichicaxtla I, 4) El Ocote I, 5) El Ocote II, 6) El Ocote III, 7) Los Corrales I, 8) Los Corrales II, 9) Los Corrales III, 10) Los Corrales IV, 11) Los Corrales V, 12) Los Corrales VI, 13) Los Corrales VII, 14) Ocojala Sierra Mojada I, 15) Ocojala Sierra Mojada II, 16) Ocojala Sierra Mojada III. Valores con la misma letra para radiación transmitida, ya sea directa o difusa, son iguales de acuerdo a la prueba de Tukey con p<0.05, unidades en Mols.m2.día-1.
Tabla 3 Coeficientes de correlación de Pearson y nivel de significancia (p-value) para diversos indicadores de estructura del rodal, dosel y distribución de luz bajo el dosel de Abies religiosa en Aquixtla, Puebla.
| G | dc | N | FAD | IAF | DAS | NSD | DPS | SF | SP | SG | CL | GLI | TransDir | TransDif | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| G | 1.00 | ||||||||||||||
| dc | 0.76 | 1.00 | |||||||||||||
| 0.001 | |||||||||||||||
| NSD | -0.47 | -0.88 | 1.00 | ||||||||||||
| 0.066 | 0.000 | ||||||||||||||
| FAD | 0.37 | 0.53 | -0.48 | 1.00 | |||||||||||
| 0.157 | 0.035 | 0.061 | |||||||||||||
| IAF | -0.37 | -0.49 | 0.42 | -0.97 | 1.00 | ||||||||||
| 0.162 | 0.056 | 0.109 | 0.000 | ||||||||||||
| DAS | 0.47 | 0.59 | -0.54 | 0.91 | -0.89 | 1.00 | |||||||||
| 0.068 | 0.017 | 0.031 | 0.000 | 0.000 | |||||||||||
| NSD | 0.29 | 0.10 | 0.12 | 0.27 | -0.32 | 0.35 | 1.00 | ||||||||
| 0.276 | 0.723 | 0.653 | 0.016 | 0.004 | 0.002 | ||||||||||
| DPS | 0.17 | 0.43 | -0.52 | 0.70 | -0.65 | 0.71 | -0.35 | 1.00 | |||||||
| 0.540 | 0.095 | 0.041 | 0.000 | 0.000 | 0.000 | 0.001 | |||||||||
| SF | -0.23 | -0.39 | 0.47 | -0.74 | 0.70 | -0.78 | 0.16 | -0.85 | 1.00 | ||||||
| 0.382 | 0.133 | 0.064 | 0.000 | 0.000 | 0.000 | 0.167 | 0.000 | ||||||||
| SP | 0.28 | 0.38 | -0.45 | 0.71 | -0.68 | 0.76 | -0.06 | 0.73 | -0.98 | 1.00 | |||||
| 0.303 | 0.141 | 0.078 | 0.000 | 0.000 | 0.000 | 0.607 | 0.000 | 0.000 | |||||||
| SG | -0.01 | 0.21 | -0.37 | 0.58 | -0.54 | 0.62 | -0.34 | 0.82 | -0.85 | 0.77 | 1.00 | ||||
| 0.956 | 0.438 | 0.153 | 0.000 | 0.000 | 0.000 | 0.002 | 0.000 | 0.000 | 0.000 | ||||||
| CL | 0.13 | 0.45 | -0.46 | 0.44 | -0.39 | 0.39 | -0.38 | 0.79 | -0.39 | 0.22 | 0.39 | 1.00 | |||
| 0.632 | 0.083 | 0.072 | 0.000 | 0.000 | 0.000 | 0.001 | 0.000 | 0.000 | 0.051 | 0.000 | |||||
| GLI | 0.36 | 0.50 | -0.47 | 0.95 | -0.94 | 0.95 | 0.26 | 0.75 | -0.76 | 0.72 | 0.63 | 0.46 | 1.00 | ||
| 0.176 | 0.050 | 0.069 | 0.000 | 0.000 | 0.000 | 0.020 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | ||||
| TransDir | 0.40 | 0.52 | -0.49 | 0.87 | -0.87 | 0.96 | 0.28 | 0.73 | -0.76 | 0.72 | 0.63 | 0.43 | 0.98 | 1.00 | |
| 0.120 | 0.040 | 0.056 | 0.000 | 0.000 | 0.000 | 0.013 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | |||
| TransDif | 0.30 | 0.47 | -0.44 | 0.98 | -0.97 | 0.89 | 0.23 | 0.72 | -0.73 | 0.68 | 0.59 | 0.47 | 0.97 | 0.90 | 1.00 |
| 0.262 | 0.069 | 0.091 | 0.000 | 0.000 | 0.000 | 0.044 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 |
G = área basal; dc = diametro cuadrático, N = densidad; FAD = fracción de apertura del dosel; IAF = índice de área foliar; DAS = duración acumulada diaria; NSD = haces de luz diarios; DPS = duración media de un haz de luz; SF = haces de luz “sunflecks”; SP = parches de luz; SG = claros de luz; CL = claros; GLI = índice de transmisión de luz; TransDiff = luz difusa transmitida; TransDir = luz directa transmitida.
Figura 1 Relación entre el diámetro cuadrático (dc) con: a) fracción de apertura del dosel (FAD); b) duración acumulada de radiación (DAS); c) radiación directa transmitida (TransDir); y d) índice de transmisión de luz (GLI); y relación entre la densidad del rodal (N) con: e) duración acumulada promedio de radiación (DAS) y f) duración promedio de un haz de luz (DPS).
La fracción de apertura del dosel (FAD) promedio (± desviación estándar) fue de 0.3059 ± 0.0946, (C.V. = 30.92%, máximo = 0.4792 y mínimo = 0.1262); mientras que para el índice de transmisión de luz (GLI) el promedio fue de 32.47% ± 10.65%, con un valor mínimo de 10.84, máximo de 53.48% y C.V. de 32.79%. La relación entre la estructura del dosel y transmisión de luz puede ser explicada por medio de una regresión lineal simple entre el grado de apertura del dosel (FAD) y la cantidad de luz transmitida (GLI) a través de este (Fig. 2). El coeficiente de determinación para el modelo de regresión lineal fue de 0.95 (Tabla 3).
La línea punteada indica la relación 1:1, y la línea continua en color rojo corresponde a la ecuación de regresión lineal GLI = 106.5676 × FAD + 0.1274 (n = 80; R² = 0.90).
Figura 2 Relación entre el índice de transmisión de luz y la fracción de apertura del dosel.
Los estadísticos para el índice de área foliar estimada a partir de las FH son los siguientes: media de 1.38 m2 m-2 ± 0.41 m2 m-2 y valores máximo y mínimo de 2.38 m2 m-2 y 0.76 m2 m-2, respectivamente (C.V. de 29.72%). El IAF se estimó a partir de características del dosel y de transmisión de luz, por lo que muestra una correlación negativa con la luz difusa transmitida (TransDif) con valor de -0.97, con la luz directa a través del dosel (TransDir) de -0.87, fracción de apertura del dosel (FAD) de -0.97 y de -0.94 con el índice de transmisión de luz (GLI) (Tabla 3). La estimación del IAF muestra que la relación existente entre los parámetros TransDif, TransDir, FAD y GLI se ajusta a un modelo de tipo exponencial (Gower y Normanm, 1991), para cada uno de los parámetros analizados (Fig. 3), lo cual indica, que la transmisión de la luz hacia la parte baja del dosel disminuirá de manera exponencial al aumentar el área de intercepción (Chazdon y Pearcy, 1991; Hardy et al., 2004; López-Sandoval et al., 2015).
Figura 3 Relación entre el índice de área foliar (IAF) con a) radiación difusa transmitida a través del dosel (TransDif); b) radiación directa transmitida (TransDir); c) fracción de apertura del dosel (FAD); y d) índice de transmisión de luz GLI (%).
El análisis estadístico mostro diferencias significativas (p < 0.05) entre sitios para la radiación transmitida a través del dosel, tanto directa como difusa; de tal manera, que la radiación directa incidente en el sotobosque bajo un dosel con una fracción de apertura de 0.19, será aproximadamente 63% menor que la radiación incidente en el sotobosque bajo un dosel con una fracción de apertura de 0.44. De igual forma, la radiación difusa incidente en el sotobosque se reducirá en alrededor de 57% en estas mismas condiciones de apertura del dosel (Tabla 2). Este patrón se explica dado que una mayor cobertura del dosel implica una mayor tasa de captura de radiación solar, expresado ello mediante el índice de transmisión de luz (GLI) (Marthews, Burslem, Phillips y Mullins, 2008; Messier et al., 1998; Promis et al., 2009).
Para la duración acumulada media de haces de luz (DAS), que describe el número de minutos de radiación directa potencial que podría recibir un sitio durante un día a lo largo de un año, se obtuvo un promedio de 224.4 min día-1 ± 76.4 min día-1, con un máximo de 377.8 min día-1 y un mínimo de 60.2 min día-1 (C.V. = 34.04%). La DAS se relacionó de manera positiva con FAD (R 2 = 0.82, p < 0.0001), con el índice de transmisión de luz (GLI) (R 2 = 0.91, p < 0.0001), y con la radiación directa transmitida (TransDir) (R 2 = 0.93, p<0.0001) (Fig. 4).
Figura 4 Relación entre la duración promedio acumulada de haces de luz por día en minutos (DAS) con: a) fracción de apertura del dosel (FAD); b) índice de transmisión de luz en porcentaje (GLI); c) radiación directa transmitida a través del dosel (TransDir); y d) índice de área foliar (IAF m2m-2).
La duración promedio de un haz de luz (DPS) es de 7.15 min ± 2.88 min, con un valor máximo de 18.26 min, mínimo de 3.38 min y C.V. de 40.27% (Tabla 2). Es conveniente señalar que el ambiente lumínico bajo diferentes tipos de dosel forestal tiende a ser diferente (Chazdon y Pearcy, 1991; Chen, 2014; Smith y Berry, 2013). El análisis de regresión muestra que DPS se relaciona de manera positiva con la FAD (R 2 = 0.49), con la transmisión de radiación a través del dosel (GLI) (R 2 = 0.43), con TransDir (R 2 = 0.54) (Fig. 5) y con TransDif (R 2 = 0.52), lo cual indica que con un incremento en cualesquiera de estos parámetros se incrementará la duración de los haces o destellos de luz bajo el dosel. La duración media de un haz de luz (DPS) tiene una relación negativa con el índice de área foliar (IAF) (R 2 = 0.43, p < 0.0001) (Fig. 5).
Figura 5 Relación existente entre la duración promedio de haces de luz (DPS) con: a) fracción de apertura del dosel (FAD); b) índice de transmisión de luz (GLI); c) radiación directa transmitida a través del dosel (TransDir); d) radiación difusa transmitida (TransDif) y e) índice de área foliar (IAF, m2m-2).
El análisis de distribución de frecuencias por clase de la duración de los haces de luz permite una mejor comprensión de la estructura del dosel (Smith, Sampson y Long, 1991) (Tabla 2). Los estadísticos obtenidos para el porcentaje del total de radiación solar directa incidente bajo el dosel (i.e. haces de luz, SF) son: media de 83.5% ± 6.26%, máximo y mínimo de 95.76% y 70.26%, respectivamente, (C.V. = 7.52%); para los parches o manchas de luz (SP) fue de 15.10% ± 5.20% con un máximo de 26.27% y mínimo de 4.23% (C.V. = 34.46%), mientras que los menores aportes al total de radiación directa presente al nivel del sotobosque fueron de las clases claros de luz (SG) y claros (CL) cuyos valores no cumplieron el supuesto de normalidad.
El análisis de regresión lineal mostró que la incidencia de haces o destellos de luz (SF), con duración menor a 8 min, se relaciona de manera negativa con FAD (R 2 = 0.55) y con la transmisión de radiación GLI (R 2 = 0.58), TransDir (R 2 = 0.58) y TransDif (R 2 = 0.53), y de manera positiva con el índice de área foliar (IAF) (R 2 = 0.50) (Fig. 6). En el caso de los parches o manchas de luz, así como de los claros de luz, los cuales se relacionan de manera negativa con el índice de área foliar (valores de r = -0.68 y r = -0.54, respectivamente) y de forma positiva con la transmisión de radiación (GLI, TranDir y TransDif) y con el FAD, los resultados indican que para un dosel en el cual se incrementa la presencia de claros, de igual forma se incrementará la transmisión de luz a través de este, estableciéndose un dosel discontinuo en el cual los eventos lumínicos incrementarán su duración (Hardy et al., 2004; Marthews et al., 2008; Sapkota y Odén, 2009). Para la radiación solar directa proveniente de eventos lumínicos de más de 120 minutos de duración que inciden en el dosel forestal (claros, CL), se obtuvo una correlación débil (r < 0.50) tanto con los parámetros de estructura del dosel como para aquellos relacionados con la transmisión de la radiación a través del dosel (Tabla 3).
Figura 6 Relación entre la fracción de haces de luz (SF) con respecto del total de radiación directa presente bajo el dosel con: a) fracción de apertura del dosel (FAD), b) índice de transmisión de luz (GLI), c) radiación directa transmitida a través del dosel (TransDir), d) radiación difusa transmitida a través del dosel (TransDif) y e) índice de área foliar.
Discusión
Bajo el supuesto de que la composición y estructura de los rodales son determinantes de las comunidades de plantas presentes en el sotobosque (Fahey y Puettmann, 2007; Macdonald y Fenniak, 2007), se ha supuesto que existe influencia directa del dosel forestal sobre la riqueza y diversidad vegetal; por lo anterior, a partir de variaciones en la densidad, expresada en términos de la apertura, transmisión de luz, índice de área foliar, o bien, del régimen lumínico, se afectará de manera determinante a las comunidades vegetales que integran al sotobosque (Macdonald y Fenniak, 2007; Roberts y Gilliam, 1995).
Los parámetros de estructura del rodal obtenidos en este estudio muestran que un incremento en el diámetro cuadrático implica una reducción en la densidad, lo que significará un incremento en la apertura del dosel (Beaudet et al., 2011) y, en consecuencia, una mayor duración e incidencia de radiación directa a nivel del sotobosque. A raíz de lo anterior, aunque el diámetro cuadrático podría ser empleado como un estimador de la capacidad del dosel de transmitir la radiación incidente hacia el sotobosque, es evidente la necesidad de ampliar las estimaciones considerando las diferentes etapas de desarrollo de los rodales y poder expresar una conclusión más sólida a este respecto. Se encontró también una relación negativa entre la densidad del rodal y los parámetros de duración de eventos lumínicos (DAS y DPS), lo que indica que un incremento en la densidad del rodal disminuirá la variabilidad espacial de aperturas en el dosel (Beaudet et al., 2011), con un efecto importante en la duración de los haces de luz incidentes en el sotobosque. Este resultado es consistente con otros estudios en bosques de coníferas y boreales, donde se han encontrado relaciones negativas entre la densidad del rodal con disponibilidad y transmisión de luz bajo el dosel (e.g. radiación difusa) (Comeau, Heineman y Newsome, 2006; Comeau, 2001).
Respecto a la relación entre fracción de apertura del dosel e índice de transmisión de luz, el resultado concuerda con otros modelos de transmisión de luz en los cuales se ha utilizado a la cobertura como una variable importante para explicar la transmisión de luz a través del dosel forestal; por ejemplo, utilizando el densitómetro esférico se han obtenido valores de R 2 en un intervalo de 0.71 a 0.99 para bosques mezclados y puros de clima templado, Fagus y Quercus (Baudry, Charmetant, Collet, y Ponette, 2014), Quercus y Pinus (Buckley, 1999) y Betula (Comeau, Gendron y Letchford, 1998; Johansson, 1996); así como para bosques de clima tropical dominados por Cedrella, Cordia, Erythrina, Eucalyptus, Hieronyma, Musa, Pouteria y Psidium (Bellow y Nair, 2003). De igual manera, en modelos de transmisión de luz en los cuales se han usado fotografías hemisféricas para estimar la estructura del dosel, se han obtenido valores del coeficiente de determinación mayores a 0.60 para bosques de Pseudotsuga, Tsuga, Thuja y Alnus (Frazer et al., 2001), Acer, Betula y Fagus (Beaudet y Messier, 2002), plantaciones de Eucalyptus, Pinus, Cupressus, Grevillea y Juniperus (Yirdaw y Luukkanen, 2004), así como Liriodendron y Liquidambar (Lhotka y Loewenstein, 2006).
Estos resultados indican que un incremento en la disponibilidad de luz bajo el dosel forestal, dependerá básicamente del incremento de claros en el dosel (Canham, 1988; Gray, Spies y Easter, 2002; Lhotka y Loewenstein, 2006). Los valores para el IAF, estimados en este trabajo, se encuentran dentro del intervalo de valores obtenido por medio de otros métodos indirectos en bosques de Abies religiosa; se han obtenido valores que van de 3.0 m2 m-2 a 4.5, para bosques densos, y de 1.0 m2 m-2 a 3.0 m2 m-2, para aquellos de menor densidad, usando datos radiométricos de imágenes de satélite ajustados mediante fotografías hemisféricas (Velasco, Champo, España y Baret, 2010), así como valores de IAF de entre 3.3 m2 m-2 ± 0.18 m2 m-2 y 2.85 m2 m-2 ± 0.35 m2 m-2 obtenidos a partir del uso de sensores de radiación fotosintéticamente activa o PAR (Pérez-Suárez, Fenn, Cetina-Alcalá y Aldrete, 2008).
Los resultados obtenidos para la duración promedio acumulada de haces de luz (DAS) sugieren que, a medida que se incrementa la presencia de claros en el dosel, se aumenta la cantidad de luz directa bajo el mismo (Hardy et al., 2004), lo cual influye en el índice de transmisión de luz (Frazer et al., 1999), como resultado de un incremento en la duración promedio de los haces de luz que podrían alcanzar algún punto en el sotobosque (Lambers, Chapin y Pons, 2008). El índice de área foliar (IAF) presentó una correlación de -0.79 (p < 0.0001) con DAS (Fig. 4) debido, principalmente, a que un incremento en el IAF representa a un dosel más denso (Chazdon y Pearcy, 1991), que da lugar a la presencia de aperturas de menor tamaño (Smith y Berry, 2013) y reduce la penetración de radiación directa a través del dosel, así como la duración de los haces de luz a nivel del sotobosque. La duración media de un haz de luz (DPS) mostró una relación negativa con el índice de área foliar (IAF) (R 2 = 0.42, p < 0.0001) (Fig. 5), manifestándose nuevamente que la estructura de las copas juega un rol importante en la dinámica del ambiente lumínico, así como en el comportamiento de los haces de luz que serán transmitidos a través del dosel y alcanzarán el sotobosque (Marthews et al., 2008). Así, cualquier incremento en la trasmisión de luz directa a través del dosel, ya sea por una mayor cantidad y tamaño de claros, o bien, por la disminución del follaje (es decir, del IAF), afectará de manera positiva a la duración de los haces de luz que alcanzarán el sotobosque.
Lo anterior es relevante porque se ha encontrado en el caso de Abies magnifica la incidencia de haces de luz de larga duración y alta densidad de flujo de fotones (i.e. mayores a 2000 µmol m-2 s-1), que ocasionan fotoinhibición (Wong, Chen, Huang y Weng, 2012), en condiciones bajo dosel y con presencia de claros grandes; los haces de luz reducen tanto la germinación como la sobrevivencia, aunque también incrementan las tasas de reclutamiento en áreas con incidencia de haces de luz de corta duración y bajo doseles densos (Ustin, Woodward, Barbour y Hatfield, 1984; Way y Pearcy, 2012). Ustin et al. (1984) sugieren que el patrón de distribución agregado para plántulas de esta especie es el resultado de la influencia de la exposición prolongada a haces de luz con efecto en los lapsos de estrés hídrico y térmico a los cuales están sujetas las plántulas. La distribución de frecuencias muestra que más de 70% de la radiación que ingresa en forma directa y se encuentra disponible en el sotobosque, se realiza a partir de haces de luz o pulsos de radiación directa de corta duración (< 8 min) (Tabla 2). Esto coincide con evaluaciones de la contribución de los haces de luz, en los cuales se determinó que estos contribuyen entre 50% y 70% del total de la luz disponible en el sotobosque en bosques tropicales (Chazdon y Pearcy, 1991). En el caso de bosques de coníferas siempre verdes, la composición, estructura, así como la etapa de sucesión es muy heterogénea, por lo cual la caracterización de los regímenes de luz es difícil (Canham y Burbank, 1994; Chazdon et al., 1996). Sin embargo, se considera que la cantidad de radiación solar disponible bajo el dosel de bosques de coníferas perennes es mayor que aquella presente bajo el dosel de bosques tropicales siempre verdes, debido a que existen menos estratos de árboles y doseles moderadamente abiertos (Canham et al., 1990; Chazdon, 1988; Chen, 2014). Esta elevada presencia de eventos lumínicos de corta duración incrementa la heterogeneidad del ambiente lumínico a nivel del sotobosque (Chazdon y Pearcy, 1991), lo que provoca variaciones rápidas en la densidad del flujo de fotones (PFD), en duración, frecuencia e intensidad, las cuales afectan la respuesta fotosintética de las plantas a nivel del sotobosque e inducen respuestas metabólicas bioquímicas y fisiológicas de algunas especies (Pearcy, Krall y Sassenrath-Cole, 1996). En otros casos, existe un efecto negativo debido a fotoinhibición y reducción de la ganancia de carbono (Le Gouallec, Cornic, y Blanc, 1990; Leakey, Press y Scholes, 2003), o marchitez y daño en hojas por temperaturas elevadas (Chazdon y Pearcy, 1991; Wagner y McGraw, 2013).
En general, la fracción de apertura del dosel, el índice de área foliar, la cantidad de luz directa (TransDir) y la cantidad de luz difusa (TransDif) en los rodales analizados son muy variables (Tabla 2), lo cual indica que tanto la estructura como la transmisión de luz a través del dosel es muy heterogénea, dando como resultado un ambiente lumínico variable en el sotobosque de los rodales analizados.
Conclusiones
El uso de FH para estimar el régimen lumínico en un rodal, o para estimar atributos estructurales del dosel, requiere de la ejecución meticulosa de una serie de pasos (i.e. selección de equipo, obtención y procesamiento de imágenes, aplicación de algoritmos), además de considerar posibles limitantes como la dependencia a condiciones meteorológicas adecuadas para la obtención de imágenes de calidad (i.e. días con nublado uniforme). También, es importante la limitada discriminación espectral, ya que los valores obtenidos para los colores rojo, verde y azul (RGB) no son considerados medidas precisas de radiación, comparadas con instrumentos especializados en espectrometría de luz. Además, se debe tener en cuenta que los valores obtenidos para el IAF a partir de FH están subestimados en comparación con métodos directos empleados para ello, y que es necesario considerar que la imagen obtenida es una representación simplificada en dos dimensiones del dosel, siendo este un complejo tridimensional. Es recomendable que el usuario se encuentre familiarizado con los modelos matemáticos empleados en el análisis, de tal manera que la aplicación así como la interpretación de los resultados del análisis de las FH sea el adecuado.
Los resultados obtenidos en este estudio indican que la transmisión de radiación solar hacia el sotobosque en los bosques de Abies religiosa estudiados se caracteriza por una gran variación espacial y temporal. Los parámetros de estructura del dosel -principalmente- y del rodal -en menor grado- explican, al menos en parte, la variación del régimen de radiación solar.
El mejor predictor para el índice de área foliar, radiación directa y radiación difusa disponible en el sotobosque y del índice de transmisión de luz del dosel fue la fracción de apertura del dosel, mientras que las estimaciones a partir de algún parámetro de la estructura del rodal mostraron una relación débil.
La estimación del índice de área foliar a partir de parámetros de transmisión de luz y estructura del dosel por medio del modelo Beer-Lambert, mostró que la relación existente entre los parámetros de radiación trasmitida (directa y difusa), apertura del dosel e índice de transmisión de luz del dosel se ajustan a un modelo de tipo exponencial, indicando que la transmisión de luz disminuirá siguiendo dicha tendencial al aumentar el área de intercepción del dosel, debido a un incremento en la cobertura y densidad de copas. El modelo de Beer-Lambert resulta ser adecuado para doseles simples y homogéneos, como el caso de rodales puros; sin embargo, investigaciones recientes indican que para doseles heterogéneos y complejos dicho modelo podría no proporcionar resultados satisfactorios y es necesario utilizar otro tipo de modelos (Forrester et. al, 2014). Los modelos de regresión ajustados para intentar explicar y estimar de manera sencilla los componentes del régimen y dinámica de la radiación solar presente en el sotobosque, a partir de parámetros de estructura del dosel, mostraron una relación moderada. Los parámetros de apertura el dosel e índice de transmisión de luz pueden ser usados para explicar buena parte de la variación en la duración y ración de la incidencia de haces de luz (e.g. duración media acumulada diaria, haces de luz “sunflecks”, parches de luz “sunpatch” y claros de luz “sungap”). La capacidad de los parámetros de estructura del rodal (área basal, densidad y diámetro cuadrático) para estimar los componentes del régimen y dinámica de la radiación solar bajo el dosel fue pobre o nula.
