Relationship between depth, sediment, latitude, and the structure of benthic infaunal assemblages on the mainland shelf of southern California

A regional benthic survey was conducted in 1994, and the data were used to assess the relationship among three habitat factors (depth, sediment grain size, and latitude) and the distribution of benthic infaunal assemblages on the southern California coastal shelf. Benthic samples were collected with a 0.1 m² Van Veen grab from 251 sites on the continental shelf (10–200 m deep) from Point Conception, California, to the United States–Mexico international border. The relationship between habitat and assemblages was investigated by conducting a Q-mode cluster analysis to define groups of stations with similar species composition and then examining whether differences were present in physical habitat attributes among those groups of stations. Analysis of data from 175 uncontaminated sites yielded four habitat-related benthic infaunal assemblages along the southern California coastal shelf: a shallow-water assemblage from 10–32 m, a mid-depth assemblage between 32 and 115 m, and two deep-water (115–200 m) assemblages, one in fine and one in coarse sediments. These empirically defined points in the depth and sediment grain size gradients can be used to define reference habitats for the development of biocriteria. Benthic abundance and diversity were greatest in the mid-depth assemblage, conforming to predictions for benthic assemblages in regions of upwelling. Within the 500 km of coastline examined, latitude was not an important factor in defining assemblages. Received: 3 December 1999 / Accepted: 9 October 2000     

Understanding multiple ecological responses to anthropogenic disturbance: rivers and potential flow regime change

Human-induced alteration of the natural flow regime is a major threat to freshwater ecosystems and biodiversity. The effects of hydrological alteration on the structural and functional attributes of riverine communities are expected to be multiple and complex, and they may not be described easily by a single model. Based on existing knowledge of key hydrological and ecological attributes, we explored potential effects of a flow-regulation scenario on macroinvertebrate assemblage composition and diversity in two river systems in Australia's relatively undeveloped wet-dry tropics. We used a single Bayesian belief network (BBN) to model potential changes in multiple assemblage attributes within each river type during dry and wet seasons given two flow scenarios: the current, near-natural flow condition, and flow regulation. We then used multidimensional scaling (MDS) ordination to visually summarize and compare the most probable attributes of assemblages and their environment under the different scenarios. The flow-regulation scenario provided less certainty in the ecological responses of one river type during the dry season, which reduced the ability to make predictions from the BBN outputs directly. However, visualizing the BBN results in an ordination highlighted similarities and differences between the scenarios that may have been otherwise difficult to ascertain. In particular, the MDS showed that flow regulation would reduce the seasonal differentiation in hydrology and assemblage characteristics that is expected under the current low level of development. Our approach may have wider application in understanding ecosystem responses to different river management practices and should be transferred easily to other ecosystems or biotic assemblages to provide researchers, managers, and decision makers an enhanced understanding of ecological responses to potential anthropogenic disturbance.     

The cost of tolerance: sensitivity of stream benthic communities to UV-B and metals.

The ability to tolerate disturbance is a defense strategy that minimizes the effects of damage to fitness and is essential for sustainability of populations, communities, and ecosystems. Despite the apparent benefits of tolerance, there may be an associated cost that results in a deficiency of a system to respond to additional disturbances. Aquatic ecosystems are often exposed to a variety of natural and anthropogenic disturbances, and the effects of these compound perturbations are not well known. In this investigation, we examine whether tolerance to one stressor, metals, results in a cost of increased sensitivity to an additional stressor, ultraviolet-B (UV-B) radiation. Heavy metal pollution is recognized as a major environmental problem in Rocky Mountain streams. These high-elevation, typically clear streams may be at particular risk to elevated UV-B levels associated with reduced levels of ozone. Microcosm experiments were conducted using natural stream benthic communities collected from a reference site and a site with a long-term history of heavy-metal pollution. Direct and interactive effects of heavy metals and UV-B radiation on structural and functional characteristics of benthic communities were evaluated among four treatments: control, UV-B, metals, and metal and UV-B. Communities from the metal-polluted site were more tolerant of metals but less tolerant to UV-B compared to reference communities. Increased mayfly drift and reduced metabolism in response to metal exposure were observed in reference communities but not in the metal-polluted communities. In contrast to these results, UV-B radiation significantly reduced community metabolism, total macroinvertebrate abundance, and abundances of mayflies, caddisflies, and dipterans from the metal-polluted site, but had no effects on benthic communities from the reference site. ANOSIM results demonstrated that community responses differed among treatments at both sites. Metals had the largest impact on community differences at both sites, while UV-B had greater impacts at the metal-polluted site. This research demonstrates the need to account for potential costs associated with tolerance and that these costs can result in behavioral, structural, and functional impacts to benthic communities.     

