Forecasting next season’s Ixodes ricinus nymphal density: the example of southern Germany 2018

Forecasting next season’s Ixodes ricinus nymphal density: the example of southern Germany 2018 Exp Appl Acarol (2018) 75:281–288 https://doi.org/10.1007/s10493-018-0267-6 Forecasting next season’s Ixodes ricinus nymphal density: the example of southern Germany 2018 1 1 2,3 Katharina Brugger  · Melanie Walter  · Lidia Chitimia‑Dobler  · 2,3,4 1 Gerhard Dobler  · Franz Rubel Received: 3 April 2018 / Accepted: 18 May 2018 / Published online: 30 May 2018 © The Author(s) 2018 Abstract The castor bean tick, Ixodes ricinus (L.) (Ixodida: Ixodidae), is the principal vector of pathogens causing tick-borne encephalitis or Lyme borreliosis in Europe. It is therefore of general interest to make an estimate of the density of I. ricinus for the whole year at the beginning of the tick season. There are two necessary conditions for making a successful prediction: a long homogeneous time series of observed tick density and a clear biological relationship between environmental predictors and tick density. A 9-year time series covering the period 2009–2017 of nymphal I. ricinus flagged at monthly intervals in southern Germany has been used. With the hypothesis that I. ricinus density is triggered by the fructification of the European beech 2 years before, the mean annual temperature of the previous year, and the current mean winter temperature (December–February), a forecast of the annual nymphal tick density has been made. Therefore, a Poisson regression model was generated resulting in an explained variance of 93.4% and an error of RMSE = 21 2 2 ticks per 100 m (annual MEAN = 260 collected ticks/100 m ). An independent verification of the forecast for the year 2017 resulted in 187 predicted versus 180 observed nymphs per 100 m . For the year 2018 a relatively high number of 443 questing I. ricinus nymphs per 100 m is forecasted, i.e., a “good” tick year. Keywords Tick-borne diseases · Tick-borne encephalitis · Lyme borreliosis · Mast seeding · Fructification Gerhard Dobler and Franz Rubel have contributed equally to this work. * Katharina Brugger katharina.brugger@vetmeduni.ac.at Institute for Veterinary Public Health, University of Veterinary Medicine Vienna, Veterinärplatz 1, 1210 Vienna, Austria Bundeswehr Institute of Microbiology, Neuherbergstraße 11, 80937 Munich, Germany German Center of Infection Research (DZIF), Partner Site Munich, Munich, Germany Parasitology Unit, University of Hohenheim, Emil-Wolff -Straße 34, 70593 Stuttgart, Germany 1 3 282 Exp Appl Acarol (2018) 75:281–288 Introduction The question whether or not the upcoming year will be a tick year in Europe, i.e., a year with many Ixodes ricinus (L.) (Ixodida: Ixodidae) ticks, comes up every spring. If reliable such a forecast would be of great interest for any nature lover, but also for dog owners. This also rises the question whether from such an annual assessment directly the risk level of getting a tick-borne disease such as tick-borne encephalitis (TBE) or Lyme borreliosis (LB) can be derived regardless of other influencing factors. With four live stages (egg, larva, nymph, and adult) and a life span of up to several years (Gray et al. 2016), a prediction on how the biological system I. ricinus will behave over the next few days, months or even years is rather complex. Each time scale has its own chal- lenges and its own accuracy, as already known from meteorology forecasts with different forecast ranges (World Meteorological Organization 2010). The weather situation and all the relevant parameters can be described in detail for the next 2 h (nowcasting) or for the next 12–72  h (short-range forecasting). Otherwise the weather situation can be described less detailed for a period beyond 3 and up to 10 days (medium-range forecast). For periods from 30 days up to 2 years (long-range forecast) averages or deviations are given. Applied for ticks, short-range forecasts of tick activity cover up to 1 week ahead and should be regularly updated. Influenced mainly by the actual weather, an estimate of tick activity including probable significant changes has to be given. Until now different meth- ods were presented to forecast tick activity in European countries. Temperature and rela- tive humidity turned out to be the main predictors for all methods. While the TICKPRO program releases a single estimate the daily tick activity for the Czech Republic (Daniel et al. 2010), the FleaTickRisk (Beugnet et al. 2009) as well as the TEKENRADAR model (Garcia-Martí et al. 2017) provide a spatial estimate of the tick activity weekly for Europe and daily for The Netherlands, respectively. Only the previous weekly forecasts for Ger- many provided by the company TICK-RADAR also considered the recent host-seeking tick activity by observing ticks in field plots (Kahl and Dautel, personal communication). Short-range forecasts of tick activity can be used to derive also the risk of acquiring a tick- borne infection for exposed people. On the other hand, long-range forecasts of tick density covering the coming months or years are still missing. Note that, for these longer time periods short-time or inter-daily fluctuations are not considered. Therefore, long-range forecasts can be interpreted as a guide for decision-making over a longer period. One important requirement for developing long-range forecasts, i.e. predictions of the annual tick density, is an existing monitoring of monthly tick activity over long periods. Although in Germany time series of regularly flagged nymphal ticks at 69 sites are available (Brugger et  al. 2016), the majority covers only 2–3 years. So far, the time series of the site Haselmühl (Bavaria, Germany) is the longest with nine consecutive years without data gaps. Time series of comparable lengths such as those from Daniel et al. (2015) or Takken et al. (2017) are seldom or not yet pub- lished (Kahl and Dautel, personal communication). The Haselmühl time series was already used to demonstrate that the seasonal inter-annual fluctuations of tick activity are most affected by time-lagged and temporal averaged variables, such as temperature and precipi- tation, rather than by contemporaneous variables (Brugger et al. 2017a). Here we use the available time series to quantify the influence of biotic and abiotic variables on the annual nymphal I. ricinus density and to forecast the next year’s density. Beside time-lagged mete- orological variables, also the fructification of European beeches (Fagus sylvatica), i.e. the annual seed production, was used. The boosting effect 2 years after a mast seeding on the 1 3 Exp Appl Acarol (2018) 75:281–288 283 nymphal density is a long-known concept, e.g for I. scapularis in the USA (Ostfeld et al. 1996). The aim of this work is to determine whether a reliable forecast of the tick density for a given location at the beginning of the tick season is possible. Which biotic and abiotic parameters can be used to forecast the next year’s nymphal tick density? Using the example of the sampling site Haselmühl, a first guess to predict the next year’s annual nymphal I. ricinus density is undertaken. Materials and methods Time series of nymphal tick densities and environmental variables For the analysis we used the data set of the site Haselmühl, a time series of monthly nym- phal ticks per 100 m for the period 2009–2017. The site as well as the flagging procedure has been described in detail by Brugger et al. (2017a). Here the main focus is on the I. rici- nus nymphs as this stage plays an important role in the epidemiology of human infec- tions (Gray 1998). Note that, the annual nymphal density is the total number of I.  ricinus nymphs monthly collected per 100 m during 1 year. Climate variables were taken from the nearest weather station Regensburg-Ober- hub (WMO No. 107760) of the German Weather Service (2018). For the period 2007–2017, time series of annual mean values as well seasonal means (i.e. Dec–Jan–Feb, Mar–April–May, June–July–Aug, Sep–Oct–Nov) for the variables temperature and precipi- tation were aggregated out of daily measurements. The density of ticks depends also on the availability of suitable hosts. Preferred hosts of I. ricinus larvae and nymphs are among others small rodents (Mihalca and Sándor 2013). As no perennial observations of small rodents are available for the study site Haselmühl, the fructification index of the European beech (Fagus sylvatica) was applied for indicating the rodent density. Beechnuts are a basic food source for small rodents resulting in popula- tion peaks 1 year after mast seeding (Ostfeld et al. 1996; Clement et al. 2009). Given the importance of forestry in Bavaria, this annual index is available back to 1954 (Kronnert et al. 2016). Fructification is defined as the annual seed production and is divided into four classes according to Eichhorn et  al. (2017): (1) absent, i.e. no fructification, (2) scarce, i.e. sporadic occurrence of fructification, but not noticeable at first sight, (3) common, i.e. clearly visible fructification, and (4) abundant, i.e. full fructification, also known as mast seeding. Statistical analysis In a first step the correlation coefficients between the annual nymphal I. ricinus density and environmental variables were calculated. As some variables are not normally distributed, the Spearman’s rank order correlation coefficient was applied. Additionally, time-lags of 1 and 2 years previous were considered to identify the highest correlation. Finally, a Pois- son regression model for the annual nymphal tick density was developed, whereas non- significant variables were removed in a stepwise procedure. Accounting the overdispersion of the time series, a quasi-Poisson GLM was used to correct the standard errors (Zuur et al. 2009). The coefficient of determination for generalised linear models R and the root mean square error (RMSE) were applied as goodness-of-fit measures. 1 3 284 Exp Appl Acarol (2018) 75:281–288 All analyses were conducted with the open-source statistical computing environment R (R Development Core Team 2018). The package rsq (Zhang 2016) was used for calculating the coefficient of determination R . Results To forecast the nymphal I.  ricinus density of the next year, a quasi-Poisson model was derived with the mean winter temperature of the months December to February T , the DJF mean annual temperature of the previous year T , as well as the beech fructification year−1 index 2 years prior F as indicated in Table 1. The model results in an explained vari- year−2 ance of 93.4% and an RMSE of 21 ticks per year. The latter is equivalent to 8% of the mean collected annual tick density of 260 nymphs/100 m . As depicted in Fig. 1, the model simulated the annual nymphal tick density with peaks almost every second year fairly accu- rately. After a year with relative low tick density below 200 nymphs per 100 m in 2017 an annual density of 443 nymphs per 100 m was estimated for 2018 (Fig.  1). This forecast Annual density simulated 237164 246196 384 202 347 378 187 443 observed 279149 260173 352 212 351 385 180 (a) 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 (b) −1 −2 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 (c) 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 (d) 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 Fig. 1 a Annual nymphal Ixodes ricinus density in Haselmühl (Germany) observed (grey bars) versus simulated (points) between 2009–2017 and forecasted for 2018 (diamond). The associated standard error is given. The explanatory variables b mean winter temperature T , c mean annual temperature of the DJF previous year T , and d the beech fructification index 2 years prior F . The dotted lines indicate year−1 year−2 the 9-year mean between 2009–2017 for the nymphal tick density and the 30-year mean of 1988–2017 for explanatory variables, respectively 1 3 F T in °C T in °C Nymphs per 100 m year−2 year−1 DJF Exp Appl Acarol (2018) 75:281–288 285 Table 1 Summary of the quasi- Estimate SE z p Poisson regression model for annual tick density Intercept 2.5404 0.7031 3.613 < 0.05 T 0.0911 0.0315 2.893 < 0.05 DJF T 0.2591 0.0706 3.672 < 0.05 year−1 F 0.2440 0.0351 6.959 < 0.001 year−2 For the explanatory variables mean winter temperature ( T ), mean DJF annual temperature of the previous year ( T ), and the beech fruc- year−1 tification index 2 years prior ( F ) the parameter estimates, the year−2 standard errors SE, the z-values (test statistics), and the p-values (sig- nificance) are given can be evaluated by independent data, when observations will be available at the end of the year 2018. In a pilot study, the same model was already used to forecast the annual nymphal tick th density for the year 2017. As presented at the 19 Annual Meeting of the International Scientific Working Group on Tick-Borne Encephalitis (Kunze and ISW-TBE 2017), a rela- tively low tick density of 187 ticks per 100 m was predicted for Haselmühl in March 2017 (Brugger et al. 2017b) and evaluated with 180 collected ticks per 100 m at the end of 2017. So far, this can be considered as the most robust model evaluation based on model-inde- pendent observations. A limitation of the model, however, is its statistical nature implying that it was fitted to the specific location of Haselmühl (Bavaria, Germany). Therefore, it is not generally applicable and needs individual adjustment when applied to other sites. Discussion Using the long-term monitoring of I.  ricinus in Haselmühl, it was demonstrated that the fructification index as well as time-lagged temperature means are appropriate pre- dictors for the annual nymphal tick density. Even more, by the end of the winter time, i.e. at the end of February, a forecast of the tick density for the upcoming year can be made. Especially, mast seeding, i.e. full fructification, turned out to be an appropriate predictor for the annual nymphal tick density. Mast seeding of beeches sets off a chain of ecological reactions (Kelly 1994; Ostfeld et  al. 1996). The effects of seeding from the perspective of ticks and tick-borne pathogens is outlined in Fig.  2. As an optimal food source, a magnitude of beech nuts attracts large mammals as roe deer (Capreo- lus capreolus) or wild boar (Sus scrofa), which are common hosts for female I.  ricinus (Gray et al. 2016). As a consequence, this boosts the density of larvae in the following summer. Coincident, the abundant beech nuts benefit an increased survival and breeding of small rodent populations mainly Myodes spp. or Apodemus spp., which are suitable hosts for I.  ricinus larvae. Together with a moderate annual temperature in the preced- ing year and warmer winter temperatures this culminates after moulting from larvae to nymphs and successful overwintering in higher nymphal densities 2 years after mast seeding. Such population peaks can result in an increased risk for humans to acquire tick-borne pathogens causing e.g. TBE or LB. The relation between mast seeding, the subsequent increase of mouse populations and disease incidence one or two years later 1 3 286 Exp Appl Acarol (2018) 75:281–288 human tick-borne infections (+) food source tick development pathogen transmission epidemiological level ecological level adult ticks larvae (+) nymphs (+) 4-6 weeks roe deer (+) beech small rodents (+) masting (+) time year year year t-2 t-1 t Fig. 2 The effect of mast seeding on Ixodes ricinus and the possible transmission of tick-borne pathogens. An increase is indicated by (+) has been already shown e.g. for human Lyme disease in the Eastern USA (Jones et  al. 1998; Schauber et al. 2005) and Poland (Bogdziewicz and Szymkowiak 2016), or Puu- mala virus infections in Germany (Reil et al. 2015). As an alternative to the fructification index, the small rodent density or hunting sta- tistics can be applied directly as biotic predictors. For Germany, several time series are available e.g. for field voles (Microtus agrestis) or bank voles (Myodes glareolus) some of which go back to 1952 (Imholt et al. 2017). However, these data sets were not used in the current study, as the observation periods available for the voles did not match those available for I. ricinus. Furthermore, the fructification index is optimal as it is synchro- nised over both large areas and climate zones (Vacchiano et al. 2017). A reliable forecasting of the annual tick density can also provide a first estimate of the risk for acquiring a tick-borne pathogen and the related human incidences. In Germany, the large majority of human TBE cases was reported each year in the south- ern federal states, especially in Bavaria (Robert Koch Institute 2018). A total of 1,361 TBE cases were registered in Bavaria between 2009–2017. The correlation between the annual tick densities in Haselmühl and the TBE cases in Bavaria is relatively high for the period 2009–2016 ( r = 0.73 , p < 0.05 , n = 8 ). However, with the increase of TBE cases in 2017 the correlation drops notably ( r = 0.46 , p = 0.21 , n = 9 ). This unexpected increase of TBE cases in 2017 not only on district level but also on federal state level is inexplicable. One possible explanation is the result of the culmination of annual cycles of different length coinciding (Zeman 2017). Nevertheless, the extraordinarly high num- ber of TBE cases reported in 2017 shows the limitations of correlative models. If a so far unidentified factor or circumstance is not considered, a reliable TBE forecast is not possible yet. One missing predictor might be hunting statistics, as already demonstrated with TBE forecasts based on statistical models for Sweden (Haemig et al. 2011) or Ger- many (Kiffner et al. 2010). On the other hand also the “good” mushroom season in 2017 could have an effect on the TBE incidences (Jaenson et al. 2012). Generally, such long-range forecasts of I. ricinus densities are essential for epidemio- logical considerations: not only for individuals at risk of acquiring a pathogen capable of causing TBE or LB, but also for public health authorities. 1 3 nymphs infections Exp Appl Acarol (2018) 75:281–288 287 Acknowledgements Open access funding provided by University of Veterinary Medicine Vienna. The tick monitoring was supported in part by a German Center of Infection Research (DZIF) Grant. The study was partly financed by Pfizer Deutschland GmbH. Please note, that both funders had no role in study design, data collection, analysis, interpretation, decision to publish, or preparation of the manuscript. The authors are very grateful to Olaf Kahl for critical reading the manuscript and his valuable hints. Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 Inter- national License (http://creat iveco mmons .org/licen ses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. References Beugnet F, Chalvet-Monfray K, Loukos H (2009) FleaTickRisk: a meteorological model developed to moni- tor and predict the activity and density of three tick species and the cat flea in Europe. Geospat Health 4:97–113. https ://doi.org/10.4081/gh.2009.213 Bogdziewicz M, Szymkowiak J (2016) Oak acorn crop and Google search volume predict Lyme disease risk in temperate Europe. Basic Appl Ecol 17(4):300–307. https ://doi.org/10.1016/j.baae.2016.01.002 Brugger K, Boehnke D, Petney T, Dobler G, Pfeffer M, Silaghi C, Schaub GA, Pinior B, Dautel H, Kahl O, Pfister K, Süss J, Rubel F (2016) A density map of the tick-borne encephalitis and Lyme borre- liosis vector Ixodes ricinus (Acari: Ixodidae) for Germany. J Med Entomol 53:1292–1302. https ://doi. org/10.1093/jme/tjw11 6 Brugger K, Walter M, Chitimia-Dobler L, Dobler G, Rubel F (2017a) Seasonal cycles of the TBE and Lyme borreliosis vector Ixodes ricinus modelled with time-lagged and interval-averaged predictors. Exp Appl Acarol 73:439–450. https ://doi.org/10.1007/s1049 3-017-0197-8 Brugger K, Walter M, Chitimia-Dobler L, Dobler G, Rubel F (2017b) Forecasting the tick-borne encephali- tis vector Ixodes ricinus in space and time, Poster presented at the 19th annual meeting of the Interna- tional Scientific Working Group on Tick-Borne Encephalitis (ISW-TBE), Vienna (Austria), 29–30 May 2017. http://iissr v1.vu-wien.ac.at/i136/BRUGG ER/pdf/ISW-TBE_2017.pdf Clement J, Vercauteren J, Verstraeten WW, Ducoffre G, Barrios JM, Vandamme AM, Maes P, Van Ranst M (2009) Relating increasing hantavirus incidences to the changing climate: the mast connection. Int J Health Geogr 8:1. https ://doi.org/10.1186/1476-072X-8-1 Daniel M, Vráblík T, Valter J, Křiž B, Danielová V (2010) The TICKPRO computer program for predicting Ixodes ricinus host-seeking activity and the warning system published on websites. Cent Eur J Public Health 18(4):230–236 Daniel M, Malỳ M, Danielová V, Kříž B, Nuttall P (2015) Abiotic predictors and annual seasonal dynam- ics of Ixodes ricinus, the major disease vector of Central Europe. Parasit Vectors 8(1):478. https ://doi. org/10.1186/s1307 1-015-1092-y Eichhorn J, Roskams P, Potočić N, Timmermann V, Ferretti M, Mues V, Szepesi A, Durrant D, Seletković I, Schröck HW, Nevalainen S, Bussotti F, Garcia P, Wulff S (2017) Part IV: Visual assessment of crown condition and damaging agents. In: UNECE ICP Forests Programme Co-ordinating Centre (ed) Manual on methods and criteria for harmonized sampling, assessment, monitoring and analysis of the effects of air pollution on forests, Thünen Institute of Forest Ecosystems, Eberswalde, Germany Garcia-Martí I, Zurita-Milla R, van Vliet AJH, Takken W (2017) Modelling and mapping tick dynamics using volunteered observations. Int J Health Geogr 16:41. https ://doi.org/10.1186/s1294 2-017-0114-8 German Weather Service (2018) Climate data Germany. http://www.dwd.de/DE/leist ungen /klima daten deuts chlan d/klima daten deuts chlan d.html. Accessed 01 Feb 2018 Gray JS (1998) The ecology of ticks transmitting Lyme borreliosis. Exp Appl Acarol 22:249–258. https :// doi.org/10.1023/A:10060 70416 135 Gray JS, Kahl O, Lane RS, Levin ML, Tsao JI (2016) Diapause in ticks of the medically important Ixodes ricinus species complex. Ticks Tick-Borne Dis 7:992–1003. https ://doi.org/10.1016/j.ttbdi s.2016.05.006 Haemig PD, Sjöstedt de Luna S, Grafström A, Lithner S, Lundkvist A, Waldenström J, Kindberg J, Stedt J, Olsén B (2011) Forecasting risk of tick-borne encephalitis (TBE): using data from wildlife and climate to predict next year’s number of human victims. Scand J Infect Dis 43:366–372. https ://doi. org/10.3109/00365 548.2011.55207 2 Imholt C, Reil D, Plašil P, Rödiger K, Jacob J (2017) Long-term population patterns of rodents and associ- ated damage in German forestry. Pest Manag Sci 73:332–340. https ://doi.org/10.1002/ps.4325 1 3 288 Exp Appl Acarol (2018) 75:281–288 Jaenson TGT, Hjertqvist M, Bergström T, Lundkvist A (2012) Why is tick-borne encephalitis increasing? A review of the key factors causing the increasing incidence of human TBE in Sweden. Parasit Vectors 5:184. https ://doi.org/10.1186/1756-3305-5-184 Jones CG, Ostfeld RS, Richard MP, Schauber EM, Wolff JO (1998) Chain reactions linking acorns to gypsy moth outbreaks and Lyme disease risk. Science 279:1023–1026 Kelly D (1994) The evolutionary ecology of mast seeding. Trends Ecol Evol 9:465–470. https ://doi. org/10.1016/0169-5347(94)90310 -7 Kiffner C, Zucchini W, Schomaker P, Vor T, Hagedorn P, Niedrig M, Rühe F (2010) Determinants of tick- borne encephalitis in counties of southern Germany, 2001–2008. Int J Health Geogr 9:42. https ://doi. org/10.1186/1476-072X-9-42 Kronnert M, Schneck D, Zollner A (2016) Blühen und Fruktifizieren unserer Waldbäume in den letzten 60 Jahren. LWF Wissen 74:37–45 Kunze U, ISW-TBE (2017) Report of the 19th Annual Meeting of the International Scientific Working Group on Tick-Borne Encephalitis (ISW-TBE)—TBE in a changing world. Ticks Tick-Borne Dis. https ://doi.org/10.1016/j.ttbdi s.2017.08.009 Mihalca AD, Sándor AD (2013) The role of rodents in the ecology of Ixodes ricinus and associated patho- gens in Central and Eastern Europe. Front Cell Infect Microbiol 3:56. https ://doi.org/10.3389/fcimb .2013.00056 Ostfeld RS, Jones CG, Wolff JO (1996) Of mice and mast. BioScience 46:323–330. https ://doi. org/10.2307/13129 46 R Development Core Team (2018) R: a language and environment for statistical computing, R Foundation for Statistical Computing, Vienna, Austria, ISBN 3-900051-07-0. Version 3.3.1. http://www.R-proje ct.org/ Reil D, Imholt C, Eccard JA, Jacob J (2015) Beech fructification and bank vole population dynamics-com- bined analyses of promoters of human Puumala virus infections in Germany. PloS ONE 10(e0134):124. https ://doi.org/10.1371/journ al.pone.01341 24 Robert Koch Institute (2018) SurvStat@RKI 2.0, https ://survs tat.rki.de. Accessed 01 Feb 2018 Schauber EM, Ostfeld RS, Evans AS Jr (2005) What is the best predictor of annual Lyme disease incidence: weather, mice, or acorns? Ecol Appl 15(2):575–586. https ://doi.org/10.1890/03-5370 Takken W, van Vliet AJH, Verhulst NO, Jacobs FHH, Gassner F, Hartemink N, Mulder S, Sprong H (2017) Acarological risk of Borrelia burgdorferi sensu lato infections across space and time in The Nether- lands. Vector Borne Zoonotic Dis 17:99–107. https ://doi.org/10.1089/vbz.2015.1933 Vacchiano G, Hacket-Pain A, Turco M, Motta R, Maringer J, Conedera M, Drobyshev I, Ascoli D (2017) Spatial patterns and broad-scale weather cues of beech mast seeding in Europe. New Phytol 215:595– 608. https ://doi.org/10.1111/nph.14600 World Meteorological Organization (2010) Manual on the global data-processing and forecasting system— Global aspects. WMO-No. 485 Zeman P (2017) Cyclic patterns in the central Europe tick-borne encephalitis incidence series. Epidemiol Infect 145:358–367. https ://doi.org/10.1017/S0950 26881 60022 23 Zhang D (2016) rsq: Coefficient of determination, R package version 0.6 Zuur AF, Ieno E, Walker N, Saveliev A, Smith G (2009) Mixed effects models and extentions in ecology with R. Statistics for biology and health. Springer, New York 1 3 http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Experimental and Applied Acarology Springer Journals

