
==== Front
Landsc Ecol
Landsc Ecol
Landscape Ecology
0921-2973
1572-9761
Springer Netherlands Dordrecht

1911
10.1007/s10980-024-01911-2
Research Article
New criteria for sustainable land use planning of metropolitan green infrastructures in the tropical Andes
LaRota-Aguilera María José mariajose.larota@uab.cat
mj.larota.aguilera@ub.edu

1
Zapata-Caldas Emmanuel 2
Buitrago-Bermúdez Oscar 3
Marull Joan 4
1 Institute of Regional and Metropolitan Studies of Barcelona, Barcelona, Spain
2 grid.8271.c 0000 0001 2295 7397 Centro Internacional de Agricultura Tropical (CIAT), University of Valle (Primary), Valle del Cauca, Colombia
3 https://ror.org/00jb9vg53 grid.8271.c 0000 0001 2295 7397 University of Valle, Valle del Cauca, Colombia
4 Minuartia, Barcelona, Spain
28 5 2024
28 5 2024
2024
39 6 1129 9 2023
21 5 2024
© The Author(s) 2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Context

Urbanization is rapidly increasing worldwide, with about 60% of the global population currently residing in cities and expected to reach 68% by 2050. In Latin America's tropical Andes region, managing these changes poses challenges, including biodiversity loss and vulnerability to climate change.

Objectives

This study assesses urban growth and agricultural intensification impacts on the ecological functionality of metropolitan green infrastructures and their capacity to provide ecosystem services using a landscape sustainability and sociometabolic approach. Specifically, it aims to identify landscape configurations promoting socio-ecological sustainability amidst rapid urbanization.

Methods

A landscape-metabolic model (IDC) was applied to evaluate the interactions between land use changes and ecosystem functions in the metropolitan region of Cali.

Results

Agricultural intensification and industrialization, coupled with uncontrolled urban growth, have significantly transformed the landscape, posing threats to its sustainability. The prevailing biocultural landscapes hold a substantial potential to provide essential ecosystem services to the metropolis. The IDC offers an approach that utilizes a land cover map and agricultural production/metabolism data to calculate an indicator closely related to ecosystem services and multifunctionality.

Conclusions

The IDC model stands out for efficiently capturing landscape dynamics, providing insights into landscape configuration and social metabolism without extensive resource requirements. This research highlights the importance of adopting a landscape-metabolic and green infrastructure framework to guide territorial policies in the tropical Andes and similar regions. It stresses the need for informed land use planning to address challenges and leverage opportunities presented by biocultural landscapes for regional sustainability amidst rapid urbanization and agricultural expansion.

Keywords

Integrated socio-ecological assessment
Social metabolism
Ecosystem services
Biocultural landscapes
Landscape sustainability
Agroforestry mosaics
http://dx.doi.org/10.13039/501100000155 Social Sciences and Humanities Research Council of Canada SSHRC 895-2011-1020 LaRota-Aguilera María José http://dx.doi.org/10.13039/501100009543 Pontificia Universidad Javeriana http://dx.doi.org/10.13039/501100007329 Universidad del Valle 4339 4339 Zapata-Caldas Emmanuel Buitrago-Bermúdez Oscar issue-copyright-statement© Springer Nature B.V. 2024
==== Body
pmcIntroduction

Urban centres and their interactions with peri-urban and rural areas have historically been pivotal in society's development, with farmland and natural ecosystems supplying essential resources (Steel 2008). However, the Industrial Revolution introduced two major socio-ecological changes: urbanization and agricultural intensification (Krausmann et al. 2008; Swyngedouw and Heynen 2003). Coupled with current economic and technological models, population growth has led to increased resource consumption in cities, impacting water, energy, materials, and food (Balatsky et al. 2015; Krausmann et al. 2009). This transformation from natural and low-intensity agriculture to urban and high-intensity agriculture has resulted in biodiversity and ecosystem service losses (Elmqvist et al. 2013), disrupted water flows (Hibbs and Sharp Jr 2012), and escalated global greenhouse gas emissions (IPCC 2014; Vermeulen et al. 2012).

At local levels, urbanization and agricultural intensification affect ecosystems' ability to provide essential services and support life (Brondizio et al. 2019). Moreover, under current climate change scenarios, biodiversity and ecosystem service loss heighten vulnerability to climate impacts (Burak Güneralp et al. 2013; IPCC 2014; McDonald et al. 2013). Despite this reality, the global population is projected to reach 9.8 billion by 2050, with 68% residing in urban areas (United Nations 2019). The impacts of this trend will disproportionately affect low- and lower-middle-income regions, where rapid urbanization is expected, potentially straining resilience and adaptation capacity (IPCC 2014).

This complex scenario presents unprecedented challenges for urban–rural relationships and land-use planning in metropolitan regions, requiring comprehensive and integrated approaches to find solutions (Yacamán-Ochoa et al. 2020). In this sense, the concept of green infrastructure has gained importance in understanding the role of metropolitan open spaces in sustaining urban environments (Benedict and McMahon 2002). However, gaps persist in how to effectively implement green infrastructure for biodiversity conservation and ecosystem service provision (Chatzimentor et al. 2020; Demuzere et al. 2014; Vásquez et al. 2019).

Notably, most advances in green infrastructure have occurred in the Global North (Chatzimentor et al. 2020; European Commission 2013; Slätmo et al. 2019), while the Global South, including Latin America's tropical Andes, lags (Pauleit et al. 2021). In these regions, land-use changes driven by agricultural intensification and resource exploitation pose unique challenges. The resulting rural-to-urban migration is exacerbated by political instability, corruption, and armed conflicts, leading to disorganized metropolitanization and inadequate urban planning (Aide et al. 2019; Canales and Canales Cerón, 2013; Angotti 1996).

Therefore, tropical Andean countries face a critical choice: sustain their natural ecosystems or pursue current economic development models. The metropolitan region of Cali, Colombia, exemplifies this dilemma, experiencing significant socioeconomic and land-use changes driven by agro-industrial practices and rural migration (Delgadillo-Vargas 2014; Marull et al. 2017; Martínez-Toro and Patiño-Gómez 2015; Centro Nacional de Memoria Histórica 2014; Uribe-Castro 2017).

To improve landscape sustainability and address contemporary challenges, it's crucial to understand how landscape patterns influence ecosystem services and human wellbeing (Wu 2006, 2013). We propose that this endeavor demands for integrated and transdisciplinary approaches that explore the relationship between different landscape-metabolic configurations within metropolitan regions and their capacity to provide ecosystem services.

This article employs a landscape-metabolism model (Marull et al. 2019, 2018) to comprehensively evaluate the ecological functionality of Cali's metropolitan green infrastructure. The study pursues three main objectives: assessing the impacts of agricultural intensification and urbanization on the green infrastructure's ecological functions and services, analysing the relationship between these configurations and their ability to supply ecosystem services, and guiding future land use policies for a functional metropolitan green infrastructure in the region. By adopting a landscape sustainability perspective, this study proposes new criteria for sustainable land use planning in metropolitan green infrastructures, ensuring that these regions can meet present and future needs without compromising the integrity of their ecosystems.

In the following sections, we provide a brief context for the case study, outline the methodologies employed, present the landscape-metabolic assessment results, and discuss the emerging opportunities and challenges posed by agricultural landscapes for the sustainability of Cali's metropolitan region. We then assess the relevance of our findings and the potential of adopting a green infrastructure framework to shape future land use policies in the tropical Andes.

Methodology

Case study

Currently, the metropolitan region of Cali is not considered an official administrative entity. There have been some initiatives to formally consolidate it, mainly based on capital criteria dictated by the predominant role of agribusiness in the region; however, they have not thrived (Martínez-Toro and Patiño-Gómez 2015; Urrea-Giraldo and Candelo-Álvarez 2017). Therefore, this study considers the Upper Cauca River Valley (henceforth UCRV) as the territorial reference for the metropolitan region and the study area.1 The UCRV limits were defined based on three criteria: i) the third metropolitan crown (Martínez-Toro and Patiño-Gómez 2015); ii) the limits of the hydrographic sub-basins in which the urban centres of the third metropolitan crown are located; and iii) the areas of influence of sugarcane cultivation, given their economic importance for the region (Fig. 1). The population in the UCRV is approximately 3′635.573 people, and 2′172,527 people live in the urban area of the city district of Cali (DANE 2018).Fig. 1 Location, administrative limits, and land cover map of the metropolitan region of Cali. Source: Corine Land Cover Map for Colombia 2012 (level 3 legend)

The UCRV is located within the geographic valley of the Cauca River (448,000 ha), between the Central and Western Mountain ranges of the Andes, in southwest Colombia (Fig. 1). The study area occupies 1′004,000 hectares and includes elevations between 800 and 5,000 m above sea level. Originally, this biogeographic unit was made up of various ecosystems, including tropical dry forests, tropical pastures, high Andean forests, and paramos. However, currently, 19.2% of the UCRV is covered by sugar cane monocultures (192,742 ha); 31,1% by mixed crops that include coffee, banana, and fruit trees, in mosaics with natural areas and pastures (312,682 ha); 13% by natural and planted open pastures (130,547 ha); 32.5% by forests and other natural areas -among which are the last enclave of tropical dry forest and paramos- (198,659 ha); and 2.4% by urban and suburban spaces built (23,954 ha) (Fig. 1).

