Modeling provincial Covid-19 epidemic data in Italy using an adjusted time-dependent SIRD model