INSTITUTO DE INVESTIGACIONES DE LA AMAZONIA PERUANA

PROGRAMA DE CAMBIO CLIMÁTICO, DESARROLLO TERRITORIAL Y AMBIENTE

(PROTERRA)

PROGRAMA PRESUPUESTAL 0068: REDUCCIÓN DE LA

VULNERABILIDAD Y ATENCIÓN DE EMERGENCIAS POR
DESASTRES
ACTIVIDAD DEL PP0068 (5005572): DESARROLLO DE INVESTIGACIÓN
APLICADA PARA LA GESTIÓN DEL RIESGO DE DESASTRES

ACTIVIDAD OPERATIVA:
EVALUACIÓN DE LA DINÁMICA DEL CAMBIO DE LA COBERTURA Y USO DE
LA TIERRA
 Instituto de Investigaciones de la Amazonía Peruana

INSTITUTO DE INVESTIGACIONES DE LA AMAZONIA PERUANA

PROGRAMA DE CAMBIO CLIMÁTICO, DESARROLLO TERRITORIAL Y AMBIENTE (PROTERRA)

RESULTADOS DE LA ACTIVIDAD OPERATIVA: EVALUACIÓN DEL CAMBIO DE COBERTURA Y

USO DE LA TIERRA, EN EL ÁREA DE INFLUENCIA DE LA PROPUESTA DE CARRETERA
BALSAPUERTO - MOYOBAMBA

Elaboración de contenidos : Programa de Cambio Climático, Desarrollo Territorial y

Ambiente (PROTERRA)

Coordinación : Lizardo Fachín Malaverri

Elaborado por : Juan Palacios Vega
Colaboración : Rocío Jarama Vilcarromero, Ricardo Zárate Gómez, Lizardo
Fachin Malaverri, Guiuseppe Torres Reyna, María Marín
Pérez, Joel Lam Rodríguez

Este libro fue revisado por pares.

Diagramación : Angel G. Pinedo Flor

Hecho el Depósito Legal en la Biblioteca Nacional del Perú Nº 2019-00000

Primera edición
Iquitos, Perú. Diciembre de 2019

© 2019 Instituto de Investigaciones de la Amazonía Peruana

Av. José Abelardo Quiñones km. 2.5
Teléfonos: (+51) (65) 265515 / 265516 Fax: (+51) (65) 265527
Iquitos, Perú, 2019
www.iiap.gob.pe / E-mail: [email protected]

El presente estudio fue financiado con fondos del Instituto de Investigaciones de la

Amazonía Peruana – IIAP, Programa Presupuestal 0068.

Cita sugerida:
Palacios, J., Jarama, A., Zárate, R., Fachín, L., Torres, G., Marín, M., Lam, J. 2019. Informe de
investigación. Evaluación del cambio de cobertura y uso de la tierra, en el área de influencia de la
propuesta de carretera Balsapuerto - Moyobamba.

La información contenida en este informe puede ser reproducida total o parcialmente

siempre y cuando se mencione la fuente de origen.

4
 Evaluación del cambio de cobertura y uso de la tierra, en el área de influencia de la
propuesta de carretera Balsapuerto - Moyobamba

CONTENIDO

I. INTRODUCCION ....................................................................................................................... 6
II. OBJETIVOS ............................................................................................................................... 7
2.1. OBJETIVO GENERAL............................................................................................................... 7
2.2. OBJETIVOS ESPECÍFICOS ........................................................................................................ 7
III. MATERIALES Y METODOS......................................................................................................... 7
3.1. MATERIALES .......................................................................................................................... 7
3.2. METODOLOGÍA ................................................................................................................... 11
IV. RESULTADOS .......................................................................................................................... 19
V. BIBLIOGRAFIA ........................................................................................................................ 34

Lista de Figuras
Figura 1. Interface de usuario del servicio web Alaska Facility, para la descarga de imágenes DEM Alos Palsar. ....... 8
Figura 2. Interface de usuario del servicio web Copernicus ESA, para la descarga de imágenes Sentinel 2. ............... 9
Figura 3. Mapa de Ubicación del área de influencia de la propuesta de carretera Balsapuerto - Moyobamba. ....... 12
Figura 4. Mapa de Uso Actual de la Tierra en el área de influencia de la propuesta de carretera
Balsapuerto - Moyobamba. ....................................................................................................................................... 21
Figura 5. Cultivo de arroz en el sector Moyobamba, San Martín. ............................................................................. 22
Figura 6. Pastos en el sector Balsapuerto, Loreto. .................................................................................................... 23
Figura 7. Cultivos de piña, papaya, plátano y café en los sectores Balsapuerto y Moyobamba. ............................... 25
Figura 8. Bosque denso alto en el sector Moyobamba.............................................................................................. 27
Figura 9. Bosque denso bajo de tierra firme en el sector Moyobamba. .................................................................... 29
Figura 10. Vegetación de Bosque denso bajo inundable en el sector Moyobamba. ................................................. 30
Figura 11. Variación de las superficies de Cobertura y Uso de la Tierra (CUT) por año, datos de CUT 2015
al 2019 y 2023 datos simulados. ............................................................................................................................... 31
Figura 12. Comparación de mapas de cobertura y uso de la tierra años 2015, 2019 y proyectado al año 2023. ...... 33

Lista de Tablas
Tabla 1. Propuesta de Leyenda de Uso y Cobertura de la Tierra - Perú (MINAM, 2012) y adaptada para el
estudio en el área de influencia de la propuesta de carretera Balsapuerto - Moyobamba. ...................................... 19
Tabla 2. Categorías por niveles de Cobertura y Uso de la Tierra adaptada al estudio del área de influencia de la
propuesta de carretera Balsapuerto – Moyobamba, leyenda de Cobertura y Uso de la Tierra - Perú
(MINAM, 2012). ........................................................................................................................................................ 20

5
 Instituto de Investigaciones de la Amazonía Peruana

I. INTRODUCCIÓN

En los últimos 50 años la preocupación sobre el propone entonces realizar reconstrucciones

cambio de cobertura y uso de la tierra (CCUT) a históricas de estos procesos y proyectar al
nivel global, regional y local ha contribuido a futuro los cambios en la cobertura y uso de la
generar mucha investigación sobre el tema, tierra.
llegando ésta a concluir, que son muchos los
impactos de estos procesos los mismos que En el presente estudio se evalúa la dinámica del
afectan al clima, a los sistemas biológicos y al cambio de cobertura y uso de la tierra en el
ser humano en particular. Son precisamente área de influencia de la propuesta de carreta
éstos los más vulnerables a estos cambios así Balsapuerto - Moyobamba, para lo cual se
como también los lugares donde se asientan, lo recopilo información histórica de la Cobertura y
que afecta las economías y actividades Uso de la Tierra (CUT) en el área de estudio, se
productivas, sociales, entre otras. generó información de CUT al año 2015
utilizando imágenes de satélite, seguidamente
En la Amazonía peruana la ocupación del se cuantifico e identifico los principales
territorio y el uso de los recursos naturales han cambios, a esto se sumó variables espaciales
generado conflictos ambientales por el mal uso explicativas con las cuales se ingresó a modelar
de la tierra, tierras con vocación forestal o de la transición potencial y finalmente simular un
protección son utilizados con fines escenario al año 2023 de CCUT.
agropecuarios, trayendo como consecuencia
pérdida de la fertilidad del suelo y problemas La importancia del presente trabajo radica en
de erosión principalmente en sitios de mayor obtener un documento que sirva en la gestión
pendiente, lo cual se manifiesta en el bajo para el ordenamiento de los usos del territorio
rendimiento de diversos cultivos y por del área estudio, es por eso que el objetivo de
consiguiente el abandono de las tierras. este estudio es analizar la dinámica del CCUT en
el periodo 2015-2019, en el área de influencia
En este contexto la evaluación de la dinámica de la propuesta de carreta Balsapuerto -
de cambio de la cobertura y uso de la tierra se Moyobamba, con la finalidad de simular un
hace necesaria, tanto para cuantificar el cambio escenario al año 2023 de CCUT que explique e
en el tiempo, así como, para comprender y identifique los principales cambios que se
predecir el impacto de estos procesos en el pudieran dar a futuro sobre el área de estudio.
territorio con enfoque de cambio climático. Se

6
 Evaluación del cambio de cobertura y uso de la tierra, en el área de influencia de la
propuesta de carretera Balsapuerto - Moyobamba

II. OBJETIVOS

2.1. OBJETIVO GENERAL

Análisis de la dinámica del cambio de cobertura y uso de la tierra en el periodo 2015-2019, en el área
de influencia de la propuesta de Carretera Balsapuerto – Moyobamba.

2.2. OBJETIVOS ESPECÍFICOS

 Recopilación de datos e información de cobertura y uso de la tierra del área del estudio
 Caracterización de la cobertura y uso actual de la tierra del área de estudio.
 Análisis del cambio de la cobertura y uso de la tierra en el periodo 2015 – 2019 del área de
