Community Perception and Communication of Volcanic Risk from the Cotopaxi Volcano in Latacunga, Ecuador

The inhabitants of Latacunga living in the surrounding of the Cotopaxi volcano (Ecuador) are exposed to several hazards and related disasters. After the last 2015 volcanic eruption, it became evident once again how important it is for the exposed population to understand their own social, physical, and systemic vulnerability. Effective risk communication is essential before the occurrence of a volcanic crisis. This study integrates quantitative risk and semi-quantitative social risk perceptions, aiming for risk-informed communities. We present the use of the RIESGOS demonstrator for interactive exploration and visualisation of risk scenarios. The development of this demonstrator through an iterative process with the local experts and potential end-users increases both the quality of the technical tool as well as its practical applicability. Moreover, the community risk perception in a focused area was investigated through online and field surveys. Geo-located interviews are used to map the social perception of volcanic risk factors. Scenario-based outcomes from quantitative risk assessment obtained by the RIESGOS demonstrator are compared with the semi-quantitative risk perceptions. We have found that further efforts are required to provide the exposed communities with a better understanding of the concepts of hazard scenario and intensity.


Introduction
An active volcanic environment is prone to produce cascading and compound natural hazards. Cascading hazards comprise a primary hazard triggering a secondary one [1], whilst compound hazards refer to events (not necessary interdependent) events whose

Framework and Objectives
Volcanic eruptions pose an enormous risk to Ecuador, because most of the exposed human settlements in the central and northern highlands are situated less than 25 km from an active volcano. Cities previously affected by volcanic eruptions include Quito, Latacunga, Salcedo, Cayambe, Ibarra-Otavalo, Ambato, Riobamba, and Baños [55]. Lahars have been among the deadliest volcanic hazards, but the emission of volcanic ash has been more frequent in the Ecuadorian Andes [56]. Ash falls do not only have direct consequences on the inhabitants' health and on the exposed infrastructure, but also on agriculture and animal husbandry, which is particularly important for the rural communities in Ecuador. Ash falls have hit the rural communities settled in the vicinity of the most active Ecuadorian volcanoes (i.e., Tungurahua, Reventador, Sangay, and Cotopaxi). Moreover, poverty, Sustainability 2021, 13, 1714 4 of 26 marginality, and high inequality of the exposed communities coexist with their physical and systemic vulnerabilities [57].

Description of the Study Area
The Cotopaxi volcano is an active stratovolcano (5897 m.a.s.l) located in the Cordillera Real of the Ecuadorian Andes ( Figure 1) and is covered by an extensive, but diminishing glacier cap. Cotopaxi is one of the most dangerous volcanoes worldwide [58] with average recurrence intervals for eruptions between 117-147 years [59]. It can produce syneruptive lahars triggered by explosive eruptions, which can travel hundreds of kilometres [60]. Three drainage systems originate on Cotopaxi (Figure 2), which have all been inundated by lahars in prehistoric times [61]. However, only the northern and southern drainage systems are densely populated: The largest urban agglomeration encountered by the northern system is "El Valle de Los Chillos" (with about 400,000 inhabitants) in the vicinity of southern Quito, whilst the southern drainage system encounters the Latacunga canton (with about 300,000 inhabitants). The last major eruption of the Cotopaxi volcano in the historical records occurred in 1877. It induced syneruptive lahars that severely affected the proximal rural communities [62], with more than 1000 deaths registered, and caused a severe economic crisis [63]. If a similar scenario occurred nowadays, the social and economic consequences would be far more catastrophic due to the high population density and the central importance of Latacunga and the Cotopaxi region for the economic development of the country [58]. Ecuador. Ash falls have hit the rural communities settled in the vicinity of the most active Ecuadorian volcanoes (i.e., Tungurahua, Reventador, Sangay, and Cotopaxi). Moreover, poverty, marginality, and high inequality of the exposed communities coexist with their physical and systemic vulnerabilities [57].

