Simulating a Watershed-Scale Strategy to Mitigate Drought, Flooding, and Sediment Transport in Drylands

: Drylands today are facing a landscape-scale water storage problem. Throughout the increasingly arid Southwest of the United States, vegetation loss in upland watersheds is leading to floods that scour soils and transport sediment that clogs downstream riparian areas and agricultural infrastructure. The resulting higher flow energies and diminished capacity to infiltrate flood flows are depleting soil water storage across the landscape, negatively impacting agriculture and ecosystems. Land and water managers face challenges to reverse the trends due to the complex interacting social and biogeophysical root causes. Presented here is an integrative system dynamics model that simulates innovative and transformative management scenarios. These scenarios include the natural and hydro-social processes and feedback dynamics critical for achieving long-term mitigation of droughts, flooding, and sediment transport. This model is a component of the Flood Flow Connectivity to the Landscape framework, which integrates spatial and hydrologic process models. Scenarios of support and collaboration for land management innovations are simulated to connect flood flow to the floodplains throughout the watershed to replenish soil storage and shallow groundwater aquifers across regional scales. The results reveal the management policy levers and trade-off balances critical for restoring management and water storage capacity to the system for long-term resilience.


Introduction
Throughout the American Southwest and many global dryland regions, drought, increased flooding, and sediment transport are exacerbating landscape-scale water challenges [1]. Less winter precipitation and higher temperatures are diminishing snowpack storage [2,3], which results in less spring runoff quantities and flow durations for downstream water users [4]. Landscape losses of deep soils and vegetative cover in this region stem from historic overgrazing and fire suppression corresponding with climate effects that began in the late nineteenth century [5][6][7]. Rain and snow melt are no longer held by the soils and released slowly throughout the year; instead, they run off immediately in floods, resulting in catastrophic flows and severe erosion [5]. As vegetation on uplands decreases, hydrologic energy in valley bottoms increases, which increases scouring of floodplain vegetation and soil along the channel bottoms, entrenching flows in the valleys. Historically, floods along the many river networks were more connected to more richly vegetated floodplains [7]. When flood flow is connected to and can overbank onto floodplains, watersheds retain more their water and soil resources and result in multiple ecosystem services such as aquifer recharge, flood mitigation, vegetation productivity, and water quality treatment [8]. Agriculture traditionally and can again support this natural dynamic through watershed-scale systems of stormwater harvesting, with a general approach to spread flood flow onto floodplains employing practices such as flood irrigation. This management approach provides critical functions for long-term resilience in that it maintains the buffering capacity to drought, flooding, and erosion challenges [9] and results in infiltration into shallow groundwater aquifers [10].

