<p>Hydrogen (H<sub>2</sub>) is one of the most abundant trace species in volcano-hydrothermal systems and is a key participant in many redox reactions occurring in the hydrothermal reservoir gas. Although H<sub>2</sub> can be produced in soils by N<sub>2</sub>-fixing and fertilizing bacteria, soils are considered nowadays as sinks of molecular hydrogen (Smith-Downey et al. 2006). Because of its chemical and physical characteristics, H<sub>2</sub> generated within the crust moves rapidly and escapes to the atmosphere. These characteristics make H<sub>2</sub> one of the best geochemical indicators of magmatic and geothermal activity at depth. Cumbre Vieja volcano (La Palma, Canary Islands) is the most active basaltic volcano in the Canaries with seven historical eruptions being Tenegu&#237;a eruption (1971) the most recent one. Cumbre Vieja volcano is characterized by a main north&#8211;south rift zone 20 km long, up to 1950 m in elevation and covering an area of 220 km<sup>2</sup> with vents located at the northwest and northeast. Cumbre Vieja does not show any visible degassing (fumaroles, etc.). For that reason, the geochemical volcano monitoring program at Cumbre Vieja volcano has been focused on soil degassing surveys. &#160;Here we show the results of soil H<sub>2</sub> emission surveys that have been carried out regularly since 2001. Soil gas samples were collected in about 600 sampling sites selected to obtain a homogeneous distribution at about 40 cm depth using a metallic probe and 60 cc hypodermic syringes and stored in 10 cc glass vials. H<sub>2</sub> content was analysed later by a VARIAN CP4900 micro-GC. A simple diffusive emission mechanism was applied to compute the emission rate of H<sub>2</sub> at each survey. Diffuse H<sub>2</sub> emission values were used to construct spatial distribution maps by using sequential Gaussian simulation (sGs) algorithm, allowing the estimation of the emission rate from the volcano. Between 2001-2003, the average diffuse H<sub>2</sub> emission rate was &#8764;2.5 kg&#183;d<sup>&#8722;1</sup> and an increase of this value was observed between 2013-2017 (&#8764;16.6 kg&#183;d<sup>&#8722;1</sup>), reaching a value of 36 kg&#183;d<sup>&#8722;1</sup> on June 2017, 4 month before the first recent seismic swarm in October, 2017 at Cumbre Vieja volcano. Six additional seismic swarms had occurred at Cumbre Vieja volcano (February 2018, July-August 2020; October 8-10, 2020; October 17-19, 2020, November 21, 2020 and December 23-26, 2020) and changes of diffuse H<sub>2</sub> emission related to this unrest had been observed reaching values up to &#8764;70 kg&#183;d<sup>&#8722;1</sup>. Diffuse H<sub>2</sub> emission surveys have demonstrated to be sensitive and excellent precursors of magmatic processes occurring at depth in Cumbre Vieja. Periodic diffuse H<sub>2</sub> emission surveys provide valuable information to improve and optimize the detection of early warning signals of volcanic unrest at Cumbre Vieja volcano.</p>