estudio
 Simulación del cambio de cobertura y uso de la tierra al año 2023 del área de estudio.

III. MATERIALES Y MÉTODOS

generadoras de datos espaciales y en otros

3.1. MATERIALES
casos el mismo IIAP posee una base de datos de
las ZEE u otros proyectos ejecutados en el área
Para el presente estudio los materiales
de estudio.
utilizados se describen como aquella
información utilizada y recopilada de las
A continuación se menciona los materiales de
fuentes o servicios oficiales encontrados libres
tipo vectorial utilizados en el desarrollo del
en la Web tanto nacionales como
proyecto:
internacionales, las cuales proveen en su
 Centros poblados, obtenidos de la data
mayoría de datos e información satelital y
del INEI del Censo del año 2017.
vectoriales (información geoespacial), esta se
 Zonificación Económica Ecológica de Alto
encuentra libre para la descarga, ordenada por
Amazonas, realizados por el Instituto de
ubicación geográfica y clasificada por día, mes y
Investigaciones de la Amazonía Peruana y
año de toma.
el Gobierno Regional de Loreto.
 Estructuras de Datos Espaciales,
3.1.1. Materiales tipo vectorial obtenidos del Gobierno Regional de San
Los materiales de tipo vectorial se refieren Martín.
específicamente a las capas de información  Red hídrica nacional, obtenidos del
geoespacial en formato estándar Shapefile Instituto Geográfico Nacional.
(shp), estos por lo general puedan ser  Carta Nacional, obtenida del Instituto
visualizados, procesados y representados en Geográfico Nacional.
plataformas o programas de tipo SIG, estos por  Red vial nacional, obtenida del Ministerio
lo general son recolectados de instituciones de Transportes y Comunicaciones.

7
 Instituto de Investigaciones de la Amazonía Peruana

 Límites departamentales, provinciales y A continuación se menciona los materiales de

distritales, obtenidos del Gobierno tipo raster utilizados en el desarrollo del
Regional de Loreto y el Gobierno Regional estudio:
de San Martín.
a. Modelo de Elevación Digital – DEM ALOS
3.1.2. Materiales tipo raster o satelital PALSAR, obtenido del portal Vertex de la
NASA
Los datos e información en formato raster
definen el mundo real por celdas ordenadas en Existen varios sistemas de teledetección
filas y columnas llamadas pixel, en las que a espacial que permiten obtener Modelos
cada una de ellas se le asigna un valor con el Digitales de Elevación (MDE) de alta resolución,
cual se puede clasificar para determinar capas entre ellos se destaca el sensor PRISM
de información temática. Este tipo de dato por (Panchromatic Remote-sensing Instrument for
lo general es recolectado de fuentes o Stereo Mapping) a bordo de la plataforma
servidores gratuitos en el internet, en otros satelital ALOS (Advanced Land Observing
casos se adquieren a proveedores. Satellite), que fue diseñado para obtener MDE
de diez metros de resolución espacial y RMS
menor a cinco metros (Díaz et al. 2010).

Figura 1. Interface de usuario del servicio web Alaska Facility, para la descarga de imágenes DEM Alos
Palsar.

El satélite ALOS durante su operación (mayo 2006 - abril 2011), colectó imágenes de Radar en escenas
de 50 km x 70 km de todo el planeta cada 45 días aproximadamente a través de la banda L de su
sensor PALSAR (Phased Array Type L-band Synthetic Aperture Radar). PALSAR es uno de los tres
instrumentos avanzados para observación terrestre, desarrollados para contribuir a los campos de la
cartografía, la observación de la tierra, la cobertura regional precisa, el seguimiento de desastres y
levantamiento de recursos. (Figura 1).

8
 Evaluación del cambio de cobertura y uso de la tierra, en el área de influencia de la
propuesta de carretera Balsapuerto - Moyobamba

b. Imágenes Sentinel 2, del Programa Copernicus de la Comisión Europea

El servidor SCHIHUB COPERNICUS, https://scihub.copernicus.eu/ servidor del cual se descargaron las 2
escenas del satélite SENTINEL 2 para el desarrollo del proyecto. (Figura 2).

Figura 2. Interface de usuario del servicio web Copernicus ESA, para la descarga de imágenes Sentinel
2.

c. Proyecto “Herramientas de mapeo reducción de la pobreza y la deforestación neta

remoto y valoración socioeconómica de la selva primaria en Perú.
para apoyar la planificación e
implementación de intervenciones de 3.1.3. Materiales de gestión para el
uso de la tierra” acopio de datos
Estos se refieren a los documentos y/o
Consiste en una capa de información raster mecanismos que se utilizan en la gestión y
correspondiente al mapeo del Inventario del solicitud a las fuentes o entidades tanto
Uso de la Tierra realizado por Remote Sensing públicas como privadas; nacionales, regionales
Applications Consultants Ltd (RSAC) en el marco o locales que puedan proveer de datos e
del proyecto “Herramientas de mapeo remoto información para la ejecución del proyecto.
y valoración socioeconómica para apoyar la
planificación e implementación de También se hace mención que Equipo
intervenciones de uso de la tierra”, como parte Funcional de Territorios Amazónicos (TERRA)
del Programa de Asociación Internacional (IPP) del Instituto de Investigaciones de la Amazonia
de la Agencia Espacial del Reino Unido (UKSA). Peruana (IIAP) es una fuente importante de
información tanto cartográfica como satelital
El Inventario de Uso de la Tierra es el principal debido al tiempo en que se vienen ejecutando
insumo del sistema IMAGE (Inventaire- proyectos para el ordenamiento del territorio
Modelisation-Alerte-Gestion-Environmental) en la Amazonia peruana y por consecuencia la
desarrollado por Vivid Economics Ltd para este generación de una cantidad considerable de
proyecto. En general, el proyecto busca información secundaria.
desarrollar herramientas para apoyar la

9
 Instituto de Investigaciones de la Amazonía Peruana

A continuación se enumera a las entidades a 30 metros y pueden ser utilizadas en diversas

que generan y gestionan información aplicaciones de la superficie del terreno.
geoespacial oficial, la cual se fue solicitada y
utilizada en el presente estudio:
3.1.5. Equipos
 Gerencia Regional de Planeamiento,
Presupuesto y Acondicionamiento Se enumera a continuación la lista de equipos
Territorial del Gobierno Regional de utilizados para realizar el ingreso y gestión de
Loreto (GOREL), a quienes se les solicito datos e información geoespacial de los
información cartográfica y satelital. proyectos:
 Municipalidad Provincial de Alto
Amazonas.  Laptop Core i7, marca Dell, modelo
 Municipalidad Distrital de Balsapuerto. VOSTRO
 Gobierno Regional de San Martín.  PC DELL modelo OPTIPLEX Core i5, con
 Municipalidad Distrital de Moyobamba. sistema operativo Windows 7
 Instituto Vial Provincial de Moyobamba.  Monitor SAMSUNG modelo SyncMAster
 Red vial nacional B2230 LED
 La presente información proporcionada  Monitor SAMSUNG modelo SyncMaster
por el Ministerio de Energía y Minas 940B
MINEM desde su base de datos para  Plotter Cannon modelo iPF810
descarga ubicada en su página web  Disco duro VERBATIM capacidad de
institucional (Figura 5), esta además almacenamiento 2 Tb
cuenta con información de Capitales de  GPS Garmin V y MAP 60
Distrito, Radiación Solar por Localidad,
Base de Datos de las Localidades a Nivel
Nacional, entre otros. 3.1.6. Software
 Para el ingreso de datos, generación de
3.1.4. Materiales interactivos de capas vectoriales y raster, gestión,
búsqueda y descarga de datos almacenamiento, análisis espacial y
en Red elaboración de mapas se utilizó el
programa en SIG ArcGIS 10.1 Desktop,
Se describen como las fuentes o servidores que
licencia del IIAP.
se encuentran en la Web tanto nacionales
 Para el ingreso de datos y análisis espacial
como internacionales las cuales proveen en su
se utilizó el programa en sistemas de
mayoría de datos e información satelital, esta
información geográfica Quantum GIS
se encuentra libre para la descarga, ordenada
(QGIS) versión 2.18. El QGIS Desktop
por ubicación geográfica y clasificada por día,
versión 2.18 Las Palmas, es un potente
mes y año de toma.
software SIG de escritorio para crear,
editar, visualizar, analizar y publicar
A continuación se enumera las fuentes de la
información geoespacial. Asimismo,
Web de las cuales se obtuvieron datos e
permite instalar complementos o
información satelital:
algoritmos especializados para realizar un
conjunto de tareas como por ejemplo el
Geoservidor del Ministerio del Ambiente, en el
complemento Semi-Automatic
cual se encuentra el catálogo de imágenes
Classification Plugin (SCP) el cual
Aster Global Digital Elevation Model (ASTER G-
proporciona adicionalmente
DEM), las cuales representan al mapa
herramientas como descarga de
topográfico mundial de resolución espacial de
imágenes satelitales disponibles sin costo,
el pre-procesamiento de imágenes y

