Fostering and Planning a Smart Governance Strategy for Evaluating the Urban...

[Full title: Fostering and Planning a Smart Governance Strategy for Evaluating the Urban Polarities of the Sardinian Island (Italy)]

sustainability Article Fostering and Planning a Smart Governance Strategy for Evaluating the Urban Polarities of the Sardinian Island (Italy) Chiara Garau * , Giulia Desogus and Mauro Coni Department of Civil and Environmental Engineering and Architecture, University of Cagliari, 09123 Cagliari, Italy * Correspondence: cgarau@unica.it; Tel.: + 39-070-675-5565 Received: 26 July 2019; Accepted: 6 September 2019; Published: 11 September 2019 Abstract: The interrelations between cities, inland areas, connecting road networks, urban, and political polarities have evolved, thereby determining economic, social, and place-based impacts. Thus, via a case study of Sardinia island (Italy), this study analyses regional transport data to evaluate the interrelations and mobility issues between the main cities and the settlement geographies of internal areas with a predominantly agricultural vocation. First, it frames the problems (common to the islands) theoretically and focuses on how the internal areas (considered marginal for a long time) have considerable material and immaterial resources to be valorised. Second, the study evaluates the internal relationship networks that characterise the island territory through the cluster and principal components analysis using origin–destination data to represent vocations and population needs A smart governance strategy is proposed for Sardinia through an assessment of the functionality of urban settlements and interconnections between the hinterlands (the small and the main cities of the case study), following the smart region paradigm. The study underlines the importance of the interconnection between urban geographical areas. Thus, given an analytical-numerical approach, the originality of this research is highlighted in how it is possible to extract social vocations of the territory, which is generally not easily quantifiable Keywords: smart governance strategy; smart island; smart region; mobility; origin–destination analysis; Sardinia 1. Introduction Literature review [ 1 – 4 ] underlines how the economic, social and territorial cohesion is one of the main objectives of the Treaty on the Functioning of the European Union (TFEU, art. 174). However, the place-based cohesion in the island territories have more structural problems that make it more problematic to identify, understand and manage the entire territory through a smart government. Therefore, this article wants to understand how the analysis of dependencies between municipalities, allows adopting a smart governance di ff erentiated according to the polarization characteristics of each area. To do this the authors start from the literature review and from analysis of urban polarities through the theory of central places (CPT), subsequently this is applied through a methodology that combines the economic, social and territorial aspects of the case study of Sardinia island (Italy) About that, the TFUE aims to reduce the gap between regions, focusing above all on rural areas a ff ected by industrial transition or those that present serious and permanent natural or demographic disadvantages, such as the northernmost regions with very low population density and the island regions [ 5 , 6 ]. Moreover, the European programme 2014 / 2020 [ 7 ] establishes an intervention strategy that can be summarised in three objectives: “(1) Protecting the territory and the safety of the Sustainability 2019 , 11 , 4962; doi:10.3390 / su 11184962 www.mdpi.com / journal / sustainability

Sustainability 2019 , 11 , 4962 2 of 24 inhabitants [ . . ]; (2) Promoting natural and cultural diversity and polycentrism through an economic and social model that is able to create models of life in the internal areas that are competitive with urban areas both in the same region and outside the region [ . . ]; and (3) Relaunching development and work through the use of potential misused resources” [ 7 ]. Hence, in Italy, the Economic and Social Cohesion Department within the EU Framework Programme Horizon 2020 highlights that 60% of the national territory in Italy consists of inland areas in which about a quarter of the Italian population lives [ 8 , 9 ]. In Italy, in 2012, the regional policy planning “methods and objectives for e ff ective use of EU funds 2014–2020” defines the inland areas as unstable agglomeration centres with strong demographic issues. However, at the same time, they are strongly polycentric; they have the potential for attraction and resources that are lacking in the central areas [ 7 ]. In this context, the Technical Committee of the “Inner Areas” and the representatives of the various Italian Regions drafted the document relating to the National Strategy of Internal Areas to adapt primary services to the needs of the population and develop the local economy through local production chains [ 10 ]. This document has, therefore, initiated a new and widespread awareness in Italy regarding the internal areas to develop a complete framework and better plan for them by considering multipolar interconnections Thus, in Italy, the National Strategy of Internal Areas (2012) developed the concept of inland areas, which are linked to regional developments and interrelations with the main cities, by employing a planning process that is applied to local economic and sustainable developments [ 11 ]. In this concept, regional peculiarities must relate to large-scale strategic planning and the internal and main areas In July 2018, the Italian o ffi ce for the evaluation of the Impact of Cohesion Policies in Europe [ 12 ] showed that, in line with the rest of Europe, despite the fact that the Italian structural funds had a significant e ff ect on the GDP per capita [ 11 ], the place-based impact has negative results for individual Italian regions. Thus, the South of Italy, including the major islands (Sicily and Sardinia), “with its twenty million inhabitants, remains the largest less-developed area” ([ 12 ], pp. 10–11) This concept becomes even more significant when it is applied to insular contexts. In fact, they have di ff erent structural disadvantages that not only create di ffi culties for economic and social development [ 13 , 14 ] but also have their own place-based, physical, and managerial characteristics [ 15 – 20 ]. Starting from these assumptions and considering a regional data set to develop new ideas on regional planning [ 21 , 22 ], this study primarily determines a smart territory approach to analyse the phenomenon of the connection between inland and coastal areas by considering a particularly sensitive context, such as the island of Sardinia Sardinia appears to be particularly emblematic as a case study not only because it is an island context but also because, within the European programming 2014–2020 ([ 6 ], pp. 34–35), it belongs (together with Abruzzo and Molise) to the Italian regions with great di ffi culty in economic and social development. Moreover, the Sardinia region participates in the National Strategy for Italian Internal Areas (NSIA, from Italian: Strategia Nazionale per le Aree Interne Italiana—SNAI), where the peculiarities of the territory are seen as a fundamental potential to counteract the phenomenon of depopulation and promote economic and social development [ 23 – 25 ]. Thus, the study begins with an in-depth bibliographic and legislative analysis on the common problems of islands with their implications on the condition of insularity. In support of these reflections, the authors reflect on theories of central places (CPT) applied to mobility as a possible concrete application to island territories Subsequently, given the latest available data on mobility regarding origin–destination, after normalisation and rendering them as indicators, the study arrived at a methodology of territorial analysis, in which a cluster and a principal component analysis (PCA) are presented Then, the articulation of homogeneous areas on which it is necessary to intervene is analysed. That is, the elaborated methodology recognises the dominant attractors of the largest urban areas within an insular system, and it can create and reorganise intermediate aggregations classified as “homogeneous areas” to provide a set of minimum services in weak internal environments, which is today not satisfied.

Sustainability 2019 , 11 , 4962 3 of 24 Finally, the results of the study are presented and discussed by showing how the methodology adopted for the Sardinia region in di ff erent homogeneous areas consents to adopting di ff erentiated policies according to the polarisation characteristics of each area 2. Literature Review on Common Issues Related to Insular Regions From a theoretical point of view, with the Treaty of Maastricht [ 26 ] in 1992, the European Union recognises the islands as ultra-peripheral regions. However, only in 1997, the “Declaration concerning the island regions” annexed to the Treaty of Amsterdam (1997) recognises, for the first time, the various structural disadvantages that hinder economic and social development, which are typical issues of island regions ([ 14 ], p. 136). Subsequently, EU legislation has increasingly recognised the presence of handicaps that lead to development shortages [ 27 ]; it introduced more specific financial and research actions to solve, at least in part, the common island problems (Table 1 ). However, the European Union has primarily shown attention to the island problems in a theoretical way Very little specific measures were produced “even though at the same time the islands have made good use of the instruments devoted to the policy of economic and social cohesion” ([ 27 ], p. 3). This situation may be partly due to the di ffi culty in finding a common classification for islands that serves as a starting point for place-based governance. Monfort, for example, argues [ 28 ] that in every island characterised by its own insularity and by political and socio-economic specificities, “the development is generally a ff ected by a complex set of constraints, which are often not all specific to the type of territory included in that category. [ . . ] The main implication of such observations is that the categories of these specific territories are far from constituting homogenous groups of regions” ([ 28 ], p. 10) The ESPON 2013 programme [ 29 ] instead underlines the intrinsic and extrinsic diversity of the island regions regarding others, thereby supporting the need to diversify the island policies “from EU and national goals for sustainable development, as well as those for economic, social and spatial cohesion. [ . . ] Thus, it becomes imperative to apply appropriate policies, which can enhance the islands’ attractiveness” ([ 29 ], p. 14). The literature also shows a discrepancy not only in the classification but also on the number of islands present in Europe [ 30 , 31 ]. Furthermore, if the instruments for economic and social cohesion and the trans-European networks are the main tools used by the European Union to reduce the island’s disadvantages [ 25 ], the proceedings of the 39 th annual general meeting of the CPMR Island Commission of 2019 [ 32 ] show the legislative failures in matters of island development. However, in the 2016 Treaty on the Functioning of the European Union, articles 174 and 175 pay much attention to insular problems; thus, Member States conducted their economic policy to reduce the gap between the levels of development of the various regions and the delay of the less favoured regions including the island regions [ 5 ] (art. 174–175) In addition, the CPMR’s Islands Commission (IC) notes, not only are “the existing regulatory framework and existing devices and measures insu ffi cient” ([ 33 ], p. 4) but also there is a need for an e ff ective solution to all the problems a ff ecting the European island territories listed on the European Union in the Opinion of the Economic and Social Committee in 2000. Table 2 shows the structural problems and permanent disadvantages common to island regions [ 16 , 34 ], highlighting a series of structural indications for evaluating policies to be implemented in the governance of islands Italy has not yet defined its own legislation on insularity issues. In 2016, with decree 497, the Senate developed a framework law for the development of the smaller islands and Sicily and Sardinia have adopted the “condition of insularity” This legislative revision shows the lack of a smart governance framework without giving general guidelines for island regions, as well as place-based strategies or policies able to solve problems of these territories with their daily habits. Namely, literature shows only specific case studies, by using, for instance, a minor island as a case study, or by giving the general characteristics of what a smart Island should be [ 35 , 36 ], or trying to define what a smart island is [ 16 , 17 , 37 ], or by applying the concept of smart governance in limited portions of territory [ 38 , 39 ]. Other studies associate the smart governance of an island with the use and impact of new information and communication technologies

