## 3 Results  
The dataset resulting from the automatic tagging step is composed of a total of 5017 parameter values for the LGS dataset and 669 for INGV dataset. Details are reported in Tables 4 and 5. By exploring the dataset, we can observe that LGS texts give evidence of SO\({}_{2}\) parameter values just starting from July 2017, whereas CO\({}_{2}\) and C/S values start from November 2019. Most of the INGV missing values are due to text report unavailability (during July-December 2015, 2016, 2017 and 2018). The INGV CO\({}_{2}\) parameter starting from 2019 is poorly present because the station was destroyed after the paroxysm of July 3 rd. Satellite information from MODIS data actually occurs in reports starting from 2019. \begin{table}
\begin{tabular}{c c c c c} \hline \hline
**Year** & **SO\({}_{2}\) (t/d)** & **CO\({}_{2}\) (t/d)** & **C/S** & **VLP (Events/h)** & **MODIS (MW)** \\ \hline
2015 & 0 & 0 & 0 & 358 & 22 \\
2016 & 0 & 0 & 0 & 180 & 3 \\
2017 & 176 & 0 & 0 & 390 & 43 \\
2018 & 396 & 0 & 0 & 383 & 34 \\
2019 & 366 & 45 & 34 & 365 & 160 \\
2020 & 383 & 291 & 259 & 415 & 145 \\
2021 & 149 & 123 & 103 & 156 & 38 \\ \hline \hline \end{tabular}
\end{table} 
Table 4: Number of parameter–value pairs for the LGS dataset (both weekly and daily). Figure 5: Validation of automatic tagging per typology (the six columns list: row number, unique ID, parameter name, and the values found in the LGS daily, LGS weekly and INGV weekly bulletins). \begin{table}
\begin{tabular}{c c} \hline \hline
**Parameter** & **Value** \\ \hline VLP & 11 events/h \\ SO\({}_{2}\) & 39 t/d \\ CO\({}_{2}\) & MEDIUM \\ CO\({}_{2}\) & 727 t/d \\ MODIS & 1 MW \\ \hline \hline \end{tabular}
\end{table} 
Table 3: Tagging of the main parameters from the daily LGS bulletin for 4 March 2021. Time series for parameters concerning the whole period 2015-2021 have been generated. LGS values for daily reports (see Figure 6a,c) are considered, because these are more numerous than weekly ones. The LGS dataset presents a gap in the range July-December 2016 due to report unavailability. INGV values regard weekly reports (see Figure 6b,d). With regard to INGV, missing values are due to text report unavailability (as explained above). \begin{table}
\begin{tabular}{c c c c c c} \hline \hline
**Year** & **SO\({}_{2}\) (t/d)** & **CO\({}_{2}\) (t/d)** & **C/S** & **VLP (Events/h)** & **MODIS (MW)** \\ \hline
2015 & 23 & 21 & 1 & 24 & 0 \\
2016 & 25 & 25 & 1 & 25 & 0 \\
2017 & 23 & 23 & 4 & 23 & 0 \\
2018 & 25 & 25 & 0 & 25 & 0 \\
2019 & 52 & 18 & 2 & 37 & 46 \\
2020 & 51 & 0 & 9 & 51 & 52 \\
2021 & 15 & 0 & 9 & 19 & 5 \\ \hline \hline \end{tabular}
\end{table} 
Table 5: Number of parameter-value pairs for the INGV dataset (both weekly and daily). Figure 6: (**a**) LGS, VLP parameter (events/h); (**b**) INGV, VLP parameter (events/h); (**c**) LGS, SO\({}_{2}\) parameter (tons/day); (**d**) INGV, SO\({}_{2}\) parameter (tons/day). In order to verify the accuracy of the automatic extraction method from a qualitative point of view, the extracted data were analyzed in the context of volcanic activity. Firstly, we performed a correlation analysis among the three main parameters: SO\({}_{2}\), CO\({}_{2}\) and VLP. We present results for the LGS dataset. The INGV correlation analysis dataset is omitted because of the poor availability of data. Figure 7 shows the SO\({}_{2}\) and CO\({}_{2}\) flux and VLP time series recorded for LGS daily in the period of November 2019-April 2021. In this time slot, a correlation analysis can be performed because the co-occurrence of values for the three parameters on the same date is highly supported. The observable quantities, i.e., SO\({}_{2}\) integrated mass flux and CO\({}_{2}\) mass flux, represented in the graph show a coherent behavior over time (Figure 8).