10
 Evaluación del cambio de cobertura y uso de la tierra, en el área de influencia de la
propuesta de carretera Balsapuerto - Moyobamba

calculadora ráster, entre otras (Cogendo,

3.2. METODOLOGÍA
2016).
 Para la colecta y manejo de datos de
3.2.1. ÁREA DE ESTUDIO
campo mediante GPS Garmin V y Map 60
se utilizó el programa Map Source 6.0. El área de estudio comprende el área de
 En la validación de las capas temáticas de influencia entre las comunidades Balsapuerto
cobertura vegetal y uso actual se utilizó el (Loreto) y Moyobamba (San Martín). Se ha
programa Google Earth Pro. considerado esta área de influencia dado que
 En los cambios de la estructura de la base existe una propuesta de construir una carretera
de datos SIG se utilizó el programa que uniría estas dos localidades. Corresponde
Microsoft Access. parte de los distritos de Balsapuerto (Loreto) y
 En la generación de reportes, tablas, Moyobamba (San Martin). La extensión
cuadros y graficas estadísticos y redacción correspondiente a Balsapuerto es de 2 839,69
de informes fueron útiles los programas km2, y contiene a los centros poblados San
Microsoft Excel y Word. Jorge, Santa Rosa, San Lorenzo, Nueva Luz,
 En la recopilación de información en Puerto Libre, Canoa Puerto, Buenos Aires y
campo se utilizó el aplicativo en tablet Nuevo Junín. El área correspondiente a
GeoODK Collet, la cual proporciona una Moyobamba tiene una extensión de 2,737.57
forma de recopilar y almacenar Km2 y contiene a los centros poblados
información georreferenciada, junto con Cordillera Andina, Sugllaquiro, Creación 2000,
un conjunto de herramientas para Cóndor, Nuevo Progreso, Nuevo Piura, Flor de
visualizar, analizar y manipular datos Mayo, Juningue, Bella Selva.
terrestres para necesidades específicas.
Permite una comprensión de los datos El área de estudio limita por el norte con los
para la toma de decisiones, investigación, distritos de Barranca y Cahuapanas, en el
negocios, gestión de desastres, departamento de Loreto; por el este con los
agricultura y más. distritos de Caynarachi, departamento de San
Martín, Yurimaguas y Jeberos, departamento
de Loreto; por el sur con los distritos de San
Roque de Cumbaza, Pinto Recodo y Jepelacio,
departamento de San Martín; y por el oeste con
los distritos de Habana, Calzada y Yantalo,
departamento de San Martín. (Figura 3)

11
 Instituto de Investigaciones de la Amazonía Peruana

Figura 3. Mapa de Ubicación del área de influencia de la propuesta de carretera Balsapuerto -

Moyobamba.

12
 Evaluación del cambio de cobertura y uso de la tierra, en el área de influencia de la
propuesta de carretera Balsapuerto - Moyobamba

3.2.2. ACONDICIONAMIENTO b. Gestión de datos geoespaciales

CARTOGRÁFICO
En el proceso de actualización de los datos
Debido a que los datos fueron procesados cartográficos del proyecto se generan errores e
mediante el uso de Sistemas de Información inconsistencias tanto de tipo geométrico como
Geográfica (SIG) se propuso un proceso de tipo tabular, los cuales deben pasar por un
metodológico de acuerdo a la funcionalidad de proceso de acondicionamiento y
los SIG para lo cual este se dividió en las etapas almacenamiento de los datos con la finalidad
de ingreso y gestión de la información de que cumplan con los estándares
geoespacial, las cuales se detallan a cartográficos y puedan ser utilizados en el
continuación. análisis espacial para la generación de
información.
a. Ingreso de datos geoespaciales
En esta etapa se considera la recopilación de Generación de estructuras de
dos tipos de datos geoespaciales, los datos de almacenamiento
tipo vectoriales y los de tipo raster.
En las etapas de ingreso, generación y
actualización de los datos geoespaciales se
Recopilación de datos de tipo vectorial generan muchos archivos y de distintos
formatos los cuales tienen que llevar un orden
Se recolecto las capas de información,
lógico y ser almacenados en carpetas y
específicamente base fundamental de la Carta
subcarpetas de forma sistemática con la
Nacional, en formato digital, elaboradas por el
finalidad de su utilización, sostenibilidad y
Instituto Geográfico Nacional (IGN) a una escala
actualización en el tiempo.
de 1:100 000 las cuales están conformada por
En esta etapa se genera la estructura de
capas que representan a ríos, lagos y lagunas
almacenamiento que soportara los datos
(polígonos y líneas), centros poblados (puntos),
obtenidos en la etapa anterior, de acuerdo a las
curvas de nivel y cotas (líneas y puntos), en
normas para el tratamiento de datos e
total se recopiló cuatro cuadrantes de la carta
información geoespacial se recomienda los
nacional actualizadas al 2016 y en formato shp.
siguientes estándares:
También se recolectó capas de información
 Escala de trabajo aplicada, de 1:50 000.
base de red la hidrográfica (polígonos y líneas),
red vial (líneas) y capas de información  Sistema de proyección cartográfica,
temática en formato shp de los diversos Universal Transversal Mercator (UTM)
proyectos de Zonificación Ecológica y  Elipsoide de Referencia, World Geodetyc
Económica que el IIAP y otras instituciones System (WGS) del año 1984.
públicas han desarrollado en el área de estudio.  Datum World Geodetyc System (WGS) del
año 1984.
 Archivos vectoriales en formato estándar
Recopilación de datos de tipo raster shp.
En esta etapa se recopiló datos e información  Archivos raster en formato Tagged Image
en formato raster, los cuales representan el File Format (TIFF) o Geo Tif.
mundo real por celdas ordenadas en filas y  Estructura de almacenamiento para el
columnas llamadas pixel; donde a cada uno de proyecto en general será en el tradicional
estos se le asigna un valor, posteriormente con en carpetas de Windows.
estos valores se pueden realizar clasificaciones  Estructura de almacenamiento para los
y generar capas de información temática. datos geográficos e información generada
será en base de datos geográfica

13
 Instituto de Investigaciones de la Amazonía Peruana

(Geodatabase de ESRI) y en archivos de comprimido jpg, que representan mapas

lenguaje de marca extensible (XML) o resultados a utilizarse en el
 Proyectos de documentos de mapas complemento de documentos
elaborados en formato (MXD) informativos.
 Documentos de mapas exportados en  SHP, carpeta en la cual se almacena
formato de documento portátil (PDF) y en archivos de capas vectoriales en formato
formato comprimido de imágenes (JPEG). estándar Shapefile.
 DBF, carpeta en la cual se almacena
Con respecto al diseño de estructura de archivos de tablas en formato dbf, xls,
almacenamiento, se generó una estructura de etc.
almacenamiento de todos los archivos
generados en el estudio basándose en la
Acondicionamiento de cartografía base
utilización de carpetas de Windows y los
nombres propuestos de acuerdo a las siglas de Esta etapa consiste en realizar la detección de
los formatos de los archivos que fueron errores e inconsistencias, así como la limpieza
almacenados en estas. de los mismos, se realiza sobre las capas que
conforman la base cartográfica, se utiliza la
técnica de generación y corrección topológica a
La estructura propuesta se conforma de las
cada uno de los tipos de objetos que
siguientes carpetas:
representan (puntos, líneas y polígonos).
 CUT, carpeta matriz o raíz en la cual se
almacenan todas las demás carpeta y Acondicionamiento de cartografía
subcarpetas que contienen toda la temática
información correspondiente al proyecto, En esta etapa se realiza el acondicionamiento
hace referencia a las siglas de Cobertura y cartográfico y alfanumérico de las capas de
Uso de la Tierra. información temática de cobertura y uso de la
 DOC, carpeta en la cual se almacenan tierra, utilizando las curvas de nivel y la red
todos los documentos de texto que se hidrográfica principalmente contenidas en el
generaron en el proyecto, soporta mapa base.
archivos en formato pdf, doc, etc.
 IMG, carpeta en la cual se almacenan El acondicionamiento cartográfico comprende
documentos de tipo imagen y está específicamente la manipulación de los datos a
conformada por las siguientes nivel de geometría y de tabla de atributos, lo
subcarpetas: cual implica la detección y corrección de
 DEM, imágenes que representan Modelos errores topológicos en la geometría y la
Digital de Elevaciones DEM, las cuales se detección y corrección de campos vacíos e
encuentran en formato tif. inconsistencia en la tabla de atributos de cada
 SAT, imágenes de satélite en bandas, capa de información temática.