In the UCRV, two key land cover changes have occurred: i) the shift from clean pastures, cropland (including sugarcane), and seminatural transition areas to low-density urban zones, and ii) the expansion of sugarcane crops across the flat region, previously a mosaic of pastures and mixed crops, dating back to the mid-nineteenth century, profoundly altering the landscape (Delgadillo-Vargas 2014). Cali's rapid growth is influenced by three factors: its designation as the Valle del Cauca department capital in 1910, concentrating public services; the growth of its industrial sector (Vásquez 2018); and over five decades of violence and internal conflict, making the city a major destination for internal refugees (140,751 individuals between 1985 and 2014; Centro Nacional de Memoria Histórica 2015). This case study is essential for understanding these complex dynamics.

Landscape ecological metrics

The landscape's functional structure evaluation was based on the adaptation of the Corine Land Cover map of land cover at a scale of 1:100,000 for Colombia, considering level 3 of the legend (Fig. 1). The reference year is 2015 since it is the most recent year with relevant information. This map was updated with supervised classification and interpretation of aerial photographs for 2013 and 2016 by the authors to correct areas with high cloudiness. The final coverage legend includes 22 classes and is shown in Table 1.Table 1 Land cover and surface distribution in the metropolitan region of Cali

Typology	Land cover	Reclassification	Area (ha)	%	
Forest and seminatural areas	Rocky outcrops	Others (rocky outcrops, glaciers)	257.81	0.0%	
Natural sandy areas	Others (rocky outcrops, glaciers)	287.51	0.0%	
Glacial and snow zones	Others (rocky outcrops, glaciers)	66.12	0.0%	
	Total others	611.44	0.1%	
Natural shrublands	Natural shrublands	39627.44	3.9%	
Natural grasslands (Paramo)	Natural grasslands (Paramo)	87665.07	8.7%	
Dense forest	Dense forest	117687.3	11.7%	
Riparian forest	Riparian forest	2487.41	0.2%	
	Total primary forest	120174.7	12.0%	
Fragmented forests	Secondary forest	29330.82	2.9%	
Secondary vegetation or in transition	Secondary vegetation or in transition	48542.91	4.8%	
	Total secondary forest	77873.73	7.8%	
		Total forest and seminatural areas		32.5%	
Agricultural land	Permanent woody crops	Permanent woody crops	3635.33	0.4%	
Permanent herbaceous crops	Sugarcane plantations[1]	192742.2	19.2%	
	Total permanent crops	196377.5	19.6%	
Pasture and forest mosaic	Pasture and forest mosaic	75165.72	7.5%	
Cropland mosaic	Cropland mosaic	1922.57	0.2%	
Cropland and forest mosaic	Cropland and forest mosaic	18478.08	1.8%	
Cropland, forest, and pasture mosaic	Cropland, forest, and pasture mosaic	152551.3	15.2%	
Cropland and pasture mosaic	Cropland and pasture mosaic	64564.77	6.4%	
	Total agricultural mosaics	312682.5	31.1%	
Wooded pastures	Wooded pastures	4673.74	0.5%	
Weeded pastures	Weeded pastures	41797.32	4.2%	
Clean pastures	Clean pastures	84076.29	8.4%	
	Total pastures	130547.3	13.0%	
Timber plantations	Timber plantations	3238.40	0.3%	
		Total agricultural land		64.0%	
Bodies of water	Rivers (50 m)	Rivers and natural water bodies	2569.59	0.3%	
Artificial water bodies	Artificial water bodies	1637.38	0.2%	
Lagoons, lakes, swamps, and natural swamps	Rivers and natural water bodies	527.56	0.1%	
	Total water	4734.54	0.5%	
Build-up areas	Continuous urban fabric (Urban areas)	Urban areas	16228.67	1.6%	
Discontinuous urban fabric (Suburban areas)	Urban areas	3447.70	0.3%	
Industrial or commercial areas (Industrial area)	Urban areas	4278.26	0.4%	
	Total urban and industrial areas	23,954.63	2.4%	
Other	Nude and degraded soils	Nude soils	854.29	0.1%	
Total			1,004,000	100.0%	
Source: Own elaboration based on the Corine Land Cover Map for Colombia

The study area was divided into 502 analysis units (hexagons), each of 20 km2 (2,000 ha) (Fig. 1), for which four indicators were calculated: i) the Shannon Index ii) the Ecological Connectivity Index iii) the Human Appropriation of Net Primary Productivity; and iv) the Intermediate Disturbance Complexity (IDC) model. Each of these landscape metrics is described below.

