Estimación de volúmenes maderables de Eucalyptus globulus en el distrito de Lajas - provincia de Chota. Estimation of timber volumes of Eucalyptus globulus in the Lajas district - Chota province. Dialeny Meliza Delgado Delgado 1; Yon Gilder Herrera Cieza 1; Roxana Jackeline Rafael Sánchez 1; Miguel Ángel Villar Cabeza 2; Fátima Elizabeth Marcelo Bazán 3; Juan Rodrigo Baselly Villanueva 3. RESUMEN El Eucalyptus globulus es la especie exótica más plantada en el país, pero se desconoce la situación actual de la producción maderable de todas estas plantaciones. En este contexto se realizó un inventario forestal para cuantificar el volumen de madera de un macizo forestal de esta especie, ubicada en el distrito de Lajas provincia de Chota departamento Cajamarca. Determinándose un Volumen total (Vt) de 36.8377 m3 ha-1 y comercial (Vc) de 18.8876 m3 ha-1. Además, se realizó los cálculos del Incremento Medio Anual total (IMAt) siendo de 3.0699 m3ha-1año-1 y del Incremento Medio Anual comercial (IMAc) que es de 1.5740 m3ha-1año-1. Finalmente se utilizó el modelo de regresión Log𝑉𝑇/𝐶=𝛽0+ 𝛽1 Log (DAP x H𝑇/𝐶) + ε para estimar VT y VC, obteniéndose buenos ajustes con Coeficientes de correlación (Ryy) de 0.989 y 0.973. Palabras claves: Eucalyptus globulus, inventario forestal, Volumen total (VT), Volumen comercial (VC), Incremento Medio Anual (IMA). ABSTRACT Eucalyptus globulus is the most planted exotic species in the country, but the current situation of the timber production of all these plantations is unknown. In this context, a forest inventory was carried out to quantify the volume of wood of a forest massif of this species, located in the district of Lajas province of Chota department of Cajamarca. Determining a Total Volume (Vt) of 36.8377 m3 ha-1 and commercial (Vc) of 18.8876 m3 ha-1. In addition, calculations were made of the Total Annual Average Increase (IMAt) being of 3.0699 m3ha- 1año-1 and the Average Annual Commercial Increase (IMAc) which is 1.5740 m3ha-1año-1. Finally, we used the regression model Log𝑉𝑇/𝐶 = β0 + β1 Log (DAP x H𝑇/𝐶) + ε to estimate VT and VC, obtaining good adjustments with correlation coefficients (Ryy) of 0.989 and 0.973. Keywords: Eucalyptus globulus, forest inventory, Total volume (Vt), Commercial volume (Vc), Annual average increase (IMA). 1 Alumnos de la carrera profesional de Ingeniería Forestal y Ambiental de la Facultad de Ciencias Agrarias y Forestales de la Universidad Nacional Autónoma de Chota. 2 Responsable del PIA (Programa de Investigación Agraria Forestal) de la Estación Experimental Agraria Baños del Inca – INIA (Instituto Nacional de Innovación Agraria). 3 Equipo técnico del PIA (Programa de Investigación Agraria Forestal) de la Estación Experimental Agraria Baños del Inca – INIA (Instituto Nacional de Innovación Agraria). INTRODUCCIÓN obteniéndose esta información mediante un inventario forestal que Cuñachi (s. f), lo define El Eucalyptus globulus es originario de Tasmania, como un método de recolección y registro de datos siendo introducido al Perú en el año 1860 con la de los diferentes árboles forestales que conforman finalidad de establecer plantaciones, encontrándose el bosque o plantación, por medio de pequeñas las primeras en los departamentos de Junín, Cusco parcelas de muestreo, SERFOR (2016) señala que y Ancash (Moreno 2015, Fernández 2008); en esta actividad constituye uno de los instrumentos Cajamarca se llevó a cabo entre 1880 y 1882, más importantes para el registro de información en siendo probablemente los mojes de la congregación la evaluación de los recursos forestales, mostrando Franciscana quienes repartieron las primeras un estado situacional cuantitativo y cualitativo. plantas (Pajares s.f.). Esta especie forestal es de mayor utilización en la reforestación del país y su El análisis de regresión es un método estadístico aceptación como madera ha ido creciendo en el para el estudio de la interdependencia entre mercado nacional, por lo que constituye una opción variables, el análisis se centra en el estudio de la muy promisoria. El IGP (2012), indica que el relación matemática entre variables, con el objetivo Eucalyptus globulus genera ingresos económicos, de predecir valores de una o más variables porque sirve para leña, construcción, delimitación dependientes a partir de un conjunto de variables de linderos, muebles, para usos medicinales, entre predictoras (Cancino s.f.). Una forma de acceder al otros. Hasta el año 2013 se contó con conocimiento del volumen de madera de una aproximadamente 1 008 000 ha de superficie plantación de forma práctica y en menor tiempo es reforestada (Samaniego 2013); por la cantidad de a través de la regresión (Da Cunha et al. 2009). En plantaciones existentes David y Cárdenas (1979), los modelos volumétricos la variable predictoria mencionan que el eucalipto ocupa el 90 % de la sólo puede ser el DAP o el DAP y alturas, el uso de extensión de plantaciones en el Perú, adaptándose estas depende de la variación en la masa forestal, y a condiciones de la sierra entre los 2000 y 3500 la variable a predecir es el volumen. metros sobre el nivel del mar. Gillis y Flores (1983) En vista de la gran importancia de esta especie realizaron un estudio situacional de plantaciones de por sus beneficios, su gran extensión en el país, y eucalipto en el departamento de Cajamarca, que en la actualidad es escasa la información plantadas entre los años 1974 y 1980, determinado sobre la producción maderable. Se realizó un que no son rentables para una industria ya que inventario de una plantación ubicada en el distrito tienen poca supervivencia, mucha heterogeneidad de Lajas provincia de Chota departamento de y poco crecimiento. Cajamarca con el fin de generar información, Para saber la cantidad de madera existente en un sobre la cantidad de madera existente y la predio se realiza una evaluación de productividad situación actual en la que se encuentra. de un sitio, una de las formas más utilizadas es la determinación del volumen ya que este valor permite apreciar la ganancia producida al combinar el diámetro y la altura (Pinilla y Gonzales 1996), MATERIAL Y MÉTODOS (IMA). Se instalaron 20 parcelas rectangulares con un tamaño de 25m x15m (375 m²), dispuesta al azar Área de la investigación en función de la fisiografía del terreno. En cada una de estas se midió el DAP, Altura total (Ht) y La investigación se desarrolló en una plantación comercial (Hc), mensurándose un total de 579 forestal de Eucalyptus globulus, ubicada en el árboles. Se utilizó forcípula para la medición del distrito de Lajas provincia de Chota departamento DAP y el hipsómetro Suunto para medir alturas. de Cajamarca. Entre las coordenadas 750918 – 751201 Este y 9275475 – 9275734 Norte UTM de Los datos obtenidos en el inventario fueron la proyección WGS 84-Zona 17S, y con un rango sistematizados y procesados en el software Excel. altitudinal entre 2186 y 2268 msnm. El macizo está Para determinar las frecuencias diamétricas se en un sistema de plantación de tres bolillos, con un usaron todos los datos. En el cálculo de los distanciamiento de 3 x 3 m. Los árboles fueron volúmenes se usaron DAPs mayores de 5 cm, aprovechados a los 21 años, y ahora sus rebrotes obteniéndose los volúmenes existentes en cada tienen 12 años, los cuales no tienen manejo y en su parcela, los cuales fueron promediados y mayoría presentan gomosis. La extensión arbórea extrapolados para estimar los volúmenes para todo actual es de 3. 36 ha, con pendientes de entre 3 y el predio. 45 %, Para el ajuste se utilizó modelo de Spurr Según, Valdivia et al. (s. f), indica que el distrito de linearizado con una modificación, donde el DAP Lajas se ubica a una latitud de 6°33´ Sur, longitud no fue elevado al cuadrado, quedando de la 78°44´Oeste y su temperatura media anual se siguiente manera: encuentra entre 15. 9 ºC a 21. 0 º C. Según Sánchez Log𝑉𝑇/𝐶=𝛽0+ 𝛽1 Log (DAP x H𝑇/𝐶) + ε y Sánchez (2010), de acuerdo a la clasificación de Javier Pulgar Vidal, Chota – Lajas se encuentra Siendo: 𝑉𝑇/𝐶=𝑉𝑜𝑙𝑢𝑚𝑒𝑛 𝑡𝑜𝑡𝑎𝑙 𝑜 𝑐𝑜𝑚𝑒𝑟𝑐𝑖𝑎𝑙, en m 3. dentro de la región Quechua baja, con un clima 𝐷𝑎𝑝=𝐷𝑖á𝑚𝑒𝑡𝑟𝑜 𝑎 𝑙𝑎 𝑎𝑙𝑡𝑢𝑟𝑎 𝑑𝑒𝑙 𝑝𝑒𝑐h𝑜, en cm. fresco. H𝑇/𝐶=𝐴𝑙𝑡𝑢𝑟𝑎 𝑡𝑜𝑡𝑎𝑙 𝑜 𝑐𝑜𝑚𝑒𝑟𝑐𝑖𝑎𝑙, en m. Metodología 𝛽0, 𝛽1, 𝛽2= Constantes propias del modelo de Usando un GPS se georreferenció el perímetro de regresión. la masa boscosa y el centroide de cada parcela. Los datos obtenidos, coordenadas, fueron procesados Usándose solamente las mediciones de los en el programa ArcGIS, calculando el área total de individuos con DAPs mayores a 10 cm para el la plantación y se verificó la adecuada distribución ajuste, los cuales fueron procesados el Software de las parcelas. Excel, mediante el método de regresión lineal. Para análisis de varianza de la regresión se usó la prueba El inventario de la plantación se realizó con la de “F” y “t” a una significancia de ∝ ≤ 0.05 y un finalidad de determinar Volumen total (Vt), 95 % de confianza. El grado de los ajustes fueron Volumen (Vc) y el Incremento Medio Anual determinados usando el Coeficiente de correlación (𝑅?̂?𝑌), Coeficiente de determinación (𝑅 2), Error m3ha-1, un IMAt de 3.0699 m3ha-1año-1 y IMAc de padrón residual 𝑆𝑌𝑋 , gráfica de distribución de 1.5740 m 3ha-1año-1. Además, todo el predio tiene residuos y su histograma. 123.77 m3 de VT y 63.46 m3 de VC. Tabla 2. Variables dasométricas por parcela. N° RESULTADOS DAP AB VT VC IMAt IMAc Parcela árbo (cm) (m2) (m3) (m3) (m3 año-1) (m3 año-1) les Frecuencias diamétrica obtenidas del inventario 1 9 8.3 0.0551 0.3596 0.2112 0.0300 0.0176 forestal 2 17 8.7 0.1155 0.7957 0.4425 0.0663 0.0369 3 13 9.6 0.1099 0.7009 0.3744 0.0584 0.0312 En la Tabla 1 se muestran 9 clases diamétricas, en 4 21 8.7 0.1509 1.2615 0.5548 0.1051 0.0462 5 40 9.4 0.3149 2.3274 1.2630 0.1939 0.1053 donde se observa que el mayor número de árboles 6 20 9.3 0.1551 1.1787 0.6779 0.0982 0.0565 se encuentra en los intervalos 4.5 – 7.5 cm con 262 7 9 7.2 0.0374 0.2050 0.0909 0.0171 0.0076 8 20 8.3 0.1254 0.7623 0.2324 0.0635 individuos, representando un 45.25 %. La clase 0.0194 9 27 7.7 0.1378 0.4977 0.1447 0.0415 0.0121 diamétrica donde se observa la menor cantidad de 10 11 8.6 0.0631 0.4273 0.2440 0.0356 0.0203 árboles se encuentra en el intervalo 32.5 – 35.5 cm 11 24 7.9 0.1288 0.8116 0.4452 0.0676 0.0371 12 16 8.1 0.0908 0.4847 0.2413 0.0404 0.0201 con un solo árbol, representado un 0.17 %. 13 20 9.3 0.1564 1.1804 0.6522 0.0984 0.0544 14 10 9.3 0.0774 0.5482 0.2900 0.0457 0.0242 15 33 10.7 0.3390 2.8634 1.6095 0.2386 Tabla 1. Frecuencias diamétricas de 0.1341 16 25 9.0 0.1752 1.2537 0.7262 0.1045 0.0605 Eucalyptus globulus. 17 24 9.4 0.1161 0.7511 0.4255 0.0626 0.0355 Intervalos de clase (cm) Marca de clase N° de Frecuencia 18 25 12.1 0.3346 3.0308 1.5105 0.2526 0.1259 (cm) árboles relativa (%) Li Ls 19 27 11.5 0.3727 3.2676 1.5202 0.2723 0.1267 20 24 13.2 0.4590 4.9207 2.5093 0.4101 0.5 3.50 2.00 77 13.30 0.2091 Prom. 21 9.32 0.1758 1.3814 0.7083 0.1151 0.0590 4.50 7.50 6.00 262 45.25 Des. 8.12 1.52 0.1200 1.247 0.6410 0.104 0.0534 8.50 11.50 10.00 141 24.35 Estand. Max. 40 13.2 0.4590 4.9207 2.5093 0.4101 0.2091 12.50 15.50 14.00 65 11.23 Min. 9 7.2 0.0374 0.205 0.0909 0.0171 0.0076 16.50 19.50 18.00 21 3.63 C.V. 39.1 16.26 68.16 90.26 90.50 90.26 90.50 20.50 23.50 22.00 6 1.04 24.50 27.50 26.00 3 0.52 Ajuste del modelo matemático para estimar 28.50 31.50 30.00 3 0.52 32.50 35.50 34.00 1 0.17 volúmenes SUMA 579 100 Al realizar las regresiones con todos los datos del inventario, estas no fueron significativas; y los Determinación del volumen maderable en las individuos con DAPs menores de 10 cm parcelas inventariadas presentaron una gran subestimación de sus En la Tabla 2 se observan los resultados obtenidos volúmenes. Por lo cual se trabajó con los rebrotes en cada una de las parcelas evaluadas, a partir de que tenían DAPs mayores de 10 cm obteniendo los cuales se obtuvieron valores promedio para regresiones estadísticamente significativas, prueba cada variable dasométrica evaluada. Extrapolando de “F” y “t”, para los dos volúmenes como se tales valores, es posible estimar que el macizo observa en la Tabla 3. Su grado de ajuste fue presenta un VT de 36.8377 m3ha-1 y VC de 18.8876 bueno, Tabla 4, pues los menores valores de 𝑅?̂?𝑌 y 𝑅2 son 0.973 y 0.946 respectivamente; además el Tabla 4. Indicador estadístico del grado del mayor 𝑆𝑌𝑋 fue de 0.024. En las gráficas de ajuste. Volumen Volumen distribución de residuos se observa que existe sobre Variable dasométricas total (Vc) comercial (Vt) y sub estimación en diámetros menores, pero con 𝑅?̂?𝑌 0.989 0.973 Indicadores 𝑅2 0.978 0.947 forme el DAP aumenta sólo hay subestimación de 𝑆 0.024 0.022 𝑌𝑋 los volúmenes (Gráfico 1-a); por otro lado, en el histograma de residuos se determina que los ajustes Tabla 5. Parámetros del modelo para los en su mayoría de estimaciones se aproximan al 0 % volúmenes. y tiene casi igual número de sobre y subestimación Volumen Volumen Variable dasométricas total comercial (Gráfico 1-b). Producto de las regresiones se (Vc) (Vt) determinaron los valores de los parámetros 𝛽 y 𝛽 β0 -4.4208 -4.0846 0 1 Parámetros β1 1.5663 1.4724 que se muestran en la Tabla 5. Tabla 3. Análisis de varianza de la regresión. DISCUSIÓN Volumen Volumen Variable dasométricas total (Vc) comercial (Vt) Campos y Leite (2013) mencionan que la Regresión 1 1 distribución diamétrica de una plantación presenta Grados de Residuos 137 137 libertad una distribución unimodal o normal, con ligera Total 138 138 frecuencia de diámetros a la izquierda y derecha en "F" 6085.980 2436.13994 Calculado β0 -98.418 -69.4762331 el eje X, y en bosque naturales la distribución de "t" β1 78.012 49.3572684 frecuencias sigue una tendencia decreciente con el "F" 3.984 Tabulado aumento del DAP. "t" 1.960 Gráfico 1. Análisis de los residuos de las estimaciones. Volumen total (Vc) Volumen comercial (Vt) 100 100 80 80 60 60 40 40 20 20 Distribución de 0 0 residuos -20 -20 (a) -40 -40 -60 -60 -80 -80 -100 -100 10 20 30 40 10 15 20 25 30 35 DAP (cm) DAP (cm) 60 40 50 30 40 30 Histograma de 20 residuos 20 (b) 10 10 0 0 Clases de error (%) Clases de error (%) N° de individuos % de error N° de individuos % de error Menez et, al. (2015) realizó un estudio en algunos o todos los factores mencionados por lo Eucalyptus urophylla, de 8 años de edad que esta plantación se encuentra en esta situación; obteniendo los mayores individuos en las clases como la falta de manejo de los rebrotes que ha diamétricas centrales; en este estudio las primeras generado una gran heterogeneidad en el desarrollo clases diamétricas tuvieron el mayor número de de los individuos y poca producción maderable. rebrotes y conforme el DAP aumenta el número de Los ajustes con todos los datos no tuvieron valides ellos disminuye, esto nos indica que esta masa estadística, pues los rebrotes con diámetros forestal tiene una mala estructura, pareciéndose a pequeños tienen mínimo volumen maderable y no la de un bosque natural y no de una plantación; al presentan una tendencia definida; Barrios et. al. presentar diámetros pequeños el valor comercial de (2014) realizó un estudio con 101 árboles de los individuos sería como leña y postes. Eucalyptus grandis con una significancia de ∝ < El eucalipto es una especie de rápido crecimiento, 0.01, y un 99 % de confianza; obteniendo un valor pero se reporta una gran variedad de cifras para su de 0.993 para el 𝑅2 y de 0.0248 para el 𝑆𝑌𝑋. Los IMA; en el suroeste de España fue de 5 a 6 m3ha- indicadores estadísticos obtenidos en el presente 1año-1, y en el noroeste 20 m3ha-1año-1, por otro estudios se encuentran cercanos a los anteriores lado, en Uruguay fue de 25 m3ha-1año-1 y en valores mencionados, el 𝑅 2?̂?𝑌 de 0.973, 𝑅 de 0.946 Portugal a la edad de 10 años fue de 20 m3ha-1año- y 𝑆𝑌𝑋 de 0.024. 1, siendo considerados estos tres últimos como muy Prodan et. al. (1997) denomina funciones generales buenos (Skolmen y Ledig 1990). A pesar de la gran a las que usan dos mediciones de los árboles para variación registradas los valores de crecimiento estimar el volumen y son recomendadas en grandes obtenidos para los rebrotes en estudio fueron muy extensiones. Los datos usados provienen de un área inferiores a los reportados en otros lugares, siendo pequeña, pero tienen alta variación (Tabla 6), a para el IMA de 3.0699 m3ha-1t año-1 y IMAc de pesar de eso se logró ajustar el modelo ya que usa 1.5740 m3ha-1año-1. dos variables independientes, DAP y Alturas, que Samaniego (2013), indica que una planta de generan más combinaciones para estimar la eucalipto con gomosis es índice de estrés por variable dependiente, volumen. Además, las factores abióticos como sequía, fuego, heladas, o estimaciones del Volumen total fueron mejores reacción frente a agentes bióticos tales como pues las Alturas totales tuvieron menor variación insectos u hongos; en la plantación la mayoría de que las comerciales. los brotes presentan está enfermedad indicando la Tabla 6. Análisis estadístico de los datos existencia de estrés. Según Gillis y Flores (1983) dasométricos usado en las regresiones. las causas para el fracaso de una plantación son el DAP Ht Hc usó de semillas y plantones de baja calidad, la mala Promedio 14.14 12.74 6.99 selección de la especie a plantar en función a la Mínimo 10.00 7.00 2.75 Máximo 32.30 21.50 15.00 calidad de sitio, pocas a nulas o mal ejecutadas Varianza 20.17 7.11 4.10 labores silviculturales en el desarrollo de los Des. estándar 4.49 2.67 2.02 individuos hasta su aprovechamiento; siendo C.V. 31.75 20.93 28.94 CONCLUSIÓN Cuñachi, EG. s. f. Manual práctico de inventarios forestales. CORPIAA Atalaya Ucayali, Perú, s.e. 3 Para obtener una buena rentabilidad maderable se p. deben usar semillas y plantones de calidad, hacer Da Cunha, TA; Guimarães Finger, CA. 2009. Modelo de una buena selección de la especie a plantar en regresión para estimar el volumen total con corteza función a la calidad de sitio y realizar labores de árboles de Pinus taeda L. en el sur de Brasil. Kurú: Revista Forestal. 6(16) 2-15. silviculturales que ayuden el desarrollo de los individuos. La plantación estudiada se encuentra en David, E; Cardenas, H. 1979. Turno Financiero en Plantaciones de Eucalyptus globulus Labill del Valle muy mal estado, teniendo la madera sólo un valor del Mantaro. Revista forestal del Perú. 9(2): 40-52. comercial como leña o postes. Por la edad de los individuos ya no se podría hacer manejo, por lo Fernández A, AF. 2008. Estudio de Prefactibilidad para la Producción y Comercialización de Papel a partir cual la solución es la tala total de la plantación y de Eucalipto. Tesis Ing. Lima- Perú. 105 p. manejar los nuevos rebrotes o la instalación total con nuevos individuos teniendo las Gillis M; Flores T, F. 1985. Las plantaciones forestales en Cajamarca. Baños del Inca, Cajamarca, Perú, consideraciones anteriormente mencionadas. CICAFOR (Centro de Investigación y Capacitación El usó de individuos con diámetros menores de 10 Forestal). 94 p. cm dificulta la regresión, el modelo Log𝑉𝑇/𝐶=𝛽0+ IGP (Instituto Geofísico del Perú). 2012. Manejo de 𝛽1 Log (DAP x H𝑇/𝐶) + ε puede ser usado en áreas Riesgos de Desastres ante Eventos Meteorológicos pequeñas con alta variación pues presenta más Extremos en el Valle del Mantaro. Perú. 197 p. combinaciones para la variable respuesta. Cuando Menez M, FC; Imaña E, J; Pereira, RS; Rezende, AV; se desea determinar el mejor modelo a usar de Eder Pereira, Miguel. 2015. Modelaje de la varios, se recomienda primero realizar la regresión Distribución Diamétrica de Eucalyptus urophylla a con el Volumen comercial pues tendrá mayor error través de las Funciones Normal, Weibull y sb de estimación a comparación del Volumen total. Johnson en la Región de Brasilia. Brasil. 15 p. Moreno R, L. 2015. Manejo e Industrialización del Eucalipto (Eucalyptus globulus). Tesis Ing. Trujillo- REFERENCIAS BIBLIOGRÁFICAS Perú. 36 p. Pajares G, U. s.f. El Eucalyptus globulus Labill sub- Barrios, A; López, A; Nieto, N. (2014). Predicción de globulus en la sierra de Cajamarca. Baños del Inca, Volúmenes Comerciales de Eucalyptus Grandis a Cajamarca, Perú, CICAFOR (Centro de Través de Modelos de Volumen Total y de Razón. Investigación y Capacitación Forestal). 13 p. Colombia. 137 p. Pinilla S, JC; Kahler G, C. 1996. Determinación de Campos, JCC; Leite, HG. 2013. Mensuração Florestal: Volúmenes para Eucalipto. Chile. 21p. Perguntas e Respostas. 4 ed. Viçosa, MG, BR. UFV. 605 p. Prodan, M; Peters, R; Cox, F; Real, P. 1997. Mensura Forestal. San José, CR. IICA. 586p. Cancino, J. s. f. Dendrometría Básica. s. l, s.e. 65. p. Samaniego M, CA. 2013. Efecto de un incendio forestal Skolmen, R.G.; LEDIG, F. T. 1990. Eucalyptus globulus en una plantación de Eucalyptus globulus Labill. Labill. Bluegum eucalyptus. En: Burns, Russell M.; subsp. Globulus en Huaraz. Tesis Ing. Lima-Perú, Honkala, Barbara H., eds. Silvics of North America: UNALM. 113 p. 2. Hardwoods. Agric. Handb. 654. Washington, DC: U.S. Department of Agriculture, Forest Service: 299- Sánchez R, A; Sánchez T, S. 2010. Mapa de Pisos 304p. Altitudinales. Cajamarca, Perú. 13 p. Valdivia V, FM; Rojas Q, AF; Comaila R, E. s. f. Estudio SERFOR (Servicio Nacional Forestal y de Fauna Definitivo para la Rehabilitacion y Mejoramiento de Silvestre, Perú). 2016. Primer Informe Parcial del la Carretera Chongoyape- Cochabamba- Inventario Nacional de Flora y Fauna Silvestre. Cajamarca Tramo Cochabamba-Chota. Cajamarca- Informe de 2016. La Molina, Lima, Perú. 268p. Perú. 143 p.