An Integrated Food, Energy, and Water Nexus, Human Well-Being, and Resilience (FEW-WISE) Framework: New Mexico

Interconnected food, energy, and water (FEW) nexus systems face many challenges to support human well-being (HWB) and maintain resilience, especially in arid and semiarid regions like New Mexico (NM), United States (US). Insufficient FEW resources, unstable economic growth due to fluctuations in prices of crude oil and natural gas, inequitable education and employment, and climate change are some of these challenges. Enhancing the resilience of such coupled socio-environmental systems depends on the efficient use of resources, improved understanding of the interlinkages across FEW system components, and adopting adaptable alternative management strategies. The goal of this study was to develop a framework that can be used to enhance the resilience of these systems. An integrated food, energy, water, well-being, and resilience (FEW-WISE) framework was developed and introduced in this study. This framework consists mainly of five steps to qualitatively and quantitatively assess FEW system relationships, identify important external drivers, integrate FEW systems using system dynamics models, develop FEW and HWB performance indices, and develop a resilience monitoring criterion using a threshold-based approach that integrates these indices. The FEW-WISE framework can be used to evaluate and predict the dynamic behavior of FEW systems in response to environmental and socioeconomic changes using resilience indicators. In conclusion, the derived resilience index can be used to inform the decision-making processes to guide the development of alternative scenario-based management strategies to enhance the resilience of ecological and socioeconomic well-being of vulnerable regions like NM.


INTRODUCTION
The availability and consumption of food, energy, and water (FEW) resources heavily rely on one another, and there are numerous ways in which these three systems overlap and interconnect. For example, food production accounts for ∼30% of global energy consumption and ∼92% of the human water footprint (Finley and Seiber, 2014); 15% of global water withdrawals are used for energy production (IEA, 2012) and 70% for food production (Ritchie, 2017); 3% of global electricity is required for water provision (Liu et al., 2015) and 6% for food production (Ritchie, 2017). However, currently, about 0.8 billion people are hungry (FAO, 2017), 1.4 billion people have no access to energy (Alstone et al., 2015), and 0.7 billion people have no access to water (FAO, 2011). The demand for FEW resources is projected to increase, further adding management, sustainability, and resilience challenges as there is a need to increase food production by ∼70% (FAO, 2014;World Bank, 2016) along with 57% more water (WWAP, 2015) and 40% more energy (OECD, 2012;World Energy Outlook, 2014;Rasul, 2016) in the next 20 years. The combination of factors that include global population growth, which is projected to grow from 7.7 billion people in 2019 to 9.7 billion in 2050 and 10.9 billion in 2100 (UNDES, 2016), increased economic development, and rapidly changing climate amplifies the scarcity of water supply, the decline in food production, and the depletion of fossil fuels that humans heavily rely on for socioeconomic development (Gerland et al., 2014;Karandish and Mousavi, 2016). Providing a framework that allows proper modeling and evaluation of FEW system components is a first step toward improving the ability to better manage and allocate these resources.
These factors which can generally be grouped into two categories-socioeconomic and climate change-act as external drivers of vulnerability of FEW systems to related disturbances and shocks (Wisner et al., 2004). Major challenges to FEW systems are those imposed by climate change-induced extreme events, such as droughts , increased temperature, and variable precipitation patterns, especially in semiarid regions with low adaptive capacity [e.g., New Mexico (NM) in the Southwest United States (US)] (Niang et al., 2014). Increased temperatures and declining precipitation can increase the demand for water supply-mostly for irrigation to ensure food security (FAO, 2012), lead to environmental degradation, and eventually deteriorate human livelihood, well-being, and economic development (Sun and Yang, 2016). Additionally, significant socioeconomic changes (e.g., food and energy prices) can influence the availability and accessibility of FEW resources, making large portions of the population unable to afford and access these basic needs (Mohtar and Daher, 2010;Ringler et al., 2016). The effects of these external drivers are complex, span the three FEW nexus components, and can drastically affect the sustainability and resilience of FEW resources.
As it is important to provide an improved understanding of the functionality of individual components of FEW systems, it is also critical to have an elaborate characterization of the interconnections and behavior of these systems as a nexus as they evolve and transform these climatic and socioeconomic setbacks (Grafton et al., 2017;Rosa et al., 2017;Scanlon et al., 2017;D'Odorico et al., 2018;Nhamo et al., 2018;Rosa et al., 2018). These interconnections are among the most important ones of nature that are essential for human well-being (IRENA, 2015;EIA, 2018;Nhamo et al., 2018;Ding et al., 2019). A distinct framework is required to investigate FEW system dynamics and assess the resilience of such coupled natural-human systems (Carpenter et al., 2001;Gunderson and Holling, 2001;Holling, 2001;Folke, 2006). The scope of this conceptual analysis was to develop a framework to assess FEW system relationships, identify external drivers, and conduct a threshold-based resilience evaluation.
The resilience challenges related to FEW systems and human well-being are particularly amplified in semiarid regions such as the state of NM in the southwest of the United States. NM, with its unique characteristics, can be considered as an informative case of FEW nexus and human well-being. NM has recently experienced frequent and prolonged drought events , and it is projected to experience extremely high-water scarcity conditions in the future (WRI, 2015). NM is a major fossil fuel (e.g., crude oil and natural gas) producer, and the energy sector significantly supports its economy. The resilience of NM's FEW nexus is challenged due to the increased pressure on its water resources, food production Zaied et al., 2019;Zaied et al., 2020), economic shocks from energy markets, and fluctuations in ranchers' and farmers' income. These NM challenges further highlight the need to appropriately conduct qualitative and quantitative analysis, and integrate and evaluate the dynamic behavior of FEW systems to enhance their resilience.
This study used the "nexus" concept because it has emerged as an effective means to address these challenges by properly describing the complex linkages and inseparable interactions between multiple distinct but interconnected systems (McGrane et al., 2019). The concept has long been used in philosophy, cell biology, and economics. It was introduced to natural resources disciplines in the Food-Energy Nexus Programme (Sachs and Silk, 1990), to provide integrated solutions to food and energy scarcity (McGrane et al., 2019). It was highlighted at the Bonn 2011 Nexus Conference to promote the understanding of FEW resources, provisioning basic needs, and ensuring security (Hoff, 2011). The nexus concept can be defined as an integrated systems approach to qualitatively describe and identify interconnected subsystems and to quantitatively account for their functionality and resource availability, enhance their synergistic use, and minimize tradeoffs with the main goal of sustaining the well-being of human societies and resilience of FEW systems. The concept has been explored to promote effective management of FEW resources qualitatively and quantitatively (Hoff, 2011;Fischer et al., 2015;Keairns et al., 2016;Scanlon et al., 2017), for example, to assess the effects of irrigation practices on water resources (de Vito et al., 2017;Cai et al., 2018), assess human livelihood under climate and socioeconomic impacts (Ding et al., 2019), and integrate FEW and human systems to promote environmental security and sustainability (Biggs et al., 2014;Biggs et al., 2015;de Grenade et al., 2016).
Most previous research efforts that focused on qualitative and (to some extent) quantitative assessments of the FEW nexus lacked modeling of its dynamic behavior (Keairns et al., 2016;Albrecht et al., 2018;Givens et al., 2018). A clear quantitative description is required to adequately integrate such coupled natural-human systems, develop effective management decisions (Ramankutty et al., 2018;Rosa et al., 2018;Mohtar and Daher, 2019), and improve resource use efficiency and synergies (Kan et al., 2016;Fang and Chen, 2017). However, most traditional modeling approaches may not be suitable for such complex systems (CE, 2010;Bazilian et al., 2011;ADB, 2013;WB, 2013;NSF, 2014) because their parametrization limits their application to local scales (Miralles-Wilhelm, 2016); hinders their ability to integrate more than two subsystems (Dubreuil et al., 2013); and lack in coupling socioeconomic systems. On the other hand, these limitations provide the opportunity to apply system dynamics (SD) (Forrester, 1994;Sterman, 2000) which allows characterizing the dynamic behavior systems. The SD approach is unique for integrating physical (e.g., FEW) and socioeconomic (e.g., employment) systems (Bagheri and Hjorth, 2007;Winz et al., 2008;Tidwell et al., 2018). An SD model can help visualizing the interrelationships between discrete sectors, the dynamic changes, and the interdependences among FEW systems (Simonovic and Fahmy, 1999). It also allows the use of scenario-based analysis to simulate past, current, and forecast future resource availability and consumption (Fiksel, 2003;Bieber et al., 2018;Laspidou et al., 2020).
It is important to monitor the dynamic behavior of the FEW nexus over time in response to disturbances. Ecological (e.g., FEW) and social (e.g., human society) systems can, mostly, adapt to external stresses over time until a resource availability threshold is surpassed, beyond which these systems can undergo significant shifts-either transform into a new equilibrium state (Angeler and Allen, 2016) or become nonresilient (Angeler and Allen, 2016). The temporal behavior of these systems can be effectively captured using the resilience concept that was introduced in the field of ecology beginning in the 1970s (Holling, 1973). Resilience thinking is a generic approach that has been increasingly adopted to understand socio-ecological systems (Carpenter et al., 2001;Folke et al., 2002;Folke, 2006). Achieving resilience across FEW and human systems signifies the amount of disturbance that these systems can withstand before shifting into a new stable state. It is important to identify the factors that allow maintaining resilience and those that can introduce risks and undesired resource scarcities.
Additionally, performance indicators that couple FEW and human well-being (HWB) systems need to be developed. HWB is a subset of social well-being and economic growth that can be measured using a number of indicators such as employment, education, and income, among others (OECD, 2012). FEW and HWB systems can seamlessly be integrated using their relevant indices-FEW Index and Human Well-being Index (HWBI), respectively. The goal of this analysis was to develop an integrated framework referred to as FEW-WISE including FEW, human Well-beIng, and reSiliEnce. The objectives were to 1) provide an improved characterization of the integrated relationships in FEW systems, 2) provide a quantitative assessment approach for the FEW nexus based on relevant indices and thresholds that can identify its dynamic equilibrium status, and 3) propose a resilience simulation framework to assess the response of FEW resources in the presence of external stresses. These objectives can be achieved using the FEW-WISE framework that consists of five steps (FEW-WISE Framework section below). In this article, the first two steps were described in detail to lay out the baseline information needed to model and conduct indicators-based analysis. The other three steps were individually conceptualized and will be described in separate modeling, indices, and resilience assessments. The SD modeling results, predicted parameters of FEW nexus, and resilience index will be reported separately in a follow-up article that shows the application of the FEW-WISE framework. This framework was developed within the context as part of an INFEWS project using NM as a case study but can be generalized for other vulnerable regions with similar conditions.

