Abstract
Background The first wave of the COVID-19 pandemic in 2020 was largely mitigated by reducing contacts in the general population. In 2022 most contact-reducing measures were lifted.
Aim We assess whether the population has reverted to pre-pandemic contact behaviour and how this would affect the transmission potential of a newly emerging pathogen.
Methods The PIENTER Corona study was held every 2-6 months in the Netherlands from April 2020, as a follow-up on the 2016-2017 PIENTER3 study. In both studies, participants (ages 1-85) reported the number and age group of all face-to-face persons contacted on the previous day. The contact behaviour during and after the COVID-19 pandemic was compared to the pre-pandemic baseline.
Results We found an average of 15.2 (13.3-16.9, 95% CI) community contacts per person per day in the post-pandemic period, which is 14% lower than the baseline value of 17.6 (16.3-18.9). Children have the highest number of contacts as before the pandemic. Mainly adults aged 20-59 have not reverted to their pre-pandemic behaviour, possibly because this age group works more often from home. Although the number of contacts is structurally lower compared to the pre-pandemic period, the effect on the potential spread of a newly emerging respiratory pathogen is limited if all age groups were equally susceptible. If younger age groups were less susceptible, as observed during the first COVID-19 wave, the transmission potential as well as the required control effort would be lower.
Conclusion Continuous monitoring of contacts is needed to be prepared for a future pandemic.
Introduction
After the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged at the end of 2019, the World Health Organization declared the outbreak a global pandemic on 11 March 2020. Most nations implemented stringent non-pharmaceutical interventions to mitigate the virus spread, encompassing physical distancing aimed at reducing exposure to SARS-CoV-2 by lowering contact rates. Face-to-face and physical contacts between members of different households were to be reduced, as each contact could be an at-risk event for transmission. Many contact surveys that were conducted during the pandemic showed that more stringent COVID-19 measures indeed lead to lower contact rates [1–8]. Most contact-reducing measures were lifted shortly after the emergence of the more transmissible Omicron variants of SARS-CoV-2 in late 2021. The question arises whether contacts in the population reverted to their pre-pandemic levels and how this would affect a future pandemic.
Answering these questions would require a contact study in the post-pandemic period, including all age groups, and a valid baseline measurement before the pandemic started. The only study we are aware of that fits these requirements is the PIENTER Corona (PiCo) contact study conducted in the Netherlands [9]. The study’s design closely mirrors that of the PIENTER3 study conducted in 2016-2017 [10], providing the pre-pandemic baseline. Both studies were part of a larger serosurveillance study encompassing a large nationwide sample of the Dutch population. The PiCo survey was conducted every 2-6 months since April 2020 and will continue until the end of 2024. An earlier analysis revealed a reduction of 76% and 41% in the number of contacts outside households during the April 2020 and June 2020 survey rounds, respectively, in comparison to the baseline survey [11].
A change in contact patterns between groups could also affect the spread of a future pandemic. During the COVID-19 pandemic, the change in behaviour varied among different groups in the population, either because they were differently affected by control measures or because they differed in ability or willingness to comply. For instance, contacts in younger age groups were greatly affected by school closures [12], while contacts in older age groups decreased due to suspended social gatherings [6,13]. Persons with comorbidities may have reduced their contacts more than persons without comorbidities [6,14,15], given their higher probability of severe outcomes [16]. Also the socio-economic status of a person may have influenced contact behaviour, as persons with a low socio-economic status often have professions that preclude working from home [17].
On 25 February 2022 all contact-reducing measures were lifted in the Netherlands. We analyze the contact behaviour in the PiCo survey rounds before and after this date, compared to the contact behaviour before the pandemic. We study the differences in contacts between age groups, medical risk groups and socio-economic groups, and we assess the potential implications for a future outbreak of a respiratory pathogen.
Methods
Study design and data
The baseline survey was conducted in a representative sample of the Dutch population aged between 0 to 89 years old [10] from February 2016 to October 2017. Most of these participants were invited for the first round of the PiCo survey conducted in April 2020 [9]. The study population was supplemented with a random selection from the Dutch population older than one year in PiCo rounds 2 (June 2020) [18] and 6 (November 2021) [19]. For comparability the baseline participants younger than one year old were excluded from the analysis.
The questionnaires were filled out either by the participant or with help from their parent or guardian if they were younger than 15 years. Participant characteristics included age and sex, as well as the age and sex of all household members. From a list of self-reported medical conditions, the medical risk status of each participant was determined. Participants with an indication for influenza vaccination [20] were classified as having a high medical risk status. Participants provided their highest obtained or current education level, that served as a proxy for socio-economic status. The education level was classified according to Dutch standards [21] as low (no education or primary education), medium (secondary school or vocational training), or high (bachelor’s degree, university). Participants up to 14 years old were assigned the highest education level of their parents or guardians. Participants were also asked whether they worked from home in the previous week.
The questionnaire included a section on contact behaviour, where contacts were defined as unique persons with whom the participant talked face-to-face, touched, kissed or played sports. Participants listed the number of contacts they had with persons outside their household on the previous day. The contacts were stratified by age group and either by gender (male/female) or by proximity (close/distant contacts). Participants that failed to provide valid contact data were excluded from the analysis.
A detailed description on the definition of the medical risk status, the validity of contact data, and handling of changes in survey questions is provided in Supplement S1. The data was reformatted to the standard format for contact surveys on socialcontactdata.org and published online [22].
Contact analysis
The age groups of participants and contacts were aggregated in 10 age groups: 0–4, 5–9, 10–19, 20–29, 30–39, 40–49, 50–59, 60–69, 70–79, and ≥ 80. In each round, the number of community contacts per participant per contact age group was truncated at 50. The household composition of the participant was used as a proxy for contacts within the household, assuming the participant would have contacted each household member on the survey day. The community and household contacts were analysed separately.
Within each participant age group, participants were weighted according to the age and sex distribution of the Dutch population [23]. Participants were also weighted according to weekday (weight 5/7) or weekend (weight 2/7), but only for community contacts, assuming the household composition is constant over the week. With these sample weights the weighted mean number of contacts per participant in that participant age group and contact age group was calculated. These numbers are the elements of a contact matrix. Assuming all contacts are reciprocal, we corrected the contact matrix for reporting errors using the age distribution of the Dutch population [24]. The uncertainty was expressed by the 95% bias-corrected bootstrap interval [25], based on 1000 bootstrap samples by participant age group.
The number of community contacts per participant per day was compared to the baseline values that were obtained from the 2016-2017 contact survey. We studied the number of contacts over time and by age group. The reported contacts by medical risk group and education level were expected to be confounded by the age of the participants. To be able to compare the subgroups, we weighted their age groups by the age distribution of the general population and calculated the population average. Such a population average of for instance the high medical risk group should therefore be interpreted as the population average if the entire population would have a high medical risk.
A total contact matrix was constructed by summing the contact matrices of the community and household contacts. If all age groups were equally infectious and susceptible to a newly emerging pathogen, this matrix can be interpreted as the Next Generation Matrix (NGM) of an epidemic process [26]. In a fully susceptible population without any control measures, the spectral radius (i.e., maximum eigenvalue) of the NGM is proportional to the basic reproduction number R0, indicating the average number of secondary cases infected by a typical primary case [27]. We compared the spectral radius over time to the baseline value. The rounds without any control measures reflect the transmission potential of a newly emerging pathogen under ‘new normal’ circumstances compared to the pre-COVID-19 period. We also determined the transmission potential if the susceptibility differed by age group, as was demonstrated for SARS-CoV-2 in the early phase of the pandemic [1].
Results
Study population
From the baseline and PiCo surveys, 387 participants younger than one year, and 1780 questionnaires without valid contact data were excluded. The final data set comprises 62490 questionnaires, filled out by 13826 unique participants, 97% of which also provided their household composition. In total 5768 persons participated in the baseline survey (round 0 in Tab. 1), 2487 of who also participated in one or more rounds of the PiCo survey. The number of participants in the PiCo survey ranged from 2594 in the first round when only baseline participants were reinvited to 8144 in the sixth round when the study population was supplemented.
Study population, by survey month, the number of participants per survey round and stratified in percentages by age group, sex, household size, medical risk group and education level. The baseline survey from 2016-2017 is indicated as survey round 0. The final column contains reference percentages for the general population by age group [23], sex [23], household size [28], medical risk status [29], and education level for 15-90 year olds [21].
Generally, age groups under 40 years were underrepresented, and age groups between 60 and 80 years were overrepresented compared to the Dutch population [23]. This imbalance seemed to increase over time, suggesting a higher drop-out rate for younger age groups. The study population contained around 56% female participants; a percentage which was constant over all rounds. Participants living in two-person households were overrepresented [28]. This was partly explained by the higher age of the participants; when the household size frequency was corrected for confounding by age and sex, two-person households were less overrepresented, but single-person households remained underrepresented (Suppl S2). The fraction of participants with a high medical risk was higher than in the general population [29], but this was completely explained by confounding with age and sex (Fig. 1A and Suppl S2). The education level of the participants was higher than in the general population [21] (Fig. 1B). This overrepresentation became even more pronounced after correcting for age and sex confounding (Suppl S2), because the abundant older participants had on average a lower education level.
Participant characteristics of the study population by age group. (A) Fraction of participants by medical risk status, for unique participants from PiCo round 4 (March 2021) onwards. Reference values (black lines) are taken from [29] (B) Fraction of participants by education level, for unique participants. For participants under 15, the highest education level of the parents or guardians is used. For comparability the age group 10-19 is split into 10-14 and 15-19. Reference values (black lines) are taken from [21].
Number of community contacts
For 1262 questionnaires, the number of community contacts was truncated as they exceeded the maximum of 50 contacts per contact age group. Because of this truncation the 95th percentile of the number of community contacts per participant decreased from 59 to 54.
The number of community contacts in the general population varied over the course of the COVID-19 pandemic (Fig. 2). Numbers were low in periods with stringent measures and high numbers of hospital admissions. Around the start of 2022 when the Omicron variant became dominant, a short period of lockdown measures did not seem to affect the average number of contacts as much, because 99% of questionnaires in survey round 6 was filled out before these measures came into effect on 19 December 2021. Since the physical distancing measures had been lifted on 25 February 2022, the number of community contacts increased from 13.6 (12.1-15.1) community contacts per person per day in survey round 7 to 15.2 (13.3-16.9) in survey rounds 9 and 10 taken together. However, these latest numbers were still 14% lower than the baseline value of 17.6 (16.3-18.9) contacts per day.
Weighted average number of community contacts per person per day (pppd) in the general population, for the baseline and PiCo survey. Shown is the weighted average number of contacts (horizontal line) and 95% bias-corrected bootstrap interval (shaded area). PiCo rounds are shown from the start to end date, with the median survey date (vertical line) with the survey round number in white. As a timeline reference the weekly number of hospital admissions [38] and the Oxford Stringency Index [39] are depicted on top. The stringency index is color coded from blue (low) to red (high).
When examining the results in more detail, we observed large differences between age groups (Fig. 3). In the baseline survey, the number of contacts was highest for the 5-9 age group and gradually decreased with age. During the pandemic, the age groups under 10 years of age were least affected in their contact behaviour. With the exception of the lockdown periods that included primary school and day care closures, their number of contacts was around the baseline value. Compared to them, the contacts of 10-19 year olds were more heavily affected, but they have reverted to pre-pandemic levels from 2022 onwards. For all other age groups, the number of contacts were well below baseline during the first two pandemic years. In the last survey rounds, the age groups of 60 years and older approached pre-pandemic behaviour, but the 20-59 year olds have settled at a below-baseline level. The difference between the modelled (rectangles in Fig. 3) and reported (circles in Fig. 3) number of contacts shows how much the reported contacts needed to be adjusted to comply with the assumption of reciprocal contacts. For instance, the 5-9 year olds reported fewer contacts with other age groups, than other age groups with them. Generally, differences between modelled and observed numbers of contacts were small, indicating consistent reporting behaviour.
Average number of community contacts per participant per day (pppd) in the general population by age group, for the baseline and PiCo survey. Shown is the average number of contacts (horizontal line) and 95% bias-corrected bootstrap interval (shaded area). PiCo rounds are shown from the start to end date, with the median survey date (vertical line). For the PiCo rounds the reported number of contacts is shown (open circle) with the 95% bias-corrected bootstrap interval (whisker) where the size of the circle is scaled by the number of participants.
In the baseline survey the average number of community contacts did not differ by medical risk group (Fig. 4). Also in the PiCo rounds no difference between the two medical risk groups was observed.
Weighted average number of community contacts per person per day (pppd) by medical risk group (left) and education level (right) for the PiCo survey. The averages are calculated by weighting the age groups in each stratum by the age distribution of the general population. Shown is the average number of contacts (point) and 95% bias-corrected bootstrap interval (whisker). For the baseline survey only the average number of contacts is indicated (dashed line).
In the baseline survey, the average number of community contacts increased with education level (Fig. 4). In the PiCo rounds with the most stringent measures (rounds 1 and 4), this ranking was reversed. In all other rounds, with less stringent measures, the average number of community contacts for the low education level is lowest, while middle and high education levels have a similar higher number of contacts.
Transmission potential
To study the impact of the post-pandemic contact behaviour on the transmission potential of a new respiratory pathogen, we summed the contact matrices of the community and household contacts and assumed values for the relative susceptibility and infectiousness per age group. When all age groups are equally susceptible and infectious, the spectral radius of the resulting matrix is around the baseline value during most of the study period, except for survey rounds 1 and 4 when measures were most stringent (Fig. 5). When all age groups are equally susceptible, transmission is driven by the age groups with most contacts, which are the youngest age groups that exhibited baseline behaviour during most of the survey rounds. For the early phase of COVID-19, Zhang et al. [1] estimated that younger age groups are less susceptible and older age groups are more susceptible to SARS-CoV-2 infection than adult age groups. When younger age groups are less susceptible, the transmission potential is generally lower than the baseline value, as transmission is now more driven by adult age groups that did structurally decrease their number of contacts. Other reported patterns of the relative susceptibility and infectiousness of COVID-19 lead to similar results (Suppl. S3).
Transmission potential, expressed as the spectral radius of the next generation matrix (NGM) compared to the baseline value, assuming equal relative susceptibility and infectiousness over age groups, or a susceptibility profile for the early phase of COVID-19 [1]. Shown is the ratio (horizontal line) and 95% bias-corrected bootstrap interval (shaded area). PiCo rounds are shown from the start to end date, with the median survey date (vertical line).
Discussion
In this study we evaluated whether the Dutch general population has reverted to its pre-pandemic contact behaviour after lifting the COVID-19 measures. Our results show that the population average of community contacts had decreased by 14% in the survey rounds in November 2022 and May 2023 when all control measures had been lifted, compared to the baseline value in 2016-2017. However, children did revert to the high number of contacts they made before the pandemic, retaining their role as potential epidemic drivers. As a consequence, the transmission potential of a newly emerging infection that would spread by a respiratory route mainly depends on how susceptible and infectious children are for this new pathogen.
The strengths of this study include the identical study design of the PiCo survey and the PIENTER3 survey, which provides a reliable baseline to compare current contact behaviour to. Both surveys are conducted with participants of all age groups, which is essential to study the spread in the entire population. Another strength is that the contact questionnaire is embedded in a larger serosurveillance study, providing a wealth of additional information on the participants. This allowed us to stratify our analysis by medical risk group and education level.
Limitations include the possibility that the study population was on average more compliant to COVID-19 measures than the general population. Evidence for this is for instance their higher COVID-19 vaccine uptake [19]. Also, because of COVID-19 information campaigns, participants may have reported contacts differently compared to the baseline survey, because they would now have a better understanding on what constitutes a contact. Finally, we have weighted participants by age and sex, but not by education level due to the lack of reference values for participants younger than 15 years. Since participants with a high education level are overrepresented, and tend to have a higher number of contacts than average with relaxed measures and a lower number of contacts than average with stringent measures, a consequence could be that the numbers of contacts are overestimated in periods with relaxed measures and underestimated in periods with stringent measures.
We found that for children under 20 years and adults over 60 years the number of contacts have reverted to their baseline levels. Only adults in the 20-59 age groups seem to have structurally fewer contacts. A possible explanation is that employees work more hours from home [30] and make fewer contacts while working from home. This trend is further substantiated by the results of the survey question whether participants worked from home last week. For all working participants in age groups 20-69 the fraction (partly) working from home has approximately doubled in the last two survey rounds, compared to a baseline value (Suppl. S4). This is in line with a different survey that reported 37% of employees worked partly from home in 2019, increasing to 45% in 2022 [31].
We found that in the baseline survey the average number of contacts per participant increased with education level. This difference was smaller or even reversed in the periods with lockdown measures (survey rounds 1 and 4). An explanation is that participants with a lower education level more often have work in essential sectors that preclude working from home [17]. This is confirmed by the results on the survey question whether participants worked from home last week stratified by education level (Suppl. S4). During the lockdown periods around 70% of the highly educated participants worked (partly) from home, against 30% of the participants with a low education level. In the last two survey rounds these percentages dropped to 50% and 25% respectively, which is still higher than the baseline values of 30% and 10% respectively. We found that the number of contacts did not differ between medical risk groups (i.e. indication for influenza vaccination) in any of the survey rounds nor the baseline survey. This could be because of our broad risk definition; persons with a very high medical risk may well have reduced their number of contacts. Participants with a high medical risk could also have had safer contacts with reduced risk of transmission by keeping distance and/or wearing protection. However, for these types of contacts we found no difference between medical risk groups (Suppl. S5).
A higher number of contacts could translate into a higher exposure to infection which, if all other factors would remain equal, could result in a higher expected infection attack rate [18,32]. For the PiCo participants with a low education level, who had the highest number of contacts in lockdown periods with high incidence, this is true for the pre-Omicron period [19]. For the PiCo participants who received an influenza vaccination a lower seroprevalence due to infection was reported [19]. However, this group does not align with our high medical risk group, as the step from being indicated to receiving an influenza vaccination also depends on risk perception and compliance.
A direct comparison of our results with other studies is difficult as, to the best of our knowledge among all contact studies only the CONNECT study [7,14] can make use of a pre-pandemic baseline. Some of our findings are in line with earlier reports. The effect of education level on the difference in number of contacts during stringent and relaxed measures was also reported for a contact study in Hungary [33]. The absence of an association between number of contacts and medical risk group (defined as indication for influenza vaccination) was found in the COVIMOD study in Germany [15]. In one of the six rounds of the CONNECT study in Canada, it was found that persons with a higher medical risk have fewer contacts, but this was also the case for their baseline survey [14]. This would suggest their risk definition could differ from ours.
Our finding that the different contact patterns have a limited impact on the spread of a new pathogen if all age groups are equally susceptible is in apparent contrast with the findings reported in the final CoMix round held in November 2022 in four European countries [34]. For the Netherlands, an average number of 9.9 (9.0-10.8) contacts per person per day was found, which is lower than the 14.8 (13.4-16.2) contacts per person per day that we found in November 2022 (PiCo round 9). Using an earlier survey from 2006-2007 as baseline [35], the CoMix study reported that the transmission potential in the Netherlands had decreased by 20%, while our study suggests no change in transmission potential when assuming equal susceptibility and infectiousness over age groups. Possible explanations could be found in the different design, study population and study period of the baseline survey, and in the fatigue effect observed for the Dutch CoMix study which leads participants to report decreasing numbers of contacts due to the high frequency of reporting [4,36].
Contact surveys such as PIENTER3 and PiCo have proven invaluable during the COVID-19 pandemic in monitoring changes in contact behaviour due to evolving measures or compliance with them. They have provided direct quantitative evidence to inform tailored infection control policies. Continuation of such contact surveys is therefore essential, even in a non-pandemic setting, not only to signal changes in contact patterns that affect the transmission potential, but also to provide a baseline for the next pandemic.
Declarations
Ethics approval and consent
This study was approved by the Medical Ethics Committee (MEC-U), the Netherlands (Clinical Trial Registration NTR8473), conformed to the principles embodied in the Declaration of Helsinki, and all participants provided written informed consent.
Availability of data and materials
The cleaned data is published online at https://doi.org/10.5281/zenodo.10370353 [22].
All code for data cleaning and analysis is publicly available at https://github.com/rivm-syso/pico_contactsurvey [37].
Competing interests
The authors have no competing interests, or other interests that might be perceived to influence the results and/or discussion reported in this paper.
Funding
This study was funded by the Ministry of Health, Welfare and Sport (VWS) in the Netherlands. The authors of this study (JB and JW) received funding from European Union’s Horizon 2020 research and innovation programme -projects EpiPose (Grant agreement number 101003688) and ESCAPE (Grant agreement number 101095619).
Authors’ contributions
F.v.d.K, G.d.H, H.d.M. and J.W. designed the study, F.v.d.K, G.d.H, C.v.d.H and E.V. collected the data, J.B. did the analyses, J.B. wrote the main manuscript text. All authors reviewed the manuscript.
Acknowledgments
We greatly acknowledge the participants of the PIENTER3 and PIENTER Corona (PICO) studies. We would also like to thank Federica Giardini (Radboud University Nijmegen), Ganna Rozhnova (Utrecht University) and Otilia Boldea (Tilburg University) for their valuable input.