Isotope Landscapes for Terrestrial Migration Research

Subarctic climate for the earliest Homo sapiens in Europe

The expansion of Homo sapiens across Eurasia marked a major milestone in human evolution that would eventually lead to our species being found across every continent. Current models propose that these expansions occurred only during episodes of warm climate, based on age correlations between archaeological and climatic records. Here, we obtain direct evidence for the temperatures faced by some of these humans through the oxygen isotope analysis of faunal remains from Bacho Kiro Cave, Bulgaria, the earliest clear record of H. sapiens in Europe. The results indicate that humans ∼45,000 years ago experienced subarctic climates with far colder climatic conditions than previously suggested. This demonstrates that the early presence of H. sapiens in Europe was not contingent on warm climates. Our results necessitate the revision of key models of human expansion and highlight the need for a less deterministic role of climate in the study of our evolutionary history.

Geographic variation in dispersal of western burrowing owl (Athene cunicularia hypugaea) populations across North America

Dispersal is one of the key elements of species’ metapopulation dynamics and, hence, influences global conservation status. Furthermore, determining the geographic variation in magnitude and direction of dispersal throughout a species’ distribution may expand our understanding of the causes of population declines in species of conservation concern. For instance, western burrowing owl (Athene cunicularia hypugaea) populations have declined at the northern and eastern edge of their breeding distribution during the 20th century. In the same period, large areas of thornscrub that did not support breeding owls were converted to irrigated agriculture in the southern edge of the subspecies’ breeding distribution in northwestern Mexico. These farmlands now support some of the highest breeding densities of owls. We tested the hypothesis that owls that colonized this recently created habitat originated from declining migratory populations from the northern portion of the subspecies’ range. We used stable isotopes 2H, 13C, and 15N in owl feathers to infer breeding dispersal patterns throughout the subspecies’ breeding range. Populations near the northern edge of the subspecies’ breeding range had immigrants that dispersed over larger distances than immigrants at low and mid latitude populations. However, agricultural populations in northwestern Mexico disrupted this latitudinal pattern, attracting owls from more distant locations. We also found immigrants originated from further distances in declining populations than increasing populations. Stable isotopes provided no evidence of contemporaneous breeding dispersal from Canadian populations to northwestern Mexico but suggest that agricultural areas in the southern edge of the subspecies’ distribution have altered the continental dispersal pattern.

Using oxygen and hydrogen stable isotopes to track the migratory movement of Sharp-shinned Hawks (Accipiter striatus) along Western Flyways of North America

The large-scale patterns of movement for the Sharp-shinned Hawk (Accipiter striatus), a small forest hawk found throughout western North America, are largely unknown. However, based on field observations we set out to test the hypothesis that juvenile migratory A. striatus caught along two distinct migration routes on opposite sides of the Sierra Nevada Mountains of North America (Pacific Coast and Intermountain Migratory Flyways) come from geographically different natal populations. We applied stable isotope analysis of hydrogen (H) and oxygen (O) of feathers, and large scale models of spatial isotopic variation (isoscapes) to formulate spatially explicit predictions of the origin of the migrant birds. Novel relationships were assessed between the measured hydrogen and oxygen isotope values of feathers from A. striatus museum specimens of known origin and the isoscape modeled hydrogen and oxygen isotope values of precipitation at those known locations. We used these relationships to predict the origin regions for birds migrating along the two flyways from the measured isotope values of migrant’s feathers and the associated hydrogen and oxygen isotopic composition of precipitation where these feathers were formed. The birds from the two migration routes had overlap in their natal/breeding origins and did not differentiate into fully separate migratory populations, with birds from the Pacific Coast Migratory Flyway showing broader natal geographic origins than those from the Intermountain Flyway. The methodology based on oxygen isotopes had, in general, less predictive power than the one based on hydrogen. There was broad agreement between the two isotope approaches in the geographic assignment of the origins of birds migrating along the Pacific Coast Flyway, but not for those migrating along the Intermountain Migratory Flyway. These results are discussed in terms of their implications for conservation efforts of A. striatus in western North America, and the use of combined hydrogen and oxygen stable isotope analysis to track the movement of birds of prey on continental scales.

Temporal and spatial trends in marine carbon isotopes in the Arctic Ocean and implications for food web studies

