Introducción

La Orinoquia colombiana es considerada como una de las regiones biogeográficas con mayor extensión en el país (9 238 277 ha), la cual alberga una gran diversidad de fauna y flora que cumple un papel fundamental dentro de la interconexión entre ecosistemas. La vegetación de la subregión conocida como altillanura está compuesta principalmente por coberturas de sabanas y pastizales, bosques de galería, morichales y plantaciones forestales (^{Rippstein, Escobar, & Motta, 2001}), que proveen una serie de servicios ecosistémicos importantes para el ambiente y la sociedad. No obstante, en las últimas tres décadas en la región se ha acelerado el cambio de las coberturas naturales de sabanas por pasturas, cultivos agrícolas y forestales, ganadería extensiva, explotación de hidrocarburos, entre otras, principalmente en la subregión denominada altillanura bien drenada (^{Romero, et al., 2009}). Cerca del 14 % de las sabanas se han transformado durante 20 años con la forestación de dichas áreas para la producción de madera, incrementar la captura de carbono y mejorar las condiciones del suelo como un modelo agroindustrial de desarrollo (^{Peñuela & Rodríguez, 2014}). Particularmente en el municipio de Villanueva-Casanare, extensiones considerables de plantaciones son establecidas principalmente con especies forestales como Pinus caribaea y Eucalyptus pellita (^{Refocosta S.A.S, 2019}), que junto con la introducción de pasturas han reemplazado las coberturas de sabana intervenidas (^{Rippstein et al., 2001}; ^{Peñuela, Fernández, Castro, & Ocampo, 2011}).

La artropofauna es considerada como elemento fundamental para la provisión de servicios ecosistémicos en un área (^{Decaëns, Jiménez, Gioia, Measey, & Lavèlle, 2006}). La fauna epiedáfica es sensible a cambios provocados en el suelo por factores ambientales y antrópicos (^{Uribe et al., 2010}) por ello son considerados como bioindicadores de la calidad de los suelos, como es el caso de la clase Collembola (^{Culik & Zeppelini, 2003}). Los organismos responden de manera diferente ante los disturbios, y aunque se tienen indicios del efecto de algunos cultivos agrícolas, pasturas y perennes que sugieren, que la fauna de lombrices se ve favorecida en las pasturas, pero los demás grupos se ven disminuidos en usos intensivos como las pasturas introducidas y los cultivos, siendo más favorables las coberturas con perennes (^{Decaëns et al., 2006}; ^{Sanabria, Dubs, Lavèlle, Fonte, & Barot, 2016}). Particularmente, las comunidades de carábidos, en la cobertura Eucaliptus pellita en la misma área de este estudio, se encuentran dominadas por especies generalistas cuya diversidad y abundancia (q0 y q1) disminuyen en plantaciones de mayor desarrollo (^{García-Suabita, Pinzón, Spence & Pinzón- Florián,2019}). La diversidad de Collembola depende de factores como la disponibilidad de materia orgánica y la humedad presente en el suelo (^{Palacios-Vargas, 2000}; ^{Socarrás, 2013}). Los colémbolos están relacionados con la disponibilidad de nutrientes para las plantas, el almacenamiento y regulación hídrica, la estabilidad de los suelos (mejorando su estructura) y el control de la humedad y pH necesarios para la fertilidad (^{Lavèlle et al., 2014}).

Los insectos de suelos en la región de la Orinoquia ha sido poco estudiados (^{Morales-Castaño y Medina, 2009}) pero particularmente se desconocen efectos de diferentes usos del suelo sobre la microfauna. De esta manera, se considera que el estudio de la colembofauna presente en el uso anterior y actual del suelo puede contribuir a identificar factores que afectan la diversidad de fauna epiedáfica como lo son la introducción de plantaciones comerciales y su relación directa con el uso del suelo, para de esta manera establecer estrategias integrales de manejo forestal sostenible.

Materiales y métodos