compuestas, corregidas y mosaicos, en
formato tif. En una primera instancia se realiza un
 RAS, imágenes o capas raster que han acondicionamiento cartográfico de cada capa
sido producto del procesamiento o temática por separado luego en una segunda
clasificación realizado a las imágenes de fase se realiza la integración de todas las capas
satélite o capas vectoriales transformados del componente físico del proyecto es aquí en
a formato raster. el cual mediante el uso de tablas de frecuencias
 JPG, carpeta en la cual se almacena los se detectó las inconsistencia a nivel de
documentos de imágenes en formato atributos.

14
 Evaluación del cambio de cobertura y uso de la tierra, en el área de influencia de la
propuesta de carretera Balsapuerto - Moyobamba

Teniendo en cuenta la estética en el diseño, el CLC se adoptó en Colombia en el marco del

programa permite elegir entre sus múltiples programa de apoyo y la colaboración prestada
paletas de colores que cuentan con índices de por el Fondo Francés para el Medio Ambiente
saturación y contraste. Para el caso de leyendas (FFEM), la Embajada de Francia en Colombia, el
con más de 10 colores se procedió a etiquetar Instituto Geográfico Nacional de Francia (ING) y
con un código numérico para una mejor lectura ONF Andina Colombia, con el fin de promover
del mapa. una metodología unificada y estandarizada para
el contexto colombiano en relación a
coberturas de la tierra. Así mismo, CLC es una
3.2.3. CARACTERIZACIÓN DE LA
COBERTURA Y USO ACTUAL DE metodología para la construcción de mapas de
LA TIERRA cobertura y uso de la tierra. Emplea una
leyenda jerárquica, que vincula distintos niveles
de detalle espacial (escala espacial) con
a. Identificación de las categorías de distintos niveles de detalle temático (niveles de
cobertura y uso de la tierra la leyenda jerárquica), y que en el 2010 paso
En esta etapa se realizó el preprocesamiento de por un proceso de discusión de la leyenda
las imágenes de satélite el cual consiste en regional ha llevado a la definición de una serie
definir y diferenciar los niveles de tonalidades de clases en los niveles 1 y 2 de la leyenda que,
mediante la modificación del histograma lo cual si bien recogen parte de los planteamientos de
generara una nueva imagen con niveles CLC, incorporan adaptaciones que se
digitales en una misma frecuencia, para esto se consideran necesarias para los países de la
utilizó la herramienta Image Analysis, la cual región (Valdivia, 2018; Arnillas et al., 2012).
viene integrada en el software ArcGIS 10.5,
luego se seleccionó la imagen a realzar y se Con los pixeles correspondientes a las áreas de
activó la herramienta Stretch, se realizó un entrenamiento se realizó el análisis de
remuestreo de los valores de intensidad, consistencia mediante una interpretación
igualando la presencia de estos en la imagen. bidimensional de los niveles digitales o
reflectividades que conforman cada una de las
Seguidamente en la imagen se generaron y clases muestreadas, seguidamente se generó
ubicaron polígonos de áreas de entrenamiento; un archivo de firma o descripción estadística
estas reúnen o agrupan pixeles de la imagen que constan del número de muestras, los
que previamente se hayan comprobado que promedios y las matrices de covarianza de las
pertenecen en su totalidad a una clase clases que derivaron de las muestras
determinada, luego estas mismas se nombran identificadas.
utilizando las clases por niveles de clasificación La identificación de categorías de cobertura y
según la leyenda CORINE Land Cover. Para una uso de la tierra, se realizó con los datos
mejor distribución de las áreas de verificados anteriormente y aplicando el
entrenamiento se generó una grilla de 1 km de método de clasificación por el algoritmo de
lado, distribuida en toda el área de estudio. Máxima Probabilidad (MLC por sus siglas en
inglés), en el cual los valores de pixeles en cada
El proyecto CORINE Land Cover (CLC) forma muestra de clases propuestas deben tener una
parte del programa CORINE, Coordination of distribución normal.
Information of the Environment, de la Agencia
Europea del Medio Ambiente, con el objetivo b. Verificación de las categorías de
fundamental de obtener una base de datos cobertura y uso de la tierra
europea de ocupación del suelo a escala
1:100.000, útil para el análisis territorial y la Para lograr una mejor consistencia y evaluación
gestión de políticas europeas. La metodología de la precisión del mapa de cobertura y uso

15
 Instituto de Investigaciones de la Amazonía Peruana

actual de la tierra se realizó la validación en evaluó todas las plantas desde 10 cm de

