Local and cross‐seasonal associations of climate and land use with abundance of monarch butterflies Danaus plexippus

Quantifying how climate and land use factors drive population dynamics at regional scales is complex because it depends on the extent of spatial and temporal synchrony among local populations, and the integration of population processes throughout a species’ annual cycle. We modeled weekly, site-specific summer abundance (1994–2013) of monarch butterflies Danaus plexippus at sites across Illinois, USA to assess relative associations of monarch abundance with climate and land use variables during the winter, spring, and summer stages of their annual cycle. We developed negative binomial regression models to estimate monarch abundance during recruitment in Illinois as a function of local climate, site-specific crop cover, and county-level herbicide (glyphosate) application. We also incorporated cross-seasonal covariates, including annual abundance of wintering monarchs in Mexico and climate conditions during spring migration and breeding in Texas, USA. We provide the first empirical evidence of a negative association between county-level glyphosate application and local abundance of adult monarchs, particularly in areas of concentrated agriculture. However, this association was only evident during the initial years of the adoption of herbicide-resistant crops (1994–2003). We also found that wetter and, to a lesser degree, cooler springs in Texas were associated with higher summer abundances in Illinois, as were relatively cool local summer temperatures in Illinois. Site-specific abundance of monarchs averaged approximately one fewer per site from 2004–2013 than during the previous decade, suggesting a recent decline in local abundance of monarch butterflies on their summer breeding grounds in Illinois. Our results demonstrate that seasonal climate and land use are associated with trends in adult monarch abundance, and our approach highlights the value of considering fine-resolution temporal fluctuations in population-level responses to environmental conditions when inferring the dynamics of migratory species.

[1]  P. Marra,et al.  SEASONAL INTERACTIONS, HABITAT QUALITY, AND POPULATION DYNAMICS IN MIGRATORY BIRDS , 2007 .

[2]  A. Peterson,et al.  Modeling current and future potential wintering distributions of eastern North American monarch butterflies , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[3]  Marc Kery,et al.  Introduction to WinBUGS for Ecologists: Bayesian approach to regression, ANOVA, mixed models and related analyses , 2010 .

[4]  K. Oberhauser,et al.  Immature Monarch Survival: Effects of Site Characteristics, Density, and Time , 2015 .

[5]  S. Gaines,et al.  Moving beyond assumptions to understand abundance distributions across the ranges of species. , 2006, Trends in ecology & evolution.

[6]  D. Slayback,et al.  Decline of monarch butterflies overwintering in Mexico: is the migratory phenomenon at risk? , 2012 .

[7]  M. Zalucki,et al.  Spatial and temporal population dynamics of monarchs downunder: Lessons for North America , 2004 .

[8]  Susan M. Haig,et al.  Links between worlds: unraveling migratory connectivity , 2002 .

[9]  Sarah P. Saunders,et al.  Evaluating confidence in climate‐based predictions of population change in a migratory species , 2016 .

[10]  M. Forister,et al.  Wherefore and Whither the Modeler: Understanding the Population Dynamics of Monarchs Will Require Integrative and Quantitative Techniques , 2016 .

[11]  Are migratory monarchs really declining in eastern North America? Examining evidence from two fall census programs , 2012 .

[12]  E. Pollard.  A method for assessing changes in the abundance of butterflies , 1977 .

[13]  I. Dworkin,et al.  Experimental evidence for within- and cross-seasonal effects of fear on survival and reproduction. , 2016, The Journal of animal ecology.

[14]  Hugh P. Possingham,et al.  Conserving mobile species , 2014 .

[15]  W. Wilhelm,et al.  Growing degree-days: one equation, two interpretations , 1997 .

[16]  Takashi Saitoh,et al.  Constraints to projecting the effects of climate change on mammals , 2006 .

[17]  D T Tyler Flockhart,et al.  Unravelling the annual cycle in a migratory animal: breeding-season habitat loss drives population declines of monarch butterflies. , 2015, The Journal of animal ecology.

[18]  P. Marra,et al.  A blind spot in climate change vulnerability assessments , 2013 .

[19]  L. Brower,et al.  Spring recolonization of eastern North America by the monarch butterfly: successive brood or single sweep migration? , 1993 .

[20]  P. Marra,et al.  Full-annual-cycle population models for migratory birds , 2015 .

[21]  L. Ries,et al.  The Disconnect Between Summer and Winter Monarch Trends for the Eastern Migratory Population: Possible Links to Differing Drivers , 2015 .

[22]  Hazel R. Parry,et al.  Monarchs in decline: a collateral landscape‐level effect of modern agriculture , 2018, Insect science.

[23]  Leslie Ries,et al.  A Citizen Army for Science: Quantifying the Contributions of Citizen Scientists to our Understanding of Monarch Butterfly Biology , 2015 .

[24]  Martyn Plummer,et al.  JAGS: A program for analysis of Bayesian graphical models using Gibbs sampling , 2003 .

[25]  Wesley W. Stone,et al.  Estimation of annual agricultural pesticide use for counties of the conterminous United States, 1992-2009 , 2013 .

[26]  Wesley W. Stone,et al.  Estimation of Annual Agricultural Pesticide Use for Counties of the Conterminous United States, 1992?2009 , 2014 .

