Research Article
Print
Research Article
Border biosecurity interceptions for air passengers – assessing intervention methods and analytic tools
expand article infoNicholas P. Moran, Anca M. Hanea, Andrew P. Robinson
‡ The University of Melbourne, Parkville, Australia
Open Access

Abstract

At-border interventions are a critical step along the biosecurity continuum, to measure and control the risks associated with the cross-border movement of people and goods. Air passengers are a high-volume pathway for a range of biosecurity risk materials, against which various interventions may be used (e.g., manual searches, detector dogs, x-rays, etc.). Tasmania is an island state, and its environment and industries have benefited from a low level of invasive pests due to their geographic isolation. Therefore, relatively strict at-border interventions are used to prevent the entry of new pests, including some serious invasive pests already present on mainland Australia (e.g., Queensland (Bactrocera tryoni) and Mediterranean (Ceratitis capitata) fruit fly). Using a large interception database for domestic air passengers entering the southern Australian state of Tasmania from mainland Australia, this study applies common statistical modelling tools to assess the efficacy of interventions (namely, dog detectors, and bag searches), and to identify pathway risk factors (e.g., flight origin/route). This analysis considered the effects of interventions on both voluntary declarations by passengers and also detections of undeclared risk material on passengers. The analysis also focused on biosecurity risk items generally (e.g., fruits and vegetables, meat products, cut flowers), and items that are specifically considered to be fruit-fly hosts. The results highlight that active at-border interventions and the presence of biosecurity inspectors capture a significant volume of biosecurity risk items at the border, and detector dogs have particularly strong positive effects on the rate of interceptions, particularly for items detected on passengers. Conducting bag searches also appears to increase interceptions, both by increasing the rate of items being detected and by encouraging voluntary declarations. Sensitivity analyses then test the robustness of results to modelling implementation methods and distributional assumptions. This study demonstrates how statistical modelling can provide robust insights into biosecurity interventions and risk factors along pathways, and further highlights the value of high-quality interception data resources for informing and improving biosecurity systems.

Key words:

Bactrocera tryoni, border biosecurity, Ceratitis capitata, detector dogs, fruit fly, invasive species, passenger screening, pathway risk analysis

Introduction

Biosecurity border interventions seek to balance the need for the cross-border movement of goods and people against their biosecurity risks. Although intervention policies inherit from a common framework of international agreements (Outhwaite 2010), there are substantial differences in the policies and implementation of biosecurity interventions across jurisdictions. These differences produce variation in the risk of pest introductions across states, countries and regions (Whattam et al. 2014; Epanchin-Niell et al. 2021). Australia’s biosecurity system tends to have a higher appropriate level of protection and to implement comparatively stricter interventions than many other jurisdictions (Whattam et al. 2014; Black and Bartlett 2020). This is partially due to the opportunity presented by Australia’s historical isolation, which has made its primary industries relatively pest-free on a global scale. However, activities such as tourism and trade have reduced and will continue to reduce this isolation, increasing the risk of pest introductions (Turner et al. 2021; Whattam et al. 2024). Recent studies estimate an aggregated cost of invasive species to Australia to be AU$389.59 billion since the 1960s (Bradshaw et al. 2021), and a net present value of AU$314 billion for the national biosecurity system in terms of the assets that it protects (Dodd et al. 2020; Stoeckl et al. 2023). This highlights the critical importance of ensuring that border interventions act as effective barriers against biosecurity threats.

This is particularly true for the Australian island state of Tasmania (Fig. 1). Due primarily to their geographic and evolutionary isolation, island ecosystems possess a disproportionate level of the earth’s endemic species and support a large percentage of its biodiversity (Kier et al. 2009; Weigelt and Kreft 2013). Tasmania has specifically been identified as a major centre of endemism for Australian flora; for example, more than half of the 30 native Eucalyptus species in Tasmania are endemic (Crisp et al. 2001; Potts et al. 2016). Local industries, communities, and natural ecosystems benefit from the state’s relative isolation and low levels of pests, including species that are present elsewhere in Australia, such as Queensland fruit fly (‘Qfly’, Bactrocera tryoni; Florec et al. 2013), Mediterranean fruit fly (‘Medfly’, Ceratitis capitata; Cook and Fraser 2015), tomato potato psyllid (Bactericera cockerelli; Moir et al. 2022), and grape phylloxera (Daktulsphaira vitifoliae; Skinner 2018). These are potential threats for Tasmania, and a 2018 incursion of Qfly in the state’s north cost millions in direct eradication costs in addition to further indirect costs (e.g., via temporary market access losses; Blake 2019). Tasmania’s low-pest status is therefore a biosecurity challenge, and islands can be particularly vulnerable to impacts from invasive pests and diseases (Keitt et al. 2011; Fraser 2016; Brettell et al. 2021). Furthermore, globalisation is expected to increase movement across borders in both goods and people and to continue to bridge the geographic barriers that once kept Australia and Tasmania isolated (Dodd et al. 2015; Seebens et al. 2017, 2021).

Figure 1.

Location of Tasmania in relation to mainland Australia (inset), including the locations of the seven arrival ports for air passengers, namely (from left to right) 1. King Island, 2. Burnie, 3. Devonport, 4. Launceston, 5. Bridport, 6. Hobart, and 7. Flinders Island. Base map produced via QGIS (v3.24.2; www.qgis.org). (Note, airports have been anonymised for the remainder of the analysis and labelled Airport_A, Airport_B. etc.)

