Skip to main content

Changes in the spatial distribution of COVID-19 incidence in Italy using GIS-based maps

Dear editor

Massive spreading of the pandemic Coronavirus Disease 2019 (COVID-19) in different continents [1, 2], have been observed [3]. Analyses mostly focused on the number of cases per country and administrative levels, multiple times without considering the relevance of the incidence rates. These help to see the concentration of disease among the population in terms of cases per 100,000 inhabitants. Even more, using geographical information systems (GIS)-based maps, stakeholders may rapidly analyze changes in the epidemiological situation [4,5,6,7]. Although the epidemic of COVID-19 caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) started in Italy on January 31, 2020, no reports on the use of GIS-based maps have been published to analyze the distinct differences in incidence rates across its regions and provinces during the last months. For these reasons, we have developed epidemiological maps of incidence rates using official populations, by regions (1st administrative level of the country) and provinces (2nd administrative level), for COVID-19 in Italy using GIS.

Surveillance cases data of the cumulative number at March 15, April 18, and June 8, 2020, officially reported by the Italian health authorities were used to estimate the cumulated incidence rates on those dates using reference population data on SARS-CoV-2 confirmed infections (cases/100,000 pop) and to develop the maps by regions and provinces, using the GIS software Kosmo® 3.1, as performed in previous related studies [6, 7]. Starting on March 8, 2020, the region of Lombardy, together with 14 additional northern and central provinces, in Piedmont, Emilia-Romagna, Veneto, and Marche, they were put under lockdown. On March 10, 2020, the government extended the lockdown measures to the whole country.

Up to March 15, 2020, after 44 days of epidemics, 24,053 cases of COVID-19 were reported in the country, for a cumulated rate of 39.6 cases/100,000 population, reaching 174,103 cases during April 18, 2020, for a rate of 286.78, and 232,855 cases during June 8, 2020, for a rate of 383.56. All the regions of the country have been affected, with rates ranging from 59.42 (Calabria) to 947.75 cases/100,000 population (Aosta Valley/Vallée d’Aoste) (June 8, 2020) (Fig. 1). Higher diversity is found in provinces, where incidence rates ranged from 28.23 (Sud Sardegna, Sardinia) to 1811.37 (Cremona, Lombardy) (June 8, 2020) (Table 1). At Lombardy are located five of the top ten provinces with higher incidence rates (Table 1, Fig. 1), with considerable increases and changes from March 15, 2020, to June 8, 2020, in approximately 3 months (Fig. 1, Table 1). Cremona (Lombardy), Piacenza (Emilia-Romagna), and Lodi (Lombardy) have become in the geographic core of the cumulated incidence rate of COVID-19 in the north of the country and Italy (Fig. 1).

Fig. 1

COVID-19 situation in Italy, on March 15, April 18, and June 8, 2020, by regions and provinces

Table 1 Top ten provinces by incidence rate (cases/100,000 inhabitants), of COVID-19, Italy, on March 15, April 18, and June 8, 2020

From the GIS-based maps, it is clear that spreading in the country is occurring from north regions and provinces such as Lombardy. On March 15, 2020, most of the southern regions were not affected (Fig. 1), but approximately a month later, all of them reported COVID-19 cases (Fig. 1), including the insular regions of Sicily and Sardinia. While the change between April 18 and June 8, has been 33%, there is still a concern in the country, mainly because, in this time, the number of deaths has reached over 34,000 deaths (14.6%).

Italy reached the top of countries with the highest number of reported COVID-19 cases, now is the ninth country in cumulated cases. It is the fourth in the European region, after Russia, the United Kingdom, and Spain. Italy represented one of the most significant sources of imported cases for other continents, as is the case of Latin America, that received their first cases from Milan, Lombardy [8,9,10].

Patient 1 of Italy (it was not possible to find patient 0) was discovered on February 20, 2020, when a 38-year-old man from the city of Codogno had shown up at the hospital. Since that date, two large clusters of outbreaks have spread first in Northern Italy, later all over the country (Fig. 1) [11]. Cases are multiplying, and the national healthcare system is collapsing [12,13,14]. Many regions are increasing intensive care beds, revolutionizing entire hospital wards. In Italy, the coordination of the swabs is managed regionally. Once the epidemic began, for example, the Veneto region started immediately with active surveillance, i.e., on asymptomatic, and this contained the spread of the virus compared to other Northern Regions.

