INTRODUCCIÓN
En los últimos siglos, como consecuencia del desarrollo económico y el crecimiento poblacional, se han presentado cambios significativos en el uso y cobertura de la tierra. El ejemplo más drástico de este proceso es la urbanización, la cual consiste en transformar y reemplazar los ecosistemas naturales por elementos artificiales como edificaciones y vías (Foley et al., 2005; Ellis y Ramankutty, 2008). Estas transformaciones, sumadas a las altas densidades poblacionales, generan unas condiciones ambientales particulares que afectan la calidad de vida de los habitantes urbanos (Patz, Campbell-Lendrum, Holloway y Foley, 2005; Basara, Basara, Illston y Crawford, 2010). Por esto, son fundamentales los servicios ecosistémicos que prestan las coberturas vegetales, al regular y mantener unas condiciones de vida sostenibles (Bolund y Hunhammar, 1999; Gómez-Baggethun et al., 2013; Lovell y Taylor, 2013).
Una de las características ambientales más contrastantes de las ciudades respecto a las áreas rurales aledañas es el aumento de la temperatura (Zhao et al., 2006; Peng et al., 2012; Livesley, McPhearson y Calfapietra, 2016). Este fenómeno conocido como isla de calor es producto principalmente de las características de los materiales empleados para la construcción de las ciudades, los cuales no permiten la evapotranspiración y absorben y retienen mayores cantidades de calor que la vegetación y otras coberturas naturales (Grimmond, 2007; Rizwan, Dennis y Liu, 2008). A su vez, la configuración espacial de las ciudades y la contaminación generada por las actividades diarias, incrementan el consumo de energía invertida para el enfriamiento de edificaciones, la demanda de combustibles fósiles y la emisión de gases de efecto invernadero (Santamouris, Cartalis, Synnefa y Kolokotsa, 2015). Por todo esto, la regulación de la temperatura es uno de los servicios ecosistémicos más relevantes que proveen las coberturas vegetales urbanas (McPhearson, 2011). A través de la evapotranspiración, el agua de las plantas absorbe calor y reduce la temperatura del aire, a la vez que se eleva la humedad relativa (Jim y Chen, 2009; Bowler, Buyung-Ali, Knight y Pullin, 2010). La provisión de sombra también regula la temperatura del aire y la temperatura superficial al reflejar la radiación solar y evitar que otras superficies absorban calor (Bolund y Hunhammar, 1999; Shashua-Bar y Hoffman, 2000; Hardin y Jensen, 2007; Bowler et al., 2010; Gómez-Baggethun et al., 2013).
Aunque la importancia del acceso a la vegetación urbana y sus servicios ecosistémicos es ampliamente reconocida, hay pruebas que demuestran que su distribución no es equitativa en muchas ciudades del mundo (Kabisch y Haase, 2014; Schwarz et al., 2015; Nesbitt y Meitner, 2016). Factores socioeconómicos como los ingresos monetarios, características étnicas y raciales, género, entre otros, han sido identificados como limitantes de la disponibilidad de coberturas vegetales urbanas tanto en cantidad como en calidad. Esto tiene implicaciones en la oferta de servicios ecosistémicos, la cual se ve disminuida en los grupos socioeconómicos menos favorecidos (Mitchell y Popham, 2008; Ernstson, 2013; Wolch, Byme y Newell, 2014). En Latinoamérica, en donde el 80 % de la población vive en ciudades y se presenta una marcada inequidad social (United Nations, 2018), se ha identificado que los residentes con mejores ingresos económicos tienen mayor acceso a espacios verdes y diversidad de especies arbóreas en ciudades como Santiago de Chile (De la Maza, Hernández, Brown, Rodríguez y Escobedo, 2002; Escobedo, Palmas-Perez, Dobbs, Gezan y Hernández, 2016) y Río de Janeiro (Pedlowski, Da Silva, Adell y Heynen, 2002). En el caso de Bogotá, se identificó menor abundancia de arbolado (Brown, 2012) y una reducción de sus parámetros estructurales en las zonas más pobres de la ciudad (Escobedo, Clerici, Staudhammer y Corzo, 2015). Esta problemática constituye una de las principales preocupaciones para la gestión urbana, en busca del bienestar generalizado y calidad de vida (Ernstson, 2013; Wolch et al., 2014).
Si bien la atención prestada a los ecosistemas urbanos es baja si se compara con la recibida por bosques y otros ecosistemas naturales (Gómez-Baggethun y Barton, 2013), el aumento en el número de publicaciones relacionadas con el mapeo de los servicios ecosistémicos indica que es un campo de investigación emergente (Ochoa y Urbina-Cardona, 2017). Su contribución como herramienta de gestión y planeación del territorio, resalta su importancia en las áreas urbanas (Egoh et al., 2008; Maes et al., 2012). Diferentes metodologías y aproximaciones se han empleado para el estudio de la distribución y el mapeo de los servicios ecosistémicos (Egoh, Drakou, Dunbar, Maes y Willemen, 2012; Maes et al., 2012). Una de las más empleadas, y que busca lidiar con la falta de información primaria, es el uso de estimaciones basadas en el tipo de cobertura (Burkhard, Kroll, Müller y Windhorst, 2009; Burkhard, Kroll, Nedkov y Müller, 2012; Schneiders, Van Daele, Van Laduyt y Van Reeth, 2012). No obstante, estas estimaciones han presentado ajustes bajos con información primaria, por lo que son recomendadas para determinar patrones de distribución a escala gruesa (Eigenbrod et al., 2010). A escalas locales como las ciudades, los datos estructurales medidos en cada árbol a través de censos forestales, son de gran valor para el estudio de los servicios ecosistémicos, no obstante, pocas ciudades cuentan con inventarios forestales de este tipo debido a sus altos costos (Nowak et al., 2008).
Pese a su relevancia, el mapeo de la regulación térmica en las ciudades ha sido poco caracterizada debido a que generalmente se relaciona con el almacenamiento y captura de carbono. Sin embargo, los ecosistemas ofrecen este servicio a diferentes niveles: en las ciudades a través de la provisión de sombra y la evapotranspiración vegetal, y a escalas espaciales mayores a través del almacenamiento y captura de carbono, contribuyendo a la regulación climática global (Maes et al., 2017). Esto implica la intervención de diferentes funciones ecosistémicas y el uso de distintos indicadores para su estudio (Egoh et al., 2012). Como antecedente del mapeo de la capacidad de regulación de temperatura urbana, se encuentran los trabajos presentados por Zulian et al. (2017) y Cortinovis y Geneletti (2018), en los que a través de una metodología basada en literatura se identificó la provisión de sombra y la evapotranspiración como las funciones principales asociadas al servicio ecosistémico; y, posteriormente, mediante el uso de indicadores se generó un mapa representando la capacidad de regulación térmica en la ciudad de Trento, Italia. De manera similar, Zardo, Geneletti, Pérez-Soba y Van Eupen (2017) desarrollaron una metodología para evaluar el potencial de regulación térmica de la infraestructura verde, identificando nuevamente la provisión de sombra y la evapotranspiración como funciones principales, y luego la aplicaron en la ciudad de Ámsterdam.
El objetivo de esta investigación fue analizar la distribución espacial de la infraestructura verde pública de Bogotá y su capacidad de regulación térmica, en relación con la estratificación socioeconómica, esto con el fin de proveer información que contribuya al mejoramiento de la provisión de este y otros servicios ecosistémicos. De acuerdo con los antecedentes, se espera identificar diferencias en la distribución de la cobertura arbórea, el verdor de los parques públicos y por consiguiente del potencial de regulación de temperatura en las áreas de la ciudad con menor estrato socioeconómico.
MATERIALES Y MÉTODOS
Área de estudio
Bogotá es la capital y ciudad más grande de Colombia. Se localiza en el centro del país sobre la cordillera oriental, a una altitud promedio de 2630m (figura 1). Cuenta con una población de alrededor de 8 millones de habitantes (Cristancho y Triana, 2018), y una superficie urbana de aproximadamente 380 km2. Durante los últimos años, factores como la contaminación y la urbanización han hecho más cálida la ciudad y acentuado la presencia del efecto isla de calor en las zonas más internas, las cuales superan en aproximadamente 3 °C la temperatura de las áreas periféricas (Ángel, Ramírez y Domínguez, 2010).
La división política y administrativa de la ciudad comprende 19 localidades que, a su vez, se subdividen en 112 UPZ (Unidades de Planeamiento Zonal), que corresponden a áreas con usos del suelo y desarrollos urbanísticos similares (Secretaría Distrital de Planeación, 2009). La división socioeconómica se representa a través de los estratos socioeconómicos, los cuales fueron establecidos por el gobierno distrital como una aproximación para capturar las diferencias socioeconómicas entre zonas de la ciudad y como una herramienta para direccionar efectivamente el gasto público hacía los ciudadanos menos favorecidos (Hernández, 2016). Existen seis estratos socioeconómicos, de los cuales el uno corresponde a las zonas más pobres de la ciudad, con carencias en infraestructura y servicios públicos, mientras que el seis representa las zonas de la ciudad donde se concentran los mayores ingresos monetarios y una gran oferta de infraestructura y servicios (Brown, 2012; Escobedo et al., 2015).
Con respecto a la infraestructura verde, la ciudad cuenta con la denominada Estructura Ecológica Principal (EEP), que agrupa los elementos fundamentales para la conservación de la biodiversidad y los procesos ecológicos esenciales para el desarrollo sostenible (Isaacs y Jaimes, 2014). Dentro de los componentes de la EEP están los parques urbanos, espacios verdes de uso público que sirven como reguladores del equilibrio ambiental y brindan espacios de recreación y aprovechamiento para los ciudadanos. Los parques urbanos se clasifican en cuatro categorías de acuerdo con su tamaño y área de influencia: parques metropolitanos (> 10 ha), parques zonales (1 a 10 ha), parques vecinales (1000 m2 a 1 ha) y parques de bolsillo (< 1000 m2) (Secretaría Distrital de Planeación, 2009). El arbolado público de la ciudad se distribuye tanto en parques como en los diferentes tipos de espacio público tales como separadores, andenes, entre otros. Algunas de las especies más abundantes son el sauco (Sambucus nigra L.), el jazmín del Cabo (Pittosporum undulatum Vent.) y la acacia japonesa (Acacia melanoxylum R. Br.) (Mahecha et al., 2010).
Estimación de la cobertura arbórea
Para estimar la cobertura arbórea se obtuvo información del Sistema de Información para la Gestión del Arbolado Urbano de Bogotá (Sigau) (tabla 1), el cual consiste en un inventario del arbolado público de la ciudad. La base de datos incluye datos morfométricos de aproximadamente 1300 000 individuos con altura igual o superior a 30 cm, y su localización en coordenadas geográficas (Dirpen, 2006). La cobertura arbórea se estimó con respecto a todos los individuos, modelando la proyección de la copa al suelo, asumiendo que el área de la copa corresponde a una circunferencia cuyo diámetro es equivalente al diámetro ecuatorial del individuo (Rodríguez-Laguna, Meza-Rangel, Vargas-Hernández y Jiménez-Pérez, 2009). Esta variable hace parte del Sigau y fue medida en terreno como la mayor longitud de la proyección ortogonal de la copa (Dirpen, 2006).
Cálculo del índice de vegetación de diferencia normalizada (NDVI)
Los índices de vegetación (IV) se fundamentan en las características físico-químicas de la vegetación que determinan su respuesta a la interacción con las diferentes longitudes de onda de la radiación electromagnética. La vegetación es altamente reflectiva en la región del espectro electromagnético correspondiente al infrarrojo cercano y a su vez presenta una baja reflectividad en la región del rojo. Por tal motivo, muchos de los IV usan las bandas del rojo y el infrarrojo cercano de las imágenes satelitales con el fin de acentuar las propiedades espectrales de la vegetación y atenuar las características de las demás coberturas (Chuvieco, 2016). El NDVI, uno de los IV más comunes, tiene numerosas aplicaciones en estudios ecológicos (Kerr y Otrovsky, 2003; Glenn, Huete, Nagler y Nelson, 2008), incluyendo el análisis y cuantificación de servicios ecosistémicos (Ayanu, Conrad, Nauss, Wengmann y Koellner, 2012). Múltiples autores han empleado este índice como una medida de verdor en zonas urbanas como resultado de la cantidad de vegetación o coberturas permeables presentes en un área (Zhou, Troy, Grove y Jenkins, 2009; Leslie, Sugiyama, Ierodiaconou y Kremer, 2010). Adicionalmente, diferentes estudios han encontrado una relación inversa entre la temperatura de la superficie y el NDVI (Weng et al., 2004; Weng, 2009; Huang y Ye, 2015).
Para calcular el NDVI, se utilizó una imagen multiespectral obtenida por el sensor OLI a bordo de la plataforma Landsat 8 (tabla 1). La cobertura de nubes alcanzó el 25 %, pero no hubo incidencia sobre el área de interés. La imagen se adquirió preprocesada, tanto geométrica (Storey, Choate y Lee, 2014) como radiométricamente (Chander, Markham y Helder, 2009) y transformada a unidades de reflectancia de la superficie por el Servicio Geológico de los Estados Unidos, mediante la aplicación del algoritmo LaSRC (Landsat Surface Reflectance Code) (United States Geological Survey, 2017). De esa forma, se obtuvo como resultado una estimación de la reflectancia espectral de la superficie que puede ser aplicada para estudios de vegetación y el cálculo más preciso de índices espectrales (Chuvieco, 2016). Sobre esta imagen se calculó el NDVI así (Rouse, Haas, Schell, Deering y Harlan, 1974):
Donde, IRC corresponde a los valores de reflectancia de la banda del infrarrojo cercano (OLI5) y ROJO a los valores de reflectancia de la banda del rojo (OLI4). El intervalo del NDVI va de -1 a 1, pero la vegetación generalmente registra valores por encima de 0.1, donde los valores más altos están asociados con mayor densidad de coberturas vegetales saludables (Jensen, 2000). Por tal motivo, se reclasificó para obtener valores entre 0 y 1, intervalo en el que se hace énfasis en el verdor.
Análisis espacial
Se usaron las capas de cobertura arbórea y NDVI previamente generadas, y otras capas del mapa de referencia para Bogotá en la versión correspondiente al mes de junio de 2017, el visor geográfico de la Secretaría Distrital de Ambiente y el portal de mapas de Bogotá (tabla 1). Los procesos y herramientas necesarios para los análisis se ejecutaron en el software QGIS 2.18 (QGIS Development Team, 2011). Para determinar la distribución de la cobertura arbórea en las zonas de la ciudad con diferente estrato socioeconómico, se realizó un análisis desde el ámbito de las UPZ. Para ello, se generaron capas de cobertura arbórea y estrato socioeconómico por UPZ que posteriormente se cruzaron para obtener el porcentaje del área total de cada estrato socioeconómico representada por cobertura arbórea (figura 2). En el caso del NDVI se obtuvo valores medios para cada UPZ y para cada parque público a partir de los valores registrados en cada píxel. Como consecuencia de la variación en la forma y tamaño de los parques, y la resolución espacial de la imagen satelital (900 m2), los polígonos no se distribuyeron de manera uniforme sobre los píxeles. Para lidiar con este problema, el NDVI se remuestreó con un tamaño de píxel de 1 m2 antes de la extracción de las estadísticas zonales, ponderando así cada píxel de tamaño original por el área intersectada por cada uno de los polígonos de los parques (Wilson et al., 2016).
Mapeo de la capacidad de regulación térmica
Inicialmente, se identificaron la provisión de sombra y la evapotranspiración como las funciones principales que influencian la capacidad de la infraestructura verde para la regulación térmica (Zardo et al., 2017; Zulian et al., 2017; Cortinovis y Geneletti, 2018). Se usó el porcentaje de cobertura arbórea como proxy de la cantidad de sombra ofrecida por el arbolado urbano (Shashua-Bar y Hoffman, 2000), la cual ha sido ampliamente reconocida como uno de los mecanismos de regulación de temperatura en las ciudades (Bolund y Hunhammar, 1999; Hardin y Jensen, 2007; Bowler et al., 2010; Gómez-Baggethun et al., 2013). La media del NDVI de los parques públicos de cada UPZ, transformado a valores entre 0 y 100, se usó como proxy del verdor o abundancia de coberturas vegetales o superficies con capacidad de evapotranspirar (Zhou et al., 2009; Leslie et al., 2010), teniendo en cuenta que el aumento en el NDVI se ha relacionado con la disminución de las temperaturas urbanas (Weng et al., 2004; Weng, 2009; Huang y Ye, 2015).
Posteriormente, se evaluó cada una por separado y finalmente se combinaron para generar un mapa coroplético que representa el potencial de regulación térmica de la infraestructura verde pública de Bogotá con respecto a las UPZ. La combinación de estás dos variables para determinar dicho potencial, corresponde a su importancia para la regulación climática urbana y lo señalado en las metodologías propuestas por Zardo et al. (2017), Zulian et al. (2017) y Cortinovis y Geneletti (2018). Respecto a la contribución relativa de las variables a la capacidad de regulación térmica, se ha identificado que varía en función del tamaño del área analizada cuando se toman parques y elementos específicos de la infraestructura verde como unidad de análisis, pero no hay suficiente información en la literatura sobre cómo combinar su contribución (Zardo et al., 2017). Por lo anterior, y teniendo en cuenta que no se consideró el área de los parques o las UPZ como un factor de análisis, se optó por seleccionar una contribución balanceada, asignando los mismos pesos a las dos variables. De esa manera, cada UPZ obtuvo un puntaje entre 0 y 200, como resultado de la adición de las variables. Usualmente en este tipo de combinaciones los valores se estandarizan entre 0 y 100, pero como los porcentajes de cobertura arbórea y los valores de NDVI fueron bajos, y ninguna combinación excedió el valor de 100, no se consideró necesaria la estandarización.
Análisis estadístico
Para evaluar el grado de autocorrelación espacial de las variables analizadas, se empleó el índice I de Moran como medida de autocorrelación global, e indicadores locales de asociación espacial (Lisa) como medida de autocorrelación local. La I de Moran, permite analizar el grado de dependencia de los valores que toma una variable en un lugar respecto a los valores de las ubicaciones vecinas. Por su parte, los análisis locales permiten identificar clústeres o agrupaciones de lugares en los que se concentran valores altos de una variable u hotspots, y valores bajos o coldspots (Celemin, 2009). Los análisis geoestadísticos fueron realizados con el software GeoDa® 1.10 (Anselin, Syabri y Kho, 2006), definiendo las vecindades de las UPZ por contigüidad y teniendo en cuenta todas las direcciones para su determinación (caso queen). Con el fin de identificar diferencias en el valor del NDVI y el porcentaje de cobertura arbórea de las UPZ según su estrato socioeconómico, y el NDVI de los parques de acuerdo con su clasificación y el estrato socioeconómico donde se localizaron, se aplicó la prueba no paramétrica de Kruskal-Wallis y la prueba post-hoc de Dunn, teniendo en cuenta que los datos no presentaron una distribución normal (Zar, 2010). Estos análisis se realizaron utilizando el programa estadístico R (R Core Team, 2016).
RESULTADOS
La distribución de la cobertura arbórea en Bogotá evidenció un comportamiento sesgado según el estrato socioeconómico (H(5) = 48.401, p < 0.001). Se identificaron 3 grupos de estratos en los que se observa que las zonas de la ciudad con estratos 1 y 2 presentaron la menor proporción de cobertura arbórea respecto a los demás (figura 3a). El verdor de las UPZ de los diferentes estratos, representado como el NDVI, presentó diferencias significativas (H(5) = 16.5, p < 0.005) pero no se evidenció un patrón ni agrupaciones claras como resultado del análisis post-hoc (figura 3b). Con respecto al estrato de la zona en donde se ubican, el verdor de los parques públicos también presentó diferencias (H(5) = 196.88, p < 0.001), aumentando en los estratos más altos, con excepción del 6 que no fue significativo (figura 3c). El análisis de los parques según su clasificación arrojó diferencias en el verdor (H(3) = 221.15, p < 0.001), presentando una tendencia de disminución en parques más pequeños y con menor área de influencia. Los parques metropolitanos mostraron los valores más altos, seguidos por los zonales y vecinales. Por su parte, los parques de bolsillo, evidenciaron el verdor más bajo (figura 3d).
![](/img/revistas/cofo/v22n2//0120-0739-cofo-22-02-00083-gf3.jpg)
Figura 3 a) Porcentaje de cobertura arbórea en las UPZ de cada estrato socioeconómico; b) NDVI de las UPZ según su estrato socioeconómico; c) NDVI de los parques públicos según el estrato socioeconómico donde se ubican; d) NDVI según tipo de parque urbano. Letras no repetidas representan valores diferentes (p < 0.05).
El estrato socioeconómico más abundante fue el 3, mientras que el más escaso fue el 6 que se presentó solo en 4 de las 112 UPZ de la ciudad (figura 4a). El porcentaje de cobertura arbórea varió entre 0 y 24 %, en donde la mayoría de las UPZ presentaron valores por debajo del 3 % y solamente 4 superaron el 10 % (figura 4b). El NDVI de los parques urbanos presentó una mayor heterogeneidad espacial; sin embargo, se pudo evidenciar al igual que en las anteriores variables una brecha marcada en la que la UPZ con parques urbanos más verdes tuvo un valor más de 6 veces mayor que aquellas con menor abundancia de coberturas vegetales en sus parques (figura 4c). Finalmente, pasando al potencial de regulación de temperatura, presentó un intervalo amplio, con valores que van desde menos de 9, hasta aproximadamente 84 (figura 4d).
![](/img/revistas/cofo/v22n2//0120-0739-cofo-22-02-00083-gf4.jpg)
Figura 4 a) Estratificación socioeconómica de la ciudad de Bogotá; b) distribución espacial de la cobertura arbórea de la ciudad; c) media ponderada del NDVI de los parques urbanos; d) puntaje de potencial de regulación térmica.
En general, los análisis del índice I de Moran registraron una tendencia débil a la autocorrelación global positiva. En el caso del estrato socioeconómico, se presentó la autocorrelación positiva más alta con un valor de I de Moran de 0.58, la cobertura arbórea registró 0.3, el potencial de regulación térmico 0.28 y el NDVI de los parques 0.24. Con un valor de p < 0.001 para todos los casos, la prueba de hipótesis estadística rechaza la hipótesis nula de aleatoriedad espacial o ausencia de autocorrelación. La figura 5 muestra los resultados de los análisis de LISA univariados para la determinación de autocorrelación espacial local. Los clústeres del estrato socioeconómico evidencian una clara división de la ciudad, en donde los estratos socioeconómicos altos, se agrupan en la zona norte y los bajos, en la zona sur de la ciudad (figura 5a). La zona suroccidental de la ciudad presentó un coldspot para el porcentaje de cobertura arbórea (figura 5b). Los parques públicos registraron valores bajos de NDVI agrupados principalmente al sur de la ciudad y los valores altos se agruparon predominantemente al norte (figura 5c). Finalmente, el potencial de regulación de temperatura registró nuevamente coldspots al sur de la ciudad y hotspots en la zona norte (figura 5d). En general, las variables analizadas presentaron intersección de algunas áreas de sus coldspots en el sur y suroccidente, y hotspots en el norte.
![](/img/revistas/cofo/v22n2//0120-0739-cofo-22-02-00083-gf5.jpg)
Figura 5 Análisis de LISA univariados (agrupaciones) para las variables analizadas en la ciudad de Bogotá: a) mapa de agrupaciones para el estrato socioeconómico; b) mapa de agrupaciones para el porcentaje de cobertura arbórea pública; c) mapa de agrupaciones para la media ponderada del NDVI de los parques públicos; d) mapa de agrupaciones para el puntaje de potencial de regulación térmica. Resultados significativos con p < 0.05.
DISCUSIÓN
El presente estudio analizó la distribución de la infraestructura verde pública de la ciudad de Bogotá con respecto a las UPZ, con énfasis en el arbolado y los parques urbanos, y su potencial para regular la temperatura con el fin de evaluar la inequidad ambiental. Aunque otros estudios han analizado la inequidad ambiental de Bogotá (Brown, 2012; Escobedo et al., 2015), ninguno ha analizado la distribución de la cobertura arbórea usando las UPZ como unidad de análisis, ni ha evaluado el verdor con datos obtenidos mediante teledetección. De igual manera, estos análisis han tenido como objetivo evaluar servicios ecosistémicos como la remoción de contaminantes del aire y el almacenamiento de carbono en biomasa, pero ninguno ha relacionado la distribución de las coberturas vegetales urbanas con el potencial de regulación térmica.
Los resultados obtenidos evidencian que la distribución de la infraestructura verde y su potencial para regular la temperatura es heterogénea en la ciudad de Bogotá respecto a los estratos socioeconómicos. En el caso del porcentaje de cobertura arbórea, se evidenció una clara tendencia a la disminución en las áreas con menor estrato socioeconómico. Esto concuerda con estudios realizados previamente en Bogotá, los cuales indicaron que tanto la cantidad de árboles (Brown, 2012) como su diversidad, área de copa y otros atributos estructurales (Escobedo et al., 2015), son mayores en los estratos más altos. De manera similar, este patrón de mayor abundancia y calidad de la infraestructura verde en las zonas con mejor condición socioeconómica, ha sido observado en otras ciudades de Latinoamérica (De la Maza et al., 2002; Pedlowski et al., 2002; Escobedo et al., 2016), Norteamérica (Schwarz et al., 2015; Nesbitt y Meitner, 2016) y Europa (Kabisch y Haase, 2014). Teniendo en cuenta que los estratos socioeconómicos buscan capturar las diferencias socioeconómicas de la población (Hernández, 2016), los resultados obtenidos indican problemas de inequidad ambiental en Bogotá, tal como han concluido previamente otros autores (Brown, 2012; Escobedo et al., 2015).
Las diferencias en la distribución de la cobertura arbórea se pueden explicar a partir del rápido aumento poblacional que ha experimentado Bogotá durante los últimos años (Cristancho y Triana, 2018) como respuesta a la tendencia global de migración hacia las áreas urbanas (Aide y Grau, 2004), particularmente acentuada en Colombia por decenios de conflicto armado (Shultz et al., 2014). Esto generó un proceso descontrolado de ocupación en las zonas periféricas de la ciudad sin ningún tipo de planificación (Romero, 2010; Muñoz y Ducón, 2016). De acuerdo con Romero (2010), a partir de la década de 1970, el crecimiento de la ciudad por la presión poblacional se direccionó hacía el norte, sur y suroccidente de la ciudad, los dos últimos casos como consecuencia de desplazamientos rurales a zonas urbanas, generando procesos de parcelación, urbanización, edificación desordenada y reducción de áreas libres hasta aproximadamente 5 %. Teniendo en cuenta que la historia de los sistemas urbanos, la interacción de sus componentes y la configuración espacial en un momento específico determinan las posibilidades futuras de manejo (Pickett et al., 2017), la ocupación desordenada de las áreas periféricas al sur y suroccidente de la ciudad se refleja actualmente en la baja cantidad de áreas verdes y cobertura arbórea en estas áreas (Brown, 2012; Escobedo et al., 2015).
Se ha reportado que la historia de urbanización en la ciudad de Bogotá no solo afecta la disponibilidad de coberturas vegetales sino también aspectos fundamentales de la flora urbana como la composición y la diversidad de especies en los jardines, y la forma en la que los habitantes de las diferentes áreas interactúan con estos elementos (Sierra-Guerrero y Amarillo-Suárez, 2017). Esta inequidad en el acceso al arbolado urbano implica también inequidad en el disfrute de los múltiples beneficios que esto genera para la salud pública (Ernstson, 2013; Wolch et al., 2014). Mitchell y Popham (2008), además, plantean que los vecindarios con menores ingresos económicos enfrentan una serie de problemáticas asociadas a la salud y mayores tasas de mortalidad, y como consecuencia estas poblaciones vulnerables recibirían un mayor beneficio neto de estas coberturas que aquellas que gozan de mejores capacidades económicas.
Las comparaciones del verdor permiten observar que los estratos 6, 5 y 1 registran los valores más altos respectivamente. Esto puede ser consecuencia de analizar no solo las áreas verdes públicas, sino toda la superficie de las UPZ, teniendo en cuenta que en el caso de los estratos más altos se presenta el mayor número de áreas verdes privadas y en el estrato 1 gran influencia de áreas verdes periféricas de la ciudad. Por otra parte, podría ser el resultado de la heterogeneidad espacial de las ciudades y la resolución espacial de los datos empleados para los análisis (Pickett et al., 2017; Zhou, Pickett y Cadenasso, 2017). Al analizar los parques según su estrato, se observan valores más altos de verdor a medida que incrementa el estrato, con un comportamiento similar a la cobertura arbórea. Nuevamente, esto tiene implicaciones en el disfrute de los beneficios generados por los parques en las diferentes zonas de la ciudad (Tzoulas et al., 2007; Bowler et al., 2010). El NDVI ha sido relacionado inversamente con la temperatura superficial en áreas urbanas (Weng et al., 2004; Weng, 2009; Huang y Ye, 2015) y con toda una serie de procesos fisiológicos tales como la capacidad de evapotranspiración, que determinan la funcionalidad de las coberturas vegetales no solo para regular la temperatura sino para la generación de bienestar en las poblaciones humanas (Kerr y Otrovsky, 2003; Glenn et al., 2008; Ayanu et al., 2012). Particularmente en las ciudades, se ha identificado la importancia del NDVI y otras medidas de verdor en los vecindarios, sobre la calidad de vida por su relación con la salud física, psicológica y la sensación de bienestar (Zhou et al., 2009; Leslie et al., 2010). Conjuntamente, lo anterior sugiere que la calidad de estas áreas verdes respecto a su capacidad de regulación de temperatura y a la provisión de otros servicios y beneficios es mayor en las zonas con mayor estrato socioeconómico.
De igual forma, el NDVI cambió de acuerdo con el tipo de parque. Si bien la clasificación de los parques se basa en su tamaño, entre otros criterios (Secretaría Distrital de Planeación, 2009), con este análisis se buscó no solo determinar la calidad de cada tipo de parque para prestar servicios ecosistémicos, sino también la variación del verdor que pueden presentar áreas con el mismo uso del suelo. Los parques metropolitanos presentaron un mayor verdor seguido por los vecinales, los zonales y los de bolsillo; esto significa que, en ese mismo orden, cada categoría tiene un mayor potencial para la oferta de servicios ecosistémicos por unidad de área. Por otra parte, la variación del verdor que existe para esta clase de uso del suelo es una evidencia de las carencias que tienen las estimaciones de servicios ecosistémicos basados en el uso del suelo a escalas pequeñas y en áreas altamente heterogéneas como las ciudades (Eigenbrod et al., 2010), ya que tienen en cuenta la variabilidad existente entre categorías de uso del suelo pero no la variabilidad al interior de cada una de ellas.
Los análisis de autocorrelación espacial indican que la distribución del estrato socioeconómico y las demás variables analizadas no es aleatoria. Esto se confirma con los resultados de los análisis Lisa, que evidencian claramente, primero, la ampliamente reconocida y marcada dualidad socioeconómica entre el norte y el sur de Bogotá (Muñoz y Ducón, 2016); y, segundo, un patrón diferencial en la distribución de la cobertura arbórea, el verdor de los parques y el potencial de regulación térmica. En pocas palabras, los hotspots y coldspots de estas variables se distribuyeron de una forma similar y refuerzan los resultados previamente analizados que indican la presencia de inequidad en la distribución del potencial de regulación de temperatura.
Si bien los resultados aportan conocimiento valioso para comprender la capacidad de regulación de temperatura en las diferentes zonas de la ciudad, la aproximación se basa en el análisis de la oferta. Un análisis más detallado que incluya la identificación de la demanda del servicio en las diferentes zonas de la ciudad en términos de la exposición y la vulnerabilidad a la temperatura, brindaría bases más sólidas para la implementación de acciones más efectivas (Norton et al., 2015). Otro aspecto a tener en cuenta sobre los resultados obtenidos es la fuente de los datos. Se ha identificado que existen diferencias importantes en la cuantificación de la infraestructura verde urbana y los servicios que presta, como consecuencia de las variaciones que presentan las diferentes fuentes de datos y las fallas institucionales en consolidar y manejar datos comunes, entorpeciendo así los procesos de planificación, gestión e investigación (Feltynowski et al., 2018). Bogotá tiene el privilegio de contar con el censo de arbolado urbano, el cual es una herramienta y una fuente de datos de gran valor en la búsqueda de una ciudad más sostenible, por su importancia para el estudio de los servicios ecosistémicos (Nowak et al., 2008; Escobedo et al., 2015). Esto recalca la importancia de contar con este tipo de datos medidos en terreno, lo más actualizados y depurados que sea posible. Aun así, es necesaria investigación adicional que contribuya al entendimiento de la prestación de este y otros servicios ecosistémicos en áreas urbanas en lo concerniente a especie, teniendo en cuenta las diferencias en los procesos fisiológicos y los rasgos funcionales, con el fin de determinar las especies más adecuadas para el bienestar del ser humano y la biodiversidad en las ciudades.
CONCLUSIONES
Se identificaron diferencias en la distribución de la cobertura arbórea y en el verdor de los parques públicos de acuerdo con el estrato socioeconómico de las UPZ. Esto es un indicador de inequidad en la distribución de estos elementos y tiene implicaciones en la provisión del servicio de regulación térmica y otros servicios ecosistémicos en las distintas zonas de la ciudad. Los análisis de autocorrelación espacial global, así como la localización de los hotspots y coldspots de las variables analizadas confirman estas diferencias, dejando en evidencia la marcada y ampliamente reconocida dualidad que existe entre el sur y el norte de la ciudad.
La heterogeneidad en la distribución de la infraestructura verde y el potencial de regulación térmica evidencian la complejidad de los procesos e interacciones que tienen lugar en la ciudad de Bogotá entre los diferentes componentes del ecosistema urbano. Esto constituye un desafío en la búsqueda de una distribución más equitativa de los beneficios que proveen las coberturas vegetales urbanas.
Un factor importante que puede estar motivando la inequidad observada es la historia del desarrollo de las diferentes zonas de la ciudad, que resulta en una configuración actual y determina la gama de posibilidades de manejo y mejora de la infraestructura verde. Lidiar con estas restricciones requiere encontrar soluciones diferentes a las tradicionales y que permitan mejorar el verdor y la calidad de vida de todos los ciudadanos.
Los parques públicos presentaron valores de verdor diferentes de acuerdo con su categoría de clasificación en la estructura ecológica principal. Esto sugiere que la capacidad de proveer servicios ecosistémicos cambia con el tipo de parque. Usualmente, en los estudios de servicios ecosistémicos y en los procesos de manejo y gestión urbana, se consideran las variaciones entre clases de usos del suelo, pero no la variabilidad al interior de estas clases. Considerar las diferencias que existen entre áreas verdes maximizaría los esfuerzos tanto en cuantificar más adecuadamente estos servicios como en garantizar una provisión más equitativa.
Este trabajo aporta información clave para el manejo de la infraestructura verde, la maximización del potencial de regulación térmica y la provisión de otros servicios ecosistémicos, ante escenarios inminentes de cambio climático.