campo tanto de la cobertura vegetal como del diámetro a la altura del pecho (DAP) a más y en
uso actual de la tierra por separado y en las parcelas de 5 x 5 m ubicadas en las esquinas
paralelo. Antes de la salida de campo se generó de la parcela mayor, se evaluaron las plantas
el material cartográfico preliminar a validar, más pequeñas.
para el caso de la cobertura y uso de la tierra en
este se muestra la ubicación de puntos de Posteriormente las muestras botánicas se
control generados a partir de un muestro herborizaron de acuerdo al siguiente
sistemático estratificado sobre cada una de las procedimiento: Se colectó los individuos de
clases generadas, para el caso de la cobertura Angiospermas, Gimnospermas y Pteridophyta
vegetal, los mapas de campo deben mostrar las (helechos). Luego se codificó con códigos
unidades de muestreo las cuales fueron de dos numéricos. Seguidamente se registró en el
tipos 1) de muestreo no detallados y 2) las formato de recopilación de información de
parcelas de evaluación de la cobertura vegetal colectas de muestras botánicas. Luego se
detalladas, todas seleccionadas en base a un prensó las muestras con la ayuda de papeles
criterio fisiográfico. periódicos. Y se conservó con alcohol en bolsas
plásticas. El secado de las muestras se realizó
Para la recolección de información de campo y en Iquitos en ambientes especializados
con la finalidad de caracterizar el uso actual de “secadores”. Las muestras se identificaron con
la tierra y confirmar la dinámica del cambio en la ayuda de bibliografía especializada.
el tiempo se generaron encuestas
estructuradas, las cuales se aplicaron a los
3.2.4. CUANTIFICACIÓN DEL
propietarios de las parcelas visitadas, para esto CAMBIO DE LA COBERTURA
se utilizó el aplicativo web Geo ODK Collect, el Y USO DE LA TIERRA
cual necesita un archivo en formato “.xls” en el
cual se programa la encuesta, este último es
cargado a un servicio web ONA a. Preparación de los datos de
(https://odk.ona.io/), en el cual se almacenan ingreso al modelo
las encuestas ya realizadas. Así mismo, a la El objetivo del presente procedimiento es
encuesta se le programo de tal manera que exportar o convertir las capas de información
permita capturar un registro fotográfico y vectoriales de cobertura y uso de la tierra de los
tomar las coordenadas del lugar visitado en el años 2015 y 2019, mediante el uso de la
Sistema Geodésico Mundial (WGS por sus siglas herramienta Polygon To Raster del módulo
en inglés). Conversion Tools/To Raster del programa
ArcGIS. Uno de los parámetros importantes
Para la caracterización de la cobertura vegetal para la aplicación de la herramienta es el de
en campo, en primer lugar se conformaron las cell_assignment el cual determina el método
brigadas por un especialista botánico, un para asignar un valor a cada celda de la nueva
asistente de botánica y cuatro obreros, capa raster, el método elegido fue
seguidamente se ubicaron las unidades de MAXIMUM_COMBINED_AREA, debido a que las
muestreo, conformadas por parcelas de 20 x 50 capas de referencia poseen más de una
m de evaluación de la flora y puntos de categoría que pueden determinar un valor por
muestro de la vegetación. En el muestreo se celda creada y en este caso el valor de la celda
incluyeron hierbas, arbustos, arbolitos, árboles, es determinado por la categoría que este
palmeras y bejucos. Adicionalmente se contenida con mayor área. El tamaño de celda
colectaron especímenes de especies de epífitos, o pixel definido para la conversión fue de 10
seleccionando un árbol al azar dentro de la metros.
parcela; así mismo, dentro de la parcela se

16
 Evaluación del cambio de cobertura y uso de la tierra, en el área de influencia de la
propuesta de carretera Balsapuerto - Moyobamba

b. Análisis del cambio de cobertura mediante otras variables espaciales explicativas

y uso de la tierra que influencien de alguna manera; tales como
Tabulación cruzada la pendiente, subpaisaje fisiográfico, altitud
relativa, distancias a las carreteras, distancia a
Con las capas de información raster centros poblados y distancia a ríos principales,
anteriormente definidas se procedió a generar a continuación se detalla el procedimiento que
la matriz de Pontius utilizando el la implicó la generación del escenario:
herramientas Tabulate Area del módulo Spatial
Analyst Tools/Zonal del programa ArcGIS, el
a. Ingreso y análisis estadísticos de
cual genera una tabla de doble entrada en la
variables explicativas
cual en la primera columna se encuentran los
datos de las categorías del año 2015 y en las En esta etapa se ingresaron los mapas inicial
siguientes columnas se distribuye los valores de (t1) y final (t2) de CUT al modelo, así como las
áreas en m2 de las categorías de la capa raster variables espaciales, se verificó la geometría y
correspondiente al 2019. Con esta tabla se coincidencia de todas las capas de información
podrá realizar el análisis del cambio y de entrada. Seguidamente se verificó la
determinar las ganancias y pérdidas por correlación existente entre las variables
categorías de cobertura y uso de la tierra. espaciales en donde se describió tanto la fuerza
como la dirección de la relación, previamente
se realizó una prueba de normalidad. Para el
Identificación de cambios de CUT
presente estudio la técnica de análisis
Con la finalidad de cuantificar y obtener la estadístico fue mediante la prueba de
localización geográfica del cambio, se realizó Spearman, para datos no paramétricos
algunos procedimientos con las capas de utilizando el programa Sigmaplot.
información generadas en las fases anteriores
de los años 2010 y 2016 las cuales fueron
denominadas t1 y t2 respectivamente. Para
determinar la posición concreta de cada tipo de
cambio o estabilidad se realizó la intersección
de capas raster utilizando la herramienta
Combine del Spatial Analyst Tools/Local del
programa ArcGIS, el cual generó una nueva
capa raster con todas las combinaciones
posibles entre categorías.

3.2.5. ESCENARIO DEL CAMBIO DE LA

COBERTURA Y USO DE LA
TIERRA
Se propuso un modelo de predicción de
escenario a futuro el cual fue definido al año
2023, debido al periodo base de información,
este escenario se generó aplicando un
algoritmo de Autómatas Celulares el cual viene
incorporado en el módulo MOLUSCE que es
ejecutado en el programa QGIS versión 2.18,
para el presente estudio el MOLUSCE fue
utilizado para detectar el cambio en el uso de la
tierra entre los años 2015 y 2019 y medirlo

17
 Instituto de Investigaciones de la Amazonía Peruana

b. Medición del cambio los conjuntos de muestras y validaciones de

entrenamiento, almacena la mejor red
A partir de las capas inicial y final de cobertura
neuronal en la memoria. El proceso de
y uso de la tierra se obtuvo, en primer lugar las
entrenamiento finaliza cuando la mejor
estadísticas del cambio, seguidamente se
precisión es alcanzada.
generó una Matriz de Transición o de relaciones
de cambio, la cual se calculó mediante una
operación de multiplicación de matrices y d. Simulación de cobertura y uso
expresa las relaciones de cambio existentes futuro
entre las diferentes coberturas y usos de la
En esta etapa se aplicó el modelo de Automatas
tierra, finalmente en esta etapa se generó un
Celulares el cual se basa en el algoritmo
mapa de cambio de la cobertura y uso de la
probabilístico Monte Carlo. El mapa de
tierra en el periodo de estudio.
potencial de transición muestra la probabilidad
o el potencial de cambiar de una clase de
c. Modelamiento de la transición cobertura o uso de la tierra a otra. Los valores
potencial del potencial de transición oscilan entre 0 (bajo
potencial de transición de cambio) 100 (alto
En el presente estudio se utilizó el modelo de
potencial de transición). Los mapas de potencial
Red Neuronal Artificial (ANN por sus siglas en
de transición se producirán a partir de los
inglés) del módulo MOLUSCE del programa
correspondientes cambios en la cobertura y uso
QGIS. Para lo cual se definió el modelo de
de la tierra.
muestreo y número de muestras. El algoritmo
de aprendizaje analiza la precisión alcanzada en

18
 Evaluación del cambio de cobertura y uso de la tierra, en el área de influencia de la
propuesta de carretera Balsapuerto - Moyobamba

IV. RESULTADOS

4.1. CARACTERIZACIÓN DE LA COBERTURA Y USO ACTUAL DE LA TIERRA

A. IDENTIFICACIÓN DE LAS CATEGORÍAS DE COBERTURA Y USO DE LA TIERRA
Con respecto a la clasificación se obtuvo como resultado las clases definidas de cobertura y uso de la
tierra, realizando posteriormente una reclasificación y adaptándola a los niveles de la leyenda de
CORINE Land Cover. (Tabla 1).

Tabla 1. Propuesta de Leyenda de Uso y Cobertura de la Tierra - Perú (MINAM, 2012) y adaptada para
el estudio en el área de influencia de la propuesta de carretera Balsapuerto - Moyobamba.

NIVEL I NIVEL II NIVEL III NIVEL IV NIVEL V

1 Áreas 11 Áreas 112 Tejido urbano
artificializadas urbanizadas discontinuo
22 Cultivos
212 Arroz
transitorios
23 Pastos
241 Mosaico de
2 Áreas agrícolas
cultivos
24 Áreas agrícolas
243 Mosaico de
heterogéneas
cultivos, pastos y
espacios naturales
3111 Bosque
311 Bosque denso
denso alto
3112 Bosque 31121 Bosque denso
31 Bosques
3 Bosques y áreas denso bajo bajo de tierra firme
mayormente 31122 Bosque denso
naturales bajo inundable
33 Áreas abiertas 333 Tierras
sin o con poca desnudas y
vegetación degradadas
5 Superficies de 51 Aguas
511 Ríos
agua continentales

19
 Instituto de Investigaciones de la Amazonía Peruana

B. CARACTERIZACIÓN DEL USO nivel II, siete al nivel III, dos al nivel IV y dos al
ACTUAL DE LA TIERRA nivel V.
Las categorías de cobertura y uso del territorio
Las categorías de Bosque denso alto y Mosaico
que fueron tomadas en cuenta en la
de cultivos, pastos y espacios naturales, poseen
clasificación, se basan en la leyenda CORINE
la mayor superficie. La primera categoría con
Land Cover propuesta por el MINAM, la cual
63459 ha, que corresponde al 59% y la segunda
pretende uniformizar los tipos de cobertura y
categoría con un total de 28311 ha equivalente
usos de la tierra en los procesos de Zonificación
al 26% del área total y la categoría que posee la
Ecológica y Económica para el Ordenamiento
menor área corresponde a las áreas de arroz
Territorial.
con 261 ha que equivalen al 0.24%, es
importante mencionar también que la categoría
Se lograron clasificar 22 categorías de CUT,
Áreas agrícolas heterogéneas contempla una
cuatro pertenecen al nivel I, siete pertenecen al
superficie de 28726.5 ha correspondiente al
26.8%. (Tabla 2).

Tabla 2. Categorías por niveles de Cobertura y Uso de la Tierra adaptada al estudio del área de
influencia de la propuesta de carretera Balsapuerto – Moyobamba, leyenda de Cobertura y Uso de la
Tierra - Perú (MINAM, 2012).
Superficie Superficie
NIVEL I NIVEL II NIVEL III NIVEL IV NIVEL V
(ha) (%)
112 Tejido
1 Áreas 11 Áreas
urbano 304.33 0.28
artificializadas urbanizadas
discontinuo

22 Cultivos
212 Arroz 261.37 0.24
transitorios
23 Pastos
3553.57 3.32
241 Mosaico de
2 Áreas agrícolas 415.35 0.39
cultivos
24 Áreas
agrícolas 243 Mosaico de
heterogéneas cultivos, pastos
28311.15 26.47
y espacios
naturales
3111
311 Bosque
Bosque 63459.37 59.34
denso
denso alto
3112
3 Bosques y
Bosque 3215.1 3.01
áreas
31 Bosques denso bajo
mayormente
31121
naturales
Bosque
denso bajo 4594.32 4.30
de tierra
firme

20
 Evaluación del cambio de cobertura y uso de la tierra, en el área de influencia de la
propuesta de carretera Balsapuerto - Moyobamba

Superficie Superficie
NIVEL I NIVEL II NIVEL III NIVEL IV NIVEL V
(ha) (%)
31122
Bosque
282.72 0.26
denso bajo
inundable
33 Áreas sin o 333 Tierras
con poca desnudas y 1767.92 1.65
vegetación degradadas
5 Superficies de 51 Aguas 511 Ríos
784.06 0.73
agua continentales
TOTAL 106949.27 100.00

A continuación se describe los tipos de categorías de cobertura y uso de la tierra para el área de
influencia de la propuesta de carretera Balsapuerto - Moyobamba. (Figura 4).

Figura 4. Mapa de Uso Actual de la Tierra en el área de influencia de la propuesta de carretera

Balsapuerto - Moyobamba.

21
 Instituto de Investigaciones de la Amazonía Peruana

C. DESCRIPCIÓN DE LAS 2. Áreas agrícolas

CATEGORÍAS DE COBERTURA Y
USO ACTUAL DE LA TIERRA 2.1. Cultivos transitorios

1. Áreas artificializadas 2.1.1. Cultivo de Arroz

Esta unidad presenta una extensión SIG
1.1. Áreas urbanizadas de 261.37 ha, equivalente al 0.24% en la
zona de Moyobamba, San Martin. Se
1.1.1. Tejido urbano discontinuo presentan fraccionadas, distribuidas en el
Comprende una superficie de 304.33 ha, distrito Moyobamba, y cercanas a los
que representa el 0,28% del ámbito del centros poblados Flor de Mayo, Santa
estudio. La superficie de esta categoría en Catalina y Nuevo Piura, las que están
la zona de Loreto es de 188.52 ha (61.95%) conectadas por trochas carrozables y a lo
y está conformado por la sumatoria de las largo de toda la cuenca del río Mayo.
áreas urbano rurales del centro poblado
Balsapuerto, los caseríos Nuevo Junín, Son unidades dispuestas en un patrón
Canoa Puerto y Puerto libre. intrincado de mosaicos geométricos que
hace difícil su separación en coberturas
En la zona de San Martin asciende a 115.81 individuales; los arreglos geométricos
ha (38.05%) ha y lo conforman los centros están relacionados con el tamaño
poblados Quilloallpa, Sugllaquiro, reducido de los predios, las condiciones
Cordillera Andina y Nuevo Progreso en la locales de los suelos, las prácticas de
zona de Moyobamba, en el departamento manejo utilizadas y las formas locales de
de San Martin. tenencia de la tierra. (Figura 5).

Figura 5. Cultivo de arroz en el sector Moyobamba, San Martín.

22
 Evaluación del cambio de cobertura y uso de la tierra, en el área de influencia de la
propuesta de carretera Balsapuerto - Moyobamba

2.2. Pastos La concentración de pastizales en la zona

de Balsapuerto está en los caseríos Nueva
Esta área corresponde a 3553.57 ha, lo
Luz, San Jorge, Balsapuerto, Buenos Aires,
que representa el 3.32 % del total del
Nuevo Jerusalén, Nuevo Cuzco y Santa
área de estudio. Comprende las áreas
Rosa, y tiene una extensión de 127.87 ha
cubiertas con hierba densa de
(3.60%).
composición florística, dedicadas al
pastoreo permanente.
En la zona de Moyobamba estas áreas se
presentan fraccionadas, siendo las
Las especies producidas son braquiaria
mayores en extensión las cercanas a los
(Brachiaria decumbens) y pasto elefante
centros poblados Quilloallpa, Sugllaquiro,
(Pennisetum purpureum). Esta es
Cordillera Andina y Nuevo Progreso con
considerada la segunda actividad
una extensión de 3425.70 ha (96.40%).
productiva del ámbito de estudio en
(Figura 6).
cuanto a extensión de área cultivada.

Figura 6. Pastos en el sector Balsapuerto, Loreto.

23
 Instituto de Investigaciones de la Amazonía Peruana

2.3. Áreas agrícolas heterogéneas

Esta unidad presenta una extensión SIG 2.3.1. Mosaico de cultivos
de 28726.50 ha, equivalente al 26.86 %
Esta unidad presenta una extensión SIG
del área de estudio. En la zona de Loreto
de 415.35 ha, lo que representa el 0.39%
la extensión es de 9779.39 ha (34.04%),
del área total del estudio. La zona más
mientras que en la zona de San Martin es
extensa de halla en el distrito de
de 18947.11 ha (65.96%). Se presentan
Balsapuerto con una extensión de 367.30
fraccionadas, distribuidas en toda la
ha (88.43%) estas áreas están
extensión de los distritos de Balsapuerto y
centralizadas cercanas a los centros
Moyobamba y cercanas a las ciudades,
poblados Balsapuerto, Buenos Aires, San
centros poblados, redes carrozables de
Jorge, Santa Rosa, Nuevo Jerusalén, San
ambos distritos y a lo largo de toda la
Lorenzo y Nuevo Cachiyacu.
cuenca de los ríos Cachiyacu, Paranapura
en Loreto y del rio Mayo en San Martin.
En la zona de Moyobamba 48.05 ha
(11.57%) cercanas a los centros poblados
Son unidades que reúnen dos o más
Santa Catalina, Bella Selva, Nuevo Sullana,
clases de coberturas agrícolas y naturales,
Juningue. Incluye las tierras ocupadas con
dispuestas en un patrón intrincado de
cultivos anuales transitorios como arroz
mosaicos geométricos que hace difícil su
(Oryza sativa), maíz (Zea mays), ají
separación en coberturas individuales; los
(Capsicum annuum), frejol (Phaseolus
arreglos geométricos están relacionados
vulgaris), y cocona (Solanum
con el tamaño reducido de los predios, las
sessiliflorum). Cultivos permanentes
condiciones locales de los suelos, las
como café (Coffea arabica), cacao
prácticas de manejo utilizadas y las
(Theobroma cacao), naranja (Citrus
formas locales de tenencia de la tierra.
sinensis), limón (Citrus limon), caña
(Saccharum officinarum), plátano (Musa
Entre las especies representativas más
paradisiaca), yuca (Manihot esculenta),
frecuente de este frente figuran una
piña (Ananas comosus), papaya (Carica
variedad de cultivos como yuca (Manihot
papaya) y coco (Cocus nucifera) en los
esculenta), plátano (Musa paradisiaca),
cuales el tamaño de las parcelas es muy
maíz (Zea mays), piña (Ananas comosus),
pequeño (inferior a 25 ha) y el patrón de
papaya (Carica papaya), limón (Citrus
distribución de los lotes es demasiado
limón), naranja (Citrus sinensis), coco
intrincado para representarlos
(Cocus nucifera), arroz (Oryza sativa), café
cartográficamente de manera individual.
(Coffea arabica) y cacao (Theobroma
(Figura 7).
cacao).

24
 Evaluación del cambio de cobertura y uso de la tierra, en el área de influencia de la
propuesta de carretera Balsapuerto - Moyobamba

Figura 7. Cultivos de piña, papaya, plátano y café en los sectores Balsapuerto y Moyobamba.

2.3.2. Mosaico de cultivos, pastos y principalmente por coberturas de cultivos y

espacios naturales pastos en combinación con espacios naturales y
dispuestas en un patrón intrincado de mosaicos
Esta unidad presenta una extensión SIG de
geométricos que hace difícil su separación en
28311.15 ha, equivalente al 26.47% del área de
coberturas individuales; los arreglos
estudio. En la zona de Loreto la extensión es de
geométricos están relacionados con el tamaño
9412.08 ha (33.25%), mientras que en la zona
reducido de los predios, las condiciones locales
de San Martin es de 18899.07 ha (66.75%). Se
de los suelos, las prácticas de manejo utilizadas
presentan en una sola extensión distribuidas
y las formas locales de tenencia de la tierra.
uniformemente en todo el distrito de
Entre las especies representativas más
Balsapuerto que recorre los centros poblados
frecuente de este frente figuran una variedad
Canoa Puerto Puerto Libre, San Lorenzo,
de cultivos como yuca (Manihot esculenta),
Buenos aires, Balsapuerto, San Jorge, Santa
plátano (Musa paradisiaca), maíz (Zea mays),
Rosa, Nuevo Cachiyacu, Nuevo Junín, Bellavista,
piña (Ananas comosus), papaya (Carica
Santa Clara, San José, Nuevo Reforma, Monte
papaya), naranja (Citrus sinensis), limón (Citrus
Alegre y Nuevo Jerusalén cercanas a las trochas
limon), coco (Cocus nucifera), arroz (Oryza
carrozables de ambos distritos y a lo largo de
sativa), café (Coffea arabica), cacao
toda la cuenca de los ríos Cachiyacu,
(Theobroma cacao), braquiaria (Brachiaria
Paranapura en Loreto y del rio Mayo en San
decumbens), pasto elefante (Pennisetum
Martin. Son unidades que comprenden las
purpureum).
superficies del territorio ocupadas

25
 Instituto de Investigaciones de la Amazonía Peruana

Bosques y áreas mayormente naturales Sapotaceae, Simaroubaceae, Meliaceae,

Polypodiaceae, Araceae, Cyclanthaceae,
2.4. Bosques Myrtaceae, Orchidaceae, Acanthaceae,
Cecropiaceae, Hippocrateaceae, Poaceae,
2.4.1. Bosque denso Burseraceae, Marantaceae, Chrysobalanaceae,
entre otras.
2.4.1.1. Bosque denso alto
Esta Unidad muestra una extensión SIG de En cuanto a las especies las más
63459.37 ha, equivalente al 59,34 % de la representativas de estos bosques son: Clavija
superficie estudiada. En el área de estudio sp., Otoba glycycarpa, Heisteria sp., Iriartea
ubicada en Loreto alcanzan una extensión deltoidea, Pourouma sp., Ocotea sp., Neea sp.,
aproximada de 50767.49 ha (80%) y en área Inga sp., Otoba parvifolia, Piper sp., Miconia
correspondiente a San Martin llegan a sp., Minquartia guianensis, Pouteria sp., Virola
aproximadamente 12691.88 ha (20%). Son sp., Simarouba sp., Iryanthera juruensis, Clusia
bosques que se desarrollan en terrazas medias, sp., Brosimum parinarioides, Tapirira sp.,
altas, colinas altas y bajas y montañas altas y Aparisthmium cordatum, Clarisia biflora,
bajas. La fisonomía de esta vegetación se Sapium marmierii, Tapirira guianensis,
caracteriza por ser bosques altamente diversos Himatanthus sucuuba, Wettinia drudei, Trichilia
de estructura alta dominada por árboles; lo laxipaniculata, Urera verrucosa, Allophylus sp.,
conforman bosques con tres estratos que se Mauria heterophylla cf., Palicourea sp.,
traslapan irregularmente. Este frente se Stylogyne cauliflora, Alchornea glandulosa,
distribuye en dos sectores, uno al sur este de la Cecropia strigosa, Condaminea corymbosa,
capital del distrito Balsapuerto que se extiende Cyclanthus bipartitus, Ladenbergia sp.,
al sur y hacia el suroeste hasta el límite Myrcianthes rhopaloides, Naucleopsis sp.,
departamental con la región San Martin y el Socratea exorrhiza, Cheiloclinium sp., Pariana
otro sector que se extiendes desde el límite sp., Graffenrieda sp., Psycotria sp., Hevea
departamental con Loreto hacia el oeste y guianensis, Brosimum utile, Iryanthera laevis,
disperso entre los centros poblados del sector Tachigali sp. Osteophloeum platyspermum,
Moyobamba. Predominan las siguientes Micropholis guyanensis, Aspidosperma
familias botánicas: Myristicaceae, schultesii, Inga sp., Tachigali, Socratea
Theophrastaceae, Lauraceae, Olacaceae, exorrhiza, Micrandra spruceana , Tapirira
Arecaceae, Euphorbiaceae, Fabaceae, guianensis, Theobroma subincanum, Matisia
Melastomataceae, Rubiaceae, Moraceae, bracteolosa, Micropholis sp., Caryocar glabrum,
Nyctaginaceae, Anacardiaceae, Piperaceae, Brosimum lactescens, Aspidosperma,
Urticaceae, Annonaceae, Apocynaceae, Cespedesia spathulata, Hevea guianensis.
Clusiaceae, Malvaceae, Sapindaceae, (Figura 8).

26
Evaluación del cambio de cobertura y uso de la tierra, en el área de influencia de la
propuesta de carretera Balsapuerto - Moyobamba

Figura 8. Bosque denso alto en el sector Moyobamba.

27
 Instituto de Investigaciones de la Amazonía Peruana

2.4.1.2. Bosque denso bajo 2.4.1.2.1. Bosque denso bajo de tierra

firme
Esta Unidad muestra una extensión SIG de
3215.1 ha, equivalente al 3,01 % de la Esta Unidad muestra una extensión SIG de
superficie estudiada. En el área de estudio 4594.32 ha, equivalente al 4.3 % de la
ubicada en Loreto alcanzan una extensión superficie estudiada. Son bosques que se
aproximada de 3202.55 ha (99.61) y en área desarrollan en terrazas altas y bajas, en áreas
correspondiente a San Martin llegan a que no se inundan o en zonas de inundación
aproximadamente 12 ha (0.39%). Son bosques temporal o permanente. En el área de estudio
que se desarrolla en colinas altas y bajas. La de la zona de Loreto, alcanzan una extensión
fisonomía de esta vegetación se caracteriza por aproximada de 4431.94 ha (99.46) y en área
ser bosques altamente diversos de estructura correspondiente a San Martin llegan a
baja a mediano que se desarrollan sobre las aproximadamente 162.38 ha (3.53%). Se ubica
colinas. Se distribuye puntualmente en el al norte de Balsapuerto capital y cerca de
sector sur-oeste del área de estudio, cercano al Nuevo Junín, Monte Alegre, Santa Mercedes,
rio Mayo y la capital del distrito Moyobamba; y Zarumilla y al noreste de Balsapuerto entre las
en el sector al nor-este, cercano a la capital del comunidades nuevo barranquita, San Vicente,
distrito Balasapuerto y hacia la cordillera. Nueva Providencia y Bellavista. Y en sector de
Moyobamba pequeñas áreas dispersas
Predominan las siguientes familias botánicas: cercanas al centro poblado Juningue, Nuevo
Myristicaceae, Arecaceae, Burseraceae, Progreso, Santa Catalina y Nuevo Milagro. La
Dioscoreaceae, Euphorbiaceae, Fabaceae, composición florística está representada por las
Melastomataceae, Sapotaceae, Urticaceae, familias: Fabaceae, Myristicaceae,
entre otras. Elaeocarpaceae, Lecythidaceae,
Chrysobalanaceae, Apocynaceae, Sapotaceae,
En cuanto a las especies las más Cecropiaceae, Vochysiaceae, Euphorbiaceae,
representativas de estos bosques son: Protium Myrtaceae, Arecaceae, Clusiaceae,
sp., Pourouma sp., Eschweilera coriacea, Bombacaceae, entre otras.
Oenocarpus batahua, Calathea sp., Socratea
exorrhiza, Lepidocaryum tenue, Eschweilera La composición florística está representada
tessmannii, Otoba glycycarpa, Tetrastylidium principalmente por: Miconia klugii, Cecropia
peruvianum, Hevea guianensis, Iryanthera ficifolia, Leonia cymosa, Eschweilera sp.,
elliptica, Protium apiculatum, Dioscorea sp., Machaerium floribundum, Pariana sp., Cecropia
Maieta sp., Iriartea deltoidea, Pouteria sp., sp., Erisma laurifolium, Inga laurina, Stylogyne
Virola pavonis, Virola calophylla, Eschweilera longifolia, Amanoa nanayensis, Aspidosperma
albiflora, Rinorea sp. X2, Iryanthera juruensis, excelsum, Brownea ucayalina, Eschweilera
Marmaroxylon ramiflorum, Marmaroxylon albiflora, Iryanthera tessmannii, Oxandra
basijugum, Eschweilera parvifolia, Virola mediocris, Tabernaemontana sananho, Aniba
decorticans, entre otras. guianensis, Inga stenoptera, entre otras. (Figura
9).

28
 Figura 9. Bosque denso bajo de tierra firme en el sector Moyobamba.

2.4.1.2.2. Bosque denso bajo inundable

Esta Unidad muestra una extensión SIG de 282.72 ha, equivalente al 0.26% de la superficie estudiada
que se encuentra ubicada en Loreto. La característica principal de este bosque es que en su
composición se halla poca diversidad de especies de palmeras las que están dominadas por aguaje
(Mauritia Flexuosa) que se desarrolla sobre suelos pantanosos. Se encuentra distribuida en el sector
norte del área de estudio, cercano a la capital del distrito Balsapuerto.

La composición florística está representada principalmente por: Mauritia flexuosa, Socratea exorrhiza,
Nephrolepis rivularis, Calathea sp., Virola pavonis, Garcinia madruno, Zygia latifolia, Connarus
wurdackii, Euterpe precatoria, Tapirira guianensis, Ampelocera edentula, Virola sp., Cecropia sp.,
Conceveiba rhytidocarpa, Hymenolobium nitidum, Brosimum utile, Cecropia sciadophylla, Ceiba
samauma, Cordia nodosa, Cyclanthus bipartitus, Iryanthera tessmannii, Otoba glycycarpa, Piper
fortalezanum, Pleurothyrium brochidodromum, Psychotria sp., entre otras. (Figura 10).
 Instituto de Investigaciones de la Amazonía Peruana

Figura 10. Vegetación de Bosque denso bajo inundable en el sector Moyobamba.

30
 Evaluación del cambio de cobertura y uso de la tierra, en el área de influencia de la
propuesta de carretera Balsapuerto - Moyobamba

3. Áreas sin o con poca vegetación De forma general los principales cambios de
CUT en el periodo 2015 al 2019 fueron de la
3.1. Tierras desnudas degradadas categoría “Bosques y áreas mayormente
Esta Unidad muestra una extensión SIG de naturales” a “Áreas agrícolas” en 13365.57 ha
1767.92 ha, equivalente al 1.65% de la lo cual indica una pérdida de la primera
superficie estudiada. Se ubican en su categoría, así mismo, se puede observar que la
totalidad en la zona de Loreto en zonas de categoría “Áreas agrícolas” cambió a “Bosques
montañas altas y bajas, de laderas y áreas mayormente naturales” en 1247.14 ha
empinadas a extremadamente empinadas, lo cual indica una recuperación de la primera
como resultado de la ocurrencia de categoría. Otros cambios menores en superficie
procesos naturales de erosión y exposición se pueden observar de la categoría “Bosques y
a condiciones climáticas extremas. áreas mayormente naturales” a “Áreas
artificializadas” y “Superficies de agua” en
13.22 ha y 77.45 ha respectivamente, lo cual
4. Superficies de agua indica que en el periodo de análisis hubo una
pérdida de bosque a raíz de la ampliación de
4.1. Aguas continentales áreas urbanizadas y cambio en el curso de los
4.1.1. Ríos ríos por erosión lateral; otro cambio es el de
“Superficies de agua” a “Áreas agrícolas” y
En el presente estudio esta categoría “Bosques y áreas mayormente naturales” en
corresponde una superficie de 784.06 ha y 142.86 ha y 20.82 ha respectivamente.
representa 0.73 % de la superficie
estudiada. Se distribuye de forma dispersa Así también, relacionando todas las categorías
en las cuencas de los ríos Cachiyacu y de clase del 2015 al 2019 se pudo cuantificar y
Mayo. representar de forma específica el cambio de la
CUT, observándose el cambio de “Bosque” a
“Mosaico de cultivos, pastos y espacios
D. CUANTIFICACIÓN DEL CAMBIO DE
naturales” en 12404.38 ha, de “Bosque” a
LA COBERTURA Y USO DE LA
“Pasto” en 940.86 ha, de “Área agrícolas” a
TIERRA
“Tejido urbano discontinuo” en 175.68 ha y de
“Bosque” a “Ríos” en 77.45 ha. (Figura 11).

Figura 11. Variación de las superficies de Cobertura y Uso de la Tierra (CUT) por año, datos de CUT
2015 al 2019 y 2023 datos simulados.

31
 Instituto de Investigaciones de la Amazonía Peruana

E. ESCENARIO DEL CAMBIO DE LA la distancia a los cuerpos de agua se relaciona

COBERTURA Y USO DE LA TIERRA directamente con la altitud relativa con un 63%,
lo cual indica que las zonas de mayor altitud en
El modelo de generación de escenario al 2023
el área de estudio están alejadas de las riveras.
requirió dos capas de información en formato
Además cabe mencionar que la variable F no se
raster de CUT de los años 2015 y 2019 las
correlaciona con las demás variables.
cuales fueron generalizadas a cuatro clases o
categorías de CUT 1) Áreas artificializadas, 2)
Las superficies modeladas de CUT para el 2023
Áreas agrícolas, 3) Bosques y áreas
se distribuyen de la siguiente manera: 271.31
mayormente naturales y 4) Superficies de agua.
ha correspondiente a la categoría “Áreas
artificializadas”, 35158.47 ha para “Áreas
Así también, se obtuvieron siete variables
agrícolas”, 70685.14 ha para la categoría
espaciales explicativas obtenidas a partir del
“Bosques y áreas mayormente naturales” y
estudio fisiográfico y del Modelo Digital de
finalmente 834.73 ha para la categoría
Elevaciones, todas estas en formato raster, 1)
“Superficies de agua”.
Distancia a los centros poblados (DCP), 2)
Distancia a la red vial (DRV), 3) Distancia a los
Teniendo en cuenta los resultados de transición
cuerpos de agua polígonos (DCAP), 4) Distancia
potencial basado en las probabilidades del
a los cuerpos de agua líneas (DCAL), 5) Altitud
cambio y los motores que influyen en el mismo
relativa (A), 6) Pendiente (P) y 7) Elemento del
(variables explicativas), los resultados del
paisaje fisiográfico (F).
escenario indican que, si la tendencia actual
continúa, para el año 2023 según este
Con la finalidad de modelar y explicar la
escenario se predice un cambio en donde la
influencia de las variables anteriormente
superficie de la categoría “Bosques y áreas
generadas en algún tipo de cambio, se aplicó un
mayormente naturales” sufriría una pérdida de
análisis de correlación del cual se puede
2635 ha aproximadamente las cuales
observar que la distancia a los centros poblados
principalmente pasarían a Áreas agrícolas. Lo
se correlaciona en un 85% a la distancia a la red
cual se puede explicar principalmente por la
vial, lo cual indica una relación directa entre
influencia de las variables espaciales de tipo
estas variables, lo que comprueba la cercanía
socioeconómicas y centros poblados, así como,
que existe entre los centros poblados y su red
aspectos de relieve o fisiográficos propicios
vial; otra relación directa importante es aquella
para ciertos cultivos principalmente pastos,
entre la pendiente y la altitud relativa,
plátano y papaya en el sector Balsapuerto y
comprobando que al aumentar la altitud en una
café y cacao en el sector Moyobamba. (Figura
determinada distancia aumenta también la
12).
pendiente; Así también, se puede observar que