STUDY AREA
NM encompasses a large geographic area with diverse interior-continental ecosystems, including mountain ranges, forests, grasslands, and deserts (NOAA, 2017). NM population has shown a sustained growth during 2001-2005 after a leveled growth around the 2000 Census (BBER, 2008). However, NM has experienced a slow growth in population during 2010-2019 when most of the counties saw a decrease in population (NMEDD, 2010); migration and economic trends provide real historical impetus toward population dynamics (USBC, 2019). NM has a considerable reserve of fossil fuel, mineral, and renewable energy resources. Crude oil and natural gas place the state among the top 10 energy producing states in the United States (EIA, 2019; EIA, 2020). NM's abundant land makes it the fifth largest state in the United States with 314,850 m 2 of land of which 175,230 m 2 are classified as farms and ranches with 88% of this area identified as rangeland ( Figure 1) (Goodwin and McDermott, 2017).
NM has an arid to semiarid climate and is ranked as the fifth driest state in the United States. NM average annual precipitation is ∼380 mm (1895-2019) and ranges from less than 254 mm over much of the southern region to more than 508 mm at higher elevations (NOAA-NCEI, 2020). NM temperatures vary widely with a monthly average temperature in the northern mountainous regions ranging from 6.7°C in January to 15.6°C in July, while in the lower elevations in the south, the range is from 4.4°C in January to 26.7°C in July. The state has a limited and variable water supply due to its normally dry conditions and frequent droughts that add challenges to explore its full potential in food and energy resources, thus affecting the sustainability of these systems and New Mexicans' livelihood and well-being.
The major industries in NM include agriculture and energy. Within the agricultural sector, NM is known for its considerable production of livestock (e.g., beef), diary (e.g., milk), and crop (pecan and hay) commodities. The livestock industry is an integral economic component for NM that adds ∼ $2.5 billion to its gross state product (GSP) and employs over 32,000 workers as of 2012 (Diemer et al., 2012). Most irrigated lands are located along the Rio Grande River corridor running from north to south through the middle of the state. The energy sector provides significant revenue to the state, which then supports New Mexicans in different ways (e.g., schools, hospitals, and state government and other services) (EMNRD, 2015). NM has plentiful energy resources (fossil fuels and renewables) that make it the eighth largest energy producer in the United States with more than 6 and 4% of the United States total proved crude oil and natural gas reserves in 2018, respectively (EIA, 2018). This wealth of energy resources also creates economic development opportunities, from attracting manufacturing to additional opportunities for energy exports. Both food and energy production are water-intensive activities.

