Estimating household contact matrices structure from easily collectable metadata.

Contact matrices are a commonly adopted data representation, used to develop compartmental models for epidemic spreading, accounting for the contact heterogeneities across age groups. Their estimation, however, is generally time and effort consuming and model-driven strategies to quantify the contacts are often needed. In this article we focus on household contact matrices, describing the contacts among the members of a family and develop a parametric model to describe them. This model combines demographic and easily quantifiable survey-based data and is tested on high resolution proximity data collected in two sites in South Africa. Given its simplicity and interpretability, we expect our method to be easily applied to other contexts as well and we identify relevant questions that need to be addressed during the data collection procedure.

The limits of human mobility traces to predict the spread of COVID-19: A transfer entropy approach.

Mobile phone data have been widely used to model the spread of COVID-19; however, quantifying and comparing their predictive value across different settings is challenging. Their quality is affected by various factors and their relationship with epidemiological indicators varies over time. Here, we adopt a model-free approach based on transfer entropy to quantify the relationship between mobile phone-derived mobility metrics and COVID-19 cases and deaths in more than 200 European subnational regions. Using multiple data sources over a one-year period, we found that past knowledge of mobility does not systematically provide statistically significant information on COVID-19 spread. Our approach allows us to determine the best metric for predicting disease incidence in a particular location, at different spatial scales. Additionally, we identify geographic and demographic factors, such as users' coverage and commuting patterns, that explain the (non)observed relationship between mobility and epidemic patterns. Our work provides epidemiologists and public health officials with a general-not limited to COVID-19-framework to evaluate the usefulness of human mobility data in responding to epidemics.

Association of close-range contact patterns with SARS-CoV-2: a household transmission study.

Background: Households are an important location for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission, especially during periods when travel and work was restricted to essential services. We aimed to assess the association of close-range contact patterns with SARS-CoV-2 transmission. Methods: We deployed proximity sensors for two weeks to measure face-to-face interactions between household members after SARS-CoV-2 was identified in the household, in South Africa, 2020-2021. We calculated the duration, frequency, and average duration of close-range proximity events with SARS-CoV-2 index cases. We assessed the association of contact parameters with SARS-CoV-2 transmission using mixed effects logistic regression accounting for index and household member characteristics. Results: We included 340 individuals (88 SARS-CoV-2 index cases and 252 household members). On multivariable analysis, factors associated with SARS-CoV-2 acquisition were index cases with minimum Ct value <30 (aOR 16.8 95% CI 3.1-93.1) vs >35, and female contacts (aOR 2.5 95% CI 1.3-5.0). No contact parameters were associated with acquisition (aOR 1.0-1.1) for any of the duration, frequency, cumulative time in contact, or average duration parameters. Conclusions: We did not find an association between close-range proximity events and SARS-CoV-2 household transmission. Our findings may be due to study limitations, that droplet-mediated transmission during close-proximity contacts plays a smaller role than airborne transmission of SARS-CoV-2 in the household, or due to high contact rates in households. Funding: Wellcome Trust (Grant number 221003/Z/20/Z) in collaboration with the Foreign, Commonwealth, and Development Office, United Kingdom.

Evidence of pandemic fatigue associated with stricter tiered COVID-19 restrictions.

Despite the availability of effective vaccines against SARS-CoV-2, non-pharmaceutical interventions remain an important part of the effort to reduce viral circulation caused by emerging variants with the capability of evading vaccine-induced immunity. With the aim of striking a balance between effective mitigation and long-term sustainability, several governments worldwide have adopted systems of tiered interventions, of increasing stringency, that are calibrated according to periodic risk assessments. A key challenge remains in quantifying temporal changes in adherence to interventions, which can decrease over time due to pandemic fatigue, under such kind of multilevel strategies. Here, we examine whether there was a reduction in adherence to tiered restrictions that were imposed in Italy from November 2020 through May 2021, and in particular we assess whether temporal trends in adherence depended on the intensity of the restrictions adopted. We analyzed daily changes in movements and in residential time, combining mobility data with the restriction tier enforced in the Italian regions. Through mixed-effects regression models, we identified a general trend of reduction in adherence and an additional effect of faster waning associated with the most stringent tier. We estimated both effects being of the same order of magnitude, suggesting that adherence decreased twice as fast during the strictest tier as in the least stringent one. Our results provide a quantitative measure of behavioral responses to tiered interventions-a metric of pandemic fatigue-that can be integrated into mathematical models to evaluate future epidemic scenarios.

