skip to main content


Title: BITS : A Bayesian Isotope Turnover and Sampling model for strontium isotopes in proboscideans and its potential utility in movement ecology
Abstract

Strontium isotope ratios (87Sr/86Sr) of incrementally grown tissues have been widely used to study movement ecology and migration of animals. However, the time scale of87Sr/86Sr incorporation from the environment into tissue and how it may influence data interpretation are still poorly understood. Using the relocation of a zoo elephant (Loxodonta africana) named Misha, we characterise and model the87Sr/86Sr turnover process using high‐resolution measurements of its tusk dentine. We seek to develop a framework that can improve quantitative interpretation of87Sr/86Sr data in tissues.

The87Sr/86Sr transition associated with the relocation is measured using laser ablation inductively coupled plasma mass spectrometry (LA‐ICP‐MS) on a prepared tusk slab. We develop a turnover model (BITS), with a rapidly exchanging central pool and a slowly exchanging peripheral pool, in a Bayesian statistical framework. The measured dentine data are first used to calibrate model parameters. The parameters are then used to estimate possible87Sr/86Sr input time series from two datasets via model inversion: a fidelity test using Misha's dentine data and a case study using published dentine measurements from an Alaskan Woolly Mammoth (Mammuthus primigenius).

The LA‐ICP‐MS data are consistent with a two‐compartment turnover process with equivalent half‐lives of 41 days for the central pool and 170 days for the peripheral pool. The model inversion shows good fidelity when estimating the intake87Sr/86Sr time series associated with Misha's relocation. In the case study, the model suggests an abrupt pattern of change in, and a much wider range of, intake87Sr/86Sr values than expressed in the woolly mammoth dentine data themselves.

