Search

We developed an endemic model of COVID-19 to assess the impact of vaccination and immunity waning on the dynamics of the disease. Our model exhibits the phenomenon of backward bifurcation and bi-stability, where a stable disease-free equilibrium coexists with a stable endemic equilibrium. The epidemiological implication of this is that the control reproduction number being less than unity is no longer sufficient to guarantee disease eradication. We showed that this phenomenon could be eliminated by either increasing the vaccine efficacy or by reducing the disease transmission rate (adhering to non-pharmaceutical interventions). Furthermore, we numerically investigated the impacts of vaccination and waning of both vaccine-induced immunity and post-recovery immunity on the disease dynamics. Our simulation results show that the waning of vaccine-induced immunity has more effect on the disease dynamics relative to post-recovery immunity waning and suggests that more emphasis should be on reducing the waning of vaccine-induced immunity to eradicate COVID-19.

The outbreak of the severe acute respiratory syndrome coronavirus 2 started in Wuhan, China, towards the end of 2019 and spread worldwide. The rapid spread of the disease can be attributed to many factors including its high infectiousness and the high rate of human mobility around the world. Although travel/movement restrictions and other non-pharmaceutical interventions aimed at controlling the disease spread were put in place during the early stages of the pandemic, these interventions did not stop COVID-19 spread. To better understand the impact of human mobility on the spread of COVID-19 between regions, we propose a hybrid gravity-metapopulation model of COVID-19. Our modeling framework has the flexibility of determining mobility between regions based on the distances between the regions or using data from mobile devices. In addition, our model explicitly incorporates time-dependent human mobility into the disease transmission rate, and has the potential to incorporate other factors that affect disease transmission such as facemasks, physical distancing, contact rates, etc. An important feature of this modeling framework is its ability to independently assess the contribution of each factor to disease transmission. Using a Bayesian hierarchical modeling framework, we calibrate our model to the weekly reported cases of COVID-19 in thirteen local health areas in Metro Vancouver, British Columbia (BC), Canada, from July 2020 to January 2021. We consider two main scenarios in our model calibration: using a fixed distance matrix and time-dependent weekly mobility matrices. We found that the distance matrix provides a better fit to the data, whilst the mobility matrices have the ability to explain the variance in transmission between regions. This result shows that the mobility data provides more information in terms of disease transmission than the distances between the regions.

Adherence to public health policies such as the non-pharmaceutical interventions implemented against COVID-19 plays a major role in reducing infections and controlling the spread of the diseases. In addition, understanding the transmission dynamics of the disease is also important in order to make and implement efficient public health policies. In this paper, we developed an SEIR-type compartmental model to assess the impact of adherence to COVID-19 non-pharmaceutical interventions and indirect transmission on the dynamics of the disease. Our model considers both direct and indirect transmission routes and stratifies the population into two groups: those that adhere to COVID-19 non-pharmaceutical interventions (NPIs) and those that do not adhere to the NPIs. We compute the control reproduction number and the final epidemic size relation for our model and study the effect of different parameters of the model on these quantities. Our results show that direct transmission has more effect on the reproduction number and final epidemic size, relative to indirect transmission. In addition, we showed that there is a significant benefit in adhering to the COVID-19 NPIs.

Monkeypox is an emerging zoonotic disease caused by the monkeypox virus, which is an infectious agent belonging to the genus Orthopoxvirus. Currently, commencing from the end of April 2022, an outbreak of monkeypox is ongoing, with more than 43,000 cases reported as of 23 August 2022, involving 99 countries and territories across all the six World Health Organization (WHO) regions. On 23 July 2022, the Director-General of the WHO declared monkeypox a global public health emergency of international concern (PHEIC), since the outbreak represents an extraordinary, unusual, and unexpected event that poses a significant risk for international spread, requiring an immediate, coordinated international response. However, the real magnitude of the burden of disease could be masked by failures in ascertainment and under-detection. As such, underestimation affects the efficiency and reliability of surveillance and notification systems and compromises the possibility of making informed and evidence-based policy decisions in terms of the adoption and implementation of ad hoc adequate preventive measures. In this review, synthesizing 53 papers, we summarize the determinants of the underestimation of sexually transmitted diseases, in general, and, in particular, monkeypox, in terms of all their various components and dimensions (under-ascertainment, underreporting, under-detection, under-diagnosis, misdiagnosis/misclassification, and under-notification).

