Ecology

Resource competition and parasite exposure both present common density-dependent fitness costs for wild animals. Because launching effective immune responses is costly in terms of resources, parasites fitness costs should be further exacerbated in high-density, resource-depleted areas. To disentangle these relationships, we related density, parasitism, and resource availability to survival and fecundity across lifespan in a long-term study of wild red deer. All fitness measures declined with a combination of parasite count, greater density, and reduced resource availability. Beyond these relationships, as expected, local density and resource scarcity exacerbated survival costs of parasitism in calves, effectively undermining tolerance of infection. However, these synergistic relationships faded in yearlings and then reversed in adults, likely through age-structured selection biases. These findings emphasize that the costs of parasites and resource scarcity can be synergistic and intertwined with density in wild populations, accentuating the value of incorporating resource competition when examining parasite-dependent population regulation.

O_LISpecies diversity potentially has a dual effect on communities: a generally positive effect on overall community biomass, reflecting the expression of species response and interaction traits, and a poorly characterised effect on mass-specific species contribution to ecosystem functions, reflecting the expression of their effect traits. Disentangling the effects of biodiversity on total biomass from those on effect trait expression would help settle a long-standing debate by clarifying how biodiversity relates to both facets of species effects on ecosystem functioning. C_LIO_LIFollowing the classical BEF approach, we calculate expected ecosystem function based on observed functioning in monoculture. We then derive a net biodiversity effect (NBE) and decompose it into four components: the classical complementarity and selection effects on total community biomass, and complementarity and selection effects on effect trait expression. The latter two reflect, respectively, a complementarity or facilitation in how effect traits influence the function, and how species with the highest potential for increasing the function become dominant in the community. C_LIO_LIWe illustrate this NBE decomposition with three ecosystem functions (nitrogen retention capacity, soil hydraulic conductivity improvement, and forage digestibility) measured in assembled communities under controlled experimental conditions of perennial grassland plants. Regarding nitrogen retention, we find a positive complementary effect via total biomass, but a negative biodiversity effect via effect trait expression. For hydraulic conductivity improvement, biodiversity effects are mostly mediated by total biomass. As for forage digestibility, we found a positive complementarity effect on trait expression, outweighed however by a negative selection effect. This analysis reveals how biodiversity may have contrasting effects on ecosystem functions via its impact on biomass and effect trait expression. C_LI SynthesisSeparating between the effect of biodiversity on plant community biomass and on effect trait expression at the community level is one important step towards understanding the pathways by which diverse plant communities drive ecosystem functioning.

Methods for estimating and modeling the long-term and short-term adult flight dynamics of the conifer silk moth Dendrolimus superans (Lepidoptera: Lasiocampidae) are examined. The analysis uses light trap adult catch data collected over 21 years, from 2005 to 2025. Three models of adult flight are considered: a flight-initiation model driven by weather factors, an autoregressive model of long-term catch dynamics, and a binary model of seasonal catch. For the flight-initiation model, we propose estimating the accumulated temperature sum ST from the date when the first derivative of the remote sensing vegetation index NDVI becomes positive until the date of the first adult capture of the season. ST is shown to be sufficiently stable across all years of observation, with flight each year beginning after this temperature sum is reached. The second model demonstrates that the long-term light trap catch time series is well described by a second-order autoregressive model AR(2), in which the catch of the current year depends on catches from the two preceding years. This long-term series is compared with a previously studied larval population density series of the Siberian silk moth; both are shown to be AR(2) series with similar coefficient values, which suggesting that adult catch data may serve as a proxy for absolute larval population density. In the third model, we describe the transition from absolute-scale seasonal catch dynamics (number of adults per day) to a binary scale (0, 1), where 0 denotes days on which no adults were attracted to the trap, and 1 denotes days on which at least one individual was captured. The seasonal absolute catch series is thereby transformed into a binary series of zeros and ones, and relationships between adjacent values in such a binary series are examined. A linear relationship between the absolute and binary seasonal dynamics series is demonstrated, making it possible to estimate absolute catches from binary catch values and to analyze seasonal flight in sparse pest populations. This potentially opens new avenues for understanding how outbreak populations function at chronically low density. Author summaryForest pests can cause catastrophic damage, yet predicting their outbreaks remains challenging. During periods of low population density, standard monitoring methods become labor-intensive and uninformative, while the transition to an outbreak often occurs unexpectedly. Using a 21-year dataset of adult Siberian silk moth (Dendrolimus superans) captures from light traps, we developed an approach combining three complementary models. First, we showed that moth flight begins upon reaching a specific temperature sum, with the starting point determined by NDVI vegetation index dynamics rather than a calendar date--making the forecast more ecologically relevant. Second, long-term adult population dynamics follow a second-order autoregressive model AR(2), matching the dynamics previously observed for larval populations. This establishes light trap data as a reliable proxy for absolute population density when ground surveys are impractical. Third, we introduced a method to analyze seasonal flight using binary data (presence/absence of moths per day), which we showed is linearly related to absolute abundance. This enables studying population dynamics during periods of extremely low density, when traditional methods fail. Our approach opens new possibilities for early warning systems to detect when a population risks transitioning from a latent state to an outbreak phase.