The impact of seabed disturbance on nematode communities: linking field and laboratory observations

Physical disturbance is a key factor in controlling the spatial and temporal composition of shallow-water benthic communities. Like shallow waters, deeper waters are increasingly subject to a range of anthropogenic disturbances, which can lead to significant alterations in sedimentation patterns. These alterations often exceed naturally occurring changes. We used a combined analysis of six independent data sets arising from large-scale field surveys and small-scale laboratory experiments to investigate the effects of seabed disturbance on nematode communities. Disturbance response was documented as a function of disturbance type (coastal development, dredged material disposal, bottom trawling, glacial fjord) and intensity (low, medium, high). Natural and man-induced seabed disturbance exerted differential effects on exposed populations, generating changes in the taxonomic (genus) and functional (feeding type) attributes of their assemblages. The genus composition of nematode assemblages from geographically separate seas converged with increased level of various types of man-made disturbance. Assemblages present along a gradient of natural disturbance in a glacial fjord followed an opposite response vector, suggesting that community changes induced by anthropogenic activities, or experimental treatments simulating the principal impacts of these, inherently differ from disturbance of natural origin. Changes in trophic diversity and structure were primarily driven by factors confounded with physical disturbance, such as metal contamination. Coupling the results of analyses at multiple scales proved a useful means of providing deeper insights into the general response of ecological communities to environmental change.     

Benthic diversity gradients and shifting baselines: implications for assessing environmental status

The increasing pressure on marine biodiversity emphasizes the importance of finding benchmarks against which to assess change. This is, however, a notoriously difficult task in estuarine ecosystems, where environmental gradients are steep, and where benthic biodiversity is highly variable in space and time. Although recent emphasis on diverse, healthy benthic communities in legislative frameworks has increased the number of indices developed for assessing benthic status, there is a lack of quantitative baselines in benthic diversity that would enable comparisons across broad spatial scales, encompassing different environmental settings and bioregions. By taking advantage of long-term monitoring data, spanning hundreds of stations over the past 40 years, we provide a comprehensive analysis of benthic a, beta, and gamma diversity, encompassing the entire' salinity gradient of the open sea areas of the large, brackish-water Baltic Sea. Using a relatively simple measure, average regional diversity, we define area-specific reference conditions and acceptable deviation against which to gauge current conditions in benthic macrofaunal diversity. Results show a severely impaired condition throughout large areas of the Baltic for the assessment period 2001-2006. All ecosystems are plagued by baselines that shift in time and space, and their definition is not trivial, but average regional diversity may offer a transparent way to deal with such changes in low-diversity systems. Identifying baselines will be of increasing importance given the potential of climatic drivers to interact with local anthropogenic stressors to affect patterns of biodiversity. Our analysis provides an evaluation of the current condition in a system that has been heavily influenced by anthropogenic impact and changing oceanographic conditions, and it provides a basis for future impact assessment and ecosystem-based management.     

Fish assemblages composition in a natural, then regulated, stream: A quantitative long-term study

This paper tests the hypothesis how effectively fish assemblage composition was shaped by local climate changes and by engineering impact in 1989. This was possible due to the monitoring protocols of the present study, which allowed estimation of the magnitude of anthropogenic changes from changes naturally occurring in nature. Fish were sampled at the end of every growing season (October) for 23 years at five contiguous sites in a stream, before (1979-1988) and after (1989-2001) regulation. In each sample, six successive electrofishing passes were used to calculate the density and mean biomass for assemblage analysis using the Zippin model. During the study, the natural, meandering stream with pools, riffles, and a moderate canopy was modified into a straight stream of uniform width and depth, stripped of all vegetation. The output layer of a self-organizing map (SOM, the artificial neural network algorithm) applied in this study for site similarity analysis was partitioned into six subclusters placed in two main clusters. Subclusters in the upper part of the SOM were occupied chiefly by regulated stream samples and those in the lower part of the SOM by natural stream samples. Subclusters in the middle position, contained both natural (19) and regulated (20) samples in nearly equal proportion. In addition, the SOM contained one subcluster with sites only from the regulated period and another with only natural sites. Differences between subclusters were attributed to differences in climate, with some differences profound. Warming of the local weather, which became most evident in the 1990s, may have resulted in changes in fish assemblages. This is shown in the SOM, in which samples from the 1980s with cold years dominate the bottom of the SOM, whereas those from the 1990s and later are at the top. Subclusters dominated by regulated or natural sites were not always significantly different when the number of species and diversity indices were considered. Clear differences between the regulated and natural samples involved qualitative characteristics and mainly concerned assemblage composition. They were also confirmed by significant IndVal values (indicator species) and neither mixed subcluster contained important species in their assemblages.     