The health system is indeed regionalized, and dispositions of the Ministry of Health are translated into multiple regional decrees and regulations, often at different timings [15]. Many regions, for example, adopted evolving criteria for testing and diagnosis, according to dispositions from the central Government, but also to test capacity, which was heavily reliant on the availability of reagents. Samples were collected either in healthcare facilities, in provisional collection points, or even in people’s houses, depending on the region and the phase of the pandemic. The epitome of these differences is the two most heavily affected regions, Lombardy and Veneto. Lombardy hospitalized, even cases with relatively modest symptoms, causing numerous nosocomial outbreaks (9% of infections were among health professionals until March). Veneto, instead, deployed widespread testing since the beginning, maintaining disease management as much as possible at the primary health care level [15].

The healthcare workers are facing COVID-19 pulling 12 h shifts in critical situations with minimal to non-existent personal protective equipment (PPE) [12,13,14]. Lacking PPE led both many healthcare workers to become COVID-19 positive (7145), and to the death of several doctors (51; about 9% of the total cases; March 27, 2020) [11].

As observed in the GIS-based maps, the COVID-19 spreading in the country has been significant and moving from north to south across the time, with provinces reaching more than 1000 cases per 100,000 inhabitants (Fig. 1) [12,13,14]. Differences in the incidence by regions would be related to different social and economic factors. Such as people who travel abroad, for whom there is a sharp difference between northern regions (about 26% of travelers) and central and southern regions (about 19%). Or net income at the household level, which is ranging from 35,000€ in the North-east to 26,000€ in the South [16]. Additionally, as has been recently suggested, climatic conditions could also influence the transmission of SARS-CoV-2 [17].

Testing capacity increased over time. While some degree of ascertainment bias is inevitable, the number of swabs performed nationally stabilized around 70,000 per day in mid-April when the peak in the number of active cases was registered (Fig. 1). From that point onwards, daily cases only decreased, consistently with the impact of the lockdown imposed in early March. Besides, and most importantly, while it is true that the epidemic might have gone undetected for some time before Case 1 was discovered in Codogno, the growth in ICU beds demand for subjects with respiratory failure that ensued in the following weeks is most likely explained by a substantial increase of cases.

Considering the limitations of diagnostics and the asymptomatic cases, these figures would be many times more. Further characterization studies should include multiple GIS-based maps with other variables at the regions and provinces levels such as deaths, hospitalizations, and ICU rates per population to understand better the critical situation of the country and its administrative levels.

Availability of data and materials

