A Qualitative Hydro-Geomorphic Prediction of the Destiny of the Mojana Region (Magdalena-Cauca Basin, Colombia), to Inform Large Scale Decision Making

Colombia is undergoing a period of rapid development. In particular, the Magdalena-Cauca Rivers basin, and the Mojana region within it, is going to experience rapid expansion in infrastructure growth, entailing hydropower development, road and navigability works along hundreds of kilometers of channels, as well as standard flood control measures. This paper argues that unexpected and undesired outcomes are very likely to occur as a consequence of the hydraulic and geomorphological reaction of river systems to such development schemes; namely, we foresee heightened hydro-morphological risks, along with the loss of environmental services and strong increases in maintenance costs. River behavior has been the subject of extensive study by diverse disciplines. As a result, key principles of fluvial dynamics have been elucidated and specific quantitative prediction tools developed. In this paper we do rely on this wealth of knowledge. However, since specific local information and interpretative tools in Colombia are either lacking or unreliable, it is inevitable that, at the moment, any basin scale analysis has to remain qualitative and must incorporate several assumptions, leaving it open to questioning and further refinement. Nonetheless, we argue that advancing such type of speculative conjectures is the “right thing to do”. The undeniably desirable but hard to achieve alternative of waiting for sufficient datasets and tools would entail excessive delay in obtaining relevant answers while large-scale development would continue to occur with potentially damaging results. Therefore, our analysis is conceived along the precautionary principle. This paper is primarily aimed at technical advisors of policy makers as it offers scientifically-based arguments for informing the political debate, hopefully guiding decision makers towards better choices. Rather than advocating specific solutions, the focus is on pointing out the likely adverse consequences of the currently planned course of action.


Introduction
Following a World Bank expert mission in the 1950s, there has been widespread political debate regarding several initiatives for the development of agricultural districts and modern transport infrastructure in the Mojana region, along with proposed policies for a complete control of its natural dynamics characterized by periodic harsh flooding and subsequent drought periods [1]; ecosystem conservation was also considered [2]. This process took on a renewed impulse after the "Ola Invernal" (winter wave) of 2008-2009 and 2010-2011, as well as two "La Niña" (i.e., cold ENSO -El Niño Southern Oscillation) periods with extreme rainfall, which caused flooding and severe damages in all of Colombia. The 2010-2011 event had a particularly notable effect in the Mojana region, where the Cauca River breached a large artificial levee, causing widespread damage within one of the most developed areas within that basically poor region.
As a result, the region received significant national and international emergency funding through the Colombia Humanitaria organization, while the government created the Fondo de Adaptación, a country-wide agency designed to promote economic recovery for the population affected by "La Niña" through sustainable and secure infrastructure work. Since then, there has been a vigorous societal and political debate around flood control concepts; although many proponents invoke the need for large defense infrastructures, this notion has been challenged, particularly within academia and the environmental movement.
We focus on the Mojana region (Figure 1), whose context within Colombia and the Magdalena-Cauca Rivers basin is briefly summarized in Table 1.   Considering the sheer size of the Mojana, it is not surprising that the decision-making process has developed in a multi-faceted, multi-channel fashion, involving several social arenas and a number of parallel leading actors with little or no coordination. The breadth of agencies involved in the process include: the National Planning Department (DNP); the Ministry of Environment and Sustainable Development and its Institute of Hydrology, Meteorology and Environmental Studies (IDEAM); the Ministry of Interior and Justice; the Ministry of Economy; the Ministry of Agriculture and Rural Development; the Ministry of Transport; the Fondo de Adaptación; the Unidad Nacional de Gestión de Riesgos y Desastres (National Unit for Risk Management and Disaster Rescue); several universities; international agencies; Non Governmental Organizations (NGOs), etc.
This has resulted in numerous, sometimes conflicting, lines of thought, with their corresponding political actions: Many participants push for the construction of new transportation infrastructure as a key element in the development of the region (part of [3,4] among others); some focus on flood control as a core component, proposing specific projects ( [5,6] among others); others emphasize the socio-economic and socio-cultural dimensions of the situation, stressing the need for empowerment of local communities, provision of basic services, and management of conflicts and illegal armed groups [7][8][9]. Finally, there are those who propose that any solution should focus on ecological capital and restoration (part of [2,3,10], some chapters of [11]; some chapters of [12,13], among others).
All of these approaches are still far from addressing the set of twelve (12) core procedural and substantive sustainability criteria proposed in [14] to be used as a guide for clarifying development purposes, identifying potentially desirable options, comparing alternatives and monitoring implementation for infrastructure.
Notwithstanding these issues, a serious coordination effort is currently underway to come up with a final Action Plan, based on the most scientifically-sound knowledge about the flood dynamics and associated risks in the Mojana region. This process is the result of a large research project conducted by a national modeling team at the IDEAM, with a select group of risk experts under the lead of the Fondo de Adaptación. Despite the very high quality of this study, two key weaknesses can be identified: The approach lacks a holistic view of the whole set of development interventions envisaged within the different planning instruments (i.e., looking at them all together).