The Arctic is undergoing unprecedented environmental change. Rapid warming, decline in sea ice extent, increase in riverine input, ocean acidification and changes in primary productivity are creating a crucible for multiple concurrent environmental stressors, with unknown consequences for the entire arctic ecosystem. Here, we synthesized 30 years of data on the stable carbon isotope (δ¹³ C) signatures in dissolved inorganic carbon (δ¹³ C-DIC; 1977-2014), marine and riverine particulate organic carbon (δ¹³ C-POC; 1986-2013) and tissues of marine mammals in the Arctic. δ¹³ C values in consumers can change as a result of environmentally driven variation in the δ¹³ C values at the base of the food web or alteration in the trophic structure, thus providing a method to assess the sensitivity of food webs to environmental change. Our synthesis reveals a spatially heterogeneous and temporally evolving δ¹³ C baseline, with spatial gradients in the δ¹³ C-POC values between arctic shelves and arctic basins likely driven by differences in productivity and riverine and coastal influence. We report a decline in δ¹³ C-DIC values (-0.011‰ per year) in the Arctic, reflecting increasing anthropogenic carbon dioxide (CO₂ ) in the Arctic Ocean (i.e. Suess effect), which is larger than predicted. The larger decline in δ¹³ C-POC values and δ¹³ C in arctic marine mammals reflects the anthropogenic CO₂ signal as well as the influence of a changing arctic environment. Combining the influence of changing sea ice conditions and isotopic fractionation by phytoplankton, we explain the decadal decline in δ¹³ C-POC values in the Arctic Ocean and partially explain the δ¹³ C values in marine mammals with consideration of time-varying integration of δ¹³ C values. The response of the arctic ecosystem to ongoing environmental change is stronger than we would predict theoretically, which has tremendous implications for the study of food webs in the rapidly changing Arctic Ocean.

The Importance of Isotopic Turnover for Understanding Key Aspects of Animal Ecology and Nutrition

Stable isotope-based methods have proved to be immensely valuable for ecological studies ranging in focus from animal movements to species interactions and community structure. Nevertheless, the use of these methods is dependent on assumptions about the incorporation and turnover of isotopes within animal tissues, which are oftentimes not explicitly acknowledged and vetted. Thus, the purpose of this review is to provide an overview of the estimation of stable isotope turnover rates in animals, and to highlight the importance of these estimates for ecological studies in terrestrial, freshwater, and marine systems that may use a wide range of stable isotopes. Specifically, we discuss 1) the factors that contribute to variation in turnover among individuals and across species, which influences the use of stable isotopes for diet reconstructions, 2) the differences in turnover among tissues that underlie so-called ‘isotopic clocks’, which are used to estimate the timing of dietary shifts, and 3) the use of turnover rates to estimate nutritional requirements and reconstruct histories of nutritional stress from tissue isotope signatures. As we discuss these topics, we highlight recent works that have effectively used estimates of turnover to design and execute informative ecological studies. Our concluding remarks suggest several steps that will improve our understanding of isotopic turnover and support its integration into a wider range of ecological studies.

Refining Stable Oxygen and Hydrogen Isoscapes for the Identification of Human Remains in Mississippi

Isoscape refinement is an essential component for accurately predicting region-of-origin in forensic investigations involving isotope analysis of unidentified human remains. Stable oxygen (δ¹⁸ O) and hydrogen (δ² H) isotopes were measured from 57 tap water samples collected across Mississippi to model refined isoscapes for the state. A tap water conversion equation, δ¹⁸ O_tw =1.64 δ¹⁸ Op-31.35, was developed for the southeastern USA to test the prediction accuracy of the δ¹⁸ O_tw isoscape using individuals with known residential histories. A local Mississippi resident (USAFA-134) was assigned with 90% probability to the correct region-of-origin reported by the participant. Assignments for Georgia residents (USAFA-118 and USAFA-205) had variable results, predicting USAFA-118 from Mississippi and USAFA-205 as a nonlocal resident. Stable isotope values often overlap geographically and a multi-isotope approach should be used when narrowing region(s)-of-origin(s). This study demonstrates the utility of refining isoscapes and the importance of tissue calibration in prediction assignments of human remains.

An Analysis of Precipitation Isotope Distributions across Namibia Using Historical Data

