Seismic liquefaction assessment at different geographical scales provides hazard maps at increasing levels of resolution and reliability. According to the considered level, the areas prone to liquefaction are identified based on specific predisposing and triggering factors, including geological and geotechnical subsoil properties and local seismicity. The robustness of the definition of these areas is strongly correlated to the availability and spatial distribution of this information. Moreover, the data type and quality considerably influence the method of analysis and the degree of uncertainty of the results. This work proposes a methodological approach to define maps in terms of liquefaction susceptibility at two different geographical scales (i.e., regional and sub-regional scale). The Calabria region in Southern Italy was studied for the scope, characterised by a relatively high level of seismicity and liquefaction occurrences recorded during the numerous historical earthquakes. At the regional scale, the level of zonation represents a preliminary assessment of areas potentially susceptible to liquefaction. The considered predisposing factors and geospatial variables were implemented from approximately 5520 available surveys, using geostatistical tools to filter inconsistent data and quantify uncertainties. A quantitative assessment of the liquefaction predisposition was also proposed at a sub-regional scale in the identified areas prone to liquefaction. A new proposal of a Territorial Liquefaction Predisposing Index was carried out by quantifying the uncertainty due to the limited availability of data with geostatistical methods. This approach could provide helpful information for land management and emergency planning to characterise the territory for liquefaction instability and optimise resources allocation to reduce the level of uncertainty at higher levels of analysis.