Postglacial viability and colonization in North America's ice-free corridor.

During the Last Glacial Maximum, continental ice sheets isolated Beringia (northeast Siberia and northwest North America) from unglaciated North America. By around 15 to 14 thousand calibrated radiocarbon years before present (cal. kyr bp), glacial retreat opened an approximately 1,500-km-long corridor between the ice sheets. It remains unclear when plants and animals colonized this corridor and it became biologically viable for human migration. We obtained radiocarbon dates, pollen, macrofossils and metagenomic DNA from lake sediment cores in a bottleneck portion of the corridor. We find evidence of steppe vegetation, bison and mammoth by approximately 12.6 cal. kyr bp, followed by open forest, with evidence of moose and elk at about 11.5 cal. kyr bp, and boreal forest approximately 10 cal. kyr bp. Our findings reveal that the first Americans, whether Clovis or earlier groups in unglaciated North America before 12.6 cal. kyr bp, are unlikely to have travelled by this route into the Americas. However, later groups may have used this north-south passageway.

Migrating bison engineer the green wave.

Newly emerging plants provide the best forage for herbivores. To exploit this fleeting resource, migrating herbivores align their movements to surf the wave of spring green-up. With new technology to track migrating animals, the Green Wave Hypothesis has steadily gained empirical support across a diversity of migratory taxa. This hypothesis assumes the green wave is controlled by variation in climate, weather, and topography, and its progression dictates the timing, pace, and extent of migrations. However, aggregate grazers that are also capable of engineering grassland ecosystems make some of the world's most impressive migrations, and it is unclear how the green wave determines their movements. Here we show that Yellowstone's bison (Bison bison) do not choreograph their migratory movements to the wave of spring green-up. Instead, bison modify the green wave as they migrate and graze. While most bison surfed during early spring, they eventually slowed and let the green wave pass them by. However, small-scale experiments indicated that feedback from grazing sustained forage quality. Most importantly, a 6-fold decadal shift in bison density revealed that intense grazing caused grasslands to green up faster, more intensely, and for a longer duration. Our finding broadens our understanding of the ways in which animal movements underpin the foraging benefit of migration. The widely accepted Green Wave Hypothesis needs to be revised to include large aggregate grazers that not only move to find forage, but also engineer plant phenology through grazing, thereby shaping their own migratory movements.

Fifty thousand years of Arctic vegetation and megafaunal diet.

Although it is generally agreed that the Arctic flora is among the youngest and least diverse on Earth, the processes that shaped it are poorly understood. Here we present 50 thousand years (kyr) of Arctic vegetation history, derived from the first large-scale ancient DNA metabarcoding study of circumpolar plant diversity. For this interval we also explore nematode diversity as a proxy for modelling vegetation cover and soil quality, and diets of herbivorous megafaunal mammals, many of which became extinct around 10 kyr bp (before present). For much of the period investigated, Arctic vegetation consisted of dry steppe-tundra dominated by forbs (non-graminoid herbaceous vascular plants). During the Last Glacial Maximum (25-15 kyr bp), diversity declined markedly, although forbs remained dominant. Much changed after 10 kyr bp, with the appearance of moist tundra dominated by woody plants and graminoids. Our analyses indicate that both graminoids and forbs would have featured in megafaunal diets. As such, our findings question the predominance of a Late Quaternary graminoid-dominated Arctic mammoth steppe.

Fossil and genomic evidence constrains the timing of bison arrival in North America.

The arrival of bison in North America marks one of the most successful large-mammal dispersals from Asia within the last million years, yet the timing and nature of this event remain poorly determined. Here, we used a combined paleontological and paleogenomic approach to provide a robust timeline for the entry and subsequent evolution of bison within North America. We characterized two fossil-rich localities in Canada's Yukon and identified the oldest well-constrained bison fossil in North America, a 130,000-y-old steppe bison, Bison cf. priscus We extracted and sequenced mitochondrial genomes from both this bison and from the remains of a recently discovered, ∼120,000-y-old giant long-horned bison, Bison latifrons, from Snowmass, Colorado. We analyzed these and 44 other bison mitogenomes with ages that span the Late Pleistocene, and identified two waves of bison dispersal into North America from Asia, the earliest of which occurred ∼195-135 thousand y ago and preceded the morphological diversification of North American bison, and the second of which occurred during the Late Pleistocene, ∼45-21 thousand y ago. This chronological arc establishes that bison first entered North America during the sea level lowstand accompanying marine isotope stage 6, rejecting earlier records of bison in North America. After their invasion, bison rapidly colonized North America during the last interglaciation, spreading from Alaska through continental North America; they have been continuously resident since then.

