The Large Volume Detector (LVD), hosted in the INFN Laboratori Nazionali del Gran Sasso, is triggered by atmospheric muons at a rate of $\sim 0.1$~Hz. The data collected over almost a quarter of century are used to study the muon intensity underground. The 50-million muon series, the longest ever exploited by an underground instrument, allows for the accurate long-term monitoring of the muon intensity underground. This is relevant as a study of the background in the Gran Sasso Laboratory, which hosts a variety of long-duration, low-background detectors. We describe the procedure to select muon-like events as well as the method used to compute the exposure. We report the value of the average muon flux measured from 1994 to 2017: $\mathrm{I_{\mu}^0 = 3.35 \pm 0.0005^{stat}\pm 0.03^{sys} \cdot 10^{-4} ~m^{-2} s^{-1}}$. We show that the intensity is modulated around this average value due to temperature variations in the stratosphere. We quantify such a correlation by using temperature data from the European Center for Medium-range Weather Forecasts: we find an effective temperature coefficient $\mathrm{\alpha_{T}} = 0.94\pm0.01^{stat} \pm0.01^{sys}$, in agreement with other measurements at the same depth. We scrutinise the spectral content of the time series of the muon intensity by means of the Lomb-Scargle analysis. This yields the evidence of a 1-year periodicity, as well as the indication of others, both shorter and longer, suggesting that the series is not a pure sinusoidal wave. Consequently, and for the first time, we characterise the observed modulation in terms of amplitude and position of maximum and minimum on a year-by-year basis.

In this work, we apply multichannel singular spectrum analysis (MSSA), a data-adaptive, multivariate, non-parametric technique that simultaneously exploits the spatial and temporal correlations of the input data to extract common modes of variability, to investigate the intermediate quasi-periodicities of the Fe XIV green coronal emission line at 530.3 nm for the period between 1944 and 2008. Our analysis reveals several significant mid-term periodicities in a range from about one to four years that are consistent with the so-called quasi-biennial oscillations (QBOs), which have been detected by several authors using different data sets and analysis methods. These QBOs display amplitudes varying significantly with time and latitude over the six solar cycles (18 to 23) covered by this study. A clear North-South asymmetry is detected both in their intensity and period distribution, with a net predominance of spectral power in the active-region belt of the northern hemisphere. On the other hand, while the QBOs with periods > 1.7 years are particularly intense around the polar regions and therefore related to the global magnetic field, the ones with shorter periods are mainly generated at mid-latitudes, in correspondence with the emergence of active regions. Our findings indicate that the North-South asymmetry manifested in the uneven latitudinal distribution of QBOs is a fundamental, albeit puzzling, characteristic of solar activity.

Typical reconstructions of historic heliospheric magnetic field (HMF) $B_{\rm HMF}$ are based on the analysis of the sunspot activity, geomagnetic data or on measurement of cosmogenic isotopes stored in terrestrial reservoirs like trees ($^{14}$C) and ice cores ($^{10}$Be). The various reconstructions of $B_{\rm HMF}$ are however discordant both in strength and trend. Cosmogenic isotopes, which are produced by galactic cosmic rays (GCRs) impacting on meteoroids and whose production rate is modulated by the varying HMF convected outward by the solar wind, may offer an alternative tool for the investigation of the HMF in the past centuries. In this work, we aim to evaluate the long-term evolution of $B_{\rm HMF}$ over a period covering the past twenty-two solar cycles by using measurements of the cosmogenic $^{44}$Ti activity ($\tau_{1/2} = 59.2 \pm 0.6$ yr) measured in 20 meteorites which fell between 1766 and 2001. Within the given uncertainties, our result is compatible with a HMF increase from $4.87^{+0.24}_{-0.30}$ nT in 1766 to $6.83^{+0.13}_{-0.11}$ nT in 2001, thus implying an overall average increment of $1.96^{+0.43}_{-0.35}$ nT over 235 years since 1766 reflecting the modern Grand maximum. The $B_{\rm HMF}$ trend thus obtained is then compared with the most recent reconstructions of the near-Earth heliospheric magnetic field strength based on geomagnetic, sunspot number and cosmogenic isotope data.