Abstract
The inventory of dung beetles in Colombia has a relatively long sampling history, represented in local checklist and occasional survey records. However, there is no diagnosis on the completeness and environmental representativeness of Scarabaeinae samplings both at regional and national levels. A study like this is essential to address questions related with the distribution of species as well as their responses to different disturbance scenarios. In this regards and based on the available taxonomic records of Scarabaeinae for the Ecoregion of the Eje Cafetero (Central Western Andes) the present study aims to estimate the reliability and thoroughness of the taxonomic inventories for this group, suggesting additionally potential areas where future samplings should be focused. Scarabaeinae surveys in the Ecoregion comprise 36 years of discontinuous samplings and the number of species increases linearly since 1990s. Nevertheless, the historical sampling covers only 4% of the ecoregion due to the marked spatial and altitudinal aggregation of the records. In addition, it has been detected a high level of uncertainty and taxonomic inconsistency between the local checklists, which reduces the reliability of the regional inventories. It is recommended that future surveys should be focused on inter-Andean valleys and on the eastern slope of Western and Central Cordillera. On the other hand, it is necessary to strengthen the protocols for the proper storage of data, also focusing efforts towards the taxonomic unification of local species checklists.
Introducción
Los escarabajos coprófagos son comúnmente usados para evaluar el impacto de la actividad humana sobre la biodiversidad (Nichols
Gran parte de la información disponible sobre diversidad a grandes escalas, deriva de la acumulación histórica de inventarios locales y regionales, así como de colectas ocasionales. Actualmente, una parte de dicha información se encuentra condensada en Bases de Datos de Biodiversidad las cuales a pesar de sus limitaciones, permiten evaluar la completitud y representatividad de los datos a lo largo de diferentes escalas espaciales (Lobo 2008). En el caso de los escarabajos coprófagos de Colombia, la construcción de la base de datos es un proceso reciente y mucha de la información secundaria (e.g. coordenadas geográficas) aún no se encuentra disponible (SIB Colombia 2012). Consecuentemente a pesar de la tradición de muestreo de Scarabaeinae en el país, no hay fiabilidad taxonómica y consenso sobre la exhaustividad del inventario nacional; falta asimismo una estimación sobre la magnitud y localización de los principales vacíos de información.
Este trabajo tiene como objetivo evaluar la completitud y determinar los sesgos espaciales de los datos disponibles sobre los escarabajos coprófagos en el ámbito geográfico de la ecorregión del Eje Cafetero (Andes Centro-Occidentales de Colombia). Esta ecorregión concentra gran cantidad de registros taxonómicos de Scarabaeinae, los cuales han sido la base para la unificación taxonómica de los distintos inventarios locales (Cultid
Materiales y métodos
Región de estudio.
Este estudio comprende la denominada Ecorregión del Eje Cafetero, la cual corresponde a gran parte de los Andes Centro Occidentales de Colombia; es una delimitación socio-económica que abarca cerca de 3,2 millones de ha e incluye los departamentos con mayor producción de café en Colombia: Risaralda, Caldas, Quindío, Valle del Cauca (norte) y Tolima (noroccidente) (Fig. 1). La Ecorregión del Eje Cafetero cubre 11% del área de la región andina colombiana, con un gradiente altitudinal que oscila entre los 250 m (valle interandino del río Magdalena) y los 5.300 m (Parque Nacional Natural Los Nevados). En esta Ecorregión se han identificado 18 zonas de vida (Espinal 1967) distribuidas en seis subregiones andinas definidas por las vertientes de cada cordillera y los valles interandinos de los ríos Cauca y Magdalena (Fig. 1) (modificado de Kattan y Franco 2004).
Además de la importancia económica de la ecorregión, por concentrar gran parte de la producción de café nacional, se destaca su articulación geo-política bajo el Sistema Regional de Área Protegidas del Eje Cafetero (SIRAP-EC 2013). En este sentido, la delimitación y consolidación del SIRAPEC ha impulsado de forma directa e indirecta el desarrollo de inventarios de diversidad biológica y de estudios de impacto ambiental. En ellos se incluyen diferentes grupos taxonómicos, algunos de ellos considerados indicadores ecológicos de perturbación antrópica (Rojas-Diaz

Ecorregión del Eje Cafetero - Sistema Regional de Áreas Protegidas del Eje Cafetero (SIRAP - EC). Se indican las subregiones andinas, definidas por la combinación de las siguientes convenciones: Vor = vertiente oriental, Voc = vertiente occidental, CO = cordillera Occidental, CC = cordillera Central VI = valles interandinos de los ríos Cauca y Magdalena (Modificado de Cultid
Registros taxonómicos.
La información usada en este estudio proviene de la base de datos actualizada, con la cual se elaboró la lista regional de escarabajos coprófagos (Cultid
Los registros en esta base de datos proceden de dos fuentes: una parte derivan de muestreos sistemáticos en los cuales se empleó un método estandarizado de colecta; otros corresponden a registros de muestreos ocasionales con diferentes niveles de información asociada (e.g. tipo de cebo, hábitat de muestreo, ubicación geográfica). Aquellos registros que no especificaban coordenadas geográficas, se georeferenciaron con base en información secundaria (e.g. informes técnicos), el gacetero GeoNames (ver http://www.geonames.org) y Google Earth. Los registros sin datos geográficos mínimos (departamento, municipio, vereda y localidad específica) no se incluyeron en el análisis.
Datos cartográficos y ambientales.
El análisis espacial de la riqueza de Scarabaeinae incluyó tres fuentes de información: i) base cartográfica de la Ecorregión (límite geográfico y límites políticos; fuente WCS-Programa Colombia), ii) zonas de vida siguiendo el sistema de clasificación de Holdridge para la Ecorregión del Eje Cafetero (fuente WCS-Programa Colombia) y iii) polígonos de las áreas protegidas que componen el SIRAP-EC. Toda la información se procesó con ArcGIS 9,3 (ESRI Inc. 1999-2008).
Análisis de datos.
Para evaluar la distribución espacial de los registros de escarabajos coprófagos, el área de la Ecorregión del Eje Cafetero se dividió en celdas de 5 x 5 km. Este tamaño de celda, aunque burdo para los requerimientos microambientales de los insectos, constituye a nuestro juicio, la resolución más detallada que es posible utilizar teniendo en cuenta la precisión geográfica de las localidades de colecta. Dicha malla se construyó teniendo en cuenta el modelo digital del terreno disponible en WorldClim (WorldClim 2013) y el programa ArcGIS (ESRI Inc. 1999-2008). Posteriormente con base en los registros taxonómicos de la base de datos, se procedió a determinar la riqueza observada en cada celda.
La completitud de muestreo en cada celda de 5 x 5 km se evaluó mediante tres estimadores no paramétricos de riqueza (ICE, Chao 2, Jack 1) recomendados para casos en donde las unidades de muestreo son pequeñas (Hortal
Los estimadores de riqueza se obtuvieron con el programa EstimateS 8,2 (Colwell 2006) y las celdas con un porcentaje de completitud ≥ 80% por cualquiera de los tres estimadores utilizados fueron consideradas "bien muestreadas" (Lobo 2008). La completitud se estimó también a escala de ecorregión y para cada una de las sub-regiones andinas. Para establecer la representatividad ambiental de la información existente, se estudió la distribución de los registros entre bandas altitudinales por subregión andina y zonas de vida. En adición, se estableció la riqueza observada y el número de visitas asociadas a las áreas protegidas del SIRAP-EC.
Debido a que los datos provenientes de colecciones entomológicas generalmente contienen un subconjunto de las listas de especies publicadas (Lobo 2008), se examinó la capacidad de la base de datos compilada en este estudio para estimar la riqueza de los escarabajos coprófagos a escala de la ecorregión. Para ello se revisaron las listas locales disponibles en artículos, informes técnicos, trabajos de grado y tesis. De estas listas se extrajeron los siguientes datos: número de localidades, riqueza observada, coberturas vegetales de muestreo, porcentaje de incertidumbre taxonómica y rango altitudinal. El porcentaje de incertidumbre taxonómica corresponde a la proporción de entidades determinadas solo hasta género o que están bajo alguna categoría de confirmación (cf., ca., aff.). El porcentaje de incertidumbre también se estimó para la lista de especies obtenida a partir de los datos de museo.
Resultados
Los datos compilados abarcan un periodo discontinuo de 36 años de muestreo (1974-2010) e información sobre 83 entidades taxonómicas (55 especies y 28 morfo-especies) pertenecientes a 19 géneros de Scarabaeinae (Tabla 1). Esto representa el 19% de las especies registradas para Colombia (Medina
El análisis subregional muestra una alta agregación espacial de los registros (Tabla 2). La vertiente occidental de la cordillera Central y la vertiente oriental de la Occidental concentraron el 54% y el 34% del total de los registros, respectivamente. Estos registros se distribuyeron en un rango altitudinal estrecho, entre los 1.500 y 2.250 m de elevación (Tabla 2). Respecto a la completitud de muestreo, las subregiones se ordenaron de las siguiente manera: Vor- Co > Voc - CO > Voc-CC>VI-Mag. > VI-Cau. > Vor-CC (Tabla 2). Sin embargo, el valle interandino del río Magdalena (VI – Mag.) y la vertiente occidental de la cordillera Central (Voc - CC) concentran, cada una, más del 50% de la riqueza de especies de la ecorregión (Tabla 2).
El muestreo de Scarabaeinae ha incluido siete de las 18 zonas de vida de la ecorregión, entre las cuales, el bosque húmedo tropical premontano (bht PM) y el bosque muy húmedo montano bajo (bmh-MB) poseen la mayor riqueza observada. Esta última zona de vida concentra la mayor cantidad de visitas (Fig. 4). Por otro lado, 11 de las 128 áreas protegidas que componen el SIRAP-EC han sido muestreadas en alguna ocasión (Fig. 5). Estas 11 áreas protegidas están distribuidas en cinco categorías de conservación y todas se encuentran dentro de la jurisdicción del departamento de Risaralda, excepto Barbas (ACM) y Barbas-Bremen (PRN) al norte del Quindío (Fig. 5). No se observó un patrón definido entre el número de visitas y la riqueza observada a lo largo de las áreas protegidas. En detalle, el Parque Municipal Natural Verdum presentó la mayor riqueza (25 especies), seguido por el Santuario de Fauna y Flora Otún Quimbaya (SFF: 16 especies), Ucumarí (PRN: 13 especies) y Campoalegre (PMN: 13 especies); las otras áreas protegidas muestran una riqueza observada menor a cinco especies (Fig. 5).
Se compilaron ocho listas locales de especies distribuidas entre publicaciones, informes técnicos y trabajos de grado (Tabla 3). En su conjunto, las listas locales generaron un catalogo regional compuesto por 136 entidades taxonómicas (42 especies y 94 morfo-especies) con un porcentaje de incertidumbre taxonómica que osciló entre el 19% y 80% (Tabla 3). Dichas listas se distribuyeron en 15 de las 64 celdas identificadas previamente, de esta forma, las listas locales han cubierto aproximadamente 1% de la Ecorregión del Eje Cafetero.
Riqueza de especies observada (Sobs) y riqueza de especies estimada por cada uno de los tres estimadores utilizados. En paréntesis el porcentaje de completitud para las 22 celdas de 5 x 5 km con más de tres registros de Scarabaeinae. Gobs = riqueza de géneros (Gobs), Nv = número de visitas y % Srobs = porcentaje de la riqueza regional representada en cada celda. TOTAL hace referencia a la riqueza observada para toda la ecorregión incluyendo la totalidad de las celdas con datos de riqueza (n = 64). Con un asterisco se señalan aquellas celdas con una completitud de muestreo ≥80% para al menos uno de los estimadores.

Acumulación de especies de Scarabaeinae a lo largo del periodo de tiempo que abarca la base de datos en la Ecorregión del Eje Cafetero (1974 y 2010).

Distribución de la riqueza de especies observada en la Ecorregión del Eje Cafetero. Se muestran las celdas de 5 x 5 km con al menos un registro, las celdas marcadas con un circulo (negro o blanco) presentaron una completitud igual o mayor al 80% (ver Tabla 1).
Discusión
La concentración espacial, y casi idiosincrática, de los muestreos realizados y la falta de articulación de los inventarios locales a iniciativas de muestreo a mayores escalas, son características comunes que aparecen cuando se examina el muestreo histórico de una región o país (Hortal
Número de visitas (Nv) para cada una de las subregiones y rangos altitudinales. Sobs es el número de especies observadas, los datos para cada estimador (ICE, Chao2 y Jack 1) corresponden al total de especies estimada y en paréntesis se muestra el porcentaje de completitud de muestreo. Subregiones andinas: Vor = vertiente oriental, Voc = vertiente occidental, CO = cordillera Occidental, CC = cordillera Central VI = valles interandinos de los ríos Cauca y Magdalena.

Número de especies observadas (Sobs) y número de visitas (Nv) por zona de vida. bs-PM = Bosque Seco Premontano; bmh-M = Bosque Muy Húmedo Montano; bh-MB = Bosque Húmedo Montano bajo; bmh-PM- Bosque Muy Húmedo Premontano; bmh-MB = Bosque Muy Húmedo Montano bajo; bh-PM Bosque Húmedo Premontano y bh-T = Bosque Húmedo Tropical.

Riqueza de especies observada (Sobs) y número de visitas (Nv) por área protegida, PMN = Parque Municipal Natural; SFF = Santuario de Fauna y Flora; PRN = Parque Regional Natural; ACM = Área de Conservación y Manejo; AME = Área de Manejo Especial.
Listas locales de especies de escarabajos coprófagos compiladas para la Ecorregión del Eje Cafetero, se muestra la altitud, el código de la celda asociada a la lista (ver Tabla 1) y la riqueza observada en cada celda, Nl = número de localidades; Gobs = géneros observados; Sobs = riqueza observada; Cvs = coberturas vegetales, se usa la nomenclatura usada por cada trabajo (B = bosque maduro, Bcr = Bosque de cerro, Bs = Bosque secundario, Bc - Bosque de cañada, Pb = Parches de bosque, Bu = Bosque de urapán, Ro robledal, BPo borde bosque potrero, Pe = Plantación de eucalipto, Po potrero); %IT = porcentaje de incertidumbre taxonómica; Tipo (P = lista publicado, Tg = trabajo de grado y IT = informe técnico).
Durante el proceso de consolidación del Sistema Regional de Áreas Protegidas del Eje Cafetero (SIRAP-EC), gran parte de los muestreos se realizaron entre Risaralda y Quindío, departamentos que concentran gran parte de las áreas protegidas del SIRAP-EC (Fig. 1). En el caso de los escarabajos coprófagos, los muestreos se desarrollaron como inventarios rápidos de diversidad que tenían como objetivo comparar áreas de bosque nativo con coberturas vegetales de origen antrópico (e.g. Rojas-Díaz
La completitud y representatividad ambiental en la ecorregión enfrenta un problema adicional, el alto grado de incertidumbre taxonómica, la cual de hecho, es una limitación de magnitud nacional para el estudio comparativo de Scarabaeinae a través de diferentes escalas espaciales. Aunque Colombia cuenta con una larga tradición de muestreo de escarabajos coprófagos (e.g. Howden y Nealis 1975; Escobar y Chacón 2000; Neita y Escobar 2012), la investigación en taxonomía o sistemática de Scarabaeinae ha avanzado de forma desigual con respecto al trabajo ecológico. En consecuencia, aún son escasas las herramientas taxonómicas (e.g. claves a especie) que permitan reducir la incoherencia entre listas locales de especies. De hecho, apenas en la pasada década se han publicado revisiones de géneros que incluyen especies presentes en los Andes colombianos (e.g. Oxysternon: Edmonds y Zídek 2004; Genieridium: Vaz de Mello 2008; Deltochilum: González
Si estamos lejos de finalizar el inventario completo regional de Scarabaeinae para la Ecorregión del Eje Cafetero, más lejos aún estamos de obtener una cobertura espacial que permita investigar las causas de las variaciones de la riqueza o las diferencias en composición entre localidades. Sólo un 0.3% de las celdas de 25 kmβ, contendría inventarios fiables. La ubicación de estas celdas, junto a los datos sobre la representatividad ambiental de la información utilizada, permite obtener algunas pautas para identificar áreas que deben ser objetivos prioritarios para futuros muestreos. En primer lugar, al igual que sucede a escala nacional, los bosques secos tropicales, ubicados en los valles inter-andinos, son las zonas más vulnerables y menos conocidas de la ecorregión (IAvH 1998; Arcila
Conclusión
Este estudio indica que el inventario de escarabajos coprófagos en el ámbito geográfico de la Ecorregión del Eje Cafetero está lejos de ser completo y representativo. En este momento, la información disponible no es fiable para la estimación de la riqueza de Scarabaeinae a través de diferentes escalas y consecuentemente, no es posible abordar preguntas ecológicas que trasciendan la escala local. En detalle, las principales debilidades del inventario a escala de la ecorregión son: sub muestreo, alta concentración espacial de los registros y baja coherencia taxonómica entre listas locales. Para suplir el vacío de información, los futuros inventarios deben dar prioridad a los valles interandinos de los ríos Cauca y Magdalena, y a la vertiente oriental de las cordilleras Occidental y Central.
El caso de la Ecorregión del Eje Cafetero sirve como base para reconocer que a escala nacional, es necesario fortalecer los siguientes aspectos en el estudio de Scarabaeinae: (1) inventarios regionales que consideren diferentes escalas espaciales bajo diseños anidados de muestreo (e.g. Numa
Footnotes
Agradecimientos
Este trabajo se realizó con apoyo del programa de Doctorados Nacionales del Departamento de Ciencia y Tecnología COLCIENCIAS, entidad que financia el estudio doctoral de Carlos A. Cultid. Agradecemos al Banco de la República (Colombia), al Instituto Alexander von Humboldt, a la Fundación Jhon & Catherine MacArthur, a la Corporación Autónoma de Risaralda (CARDER), al Centro Nacional de Investigación de Café (CENICAFÉ) y a la Federación Nacional de Cafeteros por su apoyo económico y logístico en las diferentes etapas del estudio. Este trabajo se desarrollo en el marco del proyecto "Análisis de heterogeneidad espacial y diversidad de hormigas y escarabajos coprófagos en paisajes cafeteros" financiado por COLCIENCIAS, la Universidad del Valle y Wildlife Conservation Society (WCS) - Colombia (Proyecto no.110652128706). Gracias a los evaluadores anónimos, a Vladimir Rojas y a Larry Jiménez por sus comentarios los cuales ayudaron a mejorar el documento.
