Extinctions of herbivorous mammals in the late Pleistocene of Australia in relation to their feeding ecology: No evidence for environmental change as cause of extinction

Abstract There has been debate over the cause of the extinction of ‘megafauna’ species during the late Pleistocene of Australia. One view is that environmental change, either natural or human‐induced, was the main factor in the extinctions. Some support for this comes from the observation that, among herbivores, most of the species that went extinct were apparently browsers rather than grazers. Browsers would presumably have been more dependent on shrubland and woodland habitats than grazers, and it has been argued that such habitats might have contracted in response to aridity or changed fire regimes in the late Pleistocene. Here, we test this idea by comparing extinction rates of browsers and grazers in the late Pleistocene, controlling for body mass in both groups. We show that in both browsers and grazers the probability of extinction was very strongly related to body mass, and the body mass at which extinction became likely was similar in the two groups. It is true that more browsers than grazers went extinct, but this is largely because most very large herbivores in the late Pleistocene were browsers, not because large browsers were more likely to go extinct than similarly sized grazers. This result provides evidence against some forms of environmental change as a cause of the extinctions.     

Forecasting species range dynamics with process‐explicit models: matching methods to applications

Knowing where species occur is fundamental to many ecological and environmental applications. Species distribution models (SDMs) are typically based on correlations between species occurrence data and environmental predictors, with ecological processes captured only implicitly. However, there is a growing interest in approaches that explicitly model processes such as physiology, dispersal, demography and biotic interactions. These models are believed to offer more robust predictions, particularly when extrapolating to novel conditions. Many process–explicit approaches are now available, but it is not clear how we can best draw on this expanded modelling toolbox to address ecological problems and inform management decisions. Here, we review a range of process–explicit models to determine their strengths and limitations, as well as their current use. Focusing on four common applications of SDMs – regulatory planning, extinction risk, climate refugia and invasive species – we then explore which models best meet management needs. We identify barriers to more widespread and effective use of process‐explicit models and outline how these might be overcome. As well as technical and data challenges, there is a pressing need for more thorough evaluation of model predictions to guide investment in method development and ensure the promise of these new approaches is fully realised.     

Population dynamics can be more important than physiological limits for determining range shifts under climate change

Evidence is accumulating that species' responses to climate changes are best predicted by modelling the interaction of physiological limits, biotic processes and the effects of dispersal‐limitation. Using commercially harvested blacklip (Haliotis rubra) and greenlip abalone (Haliotis laevigata) as case studies, we determine the relative importance of accounting for interactions among physiology, metapopulation dynamics and exploitation in predictions of range (geographical occupancy) and abundance (spatially explicit density) under various climate change scenarios. Traditional correlative ecological niche models (ENM) predict that climate change will benefit the commercial exploitation of abalone by promoting increased abundances without any reduction in range size. However, models that account simultaneously for demographic processes and physiological responses to climate‐related factors result in future (and present) estimates of area of occupancy (AOO) and abundance that differ from those generated by ENMs alone. Range expansion and population growth are unlikely for blacklip abalone because of important interactions between climate‐dependent mortality and metapopulation processes; in contrast, greenlip abalone should increase in abundance despite a contraction in AOO. The strongly non‐linear relationship between abalone population size and AOO has important ramifications for the use of ENM predictions that rely on metrics describing change in habitat area as proxies for extinction risk. These results show that predicting species' responses to climate change often require physiological information to understand climatic range determinants, and a metapopulation model that can make full use of this data to more realistically account for processes such as local extirpation, demographic rescue, source‐sink dynamics and dispersal‐limitation.     

How complex should models be? Comparing correlative and mechanistic range dynamics models

Criticism has been levelled at climate‐change‐induced forecasts of species range shifts that do not account explicitly for complex population dynamics. The relative importance of such dynamics under climate change is, however, undetermined because direct tests comparing the performance of demographic models vs. simpler ecological niche models are still lacking owing to difficulties in evaluating forecasts using real‐world data. We provide the first comparison of the skill of coupled ecological‐niche‐population models and ecological niche models in predicting documented shifts in the ranges of 20 British breeding bird species across a 40‐year period. Forecasts from models calibrated with data centred on 1970 were evaluated using data centred on 2010. We found that more complex coupled ecological‐niche‐population models (that account for dispersal and metapopulation dynamics) tend to have higher predictive accuracy in forecasting species range shifts than structurally simpler models that only account for variation in climate. However, these better forecasts are achieved only if ecological responses to climate change are simulated without static snapshots of historic land use, taken at a single point in time. In contrast, including both static land use and dynamic climate variables in simpler ecological niche models improve forecasts of observed range shifts. Despite being less skilful at predicting range changes at the grid‐cell level, ecological niche models do as well, or better, than more complex models at predicting the magnitude of relative change in range size. Therefore, ecological niche models can provide a reasonable first approximation of the magnitude of species' potential range shifts, especially when more detailed data are lacking on dispersal dynamics, demographic processes underpinning population performance, and change in land cover.     

