- Open Access
The influence of probabilistic volcanic hazard map properties on hazard communication
© Thompson et al.; licensee Springer. 2015
- Received: 30 June 2014
- Accepted: 6 January 2015
- Published: 11 February 2015
Probabilistic volcanic hazard analysis is becoming an increasingly popular component of volcanic risk reduction strategies worldwide. While probabilistic hazard analyses offer many advantages for decision-making, displaying the statistical results of these analyses on a map presents new hazard communication challenges. Probabilistic information is complex, difficult to interpret, and associated with uncertainties. Conveying such complicated data on a static map image without careful consideration of user perspectives or context, may result in contrasting interpretations, misunderstandings, or aversion to using the map. Here, we present the results of interviews and surveys conducted with organisational stakeholders and scientists in New Zealand which explored how probabilistic volcanic hazard map properties influence map interpretation, understanding, and preference. Our results suggest that data classification, colour scheme, content, and key expression play important roles in how users engage with and interpret probabilistic volcanic hazard maps.
Data classification was found to influence the participants’ perceived uncertainty and data reading accuracy, with isarithmic style maps reducing uncertainty and increasing accuracy best. Colour scheme had a strong influence on the type of hazard messages interpreted, with a red-yellow scheme conveying the message of a hazard distribution (high to low), and a red-yellow-blue scheme conveying the message of hazard state (present or absent) and/or risk. Multiple types of map content were found to be useful, and hazard curves were viewed as valuable supplements. The concept of “confidence” was more easily interpreted than upper and lower percentiles when expressing uncertainty on the hazard curves. Numerical and verbal expression in the key also had an influence on interpretation, with a combination of both a percent (e.g., 25%) and a natural frequency (e.g., 1 in 4) “probability” being the most inclusive and widely-understood expression. The importance of these map property choices was underscored by a high portion of participants preferring to receive maps in unalterable formats, such as PDF.
This study illustrates how engaging with users in a bottom-up approach can complement and enhance top-down approaches to volcanic hazard mapping through a collaborative and integrative design process which may help to prevent miscommunications in a future crisis when maps are likely to be drafted and disseminated rapidly.
- Hazard communication
- Hazard maps
- Probabilistic volcanic hazard analysis
- Hazard curve
Volcanic hazards are spatial phenomena, described and communicated globally using maps. In conventional hazard mapping practice, scientists analyse data about volcanic hazards, and then display the results on a map for a diverse audience of readers. However, this insular top-down approach to volcanic hazard communication does not necessarily translate into effective hazard awareness or informed decision-making, which are key contributors to risk reduction. These factors are dependent on how the reader interprets and understands the hazard map, a facet of volcanic hazard communication that has been relatively unexplored.
Wide variation in data type, hazard modelling techniques, cartographic styles, and subjective preference has led to a vast diversity in the ways volcanic hazard information is represented on maps (Calder et al. 2012), yet limited attention has been paid to the impact of these choices on reader response and comprehension. For example, does the use of a particular colour scheme (e.g., blues versus reds), data classification (e.g., stretched versus binned), or key expression (e.g., 10% versus 1 in 10) affect the way that a user interprets hazard information from the map? Geography (e.g., Robinson 1967; Bertin 1983; MacEachren 1995; Monmonier 1996; Kunz and Hurni 2011) and communication (e.g., Fisher 1991; Hoffrage et al. 2000; Severtson and Myers 2013; Doyle et al. 2014) research literature suggests that these types of properties may play an important role in hazard information processing. As capabilities and accessibility of hazard analyses continue to advance, it is important and timely to investigate how quantitative hazard map information is interpreted by those who may refer to the maps for decision-making. Initial interpretations of maps become particularly important in the context of decision-making during high-stakes, time-pressure crisis situations when reliance on affective heuristics, or intuitive feelings about the map, is likely to be high (Finucane et al. 2000).
In this New Zealand-based study, we use qualitative and quantitative methodologies to empirically investigate how certain map properties affect the way that stakeholders and scientists respond to and interpret information from probabilistic-style volcanic hazard maps, and aim to gain key map user perspectives on these hazard map properties. Understanding how map properties influence map legibility and accessibility can help build awareness in the scientific community about the impact of certain choices when designing hazard maps, and accordingly, help improve the consistency, clarity, and effectiveness of volcanic hazard communication.
Properties and users of probabilistic volcanic hazard maps
Volcanic hazard maps are thematic maps, in which the features selected as primary map content are themed on information about the potential dangers associated with volcanoes. Probabilistic volcanic hazard maps depict, in varying ways, the likelihood of a certain dangerous volcanic phenomenon(a) occurring in a particular area in a specific time frame. Computational advances and development of a number of probabilistic volcanic hazard analysis methodologies over the past decade (e.g., Mader et al. 2006) have led to a marked increase in the amount of probabilistic volcanic hazard maps being created today. There are many advantages of using a probabilistic approach to analyse volcanic hazards. Volcanoes are inherently complex systems associated with many degrees of uncertainty. Probabilistic methods allow scientists to evaluate the unknowns and to quantify and compare the likelihood of a broad range of different possible hazard scenarios, introducing more practical applications to risk assessment than traditional methods which typically analyse only one scenario (Newhall and Hoblitt 2002; Marzocchi et al. 2012). However, using probabilistic volcanic hazard data to build hazard maps also presents new challenges. Probabilities are inherently difficult to communicate and understand among experts and non-experts alike (Hoffrage et al. 2000; Spiegelhalter et al. 2011), and are interpreted and contextualised differently by different types of map users (Reyna and Brainerd 2008; Doyle et al. 2014). Conveying uncertain data using concrete map images may also introduce communication obstacles (Severtson and Myers 2013).
Many methods exist for both generalising (simplifying) and representing (depicting) thematic content on maps (Slocum 1999; Dent et al. 2009). Extensive options mean that strikingly different maps can be created from a single set of data. Such choices for representing the data may have important impacts on how the map is perceived and interpreted (Monmonier 1996). While a map may be designed to convey a single dataset of interest, it does not convey a single universal message. Each user constructs their own individual knowledge from interpreting the different symbols, colours, and expressions on the map (Bertin 1983; MacEachren 1995). Some key thematic map properties for which it is important to consider the choice of representational style are: data classification, map content, and colour scheme (Monmonier 1996; Dent et al. 2009).
Data classification refers to the characterisation and categorisation of a dataset into a particular visual symbology (Dent et al. 2009; Slocum 1999). For example, for a dataset which comprises 1,000 numerical values that range from 1 to 100, e.g., probabilistic hazard values, those values can be stretched gradationally over a gradual ramp of colours, where each value responds to a unique shade of colour on the ramp (e.g., white = 0 to black = 100 for a greyscale ramp). Isopleths or isarithms (lines connecting points of equal value) can be applied on top of the shaded classification to provide points of reference for the reader. Alternatively, the numerical values can be classed into exclusive bins which contain a specific range of values (e.g., 1 – 9, 10 – 19, and 20 – 29) and are represented by a unique colour. All values that fall together within a certain bin, or colour, are equally representative of that category. Applying different data classifications can yield markedly different patterns on the map, and can also influence how readers interpret the data (Robinson et al. 1995; Monmonier 1996; Brewer and Pickle 2002).
Map content is another important property which can influence interpretation of thematic maps (Monmonier 1996; Dent et al. 2009). Background content, such as boundaries, rivers, and landforms, are provided as points of reference for the reader, and are essential for users to be able to orient themselves with the data. Haynes et al. (2007) and Nave et al. (2010) found that selection of background content can have a significant influence on how well readers understand volcanic hazard maps, particularly among local populations living proximal to the volcano. The representation of thematic feature content also has an impact on how readers interpret the map. Probabilistic hazard datasets are typically quite robust in that they often include a range of data types (e.g., volcanic ash thickness, ash grain size) at multiple levels of uncertainty (e.g., 90th percentile, average) over a number of different time frames (e.g., in the event of an eruption, in the next 10 years). For example, some probabilistic volcanic hazard map toolsets are currently being designed to output two types of mappable content: a distribution of the probabilities of reaching and exceeding a fixed hazard intensity threshold (e.g., 10 mm of ash), or a distribution of hazard intensities expected at a fixed probability threshold (e.g., 25%) over any user-defined timeframe (Tonini et al. in press). Certain types of hazard information may have different meanings and salience to different types of map users.
Influencing more than just aesthetics, colour is a map property that has the power to provide visual contrast, unify elements, and guide attention of the map reader (Robinson 1967; Wolfe and Horowitz 2004). Blends of colour hue, saturation, and lightness can be manipulated to create a broad range of appropriate and logical colour schemes for maps. However, colour schemes also have a high margin for confusion (Brewer 1994; Robinson et al. 1995; Olson and Brewer 1997). Colours are inherently imbued with meanings for the map reader (Robinson 1967; Bertin 1983; Brewer 1994; Dent et al. 2009). While many volcanic hazard maps employ a Western green-yellow-red “stoplight” colour scheme, or a blue-yellow-red “cold to hot” colour scheme, these may have unintended meanings for some map readers (Monmonier 1996). For example, red may be associated with “hot” or “bad” and green with “vegetation” or “good”. A green-red diverging colour scheme is also likely to be very difficult to perceive for a colour blind user (Olson and Brewer 1997; Jenny and Kelso 2007; Brewer et al. 2013). Approximately 8% of males and 0.5% of females have some form of colour-vision deficiency (AOA 2014) and this can cause complications if not taken into consideration in map colour scheme design (Olson and Brewer 1997; Jenny and Kelso 2007). Although map makers strive for the highest degree of objectivity in designing probabilistic hazard maps, these fundamental choices for colour scheme, content type, and data classification are largely driven by subjective preference. Currently, there are few recommendations available for choosing appropriate styles, as the impact of these choices has not previously been tested with volcanic hazard map users.
Volcanic hazard maps have a wide range of applications and uses, and accordingly, they have many different types of users. In New Zealand, some important hazard map stakeholders are in local and regional government, health, education, agriculture, aviation, communications, emergency management, and in other groups which may use volcanic hazard information to make decisions regarding national, institutional, or local risk reduction. The detailed quantitative nature of probabilistic volcanic hazard information is well suited for organisational stakeholder decision-making, and could enhance the structure, performance, and reliability of decision-making strategies (Woo 2008, 2009; Marzocchi et al. 2012). Although each stakeholder will have different perspectives and informational needs, resource and time constraints make it impractical to develop unique hazard maps for each user group, particularly in a time of volcanic crisis (Leonard et al. 2014). It is therefore important to investigate how hazards maps can be designed to maintain a wide degree of relevance, legibility, and applicability among a diverse group of stakeholders. In order to do this, stakeholder perspectives and opinions must be taken into account, as studies show that differences exist in the way that scientists (typically the hazard map makers) and other stakeholders, such as emergency managers, understand and use probabilistic volcanic hazard information (Doyle et al. 2011, 2014). Emergency managers are a particularly important group of hazard map users in New Zealand, as they are responsible for managing many short- and long-term natural hazard and risk reduction decisions which directly affect populations at risk (MCDEM 2002).
We emphasise that people at risk are also very important stakeholders in probabilistic volcanic hazard maps and in the decisions made using such maps, but acknowledge that a survey of those at risk is outside of the scope of this study. We also acknowledge that the New Zealand stakeholder community is a generally well-educated populace, and typically has a good grasp of the concept of volcanic hazard due to regular engagement with the scientific community. Exploring probabilistic volcanic hazard map perception in other cultural, volcanic, and socio-economic settings both within New Zealand and internationally, is an important area for further research.
A pragmatic mixed-methods approach was adopted (Morgan 2007), where qualitative semi-structured in-person interviews were implemented as a pilot study to test ideas and questions which informed the development of a qualitative and quantitative online survey exploring the influence of probabilistic volcanic hazard map design properties on understanding and communicating hazard. Inductive thematic analysis of semi-structured interviews was used to identify themes regarding how scientists and stakeholders felt about certain map design properties and how they engaged with the hazard maps. Design properties which emerged as themes which have a sensitive or powerful impact on interpretation of the hazard maps were explored further in the survey among a broader sample group. Both methodological components of this study were individually approved by a human participation ethics committee.
It was emphasised to participants that the data on the maps were hypothetical, and all participants were reminded that the maps seen were not to be referred to for any type of decision-making. The datasets used were in raster format and represented the probability of accumulating a certain thickness of volcanic ash in the event of a rhyolitic eruption from Tarawera volcano in New Zealand’s North Island. Datasets were based on actual ash hazard analyses done by Thompson et al. (under review). Tarawera volcano is a well-known active volcano in New Zealand, and has had two very large explosive eruptions in the past 1 ka (Walker et al. 1984; Nairn et al. 2001). Ashfall hazard was chosen because it is one of the most common products of volcanic activity worldwide, and it is a widespread, disruptive hazard that can impact society on many different spatial, temporal, and socio-economic scales (Blong 1984). Ashfall is also a volcanic hazard that is commonly analysed using a probabilistic approach because of the need to consider variable atmospheric conditions (e.g., Scollo et al. 2008; Folch 2012). New Zealand regularly experiences small scale (~1 – 2 mm) volcanic ashfalls from its most frequently active volcanoes, with the most recent incidence of minor ashfall occurring in 2012 (Scott and Potter 2014).
Four of the maps showed the hazard data displayed using different data classification styles: gradational shaded (minimum-maximum stretch), isarithmic with solid isopleths (5% labelled intervals), isarithmic with dashed isopleths (5% labelled intervals), and binned (10% intervals). These maps were displayed using a red-yellow multi-hue sequential colour scheme similar to that used commonly in volcanic hazard maps worldwide. The fifth map displayed the data using a blue-yellow-red multi-hue diverging colour scheme similar to that used in other types of hazard maps worldwide (e.g., flooding), and the sixth map showed a deterministic scenario hazard map based on a single pre-determined eruption scenario (i.e., one wind condition and one set of eruptive parameters). A series of cards depicting numerical expressions of percent (10%), natural frequency (1 in 10), and decimal probability (0.1) and verbal expressions of “probability”, “likelihood”, and “chance” were also presented for discussion.
Participants were recruited by e-mailing informational flyers to regional and district councils near Tarawera volcano and to the national geological research institute (GNS Science), with encouragement to circulate the information among colleagues, affiliates, and peers who would be interested in participating. Initial contact was made by interested persons contacting the researchers with an expression of interest. A total of 14 people participated in the interviews, including four volcano scientists with specialties in volcano geophysics, monitoring, and hazard analysis, and ten stakeholders with specialties in emergency management, planning, infrastructure and resource management, and public health.
Full transcripts of the semi-structured interviews were analysed using inductive thematic analysis at the semantic level, based on the framework outlined by Braun and Clarke (2006). The transcripts were analysed in order to identify emergent themes regarding how map properties such as data classification, colour scheme, key expression, and content influenced the way the participant thought about or interpreted the map and its usefulness.
Summary of questions asked in the online survey
Questions in survey
Understanding of voluntary consent
Understanding of access to summary of study findings
Field of work
Primary field of work
Personal associations/relevance of volcanic hazard
Experience/memories of volcanoes
Importance of volcanic hazard maps
Data classification (Figure 2 )
Read probabilistic hazard from gradational shaded map
Read qualitative hazard from gradational shaded map
Read probabilistic hazard from binned map
Read qualitative hazard from binned map
Read probabilistic hazard from isarithmic map
Read qualitative hazard from isarithmic map
Order data classification style maps preference
Explain data classification style maps preference
Map content (Figure 3 )
Helpfulness of fixed ash/probability threshold maps
Preferred map to receive of fixed ash/probability threshold maps
Preferred key text
Influence of key text on perceived level of hazard
Hazard curves (Figure 4 )
Read hazard curve with percentile lines
Ease/difficulty of reading hazard curve with percentile lines
Read hazard curve with confidence interval area
Ease/difficulty of reading hazard curve with confidence area
Rank helpfulness of hazard curves being available with map
Map colour scheme (Figure 5 )
Preferred colour scheme set of maps
Explain reason for colour scheme preference
Influence of colour scheme on hazard perception
Explanatory text (Figure 6 )
Helpfulness of explanatory text to understanding map
Importance of explanatory text on map
Desired format for receiving volcanic hazard maps
Demographics (Table 2 )
Highest level of education
Survey participants’ occupational discipline
Participant occupational discipline
Primary field of work*
Physical/chemical/earth science (academic or scientific organisation)
Emergency management and planning
Social science (academic or scientific organisation)
Survey participant demographics
Highest level of education
Did not complete high school
Prefer not to answer
Prefer not to answer
Prefer not to answer
Descriptive and inferential statistical analyses of survey results were performed using IBM SPSS predictive analytics software. Responses to this survey generated primarily non-parametric nominal (i.e., categorical) data, and chi-square tests were used to identify the existence of significant differences within and among samples. Text responses to open-ended survey questions were qualitatively analysed using inductive thematic analysis at the semantic level to identify key patterns and themes.
Map data classification
It’s more than just communicating it [the map] to fellow scientists, or fellow emergency managers, or engineers, or whatever. It’s also how do you communicate that further down to, to members of the public. – Stakeholder 5
The advantage is you can still quite easily see the gradational nature of the probability, but at the same time you can read off the actual percentage if you feel you need to read that off. – Scientist 4
No relationship was observed between dashed (Figure 1C) and solid isopleths (Figure 1B) and interpretation of hazard information. However, the black isopleths and labels across the map were viewed as busy and distracting, detracting from the aesthetic appeal of the map. Participants overall recommended smoothing of the isopleths, citing that straight lines and right angles resulting from the square grid used in the model were visually unappealing and could falsely imply certainty. Participants also viewed the 5% intervals on the isopleths (Figure 1B,C) as misleadingly precise, and instead preferred the 10% intervals used in the binned map (Figure 1D).
Although there’s defined areas between them [bins], it’s kind of useful in terms of a planning aspect that you’ve got some boundaries to work on. – Stakeholder 6
Except, when you put a boundary on it, then people probably think if they’re on one side of the boundary or the other there’s a huge difference in probability when there isn’t. – Scientist 3
The background content was generally well liked, and viewed as helpful for orientation and location tasks, without cluttering the map. In contrast, geopolitical boundary limits on the data were not liked. Participants suggested that the boundary made the hazard appear to “stop” artificially. Participants also emphasized that managing and responding to hazards is a collaborative process that crosses boundaries, and that uncertainties or unknowns for the hazard in the surrounding regions introduced potential issues for interpretation and application.
Map colour scheme
My first impression of that [diverging map] is, blue would be safe. Whereas, where this [blue area] is yellow in the other maps, it implies there is some sort of risk that we need to consider. – Stakeholder 2
It makes it very clear, ‘that’s the area we’re worried about. We’re not worried about anything else around the district’. – Stakeholder 7
Automatically [I] relate that to weather maps… and just to complicate matters worse, we’re also doing flooding maps, and the models that we use are based very much on the same colour format. So I would look at that straight away and think, ‘Oh, it’s a flood map.’…because you see the blue. – Stakeholder 8
This was the only significant difference observed between scientists and stakeholders in the interviews. The potential difficulty in reading maps due to colour blindness was raised as an issue by several participants.
Map key expression
Depends on how accurate you want to be. Probability means yes, definitely it is 10% probable. Whereas 10% likelihood [means] well it might, it might not. - Stakeholder 9
“Chance” was perceived as untrustworthy, and akin to “slang”, invoking associations with gambling and horse-racing. Percentages (e.g., 10%) were overall the most preferred numerical expression, with participants explaining that they are a commonly used format in many different walks of life. However, some participants noted that natural frequencies (e.g., 1 in 10) were easier for people to get a “feeling” or gist for the value. Decimal values (e.g., 0.1) were universally considered as unfamiliar and too difficult for most people to understand. Many participants also commented that the gradational shaded key symbology was difficult to read with only three percent values marked (min., med., and max.).
Any ash to me is the worst case scenario. - Stakeholder 3
I think the probability is more important than the thickness because [of] the way people work. It’s, ‘Will I get ash?’ Not, ‘How much will I get?’ …They’re not sort of thinking ‘We can do this with 10 mill[imetre]s; we can’t do this with 100.’ It’s kind of ‘Oh, ok, we have to deal with volcanic ash.’ Saying that there’s a 30 or 40 percent chance…drives home a pretty strong message that they’ve got to deal with it. – Scientist 2
Many participants also suggested that addition of text explaining the possible impacts of the ash thickness would make the map more relevant and useful. It was advised that the text be placed directly onto the map face, explaining that it could easily be misplaced, truncated, or disregarded if supplied separately.
Revisions to maps for the survey
Geopolitical boundary constraint on hazard data removed (Figure 2)
Colour blind accessible colours were adopted (Figure 3)
10% bins and isopleth intervals were adopted (instead of 5%)
More intervals were labelled on the gradational stretch symbology in the key
White contour lines were used instead of black (to reduce map noise)
Isopleths were smoothed to remove the jagged artefact of the modelling grid
≥1 mm of ash was displayed instead of ≥10 mm of ash.
Out of the 110 respondents, seven participants did not fully complete the survey. Results for questions with missing values were measured using only the number of completed responses (min. 103). Results are reported for the 5% significance level (95% confidence), where χ 2 is the chi-square value for a goodness-of-fit chi-square test (one sample), Pearson χ 2 is the chi-square value for a Pearson chi-square test of independence (two or more samples), df is the degrees of freedom, and α is the level of significance.
In response to a question about whether they had ever experienced volcanic activity or had any memories associated with volcanic activity (Table 1, Question 4), 68.8% of participants responded “yes”. Of those, 64.0% described specific volcanic activity in New Zealand. When asked to rank their opinion about how important it is to have a volcanic hazard map for the potentially active volcanoes in New Zealand on a 5-point Likert scale (Table 1, Question 5), 99.1% of respondents expressed that it was important to some degree, with 73.4% ranking it “very important”.
Map data classification
More than half of respondents interpreted the qualitative level of hazard as “medium” regardless of map data classification style used (Figure 7C). The level of qualitative hazard interpreted from each map did not change for 71.8% of respondents who chose only one hazard level for all three maps. Scientists were more likely to show a change in the level of qualitative hazard interpreted from the different maps than stakeholders (Pearson χ 2 = 7.987, df = 2, α = 0.018), with 40.0% of scientists choosing two hazard levels among the three maps, and 3.3% choosing three, and 25.0%, and 0.0% of stakeholders choosing two and three levels of hazard, respectively.
Clearest because it shows that the hazard is gradational, steadily decreasing away from the source. – Scientist 30
The [actual] change is gradational, so [the binned map] is too stepwise. Gives the wrong impression. – Stakeholder 54
[It is] instantly obvious which range of values an area falls within, plus you can see where in the range it falls, so you can get quite a precise value just by quickly looking at the [isarithmic] map. – Stakeholder 87
The labels indicating the probability band removes any confusion associated with the colour symbology. – Scientist 104
Map key expression
Map colour scheme
Red always denotes hazard to me. – Stakeholder 89
The blue looks negative rather than low. The colour blue is usually associated with water. – Stakeholder 57
[The diverging map] give the impression of safety [in the blue areas] and emphasizes danger in the red areas... [The red map] gives the impression of increasing danger closer to the volcano. – Scientist 99
Red-blue seems to imply hazard and non-hazard, instead of hazard and less hazard. – Stakeholder 58
The blue de-emphasises lower likelihood areas allowing for a more risk-based focus. – Stakeholder 46
The transition between hot and cold colours should be carefully set at a standard point, as anywhere in cold colours is unlikely to receive attention when planning in government departments is carried out. – Stakeholder 91
Map explanatory text and map format
Data classification and uncertainty
Regardless of map data classification style viewed (Figure 2), more than half of survey participants considered the probability of accumulating ≥1 mm of ash in Whakatane as a “medium” level of hazard (Figure 7C). Few participants, and particularly few stakeholders, changed the qualitative descriptor used, even when their estimates of numerical value varied between maps, highlighting that qualitative descriptors of hazard may have less sensitivity for reflecting subtle perceived changes in hazard data. This demonstrates possible issues with relying on a broad verbal probability translation table for volcanic uncertainty, and emphasises the importance of including numerical values along with qualitative descriptors of probability as suggested by Doyle et al. (2011) and Budescu et al. (2009).
Estimates of hazard read from a probabilistic hazard map displaying only a gradational shaded stretch of values (Figure 2A) were less consistent than for other data classification types (Figure 7A,B), which suggests that there is low reliability in the way users may read and interpret information from this style of map. The vast majority of participants also ranked the map as least preferred (Figure 7D) or expressed dislike, which could manifest in practice as an aversion to using and applying maps with data presented in this style. However, the addition of probability isopleths dramatically improved the favourableness and readability of gradational shaded data classification. The isarithmic map (Figure 2C) was the most preferred (Figure 7D) and most accurately-read data classification style (Figure 7B). Presence of isopleths also increased the frequency of reading hazard as a range of values instead of a single value (Figure 7A), which is a more accurate estimation of hazard for the 78 km2 area of Whakatane (actual hazard ~40 – 45%, Figure 2). Expressing a range of values may have positive implications for practice, as Dieckmann et al. (2010) found that using a range of probabilities to express uncertain information is likely to be considered more credible and responsible in a decision-making context than use of a single value probability.
The isarithmic data classification (Figure 2C) also appeared to enhance the level of map engagement, by encouraging participants to critically assess the map and estimate a more precise range of probabilities. Reading the gradational shaded map (Figure 2A) required participants to spend increased effort, only to achieve a broad variation in accuracy and a low level of consistency in responses. Reading the binned map (Figure 2B) facilitated quick perception of a range of hazard with little effort, but the majority of respondents read this as “definitive” and did not elaborate beyond the 10% accuracy conveyed in the map key. When reading the isarithmic map (Figure 2C), participants were able to quickly and easily estimate a 10% range of hazard, and more participants opted to critically assess the shading within the 10% interval to increase the precision of their estimate to within 5% accuracy (Figure 7B). The non-discrete “fuzzy” appearance of gradational shading is considered an appropriate strategy to depict uncertain hazard and risk information (MacEachren 1992; Slocum 1999), and our results suggest that adding isopleths to this style effectively improves reader ability to interpolate shading for those who seek to reduce uncertainty.
Concern about precision and uncertainty in the modelled data was a primary theme regarding data classification in both the interviews and the surveys. For example, most participants found 5% interval isopleths (Figure 1B,C) unfavourable because they appeared to depict “unwarranted” precision or certainty. The prominence of uncertainty and precision in discussion of the maps suggests that it is recognised as an important component of probabilistic volcanic hazard data by both scientists and stakeholders. Donovan et al. (2012) found that effective communication and acknowledgement of uncertainty among both groups was critical for achieving useful application of probabilistic volcanic hazard and risk data in decision-making during volcanic crises at Montserrat. However, uncertainty in volcanic hazard maps is not always recognised or accepted, and it is important to explore how users perceive and approach uncertainty in maps. In the context of map properties, our results indicate that depicting 10% intervals in probability may be perceived as appropriate, adequate, and accessible for visualising and processing hazard information among a diverse group of users.
As a whole, participants expressed a general understanding that there are uncertainties associated with modelled probabilistic volcanic hazard data. However, more specific concepts of uncertainty related to the nature of the probabilistic data were less clear. The concept and meanings of lower and upper percentile values were not consistently or easily understood by the majority of participants when reading from a hazard curve (Figure 10). A hazard curve displaying the concept of confidence (Figure 4B) was better understood, with more participants also finding it easier to read than the curve with percentiles (Figure 4A). These differences have important implications for the language and graphics we use to communicate probabilistic hazard. As probabilistic tools such as volcanic hazard curves become more integrated in volcanic hazard and risk management (e.g., Tonini et al. in press), it is important to understand the contextual measures which can be implemented to improve the presentation of results and achieve clearer communication, particularly because statistical interpretation is a common challenge for many people (Reyna and Brainerd 2008).
Colour scheme and messages of hazard
The emergence of “risk and response” as a theme associated with the diverging colour scheme showed that the concept of hazard was also ambiguous in some cases. Volcanic hazard refers to the possible occurrence of a dangerous or destructive volcanic phenomenon, and can vary in frequency, magnitude, and spatial extent. In contrast, volcanic risk implies consideration of both the hazard and the vulnerabilities and capacities of people and communities in the areas exposed to the hazard (Fournier d'Albe 1979; UNISDR 2009). While the volcanology community agrees that scientists are responsible for analysing and communicating hazard, the role of scientists regarding volcanic risk and response advice is more controversial, as risk-based decision-making must take into account many social, political, and economic factors (Marzocchi et al. 2012). The diverging colour scheme (Figure 5C) was preferred by a significantly greater proportion of stakeholders than scientists, as this group may have a greater propensity for seeking risk-related information. Our results suggest that unintended messages or connotations of risk and risk-based response may be avoided by using sequential, rather than diverging, colour schemes to illustrate hazard. Diverging hues convey messages of a changing state or phenomena (Slocum 1999), such as a safe state and a dangerous state. In contrast, adjusting lightness and saturation of incremental hues on a sequential scale conveys a message of changing intensity or degree (Slocum 1999), recognised in the participants’ perception of low, but not absent, hazard in light yellow-coloured areas on the red colour scheme map (Figure 5A). Sequential hue colour schemes have also been found to convey changing intensity in health risk maps (Severtson and Vatovec 2012; Severtson and Myers 2013).
Key expression and map trust
The verbal expression “probability” was associated with greater levels of trust than “likelihood”. This suggests that verbal expression in the key can play a role in interpretation of hazard map data, as trust is a critical factor in stakeholder engagement with uncertain hazard map data (Frewer and Salter 2007; Wachinger et al. 2013). For numerical expression, preference of percentages over natural frequencies may reflect a high level of numeracy in the sampled stakeholder group which reported high education levels. While many studies suggest natural frequencies are the most widely- and well-understood numerical expression of probabilistic information (e.g., Hoffrage et al. 2000; Lipkus 2007), others (Reyna and Brainerd 2008; Keller and Siegrist 2009; Severtson and Myers 2013) advise that the effect of expression is dependent on reader numeracy skills. Using a combination of both percentages and natural frequencies on the key was the most favoured option by more than half of participants (Figure 9), and also may be the most inclusive option for reaching trust and understanding among audiences with different numeracy skills. Including natural frequency expressions may also have additional importance in the context of low probability hazard values, as users may be more likely to interpret very low probabilities as “very low hazard” when reading natural frequencies, as opposed to interpreting very low probabilities as “no hazard” when reading percentage values (Kahneman and Tversky 1979; Reyna and Brainerd 1991; Stone et al. 1994). This trend may be indicated by the ~15% of respondents who viewed the hazard as greater when presented with a natural frequency compared to a percentage.
Map content and map usefulness
Overall, participants showed an understanding for the gradational nature of volcanic ash hazard, acknowledging the “realistic” appearance of gradational shading, and dismissing a geopolitical constraint on hazard. These positive indications of understanding of volcanic ash hazard may be due to the high proportion of participants with memories or experiences of volcanoes. Equal preference of fixed ash threshold and fixed probability threshold maps indicates that there are audiences for both types of probabilistic volcanic hazard map content, and that offering different types of outputs is a valuable aspect of probabilistic volcanic hazard analysis tools. Outputs at a number of different threshold levels are likely to be useful, with results of this study suggesting even thresholds as low as 1 mm of ash can be useful information for many stakeholders. However, many participants drew upon experiences from New Zealand’s past minor ashfall events to express the perceived value in information for 1 mm of ash, and this is likely to be different in other regions, where decision-makers and communities are more frequently exposed to greater thicknesses of volcanic ash.
Map format and the importance of map-maker choices
The high preference to receive PDF, JPEG, and paper formats (Figure 13), suggest that many hazard map users do not expect to alter any of the properties on the original map (via GIS or other digital tools). Accordingly, the many representation choices for map properties made by the map makers (often volcanic hazard scientists) are important because they are likely to be carried through decision-making processes. Participants also expressed intentions to use the maps as tools for communicating with others during decision-making processes, stressing the importance of being able to clearly understand the map themselves in order to be capable of clearly and consistently explaining the hazard map to others, as recommended in a collective communication strategy by key national stakeholders (MCDEM 2010). Considering that maps are likely to be used “as is”, it is important that makers of volcanic hazard maps are aware of how choices in map properties can influence communication and understanding.
Semi-structured interviews served as a valuable pilot study methodology for acquiring rich contextual thematic data and for recognising important topics and issues, and the surveys were a valuable method for exploring and testing these topics among a broader sample group. Overall, no major discords or contrasts were recognised between the results of the two methodologies. Inclusion of open-ended text questions in the survey enhanced the richness of the survey data collected and delivered insight into motivations behind key trends. Limitations in the survey methods include possible sampling bias, question order bias, and internet access. The survey was also limited to investigation of the selective map properties and topics addressed in the questions, and other important factors or issues may exist that are not covered here. In snowball sampling, the sample is affected by the network of relationships that exist. Accordingly, more well-connected populations are more likely to be included in the sample (Patton 2002). However, we assume that the proportion of stakeholders (70%) to scientists (30%) sampled is an appropriate representation of the target population. We found this purposive sampling method was suitable for achieving a high level of participation among those who met the specific criteria related to volcanic hazard maps and decision-making, but acknowledge that potential biases could exist within the sample. For example, no statistically significant differences were recognised among age groups for all parameters. However, the 18 – 24 and 65+ year age categories may have been underrepresented for this spread of occupational disciplines (Table 3). Findings may also contrast significantly among groups with a lower rate of tertiary and postgraduate education (Table 3). This study focused on a New Zealand population sample, and important cultural differences may exist for other regions. For example, although red is globally considered a potent, emotional, active colour (Adams and Osgood 1973), it has connotations of danger in many Western cultures, but is a colour of joy in some Eastern ones, and a colour of life in New Zealand Māori culture (e.g., Jackson 1972), and these types of cultural meanings could influence the way users interpret hazard depicted in red colours (Edsall 2007).
Ash fall hazard was investigated in this study because probabilistic hazard analysis approaches are widely available and commonly employed for this type of volcanic hazard (e.g., Scollo et al. 2008; Folch 2012). Other important volcanic hazards, such as lahars (volcanic mudflows) and lava flows are becoming closer to having widely applied probabilistic analysis toolsets, but a probabilistic approach is not yet as common for mapping these gravity- and topography- driven volcanic hazards. Some of the results and findings presented here may be different for probabilistic volcanic hazard maps depicting these and other hazards.
Implications for volcanic hazard map design
Results presented here indicate that stakeholders hold strong views about certain hazard map properties, and highlight the importance of considering user perspectives in hazard map development. Opening the map design process to user contemplation, criticism, and testing in times of volcanic dormancy can help build stronger, more reliable maps and can prevent miscommunications and mistrust during a crisis, when maps are likely to be drafted, disseminated, and duplicated rapidly. The quality of visual graphics and maps presented for decision-making can play a vital and significant role in the outcomes of a crisis (Tufte 1997). A collaborative and iterative design process which considers diverse user perspectives can help identify key issues, and may help contribute to the development of a more valuable end product.
Our findings suggest that the representation of data classification, colour scheme, key expression, and content, and map presentation format all have an influence on interpretation and communication of probabilistic volcanic hazard map information. These properties have varying influence on map preference, understanding, trust, and application, which could have important effects for hazard management decision-making and incorporation into risk reduction strategies, particularly if a reader relies on intuitive gist or feeling for the information presented to some degree (Reyna and Brainerd 1991; Finucane et al. 2000). Accordingly, while research into user preference is important, as readers are more likely to use a map which appeals to them, complementary empirical measurements of map interpretation are also important because the subjectively preferred map may not always be the most effective map for decision-making (Mendonça and Delazari 2011). For example, our results show that although stakeholders preferred the diverging colour scheme (Figure 5C), a diverging colour scheme may convey unintended or false messages of risk (Figure 12C).
Our results also emphasise that there is no one-size-fits-all map presentation. One map cannot comprehensively meet the precise needs of a diverse audience of stakeholders. For example, different ash thresholds may be of concern to different stakeholder groups, and each group may use the map for different modes of communication with different audiences and communities. However, giving a voice to different stakeholder groups in an inclusive approach to hazard mapping can help map makers strive for the most integrative map possible to mediate this issue. Empirical investigation into the interplay of different map properties and how they cumulatively affect the appeal and interpretations of a map with different stakeholder groups can help map makers work towards achieving the most collaborative and effective map possible. A more consentient map can facilitate consistent messaging among stakeholder communities and stakeholder audiences, and reduce the circulation of potentially conflicting information.
Although it is important to communicate a consistent message, it may often be appropriate to create a set of maps instead of a single map, for example, a map set comprising a fixed ash threshold map (Figure 3A) and a fixed probability threshold map (Figure 3B) supplemented with available hazard curves (Figure 4). It is also important to consider that probabilistic volcanic hazard maps are built on the results of a particular input model, which uses a particular set of parameters, assumptions, and limitations. In many circumstances, it may be appropriate to present the results of more than one model. In all cases, explanatory text should be provided on the visual itself to ensure the information is contextually supported.
We further note that the challenges of displaying uncertain information on a map is not unique to volcanic hazards, but has also been discussed in fields such as hurricane hazards (Broad et al. 2007), seismic hazard (Newman et al. 2001), and health hazards (Severtson and Myers 2013). Although the target audiences and messages vary for each discipline, many of these visualisation and communication challenges are analogous or shared, and research undertaken in one field can serve as a valuable resource for guiding enquiries and applications in other hazard-related fields.
Few significant differences were recognised between stakeholder and scientist participants in this study, which may be due to limitations in sample size of the study and design of the survey questionnaire (Table 2). Recommendations for future work in this area would be to focus on identifying any key differences affecting hazard map perception among various stakeholder groups (e.g., scientists, emergency managers, and the public). People living in areas potentially exposed to volcanic hazard are very important stakeholders in hazard map information, and are a primary audience for hazard maps in the event of volcanic unrest. However, in many cases, these people are likely to be less familiar with visualisations of volcanic hazard and probabilistic data. As such, future research exploring how volcanic hazard map design influences hazard communication and understanding among these people should be undertaken in order to improve the successful transfer of important hazard information to the wider at-risk community. Further work also needs to be undertaken to identify any further properties which may influence hazard map communication and interpretation, in addition to a more probing investigation of the parameters explored here, such as testing more colour schemes and map content types. We propose that future studies in other volcanic regions or on other natural hazards could adopt a similar methodology to that presented here, and refer the reader to Additional file 1 for a full copy of the online survey.
Using an isarithmic map comprised of gradational shaded stretched data classification with smoothed and labelled 10% interval probability isopleths can be an effective format for conveying the gradual nature of ashfall hazard and for encouraging critical analysis of the data with a high level of accuracy.
Red-yellow sequential colour schemes are commonly associated with volcanoes and hazard, and reduce the potential for confusion with hydrological hazard map types. The colour scheme can effectively convey the message that areas depicted as low hazard areas still have some level of hazard present.
Diverging colour schemes may facilitate unwanted associations with volcanic risk concepts, and may enable use of the map as a reference tool for risk-based response decisions based on the transition area between hues. The colour scheme may convey messages of hazard absence in low hazard areas, and, depending on colour, may also be confused with other hazard map types (e.g., flooding hazard map). However, this scheme may be helpful for facilitating the interpretation of a high hazard “zone”.
Color Brewer (http://www.colorbrewer2.org, Brewer et al. 2013) is a helpful free online tool for selecting colour blind safe, appropriately balanced map colour schemes. Any colour scheme should be tested using a colour blind filter tool, such as those offered through Color Oracle (Jenny and Kelso 2007) or Adobe Photoshop®, before dissemination. Including labels of probability on the map face also assists vision deficient users.
When expressing uncertainty in a dataset, the concept of “confidence” may be grasped more easily and more clearly by stakeholders than lower and upper percentile data ranks, particularly in hazard curve graphics.
A percent (e.g., 25%) “probability” is a trusted numerical-verbal expression of probabilistic volcanic hazard, which is considered familiar and reliable. Including a natural frequency (e.g., 1 in 4) in addition is likely to help to increase understanding of the expression among a broader audience and may reduce the tendency to interpret very low hazard as “no hazard”.
Qualitative hazard descriptors (e.g., “low”, “medium”, and “high”) are less sensitive than numerical hazard values (e.g., 30%, 45 - 50%), and may not accurately reflect users’ perceived change in hazard probability.
Both fixed probability and fixed ash threshold maps are valuable probabilistic volcanic hazard map content types for organisational stakeholder groups.
Static PDF format is a convenient, popular, and preferred format for receiving and sharing volcanic hazard maps, emphasising the importance of map property representation choices.
Some differences exist in probabilistic volcanic hazard map interpretation between scientists and stakeholders, suggesting that engaging with users in a multi-perspective bottom-up approach to hazard mapping can complement and enhance traditional top-down approaches.
The authors would like to express thanks and appreciation to the participants of the interviews and surveys. Thanks to E. E. Hudson-Doyle and D.M. Johnston for constructive feedback and discussion, and to S. Jenkins and C. Fearnley for insightful and helpful reviews which improved this manuscript. MAT acknowledges support from the University of Auckland and the DEVORA Project (DEtermining VOlcanic Risk in Auckland) with thanks. JML gratefully acknowledges support from the New Zealand Earthquake Commission.
- Adams FM, Osgood CE (1973) A cross-cultural study of the affective meanings of color. J Cross-Cult Psychol 4(2):135–156View ArticleGoogle Scholar
- American Optometric Association (AOA) (2014) Colour vision deficiency. http://www.aoa.org/patients-and-public/eye-and-vision-problems/glossary-of-eye-and-vision-conditions/color-deficiency?sso=y. Accessed 11 June 2014
- Bertin J (1983) Semiology of Graphics (translated by William J. Berg). University of Wisconsin Press, Madison, WIGoogle Scholar
- Blong RJ (1984) Volcanic Hazards: A Sourcebook on the Effects of Eruptions. Academic Press Australia, North Ryde, AustraliaGoogle Scholar
- Braun V, Clarke V (2006) Using thematic analysis in psychology. Qual Res Psychol 3(2):77–101View ArticleGoogle Scholar
- Brewer CA (1994) Color use guidelines for mapping and visualization. In: MacEachren AM, Taylor DRF (eds) Visualization in Modern Cartography. Elsevier, Tarrytown, New York, pp 123–147View ArticleGoogle Scholar
- Brewer CA (1996) Guidelines for selecting colors for diverging schemes on maps. Cartogr J 33(2):79–86View ArticleGoogle Scholar
- Brewer CA, Harrower M, Sheesley B, Woodruff A, Heyman D (2013) ColorBrewer 2.0. The Pennyslvania State University, State College, PA, http://www.colorbrewer2.org Google Scholar
- Brewer CA, Pickle L (2002) Evaluation of methods for classifying epidemiological data on choropleth maps in series. Ann Assoc Am Geogr 92(4):662–681View ArticleGoogle Scholar
- Broad K, Leiserowitz A, Weinkle J, Steketee M (2007) Misinterpretations of the “Cone of Uncertainty” in Florida during the 2004 hurricane season. Bull Am Meteorol Soc 88(5):651–667View ArticleGoogle Scholar
- Budescu DV, Broomell S, Por H-H (2009) Improving communication of uncertainty in the reports of the intergovernmental panel on climate change. Psychol Sci 20(3):299–308View ArticleGoogle Scholar
- Calder E, Wagner K, Ogburn S (2012) One hundred and one volcanic hazard maps: what can we learn? In: Abstracts of Cities on Volcanoes 7. IAVCEI, Colima, Mexico, pp 19–23, November 2012Google Scholar
- Dent BD, Torguson J, Hodler TW (eds) (2009) Cartography: Thematic Map Design, 6th edn. McGraw-Hill Higher Education, New YorkGoogle Scholar
- Dieckmann NF, Mauro R, Slovic P (2010) The effects of presenting imprecise probabilities in intelligence forecasts. Risk Anal 30(6):987–1001View ArticleGoogle Scholar
- Donovan A, Oppenheimer C, Bravo M (2012) The use of belief-based probabilistic methods in volcanology: scientists’ views and implications for risk assessments. J Volcanol Geotherm Res 247–248:168–180View ArticleGoogle Scholar
- Doyle EEH, Johnston DM, McClure J, Paton D (2011) The communication of uncertain scientific advice during natural hazard events. N Z J Psychol 40(4):39–50Google Scholar
- Doyle EEH, McClure J, Johnston DM, Paton D (2014) Communicating likelihoods and probabilities in forecasts of volcanic eruptions. J Volcanol Geotherm Res 272:1–15View ArticleGoogle Scholar
- Edsall R (2007) Cultural factors in digital cartographic design: Implications for communication to diverse users. Cartogr Geogr Inf Sci 34:121–128View ArticleGoogle Scholar
- Finucane ML, Alhakami A, Slovic P, Johnson SM (2000) The affect heuristic in judgments of risks and benefits. J Behav Decis Mak 13(1):1–17View ArticleGoogle Scholar
- Folch A (2012) A review of tephra transport and dispersal models; evolution, current status, and future perspectives. J Volcanol Geotherm Res 235–236:96–115View ArticleGoogle Scholar
- Fournier d'Albe EM (1979) Objectives of volcanic monitoring and prediction. J Geol Soc 136(3):321–326View ArticleGoogle Scholar
- Fisher A (1991) Risk communication challenges. Risk Anal 11(2):173–179View ArticleGoogle Scholar
- Frewer L, Salter B (2007) Societal trust in risk analysis: Implications for the interface of risk assessment and risk management. In: Siegrist M, Earle T, Gutscher H (eds) Trust in Cooperative Risk Management: Uncertainty and Scepticism in the Public Mind. Earthscan, LondonGoogle Scholar
- Harrower M, Brewer CA (2003) ColorBrewer.org: An online tool for selecting colour schemes for maps. Cartogr J 40(1):27–37View ArticleGoogle Scholar
- Haynes K, Barclay J, Pidgeon N (2007) Volcanic hazard communication using maps: an evaluation of their effectiveness. Bull Volcanol 70(2):123–138View ArticleGoogle Scholar
- Hoffrage U, Lindsey S, Hertwig R, Gigerenzer G (2000) Communicating statistical information. Science 290(5500):2261–2262View ArticleGoogle Scholar
- Jackson M (1972) Aspects of symbolism and composition in Maori art. Bijdragen tot de Taal-, Land- en Volkenkunde 128(1):33–80View ArticleGoogle Scholar
- Jenny B, Kelso NV (2007) Color design for the color vision impaired. Cartographic Perspectives 57:61–67View ArticleGoogle Scholar
- Kahneman D, Tversky A (1979) Prospect theory: an analysis of decision under risk. Econometrica 47:263–291View ArticleGoogle Scholar
- Keller C, Siegrist M (2009) Effect of risk communication formats on risk perception depending on numeracy. Med Decis Making 29:483–490View ArticleGoogle Scholar
- Kunz M, Hurni L (2011) How to enhance cartographic visualisations of natural hazards assessment results. Cartogr J 48(1):60–71View ArticleGoogle Scholar
- Leonard GS, Stewart C, Wilson TM, Procter JN, Scott BJ, Keys HJ, Jolly GE, Wardman JB, Cronin SJ, McBride SK (2014) Integrating multidisciplinary science, modelling and impact data into evolving, syn-event volcanic hazard mapping and communication: A case study from the 2012 Tongariro eruption crisis, New Zealand. J Volcanol Geotherm Res 286:208–232View ArticleGoogle Scholar
- Lipkus IM (2007) Numeric, verbal, and visual formats of conveying health risks: suggested best practices and future recommendations. Med Decis Making 27(5):696–713View ArticleGoogle Scholar
- MacEachren AM (1992) Visualizing uncertain information. Cartographic Perspective 13:10–19View ArticleGoogle Scholar
- MacEachren AM (1995) How Maps Work. The Guilford Press, New YorkGoogle Scholar
- Mader HM, Coles SG, Connor CB, Connor LJ (eds) (2006) Statistics in Volcanology, vol 1. Special Publications of IAVCEI. Geological Society, LondonGoogle Scholar
- Marzocchi W, Newhall C, Woo G (2012) The scientific management of volcanic crises. J Volcanol Geotherm Res 247–248:181–189View ArticleGoogle Scholar
- MCDEM (2002) Civil Defence Emergency Management Act 2002 (CDEM Act). Ministry of Civil Defence & Emergency Management, New Zealand, Public Act 2002 No. 33Google Scholar
- MCDEM (2010) Working from the same page: consistent messages for CDEM. Ministry of Civil Defence & Emergency Management, Wellington, New ZealandGoogle Scholar
- Mendonça A, Delazari L (2011) What do People prefer and what is more effective for maps: a decision making test. In: Ruas A (ed) Advances in Cartography and GIScience. Volume 1. Lecture Notes in Geoinformation and Cartography. Springer, Berlin HeidelbergGoogle Scholar
- Monmonier M (1996) How to Lie with Maps. The University of Chicago Press, ChicagoView ArticleGoogle Scholar
- Morgan DL (2007) Paradigms lost and pragmatism regained: methodological implications of combining qualitative and quantitative methods. J Mixed Methods Res 1(1):48–76View ArticleGoogle Scholar
- Nairn IA, Self S, Cole JW, Leonard GS, Scutter C (2001) Distribution, stratigraphy, and history of proximal deposits from the c. AD 1305Kaharoa eruptive episode at Tarawera Volcano, New Zealand. N Z J Geol Geophys 44:467–484View ArticleGoogle Scholar
- Nave R, Isaia R, Vilardo G, Barclay J (2010) Re-assessing volcanic hazard maps for improving volcanic risk communication: application to Stromboli Island, Italy. J Maps 2010:260–269View ArticleGoogle Scholar
- Newhall CG, Hoblitt RP (2002) Constructing event trees for volcanic crises. Bull Volcanol 64:3–20View ArticleGoogle Scholar
- Newman A, Schneider J, Stein S, Mendez A (2001) Uncertainties in seismic hazard maps for the New Madrid Seismic Zone and implications for seismic hazard communication. Seismol Res Lett 72(6):647–663View ArticleGoogle Scholar
- Olson JM, Brewer CA (1997) An evaluation of color selections to accommodate map users with color-vision impairments. Ann Assoc Am Geogr 87(1):103–134View ArticleGoogle Scholar
- Patton MQ (ed) (2002) Qualitative Research and Evaluation Methods, 3rd edn. Sage Publications, Thousand Oaks, CAGoogle Scholar
- Reyna VF, Brainerd CJ (1991) Fuzzy-trace theory and framing effects in choice: Gist extraction, truncation, and conversion. J Behav Decis Mak 4(4):249–262View ArticleGoogle Scholar
- Reyna VF, Brainerd CJ (2008) Numeracy, ratio bias, and denominator neglect in judgments of risk and probability. Learn Individ Differ 18(1):89–107View ArticleGoogle Scholar
- Robinson AH (1967) Psychological aspects of color in cartography. International Yearbook of Cartography 7:50–59Google Scholar
- Robinson AH, Morrison JL, Muehrcke PC, Jon Kimerling A, Guptill SC (1995) Elements of Cartography, 6th edn. John Wiley & Sons, Inc, New YorkGoogle Scholar
- Scollo S, Folch A, Costa A (2008) A parametric and comparative study of different tephra fallout models. J Vol Geother Res 176:199–211View ArticleGoogle Scholar
- Scott BJ, Potter SH (2014) Aspects of historical eruptive activity and volcanic unrest at Mt. Tongariro, New Zealand: 1846–2013. J Volcanol Geotherm Res 286:263–276View ArticleGoogle Scholar
- Severtson DJ, Myers JD (2013) The influence of uncertain map features on risk beliefs and perceived ambiguity for maps of modeled cancer risk from air pollution. Risk Anal 33(5):818–837View ArticleGoogle Scholar
- Severtson DJ, Vatovec C (2012) The theory-based influence of map features on risk beliefs: Self-reports of what is seen and understood for maps depicting an environmental health hazard. J Health Commun 17:836–856View ArticleGoogle Scholar
- Slocum TA (1999) Thematic Cartography and Visualization. Prentice Hall, Upper Saddle, NJGoogle Scholar
- Spiegelhalter D, Pearson M, Short I (2011) Visualizing uncertainty about the future. Science 333(6048):1393–1400View ArticleGoogle Scholar
- Stone ER, Yates JF, Parker AM (1994) Risk communication: absolute versus relative expressions of low-probability risks. Organ Behav Hum Decis Process 60(3):387–408View ArticleGoogle Scholar
- Thompson MA, Lindsay JM, Sandri L, Biass S, Bonadonna C, Jolly G, Marzocchi W. (under review) Exploring the influence of vent location and eruption style on tephra fall hazard from the Okataina Volcanic Centre, New Zealand.Google Scholar
- Tonini R, Sandri L, Thompson MA (in press) PyBetVH: a Python tool for probabilistic volcanic hazard assessment and for generation of Bayesian hazard curves and maps, Computers and Geosciences.Google Scholar
- Tufte ER (1997) Visual and statistical thinking: displays of evidence for making decisions. In: Visual Explanations: images and Quantities, Evidence and Narrative. Graphics Press, Chesire, ConnecticutGoogle Scholar
- UNISDR (2009) UNISDR Terminology on Disaster Risk Reduction. United Nations International Strategy for Disaster Reduction (UNISDR), Geneva, SwitzerlandGoogle Scholar
- Wachinger G, Renn O, Begg C, Kuhlicke C (2013) The risk perception paradox—Implications for governance and communication of natural hazards. Risk Anal 33(6):1049–1065View ArticleGoogle Scholar
- Walker GPL, Self S, Wilson L (1984) Tarawera 1886, New Zealand - A basaltic phinian fissure eruption. J Volcanol Geotherm Res 21:61–78View ArticleGoogle Scholar
- Wolfe JM, Horowitz TS (2004) What attributes guide the deployment of visual attention and how do they do it? Nat Rev Neurosci 5:495–501View ArticleGoogle Scholar
- Woo G (2008) Probabilistic criteria for volcano evacuation decisions. Nat Hazards 45:87–97View ArticleGoogle Scholar
- Woo G (2009) A New Era of Volcano Risk Management. Risk Management Solutions, Inc, Newark, CAGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.