Why Italy First? Health, Geographical and Planning Aspects of the COVID-19 Outbreak

: COVID-19 hit Italy in February 2020 after its outbreak in China at the beginning of January. Italy the of What characteristics and di ﬀ usion patterns and hypothesized its outbreak to the of after containment including a national lockdown, we to to these questions, analyzing the issue from medical, geographical and planning of With reference to the Italian case, we observed the phenomenon in terms of the spatial di ﬀ usion process and by observing the relation between the epidemic and various environmental elements. In particular, we started from a hypothesis of the comparable economic, geographical, climatic and environmental conditions of the areas of Wuhan (in the Hubei Province in China, where the epidemic broke out) and the Po Valley area (in Italy) where most cases and deaths were registered. Via an ecological approach, we compared the spatial distribution and pattern of COVID-19-related mortality in Italy with several geographical, environmental and socio-economic variables at a Provincial level, analyzing them by means of spatial analytical techniques such as LISA (Local Indicators of Spatial Association). Possible evidence arose relating to COVID-19 cases and Nitrogen-related pollutants and land take, particularly in the Po Valley area.


Why Italy (1)? The Epidemiologic Point of View
In December 2019, in the Wuhan province of China, a new Coronavirus emerged following a spillover; this RNA virus was roughly 80% homogenous with the SARS virus, hence the name SARS-Cov-2 (Severe Acute Respiratory Syndrome Coronavirus 2) [1]. This led to a widespread epidemic of a new respiratory disease (COVID- 19), which, in three months, had crossed the borders of Asia becoming a pandemic with over 2,300,000 cases and 160,000 deaths as of 21st April 2020 [2].
The disease is spread by way of inter-human transmission, through the Flugge droplets, although it is also airborne if aerosols are generated and, as with normal influenza, it can also be transmitted Some similarities can be found in terms of the dimensions of the Province of Hubei with 58.52 million inhabitants (2015) and 158,000 km 2 , and those of Italy with 60,359,546 inhabitants and 302,072.84 km 2 . Although not an easy comparison, it is however possible to trace geo-climatic similarities, as well as those concerning human activities. In particular, both areas correspond to the Cfa subclass in the Köppen climate classification system [20] as 'humid subtropical', typical of similarities, as well as those concerning human activities. In particular, both areas correspond to the Cfa subclass in the Köppen climate classification system [20] as 'humid subtropical', typical of temperate continental areas. Both are located in an alluvial plain, the Wuhan urban agglomeration-Yangtze river, and the Greater Milan metropolitan area-Po river.
The Wuhan urban agglomeration is the most important hub in central China at the crossroads of corridors that connect Northern and Southern China [21] as well as internal China and the coast. Also, Chinese national highways cross Wuhan, as well as the Shanghai-Chengdu and Beijing-Hong Kong-Macau expressways [22]. In addition, it lays also at the centre of the Beijing-Wuhan-Guangzhou line, China's most important high-speed railway. Lastly, the Wuhan urban agglomeration hosts the international airport of Wuhan Tianhe, moving around 25 million passengers in 2018, the major hub of central China, with direct connections with mainland China, Western Europe and the USA.
The Greater Milan metropolitan area represents the most important urban and industrial agglomeration in Italy and is a connection to Central and Northern Europe. The main Italian national highways cross Milan, namely the West-East Turin-Trieste and the North-South Milan-Naples Expressways. Moreover, national and international High Speed Railways converge in Milan, connecting the area with major European Cities and linking it to the national major metropolitan areas. Lastly, the Greater Milan metropolitan area hosts the international airports of Linate, Malpensa and Bergamo, moving 49.3 million passengers in 2019, the second major hub system in Italy, with direct connections to Europe, China and the USA [23].
Both mega urbanizations have industrial and post-industrial functions, with a heavy presence of manufacturing companies, in machinery, automotive and ICT, as well as advanced and cultural services, particularly in the major centers. Both areas share a strong intermingling with agricultural activities and a wide progression of the sprawl [24][25][26][27].

Air Quality
In attempting to trace if and to what extent prolonged exposure to air pollution, in terms of peaks of concentration of fine dust and other pollutants, constitutes a detrimental factor in COVID-19 [28] cases, we paid particular attention to the relationship between climate and air quality [29].
All anthropic activities generate emission of gaseous and particulate pollutants that modify the composition of the atmosphere. Air quality and climate change are two closely related environmental issues [30]. Climate change affects the atmospheric processes and causes changes in the functioning of terrestrial and marine ecosystems which can, in turn, affect the atmospheric processes [31]. However, these two environmental emergencies are still considered separately, both at the level of the scientific community and from those responsible for environmental policies, as in the case of the recent COVID-19 emergency [28]. For this reason, policies to improve air quality and mitigate climate change must inevitably be integrated. These options favor one of the two aspects and worsen the situation of the other (win-lose policies). Coordinated actions that take due account of the connections between air quality and climate change constitute the best strategy in terms of economic and social costs (win-win policies) [32]. According to the EEA-European Environmental Agency, although air pollution [33] affects the whole population (collective health costs), only a part is more exposed [34] to individual risks [35][36][37].
In this sense, although the air pollutant containment measures [38] that derive from the numerous urban initiatives (smart building, mobility and industry 4.0) are particularly important, the megatrend of the globalization of industrial and agricultural production with related post-industrial lifestyles shows that it is not aligned with the green energy production, circular economy and ecosystem services. In particular, in the Greater Milan metropolitan area, Po Valley represents the outcome of industrial and agricultural globalization in Italy characterized by an increasingly critical quality of the air [25]. Although in the last decade in Italy there have been tax incentive measures for the purchase or improvement of the ecological performance of home-heating [39] and public and private road vehicles, the levels of air pollution for 150 days (2018) [40] exceeded the EU regulatory limits, which are much Sustainability 2020, 12, 5064 6 of 44 lower than the WHO guidelines. This situation has been protracted, as high levels of air pollution and the concentration of pollutants in the air have been constantly reported in recent years [40]. In addition, it is also necessary to remember how the climatic and geographical 'handicap effect' of the Greater Milan metropolitan area is not secondary in the quality of the air. In summary, the air pollution in big urban agglomerations [41][42][43] such as the Greater Milan metropolitan area contribute to climatic variations although there are many synergies and points of conflict between air quality and climate and urban territorial policies.