32
 Evaluación del cambio de cobertura y uso de la tierra, en el área de influencia de la
propuesta de carretera Balsapuerto - Moyobamba

Figura 12. Comparación de mapas de cobertura y uso de la tierra años 2015, 2019 y proyectado al año
2023.

33
 Instituto de Investigaciones de la Amazonía Peruana

V. BIBLIOGRAFIA

Ariza, A. (2013) Descripción y Corrección de Productos Landsat 8 LDCM (Landsat Data Continuity
Mission) versión 1.0. Cogota, Colombia.

Arnillas, C.A; Barrena, V; Llactayo, W; Ortíz, H; Regal, F; Rubín de Celis, E; Vásquez, P; Drenkhan, F;
Llerena, C.E; Espino, P; García, V; Gonzáles-Zúñiga, S; Gushiken, E; López, D. y Vásquez, R. 2012.
Informe Final del Proyecto: Análisis de las Dinámicas de Cambio de Cobertura de la Tierra en la
Comunidad Andina. Componente Nacional Perú – Primera Etapa. Centro de Datos para la
Conservación (CDC-UNALM). Laboratorio de Teledetección (LTAUNALM). Dirección General de
Ordenamiento Territorial - Ministerio del Ambiente (DGOT - MINAM). Fundación para el
Desarrollo Agrario (FDA). Lima - Perú. pp. 56.