A hybrid asymptotic-numerical method is developed to approximate the mean first passage time (MFPT) and the splitting probability for a Brownian particle in a bounded two-dimensional (2D) domain that contains absorbing disks, referred to as "traps”, of asymptotically small radii. In contrast to previous studies that required traps to be spatially well separated, we show how to readily incorporate the effect of a cluster of closely spaced traps by adapting a recently formulated least-squares approach in order to numerically solve certain local problems for the Laplacian near the cluster. We also provide new asymptotic formulae for the MFPT in 2D spatially periodic domains where a trap cluster is centred at the lattice points of an oblique Bravais lattice. Over all such lattices with fixed area for the primitive cell, and for each specific trap set, the average MFPT is smallest for a hexagonal lattice of traps. doi:10.1017/S1446181121000018

The determination of the mean first passage time (MFPT) for a Brownian particle in a bounded 2-D domain containing small absorbing traps is a fundamental problem with biophysical applications. The average MFPT is the expected capture time assuming a uniform distribution of starting points for the random walk. We develop a hybrid asymptotic-numerical approach to predict optimal configurations of $m$ small stationary circular absorbing traps that minimize the average MFPT in near-disk and elliptical domains. For a general class of near-disk domains, we illustrate through several specific examples how simple, but yet highly accurate, numerical methods can be used to implement the asymptotic theory. From the derivation of a new explicit formula for the Neumann Green's function and its regular part for the ellipse, a numerical approach based on our asymptotic theory is used to investigate how the spatial distribution of the optimal trap locations changes as the aspect ratio of an ellipse of fixed area is varied. The results from the hybrid theory for the ellipse are compared with full PDE numerical results computed from the closest point method \cite{IWWC2019}. For long and thin ellipses, it is shown that the optimal trap pattern for $m=2,\ldots,5$ identical traps is collinear along the semi-major axis of the ellipse. For such essentially 1-D patterns, a thin-domain asymptotic analysis is formulated and implemented to accurately predict the optimal locations of collinear trap patterns and the corresponding optimal average MFPT., 28 pages, 41 figures

We analyze a class of cell-bulk coupled PDE-ODE models, motivated by quorum-sensing and diffusion-mediated behavior in microbial systems, that characterize communication between localized spatially...

We extended a class of coupled PDE-ODE models for studying the spatial spread of airborne diseases by incorporating human mobility. Human populations are modeled with patches, and a Lagrangian perspective is used to keep track of individuals’ places of residence. The movement of pathogens in the air is modeled with linear diffusion and coupled to the SIR dynamics of each human population through an integral of the density of pathogen around the population patch. In the limit of fast diffusion pathogens, the method of matched asymptotic analysis is used to reduce the coupled PDE-ODE model to a nonlinear system of ODEs for the average density of pathogens in the air. The reduced system of ODEs is used to derive the basic reproduction number and the final size relation for the model. Numerical simulations of the full PDE-ODE model and the reduced system of ODEs are used to assess the impact of human mobility, together with the diffusion of pathogens on the dynamics of the disease. Results from the two models are consistent and show that human mobility significantly affects disease dynamics. In addition, we show that an increase in the diffusion rate of pathogen leads to a smaller epidemic.

PurposeClose-contact rates are thought to be a significant driving force behind the dynamics of transmission for many infectious respiratory diseases. Efforts to control such infections typically focus on the practice of strict contact-avoidance measures. Yet, contact rates and their relation to transmission, and the impact of control measures, are seldom quantified. Here, we quantify the response of contact rates, transmission and new cases of COVID-19 to public health contact-restriction orders, and the associations among these three variables, in the Canadian province of British Columbia (BC) and within its two most densely populated regional health authorities: Fraser Health Authority (FHA) and Vancouver Coastal Health Authority (VCHA).MethodsWe obtained time series for self-reported close-contact rates from the BC Mix COVID-19 Survey, new reported cases of COVID-19 from the BC Center for Disease Control, and transmission rates based on dynamic model fits to reported cases. Our study period was from September 13, 2020 to February 19, 2021, during which three public health contact-restriction orders were introduced (October 26, November 7 and November 19, 2020). We used segmented linear regression to quantify impacts of public health orders, Pearson correlation to assess the instantaneous relation between contact rates and transmission, and vector autoregressive modeling to study the lagged relations among the three variables.ResultsOverall, declines in contact rates and transmission occurred concurrently with the announcement of public health orders, whereas declines in new cases showed a reporting delay of roughly two weeks. The impact of the first public health order (October 26, 2020) on contact rates and transmission was more pronounced than that of the other two health orders. Contact rates and transmission on the same day were strongly correlated (correlation coefficients = 0.64, 0.53 and 0.34 for BC, FHA, and VCHA, respectively). Moreover, contact rates were a significant time-series driver of COVID-19 and explained roughly 30% and 18% of the variation in new cases and transmission, respectively. Interestingly, increases in transmission and new cases were followed by reduced rates of contact: overall, average daily cases explained about 10% of the variation in provincial contact rates.ConclusionWe show that close-contact rates were a significant driver of transmission of COVID-19 in British Columbia, Canada and that they varied in response to public health orders. Our results also suggest a possible feedback, by which contact rates respond to recent changes in reported cases. Our findings help to explain and validate the commonly assumed, but rarely measured, response of close contact rates to public health guidelines and their impact on the dynamics of infectious diseases.