Sustainability 2019 , 11 , 4962 4 of 24 (ICT) [ 40 ], the e ff ects of e-government development on environmental sustainability [ 41 – 43 ], or on urban sustainability [ 44 ]. Table 1. Summary of the major European initiatives on smart islands Years European Initiatives on Smart Islands 1992 1 Treaty of Maastricht, art. 129 b and “Declaration on the outermost regions of the Community” 1993 2 ISLENET is the first network of island authorities promoting sustainable energy and environmental management 1997 3 Treaty of Amsterdam, art. 158 CE (ex 130 A) and “Declaration on island regions” 2000 Guidelines for integrated actions on the island regions of the European Union (EU O ffi cial Journal C 268, 19 September 2000) 2002 2 The Committee of the Regions presents the problems of island regions in the European Union, emphasising the EU policies and regulation in the relevant areas 2003 4 Analysis of the island regions and outermost regions of the European Union. Part I: The island regions and territories 2007 2 The European Parliament adopts the Motion for Resolution tabled by MEP Francesco Musotto, which proposes tailored measures to some challenges facing European islands 2009 5 Territories with specific geographical features. A series of short papers on regional research and indicators produced by the Directorate-General for Regional Policy 2011 2 The European Parliament o ffi cially recognises the Pact of Islands as an EU initiative that runs parallel to the Covenant of Mayors 2012 2 EURELECTRIC publishes a report demonstrating islands’ potential to become test-beds for cutting-edge sustainable energy technologies 2013 2 Thirteen islands partners across Europe work together in clusters to mature sustainable energy projects through enhanced multi-level governance. The key deliverable of SMILEGOV is the Smart Islands Strategy 2013 6 The ESPON 2013 Programme. The Development of the Islands—European Islands and Cohesion Policy (EUROISLANDS) 2014 2 The Opinion calls for the implementation of smart policies and smart development initiatives at EU, national, and regional levels that cater to the specific characteristics of islands 2015 2 The 35 th Annual Conference of the CPMR Islands Commission marks a strong comeback of islands in the European Parliament. Forty-one new members join the Pact of Islands, which today counts 117 members 2015 2 The European Parliament adopts the Resolution on the special situation of islands calling for an Agenda for EU Islands accompanied by an EU Strategic Framework for Islands to link up financial instruments that can have a major territorial impact 2016 2 In Athens, on 21–22 June, DAFNI hosts the first Smart Islands Forum 2017 2 On 20 February, EURELECTRIC publishes the “Towards the Energy Transition on Europe’s Islands” report, which highlights the role of islands in Europe’s energy transition 2017 2 On 28 March, the islands make a strong comeback at the European Parliament to present the Smart Islands Initiative and communicate their potential to transform into smart, inclusive, and thriving societies for an innovative and sustainable Europe 1 Treaty on European Union. O ffi cial Journal of the European Communities No C 191 / 3 https: // eur-lex.europa.eu / legal-content / IT / TXT / ?uri = CELEX:11992 M / TXT (accessed on 13 July 2019) 2 Smart Island Initiative. The initiative: the story behind http: // www.smartislandsinitiative.eu / en / history.php (accessed on 13 July 2019) 3 Treaty of Amsterdam http: // www.europarl.europa.eu / topics / treaty / pdf / amst-en.pdf (accessed on 13 July 2019) 4 Analysis of the island regions and outermost regions of the European Union https: // ec.europa.eu / regional_policy / sources / docgener / studies / pdf / ilesrup / islands_part 1_summary_en.pdf (accessed on 13 July 2019) 5 Monfort, P. Territories with specific geographical features https: // ec.europa.eu / regional_policy / sources / docgener / work / 2009_02_geographical. pdf (accessed on 13 July 2019) 6 The ESPON 2013 Programme The Development of the Islands—European Islands and Cohesion Policy (EUROISLANDS) https: // www.espon.eu / sites / default / files / attachments / inception_report_full_ version.pdf (accessed on 13 July 2019).

Sustainability 2019 , 11 , 4962 5 of 24 Table 2. Structural problems common to island regions derived from the Opinion of the Economic and Social Committee (O ffi cial Journal of the European Union C 268, 19 September 2000) Permanent Disadvantages and Common Characteristics of Islands Contexts 1 Isolation from land 2 Higher cost of maritime and air transport, communications and infrastructure due to natural and climatic obstacles 3 Limited extension of usable land 4 Poor fishing resources 5 Limited water resources 6 Limited energy sources 7 Marine and coastal pollution 8 Particular di ffi culty in waste management 9 Population decline, particularly of young people 10 Coastal erosion 11 Lack of qualified personnel 12 Absence of an economic environment favourable to entrepreneurship 13 Di ffi cult access to school and health services For these reasons, this paper aims to show how the study of mobility, through the analysis of origin–destination shifts (O / D), can describe the vocations and needs that characterize the territory through the network of relations and attitudes of citizens However, a general governance framework for major islands is still lacking. In fact, these kind of islands by their nature have a net distinction between inland areas and coastal ones, as well as a replicable methodology that can allow not only a complete study of the island territory, starting from the usual habits and actions of population, but also a smart governance able to solve the structural and common problems of many European Islands Before describing this, the authors consider it necessary to highlight, in the next section, that one of the major problems of insularity is the clear di ff erence in planning between inland areas and coastal cities [ 45 , 46 ]; thus, in most cases, the settlements that are on the coastal areas are the main centres of the whole island [ 47 ]. This phenomenon, which the authors explain through the spatial distribution linked to the theory of central places (CPT), has important consequences on the interrelations between major cities and inland areas 3. The Analysis of Urban Polarities to Understand the Interconnection between Inner End Coastal Areas in Insular Regions The problem with inland areas is particularly evident in the literature on an island context [ 45 , 48 , 49 ], especially when (as in Sicily and Sardinia for Italy) larger islands are considered for their development, comparable with non-island regions. Therefore, it was necessary to find a solution to the problem of managing these areas through a common policy that concerns the whole Region. This section explained, through central theories, how the multipolar interconnection can provide a common smart governance to di ff erent urban poles within an insular system Moreover, “the identification of the internal areas [ . . ] starts from a polycentric interpretation of the territory [ . . ] characterised by a network of municipalities or groups of municipalities [ . . ] around which gravitate areas characterised by di ff erent levels of the spatial periphery” ([ 50 ], p. 1) This situation makes it necessary to support this study with the main models of representation of the regional structure as the planning of multipolar interconnections. In fact, the theme of island development, especially in terms of smart planning in a regional geographical area, reflects on the

Sustainability 2019 , 11 , 4962 6 of 24 theory of central places (CPT) through a synthesis of the organisation mechanisms between major cities, suburbs, inland areas and rural areas [ 51 , 52 ]. The first theories on CPT [ 53 – 55 ] studied the territorial rebalancing through multi-level urban hierarchies that established places of great economic and social development [ 53 ]. In the course of its evolution, the concept of central places was applied by di ff erent urban planning models, mainly based on the regional economy [ 56 – 58 ]. Christaller (1933) and Lösch (1941), pioneers of the concept of central place system [ 53 – 56 ] interpret central locations as hierarchical centres at multiple levels that establish the main places on the economic and social level [ 54 ]. These models also allow for evaluating the principles that determine the number, size, and distribution of cities [ 6 , 54 , 59 ], thus, including relationships between the city centre and internal areas and the e ff ects they produce Although they are old models since they were applied in di ff erent periods and contexts from those of today, they can still be valid in an insular system in their theoretical formulation, adopting more current tools and methodologies. This situation is because the characteristics of the geographical and spatial structure of an island have, by their nature, a precise boundary (the system of an island is closed and isolated) that allows the analysis of the urban system within it, stressing the concept of main and secondary hubs. Moreover, if, over the years, the literature changed the forms of use of central place theory [ 60 – 66 ] and their application to specific contexts [ 67 – 71 ], the concept of central localization is still applied to urban planning, correlating the spatial distribution of settlements with, for example, public structures [ 72 – 74 ], industries [ 75 , 76 ], and the tourism sector [ 77 – 82 ]. Furthermore, the study correlated the optimal localisation strategy with the economic demand of defined geographical areas Thus, Van Meeteren and Poorthuis [ 68 ] propose the study of settlement geographies through the analysis of social media by adopting as the basis of research the theory of central places to demonstrate how, today, through the data collected by computer analysis, it is possible to shape the geographical economic-urban areas [ 68 ]. The model was also applied to insular contexts, mainly using an economic approach to enhance the area through the location of the main tourist polarities [ 83 – 86 ]. In regional planning, these models are integrated with the analysis of mobility flows [ 87 – 90 ] as a fundamental component to understanding the movements and the population needs. Lenormand and Ramasco [ 91 ] argue on how technology can help regarding the analysis of urban systems and human mobility to better understand urban dynamics. Mul í ˇcek and Mal ý [ 89 ] argue that changes in commuting workflows show the development of the urban system hierarchy. Similarly, Mitchell [ 90 ] notes that the settlements are organised in networks of shifts, identifiable around the major hub. In fact, the importance of a representation of the dynamics of population displacement in urban hubs at di ff erent levels allows for understanding not only a local classification of the centres but also the e ff ects that strong hubs can produce in the weaker urban centres with the aim of a regional and place-based rebalancing [ 92 ]. The analysis of the urban polarities of an island acquires even more value if analysed from the point of view of mobility as the impossibility of expanding beyond island boundaries makes the analysis of population flows e ff ective in understanding the main hubs. This situation creates the opportunity to evaluate the functionality of urban settlements and their interconnections within the system and, therefore, reflects a smart governance strategy [ 93 ] capable of modifying normal flows through the analysis of the needs and of the vocations of the same population In this regard, the study uses a methodology applied to the insular context of Sardinia in Italy, to redefine urban strategic polarities using the theory of central places (CPT) applied to the mobility sector through the principal component analysis and the cluster analysis. This methodology not only demonstrated the relevance and practicality of the CPT in the study of urban strategic polarities in the Islands context but allows for responding to the three problems indicated by the National Strategy of the Italian Areas, which are key points of the regional island development (1) The polycentric articulation of the territory describes Italy in two historical subsystems (large, medium, small urban systems; and internal areas). They had divergent economic development

Sustainability 2019 , 11 , 4962 7 of 24 and industrial take-o ff , causing “the stagnation, degrowth and development of the Inner Areas” ([ 10 ], p.14) (2) The internal areas have an unused territorial capital, where the economic decrease regards the social unease, linked to the lack of basic services and the social costs generated by the production and consumption processes ([ 10 ], pp. 14–15) (3) The demographic decline in the population, with a strong decrease in births and an increase in the elderly population, remains a strategic objective to intervene to develop inland areas, reversing today’s demographic trend ([ 10 ], p. 20) However, before going into detail and describing the methodology, the authors present the case study of Sardinia, an island in the middle of the Mediterranean 4. The Case Study of the Island of Sardinia in Italy Sardinia is the largest island in Italy after Sicily (Figure 1 ). It covers 24,100.02 square kilometres and has a population of 1,639,591 inhabitants [ 94 ] with a density of 68 inhabitants per sq km, unlike Sicily which extends for 25,832.39 square kilometres with a density of 194 inhabitants per sq km [ 95 ]. Sardinia has 377 municipalities [ 95 ], which is distributed according to the five administrative territorial divisions [ 96 ]: 17 to the metropolitan city of Cagliari; 74 to the province of Nuoro; 87 to the province of Oristano; 92 to the province of Sassari; and 107 to the province of Southern Sardinia [ 97 ]. Sustainability 2019 , 11 , x FOR PEER REVIEW 7 of 24 (3) The demographic decline in the population, with a strong decrease in births and an increase in the elderly population, remains a strategic objective to intervene to develop inland areas, reversing today’s demographic trend ([10], p. 20). However, before going into detail and describing the methodology, the authors present the case study of Sardinia, an island in the middle of the Mediterranean. 4. The Case Study of the Island of Sardinia in Italy Sardinia is the largest island in Italy after Sicily (Figure 1). It covers 24,100.02 square kilometres and has a population of 1,639,591 inhabitants [95] with a density of 68 inhabitants per sq km, unlike Sicily which extends for 25,832.39 square kilometres with a density of 194 inhabitants per sq km [96]. Sardinia has 377 municipalities [96], which is distributed according to the five administrative territorial divisions [97]: 17 to the metropolitan city of Cagliari; 74 to the province of Nuoro; 87 to the province of Oristano; 92 to the province of Sassari; and 107 to the province of Southern Sardinia [98]. The morphological and geographical characteristics, which are very different between coastal landscapes and the hinterland, have, over time, created different degrees of polarity. The coastal areas, which are an unquestionable heritage for the island, are the main polarities as they constitute the greatest potential for tourism and economic development. The Sardinian hinterland is less known by mass tourism than the coastal landscapes. Here, the economy suffers due to lack of work and services, as well as the depopulation of municipalities, thus, becoming secondary polarities that need connections with the coastal area to survive. In these areas, the weakness of the general system leads to considering the need to redesign the levels of attraction by acting on the reduction of circulation impedances and medium-sized services (hospitals, school facilities, and shopping centres). Furthermore, potential polarities can be distinguished by considering those that, due to their proximity to coastal areas and average infrastructure connections, can maintain an economy and good social development. Figure 1. Sardinia island with its morphology and its five administrative territorial divisions. Figure 1. Sardinia island with its morphology and its five administrative territorial divisions The morphological and geographical characteristics, which are very di ff erent between coastal landscapes and the hinterland, have, over time, created di ff erent degrees of polarity. The coastal areas, which are an unquestionable heritage for the island, are the main polarities as they constitute the greatest potential for tourism and economic development. The Sardinian hinterland is less known by mass tourism than the coastal landscapes. Here, the economy su ff ers due to lack of work