FEW-WISE FRAMEWORK
The FEW-WISE framework follows five steps: 1) qualitatively identify the relationships between FEW nexus components and quantitatively assess resource exchanges, 2) identify the drivers of the systems, 3) integrate FEW systems using SD models to evaluate their response to the identified drivers, 4) develop FEW and HWB performance indices, and 5) develop a resilience monitoring criterion using a threshold-based approach that integrates these indices ( Figure 2). The qualitative and quantitative assessments of NM's FEW nexus along with the external drivers can be used as inputs to simulation FIGURE 1 | Study area showing county-wise distribution of food-energy-water systems and climate over New Mexico that includes livestock production (NASS, 2015), rangeland cover (USFS, 2019), total crops sold (NASS, 2015), total water withdrawal (Molly et al., 2015), Oil & Gas production in Oil Conservation Division (OCD) Districts, and30-year (1981-2010) mean annual temperature and rainfall (PRISM, 2010).
Frontiers in Environmental Science | www.frontiersin.org June 2021 | Volume 9 | Article 667018 4 and prediction of FEW systems' temporal response to identified changes and can be used to develop the corresponding FEW and HWB resilience indices at the county and state levels. These five steps are discussed and explained in more detail in NM FEW Nexus Description to Resilience section below.

NM FEW Nexus Description
To operationalize the analysis of NM FEW systems, the interactions between FEW subsystems under external drivers were first qualitatively characterized, as depicted in Figure 3. Such depiction was needed to promote cross-sector collaboration, develop coherent decisions, and enhance HWB and resilience (Keskinen et al., 2016). This qualitative description of NM FEW systems was based on resource use, for example, food commodities produced using NM land, water, and energy. Much of the food commodities produced in NM are exported, and thus, additional supplies are imported to meet local demands.

Food and Energy
In NM, energy is used in harvesting, production, processing, and transportation of agricultural (e.g., crop and livestock) commodities. Likewise, agricultural by-products (e.g., biofuels and biopower) can be used, to a limited extent, to generate energy. The important variables that link these two sectors were identified based on a few factors including the total amount and economic value of production as well as the amount of resources used or needed from the other sector to produce these variables.
NM food production systems include crops and livestock commodities. Major crops include field crops (e.g., hay, winter wheat, sorghum, corn, cotton, peanuts, and dry beans), vegetables, and nuts (e.g., chili, onions, and pecans). In 2017, the total crop production in equivalent metric tons was about 4.012 million (NASS, 2018). Corn (grain and silage) accounted for about 53%, followed by hay (all types) and onions about 28 and 6% of the total production, respectively. While pecans accounted for only 1% of the total crop production by weight, its economic value outpaced that of corn and hay ( Figure 4). NM FIGURE 2 | Integrated FEW-WISE framework that couples food-energy-water systems and drought impact, socioeconomic, and human well-being components.
FIGURE 3 | New Mexico food-energy-water (FEW) nexus describing the relationships and interconnections of its components under climate and socioeconomic drivers.
Frontiers in Environmental Science | www.frontiersin.org June 2021 | Volume 9 | Article 667018 5 was among the largest producers of pecans (a high value cash crop) as the state ranked second in 2017, with a 28% increase from 2016. In 2018, NM surpassed Georgia in pecan production for the first time and ranked as the top producer nationally. Most of the pecan acreage is in Doña Ana county, followed by Eddy and Chavez counties (NMDA, 2018).
NM livestock production systems consist of beef cattle, dairy, sheep, and lambs. Except for dairy, most of NM livestock production depends on pasture and rangelands. Based on January 2018 inventory, the total number of all cattle and calves was ∼1.51 million and that of milk cows was ∼483 thousand. The economic value of dairy production (mainly milk) was more than that of all meat animals as they accounted for 41.4 and 31.5% of the total cash receipts, respectively, of all agriculture commodities (NASS, 2018). Nationally, in 2012, NM was ranked ninth for the value of milk sold and the number of milk cows. Most of this production occurs in Chavez, Curry, and Roosevelt counties.
In terms of energy production, NM has a diverse energy portfolio, and it is considered a major producer of crude oil and natural gas as it was ranked third and eighth in the United States in 2020, respectively. NM natural gas and crude oil account for close to one-tenth and over 3% of the United States total production, respectively. Major crude oil and natural gas reserves are in the Permian Basin in the southeast and in the San Juan Basin in the northwest (see oil and gas production districts map in Figure 1). The Permian Basin has supplied more than 5.7 billion m 3 of crude oil and about 3.75 trillion m 3 of natural gas as of January 2020. As of 2018, the U.S. Energy Information Administration (EIA) estimated the remaining proven reserves in the Permian Basin exceed 11 billion barrels of crude oil and 1.38 trillion m 3 of natural gas, making it one of the largest hydrocarbon-producing basins in the United States and globally (EIA, 2019). In 2010, EIA ranked the combined San Juan Basin Gas Area in Colorado and NM as the second largest natural gas field in the United States in terms of proven reserves, with a production of 54 billion m 3 in 2019. Other sources of energy in NM include coal and renewables (e.g., geothermal, hydroelectric, biomass, wind, and solar) but have minimal contribution to the state's total energy production. NM is the seventh largest net supplier of energy in the United States, with a total annual production of 2,820 trillion Btu in 2017 (mostly crude oil and natural gas).
The total demand for energy by the agriculture sector was 1,700 trillion Btu in 2009 (all sources of energy), which decreased to 1,500 trillion Btu in 2012, and has been increasing since 2014 reaching 1,714 trillion Btu accounting for about 1.74% of the total U.S. primary (e.g., coal, natural gas, and oil) energy consumption (Hitaj and Suttles, 2016). The secondary form of energy (e.g., electricity) is also used at all the stages of food production, including pumping water for irrigation (7.5% US cropland and pastureland were irrigated using electricity in 2007), powering tractors for tillage and harvesting; transporting and distributing food products, and for heating and cooling in livestock activities (23.7% in the year 2002) (Miranowski, 2005;Nord et al., 2005).
Fossil fuels are used throughout the feed chain of the livestock production systems (production, transportation, storage, processing), farm operations (machinery, equipment, climate control), and products chain (transportation, processing, storage) (Sainz, 2003). (Patrick, 1977) estimated the energy demand for dairy, range beef, and feedlots beef productions as 57.56 million Btu, 57.17 million Btu, and 39.30 million Btu per animal unit, respectively. NM dairy production systems consistently showed the highest energy demand. NASS (2018) reported 465,000 beef cattle head and national average beef cow weight as 608 kg for the year 2017 that is equivalent to 6,28,215 animal unit (AU). Energy requirement of beef cattle in the year 2017 is approximately 60,652 billion Btu.
Energy requirements of alfalfa, corn, wheat, and sorghum were estimated using energy requirements reported by Patrick (1977) using NM crop yields of the year 2018. The energy requirements  Frontiers in Environmental Science | www.frontiersin.org June 2021 | Volume 9 | Article 667018 6 of crop (production) for alfalfa (4,230 kg), corn grain (4,712 kg), sorghum (957 kg), and wheat (405 kg) are expected to be around 15.8, 27.4, 8.0, and 4.2 million Btu, respectively. Overall, NM energy consumption varies widely by sectors defined by the EIA ( Figure 5). NM agricultural industries spent nearly $90 million for petroleum fuel and oils and another $49.6 million for electricity in 2016 to power farm equipment, manufacture fertilizer, and transport products (CEA, 2016). The transportation sector is the largest energy consumer in NM, followed closely by the industrial sector, where more energy is consumed per capita than in three-fourths of the United States (EIA, 2017) ( Figure 5).

