The NEMOH network offers to appointed researchers the unique opportunity to design and develop a personalized training-through-research route, based on research tools representing state-of-the-art in numerical modelling, experimental facilities, and hazard evaluation methods applied to volcano science.


INGV - Istituto Nazionale di Geofisica e Vulcanologia

 Name  Description


Code for computation of the fully non-ideal saturation surface of H2O-CO2-bearing silicate melts of any natural composition


Package of codes for the determination of confidence regions of total volatiles in magmas from melt inclusion data, based on Monte Carlo simulations


Code for the computation of the oxidation state, and non-ideal gas-phase homogeneous and gas-liquid heterogeneous equilibrium reactions between natural silicate melts of any composition and volatiles in the C-O-H-S system



Code for the computation of the steady, 1D, isothermal, multiphase nonequilibrium dynamics of magma ascent through volcanic conduits and fissures


Code for the computation of the transient, 2D dynamics of multi-component fluids in regimes from low-Re, incompressible, to high-Re, compressible, including transonic flow and shock waves

 PDAC-2D and PDAC-3D

Codes for the computation of the transient, 2D or 3D, multiphase multi-component dynamics of gas-particle dispersions and pyroclastic flow generation and movement


Lagrangian code for the computation of the trajectory of large volcanic particles ejected during explosive eruptions


Code for the modeling of the steady-state and transient dynamics of magma flow along geometry-variable conduits of dome systems


Code for the simulation of the transient and 3D dispersal dynamics of volcanic plume and ash cloud


Code for the simulation of the 1D quasi-steady-state dynamics of magma chamber plus conduit flow


Code for short- and long-term eruption forecasting. The code accounts for the past history of the volcano and the monitoring observations. The output consists of a probabilistic estimation and related uncertainty of any possible event chosen by the end-user

 Hydraulic press

840 tons hydraulic press (multi-anvil and piston cylinder), for exploration of material properties
and behaviour in a wide range of P-T conditions

 Jeol JXA 8200

Jeol JXA 8200 super-probe for quantitative analysis of glass and minerals, equipped with five wavelength dispersive spectrometers, one energy dispersive spectrometer and transmission illuminator
 JSM-6500F JSM-6500F high performance analytical field emission SEM integrating JEOL’s energy dispersive X-ray analyzer for simple, speedy morphological observation and elemental analysis
of volcanic product

 Uniaxial testing machine

Uniaxial testing machine (250kN, T max 1000 °C) for investigation of the physical properties of volcanic rocks in a wide range of P-T conditions


 Permeameter (up to 200 MPa effective pressure) equipped with double intensifier

 Digital high-speed cameras

Digital high-speed cameras CORDIN and NAC, for analysis of analogue experiments in volcanology and rock mechanics

 Linux Cluster

 Two local high-Performance Linux Clusters at INGV Pisa, one with 94 CPUs at 2.4 GHz with total peak performance of 460 Gflops/s, the other with 120 CPUs at 2.4 GHz with total peak performance of 580 Gflops/s

UCD - University College Dublin, National University of Ireland, Dublin

