Complex Socio-Ecological Systems: Translating Narratives into Future Land Use and Land Cover Scenarios in the Kilombero Catchment, Tanzania

The Kilombero wetland in Tanzania is affected by advancing land use and land cover changes (LULCC), where we observe a conflict between development interests and the necessity of conservation measures to maintain the functionalities of the ecosystem. Thus, assessing patterns of LULCC is crucial to foresee potential future developments and to develop sustainable future management strategies. In this study, we use a multi-method scenario approach to assess the spatial implications and underlying driving forces of potential change by (1) developing a System Dynamics Causal Loop Diagram (CLD) to disentangle the underlying socio-economic and ecologic driving forces, (2) deriving a qualitative business-as-usual (BAU) and a conservation scenario from participatory narratives elaborated during a stakeholder workshop, and (3) quantifying the spatial implications of these scenarios with the Land Change Modeler (LCM). Results indicate that under the BAU assumption only 37% of the natural vegetation is expected to persist until 2030 in the wetland. In contrast, strict enforcement of protected areas (conservation scenario) halts further conversion of the wetland. However, both scenarios pinpoint considerable expansions of cropland in the western highlands with potentially serious impacts on catchment-wide hydrological processes. The produced qualitative and quantitative outputs reveal hotspots of possible future change and starting points for advisable further research and management interventions.