Global precipitation isoscapes based on the Global Network for Isotopes in Precipitation (GNIP) network are an important toolset that aid our understanding of global hydrologic cycles. Although the GNIP database is instrumental in developing global isoscapes, data coverage in some regions of hydrological interest (e.g., drylands) is low or non-existent thus the accuracy and relevance of global isoscapes to these regions is debatable. Capitalizing on existing literature isotope data, we generated rainfall isoscapes for Namibia (dryland) using the cokriging method and compared it to a globally fitted isoscape (GFI) downscaled to country level. Results showed weak correlation between observed and predicted isotope values in the GFI model (r² < 0.20) while the cokriging isoscape showed stronger correlation (r² = 0.67). The general trend of the local cokriging isoscape is consistent with synoptic weather systems (i.e., influences from Atlantic Ocean maritime vapour, Indian Ocean maritime vapour, Zaire Air Boundary, the Intertropical Convergence Zone and Tropical Temperate Troughs) and topography affecting the region. However, because we used the unweighted approach in this method, due to data scarcity, the absolute values could be improved in future studies. A comparison of local meteoric water lines (LMWL) constructed from the cokriging and GFI suggested that the GFI model still reflects the global average even when downscaled. The cokriging LMWL was however more consistent with expectations for an arid environment. The results indicate that although not ideal, for data deficient regions such as many drylands, the unweighted cokriging approach using historical local data can be an alternative approach to modelling rainfall isoscapes that are more relevant to the local conditions compared to using downscaled global isoscapes.

Geographically Sourcing Cocaine's Origin - Delineation of the Nineteen Major Coca Growing Regions in South America

Previously, geo-sourcing to five major coca growing regions within South America was accomplished. However, the expansion of coca cultivation throughout South America made sub-regional origin determinations increasingly difficult. The former methodology was recently enhanced with additional stable isotope analyses (²H and ¹⁸O) to fully characterize cocaine due to the varying environmental conditions in which the coca was grown. An improved data analysis method was implemented with the combination of machine learning and multivariate statistical analysis methods to provide further partitioning between growing regions. Here, we show how the combination of trace cocaine alkaloids, stable isotopes, and multivariate statistical analyses can be used to classify illicit cocaine as originating from one of 19 growing regions within South America. The data obtained through this approach can be used to describe current coca cultivation and production trends, highlight trafficking routes, as well as identify new coca growing regions.

Experimental evidence shows no fractionation of strontium isotopes ((87)Sr/(86)Sr) among soil, plants, and herbivores: implications for tracking wildlife and forensic science

Strontium isotopes ((87)Sr/(86)Sr) can be useful biological markers for a wide range of forensic science applications, including wildlife tracking. However, one of the main advantages of using (87)Sr/(86)Sr values, that there is no fractionation from geological bedrock sources through the food web, also happens to be a critical assumption that has never been tested experimentally. We test this assumption by measuring (87)Sr/(86)Sr values across three trophic levels in a controlled greenhouse experiment. Adult monarch butterflies were raised on obligate larval host milkweed plants that were, in turn, grown on seven different soil types collected across Canada. We found no significant differences between (87)Sr/(86)Sr values in leachable Sr from soil minerals, organic soil, milkweed leaves, and monarch butterfly wings. Our results suggest that strontium isoscapes developed from (87)Sr/(86)Sr values in bedrock or soil may serve as a reliable biological marker in forensic science for a range of taxa and across large geographic areas.

Sulphur isotopes in animal hair track distance to sea

Stable sulphur isotope ratios ((34)S/(32)S) in animal tissues have been suggested as a tracer of coastal residency of terrestrial animals, but data are lacking that quantify the inland range of the sulphur coastal signal and the effects of seasonality. Here, we present δ(34)S measurements of sheep wool collected seasonally on eight farms across Ireland and wool samples collected opportunistically along the west and east coasts. We observed large (>10‰) δ(34)S differences across the island and we show that wool δ(34)S values were negatively correlated with distance to the west coast. We propose that this is due to the predominantly (south-)westerly airflow, possibly combined with the influence of anthropogenic sulphur deposited from the east. While essentially all the sulphur contained in west-coast wool is of marine origin, relatively high δ(34)S values were still measured >100 km inland, suggesting that marine sulphur can be carried over long distances. Seasonal variations are small at the individual level for sedentary grazing animals. We conclude that sulphur isotopes ratios measured in archival keratinous tissues can be used to describe regional δ(34)S isoscapes primarily defined by distance to coasts and thus provide a tool to detect short-term movements of domestic, feral and wild animals within such isoscapes.