Sitio de estudio: El estudio se desarrolló en el municipio de Villanueva en el departamento del Casanare, localizado a 550 m.s.n.m. en las coordenadas (4°39’10.73” N & 72°54’58.76” W). El clima está clasificado como tropical monzónico (Am) (^{García, 2004}); presenta unrégimen deprecipitación monomodal siendo los meses más lluviosos de abril a noviembre y el periodo de sequía entre diciembre y marzo. En promedio la precipitación alcanza 2 976 mm al año, la temperatura media anual es de 26.1 °C y la humedad relativa es del 76 %. La zona de estudio está ubicada en la zona de vida bosque húmedo tropical (bh-T) y presenta suelos de tipo arenoso a franco-arenoso, con fertilidad baja, alta concentración de aluminio y hierro, pobres en fósforo, nitrógeno y materia orgánica (^{Refocosta S.A.S, 2019}).

Muestreo de colémbolos: Con el objeto de comparar la diversidad de colémbolos se muestrearon tres lotes de cada edad de plantaciones jóvenes y adultas de E. pellita (2, 6, 8 y 14 años) con poco manejo (labores silviculturales y culturales), tres bosques de galería y dos lotes con pasturas, para un total de seis coberturas comparadas en 17 unidades de muestreo. El estudio se ejecutó durante los meses de noviembre 2014 (época seca) y mayo 2015 (época húmeda). No se reemplazaron bosques naturales para establecer las plantaciones, la inclusión de la diversidad de los bosques de galería es solamente de referencia. En cada uno de los lotes se delimitaron 3 transectos en banda de 50 m de largo por 3 m de ancho con una distancia de separación de 20 m entre cada transecto. Dentro de cada transecto se delimitaron 5 sub-transectos de 10 m de largo por 3 m de ancho, correspondiente a la distancia de siembra de la plantación y en cada uno se instaló una trampa de caída (^{Moreira, Huising, & Bignell, 2012}). Cada trampa consistió en un recipiente plástico con capacidad de 500 cm³ enterrado al límite de la superficie del suelo y protegido de la lluvia en el cual se utilizaron 200 ml de líquido preservante (consistente de una solución de Propilenglicol al 40 %, ^{Spence & Niemelä, 1994}). Se instalaron un total de 255 trampas de caída (por época de muestreo) y se analizaron las 3 trampas centrales de cada transecto, es decir se evaluaron únicamente un total de 153 trampas por cada época.

Con el objeto de explicar la composición de la colembofauna en función de variables ambientales, se incluyeron variables físico-químicas del suelo y variables dasométricas de las plantaciones. A cada uno de los árboles de la línea central del transecto se le determinó el diámetro a la altura del pecho (DAP) con una cinta diamétrica y la altura total con un hipsómetro Vertex IV (Haglof®, Langsele, Sweden). Se estimó el porcentaje de cobertura de dosel promedio realizando cuatro medidas sobre cada una de las trampas instaladas con un densitómetro convexo (Forestry Suppliers, Inc., Jackson, MS, USA). Además, se midió la altura de la acumulación de materia orgánica enterrando una escala métrica hasta alcanzar la superficie del horizonte A del suelo en tres puntos a cada lado de cada trampa. También se estimó el porcentaje de área sin cobertura (vegetal o por acumulación de materia orgánica) en el cuadro circundante de cada trampa, equivalente a los cuatro árboles más cercanos sembrados a distancia de 3 x 3 m (correspondiente a la distancia de siembra). Por otra parte, en campo se determinó la humedad, conductividad volumétrica y la temperatura del suelo al lado de cada trampa de caída utilizando un reflectómetro (TDR; ICT®) y termómetros de suelo (Amprobe®, Everett, WA, USA). Asimismo, se recolectaron 784 cm³ de suelo del centro de cada transecto referenciado y posteriormente, se realizó el análisis físico-químico de la siguiente manera: el porcentaje contenido de humedad se determinó mediante el cálculo de diferencia de pesos húmedo y seco de las muestras, el pH se determinó con ayuda de un potenciómetro, el potasio y el calcio mediante el método de extracción usando NH4-Acetato 1 M pH 7 y posteriormente analizado con absorción atómica, el fósforo disponible se determinó mediante el método Bray II por colorimetría. El porcentaje de carbono orgánico se procesó mediante el método Walkley Black por colorimetría, el porcentaje de nitrógeno mediante el método volumétrico Kieldahl, y finalmente la textura se determinó mediante el método de Bouyoucos (valoración densimétrica por dispersión).