Motile cryptofauna associated with live and dead coral substrates: implications for coral mortality and framework erosion

Coral reef cryptofauna are a diverse group of metazoan taxa that live within intra- and inter-skeletal voids formed by framework structures. Despite a hypothesized high biomass and numerous trophic roles, they remain uncharacterized relative to exposed reef communities. Motile cryptofauna were sampled from live coral colonies and dead frameworks typifying four successive levels of degradation on an eastern Pacific pocilloporid reef. Abundances and biomass were higher on live versus dead corals habitats. The density of cryptofauna per volume substrate was highest on dead coral frameworks of intermediate degradation, where complex eroded substrates provide abundant shelters. These data have important and far-reaching ramifications for how the diverse multispecies assemblages that are reef ecosystems will respond to anthropogenic stressors such as those associated with climate change. Extreme levels of coral mortality, bioerosion, and habitat destruction will lead to impairment and eventually loss of ecosystem functions.     

Responses of marine benthic microalgae to elevated CO2

Increasing anthropogenic CO₂ emissions to the atmosphere are causing a rise in pCO₂ concentrations in the ocean surface and lowering pH. To predict the effects of these changes, we need to improve our understanding of the responses of marine primary producers since these drive biogeochemical cycles and profoundly affect the structure and function of benthic habitats. The effects of increasing CO₂ levels on the colonisation of artificial substrata by microalgal assemblages (periphyton) were examined across a CO₂ gradient off the volcanic island of Vulcano (NE Sicily). We show that periphyton communities altered significantly as CO₂ concentrations increased. CO₂ enrichment caused significant increases in chlorophyll a concentrations and in diatom abundance although we did not detect any changes in cyanobacteria. SEM analysis revealed major shifts in diatom assemblage composition as CO₂ levels increased. The responses of benthic microalgae to rising anthropogenic CO₂ emissions are likely to have significant ecological ramifications for coastal systems.     

Spatial variability of infaunal nematode and polychaete assemblages in two muddy subtidal habitats

Spatial variation of populations and their assemblages is an important component of many aspects of ecology, including the maintenance of species diversity. In nature, organisms are generally aggregated in patches or form gradients or other spatially related patterns. This investigation quantified the degree of spatial structure in benthic invertebrate assemblages. It examined the distribution of infaunal assemblages of different size, mobility and contrasting life-history strategies (i.e. meiofaunal nematodes and macrofaunal polychaetes) in two offshore muddy habitats in the Celtic Deep and the NW Irish Sea off the west coast of the United Kingdom. The more heterogeneous habitat in the NW Irish Sea was characterised by higher tidal stress and bottom temperature while greater water depth, mean particle diameter and organic carbon content was typical for the comparatively homogeneous environment in the Celtic Deep. In both habitats, the environmental conditions became increasingly dissimilar with separation. A total of 125 nematode and 88 polychaete species were recorded with 69% of the nematode and 49% of polychaete species present at both study sites. Occurrence frequency of species, species diversity and average similarity of assemblage composition was higher in the Celtic Deep than in the NW Irish Sea. Results from correlation analyses revealed statistically robust relationships of community similarity and sample distance. Given their small size and low mobility, nematodes were more susceptible to within-habitat physical variability than larger-sized, more mobile polychaetes. This, coupled with limits to long-distance dispersal and likely restrictions in gene flow, resulted in a significant decrease in community similarity with distance at the spatial scales sampled (i.e. within 0.1–23 km). Polychaetes, in contrast, combined a higher dispersal potential (>23 km) with a relatively high tolerance to within-habitat environmental variability and these were the most likely causes for non-significant relationships between the similarity of their assemblages and sample distance. The potential mechanisms causing the observed variation and implications of results for environmental management strategies are discussed. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Are soil mite assemblages structured by the identity of native and invasive alien grasses?