We analyse oscillatory instabilities for a coupled partial-ordinary differential equation (PDE-ODE) system modelling the communication between localised spatially segregated dynamically active signalling compartments that are coupled through a passive extracellular bulk diffusion field in a bounded 2D domain. Each signalling compartment is assumed to secrete a chemical into the extracellular medium (bulk region), and it can also sense the concentration of this chemical in the region around its boundary. This feedback from the bulk region, resulting from the entire collection of cells, in turn modifies the intracellular dynamics within each cell. In the limit where the signalling compartments are circular discs with a small common radius ɛ ≪ 1 and where the bulk diffusivity is asymptotically large, a matched asymptotic analysis is used to reduce the dimensionless PDE–ODE system into a nonlinear ODE system with global coupling. For Sel’kov reaction kinetics, this ODE system for the intracellular dynamics and the spatial average of the bulk diffusion field are then used to investigate oscillatory instabilities in the dynamics of the cells that are triggered due to the global coupling. In particular, numerical bifurcation software on the ODEs is used to study the overall effect of coupling defective cells (cells that behave differently from the remaining cells) to a group of identical cells. Moreover, when the number of cells is large, the Kuramoto order parameter is computed to predict the degree of phase synchronisation of the intracellular dynamics. Quorum sensing behaviour, characterised by the onset of collective behaviour in the intracellular dynamics as the number of cells increases above a threshold, is also studied. Our analysis shows that the cell population density plays a dual role of triggering and then quenching synchronous oscillations in the intracellular dynamics.

We developed an endemic model of COVID-19 to assess the impact of vaccination and immunity waning on the dynamics of the disease. Our model exhibits the phenomenon of backward bifurcation and bi-stability, where a stable disease-free equilibrium co-exists with a stable endemic equilibrium. The epidemiological implication of this is that the control reproduction number being less than unity is no longer sufficient to guarantee disease eradication. We showed that this phenomenon could be eliminated by either increasing the vaccine efficacy or by reducing the disease transmission rate (adhering to non-pharmaceutical interventions). Furthermore, we numerically investigated the impacts of vaccination and waning of both vaccine-induced immunity and post-recovery immunity on the disease dynamics. Our simulation results show that the waning of vaccine-induced immunity has more effect on the disease dynamics relative to post-recovery immunity waning, and suggests that more emphasis should be on reducing the waning of vaccine-induced immunity to eradicate COVID-19.

The pandemic of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) took the world by surprise. Following the first outbreak of COVID-19 in December 2019, several models have been developed to study and understand its transmission dynamics. Although the spread of COVID-19 is being slowed down by vaccination and other interventions, there is still a need to have a clear understanding of the evolution of the pandemic across countries, states and communities. To this end, there is a need to have a clearer picture of the initial spread of the disease in different regions. In this project, we used a simple SEIR model and a Bayesian inference framework to estimate the basic reproduction number of COVID-19 across Africa. Our estimates vary between 1.98 (Sudan) and 9.66 (Mauritius), with a median of 3.67 (90% CrI: 3.31 - 4.12). The estimates provided in this paper will help to inform COVID-19 modeling in the respective countries/regions.