Selección, identificacióny análisis de muestras: Para la selección de los colémbolos recolectados se tomó una sub-muestra del 50 % del total de los individuos mismos que se identificaron con ayuda de un estereoscopio. Previamente, se realizó el proceso de aclarado y montaje según la metodología propuesta por ^{Palacios-Vargas y Mejía-Recamier (2007}). Adicionalmente, se llevó a cabo la clasificación morfológica con ayuda de un microscopio de contraste de fases y posteriormente con ayuda de textos de referencia y claves de identificación de géneros (^{Palacios-Vargas, 1990}; ^{Díaz Aspiazu, González Cairo, Palacios-Vargas, & Luciañez-Sánchez, 2004}; ^{Ospina, Serna, Peñaranda, & Serna, 2003}; ^{Palacios-Vargas & Peñaranda-Parada, 2005}; ^{Ávila & Jaramillo, 2009}; ^{Ospina & Palacios-Vargas., 2009}; ^{Bellinger, Christiansen, & Janssens, 1996}).

Análisis de datos: La riqueza de las especies presentes se determinó utilizando la metodología de rarefacción basada en cobertura para los órdenes de diversidad q0, q1 y q2 (^{Chao et al., 2014}), los cuales se relacionan directamente con la riqueza (se analiza bajo la importancia de especies poco comunes), el índice de Shannon (donde se relaciona el peso equivalente entre las especies teniendo en cuenta su abundancia relativa) y el índice de Simpson (mayor importancia a especies comunes). Teniendo en cuenta esto, los valores obtenidos en cada orden de diversidad representan el número de especies (o géneros) equivalentes los que se pueden comparar directamente entre sí (^{Jost, 2006}, ^Jost,2007). La metodología de rarefacción establece valores de riqueza comparables donde se determinan a partir del nivel de completitud del muestreo (el cual se analizó mediante curvas de acumulación de especies) y adicionalmente, permite establecer intervalos de confianza del 95 % para determinar si existen diferencias significativas entre los valores calculados. El análisis de los datos recolectados se realizó empleando el software R versión 3.2.5 mediante los paquetes estadísticos Vegan para el análisis de diversidad involucrando variables ambientales (^{Oksanen et al., 2016}). iNEXT para el análisis de riqueza (^{Hsieh & Chao, 2016}). Vegetarian para determinar la diversidad Alfa y Beta (^{Charney & Record, 2012}). Labdsv para la ordenación y clasificación de los conjuntos de datos de diversidad identificados (^{Roberts, 2016}) y también Indicspecies para determinar la relación entre presencia y abundancia de géneros de colémbolos encontrados (^{De Caceres & Legendre, 2009}).

Para establecer diferencias entre los promedios de las variables dasométricas y ambientales evaluadas en cada uno de los sitios y los obtenidos en la plantación se utilizó la prueba no paramétrica Kruskall-Wallis (^{García-Suabita et al., 2019}). Adicionalmente, se analizó el ensamblaje de colémbolos en función de las variables ambientales, donde se implementó un Análisis de Redundancia (RDA) el cual ajusta los datos en un modelo de regresión multivariada mediante técnicas de ordenación. Para el estudio se excluyeron aquellos géneros cuyas muestras tuvieran menos de 10 individuos, y previamente al análisis se eliminaron aquellas variables que evidenciaran una alta correlación (arena, limo, porcentaje de cobertura del dosel, contenido de humedad del suelo, DAP y CICE). Las variables restantes se analizaron mediante el factor de inflación de varianza (VIF) y se seleccionaron aquellas con VIF < 10 (^{Hair, Anderson, Tatham & Black, 1995}). Los análisis se realizaron por transecto ya que las muestras de suelo fueron representativas de cada uno de éstos. Igualmente, se analizó la diversidad Beta entre sitios mediante el cálculo de similitud de distancias de Bray-Curtis (^{Moreno, 2001}). La existencia de géneros indicadores se analizó según ^{Dufrêne y Legendre (1997}) y ^{De Caceres y Legendre (2009}), que tiene en cuenta la probabilidad de encontrar un género en una cobertura determinada y de igual manera en un hábitat específico (especificidad y fidelidad).

