Servicios Personalizados
Revista
Articulo
Indicadores
-
Citado por SciELO
-
Accesos
Links relacionados
-
Similares en SciELO
Compartir
Revista Chapingo serie ciencias forestales y del ambiente
versión On-line ISSN 2007-4018versión impresa ISSN 2007-3828
Rev. Chapingo ser. cienc. for. ambient vol.28 no.1 Chapingo ene./abr. 2022 Epub 02-Feb-2024
https://doi.org/10.5154/r.rchscfa.2020.12.077
Artículo científico
Estequiometría de la caída de acículas de Pinus hartwegii Lindl. en dos bosques alpinos del centro de México
1 Colegio de Postgraduados, Campus Montecillo, Posgrado Forestal. Carretera México-Texcoco km 36.5. C. P. 56230. Texcoco, Estado de México, México.
2 Colegio de Postgraduados, Campus Montecillo, Posgrado en Edafología. Carretera México-Texcoco km 36.5. C. P. 56230. Texcoco, Estado de México, México.
3 Instituto Nacional de Investigaciones Forestales, Agrícolas y Pecuarias (INIFAP), Centro Nacional de Investigación Disciplinaria en Conservación y Mejoramiento de Ecosistemas Forestales (CENID-COMEF). Av. Progreso núm. 5, Barrio de Santa Catarina. C. P. 04010. Coyoacán, Ciudad de México, México.
Introducción:
Es imprescindible tener líneas base sobre la dinámica de nutrientes en los bosques, debido a las alteraciones que el cambio climático pueda causar.
Objetivo:
Cuantificar la producción anual de acículas de Pinus hartwegii Lindl. y proporción de nutrientes en los bosques alpinos cerro Jocotitlán (JO) y cerro Tláloc (TL), Estado de México.
Materiales y métodos:
En cada bosque se colocaron 12 trampas circulares (30 cm de diámetro) para acícula, a nivel del suelo, distribuidas en cuatro sitios topográficamente contrastantes. Durante un año, se realizaron 228 mediciones de masa foliar y 1 140 determinaciones químicas para conocer la estequiometría de acículas. Las mediciones se sometieron a un análisis de varianza longitudinal, probando las tendencias en el tiempo (P < 0.05).
Resultados y discusión:
La producción de acículas en JO fue 67 % superior (11.2 Mg∙ha-1∙año-1) que en TL (6.7 Mg∙ha-1∙año-1); la caída fue mayor durante el verano (junio y julio, meses con mayor precipitación). En JO, el flujo de nutrientes fue 98.0, 5.2, 8.7, 24.6 y 5.6 kg∙ha-1∙año-1 para N, P, K, Ca y Mg, respectivamente; para TL fue 55.3, 3.4, 7.8, 14.4 y 4.7 kg∙ha-1∙año-1 en el mismo orden de nutrientes. Las concentraciones de nutrientes fueron más bajas de marzo a mayo. Excepto por K, las concentraciones de nutrientes y producción de acículas mostraron tendencias estacionales cuadráticas y cúbicas. La dinámica del Mg y proporciones N:Mg y N:K en TL fueron más favorables para el crecimiento arbóreo.
Conclusiones:
Los bosques Jocotitlán y Tláloc son muy productivos en acícula (comparados con otros ecosistemas) y dinámicos en la transferencia de nutrientes.
Palabras clave: masa foliar; piso forestal; flujo de nutrientes; pino; cambio climático.
Introduction:
It is essential to have baselines on nutrient dynamics in forests, due to disturbances that climate change may cause.
Objective:
To quantify the annual production of needles of Pinus hartwegii Lindl. and the proportion of nutrients in the alpine forests of Jocotitlán (JO) and Tláloc (TL) mountains, Estado de México.
Materials and methods:
A total of 12 circular needle litter traps (30 cm diameter) were placed at ground level, in each forest, distributed in four topographically contrasting sites. For one year, 228 leaf mass measurements and 1 140 chemical determinations were made to determine needle stoichiometry. Measurements were subjected to a longitudinal analysis of variance, by testing trends over time (P < 0.05).
Results and discussion:
Needle production in JO were 67 % higher (11.2 Mg∙ha-1∙year-1) than in TL (6.7 Mg∙ha-1∙year-1); needle litterfall was higher during summer (June and July, months with higher precipitation). For JO, nutrient flux was 98.0, 5.2, 8.7, 24.6, and 5.6 kg∙ha-1∙year-1 for N, P, K, Ca, and Mg, respectively; for TL it was 55.3, 3.4, 7.8, 14.4, and 4.7 kg∙ha-1∙year-1 in the same order of nutrients. Nutrient concentrations were lower from March to May. Except for K, nutrient concentrations and needle production showed quadratic and cubic seasonal trends. Mg dynamics and N:Mg and N:K ratios in TL were more positive for tree growth.
Conclusions:
Jocotitlán and Tláloc forests produce significant needle mass (compared to other ecosystems) with high dynamic in nutrient transfers.
Keywords: leaf mass; forest floor; nutrient flux; pine; climate change
Introducción
La caída de follaje en los bosques representa la principal vía de transferencia y retorno de nutrientes al suelo (Sayer et al., 2020); a su vez, la masa foliar en el piso forestal (mantillo) constituye el reservorio clave para la liberación gradual de nutrientes hacia el suelo mineral. La actividad fisiológica estacional de los árboles, longevidad de las hojas y el clima influye en la caída de la hoja (Gavinet, Ourcival, & Limousin, 2019), la cual da lugar a un flujo estacional de nutrientes, característico para cada ecosistema y su fase de sucesión. No obstante, la alteración de ciclos del carbono y del nitrógeno puede conducir a procesos de fertilización por el incremento de CO2 atmosférico en los bosques, estimulando el crecimiento de la masa foliar y modificando su composición química, pero su efecto dependerá de la limitación progresiva de los nutrientes (González de Andrés, 2019; Silva, Gómez‐Guerrero, Doane, & Horwath, 2015). Aunque localmente la fertilidad del suelo y el microclima del bosque influyen en los niveles de caída de hoja y sus nutrientes (Cheng, Lee, Lee, Chen, & Menyailo, 2020; Wood, Lawrence, & Clark, 2006), indirectamente, los eventos climáticos globales como El Niño-Oscilación del Sur (ENSO) y la Oscilación del Atlántico Norte (NOA) afectan geográficamente la masa y estequiometría de la caída de follaje (González de Andrés et al., 2019).
A nivel global, el estudio del flujo de nutrientes vía caída de follaje es limitado y requiere de más investigación (Neumann et al., 2018). En México existen pocos trabajos sobre la dinámica de la masa foliar y concentración de nutrientes en los bosques. Algunos se han realizado en la mesa central y norte de México, en bosques de pino-encino y matorrales con climas semisecos y elevaciones menores de 2 800 m (González-Rodríguez et al., 2018; Pérez-Suárez, Arredondo-Moreno, Huber-Sannwald, & Vargas-Hernández, 2009); sin embargo, la información en bosques alpinos, ecosistemas particularmente vulnerables al cambio climático, sigue pendiente. En consecuencia, se recurre a la experiencia de trabajos realizados en otras partes del mundo, para explicar el comportamiento de nuestras áreas forestales. Debido a lo anterior, y considerando que la variación climática influye de manera diferente en cada región, se requieren estudios en ecosistemas representativos del país, para saber cómo responderán los ecosistemas forestales a dichos cambios (Gómez-Guerrero & Doane, 2018).
Como resultado del cambio climático se esperan variaciones en la producción de masa foliar y en el flujo de nutrientes en los bosques (Gómez-Guerrero & Doane, 2018). El aumento pronosticado de temperatura podría resultar en estrés hídrico y cambios en la distribución de nutrientes en los reservorios de la vegetación (Neumann et al., 2018). Por tanto, es importante contar con estudios base sobre la dinámica de nutrientes en los ecosistemas forestales para valorar sus cambios futuros (Ge & Xie, 2017).
El objetivo de este estudio fue evaluar la producción anual de acículas, su concentración y masa de nutrientes (N, P, K, Ca y Mg) en bosques alpinos de Pinus hartwegii Lindl. del cerro Jocotitlán (JO) y cerro Tláloc (TL) del Estado de México. Este estudio subsana la carencia de información sobre dinámica de nutrientes en bosques alpinos y establece una línea base para medir los cambios que puedan ocurrir como respuesta de los ecosistemas ante el cambio climático. La información generada también es útil para fines de modelación de procesos biogeoquímicos en la vegetación forestal.
Materiales y métodos
Sitios de estudio
El trabajo se llevó a cabo en las montañas Jocotitlán (JO) y Tláloc (TL) en el Estado de México (Figura 1). Las montañas son parte de la Faja Volcánica Transmexicana compuesta geológicamente por un complejo de andesitas-dacitas y depósitos piroclásticos (Arce et al., 2019). El clima es tipo C(w2)(w), semifrío húmedo con lluvias en verano y menos de 5 % de lluvia invernal, con cociente precipitación/temperatura mayor de 55 (García, 2004). La precipitación media anual varía de 900 a 1 200 mm con temperatura media anual entre 5 y 10 °C. La masa forestal en JO se compone en promedio por 180 árboles∙ha-1 con un área basal de 14 m2∙ha-1 (Núñez-García, Gómez-Guerrero, Terrazas-Salgado, Vargas-Hernández, & Villanueva-Díaz, 2020), mientras que en TL, se forma por 250 árboles∙ha-1 y 27 m2∙ha-1 (Acosta Mireles, Carrillo Anzures, Delgado, & Velasco Bautista, 2014). El incremento en área basal de los árboles en TL es en promedio 2.5 mayor que en JO (Correa‐Díaz et al., 2020; Núñez et al., 2020).
Después de recorridos de campo y revisión de material cartográfico y satelital, en cada montaña se establecieron cuatro sitios de colecta periódica para la caída de acícula. Para captar la variación de cada montaña se incluyeron sitios de colecta contrastantes en elevación y exposición, pero con masas puras de P. hartwegii. La ubicación de sitios se muestra en el Cuadro 1.
![](/img/revistas/rcscfa/v28n1//2007-4018-rcscfa-28-01-57-gf5.png)
Figura 1 Localización de los sitios de estudio en las montañas de Jocotitlán (a) y Tláloc (b), Estado de México.
Cuadro 1 Sitios utilizados en la colecta periódica de acículas de Pinus hartwegii en las montañas de Jocotitlán y Tláloc, Estado de México.
Montaña | Sitio | Elevación (m) | Exposición | Latitud (°) | Longitud (°) |
---|---|---|---|---|---|
Jocotitlán | NO-3800 | 3 800 | Noroeste | 19.7423 | -99.7645 |
SO-3800 | 3 800 | Suroeste | 19.7412 | -99.7662 | |
NO-3700 | 3 700 | Noroeste | 19.7427 | -99.761 | |
SO-3700 | 3 700 | Suroeste | 19.7397 | -99.768 | |
Tláloc | NO-3900 | 3 900 | Noroeste | 19.424 | -98.7245 |
NO-3500 | 3 500 | Noroeste | 19.4161 | -98.7258 | |
SO-3900 | 3 900 | Suroeste | 19.7255 | -98.7255 | |
SO-3500 | 3 500 | Suroeste | 19.4045 | -98.7419 |
Masa de acículas
En cada sitio, la masa de acículas (kg∙ha-1) se cuantificó en un transecto de 20 m con tres trampas de 30 cm de diámetro, ubicadas a nivel del suelo y delimitadas con poliducto de 127 mm de grosor. La distancia entre trampas fue de 10 m, colocadas sobre curva de nivel. En total se establecieron 12 trampas en cada montaña (tres repeticiones x cuatro sitios de colecta).
Para identificar el material fresco depositado, se removió todo el mantillo dentro del círculo hasta exponer el suelo mineral, permitiendo que el contraste entre el color negro del suelo y el café de las acículas facilitara la identificación de material recientemente caído. Todas las trampas se ubicaron debajo del dosel del arbolado, procurando coberturas de al menos 60 %. Las acículas acumuladas se colectaron durante 9 (JO) y 10 (TL) visitas de junio del 2017 a junio del 2018. Las acículas se secaron en estufa de aire de circulación forzada (Sheldon Fx14-S, EUA) a 70 °C por 72 h y se pesaron en una balanza digital (Ohaus Scout Pro, EUA) con precisión de 0.01g.
Nutrientes del mantillo
Al inicio del experimento se tomaron tres muestras de mantillo en cada sitio (12 por montaña) y se analizaron las concentraciones de N, P, K, Ca y Mg, para estimar los reservorios de nutrientes en el piso forestal. De la misma forma, y para correlacionar el contenido de nutrientes de la caída con el del piso forestal y el suelo mineral, se determinaron las concentraciones de N y P y los cationes intercambiables K+, Ca+2 y Mg+2 del suelo superficial (0 a 10 cm), tomando 36 muestras en JO y 40 en TL durante el año de observación.
En cada sitio se instaló un registrador electrónico de temperatura y humedad relativa (Datalogger, HoBo®, EUA), para obtener los datos a intervalos de 4 h y conocer la variación en las condiciones climáticas de cada montaña.
Análisis de las muestras
Las muestras secas de mantillo y acículas se molieron para el análisis de concentración de nutrientes. El N en el tejido vegetal se determinó mediante el método semimicro-Kjeldahl (Bremmer, 1965). Las concentraciones de P, K, Ca y Mg se determinaron en el extracto resultante de la digestión con HNO3:HClO4 (2:1, v:v) de tejido vegetal seco y molido, con un espectrofotómetro de emisión atómica con plasma acoplado inductivamente (ICP-OES 725 Series, Agilent) (Alcántara & Sandoval, 1999). El pH de las muestras de suelo mineral se determinó con potenciómetro en una mezcla de suelo agua 1:2. También se determinó el nitrógeno total por el método de micro-Kjeldahl, el P por el método Bray y Kurtz (1945) y los cationes intercambiables con acetato de amonio 1 N pH 7.0.
Análisis estadístico
Se verificó que los datos cumplieran con el supuesto de normalidad con la prueba de Shapiro-Wilk y cuando no fue el caso se utilizaron transformaciones logarítmicas, raíz cuadrada o de tipo Box-Cox. Cuando el ajuste a la normalidad no fue posible mediante transformaciones, se utilizó la prueba no paramétrica de Wilcoxon (Zar, 2010), para comparaciones entre montañas. La masa de nutrientes (kg∙ha-1), la concentración (%, ppm o cmolc∙kg-1) y los flujos (kg∙ha-1∙año-1) fueron las variables respuesta para la comparación entre montañas. También se utilizaron proporciones de nutrientes con respecto al N, que es una referencia común para este tipo de estudios (Sardans et al., 2016). Para los flujos se hizo un análisis estadístico de las tendencias, ya que las muestras se repitieron en el tiempo en cada sitio; por lo anterior, las mediciones están correlacionadas en el tiempo. Con modelos lineares mixtos, se hizo un análisis de varianza longitudinal considerando el tiempo, montaña e interacción tiempo*montaña como factores fijos (Fitzmaurice, Laird, & Ware, 2004). Se probó si las tendencias temporales seguían un patrón lineal, cuadrático o cúbico (P > 0.05). Los análisis se realizaron con los paquetes estadísticos SAS (SAS Institute Inc., 2011) y R (3.6.3) (R Core Team, 2020).
Resultados
Acorde con el Cuadro 2, las características químicas del mantillo (N, P, K y Ca) entre ambas montañas fueron similares, excepto por la concentración de Mg que fue 37 % superior en TL (P > 0.0001). La concentración de Ca fue significativamente superior en 11 % en TL (P < 0.09). Aunque la concentración media de K es 54 % superior en TL, la diferencia no fue significativa entre montañas.
Cuadro 2 Concentraciones de nutrientes del mantillo (piso forestal) en dos bosques alpinos de Pinus hartwegii en el centro de México.
Sitio | NT | P | K | Ca | Mg |
---|---|---|---|---|---|
(%) | --------------------(ppm)-------------------------- | ||||
Jocotitlán | |||||
Media | 0.8 | 535.5 | 722 | 2 292.4 | 491 |
DE | 0.2 | 94.6 | 205.3 | 335.6 | 106.6 |
CV (%) | 21 | 17.7 | 28.4 | 14.6 | 21.7 |
n | 12 | 12 | 12 | 12 | 12 |
Tláloc | |||||
Media | 0.9 | 548.5 | 1 112.4 | 2 540.6 | 674.9 |
DE | 0.2 | 81.1 | 575 | 417.8 | 88.1 |
CV (%) | 27.8 | 14.8 | 51.7 | 16.4 | 13.1 |
n | 12 | 12 | 12 | 12 | 12 |
P< | ns | ns | ns | 0.0941 | 0.0001 |
NT = nitrógeno total; DE: desviación estándar de la media; CV: coeficiente de variación. El valor P< indica diferencia estadística entre montañas de acuerdo con la prueba de Wilcoxon; ns = no significativo.
En el Cuadro 3 se observa que el suelo mineral es muy ácido (pH ≈ 5.2) en ambas montañas. El N total, P extraíble y el K+ intercambiable son estadísticamente iguales entre montañas. Las diferencias estadísticas se reflejaron en Mg, Ca y en la suma de cationes (K+ + Ca+2 + Mg+2) con valores superiores en TL (P < 0.021).
Cuadro 3 Concentraciones de nutrientes del suelo superficial (0-10 cm) en dos bosques alpinos de Pinus hartwegii en el centro de México.
Montaña | pH (1:2) | NT (%) | P (ppm) | Cationes intercambiables (cmolc(+)∙kg-1) | |||
---|---|---|---|---|---|---|---|
K+ | Ca+2 | Mg+2 | Suma | ||||
Jocotitlán | |||||||
Media | 5.3 | 0.4 | 3.3 | 1.2 | 2.5 | 0.2 | 3.9 |
DE | 0.9 | 0.2 | 1.8 | 1.3 | 2 | 0.1 | 1.3 |
CV (%) | 16.5 | 42.1 | 53.1 | 106.8 | 81.2 | 55.2 | 34.1 |
n | 36 | 36 | 36 | 36 | 36 | 36 | 36 |
Tláloc | |||||||
Media | 5.1 | 0.4 | 3.4 | 1.4 | 3.6 | 0.24 | 5.2 |
DE | 0.5 | 0.1 | 2.2 | 1.6 | 2.1 | 0.3 | 1.3 |
CV (%) | 10 | 32.8 | 64.7 | 108.1 | 60.2 | 129.9 | 25.4 |
n | 40 | 40 | 40 | 40 | 40 | 40 | 40 |
P< | ns | ns | ns | ns | 0.0001 | 0.0211 | 0.0001 |
NT = nitrógeno total; DE: desviación estándar de la media; CV: coeficiente de variación. El valor P< indica diferencia estadística entre montañas de acuerdo con la prueba de Wilcoxon; ns = no significativo.
El Cuadro 4 muestra que tanto la masa del mantillo como la masa de sus nutrientes fue estadísticamente mayor en JO (P < 0.04). La masa de mantillo en JO (46.5 Mg∙ha-1) fue 175 % mayor que la de TL (16.9 Mg∙ha-1). La masa de nutrientes del mantillo fue superior en JO con valores superiores en 164, 181, 60, 164 y 118 % para N, P, K, Ca y Mg, respectivamente.
Cuadro 4 Masa del piso forestal (mantillo) y de nutrientes en dos bosques alpinos de Pinus hartwegii en el centro de México.
Sitio | Masa | N | P | K | Ca | Mg |
---|---|---|---|---|---|---|
(Mg∙ha-1) | -------------------------(kg∙ha-1)------------------------- | |||||
Jocotitlán | ||||||
Media | 46.5 | 360.6 | 24.7 | 34.5 | 106.8 | 23.8 |
DE | 21.4 | 194.4 | 11.5 | 18.3 | 53.2 | 12.7 |
CV (%) | 46 | 53.9 | 46.5 | 53.2 | 49.8 | 53.4 |
n | 12 | 12 | 12 | 12 | 12 | 12 |
Tláloc | ||||||
Media | 16.9 | 136.7 | 8.8 | 21.5 | 40.5 | 10.9 |
DE | 19.1 | 151.7 | 10.5 | 31.7 | 43.4 | 12.6 |
CV (%) | 113.2 | 111 | 119.5 | 147.8 | 107.1 | 115.3 |
n | 12 | 12 | 12 | 12 | 12 | 12 |
P< | 0.0014 | 0.0001 | 0.001 | 0.0404 | 0.0003 | 0.0141 |
DE: desviación estándar de la media; CV: coeficiente de variación. El valor P< indica diferencia estadística entre montañas de acuerdo con la prueba de Wilcoxon.
De acuerdo con el Cuadro 5, el flujo de masa de acículas y nutrientes fue estadísticamente superior en JO (P < 0.04). La transferencia de masa de acículas del dosel al piso forestal fue 67 % superior en JO con respecto a TL. En cuanto al flujo de nutrientes en masa, las diferencias fueron 77, 54, 12, 70 y 19 % superiores en JO con respecto a TL, para N, P, K, Ca y Mg, respectivamente.
Cuadro 5 Flujo anual de masa de acículas y nutrientes vía caída de follaje en dos bosques alpinos de Pinus hartwegii en el centro de México.
Sitio | Masa | N | P | K | Ca | Mg |
---|---|---|---|---|---|---|
(Mg∙ha-1∙año-1) | -----------------------(kg∙ha-1∙año-1)------------------- | |||||
Jocotitlán | ||||||
Media | 11.2 | 98 | 5.2 | 8.7 | 24.6 | 5.6 |
DE | 3 | 20.9 | 1.1 | 2.6 | 6.8 | 1.9 |
CV (%) | 26.8 | 21.3 | 23.1 | 30.7 | 27.7 | 33.4 |
n | 12 | 12 | 12 | 12 | 12 | 12 |
Tláloc | ||||||
Media | 6.7 | 55.3 | 3.4 | 7.8 | 14.4 | 4.7 |
DE | 3.4 | 28.3 | 1.8 | 6.6 | 5.8 | 2.2 |
CV (%) | 50.7 | 51.2 | 54.3 | 85.2 | 40.2 | 46.9 |
n | 12 | 12 | 12 | 12 | 12 | 12 |
P< | 0.0027 | 0.0001 | 0.001 | 0.0404 | 0.0003 | 0.0141 |
DE: desviación estándar de la media; CV: coeficiente de variación. El valor P< indica diferencia estadística entre montañas de acuerdo con la prueba de Wilcoxon.
Tendencias estacionales
La tendencia estacional más común fue cúbica. Las concentraciones de K no mostraron una tendencia definida, excepto cuando se consideró la proporción N:K. El Cuadro 6 indica que, en todos los casos en que se pudo ajustar una tendencia, el factor tiempo fue significativo (P < 0.001), mientras que el efecto de montaña no fue estadísticamente significativo para las concentraciones de P, Ca y la relación N:Ca. La interacción montaña*tiempo fue estadísticamente significativa para N y las relaciones N:K y N:Mg.
Cuadro 6 Significancia estadística del análisis de varianza longitudinal de las tendencias temporales de masa de acículas y concentración de nutrientes en en dos bosques alpinos de Pinus hartwegii en el centro de México.
Variable | Tendencia | Tiempo | Montaña | Interacción tiempo*montaña |
---|---|---|---|---|
Masa de acículas | Cuadrática | <.0001 | <.0001 | 0.0443 |
N | Cúbica | 0.0002 | 0.0002 | 0.0153 |
P | Cuadrática | <.0001 | 0.8172 | ns |
K | Sin tendencia | --- | --- | --- |
Ca | Cuadrática | <.0001 | 0.7372 | ns |
Mg | Cúbica | 0.0322 | <.0001 | ns |
N:P | Cúbica | 0.0007 | 0.0058 | ns |
N:K | Cúbica | <.0001 | <.0001 | 0.0007 |
N:Ca | Cuadrática | <.0001 | 0.2002 | ns |
N:Mg | Cúbica | <.0001 | <.0001 | 0.0024 |
Específicamente para esta investigación, durante el ciclo de un año, el hecho de que la tendencia de nutrientes o sus proporciones sea cuadrática indica que hay un punto bajo para una estación determinada, mientras que cuando es cúbica se presenta un punto máximo y un mínimo estacional; es decir, dos cambios de pendiente en la tendencia.
La Figura 2 muestra que la tendencia de caída de acícula fue mayor durante el verano y menor en invierno. Las tendencias por nutrientes y proporciones con respecto al N se muestran es las Figuras 3 y 4, respectivamente. La concentración de N fue mayor a mediados del otoño (octubre 17), mientras que para P y Ca fue mayor en el verano (junio 20 a julio 7). En el caso del Mg, las concentraciones fueron mayores en primavera (marzo 16 a mayo 5) (Figura 3). Excepto por la relación N:Ca, la mayoría de la proporciones de nutrientes (N:P, N:K, N:Mg) presentaron máximos al final del verano (agosto 28) y otoño (octubre 17) (Figura 4).
![](/img/revistas/rcscfa/v28n1//2007-4018-rcscfa-28-01-57-gf6.png)
Figura 2 Tendencias temporales en la masa de caída de acículas de Pinus hartwegii con valores estimados e intervalos de confianza del 95 % en dos bosques alpinos del centro de México. Letras minúsculas y mayúsculas indican meses del 2017 y 2018, respectivamente.
![](/img/revistas/rcscfa/v28n1//2007-4018-rcscfa-28-01-57-gf7.png)
Figura 3 Tendencias en la concentración de nutrientes de caída de acículas de Pinus hartwegii con intervalos de confianza del 95 % en dos bosques alpinos del centro de México. Letras minúsculas y mayúsculas indican meses del 2017 y 2018, respectivamente.
Discusión
Concentraciones en piso forestal y suelo
Comparado con JO, el piso forestal en TL mostró mayor concentración de Mg (P < 0.001) y Ca (P < 0.009). Tal diferencia también se reflejó en las concentraciones del suelo mineral, lo que corrobora la estrecha relación de la fertilidad del suelo y la química de las acículas (Wood et al., 2006). Una suma mayor de cationes (K + Ca + Mg) en TL indica mejor fertilidad y refuerza la relación entre la cantidad de nutrientes potencialmente disponibles y la composición química de las acículas, ya que el suelo es de origen volcánico (Andisol) en ambas montañas. En estas, la baja concentración de los principales cationes (<6 cmolc∙kg-1) es la esperada en los suelos forestales de bosques de P. hartwegii en el centro de México (Marín, Escolero-Fuentes, & Trinidad-Santos, 2002).
Masa de acículas
El hecho de que JO haya presentado mayor cantidad de masa de mantillo se tradujo también en reservorios más grandes de todos los nutrientes (P < 0.001); sin embargo, llama la atención que aun con mayores reservorios en JO, el bosque TL tenga mayor productividad. Los árboles en JO tienen incremento medio anual en área basal de 6 cm2, mientras que en TL es de 15 cm2 (Correa‐Díaz et al., 2020; Núñez-García et al., 2020).
Los resultados indican que, en JO, una proporción alta de la productividad neta se destina a la producción de masa foliar. Dependiendo de la estructura y edad del bosque, la producción de acículas puede representar más del 70 % de los detritos del dosel (Schlesinger & Bernhardt, 2013). Anualmente, fluyen mayores cantidades de nutrientes del dosel al piso forestal en JO, pero no parece haber prioridad en la productividad de componentes maderables. De hecho, la sincronía entre la actividad fotosintética a nivel dosel (a través de índices de vegetación- NDVI) y la formación de anillos de crecimiento anual es mayor en TL que en JO (Correa-Díaz, Gómez-Guerrero, Vargas-Hernández, Rozenberg, & Horwath, 2020).
Otros autores tampoco han encontrado correlación entre el crecimiento radial y la producción de acículas en bosques de Pinus silvestris L. (Lehtonen, Lindholm, Hokkanen, Salminen, & Jalkanen, 2008), pero sí han encontrado pulsos de producción de acícula cada cuatro años, posiblemente asociados a eventos de circulación climática global (González de Andrés et al., 2019). En cuanto a mayores reservorios del piso forestal, lo que se sabe es que la acumulación de mantillo está más bien correlacionada con los procesos de descomposición que con la productividad de un bosque (Frank, Pontes, & McFarlane, 2012). Los resultados de este estudio sugieren que podría haber otros factores de la calidad de acícula y descomposición, como la proporción N:lignina, que podrían explicar este resultado (Schlesinger & Bernhardt, 2013).
Es probable que en el sitio JO, por ser un volcán aislado y no parte de una cadena montañosa como TL (Figura 1), ocurran eventos más extremos en humedad y temperatura y mayor velocidad de viento que provoquen mayor caída de acícula, menor descomposición y menor crecimiento de los árboles (Zhu et al., 2021). Los datos de los dataloggers mostraron que la temperatura es ligeramente menor y la humedad relativa es mayor en JO que en TL (6.83 vs. 7.25 °C; 73.13 vs. 69.74 %).
Flujo anual de nutrientes
En ambas montañas hay tendencias definidas de la masa de acícula y nutrientes entre junio del año previo (2017) y junio del año posterior (2018), con los valores más bajos a principios de invierno. El aporte total de acículas de pino fue 11.2 Mg∙ha-1∙año-1 en JO, mientras que en TL fue 6.7 kg∙ha-1∙año-1 (Cuadro 5). Estos datos son altos comparados con los reportados en otros ecosistemas; por ejemplo, 7.6 Mg∙ha-1∙año-1 para hojarasca total (hojas, ramillas y estructuras reproductivas) en bosques de pino-encino de Chiapas (Rocha-Loredo & Ramírez-Marcial, 2009); 5.9 Mg∙ha-1∙año-1 en bosques de pino-encino de Durango (Moreno Valdez et al., 2018) y 7.4 y 10.6 Mg∙ha-1∙año-1 en bosques de encino de Nuevo León y Morelos, respectivamente (González-Rodríguez et al., 2011; Mantero-García, Gómez-Guerrero, Gavi-Reyes, Zamora-Morales, & Ramírez-Ayala, 2019). De acuerdo con una revisión global, Zhang, Yuan, Dong, y Liu (2014) determinaron que la hojarasca total en ecosistemas forestales varía de 3 a 11 Mg∙ha-1∙año-1, aunque los extremos superiores se extienden hasta los 13 Mg∙ha-1∙año-1.
En años húmedos y con vientos fuertes, la caída de hoja puede ser hasta más de dos veces mayor en comparación con años de clima promedio, lo cual afecta particularmente a los bosques de coníferas con mayor desprendimiento de acículas (Cheng et al., 2020). Por falta de estaciones meteorológicas, no fue posible saber si en el año de estudio haya habido mayor caída de hoja por factores abióticos. Los valores de caída de hoja representan sitios con buena cobertura arbórea (>60 %), pero en sitos de menor cobertura los valores serán menores. Aunque todos los sitios presentan cobertura del 100 % en sotobosque, es probable también que por la ubicación de la trampa a nivel de suelo se puedan captar algunas acículas por efecto de la gravedad.
En ambas montañas, la máxima acumulación de acículas ocurrió en los meses con mayor precipitación (Figura 2), cuando la actividad de crecimiento es mayor, los árboles reemplazan el follaje viejo y los nutrientes móviles como el N y P se transportan de tejido viejo a nuevo (Schlesinger & Bernhardt, 2013). Estos resultados contrastan con los de otros estudios, donde la mayor producción de hojarasca se asocia a la temporada más seca del año (López-Escobar et al., 2018); sin embargo, una tendencia bimodal, incluyendo temporada seca y húmeda se reportó para bosques de encino-pino (Quercus potosina Trel. y Pinus cembroides Zucc.) en San Luis Potosí, con producción total de hojarasca de 5.9 Mg∙ha-1∙año-1 (Pérez-Suárez et al., 2009). Los resultados del presente estudio indican un solo pico de producción de acícula que ocurre en verano (Figura 2). Así, a pesar de desarrollarse sobre climas fríos, los resultados evidencian que los bosques estudiados son muy productivos en acícula y dinámicos en la transferencia de nutrientes.
Eficiencia de uso de nutrientes
Utilizando la información del Cuadro 5, la eficiencia de uso de nutrientes en caída de acícula se puede estimar a nivel de ecosistema. Dicha estimación es el cociente del flujo anual de masa de acícula con respecto al flujo anual de cada nutriente (kg∙ha-1∙año-1). Las eficiencias en JO son 114, 2 153, 1 287, 485, 2 015 para N, P, K, Ca y Mg, respectivamente, en tanto que para TL son 121, 1 971, 859, 465, 1 425 en el mismo orden de nutrientes. Estos resultados indican que la eficiencia de uso de N es más eficiente en la montaña TL, ya que una misma unidad de nutriente (kg) se tradujo en cantidades mayores de masa foliar (kg). Pero en términos de K y Mg hubo mejor eficiencia en JO, lo que particularmente es consistente con una mayor suma de cationes en el suelo de TL (Cuadro 3). Independientemente de la localidad, estos niveles de eficiencia son más altos para las bases K, Ca y Mg que los encontrados en un rodal de Pinus pseudostrobus Lindl. mezclado con especies caducifolias (González-Rodríguez et al., 2011) y, excepto por la eficiencia en P, también superan la eficiencia de nutrientes del bosque tropical lluvioso (Kotowska, Leuschner, Triadiati, & Hertel, 2016).
Tendencias en el tiempo
En ambas montañas, la producción de acícula se reduce en los meses fríos, reflejando la baja actividad fisiológica en los árboles. Para bosques alpinos este es un resultado esperado, debido a que la temperatura limita un gran número de procesos ecofisiológicos (Körner, 2012). Las concentraciones de los nutrientes en la caída de acícula son resultado de las concentraciones en acícula verde y la reabsorción (Kobe, Lepczyk, & Iyer, 2005); cuanto más baja es la concentración de nutrientes en caída de acículas, mayor es la cantidad de reabsorción al interior de los árboles. Ambas montañas tienen un patrón similar con concentraciones más bajas (mayor reabsorción) en los meses de marzo a mayo (Figura 2); es decir, en los meses secos cuando la provisión de nutrientes vía suelo es baja por la menor humedad del suelo (Gómez-Guerrero & Doane, 2018).
El P mostró un aumento gradual de reutilización de julio a enero, esto es, mayor demanda en la medida que el periodo lluvias se iba alejando (Figura 3). Lo anterior podría deberse a que P es un elemento de baja movilidad que se absorbe principalmente por difusión de masas, requiriendo humedad del suelo para favorecer la presencia de sus formas solubles (Binkley & Fisher, 2020). En ambas montañas, la similitud de la proporción N:P indicó una reutilización mayor en los meses de junio y julio; posteriormente, esta reutilización disminuyó en invierno y se reactivó a finales de la primavera del siguiente año, que es una tendencia acoplada con el ritmo de crecimiento arbóreo y que denota la importancia de estos dos macroelementos (N y P) para la formación de biomasa (Schlesinger & Bernhardt, 2013).
El Ca indicó una demanda progresiva de julio a marzo, lo que refleja el proceso de formación de pared celular en las acículas; es decir, a diferencia de otros elementos, se infiere un proceso acumulativo en la medida que se concluye la expansión foliar. La variación de la proporción N:Ca indica que el mayor consumo de Ca ocurre en la época de lluvias, cuando las acículas empiezan su expansión. El hecho de que dicha relación se mantenga estable el resto del año indica que el Ca tiene baja movilidad y que se va acumulando en las estructuras de soporte como la pared celular.
Las concentraciones más bajas de Mg en JO concuerdan con la menor disponibilidad en el suelo (Cuadro 3). Los resultados indican que, en JO, las acículas que caen al piso ya van con bajas concentraciones de Mg, debido a su baja disponibilidad en el suelo, lo que condiciona a los árboles incrementar la reutilización de Mg. En cambio en TL, el suelo tiene mejor abasto de Mg, de tal forma que las acículas regresan al suelo con concentraciones más altas de Mg (Figura 3). Dado que la productividad de TL es alta comparada con JO, los resultados denotan provisión suficiente de Mg y menor necesidad de reutilización; además, la proporción N:Mg es más estable en TL a lo largo del año. La proporción más baja de N:Mg en TL indica un mejor balance de nutrientes, sustentado por la mayor fertilidad de suelo (Cuadro 2). Aunque en ambos lugares la proporción N:Mg varía durante el año, es evidente que los cambios son mayores en JO, donde las acículas caídas de agosto a octubre se empobrecen considerablemente en Mg (mayor relación N:Mg), denotando que en JO se requieren más ajustes en la proporción de nutrientes para mantener el crecimiento arbóreo.
Llama la atención que el K por sí solo no haya mostrado tendencia temporal, pero cuando se considera la proporción N:K, es notable la importancia del K y su mayor amplitud de respuesta. La proporción N:K se manifestó en dos fases; la primera de agosto a diciembre (verano-otoño) cuando proporcionalmente el K es bajo respecto a N, mientras que la segunda fase ocurre en los meses secos, de enero a mayo (invierno-primavera), cuando la proporción K respecto a N aumenta (disminuye N:K). Dado que el K está ligado a la regulación estomática, transporte de floema y fotosíntesis, este resultado podría indicar un aumento en la reabsorción de K en otoño, previo a la estación fría, mientras que en primavera la reabsorción de K disminuye, cuando la nueva estación de crecimiento y fotosíntesis está por iniciar. Por otra parte, previo al invierno, las hojas ya han acumulado celulosa y el N se diluye, resultando en proporciones N:K más bajas antes de empezar la nueva estación de crecimiento (Schlesinger & Bernhardt, 2013). A nivel de ecosistema, el papel del K es complejo, pero una revisión amplia del tema indica una fuerte relación con el N (Tripler, Kaushal, Likens, & Walter, 2006), lo que coincide completamente con este estudio. Al respecto, se ha propuesto que los pulsos paralelos de K y nitratos en los arroyos de las cuencas indican la fuerte relación de estos elementos (Tripler et al., 2006). Aunque la temperatura media de ambos sitios es similar, los días con temperaturas bajo cero son mayores en JO (Núñez-García et al., 2020). De hecho, hay un pulso más marcado de reabsorción de K, indicado por la proporción N:K en JO, lo que coincide con la mayor frecuencia de heladas en este sitio (Núñez-García et al., 2020).
Conclusiones
Los resultados de este estudio son una línea base útil para detectar cambios en la funcionalidad futura de los bosques alpinos. Los bosques de los cerros Jocotitlan y Tláloc difieren en la producción de acícula y flujo de nutrientes. El bosque Jocotitlán, menos productivo en componentes maderables (v. g. área basal), tuvo mayor caída de acículas. La relación N:K en caída de acícula es más relevante de lo que se pensaba, ya que mostró variación estacional más marcada, resaltando la importancia del balance de estos macronutrientes en el follaje. El cambio climático podría conducir al estrés nutricional e hídrico de los bosques, por lo tanto, las investigaciones sobre la producción de caída de hoja y su composición química pueden ser indicadores importantes para valorar los efectos en la vegetación y proponer estrategias de conservación. La estequiometría de la caída de follaje da información crucial para comprender la biogeoquímica de los bosques alpinos.
Agradecimientos
Al Consejo Nacional de Ciencia y Tecnología (CONACyT) por la beca de posgrado asignada a la primera autora. A la Unión de Ejidos de la Montaña y autoridades municipales de Jocotitlán del Estado de México, por el permiso de acceso a sus áreas forestales.
REFERENCIAS
Acosta Mireles, M., Carrillo Anzures, F., Delgado, D., & Velasco Bautista, E. (2014). Establecimiento de parcelas permanentes para evaluar impactos del cambio climático en el Parque Nacional Izta-Popo. Revista Mexicana de Ciencias Forestales, 5(26), 06‒29. doi: 10.29298/rmcf.v5i26.287 [ Links ]
Alcántara, G. G., & Sandoval, V. M. (1999). Manual de análisis químico de tejido vegetal. Guía de muestreo, preparación, análisis e interpretación. Chapingo, Estado de México, México: Sociedad Mexicana de la Ciencia del Suelo A. C. [ Links ]
Arce, J. L., Layer, P. W., Macías, J. L., Morales-Casique, E., García-Palomo, A., Jiménez-Domínguez, F. J., …Vásquez-Serrano, A. (2019). Geology and stratigraphy of the Mexico basin (Mexico city), central Trans-Mexican volcanic Belt. Journal of Maps, 15(2), 320‒332. doi: 10.1080/17445647.2019.1593251 [ Links ]
Binkley, D., & Fisher, R. (2020). Ecology and management of forest soils (5th ed.). UK: Wiley-Blackwell. [ Links ]
Bray, R. H., & Kurtz, L. (1945). Determination of total, organic, and available forms of phosphorus in soils. Soil Science, 59(1), 39‒46. doi: 10.1097/00010694-194501000-00006 [ Links ]
Bremmer, J. M. (1965). Total nitrogen. In C. A. Black (Ed.), Methods of soil analysis. Part 2: Chemical and microbial properties (1049‒1178). Madison, WI, USA: American Society of Agronomy. [ Links ]
Cheng, C.-H., Lee, C.-Y., Lee, H.-R., Chen, C.-P., & Menyailo, O. V. (2020). Effects of typhoon disturbances on seasonal and interannual patterns of litterfall on coniferous and broadleaf plantations in Xitou, central Taiwan. Journal of Forest Research, 25(3), 155‒162. doi: 10.1080/13416979.2020.1762026 [ Links ]
Correa-Díaz, A., Gómez-Guerrero, A., Vargas-Hernández, J. J., Rozenberg, P., & Horwath, W. (2020). Long-term wood micro-density variation in alpine forests at Central México and their spatial links with remotely sensed information. Forests, 11(4), 452. doi: 10.3390/f11040452 [ Links ]
Correa‐Díaz, A., Silva, L., Horwath, W., Gómez‐Guerrero, A., Vargas‐Hernández, J., Villanueva‐Díaz, J., …Velázquez‐Martínez, A. (2020). From trees to ecosystems: Spatiotemporal scaling of climatic impacts on montane landscapes using dendrochronological, isotopic, and remotely sensed data. Global Biogeochemical Cycles, 34(3), e2019GB006325. doi: 10.1029/2019GB006325 [ Links ]
Fitzmaurice, G. M., Laird, N. M., & Ware, J. H. (2004). Applied longitudinal analysis. USA Wiley. [ Links ]
Frank, D. A., Pontes, A. W., & McFarlane, K. J. (2012). Controls on soil organic carbon stocks and turnover among North American ecosystems. Ecosystems, 15(4), 604‒615. doi: 10.1007/s10021-012-9534-2 [ Links ]
García, E. (2004). Modificaciones al sistema de clasificación climática de Köppen (Para adaptarlo a las condiciones de la República Mexicana) (5.a ed.). México: Instituto de Geografía UNAM. [ Links ]
Gavinet, J., Ourcival, J.-M., & Limousin, J.-M. (2019). Rainfall exclusion and thinning can alter the relationships between forest functioning and drought. New Phytologist, 223(3), 1267‒1279. doi: 10.1111/nph.15860 [ Links ]
Ge, J., & Xie, Z. (2017). Leaf litter carbon, nitrogen, and phosphorus stoichiometric patterns as related to climatic factors and leaf habits across Chinese broad-leaved tree species. Plant Ecology, 218(9), 1063‒1076. doi: 10.1007/s11258-017-0752-8 [ Links ]
Gómez-Guerrero, A., & Doane, T. (2018). The response of forest ecosystems to climate change. In W. R. Horwath & K. Yakov (Eds.), Climate change impacts on soil processes and ecosystem properties (vol. 35, pp. 185‒206). Elsevier. doi: 10.1016/B978-0-444-63865-6.00007-7 [ Links ]
González-Rodríguez, H., Domínguez-Gómez, T. G., Cantú-Silva, I., Gómez-Meza, M. V., Ramírez-Lozano, R. G., Pando-Moreno, M., & Fernández, C. (2011). Litterfall deposition and leaf litter nutrient return in different locations at Northeastern Mexico. Plant Ecology, 212(10), 1747. doi: 10.1007/s11258-011-9952-9 [ Links ]
González-Rodríguez, H., Ramírez-Lozano, R. G., Cantú-Silva, I., Gómez-Meza, M. V., Estrada-Castillón, E., & Arévalo, J. R. (2018). Deposition of litter and nutrients in leaves and twigs in different plant communities of northeastern Mexico. Journal of Forestry Research, 29(5), 1307‒1314. doi: 10.1007/s11676-017-0553-x [ Links ]
González de Andrés, E. (2019). Interactions between climate and nutrient cycles on forest response to global change: The role of mixed forests. Forests, 10(8), 609. doi: 10.3390/f10080609 [ Links ]
González de Andrés, E., Blanco, J. A., Imbert, J. B., Guan, B. T., Lo, Y. H., & Castillo, F. J. (2019). ENSO and NAO affect long‐term leaf litter dynamics and stoichiometry of Scots pine and European beech mixedwoods. Global Change Biology, 25(9), 3070‒3090. doi: 10.1111/gcb.14672 [ Links ]
Kobe, R. K., Lepczyk, C. A., & Iyer, M. (2005). Resorption efficiency decreases with increasing green leaf nutrients in a global data set. Ecology, 86(10), 2780‒2792. doi: 10.1890/04-1830 [ Links ]
Körner, C. (2012). Alpine treelines: functional ecology of the global high elevation tree limits. Springer Science & Business Media. doi: 10.1007/978-3-0348-0396-0 [ Links ]
Kotowska, M. M., Leuschner, C., Triadiati, T., & Hertel, D. (2016). Conversion of tropical lowland forest reduces nutrient return through litterfall, and alters nutrient use efficiency and seasonality of net primary production. Oecologia, 180(2), 601‒618. doi: 10.1007/s00442-015-3481-5 [ Links ]
Lehtonen, A., Lindholm, M., Hokkanen, T., Salminen, H., & Jalkanen, R. (2008). Testing dependence between growth and needle litterfall in Scots pine-a case study in northern Finland. Tree Physiology, 28(11), 1741‒1749. doi: 10.1093/treephys/28.11.1741 [ Links ]
López-Escobar, N. F., Gómez-Guerrero, A., Velázquez-Martínez, A., Fierros-González, A. M., Castruita-Esparza, L. U., & Vera-Castillo, J. A. (2018). Reservorios y dinámica de nutrientes en dos rodales bajo aprovechamiento de Pinus montezumae Lamb. en Tlaxcala, México. Revista Chapingo Serie Ciencias Forestales y del Ambiente, 24(1), 115‒129. doi: 10.5154/r.rchscfa.2017.09.055 [ Links ]
Mantero-García, H. D., Gómez-Guerrero, A., Gavi-Reyes, F., Zamora-Morales, B. P., & Ramírez-Ayala, C. (2019). ¿Es sustentable el aprovechamiento de tierra de hoja en bosques de encino? Madera y Bosques, 25(3), e2531807. doi: 10.21829/myb.2019.2531807 [ Links ]
Marín, L. E., Escolero-Fuentes, O., & Trinidad-Santos, A. (2002). Physical geography, hydrogeology, and forest soils of the basin of Mexico. In M. E. Fenn, L. I. de Bauer, & T. Hernandez-Tejeda (Eds.), Urban Air Pollution and Forests (pp. 44‒67). USA: Springer. doi: 10.1007/978-0-387-22520-3 [ Links ]
Moreno Valdez, M. E., Domínguez Gómez, T. G., Alvarado, M. d. S., Colín, J. G., Corral Rivas, S., & González Rodríguez, H. (2018). Aporte y descomposición de hojarasca en bosques templados de la región de El Salto, Durango. Revista Mexicana de Ciencias Forestales, 9(47), 70‒93. doi: 10.29298/rmcf.v9i47.180 [ Links ]
Neumann, M., Ukonmaanaho, L., Johnson, J., Benham, S., Vesterdal, L., Novotný, R., …Michopoulos, P. (2018). Quantifying carbon and nutrient input from litterfall in European forests using field observations and modeling. Global Biogeochemical Cycles, 32(5), 784‒798. doi: 10.1029/2017GB005825 [ Links ]
Núñez-García, A., Gómez-Guerrero, A., Terrazas-Salgado, T. M., Vargas-Hernández, J. J., & Villanueva-Díaz, J. (2020). Analysis of basal area increment of Pinus hartwegii Lindl. at different elevations and aspects on Jocotitlán Mountain, State of Mexico. Revista Chapingo Serie Ciencias Forestales y del Ambiente, 27(1), 77‒88. doi: 10.5154/r.rchscfa.2019.10.074 [ Links ]
Pérez-Suárez, M., Arredondo-Moreno, J., Huber-Sannwald, E., & Vargas-Hernández, J. (2009). Production and quality of senesced and green litterfall in a pine-oak forest in central-northwest Mexico. Forest Ecology and Management, 258(7), 1307‒1315. doi: 10.1016/j.foreco.2009.06.031 [ Links ]
R Core Team. (2020). A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. Retrieved from https://www.R-project.org/ [ Links ]
Rocha-Loredo, A. G., & Ramírez-Marcial, N. (2009). Producción y descomposición de hojarasca en diferentes condiciones sucesionales del bosque de pino-encino en Chiapas, México. Boletín de la Sociedad Botánica de México, 84, 1‒12. doi: 10.17129/botsci.2287 [ Links ]
Sardans, J., Alonso, R., Janssens, I. A., Carnicer, J., Vereseglou, S., Rillig, M. C., …Penuelas, J. (2016). Foliar and soil concentrations and stoichiometry of nitrogen and phosphorous across European Pinus sylvestris forests: relationships with climate, N deposition and tree growth. Functional Ecology, 30(5), 676‒689. doi: 10.1111/1365-2435.12541 [ Links ]
SAS Institute Inc. (SAS). (2011). SAS/ETS 9.3 User’s Guide. Cary, NC, USA: Author. [ Links ]
Sayer, E. J., Rodtassana, C., Sheldrake, M., Bréchet, L. M., Ashford, O. S., Lopez-Sangil, L., …Wright, S. J. (2020). Revisiting nutrient cycling by litterfall-Insights from 15 years of litter manipulation in old-growth lowland tropical forest. In A. J. Dumbrell, E. C. Turner, & T. M. Fayle (Eds.), Advances in ecological research (vol. 62, pp. 173‒223): Elsevier. doi: 10.1016/bs.aecr.2020.01.002 [ Links ]
Schlesinger, W. H., & Bernhardt, E. S. (2013). Biogeochemistry: an analysis of global change (3rd ed.). San Diego, California, USA: Academic Press. [ Links ]
Silva, L. C., Gómez‐Guerrero, A., Doane, T. A., & Horwath, W. R. (2015). Isotopic and nutritional evidence for species‐and site‐specific responses to N deposition and elevated CO2 in temperate forests. Journal of Geophysical Research: Biogeosciences, 120(6), 1110‒1123. doi: 10.1002/2014JG002865 [ Links ]
Tripler, C. E., Kaushal, S. S., Likens, G. E., & Walter, M. T. (2006). Patterns in potassium dynamics in forest ecosystems. Ecology Letters, 9(4), 451‒466. doi: 10.1111/j.1461-0248.2006.00891.x [ Links ]
Wood, T. E., Lawrence, D., & Clark, D. A. (2006). Determinants of leaf litter nutrient cycling in a tropical rain forest: soil fertility versus topography. Ecosystems, 9(5), 700‒710. doi: 10.1007/s10021-005-0016-7 [ Links ]
Zar, J. H. (2010). Biostatistical analysis: Pearson new international edition (5th ed.). New Jersey, USA: Prentice Hall. [ Links ]
Zhang, H., Yuan, W., Dong, W., & Liu, S. (2014). Seasonal patterns of litterfall in forest ecosystem worldwide. Ecological Complexity, 20, 240‒247. doi: 10.1016/j.ecocom.2014.01.003 [ Links ]
Zhu, X., Zou, X., Lu, E., Deng, Y., Luo, Y., Chen, H., & Liu, W. (2021). Litterfall biomass and nutrient cycling in karst and nearby non-karst forests in tropical China: A 10-year comparison. Science of the Total Environment, 758, 143619. doi: 10.1016/j.scitotenv.2020.143619 [ Links ]
Recibido: 16 de Diciembre de 2020; Aprobado: 11 de Octubre de 2021