Sustainability 2019 , 11 , 4962 8 of 24 and services, as well as the depopulation of municipalities, thus, becoming secondary polarities that need connections with the coastal area to survive. In these areas, the weakness of the general system leads to considering the need to redesign the levels of attraction by acting on the reduction of circulation impedances and medium-sized services (hospitals, school facilities, and shopping centres). Furthermore, potential polarities can be distinguished by considering those that, due to their proximity to coastal areas and average infrastructure connections, can maintain an economy and good social development Therefore, the case study of Sardinia is emblematic for several reasons. The economic development and the social and demographic context of the island show “some signs of fragility: the low birth rate, combined with the low incoming mobility, determines a spiral of natural decrease in the population [ . . ] and the ageing process of the population” [ 98 ]. In addition to the demographic and economic decline, the Sardinian territorial system is particularly interesting as a closed system that allows an unambiguous and stable reading of the flows. It is important to study the origin-destination shifts to understand how the population moves within the island Thus, supported by the aforementioned theoretical models, the needs of the population are first identified, and then the main centres of the island are identified Therefore, given the premises of the European Community and of the National Internal Technical Committee concerning the problems of economic, social, and territorial cohesion and the issues highlighted by the analysis of central place theory and structural problems typical of the region’s insular context, an analysis of the origin-destination (O / D) shifts is proposed to describe the vocations and needs that characterise the territory through the network of relationships and attitudes of citizens 5. Methods Sardinia, given the related problems common to the islands previously described, finds a better declination in the analysis of urban polarities through the theory of central places. The CPT allows for reflecting on which connection mechanisms are useful in a closed system where the coastal area is di ff erent (in terms of economy, morphology and demography) from the internal area. The main goal is to understand the strongest urban polarities of the island system through the interactions between geographical areas analysed through the analysis of mobility data The analysis of mobility data allows for describing the network of relationships and attitudes of the citizens through the analysis of the origin–destination (O / D) shifts. They represent the result of the vocations and needs that characterise a subsystem of the entire territory. Vocations such as identity, demographic structure, history, culture, and the needs of the citizens regarding education, work, health, services, and entertainment are bound by the physical characteristics, such as the transport networks, orography, climate. The movements of citizens express their vocations and needs, and they emerge in the statistical data recorded in the O / D matrices. The analysis of the spatial distribution of the matrices can, therefore, highlight structures and dependencies between di ff erent parts of the territory The data on which the analysis is based is the recurrent daily trips of entering and leaving each municipality divided by the destination municipality. Additional information is the average travel time, the means used, and the reason for moving (study or work) The data do not, therefore, represent the universe of mobility because an important percentage of movements is subtracted from monitoring. Repetitive trips can take place for other reasons (such as health needs, leisure, shopping, tourism) in destinations and modes. Several studies used PCA and CA analysis for land-use and urban plan [ 99 – 103 ]. These techniques have been frequently used to study socio-economic parameters including social, environmental, and urban evolution and as a numerical base for land-use and urban planning. However, they have strengths and weaknesses: the strong point is the robust mathematical objective method and their weakness is the same mathematical objective method. In fact, it fails to fully consider many social aspects of planning that cannot be quantified.

Sustainability 2019 , 11 , 4962 9 of 24 While it is di ffi cult to define the part related to tourism and leisure, this rate, estimated at around 30% of the total, can be, in part, related to the previous one as far as the territorial distribution is concerned From a methodological point of view, the authors conducted an analysis in main components and later, a cluster analysis of the origin-destination (O / D) mobility data surveyed by ISTAT (Italian Institute of Statistics) of 2015 (the latest data available) The study illustrates a series of statistical calculations and territorial contouring, where di ff erent significant indexes of the characteristics of mobility and functional dependencies between the di ff erent municipalities are mapped An important consideration for the analysis developed in the study consists of observing how the depopulation phenomena involve the mountain areas to a greater extent than the coastal and flat areas. The abandonment of mountain areas is a well-known phenomenon and is well correlated with the altitude of the centre considered. Firstly, cluster analysis is conducted to group statistical units by minimising the “logical distance” within each group and maximising the external one between groups [ 104 ]. The “logical distance” is evaluated by parameters of similarity-dissimilarity between the statistical units The k-means function allows to aggregate data into k mutually exclusive clusters and returns the clustering index for each data. Each cluster is defined by its objects and by its centroid, the point to which the sum of distances from all objects is minimised. The “distance” can be evaluated di ff erently in the O / D data analysis algorithm for cluster centre initialisation and the squared Euclidean metric to determine distances (Figure 2 ). The cluster analysis (CA) allows for identifying a minimum number of groups such that the elements belonging to a group are, in some sense, more similar to each other than the elements of other groups. In summary, it is possible to search for groups of objects such that the elements that belong to them are “similar” to each other and di ff erent from those of other groups The principal component analysis (PCA) identifies the underlying structure of the observed relationships. A limited number of new variables, called main components, resulting from the linear combination of the starting variables, allow for a convenient description of the system, eliminating redundancies in the data [ 105 , 106 ]. Sustainability 2019 , 11 , x FOR PEER REVIEW 9 of 24 The study illustrates a series of statistical calculations and territorial contouring, where different significant indexes of the characteristics of mobility and functional dependencies between the different municipalities are mapped. An important consideration for the analysis developed in the study consists of observing how the depopulation phenomena involve the mountain areas to a greater extent than the coastal and flat areas. The abandonment of mountain areas is a well-known phenomenon and is well correlated with the altitude of the centre considered. Firstly, cluster analysis is conducted to group statistical units by minimising the “logical distance” within each group and maximising the external one between groups [105]. The “logical distance” is evaluated by parameters of similarity-dissimilarity between the statistical units. The k-means function allows to aggregate data into k mutually exclusive clusters and returns the clustering index for each data. Each cluster is defined by its objects and by its centroid, the point to which the sum of distances from all objects is minimised. The “distance” can be evaluated differently in the O/D data analysis algorithm for cluster centre initialisation and the squared Euclidean metric to determine distances (Figure 2). The cluster analysis (CA) allows for identifying a minimum number of groups such that the elements belonging to a group are, in some sense, more similar to each other than the elements of other groups. In summary, it is possible to search for groups of objects such that the elements that belong to them are “similar” to each other and different from those of other groups. The principal component analysis (PCA) identifies the underlying structure of the observed relationships. A limited number of new variables, called main components, resulting from the linear combination of the starting variables, allow for a convenient description of the system, eliminating redundancies in the data [106,107]. In the multivariate analysis, it is difficult to visualise data when there are more than five parameters. Thus, using colours and marker size, it is possible to describe the fourth and fifth variables in the xyz Cartesian system. Thus, when many variables are involved in the investigation, their relationships cannot be extracted by analysing the graphs. Fortunately, it is very often observed that groups of variables move together because of the same driving force governing their behaviour. In this case, it is possible to simplify the problem by replacing a group of parameters with a single new variable. ( a ) ( b ) Figure 2. Evaluation of the “logical distance” in cluster analysis: ( a ) distribution of units; ( b ) example of grouping statistical units by minimising the “logical distance” within each group. As an example, Figure 3 shows how the general behaviour of the winds can be approximated by the five vectors on the right rather than by the thousands of single values on the left. Figure 2. Evaluation of the “logical distance” in cluster analysis: ( a ) distribution of units; ( b ) example of grouping statistical units by minimising the “logical distance” within each group In the multivariate analysis, it is di ffi cult to visualise data when there are more than five parameters Thus, using colours and marker size, it is possible to describe the fourth and fifth variables in the xyz Cartesian system. Thus, when many variables are involved in the investigation, their relationships

Sustainability 2019 , 11 , 4962 10 of 24 cannot be extracted by analysing the graphs. Fortunately, it is very often observed that groups of variables move together because of the same driving force governing their behaviour. In this case, it is possible to simplify the problem by replacing a group of parameters with a single new variable As an example, Figure 3 shows how the general behaviour of the winds can be approximated by the five vectors on the right rather than by the thousands of single values on the left Sustainability 2019 , 11 , x FOR PEER REVIEW 10 of 24 ( a ) ( b ) Figure 3. Example of the principal component analysis (PCA): ( a ) general behaviour of the winds by the thousands of single values; ( b ) general behaviour of the winds approximated by the five vectors. The PCA is a quantitatively rigorous technique for this simplification, substituting old variables in a new set of parameters called principal components (PA). Each PA is a linear combination of the old variables. Among an infinite number of combinations, the first PA has the maximum variance of old variables, and it is described as a single axis in space. The second PA is another axis in space, orthogonal to the first. The new combination has the maximum variance on this second axis. It is possible to iterate the procedure until the same number of initial variables explains 100% of the original variance. However, a few PAs can explain more than 80–90% of the original variance. Thus, there is clarity on how the outgoing movements from each municipality can be mainly traced back to the presence of dominant attractors: the main urban areas. The combination of the theory of central places (CPT), CA, and PCA, which is applied to the O/D matrix, represents the original contribution of the study, as the literature does not report similar examples. Finally, the physical characteristics of the transport networks are considered through the isochrones; namely, the time needed to reach any other point in the transport network from a centroid. For this analysis, a network graph is used where each node (n) represents a point of intersection of the branches of the transport network, and the arcs (a) represents the physical connection between them. Each node is characterised by its geographical coordinates and by the time required for its crossing. The arcs are defined not only by the nodes that also connect to the travelling speeds but also the flows that interest them. In this way, it is possible to determine the crossing time of each arc. On the graph so defined, the Dijkstra’s algorithm is applied [108]. Dijkstra’s algorithm is a procedure to find the shortest paths between nodes in a graph. The original algorithm has been applied to obtain the path between any two nodes and the total travel time, adding the time of each link. The process can be iterated to calculate the total time between one Figure 3. Example of the principal component analysis (PCA): ( a ) general behaviour of the winds by the thousands of single values; ( b ) general behaviour of the winds approximated by the five vectors The PCA is a quantitatively rigorous technique for this simplification, substituting old variables in a new set of parameters called principal components (PA). Each PA is a linear combination of the old variables Among an infinite number of combinations, the first PA has the maximum variance of old variables, and it is described as a single axis in space. The second PA is another axis in space, orthogonal to the first. The new combination has the maximum variance on this second axis It is possible to iterate the procedure until the same number of initial variables explains 100% of the original variance. However, a few PAs can explain more than 80–90% of the original variance Thus, there is clarity on how the outgoing movements from each municipality can be mainly traced back to the presence of dominant attractors: the main urban areas The combination of the theory of central places (CPT), CA, and PCA, which is applied to the O / D matrix, represents the original contribution of the study, as the literature does not report similar examples Finally, the physical characteristics of the transport networks are considered through the isochrones; namely, the time needed to reach any other point in the transport network from a centroid. For this analysis, a network graph is used where each node (n) represents a point of intersection of the branches of the transport network, and the arcs (a) represents the physical connection between them. Each node