Connectivity Processes Drive Dryland Landscape Dynamics
To achieve transformation of the landscapes and the communities that rely upon them to a resilient state requires collaborations of integrated management systems and disciplines across the scales [11]. To achieve such holistic land management requires developing common understandings of the system dynamics, the functional as well as the structural attributes of a system, which is often described as restoring the natural or biophysical functions and processes [12,13], or a process approach [14]. A resilient system can then be characterized and assessed by the processes and feedback loops that propel a system across a critical threshold to a system state that achieves desired functional goals [12,15]. Despite a general recognition of the importance of function throughout natural system disciplines, most approaches still assess natural systems from only a structural perspective, and many disciplines face challenges in taking a process approach due to the difficulties in conceptualizing and measuring certain processes [12,16]. An additional obstacle to a holistic system understanding is that the boundary of the system that reveals the dynamics is often divided; for example, upper watershed managers are not commonly in collaboration with downstream community efforts to address flooding and water scarcity issues. Current approaches to flood risk mitigation commonly do not often look to upper watershed management as part of the solution. "Mono-functional" flood control infrastructure (FCI) approaches often limit considerations of the surrounding natural system [17,18]. The goal is typically to control flooding through restricting area for river flow, such as using levees to protect floodplains under use. These spatial restrictions, however, decrease the buffering capacity of flooding processes and the overall system resiliency, which increases our vulnerability to catastrophic failures [17,19,20]. As increased flooding exposes these vulnerabilities, calls are emerging for process approaches on watershed scales, often termed eco-engineering approaches or nature based solutions [21,22], to adapt the socio-natural systems to flooding processes within the flood control disciplines [23]. A few programs under development have adopted these approaches, such as the "Room for the River" in the Netherlands [21,[24][25][26].
Connectivity in watershed systems is emerging as a process that organizes dryland landscapes and is thus useful for revealing dryland dynamics [15,18,[27][28][29][30][31][32][33]. Practitioners that developed connectivity as a tool for dryland watershed restoration placed the relationship between structure and function as a central perspective [15,18]. The dryland connectivity conceptual framework defines connectivity as the extent to which materials such as water, nutrients, and organisms can move, spread, or be redistributed from one place to another within the landscape [15,34]. Quantifying the connectivity and disconnectivity attributes can determine the condition of the watersheds and characterize the flux dynamics [30,31].
A critical driver of watershed conditions are the connectivity conditions of valley bottoms where water and other flows are concentrated throughout the watershed. When hydrologists refer to hydrologic connectivity, they are often referring to longitudinal connectivity, i.e., how fast flow proceeds to the watershed outlet [30]. An overall restoration strategy is to increase lateral connectivity, whereby channel flows are connected and able to flow laterally onto the adjacent floodplains [35]. This reconnection of the flow to the floodplain can address the extensive problem of increased hydrologic energy with scouring peaks of flow, as shown in Figure 1. More of the watershed would serve as a runoff sink, and less as a runoff source, with more water retained and overall storage increased [36]. An increase in lateral connectivity throughout a watershed increases the area of runoff sinks by spreading and slowing the flow, recharging soil storage capacity, and supporting vegetation productivity. Vegetation then provides the feedback function of increasing infiltration and water holding storage capacity through increasing surface roughness, introducing soil pathways, and increasing organic matter into the soil [37]. This dynamic is well-understood in perennial river ecosystems [38], but often overlooked in the upland (often managed grazing) systems, which contain the intermittent headwater valley bottoms (here, intermittent includes all temporary waterways, both intermittent and ephemeral). Intermittent systems, however, have unique ecosystem dynamics and require distinct management, monitoring, and research approaches [39,40]. These critical conduits of water, energy, material, and organisms support high biodiversity and important ecosystem processes [41]. A valley bottom system of well-connected channels to floodplains throughout the landscape, including the intermittent systems, serves as the buffer against both flooding and the effects of growing aridity [41].
For long-term resiliency, the natural system requires management of the connectivity of flood flow to the floodplains throughout the landscape to mitigate hydrologic energy. The slowing of hydrologic energy has the potential to result in replenishing soil storage and groundwater systems, which provides critical buffers against disturbances. This was the deciding factor for the boundary of this work. Other practices could contribute to the overall objectives, such as aquifer storage and recovery (injection wells) increasing groundwater storage [42]). However, targeting the existing management systems and incorporating floodplain reconnection restoration practices into daily activities of land managers across the landscape has the potential to have the furthest reach [43]. Floodplain reconnection is still not widespread, but it gained global recognition for its ability to control floods through restoring floodplain functions that yield multifunctional social and ecological benefits [25,[44][45][46][47][48][49]. Many traditional and innovative land management practices can address these natural processes and functional goals of connecting flood flows to floodplains, with different practices appropriate for varying local conditions and dynamics, particularly for differing scales and hydrologic energy conditions. These include flood or spate irrigation [9,50,51], terraced agriculture [52], riparian buffers that filter and slow flow [53,54], and porous large wood and other debris and rock structures in both perennial and intermittent flow systems that slow but do not impound flow, mimicking beaver structure functionality [30,[55][56][57][58][59][60][61][62]. All of these ultimately modify the connectivity on the landscape to increase vegetation and capture litter and nutrients, with the term "conmods" [15] coined to describe them.

A Systems and Action Research Approach
To understand the processes and dynamics in complex systems, systems science approaches fulfill a crucial function by characterizing the key drivers and relationship dynamics that determine system states over time. System dynamics (SD) modeling approaches offer unique abilities to characterize the effects of key feedback loops and quantify their magnitudes in simulation models, and to assess potential alternative future scenarios [63][64][65][66]. Many hydrologists recognized that social system feedbacks are essential for assessment of hydrologic systems [63,[66][67][68], and hydrologists and other scientists in socio-environmental disciplines of socio-hydrology and hydro-sociology took up system dynamics approaches to understand the socio-natural system interactions, as "linear causal thinking cannot address complex challenges adequately" [69]. Hydrologists also took up SD modeling approaches to develop integrated water resources and watershed management approaches [69][70][71][72].
A systems approach enables revealing underlying dynamics that are fundamental to many systems, as well as generalizable innovative alternative management and scenarios that can transform them. Research in systems thinking identified that a core or generic set of "archetypes" can be found in many systems which are classifications of structures responsible for generic patterns of behavior over time [73,74]. Wolstenhome defines the most generic or basic archetypes as "underachievement", "out of control", "relative achievement", and "relative control" [73]. A generic structure such as the model described in this work typically results in dynamics that correspond to several archetypes over the course of the model run. While it is outside of the scope of this article to map each dynamic's corresponding generic archetype, this work did aim to explicitly address a widespread land management challenge that is one manifestation of the "tragedy of the commons" archetype. This is a commonly identified multiple-dynamic archetype where collective management tends a common pool resource system toward collapse. Each individual actor has an incentive to outcompete others by increasing his or her extraction of resources, leading to the reduction of resources below the levels of what is sustainable for the group as a whole [75]. In the system addressed by this work, interconnected resources across the landscape are being depleted. Management such as overgrazing in combination with increased occurrences of droughts resulted in over-extraction of the upland vegetation system. Addressing this root cause through spreading flood flow onto floodplains to support increased vegetation along upper watershed valley bottoms has cascading resulting effects of decreased flooding, increased infiltration, and an increase in an alternative water supply. This approach protects resources, infrastructure, and can help address surface water availability challenges that have increased groundwater reliance, a resource under the threat of collapse. Potential solutions require or can be more effective with a response archetype that fits into the category of "seeing and acting holistically" [73].
This management and policy approach is informed by collaborative action research and socioecological science (SES), approaches that have identified underlying principles that can thwart the tragedy of the commons. Current efforts to further understandings of "usable" science was long the terrain of the discipline of "action research", where the scientific efforts are defined by communities' actual problems and the need for collaborative efforts to address them [76][77][78]. A central tenet of action research that can be facilitated by SD modeling is that collaborative experimentation is critical to achieve any sustainable socioenvironmental system transformation [43,77,[79][80][81][82][83]. Enacting pilot studies and experiments can be conducted as scenario testing with a broad range of stakeholders and interdisciplinary scientists. These experiments both build the social capacity to innovate in the face of complex challenges and create models that enable society to visualize and integrate alternative management into daily life [77,83]. Scenario-testing also accelerates the refining of strategies that can be tested on the ground. A sociologist posited in the 1940s that, because these are complex systems, "if you want to know how things really work, just try to change them" [84]. As Ostrom showed, the tragedy of the commons and the incentives for individuals to compete in the depletion of common resources to catastrophic result is a dominant driver across societies, but not an inevitable one [83].

Main Aim of Work
The occurrences of droughts and a general increase in aridity are drying upland soils, decreasing vegetation, and resulting in scouring floods across the American Southwest. The complexity and scale of ecosystems and management systems impede efforts to distill process-based solutions. The main aim of this work was to create a generalizable and adaptable integrative socio-environmental system dynamics model to facilitate development of regional collaborative approaches to define thresholds for resilient states. This model simulates an innovative watershed-scale stormwater harvesting and aquifer recharge strategy to reveal the extent of management (connectivity between flood flow and floodplains) required to reverse degradation trends. Key to this approach is simulating the water budget implications of alternate management strategies, which underlies the structure of the SD model ( Figure 2). The scenarios were developed and iteratively modified to test the potential of achieving the following objectives: (i) the overall balance of system to manage flood flow/stormwater to increase downstream benefit while not negatively impacting downstream users, and (ii) to increase productivity/water availability to agriculture and reverse trends of groundwater storage declines. Managing the water budget targets the processes that connect water inputs into the storage functions and, thus, reveals the mechanisms for increasing water availability.

Integration of Other Elements of FlowCon Framework
The SD model addressed in this article is integrated with the outputs of other Flood Flow Connectivity to the Landscape (FlowCon) framework models. These include a spatially explicit landscape indicator model and semi-distributed hydrologic models. Aggregated FlowCon outputs that are used as inputs to this SD model are shown in Figure 3a,b, and more details are provided in Section 2.1.3. The initial FlowCon framework models identify optimum locations and quantify the resulting benefits and extent of management and collaborative support required for restoration of the critical landscape processes. These target processes and functional goals that reduce hydrologic energy through increasing water retention, recharge, and vegetative productivity. Figure 3a shows the locations identified for reconnecting floodplains at varying levels of priority, and Figure 3b shows the effects on a synthetic hydrograph of executing the various levels. The top two priority levels shown were found to be the optimum levels for reversing degradation trends. The FlowCon framework provides a tool for use across disciplines to identify, quantify, and locate the critical process intervention targets for employing practices which then can be fit to local conditions.
(a) (b) Figure 3. Outputs of the Flood Flow Connectivity to the Landscape (FlowCon) modeling framework that contributed to the generalizable and adaptable integrative socio-environmental system dynamics (SD) model addressed in this article. FlowCon identifies optimum locations and quantifies the resulting benefits and extent of management and collaborative support required for restoration of the critical landscape processes of reducing hydrologic energy through increasing water retention, recharge, and vegetative productivity. (a) Locations for reconnecting floodplains at varying priority levels; (b) effects on a synthetic hydrograph of executing the various levels. The top two priority levels were found to be the optimum levels for reversing degradation trends.

Materials and Methods
Presented here is the system dynamics model of the FlowCon framework, a generalizable and adaptable model structure that includes the natural and management processes and feedback loops critical for achieving transformation of watersheds, mitigating droughts, flooding, and sediment transport over the long term. Note that stormwater is usually associated with urban areas, but here the strategy is to spread and slow flood flows across the landscape; thus, in upper watershed areas, we use the term flood flow, and, in the valley areas with more developed infrastructure, we use the term stormwater.  Figure 4). Using this real-world example tests the model's outputs against our interpretation of history and the real system's behavior. This generic structure can then be adapted to specific management decisions specific to a region under question and, thus, serves as a starting point for working with communities to experiment on potential management innovations. Hatch and Mesilla Valleys, the largest areas of agricultural valleys along the New Mexican Rio Grande, rely upon snowpack to fill the reservoirs, which is becoming increasingly variable in its supply. The occurrences of droughts and a general increase in aridity are drying upland soils, decreasing vegetation, and resulting in scouring floods. Agriculture has the potential to manage these landscapes by slowing and spreading flood flow, which can recharge soil stocks and refill aquifers.

Model Water Budget Approach
The essential structure of the model uses a water budget approach (Figure 3). Critical system conditions which can be assessed through a water budget have important effects on a community's ability to mitigate drought, flooding, and sediment transport. The hydrology of the system is a major driver, both in terms of water quantity and distribution. A water budget as a management tool provides an estimate of the balance and distribution of water as it cycles through the system, and is most generally defined and calculated as the inputs minus the outputs are equal to the change in storage [86,87]. The budget estimates in this work allow for estimates of the effects on the water distribution and quantities from system innovations of interventions and the alternative land management simulated.
Hydrologists have utilized SD modeling approaches in conjunction with water budgeting to achieve greater system understandings. The New Mexico Water Resources Research Institute developed a Dynamic Statewide Water Budget system dynamics water budget model which produced outputs that this research relied upon for parameter ranges [88]. Modeling approaches that incorporate water budgeting allow for robust assessment of the competition between sectors and population clusters; for example, without policies that concretely identify agricultural needs, agriculture can be predicted to decline [89], an experience borne out over time [90][91][92]. Specific hydrologic processes modeled in system dynamics, including runoff, sediment transport, aquifer recharge, and aquifer storage and recovery, were shown to produce comparable results to other commonly used hydrologic models [42,93,94]). A critical link is to integrate with other modeling approaches, such as suites of models that include spatial modeling as this framework employs to facilitate the development of alternative watershed management strategies [94][95][96][97]. Integrating and linking spatially explicit and process-based models with integrative system dynamics frameworks was identified as critical to characterizing reference states and facilitate management planning [98]. Flooding risks motivated modeling approaches that explicitly address floodplain management as addressed in this work, although none to date included upland intermittent floodplain system interventions [95,99]. As Zischg [95] stated, floodplains are complex adaptive systems that are composed of co-evolving natural and human systems whose internal dynamics spawn emergent behavior that influences future risk pathways. These behaviors and management strategies include trade-offs between upstream and downstream users and the collective shouldering of risk. Baldassarre et al. [99] showed that, while human-water systems can be described as coupled, the two sides are not discretely distinct, and the attempt to describe the relationship between the systems is key to yielding innovative findings; in their example, the relationship between development aspirations in floodplains and the risks of flooding damage was described.

Model Boundary, Assumptions, and Inputs from Other FlowCon Framework Models
The boundary of the model is determined by the purpose that it serves, which in this work includes the extent of the system that reveals high-leverage management options for achieving long-term mitigation of droughts, flooding, and sediment transport. Table 1 summarizes the model boundary by distinguishing endogenous (calculated and affected by other variables within the model) and exogenous variables (inputs not changed by other model variables). Included as well is a list of variables that could have potentially been in the model but are excluded due to them being outside of the scope and focus of this research. The parameters used in this model were the scenarios of extent of floodplain connectivity that would result in mitigation of hydrologic energy to achieve the multifunctional ecosystem service benefits, increasing water availability, reducing sediment transport, and increasing productivity. This model enables projections of management of four main conditions: (i) the effect of water availability on agriculture that can increase or protect productivity, (ii) the balance of surface water between the upland and valleys below, (iii) the spreading of flood flow which enables the reduction of groundwater pumping (in-lieu recharge) [100], and (iv) the extent to which the interventions affect downstream users. Productivity benefits (PBu, PBv) (see Figure 14 for the graphical functions) Withdrawals change (Wc) Recovery policy ratio (Rp) Outputs from other FlowCon framework models were used as input parameters to the model (the major parameters are noted in Table 1). Utilizing other FlowCon framework analyses allows this aggregated SD model to benefit from analysis that is often best conducted spatially across the scale of the study boundary. The linking of spatial and process-based models is an important expanding area of work due to the unique qualities and management needs of particular areas of the landscape, such as the differences in infiltration rates of uplands, upland floodplains, and wide farming valleys [98]). A significant feature of the FlowCon framework is the delineation of the upland floodplains (using Topographic Wetness Index (TWI) values correlated to valley bottom and channel contours [101,102]) and the subsequent spatial and hydrologic modeling that characterized and estimated the dynamics of the flow as it courses through the upland valley bottoms and the benefits of reconnecting that flow to the floodplains. The floodplain reconnection benefit analysis served as the intervention targets for the tested scenarios for this model.
The overall FlowCon framework includes this SD model, which integrates data from the suite of the three other models. The first is a spatially explicit GIS-based landscape and hydrologic indicator model using remote sensing inputs (fine-scale 4.5 meter Digital Elevation Models (DEMs) [103] and observed radar precipitation inputs (NEXRAD) [104]), to produce a newly created expected wetness index (EWI). The EWI used here synthesizes the effect of precipitation intensity and antecedent wetness conditions, but has the potential and is currently being investigated to include vegetation density and pattern indicators using multi-spectral inputs [105]. The second model is a semi-distributed hydrological model using the SCS curve number (CN) method [106] in the HEC-HMS package, a commonly used tool by flood mitigation engineers [107]. The inputs of this model (CNs) are modified by the EWI inputs and calibrated by gauge data to quantify the locations, magnitudes, and frequencies of flow, and used to create a synthetic unit hydrograph from a series of actual rain events. The third model is a newly created hydraulic channel and floodplain network routing model in Excel. This model uses the hydrologic model results to distribute the runoff onto a finer scale, estimate the flow mitigation benefits from managing connectivity in optimum locations, and produce the spatially explicit FlowCon data. Several scenarios are tested to determine the extent of intervention required to address the functional goals. The results are then projected back onto the GIS-based landscape and hydrologic indicator model maps for use by the land managers.

Calibration, Validation, and Confidence Building Tests
The model successfully passed almost all conventional confidence building tests, including unit consistency, integration error, boundary adequacy, structural assessment, parameter assessment, behavioral anomaly, parameter sensitivity, and extreme conditions. On the topic of confidence building tests in system dynamics, refer to Sterman (Chapter 21) [108] for an exhaustive list, Barlas [109] for a methodological discussion, and Langarudi and Radzicki [110] for a practical example and guidelines. Unit consistency and integration error tests are straightforward procedures performed through Stella's built-in features. The boundary adequacy test was conducted by examining the impact of each additional structure on the model's behavior. The results indicated that all the endogenous variables were significant for the purpose of the study. Structural and parameter assessments were iterative processes which actually comprised the main model development phase. The majority of the logical structure of the model is explained in Section 1.1. The feedback structures and the parameters feeding them consist of local knowledge, literature, and technical estimations as described above throughout Section 2.1. At each stage of development, the model was tested for its behavior. It is normal that some anomalous behaviors appear at early steps, but they are fixed by modifying the structure. The same procedure was followed for parameter sensitivity and extreme condition tests which helped to identify a few flaws in the model structure. As an example, groundwater availability (GWa) and stormwater (SW) availability (Qa) were added to the model as first-order controls to prevent the water stocks from going below zero in extreme situations. These tests (especially parameter sensitivity) also helped to spot high-leverage points of the system, as discussed in detail in Section 4.
The model was also tested partially for policy sensitivity, which is further discussed in Sections 3 and 4. The only tests that the model was not exposed to were family member and system improvement. Family member tests examine generalizability of the model by applying it to real-world problems that are similar to the original study case in terms of structure and characteristics. To pass the test, a recalibrated version of the model should be able to mimic those problems with reasonable accuracy. The system improvement test examines the robustness of the interventions proposed by the model in the real world. In this regard, the recommendations should actually be implemented in the real system, which is obviously yet to happen. Considering that the model is at its earliest stage of development, none of these tests could be realistically performed. They will, nonetheless, remain a priority for future research.
It is worth mentioning that the behavior reproduction test was not performed on the model in a conventional sense mainly due to the fact that historical data for the major variables were unavailable, such as evapotranspiration. The advantage to a water budget approach is that one can use several variables as closure terms, where the variables are calibrated to the known variable values. The advantage to the system dynamics model is that the inclusion of feedback gives added confidence, and these calibrations are often finely tuned and sensitive. Another reason for not pursuing a data-fitting exercise was that exogenous factors that might be significant but out of our endogenous scope were omitted. This reduces the model's precision for point-by-point prediction but increases our ability to focus on key internal drivers. That is indeed why prominent system dynamics experts do not consider data-fitting as a strong validation test for SD models [111][112][113][114]. Therefore, the validity of our model relies heavily on the accuracy of its own internal feedback structure and data inputs that set its initial conditions. If the internal structure is deemed sound, then higher confidence can be placed in the outputs. In fact, the model's outputs were reviewed by local domain experts and their general plausibility was approved.

Upland System to Support Ranchers to Increase Vegetation
The SD model includes both the existing system and the scenarios of spreading flood flow and storm water, which are indicated in the model by variable shapes filled with dark-orange and bold connectors, as can be seen in Figures 5 and 6. The upland system is shown in Figure 5, with the corresponding downstream valley section in Figure 6. The water "into" the system, precipitation (Pu), falls onto the landscape and produces surface water (Qu), which is then "lost" through evaporation (Eu), infiltration (Iu or AIup), or runoff (Qr). Normal infiltration (Iu) from Qu is estimated, as well as additional infiltration (AIup) from surface spreading (SSr). This infiltration fills the soil moisture stocks (SMup), which is then available for recharge (Ru) to the shallow groundwater aquifer (GW) and evapotranspiration (ETu). Increased ET availability increases vegetation coverage (VCa), which is checked by the additional vegetation increasing ET. VCa then provides productivity benefits to upland managers, which directly increases (or decreases) the ratio of surface spreading (the SSa portion of SSr), as does the support to land managers (S). Potential surface spreading of the landscape is specified in a ratio (SSr, here 0.2). This spreading then reduces the stormwater runoff to the valley system (Qr).  Figure 6 shows the valley system. The water inputs are Qr from the uplands, the precipitation (Pv) which falls onto the valley itself, and compact allocations (Cl and Cd, which are filled in black). Cl and Cd determine the surface water into the system that varies according to precipitation availability, which is modeled here as a distribution from an upland reservoir, a common condition. Water is lost through the surface water that flows out of the valley (Qout) which is affected by the downstream compact allocations, evaporation (Ev), evapotranspiration (ETv), infiltration (Iv and AIv), and groundwater pumping out (GWout) of the shallow groundwater aquifer storage (GW). Here, surface spreading of stormwater (SWSr) is either based upon a ratio (DF) of the ditch irrigation infrastructure to recharges directly into the GW, or is spread onto the fields. This increases the surface water withdrawal quantity (W) that the farmers get based upon their compact allocations (Cl) and availability (Qa and affected by Pe). Farmers pump a minimum quantity of groundwater (GWp) based on their water rights, and also pump an amount to match whatever surface water is not delivered (see further discussion in Section 4.1).

Results
Variability in the biosphere is effected by behavior in various time scales, including the period directly preceding it (or spatial period adjacent), which is a random behavior categorized as pink noise [115]. The precipitation inputs into the system were, thus, characterized by pink noise, and results in trends above and below the norm of various durations (see Figure 7). The pink-noise random variation used a constant seed in order to produce consistent quantities with all model runs allowing for scenario comparison. The mean was 10.41 inches (264 mm)/year with a standard deviation of 1.8 inches (45.72 mm). (Note that a common standard in system dynamics models and employed in several graphs in this work is to characterize the behavior in a graph more clearly by not starting the y-axis at zero.)

Support
Support for managing flow connectivity to the landscape is provided by two functions, direct support to execute the new interventions and management, and support to managers to experiment with fitting the practices to their local conditions, thereby realizing more productivity benefits directly from the interventions. Three scenarios are shown in the results beginning in Figure 8: "No interventions" which represents the base case scenario with no additional flow connectivity management, "Support 1-execution only" which represents outside financial support to execute the new practices on a wider scale quickly, and "Support 2-experimentation and execution" which represents additional support for a phase of experimentation for a better fit to the local conditions to realize productivity benefits more quickly and to a greater extent, which leads to a perception of greater legitimacy of the practices. More upland interventions reduce the stormwater available for valley practices; thus, a tradeoff is that productivity benefits decrease in the valley with increased benefits in the uplands. . Three scenarios are shown in the results beginning in this Figure: "No interventions" which represents the base case scenario with no additional flow connectivity management, "Support 1-execution only" which represents outside financial support to execute the new practices on a wider scale quickly, and "Support 2experimentation and execution" which represents additional support for a phase of experimentation for a better fit to the local conditions to realize productivity benefits more quickly and to a greater extent, which leads to a perception of greater legitimacy of the practices.

Strategy 1: Overall Balance of System Management of Stormwater to Increase Downstream Benefit While Not Negatively Impacting Downstream Users
Upland system. Mitigating high-energy flows that transport sediment downstream is a central strategy, which then yields additional critical ecosystem service benefits to the upland ranching managers. An optimal target identified by the FlowCon framework for this system (which would be different in other systems) is to reduce the runoff quantity by 35%, as that target correlates to sufficient reductions in the peak flows to achieve the desired functions. Approximately 15% of the diverted flow on average is infiltrated, therefore, to result in this reduction of runoff the initial targeted surface spreading water quantity is approximately 120,000 AF (Acre Feet). As shown in Figure 9a,b, surface spreading in the uplands (SSa) reduces stormwater runoff (Qr) to the targeted amount by the Support 2 scenario. . Surface spreading in the uplands (SSa) begins with an optimal target identified by the FlowCon framework for this system (which would be different in other systems) to reduce the runoff quantity by 35%, as that target correlates to sufficient reductions in the peak flows. SSa (a) reduces stormwater runoff (Qr) (b). Approximately 15% of the diverted flow on average is infiltrated and results in reduction of runoff; therefore, the initial targeted surface spreading is approximately 120,000 AF (Acre Feet), as is shown achieved by the Support 2 scenario (a). The reduction of 35% of Qr is also shown achieved by the Support 2 scenario (b).
Valley farming system. Increasing water availability for the downstream farmers is the central valley goal. The valley system of farmers is downstream to the upland system and upstream to Qout. While stormwater runoff (Qr) is reduced by upland interventions, the combined effects are still enough to reverse the trend of shallow groundwater aquifer decline, as shown in Figure 10. While the differences might appear insignificant, typically only a ratio of the groundwater aquifer storage is usable and, as levels decline, the corresponding water quality generally also declines, decreasing that ratio. Note that each system would have different base conditions and different targets. No interventions Support 1 -execution only Support 2 -experimentation and execution Figure 10. This system is structured with a trend of modest groundwater declines to mimic likely site and common aquifer dynamics, as can be seen in the blue line and noting the fine scale of the y-axis starting at 19 million AF. While stormwater runoff (Qr) is substantially reduced, as shown in Figure 9, the combined intervention effects considered legitimate are still enough to quickly reverse the trend of shallow groundwater aquifer decline. While the differences might appear insignificant, typically only a ratio of the groundwater aquifer storage is usable and, as levels decline, the corresponding water quality generally also declines, decreasing that ratio. Each system would have different base conditions and different targets.

Downstream users.
At the most downstream point in the model, the quantity of SW out of the valley (Qout) delivers the compact downstream allocation (Cd) to users downstream of the system (Figure 11). In dryland areas, this allocation is typically supplied by surface water that is distributed into the system from an upstream reservoir and, thus, the majority of the water supply is not from stormwater. Figure 11 shows that the Qout quantity does not change substantially, and it is well above the Cd amount of 679,000 AFY.  . Stormwater (SW) out of the valley (Qout) is reduced in a scenarios resulting from management targets of reduction of higher-energy storm events that would carry sediment transport. The reduction does not impact the downstream compact allocation (Cd) quantity (679,000 AFY).

Strategy 2: Increase Productivity/Water Availability to Agriculture and Reverse Trends of Groundwater
Storage Declines.
Upland ranching system. The additional infiltration in the floodplains of the valley bottoms results in additional vegetation, as shown in Figure 12, which in turn fuels surface spreading (SSa) which produces additional infiltration (AIup) to the normal infiltration (Iu) in the valley (as shown in Figure 5).  Figure 5).
Valley farming. The valley farmers' production is directly related to water availability. As shown in Figure 13, withdrawals (W) increase with the stormwater spreading (SWa) interventions, reducing any surface water supply compact allocation gap, and thus reducing the need for groundwater pumping (GWout), known as "in-lieu recharge". The recovery policy ratio (Rp) for farmers able to capture and reuse groundwater is 0.75. This resulting effect of increasing groundwater levels and storage benefits the valley farmers with greater water availability in the form of W. Both support scenarios shown in red and orange produce similar results.