Associations between plants and animals in aboveground communities are often predictable and specific. This has been exploited for the purposes of estimating the diversity of animal species based on the diversity of plant species. The introduction of invasive alien plants into an ecosystem can result in dramatic changes in both the native plant and animal assemblages. Few data exist at the species level to determine whether belowground animal assemblages share the same degree of association to plants. The hypotheses that soil mites (Acari) form assemblages specifically associated with different native grass species in an unmanipulated natural ecosystem and that invasive alien grasses will impact soil mite assemblage composition in this setting were tested. Soil mites sampled beneath five native and two invasive alien species of grasses at the Konza Prairie Biological Station, Kansas, USA, were similarly abundant, species rich, diverse, and taxonomically distinct. No mite species had affinities for a specific grass species. There was no evidence from analysis of similarity, canonical correspondence analysis, or a nonparametric assemblage analysis that the assemblage composition of soil mites was specific to grass species. Results suggest that soil mite assemblages were more related to characteristics of the plant assemblage as a whole or prevailing soil conditions. The most recent invasive alien grass did not support a successionally younger mite fauna, based on the ratio of mesostigmatid to oribatid mites, and neither of the two invasive grasses influenced mite assemblage structure, possibly because they had not yet substantially altered the soil environment. Our results suggest that extrapolations of soil mite diversity based on assumptions of plant specificity would be invalid.     

The application of predicted habitat models to investigate the spatial ecology of demersal fish assemblages

Benthic habitats are known to influence the abundance and richness of demersal fish assemblages; however, little is known about how habitat structure and composition influences these distributions at very fine scales. We examined how the benthic environment structures marine fish assemblages using high-resolution bathymetry and accurate predicted benthic habitat maps. Areas characterised by a mosaic of habitat patches supported the highest richness of demersal fishes. A total of 37.4% of the variation in the distribution of the fish assemblage was attributed to 6 significant variables. Depth explained 23.0% of the variation, with the boulders explaining 12.6% and relief 1.4%. The remaining measures (seawhips, light/exposure and solid reef) provided a small (<1.0%) but significant contribution. Identifying components of the benthic environment important in structuring fish assemblages and understanding how they influence the spatial distribution of marine fishes is imperative for better management of demersal fish populations.     

Nutzen und Grenzen von Modellökosystemen zur Erfassung der Wirkung von Schadstoffen auf benthische Lebensgemeinschaften

The aim of this short article is to deliver insight into the prospects and constraints of model ecosystems in the assessment of effects of contaminants on benthic communities, especially the meiofauna. Numerous substances accumulate in sediments and thus a complex contamination exists which is able to impair benthic organisms. Besides direct effects, indirect effects are of importance and should be considered. Changes in competition and predator-prey-relationships lead to a variety of effects in the benthic food web that cannot be explained by direct toxicity of contaminants. In this article two microcosm studies, one with the heavy metal cadmium and one with the veterinary pharmaceutical Ivermectin, are presented which both showed direct as well as indirect effects. A significant change in abundance and composition of the meiofauna community has been observed. The prediction of field-effects of a certain substance as well as the interpretation of field data are the aims of model ecosystems. Both studies proofed that the investigation of meiofauna in microcosms is a promising tool to achieve these aims. Furthermore, microcosms with meiofauna might be able to assist the effect-directed analysis (EDA) in a risk assessment by identifying the effect-causing contaminants of complex environmental samples.     

Settlement seasonality and temporal changes in hard substrate macrozoobenthic communities of Lesina Lagoon (Apulia, Southern Adriatic Sea)

The macrofauna settling on experimental substrates was studied at two sites of the Lesina Lagoon to test its possible role in monitoring ecological variations in a brackish-water ecosystem. The community settlement was seasonally investigated on 3-month-old wooden poles; the development was monitored from 2001 to 2005. Comparisons with benthic assemblages settled on 10-yr-old poles were also performed. The main hydrological parameters were periodically measured during the study. A total of 38 species were collected. The assemblage reached the highest development in the central lagoon, showing relevant carbonate structures which supported a rich vagile fauna. Close to the sea-water inlet species richness and abundance values were lower, with the disappearance of some brackish-water species. A remarkable salinity drop during 2004 produced some faunistic changes in the assemblages, which however maintained different structures between the study sites, thus confirming macrozoobenthos as an efficient bioindicator of different environmental conditions for transition biotopes and a useful investigation tool in monitoring programmes.     

