La Zona Conurbada Colima-Villa de Álvarez (ZCCVA) se localiza en el estado de Colima, México, donde residen 265 770 habitantes sobre una superficie de 7 013.47 ha (Instituto Nacional de Estadística y Geografía [INEGI], 2020; 2023), siendo una de las 453 ciudades en el Sistema Urbano Nacional. Con un clima cálido subhúmedo, la ZCCVA alberga tres cauces fluviales relevantes: el río Colima y los arroyos Pereyra y El Manrique, con flujos norte-sur (figura 1).

Figura 1 Figura 1 Área de estudio Figura 1 Área de estudio Fuente: Elaboración propia con base en datos vectoriales de INEGI (2023).

La construcción del tejido ecológico intraurbano desarrollada en la presente investigación parte del reconocimiento de las microvarianzas espaciales, definitorias de la aptitud y potencial para la conectividad e integración entre elementos y áreas de alto valor ecológico-ambiental al interior de la ZCCVA. Para ello, el desarrollo metodológico incorporó la ponderación de criterios y escenarios diferenciados para el cálculo del tejido ecológico, así como factores de ajuste para evaluar jerarquías funcionales entre los vínculos generados, de acuerdo con su articulación y beneficios sobre el entorno urbano (figura 2).

Figura 2 Figura 2 Principios y contribuciones metodológicas Figura 2 Principios y contribuciones metodológicas Fuente: Elaboración propia.

La estructura de la investigación consideró tres etapas generales secuenciadas: 1) la identificación de áreas de alto valor ecológico-ambiental o de origen ecológico, derivados de la generación y procesamiento de datos relativos a la cobertura y uso del suelo urbano; 2) la construcción del tejido ecológico en escenarios diferenciados, a partir de la aplicación de un modelo gravitacional basado en la relación costo-eficiencia de cada vínculo y 3) el análisis de la articulación funcional y beneficios urbano-ambientales del tejido ecológico (figura 3).

Figura 3 Figura 3 Esquema metodológico Figura 3 Esquema metodológico Fuente: Elaboración propia.

Las áreas de origen ecológico representan nodos fundamentales para la articulación, integridad y seguridad ecosistémica (Huang et al., 2021; Kirk et al., 2023; Zhang et al., 2024); tienen la función de dar soporte y conducto a procesos bióticos. Para su identificación, se utilizaron bases de datos relativas a la cobertura y uso del suelo (LULC) en la ZCCVA, considerando la categorización CORINE de Esri (2024), el Índice de Vegetación de Diferencia Normalizada (NDVI) de la European Space Agency (ESA, 2024) y un mapeo poligonal elaborado con base en el Marco Geoestadístico Nacional del INEGI (2023) (tabla 1).

Tabla 1 Tabla 1 Clasificación de cobertura; fuentes de datos Tabla 1 Clasificación de cobertura; fuentes de datos

Tipo de datos	Descripción	Fuente
Clasificación de cobertura del suelo	Clasificación CORINE (ráster) de la cobertura y uso del suelo (LULC)	Esri (2024)
Índice de Vegetación de Diferencia Normalizada (NDVI)	Índice de Vegetación de Diferencia Normalizada con resolución de 15 metros	ESA (2024)
Mapeo poligonal de cobertura del suelo	Clasificación manual del territorio urbano según categorías de cobertura y uso del suelo	Elaboración de los autores con base en MGN INEGI (2023)

Fuente: Elaboración propia.

El principio básico contemplado deriva del potencial localizado para la conectividad, de acuerdo con: a) la permeabilidad del suelo; b) la extensión y continuidad de los parches permeables; c) la protección y destino del suelo urbano. Los datos LULC posibilitaron la detección de áreas generales de aplicabilidad, en donde se clasificaron zonas de acuerdo con su uso público o privado como primer filtro. El segundo paso fue el cálculo del NDVI (ecuación 1).

Ecuación 1

Dónde:

B4: Banda 4 Sentinel-2: Rojo.

B8: Banda 8 Sentinel-2: Infrarrojo cercano (NIR).

Para la detección de la aptitud del suelo urbano y la conectividad ecológica, se definió que los valores ≥0.16 del NDVI representan parches de valor ecológico-ambiental. Dichos pixeles ráster se utilizaron para el Análisis Espacial de Patrones Morfológicos (MSPA) (Soille y Vogt, 2008; Vogt et al., 2022), para determinar la conectividad y significancia de las áreas por encima del umbral de arborización. Dicho proceso permitió definir un área mínima de cobertura para seleccionar nodos significantes en la configuración del tejido ecológico (núcleos). Finalmente, se realizó un mapeo poligonal de la ZCCVA, de acuerdo con la cobertura, uso y protección de cada elemento urbano.