2.
The likely effects of the different development projects on the future geomorphic evolution of the river system have not been considered nor assessed.
In this paper, we argue that seriously considering these two issues could profoundly change the development vision-and course of actions-for the whole Mojana region (and the Magdalena-Cauca Rivers basin as well). This is a different facet of the discourse about the importance of considering fluvial geomorphology to support land-use decision-making developed in [15], along with similar challenges as in several other large rivers worldwide [16]. An exercise with comparable ambitions, but with a quantitative approach, has been conducted for the Middle-Lower Parana River [17]; it is interesting to note that, in spite of the greater rigor offered by that modeling framework, the analysis still contains-and somehow hides-several subjective hypotheses and assumptions.
We want to make clear from the onset that we are not disregarding specific models; rather, we claim that the descriptive, qualitative analysis of the type we present here is a fundamental, although preliminary, step that can then be supplemented by additional modeling exercises, provided that the underlying hypotheses and limitations are made explicit and critiqued carefully. Current geomorphic knowledge and interpretative tools [1,18] are still too weak to support a quantitative assessment for the Mojana case; nevertheless the arguments presented here is hardly refutable, being based on geomorphic knowledge and world-wide experiences.
In a more rigorous scientific approach, one would prefer to obtain all the data and tools needed before announcing the type of strong conclusion we present in this paper; but waiting until those elements are available and obtained would not provide support to the decisions that are being taken now. Hence we preferred to sacrifice some level of rigor, and strive to provide qualitative, but relevant information at this time, while the decision-making processes are still under way. This is an aspect that is certainly of interest to many other cases throughout the world: the attempt to significantly contribute to the decision process by advocating scientific arguments while relying on scarce and mainly qualitative information, particularly concerning the long-term geomorphic behavior of the system considered.
Our case can be framed within the precautionary principle [19] as it deals with a decision of national scope which may induce adverse and irreversible impacts and where consequences are highly uncertain. According to the analysis framework set up in [20], we should try to separate the uncertainty related to the possible effects (scientific side) from the degree of risk acceptance (socio-political side). Coherently, we could assess the level of scientific uncertainty we manage as "reasonable belief" or "clear indication" and label the aptitude we recommend as that of the "Cautious environmentalist" (in Weiss' scale [20] (p. 144), the lowest category of 5 is the "Environmental absolutist" who says take no action until it is proven that that it will cause no harm; the highest (5th) is the "Scientific absolutist" who says make no precautionary intervention until the danger is scientifically proven; the "Cautious environmentalist" is the 2nd).
While this knowledge level would be insufficient to permit an exploitation or development project to be undertaken when proposed by a private party, it could be enough, according to [20], in case whereby society (through its environmental scientists) were to demonstrate that counteraction needs to be undertaken due to the high likelihood that an adverse effect is certain to occur. Or, said in the words of ( [21], p. 318), " . . . uncertainty erodes the traditional positivistic model of knowledge, in which science speaks truth to power"; these authors also discuss the uncertainty paradox which is based on the fact that " . . . the role of experts was framed in terms of providing certainty about uncertain risks" and point out, citing [22] in ( [21], p. 331), that " . . . uncertain risks 'defy "normal" scientific analysis, for they are indeterminate' and " . . . inseparable from the political and social context of which they are an outcome and a contributor . . . ". Hence, a more figurative approach would be the classification proposed by [23] where the classic impulse to defeat uncertainty in scientific knowledge would be labeled as "monster exorcism" and would typically lead just to nothing, or at least to time incompatible with the need to make a decision.
Coherent with this position, this paper does not claim to be a rigorously scientific contribution; it aims to offer a scientifically consistent understanding and arguments to the community and, particularly, to the technical (not scientific) advisors of policy makers. Rather than advocating particular solutions, it addresses the likely adverse consequences of the course of actions currently planned. The methodology we adopted is reflected by the structure of this paper: first, a description of the Mojana region and a rapid River Style assessment of its main hydro-morphological actor, the Cauca River, are presented in order to characterize the system. A holistic picture of the main structural interventions planned at both scales (i.e., that of the Magdalena-Cauca Rivers basin and that of the Mojana region) is then provided: the former is needed to specify the boundary conditions of the Mojana (which may affect the processes inside it); while the latter refers specifically to that territory.
A discussion of the likely future evolution of the river basin system is then developed, based on the understanding provided by the rapid River Style assessment and a qualitative cause-effect conceptual model developed for the study. Finally, management conclusions and recommendations are provided. Note that we are often obliged to jump from one scale (Mojana region) to another (whole Magdalena-Cauca Rivers basin) depending on the step of this process: without maintaining a look at the wider scale, we would lose very important relationships that can be determinant for the long term destiny of the Mojana region.

Physiography and Peculiarities
The Mojana region is not a basin, rather it is a part of a larger region called Depresión Momposina (i.e., depression of Mompox), itself part of the Magdalena-Cauca Rivers basin (see Figure 1). It is a fluvio-deltaic plain formed by the rivers San Jorge, Cauca and the Magdalena branch called Brazo de Loba. Over the course of time, the notably flat topography and the presence of large rivers generated a complex ecosystem with ever-changing wetlands, inter-connecting channels (caños) and other components (see Table 2) according to the flood pulses and complex bidirectional temporary flows. These were well exploited by ancient Zenúes civilizations through extensive land preparation, including a number of man-made small canals (see Figure 2a).
Resources 2016, 5, 22 5 of 25 that of the Mojana region) is then provided: the former is needed to specify the boundary conditions of the Mojana (which may affect the processes inside it); while the latter refers specifically to that territory. A discussion of the likely future evolution of the river basin system is then developed, based on the understanding provided by the rapid River Style assessment and a qualitative cause-effect conceptual model developed for the study. Finally, management conclusions and recommendations are provided. Note that we are often obliged to jump from one scale (Mojana region) to another (whole Magdalena-Cauca Rivers basin) depending on the step of this process: without maintaining a look at the wider scale, we would lose very important relationships that can be determinant for the long term destiny of the Mojana region.

