COVID-19 Epidemic Modelling For Policy Decision Support in Victoria, Australia 2020 - 2021
COVID-19 Epidemic Modelling For Policy Decision Support in Victoria, Australia 2020 - 2021
COVID-19 Epidemic Modelling For Policy Decision Support in Victoria, Australia 2020 - 2021
Abstract
Background Policy responses to COVID-19 in Victoria, Australia over 2020–2021 have been supported by evidence
generated through mathematical modelling. This study describes the design, key findings, and process for policy
translation of a series of modelling studies conducted for the Victorian Department of Health COVID-19 response
team during this period.
Methods An agent-based model, Covasim, was used to simulate the impact of policy interventions on COVID-19
outbreaks and epidemic waves. The model was continually adapted to enable scenario analysis of settings or policies
being considered at the time (e.g. elimination of community transmission versus disease control). Model scenarios
were co-designed with government, to fill evidence gaps prior to key decisions.
Results Understanding outbreak risk following incursions was critical to eliminating community COVID-19 transmis‑
sion. Analyses showed risk depended on whether the first detected case was the index case, a primary contact of
the index case, or a ‘mystery case’. There were benefits of early lockdown on first case detection and gradual easing
of restrictions to minimise resurgence risk from undetected cases. As vaccination coverage increased and the focus
shifted to controlling rather than eliminating community transmission, understanding health system demand was
critical. Analyses showed that vaccines alone could not protect health systems and need to be complemented with
other public health measures.
Conclusions Model evidence offered the greatest value when decisions needed to be made pre-emptively, or for
questions that could not be answered with empiric data and data analysis alone. Co-designing scenarios with policy-
makers ensured relevance and increased policy translation.
Keywords COVID-19, Mathematical model, Outbreak analysis, Disease control
†
Nick Scott and Romesh G. Abeysuriya are equal first authors.
*Correspondence:
Nick Scott
Nick.Scott@burnet.edu.au
Full list of author information is available at the end of the article
© The Author(s) 2023. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which
permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the
original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or
other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line
to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory
regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this
licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecom‑
mons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
Scott et al. BMC Public Health (2023) 23:988 Page 2 of 12
Prospective outbreak • Probability of outbreak growing under No community cases; elimination No parameter fitting, sensitivity analyses
different theoretical circumstances strategy for key parameters
Reactive outbreak • Probability of existing outbreak reaching Small number of diagnosed cases fol‑ No parameter fitting, filter outbreaks to
different levels or being eliminated lowing period of no community cases; match observed
• Time required to contain outbreak elimination strategy
Easing restrictions • Timing and magnitude of epidemic Large number of cases; elimination Fit model to cases/hospital demand/
peak strategy deaths
• Probability of resurgence following eas‑
ing of restrictions
Managing health • Timing and magnitude of epidemic Large number of cases; epidemic control
system utilization peak(s) strategy
• Peak hospital/ICU demand
• Number of deaths
Scott et al. BMC Public Health (2023) 23:988 Page 4 of 12
As with the prospective outbreak analyses, the main impact of dynamically introducing or easing restrictions
output measure was the percentage of simulations reach- in the context of higher case numbers. The key outcome
ing different categories of 7-day average daily detected measures were hospital and ICU demand rather than the
cases and key questions focused on the duration of number of cases.
restrictions required to contain the outbreak.
Results
Scenario type 3: easing restrictions (elimination strategy An overview of overarching strategy and when principal
context) scenario types were used is shown in Fig. 1, overlaid on
For situations where a large outbreak had occurred and case numbers and vaccine coverage for epidemic context.
restrictions had already been imposed, policy questions
arose about when to ease restrictions and whether to do Easing restrictions example: wild‑type epidemic wave,
it incrementally or collectively. Model scenarios were leading to elimination (September 2020)
calibrated analogously to the constrained outbreak analy- Following a wild-type epidemic wave in Victoria (June-
ses, where the model was initialized with a small num- Nov 2020) that was contained by imposing restrictions,
ber of cases and simulations were only retained if they an analysis was conducted on 14 September to assess
were consistent with the actual outbreak, accounting for the impact of the timing of easing restrictions on resur-
restrictions imposed to date. Using only the retained sim- gence risk (Fig. 2) [8]. Model calibration involved running
ulations, scenarios were run comparing outcomes with simulations starting from zero cases, sampling over ini-
restrictions in place for different periods. tial seed infections (i.e. who in the model is the first case)
The main outcome measure was the probability of and uncertainty ranges of calibrated transmission param-
reaching > N diagnoses per day following the easing of eters, and retaining simulations within a threshold of the
restrictions (i.e., “resurgence risk”), and how this varied case data (Fig. 2, grey lines). For a set of 1000 retained
according to the timing and extent that restrictions were simulations, scenarios were applied to a short-term pro-
eased. jection to compare the impact of easing restrictions
on different dates. The red and blue lines in Fig. 2 show
Scenario type 4: health system utilization (control strategy individual simulations for each restriction scenario, high-
context) lighting the broad range of outcomes possible. For each
As the broader COVID-19 strategy transitioned from scenario, a proportion of simulations trend to zero while
elimination to control, health system utilization became a proportion result in a resurgence of infections, defined
increasingly relevant. Scenarios were calibrated similarly as reaching a threshold of daily diagnoses. This analysis
to the easing restrictions analyses, and compared the showed that an additional two weeks of restrictions could
Fig. 1 Victorian COVID-19 timeline. New daily cases and vaccine coverage are shown over time, contextualising the overall strategic phase
(elimination or suppression) and the dominant scenario type being run at different stages of the pandemic (vertical shaded regions)
Scott et al. BMC Public Health (2023) 23:988 Page 5 of 12
Fig. 2 Wild-type epidemic wave in Victoria in 2020, in the context of an elimination COVID-19 strategy, with scenarios considering resurgence risk
if restrictions were eased on 14 Sep (red) or 28 Sep (blue). Model simulations were started with random infected seed cases and randomly sampled
transmission parameters, and were retained if they were within sufficient bounds of the observed data. Throughout the simulations ‘Stage 3’, masks
and ‘Stage 4’ restrictions were imposed, with their impact derived through model calibration (see [8])
more than halve the risk of a resurgence in cases. Factor- days of the start of the simulation (e.g. >30 in Fig. 3). All
ing this and other evidence, restrictions were maintained simulations included contact tracing (including second-
and cases continued to decline, before being gradually ring tracing for some contacts).
eased from 28 September (e.g. small gatherings of five Several factors were identified that influenced out-
people outdoors). break risk including the infectiousness of the variant,
the number of generations of transmissions before the
Prospective outbreak analysis example (October 2020) first detected case, and the policy response to the ini-
After the incidence of community acquired cases tial case detection. The number of infection generations
returned to zero from 29 to 2020, outbreak risks from prior to the first diagnosis being recorded was particu-
potential incursions were assessed. Simulations started larly predictive of the outcome. In simulations where the
with one undiagnosed case. The outbreak risk was index case was the first identified, light restriction con-
defined as the percentage of simulations that reached a tained the outbreak 94% of the time (Fig. 3). However, if
particular 7-day average diagnosis threshold within 90 the index case was not the first case detected, this made
Fig. 3 Outbreak analysis in Victoria in 2021, in the context of pursuing an elimination strategy. Simulations were run starting from zero cases and a
randomly infected seed case. Table shows the proportion of simulations where an outbreak occurred (defined as reaching a 7-day average of > 30
diagnoses/day within 90 days), according to which infection generation was detected, the infectiousness of the virus (baseline = wild type; 50%
more infectious ~ Alpha variant) and what responses were implemented on detection of the first case (light = masks and limits on social gatherings;
moderate = light + density limits; heavy = lockdown)
Scott et al. BMC Public Health (2023) 23:988 Page 6 of 12
containment much more difficult. At the time of analysis, the first few diagnoses. In May 2021 an outbreak of the
the Alpha variant was estimated to be 50% more infec- Kappa variant occurred, with the first diagnosis on 24
tious than the wild type variant [10, 11], and in simula- May, and 45 and 58 cases were diagnosed cumulatively
tions where an Alpha variant infection was identified in within the next 7 and 10 days, respectively. Mandatory
the community and the index case was not the first case masks were imposed on the first day, and a lockdown was
diagnosed, then even a moderate response (mask man- imposed on the third day.
dates, limits on social gatherings, and density limits) Modelling was used to estimate the duration of the
would be unlikely to contain it (Fig. 3). These analyses lockdown needed to contain the outbreak. Simulations
highlighted the importance of interventions to increase started with a single case and were retained if transmis-
case detection (e.g., asymptomatic screening of quaran- sion occurred (i.e., the incursion did not fizzle out), the
tine hotel workers). index case was never diagnosed, and there were between
The impact of the timing of the response to an outbreak 43 and 65 diagnoses after 7 days and 53–80 diagnoses
was also assessed. For early variants (wild-type, Alpha), after 12 days. Using 1000 simulations that met these
heavy restrictions were always able to contain outbreaks, criteria (sampling over initial seeds, contact network
but a delay in introducing restrictions increased the dura- structures and transmission parameters), scenarios were
tion of restrictions to contain the outbreak (Fig. 4). The run for policy changes under consideration, where after
analysis showed that if some restrictions were already two weeks of lockdown (June 11), either the lockdown
active at the time of the incursion, the duration of restric- was maintained, or restrictions were eased by: opening
tions would also be shorter. schools only; opening schools and venues with density
This supported a “go hard go early” approach to restric- limit; opening schools, venues with density limits and
tions in the elimination context, which was applied when allowing small social gatherings; or returning to masks
lockdowns were imposed from 13 to 17 February 2021 only. The impact of these policy changes on daily diag-
after 13 Alpha variant cases were detected, from 28 May noses four weeks after implementation was recorded as
to 3 June 2021 (later extended to 10 June) after 30 Kappa the main outcome measure (Fig. 5). The analysis quanti-
variant cases were detected, and from 16 July to 20 July fied how much the risk of a resurgence increased when
(later extended to 27 July) after 18 Delta variant cases restrictions were eased more generously, and enabled
were detected. After each lockdown restrictions were resurgence risk to be incorporated into decision-making
incrementally eased, rather than lifted at once. alongside the costs associated with the lockdown.
Based on this and other evidence, from 11 June schools
Reactive outbreak analysis example: Kappa variant were re-opened, hospitality was reopened with den-
outbreak (May 2021) sity limits, small gatherings were allowed but only out-
Once an incursion was detected, outbreak risk esti- doors, and work from home and mask mandates were
mates were refined by incorporating characterisations of maintained.
Fig. 4 Outbreak time to containment. Following an outbreak, the time required in lockdown to return to < 5 diagnoses per day, according to the
delay in implementing restrictions. Lines show median and inter-quartile range (error bars) across 1000 simulations. Based on the wild type variant
Scott et al. BMC Public Health (2023) 23:988 Page 7 of 12
Fig. 5 Outbreak analysis in Victoria in 2021, when an outbreak has been detected but limited information is known, pursuing an elimination
strategy. Top: Simulations were run starting from zero cases and a randomly infected seed case. Simulations were only retained if they produced,
through stochastic variation, approximately the same initial outbreaks as occurred. In this example it was 43–65 diagnoses within 7 days of the first
diagnosis, and 53–80 diagnoses within 12 days of the first diagnosis. This is against a background of masks and work from home if possible being
implemented on day 1, and lockdown being implemented on day 3. Bottom: From the retained simulations, scenarios compared whether the
lockdown was maintained, or after 14 days either schools were reopened, venues were opened with density limits, small social gatherings were also
allowed, or a return to masks and working from home only
Prospective outbreak analysis example: context of vaccine from home) could potentially reduce COVID-19-related
rollout and coverage (March 2021) mortality to a similar rate as experienced in the 2017
When vaccination coverage increased, other public influenza season (a year with particularly high influenza-
health responses were also expected to change. In the related mortality), depending on vaccine efficacy. In
context of vaccines becoming available, outbreak risk addition, a trade-off was identified between the degree
(defined as the percentage of simulations reaching dif- of restrictions and the duration, where similar case num-
ferent diagnosis thresholds over the first 90 days after a bers could be achieved with a short duration of strict
single case) was assessed by vaccine type and coverage restrictions, or longer duration of light restrictions. This
(Fig. 6). While outbreak risk reduced with increasing vac- showed that even with high vaccine coverage, without
cine coverage, it was clear that herd immunity through NPIs the health impacts of COVID-19 would be much
vaccination was unlikely to be achieved, supporting an worse than influenza, but that intermittent light restric-
eventual move from an elimination strategy to a control tions could considerably improve outcomes.
strategy.
Managing health system utilization example: Delta variant
Managing health system utilization example: moving epidemic wave (Sept 2021)
from elimination to control (July 2021) Prior to the Delta variant, NPIs were sufficient to contain
Later we investigated the feasibility of a range of possi- outbreaks in unvaccinated populations. However, in mid-
ble control strategies, from intermittent restrictions to 2021 Victoria experienced a large Delta variant epidemic
removal of all non-pharmaceutical interventions (NPIs) wave before high coverage of vaccination was achieved
[12, 13]. Since vaccines were highly efficacious but imper- (as of 5 August only 43% of Victorians over 16 years had
fect, removing all NPIs was found to come at consider- received one dose and 21% two doses). With the Delta
able health cost (Fig. 7). This analysis showed that with variant having a shorter serial interval that reduced the
the Delta variant, intermittent low-level restrictions (for effectiveness of contact tracing, and community fatigue
example, high levels of testing, mask mandates and work leading to declining compliance with restrictions,
Scott et al. BMC Public Health (2023) 23:988 Page 8 of 12
Fig. 6 Outbreak analysis in Victoria in 2021, over 90 days starting from a single case and different levels of population vaccine coverage. Left:
AstraZeneca type vaccine, with assumed protection against infection of 50%, and protection against symptomatic disease of 65%. Right: Pfizer type
vaccines, with assumed protection against infection of 80%, and protection against symptomatic disease of 90%. TPOD = total protection against
disease; POI = protection against infection. Note that vaccine efficacy estimates were based on best available data as at March 2021
community transmission continued to increase despite a phase, community incursions from hotel quarantine reg-
lockdown being imposed. Plans to transition away from ularly occurred, and modelling was used as a critical tool
an elimination approach were accelerated, with a focus to understand outbreak risks and how they varied in dif-
on using NPIs to control the epidemic while conducting a ferent circumstances or with different responses.
mass vaccination campaign. Modelling was used to assess In the elimination phase, model outputs showed out-
the implications of different strategies (i.e., sequences break risk depended on whether the detected case was
and timings of restrictions being eased, relative to the the incursion, first ring of infections or a ‘mystery case’.
vaccination rollout and case numbers) on infections and There were benefits of early lockdown on detection of the
health system demand [14, 15]. This analysis showed that first case and exiting lockdowns slowly through gradual
maintaining testing among those vaccinated was likely to easing of restrictions to minimise resurgence risk from
be necessary for limiting transmission, and early actions undetected cases. Key metrics in this period were the
to achieve transmission reduction could have significant number of cases, the intensity of restrictions, and the
downstream implications due to rising vaccine cover- duration of restrictions, as returning to COVID-zero was
age (Fig. 8). The roadmap that was developed ultimately viewed as the most efficient way to meet other targets.
included incremental easing of restrictions in line with The emergence of COVID-19 vaccines meant that
different vaccination coverage thresholds. catastrophic health outcomes could be avoided with-
out the need for some of the stricter border control
Discussion and lockdown measures. This, combined with the fact
After Australia’s initial response successfully eliminated that the elimination of community transmission was no
community COVID-19 transmission from most jurisdic- longer realistic due to the emergence of more infectious
tions in early 2020, and as the dire health consequences of variants, led to a move from an elimination to a control
widespread transmission in the absence of a vaccine were strategy. While high vaccine coverage was an essential
observed internationally, Australia embarked on a phase element, our analyses indicated that sustainable con-
lasting until mid-2021 where similar catastrophic health trol strategies would require vaccines plus intermittent
outcomes were avoided by using tight border controls, ‘light’ restrictions, introduced early in times of case
contact tracing and restrictions to achieve and maintain number escalation, to avoid excess case and health sys-
the elimination of community transmission. During this tem demand escalation. During the critical transition
Scott et al. BMC Public Health (2023) 23:988 Page 9 of 12
Fig. 7 Projected COVID-19 deaths in Victoria over a 12-month period under different outbreak response scenarios, triggered by case numbers
(left), and average time spent under restrictions (right). Top: Grey bars show the projected number of deaths on a log scale, with the lighter shading
indicating the proportion that are among the vaccinated population. Scenarios are based on 95%/70%/70% vaccine coverage being achieved
among people aged 60+/12–59/<12 years, the vaccine having an assumed 80% protection against infection and 92% protection against death,
and 5 cases per day were seeded into the community through reduced quarantine measures. Light restrictions = mandatory masks, density limits
and work from home if possible; medium restriction = light restrictions + additional limits on gathering sizes; lockdown = mandatory masks, work
from home, schools closed, retail closed, hospitality take-away only, social gatherings up to two outdoors only. Bottom: example time series output
from one simulation, where increasing restrictions are triggered with increased hospital number thresholds, to maintain disease control
period from an elimination to a control strategy, and answered directly with data, this meant that the ration-
as vaccine coverage increased, the focus of analyses ale for the use of modelling was always clear.
shifted to health system demand, with case numbers The utility of modelling was also increased by the co-
becoming less important. design of scenarios with government and other relevant
In all analyses conducted, model evidence was most stakeholders [16]. Scenarios were only modelled if (a)
crucial for decision making when used to investigate they could not be shown to be superior / inferior based
questions that could not be answered in any other way. on existing data; (b) they were considered reasonable by
For example, where relevant data was available, direct epidemiologists and public health teams; and (c) they
analysis was significantly faster to inform decisions and were being seriously considered by government given
required fewer assumptions to be made. Model analy- other economic, social, and political implications. Many
ses were therefore typically used to support decisions policy options could be ruled out before even making it
that needed to be made pre-emptively (e.g., how long to to the modelling stage, because they were not likely to
lockdown for, when minimal information was available) be feasible for other reasons. For example, working with
[16]. Since the primary aims of each study could not be stakeholders ensured that scenarios could incorporate
Scott et al. BMC Public Health (2023) 23:988 Page 10 of 12
Fig. 8 2021 roadmap with reduced testing for people vaccinated (blue), maintained testing for people vaccinated (red), and an additional 15%
reduction in non-household transmission (green). Scenarios include schools returning to in person learning throughout October; childcare
returning and mobility restrictions easing in October; limited outdoor gatherings at 70% two-dose vaccine coverage among people 16 + years;
indoor gathering with density limits at 80% two-dose coverage among people 16 + years and mandatory vaccine requirements. Dashed vertical
lines represent estimated dates of reaching 70% and 80% two-dose coverage among people 16 + years
operational or capacity constraints in the delivery of ser- evidence. However, in the case of an outbreak, decisions
vices such as testing or vaccination. Co-designing the often need to be made quickly, with significant impacts on
scenarios also helped align model outputs with other the lives of large numbers of people. It is therefore imprac-
evidence under consideration, to make the modelling as tical to have outputs peer-reviewed in advance, which can
informative as possible. For example, running epidemic take months. Review and validation of raw model outputs
projections for similar scenarios to those separately by government analysts on standby was conducted in place
used for economic forecasts facilitates integrating both of such checks, and where feasible identical requests were
sources of evidence. In addition, the continued involve- often put to other modelling groups [17, 18] to be con-
ment of the modelling team with public health teams ducted simultaneously and independently [19, 20].
allowed ongoing refinement and improvement of the Detailing how these analyses were conducted is impor-
modelling, and validation against past accuracy. tant because the role of modelling in informing policy is
There were many challenges associated with the need not well understood. Epidemiological modelling is a spe-
for fast decision-making. Ideally, modelling analyses would cialized discipline that experienced increased media atten-
undergo peer review before being used to inform decisions. tion because of its utility to generate evidence to support
The purpose of the peer-review process is to ensure that decision-making in the COVID-19 pandemic; however, it
model design, assumptions and inputs are fit-for-purpose, is often not communicated well through the media [19].
consistent with best practice, and incorporate the latest There are common misconceptions about the differences
Scott et al. BMC Public Health (2023) 23:988 Page 11 of 12
between forecasting (trying to predict the specific course must be used in conjunction with other information
of the epidemic) and scenario analyses (comparing projec- sources to make policy decisions.
tions under different policy options to estimate the impact
of policies and inform decisions) [21]. This often results in Supplementary Information
a dichotomisation of modelling being viewed by the public The online version contains supplementary material available at https://doi.
as being right or wrong, with many of the key insights gen- org/10.1186/s12889-023-15936-w.