INTRODUCCIÓN

Según el Minam (2015), los impactos del cambio climático son hechos innegables y evidentes que afectan a los ecosistemas, diversidad biológica, la base de recursos naturales, su distribución geográfica, la estructura productiva y de servicios, así como la infraestructura económica y social, la salud y bienestar de la población, familias e individuos. En zonas altoandinas los efectos a la biodiversidad son aún mayores, ya que según ^{Castaño (2002)} la fragilidad del espacio reducido y específico se vuelve un agravante para la flora y la fauna. ^{Kappelle y Horn (2005)} indican que cuando las condiciones no permiten la migración altitudinal de las especies, el ecosistema es igual de vulnerable que una isla pequeña incrementando el riesgo a la extinción.

Se considera que una especie de flora se encuentra amenazada cuando sus poblaciones silvestres se reducen extremadamente por el desarrollo de actividades humanas o por efecto del cambio climático, y son protegidas estatalmente con base en las clasificaciones en el que se incluyen criterios de la UICN (2012) como son: peligro crítico (CR), en peligro (EN), vulnerable (VU) y casi amenazado (NT) como la lista oficial de especies de flora categorizada del Perú (Decreto Supremo N.º043-2006-AG, 2006). Además, existe el Listado de especies de flora silvestre Cites que tiene la finalidad de velar por que el comercio internacional de especímenes silvestres no constituya una amenaza para su supervivencia, en el cual se ofrecen diversos grados de protección a más de 30 000 especies en el ámbito mundial (Minam, 2018).

Este grado de amenaza puede incrementarse debido al aumento de emisiones de GEI que el IPCC (2013) las define como las cuatro trayectorias de concentración representativas (RCP por sus siglas en inglés) en función de la concentración de CO₂ en la atmosfera (forzamiento radiactivo total para el 2100 que oscilará entre 2.6 y 8.5 W.m^-2 y será mayor que en la actualidad). Las RCP comprenden un escenario en el que los esfuerzos en mitigación conducen a un nivel de forzamiento muy bajo (RCP2.6), dos escenarios de estabilización (RCP4.5 y RCP6.0) y un escenario con un nivel muy alto de emisiones de GEI (RCP8.5). Las proyecciones de temperatura y precipitación en estos RCP han sido publicadas en el Worldclim (2016) útiles para el mapeo y modelado espacial.

Entonces, siendo el departamento de Tacna de poca superficie territorial, en el cual se presentan variaciones de temperaturas en diversas zonas, p. ej. La Yarada, (Senamhi, 2010) con variaciones extremas máximas y mínimas de entre 1 °C y 1.3 °C, es que se le ha considerado como uno de los primeros departamentos en donde ya se evidencian los efectos del cambio climático. Por lo anterior, el gobierno regional de Tacna (Gore Tacna), elaboró su Estrategia regional de cambio climático (ERCC) aprobada mediante Ordenanza Regional N.º 005-2015-CR/GOB. REG. TACNA (2015).

En este contexto es válido preguntarse ¿cuál será la distribución geográfica actual de los nichos ecológicos de las especies de flora silvestre categorizada CR del departamento de Tacna y como cambiarán frente a las variaciones climáticas en las cuatro RCP del año 2050 y 2070?

Como respuesta se realizó el modelamiento actual y proyectado de los nichos ecológicos de las especies categorizadas CR utilizando el software MaxEnt, los datos de ocurrencia de campo y del estudio Zonificación ecológica y económica (ZEE) Tacna (Ordenanza Regional N.º 016-2012-CR/GOB. REG. TACNA, 2013) y un conjunto de variables ambientales (bioclimáticas, altitud, radiación solar, etc.) actuales y proyectadas a los RCP futuros. La información científica obtenida servirá para que se diseñen, a nivel departamental, adecuados planes de gestión territorial con eficientes medidas de adaptación y mitigación al cambio climático, que estén estrechamente vinculadas con la conservación.