Physiography and Peculiarities
The Mojana region is not a basin, rather it is a part of a larger region called Depresión Momposina (i.e., depression of Mompox), itself part of the Magdalena-Cauca Rivers basin (see Figure 1). It is a fluvio-deltaic plain formed by the rivers San Jorge, Cauca and the Magdalena branch called Brazo de Loba. Over the course of time, the notably flat topography and the presence of large rivers generated a complex ecosystem with ever-changing wetlands, inter-connecting channels (caños) and other components (see Table 2) according to the flood pulses and complex bidirectional temporary flows. These were well exploited by ancient Zenúes civilizations through extensive land preparation, including a number of man-made small canals (see Figure 2a). National and regional roads do exist around the region, but they are scarce and in bad condition, therefore access to the internal areas is difficult; communication is mainly by boat in the Cauca, San Jorge and Magdalena Rivers and the swamps. Some paths (e.g., San Marcos-Majagual-Achí) were built that had profound disturbance of the ecological equilibrium of the area, cutting through swamps and in some cases restricting hydraulic interconnection-resulting in more frequent floods in some areas  National and regional roads do exist around the region, but they are scarce and in bad condition, therefore access to the internal areas is difficult; communication is mainly by boat in the Cauca, San Jorge and Magdalena Rivers and the swamps. Some paths (e.g., San Marcos-Majagual-Achí) were built that had profound disturbance of the ecological equilibrium of the area, cutting through swamps and in some cases restricting hydraulic interconnection-resulting in more frequent floods in some areas [1].

Rapid River Style Assesment
Three main rivers act within the Mojana region: the great Magdalena River, which borders its northeastern side; the San Jorge River, marking the main drainage axis of the large system of pulsing wetlands (humedales or ciénagas); and the Cauca River in the south-eastern part (see Figure 3). In this paragraph we perform a rapid assessment of the Cauca River only. This analysis is followed later on by a more comprehensive examination of the whole Magdalena-Cauca basin and its other rivers, with the purpose of seeing how "the rest of the world" surrounding the Mojana region may affect it in terms of boundary conditions.
Starting with the Cauca River is justified in our view, because it is in some sense the key actor for a number of reasons:

1.
It is this river that historically produced the harshest damage due to its geomorphic behavior.

2.
It greatly influences the hydrological dynamics of the humedales (wetlands) as a result of its periodic overflowing events.

3.
It is in its basin that most of the future dams are envisaged.
Moreover, looking at the Cauca in some detail as a first step will help us to develop an understanding and basic (qualitative) predictive capability.
The key idea of the River Style framework [26] is to look at the river on different scales at the same time: from the landscape and river basin scales, down to corridor and reach scales, until that of single geomorphological units within either the fluvial corridor or the river channel itself. From this analysis, an effort is undertaken to understand the functioning of the river in different hydrological conditions, particularly during bankfull stage and large floods events. It is hence not just a classification exercise (like the approach in [27]), but a real framework for a thorough understanding that opens the door to the analyst to predict future behavior in a rather qualitative, but very illuminating fashion.  Figure 1). In the floodplain between the rivers, light blue marks the wetlands. The grey base is a 10 m pixel Digital Terrain Model (DTM 10 m). On the left, is the highly sinuous San Jorge River, partly confusing with the wetlands it feeds. The major splays (track visible in black arrows in the inset) from left bank natural levees crevasses of Cauca River-named Nuevo Mundo and Santa Ana-are located within reach 1 and 2, pouring their water and sediments northwards. Table 3 presents a set of attributes adopted to characterize the reaches under consideration according to the River Styles approach of [26]. Their assessment is still mainly qualitative and preliminary, but sufficient to capture some of the key features of the Cauca River and its territory. The following Figure 4 presents in schematic form some planform features of the river to underline its setting in terms of "freedom of movement", i.e., confinement and alignment. Figure 5 shows a representative detail of river morphology to offer a multi-scale view. Figure 6 then shows some long topographic cross sections to point out the relationship between the river and its floodplain.  Figure 1). In the floodplain between the rivers, light blue marks the wetlands. The grey base is a 10 m pixel Digital Terrain Model (DTM 10 m). On the left, is the highly sinuous San Jorge River, partly confusing with the wetlands it feeds. The major splays (track visible in black arrows in the inset) from left bank natural levees crevasses of Cauca River-named Nuevo Mundo and Santa Ana-are located within reach 1 and 2, pouring their water and sediments northwards. Table 3 presents a set of attributes adopted to characterize the reaches under consideration according to the River Styles approach of [26]. Their assessment is still mainly qualitative and preliminary, but sufficient to capture some of the key features of the Cauca River and its territory. The following Figure 4 presents in schematic form some planform features of the river to underline its setting in terms of "freedom of movement", i.e., confinement and alignment. Figure 5 shows a representative detail of river morphology to offer a multi-scale view. Figure 6 then shows some long topographic cross sections to point out the relationship between the river and its floodplain.      Table 3 (the reach identifier is reported in the horizontal colored bar at the top), which points out the degree of confinement, the alignment of the channel within its valley (dashed zone) and secondarily, its sinuosity. Dashed vertical reddish segments indicate the location of the long topographic cross sections shown in the following Figure 6. Notice the presence of a clear confined reach (4) where Caribona River joins the Cauca. Along reaches 3, 5 and 6 (where sections C-C 1 , E-E 1 and F-F 1 are located), Cauca River can move in a very wide floodplain; although the major crevasses are located along reach 1 and 2 (see Figure 2), their splays pour into reach 3 valley domain.    On the basis of this information, the quite peculiar character of Cauca River within the Mojana region can be inferred: it is a river with a very low slope (0.15 per thousand or less) and low or very low specific stream power. Its sinuosity is low to moderate (the exception being the last reach where it is higher); mainly with a single main channel, but sometimes with up to three (rarely four) channels, seemingly denoting a wandering character, compatible with its high w/h ration (even if the bankfull depth "h" data are very imprecise); accordingly, higher flow rates increase velocity and shear stress more than water depth.
The Cauca River displays very significant (several km) lateral displacements in several reaches and even avulsions (witnessed by the presence of marked longitudinal paleo-channels), facilitated by the generally very low confinement (see Table 3). Its low energy is compatible with the presence of anastomosed channels, while the seemingly wandering behavior of some reaches is indeed witnessing an avulsive character associated with the dynamics of the natural levees that accompany most of its length. Natural levees are created by deposition of coarser sediments (usually transported as bedload) during frequent, lateral overflows, while suspended load is transported further away creating backswamp deposits. Indeed, during significantly high flow events (with flow rates up to even three times as large as the bankfull flow), the occurrence of crevasses in the natural levees of Cauca River is not unusual.
This generally occurs in known preferential sites (called rompederos), but sometimes opening new ways. The crevasse flow then feeds an irregular and ever changing, very wide, splay zone, with significant water and sediment volumes, part of which is conveyed by a complex system of more stable natural channels (called caños); this whole draining system carries water and sediments to the main drainage axis where the great swamps and wetlands (humedales and ciénagas) lie. The associated deposition of sediments contributes to contrast the natural subsidence process [28,29].
The long topographic cross sections (see Figure 6) show the peculiar setting of the river on a kind of crest, with a clear West-Northwards slope of the same magnitude as its longitudinal slope (exception made for the confined reach 4). This setting indicates a river in the process of building-up its alluvial valley mainly by vertical deposition, but also by a re-working due to shifting meanders (ridge and swale topography) and partly because of splays and deposition associated with crevasses. It is exactly on the slightly more elevated spots of land left by natural levees or meanders' convex bank ridges that people settled and created the towns.
Perhaps other factors, including tectonic or differential subsidence, contribute to the creation of this topographic setting. In practice, the river channel keeps its current alignment mainly because of the presence of its own natural levees which perform a kind of guiding function; it cannot be excluded, however, that in a particular flood event, the river might suddenly change its trajectory, "plunging" into its wide floodplain towards the draining axis of the San Jorge River and its wetlands (humedales and ciénagas) and giving life to a completely new main river bed (coherently with the generally very low value of the entrenchment and pseudo-entrenchment indices of Table 3).