How to understand species’ niches and range dynamics: a demographic research agenda for biogeography

Range dynamics causes mismatches between a species’ geographical distribution and the set of suitable environments in which population growth is positive (the Hutchinsonian niche). This is because source–sink population dynamics cause species to occupy unsuitable environments, and because environmental change creates non‐equilibrium situations in which species may be absent from suitable environments (due to migration limitation) or present in unsuitable environments that were previously suitable (due to time‐delayed extinction). Because correlative species distribution models do not account for these processes, they are likely to produce biased niche estimates and biased forecasts of future range dynamics. Recently developed dynamic range models (DRMs) overcome this problem: they statistically estimate both range dynamics and the underlying environmental response of demographic rates from species distribution data. This process‐based statistical approach qualitatively advances biogeographical analyses. Yet, the application of DRMs to a broad range of species and study systems requires substantial research efforts in statistical modelling, empirical data collection and ecological theory. Here we review current and potential contributions of these fields to a demographic understanding of niches and range dynamics. Our review serves to formulate a demographic research agenda that entails: (1) advances in incorporating process‐based models of demographic responses and range dynamics into a statistical framework, (2) systematic collection of data on temporal changes in distribution and abundance and on the response of demographic rates to environmental variation, and (3) improved theoretical understanding of the scaling of demographic rates and the dynamics of spatially coupled populations. This demographic research agenda is challenging but necessary for improved comprehension and quantification of niches and range dynamics. It also forms the basis for understanding how niches and range dynamics are shaped by evolutionary dynamics and biotic interactions. Ultimately, the demographic research agenda should lead to deeper integration of biogeography with empirical and theoretical ecology.     

Colonization and persistence ability explain the extent to which plant species fill their potential range

Aim How species traits and environmental conditions affect biogeographical dynamics is poorly understood. Here we test whether estimates of a species’ evolutionary age, colonization and persistence ability can explain its current ‘range filling’ (the ratio between realized and potential range size). Location Fynbos biome (Cape Floristic Region, South Africa). Methods For 37 species of woody plants (Proteaceae), we estimate range filling using atlas data and distribution models, evolutionary age using molecular phylogenies, and persistence ability using estimates of individual longevity (which determines the probability of extinction of local populations). Colonization ability is estimated from validated process‐based seed dispersal models, the arrangement of potential habitat, and data on local abundance. To relate interspecific variation in range filling to evolutionary age, colonization and persistence ability, we use two complementary model types: phenomenological linear models and the process‐based metapopulation model of Levins. Results Linear model analyses show that range filling increases with a species’ colonization and persistence ability but is not affected by species age. Moreover, colonization ability is a better predictor of range filling than its component variables (local abundance and dispersal ability). The phylogenetically independent interaction between colonization and persistence ability is significant (P < 0.05) for 97% of 180 alternative phylogenies. While the selected linear model explains 42% of the variance in arcsine transformed range filling, the Levins model performs more poorly. It overestimates range filling for realistic parameter values and produces unrealistic parameter estimates when fitted statistically. Main conclusions Colonization and local extinction seem to shape Proteaceae range dynamics on ecological rather than macroevolutionary time‐scales. Our results suggest that the positive abundance–range size relationship in this group is due primarily to the effect of abundance on colonization. In summary, this study contributes to a process‐based understanding of range dynamics and highlights the importance of colonization for the future survival of Fynbos Proteaceae.     

Estimating indices of range shifts in birds using dynamic models when detection is imperfect

There is intense interest in basic and applied ecology about the effect of global change on current and future species distributions. Projections based on widely used static modeling methods implicitly assume that species are in equilibrium with the environment and that detection during surveys is perfect. We used multiseason correlated detection occupancy models, which avoid these assumptions, to relate climate data to distributional shifts of Louisiana Waterthrush in the North American Breeding Bird Survey (BBS) data. We summarized these shifts with indices of range size and position and compared them to the same indices obtained using more basic modeling approaches. Detection rates during point counts in BBS surveys were low, and models that ignored imperfect detection severely underestimated the proportion of area occupied and slightly overestimated mean latitude. Static models indicated Louisiana Waterthrush distribution was most closely associated with moderate temperatures, while dynamic occupancy models indicated that initial occupancy was associated with diurnal temperature ranges and colonization of sites was associated with moderate precipitation. Overall, the proportion of area occupied and mean latitude changed little during the 1997–2013 study period. Near‐term forecasts of species distribution generated by dynamic models were more similar to subsequently observed distributions than forecasts from static models. Occupancy models incorporating a finite mixture model on detection – a new extension to correlated detection occupancy models – were better supported and may reduce bias associated with detection heterogeneity. We argue that replacing phenomenological static models with more mechanistic dynamic models can improve projections of future species distributions. In turn, better projections can improve biodiversity forecasts, management decisions, and understanding of global change biology.     