Resultados

En las áreas muestreadas en este estudio, incluyendo el uso anterior (pasturas) y el actual con plantaciones de E. pellita, se recolectaron 1 157 colémbolos, pertenecientes al orden Entomobryomorpha (61.7 %), Poduromorpha (17 %) y Symphypleona (21.2 %); 12 familias y 26 géneros en total (Apéndice 1). La mayor abundancia absoluta se encontró en la plantación de 8 años (22.2 %) y en la cobertura de pastos (20.1 %). Como era de esperarse, se observó mayor abundancia (57 %) en la época húmeda comparado con la época seca (42.9 %). Por otra parte, los muestreos representaron apropiadamente los géneros esperados en el área de estudio (Apéndice 2 y Apéndice 3).

Riqueza de colémbolos epiedáficos: Durante la época húmeda, la riqueza de colémbolos expresada en géneros poco comunes (q0) y la abundancia (q1) fue significativamente mayor en los pastos y en las plantaciones de 2 años en comparación con las otras edades de plantación (Fig. 1A). La época seca afectó la riqueza y abundancia de colémbolos en las pasturas en comparación con la misma cobertura en la época húmeda (de 26 géneros encontrados en la época húmeda se redujo a 12 en la época seca, en un 53.8 %) y en las plantaciones hasta 6 años (en las plantaciones de 2 y 6 años la reducción fue de 64.7 % y 25 % respectivamente). Mientras que en las plantaciones de 8 años en adelante no se observó una fluctuación marcada de número de especies o de abundancia, y contrario a lo esperado, en la época seca se observó mayor número efectivo de géneros (Fig. 1B).

Por otro lado, cuando se consideró la riqueza equivalente de géneros del orden q2 (con más peso los géneros comunes), en las plantaciones de 2 y 6 años disminuyó sustancialmente el número de géneros identificados para la época seca en comparación con la misma edad para la época húmeda (60 y 37.5 % respectivamente, Fig. 1C). Durante la época húmeda la mayor abundancia se observó en las plantaciones de 8 y 14 años (137 y 144 respectivamente), mientras que durante la época seca se observó en los bosques de galería.

Fig. 1 Fig. 1 Número efectivo de géneros de colémbolosepiedáficos por tipo de cobertura y época de muestreo en Villanueva-Casanare. A. Orden de riqueza q0. B. Orden de riqueza q1. C. Orden de riqueza q2. Effective number of epiedaphic Collembola genera by type of coverage and sampling season in Villanueva-Casanare. A. Richness order q0. B. Richness order q1. C. Richness order q2. 

Géneros indicadores: El análisis de especificidad y fidelidad permitió establecer los géneros Rastriopes (Börner, 1906), Folsomides (Stach, 1922)y Salina (MacGillivray, 1894) como afines a las coberturas de plantaciones de E. pellita de 8 y 14 años mientras que no se encontraron géneros afines a las edades de 2 y 6 años. Asimismo, Seira (Lubbock, 1870), Entomobrya (Rondani, 1861), Dicyrtomina (Börner, 1903) y Brachystomella (Ågren, 1904) son afines a las pasturas y Ptenothrix, Sphyrotheca (Börner, 1906) e Isotoma (Bourlet, 1839) fueron afines a la cobertura de bosque de galería (Tabla 1).

Similitud de géneros de colémbolos encontrados en las coberturas analizadas: La similitud de colembofauna epiedáfica de las pasturas estuvo representada parcialmente en las plantaciones. La similitud varió entre 41 y 66 % en la época húmeda y entre 51 y 67 % en la época seca. La mayor similaridad se observó en las plantaciones de mayor edad, E. pellita de 14 años (66 %) y E. pellita de 8 años (67 %). Entre edades de plantación también se observó variación en la similitud, los mayores valores (79 %) se observaron entre las plantaciones de E. pellita de 8 y 14 años de edad para la época húmeda, y 74 % entre E. pellita de 6 y 14 años de edad en la época seca. La similitud de géneros de colémbolos entre las plantaciones de E. pellita y el bosque de galería varió entre un 41 y 49 % para la época húmeda, mientras que para la época seca entre un 53 y 64 % (Tabla 2).