MATERIALES Y MÉTODOS

Área de estudio

El departamento de Tacna está localizado en el extremo sur del Perú, limita por el noroeste con Moquegua, por el noreste con Puno, por el sur con la República de Chile, por el este con la República de Bolivia y por el oeste con el océano Pacífico (figura 1). Tiene una superficie territorial de 16 075.73 km² que representa el 1.25 % del país y está dividida en cuatro provincias: Tacna, Tarata, Candarave y Jorge Basadre (Gobierno Regional de Tacna, 2014).

Figura 1 Mapa de localización estratégica del departamento de Tacna en Sudamérica y de distribución de las coordenadas geográficas en las que fueron reportadas las especies de flora categorizada (CR) evaluadas al interior del área de estudio 

Especies de flora silvestre categorizada evaluadas

De acuerdo con la lista oficial de especies de flora categorizada del Perú (Decreto Supremo N.º043-2006-AG, 2006), en el departamento de Tacna se encuentran dentro de la categoría en peligro crítico (CR) a Buddleja coriacea J. Remy, Carica candicans A. Gray, Haplorhus peruviana Engl, Kageneckia lanceolata Ruiz & Pav y Weberbauerella brongniartioides Ulbr. La cantidad de ocurrencias o coordenadas de ubicación geográfica para cada especie categorizada CR que se reportaron en la ZEE Tacna se presenta en la tabla 1, y estas fueron contrastadas con información técnica no publicada del Herbario Takana, informes técnicos de estudios poblacionales de biodiversidad realizados por el Gore Tacna, así como el estudio poblacional de C. candicans realizado por ^{Franco (2013)}. Del mismo modo, se realizaron 10 salidas de campo con 60 puntos de búsqueda, observación y fotografiado de las especies (sin colecta) para confirmar su presencia/ausencia en las coordenadas referidas, las cuales se presentan en la figura 1. Respecto a la importancia de las especies categorizadas evaluadas, B. coriacea y K. lanceolata son especies forestales proveedoras de madera de excelente calidad, utilizadas en construcción de vigas, puertas, ventanas, dinteles y herramientas agrícolas (yuntas, arados, tacllas) así como artesanías y utilería; igualmente, es común verlo en compuertas y partes de los canales de regadío, además su leña y carbón son de excelente calidad (^{Reynel y León, 1990}; Reynel y Marcelo, 2009). Asimismo, C. candicans es una planta frondosa y foliosa y juega un papel importante como captadores naturales de niebla, son productores primarios en conjunto con el estrato arbóreo arbustivo y sirven como hábitat o nicho ecológico de organismos que cumplen similares funciones dentro del ecosistema (^{Mendoza, 2006}; ^{Sagástegui, Rodríguez y Arroyo, 2007}; Weberbauer 1945). Respecto a H. peruviana, Pumasupa (2017) determinó un stock de carbono de 4.5244 tC/ha, demostrando así su importancia como proveedores de servicios ecosistémicos de regulación, mientras que W. brongniartioides es endémica para las costas del sur del Perú del que aún no se tienen muchos datos de utilidad de la especie.

Tabla 1 Nombre científico, nombre común, usos e importancia y número de ocurrencias o coordenadas geográficas de reporte de cada una de las cinco especies de flora categorizada (CR) evaluadas en el modelamiento de nicho ecológico actual y proyectado 

RESULTADOS

Todos los modelos de nicho ecológico obtenidos para las cinco especies de flora categorizada CR elaborados en MaxEnt dieron como resultado valores estadísticos AUC mayores a 0.9, siendo el menor valor 0.905 y el mayor 0.999, estos se consideran modelos predictivos de alta confianza.

