Data and code from: Mapping italy's sacred natural sites to support the establishment of other effective area-based conservation measures (OECMS)
Data files
Jun 18, 2026 version files 5.65 MB
-
0.1_Compute_road_function.R
4.49 KB
-
0.1_Road_distances.R
1.18 KB
-
0.2_Extract_elev_PAs.R
4.32 KB
-
0.3_CorineBuffers_FigS3.R
13.32 KB
-
0.3_Extract_corine_function.R
1.63 KB
-
0.4_Extract_bioclim.R
1.39 KB
-
0.5_Join_data.R
1.15 KB
-
Figure_1.R
6.56 KB
-
Figure_3.R
7.66 KB
-
Figure2_FigureS1-2.R
10.12 KB
-
README.md
11.86 KB
-
SNS_all_data.csv
1.24 MB
-
sns_all_data.gpkg
1.36 MB
-
SNS_bioclim.gpkg
909.31 KB
-
SNS_corine.gpkg
692.22 KB
-
SNS_dem_distPA.gpkg
622.59 KB
-
SNS_raw.csv
226.64 KB
-
SNS_roads.gpkg
532.48 KB
Abstract
Aim
Sacred natural sites (SNS) are culturally governed places that often retain remnants of natural vegetation, safeguard species of conservation concern, and maintain landscape continuity, positioning them as potential Other Effective Area-Based Conservation Measures (OECMs). Despite increasing ecological interest, national datasets remain unavailable, limiting reproducibility of studies and preventing systematic evaluation of SNS as OECMs. Here we release the first harmonized, openly accessible, geospatial inventory of SNS in Italy, comprising 2,332 records. The dataset provides detailed information on the types of SNS, their ecological attributes, and their landscape context, offering a valuable resource for biodiversity conservation as well as research in biogeography, landscape ecology, and biocultural science.
Location
Italy.
Feature
SNS of Catholic tradition (e.g., shrines, hermitages, monasteries, chapels, pilgrimage sites, groves, springs).
Results
The dataset contains 2,332 georeferenced SNS distributed across Italy, mainly concentrated in central and northern regions. Only a small proportion of SNS occur within protected areas (17%), with the majority corresponding to sites lacking an assigned IUCN category (38%). SNS are predominantly situated at low to mid elevations (mean = 458 m, SD = 343 m) and in landscapes dominated by agricultural (45%) and natural land covers (38%). Their climatic space overlaps with, but extends beyond, that encompassed by the protected area network, indicating potential complementarity in conservation representation.
Main Conclusions
This dataset provides a foundational resource for assessing SNS as potential OECMs, quantifying their ecological contexts, and integrating biocultural heritage into conservation planning, including pathways toward 30 x 30 implementation.
Description of the data and file structure
0.1_Compute_road_function.R
Function to compute road distance and category for SNS points.
0.1_Road_distances.R
Script to calculate distances from SNS points to nearest roads and classify road types.
0.2_Extract_elev_PAs.R
Script to extract elevation data and protected area status for SNS points; and protected areas distances for unprotected SNS.
0.3_Extract_corine_function.R
Function to extract land cover areas from CORINE dataset within specified buffer distances around SNS points.
0.3_CorineBuffers_FigS3.R
Script to calculate land cover areas from CORINE dataset within multiple buffer distances around SNS points and create Figure S3 of the Supplementary Material.
0.4_Extract_bioclim.R
Script to extract bioclimatic variables (BIO1-19) for SNS points from CHELSA dataset.
0.5_Join_data.R
Script to join all the output of the different script and create the final dataset in .csv and .gpkg format
Figure_1.R
Script to create Figure 1 of the manuscript
Figure2_FigureS1-2.R
Script to create Figure 2, S1 and S2 of the manuscript and the Supplementary Material
Figure_3.R
Script to create Figure 3 of the manuscript
SNS_roads.gpkg
Geopackage file containing SNS points with computed road distances and categories (generated from 0.1_Road_distances.R).
SNS_dem_distPA.gpkg
Geopackage file containing SNS points with extracted elevation, information about protected areas and protected area distances (generated from 0.2_Extract_elev_PAs.R).
SNS_corine.gpkg
Geopackage file containing SNS points with extracted CORINE land cover areas within multiple buffer distances (generated from 0.3_CorineBuffers_FigS3.R).
SNS_bioclim.gpkg
Geopackage file containing SNS points with extracted bioclimatic variables (BIO1-19) (generated from 0.4_Extract_bioclim.R).
sns_all_data.gpkg
Geopackage file containing the final joined dataset with all variables for SNS points (generated from 0.5_Join_data.R). See SNS_all_data.csv for column descriptions.
SNS_raw.csv
Description
This file contains the raw data for 2,332 SNS across Italy, including their spatial coordinates and basic attributes.
Columns
The dataset contains the following columns:
| Column name | Description |
|---|---|
prin |
Unique identifier for each SNS unit |
name |
Name of each SNS site |
region |
Region where the SNS is located |
district |
District where the SNS is located |
municipality |
Municipality where the SNS is located |
Ycoord |
Latitude coordinate of the SNS |
Xcoord |
Longitude coordinate of the SNS |
SNS_all_data.csv
Description
This file contains tabular data used for 2,332 SNS across Italy.
The dataset includes spatial attributes, environmental variables, protection status,
and proximity metrics related to protected areas, and is intended for use in spatial and ecological analyses.
The README documents only the column structure, not interpretations or derived meaning.
Columns
The dataset contains the following columns:
| Column name | Description |
|---|---|
prin |
Unique identifier for each SNS unit |
name |
Name of each SNS site |
region |
Region where the SNS is located |
district |
District where the SNS is located |
municipality |
Municipality where the SNS is located |
distance_to_road_m |
Distance in meters to the nearest road |
nearest_road_type |
Category of the nearest road |
road_category |
Macro-level classification of the nearest road |
X50m_Agricultural.areas |
Area (m²) of the specified CORINE Land Cover class within a 50 m buffer around the SNS point |
X100m_Agricultural.areas |
Area (m²) of the specified CORINE Land Cover class within a 100 m buffer around the SNS point |
X150m_Agricultural.areas |
Area (m²) of the specified CORINE Land Cover class within a 150 m buffer around the SNS point |
X200m_Agricultural.areas |
Area (m²) of the specified CORINE Land Cover class within a 200 m buffer around the SNS point |
X250m_Agricultural.areas |
Area (m²) of the specified CORINE Land Cover class within a 250 m buffer around the SNS point |
X50m_Natural.areas |
Area (m²) of the specified CORINE Land Cover class within a 50 m buffer around the SNS point |
X100m_Natural.areas |
Area (m²) of the specified CORINE Land Cover class within a 100 m buffer around the SNS point |
X150m_Natural.areas |
Area (m²) of the specified CORINE Land Cover class within a 150 m buffer around the SNS point |
X200m_Natural.areas |
Area (m²) of the specified CORINE Land Cover class within a 200 m buffer around the SNS point |
X250m_Natural.areas |
Area (m²) of the specified CORINE Land Cover class within a 250 m buffer around the SNS point |
X50m_Urban.surfaces |
Area (m²) of the specified CORINE Land Cover class within a 50 m buffer around the SNS point |
X100m_Urban.surfaces |
Area (m²) of the specified CORINE Land Cover class within a 100 m buffer around the SNS point |
X150m_Urban.surfaces |
Area (m²) of the specified CORINE Land Cover class within a 150 m buffer around the SNS point |
X200m_Urban.surfaces |
Area (m²) of the specified CORINE Land Cover class within a 200 m buffer around the SNS point |
X250m_Urban.surfaces |
Area (m²) of the specified CORINE Land Cover class within a 250 m buffer around the SNS point |
X50m_Water.bodies |
Area (m²) of the specified CORINE Land Cover class within a 50 m buffer around the SNS point |
X100m_Water.bodies |
Area (m²) of the specified CORINE Land Cover class within a 100 m buffer around the SNS point |
X150m_Water.bodies |
Area (m²) of the specified CORINE Land Cover class within a 150 m buffer around the SNS point |
X200m_Water.bodies |
Area (m²) of the specified CORINE Land Cover class within a 200 m buffer around the SNS point |
X250m_Water.bodies |
Area (m²) of the specified CORINE Land Cover class within a 250 m buffer around the SNS point |
X50m_Wetlands |
Area (m²) of the specified CORINE Land Cover class within a 50 m buffer around the SNS point |
X100m_Wetlands |
Area (m²) of the specified CORINE Land Cover class within a 100 m buffer around the SNS point |
X150m_Wetlands |
Area (m²) of the specified CORINE Land Cover class within a 150 m buffer around the SNS point |
X200m_Wetlands |
Area (m²) of the specified CORINE Land Cover class within a 200 m buffer around the SNS point |
X250m_Wetlands |
Area (m²) of the specified CORINE Land Cover class within a 250 m buffer around the SNS point |
bio1 |
Annual Mean Temperature (°C × 10) |
bio2 |
Mean Diurnal Range (Mean of monthly (max temp − min temp)) (°C × 10) |
bio3 |
Isothermality (BIO2 / BIO7 × 100) |
bio4 |
Temperature Seasonality (standard deviation × 100) |
bio5 |
Max Temperature of Warmest Month (°C × 10) |
bio6 |
Min Temperature of Coldest Month (°C × 10) |
bio7 |
Temperature Annual Range (BIO5 − BIO6) (°C × 10) |
bio8 |
Mean Temperature of Wettest Quarter (°C × 10) |
bio9 |
Mean Temperature of Driest Quarter (°C × 10) |
bio10 |
Mean Temperature of Warmest Quarter (°C × 10) |
bio11 |
Mean Temperature of Coldest Quarter (°C × 10) |
bio12 |
Annual Precipitation (mm) |
bio13 |
Precipitation of Wettest Month (mm) |
bio14 |
Precipitation of Driest Month (mm) |
bio15 |
Precipitation Seasonality (Coefficient of Variation, %) |
bio16 |
Precipitation of Wettest Quarter (mm) |
bio17 |
Precipitation of Driest Quarter (mm) |
bio18 |
Precipitation of Warmest Quarter (mm) |
bio19 |
Precipitation of Coldest Quarter (mm) |
elevation |
Elevation of each SNS (m) |
PAs |
Protection status of the site (inside or outside a protected area) |
IUCN_Cat |
IUCN management category of the protected area, if applicable |
nearest_PA_km |
Distance (in kilometers) from an unprotected SNS to the nearest protected area; NA for protected sites |
nearest_PA_IUCN |
IUCN management category of the nearest protected area |
dist_class |
Categorical distance class to the nearest protected area (in km) |
lon |
Longitude coordinate of the SNS |
lat |
Latitude coordinate of the SNS |
The dataset derives from a standardized nationwide identification protocol combining documentary sources, remote-sensing geolocation, and land cover filtering to identify sites embedded within natural or semi-natural landscapes. We provide a georeferenced point vector dataset, with each point representing the location of a SNS. For each site we report geographic coordinates, elevation, bioclimatic variables, proximity to road infrastructure, surrounding land cover composition, and indicators of formal protection status.
