Climatically driven changes in primary production propagate through trophic levels

Relationships of primary productivity with anuran abundance, richness, and community composition in tropical streams

The relationship between primary productivity and diversity has been demonstrated across taxa and spatial scales, but for organisms with biphasic life cycles, little research has examined whether productivity of larval and adult environments influence each life stage independently, or whether productivity of one life stage's environment outweighs the influence of the other. Experimental work demonstrates that tadpoles of stream-breeding anurans can exhibit a top-down influence on aquatic primary productivity (APP), but few studies have sought evidence of a bottom-up influence of primary productivity on anuran abundance, species richness and community composition, as seen in other organisms. We examined aquatic and terrestrial primary productivity in two forest types in Borneo, along with amphibian abundance, species richness, and community composition at larval and adult stages, to determine whether there is evidence for a bottom-up influence of APP on tadpole abundance and species richness across streams, and the relative importance of aquatic and terrestrial primary productivity on larval and adult phases of anurans. We predicted that adult richness, abundance, and community composition would be influenced by terrestrial primary productivity, but that tadpole richness, abundance, and community composition would be influenced by APP. Contrary to expectations, we did not find evidence that primary productivity, or variation thereof, predicts anuran richness at larval or adult stages. Further, no measure of primary productivity or its variation was a significant predictor of adult abundance, or of adult or tadpole community composition. For tadpoles, we found that in areas with low terrestrial primary productivity, abundance was positively related to APP, but in areas with high terrestrial primary productivity, abundance was negatively related to APP, suggesting a bottom-up influence of primary productivity on abundance in secondary forest, and a top-down influence of tadpoles on primary productivity in primary forest. Additional data are needed to better understand the ecological interactions between terrestrial primary productivity, aquatic primary productivity, and tadpole abundance.

Early post-settlement events, rather than settlement, drive recruitment and coral recovery at Moorea, French Polynesia

Understanding population dynamics is a long-standing objective of ecology, but the need for progress in this area has become urgent. For coral reefs, achieving this objective is impeded by a lack of information on settlement versus post-settlement events in determining recruitment and population size. Declines in coral abundance are often inferred to be associated with reduced densities of recruits, which could arise from mechanisms occurring at larval settlement, or throughout post-settlement stages. This study uses annual measurements from 2008 to 2021 of coral cover, the density of coral settlers (S), the density of small corals (SC), and environmental conditions, to evaluate the roles of settlement versus post-settlement events in determining rates of coral recruitment and changes in coral cover at Moorea, French Polynesia. Coral cover, S, SC, and the SC:S ratio (a proxy for post-settlement success), and environmental conditions, were used in generalized additive models (GAMs) to show that: (a) coral cover was more strongly related to SC and SC:S than S, and (b) SC:S was highest when preceded by cool seawater, low concentrations of Chlorophyll a, and low flow speeds, and S showed evidence of declining with elevated temperature. Together, these results suggest that changes in coral cover in Moorea are more strongly influenced by post-settlement events than settlement. The key to understanding coral community resilience may lie in elucidating the factors attenuating the bottleneck between settlers and small corals.

Stepping stones to extirpation: Puma patch occupancy thresholds in an urban-wildland matrix

Habitat loss and fragmentation are the leading causes of species range contraction and extirpation, worldwide. Factors that predict sensitivity to fragmentation include high trophic level, large body size, and extensive spatial requirements. Pumas (Puma concolor) exemplify these qualities, making them particularly susceptible to fragmentation and subsequent reductions in demographic connectivity. The chaparral-dominated ecosystems surrounding the greater San Francisco Bay Area encompass over 10,000 km2 of suitable puma habitat, but inland waterways, croplands, urban land uses, and extensive transportation infrastructure have resulted in widespread habitat fragmentation. Pumas in this region now exist as a metapopulation marked by loss of genetic diversity, collisions with vehicles, and extensive human-puma conflict. Given these trends, we conducted a photo survey from 2017 to 2021 across 19 patches of predicted habitat and compiled a dataset of >6584 puma images. We used a logistic regression analytical framework to evaluate the hypothesis that puma patch occupancy would exhibit a threshold response explained by patch size, isolation, and habitat quality. Contrary to predictions, only variables related to patch size demonstrated any power to explain occupancy. On average, occupied patches were 18× larger than those where they were not detected (825 ± 1238 vs. 46 ± 101 km2). Although we observed pumas in patches as small as 1 km2, logistic regression models indicated a threshold occupancy probability between 300 and 400 km2, which is remarkably close to the mean male puma home range size in coastal California (~381 km2). Puma populations dependent on habitats below this value may be susceptible to inbreeding depression and human-wildlife conflict, and therefore vulnerable to extirpation. For species conservation, we suggest conflicts might be ameliorated by identifying the largest, isolated patches for public education campaigns with respect to management of domestic animals, and remaining connective parcels be identified, mapped, and prioritized for targeted mitigation.

