Fitzpatrick, Preisser, Porter, Elkinton, Waller, Carlin, and Ellison adapt Lu and Carlin's Bayesian areal wombling β previously used for disease and rabies mapping β to ecological boundaries, using county-level dates of first reported hemlock woolly adelgid (HWA, Adelges tsugae) infestation across eastern North America . Boundary likelihood values Ξij = |Yi β Yj| between neighboring counties yield posterior probabilities P(Ξij > c | y) via G = 2000 MCMC draws, with 100,000 burn-in plus 100,000 sampling iterations and GelmanβRubin convergence checks .
A validation simulation correctly recovered designed boundaries around York County, PA for both ΞΌ- and Ο-based boundary maps β the strongest internal evidence the method works as intended . Applied to real HWA data, the areal model placed highest boundary probabilities near initial establishment sites, along Appalachian ridges, and at the northern range edge; however, winter temperature and hemlock abundance were not significant, and residual-based maps showed covariates explained few boundaries . The authors' own local edge wombling extension (modeling border-level differences) improved explanatory power β hemlock abundance and population density became significant β an honest self-correction showing global county-mean covariates poorly couple to local barriers .
Reported strengths: direct posterior boundary probabilities (impossible in classical null-based wombling), explicit spatial-dependence handling, and missing-data estimation capability . Limitations the authors acknowledge: CAR smoothing can shift an early-colonized county's modeled infestation date later (asymmetric error, since reported dates cannot be earlier than actual), and spatially correlated errors may distort fixed-effect estimates and residual-based maps . Limitations underemphasized: ad hoc county-level HWA reporting is itself observation-biased (low-density populations are detected late, per the authors' own earlier work), the exact county count is not stated in the text, single threshold choices (60 vs. 36 months) drive map appearance, and no external ground-truth boundary data validate the maps beyond the simulation . The method's lineage traces to Womble's 1951 differential systematics and public-health areal wombling, so the novelty is in ecological transfer rather than statistical invention β a later methodological review confirms Bayesian wombling matured as a distinct framework with both gradient and boundary measures .
What would change this assessment: independent validation against field-surveyed HWA fronts, sensitivity analysis over thresholds c and neighbor definitions (the authors note inverse-distance weighting was untested), and demonstration on a second invasion system. Confidence in the methodological demonstration is high; confidence in the specific ecological conclusions about HWA barriers is moderate given covariate failures and reporting bias.
Know what changed, what holds up, and what remains uncertain. Every Friday. No ads.