Using multiple data sets to populate probabilistic volcanic event trees

USGS Publications Warehouse

Newhall, C.G.; Pallister, John S.

2014-01-01

The key parameters one needs to forecast outcomes of volcanic unrest are hidden kilometers beneath the Earth’s surface, and volcanic systems are so complex that there will invariably be stochastic elements in the evolution of any unrest. Fortunately, there is sufficient regularity in behaviour that some, perhaps many, eruptions can be forecast with enough certainty for populations to be evacuated and kept safe. Volcanologists charged with forecasting eruptions must try to understand each volcanic system well enough that unrest can be interpreted in terms of pre-eruptive process, but must simultaneously recognize and convey uncertainties in their assessment. We have found that use of event trees helps to focus discussion, integrate data from multiple sources, reach consensus among scientists about both pre-eruptive process and uncertainties and, in some cases, to explain all of this to officials. Figure 1 shows a generic volcanic event tree from Newhall and Hoblitt (2002) that can be modified as needed for each specific volcano. This paper reviews how we and our colleagues have used such trees during a number of volcanic crises worldwide, for rapid hazard assessments in situations in which more formal expert elicitations could not be conducted. We describe how Multiple Data Sets can be used to estimate probabilities at each node and branch. We also present case histories of probability estimation during crises, how the estimates were used by public officials, and some suggestions for future improvements.

Scientists vs. Vesuvius: limits of volcanology

NASA Astrophysics Data System (ADS)

Carlino, Stefano; Somma, Renato

2014-05-01

Recently, Italian newspapers reported the statements of Japanese and American volcanologists which declared the high hazard related to the future occurrence of catastrophic eruption at Vesuvius. Is this a reliable picture from scientific point of view? The evaluation of volcanic hazard is based on a general statistical law for which the chances of an eruptive event increase when energy decreases. This law is constructed on the basis of empirical data. Thus, the possibility that a plinian-like eruption occurs, for each volcano, is rare and further reduced for worst-case scenario. However, empirical data are not supported by a robust scientific theory, experimentally verifiable through an exact forecast of a long-term eruption, both in time limits and in energy. Today, the lack of paradigms able to predict in a deterministic way such a complex phenomena, limit the field of the scientists that cannot go further evaluations of a purely probabilistic nature. From this point of view volcanology cannot be considered an hard quantitative Science. The declaration according to which Vesuvius, sooner or later, will produce a catastrophic eruption, yet apparently obvious if we consider the very high degree of urbanization, is not supported by any experimentally verifiable theory. Therefore, the statement according to which Vesuvius next eruptive event will be catastrophic is false. In probabilistic terms, it is actually the least possible scenario. Recognizing the cognitive limits in this research field means to encourage research itself towards the determination of more solid paradigms, in order to get more exact forecasts about such complex phenomena. The scientific compromise of defining risk scenarios, rather than deterministic evaluations about future eruptive events, precisely reflects the limits of research that have to be contemplated even by Civil Protection. Having considered these limits, every risk scenario, even the most conservative, will be ineffective in absence of an adequate political program about the reduction of the exposed value of the area and the systemic risk. In such a context, the Vesuvius area, the recent enlargement of the red zone could not represent an effective method of defence from natural disasters.

Natural disasters: forecasting economic and life losses

USGS Publications Warehouse

Nishenko, Stuart P.; Barton, Christopher C.

1997-01-01

Events such as hurricanes, earthquakes, floods, tsunamis, volcanic eruptions, and tornadoes are natural disasters because they negatively impact society, and so they must be measured and understood in human-related terms. At the U.S. Geological Survey, we have developed a new method to examine fatality and dollar-loss data, and to make probabilistic estimates of the frequency and magnitude of future events. This information is vital to large sectors of society including disaster relief agencies and insurance companies.

Learning to recognize volcanic non-eruptions

USGS Publications Warehouse

Poland, Michael P.

2010-01-01

An important goal of volcanology is to answer the questions of when, where, and how a volcano will erupt—in other words, eruption prediction. Generally, eruption predictions are based on insights from monitoring data combined with the history of the volcano. An outstanding example is the A.D. 1980–1986 lava dome growth at Mount St. Helens, Washington (United States). Recognition of a consistent pattern of precursors revealed by geophysical, geological, and geochemical monitoring enabled successful predictions of more than 12 dome-building episodes (Swanson et al., 1983). At volcanic systems that are more complex or poorly understood, probabilistic forecasts can be useful (e.g., Newhall and Hoblitt, 2002; Marzocchi and Woo, 2009). In such cases, the probabilities of different types of volcanic events are quantified, using historical accounts and geological studies of a volcano's past activity, supplemented by information from similar volcanoes elsewhere, combined with contemporary monitoring information.

A review of tephra transport and dispersal models: Evolution, current status, and future perspectives

NASA Astrophysics Data System (ADS)

Folch, A.

2012-08-01

Tephra transport models try to predict atmospheric dispersion and sedimentation of tephra depending on meteorology, particle properties, and eruption characteristics, defined by eruption column height, mass eruption rate, and vertical distribution of mass. Models are used for different purposes, from operational forecast of volcanic ash clouds to hazard assessment of tephra dispersion and fallout. The size of the erupted particles, a key parameter controlling the dynamics of particle sedimentation in the atmosphere, varies within a wide range. Largest centimetric to millimetric particles fallout at proximal to medial distances from the volcano and sediment by gravitational settling. On the other extreme, smallest micrometric to sub-micrometric particles can be transported at continental or even at global scales and are affected by other deposition and aggregation mechanisms. Different scientific communities had traditionally modeled the dispersion of these two end members. Volcanologists developed families of models suitable for lapilli and coarse ash and aimed at computing fallout deposits and for hazard assessment. In contrast, meteorologists and atmospheric scientists have traditionally used other atmospheric transport models, dealing with finer particles, for tracking motion of volcanic ash clouds and, eventually, for computing airborne ash concentrations. During the last decade, the increasing demand for model accuracy and forecast reliability has pushed on two fronts. First, the original gap between these different families of models has been filled with the emergence of multi-scale and multi-purpose models. Second, new modeling strategies including, for example, ensemble and probabilistic forecast or model data assimilation are being investigated for future implementation in models and or modeling strategies. This paper reviews the evolution of tephra transport and dispersal models during the last two decades, presents the status and limitations of the current modeling strategies, and discusses some emergent perspectives expected to be implemented at operational level during the next few years. Improvements in both real-time forecasting and long-term hazard assessment are necessary to loss prevention programs on a local, regional, national and international level.

Forecasting eruptions of Mauna Loa Volcano, Hawaii

NASA Astrophysics Data System (ADS)

Decker, Robert W.; Klein, Fred W.; Okamura, Arnold T.; Okubo, Paul G.

Past eruption patterns and various kinds of precursors are the two basic ingredients of eruption forecasts. The 39 historical eruptions of Mauna Loa from 1832 to 1984 have intervals as short as 104 days and as long as 9,165 days between the beginning of an eruption and the beginning of the next one. These recurrence times roughly fit a Poisson distribution pattern with a mean recurrence time of 1,459 days, yielding a probability of 22% (P=.22) for an eruption of Mauna Loa during any next year. The long recurrence times since 1950, however, suggest that the probability is not random, and that the current probability for an eruption during the next year may be as low as 6%. Seismicity beneath Mauna Loa increased for about two years prior to the 1975 and 1984 eruptions. Inflation of the summit area took place between eruptions with the highest rates occurring for a year or two before and after the 1975 and 1984 eruptions. Volcanic tremor beneath Mauna Loa began 51 minutes prior to the 1975 eruption and 115 minutes prior to the 1984 eruption. Eruption forecasts were published in 1975, 1976, and 1983. The 1975 and 1983 forecasts, though vaguely worded, were qualitatively correct regarding the timing of the next eruption. The 1976 forecast was more quantitative; it was wrong on timing but accurate on forecasting the location of the 1984 eruption. This paper urges that future forecasts be specific so they can be evaluated quantitatively.

Volcanic risk metrics at Mt Ruapehu, New Zealand: some background to a probabilistic eruption forecasting scheme and a cost/benefit analysis at an open conduit volcano

NASA Astrophysics Data System (ADS)

Jolly, Gill; Sandri, Laura; Lindsay, Jan; Scott, Brad; Sherburn, Steve; Jolly, Art; Fournier, Nico; Keys, Harry; Marzocchi, Warner

2010-05-01

The Bayesian Event Tree for Eruption Forecasting software (BET_EF) is a probabilistic model based on an event tree scheme that was created specifically to compute long- and short-term probabilities of different outcomes (volcanic unrest, magmatic unrest, eruption, vent location and eruption size) at long-time dormant and routinely monitored volcanoes. It is based on the assumption that upward movements of magma in a closed conduit volcano will produce detectable changes in the monitored parameters at the surface. In this perspective, the goal of BET_EF is to compute probabilities by merging information from geology, models, past data and present monitoring measurements, through a Bayesian inferential method. In the present study, we attempt to apply BET_EF to Mt Ruapehu, a very active and well-monitored volcano exhibiting the typical features of open conduit volcanoes. In such conditions, current monitoring at the surface is not necessarily able to detect short term changes at depth that may occur only seconds to minutes before an eruption. This results in so-called "blue sky eruptions" of Mt Ruapehu (for example in September 2007), that are volcanic eruptions apparently not preceded by any presently detectable signal in the current monitoring. A further complication at Mt Ruapehu arises from the well-developed hydrothermal system and the permanent crater lake sitting on top of the magmatic conduit. Both the hydrothermal system and crater lake may act to mask or change monitoring signals (if present) that magma produces deeper in the edifice. Notwithstanding these potential drawbacks, we think that an attempt to apply BET_EF at Ruapehu is worthwhile, for several reasons. First, with the exception of a few "blue sky" events, monitoring data at Mt Ruapehu can be helpful in forecasting major events, especially if a large amount of magma is intruded into the edifice and becomes available for phreatomagmatic or magmatic eruptions, as for example in 1995-96. Secondly, in setting up BET_EF for Mt Ruapehu we are forced to define quantitatively what the background activity is. This will result in a quantitative evaluation of what changes in long time monitored parameters may influence the probability of future eruptions. The slopes of Mt Ruapehu host the largest ski area in North Island, New Zealand. Lahars have been generated as a result of several eruptions in the last 50 years, and some of these have reached the ski runs in a very short time frame (around 90 seconds from the beginning of the eruption). In the light of these potentially hazardous lahars, we use the output probabilities provided by BET_EF in a practical and rational decision scheme recently proposed by Marzocchi and Woo (2009) based on a cost/benefit analysis (CBA). In such scheme, a C/L ratio is computed, based on the costs (C) of practical mitigation actions to reduce risk (e.g., a public warning scheme and other means of raising awareness, and a call for a temporary and/or partial closure of the ski area) and on the potential loss (L) if no mitigation action is taken and an eruption occurs causing lahars down the ski fields. By comparing the probability of eruption-driven lahars and the C/L ratio, it is possible to define the most rational mitigation actions that can be taken to reduce the risk to skiers, snowboarders and staff on skifield. As BET_EF probability of eruption changes dynamically as updated monitoring data are received, the authorities can decide, at any specific point in time, what is the best action according to the current monitoring of the volcano. In this respect, CBA represents a bridge linking scientific output (probabilities) and Decision Makers (practical mitigation actions).

Probabilistic short-term volcanic hazard in phases of unrest: A case study for tephra fallout

NASA Astrophysics Data System (ADS)

Selva, Jacopo; Costa, Antonio; Sandri, Laura; Macedonio, Giovanni; Marzocchi, Warner

2014-12-01

During volcanic crises, volcanologists estimate the impact of possible imminent eruptions usually through deterministic modeling of the effects of one or a few preestablished scenarios. Despite such an approach may bring an important information to the decision makers, the sole use of deterministic scenarios does not allow scientists to properly take into consideration all uncertainties, and it cannot be used to assess quantitatively the risk because the latter unavoidably requires a probabilistic approach. We present a model based on the concept of Bayesian event tree (hereinafter named BET_VH_ST, standing for Bayesian event tree for short-term volcanic hazard), for short-term near-real-time probabilistic volcanic hazard analysis formulated for any potential hazardous phenomenon accompanying an eruption. The specific goal of BET_VH_ST is to produce a quantitative assessment of the probability of exceedance of any potential level of intensity for a given volcanic hazard due to eruptions within restricted time windows (hours to days) in any area surrounding the volcano, accounting for all natural and epistemic uncertainties. BET_VH_ST properly assesses the conditional probability at each level of the event tree accounting for any relevant information derived from the monitoring system, theoretical models, and the past history of the volcano, propagating any relevant epistemic uncertainty underlying these assessments. As an application example of the model, we apply BET_VH_ST to assess short-term volcanic hazard related to tephra loading during Major Emergency Simulation Exercise, a major exercise at Mount Vesuvius that took place from 19 to 23 October 2006, consisting in a blind simulation of Vesuvius reactivation, from the early warning phase up to the final eruption, including the evacuation of a sample of about 2000 people from the area at risk. The results show that BET_VH_ST is able to produce short-term forecasts of the impact of tephra fall during a rapidly evolving crisis, accurately accounting for and propagating all uncertainties and enabling rational decision making under uncertainty.

Eruption Forecasting in Alaska: A Retrospective and Test of the Distal VT Model

NASA Astrophysics Data System (ADS)

Prejean, S. G.; Pesicek, J. D.; Wellik, J.; Cameron, C.; White, R. A.; McCausland, W. A.; Buurman, H.

2015-12-01

United States volcano observatories have successfully forecast most significant US eruptions in the past decade. However, eruptions of some volcanoes remain stubbornly difficult to forecast effectively using seismic data alone. The Alaska Volcano Observatory (AVO) has responded to 28 eruptions from 10 volcanoes since 2005. Eruptions that were not forecast include those of frequently active volcanoes with basaltic-andesite magmas, like Pavlof, Veniaminof, and Okmok volcanoes. In this study we quantify the success rate of eruption forecasting in Alaska and explore common characteristics of eruptions not forecast. In an effort to improve future forecasts, we re-examine seismic data from eruptions and known intrusive episodes in Alaska to test the effectiveness of the distal VT model commonly employed by the USGS-USAID Volcano Disaster Assistance Program (VDAP). In the distal VT model, anomalous brittle failure or volcano-tectonic (VT) earthquake swarms in the shallow crust surrounding the volcano occur as a secondary response to crustal strain induced by magma intrusion. Because the Aleutian volcanic arc is among the most seismically active regions on Earth, distinguishing distal VT earthquake swarms for eruption forecasting purposes from tectonic seismicity unrelated to volcanic processes poses a distinct challenge. In this study, we use a modified beta-statistic to identify pre-eruptive distal VT swarms and establish their statistical significance with respect to long-term background seismicity. This analysis allows us to explore the general applicability of the distal VT model and quantify the likelihood of encountering false positives in eruption forecasting using this model alone.

The meta-Gaussian Bayesian Processor of forecasts and associated preliminary experiments

NASA Astrophysics Data System (ADS)

Chen, Fajing; Jiao, Meiyan; Chen, Jing

2013-04-01

Public weather services are trending toward providing users with probabilistic weather forecasts, in place of traditional deterministic forecasts. Probabilistic forecasting techniques are continually being improved to optimize available forecasting information. The Bayesian Processor of Forecast (BPF), a new statistical method for probabilistic forecast, can transform a deterministic forecast into a probabilistic forecast according to the historical statistical relationship between observations and forecasts generated by that forecasting system. This technique accounts for the typical forecasting performance of a deterministic forecasting system in quantifying the forecast uncertainty. The meta-Gaussian likelihood model is suitable for a variety of stochastic dependence structures with monotone likelihood ratios. The meta-Gaussian BPF adopting this kind of likelihood model can therefore be applied across many fields, including meteorology and hydrology. The Bayes theorem with two continuous random variables and the normal-linear BPF are briefly introduced. The meta-Gaussian BPF for a continuous predictand using a single predictor is then presented and discussed. The performance of the meta-Gaussian BPF is tested in a preliminary experiment. Control forecasts of daily surface temperature at 0000 UTC at Changsha and Wuhan stations are used as the deterministic forecast data. These control forecasts are taken from ensemble predictions with a 96-h lead time generated by the National Meteorological Center of the China Meteorological Administration, the European Centre for Medium-Range Weather Forecasts, and the US National Centers for Environmental Prediction during January 2008. The results of the experiment show that the meta-Gaussian BPF can transform a deterministic control forecast of surface temperature from any one of the three ensemble predictions into a useful probabilistic forecast of surface temperature. These probabilistic forecasts quantify the uncertainty of the control forecast; accordingly, the performance of the probabilistic forecasts differs based on the source of the underlying deterministic control forecasts.

A volcanic event forecasting model for multiple tephra records, demonstrated on Mt. Taranaki, New Zealand

NASA Astrophysics Data System (ADS)

Damaschke, Magret; Cronin, Shane J.; Bebbington, Mark S.

2018-01-01

Robust time-varying volcanic hazard assessments are difficult to develop, because they depend upon having a complete and extensive eruptive activity record. Missing events in eruption records are endemic, due to poor preservation or erosion of tephra and other volcanic deposits. Even with many stratigraphic studies, underestimation or overestimation of eruption numbers is possible due to mis-matching tephras with similar chemical compositions or problematic age models. It is also common to have gaps in event coverage due to sedimentary records not being available in all directions from the volcano, especially downwind. Here, we examine the sensitivity of probabilistic hazard estimates using a suite of four new and two existing high-resolution tephra records located around Mt. Taranaki, New Zealand. Previous estimates were made using only single, or two correlated, tephra records. In this study, tephra data from six individual sites in lake and peat bogs covering an arc of 120° downwind of the volcano provided an excellent temporal high-resolution event record. The new data confirm a previously identified semi-regular pattern of variable eruption frequency at Mt. Taranaki. Eruption intervals exhibit a bimodal distribution, with eruptions being an average of 65 years apart, and in 2% of cases, centuries separate eruptions. The long intervals are less common than seen in earlier studies, but they have not disappeared with the inclusion of our comprehensive new dataset. Hence, the latest long interval of quiescence, since AD 1800, is unusual, but not out of character with the volcano. The new data also suggest that one of the tephra records (Lake Rotokare) used in earlier work had an old carbon effect on age determinations. This shifted ages of the affected tephras so that they were not correlated to other sites, leading to an artificially high eruption frequency in the previous combined record. New modelled time-varying frequency estimates suggest a 33-42% probability of an explosive eruption from Mt. Taranaki in the next 50 years, which is significantly lower than suggested by previous studies. This work also demonstrates some of the pitfalls to be avoided in combining stratigraphic records for eruption forecasting.

Combining probabilistic hazard assessment with cost-benefit analysis to support decision making in a volcanic crisis from the Auckland Volcanic Field, New Zealand

NASA Astrophysics Data System (ADS)

Sandri, Laura; Jolly, Gill; Lindsay, Jan; Howe, Tracy; Marzocchi, Warner

2010-05-01

One of the main challenges of modern volcanology is to provide the public with robust and useful information for decision-making in land-use planning and in emergency management. From the scientific point of view, this translates into reliable and quantitative long- and short-term volcanic hazard assessment and eruption forecasting. Because of the complexity in characterizing volcanic events, and of the natural variability of volcanic processes, a probabilistic approach is more suitable than deterministic modeling. In recent years, two probabilistic codes have been developed for quantitative short- and long-term eruption forecasting (BET_EF) and volcanic hazard assessment (BET_VH). Both of them are based on a Bayesian Event Tree, in which volcanic events are seen as a chain of logical steps of increasing detail. At each node of the tree, the probability is computed by taking into account different sources of information, such as geological and volcanological models, past occurrences, expert opinion and numerical modeling of volcanic phenomena. Since it is a Bayesian tool, the output probability is not a single number, but a probability distribution accounting for aleatory and epistemic uncertainty. In this study, we apply BET_VH in order to quantify the long-term volcanic hazard due to base surge invasion in the region around Auckland, New Zealand's most populous city. Here, small basaltic eruptions from monogenetic cones pose a considerable risk to the city in case of phreatomagmatic activity: evidence for base surges are not uncommon in deposits from past events. Currently, we are particularly focussing on the scenario simulated during Exercise Ruaumoko, a national disaster exercise based on the build-up to an eruption in the Auckland Volcanic Field. Based on recent papers by Marzocchi and Woo, we suggest a possible quantitative strategy to link probabilistic scientific output and Boolean decision making. It is based on cost-benefit analysis, in which all costs and benefits of mitigation actions have to be evaluated and compared, weighting them with the probability of occurrence of a specific threatening volcanic event. An action should be taken when the benefit of that action outweighs the costs. It is worth remarking that this strategy does not guarantee to recommend a decision that we would have taken with the benefit of hindsight. However, this strategy will be successful over the long-tem. Furthermore, it has the overwhelming advantage of providing a quantitative decision rule that is set before any emergency, and thus it will be justifiable at any stage of the process. In our present application, we are trying to set up a cost-benefit scheme for the call of an evacuation to protect people in the Auckland Volcanic Field against base surge invasion. Considering the heterogeneity of the urban environment and the size of the region at risk, we propose a cost-benefit scheme that is space dependent, to take into account higher costs when an eruption threatens sensible sites for the city and/or the nation, such as the international airport or the harbour. Finally, we compare our findings with the present Contingency Plan for Auckland.

Developing International Guidelines on Volcanic Hazard Assessments for Nuclear Facilities

NASA Astrophysics Data System (ADS)

Connor, Charles

2014-05-01

Worldwide, tremendous progress has been made in recent decades in forecasting volcanic events, such as episodes of volcanic unrest, eruptions, and the potential impacts of eruptions. Generally these forecasts are divided into two categories. Short-term forecasts are prepared in response to unrest at volcanoes, rely on geophysical monitoring and related observations, and have the goal of forecasting events on timescales of hours to weeks to provide time for evacuation of people, shutdown of facilities, and implementation of related safety measures. Long-term forecasts are prepared to better understand the potential impacts of volcanism in the future and to plan for potential volcanic activity. Long-term forecasts are particularly useful to better understand and communicate the potential consequences of volcanic events for populated areas around volcanoes and for siting critical infrastructure, such as nuclear facilities. Recent work by an international team, through the auspices of the International Atomic Energy Agency, has focused on developing guidelines for long-term volcanic hazard assessments. These guidelines have now been implemented for hazard assessment for nuclear facilities in nations including Indonesia, the Philippines, Armenia, Chile, and the United States. One any time scale, all volcanic hazard assessments rely on a geologically reasonable conceptual model of volcanism. Such conceptual models are usually built upon years or decades of geological studies of specific volcanic systems, analogous systems, and development of a process-level understanding of volcanic activity. Conceptual models are used to bound potential rates of volcanic activity, potential magnitudes of eruptions, and to understand temporal and spatial trends in volcanic activity. It is these conceptual models that provide essential justification for assumptions made in statistical model development and the application of numerical models to generate quantitative forecasts. It is a tremendous challenge in quantitative volcanic hazard assessments to encompass alternative conceptual models, and to create models that are robust to evolving understanding of specific volcanic systems by the scientific community. A central question in volcanic hazards forecasts is quantifying rates of volcanic activity. Especially for long-dormant volcanic systems, data from the geologic record may be sparse, individual events may be missing or unrecognized in the geologic record, patterns of activity may be episodic or otherwise nonstationary. This leads to uncertainty in forecasting long-term rates of activity. Hazard assessments strive to quantify such uncertainty, for example by comparing observed rates of activity with alternative parametric and nonparametric models. Numerical models are presented that characterize the spatial distribution of potential volcanic events. These spatial density models serve as the basis for application of numerical models of specific phenomena such as development of lava flow, tephra fallout, and a host of other volcanic phenomena. Monte Carlo techniques (random sampling, stratified sampling, importance sampling) are methods used to sample vent location and other key eruption parameters, such as eruption volume, magma rheology, and eruption column height for probabilistic models. The development of coupled scenarios (e.g., the probability of tephra accumulation on a slope resulting in subsequent debris flows) is also assessed through these methods, usually with the aid of event trees. The primary products of long-term forecasts are a statistical model of the conditional probability of the potential effects of volcanism, should an eruption occur, and the probability of such activity occurring. It is emphasized that hazard forecasting is an iterative process, and board consideration must be given to alternative conceptual models of volcanism, weighting of volcanological data in the analyses, and alternative statistical and numerical models. This structure is amenable to expert elicitation in order to weight alternative models and to explore alternative scenarios.

Scientific assessment of accuracy, skill and reliability of ocean probabilistic forecast products.

NASA Astrophysics Data System (ADS)

Wei, M.; Rowley, C. D.; Barron, C. N.; Hogan, P. J.

2016-02-01

As ocean operational centers are increasingly adopting and generating probabilistic forecast products for their customers with valuable forecast uncertainties, how to assess and measure these complicated probabilistic forecast products objectively is challenging. The first challenge is how to deal with the huge amount of the data from the ensemble forecasts. The second one is how to describe the scientific quality of probabilistic products. In fact, probabilistic forecast accuracy, skills, reliability, resolutions are different attributes of a forecast system. We briefly introduce some of the fundamental metrics such as the Reliability Diagram, Reliability, Resolution, Brier Score (BS), Brier Skill Score (BSS), Ranked Probability Score (RPS), Ranked Probability Skill Score (RPSS), Continuous Ranked Probability Score (CRPS), and Continuous Ranked Probability Skill Score (CRPSS). The values and significance of these metrics are demonstrated for the forecasts from the US Navy's regional ensemble system with different ensemble members. The advantages and differences of these metrics are studied and clarified.

The case for probabilistic forecasting in hydrology

NASA Astrophysics Data System (ADS)

Krzysztofowicz, Roman

2001-08-01

That forecasts should be stated in probabilistic, rather than deterministic, terms has been argued from common sense and decision-theoretic perspectives for almost a century. Yet most operational hydrological forecasting systems produce deterministic forecasts and most research in operational hydrology has been devoted to finding the 'best' estimates rather than quantifying the predictive uncertainty. This essay presents a compendium of reasons for probabilistic forecasting of hydrological variates. Probabilistic forecasts are scientifically more honest, enable risk-based warnings of floods, enable rational decision making, and offer additional economic benefits. The growing demand for information about risk and the rising capability to quantify predictive uncertainties create an unparalleled opportunity for the hydrological profession to dramatically enhance the forecasting paradigm.

Constructing event trees for volcanic crises

USGS Publications Warehouse

Newhall, C.; Hoblitt, R.

2002-01-01

Event trees are useful frameworks for discussing probabilities of possible outcomes of volcanic unrest. Each branch of the tree leads from a necessary prior event to a more specific outcome, e.g., from an eruption to a pyroclastic flow. Where volcanic processes are poorly understood, probability estimates might be purely empirical - utilizing observations of past and current activity and an assumption that the future will mimic the past or follow a present trend. If processes are better understood, probabilities might be estimated from a theoritical model, either subjectively or by numerical simulations. Use of Bayes' theorem aids in the estimation of how fresh unrest raises (or lowers) the probabilities of eruptions. Use of event trees during volcanic crises can help volcanologists to critically review their analysis of hazard, and help officials and individuals to compare volcanic risks with more familiar risks. Trees also emphasize the inherently probabilistic nature of volcano forecasts, with multiple possible outcomes.

Weighing costs and losses: A decision making game using probabilistic forecasts

NASA Astrophysics Data System (ADS)

Werner, Micha; Ramos, Maria-Helena; Wetterhall, Frederik; Cranston, Michael; van Andel, Schalk-Jan; Pappenberger, Florian; Verkade, Jan

2017-04-01

Probabilistic forecasts are increasingly recognised as an effective and reliable tool to communicate uncertainties. The economic value of probabilistic forecasts has been demonstrated by several authors, showing the benefit to using probabilistic forecasts over deterministic forecasts in several sectors, including flood and drought warning, hydropower, and agriculture. Probabilistic forecasting is also central to the emerging concept of risk-based decision making, and underlies emerging paradigms such as impact-based forecasting. Although the economic value of probabilistic forecasts is easily demonstrated in academic works, its evaluation in practice is more complex. The practical use of probabilistic forecasts requires decision makers to weigh the cost of an appropriate response to a probabilistic warning against the projected loss that would occur if the event forecast becomes reality. In this paper, we present the results of a simple game that aims to explore how decision makers are influenced by the costs required for taking a response and the potential losses they face in case the forecast flood event occurs. Participants play the role of one of three possible different shop owners. Each type of shop has losses of quite different magnitude, should a flood event occur. The shop owners are presented with several forecasts, each with a probability of a flood event occurring, which would inundate their shop and lead to those losses. In response, they have to decide if they want to do nothing, raise temporary defences, or relocate their inventory. Each action comes at a cost; and the different shop owners therefore have quite different cost/loss ratios. The game was played on four occasions. Players were attendees of the ensemble hydro-meteorological forecasting session of the 2016 EGU Assembly, professionals participating at two other conferences related to hydrometeorology, and a group of students. All audiences were familiar with the principles of forecasting and water-related risks, and one of the audiences comprised a group of experts in probabilistic forecasting. Results show that the different shop owners do take the costs of taking action and the potential losses into account in their decisions. Shop owners with a low cost/loss ratio were found to be more inclined to take actions based on the forecasts, though the absolute value of the losses also increased the willingness to take action. Little differentiation was found between the different groups of players.

Using volcanic tremor for eruption forecasting at White Island volcano (Whakaari), New Zealand

NASA Astrophysics Data System (ADS)

Chardot, Lauriane; Jolly, Arthur D.; Kennedy, Ben M.; Fournier, Nicolas; Sherburn, Steven

2015-09-01

Eruption forecasting is a challenging task because of the inherent complexity of volcanic systems. Despite remarkable efforts to develop complex models in order to explain volcanic processes prior to eruptions, the material Failure Forecast Method (FFM) is one of the very few techniques that can provide a forecast time for an eruption. However, the method requires testing and automation before being used as a real-time eruption forecasting tool at a volcano. We developed an automatic algorithm to issue forecasts from volcanic tremor increase episodes recorded by Real-time Seismic Amplitude Measurement (RSAM) at one station and optimised this algorithm for the period August 2011-January 2014 which comprises the recent unrest period at White Island volcano (Whakaari), New Zealand. A detailed residual analysis was paramount to select the most appropriate model explaining the RSAM time evolutions. In a hindsight simulation, four out of the five small eruptions reported during this period occurred within a failure window forecast by our optimised algorithm and the probability of an eruption on a day within a failure window was 0.21, which is 37 times higher than the probability of having an eruption on any day during the same period (0.0057). Moreover, the forecasts were issued prior to the eruptions by a few hours which is important from an emergency management point of view. Whereas the RSAM time evolutions preceding these four eruptions have a similar goodness-of-fit with the FFM, their spectral characteristics are different. The duration-amplitude distributions of the precursory tremor episodes support the hypothesis that several processes were likely occurring prior to these eruptions. We propose that slow rock failure and fluid flow processes are plausible candidates for the tremor source of these episodes. This hindsight exercise can be useful for future real-time implementation of the FFM at White Island. A similar methodology could also be tested at other volcanoes even if only a limited network is available.

Real-time prediction of rain-triggered lahars: incorporating seasonality and catchment recovery

NASA Astrophysics Data System (ADS)

Jones, Robbie; Manville, Vern; Peakall, Jeff; Froude, Melanie J.; Odbert, Henry M.