Sediment characteristics and colonization of soft-bottom benthos: a field manipulation experiment

A field manipulation experiment was carried out to test the effects of sediment characteristics (particle size and organic content) on colonization of soft-bottom benthos. Defaunated mud and sand in experimental trays were exposed at sub-tidal levels of two sites with sandy and muddy sea bottom, and retrieved monthly for examination of benthic composition. Seventy-eight out of the 107 species recorded, and 50% of the dominant species identified in the present study, were common in both the sand and mud trays. Classification analysis revealed higher faunal similarity of benthic composition in both sediment types within each of the two study sites. Results of three-way ANOVA (Sediment, Site and Time) showed that sediment effect was only significant for colonization of bivalves and gastropods, but not for polychaetes, amphipods, total species and individual numbers, or overall benthic diversity and evenness. This suggested that many of the benthic species can colonize both sediment types, and that factors other than sediment particle size and organic content may play a significant role in determining colonization of these species. The faunal composition of the same sediment type between sites was found to be different, indicating that the availability and abundance of benthic larvae/adults may be more important than sediment characteristics in determining benthic settlement. An increase in faunal similarity between mud and sand was observed over the study period, implying that the importance of sediment characteristics to benthic colonization decreased with exposure time. Received: 29 November 1996 / Accepted: 16 January 1997     

Lionfish alter benthic invertebrate assemblages in patch habitats of a subtropical estuary

Invasive lionfish (primarily Pterois volitans) have spread throughout the Caribbean region, Gulf of Mexico, and the Southeast US coast, and resulting impacts on reef fish populations have been well documented. We examined whether lionfish can likewise affect benthic invertebrate communities, using an in situ caging experiment, in the Loxahatchee Estuary, Florida, USA. We found that lionfish caused significant declines in the three most abundant benthic invertebrate species, driving an overall shift in assemblage composition. For example, grass shrimp (Palaemonidae) abundance was reduced by nearly 90 % in the presence of lionfish. Species richness of benthic organisms was significantly higher when lionfish were present, suggesting potentially complex emergent effects of lionfish predation on benthic assemblages. Despite the fact that this experiment was conducted in just a single location using relatively small experimental units, we show altered benthic invertebrate communities could well be an additional outcome of the lionfish invasion.     

Sewage discharges in oceanic islands: effects and recovery of eulittoral macrofauna assemblages

Sewage discharges are among the most common anthropogenic stressors on rocky shores. However, studies on benthic invertebrates’ responses to sewage pollution on insular assemblages are rare. In order to fulfil this gap, (i) the effects of sewage pollution and (ii) the recovery capacity of eulittoral insular assemblages were examined. The study was conducted on Terceira Island (Azorean Archipelago, Portugal), in four distinct areas: one impacted area (outfall in operation), one post-impacted area (decommissioned outfall) and two undisturbed reference areas. Results showed that the abundance of the dominant species changed between reference and impacted areas. More specifically, the abundance of limpets Patella candei, and barnacles Chthamalus stellatus decreased, while the abundance of littorinids Tectarius striatus and Melarhaphe neritoides increased near the sewage discharges. Seventeen months after removal of the outfall, limpet and barnacle populations had partially recovered, but littorinid populations had still not fully recovered. This study confirms the negative effects of sewage discharges on eulittoral assemblages, and provides information on the resilience of benthic species following the cessation of sewage disposal.     

Trophic structure of coastal Antarctic food webs associated with changes in sea ice and food supply