Sustainability 2019 , 11 , 4962 11 of 24 is characterised by its geographical coordinates and by the time required for its crossing. The arcs are defined not only by the nodes that also connect to the travelling speeds but also the flows that interest them. In this way, it is possible to determine the crossing time of each arc. On the graph so defined, the Dijkstra’s algorithm is applied [ 107 ]. Dijkstra’s algorithm is a procedure to find the shortest paths between nodes in a graph. The original algorithm has been applied to obtain the path between any two nodes and the total travel time, adding the time of each link. The process can be iterated to calculate the total time between one node (the source) and all other nodes of the graph. The algorithm explores the node starting from the source, considering the infinite time to reach other points. During the first iteration, the nodes connected to the source (border nodes) are considered by assigning the minimum value between infinity and the travel time of each link. At the first step, all nodes linked to the source have the correct total travel time. The source is marked as a visited node. In the second step, all border nodes are selected, and the previous procedure is repeated for each. At the end of the second step, the neighbouring border nodes have the correct total travel time. The border nodes are marked as visited nodes, and the neighbouring of border nodes became new border nodes Once the distribution of the time needed to reach all the other points of the network from a centroid is known, it is possible to evaluate the overall time and the average time of the centroid Finally, the trend line is investigated on the available population data; that is, its density and variation over time, the number of displacements, the average travel time from each centre, and altitude. This trend line explains the variance observed regarding maximum correlation ( R and R 2) with depopulation The correlation expresses any statistical relationship between two variables (causal or not) and can be understood as their dependence or association. Correlations can be useful tools in predicting relationships or revealing hidden associations. Among di ff erent approaches, Pearson’s correlation was applied, where R is evaluated as: R = σ xy σ x σ y = P n i = 1 x i − µ y y i − µ y r P n i = 1 x i − µ y 2 q P n i = 1 x i − µ y 2 If the variables are independent, then R = 0. For a perfect increasing linear dependence, R assumes the maximum value of + 1, while R = − 1 express the perfect opposite linear relationship dependence (anticorrelation) Moreover, a coe ffi cient of determination ( R 2) is evaluated to provide an estimation of how the data are well simulated by the model. There are several definitions of R 2. In the simple linear regression (when an intercept is included) R 2 is simply the square of the sample correlation coe ffi cient R 6. Application to the Case Study of Sardinia The percentage of internal displacements in Sardinia is equal to 98.4% of the total. The island’s daily mobility takes place predominantly as a closed system; movements are prevalently directed to and from areas of the same region. Considering the theory of central places that, starting from the concept of urban hierarchy, interpreted the size, frequency, and distance of urban centres, the territory of the island is analyzed to highlight its structure, subsystems, and related relationships.

Sustainability 2019 , 11 , 4962 12 of 24 Therefore, for the cluster analysis, a Euclidean metric is used, and a hierarchical aggregative method is applied to 35 clusters compared to 377 initial objects (the 377 Sardinian municipal administrations) The O / D matrix (377 × 377) in which each row represents the number of shifts facing the remaining 376 municipalities and the municipality itself. The diagonal represents intra-municipal movements (exiting from i and facing i ). It is possible to interpret these shifts as the intensity of the relationships or attraction of i towards other gravitational areas. The cluster analysis was conducted using the “Statistics and Machine Learning” Toolbox in the Matlab environment [ 108 ] with the k-means technique (Table 3 ). Given the assessment of accessibility, Figure 4 shows how the network graph was outlined, starting from the geometric and velocity characteristics of each stretch of the transport networks The application of Dijkstra’s algorithm on the graph allows for determining the travel times Sustainability 2019 , 11 , x FOR PEER REVIEW 12 of 24 ( a ) ( b ) ( c ) Figure 4. ( a ) The mean speed of road network; ( b ) the network graph; ( c ) the isochronous with the centre source of Cagliari. Table 4 shows the multiplicity of data on which a systematic statistical analysis is performed, in order to identify the maximum correlations between the physical and functional characteristics of the transport network (isochrones, altitude, average times), displacement indices (generation and attraction, inter-municipal and intra-municipal ones) and population changes. Figure 4. ( a ) The mean speed of road network; ( b ) the network graph; ( c ) the isochronous with the centre source of Cagliari Table 4 shows the multiplicity of data on which a systematic statistical analysis is performed, in order to identify the maximum correlations between the physical and functional characteristics of the transport network (isochrones, altitude, average times), displacement indices (generation and attraction, inter-municipal and intra-municipal ones) and population changes.

Sustainability 2019 , 11 , 4962 13 of 24 Table 3. The daily matrix origin–destination (O / D) with the 377 municipalities of Sardinia Municipalities Abbasanta Aggius Aglientu Aidomaggiore Al à dei Sardi Albagiara Ales Alghero Allai Anela Arborea Arbus Ardara V illanova T ruschedu V illanova T ulo V illanovaforru V illanovafranca V illaperuccio V illaputzu V illasalto V illasimius V illasor V illaspeciosa V illaurbana Zeddiani Zerfaliu Abbasanta 1190 0 0 6 0 0 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 2 0 2 Aggius 0 780 12 0 0 0 0 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Aglientu 0 8 504 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Aidomaggiore 12 0 0 110 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Al à dei Sardi 0 0 0 0 1048 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Albagiara 0 0 0 0 0 44 10 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 2 0 Ales 0 0 0 0 0 6 744 0 0 0 0 2 0 0 0 2 0 0 0 0 0 0 0 4 0 0 Alghero 0 2 0 0 2 0 0 29,024 0 0 6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Allai 2 0 0 0 0 0 0 0 166 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 0 0 Anela 0 0 0 0 0 0 0 2 0 246 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Arborea 4 0 0 0 0 0 2 0 0 0 2878 6 0 0 0 0 0 0 0 0 0 2 0 2 0 0 Arbus 0 0 0 0 0 0 0 0 0 0 6 2978 0 0 0 2 0 2 2 0 0 10 2 0 0 0 Ardara 0 0 0 0 0 0 0 10 0 0 0 0 358 0 0 0 0 0 0 0 0 0 0 0 0 0 Villanova Truschedu 6 0 0 0 0 0 0 0 4 0 0 0 0 64 0 0 0 0 0 0 0 0 0 0 0 0 Villanova Tulo 0 0 0 0 0 0 0 0 0 0 2 0 0 0 612 0 0 0 0 0 0 0 0 0 0 0 Villanovaforru 0 0 0 0 0 2 24 0 0 0 0 0 0 0 0 218 0 0 0 0 0 0 0 0 0 0 Villanovafranca 0 0 0 0 0 0 2 0 0 0 0 0 0 0 0 4 564 0 0 0 0 4 0 0 0 0 Villaperuccio 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 344 0 0 0 4 0 0 0 0 Villaputzu 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2496 4 62 0 0 0 0 0 Villasalto 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 18 582 2 0 0 0 0 0 Villasimius 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 14 0 2652 2 2 0 0 0 Villasor 0 0 0 0 0 0 0 0 0 0 0 4 0 0 0 0 2 0 4 0 2 2926 14 0 0 0 Villaspeciosa 0 0 0 0 0 0 0 0 0 0 2 0 0 0 0 0 0 0 0 0 0 30 870 0 0 0 Villaurbana 0 0 0 0 0 0 6 0 4 0 12 0 0 0 0 0 0 0 0 0 0 0 0 626 0 0 Zeddiani 8 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 354 0 Zerfaliu 4 0 0 0 0 0 0 0 0 0 14 0 0 2 0 0 0 0 0 0 0 0 0 4 8 386