The Planned Future
In Section 2, we provided some insight into how the Cauca River behaves as the main actor within the Mojana region; however, it is impossible to predict what will happen there, without specifying what will happen at its border, owing to the evolution of the rest of the whole Magdalena-Cauca Rivers basin. The Magdalena River actually marks the northeast border of the Mojana; but, more importantly, even if that were a short reach, some changes in the Magdalena River may deeply transform the whole territory. For instance, incision: By regressive erosion it can lead to incision of the whole hydrological network and that can have a very significant effect.
This vast basin is going to be affected by a number of significant processes: Land use change in the river basin: this is a process which possibly already touched its apex, but it is still going on (particularly deforestation and urbanization more or less un-controlled) and whose consequences-particularly the increment of sediment load to rivers-are on-going [31][32][33]; this is accompanied by industrial and illegal mining development and its proved toxic consequences [34][35][36][37][38][39][40][41]; 3.
Population growth and its associated demand for goods and production of wastes.
In addition, a set of very heavy interventions is on the political agenda (see 2. Land use change in the river basin: this is a process which possibly already touched its apex, but it is still going on (particularly deforestation and urbanization more or less un-controlled) and whose consequences-particularly the increment of sediment load to rivers-are on-going [31][32][33]; this is accompanied by industrial and illegal mining development and its proved toxic consequences [34][35][36][37][38][39][40][41]; 3. Population growth and its associated demand for goods and production of wastes.
In addition, a set of very heavy interventions is on the political agenda (see Figures 7 and 8), including:  Hydraulic works to control floods (levees) and fluvial dynamics (groynes, rip-rap protections, gabions, etc.) in several river sectors (which, as experience demonstrates-as partly shown in the Supplementary Material-very often activate a never-ending process that affects the whole river network, while implying much higher expenditures than initially planned);

3.
Morphological re-sectioning of several caños in the region designed to ease the evacuation of flood waters and the reduction of flooding in the exploited region; 4.
Mega interventions to ensure the commercial navigability along 900 km of the Magdalena River and in some caños of the Mojana; this includes very extensive and continuous dredging as well as bank defenses (rip-rap, gabions, walls, etc.).    Figures 1 and 7).

The Foreseen Future: A Different Picture
In this paragraph, we put forward a large scale, long term prediction of the likely geomorphic effects associated with the whole array of interventions and boundary conditions portrayed above. As the current political debate is concentrated on the Mojana region, that is our focus of discussion; however, we have to consider the whole basin because interventions there may affect the boundary conditions of the Mojana, and because the basin itself (particularly the lower basin) is likely to suffer consequences, including those due to the Mojana region dynamics.
The "model" adopted to carry out such a prediction is merely conceptual and is based on known physics of fluvial dynamics (see for instance [42]) as well as a number of world-wide experiences of similar cases (the "conceptual cause-effect model", see supplementary material to this paper; see Figure 9). The model addresses three main components: hydro-morphological risk, environmental consequences and economics; we describe synthetically here just the most relevant aspects.