2017-12-01

Rain-triggered lahars are a significant secondary hydrological and geomorphic hazard at volcanoes where unconsolidated pyroclastic material produced by explosive eruptions is exposed to intense rainfall, often occurring for years to decades after the initial eruptive activity. Previous studies have shown that secondary lahar initiation is a function of rainfall parameters, source material characteristics and time since eruptive activity. In this study, probabilistic rain-triggered lahar forecasting models are developed using the lahar occurrence and rainfall record of the Belham River valley at the Soufrière Hills volcano (SHV), Montserrat, collected between April 2010 and April 2012. In addition to the use of peak rainfall intensity (PRI) as a base forecasting parameter, considerations for the effects of rainfall seasonality and catchment evolution upon the initiation of rain-triggered lahars and the predictability of lahar generation are also incorporated into these models. Lahar probability increases with peak 1 h rainfall intensity throughout the 2-year dataset and is higher under given rainfall conditions in year 1 than year 2. The probability of lahars is also enhanced during the wet season, when large-scale synoptic weather systems (including tropical cyclones) are more common and antecedent rainfall and thus levels of deposit saturation are typically increased. The incorporation of antecedent conditions and catchment evolution into logistic-regression-based rain-triggered lahar probability estimation models is shown to enhance model performance and displays the potential for successful real-time prediction of lahars, even in areas featuring strongly seasonal climates and temporal catchment recovery.

Do probabilistic forecasts lead to better decisions?

NASA Astrophysics Data System (ADS)

Ramos, M. H.; van Andel, S. J.; Pappenberger, F.

2012-12-01

The last decade has seen growing research in producing probabilistic hydro-meteorological forecasts and increasing their reliability. This followed the promise that, supplied with information about uncertainty, people would take better risk-based decisions. In recent years, therefore, research and operational developments have also start putting attention to ways of communicating the probabilistic forecasts to decision makers. Communicating probabilistic forecasts includes preparing tools and products for visualization, but also requires understanding how decision makers perceive and use uncertainty information in real-time. At the EGU General Assembly 2012, we conducted a laboratory-style experiment in which several cases of flood forecasts and a choice of actions to take were presented as part of a game to participants, who acted as decision makers. Answers were collected and analyzed. In this paper, we present the results of this exercise and discuss if indeed we make better decisions on the basis of probabilistic forecasts.

Do probabilistic forecasts lead to better decisions?

NASA Astrophysics Data System (ADS)

Ramos, M. H.; van Andel, S. J.; Pappenberger, F.

2013-06-01

The last decade has seen growing research in producing probabilistic hydro-meteorological forecasts and increasing their reliability. This followed the promise that, supplied with information about uncertainty, people would take better risk-based decisions. In recent years, therefore, research and operational developments have also started focusing attention on ways of communicating the probabilistic forecasts to decision-makers. Communicating probabilistic forecasts includes preparing tools and products for visualisation, but also requires understanding how decision-makers perceive and use uncertainty information in real time. At the EGU General Assembly 2012, we conducted a laboratory-style experiment in which several cases of flood forecasts and a choice of actions to take were presented as part of a game to participants, who acted as decision-makers. Answers were collected and analysed. In this paper, we present the results of this exercise and discuss if we indeed make better decisions on the basis of probabilistic forecasts.

Probabilistically modeling lava flows with MOLASSES

NASA Astrophysics Data System (ADS)

Richardson, J. A.; Connor, L.; Connor, C.; Gallant, E.

2017-12-01

Modeling lava flows through Cellular Automata methods enables a computationally inexpensive means to quickly forecast lava flow paths and ultimate areal extents. We have developed a lava flow simulator, MOLASSES, that forecasts lava flow inundation over an elevation model from a point source eruption. This modular code can be implemented in a deterministic fashion with given user inputs that will produce a single lava flow simulation. MOLASSES can also be implemented in a probabilistic fashion where given user inputs define parameter distributions that are randomly sampled to create many lava flow simulations. This probabilistic approach enables uncertainty in input data to be expressed in the model results and MOLASSES outputs a probability map of inundation instead of a determined lava flow extent. Since the code is comparatively fast, we use it probabilistically to investigate where potential vents are located that may impact specific sites and areas, as well as the unconditional probability of lava flow inundation of sites or areas from any vent. We have validated the MOLASSES code to community-defined benchmark tests and to the real world lava flows at Tolbachik (2012-2013) and Pico do Fogo (2014-2015). To determine the efficacy of the MOLASSES simulator at accurately and precisely mimicking the inundation area of real flows, we report goodness of fit using both model sensitivity and the Positive Predictive Value, the latter of which is a Bayesian posterior statistic. Model sensitivity is often used in evaluating lava flow simulators, as it describes how much of the lava flow was successfully modeled by the simulation. We argue that the positive predictive value is equally important in determining how good a simulator is, as it describes the percentage of the simulation space that was actually inundated by lava.

Collaborative Cyber-infrastructures for the Management of the UNESCO-IGCP Research Project "Forecast of tephra fallout"

NASA Astrophysics Data System (ADS)

Folch, A.; Costa, A.; Cordoba, G.

2009-04-01

Tephra fallout following explosive volcanic eruptions produces several hazardous effects on inhabitants, infrastructure, and property and represents a serious threat for communities located around active volcanoes. In order to mitigate the effects on the surrounding areas, scientists and civil decision-making authorities need reliable short-term forecasts during episodes of eruptive crisis and long-term probabilistic maps to plan territorial policies and land use. Modelling, together with field studies and volcano monitoring, constitutes an indispensable tool to achieve these objectives. The UNESCO-IGCP research project proposal "Forecast of tephra fallout" has the aim to produce a series of tools capable to elaborate both short-term forecasts and long-term hazard assessments using the cutting-edge models for tephra transport and sedimentation. A special project website will be designed to supply a set of models, procedures and expertise to several Latino-American Institutes based in countries seriously threatened by this geo-hazard (Argentina, Chile, Colombia, Ecuador, Mexico, and Nicaragua). This will proportionate to the final users a tool to elaborate short-term forecasts of tephra deposition on the ground, and determine airborne ash concentrations (a quantity of special relevance for aerial navigation safety) during eruptions and emergencies. The project web-site will have a public section and a password-protected area to exchange information and data among participants and, eventually, to allow remote execution of high-resolution mesoscale meteorological forecasts at the BSC facilities. The public website section will be updated periodically and will include sections describing the project objectives and achievements as well as the hazard maps for the investigated volcanoes, and will be linked to other relevant websites such as IAVCEI, IGCP, IUGS and UNESCO homepages. A part of the public section of the website will be devoted to disseminate achieved scientific results, provide general advice, and display hazard maps to a larger public beyond the scientific community. The website private section will include a software and documentation download section as well as a gateway to run the WRF mesoscale meteorological model and the parallel version of the FALL3D model at the BSC facilities. It will be invaluable during an eventual emergency if the affected institution does not yet have an agreement with its national weather service.

Comparison of the economic impact of different wind power forecast systems for producers

NASA Astrophysics Data System (ADS)

Alessandrini, S.; Davò, F.; Sperati, S.; Benini, M.; Delle Monache, L.

2014-05-01

Deterministic forecasts of wind production for the next 72 h at a single wind farm or at the regional level are among the main end-users requirement. However, for an optimal management of wind power production and distribution it is important to provide, together with a deterministic prediction, a probabilistic one. A deterministic forecast consists of a single value for each time in the future for the variable to be predicted, while probabilistic forecasting informs on probabilities for potential future events. This means providing information about uncertainty (i.e. a forecast of the PDF of power) in addition to the commonly provided single-valued power prediction. A significant probabilistic application is related to the trading of energy in day-ahead electricity markets. It has been shown that, when trading future wind energy production, using probabilistic wind power predictions can lead to higher benefits than those obtained by using deterministic forecasts alone. In fact, by using probabilistic forecasting it is possible to solve economic model equations trying to optimize the revenue for the producer depending, for example, on the specific penalties for forecast errors valid in that market. In this work we have applied a probabilistic wind power forecast systems based on the "analog ensemble" method for bidding wind energy during the day-ahead market in the case of a wind farm located in Italy. The actual hourly income for the plant is computed considering the actual selling energy prices and penalties proportional to the unbalancing, defined as the difference between the day-ahead offered energy and the actual production. The economic benefit of using a probabilistic approach for the day-ahead energy bidding are evaluated, resulting in an increase of 23% of the annual income for a wind farm owner in the case of knowing "a priori" the future energy prices. The uncertainty on price forecasting partly reduces the economic benefit gained by using a probabilistic energy forecast system.

Probabilistic approach to decision making under uncertainty during volcanic crises. Retrospective analysis of the 2011 eruption of El Hierro, in the Canary Islands

NASA Astrophysics Data System (ADS)

Sobradelo, Rosa; Martí, Joan; Kilburn, Christopher; López, Carmen

2014-05-01

Understanding the potential evolution of a volcanic crisis is crucial to improving the design of effective mitigation strategies. This is especially the case for volcanoes close to densely-populated regions, where inappropriate decisions may trigger widespread loss of life, economic disruption and public distress. An outstanding goal for improving the management of volcanic crises, therefore, is to develop objective, real-time methodologies for evaluating how an emergency will develop and how scientists communicate with decision makers. Here we present a new model BADEMO (Bayesian Decision Model) that applies a general and flexible, probabilistic approach to managing volcanic crises. The model combines the hazard and risk factors that decision makers need for a holistic analysis of a volcanic crisis. These factors include eruption scenarios and their probabilities of occurrence, the vulnerability of populations and their activities, and the costs of false alarms and failed forecasts. The model can be implemented before an emergency, to identify actions for reducing the vulnerability of a district; during an emergency, to identify the optimum mitigating actions and how these may change as new information is obtained; and after an emergency, to assess the effectiveness of a mitigating response and, from the results, to improve strategies before another crisis occurs. As illustrated by a retrospective analysis of the 2011 eruption of El Hierro, in the Canary Islands, BADEMO provides the basis for quantifying the uncertainty associated with each recommended action as an emergency evolves, and serves as a mechanism for improving communications between scientists and decision makers.

Global forecasting of thermal health hazards: the skill of probabilistic predictions of the Universal Thermal Climate Index (UTCI).

PubMed

Pappenberger, F; Jendritzky, G; Staiger, H; Dutra, E; Di Giuseppe, F; Richardson, D S; Cloke, H L

2015-03-01

Although over a hundred thermal indices can be used for assessing thermal health hazards, many ignore the human heat budget, physiology and clothing. The Universal Thermal Climate Index (UTCI) addresses these shortcomings by using an advanced thermo-physiological model. This paper assesses the potential of using the UTCI for forecasting thermal health hazards. Traditionally, such hazard forecasting has had two further limitations: it has been narrowly focused on a particular region or nation and has relied on the use of single 'deterministic' forecasts. Here, the UTCI is computed on a global scale, which is essential for international health-hazard warnings and disaster preparedness, and it is provided as a probabilistic forecast. It is shown that probabilistic UTCI forecasts are superior in skill to deterministic forecasts and that despite global variations, the UTCI forecast is skilful for lead times up to 10 days. The paper also demonstrates the utility of probabilistic UTCI forecasts on the example of the 2010 heat wave in Russia.

Short-term volcanic hazard assessment through Bayesian inference: retrospective application to the Pinatubo 1991 volcanic crisis

NASA Astrophysics Data System (ADS)

Sobradelo, Rosa; Martí, Joan

2015-01-01

One of the most challenging aspects of managing a volcanic crisis is the interpretation of the monitoring data, so as to anticipate to the evolution of the unrest and implement timely mitigation actions. An unrest episode may include different stages or time intervals of increasing activity that may or may not precede a volcanic eruption, depending on the causes of the unrest (magmatic, geothermal or tectonic). Therefore, one of the main goals in monitoring volcanic unrest is to forecast whether or not such increase of activity will end up with an eruption, and if this is the case, how, when, and where this eruption will take place. As an alternative method to expert elicitation for assessing and merging monitoring data and relevant past information, we present a probabilistic method to transform precursory activity into the probability of experiencing a significant variation by the next time interval (i.e. the next step in the unrest), given its preceding evolution, and by further estimating the probability of the occurrence of a particular eruptive scenario combining monitoring and past data. With the 1991 Pinatubo volcanic crisis as a reference, we have developed such a method to assess short-term volcanic hazard using Bayesian inference.

Clustering and Hazard Estimation in the Auckland Volcanic Field, New Zealand

NASA Astrophysics Data System (ADS)

Cronin, S. J.; Bebbington, M. S.

2009-12-01

The Auckland Volcanic Field (AVF) with its 49 eruptive centres formed over the last c. 250 ka presents several unique challenges to our understanding of distributed volcanic field construction and evolution. Due to the youth of the field, high-resolution stratigraphy of eruption centres and ash-fall sequences is possible, allowing time-breaks, soil and peat formation between eruption units to be identified. Radiocarbon dating of sediments between volcanic deposits shows that at least five of the centres have erupted on more than one occasion, with time breaks of 50-100 years between episodes. In addition, paleomagnetic and ash fall evidence implies that there has been strong clustering of eruption events over time, with a specific “flare-up” event involving over possibly up to 19 eruptions occurring between 35-25 ka, in spatially disparate locations. An additional complicating factor is that the only centre that shows any major evidence for evolution out of standard alkali basaltic compositions is also the youngest and largest in volume by several orders of magnitude. All of these features of the AVF, along with relatively poor age-control for many of the vents make spatio-temporal hazard forecasting for the field based on assumptions of past behaviour extremely difficult. Any relationships that take volumetric considerations into account are particularly difficult, since any trend analysis produces unreasonably large future eruptions. The most reasonable model is spatial, via eruption location. We have re-examined the age progression of eruptive events in the AVF, incorporating the most reliable sources of age and stratigraphic data, including developing new correlations between ashfall records in cores and likely vent locations via a probabilistic model of tephra dispersal. A Monte Carlo procedure using the age-progression, stratigraphy and dating constraints can then randomly reproduce likely orderings of events in the field. These were fitted by a clustering-based model of vent locations as originally applied by Magill et al (2005: Mathematical Geol. 37: 227-242) to the Allen and Smith (1994; Geosci. Report Shizuoka Univ 20: 5-14) age ordering of volcanism at AVF. Applying this model, modified by allowing continuation of activity at or around the youngest event, to sampled age orderings from the Monte Carlo procedure shows a very different spatial forecast to the earlier analysis. It is also different to the distribution from randomly ordered events, implying there is at least some clustering control on the location of eruptions in the field. Further iterations of this modelling approach will be tested in relation to eruptive volume and applied to other comparative volcanic fields.

Forecasting volcanic ash dispersal and coeval resuspension during the April-May 2015 Calbuco eruption

NASA Astrophysics Data System (ADS)

Reckziegel, F.; Bustos, E.; Mingari, L.; Báez, W.; Villarosa, G.; Folch, A.; Collini, E.; Viramonte, J.; Romero, J.; Osores, S.

2016-07-01

Atmospheric dispersion of volcanic ash from explosive eruptions or from subsequent fallout deposit resuspension causes a range of impacts and disruptions on human activities and ecosystems. The April-May 2015 Calbuco eruption in Chile involved eruption and resuspension activities. We overview the chronology, effects, and products resulting from these events, in order to validate an operational forecast strategy for tephra dispersal. The modelling strategy builds on coupling the meteorological Weather Research and Forecasting (WRF/ARW) model with the FALL3D dispersal model for eruptive and resuspension processes. The eruption modelling considers two distinct particle granulometries, a preliminary first guess distribution used operationally when no field data was available yet, and a refined distribution based on field measurements. Volcanological inputs were inferred from eruption reports and results from an Argentina-Chilean ash sample data network, which performed in-situ sampling during the eruption. In order to validate the modelling strategy, results were compared with satellite retrievals and ground deposit measurements. Results indicate that the WRF-FALL3D modelling system can provide reasonable forecasts in both eruption and resuspension modes, particularly when the adjusted granulometry is considered. The study also highlights the importance of having dedicated datasets of active volcanoes furnishing first-guess model inputs during the early stages of an eruption.

Benchmarking computational fluid dynamics models of lava flow simulation for hazard assessment, forecasting, and risk management

USGS Publications Warehouse

Dietterich, Hannah; Lev, Einat; Chen, Jiangzhi; Richardson, Jacob A.; Cashman, Katharine V.

2017-01-01

Numerical simulations of lava flow emplacement are valuable for assessing lava flow hazards, forecasting active flows, designing flow mitigation measures, interpreting past eruptions, and understanding the controls on lava flow behavior. Existing lava flow models vary in simplifying assumptions, physics, dimensionality, and the degree to which they have been validated against analytical solutions, experiments, and natural observations. In order to assess existing models and guide the development of new codes, we conduct a benchmarking study of computational fluid dynamics (CFD) models for lava flow emplacement, including VolcFlow, OpenFOAM, FLOW-3D, COMSOL, and MOLASSES. We model viscous, cooling, and solidifying flows over horizontal planes, sloping surfaces, and into topographic obstacles. We compare model results to physical observations made during well-controlled analogue and molten basalt experiments, and to analytical theory when available. Overall, the models accurately simulate viscous flow with some variability in flow thickness where flows intersect obstacles. OpenFOAM, COMSOL, and FLOW-3D can each reproduce experimental measurements of cooling viscous flows, and OpenFOAM and FLOW-3D simulations with temperature-dependent rheology match results from molten basalt experiments. We assess the goodness-of-fit of the simulation results and the computational cost. Our results guide the selection of numerical simulation codes for different applications, including inferring emplacement conditions of past lava flows, modeling the temporal evolution of ongoing flows during eruption, and probabilistic assessment of lava flow hazard prior to eruption. Finally, we outline potential experiments and desired key observational data from future flows that would extend existing benchmarking data sets.

Probabilistic Wind Power Ramp Forecasting Based on a Scenario Generation Method

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wang, Qin; Florita, Anthony R; Krishnan, Venkat K

Wind power ramps (WPRs) are particularly important in the management and dispatch of wind power and currently drawing the attention of balancing authorities. With the aim to reduce the impact of WPRs for power system operations, this paper develops a probabilistic ramp forecasting method based on a large number of simulated scenarios. An ensemble machine learning technique is first adopted to forecast the basic wind power forecasting scenario and calculate the historical forecasting errors. A continuous Gaussian mixture model (GMM) is used to fit the probability distribution function (PDF) of forecasting errors. The cumulative distribution function (CDF) is analytically deduced.more » The inverse transform method based on Monte Carlo sampling and the CDF is used to generate a massive number of forecasting error scenarios. An optimized swinging door algorithm is adopted to extract all the WPRs from the complete set of wind power forecasting scenarios. The probabilistic forecasting results of ramp duration and start-time are generated based on all scenarios. Numerical simulations on publicly available wind power data show that within a predefined tolerance level, the developed probabilistic wind power ramp forecasting method is able to predict WPRs with a high level of sharpness and accuracy.« less

Probabilistic Wind Power Ramp Forecasting Based on a Scenario Generation Method: Preprint

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wang, Qin; Florita, Anthony R; Krishnan, Venkat K

2017-08-31

Wind power ramps (WPRs) are particularly important in the management and dispatch of wind power, and they are currently drawing the attention of balancing authorities. With the aim to reduce the impact of WPRs for power system operations, this paper develops a probabilistic ramp forecasting method based on a large number of simulated scenarios. An ensemble machine learning technique is first adopted to forecast the basic wind power forecasting scenario and calculate the historical forecasting errors. A continuous Gaussian mixture model (GMM) is used to fit the probability distribution function (PDF) of forecasting errors. The cumulative distribution function (CDF) ismore » analytically deduced. The inverse transform method based on Monte Carlo sampling and the CDF is used to generate a massive number of forecasting error scenarios. An optimized swinging door algorithm is adopted to extract all the WPRs from the complete set of wind power forecasting scenarios. The probabilistic forecasting results of ramp duration and start time are generated based on all scenarios. Numerical simulations on publicly available wind power data show that within a predefined tolerance level, the developed probabilistic wind power ramp forecasting method is able to predict WPRs with a high level of sharpness and accuracy.« less

Volcanic Eruption Forecasts From Accelerating Rates of Drumbeat Long-Period Earthquakes

NASA Astrophysics Data System (ADS)

Bell, Andrew F.; Naylor, Mark; Hernandez, Stephen; Main, Ian G.; Gaunt, H. Elizabeth; Mothes, Patricia; Ruiz, Mario

2018-02-01

Accelerating rates of quasiperiodic "drumbeat" long-period earthquakes (LPs) are commonly reported before eruptions at andesite and dacite volcanoes, and promise insights into the nature of fundamental preeruptive processes and improved eruption forecasts. Here we apply a new Bayesian Markov chain Monte Carlo gamma point process methodology to investigate an exceptionally well-developed sequence of drumbeat LPs preceding a recent large vulcanian explosion at Tungurahua volcano, Ecuador. For more than 24 hr, LP rates increased according to the inverse power law trend predicted by material failure theory, and with a retrospectively forecast failure time that agrees with the eruption onset within error. LPs resulted from repeated activation of a single characteristic source driven by accelerating loading, rather than a distributed failure process, showing that similar precursory trends can emerge from quite different underlying physics. Nevertheless, such sequences have clear potential for improving forecasts of eruptions at Tungurahua and analogous volcanoes.

Visualising probabilistic flood forecast information: expert preferences and perceptions of best practice in uncertainty communication

NASA Astrophysics Data System (ADS)

Pappenberger, F.; Stephens, E. M.; Thielen, J.; Salomon, P.; Demeritt, D.; van Andel, S.; Wetterhall, F.; Alfieri, L.

2011-12-01

The aim of this paper is to understand and to contribute to improved communication of the probabilistic flood forecasts generated by Hydrological Ensemble Prediction Systems (HEPS) with particular focus on the inter expert communication. Different users are likely to require different kinds of information from HEPS and thus different visualizations. The perceptions of this expert group are important both because they are the designers and primary users of existing HEPS. Nevertheless, they have sometimes resisted the release of uncertainty information to the general public because of doubts about whether it can be successfully communicated in ways that would be readily understood to non-experts. In this paper we explore the strengths and weaknesses of existing HEPS visualization methods and thereby formulate some wider recommendations about best practice for HEPS visualization and communication. We suggest that specific training on probabilistic forecasting would foster use of probabilistic forecasts with a wider range of applications. The result of a case study exercise showed that there is no overarching agreement between experts on how to display probabilistic forecasts and what they consider essential information that should accompany plots and diagrams. In this paper we propose a list of minimum properties that, if consistently displayed with probabilistic forecasts, would make the products more easily understandable.

Improvements on the relationship between plume height and mass eruption rate: Implications for volcanic ash cloud forecasting

NASA Astrophysics Data System (ADS)

Webley, P. W.; Dehn, J.; Mastin, L. G.; Steensen, T. S.

2011-12-01

Volcanic ash plumes and the dispersing clouds into the atmosphere are a hazard for local populations as well as for the aviation industry. Volcanic ash transport and dispersion (VATD) models, used to forecast the movement of these hazardous ash emissions, require eruption source parameters (ESP) such as plume height, eruption rate and duration. To estimate mass eruption rate, empirical relationships with observed plume height have been applied. Theoretical relationships defined by Morton et al. (1956) and Wilson et al. (1976) use default values for the environmental lapse rate (ELR), thermal efficiency, density of ash, specific heat capacity, initial temperature of the erupted material and final temperature of the material. Each volcano, based on its magma type, has a different density, specific heat capacity and initial eruptive temperature compared to these default parameters, and local atmospheric conditions can produce a very different ELR. Our research shows that a relationship between plume height and mass eruption rate can be defined for each eruptive event for each volcano. Additionally, using the one-dimensional modeling program, Plumeria, our analysis assesses the importance of factors such as vent diameter and eruption velocity on the relationship between the eruption rate and measured plume height. Coupling such a tool with a VATD model should improve pre-eruptive forecasts of ash emissions downwind and lead to improvements in ESP data that VATD models use for operational volcanic ash cloud forecasting.

A novel visualisation tool for climate services: a case study of temperature extremes and human mortality in Europe

NASA Astrophysics Data System (ADS)

Lowe, R.; Ballester, J.; Robine, J.; Herrmann, F. R.; Jupp, T. E.; Stephenson, D.; Rodó, X.

2013-12-01

Users of climate information often require probabilistic information on which to base their decisions. However, communicating information contained within a probabilistic forecast presents a challenge. In this paper we demonstrate a novel visualisation technique to display ternary probabilistic forecasts on a map in order to inform decision making. In this method, ternary probabilistic forecasts, which assign probabilities to a set of three outcomes (e.g. low, medium, and high risk), are considered as a point in a triangle of barycentric coordinates. This allows a unique colour to be assigned to each forecast from a continuum of colours defined on the triangle. Colour saturation increases with information gain relative to the reference forecast (i.e. the long term average). This provides additional information to decision makers compared with conventional methods used in seasonal climate forecasting, where one colour is used to represent one forecast category on a forecast map (e.g. red = ';dry'). We use the tool to present climate-related mortality projections across Europe. Temperature and humidity are related to human mortality via location-specific transfer functions, calculated using historical data. Daily mortality data at the NUTS2 level for 16 countries in Europe were obtain from 1998-2005. Transfer functions were calculated for 54 aggregations in Europe, defined using criteria related to population and climatological similarities. Aggregations are restricted to fall within political boundaries to avoid problems related to varying adaptation policies between countries. A statistical model is fit to cold and warm tails to estimate future mortality using forecast temperatures, in a Bayesian probabilistic framework. Using predefined categories of temperature-related mortality risk, we present maps of probabilistic projections for human mortality at seasonal to decadal time scales. We demonstrate the information gained from using this technique compared to more traditional methods to display ternary probabilistic forecasts. This technique allows decision makers to identify areas where the model predicts with certainty area-specific heat waves or cold snaps, in order to effectively target resources to those areas most at risk, for a given season or year. It is hoped that this visualisation tool will facilitate the interpretation of the probabilistic forecasts not only for public health decision makers but also within a multi-sectoral climate service framework.

Performance of the 'material Failure Forecast Method' in real-time situations: A Bayesian approach applied on effusive and explosive eruptions

NASA Astrophysics Data System (ADS)

Boué, A.; Lesage, P.; Cortés, G.; Valette, B.; Reyes-Dávila, G.; Arámbula-Mendoza, R.; Budi-Santoso, A.

2016-11-01

Most attempts of deterministic eruption forecasting are based on the material Failure Forecast Method (FFM). This method assumes that a precursory observable, such as the rate of seismic activity, can be described by a simple power law which presents a singularity at a time close to the eruption onset. Until now, this method has been applied only in a small number of cases, generally for forecasts in hindsight. In this paper, a rigorous Bayesian approach of the FFM designed for real-time applications is applied. Using an automatic recognition system, seismo-volcanic events are detected and classified according to their physical mechanism and time series of probability distributions of the rates of events are calculated. At each time of observation, a Bayesian inversion provides estimations of the exponent of the power law and of the time of eruption, together with their probability density functions. Two criteria are defined in order to evaluate the quality and reliability of the forecasts. Our automated procedure has allowed the analysis of long, continuous seismic time series: 13 years from Volcán de Colima, Mexico, 10 years from Piton de la Fournaise, Reunion Island, France, and several months from Merapi volcano, Java, Indonesia. The new forecasting approach has been applied to 64 pre-eruptive sequences which present various types of dominant seismic activity (volcano-tectonic or long-period events) and patterns of seismicity with different level of complexity. This has allowed us to test the FFM assumptions, to determine in which conditions the method can be applied, and to quantify the success rate of the forecasts. 62% of the precursory sequences analysed are suitable for the application of FFM and half of the total number of eruptions are successfully forecast in hindsight. In real-time, the method allows for the successful forecast of 36% of all the eruptions considered. Nevertheless, real-time forecasts are successful for 83% of the cases that fulfil the reliability criteria. Therefore, good confidence on the method is obtained when the reliability criteria are met.

Willingness-to-pay for a probabilistic flood forecast: a risk-based decision-making game

NASA Astrophysics Data System (ADS)

Arnal, Louise; Ramos, Maria-Helena; Coughlan de Perez, Erin; Cloke, Hannah Louise; Stephens, Elisabeth; Wetterhall, Fredrik; van Andel, Schalk Jan; Pappenberger, Florian

2016-08-01

Probabilistic hydro-meteorological forecasts have over the last decades been used more frequently to communicate forecast uncertainty. This uncertainty is twofold, as it constitutes both an added value and a challenge for the forecaster and the user of the forecasts. Many authors have demonstrated the added (economic) value of probabilistic over deterministic forecasts across the water sector (e.g. flood protection, hydroelectric power management and navigation). However, the richness of the information is also a source of challenges for operational uses, due partially to the difficulty in transforming the probability of occurrence of an event into a binary decision. This paper presents the results of a risk-based decision-making game on the topic of flood protection mitigation, called "How much are you prepared to pay for a forecast?". The game was played at several workshops in 2015, which were attended by operational forecasters and academics working in the field of hydro-meteorology. The aim of this game was to better understand the role of probabilistic forecasts in decision-making processes and their perceived value by decision-makers. Based on the participants' willingness-to-pay for a forecast, the results of the game show that the value (or the usefulness) of a forecast depends on several factors, including the way users perceive the quality of their forecasts and link it to the perception of their own performances as decision-makers.

Visualizing Uncertainty for Probabilistic Weather Forecasting based on Reforecast Analogs

NASA Astrophysics Data System (ADS)

Pelorosso, Leandro; Diehl, Alexandra; Matković, Krešimir; Delrieux, Claudio; Ruiz, Juan; Gröeller, M. Eduard; Bruckner, Stefan

2016-04-01

Numerical weather forecasts are prone to uncertainty coming from inaccuracies in the initial and boundary conditions and lack of precision in numerical models. Ensemble of forecasts partially addresses these problems by considering several runs of the numerical model. Each forecast is generated with different initial and boundary conditions and different model configurations [GR05]. The ensembles can be expressed as probabilistic forecasts, which have proven to be very effective in the decision-making processes [DE06]. The ensemble of forecasts represents only some of the possible future atmospheric states, usually underestimating the degree of uncertainty in the predictions [KAL03, PH06]. Hamill and Whitaker [HW06] introduced the "Reforecast Analog Regression" (RAR) technique to overcome the limitations of ensemble forecasting. This technique produces probabilistic predictions based on the analysis of historical forecasts and observations. Visual analytics provides tools for processing, visualizing, and exploring data to get new insights and discover hidden information patterns in an interactive exchange between the user and the application [KMS08]. In this work, we introduce Albero, a visual analytics solution for probabilistic weather forecasting based on the RAR technique. Albero targets at least two different type of users: "forecasters", who are meteorologists working in operational weather forecasting and "researchers", who work in the construction of numerical prediction models. Albero is an efficient tool for analyzing precipitation forecasts, allowing forecasters to make and communicate quick decisions. Our solution facilitates the analysis of a set of probabilistic forecasts, associated statistical data, observations and uncertainty. A dashboard with small-multiples of probabilistic forecasts allows the forecasters to analyze at a glance the distribution of probabilities as a function of time, space, and magnitude. It provides the user with a more accurate measure of forecast uncertainty that could result in better decision-making. It offers different level of abstractions to help with the recalibration of the RAR method. It also has an inspection tool that displays the selected analogs, their observations and statistical data. It gives the users access to inner parts of the method, unveiling hidden information. References [GR05] GNEITING T., RAFTERY A. E.: Weather forecasting with ensemble methods. Science 310, 5746, 248-249, 2005. [KAL03] KALNAY E.: Atmospheric modeling, data assimilation and predictability. Cambridge University Press, 2003. [PH06] PALMER T., HAGEDORN R.: Predictability of weather and climate. Cambridge University Press, 2006. [HW06] HAMILL T. M., WHITAKER J. S.: Probabilistic quantitative precipitation forecasts based on reforecast analogs: Theory and application. Monthly Weather Review 134, 11, 3209-3229, 2006. [DE06] DEITRICK S., EDSALL R.: The influence of uncertainty visualization on decision making: An empirical evaluation. Springer, 2006. [KMS08] KEIM D. A., MANSMANN F., SCHNEIDEWIND J., THOMAS J., ZIEGLER H.: Visual analytics: Scope and challenges. Springer, 2008.