Effect of season and superstimulatory treatment on in vivo and in vitro embryo production in wood bison (Bison bison athabascae).

Two experiments were done using a two-by-two design to determine the effects of season and superstimulatory protocol on embryo production in wood bison. In Experiment 1 (in vivo-derived embryos), ovarian superstimulation was induced in female bison during the ovulatory and anovulatory seasons with either two or three doses of FSH given every-other-day (FSH × 2 vs. FSH × 3, respectively). Bison were given hCG to induce ovulation, inseminated 12 and 24 hr after hCG, and embryos were collected 8 days after hCG (n = 10 bison/group). In Experiment 2 (in vitro embryo production), ovarian superstimulation was induced in female bison during the ovulatory and anovulatory seasons with two doses of FSH, and in vivo maturation of the cumulus-oocyte complexes (COC) was induced with hCG at either 48 or 72 hr after the last dose of FSH. COC were collected 34 hr after hCG, and expanded COC were used for in vitro fertilization and culture. In Experiment 1, the number of follicles ≥9 mm, the proportion of follicles that ovulated, the number of CL, and the total number of ova/embryos collected did not differ between seasons or treatment groups, but the number of transferable embryos was greater (p < .05) in the ovulatory season. In Experiment 2, no differences were detected between seasons or treatment groups for any end point. The number of transferable embryos produced per bison was greatest (p < .05) using in vitro fertilization and was unaffected by season (1.5 ± 0.2 and 1.1 ± 0.3 during anovulatory and ovulatory seasons, respectively), in contrast to in vivo embryo production which was affected by season (0.1 ± 0.01 and 0.7 ± 0.2 during anovulatory and ovulatory seasons, respectively). Results demonstrate that transferable embryos can be produced throughout the year in wood bison by both in vivo and in vitro techniques, but the efficiency of embryo production of in vivo-derived embryos is significantly lower during the anovulatory season.

Obtaining Wisent early blastocyst in vitro is a basic for protection and creation of biodiversity for this threatened species.

Wisent, or European bison (Bison bonasus), is listed as "vulnerable" on the IUCN Red List of Threatened Species and is therefore protected by international law. For the first time, a Wisent embryo has been obtained in vitro. This procedure creates a new opportunity to protect and increase Wisent reproductive potential and thereby opens new possibilities for the establishment of a controlled and broad reserve of the gene pool.

Collective decision-making promotes fitness loss in a fusion-fission society.

While collective decision-making is recognised as a significant contributor to fitness in social species, the opposite outcome is also logically possible. We show that collective movement decisions guided by individual bison sharing faulty information about habitat quality promoted the use of ecological traps. The frequent, but short-lived, associations of bison with different spatial knowledge led to a population-wide shift from avoidance to selection of agricultural patches over 9 years in and around Prince Albert National Park, Canada. Bison were more likely to travel to an agricultural patch for the first time by following conspecifics already familiar with agricultural patches. Annual adult mortality increased by 12% due to hunting of bison on agricultural lands. Maladaptive social behaviour accordingly was a major force that contributed to a ~50% population decline in less than a decade. In human-altered landscapes, social learning by group-living species can lead to fitness losses, particularly in fusion-fission societies.

Diet segregation in American bison (Bison bison) of Yellowstone National Park (Wyoming, USA).

BACKGROUND: Body size is a major factor in the nutritional ecology of ruminant mammals. Females, due to their smaller size and smaller rumen, have more rapid food-passage times than males and thereby require higher quality forage. Males are more efficient at converting high-fiber forage into usable energy and thus, are more concerned with quantity. American bison are sexually dimorphic and sexually segregate for the majority of their adult lives, and in Yellowstone National Park, they occur in two distinct subpopulations within the Northern and Central ranges. We used fecal nitrogen and stable isotopes of carbon and nitrogen from American bison to investigate sex-specific differences in diet composition, diet quality, and dietary breadth between the mating season and a time period spanning multiple years, and compared diet indicators for these different time periods between the Northern and Central ranges. RESULTS: During mating season, diet composition of male and female American bison differed significantly; females had higher quality diets, and males had greater dietary breadth. Over the multi-year period, females had higher quality diets and males, greater dietary breadth. Diet segregation for bison in the Central Range was more pronounced during the mating season than for the multi-year period and females had higher quality diets than males. Finally, diet segregation in the Northern Range was more pronounced during the multi-year period than during the mating season, and males had greater dietary breadth. CONCLUSIONS: Female bison in Yellowstone National Park have higher quality diets than males, whereas males ingest a greater diversity of plants or plants parts, and bison from different ranges exhibited more pronounced diet segregation during different times. Collectively, our results suggest that diet segregation in bison of Yellowstone National Park is associated with sex-specific differences in nutritional demands. Altogether, our results highlight the importance of accounting for spatial and temporal heterogeneity when conducting dietary studies on wild ungulates.

