Ensemble-Based Forecast of Volcanic Clouds Using FALL3D-8.1

Operational forecasting of volcanic ash and SO2 clouds is challenging due to the large uncertainties that typically exist on the eruption source term and the mass removal mechanisms occurring downwind. Current operational forecast systems build on single-run deterministic scenarios that do not account for model input uncertainties and their propagation in time during transport. An ensemble-based forecast strategy has been implemented in the FALL3D-8.1 atmospheric dispersal model to configure, execute, and post-process an arbitrary number of ensemble members in a parallel workflow. In addition to intra-member model domain decomposition, a set of inter-member communicators defines a higher level of code parallelism to enable future incorporation of model data assimilation cycles. Two types of standard products are automatically generated by the ensemble post-process task. On one hand, deterministic forecast products result from some combination of the ensemble members (e.g., ensemble mean, ensemble median, etc.) with an associated quantification of forecast uncertainty given by the ensemble spread. On the other hand, probabilistic products can also be built based on the percentage of members that verify a certain threshold condition. The novel aspect of FALL3D-8.1 is the automatisation of the ensemble-based workflow, including an eventual model validation. To this purpose, novel categorical forecast diagnostic metrics, originally defined in deterministic forecast contexts, are generalised here to probabilistic forecasts in order to have a unique set of skill scores valid to both deterministic and probabilistic forecast contexts. Ensemble-based deterministic and probabilistic approaches are compared using different types of observation datasets (satellite cloud detection and retrieval and deposit thickness observations) for the July 2018 Ambae eruption in the Vanuatu archipelago and the April 2015 Calbuco eruption in Chile. Both ensemble-based approaches outperform single-run simulations in all categorical metrics but no clear conclusion can be extracted on which is the best option between these two.

Uncertainty in soil hydraulic parameters: A review of basics and methods

The paper is a review article on the basics of uncertainty, necessity of its quantification and a comparative study of various methods of uncertainty estimation. The paper primarily focusses on uncertainty estimation of soil hydraulic parameters as of their pivotal importance in groundwater flow and transport simulations, soil moisture modelling techniques etc. The deterministic and probabilistic approaches of uncertainty quantification are studied and an understanding of uncertainty based on field scale measurements, empirical methods and pedotransfer functions is established. A comparative analysis of the basic methods of uncertainty analysis Monte Carlo, Bayesian, FORM/SORM and GLUE is done and the preferential use based on the importance is suggested. Bayesian approach was most suitable for evaluating parametric uncertainty, Monte Carlo was one of the most powerful tools but computationally expensive, FORM was applicable to both numerical and analytical solutions but didn’t guarantee a global convergence and GLUE was conceptually simple but gave only a statistical measure.

Some Remarks on the Use of Deterministic and Probabilistic Approaches in the Evaluation of Rock Slope Stability

The rock slope stability assessment can be performed by means of deterministic and probabilistic approaches. As the deterministic analysis needs only representative values (generally, the mean value) for each physical and geo-mechanical parameter involved, it does not take into account the variability and uncertainty of geo-structural and geo-mechanical properties of joints. This analysis can be usually carried out using different methods, such as the Limit Equilibrium method or numerical modeling techniques sometimes implemented in graphical tests to identify different failure mechanisms (kinematic approach). Probabilistic methods (kinetic approach) aimed to calculate the slope failure probability, consider all orientations, physical characters and shear strength of joints and not only those recognized as kinematically possible. Consequently, the failure probability can be overestimated. It is, therefore, considered more realistic to perform both kinematic and kinetic analyses and to calculate a conditional probability given by the product of the kinematic and kinetic probabilities assuming that they are statistically independent variables. These approaches have been tested on two rock slopes in the Campanian region of Southern Italy affected by possible plane and wedge failures, respectively. Kinematic and kinetic probabilities have been evaluated both by means of the Markland’s test and the Monte Carlo simulation. Using the Eurocode 7, also a deterministic limit equilibrium analysis was performed. The obtained results were compared and commented on.