Conservation management and sustainable harvest quotas are sensitive to choice of climate modelling approach for two marine gastropods

Aim

To establish the robustness of two alternative methods for predicting the future ranges and abundances for two wild‐harvested abalone species (Haliotis rubra Donovan 1808 and H. laevigata Leach 1814): single atmosphere–ocean general circulation model (GCM) or ensemble‐averaged GCM forecasts.

Location

South Australia.

Methods

We assessed the ability of 20 GCMs to simulate observed seasonal sea surface temperature (SST) between 1980–1999, globally, and regionally for the Indian and Pacific Oceans south of the Equator. We used model rankings to characterize a set of representative climate futures, using three different‐sized GCM ensembles and two individual GCMs (the Parallel Climate Model and the Community Climate System Model, version 3.0). Ecological niche models were then coupled to physiological information to compare forecast changes in area of occupancy, population size and harvest area based on forecasts using the various GCM selection methods, as well as different greenhouse gas emission scenarios and climate sensitivities.

Results

We show that: (1) the skill with which climate models reproduce recent SST records varies considerably amongst GCMs, with multimodel ensemble averages showing closer agreement to observations than single models; (2) choice of GCM, and the decision on whether or not to use ensemble‐averaged climate forecasts, can strongly influence spatiotemporal predictions of range, abundance and fishing potential; and (3) comparable hindcasting skill does not necessarily guarantee that GCM forecasts and ecological and evolutionary responses to these forecast changes, will be similar amongst closely ranked models.

Conclusion

By averaging across an ensemble of seven highly ranked skilful GCMs, inherent uncertainties stemming from GCM differences are incorporated into forecasts of change in species range, abundance and sustainable fishing area. Our results highlight the need to make informed and explicit decisions on GCM choice, model sensitivity and emission scenarios when exploring conservation options for marine species and the sustainability of future harvests using ecological niche models.

     

Informing management decisions for ecological networks,using dynamic models calibrated to noisy time‐series data

Well‐intentioned environmental management can backfire, causing unforeseen damage. To avoid this, managers and ecologists seek accurate predictions of the ecosystem‐wide impacts of interventions, given small and imprecise datasets, which is an incredibly difficult task. We generated and analysed thousands of ecosystem population time series to investigate whether fitted models can aid decision‐makers to select interventions. Using these time‐series data (sparse and noisy datasets drawn from deterministic Lotka‐Volterra systems with two to nine species, of known network structure), dynamic model forecasts of whether a species’ future population will be positively or negatively affected by rapid eradication of another species were correct > 70% of the time. Although 70% correct classifications is only slightly better than an uninformative prediction (50%), this classification accuracy can be feasibly improved by increasing monitoring accuracy and frequency. Our findings suggest that models may not need to produce well‐constrained predictions before they can inform decisions that improve environmental outcomes.     

Responses to atmospheric CO2 concentrations in crop simulation models: a review of current simple and semicomplex representations and options for model development

Elevated atmospheric CO₂ concentrations ([CO₂]) cause direct changes in crop physiological processes (e.g. photosynthesis and stomatal conductance). To represent these CO₂ responses, commonly used crop simulation models have been amended, using simple and semicomplex representations of the processes involved. Yet, there is no standard approach to and often poor documentation of these developments. This study used a bottom‐up approach (starting with the APSIM framework as case study) to evaluate modelled responses in a consortium of commonly used crop models and illuminate whether variation in responses reflects true uncertainty in our understanding compared to arbitrary choices of model developers. Diversity in simulated CO₂ responses and limited validation were common among models, both within the APSIM framework and more generally. Whereas production responses show some consistency up to moderately high [CO₂] (around 700 ppm), transpiration and stomatal responses vary more widely in nature and magnitude (e.g. a decrease in stomatal conductance varying between 35% and 90% among models was found for [CO₂] doubling to 700 ppm). Most notably, nitrogen responses were found to be included in few crop models despite being commonly observed and critical for the simulation of photosynthetic acclimation, crop nutritional quality and carbon allocation. We suggest harmonization and consideration of more mechanistic concepts in particular subroutines, for example, for the simulation of N dynamics, as a way to improve our predictive understanding of CO₂ responses and capture secondary processes. Intercomparison studies could assist in this aim, provided that they go beyond simple output comparison and explicitly identify the representations and assumptions that are causal for intermodel differences. Additionally, validation and proper documentation of the representation of CO₂ responses within models should be prioritized.