Ecosystem productivity affected the spatiotemporal disappearance of Neanderthals in Iberia

What role did fluctuations play in biomass availability for secondary consumers in the disappearance of Neanderthals and the survival of modern humans? To answer this, we quantify the effects of stadial and interstadial conditions on ecosystem productivity and human spatiotemporal distribution patterns during the Middle to Upper Palaeolithic transition (50,000-30,000 calibrated years before the present) in Iberia. First, we used summed probability distribution, optimal linear estimation and Bayesian age modelling to reconstruct an updated timescale for the transition. Next, we executed a generalized dynamic vegetation model to estimate the net primary productivity. Finally, we developed a macroecological model validated with present-day observations to calculate herbivore abundance. The results indicate that, in the Eurosiberian region, the disappearance of Neanderthal groups was contemporaneous with a significant decrease in the available biomass for secondary consumers, and the arrival of the first Homo sapiens populations coincided with an increase in herbivore carrying capacity. During stadials, the Mediterranean region had the most stable conditions and the highest biomass of medium and medium-large herbivores. These outcomes support an ecological cause for the hiatus between the Mousterian and Aurignacian technocomplexes in Northern Iberia and the longer persistence of Neanderthals in southern latitudes.

Net effect of environmental fluctuations in multiple global-change drivers across the tree of life

Jensen's inequality predicts that the response of any given system to average constant conditions is different from its average response to varying ones. Environmental fluctuations in abiotic conditions are pervasive on Earth; yet until recently, most ecological research has addressed the effects of multiple environmental drivers by assuming constant conditions. One could thus expect to find significant deviations in the magnitude of their effects on ecosystems when environmental fluctuations are considered. Drawing on experimental studies published during the last 30 years reporting more than 950 response ratios (n = 5,700), we present a comprehensive analysis of the role that environmental fluctuations play across the tree of life. In contrast to the predominance of interactive effects of global-change drivers reported in the literature, our results show that their cumulative effects were additive (58%), synergistic (26%), and antagonistic (16%) when environmental fluctuations were present. However, the dominant type of interaction varied by trophic level (autotrophs: interactive; heterotrophs: additive) and phylogenetic group (additive in Animalia; additive and positive antagonism in Chromista; negative antagonism and synergism in Plantae). In addition, we identify the need to tackle how complex communities respond to fluctuating environments, widening the phylogenetic and biogeographic ranges considered, and to consider other drivers beyond warming and acidification as well as longer timescales. Environmental fluctuations must be taken into account in experimental and modeling studies as well as conservation plans to better predict the nature, magnitude, and direction of the impacts of global change on organisms and ecosystems.

Productivity, niche availability, species richness, and extinction risk: Untangling relationships using individual-based simulations

It has often been suggested that the productivity of an ecosystem affects the number of species that it can support. Despite decades of study, the nature, extent, and underlying mechanisms of this relationship are unclear. One suggested mechanism is the "more individuals" hypothesis (MIH). This proposes that productivity controls the number of individuals in the ecosystem, and that more individuals can be divided into a greater number of species before their population size is sufficiently small for each to be at substantial risk of extinction. Here, we test this hypothesis using REvoSim: an individual-based eco-evolutionary system that simulates the evolution and speciation of populations over geological time, allowing phenomena occurring over timescales that cannot be easily observed in the real world to be evaluated. The individual-based nature of this system allows us to remove assumptions about the nature of speciation and extinction that previous models have had to make. Many of the predictions of the MIH are supported in our simulations: Rare species are more likely to undergo extinction than common species, and species richness scales with productivity. However, we also find support for relationships that contradict the predictions of the strict MIH: species population size scales with productivity, and species extinction risk is better predicted by relative than absolute species population size, apparently due to increased competition when total community abundance is higher. Furthermore, we show that the scaling of species richness with productivity depends upon the ability of species to partition niche space. Consequently, we suggest that the MIH is applicable only to ecosystems in which niche partitioning has not been halted by species saturation. Some hypotheses regarding patterns of biodiversity implicitly or explicitly overlook niche theory in favor of neutral explanations, as has historically been the case with the MIH. Our simulations demonstrate that niche theory exerts a control on the applicability of the MIH and thus needs to be accounted for in macroecology.

