Epidemic spread on networks is one of the most studied dynamics in network science and has important implications in real epidemic scenarios. Nonetheless, the dynamics of real epidemics and how it is affected by the underline structure of the infection channels are still not fully understood. Here we apply the susceptible-infected-recovered model and study analytically and numerically the epidemic spread on a recently developed spatial modular model imitating the structure of cities in a country. The model assumes that inside a city the infection channels connect many different locations, while the infection channels between cities are less and usually directly connect only a few nearest neighbor cities in a two-dimensional plane. We find that the model experience two epidemic transitions. The first lower threshold represents a local epidemic spread within a city but not to the entire country and the second higher threshold represents a global epidemic in the entire country. Based on our analytical solution we proposed several control strategies and how to optimize them. We also show that while control strategies can successfully control the disease, early actions are essentials to prevent the disease global spread.

We study the critical features of the order parameter’s fluctuations near the threshold of mixed-order phase transitions in randomly interdependent spatial networks. Remarkably, we find that although the structure of the order parameter is not scale invariant, its fluctuations are fractal up to a well-defined correlation length ${\xi }^{\prime }$ that diverges when approaching the mixed-order transition threshold. We characterize the self-similar nature of these critical fluctuations through their effective fractal dimension $d_f^{\prime }=3d/4$, and correlation length exponent ${\nu }^{\prime }=2/d$, where $d$ is the dimension of the system. By analyzing percolation and magnetization, we demonstrate that $d_f^{\prime }$ and ${\nu }^{\prime }$ are the same for both, i.e., independent of the symmetry of the process for any $d$ of the underlying networks.

The new coronavirus known as COVID-19 is spread world-wide since December 2019. Without any vaccination or medicine, the means of controlling it are limited to quarantine and social distancing. Here we study the spatio-temporal propagation of the first wave of the COVID-19 virus in China and compare it to other global locations. We provide a comprehensive picture of the spatial propagation from Hubei to other provinces in China in terms of distance, population size, and human mobility and their scaling relations. Since strict quarantine has been usually applied between cities, more insight about the temporal evolution of the disease can be obtained by analyzing the epidemic within cities, especially the time evolution of the infection, death, and recovery rates which affected by policies. We study and compare the infection rate in different cities in China and provinces in Italy and find that the disease spread is characterized by a two-stages process. At early times, at order of few days, the infection rate is close to a constant probably due to the lack of means to detect infected individuals before infection symptoms are observed. Then at later times it decays approximately exponentially due to quarantines. The time evolution of the death and recovery rates also distinguish between these two stages and reflect the health system situation which could be overloaded.

Many interdependent, real-world infrastructures involve interconnections between different communities or cities. Here we show how the effects of such interconnections can be described as an external field for interdependent networks experiencing a first-order percolation transition. We find that the critical exponents $\gamma$ and $\delta$, related to the external field, can also be defined for first-order transitions but that they have different values than those found for second-order transitions. Surprisingly, we find that both sets of different exponents (for first and second order) can even be found within a single model of interdependent networks, depending on the dependency coupling strength. Nevertheless, in both cases both sets satisfy Widom's identity, $\delta -1=\gamma /\beta$, which further supports the validity of their definitions. Furthermore, we find that both Erdős-Rényi and scale-free networks have the same values of the exponents in the first-order regime, implying that these models are in the same universality class. In addition, we find that in $k$-core percolation the values of the critical exponents related to the field are the same as for interdependent networks, suggesting that these systems also belong to the same universality class.