Figures & data
Figure 1. Bifurcation diagram of the coexistence state for different system parameters. The black line indicates the sign change of the eigenvalues.
![Figure 1. Bifurcation diagram of the coexistence state for different system parameters. The black line indicates the sign change of the eigenvalues.](/cms/asset/fcdf8ec0-d73b-4e3d-98bb-1f328e301bd1/tizo_a_2134477_f0001_oc.jpg)
Figure 2. NOAA Climate.gov graph of annual surface temperature from 1880 to 2020, based on National Centers for Environmental Information (NCEI) data (NOAA National Centers for Environmental information Citation2022). Using the Matlab curve fitting toolbox’s linear approximation, the thick red line indicates a linear approximation to field data determined between the relevant years.
![Figure 2. NOAA Climate.gov graph of annual surface temperature from 1880 to 2020, based on National Centers for Environmental Information (NCEI) data (NOAA National Centers for Environmental information Citation2022). Using the Matlab curve fitting toolbox’s linear approximation, the thick red line indicates a linear approximation to field data determined between the relevant years.](/cms/asset/3b7f4261-3296-457d-bc65-f7542487d290/tizo_a_2134477_f0002_oc.jpg)
Figure 3. (a) Sketch of brood temperature versus time (red), where the lower and upper temperature are assumed to be 32 and 36°C, respectively; (b) sketch of SHB growth (green) over time.
![Figure 3. (a) Sketch of brood temperature versus time (red), where the lower and upper temperature are assumed to be 32 and 36°C, respectively; (b) sketch of SHB growth (green) over time.](/cms/asset/ced70b65-1e5c-4e28-884c-e5ba3ea7d171/tizo_a_2134477_f0003_oc.jpg)
Figure 4. Honey bee population (blue) and SHB population (black) against time obtained for given temperature function as in ). The rest of the parameters are the same as in the text. The figure shows the effect of the maximum developmental rate of SHB on the system dynamics for (a) , (b)
, (c)
and (d)
.
![Figure 4. Honey bee population (blue) and SHB population (black) against time obtained for given temperature function as in Figure 2(a). The rest of the parameters are the same as in the text. The figure shows the effect of the maximum developmental rate of SHB on the system dynamics for (a) rm=0.01, (b) rm=0.02, (c) rm=0.03 and (d) rm=0.08..](/cms/asset/a22f02c5-e775-4a0c-bfc0-9d85ebda97c6/tizo_a_2134477_f0004_oc.jpg)
Figure 5. Honey bee population (blue) and SHB population (black) against time obtained for a given temperature function, as in ). The rest of the parameters are the same as in the text. The effect of the maximum developmental rate of SHB on the system dynamics is fixed and the different SHB invasion rates γ are compared as (a) (b)
(c)
and (d)
![Figure 5. Honey bee population (blue) and SHB population (black) against time obtained for a given temperature function, as in Figure 3(a). The rest of the parameters are the same as in the text. The effect of the maximum developmental rate of SHB on the system dynamics is fixed and the different SHB invasion rates γ are compared as (a) γ=10−8 (b) γ=10−7, (c) γ=10−6 and (d) γ=10−5.](/cms/asset/df9b0cf9-d1ab-434b-a68c-ccb07bceae06/tizo_a_2134477_f0005_oc.jpg)
Figure 6. Honey bee population (blue) and SHB population (black) against time obtained for a given temperature function, as in ). The rest of the parameters are the same as in the text. The effect of the SHB death rate increases with an increase in temperature as (a) (b)
(c)
and (d)
![Figure 6. Honey bee population (blue) and SHB population (black) against time obtained for a given temperature function, as in Figure 2(a). The rest of the parameters are the same as in the text. The effect of the SHB death rate increases with an increase in temperature as (a) δ2=0.01, (b) δ2=0.02 , (c) δ2=0.035 and (d) δ2=0.08 .](/cms/asset/a134bbbf-603f-4987-b6de-2e629b3cd9bd/tizo_a_2134477_f0006_oc.jpg)
Data availability statement
The data that support the findings reported in this study are available from the corresponding author upon reasonable request.