Characterization of biological responses under different environmental conditions: A hierarchical modeling approach

A natural river system is organized as a nested hierarchy of interconnected habitats with specific environmental conditions to which the biological community has adapted. Due to this hierarchical structure, identifying the role of different stressors on the biological community is a formidable task. Efforts trying to link stressors to biological integrity have always been bound to the geographic scale of the selected study area, leading to scale-specific results. In this research, an attempt is made to lift this limitation and develop a hierarchical, scale-sensitive methodology that can identify the significant environmental stressors to the biological community at different scales. Sites with similar background environmental conditions are clustered using self-organizing maps (SOM). This is used to identify stressors which affect the biological community throughout the area of study - called environmental gradients or large-scale stressors. Subsequently, these clusters of similar observations (sampling sites) are progressively sub-divided using environmental variables with a significant but localized effect on the biological community - called small-scale stressors. A parent group of sites is split only when the resulting sub-groups have significantly different biological responses. At the end of this recursive sites decomposition procedure, the original set of observations is organized as a tree of environmentally homogeneous groups of observations characterized by unique biological responses to multiple stressors with different geographic extents. The developed hierarchical analysis methodology has been validated using a large-size dataset of environmental observations from the State of Ohio. Our results show that habitat degradation and increased nutrient loading are the large-scale stressors with a widespread impact in Ohio. Other stressors, such as heavy metals, pH or nitrate concentrations have significant albeit localized effects on biological integrity.     

Modelling landscape effects on density-contact rate relationships of deer in eastern Alberta: Implications for chronic wasting disease

Managing wildlife diseases requires an understanding of disease transmission, which may be strongly affected by host population density and landscape features. Transmission models are typically fit from time-series disease prevalence data and modelled based on how the contact rate among hosts is affected by density, which is often assumed to be a linear (density-dependent transmission) or constant (frequency-dependent transmission) relationship. However, long-term time-series data is unavailable for emerging diseases, and this approach cannot account for independent effects of landscape. We developed a mechanistic model based on ecological data to empirically derive the contact rate-density relationship in white-tailed and mule deer in an enzootic region of chronic wasting disease (CWD) in Alberta, Canada and to determine whether it was affected by landscape. Using data collected from aerial surveys and GPS-telemetry, we developed empirical relationships predicting deer group size, home range size, and habitat selection to iteratively simulate deer distributions across a range of densities and landscapes. We calculated a relative measure of total per-capita contact rate, which is proportional to the number of other deer contacted per individual per unit time, for each distribution as the sum of pairwise contact rates between a target deer and all other individuals. Each pairwise contact rate was estimated from an empirical relationship developed from GPS-telemetry data predicting pairwise contact rates as a function of home range overlap and landscape structure. Total per-capita contact rates increased as a saturating function of density, supporting a transmission model intermediate between density- and frequency-dependent transmission. This pattern resulted from group sizes that reached an asymptote with increasing deer density, although this relationship was mediated by tree and shrub coverage in the landscape, such that in heavily wooded areas, the contact rate saturated at much lower densities. These results suggest that CWD management based on herd reductions, which require a density-dependent contact rate to be effective, may have variable effects on disease across a single management region. The novel mechanistic approach we employed for estimating effects of density and landscape on transmission is a powerful complement to typical data-fitting approaches for modelling disease transmission.     

Lake to land subsidies: experimental addition of aquatic insects increases terrestrial arthropod densities

