Appendix A. Sampling localities, equilibrium simulations, and simulations with varying r.
Site Name |
Latitude |
Longitude |
Sampling |
Years |
Total n |
Mean prevalence |
2 Mile Creek, 48 |
-14.4289 |
132.3094 |
1 |
2007 |
7 |
100.0 |
46 |
-15.4311 |
131.4753 |
1 |
2007 |
2 |
100.0 |
47 |
-15.4922 |
131.3586 |
1 |
2007 |
2 |
100.0 |
51 |
-13.2786 |
131.1644 |
1 |
2007 |
10 |
0.0 |
53 |
-14.3349 |
132.1220 |
1 |
2007 |
1 |
100.0 |
55 |
-14.0655 |
131.9728 |
1 |
2007 |
2 |
50.0 |
57 |
-13.7295 |
131.7179 |
1 |
2007 |
5 |
0.0 |
Abergowrie |
-18.4667 |
145.8833 |
1 |
1989 |
17 |
5.9 |
Adelaide River Township, ART |
-13.1770 |
131.1700 |
1 |
2008 |
21 |
33.3 |
Bam Bam Springs, HSBBSHS |
-13.3878 |
131.5073 |
1 |
2008 |
10 |
20.0 |
Beerboom Crossing, BBX |
-13.8627 |
131.0754 |
1 |
2008 |
16 |
37.5 |
Bentley |
-19.2610 |
146.8260 |
10 |
19901992 |
186 |
73.1 |
Bloomfield |
-15.8500 |
145.3167 |
2 |
1989 |
52 |
1.9 |
Borroloola |
-16.0680 |
136.3070 |
1 |
2006 |
31 |
87.1 |
Boyne island |
-23.9333 |
151.3500 |
2 |
19891991 |
26 |
0.0 |
Brisbane |
-27.4667 |
153.0167 |
1 |
1989 |
17 |
0.0 |
Bull Ck, O4 |
-13.8835 |
131.2777 |
1 |
2008 |
5 |
60.0 |
Cairns |
-16.9220 |
145.7760 |
1 |
2006 |
30 |
80.0 |
Calvert Hills |
-17.2167 |
137.3167 |
1 |
1989 |
72 |
0.0 |
CapeTribulation |
-16.0167 |
145.4333 |
2 |
1990 |
13 |
0.0 |
CapeWeymouth |
-12.6000 |
142.4333 |
2 |
1989 |
42 |
0.0 |
ChartersTowers |
-20.0667 |
146.2500 |
1 |
1989 |
19 |
0.0 |
Closest swamp to Ooloo Crossing, O5 |
-14.0594 |
131.2565 |
1 |
2008 |
10 |
70.0 |
Coen |
-13.9333 |
143.1833 |
1 |
1989 |
4 |
0.0 |
Copperfield Dam |
-14.0400 |
131.9380 |
3 |
2008 |
31 |
45.2 |
Corroborree |
-12.7692 |
131.4706 |
1 |
2008 |
16 |
37.5 |
Creek Crossing, Ooloo Rd, O6 |
-14.0058 |
131.0799 |
1 |
2008 |
10 |
60.0 |
Cullen, 56 |
-14.0040 |
131.9350 |
2 |
20072008 |
8 |
0.0 |
Darwin |
-12.4620 |
130.8420 |
1 |
2008 |
33 |
0.0 |
Depot Ck |
-13.6479 |
131.5476 |
1 |
2008 |
14 |
21.4 |
Dorat Road drainage line, DR3 |
-13.4808 |
132.4367 |
1 |
2008 |
9 |
0.0 |
Edith Falls |
-14.1830 |
132.1830 |
1 |
2008 |
1 |
0.0 |
Edith Falls Rd |
-14.1830 |
132.1830 |
1 |
2008 |
3 |
100.0 |
Elsey Creek, 24 |
-15.0107 |
133.3537 |
1 |
2007 |
5 |
100.0 |
First Ck inside Hot Springs National Park, O2 |
-13.7571 |
131.4331 |
1 |
2008 |
10 |
0.0 |
Flora River Rd, 43 |
-14.8519 |
131.5778 |
1 |
2007 |
6 |
50.0 |
From Edith Falls to Cullen |
-14.0935 |
132.0590 |
1 |
2008 |
2 |
0.0 |
Greenant Ck, DR2 |
-13.5400 |
131.2493 |
1 |
2008 |
6 |
16.7 |
Hayes Ck Caravan Park, HCCP |
-13.5838 |
131.4600 |
1 |
2008 |
12 |
33.3 |
Hayes Ck Pond, O1 |
-13.6601 |
131.4002 |
1 |
2008 |
12 |
0.0 |
Hayes Creek, 50 |
-13.5692 |
131.4817 |
2 |
20072008 |
6 |
0.0 |
Humpty Doo |
-12.5840 |
131.1260 |
1 |
2008 |
25 |
0.0 |
Katherine Gorge campground |
-14.3163 |
132.4211 |
1 |
2007 |
3 |
0.0 |
katherine Town |
-14.4667 |
132.2667 |
1 |
2007 |
2 |
50.0 |
King River, 128B |
-14.6219 |
132.5956 |
1 |
2007 |
5 |
100.0 |
Leaning Tree Lagoon |
-12.7140 |
131.4200 |
1 |
2008 |
26 |
23.1 |
Limestone Creek dam, 41 |
-14.8122 |
131.9183 |
1 |
2007 |
5 |
80.0 |
Mackay |
-21.1333 |
149.1833 |
2 |
1991 |
25 |
24.0 |
Manton Dam, 52 |
-12.8769 |
131.1408 |
1 |
2007 |
5 |
0.0 |
Mareeba |
-16.9833 |
145.5167 |
1 |
1991 |
8 |
25.0 |
Mary R Station Acess Road Dam, KH2 |
-13.6747 |
131.9734 |
1 |
2008 |
10 |
10.0 |
Mary River Roadhouse, KH3 |
-13.6037 |
132.2190 |
1 |
2008 |
10 |
90.0 |
Mary River Statio HOmestead, MRS |
-13.4863 |
132.0195 |
1 |
2008 |
10 |
70.0 |
Mataranka |
-14.9167 |
133.0667 |
1 |
2007 |
1 |
100.0 |
Middle Ck (Liddy Bridge), O7 |
-13.8098 |
131.3419 |
1 |
2008 |
5 |
80.0 |
Middle Pt |
-12.5990 |
131.3260 |
23 |
20052008 |
160 |
5.6 |
Mt View Rd |
-19.2680 |
146.7960 |
4 |
19891990 |
92 |
2.2 |
Normanton |
-17.6700 |
141.0790 |
1 |
2006 |
30 |
93.3 |
Palmerston |
-12.4810 |
130.9840 |
1 |
2008 |
30 |
0.0 |
Paluma |
-19.000 |
146.2000 |
2 |
19891990 |
2 |
0.0 |
Pine Ck, 49 |
-13.8186 |
131.8261 |
5 |
20072008 |
23 |
30.4 |
Pond on edge of highway, SH1 |
-13.7294 |
131.7179 |
1 |
2008 |
11 |
9.1 |
Port Douglas |
-16.4833 |
145.4500 |
1 |
1991 |
11 |
100.0 |
Sayle Creek, 21c |
-14.8928 |
133.8664 |
1 |
2007 |
5 |
40.0 |
Scott Creek, 42 |
-14.9258 |
131.8764 |
1 |
2007 |
5 |
60.0 |
Dorat Road |
-13.4809 |
131.1780 |
1 |
2008 |
8 |
25.0 |
Edith Falls turn-off, 54 |
-14.1836 |
132.0358 |
1 |
2007 |
2 |
100.0 |
Timber Creek |
-15.6460 |
130.4770 |
2 |
20062008 |
25 |
8.0 |
Tipperary Homestead, TIPPY |
-13.7348 |
131.0437 |
1 |
2008 |
15 |
93.3 |
Townsville |
-19.2600 |
146.8170 |
2 |
2006 |
31 |
54.8 |
Tully |
-17.9310 |
145.9240 |
1 |
2006 |
17 |
64.7 |
Yungaburra |
-17.2667 |
145.5833 |
1 |
1989 |
10 |
0.0 |
|
FIG. A1. At equilibrium, the choice of cell-size and carrying capacity can affect parasite persistence. To test our model and to provide baseline simulations for comparison with scenarios entailing range-expansion, we ran numerous simulations with the population at spatial equilibrium. We ran equilibrium simulations (with ten replicates each over fifty generations) for all combinations of global carrying capacity Kglobal= 8000, 16000, 24000, 32000; grid-cell size = 2, 4, 5, 10 units; and global (frequency independent) transmission rates T = 0.1, 0.2, …, 0.9. Ideally, we wanted a grid-cell size that was small relative to the scale of dispersal, but large enough to accommodate a sufficient number of individuals that parasite prevalence was maintained at the transmission rate at equilibrium. Large carrying-capacity simulations also carry a computational cost (they are slow to run), so we searched for the smallest grid-cell size that gave us stable parasite populations over a wide range of transmission rates at equlibrium, but which was computationally tractable. At a constant dispersal scale, the cell size and carrying capacity we chose had a strong effect on parasite persistence. Carrying capacities below 32000 tended to result in parasite extinction, particularly at low transmission rates, unless cell size was large (ten units). This figure shows observed prevalence versus transmission rate at equilibrium for global carrying capacity K=32000. Cell sizes vary as per the legend. Lines are drawn through mean values from ten simulations. We chose this carrying capacity and a cell size of five units for all following simulations, because this combination allowed parasite persistence across a wide range of transmission rates at spatial equilibrium, but had a cell size only one-third of the maximum dispersal distance. |
|
FIG. A2. The lag between host and parasite with density dependent transmission under different values of r and parasite cost. Note, for r = 3, the spatial extent was increased to 300, as the lag became greater than the spatial extent (200) used elsewhere. In all panels error bars (1 x standard deviation across 20 replicate runs) are provided for only one of the series. Truncated series in the lower panels occur where the lag exceeded the spatial extent of simulation. |
|
FIG. A3. The lag between host and parasite with density independent transmission under different values of r and parasite cost. In each panel, results are from twenty replicate simulations at each level and error bars represent ± 1 x standard deviation across replicate simulations. |