Air passengers are an important high-volume pathway for pest introductions, with around 20 million passengers arriving annually in Australia in 2023–24 (BITRE 2024). Pests may be introduced via luggage (e.g., via infested fruit) or by being attached to passenger clothing or belongings (e.g., soil on shoes or sports equipment; McNeill et al. 2011; Pace et al. 2022; Robinson and McNeill 2022). Air passenger pathways may be particularly important for pest insects, including Medfly (Liebhold et al. 2006; McCullough et al. 2006). Passengers are also a pathway for animal diseases, for example, measures targeting African swine fever detected a yearly average of 33,684 pork products from 2% of screened passengers entering Australia from 2021/22–2023/24 (DAFF 2024). Interventions on these pathways must therefore be able to mitigate a high volume and a diverse range of biosecurity threats.

A range of risk mitigation tools can be employed at multiple points on the air passenger pathway, from pre-departure, in-transit and on-arrival screening/inspection phases of the biosecurity continuum (Sequeira and Griffin 2014; Whattam et al. 2014). For international arrivals into Australia, travellers are subject to pre-arrival risk profiling, and a subset of arrivals are subject to active interventions/screening at the border (Inspector-General of Biosecurity 2019). Common at-border interventions include manual examination, dog detector teams, and x-rays (Inspector-General of Biosecurity 2022). Tasmania also applies similar active interventions for domestic arrivals, including manual bag searches and detector dogs.

Detector dogs may be particularly valuable in air passenger screening, being able to screen large volumes of passengers and luggage efficiently, and able to be trained to target general biosecurity materials as well as specific pests and diseases (Whattam et al. 2014; Moser et al. 2020). Nonetheless, there are limited studies assessing their efficacy relative to other intervention methods. Furthermore, biosecurity interventions are often targeted towards specific flights to maximise the utility of limited resources, e.g., flight-based-traveller profiles used to target international arrivals into Australia (Inspector-General of Biosecurity 2019). Therefore, further empirical evidence about the relative efficacy of different interventions may help target intervention resources towards the highest-risk arrivals.

The analysis of biosecurity interception/surveillance data is prone to some common issues in statistical modelling, including zero-inflation (i.e., where data includes a large proportion of zeros, for example where detections of targeted items are rare), overdispersion (e.g., where variance is much higher than predicted), and censoring (e.g., if data is only recorded where contamination is detected; Kachigunda 2020; Turner et al. 2020; Trouvé and Robinson 2021; Kachigunda et al. 2022). Failing to account for overdispersion or zero-inflation in data can lead to biased or inaccurate parameter or error estimates (Harrison 2014; Campbell 2021; Feng 2021). Although some studies suggest that the outputs of mixed-effects models can be robust to violations of distributional assumptions (e.g., Schielzeth et al. 2020; Knief and Forstmeier 2021), exploring the potential effects of model design and implementation may be important considerations when using interception data to inform biosecurity decision-making.

Focusing on domestic flight arrivals into Tasmania from mainland Australia, the goal of this study is to assess the efficacy of passenger interventions and pathway risk factors on biosecurity interceptions. This focuses on both general biosecurity risk material (‘BRM’) interceptions, and interceptions of BRM material specifically relevant to Qfly and Medfly (collectively referred to as fruit fly, ‘FF’). Preventing FF incursions has been a focus of Tasmania’s border biosecurity system, particularly following the 2018 incursion. This study uses recent air pathway intervention data for Tasmania. These data are rich resources for our study, as they include records for all commercial arrivals, as well as relevant data on the types/amount of BRM intercepted. This provides a valuable opportunity to apply statistical modelling approaches and assess their sensitivity to implementation methods. The specific aims of this analysis were:

  1. To determine the relative effects of different interventions (namely, dog detector teams and luggage searches) on the rates of BRM and FF host interceptions, including voluntary declarations by passengers and involuntary detections of items by biosecurity officers. We did not make any specific directional predictions about the effects of searches and detector dogs on interceptions.
  2. To identify pathway-risk heterogeneity based on the origin and specific routes of flights. We expected substantial variation in interception rates related to flight origin and route, which may be used to identify high-risk arrivals.
  3. To test whether our results are sensitive to overdispersion and zero-inflation by implementing Bayesian mixed models with zero-inflated Poisson and negative binomial distributions. We expected the outputs of models and the estimated effects of intervention methods to be robust to different implementation approaches.

Methods

Data context and overview

Tasmania is an island state (see Fig. 1), with a cool temperate climate, unique natural ecosystems characterised by high endemism (Crisp et al. 2001; Potts et al. 2016), and a large primary industry sector with an income from agriculture, forestry and fishing industries worth around AU$3.5 billion in 2022–23, or ~9% of the Gross State Product (ABS 2023).

The main entry pathway for domestic passengers is from flights originating in six Australian mainland states/territories, with a significant but smaller volume of maritime arrivals (e.g., ferries, cruise vessels, private vessels, etc.). Interceptions from air passengers are recorded in the Biosecurity Activity Database System (referred to as ‘BAS data’), from which data from 1 January 2019–1 September 2023 were available. There are BAS interception records for 59,917 domestic interstate flight arrivals, carrying over 6.5 M passengers (~1.4 M/year on average), from which 66,675 BRM interceptions were made. This study focuses specifically on domestic interstate arrivals, as these are routinely recorded in BAS data, and domestic flights represent a large majority of air arrivals in the state (~99% based on 2023/24 volumes; BITRE 2024).