The Foreseen Future: A Different Picture
In this paragraph, we put forward a large scale, long term prediction of the likely geomorphic effects associated with the whole array of interventions and boundary conditions portrayed above. As the current political debate is concentrated on the Mojana region, that is our focus of discussion; however, we have to consider the whole basin because interventions there may affect the boundary conditions of the Mojana, and because the basin itself (particularly the lower basin) is likely to suffer consequences, including those due to the Mojana region dynamics.
The "model" adopted to carry out such a prediction is merely conceptual and is based on known physics of fluvial dynamics (see for instance [42]) as well as a number of world-wide experiences of similar cases (the "conceptual cause-effect model", see supplementary material to this paper; see Figure 9). The model addresses three main components: hydro-morphological risk, environmental consequences and economics; we describe synthetically here just the most relevant aspects.
As stated at the beginning of the paper, a scientifically prudent approach would prefer not to put black on white the following statements as they involve several assumptions; but the desire and need to support as far as possible the decision making process before it is too late motivated us to undertake this daring and uncertain exercise (along the line of what has been carried out by [43]). We are aware that some phenomena counteract others; at the moment it is virtually impossible to state which one would prevail. The most probable future is that they will all take place, but in different reaches and at different times. Notice that the figure is somehow aggregated: although it presents the set of dikes, weirs, and dams as a unified element, their effects are of course not the same and sometimes the arrow only carries a partial relationship. For instance, only dams and partly weirs have a sediment retention effect which hence reduces solid load; on the contrary, it is mainly levees which reduce overbank flows; and so on. Red boxes denote "pivot effects" that are key to other components of the cause-effect network, not represented in this scheme, but in others. Notice that the figure is somehow aggregated: although it presents the set of dikes, weirs, and dams as a unified element, their effects are of course not the same and sometimes the arrow only carries a partial relationship. For instance, only dams and partly weirs have a sediment retention effect which hence reduces solid load; on the contrary, it is mainly levees which reduce overbank flows; and so on. Red boxes denote "pivot effects" that are key to other components of the cause-effect network, not represented in this scheme, but in others.
As stated at the beginning of the paper, a scientifically prudent approach would prefer not to put black on white the following statements as they involve several assumptions; but the desire and need to support as far as possible the decision making process before it is too late motivated us to undertake this daring and uncertain exercise (along the line of what has been carried out by [43]). We are aware that some phenomena counteract others; at the moment it is virtually impossible to state which one would prevail. The most probable future is that they will all take place, but in different reaches and at different times.
To give an idea of why we are obliged to accept a qualitative, semi-subjective approach such as the one here adopted, consider that overcoming these kind of limitations would require a very detailed, broad scale, time consuming and costly exercise to fill several information gaps at the scale of the Magdalena-Cauca Rivers basin (271,249 km 2 ). Namely, a thorough process would probably include the following: (i) The available official time-series of flow rates and suspended solid concentrations at different gauging stations should be verified to eliminate the likely effects of subsidence and morphological changes (particularly aggradation/incision) occurred during time so far; they should be depurated by eliminating wrong or doubtful data, and-where needed-completed by reconstructing missing data; finally-and so adjusted-they should be shared and recognized as an accepted database within the Colombian scientific community. (ii) Data on granulometry should be taken regularly at several monitoring points (un-existent today). (iii) Also topographic sections in the rivers and topography of the floodplain should be surveyed systematically to create a complete database, detect changes and verify how subsidence (and tectonics or orogeny) works. (iv) A systematic, complete and updated mapping of relevant infrastructure and its status (bridges, dams, intakes, bank protection works, levees, etc.) should be made available. (v) Identification and measurement of morphological characteristics at different scales (e.g., the river valley, the floodplain, the bankfull, the islands, the bars, etc.) should be made at different time instants with an homogeneous approach-a process that requires a wide set of good quality, whole coverage, aerial or satellite photographs (a very interesting example is given in [44]). (vi) A geomorphic classification possibly based on the River Style approach (similar, but more refined than the Rapid River Assessment of Section 2.2, Table 3) would usefully identify homogeneous geomorphic reaches. (vii) Sediment transport capacity should be assessed (and suitable formulas calibrated) and bankfull and effective flow rates should be assessed and compared. (viii) Mapping of sediment sources should cover the whole region and models estimating supply-based on empirical formulas like for instance the RUSLE [45]-should be calibrated by performing historical sediment budgets of relevant reaches and reservoirs, which would take into account morphological changes and all inputs/losses. (ix) A "story of the river" for relevant, representative river reaches would be elaborated which would describe synthetically and systematically the changes experienced so far, capturing amongst others those changes that are due to damming (the idea is introduced in [46], a very nice example is provided by [47]). (x) The historical mobility space of the rivers should be determined as a reference [48]. (xi) Geomorphic dynamic equilibrium should be assessed by integrating several elements like those already included within the River Style classification, the comparison of transport capacity amongst reaches (to identify sedimentation/incision propensity), the historical sediment budgets, the "story of the river" itself, the evidence of armoring/paving, the condition of bars and riparian vegetation, and of course modeling. (xii) An "interpretative theory", trying to find credible explanations that show relationships between causes and effects should then be conceived and described. (xiii) Suitable mathematical models should be developed (and calibrated), and then coupled, to describe all relevant processes, amongst which in particular: hydrological rainfall-runoff behavior (to capture the effects of climate change and land use changes); reservoir operation and sediment trapping; rivers hydraulic simulation coupled to sediment transport/balance and morpho-dynamics, to eventually determine the characteristics of overbank flooding as well as channel changes (more aggregated models, despite the clever ideas they incorporate (e.g., [33]), cannot provide the sought answers as they always miss part of the problem/phenomena, as for instance morphological changes of river beds).

