Author contributions
M.d.B. with support from F.D.R. and P.D.L., D.L., G.R. and F.L.I., who contributed equally to this work, designed the research; F.D.R., A.D.M, P.D.L. and D.S. carried out the model parameterization; F.L.I., C.C. and A.G. analysed the data and estimated the inter-regional flux matrix; F.D.L., M.C. and R.C. wrote and checked the numerical code used for all simulations; D.L., M.C. and R.C. investigated feedback strategies for mitigation and containment; F.D.L., D.S., M.C., C.C., A.G. and R.C. carried out the numerical simulations with inputs from the rest of the authors; M.d.B. wrote the manuscript with inputs from F.D.R., P.D.L., A.D.M., D.L., F.L.I. and G.R.