Water and Energy
Generally, energy production relies on water availability, and the supply and distribution of water require energy (King et al., 2008;Sanders and Webber, 2012). In NM, extraction of energy resources including crude oil and natural gas involves the use of hydraulic fracturing technology-a process of high-pressure injection of a fluid mixture that consists of water, sand, and chemicals into bedrock formation for increased production. Also, power generation is a major source of water consumption (e.g., mainly electricity) for cooling (specifically coal and natural gas). There are no nuclear power generation plants in NM. Minimal amounts of water are used in other mining activities. On the other hand, energy is needed for pumping, conveyance of water for agriculture (from surface and/or groundwater resources), purification, treatment, and distribution of freshwater and wastewater, and other domestic uses at the household level (e.g., water heating and laundry).
To account for how much water is used for energy production, it is important to evaluate the exchange of quantities between these two systems. NM water supply depends on surface and groundwater resources relatively equally to meet its demands. During a typical year, almost half of NM water comes from groundwater aquifers (Longworth et al., 2019)-about 2.1 billion m 3 per year, pumped from the five major freshwater aquifers underlying NM. More groundwater pumping may occur during drought years to make up for the deficit in surface water supply. Groundwater withdrawal is more energy-intensive than that of surface water. NM uses about 12.87 billion m 3 per day (Maupin et al., 2010) as mining and irrigation account for 1 and 76%, respectively, of total groundwater use (Longworth et al., 2013) (Figure 6). In NM, about 1.65 billion m 3 of water were applied to irrigate 2.73 billion m 2 in 2018 (US DOE, 2014).
The growing practice of multistage hydraulic fracturing has increased the amount of water used per well in recent years, which, in NM, can range from 1.9 to 30 million liters per well. While this amount of water used for energy extraction (crude oil and natural gas) represents a small fraction of the total water use, it can trigger local stress on freshwater supply in high-production areas. Withdrawal in all water use categories combined was about 3.84 billion m 3 , of which surface water and groundwater accounted for 52.34 and 47.66%, respectively. In 2015, water use under the mining category accounted for 52.  (Longworth et al., 2013). However, in the United States, irrigation accounted for 42% of freshwater withdrawal in 2015 as compared to 84% in NM. Therefore, the data about the nature and volume of the aquifers are important to manage NM water resources, economics, and environment.
To represent the variation in water use for energy production, water intensity of the power generation variable can be developed and used to account for the combined water use in mining and power generation. For the United States, water intensity of the total power generation was about 57.2 L per kilowatt-hour in 2014 (49.2 L per kilowatt-hour in 2017). Water intensity referred to here is the average amount of water withdrawn per unit of total net electricity generated. The total energy used for irrigation and sprinkler operation accounted for 1,160 GWh or approximately 3% of NM's total energy use in 2011. This estimate was based on acres of land using groundwater and sprinklers. In 2010, the total water used for energy production accounted for about 3% of NM's overall withdrawal from surface and groundwater sources (Tatro, 2018). The amount of energy needed for public supply was about 26 MWh of electric power for groundwater pumping from an average depth of 150 m.

Water and Food
NM's water, which supports its crop and livestock production and food processing, is limited due to variable precipitation and recurring drought conditions that consequently create production uncertainties, inhibit farmer and rancher livelihood, and can result in more fallow land or land transitions out of agriculture. For the last half century, water use by agriculture was roughly evenly split between groundwater and surface water. However, in 2013, out of 1.62 billion m 3 irrigation water, 0.99 billion m 3 ground water from wells, 0.18 billion m 3 on-farm surface water, and 0.45 billion m 3 off-farm water from other sources were used (USDA, 2013). In addition to the water for food needs for growing crops, water is also essential for livestock production-beef cattle and dairy are an important source of food in NM. Globally, irrigated crops and raising animals consume ∼70% of the total freshwater. Livestock production needs water for animal watering, growing feeds, and on-farm needs such as cleaning, sanitation, cooling, and waste disposal systems; however, direct water consumption for drinking accounts only for 0.5% of the total livestock water footprint (Sawalhah et al., 2021). Most water requirement in livestock production is to grow feeds (95% of the total water footprint). In NM, 28,203 L of green (precipitation) and blue (surface and ground) water are required to produce one kilogram of beef, where blue water accounted for only 18% of the total water footprint (Sawalhah et al., 2021). In 2015, livestock production accounted for 1.1% from NM total freshwater withdrawals compared to 82.5% withdrawals by irrigation (USGS, 2015). The highest annual blue water footprint in NM's livestock production is associated with dairy cows. Marston et al. (2018) reported that dairy cows, beef cattle, sheep, goats, and hogs consumed around 23, 597; 20, 610; 247; 78; and 7 thousand m 3 per year, respectively. Additionally, more water is being used to produce beef than anything else humans eat, requiring 15, 415 liters of water (primarily to grow feed for animals).
Crops also affect natural flow regimes, yet these effects are not well-understood. Return flow from irrigation systems (nonconsumptive portion of water withdrawals) affects flow balance and water use accounting at a basin scale (Cai et al., 2003;Ochoa et al., 2020). Determining return flow, especially the utilizable return flow volume, is important for not only understanding water balance in streams and aquifers but also determining water availability for the development of more reasonable (and sustainable) water rights at the river basin scale (Grafton et al., 2012). For instance, the conversion of water-intensive crops to low-water crops (e.g., sorghum) has been identified as one of the effective methods to ease pressure on an increasingly limited water supply. Water managers face a continual challenge to meet the needs of multiple users.

NM FEW Nexus Drivers
FEW nexus components act as endogenous factors influencing one another (Chang et al., 2016). External (exogenous) drivers can simultaneously affect their behavior; thus, they need to be considered to address resilience challenges in terms of opportunities and trade-offs. External drivers can include climate variables and related extreme events-that is, temperature, precipitation, and drought-and socioeconomic stresses such as economic growth, population growth, poverty, political stability (Hameed et al., 2019). They can directly be linked in the FEW-WISE framework through multiple interactions and feedbacks across spatial and temporal scales (Chenoweth et al., 2011;King and Jaafar, 2015;Scanlon et al., 2017), thus creating risks and management challenges.