Shannon index (H') (Shannon and Weaver 1948)

Characterises the landscape structure as a function of the land cover heterogeneityH=∑i=1j=kpj∙Logpi

For its calculation, eight land covers (j) considered potential habitats for biodiversity (forest, shrubland, grassland, heterogeneous crops, sugarcane, pastures, bodies of water, rocky outcrops, sandy areas) and one "no habitat" category grouping urban and industrial areas, degraded land, and road infrastructure, were defined. Thus, Pji is the proportion of land cover j in hexagon i. The H' values range from 0 to 1, with 0 being a homogeneous landscape with a single predominant cover and 1 being a theoretical landscape with many land cover classes distributed equitably.

Ecological connectivity index (ECI) (Marull and Mallarach 2005)

Assesses the functionality of the landscape in terms of the ecological connectivity between related land covers (Pino and Marull 2012). ECI was calculated through a cost-distance model based on an affinity matrix considering 7 types of 'functional ecological areas' (i.e., forest, shrubland, grassland or paramo vegetation, agroforestry mosaics, crops, pastures, and sugar cane), and an impact matrix considering the 'anthropogenic barriers' (i.e., urban areas, infrastructures). The selection criteria for functional ecological areas, the coefficients, and the type of anthropogenic barriers were obtained from a review of the literature and local experts' knowledge. Ecological connectivity is calculated for each of the different functional ecological areas:ECIb=10-9In1+Xi/In1+Xt3

where Xi is the value of the sum of the cost distance per pixel and Xt is the maximum theoretical cost distance.

The total ecological connectivity values are calculated from the values obtained for each type of functional ecological area:ECI=∑ECIb/m

where m is the absolute number of functional ecological areas considered. The highest values, in a range of 0 to 10, represent high ecological connectivity.

Human appropriation of net primary production (HANPP)

Measures the disturbance exerted by society on a particular ecosystem (Haberl et al. 2007) as a function of the degree to which humans modify the amount of NPP available to other species, fundamentally through two processes: the removal of a portion of the NPP as food., fibre and material of use for society (NPPh) and the change in land cover (∆NPPLu) (Haberl et al. 2007; Krausmann et al. 2013). Recent studies suggest that the indicators associated with HANPP provide key information for planning and evaluating ecosystem services (Mayer et al. 2021). HANPP considers NPP as the net amount of biomass produced by autotrophic organisms, in this case, plants, which constitute the primary energy source for the rest of the food chain for one year. In this sense, HANPP measures the:HANPPi=ΔNPPIn+NPPhΔNPPLu=NPP0-NPPact

In turn, ∆NPPLu is the difference between the potential NPP (Krausmann et al. 2013; available at https://www.aau.at/blog/global-hanpp-2000/) and the actual NPP (NPPact) based on disaggregated agricultural production data for the region obtained from the 3rd National Agricultural Census (2014) and based on (Guzmán-Casado et al. 2014).

To obtain the HANPP value per unit of analysis (hexagon), HANPP values (P) for each land cover i were multiplied by a wi coefficient representing the proportion of land cover i in each hexagon. HANPP units are presented in Tons of C / ha.HANPPhex=∑i=1kwiPi

Intermediate disturbance complexity (IDC)

The IDC model proposed by Marull et al. (2016a, b) transfers the concept of intermediate disturbance in natural ecosystems (Cornell 1978) to human-transformed landscapes (e.g., agroecosystems). The IDC argues that heterogeneous and well-connected land covers, with intermediate levels of agricultural activity, reflect an interaction between landscape complexity and energy availability that constitutes an agroecological matrix capable of harbouring great biodiversity (Loreau 2000; Tscharntke et al. 2012). Therefore, the IDC assesses the landscapes' capacity to host biodiversity and provide ecosystem services (Marull et al. 2019, 2018, 2016b).

The IDC model is calculated from the biomass available for other species (1-HANPP / 100) and the complexity of the landscape (Le). Le describes in a combined way the patterns (L) and processes (ECI) of the landscape (Marull et al. 2018).Le=aL+bECI101/a+b

where a and b are the canonical coefficients for the ortorthogonalization of the indices.IDC=Le1-HANPP/100

Ecosystem services

The UCRV landscapes' supply and demand of ecosystem services were assessed by Tabares-Mosquera et al. (2020) and based on the Common International Classification of Ecosystem Services (CICES) v.4.3 typology (Haines-young and Potschin 2012), which is associated with the defined ecosystem service categories: provisioning, regulation and cultural, of the Millennium Ecosystem Assessment (2005). Twenty-one ecosystem services were selected that appropriately fit the spatial scale of the Cali metropolitan phenomenon (Table 2).Table 2 Ecosystem services considered in the analysis

Section	Division	Code	Group	
Provisioning	Nutrition	1.1.1	Biomass	
1.1.2	Water	
Materials	1.2.1	Biomass	
1.2.2	Water	
1.2.3	Metallic and non-metallic abiotic materials	
Energy	1.3.1	Biomass	
1.3.2	Renewable abiotic source	
1.3.3	Non-renewable abiotic source	
Regulation and support	Regulation of waste, toxic substances, and other nuisances	2.1.1	Mediation by living systems	
Flow regulation	2.2.1	Mass flows	
2.2.2	Liquid flows	
2.2.3	Gas/air flows	
Maintenance of physical, chemical, and biological conditions	2.3.1	Maintenance of the life cycle, habitat, and protection of the gene pool	
2.3.2	Control of pests and diseases	
2.3.3	Soil formation and composition	
2.3.4	Maintenance of the chemical composition of water	
2.3.5	Atmospheric composition and climate regulation	
Cultural	Physical and intellectual interaction	3.1.1	Physical experience	
3.1.2	Intellectual and representative	
Symbolic and spiritual interaction	3.2.1	Spiritual or emblematic	
3.2.2	Existence and natural intrinsic value	
Source: Tabares-Mosquera et al. 2020, adapted from Haines-young and Potschin 2012

Since the land cover pattern is one of the most critical factors affecting the ability of a landscape to provide ecosystem services (Burkhard et al. 2009), the land cover map was used as the basis for quantifying multifunctionality and the capacity to provide ecosystem services of the landscapes that make up the UCRV. Ecosystem services were assessed using an expert-knowledge approach, given the region's incipient developments of ecosystem service mapping (Jacobs and Burkhard 2017). An interdisciplinary group of Twenty-seven experts was selected to evaluate the capacity to supply and demand ecosystem services for each land cover type (Tabares-Mosquera et al. 2020). The experts had to: i) be specialists in at least one of the following groups of land use: agricultural production areas, forests, seminatural areas, humid areas, or artificial areas; ii) be knowledgeable about the study area; and iii) be affiliated with public administrative institutions (mayors and governments), universities or private research centres.

The experts qualitatively evaluated the capacity of land covers to supply or demand ecosystem services, based on a six-class Likert scale: not relevant (0); very low (1); low (2); medium (3); high (4); very high (5) (Albert et al. 2016; Koschke et al. 2012). Experts' responses were averaged to obtain the supply and demand matrices. Finally, two criteria were established to evaluate the ecosystem function of each land cover. The first criterion, 'capacity', is defined as the long-term ability to provide different ecosystem services (Jacobs and Burkhard 2017). Capacity was calculated from the difference between ecosystem services supply and demand and can take values between -5 (very low capacity) and 5 (very high capacity) (Burkhard et al. 2014, 2009). The second criterion, 'multifunctionality', accounts for the number of ecosystem services of various categories (i.e., provisioning, regulation, and cultural) offered by each land cover class(Tabares-Mosquera et al. 2020). The Multifunctionality can take values from 0 (minimum) to 5 (maximum).

The Ecosystem Services Capacity (ESCj) and Multifunctionality of Ecosystem Services (MFj) for each hexagon (j) were calculated as the sum of each land cover capacity or multifunctionality within j, and weighted by the proportion of land cover I in j (Pij).ESCj=∑j=1kESCi∙PijMFj=∑j=1kESCi∙Pij

Statistical analyses

To assess the contribution of the different land covers to the expression of the IDC, a step-wise multiple regression model (MRM) was performed where IDC was the dependent variable, and the 22 land cover classes of the area were the predictor variables. Additionally, a linear regression analysis was performed to explore the relationship between IDC and the ESC based on the results obtained for each unit of analysis (hexagons; n = 502).

Results and discussion

Ecological functions and services of the metropolitan green infrastructure

Landscape metrics

The results indicate the highest levels of anthropogenic disturbance in the flat zone of the UCRV, where the main urban and industrial centres, road infrastructure, and sugarcane monoculture are located (HANPP > 61%; Fig. 2c). HANPP shows the existence of a gradient in the intensity of agricultural land use, which decreases with elevation as we move to the southern zone of the metropolitan region. This spatial pattern of anthropic disturbance would be associated with a disruption of ecological connectivity (ECI) in the flat zone of the study area (Fig. 2a). This area is also characterised by low landscape heterogeneity (H') (Fig. 2b). On the contrary, a strip of higher connectivity (Fig. 2a) and heterogeneity (Fig. 2b) stands out on the slopes of both mountain ranges, which also coincide with intermediate levels of anthropic disturbance (Fig. 2c).Fig. 2 Ecological functions and services of the metropolitan region of Cali. Ecological Connectivity Index -ECI (a); Shannon Index -H' (b); Human Appropriation of Net Primary Production -HANPP (c); Ecosystem Services Capacity -ESC (d); Multifunctionality of Ecosystem Services -MF (e), and Intermediate Disturbance Complexity -IDC (f)

The ECI levels reflect an overall ecological disconnection between the valley and key natural areas (i.e., forests, natural shrubs, natural grasslands; Fig. 1) mainly located at high elevations of both Andean mountains (Fig. 3). Even though the ECI model considers the hydrological network a fundamental element for ecological connectivity, the connective function exerted by rivers, especially the Cauca River, is very subtle. This situation contrasts with the one found in other metropolitan regions of the world, where riparian ecosystems play a key role in maintaining ecological connectivity even within highly antrhopogenized areas (Dupras et al. 2016; Padró et al. 2020a).Fig. 3 Ecological Connectivity Index (ECI) of the metropolitan region of Cali

The low connectivity associated with the hydrological network (Fig. 3) could be explained by the low ecological affinity between forests and pastureland covers. In this region, pastures play a critical and controversial role in the landscape functioning of this region and the socioecological sustainability. For instance, many clean pastures adjacent to rivers result from deforestation driven by extensive cattle ranching and land speculation (Garcia Corrales et al. 2019; Rodríguez Eraso et al. 2013; Zuluaga et al. 2021). The land cover changes associated with pastures establishment is a phenomenon risking homogenising the UCRV and overall mountain landscapes in the tropical Andes. Similarly, riparian natural land covers have also been pushed to the limits by sugarcane monoculture (Ayala-Osorio 2019; Delgadillo-Vargas 2014; Pérez-Rincón 2018). Additionally, although beyond the scope of this work, illegal human settlements and resulting water pollution are critical issues for the Cauca River (Holguin Gonzalez and Goethals 2010). All these conditions have limited the river's ecological connectivity potential. These findings are aligned with the water use conflict between agricultural demands and urban water consumption (Pérez-Rincón et al. 2011) and offer a comprehensive analysis of the larger river basin, complementing prior research conducted at local and sub-basin scales which identified desertification, high scarcity indexes (Pérez-Rincón et al. 2011), and river reversion (Delgadillo-Vargas 2014; Marull et al. 2017) in the region.

Landscape-metabolism model (IDC)

The IDC model serves as a valuable tool for exploring landscape sustainability by investigating the regional biophysical connections between humans and nature (Wu 2021). It does so by elucidating the landscape-metabolic relationships that reflect the complex socio-ecological processes occurring in this rapidly evolving metropolitan region. The relationship between land cover types and the Intermediate Disturbance Complexity (IDC) within the study area showed that the overall model fit statistics, including the multiple R-squared (0.7433) and adjusted R-squared (0.7381), indicating that approximately 74.33% of the variance in IDC values could be explained by the land cover variables included in the analysis. Similarly, the results revealed notable associations between specific land cover types and IDC values (Table 3). Specifically, sugarcane plantations, urban areas, and natural areas (dense forest, natural grasslands) have a significant (p < 0.050) negative relationship with the IDC, while Natural shrublands, cropland, and pasture mosaics, rivers and water bodies, and clean pastures significantly contribute to higher IDC values (p < 0.050) (Table 3).Table 3 Step-wise multiple regression model

Residuals:	Min	1Q	Median	3Q	Max	
-0.17760	-0.04293	-0.00521	0.02989	0.21350	
Variables/Coefficients	Estimate	Std. Error	t value	Pr( >|t|)	Significance

(< 0.050)

	
Intercept	0.43586	0.007971	54.683	 < 2e−16	*	
Sugarcane plantations	-0.26907	0.010776	-24.970	 < 2e−16	*	
Dense forest	-0.24005	0.014356	-16.721	 < 2e−16	*	
Urban areas	-0.40903	0.027881	-14.671	 < 2e−16	*	
Natural grasslands (Paramo)	-0.14953	0.016800	-8.901	 < 2e−16	*	
Natural shrublands	0.24755	0.031203	7.934	 < 1.45e−14	*	
Cropland and pasture mosaic	0.13680	0.031376	4.360	 < 1.58e−05	*	
Rivers and natural water bodies	0.90214	0.256902	3.512	0.000486	*	
Clean pastures	0.06582	0.027501	2.394	0.017062	*	
Weeded pastures	-0.07655	0.042348	-1.808	0.071252		
Cropland mosaic	0.25399	0.175123	1.450	0.147594		
Dependent variable: IDC. Predictive variables: 22 land covers. n: 502. Multiple R-squared: 0.7433. Adjusted R-squared: 0.7381. Residual standard error: 0.06372 on 491 degrees of freedom. F-statistic: 142.2 on 10 and 491 DF. p-value: < 2.2e−16

Significance codes: * < 0.050

These relationships reflect the distinct levels of ecological disturbances exerted by various agricultural practices of the region resulting in a land cover and land-use intensity spatial gradient (Fig. 2f). In this sense, the IDC revealed the existence of at least three types of landscape-metabolic configurations of the UCRV landscapes (Fig. 2f). The first type (anthropic) is defined by a very low IDC (0.05 < IDC < 0.3). This type depicts landscapes resulting from an industrial agriculture metabolism and is mainly concentrated along the river valley's flat area. The second type (natural areas) is defined by low-to-moderate IDC, low anthropogenic activity, low land cover heterogeneity, predominantly dense forests, and paramos. These areas are located mainly above the 3000 m a.s.l, where population density and activity are very low, usually associated with subsistence agriculture. A high IDC defines the third type of landscape-metabolic configuration (mosaics). It reflects heterogeneous landscapes with less intensive agricultural activities (intermediate disturbance levels). These landscapes are found at mid-elevations (1,200 m a.s.l. to 2,800 m a.s.l.) on the slopes of the Andean Mountain ranges. This area comprises agricultural, agroforestry, and agropastoral mosaics (Fig. 1), reflecting different agricultural practices, including traditional peasant, Afro, and Indigenous agroecological models (Duarte Torres et al. 2018). The findings build upon those of Marull et al. (2017) by extending the study beyond the administrative boundaries of a single municipality to encompass the Cauca River Valley region. This approach highlights the importance of adopting a green infrastructure approach within the metropolitan region of Cali and incorporating a broader spectrum of biocultural landscapes representative of the entire biogeographic region.

The significance of integrating findings from this landscape-metabolism approach into the strategic planning of sustainable metropolises, particularly within socioeconomically intricate regions like the tropical Andes, is underscored by its contextualization within the land-sharing and land-sparing discourses (Fischer et al. 2014, 2008; Grass et al. 2019). The proposed classification of landscape metabolic configurations resulting from the IDC analysis, identifies, on one side, anthropic landscape-metabolic configuration is characteristic of a productive paradigm that has often been associated with land-sparing strategies to balance biodiversity conservation and agricultural production (mainly for food but see Anderson-Teixeira et al. (2012)). Secondly, the natural landscape-metabolic configurations concentrated in protected areas host critical ecosystems and perform essential functions and services for natural and societal communities. However, these areas are highly disconnected, possibly affecting the effectiveness of a land-sparing strategy (Cannon et al. 2019; Edwards et al. 2021). Finally, the areas with landscape-metabolic configurations of mosaics would reflect a land-sharing strategy for food production and biodiversity conservation (Perfecto et al. 2009; Perfecto and Vandermeer 2008). These areas increasingly offer suitable habitats, higher biodiversity, and multiple ecosystem services (Loreau et al. 2003; Marull et al. 2019; Margalef 1973; Tscharntke et al. 2005), and play a crucial role in ecological connectivity (Fig. 3). To illustrate this, we explore the relationship between the IDC and the capacity of metropolitan landscapes to supply ecosystem services (ESC) (Fig. 4). Landscape-metabolic configurations related to agro-industrial activity are associated with a lower capacity to supply ecosystem services for the metropolitan population (yellow dots in Fig. 4). On the contrary, the agricultural mosaic revealed a higher capacity to supply ecosystem services such as water supply, food production, and flood regulation among others (orange dots in Fig. 4). Since the IDC is based on the theoretical assumption that agricultural landscapes can retain more farm-associated biodiversity at intermediate levels of human net primary production appropriation (Loreau et al. 2003; Marull et al. 2015; Montero et al. 2021), the IDC predictive power decreases as non-anthropogenic land covers (i.e., natural forest, shrublands, pasturelands, and paramos) increases, as seen in Fig. 3. Accordingly, metropolitan green infrastructure must contain different interrelated and connected elements to provide a structure that provides distinct functions and services (Basnou et al. 2020; Hansen and Pauleit 2014). The ESC of landscapes is expressed along the land-use intensity gradient, with the lower capacity index where urban and sugarcane land covers are predominant (Fig. 2d). This is expected, as these land covers have a high demand for ecosystem services while offering none (i.e., urban) or very few (i.e., sugarcane monocultures) ecosystem services, for instance, the capacity to provide energy in the form of biomass. In contrast, the strip of intermediately disturbed well-connected and heterogeneous landscapes (higher IDC values) shows higher capacity and multifunctionality to supply ecosystem services.Fig. 4 Relationship between the Intermediate Disturbance Complexity (IDC) model and the Ecosystem Services Capacity (ESC) in the metropolitan region of Cali. Note: The points represent each of the 502 landscape-scale analysis units. Anthropogenized areas: pastures, agricultural mosaics, sugarcane plantations, built-up areas. Natural areas: forests, grasslands, and shrubs. Colours indicate the predominant land cover in each unit of analysis (i.e., those occupying more than 50% of the total hexagon area): Natural areas (i.e., forests, grasslands, and shrublands), Agricultural mosaics, sugarcane plantations, and build-up areas (i.e., urban, suburban, and industrial)

We propose that the role of these biocultural agricultural mosaics is increasingly pertinent to the ongoing process of metropolitan growth in this region. While the nuances of the effectiveness of the land-sharing or land-sparing strategies go beyond the purpose and scope of this study and have already been widely revised in the literature (Fischer et al. 2017, 2014; Grau et al. 2013; Kremen 2015; Scariot 2013), their postulates can complement our analysis. The UCRV reflects the necessity of combined land-sharing and land-sparing strategies, the integrative assessment exposes that the effectiveness of these two strategies would be highly compromised by the poor ecological connectivity of the area. Our assessment reveals an important emerging category of metropolitan open spaces: the agrosilvopastoral mosaics. These biocultural landscapes hold a high capacity to supply ecosystem services and high ecological connectivity and help structure a well-connected, multifunctional green infrastructure (Table 4 and Table 5). Therefore, land use planning policies should consider the socioecological impacts of maintaining large agricultural areas with low levels of provision of ecosystem services (i.e., intensive sugarcane cultivation) at the cost of losing the ecological quality of the Cauca Valley by decreasing the provision of other essential ecosystem services and condemning the most populated area of the UCRV to ecological isolation and environmental degradation (Kremen and Miles 2012).Table 4 Ecosystem service capacity of the land covers of the metropolitan region of Cali

*These land cover classes were not mapped in Fig. 1 given their low representation in the study area

Table 5 Ecosystem services multifunctionality of the land covers of the metropolitan region of Cali

*These classes were not mapped given their low representation in the study area

Our results show that the ecological rupture affects essential ecosystem flows between the paramos and high Andean Forest. This situation jeopardises the delivery of vital ecosystem services from these areas to the metropolis, including collecting, regulating, and providing water for human consumption and agriculture in the entire region (Table 4), especially, in a tropical Andes region where climatic phenomena, such as "El Niño," exacerbated by climate change trends, can bring extreme drought conditions, amplifying the risks of water shortages, crop failures, and floodings. These challenges highlight the urgent need for adaptive strategies and sustainable water resource management practices to mitigate the socio-economic and environmental impacts of these events and other manifestations of climate change in the region. Therefore, a systemic green infrastructure approach can play a pivotal role in both mitigation and adaptation efforts, offering nature-based solutions to enhance water retention, reduce flood risks, and promote ecosystem resilience in the face of changing climatic conditions.

Our results underscore the significance of preserving a diverse metropolitan green infrastructure and comprehending it not merely as the sum of its parts, but as a cohesive functional unit—a complex system wherein each component interacts through landscape-metabolic processes, giving rise to emergent properties that support various functions and services crucial for the metropolitan region. These functions include connectivity, complementarity, and the provision of ecosystem services. Furthermore, the proposed systemic perspective interprets Cali's process of metropolization as the interaction between the built environment and green infrastructure, from which humans derive benefits. This perspective integrates aspects of green infrastructure as a subsystem into metropolitan urban analyses, recognizing its capacity to provide benefits to society (Tabares-Mosquera et al. 2020).

These results confirm the importance of polyculture landscapes and different types of forests (dense and riparian), which not only provide resources to metropolitan areas but also regulate ecosystem functions like supporting life cycles, protecting habitats and gene pools, pest control, and soil health (Bennett and Radford 2008; Grass et al. 2019; Kennedy et al. 2013; Tscharntke et al. 2012). These findings are pertinent to land scarcity discussions and can inform analyses of conflicts related to biofuel, food, and biodiversity production (Fischer et al. 2014; Lambin and Meyfroidt 2011). Thus, it's crucial to explore various agricultural practices (conventional, organic, agroecological) at both farm and landscape scales, as they can significantly influence agroecosystem functioning and contribute to landscape sustainability (Guzmán et al. 2011) and ecosystem services (Padró et al. 2020a).

Importance of the metropolitan green infrastructure for land use planning

In this study, metropolitan green infrastructure is defined as the collection of natural and semi-natural land covers whose condition, extent, and connectivity allow them to provide ecosystem services (Tabares-Mosquera et al. 2020). The results derived from the IDC model reveal a gradient of landscape-metabolic configurations within the metropolitan green infrastructure. These configurations, shaped by the intricate interactions among different green infrastructure elements, align with the IDC principles, where higher IDC values correlate with increased capacity and multifunctionality in providing ecosystem services (Marull et al. 2015). This gradient-based approach emphasizes the varied roles that landscape components play in supplying ecosystem services (Jones et al. 2013), shedding light on the complex relationship between these elements within the metropolitan area. For instance, the COVID-19 pandemic made evident the food sovereignty challenges of a metropolitan region (Rankin et al. 2021) isolated in landscape-metabolic configurations resulting from sugarcane industrialized monocultures, and it became evident that this homogeneous landscape-metabolic territorial set-up reduces the ecosystem service supply, and ultimately, the resilience of the entire metropolitan system.

Lastly, our research provides key insights into sustainable landscape patterns (Peng et al. 2021). The results support the idea that complex and heterogeneous landscape configurations can enhance ecosystem services while promoting landscape connectivity and multifunctionality, underscoring the intricate interactions among green infrastructure elements, landscape sustainability and resilience in metropolitan regions.

To offer valuable guidance for policy-makers and land planners aiming to develop strategies that enhance long-term sustainability at the landscape level, with a focus on promoting landscapes that are both ecologically robust and socially beneficial, we propose five reflections to guide land-use planning discussions in the Upper Cauca River Valley (UCRV) and other metropolitan regions in the tropical Andes, particularly those facing challenges from land-use change and intensification.

Adopting a landscape approach

Adopting the landscape scale for urban and territorial planning in metropolitan systems is increasingly relevant for land planning. In the first place, many of the ecosystem functions and services of the metropolises occur at this scale, which is key to overcoming the methodological difficulties of their evaluation, monitoring, and incorporation into public policy. The landscape scale allows closing metabolic cycles (water, energy, food, waste,) necessary to move towards a more circular economic model, which implies fostering interactions between built and open spaces, relationships of great relevance to society (Bennett and Radford 2008; Tello et al. 2017). Finally, considering the landscape as one of the essential background elements for the sustainability of the metropolitan system allows, under the right conditions, the reproduction of biophysical flows essential for its sustainability, such as those related to the agri-food system (Cattaneo et al. 2018; Marull et al. 2016a).

Re-establish the functional structure of the landscape

Restoring landscape complexity (i.e., heterogeneous, and well-connected land covers) would imply a substantial and long-term productive transformation of the valley. However, interstitial spaces in the flat plain can be a starting point for reconfiguring the territory and promoting an agricultural mosaic with a high capacity to provide ecosystem services. Improving the ecological structure of pastureland covers (i.e., promoting wooded pastures) would be crucial because of their affinity with agroforestry mosaics and heterogeneous crops, serving as stepping stones for ecological processes.

Enhancing ecological connectivity and protecting high mountain nature

Altitudinal gradients are associated with highly endemic biodiversity (Larsen et al. 2009). Therefore, a well-connected network of natural areas should facilitate the altitudinal migration of species threatened by global warming, anthropogenic habitat loss, and degradation (Balthazar et al. 2015; Cresso et al. 2020; Lambin et al. 2003). For instance, despite occupying less than 8% of the study area, the land covers associated with the Páramos ecosystem (i.e., grasslands and natural grasslands) have a high capacity to supply essential ecosystem services. A region with good connectivity at the macro-basin scale will likely counteract the impact suffered by the hydrological system of the UCRV.

Configure mosaic landscapes following agroecological management

The results show the great weight that the agrosilvopastoral mosaic area exerts on the provision of ecosystem services and the ecological connectivity of the metropolis. However, these values may vary depending on the type of agriculture practised (e.g., conventional, organic, agroecological) (Font et al. 2020; Marull et al., 2020; Padró et al. 2020b). Given this condition, the international consensus points to a necessary global agroecological transition, for which Latin America plays a fundamental role (Altieri and Nicholls 2012; Altieri and Toledo 2011; Jeanneret et al. 2021; Perfecto and Vandermeer 2008). The priorities should focus on improving the metropolitan areas' capacity to close metabolic cycles (Billen et al. 2021; Cattaneo et al. 2018) and providing multiple ecosystem services. This can be achieved in combination with promoting highly multifunctional land uses, where crop and livestock systems are integrated. The potential of agricultural mosaics for sustainable land use planning and management in metropolitan regions is critical (Tscharntke et al. 2021).

Define and implement a metropolitan green infrastructure

Based on the five previous elements, both the needs and the opportunities to adopt a conceptual and methodological green infrastructure framework to face the sustainability challenges of the UCRV are evident. The results are in line with cutting-edge studies on urban and territorial planning of metropolitan areas, which increasingly reinforce the need to include agricultural open spaces as fundamental elements for the sustainability of the metropolitan system given its multifunctional character (Basnou et al. 2020; Slätmo et al. 2019; Yacamán-Ochoa et al. 2020). This implies defining a network of interconnected open spaces, including peri-urban and rural spaces and natural spaces, capable of providing diverse ecological services, goods, and functions for society. The agrarian elements of this green infrastructure are fundamental, and, therefore, it is recommended to include them in the landscape and urban planning of metropolitan regions.

Conclusions

This article presents a territorial-metabolic approach for studying landscape sustainability and offers a framework that explores the interconnectedness of landscape configuration, ecosystem services, and, ultimately, human well-being. The approach integrates landscape-metabolism to assess how urban growth and agricultural intensification impact the socioecological functionality of a tropical Andean metropolitan region. Our results elaborate on Marull et al. 2017, by providing a classification of landscape metabolic configurations, and relating them with the provision of ecosystem services (i.e., water supply, food production, flood regulation among others), increasingly pertinent to the ongoing process of metropolitan growth in this region and the occurrence of regular yet increasingly severe climatic phenomena. Different landscape-metabolic configurations along a land-use gradient are linked to ecosystem service provisioning. The findings reveal that the current land use planning model, driven by agricultural intensification and industrialization, has degraded the metropolitan socioecological quality, undermining sustainable urban and rural progress. The proposed landscape-metabolism model aims to guide decisions in land use planning and territorial management to enhance sustainability and the multifunctionality of biocultural landscapes for ecosystem service provision.

The Intermediate Disturbance Complexity (IDC) model, originally validated in Mediterranean ecosystems, was applied in a tropical Andean region. Evaluating ESC in such intricate landscapes presents challenges, especially in economically developing nations constrained by resources. However, the IDC model offers a methodologically straightforward solution, utilizing a land cover map and data on metropolitan green infrastructure metabolism (e.g. agricultural production and management) to compute an indicator closely linked to ESC and multifunctionality. The IDC model distinguishes itself for its efficiency and efficacy in capturing not only landscape dynamics but also furnishing insights into landscape configuration and social metabolism, allowing researchers and land planners to understand the importance of agroecological landscapes for socioecological sustainability (Tello & González de Molina 2023).

Future research should prioritize assessing synergies and trade-offs in territorial planning to address multiple socio-ecological challenges, such as enhancing ecological connectivity in the valley's hydrological network and mitigating conflicts over water for human consumption and agriculture. Given the prevalence of inequality and poverty in tropical Andes metropolitan regions and Latin America, future socio-ecological studies should also incorporate environmental justice considerations into their assessments and objectives.

Authors' contributions

MJLRA: Conceptualization, Methodology, Formal analysis, Investigation, Data Curation, Writing—Original Draft, Writing—Review & Editing, Visualization. EZC: Conceptualization, Methodology, Investigation, Data Curation, Writing—Original Draft, Writing—Review & Editing. OB: Methodology, Investigation, Writing—Review & Editing, Supervision, Funding acquisition. JM: Conceptualization, Methodology, Writing—Review & Editing, Supervision.

Funding

LaRota-Aguilera received pre-doctoral funding from the’Sustainable Farm Systems: Long-Term Socioecological Metabolism in Western Agriculture’ (2012–2018) funded by the Social Sciences and Humanities Research Council of Canada (SSHRC 895–2011-1020) and the Javeriana University of Cali Colombia (2019). Zapata-Caldas and Buitrago-Bermudez received funding from Universidad del Valle (2020), Colombia (project 4339).

Data availability

The data used in this study is available upon request by contacting the authors directly.

Declarations

Ethical approval

Not applicable.

Competing interests

We declare that there are no conflicting interests associated with this research paper. We have no financial, personal, or professional affiliations that could be perceived as conflicting with the objectives, methods, or outcomes of this research. This includes any financial interests, employment, consultancies, patent applications, or any other potential conflicts of interest. We are committed to upholding the highest standards of transparency and integrity in the presentation of our research findings.

1 The population figures provided for the Upper Cauca River Valley in this study correspond to the aggregated municipal population reported by the 2018 National Census for the 25 municipalities considered part of the UCRV region in this study. Municipalities: Buenos Aires, Caldono, Cali, Caloto, Candelaria, Corinto, El Cerrito, Florida, Ginebra, Guacari, Guachene, Jambaló, Jamundí, Miranda, Morales, Padilla, Palmira, Pradera, Puerto Tejada, Santander de Quilichao, Suarez, Toribío, Vijes, Villa Rica, Yumbo,

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

Aide M Grau R Graesser J Andrade-Nuñez MJ Aráoz E Barros AP Campos-Cerqueira M Chacon-Moreno E Cuesta F Espinoza R Peralvo M Polk MH Rueda X Sanchez A Young KR Zarbá L Zimmerer KS Woody vegetation dynamics in the tropical and subtropical Andes from 2001 to 2014: Satellite image interpretation and expert validation Glob Change Biol 2019 25 2112 2126 10.1111/gcb.14618
Albert C Galler C Hermes J Neuendorf F von Haaren C Lovett A Applying ecosystem services indicators in landscape planning and management: The ES-in-Planning framework Ecol Ind 2016 61 100 113 10.1016/j.ecolind.2015.03.029
Altieri, M.A., Nicholls, C.I., 2012. Agroecology Scaling Up for Food Sovereignty and Resiliency 1–29. 10.1007/978-94-007-5449-2
Altieri MA Toledo VM The agroecological revolution in Latin America: Rescuing nature, ensuring food sovereignty, and empowering peasants J Peasant Stud 2011 38 587 612 10.1080/03066150.2011.582947
Anderson-Teixeira KJ Duval BD Long SP Delucia EH Biofuels on the landscape: Is “land sharing” preferable to “land sparing”? Ecol Appl 2012 22 2035 2048 10.1890/12-0711.1 23387108
Angotti T Latin American urbanization and planning: inequality and unsustainability in North and South Latin Am Perspect 1996 23 12 34 10.1177/0094582X9602300403
Ayala-Osorio G El monocultivo de la caña de azúcar en el valle geográfi co del río Cauca (Valle del Cauca, Colombia): un enclave que desnaturaliza la vida ecosistémica Forum. Revista Departamento de Ciencia Política 2019 15 37 66 10.15446/frdcp.n15.72452
Balatsky AV Balatsky GI Borysov SS Resource demand growth and sustainability due to increased world consumption Sustainability 2015 7 3430 3440 10.3390/su7033430
Balthazar V Vanacker V Molina A Lambin EF Impacts of forest cover change on ecosystem services in high Andean mountains Ecol Ind 2015 48 63 75 10.1016/j.ecolind.2014.07.043
Basnou C Baró F Langemeyer J Castell C Dalmases C Pino J Advancing the green infrastructure approach in the Province of Barcelona: integrating biodiversity, ecosystem functions, and services into landscape planning Urban Forest Urban Green 2020 55 126797 10.1016/j.ufug.2020.126797
Benedict MA Mcmahon ET Green Infrastructure : Smart Conservation for the 21 Century Renew Resourc J 2002 20 12 18
Bennett, A.F., Radford, J.Q., 2008. Emergent Properties of Land Mosaics: Implications for Land Management and Biodiversity Conservation, in: Managing and Designing Landscapes for Conservation: Moving from Perspectives to Principles. 199–214. 10.1002/9780470692400.ch18
Billen G Aguilera E Einarsson R Garnier J Gingrich S Grizzetti B Lassaletta L le Noë J Sanz-Cobena A Reshaping the European agro-food system and closing its nitrogen cycle: The potential of combining dietary change, agroecology, and circularity One Earth 2021 4 839 850 10.1016/j.oneear.2021.05.008
Brondizio, E.S., Settele, J., Díaz, S., Ngo, H.T., 2019. Global assessment report on biodiversity and ecosystem services of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services. Bonn, Germany. 10.5281/zenodo.3831673
Burkhard B Kroll F Müller F Windhorst W Landscapes’ capacities to provide ecosystem services - A concept for land-cover based assessments Landscape Online 2009 15 1 22 10.3097/LO.200915
Burkhard B Kandziora M Hou Y Müller F Ecosystem service potentials, flows and demands-concepts for spatial localisation, indication and quantification Landscape Online 2014 34 1 32 10.3097/LO.201434
Canales AI Canales Cerón M De la metropolización a las agrópolis: El nuevo poblamiento urbano en el Chile actual Polis 2013 12 31 56 10.4067/s0718-65682013000100003
Cannon PG Gilroy JJ Tobias JA Anderson A Haugaasen T Edwards DP Land-sparing agriculture sustains higher levels of avian functional diversity than land sharing Glob Change Biol 2019 25 1576 1590 10.1111/gcb.14601
Cattaneo C Marull J Tello E Landscape Agroecology. The dysfunctionalities of industrial agriculture and the loss of the circular bioeconomy in the Barcelona Region Sustainability 2018 10 12 4722 10.3390/su10124722
Chatzimentor A Apostolopoulou E Mazaris AD A review of green infrastructure research in Europe: Challenges and opportunities Landsc Urban Plan 2020 198 103775 10.1016/j.landurbplan.2020.103775
Centro Nacional de Memoria Histórica (2015) Una nación desplazada: informe nacional del desplazamiento forzado en Colombia, Bogotá, CNMH - UARIV
Cornell JH Diversity in tropical rain forest and coral reefs Science 1978 199 4335 1302 1310 10.1126/science.199.4335.1302 17840770
Cresso M Clerici N Sanchez A Jaramillo F Future climate change renders unsuitable conditions for paramo ecosystems in Colombia Sustainability 2020 12 1 13 10.3390/su12208373 35136666
DANE (2018) Censo Nacional de Población y Vivienda 2018: Población censal ajustada por cobertura y porcentajes de omisión municipal por área. Retrieved from http://www.dane.gov.co/files/censo2018/poblacion_ajustada.pdf
Delgadillo-Vargas O (2014) La caña de azúcar en la historia ambiental del Valle geográfico del río Cauca (Doctoral dissertation). Pontificia Universidad Javeriana. Bogotá, Colombia. Retrieved from https://repository.javeriana.edu.co/handle/10554/15735
Demuzere M Orru K Heidrich O Olazabal E Geneletti D Orru H Bhave AG Mittal N Feliú E Faehnle M Mitigating and adapting to climate change: Multi-functional and multi-scale assessment of green urban infrastructure J Environ Manag 2014 146 107 115 10.1016/j.jenvman.2014.07.025
Duarte Torres CA LaRota-Aguilera MJ Gómez MC Ducká E Galindo PA Marulanda JL Valencia D Trujillo D La estructura de la propiedad rural en el Cauca : perspectivas sobre las necesidades de tierra en contextos interculturales 2018 Cali Universidad Javeriana de Cali
Dupras J Marull J Parcerisas L Coll F Gonzalez A Girard M Tello E González A Girard M Tello E Gonzalez A Girard M Tello E González A Girard M Tello E Gonzalez A Girard M Tello E The impacts of urban sprawl on ecological connectivity in the Montreal Metropolitan Region Environ Sci Policy 2016 58 61 73 10.1016/j.envsci.2016.01.005
Edwards FA Massam MR Cosset CCP Cannon PG Haugaasen T Gilroy JJ Edwards DP Sparing land for secondary forest regeneration protects more tropical biodiversity than land sharing in cattle farming landscapes Curr Biol 2021 31 1284 1293.e4 10.1016/j.cub.2020.12.030 33482111
Elmqvist T Fragkias M Goodness J Güneralp B Marcotullio PJ McDonald RI Parnell S Schewenius M Sendstad M Seto KC Wilkinson C Urbanization, Biodiversity and Ecosystem Services: Challenges and Opportunities, Urbanization, Biodiversity and Ecosystem Services: Challenges and Opportunities: A Global Assessment 2013 10.1007/978-94-007-7088-1_23
European Commission (2013) Green Infrastructure (GI) — Enhancing Europe’s natural capital (Communication from the Commission to the European Parliament, the Council, the European Economic and Social Committee and the Committee of the Regions). Retrieved from https://www.eea.europa.eu/policy-documents/green-infrastructure-gi-2014-enhancing
Fischer J Brosi B Daily GC Ehrlich PR Goldman R Goldstein J Lindenmayer DB Manning AD Mooney HA Pejchar L Ranganathan J Tallis H Should agricultural policies encourage land sparing or wildlife-friendly farming? Front Ecol Environ 2008 6 380 385 10.1890/070019
Fischer J Abson DJ Butsic V Chappell MJ Ekroos J Hanspach J Kuemmerle T Smith HG von Wehrden H Land sparing versus land sharing: Moving forward Conserv Lett 2014 7 149 157 10.1111/conl.12084
Fischer J Abson DJ Bergsten A French Collier N Dorresteijn I Hanspach J Hylander K Schultner J Senbeta F Reframing the Food-Biodiversity Challenge Trends Ecol Evol 2017 32 335 345 10.1016/j.tree.2017.02.009 28284373
Font C Padró R Cattaneo C Marull J Tello E Alabert A Farré M How farmers shape cultural landscapes. Dealing with information in farm systems (Vallès County, Catalonia, 1860) Ecological Indicators 2020 112 106104 10.1016/j.ecolind.2020.106104
Garcia Corrales LM Avila Rangel H Gutierrez Llantoy R Land-use and socioeconomic changes related to armed conflicts: A Colombian regional case study Environ Sci Policy 2019 97 116 124 10.1016/j.envsci.2019.04.012
Grass I Loos J Baensch S Batáry P Librán-Embid F Ficiciyan A Klaus F Riechers M Rosa J Tiede J Udy K Westphal C Wurz A Tscharntke T Land-sharing/-sparing connectivity landscapes for ecosystem services and biodiversity conservation People and Nature 2019 1 262 272 10.1002/pan3.21
Grau R Kuemmerle T Macchi L Beyond “land sparing versus land sharing”: Environmental heterogeneity, globalization and the balance between agricultural production and nature conservation Curr Opin Environ Sustain 2013 5 477 483 10.1016/j.cosust.2013.06.001
Güneralp B McDonald RI Fragkias M Goodness J Marcotullio PJ Seto KC Urbanization Forecasts, Effects on Land Use, Biodiversity, and Ecosystem Services Urbanization 2013 Challenges and Opportunities Biodiversity and Ecosystem Services 31 47
Guzmán GI de Molina MG Alonso AM The land cost of agrarian sustainability An Assess Land Use Policy 2011 28 4 825 835 10.1016/j.landusepol.2011.01.010
Guzmán-Casado GI, Aguilera E, Soto D, Cid A, Ruiz RG, Herrera A, Villa I, Molina MG de (2014) Methodology and conversion factors to estimate the net primary productivity of historical and contemporary agroecosystems. Working Papers Sociedad Española de Historia Agraria. Retrieved from https://repositori.uji.es/xmlui/handle/10234/91670
Haberl H Erb KH Krausmann F Gaube V Bondeau A Plutzar C Gingrich S Lucht W Fischer-kowalski M Quantifying and mapping the human appropriation of net primary production in earth ’ s terrestrial ecosystems Proc Natl Acad Sci USA 2007 104 12942 12947 10.1073/pnas.0704243104 17616580
Haines-Young R, Potschin M (2012) Common international classification of ecosystem services (CICES, Version 4.1). Eur Environ Agency 33:107
Hansen R Pauleit S From multifunctionality to multiple ecosystem services? A conceptual framework for multifunctionality in green infrastructure planning for Urban Areas Ambio 2014 43 516 529 10.1007/s13280-014-0510-2 24740622
Hibbs BJ Sharp JM Jr Hydrogeological impacts of urbanization Environ Eng Geosci 2012 18 3 24 10.2113/gseegeosci.18.1.3
Holguin GJ, Goethals Peter LM (2010) Modelling the ecological impact of discharged urban waters upon receiving aquatic ecosystems. A tropical lowland river case study: city Cali and the Cauca river in Colombia. Int Congr Environ Model Softw 539. https://scholarsarchive.byu.edu/iemssconference/2010/all/539
IPCC (2014) Cambio climático 2014: Informe de síntesis. Contribución de los Grupos de trabajo I, II y III al Quinto Informe de Evaluación del Grupo Intergubernamental de Expertos sobre el Cambio Climático, Contribución de los Grupos de trabajo I,II y III al Quinto Informe de Evaluación del Grupo Intergubernamental de Expertos sobre el Cambio Climático.
Jacobs S, Burkhard B (2017) Applying expert knowledge for ecosystem services quantification. In B. Burkhard & J. Maes (Eds.), Mapping Ecosystem Services (pp. 144-148). Pensoft Publishers. 10.3897/ab.e12837
Jeanneret Ph Aviron S Alignier A Lavigne C Helfenstein J Herzog F Kay S Petit S Agroecology landscapes Landscape Ecol 2021 36 2235 2257 10.1007/s10980-021-01248-0
Jones KB, Zurlini G, Kienast F, Petrosillo I, Edwards T, Wade TG,...& Zaccarelli N (2013) Informing landscape planning and design for sustaining ecosystem services from existing spatial patterns and knowledge. Landscape ecology, 28, 1175–1192
Kennedy CM Lonsdorf E Neel MC Williams NM Ricketts T Winfree R Bommarco R Brittain C Burley AL Cariveau D A global quantitative synthesis of local and landscape effects on wild bee pollinators in agroecosystems Ecology letters 2013 16 584 599 10.1111/ele.12082 23489285
Koschke L Fürst C Frank S Makeschin F A multi-criteria approach for an integrated land-cover-based assessment of ecosystem services provision to support landscape planning Ecol Indic 2012 21 54 66 10.1016/j.ecolind.2011.12.010
Krausmann F Schandl H Sieferle RP Socio-ecological regime transitions in Austria and the United Kingdom Ecol Econ 2008 65 187 201 10.1016/j.ecolecon.2007.06.009
Krausmann F Gingrich S Eisenmenger N Erb KH Haberl H Fischer-Kowalski M Growth in global materials use, GDP and population during the 20th century Ecol Econ 2009 68 2696 2705 10.1016/j.ecolecon.2009.05.007
Krausmann F Erb K-H Gingrich S Haberl H Bondeau A Gaube V Lauk C Plutzar C Searchinger TD Global human appropriation of net primary production doubled in the 20th century Proc Natl Acad Sci 2013 110 10324 10329 10.1073/pnas.1211349110 23733940
Kremen C Reframing the land-sparing/land-sharing debate for biodiversity conservation Ann N Y Acad Sci 2015 1355 52 76 10.1111/nyas.12845 26213864
Kremen, C., Miles, A., 2012. Ecosystem services in biologically diversified versus conventional farming systems: Benefits, externalities, and trade-offs. Ecology and Society 17. 10.5751/ES-05035-170440
Lambin EF Meyfroidt P Global land use change, economic globalization, and the looming land scarcity Proc Natl Acad Sci USA 2011 108 3465 3472 10.1073/pnas.1100480108 21321211
Lambin EF Geist HJ Lepers E Dynamics of land-use and land-cover change in tropical regions Annu Rev Environ Resour 2003 28 205 241 10.1146/annurev.energy.28.050302.105459
Larsen, T.H., Brehm, G., Navarrete, H., Franco, P., Gómez, H., Mena, J.L., Morales, V., Argollo, J., Blacutt, L., Canhos, V., 2009. Desplazamientos de los Rangos de Distribución y Extinciones Impulsados por el Cambio Climático en los Andes Tropicales : Síntesis y Orientaciones. Cambio Climàtico en los Andes Tropicales 57–82.
Loreau M Biodiversity and ecosystem functioning: recent theoretical advances Oikos 2000 91 3 17 10.1034/j.1600-0706.2000.910101.x
Loreau M Mouquet N Gonzalez A Biodiversity as spatial insurance in heterogeneous landscapes Proc Natl Acad Sci USA 2003 100 12765 12770 10.1073/pnas.2235465100 14569008
Margalef R Ecological theory and prediction in the study of the interaction between man and the rest of the biosphere 1973 Tecnologia i Cultura Medi Ambient 38
Martínez-Toro P Patiño-Gómez ZL Área metropolitana del sur del Valle y norte del Cauca (AMVACA). Definición, caracterización y propuesta de conformación 2015 Primera Programa Editorial, Santiago de Cali Universidad del Valle
Marull J Mallarach JM A GIS methodology for assessing ecological connectivity: application to the Barcelona Metropolitan Area Landsc Urban Plan 2005 71 243 262 10.1016/j.landurbplan.2004.03.007
Marull J Tello E Fullana N Murray I Jover G Font C Coll F Domene E Leoni V Decolli T Long-term bio-cultural heritage: exploring the intermediate disturbance hypothesis in agro-ecological landscapes (Mallorca, c. 1850–2012) Biodivers Conserv 2015 24 3217 3251 10.1007/s10531-015-0955-z
Marull J Font C Padró R Tello E Panazzolo A Energy-Landscape Integrated Analysis: A proposal for measuring complexity in internal agroecosystem processes (Barcelona Metropolitan Region, 1860–2000) Ecol Ind 2016 66 30 46 10.1016/j.ecolind.2016.01.015
Marull J Font C Tello E Fullana N Domene E Pons M Galán E Towards an energy–landscape integrated analysis? Exploring the links between socio-metabolic disturbance and landscape ecology performance (Mallorca, Spain, 1956–2011) Landscape Ecol 2016 31 317 336 10.1007/s10980-015-0245-x
Marull J Delgadillo-Vargas O Cattaneo C LaRota-Aguilera MJ Krausmann F Socioecological transition in the Cauca river valley, Colombia (1943–2010): towards an energy – landscape integrated analysis Reg Environ Change 2017 10.1007/s10113-017-1128-2
Marull J Tello E Bagaria G Font X Cattaneo C Pino J Exploring the links between social metabolism and biodiversity distribution across landscape gradients: A regional-scale contribution to the land-sharing versus land-sparing debate Sci Total Environ 2018 619–620 1272 1285 10.1016/j.scitotenv.2017.11.196 29734605
Marull J Herrando S Brotons L Melero Y Pino J Cattaneo C Pons M Llobet J Tello E Building on Margalef: Testing the links between landscape structure, energy and information flows driven by farming and biodiversity Sci Total Environ 2019 674 603 614 10.1016/j.scitotenv.2019.04.129 31026792
Mayer A Kaufmann L Kalt G Matej S Theurl MC Morais TG Leip A Erb K-H Applying the Human Appropriation of Net Primary Production framework to map provisioning ecosystem services and their relation to ecosystem functioning across the European Union Ecosyst Serv 2021 51 101344 10.1016/j.ecoser.2021.101344 34631401
Mcdonald R, Marcotullio PJ, Güneralp B (2013) Urbanization and global trends in biodiversity and ecosystem services. In: Elmqvist T, Fragkias M, Goodness J, Güneralp B, Marcotullio PJ, McDonald RI et al (eds) Urbanization, biodiversity and ecosystem services: Challenges and opportunities. Springer,  pp 31–52. 10.1007/978-94-007-7088-1_3
Montero A Marull J Tello E Cattaneo C Coll F Pons M Infante-Amate J Urrego-Mesa A Fernández-Landa A Vargas M The impacts of agricultural and urban land-use changes on plant and bird biodiversity in Costa Rica (1986–2014) Reg Environ Change 2021 21 1 19 10.1007/s10113-021-01767-1 33362432
Padró R la Rota-Aguilera MJ Giocoli A Cirera J Coll F Pons M Pino J Pili S Serrano T Villalba G Marull J Assessing the sustainability of contrasting land use scenarios through the Socioecological Integrated Analysis (SIA) of the metropolitan green infrastructure in Barcelona Landsc Urban Plan 2020 203 103905 10.1016/j.landurbplan.2020.103905 32834265
Padró R Tello E Marco I Olarieta JR Grasa MM Font C Modelling the scaling up of sustainable farming into Agroecology Territories: Potentials and bottlenecks at the landscape level in a Mediterranean case study J Clean Prod 2020 275 124043 10.1016/j.jclepro.2020.124043
Pauleit S, Vasquéz A, Maruthaveeran, S., Liu, L., Cilliers, S.S., 2021. Urban green infrastructure in the Global South. Urban ecology in the Global South. Springer, Cham 107–143.
Peng J Liu Y Corstanje R Meersmans J Promoting sustainable landscape pattern for landscape sustainability Landscape Ecol 2021 36 7 1839 1844 10.1007/s10980-021-01271-1
Pérez-Rincón MA Peña MR Alvarez P Agro-industria cañera y uso del agua: Análisis crítico en el contexto de la política de agrocombustibles en colombia Ambiente e Sociedade 2011 14 153 178 10.1590/S1414-753X2011000200011
Perfecto I Vandermeer J Biodiversity conservation in tropical agroecosystems: A new conservation paradigm Ann N Y Acad Sci 2008 1134 173 200 10.1196/annals.1439.011 18566094
Perfecto I Vandermeer J Wright A Nature’s matrix: linking agriculture, conservation and food sovereignty 2009 Routledge
Pino J Marull J Ecological networks: Are they enough for connectivity conservation? A case study in the Barcelona Metropolitan Region (NE Spain) Land Use Policy 2012 29 684 690 10.1016/j.landusepol.2011.11.004
Rankin S Hurtado-Bermudez LJ Bonilla Findji O Mosquera EE Lundy M Perfil del sistema alimentario de Cali, ciudad-región: Construyendo un entendimiento común que articule e impulse la acción 2021 Cali (Colombia) Alliance of Bioversity International and CIAT 27
Centro Nacional de Memoria Histórica (2014) “Patrones” y campesinos: tierra, poder y violencia en el Valle del Cauca (1960 – 2012). CNMH. Bogotá, Colombia. Retrived from: https://www.centrodememoriahistorica.gov.co/descargas/informes2014/patronesyCampesinos/patrones-y-campesinos-tierra-poder-violencia-valle-del-cauca-insertos-baja.pdf
Pérez-Rincón M (2018) Impactos y pasivos ambientales de la caña de azúcar en el Valle del río Cauca en Colombia. Universidad del Valle. Cali, Colombia
Rodríguez Eraso N Armenteras-Pascual D Alumbreros JR Land use and land cover change in the Colombian Andes: Dynamics and future scenarios J Land Use Sci 2013 8 154 174 10.1080/1747423X.2011.650228
Scariot A Land sparing or land sharing: The missing link Front Ecol Environ 2013 11 177 178 10.1890/13.WB.008
Shannon CE Weaver W A mathematical model of communication Bell Syst Tech J 1948 27 379 423 10.1002/j.1538-7305.1948.tb01338.x
Slätmo E, Nilsson K, Turunen E (2019) Implementing green infrastructure in spatial planning in Europe. Land 8. 10.3390/land8040062
Steel C (2008) Hungry city: How food shapes our lives. Vintage Books, London
Swyngedouw E Heynen NC Urban political ecology, justice and the politics of scale Antipode 2003 35 898 918 10.1111/j.1467-8330.2003.00364.x
Tabares-Mosquera M Zapata-Caldas E Buitrago-Bermúdez O Valuation of ecosystem services for the identification of metropolitan ecological structures: The case of Cali, Colombia Geograph Res Lett 2020 46 603 631 10.18172/cig.3952
Tello E González de Molina M Agrarian Metabolism and Socio-ecological Transitions to Agroecology Landscapes The Barcelona School of Ecological Economics and Political Ecology: A Companion in Honour of Joan Martinez-Alier 2023 Cham Springer International Publishing 93 107
Tello E González de Molina M de Molina MG Fraňková E Haas W Singh SJ Methodological Challenges and General Criteria for Assessing and Designing Local Sustainable Agri-Food Systems: A Socio-Ecological Approach at Landscape Level Socio-Metabolic Perspectives on the Sustainability of Local Food Systems: Insights for Science 2017 Policy and Practice Springer 27 68
Tscharntke T Tylianakis JM Rand TA Didham RK Fahrig L Batáry P Bengtsson J Clough Y Crist TO Dormann CF Ewers RM Fründ J Holt RD Holzschuh A Klein AM Kleijn D Kremen C Landis DA Laurance WF Lindenmayer DL Scherber C Sodhi N Steffan-Dewenter I Thies C van der Putten WH Westphal C Landscape moderation of biodiversity patterns and processes - eight hypotheses Biol Rev 2012 87 661 685 10.1111/j.1469-185X.2011.00216.x 22272640
Tscharntke T, Klein AM, Kruess I, Steffan-Dewenter I, Thies C (2005) Landscape perspectives on agricultural intensification and biodiversity: ecosystem service management. Ecol Lett 857–874. 10.1111/j.1461-0248.2005.00782.x
Tscharntke T, Grass I, Wanger TC, Westphal C, Batáry P (2021) Beyond organic farming – harnessing biodiversity-friendly landscapes. Trends in Ecology & Evolution 1–12. 10.1016/j.tree.2021.06.010
United Nations World population prospects 2019 2019 World Population Prospects Department of Economic and Social Affairs 2019
Uribe H (2017) Capítulo 20: El Valle geográfico del río Cauca: Un espacio transformado por el capital agroindustrial. Repositorio Institucional UNAD. Retrived from: https://repository.unad.edu.co/handle/10596/24849
Urrea-Giraldo F, Candelo-Álvarez AF (2017) Cali, ciudad región ampliada: una aproximación desde la dimensión étnica-racial y los flujos poblacionales. Sociedad y Economía 145–174. 10.25100/sye.v0i33.5628
Vásquez A, Giannotti E, Galdámez E, Velásquez P, Devoto C (2019) Green Infrastructure Planning to Tackle Climate Change in Latin American Cities BT - Urban Climates in Latin America, in: Henríquez, C., Romero, H. (Eds.), . Springer International Publishing, Cham, 329–354. 10.1007/978-3-319-97013-4_13
Vásquez É (2018) Historia de Cali en el siglo 20 (Sociedad, economía, cultura y espacio), Historia y Espacio. 10.25100/hye.v0i18.7017
Vermeulen SJ Campbell BM Ingram JSI Climate change and food systems Annu Rev Environ Resour 2012 37 195 222 10.1146/annurev-environ-020411-130608
Wu J Landscape ecology, cross-disciplinarity, and sustainability science Landscape Ecol 2006 21 1 4 10.1007/s10980-006-7195-2
Wu J Landscape sustainability science: ecosystem services and human well-being in changing landscapes Landsc Ecol 2013 28 6 999 1023 10.1007/s10980-013-9894-9
Wu J Landscape sustainability science (II): core questions and key approaches Landsc Ecol 2021 36 2453 2485 10.1007/s10980-021-01245-3
Yacamán-Ochoa C, Ferrer-Jiménez D, Mata-Olmo R (2020) Green Infrastructure Planning in Metropolitan Landscapes and Food Security. Land. 10.3390/land9110414
Zuluaga A, Etter A, Nepstad D, Chará J, Stickler C, Warren M (2021) Colombia’s pathway to a more sustainable cattle sector: A spatial multi-criteria analysis. Land Use Policy 109. 10.1016/j.landusepol.2021.105596