Composición de colémbolos en función de las variables de la estructura del bosque y las propiedades del suelo: El modelo de redundancia completo (RDA), explicó un 45 % de la varianza observada en la abundancia de colémbolos (R² = 0.26, F = 5.31, P = 0.001). De las 14 variables ambientales analizadas (porcentaje exposición del suelo alrededor de la trampa, acumulación de materia orgánica, altura de los árboles, temperatura del suelo, pH, humedad volumétrica, conductividad volumétrica, fósforo, carbono, nitrógeno, aluminio, calcio, potasio y textura), tres se relacionan con la diversidad de colémbolos observada: la exposición del suelo alrededor de la trampa (ESAT), la altura de la vegetación (ALTL) y el contenido de Nitrógeno (N) en el suelo (Apéndice 4). En el primer eje (RDA₁) se observan dos gradientes representados por el Nitrógeno (F = 8.67, P = 0.001), donde Ptenothrix está relacionado con altos contendidos de N en el suelo, lo cual diferencia al bosque de galería de otras coberturas, y para el caso de los pastos, la altura (F = 4.63, P = 0.001) se muestra como una variable inversa. Trogolaphysa y Ptenothrix son los de mayor dominancia para los bosques de galería, mientras Dicyrtomina y Brachystomella son los de mayor dominancia para los pastos. No obstante, en el segundo eje (RDA₂) se observa la representación del gradiente exposición del suelo alrededor de la trampa (F = 2.63, P = 0.010) (Fig. 2) relacionados con Lepidocyrtus y Collophora.

Fig. 2 Fig. 2 Diagrama RDA que muestra la composición del ensamblaje de colémbolos en función de las variables de la estructura del bosque y las propiedades del suelo. A. Diagrama de ordenación que muestra la distribución de los géneros de colémbolos en relación con las variables explicativas. B. El mismo modelo con intervalos de confianza del 95 % alrededor de los centroides del grupo superpuestos en el diagrama de ordenación que representa los diferentes hábitats (las elipses sólo se pueden usar como referencia para la interpretación, ya que el tipo de hábitat no se incluyó en el modelo). ESAT: Exposición del suelo alrededor de la trampa, ALTL: Altura del árbol, N: Nitrógeno. 1: Seira, 2: Entomobrya, 3: Cyphoderus, 4: Arlesia, 5: Sphaeridia, 6: Salina, 7: Collophora, 8: Lepidocyrtus, 9: Troglolaphysa, 10: Ptenothrix, 11: Dicyrtomina, 12: Brachystomella. RDA ordinations showing assemblage composition of Collembola as a function of forest structure variables and soil properties A. Ordination plot showing distribution of Collembola genera in relation to explanatory variables. B. Same model with 95 % confidence intervals around group centroids overlaid onto the ordination plot representing the different habitats (ellipses may be used only as a reference for interpretation as habitat type was not included in the model). ESAT: Exposure of the soil around the trap, ALTL: Height of trees, N: Nitrogen. 1: Seira, 2: Entomobrya, 3: Cyphoderus, 4: Arlesia, 5: Sphaeridia, 6: Salina, 7: Collophora, 8: Lepidocyrtus, 9: Troglolaphysa, 10: Ptenothrix, 11: Dicyrtomina, 12: Brachystomella. 

Discusión

Este estudio reporta la comparación de la colembofauna observada en cultivos con especies exóticas establecidos en áreas de la altillanura colombiana: pasturas y E. pellita. En general, el número de géneros de colémbolos epiedáficos observado en las plantaciones de E. pellita fue menor y el ensamblaje parcialmente similar al encontrado en las pasturas; sin embargo, sólo se observa menor abundancia de colémbolos en las plantaciones más jóvenes durante la época seca. Por otra parte, se observa un efecto importante de la edad de la plantación y de la época seca. Los ensamblajes de colémbolos encontrados estuvieron relacionados con el contenido de Nitrógeno en el suelo y por el desarrollo de cobertura arbórea (altura de los árboles). Aun cuando los cultivos comparados no reemplazaron los bosques de galería, la colembofauna observada en estas áreas se incluyó para tratar de entender el origen de la composición de la colembofauna encontrada en las plantaciones.