Temperature and Precipitation
NM annual precipitation showed a slight decreasing trend, while an increasing trend in mean annual temperature was observed for the period between 1895 and 2019 ( Figure 7A). This combination can have a profound negative impact on mountain snowpack that feeds water supply reservoirs by reducing water flow to the river basins and thus affecting water availability during the growing season. Even if snowpack accumulation was not to decrease, the projected higher temperatures will lead to an earlier initiation and end of snowmelt, potentially necessitating changes in water management.
Multiyear periods of high and low precipitation have resulted in very large swings in reservoir water supplies for agriculture. For example, the water levels in the Elephant Butte Reservoir were high from the 1920s to the 1940s before dropping to low levels during the mega drought of the 1950s until the 1980s. High levels remained throughout the 1980s and 1990s until falling again in the first part of the 21st century to about 10-15% of its capacity ( Figure 7B).
NM's precipitation and temperature trends ( Figure 7A) demonstrate that the deficits between potential evapotranspiration from freshwater bodies or vegetated surfaces and precipitation are increasing. This increased deficit can reduce streamflow and reservoir storage ( Figure 7B) and increase drought severity (Figure 8). With persisting drought conditions, farmers are increasingly dependent on pumping groundwater, rather than surface water sources (e.g., reservoirs and streams) to irrigate their crops and make up for the deficit that results in high production costs. Subsequently, the variable water supply threat has led farmers to plant more droughtresistant crops like beans, or to abandon their fields altogether. In NM, the most extensively pursued in terms of agriculture land area used is livestock grazing. Due to increased variability in precipitation, the amount of water is not sufficient for the optimum growth of forage on rangeland. Thus, this limited and variable rangeland productivity has resulted in reduced feed for livestock, and some ranchers adopt risk-averse strategies of selling their cattle to reduce the cost of buying additional feed supplements (e.g., hay) Gedefaw et al., 2021). At the same time, ranchers significantly reduce cattle herd sizes, allowing grasslands to recover from drought conditions (Uyttebrouck, 2013).

Drought Impacts
The southwestern United States is particularly vulnerable to drought, and even a small decrease in water availability in this already arid region can stress natural systems and further threaten water supplies. Drought poses a persistent risk to NM, added to its normally arid to semiarid climatic conditions. Drought events have broken historical records in recent years (NOAA, 2017). The extended record indicates that droughts were frequent in NM and more severe in recent years ( Figure 8) (NOAA, 2017;Johnson et al., 2020), and this trend is projected to continue (Schwalm et al., 2012).
Drought can directly and indirectly disrupt the state's most vulnerable economic activities including farming, ranching, and other sectors linked to agriculture. Recent extreme droughts have negatively impacted NM ecosystems such as the Chihuahuan Desert, causing grassland degradation and resulting in reduced grazing capacity for livestock (Gedefaw et al., 2021). Indirect impacts include, for example, reduced employment opportunities and increased agricultural input purchases. Drought also has indirect economic impacts on other sectors including public water supply, industry (e.g., mining, crude oil, and natural gas production), tourism and recreation, and tribal economies. Recurrence of a multiyear severe drought like that of the 1950s would have greater impacts on the food and water resources and the economy of the state than in the 1950s Frontiers in Environmental Science | www.frontiersin.org June 2021 | Volume 9 | Article 667018 because of warmer temperatures, population growth, and increased demand for water since the 1950s. Nevertheless, climate change will affect NM food systems and could undermine food security by reducing farmlands by 20-25% and irrigation water by more than 90% (USEPA, 1998;Voiland, 2013). NM experiences lack of funding for short-and long-term water planning, which can increase water sustainability risks compounded by the general lack of awareness of predicted future declines in water availability . NM is not alone as it is facing the same challenges experienced in the Western United States with increasing demand, reduced supply, and an inadequate legal and regulatory framework that was set up over 115 years ago when the conditions were different. There is a concern that the lack of funding and appropriate efforts will result in prolonged periods of extremely low water flow as warmer temperatures increase water losses due to increased evaporation. With increased reliance on groundwater, deeper groundwater wells would be needed to accommodate declining groundwater table, posing additional accessibility challenges and risks as groundwater often contains higher levels of salinity and other minerals that can reduce crop production and increase soil contamination (Frisvold et al., 2013).

Socioeconomics
NM FEW nexus functionality has direct and indirect impacts on the state's revenue and expenses, employment, and gross development product (GDP). Specifically, NM energy and food sectors have wide-ranging benefits to the state economy and New Mexican's livelihood and well-being. Together, agriculture and food processing support almost 50,000 jobs and contribute $10.6 billion to NM gross state product (GSP) in 2012 (Diemer et al., 2012). The agriculture and food processing industries directly created 32,578 jobs and 18,308 jobs in related support activities for a total of 50,886 jobs statewide (Diemer et al., 2012). These socioeconomic indicators have varied over the years based on several factors that include price fluctuations (crude oil and natural gas) and natural hazards such as drought. For example, the economic well-being of farmers and ranchers, which can be evaluated using net gains and losses in cash income, decreased by 36% between 2007 and 2012. Also, from 2002 to 2012, the number of farms in the United States with net losses increased by 0.2%, while that for NM increased by about 81.4% (New Mexico First and New Mexico State University, 2016).
The total value of the NM agriculture sector in 2018 was ∼$3.17 billion, while that of livestock production was $2.18 billion (NASS, 2018). Beef cattle and dairy are the most important agricultural activities in NM, contributing ∼42 and 56%, respectively, in cash receipts (NASS, 2016). Major portions of NM agricultural and processed food products (∼97%) and cattle (∼99%) are processed out of the state (NMDA, 2010). To meet NM consumers domestic demands, locally produced agricultural products of $13 million were sold directly to consumers in 2010. The prices of these food commodities coincided with higher energy consumption required for their production and processing.
Crop production requires significant direct (fuel and electricity) and indirect (fertilizers and pesticides) energy inputs that are mostly based on crude oil and natural gas. Price fluctuations of these energy sources directly impact agricultural production. Energy prices during 2011-2014 were 20 times higher than those of the 1970s and can be largely attributed to a combination of increasing global demand and increasing market fluctuations ( Figure 9A). On the other hand, a huge drop in crude oil prices from June 2014 to March 2016 negatively impacted NM's economy through declines in revenue that ultimately affected services provided to New Mexicans and resulted in a state fiscal crisis ( Figure 9B).
Regarding the energy sector, price fluctuations have a direct impact on the state's economic growth and stability. During 2012-2014, there was considerable economic growth that was attributed to increased crude oil and natural gas revenues primarily because of increased prices and production (EMNRD, 2015). Depending on the demand and prices, NM generally receives over $2 billion annually in direct revenue from crude oil and natural gas industry through severance taxes, property taxes, and royalty and rental income. Additional indirect income comes from sales and income taxes on crude oil and natural gas drilling and services, which can generate ∼$300 million. In NM, the revenue generated by the crude oil and natural gas industry can directly and indirectly impact New Mexicans in multiple ways as it contributes about 35% of funding   (Moskowitz, 2017;BEA, 2018) and by 3.7% in the year 2019 from that of 2018. On average, the total impact of agricultural production is about 7.4% of NM GDP and that of food processing is about 5.7%. These two broad industries accounted for $10.6 billion (around 12.3%) toward NM GDP by state of $86.5 billion in the year 2012. The agriculture sector provided 41,961 jobs in NM in 2012 (Diemer et al., 2012). Figure 10 shows the contribution of different industries to NM GDP 2019.
In 2017, the largest contributor (∼7% of $11.3 billion) to real GDP growth was natural gas and crude oil production (BEA, 2018). Mining, quarrying, and crude oil and natural gas extraction were the state's largest growth sectors, contributing about 2.2% to NM GDP and 1.7% to GDP change from 2018 to 2019, while this only contributed about 0.33% to the GDP of the entire United States. One-third of the annual revenues and contributions from crude oil and natural gas industry are used to provide services in NM. Energy-related high-paying jobs are estimated as 6.4% of employment in the state (EMNRD, 2019). Accounting for the impacts of FEW nexus on New Mexicans' well-being needs to consider the abovementioned factors to develop resilient social life before, during, and after climate extreme events and socioeconomic stresses.