BRM items are generally defined to include fresh produce (i.e., fruits and vegetables), animal products including seafood, live animals, plant material (e.g., nursery stock, seeds), and soil attached to sports equipment or clothing (DNRET 2023). Biosecurity interventions for Tasmania have a particular focus on preventing incursions of FF into Tasmania, and a large subset of BRM interceptions (43,803, or approximately 2/3) are of items considered to be FF hosts. For this analysis, FF hosts include 130 taxa listed as Medfly and/or Qfly hosts in the Plant Biosecurity Manual Tasmania 2023 (Biosecurity Tasmania 2023). For details of BRM and FF host item definitions and of intercepted BRM items, see Suppl. material 1: A.

Data processing

Data for 59,917 flight arrivals was found to be within the scope of analysis, which excludes flights from international origins or within the state, and flights with no data (i.e., cancelled, diverted, missed, or cleared remotely; ~14% of all records). A further subset of 27 arrivals was excluded because of apparent data entry issues, and 25 arrivals into one airport were excluded as no commercial flights arrive at this location, and interceptions for the remaining private arrivals were extremely low, causing computational issues with model implementation. As only a small fraction of actual arrivals is excluded, and missed arrivals do not appear to be targeted/biased towards specific arrivals, we are confident that these exclusions do not reduce the operational relevance of the analyses.

Six count variables were used as response variables, namely:

  1. the total number of BRM interceptions per flight (N_Total);
  2. the number of BRM declarations by passengers (N_Declarations);
  3. the number of undeclared BRM interceptions (N_Detections);
  4. the total number of FF host interceptions per flight (N_Total_FF);
  5. the number of FF host declarations by passengers (N_Declarations_FF); and,
  6. the number of undeclared FF host interceptions (N_Detections_FF).

Total BRM and FF host interceptions are the sum of their corresponding declared and undeclared detection counts. Both BRM and FF variables were used to explore how interventions perform against both general biosecurity threats as well as high-priority/high-risk biosecurity materials, respectively. Detections may occur through manual searches or via detector dogs, and biosecurity staff also ask for passengers to voluntarily declare any BRM items, both of which may occur at several stages of the arrival process (e.g., as passengers enter terminals, or in luggage collection areas). Both detections and declarations were included to explore how interventions influence both voluntary and involuntary compliance behaviour in passengers (e.g., whether detector dogs primarily increase interception through direct detections, or whether their visual presence also encourages voluntary declarations).

The number of interceptions was calculated as the sum of each distinct type of BRM or FF host, separated by the passenger (e.g., if 2 passengers are intercepted each carrying 3 types of BRM, N_Total = 6). The rationale is that each commodity type may represent a distinct biosecurity threat, as may the same kind of commodity being carried by two separate passengers.

Statistical analysis A: Intervention and pathway risk effects

Generalized linear mixed effects (‘glm’) models with a Poisson distribution were implemented via package ‘lme4’ (v1.1-33, Bates et al. 2015), in the R statistical environment (v4.2.3, R Core Team 2013). This was chosen for the primary analysis, as lme4 is an accessible package that can implement models using common distributions, relative to more complex Bayesian implementation methods that may be required for more advanced model types. Therefore, this approach may be more relevant for use by non-academic users such as biosecurity managers.

Four fixed effects were included, namely: arrival airport ("Location"), intervention regime ("Regime"), number of bag searches ("BagSearchCount"), and number of passengers per flight ("PassengerCount"). "Regime" includes five combinations of one or more dog detector team (‘DDT’) and biosecurity inspector (‘BI’), i.e. one BI, two BIs, one DDT, one DDT with one BI, and two DDTs. Both DDTs and BIs have been deployed across all airports. Airports were included as fixed effects, as they may differ both in their interception efficacy and in the underlying rates of contamination on flights arriving at each location. Count predictor variables (i.e., passenger and bag search counts) were square-root transformed and Z-scaled to reduce skewness, to improve both model performance and the interpretability of effect estimates (per Schielzeth 2010).

Models included two random effects, to assess the level of variance associated with the flight’s Australian state/territory of origin ("FlightOrigin") and specific flight route ("FlightNumber"; nested within origin). Flights without a number recorded were categorised as ‘Itinerant/Other’, with a large majority considered to be private non-commercial arrivals, but also likely to include a small percentage of commercial flights for which their numbers were not entered into the database. For further details of model structure see Suppl. material 1: B.

Unless otherwise stated, all values in square brackets below represent 95% confidence intervals (or credibility intervals for Bayesian models below; ‘95CI’) for the estimated effects. The statistical significance of any fixed effects is inferred from whether their 95CIs include zero. Random effects are assessed based on how much variance is explained in models, and whether 95CIs for any specific random intercept predictions include zero. Where appropriate, parameter/effect estimates below have been converted to percentage changes in the expected number of interceptions for ease of interpretability. Marginal means were extracted from models using the package ‘emmeans’ (v1.8.7, Lenth 2023), to estimate expected interception rates under different intervention regimes.

Statistical analysis B: Model sensitivity

Sensitivity to overdispersion and zero-inflation was tested by re-fitting a subset of models in a Bayesian framework via the package ‘brms’ (v2.19.0, Bürkner 2017). From the six response variables used in the main analysis, two were selected for sensitivity analyses. These were the total BRM interceptions (N_Total, i.e., the most inclusive aggregation of interception data), and the number of FF host declarations (N_Declarations_FF, i.e., the most sparse response variable).