Se obtuvieron nueve modelos de nicho ecológico para cada especie (el modelo actual más cuatro proyectados al año 2050 y para el 2070), conformando así un total 45 archivos raster del tipo ASCI que contenían la distribución geográfica de los nichos ecológicos de las cinco especies, tanto para la actualidad como para los cuatro RCP proyectados. En la tabla 2 se presentan los resultados de cobertura actual del nicho ecológico y el porcentaje de pérdida o ganancia de cobertura en el RCP proyectado. Allí se observa que solo la especie B. coriacea tendrá pérdidas de superficie de nicho ecológico en todos los RCP proyectados al 2050 y 2070; mientras que C. candicans y W. brongniartioides presentan ganancia de superficie en todos los escenarios evaluados. H. peruviana demuestra pérdida de superficie de su nicho ecológico en los tres primeros escenarios del 2050, pero a partir del RCP8.5 estos tienden a volverse positivos (ganancia), mientras que K. lanceolata también presenta ganancia en todas las proyecciones RCP. En la misma tabla 2 se presentan los resultados de la prueba de Jakknife, resaltándose a aquella variable que aporto mayor información al modelamiento de los nichos ecológicos por especie. Se evidencia que predominan las variables bioclimáticas, en especial las de precipitación.

Tabla 2 Resultados obtenidos en la predicción de los nichos ecológico para las especies de flora categorizada (CR). La fila km² representa la superficie actual de cada especie. Las filas % representan la ganancia o pérdida (según signo) de superficie de los nichos ecológicos de las especies en los RCP del año 2050 y 2070. Las filas Var indica el nombre de las variables que contribuyeron con mayor información en el diseño del modelo de nicho ecológico, según la prueba de Jakknife realizado con el MaxEnt 

DISCUSIÓN

Los resultados de AUC obtenidos durante el modelamiento de nicho ecológico de las cinco especies categorizadas CR del departamento de Tacna (mayores a 0.9) demuestran alta confiabilidad (^{Phillips et al., 2006}) como los resultados de ^{Reynoso et al. (2018)} de 0.882 para Pinus oocarpa y 0.947 para Pinus pseudostrobus; o como los de ^{Wan et al. (2015)} entre 0.7 y 0.9 para Pinus koraiensis y Eleutherococcus senticosus; o los de ^{Ibarra et al. (2016)} para Coryphantha chihuahuensis de 0.9. Todo esto comprueba que las variables ambientales y el software MaxEnt son las principales herramientas para predicciones de nicho ecológico (^{Elith et al., 2006}) demostrado por ^{Ortega y Townsend (2008)} al compararlo entre seis diferentes métodos (algoritmos) resultando predicciones más significativas, sin fallas graves y con el análisis de participación de cada variable.

En general, el modelamiento de nichos ecológicos viene siendo aplicado a diferentes especies silvestres (^{Morales, 2012}), ya sea para predecir su distribución o para gestionar oportunidades de aprovechamiento; como la publicación de Osinfor (2013) que presenta modelos de nicho ecológico para 29 especies forestales utilizando el software MaxEnt. También ha sido utilizado para determinar zonas potenciales para especies cultivables de flora con elevado valor comercial, como la Vanilla planifolia (^{Hernández et al., 2016}) a manera de oportunidad de econegocio. Sin embargo, es necesario considerar el cuestionamiento a la utilidad de los análisis ROC al ponderar igual los errores de omisión que los de comisión (^{Lobo et al., 2007}; ^{Peterson, Papes y Soberón, 2008}) por lo que ^{Martínez, Aguirre, Eguiarte y Jaramillo (2016)} utilizaron el programa Tool for Partial-ROC (^{Narayani, 2008}) para incrementar el nivel de confianza de los modelos motivo de su investigación.