Effects at the Border of the Mojana Region: Hydro-Morphological Risk
The whole picture of the causal factors/processes that are likely to affect the borders of the Mojana region in the next future can be summarized as follows: ‚ "From the sky" (climate change): Higher and/or more frequent flood peaks and longer and harsher droughts (deeper hydrological variability).
‚ "From the sea side" (climate change): Rising of sea level, which controls the level of the Magdalena River at its outlet at Barranquilla.
‚ "From the mountains": Reduction of sediment load to the rivers because of dams (see following Box 4.2), of the control works and of the dredging (we ignore, however, when these factors will overcome the experienced increment of solid load due to land use change occurred during the last decades).
‚ "From the rivers": Geomorphic adaptation reaction (change of the transport capacity of liquid and solid flows) because of anthropogenic interventions as well as the change in the hydrological regime due to climate change. This implies a change of longitudinal profile and slope (incision and aggradation); bankfull width; sinuosity; number and planform of channels; riparian vegetation, fish, fauna, etc.; and increment of instability (divagation with bank erosion, avulsion, incision/aggradation, etc.), eventually affecting river dynamics also within the Mojana region.
In particular, the middle reach of the Magdalena river (possibly between La Dorada and Barrancabermeja-see Figure 7-or even more) might experience an incision process (once the increment of solid load due to land use change upstream will have come to an end: when?) owing to the construction of dams upstream which will capture a significant portion of the sediment load (see Box 4.2 impact of dams), as well as owing to the extensive dredging and the control works against river dynamics envisaged for the sake of fluvial navigability (bank protections, re-sectioning) (see Figure 8).

Evolution of Cauca River
Without new dams upstream, but with the "dique marginal Nechí-Astilleros" in the Cauca River (the artificial levee on the left bank, between Colorado and La Raya settlements, see Figure 8), the riverbed is likely to suffer a significant aggradation. In general, the effect of embankments depends on the balance of two opposite phenomena: (i) if the cross section is significantly narrowed, velocity-under the same flow rate-increases and so does scour and the river incises trying to lower its slope and hence its transport capacity towards a new balance (one of such cases is described for instance in [49] for a gravel bed river); and (ii) when naturally a significant fraction of (suspended) sediment load is frequently lost by overbank flows to the floodplain (as in the Cauca River), embankments prevent this loss, so increasing sediment load; hence, sediment start accumulating within the river bed itself and aggradation occurs as the river requires a higher slope in order to increase its transport capacity and get rid of this sediment excess towards a new balance (this is a case, for instance, of the lower Po River in Italy). The Cauca River, according to the rapid River Style assessment (see Section 2.2), is likely to exhibit this latter process.
The available information does not allow us to determine whether, with the construction of the new dams foreseen upstream (see Figure 7), the process of sediment capture by reservoirs would compensate or even overwhelm the aggradation process; but some evidence speaks in favor of this hypothesis (see Figure 10).

Preliminary Estimation of the Impact of Dams on Solid Transport
We introduce here a very simplified approach by assuming that the progressive filling of the dead storage of reservoirs can be an acceptable estimator of the solid load captured and subtracted from the solid load transported by the river.
Supposing then that the dead storage (Vd) of each reservoir be approximately one-tenth of its capacity (we have no more detailed information) and supposing that all reservoirs will take T years to fill it with sediments without significant artificial emptying, it is possible to estimate the fraction  of the solid flow (QS Magda) of Magdalena River captured, as follows: where η: porosity of wet sediments; γS= g × S, with g acceleration of gravity (m/s 2 ) and S density of dry sediments (kg/m 3 ). According to what is shown in the following Table 4, the upper limit of β would result significantly high (around 50%) in the scenario considering the dead storage of all new foreseen reservoirs (including those being currently constructed: Ituango and Quimbo) plus half of that of existing reservoirs (approximately 1200 × 10 6 m 3 ) and that filling would take only 30 years with a porosity of 58% (corresponding to a bulk density of 1080 kg/m 3 ): the effect of reservoirs on the subtraction of sediments load would be extremely high (most probably exaggerated). However, the lower limit (around 8%) would still be quite significant-considering only the new reservoirs, a longer horizon of 50 years and a porosity of 85% (corresponding to a bulk density of 380 kg/m 3 ). Results more alarming than these have been found in the Mekong River [51].
An efficiency of sediment capturing around 20%-40% for existing reservoirs in the Magdalena-Cauca Rivers basin was estimated by [52]; this datum, however, is not immediately comparable with ours-because there is a whole river network in between-but speaks anyway of significant impacts.
In summary, the problem of sediments capture seems to be quite important.  [50].

