This database includes a large collection of quality-controlled and homogenized historical snow records measured in the 1951-2001 period in the Central and Southern Apennine Mountains (Italy). Such data have been manually digitized from the Hydrological Yearbooks of the Italian National Hydrological and Mareographic Service (hereafter, NHMS), the institution that managed the hydro-meteorological data collection in Italy from 1917 to 2002. More specifically, the rescued dataset includes the monthly observations of three different variables: · The snow cover duration (SCD), which is defined as total number of days in a given month with snow depth on the ground >=1 cm. This variable is available for 110 stations between 288 and 1430 m above the sea level (ASL). · The number of days with snowfall (NDS), which is total number of days in a given month on which the accumulated snowfall (i.e. the amount of fresh snow with respect to the previous observations) is at least 1 cm. This variable is available for 114 stations between 288 and 1430 m ASL. · The height of new snow (HN), which is defined as the monthly amount of fresh snow (expressed in cm). The monthly value is intended as the sum of daily HN data observed in a determined month. This variable is available for 120 stations between 288 and 1750 m ASL. Note that for HN variable, the data availability is restricted to the period 1971-2001. The considered dataset has been subjected to an accurate quality control consisting of several statistical tests: the gross error test, which flags the data that are above or below acceptable physical limits, the consistency test, which involves an inter-variable check, and the tolerance test, which is focused on the outlier detection. In addition, the homogeneity of the rescued time series has been checked using Climatol method (Guijarro, 2018). The latter is based on the Standard Normal Homogeneity Test (Alexandersson, 1986) for the identification of the breaks and on a linear regression approach for the adjustments (Easterling and Peterson, 1995). Climatol has been also employed for the filling of missing values. The database is structured into three different folders (one for each variable). In a determined folder, the user finds two files, one containing the main information regarding the available stations (code, station name, latitude and longitude (in decimal degrees) and altitude ASL (in m)), the other one the monthly time series for the considered variable. Note that the original data sources of this database, the Hydrological Yearbooks of the NHMS, are freely accessible in printed version (i.e. as scanned images in portable document format) through the Italian Institute for Environmental Protection and Research (ISPRA) website (http://www.bio.isprambiente.it/annalipdf). Additional information about the data rescue processing can be found in the preprint “Historical snowfall measurements in the Central and Southern Apennine Mountains: climatology, variability and trend”, open for discussion in The Cryosphere journal (https://doi.org/10.5194/egusphere-2024-1056). References Alexandersson, H.: A homogeneity test applied to precipitation data, J. Climatol., 6, 661–675, 1986. Easterling, D. R. and Peterson, T.C.: A new method for detecting and adjusting for undocumented discontinuities in climatological time series, International Journal Climatol.,15, 369–377, https://doi.org/10.1002/joc.3370150403, 1995. Guijarro, J. A.: Homogenization of climatic series with Climatol, Climatol manual, https://www.climatol.eu/homog_climatolen.pdf (last access: 15 February 2024), 2018.