A PRELIMINARY ASSESSMENT OF THE EFFECTS OF REPRODUCTIVE SEASONALITY ON SEMINAL PLASMA CHARACTERISTICS IN NORTH AMERICAN BISON (BISON BISON).

The two subspecies of the North American bison, plains ( Bison bison bison) and wood ( Bison bison athabascae) bison, are seasonal breeders. The objective of this study was to conduct a preliminary investigation into the effects of season on semen. To test the hypothesis that there are seasonal effects on seminal plasma, protein profiles of seminal plasma from plains and wood bison (n = 2 of each subspecies) were compared between breeding and nonbreeding seasons. Using two-dimensional gel electrophoresis of seminal plasma proteins, 54 of 170 spots (plains bison) and 19 of 153 spots (wood bison) had differential expression (≥2-fold; P < 0.01) between seasons. Based on immunoblotting, BSP5 and TIMP-2 (two fertility-associated proteins in cattle) were higher during the breeding vs. nonbreeding season. Furthermore, epididymal sperm incubated with seminal plasma from the nonbreeding season had lower postthaw progressive motility (17.33 ± 7.47 vs. 22.09 ± 6.67%; mean ± SD) and an increased ability to undergo a lysophosphatidylcholine-induced acrosome reaction (77.83 ± 8.47 vs. 52.67 ± 7.76%; mean ± SD) as compared to epididymal sperm incubated with seminal plasma from the breeding season. In a heterologous in vitro fertilization system (using bovine oocytes), cleavage rate was higher for sperm exposed to seminal plasma from the breeding vs. the nonbreeding season (75.35 ± 16.55 vs. 33.63 ± 12.44%; mean ± SD). This study suggested that differential expression of seminal plasma characteristics modulating sperm function is one of the mechanisms by which reproductive seasonality affects sperm function in the North American bison.

To follow or not? How animals in fusion-fission societies handle conflicting information during group decision-making.

When group members possess differing information about the environment, they may disagree on the best movement decision. Such conflicts result in group break-ups, and are therefore a fundamental driver of fusion-fission group dynamics. Yet, a paucity of empirical work hampers our understanding of how adaptive evolution has shaped plasticity in collective behaviours that promote and maintain fusion-fission dynamics. Using movement data from GPS-collared bison, we found that individuals constantly associated with other animals possessing different spatial knowledge, and both personal and conspecific information influenced an individual's patch choice decisions. During conflict situations, bison used group familiarity coupled with their knowledge of local foraging options and recently sampled resource quality when deciding to follow or leave a group - a tactic that led to energy-rewarding movements. Natural selection has shaped collective behaviours for coping with social conflicts and resource heterogeneity, which maintain fusion-fission dynamics and play an essential role in animal distribution.

Bison distribution under conflicting foraging strategies: site fidelity vs. energy maximization.

Foraging strategies based on site fidelity and maximization of energy intake rate are two adaptive forces shaping animal behavior. Whereas these strategies can both be evolutionarily stable, they predict conflicting optimal behaviors when population abundance is in decline. In such a case, foragers employing an energy-maximizing strategy should reduce their use of low-quality patches as interference competition becomes less intense for high-quality patches. Foragers using a site fidelity strategy, however, should continue to use familiar patches. Because natural fluctuations in population abundance provide the only non-manipulative opportunity to evaluate adaptation to these evolutionary forces, few studies have examined these foraging strategies simultaneously. Using abundance and space use data from a free-ranging bison (Bison bison) population living in a meadow-forest matrix in Prince Albert National Park, Canada, we determined how individuals balance the trade-off between site fidelity and energy-maximizing patch choice strategies with respect to changes in population abundance. From 1996 to 2005, bison abundance increased from 225 to 475 and then decreased to 225 by 2013. During the period of population increase, population range size increased. This expansion involved the addition of relatively less profitable areas and patches, leading to a decrease in the mean expected profitability of the range. Yet, during the period of population decline, we detected neither a subsequent retraction in population range size nor an increase in mean expected profitability of the range. Further, patch selection models. during the population decline indicated that, as density decreased, bison portrayed stronger fidelity to previously visited meadows, but no increase in selection strength for profitable meadows. Our analysis reveals that an energy-maximizing patch choice strategy alone cannot explain the distribution ofindividuals and populations, and site fidelity is an important evolutionary force shaping animal distribution. Animals may not always forage in the richest patches available, as ecological theory would often predict, but their use of profitable patches is dependent on population dynamics and the strength of site fidelity. Our findings are likewise relevant to applied inquiries such as forecasting species range shifts and reducing human-wildlife conflicts.