If required.


  1. 1.

    Rodriguez-Morales AJ, MacGregor K, Kanagarajah S, Patel D, Schlagenhauf P. Going global—travel and the 2019 novel coronavirus. Travel Med Infect Dis. 2020;33:101578.

    Article  Google Scholar 

  2. 2.

    Rodriguez-Morales AJ, Bonilla-Aldana DK, Balbin-Ramon GJ, Paniz-Mondolfi A, Rabaan A, Sah R, Pagliano P, Esposito S. History is repeating itself, a probable zoonotic spillover as a cause of an epidemic: the case of 2019 novel Coronavirus. Infez Med. 2020;28(1):3–5.

    CAS  PubMed  Google Scholar 

  3. 3.

    Dhama K, Khan S, Tiwari R, Sircar S, Bhat S, Malik YS, Singh KP, Chaicumpa W, Bonilla-Aldana DK, Rodriguez-Morales AJ. Coronavirus disease 2019-COVID-19. Clin Microbiol Rev. 2020;33(4):e00028-20.

    Article  Google Scholar 

  4. 4.

    Bastola A, Sah R, Rodriguez-Morales AJ, Lal BK, Jha R, Ojha HC, Shrestha B, Chu DKW, Poon LLM, Costello A, et al. The first 2019 novel coronavirus case in Nepal. Lancet Infect Dis. 2020;20(3):279–80.

    CAS  Article  Google Scholar 

  5. 5.

    Dong E, Du H, Gardner L. An interactive web-based dashboard to track COVID-19 in real time. Lancet Infect Dis. 2020;20:533–4.

    CAS  Article  Google Scholar 

  6. 6.

    Mousavi SH, Zahid SU, Wardak K, Azimi KA, Reza Hosseini SM, Wafaee M, Dhama K, Sah R, Rabaan AA, Arteaga-livias K, et al. Mapping the changes on incidence, case fatality rates and recovery proportion of COVID-19 in Afghanistan using geographical information systems. Arch Med Res. 2020.

    Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Arab-Mazar Z, Sah R, Rabaan AA, Dhama K, Rodriguez-Morales AJ. Mapping the incidence of the COVID-19 hotspot in Iran—implications for travellers. Travel Med Infect Dis. 2020;34:101630.

    Article  Google Scholar 

  8. 8.

    Rodriguez-Morales AJ, Gallego V, Escalera-Antezana JP, Mendez CA, Zambrano LI, Franco-Paredes C, Suarez JA, Rodriguez-Enciso HD, Balbin-Ramon GJ, Savio-Larriera E, et al. COVID-19 in Latin America: the implications of the first confirmed case in Brazil. Travel Med Infect Dis. 2020;35:101613.

    Article  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Escalera-Antezana JP, Lizon-Ferrufino NF, Maldonado-Alanoca A, Alarcon-De-la-Vega G, Alvarado-Arnez LE, Balderrama-Saavedra MA, Bonilla-Aldana DK, Rodriguez-Morales AJ. Lancovid: clinical features of the first cases and a cluster of coronavirus disease 2019 (COVID-19) in Bolivia imported from Italy and Spain. Travel Med Infect Dis. 2020;35:101653.

    Article  Google Scholar 

  10. 10.

    Millan-Onate J, Millan W, Mendoza LA, Sanchez CG, Fernandez-Suarez H, Bonilla-Aldana DK, Rodriguez-Morales AJ. Successful recovery of COVID-19 pneumonia in a patient from Colombia after receiving chloroquine and clarithromycin. Ann Clin Microbiol Antimicrob. 2020;19(1):16.

    CAS  Article  Google Scholar 

  11. 11.

    De Giorgio A. COVID-19 is not just a flu Learn from Italy and act now. Travel Med Infect Dis. 2020;35:101655.

    Article  Google Scholar 

  12. 12.

    Livingston E, Bucher K. Coronavirus disease 2019 (COVID-19) in Italy. JAMA. 2020;323:1335.

    Article  Google Scholar 

  13. 13.

    Onder G, Rezza G, Brusaferro S. Case-fatality rate and characteristics of patients dying in relation to COVID-19 in Italy. JAMA. 2020;323:1775–6.

    CAS  Google Scholar 

  14. 14.

    Remuzzi A, Remuzzi G. COVID-19 and Italy: what next? Lancet. 2020;395(10231):1225–8.

    CAS  Article  Google Scholar 

  15. 15.

    Boccia S, Ricciardi W, Ioannidis JPA. What other countries can learn from Italy during the COVID-19 pandemic. JAMA Intern Med. 2020;180:927–8.

    Article  Google Scholar 

  16. 16.

    ISTAT data warehouse. 2020. Accessed 19 June 2019.

  17. 17.

    Mendez-Arriaga F. The temperature and regional climate effects on communitarian COVID-19 contagion in Mexico throughout phase 1. Sci Total Environ. 2020;735:139560.

    CAS  Article  Google Scholar 

Download references





Author information




AJRM, conceptualization; data curation; formal analysis; methodology; software; writing-original draft; writing-review and editing. CAM, data curation; formal analysis; methodology; writing-review and editing. RS, AAR, KD, CC, KAL, TS, AO, DB, AH, YK, ZF, NR, MAUR, TT, writing-review and editing. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Alfonso J. Rodriguez-Morales.

Ethics declarations

Ethics approval and consent to participate

Not required.

Consent for publication

Not required.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Martellucci, C.A., Sah, R., Rabaan, A.A. et al. Changes in the spatial distribution of COVID-19 incidence in Italy using GIS-based maps. Ann Clin Microbiol Antimicrob 19, 30 (2020).

Download citation


  • SARS-CoV-2
  • COVID-19
  • Geographical information systems
  • Pandemic
  • Italy