Global link between deformation and volcanic eruption quantified by satellite imagery.

PubMed

Biggs, J; Ebmeier, S K; Aspinall, W P; Lu, Z; Pritchard, M E; Sparks, R S J; Mather, T A

2014-04-03

A key challenge for volcanological science and hazard management is that few of the world's volcanoes are effectively monitored. Satellite imagery covers volcanoes globally throughout their eruptive cycles, independent of ground-based monitoring, providing a multidecadal archive suitable for probabilistic analysis linking deformation with eruption. Here we show that, of the 198 volcanoes systematically observed for the past 18 years, 54 deformed, of which 25 also erupted. For assessing eruption potential, this high proportion of deforming volcanoes that also erupted (46%), together with the proportion of non-deforming volcanoes that did not erupt (94%), jointly represent indicators with 'strong' evidential worth. Using a larger catalogue of 540 volcanoes observed for 3 years, we demonstrate how this eruption-deformation relationship is influenced by tectonic, petrological and volcanic factors. Satellite technology is rapidly evolving and routine monitoring of the deformation status of all volcanoes from space is anticipated, meaning probabilistic approaches will increasingly inform hazard decisions and strategic development.

How much are you prepared to PAY for a forecast?

NASA Astrophysics Data System (ADS)

Arnal, Louise; Coughlan, Erin; Ramos, Maria-Helena; Pappenberger, Florian; Wetterhall, Fredrik; Bachofen, Carina; van Andel, Schalk Jan

2015-04-01

Probabilistic hydro-meteorological forecasts are a crucial element of the decision-making chain in the field of flood prevention. The operational use of probabilistic forecasts is increasingly promoted through the development of new novel state-of-the-art forecast methods and numerical skill is continuously increasing. However, the value of such forecasts for flood early-warning systems is a topic of diverging opinions. Indeed, the word value, when applied to flood forecasting, is multifaceted. It refers, not only to the raw cost of acquiring and maintaining a probabilistic forecasting system (in terms of human and financial resources, data volume and computational time), but also and most importantly perhaps, to the use of such products. This game aims at investigating this point. It is a willingness to pay game, embedded in a risk-based decision-making experiment. Based on a ``Red Cross/Red Crescent, Climate Centre'' game, it is a contribution to the international Hydrologic Ensemble Prediction Experiment (HEPEX). A limited number of probabilistic forecasts will be auctioned to the participants; the price of these forecasts being market driven. All participants (irrespective of having bought or not a forecast set) will then be taken through a decision-making process to issue warnings for extreme rainfall. This game will promote discussions around the topic of the value of forecasts for decision-making in the field of flood prevention.

Flash-flood early warning using weather radar data: from nowcasting to forecasting

NASA Astrophysics Data System (ADS)

Liechti, Katharina; Panziera, Luca; Germann, Urs; Zappa, Massimiliano

2013-04-01

In our study we explore the limits of radar-based forecasting for hydrological runoff prediction. Two novel probabilistic radar-based forecasting chains for flash-flood early warning are investigated in three catchments in the Southern Swiss Alps and set in relation to deterministic discharge forecast for the same catchments. The first probabilistic radar-based forecasting chain is driven by NORA (Nowcasting of Orographic Rainfall by means of Analogues), an analogue-based heuristic nowcasting system to predict orographic rainfall for the following eight hours. The second probabilistic forecasting system evaluated is REAL-C2, where the numerical weather prediction COSMO-2 is initialized with 25 different initial conditions derived from a four-day nowcast with the radar ensemble REAL. Additionally, three deterministic forecasting chains were analysed. The performance of these five flash-flood forecasting systems was analysed for 1389 hours between June 2007 and December 2010 for which NORA forecasts were issued, due to the presence of orographic forcing. We found a clear preference for the probabilistic approach. Discharge forecasts perform better when forced by NORA rather than by a persistent radar QPE for lead times up to eight hours and for all discharge thresholds analysed. The best results were, however, obtained with the REAL-C2 forecasting chain, which was also remarkably skilful even with the highest thresholds. However, for regions where REAL cannot be produced, NORA might be an option for forecasting events triggered by orographic forcing.

Flash-flood early warning using weather radar data: from nowcasting to forecasting

NASA Astrophysics Data System (ADS)

Liechti, K.; Panziera, L.; Germann, U.; Zappa, M.

2013-01-01

This study explores the limits of radar-based forecasting for hydrological runoff prediction. Two novel probabilistic radar-based forecasting chains for flash-flood early warning are investigated in three catchments in the Southern Swiss Alps and set in relation to deterministic discharge forecast for the same catchments. The first probabilistic radar-based forecasting chain is driven by NORA (Nowcasting of Orographic Rainfall by means of Analogues), an analogue-based heuristic nowcasting system to predict orographic rainfall for the following eight hours. The second probabilistic forecasting system evaluated is REAL-C2, where the numerical weather prediction COSMO-2 is initialized with 25 different initial conditions derived from a four-day nowcast with the radar ensemble REAL. Additionally, three deterministic forecasting chains were analysed. The performance of these five flash-flood forecasting systems was analysed for 1389 h between June 2007 and December 2010 for which NORA forecasts were issued, due to the presence of orographic forcing. We found a clear preference for the probabilistic approach. Discharge forecasts perform better when forced by NORA rather than by a persistent radar QPE for lead times up to eight hours and for all discharge thresholds analysed. The best results were, however, obtained with the REAL-C2 forecasting chain, which was also remarkably skilful even with the highest thresholds. However, for regions where REAL cannot be produced, NORA might be an option for forecasting events triggered by orographic precipitation.

Understanding causality and uncertainty in volcanic observations: An example of forecasting eruptive activity on Soufrière Hills Volcano, Montserrat

NASA Astrophysics Data System (ADS)

Sheldrake, T. E.; Aspinall, W. P.; Odbert, H. M.; Wadge, G.; Sparks, R. S. J.

2017-07-01

Following a cessation in eruptive activity it is important to understand how a volcano will behave in the future and when it may next erupt. Such an assessment can be based on the volcano's long-term pattern of behaviour and insights into its current state via monitoring observations. We present a Bayesian network that integrates these two strands of evidence to forecast future eruptive scenarios using expert elicitation. The Bayesian approach provides a framework to quantify the magmatic causes in terms of volcanic effects (i.e., eruption and unrest). In October 2013, an expert elicitation was performed to populate a Bayesian network designed to help forecast future eruptive (in-)activity at Soufrière Hills Volcano. The Bayesian network was devised to assess the state of the shallow magmatic system, as a means to forecast the future eruptive activity in the context of the long-term behaviour at similar dome-building volcanoes. The findings highlight coherence amongst experts when interpreting the current behaviour of the volcano, but reveal considerable ambiguity when relating this to longer patterns of volcanism at dome-building volcanoes, as a class. By asking questions in terms of magmatic causes, the Bayesian approach highlights the importance of using short-term unrest indicators from monitoring data as evidence in long-term forecasts at volcanoes. Furthermore, it highlights potential biases in the judgements of volcanologists and identifies sources of uncertainty in terms of magmatic causes rather than scenario-based outcomes.

Improvement forecasting of volcanic activity by applying a Kalman filter to the SSEM signal. The case of the El Hierro Island eruption (October 2011)

NASA Astrophysics Data System (ADS)

Garcia, A.; Berrocoso, M.; Marrero, J. M.; Ortiz, R.

2012-04-01

The FFM (Failure Forecast Method) is developed from the eruption of St. Helens, being repeatedly applied to forecast eruptions and recently to the prediction of seismic activity in active volcanic areas. The underwater eruption of El Hierro Island has been monitored from three months before starting (October 10, 2011). This allowed a large catalogue of seismic events (over 11000) and continuous recording seismic signals that cover the entire period. Since the beginning of the seismic-volcanic crisis (July 2011), the FFM was applied to the SSEM signal of seismic records. Mainly because El Hierro is a very small island, the SSEM has a high noise (traffic and oceanic noise). To improve the signal / noise ratio has been used a Kalman filter. The Kalman filter coefficients are adjusted using an inversion process based on forecasting errors occurred in the twenty days preceding. The application of this filter has been a significant improvement in the reliability of forecasts. The analysis of the results shows, before the start of the eruption, that 90% of the forecasts are obtained with errors less than 10 minutes with more than 24 hours in advance. It is noteworthy that the method predicts the events of greater magnitude and especially the beginning of each swarm of seismic events. At the time the eruption starts reducing the efficiency of the forecast 50% with a dispersion of more than one hour. This fact is probably due to decreased detectability by saturation of some of the seismic stations and decreased the average magnitude. However, the events of magnitude greater than 4 were predicted with an error less than 20 minutes.

Global link between deformation and volcanic eruption quantified by satellite imagery

PubMed Central

Biggs, J.; Ebmeier, S. K.; Aspinall, W. P.; Lu, Z.; Pritchard, M. E.; Sparks, R. S. J.; Mather, T. A.

2014-01-01

A key challenge for volcanological science and hazard management is that few of the world’s volcanoes are effectively monitored. Satellite imagery covers volcanoes globally throughout their eruptive cycles, independent of ground-based monitoring, providing a multidecadal archive suitable for probabilistic analysis linking deformation with eruption. Here we show that, of the 198 volcanoes systematically observed for the past 18 years, 54 deformed, of which 25 also erupted. For assessing eruption potential, this high proportion of deforming volcanoes that also erupted (46%), together with the proportion of non-deforming volcanoes that did not erupt (94%), jointly represent indicators with ‘strong’ evidential worth. Using a larger catalogue of 540 volcanoes observed for 3 years, we demonstrate how this eruption–deformation relationship is influenced by tectonic, petrological and volcanic factors. Satellite technology is rapidly evolving and routine monitoring of the deformation status of all volcanoes from space is anticipated, meaning probabilistic approaches will increasingly inform hazard decisions and strategic development. PMID:24699342

Sensibility analysis of VORIS lava-flow simulations: application to Nyamulagira volcano, Democratic Republic of Congo

NASA Astrophysics Data System (ADS)

Syavulisembo, A. M.; Havenith, H.-B.; Smets, B.; d'Oreye, N.; Marti, J.

2015-03-01

Assessment and management of volcanic risk are important scientific, economic, and political issues, especially in densely populated areas threatened by volcanoes. The Virunga area in the Democratic Republic of Congo, with over 1 million inhabitants, has to cope permanently with the threat posed by the active Nyamulagira and Nyiragongo volcanoes. During the past century, Nyamulagira erupted at intervals of 1-4 years - mostly in the form of lava flows - at least 30 times. Its summit and flank eruptions lasted for periods of a few days up to more than two years, and produced lava flows sometimes reaching distances of over 20 km from the volcano, thereby affecting very large areas and having a serious impact on the region of Virunga. In order to identify a useful tool for lava flow hazard assessment at the Goma Volcano Observatory (GVO), we tested VORIS 2.0.1 (Felpeto et al., 2007), a freely available software (http://www.gvb-csic.es) based on a probabilistic model that considers topography as the main parameter controlling lava flow propagation. We tested different Digital Elevation Models (DEM) - SRTM1, SRTM3, and ASTER GDEM - to analyze the sensibility of the input parameters of VORIS 2.0.1 in simulation of recent historical lava-flow for which the pre-eruption topography is known. The results obtained show that VORIS 2.0.1 is a quick, easy-to-use tool for simulating lava-flow eruptions and replicates to a high degree of accuracy the eruptions tested. In practice, these results will be used by GVO to calibrate VORIS model for lava flow path forecasting during new eruptions, hence contributing to a better volcanic crisis management.

Modelling the dynamics and hazards of explosive eruptions: Where we are now, and confronting the next challenges (Sergey Soloviev Medal Lecture)

NASA Astrophysics Data System (ADS)

Neri, Augusto

2017-04-01

Understanding of explosive eruption dynamics and assessment of their hazards continue to represent challenging issues to the present-day volcanology community. This is largely due to the complex and diverse nature of the phenomena, and the variability and unpredictability of volcanic processes. Nevertheless, important and continuing progress has been made in the last few decades in understanding fundamental processes and in forecasting the occurrences of these phenomena, thanks to significant advances in field, experimental and theoretical modeling investigations. For over four decades, for example, volcanologists have made major progress in the description of the nature of explosive eruptions, considerably aided by the development, improvement, and application of physical-mathematical models. Integral steady-state homogeneous flow models were first used to investigate the different controlling mechanisms and to infer the genesis and evolution of the phenomena. Through continuous improvements and quantum-leap developments, a variety of transient, 3D, multiphase flow models of volcanic phenomena now can implement state-of-the-art formulations of the underlying physics, new-generation analytical and experimental data, as well as high-performance computational techniques. These numerical models have proved to be able to provide key insights in the understanding of the dynamics of explosive eruptions (e.g. convective plumes, collapsing columns, pyroclastic density currents, short-lived explosions, etc.), as well as to represent a valuable tool in the quantification of potential eruptive scenarios and associated hazards. Simplified models based on a reduction of the system complexity have been also proved useful, combined with Monte Carlo and statistical methods, to generate quantitative probabilistic hazard maps at different space and time scales, some including the quantification of important sources of uncertainty. Nevertheless, the development of physical models able to accurately replicate, within acceptable statistical uncertainty, the evolution of explosive eruptions remains a challenging goal still to be achieved. Testing of the developed models versus large-scale experimental data and well-measured real events, real-time assimilation of observational data to forecast the process nature and evolution, as well as the quantification of the uncertainties affecting our system and modelling representations appear key next steps to further progress volcanological research and its essential contribution to the mitigation of volcanic risk.

Beyond eruptive scenarios: assessing tephra fallout hazard from Neapolitan volcanoes.

PubMed

Sandri, Laura; Costa, Antonio; Selva, Jacopo; Tonini, Roberto; Macedonio, Giovanni; Folch, Arnau; Sulpizio, Roberto

2016-04-12

Assessment of volcanic hazards is necessary for risk mitigation. Typically, hazard assessment is based on one or a few, subjectively chosen representative eruptive scenarios, which use a specific combination of eruptive sizes and intensities to represent a particular size class of eruption. While such eruptive scenarios use a range of representative members to capture a range of eruptive sizes and intensities in order to reflect a wider size class, a scenario approach neglects to account for the intrinsic variability of volcanic eruptions, and implicitly assumes that inter-class size variability (i.e. size difference between different eruptive size classes) dominates over intra-class size variability (i.e. size difference within an eruptive size class), the latter of which is treated as negligible. So far, no quantitative study has been undertaken to verify such an assumption. Here, we adopt a novel Probabilistic Volcanic Hazard Analysis (PVHA) strategy, which accounts for intrinsic eruptive variabilities, to quantify the tephra fallout hazard in the Campania area. We compare the results of the new probabilistic approach with the classical scenario approach. The results allow for determining whether a simplified scenario approach can be considered valid, and for quantifying the bias which arises when full variability is not accounted for.

Online probabilistic learning with an ensemble of forecasts

NASA Astrophysics Data System (ADS)

Thorey, Jean; Mallet, Vivien; Chaussin, Christophe

2016-04-01

Our objective is to produce a calibrated weighted ensemble to forecast a univariate time series. In addition to a meteorological ensemble of forecasts, we rely on observations or analyses of the target variable. The celebrated Continuous Ranked Probability Score (CRPS) is used to evaluate the probabilistic forecasts. However applying the CRPS on weighted empirical distribution functions (deriving from the weighted ensemble) may introduce a bias because of which minimizing the CRPS does not produce the optimal weights. Thus we propose an unbiased version of the CRPS which relies on clusters of members and is strictly proper. We adapt online learning methods for the minimization of the CRPS. These methods generate the weights associated to the members in the forecasted empirical distribution function. The weights are updated before each forecast step using only past observations and forecasts. Our learning algorithms provide the theoretical guarantee that, in the long run, the CRPS of the weighted forecasts is at least as good as the CRPS of any weighted ensemble with weights constant in time. In particular, the performance of our forecast is better than that of any subset ensemble with uniform weights. A noteworthy advantage of our algorithm is that it does not require any assumption on the distributions of the observations and forecasts, both for the application and for the theoretical guarantee to hold. As application example on meteorological forecasts for photovoltaic production integration, we show that our algorithm generates a calibrated probabilistic forecast, with significant performance improvements on probabilistic diagnostic tools (the CRPS, the reliability diagram and the rank histogram).

Maximizing Statistical Power When Verifying Probabilistic Forecasts of Hydrometeorological Events

NASA Astrophysics Data System (ADS)

DeChant, C. M.; Moradkhani, H.

2014-12-01

Hydrometeorological events (i.e. floods, droughts, precipitation) are increasingly being forecasted probabilistically, owing to the uncertainties in the underlying causes of the phenomenon. In these forecasts, the probability of the event, over some lead time, is estimated based on some model simulations or predictive indicators. By issuing probabilistic forecasts, agencies may communicate the uncertainty in the event occurring. Assuming that the assigned probability of the event is correct, which is referred to as a reliable forecast, the end user may perform some risk management based on the potential damages resulting from the event. Alternatively, an unreliable forecast may give false impressions of the actual risk, leading to improper decision making when protecting resources from extreme events. Due to this requisite for reliable forecasts to perform effective risk management, this study takes a renewed look at reliability assessment in event forecasts. Illustrative experiments will be presented, showing deficiencies in the commonly available approaches (Brier Score, Reliability Diagram). Overall, it is shown that the conventional reliability assessment techniques do not maximize the ability to distinguish between a reliable and unreliable forecast. In this regard, a theoretical formulation of the probabilistic event forecast verification framework will be presented. From this analysis, hypothesis testing with the Poisson-Binomial distribution is the most exact model available for the verification framework, and therefore maximizes one's ability to distinguish between a reliable and unreliable forecast. Application of this verification system was also examined within a real forecasting case study, highlighting the additional statistical power provided with the use of the Poisson-Binomial distribution.

Tidal Triggering of Microearthquakes Over an Eruption Cycle at 9°50'N East Pacific Rise

NASA Astrophysics Data System (ADS)

Tan, Yen Joe; Tolstoy, Maya; Waldhauser, Felix; Bohnenstiehl, DelWayne R.

2018-02-01

Studies have found that earthquake timing often correlates with tides at mid-ocean ridges and some terrestrial settings. Studies have also suggested that tidal triggering may preferentially happen when a region is critically stressed, making it a potential tool to forecast earthquakes and volcanic eruptions. We examine tidal triggering of ˜100,000 microearthquakes near 9°50'N East Pacific Rise recorded between October 2003 and January 2007, which encompasses an eruption in January 2006. This allows us to look at how tidal triggering signal varies over an eruption cycle to examine its utility as a forecasting tool. We find that tidal triggering signal is strong but does not vary systematically in the 2+ years leading up to the eruption. However, tidal triggering signal disappears immediately posteruption. Our findings suggest that tidal triggering variation may not be useful for forecasting mid-ocean ridge eruptions over a 2+ year timescale but might be useful over a longer timescale.

Using Bayesian Model Averaging (BMA) to calibrate probabilistic surface temperature forecasts over Iran

NASA Astrophysics Data System (ADS)

Soltanzadeh, I.; Azadi, M.; Vakili, G. A.

2011-07-01

Using Bayesian Model Averaging (BMA), an attempt was made to obtain calibrated probabilistic numerical forecasts of 2-m temperature over Iran. The ensemble employs three limited area models (WRF, MM5 and HRM), with WRF used with five different configurations. Initial and boundary conditions for MM5 and WRF are obtained from the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS) and for HRM the initial and boundary conditions come from analysis of Global Model Europe (GME) of the German Weather Service. The resulting ensemble of seven members was run for a period of 6 months (from December 2008 to May 2009) over Iran. The 48-h raw ensemble outputs were calibrated using BMA technique for 120 days using a 40 days training sample of forecasts and relative verification data. The calibrated probabilistic forecasts were assessed using rank histogram and attribute diagrams. Results showed that application of BMA improved the reliability of the raw ensemble. Using the weighted ensemble mean forecast as a deterministic forecast it was found that the deterministic-style BMA forecasts performed usually better than the best member's deterministic forecast.

Communicating uncertainty in hydrological forecasts: mission impossible?

NASA Astrophysics Data System (ADS)

Ramos, Maria-Helena; Mathevet, Thibault; Thielen, Jutta; Pappenberger, Florian

2010-05-01

Cascading uncertainty in meteo-hydrological modelling chains for forecasting and integrated flood risk assessment is an essential step to improve the quality of hydrological forecasts. Although the best methodology to quantify the total predictive uncertainty in hydrology is still debated, there is a common agreement that one must avoid uncertainty misrepresentation and miscommunication, as well as misinterpretation of information by users. Several recent studies point out that uncertainty, when properly explained and defined, is no longer unwelcome among emergence response organizations, users of flood risk information and the general public. However, efficient communication of uncertain hydro-meteorological forecasts is far from being a resolved issue. This study focuses on the interpretation and communication of uncertain hydrological forecasts based on (uncertain) meteorological forecasts and (uncertain) rainfall-runoff modelling approaches to decision-makers such as operational hydrologists and water managers in charge of flood warning and scenario-based reservoir operation. An overview of the typical flow of uncertainties and risk-based decisions in hydrological forecasting systems is presented. The challenges related to the extraction of meaningful information from probabilistic forecasts and the test of its usefulness in assisting operational flood forecasting are illustrated with the help of two case-studies: 1) a study on the use and communication of probabilistic flood forecasting within the European Flood Alert System; 2) a case-study on the use of probabilistic forecasts by operational forecasters from the hydroelectricity company EDF in France. These examples show that attention must be paid to initiatives that promote or reinforce the active participation of expert forecasters in the forecasting chain. The practice of face-to-face forecast briefings, focusing on sharing how forecasters interpret, describe and perceive the model output forecasted scenarios, is essential. We believe that the efficient communication of uncertainty in hydro-meteorological forecasts is not a mission impossible. Questions remaining unanswered in probabilistic hydrological forecasting should not neutralize the goal of such a mission, and the suspense kept should instead act as a catalyst for overcoming the remaining challenges.

A probabilistic drought forecasting framework: A combined dynamical and statistical approach

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yan, Hongxiang; Moradkhani, Hamid; Zarekarizi, Mahkameh

In order to improve drought forecasting skill, this study develops a probabilistic drought forecasting framework comprised of dynamical and statistical modeling components. The novelty of this study is to seek the use of data assimilation to quantify initial condition uncertainty with the Monte Carlo ensemble members, rather than relying entirely on the hydrologic model or land surface model to generate a single deterministic initial condition, as currently implemented in the operational drought forecasting systems. Next, the initial condition uncertainty is quantified through data assimilation and coupled with a newly developed probabilistic drought forecasting model using a copula function. The initialmore » condition at each forecast start date are sampled from the data assimilation ensembles for forecast initialization. Finally, seasonal drought forecasting products are generated with the updated initial conditions. This study introduces the theory behind the proposed drought forecasting system, with an application in Columbia River Basin, Pacific Northwest, United States. Results from both synthetic and real case studies suggest that the proposed drought forecasting system significantly improves the seasonal drought forecasting skills and can facilitate the state drought preparation and declaration, at least three months before the official state drought declaration.« less

Understanding and forecasting phreatic eruptions driven by magmatic degassing

NASA Astrophysics Data System (ADS)

Stix, John; de Moor, J. Maarten

2018-05-01