Following the emergence of COVID-19 at the end of 2019, several mathematical models have been developed to study the transmission dynamics of this disease. Many of these models assume homogeneous mixing in the underlying population. However, contact rates and mixing patterns can vary dramatically among individuals depending on their age and activity level. Variation in contact rates among age groups and over time can significantly impact how well a model captures observed trends. To properly model the age-dependent dynamics of COVID-19 and understand the impacts of interventions, it is essential to consider heterogeneity arising from contact rates and mixing patterns. We developed an age-structured model that incorporates time-varying contact rates and population mixing computed from the ongoing BC Mix COVID-19 survey to study transmission dynamics of COVID-19 in British Columbia (BC), Canada. Using a Bayesian inference framework, we fit four versions of our model to weekly reported cases of COVID-19 in BC, with each version allowing different assumptions of contact rates. We show that in addition to incorporating age-specific contact rates and mixing patterns, time-dependent (weekly) contact rates are needed to adequately capture the observed transmission dynamics of COVID-19. Our approach provides a framework for explicitly including empirical contact rates in a transmission model, which removes the need to otherwise model the impact of many non-pharmaceutical interventions. Further, this approach allows projection of future cases based on clear assumptions of age-specific contact rates, as opposed to less tractable assumptions regarding transmission rates.

BackgroundWith the return of in-person classes, an understanding of COVID-19 transmission in vaccinated university campuses is essential. Given the context of high anticipated vaccination rates and other measures, there are outstanding questions of the potential impact of campus-based asymptomatic screening.MethodsWe estimated the expected number of cases and hospitalizations in one semester using rates derived for British Columbia (BC), Canada up to September 15th, 2021 and age-standardizing to a University population. To estimate the expected number of secondary cases averted due to routine tests of unvaccinated individuals in a BC post-secondary institution, we used a probabilistic model based on the incidence, vaccination effectiveness, vaccination coverage and R0. We examined multiple scenarios of vaccine coverage, screening frequency, and pre-vaccination R0.ResultsFor one 12 week semester, the expected number of cases is 67 per 50,000 for 80% vaccination coverage and 37 per 50,000 for 95% vaccination coverage. Screening of the unvaccinated population averts an expected 6-16 cases per 50,000 at 80% decreasing to 1-2 averted cases per 50,000 at 95% vaccination coverage for weekly to daily screening. Further scenarios can be explored using a web-based application.InterpretationRoutine screening of unvaccinated individuals may be of limited benefit if vaccination coverage is 80% or greater within a university setting.

IntroductionSeveral non-pharmaceutical interventions such as physical distancing, hand washing, self-isolation, and schools and business closures, were implemented in British Columbia (BC) following the first laboratory-confirmed case of coronavirus disease 2019 (COVID-19) on January 26, 2020, to minimize in-person contacts that could spread infections. The BC COVID-19 Population Mixing Patterns survey (BC-Mix) was established as a surveillance system to measure behaviour and contact patterns in BC over time to inform the timing of the easing/re-imposition of control measures. In this paper, we describe the BC-Mix survey design and the demographic characteristics of respondents.MethodsThe ongoing repeated online survey was launched in September 2020. Participants are mainly recruited through social media platforms (including Instagram, Facebook, YouTube, WhatsApp). A follow up survey is sent to participants two to four weeks after completing the baseline survey. Survey responses are weighted to BC’s population by age, sex, geography, and ethnicity to obtain generalizable estimates. Additional indices such as the material and social deprivation index, residential instability, economic dependency, and others are generated using census and location data.ResultsAs of July 26, 2021, over 61,000 baseline survey responses were received of which 41,375 were eligible for analysis. Of the eligible participants, about 60% consented to follow up and about 27% provided their personal health numbers for linkage with healthcare databases. Approximately 50% of respondents were female, 39% were 55 years or older, 65% identified as white and 50% had at least a university degree.ConclusionThe pandemic response is best informed by surveillance systems capable of timely assessment of behaviour patterns. BC-Mix survey respondents represent a large cohort of British Columbians providing near real-time information on behavioural and contact patterns in BC. Data from the BC-Mix survey would inform provincial COVID-19-related control measures.

In many micro- and macro-scale systems, collective dynamics occurs from the coupling of small spatially segregated, but dynamically active, units through a bulk diffusion field. This bulk diffusion field, which is both produced and sensed by the active units, can trigger and then synchronize oscillatory dynamics associated with the individual units. In this context, we analyze diffusion-induced synchrony for a class of cell-bulk ODE–PDE system in $${\mathbb {R}}^2$$ that has two spatially segregated dynamically active circular cells of small radius. By using strong localized perturbation theory in the limit of small cell radius, we calculate the steady-state solution and formulate the linear stability problem. For Sel’kov intracellular reaction kinetics, we analyze how the effect of bulk diffusion can trigger, via a Hopf bifurcation, either in-phase or anti-phase intracellular oscillations that would otherwise not occur for cells that are uncoupled from the bulk medium. Phase diagrams in parameter space where these oscillations occur are presented, and the theoretical results from the linear stability theory are validated from full numerical simulations of the ODE–PDE system. In addition, the two-cell case is extended to study the onset of synchronous oscillatory instabilities associated with an infinite hexagonal arrangement of small identical cells in $${\mathbb {R}}^2$$ with Sel’kov intracellular kinetics. This analysis for the hexagonal cell pattern relies on determining a new, computationally efficient, explicit formula for the regular part of a certain periodic reduced-wave Green’s function.