Habitat complexity (HC) in part determines the diversity, stability, and behavior of food webs and can influence predation according to a wide variety of functional relationships. Many aquatic species provide habitat complexity and are also consumed by other species (e.g., macrophytes, corals, mussels). However, food web theory does not readily account for these species that act as edible habitat complexity (EHC). Here, we combine existing theory on predator-prey interactions, HC, and prey switching to describe the role of EHC in benthic food web models. We dissect feedback loops in each model to demonstrate how self-regulation of the prey species, mediated by species densities and HC, drives that food webs behavior. HC can stabilize predator-prey interactions by coupling prey self-regulation with HC self-regulation. EHC can further stabilize predator-prey interactions across a wide variety of "HC functions" that relate HC to predation rates. Significance StatementHabitat complexity (HC) plays a critical role in trophic interactions, population dynamics, and food web stability. However, little theory exists to describe edible habitat complexity (EHC), where a species is both consumed and confers habitat complexity for other species. We provide a series of models demonstrating how HC and EHC alter the population dynamics and stability of simple aquatic food webs. HC is strongly stabilizing in food webs by providing safety in rarity for prey. EHC provides safety in rarity for both prey and the EHC species because their predators are omnivorous. Given the prevalence of EHC species in aquatic systems (e.g., macrophytes, corals, mussels), our models demonstrate the importance of maintaining EHC species in aquatic systems for stable food webs.

Deadwood is a key habitat for forest biodiversity, yet how tree species and deadwood type shape linked fungal-beetle communities remain poorly understood. We explored saproxylic fungi and beetles in a large-scale restoration experiment on birch, pine, and spruce deadwood created as burned standing trees, felled logs, girdled trees, high stumps, and uprooted trees. As expected, we found that tree species was the main driver of both fungal and beetle community composition, while deadwood type was the second most important driver. Fungal-beetle community correlations were context dependent: significant multivariate correlations were detected for pine and spruce, but not birch, and were strongest in burned standing pine, burned standing spruce, and girdled spruce. Across all tree species and deadwood types, fungal-beetle co-occurrence networks were consistently less nested and more modular than expected by chance, indicating structured, compartmentalized associations of fungi and beetles even within single deadwood units. SynthesisThese results show that maintaining diverse tree species and deadwood types is essential to retain specialized multitrophic communities and the ecological processes they support.