De esta manera fue que se calculó el nicho ecológico actual de la especie forestal B. coriacea cuya superficie es de 540.98 km² como se muestra en la tabla 2, pero esta tenderá a disminuir en las proyecciones RCP8.5 del año 2050, con una pérdida de 319.04 km² (-58.97 %) y de 370.80 km² (-68.54 %) para el año 2070, como se observa en la figura 2. Pero, la pérdida más significativa se proyecta para el RCP6.0 del año 2050, con un -86.19 %, lo que podría considerarse como una futura extinción regional de la especie por cambio climático (figura 3), lo que difiere de lo que ^{Quesada, Acosta, Arias y Rodríguez (2017)} obtuvieron al modelar nichos ecológicos para cinco especies de plantas en zonas altas de Costa Rica (1500 a 3820 metros de altitud), en el que el valor más alto de pérdida de superficie de nicho ecológico en el escenario más pesimista (RCP8.5) sería de -69 %. Está pérdida de nicho ecológico de B. coriacea podría deberse a la distribución altitudinal natural de la especie, entre los 3400 a 4500 metros (^{Reynel y Marcelo, 2009}), ya que en los gradientes altitudinales el clima en combinación con otros factores ambientales ejercen influencia en la sobrevivencia de la flora a variaciones extremas de temperatura y precipitación (Guariguata y Matan, 2003). Esto permite explicar patrones de distribución de plantas (Wood Ward, 1987) al registrar cambios en el tamaño y localización de las zonas de idoneidad ambiental para muchas especies en el mundo con cambios más severos en zonas montañosas, donde muchas especies tienden a elevar su límite altitudinal (Wilson et al., 2005).

Figura 2 Gráficos de barras que muestran las tendencias en los cambios de la distribución actual de los nichos ecológicos de las especies categorizadas CR del departamento de Tacna en los escenarios RCP2.6, RCP4.5, RCP6.0 y RCP8.5 proyectados al 2050 y 2070 con MaxEnt. 

Igualmente, se determinó un nicho ecológico actual de 1135.55 km² para la especie K. lanceolata, especie forestal que ^{Brako y Zarucchi (1993)} describen su presencia en zonas con altitudes entre los 2000 y 4100 metros, mientras que ^{Montesinos, Sýkora, Quipuscoa y Cleef (2015)} determinaron su distribución en el departamento de Arequipa entre los 2800 y 3040 metros de altitud. Es posible que este amplio rango altitudinal de desarrollo natural de la especie haga que se adapte rápidamente a diferentes condiciones ambientales y, por ello, su nicho ecológico tienda a ganar superficie a pesar de los cambios ambientales futuros, como se observa en la figura 3. De manera similar a la especie Ilex pallida en la investigación de ^{Quesada et al. (2017)}, la cual tiene mayor distribución en el gradiente altitudinal de las cinco especies seleccionadas y que para los escenarios RCP 2.6 y RCP 4.5 el área aumenta con respecto al nicho ecológico actual, asegurando nuevos terrenos de zonas altas con características similares a las actuales. Esto hace suponer que la categoría de amenaza CR de K. lanceolata se deba principalmente al desarrollo de actividades antropogénicas como la agricultura, más que en los efectos del cambio climático.

Figura 3 Distribución potencial de los nichos ecológicos actuales y proyectadas para los años 2050 y 2070, en cuatro escenarios de cambio climático: RCP2.6, RCP4.5, RCP6.0 y RCP8.5 de las especies a) B. coriacea b) C. candicans c) H. peruviana d) K. lanceolata e) W. brongniartioides. 

Con respecto a los nichos ecológicos obtenidos para las especies de flora silvestre categorizada que se distribuyen a menor altitud (de 0 a 2000 metros) que habitan zonas desérticas, costeras y de lomas, como C. candicans (75.47 km²), H. peruviana (254.05 km²) y W. brongniartioides (185.94 km²), sus superficies tienden a mantenerse e incrementarse incluso en los RCP proyectados al año 2050 y 2070 como se muestra en la figura 3. En estas se observa que no se verán afectadas por las nuevas condiciones climáticas de los RCP futuros, lo que estaría asociada a sus ciclos vegetativos. Aunque ^{Leiva, Galloso y Chang (2016)} describieron a C. candicans como especie psicrófila (crece en lugares fríos), heliófila (requiere abundante luz del sol), higrófila (prefiere suelos húmedos, tierras negras con abundante humus), psamófila (puede desarrollarse en suelos arenosos, pedregosos) y argilícola (prospera en suelos arcillosos); aunque, por lo general, es una especie eutrofa que prefiere suelos con abundantes nutrientes. Esto también hace suponer que su categoría de amenaza CR se debe principalmente a la frecuente actividad antropogénica establecida en esta parte del país que el PNUD (2018) identifica como minería desregulada, invasiones de terrenos para viviendas, el tráfico de terrenos y el pastoreo no regulado de ganado.