System Dynamics Modeling
The SD approach is defined as "the study of the information-feedback characteristic of industrial activity to show how organization structure, amplification, and time delays interact to influence the success of the enterprises" (Forrester, 1958;Forrester, 1961). It has been used since the 1950s in diverse applications including economics, sociology, ecology, and engineering. It is based on the notion of system thinking (Forrester, 1961;Sterman, 2000) and allows to comprehensively review the structure and dynamics of complex systems (Tenza et al., 2017). It emerged as an innovative approach to facilitate holistic analysis of coupled human-environmental systems such as FEW systems (Kotir et al., 2016;Turner et al., 2016;Tenza et al., 2017;Xu and Szmerekovsky, 2017).
SD has the ability to integrate disparate systems such as FEW and HWB (Green et al., 2011;Qin et al., 2011;Shannak et al., 2018;Sušnik et al., 2018). It was therefore proposed to develop the FEW-WISE framework to identify the interactions between FEW's drivers at multiple spatial and temporal scales (Dietz et al., 2003). Specifically, development of SD models uses causal loop diagrams (CLD) to qualitatively represent systems' interactions. CLD can visually represent the links between the system variables using arrows and describe their effects. For example, a positive relationship indicates that an increase in a variable brings about an increase in the other (i.e., both variables change in the same direction). In contrast, a negative relationship indicates that an increase in a variable brings about a decrease in the other (i.e., the variables change in opposite directions). Another feature of SD is its ability to model feedbacks and delays in systems' response over time.
A generalized CLD for NM FEW nexus ( Figure 11) highlights multiple loops that describe how food production and prices; energy production and prices; and state revenue and severance tax incentives interact. Qualitative CLD can be translated to quantitative estimation of FEW nexus parameters. An SD simulation and prediction can offer a "virtual world" to analyze the influence of interconnected variables on FEW systems' behavior through scenario-based sensitivity analyses (Richardson, 1995;Ford, 1999). This SD approach will be further used to develop livestock, forage, dairy, energy, and water interconnected relationships for NM.

FEW Nexus Assessment Using Indices
As a part of FEW-WISE framework, the qualitative description of FEW systems' component linkages (NM FEW Nexus Description and NM FEW Nexus Drivers section) will inform the flows and connections of stocks and drivers in the SD model. The time series output of SD modeling will further be utilized to derive indices to monitor the FEW nexus, assess its adaptive capacity and resilience, assess potential resource management scenarios, and guide the decision-making process to develop effective policies to enhance resilience. The indices considered in this analysis represent and integrate two elements: FEW resource use (and availability) and HWB. The availability and use of FEW resources will be quantified based on the qualitative description of the interactions among different FEW system components (Food and Energy to Water and Food sections). Enhancing HWB can be achieved mostly through economic development but at the cost of degrading FEW resources (Foreman et al., 2003;Overpeck et al., 2005;Carpenter et al., 2006)-a process that acts as a barrier to resilience. A limited number of indices in nexus research have recently emerged (Willis et al., 2016) such as the Notre Dame Global Adaptation Index (Notre Dame Global Adaptation Index (ND-GAIN, 2015), the Environmental Performance Index (Wendling et al., 2020), the Human Insecurity Index (Werthes However, these indices focus on resource vulnerability, use, or accessibility but lacked in integrating HWB. For instance, the FEW Security Index is based on the principle that sufficient, widely accessible FEW systems are required for human development based on the fact that high security in one or two of these components cannot completely compensate for low security in others. This research focuses on developing an adapted index that couples FEW and HWB systems based on the principle that measuring only the aggregated amount of available resources to support human development is an insufficient measure of security and that the distribution of FEW resources is also an important consideration. The linkage between resources and human development can directly be determined by whether an adequate amount and good quality of resources are provided to meet the needs of growing population.