En las pasturas y plantaciones más jóvenes se observó un mayor número de géneros que en las plantaciones de edad media y adulta; sin embargo, al tomar en cuenta la abundancia general y los géneros más comunes (q1 y q2) las diferencias se disiparon, especialmente durante la época húmeda.La mayor diversidad de géneros observados en las pasturas puede obedecer a que en éstas se conserven y adapten mejor las especies adaptadas a áreas abiertas propias de las sabanas (uso anterior a los pastos). Teniendo en cuenta que las plantaciones más jóvenes tienen baja cobertura de dosel, la diversidad de áreas abiertas se ve favorecida por un tiempo dentro de la plantación. De manera contraria, en las plantaciones de desarrollo intermedio, escasea la vegetación herbácea y predomina la acumulación de hojarasca del eucalipto, aunque hacia el final de turno y debido a las entresacas, se pueden nuevamente desarrollar algunas herbáceas y arbustos de las familias Melastomataceae, Hypericaceae, Gentianaceae y Lauraceae en plantaciones de E. pellita, las cuales pueden favorecer el incremento de la diversidad de entomofauna al existir una mayor disponibilidad del recurso alimenticio (^{Fernández Méndez, Camargo Martínez, & Sarmiento, 2012}).

En este estudio, el régimen climático afectó de manera diferente a los colémbolos en las coberturas estudiadas. En general en la época seca, la región alcanza temperaturas extremas (alrededor de los 35 °C, ^{Weather Spark, 2016}) y baja humedad en el suelo (^{Camara, Correia, & Villela, 2012}).Durante la época secadisminuyó significativamente la diversidad y la abundancia de los colémbolos particularmente en las pasturas y las plantaciones de dos y seisaños (q0, q1 y q2), pero no en las plantaciones de mayor desarrollo. Lo anterior sugiere un efecto benéfico de la cobertura de dosel, posiblemente debido a la regulación de la temperatura del suelo y debido a que dichas coberturas presentan un conjunto heterogéneo de vegetación que satisface las necesidades de la fauna epiedáfica (^{CONABIO, 2011}), que también se observa en el bosque de galería. Para la época seca se registraron hallazgos de colémbolos que han desarrollado estrategias de adaptación para sobrevivir, como es el caso de Pseudachorutes, Brachystomella y Folsomides, y en general algunos pertenecientes al orden Entomobryomorpha y Poduromorpha (^{Arbea & Blasco-Zumeta, 2001}; ^{Guzmán, Melo, Lozano, & Rivera, 2010}).

En las pasturas, las plantaciones y los bosques de galería se observaron géneros indicadores, diferentes para cada tipo de cobertura. Los génerosSeira, Entomobrya, Dicyrtomina y Brachystomella, aunque presentes en todas las coberturas, fueron indicadores de pasturas.La cobertura de pastosfavorece la retención de humedad y controla la temperatura del suelo, variables que influyen significativamente sobre la tasa de reproducción de los colémbolos, el crecimiento de los individuos y su distribución a lo largo del perfil (^{Guzmán et al., 2010}), ya que de no existir cobertura vegetal sobre el suelo aumenta la evaporación del agua e incrementa la temperatura, con lo cual se ve afectada directamente la densidad poblacional de los colémbolos, sus ciclos de reproducción y su dinámica poblacional (^{Cole, Bradford, Shaw, & Bardgett, 2006}). Seira y Entomobrya presentan especies epigeas, adaptadas a áreas abiertas (e.g. sabanas extensivas, ^{RamosDíaz, Palacios-Vargas, & Pinzón, en prensa}), con adaptaciones morfológicas que los protegen de la radiación solar al reducir la transpiración de la superficie del cuerpo (e.g. escamas, sedas e incluso una pigmentación oscura), así como la existencia frecuente de mudas que les favorecen de condiciones extremas de humedad y temperatura (^{Guzmán et al., 2010}); o incluso a la exposición directa de la lluvia. Dichas adaptaciones resultan favorables en pasturas donde la radiación solar se da de manera directa, al no existir cobertura forestal que proteja al suelo del impacto que pueda ocasionar dicha radiación (^{Arbea & Blasco-Zumeta, 2001}). El género Brachystomella por su parte tiene hábitos fitófagos micrófagos (son capaces de alimentarse de esporas) y se les encuentra sobre suelos húmedos asociados a hongos y musgos (^{Ospina et al., 2003}). Los géneros observados en las pasturas en el área de estudio coinciden con reportes en pasturas de otras regiones de Colombia como lo registrado por ^{Ospina et al. (2003}), donde se identificaron familias como Hypogastruridae, Brachystomellidae, Neanuridae, Isotomidae, Entomobryidae, Paronellidae, Cyphoderidae, Sminthurididae, Katiannidae, Dicyrtomidae, Bourletiellidae, entre otras.