Es así como al utilizar las mismas 19 variables bioclimáticas y otras de orientación, pendiente e índice topográfico de humedad y con la misma resolución espacial, pero para el área correspondiente a México, Estados Unidos y Centroamérica. ^{Martínez et al. (2016)} modelaron el nicho ecológico de las especies del género Abies (pinaceae) en México, obteniendo que las variables que más información aportaron en la mayor cantidad de especies fueron la Bio_19 (precipitación del cuarto más frío), común para seis especies, seguida de elevación, Bio_2 (intervalo promedio de temperaturas diurnas), Bio_18 (precipitación del cuarto más cálido) y pendiente, las cuales fueron compartidas por cinco especies. De igual manera, con las mismas variables bioclimáticas ^{Quesada et al. (2017)} concluyeron que las variables más determinantes para el modelamiento de los nichos ecológicos de las especies utilizadas en su investigación son dependientes de las variables de temperatura, mientras que las de precipitación tienen aportes mínimos. En el caso del nicho ecológico actual de B. coriacea la variable predominante fue la Bio_1 (temperatura media anual) y las proyecciones de los RCP por la Bio_14 (precipitación del mes más seco) a excepción del RCP85 del año 2070, donde vuelve a influir la Bio_1. Para C. candicans todos los modelos obtenidos fueron influenciados por la Bio_12 (precipitación anual). Para H. peruviana la variable más influyente fue la Bio_7 (rango anual de temperatura), para K. lanceolata la Bio_18 (precipitación de cuarto mes más caluroso), al igual que para W. brongniartioides. Algunas otras variables de temperatura como Bio_9 (promedio en el trimestre más seco) obtuvieron valores altos solo para dos especies: C. candicans y H. peruviana.

Asimismo, en el Perú Osinfor (2016) realizó el modelamiento de la distribución potencial de nicho ecológico para 18 especies forestales en el departamento de Loreto, utilizando las mismas variables bioclimáticas y MaxEnt. En esta se concluye que la variable Bio_12 (precipitación total anual), es la que aporta mayor información a todos los modelos. Esta misma variable influyó principalmente en el modelo de C. candicans (tabla 2) y significa que para el ciclo vegetativo de estas especies es más importante el volumen anual acumulado de precipitación y no tanto las variaciones que por factores externos podrían alterar el régimen de lluvias en determinados meses del año.

Igualmente, durante el ciclo vegetativo de las cactáceas existe un predominio de la temperatura sobre la precipitación, por lo que se espera que en el modelamiento de sus nichos ecológicos estas variables influyan prioritariamente, como quedó evidenciado en el trabajo de ^{Ibarra et al. (2016)} en el que las variables bioclimáticas de temperatura fueron importantes en la distribución y persistencia de sus poblaciones y por sí solas fueron capaces de delimitar la distribución actual de Coryphantha chihuahuensis. Este mismo patrón sucedió en la especie forestal H. peruviana, cuyos nichos ecológicos derivaron exclusivamente de variables de temperatura, por lo que a pesar de no ser una cactácea, demuestra un similar ciclo vegetativo.

Con respecto a las variables de radiación solar, estas influyeron solo en el modelo actual del nicho ecológico de la especie desértica W. brongniartioides, pero no tuvieron mayor influencia en otras especies ni en los RCP proyectados. La variable ELE que representa a los metros de altitud del terreno tuvo valores medios de contribución en todos los modelos, pero sin ser predominantes.