Interplay between mobility, multi-seeding and lockdowns shapes COVID-19 local impact.

Assessing the impact of mobility on epidemic spreading is of crucial importance for understanding the effect of policies like mass quarantines and selective re-openings. While many factors affect disease incidence at a local level, making it more or less homogeneous with respect to other areas, the importance of multi-seeding has often been overlooked. Multi-seeding occurs when several independent (non-clustered) infected individuals arrive at a susceptible population. This can lead to independent outbreaks that spark from distinct areas of the local contact (social) network. Such mechanism has the potential to boost incidence, making control efforts and contact tracing less effective. Here, through a modeling approach we show that the effect produced by the number of initial infections is non-linear on the incidence peak and peak time. When case importations are carried by mobility from an already infected area, this effect is further enhanced by the local demography and underlying mixing patterns: the impact of every seed is larger in smaller populations. Finally, both in the model simulations and the analysis, we show that a multi-seeding effect combined with mobility restrictions can explain the observed spatial heterogeneities in the first wave of COVID-19 incidence and mortality in five European countries. Our results allow us for identifying what we have called epidemic epicenter: an area that shapes incidence and mortality peaks in the entire country. The present work further clarifies the nonlinear effects that mobility can have on the evolution of an epidemic and highlight their relevance for epidemic control.

Spatial and temporal variation in proximity networks of commercial dairy cattle in Great Britain.

The nature of contacts between hosts can be important in facilitating or impeding the spread of pathogens within a population. Networks constructed from contacts between hosts allow examination of how individual variation might influence the spread of infections. Studying the contact networks of livestock species managed under different conditions can additionally provide insight into their influence on these contact structures. We collected high-resolution proximity and GPS location data from nine groups of domestic cattle (mean group size = 85) in seven dairy herds employing a range of grazing and housing regimes. Networks were constructed from cattle contacts (defined by proximity) aggregated by different temporal windows (2 h, 24 h, and approximately 1 week) and by location within the farm. Networks of contacts aggregated over the whole study were highly saturated but dividing contacts by space and time revealed substantial variation in cattle interactions. Cows showed statistically significant variation in the frequency of their contacts and in the number of cows with which they were in contact. When cows were in buildings, compared to being on pasture, contact durations were longer and cows contacted more other cows. A small number of cows showed evidence of consistent relationships but the majority of cattle did not. In one group where management allowed free access to all farm areas, cows showed asynchronous space use and, while at pasture, contacted fewer other cows and showed substantially greater between-individual variation in contacts than other groups. We highlight the degree to which variations in management (e.g. grazing access, milking routine) substantially alter cattle contact patterns, with potentially major implications for infection transmission and social interactions. In particular, where individual cows have free choice of their environment, the resulting contact networks may have a less-risky structure that could reduce the likelihood of direct transmission of infections.

Socio-economic determinants of mobility responses during the first wave of COVID-19 in Italy: from provinces to neighbourhoods.

After more than 1 year into the COVID-19 pandemic, governments worldwide still face the challenge of adopting non-pharmaceutical interventions to mitigate the risks posed by the emergence of new SARS-CoV-2 variants and the lack of a worldwide equitable vaccine allocation. Thus, it becomes crucial to identify the drivers of mobility responses to mitigation efforts during different restriction regimes, for planning interventions that are both economically and socially sustainable while effective in controlling an outbreak. Here, using anonymous and privacy-enhanced cell phone data from Italy, we investigate the determinants of spatial variations of reductions in mobility and co-location in response to the adoption and the lift of restrictions, considering both provinces and city neighbourhoods. In large urban areas, our analysis uncovers the desertification of historic city centres, which persisted after the end of the lockdown. Such centre-periphery gradient was mainly associated with differences in educational attainment. At the province level, the local structure of the labour market mainly explained the variations in mobility responses, together with other demographic factors, such as the population's age and sex composition. In the future, targeted interventions should take into account how the ability to comply with restrictions varies across geographical areas and socio-demographic groups.

Early social distancing policies in Europe, changes in mobility & COVID-19 case trajectories: Insights from Spring 2020.