Disentangling the effects of climate, density dependence, and harvest on an iconic large herbivore's population dynamics.

Understanding the relative effects of climate, harvest, and density dependence on population dynamics is critical for guiding sound population management, especially for ungulates in arid and semiarid environments experiencing climate change. To address these issues for bison in southern Utah, USA, we applied a Bayesian state-space model to a 72-yr time series of abundance counts. While accounting for known harvest (as well as live removal) from the population, we found that the bison population in southern Utah exhibited a strong potential to grow from low density (ß0 = 0.26; Bayesian credible interval based on 95% of the highest posterior density [BCI] = 0.19-0.33), and weak but statistically significant density dependence (ß1 = -0.02, BCI = -0.04 to -0.004). Early spring temperatures also had strong positive effects on population growth (Pfat1 = 0.09, BCI = 0.04-0.14), much more so than precipitation and other temperature-related variables (model weight > three times more than that for other climate variables). Although we hypothesized that harvest is the primary driving force of bison population dynamics in southern Utah, our elasticity analysis indicated that changes in early spring temperature could have a greater relative effect on equilibrium abundance than either harvest or. the strength of density dependence. Our findings highlight the utility of incorporating elasticity analyses into state-space population models, and the need to include climatic processes in wildlife management policies and planning.

Effect of reproductive seasonality on gamete quality in the North American bison (Bison bison bison).

The objective was to investigate the effects of reproductive seasonality on gamete quality in plains bison (Bison bison bison). Epididymal sperm (n = 61 per season), collected during the breeding season (July-September), had significantly higher post-thaw total motility (36.76 ± 14.18 vs 31.24 ± 12.74%), and lower linearity (0.36 ± 0.06 vs 0.39 ± 0.04) and wobbliness (0.49 ± 0.04 vs 0.51 ± 0.03; mean ± SD) compared to non-breeding season (January-March) samples. Representative samples (n = 4) from each season were used in heterologous IVF trials using cattle oocytes. Cleavage, morulae and blastocyst percentage were higher for breeding vs non-breeding season sperm samples (81.88 ± 6.8 vs 49.94 ± 6.77; 41.89 ± 13.40 vs 27.08 ± 23.21; and 30.49 ± 17.87 vs 13.72 ± 18.98%, respectively). Plains bison ovaries collected during the breeding (n = 97 pairs) and non-breeding (n = 100 pairs) seasons were classified as luteal or follicular. Oocytes recovered from these ovaries were classified into five grades based on morphology. There was no significant difference in the number of luteal ovaries or grades of oocytes recovered. Oocytes were matured, fertilized (with frozen sperm from three bison bulls) and cultured in vitro. Cleavage percentage was higher for oocytes collected during breeding vs non-breeding season (83.72 ± 6.42 vs 73.98 ± 6.43), with no significant difference in subsequent development to blastocysts. In summary, epididymal sperm from non-breeding season had decreased total motility and resulted in reduced embryo production in vitro. Oocytes collected during non-breeding season had reduced ability to be matured, fertilized and/or undergo cleavage in vitro. Data suggested that season influenced gamete quality in plains bison.

A memory-based foraging tactic reveals an adaptive mechanism for restricted space use.

The restricted area of space used by most mobile animals is thought to result from fitness-rewarding decisions derived from gaining information about the environment. Yet, assessments of how animals deal with uncertainty using memory have been largely theoretical, and an empirically derived mechanism explaining restricted space use in animals is still lacking. Using a patch-to-patch movement analysis, we investigated predictions of how free-ranging bison (Bison bison) living in a meadow-forest matrix use memory to reduce uncertainty in energy intake rate. Results indicate that bison remembered pertinent information about location and quality of meadows, and they used this information to selectively move to meadows of higher profitability. Moreover, bison chose profitable meadows they had previously visited, and this choice was stronger after visiting a relatively poor quality meadow. Our work demonstrates a link between memory, energy gains and restricted space use while establishing a fitness-based integration of movement, cognitive and spatial ecology.