* these codes codes developed by UCD, along with user instructions are available on request to all the partners of NEMOH Consortium

 Name  Description


 Code for full moment-tensor inversion in frequency domain + grid search for the source location (


 Elastic and visco-elastic lattice method codes for computing wave propagation in 2D and 3D heterogeneous Earth models. Theory, performance and numerical banchmarking of these codes can be found in O'Brien et al., 2004, O'Brien and Bean 2004, O'Brien 2008, O'Brien et al., 2009

Linux Claster

 High Performance Linux Clusters at UCD with 96 CPU's at 2.4 Ghz with 2Gb RAM per CPU

Finite diff. scheme

 Finite-difference scheme for computing seismic wave propagation in 2D and 3D visco-elastic heterogeneous earth models (velocity-stress, displacement stress and rotated staggered grid formulations)

Field Equipment

 The core of the UCD Geophysics equipment pool comprises 56 broadband seismometers which are mainly used for temporay installations related to volcanological and microseism studies

- O'Brien, Gareth S., Bean, Christopher J., A 3D discrete numerical elastic lattice method for seismic wave propagation in heterogeneous media with topography, Geophys. Res. Lett., Vol. 31, No. 14, L14608 10.1029/2004GL020069 24 July 2004.

- O'Brien, G. S., Bean, C. J., A discrete numerical method for modeling volcanic earthquake source mechanisms, J. Geophys. Res., Vol. 109, No. B9, B09301 10.1029/2004JB003023 02 September 2004.

- O'Brien, G. S. (2008), Discrete visco-elastic lattice methods for seismic wave propagation, Geophys. Res. Lett., 35, L02302, doi:10.1029/2007GL032214.

- O'Brien G.S., Bean C.J. and Tapamo H. (2009), Dispersion analysis and computational efficiency of Elastic Lattice methods for seismic wave propagation, Computers & Geosciences, Vol. 35, 1768-1775.

DEES - Ludwig-Maximilians Universitaet Muenchen

 Name  Description

 Magma rheology

High temperature magma rheology facilities for the experimental construction of constitutive flow laws for magma. (3 x concentric cylinder, parallel plate, micropenetration, acoustic emissions sensoring)

 Fragmentation Lab

High-temperature, high-pressure magma fragmentation facility for mapping the dynamics and mechanisms of explosive volcanic eruptions (4 shock tubes, capabilities for complete particle characterisation and dynamic monitoring)

 High-speed camera

Digital high-speed camera from Vision Research to analyze fragmentation & ejection behavior of rocks and magma analogues

 Thermal Analysis Lab

Thermal analysis laboratory for obtaining thermal histories and high temperature enthalpic and volumetric properties of magma systems to magmatic temperatures (full scanning calorimetric and scanning dilatometric (horizontal and vertical; thermogravimetry)

 Computed Tomography

Tomographic texture analysis via micro-CT and neutron-CT

 High Temperature - Melt Lab

Dynamic high temperature facilities for the forced convective mingling of magmas and centrifuge for the high temperature separation of unmixed liquids

DESB - University of Bristol

 Name  Description

 Rotational rheometer

Haake Mars rotational rheometer with integrated microscope and digital camera for imaging fluid microstructures


Aerator and mixer for generation of multiphase analogue materials

 Shock tube apparatus

Decompression and shock tube apparatus (including high speed camera) for multiphase
analogue volcanological flows
  •  Experiment benchmarks that can be done at the University of Bristol for comparison with theoretical and numerical modelling

ULANC - Lancaster University

 Name  Description

 Laboratory flow tubes

A suite of laboratory tube systems for fluid-dynamic experiments, with pressure measurement from zero to one megasample per second, combined with imaging, displacement, force and acceleration sensing

 X-ray tomography

X-ray tomography facilities for study of solid material resulting from flow experiments as well as characterisation facilities for analogue fluids


Computational fluid dynamics simulations of conduit flow with Flow3D

 Laser scanner

Very-long-range laser scanner (Riegl LPM-321) for topographic measurements at up to ca. 4 km


Mass spectrometry

Thermogravimetry-differential scanning calorimetry facilities for measurement of volatile contents of volcanic rock samples and characterisation of thermal transformations

 Microscope hotstage

Hotstage microscopy facility (up to 1500 °C) for experimental melting, vesiculation and crystallization of volcanic rock samples and high-quality imaging of the products

BSC - Barcelona Supercomputing Center – Centro Nacional de Supercomputacion

 Name  Description



Code for the computation of ash-cloud transport and fallout


Code for 3D computation of ground deformation in elastic and thermo-elastic media. FEMFES is a program to solve 2D, axisymmetric or 3D equations of linear elasticity, heat transfer, linear thermoelasticity or viscosity for both homogeneous and nonhomogeneous materials. In any case, the system of governing equations is solved by means of a Finite Element Method (FEM)


Code to solve a wide variety of Computational Mechanics problems (fluid, structure,
multiphase, etc.) of HPC

 IBM cluster

High-Performance IBM cluster at BSC Barcelona, with 10000 CPUs and 94 Tflops/s of peak performance. This capacity is foreseen to increase to 1 Pflop/s by the end of 2011

 Bull CUDA cluster

Bull CUDA cluster composed by 128 Bull Blade B505 and 2 NVIDIA M2090

IES - Haskoli Islands

 Name  Description

 GPS analysis

Facilities for analysis GPS ground deformation data


Facilities for interferometric analysis of synthetic aperture radar data (InSAR)

Facilities for interpretation of deformation data (GPS, InSAR, tilt, etc.) including 2D and 3D
finite element modeling (FEM)


Facilities for simulation of multiphase fluid flow and heat transfer in geothermal reservoirs as fractured porous media, including TOUGH2


A full Stokes finite element flow model for glaciers, applicable to volcano-ice interaction

BRGM - Bureau de Recherches Geologiques et Minieres

 Name  Description

 Linux Clusters

High-Performance Linux Clusters at BRGM Orléans, 320 cores and 1.6 Tflops (64 AMD Opteron Quad core processors + 32 AMD Opteron Dual core processors)


Spectral-element method code for modeling seismic wave propagation in 3D complex media (viscoelastic + topography + large velocity gradient)


3D full waveform tomography by spectral-element method (code EFISPEC3D) and genetic algorithm optimization


Interferometric processor for InSAR analysis


Interferometric Point Target Analysis for persistent scatterers InSAR

UNIGE - Universite de Geneve

 Name  Description


Code for probabilistic prediction of tephra dispersion


GIS-Matlab based tool for the hazard and risk assesment of tephra fall

 Fluid dynamics lab

Sedimentology and geophysical fluid dynamics lab including a vertical wind tunnel for the study of particle terminal velocity and aggregation and a long flume for the study of granularflow propagation. Measuring tools include a Phantom v10 high speed camera (14-bit image depth, and 480 frames per second at a full resolution of 2400 x 1800 pixels)