Once the phenomena have been identified and a perceptual model of the mechanisms involved has been developed, an important modelling step is the definition of the relevant variables. In the models reviewed here, the hydrology variable is always related to floods in some way while the societal aspects of the coupled system may be represented by different variables such as community sensitivity [8] and flood risk awareness [17] . Fig. 3 summarises the system components used by seven coupled human-flood models from the literature. The components are hydrology, economy, technology, politics, environment and society. For each of these components, specific variables have been chosen in the studies (shown in brackets in Fig. 3 ), depending on the specific focus and the interactions modelled. For example, in most of the cases described in Section 2 of this paper, floods affect society through damage, which is an economical variable, and in some cases society affects hydrology through building levees, which is a technological variable.
Fig. 3 . System components, variables (in brackets) and feedback loops (arrows) for seven coupled human-flood models from the literature. All six system components are always shown to facilitate comparison. Positive (negative) connections are indicated with ‘+’ (‘−’). Feedback loops have been presented by the authors of the original papers or, if the information was not available, interpreted by the authors of this paper. a. Hazarika et al. [24] , b. Di Baldassarre et al. [15] , c. Grelot and Barreteau [22] , d. Grames et al. [21] , e. Liao [27] , f. Chen et al. [8] , g. Di Baldassarre et al. [17] , h. Diagram with all possible connections.
It is important to clearly specify the meaning of each variable in terms of process, spatial and temporal scales and unit. Damage, for example, can be expressed in monetary terms (Euros) and structural protection level in terms of levee height (m). In many instances in the literature, the variables represent stocks (such as capital, levee height, flood risk awareness) rather than fluxes. Choosing variables that are measurable is usually an advantage as these facilitate model parameter calibration and model testing when confronting them with measured data. Variables such as flood level are relatively easy to measure while variables representing society may be less easily measurable. Monetary variables may be more accessible than social variables such as flood awareness which require questionnaires to elicit response from people that have experienced floods [26] .
Once the variables of the system have been specified in an unambiguous way, the interactions between these variables can be identified. The arrows in Fig. 3 show loop diagrams representing the interactions of variables of models from the literature, as presented by the authors of the original papers or, if the information was not available, interpreted by the authors of this paper. The arrows represent fluxes between the variables (stocks), as most equations are balance equations. The fluxes can be positive (i.e. increasing the influenced stock) or negative (i.e. decreasing the influenced stock). When two variables influence each other through fluxes, this results in a ‘feedback loop’. In Fig. 3 a, for example increasing vulnerability increases damage in case a flood occurs and this leads people to implement adaptation strategies to reduce their vulnerability [24] . Feedback loops can also involve more than two variables. In Fig. 3 b, floods increase the awareness of flood risk, which results in a change of the dam operation to mitigate floods, which consequently affects flood discharges [15] . Also competing feedback loops may occur, like in Fig. 3 f, where two loops connect hydrology and society, one additionally including economy and the other environment [8] . If both fluxes between two variables are negative, a negative feedback loop (or dampening feedback loop) will occur, such as between economy and society in Fig. 3 c, where an increase in damage discourages people from settling in the damaged area which, in turn, decreases the potential damage [22] . Without other controls, the floodplain would eventually be abandoned. If both fluxes between two variables are positive, a positive feedback loop will occur, such as between economy and society in Fig. 3 d [21] . Without other controls this behaviour will lead to exponential economic growth.
The different strategies to flood risk management represented in the models are reflected in the diagrams in Fig. 3 . According to these models, society tries to minimise flood risk through decreasing hazard potential (e.g., through structural protection measures, Fig. 3 b,d,e,f,g) or reducing exposure (e.g., through insurance policy, Fig. 3 c), or reducing vulnerability (e.g., through individual measures, Fig. 3 a). These actions are part of conventional flood risk assessment and management but they only reflect part of the possible connections in Fig. 3 . In reality, human-flood systems are more complex and there are other connections such as public sensitivity to environmental issues (Fig. 3 f) or the influence of regular flooding on the risk awareness of society (Fig. 3 b,g). The coupled human-flood models represented in Fig. 3 , analyse the system as a whole and may lead to the unexpected phenomena discussed in Section 2 (e.g., call effect in Fig. 3 c, levee effect in Fig. 3 e,g, adaptation effect in Fig. 3 a, change in values in Fig. 3 b,f). Ideally, scientists or flood risk managers should examine all possible connections between the system components (Fig. 3 h) and implement the dominant ones for the real world situation in question. For example, when local managers are considering building a dam for flood protection and hydro-power generation, the hydrological, economic and technological aspects will be on their minds. However, they should also take into account the environment, society and politics. Otherwise, unexpected consequences may occur in the near or far future.
4.Models of feedbacks
Models that represent the dynamic phenomena and feedbacks of the previous sections may be classified, according to their complexity, into (a) system-of-systems models, and (b) stylised models [35] .
System-of-Systems models (or coupled component models, [1] ) consist of coupled, detailed models of the hydrologic and societal parts of the system. Typically, well established and tested models for each discipline are used, which enhances the credibility of the overall model. They are usually spatially explicit and can be tailored to specific, real-world management problems. Falter et al. [20] , for example, coupled a set of models representing the complete flood risk chain: a weather generator, a hydrological model, a coupled 1D-2D hydrodynamic model and a flood loss model. While Falter et al. [20] did not include the impact of flood effects on societal processes, a more detailed model of social behaviour could be implemented as suggested by O’Connell and O’Donnell [30] . System-of-Systems models are amenable to practical problems, because of the local detail included. However, representing feedbacks may not be straightforward, because the disciplinary models may not have been designed for being integrated with other models. Also, understanding the overall system behaviour may be difficult because of the system complexity.
Stylised models, on the other hand, are simple/parsimonious models that directly couple the hydrologic and societal parts of the system. Feedbacks are built into the models during their development, e.g., starting from the loops in Fig. 3 , and the models are used to understand how the feedbacks lead to the phenomena described in Section 2 . Because of their relative simplicity, stylised models usually allow a clear understanding of the system behaviour but are not ideal for describing real-world management problems in specific cases. They are used for generalising the understanding of phenomena beyond individual case studies. These stylised models could be used to inform the development of more complex system-of-systems models. The outcomes of stylised models can inform the modeler which components should be represented in system-of-systems models and how these components should be connected. The system-of-systems models should be able to reproduce the same kind of feedbacks as the stylized models but can include more detail about the specific characteristics of the case study.
Most of the coupled human-flood models proposed in the literature are based on systems of a few coupled ordinary differential equations (termed ‘system dynamic’ models by [1] ). For example, the model of Di Baldassarre et al. [17] and Viglione et al. [42] represents the evolution of four variables by four coupled non-linear differential equations. The equations are coupled in the sense that the change of one variable with time depends on some of the other variables. They are non-linear mainly due to threshold processes such as when a flood overtops a levee or the discharge exceeds a certain value like in the Kissimmee system of Chen et al. [8] . In system dynamic models, most differential equations are conceptual representations of the lumped system behaviour, i.e. one variable such as flood awareness is an aggregated, representative value of the entire domain.
In contrast, agent based models start from the behaviour of individuals (or small groups of people) termed agents. The modeller describes the interactions between agents through decision rules, and these interactions alter the agents’ state [1] . In the model of Grelot and Barreteau [22] , households decide to move into a floodplain based on the attractiveness of being close to the river and living in a populated area, they are discouraged if other houses in the area are dilapidated, and a percentage of the population leaves the floodplain in every time step. There is also a regulator (also an agent) determining the insurance fee based on the number of households and the potential flood damage. Even with very simple rules at the agent level, these models are able to simulate complex macro-scale consequences (‘emergent patterns’, [36] because of the number and heterogeneity of the agents. In Grelot and Barreteau [22] , the introduction of an insurance system generates a call effect resulting in an increase of the risk taken by the households. While the upscaling through the interaction of agents is attractive, the output of agent based models and the aggregate system behaviour is sometimes difficult to understand because the connection between the variables is less clear at the macro-scale than for system dynamic models.
Both system dynamic and agent based models can be completely deterministic or may allow for stochasticity. Most of the coupled human-flood models reviewed in this paper are internally deterministic with external forcing that may be stochastic (i.e., flood occurrence and magnitude). However, stochasticity could be incorporated into the models. In system dynamic models, noise terms could be added to the ordinary differential equations making them stochastic, similar to models of eco-hydrology [33] . Examples of stochastic agent based models include the model of Dawson et al. [12] where the daily routine of individuals may or may not be interrupted (with a given probability) when they become aware of a flood, and the model of Grelot and Barreteau [22] , where the population leaving the floodplain is selected randomly. The advantage of stochastic over deterministic models is their ability to represent the uncertainty associated with the connections between variables (i.e., with the feedbacks) as well as ‘random’ modifications of state variables not determined by the other variables in the model. The uncertainty in the outputs of the models then may be evaluated by uncertainty propagation methods.
Almost all coupled human-flood models published so far in the literature are descriptive models, in that they aim at describing the human – flood interactions as they are observed (or could be observed). The purpose is to understand the interactions, rather than searching for optimal policies/investments to minimise flood risk and/or maximise economic income. An exception is Grames et al. [21] who used a model in a prescriptive way. In this model, an economic decision maker interacts with the hydrological system by choosing how much to invest in flood defence in order to maximise the total economic income. O’Connell and O’Donnell [30] suggest that prescriptive coupled human-flood models may be more useful in determining appropriate flood investments than the traditional cost-benefit analysis, in particular if flood rich and flood poor periods occur.
All models reviewed in this paper have been developed starting from concepts derived from personal experience/observation of how the system works and comparisons with data have been performed in a qualitative way. For example, Ciullo et al. [10] compared the outputs of the model of Di Baldassarre et al. [18] with data of population density, flood losses and levee heights in the city of Rome and in Bangladesh. In a similar way, Chen et al. [8] used flood series and data about wetland area to assess whether the model outputs are realistic. The next step for coupled human-flood modelling is the development of more formal methodologies to test and calibrate coupled human-flood models, perhaps based on Bayesian methods [43,41] . However, this poses a challenge, because the usual time-series fitting of hydrological models may not be the best option, since phenomena and feedbacks are of interest. New ways of testing how well the models reproduce observed patterns need to be developed, such as matching trajectories in the phase space of the model variables.
Starting from concepts and testing model outputs against data is one way of model development (i.e., a hypothesis testing approach). Alternatively, investigative data-based modelling could be performed starting from the data, without assuming a priori model structures. Data analyses that attempt to find relationships between the system components from data correlations (e.g., [16,11] ) could be starting points for such models and may help bring coupled human-flood modelling to the next level. However, the lack of long-term data on individual system components poses a challenge for this type of methods. Comparative analyses could be of help, in the spirit of the PUB initiative [3] , where the analysis in time is replaced by analyses in space [16] , although the validity of this assumption in the context of co-evolutionary systems needs to be tested [31] . Comparative analyses may also become more wide spread thanks to the availability of new types of data [7] . Especially for coupled human-water systems this is more difficult than for hydrologic systems alone [35] , but has the potential of boosting new insights and the development of new methodologies.
5.Implications for flood risk management
The European Flood Directive [19] explicitly requires the establishment of flood risk management plans that account for the complexity of the system as a whole. The underlying framework is termed ‘Integrated Flood Risk Management’ (IFRM), which implies integration between sectors (e.g., water management, transport, regional planning and tourism) and between upstream and downstream reaches in a river basin (i.e., the so called ‘solidarity principle’) [6] . IFRM seeks to reduce hazard and exposure at a regional scale by implementing a combination of various measures, such as structural (e.g., levees, polders) and non-structural (e.g., insurances, land-use zoning) measures. IFRM also seeks to reduce vulnerability through increasing people’s awareness of flood risk [39] . Implementation of the Flood Directive is well underway in most European countries. For the specific case of Vienna, as discussed in Section 1 of this paper, the main management measure is the relief channel, supported by non-structural measures such as warning, evacuation plans, and insurance.
Traditionally, decisions on what flood management measures to implement and to what level to protect have been highly political, since citizens may be directly affected and management measures can be very expensive. More recently, cost-benefit analyses are being performed, as stipulated in the Flood Directive. Typically, cost-benefit analyses compare a small number of scenarios to ascertain which scenarios are economically more favourable. However such scenarios do not usually take into account that the measures implemented will have an effect on the future flood risk because of the feedback loops in the system. Therefore the optimum combinations of measures chosen now may not be robust against the future evolution of the risk [5] . An alternative consists of adopting measures that are not optimal for a single future situation but are reasonable for a wide variety of different situations (termed robust, ‘low-regret’, or resilient approaches) [46] . For example, allowing floods to inundate into flood plains may be a more robust strategy than levees that will be overtopped once the protection level is exceeded. In practice, there is often a combination of a number of different measures and it is difficult to understand their interplay. A dynamic approach with coupled human-flood models can help explore the possible future situations that may arise because of the interplay between the system components.
Consider the flood damage simulations for a hypothetical city in Fig. 4 . The simulations show the effect of two alternative flood risk management approaches over a period of two centuries using the model of Di Baldassarre et al. [17] . In the top right panel of Fig. 4 the options are to build close or far away from the river but there is no option to build levees. In the bottom right panel, the flood management options also involve the construction of levees. In this second case, the hazard has been decreased substantially after the construction of levees around year 50, however exposure has increased because of the levee effect and the consequences of the events around year 150 are catastrophic. In the top panel, instead, more floods have produced damage in the years before 150 and therefore exposure has not increased as much, which results in a lower damage for the big events around year 150. The levee effect is triggered by the flood awareness of the community which may be low during a period during which few floods occur [23] . When a flood rich period ensues (around year 130 in Fig. 4 ) the damage may be much larger than expected. Simulations such as those in Fig. 4 explore what could possibly happen and may therefore be used to advise decision making in the implementation of IFRM. They can show the consequences of a particular flood management decision and may help in comparing the measures that are being considered or show a need for additional measures. Simulations like in Fig. 4 could have informed flood risk management in the case of Vienna: additional measures (e.g. land use planning or flood zoning) could have been taken to reduce the increase in exposure that occurred as a result of the implemented structural measures.
Fig. 4 . Two scenarios of flood damage for a hypothetical city. Top right: Damage resulting from flood management options involving the choice of building close to or far away from a flood prone river but no levees. Bottom right: Damage resulting from flood management options also involving the construction of levees. For comparison, the top left panel shows the input flood water levels (for both cases) and the modelled flood risk awareness of the community for the case involving levees. Light blue areas indicate flood-rich periods, white areas flood-poor periods. Results from the model of Di Baldassarre et al. (2013b) and Viglione et al. (2014).
Coupled human-flood modelling of many different systems around the world, such as those reviewed in this paper, may also be useful to inform IFRM at a particular location because they may provide a wider overview of the possible phenomena that have not occurred yet, but may occur. For example, for the Danube case, the levee effect has already occurred but other phenomena not yet observed such as the adaptation effect or a shift from flood risk concern to environmental concern may occur in the future. Stylised models, such as those reviewed here, could be used to investigate the effects of different feedback loops. The outcome of these models could then be the basis for building comprehensive and spatially explicit system-of-systems models to investigate the effects of alternative flood management decisions in more detail and support IFRM through a dynamic framework.
Acknowledgements
We would like to acknowledge funding from the Marie Sklodowska-Curie Innovative Training Network “A Large-Scale Systems Approach to Flood Risk Assessment and Management – SYSTEM-RISK“ (Grant Agreement number: 676027); the ERC FloodChange project (ERC Advanced Grant FP7-IDEAS-ERC-AG-PE10 291152); and the FWF Vienna Doctoral Programme on Water Resource Systems (W1219-N22).
We would like to thank two anonymous reviewers and Bruno Merz for their useful comments to the original version of the paper.