An Integrated FEW Index (Idx FEW )
The concept that will be followed in developing these indices is based on the availability, sustainability, and resiliency of FEW resource use to account for short-and long-term changes and shocks. To calculate these indices, it is important to describe the quality and quantity of FEW resources in terms of use, production, and availability through indicators. The qualitative description of the interlinkages (NM FEW Nexus Description section) can be used to guide the selection of the most important indicators to represent the FEW systems under investigation, expressing the linkages between variables using a process-based approach as exemplified in a CLD (Figure 11). This study, for example, considered the amount of food commodities produced (regardless of their respective nutritive values and ability to meet the demand) and required, in one way or the other, a form of resources (water and energy). For each subsystem, various indicators will be chosen to provide information about the behavior of the integrated system for which high-quality historical data are needed at the NM county and state levels ( Table 1). These indicators of FEW index will be contextdependent as they are developed considering the qualitative and quantitative description of FEW system components of NM. However, the calculation and derivation of this index are more generic that can be applied and used in other regions and other contexts like food security.
The FEW indices can be calculated as follows: where Idx Food is based on food production amounts for water and energy to the total food production (including imported) multiplied by the fraction of the total food production of a year to the long-term average of total food production. The total food production refers to the amount produced in a specific year, and the long-term average refers to the average production for the entire period of the data, for example, there is total production data for NM since 1950s. Thus, the food index can be directly (but partially) linked to ET (which can depend on climate variables-temperature and other variables), cropland area, number of animals, and animal forage consumption among others. The Idx Energy is based on the ratio of energy used for food and water to the total energy production multiplied by the FIGURE 11 | Casual loop diagram for NM's FEW systems. The arrows describe cause-and-effect behavior between variables. For example, an arrow from A to B indicates that A causes B. Pair-wise variable polarities represented as positive (+) (i.e., an increase or decrease in one factor causes an increase or decrease in the other factor) or (-), which is the opposite of a positive influence (i.e., an increase or decrease in one factor causes a decrease or increase in the other).
fraction of the total energy production (from oil and gas) to the long-term average of the total energy production from all sources. Idx Water is water used for food and energy production relative to total water withdrawal multiplied by the fraction of the total withdrawal to the long-term average withdrawal. For example, the amount of water consumed by crops can be estimated from evapotranspiration (ET) and the area of cropland. The amount of water consumed by livestock can be based on ET of natural vegetation and animal feed consumption, animal numbers, and other variables. These two water estimates can be estimated based on physical and/or empirical (statistical) models. The total amount of water used to produce food commodities is considered in the calculation of water index (Idx Water ). The above three indices will be combined using the weighted average to derive a single indicator explaining the status of FEW systems as follows: where i is the number of FEW indicators. The selection of an appropriate time scale based on the indicators is required to develop a consistent time series of FEW indices. The availability of such time series is key in monitoring and evaluating the behavior FEW nexus in response to the abovementioned stresses. In other words, these indices can be compared with drought indices and energy prices.

Human Wellbeing Index (Idx HWB )
Humans benefit from their interaction with the Earth's natural resources base as it provides environmental, economic, and social capital to enhance their well-being. HWB is a broad concept that cannot be observed directly from surveys or measured independent of various social factors. It could possibly have multiple domains, and social scientists have developed broad categories to draw general distinctions between them. Within each domain, there is a set of subcategories (or indicators) that identify specific components of HWB. There is no general agreement on a set of domains and indicators to describe HWB. In the literature, seven domains have been identified to be broad enough to encompass most research frameworks such as (Hagerty et al., 2001) relationships with family and friends, emotional well-being, material well-being, health, work and productivity, feeling part of one's community, and personal safety (Cummins et al., 1994;Cummins, 1996). The list of potential indicators is even longer, and no comprehensive list exists (e.g., education, employment, energy, population, and leisure activities.) (Boelhouwer and Stoop, 1999;Diener et al., 1999;Marks et al., 2006;Costanza et al., 2007). These domains and indicators are used to understand and categorize the concepts, status, and trends of HWB based on a set of indices. One such index is the Human Wellbeing Index (HWBI) of the U.S. Environmental Protection Agency that characterizes HWB at multiple scales (e.g., national, regional, or local) (Summers et al., 2017). The HWBI can use a substantial group of indicators to evaluate the influence of social, economic, and environmental domains in an integrated fashion based on well-being, applicable to communities at multiple scales. In this study, a subset of social and economic domains defined by Summers et al. (2017) will be modified to estimate the HWB index for NM (i.e., Idx HWB ) ( Table 2).
Using the following equation, three indicators associated with the economic domain of HWB will be combined to calculate Idx HWB according to the following equation: where and i w HWB i 1. where i is the number of HWB indicators. The indicators (e.g., adults working long hours spending less leisure time, degradation of natural FEW systems, and educated population associated with specific age-groups) associated with the social domain of the HWB system will be defined objectively to determine the time series HWB index. A qualitative nonmeasurable criterion or scores from low to high will be given due to their related effectiveness to the quality of life. The selected economic and social domains of HWB are reflective of NM economy and social conditions due to contributions of energy (e.g., crude oil and natural gas) and food (e.g., livestock)

System
Measures/Indicators Food Amount of food products (metric tons) used to generate energy (e.g., biofuels) (Btu); food imports over total exports (%); total food production (metric tons) Energy Electricity (KWh) and crude oil and natural gas consumption (Btu) needed to produce food (metric tons) and water (m 3 ); total energy production (Btu) Water Total water withdrawals (m 3 ), water for human consumption (m 3 ), water withdrawal to produce food and energy (liters per kg/KWh) Frontiers in Environmental Science | www.frontiersin.org June 2021 | Volume 9 | Article 667018 industries toward revenues, degradation of FEW, or natural resources (Gedefaw et al., 2020) that contributes indirectly to the well-being of the human society by promoting and attracting tourism-based economic and social welfare. Both FEW and HWB indices are generic, and the weights are dimensionless, varying from 0 -1, but they need to be evaluated regionally and contextually. The units of different indicators of the two indices (i.e., Idx FEW and Idx HWB ) are described in Table 1 and Table 2, respectively. The approach of estimating indices would allow the FEW-WISE framework to be adaptable and replicable in another contexts.

Resilience
NM FEW nexus is strained by some of the most challenging risks toward achieving resilience. To address these risks, resource managers heavily rely on these socioecological systems' capacities to absorb disturbances while maintaining essential functions (Holling, 1973;Folke et al., 2002). But when a critical threshold is surpassed, a system under stress can undergo catastrophic changes and reorganize into a different state (Angeler and Allen, 2016). An improved understanding of the boundaries of a system resilience (i.e., thresholds) that separate one state from another alternative potentially undesirable state can help in developing adequate management practices.