Scenarios for Achieving Transformation of Trends of Resource Declines
Two main critical policy/intervention dynamics in the uplands and valley management systems are addressed through highly aggregated support structures to land managers. As can be seen in the uplands structure in Figure 5, the benefits that the local managers would yield (BLu) is the more developed structure in the model. Those benefits ultimately support mainstreaming, the adoption of resilient practices into daily management [43], in these primarily range management systems. An increase in actual vegetation coverage represents forage increases, which is both the most critical productivity objective and the most challenging to achieve. Here, the innovation is to manage the range landscape through the critical driver of flood flow dynamics. Reconnecting floodplains would allow for two management improvements. Firstly, the larger floodplains can be developed as pastures are in irrigated areas providing direct productivity gains of grass responses, which can begin to be realized in the first season. Secondly, improved productivity in the floodplains would allow for more opportunity for longer recovery periods of upland rest from grazing. As raised earlier, floodplain reconnection is not a widespread resource management strategy, particularly in uplands, as the hydrologic energy behind flood flows is extensive, and ways to harness the water quantity without incurring damage from the energy represent a challenge. Support for land managers to experiment with fitting the practices to their local conditions is critical, and the details of the scenarios are discussed in further detail in section 3.1. The complexity of every ecosystem and the management approaches and opportunities provide unique conditions and needs. Here, the support to experiment is represented by the reference mode of behavior of productivity benefit uplands (PBu) and the productivity benefit valley (PBv), as shown in Figure 14. The support scenario 2 experimentation to fit the practice to the local conditions would yield productivity benefits more quickly and to a greater extent, leading to a perception of greater legitimacy of the practices. Additionally, it would provide local collaboration, which would lead to increased innovation and greater adoption. The support scenario 1 of outside support to simply shoulder the management burden would experience a longer learning curve and greater adaptive management inputs to correct unintended consequences, yielding slower and less robust benefit yields over time. An important strategy is for the land managers to experience benefits as quickly as possible.
In the valley farming systems, the benefit modeled is the result of scenarios of an increase in water availability through stormwater spreading (SWa), resulting in both the shallow groundwater aquifer recharge and additional water supply directly to fields. The spreading results in additional infiltration (AIv), which then reduces the need for surface water withdrawals, which then reduces any potential water delivery gap (Cg), resulting in "in-lieu recharge", i.e., less groundwater pumping (GWout).
(a) (b) Figure 14. The reference modes of behavior of productivity benefit uplands (PBu) (a) and productivity benefit valley (PBv) (b). For both figures, the support scenario 2 experimentation to fit the practice to the local conditions would yield productivity benefits more quickly and to a greater extent, leading to a perception of greater legitimacy of the practices. Additionally, it would provide local collaboration, which would lead to increased innovation and greater adoption. The support scenario 1 of outside support to simply shoulder the management burden would experience a longer learning curve and greater adaptive management inputs to correct unintended consequences, yielding slower and less robust benefit yields over time.
Support to land managers (S) from outside the system represents the result of collaboration across scales and can represent a range of initial injections to sustained support. Initial injections are commonly in the form of financial and technical support through vehicles such as grants. Sustained support generally requires an ecosystem service payments approach where beneficiaries support, either directly or through the state, the land managers in shouldering the management burden over the long term. The structure in this model is a constant ratio of direct support to accomplish a ratio (20%) of the potential maximum amount of floodplain area that can be reconnected and the surface spreading potential ratio (SSr), which is 20% of the landscape; thus, the floodplains reconnected in this model represent 4% of the upland landscape.
Working landscape communities of rural farmers and ranchers likely require sustained support for planning, testing, installation, monitoring, and adaptive management phases. The initial responses of the system to the interventions in this model take into account that the interventions start on smaller scales and expand. It takes time, experimentation, and expertise to build innovative approaches. The support included here assumes a strategy of ecosystem service payments from downstream beneficiaries to addresses this issue [116]. This modeling framework can be used to identify the beneficiaries of an alternate scenario and develop a collaboration or negotiations to fund or incentivize the needed management changes. The typical application is that (in this case, downstream) beneficiaries of (upstream) restoration would pay a significantly lesser amount to address the problem before it develops. A prominent example is when farmers in upstate New York felt significant economic pressure to survive and responded with intensifying their farming and livestock management. By the end of the 1980s, the increased fertilizer use and general production intensification started contaminating the New York City (NYC) water supply. NYC's Commissioner of Environmental Protection and Director of the Water and Sewer System, Albert Appleton, created an urban-rural collaboration to support farmers to voluntarily maintain sustainable practices in combination with purchasing and protecting key upstream properties [117]. As Appleton and others identified, the key for communities of any resource level, and particularly those with less than NYC, is to identify the costs of the problem and compare the costs of solutions to show evidence that they can legitimately address the problem over the long term.
In communities without significant economic resources, non-cash rewards show promise to have greater effect [118][119][120]. In New Mexico, as in many states across the West, water policy is driven by managers' need to meet compact requirements, and allowing upstream managers such as ranchers to divert and infiltrate water unilaterally could threaten the ability to meet those requirements. Evidence of the application of these approaches where no downstream impairment occurs could open up opportunities for support for water management that increases landscape-scale soil and groundwater storage over the long term. Upon such evidence, water rights or water use through a water-banking program could be used as a non-cash reward.
Land restoration on a large scale could also be potentially fueled by carbon credit or tax systems. These systems would pay land managers to increase their vegetation productivity to result in increased carbon storage. Carbon is a measure of ecological conditions increasingly used by scientists and regulators globally. In California, the carbon credit market is currently funding restoration projects and, historically, the Chicago carbon market did as well. Studies showed how carbon markets could become significant incentives to fuel restoration [121]. A 2015 study determined that, while carbon markets alone would not likely incentivize farmers to convert to their fields to riparian floodplains, it could repay 100% of the costs to do so, showing the potential for a hybrid solution of riparian buffers and a significant fund source [122].
A basin-scale stormwater harvesting and recharge strategy is likely far less expensive than the "gray" or "hard" engineering approaches typically used for flood management infrastructure [123,124]), as well as yielding far more economically advantageous ecosystem services [123,125]. However, the management burden on the agricultural community would likely be significant, particularly in early experimentation and adoption. To realize innovations in management structures to the extent that can transform the natural systems to a resilient condition, support is critical. The action research approach of building collaboration within and across scales garners outside support for both shouldering the management burden and increasing innovation capacity to yield local productivity increases considered legitimate [11,77,126].