Aquatic insects are a common and important subsidy to terrestrial systems, yet little is known about how these inputs affect terrestrial food webs, especially around lakes. Myvatn, a lake in northern Iceland, has extraordinary midge (Chironomidae) emergences that result in large inputs of biomass and nutrients to terrestrial arthropod communities. We simulated this lake-to-land resource pulse by collecting midges from Myvatn and spreading their dried carcasses on 1-m2 plots at a nearby site that receives very little midge deposition. We hypothesized a positive bottom-up response of detritivores that would be transmitted to their predators and would persist into the following year. We sampled the arthropod community once per month for two consecutive summers. Midge addition resulted in significantly different arthropod communities and increased densities of some taxa in both years. Detritivores, specifically Diptera larvae, Collembola, and Acari increased in midge-addition plots, and so did some predators and parasitoids. Arthropod densities were still elevated a year after midge addition, and two years of midge addition further increased the density of higher-order consumers (e.g., Coleoptera and Hymenoptera). Midge addition increased arthropod biomass by 68% after one year and 108% after two years. By manipulating the nutrient pulse delivered by midges we were able to elucidate food web consequences of midge deposition and spatial and temporal dynamics that are difficult to determine based on comparative approaches alone. Resources cross ecosystem boundaries and are assimilated over time because of life-history strategies that connect aquatic and terrestrial food webs and these systems cannot be fully understood in isolation from each other.     

Scale-dependent response diversity of seabirds to prey in the North Sea

Functional response diversity is defined as the diversity of responses to environmental change among species that contribute to the same ecosystem function. Because different ecological processes dominate on different spatial and temporal scales, response diversity is likely to be scale dependent. Using three extensive data sets on seabirds, pelagic fish, and zooplankton, we investigate the strength and diversity in the response of seabirds to prey in the North Sea over three scales of ecological organization. Two-stage analyses were used to partition the variance in the abundance of predators and prey among the different scales of investigation: variation from year to year, variation among habitats, and variation on the local patch scale. On the year-to-year scale, we found a strong and synchronous response of seabirds to the abundance of prey, resulting in low response diversity. Conversely, as different seabird species were found in habitats dominated by different prey species, we found a high diversity in the response of seabirds to prey on the habitat scale. Finally, on the local patch scale, seabirds were organized in multispecies patches. These patches were weakly associated with patches of prey, resulting in a weak response strength and a low response diversity. We suggest that ecological similarities among seabird species resulted in low response diversity on the year-to-year scale. On the habitat scale, we suggest that high response diversity was due to interspecific competition and niche segregation among seabird species. On the local patch scale, we suggest that facilitation with respect to the detection and accessibility of prey patches resulted in overlapping distribution of seabirds but weak associations with prey. The observed scale dependencies in response strength and diversity have implications for how the seabird community will respond to different environmental disturbances.     

Probabilistic movement model with emigration simulates movements of deer in Nebraska, 1990–2006

Movements of deer can affect population dynamics, spatial redistribution, and transmission and spread of diseases. Our goal was to model the movement of deer in Nebraska in an attempt to predict the potential for spread of chronic wasting disease (CWD) into eastern Nebraska. We collared and radio-tracked >600 white-tailed deer (Odocoileus virginianus) and mule deer (Odocoileus hemionus) in Nebraska during 1990–2006. We observed large displacements (>10 km) for both species and sexes of deer, including migrations up to 100 km and dispersals up to 50 km. Average distance traveled between successive daily locations was 166 m for male and 173 for female deer in eastern Nebraska, and 427 m for male and 459 for female deer in western Nebraska. Average daily displacement from initial capture point was 10 m for male and 14 m for female deer in eastern Nebraska, and 27 m for male and 28 m for female deer in western Nebraska. We used these data on naturally occurring movements to create and test 6 individual-based models of movement for white-tailed deer and mule deer in Nebraska, including models that incorporated sampling from empirical distributions of movement lengths and turn angles (DIST), correlated random walks (CRW), home point fidelity (FOCUS), shifting home point (SHIFT), probabilistic movement acceptance (MOVE), and probabilistic movement with emigration (MOVEwEMI). We created models in sequence in an attempt to account for the shortcomings of the previous model(s). We used the Kolmogrov–Smirnov goodness-of-fit test to verify improvement of simulated annual displacement distributions to empirical displacement distributions. The best-fit model (D = 0.07 and 0.08 for eastern and western Nebraska, respectively) included a probabilistic movement chance with emigration (MOVEwEMI) and resulted in an optimal daily movement length of 350 m (maximum daily movement length of 2800 m for emigrators) for eastern Nebraska and 370 m (maximum of 2960 m) for western Nebraska. The proportion of deer that moved as emigrators was 0.10 and 0.13 for eastern and western Nebraska, respectively. We propose that the observed spread of CWD may be driven by large movements of a small proportion of deer that help to establish a low prevalence of the disease in areas east of the current endemic area. Our movement models will be used in a larger individual-based simulation of movement, survival, and transmission of CWD to help determine future surveillance and management actions.     

