<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Insect Sci.</journal-id>
<journal-title>Frontiers in Insect Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Insect Sci.</abbrev-journal-title>
<issn pub-type="epub">2673-8600</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/finsc.2022.783285</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Insect Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Quantifying Spillover of an Urban Invasive Vector of Plant Disease: Asian Citrus Psyllid (<italic>Diaphorina citri</italic>) in California Citrus</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name><surname>Bayles</surname> <given-names>Brett R.</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1353425/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Thomas</surname> <given-names>Shyam M.</given-names></name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1500519/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Simmons</surname> <given-names>Gregory S.</given-names></name>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Daugherty</surname> <given-names>Matthew P.</given-names></name>
<xref ref-type="aff" rid="aff5"><sup>5</sup></xref>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Department of Global Public Health, Dominican University of California</institution>, <addr-line>San Rafael, CA</addr-line>, <country>United States</country></aff>
<aff id="aff2"><sup>2</sup><institution>Department of Natural Sciences and Mathematics, Dominican University of California</institution>, <addr-line>San Rafael, CA</addr-line>, <country>United States</country></aff>
<aff id="aff3"><sup>3</sup><institution>Department of Fisheries, Wildlife and Conservation Biology, University of Minnesota</institution>, <addr-line>St. Paul, MN</addr-line>, <country>United States</country></aff>
<aff id="aff4"><sup>4</sup><institution>United States Department of Agriculture, Animal and Plant Health Inspection Service</institution>, <addr-line>Salinas, CA</addr-line>, <country>United States</country></aff>
<aff id="aff5"><sup>5</sup><institution>Department of Entomology, University of California, Riverside</institution>, <addr-line>Riverside, CA</addr-line>, <country>United States</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Robert C. Venette, Forest Service, United States Department of Agriculture (USDA), United States</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Mamoudou Setamou, Texas A&#x00026;M University Kingsville, United States; Bob Douma, Wageningen University and Research, Netherlands</p></fn>
<corresp id="c001">&#x0002A;Correspondence: Brett R. Bayles <email>brett.bayles&#x00040;dominican.edu</email></corresp>
<fn fn-type="other" id="fn001"><p>This article was submitted to Invasive Insect Species, a section of the journal Frontiers in Insect Science</p></fn></author-notes>
<pub-date pub-type="epub">
<day>16</day>
<month>02</month>
<year>2022</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>2</volume>
<elocation-id>783285</elocation-id>
<history>
<date date-type="received">
<day>25</day>
<month>09</month>
<year>2021</year>
</date>
<date date-type="accepted">
<day>12</day>
<month>01</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2022 Bayles, Thomas, Simmons and Daugherty.</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>Bayles, Thomas, Simmons and Daugherty</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license> 
</permissions>
<abstract>
<p>Urban environments frequently play an important role in the initial stages of biological invasions, often serving as gateways for non-native species, which may propagate to nearby natural and agricultural ecosystems in the event of spillover. In California, citrus trees are a dominant ornamental and food plant in urban and peri-urban environments. We studied the invasion dynamics of the Asian citrus psyllid (<italic>Diaphorina citri</italic>), which became widespread in urban areas of southern California starting in 2008, to understand the factors driving its more recent invasion in commercial citrus groves. Using a multi-year monitoring database, we applied a suite of models to evaluate the rate at which groves accrued their first <italic>D. citri</italic> detection and the cumulative number of detections thereafter. Grove characteristics and landscape context proved to be important, with generally higher invasion rates and more cumulative detections in groves that were larger, had more edge, or had more perforated shapes, with greater urbanization intensity favoring more rapid invasion, but with inconsistent effects of distance to roads among models. Notably, distance to urban or other grove occurrences proved to be among the most important variables. During the early phase of <italic>D. citri</italic> invasion in the region, groves closer to urban occurrences were invaded more rapidly, whereas more recently, invasion rate depended primarily on proximity to grove occurrences. Yet, proximity to urban and grove occurrences contributed positively to cumulative <italic>D. citri</italic> detections, suggesting a continued influx from both sources. These results suggest that inherent features of agroecosystems and spatial coupling with urban ecosystems can be important, temporally dynamic, drivers of biological invasions. Further consideration of these issues may guide the development of strategic responses to <italic>D. citri</italic>&#x00027;s ongoing invasion.</p></abstract>
<kwd-group>
<kwd>biological invasion</kwd>
<kwd>invasion dynamics</kwd>
<kwd>landscape ecology</kwd>
<kwd>urban-rural interface</kwd>
<kwd>herbivore spillover</kwd>
<kwd>citrus greening disease</kwd>
<kwd>huanglongbing</kwd>
</kwd-group>
<counts>
<fig-count count="7"/>
<table-count count="0"/>
<equation-count count="1"/>
<ref-count count="49"/>
<page-count count="10"/>
<word-count count="6566"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>Introduction</title>
<p>Urban ecosystems provide a suite of services for human health and well-being, including food production from urban gardens and commercial agriculture (<xref ref-type="bibr" rid="B1">1</xref>&#x02013;<xref ref-type="bibr" rid="B3">3</xref>). However, such benefits should be weighed against potential disservices created as a result of these human-altered landscapes (<xref ref-type="bibr" rid="B4">4</xref>, <xref ref-type="bibr" rid="B5">5</xref>). In particular, urban areas contain unique sets of social and ecological processes that may create focal points for the introduction, establishment and spread of non-native species (<xref ref-type="bibr" rid="B6">6</xref>&#x02013;<xref ref-type="bibr" rid="B10">10</xref>).</p>
<p>Anthropogenic landscapes often exist along gradients of different land-use types, including urban and peri-urban cities, commercial agriculture, and a range of natural ecosystems (<xref ref-type="bibr" rid="B8">8</xref>, <xref ref-type="bibr" rid="B11">11</xref>, <xref ref-type="bibr" rid="B12">12</xref>). The particular spatial configurations of these landscapes in a given area are idiosyncratic, which suggests that a diverse interplay of social and ecological processes determine the initial risk of invasion, as well as if and when spillover into surrounding urban forest landscapes will occur (<xref ref-type="bibr" rid="B4">4</xref>, <xref ref-type="bibr" rid="B10">10</xref>, <xref ref-type="bibr" rid="B13">13</xref>, <xref ref-type="bibr" rid="B14">14</xref>).</p>
<p>Cross-habitat spillover involving the movement of invasive species between agroecosystems and natural habitats have been documented (<xref ref-type="bibr" rid="B15">15</xref>, <xref ref-type="bibr" rid="B16">16</xref>). Most studies have focused on asymmetric movement of insects between agricultural areas and surrounding natural habitats (<xref ref-type="bibr" rid="B13">13</xref>, <xref ref-type="bibr" rid="B15">15</xref>, <xref ref-type="bibr" rid="B17">17</xref>). Landscape features, such as the extent and characteristics of habitat edges and the relative proximity to source habitats, have been shown to affect spillover across habitats (<xref ref-type="bibr" rid="B13">13</xref>). Such spillover events can further influence community dynamics of the recipient habitat through processes such as parasitism, predation, and disease transmission (<xref ref-type="bibr" rid="B13">13</xref>, <xref ref-type="bibr" rid="B14">14</xref>, <xref ref-type="bibr" rid="B18">18</xref>, <xref ref-type="bibr" rid="B19">19</xref>).</p>
<p>Southern California, a region characterized by extensive urban development, has recently experienced the introduction of numerous phytophagous insects (<xref ref-type="bibr" rid="B20">20</xref>). Of the &#x0007E;10 exotic arthropod introduction events per year in the state, an estimated 20% become invasive pests (<xref ref-type="bibr" rid="B21">21</xref>). One invader of particular concern is the Asian citrus psyllid (<italic>Diaphorina citri</italic>) because of its ability to transmit the bacterium <italic>Candidatus</italic> Liberibacter asiaticus (<italic>C</italic>Las) associated with huanglongbing (HLB; citrus greening disease) [Grafton-Cardwell et al., (<xref ref-type="bibr" rid="B22">22</xref>)]. Symptoms of HLB, to which all commercial varieties of citrus are considered susceptible, include progressive yellowing of leaves, stunting of shoots, defoliation, reduced yield and fruit quality, and tree mortality (<xref ref-type="bibr" rid="B23">23</xref>). Citrus growing regions where both the pathogen and vector successfully invaded, such as Florida, have faced severe economic impacts (<xref ref-type="bibr" rid="B24">24</xref>).</p>
<p>Citrus has unique historical and cultural significance in California, following its likely introduction into the region by Spanish missionaries in the late 18th century (<xref ref-type="bibr" rid="B25">25</xref>). The first commercial production began in southern California in the mid-19th century, expanded to other areas of the state in subsequent decades, and eventually developed into what is currently estimated to be approximately 270,000 ac of citrus farmed at a yearly production value of more than $3B (<xref ref-type="bibr" rid="B25">25</xref>&#x02013;<xref ref-type="bibr" rid="B27">27</xref>). The development of the citrus industry is credited with contributing to California&#x00027;s population expansion throughout the 19th century, the creation or evaluation of novel citrus varieties that are commonly grown today, and the widespread acceptance of biological control as an element of integrated pest management (<xref ref-type="bibr" rid="B25">25</xref>, <xref ref-type="bibr" rid="B28">28</xref>). Enthusiasm for citrus carries over to urban and other residential areas of California where it is frequently grown as an ornamental and food plant. Surveys conducted in select neighborhoods in southern California estimated that more than 60% of properties had at least one citrus tree present (<xref ref-type="bibr" rid="B29">29</xref>). Thus, residential and commercial trees are both at risk to invasive pests and diseases, such as <italic>D. citri</italic> and huanglongbing, and indeed may exacerbate damage to each other in the event of cross-habitat spillover (<xref ref-type="fig" rid="F1">Figure 1</xref>).</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p>Map of Southern California showing the distribution of commercial citrus groves and surrounding urban land cover.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="finsc-02-783285-g0001.tif"/>
</fig>
<p><italic>D. citri</italic> was first detected in California in a residential area of San Diego County in 2008 and in Los Angeles by 2009 (<xref ref-type="bibr" rid="B30">30</xref>). In response, a strategic plan was implemented to mitigate the impact of this invasive pathosystem in California (<xref ref-type="bibr" rid="B27">27</xref>). Over the next few years, the insect spread throughout urban and peri-urban areas, apparently driven by favorable climatic conditions and landscape characteristics in these environments (<xref ref-type="bibr" rid="B30">30</xref>, <xref ref-type="bibr" rid="B31">31</xref>). <italic>D. citri</italic> was first detected in a handful of commercial citrus groves in 2011, and over the next few years became increasingly prevalent in these settings (<xref ref-type="bibr" rid="B30">30</xref>). Finally, the first case of HLB was detected in a residential area of Los Angeles County in 2012 (<xref ref-type="bibr" rid="B32">32</xref>). Disease surveys have subsequently identified upward of 2500 infected trees that have been removed to limit pathogen spread (<xref ref-type="bibr" rid="B27">27</xref>, <xref ref-type="bibr" rid="B33">33</xref>). Thus far, all cases of the disease have occurred in residential areas, but there was one detection of a <italic>C</italic>Las positive psyllid from a commercial grove, and nearby groves are considered at risk.</p>
<p>While invasive species spillover between agricultural and natural ecosystems has been quantified, less is known about how common spillover is from urban to agricultural ecosystems, or the processes governing it across space and time. The general sequence of <italic>D. citri</italic> occurrence and patterns of <italic>C</italic>Las spread in Southern California are consistent with urban environments playing a potentially important role in facilitating the establishment of <italic>D. citri</italic> (<xref ref-type="bibr" rid="B30">30</xref>, <xref ref-type="bibr" rid="B31">31</xref>, <xref ref-type="bibr" rid="B33">33</xref>). Yet, formal analyses are lacking that identify the factors driving <italic>D. citri</italic> invasion of commercial citrus groves, which is critically important for refining disease detection and management strategies. Here, we quantify the magnitude of urban-to-agriculture <italic>D. citri</italic> spillover, and the underlying role of grove-level characteristics, urban landscape-context, and neighborhood proximity in governing the spatial and temporal distribution of <italic>D. citri</italic> activity in commercial groves.</p>
</sec>
<sec sec-type="materials and methods" id="s2">
<title>Materials and Methods</title>
<sec>
<title>Data</title>
<p>We obtained a georeferenced database of 78,000 <italic>D. citri</italic> panel trap collection sites in Southern California from the California Department of Food and Agriculture (CDFA) and the Citrus Pest and Disease Prevention Program (CPDPP) via the Citrus Research Board (CRB) from 2008 through 2014. Detection records were obtained from 14 x 23 cm double-sided yellow panel traps that were positioned in parts of 10 counties (Imperial, Kern, Los Angeles, Orange, Riverside, Santa Barbara, San Bernardino, San Diego, and San Luis Obispo) [see (<xref ref-type="bibr" rid="B30">30</xref>) for a detailed description of the trapping database].</p>
<p>We developed three indicators of &#x0201C;grove-level&#x0201D; characteristics, including: (1) grove area, (2) grove edge, and (3) grove perforation. Grove area was calculated as the size of each commercial citrus grove in square meters using a commercial citrus grove (definition from 26) layer (obtained from the CRB) in ArcGIS v. 10.2. Given that the number of traps varied with the size of the grove, we expect to find an area effect in the analyses, and therefore use it in part to minimize confounding effects on the estimation of other variables in the models.</p>
<p>Because there is a biologically plausible mechanism by which an &#x0201C;edge effect&#x0201D; in citrus groves can impact insect abundance (<xref ref-type="bibr" rid="B34">34</xref>) and disease transmission risk (<xref ref-type="bibr" rid="B23">23</xref>), we used the Landscape Fragmentation Tool (LFT v2.0) to quantify grove edge and grove perforation (<xref ref-type="bibr" rid="B35">35</xref>). Grove edge was calculated as the proportion of citrus tree pixels within 100 meters of the citrus grove/non-citrus grove outside edge. Grove perforation was calculated as the proportion of citrus tree pixels along the edge of small clearings found within groves. Together, these landscape metrics characterize the extent to which citrus groves are fragmented and are therefore an effective proxy for quantifying the edge effect (<xref ref-type="bibr" rid="B36">36</xref>).</p>
<p>We developed two indicators of &#x0201C;urban landscape-level&#x0201D; characteristics surrounding groves, including: (1) urbanization intensity, and (2) proximity to transportation corridors (i.e., distance to roads). Urbanization was calculated as the proportion of both treeless cover and impervious surface within 100 m of each grove:</p>
<disp-formula id="E1"><mml:math id="M1"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:mtext>&#x000A0;</mml:mtext><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mn>100</mml:mn><mml:mi>%</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mo>-</mml:mo><mml:mtext>&#x000A0;</mml:mtext><mml:mi>%</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>T</mml:mi><mml:mi>r</mml:mi><mml:mi>e</mml:mi><mml:mi>e</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>C</mml:mi><mml:mi>o</mml:mi><mml:mi>v</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mtext>&#x000A0;</mml:mtext><mml:mo>&#x0002B;</mml:mo><mml:mtext>&#x000A0;</mml:mtext><mml:mi>%</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>I</mml:mi><mml:mi>m</mml:mi><mml:mi>p</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mi>v</mml:mi><mml:mi>i</mml:mi><mml:mi>o</mml:mi><mml:mi>u</mml:mi><mml:mi>s</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>S</mml:mi><mml:mi>u</mml:mi><mml:mi>r</mml:mi><mml:mi>f</mml:mi><mml:mi>a</mml:mi><mml:mi>c</mml:mi><mml:mi>e</mml:mi></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:mfrac></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>This produced an index of urbanization intensity ranging from 0 (low urbanization) to 100 (high urbanization). Proximity to transportation corridors was calculated as the minimum distance between groves and the nearest road designated for commercial truck traffic by the California Department of Transportation (available from <ext-link ext-link-type="uri" xlink:href="http://www.dot.ca.gov/hq/tsip/gis/datalibrary/">http://www.dot.ca.gov/hq/tsip/gis/datalibrary/</ext-link>).</p>
<p>Finally, we developed a set of indicators to measure the &#x0201C;neighbor effect&#x0201D; to more mechanistically capture spatial linkages stemming from urban influences. We included measures of distance to urban or peri-urban <italic>D. citri</italic> detections to capture any spillover from urban environments, and distance to <italic>D. citri</italic> detections in other citrus groves to capture grove-to-grove spillover. Due to strong year-to-year collinearity among distances to urban <italic>D. citri</italic> detections and among distances to <italic>D. citri</italic> grove detections, we grouped neighbor effect indicators into two temporal phases relative to when <italic>D. citri</italic> detections in urban areas and in groves started to become common (<xref ref-type="bibr" rid="B30">30</xref>). Our indicators of neighbor effect included: (1) early-phase (2008&#x02013;2011) distance between citrus groves and nearest residential citrus detection (i.e., &#x0201C;urban early-phase&#x0201D;), (2) late-phase (2012&#x02013;2014) distance between citrus groves and nearest residential citrus detection (i.e., &#x0201C;urban late-phase&#x0201D;), (3) early-phase (2011&#x02013;2012) distance between commercial citrus groves and nearest commercial grove detection (i.e., &#x0201C;grove early-phase&#x0201D;), and (4) late-phase (2013&#x02013;2014) distance between commercial citrus groves and nearest commercial grove detection (i.e., &#x0201C;grove late-phase&#x0201D;). An initial principal components analysis supported grouping the neighbor effects in this manner (MP Daugherty, <italic>unpublished data</italic>).</p>
</sec>
<sec>
<title>Data Analysis</title>
<p>We used two sets of analyses to identify the factors related to <italic>D. citri:</italic> (1) the invasion rate of commercial citrus groves (i.e., why some groves were invaded more quickly than others), and (2) the cumulative number of <italic>D. citri</italic> detections in all invaded groves as a proxy for the potential impact of this invasive vector (i.e., why some groves have more repeat detections than others). Cumulative detections were calculated as the total number of traps in a grove, to date, that collected at least one <italic>D. citri</italic> (i.e., total number of traps with <italic>D. citri</italic> present). We adopted this approach, instead of using explicit <italic>D. citri</italic> counts on each trap, due to concerns about inconsistency in reporting counts versus simply reporting <italic>D. citri</italic> presence absence. Prior to analysis, all grove, landscape, and neighbor effects were standardized by converting them to z-scores. This standardization is important to ensure that the different units in which the variables were measured (e.g., grove area [m<sup>2</sup>] vs. urbanization [0&#x02013;100%]) did not bias the model fitting and selection process in favor of particular variables (<xref ref-type="bibr" rid="B37">37</xref>). All analyses were conducted using the R programming language (v. 3.6.1).</p>
<p>For the analysis of invasion rate, we developed a pair of Cox proportional hazard models, which are semi-parametric tests often used to test &#x0201C;time-to-event&#x0201D; data that includes multiple variables (<xref ref-type="bibr" rid="B38">38</xref>), to discern the effects of grove-level variables in an urban landscape-level context from that of neighbor effects that capture spillover pressure [(<xref ref-type="bibr" rid="B39">39</xref>); coxph() in R]. The first model evaluated the effects of grove-level characteristics (i.e., grove area, grove edge, grove perforation) and urban landscape-level (i.e., urbanization intensity, distance to roads) variables on the rate at which commercial groves incurred their first <italic>D. citri</italic> detection. Because the distance to roads variable did not meet the proportional hazards assumption, an interaction with time was included in the model (<xref ref-type="bibr" rid="B40">40</xref>). Next, to capture more mechanistically the effects of potential spatially dependent drivers of invasion rate, a second model evaluated effects of grove-level characteristics and neighbor effects [i.e., grove early-phase (2011&#x02013;2012); urban early-phase (2008&#x02013;2011); grove late-phase (2013&#x02013;2014); urban late-phase (2012&#x02013;2014)] on invasion rate. Due to collinearity with some of the neighbor effects, the urban landscape-level effects were not included in the second model. For both models, we fit different formulations with all possible combinations of the 5 or 7 unique variables. The different model formulations were then ranked via Akaike Information Criterion [AIC; (<xref ref-type="bibr" rid="B41">41</xref>)] to identify one or more best-supported models, which we define as falling within five AIC units (i.e. &#x00394;AIC &#x0003C; 5) of the highest ranked model (i.e. lowest AIC value). Finally, we used model averaging with shrinkage [model.avg() in package MuMLn; Barto&#x00144;, (<xref ref-type="bibr" rid="B42">42</xref>)] of the best-supported models to estimate the effect size and 95% confidence interval for each variable represented in at least one of the best-supported models. In the event that a single model substantially outperformed the others, we evaluated the significance of individual terms with a series of likelihood-ratio tests (<xref ref-type="bibr" rid="B38">38</xref>).</p>
<p>Unlike invasion rate, cumulative number of <italic>D. citri</italic> detections in a citrus grove is a proxy measure of vector impact, stemming from effects of vector density or frequency relative to hosts on disease outbreaks (<xref ref-type="bibr" rid="B43">43</xref>). We used a similar approach to analyze the effects of grove-level characteristics, urban landscape-level context, and neighbor effects on the cumulative number of grove detections (i.e. number of <italic>D. citri</italic> positive traps) for the approximately 500 groves that had at least one detection. Due the count nature of cumulative detections, and high frequency of groves with a low number of observed detections, they were analyzed with a pair of generalized linear models with negative binomial error (<xref ref-type="bibr" rid="B38">38</xref>). As before, one set of models analyzed grove-level characteristics and urban landscape-level context variables and a second analyzed grove-level characteristics and neighbor effects. In addition, each set of models included an effect of time since infestation (i.e., number of months since the first <italic>D. citri</italic> detection in that grove) to control for some groves being first invaded much (i.e., up to 3 years) earlier than others. Again, we fit all possible combinations of variables within each set of models, ranked those models using an information-theoretic approach, and used model averaging to estimate the effect size of coefficients in the best-supported models.</p>
</sec>
</sec>
<sec sec-type="results" id="s3">
<title>Results</title>
<p>Of the nearly 2,500 unique citrus groves for which <italic>D. citri</italic> trapping data was available, 539 had at least one detection by the end of 2014. The earliest commercial citrus detection occurred in April 2011, with six additional groves having detections in that year, an additional 208 groves in 2012, and the remainder in 2013 and 2014. Of those groves in which <italic>D. citri</italic> was detected, 6.5% (35/538) had a single trap detection, but with a maximum of 359 cumulative detections (mean &#x000B1; SD: 18.76 &#x000B1; 33.46).</p>
<sec>
<title>Invasion Rate</title>
<p>The analysis of grove characteristics and landscape context on rate of invasion revealed the full model with all five variables and an interaction between distance to roads and time provided the best description of the data (&#x0003E; 5 AIC units lower than all other models). As a result, we did not use model averaging for this part of the analysis. The full model included significant effects of grove area (&#x003C7;<sup>2</sup> = 5.637, <italic>df</italic> = 1, <italic>p</italic> = 0.0176), grove edge (&#x003C7;<sup>2</sup> = 70.664, <italic>df</italic> = 1, <italic>p</italic> &#x0003C; 0.0001), grove perforation (&#x003C7;<sup>2</sup> = 23.423, <italic>df</italic> = 1, <italic>p</italic> &#x0003C; 0.0001), distance to roads (&#x003C7;<sup>2</sup> = 13.9622, <italic>df</italic> = 1, <italic>p</italic> = 0.0002), and urbanization intensity (&#x003C7;<sup>2</sup> = 17.672, <italic>df</italic> = 1, <italic>p</italic> &#x0003C; 0.0001), and a road distance-by-time interaction (&#x003C7;<sup>2</sup> = 19.558, <italic>df</italic> = 1, <italic>p</italic> &#x0003C; 0.0001). With respect to the grove-level characteristics, grove area [hazard ratio (SE) = 1.063 (0.017)], amount of grove edge [1.391 (0.04)], and the extent of grove perforation [1.15 (0.026)] were all positively related to the rate at which individual citrus groves incurred their first <italic>D. citri</italic> detection (i.e., hazard ratio &#x0003E; 1; <xref ref-type="fig" rid="F2">Figure 2</xref>). Both of the urban landscape-context variables were also positively related to invasion rate [urbanization intensity: 1.197 (0.042); distance to roads: 2.167 (0.188)]. However, the road-by-time interaction with a hazard ratio &#x0003C;one [0.965 (0.003)] indicates that early on groves further from roadways were invaded at equivalent or lower rates than groves distant from roadways (<xref ref-type="fig" rid="F3">Figure 3</xref>).</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p>Effect of <bold>(A)</bold> grove area, <bold>(B)</bold> amount of grove edge, or <bold>(C)</bold> grove perforation on the time to first <italic>D. citri</italic> detection in a grove. Kaplan-Meir plots were created by grouping based on natural breaks in the independent variables, resulting in three relative grove sizes and two relative values for edge and perforation. The x-axis reflects the amount of time since the first detection of <italic>D. citri</italic> in California in 2008.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="finsc-02-783285-g0002.tif"/>
</fig>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p>Effect of <bold>(A)</bold> distance to major roads, or <bold>(B)</bold> surrounding urban intensity on the time to first <italic>D. citri</italic> detection in a grove. Groupings for plots were created based on natural breaks in the independent variables, resulting in two relative values for roads and three for urban intensity.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="finsc-02-783285-g0003.tif"/>
</fig>
<p>For the analysis of grove-level characteristics and neighbor proximity effects on invasion rate, of the dozens of models evaluated, four performed similarly well, with all other model formulations having AIC values at least 15 units greater (<xref ref-type="supplementary-material" rid="SM1">Supplementary Table S1</xref>). Based on the model-averaged coefficients for these preferred models, grove area, grove edge, and grove perforation were significantly positively related to invasion rate (<xref ref-type="fig" rid="F4">Figure 4</xref>). The model-averaged coefficients indicated that both early-phase distance to nearest grove occurrence (i.e., grove early-phase) and late-phase distance to nearest urban occurrence (i.e., urban late-phase) were not significantly related to invasion rate. However, the late-phase distance to nearest grove occurrence (i.e., grove late-phase) and early-phase distance to nearest urban occurrence (i.e., urban early-phase) are strongly negatively related to invasion rate (<xref ref-type="fig" rid="F4">Figure 4</xref>; <xref ref-type="supplementary-material" rid="SM1">Supplementary Figure S1</xref>). With respect to these last effects, the overall pattern of new <italic>D. citri</italic> grove detections suggest that potential spillover (i.e., the invasion kernel) from urban early-phase occurrences is strongest within 25 km and spillover from grove late-phase occurrences is strongest within 10 km (<xref ref-type="fig" rid="F5">Figure 5</xref>).</p>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p>Model-averaged coefficient values (&#x000B1;95% confidence intervals) across the four best models for grove characteristics (grove area [&#x0201C;area&#x0201D;], grove edge [&#x0201C;edge&#x0201D;], grove perforation [&#x0201C;perf&#x0201D;]) and neighbor effects (grove early-phase distance [&#x0201C;grove.e&#x0201D;], grove late-phase [&#x0201C;grove.l&#x0201D;], urban early-phase [urban.e], or urban late-phase [&#x0201C;urban.l&#x0201D;] occurrence of <italic>D. citri</italic>) on time to <italic>D. citri</italic> first detection in a given citrus grove. Variables with error bars overlapping 0 are appropriately viewed as having minimal effects on invasion rate.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="finsc-02-783285-g0004.tif"/>
</fig>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p>Overall proportion of groves with at least one <italic>D. citri</italic> detection as a function of distance (km) from the nearest <bold>(A)</bold> urban early-phase, or <bold>(B)</bold> grove late-phase occurrence of <italic>D. citri</italic>.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="finsc-02-783285-g0005.tif"/>
</fig>
</sec>
<sec>
<title>Cumulative <italic>D. citri</italic> Detections</title>
<p>For the analysis of grove-level and urban landscape-level effects, seven models performed similarly well with respect to AIC rankings (<xref ref-type="supplementary-material" rid="SM1">Supplementary Table S2</xref>). Based on the model-averaged coefficients grove area and time since first detection were strongly and positively related to more <italic>D. citri</italic> detections (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure S2</xref>), distance to roads was slightly negatively related to the number of detections (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure S2C</xref>), and the other three variables had equivocal effects (<xref ref-type="fig" rid="F6">Figure 6</xref>).</p>
<fig id="F6" position="float">
<label>Figure 6</label>
<caption><p>Model-averaged coefficient values (&#x000B1; 95% confidence intervals) across the seven best models for grove characteristics, landscape context (distance to roads [&#x0201C;road&#x0201D;], urbanization intensity [&#x0201C;urb&#x0201D;]), and temporal effects (number of months after first <italic>D. citri</italic> detection in the grove) on cumulative number of <italic>D. citri</italic> detections.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="finsc-02-783285-g0006.tif"/>
</fig>
<p>For the analysis of grove-level characteristics and neighbor effects, of the dozens of models evaluated, one substantially outperformed the others (&#x0003E;10 AIC units lower). That model included significant effects of grove area (&#x003C7;<sup>2</sup> = 55.991, <italic>df</italic> = 1, <italic>p</italic> &#x0003C; 0.0001), time since first detection (&#x003C7;<sup>2</sup> = 255.63, <italic>df</italic> = 1, <italic>p</italic> &#x0003C; 0.0001), urban late-phase distance (&#x003C7;<sup>2</sup> = 23.784, <italic>df</italic> = 1, <italic>p</italic> &#x0003C; 0.0001), and grove late-phase distance (&#x003C7;<sup>2</sup> = 24.470, <italic>df</italic> = 1, <italic>p</italic> &#x0003C; 0.0001). Both grove area (slope &#x000B1; SE = 0.677 &#x000B1; 0.09) and time since first detection (0.634 &#x000B1; 0.04) were positively associated with <italic>D. citri</italic> detections (<xref ref-type="fig" rid="F7">Figures 7A,B</xref>). Conversely, greater urban late-phase distances (&#x02212;1.136 &#x000B1; 0.233) and nearest grove detections (&#x02212;1.737 &#x000B1; 0.351) were associated with fewer cumulative <italic>D. citri</italic> detections (<xref ref-type="fig" rid="F7">Figures 7C,D</xref>).</p>
<fig id="F7" position="float">
<label>Figure 7</label>
<caption><p>Effect of <bold>(A)</bold> grove area, <bold>(B)</bold> time since first <italic>D. citri</italic> detection, and distance to the nearest <bold>(C)</bold> grove late-phase <italic>D. citri</italic> occurrence or <bold>(D)</bold> urban late-phase occurrence, on the cumulative number of <italic>D. citri</italic> detections in a citrus grove. Independent variables were all standardized via z-scores; smaller or more negative values equate to smaller groves, more recent first detections, or shorter distances.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="finsc-02-783285-g0007.tif"/>
</fig>
</sec>
</sec>
<sec sec-type="discussion" id="s4">
<title>Discussion</title>
<p>Cross-habitat spillover of organisms is increasingly gaining recognition, particularly in agroecosystems and urban forests. Many studies have focused on the movement of organisms between agricultural habitats and natural habitats (<xref ref-type="bibr" rid="B13">13</xref>, <xref ref-type="bibr" rid="B15">15</xref>, <xref ref-type="bibr" rid="B17">17</xref>). Unlike these studies, our research addresses the less familiar spillover between urban and agricultural habitats within an anthropogenic-dominated landscape, and its importance for the spread and impact of an invasive insect pest of high economic significance. In our study area, citrus constitutes a substantial fraction of urban tree plantings (<xref ref-type="bibr" rid="B29">29</xref>). <italic>D. citri</italic> was first recorded in these urban micro-habitats (such as house backyards and gardens) of southern California in 2008. Within three years, <italic>D. citri</italic> had invaded commercial citrus groves (<xref ref-type="bibr" rid="B30">30</xref>). Overall, our results highlight the role of grove-level factors in combination with two spatially-explicit factors, namely urban landscape-level context and neighbor effects in modulating urban-agricultural spillover.</p>
<p>At the grove-level, grove area stood out as a common factor affecting both rate of invasion and magnitude of <italic>D. citri</italic> infestation such that large groves were invaded faster and also contained greater cumulative detections of <italic>D. citri</italic>. These results are not surprising given that, as noted earlier, trap number varied grove size. Yet, part of the observed effect may also be attributable to larger plantings of citrus trees acting as larger targets for <italic>D. citri</italic> invasion. A similar &#x02018;magnetic effect&#x02019; has been reported among mass flowering crops wherein pollinators tend to spillover in large numbers during blooming periods (<xref ref-type="bibr" rid="B13">13</xref>, <xref ref-type="bibr" rid="B44">44</xref>). After controlling for grove area, the proportion of citrus along grove edges and the extent of grove perforation only affected the rate of invasion. In both cases, a higher proportion of grove edge and grove perforation increased the rate of invasion. Studies have similarly shown that the amount and ecological nature of the edge, and consequent edge-effects, are often important drivers of spillover from groves to natural habitats (<xref ref-type="bibr" rid="B15">15</xref>). Citrus trees at the outer edge of groves have been shown to contain a higher abundance of <italic>D. citri</italic> compared to trees in the interior, particularly those groves on the exterior edge of large farms (<xref ref-type="bibr" rid="B34">34</xref>). Interestingly, grove edge and grove perforation had no significant effect on the cumulative number of <italic>D. citri</italic> detections in a citrus grove, perhaps indicating that once a grove is invaded local recruitment is an important element of <italic>D. citri</italic> population dynamics. Studies have shown that resource availability and abundance can affect population growth and eventual abundance of organisms that spillover into high quality recipient habitats such as mass flowering crops plants (<xref ref-type="bibr" rid="B45">45</xref>, <xref ref-type="bibr" rid="B46">46</xref>).</p>
<p>Both urban landscape-level variables, distance to major roads and urbanization intensity, had significant effects on the rate of <italic>D. citri</italic> invasion. The distance to roads had an overall negative effect on the rate of spillover during early stages of invasion, but then shifted toward a positive effect in more recent years. Urbanization intensity increased the rate of invasion, implying that citrus groves surrounded by higher levels of urbanization are more quickly invaded by <italic>D. citri</italic>. Taken together, the positive effect of urbanization intensity on rate of invasion suggests heavily disturbed urban-grove edges, along with the distance to roads, can facilitate introduction pathways and hasten <italic>D. citri</italic> invasion. Distance to major roads tends to result in greater <italic>D. citri</italic> abundance within groves, which again suggests the role of roads in serving as introduction pathways and potentially affecting multiple introductions of <italic>D. citri</italic> individuals from the surrounding landscape. An earlier study on the initial <italic>D. citri</italic> detections within urban Los Angeles identified major roads entering the city as a likely introduction pathway (<xref ref-type="bibr" rid="B31">31</xref>). In the only other study known to us, proximity to anthropogenic edges such as roads was found to enhance parasitic attacks of a leaf miner in suburban woodland fragments (<xref ref-type="bibr" rid="B15">15</xref>, <xref ref-type="bibr" rid="B47">47</xref>). Our findings highlight the importance of incorporating urban landscape-level factors in spillover studies.</p>
<p>Studies have shown that cross-habitat spillover can influence community dynamics of the recipient habitat through processes such as predation, parasitism, and disease transmission (<xref ref-type="bibr" rid="B13">13</xref>, <xref ref-type="bibr" rid="B14">14</xref>, <xref ref-type="bibr" rid="B18">18</xref>, <xref ref-type="bibr" rid="B19">19</xref>). While grove-level characteristics set the risk threshold for invasion, neighbor effects suggest that proximity to <italic>D. citri</italic> invaded neighborhoods are the eventual deterministic factor driving the spillover. In this respect, it is important to note that the early-phase distance to nearest urban occurrence (urban early-phase) and late-phase distance to nearest grove occurrence (grove late-phase) are negatively related to the rate of invasion. In other words, there appear to be strong spatial linkages between urban occurrences of <italic>D. citri</italic> and its arrival in nearby groves, and among groves once one of them is invaded. The relatively sporadic presence of citrus host plants in urban micro-habitats can initiate the urban-grove spillover, which suggests the need for more spatially comprehensive management plans that take into consideration the potential role of surrounding landscape.</p>
<p>The strong effects of early-phase urban and late-phase grove invasion implies that the spillover pressure from <italic>D. citri</italic> invasion is most pronounced after spreading substantially in the surrounding landscape. In a recent spillover study, the magnitude of bee spillover from mass flowering citrus groves to nearby woodland patches was similarly found to vary in both space and time, such that spillover was not only higher in woodland patches in landscapes with high cover of mass flowering citrus groves but was particularly higher after the citrus blooming period (<xref ref-type="bibr" rid="B44">44</xref>). The study further demonstrated that the observed spillover dynamics are largely driven by variation in flower cover across the landscape, and hence is considered to be resource-mediated. Our findings point to a similar resource-mediated mechanism behind the spillover but with the striking difference that the movement of <italic>D. citri</italic> individuals were initially from relatively resource-poor urban conditions to resource-rich groves. Overall, the neighbor effects suggest that spillover of <italic>D. citri</italic> is not merely a function of spatial proximity to nearest invaded grove and/or urban location but is also shaped by the temporal trajectory of <italic>D. citri</italic> invasion.</p>
<p>To mitigate the impact of the <italic>D. citri</italic> invasion in Southern California, a strategic plan was implemented that included monitoring for <italic>D. citri</italic> and <italic>C</italic>Las, removal of infected trees, the release of biological control agents, a residential insecticide treatment program, restrictions on the movement of plant material, and community outreach (<xref ref-type="bibr" rid="B27">27</xref>, <xref ref-type="bibr" rid="B48">48</xref>). Despite this, urban-to-agriculture spillover occurred within the first few years of the southern California invasion. This observed early-phase effect of urban invasion suggests <italic>D. citri</italic> spillover into citrus groves was initially driven by <italic>D. citri</italic> from citrus plantings in urban and peri-urban locations in close proximity to groves. Later citrus grove invasions were driven largely by spillover from nearby invaded groves, which suggests that the source of <italic>D. citri</italic> spillover into citrus groves shifts over time. This apparent shift in spillover dynamics has implications for management of <italic>D. citri</italic> and similar insect pests that are not restricted to a single habitat type. Specifically, for areas where the invasion is at an advanced stage, such as in Southern California, attempts to mitigate further <italic>D. citri</italic> spread into commercial citrus will likely benefit more from management of invaded groves than from measures focused in residential areas (<xref ref-type="bibr" rid="B27">27</xref>, <xref ref-type="bibr" rid="B48">48</xref>). Conversely, in areas where <italic>D. citri</italic> is not as well-established, such as the Central Valley of California (<xref ref-type="bibr" rid="B49">49</xref>), given competing demands on resources, our results indicate mitigating future <italic>D. citri</italic> invasion of groves is likely to benefit most from management in urban areas. Meanwhile, attempts to mitigate further spread of <italic>C</italic>Las and huanglongbing may require different strategies. Given that to date all of the &#x0007E;1,500 documented cases of huanglongbing in California have occurred in residential citrus (<xref ref-type="bibr" rid="B33">33</xref>), and within cumulative counts of <italic>D. citri</italic> in groves indicating a continuing influx from urban areas, concentrating vector control and inoculum reduction programs (<xref ref-type="bibr" rid="B27">27</xref>) in urban areas is likely to be more impactful than in commercial citrus. In other words, development of appropriate management strategies requires a recognition of the fluid nature of invasion dynamics and careful consideration of the specific management objectives.</p>
</sec>
<sec sec-type="data-availability" id="s5">
<title>Data Availability Statement</title>
<p>The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec id="s6">
<title>Author Contributions</title>
<p>BB and MD: conceptualization and study design. BB, MD, and ST: statistical analysis and writing. MD: funding acquisition and supervision. GS: data acquisition. All authors contributed to the article and approved the submitted version.</p>
</sec>
<sec sec-type="funding-information" id="s7">
<title>Funding</title>
<p>This project was supported by funding from USDA-APHIS-PPQ to MD.</p>
</sec>
<sec sec-type="COI-statement" id="conf1">
<title>Conflict of Interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s8">
<title>Publisher&#x00027;s Note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec> 
</body>
<back>
<ack><p>Thanks to R. Johnson and R. Dunn for helping with access to the ACP trapping geodatabase and to J. Morse, E. Grafton-Cardwell, R. Stouthamer, T. Schartel, and E. Conlisk for helpful discussion.</p>
</ack>
<sec sec-type="supplementary-material" id="s9">
<title>Supplementary Material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/finsc.2022.783285/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/finsc.2022.783285/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Data_Sheet_1.docx" id="SM1" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<label>1.</label>
<citation citation-type="book"><person-group person-group-type="author"><name><surname>G&#x000F3;mez-Baggethun</surname> <given-names>E</given-names></name> <name><surname>Gren</surname> <given-names>&#x000C5;</given-names></name> <name><surname>Barton</surname> <given-names>DN</given-names></name> <name><surname>Langemeyer</surname> <given-names>J</given-names></name> <name><surname>McPhearson</surname> <given-names>T</given-names></name> <name><surname>O&#x00027;Farrell</surname> <given-names>P</given-names></name> <etal/></person-group>. <article-title>Urban ecosystem services</article-title>. In: <source>Urbanization, Biodiversity and Ecosystem Services: Challenges and Opportunities</source>. <publisher-loc>Dordrecht</publisher-loc>: <publisher-name>Springer Nature.</publisher-name> (<year>2013</year>). <pub-id pub-id-type="doi">10.1007/978-94-007-7088-1_11</pub-id></citation>
</ref>
<ref id="B2">
<label>2.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lovell</surname> <given-names>ST</given-names></name> <name><surname>Taylor</surname> <given-names>JR</given-names></name></person-group>. <article-title>Supplying urban ecosystem services through multifunctional green infrastructure in the United States</article-title>. <source>Landsc Ecol.</source> (<year>2013</year>) <volume>28</volume>:<fpage>1447</fpage>&#x02013;<lpage>63</lpage>. <pub-id pub-id-type="doi">10.1007/s10980-013-9912-y</pub-id></citation>
</ref>
<ref id="B3">
<label>3.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Derkzen</surname> <given-names>ML</given-names></name> <name><surname>van Teeffelen</surname> <given-names>AJ</given-names></name> <name><surname>Verburg</surname> <given-names>PH</given-names></name></person-group>. <article-title>Quantifying urban ecosystem services based on high-resolution data of urban green space: an assessment for Rotterdam, the Netherlands</article-title>. <source>J Appl Ecol.</source> (<year>2015</year>) <volume>52</volume>:<fpage>1020</fpage>&#x02013;<lpage>32</lpage>. <pub-id pub-id-type="doi">10.1111/1365-2664.12469</pub-id></citation>
</ref>
<ref id="B4">
<label>4.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cumming</surname> <given-names>GS</given-names></name> <name><surname>Buerkert</surname> <given-names>A</given-names></name> <name><surname>Hoffmann</surname> <given-names>EM</given-names></name> <name><surname>Schlecht</surname> <given-names>E</given-names></name> <name><surname>von Cramon-Taubadel</surname> <given-names>S</given-names></name> <name><surname>Tscharntke</surname> <given-names>T</given-names></name></person-group>. <article-title>Implications of agricultural transitions and urbanization for ecosystem services</article-title>. <source>Nature.</source> (<year>2014</year>) <volume>515</volume>:<fpage>50</fpage>&#x02013;<lpage>7</lpage>. <pub-id pub-id-type="doi">10.1038/nature13945</pub-id><pub-id pub-id-type="pmid">25373674</pub-id></citation></ref>
<ref id="B5">
<label>5.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Power</surname> <given-names>AG</given-names></name></person-group>. <article-title>Ecosystem services and agriculture: tradeoffs and synergies</article-title>. <source>Philos Trans R Soc B.</source> (<year>2010</year>) <volume>365</volume>:<fpage>2959</fpage>&#x02013;<lpage>71</lpage>. <pub-id pub-id-type="doi">10.1098/rstb.2010.0143</pub-id><pub-id pub-id-type="pmid">20713396</pub-id></citation></ref>
<ref id="B6">
<label>6.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cadotte</surname> <given-names>MW</given-names></name> <name><surname>Yasui</surname> <given-names>SLE</given-names></name> <name><surname>Livingstone</surname> <given-names>S</given-names></name> <name><surname>MacIvor</surname> <given-names>JS</given-names></name></person-group>. <article-title>Are urban systems beneficial, detrimental, or indifferent for biological invasion?</article-title> <source>Biol. Invasions.</source> (<year>2017</year>) <volume>19</volume>:<fpage>3489</fpage>&#x02013;<lpage>503</lpage>. <pub-id pub-id-type="doi">10.1007/s10530-017-1586-y</pub-id></citation>
</ref>
<ref id="B7">
<label>7.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gaertner</surname> <given-names>M</given-names></name> <name><surname>Wilson</surname> <given-names>JR</given-names></name> <name><surname>Cadotte</surname> <given-names>MW</given-names></name> <name><surname>MacIvor</surname> <given-names>JS</given-names></name> <name><surname>Zenni</surname> <given-names>RD</given-names></name> <name><surname>Richardson</surname> <given-names>DM</given-names></name></person-group>. <article-title>Non-native species in urban environments: patterns, processes, impacts and challenges</article-title>. <source>Biol Invasions.</source> (<year>2017</year>) <volume>19</volume>:<fpage>3461</fpage>&#x02013;<lpage>9</lpage>. <pub-id pub-id-type="doi">10.1007/s10530-017-1598-7</pub-id></citation>
</ref>
<ref id="B8">
<label>8.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Grimm</surname> <given-names>NB</given-names></name> <name><surname>Foster</surname> <given-names>D</given-names></name> <name><surname>Groffman</surname> <given-names>P</given-names></name> <name><surname>Grove</surname> <given-names>JM</given-names></name> <name><surname>Hopkinson</surname> <given-names>CS</given-names></name> <name><surname>Nadelhoffer</surname> <given-names>KJ</given-names></name> <etal/></person-group>. <article-title>The changing landscape: ecosystem responses to urbanization and pollution across climatic and societal gradients</article-title>. <source>Front Ecol Environ.</source> (<year>2008</year>) <volume>6</volume>:<fpage>264</fpage>&#x02013;<lpage>72</lpage>. <pub-id pub-id-type="doi">10.1890/070147</pub-id></citation>
</ref>
<ref id="B9">
<label>9.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hui</surname> <given-names>C</given-names></name> <name><surname>Richardson</surname> <given-names>DM</given-names></name> <name><surname>Visser</surname> <given-names>V</given-names></name></person-group>. <article-title>Ranking of invasive spread through urban green areas in the world&#x00027;s 100 most populous cities</article-title>. <source>Biol Invasions.</source> (<year>2017</year>) <volume>19</volume>:<fpage>3527</fpage>&#x02013;<lpage>39</lpage>. <pub-id pub-id-type="doi">10.1007/s10530-017-1584-0</pub-id></citation>
</ref>
<ref id="B10">
<label>10.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Padayachee</surname> <given-names>AL</given-names></name> <name><surname>Irlich</surname> <given-names>UM</given-names></name> <name><surname>Faulkner</surname> <given-names>KT</given-names></name> <name><surname>Gaertner</surname> <given-names>M</given-names></name> <name><surname>Proche&#x0015F;</surname> <given-names>S</given-names></name> <name><surname>Wilson</surname> <given-names>JR</given-names></name> <etal/></person-group>. <article-title>How do invasive species travel to and through urban environments?</article-title> <source>Biol. Invasions.</source> (<year>2017</year>) <volume>19</volume>:<fpage>3557</fpage>&#x02013;<lpage>70</lpage>. <pub-id pub-id-type="doi">10.1007/s10530-017-1596-9</pub-id></citation>
</ref>
<ref id="B11">
<label>11.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Grimm</surname> <given-names>NB</given-names></name> <name><surname>Faeth</surname> <given-names>SH</given-names></name> <name><surname>Golubiewski</surname> <given-names>NE</given-names></name> <name><surname>Redman</surname> <given-names>CL</given-names></name> <name><surname>Wu</surname> <given-names>J</given-names></name> <name><surname>Bai</surname> <given-names>X</given-names></name> <etal/></person-group>. <article-title>Global change and the ecology of cities</article-title>. <source>Science.</source> (<year>2008</year>) <volume>319</volume>:<fpage>756</fpage>&#x02013;<lpage>60</lpage>. <pub-id pub-id-type="doi">10.1126/science.1150195</pub-id><pub-id pub-id-type="pmid">18258902</pub-id></citation></ref>
<ref id="B12">
<label>12.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Groffman</surname> <given-names>PM</given-names></name> <name><surname>Avolio</surname> <given-names>M</given-names></name> <name><surname>Cavender-Bares</surname> <given-names>J</given-names></name> <name><surname>Bettez</surname> <given-names>ND</given-names></name> <name><surname>Grove</surname> <given-names>JM</given-names></name> <name><surname>Hall</surname> <given-names>SJ</given-names></name> <etal/></person-group>. <article-title>Ecological homogenization of residential macrosystems</article-title>. <source>Nat Ecol Evol.</source> (<year>2017</year>) <volume>1</volume>:<fpage>1</fpage>&#x02013;<lpage>3</lpage>. <pub-id pub-id-type="doi">10.1038/s41559-017-0191</pub-id><pub-id pub-id-type="pmid">28812594</pub-id></citation></ref>
<ref id="B13">
<label>13.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Blitzer</surname> <given-names>EJ</given-names></name> <name><surname>Dormann</surname> <given-names>CF</given-names></name> <name><surname>Holzschuh</surname> <given-names>A</given-names></name> <name><surname>Klein</surname> <given-names>AM</given-names></name> <name><surname>Rand</surname> <given-names>TA</given-names></name> <name><surname>Tscharntke</surname> <given-names>T</given-names></name></person-group>. <article-title>Spillover of functionally important organisms between managed and natural habitats</article-title>. <source>Agric Ecosyst Environ.</source> (<year>2012</year>) <volume>146</volume>:<fpage>34</fpage>&#x02013;<lpage>43</lpage>. <pub-id pub-id-type="doi">10.1016/j.agee.2011.09.005</pub-id></citation>
</ref>
<ref id="B14">
<label>14.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tscharntke</surname> <given-names>T</given-names></name> <name><surname>Rand</surname> <given-names>TA</given-names></name> <name><surname>Bianchi</surname> <given-names>FJ</given-names></name></person-group>. <article-title>The landscape context of trophic interactions: insect spillover across the crop&#x02014;noncrop interface</article-title>. In: <source>Annales Zoologici Fennici.</source> Finnish Zoological and Botanical Publishing Board (<year>2005</year>). p. <fpage>421</fpage>&#x02013;<lpage>32</lpage>.</citation>
</ref>
<ref id="B15">
<label>15.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rand</surname> <given-names>TA</given-names></name> <name><surname>Tylianakis</surname> <given-names>JM</given-names></name> <name><surname>Tscharntke</surname> <given-names>T</given-names></name></person-group>. <article-title>Spillover edge effects: the dispersal of agriculturally subsidized insect natural enemies into adjacent natural habitats</article-title>. <source>Ecol Lett.</source> (<year>2006</year>) <volume>9</volume>:<fpage>603</fpage>&#x02013;<lpage>14</lpage>. <pub-id pub-id-type="doi">10.1111/j.1461-0248.2006.00911.x</pub-id><pub-id pub-id-type="pmid">16643305</pub-id></citation></ref>
<ref id="B16">
<label>16.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rand</surname> <given-names>TA</given-names></name> <name><surname>Louda</surname> <given-names>SM</given-names></name></person-group>. <article-title>Spillover of agriculturally subsidized predators as a potential threat to native insect herbivores in fragmented landscapes</article-title>. <source>Conserv Biol.</source> (<year>2006</year>) <volume>20</volume>:<fpage>1720</fpage>&#x02013;<lpage>9</lpage>. <pub-id pub-id-type="doi">10.1111/j.1523-1739.2006.00507.x</pub-id><pub-id pub-id-type="pmid">17181807</pub-id></citation></ref>
<ref id="B17">
<label>17.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Frost</surname> <given-names>CM</given-names></name> <name><surname>Didham</surname> <given-names>RK</given-names></name> <name><surname>Rand</surname> <given-names>TA</given-names></name> <name><surname>Peralta</surname> <given-names>G</given-names></name> <name><surname>Tylianakis</surname> <given-names>JM</given-names></name></person-group>. <article-title>Community-level net spillover of natural enemies from managed to natural forest</article-title>. <source>Ecology.</source> (<year>2015</year>) <volume>96</volume>:<fpage>193</fpage>&#x02013;<lpage>202</lpage>. <pub-id pub-id-type="doi">10.1890/14-0696.1</pub-id><pub-id pub-id-type="pmid">26236904</pub-id></citation></ref>
<ref id="B18">
<label>18.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Power</surname> <given-names>AG</given-names></name> <name><surname>Mitchell</surname> <given-names>CE</given-names></name></person-group>. <article-title>Pathogen spillover in disease epidemics</article-title>. <source>Am Nat.</source> (<year>2004</year>) <volume>164</volume>:<fpage>S79</fpage>&#x02013;<lpage>89</lpage>. <pub-id pub-id-type="doi">10.1086/424610</pub-id><pub-id pub-id-type="pmid">15540144</pub-id></citation></ref>
<ref id="B19">
<label>19.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pereira-Peixoto</surname> <given-names>MH</given-names></name> <name><surname>Pufal</surname> <given-names>G</given-names></name> <name><surname>Martins</surname> <given-names>CF</given-names></name> <name><surname>Klein</surname> <given-names>AM</given-names></name></person-group>. <article-title>Spillover of trap-nesting bees and wasps in an urban&#x02013;rural interface</article-title>. <source>J Insect Conserv.</source> (<year>2014</year>) <volume>18</volume>:<fpage>815</fpage>&#x02013;<lpage>26</lpage>. <pub-id pub-id-type="doi">10.1007/s10841-014-9688-7</pub-id></citation>
</ref>
<ref id="B20">
<label>20.</label>
<citation citation-type="web"><person-group person-group-type="author"><collab>University University of California Riverside</collab></person-group>. <article-title>Center for Invasive Species Research</article-title>. (<year>2021</year>). Available online at: <ext-link ext-link-type="uri" xlink:href="https://cisr.ucr.edu/invasive-species/non-native-invasive-species-california">https://cisr.ucr.edu/invasive-species/non-native-invasive-species-california</ext-link> (accessed September 22, 2021).</citation>
</ref>
<ref id="B21">
<label>21.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dowell</surname> <given-names>RV</given-names></name> <name><surname>Gill</surname> <given-names>RJ</given-names></name> <name><surname>Jeske</surname> <given-names>DR</given-names></name> <name><surname>Hoddle</surname> <given-names>MS</given-names></name></person-group>. <article-title>Exotic terrestrial macro-invertebrate invaders in California from 1700 to 2010: an analysis of records</article-title>. <source>CA Acad Sci Ser 4</source>. (<year>2016</year>) <volume>63</volume>: <fpage>63</fpage>&#x02013;<lpage>157</lpage>.</citation>
</ref>
<ref id="B22">
<label>22.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Grafton-Cardwell</surname> <given-names>EE</given-names></name> <name><surname>Stelinski</surname> <given-names>LL</given-names></name> <name><surname>Stansly</surname> <given-names>PA</given-names></name></person-group>. <article-title>Biology and management of Asian citrus psyllid, vector of the huanglongbing pathogens</article-title>. <source>Annu Rev Entomol.</source> (<year>2013</year>) <volume>58</volume>:<fpage>413</fpage>&#x02013;<lpage>32</lpage>. <pub-id pub-id-type="doi">10.1146/annurev-ento-120811-153542</pub-id><pub-id pub-id-type="pmid">23317046</pub-id></citation></ref>
<ref id="B23">
<label>23.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gottwald</surname> <given-names>TR</given-names></name></person-group>. <article-title>Current epidemiological understanding of Citrus Huanglongbing<sup>&#x0002A;</sup></article-title>. <source>Annu Rev Phytopathol.</source> (<year>2010</year>) <volume>48</volume>:<fpage>119</fpage>&#x02013;<lpage>39</lpage>. <pub-id pub-id-type="doi">10.1146/annurev-phyto-073009-114418</pub-id><pub-id pub-id-type="pmid">20415578</pub-id></citation></ref>
<ref id="B24">
<label>24.</label>
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Hodges</surname> <given-names>AW</given-names></name> <name><surname>Spreen</surname> <given-names>TH</given-names></name></person-group>. <source>Economic Impacts of Citrus Greening (HLB) in Florida, 2006/07&#x02013;2010/11</source>. <publisher-loc>Gainesville, FL</publisher-loc>: <publisher-name>University of Florida, Institute of Food and Agricultural Sciences</publisher-name>. (<year>2012</year>).</citation>
</ref>
<ref id="B25">
<label>25.</label>
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Webber</surname> <given-names>HJ</given-names></name> <name><surname>Barker</surname> <given-names>R</given-names></name> <name><surname>Ferguson</surname> <given-names>L</given-names></name></person-group>. <article-title>History and development of the California citrus industry</article-title>. In: L. Ferguson and E.E. Grafton-Cardwell, editors. <source>Citrus Production Manual</source>. <publisher-loc>Oakland, CA</publisher-loc>: <publisher-name>University of California, Agriculture and Natural Resources</publisher-name>. (<year>2014</year>). p. <fpage>3</fpage>&#x02013;<lpage>20</lpage>.</citation>
</ref>
<ref id="B26">
<label>26.</label>
<citation citation-type="web"><person-group person-group-type="author"><name><surname>Babcock</surname> <given-names>BA</given-names></name></person-group>. <source>Economic Impact of California&#x00027;s Citrus Industry</source>. (<year>2018</year>). Available online at: <ext-link ext-link-type="uri" xlink:href="https://ccqc.org/wp-content/uploads/Economic-Contribution-of-California-Citrus-Industry21_rev_060418.pdf">https://ccqc.org/wp-content/uploads/Economic-Contribution-of-California-Citrus-Industry21_rev_060418.pdf</ext-link></citation>
</ref>
<ref id="B27">
<label>27.</label>
<citation citation-type="web"><person-group person-group-type="author"><collab>California Department of Food Agriculture (CDFA)</collab></person-group>. <source>Action Plan for Asian Citrus Psyllid and Huanglongbing (citrus greening) in California</source>. (<year>2021</year>). Available online at: <ext-link ext-link-type="uri" xlink:href="https://www.cdfa.ca.gov/citruscommittee/docs/ActionPlan.pdf">https://www.cdfa.ca.gov/citruscommittee/docs/ActionPlan.pdf</ext-link></citation>
</ref>
<ref id="B28">
<label>28.</label>
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Sawyer</surname> <given-names>RC</given-names></name></person-group>. <source>To Make a Spotless Orange</source>. <publisher-loc>West Lafayette, IN</publisher-loc>: <publisher-name>Purdue University Press</publisher-name> (<year>2002</year>).<fpage>22</fpage><pub-id pub-id-type="pmid">30313346</pub-id></citation></ref>
<ref id="B29">
<label>29.</label>
<citation citation-type="web"><person-group person-group-type="author"><name><surname>Hoddle</surname> <given-names>M</given-names></name></person-group>. <article-title>Biological control of Asian citrus psyllid takes big steps forward in Southern California</article-title>. <source>CAPCA Adviser.</source> (<year>2012</year>) <volume>15</volume>:<fpage>18</fpage>&#x02013;<lpage>20</lpage>. Available online at: <ext-link ext-link-type="uri" xlink:href="https://ucanr.edu/sites/scienceforcitrushealth/Research_Snapshots/Psyllid_Management/Hoddle/">https://ucanr.edu/sites/scienceforcitrushealth/Research_Snapshots/Psyllid_Management/Hoddle/</ext-link> (accessed September 25, 2021).</citation>
</ref>
<ref id="B30">
<label>30.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bayles</surname> <given-names>BR</given-names></name> <name><surname>Thomas</surname> <given-names>SM</given-names></name> <name><surname>Simmons</surname> <given-names>GS</given-names></name> <name><surname>Grafton-Cardwell</surname> <given-names>EE</given-names></name> <name><surname>Daugherty</surname> <given-names>MP</given-names></name></person-group>. <article-title>Spatiotemporal dynamics of the Southern California Asian Citrus Psyllid (<italic>Diaphorina citri</italic>) invasion</article-title>. <source>PLoS ONE.</source> (<year>2017</year>) <volume>12</volume>:<fpage>e0173226</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0173226</pub-id><pub-id pub-id-type="pmid">28278188</pub-id></citation></ref>
<ref id="B31">
<label>31.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Thomas</surname> <given-names>SM</given-names></name> <name><surname>Simmons</surname> <given-names>GS</given-names></name> <name><surname>Daugherty</surname> <given-names>MP</given-names></name></person-group>. <article-title>Spatiotemporal distribution of an invasive insect in an urban landscape: introduction, establishment and impact</article-title>. <source>Landsc Ecol.</source> (<year>2017</year>) <volume>32</volume>:<fpage>2041</fpage>&#x02013;<lpage>57</lpage>. <pub-id pub-id-type="doi">10.1007/s10980-017-0565-0</pub-id></citation>
</ref>
<ref id="B32">
<label>32.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kumagai</surname> <given-names>LB</given-names></name> <name><surname>LeVesque</surname> <given-names>CS</given-names></name> <name><surname>Blomquist</surname> <given-names>CL</given-names></name> <name><surname>Madishetty</surname> <given-names>K</given-names></name> <name><surname>Guo</surname> <given-names>Y</given-names></name> <name><surname>Woods</surname> <given-names>PW</given-names></name> <etal/></person-group>. <article-title>First report of Candidatus <italic>Liberibacter asiacticus</italic> associated with citrus huanglongbing in California</article-title>. <source>Plant Dis.</source> (<year>2013</year>) <volume>97</volume>:<fpage>283</fpage>. <pub-id pub-id-type="doi">10.1094/PDIS-09-12-0845-PDN</pub-id><pub-id pub-id-type="pmid">30722341</pub-id></citation></ref>
<ref id="B33">
<label>33.</label>
<citation citation-type="web"><person-group person-group-type="author"><collab>Data Analysis and Tactical Operations Cell (DATOC)</collab></person-group>. <article-title>The HLB epidemic in California</article-title>. (<year>2021</year>). Available online at: <ext-link ext-link-type="uri" xlink:href="https://www.datoc.us/the-hlb-epidemic">https://www.datoc.us/the-hlb-epidemic</ext-link> (accessed November 24, 2021).</citation>
</ref>
<ref id="B34">
<label>34.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>S&#x000E9;tamou</surname> <given-names>M</given-names></name> <name><surname>Bartels</surname> <given-names>DW</given-names></name></person-group>. <article-title>Living on the edges: spatial niche occupation of Asian citrus psyllid, <italic>Diaphorina citri</italic> Kuwayama (Hemiptera: Liviidae), in citrus groves</article-title>. <source>PLoS ONE.</source> (<year>2015</year>) <volume>10</volume>:<fpage>e0131917</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0131917</pub-id><pub-id pub-id-type="pmid">26193111</pub-id></citation></ref>
<ref id="B35">
<label>35.</label>
<citation citation-type="web"><person-group person-group-type="author"><name><surname>Parent</surname> <given-names>J</given-names></name> <name><surname>Hurd</surname> <given-names>J</given-names></name></person-group>. <source>Landscape Fragmentation Tool (LFT v2.0). Center for Land Use Education Research</source>. (<year>2010</year>). Available online at: <ext-link ext-link-type="uri" xlink:href="http://clear.uconn.edu/tools/lft/lft2/index.htm">http://clear.uconn.edu/tools/lft/lft2/index.htm</ext-link></citation>
</ref>
<ref id="B36">
<label>36.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fischer</surname> <given-names>J</given-names></name> <name><surname>Lindenmayer</surname> <given-names>DB</given-names></name></person-group>. <article-title>Landscape modification and habitat fragmentation: a synthesis</article-title>. <source>Glob Ecol Biogeogr.</source> (<year>2007</year>) <volume>16</volume>:<fpage>265</fpage>&#x02013;<lpage>80</lpage>. <pub-id pub-id-type="doi">10.1111/j.1466-8238.2007.00287.x</pub-id></citation>
</ref>
<ref id="B37">
<label>37.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gelman</surname> <given-names>A</given-names></name></person-group>. <article-title>Scaling regression inputs by dividing by two standard deviations</article-title>. <source>Stat Med.</source> (<year>2008</year>) <volume>27</volume>:<fpage>2865</fpage>&#x02013;<lpage>73</lpage>. <pub-id pub-id-type="doi">10.1002/sim.3107</pub-id><pub-id pub-id-type="pmid">17960576</pub-id></citation></ref>
<ref id="B38">
<label>38.</label>
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Crawley</surname> <given-names>MJ</given-names></name></person-group>. <source>The R Book</source>. <publisher-loc>Chichester</publisher-loc>: <publisher-name>John Wiley and Sons, Inc.</publisher-name> (<year>2009</year>).</citation>
</ref>
<ref id="B39">
<label>39.</label>
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Hosmer DW</surname> <given-names>Jr</given-names></name> <name><surname>Lemeshow</surname> <given-names>S</given-names></name></person-group>. <source>Applied Survival Analysis: Regression Modeling of Time to Event Data</source>. <publisher-loc>New York, NY</publisher-loc>: <publisher-name>Wiley</publisher-name>. (<year>1999</year>).</citation>
</ref>
<ref id="B40">
<label>40.</label>
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Fox</surname> <given-names>J</given-names></name> <name><surname>Weisberg</surname> <given-names>S</given-names></name></person-group>. <source>An R Companion to Applied Regression</source>. <publisher-loc>Thousand Oaks, CA</publisher-loc>: <publisher-name>Sage Publishing, Inc</publisher-name>. (<year>2018</year>).</citation>
</ref>
<ref id="B41">
<label>41.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Burnham</surname> <given-names>KP</given-names></name> <name><surname>Anderson</surname> <given-names>DR</given-names></name></person-group>. <source>Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach. New York, NY: Springer-Verlag.</source> (<year>2002</year>).</citation>
</ref>
<ref id="B42">
<label>42.</label>
<citation citation-type="web"><person-group person-group-type="author"><name><surname>Barton</surname> <given-names>K.</given-names></name></person-group> <article-title>MuMIn: <italic>Multi-Model Inference</italic></article-title>. <source>R Package vs.1.9.0</source>. (<year>2013</year>). Available online at: <ext-link ext-link-type="uri" xlink:href="http://cran.r-project.org/web/packages/MuMIn/">http://cran.r-project.org/web/packages/MuMIn/</ext-link></citation>
</ref>
<ref id="B43">
<label>43.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wonham</surname> <given-names>MJ</given-names></name> <name><surname>Lewis</surname> <given-names>MA</given-names></name> <name><surname>Renclawowicz</surname> <given-names>J</given-names></name> <name><surname>van den Driessche</surname> <given-names>P</given-names></name></person-group>. <article-title>Transmission assumptions generate conflicting predictions in host-vector disease models: a case study in West Nile virus</article-title>. <source>Ecol Lett.</source> (<year>2000</year>) <volume>9</volume>:<fpage>706</fpage>&#x02013;<lpage>25</lpage>. <pub-id pub-id-type="doi">10.1111/j.1461-0248.2006.00912.x</pub-id><pub-id pub-id-type="pmid">16706915</pub-id></citation></ref>
<ref id="B44">
<label>44.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gonz&#x000E1;lez-Varo</surname> <given-names>JP</given-names></name> <name><surname>Vil&#x000E0;</surname> <given-names>M</given-names></name></person-group>. <article-title>Spillover of managed honeybees from mass-flowering crops into natural habitats</article-title>. <source>Biol Conserv.</source> (<year>2017</year>) <volume>212</volume>:<fpage>376</fpage>&#x02013;<lpage>82</lpage>. <pub-id pub-id-type="doi">10.1016/j.biocon.2017.06.018</pub-id></citation>
</ref>
<ref id="B45">
<label>45.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Westphal</surname> <given-names>C</given-names></name> <name><surname>Steffan-Dewenter</surname> <given-names>I</given-names></name> <name><surname>Tscharntke</surname> <given-names>T</given-names></name></person-group>. <article-title>Mass flowering crops enhance pollinator densities at a landscape scale</article-title>. <source>Ecol Lett.</source> (<year>2003</year>) <volume>6</volume>:<fpage>961</fpage>&#x02013;<lpage>5</lpage>. <pub-id pub-id-type="doi">10.1046/j.1461-0248.2003.00523.x</pub-id></citation>
</ref>
<ref id="B46">
<label>46.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Westphal</surname> <given-names>C</given-names></name> <name><surname>Steffan-Dewenter</surname> <given-names>I</given-names></name> <name><surname>Tscharntke</surname> <given-names>T</given-names></name></person-group>. <article-title>Mass flowering oilseed rape improves early colony growth but not sexual reproduction of bumblebees</article-title>. <source>J Appl Ecol.</source> (<year>2009</year>) <volume>46</volume>:<fpage>187</fpage>&#x02013;<lpage>93</lpage>. <pub-id pub-id-type="doi">10.1111/j.1365-2664.2008.01580.x</pub-id></citation>
</ref>
<ref id="B47">
<label>47.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>McGeoch</surname> <given-names>MA</given-names></name> <name><surname>Gaston</surname> <given-names>KJ</given-names></name></person-group>. <article-title>Edge effects on the prevalence and mortality factors of Phytomyza ilicis (Diptera, Agromyzidae) in a suburban woodland</article-title>. <source>Ecol Lett.</source> (<year>2000</year>) <volume>3</volume>:<fpage>23</fpage>&#x02013;<lpage>9</lpage>. <pub-id pub-id-type="doi">10.1046/j.1461-0248.2000.00114.x</pub-id></citation>
</ref>
<ref id="B48">
<label>48.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hoddle</surname> <given-names>M</given-names></name> <name><surname>Mace</surname> <given-names>K</given-names></name> <name><surname>Steggall</surname> <given-names>J</given-names></name></person-group>. <article-title>Proactive biological control: a cost-effective management option for invasive pests</article-title>. <source>Calif Agric.</source> (<year>2018</year>) <volume>72</volume>:<fpage>148</fpage>&#x02013;<lpage>50</lpage>. <pub-id pub-id-type="doi">10.3733/ca.2018a0027</pub-id></citation>
</ref>
<ref id="B49">
<label>49.</label>
<citation citation-type="web"><person-group person-group-type="author"><collab>University University of California Agriculture Natural Resources.</collab></person-group> (<year>2021</year>). <source>Asian Citrus Psyllid Distribution and Management</source>. Available online at: <ext-link ext-link-type="uri" xlink:href="https://ucanr.edu/sites/ACP/">https://ucanr.edu/sites/ACP/</ext-link> (accessed November 24, 2021).</citation>
</ref>
</ref-list> 
</back>
</article> 