skip to main content


Title: Plant species richness at archaeological sites suggests ecological legacy of Indigenous subsistence on the Colorado Plateau
Humans have both intentional and unintentional impacts on their environment, yet identifying the enduring ecological legacies of past small-scale societies remains difficult, and as such, evidence is sparse. The present study found evidence of an ecological legacy that persists today within an semiarid ecosystem of western North America. Specifically, the richness of ethnographically important plant species is strongly associated with archaeological complexity and ecological diversity at Puebloan sites in a region known as Bears Ears on the Colorado Plateau. A multivariate model including both environmental and archaeological predictors explains 88% of the variation in ethnographic species richness (ESR), with growing degree days and archaeological site complexity having the strongest effects. At least 31 plant species important to five tribal groups (Navajo, Hopi, Zuni, Ute Mountain Ute, and Apache), including the Four Corners potato ( Solanum jamesii ), goosefoot ( Chenopodium sp.), wolfberry ( Lycium pallidum ), and sumac ( Rhus trilobata ), occurred at archaeological sites, despite being uncommon across the wider landscape. Our results reveal a clear ecological legacy of past human behavior: even when holding environmental variables constant, ESR increases significantly as a function of past investment in habitation and subsistence. Consequently, we suggest that propagules of some species were transported and cultivated, intentionally or not, establishing populations that persist to this day. Ensuring persistence will require tribal input for conserving and restoring archaeo-ecosystems containing “high-priority” plant species, especially those held sacred as lifeway medicines. This transdisciplinary approach has important implications for resource management planning, especially in areas such as Bears Ears that will experience greater visitation and associated impacts in the near future.  more » « less
Award ID(s):
1714972
NSF-PAR ID:
10297926
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Volume:
118
Issue:
21
ISSN:
0027-8424
Page Range / eLocation ID:
e2025047118
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Motivation

    Rapid climate change is altering plant communities around the globe fundamentally. Despite progress in understanding how plants respond to these climate shifts, accumulating evidence suggests that disturbance could not only modify expected plant responses but, in some cases, have larger impacts on compositional shifts than climate change. Climate‐driven disturbances are becoming increasingly common in many biomes and are key drivers of vegetation dynamics at both species and community levels. Palaeoecological records provide valuable observational windows for elucidating the long‐term impacts of these disturbances on plant dynamics; however, sparse resolution and difficulty in disentangling drivers of change limit our ability to understand the impact of disturbance on plant communities. In this targeted review, we highlight emerging opportunities in palaeoecology to advance our understanding about how disturbance, especially fire, impacts the ecological and evolutionary dynamics of terrestrial plant communities.

    Location

    Global examples, with many from North America.

    Conclusions

    We propose a set of palaeoecological and integrative approaches that could greatly enhance our understanding of how disturbance regimes influence global plant dynamics. Specifically, we identify four future study areas: (1) focus on palaeoecological disturbance proxies beyond fire and leverage multi proxy research to examine the influence of interacting disturbances on plant community dynamics; (2) use advances in disturbance and vegetation reconstructions, including ancient sedimentary DNA, to provide the spatial, temporal and taxonomic resolution needed to resolve the relationship between changing disturbance regimes and corresponding shifts in plant community composition; (3) integrate palaeoecological, archaeological and Indigenous knowledge to disentangle the complex interplay between climate, human land use, fire and vegetation structure; and (4) apply “functional palaeoecology” and the synergy between palaeoecology and genetics to understand how fire disturbance has served as a long‐standing selective agent on plants. These frameworks could increase the resolution of disturbance‐driven plant dynamics, potentially providing valuable information for future management.

     
    more » « less
  2. Abstract

    Human activities are altering ecological communities around the globe. Understanding the implications of these changes requires that we consider the composition of those communities. However, composition can be summarized by many metrics which in turn are influenced by different ecological processes. For example, incidence‐based metrics strongly reflect species gains or losses, while abundance‐based metrics are minimally affected by changes in the abundance of small or uncommon species. Furthermore, metrics might be correlated with different predictors. We used a globally distributed experiment to examine variation in species composition within 60 grasslands on six continents. Each site had an identical experimental and sampling design: 24 plots × 4 years. We expressed compositional variation within each site—not across sites—using abundance‐ and incidence‐based metrics of the magnitude of dissimilarity (Bray–Curtis and Sorensen, respectively), abundance‐ and incidence‐based measures of the relative importance of replacement (balanced variation and species turnover, respectively), and species richness at two scales (per plot‐year [alpha] and per site [gamma]). Average compositional variation among all plot‐years at a site was high and similar to spatial variation among plots in the pretreatment year, but lower among years in untreated plots. For both types of metrics, most variation was due to replacement rather than nestedness. Differences among sites in overall within‐site compositional variation were related to several predictors. Environmental heterogeneity (expressed as the CV of total aboveground plant biomass in unfertilized plots of the site) was an important predictor for most metrics. Biomass production was a predictor of species turnover and of alpha diversity but not of other metrics. Continentality (measured as annual temperature range) was a strong predictor of Sorensen dissimilarity. Metrics of compositional variation are moderately correlated: knowing the magnitude of dissimilarity at a site provides little insight into whether the variation is driven by replacement processes. Overall, our understanding of compositional variation at a site is enhanced by considering multiple metrics simultaneously. Monitoring programs that explicitly incorporate these implications, both when designing sampling strategies and analyzing data, will have a stronger ability to understand the compositional variation of systems and to quantify the impacts of human activities.

     
    more » « less
  3. Abstract Questions

    A recently introduced non‐native annual grass,Ventenata dubia, is challenging previous conceptions of community resistance in forest mosaic communities in the Inland Northwest. However, little is known of the drivers and potential ecological impacts of this rapidly expanding species. Here we (1) identify abiotic and biotic habitat characteristics associated with theV. dubiainvasion and examine how these differ betweenV. dubiaand other problematic non‐native annual grasses,Bromus tectorumandTaeniatherum caput‐medusae; and (2) determine how burning influences relationships betweenV. dubiaand plant community composition and structure to address potential impacts on Inland Northwest forest mosaic communities.

    Location

    Blue Mountains of the Inland Northwest, USA.

    Methods

    We measured environmental and plant community characteristics in 110 recently burned and nearby unburned plots. Plots were stratified to capture a range ofV. dubiacover, elevations, biophysical classes, and fire severities. We investigated relationships betweenV. dubia, wildfire, environmental, and plant community characteristics using non‐metric multidimensional scaling and linear regressions.

    Results

    Ventenata dubiawas most abundant in sparsely vegetated, basalt‐derived rocky scablands interspersed throughout the forested landscape. Plant communities most heavily invaded byV. dubiawere largely uninvaded by other non‐native annual grasses.Ventenata dubiawas abundant in both unburned and burned areas, but negative relationships betweenV. dubiacover and community diversity were stronger in burned plots, where keystone sagebrush species were largely absent after fire.

    Conclusions

    Ventenata dubiais expanding the overall invasion footprint into previously uninvaded communities. Burning may exacerbate negative relationships betweenV. dubiaand species richness, evenness, and functional diversity, including in communities that historically rarely burned. Understanding the drivers and impacts of theV. dubiainvasion and recognizing how these differ from other annual grass invasions may provide insight into mechanisms of community invasibility, grass‐fire feedbacks, and aid the development of species‐specific management plans.

     
    more » « less
  4. Stable isotope proxies measured in the proteinaceous fraction of archaeological mollusc shell represents an increasingly important archive for reconstructing past ecological and biogeochemical conditions of nearshore environments. A major issue, however, is understanding the impact of diagenetic alteration in sub-fossil shell isotope values. “Bulk” stable isotope values of nitrogen (δ15N), and especially carbon (δ13C) often shift strongly with increasing C/N ratios in degraded shell, resulting in unreliable data. Here, we examine preservation of an entirely new set of shell paleo-proxies, compound-specific isotopes of amino acids (CSI-AA). We examine carbon (δ13CAA) and nitrogen (δ15NAA) patterns and values from the organic fraction of California mussel (Mytilus californianus) shells from the California Channel Islands. Archaeological shell samples ranging in age from ca. 6,100 to 250 cal BP exhibiting a wide range of degradation states were collected from varied depositional environments (e.g., exposed coastal bluff, buried strata, etc.), and were directly compared to modern shells of the same species and region. Our results indicate organic matter C/N ratios as the best bulk diagnostic indicator of the relative degradation state of shell organic fraction, including changes at the molecular level. Modern shell organic C/N ratios ranged from 2.8 to 3.5, while those in archaeological shell were substantially elevated (3.4–9.5), exhibiting strong and significant negative correlations with bulk δ13C values, weight %C, and weight %N, and a significant but weaker correlation with δ15 N values. An additional “cleaning” step using weak NaOH helped to remove possible exogenous contaminants and improved bulk values of some samples. However, relative molar AA abundances revealed that some AAs, especially the two most abundant, Glycine and Alanine, progressively decreased with increasing C/N ratio. The loss of these amino acids permanently alters bulk isotope values regardless of removal of contaminants. Modeling the bulk isotope change expected due to amino acid molar composition showed major and predictable shifts in bulk δ13C values from selected AA loss, and similarly large but far more variable impacts from exogenous contaminants. In contrast to bulk data, key CSI-AA values and patterns remained almost entirely unaltered, even in the most degraded shell samples, closely matching expected biosynthetic isotope patterns in modern mussel shell. AA isotope proxies for “baseline” (δ15N-Phenylalanine and average δ13C-Essential AAs) and planktonic trophic structure (δ15N-Glutamic Acid and δ15N-Phenylalanine) were not statistically altered with degradation in any sample. Overall, we conclude that while bulk isotopes, particularly δ13C, are very likely to be unreliable in archaeological or subfossil shell with C/N ratios higher than ∼4.0, CSI-AA proxies can still be used to reconstruct past climatic and ecological conditions of the nearshore marine environment. 
    more » « less
  5. Abstract

    In 1949–1951, ecologist Robert H. Whittaker sampled plant community composition at 470 sites in the Siskiyou Mountains (Oregon and California; also known as Klamath or Klamath‐Siskiyou Mountains). His primary goal was to develop methods to quantify plant community variation across environmental gradients, following on his seminal work challenging communities as discrete entities. He selected the Siskiyous because of their diverse and endemic‐rich flora, which he attributed to geological complexity and an ancient stable climate. He chose sites to span gradients of topography, elevation, geologic substrate, and distance from the coast. He used the frequencies of indicator species in his data to assign sampling locations to positions on the topographic gradient, nested within the elevational and substrate gradients. He originated in this study the concept of diversity partitioning, in which gamma diversity (species richness of a community) equals alpha diversity (species richness in homogeneous sites) times beta diversity (species turnover among sites along gradients). Diversity partitioning subsequently became highly influential and new developments on it continue. Whittaker published his Siskiyou work covering paleohistory, biogeography, floristics, vegetation, gradient analysis, and diversity partitioning inEcological Monographsin 1960. Discussed in 2 pages of his 60‐page monograph, diversity partitioning accounts for >95% of its current >4300 citations. In 2006, we retrieved Whittaker's Siskiyou data in hard copy from the Cornell University archives and entered them in a database. We used these data for multiple published analyses, including some based on (re)sampling the approximate locations of a subset of his sites. Because of the continued interest in diversity partitioning and in historic data sets, here we present his data, including 359 sampling locations and their descriptors and, for each sample, a list of species with their estimated percent cover (herbs and shrubs) and numbers by diameter at breast height (DBH) category (trees). Site descriptors include the approximate location (road, trail, or stream), elevation, topographic aspect, geologic substrate (serpentine, gabbro, or diorite), and dominant woody vegetation of each location. For 111 sites, including the small number chosen to represent the distance‐to‐coast gradient, we could not locate his data. There are no copyright restrictions and users of these data should cite this data paper in any publications that result from its use. The authors are available for consultations about and collaborations involving the data.

     
    more » « less