Introduction
Anthropogenic land use and land cover change (LULCC) is a fundamental component of global change [1]. Primarily driven by the accelerating societal demand for ecosystem goods and services, more than 50% of the terrestrial surfaces are already anthropogenically transformed [2]. This has profound implications on the climate [3], soils [4], water resources [5] as well as biodiversity and ecosystems [6]. On the one hand, the production of food and the use of natural resources are fundamental to human development and the achievement of the Sustainable Development Goals (SDGs). On the other hand, the ongoing conversion and degradation of natural vegetation are at risk to ultimately erode the ability of the biosphere to permanently provide ecosystem goods and services, thus risking the long-term sustainability of human livelihoods and well-being [7]. The assessment of potential future LULCC and its underlying dynamics is therefore crucial for ensuring sustainable land management, development strategies, and conservation measures [8]. LULCC arises from dynamic and multi-scalar interactions of biophysical, socio-economic, and political systems [9]. Scenario analysis is widely applied to investigate possible future pathways of such coupled socio-ecologic systems [9][10][11]. Rather than giving accurate forecasts or predictions, the objective of scenario approaches is to investigate a variety of alternative future trajectories based on coherent and internally consistent assumptions about the development of key driving forces [8,10]. Diverse approaches, frameworks, and methods for scenario analysis exist [10,[12][13][14]. Commonly, a distinction is made between qualitative scenarios expressed in narrative storylines and quantitative scenarios usually based on numerical computer models. Qualitative LULCC scenarios provide the benefit that they enable the investigation of complex interrelationships and interdependencies [13,15]. Moreover, qualitative scenarios easily facilitate the communication between scientists, decision-makers, and the public both during the scenario development process and in the presentation of results [15]. Participatory scenario building is increasingly applied since the inclusion of different perspectives and insider knowledge considerably enhances the relevance and legitimacy of scenarios [8,12,15]. At the same time, however, qualitative scenarios are criticized firstly for their difficult reproducibility and secondly for their lack of numerical information, which is considered to prevent further scientific and decision-making processes [12,15]. These criteria are in turn met by quantitative scenarios. Quantitative LULCC scenarios examine potential futures through numerical simulations of the rates, types, and tempo-spatial allocations of land use/land cover (LULC) related to the investigated development paths. The comprehensive modeling process produces detailed and reproducible numerical information, but risks the exclusion of stakeholder knowledge and the loss of communicability [9,13]. Therefore, a combination of both well-elaborated and described qualitative scenarios and quantitative scenarios is well recognized as it makes studies more holistic and robust: The contextualization and inclusiveness of a quantitative scenario can be enriched through qualitative information while a qualitative scenario can become more plausible and meaningful through substantiation with quantitative information [10,[13][14][15][16]. Yet, the integration of both qualitative and quantitative scenarios poses a methodological challenge which Alcamo (2008) [15] denotes the reproducibility problem and the conversion problem. The reproducibility problem relates to the fact that the assumptions and logics behind the qualitative scenario storylines are often based on specific participatory processes and therefore difficult to replicate. The conversion problem arises because scenario storylines cannot be directly transformed into numerical model inputs. In order to improve scientific credibility, the transformation process of the underlying assumptions and scenario logics therefore requires a transparent translation process [10,11,13].
This study presents an approach for transparently combining qualitative and quantitative scenario building methods using the Kilombero catchment in Tanzania as an example. By combining the contextual strengths of the qualitative System Dynamics method with participatory narratives on potential future pathways in the Kilombero catchment and a spatial LULCC modeling method, two LULCC scenarios are scrutinized. Hence, the objective of this study is to investigate how participatory narratives on potential future development trajectories in the Kilombero catchment can transparently be translated into quantitative scenarios and made spatially explicit. In doing so, we are able to assess the different impacts of the scenarios on future LULC in the Kilombero catchment, a highly contested area where conflicting land-use interests and pressures collide with the need for more conservation measures to maintain the functionalities of the ecosystem. Thus, our findings build the basis for further research on the consequences of the different LULCC scenarios for biodiversity, hydrological processes, or livelihoods, hereby supporting the identification of intervention priorities and the development of adapted and sustainable future strategies.

Study Area
The Kilombero catchment is situated in south central Tanzania in the Morogoro Region ( Figure 1). Enclosed by the Mbarika Mountains and Mahenge Highlands in the south and southeast as well as the Udzungwa Mountains in the north and northeast, the catchment covers an area of 40,240 km². The mountain ranges constitute important habitats for flora and fauna, including several endemic species [17], which is why considerable parts of the catchment area are put under protection ( Figure 1). In its wide valley floor, the catchment comprises one of East Africa's largest freshwater wetlands [18]. The wetland is regularly flooded by the perennial Kilombero River with water levels strongly influenced by the inflow from the upland areas [19]. Due to its outstanding global importance in terms of biodiversity, flow regulation, and nutrient provision, the wetland has been designated as a Ramsar Site since 2002 [20,21]. Concurrently, the valley is often also referred to as the 'Breadbasket of Tanzania', as more than 80% of its steadily growing population is engaged in farming [27]. Especially the production of rice in the floodplain is of great national importance [28]. From the 2012 census, it can be deduced that 670,000 people lived in the area in 2012 with a population growth rate exceeding 3% [21]. The high pace of population growth coupled with weak land management led to an unregulated expansion of settlement, livestock, and farming areas [20,21]. As a result, the average farm size declined to only one hectare, resulting in low food security [29,30].
To improve food security, the Kilombero floodplain was targeted as one of the central clusters of the Southern Agricultural Growth Corridor of Tanzania (SAGCOT). As part of the 'Kilimo Kwanza' (agriculture first) initiative, introduced by the former president Kikwete in 2010, the objective of the SAGCOT corridor is the reduction of poverty and the improvement of food security through public-private partnerships and the commercialization of agriculture [31]. Following the Green Growth paradigm, a triple-win situation was envisioned that benefited investors, local communities, and nature conservation [32,33]. Yet, the associated expanding large-scale investments (e.g., in rice, sugarcane, or plantations) constitute further competitors for land.
As a consequence of the high land pressure, the LULC of the Kilombero catchment has considerably changed in recent decades [5,19,27,[34][35][36]. By 2014, around 18% of the catchment's natural vegetation had been converted to cropland [36]. Wildlife populations declined considerably [20,37] and former wildlife corridors ( Figure 1) are inactive [18,19,23,34,37,38]. This applies, to varying degrees, also to the protected areas found in the Kilombero catchment [18]. In total, the protected area categories considered in this study cover an area of 1,670,296 ha (41% of the catchment) and are characterized by differing conservation goals and levels of enforcement (Table 1) [18,21,[39][40][41]. In the Ramsar wetland, for instance, rice now constitutes the dominant LULC category [36] although agriculture and grazing are legally prohibited [40] (section 32 (2) and section 34 (1b)). Facing these ongoing anthropogenic disturbances, the long-term viability of the ecosystem and the Ramsar status of the wetland are endangered [21,42]. A recent study indicates that the ecosystem service values provided by the Kilombero Valley floodplain have already decreased by more than 25% over the past decades [43].
Accordingly, the Kilombero wetland represents an illustrative example of the dichotomy between investment interests and the urgent need for food production on the one hand and the dependency of the population on an intact ecosystem on the other. Since the wetland is embedded in a larger socio-ecological system, a catchment-wide approach is recommended for investigating the degradation of the wetland [44][45][46]. Our multimethod approach for investigating future LULCC scenarios in the catchment is therefore of high relevance not only for the Kilombero Valley itself but can also be transferred to other areas globally.

Conceptual Framework
To assess the spatial implications and underlying driving forces of participatory scenario narratives, a multi-method scenario approach is applied which enables the combination of the strengths of qualitative and quantitative scenario building [13]. First, in order to construct consistent and plausible scenarios, an understanding of the complex driving forces of LULCC in the Kilombero catchment is developed by creating a System Dynamics Causal Loop Diagram (CLD). Second, qualitative scenario logics are constructed based on participatory narratives from a local stakeholder workshop and the CLD. Finally, the qualitative scenarios are translated into quantitative and spatially explicit LULC simulations using the Land Change Modeler (LCM). With this step, we are able to analyze future LULC configurations under the respective scenario conditions (Figure 2).

Data
Our research forms part of the collaborative research center 'Future Rural Africa: Future-making and social-ecological transformation' (sub-project A03: Agro Futures, www.crc228.de). Within this framework, a two-day stakeholder workshop was held at Mzumbe University, Tanzania, in February 2019 [49]. It dealt with the stakeholders' perceptions of human-water interactions as well as their expectations for future development in the Kilombero catchment, following the guiding question: How is future made in the Kilombero Valley? The workshop included three interactive sessions. First, applying a systems approach, the participants determined prevailing problems, their driving forces and consequences. Second, in a participatory mapping exercise, the stakeholders localized hotspots of current and future LULCC. Third, employing a scenario-building technique, narratives on potential future development pathways for the Kilombero Valley were identified. The 26 Tanzanian workshop participants were purposely selected from different disciplines, targeting representatives engaged in research and academic institutions, national and regional government, local authorities, the SAGCOT initiative, and non-governmental organizations (NGOs). Thereby, different perspectives from conservation, ecology, hydrology, economy, and the social sciences were included.
The central baseline data set used in this study consists of two LULC maps derived from multitemporal metrics of 2004 and 2014 Landsat imagery by Thonfeld et al. 2020 [36], who also provide a detailed description of the methodological processing. The LULC maps provide satisfying overall classification accuracies of 73% (2004) and 71% (2014) respectively [36]. Table 2 provides an overview of the data utilized in this study. Wildlife corridors Location and extent of the Nyanganje and Rupia corridors were retrieved from data provided by [18] and supplemented by literature [22,23,37]. The Mngenta Corridor was digitized from [22] (p. 49) and [23] (p. 58). Planned irrigation schemes, planned dams, planned road Results from the participatory mapping exercise conducted during the stakeholder workshop [49]

System Definition-Causal Loop Diagram for Understanding the Complex System of LULCC
For constructing internally consistent and realistic scenario logics, prior identification of driving forces and important factors shaping the investigated developments is essential [10]. Therefore, an approach is required that enables the systematic portrayal and analysis of the cause and effect relationships and feedbacks between the socio-economic, political, and ecological factors contributing to LULCC in the Kilombero catchment [8,45]. Systems thinking allows the reduction of complex phenomena to their main components in order to understand their functioning and interrelations, predict their behavior, and identify opportunities for intervention [52]. Mallampalli et al. (2016) [13] and Kelly et al. (2013) [53] give an overview on different approaches for system analysis. Among these, qualitative System Dynamics was identified as most suitable for this research as it enables an in-depth understanding and description of complex systems through the identification of reciprocal interconnections and feedback relationships between system factors. System dynamics was initially developed by Forrester (1961) [54] in the 1960s for the simulation of industrial and urban dynamics and has been further elaborated ever since [55][56][57]. Today, due to its interdisciplinary scope, it is widely applied in numerous research fields, especially in hydrology and water resource management [58,59], but also in the context of LULCC [60,61] and scenario analysis [13,53,57]. The key strength of qualitative System Dynamics is the graphic representation of the system components, their causal relationships, balancing and self-reinforcing structures as well as central system determinants in Causal Loop Diagrams (CLDs) [58]. A CLD consists of several internal and external factors (or system components). Internal system components are interactive within the system boundaries whereas external factors influence the system though they are not interactive and outside the system boundaries [57]. Within the CLD, causal links between the system components are marked by arrows. The arrows express positive or negative correlations. For instance, if an increase in factor A leads to an increase in factor B, there is a positive relationship between the two factors (change in the same direction). If an increase in factor A leads to a decrease in factor B, the relationship is negative (change in the opposite direction) [56]. In this study, a CLD was developed to analyze and more profoundly understand interconnections and feedback mechanisms leading to LULCC in the Kilombero Valley. Following Coyle (1996) [56], a three-step approach was applied. First, the problem setting, as well as the scope and boundaries, of the study were defined by a review of the literature and an evaluation of the workshop results. Second, the system components of the CLD and their connections were determined using a mixed-method approach combining participatory methods and literature research [13]. In the framework of the workshop, the stakeholders interactively developed CLD outlines to identify current problems, their causes, and consequences [49]. These insights were conceptualized and aggregated to a broader CLD on LULCC in the region, thus representing the experts' perceptions in a quasi-participatory manner [59]. On the other hand, extensive literature research of both peer-reviewed, as well as grey literature and reports, allowed for substantiation and extension of the CLD [13]. Particular emphasis was put on ensuring that the evaluated literature had an explicit spatial focus on the study area. Using the software Vensim PLE [62], the CLD was subsequently developed in an iterative process; hence several versions of the CLD were created and continuously developed. The third step regarding the qualitative System Dynamics method comprises the analysis of the generated CLD. First, system variables with the most incoming and outgoing arrows were statistically determined to identify key control factors influencing LULCC in the Kilombero catchment. By calculating 90% quantile of all factors, the upper 10% of the system variables with the most incoming and/or outgoing arrows were detected. Second, underlying causal structures and feedback loops leading to and resulting from LULCC were identified. "A loop exists when, starting from a given variable and following arrows in the direction they lead, it is possible to get back to the start, without going through any variable more than once" [56] (p. 21). There are positive and negative feedback loops. Negative feedback loops (B) contain an odd number of negative and positive links [56]. They are self-limiting and have a balancing effect on the system [57]. Positive feedback loops (R), on the other hand, contain either solely positive links or an even number of negative links [56]. They are also referred to as self-reinforcing or exponential loops [57]. If these self-reinforcing loops set desired processes in motion, targeted stimulation of such loops is possible. However, if they increasingly aggravate the situation in a system, corrective measures can be developed to interrupt or balance these feedback processes [56]. Complex CLDs may contain several hundreds of feedback loops. Therefore, the analysis must focus on the most relevant causal loops only [63]. Based on the preceding statistical analysis of central system components, the most relevant feedback loops representing key system mechanisms were identified through a combination of both the automatic loop detection provided by the software Vensim PLE together with a visual examination performed by the researcher.

Qualitative Scenarios-Deriving Scenarios from Participatory Narratives and the CLD
In this study, scenarios were derived from participatory narratives identified in the stakeholder workshop [14,49]. Stakeholder involvement is a crucial factor to make scenarios more legitimate and plausible as it enables the inclusion of different points of view as well as local knowledge [11]. Applying a participatory scenario building technique, the 26 stakeholders discussed their assumptions and expectations for the future development in the next 10-20 years. In the course of the discussion, the participants concluded on four future narratives, namely the conservation narrative, the agricultural intensification narrative, the livestock intensification narrative as well as the hydropower and dams narrative. Combining these narratives with the central findings of the CLD, two spatially representable scenarios were derived and further investigated (Sections 3.1 and 3.2).
To construct the underlying assumptions and scenario logics, the preliminarily elaborated CLD, in which driving forces, feedback structures, and key elements were determined, served as an important basis [11,13]. Information from the workshop was combined with literature analysis [13,14]. References of particular relevance for the construction of the scenario logics included studies that have likewise developed scenarios for the Kilombero catchment [5,20,64,65]. In essence, by combining researcher-driven and participatory scenario building approaches, both normative and explorative approaches could be included in the scenario logics (see also Section 3.2, especially Table 4).

Quantitative Scenarios-Spatial LULCC Scenario Modeling
For locating and quantifying the implications of alternative development trajectories for LULC, the narrative-based qualitative scenarios were translated into spatially explicit scenarios. Thereby, the previously developed CLD and the workshop results of the participatory mapping exercise supported the contextualization [13,49]. The simulation of the LULC scenarios was performed using the Land Change Modeler (LCM) [66][67][68] integrated into the IDRISI GIS-Software TerrSet (version 18.31) [67]. The LCM is widely applied in LULC scenario modeling such as deforestation [69], urban growth [70], conservation [71,72], and water resources modeling [5].
The LCM modeling process is divided into three consecutive steps. First, dominant LULC transitions and patterns of change were identified by performing a change analysis between two historic LULC maps of 2004 and 2014 [66]. Given the available data, those years represent the current changes in the Kilombero catchment most appropriately [5].
Second, two transition submodels representing the identified major LULC transformations were calibrated deploying the multi-layer perceptron (MLP) neural network [68,70,71]. The submodel 'cropland' included all LULC class transitions from natural vegetation to cropland while the submodel 'rice' included all transitions from natural vegetation to rice. A distinction was made between rice and cropland due to the special characteristics of rice and its outstanding role in the floodplain area [5]. Submodel calibration was performed through an empirical evaluation of spatial explanatory variables considered to influence LULCC [72]. The explanatory variables were determined through the application of both qualitative and quantitative methods [9]. First, the preliminarily developed CLD, which depicts factors contributing to LULCC, was examined [61]. Spatially representable system variables were identified, and available region-specific GIS data were collected. As the CLD was developed based on causal loops identified within the stakeholder workshop, the participatory insights were thereby indirectly incorporated into the model [13,49]. Second, literature on driver variables frequently used in the LCM modeling process was reviewed [5,61,[69][70][71][72]. Consequently, eight explanatory variables were included for calibrating both submodels: (i) A digital elevation model, (ii) slope gradients, (iii) soil types, and (iv) the evidence likelihood of change for each LULC category as well as distance variables such as (v) the proximity to rivers, (vi) the proximity to roads, (vii) to settlement structures and (viii) to already established rice fields (submodel 'rice') or cropland (submodel 'cropland') as of 2004. The MLP was run in the automatic training mode including dynamic learning rates and a selection of 10,000 sample pixels per LULC category. The process was repeated multiple times, continuously testing combinations of explanatory driver variables and changing learning rates. The MLP predicts the potential of a pixel to transition derived from the values of the explanatory variables for the respective pixel [71]. It models based on samples of pixels that went through each transition investigated in the submodel as well as on samples that were eligible for the transition but did not go through it. After, the sample pixels are randomly assigned to two groups: Half of the selected pixels are used to calibrate the model; the other half is used to validate how well the model predicts changes expressed through the accuracy rate [68,71]. Since some of the pixels are correctly assigned by pure coincidence, additionally to the accuracy rate also the skill statistic of the model is provided. Potential values of the skill measure can range between −1 and +1 where a skill less than 0 means that the prediction performs worse than chance, while a skill of 1 indicates an entirely correct prediction [68]. After achieving an accuracy rate of 82.62% for the submodel rice and 77.07% for the submodel cropland (skill measures: 0.8 and 0.72), the submodels were simulated to predict the most probable LULC allocation in 2030. To incorporate feedback, the dynamic explanatory variables 'proximity to settlements' and 'proximity to rice/cropland' were recalculated three times during the simulation process. The quantity of expected LULCC was determined by applying a Markov Chain analysis, extrapolating the observed trends from the 2004 and 2014 LULC maps to the year 2030 [66].
The different scenario assumptions were included in the modeling process by preparing constraints and incentives maps. These exert an influence on the prediction outputs by being multiplied with the individual transition potential maps in the course of the change prediction process. Maps with a value of 0 impose absolute restrictions, whereas a value of 1 has no effect when being multiplied. Yet, a value greater than 1 increases the transition potential thus acting as an incentive [67]. Consequently, according to the respective scenario assumptions, areas that were considered as incentives were classified with a value above 1, while areas that were considered as constraints were in turn classified with a value below 1.

Results
The results of the consecutive methods are presented following the conceptual framework ( Figure 2). First, the Causal Loop Diagram (Figure 3) with its identified key factors and feedback structures is introduced. Second, the narrative-based scenarios are outlined. Third, the quantitative LULC scenario modeling results are presented.

Socio-Ecological System of LULCC
The CLD visualizes the complex socio-ecological system of LULCC in the Kilombero catchment. To adequately account for the wetland-catchment interactions, the boundary for the diagram was drawn at the catchment level. In total, the CLD ( Figure 3) highlights four external factors (rectangular shape) and 32 internal factors. The factors are connected by 183 links (93 outgoing and 90 incoming arrows). All system variables and connections within the CLD are proven by expert knowledge from the stakeholder workshop and/or the evaluation of more than 90 region-specific publications (c). The CLD contains more than 500 feedback loops. The statistical determination of the 90% quantile revealed that there are three system main factors with more than nine incoming and outgoing arrows within the CLD, namely ecosystem health (14 arrows in total), agricultural area (12 arrows) and well-established and effective conservation areas (10 arrows). Based on these identified factors the following feedback loops were selected, hereby genuinely describing the tension of LULCC in Kilombero catchment, where there is a sensitive interplay between the integrity of the ecosystem, the growing agricultural area and the effort to conserve pristine vegetation. The CLD shows that those factors are embedded within a web of socio-ecological factors fine-tuning the main processes, for example, food security, social conflicts about land availability, or natural vegetation. By focusing on the main factors, we were able to identify balancing (B) and reinforcing loops (R) demonstrating the potential of some factors to redirect change.

Consequences of Agricultural Expansion
The first pair of positive, self-reinforcing feedback loops R1-R2 (Figure 4) concerns the negative consequences agricultural expansion may have for the environment and food security. Starting from the agricultural area, an increase in cultivated land leads to a decrease in natural vegetation. A decline in natural vegetation affects ecosystem health negatively. The worse the condition of ecosystem health, the less food security prevails. The lower food security, the higher the need to further expand agricultural areas whereby the loop closes (R1). In the long term, this leads to a self-reinforcing deterioration of the situation. The second feedback loop (R2) constitutes an extension of the first feedback loop: The more agricultural area increases, the less natural vegetation is left and the more surface run-off is generated. This in turn results in an acceleration of high and low flows. The longer and more intense high and low flows occur, the further food security is jeopardized. Hence, more agricultural areas are cultivated. The two feedback loops reinforce LULCC through an expansion of agricultural areas at the expense of natural vegetation. This leads to continuing environmental degradation and food insecurity. The references for the drawn connections within the feedback loops are provided in Table 3.

Contrasting Consequences of Conservation Measures
The second pair of feedback loops concern positive and negative consequences conservation measures may have. Feedback loop R3 ( Figure 5) demonstrates that conservation areas positively affect the environment and might be self-sustaining through the income generated from tourism. However, these effects might be constrained as conservation measures potentially stimulate LULCC outside protected areas (loop B1, Figure 5). The more well-established and effective conservation areas there are, the more natural vegetation is resulting in an increase in ecosystem health. The more ecosystem health, the more wildlife exists. The more wildlife is to watch, the more tourists are attracted. The more tourism generates income through park fees, the more money is available for perpetuating well-established and effective conservation areas. Hence, nature conservation efforts benefit the environment and can be continued through tourism (R3). On the contrary, if the situation is the reverse and protected areas are less well-established and effective, this can have negative consequences for both the environment and tourism, leading to a downward spiral.
However, isolated analysis of this feedback loop may be misleading since other interlinkages also influence the variables concerned. The negative feedback loop B1 has a balancing effect on the situation, for instance. The more well-established and effective conservation areas there are, the less land is available for agriculture. Moreover, relocating farmers from their already existing fields may be necessary in the course of enforcing protected areas. The literature review revealed that several studies report on evictions and resettlements of the inhabitants in order to create protected areas [29,73,74]. As a consequence, the population is required to cultivate new areas for their livelihoods, leading to an expansion of agricultural land outside protected areas. This further depletes natural vegetation, which has a balancing effect on feedback loop R3. Additionally, the concerned feedback loops are also interconnected with the feedback loops R1 and R2 described above via the factors agricultural area, natural vegetation, and ecosystem health. The two loops, R3 and B1, highlight the fact that feedback loops are not isolated sub-systems but rather influenced by the other factors interlinked within the CLD. Measures such as nature conservation, which have a positive impact on the environment in one place, can lead to negative repercussions elsewhere. Spatial analysis can clarify these interrelationships.

Narratives about the Future of the Kilombero Catchment-Qualitative Scenario Development
The stakeholders' expectations of potential future developments in the Kilombero catchment are clustered in four visions: the conservation narrative, the agricultural intensification narrative, the livestock intensification narrative, as well as the hydropower and dam narrative. The scope of this study required focusing on two spatially representable scenarios. A business-as-usual (BAU) scenario was created, expecting existing trends and policy decisions to continue without any interventions or alterations. It anticipates that the intensification and expansion of agricultural land and the respective degradation of the ecosystem continue unrestrictedly as envisioned in the stakeholders' intensification narrative.
The second scenario was deduced from the stakeholders' conservation narrative which was also highlighted by the CLD, where ecosystem health and well-established and effective conservation areas were identified as system main factors (Figure 3). The conservation scenario envisages that protected areas are effectively managed and well protected. It assumes that protected area regulations are enforced and accordingly, where it is legally prohibited, no further land is converted into cropland. Existing anthropogenic land uses however persist. Anthropogenic LULC conversions into wildlife corridors are expected to be reduced. In the CLD, the conservation scenario corresponds to a strong increase of the factor well-established and sustainable conservation areas. An increase in this factor directly and indirectly affects 29 of the 36 CLD system components. The concerned factors reveal the manifold and intertwined consequences a strengthening of the protection efforts in the Kilombero catchment would have. Moreover, the portrayed feedback loops ( Figures  4 and 5) outline some of the potential repercussions of the conservation scenario.
Both scenarios further indirectly comprise the hydropower and dam narrative since three planned dams (namely Mpaga Dam, Ruhudji Dam, and Mnyer Dam) indicated by the stakeholders during the participatory mapping exercise were also included in the scenario logics. Due to data unavailability, the livestock intensification narrative was excluded from the analysis.
The time horizon of the scenarios was set to the year 2030 since the narratives of the workshop participants referred to the next 10 to 20 years. Thereby, a 10-year period could be covered, and the outcomes are consistent with the Agenda 2030 Sustainable Development Goals timeline. The detailed underlying assumptions of the two scenarios are presented in Table 4. Table 4. Scenario assumptions for the BAU and conservation scenarios. Each category in the left column as well as words in italics refer to system factors of the CLD.

Conservation
No enforcement or expansion of existing protected areas. The encroachment of the Kilombero Ramsar Site and GCA as well as other protected area categories through cultivation and livestock keeping continues.
All conservation areas are managed and protected effectively according to their legal protection status. Anthropogenic uses are limited to the activities permitted in the respective protected area categories.

Population growth and settlement areas
Following the current trend, the expected high rate of population growth leads to an expansion of settlement areas and a rising demand for natural resources.
Population grows similarly to the BAU scenario. However, settlement areas are successfully prevented to expand into protected areas where the regulations prohibit such activities.

Agricultural areas
Agricultural food production continues to increase consistently with population growth. The extension of agricultural areas follows the observed trends of the last decades, thereby continuing the unrestricted expansion of cultivated land at the expense of natural vegetation.
The demand for agricultural land corresponds to that in the BAU scenario. However, the enforcement of conservation measures prevents a further expansion of cultivated land in protected areas. Existing agricultural fields remain in place and are not restored, but no new areas are converted.

Ecosystem condition and wildlife
The continuation of current activities leads to the ongoing degradation of the ecosystem and loss of ecosystem services. Consequently, food insecurity rises, habitats degrade, and wildlife declines.
Effective conservation measures mitigate ecosystem degradation and preserve important ecosystem services. Wildlife populations recover.

Tourism
The ongoing conversion of the landscape to agricultural land and the degradation of the ecosystem lead to decreasing attractiveness of the catchment and hamper its potential for tourism.
Well-managed protected areas with high wildlife populations attract tourism. Tourism in turn has a positive impact on the macroeconomy thus lowering the general poverty level. Entrance fees support the maintenance of effective protected areas.

Agricultural intensification
Agricultural intensification follows recent trends. Agribusiness investments fuel the establishment of, for example, new irrigation schemes. In the participatory mapping exercise, four planned irrigation schemes were located.
Same as BAU, but in protected areas where no agricultural activities are allowed, agricultural intensification measures such as the construction of planned irrigation schemes are prohibited to be realized.

Social conflicts
Social conflicts about land availability between largescale farmers, small-scale farmers, conservationists, and pastoralists continue.
The enforcement of protected areas intensifies land pressure and social conflicts about land availability even more than the BAU scenario.

Wildlife corridors
The expansion of anthropogenic land uses (e.g. cultivated land, grazing, settlements) continues without restrictions. Wildlife connectivity is further degraded. Protected areas remain as isolated patches in a fragmented landscape.
A further expansion of agricultural land in the important wildlife corridors Mngeta, Nyanganje, and Rupia is minimized to facilitate human-wildlife co-existence. This maintains wildlife connectivity at the current, albeit low, level.

Infrastructure development
Infrastructure and the associated accessibility of new areas act as a nucleus for anthropogenic activities such as agricultural expansion as well as tourism. In the participatory mapping exercise, a planned trunk road connecting the southern part of the catchment was expected to have a significant impact on the surrounding areas.
Same as BAU, but LULCC is restricted in protected areas.

Dams
New dams are constructed for hydropower generation and flood regulation. Dams also lead to more water resources availability which in turn leads to a further expansion of agricultural areas. In the participatory mapping exercise, the location of three planned dams (namely Mpaga Dam, Ruhudji Dam, and Mnyer Dam) was indicated.
Same as BAU.

Translated Quantitative Scenarios
Both the BAU and the conservation scenario assume that the planned road, dams, and irrigation schemes, which were identified by the stakeholders in the participatory mapping exercise, will act as a nucleus for anthropogenic LULCC in the future. The CLD ( Figure 3) visualizes these expected developments: An increase in the factor infrastructure leads to the colonization of new agricultural areas. An increase in the factor agricultural intensification which includes irrigation schemes also unlocks new agricultural areas. The construction of dams enhances water availability for cultivation. Creating an incentives map, the value 1.2 was selected to represent these expected developments. In order to account for the nucleus effect, buffers with a radius of 1 km were drawn around the respective planned measures ( Table 5).
The conservation scenario presumes that all protected areas are effectively implemented and that only anthropogenic LULCC in accordance with their legal protection status (Table 1) are allowed. Table 5 provides the constraints values assigned to the respective protected area categories. Although no GIS data exist for the Nyerere National Park, it was indirectly included in the analysis as it is located within the boundaries of the larger Selous Game Reserve (Figure 1). In addition, wildlife corridors were included with the constraints value 0.8, thus not completely preventing but reducing anthropogenic LULCC. Table 5. Results of the scenario quantification. To include the scenario assumptions in the modeling process, the constraints and incentives maps were multiplied with the individual transition potential maps in the course of the change prediction process. Consequently, a value above 1 creates an incentive while a value below 1 creates a constraint.   (Figure 6a,b).

Constraints
The scenario maps highlight two hot spots of change. First, the extent and density of cropland increases in the western, central northern, and southern parts of the catchment. This applies to both scenarios, although the expansion is more pronounced in the conservation scenario. Overall, in both scenarios, cropland expands further in almost all areas where initial conversions were already visible in 2014. Second, in the BAU scenario, rice encroaches nearly the entire floodplain area. Almost 65% of the Ramsar wetland is cultivated with rice and cropland in this scenario. Solely the areas immediately adjacent to the river banks remain covered with grassland. In the conservation scenario, in contrast, the agricultural uses in the Ramsar site and in the other strictly protected areas persist on the level of 2014. For the Ramsar wetland, this corresponds to a share of almost 45% for rice and cropland. Beyond that, areas of rice expansion are noticeable in both scenarios in the northern and western parts of the catchment, partially also in the highlands. Moreover, areas in the east, located along the Kilombero River in the Selous Game Reserve, are converted to cultivated land in the BAU scenario. In the conservation scenario, in contrast, the natural vegetation in this area is not altered.
The Markov Chain analysis revealed that in 2030, land demand for cropland is expected to amount to 8921 km². This represents an increase of 70% compared to 2014. For rice, land demand is expected to amount to 5401 km² which corresponds to an increase of 64%. Table 6 compares the calculated land demand for the two LULC classes with the areas actually distributed in the individual scenarios. The table demonstrates that for cropland, the effectively distributed area is almost congruent with the calculated land demand. Deviations account for only −0.45% in the BAU and −2.4% in the conservation scenario. For rice, however, the areas allocated in the conservation scenario differ considerably from the calculated land demand (−14%).  The impact of LULCC on the condition of the ecosystem, biodiversity, or ecosystem services is largely determined by both the intensity of change and the resulting LULCC [18]. While changes within natural LULC classes such as transformations from closed to open woodland may have inferior impacts, conversions from natural land covers to anthropogenic land uses normally have pronounced effects. Figure 7 compares the proportion of anthropogenic and natural LULC in the different scenarios within protected areas, outside protected areas, and in wildlife corridors. Anthropogenic LULC comprises the LULC classes cropland, rice, plantations, and built-up. Montane forest, swamp, water, grassland, savanna, open woodland, and closed woodland were combined and placed in the natural LULC category.
Notably, in the conservation scenario, the share of anthropogenic LULC within protected areas remains at the 2014 level of around 23%. As a counterpart, anthropogenic LULCC in unprotected areas rises by 20% and is even more pronounced than in the BAU scenario. In the BAU scenario, in turn, the share of anthropogenic LULC increases significantly within protected areas. Figure 6c indicates that these conversions particularly concentrate on the Ramsar floodplain where the natural vegetation declines to 37% ( Figure  7).
In the conservation scenario, a pixel line classified as rice and to a minor extent as cropland extends horizontally through the map. Furthermore, salt and pepper noise allocations of pixels classified as rice and cropland were observed in the background values of this scenario. An interpretation of this artifact is given in the discussion. Both, the pixel line and the scattered background pixels were excluded from further analysis of the results.

Complexity of LULCC Processes
With its 36 system components and 183 links, the developed CLD highlights the complexity of LULCC processes in the Kilombero Valley. Thereby, the catchment scale, which sets the system boundary for the diagram, enables an investigation not only of the processes within the Ramsar wetland itself but also within the surrounding and interconnected socio-ecological system. This wetland-catchment approach is requested by several studies [20,44,46,92] and makes this research a new contribution to the understanding of the dynamic and non-linear driving forces of LULCC in the overall catchment. Notwithstanding the wetland-catchment approach, it must be noted that the processes portrayed do not unfold in isolation, but are in turn influenced by other multi-scalar processes [1]. For this reason, particularly influential exogenous factors are addressed by the inclusion of four 'external factors' in the CLD. The development of these external factors lies outside the scope of the considered system boundaries and is therefore subject to great uncertainty [56,58]. Also within the CLD's system boundaries, many factors and interrelationships could not be included or could only be included in a simplified form. These factors and connections form subsystems for other research interests, such as market price logics [27], detailed influences of climate change [75], or electrification [30]. A compromise between clarity and level of detail had to be found, whereby only the most relevant and scientifically provable connections for this study were included. Therefore, CLDs are never complete representations of actual processes, but always a simplification of reality [57]. The investigated feedback loops R1-R2 (Figure 4) highlight the detrimental and selfreinforcing consequences agricultural expansion may have for the condition of the ecosystem and for food security in the long term. Consequently, these loops underline the necessity of analyzing LULCC trajectories for enabling environmental monitoring and supporting sustainable local and national decision-making [19,60]. Further, the CLD points out that factors and feedback processes are not isolated subsystems but rather connected to and influenced by the overall system [52]. When designing measures aimed at supporting or interrupting loops, the connections within the entire CLD have to be considered. In line with systems thinking for complexity, the dynamic effects of interventions in subsystems must be analyzed in the overall context.
The pair of feedback loops R3 and B1 ( Figure 5) indicates that conservation measures may have positive and growth-producing impacts in one place, but negative or balancing effects in another. Yet it is not possible to deduce where and to what extent these effects will materialize. Consequently, this pair of loops highlights the important spatial component of LULCC processes [61] which cannot be covered by the qualitative CLD [13,53,58]. Combining the contextual strengths of the qualitative System Dynamics CLD on the one hand with a spatial LULCC modeling method on the other hand proved to be a suitable approach to build and examine coherent and robust scenarios.

LULCC Scenarios
Combining qualitative and quantitative methods, two opposing scenarios envisioning future socio-economic and ecologic development trajectories in the Kilombero catchment were analyzed. The reproducibility problem and the conversion problem as described by Alcamo (2008) [15] in integrating qualitative and quantitative scenarios were successfully addressed in this study by contextualizing the underlying scenario assumptions with the CLD. Thereby, both the reproducibility of the qualitative scenarios [15] and the conversion to a quantitative model [13,53] could be enhanced.
The modeled LULC maps for 2030 provide reasonable results. Both submodels yielded good skill measures (rice: 0.8; cropland: 0.72). The western parts of the catchment and the central northern parts constitute hot spots of future cropland expansion. In a recent study, Thonfeld et al. (2020) [18] assessed the potential suitability for anthropogenic use in the catchment area based on morphometric, topographic, and bioclimatic characteristics. They also identified the western and central northern parts as especially suitable.
The scenario LULC maps show a successful differentiation between areas of future cropland and rice production. However, in both scenarios, some rice-growing areas were allocated within the western highlands. This is rather unlikely due to unsuitable climatic conditions in the highlands and since market availability and processing tools there tend to focus on more profitable products rather than land-intensive rice cultivation [5].
Another outstanding hotspot of change in the BAU scenario is the expansion of rice in the Kilombero floodplain. This is confirmed by recent studies [5,18,64] which also highlight the Kilombero floodplain as a hotspot of future anthropogenic pressure. The simulated results suggest that under BAU conditions only 37% of the natural vegetation will be remaining in the Ramsar site until 2030 (Figure 7). The future integrity of the sensitive ecosystem and the Ramsar status are highly questionable given this setup [21]. Yet the modeled extent of the expansion of the rice cultivation area almost up to the riverbed raises questions as crops too close to the riverbed are at risk of being flooded. Cultivation in the central floodplain area would require the floods of the Kilombero River to be regulated [21,36]. This would in turn have serious negative effects on the ecosystem [93].
In the conservation scenario, it is assumed that no further anthropogenic LULCC occurs in protected areas in which it is legally prohibited. The results reveal that this assumption has been successfully implemented in the model as the shares of anthropogenic land uses within protected areas do not increase in the conservation scenario (Figure 7). Outside protected areas, however, the conversion of LULC is even more pronounced as compared to the BAU scenario. Referring to feedback loop B1 (Figure 5), here the dislocation of agricultural conversion is demonstrated. Especially in the western part of the catchment, the density and extent of cropland increase markedly more than in the BAU scenario (Figure 6c,d). This is to be considered critical since the upland areas are of great importance for the water supply of the wetland. Näschen et al. (2018) [44] (p. 20) warn that "(t)he increased share of cropland, which results in a reduced retention capacity, will influence the flow regime, with declining low flows and aggravated flooding". Conservation areas such as the Ramsar site therefore cannot be considered as isolated systems, but are connected within the overall socio-ecologic system [52]. Hence, these observations justify the catchment-wide approach of this analysis.
With respect to the catchment-wide approach of this study, it has to be acknowledged that the observed changes between 2004 and 2014 were linearly extrapolated to 2030. Consequently, neither external forces such as changes in policies or (inter-)national market demand ( Figure 3) nor tipping points or irreversible changes were included. Additionally, not all influencing driver variables of LULCC as identified within the CLD could be included in the LCM submodel calibration due to either data unavailability or because the variables could not be represented spatially. Agribusiness investments in large-scale ventures or the interactions between agricultural expansion and ecosystem health, for instance, therefore, constitute sources of uncertainty in the model. The difficulty to include socio-economic data was also highlighted by other studies [5,16,72]. However, the objective of the quantitative LULCC scenario simulation was not to make pixel-precise predictions but rather to reveal trends and future hotspots of LULCC. Additionally, even though some qualitative data cannot be transferred into the model, these qualitative data are not neglected but build the context and the quantitative results are reflected within this broader qualitative setting. Accordingly, we argue that here the relevance of combining qualitative and quantitative methods for a holistic investigation of scenarios is evident.
Although the two future scenarios differ in terms of their respective spatial distribution of LULC, their absolute amount of progressing LULCC is roughly similar ( Table 6). The reason for this is that the calculated land demand (Table 6) in both scenarios is based on the extrapolated trends from 2004 to 2014. Nonetheless, in the conservation scenario, less rice (−14%) and slightly less cropland (−2.4%) was distributed than the actual land demand was calculated. This can be explained by the fact that in view of the constraints in the conservation scenario, a further expansion of agricultural areas in protected areas is not possible. Outside protected areas, however, there appears not to be enough land suitable for rice and crop cultivation. Insufficient pixels with suitable characteristics for rice or cropland are available (Table 6). This also explains the horizontal pixel line classified as rice and subordinately as cropland as well as the salt and pepper noise allocations that are spread over the LULC map of the conservation scenario. Through this over-allocation, the LCM indicates that no more applicable pixels are available [68]. Accordingly, the conservation scenario implies insufficient usable land for agriculture, even though existing farmland is not even reduced. In view of the increasing demand for agricultural products both for the growing population and for the national and international market, this raises questions regarding the future supply of the population with both food and income. The enforcement of strictly protected areas and the halting of LULCC can only be feasible if food production becomes more efficient and/or new sources of income are created in parallel. Agricultural intensification such as the application of fertilizers, pesti-cides, irrigation, recession cropping, or agricultural machinery has the potential to increase yields and thus enhance agricultural productivity [28,65]. Yet, these measures require careful consideration. The CLD (Figure 3) indicates that even without areal expansion, agricultural intensification may profoundly degrade the ecosystem through surface and groundwater contamination or soil degradation, for instance. Therefore intensification measures must be balanced and should follow eco-efficient and site-specific strategies [28]. For this purpose, the CLD is an excellent tool to evaluate the interdependencies and potential consequences of such measures. Targeted activation of desired self-reinforcing feedback loops on the one hand and the interruption of undesired feedback effects on the other can support the identification of appropriate strategies.

Conclusions
The study shows the relevance of understanding the interconnectedness of the coupled socio-economic and ecological system of LULCC and its dynamic and diverging future consequences for designing and implementing sustainable development measures and policies. In line with systems thinking for complexity, the developed Causal Loop Diagram provides a unique level of insight into the non-linear interplay of the manifold interdependent push and pull factors and feedback structures leading to and affected by LULC transformations in the Kilombero catchment. Combining the contextual strengths of the qualitative System Dynamics CLD with participatory narratives on future trajectories in the Kilombero catchment and a spatial LULCC modeling method proved to be a suitable approach to transparently build and examine coherent and robust scenarios. Thereby, the integration of participatory inputs and normative visions in the CLD and scenario building process enabled the inclusion of both bottom-up and top-down perspectives. The two scenarios examined represent diverging paths into the future. Given the BAU assumption of a continuation of all current trends in the future, a collapse of the Kilombero wetland is likely. In this setup, rice encroaches almost the entire Ramsar site, leaving only 37% of natural vegetation. Strict enforcement of protected areas as envisioned in the conservation scenario could halt the degradation of the wetland but exacerbate land scarcity and leave insufficient suitable land for agricultural production. Apart from the expansion of rice, the western part of the catchment is confronted with an immense expansion of cropland areas in both scenarios which is expected to substantially alter hydrological processes in the overall catchment. Both scenarios represent two extremes among a multitude of possible development paths. Which hypotheses will translate into reality depends on the complex interplay of the numerous interlinked factors at play in the Kilombero catchment. However, targeted policy interventions and initiatives can strategically influence future development. The study produced qualitative, quantitative, and spatial outputs that can be used as a basis for such decision-making. Our results build a starting point for further research on the impacts of the different LULCC scenarios on biodiversity, hydrological processes, or livelihoods, thereby supporting the identification of intervention priorities and the development of adapted and sustainable future strategies.

Data Availability Statement:
The data presented in this study are available in the article and Appendix A Table A1.

Acknowledgments:
We gratefully acknowledge the participants of our workshop for their contributions. We appreciate that they shared their time, expertise, and experience during the collaborative workshop exercises.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.
Appendix A Table A1. The table explains each system factor of the CLD (Figure 3) and provides the references for the outgoing arrows of each factor. Thus, all arrows are substantiated with literature.  [21,30,80,82]