El cálculo de las trayectorias más eficientes para la vinculación espacial del tejido ecológico intraurbano se derivó de la herramienta Gravity Model Toolbox en ArcMap 10.8 (Wanghe et al., 2020), tomando la base poligonal sobre la cobertura del suelo como datos de entrada (input). La definición del peso, impedancia y costo para la conectividad, se establecieron con base en dos escenarios, en donde se adaptaron los valores propuestos por Kong et al. (2010) y Wanghe et al., (2020) al contexto urbano de la ZCCVA. El primer escenario otorga una preponderancia a la conectividad espacial por medio de la infraestructura verde urbana, principalmente en camellones arbolados, mientras que el segundo promueve la conectividad por medio de los cauces fluviales locales; asimismo, se definió un rango de valores (mínimos/máximos) ponderados por el NDVI promedio sobre cada polígono, considerando valores máximos y mínimos en cada clasificación (tabla 2).

Tabla 2 Tabla 2 Parametrización de valores de cobertura según escenarios Tabla 2 Parametrización de valores de cobertura según escenarios

Componente	Peso	Escenario 1		Escenario 2
		Impedancia	Costo	Impedancia	Costo
Espacios públicos	0.66	1.00-10.00	47.20-73.60	1.00-10.00	47.20-73.60
Áreas residuales	0.40	0.20-0.50	68.00-84.00	0.20-0.500	68.00-84.00
Camellones	0.10	0.125-0.50	47.20-60.40	0.125-0.50	47.20-60.40
Cauces fluviales	0.10	1.00	99.00	0.10	20.00
Manzanas urbanas	0.14	1000.00	1000.00	1000.00	1000.00
Vialidades	0.14	60.00-100.00	86.00-88.00	60.00-100.00	86.00-88.00

Fuente: Elaboración propia.

Las ecuaciones 2.1, 2.2 y 2.3 definen el sustento estadístico aplicado por el modelo gravitacional, ajustadas por Wanghe et al. (2020) a partir del trabajo de Kong et al. (2010). Su aplicación en el área de estudio consideró una microescala para el procesamiento unitario básico, con una precisión de 2.5 metros [1.4*10^-6 grados] para el cálculo y definición de vínculos en el tejido ecológico intraurbano. Con base en los datos [multipuntos] resultantes, se definieron segmentos vectoriales a mediante los cuales se evaluó la articulación funcional de cada escenario y componente.

Ecuación 2.1

Ecuación 2.2

Ecuación 2.3

Dónde:

G_ab = Interacción entre los nodos a-b

(a,b) = Nodos a y b (origen ecológico-ambiental)

N_a/b = Peso adscrito al nodo a o b

S_a/b = Área adscrita al nodo a o b

P_af/bf = Impedancia/resistencia del nodo a o b

L_ab = Impedancia acumulada entre los nodos a-b

De manera complementaria, se utilizaron datos resultantes de la impedancia del suelo urbano, para la conectividad en relación con el análisis de la concentración espacial de los trayectos más eficientes (McRae et al., 2008) Se consideró el parámetro de corriente acumulada obtenido en Circuitscape (McRae et al., 2016), con el propósito de la definición de la preponderancia y jerarquía de tramos/nodos viales que resulten críticos en la conectividad.

La definición de los vínculos con mayor articulación funcional en el entorno urbano se definió por medio de la ponderación y ajuste geoestadístico de acuerdo con los servicios/beneficios brindados por el tejido ecológico, considerando indicadores de la población beneficiada, su conectividad comercial e incidencia sobre el fortalecimiento de la red fluvial local (tabla 3 y figura 4).

Tabla 3 Tabla 3 Vinculación funcional del tejido ecológico; fuentes de datos Tabla 3 Vinculación funcional del tejido ecológico; fuentes de datos

Tipo de datos		Descripción	Fuente
1	Densidad poblacional	Población total por manzana urbana	Censo de Población y Vivienda INEGI (2020)
2	Conectividad comercial	Unidades económicas en el entorno urbano	Directorio Estadístico Nacional de Unidades Económicas (DENUE) INEGI (2024)
3	Conectividad fluvial	Centro geométrico (vectores) de cauces fluviales locales	Marco Geoestadístico Nacional (MGN) INEGI (2023)