Preliminary Estimation of the Impact of Dams on Solid Transport
We introduce here a very simplified approach by assuming that the progressive filling of the dead storage of reservoirs can be an acceptable estimator of the solid load captured and subtracted from the solid load transported by the river.
Supposing then that the dead storage (V d ) of each reservoir be approximately one-tenth of its capacity (we have no more detailed information) and supposing that all reservoirs will take T years to fill it with sediments without significant artificial emptying, it is possible to estimate the fraction β of the solid flow (Q S Magda ) of Magdalena River captured, as follows: According to what is shown in the following Table 4, the upper limit of β would result significantly high (around 50%) in the scenario considering the dead storage of all new foreseen reservoirs (including those being currently constructed: Ituango and Quimbo) plus half of that of existing reservoirs (approximately 1200ˆ10 6 m 3 ) and that filling would take only 30 years with a porosity of 58% (corresponding to a bulk density of 1080 kg/m 3 ): the effect of reservoirs on the subtraction of sediments load would be extremely high (most probably exaggerated). However, the lower limit (around 8%) would still be quite significant-considering only the new reservoirs, a longer horizon of 50 years and a porosity of 85% (corresponding to a bulk density of 380 kg/m 3 ). Results more alarming than these have been found in the Mekong River [51].
An efficiency of sediment capturing around 20%-40% for existing reservoirs in the Magdalena-Cauca Rivers basin was estimated by [52]; this datum, however, is not immediately comparable with ours-because there is a whole river network in between-but speaks anyway of significant impacts.
In summary, the problem of sediments capture seems to be quite important. In this latter case, the river is likely to change style (possibly downstream of Caucasia-see Figure 8, lowest part) because of the solid load deficit: we would hence witness a progressive incision and/or an increment of its sinuosity, in order to reduce its slope and hence its transport capacity and, at the same time, trying to get some material from its own banks and bed (see Figure 11). In any case, this process will also lead to a more frequent/likely de-stabilization of works (bridges, bank protections, levees, etc.) and eventually again to an (important) increment of OMR costs (Operation, Maintenance and Replacement). All these phenomena will be aggravated by climate change (more frequent and intense floods).  In this latter case, the river is likely to change style (possibly downstream of Caucasia-see Figure 8, lowest part) because of the solid load deficit: we would hence witness a progressive incision and/or an increment of its sinuosity, in order to reduce its slope and hence its transport capacity and, at the same time, trying to get some material from its own banks and bed (see Figure 11). In any case, this process will also lead to a more frequent/likely de-stabilization of works (bridges, bank protections, levees, etc.) and eventually again to an (important) increment of OMR costs (Operation, Maintenance and Replacement). All these phenomena will be aggravated by climate change (more frequent and intense floods). Figure 11. The capture of sediments by existing and future reservoirs (mainly on Cauca River, upstream of Caucasia and on Nechí River basin) will lead to a reduction of sediment load downstream; this will trigger an adaptation of the river to reduce its transport capacity, basically by lowering its slope and getting to a new river style: it probably will widen its meanders belt downstream of Caucasia, and increase its water depth owing to the lower slope. Figure 11. The capture of sediments by existing and future reservoirs (mainly on Cauca River, upstream of Caucasia and on Nechí River basin) will lead to a reduction of sediment load downstream; this will trigger an adaptation of the river to reduce its transport capacity, basically by lowering its slope and getting to a new river style: it probably will widen its meanders belt downstream of Caucasia, and increase its water depth owing to the lower slope.

Flooding Hazard Increases in the Lower Basin
The phenomenon of incision of the middle Magdalena River reach (which would take place probably in few decades, after dams construction and in parallel to dredging) might lead to increased peak flows downstream as the water retention of current floodplain would be partially inhibited, so concentrating enormous water volumes downstream and elevating water levels there, around Magangué, and-owing to backwatering-even in the zone around the San Jorge River and possibly affecting Plato, Calamar and even the very important city of Barranquilla and along the Canal del Dique towards Cartagena (see Figure 12). This phenomenon would be exacerbated by the presence of the "dique marginal" (Cauca River left bank levee) as this work is exactly meant to reduce flooding towards the North Western humedales zone (furthermore, it would definitely worsen the flood risk in stretches of the Cauca River itself downstream of the dike itself).
An increased flooding frequency/intensity of the lower Magdalena-Cauca Rivers basin will occur also-and without any possibility of countermeasures-because of the expected sea level rise which will affect hundreds of km upstream owing to the very low river slope (virtually zero or few units per thousand). In addition, the harsher flood events due to climate change will contribute to worsen the flooding hazard. The problem will also worsen because of the exacerbated subsidence process that the Mojana region will very probably experience as a result of the reduction of solid load to the floodplain caused by the "dique marginal", the interventions (dredging and re-sectioning) for increasing the conveyance of some key caños to reduce their overflows (while others are going to be closed to avoid flooding of nearby properties dedicated to cattle farming), and the embankments that will accompany the extensive system of new roads foreseen (e.g., the San Marcos-Majagual-Achí road) (see Figure 8) (an additional point refers to the effect of tectonics and orogeny, which we did not explore), while the floodplain will plunge, the riverbed elevation in lower reaches protected by levees is, on the contrary, likely to rise (as lateral sediment dispersion is impeded) and with it also the flood peak levels.  An increased flooding frequency/intensity of the lower Magdalena-Cauca Rivers basin will occur also-and without any possibility of countermeasures-because of the expected sea level rise which will affect hundreds of km upstream owing to the very low river slope (virtually zero or few units per thousand). In addition, the harsher flood events due to climate change will contribute to worsen the flooding hazard. The problem will also worsen because of the exacerbated subsidence process that the Mojana region will very probably experience as a result of the reduction of solid load to the floodplain caused by the "dique marginal", the interventions (dredging and re-sectioning) for increasing the conveyance of some key caños to reduce their overflows (while others are going to be closed to avoid flooding of nearby properties dedicated to cattle farming), and the embankments that will accompany the extensive system of new roads foreseen (e.g., the San Marcos-Majagual-Achí road) (see Figure 8) (an additional point refers to the effect of tectonics and orogeny, which we did not explore), while the floodplain will plunge, the riverbed elevation in lower reaches protected by levees is, on the contrary, likely to rise (as lateral sediment dispersion is impeded) and with it also the flood peak levels.
In summary, higher flood levels, higher elevation difference between the water level within the river and the surrounding terrain because of subsidence (where human activities lie) and river bed aggradation and flood concentration due to levees, together with an increase of the exposed value-as settlements and activities will assuredly grow towards hazard-prone areas (as the international experience demonstrates)-and together with harsher climatic events, will inevitably increase flood risk and eventually provoke more damages.

Both Drought Hazard and Vulnerability Increase
More prolonged dry periods due to climate change, together with the incision of the Magdalena River, and its reduced frequency of overflows on its floodplain, will lower the water table. The resulting lower water availability will provoke harsher water crisis in agricultural activities, increase human water supply costs and affect ecosystems.