BACKGROUND: Social distancing have been widely used to mitigate community spread of SARS-CoV-2. We sought to quantify the impact of COVID-19 social distancing policies across 27 European counties in spring 2020 on population mobility and the subsequent trajectory of disease. METHODS: We obtained data on national social distancing policies from the Oxford COVID-19 Government Response Tracker and aggregated and anonymized mobility data from Google. We used a pre-post comparison and two linear mixed-effects models to first assess the relationship between implementation of national policies and observed changes in mobility, and then to assess the relationship between changes in mobility and rates of COVID-19 infections in subsequent weeks. RESULTS: Compared to a pre-COVID baseline, Spain saw the largest decrease in aggregate population mobility (~70%), as measured by the time spent away from residence, while Sweden saw the smallest decrease (~20%). The largest declines in mobility were associated with mandatory stay-at-home orders, followed by mandatory workplace closures, school closures, and non-mandatory workplace closures. While mandatory shelter-in-place orders were associated with 16.7% less mobility (95% CI: -23.7% to -9.7%), non-mandatory orders were only associated with an 8.4% decrease (95% CI: -14.9% to -1.8%). Large-gathering bans were associated with the smallest change in mobility compared with other policy types. Changes in mobility were in turn associated with changes in COVID-19 case growth. For example, a 10% decrease in time spent away from places of residence was associated with 11.8% (95% CI: 3.8%, 19.1%) fewer new COVID-19 cases. DISCUSSION: This comprehensive evaluation across Europe suggests that mandatory stay-at-home orders and workplace closures had the largest impacts on population mobility and subsequent COVID-19 cases at the onset of the pandemic. With a better understanding of policies' relative performance, countries can more effectively invest in, and target, early nonpharmacological interventions.

Detecting adherence to the recommended childhood vaccination schedule from user-generated content in a US parenting forum.

Vaccine hesitancy is considered as one of the leading causes for the resurgence of vaccine preventable diseases. A non-negligible minority of parents does not fully adhere to the recommended vaccination schedule, leading their children to be partially immunized and at higher risk of contracting vaccine preventable diseases. Here, we leverage more than one million comments of 201,986 users posted from March 2008 to April 2019 on the public online forum BabyCenter US to learn more about such parents. For 32% with geographic location, we find the number of mapped users for each US state resembling the census population distribution with good agreement. We employ Natural Language Processing to identify 6884 and 10,131 users expressing their intention of following the recommended and alternative vaccination schedule, respectively RSUs and ASUs. From the analysis of their activity on the forum we find that ASUs have distinctly different interests and previous experiences with vaccination than RSUs. In particular, ASUs are more likely to follow groups focused on alternative medicine, are two times more likely to have experienced adverse events following immunization, and to mention more serious adverse reactions such as seizure or developmental regression. Content analysis of comments shows that the resources most frequently shared by both groups point to governmental domains (.gov). Finally, network analysis shows that RSUs and ASUs communicate between each other (indicating the absence of echo chambers), however with the latter group being more endogamic and favoring interactions with other ASUs. While our findings are limited to the specific platform analyzed, our approach may provide additional insights for the development of campaigns targeting parents on digital platforms.

COVID-19 outbreak response, a dataset to assess mobility changes in Italy following national lockdown.

Italy has been severely affected by the COVID-19 pandemic, reporting the highest death toll in Europe as of April 2020. Following the identification of the first infections, on February 21, 2020, national authorities have put in place an increasing number of restrictions aimed at containing the outbreak and delaying the epidemic peak. On March 12, the government imposed a national lockdown. To aid the evaluation of the impact of interventions, we present daily time-series of three different aggregated mobility metrics: the origin-destination movements between Italian provinces, the radius of gyration, and the average degree of a spatial proximity network. All metrics were computed by processing a large-scale dataset of anonymously shared positions of about 170,000 de-identified smartphone users before and during the outbreak, at the sub-national scale. This dataset can help to monitor the impact of the lockdown on the epidemic trajectory and inform future public health decision making.

weg2vec: Event embedding for temporal networks.

Network embedding techniques are powerful to capture structural regularities in networks and to identify similarities between their local fabrics. However, conventional network embedding models are developed for static structures, commonly consider nodes only and they are seriously challenged when the network is varying in time. Temporal networks may provide an advantage in the description of real systems, but they code more complex information, which could be effectively represented only by a handful of methods so far. Here, we propose a new method of event embedding of temporal networks, called weg2vec, which builds on temporal and structural similarities of events to learn a low dimensional representation of a temporal network. This projection successfully captures latent structures and similarities between events involving different nodes at different times and provides ways to predict the final outcome of spreading processes unfolding on the temporal structure.