Sustainability 2019 , 11 , 4962 14 of 24 Table 4. The set and organization of available data Ord. Isocr . Number Name East North Land Area Quote 1999 Inhabitants 2001 Population 2008 Population 2009 Population 2001 Intra 2011 Intra 2001 Gen Inter 2001 Att Inter 2011 Gen Inter 2011 Att Inter T ot 2001 Gen Inter T ot 2001 Att Inter T ot 2011 Gen Inter T ot 2011 Att Inter 2001 Population Density 2011 Population Density 2001 Gen. Interc. / Intrac. 2001 Attr . Interc. / Intrac. 2001 Intrac. / ab 2001 Gen. Interc. / ab 2001 Gen. Interc. / Att Interc 2001 Gen. T ot / Att tot (ab 2001–ab 1999) / ab 1999 2011 Gen. Interc. / Intrac. 2011 Attr . Interc. / Intrac. 2011 Intrac. / ab 2011 Gen. Interc. / ab 2011 Gen. Interc. / Att Interc 2011 Gen. T ot / Att tot (ab 2009–ab 2001) / ab 2001 (ab 2009–ab 2001) / ab 2001 for 9 Y ears Accessibility Index 1 Abbasanta 1,484,174 4,440,529 39.9 315 2640 2815 2873 2877 1282 1190 2338 2342 1230 1050 3620 3624 2420 2240 70.64 72.20 0.55 1.83 0.46 0.83 1.00 1.00 6.6% 1.034 0.882 0.414 0.428 1.171 1.08 2.2% 0.02 1.03 2 Aggius 1,504,417 4,538,322 83.6 514 1787 1686 1642 1629 812 780 1444 1026 660 213 2256 1838 1440 993 20.18 19.49 0.56 1.26 0.48 0.86 1.41 1.23 − 5.7% 0.846 0.273 0.479 0.405 3.099 1.45 − 3.4% − 0.03 0.73 3 Aglientu 1,512,856 4,550,908 148.6 420 1102 1093 1198 1207 526 504 720 646 200 197 1246 1172 704 701 7.36 8.12 0.73 1.23 0.48 0.66 1.11 1.06 − 0.8% 0.397 0.391 0.418 0.166 1.015 1.004 10.4% 0.10 1.11 4 Aidomaggiore 1,484,703 4,447,784 41.3 250 617 542 473 478 176 110 376 238 182 78 552 414 292 188 13.11 11.57 0.47 1.35 0.32 0.69 1.58 1.33 − 12.2% 1.655 0.709 0.23 0.381 2.333 1.553 − 11.8% − 0.12 2.85 5 Al à dei Sardi 1,524,459 4,501,178 188.6 663 2052 1949 1949 1963 842 1048 1136 906 364 88 1978 1748 1412 1136 10.33 10.41 0.74 1.08 0.43 0.58 1.25 1.13 − 5.0% 0.347 0.084 0.534 0.185 4.136 1.243 0.7% 0.01 1.04 7 Ales 1,480,676 4,403,508 21.7 194 1691 1628 1550 1533 800 744 1248 1974 376 1170 2048 2774 1120 1914 75.20 70.81 0.64 2.47 0.49 0.77 0.63 0.74 − 3.7% 0.505 1.573 0.485 0.245 0.321 0.585 − 5.8% − 0.06 1.88 10 Anela 1,508,106 4,483,800 37.0 446 939 817 716 712 318 246 592 378 246 96 910 696 492 342 22.10 19.26 0.54 1.19 0.39 0.72 1.57 1.31 − 13.0% 1 0.39 0.346 0.346 2.563 1.439 − 12.9% − 0.13 3.21 11 Arborea 1,465,970 4,403,756 115.5 7 3785 3927 3991 4028 2434 2878 3544 3576 1204 1686 5978 6010 4082 4564 34.00 34.87 0.69 1.47 0.62 0.90 0.99 0.99 3.8% 0.418 0.586 0.714 0.299 0.714 0.894 2.6% 0.03 0.61 12 Arbus 1,465,374 4,372,350 267.2 311 7596 7021 6677 6659 3198 2978 4950 3994 1672 891 8148 7192 4650 3869 26.28 24.93 0.65 1.25 0.46 0.71 1.24 1.13 − 7.6% 0.561 0.299 0.447 0.251 1.877 1.202 − 5.2% − 0.05 1.02 13 Ardara 1,485,613 4,494,413 38.1 296 839 847 820 822 452 358 722 530 386 74 1174 982 744 432 22.25 21.59 0.63 1.17 0.53 0.85 1.36 1.20 1.0% 1.078 0.207 0.436 0.47 5.216 1.722 − 3.0% − 0.03 1.34 14 Ardauli 1,499,141 4,435,869 20.6 421 1352 1158 1010 994 366 278 668 514 300 128 1034 880 578 406 56.35 48.37 0.55 1.40 0.32 0.58 1.30 1.18 − 14.3% 1.079 0.46 0.28 0.302 2.344 1.424 − 14.2% − 0.14 2.24 16 Armungia 1,543,460 4,380,741 54.8 366 668 584 521 510 226 178 340 278 130 50 566 504 308 228 10.66 9.31 0.66 1.23 0.39 0.58 1.22 1.12 − 12.6% 0.73 0.281 0.349 0.255 2.6 1.351 − 12.7% − 0.13 3.24 20 Assolo 1,495,160 4,405,686 16.3 255 544 485 455 447 88 88 294 162 148 18 382 250 236 106 29.72 27.39 0.30 1.84 0.18 0.61 1.81 1.53 − 10.8% 1.682 0.205 0.197 0.331 8.222 2.226 − 7.8% − 0.08 3.63 22 Atzara 1,504,417 4,423,283 35.8 540 1332 1310 1238 1230 546 456 922 644 296 78 1468 1190 752 534 36.58 34.35 0.59 1.18 0.42 0.70 1.43 1.23 − 1.7% 0.649 0.171 0.371 0.241 3.795 1.408 − 6.1% − 0.06 1.89 23 Austis 1,505,993 4,439,459 50.7 737 1054 959 897 891 422 392 670 466 200 88 1092 888 592 480 18.91 17.57 0.63 1.10 0.44 0.70 1.44 1.23 − 9.0% 0.51 0.224 0.44 0.224 2.273 1.233 − 7.1% − 0.07 2.70 24 Badesi 1,492,910 4,533,934 30.7 102 1860 1862 1892 1883 954 1022 1368 1272 480 342 2322 2226 1502 1364 60.63 61.32 0.70 1.33 0.51 0.73 1.08 1.04 0.1% 0.47 0.335 0.543 0.255 1.404 1.101 1.1% 0.01 0.91 25 Ballao 1,530,132 4,379,718 46.7 90 1086 971 930 907 374 328 640 444 224 86 1014 818 552 414 20.80 19.43 0.58 1.19 0.39 0.66 1.44 1.24 − 10.6% 0.683 0.262 0.362 0.247 2.605 1.333 − 6.6% − 0.07 2.47 26 Banari 1,472,237 4,495,780 21.3 419 756 677 620 617 224 186 438 256 244 64 662 480 430 250 31.83 29.01 0.51 1.14 0.33 0.65 1.71 1.38 − 10.4% 1.312 0.344 0.301 0.395 3.813 1.72 − 8.9% − 0.09 1.95

Sustainability 2019 , 11 , 4962 15 of 24 7. Results The cluster analysis and the PCA highlighted a composite system, divided into subsystems, within which the hubs that provide services to the insular territory emerge. The dependencies are not well-defined; they are blurred, especially in the border municipalities. Unlike the Christallerian theory that previously indicated that urban centres are arranged at regular distances for perfect competition, the homogeneity of the territorial areas in the groups formed through the cluster analysis is variable This observation is evident in Figure 5 where, in addition to the analysis in main components, the level of correlation that each municipality has with the overall tendency of the group to which it belongs is represented Moreover, the analysis correctly captures the perimeter of the provincial administrative limits except for some border municipalities which, as already mentioned, are not strongly polarised on a single subsystem. Figure 5 , in the three images, shows how each area has a more complex articulation than its administrative boundaries, almost coinciding with the historical areas of Sardinia. It is remarkable to observe how the mathematical analysis, conducted aseptically on the O / D matrix of Istat 2015 mobility, can capture the design of the socio-economic relationships generated by history and orography. Figure 5 shows the results of the cluster analysis, PCA, and the historical regions of Sardinia Sustainability 2019 , 11 , x FOR PEER REVIEW 15 of 24 7. Results The cluster analysis and the PCA highlighted a composite system, divided into subsystems, within which the hubs that provide services to the insular territory emerge. The dependencies are not well-defined; they are blurred, especially in the border municipalities. Unlike the Christallerian theory that previously indicated that urban centres are arranged at regular distances for perfect competition, the homogeneity of the territorial areas in the groups formed through the cluster analysis is variable. This observation is evident in Figure 5 where, in addition to the analysis in main components, the level of correlation that each municipality has with the overall tendency of the group to which it belongs is represented. Moreover, the analysis correctly captures the perimeter of the provincial administrative limits except for some border municipalities which, as already mentioned, are not strongly polarised on a single subsystem. Figure 5, in the three images, shows how each area has a more complex articulation than its administrative boundaries, almost coinciding with the historical areas of Sardinia. It is remarkable to observe how the mathematical analysis, conducted aseptically on the O/D matrix of Istat 2015 mobility, can capture the design of the socio-economic relationships generated by history and orography. Figure 5 shows the results of the cluster analysis, PCA, and the historical regions of Sardinia. ( a ) ( b ) ( c ) Figure 5. Results of ( a ) the cluster analysis; ( b ) PCA; ( c ) the historical regions of Sardinia. The analyses only partially frame the functional dependencies and the coherent behaviour of urban centres, as they refer to the cadenced mobility for study and work reasons. However, this has led to further analysis of the density of the infrastructural endowment to understand how to make inland areas more accessible to coastal areas, which are the areas with the most services. Starting from one of the initial postulates of Chirstaller’s central location theory that saw the increase in distance as a real problem for the demand for a service [110], the analysis of the territorial density of the infrastructural endowment shows a limited territorial coverage compared to the Italian average. In particular, the state network is characterised by an infrastructure density that is less than about 20% of the national one, while the provincial one, with a development of 22.63 km per 100 kmq compared to the national value of 37.81 km per 100 kmq, has a negative gap of 40%. Figure 5. Results of ( a ) the cluster analysis; ( b ) PCA; ( c ) the historical regions of Sardinia The analyses only partially frame the functional dependencies and the coherent behaviour of urban centres, as they refer to the cadenced mobility for study and work reasons. However, this has led to further analysis of the density of the infrastructural endowment to understand how to make inland areas more accessible to coastal areas, which are the areas with the most services. Starting from one of the initial postulates of Chirstaller’s central location theory that saw the increase in distance as a real problem for the demand for a service [ 109 ], the analysis of the territorial density of the infrastructural endowment shows a limited territorial coverage compared to the Italian average. In particular, the state network is characterised by an infrastructure density that is less than about 20% of the national one,

Sustainability 2019 , 11 , 4962 16 of 24 while the provincial one, with a development of 22.63 km per 100 kmq compared to the national value of 37.81 km per 100 kmq, has a negative gap of 40% The part of the territory classified as “mountain” is not considered; the data significantly worsens In flat and hilly areas, the consistency index of the provincial network is less than half of the national data, while the state network has an index with a negative gap of around 40%. The traditional indexes (network km per kmq of territory or network km per inhabitants) used by ISTAT do not allow for highlighting the real criticalities of the regional road network. What makes the situation serious is not so much the overall development but, above all, the standards that it possesses. In the Sardinia region, over 35% of the road network has sections of less than 7 m while the national data is 13.5%. Furthermore, road infrastructures with a section of fewer than 9 m on the island represent 73%, while the national average is 59%. Figure 6 a shows the territorial distribution of infrastructures according to their width (including the carriageway, the platforms and the marginal elements). It is significant to observe how the province of Nuoro, historically characterised by minimum accessibility values, is endowed with a considerable development of state and provincial roads; however, these are widely characterized by sections of less than 8 m. The speeds allowed by the current road network are also shown (Figure 6 b). It exemplifies, even more clearly, not only the areas of the regional territory in which the most disadvantageous accessibility conditions are verified (Figure 6 c) but also the lack of spatial balance, in terms of time and transport costs, which the theory of central places sets as a fundamental principle [ 110 ]. Sustainability 2019 , 11 , x FOR PEER REVIEW 16 of 24 The part of the territory classified as “mountain” is not considered; the data significantly worsens. In flat and hilly areas, the consistency index of the provincial network is less than half of the national data, while the state network has an index with a negative gap of around 40%. The traditional indexes (network km per kmq of territory or network km per inhabitants) used by ISTAT do not allow for highlighting the real criticalities of the regional road network. What makes the situation serious is not so much the overall development but, above all, the standards that it possesses. In the Sardinia region, over 35% of the road network has sections of less than 7 m while the national data is 13.5%. Furthermore, road infrastructures with a section of fewer than 9 m on the island represent 73%, while the national average is 59%. Figure 6 a shows the territorial distribution of infrastructures according to their width (including the carriageway, the platforms and the marginal elements). It is significant to observe how the province of Nuoro, historically characterised by minimum accessibility values, is endowed with a considerable development of state and provincial roads; however, these are widely characterized by sections of less than 8 m. The speeds allowed by the current road network are also shown (Figure 6 b). It exemplifies, even more clearly, not only the areas of the regional territory in which the most disadvantageous accessibility conditions are verified (Figure 6 c) but also the lack of spatial balance, in terms of time and transport costs, which the theory of central places sets as a fundamental principle [111]. ( a ) ( b ) ( c ) Figure 6. ( a ) Distribution of the road network in Sardinia; ( b ) speed distribution on the road network; ( c ) areas with inadequate speeds. Once these analyses were defined, the study analysed the accessibility of the newly defined homogeneous areas regarding the depopulation that particularly affects the internal areas of Sardinia. Figure 7 summarises the results of the analyses. The image shows the depopulation values in percentage terms on the left and the accessibility index on the right. This result was defined as the product of the average times required by each centre to reach the other 377, and the average altitude above sea level normalised according to the extreme values of depopulation. According to this analysis, 68.2% of the centres are in a phase of depopulation, and about 13.3% show strong abandonment values, higher than 15%. The most affected centres are in the inner and mountainous part of the island. Moreover, 94% of the municipalities included in a 10 km coastal strip have a growing population, and the top 10 coastal centres have growth rates of over 23%. Figure 6. ( a ) Distribution of the road network in Sardinia; ( b ) speed distribution on the road network; ( c ) areas with inadequate speeds Once these analyses were defined, the study analysed the accessibility of the newly defined homogeneous areas regarding the depopulation that particularly a ff ects the internal areas of Sardinia Figure 7 summarises the results of the analyses. The image shows the depopulation values in percentage terms on the left and the accessibility index on the right. This result was defined as the product of the average times required by each centre to reach the other 377, and the average altitude above sea level normalised according to the extreme values of depopulation.