Briceño Valera, Francisco Las cadenas de Markov en el análisis de cambios y asignación de usos de la
tierra Revista Geográfica Venezolana, vol. 46, núm. 1, enero-junio, 2005, pp. 35-45 Universidad
de los Andes Mérida, Venezuela.

Dirección General de Ordenamiento Territorial (2013) Guía Técnica para el Uso de las Imágenes
RapidEye. Lima, Perú.

Labrador Garcia, M.; Evora Brondo J. A.; Arbelo Perez, M. (2012) Satelites de Teledetcción para la
Gestion del Territorio. Canarias, España.

Paegelow, M, M. T. Camacho Olmedo y J. Menor T. 2003. Cadenas de Markov, evaluación

multicriterio y evaluación multiobjetivo para la modelización prospectiva del paisaje.
GeoFocus. Vol. (3): 22-44.

Lucchese L and Mitra S (133-136), "A New Method for Color Image Equalization", In Proc. of 2001 Int'l
Conference on Image Processing (ICIP 2001), Thessaloniki, Greece, Sept. 2001, Vol. I.

MacEachren A (1994), "Visualizacion in modern cartography", pp. 1-12. Pergamon Press, Inc.

Maguire D and Batty Goodchild M (2005), "GIS, spatial analysis and modeling" ESRI Press.