Wearable Proximity Sensors for Monitoring a Mass Casualty Incident Exercise: Feasibility Study.

BACKGROUND: Over the past several decades, naturally occurring and man-made mass casualty incidents (MCIs) have increased in frequency and number worldwide. To test the impact of such events on medical resources, simulations can provide a safe, controlled setting while replicating the chaotic environment typical of an actual disaster. A standardized method to collect and analyze data from mass casualty exercises is needed to assess preparedness and performance of the health care staff involved. OBJECTIVE: In this study, we aimed to assess the feasibility of using wearable proximity sensors to measure proximity events during an MCI simulation. In the first instance, our objective was to demonstrate how proximity sensors can collect spatial and temporal information about the interactions between medical staff and patients during an MCI exercise in a quasi-autonomous way. In addition, we assessed how the deployment of this technology could help improve future simulations by analyzing the flow of patients in the hospital. METHODS: Data were obtained and collected through the deployment of wearable proximity sensors during an MCI functional exercise. The scenario included 2 areas: the accident site and the Advanced Medical Post, and the exercise lasted 3 hours. A total of 238 participants were involved in the exercise and classified in categories according to their role: 14 medical doctors, 16 nurses, 134 victims, 47 Emergency Medical Services staff members, and 27 health care assistants and other hospital support staff. Each victim was assigned a score related to the severity of his/her injury. Each participant wore a proximity sensor, and in addition, 30 fixed devices were placed in the field hospital. RESULTS: The contact networks show a heterogeneous distribution of the cumulative time spent in proximity by the participants. We obtained contact matrices based on the cumulative time spent in proximity between the victims and rescuers. Our results showed that the time spent in proximity by the health care teams with the victims is related to the severity of the patient's injury. The analysis of patients' flow showed that the presence of patients in the rooms of the hospital is consistent with the triage code and diagnosis, and no obvious bottlenecks were found. CONCLUSIONS: Our study shows the feasibility of the use of wearable sensors for tracking close contacts among individuals during an MCI simulation. It represents, to our knowledge, the first example of unsupervised data collection-ie, without the need for the involvement of observers, which could compromise the realism of the exercise-of face-to-face contacts during an MCI exercise. Moreover, by permitting detailed data collection about the simulation, such as data related to the flow of patients in the hospital, such deployment provides highly relevant input for the improvement of MCI resource allocation and management.

Change points, memory and epidemic spreading in temporal networks.

Dynamic networks exhibit temporal patterns that vary across different time scales, all of which can potentially affect processes that take place on the network. However, most data-driven approaches used to model time-varying networks attempt to capture only a single characteristic time scale in isolation - typically associated with the short-time memory of a Markov chain or with long-time abrupt changes caused by external or systemic events. Here we propose a unified approach to model both aspects simultaneously, detecting short and long-time behaviors of temporal networks. We do so by developing an arbitrary-order mixed Markov model with change points, and using a nonparametric Bayesian formulation that allows the Markov order and the position of change points to be determined from data without overfitting. In addition, we evaluate the quality of the multiscale model in its capacity to reproduce the spreading of epidemics on the temporal network, and we show that describing multiple time scales simultaneously has a synergistic effect, where statistically significant features are uncovered that otherwise would remain hidden by treating each time scale independently.

Estimating the outcome of spreading processes on networks with incomplete information: A dimensionality reduction approach.

Recent advances in data collection have facilitated the access to time-resolved human proximity data that can conveniently be represented as temporal networks of contacts between individuals. While the structural and dynamical information revealed by this type of data is fundamental to investigate how information or diseases propagate in a population, data often suffer from incompleteness, which possibly leads to biased estimations in data-driven models. A major challenge is thus to estimate the outcome of spreading processes occurring on temporal networks built from partial information. To cope with this problem, we devise an approach based on non-negative tensor factorization, a dimensionality reduction technique from multilinear algebra. The key idea is to learn a low-dimensional representation of the temporal network built from partial information and to use it to construct a surrogate network similar to the complete original network. To test our method, we consider several human-proximity networks, on which we perform resampling experiments to simulate a loss of data. Using our approach on the resulting partial networks, we build a surrogate version of the complete network for each. We then compare the outcome of a spreading process on the complete networks (nonaltered by a loss of data) and on the surrogate networks. We observe that the epidemic sizes obtained using the surrogate networks are in good agreement with those measured on the complete networks. Finally, we propose an extension of our framework that can leverage additional data, when available, to improve the surrogate network when the data loss is particularly large.