Four alternative distributions were tested for each response variable, namely Poisson (as in the main analysis but implemented in a Bayesian framework), zero-inflated Poisson, negative binomial, and zero-inflated negative binomial, respectively. These were chosen as common alternatives to account for cases with excess zeros and overdispersion in ecology and other fields where count data is common (Lindén and Mäntyniemi 2011; Campbell 2021; Pittman et al. 2022). Models used the same fixed and random effects specifications as in the glm models, with default noninformative priors to reflect our lack of prior knowledge for parameter estimates (chains = 3, iterations = 3000, warmup = 1000). The outputs for the fixed effects of intervention regimes and bag searches and the random intercepts associated with flight origins were estimated and qualitatively compared between models.

Measures of model fit were also estimated for all models, i.e., Akaike/Watanabe–Akaike information criterion ‘AIC’/‘WAIC’ as a measure of the quality of model fit for the dataset. Marginal and conditional R2 values as measures of the proportion of variance were explained by fixed effects and both fixed and random effects respectively (via package ‘performance’, v0.10.3, Nakagawa and Schielzeth 2013; Lüdecke et al. 2021). Overdispersion and zero-inflation tests were also conducted (also via ‘performance’).

Results

Intervention and pathway risk effects

Models identified significant effects of biosecurity interventions upon interceptions of both BRM and FF host items. Estimated BRM and FF host interception rates were significantly higher when detector dogs were present. For example, the estimated total BRM interceptions (N_Total) per flight with one DDT was 0.88 [95CI: 0.80, 0.96], compared to 0.50 [0.46, 0.55] for one BI when using our default Poisson model. This effect appeared to primarily be driven by increases in the number of detections, with DDTs having strong positive effects on BRM and FF host item detections (Fig. 2). Full results, code, models and outputs are available via Open Science Framework (osf.io/78tv9/; doi: 10.17605/OSF.IO/78TV9), and detailed model outputs are available in the Suppl. material 1: B.

Figure 2.

Estimated interception rates for air passengers under different border intervention regimes, for (A) biosecurity risk material (BRM) interceptions and (B) the subset of BRM that are fruit fly (FF) host items. Regimes include combinations of detector dog teams (DDT) and biosecurity inspectors (BI). Note, that estimated rates are the predicted number of interceptions per flight, and are independent of other factors included in the models (i.e., are estimated based on a flight with a mean number of passengers, and a mean number of bag searches, and averaged across arrival airports). Total, declared and detected rates come from distinct models, so estimates are not expected to be additive.

