Effects of eruption source parameter variation and meteorological dataset on tephra fallout hazard assessment: example from Vesuvius (Italy)

In this study, using the tephra dispersal model HAZMAP, we investigate the effect of using different meteorological datasets and eruption source parameters on tephra fallout hazard assessment for a sub-Plinian eruption of Vesuvius, which is considered as a reference case for hazard assessment analysis. We analyze the effect of using different meteorological data, from: i) radio-sounding carried out at the meteorological station of Brindisi (Italy) between 1962 and 1976 and between 1996 and 2012, and at Pratica di Mare (Rome, Italy) between 1995 and 2013; ii) meteorological models of the National Oceanic and Atmospheric Administration (NOAA), and of the European Centre for Medium-Range Weather Forecasts (ECMWF). Furthermore, we consider the effects of perturbing reference eruptive source parameters. In particular, we vary the total mass, the total grain-size distribution, the column height, and the effective atmospheric diffusion coefficient to evaluate how these parameters affect the hazard probability maps. Moreover, the effect of the seasonal variation of the wind field and the effect of the rain on the deposit loading are considered. Results show that the parameter that mostly affects hazard maps is, as expected, the total erupted mass; furthermore, keeping constant the erupted mass, the most important control on hazard is due to the particle terminal settling velocity distribution which is a function of the total grain-size distribution, particle density and shape. Within the considered range variations, the hazard depends less on the use of different meteorological datasets, column height and effective diffusion coefficient

Grain size distribution uncertainty quantification in volcanic ash dispersal and deposition from weak plumes

We present the results of uncertainty quantification and sensitivity analysis applied to volcanic ash dispersal from weak plumes with focus on the uncertainties associated to the original grain size distribution of the mixture. The Lagrangian particle model Lagrangian Particles Advection Code is used to simulate the transport of inertial particles under the action of realistic atmospheric conditions. The particle motion equations are derived by expressing the particle acceleration as the sum of forces acting along its trajectory, with the drag force calculated as a function of particle diameter, density, shape, and Reynolds number. Simulations are representative of a weak plume event of Mount Etna (Italy) and aimed at quantifying the effect on the dispersal process of the uncertainty in the mean and standard deviation of a lognormal function describing the initial grain size distribution and in particle sphericity. In order to analyze the sensitivity of particle dispersal to these uncertain variables with a reasonable number of simulations, response surfaces in the parameter space are built by using the generalized polynomial chaos expansion technique. The mean diameter and standard deviation of particle size distribution, and their probability density functions, at various distances from the source, both airborne and on ground, are quantified. Results highlight that uncertainty ranges in these quantities are drastically reduced with distance from source, making them largely dependent just on the location. Moreover, at a given distance from source, the distribution is mostly controlled by particle sphericity, particularly on the ground, whereas in air also mean diameter and sorting play a main role

Depressive Symptoms in Fathers during the First Postpartum Year: The Influence of Severity of Preterm Birth, Parenting Stress and Partners' Depression

Although preterm birth constitutes a risk factor for postpartum depressive symptomatology, perinatal depression (PND) has not been investigated extensively in fathers of very low (VLBW) and extremely low birth weight (ELBW) infants. This study explored paternal depression levels at 3, 9, and 12 months of infant corrected age, investigating also the predictive role played by the severity of prematurity, maternal and paternal PND levels, and parenting stress. We recruited 153 fathers of 33 ELBW, 42 VLBW, and 78 full-term (FT) infants, respectively. Depression was investigated by the Edinburgh Postnatal Depression Scale (EPDS) and distress by the Parenting Stress Index-Short Form-PSI-SF (Total and subscales: Parental Distress, Parent-Child Dysfunctional Interaction, and Difficult Child). ELBW fathers showed a significant decrease (improvement) in EPDS, total PSI-SF, and Parental Distress mean scores after 3 months. Paternal EPDS scores at 12 months were significantly predicted by VLBW and FT infants' birth weight categories, fathers' EPDS scores at 3 and 9 months, Parent-Child Dysfunctional Interaction subscale at 3 months, and Difficult Child subscale at 9 months. This study strengthens the relevance of including early routine screening and parenting support for fathers in perinatal health services, with particular attention to fathers who might be more vulnerable to mental health difficulties due to severely preterm birth

Rapid emergency assessment of ash and gas hazard for future eruptions at Santorini Volcano, Greece

Hazard assessments for long-dormant volcanoes, where information is rarely available, typically have to be made rapidly and in the face of considerable uncertainty and often poor information. A conditional (assuming an eruption), scenario-based probabilistic approach to such an assessment is presented here for Santorini volcano (Greece). The rapid assessment was developed and implemented in response to the 2011-2012 unrest crisis in order to inform emergency management and planning. This paper synthesises the results presented to the Greek National Committee and scientific community involved. Two plausible eruptions at Santorini were investigated, using multiple inputs and dispersal models, based on observations of historic eruptions and expert judgement. For ash hazard, a 'most likely' eruption scenario was developed, characterised by slow lava extrusion over periods of one to two years with weak but persistent explosions and ash venting up to 3 km. A second 'largest considered' sub-Plinian explosive scenario assumed a 12 km high column of 4-h duration. For gas hazard, constant fluxes of 200 and 800 tons/day SO2 were assumed for the duration of the eruption scenarios, noting that there is very little evidence to constrain SO2 flux from Santorini eruptions. Statistical models of likely wind conditions with height and season were developed from decadal reanalysis time series showing that consistent low-altitude winds were rarely maintained for more than a few days. Stochastic models of ash (TEPHRA2, VOL-CALPUFF) and gas (AERMOD) dispersal provided outputs in the form of probability maps and exceedance probability curves for key loading and concentration thresholds at important locations on the island. The results from the rapid assessments presented in this paper confirm that ash and gas hazard is likely to be of concern if an eruption of Santorini occurs. Higher hazard may be expected to the south and east of the volcano, notably at important tourist and transport hubs. Low hazard to the north and northwest suggests that these may be suitable locations for emergency response centres and emergency critical infrastructure. This approach may provide a blueprint for rapid ash and gas assessment for other long-dormant volcanoes and we provide suggestions for refining the methods used.</p