Transformation of Wetlands
The destiny of the vast, characteristic wetlands of the Mojana region is uncertain. The "dique marginal" along the left bank of the Cauca River, together with the other interventions, will reduce water supply to the humedales zone and, what is worse, bring about a reduction of the water "pulses" which characterize so profoundly their functioning; governing fish reproduction, migration and hence fishery, amongst other aspects. However, in case the net subsidence and the flood peaks concentration downstream prove very significant, the humedales will experience, on the contrary, a progressive transformation into permanent water bodies. The sad part of it is that in both cases an impact on ecosystems will take place, particularly on the more delicate "transition components" (playones: temporary emerged bars; zápales: wet woods; etc.) together with their associated processes and environmental services [1]. This does not prevent the wetlands from drying up in an even more radical fashion than today, when the El Niño (i.e., hot ENSO) phenomenon occurs, particularly if the current tendency to close some natural "caños" which feed them (in order to widen dry lands for cattle raising (as already noted, the situation is articulated, because some "caños" are going to be modified to increase conveyance and/or navigability, while others are being closed)) and to destroy the riparian vegetation (which contributes to capture atmospheric water and reduce evaporation through shadow effect) will keep going on.

Economics: Increased Operation, Maintenance and Replacement Costs (OMR)
Any work, once it has been built, implies ongoing expenditure virtually forever-or at least until a new configuration is implemented-in order to operate, maintain and, periodically, rebuild it (depending on its technical life span). Our argument is that such costs are usually estimated on the basis of engineering criteria, which do not include the effects of such a wide infrastructural change throughout the whole river basin and its associated geomorphological consequences.
It is very difficult to estimate the increase in OMR costs due to fluvial dynamics and the impact of interventions; attempts based on historical data usually fail because: (i) It is extremely hard to get the whole list of interventions (and associated dates and costs). (ii) It is very hard to discern based on denominations/descriptions whether it was actually an OMR cost of an existing work or rather a different, new work. (iii) The behavior of the system was certainly (very) different from the one that will take place in the future once the new interventions will be implemented.
As an indicative figure, we attempted an estimation based on the official data (available online in a Excel ® file that can be consulted as Supplementary Material), which presents the investment figures for risk management in the Mojana region between 2003-2013 spent by governmental entities. By performing the most credible and prudent selection of the items presented, we obtained (in Colombian pesos with no actualization of money value) a total expense in the period of about 338 mil M$, a total OMR cost of 121 mil M$ and a net total investment C of 217 mil M$; with these data, the ratio OMR/C is about 56%, a figure which, in just 10 years, is undoubtedly an extremely high one: in a planning time horizon of 50 years, a similar OMR cost would triple the net investment cost. In addition, because of the hydro-geomorphological processes described above, this figure is very likely to become much higher than this in the future.

Conclusions
This paper is not intended to propose a particular solution; it rather makes a point on the likely adverse consequences of the course of actions currently planned. We took a delicate position according to a "preventive approach" [19][20][21][22][23] in order to provide information useful to the decision making process even if scientifically questionable, as several assumptions need to be supported by evidence yet lacking. We think this is a sensible choice, as a more scientifically prudent approach requiring more data (which would be available possibly in 10 or 20 years from now if a coordinated research effort were put into place immediately) and more scientific evidence would miss the purpose of a timely contribution to a balanced country development.
We conclude that the assessment of flood hazard recently conducted by the modeling group at IDEAM through mathematical modeling-although very refined and consistent with its hypotheses-will inevitably lead to an under-estimation of risk, a fact which, in turn, will bias the cost-benefit evaluation and, in turn, decisions.
Indeed, actual hazard will definitely be higher owing to the dynamic, geomorphic behavior of the system, very much linked to sediment transport, an aspect (as already noted by [53]) not considered in the hydrological-hydraulic modeling which performs all of its simulations on a fixed morphological context coinciding with a picture of the current system. We argued in this paper that, to the contrary, morphological evolution will conduct to a very distinct behavior with undesired increases of flood levels in vulnerable zones, increasing gradients between fluvial water level and terrain (subject to subsidence) and higher morphological risk due to more frequent and stronger destabilization action of fluvial dynamics on existing and foreseen works and infrastructures. The benefit-cost balance may seem now more positive than it actually should be because of risk under-estimation. However, another reason is that actual costs will be significantly higher simply because OMR costs of works are very likely to grossly exceed what can be estimated by ignoring the above aspect. This very fact can, by itself, drastically change the economic judgment of the whole plan.
As well, risks due to droughts would increase, while a broad array of environmental services associated with the characteristic pseudo-periodic "water pulses" would be progressively, but inevitably, lost or diminished with enormous adverse consequences.
These observations-together with socio-economic and cultural considerations not discussed here-should foster a deep re-thinking of the development model envisaged for the Mojana leading to the recognition of its very peculiar nature and to a strategy aiming at maintaining the "water pulses" peculiarity, while helping the anthropogenic system to adapt to it and take more advantage from it than was previously done. In other words, the recommendation we give is to not interfere, within the limits of possibility, with natural processes and, as such, build as few infrastructures as possible or, at least, ensure their compatibility with such processes.
In line with several other studies, like for instance [49], we can state that the geomorphic view of the consequences of development actions opens an illuminating perspective that adds important arguments for the decision process and can help avoiding choices that might dump heavy impacts on future generations.
To overcome the limitations of the qualitative approach adopted here, several information gaps should be filled at the scale of the Magdalena-Cauca Rivers basin, a task very hard to accomplish and which would require definitely several years or decades. Undertaking this exercise is certainly advisable, but not for the decisions to be made now.