O_LIUnderstanding functional community structure and the niche-based mechanisms that enable coexistence among sympatric species is essential for explaining how biodiversity is maintained in natural systems, and for anticipating how ecological communities will respond to ongoing environmental change. Stable isotope analysis provides a process-oriented perspective on resource use by integrating information across time and space, thereby allowing reconstruction of realised isotopic niches that reflect multiple dimensions of ecological differentiation. C_LIO_LIWe applied this framework to a community of ungulates in the Central-Eastern Italian Alps, including red deer (Cervus elaphus), roe deer (Capreolus capreolus), and Alpine chamois (Rupicapra rupicapra). Using stable isotope ratios in summer-grown hair segments ({delta}13C, {delta}15N, {delta}34S, {delta}18O, {delta}2H), we quantified species-specific n-dimensional niche hypervolumes within a Bayesian framework and estimated niche regions, overlap probabilities, univariate differentiation and multivariate structure. C_LIO_LIDespite broad dietary overlap typically observed among these ungulates, we found clear isotopic niche segregation, with mean pairwise overlap consistently remaining below 40%. Three dimensions emerged as primary drivers of differentiation: water sourcing ({delta}18O), diet quality ({delta}15N), and habitat openness ({delta}13C). Specifically, chamois appeared to derive more water from plants in their diet rather than from drinking, and to consume a higher-quality diet compared to Cervids. Red deer relied more heavily on forested habitats for resource use compared to roe deer and chamois, and additional isotopic differences between red deer and roe deer may stem from fine-scale abiotic conditions like microclimate and topography. We found no isotopic evidence for differential niche breadth among the three ungulate species. C_LIO_LITogether, these patterns highlight functional differentiation across multiple ecological axes, offering mechanistic insight into how these ungulates segregate realised niche space despite substantial potential for resource overlap. This multi-element isotope perspective underscores the value of integrative, process-based approaches for understanding current coexistence as well as improving predictions of how mammal communities may reorganise under accelerating environmental change. C_LI

Resource availability is a central driver of ecological and evolutionary processes, yet its effects on infectious disease and virulence are not fully understood. A key limitation is that many studies consider only a narrow range of resource conditions or a limited subset of host and pathogen traits, potentially obscuring non-linear relationships. Here, we quantify how a gradient of six food levels simultaneously shapes host fitness and pathogen performance in the Daphnia magna- Ordospora colligata system. Across two laboratory experiments, we measured infection rates, pathogen burden, host fecundity, survival, and filtration rates. Increased food availability enhanced pathogen fitness, with both infection rates and spore burden increasing with provisioning. In contrast, host responses were trait-specific: while fecundity increased with food availability, pathogen-induced reductions in fecundity (i.e., virulence) peaked at intermediate resource levels, despite continued increases in pathogen load. This pattern indicates that resource availability alters host tolerance as well as pathogen growth, generating non-linear disease outcomes. Host survival was unaffected by either food provisioning or infection, further demonstrating that resource availability can simultaneously influence host and pathogen traits in different directions. Our results highlight the importance of integrating multiple fitness components across provisioning levels to understand disease dynamics and suggest that ongoing anthropogenic changes in resource availability may alter host-pathogen interactions.

The phenology of organisms worldwide is shifting in response to changes in environmental conditions. There is growing concern that resulting timing mismatches among interacting species will negatively impact system-level properties, yet there is no general framework for evaluating community responses to changes in phenology. To address this gap, we developed a mathematical framework based on local stability analysis and used it to assess the resilience implications of phenological perturbations with a multi-year, highly time-resolved empirical dataset on subalpine plant-pollinator communities. The forecasted effects of phenological perturbations were largely independent of perturbations to species densities, indicating the potential for even small changes in phenology to disrupt the functioning of ecosystems that are otherwise highly stable.

Human activities are driving unprecedented environmental change, yet assessments of ecosystem resilience often overlook the rapid pace of change in the Anthropocene. Predator-prey systems are sensitive to the rate of environmental change and the whole system can collapse if predator population fail to promptly adjust to environmentally-driven shifts in resource population. Here, we investigate how different combinations of predator responsiveness and rates of environmental change influence the system vulnerability to critical transitions, explicitly addressing its interplay with magnitude of change. We found that, as predator responsiveness decreases, relatively slower rates and smaller magnitudes of environmental change leads to system collapse. Hence, even low and seemingly inoffensive total magnitudes of environmental change can be catastrophic if the rate of change is beyond a critical threshold. We propose considering predator responsiveness and current rates of environmental change as crucial factors in predicting the Anthropocenes impact on ecosystems.

