- On the setting of environmental noise and the performance of population dynamical models (2010)
- Background: Environmental noise is ubiquitous in population growth processes, with a well acknowledged potential to affect populations regardless of their sizes. It therefore deserves consideration in population dynamics modelling. The usual approach to incorporating noise into population dynamical models is to make some model parameter(s) (typically the growth rate, the carrying capacity, or both) stochastic and responsive to environment fluctuations. It is however still unclear whether including noise in one or/and another parameter makes a difference to the model performance. Here we investigated this issue with a focus on model fit and predictive accuracy. To do this, we developed three population dynamical models of the Ricker type with the noise included in the growth rate (Model 1), in the carrying capacity (Model 2), and in both (Model 3). We generated several population time series under each model, and used a Bayesian approach to fit the three models to the simulated data. We then compared the model performances in fitting to the data and in forecasting future observations. Results: When the mean intrinsic growth rate, r, in the data was low, the three models had roughly comparable performances, irrespective of the true model and the level of noise. As r increased, Models 1 performed best on data generated from it, and Model 3 tended to perform best on data generated from either Models 2 or Model 3. Model 2 was uniformly outcompeted by the other two models, regardless of the true model and the level of noise. The correlation between the deviance information criterion (DIC) and the mean square error (MSE) used respectively as measure of fit and predictive accuracy was broadly positive. Conclusion: Our results suggested that the way environmental noise is incorporated into a population dynamical model may profoundly affect its performance. Overall, we found that including noise in one or/and another parameter does not matter as long as the mean intrinsic growth rate, r, is low. As r increased, however, the three models performed differently. Models 1 and 3 broadly outperformed Model 2, the first having the advantage of being simple and more computationally tractable. A comforting result emerging from our analysis is the broad positive correlation between MSEs and DICs, suggesting that the latter may also be informative about the predictive performance of a model.
- Hierarchical modelling of temperature and habitat size effects on population dynamics of North Atlantic cod (2010)
- Understanding how temperature affects cod (Gadus morhua) ecology is important for forecasting how populations will develop as climate changes in future. The effects of spawning-season temperature and habitat size on cod recruitment dynamics have been investigated across the North Atlantic. Ricker and Beverton and Holt stock–recruitment (SR) models were extended by applying hierarchical methods, mixed-effects models, and Bayesian inference to incorporate the influence of these ecosystem factors on model parameters representing cod maximum reproductive rate and carrying capacity. We identified the pattern of temperature effects on cod productivity at the species level and estimated SR model parameters with increased precision. Temperature impacts vary geographically, being positive in areas where temperatures are <5°C, and negative for higher temperatures. Using the relationship derived, it is possible to predict expected changes in population-specific reproductive rates and carrying capacities resulting from temperature increases. Further, carrying capacity covaries with available habitat size, explaining at least half its variability across stocks. These patterns improve our understanding of environmental impacts on key population parameters, which is required for an ecosystem approach to cod management, particularly under ocean-warming scenarios. Key words: carrying capacity , cod , hierarchical models , North Atlantic , temperature , uncertainty
- Host genotype shapes the foliar fungal microbiome of balsam poplar (Populus balsamifera) (2013)
- Foliar fungal communities of plants are diverse and ubiquitous. In grasses endophytes may increase host fitness; in trees, their ecological roles are poorly understood. We investigated whether the genotype of the host tree influences community structure of foliar fungi. We sampled leaves from genotyped balsam poplars from across the species' range, and applied 454 amplicon sequencing to characterize foliar fungal communities. At the time of the sampling the poplars had been growing in a common garden for two years. We found diverse fungal communities associated with the poplar leaves. Linear discriminant analysis and generalized linear models showed that host genotypes had a structuring effect on the composition of foliar fungal communities. The observed patterns may be explained by a filtering mechanism which allows the trees to selectively recruit fungal strains from the environment. Alternatively, host genotype-specific fungal communities may be present in the tree systemically, and persist in the host even after two clonal reproductions. Both scenarios are consistent with host tree adaptation to specific foliar fungal communities and suggest that there is a functional basis for the strong biotic interaction.
- Variable strength of forest stand attributes and weather conditions on the questing activity of Ixodes ricinus ticks over years in managed forests (2013)
- Given the ever-increasing human impact through land use and climate change on the environment, we crucially need to achieve a better understanding of those factors that influence the questing activity of ixodid ticks, a major disease-transmitting vector in temperate forests. We investigated variation in the relative questing nymph densities of Ixodes ricinus in differently managed forest types for three years (2008–2010) in SW Germany by drag sampling. We used a hierarchical Bayesian modeling approach to examine the relative effects of habitat and weather and to consider possible nested structures of habitat and climate forces. The questing activity of nymphs was considerably larger in young forest successional stages of thicket compared with pole wood and timber stages. Questing nymph density increased markedly with milder winter temperatures. Generally, the relative strength of the various environmental forces on questing nymph density differed across years. In particular, winter temperature had a negative effect on tick activity across sites in 2008 in contrast to the overall effect of temperature across years. Our results suggest that forest management practices have important impacts on questing nymph density. Variable weather conditions, however, might override the effects of forest management practices on the fluctuations and dynamics of tick populations and activity over years, in particular, the preceding winter temperatures. Therefore, robust predictions and the detection of possible interactions and nested structures of habitat and climate forces can only be quantified through the collection of long-term data. Such data are particularly important with regard to future scenarios of forest management and climate warming.
- Dealing with varying detection probability, unequal sample sizes and clumped distributions in count data (2012)
- Temporal variation in the detectability of a species can bias estimates of relative abundance if not handled correctly. For example, when effort varies in space and/or time it becomes necessary to take variation in detectability into account when data are analyzed. We demonstrate the importance of incorporating seasonality into the analysis of data with unequal sample sizes due to lost traps at a particular density of a species. A case study of count data was simulated using a spring-active carabid beetle. Traps were ‘lost’ randomly during high beetle activity in high abundance sites and during low beetle activity in low abundance sites. Five different models were fitted to datasets with different levels of loss. If sample sizes were unequal and a seasonality variable was not included in models that assumed the number of individuals was log-normally distributed, the models severely under- or overestimated the true effect size. Results did not improve when seasonality and number of trapping days were included in these models as offset terms, but only performed well when the response variable was specified as following a negative binomial distribution. Finally, if seasonal variation of a species is unknown, which is often the case, seasonality can be added as a free factor, resulting in well-performing negative binomial models. Based on these results we recommend (a) add sampling effort (number of trapping days in our example) to the models as an offset term, (b) if precise information is available on seasonal variation in detectability of a study object, add seasonality to the models as an offset term; (c) if information on seasonal variation in detectability is inadequate, add seasonality as a free factor; and (d) specify the response variable of count data as following a negative binomial or over-dispersed Poisson distribution.
- Southern high-latitude terrestrial climate change during the Paleocene–Eocene derived from a marine pollen record (ODP Site 1172, East Tasman Plateau) (2014)
- Reconstructing the early Paleogene climate dynamics of terrestrial settings in the high southern latitudes is important to assess the role of high-latitude physical and biogeochemical processes in the global climate system. However, whereas a number of high-quality Paleogene climate records has become available for the marine realm of the high southern latitudes over the recent past, the long-term evolution of coeval terrestrial climates and ecosystems is yet poorly known. We here explore the climate and vegetation dynamics on Tasmania from the middle Paleocene to the early Eocene (60.7–54.2 Ma) based on a sporomorph record from Ocean Drilling Program (ODP) Site 1172 on the East Tasman Plateau. Our results show that three distinctly different vegetation types thrived on Tasmania under a high-precipitation regime during the middle Paleocene to early Eocene, with each type representing different temperature conditions: (i) warm-temperate forests dominated by gymnosperms that were dominant during the middle and late Paleocene; (ii) cool-temperate forests dominated by southern beech (Nothofagus) and araucarians across the middle/late Paleocene transition interval (~59.5 to ~59.0 Ma); and (iii) paratropical forests rich in ferns that were established during and in the wake of the Paleocene–Eocene Thermal Maximum (PETM). The transient establishment of cool-temperate forests lacking any frost-sensitive elements (i.e., palms and cycads) across the middle/late Paleocene transition interval indicates markedly cooler conditions, with the occurrence of frosts in winter, on Tasmania during that time. The integration of our sporomorph data with previously published TEX86-based sea-surface temperatures from ODP Site 1172 documents that the vegetation dynamics on Tasmania were closely linked with the temperature evolution in the Tasman sector of the Southwest Pacific region. Moreover, the comparison of our season-specific climate estimates for the sporomorph assemblages from ODP Site 1172 with the TEX86L- and TEX86H-based temperature data suggests a warm-season bias of both calibrations for the early Paleogene of the high southern latitudes.
- Southern high-latitude terrestrial climate change during the Palaeocene–Eocene derived from a marine pollen record (ODP Site 1172, East Tasman Plateau) (2014)
- Global warming, changes in the hydrological cycle and enhanced marine primary productivity all have been invoked to have contributed to the occurrence of widespread ocean anoxia during the Cenomanian-Turonian Oceanic Anoxic Event (OAE2; ~ 94 Ma), but disentangling these factors on a regional scale has remained problematic. We generated palynological and organic geochemical records that allow the separation of these forcing factors in a core spanning the OAE2 from Wunstorf, Lower Saxony Basin (LSB; North Gemany), which exhibits cyclic black shale–marl alternations related to the orbital precession cycle. Despite the widely varying depositional conditions complicating the interpretation of the obtained records, TEX86H indicates that sea-surface temperature (SST) evolution in the LSB during OAE2 resembles that of previously studied sites throughout the proto-North Atlantic. Cooling during the so-called Plenus Cold Event interrupted black shale deposition during the early stages of OAE2. However, TEX86 does not vary significantly across marl–black shale alternations, suggesting that temperature variations did not force the formation of the cyclic black shale horizons. Relative (i.e., with respect to marine palynomorphs) and absolute abundances of pollen and spores are elevated during phases of black shale deposition, indicative of enhanced precipitation and run-off. High abundances of cysts from inferred heterotrophic and euryhaline dinoflagellates supports high run-off, which likely introduced additional nutrients to the epicontinental shelf resulting in elevated marine primary productivity. We conclude that orbitally-forced enhanced precipitation and run-off, in tandem with elevated marine primary productivity, were critical in cyclic black shale formation on the northwest European epicontinental shelf and potentially for other OAE2 sections in the proto-Atlantic and Western Interior Seaway at similar latitudes as well.
- Risk Factors for the Presence of Chikungunya and Dengue Vectors (Aedes aegypti and Aedes albopictus), Their Altitudinal Distribution and Climatic Determinants of Their Abundance in Central Nepal (2015)
- Background:The presence of the recently introduced primary dengue virus vector mosquito Aedes aegypti in Nepal, in association with the likely indigenous secondary vector Aedes albopictus, raises public health concerns. Chikungunya fever cases have also been reported in Nepal, and the virus causing this disease is also transmitted by these mosquito species. Here we report the results of a study on the risk factors for the presence of chikungunya and dengue virus vectors, their elevational ceiling of distribution, and climatic determinants of their abundance in central Nepal. Methodology/Principal Findings: We collected immature stages of mosquitoes during six monthly cross-sectional surveys covering six administrative districts along an altitudinal transect in central Nepal that extended from Birgunj (80 m above sea level [asl]) to Dhunche (highest altitude sampled: 2,100 m asl). The dengue vectors Ae. aegypti and Ae. albopictus were commonly found up to 1,350 m asl in Kathmandu valley and were present but rarely found from 1,750 to 2,100 m asl in Dhunche. The lymphatic filariasis vector Culex quinquefasciatus was commonly found throughout the study transect. Physiographic region, month of collection, collection station and container type were significant predictors of the occurrence and co-occurrence of Ae. aegypti and Ae. albopictus. The climatic variables rainfall, temperature, and relative humidity were significant predictors of chikungunya and dengue virus vectors abundance. Conclusions/Significance: We conclude that chikungunya and dengue virus vectors have already established their populations up to the High Mountain region of Nepal and that this may be attributed to the environmental and climate change that has been observed over the decades in Nepal. The rapid expansion of the distribution of these important disease vectors in the High Mountain region, previously considered to be non-endemic for dengue and chikungunya fever, calls for urgent actions to protect the health of local people and tourists travelling in the central Himalayas.