Conduit geometry and evolution of effusion rate during basaltic effusive events: Insights from numerical modeling

The dynamics of effusive events is controlled by the interplay between conduit geometry and source conditions. Dyke-like geometries have been commonly employed for describing conduits during effusive eruptions, but their depth-dependent and temporal modifications are largely unknown. Here we present a novel model that describes the evolution of conduit geometry during effusive eruptions by using a quasisteady state approach based on a 1D conduit model and appropriate criteria to model the geometric evolution of the conduit due to fluid shear stress and elastic deformation. Such approach provides time-dependent trends for effusion rate, conduit geometry, exit velocity and gas flow, among other output variables. Fluid shear stress leads to upward widening conduits, whereas elastic deformation becomes relevant only during final phases of the eruptions. Since the model is able to reproduce different trends of effusion rate, it was employed for addressing the effects of magma source conditions and conduit properties on the main characteristics of the resulting effusive eruptions (e.g. duration, erupted mass, maximum effusion rate). We show that the total erupted mass is mainly controlled by magma reservoir dimensions and their conditions before the eruption (i.e., initial overpressure), whereas conduit processes and geometry are able to affect the magma withdrawal rate and thus the eruption duration and effusion rate. The resulting effusion rate trends were classified in different types, and associated to the curves described in the literature for different volcanic events. Results well reproduce these trends and provide new insights for interpreting them, highlighting the importance of reservoir overpressure and the initial dimensions of the feeding dyke on the resulting effusion rate curve

Radial interpolation of GPS and leveling data of ground deformation in a resurgent caldera: application to Campi Flegrei (Italy)

This study presents a new method, called the Radial Interpolation Method, to interpolate data characterized by an approximately radial pattern around a relatively constrained central zone, such as the ground deformation patterns shown in many active volcanic areas. The method enables the fast production of short-term deformation maps on the base of spatially sparse ground deformation measurements and can provide uncertainty quantification on the interpolated values, fundamental for hazard assessment purposes and deformation source reconstruction. The presented approach is not dependent on a priori assumptions about the geometry, location and physical properties of the source, except for the requirement of a locally radial pattern, i.e., allowing multiple centers of symmetry. We test the new method on a synthetic point source example, and then, we apply the method to selected time intervals of real geodetic data collected at the Campi Flegrei caldera during the last 39 years, including examples of leveling, Geodetic Precise Traversing measurements and Global Positioning System. The maps of horizontal displacement, calculated inland, show maximum values lying along a semicircular annular region with a radius of about 2–3 km in size. This semi-annular area is marked by mesoscale structures such as faults, sand dikes and fractures. The maps of vertical displacement describe a linear relation between the maximum vertical uplift measured and the volume variation. The multiplicative factor in the linear relation is about 0.3 × 106 m3/cm if we estimate the proportion of the ΔV that is captured by the GPS network onland and we use this to estimate the full ΔV. In this case, the 95% confidence interval on K because of linear regression is ± 5%. Finally, we briefly discuss how the new method could be used for the production of short-term vent opening maps on the base of real-time geodetic measurements of the horizontal and vertical displacements

The effects of vent location, event scale and time forecasts on pyroclastic density current hazard maps at Campi Flegrei caldera (Italy)

This study presents a new method for producing long-term hazard maps for pyroclastic density currents (PDC) originating at Campi Flegrei caldera. Such method is based on a doubly stochastic approach and is able to combine the uncertainty assessments on the spatial location of the volcanic vent, the size of the flow and the expected time of such an event. The results are obtained by using a Monte Carlo approach and adopting a simplified invasion model based on the box model integral approximation. Temporal assessments are modeled through a Cox-type process including self-excitement effects, based on the eruptive record of the last 15 kyr.Mean and percentilemaps of PDC invasion probability are produced, exploring their sensitivity to some sources of uncertainty and to the effects of the dependence between PDC scales and the caldera sector where they originated. Conditional maps representative of PDC originating inside limited zones of the caldera, or of PDC with a limited range of scales are also produced. Finally, the effect of assuming different time windows for the hazard estimates is explored, also including the potential occurrence of a sequence of multiple events. Assuming that the last eruption of Monte Nuovo (A.D. 1538) marked the beginning of a new epoch of activity similar to the previous ones, results of the statistical analysis indicate a mean probability of PDC invasion above 5% in the next 50 years on almost the entire caldera (with a probability peak of 25% in the central part of the caldera). In contrast, probability values reduce by a factor of about 3 if the entire eruptive record is considered over the last 15 kyr, i.e., including both eruptive epochs and quiescent periods