Geology ReportsSearch

USGS · 70272766

Do bighorn sheep use desert tortoise burrow spoil piles as mineral licks in southern California?

Abstract

In a previous study, bighorn sheep ( Ovis canadensis ) were periodically photographed by trail cameras when they visited desert tortoise ( Gopherus agassizii ) burrows in the San Bernardino Mountains near Palm Springs, California, USA. The authors suggested that bighorn sheep may utilize the excavated spoil material from burrows as a mineral lick. To test that supposition, we collected soil samples from the spoil piles and undisturbed surface soil up to 1 m away at two desert tortoise study sites with bighorn sheep at the aforementioned site and another nearby site. We hypothesized that sodium concentrations would be greater in spoil than in nearby surface soil samples. Common major cations were quantified in soil samples using inductively coupled plasma-mass spectrometry following microwave-assisted digestion to solubilize the fraction of elements that could become biologically or environmentally available. Mean sodium concentrations were greater in burrow spoil samples than in undisturbed soil samples at both sites, likely due to aqueous leaching of sodium from surface to subsurface soil horizons. Although we did not collect evidence of geophagy by bighorn sheep, the behavior is well known. However, trail camera photographs of bighorn sheep with their head next to the entrance of tortoise burrows suggests that they may periodically utilize sodium and other elements brought to the surface by desert tortoises that are otherwise unavailable to bighorn sheep on the surface.

Explore related subjects

Keep this discovery

BibTeXRIS

Jeffrey E. Lovich, Danielle M. Cleveland, Charles B. Yackulic, Kristy L. Cummings, Michele R. Puffer, Molly J. Bechtel, Christopher R. Tracy, William J. Hoese, Gerardo Avila, Megan Peukert, Samantha Hubbard. 2025. Do bighorn sheep use desert tortoise burrow spoil piles as mineral licks in southern California?. https://doi.org/10.3160/0038-3872-124.2.55

Cite the original work for its findings. Save a collection to share your selection of sources.

Discover connections

Connections use source metadata and explicit phrase matches, not verified experimental comparisons.

KEEP EXPLORING

Related discoveries

Widespread terrestrial ecosystem disruption at the onset of the Paleocene–Eocene Thermal Maximum

The Paleocene–Eocene Thermal Maximum (PETM, ~56 Mya) interval was marked by massive 13 C-depleted carbon emissions into the ocean/atmosphere system, manifested as a negative carbon isotope excursion (CIE) in sedimentary components, and ~5 °C global average warming. Episodes of hydrological perturbations and soil-erosion have been widely documented for the PETM but their link with vegetation- and carbon cycle changes remain poorly constrained. Here, we present organic microfossil evidence showing a strong increase in fern-dominated pioneer vegetation that replaced coniferous forests on the margin of the Norwegian Sea during the first millennia of the CIE. With the present stratigraphic constraints, the “fern spike” occurred simultaneously in terrestrial settings along the North Sea, Arctic Ocean, the US east coast and in southern Australia, indicating that pioneer vegetation persisted for several millennia following a partial collapse of previously stable terrestrial ecosystems. Both the ferns and influx of microcharcoal imply recurrent physical disturbance, including soil destabilization and erosion, potentially linked to droughts, wildfires, and strong hydrological forcing resulting from extreme climate change. Together with evidence for reworked clay minerals and ancient organic matter (kerogen), these findings show that highly disturbed terrestrial ecosystems were widespread across mid- and high-latitude regions globally. Carbon cycle model simulations suggest that a substantial loss of standing and buried biomass, along with oxidation of soil organic matter, acted as important positive feedbacks during the onset of the CIE. Additionally, enhanced kerogen weathering likely contributed as another major positive feedback throughout both the onset and main phase of the CIE.

Proceedings of the National Academy of Sciences

Efficient physics‐informed ground‐motion simulations with reduced‐order models: CyberShake implications and high‐resolution site terms for southern San Andreas fault earthquakes

Recent advances in Probabilistic Seismic Hazard Analysis (PSHA) leverage physics‐based ground‐motion simulations to estimate seismic hazard, such as the CyberShake project. However, computational costs quickly escalate when performing PSHA for numerous faults or sites and can become prohibitively expensive. To reduce computational demands, CyberShake uses reciprocity and interpolates physics‐informed corrections from simulations conducted at fewer locations, but the accuracy of these interpolations remains poorly quantified. To quantify the interpolation accuracy, we derive high‐resolution, frequency‐dependent site terms for southern California and compare them with interpolated site terms using the CyberShake approach. We accomplish this by performing a set of earthquake point‐source simulations distributed along the nonplanar fault geometry for the southern San Andreas fault (SSAF) extending from Bombay Beach to Lake Hughes. Using SeisSol, we simulate three minutes of viscoelastic seismic wave propagation for these sources and store the horizontal‐component Green’s functions for 480,000 sites. We then use a scientific machine learning approach based on interpolated proper orthogonal decomposition to construct an accurate reduced‐order model of the Green’s functions to efficiently predict effective amplitude spectra (EAS) for finite‐source rupture models of SSAF earthquakes. Using minimum curvature interpolation with tension, as used in CyberShake, we compare the interpolated site terms against our high‐resolution site terms. We identify local discrepancies with EAS differing by up to a factor of approximately three. Furthermore, we identify locations where unexpectedly high or low ground motions are missed when using the interpolated dataset for these earthquakes. We estimate that our approach may be used within CyberShake to reduce the time‐to‐solution by a factor of 336 for the entire earthquake rupture forecast. Our analysis of physics‐based site terms provides more insight into the seismic hazard due to SSAF ruptures and guides future developments by combining high‐performance computing and reduced‐order modeling techniques for PSHA.

California

Orientation dependence of probabilistic seismic hazard estimates from CyberShake physics-based simulations

Earthquake ground‐motion intensities, such as pseudospectral accelerations (SAs), can vary significantly with horizontal orientation. However, conducting probabilistic seismic hazard analysis (PSHA) for each horizontal orientation is challenging because current ground‐motion models used in PSHA consider only a single horizontal intensity value, usually the median across all orientations, known as RotD50. To address this limitation, we employ physics‐based simulations for PSHA, which contain full waveforms from which ground‐motion intensities can be computed for all horizontal orientations to study directional seismic hazard. We apply our approach to the latest CyberShake study of the Greater Los Angeles metropolitan area, developed by the Statewide California Earthquake Center, finding that seismic hazard at a 2475‐yr return period, a common value used for earthquake‐resistant design, varies significantly with horizontal orientation. For instance, for SAs at 3 s, the maximum seismic hazard across all horizontal orientations is, on average, 15% higher than the median RotD50 hazard, with these differences becoming more pronounced at longer periods. These observed variations can generally be attributed to physical mechanisms that polarize seismic waves, such as the radiation pattern of the earthquake source and the influence of the subsurface structure. These results may have important implications for earthquake engineering applications, particularly for long‐period structures in areas with substantial horizontal variations in seismic hazard.

California