The spread of SARS-CoV-2 through direct transmission (person-to-person) has been the focus of most studies on the dynamics of COVID-19. The efficacy of social distancing and mask usage at reducing the risk of direct transmission of COVID-19 has been studied by many researchers. Little or no attention is given to indirect transmission of the virus through shared items, commonly touch surfaces and door handles. The impact of the persistence of SARS-CoV-2 on hard surfaces and in the environment, on the dynamics of COVID-19 remain largely unknown. Also, the current increase in the number of cases despite the strict non-pharmaceutical interventions suggests a need to study the indirect transmission of COVID-19 while incorporating testing of infected individuals as a preventive measure. Assessing the impact of indirect transmission of the virus may improve our understanding of the overall dynamics of COVID-19. We developed a novel deterministic susceptible-exposed-infected-removed-virus-death compartmental model to study the impact of indirect transmission pathway on the spread of COVID-19, the sources of infection, and prevention/control. We fitted the model to the cumulative number of confirmed cases at episode date in Toronto, Canada using a Markov Chain Monte Carlo optimization algorithm. We studied the effect of indirect transmission on the epidemic peak, peak time, epidemic final size and the effective reproduction number, based on different initial conditions and at different stages. Our findings revealed an increase in cases with indirect transmission. Our work highlights the importance of implementing additional preventive and control measures involving cleaning of surfaces, fumigation, and disinfection to lower the spread of COVID-19, especially in public areas like the grocery stores, malls and so on. We conclude that indirect transmission of SARS-CoV-2 has a significant effect on the dynamics of COVID-19, and there is need to consider this transmission route for effective mitigation, prevention and control of COVID-19 epidemic.

Extensive non-pharmaceutical and physical distancing measures are currently the primary interventions against coronavirus disease 2019 (COVID-19) worldwide. It is therefore urgent to estimate the impact such measures are having. We introduce a Bayesian epidemiological model in which a proportion of individuals are willing and able to participate in distancing, with the timing of distancing measures informed by survey data on attitudes to distancing and COVID-19. We fit our model to reported COVID-19 cases in British Columbia (BC), Canada, and five other jurisdictions, using an observation model that accounts for both underestimation and the delay between symptom onset and reporting. We estimated the impact that physical distancing (social distancing) has had on the contact rate and examined the projected impact of relaxing distancing measures. We found that, as of April 11 2020, distancing had a strong impact in BC, consistent with declines in reported cases and in hospitalization and intensive care unit numbers; individuals practising physical distancing experienced approximately 0.22 (0.11–0.34 90% CI [credible interval]) of their normal contact rate. The threshold above which prevalence was expected to grow was 0.55. We define the “contact ratio” to be the ratio of the estimated contact rate to the threshold rate at which cases are expected to grow; we estimated this contact ratio to be 0.40 (0.19–0.60) in BC. We developed an R package ‘covidseir’ to make our model available, and used it to quantify the impact of distancing in five additional jurisdictions. As of May 7, 2020, we estimated that New Zealand was well below its threshold value (contact ratio of 0.22 [0.11–0.34]), New York (0.60 [0.43–0.74]), Washington (0.84 [0.79–0.90]) and Florida (0.86 [0.76–0.96]) were progressively closer to theirs yet still below, but California (1.15 [1.07–1.23]) was above its threshold overall, with cases still rising. Accordingly, we found that BC, New Zealand, and New York may have had more room to relax distancing measures than the other jurisdictions, though this would need to be done cautiously and with total case volumes in mind. Our projections indicate that intermittent distancing measures—if sufficiently strong and robustly followed—could control COVID-19 transmission. This approach provides a useful tool for jurisdictions to monitor and assess current levels of distancing relative to their threshold, which will continue to be essential through subsequent waves of this pandemic.