Sustainability 2019 , 11 , 4962 17 of 24 According to this analysis, 68.2% of the centres are in a phase of depopulation, and about 13.3% show strong abandonment values, higher than 15%. The most a ff ected centres are in the inner and mountainous part of the island. Moreover, 94% of the municipalities included in a 10 km coastal strip have a growing population, and the top 10 coastal centres have growth rates of over 23% The comparison between the two images of Figure 7 clearly shows the correlation between the poor conditions of accessibility and high values of decrease. The link between the two indices is non-linear with a correlation coe ffi cient R = 66.5% and an explained variance of R 2 = 44.3%. In other words, the causes of depopulation can be significantly traced to the lack of accessibility Sustainability 2019 , 11 , x FOR PEER REVIEW 17 of 24 The comparison between the two images of Figure 7 clearly shows the correlation between the poor conditions of accessibility and high values of decrease. The link between the two indices is nonlinear with a correlation coefficient R = 66.5% and an explained variance of R 2 = 44.3%. In other words, the causes of depopulation can be significantly traced to the lack of accessibility. ( a ) ( b ) Figure 7. Results of the conducted analyzes: ( a ) depopulation values in percentage terms; ( b ) the accessibility index. 8. Discussion and Conclusions A Smart Governance Strategy is put forward in this study through the union between structural problems common to the islands and the theory of central places; this has allowed for observing how the study of mobility within a closed system (as an island associated with the vocation and needs of the population) is necessary to implement smart territorial policies. First, it was presented theoretically and, subsequently, it was demonstrated through the case study of Sardinia. The CPT (a means to understand the organisational dynamics based on urban hierarchies), which is linked to insularity, leads to reflections on which connection mechanisms are useful in a closed system where the coastal area is different (regarding the economy, morphology, and demography) from the internal area. Furthermore, the case study analysed leads to make some considerations for the major islands. Firstly, this study can make possible the definition of guidelines for the key points of regional island development (as shown in paragraph 3 under the National Strategy of the Italian Areas). In particular, the authors respond to the three points indicated by the National Strategy of the Italian Areas as the basis of development: (1) By understanding the attractive force that each region has within it, it is possible to reverse the negative economic and social trend that blocks the development of internal areas. This action can be done by modifying the primary services to change the polycentric articulation of the island territory; (2) By exploiting the unused internal capital and understanding how to modify the internal attractive force through basic services and the social costs generated by the production and consumption processes; Figure 7. Results of the conducted analyzes: ( a ) depopulation values in percentage terms; ( b ) the accessibility index 8. Discussion and Conclusions A Smart Governance Strategy is put forward in this study through the union between structural problems common to the islands and the theory of central places; this has allowed for observing how the study of mobility within a closed system (as an island associated with the vocation and needs of the population) is necessary to implement smart territorial policies. First, it was presented theoretically and, subsequently, it was demonstrated through the case study of Sardinia. The CPT (a means to understand the organisational dynamics based on urban hierarchies), which is linked to insularity, leads to reflections on which connection mechanisms are useful in a closed system where the coastal area is di ff erent (regarding the economy, morphology, and demography) from the internal area. Furthermore, the case study analysed leads to make some considerations for the major islands. Firstly, this study can make possible the definition of guidelines for the key points of regional island development (as shown in paragraph 3 under the National Strategy of the Italian Areas). In particular, the authors respond to the three points indicated by the National Strategy of the Italian Areas as the basis of development:

Sustainability 2019 , 11 , 4962 18 of 24 (1) By understanding the attractive force that each region has within it, it is possible to reverse the negative economic and social trend that blocks the development of internal areas. This action can be done by modifying the primary services to change the polycentric articulation of the island territory; (2) By exploiting the unused internal capital and understanding how to modify the internal attractive force through basic services and the social costs generated by the production and consumption processes; (3) By understanding that the depopulation phenomena remain one of the most important problems in Sardinia. Reversing the economic trend through the use of inland areas is possible to avoid the depopulation phenomena over time The challenge was to understand the strongest urban polarities and the interactions between the geographical areas of the system, which isolates through urban and regional planning data that allowed for developing a smart governance approach to the territory Thus, the design of the relationships, which the authors established through the analysis of the main components and the cluster analysis allowed for establishing the presence of stronger areas within the Sardinian insular system, highlighting a spontaneous creation (elaborated from population movements) to groupings of municipalities (of homogeneous areas) This, at first, allowed for observing how (1) the depopulation phenomena involve the mountain areas to a greater extent than the coastal and flat areas, (2) the outgoing displacements from each municipality can be mainly attributed to the presence of dominant attractors (main urban areas), and (3) the existence of insular subsystems comprises reference hubs and secondary zones Furthermore, the analysis, repeatable in other national and international insular contexts, represented the functional dependencies between municipalities; therefore, it is the strength with which the municipalities attract each other. In particular, through the overlap of the Historical Regions of Sardinia and the analysis of the main components, how the homogeneity of the groupings has proved variable and pertinent both to the historical regions and the administrative limits of today became evident The result is an organic system composed of the set of identified subsystems that is capable of achieving adequate accessibility conditions throughout the entire regional territory, thereby, making maximum integration between the infrastructural systems, the main vocations of the Sardinian territory (identity, demographic structure, history, culture), and population needs (work, health, services, leisure) Within each area, the hubs that need to be equipped with services are identified, and, at the same time, how to enhance accessibility within the homogeneous areas and, subsequently, how to improve accessibility between the di ff erent areas are studied Finally, for further research, this study allowed for having a starting point for developing a set of indicators to identify the polarities that have more services and fewer services, and to define a methodology that can solve similar problems in other insular contexts, considering the dynamics that emerged Author Contributions: This paper is the result of the joint work of the authors. In particular, “The analysis of urban polarities to understand the interconnection between inner end coastal areas in insular regions”, “The case study of the island of Sardinia in Italy”, “Methods”, “Application to the case study of Sardinia” and “Results” were written jointly by the authors. C.G. wrote “Introduction”, G.D. wrote “Literature review on common issues related to insular regions”. M.C. wrote “Discussion and Conclusions”. Conceptualization, C.G. and G.D.; data curation, C.G., G.D., and M.C.; formal analysis, C.G.; investigation, M.C.; methodology, C.G., G.D. and M.C.; supervision, C.G. and M.C.; validation, M.C.; writing—original draft, C.G., M.C. and M.C.; writing—review and editing, C.G Funding: This study was funded by the MIUR (Ministry of Education, Universities and Research [Italy]) through two projects: (1) Governing the smart city: a governance-centred approach to smart urbanism—GHOST (Project code: RBSI 14 FDPF; CUP Code: F 22 I 15000070008), financed with the SIR (Scientific Independence of Young Researchers) programme. (2) “CAGLIARI 2020” (Project code: PON 04 a 2_00381). We authorize the MIUR to reproduce and distribute reprints for governmental purposes, notwithstanding any copyright notations thereon Any opinions, findings, and conclusions or recommendations expressed in this material are those of the authors,

Sustainability 2019 , 11 , 4962 19 of 24 and do not necessarily reflect the views of the MIUR. This study was also supported by the project “Healthy Cities and Smart Territories,” founded by the Foundation of Sardinia and Autonomous Region of Sardinia (Fondazione di Sardegna—Convenzione triennale tra la Fondazione di Sardegna e gli Atenei Sardi Regione Sardegna 2016) Acknowledgments: The authors would like to thank the anonymous reviewers of Sustainability for their valuable contributions and suggestions, which have made possible significant improvements to the text Conflicts of Interest: The authors declare no conflict of interest References and Notes 1 Davoudi, S. Understanding territorial cohesion J. Plan. Pract. Res 2006 , 20 , 433–441. [ CrossRef ] 2 Schön, K.P. Territorial Cohesion—Current Views of the Commission and the Member States of the European Union Ger. Annu. Spat. Res. Policy 2009 , 7–17. [ CrossRef ] 3 Zaucha, J.; Böhme, K. Measuring territorial cohesion is not a mission impossible J. Eur. Plan. Stud 2019 [ CrossRef ] 4 Medeiros, E. Debating the Urban Dimension of Territorial Cohesion Territ. Cohes 2019 , 3–22. [ CrossRef ] 5 Consolidated Version of the Treaty on the Functioning of the European Union, O ffi cial EN Journal of the European Union C 326 / 47. 2012. Available online: https: // eur-lex.europa.eu / legal-content / EN / TXT / PDF / ?uri = CELEX:12012 E / TXT&from = IT (accessed on 17 June 2019) 6 European Structural and Investment Funds Regulations 2014–2020. 2015. Available online: https: // ec.europa. eu / regional_policy / en / information / legislation / regulations / (accessed on 17 June 2019) 7 Methods and Objectives for E ff ective Use of Community Funds 2014–2020. 2012. Available online: https: // www.sicilia-fse.it / Uploads / Temporaneo / Sintesi%20 Doc%20 Barca%20 Metodi-e-obiettivi-per-unuso-e ffi cace-dei-fondi-comunitari-2014-20.pdf (accessed on 17 June 2019) 8 Partnership Agreement 2014–2020 National Strategy for Internal Areas: Definition, Objectives, Tools and Governance Department for the Development of Social Coesiene 2013 Available online: file: /// C: / Users / user / Downloads / Accordo-P-Strategia_nazionale_per_le_Aree_interne_definizione_ obiettivi_strumenti_e_governance_2014%20 (accessed on 17 June 2019) 9 Coronato, M. The European Instrument to internal areas Boll. Della Assoc. Ital. Cartogr 2016 , 157 , 53–59 [ CrossRef ] 10 National Strategy for Internal Areas, Department for Planning and Coordination of Economic Policy. Available online: http: // www.programmazioneeconomica.gov.it / 2019 / 05 / 23 / strategia-nazionale-delle-aree-interne / (accessed on 17 June 2019) 11 Scorza, F.; Murgante, B.; Las Casas, G.; Fortino, Y.; Pilogallo, A. Investigating Territorial Specialization in Tourism Sector by Ecosystem Services Approach Mediterr. Cities Isl. Communities 2019 , 161–179. [ CrossRef ] 12 The Impact of Cohesion Policies in Europe and Italy. 2018. Available online: http: // www.senato.it / service / PDF / PDFServer / BGT / 01083823.pdf (accessed on 15 July 2019) 13 Access to European Union Law. O ffi cial Journal of the European Communities C 268 Volume 43. 2000 Available online: https: // eur-lex.europa.eu / legal-content / IT / TXT / ?uri = OJ%3 AC%3 A 2000%3 A 268%3 ATOC (accessed on 17 June 2019) 14 Amsterdam Treaty Available online: http: // www.eurotreaties.com / amsterdam-treaty / (accessed on 17 June 2019) 15 Pan, Y.; Zhai, M.; Lin, L.; Lin, Y. Characterizing the spatiotemporal evolutions and impact of rapid urbanization on island sustainable development Habitat Int 2016 , 53 , 215–227. [ CrossRef ] 16 Garau, C.; Desogus, G. A preliminary survey for a smart framework for the island contexts. In Proceedings of the IFAU 2018—2 nd International Forum on Architecture and Urbanism, Pescara, Italy, 8–10 November 2018 (accepted; in press) 17 Priano, F.H.; Armas, R.L.; Guerra, C.F. Developing smart regions: Proposal and application of a model for island territories Int. J. E Plan. Res 2018 , 7 , 89–114. [ CrossRef ] 18 Smart Islands Declaration, New Pathways for European Islands to Create Smart, Inclusive and Thriving Island Societies for an Innovative and Sustainable Europe. Available online: http: // www.smartislandsinitiative.eu / pdf / Smart_Islands_Declaration.pdf (accessed on 17 June 2019) 19 Smart Islands Projects and Strategies Issued from the 1 st European Smart Islands Forum, June 2016, Athens, Greece. Available online: http: // library.fes.de / pdf-files / bueros / athen / 12860.pdf (accessed on 17 June 2019).