Land Take
The environmental pressure is characterized not only by atmospheric emissions, but also land taken by human activities is important in terms of the impact on the environment. As in Section 2.1.2, land take is a relevant and spreading phenomenon, that is not only related to the competition for land among different sectors and activities, but deals also with the capacity of the environment to sustain the outcome of human activities. Specifically, land take can be considered as the increase in the amount of agriculture, forest, and other semi-natural and natural land taken by urban and other artificial land development. It includes areas sealed by construction and urban infrastructure as well as urban green areas and sport and leisure facilities. The main drivers of land take are grouped in processes, resulting in the extension of housing, services, and recreation, industrial, and commercial sites, transport networks and infrastructures mines, as well as quarries and waste dumpsites. In the literature examined in this paper, a similar concept encountered is soil consumption (ISPRA).
In the present research, we were interested in understanding the causes and origins of such a violent and massive COVID-19 outbreak in Northern Italy and what could be the special-and spatial -factors reinforcing this. To accomplish this, we considered COVID-19 mortality as the most usable data for giving an idea of the gravity and strength of COVID-19 in different parts of Italy. We analyzed such aspects related to COVID-19 by means of their relation and consideration of the spatial proximity, by means of LISA (local indicators of spatial autocorrelation). Such aspects were compared with a selection of variables from a wide database of economic, geographical, climatic and environmental elements we collected during the research. In particular, the subset of variables in four major groups, those in our view more likely to be related to the COVID-19 outbreak, as Land use, air quality, climate and weather, population, health and life expectancy. Our study, which does not claim to be exhaustive but only intends to demonstrate the initial results, develops in a complex international framework.
The rest of the paper is organized as it follows. Section 2 deals with Materials, Data and Methods, with an overview, in Section 2.1, of the study area, the issues of land take phenomenon and the basics of spatial diffusion process related to pandemics, followed by a broad analysis on the data collected and prepared for the analysis in Section 2.2. We present Methods in Section 2.3, with particular reference to the ecological approach adopted and an explanation on LISA-Local Indicators of Spatial Autocorrelation. Results follow in Section 3, where the analysis of spatial autocorrelation is performed over the main groups of variables considered in the research. Section 4 presents a discussion concerning the local and global geographical aspects of the COVID-19 and issues related to pollutants, land take, and some suggestions for policies. Finally, in Section 5 are the conclusions and future development of the research. Most of the population is concentrated in the Po Valley geographical region, surrounded by the Alpine and Apennine mountains and by the Adriatic Sea, eastwards towards the Po Delta. The Po Valley alone represents Italy's economic 'core'. In an area of approximately 55,000 km 2 nearly 22 million people live, with a density (400 inhabitants per km 2 ) double that of the rest of the peninsula, reaching different peaks in the main urban areas of the Greater Milan metropolitan area (the neighboring Milan and Monza Provinces exceed 2000 inhabitants per square kilometer). It is noticeable that in 2015 the former 'Province of Milan' become the Metropolitan City of Milan, actually covering the same area and therefore the same municipalities of the former qualification. From a functional point of view, as generally happens with metropolitan areas worldwide, the 'Greater Milan Metropolitan', can be considered a wider area, covering Milan's neighboring provinces (Varese, Como, Lecco, Monza-Brianza, Pavia, Lodi and Cremona). Extended interpretations of the concept involve also other provinces such as Bergamo, in Brescia, Lombardy, and others belonging administratively to other regions, like Novara (hosting Milan-Malpensa International Airport), Alessandria in Piedmont, and Piacenza in Emilia Romagna (Figure 3).  The Milan area is the offset center of the Po Valley Megalopolis, an urbanized and industrialized area including the major Northern Italian cities, from the Western former 'industrial triangle' of Milan, Genoa and Turin, and moving eastwards towards Venice and South East to Bologna and beyond, on the Adriatic Coast [44]. A high level of mobility is present in the area, both within the same area, as commuters crowd roads and railway lines, and in terms of the national and international connections, due to the presence of international airports, major motorways and HSR with fast services towards major national urban destinations. Such an area is also the Southern part of the image of the 'Blue Banana' [45]: such area is known as the Po Valley Megalopolis, and part of a wider European Megalopolis, gathering the major metropolitan areas and major cities in Europe. This is an area that stretches in the fruit-like shape from Northern Italy through Germany, North Eastern France, Benelux, and Southern England, as the urban, economic and industrial core of Europe (References on Italian landscape and Po Valley area: [46][47][48][49][50][51][52][53]. In this complex of urban and territorial dynamics of the Po Valley Megalopolis, it is possible to distinguish a geo-spatial correspondence between population density and emissions of NO2, which compromises air quality (Figure 4; Figure  A1). The Milan area is the offset center of the Po Valley Megalopolis, an urbanized and industrialized area including the major Northern Italian cities, from the Western former 'industrial triangle' of Milan, Genoa and Turin, and moving eastwards towards Venice and South East to Bologna and beyond, on the Adriatic Coast [44]. A high level of mobility is present in the area, both within the same area, as commuters crowd roads and railway lines, and in terms of the national and international connections, due to the presence of international airports, major motorways and HSR with fast services towards major national urban destinations. Such an area is also the Southern part of the image of the 'Blue Banana' [45]: such area is known as the Po Valley Megalopolis, and part of a wider European Megalopolis, gathering the major metropolitan areas and major cities in Europe. This is an area that stretches in the fruit-like shape from Northern Italy through Germany, North Eastern France, Benelux, and Southern England, as the urban, economic and industrial core of Europe (References on Italian landscape and Po Valley area: [46][47][48][49][50][51][52][53]. In this complex of urban and territorial dynamics of the Po Valley Megalopolis, it is possible to distinguish a geo-spatial correspondence between population density and emissions of NO 2 , which compromises air quality (Figure 4; Figure A1).  However, the issue is more complex as the variables involved are manifold and aspects such as the scale of investigation (wind, relative humidity, more generally weather and climate conditions) and all air pollutants must be considered, for example, the PM10 which during the lockdown exceeded the limits in numerous units of the Po Valley megalopolis [54,55].

Land Take Phenomenon in Italy
In the last fifty years, the phenomenon of land take [56] heavily occurred in Italy in different forms and in various areas [57][58][59]24]. In particular, with proximity to metropolitan and production areas, the phenomenon is more intense and takes the form of urban sprawl: new low density settlements, with poor services and connections, close to the city, where people have the feeling of living in a more natural context. The main driving forces of this trend are the spatial configuration and the appeal of the areas. When space is characterized by a homogeneous and isotropic form, the phenomenon of sprawl is more elevated. Furthermore, socio-economic indicators play an important role in attracting people, such as new residents or commuters, generating the building of new neighborhoods or transport infrastructures. Following this model of urban growth, the soil loses its biological value, becoming unable to absorb and filter rainwater, producing negative effects on biodiversity as well as on agricultural production [60][61][62]. This sealing process produces the loss of natural ecosystem functions generating complete soil degradation [63][64][65][66]. Soil is where energy and substances exchange with other environmental elements. The role of soil in the hydrogeological cycle is very important. Solar radiation causes the evaporation of water from accumulation areas towards the atmosphere. The steam rises to a high altitude, cools and condenses forming clouds. The water then returns to the emerged lands in form of precipitation, part of it falls into the rivers and the surface water network, another is absorbed by the soil reaching the groundwater. The soil controls the flow of surface water and regulates its absorption by filtering polluting substances. Infiltration also depends on the permeability and porosity of soil.
The soil is also fundamental in the carbon cycle. Carbon is everywhere in nature and it is transformed into oxygen through photosynthesis in the carbon cycle. Through the plants, soil absorbs carbon dioxide, which can remain underground for thousands of years, feeding the soil microorganisms. Consequently, the soil is a sort of absorption well where CO2 sequestration and storage are possible. Poor soil management can generate a loss of these properties, thereby producing negative effects [67]. Soil and related ecosystem services are important elements in the improvement of air quality reducing PM10 and O3 [68,69].
Another model of land take occurred in Italy in more remote zones that were less accessible and mostly located in mountainous or hilly areas [70]. This uncontrolled urban development has been defined as sprinkling [71,72]. Unlike urban sprawl, sprinkling is characterized by low density settlements and a more spontaneous, dispersed and chaotic development. The cost of sprinkling is higher than that of sprawl [73,74] because this model requires a lot of small infrastructures for transport, electricity, water distribution etc. If, on the one hand, sprawl is denser than sprinkling However, the issue is more complex as the variables involved are manifold and aspects such as the scale of investigation (wind, relative humidity, more generally weather and climate conditions) and all air pollutants must be considered, for example, the PM 10 which during the lockdown exceeded the limits in numerous units of the Po Valley megalopolis [54,55].

Land Take Phenomenon in Italy
In the last fifty years, the phenomenon of land take [56] heavily occurred in Italy in different forms and in various areas [24,[57][58][59]. In particular, with proximity to metropolitan and production areas, the phenomenon is more intense and takes the form of urban sprawl: new low density settlements, with poor services and connections, close to the city, where people have the feeling of living in a more natural context. The main driving forces of this trend are the spatial configuration and the appeal of the areas. When space is characterized by a homogeneous and isotropic form, the phenomenon of sprawl is more elevated. Furthermore, socio-economic indicators play an important role in attracting people, such as new residents or commuters, generating the building of new neighborhoods or transport infrastructures. Following this model of urban growth, the soil loses its biological value, becoming unable to absorb and filter rainwater, producing negative effects on biodiversity as well as on agricultural production [60][61][62]. This sealing process produces the loss of natural ecosystem functions generating complete soil degradation [63][64][65][66]. Soil is where energy and substances exchange with other environmental elements. The role of soil in the hydrogeological cycle is very important. Solar radiation causes the evaporation of water from accumulation areas towards the atmosphere. The steam rises to a high altitude, cools and condenses forming clouds. The water then returns to the emerged lands in form of precipitation, part of it falls into the rivers and the surface water network, another is absorbed by the soil reaching the groundwater. The soil controls the flow of surface water and regulates its absorption by filtering polluting substances. Infiltration also depends on the permeability and porosity of soil.
The soil is also fundamental in the carbon cycle. Carbon is everywhere in nature and it is transformed into oxygen through photosynthesis in the carbon cycle. Through the plants, soil absorbs carbon dioxide, which can remain underground for thousands of years, feeding the soil microorganisms. Consequently, the soil is a sort of absorption well where CO 2 sequestration and storage are possible. Poor soil management can generate a loss of these properties, thereby producing negative effects [67]. Soil and related ecosystem services are important elements in the improvement of air quality reducing PM 10 and O 3 [68,69].
Another model of land take occurred in Italy in more remote zones that were less accessible and mostly located in mountainous or hilly areas [70]. This uncontrolled urban development has been defined as sprinkling [71,72]. Unlike urban sprawl, sprinkling is characterized by low density settlements and a more spontaneous, dispersed and chaotic development. The cost of sprinkling is higher than that of sprawl [73,74] because this model requires a lot of small infrastructures for transport, electricity, water distribution etc. If, on the one hand, sprawl is denser than sprinkling educing a depauperation of the landscape, on the other, this development creates more sealing of soils, generating the total loss of their natural properties.

Geography of Diffusion
Diffusion Processes in Geography: Some Theories A virus outbreak is a typical, still dramatic and frightening case of spatial diffusion, a topic that is well-known and studied in geography and its fundamentals. Diffusion in geography implies the movement of an event or set of events in space and time, and brings with it the idea of a process and the drawing of a pattern, as the outcome of the event's movement in space and time [75][76][77][78]. Diffusion has been studied in geography with reference to very different sets of cases and situations, from plagues to financial crises, from migration to music styles, from physical geography to human and economic geography. The analysis of these phenomena, performed by authors in several different contexts, introduced some basic elements to be reprised.
A first classification of spatial diffusion can distinguish cases between relocation and expansion. Relocation implies the physical movement and abandonment of the site of origin of the event, towards a new one. Expansion implies the spatial and temporal extension of a given state, or event, to cover and fill (all of) the available space ( Figure 5).
Sustainability 2020, 12, x FOR PEER REVIEW 10 of 48 educing a depauperation of the landscape, on the other, this development creates more sealing of soils, generating the total loss of their natural properties.

Geography of Diffusion
(1) Diffusion Processes in Geography: Some Theories A virus outbreak is a typical, still dramatic and frightening case of spatial diffusion, a topic that is well-known and studied in geography and its fundamentals. Diffusion in geography implies the movement of an event or set of events in space and time, and brings with it the idea of a process and the drawing of a pattern, as the outcome of the event's movement in space and time [75][76][77][78]. Diffusion has been studied in geography with reference to very different sets of cases and situations, from plagues to financial crises, from migration to music styles, from physical geography to human and economic geography. The analysis of these phenomena, performed by authors in several different contexts, introduced some basic elements to be reprised.
A first classification of spatial diffusion can distinguish cases between relocation and expansion. Relocation implies the physical movement and abandonment of the site of origin of the event, towards a new one. Expansion implies the spatial and temporal extension of a given state, or event, to cover and fill (all of) the available space ( Figure 5). The expansion process can proceed in different ways and formats, and follows different rules: contagion, network, hierarchical, waterfall. Contagion is the typical 'local' process that implies a contact between the event carrying the 'innovation' and those not yet affected. The network diffusion process deals with the network structure of contact between subjects (involved in the diffusion) at local and global levels. It also implies the existence of social networks between people, both locally and globally and involves the presence of major transport infrastructures and networks, i.e., local transit systems and global major air transport routes [79]. The hierarchical expansion process comes about when innovation is spread by means of privileged channels of communication and between centers of higher importance. Major transport and communication routes help in channeling the spread of the innovation in space and time. Waterfall implies the direction and speed of the diffusion in the following manner: it is generally fast, following a top-down approach , i.e., from major centers to minor ones and slow when it moves from minor centers towards major ones, in a bottom-up approach ( Figure 6). The expansion process can proceed in different ways and formats, and follows different rules: contagion, network, hierarchical, waterfall. Contagion is the typical 'local' process that implies a contact between the event carrying the 'innovation' and those not yet affected. The network diffusion process deals with the network structure of contact between subjects (involved in the diffusion) at local and global levels. It also implies the existence of social networks between people, both locally and globally and involves the presence of major transport infrastructures and networks, i.e., local transit systems and global major air transport routes [79]. The hierarchical expansion process comes about when innovation is spread by means of privileged channels of communication and between centers of higher importance. Major transport and communication routes help in channeling the spread of the innovation in space and time. Waterfall implies the direction and speed of the diffusion in the following manner: it is generally fast, following a top-down approach, i.e., from major centers to minor ones and slow when it moves from minor centers towards major ones, in a bottom-up approach ( Figure 6). Obviously, when the innovation reaches a higher centre, it will start the fast top-down diffusion. (Figure 7). (2) COVID-19 Theory and Practice The virus and disease diffusion possibly act according to a mix of the above-mentioned ways of diffusion. Furthermore, Haggett and Cliff [80] recall that diffusion processes come about as spatial diffusion waves, starting either in a single or a set of locations and then spreading out through different processes, covering wider areas. Haggett and Cliff [81][82][83] along with Smallman-Raynor [84,83] [85] among other geographers have modelled such diffusion modes, examining also the relation of epidemics in space and time and the wave-nature of epidemics.
The diffusion process is a mix of expansion and relocation processes: usually an epidemic starts in a given region expanding in space, and relocation occurs when its footprint fades in the place of origin and continues to grow in newly affected areas. Therefore, the process is contagious, when the virus spreads through direct contacts; network, when it follows networks of relations and flows between individuals and places; hierarchical, because major centres affect a higher number of lower order centres; and waterfall because the movement is generally stronger between centres following a top-down approach. The 'wave' also reverses its direction as the population recovers and early infected regions go back to a cleared situation [86,87]. In geographical terms, a diffusion wave follows, generally and theoretically, in 5 steps, (Figure 8). Obviously, when the innovation reaches a higher centre, it will start the fast top-down diffusion. (Figure 7). Obviously, when the innovation reaches a higher centre, it will start the fast top-down diffusion. (Figure 7). (2) COVID-19 Theory and Practice The virus and disease diffusion possibly act according to a mix of the above-mentioned ways of diffusion. Furthermore, Haggett and Cliff [80] recall that diffusion processes come about as spatial diffusion waves, starting either in a single or a set of locations and then spreading out through different processes, covering wider areas. Haggett and Cliff [81][82][83] along with Smallman-Raynor [84,83] [85] among other geographers have modelled such diffusion modes, examining also the relation of epidemics in space and time and the wave-nature of epidemics.
The diffusion process is a mix of expansion and relocation processes: usually an epidemic starts in a given region expanding in space, and relocation occurs when its footprint fades in the place of origin and continues to grow in newly affected areas. Therefore, the process is contagious, when the virus spreads through direct contacts; network, when it follows networks of relations and flows between individuals and places; hierarchical, because major centres affect a higher number of lower order centres; and waterfall because the movement is generally stronger between centres following a top-down approach. The 'wave' also reverses its direction as the population recovers and early infected regions go back to a cleared situation [86,87]. In geographical terms, a diffusion wave follows, generally and theoretically, in 5 steps, (Figure 8).

COVID-19 Theory and Practice
The virus and disease diffusion possibly act according to a mix of the above-mentioned ways of diffusion. Furthermore, Haggett and Cliff [80] recall that diffusion processes come about as spatial diffusion waves, starting either in a single or a set of locations and then spreading out through different processes, covering wider areas. Haggett and Cliff [81][82][83] along with Smallman-Raynor [83][84][85] among other geographers have modelled such diffusion modes, examining also the relation of epidemics in space and time and the wave-nature of epidemics.
The diffusion process is a mix of expansion and relocation processes: usually an epidemic starts in a given region expanding in space, and relocation occurs when its footprint fades in the place of origin and continues to grow in newly affected areas. Therefore, the process is contagious, when the virus spreads through direct contacts; network, when it follows networks of relations and flows between individuals and places; hierarchical, because major centres affect a higher number of lower order centres; and waterfall because the movement is generally stronger between centres following a top-down approach. The 'wave' also reverses its direction as the population recovers and early infected regions go back to a cleared situation [86,87]. In geographical terms, a diffusion wave follows, generally and theoretically, in 5 steps, (Figure 8). Sustainability 2020, 12, x FOR PEER REVIEW 12 of 48 Evidences from past outbreaks lead to highlighting both local diffusion (contagion) and longrange ones (hierarchical, cascade). C. Maturity: The highest intensity is reached with clusters spread over a vulnerable populationthe entire areas involved in the epidemic. The intensity is at its maximum with contrasts in infection density in different sub-regions. D. Decay: Fewer reported cases and decline is registered, with a slower spatial contraction than the proper diffusion steps. Low intensity infected areas appear as scattered. E. Extinction: The tail of the epidemic wave can be spotted through few and scattered cases, that can be found mostly in less accessible areas.
The data available so far related to the COVID- 19 virus is yet to be fully validated and understood and, at international and national levels allows only for a limited possibility of analysis and understanding of the phenomena at stake and under the process of change in time and space. From beginning of the outbreak, we can observe that the contagion, as a diffusion process characterizes the local ones, that can be observed both in the place of origin, i.e., Wuhan city and Hubei region, as well as taking place at the different local scales, i.e., South Korea and other neighboring countries at the beginning of the phenomenon, and Italy and other Western countries in the subsequent stages. With medium and long distances, but also shorter ones, as will be seen in the next examples, diffusion is based on long transport networks, namely rail (high speed trains) and air routes, and on shorter ones, such as rail (regional) and maritime (ferry) routes.
While the diffusion that occurs locally follows a contagion model, a hierarchical diffusion is responsible for the regional and international diffusion of the disease. Not surprisingly, the neighboring South Korea was the first major country affected, followed, after some weeks, by a western country, i.e., Italy, and then other European and North American countries the weeks after that. Recent studies by Tatem et al. [88], or Ben-Zion et al. [89], Bowen and Laroe [90] and, particularly, that of Brokmann and Helbing [91] on the SARS and swine flu diseases, show the network geometry of the transport system, namely air transport, as the backbone for human interactions at a global scale and also as the backbone for a virus outbreak diffusion outside its region of origin. The simulations presented by the authors help in understanding and explaining the geographical outbreaks of the major epidemics that took place during the first decade of this century and provide an estimate of the possible evolution in terms of the temporal order of the areas hit. By means of an example, we can theoretically dispute what has been presented by the abovementioned scholars, that after an outbreak in China where infectious diseases of different kinds have often originated in the past, even if at a slower pace, from the 1960s and characterized particularly as respiratory syndromes [92], major air connections facilitated spreading on the mainland and towards neighboring countries such as South Korea and Japan, Europe and the United States, to cite a few examples of major destination areas. Major air transport routes to and from China connect European destinations-count for 9.8% of EU air traffic, with Amsterdam Schiphol, Frankfurt, London and Evidences from past outbreaks lead to highlighting both local diffusion (contagion) and long-range ones (hierarchical, cascade). C. Maturity: The highest intensity is reached with clusters spread over a vulnerable population-the entire areas involved in the epidemic. The intensity is at its maximum with contrasts in infection density in different sub-regions. D. Decay: Fewer reported cases and decline is registered, with a slower spatial contraction than the proper diffusion steps. Low intensity infected areas appear as scattered. E. Extinction: The tail of the epidemic wave can be spotted through few and scattered cases, that can be found mostly in less accessible areas.
The data available so far related to the COVID- 19 virus is yet to be fully validated and understood and, at international and national levels allows only for a limited possibility of analysis and understanding of the phenomena at stake and under the process of change in time and space. From beginning of the outbreak, we can observe that the contagion, as a diffusion process characterizes the local ones, that can be observed both in the place of origin, i.e., Wuhan city and Hubei region, as well as taking place at the different local scales, i.e., South Korea and other neighboring countries at the beginning of the phenomenon, and Italy and other Western countries in the subsequent stages. With medium and long distances, but also shorter ones, as will be seen in the next examples, diffusion is based on long transport networks, namely rail (high speed trains) and air routes, and on shorter ones, such as rail (regional) and maritime (ferry) routes.
While the diffusion that occurs locally follows a contagion model, a hierarchical diffusion is responsible for the regional and international diffusion of the disease. Not surprisingly, the neighboring South Korea was the first major country affected, followed, after some weeks, by a western country, i.e., Italy, and then other European and North American countries the weeks after that. Recent studies by Tatem et al. [88], or Ben-Zion et al. [89], Bowen and Laroe [90] and, particularly, that of Brokmann and Helbing [91] on the SARS and swine flu diseases, show the network geometry of the transport system, namely air transport, as the backbone for human interactions at a global scale and also as the backbone for a virus outbreak diffusion outside its region of origin. The simulations presented by the authors help in understanding and explaining the geographical outbreaks of the major epidemics that took place during the first decade of this century and provide an estimate of the possible evolution in terms of the temporal order of the areas hit. By means of an example, we can theoretically dispute what has been presented by the abovementioned scholars, that after an outbreak in China where infectious diseases of different kinds have often originated in the past, even if at a slower pace, from the 1960s and characterized particularly as respiratory syndromes [92], major air connections facilitated spreading on the mainland and towards neighboring countries such as South Korea and Japan, Europe and the United States, to cite a few examples of major destination areas. Major air transport routes to and from China connect European destinations-count for 9.8% of EU air traffic, with Amsterdam Schiphol, Frankfurt, London and Paris among the major airports with the highest number of international connections (Rome Fiumicino Airport, however, also offered a direct connection to Wuhan Tianhe International airport) [93][94][95].
Recent studies seem to show that patient zero in Europe, although asymptomatic, was identified in Germany, in January [12]. Furthermore, a particularly high seasonal flu peak was registered in Germany in the early weeks of the year 2020 [96]. The outbreak conditions, however, were found in Italy, that as a consequence was hit first and in a very aggressive way.
Diffusion as a Local Spatial Process Before and after the Italian Lockdown With reference to diffusion processes at a regional/Italian level, we can recall some possible diffusion dynamics that occurred after the first two outbreaks in Vo (Veneto) and Codogno (Lombardy) around 20th February 2020. The local diffusion processes led to Lombardy and other provinces in Piedmont, Emilia-Romagna and Veneto being inserted in a Red Zone at the beginning of March (8th March 2020), shortly before the decision of placing the whole country into a single red zone (10th March 2020), in a lock down with severe limitations to individual movement and industrial production. Italian internal migrations still see many people from Southern Italy moving to the cities in the north for work purposes. Furthermore, many southern Italian locations host holiday homes for tourism. Journeys back home and towards holiday homes, moved people from the north to the south in the days before the full lockdown became fully operational. Fear arose of a spatial diffusion of the virus towards southern regions via long-distance means of transport (high speed trains; air connections).
The containment policies evolved from mild forms concerning early affected provinces to stricter ones limiting mobility and actions. The picture we drew up at the end of March 2020 can be considered a reasonable representation of the 'natural' diffusion of the phenomenon, before the effects of the lock down policies and given the two-week maximum incubation time [97], although such a duration can change in time and space [98][99][100].

The Data
Research was carried out using different datasets mainly referred to Italy and related to the COVID-19 outbreak, as well as socio-economic and environmental data, considered useful for the examination of the territorial aspects of the virus outbreak in Italy.
We made a preparatory work of selection of variables, of which only a part was used in the present research, while other are currently undergoing further analyses in other research lines.
COVID-19 data considered the number of total infected people as at 31 March 2019 at a provincial level, as reported by the Italian Ministry of Health, and as collected by the Civil Protection. It is to be outlined that such data is considered as raw data, including those cases were the virus was not validated by the Italian Higher Institute of Health (ISS-Istituto Superiore di Sanità) that checks the actual cases, even after a one-to-one evaluation of the causes of death. We considered data as at that day to 'close an important month' in terms of the virus outbreak, and in order to have a picture of the situation after the austere choices of locking down, on a national level, most of the activities and individual mobility. Data after that moment was difficult to relate to diffusion processes in strict terms and more related to the regional policies taken after the national lockdown.
An important novel dataset, originally built from scratch by the research group, is the number of deaths at a provincial level. Such data was collected from different sources, reason being that such data was not always available from the same sources. In many cases data was provided by regional administrations, while in other cases the research required counting and referring data to the provinces from the local health agencies, or even from other sources such as newspapers that provided data locally. Among others, major difficulties were found in locating data at a provincial level for important regions in terms of the COVID-19 outbreak, such as Lombardy and Piedmont as well as the big regions like Liguria, Lazio, Campania, and Sicily, which required an extra effort to tracing the deaths at a provincial level. We succeeded in localizing, provincially, 11,336 deaths over 12,428 nationwide, and 102,440 infected people over 105,792.
The socio-economic and environmental data taken into consideration comes from different official sources. Socio-economic and demographic data comes from the ISTAT (Italian Statistical Institute-Istituto Nazionale di Statistica), such as population (total and organized in age groups), as well as mortality, differentiated by causes, as at 2019.
Environmental data and indicators come from ISPRA (Higher Institute for Environmental Protection and Research-Istituto superiore per la protezione e la ricerca ambientale), WHO (World Health Organization), ISS, EEA (European Environmental Agency), Il Sole 24 Ore (economic and business newspaper, which provides constant reports on economical facts), Legambiente (no-profit association for environmental protection), ACI (Italian Automobile Club), ilmeteo.com and windfinder.com (weather and wind data). Furthermore, the air quality (PM 2.5 , PM 10 , NH 3 , CO, CO 2 , NO x ) and the weather conditions (humidity, wind, rain) were also monitored promptly through specific dashboards. In this sense, we have elaborated Figure 9 that summarizes the data set (dynamic and static) in reference to the ecological approach used. The complete data set (open data) is shown at the end of the paper (The full dataset is resumed in Supplementary Table S1; the subset, as the dataset used for the analysis is presented in Appendix A, Table A1). The socio-economic and environmental data taken into consideration comes from different official sources. Socio-economic and demographic data comes from the ISTAT (Italian Statistical Institute-Istituto Nazionale di Statistica), such as population (total and organized in age groups), as well as mortality, differentiated by causes, as at 2019.
Environmental data and indicators come from ISPRA (Higher Institute for Environmental Protection and Research-Istituto superiore per la protezione e la ricerca ambientale), WHO (World Health Organization), ISS, EEA (European Environmental Agency), Il Sole 24 Ore (economic and business newspaper, which provides constant reports on economical facts), Legambiente (no-profit association for environmental protection), ACI (Italian Automobile Club), ilmeteo.com and windfinder.com (weather and wind data). Furthermore, the air quality (PM2.5, PM10, NH3, CO, CO2, NOx) and the weather conditions (humidity, wind, rain) were also monitored promptly through specific dashboards. In this sense, we have elaborated Figure 9 that summarizes the data set (dynamic and static) in reference to the ecological approach used. The complete data set (open data) is shown at the end of the paper (The full dataset is resumed in Supplementary Table S1; the subset, as the dataset used for the analysis is presented in Appendix A, Table A1). The spatial units selected are those once known as Provinces, the intermediate levels between Municipalities and Regions, still used as statistical units also by ISTAT, as many of them lose or change their administrative role for a total of 107 units. Spatial units are provided by ISTAT as at 2019. The choice of these units, from the geographical, cartographical and spatial analytical points of view, holds several limitations. This happens as the phenomena referring to such units tend to be diluted over irregular and non-homogeneous areas, in terms of both spatial dimension, number of inhabitants, and population densities, while incorporating several differences, not only among each other, but in terms of the spatial variation within the same area. The risk, as is often outlined, is confusing the spatial pattern drawn by the geographical units as well as that of the underlying population, instead of the phenomena per se [101][102][103][104]. Such issues are however well known. The The spatial units selected are those once known as Provinces, the intermediate levels between Municipalities and Regions, still used as statistical units also by ISTAT, as many of them lose or change their administrative role for a total of 107 units. Spatial units are provided by ISTAT as at 2019. The choice of these units, from the geographical, cartographical and spatial analytical points of view, holds several limitations. This happens as the phenomena referring to such units tend to be diluted over irregular and non-homogeneous areas, in terms of both spatial dimension, number of inhabitants, and population densities, while incorporating several differences, not only among each other, but in terms of the spatial variation within the same area. The risk, as is often outlined, is confusing the spatial pattern drawn by the geographical units as well as that of the underlying population, instead of the phenomena per se [101][102][103][104]. Such issues are however well known. The choice of provinces as spatial units, instead of regions, was the only one to allow a finer and disaggregate analysis at a local level. It has to be pointed out however, that with particular reference to obvious areas at risk, the Po Valley area holds some characteristics of an area that is nearly homogeneous and isotropic as per Christaller's terms, because the provinces in part of Piedmont, Lombardy, Veneto, and Emilia Romagna present quite comparable spatial dimensions.
Moreover, as will be more evident in the rest of the paper, the consideration of neighboring provinces will help to understand the spatial variation of the phenomena that are not dependent from the regional differences.
Similarly, that would have also needed reasoning at a nearly 'international' level, such as the high level of autonomy granted in Italy to the 19 Regions and 2 Autonomous Provinces of Trento and Bolzano, particularly in terms of Health system and Mobility, would suggest approaching the issue as a comparison between independent states. In this sense, the retrieval of data, specifically open data, the cataloging, representation and geospatial correlation thereof, has always been consistent with the ecological approach, specifically in order to evaluate the phenomena in their complexity and entirety.

The Ecological Approach
In this paper the ecological approach was adopted, because the physiological traits of the virus are combined with a wide set of selected relevant environmental variables. In this sense, occurrences of the outbreak (as infected cases and deaths), were examined and referred to several variables. To do that, throughout the paper, a regional focus on the characters of the study area-the Po Valley in the Italian context, from different points of view-is applied, in terms of the physical and human-economic geographical characters, describing them and observing them from a qualitative and quantitative point of view. The Po Valley area analysis-especially the Greater Milan metropolitan area-is performed, recalling the similarities with Hubei (Wuhan areas) in order to explore the potential analogies with the conditions of COVID-19 outbreaks. For this to come about, we concentrated on particular elements related to aspects that, in an integrated manner, can be considered important in understanding the human-environmental relations between human activities, geographic and climatic conditions and virus outbreaks. Focus on air and climate characteristics and on soil consumption, being the common features concerning several human-related behaviors that affect the environmental balance, is also applied. In fact, as soil consumption increases, the storage capacity for air purification decreases and more generally a deterioration of the physical, chemical, biological and economic characteristics connected to it [105]. In other words, land take interferes with collective well-being [106]. In this sense, ecosystem services play a key role for the collective well-being, namely ecological, social and productive [107,108]. Our research effort is part of this interdisciplinary ecological approach to investigate the spread of COVID-19 in Italy, based on a theoretical and quantitative analysis on a large dataset of environmental variables.

Calculation of Case Fatality Rate
The absolute frequencies of deaths per province, obtained through the reconstruction of data from distinct information sources, have been compared with the number of positive COVID-19 cases published by the ISS [6], according to the formula: Standardized Mortality Ratio is a standardization method used to make comparisons of death rates between different regions, considering that a given region can have a population that is older than another, and considering that younger people are less likely to die than older people. To do that, it was necessary to investigate the pattern of deaths and the dependency on the age composition. For each areal unit, based on the distribution of population per age group, and the age-specific rates of deaths in some larger populations, the expectation of number of deaths was calculated. The ratio of observed over expected deaths was calculated. A value of 1 indicates that the area considered was reacting 'as expected' in terms of mortality, in line with that of a wider reference area. Values higher than 1 show a mortality rate that was higher than expected, even in terms of population structure, while values lower than 1 suggested mortality was reduced and lower than expected [109].
Specific mortality from COVID-19 has been standardized for each Italian province and per age groups (10 groups) where the first group was 0-9 years and the last group 90-∞, with reference to the national population figures in the year 2019 [110]. The indirect standardization process initially provided for the calculation of national specific mortality by age group, obtained by dividing the number of COVID-19 deaths confirmed by ISS [6] with the 10 defined age groups. Thus, the number of deaths expected in the Italian provinces for the age groups previously identified and based on the 2019 provincial populations [110], was calculated according to the formula: where n i is the specific age group population in each observed area (province); R i is the national mortality rate for the specific age group. The standardized mortality ratio (SMR) was obtained by comparing the number of events observed in each province with the respective number of expected events: where d is the number of observed deaths; e the number of expected deaths. Finally, the 95% confidence intervals (95% CI) were calculated as proposed by Vandenbroucke [111].

Spatial Autocorrelation
When dealing with data referred to spatial units, the characters involved are locational information and the properties or 'attribute' data. The pattern drawn by geographical features and the data referred thereto can be mutually influenced, resulting in spatial autocorrelation. In geographical analytical terms, it is the capability of analyzing locational and attribute information at the same time [112]. The study of spatial autocorrelation can be very effective in analyzing the spatial distribution of objects, examining simultaneously the influence of neighboring objects, in a concept anticipated by Waldo Tobler [113] in the first law of geography, which states "All things are related, but nearby things are more related than distant things". Although this approach is very simple and intuitive [114] and very important in a huge variety of application domains, it has not been applied for more than twenty years [115].

1.
i and j are two objects;

2.
N is the number of objects; 3.
C ij is a degree of similarity of attributes i and j; 4.
W ij is a degree of similarity of location i and j; This general formula gave origin to two indices widely used in spatial analysis, namely the Geary C Ration [117] and the Moran Index I [118].
Defining x i as the value of object i attribute; if c ij = (x i − x j ) 2 , Geary C Ratio can be defined as follows: Index I can be defined as follows: As recalled in Murgante and Borruso [119], these indices are quite similar, with the main difference in the cross-product term in the numerator, calculated using the deviations from the mean in Moran, while in Geary it is directly computed.
The indices are useful for highlighting the presence-or absence-of spatial autocorrelation at a global level in the overall distribution, while the local presence of autocorrelation can be highlighted by the so-called LISA (Local Indicators of Spatial Association). Anselin [120,121] considered LISA as a local Moran index. The sum of all local indices is proportional to the value of the Moran one: The index is calculated as follows: For each location it allows to assess the similarity of each observation with its surrounding elements. Five scenarios emerge: • locations with high values of the phenomenon and a high level of similarity with its surroundings (high-high H-H), defined as hot spots; • locations with low values of the phenomenon and a low level of similarity with its surroundings (low-low L-L), defined as cold spots; • locations with high values of the phenomenon and a low level of similarity with its surroundings (high-low H-L), defined as potentially spatial outliers; • locations with low values of the phenomenon and a high level of similarity with its surroundings (low-high L-H), defined as potentially spatial outliers; • locations completely lacking of significant autocorrelations.
The interesting characteristic of LISA is in providing an effective measure of the degree of relative spatial association between each territorial unit and its surrounding elements, thereby highlighting the type of spatial concentration and clustering.
An important element to be considered in the above-mentioned equations is the parameter weight, w ij , related to the neighborhood property. Generally, w ij [122] values indicate the presence, or absence, of neighboring spatial units to a given one. A spatial weight matrix is realized, with w ij assuming values of 0 in cases in which i and j are not neighbors, or 1 when i and j are neighbors. Neighborhood is computed in terms of contiguity such as, in the case of areal units, sharing a common border of non-zero length [104]. If one were to adopt a chess game metaphor [104], contiguity can be considered as allowed along the paths of rook, bishop and queen.

Results
Through the evaluation of the elaborations from the data collected for this research, according to an interdisciplinary ecological approach and referring to different scales from global to provincial, we have been able to obtain some important results, which we believe will also be developed in a subsequent research.

Mortality rates, SMRs, Population Density, Commuting
The death toll at a provincial level required an in-depth search to correctly attribute the regional values to such an intermediate level and integrate different data sources that are generally not communicated officially. Some first important results can be observed. The fatalities are observed in relative terms, based on population, and as standardized data. Starting from the data collected at a provincial level and related to the death toll, we could draw up a map of the number of COVID-19 related deaths connected over the population of the area and expressed as deaths per 50,000 inhabitants (Figures 10a and 11).

Mortality rates, SMRs, Population Density, Commuting
The death toll at a provincial level required an in-depth search to correctly attribute the regional values to such an intermediate level and integrate different data sources that are generally not communicated officially. Some first important results can be observed. The fatalities are observed in relative terms, based on population, and as standardized data. Starting from the data collected at a provincial level and related to the death toll, we could draw up a map of the number of COVID-19 related deaths connected over the population of the area and expressed as deaths per 50,000 inhabitants (Figures 10a and 11). The five-class classification shows that a majority of provinces, mainly distributed in Central and Southern Italy (major islands included), present less than 5deaths per 50,000 inhabitants, confirming a relatively low diffusion of the virus in such parts of the country. The highest values can be found in the five provinces in the northeast, east and southeast of Milan, (Bergamo, Brescia, Cremona, and Lodi in the Lombardy region), and Piacenza in Emilia Romagna ( Figure A3, Appendix  B).
The death 'density' seems to be decreasing in such a core area. Central and Southern Italian provinces seem to have very low values, apart from the peaks that can be found in Southern Emilia Romagna (Province of Rimini) and the Marche region (Province of Pesaro in particular).
The standardized mortality ratio (SMR) represents a further relevant analysis, which is fundamental at this stage of the research on the COVID-19 Italian outbreak in order to manage age as confounder. As introduced above, it compares the COVID-19 mortalities with that which is expected, based on the 2019 data. The thematic map produced (Figure 10b) shows quite a neat separation of the value around unity (below unity values) that represents the Italian provinces where, as at 31 st March 2019, the mortality rate is in line with the expected forecasts or even lower. In other provinces, namely the Western Po Valley provinces, including the mountainous ones, and on the Adriatic Coast of the Emilia Romagna and Marche Regions, the Standardized Mortality Ratio is much The five-class classification shows that a majority of provinces, mainly distributed in Central and Southern Italy (major islands included), present less than 5deaths per 50,000 inhabitants, confirming a relatively low diffusion of the virus in such parts of the country. The highest values can be found in the five provinces in the northeast, east and southeast of Milan, (Bergamo, Brescia, Cremona, and Lodi in the Lombardy region), and Piacenza in Emilia Romagna ( Figure A3, Appendix B).
The death 'density' seems to be decreasing in such a core area. Central and Southern Italian provinces seem to have very low values, apart from the peaks that can be found in Southern Emilia Romagna (Province of Rimini) and the Marche region (Province of Pesaro in particular). both for the absolute values (infected and deaths, Cov_57 and Cov_58; Figure 11) and for the relative ones (lethality and SMR, Cov_65, and Cov_66; Figure 12).  This data is also compared to some other indicators that are typical of the 'human geography' of the area, such as population density and commuting. The first element, population density, inhabitants over unit's surface in km 2 , particularly characterizes the Po Valley area, namely in two The standardized mortality ratio (SMR) represents a further relevant analysis, which is fundamental at this stage of the research on the COVID-19 Italian outbreak in order to manage age as confounder. As introduced above, it compares the COVID-19 mortalities with that which is expected, based on the 2019 data. The thematic map produced (Figure 10b) shows quite a neat separation of the value around unity (below unity values) that represents the Italian provinces where, as at 31st March 2019, the mortality rate is in line with the expected forecasts or even lower. In other provinces, namely the Western Po Valley provinces, including the mountainous ones, and on the Adriatic Coast of the Emilia Romagna and Marche Regions, the Standardized Mortality Ratio is much higher than expected. It is interesting to notice that also major urban areas (Turin, Verona and Bologna) present lower values, thereby showing a mortality rate that is apparently less affected by the COVID-19 outbreak. LISA maps on the indicators related to such phenomena confirm these first evaluations, displaying a strong autocorrelation of areas in terms of values and spatial proximity, both for the absolute values (infected and deaths, Cov_57 and Cov_58; Figure 11) and for the relative ones (lethality and SMR, Cov_65, and Cov_66; Figure 12). both for the absolute values (infected and deaths, Cov_57 and Cov_58; Figure 11) and for the relative ones (lethality and SMR, Cov_65, and Cov_66; Figure 12).  This data is also compared to some other indicators that are typical of the 'human geography' of the area, such as population density and commuting. The first element, population density, This data is also compared to some other indicators that are typical of the 'human geography' of the area, such as population density and commuting. The first element, population density, inhabitants over unit's surface in km 2 , particularly characterizes the Po Valley area, namely in two distinct zones, the foothills of the Alps in the north, and the Apennines in the south, although this latter area with much lower values. The Greater Milan area, including Bergamo, one of the provinces most affected by the outbreak, appears as the most uniform and densest, represented by with darker areas ( Figure A4,  Appendix B). This is also confirmed by the LISA analysis (Figure 13a), where provinces around Milan metropolitan city appears positively autocorrelated with reference to density.
those areas that act either as attractor (with values higher than 1, meaning that incoming flows are higher than the outgoing ones) or tributaries (with values lower than 1, meaning that outgoing flows are higher than the incoming ones). The magnitude of the index indicates the 'weight' of the area as the capacity to attract people from neighboring areas, emphasizing metropolitan areas and central cities in a given system. It was interesting to consider the latter index to observe the overall degree of mobility of the area, proposing such a comparison between extra-provincial flows and internal ones ( Figure 14). The Commuting Index (1) illustrates, in the darkest color that represents the highest values, the metropolitan centers, those attracting the higher flows of commuting divided by the population. Major cities and metropolises are highlighted (Figure 14a). The Commuting Index (2) illustrates the overall mobility, expressed as the sum of inbound and outbound flows divided by the flows internal to the area. It can therefore depict both areas of wide mobility and metropolitan roles, as well as a high level of 'openness' in terms of flows (Figure 14b). This analysis was confirmed through LISA (Figure 13b) Cov_83, that displays, through high level of spatial autocorrelation, a generalized self-containment, in terms of inbound and outbound flows, of the Greater Milan area in the extension, as it was introduced in the study area paragraph.  A second element is related to mobility: commuting-although data taken from the 2011 census and available for the whole country, is considered here, different indicators were considered. Two Commuting Indexes were used. The first, a standard one, considers the difference between incoming and outgoing flows over the population of the area. A second index we used here considers the sum of incoming and outgoing flows over the internal flows of the area. The first index tends to emphasize those areas that act either as attractor (with values higher than 1, meaning that incoming flows are higher than the outgoing ones) or tributaries (with values lower than 1, meaning that outgoing flows are higher than the incoming ones). The magnitude of the index indicates the 'weight' of the area as the capacity to attract people from neighboring areas, emphasizing metropolitan areas and central cities in a given system. It was interesting to consider the latter index to observe the overall degree of mobility of the area, proposing such a comparison between extra-provincial flows and internal ones ( Figure 14). The Commuting Index (1) illustrates, in the darkest color that represents the highest values, the metropolitan centers, those attracting the higher flows of commuting divided by the population. Major cities and metropolises are highlighted (Figure 14a). The Commuting Index (2) illustrates the overall mobility, expressed as the sum of inbound and outbound flows divided by the flows internal to the area. It can therefore depict both areas of wide mobility and metropolitan roles, as well as a high level of 'openness' in terms of flows (Figure 14b).

Local Climate Change and Air Quality
A first comment regards the local climate changes affecting the Po Valley megalopolis ( Figure  15) and the significant change of the relative humidity and air quality [123]. These are apparently unconnected phenomena, which in fact, besides being deeply dependent on each other, do not act as a simple sum on the environmental ecosystems and on the community, but in the form of combinations, which are, in turn, correlated with urban geography-land use: efficient land use, sprawl and ecosystem services.
In particular, we have found different regulations in relation to different states, so air quality is not an equal concept for everyone. In this case, to better represent the PM2.5 situation, we used the WHO 2005 guidelines: daily limit 25 µ g/m 3 ; year limit 10 µ g/m 3 ; Italian legislative decree 13 August 2010, no. 255: year limit 25 µ g/m 3 lowered to 20 µ g/m 3 from 1 January 2020-we elaborated the following Figure 16. This analysis was confirmed through LISA (Figure 13b) Cov_83, that displays, through high level of spatial autocorrelation, a generalized self-containment, in terms of inbound and outbound flows, of the Greater Milan area in the extension, as it was introduced in the study area paragraph.

Local Climate Change and Air Quality
A first comment regards the local climate changes affecting the Po Valley megalopolis ( Figure 15) and the significant change of the relative humidity and air quality [123].

Local Climate Change and Air Quality
A first comment regards the local climate changes affecting the Po Valley megalopolis ( Figure  15) and the significant change of the relative humidity and air quality [123]. These are apparently unconnected phenomena, which in fact, besides being deeply dependent on each other, do not act as a simple sum on the environmental ecosystems and on the community, but in the form of combinations, which are, in turn, correlated with urban geography-land use: efficient land use, sprawl and ecosystem services.
In particular, we have found different regulations in relation to different states, so air quality is not an equal concept for everyone. In this case, to better represent the PM2.5 situation, we used the WHO 2005 guidelines: daily limit 25 µ g/m 3 ; year limit 10 µ g/m 3 ; Italian legislative decree 13 August 2010, no. 255: year limit 25 µ g/m 3 lowered to 20 µ g/m 3 from 1 January 2020-we elaborated the following Figure 16. These are apparently unconnected phenomena, which in fact, besides being deeply dependent on each other, do not act as a simple sum on the environmental ecosystems and on the community, but in the form of combinations, which are, in turn, correlated with urban geography-land use: efficient land use, sprawl and ecosystem services.
In particular, we have found different regulations in relation to different states, so air quality is not an equal concept for everyone. In this case, to better represent the PM 2.5 situation, we used the WHO 2005 guidelines: daily limit 25 µg/m 3 ; year limit 10 µg/m 3 ; Italian legislative decree 13 August 2010, Sustainability 2020, 12, 5064 22 of 44 no. 255: year limit 25 µg/m 3 lowered to 20 µg/m 3 from 1 January 2020-we elaborated the following Figure 16.  Figure 16 highlights how the EU represents about 40% of the cities that comply with the WHO 2005 guidelines, the critical learning concentration for the delicate part of the population, particularly the general public and those sensitive individuals who risk having irritation and respiratory problems [124].
Among the urban areas that do not respect the WHO target is the Po Valley Megalopolis, an Italian economic core, with a high population density and a strong commuting aspect (interregional on the EAST side).
Furthermore, according to EEA (2018) [125] the sectors that continue to contribute to the formation of PM2.5 are the commercial, public institutions and households. Similarly, this also occurs with PM10, found in the Po Valley and remained high throughout the Italian lockdown period [126], mainly produced by heating. This failure to reduce PM10 is in fact attributable to the limited progress recorded in Italy on savings on final domestic energy consumption. In fact, since the beginning of the century, our country has shown an improvement performance of 11.6%, well below the European average value of above 30%. According to the European classification for the energy efficiency of buildings (Energy Performance Class, EPC) over 70% of Italian homes falls into a class higher than D [127]. Furthermore, again in Italy, due to the non-renewal of diesel cars with electric and hybrid vehicles (16% registration in 2019), diesel cars have been replaced with petrol ones. Therefore, after several years of improvement, since 2017, CO2 emissions have increased every year [128]. In this ecological approach aimed at evaluating why Northern Italy was marked by COVID-19, the data set (selected from different sources and open data) used for the development of the Lisa Maps, played an important role, which supports and confirms our interdisciplinary assessment. In particular, the LISA maps on indicators relating to these phenomena: Cov_14 (PM2.5), Cov_15 (PM10), Cov_19 (O3), and Cov_72 (O3+PM10), confirm these first evaluations of the air pollution in the Po Valley megalopolis (Figures 17 and 18).  Figure 16 highlights how the EU represents about 40% of the cities that comply with the WHO 2005 guidelines, the critical learning concentration for the delicate part of the population, particularly the general public and those sensitive individuals who risk having irritation and respiratory problems [124].
Among the urban areas that do not respect the WHO target is the Po Valley Megalopolis, an Italian economic core, with a high population density and a strong commuting aspect (interregional on the EAST side).
Furthermore, according to EEA (2018) [125] the sectors that continue to contribute to the formation of PM 2.5 are the commercial, public institutions and households. Similarly, this also occurs with PM 10 , found in the Po Valley and remained high throughout the Italian lockdown period [126], mainly produced by heating. This failure to reduce PM 10 is in fact attributable to the limited progress recorded in Italy on savings on final domestic energy consumption. In fact, since the beginning of the century, our country has shown an improvement performance of 11.6%, well below the European average value of above 30%. According to the European classification for the energy efficiency of buildings (Energy Performance Class, EPC) over 70% of Italian homes falls into a class higher than D [127]. Furthermore, again in Italy, due to the non-renewal of diesel cars with electric and hybrid vehicles (16% registration in 2019), diesel cars have been replaced with petrol ones. Therefore, after several years of improvement, since 2017, CO 2 emissions have increased every year [128]. In this ecological approach aimed at evaluating why Northern Italy was marked by COVID-19, the data set (selected from different sources and open data) used for the development of the Lisa Maps, played an important role, which supports and confirms our interdisciplinary assessment. In particular, the LISA maps on indicators relating to these phenomena: Cov_14 (PM 2.5 ), Cov_15 (PM 10 ), Cov_19 (O 3 ), and Cov_72 (O 3 +PM 10 ), confirm these first evaluations of the air pollution in the Po Valley megalopolis (Figures 17 and 18).  In this complex picture of air pollutants in the air, from the diversity of the assessment of international targets and the relative multiple production sectors, in Italy and in particular in the Po valley, further obstacles are added in terms of physical geography (river valley) and climatic conditions, expressed according to the climatic well-being index (sunshine, perceived temperature,  In this complex picture of air pollutants in the air, from the diversity of the assessment of international targets and the relative multiple production sectors, in Italy and in particular in the Po valley, further obstacles are added in terms of physical geography (river valley) and climatic conditions, expressed according to the climatic well-being index (sunshine, perceived temperature, In this complex picture of air pollutants in the air, from the diversity of the assessment of international targets and the relative multiple production sectors, in Italy and in particular in the Po valley, further obstacles are added in terms of physical geography (river valley) and climatic conditions, expressed according to the climatic well-being index (sunshine, perceived temperature, heat waves, extreme events, breeze, relative humidity, gusts of wind, rain, fog) which significantly influences the air quality [129].
In this sense, the LISA maps on indicators relating to these phenomena: Cov_37 (climate wellbeing), Cov_39 (wind gusts), Cov_41 (fog), Cov_55 (wind speed) confirm these evaluations of the climate conditions (Figures 19 and 20). heat waves, extreme events, breeze, relative humidity, gusts of wind, rain, fog) which significantly influences the air quality [129].
In this sense, the LISA maps on indicators relating to these phenomena: Cov_37 (climate wellbeing), Cov_39 (wind gusts), Cov_41 (fog), Cov_55 (wind speed) confirm these evaluations of the climate conditions (Figures 19 and 20).  Recalling that climatic variations affect air quality and that air pollution induces climatic variations, we can dispute that the LISA Maps have constituted a valid support for the evaluation of autocorrelation, which will be the subject of further research developments. In this sense, the LISA maps on indicators relating to these phenomena: Cov_37 (climate wellbeing), Cov_39 (wind gusts), Cov_41 (fog), Cov_55 (wind speed) confirm these evaluations of the climate conditions (Figures 19 and 20).  Recalling that climatic variations affect air quality and that air pollution induces climatic variations, we can dispute that the LISA Maps have constituted a valid support for the evaluation of autocorrelation, which will be the subject of further research developments. Recalling that climatic variations affect air quality and that air pollution induces climatic variations, we can dispute that the LISA Maps have constituted a valid support for the evaluation of autocorrelation, which will be the subject of further research developments.

Land Take and COVID-19
As previously explained, the physical conformation of the northern part of the Italian peninsula is influenced by the Alpine Chain. It represents a sort of barrier against the winds, which influence air circulation and distribution. This hypothesis has been confirmed by the LISA analysis Cov_39 (Figure 19), which highlights a strong spatial autocorrelation of a low level of data concerning the annual days with gusts of wind stronger than 25 knots.
Considering also that most of the productive activities are concentrated in this area, it should be important to have a lot of non-urbanized areas which enable CO 2 storage. Cov_82 (Figure 21) shows high values of spatial autocorrelation of CO 2 related to non-urbanized areas. In particular, the areas with the greatest number of positive cases and deaths do not have enough spaces with natural functions able to store the produced CO 2 . decrease of PM2.5 PM10 and O3, is highlighted in Figures 17 and 18 where the highest values of spatial autocorrelation occur in COVID-19 red zones.
The Italian Institute for Environmental Protection and Research (ISPRA) produces an annual report on Land Take, Territorial Dynamics, and Ecosystem Services which highlights all the critical issues [130]. Unfortunately, according to this report, the northern part of Italy is the area where the phenomenon of sprawl is more concentrated.
This trend is also confirmed by a spatial autocorrelation analysis. Figures 21 and 22    Considering these structural features of the territory, unsealed soils can play a fundamental role in CO 2 storage and, at the same time, green areas can be very useful in improving air quality, reducing PM 2.5 PM 10 and O 3 . The lack of forest and green areas able to raise air quality, thereby determining a decrease of PM 2.5 PM 10 and O 3 , is highlighted in Figures 17 and 18 where the highest values of spatial autocorrelation occur in COVID-19 red zones.
The Italian Institute for Environmental Protection and Research (ISPRA) produces an annual report on Land Take, Territorial Dynamics, and Ecosystem Services which highlights all the critical issues [130]. Unfortunately, according to this report, the northern part of Italy is the area where the phenomenon of sprawl is more concentrated.
This trend is also confirmed by a spatial autocorrelation analysis. Martellozzo et al. [131] analysed on a national level, how an approach based on laissez-faire generated an uncontrolled growth in most of the country. In particular, this trend is more Martellozzo et al. [131] analysed on a national level, how an approach based on laissez-faire generated an uncontrolled growth in most of the country. In particular, this trend is more concentrated in the northern part of the country.
In this study, two simulations of change in land use for 2030 have been developed, the first adopts a more sustainable approach and the second considers a "business-as-usual" development. Despite big differences in urban sprawl between the two scenarios, in the northern part of Italy land take reaches important values, also when analysing the sustainable scenario. This trend has been described in the Lombardy region, with a detailed analysis on Milan and Brescia by Pileri [132].
This phenomenon is mainly due to the implementation of investments in total absence of planning or with very old tools. Even in the case of recent planning, the new tools have been applied with old laws thereby generating a paradoxical gap between the plan contents and territorial reality. For this reason, some authors adopt the term "Vintage Urban Planning" [133] or "Ghost Planning" [134]. Many of the interventions are carried out without an overall strategic framework able to relate investments to territorial features [135]. This is possible because, in a lot of cases, the municipalities possess old plans, with demographic projections calculated in periods of great population growth, producing an impressive availability of development areas with a high chance of never completely using them to this aim. In the more competitive areas of the country, e.g., Northern Italy, these lands are owned by big companies that decide to transform them, not just for selling, but also for the advantage that could be attained in inflating their accounting budget, assigning real estate assessments higher than market values [136,137]. The continuous urbanization of northern Italy is also due to the strong demand that derives from migratory flows from Southern Italy and from abroad [138].

Diffusion Processes, Local and Global Effects
In this framework, we hypothesized some characters of the diffusion process that occurred in Italy in the early stages of the virus outbreak among the different areas. Considering that a quite limited set of data is available to date and the process appears to remain ongoing, a conclusion cannot therefore be drawn.
Besides the debate concerning the existing relation between the two virus hotbeds, Codogno and Vo, we could dispute that, being minor centers in their areas, we faced a hierarchical, bottom-up diffusion process towards medium and major centers in the Lombardy, Piedmont, Veneto, and Emilia Romagna regions. Codogno in particular is located in the center of a triangle of three medium-sized cities, such as Piacenza, Cremona and Lodi, at the top of the particulate emissions index in Italy for several years, as in Table 1, and very well connected to the Milan area and the industrial region on its eastern outskirts. The connected motorway and state road system is the backbone of commuting in an area characterized by a high level of accessibility, although often saturated in terms of heavy trucks and car congestion. The Brescia Provinces, as well as a by-pass to enter the city of Milan from the south, avoiding congestion of the A1 in accessing the orbital route around the city (Figure 23).    The bottom-up hierarchical diffusion process possibly reached such medium and bigger centers, spreading then, as a waterfall diffusion process, top-down, towards minor centers, and again activating local contagious processes. Such higher order centers can be considered to be those cities that host specialized services and facilities, such as hospitals, health care structures and retirement homes.
However, hierarchical, top-down diffusion presumably occurred mainly from medium-density and sized cities: Bergamo, Brescia, Cremona, Parma Piacenza, Rimini, and Pesaro, just to cite among the most involved ones which were severely hit, but proportionally less than major centers such as Milan, Turin, and Bologna. These latter, major cities could be primarily involved, in their numbers, thanks to the presence of hospitals and retirement homes as hotbeds for further infection, rather than as dense environments which favor stronger diffusion processes.
At this stage of data availability and knowledge, we cannot make particular considerations, regarding the pre-and post-lock down periods, about the spread of COVID-19 by hierarchical diffusion processes related to medium and long connections, such as air and rail transport, including high speed trains particularly in the north-south corridors, hence it is worthy of attention in further research.
The above-mentioned considerations regarding the local spatial diffusion of COVID-19, as well as the issues related to air quality and land take, can be observed not only in the analysis of the results from the different methods adopted, but also in a synthetic table (Table 1). Here, we ranked the Italian provinces where the SMR as at 31 March is higher than 1, which indicates a higher increase of mortality than expected. It is a set of 29 provinces, where, by means of example in the Po Valley area, a city like Bergamo, with a value of 12.356, presents an increase 12 times more than expected. We can notice that the most affected areas are localized in the Po Valley and are particularly characterized by a high number of days in which the limits of air particulate emissions and soil consumption have been exceeded. One can also notice a similarity in the density classes, which mainly include between 300 and 400 people per square kilometer. As other research pointed out [25], they host their provincial capital spanning from 100,000 to 200,000 inhabitants ( Figure 23, Figure A2). Furthermore, a similar level of high, interprovincial mobility is reached, seeing that most of the cities are tributaries of other major metropolitan areas, or traffic collectors because of the presence of major industrial activities, i.e., Bergamo and Brescia. A negative commuting index (1) generally suggests major outgoing flows towards other centers; a high commuting index (2) suggests a high level of overall external mobility and therefore openness to other provinces and areas. Such elements appear as interesting in helping to understand both the potential local diffusion of COVID-19, as well as the patterns and consequences on the air quality and land take.

Air Issues and Policies
In this context of diffusion from top to bottom, the geographical, climatic, and air quality conditions played an important role and contributed to increasing the effects.
In particular, local climate changes such as temperature and humidity, poor air quality and the persistent absence of wind, make the Po Valley a one-of-a-kind area, both at a national and international level. Furthermore, frequent and persistent thermal inversion phenomena in the winter months, especially in periods of high atmospheric pressure, trap the cold air to the ground, together with the pollutants.
In the Po valley, urban industrialization and intensive agriculture intertwine and more than 50% of the national GDP is produced, and almost 50% of the national energy is consumed [139]. The transition to compatible solutions related to well-being does not exactly seem close, despite the several regional and national plans to monitor and improve air quality. In fact, with the COVID-19 outbreak, PM 10 emissions in the Po Valley were high and sometimes exceeded the limits and are to be attributed to a combined climate action (wind, winter thermal inversion) and activities, as shown ( Figure 24). of the national GDP is produced, and almost 50% of the national energy is consumed [139]. The transition to compatible solutions related to well-being does not exactly seem close, despite the several regional and national plans to monitor and improve air quality. In fact, with the COVID-19 outbreak, PM10 emissions in the Po Valley were high and sometimes exceeded the limits and are to be attributed to a combined climate action (wind, winter thermal inversion) and activities, as shown ( Figure 24). In this context, which is certainly not simple for human, environmental, and anthropic geography, the ecological approach has allowed us to obtain the first results and the first policy proposals supported by the LISA analysis ( Figure 25). In this context, which is certainly not simple for human, environmental, and anthropic geography, the ecological approach has allowed us to obtain the first results and the first policy proposals supported by the LISA analysis ( Figure 25). The results produced are in line with the most recent scientific literature of the Po Valley and COVID-19 [140][141][142] and its guidelines: thinking globally and acting locally.

New Approaches to Planning and Policies
Generally, plans are the results of long negotiations, which produce a long list of strict norms. Consequently, the main aim of this approach to planning is to enforce such rules. In several cases, The results produced are in line with the most recent scientific literature of the Po Valley and COVID-19 [140][141][142] and its guidelines: thinking globally and acting locally.

New Approaches to Planning and Policies
Generally, plans are the results of long negotiations, which produce a long list of strict norms. Consequently, the main aim of this approach to planning is to enforce such rules. In several cases, plans are old, very far from the current reality, or based on old laws, which do not allow for the production of tools to solve current problems. The result is that the main planning goals are very far from providing a serious response to the transformation demands that arise daily.
This approach, based on vintage planning [133] or ghost planning [134], leads to a situation synthetized on the left of Figure 26. The results are different: a city with an uncontrolled real estate market, which produces urban sprawl; several portions of the city inhabited only by tourists due to the effects of Airbnb; cities dominated by cars with serious pollution problems. This scenario leads to a consumption of resources greater than the capacity of the planet. In 2019, Overshoot Day occurred in late July. Consequently, for the following five months, humanity used resources that the planet could not provide. An alternative is an approach based on simulations in assessing transformation impacts, allowing planners to take into account several land use scenarios, choosing the more suitable solutions for the transformations. This approach to planning also considers possible losses of the ecosystem services in simulations [61,[143][144][145]. can create a more sustainable scenario able to flatten the curve under the earth's carrying capacity [160][161][162].

Conclusions
In this paper, we tried to find some possible answers to questions related to the severe COVID-19 outbreak, which saw Italy as the first country in Europe after the Southeastern Asian ones. In particular, we tried to highlight some elements related to the causes of the development of virus hotbeds in Northern Italy, the Milan and Po Valley areas in particular, and their diffusion in other parts of the country. In this sense, we analysed COVID-19 related data, in terms of cases and deaths until 31 March 2020, in order to consider the situation at a suitable time after the hard lockdown policies set up by the Italian government-from 10 March 2020-and which, in posing severe limitations to movement and production, should have stopped the spread of the disease. Therefore, our observations can be considered in a "frozen" situation of the diffusion process. The analysis was involved observing the spatial diffusion and distribution of the phenomenon and observing the COVID-19 outbreak with reference to some major groups of variables to elucidate the possible causes of such a violent phenomenon: land use, air quality, climate and weather, population, health and life expectancy. These have been analysed considering their spatial autocorrelation by means of LISA.
As recent studies anticipated, multiple origins of the virus in Italy, through Germany and Several authors adopted the term performance-based planning [146][147][148][149][150][151][152][153] to synthetize this approach. Consequently, this "umbrella" can contain all simulation models and tools. Due to data availability, all models based on Cellular automata or Multiagent Systems, Space Syntax, Geodesign [154][155][156], etc. can take into account a lot of components in detailed simulations. Therefore, goals for the protection of natural areas will be pursued more easily. Furthermore, adopting urban policies based on urban regeneration, sustainable mobility and the creation of green infrastructures [157][158][159] can create a more sustainable scenario able to flatten the curve under the earth's carrying capacity [160][161][162].

Conclusions
In this paper, we tried to find some possible answers to questions related to the severe COVID-19 outbreak, which saw Italy as the first country in Europe after the Southeastern Asian ones. In particular, we tried to highlight some elements related to the causes of the development of virus hotbeds in Northern Italy, the Milan and Po Valley areas in particular, and their diffusion in other parts of the country. In this sense, we analysed COVID-19 related data, in terms of cases and deaths until 31 March 2020, in order to consider the situation at a suitable time after the hard lockdown policies set up by the Italian government-from 10 March 2020-and which, in posing severe limitations to movement and production, should have stopped the spread of the disease. Therefore, our observations can be considered in a "frozen" situation of the diffusion process. The analysis was involved observing the spatial diffusion and distribution of the phenomenon and observing the COVID-19 outbreak with reference to some major groups of variables to elucidate the possible causes of such a violent phenomenon: land use, air quality, climate and weather, population, health and life expectancy. These have been analysed considering their spatial autocorrelation by means of LISA.
As recent studies anticipated, multiple origins of the virus in Italy, through Germany and directly in Italy from China (particularly in the Lazio region), could be traced, but it is in the Po Valley area that the most severe outbreak took place, and then spread to other parts of Italy.
We observed similarities in the conditions of the Wuhan area in the Hubei Province, with those in the Po Valley metropolis. These seemed particularly related to the geographical, climatic-(presence of water bodies, flat lands, limited air circulation, similar climate zones) and socio economic ones, industrial production, transport infrastructures and mobility, population distribution and density, as well as similarities in terms of presence, concentration and persistence of pollutants in the atmosphere.
We hypothesized the existence of relationships between pollution and both the spread of the virus in generating outbreaks and the lethality of the disease. In particular, we took into consideration land take and air pollution, this latter referred to particulates (PM 2.5 and PM 10 ) and Nitrogen based components, such as NO x and NH 3 deriving from human activities (industry, traffic, home heating); agriculture and cattle breeding. The idea was that the presence of air-related pollutants can generate pressure on the health conditions of the populations at risk and offer preconditions for both the development of respiratory related diseases and complications, including some that are life-threatening, which may explain the high death rate observed in the area under study. Moreover, particular atmospheric conditions in the early weeks of 2020 could have worsened the environmental situation in the Wuhan and Po Valley areas, namely scarce precipitation and a warmer winter than the previous ones.
The analysis on COVID-19 related mortality showed quite a neat divider between Northern Italy on one side and Central and Southern Italy on the other, along the Apennine mountain chain, with values higher (much higher than expected) in the north, and values in line the with expected mortality in other Italian regions, particularly in the south. Starting from this, we observed the spatial distribution of COVID-19 related deaths through SMR with climatic, pollutant, and human-pressure indicators, finding evidences of similarities in the spatial patterns drawn particularly with PM 2.5 and other particulates, Nitrogen-related pollutants, land take, population density and commuting, mainly in the areas of the Western Po Valley, the Milan and neighboring areas towards the Provinces of Bergamo and Brescia, and in the south and southeastern part, towards the Emilia Romagna Region (Cremona, Piacenza, Parma, etc.).
In geographical terms of spatial diffusion processes, at a global level we would have expected, hierarchically, a diffusion following the major air connections, and therefore outbreaks in areas of the world that are more connected to China. These can be the cases of the United States (West Coast in particular) and of Central and Northern Europe, or, in the Italian case, an outbreak following a diffusion process could have started in the Lazio region: a direct flight existed between Wuhan Tianhe and Rome Fiumicino airport [163,164]. As geographers point out [80], an outbreak takes place where a susceptible population can be detected, and this evidently took place starting from the two hotbeds of Vo and Codogno, two municipalities in the Po Valley area. Codogno, more specifically, in the southern outskirts of the Metropolitan area of Milan is located at the center of three of the cities in Italy-Lodi, Cremona and Piacenza-that have registered (for consecutive years) the most days per year on which the limits of particulate presence in the atmosphere were exceeded.
Local, contagion, and hierarchical diffusion processes could be hypothesized starting from this area, with more severe cases in the Lombardy, Piedmont, Emilia Romagna, and Marche Regions in particular, confirmed by the COVID-19 related Standardized Mortality Ratio computed over all the Italian provinces. The main branches of the diffusion appear related to the main gravitation axis of the Milan metropolitan areas: westwards towards Piedmont and bordering provinces, eastwards towards the industrial cities of Bergamo and Brescia, south-eastwards towards the Lombardy and Emilia Romagna regions, and through the Marche region, along the 'Via Emilia'. Generally, the orientation seems to be compatible with the major communication routes, namely the A4 and A1 motorways as well as the state road system. However, when observing the COVID-19 outbreak in terms of mortality and case diffusion, a major incidence can be noticed particularly in the centers of average density and dimension, more than in the major, denser cities. Denser and bigger regional capital cities, namely Milan, Turin, Verona, and Bologna appear to be relatively less affected than their neighboring, surrounding provinces. This element is also apparently in contrast with a model of a hierarchical, top-down diffusion from major centers towards minor ones. The deindustrialization of these capital cities, and their orientation towards an economy of services, could be related to that, as well as the different set up of social and human relations in the medium-sized cities if compared with the bigger ones. However, this was not analyzed in this research and needs further elaboration and analysis.
What has emerged therefore is that a profound reflection on the monitoring of air emissions is required, in particular of the PM 10 , which did not substantially decrease during the lockdown. In fact, through monitoring it is possible to verify the effectiveness of the measures implemented at a local level to reduce air pollution. This precaution should be included in the agreement of the Po Valley for the improvement of air quality, signed in Bologna during the G7 Environment Ministers' Meeting of 9 June 2017, by the Minister of the Environment and the presidents of Lombardy, Piedmont, Veneto and Emilia-Romagna.
As regards to suggestions in terms of policies, we can strengthen the need to act in favor of policies aimed at the reduction of pollutants in the atmosphere, by means of speeding up the already existing plans and policies, aimed at all the sources of atmospheric pollution: industries, home heating and traffic. That, however, does not appear sufficient if (as noticed through the satellite images and related pollution maps-ESA before and after the lock down), as the reduction of pollutants in the Po Valley area was not as intense as expected, demonstrating an inertia of pollutant in the air due to the climate and geographical conditions of the area. Investment in clean transport and building should therefore be reinforced, beginning with rapidly applicable measures, e.g., road washing, pollution eating paints, façades, and plants, etc.
The results so far obtained, as well as the conclusions hypothesized so far, need to be better studied and understood with further data, which at present is still not validated and available. With reference to mortality, our analysis on SMR with data at 31 March 2020 were recently compared with the situation at 30 April 2020, which analysis is ongoing in another research. That showed how the absolute pattern did not change among the different provinces, and therefore validated the first results obtained. In particular, more observational studies based on individual complete data related to the COVID-19 cases and deaths, registered also at a finer geographical scale than the regional and provincial ones, would allow overcoming the intrinsic limit of the use of provincial units for data collection, in terms of heterogeneity of their shape and population. We also need to better analyze the overall amount of data and indicators collected so far, also in terms of the evolution of the virus outbreak and its hoped-for reduction, in light of the other results obtained in medical and epidemiological research by other research groups. In this paper, the analysis was carried on considering the Italian situation, although in more recent times other countries experimented the COVID-19 outbreak on their territories, also with more severe consequences than in the Italian case. Representing Italy as a sort of pioneer in Western countries, the analyses appeared quite relevant as no previous case studies could be observed and therefore compared.
Future developments of the present research include an analysis of the full dataset realized with further spatial analytical techniques for a COVID-19 risk map for Italy.