This paper examines phreatic eruptions which are driven by inputs of magma and magmatic gas. We synthesize data from several significant phreatic systems, including two in Costa Rica (Turrialba and Poás) which are currently highly active and hazardous. We define two endmember types of phreatic eruptions, the first (type 1) in which a deeper hydrothermal system fed by magmatic gases is sealed and produces overpressure sufficient to drive explosive eruptions, and the second (type 2) where magmatic gases are supplied via open-vent degassing to a near-surface hydrothermal system, vaporizing liquid water which drives the phreatic eruptions. The surficial source of type 2 eruptions is characteristic, while the source depth of type 1 eruptions is commonly greater. Hence, type 1 eruptions tend to be more energetic than type 2 eruptions. The first type of eruption we term "phreato-vulcanian", and the second we term "phreato-surtseyan". Some systems (e.g., Ruapehu, Poás) can produce both type 1 and type 2 eruptions, and all systems can undergo sealing at various timescales. We examine a number of precursory signals which appear to be important in understanding and forecasting phreatic eruptions; these include very long period events, banded tremor, and gas ratios, in particular H2S/SO2 and CO2/SO2. We propose that if these datasets are carefully integrated during a monitoring program, it may be possible to accurately forecast phreatic eruptions.[Figure not available: see fulltext.

Exploring the calibration of a wind forecast ensemble for energy applications

NASA Astrophysics Data System (ADS)

Heppelmann, Tobias; Ben Bouallegue, Zied; Theis, Susanne

2015-04-01

In the German research project EWeLiNE, Deutscher Wetterdienst (DWD) and Fraunhofer Institute for Wind Energy and Energy System Technology (IWES) are collaborating with three German Transmission System Operators (TSO) in order to provide the TSOs with improved probabilistic power forecasts. Probabilistic power forecasts are derived from probabilistic weather forecasts, themselves derived from ensemble prediction systems (EPS). Since the considered raw ensemble wind forecasts suffer from underdispersiveness and bias, calibration methods are developed for the correction of the model bias and the ensemble spread bias. The overall aim is to improve the ensemble forecasts such that the uncertainty of the possible weather deployment is depicted by the ensemble spread from the first forecast hours. Additionally, the ensemble members after calibration should remain physically consistent scenarios. We focus on probabilistic hourly wind forecasts with horizon of 21 h delivered by the convection permitting high-resolution ensemble system COSMO-DE-EPS which has become operational in 2012 at DWD. The ensemble consists of 20 ensemble members driven by four different global models. The model area includes whole Germany and parts of Central Europe with a horizontal resolution of 2.8 km and a vertical resolution of 50 model levels. For verification we use wind mast measurements around 100 m height that corresponds to the hub height of wind energy plants that belong to wind farms within the model area. Calibration of the ensemble forecasts can be performed by different statistical methods applied to the raw ensemble output. Here, we explore local bivariate Ensemble Model Output Statistics at individual sites and quantile regression with different predictors. Applying different methods, we already show an improvement of ensemble wind forecasts from COSMO-DE-EPS for energy applications. In addition, an ensemble copula coupling approach transfers the time-dependencies of the raw ensemble to the calibrated ensemble. The calibrated wind forecasts are evaluated first with univariate probabilistic scores and additionally with diagnostics of wind ramps in order to assess the time-consistency of the calibrated ensemble members.

Against all odds -- Probabilistic forecasts and decision making

NASA Astrophysics Data System (ADS)

Liechti, Katharina; Zappa, Massimiliano

2015-04-01

In the city of Zurich (Switzerland) the setting is such that the damage potential due to flooding of the river Sihl is estimated to about 5 billion US dollars. The flood forecasting system that is used by the administration for decision making runs continuously since 2007. It has a time horizon of max. five days and operates at hourly time steps. The flood forecasting system includes three different model chains. Two of those are run by the deterministic NWP models COSMO-2 and COSMO-7 and one is driven by the probabilistic NWP COSMO-Leps. The model chains are consistent since February 2010, so five full years are available for the evaluation for the system. The system was evaluated continuously and is a very nice example to present the added value that lies in probabilistic forecasts. The forecasts are available on an online-platform to the decision makers. Several graphical representations of the forecasts and forecast-history are available to support decision making and to rate the current situation. The communication between forecasters and decision-makers is quite close. To put it short, an ideal situation. However, an event or better put a non-event in summer 2014 showed that the knowledge about the general superiority of probabilistic forecasts doesn't necessarily mean that the decisions taken in a specific situation will be based on that probabilistic forecast. Some years of experience allow gaining confidence in the system, both for the forecasters and for the decision-makers. Even if from the theoretical point of view the handling during crisis situation is well designed, a first event demonstrated that the dialog with the decision-makers still lacks of exercise during such situations. We argue, that a false alarm is a needed experience to consolidate real-time emergency procedures relying on ensemble predictions. A missed event would probably also fit, but, in our case, we are very happy not to report about this option.

Probability hazard map for future vent opening at Etna volcano (Sicily, Italy).

NASA Astrophysics Data System (ADS)

Brancato, Alfonso; Tusa, Giuseppina; Coltelli, Mauro; Proietti, Cristina

2014-05-01

Mount Etna is a composite stratovolcano located along the Ionian coast of eastern Sicily. The frequent flank eruptions occurrence (at an interval of years, mostly concentrated along the NE, S and W rift zones) lead to a high volcanic hazard that, linked with intense urbanization, poses a high volcanic risk. A long-term volcanic hazard assessment, mainly based on the past behaviour of the Etna volcano, is the basic tool for the evaluation of this risk. Then, a reliable forecast where the next eruption will occur is needed. A computer-assisted analysis and probabilistic evaluations will provide the relative map, thus allowing identification of the areas prone to the highest hazard. Based on these grounds, the use of a code such BET_EF (Bayesian Event Tree_Eruption Forecasting) showed that a suitable analysis can be explored (Selva et al., 2012). Following an analysis we are performing, a total of 6886 point-vents referring to the last 4.0 ka of Etna flank activity, and spread over an area of 744 km2 (divided into N=2976 squared cell, with side of 500 m), allowed us to estimate a pdf by applying a Gaussian kernel. The probability values represent a complete set of outcomes mutually exclusive and the relative sum is normalized to one over the investigated area; then, the basic assumptions of a Dirichlet distribution (the prior distribution set in the BET_EF code (Marzocchi et al., 2004, 2008)) still hold. One fundamental parameter is the the equivalent number of data, that depicts our confidence on the best guess probability. The BET_EF code also works with a likelihood function. This is modelled by a Multinomial distribution, with parameters representing the number of vents in each cell and the total number of past data (i.e. the 6886 point-vents). Given the grid of N cells, the final posterior distribution will be evaluated by multiplying the a priori Dirichlet probability distribution with the past data in each cell through the likelihood. The probability hazard map shows a tendency to concentrate along the NE and S rifts, as well as Valle del Bove, increasing the difference in probability between these areas and the rest of the volcano edifice. It is worthy notice that a higher significance is still evident along the W rift, even if not comparable with the ones of the above mentioned areas. References Marzocchi W., Sandri L., Gasparini P., Newhall C. and Boschi E.; 2004: Quantifying probabilities of volcanic events: The example of volcanic hazard at Mount Vesuvius, J. Geophys. Res., 109, B11201, doi:10.1029/2004JB00315U. Marzocchi W., Sandri, L. and Selva, J.; 2008: BET_EF: a probabilistic tool for long- and short-term eruption forecasting, Bull. Volcanol., 70, 623 - 632, doi: 10.1007/s00445-007-0157-y. Selva J., Orsi G., Di Vito M.A., Marzocchi W. And Sandri L.; 2012: Probability hazard mapfor future vent opening atthe Campi Flegrei caldera, Italy, Bull. Volcanol., 74, 497 - 510, doi: 10.1007/s00445-011-0528-2.

A parimutuel gambling perspective to compare probabilistic seismicity forecasts

NASA Astrophysics Data System (ADS)

Zechar, J. Douglas; Zhuang, Jiancang

2014-10-01

Using analogies to gaming, we consider the problem of comparing multiple probabilistic seismicity forecasts. To measure relative model performance, we suggest a parimutuel gambling perspective which addresses shortcomings of other methods such as likelihood ratio, information gain and Molchan diagrams. We describe two variants of the parimutuel approach for a set of forecasts: head-to-head, in which forecasts are compared in pairs, and round table, in which all forecasts are compared simultaneously. For illustration, we compare the 5-yr forecasts of the Regional Earthquake Likelihood Models experiment for M4.95+ seismicity in California.

Towards real-time eruption forecasting in the Auckland Volcanic Field: application of BET_EF during the New Zealand National Disaster Exercise `Ruaumoko'

NASA Astrophysics Data System (ADS)

Lindsay, Jan; Marzocchi, Warner; Jolly, Gill; Constantinescu, Robert; Selva, Jacopo; Sandri, Laura

2010-03-01

The Auckland Volcanic Field (AVF) is a young basaltic field that lies beneath the urban area of Auckland, New Zealand’s largest city. Over the past 250,000 years the AVF has produced at least 49 basaltic centers; the last eruption was only 600 years ago. In recognition of the high risk associated with a possible future eruption in Auckland, the New Zealand government ran Exercise Ruaumoko in March 2008, a test of New Zealand’s nation-wide preparedness for responding to a major disaster resulting from a volcanic eruption in Auckland City. The exercise scenario was developed in secret, and covered the period of precursory activity up until the eruption. During Exercise Ruaumoko we adapted a recently developed statistical code for eruption forecasting, namely BET_EF (Bayesian Event Tree for Eruption Forecasting), to independently track the unrest evolution and to forecast the most likely onset time, location and style of the initial phase of the simulated eruption. The code was set up before the start of the exercise by entering reliable information on the past history of the AVF as well as the monitoring signals expected in the event of magmatic unrest and an impending eruption. The average probabilities calculated by BET_EF during Exercise Ruaumoko corresponded well to the probabilities subjectively (and independently) estimated by the advising scientists (differences of few percentage units), and provided a sound forecast of the timing (before the event, the eruption probability reached 90%) and location of the eruption. This application of BET_EF to a volcanic field that has experienced no historical activity and for which otherwise limited prior information is available shows its versatility and potential usefulness as a tool to aid decision-making for a wide range of volcano types. Our near real-time application of BET_EF during Exercise Ruaumoko highlighted its potential to clarify and possibly optimize decision-making procedures in a future AVF eruption crisis, and as a rational starting point for discussions in a scientific advisory group. It also stimulated valuable scientific discussion around how a future AVF eruption might progress, and highlighted areas of future volcanological research that would reduce epistemic uncertainties through the development of better input models.

Statistical forecasting of repetitious dome failures during the waning eruption of Redoubt Volcano, Alaska, February-April 1990

USGS Publications Warehouse

Page, R.A.; Lahr, J.C.; Chouet, B.A.; Power, J.A.; Stephens, C.D.

1994-01-01

The waning phase of the 1989-1990 eruption of Redoubt Volcano in the Cook Inlet region of south-central Alaska comprised a quasi-regular pattern of repetitious dome growth and destruction that lasted from February 15 to late April 1990. The dome failures produced ash plumes hazardous to airline traffic. In response to this hazard, the Alaska Volcano Observatory sought to forecast these ash-producing events using two approaches. One approach built on early successes in issuing warnings before major eruptions on December 14, 1989 and January 2, 1990. These warnings were based largely on changes in seismic activity related to the occurrence of precursory swarms of long-period seismic events. The search for precursory swarms of long-period seismicity was continued through the waning phase of the eruption and led to warnings before tephra eruptions on March 23 and April 6. The observed regularity of dome failures after February 15 suggested that a statistical forecasting method based on a constant-rate failure model might also be successful. The first statistical forecast was issued on March 16 after seven events had occurred, at an average interval of 4.5 days. At this time, the interval between dome failures abruptly lengthened. Accordingly, the forecast was unsuccessful and further forecasting was suspended until the regularity of subsequent failures could be confirmed. Statistical forecasting resumed on April 12, after four dome failure episodes separated by an average of 7.8 days. One dome failure (April 15) was successfully forecast using a 70% confidence window, and a second event (April 21) was narrowly missed before the end of the activity. The cessation of dome failures after April 21 resulted in a concluding false alarm. Although forecasting success during the eruption was limited, retrospective analysis shows that early and consistent application of the statistical method using a constant-rate failure model and a 90% confidence window could have yielded five successful forecasts and two false alarms; no events would have been missed. On closer examination, the intervals between successive dome failures are not uniform but tend to increase with time. This increase attests to the continuous, slowly decreasing supply of magma to the surface vent during the waning phase of the eruption. The domes formed in a precarious position in a breach in the summit crater rim where they were susceptible to gravitational collapse. The instability of the February 15-April 21 domes relative to the earlier domes is attributed to reaming the lip of the vent by a laterally directed explosion during the major dome-destroying eruption of February 15, a process which would leave a less secure foundation for subsequent domes. ?? 1994.

Integrating geological and geophysical data to improve probabilistic hazard forecasting of Arabian Shield volcanism

NASA Astrophysics Data System (ADS)

Runge, Melody G.; Bebbington, Mark S.; Cronin, Shane J.; Lindsay, Jan M.; Moufti, Mohammed R.

2016-02-01

During probabilistic volcanic hazard analysis of volcanic fields, a greater variety of spatial data on crustal features should help improve forecasts of future vent locations. Without further examination, however, geophysical estimations of crustal or other features may be non-informative. Here, we present a new, robust, non-parametric method to quantitatively determine the existence of any relationship between natural phenomena (e.g., volcanic eruptions) and a variety of geophysical data. This provides a new validation tool for incorporating a range of potentially hazard-diagnostic observable data into recurrence rate estimates and hazard analyses. Through this study it is shown that the location of Cenozoic volcanic fields across the Arabian Shield appear to be related to locations of major and minor faults, at higher elevations, and regions where gravity anomaly values were between - 125 mGal and 0 mGal. These findings support earlier hypotheses that the western shield uplift was related to Cenozoic volcanism. At the harrat (volcanic field)-scale, higher vent density regions are related to both elevation and gravity anomaly values. A by-product of this work is the collection of existing data on the volcanism across Saudi Arabia, with all vent locations provided herein, as well as updated maps for Harrats Kura, Khaybar, Ithnayn, Kishb, and Rahat. This work also highlights the potential dangers of assuming relationships between observed data and the occurrence of a natural phenomenon without quantitative assessment or proper consideration of the effects of data resolution.

Forecasting deflation, intrusion and eruption at inflating volcanoes

NASA Astrophysics Data System (ADS)

Blake, Stephen; Cortés, Joaquín A.

2018-01-01

A principal goal of volcanology is to successfully forecast the start of volcanic eruptions. This paper introduces a general forecasting method, which relies on a stream of monitoring data and a statistical description of a given threshold criterion for an eruption to start. Specifically we investigate the timing of intrusive and eruptive events at inflating volcanoes. The gradual inflation of the ground surface is a well-known phenomenon at many volcanoes and is attributable to pressurised magma accumulating within a shallow chamber. Inflation usually culminates in a rapid deflation event caused by magma escaping from the chamber to produce a shallow intrusion and, in some cases, a volcanic eruption. We show that the ground elevation during 15 inflation periods at Krafla volcano, Iceland, increased with time towards a limiting value by following a decaying exponential with characteristic timescale τ. The available data for Krafla, Kilauea and Mauna Loa volcanoes show that the duration of inflation (t*) is approximately equal to τ. The distribution of t* / τ values follows a log-logistic distribution in which the central 60% of the data lie between 0.99

Probabilistic-numerical assessment of pyroclastic current hazard at Campi Flegrei and Naples city: Multi-VEI scenarios as a tool for "full-scale" risk management.

PubMed

Mastrolorenzo, Giuseppe; Palladino, Danilo M; Pappalardo, Lucia; Rossano, Sergio

2017-01-01

The Campi Flegrei volcanic field (Italy) poses very high risk to the highly urbanized Neapolitan area. Eruptive history was dominated by explosive activity producing pyroclastic currents (hereon PCs; acronym for Pyroclastic Currents) ranging in scale from localized base surges to regional flows. Here we apply probabilistic numerical simulation approaches to produce PC hazard maps, based on a comprehensive spectrum of flow properties and vent locations. These maps are incorporated in a Geographic Information System (GIS) and provide all probable Volcanic Explosivity Index (VEI) scenarios from different source vents in the caldera, relevant for risk management planning. For each VEI scenario, we report the conditional probability for PCs (i.e., the probability for a given area to be affected by the passage of PCs in case of a PC-forming explosive event) and related dynamic pressure. Model results indicate that PCs from VEI<4 events would be confined within the Campi Flegrei caldera, PC propagation being impeded by the northern and eastern caldera walls. Conversely, PCs from VEI 4-5 events could invade a wide area beyond the northern caldera rim, as well as part of the Naples metropolitan area to the east. A major controlling factor of PC dispersal is represented by the location of the vent area. PCs from the potentially largest eruption scenarios (analogous to the ~15 ka, VEI 6 Neapolitan Yellow Tuff or even the ~39 ka, VEI 7 Campanian Ignimbrite extreme event) would affect a large part of the Campanian Plain to the north and the city of Naples to the east. Thus, in case of renewal of eruptive activity at Campi Flegrei, up to 3 million people will be potentially exposed to volcanic hazard, pointing out the urgency of an emergency plan. Considering the present level of uncertainty in forecasting the future eruption type, size and location (essentially based on statistical analysis of previous activity), we suggest that appropriate planning measures should face at least the VEI 5 reference scenario (at least 2 occurrences documented in the last 10 ka).

Probabilistic-numerical assessment of pyroclastic current hazard at Campi Flegrei and Naples city: Multi-VEI scenarios as a tool for “full-scale” risk management

PubMed Central

Mastrolorenzo, Giuseppe; Palladino, Danilo M.; Pappalardo, Lucia; Rossano, Sergio

2017-01-01

The Campi Flegrei volcanic field (Italy) poses very high risk to the highly urbanized Neapolitan area. Eruptive history was dominated by explosive activity producing pyroclastic currents (hereon PCs; acronym for Pyroclastic Currents) ranging in scale from localized base surges to regional flows. Here we apply probabilistic numerical simulation approaches to produce PC hazard maps, based on a comprehensive spectrum of flow properties and vent locations. These maps are incorporated in a Geographic Information System (GIS) and provide all probable Volcanic Explosivity Index (VEI) scenarios from different source vents in the caldera, relevant for risk management planning. For each VEI scenario, we report the conditional probability for PCs (i.e., the probability for a given area to be affected by the passage of PCs in case of a PC-forming explosive event) and related dynamic pressure. Model results indicate that PCs from VEI<4 events would be confined within the Campi Flegrei caldera, PC propagation being impeded by the northern and eastern caldera walls. Conversely, PCs from VEI 4–5 events could invade a wide area beyond the northern caldera rim, as well as part of the Naples metropolitan area to the east. A major controlling factor of PC dispersal is represented by the location of the vent area. PCs from the potentially largest eruption scenarios (analogous to the ~15 ka, VEI 6 Neapolitan Yellow Tuff or even the ~39 ka, VEI 7 Campanian Ignimbrite extreme event) would affect a large part of the Campanian Plain to the north and the city of Naples to the east. Thus, in case of renewal of eruptive activity at Campi Flegrei, up to 3 million people will be potentially exposed to volcanic hazard, pointing out the urgency of an emergency plan. Considering the present level of uncertainty in forecasting the future eruption type, size and location (essentially based on statistical analysis of previous activity), we suggest that appropriate planning measures should face at least the VEI 5 reference scenario (at least 2 occurrences documented in the last 10 ka). PMID:29020018

Probabilistic Solar Wind and Geomagnetic Forecasting Using an Analogue Ensemble or "Similar Day" Approach

NASA Astrophysics Data System (ADS)

Owens, M. J.; Riley, P.; Horbury, T. S.

2017-05-01

Effective space-weather prediction and mitigation requires accurate forecasting of near-Earth solar-wind conditions. Numerical magnetohydrodynamic models of the solar wind, driven by remote solar observations, are gaining skill at forecasting the large-scale solar-wind features that give rise to near-Earth variations over days and weeks. There remains a need for accurate short-term (hours to days) solar-wind forecasts, however. In this study we investigate the analogue ensemble (AnEn), or "similar day", approach that was developed for atmospheric weather forecasting. The central premise of the AnEn is that past variations that are analogous or similar to current conditions can be used to provide a good estimate of future variations. By considering an ensemble of past analogues, the AnEn forecast is inherently probabilistic and provides a measure of the forecast uncertainty. We show that forecasts of solar-wind speed can be improved by considering both speed and density when determining past analogues, whereas forecasts of the out-of-ecliptic magnetic field [BN] are improved by also considering the in-ecliptic magnetic-field components. In general, the best forecasts are found by considering only the previous 6 - 12 hours of observations. Using these parameters, the AnEn provides a valuable probabilistic forecast for solar-wind speed, density, and in-ecliptic magnetic field over lead times from a few hours to around four days. For BN, which is central to space-weather disturbance, the AnEn only provides a valuable forecast out to around six to seven hours. As the inherent predictability of this parameter is low, this is still likely a marked improvement over other forecast methods. We also investigate the use of the AnEn in forecasting geomagnetic indices Dst and Kp. The AnEn provides a valuable probabilistic forecast of both indices out to around four days. We outline a number of future improvements to AnEn forecasts of near-Earth solar-wind and geomagnetic conditions.

Generating short-term probabilistic wind power scenarios via nonparametric forecast error density estimators: Generating short-term probabilistic wind power scenarios via nonparametric forecast error density estimators

DOE PAGES

Staid, Andrea; Watson, Jean -Paul; Wets, Roger J. -B.; ...

2017-07-11

Forecasts of available wind power are critical in key electric power systems operations planning problems, including economic dispatch and unit commitment. Such forecasts are necessarily uncertain, limiting the reliability and cost effectiveness of operations planning models based on a single deterministic or “point” forecast. A common approach to address this limitation involves the use of a number of probabilistic scenarios, each specifying a possible trajectory of wind power production, with associated probability. We present and analyze a novel method for generating probabilistic wind power scenarios, leveraging available historical information in the form of forecasted and corresponding observed wind power timemore » series. We estimate non-parametric forecast error densities, specifically using epi-spline basis functions, allowing us to capture the skewed and non-parametric nature of error densities observed in real-world data. We then describe a method to generate probabilistic scenarios from these basis functions that allows users to control for the degree to which extreme errors are captured.We compare the performance of our approach to the current state-of-the-art considering publicly available data associated with the Bonneville Power Administration, analyzing aggregate production of a number of wind farms over a large geographic region. Finally, we discuss the advantages of our approach in the context of specific power systems operations planning problems: stochastic unit commitment and economic dispatch. Here, our methodology is embodied in the joint Sandia – University of California Davis Prescient software package for assessing and analyzing stochastic operations strategies.« less

Probabilistic versus deterministic skill in predicting the western North Pacific-East Asian summer monsoon variability with multimodel ensembles

NASA Astrophysics Data System (ADS)

Yang, Xiu-Qun; Yang, Dejian; Xie, Qian; Zhang, Yaocun; Ren, Xuejuan; Tang, Youmin

2017-04-01

Based on historical forecasts of three quasi-operational multi-model ensemble (MME) systems, this study assesses the superiority of coupled MME over contributing single-model ensembles (SMEs) and over uncoupled atmospheric MME in predicting the Western North Pacific-East Asian summer monsoon variability. The probabilistic and deterministic forecast skills are measured by Brier skill score (BSS) and anomaly correlation (AC), respectively. A forecast-format dependent MME superiority over SMEs is found. The probabilistic forecast skill of the MME is always significantly better than that of each SME, while the deterministic forecast skill of the MME can be lower than that of some SMEs. The MME superiority arises from both the model diversity and the ensemble size increase in the tropics, and primarily from the ensemble size increase in the subtropics. The BSS is composed of reliability and resolution, two attributes characterizing probabilistic forecast skill. The probabilistic skill increase of the MME is dominated by the dramatic improvement in reliability, while resolution is not always improved, similar to AC. A monotonic resolution-AC relationship is further found and qualitatively explained, whereas little relationship can be identified between reliability and AC. It is argued that the MME's success in improving the reliability arises from an effective reduction of the overconfidence in forecast distributions. Moreover, it is examined that the seasonal predictions with coupled MME are more skillful than those with the uncoupled atmospheric MME forced by persisting sea surface temperature (SST) anomalies, since the coupled MME has better predicted the SST anomaly evolution in three key regions.

Generating short-term probabilistic wind power scenarios via nonparametric forecast error density estimators: Generating short-term probabilistic wind power scenarios via nonparametric forecast error density estimators

DOE Office of Scientific and Technical Information (OSTI.GOV)

Staid, Andrea; Watson, Jean -Paul; Wets, Roger J. -B.

Forecasts of available wind power are critical in key electric power systems operations planning problems, including economic dispatch and unit commitment. Such forecasts are necessarily uncertain, limiting the reliability and cost effectiveness of operations planning models based on a single deterministic or “point” forecast. A common approach to address this limitation involves the use of a number of probabilistic scenarios, each specifying a possible trajectory of wind power production, with associated probability. We present and analyze a novel method for generating probabilistic wind power scenarios, leveraging available historical information in the form of forecasted and corresponding observed wind power timemore » series. We estimate non-parametric forecast error densities, specifically using epi-spline basis functions, allowing us to capture the skewed and non-parametric nature of error densities observed in real-world data. We then describe a method to generate probabilistic scenarios from these basis functions that allows users to control for the degree to which extreme errors are captured.We compare the performance of our approach to the current state-of-the-art considering publicly available data associated with the Bonneville Power Administration, analyzing aggregate production of a number of wind farms over a large geographic region. Finally, we discuss the advantages of our approach in the context of specific power systems operations planning problems: stochastic unit commitment and economic dispatch. Here, our methodology is embodied in the joint Sandia – University of California Davis Prescient software package for assessing and analyzing stochastic operations strategies.« less

Beyond eruptive scenarios: assessing tephra fallout hazard from Neapolitan volcanoes

PubMed Central

Sandri, Laura; Costa, Antonio; Selva, Jacopo; Tonini, Roberto; Macedonio, Giovanni; Folch, Arnau; Sulpizio, Roberto

2016-01-01

Assessment of volcanic hazards is necessary for risk mitigation. Typically, hazard assessment is based on one or a few, subjectively chosen representative eruptive scenarios, which use a specific combination of eruptive sizes and intensities to represent a particular size class of eruption. While such eruptive scenarios use a range of representative members to capture a range of eruptive sizes and intensities in order to reflect a wider size class, a scenario approach neglects to account for the intrinsic variability of volcanic eruptions, and implicitly assumes that inter-class size variability (i.e. size difference between different eruptive size classes) dominates over intra-class size variability (i.e. size difference within an eruptive size class), the latter of which is treated as negligible. So far, no quantitative study has been undertaken to verify such an assumption. Here, we adopt a novel Probabilistic Volcanic Hazard Analysis (PVHA) strategy, which accounts for intrinsic eruptive variabilities, to quantify the tephra fallout hazard in the Campania area. We compare the results of the new probabilistic approach with the classical scenario approach. The results allow for determining whether a simplified scenario approach can be considered valid, and for quantifying the bias which arises when full variability is not accounted for. PMID:27067389

Detecting and Characterizing Repeating Earthquake Sequences During Volcanic Eruptions

NASA Astrophysics Data System (ADS)

Tepp, G.; Haney, M. M.; Wech, A.

2017-12-01

A major challenge in volcano seismology is forecasting eruptions. Repeating earthquake sequences often precede volcanic eruptions or lava dome activity, providing an opportunity for short-term eruption forecasting. Automatic detection of these sequences can lead to timely eruption notification and aid in continuous monitoring of volcanic systems. However, repeating earthquake sequences may also occur after eruptions or along with magma intrusions that do not immediately lead to an eruption. This additional challenge requires a better understanding of the processes involved in producing these sequences to distinguish those that are precursory. Calculation of the inverse moment rate and concepts from the material failure forecast method can lead to such insights. The temporal evolution of the inverse moment rate is observed to differ for precursory and non-precursory sequences, and multiple earthquake sequences may occur concurrently. These observations suggest that sequences may occur in different locations or through different processes. We developed an automated repeating earthquake sequence detector and near real-time alarm to send alerts when an in-progress sequence is identified. Near real-time inverse moment rate measurements can further improve our ability to forecast eruptions by allowing for characterization of sequences. We apply the detector to eruptions of two Alaskan volcanoes: Bogoslof in 2016-2017 and Redoubt Volcano in 2009. The Bogoslof eruption produced almost 40 repeating earthquake sequences between its start in mid-December 2016 and early June 2017, 21 of which preceded an explosive eruption, and 2 sequences in the months before eruptive activity. Three of the sequences occurred after the implementation of the alarm in late March 2017 and successfully triggered alerts. The nearest seismometers to Bogoslof are over 45 km away, requiring a detector that can work with few stations and a relatively low signal-to-noise ratio. During the Redoubt eruption, earthquake sequences were observed in the months leading up to the eruptive activity beginning in March 2009 as well as immediately preceding 7 of the 19 explosive events. In contrast to Bogoslof, Redoubt has a local monitoring network which allows for better detection and more detailed analysis of the repeating earthquake sequences.

WIPCast: Probabilistic Forecasting for Aviation Decision Aid Applications

DTIC Science & Technology

2011-06-01

traders, or families planning an outing – manage weather-related risk. By quantifying risk , probabilistic forecasting enables optimization of actions via...confidence interval to the user’s risk tolerance helps drive highly effective and innovative decision support mechanisms for visually quantifying risk for

Probabilistic forecasting for extreme NO2 pollution episodes.

PubMed

Aznarte, José L

2017-10-01

In this study, we investigate the convenience of quantile regression to predict extreme concentrations of NO 2 . Contrarily to the usual point-forecasting, where a single value is forecast for each horizon, probabilistic forecasting through quantile regression allows for the prediction of the full probability distribution, which in turn allows to build models specifically fit for the tails of this distribution. Using data from the city of Madrid, including NO 2 concentrations as well as meteorological measures, we build models that predict extreme NO 2 concentrations, outperforming point-forecasting alternatives, and we prove that the predictions are accurate, reliable and sharp. Besides, we study the relative importance of the independent variables involved, and show how the important variables for the median quantile are different than those important for the upper quantiles. Furthermore, we present a method to compute the probability of exceedance of thresholds, which is a simple and comprehensible manner to present probabilistic forecasts maximizing their usefulness. Copyright © 2017 Elsevier Ltd. All rights reserved.

Forecasting the Onset Time of Volcanic Eruptions Using Ground Deformation Data

NASA Astrophysics Data System (ADS)

Blake, S.; Cortes, J. A.

2016-12-01

The pre-eruptive inflation of the ground surface is a well-known phenomenon at many volcanoes. In a number of intensively studied cases, elevation and/or radial tilt increase with time (t) towards a limiting value by following a decaying exponential with characteristic timescale τ (Kilauea and Mauna Loa: Dvorak and Okamura 1987, Lengliné et al., 2008) or, after sufficiently long times, by following the sum of two such functions such that two timescales, τ1 and τ2, are required to describe the temporal pattern of inflation (Axial Seamount: Nooner and Chadwick, 2009). We have used the Levenberg-Marquardt non-linear fit algorithm to analyse data for 18 inflation periods at Krafla volcano, Iceland, (Björnsson and Eysteinsson, 1998) and found the same functional relationship. Pooling all of the available data from 25 eruptions at 4 volcanoes shows that the duration of inflation before an eruption or shallow intrusion (t*) is comparable to τ (or the longer of τ1 and τ2) and follows an almost 1:1 linear relationship (r2 0.8). We also find that this scaling is replicated by Monte Carlo simulations of physics-based forward models of hydraulically connected dual magma chamber systems which erupt when the chamber pressure reaches a threshold value. These results lead to a new forecasting method which we describe and assess here: if τ can be constrained during an on-going inflation period, then the statistical distribution of t*/τ values calibrated from other pre-eruptive inflation periods allows the probability of an eruption starting before (or after) a specified time to be estimated. The time at which there is a specified probability of an eruption starting can also be forecast. These approaches rely on fitting deformation data up to time t in order to obtain τ(t) which is then used to forecast t*. Forecasts can be updated after each new deformation measurement.

Real Time Volcanic Cloud Products and Predictions for Aviation Alerts

NASA Technical Reports Server (NTRS)

Krotkov, Nickolay A.; Habib, Shahid; da Silva, Arlindo; Hughes, Eric; Yang, Kai; Brentzel, Kelvin; Seftor, Colin; Li, Jason Y.; Schneider, David; Guffanti, Marianne;

Forecasting eruption size: what we know, what we don't know

NASA Astrophysics Data System (ADS)

Papale, Paolo

2017-04-01

Any eruption forecast includes an evaluation of the expected size of the forthcoming eruption, usually expressed as the probability associated to given size classes. Such evaluation is mostly based on the previous volcanic history at the specific volcano, or it is referred to a broader class of volcanoes constituting "analogues" of the one under specific consideration. In any case, use of knowledge from past eruptions implies considering the completeness of the reference catalogue, and most importantly, the existence of systematic biases in the catalogue, that may affect probability estimates and translate into biased volcanic hazard forecasts. An analysis of existing catalogues, with major reference to the catalogue from the Smithsonian Global Volcanism Program, suggests that systematic biases largely dominate at global, regional and local scale: volcanic histories reconstructed at individual volcanoes, often used as a reference for volcanic hazard forecasts, are the result of systematic loss of information with time and poor sample representativeness. That situation strictly requires the use of techniques to complete existing catalogues, as well as careful consideration of the uncertainties deriving from inadequate knowledge and model-dependent data elaboration. A reconstructed global eruption size distribution, obtained by merging information from different existing catalogues, shows a mode in the VEI 1-2 range, <0.1% incidence of eruptions with VEI 7 or larger, and substantial uncertainties associated with individual VEI frequencies. Even larger uncertainties are expected to derive from application to individual volcanoes or classes of analogue volcanoes, suggesting large to very large uncertainties associated to volcanic hazard forecasts virtually at any individual volcano worldwide.

Probabilistic forecasting of extreme weather events based on extreme value theory

NASA Astrophysics Data System (ADS)

Van De Vyver, Hans; Van Schaeybroeck, Bert

2016-04-01

Extreme events in weather and climate such as high wind gusts, heavy precipitation or extreme temperatures are commonly associated with high impacts on both environment and society. Forecasting extreme weather events is difficult, and very high-resolution models are needed to describe explicitly extreme weather phenomena. A prediction system for such events should therefore preferably be probabilistic in nature. Probabilistic forecasts and state estimations are nowadays common in the numerical weather prediction community. In this work, we develop a new probabilistic framework based on extreme value theory that aims to provide early warnings up to several days in advance. We consider the combined events when an observation variable Y (for instance wind speed) exceeds a high threshold y and its corresponding deterministic forecasts X also exceeds a high forecast threshold y. More specifically two problems are addressed:} We consider pairs (X,Y) of extreme events where X represents a deterministic forecast, and Y the observation variable (for instance wind speed). More specifically two problems are addressed: Given a high forecast X=x_0, what is the probability that Y>y? In other words: provide inference on the conditional probability: [ Pr{Y>y|X=x_0}. ] Given a probabilistic model for Problem 1, what is the impact on the verification analysis of extreme events. These problems can be solved with bivariate extremes (Coles, 2001), and the verification analysis in (Ferro, 2007). We apply the Ramos and Ledford (2009) parametric model for bivariate tail estimation of the pair (X,Y). The model accommodates different types of extremal dependence and asymmetry within a parsimonious representation. Results are presented using the ensemble reforecast system of the European Centre of Weather Forecasts (Hagedorn, 2008). Coles, S. (2001) An Introduction to Statistical modelling of Extreme Values. Springer-Verlag.Ferro, C.A.T. (2007) A probability model for verifying deterministic forecasts of extreme events. Wea. Forecasting {22}, 1089-1100.Hagedorn, R. (2008) Using the ECMWF reforecast dataset to calibrate EPS forecasts. ECMWF Newsletter, {117}, 8-13.Ramos, A., Ledford, A. (2009) A new class of models for bivariate joint tails. J.R. Statist. Soc. B {71}, 219-241.

Evaluation of ensemble forecast uncertainty using a new proper score: application to medium-range and seasonal forecasts

NASA Astrophysics Data System (ADS)

Christensen, Hannah; Moroz, Irene; Palmer, Tim

2015-04-01

Forecast verification is important across scientific disciplines as it provides a framework for evaluating the performance of a forecasting system. In the atmospheric sciences, probabilistic skill scores are often used for verification as they provide a way of unambiguously ranking the performance of different probabilistic forecasts. In order to be useful, a skill score must be proper -- it must encourage honesty in the forecaster, and reward forecasts which are reliable and which have good resolution. A new score, the Error-spread Score (ES), is proposed which is particularly suitable for evaluation of ensemble forecasts. It is formulated with respect to the moments of the forecast. The ES is confirmed to be a proper score, and is therefore sensitive to both resolution and reliability. The ES is tested on forecasts made using the Lorenz '96 system, and found to be useful for summarising the skill of the forecasts. The European Centre for Medium-Range Weather Forecasts (ECMWF) ensemble prediction system (EPS) is evaluated using the ES. Its performance is compared to a perfect statistical probabilistic forecast -- the ECMWF high resolution deterministic forecast dressed with the observed error distribution. This generates a forecast that is perfectly reliable if considered over all time, but which does not vary from day to day with the predictability of the atmospheric flow. The ES distinguishes between the dynamically reliable EPS forecasts and the statically reliable dressed deterministic forecasts. Other skill scores are tested and found to be comparatively insensitive to this desirable forecast quality. The ES is used to evaluate seasonal range ensemble forecasts made with the ECMWF System 4. The ensemble forecasts are found to be skilful when compared with climatological or persistence forecasts, though this skill is dependent on region and time of year.

The longevity of lava dome eruptions

NASA Astrophysics Data System (ADS)

Wolpert, Robert L.; Ogburn, Sarah E.; Calder, Eliza S.

2016-02-01

Understanding the duration of past, ongoing, and future volcanic eruptions is an important scientific goal and a key societal need. We present a new methodology for forecasting the duration of ongoing and future lava dome eruptions based on a database (DomeHaz) recently compiled by the authors. The database includes duration and composition for 177 such eruptions, with "eruption" defined as the period encompassing individual episodes of dome growth along with associated quiescent periods during which extrusion pauses but unrest continues. In a key finding, we show that probability distributions for dome eruption durations are both heavy tailed and composition dependent. We construct objective Bayesian statistical models featuring heavy-tailed Generalized Pareto distributions with composition-specific parameters to make forecasts about the durations of new and ongoing eruptions that depend on both eruption duration to date and composition. Our Bayesian predictive distributions reflect both uncertainty about model parameter values (epistemic uncertainty) and the natural variability of the geologic processes (aleatoric uncertainty). The results are illustrated by presenting likely trajectories for 14 dome-building eruptions ongoing in 2015. Full representation of the uncertainty is presented for two key eruptions, Soufriére Hills Volcano in Montserrat (10-139 years, median 35 years) and Sinabung, Indonesia (1-17 years, median 4 years). Uncertainties are high but, importantly, quantifiable. This work provides for the first time a quantitative and transferable method and rationale on which to base long-term planning decisions for lava dome-forming volcanoes, with wide potential use and transferability to forecasts of other types of eruptions and other adverse events across the geohazard spectrum.

Climatological Observations for Maritime Prediction and Analysis Support Service (COMPASS)

NASA Astrophysics Data System (ADS)

OConnor, A.; Kirtman, B. P.; Harrison, S.; Gorman, J.

2016-02-01

Current US Navy forecasting systems cannot easily incorporate extended-range forecasts that can improve mission readiness and effectiveness; ensure safety; and reduce cost, labor, and resource requirements. If Navy operational planners had systems that incorporated these forecasts, they could plan missions using more reliable and longer-term weather and climate predictions. Further, using multi-model forecast ensembles instead of single forecasts would produce higher predictive performance. Extended-range multi-model forecast ensembles, such as those available in the North American Multi-Model Ensemble (NMME), are ideal for system integration because of their high skill predictions; however, even higher skill predictions can be produced if forecast model ensembles are combined correctly. While many methods for weighting models exist, the best method in a given environment requires expert knowledge of the models and combination methods.We present an innovative approach that uses machine learning to combine extended-range predictions from multi-model forecast ensembles and generate a probabilistic forecast for any region of the globe up to 12 months in advance. Our machine-learning approach uses 30 years of hindcast predictions to learn patterns of forecast model successes and failures. Each model is assigned a weight for each environmental condition, 100 km2 region, and day given any expected environmental information. These weights are then applied to the respective predictions for the region and time of interest to effectively stitch together a single, coherent probabilistic forecast. Our experimental results demonstrate the benefits of our approach to produce extended-range probabilistic forecasts for regions and time periods of interest that are superior, in terms of skill, to individual NMME forecast models and commonly weighted models. The probabilistic forecast leverages the strengths of three NMME forecast models to predict environmental conditions for an area spanning from San Diego, CA to Honolulu, HI, seven months in-advance. Key findings include: weighted combinations of models are strictly better than individual models; machine-learned combinations are especially better; and forecasts produced using our approach have the highest rank probability skill score most often.

Improved water allocation utilizing probabilistic climate forecasts: Short-term water contracts in a risk management framework

NASA Astrophysics Data System (ADS)

Sankarasubramanian, A.; Lall, Upmanu; Souza Filho, Francisco Assis; Sharma, Ashish

2009-11-01

Probabilistic, seasonal to interannual streamflow forecasts are becoming increasingly available as the ability to model climate teleconnections is improving. However, water managers and practitioners have been slow to adopt such products, citing concerns with forecast skill. Essentially, a management risk is perceived in "gambling" with operations using a probabilistic forecast, while a system failure upon following existing operating policies is "protected" by the official rules or guidebook. In the presence of a prescribed system of prior allocation of releases under different storage or water availability conditions, the manager has little incentive to change. Innovation in allocation and operation is hence key to improved risk management using such forecasts. A participatory water allocation process that can effectively use probabilistic forecasts as part of an adaptive management strategy is introduced here. Users can express their demand for water through statements that cover the quantity needed at a particular reliability, the temporal distribution of the "allocation," the associated willingness to pay, and compensation in the event of contract nonperformance. The water manager then assesses feasible allocations using the probabilistic forecast that try to meet these criteria across all users. An iterative process between users and water manager could be used to formalize a set of short-term contracts that represent the resulting prioritized water allocation strategy over the operating period for which the forecast was issued. These contracts can be used to allocate water each year/season beyond long-term contracts that may have precedence. Thus, integrated supply and demand management can be achieved. In this paper, a single period multiuser optimization model that can support such an allocation process is presented. The application of this conceptual model is explored using data for the Jaguaribe Metropolitan Hydro System in Ceara, Brazil. The performance relative to the current allocation process is assessed in the context of whether such a model could support the proposed short-term contract based participatory process. A synthetic forecasting example is also used to explore the relative roles of forecast skill and reservoir storage in this framework.

A temporal-spatial postprocessing model for probabilistic run-off forecast. With a case study from Ulla-Førre with five catchments and ten lead times

NASA Astrophysics Data System (ADS)

Engeland, K.; Steinsland, I.

2012-04-01

This work is driven by the needs of next generation short term optimization methodology for hydro power production. Stochastic optimization are about to be introduced; i.e. optimizing when available resources (water) and utility (prices) are uncertain. In this paper we focus on the available resources, i.e. water, where uncertainty mainly comes from uncertainty in future runoff. When optimizing a water system all catchments and several lead times have to be considered simultaneously. Depending on the system of hydropower reservoirs, it might be a set of headwater catchments, a system of upstream /downstream reservoirs where water used from one catchment /dam arrives in a lower catchment maybe days later, or a combination of both. The aim of this paper is therefore to construct a simultaneous probabilistic forecast for several catchments and lead times, i.e. to provide a predictive distribution for the forecasts. Stochastic optimization methods need samples/ensembles of run-off forecasts as input. Hence, it should also be possible to sample from our probabilistic forecast. A post-processing approach is taken, and an error model based on Box- Cox transformation, power transform and a temporal-spatial copula model is used. It accounts for both between catchment and between lead time dependencies. In operational use it is strait forward to sample run-off ensembles from this models that inherits the catchment and lead time dependencies. The methodology is tested and demonstrated in the Ulla-Førre river system, and simultaneous probabilistic forecasts for five catchments and ten lead times are constructed. The methodology has enough flexibility to model operationally important features in this case study such as hetroscadasety, lead-time varying temporal dependency and lead-time varying inter-catchment dependency. Our model is evaluated using CRPS for marginal predictive distributions and energy score for joint predictive distribution. It is tested against deterministic run-off forecast, climatology forecast and a persistent forecast, and is found to be the better probabilistic forecast for lead time grater then two. From an operational point of view the results are interesting as the between catchment dependency gets stronger with longer lead-times.

Multi-parametric variational data assimilation for hydrological forecasting

NASA Astrophysics Data System (ADS)

Alvarado-Montero, R.; Schwanenberg, D.; Krahe, P.; Helmke, P.; Klein, B.

2017-12-01

Ensemble forecasting is increasingly applied in flow forecasting systems to provide users with a better understanding of forecast uncertainty and consequently to take better-informed decisions. A common practice in probabilistic streamflow forecasting is to force deterministic hydrological model with an ensemble of numerical weather predictions. This approach aims at the representation of meteorological uncertainty but neglects uncertainty of the hydrological model as well as its initial conditions. Complementary approaches use probabilistic data assimilation techniques to receive a variety of initial states or represent model uncertainty by model pools instead of single deterministic models. This paper introduces a novel approach that extends a variational data assimilation based on Moving Horizon Estimation to enable the assimilation of observations into multi-parametric model pools. It results in a probabilistic estimate of initial model states that takes into account the parametric model uncertainty in the data assimilation. The assimilation technique is applied to the uppermost area of River Main in Germany. We use different parametric pools, each of them with five parameter sets, to assimilate streamflow data, as well as remotely sensed data from the H-SAF project. We assess the impact of the assimilation in the lead time performance of perfect forecasts (i.e. observed data as forcing variables) as well as deterministic and probabilistic forecasts from ECMWF. The multi-parametric assimilation shows an improvement of up to 23% for CRPS performance and approximately 20% in Brier Skill Scores with respect to the deterministic approach. It also improves the skill of the forecast in terms of rank histogram and produces a narrower ensemble spread.

Hazard Monitoring of Growing Lava Flow Fields Using Seismic Tremor

NASA Astrophysics Data System (ADS)

Eibl, E. P. S.; Bean, C. J.; Jónsdottir, I.; Hoskuldsson, A.; Thordarson, T.; Coppola, D.; Witt, T.; Walter, T. R.

2017-12-01

An effusive eruption in 2014/15 created a 85 km2 large lava flow field in a remote location in the Icelandic highlands. The lava flows did not threaten any settlements or paved roads but they were nevertheless interdisciplinarily monitored in detail. Images from satellites and aircraft, ground based video monitoring, GPS and seismic recordings allowed the monitoring and reconstruction of a detailed time series of the growing lava flow field. While the use of satellite images and probabilistic modelling of lava flows are quite common tools to monitor the current and forecast the future growth direction, here we show that seismic recordings can be of use too. We installed a cluster of seismometers at 15 km from the vents and recorded the ground vibrations associated with the eruption. This seismic tremor was not only generated below the vents, but also at the edges of the growing lava flow field and indicated the parts of the lava flow field that were most actively growing. Whilst the time resolution is in the range of days for satellites, seismic stations easily sample continuously at 100 Hz and could therefore provide a much better resolution and estimate of the lava flow hazard in real-time.

An operational hydrological ensemble prediction system for the city of Zurich (Switzerland): assessing the added value of probabilistic forecasts

NASA Astrophysics Data System (ADS)

Addor, N.; Jaun, S.; Fundel, F.; Zappa, M.

2012-04-01

The Sihl River flows through Zurich, Switzerland's most populated city, for which it represents the largest flood threat. To anticipate extreme discharge events and provide decision support in case of flood risk, a hydrometeorological ensemble prediction system (HEPS) was launched operationally in 2008. This model chain relies on deterministic (COSMO-7) and probabilistic (COSMO-LEPS) atmospheric forecasts, which are used to force a semi-distributed hydrological model (PREVAH) coupled to a hydraulic model (FLORIS). The resulting hydrological forecasts are eventually communicated to the stakeholders involved in the Sihl discharge management. This fully operational setting provides a real framework with which we assessed the potential of deterministic and probabilistic discharge forecasts for flood mitigation. To study the suitability of HEPS for small-scale basins and to quantify the added value conveyed by the probability information, a 31-month reforecast was produced for the Sihl catchment (336 km2). Several metrics support the conclusion that the performance gain is of up to 2 days lead time for the catchment considered. Brier skill scores show that probabilistic hydrological forecasts outperform their deterministic counterparts for all the lead times and event intensities considered. The small size of the Sihl catchment does not prevent skillful discharge forecasts, but makes them particularly dependent on correct precipitation forecasts. Our evaluation stresses that the capacity of the model to provide confident and reliable mid-term probability forecasts for high discharges is limited. We finally highlight challenges for making decisions on the basis of hydrological predictions, and discuss the need for a tool to be used in addition to forecasts to compare the different mitigation actions possible in the Sihl catchment.

Visualisation and communication of probabilistic climate forecasts to renewable-energy policy makers

NASA Astrophysics Data System (ADS)

Steffen, Sophie; Lowe, Rachel; Davis, Melanie; Doblas-Reyes, Francisco J.; Rodó, Xavier

2014-05-01

Despite the strong dependence on weather and climate variability of the renewable-energy industry, and the existence of several initiatives towards demonstrating the added benefits of integrating probabilistic forecasts into energy decision-making processes, weather and climate forecasts are still under-utilised within the sector. Improved communication is fundamental to stimulate the use of climate forecast information within decision-making processes, in order to adapt to a highly climate dependent renewable-energy industry. This work focuses on improving the visualisation of climate forecast information, paying special attention to seasonal time scales. This activity is central to enhance climate services for renewable energy and to optimise the usefulness and usability of inherently complex climate information. In the realm of the Global Framework for Climate Services (GFCS) initiative, and subsequent European projects: Seasonal-to-Decadal Climate Prediction for the Improvement of European Climate Service (SPECS) and the European Provision of Regional Impacts Assessment in Seasonal and Decadal Timescales (EUPORIAS), this paper investigates the visualisation and communication of seasonal forecasts with regards to their usefulness and usability, to enable the development of a European climate service. The target end user is the group of renewable-energy policy makers, who are central to enhance climate services for the energy industry. The overall objective is to promote the wide-range dissemination and exchange of actionable climate information based on seasonal forecasts from Global Producing Centres (GPCs). It examines the existing main barriers and deficits. Examples of probabilistic climate forecasts from different GPC's are used to make a catalogue of current approaches, to assess their advantages and limitations and, finally, to recommend better alternatives. Interviews have been conducted with renewable-energy stakeholders to receive feedback for the improvement of existing visualisation techniques of forecasts. The overall aim is to establish a communication protocol for the visualisation of probabilistic climate forecasts, which does not currently exist. GPCs show their own probabilistic forecasts with limited consistency in their communication across different centres, which complicates the understanding for the end user. The recommended communication protocol for both the visualisation and description of climate forecasts can help to introduce a standard format and message to end users from several climate-sensitive sectors, such as energy, tourism, agriculture and health.

New Aspects of Probabilistic Forecast Verification Using Information Theory

NASA Astrophysics Data System (ADS)

Tödter, Julian; Ahrens, Bodo

2013-04-01

This work deals with information-theoretical methods in probabilistic forecast verification, particularly concerning ensemble forecasts. Recent findings concerning the "Ignorance Score" are shortly reviewed, then a consistent generalization to continuous forecasts is motivated. For ensemble-generated forecasts, the presented measures can be calculated exactly. The Brier Score (BS) and its generalizations to the multi-categorical Ranked Probability Score (RPS) and to the Continuous Ranked Probability Score (CRPS) are prominent verification measures for probabilistic forecasts. Particularly, their decompositions into measures quantifying the reliability, resolution and uncertainty of the forecasts are attractive. Information theory sets up a natural framework for forecast verification. Recently, it has been shown that the BS is a second-order approximation of the information-based Ignorance Score (IGN), which also contains easily interpretable components and can also be generalized to a ranked version (RIGN). Here, the IGN, its generalizations and decompositions are systematically discussed in analogy to the variants of the BS. Additionally, a Continuous Ranked IGN (CRIGN) is introduced in analogy to the CRPS. The useful properties of the conceptually appealing CRIGN are illustrated, together with an algorithm to evaluate its components reliability, resolution, and uncertainty for ensemble-generated forecasts. This algorithm can also be used to calculate the decomposition of the more traditional CRPS exactly. The applicability of the "new" measures is demonstrated in a small evaluation study of ensemble-based precipitation forecasts.

Communicating weather forecast uncertainty: Do individual differences matter?

PubMed

Grounds, Margaret A; Joslyn, Susan L

2018-03-01

Research suggests that people make better weather-related decisions when they are given numeric probabilities for critical outcomes (Joslyn & Leclerc, 2012, 2013). However, it is unclear whether all users can take advantage of probabilistic forecasts to the same extent. The research reported here assessed key cognitive and demographic factors to determine their relationship to the use of probabilistic forecasts to improve decision quality. In two studies, participants decided between spending resources to prevent icy conditions on roadways or risk a larger penalty when freezing temperatures occurred. Several forecast formats were tested, including a control condition with the night-time low temperature alone and experimental conditions that also included the probability of freezing and advice based on expected value. All but those with extremely low numeracy scores made better decisions with probabilistic forecasts. Importantly, no groups made worse decisions when probabilities were included. Moreover, numeracy was the best predictor of decision quality, regardless of forecast format, suggesting that the advantage may extend beyond understanding the forecast to general decision strategy issues. This research adds to a growing body of evidence that numerical uncertainty estimates may be an effective way to communicate weather danger to general public end users. (PsycINFO Database Record (c) 2018 APA, all rights reserved).

Anticipatory Competence and Ability to Probabilistic Forecasting in Adolescents: Research Results

ERIC Educational Resources Information Center

Akhmetzyanova, Anna I.

2016-01-01

The relevance of this problem is related to the urgent need to explain peculiarities of anticipation and probabilistic forecasting in adolescence. It has revealed a contradiction: on the one hand, the problem of anticipation in ontogenesis is well developed, and, on the other hand, there remain understudied mechanisms of anticipation in…

Probabilistic accounting of uncertainty in forecasts of species distributions under climate change

Treesearch

Seth J. Wenger; Nicholas A. Som; Daniel C. Dauwalter; Daniel J. Isaak; Helen M. Neville; Charles H. Luce; Jason B. Dunham; Michael K. Young; Kurt D. Fausch; Bruce E. Rieman

2013-01-01

Forecasts of species distributions under future climates are inherently uncertain, but there have been few attempts to describe this uncertainty comprehensively in a probabilistic manner. We developed a Monte Carlo approach that accounts for uncertainty within generalized linear regression models (parameter uncertainty and residual error), uncertainty among competing...

Forecasting volcanic unrest using seismicity: The good, the bad and the time consuming

NASA Astrophysics Data System (ADS)

Salvage, Rebecca; Neuberg, Jurgen W.

2013-04-01

Volcanic eruptions are inherently unpredictable in nature, with scientists struggling to forecast the type and timing of events, in particular in real time scenarios. Current understanding suggests that the use of statistical patterns within precursory datasets of seismicity prior to eruptive events could hold the potential to be used as real time forecasting tools. They allow us to determine times of clear deviation in data, which might be indicative of volcanic unrest. The identification of low frequency seismic swarms and the acceleration of this seismicity prior to observed volcanic unrest may be key in developing forecasting tools. The development of these real time forecasting models which can be implemented at volcano observatories is of particular importance since the identification of early warning signals allows danger to the proximal population to be minimized. We concentrate on understanding the significance and development of these seismic swarms as unrest develops at the volcano. In particular, analysis of accelerations in event rate, amplitude and energy rates released by seismicity prior to eruption suggests that these are important indicators of developing unrest. Real time analysis of these parameters simultaneously allows possible improvements to forecasting models. Although more time and computationally intense, cross correlation techniques applied to continuous seismicity prior to volcanic unrest scenarios allows all significant seismic events to be analysed, rather than only those which can be detected by an automated identification system. This may allow a more accurate forecast since all precursory seismicity can be taken into account. In addition, the classification of seismic events based on spectral characteristics may allow us to isolate individual types of signals which are responsible for certain types of unrest. In this way, we may be able to better forecast the type of eruption that may ensue, or at least some of its prevailing characteristics.

Is there a basis for preferring characteristic earthquakes over a Gutenberg–Richter distribution in probabilistic earthquake forecasting?

USGS Publications Warehouse

Parsons, Thomas E.; Geist, Eric L.

2009-01-01

The idea that faults rupture in repeated, characteristic earthquakes is central to most probabilistic earthquake forecasts. The concept is elegant in its simplicity, and if the same event has repeated itself multiple times in the past, we might anticipate the next. In practice however, assembling a fault-segmented characteristic earthquake rupture model can grow into a complex task laden with unquantified uncertainty. We weigh the evidence that supports characteristic earthquakes against a potentially simpler model made from extrapolation of a Gutenberg–Richter magnitude-frequency law to individual fault zones. We find that the Gutenberg–Richter model satisfies key data constraints used for earthquake forecasting equally well as a characteristic model. Therefore, judicious use of instrumental and historical earthquake catalogs enables large-earthquake-rate calculations with quantifiable uncertainty that should get at least equal weighting in probabilistic forecasting.

Forecasting the climate response to volcanic eruptions: prediction skill related to stratospheric aerosol forcing

NASA Astrophysics Data System (ADS)

Ménégoz, M.; Bilbao, R.; Bellprat, O.; Guemas, V.; Doblas-Reyes, F. J.

2018-06-01

The last major volcanic eruptions, the Agung in 1963, El Chichon in 1982 and Pinatubo in 1991, were each associated with a cooling of the troposphere that has been observed over large continental areas and over the Western Pacific, the Indian Ocean and the Southern Atlantic. Simultaneously, Eastern tropical Pacific temperatures increased due to prevailing El Niño conditions. Here we show that the pattern of these near-surface temperature anomalies is partly reproduced with decadal simulations of the EC-Earth model initialised with climate observations and forced with an estimate of the observed volcanic aerosol optical thickness. Sensitivity experiments highlight a cooling induced by the volcanic forcing, whereas El Niño events following the eruptions would have occurred even without volcanic eruptions. Focusing on the period 1961–2001, the main source of skill of this decadal forecast system during the first 2 years is related to the initialisation of the model. The contribution of the initialisation to the skill becomes smaller than the contribution of the volcanic forcing after two years, the latter being substantial in the Western Pacific, the Indian Ocean and the Western Atlantic. Two simple protocols for real time forecasts are investigated: using the forcing of a past volcanic eruption to simulate the forcing of a new one, and applying a two-year exponential decay to the initial stratospheric aerosol load observed at the beginning of the forecast. This second protocol applied in retrospective forecasts allows a partial reproduction of the skill attained with observed forcing.

Sensitivity analysis and uncertainty estimation in ash concentration simulations and tephra deposit daily forecasted at Mt. Etna, in Italy

NASA Astrophysics Data System (ADS)

Prestifilippo, Michele; Scollo, Simona; Tarantola, Stefano

2015-04-01

The uncertainty in volcanic ash forecasts may depend on our knowledge of the model input parameters and our capability to represent the dynamic of an incoming eruption. Forecasts help governments to reduce risks associated with volcanic eruptions and for this reason different kinds of analysis that help to understand the effect that each input parameter has on model outputs are necessary. We present an iterative approach based on the sequential combination of sensitivity analysis, parameter estimation procedure and Monte Carlo-based uncertainty analysis, applied to the lagrangian volcanic ash dispersal model PUFF. We modify the main input parameters as the total mass, the total grain-size distribution, the plume thickness, the shape of the eruption column, the sedimentation models and the diffusion coefficient, perform thousands of simulations and analyze the results. The study is carried out on two different Etna scenarios: the sub-plinian eruption of 22 July 1998 that formed an eruption column rising 12 km above sea level and lasted some minutes and the lava fountain eruption having features similar to the 2011-2013 events that produced eruption column high up to several kilometers above sea level and lasted some hours. Sensitivity analyses and uncertainty estimation results help us to address the measurements that volcanologists should perform during volcanic crisis to reduce the model uncertainty.

A global empirical system for probabilistic seasonal climate prediction

NASA Astrophysics Data System (ADS)

Eden, J. M.; van Oldenborgh, G. J.; Hawkins, E.; Suckling, E. B.

2015-12-01

Preparing for episodes with risks of anomalous weather a month to a year ahead is an important challenge for governments, non-governmental organisations, and private companies and is dependent on the availability of reliable forecasts. The majority of operational seasonal forecasts are made using process-based dynamical models, which are complex, computationally challenging and prone to biases. Empirical forecast approaches built on statistical models to represent physical processes offer an alternative to dynamical systems and can provide either a benchmark for comparison or independent supplementary forecasts. Here, we present a simple empirical system based on multiple linear regression for producing probabilistic forecasts of seasonal surface air temperature and precipitation across the globe. The global CO2-equivalent concentration is taken as the primary predictor; subsequent predictors, including large-scale modes of variability in the climate system and local-scale information, are selected on the basis of their physical relationship with the predictand. The focus given to the climate change signal as a source of skill and the probabilistic nature of the forecasts produced constitute a novel approach to global empirical prediction. Hindcasts for the period 1961-2013 are validated against observations using deterministic (correlation of seasonal means) and probabilistic (continuous rank probability skill scores) metrics. Good skill is found in many regions, particularly for surface air temperature and most notably in much of Europe during the spring and summer seasons. For precipitation, skill is generally limited to regions with known El Niño-Southern Oscillation (ENSO) teleconnections. The system is used in a quasi-operational framework to generate empirical seasonal forecasts on a monthly basis.

An empirical system for probabilistic seasonal climate prediction

NASA Astrophysics Data System (ADS)

Eden, Jonathan; van Oldenborgh, Geert Jan; Hawkins, Ed; Suckling, Emma

2016-04-01

Preparing for episodes with risks of anomalous weather a month to a year ahead is an important challenge for governments, non-governmental organisations, and private companies and is dependent on the availability of reliable forecasts. The majority of operational seasonal forecasts are made using process-based dynamical models, which are complex, computationally challenging and prone to biases. Empirical forecast approaches built on statistical models to represent physical processes offer an alternative to dynamical systems and can provide either a benchmark for comparison or independent supplementary forecasts. Here, we present a simple empirical system based on multiple linear regression for producing probabilistic forecasts of seasonal surface air temperature and precipitation across the globe. The global CO2-equivalent concentration is taken as the primary predictor; subsequent predictors, including large-scale modes of variability in the climate system and local-scale information, are selected on the basis of their physical relationship with the predictand. The focus given to the climate change signal as a source of skill and the probabilistic nature of the forecasts produced constitute a novel approach to global empirical prediction. Hindcasts for the period 1961-2013 are validated against observations using deterministic (correlation of seasonal means) and probabilistic (continuous rank probability skill scores) metrics. Good skill is found in many regions, particularly for surface air temperature and most notably in much of Europe during the spring and summer seasons. For precipitation, skill is generally limited to regions with known El Niño-Southern Oscillation (ENSO) teleconnections. The system is used in a quasi-operational framework to generate empirical seasonal forecasts on a monthly basis.

Decision Analysis Tools for Volcano Observatories

NASA Astrophysics Data System (ADS)

Hincks, T. H.; Aspinall, W.; Woo, G.

2005-12-01

Staff at volcano observatories are predominantly engaged in scientific activities related to volcano monitoring and instrumentation, data acquisition and analysis. Accordingly, the academic education and professional training of observatory staff tend to focus on these scientific functions. From time to time, however, staff may be called upon to provide decision support to government officials responsible for civil protection. Recognizing that Earth scientists may have limited technical familiarity with formal decision analysis methods, specialist software tools that assist decision support in a crisis should be welcome. A review is given of two software tools that have been under development recently. The first is for probabilistic risk assessment of human and economic loss from volcanic eruptions, and is of practical use in short and medium-term risk-informed planning of exclusion zones, post-disaster response, etc. A multiple branch event-tree architecture for the software, together with a formalism for ascribing probabilities to branches, have been developed within the context of the European Community EXPLORIS project. The second software tool utilizes the principles of the Bayesian Belief Network (BBN) for evidence-based assessment of volcanic state and probabilistic threat evaluation. This is of practical application in short-term volcano hazard forecasting and real-time crisis management, including the difficult challenge of deciding when an eruption is over. An open-source BBN library is the software foundation for this tool, which is capable of combining synoptically different strands of observational data from diverse monitoring sources. A conceptual vision is presented of the practical deployment of these decision analysis tools in a future volcano observatory environment. Summary retrospective analyses are given of previous volcanic crises to illustrate the hazard and risk insights gained from use of these tools.

The Hawaiian Volcano Observatory's current approach to forecasting lava flow hazards (Invited)

NASA Astrophysics Data System (ADS)

Kauahikaua, J. P.

2013-12-01

Hawaiian Volcanoes are best known for their frequent basaltic eruptions, which typically start with fast-moving channelized `a`a flows fed by high eruptions rates. If the flows continue, they generally transition into pahoehoe flows, fed by lower eruption rates, after a few days to weeks. Kilauea Volcano's ongoing eruption illustrates this--since 1986, effusion at Kilauea has mostly produced pahoehoe. The current state of lava flow simulation is quite advanced, but the simplicity of the models mean that they are most appropriately used during the first, most vigorous, days to weeks of an eruption - during the effusion of `a`a flows. Colleagues at INGV in Catania have shown decisively that MAGFLOW simulations utilizing satellite-derived eruption rates can be effective at estimating hazards during the initial periods of an eruption crisis. However, the algorithms do not simulate the complexity of pahoehoe flows. Forecasts of lava flow hazards are the most common form of volcanic hazard assessments made in Hawai`i. Communications with emergency managers over the last decade have relied on simple steepest-descent line maps, coupled with empirical lava flow advance rate information, to portray the imminence of lava flow hazard to nearby communities. Lavasheds, calculated as watersheds, are used as a broader context for the future flow paths and to advise on the utility of diversion efforts, should they be contemplated. The key is to communicate the uncertainty of any approach used to formulate a forecast and, if the forecast uses simple tools, these communications can be fairly straightforward. The calculation of steepest-descent paths and lavasheds relies on the accuracy of the digital elevation model (DEM) used, so the choice of DEM is critical. In Hawai`i, the best choice is not the most recent but is a 1980s-vintage 10-m DEM--more recent LIDAR and satellite radar DEM are referenced to the ellipsoid and include vegetation effects. On low-slope terrain, steepest descent lines calculated on a geoid-based DEM may differ significantly from those calculated on an ellipsoid-based DEM. Good estimates of lava flow advance rates can be obtained from empirical compilations of historical advance rates of Hawaiian lava flows. In this way, rates appropriate for observed flow types (`a`a or pahoehoe, channelized or not) can be applied. Eruption rate is arguably the most important factor, while slope is also significant for low eruption rates. Eruption rate, however, remains the most difficult parameter to estimate during an active eruption. The simplicity of the HVO approach is its major benefit. How much better can lava-flow advance be forecast for all types of lava flows? Will the improvements outweigh the increased uncertainty propagated through the simulation calculations? HVO continues to improve and evaluate its lava flow forecasting tools to provide better hazard assessments to emergency personnel.

Evaluating sub-seasonal skill in probabilistic forecasts of Atmospheric Rivers and associated extreme events

NASA Astrophysics Data System (ADS)

Subramanian, A. C.; Lavers, D.; Matsueda, M.; Shukla, S.; Cayan, D. R.; Ralph, M.

2017-12-01

Atmospheric rivers (ARs) - elongated plumes of intense moisture transport - are a primary source of hydrological extremes, water resources and impactful weather along the West Coast of North America and Europe. There is strong demand in the water management, societal infrastructure and humanitarian sectors for reliable sub-seasonal forecasts, particularly of extreme events, such as floods and droughts so that actions to mitigate disastrous impacts can be taken with sufficient lead-time. Many recent studies have shown that ARs in the Pacific and the Atlantic are modulated by large-scale modes of climate variability. Leveraging the improved understanding of how these large-scale climate modes modulate the ARs in these two basins, we use the state-of-the-art multi-model forecast systems such as the North American Multi-Model Ensemble (NMME) and the Subseasonal-to-Seasonal (S2S) database to help inform and assess the probabilistic prediction of ARs and related extreme weather events over the North American and European West Coasts. We will present results from evaluating probabilistic forecasts of extreme precipitation and AR activity at the sub-seasonal scale. In particular, results from the comparison of two winters (2015-16 and 2016-17) will be shown, winters which defied canonical El Niño teleconnection patterns over North America and Europe. We further extend this study to analyze probabilistic forecast skill of AR events in these two basins and the variability in forecast skill during certain regimes of large-scale climate modes.

The Ongoing 2011 Eruption of Cordón Caulle (Southern Andes) and its Related Hazards

NASA Astrophysics Data System (ADS)

Amigo, A.; Lara, L. E.; Silva, C.; Orozco, G.; Bertin, D.

2011-12-01

On June 4, 2011, at 18:45 UTC, Cordón Caulle volcano (Southern Andes, 40.52S, 72.14W) erupted explosively after 51 years of quiescence. The last eruption occurred in 1960 and was triggered by the great Mw 9.5 Chile earthquake. The ongoing eruption started after 2 months of increased shallow seismicity as recorded by OVDAS (the volcano observatory at Sernageomin). This close monitoring effort allowed a timely eruption forecast with at least 3 hours of warning, which facilitated the crisis response. In addition to this successful performance, for the first time in Chile volcanic hazards were assessed in advance supporting the emergency management. In particular, tephra dispersal was daily forecasted using the ASHFALL advection-diffusion model and potential lahars and PDC impact zones were delineated according to numerical approaches. The first eruptive stage lasted 27 hours. It was characterized by ca. 15-km strong Plinian-like column, associated with the emission of 0.2 - 0.4 km3 of magma (DRE). Tephra fallout mostly occurred in Chile and Argentina, although fine particles and aerosols circumnavigated the globe twice, causing disruptions on air navigation across the Southern Hemisphere. The second ongoing eruptive stage has been characterized by persistent weak plumes and lava emission at effusion rates in the range of 20 and 60 m3/s, which total volume is estimated <0.20 km3 (at the end of July 2011). Eruptive products have virtually the same bulk composition as those of the historical 1921 and 1960 eruptions, corresponding to phenocryst-poor rhyodacites (67 - 70% SiO2) for what a pre-eruptive temperature of ca. 920C could be inferred. In contrast to the previous eruptive cycles, the ongoing eruption has not evolved (at the time of writing) as a fissure eruption although the vent is atop of fault scarp that borders the Pleistocene-Holocene extensional graben of the Cordón Caulle. This episode is a good case of successful eruption forecast and hazards assessment but it is also an important case-study of silicic eruptions in an arc segment where mostly mafic magmas have been erupted during the Holocene.

Real-time eruption forecasting using the material Failure Forecast Method with a Bayesian approach

NASA Astrophysics Data System (ADS)

Boué, A.; Lesage, P.; Cortés, G.; Valette, B.; Reyes-Dávila, G.

2015-04-01

Many attempts for deterministic forecasting of eruptions and landslides have been performed using the material Failure Forecast Method (FFM). This method consists in adjusting an empirical power law on precursory patterns of seismicity or deformation. Until now, most of the studies have presented hindsight forecasts based on complete time series of precursors and do not evaluate the ability of the method for carrying out real-time forecasting with partial precursory sequences. In this study, we present a rigorous approach of the FFM designed for real-time applications on volcano-seismic precursors. We use a Bayesian approach based on the FFM theory and an automatic classification of seismic events. The probability distributions of the data deduced from the performance of this classification are used as input. As output, it provides the probability of the forecast time at each observation time before the eruption. The spread of the a posteriori probability density function of the prediction time and its stability with respect to the observation time are used as criteria to evaluate the reliability of the forecast. We test the method on precursory accelerations of long-period seismicity prior to vulcanian explosions at Volcán de Colima (Mexico). For explosions preceded by a single phase of seismic acceleration, we obtain accurate and reliable forecasts using approximately 80% of the whole precursory sequence. It is, however, more difficult to apply the method to multiple acceleration patterns.

Improving medium-range ensemble streamflow forecasts through statistical post-processing

NASA Astrophysics Data System (ADS)

Mendoza, Pablo; Wood, Andy; Clark, Elizabeth; Nijssen, Bart; Clark, Martyn; Ramos, Maria-Helena; Nowak, Kenneth; Arnold, Jeffrey

2017-04-01

Probabilistic hydrologic forecasts are a powerful source of information for decision-making in water resources operations. A common approach is the hydrologic model-based generation of streamflow forecast ensembles, which can be implemented to account for different sources of uncertainties - e.g., from initial hydrologic conditions (IHCs), weather forecasts, and hydrologic model structure and parameters. In practice, hydrologic ensemble forecasts typically have biases and spread errors stemming from errors in the aforementioned elements, resulting in a degradation of probabilistic properties. In this work, we compare several statistical post-processing techniques applied to medium-range ensemble streamflow forecasts obtained with the System for Hydromet Applications, Research and Prediction (SHARP). SHARP is a fully automated prediction system for the assessment and demonstration of short-term to seasonal streamflow forecasting applications, developed by the National Center for Atmospheric Research, University of Washington, U.S. Army Corps of Engineers, and U.S. Bureau of Reclamation. The suite of post-processing techniques includes linear blending, quantile mapping, extended logistic regression, quantile regression, ensemble analogs, and the generalized linear model post-processor (GLMPP). We assess and compare these techniques using multi-year hindcasts in several river basins in the western US. This presentation discusses preliminary findings about the effectiveness of the techniques for improving probabilistic skill, reliability, discrimination, sharpness and resolution.

The visualisation and communication of probabilistic climate forecasts to renewable energy policy makers

NASA Astrophysics Data System (ADS)

Doblas-Reyes, F.; Steffen, S.; Lowe, R.; Davis, M.; Rodó, X.

2013-12-01

Despite the strong dependence of weather and climate variability on the renewable energy industry, and several initiatives towards demonstrating the added benefits of integrating probabilistic forecasts into energy decision making process, they are still under-utilised within the sector. Improved communication is fundamental to stimulate the use of climate forecast information within decision-making processes, in order to adapt to a highly climate dependent renewable energy industry. This paper focuses on improving the visualisation of climate forecast information, paying special attention to seasonal to decadal (s2d) timescales. This is central to enhance climate services for renewable energy, and optimise the usefulness and usability of inherently complex climate information. In the realm of the Global Framework for Climate Services (GFCS) initiative, and subsequent European projects: Seasonal-to-Decadal Climate Prediction for the Improvement of European Climate Service (SPECS) and the European Provision of Regional Impacts Assessment in Seasonal and Decadal Timescales (EUPORIAS), this paper investigates the visualisation and communication of s2d forecasts with regards to their usefulness and usability, to enable the development of a European climate service. The target end user will be renewable energy policy makers, who are central to enhance climate services for the energy industry. The overall objective is to promote the wide-range dissemination and exchange of actionable climate information based on s2d forecasts from Global Producing Centres (GPC's). Therefore, it is crucial to examine the existing main barriers and deficits. Examples of probabilistic climate forecasts from different GPC's were used to prepare a catalogue of current approaches, to assess their advantages and limitations and finally to recommend better alternatives. In parallel, interviews were conducted with renewable energy stakeholders to receive feedback for the improvement of existing visualisation techniques of forecasts. The overall aim is to establish a communication protocol for the visualisation of probabilistic climate forecasts, which does not currently exist. Global Producing Centres show their own probabilistic forecasts with limited consistency in their communication across different centres, which complicates the understanding for the end user. A communication protocol for both the visualisation and description of climate forecasts can help to introduce a standard format and message to end users from several climate-sensitive sectors, such as energy, tourism, agriculture and health. It is hoped that this work will facilitate the improvement of decision-making processes relying on forecast information and enable their wide-range dissemination based on a standardised approach.

Added value of non-calibrated and BMA calibrated AEMET-SREPS probabilistic forecasts: the 24 January 2009 extreme wind event over Catalonia

NASA Astrophysics Data System (ADS)

Escriba, P. A.; Callado, A.; Santos, D.; Santos, C.; Simarro, J.; García-Moya, J. A.

2009-09-01

At 00 UTC 24 January 2009 an explosive ciclogenesis originated over the Atlantic Ocean reached its maximum intensity with observed surface pressures lower than 970 hPa on its center and placed at Gulf of Vizcaya. During its path through southern France this low caused strong westerly and north-westerly winds over the Iberian Peninsula higher than 150 km/h at some places. These extreme winds leaved 10 casualties in Spain, 8 of them in Catalonia. The aim of this work is to show whether exists an added value in the short range prediction of the 24 January 2009 strong winds when using the Short Range Ensemble Prediction System (SREPS) of the Spanish Meteorological Agency (AEMET), with respect to the operational forecasting tools. This study emphasizes two aspects of probabilistic forecasting: the ability of a 3-day forecast of warn an extreme windy event and the ability of quantifying the predictability of the event so that giving value to deterministic forecast. Two type of probabilistic forecasts of wind are carried out, a non-calibrated and a calibrated one using Bayesian Model Averaging (BMA). AEMET runs daily experimentally SREPS twice a day (00 and 12 UTC). This system consists of 20 members that are constructed by integrating 5 local area models, COSMO (COSMO), HIRLAM (HIRLAM Consortium), HRM (DWD), MM5 (NOAA) and UM (UKMO), at 25 km of horizontal resolution. Each model uses 4 different initial and boundary conditions, the global models GFS (NCEP), GME (DWD), IFS (ECMWF) and UM. By this way it is obtained a probabilistic forecast that takes into account the initial, the contour and the model errors. BMA is a statistical tool for combining predictive probability functions from different sources. The BMA predictive probability density function (PDF) is a weighted average of PDFs centered on the individual bias-corrected forecasts. The weights are equal to posterior probabilities of the models generating the forecasts and reflect the skill of the ensemble members. Here BMA is applied to provide probabilistic forecasts of wind speed. In this work several forecasts for different time ranges (H+72, H+48 and H+24) of 10 meters wind speed over Catalonia are verified subjectively at one of the instants of maximum intensity, 12 UTC 24 January 2009. On one hand, three probabilistic forecasts are compared, ECMWF EPS, non-calibrated SREPS and calibrated SREPS. On the other hand, the relationship between predictability and skill of deterministic forecast is studied by looking at HIRLAM 0.16 deterministic forecasts of the event. Verification is focused on location and intensity of 10 meters wind speed and 10-minutal measures from AEMET automatic ground stations are used as observations. The results indicate that SREPS is able to forecast three days ahead mean winds higher than 36 km/h and that correctly localizes them with a significant probability of ocurrence in the affected area. The probability is higher after BMA calibration of the ensemble. The fact that probability of strong winds is high allows us to state that the predictability of the event is also high and, as a consequence, deterministic forecasts are more reliable. This is confirmed when verifying HIRLAM deterministic forecasts against observed values.

Bayesian Probabilistic Projections of Life Expectancy for All Countries

PubMed Central

Raftery, Adrian E.; Chunn, Jennifer L.; Gerland, Patrick; Ševčíková, Hana

2014-01-01

We propose a Bayesian hierarchical model for producing probabilistic forecasts of male period life expectancy at birth for all the countries of the world from the present to 2100. Such forecasts would be an input to the production of probabilistic population projections for all countries, which is currently being considered by the United Nations. To evaluate the method, we did an out-of-sample cross-validation experiment, fitting the model to the data from 1950–1995, and using the estimated model to forecast for the subsequent ten years. The ten-year predictions had a mean absolute error of about 1 year, about 40% less than the current UN methodology. The probabilistic forecasts were calibrated, in the sense that (for example) the 80% prediction intervals contained the truth about 80% of the time. We illustrate our method with results from Madagascar (a typical country with steadily improving life expectancy), Latvia (a country that has had a mortality crisis), and Japan (a leading country). We also show aggregated results for South Asia, a region with eight countries. Free publicly available R software packages called bayesLife and bayesDem are available to implement the method. PMID:23494599

On the skill of various ensemble spread estimators for probabilistic short range wind forecasting

NASA Astrophysics Data System (ADS)

Kann, A.

2012-05-01

A variety of applications ranging from civil protection associated with severe weather to economical interests are heavily dependent on meteorological information. For example, a precise planning of the energy supply with a high share of renewables requires detailed meteorological information on high temporal and spatial resolution. With respect to wind power, detailed analyses and forecasts of wind speed are of crucial interest for the energy management. Although the applicability and the current skill of state-of-the-art probabilistic short range forecasts has increased during the last years, ensemble systems still show systematic deficiencies which limit its practical use. This paper presents methods to improve the ensemble skill of 10-m wind speed forecasts by combining deterministic information from a nowcasting system on very high horizontal resolution with uncertainty estimates from a limited area ensemble system. It is shown for a one month validation period that a statistical post-processing procedure (a modified non-homogeneous Gaussian regression) adds further skill to the probabilistic forecasts, especially beyond the nowcasting range after +6 h.

Willingness-to-pay for a probabilistic flood forecast: a risk-based decision-making game

NASA Astrophysics Data System (ADS)

Arnal, Louise; Ramos, Maria-Helena; Coughlan, Erin; Cloke, Hannah L.; Stephens, Elisabeth; Wetterhall, Fredrik; van Andel, Schalk-Jan; Pappenberger, Florian

2016-04-01

Forecast uncertainty is a twofold issue, as it constitutes both an added value and a challenge for the forecaster and the user of the forecasts. Many authors have demonstrated the added (economic) value of probabilistic forecasts over deterministic forecasts for a diversity of activities in the water sector (e.g. flood protection, hydroelectric power management and navigation). However, the richness of the information is also a source of challenges for operational uses, due partially to the difficulty to transform the probability of occurrence of an event into a binary decision. The setup and the results of a risk-based decision-making experiment, designed as a game on the topic of flood protection mitigation, called ``How much are you prepared to pay for a forecast?'', will be presented. The game was played at several workshops in 2015, including during this session at the EGU conference in 2015, and a total of 129 worksheets were collected and analysed. The aim of this experiment was to contribute to the understanding of the role of probabilistic forecasts in decision-making processes and their perceived value by decision-makers. Based on the participants' willingness-to-pay for a forecast, the results of the game showed that the value (or the usefulness) of a forecast depends on several factors, including the way users perceive the quality of their forecasts and link it to the perception of their own performances as decision-makers. Balancing avoided costs and the cost (or the benefit) of having forecasts available for making decisions is not straightforward, even in a simplified game situation, and is a topic that deserves more attention from the hydrological forecasting community in the future.

The Forecast Interpretation Tool—a Monte Carlo technique for blending climatic distributions with probabilistic forecasts

USGS Publications Warehouse

Husak, Gregory J.; Michaelsen, Joel; Kyriakidis, P.; Verdin, James P.; Funk, Chris; Galu, Gideon

2011-01-01

Probabilistic forecasts are produced from a variety of outlets to help predict rainfall, and other meteorological events, for periods of 1 month or more. Such forecasts are expressed as probabilities of a rainfall event, e.g. being in the upper, middle, or lower third of the relevant distribution of rainfall in the region. The impact of these forecasts on the expectation for the event is not always clear or easily conveyed. This article proposes a technique based on Monte Carlo simulation for adjusting existing climatologic statistical parameters to match forecast information, resulting in new parameters defining the probability of events for the forecast interval. The resulting parameters are shown to approximate the forecasts with reasonable accuracy. To show the value of the technique as an application for seasonal rainfall, it is used with consensus forecast developed for the Greater Horn of Africa for the 2009 March-April-May season. An alternative, analytical approach is also proposed, and discussed in comparison to the first simulation-based technique.

Benefits of volcano monitoring far outweigh costs - the case of Mount Pinatubo

USGS Publications Warehouse

Newhall, Chris G.; Hendley, James W.; Stauffer, Peter H.

1997-01-01

The climactic June 1991 eruption of Mount Pinatubo, Philippines, was the largest volcanic eruption in this century to affect a heavily populated area. Because it was forecast by scientists from the Philippine Institute of Volcanology and Seismology and the U.S. Geological Survey, civil and military leaders were able to order massive evacuations and take measures to protect property before the eruption. Thousands of lives were saved and hundreds of millions of dollars in property losses averted. The savings in property alone were many times the total costs of the forecasting and evacuations.

Evaluation of the Plant-Craig stochastic convection scheme in an ensemble forecasting system

NASA Astrophysics Data System (ADS)

Keane, R. J.; Plant, R. S.; Tennant, W. J.

2015-12-01

The Plant-Craig stochastic convection parameterization (version 2.0) is implemented in the Met Office Regional Ensemble Prediction System (MOGREPS-R) and is assessed in comparison with the standard convection scheme with a simple stochastic element only, from random parameter variation. A set of 34 ensemble forecasts, each with 24 members, is considered, over the month of July 2009. Deterministic and probabilistic measures of the precipitation forecasts are assessed. The Plant-Craig parameterization is found to improve probabilistic forecast measures, particularly the results for lower precipitation thresholds. The impact on deterministic forecasts at the grid scale is neutral, although the Plant-Craig scheme does deliver improvements when forecasts are made over larger areas. The improvements found are greater in conditions of relatively weak synoptic forcing, for which convective precipitation is likely to be less predictable.

Forecasts and predictions of eruptive activity at Mount St. Helens, USA: 1975-1984

USGS Publications Warehouse

Swanson, D.A.; Casadevall, T.J.; Dzurisin, D.; Holcomb, R.T.; Newhall, C.G.; Malone, S.D.; Weaver, C.S.

1985-01-01

Public statements about volcanic activity at Mount St. Helens include factual statements, forecasts, and predictions. A factual statement describes current conditions but does not anticipate future events. A forecast is a comparatively imprecise statement of the time, place, and nature of expected activity. A prediction is a comparatively precise statement of the time, place, and ideally, the nature and size of impending activity. A prediction usually covers a shorter time period than a forecast and is generally based dominantly on interpretations and measurements of ongoing processes and secondarily on a projection of past history. The three types of statements grade from one to another, and distinctions are sometimes arbitrary. Forecasts and predictions at Mount St. Helens became increasingly precise from 1975 to 1982. Stratigraphic studies led to a long-range forecast in 1975 of renewed eruptive activity at Mount St. Helens, possibly before the end of the century. On the basis of seismic, geodetic and geologic data, general forecasts for a landslide and eruption were issued in April 1980, before the catastrophic blast and landslide on 18 May 1980. All extrusions except two from June 1980 to the end of 1984 were predicted on the basis of integrated geophysical, geochemical, and geologic monitoring. The two extrusions that were not predicted were preceded by explosions that removed a substantial part of the dome, reducing confining pressure and essentially short-circuiting the normal precursors. ?? 1985.

The VUELCO project consortium: new interdisciplinary research for improved risk mitigation and management during volcanic unrest

NASA Astrophysics Data System (ADS)

Gottsmann, J.

2012-04-01

Volcanic unrest is a complex multi-hazard phenomenon of volcanism. The fact that unrest may, but not necessarily must lead to an imminent eruption contributes significant uncertainty to short-term hazard assessment of volcanic activity world-wide. Although it is reasonable to assume that all eruptions are associated with precursory activity of some sort, the knowledge of the causative links between subsurface processes, resulting unrest signals and imminent eruption is, today, inadequate to deal effectively with crises of volcanic unrest. This results predominantly from the uncertainties in identifying the causative processes of unrest and as a consequence in forecasting its short-term evolution. However, key for effective risk mitigation and management during unrest is the early and reliable identification of changes in the subsurface dynamics of a volcano and their assessment as precursors to an impending eruption. The VUELCO project consortium has come together for a multi-disciplinary attack on the origin, nature and significance of volcanic unrest from the scientific contributions generated by collaboration of ten partners in Europe and Latin America. Dissecting the science of monitoring data from unrest periods at six type volcanoes in Italy, Spain, the West Indies, Mexico and Ecuador the consortium will create global strategies for 1) enhanced monitoring capacity and value, 2) mechanistic data interpretation and 3) identification of reliable eruption precursors; all from the geophysical, geochemical and geodetic fingerprints of unrest episodes. Experiments will establish a mechanistic understanding of subsurface processes capable of inducing unrest and aid in identifying key volcano monitoring parameters indicative of the nature of unrest processes. Numerical models will help establish a link between the processes and volcano monitoring data to inform on the causes of unrest and its short-term evolution. Using uncertainty assessment and new short-term probabilistic hazard forecasting tools the scientific knowledge base will provide the crucial parameters for a comprehensive and best-practice approach to 1) risk mitigation, 2) communication, 3) decision-making and 4) crisis management during unrest periods. The VUELCO project consortium efforts will generate guidance in the definition and implementation of strategic options for effective risk mitigation, management and governance during unrest episodes. Such a mechanistic platform of understanding, impacting on the synergy of scientists, policy-makers, civil protection authorities, decision-makers, and the public, will place volcanic unrest management on a new basis, with European expertise at its peak. The project is financed by the European Commission under the 7th Framework Programme for Research and Technological Development, Area "Environment".

Improving Global Forecast System of extreme precipitation events with regional statistical model: Application of quantile-based probabilistic forecasts

NASA Astrophysics Data System (ADS)

Shastri, Hiteshri; Ghosh, Subimal; Karmakar, Subhankar

2017-02-01

Forecasting of extreme precipitation events at a regional scale is of high importance due to their severe impacts on society. The impacts are stronger in urban regions due to high flood potential as well high population density leading to high vulnerability. Although significant scientific improvements took place in the global models for weather forecasting, they are still not adequate at a regional scale (e.g., for an urban region) with high false alarms and low detection. There has been a need to improve the weather forecast skill at a local scale with probabilistic outcome. Here we develop a methodology with quantile regression, where the reliably simulated variables from Global Forecast System are used as predictors and different quantiles of rainfall are generated corresponding to that set of predictors. We apply this method to a flood-prone coastal city of India, Mumbai, which has experienced severe floods in recent years. We find significant improvements in the forecast with high detection and skill scores. We apply the methodology to 10 ensemble members of Global Ensemble Forecast System and find a reduction in ensemble uncertainty of precipitation across realizations with respect to that of original precipitation forecasts. We validate our model for the monsoon season of 2006 and 2007, which are independent of the training/calibration data set used in the study. We find promising results and emphasize to implement such data-driven methods for a better probabilistic forecast at an urban scale primarily for an early flood warning.

IEA Wind Task 36 Forecasting

NASA Astrophysics Data System (ADS)

Giebel, Gregor; Cline, Joel; Frank, Helmut; Shaw, Will; Pinson, Pierre; Hodge, Bri-Mathias; Kariniotakis, Georges; Sempreviva, Anna Maria; Draxl, Caroline

2017-04-01

Wind power forecasts have been used operatively for over 20 years. Despite this fact, there are still several possibilities to improve the forecasts, both from the weather prediction side and from the usage of the forecasts. The new International Energy Agency (IEA) Task on Wind Power Forecasting tries to organise international collaboration, among national weather centres with an interest and/or large projects on wind forecast improvements (NOAA, DWD, UK MetOffice, …) and operational forecaster and forecast users. The Task is divided in three work packages: Firstly, a collaboration on the improvement of the scientific basis for the wind predictions themselves. This includes numerical weather prediction model physics, but also widely distributed information on accessible datasets for verification. Secondly, we will be aiming at an international pre-standard (an IEA Recommended Practice) on benchmarking and comparing wind power forecasts, including probabilistic forecasts aiming at industry and forecasters alike. This WP will also organise benchmarks, in cooperation with the IEA Task WakeBench. Thirdly, we will be engaging end users aiming at dissemination of the best practice in the usage of wind power predictions, especially probabilistic ones. The Operating Agent is Gregor Giebel of DTU, Co-Operating Agent is Joel Cline of the US Department of Energy. Collaboration in the task is solicited from everyone interested in the forecasting business. We will collaborate with IEA Task 31 Wakebench, which developed the Windbench benchmarking platform, which this task will use for forecasting benchmarks. The task runs for three years, 2016-2018. Main deliverables are an up-to-date list of current projects and main project results, including datasets which can be used by researchers around the world to improve their own models, an IEA Recommended Practice on performance evaluation of probabilistic forecasts, a position paper regarding the use of probabilistic forecasts, and one or more benchmark studies implemented on the Windbench platform hosted at CENER. Additionally, spreading of relevant information in both the forecasters and the users community is paramount. The poster also shows the work done in the first half of the Task, e.g. the collection of available datasets and the learnings from a public workshop on 9 June in Barcelona on Experiences with the Use of Forecasts and Gaps in Research. Participation is open for all interested parties in member states of the IEA Annex on Wind Power, see ieawind.org for the up-to-date list. For collaboration, please contact the author grgi@dtu.dk).

Probabilistic evaluation of the physical impact of future tephra fallout events for the Island of Vulcano, Italy

NASA Astrophysics Data System (ADS)

Biass, Sebastien; Bonadonna, Costanza; di Traglia, Federico; Pistolesi, Marco; Rosi, Mauro; Lestuzzi, Pierino

2016-05-01

A first probabilistic scenario-based hazard assessment for tephra fallout is presented for La Fossa volcano (Vulcano Island, Italy) and subsequently used to assess the impact on the built environment. Eruption scenarios are based upon the stratigraphy produced by the last 1000 years of activity at Vulcano and include long-lasting Vulcanian and sub-Plinian eruptions. A new method is proposed to quantify the evolution through time of the hazard associated with pulsatory Vulcanian eruptions lasting from weeks to years, and the increase in hazard related to typical rainfall events around Sicily is also accounted for. The impact assessment on the roofs is performed by combining a field characterization of the buildings with the composite European vulnerability curves for typical roofing stocks. Results show that a sub-Plinian eruption of VEI 2 is not likely to affect buildings, whereas a sub-Plinian eruption of VEI 3 results in 90 % of the building stock having a ≥12 % probability of collapse. The hazard related to long-lasting Vulcanian eruptions evolves through time, and our analysis shows that the town of Il Piano, located downwind of the preferential wind patterns, is likely to reach critical tephra accumulations for roof collapse 5-9 months after the onset of the eruption. If no cleaning measures are taken, half of the building stock has a probability >20 % of suffering roof collapse.

Analog-Based Postprocessing of Navigation-Related Hydrological Ensemble Forecasts

NASA Astrophysics Data System (ADS)

Hemri, S.; Klein, B.

2017-11-01

Inland waterway transport benefits from probabilistic forecasts of water levels as they allow to optimize the ship load and, hence, to minimize the transport costs. Probabilistic state-of-the-art hydrologic ensemble forecasts inherit biases and dispersion errors from the atmospheric ensemble forecasts they are driven with. The use of statistical postprocessing techniques like ensemble model output statistics (EMOS) allows for a reduction of these systematic errors by fitting a statistical model based on training data. In this study, training periods for EMOS are selected based on forecast analogs, i.e., historical forecasts that are similar to the forecast to be verified. Due to the strong autocorrelation of water levels, forecast analogs have to be selected based on entire forecast hydrographs in order to guarantee similar hydrograph shapes. Custom-tailored measures of similarity for forecast hydrographs comprise hydrological series distance (SD), the hydrological matching algorithm (HMA), and dynamic time warping (DTW). Verification against observations reveals that EMOS forecasts for water level at three gauges along the river Rhine with training periods selected based on SD, HMA, and DTW compare favorably with reference EMOS forecasts, which are based on either seasonal training periods or on training periods obtained by dividing the hydrological forecast trajectories into runoff regimes.

A probabilistic verification score for contours demonstrated with idealized ice-edge forecasts

NASA Astrophysics Data System (ADS)

Goessling, Helge; Jung, Thomas

2017-04-01

We introduce a probabilistic verification score for ensemble-based forecasts of contours: the Spatial Probability Score (SPS). Defined as the spatial integral of local (Half) Brier Scores, the SPS can be considered the spatial analog of the Continuous Ranked Probability Score (CRPS). Applying the SPS to idealized seasonal ensemble forecasts of the Arctic sea-ice edge in a global coupled climate model, we demonstrate that the SPS responds properly to ensemble size, bias, and spread. When applied to individual forecasts or ensemble means (or quantiles), the SPS is reduced to the 'volume' of mismatch, in case of the ice edge corresponding to the Integrated Ice Edge Error (IIEE).

Probabilistic Volcanic Multi-Hazard Assessment at Somma-Vesuvius (Italy): coupling Bayesian Belief Networks with a physical model for lahar propagation

NASA Astrophysics Data System (ADS)

Tierz, Pablo; Woodhouse, Mark; Phillips, Jeremy; Sandri, Laura; Selva, Jacopo; Marzocchi, Warner; Odbert, Henry

2017-04-01

Volcanoes are extremely complex physico-chemical systems where magma formed at depth breaks into the planet's surface resulting in major hazards from local to global scales. Volcano physics are dominated by non-linearities, and complicated spatio-temporal interrelationships which make volcanic hazards stochastic (i.e. not deterministic) by nature. In this context, probabilistic assessments are required to quantify the large uncertainties related to volcanic hazards. Moreover, volcanoes are typically multi-hazard environments where different hazardous processes can occur whether simultaneously or in succession. In particular, explosive volcanoes are able to accumulate, through tephra fallout and Pyroclastic Density Currents (PDCs), large amounts of pyroclastic material into the drainage basins surrounding the volcano. This addition of fresh particulate material alters the local/regional hydrogeological equilibrium and increases the frequency and magnitude of sediment-rich aqueous flows, commonly known as lahars. The initiation and volume of rain-triggered lahars may depend on: rainfall intensity and duration; antecedent rainfall; terrain slope; thickness, permeability and hydraulic diffusivity of the tephra deposit; etc. Quantifying these complex interrelationships (and their uncertainties), in a tractable manner, requires a structured but flexible probabilistic approach. A Bayesian Belief Network (BBN) is a directed acyclic graph that allows the representation of the joint probability distribution for a set of uncertain variables in a compact and efficient way, by exploiting unconditional and conditional independences between these variables. Once constructed and parametrized, the BBN uses Bayesian inference to perform causal (e.g. forecast) and/or evidential reasoning (e.g. explanation) about query variables, given some evidence. In this work, we illustrate how BBNs can be used to model the influence of several variables on the generation of rain-triggered lahars and, finally, assess the probability of occurrence of lahars of different volumes. The information utilized to parametrize the BBNs includes: (1) datasets of lahar observations; (2) numerical modelling of tephra fallout and PDCs; and (3) literature data. The BBN framework provides an opportunity to quantitatively combine these different types of evidence and use them to derive a rational approach to lahar forecasting. Lastly, we couple the BBN assessments with a shallow-water physical model for lahar propagation in order to attach probabilities to the simulated hazard footprints. We develop our methodology at Somma-Vesuvius (Italy), an explosive volcano prone to rain-triggered lahars or debris flows whether right after an eruption or during inter-eruptive periods. Accounting for the variability in tephra-fallout and dense-PDC propagation and the main geomorphological features of the catchments around Somma-Vesuvius, the areas most likely of forming medium-large lahars are the flanks of the volcano and the Sarno mountains towards the east.

Exploring the influence of vent location and eruption style on tephra fall hazard from the Okataina Volcanic Centre, New Zealand

NASA Astrophysics Data System (ADS)

Thompson, Mary Anne; Lindsay, Jan M.; Sandri, Laura; Biass, Sébastien; Bonadonna, Costanza; Jolly, Gill; Marzocchi, Warner

2015-05-01

Uncertainties in modelling volcanic hazards are often amplified in geographically large systems which have a diverse eruption history that comprises variable eruption styles from many different vent locations. The ~700 km2 Okataina Volcanic Centre (OVC) is a caldera complex in New Zealand which has displayed a range of eruption styles and compositions over its current phase of activity (26 ka-present), including one basaltic maar-forming eruption, one basaltic Plinian eruption and nine rhyolitic Plinian eruptions. All three of these eruption styles occurred within the past 3.5 ky, and any of these styles could occur in the event of a future eruption. The location of a future eruption is also unknown. Future vents could potentially open in one of three different areas which have been activated in the past 26 ky at the OVC: the Tarawera linear vent zone (LVZ) (five eruptions), the Haroharo LVZ (five eruptions) or outside of these LVZs (one eruption). A future rhyolitic or basaltic Plinian eruption from the OVC is likely to generate widespread tephra fall in loads that will cause significant disruption and have severe socio-economic impacts. Past OVC tephra hazard studies have focused on evaluating hazard from a rhyolitic Plinian eruption at select vent locations in the OVC's Tarawera LVZ. Here, we expand upon past studies by evaluating tephra hazard for all possible OVC eruption vent areas and for both rhyolitic and basaltic Plinian eruption styles, and explore how these parameters influence tephra hazard forecasts. Probabilistic volcanic hazard model BET_VH and advection-diffusion model TEPHRA2 were used to assess the hazard of accumulating ≥10 kg m-2 of tephra from both basaltic Plinian and rhyolitic Plinian eruption styles, occurring from within the Tarawera LVZ, the Haroharo LVZ or other potential vent areas within the caldera. Our results highlight the importance of considering all the potential vent locations of a volcanic system, in order to capture the full eruption catalogue in analyses (e.g. 11 eruptions over 26 ky for the OVC versus only five eruptions over 26 ky for the Tarawera LVZ), as well as the full spatial distribution of tephra hazard. Although the Tarawera LVZ has been prominently discussed in studies of OVC hazard because of its recent activity (1886 and ~1315 ad), we find that in the event of a future eruption, the estimated likelihood of a vent opening within the Haroharo LVZ (last eruption 5.6 ka) is equivalent (<1 % difference) to that for the Tarawera LVZ (31.8 compared to 32.5 %). Including both the Haroharo LVZ and the Tarawera LVZ as possible source areas in the hazard analysis allows us to assess the full spatial extent of OVC tephra fall hazard. By considering both basaltic Plinian and rhyolitic Plinian eruption styles, as well as multiple vent location areas, we present a hazard assessment which aims to reduce bias through incorporating a greater range of eruption variables.

Meta-heuristic CRPS minimization for the calibration of short-range probabilistic forecasts

NASA Astrophysics Data System (ADS)

Mohammadi, Seyedeh Atefeh; Rahmani, Morteza; Azadi, Majid

2016-08-01

This paper deals with the probabilistic short-range temperature forecasts over synoptic meteorological stations across Iran using non-homogeneous Gaussian regression (NGR). NGR creates a Gaussian forecast probability density function (PDF) from the ensemble output. The mean of the normal predictive PDF is a bias-corrected weighted average of the ensemble members and its variance is a linear function of the raw ensemble variance. The coefficients for the mean and variance are estimated by minimizing the continuous ranked probability score (CRPS) during a training period. CRPS is a scoring rule for distributional forecasts. In the paper of Gneiting et al. (Mon Weather Rev 133:1098-1118, 2005), Broyden-Fletcher-Goldfarb-Shanno (BFGS) method is used to minimize the CRPS. Since BFGS is a conventional optimization method with its own limitations, we suggest using the particle swarm optimization (PSO), a robust meta-heuristic method, to minimize the CRPS. The ensemble prediction system used in this study consists of nine different configurations of the weather research and forecasting model for 48-h forecasts of temperature during autumn and winter 2011 and 2012. The probabilistic forecasts were evaluated using several common verification scores including Brier score, attribute diagram and rank histogram. Results show that both BFGS and PSO find the optimal solution and show the same evaluation scores, but PSO can do this with a feasible random first guess and much less computational complexity.

Forecasting volcanic eruptions and other material failure phenomena: An evaluation of the failure forecast method

NASA Astrophysics Data System (ADS)

Bell, Andrew F.; Naylor, Mark; Heap, Michael J.; Main, Ian G.

2011-08-01

Power-law accelerations in the mean rate of strain, earthquakes and other precursors have been widely reported prior to material failure phenomena, including volcanic eruptions, landslides and laboratory deformation experiments, as predicted by several theoretical models. The Failure Forecast Method (FFM), which linearizes the power-law trend, has been routinely used to forecast the failure time in retrospective analyses; however, its performance has never been formally evaluated. Here we use synthetic and real data, recorded in laboratory brittle creep experiments and at volcanoes, to show that the assumptions of the FFM are inconsistent with the error structure of the data, leading to biased and imprecise forecasts. We show that a Generalized Linear Model method provides higher-quality forecasts that converge more accurately to the eventual failure time, accounting for the appropriate error distributions. This approach should be employed in place of the FFM to provide reliable quantitative forecasts and estimate their associated uncertainties.

The quality and value of seasonal precipitation forecasts for an early warning of large-scale droughts and floods in West Africa

NASA Astrophysics Data System (ADS)

Bliefernicht, Jan; Seidel, Jochen; Salack, Seyni; Waongo, Moussa; Laux, Patrick; Kunstmann, Harald

2017-04-01

Seasonal precipitation forecasts are a crucial source of information for an early warning of hydro-meteorological extremes in West Africa. However, the current seasonal forecasting system used by the West African weather services in the framework of the West African Climate Outlook forum (PRESAO) is limited to probabilistic precipitation forecasts of 1-month lead time. To improve this provision, we use an ensemble-based quantile-quantile transformation for bias correction of precipitation forecasts provided by a global seasonal ensemble prediction system, the Climate Forecast System Version 2 (CFS2). The statistical technique eliminates systematic differences between global forecasts and observations with the potential to preserve the signal from the model. The technique has also the advantage that it can be easily implemented at national weather services with low capacities. The statistical technique is used to generate probabilistic forecasts of monthly and seasonal precipitation amount and other precipitation indices useful for an early warning of large-scale drought and floods in West Africa. The evaluation of the statistical technique is done using CFS hindcasts (1982 to 2009) in a cross-validation mode to determine the performance of the precipitation forecasts for several lead times focusing on drought and flood events depicted over the Volta and Niger basins. In addition, operational forecasts provided by PRESAO are analyzed from 1998 to 2015. The precipitation forecasts are compared to low-skill reference forecasts generated from gridded observations (i.e. GPCC, CHIRPS) and a novel in-situ gauge database from national observation networks (see Poster EGU2017-10271). The forecasts are evaluated using state-of-the-art verification techniques to determine specific quality attributes of probabilistic forecasts such as reliability, accuracy and skill. In addition, cost-loss approaches are used to determine the value of probabilistic forecasts for multiple users in warning situations. The outcomes of the hindcasts experiment for the Volta basin illustrate that the statistical technique can clearly improve the CFS precipitation forecasts with the potential to provide skillful and valuable early precipitation warnings for large-scale drought and flood situations several months in ahead. In this presentation we give a detailed overview about the ensemble-based quantile-quantile-transformation, its validation and verification and the possibilities of this technique to complement PRESAO. We also highlight the performance of this technique for extremes such as the Sahel drought in the 80ties and in comparison to the various reference data sets (e.g. CFS2, PRESAO, observational data sets) used in this study.

An online tool for Operational Probabilistic Drought Forecasting System (OPDFS): a Statistical-Dynamical Framework

NASA Astrophysics Data System (ADS)

Zarekarizi, M.; Moradkhani, H.; Yan, H.

2017-12-01

The Operational Probabilistic Drought Forecasting System (OPDFS) is an online tool recently developed at Portland State University for operational agricultural drought forecasting. This is an integrated statistical-dynamical framework issuing probabilistic drought forecasts monthly for the lead times of 1, 2, and 3 months. The statistical drought forecasting method utilizes copula functions in order to condition the future soil moisture values on the antecedent states. Due to stochastic nature of land surface properties, the antecedent soil moisture states are uncertain; therefore, data assimilation system based on Particle Filtering (PF) is employed to quantify the uncertainties associated with the initial condition of the land state, i.e. soil moisture. PF assimilates the satellite soil moisture data to Variable Infiltration Capacity (VIC) land surface model and ultimately updates the simulated soil moisture. The OPDFS builds on the NOAA's seasonal drought outlook by offering drought probabilities instead of qualitative ordinal categories and provides the user with the probability maps associated with a particular drought category. A retrospective assessment of the OPDFS showed that the forecasting of the 2012 Great Plains and 2014 California droughts were possible at least one month in advance. The OPDFS offers a timely assistance to water managers, stakeholders and decision-makers to develop resilience against uncertain upcoming droughts.

Naples between two fires: eruptive scenarios for the next eruptions by an integrated volcanological-probabilistic approach.

NASA Astrophysics Data System (ADS)

Mastrolorenzo, G.; Pappalardo, L.; de Natale, G.; Troise, C.; Rossano, S.; Panizza, A.

2009-04-01

Probabilistic approaches based on available volcanological data from real eruptions of Campi Flegrei and Somma-Vesuvius, are assembled in a comprehensive assessment of volcanic hazards at the Neapolitan area. This allows to compare the volcanic hazards related to the different types of events, which can be used for evaluating the conditional probability of flows and falls hazard in case of a volcanic crisis. Hazard maps are presented, based on a rather complete set of numerical simulations, produced using field and laboratory data as input parameters relative to a large range (VEI 1 to 5) of fallout and pyroclastic-flow events and their relative occurrence. The results allow us to quantitatively evaluate and compare the hazard related to pyroclastic fallout and density currents (PDCs) at the Neapolitan volcanoes and their surroundings, including the city of Naples. Due to its position between the two volcanic areas, the city of Naples is particularly exposed to volcanic risk from VEI>2 eruptions, as recorded in the local volcanic succession. Because dominant wind directions, the area of Naples is particularly prone to fallout hazard from Campi Flegrei caldera eruptions in the VEI range 2-5. The hazard from PDCs decreases roughly radially with distance from the eruptive vents and is strongly controlled by the topographic heights. Campi Flegrei eruptions are particularly hazardous for Naples, although the Camaldoli and Posillipo hills produce an effective barrier to propagation to the very central part of Naples. PDCs from Vesuvius eruptions with VEI>4 can cover the city of Naples, whereas even VEI>3 eruptions have a moderate fallout hazard there.

Evaluation of the Plant-Craig stochastic convection scheme (v2.0) in the ensemble forecasting system MOGREPS-R (24 km) based on the Unified Model (v7.3)

NASA Astrophysics Data System (ADS)

Keane, Richard J.; Plant, Robert S.; Tennant, Warren J.

2016-05-01

The Plant-Craig stochastic convection parameterization (version 2.0) is implemented in the Met Office Regional Ensemble Prediction System (MOGREPS-R) and is assessed in comparison with the standard convection scheme with a simple stochastic scheme only, from random parameter variation. A set of 34 ensemble forecasts, each with 24 members, is considered, over the month of July 2009. Deterministic and probabilistic measures of the precipitation forecasts are assessed. The Plant-Craig parameterization is found to improve probabilistic forecast measures, particularly the results for lower precipitation thresholds. The impact on deterministic forecasts at the grid scale is neutral, although the Plant-Craig scheme does deliver improvements when forecasts are made over larger areas. The improvements found are greater in conditions of relatively weak synoptic forcing, for which convective precipitation is likely to be less predictable.

Evaluating probabilistic dengue risk forecasts from a prototype early warning system for Brazil.

PubMed

Lowe, Rachel; Coelho, Caio As; Barcellos, Christovam; Carvalho, Marilia Sá; Catão, Rafael De Castro; Coelho, Giovanini E; Ramalho, Walter Massa; Bailey, Trevor C; Stephenson, David B; Rodó, Xavier

2016-02-24

Recently, a prototype dengue early warning system was developed to produce probabilistic forecasts of dengue risk three months ahead of the 2014 World Cup in Brazil. Here, we evaluate the categorical dengue forecasts across all microregions in Brazil, using dengue cases reported in June 2014 to validate the model. We also compare the forecast model framework to a null model, based on seasonal averages of previously observed dengue incidence. When considering the ability of the two models to predict high dengue risk across Brazil, the forecast model produced more hits and fewer missed events than the null model, with a hit rate of 57% for the forecast model compared to 33% for the null model. This early warning model framework may be useful to public health services, not only ahead of mass gatherings, but also before the peak dengue season each year, to control potentially explosive dengue epidemics.

On a Possible Unified Scaling Law for Volcanic Eruption Durations

PubMed Central

Cannavò, Flavio; Nunnari, Giuseppe

2016-01-01

Volcanoes constitute dissipative systems with many degrees of freedom. Their eruptions are the result of complex processes that involve interacting chemical-physical systems. At present, due to the complexity of involved phenomena and to the lack of precise measurements, both analytical and numerical models are unable to simultaneously include the main processes involved in eruptions thus making forecasts of volcanic dynamics rather unreliable. On the other hand, accurate forecasts of some eruption parameters, such as the duration, could be a key factor in natural hazard estimation and mitigation. Analyzing a large database with most of all the known volcanic eruptions, we have determined that the duration of eruptions seems to be described by a universal distribution which characterizes eruption duration dynamics. In particular, this paper presents a plausible global power-law distribution of durations of volcanic eruptions that holds worldwide for different volcanic environments. We also introduce a new, simple and realistic pipe model that can follow the same found empirical distribution. Since the proposed model belongs to the family of the self-organized systems it may support the hypothesis that simple mechanisms can lead naturally to the emergent complexity in volcanic behaviour. PMID:26926425

On a Possible Unified Scaling Law for Volcanic Eruption Durations.

PubMed

Cannavò, Flavio; Nunnari, Giuseppe

2016-03-01

Volcanoes constitute dissipative systems with many degrees of freedom. Their eruptions are the result of complex processes that involve interacting chemical-physical systems. At present, due to the complexity of involved phenomena and to the lack of precise measurements, both analytical and numerical models are unable to simultaneously include the main processes involved in eruptions thus making forecasts of volcanic dynamics rather unreliable. On the other hand, accurate forecasts of some eruption parameters, such as the duration, could be a key factor in natural hazard estimation and mitigation. Analyzing a large database with most of all the known volcanic eruptions, we have determined that the duration of eruptions seems to be described by a universal distribution which characterizes eruption duration dynamics. In particular, this paper presents a plausible global power-law distribution of durations of volcanic eruptions that holds worldwide for different volcanic environments. We also introduce a new, simple and realistic pipe model that can follow the same found empirical distribution. Since the proposed model belongs to the family of the self-organized systems it may support the hypothesis that simple mechanisms can lead naturally to the emergent complexity in volcanic behaviour.

The climatic effect of explosive volcanic activity: Analysis of the historical data

NASA Technical Reports Server (NTRS)

Bryson, R. A.; Goodman, B. M.

1982-01-01

By using the most complete available records of direct beam radiation and volcanic eruptions, an historical analysis of the role of the latter in modulating the former was made. A very simple fallout and dispersion model was applied to the historical chronology of explosive eruptions. The resulting time series explains about 77 percent of the radiation variance, as well as suggests that tropical and subpolar eruptions are more important than mid-latitude eruptions in their impact on the stratospheric aerosol optical depth. The simpler climatic models indicate that past hemispheric temperature can be stimulated very well with volcanic and CO2 inputs and suggest that climate forecasting will also require volcano forecasting. There is some evidence that this is possible some years in advance.

New insights into the Kawah Ijen hydrothermal system from geophysical data

USGS Publications Warehouse

Caudron, Corentin; Mauri, G.; Williams-Jones, Glyn; Lecocq, Thomas; Syahbana, Devy Kamil; de Plaen, Raphael; Peiffer, Loic; Bernard, Alain; Saracco, Ginette

2017-01-01

Volcanoes with crater lakes and/or extensive hydrothermal systems pose significant challenges with respect to monitoring and forecasting eruptions, but they also provide new opportunities to enhance our understanding of magmatic–hydrothermal processes. Their lakes and hydrothermal systems serve as reservoirs for magmatic heat and fluid emissions, filtering and delaying the surface expressions of magmatic unrest and eruption, yet they also enable sampling and monitoring of geochemical tracers. Here, we describe the outcomes of a highly focused international experimental campaign and workshop carried out at Kawah Ijen volcano, Indonesia, in September 2014, designed to answer fundamental questions about how to improve monitoring and eruption forecasting at wet volcanoes.

The game of making decisions under uncertainty: How sure must one be?

NASA Astrophysics Data System (ADS)

Werner, Micha; Verkade, Jan; Wetterhall, Fredrik; van Andel, Schalk-Jan; Ramos, Maria-Helena

2016-04-01

Probabilistic hydrometeorological forecasting is now widely accepted to be more skillful than deterministic forecasts, and is increasingly being integrated into operational practice. Provided they are reliable and unbiased, probabilistic forecasts have the advantage that they give decision maker not only the forecast value, but also the uncertainty associated to that prediction. Though that information provides more insight, it does now leave the forecaster/decision maker with the challenge of deciding at what level of probability of a threshold being exceeded the decision to act should be taken. According to the cost-loss theory, that probability should be related to the impact of the threshold being exceeded. However, it is not entirely clear how easy it is for decision makers to follow that rule, even when the impact of a threshold being exceeded, and the actions to choose from are known. To continue the tradition in the "Ensemble Hydrometeorological Forecast" session, we will address the challenge of making decisions based on probabilistic forecasts through a game to be played with the audience. We will explore how decisions made differ depending on the known impacts of the forecasted events. Participants will be divided into a number of groups with differing levels of impact, and will be faced with a number of forecast situations. They will be asked to make decisions and record the consequence of those decisions. A discussion of the differences in the decisions made will be presented at the end of the game, with a fuller analysis later posted on the HEPEX web site blog (www.hepex.org).

Towards an improved ensemble precipitation forecast: A probabilistic post-processing approach

NASA Astrophysics Data System (ADS)

Khajehei, Sepideh; Moradkhani, Hamid

2017-03-01

Recently, ensemble post-processing (EPP) has become a commonly used approach for reducing the uncertainty in forcing data and hence hydrologic simulation. The procedure was introduced to build ensemble precipitation forecasts based on the statistical relationship between observations and forecasts. More specifically, the approach relies on a transfer function that is developed based on a bivariate joint distribution between the observations and the simulations in the historical period. The transfer function is used to post-process the forecast. In this study, we propose a Bayesian EPP approach based on copula functions (COP-EPP) to improve the reliability of the precipitation ensemble forecast. Evaluation of the copula-based method is carried out by comparing the performance of the generated ensemble precipitation with the outputs from an existing procedure, i.e. mixed type meta-Gaussian distribution. Monthly precipitation from Climate Forecast System Reanalysis (CFS) and gridded observation from Parameter-Elevation Relationships on Independent Slopes Model (PRISM) have been employed to generate the post-processed ensemble precipitation. Deterministic and probabilistic verification frameworks are utilized in order to evaluate the outputs from the proposed technique. Distribution of seasonal precipitation for the generated ensemble from the copula-based technique is compared to the observation and raw forecasts for three sub-basins located in the Western United States. Results show that both techniques are successful in producing reliable and unbiased ensemble forecast, however, the COP-EPP demonstrates considerable improvement in the ensemble forecast in both deterministic and probabilistic verification, in particular in characterizing the extreme events in wet seasons.

Predictability of short-range forecasting: a multimodel approach

NASA Astrophysics Data System (ADS)

García-Moya, Jose-Antonio; Callado, Alfons; Escribà, Pau; Santos, Carlos; Santos-Muñoz, Daniel; Simarro, Juan

2011-05-01

Numerical weather prediction (NWP) models (including mesoscale) have limitations when it comes to dealing with severe weather events because extreme weather is highly unpredictable, even in the short range. A probabilistic forecast based on an ensemble of slightly different model runs may help to address this issue. Among other ensemble techniques, Multimodel ensemble prediction systems (EPSs) are proving to be useful for adding probabilistic value to mesoscale deterministic models. A Multimodel Short Range Ensemble Prediction System (SREPS) focused on forecasting the weather up to 72 h has been developed at the Spanish Meteorological Service (AEMET). The system uses five different limited area models (LAMs), namely HIRLAM (HIRLAM Consortium), HRM (DWD), the UM (UKMO), MM5 (PSU/NCAR) and COSMO (COSMO Consortium). These models run with initial and boundary conditions provided by five different global deterministic models, namely IFS (ECMWF), UM (UKMO), GME (DWD), GFS (NCEP) and CMC (MSC). AEMET-SREPS (AE) validation on the large-scale flow, using ECMWF analysis, shows a consistent and slightly underdispersive system. For surface parameters, the system shows high skill forecasting binary events. 24-h precipitation probabilistic forecasts are verified using an up-scaling grid of observations from European high-resolution precipitation networks, and compared with ECMWF-EPS (EC).

Seismic energy data analysis of Merapi volcano to test the eruption time prediction using materials failure forecast method (FFM)

NASA Astrophysics Data System (ADS)

Anggraeni, Novia Antika

2015-04-01

The test of eruption time prediction is an effort to prepare volcanic disaster mitigation, especially in the volcano's inhabited slope area, such as Merapi Volcano. The test can be conducted by observing the increase of volcanic activity, such as seismicity degree, deformation and SO2 gas emission. One of methods that can be used to predict the time of eruption is Materials Failure Forecast Method (FFM). Materials Failure Forecast Method (FFM) is a predictive method to determine the time of volcanic eruption which was introduced by Voight (1988). This method requires an increase in the rate of change, or acceleration of the observed volcanic activity parameters. The parameter used in this study is the seismic energy value of Merapi Volcano from 1990 - 2012. The data was plotted in form of graphs of seismic energy rate inverse versus time with FFM graphical technique approach uses simple linear regression. The data quality control used to increase the time precision employs the data correlation coefficient value of the seismic energy rate inverse versus time. From the results of graph analysis, the precision of prediction time toward the real time of eruption vary between -2.86 up to 5.49 days.

Seismic energy data analysis of Merapi volcano to test the eruption time prediction using materials failure forecast method (FFM)

DOE Office of Scientific and Technical Information (OSTI.GOV)

Anggraeni, Novia Antika, E-mail: novia.antika.a@gmail.com

The test of eruption time prediction is an effort to prepare volcanic disaster mitigation, especially in the volcano’s inhabited slope area, such as Merapi Volcano. The test can be conducted by observing the increase of volcanic activity, such as seismicity degree, deformation and SO2 gas emission. One of methods that can be used to predict the time of eruption is Materials Failure Forecast Method (FFM). Materials Failure Forecast Method (FFM) is a predictive method to determine the time of volcanic eruption which was introduced by Voight (1988). This method requires an increase in the rate of change, or acceleration ofmore » the observed volcanic activity parameters. The parameter used in this study is the seismic energy value of Merapi Volcano from 1990 – 2012. The data was plotted in form of graphs of seismic energy rate inverse versus time with FFM graphical technique approach uses simple linear regression. The data quality control used to increase the time precision employs the data correlation coefficient value of the seismic energy rate inverse versus time. From the results of graph analysis, the precision of prediction time toward the real time of eruption vary between −2.86 up to 5.49 days.« less

Eruption-induced modifications to volcanic seismicity at Ruapehu, New Zealand, and its implications for eruption forecasting

USGS Publications Warehouse

Bryan, C.J.; Sherburn, S.

2003-01-01

Broadband seismic data collected on Ruapehu volcano, New Zealand, in 1994 and 1998 show that the 1995-1996 eruptions of Ruapehu resulted in a significant change in the frequency content of tremor and volcanic earthquakes at the volcano. The pre-eruption volcanic seismicity was characterized by several independent dominant frequencies, with a 2 Hz spectral peak dominating the strongest tremor and volcanic earthquakes and higher frequencies forming the background signal. The post-eruption volcanic seismicity was dominated by a 0.8-1.4 Hz spectral peak not seen before the eruptions. The 2 Hz and higher frequency signals remained, but were subordinate to the 0.8-1.4 Hz energy. That the dominant frequencies of volcanic tremor and volcanic earthquakes were identical during the individual time periods prior to and following the 1995-1996 eruptions suggests that during each of these time periods the volcanic tremor and earthquakes were generated by the same source process. The overall change in the frequency content, which occurred during the 1995-1996 eruptions and remains as of the time of the writing of this paper, most likely resulted from changes in the volcanic plumbing system and has significant implications for forecasting and real-time assessment of future eruptive activity at Ruapehu.

Changes in shear-wave splitting before volcanic eruptions

NASA Astrophysics Data System (ADS)

Liu, Sha; Crampin, Stuart

2015-04-01

We have shown that observations of shear-wave splitting (SWS) monitor stress-accumulation and stress-relaxation before earthquakes which allows the time, magnitude, and in some circumstances fault-plane of impending earthquakes to be stress-forecast. (We call this procedure stress-forecasting rather than predicting or forecasting to emphasise the different formalism.) We have stress-forecast these parameters successfully three-days before a 1988 M5 earthquake in SW Iceland, and identified characteristic anomalies retrospectively before ~16 other earthquakes in Iceland and elsewhere. SWS monitors microcrack geometry and shows that microcracks are so closely spaced that they verge on fracturing and earthquakes. Phenomena verging on failure in this way are critical-systems with 'butterfly wings' sensitivity. Such critical-systems are very common. The Earth is an archetypal complex heterogeneous interactive phenomenon and must be expected to be a critical-system. We claim this critical system as a New Geophysics of a critically-microcracked rock mass. Such critical systems impose a range of fundamentally-new properties on conventional sub-critical physics/geophysics, one of which is universality. Consequently it is expected that we observe similar stress-accumulation and stress-relaxation before volcanic eruptions to those before earthquakes. There are three eruptions where appropriate changes in SWS have been observed similar to those observed before earthquakes. These are: the 1996 Gjálp fissure eruption, Vatnajökull, Iceland; a 2001 flank eruption on Mount Etna, Sicily (reported by Francesca Bianco, INGV, Naples); and the 2010 Eyjafjajökull ash-cloud eruption, SW Iceland. These will be presented in the same normalised format as is used before earthquakes. The 1996 Gjálp eruption showed a 2½-month stress-accumulation, and a ~1-year stress-relaxation (attributed to the North Atlantic Ridge adjusting to the magma injection beneath the Vatnajökull Ice Cap). The 2001 flank eruption of Etna showed stress-accumulation and stress-relaxation typical of a small earthquake. However, the changes in SWS before the 2010 Eyjafjajökull Eruption, SW Iceland, showed the most distinctive correlations with earthquakes, as it was only ~90km-west of the 1988 M5 in SW Iceland, which was successfully stress-forecast. The behaviour of SWS before the M5 earthquake and the Eyjafjajökull flank (ash-cloud) eruption is almost identical both showing linear stress-accumulation increases, and linear stress-relaxation decreases to the earthquake and the onset of the flank eruption, respectively. There are comparable slopes and durations. We consider this strong confirmation of the universality property of the New Geophysics of a critically-microcracked Earth. Papers referring to these developments can be found in geos.ed.ac.uk/home/scrampin/opinion. Also see abstracts in EGU2015 Sessions: Crampin & Gao (SM1.1), Gao & Crampin (SM3.1), and Crampin & Gao (GD.1).

Identifying open and closed system behaviors at Tungurahua volcano (Ecuador) using SO2 and seismo-acoustic measurements

NASA Astrophysics Data System (ADS)

Hidalgo, Silvana; Battaglia, Jean; Bernard, Benjamin; Steele, Alexander; Arellano, Santiago; Galle, Bo

2014-05-01

Tungurahua is one of the most active volcanoes in Ecuador. It is located in Central Ecuador, 160 km South of Quito and 8 km South of the touristic town of Baños. Tungurahua had one eruption every century since 1500, with an activity characterized by ash fallouts and pyroclastic and lava flows. The current eruptive period of Tungurahua began in 1999 with multiple episodes of explosive activity that have threatened the local population. The monitoring network is constituted by 5 short period and 5 broadband seismic stations, 4 DOAS permanent instruments, 4 tiltmeters, 2 permanent high resolution GPS, 4 digital cameras and 10 acoustic flow monitors. The correct interpretation of the different data acquired by this network allows a better understanding of the eruptive behavior of Tungurahua in order to provide early warning to the local population. Tungurahua changed its behavior from a continuously erupting volcano, as it was until 2008, to a sporadically erupting one, showing clear quiescence phases lasting from 40 to 184 days, and intense activity phases lasting from 15 to 70 days. Activity phases are characterized by Strombolian and Vulcanian eruptive styles, producing ash fallouts and in a few occasions pyroclastic flows. In terms of hazard to the local population, one of the goals of monitoring Tungurahura is to forecast the onset and evolution of eruptive phases. In particular the occurrence of large Vulcanian explosions which occur when the conduit is closed is a major issue. Since 2010 we focused our study on the relation between SO2 gas emissions, the seismic and acoustic energies of explosions and the tremor amplitudes. The first observation of comparing these different datasets is that the correlation between seismic and SO2 degassing is not straightforward, and actually the relation reflects the conditions at the vent: open or closed. The onset of eruptive phases in open conduit conditions can be identified which leads to an effective eruption forecasting. An example of this behavior is the eruptive phase between December 2009 and March 2010 when SO2 measurements increased 4 days before the amplitude of tremor and 9 days before the occurrence of the first explosions. Conversely, if the vent is closed at the beginning of a phase and no evident seismic precursors are observed forecasting is hardly possible. During an ongoing eruptive phase, the relation between these parameters allows to identify periods when the conduit is totally open as degassing may occur almost without generating any seismicity. Therefore the forecasting of escalating open conduit activity or a partial closing of the system is possible. Such a case was observed and forecasted on December 2011. In this work, we present observational evidence of these mechanisms which are used to identify possible patterns of evolution of the activity, contributing to a more effective volcanic hazard assessment.

Leveraging Past and Current Measurements to Probabilistically Nowcast Low Visibility Procedures at an Airport

NASA Astrophysics Data System (ADS)

Mayr, G. J.; Kneringer, P.; Dietz, S. J.; Zeileis, A.

2016-12-01

Low visibility or low cloud ceiling reduce the capacity of airports by requiring special low visibility procedures (LVP) for incoming/departing aircraft. Probabilistic forecasts when such procedures will become necessary help to mitigate delays and economic losses.We compare the performance of probabilistic nowcasts with two statistical methods: ordered logistic regression, and trees and random forests. These models harness historic and current meteorological measurements in the vicinity of the airport and LVP states, and incorporate diurnal and seasonal climatological information via generalized additive models (GAM). The methods are applied at Vienna International Airport (Austria). The performance is benchmarked against climatology, persistence and human forecasters.

Topics in Probabilistic Judgment Aggregation

ERIC Educational Resources Information Center

Wang, Guanchun

2011-01-01

This dissertation is a compilation of several studies that are united by their relevance to probabilistic judgment aggregation. In the face of complex and uncertain events, panels of judges are frequently consulted to provide probabilistic forecasts, and aggregation of such estimates in groups often yield better results than could have been made…

Probabilistic Hazard Estimation at a Densely Urbanised Area: the Neaples Volcanoes

NASA Astrophysics Data System (ADS)

de Natale, G.; Mastrolorenzo, G.; Panizza, A.; Pappalardo, L.; Claudia, T.

2005-12-01

The Neaples volcanic area (Southern Italy), including Vesuvius, Campi Flegrei caldera and Ischia island, is the highest risk one in the World, where more than 2 million people live within about 10 km from an active volcanic vent. Such an extreme risk calls for accurate methodologies aimed to quantify it, in a probabilistic way, considering all the available volcanological information as well as modelling results. In fact, simple hazard maps based on the observation of deposits from past eruptions have the major problem that eruptive history generally samples a very limited number of possible outcomes, thus resulting almost meaningless to get the event probability in the area. This work describes a methodology making the best use (from a Bayesian point of view) of volcanological data and modelling results, to compute probabilistic hazard maps from multi-vent explosive eruptions. The method, which follows an approach recently developed by the same authors for pyroclastic flows hazard, has been here improved and extended to compute also fall-out hazard. The application of the method to the Neapolitan volcanic area, including the densely populated city of Naples, allows, for the first time, to get a global picture of the areal distribution for the main hazards from multi-vent explosive eruptions. From a joint consideration of the hazard contributions from all the three volcanic areas, new insight on the volcanic hazard distribution emerges, which will have strong implications for urban and emergency planning in the area.

Enhancing Community Based Early Warning Systems in Nepal with Flood Forecasting Using Local and Global Models

NASA Astrophysics Data System (ADS)

Dugar, Sumit; Smith, Paul; Parajuli, Binod; Khanal, Sonu; Brown, Sarah; Gautam, Dilip; Bhandari, Dinanath; Gurung, Gehendra; Shakya, Puja; Kharbuja, RamGopal; Uprety, Madhab

2017-04-01

Operationalising effective Flood Early Warning Systems (EWS) in developing countries like Nepal poses numerous challenges, with complex topography and geology, sparse network of river and rainfall gauging stations and diverse socio-economic conditions. Despite these challenges, simple real-time monitoring based EWSs have been in place for the past decade. A key constraint of these simple systems is the very limited lead time for response - as little as 2-3 hours, especially for rivers originating from steep mountainous catchments. Efforts to increase lead time for early warning are focusing on imbedding forecasts into the existing early warning systems. In 2016, the Nepal Department of Hydrology and Meteorology (DHM) piloted an operational Probabilistic Flood Forecasting Model in major river basins across Nepal. This comprised a low data approach to forecast water levels, developed jointly through a research/practitioner partnership with Lancaster University and WaterNumbers (UK) and the International NGO Practical Action. Using Data-Based Mechanistic Modelling (DBM) techniques, the model assimilated rainfall and water levels to generate localised hourly flood predictions, which are presented as probabilistic forecasts, increasing lead times from 2-3 hours to 7-8 hours. The Nepal DHM has simultaneously started utilizing forecasts from the Global Flood Awareness System (GLoFAS) that provides streamflow predictions at the global scale based upon distributed hydrological simulations using numerical ensemble weather forecasts from the ECMWF (European Centre for Medium-Range Weather Forecasts). The aforementioned global and local models have already affected the approach to early warning in Nepal, being operational during the 2016 monsoon in the West Rapti basin in Western Nepal. On 24 July 2016, GLoFAS hydrological forecasts for the West Rapti indicated a sharp rise in river discharge above 1500 m3/sec (equivalent to the river warning level at 5 meters) with 53% probability of exceeding the Medium Level Alert in two days. Rainfall stations upstream of the West Rapti catchment recorded heavy rainfall on 26 July, and localized forecasts from the probabilistic model at 8 am suggested that the water level would cross a pre-determined warning level in the next 3 hours. The Flood Forecasting Section at DHM issued a flood advisory, and disseminated SMS flood alerts to more than 13,000 at-risk people residing along the floodplains. Water levels crossed the danger threshold (5.4 meters) at 11 am, peaking at 8.15 meters at 10 pm. Extension of the warning lead time from probabilistic forecasts was significant in minimising the risk to lives and livelihoods as communities gained extra time to prepare, evacuate and respond. Likewise, longer timescale forecasts from GLoFAS could be potentially linked with no-regret early actions leading to improved preparedness and emergency response. These forecasting tools have contributed to enhance the effectiveness and efficiency of existing community based systems, increasing the lead time for response. Nevertheless, extensive work is required on appropriate ways to interpret and disseminate probabilistic forecasts having longer (2-14 days) and shorter (3-5 hours) time horizon for operational deployment as there are numerous uncertainties associated with predictions.

Towards a Proactive Risk Mitigation Strategy at La Fossa Volcano, Vulcano Island

NASA Astrophysics Data System (ADS)

Biass, S.; Gregg, C. E.; Frischknecht, C.; Falcone, J. L.; Lestuzzi, P.; di Traglia, F.; Rosi, M.; Bonadonna, C.

2014-12-01

A comprehensive risk assessment framework was built to develop proactive risk reduction measures for Vulcano Island, Italy. This framework includes identification of eruption scenarios; probabilistic hazard assessment, quantification of hazard impacts on the built environment, accessibility assessment on the island and risk perception study. Vulcano, a 21 km2 island with two primary communities host to 900 permanent residents and up to 10,000 visitors during summer, shows a strong dependency on the mainland for basic needs (water, energy) and relies on a ~2 month tourism season for its economy. The recent stratigraphy reveals a dominance of vulcanian and subplinian eruptions, producing a range of hazards acting at different time scales. We developed new methods to probabilistically quantify the hazard related to ballistics, lahars and tephra for all eruption styles. We also elaborated field- and GIS- based methods to assess the physical vulnerability of the built environment and created dynamic models of accessibility. Results outline the difference of hazard between short and long-lasting eruptions. A subplinian eruption has a 50% probability of impacting ~30% of the buildings within days after the eruption, but the year-long damage resulting from a long-lasting vulcanian eruption is similar if tephra is not removed from rooftops. Similarly, a subplinian eruption results in a volume of 7x105 m3 of material potentially remobilized into lahars soon after the eruption. Similar volumes are expected for a vulcanian activity over years, increasing the hazard of small lahars. Preferential lahar paths affect critical infrastructures lacking redundancy, such as the road network, communications systems, the island's only gas station, and access to the island's two evacuation ports. Such results from hazard, physical and systemic vulnerability help establish proactive volcanic risk mitigation strategies and may be applicable in other island settings.

Calibration of decadal ensemble predictions

NASA Astrophysics Data System (ADS)

Pasternack, Alexander; Rust, Henning W.; Bhend, Jonas; Liniger, Mark; Grieger, Jens; Müller, Wolfgang; Ulbrich, Uwe

2017-04-01

Decadal climate predictions are of great socio-economic interest due to the corresponding planning horizons of several political and economic decisions. Due to uncertainties of weather and climate, forecasts (e.g. due to initial condition uncertainty), they are issued in a probabilistic way. One issue frequently observed for probabilistic forecasts is that they tend to be not reliable, i.e. the forecasted probabilities are not consistent with the relative frequency of the associated observed events. Thus, these kind of forecasts need to be re-calibrated. While re-calibration methods for seasonal time scales are available and frequently applied, these methods still have to be adapted for decadal time scales and its characteristic problems like climate trend and lead time dependent bias. Regarding this, we propose a method to re-calibrate decadal ensemble predictions that takes the above mentioned characteristics into account. Finally, this method will be applied and validated to decadal forecasts from the MiKlip system (Germany's initiative for decadal prediction).

Multidimensional analysis and probabilistic model of volcanic and seismic activities

NASA Astrophysics Data System (ADS)

Fedorov, V.

2009-04-01

A search for space and time regularities in volcanic and seismic events for the purpose of forecast method development seems to be of current concern, both scientifically and practically. The seismic and volcanic processes take place in the Earth's field of gravity which in turn is closely related to gravitational fields of the Moon, the Sun, and the planets of the Solar System. It is mostly gravity and tidal forces that exercise control over the Earth's configuration and relief. Dynamic gravitational interaction between the Earth and other celestial bodies makes itself evident in tidal phenomena and other effects in the geospheres (including the Earth's crust). Dynamics of the tidal and attractive forces is responsible for periodical changes in gravity force, both in value and direction [Darwin, 1965], in the rate of rotation and orbital speed; that implies related changes in the endogenic activity of the Earth. The Earth's rotation in the alternating gravitational field accounts to a considerable extent for regular pattern of crustal deformations and dislocations; it is among principal factors that control the Earth's form and structure, distribution of oceans and continents and, probably, continental drift [Peive, 1969; Khain, 1973; Kosygin, 1983]. The energy of gravitational interaction is transmitted through the tidal energy to planetary spheres and feeds various processes there, including volcanic and seismic ones. To determine degree, character and special features of tidal force contribution to the volcanic and seismic processes is of primary importance for understanding of genetic and dynamic aspects of volcanism and seismicity. Both volcanic and seismic processes are involved in evolution of celestial bodies; they are operative on the planets of the Earth group and many satellites [Essays…, 1981; Lukashov, 1996]. From this standpoint, studies of those processes are essential with a view to development of scenarios of the Earth's evolution as a celestial body, as well as to forecast of changes in its relief. As the volcanic and seismic processes are of cosmic nature and occurrence, it seems logical to investigate their chronological structure in terms of astronomical time reference system or in parameters of the Earth orbital movement. Gravitational interaction of the Earth with the moon, the Sun and planets of the Solar system forms the physical basis of this multidimensional system; it manifests itself in tidal deformations of the Earth's lithosphere and in periodical changes in the planet rotation and orbital speed. A search for chronological correlation between the Earth's volcanism and seismicity on one hand and the orbital parameters dynamic on the other shows a certain promise in relation to prognostic decisions. It should be kept in mind that the calculation of astronomical characteristics (Ephemerides), which is one of the main lines in theoretical astronomy, spans many years both in the past and in future. It seems appropriate therefore to apply the astronomical time reference system to investigations of chronological structure of volcanic and seismic processes from the methodical viewpoint, as well as for retrospective and prognostic analyses. To investigate temporal pattern of the volcanic and seismic processes and to find a degree of their dependence on tidal forces, we used the astronomical time reference system as related to the Earth's orbital movement. The system is based on substitution of calendar dates of eruption and earthquakes for corresponding values of known astronomical characteristics, such as the Earth to Sun and Earth to Moon distances, ecliptic latitude of the Moon, etc. In coordinates of astronomical parameters (JPL Planetary and Lunar Efemerides, 1997, as compiled by the Jet Propulsion Laboratory, California Institute of Technology, on the basis of DE 406 block developed by NASA), we analyzed massifs of information, both volcanological (Catalogue of the World volcanic eruptions by I.I. Gushchenko, 1979) and seismological (database of USGS/NEIC Significant Worldwide Earthquakes, 2150 B.C.- 1994 A.D.) information which displays dynamics of endogenic relief-forming processes over a period of 1900 to 1994. In the course of the analysis, a substitution of calendar variable by a corresponding astronomical one has been performed and the epoch superposition method was applied. In essence, the method consists in that the massifs of information on volcanic eruptions (over a period of 1900 to 1977) and seismic events (1900-1994) are differentiated with respect to value of astronomical parameters which correspond to the calendar dates of the known eruptions and earthquakes, regardless of the calendar year. The obtained spectra of volcanic eruptions and violent earthquake distribution in the fields of the Earth orbital movement parameters were used as a basis for calculation of frequency spectra and diurnal probability of volcanic and seismic activity. The objective of the proposed investigations is a probabilistic model development of the volcanic and seismic events, as well as GIS designing for monitoring and forecast of volcanic and seismic activities. In accordance with the stated objective, three probability parameters have been found in the course of preliminary studies; they form the basis for GIS-monitoring and forecast development. 1. A multidimensional analysis of volcanic eruption and earthquakes (of magnitude 7) have been performed in terms of the Earth orbital movement. Probability characteristics of volcanism and seismicity have been defined for the Earth as a whole. Time intervals have been identified with a diurnal probability twice as great as the mean value. Diurnal probability of volcanic and seismic events has been calculated up to 2020. 2. A regularity is found in duration of dormant (repose) periods has been established. A relationship has been found between the distribution of the repose period probability density and duration of the period. 3. Features of spatial distribution of volcanic eruptions and earthquakes of magnitude 7 were analyzed, and those related to the Earth rotation identified. Frequencies of their spatial distribution are calculated. Using those parameters as the base, a scheme (algorithm) of probabilistic monitoring (long-range forecast) has been developed for volcanic and seismic events. Refereces (in Russian): 1. Fedorov V.M. Gravitational factors and astronomy-based chronology of processes in geospheres. Moscow University Publishing House, 2000. 368 p. 2. Fedorov V.M. Comparison between chronology of the Earth volcanic activity and characteristics of its orbital motion // Vulkanologiya i seismologiya, № 5, 2001, p. 65-67. 3. Fedorov V.M. Specific features of latitudinal distribution of volcanic eruptions// Vulkanologiya i seismologiya, № 4, 2002, p.39-43. 4. Fedorov V.M. Specific features of latitudinal distribution of endogenic relief-forming processes and the rotation of the Earth // Geomorphologiya, № 1, 2003, p.3-9. 5. Fedorov V.M. Comparison between chronology of the Earth volcanic and seismic activity and characteristics of its orbital motion // Izvestiya RAS. Ser. Geogr. № 5, 2003, p.16-20. 6. Fedorov V.M. Chronological structure and probability of volcanic events as related to tidal deformation of lithosphere // Vulkanologiya i seismologiya, № 1, 2005, p.44-50. 7. Fedorov V.M. Multidimensional analysis and a probabilistic model of the activity of endogenic relief-forming processes // Geomorphology, № 2, 2007, p. 37 - 48. 8. Fedorov V.M. Multidimensional analysis - is a spatiotemporal structure of the geodynamic activity of Earth// Vestnik Moskovskogo Universiteta; Ser. 4. Geology, № 4, 2007, p. 24-31.

Testing for ontological errors in probabilistic forecasting models of natural systems

PubMed Central

Marzocchi, Warner; Jordan, Thomas H.

2014-01-01

Probabilistic forecasting models describe the aleatory variability of natural systems as well as our epistemic uncertainty about how the systems work. Testing a model against observations exposes ontological errors in the representation of a system and its uncertainties. We clarify several conceptual issues regarding the testing of probabilistic forecasting models for ontological errors: the ambiguity of the aleatory/epistemic dichotomy, the quantification of uncertainties as degrees of belief, the interplay between Bayesian and frequentist methods, and the scientific pathway for capturing predictability. We show that testability of the ontological null hypothesis derives from an experimental concept, external to the model, that identifies collections of data, observed and not yet observed, that are judged to be exchangeable when conditioned on a set of explanatory variables. These conditional exchangeability judgments specify observations with well-defined frequencies. Any model predicting these behaviors can thus be tested for ontological error by frequentist methods; e.g., using P values. In the forecasting problem, prior predictive model checking, rather than posterior predictive checking, is desirable because it provides more severe tests. We illustrate experimental concepts using examples from probabilistic seismic hazard analysis. Severe testing of a model under an appropriate set of experimental concepts is the key to model validation, in which we seek to know whether a model replicates the data-generating process well enough to be sufficiently reliable for some useful purpose, such as long-term seismic forecasting. Pessimistic views of system predictability fail to recognize the power of this methodology in separating predictable behaviors from those that are not. PMID:25097265

Forecasting magma-chamber rupture at Santorini volcano, Greece.

PubMed

Browning, John; Drymoni, Kyriaki; Gudmundsson, Agust

2015-10-28

How much magma needs to be added to a shallow magma chamber to cause rupture, dyke injection, and a potential eruption? Models that yield reliable answers to this question are needed in order to facilitate eruption forecasting. Development of a long-lived shallow magma chamber requires periodic influx of magmas from a parental body at depth. This redistribution process does not necessarily cause an eruption but produces a net volume change that can be measured geodetically by inversion techniques. Using continuum-mechanics and fracture-mechanics principles, we calculate the amount of magma contained at shallow depth beneath Santorini volcano, Greece. We demonstrate through structural analysis of dykes exposed within the Santorini caldera, previously published data on the volume of recent eruptions, and geodetic measurements of the 2011-2012 unrest period, that the measured 0.02% increase in volume of Santorini's shallow magma chamber was associated with magmatic excess pressure increase of around 1.1 MPa. This excess pressure was high enough to bring the chamber roof close to rupture and dyke injection. For volcanoes with known typical extrusion and intrusion (dyke) volumes, the new methodology presented here makes it possible to forecast the conditions for magma-chamber failure and dyke injection at any geodetically well-monitored volcano.

Evaluating probabilistic dengue risk forecasts from a prototype early warning system for Brazil

PubMed Central

Lowe, Rachel; Coelho, Caio AS; Barcellos, Christovam; Carvalho, Marilia Sá; Catão, Rafael De Castro; Coelho, Giovanini E; Ramalho, Walter Massa; Bailey, Trevor C; Stephenson, David B; Rodó, Xavier

2016-01-01

Recently, a prototype dengue early warning system was developed to produce probabilistic forecasts of dengue risk three months ahead of the 2014 World Cup in Brazil. Here, we evaluate the categorical dengue forecasts across all microregions in Brazil, using dengue cases reported in June 2014 to validate the model. We also compare the forecast model framework to a null model, based on seasonal averages of previously observed dengue incidence. When considering the ability of the two models to predict high dengue risk across Brazil, the forecast model produced more hits and fewer missed events than the null model, with a hit rate of 57% for the forecast model compared to 33% for the null model. This early warning model framework may be useful to public health services, not only ahead of mass gatherings, but also before the peak dengue season each year, to control potentially explosive dengue epidemics. DOI: http://dx.doi.org/10.7554/eLife.11285.001 PMID:26910315

A Local Forecast of Land Surface Wetness Conditions, Drought, and St. Louis Encephalitis Virus Transmission Derived from Seasonal Climate Predictions

NASA Astrophysics Data System (ADS)

Shaman, J.; Stieglitz, M.; Zebiak, S.; Cane, M.; Day, J. F.

2002-12-01

We present an ensemble local hydrologic forecast derived from the seasonal forecasts of the International Research Institute (IRI) for Climate Prediction. Three- month seasonal forecasts were used to resample historical meteorological conditions and generate ensemble forcing datasets for a TOPMODEL-based hydrology model. Eleven retrospective forecasts were run at a Florida and New York site. Forecast skill was assessed for mean area modeled water table depth (WTD), i.e. near surface soil wetness conditions, and compared with WTD simulated with observed data. Hydrology model forecast skill was evident at the Florida site but not at the New York site. At the Florida site, persistence of hydrologic conditions and local skill of the IRI seasonal forecast contributed to the local hydrologic forecast skill. This forecast will permit probabilistic prediction of future hydrologic conditions. At the Florida site, we have also quantified the link between modeled WTD (i.e. drought) and the amplification and transmission of St. Louis Encephalitis virus (SLEV). We derive an empirical relationship between modeled land surface wetness and levels of SLEV transmission associated with human clinical cases. We then combine the seasonal forecasts of local, modeled WTD with this empirical relationship and produce retrospective probabilistic seasonal forecasts of epidemic SLEV transmission in Florida. Epidemic SLEV transmission forecast skill is demonstrated. These findings will permit real-time forecast of drought and resultant SLEV transmission in Florida.

Parametric decadal climate forecast recalibration (DeFoReSt 1.0)

NASA Astrophysics Data System (ADS)

Pasternack, Alexander; Bhend, Jonas; Liniger, Mark A.; Rust, Henning W.; Müller, Wolfgang A.; Ulbrich, Uwe

2018-01-01

Near-term climate predictions such as decadal climate forecasts are increasingly being used to guide adaptation measures. For near-term probabilistic predictions to be useful, systematic errors of the forecasting systems have to be corrected. While methods for the calibration of probabilistic forecasts are readily available, these have to be adapted to the specifics of decadal climate forecasts including the long time horizon of decadal climate forecasts, lead-time-dependent systematic errors (drift) and the errors in the representation of long-term changes and variability. These features are compounded by small ensemble sizes to describe forecast uncertainty and a relatively short period for which typically pairs of reforecasts and observations are available to estimate calibration parameters. We introduce the Decadal Climate Forecast Recalibration Strategy (DeFoReSt), a parametric approach to recalibrate decadal ensemble forecasts that takes the above specifics into account. DeFoReSt optimizes forecast quality as measured by the continuous ranked probability score (CRPS). Using a toy model to generate synthetic forecast observation pairs, we demonstrate the positive effect on forecast quality in situations with pronounced and limited predictability. Finally, we apply DeFoReSt to decadal surface temperature forecasts from the MiKlip prototype system and find consistent, and sometimes considerable, improvements in forecast quality compared with a simple calibration of the lead-time-dependent systematic errors.

Probabilistic Space Weather Forecasting: a Bayesian Perspective

NASA Astrophysics Data System (ADS)

Camporeale, E.; Chandorkar, M.; Borovsky, J.; Care', A.

2017-12-01

Most of the Space Weather forecasts, both at operational and research level, are not probabilistic in nature. Unfortunately, a prediction that does not provide a confidence level is not very useful in a decision-making scenario. Nowadays, forecast models range from purely data-driven, machine learning algorithms, to physics-based approximation of first-principle equations (and everything that sits in between). Uncertainties pervade all such models, at every level: from the raw data to finite-precision implementation of numerical methods. The most rigorous way of quantifying the propagation of uncertainties is by embracing a Bayesian probabilistic approach. One of the simplest and most robust machine learning technique in the Bayesian framework is Gaussian Process regression and classification. Here, we present the application of Gaussian Processes to the problems of the DST geomagnetic index forecast, the solar wind type classification, and the estimation of diffusion parameters in radiation belt modeling. In each of these very diverse problems, the GP approach rigorously provide forecasts in the form of predictive distributions. In turn, these distributions can be used as input for ensemble simulations in order to quantify the amplification of uncertainties. We show that we have achieved excellent results in all of the standard metrics to evaluate our models, with very modest computational cost.

Long-range hazard assessment of volcanic ash dispersal for a Plinian eruptive scenario at Popocatépetl volcano (Mexico): implications for civil aviation safety

USGS Publications Warehouse

Bonasia, Rosanna; Scaini, Chirara; Capra, Lucia; Nathenson, Manuel; Siebe, Claus; Arana-Salinas, Lilia; Folch, Arnau

2013-01-01

Popocatépetl is one of Mexico’s most active volcanoes threatening a densely populated area that includes Mexico City with more than 20 million inhabitants. The destructive potential of this volcano is demonstrated by its Late Pleistocene–Holocene eruptive activity, which has been characterized by recurrent Plinian eruptions of large magnitude, the last two of which destroyed human settlements in pre-Hispanic times. Popocatépetl’s reawakening in 1994 produced a crisis that culminated with the evacuation of two villages on the northeastern flank of the volcano. Shortly after, a monitoring system and a civil protection contingency plan based on a hazard zone map were implemented. The current volcanic hazards map considers the potential occurrence of different volcanic phenomena, including pyroclastic density currents and lahars. However, no quantitative assessment of the tephra hazard, especially related to atmospheric dispersal, has been performed. The presence of airborne volcanic ash at low and jet-cruise atmospheric levels compromises the safety of aircraft operations and forces re-routing of aircraft to prevent encounters with volcanic ash clouds. Given the high number of important airports in the surroundings of Popocatépetl volcano and considering the potential threat posed to civil aviation in Mexico and adjacent regions in case of a Plinian eruption, a hazard assessment for tephra dispersal is required. In this work, we present the first probabilistic tephra dispersal hazard assessment for Popocatépetl volcano. We compute probabilistic hazard maps for critical thresholds of airborne ash concentrations at different flight levels, corresponding to the situation defined in Europe during 2010, and still under discussion. Tephra dispersal mode is performed using the FALL3D numerical model. Probabilistic hazard maps are built for a Plinian eruptive scenario defined on the basis of geological field data for the “Ochre Pumice” Plinian eruption (4965 14C yr BP). FALL3D model input eruptive parameters are constrained through an inversion method carried out with the semi-analytical HAZMAP model and are varied by sampling them using probability density functions. We analyze the influence of seasonal variations on ash dispersal and estimate the average persistence of critical ash concentrations at relevant locations and airports. This study assesses the impact that a Plinian eruption similar to the Ochre Pumice eruption would have on the main airports of Mexico and adjacent areas. The hazard maps presented here can support long-term planning that would help minimize the impacts of such an eruption on civil aviation.

Long-range hazard assessment of volcanic ash dispersal for a Plinian eruptive scenario at Popocatépetl volcano (Mexico): implications for civil aviation safety

NASA Astrophysics Data System (ADS)

Bonasia, Rosanna; Scaini, Chiara; Capra, Lucia; Nathenson, Manuel; Siebe, Claus; Arana-Salinas, Lilia; Folch, Arnau

2014-01-01

Popocatépetl is one of Mexico's most active volcanoes threatening a densely populated area that includes Mexico City with more than 20 million inhabitants. The destructive potential of this volcano is demonstrated by its Late Pleistocene-Holocene eruptive activity, which has been characterized by recurrent Plinian eruptions of large magnitude, the last two of which destroyed human settlements in pre-Hispanic times. Popocatépetl's reawakening in 1994 produced a crisis that culminated with the evacuation of two villages on the northeastern flank of the volcano. Shortly after, a monitoring system and a civil protection contingency plan based on a hazard zone map were implemented. The current volcanic hazards map considers the potential occurrence of different volcanic phenomena, including pyroclastic density currents and lahars. However, no quantitative assessment of the tephra hazard, especially related to atmospheric dispersal, has been performed. The presence of airborne volcanic ash at low and jet-cruise atmospheric levels compromises the safety of aircraft operations and forces re-routing of aircraft to prevent encounters with volcanic ash clouds. Given the high number of important airports in the surroundings of Popocatépetl volcano and considering the potential threat posed to civil aviation in Mexico and adjacent regions in case of a Plinian eruption, a hazard assessment for tephra dispersal is required. In this work, we present the first probabilistic tephra dispersal hazard assessment for Popocatépetl volcano. We compute probabilistic hazard maps for critical thresholds of airborne ash concentrations at different flight levels, corresponding to the situation defined in Europe during 2010, and still under discussion. Tephra dispersal mode is performed using the FALL3D numerical model. Probabilistic hazard maps are built for a Plinian eruptive scenario defined on the basis of geological field data for the "Ochre Pumice" Plinian eruption (4965 14C yr BP). FALL3D model input eruptive parameters are constrained through an inversion method carried out with the semi-analytical HAZMAP model and are varied by sampling them using probability density functions. We analyze the influence of seasonal variations on ash dispersal and estimate the average persistence of critical ash concentrations at relevant locations and airports. This study assesses the impact that a Plinian eruption similar to the Ochre Pumice eruption would have on the main airports of Mexico and adjacent areas. The hazard maps presented here can support long-term planning that would help minimize the impacts of such an eruption on civil aviation.

Near-term probabilistic forecast of significant wildfire events for the Western United States

Treesearch

Haiganoush K. Preisler; Karin L. Riley; Crystal S. Stonesifer; Dave E. Calkin; Matt Jolly

2016-01-01

Fire danger and potential for large fires in the United States (US) is currently indicated via several forecasted qualitative indices. However, landscape-level quantitative forecasts of the probability of a large fire are currently lacking. In this study, we present a framework for forecasting large fire occurrence - an extreme value event - and evaluating...

The longevity of lava dome eruptions: analysis of the global DomeHaz database

NASA Astrophysics Data System (ADS)

Ogburn, S. E.; Wolpert, R.; Calder, E.; Pallister, J. S.; Wright, H. M. N.

2015-12-01

The likely duration of ongoing volcanic eruptions is a topic of great interest to volcanologists, volcano observatories, and communities near volcanoes. Lava dome forming eruptions can last from days to centuries, and can produce violent, difficult-to-forecast activity including vulcanian to plinian explosions and pyroclastic density currents. Periods of active dome extrusion are often interspersed with periods of relative quiescence, during which extrusion may slow or pause altogether, but persistent volcanic unrest continues. This contribution focuses on the durations of these longer-term unrest phases, hereafter eruptions, that include periods of both lava extrusion and quiescence. A new database of lava dome eruptions, DomeHaz, provides characteristics of 228 eruptions at 127 volcanoes; for which 177 have duration information. We find that while 78% of dome-forming eruptions do not continue for more than 5 years, the remainder can be very long-lived. The probability distributions of eruption durations are shown to be heavy-tailed and vary by magma composition. For this reason, eruption durations are modeled with generalized Pareto distributions whose governing parameters depend on each volcano's composition and eruption duration to date. Bayesian predictive distributions and associated uncertainties are presented for the remaining duration of ongoing eruptions of specified composition and duration to date. Forecasts of such natural events will always have large uncertainties, but the ability to quantify such uncertainty is key to effective communication with stakeholders and to mitigation of hazards. Projections are made for the remaining eruption durations of ongoing eruptions, including those at Soufrière Hills Volcano, Montserrat and Sinabung, Indonesia. This work provides a quantitative, transferable method and rationale on which to base long-term planning decisions for dome forming volcanoes of different compositions, regardless of the quality of an individual volcano's eruptive record, by leveraging a global database.

Satellite-driven modeling approach for monitoring lava flow hazards during the 2017 Etna eruption

NASA Astrophysics Data System (ADS)

Del Negro, C.; Bilotta, G.; Cappello, A.; Ganci, G.; Herault, A.; Zago, V.

2017-12-01

The integration of satellite data and modeling represents an efficient strategy that may provide immediate answers to the main issues raised at the onset of a new effusive eruption. Satellite-based thermal remote sensing of hotspots related to effusive activity can effectively provide a variety of products suited to timing, locating, and tracking the radiant character of lava flows. Hotspots show the location and occurrence of eruptive events (vents). Discharge rate estimates may indicate the current intensity (effusion rate) and potential magnitude (volume). High-spatial resolution multispectral satellite data can complement field observations for monitoring the front position (length) and extension of flows (area). Physics-based models driven, or validated, by satellite-derived parameters are now capable of fast and accurate forecast of lava flow inundation scenarios (hazard). Here, we demonstrate the potential of the integrated application of satellite remote-sensing techniques and lava flow models during the 2017 effusive eruption at Mount Etna in Italy. This combined approach provided insights into lava flow field evolution by supplying detailed views of flow field construction (e.g., the opening of ephemeral vents) that were useful for more accurate and reliable forecasts of eruptive activity. Moreover, we gave a detailed chronology of the lava flow activity based on field observations and satellite images, assessed the potential extent of impacted areas, mapped the evolution of lava flow field, and executed hazard projections. The underside of this combination is the high sensitivity of lava flow inundation scenarios to uncertainties in vent location, discharge rate, and other parameters, which can make interpreting hazard forecasts difficult during an effusive crisis. However, such integration at last makes timely forecasts of lava flow hazards during effusive crises possible at the great majority of volcanoes for which no monitoring exists.

Ensemble superparameterization versus stochastic parameterization: A comparison of model uncertainty representation in tropical weather prediction

NASA Astrophysics Data System (ADS)

Subramanian, Aneesh C.; Palmer, Tim N.

2017-06-01

Stochastic schemes to represent model uncertainty in the European Centre for Medium-Range Weather Forecasts (ECMWF) ensemble prediction system has helped improve its probabilistic forecast skill over the past decade by both improving its reliability and reducing the ensemble mean error. The largest uncertainties in the model arise from the model physics parameterizations. In the tropics, the parameterization of moist convection presents a major challenge for the accurate prediction of weather and climate. Superparameterization is a promising alternative strategy for including the effects of moist convection through explicit turbulent fluxes calculated from a cloud-resolving model (CRM) embedded within a global climate model (GCM). In this paper, we compare the impact of initial random perturbations in embedded CRMs, within the ECMWF ensemble prediction system, with stochastically perturbed physical tendency (SPPT) scheme as a way to represent model uncertainty in medium-range tropical weather forecasts. We especially focus on forecasts of tropical convection and dynamics during MJO events in October-November 2011. These are well-studied events for MJO dynamics as they were also heavily observed during the DYNAMO field campaign. We show that a multiscale ensemble modeling approach helps improve forecasts of certain aspects of tropical convection during the MJO events, while it also tends to deteriorate certain large-scale dynamic fields with respect to stochastically perturbed physical tendencies approach that is used operationally at ECMWF.