Estos resultados permitieron conocer y analizar la distribución actual y potencial de los nichos ecológicos de las especies categorizadas CR de Tacna en los escenarios de cambio climático (IPCC, 2007), demostrándose que la altitud de distribución de las especies será influyente en la reducción del nicho ecológico frente a las variaciones de precipitación y temperatura proyectados para los RCP de los años 2050 y 2070, como se observa en la figura 3; mientras que habrá incremento del nicho ecológico de las especies de mayor rango altitudinal y que habitan zonas costeras. Sin embargo, considerando que los cambios de temperatura y precipitación en Tacna son muy variables en distancias muy cortas (Senamhi, 2010), existe una mayor probabilidad de error en las bases de datos de Worldclim en zonas de mayor altitud. Pero esto no significa que la información obtenida pueda ser útil para que las autoridades competentes encaren el cambio climático con mejores herramientas y se trabaje en lograr lo que ^{Vegas (2010)} señala en cuanto que las variaciones climáticas son una oportunidad para movilizar a la sociedad a enfrentar los problemas del subdesarrollo, capitalizar la economía, superar la pobreza y la desigualdad social en el ámbito local, regional y nacional, en especial con las especies nativas y endémicas.

En el departamento de Tacna ^{Franco, Cáceres y Sulca (2019)} determinaron 708 especies de flora vascular, siendo 92 endémicas. Los endemismos peruanos están principalmente vinculados a los Andes (Young, Ulloa, Luteyn y Knapp, 2002), destacando a B. coriacea como la especie más andina de la presente investigación. ^{León, Pitman y Roque (2013)} señalan que a través del conocimiento que se tenga sobre las especies endémicas de flora será posible discutir sobre el estado de conservación y el real futuro de la diversidad biológica, señalando que cuatro de las cinco especies categorizadas evaluadas en la presente investigación son endémicas (C. candicans, H. peruviana, B. coriacea y W. brongniartioides) mientras que se comparte la distribución de K. lanceolata con otros países de Sudamérica. Sin duda, el conocimiento generado sobre sus nichos ecológicos actuales y futuros de las especies endémicas evaluadas es un gran aporte para la conservación de la diversidad biológica. Por ello, se espera que a partir de los resultados obtenidos se empiece a discutir y promover el establecimiento estratégico de áreas naturales protegidas de administración pública o privadas en la superficie de los nichos ecológicos modelados, con fines de aprovechamiento, ecoturismo y conservación, como lo mencionan ^{Ranganathan y Daily (2008)}, y con ella mejorar las oportunidades de desarrollo económico mediante el manejo de especies forestales o arbustivas en zonas de distribución silvestre o en plantaciones forestales. Estas pueden ser instaladas en zonas eriazas o agrícolas con las condiciones adecuadas a los nichos ecológicos de las especies evaluadas, a manera de recambio, como estrategia para la adaptación ante el cambio climático. Puesto que el IPCC (2013) señala que las tierras agrícolas serán afectadas a nivel productivo y, por ende, la seguridad alimentaria y el desarrollo económico deberán orientarse hacia la nueva estrategia denominada “Adaptación basada en ecosistemas” que propone la UICN en la publicación de ^{Lhumeau y Cordero (2012)}).

Es por esto que las plantaciones forestales, en especial aquellas que trabajan con especies nativas, son una importante herramienta de mitigación y adaptación al cambio climático (^{Montesinos et al., 2019}). De esta manera, es necesario lograr eficiencia en los proyectos de forestación y reforestación, por lo que es fundamental conocer las condiciones ambientales donde naturalmente habitan y que ^{Mindreau y Zúñiga (2010)} señalan como la exposición del terreno, humedad, altitud, pendiente, suelos, climas, flora acompañante y de conectividad; insumos que lograron obtenerse de las especies forestales evaluadas y que también pueden representar áreas potenciales para forestación.