Predicting the dynamics of ecosystems requires an understanding of how trophic interactions respond to environmental change. In Antarctic marine ecosystems, food web dynamics are inextricably linked to sea ice conditions that affect the nature and magnitude of primary food sources available to higher trophic levels. Recent attention on the changing sea ice conditions in polar seas highlights the need to better understand how marine food webs respond to changes in such broad-scale environmental drivers. This study investigated the importance of sea ice and advected primary food sources to the structure of benthic food webs in coastal Antarctica. We compared the isotopic composition of several seafloor taxa (including primary producers and invertebrates with a variety of feeding modes) that are widely distributed in the Antarctic. We assessed shifts in the trophic role of numerically dominant benthic omnivores at five coastal Ross Sea locations. These locations vary in primary productivity and food availability, due to their different levels of sea ice cover, and proximity to polynyas and advected primary production. The delta15N signatures and isotope mixing model results for the bivalves Laternula elliptica and Adamussium colbecki and the urchin Sterechinus neumeyeri indicate a shift from consumption of a higher proportion of detritus at locations with more permanent sea ice in the south to more freshly produced algal material associated with proximity to ice-free water in the north and east. The detrital pathways utilized by many benthic species may act to dampen the impacts of large seasonal fluctuations in the availability of primary production. The limiting relationship between sea ice distribution and in situ primary productivity emphasizes the role of connectivity and spatial subsidies of organic matter in fueling the food web. Our results begin to provide a basis for predicting how benthic ecosystems will respond to changes in sea ice persistence and extent along environmental gradients in the high Antarctic.     

Temporal patterns in diversity and species composition of deep-sea benthic copepods in bathyal Sagami Bay,central Japan

Little is known about temporal changes in the diversity and species composition of deep-sea metazoan meiofauna and their relationships with changes in the food supply. Those changes were studied for benthic copepod assemblages based on 2-year time-series data at a bathyal site in Sagami Bay (1430 m depth), central Japan, where annual fluctuation in the abundance of benthic foraminiferans was previously observed. Species diversity of benthic copepods at the site was as high as, or slightly higher than, that observed at other deep-sea sites, but did not fluctuate temporally through the study period. Multivariate analyses did not reveal any clear seasonal or directional change occurring over the longer term in their species composition, although there was some consistent pattern. These results indicate a lack of, or only weak, seasonality in the diversity and species structure of the deep-sea benthic copepod assemblages, even though the fresh organic food supply fluctuates seasonally. They also suggest that there are differences between copepods and foraminiferans in the response to changes in environmental factors, and that spatial differences in the composition of copepod communities are greater than temporal ones at this deep-sea site.Communicated by T. Ikeda, Hakodate     

Testing a new classification of morphological functional groups of marine macroalgae for the detection of responses to stress

Morpho-functional group classifications have been widely used in studies concerning benthic macroalgal communities, due to their easy application and the reduced cost and workload compared with species-level identifications. However, their effectiveness in representing spatio-temporal patterns and describing the effect of disturbances has been disputed. We propose a new, expanded classification of morphological functional groups based on thallus structure, growth form, branching pattern, and taxonomic affinities. This approach was tested on subtidal macroalgal assemblages of the western Mediterranean. We examined the differences in composition of assemblages subjected to different types of stressors using species, traditional morphological groups, and new morphological groups. Results showed that the new morphological groups allowed to detect patterns highlighted through species analyses generally better than the use of traditional morpho-functional groups. The new morphological groups may represent a valid tool in monitoring programs and environmental impact assessments, where it is important to detect early stages of environmental changes using methods that allow to examine a large number of samples in a limited time.     

A landscape-scale framework to identify refugia from multiple stressors

From a conservation perspective, quantifying potential refugial capacity has been predominantly focused on climate refugia, which is critical for maintaining the persistence of species and ecosystems. However, protection from other stressors, such as human-induced changes in fire and hydrology, that cause habitat loss, degradation, and fragmentation is also necessary to ensure that conservation efforts focused on climate are not undermined by other threats. Thus, conceptual and methodological advances for quantifying potential refugia from multiple anthropogenic stressors are important to support conservation efforts. We devised a new conceptual approach, the domains of refugia, for assessing refugial capacity that identifies areas where exposure to multiple stressors is low. In our framework, patterns of environmental variability (e.g., increased frequency of warm summers), thresholds of resilience, and extent and intensity of stressors are used to identify areas of potential refugia from a suite of ongoing anthropogenic stressors (e.g., changes in fire regime). To demonstrate its utility, we applied the framework to a Southern California landscape. Sites with high refugial capacity (super-refugia sites) had on average 30% fewer extremely warm summers, 20% fewer fire events, 10% less exposure to altered river channels and riparian areas, and 50% fewer recreational trails than the surrounding landscape. Our results suggest that super-refugia sites (∼8200 km²) for some natural communities are underrepresented in the existing protected area network, a finding that can inform efforts to expand protected areas. Our case study highlights how considering exposure to multiple stressors can inform planning and practice to conserve biodiversity in a changing world.