1. Introduction
Urbanization is a complex process that involves various social, economic, and environmental processes that are continuously interacting and shaping the landscape. Economic development can be viewed, under a lens of increasing urbanization, as one metric of progress [
1,
2]. Industrialization has fostered the transition and expansion of economic activities, replacing primary economic sectors (e.g., agriculture and forestry) with modern industries that provide services and create new technologies that advance the human condition [
3,
4]. Consequently, urban areas have become the main human habitat. Urbanization is often problematic, since it is accompanied by rapid landscape transformations and socio-economic and demographic pressures that challenge the availability, access, and distribution of food. Additionally, rapid urbanization has displaced agricultural and forest land, thus impacting soils and food production. Local populations living in urbanized areas face multiple stressors such as environmental degradation, overcrowding, and lack of access to locally grown healthy food. Meeting demands for food with a region’s limited resources is challenging and requires other strategies to satisfy an increasing population’s needs, which may include new agricultural production technologies, importing foods, and changes in traditional practices. As agricultural land decreases due to land cover change, more production, or higher quality foods, are needed. High-intensity agricultural systems are required to increase crop production per unit area. However, intensive irrigation, the use of fertilizers and pesticides, and high-yielding crop varieties (HYV) have resulted in increased water pollution, soil erosion, and have decreased the fertility of soil [
5,
6]. Alternatively, crop-diversification production systems that can maximize specific nutrients but reduce high-intensity agricultural options, could be utilized. In this regard, environmental policies and strategies are crucial to conserve resources and ecosystem services. Additionally, sustainable land use planning and decision making are critical to balance the supply and demand of the resources and services needed to ensure food security for the population. While various approaches have been developed to foster sustainable land management and agricultural practices, these are often context-specific and mainly focus on the prevention of environmental degradation [
7,
8,
9]. There is a need for an approach that considers environmental sustainability but focuses on agricultural production for local nutrition.
Diverse conceptual frameworks have been developed to link the urban system, the agricultural system, the natural system, and human wellbeing, including nutritional aspects. Examples of these frameworks are the social–ecological system, the urban social-ecological–technological systems, the FEW nexus framework, the landscape sustainability science framework, and the urban ecological economic systems framework [
10,
11,
12,
13,
14,
15]. In this paper, we propose an urban nutrition (UN) extension (
Figure 1) to the urban ecological economic system framework, by developing a focused component that simulates how existing agricultural systems can be modified to achieve nutritional food security through the use of sustainable practices that yield increased crop production and better nutritional support. The urban ecological economic systems framework proposed by Huang and others in 1988 attempts to study urban societies using a holistic view [
15]. We use the UN extension to narrow the focus of the four main components of the framework: natural system, agricultural system, urban system, and life support service. The natural system is represented by ecosystem services, the agricultural system is represented by agricultural practices (monoculture versus crop diversity), the urban system is represented by land use change, and the life support service is represented by human nutrition. This research examines these complex relationships to help decision makers with policy decisions concerning land use choices that might impact the interrelationships between urbanization and food production. Details of this framework are described in the Methods section. Here, we use three hypothetical scenarios to evaluate how land use choices affect crop production, leading to impacts on human nutrition and the environment measured through soil erosion. The scenarios include three simulated landscapes with different levels of urban, agriculture, and forest-land use. Each scenario was analyzed for (1) soil erosion and (2) crop production under both high-intensity agricultural methods (monoculture cultivation) and diverse crop cultivation, to explore the minimum population size that can be supported given a fixed amount of agricultural land and identify the limited nutrients that could be supported under each strategy. This framework promotes local sustainability practices, such as the “farm-to-table” (or farm-to-fork) movement, by enabling decision makers to evaluate options under various land use constraints. This is important because as cities grow, populations shift to urban landscapes and disengage from healthy locally grown food.
Urbanization Impacts on Food and Health
Food production, influenced by urbanization and technology, has direct impacts on human wellbeing. Production systems such as crop diversification—the introduction of new crop species to diversify crop production—can be used to increase food security in communities, increase income on small farm holdings, reduce environmental degradation such as soil erosion, and mitigate climate change [
16]. This paper compares crop production systems that are driven by technologies and urbanization, such as monoculture production systems, where only one genetically identical crop is continuously grown, and crop diversification systems, where multiple crops are cultivated based on their nutrient content.
Monocultures are widely used in industrial agricultural systems since they allow the overproduction of a particular commodity that is valuable for society (i.e., corn, rice, wheat). For example, maize is used for several purposes, including fuel ethanol, industrial products, and livestock feed. Global dietary trends of increasing demand for foods of animal origin have resulted in a high percentage of agricultural land being allocated for growing feed for livestock rather than crops for direct human consumption. Growing a diversity of crops is a straightforward strategy for providing sufficient nutrients for a healthy urban population. A perspective of sustainability is needed to produce enough food for the current population, while conserving the land for future generations and ensuring human welfare. This perspective can center on the study and management of ecosystem services that directly or indirectly affect the supply of crops, including supportive (e.g., soil fertility), regulating (e.g., erosion control), and provisioning (e.g., food production) ecosystem services. Bommarco and others illustrate how ecological intensification through the management of ecosystem services can sustainably enhance food security globally [
17].
Linking sustainability with ecosystem services is complex, since the definitions of both are still being developed [
18,
19,
20]. For instance, some authors use “environmental sustainability” to prioritize conservation and to highlight the importance of maintaining natural capital [
21]. Other definitions attempt to capture the relationship between economic development and global ecology [
22]. Here, we use one of the most common definitions of sustainability: “meeting the needs of the present without compromising the ability of future generations to meet their own needs”, because it addresses the future of human societies and their environments [
23].
The overarching hypothesis of the work presented here is that, given a limited amount of cropland, agricultural strategies that diversify crop types nutritionally support larger populations, promote a healthier diet, and synergize ecosystem services, leading to more sustainable urban societies as compared to strategies that prioritize monocultures. This hypothesis was explored using synthetic landscapes where agricultural area and crop production systems were specified, and nutrient production and soil erosion were modeled for each landscape.
2. Data and Methods
2.1. Urban Nutrition Extension
The UN extension illustrates how the urban system, the natural system, and the agricultural system interact and affect nutritional security in local populations. The extension could help guide decision makers and stakeholders to differentially prioritize between land uses (urban surface versus cropland) and agricultural-production techniques (monoculture versus crop diversity), which affect ecosystem services (crop production and soil erosion) that promote urban population’s nutritional needs. Decisions are affected by many drivers, including the food system, agricultural technologies, science, education, economic growth, governance, and policies.
2.2. Coupling UN Extension with InVEST
The Integrated Valuation of Ecosystem Services and Tradeoffs (InVEST) software, version 3.9.0, was used to explore this interplay. The following steps describe the construction and assessment of scenarios of land use change, crop production, soil erosion, and nutrition: (1) development of hypothetical scenarios, (2) data collection and preparation, (3) modeling and analysis of the effects of each scenario on ecosystem services, and (4) comparison between nutrition and calories for each scenario.
This paper uses two InVEST models to evaluate the impacts on ecosystem services: the Crop Production-Percentile (CPP) model, and the Sediment Delivery Ratio (SDR) model. The InVEST software was developed by the Natural Capital Project in collaboration with Stanford University and several other academic and non-profit institutions to assist and facilitate the inclusion of ecosystem services in decision making [
24]. The different models can be applied in developed and developing countries and for different types of ecosystem services (e.g., provisioning, regulating, cultural) or multiple ecosystem services under different scenarios [
25]. The suite of models can be operated by non-experts, making it a relatively simple and valuable decision support tool. However, the complexity of the model can vary, and availability of data can be problematic [
26].
The CPP model evaluates different strategies and scenarios of crop production impacted by land use/land cover change across the globe depending on climate. The model uses the United Nations Food and Agriculture Organization (FAO) database comprised of 175 kinds of crops. This database is supplemented by national and regional datasets [
27,
28]. The SDR model uses global data to determine annual soil loss in a determined area. This model is based on the Universal Soil Loss Equation and requires spatially explicit data including a land use/land cover raster, soil type, topography, and rainfall intensity data [
26,
28,
29].
2.3. Development of Hypothetical Scenarios
Three land use scenarios were developed to assess crop production, soil erosion, and nutritional yield. Each scenario utilizes a synthetic landscape (after Smith and Atkinson [
30]) that emphasizes one of three planning preferences: (1) economic and urban growth, (2) forest conservation, and (3) integrated (
Table S1 in the Supplementary Materials). These planning preferences were selected to enable the comparison of the economic and environmental components of sustainability with an integrated perspective that seeks to balance land uses that influence components of economy, environment, and society. Further, the integrated scenario follows a current and realistic pattern of urbanization. These scenarios were based on a hypothetical metropolitan region that occupies approximately 220,000 ha with specific land uses distributed across the area based on a desired composition that emphasizes a planning preference. The first scenario was composed of 70% urban surface, 20% agriculture, and 10% forest. This scenario is based on the current perspective of countries where economic and urban growth is a typical political priority (e.g., USA, China, and Mexico). Consequently, there is a reduction in natural capital and environmental degradation is more likely to occur. The second scenario was composed of 30% urban surface, 25% agriculture, and 45% forest. This scenario has an environmental emphasis, where conservation of forest and natural resources are the priority. For example, northern European countries (e.g., Denmark and Sweden) encourage strategies to prioritize natural capital. However, the proposed distribution of land use types in this scenario may be unrealistic in an already highly urbanized area where many stakeholders may not support halting economic development due to common perceptions that link economic progress to the advancement of human welfare. The third scenario was composed of 50% urban surface, 30% agriculture, and 20% forest. This scenario tries to integrate the three components of sustainability but considers the inevitability of urbanization while global population continues to increase. We expect that this option is the most suitable for present and future human societies, at least for the foreseeable future. Additionally, we expect that, despite a continuing increase in urban surfaces, every scenario can improve human wellbeing alongside progress towards sustainability. Therefore, three additional scenarios were created using the same three land use-composition scenarios described above but examining the effect of agriculture lands primarily driven by monoculture crop cultivation versus a diversity of crop cultivation that emphasizes nutrition through macronutrients and micronutrients for optimal health. These three additional scenarios used the same set of ten crops but the percentages for each crop varied according to the planning preference. For example, the first scenario used 2% of agricultural land for each crop (total of 20% agricultural lands), the second scenario used 2.5% per crop (25% agricultural lands), and the third scenario used 3% for each crop (30% agricultural lands).
2.4. Data Collection and Preparation
Boundaries of the synthetic landscapes were based on the Toluca Metropolitan Zone (TMZ), which is a large and rapidly urbanizing metropolitan region in the center of Mexico. This region was selected because it is within Latin America, a part of the Global South [
31], and due to its similarity to many other rapidly urbanizing contexts in the region. Additionally, this region is characterized by heterogenous landscapes, soil types, climate, and elevation, mixed with disparate land use patterns and agricultural priorities. This complex system, representing multiple tradeoffs, was used as a foundation for the synthetic landscape research. These synthetic landscapes were constructed using a classified 2015 Landsat image of this area. The original classified image contains 5 main classes: (1) forest, (2) water, (3) agriculture, (4) bare surface, and (5) urban surface. Each synthetic landscape is a modified version of the original classified image, where the number of pixels contained within each class, i.e., the percentage coverage of each land use category, varies in accordance with the sustainability component emphasis (economic and urban growth, forest conservation, or integrated). The algorithm to develop each synthetic landscape starts with the original classified image. A frequency distribution of pixels classified as forest, grass, agriculture, and urban is then computed. The number of pixels that need to be converted from one class to another (e.g., agriculture to urban) is then computed using the desired distribution of urban, agriculture, and forest land uses (e.g., 70% urban surface, 20% agriculture, and 10% forest). The process of converting pixels from one class to another cannot be performed sequentially starting at one corner of the image, as this leads to linear land use classification that does not represent realistic geographic patterns. To better control the placement of converted land use classes across the study area, the algorithm randomly chooses starting points in the image and converts a user-specified number of pixels from one land use class to another. This number ranged from 500 pixels to 10,000 pixels at a time. This procedure of selecting a random starting point is repeated until all conversions are complete. Pixels classified as bare surface or water are excluded from modification and remain consistent across all the synthetic landscapes.
The forest class was dominated by coniferous (mostly pine trees) and deciduous forest (mostly oak trees). The water class comprised small bodies of water, peatlands, and dams. The agriculture class included both irrigated and rain-fed crops, but mostly dominated by irrigated maize. The urban surface consisted of built-up area, roads, and pavement.
Agricultural land has been increasingly dominated by monocultures (e.g., maize) in Mexico; therefore, we developed similar synthetic monoculture landscapes, but also developed diverse crop landscapes with an emphasis on nutrition, as explained in the section above. The crops included 1 fruit type; 5 vegetable types, including a leafy green and a brassica; 2 types of legumes; 1 type of nut; and 1 type of grain. Specifically, we used the following crops: apple, barely, cabbage, cucumber, fava beans, maize, peas, pecan nuts, spinach, and tomatoes (
Table S1 in Supplementary Materials). These 10 crops were selected based on climate and geography and following the food groups recommended in a whole food plant-based diet (WFPBD). The WFPBD has been demonstrated to be a health-promoting diet that can meet most of the nutritional requirements of the population. Additionally, evidence shows that a WFPBD may help with the prevention and treatment of chronic diseases, such as type-2 diabetes and cardiovascular diseases [
32,
33,
34].
2.5. Modelling Ecosystem Services Using InVEST
Two InVEST models were used to evaluate the impacts of planning preferences on ecosystem services: crop production percentile (CPP) and sediment delivery ratio (SDR). Both models require spatial data projected in a coordinate system with units of meters. A spatial reference is important in the CPP model because crop yield is primarily driven by climate, while the SDR model uses global sediment database [
28]. We selected the latitude of the Toluca metropolitan area in Mexico. Spatial files were projected in UTM coordinates (WGS 1984 Zone 14).
The main data needs of the CPP model include projected land use or land cover raster, and a comma-separated values file containing the name of the crop as specified by the model with the respective land use code. The main outputs include total observed production, total area covered by each crop, and nutrient information for each crop (e.g., energy, vitamins, minerals). The main data needs for the SDR model include a projected land use or land cover raster, topography, rainfall erosivity index, soil erodibility, a comma-separated values file containing management factor and support practice factor associated with the land use, and a watershed polygon obtained with the InVEST “DelineatIT” tool. This tool uses the deterministic eight-neighbor (D8) method to route flow directions using the digital elevation model (DEM). The D8 approach determines the boundary of the watershed, river networks, and sub-watersheds associated with each river segment [
35]. The main outputs for the SDR model include total amount of potential soil loss, total amount of sediment exported to the stream, and total amount of sediment deposited to the landscape. For the full list of data needs and main outputs for each InVEST model, see
Table S2 in the Supplementary Materials.
Spatial data were obtained from “Instituto Nacional de Geografía, Estadística e Informática” INEGI and “Consejo Nacional de Biodiversidad” CONABIO [
36,
37]. Data included soil types, digital elevation model, and land use. Soils in Mexico use the 1974 FAO-UNESCO classification. As erodibility was not calculated for the area, we conducted a literature review to obtain characteristics and properties for each type of soil. Percentage of sand, silt, clay, and organic matter for each soil was obtained from INEGI [
36]. The K factor was calculated following the InVEST 3.9.0 user’s guide [
28]. The R factor was obtained from European Soil Data Centre, which provides global rainfall erosivity data [
38].
2.6. Nutritional Needs Assessment
Data generated by the crop production model was used to compare nutrition provided by monoculture cultivation (generated in the initial land use scenarios) and crop diversity cultivation (generated in the land use scenarios recommended for nutrition). The model’s output included the total area covered by the crop, annual crop production and resultant nutrients. The nutritional data were used to calculate the population size that could be supported by each scenario on an annual basis. Those values were divided by 365 to obtain daily nutrition production, which was subsequently divided by the Recommended Dietary Allowance (RDA) or Adequate Intake (AI) for each macro- and micronutrient to estimate the population size that could be supported under each scenario. Both RDA and AI set reference values to meet nutrient requirements of healthy people (NIH, 2021), but because nutrient requirements vary for children, women, men, and elderly people, we used the upper limit RDA or AI for each nutrient to ensure all healthy people were supported despite gender and age, yielding a conservative estimate of the population size supported for each scenario.
4. Discussion
This paper proposes an urban nutrition simulation extension of the urban ecological economic system proposed by Huang and others in 1988 and describes how the InVEST software model was used to assess the interactions between urbanization, agricultural production techniques, and ecosystem services, to better understand their impacts on local nutritional security. These components were analyzed using synthetic landscapes and crop production systems on (1) provisioning and regulating ecosystem services, such as food production and erosion control, and (2) calories and nutrients to support the population. This research contributes to different themes in the ecosystem services literature including human welfare, land use, agroecosystems and food security, and landscape planning [
20]. Furthermore, this paper explicitly integrates the concept of sustainability in a world where urbanization, land use change, economic development, and population growth are inexorable and impact nutritional security. The synthetic landscapes along with the InVEST models illustrated how different planning preferences and agricultural choices impact ecosystem services and nutritional output that consequently affect urban sustainability. Results derived from these simulations showed that: (a) provisioning ecosystem services and human nutrition can be improved when land use/land cover (LU/LC) are balanced as shown in the integrated scenario; (b) crop diversity techniques yield a higher production compared to monocultures; (c) larger populations can be supported when adding a variety of crop types, even in scenarios where there is limited cropland area and urban surface dominates; (d) nutritional benefit, at minimum, is at least two times higher under the crop-diversity priority compared to monoculture; and (e) limiting nutrients can be identified and prioritized to ensure the population meets all nutritional requirements for optimal human wellbeing.
The results obtained from the InVEST models demonstrated that the integrated scenario was the most suitable land use-planning preference to increase crop production and enhance nutrition. However, it is important to note that the forest conservation scenario had the least potential soil loss compared to the other scenarios and, therefore, was most likely to preserve ecosystem services. Crop production under the economic and urban growth scenario supported the smallest population and induced a large amount of erosion. This could be considered the least favorable approach for optimal nutrition and the conservation of ecosystem services from a sustainability perspective. It is important to note that the integrated scenario showed the least sediment retention and highest soil export, even compared to the economic and urban growth scenario. This can be due the high proportion of impermeable surface found in the urban growth scenario. After other types of land use (e.g., agriculture and forest) are converted into concrete, soil is no longer lost. However, these land transformations require construction efforts that temporarily lead to high rates of sediment loss and impact other ecosystem services, including habitat quality. Additionally, cover management and support practice factors impact soil loss and greatly vary from one jurisdiction to another. It is likely that C and P factors are higher in integrated and conservation scenarios, resulting in lower levels of soil export and loss; however, these values were held constant across scenarios. Sustainable agricultural practices (i.e., crop diversification) are necessary to decrease soil erosion and support other ecosystem services.
In addition to assessing crop production and soil erosion in landscapes that prioritize corn monoculture, crop production and derived nutrition was also assessed in landscapes with diversified crops. The integrated scenario demonstrates that locally grown crops provide adequate nutritional support to urban communities, thereby reducing reliance on food imports. This option is relevant at the local scale when desired or mandated changes to promote agricultural land use are not feasible. In response, current practices encourage food imports, typically resulting in lower food prices. However, food imports that prioritize lower costs are often limited in variety and poor quality, which in turn can adversely affect human wellbeing, particularly for vulnerable groups. Further, there are many costs associated with poor quality food, i.e., malnutrition, cardiovascular diseases, and healthcare costs. Current food systems also encourage overproduction, which has led to a drop in commodity prices. Hence, commercial large-scale agribusinesses control the market and local small farms suffer the consequences of this unsustainable activity [
39,
40]. Monoculture production systems, driven by industrialization and technological developments, enable the overproduction of valuable commodities. Corn, for instance, is a versatile commodity used for livestock feed, food, seed, biofuel, and industrial purposes [
41]. The current demand and profitability of corn production has resulted in a lack of diversity of crops. The ecological footprint of corn monocultures includes food waste; soil erosion and nutrient losses; reduced biodiversity and increased vulnerabilities to pests; the use of large amounts of herbicides, particularly atrazine, which is known to be an endocrine disruptor; the use of large amounts of nitrogenous fertilizer, which leads to surface and groundwater pollution; as well as air pollution that contributes to climate change [
42,
43].
An agricultural system that encourages crop diversity not only reduces these environmental impacts, but substantially increases the nutrition in communities, supports local farmers and reduces overproduction and food waste.
Figure 3 showed the different levels of population size supported across the three planning scenarios, and the type of production system. Nutritional output, and therefore the population size supported, was the highest under the integrated scenario when prioritizing crop diversity. In fact, crop diversity always improves nutritional output regardless of land use distributions. The graphs also revealed that, as you transition across scenarios, the percentage increase associated with crop diversification is always greater than the percentage increase associated with monoculture, i.e., the slopes associated with crop diversification are greater in comparison with monoculture. This suggests greater incentive to transition from the urban and economic growth scenario to the integrated scenario under crop diversification strategies.
The crop selection used in this research (
Tables S3–S5 in Supplementary Materials) showed that communities can obtain most required nutrients from plants. Additionally, there are other phytonutrients obtained from plants that improve health. Nutritional transition has impacted the health of the population by contributing to non-communicable diseases such as heart disease, obesity, diabetes, and cancer, which have become leading causes of death globally [
44,
45]. This is important in an urbanized and industrialized world where processed and animal-based foods are predominant [
46]. There are several studies that show that whole-food plant-based diets protect against many types of cancers, such as breast, prostate, colorectal, and gastrointestinal cancers [
47,
48]. They also are a safe, low-cost, and effective option to lose weight, reduce cholesterol, and prevent and treat cardiovascular diseases and type-2 diabetes [
33,
34,
49,
50]. Lycopene, for instance, a carotenoid present in tomatoes, is associated with a reduced risk of cardiovascular, and degenerative, diseases, and has shown to defend against cancer development and progression [
51,
52,
53]. Hever suggests that a health-promoting plant-based diet should be composed mainly of vegetables, fruits, whole grains, legumes, herbs and spices, and in smaller amounts, nuts and seeds [
54]. These food groups contain all the nutrients necessary for optimal wellbeing. The only exception is vitamin B12; animals are the exclusive natural sources of active cobalamin [
55].
Decision makers can use the InVEST CPP model and simulation landscapes to explore differences in nutrients resulting from distinct land uses and crop production systems. The NII developed in this research can determine how much benefit can be obtained when switching scenarios of crop production and land use. It is also possible to determine the limiting nutrient affecting a community. For example, if fat is insufficient, decision makers can focus on this limiting nutrient and encourage healthful whole food sources of fat such as seeds and nuts. Additionally, the NII can provide guidance on how much of this nutrient is needed to support a larger population.
The UN extension along with the InVEST modeling guided the exploration of the effects of land use change, crop techniques, and ecosystem services on nutrition. Decision makers can use this guide to identify nutritional needs and pathways that can be taken (e.g., land use planning and agricultural production techniques) to improve nutritional yield in a sustainable manner. However, there are some limitations when using the InVEST model. The CPP model only provides estimates from the 175 crops (based on FAO national datasets) and many crops have incomplete nutritional information. Research presented in this paper utilized the best combination of crops in the Toluca, Mexico region according to the climate bin maps provided by the model. Additionally, spatial resolution might be another factor that affects overall results. This research took a conservative approach to model the population size that can be supported under each scenario by using the highest recommended dietary allowance, which varies among men, women, and children. If average RDA were used as the basis, or if a representative population distribution of sex- and age-based RDAs were modeled, the population size that could be supported would be greater. Despite these limitations, the framework should be valuable to stakeholders and decision makers since it provides guidance on how to preserve ecosystem services while increasing human wellbeing. Protecting ecosystem services and including them in planning is crucial to promote sustainable cities [
56]. This paper focuses on soil erosion and crop production; however, other types of ecosystem services can be considered. For example, urban ecosystem services can be used in the framework to foster resilience in urban settings [
57]. Sustainability is a work in progress, an evolving concept that has different emphases in different countries and at different points in time. However, one of the shared challenges is to balance different components (environment, society, economy) in a world that is continually transforming and becoming increasingly urban.