Mancebo Quintana, S.; Ortega Pérez, E.; Valentín Criado, A. C.; Martin Ramos, B.; Martin Fernández, L.
(2008) Libro SIG: aprendiendo a manejar los SIG en la gestión ambiental. Madrid, España, los
autores.

Pontius Jr., R. G.; Shusas, E.; McEachern, M., Detecting important categorical land changes while
accounting for persistence. Agriculture, Ecosystems and Environment 101 (2004) 251–268.

Ramírez-Sánchez H.U., García-Guadalupe M.E. y Villalpando-Piña G. EVALUACIÓN ESPACIO-TEMPORAL

DE VEGETACIÓN Y USO DE SUELO DE LA CUENCA RÍO CALIENTE, JALISCO, MÉXICO. 2013, pp. 43 -
55.

34
 Evaluación del cambio de cobertura y uso de la tierra, en el área de influencia de la
propuesta de carretera Balsapuerto - Moyobamba

Valdivia Márquez, L. N. (2018). Clasificación de la cobertura y uso del territorio para el año 2014 de la
provincia de Padre Abad, departamento de Ucayali.

Velázquez, A.; Mas, J. F.; Díaz Gallegos, J. R.; Mayorga Saucedo, R.; Alcántara, P. C.; Castro, R.;
Fernández, T.; Bocco, G.; Ezcurra, E.; Palacio, J. L. Patrones y tasas de cambio de uso del suelo en
México Gaceta Ecológica, núm. 62, 2002, pp. 21-37.

35