Sustainability 2019 , 11 , 4962 20 of 24 20 Desogus, G.; Garau, C.; Mistretta, P. Smart Islands: A Systematic Review on Urban Policies and Smart Governance. In Proceedings of the ICCSA 2019—The International conference on computational science and its application, Saint Petersburg, Russia, 1–4 July 2019. (in press) 21 Las Casas, G.; Murgante, B.; Scorza, F. Regional Local Development Strategies Benefiting from Open Data and Open Tools and an Outlook on the Renewable Energy Sources Contribution Smart Energy Smart City 2016 , 275–290. [ CrossRef ] 22 Stratigea, A.; Papadopoulou, C.A.; Panagiotopoulou, M. Tools and Technologies for Planning the Development of Smart Cities J. Urban Technol 2015 , 22 , 43–62. [ CrossRef ] 23 State of Implementation of the Public Preliminary Investigation State-Regions / Autonomous Provinces Available online: http: // old 2018.agenziacoesione.gov.it / it / arint / Stato_di_attuazione / stato_di_attuazione.html (accessed on 17 June 2019) 24 Sardinia Programming The Strategy for Internal Areas Available online: https: // www. sardegnaprogrammazione.it / index.php?xsl = 1384&s = 278011&v = 2&c = 12950 (accessed on 17 June 2019) 25 Autonomous Region of Sardinia, Resolution n. 6 / 13 of 10.2.2015. Available online: http: // www.regione. sardegna.it / documenti / 1_274_20150211164206.pdf (accessed on 17 June 2019) 26 Treaty on European Union. O ffi cial Journal of the European Communities No C 191 / 3. Available online: https: // eur-lex.europa.eu / legal-content / IT / TXT / ?uri = CELEX:11992 M / TXT (accessed on 13 July 2019) 27 Analysis of the Island Regions and Outermost Regions of the European Union. Available online: https: // ec. europa.eu / regional_policy / sources / docgener / studies / pdf / ilesrup / islands_part 1_summary_en.pdf (accessed on 13 July 2019) 28 Monfort, P. Territories with Specific Geographical Features. Available online: https: // ec.europa.eu / regional_ policy / sources / docgener / work / 2009_02_geographical.pdf (accessed on 13 July 2019) 29 The ESPON 2013 Programme the Development of the Islands—European Islands and Cohesion Policy (EUROISLANDS). Available online: https: // www.espon.eu / sites / default / files / attachments / inception_report_ full_version.pdf (accessed on 13 July 2019) 30 The Statistical Data of EUROSTAT Is Based on the Definition of Island of the European Community Which Catalogs These as NUTS 3: “Islands Are Defined as Territories Having: A Minimum Surface of 1 km 2 ; a Minimum Distance between the Island and the Mainland of 1 km; a Resident Population of More Than 50 Inhabitants; and No Fixed Link between the Island and the Mainland” Based on the above Definitions, in Europe, There Are Just 76 Insular Regions. EUROSTAT Statistic Explained: Territorial Typologies Manual—Island Regions. Available online: https: // ec.europa.eu / eurostat / statistics-explained / index.php?title = Territorial_typologies_manual_ -_island_regions#Classes_for_the_typology_and_their_conditions (accessed on 16 July 2019). European Parliament, Parliamentary Questions. Available online: http: // www.europarl.europa.eu / doceo / document / E- 8-2018-000443-ASW_EN.html (accessed on 16 July 2019) 31 According to the Final Report on Analysis of the Island Regions of 2000, Europe’s Islands are 286 are Home to almost 10 million People Who Occupy an Area of 100,000 Km 2 . Analysis of the Island Regions and Outermost Regions of the European Union. Available online: https: // ec.europa.eu / regional_policy / sources / docgener / studies / pdf / ilesrup / islands_part 1_summary_en.pdf (accessed on 13 July 2019) 32 Annual General Meeting of the CPMR Islands Commission, Corfu (Ionian Islands, Greece). 21 March 2019 Available online: https: // cpmr-islands.org / download / final-declaration-adopted-by-the-islands-commissionin-corfu-on-march-2019 / ?wpdmdl = 4650&ind = 1554795693625 (accessed on 16 July 2019) 33 Final declaration approved by the Islands Commission March 21, 2019, 39 th Annual General Assembly of the CPMR Islands Commission. 2019. Available online: https: // cpmr-islands.org / download / final-declarationadopted-by-the-islands-commission-in-corfu-on-march-2019 / ?wpdmdl = 4650&ind = 1554795700924 (accessed on 15 July 2019) 34 Opinion of the Economic and Social Committee on ‘Guidelines for integrated actions on the island regions of the European Union following the Amsterdam Treaty (Article 158), Gazzetta u ffi ciale dell’Unione europea C 268, il 19 / 09 / 2000. Available online: https: // eur-lex.europa.eu / legal-content / EN / TXT / ?uri = CELEX: 52000 IE 0805 (accessed on 16 July 2019) 35 Moncada, S.; Camilleri, M.; Formosa, S.; Galea, R. Islands at the Periphery: Integrating the Challenges of Island Sustainability into European Policy. European FP 6 Project SENSOR. 2009. Available online: https: // www.um.edu.mt / europeanstudies / books / CD_CSP 5 / pdf / smoncada.pdf (accessed on 28 August 2019).

Sustainability 2019 , 11 , 4962 21 of 24 36 Pantazis, D.N.; Moussas, V.C.; Murgante, B.; Daverona, A.C.; Stratakis, P.; Nikolaos, V.; Kavadias, A.; Economou, D.; Santimpantakis, K.; Karathanasis, B.; et al. Smart Sustainable Islands vs Smart Sustainable Cities. In Proceedings of the ISPRS Annals of the Photogrammetry, Remote Sensing and Spatial Information Sciences, 2 nd International Conference on Smart Data and Smart Cities, Puebla, Mexico, 4–6 October 2017 37 Dominguez, C.D.; Hern á ndez, M.R.; Talavera, A.S.; L ó pez, E.P. Smart Island Tourism and Strategic Marketing: The Case of the Island of El Hierro. Available online: https: // pdfs.semanticscholar.org / c 411 / 6406 d 3 e 277780 f 018 e 7 d 4 adee 8 f 27 f 5 c 2106.pdf (accessed on 28 August 2019) 38 Banos-Gonz á lez, I.; Mart í nez-Fern á ndez, J.; Esteve, M.A. Tools for sustainability assessment in island socio-ecological systems: An application to the Canary Islands Isl. Stud. J 2016 , 11 , 9–34. Available online: https: // islandstudies.ca / sites / default / files / ISJ-11-1-B-BanosGonz%C 3%A 1 lez-et-al.pdf (accessed on 28 August 2019) 39 Abis, E.; Garau, C. An assessment of the e ff ectiveness of strategic spatial planning: A study of Sardinian municipalities Eur. Plan. Stud 2016 , 24 , 139–162. [ CrossRef ] 40 Tuugalei, I.; Mow, C. Issues and challenges, strategies and recommendations, in the development of ICT in a Small Island Developing State: The case of Samoa EJISDC 2014 , 63 , 1–24. [ CrossRef ] 41 Kurniawan, F.; Adrianto, L.; Bengen, D.G.; Prasetyo, L.B. Patterns of landscape change on small islands: A case of Gili Matra Islands, Marine Tourism Park, Indonesia Procedia Soc. Behav. Sci 2016 , 227 , 553–559 [ CrossRef ] 42 Lee, Y.B. Exploring the Relationship between E-Government Development and Environmental Sustainability: A Study of Small Island Developing States Sustainability 2017 , 9 , 732. [ CrossRef ] 43 Castanho, R.A. Dynamics of the Land Use Changes and the Associated Barriers and Opportunities for Sustainable Development on Peripheral and Insular Territories: The Madeira Island, Portugal. 2018. Available online: https: // www.intechopen.com / online-first / dynamics-of-the-land-use-changes-and-the-associatedbarriers-and-opportunities-for-sustainable-devel (accessed on 28 August 2019) 44 S á nchez, A.F.; Gil, F.S.; Sabater, L.A.; Dentinho, T.P. AQ Methodology approach to define urban sustainability challenges in a small insular city. In Proceedings of the 51 st European Congress of the Regional Association International, Barcelona, Spain, 30 August–3 September 2011 45 Andriotis, K. Researching the development gap between the hinterland and the coast—Evidence from the island of Crete Tour. Manag 2006 , 27 , 629–639. [ CrossRef ] 46 Papatheochari, T.; Coccossis, H. Development of a waterfront regeneration tool to support local decision making in the context of integrated coastal zone management Ocean Coast. Manag 2019 , 169 , 284–295 [ CrossRef ] 47 Vallega, A. Urban waterfront facing integrated coastal management Ocean Coast. Manag 2001 , 44 , 379–410 [ CrossRef ] 48 Novembre, C. Le aree interne della sicilia tra problemi di sviluppo e ricomposizione territorial Geogr. Ital 2015 , 122 , 235–253. Available online: https: // dialnet.unirioja.es / servlet / articulo?codigo = 5173381 (accessed on 26 July 2019) 49 Core, B. Spatial dynamics of Mediterranean coastal regions J. Coast. Conserv 1999 , 5 , 105–112. Available online: https: // link.springer.com / content / pdf / 10.1007%2 FBF 02802747.pdf (accessed on 26 July 2019) 50 Internal Areas: Which Territories Are We Talking About? Explanatory Note on the Method of Classifying Areas. Available online: http: // old 2018.agenziacoesione.gov.it / opencms / export / sites / dps / it / documentazione / Aree_interne / Nota_metodologica_Aree_interne.pdf (accessed on 15 July 2019) 51 Mulligan, G.F.; Partridge, M.D.; Carruthers, J.I. Central place theory and its reemergence in regional science Ann. Reg. Sci 2012 , 48 , 405–431. [ CrossRef ] 52 Doran, D.; Fox, A. Operationalizing Central Place and Central Flow Theory with Mobile Phone Data Ann. Data Sci 2015 , 3 , 1–24. [ CrossRef ] 53 Berry, B.J.I.; Baskin, C.W.; Christaller, W. Central Places in Southern Germany Econ. Geogr 1967 , 43 , 275–276 [ CrossRef ] 54 Woglom, H.W.; Stolper, W.F. The Economics of Location by August Lösch Am. Sociol. Rev 1995 , 20 , 370–361 [ CrossRef ] 55 Getis, A.; Getis, J. Christaller’s Central Place Theory J. Geogr 1966 . [ CrossRef ] 56 Fujita, M.; Ogawa, H.; Thisse, J.F. A spatial competition approach to central place theory: Some basic principles J. Reg. Sci 1988 , 28 , 477–494. [ CrossRef ]

