Before revision:  We start with an ODE version of a SEIHARD model that we optimize to obtain relevant parameters for these regions . Subsequently, spatially distributed versions of the corresponding models in the form of convection-reaction-diffusion equations are considered. We show that these models bear significant potential towards capturing the spatial expansion of the pandemic and veins of future work are considered to that effect .
After revision:  We start with an ODE version of a SEIHARD model that we optimize to obtain relevant parameters for these regions . Subsequently, spatially distributed versions of the corresponding models in the form of reaction-diffusion equations is developed. We consider that, after an initial localized seeding of the infection, its spread is governed by the diffusion (and 0D model "reactions") of the asymptomatic and symptomatically infected populations, which decrease with the imposed restrictive measures. We inserted the maps of the two regions, and we imported population-density data into COMSOL, which was subsequently used to solve numerically the model PDEs. Upon discussing how to adapt the 0D model to this spatial setting, we show that these models bear significant potential towards capturing the spatial expansion of the pandemic and veins of future work are considered to that effect .