Articles | Volume 16, issue 12
https://doi.org/10.5194/gmd-16-3459-2023
© Author(s) 2023. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/gmd-16-3459-2023
© Author(s) 2023. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Development and technical paper
| 
27 Jun 2023
Development and technical paper |
| 27 Jun 2023
| 27 Jun 2023
Reconstructing tephra fall deposits via ensemble-based data assimilation techniques
Geosciences Barcelona (GEO3BCN-CSIC), Barcelona, Spain
Barcelona Supercomputing Center, Barcelona, Spain
Antonio Costa
Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Bologna, Bologna, Italy
Giovanni Macedonio
Istituto Nazionale di Geofisica e Vulcanologia, Osservatorio Vesuviano, Naples, Italy
Arnau Folch
Geosciences Barcelona (GEO3BCN-CSIC), Barcelona, Spain
Related authors
Leonardo Mingari, Arnau Folch, Andrew T. Prata, Federica Pardini, Giovanni Macedonio, and Antonio Costa
Atmos. Chem. Phys., 22, 1773–1792, https://doi.org/10.5194/acp-22-1773-2022, https://doi.org/10.5194/acp-22-1773-2022, 2022
Short summary
Short summary
We present a new implementation of an ensemble-based data assimilation method to improve forecasting of volcanic aerosols. This system can be efficiently integrated into operational workflows by exploiting high-performance computing resources. We found a dramatic improvement of forecast quality when satellite retrievals are continuously assimilated. Management of volcanic risk and reduction of aviation impacts can strongly benefit from this research.
Manuel Titos, Beatriz Martínez Montesinos, Sara Barsotti, Laura Sandri, Arnau Folch, Leonardo Mingari, Giovanni Macedonio, and Antonio Costa
Nat. Hazards Earth Syst. Sci., 22, 139–163, https://doi.org/10.5194/nhess-22-139-2022, https://doi.org/10.5194/nhess-22-139-2022, 2022
Short summary
Short summary
This work addresses a quantitative hazard assessment on the possible impact on air traffic of a future ash-forming eruption on the island of Jan Mayen. Through high-performance computing resources, we numerically simulate the transport of ash clouds and ash concentration at different flight levels over an area covering Iceland and the UK using the FALL3D model. This approach allows us to derive a set of probability maps explaining the extent and persisting concentration conditions of ash clouds.
Andrew T. Prata, Leonardo Mingari, Arnau Folch, Giovanni Macedonio, and Antonio Costa
Geosci. Model Dev., 14, 409–436, https://doi.org/10.5194/gmd-14-409-2021, https://doi.org/10.5194/gmd-14-409-2021, 2021
Short summary
Short summary
This paper presents FALL3D-8.0, the latest version release of an open-source code with a track record of 15+ years and a growing number of users in the volcanological and atmospheric communities. The code, originally conceived for atmospheric dispersal and deposition of tephra particles, has been extended to model other types of particles, aerosols and radionuclides. This paper details new model applications and validation of FALL3D-8.0 using satellite, ground-deposit load and radionuclide data.
Arnau Folch, Leonardo Mingari, Natalia Gutierrez, Mauricio Hanzich, Giovanni Macedonio, and Antonio Costa
Geosci. Model Dev., 13, 1431–1458, https://doi.org/10.5194/gmd-13-1431-2020, https://doi.org/10.5194/gmd-13-1431-2020, 2020
Short summary
Short summary
This paper presents FALL3D-8.0, the latest version release of an open-source code with a track record of 15+ years and a growing number of users in the volcanological and atmospheric communities. The code, originally conceived for atmospheric dispersal and deposition of tephra particles, has been extended to model other types of particles, aerosols and radionuclides. This paper details the FALL3D-8.0 model physics and the numerical implementation of the code.
Leonardo A. Mingari, Estela A. Collini, Arnau Folch, Walter Báez, Emilce Bustos, María Soledad Osores, Florencia Reckziegel, Peter Alexander, and José G. Viramonte
Atmos. Chem. Phys., 17, 6759–6778, https://doi.org/10.5194/acp-17-6759-2017, https://doi.org/10.5194/acp-17-6759-2017, 2017
Short summary
Short summary
In this paper, we provide the first comprehensive description of
a dust episode occurred in South America in June 2015 through
observations and numerical simulations. We have investigated
the spatiotemporal distribution of aerosols and the emission
process over complex terrain to gain insight into the key role
played by the orography and the condition that triggered the
long-range transport episode.
A. Folch, L. Mingari, M. S. Osores, and E. Collini
Nat. Hazards Earth Syst. Sci., 14, 119–133, https://doi.org/10.5194/nhess-14-119-2014, https://doi.org/10.5194/nhess-14-119-2014, 2014
Laura Sandri, Mattia de' Michieli Vitturi, Antonio Costa, Mauro Antonio Di Vito, Ilaria Rucco, Domenico Maria Doronzo, Marina Bisson, Roberto Gianardi, Sandro de Vita, and Roberto Sulpizio
Solid Earth, 15, 459–476, https://doi.org/10.5194/se-15-459-2024, https://doi.org/10.5194/se-15-459-2024, 2024
Short summary
Short summary
We study the lahar hazard due to the remobilization of tephra deposits from reference eruptions at Somma–Vesuvius. To this end, we rely on the results of two companion papers dealing with field data and model calibration and run hundreds of simulations from the catchments around the target area to capture the uncertainty in the initial parameters. We process the simulations to draw maps of the probability of overcoming thresholds in lahar flow thickness and dynamic pressure relevant for risk.
Mattia de' Michieli Vitturi, Antonio Costa, Mauro A. Di Vito, Laura Sandri, and Domenico M. Doronzo
Solid Earth, 15, 437–458, https://doi.org/10.5194/se-15-437-2024, https://doi.org/10.5194/se-15-437-2024, 2024
Short summary
Short summary
We present a numerical model for lahars generated by the mobilization of tephra deposits from a reference size eruption at Somma–Vesuvius. The paper presents the model (pyhsics and numerics) and a sensitivity analysis of the processes modelled, numerical schemes, and grid resolution. This work provides the basis for application to hazard quantification for lahars in the Vesuvius area. To this end, we rely on results of the two companion papers (Part 1 on field data, Part 3 on hazard maps).
Mauro Antonio Di Vito, Ilaria Rucco, Sandro de Vita, Domenico Maria Doronzo, Marina Bisson, Mattia de' Michieli Vitturi, Mauro Rosi, Laura Sandri, Giovanni Zanchetta, Elena Zanella, and Antonio Costa
Solid Earth, 15, 405–436, https://doi.org/10.5194/se-15-405-2024, https://doi.org/10.5194/se-15-405-2024, 2024
Short summary
Short summary
We study the distribution of two historical pyroclastic fall–flow and lahar deposits from the sub-Plinian Vesuvius eruptions of 472 CE Pollena and 1631. The motivation comes directly from the widely distributed impact that both the eruptions and lahar phenomena had on the Campanian territory, not only around the volcano but also down the nearby Apennine valleys. Data on about 500 stratigraphic sections and modeling allowed us to evaluate the physical and dynamical impact of these phenomena.
Fabio Dioguardi, Giovanni Chiodini, and Antonio Costa
EGUsphere, https://doi.org/10.5194/egusphere-2023-2867, https://doi.org/10.5194/egusphere-2023-2867, 2023
Short summary
Short summary
We present results of non-volcanic gas (CO2) hazard assessment at the Mefite d’Ansanto area (Italy) where a cold gas stream, which had already been lethal for humans and animals, forms in the valleys surrounding the emission zone. We took the uncertainty related to the gas emission and meteorological conditions into account. Results include maps of CO2 concentration at defined probability levels and of the probability to overcome specified CO2 concentrations over specified time intervals.
Silvia Massaro, Manuel Stocchi, Beatriz Martínez Montesinos, Laura Sandri, Jacopo Selva, Roberto Sulpizio, Biagio Giaccio, Massimiliano Moscatelli, Edoardo Peronace, Marco Nocentini, Roberto Isaia, Manuel Titos Luzón, Pierfrancesco Dellino, Giuseppe Naso, and Antonio Costa
Nat. Hazards Earth Syst. Sci., 23, 2289–2311, https://doi.org/10.5194/nhess-23-2289-2023, https://doi.org/10.5194/nhess-23-2289-2023, 2023
Short summary
Short summary
A new methodology to calculate a probabilistic long-term tephra fallout hazard assessment in southern Italy from the Neapolitan volcanoes is provided. By means of thousands of numerical simulations we quantify the mean annual frequency with which the tephra load at the ground exceeds critical thresholds in 50 years. The output hazard maps account for changes in eruptive regimes of each volcano and are also comparable with those of other natural disasters in which more sources are integrated.
Andrea Bevilacqua, Alvaro Aravena, Willy Aspinall, Antonio Costa, Sue Mahony, Augusto Neri, Stephen Sparks, and Brittain Hill
Nat. Hazards Earth Syst. Sci., 22, 3329–3348, https://doi.org/10.5194/nhess-22-3329-2022, https://doi.org/10.5194/nhess-22-3329-2022, 2022
Short summary
Short summary
We evaluate through first-order kinetic energy models, the minimum volume and mass of a pyroclastic density current generated at the Aso caldera that might affect any of five distal infrastructure sites. These target sites are all located 115–145 km from the caldera, but in well-separated directions. Our constraints of volume and mass are then compared with the scale of Aso-4, the largest caldera-forming eruption of Aso.
Leonardo Mingari, Arnau Folch, Andrew T. Prata, Federica Pardini, Giovanni Macedonio, and Antonio Costa
Atmos. Chem. Phys., 22, 1773–1792, https://doi.org/10.5194/acp-22-1773-2022, https://doi.org/10.5194/acp-22-1773-2022, 2022
Short summary
Short summary
We present a new implementation of an ensemble-based data assimilation method to improve forecasting of volcanic aerosols. This system can be efficiently integrated into operational workflows by exploiting high-performance computing resources. We found a dramatic improvement of forecast quality when satellite retrievals are continuously assimilated. Management of volcanic risk and reduction of aviation impacts can strongly benefit from this research.
Manuel Titos, Beatriz Martínez Montesinos, Sara Barsotti, Laura Sandri, Arnau Folch, Leonardo Mingari, Giovanni Macedonio, and Antonio Costa
Nat. Hazards Earth Syst. Sci., 22, 139–163, https://doi.org/10.5194/nhess-22-139-2022, https://doi.org/10.5194/nhess-22-139-2022, 2022
Short summary
Short summary
This work addresses a quantitative hazard assessment on the possible impact on air traffic of a future ash-forming eruption on the island of Jan Mayen. Through high-performance computing resources, we numerically simulate the transport of ash clouds and ash concentration at different flight levels over an area covering Iceland and the UK using the FALL3D model. This approach allows us to derive a set of probability maps explaining the extent and persisting concentration conditions of ash clouds.
Andrew T. Prata, Leonardo Mingari, Arnau Folch, Giovanni Macedonio, and Antonio Costa
Geosci. Model Dev., 14, 409–436, https://doi.org/10.5194/gmd-14-409-2021, https://doi.org/10.5194/gmd-14-409-2021, 2021
Short summary
Short summary
This paper presents FALL3D-8.0, the latest version release of an open-source code with a track record of 15+ years and a growing number of users in the volcanological and atmospheric communities. The code, originally conceived for atmospheric dispersal and deposition of tephra particles, has been extended to model other types of particles, aerosols and radionuclides. This paper details new model applications and validation of FALL3D-8.0 using satellite, ground-deposit load and radionuclide data.
Silvia Massaro, Roberto Sulpizio, Gianluca Norini, Gianluca Groppelli, Antonio Costa, Lucia Capra, Giacomo Lo Zupone, Michele Porfido, and Andrea Gabrieli
Solid Earth, 11, 2515–2533, https://doi.org/10.5194/se-11-2515-2020, https://doi.org/10.5194/se-11-2515-2020, 2020
Short summary
Short summary
In this work we provide a 2D finite-element modelling of the stress field conditions around the Fuego de Colima volcano (Mexico) in order to test the response of the commercial Linear Static Analysis software to increasingly different geological constraints. Results suggest that an appropriate set of geological and geophysical data improves the mesh generation procedures and the degree of accuracy of numerical outputs, aimed at more reliable physics-based representations of the natural system.
Arnau Folch, Leonardo Mingari, Natalia Gutierrez, Mauricio Hanzich, Giovanni Macedonio, and Antonio Costa
Geosci. Model Dev., 13, 1431–1458, https://doi.org/10.5194/gmd-13-1431-2020, https://doi.org/10.5194/gmd-13-1431-2020, 2020
Short summary
Short summary
This paper presents FALL3D-8.0, the latest version release of an open-source code with a track record of 15+ years and a growing number of users in the volcanological and atmospheric communities. The code, originally conceived for atmospheric dispersal and deposition of tephra particles, has been extended to model other types of particles, aerosols and radionuclides. This paper details the FALL3D-8.0 model physics and the numerical implementation of the code.
Soledad Osores, Juan Ruiz, Arnau Folch, and Estela Collini
Geosci. Model Dev., 13, 1–22, https://doi.org/10.5194/gmd-13-1-2020, https://doi.org/10.5194/gmd-13-1-2020, 2020
Short summary
Short summary
Volcanic ash dispersal forecasts are routinely used to avoid aircraft encounters with volcanic ash. However, the accuracy of these forecasts depends on the knowledge of key factors that are usually difficult to observe directly. In this work we apply an inverse methodology to improve ash concentration forecasts. Results are encouraging, showing that accurate estimations of ash emissions can be performed using the proposed approach, leading to an improvement in ash concentration forecasts.
Silvia Massaro, Antonio Costa, Roberto Sulpizio, Diego Coppola, and Lucia Capra
Solid Earth, 10, 1429–1450, https://doi.org/10.5194/se-10-1429-2019, https://doi.org/10.5194/se-10-1429-2019, 2019
Short summary
Short summary
The Fuego de Colima volcano (Mexico) shows a complex eruptive history, with periods of rapid and slow lava dome growth punctuated by explosive activity. Here we reconstructed the 1998–2018 average discharge rate by means of satellite thermal data and the literature. Using spectral and wavelet analysis, we found a multi-term cyclic behavior that is in good agreement with numerical modeling, accounting for a variable magmatic feeding system composed of a single or double magma chamber system.
Matthieu Poret, Stefano Corradini, Luca Merucci, Antonio Costa, Daniele Andronico, Mario Montopoli, Gianfranco Vulpiani, and Valentin Freret-Lorgeril
Atmos. Chem. Phys., 18, 4695–4714, https://doi.org/10.5194/acp-18-4695-2018, https://doi.org/10.5194/acp-18-4695-2018, 2018
Short summary
Short summary
This study aims at proposing a method to better assess the initial magma fragmentation produced during explosive volcanic eruptions. We worked on merging field, radar, and satellite data to estimate the total grain-size distribution, which is used within simulations to reconstruct the tephra loading and far-travelling airborne ash dispersal. This approach is applied to 23 November 2013, giving the very fine ash fraction related to volcanic hazards (e.g. air traffic safety).
Alejandro Marti and Arnau Folch
Atmos. Chem. Phys., 18, 4019–4038, https://doi.org/10.5194/acp-18-4019-2018, https://doi.org/10.5194/acp-18-4019-2018, 2018
Short summary
Short summary
We use the NMMB-MONARCH-ASH model to quantify the systematic errors associated with traditional offline modeling systems used for operational volcanic ash forecast. Evaluation scores indicate that uncertainties credited to offline modeling are of the same order of magnitude as those associated with the source term, failing to reproduce up to 45–70 % of the ash cloud of an online forecast. This work encourages operational groups to consider online dispersal models for real-time aviation advisory.
Arnau Folch, Jordi Barcons, Tomofumi Kozono, and Antonio Costa
Nat. Hazards Earth Syst. Sci., 17, 861–879, https://doi.org/10.5194/nhess-17-861-2017, https://doi.org/10.5194/nhess-17-861-2017, 2017
Short summary
Short summary
Atmospheric dispersal of a gas denser than air can threat the environment and surrounding communities. In complex terrains, microscale winds and local orographic features can have a strong influence on the gas cloud behavior, potentially leading to inaccurate model results if not captured by coarser-scale simulations. We introduce a methodology for microscale wind field characterization and validate it using, as a test case, the CO2 gas dispersal from 1986 Lake Nyos eruption.
Leonardo A. Mingari, Estela A. Collini, Arnau Folch, Walter Báez, Emilce Bustos, María Soledad Osores, Florencia Reckziegel, Peter Alexander, and José G. Viramonte
Atmos. Chem. Phys., 17, 6759–6778, https://doi.org/10.5194/acp-17-6759-2017, https://doi.org/10.5194/acp-17-6759-2017, 2017
Short summary
Short summary
In this paper, we provide the first comprehensive description of
a dust episode occurred in South America in June 2015 through
observations and numerical simulations. We have investigated
the spatiotemporal distribution of aerosols and the emission
process over complex terrain to gain insight into the key role
played by the orography and the condition that triggered the
long-range transport episode.
Alejandro Marti, Arnau Folch, Oriol Jorba, and Zavisa Janjic
Atmos. Chem. Phys., 17, 4005–4030, https://doi.org/10.5194/acp-17-4005-2017, https://doi.org/10.5194/acp-17-4005-2017, 2017
Short summary
Short summary
We describe and evaluate NMMB-MONARCH-ASH, a novel online multi-scale meteorological and transport model developed at the BSC-CNS capable of forecasting the dispersal and deposition of volcanic ash. The forecast skills of the model have been validated and they improve on those from traditional operational offline (decoupled) models. The results support the use of online coupled models to aid civil aviation and emergency management during a crisis such as the 2010 eruption of Eyjafjallajökull.
A. Folch, A. Costa, and G. Macedonio
Geosci. Model Dev., 9, 431–450, https://doi.org/10.5194/gmd-9-431-2016, https://doi.org/10.5194/gmd-9-431-2016, 2016
Short summary
Short summary
We present FPLUME-1.0, a steady-state 1-D cross-section-averaged eruption column model based on the buoyant plume theory (BPT). The model accounts for plume bending by wind, entrainment of ambient moisture, effects of water phase changes, particle fallout and re-entrainment, a new parameterization for the air entrainment coefficients and a model for wet aggregation of ash particles in presence of liquid water or ice.
R. Tonini, L. Sandri, A. Costa, and J. Selva
Nat. Hazards Earth Syst. Sci., 15, 409–415, https://doi.org/10.5194/nhess-15-409-2015, https://doi.org/10.5194/nhess-15-409-2015, 2015
S. Biass, C. Scaini, C. Bonadonna, A. Folch, K. Smith, and A. Höskuldsson
Nat. Hazards Earth Syst. Sci., 14, 2265–2287, https://doi.org/10.5194/nhess-14-2265-2014, https://doi.org/10.5194/nhess-14-2265-2014, 2014
C. Scaini, S. Biass, A. Galderisi, C. Bonadonna, A. Folch, K. Smith, and A. Höskuldsson
Nat. Hazards Earth Syst. Sci., 14, 2289–2312, https://doi.org/10.5194/nhess-14-2289-2014, https://doi.org/10.5194/nhess-14-2289-2014, 2014
A. Folch, L. Mingari, M. S. Osores, and E. Collini
Nat. Hazards Earth Syst. Sci., 14, 119–133, https://doi.org/10.5194/nhess-14-119-2014, https://doi.org/10.5194/nhess-14-119-2014, 2014
Related subject area
Solid Earth
ShellSet v1.1.0 parallel dynamic neotectonic modelling: a case study using Earth5-049
FastIsostasy v1.0 – a regional, accelerated 2D glacial isostatic adjustment (GIA) model accounting for the lateral variability of the solid Earth
Automatic adjoint-based inversion schemes for geodynamics: reconstructing the evolution of Earth's mantle in space and time
Reconciling Surface Deflections From Simulations of Global Mantle Convection
Benchmarking the accuracy of higher-order particle methods in geodynamic models of transient flow
REHEATFUNQ (REgional HEAT-Flow Uncertainty and aNomaly Quantification) 2.0.1: a model for regional aggregate heat flow distributions and anomaly quantification
A new temperature–photoperiod coupled phenology module in LPJ-GUESS model v4.1: optimizing estimation of terrestrial carbon and water processes
High-precision 1′ × 1′ bathymetric model of Philippine Sea inversed from marine gravity anomalies
Deciphering past earthquakes from the probabilistic modeling of paleoseismic records – the Paleoseismic EArthquake CHronologies code (PEACH, version 1)
Modelling detrital cosmogenic nuclide concentrations during landscape evolution in Cidre v2.0
IMEX_SfloW2D v2: a depth-averaged numerical flow model for volcanic gas–particle flows over complex topographies and water
A Fast Surrogate Model for 3D-Earth Glacial Isostatic Adjustment using Tensorflow (v2.8.10) Artificial Neural Networks
Simulation of a fully coupled 3D glacial isostatic adjustment – ice sheet model for the Antarctic ice sheet over a glacial cycle
Three-Dimensional Analytical Solution of Self-potential from Regularly Polarized Bodies in Layered Seafloor Model
AdaHRBF v1.0: gradient-adaptive Hermite–Birkhoff radial basis function interpolants for three-dimensional stratigraphic implicit modeling
PySubdiv 1.0: open-source geological modeling and reconstruction by non-manifold subdivision surfaces
ClinoformNet-1.0: stratigraphic forward modeling and deep learning for seismic clinoform delineation
Addressing challenges in uncertainty quantification: the case of geohazard assessments
DeepISMNet: three-dimensional implicit structural modeling with convolutional neural network
Towards automatic finite-element methods for geodynamics via Firedrake
MagmaFOAM-1.0: a modular framework for the simulation of magmatic systems
A global, spherical finite-element model for post-seismic deformation using Abaqus
SMAUG v1.0 – a user-friendly muon simulator for the imaging of geological objects in 3-D
CliffDelineaTool v1.2.0: an algorithm for identifying coastal cliff base and top positions
Capturing the interactions between ice sheets, sea level and the solid Earth on a range of timescales: a new “time window” algorithm
Structural, petrophysical, and geological constraints in potential field inversion using the Tomofast-x v1.0 open-source code
Spatial agents for geological surface modelling
RHEA v1.0: Enabling fully coupled simulations with hydro-geomechanical heterogeneity
Modelling of faults in LoopStructural 1.0
PALEOSTRIPv1.0 – a user-friendly 3D backtracking software to reconstruct paleo-bathymetries
LoopStructural 1.0: time-aware geological modelling
Sub3DNet1.0: a deep-learning model for regional-scale 3D subsurface structure mapping
Analytical solutions for mantle flow in cylindrical and spherical shells
Towards a model for structured mass movements: the OpenLISEM hazard model 2.0a
GO_3D_OBS: the multi-parameter benchmark geomodel for seismic imaging method assessment and next-generation 3D survey design (version 1.0)
PLUME-MoM-TSM 1.0.0: a volcanic column and umbrella cloud spreading model
HydrothermalFoam v1.0: a 3-D hydrothermal transport model for natural submarine hydrothermal systems
Synthetic seismicity distribution in Guerrero–Oaxaca subduction zone, Mexico, and its implications on the role of asperities in Gutenberg–Richter law
A new open-source viscoelastic solid earth deformation module implemented in Elmer (v8.4)
CobWeb 1.0: machine learning toolbox for tomographic imaging
pygeodyn 1.1.0: a Python package for geomagnetic data assimilation
IMEX_SfloW2D 1.0: a depth-averaged numerical flow model for pyroclastic avalanches
A multilayer approach and its application to model a local gravimetric quasi-geoid model over the North Sea: QGNSea V1.0
Development of an automatic delineation of cliff top and toe on very irregular planform coastlines (CliffMetrics v1.0)
Bayesian inference of earthquake rupture models using polynomial chaos expansion
Geodynamic diagnostics, scientific visualisation and StagLab 3.0
SaLEM (v1.0) – the Soil and Landscape Evolution Model (SaLEM) for simulation of regolith depth in periglacial environments
SILLi 1.0: a 1-D numerical tool quantifying the thermal effects of sill intrusions
The SPACE 1.0 model: a Landlab component for 2-D calculation of sediment transport, bedrock erosion, and landscape evolution
Ellipsoids (v1.0): 3-D magnetic modelling of ellipsoidal bodies
Jon B. May, Peter Bird, and Michele M. C. Carafa
Geosci. Model Dev., 17, 6153–6171, https://doi.org/10.5194/gmd-17-6153-2024, https://doi.org/10.5194/gmd-17-6153-2024, 2024
Short summary
Short summary
ShellSet is a combination of well-known geoscience software packages. It features a simple user interface and is optimised through the addition of a grid search input option (automatically searching for optimal models within a defined N-dimensional parameter space) and the ability to run multiple models in parallel. We show that for each number of models tested there is a performance benefit to parallel running, while two examples demonstrate a use case by improving an existing global model.
Jan Swierczek-Jereczek, Marisa Montoya, Konstantin Latychev, Alexander Robinson, Jorge Alvarez-Solas, and Jerry Mitrovica
Geosci. Model Dev., 17, 5263–5290, https://doi.org/10.5194/gmd-17-5263-2024, https://doi.org/10.5194/gmd-17-5263-2024, 2024
Short summary
Short summary
Ice sheets present a thickness of a few kilometres, leading to a vertical deformation of the crust of up to a kilometre. This process depends on properties of the solid Earth, which can be regionally very different. We propose a model that accounts for this often-ignored heterogeneity and run 100 000 simulation years in minutes. Thus, the evolution of ice sheets is modeled with better accuracy, which is critical for a good mitigation of climate change and, in particular, sea-level rise.
Sia Ghelichkhan, Angus Gibson, D. Rhodri Davies, Stephan C. Kramer, and David A. Ham
Geosci. Model Dev., 17, 5057–5086, https://doi.org/10.5194/gmd-17-5057-2024, https://doi.org/10.5194/gmd-17-5057-2024, 2024
Short summary
Short summary
We introduce the Geoscientific ADjoint Optimisation PlaTform (G-ADOPT), designed for inverse modelling of Earth system processes, with an initial focus on mantle dynamics. G-ADOPT is built upon Firedrake, Dolfin-Adjoint and the Rapid Optimisation Library, which work together to optimise models using an adjoint method, aligning them with seismic and geologic datasets. We demonstrate G-ADOPT's ability to reconstruct mantle evolution and thus be a powerful tool in geosciences.
Conor P. B. O'Malley, Gareth G. Roberts, James Panton, Fred D. Richards, J. Huw Davies, Victoria M. Fernandes, and Sia Ghelichkhan
EGUsphere, https://doi.org/10.5194/egusphere-2024-1893, https://doi.org/10.5194/egusphere-2024-1893, 2024
Short summary
Short summary
We wish to understand how the history of flowing rock within Earth's interior impacts deflection of its surface. Observations exist to address this problem, and mathematics and different computing tools can be used to predict histories of flow. We explore how modelling choices impact calculated vertical deflections. The sensitivity of vertical motions at Earth's surface to deep flow is assessed, demonstrating how surface observations can enlighten flow histories.
Rene Gassmöller, Juliane Dannberg, Wolfgang Bangerth, Elbridge Gerry Puckett, and Cedric Thieulot
Geosci. Model Dev., 17, 4115–4134, https://doi.org/10.5194/gmd-17-4115-2024, https://doi.org/10.5194/gmd-17-4115-2024, 2024
Short summary
Short summary
Numerical models that use simulated particles are a powerful tool for investigating flow in the interior of the Earth, but the accuracy of these models is not fully understood. Here we present two new benchmarks that allow measurement of model accuracy. We then document that better accuracy matters for applications like convection beneath an oceanic plate. Our benchmarks and methods are freely available to help the community develop better models.
Malte Jörn Ziebarth and Sebastian von Specht
Geosci. Model Dev., 17, 2783–2828, https://doi.org/10.5194/gmd-17-2783-2024, https://doi.org/10.5194/gmd-17-2783-2024, 2024
Short summary
Short summary
Thermal energy from Earth’s active interior constantly dissipates through Earth’s surface. This heat flow is not spatially uniform, and its exact pattern is hard to predict since it depends on crustal and mantle properties, both varying across scales. Our new model REHEATFUNQ addresses this difficulty by treating the fluctuations of heat flow within a region statistically. REHEATFUNQ estimates the regional distribution of heat flow and quantifies known structural signals therein.
Shouzhi Chen, Yongshuo H. Fu, Mingwei Li, Zitong Jia, Yishuo Cui, and Jing Tang
Geosci. Model Dev., 17, 2509–2523, https://doi.org/10.5194/gmd-17-2509-2024, https://doi.org/10.5194/gmd-17-2509-2024, 2024
Short summary
Short summary
It is still a challenge to achieve an accurate simulation of vegetation phenology in the dynamic global vegetation models (DGVMs). We implemented and coupled the spring and autumn phenology models into one of the DGVMs, LPJ-GUESS, and substantially improved the accuracy in capturing the start and end dates of growing seasons. Our study highlights the importance of getting accurate phenology estimations to reduce the uncertainties in plant distribution and terrestrial carbon and water cycling.
Dechao An, Jinyun Guo, Xiaotao Chang, Zhenming Wang, Yongjun Jia, Xin Liu, Valery Bondur, and Heping Sun
Geosci. Model Dev., 17, 2039–2052, https://doi.org/10.5194/gmd-17-2039-2024, https://doi.org/10.5194/gmd-17-2039-2024, 2024
Short summary
Short summary
Seafloor topography, as fundamental geoinformation in marine surveying and mapping, plays a crucial role in numerous scientific studies. In this paper, we focus on constructing a high-precision seafloor topography and bathymetry model for the Philippine Sea (5° N–35° N, 120° E–150° E), based on shipborne bathymetric data and marine gravity anomalies, and evaluate the reliability of the model's accuracy.
Octavi Gómez-Novell, Bruno Pace, Francesco Visini, Joanna Faure Walker, and Oona Scotti
Geosci. Model Dev., 16, 7339–7355, https://doi.org/10.5194/gmd-16-7339-2023, https://doi.org/10.5194/gmd-16-7339-2023, 2023
Short summary
Short summary
Knowing the rate at which earthquakes happen along active faults is crucial to characterize the hazard that they pose. We present an approach (Paleoseismic EArthquake CHronologies, PEACH) to correlate and compute seismic histories using paleoseismic data, a type of data that characterizes past seismic activity from the geological record. Our approach reduces the uncertainties of the seismic histories and overall can improve the knowledge on fault rupture behavior for the seismic hazard.
Sébastien Carretier, Vincent Regard, Youssouf Abdelhafiz, and Bastien Plazolles
Geosci. Model Dev., 16, 6741–6755, https://doi.org/10.5194/gmd-16-6741-2023, https://doi.org/10.5194/gmd-16-6741-2023, 2023
Short summary
Short summary
We present the development of a code to simulate simultaneously the dynamics of landscapes over geological time and the evolution of the concentration of cosmogenic isotopes in grains throughout their transport from the slopes to the river outlets. This new model makes it possible to study the relationship between the detrital signal of cosmogenic isotope concentration measured in sediment and the erosion--deposition processes in watersheds.
Mattia de' Michieli Vitturi, Tomaso Esposti Ongaro, and Samantha Engwell
Geosci. Model Dev., 16, 6309–6336, https://doi.org/10.5194/gmd-16-6309-2023, https://doi.org/10.5194/gmd-16-6309-2023, 2023
Short summary
Short summary
We present version 2 of the numerical code IMEX-SfloW2D. With this version it is possible to simulate a wide range of volcanic mass flows (pyroclastic avalanches, lahars, pyroclastic surges), and here we present its application to transient dilute pyroclastic density currents (PDCs). A simulation of the 1883 Krakatau eruption demonstrates the capability of the numerical model to face a complex natural case involving the propagation of PDCs over the sea surface and across topographic obstacles.
Ryan Love, Glenn A. Milne, Parviz Ajourlou, Soran Parang, Lev Tarasov, and Konstantin Latychev
EGUsphere, https://doi.org/10.5194/egusphere-2023-2491, https://doi.org/10.5194/egusphere-2023-2491, 2023
Short summary
Short summary
A relatively recent advance in glacial isostatic adjustment modelling has been the development of models that include 3D Earth structure, as opposed to 1D structure. However, a major limitation is the computational expense. We have developed a method using artificial neural networks to emulate the influence of 3D Earth models to affordably constrain the viscosity parameter space. Our results indicate that the misfits are of a scale such that useful predictions of relative sea level can be made.
Caroline J. van Calcar, Roderik S. W. van de Wal, Bas Blank, Bas de Boer, and Wouter van der Wal
Geosci. Model Dev., 16, 5473–5492, https://doi.org/10.5194/gmd-16-5473-2023, https://doi.org/10.5194/gmd-16-5473-2023, 2023
Short summary
Short summary
The waxing and waning of the Antarctic ice sheet caused the Earth’s surface to deform, which is stabilizing the ice sheet and mainly determined by the spatially variable viscosity of the mantle. Including this feedback in model simulations led to significant differences in ice sheet extent and ice thickness over the last glacial cycle. The results underline and quantify the importance of including this local feedback effect in ice sheet models when simulating the Antarctic ice sheet evolution.
Pengfei Zhang, Yi-an Cui, Jing Xie, Youjun Guo, Jianxin Liu, and Jieran Liu
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2023-94, https://doi.org/10.5194/gmd-2023-94, 2023
Revised manuscript accepted for GMD
Short summary
Short summary
A reasonable self-potential (SP) forward modeling is fundamental for mineral exploration. In this paper, we present a method to obtain the theoretical solution of SP generated by regularly polarized bodies in layered media. The results demonstrate that the measured SP data is consistent with the analytical solution, validating the proposed method and corresponding analytical solution.
Baoyi Zhang, Linze Du, Umair Khan, Yongqiang Tong, Lifang Wang, and Hao Deng
Geosci. Model Dev., 16, 3651–3674, https://doi.org/10.5194/gmd-16-3651-2023, https://doi.org/10.5194/gmd-16-3651-2023, 2023
Short summary
Short summary
We propose a Hermite–Birkhoff radial basis function (HRBF) formulation, AdaHRBF, with an adaptive gradient magnitude for continuous 3D stratigraphic potential field (SPF) modeling of multiple stratigraphic interfaces. In the linear system of HRBF interpolants constrained by the scattered on-contact attribute points and off-contact attitude points of a set of strata in 3D space, we add a novel optimization term to iteratively obtain the true gradient magnitude.
Mohammad Moulaeifard, Simon Bernard, and Florian Wellmann
Geosci. Model Dev., 16, 3565–3579, https://doi.org/10.5194/gmd-16-3565-2023, https://doi.org/10.5194/gmd-16-3565-2023, 2023
Short summary
Short summary
In this work, we propose a flexible framework to generate and interact with geological models using explicit surface representations. The essence of the work lies in the determination of the flexible control mesh, topologically similar to the main geological structure, watertight and controllable with few control points, to manage the geological structures. We exploited the subdivision surface method in our work, which is commonly used in the animation and gaming industry.
Hui Gao, Xinming Wu, Jinyu Zhang, Xiaoming Sun, and Zhengfa Bi
Geosci. Model Dev., 16, 2495–2513, https://doi.org/10.5194/gmd-16-2495-2023, https://doi.org/10.5194/gmd-16-2495-2023, 2023
Short summary
Short summary
We propose a workflow to automatically generate synthetic seismic data and corresponding stratigraphic labels (e.g., clinoform facies, relative geologic time, and synchronous horizons) by geological and geophysical forward modeling. Trained with only synthetic datasets, our network works well to accurately and efficiently predict clinoform facies in 2D and 3D field seismic data. Such a workflow can be easily extended for other geological and geophysical scenarios in the future.
Ibsen Chivata Cardenas, Terje Aven, and Roger Flage
Geosci. Model Dev., 16, 1601–1615, https://doi.org/10.5194/gmd-16-1601-2023, https://doi.org/10.5194/gmd-16-1601-2023, 2023
Short summary
Short summary
We discuss challenges in uncertainty quantification for geohazard assessments. The challenges arise from limited data and the one-off nature of geohazard features. The challenges include the credibility of predictions, input uncertainty, and assumptions’ impact. Considerations to increase credibility of the quantification are provided. Crucial tasks in the quantification are the exhaustive scrutiny of the background knowledge coupled with the assessment of deviations of assumptions made.
Zhengfa Bi, Xinming Wu, Zhaoliang Li, Dekuan Chang, and Xueshan Yong
Geosci. Model Dev., 15, 6841–6861, https://doi.org/10.5194/gmd-15-6841-2022, https://doi.org/10.5194/gmd-15-6841-2022, 2022
Short summary
Short summary
We present an implicit modeling method based on deep learning to produce a geologically valid and structurally compatible model from unevenly sampled structural data. Trained with automatically generated synthetic data with realistic features, our network can efficiently model geological structures without the need to solve large systems of mathematical equations, opening new opportunities for further leveraging deep learning to improve modeling capacity in many Earth science applications.
D. Rhodri Davies, Stephan C. Kramer, Sia Ghelichkhan, and Angus Gibson
Geosci. Model Dev., 15, 5127–5166, https://doi.org/10.5194/gmd-15-5127-2022, https://doi.org/10.5194/gmd-15-5127-2022, 2022
Short summary
Short summary
Firedrake is a state-of-the-art system that automatically generates highly optimised code for simulating finite-element (FE) problems in geophysical fluid dynamics. It creates a separation of concerns between employing the FE method and implementing it. Here, we demonstrate the applicability and benefits of Firedrake for simulating geodynamical flows, with a focus on the slow creeping motion of Earth's mantle over geological timescales, which is ultimately the engine driving our dynamic Earth.
Federico Brogi, Simone Colucci, Jacopo Matrone, Chiara Paola Montagna, Mattia De' Michieli Vitturi, and Paolo Papale
Geosci. Model Dev., 15, 3773–3796, https://doi.org/10.5194/gmd-15-3773-2022, https://doi.org/10.5194/gmd-15-3773-2022, 2022
Short summary
Short summary
Computer simulations play a fundamental role in understanding volcanic phenomena. The growing complexity of these simulations requires the development of flexible computational tools that can easily switch between sub-models and solution techniques as well as optimizations. MagmaFOAM is a newly developed library that allows for maximum flexibility for solving multiphase volcanic flows and promotes collaborative work for in-house and community model development, testing, and comparison.
Grace A. Nield, Matt A. King, Rebekka Steffen, and Bas Blank
Geosci. Model Dev., 15, 2489–2503, https://doi.org/10.5194/gmd-15-2489-2022, https://doi.org/10.5194/gmd-15-2489-2022, 2022
Short summary
Short summary
We present a finite-element model of post-seismic solid Earth deformation built in the software package Abaqus for the purpose of calculating post-seismic deformation in the far field of major earthquakes. The model is benchmarked against an existing open-source post-seismic model demonstrating good agreement. The advantage over existing models is the potential for simple modification to include 3-D Earth structure, non-linear rheologies and alternative or multiple sources of stress change.
Alessandro Lechmann, David Mair, Akitaka Ariga, Tomoko Ariga, Antonio Ereditato, Ryuichi Nishiyama, Ciro Pistillo, Paola Scampoli, Mykhailo Vladymyrov, and Fritz Schlunegger
Geosci. Model Dev., 15, 2441–2473, https://doi.org/10.5194/gmd-15-2441-2022, https://doi.org/10.5194/gmd-15-2441-2022, 2022
Short summary
Short summary
Muon tomography is a technology that is used often in geoscientific research. The know-how of data analysis is, however, still possessed by physicists who developed this technology. This article aims at providing geoscientists with the necessary tools to perform their own analyses. We hope that a lower threshold to enter the field of muon tomography will allow more geoscientists to engage with muon tomography. SMAUG is set up in a modular way to allow for its own modules to work in between.
Zuzanna M. Swirad and Adam P. Young
Geosci. Model Dev., 15, 1499–1512, https://doi.org/10.5194/gmd-15-1499-2022, https://doi.org/10.5194/gmd-15-1499-2022, 2022
Short summary
Short summary
Cliff base and top lines that delimit coastal cliff faces are usually manually digitized based on maps, aerial photographs, terrain models, etc. However, manual mapping is time consuming and depends on the mapper's decisions and skills. To increase the objectivity and efficiency of cliff mapping, we developed CliffDelineaTool, an algorithm that identifies cliff base and top positions along cross-shore transects using elevation and slope characteristics.
Holly Kyeore Han, Natalya Gomez, and Jeannette Xiu Wen Wan
Geosci. Model Dev., 15, 1355–1373, https://doi.org/10.5194/gmd-15-1355-2022, https://doi.org/10.5194/gmd-15-1355-2022, 2022
Short summary
Short summary
Interactions between ice sheets, sea level and the solid Earth occur over a range of timescales from years to tens of thousands of years. This requires coupled ice-sheet–sea-level models to exchange information frequently, leading to a quadratic increase in computation time with the number of model timesteps. We present a new sea-level model algorithm that allows coupled models to improve the computational feasibility and precisely capture short-term interactions within longer simulations.
Jérémie Giraud, Vitaliy Ogarko, Roland Martin, Mark Jessell, and Mark Lindsay
Geosci. Model Dev., 14, 6681–6709, https://doi.org/10.5194/gmd-14-6681-2021, https://doi.org/10.5194/gmd-14-6681-2021, 2021
Short summary
Short summary
We review different techniques to model the Earth's subsurface from geophysical data (gravity field anomaly, magnetic field anomaly) using geological models and measurements of the rocks' properties. We show examples of application using idealised examples reproducing realistic features and provide theoretical details of the open-source algorithm we use.
Eric A. de Kemp
Geosci. Model Dev., 14, 6661–6680, https://doi.org/10.5194/gmd-14-6661-2021, https://doi.org/10.5194/gmd-14-6661-2021, 2021
Short summary
Short summary
This is a proof of concept and review paper of spatial agents, with initial research focusing on geomodelling. The results may be of interest to others working on complex regional geological modelling with sparse data. Structural agent-based swarming behaviour is key to advancing this field. The study provides groundwork for research in structural geology 3D modelling with spatial agents. This work was done with NetLogo, a free agent modelling platform used mostly for teaching complex systems.
José M. Bastías Espejo, Andy Wilkins, Gabriel C. Rau, and Philipp Blum
Geosci. Model Dev., 14, 6257–6272, https://doi.org/10.5194/gmd-14-6257-2021, https://doi.org/10.5194/gmd-14-6257-2021, 2021
Short summary
Short summary
The hydraulic and mechanical properties of the subsurface are inherently heterogeneous. RHEA is a simulator that can perform couple hydro-geomechanical processes in heterogeneous porous media with steep gradients. RHEA is able to fully integrate spatial heterogeneity, allowing allocation of distributed hydraulic and geomechanical properties at mesh element level. RHEA is a valuable tool that can simulate problems considering realistic heterogeneity inherent to geologic formations.
Lachlan Grose, Laurent Ailleres, Gautier Laurent, Guillaume Caumon, Mark Jessell, and Robin Armit
Geosci. Model Dev., 14, 6197–6213, https://doi.org/10.5194/gmd-14-6197-2021, https://doi.org/10.5194/gmd-14-6197-2021, 2021
Short summary
Short summary
Fault discontinuities in rock packages represent the plane where two blocks of rock have moved. They are challenging to incorporate into geological models because the geometry of the faulted rock units are defined by not only the location of the discontinuity but also the kinematics of the fault. In this paper, we outline a structural geology framework for incorporating faults into geological models by directly incorporating kinematics into the mathematical framework of the model.
Florence Colleoni, Laura De Santis, Enrico Pochini, Edy Forlin, Riccardo Geletti, Giuseppe Brancatelli, Magdala Tesauro, Martina Busetti, and Carla Braitenberg
Geosci. Model Dev., 14, 5285–5305, https://doi.org/10.5194/gmd-14-5285-2021, https://doi.org/10.5194/gmd-14-5285-2021, 2021
Short summary
Short summary
PALEOSTRIP has been developed in the framework of past Antarctic ice sheet reconstructions for periods when bathymetry around Antarctica differed substantially from today. It has been designed for users with no knowledge of numerical modelling and allows users to switch on and off the processes involved in backtracking and backstripping. Applications are broad, and it can be used to restore any continental margin bathymetry or sediment thickness and to perform basin analysis.
Lachlan Grose, Laurent Ailleres, Gautier Laurent, and Mark Jessell
Geosci. Model Dev., 14, 3915–3937, https://doi.org/10.5194/gmd-14-3915-2021, https://doi.org/10.5194/gmd-14-3915-2021, 2021
Short summary
Short summary
LoopStructural is an open-source 3D geological modelling library with a model design allowing for multiple different algorithms to be used for comparison for the same geology. Geological structures are modelled using structural geology concepts and techniques, allowing for complex structures such as overprinted folds and faults to be modelled. In the paper, we demonstrate automatically generating a 3-D model from map2loop-processed geological survey data of the Flinders Ranges, South Australia.
Zhenjiao Jiang, Dirk Mallants, Lei Gao, Tim Munday, Gregoire Mariethoz, and Luk Peeters
Geosci. Model Dev., 14, 3421–3435, https://doi.org/10.5194/gmd-14-3421-2021, https://doi.org/10.5194/gmd-14-3421-2021, 2021
Short summary
Short summary
Fast and reliable tools are required to extract hidden information from big geophysical and remote sensing data. A deep-learning model in 3D image construction from 2D image(s) is here developed for paleovalley mapping from globally available digital elevation data. The outstanding performance for 3D subsurface imaging gives confidence that this generic novel tool will make better use of existing geophysical and remote sensing data for improved management of limited earth resources.
Stephan C. Kramer, D. Rhodri Davies, and Cian R. Wilson
Geosci. Model Dev., 14, 1899–1919, https://doi.org/10.5194/gmd-14-1899-2021, https://doi.org/10.5194/gmd-14-1899-2021, 2021
Short summary
Short summary
Computational models of Earth's mantle require rigorous verification and validation. Analytical solutions of the underlying Stokes equations provide a method to verify that these equations are accurately solved for. However, their derivation in spherical and cylindrical shell domains with physically relevant boundary conditions is involved. This paper provides a number of solutions. They are provided in a Python package (Assess) and their use is demonstrated in a convergence study with Fluidity.
Bastian van den Bout, Theo van Asch, Wei Hu, Chenxiao X. Tang, Olga Mavrouli, Victor G. Jetten, and Cees J. van Westen
Geosci. Model Dev., 14, 1841–1864, https://doi.org/10.5194/gmd-14-1841-2021, https://doi.org/10.5194/gmd-14-1841-2021, 2021
Short summary
Short summary
Landslides, debris flows and other types of dense gravity-driven flows threaten livelihoods around the globe. Understanding the mechanics of these flows can be crucial for predicting their behaviour and reducing disaster risk. Numerical models assume that the solids and fluids of the flow are unstructured. The newly presented model captures the internal structure during movement. This important step can lead to more accurate predictions of landslide movement.
Andrzej Górszczyk and Stéphane Operto
Geosci. Model Dev., 14, 1773–1799, https://doi.org/10.5194/gmd-14-1773-2021, https://doi.org/10.5194/gmd-14-1773-2021, 2021
Short summary
Short summary
We present the 3D multi-parameter synthetic geomodel of the subduction zone, as well as the workflow designed to implement all of its components. The model contains different geological structures of various scales and complexities. It is intended to serve as a tool for the geophysical community to validate imaging approaches, design acquisition techniques, estimate uncertainties, benchmark computing approaches, etc.
Mattia de' Michieli Vitturi and Federica Pardini
Geosci. Model Dev., 14, 1345–1377, https://doi.org/10.5194/gmd-14-1345-2021, https://doi.org/10.5194/gmd-14-1345-2021, 2021
Short summary
Short summary
Here, we present PLUME-MoM-TSM, a volcanic plume model that allows us to quantify the formation of aggregates during the rise of the plume, model the phase change of water, and include the possibility to simulate the initial spreading of the tephra umbrella cloud intruding from the volcanic column into the atmosphere. The model is first applied to the 2015 Calbuco eruption (Chile) and provides an analytical relationship between the upwind spreading and some characteristic of the volcanic column.
Zhikui Guo, Lars Rüpke, and Chunhui Tao
Geosci. Model Dev., 13, 6547–6565, https://doi.org/10.5194/gmd-13-6547-2020, https://doi.org/10.5194/gmd-13-6547-2020, 2020
Short summary
Short summary
We present the 3-D hydro-thermo-transport model HydrothermalFoam v1.0, which we designed to provide the marine geosciences community with an easy-to-use and state-of-the-art tool for simulating mass and energy transport in submarine hydrothermal systems. HydrothermalFoam is based on the popular open-source platform OpenFOAM, comes with a number of tutorials, and is published under the GNU General Public License v3.0.
Marisol Monterrubio-Velasco, F. Ramón Zúñiga, Quetzalcoatl Rodríguez-Pérez, Otilio Rojas, Armando Aguilar-Meléndez, and Josep de la Puente
Geosci. Model Dev., 13, 6361–6381, https://doi.org/10.5194/gmd-13-6361-2020, https://doi.org/10.5194/gmd-13-6361-2020, 2020
Short summary
Short summary
The Mexican subduction zone along the Pacific coast is one of the most active seismic zones in the world, where every year larger-magnitude earthquakes shake huge inland cities such as Mexico City. In this work, we use TREMOL (sThochastic Rupture Earthquake ModeL) to simulate the seismicity observed in this zone. Our numerical results reinforce the hypothesis that in some subduction regions single asperities are responsible for producing the observed seismicity.
Thomas Zwinger, Grace A. Nield, Juha Ruokolainen, and Matt A. King
Geosci. Model Dev., 13, 1155–1164, https://doi.org/10.5194/gmd-13-1155-2020, https://doi.org/10.5194/gmd-13-1155-2020, 2020
Short summary
Short summary
We present a newly developed flat-earth model, Elmer/Earth, for viscoelastic treatment of solid earth deformation under ice loads. Unlike many previous approaches with proprietary software, this model is based on the open-source FEM code Elmer, with the advantage for scientists to apply and alter the model without license constraints. The new-generation full-stress ice-sheet model Elmer/Ice shares the same code base, enabling future coupled ice-sheet–glacial-isostatic-adjustment simulations.
Swarup Chauhan, Kathleen Sell, Wolfram Rühaak, Thorsten Wille, and Ingo Sass
Geosci. Model Dev., 13, 315–334, https://doi.org/10.5194/gmd-13-315-2020, https://doi.org/10.5194/gmd-13-315-2020, 2020
Short summary
Short summary
We present CobWeb 1.0, a graphical user interface for analysing tomographic images of geomaterials. CobWeb offers different machine learning techniques for accurate multiphase image segmentation and visualizing material specific parameters such as pore size distribution, relative porosity and volume fraction. We demonstrate a novel approach of dual filtration and dual segmentation to eliminate edge enhancement artefact in synchrotron-tomographic datasets and provide the computational code.
Loïc Huder, Nicolas Gillet, and Franck Thollard
Geosci. Model Dev., 12, 3795–3803, https://doi.org/10.5194/gmd-12-3795-2019, https://doi.org/10.5194/gmd-12-3795-2019, 2019
Short summary
Short summary
The pygeodyn package is a geomagnetic data assimilation tool written in Python. It gives access to the Earth's core flow dynamics, controlled by geomagnetic observations, by means of a reduced numerical model anchored to geodynamo simulation statistics. It aims to provide the community with a user-friendly and tunable data assimilation algorithm. It can be used for education, geomagnetic model production or tests in conjunction with webgeodyn, a set of visualization tools for geomagnetic models.
Mattia de' Michieli Vitturi, Tomaso Esposti Ongaro, Giacomo Lari, and Alvaro Aravena
Geosci. Model Dev., 12, 581–595, https://doi.org/10.5194/gmd-12-581-2019, https://doi.org/10.5194/gmd-12-581-2019, 2019
Short summary
Short summary
Pyroclastic avalanches are a type of granular flow generated at active volcanoes by different mechanisms, including the collapse of steep pyroclastic deposits (e.g., scoria and ash cones) and fountaining during moderately explosive eruptions. We present IMEX_SfloW2D, a depth-averaged flow model describing the granular mixture as a single-phase granular fluid. Benchmark cases and preliminary application to the simulation of the 11 February pyroclastic avalanche at Mt. Etna (Italy) are shown.
Yihao Wu, Zhicai Luo, Bo Zhong, and Chuang Xu
Geosci. Model Dev., 11, 4797–4815, https://doi.org/10.5194/gmd-11-4797-2018, https://doi.org/10.5194/gmd-11-4797-2018, 2018
Short summary
Short summary
A multilayer approach is parameterized for model development, and the multiple layers are located at different depths beneath the Earth’s surface. This method may be beneficial for gravity/manget field modeling, which may outperform the traditional single-layer approach.
Andres Payo, Bismarck Jigena Antelo, Martin Hurst, Monica Palaseanu-Lovejoy, Chris Williams, Gareth Jenkins, Kathryn Lee, David Favis-Mortlock, Andrew Barkwith, and Michael A. Ellis
Geosci. Model Dev., 11, 4317–4337, https://doi.org/10.5194/gmd-11-4317-2018, https://doi.org/10.5194/gmd-11-4317-2018, 2018
Short summary
Short summary
We describe a new algorithm that automatically delineates the cliff top and toe of a cliffed coastline from a digital elevation model (DEM). The algorithm builds upon existing methods but is specifically designed to resolve very irregular planform coastlines with many bays and capes, such as parts of the coastline of Great Britain.
Hugo Cruz-Jiménez, Guotu Li, Paul Martin Mai, Ibrahim Hoteit, and Omar M. Knio
Geosci. Model Dev., 11, 3071–3088, https://doi.org/10.5194/gmd-11-3071-2018, https://doi.org/10.5194/gmd-11-3071-2018, 2018
Short summary
Short summary
One of the most important challenges seismologists and earthquake engineers face is reliably estimating ground motion in an area prone to large damaging earthquakes. This study aimed at better understanding the relationship between characteristics of geological faults (e.g., hypocenter location, rupture size/location, etc.) and resulting ground motion, via statistical analysis of a rupture simulation model. This study provides important insight on ground-motion responses to geological faults.
Fabio Crameri
Geosci. Model Dev., 11, 2541–2562, https://doi.org/10.5194/gmd-11-2541-2018, https://doi.org/10.5194/gmd-11-2541-2018, 2018
Short summary
Short summary
Firstly, this study acts as a compilation of key geodynamic diagnostics and describes how to automatise them for a more efficient scientific procedure. Secondly, it outlines today's key pitfalls of scientific visualisation and provides means to circumvent them with, for example, a novel set of fully scientific colour maps. Thirdly, it introduces StagLab 3.0, a software that applies such fully automated diagnostics and state-of-the-art visualisation in the blink of an eye.
Michael Bock, Olaf Conrad, Andreas Günther, Ernst Gehrt, Rainer Baritz, and Jürgen Böhner
Geosci. Model Dev., 11, 1641–1652, https://doi.org/10.5194/gmd-11-1641-2018, https://doi.org/10.5194/gmd-11-1641-2018, 2018
Short summary
Short summary
We introduce the Soil and
Landscape Evolution Model (SaLEM) for the prediction of soil parent material evolution following a lithologically differentiated approach. The GIS tool is working within the software framework SAGA GIS. Weathering, erosion and transport functions are calibrated using extrinsic and intrinsic parameter data. First results indicate that our approach shows evidence for the spatiotemporal prediction of soil parental material properties.
Karthik Iyer, Henrik Svensen, and Daniel W. Schmid
Geosci. Model Dev., 11, 43–60, https://doi.org/10.5194/gmd-11-43-2018, https://doi.org/10.5194/gmd-11-43-2018, 2018
Short summary
Short summary
Igneous intrusions in sedimentary basins have a profound effect on the thermal structure of the hosting sedimentary rocks. In this paper, we present a user-friendly 1-D FEM-based tool, SILLi, that calculates the thermal effects of sill intrusions on the enclosing sedimentary stratigraphy. The motivation is to make a standardized numerical toolkit openly available that can be widely used by scientists with different backgrounds to test the effects of magmatic bodies in a wide variety of settings.
Charles M. Shobe, Gregory E. Tucker, and Katherine R. Barnhart
Geosci. Model Dev., 10, 4577–4604, https://doi.org/10.5194/gmd-10-4577-2017, https://doi.org/10.5194/gmd-10-4577-2017, 2017
Short summary
Short summary
Rivers control the movement of sediment and nutrients across Earth's surface. Understanding how rivers change through time is important for mitigating natural hazards and predicting Earth's response to climate change. We develop a new computer model for predicting how rivers cut through sediment and rock. Our model is designed to be joined with models of flooding, landslides, vegetation change, and other factors to provide a comprehensive toolbox for predicting changes to the landscape.
Diego Takahashi and Vanderlei C. Oliveira Jr.
Geosci. Model Dev., 10, 3591–3608, https://doi.org/10.5194/gmd-10-3591-2017, https://doi.org/10.5194/gmd-10-3591-2017, 2017
Short summary
Short summary
Ellipsoids are the only bodies for which the self-demagnetization can be treated analytically. This property is useful for modelling compact orebodies having high susceptibility. We present a review of the magnetic modelling of ellipsoids, propose a way of determining the isotropic susceptibility above which the self-demagnetization must be considered, and discuss the ambiguity between confocal ellipsoids, as well as provide a set of routines to model the magnetic field produced by ellipsoids.
Cited articles
Amezcua, J. and Leeuwen, P. J. V.:
Gaussian anamorphosis in the analysis step of the EnKF: a joint state-variable/observation approach, Tellus A, 66, 23493, https://doi.org/10.3402/tellusa.v66.23493, 2014. a
Anderson, J. L.:
A Local Least Squares Framework for Ensemble Filtering, Mon. Weather Rev., 131, 634–642, https://doi.org/10.1175/1520-0493(2003)131<0634:ALLSFF>2.0.CO;2, 2003. a, b
Beckett, F. M., Witham, C. S., Leadbetter, S. J., Crocker, R., Webster, H. N., Hort, M. C., Jones, A. R., Devenish, B. J., and Thomson, D. J.:
Atmospheric Dispersion Modelling at the London VAAC: A Review of Developments since the 2010 Eyjafjallajökull Volcano Ash Cloud, Atmosphere, 11, 352, https://doi.org/10.3390/atmos11040352, 2020. a
Bonadonna, C. and Costa, A.:
Plume height, volume, and classification of volcanic explosive eruptions based on the Weibull function, B. Volcanol., 75, 1–19, https://doi.org/10.1007/s00445-013-0742-1, 2013. a
Bonadonna, C., Ernst, G., and Sparks, R.:
Thickness variations and volume estimates of tephra fall deposits: the importance of particle Reynolds number, J. Volcanol. Geoth. Res., 81, 173–187, https://doi.org/10.1016/S0377-0273(98)00007-9, 1998. a
Bonadonna, C., Folch, A., Loughlin, S., and Puempel, H.:
Future developments in modelling and monitoring of volcanic ash clouds: outcomes from the first IAVCEI-WMO workshop on Ash Dispersal Forecast and Civil Aviation, B. Volcanol., 74, 1–10, 2012. a
Bonadonna, C., Biass, S., and Costa, A.:
Physical characterization of explosive volcanic eruptions based on tephra deposits: Propagation of uncertainties and sensitivity analysis, J. Volcanol. Geoth. Res., 296, 80–100, https://doi.org/10.1016/j.jvolgeores.2015.03.009, 2015. a
Carrassi, A., Bocquet, M., Bertino, L., and Evensen, G.:
Data assimilation in the geosciences: An overview of methods, issues, and perspectives, WIREs Clim. Change, 9, e535, https://doi.org/10.1002/wcc.535, 2018. a, b
Clarkson, R. J., Majewicz, E. J., and Mack, P.:
A re-evaluation of the 2010 quantitative understanding of the effects volcanic ash has on gas turbine engines, P. I. Mech. Eng. G-J. Aer., 230, 2274–2291, https://doi.org/10.1177/0954410015623372, 2016. a
Constantinescu, R., White, J. T., Connor, C. B., Hopulele-Gligor, A., Charbonnier, S., Thouret, J.-C., Lindsay, J. M., and Bertin, D.:
Uncertainty Quantification of Eruption Source Parameters Estimated From Tephra Fall Deposits, Geophys. Res. Lett., 49, e2021GL097425, https://doi.org/10.1029/2021GL097425, 2022. a
Cornell, W., Carey, S., and Sigurdsson, H.:
Computer simulation and transport of the Campanian Y-5 ash, J. Volcanol. Geoth. Res., 17, 89–109, https://doi.org/10.1016/0377-0273(83)90063-X, 1983. a
Costa, A., Folch, A., and Macedonio, G.:
A model for wet aggregation of ash particles in volcanic plumes and clouds: I. Theoretical formulation, J. Geophys. Res.-Sol. Ea., 115, B09201, https://doi.org/10.1029/2009JB007175, 2010. a
Costa, A., Pioli, L., and Bonadonna, C.:
Assessing tephra total grain-size distribution: Insights from field data analysis, Earth Planet. Sc. Lett., 443, 90–107, https://doi.org/10.1016/j.epsl.2016.02.040, 2016. a
Dacre, H. F. and Harvey, N. J.:
Characterizing the Atmospheric Conditions Leading to Large Error Growth in Volcanic Ash Cloud Forecasts, J. Appl. Meteorol. Clim., 57, 1011–1019, https://doi.org/10.1175/JAMC-D-17-0298.1, 2018. a
Dare, R. A., Smith, D. H., and Naughton, M. J.:
Ensemble Prediction of the Dispersion of Volcanic Ash from the 13 February 2014 Eruption of Kelut, Indonesia, J. Appl. Meteorol. Clim., 55, 61–78, https://doi.org/10.1175/JAMC-D-15-0079.1, 2016. a
Degruyter, W. and Bonadonna, C.:
Improving on mass flow rate estimates of volcanic eruptions, Geophys. Res. Lett., 39, L16308, https://doi.org/10.1029/2012GL052566, 2012. a, b
Denlinger, R. P., Pavolonis, M., and Sieglaff, J.:
A robust method to forecast volcanic ash clouds, J. Geophys. Res.-Atmos., 117, D13208, https://doi.org/10.1029/2012JD017732, 2012. a
Dominguez, L., Bonadonna, C., Forte, P., Jarvis, P. A., Cioni, R., Mingari, L., Bran, D., and Panebianco, J. E.:
Aeolian remobilisation of the 2011-Cordón Caulle Tephra-Fallout Deposit: example of an important process in the life cycle of Volcanic Ash, Front. Earth Sci., 7, 343, https://doi.org/10.3389/feart.2019.00343, 2020. a
Durant, A. J., Rose, W. I., Sarna-Wojcicki, A. M., Carey, S., and Volentik, A. C. M.:
Hydrometeor-enhanced tephra sedimentation: Constraints from the 18 May 1980 eruption of Mount St. Helens, J. Geophys. Res.-Sol. Ea., 114, B03204, https://doi.org/10.1029/2008JB005756, 2009. a
Evensen, G.:
Sequential data assimilation with a nonlinear quasi-geostrophic model using Monte Carlo methods to forecast error statistics, J. Geophys. Res.-Oceans, 99, 10143–10162, https://doi.org/10.1029/94JC00572, 1994. a, b
Folch, A. and Mingari, L.: Data Assimilation of Volcanic Clouds: Recent Advances and Implications on Operational Forecasts, in: Applications of Data Assimilation and Inverse Problems in the Earth Sciences (Special Publications of the International Union of Geodesy and Geophysics, edited by: Ismail-Zadeh, A., Castelli, F., Jones, D., and Sanchez, S., 144–156, Cambridge University Press, Cambridge, https://doi.org/10.1017/9781009180412.010, 2023. a, b, c
Folch, A., Costa, A., Durant, A., and Macedonio, G.:
A model for wet aggregation of ash particles in volcanic plumes and clouds: II. Model application, J. Geophys. Res.-Sol. Ea., 115, B09202, https://doi.org/10.1029/2009JB007176, 2010. a
Folch, A., Mingari, L., Osores, M. S., and Collini, E.:
Modeling volcanic ash resuspension – application to the 14–18 October 2011 outbreak episode in central Patagonia, Argentina, Nat. Hazards Earth Syst. Sci., 14, 119–133, https://doi.org/10.5194/nhess-14-119-2014, 2014. a
Folch, A., Mingari, L., Gutierrez, N., Hanzich, M., Macedonio, G., and Costa, A.:
FALL3D-8.0: a computational model for atmospheric transport and deposition of particles, aerosols and radionuclides – Part 1: Model physics and numerics, Geosci. Model Dev., 13, 1431–1458, https://doi.org/10.5194/gmd-13-1431-2020, 2020. a, b
Folch, A., Costa, A., Macedonio, G. and Mingari, L.: FALL3D (8.2.0), Zenodo [code], https://doi.org/10.5281/zenodo.7267808, 2022a. a
Folch, A., Mingari, L., and Prata, A. T.:
Ensemble-Based Forecast of Volcanic Clouds Using FALL3D-8.1, Front. Earth Sci., 9, 741841, https://doi.org/10.3389/feart.2021.741841, 2022b. a, b
Folch, A., Abril, C., Afanasiev, M., Amati, G., Bader, M., Badia, R. M., Bayraktar, H. B., Barsotti, S., Basili, R., Bernardi, F., Boehm, C., Brizuela, B., Brogi, F., Cabrera, E., Casarotti, E., Castro, M. J., Cerminara, M., Cirella, A., Cheptsov, A., Conejero, J., Costa, A., de la Asunción, M., de la Puente, J., Djuric, M., Dorozhinskii, R., Espinosa, G., Esposti-Ongaro, T., Farnós, J., Favretto-Cristini, N., Fichtner, A., Fournier, A., Gabriel, A.-A., Gallard, J.-M., Gibbons, S. J., Glimsdal, S., González-Vida, J. M., Gracia, J., Gregorio, R., Gutierrez, N., Halldorsson, B., Hamitou, O., Houzeaux, G., Jaure, S., Kessar, M., Krenz, L.,Krischer, L., Laforet, S., Lanucara, P., Li, B., Lorenzino, M. C., Lorito, S., Løvholt, F., Macedonio, G., Macías, J., Marín, G., Martínez Montesinos, B., Mingari, L., Moguilny, G., Montellier, V., Monterrubio-Velasco, M., Moulard, G. E., Nagaso, M., Nazaria, M., Niethammer, C., Pardini, F., Pienkowska, M., Pizzimenti, L., Poiata, N., Rannabauer, L., Rojas, O., Rodriguez, J. E., Romano, F., Rudyy, O., Ruggiero, V., Samfass, P., Sánchez-Linares, C., Sanchez, S., Sandri, L., Scala, A., Schaeffer, N., Schuchart, J., Selva, J., Sergeant, A., Stallone, A., Taroni, M., Thrastarson, S., Titos, M., Tonelllo, N., Tonini, R., Ulrich, T., Vilotte, J.-P., Vöge, M., Volpe, M., Wirp, S. A., and Wössner, U.:
The EU Center of Excellence for Exascale in Solid Earth (ChEESE): Implementation, results, and roadmap for the second phase, Future Gener. Comp. Sy., 146, 47–61, https://doi.org/10.1016/j.future.2023.04.006, 2023. a
Fu, G., Heemink, A., Lu, S., Segers, A., Weber, K., and Lin, H.-X.:
Model-based aviation advice on distal volcanic ash clouds by assimilating aircraft in situ measurements, Atmos. Chem. Phys., 16, 9189–9200, https://doi.org/10.5194/acp-16-9189-2016, 2016. a
Fu, G., Prata, F., Lin, H. X., Heemink, A., Segers, A., and Lu, S.:
Data assimilation for volcanic ash plumes using a satellite observational operator: a case study on the 2010 Eyjafjallajökull volcanic eruption, Atmos. Chem. Phys., 17, 1187–1205, https://doi.org/10.5194/acp-17-1187-2017, 2017. a
Harvey, N. J., Dacre, H. F., Webster, H. N., Taylor, I. A., Khanal, S., Grainger, R. G., and Cooke, M. C.:
The impact of ensemble meteorology on inverse modeling estimates of volcano emissions and ash dispersion forecasts: Grímsvötn 2011, Atmosphere, 11, 1022, https://doi.org/10.3390/atmos11101022, 2020. a
Husak, G. J., Michaelsen, J., and Funk, C.:
Use of the gamma distribution to represent monthly rainfall in Africa for drought monitoring applications, Int. J. Climatol., 27, 935–944, https://doi.org/10.1002/joc.1441, 2007. a
Jazwinski, A. H.:
Stochastic processes and filtering theory, Academic Press, Inc., New York, ISBN 9780080960906, 1970. a
Kalman, R. E.:
A New Approach to Linear Filtering and Prediction Problems, J. Basic Eng.-T. ASME, 82, 35–45, https://doi.org/10.1115/1.3662552, 1960. a, b
Kliewer, A. J., Fletcher, S. J., Jones, A. S., and Forsythe, J. M.:
Comparison of Gaussian, logarithmic transform and mixed Gaussian–log-normal distribution based 1DVAR microwave temperature–water-vapour mixing ratio retrievals, Q.+J. Roy. Meteor. Soc., 142, 274–286, https://doi.org/10.1002/qj.2651, 2016. a
Madankan, R., Pouget, S., Singla, P., Bursik, M., Dehn, J., Jones, M., Patra, A., Pavolonis, M., Pitman, E., Singh, T., and Webley, P.:
Computation of probabilistic hazard maps and source parameter estimation for volcanic ash transport and dispersion, J. Comput. Phys., 271, 39–59, https://doi.org/10.1016/j.jcp.2013.11.032, 2014. a
Martí, A., Folch, A., Costa, A., and Engwell, S.:
Reconstructing the plinian and co-ignimbrite sources of large volcanic eruptions: a novel approach for the Campanian Ignimbrite, Sci. Rep., 6, 1–11, https://doi.org/10.1038/srep21220, 2016. a, b
Martinez, B. M., Luzón, M. T., Sandri, L., Rudyy, O., Cheptsov, A., Macedonio, G., Folch, A., Barsotti, S., Selva, J., and Costa, A.:
On the feasibility and usefulness of high-performance computing in probabilistic volcanic hazard assessment: An application to tephra hazard from Campi Flegrei, Front. Earth Sci., 10, 941789, https://doi.org/10.3389/feart.2022.941789, 2022. a
McKay, M. D., Beckman, R. J., and Conover, W. J.:
A Comparison of Three Methods for Selecting Values of Input Variables in the Analysis of Output from a Computer Code, Technometrics, 21, 239–245, 1979. a
Mingari, L., Folch, A., Dominguez, L., and Bonadonna, C.:
Volcanic Ash Resuspension in Patagonia: Numerical Simulations and Observations, Atmosphere, 11, 977, https://doi.org/10.3390/atmos11090977, 2020. a
Mingari, L., Costa, A., Macedonio, G., and Folch, A.: gmd-2022-246 (2.0), Zenodo [data set], https://doi.org/10.5281/zenodo.7259531, 2022a. a
Mingari, L., Folch, A., Prata, A. T., Pardini, F., Macedonio, G., and Costa, A.:
Data assimilation of volcanic aerosol observations using FALL3D+PDAF, Atmos. Chem. Phys., 22, 1773–1792, https://doi.org/10.5194/acp-22-1773-2022, 2022b. a, b, c
O'Neill, N. T., Ignatov, A., Holben, B. N., and Eck, T. F.:
The lognormal distribution as a reference for reporting aerosol optical depth statistics; Empirical tests using multi-year, multi-site AERONET Sunphotometer data, Geophys. Res. Lett., 27, 3333–3336, https://doi.org/10.1029/2000GL011581, 2000. a
Osores, S., Ruiz, J., Folch, A., and Collini, E.:
Volcanic ash forecast using ensemble-based data assimilation: an ensemble transform Kalman filter coupled with the FALL3D-7.2 model (ETKF–FALL3D version 1.0), Geosci. Model Dev., 13, 1–22, https://doi.org/10.5194/gmd-13-1-2020, 2020. a
Pardini, F., Corradini, S., Costa, A., Esposti Ongaro, T., Merucci, L., Neri, A., Stelitano, D., de' Michieli Vitturi, M.: Ensemble-Based Data Assimilation of Volcanic Ash Clouds from Satellite Observations: Application to the 24 December 2018 Mt. Etna Explosive Eruption, Atmosphere, 11, 359, https://doi.org/10.3390/atmos11040359, 2020. a
Pelley, R., Cooke, M., Manning, A., Thomson, D., Witham, C., and Hort, M.:
Initial implementation of an inversion technique for estimating volcanic ash source parameters in near real time using satellite retrievals, Tech. rep., Forecasting Research Technical Report, Met Office, Exeter, 2015. a
Pfeiffer, T., Costa, A., and Macedonio, G.:
A model for the numerical simulation of tephra fall deposits, J. Volcanol. Geoth. Res., 140, 273–294, https://doi.org/10.1016/j.jvolgeores.2004.09.001, 2005. a, b, c
Posselt, D. J. and Bishop, C. H.:
Nonlinear data assimilation for clouds and precipitation using a gamma inverse-gamma ensemble filter, Q. J. Roy. Meteor. Soc., 144, 2331–2349, https://doi.org/10.1002/qj.3374, 2018. a
Prata, A. T., Mingari, L., Folch, A., Macedonio, G., and Costa, A.:
FALL3D-8.0: a computational model for atmospheric transport and deposition of particles, aerosols and radionuclides – Part 2: Model validation, Geosci. Model Dev., 14, 409–436, https://doi.org/10.5194/gmd-14-409-2021, 2021. a, b
Pyle, D. M.:
The thickness, volume and grainsize of tephra fall deposits, B. Volcanol., 51, 1–15, 1989. a
Romero, J., Morgavi, D., Arzilli, F., Daga, R., Caselli, A., Reckziegel, F., Viramonte, J., Díaz-Alvarado, J., Polacci, M., Burton, M., and Perugini, D.:
Eruption dynamics of the 22–23 April 2015 Calbuco Volcano (Southern Chile): Analyses of tephra fall deposits, J. Volcanol. Geoth. Res., 317, 15–29, https://doi.org/10.1016/j.jvolgeores.2016.02.027, 2016. a, b, c, d, e, f, g, h
Sandri, L., Costa, A., Selva, J., Tonini, R., Macedonio, G., Folch, A., and Sulpizio, R.:
Beyond eruptive scenarios: Assessing tephra fallout hazard from Neapolitan volcanoes, Sci. Rep., 6, 1–13, https://doi.org/10.1038/srep24271, 2016. a
Scasso, R. A., Corbella, H., and Tiberi, P.:
Sedimentological analysis of the tephra from the 12–15 August 1991 eruption of Hudson volcano, B. Volcanol., 56, 121–132, 1994. a
Sha, F., Lin, Y., Saul, L. K., and Lee, D. D.:
Multiplicative Updates for Nonnegative Quadratic Programming, Neural Comput., 19, 2004–2031, https://doi.org/10.1162/neco.2007.19.8.2004, 2007. a, b
Stefanescu, E. R., Patra, A. K., Bursik, M. I., Madankan, R., Pouget, S., Jones, M., Singla, P., Singh, T., Pitman, E. B., Pavolonis, M., Morton, D., Webley, P., and Dehn, J.:
Temporal, probabilistic mapping of ash clouds using wind field stochastic variability and uncertain eruption source parameters: Example of the 14 April 2010 Eyjafjallajökull eruption, J. Adv. Model. Earth Sy., 6, 1173–1184, https://doi.org/10.1002/2014MS000332, 2014. a
Suzuki, T.:
A theoretical model for dispersion of tephra, in: Arc Volcanism: Physics and Tectonics, edited by: Shimozuru, D. and Yokoyama, I., Terra Scientific Publishing Company (TERRAPUB), Tokyo, 93–113, ISBN 90-277-1612-9, 1983. a
Van Eaton, A. R., Amigo, A., Bertin, D., Mastin, L. G., Giacosa, R. E., González, J., Valderrama, O., Fontijn, K., and Behnke, S. A.:
Volcanic lightning and plume behavior reveal evolving hazards during the April 2015 eruption of Calbuco volcano, Chile, Geophys. Res. Lett., 43, 3563–3571, https://doi.org/10.1002/2016GL068076, 2016.
a, b, c, d, e
Vogel, H., Förstner, J., Vogel, B., Hanisch, T., Mühr, B., Schättler, U., and Schad, T.:
Time-lagged ensemble simulations of the dispersion of the Eyjafjallajökull plume over Europe with COSMO-ART, Atmos. Chem. Phys., 14, 7837–7845, https://doi.org/10.5194/acp-14-7837-2014, 2014. a
Wilson, G., Wilson, T., Deligne, N., and Cole, J.:
Volcanic hazard impacts to critical infrastructure: A review, J. Volcanol. Geoth. Res., 286, 148–182, https://doi.org/10.1016/j.jvolgeores.2014.08.030, 2014. a
Zidikheri, M. J. and Lucas, C.:
A Computationally Efficient Ensemble Filtering Scheme for Quantitative Volcanic Ash Forecasts, J. Geophys. Res.-Atmos., 126, e2020JD033094, https://doi.org/10.1029/2020JD033094, 2021a. a
Zidikheri, M. J. and Lucas, C.:
Improving Ensemble Volcanic Ash Forecasts by Direct Insertion of Satellite Data and Ensemble Filtering, Atmosphere, 12, 1215, https://doi.org/10.3390/atmos12091215, 2021b. a
Zidikheri, M. J., Lucas, C., and Potts, R. J.:
Estimation of optimal dispersion model source parameters using satellite detections of volcanic ash, J. Geophys. Res.-Atmos., 122, 8207–8232, https://doi.org/10.1002/2017JD026676, 2017. a
Zidikheri, M. J., Lucas, C., and Potts, R. J.:
Quantitative Verification and Calibration of Volcanic Ash Ensemble Forecasts Using Satellite Data, J. Geophys. Res.-Atmos., 123, 4135–4156, https://doi.org/10.1002/2017JD027740, 2018. a
Short summary
Two novel techniques for ensemble-based data assimilation, suitable for semi-positive-definite variables with highly skewed uncertainty distributions such as tephra deposit mass loading, are applied to reconstruct the tephra fallout deposit resulting from the 2015 Calbuco eruption in Chile. The deposit spatial distribution and the ashfall volume according to the analyses are in good agreement with estimations based on field measurements and isopach maps reported in previous studies.
Two novel techniques for ensemble-based data assimilation, suitable for semi-positive-definite...