Does the risk of encountering lions influence African herbivore behaviour at waterholes?

A central question in the study of predator–prey relationships is to what extent prey behaviour is determined by avoidance of predators. Here, we test whether the long-term risk of encountering lions and the presence of lions in the vicinity influence the behaviour of large African herbivores at waterholes through avoidance of high-risk areas, increases in group size, changes in temporal niche or changes in the time spent in waterhole areas. In Hwange National Park, Zimbabwe, we monitored waterholes to study the behaviour of nine herbivore species under different risks of encountering lions. We radio-collared 26 lions in the study area which provided the opportunity to monitor whether lions were present during observation sessions and to map longer-term seasonal landscapes of risk of encountering lions. Our results show that the preferred prey species for lions (buffalo, kudu and giraffe) avoided risky waterholes. Group size increased as encounter risk increased for only two species (wildebeest and zebra), but this effect was not strong. Interestingly, buffalo avoided the hours of the day which are dangerous when the long-term and short-term risks of encountering lions were high, and all species showed avoidance of waterhole use at night times when lions were in the vicinity. This illustrates well how prey can make temporal adjustments to avoid dangerous periods coinciding with predator hunting. Additionally, many herbivores spent more time accessing water to drink when the long- and short-term risks of encountering lions were high, and they showed longer potential drinking time when the long-term risk of encountering lions was high, suggesting higher levels of vigilance. This study illustrates the diversity of behavioural adjustments to the risk of encountering a predator and how prey respond differently to temporal variations in this risk.     

Seasonal variation in the diversity and abundance of pelagic larvae of Antarctic marine invertebrates

Most marine benthic macroinvertebrate species reproduce via a larval phase but attempts to explain the occurrence of different larval strategies (feeding or non-feeding, pelagic or benthic) in different habitats have been largely inconclusive. There have been very few year-round surveys of meroplankton at any latitude and in consequence fundamental data on the diversity, abundance, and timings of larval life history phases are lacking. There has been considerable debate regarding the viability of pelagic larvae in cold waters with highly seasonal primary production but there has been only one year-round study of meroplankton in the Southern Ocean, and that was outside of the Antarctic Circle. We present data from the first year-round survey of meroplankton assemblages at a location within the Antarctic Circle. We surveyed abundances of meroplanktonic larvae over 1.5 year at Rothera Point, West Antarctic Peninsula (67°34′S, 68°07′W). Larvae were collected in monthly diver-towed net samples close to the seabed at 20 and 6 m total water depths at each of three locations and were identified and counted live immediately after sampling. A total of 99 operationally defined taxonomic types representing 11 phyla were recorded but this is likely to be an underestimate of true diversity because of inherent difficulties of identification. Larvae were present in all months of the year and although planktotrophic larvae were more abundant in summer, both feeding and non-feeding types were present in all months. Comparisons of seasonal larval abundances with data from a settlement study at the same sites and from the literature show that larvae of mobile adults settle in summer regardless of developmental type, whereas sessile taxa settle in all seasons. We suggest that this is a consequence of differences in the food requirements of mobile and sessile fauna and that the availability of food for post-larval juveniles is more critical for survival than factors affecting the larval stage itself.