Fuente: Elaboración propia.

Figura 4 Figura 4 Distribución de indicadores socioambientales Figura 4 Distribución de indicadores socioambientales Fuente: Elaboración propia con base en datos de INEGI (2020; 2023; 2024)

Su cálculo se derivó del desempeño diferenciado de los vínculos que configuran el tejido ecológico en cada escenario, en donde se consideró un área de influencia efectiva de 400 metros para el procesamiento geoestadístico de los indicadores contemplados, mismo que se define por las ecuaciones 3.1 y 3.2.

Ecuación 3.1

Ecuación 3.2

Dónde:

AF_{soc,eco,amb i} = Articulación funcional (social/económica/ambiental) del vínculo i.

UG_Pol(Pob)= Población total en la Unidad Geoestadística poligonal (manzanas).

UG_Pun = Unidades Geoestadísticas puntuales (Unidades económicas, puntos a lo largo de cauces fluviales)

UϵA(V; R) _i = Unidades Geoestadísticas U dentro del área A definida por un vector V y un radio R del vínculo i.

Los valores resultantes fueron normalizados (0-1) para la integración geoestadística de la articulación funcional de cada vínculo, tomando los datos del modelo gravitacional, incorporando un ajuste ponderado por la distancia. Esto posibilitó la identificación de la estructura jerárquica del tejido ecológico en la ZCCVA; se define por las ecuaciones 4.1 y 4.2.

Ecuación 4.1

Ecuación 4.2

Dónde:

AF ^norm = Valores normalizados de la articulación funcional (social/económica/ambiental) del vínculo i.

G_ab = Interacción entre los nodos a-b

G_ab(int) _i = Interacción integral entre los nodos a-b

La tabla 4 detalla los puntos de origen ecológico-ambiental identificados en la ZCCVA, resultado del análisis del uso y cobertura del suelo, junto con la aplicación del MSPA basada en el NDVI (figura 5 y tabla 5). Para la definición de escenarios, el criterio principal fue el uso público del suelo urbano y la protección institucional de las áreas, según su estatus —zonas arqueológicas, parques ecológicos, cauces fluviales—, factores que determinan su vinculación y viabilidad como nodos del tejido ecológico. Se encontró que solo el 4.69 % del área de estudio (7 019.49 ha) es apta como zona de origen, destacando “La Campana” (112.85 ha) como la mayor reserva ecológica.

Tabla 4 Tabla 4 Puntos de origen ecológico-ambiental Tabla 4 Puntos de origen ecológico-ambiental

Parque Arqueológico La Campana	1	112.85	34.08 %	Federal
Zona Arqueológica El Chanal	1	3.60	1.09 %	Federal
Parque Regional Metropolitano Griselda Álvarez	2	21.39	6.46 %	Estatal
Parque Huertas del Cura Arzac	2	1.31	0.40 %	Municipal
Parque El Rodeo	1	12.16	3.67 %	Estatal
Parque Ecológico El Topocharco	1	1.33	0.40 %	Municipal
Ex Zona Militar del municipio de Colima	1	14.90	4.50 %	Estatal
Parque Hidalgo	1	5.60	1.69 %	Municipal
Jardín Núñez	1	1.99	0.60 %	Estatal
Arroyo Pereyra	3	25.61	7.73 %	Federal
Río Colima	5	44.90	13.56 %	Federal
Arroyo El Manrique	2	29.32	8.86 %	Federal
Cauces fluviales secundarios	11	56.17	16.96 %	Federal

Fuente: Elaboración propia. Nota

Porcentaje del total de áreas con valor ecológico-ambiental en la ZCCVA (297.67ha).

Tabla 5 Tabla 5 Síntesis estadística: cobertura del suelo urbano según fuentes Tabla 5 Síntesis estadística: cobertura del suelo urbano según fuentes

Clasificación CORINE de la cobertura del suelo	Arbolado	7 620	76.2	1.08 %
	Cultivos	18 179	181.79	2.58 %
	Área construida	620 169	6 201.69	87.98 %
	Suelo desnudo	568	5.68	0.08 %
	Pastizales	58 383	583.83	8.28 %
NDVI (Sentinel-2)	Área permeable [>0.16]	384 256	3 842.56	54.51 %
	Entorno urbano [<0.16]	320 663	3 206.63	45.49 %