Ecosystem response to earlier ice break-up date: Climate-driven changes to water temperature, lake-habitat-specific production, and trout habitat and resource use

Climate warming has yielded earlier ice break-up dates in recent decades for lakes leading to water temperature increases, altered habitat, and both increases and decreases to ecosystem productivity. Within lakes, the effect of climate warming on secondary production in littoral and pelagic habitats remains unclear. The intersection of changing habitat productivity and warming water temperatures on salmonids is important for understanding how climate warming will impact mountain ecosystems. We develop and test a conceptual model that expresses how earlier ice break-up dates influence within lake habitat production, water temperatures and the habitat utilized by, resources obtained and behavior of salmonids in a mountain lake. We measured zoobenthic and zooplankton production from the littoral and pelagic habitats, thermal conditions, and the habitat use, resource use, and fitness of Brook Trout (Salvelinus fontinalis). We show that earlier ice break-up conditions created a "resource-rich" littoral-benthic habitat with increases in zoobenthic production compared to the pelagic habitat which decreased in zooplankton production. Despite the increases in littoral-benthic food resources, trout did not utilize littoral habitat or zoobenthic resources due to longer durations of warm water temperatures in the littoral zone. In addition, 87% of their resources were supported by the pelagic habitat during periods with earlier ice break-up when pelagic resources were least abundant. The decreased reliance on littoral-benthic resources during earlier ice break-up caused reduced fitness (mean reduction of 12 g) to trout. Our data show that changes to ice break-up drive multi-directional results for resource production within lake habitats and increase the duration of warmer water temperatures in food-rich littoral habitats. The increased duration of warmer littoral water temperatures reduces the use of energetically efficient habitats culminating in decreased trout fitness.

Comparative Quality and Trend of Remotely Sensed Phenology and Productivity Metrics across the Western United States

Vegetation phenology and productivity play a crucial role in surface energy balance, plant and animal distribution, and animal movement and habitat use and can be measured with remote sensing metrics including start of season (SOS), peak instantaneous rate of green-up date (PIRGd), peak of season (POS), end of season (EOS), and integrated vegetation indices. However, for most metrics, we do not yet understand the agreement of remotely sensed data products with near-surface observations. We also need summaries of changes over time, spatial distribution, variability, and consistency in remote sensing dataset metrics for vegetation timing and quality. We compare metrics from 10 leading remote sensing datasets against a network of PhenoCam near-surface cameras throughout the western United States from 2002 to 2014. Most phenology metrics representing a date (SOS, PIRGd, POS, and EOS), rather than a duration (length of spring, length of growing season), better agreed with near-surface metrics but results varied by dataset, metric, and land cover, with absolute value of mean bias ranging from 0.38 (PIRGd) to 37.92 days (EOS). Datasets had higher agreement with PhenoCam metrics in shrublands, grasslands, and deciduous forests than in evergreen forests. Phenology metrics had higher agreement than productivity metrics, aside from a few datasets in deciduous forests. Using two datasets covering the period 1982–2016 that best agreed with PhenoCam metrics, we analyzed changes over time to growing seasons. Both datasets exhibited substantial spatial heterogeneity in the direction of phenology trends. Variability of metrics increased over time in some areas, particularly in the Southwest. Approximately 60% of pixels had consistent trend direction between datasets for SOS, POS, and EOS, with the direction varying by location. In all ecoregions except Mediterranean California, EOS has become later. This study comprehensively compares remote sensing datasets across multiple growing season metrics and discusses considerations for applied users to inform their data choices.

Competition for safe real estate, not food, drives density-dependent juvenile survival in a large herbivore

