Research Article |
Corresponding author: Matthew L. Bowser ( matthewnkim@gmail.com ) Academic editor: Jianghua Sun
© 2016 Deanna Marie Saltmarsh, Matthew L. Bowser, John M. Morton, Shirley Lang, Daniel Shain, Roman Dial.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Saltmarsh DM, Bowser ML, Morton JM, Sirley Lang S, Shain D, Dial R (2016) Distribution and abundance of exotic earthworms within a boreal forest system in southcentral Alaska. NeoBiota 28: 67-86. doi: 10.3897/neobiota.28.5503
|
![]() |
Little is known about exotic earthworms (Oligochaeta: Lumbricidae) in Alaska outside its southeastern panhandle. This study documents the distribution of exotic earthworms in the relatively undisturbed Kenai National Wildlife Refuge (KNWR), a large, primarily wilderness refuge in southcentral Alaska. We sampled 69 sites near boat launches, along road corridors, and in low human impact areas > 5 km from the road, finding three species of earthworms (Dendrobaena octaedra, Dendrodrilus rubidus, and Lumbricus terrestris). Most road sites (90%) and boat launches (80%) contained earthworms; half (50%) of low human impact sites contained earthworms. Distance to roads was the only significant factor in predicting earthworm occurrence; soil pH, soil moisture, leaf litter depth, and vegetation cover were not. The disparate distributions of these three species suggest that within the KNWR road construction and vehicle traffic played a role in dispersal of the widespread, abundant Dendrobaena octaedra and uncommon Dendrodrilus rubidus; bait abandonment appeared to be the primary method of introduction of Lumbricus terrestris. While the distribution of harmful anecic earthworms in KNWR is currently limited, the prohibition of Lumbricus spp. as bait within conservation units in Alaska may be warranted.
Lumbricidae, earthworm invasion, taiga, bait abandonment, non-native species
Pleistocene glaciations extirpated native earthworms from much of North America, leaving landscapes devoid of earthworms until the introduction of exotic earthworms (Oligochaeta; Lumbricidae) during European settlement (
Material transport by anecic worms, their large adult size, and dense populations have led to substantial ecosystem changes in some parts of North America (
Most studies of exotic earthworms have occurred in temperate regions (
As for most invasive species, human activities, particularly road construction and unintentional transport, likely increase the rate of spread for exotic earthworms above their natural dispersal rate of 5-10 meters a year (
Seventeen species of earthworms are known to occur in Alaska (see records in
Factors such as pH and temperature likely limit earthworm distribution, especially in boreal regions like Alaska (
A rapidly warming climate in Alaska is likely improving environmental conditions for earthworms. Wetlands in Alaska are warming and drying (
The present study documents species composition, distribution, and habitat correlates for earthworms in the Kenai National Wildlife Refuge, a conservation area in southcentral Alaska. A secondary goal is to examine relationships between earthworm occurrence and distance from human-disturbed areas, such as roads and popular fishing areas. The final goal is to identify factors potentially limiting earthworm distribution, such as pH and soil moisture, which are likely to change as the climate continues to warm on the Kenai Peninsula.
Located on the Kenai Peninsula in southcentral Alaska, USA (60°N, 150°W), the Kenai National Wildlife Refuge (KNWR) covers 777,000 ha. Mountains and glaciers characterize the southeastern KNWR (Figure
While most of the KNWR is currently managed as congressionally designated Wilderness, over 130 historic cabins have been inventoried on the Refuge along with other historic resources associated with mining, trapping, oil development, and road construction, mostly in the north of the refuge. Commercial mining and fishing in the area occurred in the late 1800s and early 1900s. Oil exploration began in the northern part of KNWR in the Swanson River area in the mid-1950s, resulting in 2,900 km of bulldozed seismic lines. There are also 240 km of utility and transmission lines and 180 km of established trails within the Refuge (U.S. Fish & Wildlife Service 2010). Fire and associated suppression activities have also been prevalent within the refuge in the past. Major fires of unknown origin occurred in 1871, 1883, 1891, and 1910. Two large, human-caused fires (1947: 125,000 ha and 1969: 35,000 ha) resulted in replacement of mature spruce forests by a mosaic of young mixed conifer-deciduous forest in various stages of succession. Many historic disturbances provided opportunities for earthworm introduction and establishment on the KNWR.
Earthworms were sampled throughout the KNWR during July and August 2011 at 69 total sites representing three levels of human impact. These levels of human impact were characterized by explicit vectors of introduction: (1) boat launches (n = 20), (2) road corridors (n = 20); and (3) low impact areas (> 600 m from any road or facility and 50 m from any trail or river; n = 29). The sampling site locations were chosen within a GIS (ArcGIS v.10.1) by first identifying suitable areas for each impact level and then randomly selecting sample sites. Three 0.09 m2 plots were established at each site to sample earthworm occurrence. At road and boat launch sites, plots were placed two, three, and five meters from the road or edge of a boat launch. The plots were located approximately 10 m apart when possible, establishing a wider area for detecting earthworm occurrence. This protocol was modified at four of the low impact sites accessed by float plane, where only a single 0.09 m2 quadrat was sampled at each site.
We sampled earthworm abundance at each plot using a 30 cm × 30 cm quadrat. Within each quadrat we removed and hand-sorted surface organic material for earthworms. We extracted additional earthworms with a liquid mustard solution of 40 g ground mustard powder in 3.8 L water (
We identified earthworm specimens to species level when possible based on visual observations of external morphology (
We measured leaf litter depth by clearing a small area and measuring the vertical depth of the leaf layer with a 30 cm ruler within each plot. Soil pH, (Soil pH Meter, HANNA, RI), and moisture (Digital Moisture Meter, General, NY) were measured in the field. In each plot, we estimated the percentage cover of litter, grass, forbs, moss, and lichen. General forest type (deciduous, mixed, conifer) of the area was determined from field observations and a GIS land cover layer.
We estimated ash-free dry biomass (g) from the length (mm) of each preserved specimen with the allometric equation of Hale (
g= e2.2853ln(mm)–11.9047
Hale (2004) found that the allometric equations for Octolasion tyrtaeum, Lumbricus spp., and Dendrobaena octaedra were not significantly different from one another, allowing one equation for all species.
Anecic Lumbricus species are potentially more damaging than other genera (
Individual plots served as replicates for each site, but we used site level data for most analyses by averaging plot level data. A site was considered to contain earthworms if individuals were detected in at least one plot. Soil moisture and leaf litter depth were square-root transformed for normality. We calculated remoteness for each site as a measure of distance from the nearest road. This distance, Y (m), was calculated in GIS by using true surface distance, as it accounted for elevation changes and also masked out lakes. To approximate a normal distribution for analyses, we transformed this distance using ln (Y + 1).
All analyses were performed using R, version 3.1.2 (
Prior to occupancy modeling, principal components analysis (PCA) was used to reduce the dimensionality of the eight habitat variables soil pH, soil moisture, leaf litter depth, and moss, grass, lichen, litter, and forb cover percentages. PCA was performed using the PCA function from the FactoMineR package (
Results of the PCA were included in occupancy models (
We used MANOVA to determine the effects of the impact level and earthworm occurrence on the three dependent variables of soil pH, soil moisture, and leaf litter depth.
We found three exotic earthworm species, Dendrobaena octaedra, Dendrodrilus rubidus, and Lumbricus terrestris, on KNWR. We failed to detect Lumbricus rubellus, known from one location on KNWR (http://arctos.database.museum/guid/KNWR:Ento:6755), even though one of our sampling sites was only ~ 48 meters from this known occurrence. Specimen records are included in Suppl. material
The six individuals that were genetically identified using the COI gene showed > 96% probability of identity to their respective species based on a BLAST results. BOLD process ID’s for sequenced specimens are also included in Suppl. material
No site had more than two species confirmed as present. Only four sites (three boat launch sites one road site) contained two species of earthworms, while the majority of sites contained only one species. Dendrobaena octaedra was the most widespread, occurring at 48 (70%) of 69 sites. Most immature earthworms appeared to be D. octaedra based on morphology. Dendrodrilus rubidus occurred at two sites geographically distant from one another. Lumbricus terrestris occurred at three sites, all of which were boat launches located adjacent to one another in the northern part of the Refuge (Figure
Overall, earthworms occurred at 49 (71%) of the 69 sampled sites. Nearly all road sites (18 of 20 total sites, 90%) had earthworms in at least one plot, as did most boat launches (17 of 20 total sites, 85%). In contrast, only half (14 of 29 total sites, 48%) of the low-impact sites contained earthworms.
Earthworms occurred more frequently at roads and boat launch sites than expected, but much less frequently at low impact sites than expected (Table
Observed and expected values for earthworm occurrence in boat launch, road, and low impact sites from a chi-square test of independence.
Boat launch | Road | Low impact | |
---|---|---|---|
Earthworms present | 16 (13.9) | 18 (13.9) | 14 (20.2) |
Earthworms absent | 4 (6.1) | 2 (6.1) | 15 (8.8) |
Observed and expected values for earthworm occurrence at sites in conifer forests and deciduous trees/shrubs from a chi-square test of independence.
Conifer | Deciduous | |
---|---|---|
Earthworms present | 7 (13.3) | 41 (34.7) |
Earthworms absent | 11(4.7) | 6 (12.3) |
Where earthworms were present, the mean density of earthworms was (28.1 ± 4.4 individuals/m2), with mean densities ranging from 23.9 ± 4.5 at road sites to 33.1 ± 6.1 at boat launches (Table
Mean densities (individuals extracted/m2 ± SE) by species and human impact level.
Species | Boat launch | Road | Low impact | Total density |
---|---|---|---|---|
Dendrobaena octaedra | 26.4 ± 6.4 | 23.9 ± 4.5 | 26.9 ± 11.9 | 25.6 ± 4.4 |
Dendrodrilus rubidus | 3.7 ± 3.7 | 7.4 ± 7.4 | - | 5.6 ± 1.9 |
Lumbricus terrestris | 7.4 ± 7.4 | - | - | 7.4 ± 7.4 |
Unidentified immatures | 43.2 ± 3.2 | - | - | 43.2 ± 3.2 |
Total density | 33.1 ± 6.1 | 24.3 ± 4.7 | 26.9 ± 11.9 | 28.1 ± 4.4 |
Mean biomass (mg/m2 ± SE) of earthworms by species and human impact level.
Species | Boat launch | Road | Low impact | Total |
---|---|---|---|---|
Dendrobaena octaedra | 105 ± 30 | 114 ± 19 | 193 ± 61.4 | 135 ± 23 |
Dendrodrilus rubidus | 32.4 | 52.8 | - | 42.6 ± 10.2 |
Lumbricus terrestris | 5651 | - | - | 5651 |
Unidentified immatures | 1891 | - | - | 1891 |
All species | 652 ± 353 | 105 ± 36 | 193 ± 61.4 | 361 ± 144 |
Where earthworms were found, ash-free dry biomass (AFD) of earthworms showed moderate variation (0.36 ± 0.14 AFD g/m2, n = 49; Table
There were many significant correlations among the habitat variables (Table
Correlation matrix for variables used in principal component analysis. *Correlation is significant at the 0.05 level (2-tailed).**Correlation is significant at the 0.01 level (2-tailed).
litter cover | moss cover | grass cover | forbs cover | lichen cover | soil pH | soil moisture | |
---|---|---|---|---|---|---|---|
moss cover | -0.39** | ||||||
grass cover | -0.54** | -0.45** | |||||
forbs cover | -0.12 | 0.24* | -0.28* | ||||
lichen cover | -0.22 | 0.51** | -0.28* | 0.01 | |||
soil pH | 0.15 | -0.62** | 0.51** | -0.33** | -0.47** | ||
soil moisture | -0.38** | -0.17 | 0.53** | -0.17 | -0.30* | 0.26* | |
litter depth | 0.53** | -0.43** | -0.24* | -0.18 | -0.19 | 0.00 | -0.26* |
Factor analysis loadings for components: (n = 65).
Variable | dim1 | dim2 |
litter cover | 0.075 | -0.905 |
moss cover | -0.816 | 0.380 |
grass cover | 0.719 | 0.566 |
forbs cover | -0.445 | 0.082 |
lichen cover | -0.672 | 0.185 |
soil pH | 0.823 | -0.042 |
soil moisture | 0.517 | 0.564 |
litter depth | 0.161 | -0.780 |
The model-averaged overall estimates of occupancy and detection probability of earthworms were, respectively, 0.83 ± 0.07 and 0.68 ± 0.04. Among the impact levels, the occupancy estimate was highest at road sites (0.90 ± 0.09) and lowest at remote sites (0.73 ± 0.16), although confidence intervals of occupancy at all three human impact levels overlapped considerably (Table
Model-averaged estimates of occupancy (Ψ) and detection probability (p). Uncond. SE: unconditional stand error.
Parameter | Estimate | Uncond. SE | 95% CI |
---|---|---|---|
Ψ | 0.83 | 0.07 | 0.63, 0.93 |
Ψboat launch | 0.84 | 0.08 | 0.59, 0.95 |
Ψroad | 0.90 | 0.09 | 0.42, 0.99 |
Ψremote | 0.73 | 0.16 | 0.38, 0.94 |
Ψconifer | 0.47 | 0.16 | 0.19, 0.76 |
Ψhardwood | 0.91 | 0.06 | 0.72, 0.97 |
p | 0.68 | 0.04 | 0.59, 0.76 |
The top-ranked occupancy model had a weight of 0.31 and included terms for forest type and distance to roads (Table
Top five occupancy models for earthworm occurrence based on the AICc. K: the number of estimated parameters.
Model (occupancy) | Log-likelihood | K | AICc | Δ AICc | Akaike weight |
hardwood + distance | -111.86 | 4 | 232.39 | 0 | 0.31 |
hardwood + boat launch + road | -111.38 | 5 | 233.78 | 1.39 | 0.16 |
hardwood + distance + Dim2 | -111.7 | 5 | 234.42 | 2.03 | 0.11 |
hardwood + distance + Dim1 | -111.71 | 5 | 234.44 | 2.05 | 0.11 |
hardwood + distance + boat launch + road | -111.22 | 6 | 235.88 | 3.5 | 0.05 |
Forest type was included in all highly-ranked models (importance value of 0.98, Table
Model-averaged parameters on logit scale from models of earthworm occurrence. Estimate: model-averaged parameter estimates. SE: Unconditional standard errors.
Parameter | Importance | Estimate | SE | 95% CI |
---|---|---|---|---|
hardwood | 0.98 | 2.45 | 0.86 | 0.76, 4.14 |
distance | 0.69 | -0.30 | 0.16 | -0.62, 0.02 |
boat launch | 0.37 | 1.54 | 1.10 | -0.61, 3.69 |
road | 0.37 | 2.86 | 1.78 | -0.62, 6.34 |
Dim1 | 0.26 | 0.08 | 0.24 | -0.39, 0.55 |
Dim2 | 0.26 | 0.16 | 0.33 | -0.48, 0.8 |
The presence of earthworms did not affect leaf litter depth or other measured soil properties. In a two-way MANOVA of the three factors impact level, forest type, and earthworm occurrence on the three dependent variables of soil pH, soil moisture, and leaf litter depth, the combined dependent variables differed among impact levels (Pilliai’s trace = 0.804, F= 11.7, p < 0.001) and between forest types (Pilliai’s trace = 0.243, F = 5.47, p = 0.002), but the combined variables did not differ between sites where earthworms were present or absent (Pilliai’s trace = 0.037, F = 0.648, p = 0.588). Follow-up univariate ANOVA tests and Bayesian model averaging confirmed that neither soil pH, soil moisture, nor leaf litter depth were affected by the presence or biomass of earthworms.
Exotic earthworms were found to inhabit 90% of road corridors and 85% of boat launch sites, but only 48% of low impact sites. These results suggest that human traffic influences earthworm presence in the KNWR. Similarly,
The road system in the KNWR, while poorly developed compared to conservation units in the contiguous U.S., is fairly extensive compared to other Federal conservation units in Alaska, constituting ~1% of the refuge (100 m buffer either side of all refuge roads gives 6,420 ha). The paved 35 km Sterling Highway and graveled 31 km Skilak Lake Road together bisect the KNWR. These two unpaved roads provide connectivity to many of the 2,900 km of seismic lines (
We found Dendrobaena octaedra to be the most widespread (adults at 70% of study sites) and abundant species (25.6 ± 4.4 ind./m2) of earthworm on KNWR. This species is most likely introduced and spread by vehicles because its small size and epigeic habits (i.e., inhabit near-surface of the leaf litter) likely increase its chances of dispersal by human activities. Dendrobaena octaedra is a prominent invader throughout North America, often both the most widespread and densest exotic earthworm (
Dendrodrilus rubidus was found at only two locations geographically distant from each other, suggesting independent introduction events and perhaps multiple vectors. In an unrelated sampling effort, we have also collected Dd. rubidus in the subalpine zone on the southern portion of the refuge at a site accessible only by floatplane or by foot (http://arctos.database.museum/guid/KNWR:Ento:7100), again suggesting another independent introduction.
Dendrodrilus rubidus
is a fairly common earthworm that appears to be present more in northern hardwood and coniferous forests throughout Alaska and Canada (
In contrast, L. terrestris is an anecic species that lives deep in the soil (
Though we did not detect L. rubellus in our sampling effort, it is presently known to occur on KNWR at only one site, a boat launch on a popular fishing lake. As with L. terrestris, it was most likely introduced by bait abandonment.
In this study, the most important factors determining earthworm occurrence appeared to be forest type (conifer versus hardwoods) followed by distance from roads. We found that earthworms were more likely to be found at sites dominated by deciduous trees and shrubs than at sites dominated by conifers. On KNWR, conifer dominated sites tend to have acidic soils covered by a thick moss carpet, conditions unfavorable to most earthworm species. In contrast, hardwood sites tend to have less acidic soils covered by deciduous leaf litter, providing more ideal conditions for most earthworm species.
Our finding that half (48%) of the low impact sites (> 600 m from any road or facility, and 50 m from any trail or river) contained earthworms was relatively higher than
Besides surface vehicle access, other anthropogenic influences likely contribute to earthworm presence in KNWR, especially in more remote areas. These remote site invasion vectors are not easily identified. Dendrobaena octaedra was found throughout the study area, most likely introduced by road construction, but also possibly by seismic exploration, fire suppression activities, and mechanical tree crushing for moose browse in the northern part of the KNWR during the 1970s. There, extensive seismic lines, mostly in the northern part of the Refuge, have been in place since as early as the 1950s, and many remain visible on the landscape today as animal, hiking, and snowmobile trails, as well as illegal access routes for all-terrain vehicles (
There are also non-anthropogenic vectors that can spread earthworms such as birds (D. Saltmarsh, pers obs.) and streams. In southeast Alaska,
Factors such as soil pH likely also limit earthworm distribution. Most earthworms prefer soil pH of 5–7.4 (
Earthworm densities showed substantial variation with a mean value (28 earthworms/m2) comparable to other studies.
We found no evidence that earthworms were affecting the soil properties pH, soil moisture, and leaf litter depth on KNWR. Likely explanations are the dominance of the epigeic D. octaedra, moderate densities of earthworms, and that these may have been young infestations. In Minnesota,
In temperate studies, earthworm invasions appear to follow a predictable successional sequence, beginning with early invasion by epigeic species, such as D. octaedra, and epi-endogeic species, like L. rubellus. Subsequently endogeic and anecic species like L. terrestris colonize (
Even though differences in cold tolerance do limit the distributions of earthworm species at larger scales (
Regionally in Alaska, the distribution of permafrost and cold winter temperatures, as well as soil moisture and pH, likely limit earthworms’ potential Alaskan distribution. Where earthworms can survive, historic and current human activity and land use practices, and the composition of particular source populations, likely determine earthworm occurrence. The fact that all earthworm records in Alaska up to the present time have been from southern Alaska (see Suppl. material
Both D. octaedra and Dd. rubidus are parthenogenic, frost-hardy species, traits that, combined with their ability to tolerate acidic soils and exploit poor litter quality, contribute to their success in colonizing large areas. As with many “weedy” species, parthenogenesis facilitates rapid reproduction from very low densities, characteristic of rare dispersal events, where a single individual can establish an entire population (
As there are no effective strategies for removing exotic earthworms once they are established, preventing invasion and slowing their dispersal are the only viable ways to reduce their overall impacts. While we found that D. octaedra was widespread, most of the KNWR was free of the more damaging Lumbricus species. Most of the KNWR can be kept free of Lumbricus species for many years due to the extremely slow natural dispersal rate of these worms. Because the main vectors of earthworms on the KNWR appeared to be vehicles and bait abandonment, logical methods for slowing the spread of species already on the landscape and preventing the introduction of additional earthworm species would be to minimize vehicular activity in areas currently devoid of earthworms and to explicitly prohibit the use of earthworms as live bait.
As the first study of earthworm diversity and distribution in the southcentral Alaska region, we established patterns of distribution likely to hold true regionally and we set the stage for considering ways to limit the further introductions of exotic earthworms in Alaska.
We thank Paige Richardson and Jennifer Gregory for their assistance with field and laboratory work, and housing and logistical support from the Kenai National Wildlife Refuge, Jason Geck for GIS assistance, and Carl Tobin for manuscript suggestions. Daniel Fernández Marchán provided helpful comments that improved the manuscript. This research was funded by the NASA Alaska Space Grant Program, the Alaska Pacific University GIS and GPS Grant, and a grant from the USFWS-NWRS Invasive Species Program.