Integrating multiple lines of evidence into historical biogeography hypothesis testing: a Bison bison case study.

One of the grand goals of historical biogeography is to understand how and why species' population sizes and distributions change over time. Multiple types of data drawn from disparate fields, combined into a single modelling framework, are necessary to document changes in a species's demography and distribution, and to determine the drivers responsible for change. Yet truly integrated approaches are challenging and rarely performed. Here, we discuss a modelling framework that integrates spatio-temporal fossil data, ancient DNA, palaeoclimatological reconstructions, bioclimatic envelope modelling and coalescence models in order to statistically test alternative hypotheses of demographic and potential distributional changes for the iconic American bison (Bison bison). Using different assumptions about the evolution of the bioclimatic niche, we generate hypothetical distributional and demographic histories of the species. We then test these demographic models by comparing the genetic signature predicted by serial coalescence against sequence data derived from subfossils and modern populations. Our results supported demographic models that include both climate and human-associated drivers of population declines. This synthetic approach, integrating palaeoclimatology, bioclimatic envelopes, serial coalescence, spatio-temporal fossil data and heterochronous DNA sequences, improves understanding of species' historical biogeography by allowing consideration of both abiotic and biotic interactions at the population level.

Integrating population- and individual-level information in a movement model of Yellowstone bison.

Throughout the world, fragmentation of landscapes by human activities has constrained the opportunity for large herbivores to migrate. Conflict between people and wildlife results when migrating animals transmit disease to livestock, damage property, and threaten human safety. Mitigating this conflict requires understanding the forces that shape migration patterns. Bison Bos bison migrating from Yellowstone National Park into the state of Montana during winter and spring concern ranchers on lands surrounding the park because bison can transmit brucellosis (Brucella abortus) to cattle. Migrations have been constrained, with bison being lethally removed or moved back into the park. We developed a state-space model to support decisions on bison management aimed at mitigating conflict with landowners outside the park. The model integrated recent GPS observations with 22 years (1990-2012) of aerial counts to forecast monthly distributions and identify factors driving migration. Wintering areas were located along decreasing elevation gradients, and bison accumulated in wintering areas prior to moving to areas progressively lower in elevation. Bison movements were affected by time since the onset of snowpack, snowpack magnitude, standing crop, and herd size. Migration pathways were increasingly used over time, suggesting that experience or learning influenced movements. To support adaptive management of Yellowstone bison, we forecast future movements to evaluate alternatives. Our approach of developing models capable of making explicit probabilistic forecasts of large herbivore movements and seasonal distributions is applicable to managing the migratory movements of large herbivores worldwide. These forecasts allow managers to develop and refine strategies in advance, and promote sound decision-making that reduces conflict as migratory animals come into contact with people.

Ovarian synchronisation in wood bison (Bison bison athabascae) during the anovulatory season.

Two experiments were performed in wood bison during the anovulatory season to establish an effective protocol for ovarian synchronisation. In an untreated control phase, bison cows (n=19) were examined daily to establish the interval to new follicular wave emergence (4.9±0.7 days) for the purposes of comparison with the experimental treatments. In Experiment 1, bison were treated by transvaginal ultrasound-guided follicular ablation (n=9) or with 2mg, i.m., 17ß-oestradiol (n=10). In Experiment 2, bison were treated by follicular ablation (n=9) or with 2mg, i.m., 17ß-oestradiol +100mg, i.m., progesterone (n=10). In Experiment 1, the interval to new wave emergence for control, follicular ablation and 17ß-oestradiol-treated groups was 4.9±0.7, 1.1±0.1 and 3.1±0.4 days, respectively (P<0.05). The degree of synchrony was 2.4±0.4, 0.2±0.1 and 0.8±0.2 days, respectively (P<0.05). In Experiment 2, the interval to new wave emergence for control, follicular ablation and 17ß-oestradiol + progesterone-treated groups was 4.9±0.7, 1.2±0.2 and 3.3±0.3 days, respectively (P<0.05), and the degree of synchrony was 2.4±0.4, 0.2±0.1, and 0.8±0.2 days, respectively (P<0.05). The degree of synchrony did not differ between ablation and hormone treatment groups in either experiment, but was greater in treatment groups than in the untreated control phase. Both follicular ablation and hormone treatment shortened and decreased the variability in the interval to follicular wave emergence in bison, but wave emergence occurred earlier after follicular ablation.