IMPORTANCE: An effective vaccine against SARS-CoV-2 will reduce morbidity and mortality and allow substantial relaxation of physical distancing policies. However, the ability of a vaccine to prevent infection or disease depends critically on protecting older individuals, who are at highest risk of severe disease. OBJECTIVE: We quantitatively estimated the relative benefits of COVID-19 vaccines, in terms of preventing infection and death, with a particular focus on effectiveness in elderly people. DESIGN: We applied compartmental mathematical modelling to determine the relative effects of vaccines that block infection and onward transmission, and those that prevent severe disease. We assumed that vaccines showing high efficacy in adults would be deployed, and examined the effects of lower vaccine efficacy among the elderly population. SETTING AND PARTICIPANTS: Our mathematical model was calibrated to simulate the course of an epidemic among the entire population of British Columbia, Canada. Within our model, the population was structured by age and levels of contact. MAIN OUTCOME(S) AND MEASURE(S): We assessed the effectiveness of possible vaccines in terms of the predicted number of infections within the entire population, and deaths among people aged 65 years and over. RESULTS: In order to reduce the overall rate of infections in the population, high rates of deployment to all age groups will be critical. However, to substantially reduce mortality among people aged 65 years and over, a vaccine must directly protect a high proportion of people in that group. CONCLUSIONS AND RELEVANCE: Effective vaccines deployed to a large fraction of the population are projected to substantially reduce infection in an otherwise susceptible population. However, even if transmission were blocked highly effectively by vaccination of children and younger adults, overall mortality would not be substantially reduced unless the vaccine is also directly protective in elderly people. We strongly recommend: (i) the inclusion of people aged 65 years and over in future trials of COVID-19 vaccine candidates; (ii) careful monitoring of vaccine efficacy in older age groups following vaccination.

Following successful non-pharmaceutical interventions (NPI) aiming to control COVID-19, many jurisdictions reopened their economies and borders. As little immunity had developed in most populations, re-establishing higher contact carried substantial risks, and therefore many locations began to see resurgence in COVID-19 cases. We present a Bayesian method to estimate the leeway to reopen, or alternatively the strength of change required to re-establish COVID-19 control, in a range of jurisdictions experiencing different COVID-19 epidemics. We estimated the timing and strength of initial control measures such as widespread distancing and compared the leeway jurisdictions had to reopen immediately after NPI measures to later estimates of leeway. Finally, we quantified risks associated with reopening and the likely burden of new cases due to introductions from other jurisdictions. We found widely varying leeway to reopen. After initial NPI measures took effect, some jurisdictions had substantial leeway (e.g., Japan, New Zealand, Germany) with > 0.99 probability that contact rates were below 80% of the threshold for epidemic growth. Others had little leeway (e.g., the United Kingdom, Washington State) and some had none (e.g., Sweden, California). For most such regions, increases in contact rate of 1.5–2 fold would have had high (> 0.7) probability of exceeding past peak sizes. Most jurisdictions experienced June–August trajectories consistent with our projections of contact rate increases of 1–2-fold. Under such relaxation scenarios for some regions, we projected up to ∼ 100 additional cases if just one case were imported per week over six weeks, even between jurisdictions with comparable COVID-19 risk. We provide an R package covidseir to enable jurisdictions to estimate leeway and forecast cases under different future contact patterns. Estimates of leeway can establish a quantitative basis for decisions about reopening. We recommend a cautious approach to reopening economies and borders, coupled with strong monitoring for changes in transmission.

Extensive physical distancing measures are currently the primary intervention against coronavirus disease 2019 (COVID-19) worldwide. It is therefore urgent to estimate the impact such measures are having. We introduce a Bayesian epidemiological model in which a proportion of individuals are willing and able to participate in distancing measures, with the timing of these measures informed by survey data on attitudes to distancing and COVID-19. We fit our model to reported COVID-19 cases in British Columbia, Canada, using an observation model that accounts for both underestimation and the delay between symptom onset and reporting. We estimate the impact that physical distancing (also known as social distancing) has had on the contact rate and examine the projected impact of relaxing distancing measures. We find that distancing has had a strong impact, consistent with declines in reported cases and in hospitalization and intensive care unit numbers. We estimate that approximately 0.78 (0.66–0.89 90% CI) of contacts have been removed for individuals in British Columbia practising physical distancing and that this fraction is above the threshold of 0.45 at which prevalence is expected to grow. However, relaxing distancing measures beyond this threshold re-starts rapid exponential growth. Because the extent of underestimation is unknown, the data are consistent with a wide range in the prevalence of COVID-19 in the population; changes to testing criteria over time introduce additional uncertainty. Our projections indicate that intermittent distancing measures—if sufficiently strong and robustly followed— could control COVID-19 transmission, but that if distancing measures are relaxed too much, the epidemic curve would grow to high prevalence.