Description of the Study Area
The Cotopaxi volcano is an active stratovolcano (5897 m.a.s.l) located in the Cordillera Real of the Ecuadorian Andes ( Figure 1) and is covered by an extensive, but diminishing glacier cap. Cotopaxi is one of the most dangerous volcanoes worldwide [58] with average recurrence intervals for eruptions between 117-147 years [59]. It can produce syneruptive lahars triggered by explosive eruptions, which can travel hundreds of kilometres [60]. Three drainage systems originate on Cotopaxi (Figure 2), which have all been inundated by lahars in prehistoric times [61]. However, only the northern and southern drainage systems are densely populated: The largest urban agglomeration encountered by the northern system is "El Valle de Los Chillos" (with about 400,000 inhabitants) in the vicinity of southern Quito, whilst the southern drainage system encounters the Latacunga canton (with about 300,000 inhabitants). The last major eruption of the Cotopaxi volcano in the historical records occurred in 1877. It induced syneruptive lahars that severely affected the proximal rural communities [62], with more than 1000 deaths registered, and caused a severe economic crisis [63]. If a similar scenario occurred nowadays, the social and economic consequences would be far more catastrophic due to the high population density and the central importance of Latacunga and the Cotopaxi region for the economic development of the country [58]. Latacunga is the largest city of the Latacunga canton (second Ecuadorian administrative division) and it is the capital of the Cotopaxi province (first division). It is located at a 14 km distance from the Cotopaxi volcano. For the year 2020, and based on the population projections of the National Institute of Statistics and Censuses [65], the city has an inferred population of approximately 205,600 inhabitants, with a major rural composition (59.8%). The last peak of volcanic activity of the Cotopaxi volcano occurred in mid-April Latacunga is the largest city of the Latacunga canton (second Ecuadorian administrative division) and it is the capital of the Cotopaxi province (first division). It is located at a 14 km distance from the Cotopaxi volcano. For the year 2020, and based on the population projections of the National Institute of Statistics and Censuses [65], the city has an inferred population of approximately 205,600 inhabitants, with a major rural composition (59.8%). The last peak of volcanic activity of the Cotopaxi volcano occurred in mid-April 2015 and lead to a crisis in risk management in Latacunga and neighbouring municipalities [66]. Firstly, an increase in the seismic activity of the volcano was accompanied by the emission of sulphur dioxide and ash fall for some weeks [56]. Subsequently, authorities and local press communicated to the inhabitants of the communities in the vicinity of the Cotopaxi volcano that it was necessary to evacuate their homes promptly due to the imminent occurrence of lahars [66]. This generated social chaos due to the ignorance of the evacuation routes, the uncontrolled behaviour of the citizens (due to generalised fear of looting), as well as a very low level of trust in government representatives [67]. Eventually, the 2015 activity never surpassed a magnitude VEI 2 and no large syneruptive lahar flows occurred [68]. The lesson learned from this experience was the need for adequate evacuation protocols and local authorities with an understanding of the complexity of the risk in the area. Moreover, it was realized how important it is for citizens to understand their own social, physical, and systemic vulnerability [68]. 2015 and lead to a crisis in risk management in Latacunga and neighbouring municipalities [66]. Firstly, an increase in the seismic activity of the volcano was accompanied by the emission of sulphur dioxide and ash fall for some weeks [56]. Subsequently, authorities and local press communicated to the inhabitants of the communities in the vicinity of the Cotopaxi volcano that it was necessary to evacuate their homes promptly due to the imminent occurrence of lahars [66]. This generated social chaos due to the ignorance of the evacuation routes, the uncontrolled behaviour of the citizens (due to generalised fear of looting), as well as a very low level of trust in government representatives [67]. Eventually, the 2015 activity never surpassed a magnitude VEI 2 and no large syneruptive lahar flows occurred [68]. The lesson learned from this experience was the need for adequate evacuation protocols and local authorities with an understanding of the complexity of the risk in the area. Moreover, it was realized how important it is for citizens to understand their own social, physical, and systemic vulnerability [68].  Latacunga is settled on ancient and recent geological materials formed by volcanic material. Some of the most representative and better-exposed stratigraphic formations Sustainability 2021, 13, 1714 6 of 26 of ancient ashes and lahar deposits originated from the previous volcanic activity of the Cotopaxi volcano were visited ( Figure 3) with the guidance of experts from the Geophysical Institute of the National Polytechnic School (IG-EPN (Instituto Geofísico de la Escuela Politécnica Nacional (Quito, Ecuador))) and the Decentralized Autonomous Government of the Province of Cotopaxi (GADPC (GADPC, Gobierno Autónomo Descentralizado Provincial de Cotopaxi, Latacunga, Ecuador)). Some of these deposits are from pre-historical times, whilst the shallower ones date from the 1877 event that destroyed Latacunga [70]. Official maps of the Geological and Energy Research Institute (IIGE) and IG-EPN [71] were used during the field reconnaissance. This field trip was relevant to visualize the geological characteristics of the study area, as well as to strengthen the cooperation and idea exchanges with the local experts. Latacunga is settled on ancient and recent geological materials formed by volcanic material. Some of the most representative and better-exposed stratigraphic formations of ancient ashes and lahar deposits originated from the previous volcanic activity of the Cotopaxi volcano were visited ( Figure 3) with the guidance of experts from the Geophysical Institute of the National Polytechnic School (IG-EPN (Instituto Geofísico de la Escuela Politécnica Nacional (Quito, Ecuador))) and the Decentralized Autonomous Government of the Province of Cotopaxi (GADPC (GADPC, Gobierno Autónomo Descentralizado Provincial de Cotopaxi, Latacunga, Ecuador)). Some of these deposits are from pre-historical times, whilst the shallower ones date from the 1877 event that destroyed Latacunga [70]. Official maps of the Geological and Energy Research Institute (IIGE) and IG-EPN [71] were used during the field reconnaissance. This field trip was relevant to visualize the geological characteristics of the study area, as well as to strengthen the cooperation and idea exchanges with the local experts. Latacunga is not only exposed to the natural hazards imposed by the Cotopaxi volcano, but also to other geodynamic (e.g., landslides and earthquakes) and hydro-climatologic hazards (e.g., frosts and droughts). As reported in [72], there has been an intensification in the variability of precipitations, droughts, and frosts in Latacunga. This has been evidenced in the period between the years 1981-2014, during which the average air temperature has increased about 0.8°C. These ongoing phenomena related to climate change have generated negative consequences mainly in the rural area and in agriculture areas [72].

Objectives
The understanding of disaster risk based on the independent investigation of their dimensions, hazards, exposure, and vulnerability, guided by a multi-hazard risk approach with risk-informed decision-makers, is the key advice of the Sendai Framework for Disaster Risk reduction (2015-2030) [73]. Having in mind the aforementioned limitations on volcanic risk assessment as well as the lack of exploration tools for risk communication, two initiatives, namely the programme "Sustainable Intermediate Cities-CIS" and the research project "Multi-Risk Analysis and Information System Components for the Andes Region-RIESGOS" have been working in Latacunga, Ecuador, with the aim of increasing awareness and preparedness and enhancing the coping capacities of the communities exposed to the Cotopaxi volcano. The particular objectives of this integrative study are: Latacunga is not only exposed to the natural hazards imposed by the Cotopaxi volcano, but also to other geodynamic (e.g., landslides and earthquakes) and hydro-climatologic hazards (e.g., frosts and droughts). As reported in [72], there has been an intensification in the variability of precipitations, droughts, and frosts in Latacunga. This has been evidenced in the period between the years 1981-2014, during which the average air temperature has increased about 0.8 • C. These ongoing phenomena related to climate change have generated negative consequences mainly in the rural area and in agriculture areas [72].

Objectives
The understanding of disaster risk based on the independent investigation of their dimensions, hazards, exposure, and vulnerability, guided by a multi-hazard risk approach with risk-informed decision-makers, is the key advice of the Sendai Framework for Disaster Risk reduction (2015-2030) [73]. Having in mind the aforementioned limitations on volcanic risk assessment as well as the lack of exploration tools for risk communication, two initiatives, namely the programme "Sustainable Intermediate Cities-CIS" and the research project "Multi-Risk Analysis and Information System Components for the Andes Region-RIESGOS" have been working in Latacunga, Ecuador, with the aim of increasing awareness and preparedness and enhancing the coping capacities of the communities exposed to the Cotopaxi volcano. The particular objectives of this integrative study are:

1.
Presenting a comprehensive risk communication process, from scenario-based volcanic risk analysis along with active participation of the exposed communities, while also investigating the risk perception of the exposed communities.

2.
Providing a recent measurement of spatially distributed risk perception in the Cotopaxi area (results from the CIS questionnaire)

3.
Testing the applicability of the RIESGOS demonstrator, a decentralized web-service architecture that allows for integrating local expert knowledge and locally designed models in a scenario-based multi-risk analysis, for the purpose of interactive communication 4.
Merging the results of 2 and 3 to investigate how well the simulated quantitative risk matches the subjectively perceived risk in a common area.