El no haber considerado durante el modelamiento de nichos ecológicos variables de cambio de uso de suelo, es decir, el remplazo del nicho ecológico actual con infraestructura de la actividad antropogénica presente, ha hecho que el MaxEnt calcule superficies entre 254.05 a 1135.55 km² (a excepción de C. candicans cuya superficie resultó en 75.47 km²), como se observa en la tabla 2. Lo anterior difiere del concepto de la categoría de amenaza de la especie “En peligro crítico”, que establece un límite de superficie de distribución de 100 km² (Decreto Supremo N.º 043-2006-AG, 2006), lo que fue realizado por ^{Kurpis, Serrato y Arroyo (2019)}, quienes usando un mapa de usos de la tierra concluyeron que las capacidades de dispersión de Tagetes lucida están mayormente amenazados por actividad humana.

Lo relevante de la presente investigación, y de los resultados obtenidos, es que se enmarcan en lo señalado por la Secretaría del Convenio de Diversidad Biológica (2004), la cual indica que las fuerzas motrices de los ecosistemas varían en el espacio y el tiempo, por lo que deben gestionarse y estudiarse a más de una escala. Esto evidencia la urgencia por analizar el territorio en los diferentes escenarios climáticos para predecir el futuro de las especies de flora y fauna que estarían en riesgo de extinción; como lo realizado por ^{Ortíz, Restrepo y Páez (2014)} al modelar la posible fluctuación de la distribución potencial de la tortuga Podocnemis lewyana Duméril en escenarios de cambio climático en Colombia y ^{Pliscoff y Fuentes (2011)} al realizar una revisión de las nuevas herramientas y enfoques empleados al modelar la distribución de las especies y ecosistemas. Asimismo, se produce información científica de calidad sobre la distribución actual y futura de los nichos ecológicos de especies de flora con alto grado de amenaza y pocos antecedentes en investigación, que resulta pertinente y ayudará a las autoridades a planificar el territorio de acuerdo con los intereses globales de conservación, protección y manejo, orientado al desarrollo sostenible.

CONCLUSIONES

Los nichos ecológicos modelados para las cinco especies de flora silvestre categorizada del departamento Tacna tienen alto nivel de confianza. De igual manera, las variables ambientales más determinantes para cuatro de las cinco especies en la elaboración de sus modelos actuales y proyectados (2050 y 2070) fueron las que derivan de la precipitación, mientras que las variables de temperatura fueron altamente influyentes solo en el modelamiento del nicho ecológico de H. peruviana.

Es por ello que las políticas, actividades e iniciativas de conservación y ordenamiento territorial que se deriven de la información científica obtenida deberán ir orientadas a lograr resultados concretos en el escenario de mitigación RCP2.6 para garantizar el menor impacto del cambio climático a las especies. Siendo el nicho ecológico de B. coriacea el más variable en el tiempo, es necesario desarrollar programas de forestación y reforestación basadas en investigaciones científicas similares y aplicadas socioeconómicamente bajo la estrategia de adaptación basada en ecosistemas propuesta por la UICN mundialmente. Asimismo, para las especies evaluadas que habitan lomas costeras y desérticas, cuyo desarrollo futuro no se ve influenciada por las RCP, será necesario complementar este estudio con información de las actividades antropogénica que se desarrollan y que aproximarían aún más la distribución actual y proyectada de los nichos ecológicos con la finalidad de establecer mecanismos de ordenamiento territorial y de modalidades de conservación y reforestación.

Finalmente, la generación de escenarios como los calculados en el presente trabajo de investigación nos permite mejorar significativamente la comprensión de la dinámica de los nichos ecológicos de las especies de flora categorizada que habitan ecosistemas costeros, andinos y altoandinos a través del monitoreo de la biodiversidad, mostrándonos una alternativa más para la adecuada gestión y planificación territorial.