Sustainability 2019 , 11 , 4962 22 of 24 57 Ye, X.; Liu, X Cities as Spatial and Social Networks , 1 st ed.; Springer: Cham, Switzerland, 2019; p. 41, ISBN 978-3-319-95350-2 58 Clark, J.; Harrison, J.; Miguelez, E. Connecting cities, revitalizing regions: The centrality of cities to regional development Reg. Stud 2018 , 52 , 1025–1028. [ CrossRef ] 59 Knitter, D.; Nakoinz, O. The Relative Concentration of Interaction-A Proposal for an Integrated Understanding of Centrality and Central Places Land 2018 , 22–39 , 86. [ CrossRef ] 60 Hsu, W.T.; Zou, X. Central Place Theory and the Power Law for Cities Math Urban Morphol 2019 , 55–75 [ CrossRef ] 61 Derudder, B.; Taylor, P.J. Central flow theory: Comparative connectivities in the world-city network Reg. Stud 2018 , 52 , 1029–1040. [ CrossRef ] 62 Mal ý , J. Polycentric Urban Systems and Territorial Cohesion. In Territorial Choesion the Urban Dimension , 1 st ed.; Springer: Cham, Switzerland, 2019; p. 41, ISBN 978-3-030-03386-6 63 Camagni, R. From city hierarchy to city network: Reflections about an emerging paradigm. In Seminal Studies in Regional and Urban Economics: Contributions from an Impressive Mind ; Springer: Cham, Switzerland, 2017 [ CrossRef ] 64 Krugman, P. Increasing Returns and Economic Geography J. Polit. Econ 1991 , 99. [ CrossRef ] 65 Von Böventer, E. Toward a united theory of spatial economic structure Pap. Reg. Sci 1962 , 10 , 163–187 [ CrossRef ] 66 Pred, A. The Interurban Transmission of Growth in Advanced Economies: Empirical Findings Versus Regional-Planning Assumptions Reg. Stud 1976 , 10 , 151–171. [ CrossRef ] 67 Davoudi, S. Polycentricity in European spatial planning: From an analytical tool to a normative agenda Eur. Plan. Stud 2003 , 11 , 979–999. [ CrossRef ] 68 Van Meeteren, M.; Poorthuis, A. Christaller and “big data”: Recalibrating central place theory via the geoweb Urban Geogr 2018 , 39 , 122–148. [ CrossRef ] 69 Chapman, G.P. The Application of Information Theory to the Analysis of Population Distributions in Space Econ. Geogr 2006 , 46 , 317–331. [ CrossRef ] 70 Mulligan, G.F. Agglomeration and Central Place Theory: A Review of the Literature Int. Reg. Sci. Rev 1984 , 9 , 1–42. [ CrossRef ] 71 Parr, J.B. Competition, supply areas and industrial location: An equilibrium perspective Ann. Reg. Sci 1993 , 27 , 191–210. [ CrossRef ] 72 Loumeau, G Locating Public Facilities: Theory and Micro Evidence from Paris ; KOF Working Papers; KOF Swiss Economic Institute: Zurich, Switzerland, 2019; Volume 452. [ CrossRef ] 73 Godwin Onokerhoraye, A.G. A conceptual framework for the location of public facilities in the urban areas of developing countries: The Nigerian case Socio Econ. Plan. Sci 1976 , 10 , 237–240. [ CrossRef ] 74 Teitz, M.B. Toward a theory of urban public facility location Pap. Reg. Sci. Assoc 1968 , 21 , 35–51. [ CrossRef ] 75 Beckmann, M.J.; Thisse, J.F. The location of production activities Handb. Reg. Urban Econ 1987 , 1 , 21–95 [ CrossRef ] 76 Jovanovic, M.N. Spatial Location of Firms and Industries: An Overview of Theory Econ. Internazionale 2003 , 56 , 23–82. [ CrossRef ] 77 Xu, F.; Zhen, F.; Qin, X.; Wang, X.; Wang, F. From central place to central flow theory: An exploration of urban catering Tour. Geogr 2019 , 21 , 121–142. [ CrossRef ] 78 Butowski, L. Tourist sustainability of destination as a measure of its development Curr. Issues Tour 2017 , 22 , 1043–1061. [ CrossRef ] 79 Daniels, M.J. Central place theory and sport tourism impacts Ann. Tour. Res 2007 , 34 , 332–347. [ CrossRef ] 80 Brown, S. Christaller Knew My Father: Recycling Central Place Theory J. Macromarketing 1995 , 15 , 60–72 [ CrossRef ] 81 Dennis, C.; Marsland, D.; Cockett, T. Central place practice: Shopping centre attractiveness measures, hinterland boundaries and the UK retail hierarchy J. Retail. Consum. Serv 2002 , 9 , 185–199. [ CrossRef ] 82 Stratigea, A.; Katsoni, V. A strategic policy scenario analysis framework for the sustainable tourist development of peripheral small island areas—The case of Lefkada-Greece Island Eur. J. Futures Res 2015 , 3. [ CrossRef ] 83 Choy, D.J.L. Life Cycle Models for Pacific Island Destinations J. Travel Res 1992 , 30 , 26–31. [ CrossRef ] 84 Fennell, A.F. A tourist space-time budget in the Shetland Islands Ann. Tour. Res 1996 , 24 , 811–829 [ CrossRef ]

Sustainability 2019 , 11 , 4962 23 of 24 85 Rodr í guez, J.R.O.; Parra-L ó pez, E.; Yanes-Est é vez, V. The sustainability of island destinations: Tourism area life cycle and teleological perspectives. The case of Tenerife Tour. Manag 2008 , 29 , 53–65. [ CrossRef ] 86 Nguyen, D.N.; Imamura, F.; Iuchi, K. Barriers towards hotel disaster preparedness: Case studies of post 2011 Tsunami, Japan Int. J. Disaster Risk Reduct 2018 , 28 , 585–594. [ CrossRef ] 87 Volgmann, K.; Rusche, K. The Geography of Borrowing Size: Exploring Spatial Distributions for German Urban Regions Early View 2019 . [ CrossRef ] 88 Van Meeteren, M. Urban System. In The Wiley Blackwell Encyclopedia of Urban and Regional Studies ; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2019. [ CrossRef ] 89 Mul í ˇcek, O.; Mal ý , J. Moving towards more cohesive and polycentric spatial patterns? Evidence from the Czech Republic Reg. Sciance 2019 , 98 , 1177–1194. [ CrossRef ] 90 Mitchell, R.E. Individual Settlements Are Members of Larger Settlement Systems. In Human Geographies within the Pale of Settlement , 1 st ed.; Mitchell, R.E., Ed.; Palgrave Macmillan: Cham, Switzerland, 2018; pp. 71–91, ISBN 978-3-319-99144-3 91 Lenormand, M.; Ramasco, J.J. Towards a Better Understanding of Cities Using Mobility Data Built Environ 2016 , 42 , 356–364. [ CrossRef ] 92 Hansen, N. Border regions: A critique of spatial theory and a European case study Ann. Reg. Sci 1977 , 11 , 1–14. [ CrossRef ] 93 Azzari, M.; Garau, C.; Nesi, P.; Paolucci, M.; Zamperlin, P. Smart City Governance Strategies to better move towards a Smart Urbanism Comput. Sci. Its Appl. ICCSA 2018 , 10962 , 639–653. [ CrossRef ] 94 Demographic Statistics, Population Sardinia 2001–2018. Available online: https: // www.tuttitalia.it / sardegna / statistiche / popolazione-andamento-demografico / (accessed on 23 June 2019) 95 Italian Regions Classified by Area. Available online: https: // www.tuttitalia.it / regioni / superficie / (accessed on 23 June 2019) 96 Regional law 4 February 2016, n. 2. Reorganization of the system of local autonomies of Sardinia. O ffi cial bulletin of the Sardinia region n. 6 of 11 February 2016 97 ISTAT Statistical Codes of Territorial Administrative Units: Municipalities, Metropolitan Cities, Provinces and Regions Available online: https: // www.istat.it / it / archivio / 6789# Elencodeicodiciedelledenominazionidelleunitterritoriali-0 (accessed on 23 June 2019) 98 CRENOS. Economy of Sardinia, 25 th Report 2018. Available online: http: // crenos.unica.it / crenosterritorio / sites / default / files / allegati-pubblicazioni-tes / Economia%20 della%20 Sardegna%20-%2025%C 2%B 0% 20 Rapporto%20 CRENoS%202018.pdf (accessed on 15 July 2019) 99 Fisale, K.; Shaker, A. An Investigation of GIS Overlay and PCA Techniques for Urban Environmental Quality Assessment: A Case Study in Toronto, Ontario, Canada Sustainability 2017 , 9 , 380. [ CrossRef ] 100. Li, X.; Gar-On Yeh, A. Urban Simulation Using Principal Components Analysis and Cellular Automata for Land-Use Planning J. Photogramm. Eng. Remote Sens 2002 , 68 , 341–351 101. Li, X.; Gar-On Yeh, A. Principal component analysis of stacked multi-temporal images for the monitoring of rapid urban expansion in the Pearl River Delta Int. J. Remote Sens. Taylor Fr 1998 , 19 , 1501–1518. [ CrossRef ] 102. Pacione, M. Urban environmental quality and human wellbeing-a social geographical perspective J. Landsc Urban Plan 2003 , 65 , 19–30. [ CrossRef ] 103. Bonaiuto, M.; Fornara, F.; Bonnes, M. Indexes of perceived residential environment quality and neighbourhood attachment in urban environments: A confirmation study on the city of Rome J. Landsc. Urban Plan 2003 , 65 , 41–52. [ CrossRef ] 104. Everitt, B.S.; Landau, S.; Leese, M.; Stahl, D Cluster Analysis , 5 th ed.; John Wiley & Sons: Hoboken, NJ, USA, 2011; ISBN 9780470749913 105. Jolli ff e, I. Principal Component Analysis In International Encyclopedia of Statistical Science , 1 st ed.; Lovric, M., Ed.; Springer: Berlin / Heidelberg, Germany, 2011. Available online: https: // link.springer. com / referenceworkentry / 10.1007%2 F 978-3-642-04898-2_455 (accessed on 15 July 2019) 106. Krzanowski, W Principles of Multivariate Analysis , 2 nd ed.; Oxford University Press: Oxford, UK, 2000; ISBN 978-0-19-850708-6 107. Cormen, T.H.; Leiserson, C.E.; Rivest, R.L.; Stein, C Introduction to Algorithms , 2 nd ed.; The MIT Press: Cambridge, MA, USA, 2001; ISBN 0-262-03293-7.

Sustainability 2019 , 11 , 4962 24 of 24 108. Statistics and Machine Learning Toolbox ™ User’s Guide. 2018. Available online: https: // it.mathworks.com / help / pdf_doc / stats / stats.pdf (accessed on 26 July 2019) 109. Parr, J.B. Central place theory: An evaluation Rev. Urban Reg. Dev. Stud 2017 , 11 , 256–14. [ CrossRef ] 110. Spatial Balance and Central Locations, Lesson 4. Available online: http: // www.gpt.adhoc.ch / Geoeconomica / Modulo 2 / Lezione 4.pdf (accessed on 23 June 2019) © 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http: // creativecommons.org / licenses / by / 4.0 / ).