Species tolerating the same environmental conditions can potentially colonize and thrive in the same habitats and eco-regions. Are any pair of those species equally probable to co-occur in the same community? Can we quantify the propensity of two species to co-occur together? Here, we focus on a simple but largely overlooked community-level pattern: the co-occurrence-occupancy curve, which relates the tendency of species to co-occur with others to their total occupancy across sites. We first define this empirical curve and then derive its expected shape under a random null model that assumes site equivalence and species independence. Building on these results, we introduce the Species Association Index (SAI), an occupancy-standardized measure that quantifies the tendency of a species to associate with others independently of its overall frequency of occurrence. The SAI enables meaningful comparisons among species with contrasting occupancies and provides a transparent benchmark against which departures from neutrality can be assessed. We illustrate the approach using two contrasting systems--tropical rain forest trees on Barro Colorado Island and organisms from Mediterranean rocky shores--highlighting both the generality of the co-occurrence-occupancy framework and its limitations.

Theory suggests that different components of environmental fluctuations, from daily and seasonal cycles to multidecadal trends, can have distinct and even opposing effects on species abundances and community dynamics, depending on their specific adaptations. But empirical research that deconstructs the influence of these different cycles on communities is lacking. Here, we used long-term biological monitoring data together with flow records of rivers across New Zealand to (i) investigate the role of fast, slow, and seasonal river-flow fluctuations in structuring macroinvertebrate communities; and (ii) to assess whether life-history and mobility traits mediate the response. Using joint species distribution models, we found striking differences in taxon and community responses to the different components of river flow variation. Responses to slow fluctuations were generally stronger and better predicted by traits, while responses to seasonal fluctuations were highly heterogeneous. Fast increases in flow, typical of flooding events, had pervasive negative effects on species abundances, but the severity of impact partly depended on mobility traits. Our results suggest that different ecological mechanisms underpin the response to distinct environmental fluctuations, highlighting the value of jointly considering multiple temporal scales of variation and species functional traits to understand and predict how communities reorganise under fluctuating environmental regimes.

1. Biotic insular systems differ from conventional islands because patch attributes change dynamically as patch-forming organisms develop. It therefore remains unclear whether the assembly mechanisms predicted by island biogeography theory (IBT) operate in such systems. Here, using epiphytic birds nest ferns (BNFs, Asplenium nidus) as a model biotic island system, we tested whether fungal and bacterial community diversity conform to species-area relationships predicted by IBT. With a stratified sampling scheme, we further evaluated the underlying mechanisms (passive sampling, disproportionate effects, and environmental heterogeneity) of species-area relationships, and assessed isolation effects using distance-decay patterns in community similarity. 2. We treated each BNF individual as a microbial island and categorized 24 BNFs into three size classes. Microbial and humus samples from multiple litter layers within each BNF individual were collected; microbial communities were characterized using next-generation sequencing, and humus chemical properties (pH and C:N ratio) were measured to characterize microhabitat conditions. To investigate mechanisms underlying species-area relationships, we applied a multi-scale rarefaction framework to partition diversity components. Spatial distances among BNFs were quantified to evaluate isolation effects. 3. Consistent with IBT predictions, both fungal and bacterial communities exhibited positive species-area relationships, indicating that larger BNFs harbored greater microbial richness. Diversity partitioning suggested that fungal richness increased through both disproportionate effects and environmental heterogeneity, whereas bacterial richness was primarily driven by environmental heterogeneity. Within larger ferns, greater heterogeneity in litter pH was associated with increased species turnover across litter layers, suggesting that decomposition-driven pH gradients create diverse microhabitats that promote microbial diversity. In addition, both microbial communities exhibited distance-decay patterns, indicating that isolation contributes to community assembly through dispersal limitation. 4. Synthesis. Our results demonstrate that BNFs function as a biotic insular system, in which both patch size and spatial isolation structure microbial diversity, consistent with predictions from IBT. Furthermore, we show that environmental heterogeneity generated by the growth of the habitatforming BNF mechanistically links island area to microbial diversity. Our study integrates both local habitat heterogeneity and regional spatial structure, highlighting the potential to extend IBT and metacommunity theory to organism-formed habitats.

