Background Healthcare factors have strongly influenced the propagation of COVID-19. The present study aims to examine whether excess mortality during the first phase of the COVID-19 outbreak in Italy was associated with health, healthcare, demographic, and socioeconomic indicators measured at a provincial level.

Methods The present ecological study concerns the raw number of deaths from Jan. 1 to Apr. 30, 2020 and the mean number of deaths in the same months of 2015 to 2019, per province. Information on socioeconomic factors and healthcare settings were extracted from the most recently updated databases on the ISTAT website. Two multilevel, multivariate models were constructed to test whether excess mortality was associated with the indicators across 107 provinces in Italy.

Results On linear multilevel, multivariate analysis, AIDS mortality rate (p-value <0.05) correlates positively with excess mortality, while a higher density of General Practitioners (number of GPs per 1,000 population) is associated with lower excess mortality (p-value <0.05). After controlling for the diffusion of COVID-19 in each province, the significance of GP density increases (p-value <0.001) and the rate of hospitalization in long-term care wards is positively associated (p-value <0.05) with excess mortality.

Conclusion Some health and healthcare variables are strongly associated with excess mortality caused by COVID-19 in Italy and should be considered to implement mitigation policies and increase healthcare resilience.