A hydrogen isoscape for tracing the migration of terrestrial herbivorous
insects across the Afro-Palearctic range
Abstract
Rationale: Many insect species undertake multi-generational migrations
in the Afro-tropical and Palearctic ranges, and understanding their
migratory connectivity remains challenging due to their small size,
short life span and large population sizes. Hydrogen isotope ( δ
2H) can be used to reconstruct the movement of
dispersing or migrating insects, but applying δ
2H for provenance requires a robust isotope baseline
map (i.e., isoscape) for the Afro-Palearctic. Methods: We analysed the
δ 2H in the wings ( δ
2H wing) of 142 resident butterflies
from 56 sites across the Afro-Palearctic. The δ
2H wing values were compared to the
predicted local growing-season precipitation δ
2H values ( δ 2H
GSP) using a linear regression model to develop an
insect wing δ 2H isoscape. We used multivariate
linear mixed models and high-resolution and time-specific remote sensing
climate and environmental data to explore the controls of the residual
δ 2H wing variability. Results:
A strong linear relationship was found between δ
2H wing and δ
2H GSP values (r
2=0.53). The resulting isoscape showed strong patterns
across the Palearctic but limited variation and high uncertainty for the
Afro-tropics. Positive residuals of this relationship were correlated
with dry conditions for the month preceding sampling whereas negative
residuals were correlated with more wet days for the month preceding
sampling. High intra-site δ 2H
wing variance was associated with lower relative
humidity for the month preceding sampling and higher elevation.
Conclusion: The δ 2H wing
isoscape is applicable to trace butterflies, moths and other terrestrial
herbivorous insects that migrate across the Afro-Palearctic range but
has limited geolocation potential in the Afro-tropics. The spatial
analysis of uncertainty using high-resolution climatic data demonstrated
that many African regions with highly variable evaporation rates and
relative humidity have δ 2H
wing values that are less related to δ
2H GSP values. Increasing geolocation
precision will require new modeling approaches using more time-specific
environmental data and/or independent geolocation tools.