Abstract

Abstract. This paper presents model validation results for the latest version release of the FALL3D atmospheric transport model. The code has been redesigned from scratch to incorporate different categories of species and to overcome legacy issues that precluded its preparation towards extreme-scale computing. The model validation is based on the new FALL3D-8.0 test suite, which comprises a set of four real case studies that encapsulate the major features of the model; namely, the simulation of long-range fine volcanic ash dispersal, volcanic SO2 dispersal, tephra fallout deposits and the dispersal and deposition of radionuclides. The first two test suite cases (i.e. the June 2011 Puyehue-Cordón Caulle ash cloud and the June 2019 Raikoke SO2 cloud) are validated against geostationary satellite retrievals and demonstrate the new FALL3D data insertion scheme. The metrics used to validate the volcanic ash and SO2 simulations are the structure, amplitude and location (SAL) metric and the figure of merit in space (FMS). The other two test suite cases (i.e. the February 2013 Mt. Etna ash cloud and associated tephra fallout deposit, and the dispersal of radionuclides resulting from the 1986 Chernobyl nuclear accident) are validated with scattered ground-based observations of deposit load and local particle grain size distributions and with measurements from the Radioactivity Environmental Monitoring database. For validation of tephra deposit loads and radionuclides, we use two variants of the normalised root-mean-square error metric. We find that FALL3D-8.0 simulations initialised with data insertion consistently improve agreement with satellite retrievals at all lead times up to 48 h for both volcanic ash and SO2 simulations. In general, SAL scores lower than 1.5 and FMS scores greater than 0.40 indicate acceptable agreement with satellite retrievals of volcanic ash and SO2. In addition, we show very good agreement, across several orders of magnitude, between the model and observations for the 2013 Mt. Etna and 1986 Chernobyl case studies. Our results, along with the validation datasets provided in the publicly available test suite, form the basis for future improvements to FALL3D (version 8 or later) and also allow for model intercomparison studies.

Highlights

  • FALL3D-8.0 is the latest major version release of FALL3D (Costa et al, 2006; Folch et al, 2009), an open-source code with a 15-year+ track record and a growing number of users in the volcanological and atmospheric science communities

  • In order to validate FALL3D-8.0 simulations of dispersal of fine ash and to test the new volcanic ash data insertion scheme, we use the retrieval method of Prata and Prata (2012) to derive fine ash mass loading estimates based on infrared (IR) measurements made by the SEVIRI (Spinning Enhanced Visible and Infrared Imager; Schmetz et al, 2002) instrument during the 2011 PuyehueCordón Caulle eruption in Chile

  • The Chernobyl-1986 case considers the computational domain shown in Fig. 8 for the period from 24 April to 10 May 1986, considering the input values reported in Table 5 and the meteorological fields obtained from ERA5 reanalysis, which accounts for atmospheric diffusion, wet deposition and radioactive decay

Read more

Summary

Introduction

FALL3D-8.0 is the latest major version release of FALL3D (Costa et al, 2006; Folch et al, 2009), an open-source code with a 15-year+ track record and a growing number of users in the volcanological and atmospheric science communities. This subfolder contains all the necessary files to obtain meteorological data depending on the meteorological driver (for possible options, see Table 12 in Folch et al, 2020). For global datasets such as the ERA5 dataset (Copernicus Climate Change Service, 2017) used in the Puyehue-2011 case, Python and shell scripts are provided so that the user can download and merge meteorological data consistently with the SetDbs pre-process task Python script to validate model results; writes SAL metrics on validation_metrics_puyehue.txt Python library needed (imported) by validate_puyehue.py

Puyehue-2011
Validation dataset
Model setup
Raikoke-2019
Etna-2013
Chernobyl-1986
Validation metrics
Amplitude
Location
Structure
Figure of merit in space
Root-mean-square error
Findings
Conclusions
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call