MSPA	Núcleos	222 279	2 222.79	31.53 %
	Bordes	47 581	475.81	6.75 %
	Islas	20 463	204.63	2.90 %
	Puentes	21 453	214.53	3.04 %
	Ramificación	31 988	319.88	4.54 %
	Circuitos	8 286	82.86	1.18 %
	Perforaciones	32 206	322.06	4.57 %
	Fondo	320 663	3 206.63	45.49 %
Mapeo poligonal de la cobertura del suelo	Manzanas urbanas	6 699	3 951.65	56.06 %
	Espacios públicos	451	174.826	2.48 %
	Áreas naturales	11	175.13	2.48 %
	Áreas permeables	389	1 054.1654	14.95 %
	Cauces fluviales	135	156.003	2.21 %
	Camellones	1 220	97.9856	1.39 %
	Vialidades	7 622	1 439.43	20.42 %

Fuente: Elaboración propia

Figura 5 Figura 5 Coberturas y usos del suelo urbano Figura 5 Coberturas y usos del suelo urbano Fuente: Elaboración propia con base en datos de Esri (2024), ESA (2024) e INEGI (2023).

A pesar de la variabilidad en los índices de permeabilidad del suelo, que oscilan entre el 12.02 % (CORINE), el 54.51 % (NDVI), el 31.51 % (MSPA) y el 23.52 % (mapeo poligonal), la fragmentación de los parches arborizados en el área de estudio señala un aislamiento espacial y funcional del tejido ecológico en la ZCCVA. Además, se detecta un potencial de conectividad del 7.58 % en la morfología (534.41 ha, correspondiente a puentes y ramificaciones según el MSPA) y del 3.87 % en la infraestructura verde de la ZCCVA (272.81 ha, incluyendo espacios públicos y camellones arbolados). Estas áreas fungen, parcial o totalmente, como vínculos espaciales o corredores verdes dentro del tejido ecológico.

La figura 6 muestra la ubicación de los puntos de origen ecológico, distinguiéndolos por su carácter urbano o fluvial. A partir del análisis de la impedancia del suelo urbano para la conectividad ecológica, basado en los valores de la tabla 2, se identificaron los tramos viales y camellones arbolados en la ZCCVA, que son críticos para determinar la articulación espacial y la jerarquía de la infraestructura vial y los corredores verdes. Este análisis incluyó una muestra del 5.75 % de los polígonos (439 de los 7 622 tramos) con la mayor corriente acumulada.

Figura 6 Figura 6 Parámetros para el modelado de la conectividad ecológica intraurbana Figura 6 Parámetros para el modelado de la conectividad ecológica intraurbana Fuente: Elaboración propia.

La interacción espacial entre dos puntos de origen ecológico (G_ab), derivada del modelo gravitacional, se calculó en tres escenarios:

Puntos formales, con preponderancia vial en la parametrización inicial.

Los resultados del modelo, basados en su desempeño sistémico, se ilustran en la figura 7. En esta, se observa una varianza significativa en el trazo de los corredores generados. Específicamente, el escenario A se estructura a partir de los camellones arbolados, mientras que los escenarios B y C toman el río Colima como eje rector para la vinculación espacial.

Figura 7 Figura 7 Resultados del modelo gravitacional Figura 7 Resultados del modelo gravitacional Fuente: Elaboración propia con base en resultados de Gravity Model Toolbox (Wanghe et al. 2020).

La vinculación funcional de los nexos en cada escenario, calculada mediante las ecuaciones 3.1 y 3.2, se ilustra en la figura 8. Se observa una tendencia de centralización espacial en los indicadores socioeconómicos (1 y 3), asociada a los mayores índices de consolidación en el núcleo urbano en comparación con las periferias. En cuanto a la vinculación de los vínculos con la red fluvial, se destaca la preponderancia del Río Colima y el Arroyo Pereyra como nexos efectivos entre las zonas norte y sur de la ZCCVA. Además, se registra un aislamiento funcional de los sectores sur y poniente.

Figura 8 Figura 8 Vinculación funcional del tejido ecológico según escenarios Figura 8 Vinculación funcional del tejido ecológico según escenarios Fuente: Elaboración propia