User-based representation of time-resolved multimodal public transportation networks.

Multimodal transportation systems, with several coexisting services like bus, tram and metro, can be represented as time-resolved multilayer networks where the different transportation modes connecting the same set of nodes are associated with distinct network layers. Their quantitative description became possible recently due to openly accessible datasets describing the geo-localized transportation dynamics of large urban areas. Advancements call for novel analytics, which combines earlier established methods and exploits the inherent complexity of the data. Here, we provide a novel user-based representation of public transportation systems, which combines representations, accounting for the presence of multiple lines and reducing the effect of spatial embeddedness, while considering the total travel time, its variability across the schedule, and taking into account the number of transfers necessary. After the adjustment of earlier techniques to the novel representation framework, we analyse the public transportation systems of several French municipal areas and identify hidden patterns of privileged connections. Furthermore, we study their efficiency as compared to the commuting flow. The proposed representation could help to enhance resilience of local transportation systems to provide better design policies for future developments.

Detecting the community structure and activity patterns of temporal networks: a non-negative tensor factorization approach.

The increasing availability of temporal network data is calling for more research on extracting and characterizing mesoscopic structures in temporal networks and on relating such structure to specific functions or properties of the system. An outstanding challenge is the extension of the results achieved for static networks to time-varying networks, where the topological structure of the system and the temporal activity patterns of its components are intertwined. Here we investigate the use of a latent factor decomposition technique, non-negative tensor factorization, to extract the community-activity structure of temporal networks. The method is intrinsically temporal and allows to simultaneously identify communities and to track their activity over time. We represent the time-varying adjacency matrix of a temporal network as a three-way tensor and approximate this tensor as a sum of terms that can be interpreted as communities of nodes with an associated activity time series. We summarize known computational techniques for tensor decomposition and discuss some quality metrics that can be used to tune the complexity of the factorized representation. We subsequently apply tensor factorization to a temporal network for which a ground truth is available for both the community structure and the temporal activity patterns. The data we use describe the social interactions of students in a school, the associations between students and school classes, and the spatio-temporal trajectories of students over time. We show that non-negative tensor factorization is capable of recovering the class structure with high accuracy. In particular, the extracted tensor components can be validated either as known school classes, or in terms of correlated activity patterns, i.e., of spatial and temporal coincidences that are determined by the known school activity schedule.

Activity clocks: spreading dynamics on temporal networks of human contact.

Dynamical processes on time-varying complex networks are key to understanding and modeling a broad variety of processes in socio-technical systems. Here we focus on empirical temporal networks of human proximity and we aim at understanding the factors that, in simulation, shape the arrival time distribution of simple spreading processes. Abandoning the notion of wall-clock time in favour of node-specific clocks based on activity exposes robust statistical patterns in the arrival times across different social contexts. Using randomization strategies and generative models constrained by data, we show that these patterns can be understood in terms of heterogeneous inter-event time distributions coupled with heterogeneous numbers of events per edge. We also show, both empirically and by using a synthetic dataset, that significant deviations from the above behavior can be caused by the presence of edge classes with strong activity correlations.

Schelling segregation in an open city: a kinetically constrained Blume-Emery-Griffiths spin-1 system.

In the 70s Schelling introduced a multiagent model to describe the segregation dynamics that may occur with individuals having only weak preferences for "similar" neighbors. Recently variants of this model have been discussed, in particular, with emphasis on the links with statistical physics models. Whereas these models consider a fixed number of agents moving on a lattice, here, we present a version allowing for exchanges with an external reservoir of agents. The density of agents is controlled by a parameter which can be viewed as measuring the attractiveness of the city lattice. This model is directly related to the zero-temperature dynamics of the Blume-Emery-Griffiths spin-1 model, with kinetic constraints. With a varying vacancy density, the dynamics with agents making deterministic decisions leads to a variety of "phases" whose main features are the characteristics of the interfaces between clusters of agents of different types. The domains of existence of each type of interface are obtained analytically as well as numerically. These interfaces may completely isolate the agents leading to another type of segregation as compared to what is observed in the original Schelling model, and we discuss its possible socioeconomic correlates.