Fuel availability and fuel dryness are consistently shown to be the primary drivers of wildfire intensity and burnt area. Here we hypothesise that differences in the timing of fuel build up and drying determine the optimal time for wildfire occurrence. We use gross primary production (GPP) as a measure of biomass production and hence fuel availability, and vapour pressure deficit (VPD) as a measure of fuel drying. We use the phase difference in the seasonal time course and magnitude of GPP and VPD to cluster regions that should therefore have distinct wildfire behaviour. We then show that each of the resultant clusters is distinctive in terms of one or more fire properties, specifically number of ignitions, burnt area, size, speed, duration, intensity, and length of the wildfire season. The emergence of distinct regimes as a function of two biophysical drivers reflects the fact that both vegetation and wildfire properties are a consequence of eco-evolutionary adaptions to environmental conditions. We then examine the degree to which human activities or vegetation properties modify these fire regimes within each of these clusters. Variability in GPP and VPD largely explains the within-cluster variation in fire properties. The type of vegetation cover has an influence on burnt area and carbon emissions in particular, while human activities are more important for fire properties such as size, rate of spread and duration largely through their influence of landscape fragmentation. Although both human activities and vegetation properties modify wildfire regimes, the ability to distinguish wildfire regimes using GPP and VPD alone emphasizes that land management, fire use and fire suppression are constrained by environmental conditions. This eco-evolutionary optimality approach to characterising wildfire regimes provides a basis for designing a simple fire model for Earth System modelling.

O_LIAlthough ecological research has long focused on the effects of temperature on population growth, arthropod pests are exposed to a wide variety of environmental factors that affect their performance, such as chemical pesticides targeted against them. Moreover, these environmental factors likely do not act in isolation. Identifying the extent to which abiotic factors interact to affect pest population dynamics can strengthen current and future pest management programs. C_LIO_LIHere, we investigated the extent to which temephos, a common pesticide applied to aquatic environments for mosquito control, influences the thermal performance of juvenile survival and development rate, as well as the intrinsic population growth rate, of the invasive mosquito pest, Aedes aegypti. We implemented a response surface experimental design to measure these traits across seven temperatures and five temephos concentrations and fit temperature- and insecticide-dependent performance curves to assess impacts on the overall performance and the thermal optimum, minimum, and maximum. C_LIO_LITemephos exposure profoundly altered the thermal performance of juvenile survival by reducing survival across all temperatures, shrinking the thermal breadth, and shifting the thermal optimum to warmer temperatures. Through this, temephos also altered the thermal performance of population growth primarily by reducing its thermal breadth. C_LIO_LISynthesis and applications: Our findings demonstrate that interactions between temperature and insecticide exposure can fundamentally reshape pest population dynamics, rather than acting as independent stressors. By quantifying this interaction, we showed that temphos is most effective below the pests thermal optimum, suggesting that larvicides may yield the greatest population suppression in cooler regions or during cooler periods of the year. Incorporating such temperature-dependent efficacy into pest management strategies could improve the timing and spatial targeting of control efforts. More broadly, these results highlight the need to integrate anthropogenic stressors with climatic drivers when predicting pest risk and optimizing management under ongoing environmental change. C_LI