The number of bag searches conducted (performed both by DDTs and BIs) had a positive effect on the estimated total, declared, and detected BRM counts. Unsurprisingly, the square root number of bag searches per arrival was associated with a per-unit increase in total BRM interceptions of 24.0% [23.3%, 24.6%]. A similar effect was found for BRM declarations (i.e., 23.8% [22.8%, 24.9%]) and detections (i.e., 24.1% [23.3%, 24.9%]). Similar effects were observed for FF host interceptions in total (23.8% [23.0%, 24.6%], declarations (23.2% [22.1%, 24.4%]), and detections (24.3% [23.2%, 25.4%]). In both cases, the effect was similar for detections and declarations, suggesting that conducting more bag searches increases the rate of BRM being detected and encourages more declarations. As expected, increased passenger counts were also associated with increased interception rates across all response variables. Finally, there were also some differences between arrival airports in their estimated interception rates (see Suppl. material 1: B).

Random factors, flight number and flight origin, both explained some variance in interception rates (e.g., for N_Total, VFlightOrigin = 0.004, VFlightNumber = 0.042), although the random effects only explained a small proportion of variation relative to fixed effects (i.e., for N_Total, R2marginal = 0.543; R2conditional = 0.570, proportional VFlightOrigin = 0.002, and proportional VFlightNumber = 0.025). It should also be noted that the overdispersion may lead to overestimates of R2 values (e.g., Harrison 2014), so these values should be interpreted cautiously. Nonetheless, random intercept predictions for flight origin show how pathway factors may be used to identify and target interceptions towards higher risk arrivals (Fig. 3). Similar pathway heterogeneity can also be identified for flight number (see Suppl. material 1: figs B.1, B.2).

Figure 3.

Predicted random intercepts by flight origin for (A–C) biosecurity risk material (BRM) interceptions, declarations, and detections; and (B) fruit fly (FF) host item interceptions, declarations, and detections. Error bars represent 95CIs, and intervals that do not include zero are considered to have significantly higher or lower levels of BRM or FF host interceptions than an average flight. Intercept estimates are in the modelled unit, i.e., the log of the proportional difference between the group and the overall expected BRM or FF host count.

Model sensitivity

Tests showed that models used in the main analysis for N_Total and N_Declarations_FF both had probable zero-inflation (ratio of predicted to observed zeros: 0.85 and 0.92, respectively), and overdispersion was present in both cases (N_Total: dispersion ratio = 1.634, χ2 = 97786.221, P < 0.001; N_Declarations_FF: dispersion ratio = 1.709, χ2 = 102303.640, P < 0.001). Notably, overdispersion can be a common consequence of zero-inflation, in which case a zero-inflated Poisson approach may be sufficient to account for both issues (see Yang et al. 2009). Measures of model fit also showed that all models accounting for zero-inflation had lower WAIC scores than those that did not, whereas negative binomial models had the lowest scores and appeared to be the best-fitting models tested (see Suppl. material 1: table B.1).

Nonetheless, sensitivity analysis showed that the outputs were relatively robust to implementation methods, with the patterns identified qualitatively similar between implementation types, but with some variation in the magnitude and uncertainty of effects. Estimated interception rates under differing regimes showed similar patterns when using a Bayesian implementation, although with slightly greater uncertainty (e.g., the estimated rate with one BI was 0.50 [0.44, 0.57] compared to 0.50 [0.46, 0.55] in the main model; see Fig. 4). Incorporating zero-inflation into Poisson models led to higher rate estimates (see also Suppl. material 1: fig. B.3). Negative binomial models, which were the best-supported models based on WAIC scores, instead produced, considerably lower estimates of BRM interceptions, while still showing qualitatively similar differences when comparing interception rates between regimes.

Figure 4.

Estimated biosecurity risk material (BRM) interception rates for air passenger intervention regimes using different model implementations. Regimes include combinations of detector dog teams (DDT) and biosecurity inspectors (BI). Estimates are for five different model implementations, based on the modelled distribution and on a frequentist (lme4) vs Bayesian (brms) framework (see further details under Suppl. material 1: table B.1).

Random intercept predictions also appeared to show qualitatively similar patterns for FF host detections (Fig. 5) and BRM interceptions (Suppl. material 1: fig. B.4). Comparing outputs based on a Poisson distribution, the Bayesian (brms) approach produced greater uncertainty in the mean intercept estimates/predictions relative to the frequentist (lme4) implementation (Figs 4, 5), although the means themselves remained relatively consistent.

Figure 5.

Predicted random intercepts for fruit fly (FF) host item detections, by flight origin. Estimates are included from five different model implementations, which from the top include the following; (green, circle) Poisson-lme4; (purple, triangle) Poisson-brms; (blue, square) zero-inflated Poisson-brms; (gold, diamond) negative binomial-brms; and, (orange, upside-down triangle) zero-inflated negative binomial-brms. Intercept estimates are in the modelled unit, i.e., the log of the proportional difference between the group and the overall expected FF host count.

Discussion

Biosecurity interceptions at Tasmanian airports were strongly influenced by the methods used. For example, the number of bag searches conducted increased both detections and declarations from passengers. This suggests that increased effort in active at-border surveillance by officers will increase the efficacy of interventions by promoting both voluntary compliance by passengers and detections of undeclared risk items that may otherwise have been missed. Similarly, dog detector teams are increasingly deployed for border interventions along high-volume phytosanitary risk pathways, both in Australia and many other countries (Whattam et al. 2014; Inspector-General of Biosecurity 2022). Few studies have quantitatively assessed the efficacy of detector dogs, although a recent study from Williams and Sharp (2023) showed how the presence of a dog versus an officer alone can alter passenger behaviour including eye contact, gestures or interactions with the officer/dog. Our study provides further insights into the effects dogs may have on passenger behaviour and the rate of interceptions at airports. The increase also appears to be largely driven by detections of undeclared items instead of voluntary declarations, suggesting that dogs may be particularly useful for capturing a component of the biosecurity risk material that may otherwise not be found through more passive, voluntary compliance-based methods (e.g., public awareness and education campaigns, biosecurity signage and announcements). Although notably, one inspector alone still appears to capture about half of the biosecurity risk material entering the state relative to a detector dog team. Therefore, although the efficacy of a single person is lower than the detector dog, active human surveillance still effectively mitigates a proportion of risk at the border.

The ability to analyse pathway risk heterogeneity is limited by the type of data collected for arrivals along a biosecurity pathway. Despite the relatively limited set of pathway factors included in models (i.e., origin and route) and the relatively small proportion of total variance explained by these factors, models were able to identify specific flight origins and routes as potentially either high- or low-risk arrivals. This shows how interception data may be useful for supporting risk-based approaches to interventions (e.g., Australia’s flight-based-traveller profiles; Inspector-General of Biosecurity 2019) by identifying higher and lower risk arrivals to allocate limited resources to these arrivals (Trouvé et al. 2024). These analyses can provide important quantitative evidence supporting targeted resource allocations at the border, particularly when combined with further contextual information such as pre-border pest prevalence data or estimates of potential post-border impacts. While Tasmania currently targets 100% of air arrivals, risk-based approaches could include decisions about where and when to allocate their most effective methods (e.g., detector dogs), or to potentially identify a subset of low-risk arrivals based on factors such as origin or route that can be met with less resource-demanding methods (e.g., passive interventions, signage, amnesty bins).

While these data are valuable for identifying how interventions or pathway risk influence the actual interception rates, which is critical information for implementing biosecurity interventions to reduce risk at the border, further information is required to fully quantify the risk of incursions along this pathway. For example, the risk of fruit fly establishing through the air passenger pathway would require us to estimate the actual volume of biosecurity risk material on flights and the proportion of those items infested with fruit fly (i.e., contamination/infestation rates), the proportion of risk material missed (i.e., leakage), or the viability of any individuals or larvae that may infest any of the risk material. These parameters may be estimated using complementary methods, such as endpoint surveys that target a subset of passengers as manual or X-ray searches to estimate contamination rates (Mannix et al. 2024). Samples of intercepted biosecurity risk material may also be further tested to measure their pest contamination/infestation rates. In many cases, particularly in biosecurity, empirical data is lacking. Structured expert judgement may then be used to elicit unknown parameters from relevant experts (e.g., biosecurity managers, entomologists, academics, etc.), using advanced methods to directly elicit uncertainty in parameter estimates and incorporate this uncertainty into the decision-making process (Hemming et al. 2018; Bau et al. 2024). Therefore, although this study highlights the value of interception data for informing biosecurity practices, additional knowledge is required to more completely assess and quantify risk across a biosecurity continuum.

This study used a simple and common approach for count data (i.e., regression modelling based on a Poisson distribution), performed with modelling tools that are accessible, and relatively easy to implement. Sensitivity analysis suggested that the outputs of this approach were qualitatively similar to approaches using more advanced tools (i.e., Bayesian modelling methods), or distributions (e.g., that account for zero-inflation). However, there were some notable differences in outputs. For example, interception rates estimated from negative binomial models tended to be lower, suggesting that failing to account for over-dispersion may lead to overestimates. Furthermore, these models accounting for both zero-inflation and overdispersion (i.e., zero-inflated negative binomial models) were the best-fitting models tested. Also, Bayesian methods tended to lead to higher uncertainty estimates around fixed- and random-effect parameters, so may represent a more conservative approach to modelling pathway risk factors. Therefore, while our conclusions were generally robust to implementation methods, sensitivity analysis may also be a valuable step for providing additional information for decision-makers about the robustness of any conclusions drawn from modelling and future researchers should consider their implementation. In this case, however, the operational interpretation of the simpler models was borne out by the more complex models, i.e. that detector dogs were more effective than biosecurity inspectors.

Finally, effective at-border interventions are a key step in the biosecurity continuum. As at-border interventions become more sophisticated and widely implemented, large interception data sets will inevitably become more available to researchers and biosecurity decision-makers. This study highlights how this data can be a valuable resource for informing management decisions for Tasmania and can provide empirical evidence to support the implementation of risk-based approaches or the use of specific methods such as detector dogs, which can improve resource allocations and lead to more effective interventions at borders.

Acknowledgments

We acknowledge Biosecurity Tasmania and its staff for supporting this project and specifically thank Rae Burrows, Andrew Bishop, Ryan Wilkinson, Meg Flanagan, Renee Tiggs, Guy Westmore, Emma Hutchinson, Cameron Healey, Jason Paul, Alex Matthews, Liam Harrex, Jesse Siebler, Christian Knapp, Andrew Carter, Deepa Adhikari, Elise Baniowski, Marlene Quinn, and Chris Jones.

Additional information

Conflict of interest

The authors have declared that no competing interests exist.

Ethical statement

No ethical statement was reported.

Funding

The project received funding from Biosecurity Tasmania, Department of Natural Resources and Environment Tasmania.

Author contributions

NPM: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Resources, Validation, Visualization, Writing – original draft. AMH: Conceptualization, Methodology, Supervision, Writing – review & editing. APR: Conceptualization, Methodology, Funding acquisition, Supervision, Writing – review & editing.

Author ORCIDs

Nicholas. P Moran https://orcid.org/0000-0002-7331-0400

Anca M. Hanea https://orcid.org/0000-0003-3870-5949

Andrew P. Robinson https://orcid.org/0000-0002-0509-6043

Data availability

For privacy and operational purposes, all identifying information has been anonymised from datasets, including airport names, flight numbers, etc. Fully anonymised datasets, analysis code, models and outputs are all available at the Open Science Framework (osf.io/78tv9/; doi: 10.17605/OSF.IO/78TV9).

References

  • Bau S, Hanea A, Robinson AP, Burgman M (2024) Elicit: Using Structured Elicitation in Biosecurity. In: Hester SM, Bland LM, Arndt E, Bau S, Camac JS, Mannix E, Trouve R, Robinson AP (Eds) Biosecurity: A Systems Perspective. CRC Press, Boca Raton, 191–210. https://www.doi.org/10.1201/9781003253204-17
  • Bradshaw CJA, Hoskins AJ, Haubrock PJ, Cuthbert RN, Diagne C, Leroy B, Andrews L, Page B, Cassey P, Sheppard AW, Courchamp F (2021) Detailed assessment of the reported economic costs of invasive species in Australia. NeoBiota 67: 511–550. https://doi.org/10.3897/neobiota.67.58834
  • Brettell LE, Martin SJ, Riegler M, Cook JM (2021) Vulnerability of island insect pollinator communities to pathogens. Journal of Invertebrate Pathology 186: 107670. https://doi.org/10.1016/j.jip.2021.107670
  • Campbell H (2021) The consequences of checking for zero-inflation and overdispersion in the analysis of count data. Methods in Ecology and Evolution 12(4): 665–680. https://doi.org/10.1111/2041-210X.13559
  • Cook DC, Fraser RW (2015) Eradication versus control of Mediterranean fruit fly in Western Australia. Agricultural and Forest Entomology 17(2): 173–180. https://doi.org/10.1111/afe.12093
  • Dodd AJ, Burgman MA, McCarthy MA, Ainsworth N (2015) The changing patterns of plant naturalization in Australia. Diversity & Distributions 21(9): 1038–1050. https://doi.org/10.1111/ddi.12351
  • Epanchin-Niell R, McAusland C, Liebhold A, Mwebaze P, Springborn MR (2021) Biological Invasions and International Trade: Managing a Moving Target. Review of Environmental Economics and Policy 15(1): 180–190. https://doi.org/10.1086/713025
  • Harrison XA (2014) Using observation-level random effects to model overdispersion in count data in ecology and evolution. PeerJ 2: e616. https://doi.org/10.7717/peerj.616
  • Hemming V, Burgman MA, Hanea AM, McBride MF, Wintle BC (2018) A practical guide to structured expert elicitation using the IDEA protocol. Methods in Ecology and Evolution 9(1): 169–180. https://doi.org/10.1111/2041-210X.12857
  • Inspector-General of Biosecurity (2022) Efficacy and adequacy of department’s X-ray scanning and detector dog screening techniques to prevent entry of biosecurity risk material into Australia. Department of Agriculture and Water Resources, Canberra, Australia. https://www.igb.gov.au/current-and-completed-reviews [Accessed on 12.09.2023]
  • Kachigunda B, Mengersen K, Perera DI, Coupland GT, Van der Merwe J, McKirdy S (2022) Use of mixed-type data clustering algorithm for characterizing temporal and spatial distribution of biosecurity border detections of terrestrial non-indigenous species. PLoS ONE 17(8): e0272413. https://doi.org/10.1371/journal.pone.0272413
  • Keitt B, Campbell K, Saunders A, Clout M, Wang Y, Heinz R, Newton K, Tershy B (2011) The global islands invasive vertebrate eradication database: a tool to improve and facilitate restoration of island ecosystems. In: Veitch CR, Clout MN, Towns DR (Eds) Island Invasives: Eradication and Management: Proceedings of the International Conference on Island Invasives. International Union for Conservation of Nature (IUCN) Gland, Switzerland, 74–77.
  • Kier G, Kreft H, Lee TM, Jetz W, Ibisch PL, Nowicki C, Mutke J, Barthlott W (2009) A global assessment of endemism and species richness across island and mainland regions. Proceedings of the National Academy of Sciences of the United States of America 106(23): 9322–9327. https://doi.org/10.1073/pnas.0810306106
  • Liebhold AM, Work TT, McCullough DG, Cavey JF (2006) Airline Baggage as a Pathway for Alien Insect Species Invading the United States. American Entomologist (Lanham, Md. ) 52(1): 48–54. https://doi.org/10.1093/ae/52.1.48
  • Lindén A, Mäntyniemi S (2011) Using the negative binomial distribution to model overdispersion in ecological count data. Ecology 92(7): 1414–1421. https://doi.org/10.1890/10-1831.1
  • Lüdecke D, Ben-Shachar MS, Patil I, Waggoner P, Makowski D (2021) performance: An R package for assessment, comparison and testing of statistical models. Journal of Open Source Software 6(60): 3139. https://doi.org/10.21105/joss.03139
  • Mannix E, Baumgartner JB, Bau S, Bland LM, Robinson AP, Page N (2024) Profiling and Automation. In: Hester SM, Bland LM, Arndt E, Bau S, Camac JS, Mannix E, Trouve R, Robinson AP (Eds) Biosecurity: A Systems Perspective. CRC Press, Boca Raton, 211–230. https://doi.org/10.1201/9781003253204-18
  • McCullough DG, Work TT, Cavey JF, Liebhold AM, Marshall D (2006) Interceptions of Nonindigenous Plant Pests at US Ports of Entry and Border Crossings Over a 17-year Period. Biological Invasions 8(4): 611–630. https://doi.org/10.1007/s10530-005-1798-4
  • McNeill M, Phillips C, Young S, Shah F, Aalders L, Bell N, Gerard E, Littlejohn R (2011) Transportation of nonindigenous species via soil on international aircraft passengers’ footwear. Biological Invasions 13(12): 2799–2815. https://doi.org/10.1007/s10530-011-9964-3
  • Moir ML, Croeser L, Telfer D, Fenner C, McCauley R (2022) Value-adding in biosecurity surveillance and monitoring: Testing colour and non-target semiochemical lures on Psylloidea and Pentatomoidea. Journal of Applied Entomology 146(10): 1333–1342. https://doi.org/10.1111/jen.13074
  • Moser AY, Brown WY, Bizo LA, Andrew NR, Taylor MK (2020) Biosecurity Dogs Detect Live Insects after Training with Odor-Proxy Training Aids: Scent Extract and Dead Specimens. Chemical Senses 45(3): 179–186. https://doi.org/10.1093/chemse/bjaa001
  • Pace R, Ascolese R, Miele F, Russo E, Griffo RV, Bernardo U, Nugnes F (2022) The Bugs in the Bags: The Risk Associated with the Introduction of Small Quantities of Fruit and Plants by Airline Passengers. Insects 13(7): 617. https://doi.org/10.3390/insects13070617
  • Pittman B, Buta E, Garrison K, Gueorguieva R (2022) Models for Zero-Inflated and Overdispersed Correlated Count Data: An Application to Cigarette Use. Nicotine & Tobacco Research: Official Journal of the Society for Research on Nicotine and Tobacco 25(5): 996–1003. https://doi.org/10.1093/ntr/ntac253
  • Potts BM, Sandhu KS, Wardlaw T, Freeman J, Li H, Tilyard P, Park RF (2016) Evolutionary history shapes the susceptibility of an island tree flora to an exotic pathogen. Forest Ecology and Management 368: 183–193. https://doi.org/10.1016/j.foreco.2016.02.027
  • Robinson AP, McNeill MR (2022) Biosecurity and post-arrival pathways in New Zealand: Relating alien organism detections to tourism indicators. NeoBiota 71: 51–69. https://doi.org/10.3897/neobiota.71.64618
  • Schielzeth H, Dingemanse NJ, Nakagawa S, Westneat DF, Allegue H, Teplitsky C, Réale D, Dochtermann NA, Garamszegi LZ, Araya-Ajoy YG (2020) Robustness of linear mixed-effects models to violations of distributional assumptions. Methods in Ecology and Evolution 11(9): 1141–1152. https://doi.org/10.1111/2041-210X.13434
  • Seebens H, Blackburn TM, Dyer EE, Genovesi P, Hulme PE, Jeschke JM, Pagad S, Pyšek P, Winter M, Arianoutsou M, Bacher S, Blasius B, Brundu G, Capinha C, Celesti-Grapow L, Dawson W, Dullinger S, Fuentes N, Jäger H, Kartesz J, Kenis M, Kreft H, Kühn I, Lenzner B, Liebhold A, Mosena A, Moser D, Nishino M, Pearman D, Pergl J, Rabitsch W, Rojas-Sandoval J, Roques A, Rorke S, Rossinelli S, Roy HE, Scalera R, Schindler S, Štajerová K, Tokarska-Guzik B, van Kleunen M, Walker K, Weigelt P, Yamanaka T, Essl F (2017) No saturation in the accumulation of alien species worldwide. Nature Communications 8(1): 14435. https://doi.org/10.1038/ncomms14435
  • Seebens H, Bacher S, Blackburn TM, Capinha C, Dawson W, Dullinger S, Genovesi P, Hulme PE, van Kleunen M, Kühn I, Jeschke JM, Lenzner B, Liebhold AM, Pattison Z, Pergl J, Pyšek P, Winter M, Essl F (2021) Projecting the continental accumulation of alien species through to 2050. Global Change Biology 27(5): 970–982. https://doi.org/10.1111/gcb.15333
  • Sequeira R, Griffin R (2014) The biosecurity continuum and trade: Pre-border operations. In: Gordh G, McKirdy S (Eds) The Handbook of Plant Biosecurity: Principles and Practices for the Identification, Containment and Control of Organisms that Threaten Agriculture and the Environment Globally. Springer, Dordrecht, 119–148. https://doi.org/10.1007/978-94-007-7365-3_5
  • Skinner W (2018) Presence Through Absence: Phylloxera and the Viticultural Imagination in McLaren Vale, South Australia. The Asia Pacific Journal of Anthropology 19(3): 245–263. https://doi.org/10.1080/14442213.2018.1461916
  • Stoeckl N, Dodd A, Kompas T (2023) The monetary value of 16 services protected by the Australian National Biosecurity System: Spatially explicit estimates and vulnerability to incursions. Ecosystem Services 60: 101509. https://doi.org/10.1016/j.ecoser.2023.101509
  • Trouvé R, Robinson AP (2021) Estimating Consignment-Level Infestation Rates from the Proportion of Consignment that Failed Border Inspection: Possibilities and Limitations in the Presence of Overdispersed Data. Risk Analysis 41(6): 992–1003. https://doi.org/10.1111/risa.13592
  • Trouvé R, Bland LM, Robinson AP, Ducey MJ, Hester SM (2024) Screen: Designing Sampling Schemes for Border Inspection. In: Hester SM, Bland LM, Arndt E, Bau S, Camac JS, Mannix E, Trouve R, Robinson AP (Eds) Biosecurity: A Systems Perspective. CRC Press, Boca Raton, 69–82. https://doi.org/10.1201/9781003253204-7
  • Turner RM, Plank MJ, Brockerhoff EG, Pawson S, Liebhold A, James A (2020) Considering unseen arrivals in predictions of establishment risk based on border biosecurity interceptions. Ecological Applications 30(8): e02194. https://doi.org/10.1002/eap.2194
  • Turner RM, Brockerhoff EG, Bertelsmeier C, Blake RE, Caton B, James A, MacLeod A, Nahrung HF, Pawson SM, Plank MJ, Pureswaran DS, Seebens H, Yamanaka T, Liebhold AM (2021) Worldwide border interceptions provide a window into human-mediated global insect movement. Ecological Applications 31(7): e02412. https://doi.org/10.1002/eap.2412
  • Whattam M, Clover G, Firko M, Kalaris T (2014) The biosecurity continuum and trade: Border operations. In: Gordh G, McKirdy S (Eds) The Handbook of Plant Biosecurity: Principles and Practices for the Identification, Containment and Control of Organisms that Threaten Agriculture and the Environment Globally. Springer, Dordrecht, 149–188. https://doi.org/10.1007/978-94-007-7365-3_6
  • Whattam M, Azzopardi S, Nehl D, Maxwell A, Davis K (2024) Protecting Australia’s plant health: Plant quarantine in an evolving biosecurity system. Historical Records of Australian Science 36. https://doi.org/10.1071/HR24012
  • Williams EEM, Sharp RA (2023) Some effects of detection dogs on passenger behavior at border control ports. Journal of Applied Behavior Analysis 56(2): 377–387. https://doi.org/10.1002/jaba.985

Supplementary material

Supplementary material 1 

Additional supporting materials include the following data and modelling outputs

Nicholas. P Moran Anca M. Hanea, Andrew P. Robinson

Data type: pdf

Explanation note: (A) Composition of BRM and FF host interceptions; (B) Supplementary model design and output details.

This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.
Download file (5.42 MB)
login to comment