Por otro lado, los géneros Folsomides y Salina presentaron un alto valor de indicación en las plantaciones adultas.Estosgénerossuelen ser abundantes en la vegetación boscosa (^{Palacios-Vargas & Castaño-Meneses, 2003}), donde sobre las coberturas de dosel maduras se acumula mayor cantidad de materia orgánica y se desarrollan con mayor facilidad hongos y bacterias que hacen parte de la dieta de estosgéneros. El tipo de cobertura vegetal presente en un lugar define la abundancia de géneros de colémbolos afines, y específicamente para el caso de plantaciones forestales, la mayor influencia está dada por la edad de los árboles y la profundidad de sus raíces (^{Palacios-Vargas & Mejía-Recamier, 2007}), además de factores como la calidad del suelo en términos de fertilidad y las condiciones climáticas presentes (^{Rodríguez Serrano, Acevedo, Álvarez Sánchez, & Uribe Gómez, 2016}). La cantidad y calidad de hojarasca que se acumula en plantaciones adultas es un recurso favorable en la cadena alimenticia para los colémbolos, ya que, durante el proceso de descomposición de la materia orgánica, se favorece la diversidad de hongos y bacterias los cuales son su principal fuente de alimento para los éstos (^{Guillén, Soto-Adames, & Springer, 2006}). La oferta de materia orgánica disponible produce altos contenidos de ácidos húmicos, disminuyendo la concentración del pH, y generando un ambiente ideal para la reproducción de hongos. Cuando existe una mayor diversidad de hongos se genera un aumento en las asociaciones entre las plantas y las micorrizas y también se incrementa la generación de exudados radicales, los cuales también son una importante fuente de alimento para los colémbolos, los hongos y las bacterias (^{Guillén et al., 2006}). Por esta razón, al encontrar una mayor cantidad de materia orgánica en plantaciones adultas se traduce a una mayor oferta de alimento para los colémbolos y su relación con laabundancia presente en dichas coberturas.

Un mayor número de géneros fueron indicadores de bosques naturales, siendo Ptenotrix, Sphyrotheca e Isotoma exclusivos de bosques de galería.Según la literatura (^{Palacios-Vargas & Castaño-Meneses, 2003}), estos génerosse registran sobre coberturas de dosel maduros y presentan adaptaciones morfológicas para facilitar su movilidad (como la forma del ungues y el apéndice empodial de sus patas), ya que se alimentan de hifas de hongos que crecen sobre la corteza de árboles e incluso se alimentan de musgo y líquenes. Adicionalmente, el desarrollo de dosel presente en el bosque de galería disminuye la radiación solar directa sobre el suelo, manteniendo las condiciones de temperatura bajas y a su vez beneficiando la composición de colémbolos presentes en el suelo (^{Camara et al., 2012}).