O_LIPhenological shifts are a major ecological consequence of climate change, yet studies often focus on single life stages meaning that the potential for carryover effects between life stages remains poorly understood. Failing to account for these effects may lead to inaccurate estimates of phenological shifts, with consequences for predicted synchrony among interacting species. This is especially relevant for temperate systems where climate warming is occurring unevenly across the year. C_LIO_LIHere, we investigated how temperature experienced the previous autumn and winter (during the pupal and egg stage) influences spring phenology in the winter moth (Operophtera brumata), a herbivorous insect with distinct life stages. Using 50 years of local climate data to create five experimental temperature regimes, we first quantified phenotypic plasticity in the duration and temporal variability of pupal and egg development. We then examined how timing of adult moth emergence affects timing of offspring hatching. C_LIO_LIWe found divergent effects of temperature on different life stages; pupal development time was shortest at intermediate temperatures while egg development time decreased linearly with increasing temperature. Furthermore, phenological shifts due to the conditions experienced by the mother were carried over to influence the phenology of her offspring. While this carryover effect was partially compensated during subsequent stages, compensation decreased under warming conditions. C_LIO_LIThese results refine our understanding of the sensitivity of the annual cycle of winter moth phenology to variation in temperature with potential implications for population dynamics and interspecific interactions. Overall, our findings highlight the need to consider the impacts of warming across multiple life stages so that carryover effects can be properly accounted for. Doing so will improve predictions of phenological shifts under future climates. C_LI

Migration is widely recognized as a strategy for animals to track seasonally shifting resources. Yet, seasonal and spatial dynamics of migration are challenging to study, particularly for difficult-to-track insects. Among insects, monarch butterflies (Danaus plexippus) have a well-documented fall migration, but spring breeding recolonization remains poorly understood, particularly for the western population. We conducted multi-year surveys across six regions in the western United States to characterize monarch breeding phenology and evaluate three related hypotheses: (i) the successive broods model, with discrete generations shifting activity across the breeding range, (ii) a diffusion-like expansion model with overlapping breeding periods, and (iii) a mid-summer lull model with temporary summer declines in breeding for areas near the overwintering habitat. Monarch immature presence served as an indicator of local breeding activity. Our results do not support the successive broods or mid-summer lull hypotheses. Breeding onset occurred earlier near overwintering areas and gradually expanded north-and eastward, with sustained activity in many regions throughout the season. Termination of breeding also occurred earlier at more distant sites, resulting in longer breeding activity nearer to overwintering habitat. Immature monarch density declined with distance from overwintering areas at onset and termination, suggesting delayed colonization of peripheral regions. Together, these results support a diffusion-like expansion of breeding rather than sequential generational replacement. Western monarchs also do not initiate or terminate breeding in close synchrony with host plant availability, contrary to predictions from resource-tracking theory. These findings highlight fundamental differences between western monarch breeding dynamics and paradigms for eastern monarchs, demonstrating that a single species can employ fundamentally different spatial strategies for recolonizing its breeding range in different regions. More generally, these results distinguish insect migration from systems with direct movements between wintering and breeding habitats, and underscore the value of long-term, landscape-scale monitoring for resolving habitat use across heterogeneous environments.

Due to long co-evolutionary histories, many zoonotic pathogens are thought to exert little or no negative effects on their wildlife reservoir hosts. However, there remains a lack of rigorous investigations in natural settings. We conducted a 3-year factorial field experiment to investigate how survival of the Puumala hantavirus (PUUV) reservoir, the bank vole, is impacted by PUUV infection, nematode infections, and food availability. We hypothesized that PUUV would not impact survival, but that coinfection with nematodes would negatively impact survival, and that increased food availability would mitigate the negative effects of coinfection. Surprisingly, we demonstrated that PUUV infected voles had substantially reduced survival when compared to uninfected voles, and this strong negative effect manifested in young voles. Nematode removal increased survival of young voles and food supplementation interacted with movement rather than survival. Our results provide empirical evidence in a natural system for infection reducing survival of its reservoir host.