Density-dependent competition for food reduces vital rates, with juvenile survival often the first to decline. A clear prediction of food-based, density-dependent competition for large herbivores is decreasing juvenile survival with increasing density. However, competition for enemy-free space could also be a significant mechanism for density dependence in territorial species. How juvenile survival is predicted to change across density depends critically on the nature of predator-prey dynamics and spatial overlap among predator and prey, especially in multiple-predator systems. Here, we used a management experiment that reduced densities of a generalist predator, coyotes, and specialist predator, mountain lions, over a 5-year period to test for spatial density dependence mediated by predation on juvenile mule deer in Idaho, USA. We tested the spatial density-dependence hypothesis by tracking the fate of 251 juvenile mule deer, estimating cause-specific mortality, and testing responses to changes in deer density and predator abundance. Overall juvenile mortality did not increase with deer density, but generalist coyote-caused mortality did, but not when coyote density was reduced experimentally. Mountain lion-caused mortality did not change with deer density in the reference area in contradiction of the food-based competition hypothesis, but declined in the treatment area, opposite to the pattern of coyotes. These observations clearly reject the food-based density-dependence hypothesis for juvenile mule deer. Instead, our results provide support for the spatial density-dependence hypothesis that competition for enemy-free space increases predation by generalist predators on juvenile large herbivores.

Bottom-up when it is not top-down: Predators and plants control biomass of grassland arthropods

We investigate where bottom-up and top-down control regulates ecological communities as a mechanism linking ecological gradients to the geography of consumer abundance and biomass. We use standardized surveys of 54 North American grasslands to test alternate hypotheses predicting 100-fold shifts in the biomass of four common grassland arthropod taxa-Auchenorrhyncha, sucking herbivores, Acrididae, chewing herbivores, Tettigoniidae, omnivores, and Araneae, predators. Bottom-up models predict that consumer biomass tracks plant quantity (e.g. productivity and standing biomass) and quality (nutrient content) and that ectotherm access to food increases with temperature. Each of the focal trophic groups responded differently to these drivers: the biomass of sucking herbivores and omnivores increased with plant biomass; that of chewing herbivores tracked plant quality; and predator biomass did not depend on plant quality, plant quantity or temperature. The Exploitation Ecosystem Hypothesis is a top-down hypothesis that predicts a shift from resource limitation of herbivores when plant production is low, to predator limitation when plant production is high. In grasslands where spider biomass was low, herbivore biomass increased with plant biomass, whereas bottom-up structuring was not evident when spiders were abundant. Furthermore, neither predator biomass nor trophic position (via stable isotope analysis) increased with plant biomass, suggesting predators themselves are top-down limited. Stable isotope analysis revealed that trophic position of the chewing herbivore and omnivore increased significantly with plant biomass, suggesting these groups increased scavenging and meat consumption in grasslands with higher carbohydrate availability. Taken together, our snapshot sampling documents gradients of food web structure across 54 grasslands, consistent with multiple hypotheses of bottom-up and top-down regulation.

Environmental Differences between Migratory and Resident Ungulates—Predicting Movement Strategies in Rocky Mountain Mule Deer (Odocoileus hemionus) with Remotely Sensed Plant Phenology, Snow, and Land Cover

Migration is a valuable life history strategy for many species because it enables individuals to exploit spatially and temporally variable resources. Globally, the prevalence of species’ migratory behavior is decreasing as individuals forgo migration to remain resident year-round, an effect hypothesized to result from anthropogenic changes to landscape dynamics. Efforts to conserve and restore migrations require an understanding of the ecological characteristics driving the behavioral tradeoff between migration and residence. We identified migratory and resident behaviors of 42 mule deer (Odocoileus hemionus) based on GPS locations and correlated their locations to remotely sensed indicators of forage quality, land cover, snow cover, and human land use. The model classified mule deer seasonal migratory and resident niches with an overall accuracy of 97.8% and cross-validated accuracy of 81.2%. The distance to development was the most important variable in discriminating in which environments these behaviors occur, with resident niche space most often closer to developed areas than migratory niches. Additionally, snow cover in December was important for discriminating summer migratory niches. This approach demonstrates the utility of niche analysis based on remotely sensed environmental datasets and provides empirical evidence of human land use impacts on large-scale wildlife migrations.