[27]  Tara L. Crewe,et al.  Long-Term Trends in the Number of Monarch Butterflies (Lepidoptera: Nymphalidae) Counted on Fall Migration at Long Point, Ontario, Canada (1995–2014) , 2015 .

[28]  Leslie Ries,et al.  Tracking climate impacts on the migratory monarch butterfly , 2012, Global change biology.

[29]  David R. Anderson,et al.  Model Selection and Multimodel Inference , 2003 .

[30]  S. Ellner,et al.  Linking the continental migratory cycle of the monarch butterfly to understand its population decline , 2016 .

[31]  Yaxing Wei,et al.  Daymet: Daily Surface Weather Data on a 1-km Grid for North America, Version 2 , 2014 .

[32]  Nancy T. Baker.  Estimated annual agricultural pesticide use by crop group for states of the conterminous United States, 1992-2014 , 2017 .

[33]  J. Kerr,et al.  Predicting the impacts of global change on species, communities and ecosystems: it takes time , 2013 .

[34]  N. G. Best,et al.  The deviance information criterion: 12 years on , 2014 .

[35]  S. R. Stevens,et al.  Host plant pattern and variation in climate predict the location of natal grounds for migratory monarch butterflies in western North America , 2010, Journal of Insect Conservation.

[36]  Robert G. Hartzler,et al.  Reduction in common milkweed (Asclepias syriaca) occurrence in Iowa cropland from 1999 to 2009 , 2010 .

[37]  A Townsend Peterson,et al.  Ecological Niches in Sequential Generations of Eastern North American Monarch Butterflies (Lepidoptera: Danaidae): The Ecology of Migration and Likely Climate Change Implications , 2007, Environmental entomology.

[38]  P. Marra,et al.  Quantifying drivers of population dynamics for a migratory bird throughout the annual cycle , 2016, Proceedings of the Royal Society B: Biological Sciences.

[39]  Keith A. Hobson,et al.  Regional climate on the breeding grounds predicts variation in the natal origin of monarch butterflies overwintering in Mexico over 38 years , 2017, Global change biology.

[40]  Alison J. Stattersfield,et al.  Key conservation issues for migratory land- and waterbird species on the world's major flyways , 2008, Bird Conservation International.

[41]  Emily B. Cohen,et al.  A call for full annual cycle research in animal ecology , 2015, Biology Letters.

[42]  M. Zalucki.  TEMPERATURE AND RATE OF DEVELOPMENT IN DANAUS PLEXIPPUS L. AND D. CHRYSIPPUS L. (LEPIDOPTERA:NYMPHALIDAE) , 1982 .

[43]  Karen S. Oberhauser,et al.  Temporal and geographical variation in monarch densities: Citizen scientists document monarch population patterns , 2004 .

[44]  P. Marra,et al.  Range-wide effects of breeding- and nonbreeding-season climate on the abundance of a Neotropical migrant songbird. , 2011, Ecology.

[45]  J. Pleasants,et al.  Conclusion of No Decline in Summer Monarch Population Not Supported , 2016 .

[46]  K. Hobson,et al.  Natal origins of migratory monarch butterflies at wintering colonies in Mexico: new isotopic evidence. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[47]  J. Pleasants,et al.  Milkweed loss in agricultural fields because of herbicide use: effect on the monarch butterfly population , 2013 .

[48]  Ian S. Pearse,et al.  Long-Term Trends in Midwestern Milkweed Abundances and Their Relevance to Monarch Butterfly Declines , 2017 .

[49]  William A. Link,et al.  Analysis of the North American Breeding Bird Survey Using Hierarchical Models , 2011 .

[50]  K. Oberhauser,et al.  Habitat Productivity and Temporal Patterns of Monarch Butterfly Egg Densities in the Eastern United States , 2015 .

[51]  J. López-García,et al.  Trends in Deforestation and Forest Degradation after a Decade of Monitoring in the Monarch Butterfly Biosphere Reserve in Mexico , 2014, Conservation biology : the journal of the Society for Conservation Biology.

[52]  Omar Vidal,et al.  Dynamics and trends of overwintering colonies of the monarch butterfly in Mexico , 2014 .

[53]  Michael W. Gerling,et al.  National Agricultural Statistics Service: Research and Perspective on the Use of Previously Reported Data During Data Collection , 2011 .

[54]  Martin Wikelski,et al.  Going, Going, Gone: Is Animal Migration Disappearing , 2008, PLoS biology.

[55]  J. Pleasants.  Milkweed restoration in the Midwest for monarch butterfly recovery: estimates of milkweeds lost, milkweeds remaining and milkweeds that must be added to increase the monarch population , 2017 .

[56]  K. Oberhauser,et al.  A trans‐national monarch butterfly population model and implications for regional conservation priorities , 2017 .

[57]  Andrew Gelman,et al.  Data Analysis Using Regression and Multilevel/Hierarchical Models , 2006 .

[58]  H. Günthard,et al.  Quantifying the drivers of HIV transmission and prevention in men who have sex with men: a population model‐based analysis in Switzerland , 2018, HIV medicine.