The Hydrology-Vegetation Feedback Dynamics
A system dynamics modeling approach and its ability to capture feedback dynamics is particularly useful for characterizing the strong link between hydrology and vegetation. In this model, the interventions result in changes in the landscape vegetation conditions over time, which then changes the landscape response to flood flows. Researchers linked threshold densities of dryland vegetation occurring in bands and patches to land productivity, carrying capacity, and resource health [15,[127][128][129][130][131][132]. Reference conditions have greater vegetation density and diversity on connected floodplains than the adjacent hillsides [133,134]. The extent of floodplain vegetation density is, therefore, a key indicator of connectivity that has the potential to define thresholds for resilience in dryland systems. The hydrology-vegetation interactions are captured here with two main feedback dynamics. Firstly, additional water infiltrated resulting from the additional spreading of flood flows across floodplains (AIup) supports increased vegetation growth (VCa), which in turn supports additional infiltration (Ifp), i.e., a positive feedback loop (as can be seen in Figure 5 or A1). Secondly, the increased vegetation growth is balanced by the function that the additional soil moisture available for evapotranspiration (ETup and Eta) is inevitably consumed by additional vegetation, providing a limit to growth. The spreading of flood flows occurs in the floodplains of the upland landscape, which has higher infiltration rates (Ifp) than the average of the landscape (If), allowing for separate functions to recognize discrete spatial dynamics.

Managed Aquifer Recharge to Reverse Groundwater Storage Declines
The storage and distribution of the water through river systems around the world are managed for downstream users. Agriculture has the greatest quantity needs and the greatest controls over what gets recharged back into the groundwater aquifer systems. The approach in this work assumes that providing agriculture with the ability to conduct the innovative practices of managed aquifer recharge (MAR) is the best approach to reverse the trends of declining groundwater storage in aquifers globally over the long term.
MAR refers most generally to the concept of land management with the purpose of recharging aquifers and includes a variety of practices, with the most common application being surface spreading [135]. Aquifer storage and recovery (ASR) is the term most commonly associated with MAR, yet this refers to a specific application which involves injection wells and rights for users to then recover the water through pumping and use. The approach here primarily looks at agricultural systems as the actors of MAR, although this model could easily adapt to incorporate other MAR techniques, such as ASR.
Generally, agriculture does not take advantage of stormwater flows, often due to water rights limitations, which are typically limited to the quantities of water distributed from reservoirs through water compact agreements. This model includes a dual system of stormwater and traditional use, whereby compact allocations drive the distribution of surface water upstream from the system, as shown in Figure  15. These allocations are affected by upstream precipitation/snowpack variability, where any quantity below a mean level set by the compact agreement reduces the allocation, as shown in Figure 15b. Stormwater spread in this management and policy approach serves to balance the tendency to decline the groundwater levels. In the current system, supply to farmers are the withdrawals from surface water, which are initially the compact local allocation (CI) as reduced by the effect of precipitation (Pe) (as can be seen in Figure 6). The reductions from precipitation variability create a gap in the water right and what is allocated (Cg). This gap is generally filled by additional groundwater pumping (GWout). The stormwater actual spread (SWa) is the policy intervention, and it provides two functions for reversing shallow groundwater aquifer (GW) levels: recharge directly into the aquifer, and additional supplies of stormwater used for fields in lieu of pumping. Future work on the model intends to further develop the characterization of the contribution of recharge in the upland floodplains as subsurface flow to the valley shallow groundwater aquifer.

Conclusions
Many drylands globally are experiencing the tragedy of the commons dynamics found in the American Southwest, where management in the face of the disturbances (drought, flooding, and sediment transport) results in diminishing of natural resources and system "underachievement." While this common system behavior is indisputably a powerful phenomenon, Ostrom proved that the dynamic can be thwarted through the collaborative actions of land managers across scales. This model proposes a land management archetype that fits into the category of "seeing and acting holistically" to restore socioenvironmental capacity for resilience to the system. Flood flow is a powerful and underutilized leverage point in water-limited dryland systems. Transforming the management system to harness this opportunity on a large scale faces several challenges: what practices can successfully spread flood flow without incurring damage from the energy, who will do the work of installing stormwater harvesting practices across the scale of the landscape, and what will the incentives be? Ultimately, how could large-scale restoration realistically be achieved? This model provides a water management framework to estimate mitigating the trade-offs of upland ranching management using just enough upland flood flow to mitigate flood energy, while allowing enough to flow to the valleys to use for aquifer and in-lieu recharge. The ultimate model objective is the goal-seeking behavior of building both the natural system resource capacity and the social management system capacity to innovate, thereby achieving a resilient state. Managing flood flow connections to the landscape by slowing and reconnecting flow to floodplains mitigates hydrologic energy and results in replenishing soil storage and recharging groundwater systems, which provide critical buffers against disturbances, essential for long-term resiliency. The model is generalizable and adaptable to dryland working landscapes globally, and it serves as a collaborative starting point for working with land managers to both characterize the system under study and project scenarios of what leverage points will achieve the desired capacity-building goals. Nearly 40% of the global land surface is managed in agriculture, and, with innovative management of system functions, agriculture can revive traditional functions and once again become a system for recharging our aquifers and restoring our watersheds.