Table S13. Occupancy model selection results. ψ[S] = initial occupancy probability; γ[U,S] = probability of colonization if site was unsuitable the previous year; γ[S,S] = probability of colonization if site was suitable the previous year; ε = probability of extinction if site was suitable the previous year; π[S] = initial probability of suitable habitat; η[S,S] = probability of habitat remaining suitable if site was suitable the previous year. η [U,S] = probability of habitat becoming suitable if site was unsuitable the previous year. p = probability of species detection. TRIB and URB reflect site groupings by tributary and urbanization category, respectively. Time variation is indicated by t while parameters that were modeled as constant across years are indicated by a dot, “(.)”. TS = temperature stability; WD = water depth; RC = rock cover; MH = presence of maidenhair ferns; CA = calcium carbonate deposition. Model ψ[S] γ[U,S] γ[S,S] *(TSt + WDt + MH), ε(t), π[S](URB), η[S,S](t), η [U,S](t), p(TRIB) ψ
γ
γ
(TSt + WDt + CA + MH), ε(t), π (URB), η
[S] [U,S] [S,S] *
[S]
[S,S]
(t), η
[U,S]
(t), p(TRIB)
ψ[S] γ[U,S] γ[S,S] *(TSt + WDt), ε(t), π[S](URB), η[S,S](t), η [U,S](t), p(TRIB) ψ[S] γ[U,S] γ[S,S] *(CA + TSt + WDt), ε(t), π[S](URB), η[S,S](t), η [U,S](t), p(TRIB) ψ[S] γ[U,S] γ[S,S] *(RCt + TSt + WDt), ε(t), π[S](URB), η[S,S](t), η [U,S](t), p(TRIB) ψ[S] γ[U,S] γ[S,S]*(TSt), ε(t), π[S](URB), η[S,S](t), η [U,S](t), p(TRIB) ψ[S] γ[U,S] γ[S,S] *(RCt), ε(t), π[S](URB), η[S,S](t), η [U,S](t), p(TRIB) ψ[S] γ[U,S] γ[S,S] *(CA + MH), ε(t), π[S](URB), η[S,S](t), η [U,S](t), p(TRIB) ψ[S] γ[U,S] γ[S,S] *(RCt + TSt), ε(t), π[S](URB), η[S,S](t), η [U,S](t), p(TRIB) ψ[S] (TRIB) γ[U,S](t), γ[S,S] (t), ε(t), π[S](TRIB), η[S,S](t), η [U,S](t), p(TRIB) ψ[S] γ[U,S] γ[S,S] *(MH), ε(t), π[S](URB), η[S,S](t), η [U,S](t), p(TRIB) ψ[S] γ[U,S] γ[S,S] *(CA), ε(t), π[S](URB), η[S,S](t), η [U,S](t), p(TRIB) ψ[S] γ[U,S] γ[S,S] *(WDt), ε(t), π[S](URB), η[S,S](t), η [U,S](t), p(TRIB) ψ[S] γ[U,S] γ[S,S] *(RCt + WDt), ε(t), π[S](URB), η[S,S](t), η [U,S](t), p(TRIB) ψ[S] (URB) γ[U,S](t), γ[S,S] (t), ε(t), π[S](URB), η[S,S](t), η [U,S](t), p(TRIB) ψ[S], γ[U,S](t), γ[S,S](t), ε(t), π[S](URB), η[S,S](t), η [U,S](t), p(TRIB) ψ[S], γ[U,S](t), γ[S,S](t), ε(t), π[S](TRIB), η[S,S](t), η [U,S](t), p(TRIB) ψ[S], γ[U,S]=γ[S,S](t), ε(t), π[S](URB), η[S,S](t), η [U,S](t), p(TRIB) ψ[S], γ[U,S](t), γ[S,S](t), ε (t), π[S](.), η[S,S](t), η [U,S](t), p(TRIB) ψ[S], γ[U,S](t), γ[S,S](t), ε (t), π[S](.) η[S,S](t) η [U,S](t) p(.) ψ[S], γ[U,S](t), γ[S,S](t), ε (t), π[S](.), η[S,S](t), η [U,S](t), p(t) ψ[S], γ[U,S](t), γ[S,S](t), ε (t), π[S](.), η[S,S](t), η [U,S](t), p(TRIB + t)
AICc
ΔAICc
wi
Model Likelihood
K
-2l
753.76 753.87 762.62 766.14 769.74 773.3 777.98 778.15 779.38 779.58 780.64 783.53 784.89 785.45 785.88 787.87 792.7 794.58 794.75 795.82 799.44 800.23
0 0.11 8.86 12.38 15.98 19.54 24.22 24.39 25.62 25.82 26.88 29.77 31.13 31.69 32.12 34.11 38.94 40.82 40.99 42.06 45.68 46.47
0.51 0.4827 0.0061 0.001 0.0002 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
1 0.947 0.0119 0.002 0.0003 0.0001 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
24 26 22 24 25 19 19 20 22 22 18 18 19 22 19 18 21 16 17 13 15 20
703.74 699.5 716.92 716.12 717.55 734.03 738.71 736.75 733.68 733.88 743.5 746.39 745.62 739.75 746.61 750.73 749.15 761.68 759.73 769.22 768.65 758.83
Table S13 continued. Model ψ[S], γ[U,S](t), γ[S,S](t), ε (t), π[S](.), η[S,S](.) η [U,S](.) p(TRIB) [S]
ψ ,γ
[U,S]
(t), γ
(t), ε (t), π (.), η
ψ ,γ
[U,S]
(t), γ
(t), ε (t), π (.), η
[S]
ψ ,γ [S]
[S,S]
[U,S]
[S,S]
=γ
[S]
[S]
[S,S]
[S]
(t), p(TRIB)
=η
[U,S]
[S,S]
(t) η [U,S](t) p(TRIB)
=(1-ε)(t), π (URB) η
[S,S]
=η
[S,S]
[U,S]
(.), p(TRIB)
AICc
ΔAICc
wi
Model Likelihood
K
-2l
828.51 829.32 859.91 1303.23
74.75 75.56 106.15 549.47
0 0 0 0
0 0 0 0
15 15 14 14
797.72 798.53 831.22 1274.54
*Indicates covariates indicated are the same for both initial occupancy and colonization parameters. The effect of covariates on colonization was assumed to be independent of the prior habitat state. Effective sample size was assumed to be 620.