Materials and Methods
An integrative framework between scientific approaches and risk communication practices with the exposed society has been set up in Latacunga (Ecuador) by two different initiatives, (1) the CIS (Sustainable Intermediate Cities) programme and (2) the RIESGOS project (Multi-risk analysis and information system components for the Andes region).
3.1. The CIS Programme: The Creation of a Local Laboratory to Evaluate the Social Perception of Risk and Resilience "The Latacunga Laboratory: Risk management, resilience, and adaptation to climate change (Laboratorio Urbano de Latacunga: Gestión de riesgos, resiliencia y adaptación al cambio climático)" has been created within the CIS programme, as part of the joint initiatives of GIZ ("Deutsche Gesellschaft für Internationale Zusammenarbeit") and Grupo FAR (Ecuadorian NGO (https://grupofaro.org/ (accessed on 1 January 2021)). The creation of so-called resilience observatories for exposed communities to natural hazards is a relatively new trend [45,74]. Similarly, the Latacunga Laboratory seeks to contribute to the risk management to natural hazards that are likely to occur in the territory, while aiming to contribute in the long-term to the development of the city embracing its urban-rural ties. With that goal, initial contributions related to social risk perceptions have been documented in [75] as a joint effort between the Latacunga Laboratory, the local government, academic institutions, and local actors. We conducted a survey by means of a custom-designed questionnaire, a fundamental tool for acquiring information on public knowledge of the community [76]. It is composed of a series of multiple-choice questions in Spanish. The survey was carried out in the field and online to collect data about the individual knowledge, attitudes, and risk perceptions of the inhabitants of Latacunga. The online survey was promoted on social media and was available on the official website of the CIS Latacunga Laboratory (https://latacungaresiliente.com/ (accessed on 1 January 2021)) for a month. In the meantime, the field survey was carried out only in the urban agglomeration of Latacunga. The collected data is used for two main objectives: (1) as input to perform the semiquantitative method proposed in [77] that ranks the social perception of volcanic risk factors (i.e., hazard recurrence, exposure, vulnerability, and resilience) among other natural hazards likely to occur in the study area (i.e., earthquakes, drought, frost, floods, and landslides), and (2) to map the spatial distribution of volcano-related risk perception into comprehensive categories (i.e., easily understandable by the exposed communities).
A design of the field surveying site was carried out. According to the last official census available [65] and population projections by the survey elaboration date (September 2019), 50,442 inhabitants over the age of 18 years were considered as qualified informants. In order to get a statistically representative sample, a confidence level of 95% and a margin of error of 5% were selected. On this basis, we estimated that a sample for the field surveys not smaller than 380 inhabitants had to be selected. Considering 10% additional surveys, a final sample size of 420 people was chosen. The population density (Figure 4a) was used to constrain the spatial distribution of the field surveys within the urban blocks with a residential occupancy (Figure 4b). The surveys were carried out by 55 students of the ISTC (Instituto Superior Tecnológico Cotopaxi) in September 2019. surveys not smaller than 380 inhabitants had to be selected. Considering 10% additional surveys, a final sample size of 420 people was chosen. The population density (Figure 4a) was used to constrain the spatial distribution of the field surveys within the urban blocks with a residential occupancy (Figure 4b). The surveys were carried out by 55 students of the ISTC (Instituto Superior Tecnológico Cotopaxi) in September 2019. It is worth mentioning that, on the one hand, some drawbacks have been found when the community perception of exposure, vulnerability, and resilience are independently addressed for large-scale studies [78,79]. On the other hand, there have been also reported benefits of this separation for mapping the social risk perception to natural hazards (e.g., [80,81]) when bottom-up approaches are carried out. Therefore, we have decided to independently investigate the social perceptions towards these components through separated questions. The Likert scale is used in this context to obtain a quantifiable level of perception of each risk factor. An integer numerical score (1, 2, or 3) is assigned to every possible answer. Although the passage from a qualitative perception to an index can be questioned, several recent studies have shown the usefulness of the Likert scale [82][83][84][85][86]. Notably, in [87], it was found to provide a good compromise between the quality of the information collected and the accessibility to respondents, while the bias in responses decreases, and there is consistency across different measurements and research domains of disaster risk reduction.
Subsequently, the average is computed for every question to obtain the perception of every component. These values are inputs to the computation of the risk perception preindex through the use of Equation (1), where P stands for "perception". An example subset of the questions is presented in Table A1 (Appendix A). The questions and answers were validated by local risk management experts from the Association of Risk Management Professionals of Ecuador (Asociación de Profesionales de Gestión de Riesgos de Ecuador, APGR). It is worth mentioning that, on the one hand, some drawbacks have been found when the community perception of exposure, vulnerability, and resilience are independently addressed for large-scale studies [78,79]. On the other hand, there have been also reported benefits of this separation for mapping the social risk perception to natural hazards (e.g., [80,81]) when bottom-up approaches are carried out. Therefore, we have decided to independently investigate the social perceptions towards these components through separated questions. The Likert scale is used in this context to obtain a quantifiable level of perception of each risk factor. An integer numerical score (1, 2, or 3) is assigned to every possible answer. Although the passage from a qualitative perception to an index can be questioned, several recent studies have shown the usefulness of the Likert scale [82][83][84][85][86]. Notably, in [87], it was found to provide a good compromise between the quality of the information collected and the accessibility to respondents, while the bias in responses decreases, and there is consistency across different measurements and research domains of disaster risk reduction.
Subsequently, the average is computed for every question to obtain the perception of every component. These values are inputs to the computation of the risk perception pre-index through the use of Equation (1), where P stands for "perception". An example subset of the questions is presented in Table A1 (Appendix A). The questions and answers were validated by local risk management experts from the Association of Risk Management Professionals of Ecuador (Asociación de Profesionales de Gestión de Riesgos de Ecuador, APGR).
The numerator of Equation (1) can have a maximum possible value of 27, whilst the minimum for the resilience term in the denominator is 1. Therefore the maximum risk perception value that this method admits is 27. The values in the range from 1-27 form a Sustainability 2021, 13, 1714 9 of 26 "pre-index". To obtain a more comprehensive numerical value, a "reduced index" in the 0-3 range is obtained through the application of Equation (2).
The relations between the "pre-index" and the "reduced index" is shown in Figure 5. For mapping purposes an "equal interval" classification for the reduced index scale is introduced with five classes of length 0.6 for finally presenting the spatialized perception of every risk factor in a compressive manner to the community. The calculated results for every answered question at each survey location (Figure 4b) are used to map the spatial distribution of the perception of hazard recurrence, exposure, vulnerability and resilience, and the risk index (computed with Equation (1). Subsequently, they were interpolated through the use of the ordinary kriging geostatistical algorithm [88].
The numerator of Equation (1) can have a maximum possible value of 27, whilst the minimum for the resilience term in the denominator is 1. Therefore the maximum risk perception value that this method admits is 27. The values in the range from 1-27 form a "pre-index". To obtain a more comprehensive numerical value, a "reduced index" in the 0-3 range is obtained through the application of Equation (2).
The relations between the "pre-index" and the "reduced index" is shown in Figure 5. For mapping purposes an "equal interval" classification for the reduced index scale is introduced with five classes of length 0.6 for finally presenting the spatialized perception of every risk factor in a compressive manner to the community. The calculated results for every answered question at each survey location (Figure 4b) are used to map the spatial distribution of the perception of hazard recurrence, exposure, vulnerability and resilience, and the risk index (computed with Equation (1). Subsequently, they were interpolated through the use of the ordinary kriging geostatistical algorithm [88].

The RIESGOS Project: Iterative Simulation Improvement and Enhanced Communication
The idea of constructing a web-tool, the RIESGOS demonstrator, as a decentralised and intraoperative environment for the exploration of the consequences in Latacunga from different volcanic hazard scenarios was proposed to the local stakeholders who participated in four participative workshops. Two of them were held in Latacunga (7 December 2018; 25 November 2019) in the headquarters of GADPC (Decentralized Autonomous Government of the Cotopaxi Province) and two workshops took place in Quito on 11 December 2018, and on 27 November 2019, respectively. The participants ranged from research partners, representatives of the rural municipalities (parishes) of the Cotopaxi province, public authorities, environment secretaries, actors of the civil society such as local representatives of agriculture associations, and urban and rural leaders. Similarly, as recently presented in [46], the workshops were used as a means to implement a user-centred iterative approach, seeking a continuous redesign of the RIESGOS demonstrator that has been guided by the needs of potential users and practical applicability. This has been ensured by a comprehensive analysis of user requirements (e.g., open-source, userfriendly graphical user interface and transferability).

The RIESGOS Demonstrator Tool for Quantitative Multi-Risk Analysis
The iteratively constructed RIESGOS demonstrator for a multi-risk information system is based on a modular and scalable concept in which the different hazards, the related exposure models, and vulnerability schemas are each represented by one individual web service. These independent and distributed web-services (managed and maintained by individual research institutions) are based on the quantitative methodologies developed within the RIESGOS framework for multi-risk analysis (i.e., [69,[89][90][91][92][93][94]

The RIESGOS Project: Iterative Simulation Improvement and Enhanced Communication
The idea of constructing a web-tool, the RIESGOS demonstrator, as a decentralised and intraoperative environment for the exploration of the consequences in Latacunga from different volcanic hazard scenarios was proposed to the local stakeholders who participated in four participative workshops. Two of them were held in Latacunga (7 December 2018; 25 November 2019) in the headquarters of GADPC (Decentralized Autonomous Government of the Cotopaxi Province) and two workshops took place in Quito on 11 December 2018, and on 27 November 2019, respectively. The participants ranged from research partners, representatives of the rural municipalities (parishes) of the Cotopaxi province, public authorities, environment secretaries, actors of the civil society such as local representatives of agriculture associations, and urban and rural leaders. Similarly, as recently presented in [46], the workshops were used as a means to implement a user-centred iterative approach, seeking a continuous redesign of the RIESGOS demonstrator that has been guided by the needs of potential users and practical applicability. This has been ensured by a comprehensive analysis of user requirements (e.g., open-source, user-friendly graphical user interface and transferability).

The RIESGOS Demonstrator Tool for Quantitative Multi-Risk Analysis
The iteratively constructed RIESGOS demonstrator for a multi-risk information system is based on a modular and scalable concept in which the different hazards, the related exposure models, and vulnerability schemas are each represented by one individual web service. These independent and distributed web-services (managed and maintained by individual research institutions) are based on the quantitative methodologies developed within the RIESGOS framework for multi-risk analysis (i.e., [69,[89][90][91][92][93][94]). Therefore, their integration into the RIESGOS demonstrator simulates the multi-risk environment of Latacunga. This modular approach offers the possibility to integrate different web services into already existing system environments.
Currently, the graphical user interface of the demonstrator can be accessed from a web browser only by users with special rights. The main screen of the graphical user interface is divided into three main display areas: the central map window, the configuration wizard for the control of each web service to the left, and the results panel to the right (e.g., see Figure 6). The code of the graphical user interface (RIESGOS frontend) is openly published on GitHub (https://github.com/riesgos/dlr-riesgos-frontend (accessed on 1 January 2021)). The use of standardized web services such as geospatial web services defined by the Open Geospatial Consortium (OGC) allows users accessing open and flexible multi-risk information and data products. Web-services and exposed data resources can be accessed using a variety of means from a simple command-line tool, over a web browser, to existing graphical user interfaces of public authorities and companies, which are equipped with a map user. OGC web services allow all kinds of geospatial functionality out-of-the-box including data access, data display, styling, and processing. Web services can easily be integrated into existing clients. The providers of web services define their products, display options, and configuration items. More details of this integrative process are reported in [94]. Through the clear separation in competencies between web services and user-interface, modularity and scalability are increased.
integration into the RIESGOS demonstrator simulates the multi-risk environment of Latacunga. This modular approach offers the possibility to integrate different web services into already existing system environments.
Currently, the graphical user interface of the demonstrator can be accessed from a web browser only by users with special rights. The main screen of the graphical user interface is divided into three main display areas: the central map window, the configuration wizard for the control of each web service to the left, and the results panel to the right (e.g., see Figure 6). The code of the graphical user interface (RIESGOS frontend) is openly published on GitHub (https://github.com/riesgos/dlr-riesgos-frontend). The use of standardized web services such as geospatial web services defined by the Open Geospatial Consortium (OGC) allows users accessing open and flexible multi-risk information and data products. Web-services and exposed data resources can be accessed using a variety of means from a simple command-line tool, over a web browser, to existing graphical user interfaces of public authorities and companies, which are equipped with a map user. OGC web services allow all kinds of geospatial functionality out-of-the-box including data access, data display, styling, and processing. Web services can easily be integrated into existing clients. The providers of web services define their products, display options, and configuration items. More details of this integrative process are reported in [94]. Through the clear separation in competencies between web services and user-interface, modularity and scalability are increased. Figure 6. Example of the graphical representation of loss distribution due to ash fall scenario in the RIESGOS demonstrator (as of December 2020) from a previously selected VEI. Reddish and greenish aggregation areas representing higher and lower values, respectively. On top of these results, the lahar model (with the same VEI) is displayed as input to calculate the cumulative damage over the same geo-cells exposed to both perils.
Precomputed hazard models of ash-falls and lahars are displayed by the RIESGOS demonstrator after the selection of a scenario in terms of the expected for an eruption of the Cotopaxi volcano. Local probabilistic ashfall models for the Cotopaxi volcano generated by the IG-EPN (following the method of [95] with 20-year observation of wind flow directions) are currently integrated as twelve explorative scenarios. They are represented by isolines ( Figure A1). The lahar models described in [69] are incorporated, showing the maximum possible values of five physical properties (i.e., flow velocity, flow depth, pressure, erosion, and deposition (see Figure A2). Figure 6. Example of the graphical representation of loss distribution due to ash fall scenario in the RIESGOS demonstrator (as of December 2020) from a previously selected VEI. Reddish and greenish aggregation areas representing higher and lower values, respectively. On top of these results, the lahar model (with the same VEI) is displayed as input to calculate the cumulative damage over the same geo-cells exposed to both perils.
Precomputed hazard models of ash-falls and lahars are displayed by the RIESGOS demonstrator after the selection of a scenario in terms of the expected for an eruption of the Cotopaxi volcano. Local probabilistic ashfall models for the Cotopaxi volcano generated by the IG-EPN (following the method of [95] with 20-year observation of wind flow directions) are currently integrated as twelve explorative scenarios. They are represented by isolines ( Figure A1). The lahar models described in [69] are incorporated, showing the maximum possible values of five physical properties (i.e., flow velocity, flow depth, pressure, erosion, and deposition (see Figure A2).
The exposure model provides the input to calculate the direct losses over residential building portfolios classified in specific building classes for every hazard. An example for lahar-building classes is depicted in Figure A3. These models were constrained through the use of taxonomic characteristics available in the official cadastral dataset of the GADPC (Gobierno Autónomo Descentralizado Provincial de Cotopaxi, Latacunga, Ecuador), such as roof and wall materials, and the proportions of the predominant building materials sug-gested for Latacunga in [96]. No further details are provided on the manner the building exposure models were constructed, since this is out of the scope of this paper.
The vulnerability analysis of the typical residential buildings is performed using representative building exposure models with their respective fragility functions and suitable economical consequence models. Specifically, this approach is an extension of the Performance-Based Earthquake-Engineering (PBEE) method developed by [97], which has more recently been adapted to other kinds of natural hazards. The fragility model proposed in [98] is used in lahar fragility, whilst the one in [99] is used in ash fall fragility for typical residential buildings that can be encountered in the study area. The demonstrator ultimately obtains the spatial distribution of damage and losses per individual hazard, plus the option of obtaining the cumulative damage and losses due to the action of both hazardous events using the novel method outlined in [90]. Some examples are depicted in Figures 6 and A4. No further technical details are provided, because it is out of the scope of this work. Furthermore, the demonstrator enables the visualization of the areas that might potentially get disconnected from different networks and thus the identification of cascading effects on the economic activity. The method of implementation in the systemic vulnerability analysis applied in this case is similar to the one proposed in [100]. This information can be related with census data for estimating the population that might be affected by a blackout [101]. One example of this process is depicted in Figure A5 for the interruption probabilities of the electrical power network due to the impact of a lahar.

The Recognition of the Latacunga Local Laboratory by the Local Actors of the Community
"The Latacunga Laboratory: Risk management, resilience, and adaptation to climate change" has strengthened its presence in the territory through several continuous participative activities that are aligned with the objectives mentioned in Section 3.1. For instance, the Laboratory has been recently working in materialising initiatives that were proposed by local entrepreneurs. One of them is currently working on the recovery of "Relatos de una erupcion" (Tales of an eruption), which works on rescuing the historical memory of what happened in the eruption of the Cotopaxi volcano in 1877. This has been carried out through audio-visual stories that are told by direct descendants who survived this event. This initiative enhances co-responsibility and respect for historical memory. The oral transmission of this information is an important input to generate awareness. Details about these initiatives can be found in the Latacunga Laboratory website (https://latacungaresiliente. com/rescate-de-la-memoria-historica-de-la-erupcion-del-volcan-cotopaxi/ (accessed on 1 January 2021)).

Comparative Analysis of the Social Risk Perception Factors to Natural Hazards and the Spatial Distribution of Volcanic-Related Risk Factors
The method described in Section 3.1.1 was applied to rank the volcanic risk perception for the most densely populated area in Latacunga conurbation. Making use of the 420 processed surveys as input data, the social perception to the recurrence of hazards, exposure, vulnerability, and resilience for six natural hazards likely to occur in Latacunga (i.e., earthquakes, volcanic eruptions, droughts, frosts, landslides, and floods) has been investigated. This is presented in the form of the comparative matrix shown in Table 1, which reports the mean values (for all the surveys) related to the perception of every component, as well as the computed risk index for every considered hazard. The higher the value, the greater the perception of risk. In the case of resilience, the interpretation is the opposite: the higher the value, the higher perception of resilience after a hazardous event.
The greatest concern among the inhabitants of Latacunga is their own perceived vulnerability to volcanic hazards. Remarkably, their resilience after a volcanic event scores the lowest value. This implies the community is aware that they would have great difficulty (or impossibility) to recover from the related damages. It is worth noticing that despite the fact that in the questionnaires there was no distinction made in terms of the type of volcanic hazards (lahar, ash fall/tephra fall, or ballistics) or in terms of their intensity, the collective imaginary always tended to associate the occurrence of a destructive lahar as "the volcanic hazard". Most likely, the oral transmission of the experiences of the survivors from the 1877 event has permeated the mental construction of their descendants. The mean results in terms of the percentage for the answered questionnaire that makes up the field and online surveys are depicted in Table A2. Contrary to the field survey, the online surveys score large values in the basic knowledge and reconnaissance of their exposed environment (i.e., evacuation routes, emergency committee, the existence of initiatives for risk reduction). The field surveys express that 64% of the inhabitants consider the volcanic related hazards as events that are likely to happen in the city within their lifetimes. Furthermore, 86% answered that they believe an eventual eruption of the Cotopaxi volcano would cause very serious damaging effects to the city. Likewise, 75% considered they will have very serious impacts directly on their families and themselves. 25% of the population considers that recovery from a serious volcanic event would be impossible, whilst 57% think it would be difficult to overcome. Regarding knowledge, 67% of the population know safe places in the event of a possible disaster, while 61% know evacuation routes. However, only 34% ensure there are emergency plans in their neighbourhood. Half of the respondents do not even know if they live in a volcanic hazard zone. Additionally, 56% of the field-surveyed inhabitants and 69% of the online-respondents consider they would have rapid reaction capacities. Finally,~42% of the population talks about how to act in case of emergency with their families.
Every answer of the 420 field surveys was spatially distributed onto the survey locations ( Figure 4). Their associated numerical values of the Likert scale were interpolated through the use of the ordinary kriging geostatistical algorithm [88]. Subsequently, every numerical value is converted to the equivalent categories presented in Figure 5. The spatially explicit categories represent the social perception of volcanic hazard recurrence, exposure, vulnerability, and resilience in the study area. They are respectively depicted in Figure 7a-d. The former factors are integrated through Equation (1) to generate Figure 7e, which represents the semi-quantitative volcanic risk perception index proposed in [77]. In general terms, the perceptions of hazard recurrence, exposure, and vulnerability are quite similar. However, in the central-easternmost and northernmost zones, there is a high perception of hazard recurrence, a very low perception of resilience, and a moderate perception of vulnerability, whilst in the southernmost part (where the Cutuchi River flows), the four assessed factors show high and very high values that ultimately lead to a generalized "very high" category in the volcanic risk index. This is contrary to what is observed in the central-western and northern parts. Due to the increasing distances from the main drainages, there is a strong anti-correlation between the higher resilience levels (reddish areas) and the other risk factors. Hence, despite the last volcanic crisis in 2015, there is still a generalised very low perception among the inhabitants that they cannot suffer any direct impact or damaging effect after increasing volcanic activity, because they consider the occurrence of lahars (within their lifetimes) to be impossible. Clearly, the inhabitants of that sector are not aware of the large intensities the Cotopaxi volcano can achieve (e.g., a Plinian activity (VEI > 4) in Figure 2).
four assessed factors show high and very high values that ultimately lead to a generalized "very high" category in the volcanic risk index. This is contrary to what is observed in the central-western and northern parts. Due to the increasing distances from the main drainages, there is a strong anti-correlation between the higher resilience levels (reddish areas) and the other risk factors. Hence, despite the last volcanic crisis in 2015, there is still a generalised very low perception among the inhabitants that they cannot suffer any direct impact or damaging effect after increasing volcanic activity, because they consider the occurrence of lahars (within their lifetimes) to be impossible. Clearly, the inhabitants of that sector are not aware of the large intensities the Cotopaxi volcano can achieve (e.g., a Plinian activity (VEI > 4) in Figure 2).

The Commnuication of the Scenario-Based Risk Assessment Concept with Local Stakeholders
During the four RIESGOS participative workshops, the invited stakeholders expressed interest in understanding the impacts of an extreme volcanic eruptions on the exposed elements such as buildings and critical infrastructure. Brainstorming exercises were carried out during the two first workshops. The participants were invited to imagine a future potential eruption with the emission of ash fall and occurrence of lahars. Thereafter, based on their perspectives and local knowledge, it was asked which physical, systemic, and cascading damaging effects they would expect on their built environment, infrastructure systems, and socioeconomic activities.
Some basic concepts of the probabilistic method, as an open-source web-service assesses the vulnerability of the exposed residential buildings (see Section 3.2.1), were presented to the local stakeholders. Due to the iterative approach used in constructing the demonstrator, some of the details that have been presented in this work as methods are actually the initial outputs of the first participative workshops. In this regard, the adaptability of "foreign" lahar vulnerability models (i.e., not developed for Ecuador) in the study area (e.g., [21,31,102]) was initially discussed with the representatives of the scientific local institutions. Due to the absence of locally developed ash fall vulnerability models for the residential buildings in the surroundings of the Cotopaxi volcano, the use of vulnerability models for the southern Colombian Galeras volcano [99] was perceived suitable to be implemented in the risk calculations rather than the fragility functions frequently developed for other areas (e.g., Italy [103]). With this feedback, the web-tool was redesigned. This is an example of how the engagement of local participants can improve both the technical development of quantitative methods (by agreeing on a proper model) as well as the understanding of such methods by the community.
Possible cascading effects that would occur in the case of critical infrastructure failure were debated. For instance, the participants realised that assessing the vulnerability of

The Commnuication of the Scenario-Based Risk Assessment Concept with Local Stakeholders
During the four RIESGOS participative workshops, the invited stakeholders expressed interest in understanding the impacts of an extreme volcanic eruptions on the exposed elements such as buildings and critical infrastructure. Brainstorming exercises were carried out during the two first workshops. The participants were invited to imagine a future potential eruption with the emission of ash fall and occurrence of lahars. Thereafter, based on their perspectives and local knowledge, it was asked which physical, systemic, and cascading damaging effects they would expect on their built environment, infrastructure systems, and socioeconomic activities.
Some basic concepts of the probabilistic method, as an open-source web-service assesses the vulnerability of the exposed residential buildings (see Section 3.2.1), were presented to the local stakeholders. Due to the iterative approach used in constructing the demonstrator, some of the details that have been presented in this work as methods are actually the initial outputs of the first participative workshops. In this regard, the adaptability of "foreign" lahar vulnerability models (i.e., not developed for Ecuador) in the study area (e.g., [21,31,102]) was initially discussed with the representatives of the scientific local institutions. Due to the absence of locally developed ash fall vulnerability models for the residential buildings in the surroundings of the Cotopaxi volcano, the use of vulnerability models for the southern Colombian Galeras volcano [99] was perceived suitable to be implemented in the risk calculations rather than the fragility functions frequently developed for other areas (e.g., Italy [103]). With this feedback, the web-tool was redesigned. This is an example of how the engagement of local participants can improve both the technical development of quantitative methods (by agreeing on a proper model) as well as the understanding of such methods by the community.
Possible cascading effects that would occur in the case of critical infrastructure failure were debated. For instance, the participants realised that assessing the vulnerability of electric networks to ash falls and lahars is fundamental because of the further consequences on daily social and economic activities. However, the most debated topic was the reliability of the road system that, in the case of failure, may induce physical disruption and affect evacuation and emergency response during a volcanic crisis. Other public infrastructures that would be affected by Cotopaxi's lahars include the Army headquarters "Brigada Patria", Latacunga hospital, and the new penitentiary [58]. The interest in relocating some of the exposed assets was discussed.
"Hands-on" sessions took place during the two last workshops. The participants could experience on their own the use of the RIESGOS demonstrator. They selected different scenarios to visually compare every hazard footprint and intensity (i.e., for ash fall and lahars) as well as their associated risk outcomes on residential buildings and electric power networks. This was done through the selection of individual and successive hazard scenarios addressing cumulative damage. During the "hands-on" session, the participants recognized the potential of the demonstrator as an exploration tool for risk communication.

Discussion
The CIS and RIESGOS projects have independently addressed the domain of risk communication in Latacunga (Ecuador) at different geographical scales. The investigation and mapping of the perception of volcanic risk factors led by the Latacunga Laboratory (created by CIS) was carried out in a focused area (urban area) due to the necessity of having control points (where field-surveys were carried out) for a further geostatistical interpolation process, whilst, in the framework of the RIESGOS project, the construction of the hazard, exposure, and vulnerability approaches for scenario-based multi-risk calculations have been carried out at the canton level. Despite that, the community perceptions of the entire canton and province can be assessed in the future through field surveys for other urban centres (e.g., Pujili, Saquisili, and Salcedo), a meaningful spatially explicit perception of volcanic risk factors could only be mapped for the urban centres. This is because, due to the scattered location of the residential buildings in the rural areas, conventional geostatistical interpolation algorithms would carry significant bias in the results. For the commonly investigated area by RIESGOS and CIS, we can see that the exposed community recognise that they are under a variable level of risk regarding volcanic events depending on their location. These perceptions match the lahar footprints from the scenarios with higher probabilities of occurrence (VEI < 3). However, for larger intensities, (e.g., lahar footprints from a VEI > 4 scenario, see Figure 2b), we observe a mismatch with the spatially explicit community perceptions of volcanic risk factors (Figure 7). For instance, the easternmost areas of the urban centre of Latacunga show low and very low reconnaissance of volcanic risk factors due to their increasing distance respect to the main drainages. The inhabitants of that particular sector have perceived as impossible the occurrence of and suffering from consequences of lahars. The ignorance of the lahar footprints expected from these large intensity scenarios means that the concepts of "safe place" and evacuation routes are not applicable for either. These results should not be interpreted as fixed or permanent, but they rather constitute a temporal reading of the collective mental construction of the inhabitants at the time the surveys were carried out. Nevertheless, considering that the community is placed in ancient lahar deposits, as well as the relatively short time since the last 2015 volcanic crisis, one can realise that from the comparison of the respective outcomes arises the need to prioritize some zones where further divulgation activities should be made in the future regarding the possible scenarios and intensities that the Cotopaxi volcano can actually produce.
The formulated questions comprised in the survey forms are locally revised by experts from the APGR while paying attention to the use of collectively known terminology and the cultural characteristics of the community. In this work, we have implemented a simple numerical expression (Equation (1)) that equally ranks the risks factors of the different volcanic risk factors. This selection carries epistemic uncertainties. For instance, a customisation weighting schema to each factor, the selection of the median or mode instead of the mean value (herein adopted), together with a broader range in the Likert scale (e.g., 1 to 7 as explored in [104]) could be alternative approaches to be compared or even integrating each other into condition trees as proposed in [105]. The selection of the Likert scale to rank the answers and ultimately map the community perceptions implied an ordinal scale that was further converted into a nominal one based on the "equalscale" (Equation (2)). This decision was made because, since the methods and results are aimed to be divulged, the categories have been found to be comprehensive, easily understandable, and culturally accepted by the community. Although the Likert scale has been extensively and recently used to successfully assess the community perception (e.g., [82][83][84][85][86]), there are several limitations in its adoption. For instance, as stated in [87], this kind of scale, despite maximizing the reliability of answers, also sacrifices the level of detail. However, it should be noted that through the simple possible answers related to the vulnerability perception and the nominal categories, we are only proposing a very simple categorization. More robust approaches that have addressed spatial multi-criteria analysis (as presented in [106]) have shown the impact of addressing diverse socioeconomic variables that we have not addressed in in our approach. A similar situation occurs with the resilience perception, which as discussed in [24], can be decomposed into very heterogeneous variables in economically developed countries.
Therefore, we are not claiming that our results related to the community perception of risk factors are exhaustive, but instead, they should be used as a basis for developing more complex analyses in future stages. For instance, even though we have already observed clear behaviour differences between the responses from online and field surveys, with explicitly designed survey and accounting variables such as work location, alphabetisation level, and economic activity, we could in the future classify the population into different social groups and find similarities and differences in their behaviour within a social environment to carry out more sophisticated methods, as proposed in [107]. Thereby, for each group, we could expect different reactions to a future volcanic crisis and then propose particular resilience practices. However, these kinds of approaches will largely depend on the data availability, which is particularly difficult in the rural tropics [25,57].
As described in recent participative experiences to assess the community perception to natural hazards (e.g., [45,47]) we have also experienced that the workshops carried out were allowed to go beyond a simple exchange of information. They paved the way for a better divulgation of concepts such as triggering and cascading hazards, dynamic vulnerability, cumulative damage, and cascading effects. These understandings in turn facilitated the knowledge flows and feedback acquisition to continuously design the RIESGOS demonstrator guided by increasingly risk-informed decision-makers. With this bottom-up iterative approach in the web-tool design, we are following the suggestions of the Sendai Framework for Disaster Risk Reduction (2015-2030) [73]. The outcomes of the demonstrator are not static hazard maps that are delivered to the exposed population from top-down approaches (e.g., [48][49][50]), but rather scenario-based online computations that can dynamically change based upon the continuous integration of local datasets and models.
During the "hands-on" sessions, the potential users perceived the RIESGOS demonstrator to intended prompt risk communication processes. For the study area, only hazard models have been typically available, and the few risk outcomes obtained in the past have been reported in tables and not in a spatially explicit manner [58]. Therefore, this work is providing the community with the availability of scenario-based risk models based on the vulnerability of the exposed elements in graphical and user-friendly interphase, which is an added value for the local community. The integrated scenario-based lahar footprints per VEI [69] and the locally developed probabilistic ash falls models [95] are themselves useful outcomes for civil protection and local-planners. They can be used to identify which human settlements and agricultural plantations might be affected or even discuss the relocation of some of the exposed components of critical infrastructure. Although we have not accounted for the conditional probabilities between triggering and cascading hazards as proposed in [17], we have instead presented fixed risk scenarios. For such a purpose, the demonstrator is served by a novel method that calculates and disaggregates the cumulative damage when there are interactions at the vulnerability level. In the specific volcanic context, although the concept of dynamic vulnerability had been already theoretically sketched in the work of [108], to the best of the authors' knowledge, we have first presented an example case of cumulative damage for risk-informed communities exposed to compound and cascading volcanic hazards. This is an innovative approach that not only contributes to reducing the generalized gap in the interactions at the vulnerability level [35], but also to communicating the results to the local stakeholders. With these contributions, the potential users could identify the most vulnerable areas for further mitigation strategies. It is worth mentioning that, since the RIESGOS demonstrator is currently not an operational tool, but rather shows the scientific and technological capabilities, the economic loss estimations for every exposure geo-cell (where residential buildings are aggregated) should not be used as definitive results. Therefore, due to the underlying uncertainties in these results, there is still the permanent necessity pointed out in [39,52] of having expert local users and scientists who can analyse and effectively communicate this information.
The technology transfer of the activities included in the CIS and RIESGOS programmes is highly relevant. The modular software architecture is particularly relevant for this aspect, for which the databases and methodologies of local Ecuadorian institutions may be ultimately integrated. However, the applicability of the demonstrator in the long-term will depend on how the local authorities will "give life" to the initiative, considering the local legal aspects. For future communication initiatives, due to the intrinsic interoperative sequence of inputs and outputs, the demonstrator can be a pedagogic tool to divulge multi-risk situations as similarly carried out by audio-visual approaches (e.g., [53,109]). Nevertheless, these kinds of local actors should be the first ones to understand the aforementioned concepts of "scenario" and "intensity" within the multi-risk chain, and most importantly, that they can be further contrasted with future and continuous spatially explicit social risk perceptions monitoring initiatives.

Conclusions
We have presented an integrative framework of qualitative community risk perceptions (carried out by the CIS Latacunga Laboratory) and scenario-based quantitative multi-hazard risk assessment (developed by the RIESGOS project). These initiatives have jointly worked on comprehensive volcanic risk communication processes in Latacunga, a city with a mainly rurally composed population, exposed to volcanic hazards from the Cotopaxi volcano.
Online and field surveys were carried out to rank the volcanic risk factors to investigate the individual knowledge and attitudes in Latacunga. Only the geo-located interviews in the field were used to map the community risk perceptions and to calculate a spatially explicit risk perception index through a semi-quantitative approach.
The participative workshops allowed the potentially affected communities to identify how their exposed assets, depending on their physical and systemic vulnerabilities, would be differently affected by several volcanic hazard scenarios. The iteratively customised RIESGOS demonstrator proved to be a useful tool for the communication of quantitative risk scenarios, raising the awareness of potentially affected population for the concept of scenarios and intensity. Its outcomes facilitate discussions among the participants on topics such as relocation of critical infrastructure elements. The demonstrator is not only enhancing the awareness of the communities, but also the user involvement in its development, improving the quality of the software. Although the development of the CIS and RIESGOS methodologies started independently, the respective outcomes of this collaborative work has allowed identifying areas where risk perception and scenariobased risk models are in disagreement. Thus, the need to continue assessing the social risk perception along with future risk communication efforts in the Cotopaxi region is highlighted.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the continuous development of the RIESGOS demonstrator and community risk perception carried out by GIZ. to Dorothea Kallenberger (GIZ), Ana Patricia (Grupo FARO), and Cristopher Velasco (president of the APGR) for having supported the research of the community risk perception study in Latacunga. Thanks to Luis Chasi for leading the initiative "Relatos de una erupcion". Thanks to Diego Molina (GADPC) for promoting discussions about the social risk perception of the local actors of the Cotopaxi province during the workshops held in Latacunga. Thanks to all SENESCYT and SNGER for the co-organization of the workshops held in Quito. Special thanks to Nils Brinkmann and Matthias Rüster (GFZ) for having participated in the construction of some of the computer codes implemented in this work.

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.  Figure A1. Example of the graphical representation of the spatially distributed ash fall intensities as isolines in the RIES-GOS demonstrator (as of December 2020) from a previously selected VEI. The thickness values are displayed. Once a point within the isolines is clicked, the expected load (kPa) value is also shown. Figure A2. Example of the graphical representation of the footprint and intensities of the lahars in the RIESGOS demonstrator (as of December 2020) from a previously selected VEI. On the top-right side of the window, the outputs of the lahar simulation are listed (i.e., lahar flow velocity, flow depth, pressure, erosion, and deposition). Figure A1. Example of the graphical representation of the spatially distributed ash fall intensities as isolines in the RIESGOS demonstrator (as of December 2020) from a previously selected VEI. The thickness values are displayed. Once a point within the isolines is clicked, the expected load (kPa) value is also shown.  Figure A1. Example of the graphical representation of the spatially distributed ash fall intensities as isolines in the RIES-GOS demonstrator (as of December 2020) from a previously selected VEI. The thickness values are displayed. Once a point within the isolines is clicked, the expected load (kPa) value is also shown.   There are displayed the quantities of every ash fall risk oriented building class proposed in [99] within a selected area. Figure A4. Example of the graphical representation of damage state distribution due to the combined effect of ash falls and lahar scenarios in the RIESGOS demonstrator (as of December 2020) calculated using the method proposed in [90], [92].  Figure A3. Example of the graphical representation of the residential building exposure model in the RIESGOS demonstrator (as of December 2020). It is represented into aggregation areas based on the official rural and urban administrative divisions of Latacunga. There are displayed the quantities of every ash fall risk oriented building class proposed in [99] within a selected area. Figure A4. Example of the graphical representation of damage state distribution due to the combined effect of ash falls and lahar scenarios in the RIESGOS demonstrator (as of December 2020) calculated using the method proposed in [90], [92]. Figure A4. Example of the graphical representation of damage state distribution due to the combined effect of ash falls and lahar scenarios in the RIESGOS demonstrator (as of December 2020) calculated using the method proposed in [90], [92]. Sustainability 2021, 13, x FOR PEER REVIEW 22 of 26 Figure A5. Example of the visualization of the expected interruption probabilities in the RIESGOS demonstrator (as of December 2020) of the electrical power network due to the action of a lahar scenario.