Functional consequences of animal community changes in managed grasslands: An application of the CAFE approach.

In the midst of an ongoing biodiversity crisis, much research has focused on species losses and their impacts on ecosystem functioning. The functional consequences (ecosystem response) of shifts in communities are shaped not only by changes in species richness, but also by compositional shifts that result from species losses and gains. Species differ in their contribution to ecosystem functioning, so species identity underlies the consequences of species losses and gains on ecosystem functions. Such research is critical to better predict the impact of disturbances on communities and ecosystems. We used the "Community Assembly and the Functioning of Ecosystems" (CAFE) approach, a modification of the Price equation to understand the functional consequences and relative effects of richness and composition changes in small nonvolant mammal and dung beetle communities as a result of two common disturbances in North American prairie restorations, prescribed fire and the reintroduction of large grazing mammals. Previous research in this system has shown dung beetles are critically important decomposers, while small mammals modulate much energy in prairie food webs. We found that dung beetle communities were more responsive to bison reintroduction and prescribed fires than small nonvolant mammals. Dung beetle richness increased after bison reintroduction, with higher dung beetle community biomass resulting from changes in remaining species (context-dependent component) rather than species turnover (richness components); prescribed fire caused a minor increase in dung beetle biomass for the same reason. For small mammals, bison reintroduction reduced energy transfer through the loss of species, while prescribed fire had little impact on either small mammal richness or energy transfer. The CAFE approach demonstrates how bison reintroduction controls small nonvolant mammal communities by increasing prairie food web complexity, and increases dung beetle populations with possible benefits for soil health through dung mineralization and soil bioturbation. Prescribed fires, however, have little effect on small mammals and dung beetles, suggesting a resilience to fire. These findings illustrate the key role of re-establishing historical disturbance regimes when restoring endangered prairie ecosystems and their ecological function.

Spatial-social familiarity complements the spatial-social interface: evidence from Yellowstone bison.

Social animals make behavioural decisions based on local habitat and conspecifics, as well as memorized past experience (i.e. 'familiarity') with habitat and conspecifics. Here, we develop a conceptual and empirical understanding of how spatial and social familiarity fit within the spatial-social interface-a novel framework integrating the spatial and social components of animal behaviour. We conducted a multi-scale analysis of the movements of GPS-collared plains bison (Bison bison, n = 66) residing in and around Yellowstone National Park, USA. We found that both spatial and social familiarity mediate how individuals respond to their spatial and social environments. For instance, individuals with high spatial familiarity rely on their own knowledge as opposed to their conspecifics, and individuals with high social familiarity rely more strongly on the movement of conspecifics to guide their own movement. We also found that fine-scale spatial and social phenotypes often scale up to broad-scale phenotypes. For instance, bison that select more strongly to align with their nearest neighbour have larger home ranges. By integrating spatial and social familiarity into the spatial-social interface, we demonstrate the utility of the interface for testing hypotheses, while also highlighting the pervasive importance of cognitive mechanisms in animal behaviour. This article is part of the theme issue 'The spatial-social interface: a theoretical and empirical integration'.

Theoretical model of impact mitigation mechanisms inherent to the North American bison skull.

North American bison (Bovidae: Bison bison) incur blunt impacts to the interparietal and frontal bones when they engage in head-to-head fights. To investigate the impact mitigation of these bones, a finite element analysis (FEA) of the skull under loading conditions was performed. Based on anatomical and histological studies, the interparietal and frontal bones are both comprised of a combination of haversian and plexiform bone and are both underlain by bony septa. Additionally, the interparietal bone is thicker than the frontal bone. Data regarding the mechanical properties of bison bone are scarce, but the results of a phylogenetic analysis infer that the material properties of the closely related domestic cow bone are a suitable proxy for use in the FEA. Results of the FEA suggest that the thickness of the interparietal bone in conjunction with the bony septa may prevent fracture stresses by helping to absorb and disperse the blunt impact energy throughout the skull. Monotonic stress levels of 294 MPa, which are below the compressive strength of bone were exhibited in the simulated bison head impacts indicating no fracture of the bones.