La ponderación integrada del tejido ecológico, calculada mediante las ecuaciones 4.1 y 4.2, se ilustra en la figura 9, detallando el desempeño de los extremos estadísticos (<10%, >90%) en la tabla 6. Los datos demuestran una varianza significativa en la cobertura territorial de los vínculos con mayor desempeño, así como en la estructura y geometría de sus componentes (G_ab e indicadores 1-3). Los vínculos más eficientes del escenario A presentan una mayor articulación socioeconómica, con un incremento del 36.90 % en la población cubierta, en vínculos 54.39 % más cortos, mientras que los escenarios B y C incorporan una longitud 48.33 % mayor de red fluvial local.

Figura 9 Figura 9 Escenarios integrados: vinculación funcional del tejido ecológico Figura 9 Escenarios integrados: vinculación funcional del tejido ecológico Fuente: Elaboración propia.

Tabla 6 Tabla 6 Comparativa estadística: desempeño funcional de corredores ecológicos Tabla 6 Comparativa estadística: desempeño funcional de corredores ecológicos

Escenario		Gab(Int) (#)	Unidades económicas (#)	Población (habitantes)	Cauces fluviales (km)	Longitud (km)
A	Máx. >90 %	286.59	2 168.67	29 884.70	3.02	3.21
	Min. <10 %	0.04	1 512.00	6 755.17	1.28	0.59
B	Máx. >90 %	240.36	2 088.67	18 858.50	6.25	5.89
	Min. <10 %	0.02	1 567.33	7 483.17	1.58	0.69
C	Máx. >90 %	451.42	2 137.70	20 059.70	6.14	5.63
	Min. <10 %	0.66	381.70	7 247.30	1.80	0.92

Fuente: Elaboración propia.

El análisis del 10 % más eficiente de la influencia territorial-urbana, según cada escenario, se detalla en la tabla 7. En esta tabla, el escenario A destaca como el más eficiente en términos de población beneficiada, mostrando un incremento del 25.12 % (19 905 habitantes) al 35.43 % (28 077 habitantes) en comparación con los escenarios B y C. Además, se observa una mayor extensión territorial en el escenario A, con diferencias del 24.99 % (300.11 ha) al 29.80 % (357.85 ha) en la misma comparación. Estas condiciones evidencian una clara diferenciación en la consolidación del tejido ecológico intraurbano: corredores verdes en el escenario A y recuperación de cauces fluviales en los escenarios B y C.

Tabla 7 Tabla 7 Cobertura poblacional y territorial de vínculos con mayor desempeño Tabla 7 Cobertura poblacional y territorial de vínculos con mayor desempeño

Escenario		Cobertura poblacional		Cobertura territorial
		(habitantes)	(% área de estudio)	(ha)	(% área de estudio)
A	Valores máximos (>90%)	79 229	29.81 %	1 200.72	17.12 %
B		59 324	22.32 %	900.601	12.84 %
C		51 152	19.24 %	842.864	12.02 %

Fuente: Elaboración propia.

Esta investigación abordó la construcción de escenarios para la conectividad y eficiencia del tejido ecológico en la ZCCVA, considerando criterios y parámetros para su optimización estructural a microescala, así como para ponderar los beneficios socioambientales y su cobertura territorial-urbana. La metodología propuesta contribuye al reconocimiento de la función integral de los servicios ecosistémicos derivados de la consolidación, el fortalecimiento o la preservación de los vínculos que configuran el tejido ecológico intraurbano. Por lo tanto, resulta de gran utilidad y aplicación para la gestión de los recursos ambientales, así como para la definición de criterios en la planificación de subsistemas de infraestructura verde.

La identificación de zonas de origen para la generación de escenarios diferenciados, mediante el uso de datos del NDVI y la aplicación del MSPA, son técnicas ampliamente utilizadas en los estudios de conectividad ecológica (Jiang et al., 2024; Hu et al., 2022; Zhang et al., 2024). No obstante, para garantizar la precisión del modelo en sistemas urbanos, resulta necesario incorporar criterios de cobertura y uso del suelo (LULC) (He et al., 2023; Jamil et al., 2024; Wanghe et al., 2020). Consecuentemente, el presente estudio consideró micro-varianzas espaciales (≤3 metros) para el mapeo poligonal requerido por el modelo gravitacional (figura 2), mismo que permitió un alto grado de precisión y fiabilidad en el cálculo de cada uno de los vínculos en los escenarios contemplados.

El reconocimiento de los elementos urbanos que estructuran el tejido ecológico, a través de corredores verdes urbanos (como camellones arbolados y espacios públicos), es de especial importancia para mitigar el aislamiento y fragmentación del paisaje ecológico urbano y sus áreas de origen o conservación ambiental (Huang et al., 2021; Jamil et al., 2024). Por lo tanto, es indispensable garantizar su continuidad y eficiencia a través del territorio urbano, con la finalidad de optimizar la integración de zonas y sectores. Esto contrasta con lo registrado en la ZCCVA (figura 9.C), donde se observa una desvinculación derivada de la deficiencia estructural/espacial de la infraestructura verde en patrones de urbanización con alta densidad y carácter periférico; ello, denota una desigualdad distributiva de los servicios ecosistémicos a través del territorio urbano.

La literatura del análisis del paisaje y la ecología urbana ha registrado esta inequidad en la distribución, acceso y calidad de las áreas verdes de significancia ecológico-ambiental. Principalmente, esto es atribuible a particularidades de la geografía humana, donde las ciudades, en sus diversas escalas, se disponen y adaptan al territorio en una diversidad de patrones y relaciones con respecto al suelo de conservación o protección ambiental.

Los trabajos de Hu et al. (2022), Huang et al. (2021), Shi et al. (2021) y Wanghe et al. (2020) coinciden en la preponderancia jerárquica de los cauces fluviales en la configuración ecológica del paisaje, así como en su fragmentación derivada de la producción intensiva del suelo urbano. En el caso de la ZCCVA, se observa una centralización en las áreas de origen ecológico, así como una marcada asimetría en la presencia y continuidad del sistema fluvial local, con mayor cobertura en el Municipio de Colima (figuras 4.1 y 5). Esto indica la necesidad de acciones compensatorias para el ordenamiento ecológico del territorio, enfocadas en la generación de áreas de reserva y conservación como recursos clave para la sustentabilidad en la planeación y el desarrollo urbano a mediano y largo plazo.

La definición formal de la estructura espacial de subsistemas ecológicos y del paisaje natural permitió demostrar la relevancia de la red fluvial, principalmente del río Colima y el arroyo Pereyra, como los elementos naturales con mayor jerarquía y centralidad en el territorio urbano. Dentro de los estudios de la ecología del paisaje, los cauces fluviales son preponderantes en estrategias para la protección y restauración de la biodiversidad, así como en la provisión de servicios ecosistémicos y el sostenimiento de los asentamientos humanos (Jiang et al., 2024; Li et al., 2024; Staccione et al., 2022). Pese a ello, en la ZCCVA se observan altos índices de fragmentación, aunados a la discontinuidad y la reducción en las dimensiones de sus cauces fluviales; lo cual, torna necesario plantear la importancia de estrategias para la recuperación, saneamiento y renaturalización de los ecosistemas riparios urbanos para fortalecer la salubridad medioambiental, promover la biodiversidad y mejorar los servicios ecosistémicos asociados a las redes fluviales (Kirk et al., 2023; Liu et al., 2024; Shi et al., 2021).

Asimismo, la investigación destaca la relevancia de articular y complementar las reservas ecológicas del paisaje urbano respecto a dos subsistemas de origen antrópico: el conjunto de espacios públicos y la infraestructura verde vial. Si bien, las características inherentes al origen, las dimensiones y la función de dichos subsistemas imposibilitan la consideración de sus componentes como áreas de origen ecológico, su disposición continua e integrada a través del territorio urbano posibilita el fortalecimiento de la estructura y funcionalidad del tejido ecológico. En ese sentido, la optimización de los subsistemas que dan sostenimiento y conducción a los flujos naturales y humanos establece parámetros para el ordenamiento y gestión sustentable del paisaje urbano, así como la profundización de la relación humano-naturaleza en las ciudades (He et al., 2023; Huang et al., 2021; Jato-Espino et al., 2023; Jiang et al., 2024; Veerkamp et al., 2021).

Las implicaciones sobre las prácticas de planificación urbana e implementación de políticas públicas asociadas a gestionar los recursos ecológicos en contacto con los asentamientos humanos, poseen dos directrices relacionadas con el método desarrollado y resultados descritos. Primeramente, el pleno reconocimiento de la interdependencia entre sistemas ecológicos para la salubridad y seguridad medioambiental, y su vinculación efectiva con la estructura socioeconómica de las ciudades (Jiang et al., 2024; Zhang et al., 2024). Por otro lado, la incorporación efectiva de parámetros para garantizar la equidad social del acceso y contacto con espacios verdes urbanos, con el objetivo de mejorar el bienestar psicofísico de sus poblaciones (He et al., 2023; Syamili et al., 2023).