Las comunidades de colémbolos en las plantaciones estuvieron compuestas tanto de géneros de áreas abiertas (pasturas), como de bosques naturales.De esta forma, géneros como Seira, Brachystomella y Arlesia pueden encontrarse tanto en pasturas como en plantaciones de E. pellita. Asimismo,se observó similitud entre las coberturas de E. pellita de 14 años y los pastos (60 %) mayor a la similitud entre el bosque de galería y la plantación (40 %), posiblemente debido a que a los 14 años de edad de la plantación la densidad de árboles es baja y permite más especies de áreas abiertas, mientras que algunas especies son especialistas del bosque natural y no se adaptan a las plantaciones. Por otra parte, la similitud de géneros entre diferentes edades en las plantaciones es significativa (70-80 %); es decir se encuentran aquellas especies generalistas de áreas abiertas que toleran áreas con dosel arbóreo, y otras generalistas de áreas de bosque que se adaptan a las plantaciones donde se mantienen a lo largo del turno de la plantación.

El modelo RDA, sugiere que la abundancia de Entomobrya responde positivamente al desarrollo de cubierta forestal (altura de los árboles), y se identificó al género Ptenotrix como el de mayor relación con el contenido de Nitrógeno que se encuentra disponible en el suelo en los bosques de galería (es decir la abundancia del género se ve afectada directamente por la disponibilidad de nitrógeno) (^{Deharveng, 1996}), mientras que Lepidocyrtus se vio favorecido con la menor cobertura del suelo (ESAT). En plantaciones jóvenes como es el caso de E. pellita de 2 y 6 años, existe una mayor exposición de suelo descubierto la cual se relaciona fuertemente con los géneros más abundantes, a medida que disminuye la exposición del suelo la abundancia aumenta.

Contrariamente a lo observado en este estudio, en otra región de la altillanura colombiana RamosDíaz et al. (en prensa), la diversidad de colémbolos aumentó en sabanas forestadas con Acacia mangium, y la abundancia disminuyó con la edad de plantación. Los resultados de ambos estudios coinciden en que la época seca es un factor limitante en la diversidad de colémbolos independientemente de la cobertura vegetal. Lepidocyrtus coincidió como uno de los géneros más abundantes en la época seca en ambos estudios. De la misma forma, según lo encontrado por Ramos Díaz et al. (en prensa), los géneros Salina y Trogolaphysa se afectan por la remoción de la materia orgánica, cambios en el microclima o la reducción de la cantidad y calidad de la hojarasca, lo cual coincide con lo observado para Lepidocyrtus cuando disminuye la disponibilidad de materia orgánica presente en el suelo en este estudio.

La forestación con E. pellita genera un hábitat favorable para algunas especies de colémbolos de áreas abiertas y otras de áreas boscosas. Por lo anterior se pueden llegar a considerar como una opción productiva intermedia entre pasturas y bosques naturales y a su vez para mantener, recuperar o incrementar la biodiversidad del lugar (^{Real de Abreu, Bessão de Assis, Frison, & Durigan, 2011}). Por ser cultivos a largo plazo y por el desarrollo de cobertura de dosel, se considera que las plantaciones forestales pueden ofrecer mayores beneficios que otros cultivos en la sostenibilidad del suelo, como ha sido observado en estudios de macrofauna en esta región (^{Decaëns et al., 2006}; ^{Lavèlle et al., 2014}; ^{Sanabria et al., 2016}), pero para el caso de la microfauna, se requiere optimizar aspectos del manejo forestal, involucrar diferencias entre especies forestales, técnicas de manejo, entre otras (e. g. densidad optima, podas). Finalmente, desde el punto de vista de la diversidad encontrada es importante mencionar que este estudio reporta 27 géneros pertenecientes a 12 familias de la clase Collembola, dentro de los cuales Cyphoderus, Arlesia, Sphaeridia, Collophora, Ptenothrix y Dicyrtomina son considerados como nuevos registros genéricos en la zona de estudio evaluada.

Declaración de ética: los autores declaran que todos están de acuerdo con esta publicación y que han hecho aportes que justifican su autoría; que no hay conflicto de interés de ningún tipo; y que han cumplido con todos los requisitos y procedimientos éticos y legales pertinentes. Todas las fuentes de financiamiento se detallan plena y claramente en la sección de agradecimientos. El respectivo documento legal firmado se encuentra en los archivos de la revista.