Our framework bridges the gap between environmental87Sr/86Sr variation and data measured in tusk dentine or other incrementally grown tissues. It could be coupled with movement models and additional isotope tracers to study seasonal residency or the spatial and temporal patterns of movement/migration. The generic turnover processes can be adapted to other isotope systems, additional incremental tissues, or other organisms, thus expanding our modelling toolkit to investigate niche partitioning, life history traits and behavioural patterns in conservation biology, archaeology and paleoecology.

 
more » « less
NSF-PAR ID:
10472215
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Methods in Ecology and Evolution
Volume:
14
Issue:
11
ISSN:
2041-210X
Format(s):
Medium: X Size: p. 2800-2813
Size(s):
["p. 2800-2813"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Little is known about woolly mammoth (Mammuthus primigenius) mobility and range. Here we use high temporal resolution sequential analyses of strontium isotope ratios along an entire 1.7-meter-long tusk to reconstruct the movements of an Arctic woolly mammoth that lived 17,100 years ago, during the last ice age. We use an isotope-guided random walk approach to compare the tusk’s strontium and oxygen isotope profiles to isotopic maps. Our modeling reveals patterns of movement across a geographically extensive range during the animal’s ~28-year life span that varied with life stages. Maintenance of this level of mobility by megafaunal species such as mammoth would have been increasingly difficult as the ice age ended and the environment changed at high latitudes.

     
    more » « less
  2. Laser ablation tandem mass spectrometry is a burgeoning field forin situRb‐Sr geochronology. Here, we determined simultaneous isotope ratios of87Sr/86Sr and87Rb/86Sr in metamorphic biotite from western Maine, using an ESL™ imageGEO™193 excimer laser ablation system coupled to a Thermo Scientific™ Neoma™ MC‐ICP‐MS/MS. Measurements were made on Faraday cups with Rb+at mass 87; Sr isotopes were reacted with SF6gas and measured as SrF+at masses 103–107. Twenty‐two laser spots in biotite from a single sample yield a traditional Rb‐Sr isochron date of 289 ± 6 Ma. Time‐resolved signals reveal significant zoning in87Sr/86Sr and87Rb/86Sr within single spot analyses, which were used to construct single spot isochrons. Individual laser spots contain multiple isochronous subpopulations; some spots contain up to three distinct Rb‐Sr isochrons that are decoupled from variations in Rb/Sr. Thirty‐five isochron dates were determined using this sub‐spot approach, with87Sr/86Sr intercepts that systematically vary with Rb‐Sr date; two‐point isochrons were calculated for individual integrations (n= 780) based on these variable intercepts. Both methods yield age peaks at 303, 270 and 240 Ma. These data suggest that the Rb‐Sr system has the potential to record multiple heating, cooling or fluid‐alteration events spanning ~ 100 My within small domains in single biotite crystals.

     
    more » « less
  3. Abstract

    Young mafic lavas from the East African Western Rift record melting of subcontinental lithospheric mantle that was metasomatically modified by multiple tectonic events. We report new isotope data from monogenetic cinder cones near Bufumbira, Uganda, in the Virunga Volcanic Field:87Sr/86Sr = 0.7059–0.7079,εNd = −6.5 to −1.3,εHf = −6.3 to +0.9,208Pb/204Pb = 40.1–40.7,207Pb/204Pb = 15.68–15.75, and206Pb/204Pb = 19.27–19.45. Olivine phenocrysts from the Bufumbira lavas have3He/4He = 6.0–7.4RA. The isotopic data, in conjunction with major and trace element systematics, indicate that primitive Bufumbira magmas are derived from two different metasomatized lithospheric source domains. Melts generated by lower degrees of melting record greater contributions from ∼1 to 2 Ga isotopically enriched garnet‐amphibole‐phlogopite pyroxenite veins within the lithosphere. As melting progresses, these vein melts become increasingly diluted by melts that originate near the lithosphere/asthenosphere boundary, shifting the isotopic compositions toward the common lithospheric mantle (CLM) proposed by Furman and Graham (1999,https://doi.org/10.1016/s0024-4937(99)00031-6). This ∼450–500 Ma source domain appears to underlie all Western Rift volcanic provinces and is characterized by87Sr/86Sr ∼ 0.705,εNd∼ 0,εHf∼ +1 to +3,206Pb/204Pb ∼ 19.0–19.2,208Pb/204Pb ∼ 39.7, and3He/4He ∼ 7RA. Basal portions of the dense subcontinental lithospheric mantle may become gravitationally unstable and founder into underlying warmer asthenosphere, exposing surfaces where melting of locally heterogeneous veins produces small‐volume, alkaline mafic melts. Mafic lavas from all Western Rift volcanic provinces record mixing between the CLM and locally variable metasomatized source domains, suggesting this style of melt generation is fundamental to the development of magma‐poor rifts.

     
    more » « less
  4. Abstract

    Blood Falls is a hypersaline, iron‐rich discharge at the terminus of the Taylor Glacier in the McMurdo Dry Valleys, Antarctica. In November 2014, brine in a conduit within the glacier was penetrated and sampled using clean‐entry techniques and a thermoelectric melting probe called the IceMole. We analyzed the englacial brine sample for filterable iron (fFe), total Fe, major cations and anions, nutrients, organic carbon, and perchlorate. In addition, aliquots were analyzed for minor and trace elements and isotopes including δD and δ18O of water, δ34S and δ18O of sulfate,234U,238U, δ11B,87Sr/86Sr, and δ81Br. These measurements were made in order to (1) determine the source and geochemical evolution of the brine and (2) compare the chemistry of the brine to that of nearby hypersaline lake waters and previous supraglacially sampled collections of Blood Falls outflow that were interpreted asendmemberbrines. The englacial brine had higher Clconcentrations than the Blood Falls end‐member outflow; however, other constituents were similar. The isotope data indicate that the water in the brine is derived from glacier melt. The H4SiO4concentrations and U and Sr isotope suggest a high degree of chemical weathering products. The brine has a low N:P ratio of ~7.2 with most of the dissolved inorganic nitrogen in the form of NH4+. Dissolved organic carbon concentrations are similar to end‐member outflow values. Our results provide strong evidence that the original source of solutes in the brine was ancient seawater, which has been modified with the addition of chemical weathering products.

     
    more » « less
  5. The Gejiu alkaline complex (GAC) within the western part of the Youjiang Basin provides a window to investigate the evolution in the junction of Cathaysia, Yangtze and Indochina blocks. Here, we investigate the GAC in terms of their petrology, zircon U–Pb geochronology, whole‐rock geochemistry, and Sr–Nd isotopic data to gain insights into the origin and evolution of the alkaline magma. The GAC is lithologically composed of alkali syenites and feldspathoid syenites, in which some were altered into sericite syenites. Zircon U–Pb dating of the alkali syenites yielded an age of 85.03 ± 0.47 Ma, which is slightly older than the feldspathoid syenites. The alkali syenites and feldspathoid syenites are silica‐saturated and silica‐undersaturated, respectively, and are characterized by high alkalinity with K2O + Na2O of 11.55–17.08 wt% and Al2O3of 18.57–22.49 wt%, low MgO of 0.11–1.39 wt%, weakly negative Eu anomalies, enrichments of LILEs such as Th and U, HFSEs like Zr and Hf but depletion of Ba, Sr, Nb, Ta, P, Ti, strongly fractionated LREEs to HREEs. Their uniform Sr–Nd isotope composition with initial87Sr/86Sr = 0.708802–0.710571 andεNd(t) = −7.1 to −6.6 indicates that they were products of a homologous magma. They crystallized atc. 810–956°C and have a relatively high magmatic oxygen fugacity. Our geochemical and isotopic data proved that the GAC magma was derived from the low‐degree partial melting (<10%) of a phlogopite‐bearing‐enriched mantle that was metasomatized by subducting sediments and originated possibly from the spinel and garnet transition zone at a depth of 60–80 km and a pressure of about 1.8–2.4 Gpa. The primary magma experienced protracted two‐stage crystal fractionation of clinopyroxene+amphibole and biotite+K‐feldspar, forming alkali syenites and feldspathoid syenites. Crustal contamination plays a negligible role in their formation. Considering previous tectonic studies, it was therefore proposed that the GAC formed in an extensional tectonic setting related to the Neo‐Tethyan tectonic regions during the Late Cretaceous.

     
    more » « less