Predicting the stability of ecological communities in changing environments is challenging. Classical theory posits that community stability cannot be understood without considering interspecific interactions. A contrasting view is that species environmental responses and their variation (response diversity) influence stability to the extent that effects of interspecific interactions can be ignored. Surprisingly, few studies have evaluated the relative importance of interactions versus species responses. Moreover, trait-based measures of response diversity often show limited predictability. Here, we introduce community performance curves, the aggregate of species performance curves, as a powerful mechanistic link between community composition and stability. This approach reveals that species responses predict most of the variation in community stability in simulated communities, even when the strength of interspecific interactions varies. An experiment with ciliate communities corroborates these findings, while a literature review reveals how rarely both mechanisms are assessed jointly. By moving from summary traits to community performance curves, we reconcile the two perspectives: while species interactions undeniably shape community dynamics, community performance curves are sufficient to predict stability. This provides the opportunity to predict community stability, even when information about the multitude and diversity of interspecific interactions is unavailable.

O_LIChanges in habitat structure across species distributions may contribute to the generation and maintenance of range limits, but few studies have evaluated this by directly measuring habitat availability across relevant spatial scales. C_LIO_LIHere, we test the predictions that coarse-scale and patch-level habitat availability decline towards and beyond the northern range limit of Pacific coastal dune endemic Camissoniopsis cheiranthifolia. We used aerial imagery and geographic information system (GIS) tools to measure the coarse-scale availability of coastal dune habitat in California and Oregon. The availability of finer-scale habitat patches specifically suitable for C. cheiranthifolia was measured in a 2-generation field survey of > 4,200 5m x 5m plots randomly distributed across 1100 km of coastal dune habitat transcending the species northern range limit. At each plot, we estimated the proportion of area that contained suitable habitat as well as recorded occupancy by C. cheiranthifolia. As an alternative approach to visually estimating habitat suitability, we recorded plant community composition at each plot to predict beyond-range habitat suitability using a random forest model. C_LIO_LIContrary to our predictions, we found that coastal dune habitat, measured coarsely from aerial imagery, was more abundant and continuous towards and beyond the northern range limit. At the fine scale, however, the proportion of plots with suitable habitat (patch suitability) and the proportion of habitat within plots that was suitable (patch size) declined across the range limit. Moreover, patches were more isolated from one another and, in one survey year, less temporally stable towards and beyond the range limit. Finally, occupancy by C. cheiranthifolia was less likely in smaller, more isolated, and temporally unstable patches, providing mechanistic insight to the previously observed decline in occupancy towards the range limit. C_LIO_LISynthesis: Taken together, our results suggest that fine-scale habitat patch configuration changes in ways that likely impede patch colonization, thereby reducing occupancy and limiting the species northern distribution. Thus, consideration of geographic variation in patch and landscape structure, rather than only coarse-scale habitat availability, may be essential for understanding the processes that limit species ranges. C_LI

Microbial communities can shift under drought in ways that enhance plant performance during drought ("microbe-mediated acclimation"). However, it is also possible for microbial communities to shift in ways that worsen the effects of drought ("mal-acclimation"). It is unclear how and where microbe-mediated acclimation vs. mal-acclimation occurs, or if there are types of soils or microbial communities that are more likely to harbor microbes that enhance plant acclimation and limit mal-acclimation. We tested for microbe-mediated plant acclimation/mal-acclimation to drought in soils from 21 maize farms in the midwestern United States, spanning a range of climate, soil types, and management practices. We first conditioned soil microbial communities to drought vs. well-watered conditions in a greenhouse and then tested for microbe-mediated acclimation by growing maize in soils inoculated with the conditioned microbial communities under drought and well-watered conditions. Drought-conditioned soils did not enhance plant performance under drought. In fact, one third of the farms exhibited mal-acclimation, especially under well-watered conditions where wet-conditioned soils reduced plant performance in well-watered contemporary conditions. Farm management practices, climate, soil texture, and microbial diversity generally did not predict when this microbe-mediated mal-acclimation occurred. Overall, these results suggest that in agricultural soils, microbes may frequently impede-rather than facilitate-plant acclimation to soil moisture levels. Open research statementThe plant and soil data used in this study are available via the Environmental Data Initiative repository at https://doi.org/10.6073/pasta/f4a0db3a076cf6d8cef908947f82736e. The bacterial and fungal amplicon sequence data are available via the European Nucleotide Archive under accessions PRJEB110071 and PRJEB109827, respectively.