Concept
The concept of resilience was introduced more than 40 years ago in ecological sciences and was defined as the amount of disturbance that a system can withstand before shifting into an alternative stable state (Holling, 1973). The concept has been used in human development (Brown and Westaway, 2011) and socioecological systems (Folke, 2006;Curtin and Parker, 2014;Desjardins et al., 2015) among others. It emphasizes that socioecological systems need to be managed and governed for flexibility and emergence, rather than for maintaining stability (Peterson et al., 2003;Carpenter et al., 2015). Hence, resilience is concerned with navigating complexity, uncertainty, and changes across scales (Berkes et al., 2003;Cash et al., 2006;Cumming et al., 2013) on a human-dominated planet (Lubchenco, 1998;Steffen et al., 2007) by combining the concepts of adaptability and transformation. Adaptation is a process of deliberate change in anticipation of or in reaction to external stimuli and stress (Nelson et al., 2007). Transformability is shifting development into new pathways and even creating novel ones.
Given this concept, the FEW-WISE framework adopted the ecological resilience concept (Holling and Meffe, 1996;Quinlan et al., 2016), which refers to the magnitude of disturbance (i.e., distance to critical transition) that a system can absorb before shifting to an alternate system state (Holling and Meffe, 1996). It assumes that a system has multiple alternate equilibria and the capacity to maintain its essential structure and functions through reorganization-a key property of complex adaptive systems. Disturbances slow down system's processes through abrupt, gradual turns, or bifurcate toward an alternate state. These turning points are referred to as tipping points as the system passes critical thresholds (Dakos et al., 2015). A system becomes unstable once reaching such thresholds but afterward should attain a new stable, adaptive, and transformed state (Kuehn, 2011). A dynamic behavior of a system consisting of two different states, A and B, before and after a disturbance, respectively, is depicted in Figure 12. The recovery rate (resilience) of a system is high when the distance to a critical transition (h) is higher and low when the distance to critical threshold (ds) is shorter. An FEW nexus can reorganize (Ullah et al., 2015) and attain multiple resilience states, while undergoing abrupt shifts (Scheffer and Carpenter, 2003;Scheffer et al., 2009) through human interventions that can introduce different adaptive measures (Folke et al., 2004;Biggs et al., 2012;Schoon and Cox, 2012). The concept described in Figure 12 can be used to assess NM FEW systems' resilience by quantitatively characterizing and predicting their response to the current and future climatic and socioeconomic changes.
Resilience as a system property should not be reduced to a simple index, but different types of indicators need to be combined to capture facets of resilience. A Resilience Index (RI) needs to reflect the complex adaptive behavior of socioecological systems-combining biophysical (e.g., FEW) feedbacks and trade-offs with socioeconomic (HWB) drivers and outcomes (Givens et al., 2018).

Modeling
SD can be the most appropriate approach to model resilience with its structured framework that allows simulation of simulate complex feedbacks of key resource variables (Luna-Reyes and Andersen, 2003;Tenza et al., 2017), can integrate numerous interactions into a set of nonlinear expressions similar to those governing coupled FEW-HWB systems, and most importantly, focuses on the evolution of a process of interest, rather than achieving a specific equilibrium or optimal solution. Particularly, what is more pertaining to resilience is that SD models are essentially systems of equations with no "closed-form" solution, where comparison of the evolution of different scenarios is the primary mode of analysis. A critical aspect of resilience is to monitor the evolution of processes and feedback loops. The analysis of CLD can provide insights into the root causes of a system behavior; help in identifying changes in variables, processes, and drivers; and indicate threshold of a state variable that is prone to be overpassed. For example, the occurrence of trade-offs, synergies, and cross-sector dependencies in a system can indicate new technology (oil and gas hydraulic fracturing), policy, or environmental changes that transform it into a new state, demonstrating the adaptive management of FEW resources. Therefore, an RI needs to couple multiple indices (e.g., Idx HWB and Idx FEW ) (Heckbert et al., 2014).

Indices
The FEW-WISE framework ( Figure 2) involves resilience assessment based on an RI (Eq. 9) that integrates and evaluates coupled FEW-HWB systems. The proposed RI (Idx Reslience ) will be developed to monitor changes in the systems at the county level based on a criterion related to potential risks above and below a threshold value as moderate, high, and extreme. The objective of the RI is to determine which NM counties are most and least vulnerable to climate and socioeconomic risks. The Idx Reslience integrated four indices: Idx Food , Idx Energy , Idx Water , and Idx HWB .
where the weights w RI 1 , w RI 2 , w RI 3 , and w RI 4 represent the contribution of each index. The weights, which are dimensionless, vary from 0 to 1 and sum to 1. They can vary according to the focus of the specific FEW systems. If all four elements are equally averaged, then weight equals 0.25. More representative weights will be developed based on questionnaire and survey analysis with different stakeholders in NM. The RI will be considered as a risk indicator that ranges from 0 (lowest resilience) to 1 (highest resilience).
Based on Idx Reslience , the counties with low resilience will be identified and prioritized for conducting a scenario-based analysis to develop alternate management strategies to enhance their resilience (Figure 2), allowing development of potentially effective resource management options along with their trade-offs and synergies (e.g., improving energy plans, changing the livestock and crop production portfolio, and identifying alternative water sources).

CONCLUDING REMARKS
The goal of this analysis was to develop the FEW-WISE framework for an improved assessment of FEW systems resilience in response to environmental and socioeconomic shocks. The state of New Mexico was used as a case study to develop this framework in terms of its unique characteristics and challenges in terms of climate change impacts that affect its FEW resources and New Mexicans' well-being. To achieve this goal, FEW-WISE followed five steps: 1) qualitatively identify interconnections between the FEW system components and quantitatively assess resource exchanges, 2) identify important systems drivers, 3) integrate FEW systems using system dynamics (SD) models to evaluate their response to the identified drivers, 4) develop FEW-HWB performance indices, and 5) develop a resilience monitoring criterion using a threshold-based approach that integrates these indices. The first two steps have been described in detail here. The other three steps were individually conceptualized and will be carried out and presented in separate modeling, indices, and resilience assessments. With the application of the framework, the expected results will include time series of predicted values of FEW components, FEW-HWB, and resilience indices ranging from 0 to 1.
The SD modeling is a key feature of the framework as it simulates the nexus as a holistic multi-sectoral system, providing insights into the vulnerability of resources to stresses, demonstrating how the nexus will respond to changes and transition to absorb the effects from these stresses. The SD modeling along with the resilience indices potentially provides the needed linkages to the decision-making processes. It should be noted that the HWB and FEW indices can be developed based on a wide range of variables or indicators. However, data availability challenges exist, especially those related to the links FIGURE 12 | (A) Stages in the performance of FEWS toward resilience in response to disturbances with (A) tipping points t1 and t2 and critical thresholds f1 and f2 which are overcome by changing the state of a system in response to drivers, and (B) and (C) show the systems with low and high resilience and recovery rate, respectively.
to food resources in water and energy systems, and water used for crop production, and human well-being. While the use of variables with most available data is important, some variables may not be the best ones to represent the systems or regions under consideration.
In conclusion, the FEW-WISE framework with its indices can be used to develop alternative scenario-based management strategies for FEW resources to enhance resilience and sustainability of ecologically and socioeconomically vulnerable regions. As NM has FEW systems that are sensitive to drought and fluctuations in energy prices, this framework would allow the state to effectively manage its resources. The framework is currently proposed at county and state scales over a semiarid region and context; however, its methodology can effectively be adapted, redefined, and transferred to operate over different regions and contexts in other parts of the world as needed.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
KY and HG conceived the idea for this study. The FEW nexus components and drivers were described by KY and HG. The indices for human well-being and resilience were developed by HG. MS, LJ, and AZ, and MG provided datasets related to FEW components. Tables and figures were generated by KY along with the first draft of the writing. The final version of the manuscript was edited by all the authors. The final edits were compiled by KY who converted the manuscript to the final format for this journal.