Forecasting next season’s Ixodes ricinus nymphal density: the example of southern Germany 2018

Free
8 pages
Loading next page...
 
/lp/springer_journal/forecasting-next-season-s-ixodes-ricinus-nymphal-density-the-example-engEVxJ0fo
Publisher
Springer International Publishing
Copyright
Copyright © 2018 by The Author(s)
Subject
Life Sciences; Entomology; Animal Systematics/Taxonomy/Biogeography; Animal Genetics and Genomics; Animal Ecology; Life Sciences, general
ISSN
0168-8162
eISSN
1572-9702
D.O.I.
10.1007/s10493-018-0267-6
Publisher site
See Article on Publisher Site

Abstract

Exp Appl Acarol (2018) 75:281–288 https://doi.org/10.1007/s10493-018-0267-6 Forecasting next season’s Ixodes ricinus nymphal density: the example of southern Germany 2018 1 1 2,3 Katharina Brugger  · Melanie Walter  · Lidia Chitimia‑Dobler  · 2,3,4 1 Gerhard Dobler  · Franz Rubel Received: 3 April 2018 / Accepted: 18 May 2018 / Published online: 30 May 2018 © The Author(s) 2018 Abstract The castor bean tick, Ixodes ricinus (L.) (Ixodida: Ixodidae), is the principal vector of pathogens causing tick-borne encephalitis or Lyme borreliosis in Europe. It is therefore of general interest to make an estimate of the density of I. ricinus for the whole year at the beginning of the tick season. There are two necessary conditions for making a successful prediction: a long homogeneous time series of observed tick density and a clear biological relationship between environmental predictors and tick density. A 9-year time series covering the period 2009–2017 of nymphal I. ricinus flagged at monthly intervals in southern Germany has been used. With the hypothesis that I. ricinus density is triggered by the fructification of the European beech 2 years before, the mean annual temperature of the previous year, and the current mean winter temperature (December–February), a forecast of the annual nymphal tick density has been made. Therefore, a Poisson regression model was generated resulting in an explained variance of 93.4% and an error of RMSE = 21 2 2 ticks per 100 m (annual MEAN = 260 collected ticks/100 m ). An independent verification of the forecast for the year 2017 resulted in 187 predicted versus 180 observed nymphs per 100 m . For the year 2018 a relatively high number of 443 questing I. ricinus nymphs per 100 m is forecasted, i.e., a “good” tick year. Keywords Tick-borne diseases · Tick-borne encephalitis · Lyme borreliosis · Mast seeding · Fructification Gerhard Dobler and Franz Rubel have contributed equally to this work. * Katharina Brugger katharina.brugger@vetmeduni.ac.at Institute for Veterinary Public Health, University of Veterinary Medicine Vienna, Veterinärplatz 1, 1210 Vienna, Austria Bundeswehr Institute of Microbiology, Neuherbergstraße 11, 80937 Munich, Germany German Center of Infection Research (DZIF), Partner Site Munich, Munich, Germany Parasitology Unit, University of Hohenheim, Emil-Wolff -Straße 34, 70593 Stuttgart, Germany 1 3 282 Exp Appl Acarol (2018) 75:281–288 Introduction The question whether or not the upcoming year will be a tick year in Europe, i.e., a year with many Ixodes ricinus (L.) (Ixodida: Ixodidae) ticks, comes up every spring. If reliable such a forecast would be of great interest for any nature lover, but also for dog owners. This also rises the question whether from such an annual assessment directly the risk level of getting a tick-borne disease such as tick-borne encephalitis (TBE) or Lyme borreliosis (LB) can be derived regardless of other influencing factors. With four live stages (egg, larva, nymph, and adult) and a life span of up to several years (Gray et al. 2016), a prediction on how the biological system I. ricinus will behave over the next few days, months or even years is rather complex. Each time scale has its own chal- lenges and its own accuracy, as already known from meteorology forecasts with different forecast ranges (World Meteorological Organization 2010). The weather situation and all the relevant parameters can be described in detail for the next 2 h (nowcasting) or for the next 12–72  h (short-range forecasting). Otherwise the weather situation can be described less detailed for a period beyond 3 and up to 10 days (medium-range forecast). For periods from 30 days up to 2 years (long-range forecast) averages or deviations are given. Applied for ticks, short-range forecasts of tick activity cover up to 1 week ahead and should be regularly updated. Influenced mainly by the actual weather, an estimate of tick activity including probable significant changes has to be given. Until now different meth- ods were presented to forecast tick activity in European countries. Temperature and rela- tive humidity turned out to be the main predictors for all methods. While the TICKPRO program releases a single estimate the daily tick activity for the Czech Republic (Daniel et al. 2010), the FleaTickRisk (Beugnet et al. 2009) as well as the TEKENRADAR model (Garcia-Martí et al. 2017) provide a spatial estimate of the tick activity weekly for Europe and daily for The Netherlands, respectively. Only the previous weekly forecasts for Ger- many provided by the company TICK-RADAR also considered the recent host-seeking tick activity by observing ticks in field plots (Kahl and Dautel, personal communication). Short-range forecasts of tick activity can be used to derive also the risk of acquiring a tick- borne infection for exposed people. On the other hand, long-range forecasts of tick density covering the coming months or years are still missing. Note that, for these longer time periods short-time or inter-daily fluctuations are not considered. Therefore, long-range forecasts can be interpreted as a guide for decision-making over a longer period. One important requirement for developing long-range forecasts, i.e. predictions of the annual tick density, is an existing monitoring of monthly tick activity over long periods. Although in Germany time series of regularly flagged nymphal ticks at 69 sites are available (Brugger et  al. 2016), the majority covers only 2–3 years. So far, the time series of the site Haselmühl (Bavaria, Germany) is the longest with nine consecutive years without data gaps. Time series of comparable lengths such as those from Daniel et al. (2015) or Takken et al. (2017) are seldom or not yet pub- lished (Kahl and Dautel, personal communication). The Haselmühl time series was already used to demonstrate that the seasonal inter-annual fluctuations of tick activity are most affected by time-lagged and temporal averaged variables, such as temperature and precipi- tation, rather than by contemporaneous variables (Brugger et al. 2017a). Here we use the available time series to quantify the influence of biotic and abiotic variables on the annual nymphal I. ricinus density and to forecast the next year’s density. Beside time-lagged mete- orological variables, also the fructification of European beeches (Fagus sylvatica), i.e. the annual seed production, was used. The boosting effect 2 years after a mast seeding on the 1 3 Exp Appl Acarol (2018) 75:281–288 283 nymphal density is a long-known concept, e.g for I. scapularis in the USA (Ostfeld et al. 1996). The aim of this work is to determine whether a reliable forecast of the tick density for a given location at the beginning of the tick season is possible. Which biotic and abiotic parameters can be used to forecast the next year’s nymphal tick density? Using the example of the sampling site Haselmühl, a first guess to predict the next year’s annual nymphal I. ricinus density is undertaken. Materials and methods Time series of nymphal tick densities and environmental variables For the analysis we used the data set of the site Haselmühl, a time series of monthly nym- phal ticks per 100 m for the period 2009–2017. The site as well as the flagging procedure has been described in detail by Brugger et al. (2017a). Here the main focus is on the I. rici- nus nymphs as this stage plays an important role in the epidemiology of human infec- tions (Gray 1998). Note that, the annual nymphal density is the total number of I.  ricinus nymphs monthly collected per 100 m during 1 year. Climate variables were taken from the nearest weather station Regensburg-Ober- hub (WMO No. 107760) of the German Weather Service (2018). For the period 2007–2017, time series of annual mean values as well seasonal means (i.e. Dec–Jan–Feb, Mar–April–May, June–July–Aug, Sep–Oct–Nov) for the variables temperature and precipi- tation were aggregated out of daily measurements. The density of ticks depends also on the availability of suitable hosts. Preferred hosts of I. ricinus larvae and nymphs are among others small rodents (Mihalca and Sándor 2013). As no perennial observations of small rodents are available for the study site Haselmühl, the fructification index of the European beech (Fagus sylvatica) was applied for indicating the rodent density. Beechnuts are a basic food source for small rodents resulting in popula- tion peaks 1 year after mast seeding (Ostfeld et al. 1996; Clement et al. 2009). Given the importance of forestry in Bavaria, this annual index is available back to 1954 (Kronnert et al. 2016). Fructification is defined as the annual seed production and is divided into four classes according to Eichhorn et  al. (2017): (1) absent, i.e. no fructification, (2) scarce, i.e. sporadic occurrence of fructification, but not noticeable at first sight, (3) common, i.e. clearly visible fructification, and (4) abundant, i.e. full fructification, also known as mast seeding. Statistical analysis In a first step the correlation coefficients between the annual nymphal I. ricinus density and environmental variables were calculated. As some variables are not normally distributed, the Spearman’s rank order correlation coefficient was applied. Additionally, time-lags of 1 and 2 years previous were considered to identify the highest correlation. Finally, a Pois- son regression model for the annual nymphal tick density was developed, whereas non- significant variables were removed in a stepwise procedure. Accounting the overdispersion of the time series, a quasi-Poisson GLM was used to correct the standard errors (Zuur et al. 2009). The coefficient of determination for generalised linear models R and the root mean square error (RMSE) were applied as goodness-of-fit measures. 1 3 284 Exp Appl Acarol (2018) 75:281–288 All analyses were conducted with the open-source statistical computing environment R (R Development Core Team 2018). The package rsq (Zhang 2016) was used for calculating the coefficient of determination R . Results To forecast the nymphal I.  ricinus density of the next year, a quasi-Poisson model was derived with the mean winter temperature of the months December to February T , the DJF mean annual temperature of the previous year T , as well as the beech fructification year−1 index 2 years prior F as indicated in Table 1. The model results in an explained vari- year−2 ance of 93.4% and an RMSE of 21 ticks per year. The latter is equivalent to 8% of the mean collected annual tick density of 260 nymphs/100 m . As depicted in Fig. 1, the model simulated the annual nymphal tick density with peaks almost every second year fairly accu- rately. After a year with relative low tick density below 200 nymphs per 100 m in 2017 an annual density of 443 nymphs per 100 m was estimated for 2018 (Fig.  1). This forecast Annual density simulated 237164 246196 384 202 347 378 187 443 observed 279149 260173 352 212 351 385 180 (a) 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 (b) −1 −2 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 (c) 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 (d) 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 Fig. 1 a Annual nymphal Ixodes ricinus density in Haselmühl (Germany) observed (grey bars) versus simulated (points) between 2009–2017 and forecasted for 2018 (diamond). The associated standard error is given. The explanatory variables b mean winter temperature T , c mean annual temperature of the DJF previous year T , and d the beech fructification index 2 years prior F . The dotted lines indicate year−1 year−2 the 9-year mean between 2009–2017 for the nymphal tick density and the 30-year mean of 1988–2017 for explanatory variables, respectively 1 3 F T in °C T in °C Nymphs per 100 m year−2 year−1 DJF Exp Appl Acarol (2018) 75:281–288 285 Table 1 Summary of the quasi- Estimate SE z p Poisson regression model for annual tick density Intercept 2.5404 0.7031 3.613 < 0.05 T 0.0911 0.0315 2.893 < 0.05 DJF T 0.2591 0.0706 3.672 < 0.05 year−1 F 0.2440 0.0351 6.959 < 0.001 year−2 For the explanatory variables mean winter temperature ( T ), mean DJF annual temperature of the previous year ( T ), and the beech fruc- year−1 tification index 2 years prior ( F ) the parameter estimates, the year−2 standard errors SE, the z-values (test statistics), and the p-values (sig- nificance) are given can be evaluated by independent data, when observations will be available at the end of the year 2018. In a pilot study, the same model was already used to forecast the annual nymphal tick th density for the year 2017. As presented at the 19 Annual Meeting of the International Scientific Working Group on Tick-Borne Encephalitis (Kunze and ISW-TBE 2017), a rela- tively low tick density of 187 ticks per 100 m was predicted for Haselmühl in March 2017 (Brugger et al. 2017b) and evaluated with 180 collected ticks per 100 m at the end of 2017. So far, this can be considered as the most robust model evaluation based on model-inde- pendent observations. A limitation of the model, however, is its statistical nature implying that it was fitted to the specific location of Haselmühl (Bavaria, Germany). Therefore, it is not generally applicable and needs individual adjustment when applied to other sites. Discussion Using the long-term monitoring of I.  ricinus in Haselmühl, it was demonstrated that the fructification index as well as time-lagged temperature means are appropriate pre- dictors for the annual nymphal tick density. Even more, by the end of the winter time, i.e. at the end of February, a forecast of the tick density for the upcoming year can be made. Especially, mast seeding, i.e. full fructification, turned out to be an appropriate predictor for the annual nymphal tick density. Mast seeding of beeches sets off a chain of ecological reactions (Kelly 1994; Ostfeld et  al. 1996). The effects of seeding from the perspective of ticks and tick-borne pathogens is outlined in Fig.  2. As an optimal food source, a magnitude of beech nuts attracts large mammals as roe deer (Capreo- lus capreolus) or wild boar (Sus scrofa), which are common hosts for female I.  ricinus (Gray et al. 2016). As a consequence, this boosts the density of larvae in the following summer. Coincident, the abundant beech nuts benefit an increased survival and breeding of small rodent populations mainly Myodes spp. or Apodemus spp., which are suitable hosts for I.  ricinus larvae. Together with a moderate annual temperature in the preced- ing year and warmer winter temperatures this culminates after moulting from larvae to nymphs and successful overwintering in higher nymphal densities 2 years after mast seeding. Such population peaks can result in an increased risk for humans to acquire tick-borne pathogens causing e.g. TBE or LB. The relation between mast seeding, the subsequent increase of mouse populations and disease incidence one or two years later 1 3 286 Exp Appl Acarol (2018) 75:281–288 human tick-borne infections (+) food source tick development pathogen transmission epidemiological level ecological level adult ticks larvae (+) nymphs (+) 4-6 weeks roe deer (+) beech small rodents (+) masting (+) time year year year t-2 t-1 t Fig. 2 The effect of mast seeding on Ixodes ricinus and the possible transmission of tick-borne pathogens. An increase is indicated by (+) has been already shown e.g. for human Lyme disease in the Eastern USA (Jones et  al. 1998; Schauber et al. 2005) and Poland (Bogdziewicz and Szymkowiak 2016), or Puu- mala virus infections in Germany (Reil et al. 2015). As an alternative to the fructification index, the small rodent density or hunting sta- tistics can be applied directly as biotic predictors. For Germany, several time series are available e.g. for field voles (Microtus agrestis) or bank voles (Myodes glareolus) some of which go back to 1952 (Imholt et al. 2017). However, these data sets were not used in the current study, as the observation periods available for the voles did not match those available for I. ricinus. Furthermore, the fructification index is optimal as it is synchro- nised over both large areas and climate zones (Vacchiano et al. 2017). A reliable forecasting of the annual tick density can also provide a first estimate of the risk for acquiring a tick-borne pathogen and the related human incidences. In Germany, the large majority of human TBE cases was reported each year in the south- ern federal states, especially in Bavaria (Robert Koch Institute 2018). A total of 1,361 TBE cases were registered in Bavaria between 2009–2017. The correlation between the annual tick densities in Haselmühl and the TBE cases in Bavaria is relatively high for the period 2009–2016 ( r = 0.73 , p < 0.05 , n = 8 ). However, with the increase of TBE cases in 2017 the correlation drops notably ( r = 0.46 , p = 0.21 , n = 9 ). This unexpected increase of TBE cases in 2017 not only on district level but also on federal state level is inexplicable. One possible explanation is the result of the culmination of annual cycles of different length coinciding (Zeman 2017). Nevertheless, the extraordinarly high num- ber of TBE cases reported in 2017 shows the limitations of correlative models. If a so far unidentified factor or circumstance is not considered, a reliable TBE forecast is not possible yet. One missing predictor might be hunting statistics, as already demonstrated with TBE forecasts based on statistical models for Sweden (Haemig et al. 2011) or Ger- many (Kiffner et al. 2010). On the other hand also the “good” mushroom season in 2017 could have an effect on the TBE incidences (Jaenson et al. 2012). Generally, such long-range forecasts of I. ricinus densities are essential for epidemio- logical considerations: not only for individuals at risk of acquiring a pathogen capable of causing TBE or LB, but also for public health authorities. 1 3 nymphs infections Exp Appl Acarol (2018) 75:281–288 287 Acknowledgements Open access funding provided by University of Veterinary Medicine Vienna. The tick monitoring was supported in part by a German Center of Infection Research (DZIF) Grant. The study was partly financed by Pfizer Deutschland GmbH. Please note, that both funders had no role in study design, data collection, analysis, interpretation, decision to publish, or preparation of the manuscript. The authors are very grateful to Olaf Kahl for critical reading the manuscript and his valuable hints. Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 Inter- national License (http://creat iveco mmons .org/licen ses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. References Beugnet F, Chalvet-Monfray K, Loukos H (2009) FleaTickRisk: a meteorological model developed to moni- tor and predict the activity and density of three tick species and the cat flea in Europe. Geospat Health 4:97–113. https ://doi.org/10.4081/gh.2009.213 Bogdziewicz M, Szymkowiak J (2016) Oak acorn crop and Google search volume predict Lyme disease risk in temperate Europe. Basic Appl Ecol 17(4):300–307. https ://doi.org/10.1016/j.baae.2016.01.002 Brugger K, Boehnke D, Petney T, Dobler G, Pfeffer M, Silaghi C, Schaub GA, Pinior B, Dautel H, Kahl O, Pfister K, Süss J, Rubel F (2016) A density map of the tick-borne encephalitis and Lyme borre- liosis vector Ixodes ricinus (Acari: Ixodidae) for Germany. J Med Entomol 53:1292–1302. https ://doi. org/10.1093/jme/tjw11 6 Brugger K, Walter M, Chitimia-Dobler L, Dobler G, Rubel F (2017a) Seasonal cycles of the TBE and Lyme borreliosis vector Ixodes ricinus modelled with time-lagged and interval-averaged predictors. Exp Appl Acarol 73:439–450. https ://doi.org/10.1007/s1049 3-017-0197-8 Brugger K, Walter M, Chitimia-Dobler L, Dobler G, Rubel F (2017b) Forecasting the tick-borne encephali- tis vector Ixodes ricinus in space and time, Poster presented at the 19th annual meeting of the Interna- tional Scientific Working Group on Tick-Borne Encephalitis (ISW-TBE), Vienna (Austria), 29–30 May 2017. http://iissr v1.vu-wien.ac.at/i136/BRUGG ER/pdf/ISW-TBE_2017.pdf Clement J, Vercauteren J, Verstraeten WW, Ducoffre G, Barrios JM, Vandamme AM, Maes P, Van Ranst M (2009) Relating increasing hantavirus incidences to the changing climate: the mast connection. Int J Health Geogr 8:1. https ://doi.org/10.1186/1476-072X-8-1 Daniel M, Vráblík T, Valter J, Křiž B, Danielová V (2010) The TICKPRO computer program for predicting Ixodes ricinus host-seeking activity and the warning system published on websites. Cent Eur J Public Health 18(4):230–236 Daniel M, Malỳ M, Danielová V, Kříž B, Nuttall P (2015) Abiotic predictors and annual seasonal dynam- ics of Ixodes ricinus, the major disease vector of Central Europe. Parasit Vectors 8(1):478. https ://doi. org/10.1186/s1307 1-015-1092-y Eichhorn J, Roskams P, Potočić N, Timmermann V, Ferretti M, Mues V, Szepesi A, Durrant D, Seletković I, Schröck HW, Nevalainen S, Bussotti F, Garcia P, Wulff S (2017) Part IV: Visual assessment of crown condition and damaging agents. In: UNECE ICP Forests Programme Co-ordinating Centre (ed) Manual on methods and criteria for harmonized sampling, assessment, monitoring and analysis of the effects of air pollution on forests, Thünen Institute of Forest Ecosystems, Eberswalde, Germany Garcia-Martí I, Zurita-Milla R, van Vliet AJH, Takken W (2017) Modelling and mapping tick dynamics using volunteered observations. Int J Health Geogr 16:41. https ://doi.org/10.1186/s1294 2-017-0114-8 German Weather Service (2018) Climate data Germany. http://www.dwd.de/DE/leist ungen /klima daten deuts chlan d/klima daten deuts chlan d.html. Accessed 01 Feb 2018 Gray JS (1998) The ecology of ticks transmitting Lyme borreliosis. Exp Appl Acarol 22:249–258. https :// doi.org/10.1023/A:10060 70416 135 Gray JS, Kahl O, Lane RS, Levin ML, Tsao JI (2016) Diapause in ticks of the medically important Ixodes ricinus species complex. Ticks Tick-Borne Dis 7:992–1003. https ://doi.org/10.1016/j.ttbdi s.2016.05.006 Haemig PD, Sjöstedt de Luna S, Grafström A, Lithner S, Lundkvist A, Waldenström J, Kindberg J, Stedt J, Olsén B (2011) Forecasting risk of tick-borne encephalitis (TBE): using data from wildlife and climate to predict next year’s number of human victims. Scand J Infect Dis 43:366–372. https ://doi. org/10.3109/00365 548.2011.55207 2 Imholt C, Reil D, Plašil P, Rödiger K, Jacob J (2017) Long-term population patterns of rodents and associ- ated damage in German forestry. Pest Manag Sci 73:332–340. https ://doi.org/10.1002/ps.4325 1 3 288 Exp Appl Acarol (2018) 75:281–288 Jaenson TGT, Hjertqvist M, Bergström T, Lundkvist A (2012) Why is tick-borne encephalitis increasing? A review of the key factors causing the increasing incidence of human TBE in Sweden. Parasit Vectors 5:184. https ://doi.org/10.1186/1756-3305-5-184 Jones CG, Ostfeld RS, Richard MP, Schauber EM, Wolff JO (1998) Chain reactions linking acorns to gypsy moth outbreaks and Lyme disease risk. Science 279:1023–1026 Kelly D (1994) The evolutionary ecology of mast seeding. Trends Ecol Evol 9:465–470. https ://doi. org/10.1016/0169-5347(94)90310 -7 Kiffner C, Zucchini W, Schomaker P, Vor T, Hagedorn P, Niedrig M, Rühe F (2010) Determinants of tick- borne encephalitis in counties of southern Germany, 2001–2008. Int J Health Geogr 9:42. https ://doi. org/10.1186/1476-072X-9-42 Kronnert M, Schneck D, Zollner A (2016) Blühen und Fruktifizieren unserer Waldbäume in den letzten 60 Jahren. LWF Wissen 74:37–45 Kunze U, ISW-TBE (2017) Report of the 19th Annual Meeting of the International Scientific Working Group on Tick-Borne Encephalitis (ISW-TBE)—TBE in a changing world. Ticks Tick-Borne Dis. https ://doi.org/10.1016/j.ttbdi s.2017.08.009 Mihalca AD, Sándor AD (2013) The role of rodents in the ecology of Ixodes ricinus and associated patho- gens in Central and Eastern Europe. Front Cell Infect Microbiol 3:56. https ://doi.org/10.3389/fcimb .2013.00056 Ostfeld RS, Jones CG, Wolff JO (1996) Of mice and mast. BioScience 46:323–330. https ://doi. org/10.2307/13129 46 R Development Core Team (2018) R: a language and environment for statistical computing, R Foundation for Statistical Computing, Vienna, Austria, ISBN 3-900051-07-0. Version 3.3.1. http://www.R-proje ct.org/ Reil D, Imholt C, Eccard JA, Jacob J (2015) Beech fructification and bank vole population dynamics-com- bined analyses of promoters of human Puumala virus infections in Germany. PloS ONE 10(e0134):124. https ://doi.org/10.1371/journ al.pone.01341 24 Robert Koch Institute (2018) SurvStat@RKI 2.0, https ://survs tat.rki.de. Accessed 01 Feb 2018 Schauber EM, Ostfeld RS, Evans AS Jr (2005) What is the best predictor of annual Lyme disease incidence: weather, mice, or acorns? Ecol Appl 15(2):575–586. https ://doi.org/10.1890/03-5370 Takken W, van Vliet AJH, Verhulst NO, Jacobs FHH, Gassner F, Hartemink N, Mulder S, Sprong H (2017) Acarological risk of Borrelia burgdorferi sensu lato infections across space and time in The Nether- lands. Vector Borne Zoonotic Dis 17:99–107. https ://doi.org/10.1089/vbz.2015.1933 Vacchiano G, Hacket-Pain A, Turco M, Motta R, Maringer J, Conedera M, Drobyshev I, Ascoli D (2017) Spatial patterns and broad-scale weather cues of beech mast seeding in Europe. New Phytol 215:595– 608. https ://doi.org/10.1111/nph.14600 World Meteorological Organization (2010) Manual on the global data-processing and forecasting system— Global aspects. WMO-No. 485 Zeman P (2017) Cyclic patterns in the central Europe tick-borne encephalitis incidence series. Epidemiol Infect 145:358–367. https ://doi.org/10.1017/S0950 26881 60022 23 Zhang D (2016) rsq: Coefficient of determination, R package version 0.6 Zuur AF, Ieno E, Walker N, Saveliev A, Smith G (2009) Mixed effects models and extentions in ecology with R. Statistics for biology and health. Springer, New York 1 3

Journal

Experimental and Applied AcarologySpringer Journals

Published: May 30, 2018

References

You’re reading a free preview. Subscribe to read the entire article.


DeepDyve is your
personal research library

It’s your single place to instantly
discover and read the research
that matters to you.

Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.

All for just $49/month

Explore the DeepDyve Library

Search

Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly

Organize

Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.

Access

Get unlimited, online access to over 18 million full-text articles from more than 15,000 scientific journals.

Your journals are on DeepDyve

Read from thousands of the leading scholarly journals from SpringerNature, Elsevier, Wiley-Blackwell, Oxford University Press and more.

All the latest content is available, no embargo periods.

See the journals in your area

DeepDyve

Freelancer

DeepDyve

Pro

Price

FREE

$49/month
$360/year

Save searches from
Google Scholar,
PubMed

Create lists to
organize your research

Export lists, citations

Read DeepDyve articles

Abstract access only

Unlimited access to over
18 million full-text articles

Print

20 pages / month

PDF Discount

20% off