More stories

  • in

    Towards an absolute light pollution indicator

    DefinitionWe present here a new statistical approach to measure and characterize light pollution. The objective is to define an indicator which is not limited to clear sky measurements and does not require a precise calibration of a photometer. The key attributes of the indicator are the following:

    It requires the automated acquisition of a large number of zenithal NSB measures when the Sun is below (-18^{circ }) and the Moon below (-5^{circ });

    The acquisitions must at least cover a period of 6 months in order to record a wide range of possible weather conditions from perfectly clear to totally overcast skies. The objective is to obtain a significant sample of every type of cloud conditions (e.g. cloud density and ceiling altitude) as well as a good characterization of the average clear sky ;

    It is based on the analysis of the zenithal NSB measure dispersion which is directly linked to the level of light pollution a site experiences.

    As presented above in Fig. 2, the NSB density histograms, which are assembled from a large number of NSB measures, display a higher density zone which denotes a characteristic clear sky level that we name nominal NSB in the scope of the indicator calculation. On both sides of the clear sky level (above and below), NSB measures are distributed in a way that reflect the zenithal night sky luminance in cloudy conditions: NSB measures above the clear sky level mean that the light pollution is amplified by clouds while those below the clear sky level indicate a darker environment where clouds mask light pollution from distant sources as well as natural light sources. The calculation of the indicator is based on the evaluation of the NSB measure dispersion on both sides of the nominal NSB (i.e. characteristic clear sky level). Since there can have strong variations of artificial light emitted into the environment at the beginning and end of the night (decrease then increase of human activity, extinction of public lighting, etc.), the range of NSB measures retained for calculating the indicator is restricted to a portion in the middle of the night, typically 2 h.Figure 8 shows a typical NSB density histogram for a site which is quite severely impacted by light pollution. It covers a 2 h time range between 23:00 UTC and 01:00 UTC and one can easily see that the zone above the nominal NSB is much higher and denser than the one below, i.e. cloud conditions create more often a brighter environment than a darker one and with a greater amplitude.Figure 8NSB density histogram where the nominal NSB that represents the most common clear sky conditions is identified. It delimits two areas, the NSB bright dispersion above the nominal NSB and the NSB dark dispersion below.Full size imageBased on the determination of the nominal NSB, a quantitative indicator, called NDR for NSB Dispersion Ratio, is calculated in the following way:$$begin{aligned} NDR = (N_b cdot MAD_b) / (N_d cdot MAD_d) end{aligned}$$where (N_b) is the number of measures above the nominal NSB (brighter sky), (N_d) is the number of measures below the nominal NSB (darker sky), (MAD_b) is the median absolute deviation of the measures in the bright dispersion zone (above the nominal NSB) and (MAD_d) is the median absolute deviation of the measures in the dark dispersion zone (below the nominal NSB). The median absolute deviation is a statistical tool used to measure the variability of a data set, which is exactly what we try to achieve with the two NSB extensions above and below the nominal NSB. It is formally defined as (MAD = median(|X_i – tilde{{mathbf {X}}}|)) where (X_i) in our case represents an NSB value and (tilde{{mathbf {X}}}) is (median(X_i)). The median absolute deviation is a better choice than the usual standard deviation to measure the spread of NSB measures since the data does not follow a normal distribution.In order to make the determination of the NSB Dispersion Ratio stronger from a statistical standpoint, we use a bootstrapping with replacement resampling method on the set of night portions used to compute the indicator. Assuming we have N night portions at our disposal, we randomly select a sample of N items in this set of night portions knowing that a given item can appear multiple times in the sample (hence the bootstrapping with replacement). The NDR value is then computed for the considered sample. This process is repeated 1000 times and the average NDR value if eventually computed. This average value represents the actual NDR indicator of the considered site.The NDR indicator takes into account both the number of NSB values on each side of the nominal NSB and the dispersion of these values. This is what makes it relevant as an indicator of light pollution which encompasses all kinds of meteorological conditions experienced at a particular site. On that aspect, it is therefore not an astronomical light pollution indicator since it is not focused on clear sky conditions. On the opposite, it requires to have a significant number of NSB measures in all sorts of cloudy conditions so that a valid NDR indicator can be derived.A key aspect of the NDR calculation methodology is to determine the level of the nominal NSB, i.e. the typical clear sky level, since it will be used to differentiate the NSB measures that go in each of the two sets to calculate the bright and dark dispersions. As we have seen earlier in the article, such a determination can be biased by natural light sources that raise or lower clear sky NSB at different times of the night. This can result into a “blurry” high density zone which makes the determination of the nominal NSB difficult or even impossible depending on the observation period. Based on the quantitative estimate of the different natural light sources presented above, the most important bias to address is the contribution of the Galactic plane. This contribution must be eliminated for all the NSB measures which are used to calculate the NDR indicator. In order to do that, Noxi, the Ninox processing software developed by DarkSkyLab, calculates for each NSB measure the corresponding Galactic plane and star fields contribution using the galactic coordinates of the zenith and integrating the combined flux of all stars in the field of view using the UCAC4 astrometry and photometry star catalogue. It is not possible to cancel the contribution of the airglow due to its unpredictable nature, but since it only appears in rare occasions, it is not seen as a problem and is ignored. Regarding the contribution of the zodiacal light, it is considered as minimal at the zenith and it is also ignored.As an example, Fig. 9 shows on the left an NSB density histogram where the Galactic plane bias has not been corrected in the data, and on the right the same data but with the Galactic plane bias corrected. It is easy to see that in the latter the nominal NSB is much easier to determine, providing a more accurate reference level to calculate the NDR indicator. Once the Galactic plane bias has been corrected, the nominal NSB is determined as the highest density zone of the NSB histogram. It must be noted that, as of today, all the NSB measures are corrected from the Galactic plane bias without regards to the presence of clouds or high levels of light pollution. This results into an additional source of inaccuracy that will be addressed in the future through the implementation of two heuristics within the Noxi software:

    1.

    A first heuristic will determine if a night portion is considered as having a clear sky or not so that the Galactic bias correction is applied only if the sky is clear. In order to do that, we have developed an indicator called the NSS (for Night Sky Stability). To determine the NSS for a full night of measures or just a night portion, we fit the NSB curve with a degree 10 polynomial and we then compute the difference between each NSB measure and it polynomial counterpart. As a result, we obtain a set of residuals. The variance of all the residuals defines the NSS for the considered NSB dataset. Below a given value, the sky is considered as clear knowing that the NSS indicator has been calibrated on several NSB data sets for which the corresponding weather conditions are known;

    2.

    A second heuristic will allow us to weight the Galactic bias correction to be applied to NSB measures according to their value. For non-polluted skies with high values of NSB, the full Galactic bias correction will be applied while below a certain NSB threshold (for instance 21 mag(_{mathrm{SQM}})/arcsec(^{2}) which corresponds roughly to the brightest parts of the Milky Way) no correction will be applied.

    Figure 9NSB density histograms of the same data set with no correction of the Galactic plane bias applied on the left and a full correction applied on the right.Full size imageThe NDR indicator is unitless since it is the ratio of two quantities with the same unit (mag(_{mathrm{SQM}})/arcsec(^{2})). For the data set presented in Fig. 9, the NDR value which is obtained is 25 (which is justified by the fact that the bright extension in the density histogram is much higher and denser than the dark extension). This denotes a quite high level of light pollution despite the fact that the nominal NSB is at a level of 21.6 mag(_{mathrm{SQM}})/arcsec(^{2}). This highlights the fact that there is not always a strict correlation between the typical clear sky NSB obtained for a given site and its NDR indicator, i.e. the presence of clouds decreases the NSB more than we could have expected just by knowing the clear sky NSB. On that respect, the NDR ratio brings more information that the clear sky NSB alone.In addition to provide an indicator which is representative of light pollution in all possible atmospheric conditions, the NDR provides a tool to compare locations in a more meaningful way than just using a set of standalone NSB evaluations. First it is not dependent of an inter-calibration between different systems and second its statistical nature makes it more robust when it comes to perform comparisons.NDR into practiceThe NDR indicator has been calculated for several different sites by DarkSkyLab during various projects in France that involved NSB measuring sessions in the field. To demonstrate some of the results that have been obtained, Fig. 10 provides the density histograms of 4 different sites which have quite different light pollution profiles.Figure 10NSB density histograms of 4 different sites used to compute the NDR indicator. The nominal NSB (which corresponds to the most common clear sky conditions) is noted with a white tick mark next to the vertical axis. Relative levels of the bright and dark dispersion terms (((N_b cdot MAD_b)) and ((N_d cdot MAD_d))) are noted respectively with an orange tick mark and a green tick mark. The computed values of the NDR indicator and nominal NSB are provided in the top-left corner of each figure.Full size imageTo build these diagrams, only the measures acquired during a few hours in the middle of the nights have been used to ensure the maximum stability of the NSB curves and avoid lighting extinctions that create large gaps in NSB profiles. The Galactic plane bias is corrected on all plots and the same NSB scale is used in order to perform comparisons between the 4 sites. One can notice that the number of measures and nights for the 4 sites are quite different. However, they are all sufficient to derive a meaningful value of the NDR indicator using the bootstrapping with replacement resampling method described above, but it is clear that the more NSB measures used, the more accurate the NDR indicator.Table 2 summarizes the NDR indicators as well as the nominal NSB for the 4 sites which are sorted in the order of decreasing NDR indicator values.Table 2 Summary of the nominal NSB and NDR indicators of the 4 different sites.Full size tableOne can see that the NDR indicator values are not strictly correlated to the nominal NSB values, e.g. despite the fact that the nominal NSB of site (a) is slightly better than the one of site (b), the NDR indicator value is much larger for site (a) than for site (b). This can be explained if we consider the specificities of each site:

    Cervières (a) is a small village in the Haut-Forez area, France, which is surrounded by large cities (Lyon, Saint-Etienne and Clermont-Ferrand at a distance between 50 to 80 km) and a closer mid-size city (Roanne at 30 km). At the top of that, the town of Noirétable and a large highway rest area are just 2 km away without any nocturnal extinction applied (as opposed to the village of Cervières itself for which public lighting is turned off from 23:00 to 05:00 local time). These conditions are favourable to the presence of a constant light pollution background which has a negative impact on the zenithal NSB measures in most cloudy conditions (distant large cities for high elevation clouds and Noirétable and the highway rest area for lower elevation clouds). Only rare cloud conditions actually protect the site from the effect of mid-distance light sources. In clear sky conditions, however, the fact that there is no close light sources provides reasonably good NSB levels;

    The Copernic Association Observatory (b) is located 6 km from the large town of Gap in the mountain area of Hautes-Alpes in the south of France. There is no significant short distant light sources but in many cloud conditions the contribution of Gap has a very negative impact on the zenithal luminance. However, due to the fact that the observatory is at a higher altitude on the hills surrounding the city of Gap, there are cloud conditions that make the site darker. In clear sky conditions, the proximity of Gap does not permit a quality better than that of a rural sky;

    The Astrièves Observatory (c) is located near the center of the small town of Gresse-en-Vercors in the Parc Naturel Régional du Vercors. There is a full nocturnal extinction of the village for a large part of the night resulting in a good sky quality in clear sky conditions. The large city of Grenoble is at a distance of 30 km in a valley at the north-east, and the two locations are separated by a few mountains which efficiently help masking the light pollution as soon as the cloud ceiling is below a certain altitude, resulting into a dark environment. On the opposite, high elevation clouds reflect the light from Grenoble and increase the zenithal luminance;

    Eourres (d) is a small and isolated village located 20 km west of Sisteron in the department of Hautes-Alpes, France, which is surrounded by mountains. There is no significant light sources closer than those of Sisteron and this results into a very good night sky quality with, most of the times, a very dark environment in cloudy conditions.

    Figure 11 provides a graphical representation of the NDR indicator values for the 4 sites. On the NDR scale, the value 1 indicates that the bright and dark dispersion terms (respectively ((N_b cdot MAD_b)) and ((N_d cdot MAD_d))) are equal, which means there is a balance between dark and bright conditions at the zenith on the considered site with reference to the most common clear sky level.Figure 11Summary of the NDR indicators obtained for the 4 sites. The diagram uses 1 as the pivotal value to delineate sites according to the two bright and dark dispersion terms ((N_b cdot MAD_b)) and ((N_d cdot MAD_d)).Full size imageThe NDR can theoretically vary between 0 (totally dark site) and several hundreds (extremely bright site) but in practice the best sites can reach NDR indicator values down to 0.3 in the best preserved locations and up to 200 for very large and polluted cities.Robustness of the NDR indicatorIt is important to evaluate how the NDR indicator is dependant on the number of measures used to compute it and to figure out what would be the minimum number of night sessions required to obtain a meaningful NDR indicator value at a given site. To achieve that, we have used the data from two of the four sites presented above (the two which have the largest number or recorded nights: Cervières with 424 nights and the Astrièves Observatory with 373 nights). The 1000-step bootstrapping procedure has been repeatedly executed on each data set with a regularly decreasing sample of nights: starting from the full number of nights, a decrement of 10 nights is applied at each step until only 20 nights are remaining. At every bootstrap step, each sample is composed of n nights randomly chosen among the N available ones knowing that any night can be selected several times.Figure 12 shows the NDR indicator values that have been obtained for each of the two sites as a function of the night sample considered. The 95% confidence interval is plotted against each NDR indicator value (it is preferred to the standard deviation since the NSB distribution in the data sets is not normal). In the right plot of Fig. 12, the last confidence interval for the 24 night sample is too wide to fit in the y-axis NDR range (the top value is 195).Figure 12Results of the NDR resampling on the two data sets of Cervières and Astrièves Observatory. The horizontal axis is the number of nights considered into the night sample and the vertical axis provides the NDR indicator obtained for each sampling set through a 1000-iteration bootstrapping with replacement calculation.Full size imageDiscussion on the required number of nightsWe can see in Fig. 12 that the NDR indicator and the confidence interval remain stable down to 200 nights. Below this threshold, the NDR starts to become unstable with growing confidence intervals. Based on this data, we can estimate that the minimum number of nights required to compute a robust NDR indicator is 200 (therefore between 7 and 8 months since there are periods around the full moons where there is no night portions recorded).However, depending on the measuring session objectives, the NDR indicator can be considered as accurate enough even when using a smaller number of nights. If the goal is simply to get a first estimate of the light pollution level at a given site, we can consider that 90 nights (a little more than 3 months of measures) are enough. On the opposite, if we want to perform a comparison between several sites for evaluating the impact of light pollution on a particular species, we might want to perform at least 200 nights of measurement to get a better accuracy for the NDR indicator. The experience from DarkSkyLab through many NSB measuring sessions is that 3 to 4 months of measures are required to get a meaningful density histogram, hence an accurate enough NDR indicator, so that a site can be sufficiently characterized from a light pollution perspective. Such a measuring period usually guarantees that the clear sky nominal NSB is well defined and that various cloud conditions have been observed. This estimate is sustained by the results obtained in Fig. 12.Value of the NDR indicator for ecological researchThe study of the impact of light pollution on biodiversity is currently in full expansion, amplifying a political and citizen demand for the reclamation of the night2,32,33.We identify three main contributions of the NDR indicator for ecological research. First, it overcomes the limits of an old problem of communication in terms of measurement units between disciplines and potentially limits the use of units without real meaning from a biodiversity point of view34,35. Secondly, the use of the NDR indicator limits the common biases linked to a characterization of the effects of anthropogenic light which is too limited in time and space35. Indeed, the life history traits of species are not only shaped by the intensity of light emitted into the nocturnal environment but also by its variation over time34,35,36. Currently, the characterization of light pollution is too often limited in time and space, which can lead to misinterpretation37. Thirdly, the NDR indicator provides ecological researchers with a unit of measurement that integrates a sufficiently long time step to study the impact of light pollution on the evolutionary processes at work in the life of species and particularly on population dynamics and animal behavior36,38,39.Limitations and future improvements of the NDR indicatorThe main limitation of the NDR indicator resides in the possible difficulty to identify a well defined value for the nominal NSB, i.e. the NSB value that represents the most common clear sky conditions of a given site. For the most part, this is due to the contribution of the Galactic plane to the zenithal sky brightness and, to a lower extent, to the contribution of other natural light sources (dense star fields, airglow and zodiacal light). The residual spread of NSB measures is due to changing atmospheric conditions at various time scales, but, for this particular contribution, we can expect a statistical compensation to eliminate a systematic associated bias.At the moment, the contribution of the Galactic plane and star fields is canceled into the NSB measures by calculating in the Noxi software the integrated flux of all the stars that belong to the field of view (using the UCAC4 star catalogue). However, this approach has proven some limitations, especially in the southern hemisphere where the Galactic center goes through the zenith and is particularly bright. A probable explanation for that lack of predictability is the fact that the Galactic plane contains diffuse sources such as nebulae which are not accounted for into the star catalogues and which actually cannot be ignored. To address this issue, DarkSkyLab has the project to create a brightness map of the Galactic plane with a square degree resolution or better so that the contribution of all sources can be correctly accounted for.In addition to better correcting the Galactic plane bias, an improvement must be made with regards to the NSB measures that need to be corrected. At the moment, all NSB measures are corrected from the Galactic plane bias without regards to the presence of clouds or high levels of light pollution. So a first heuristic must be implemented to only apply the bias correction to clear sky NSB measures. An other heuristic must also be developed to reduce the correction applied as a function of the NSB level.A third limitation of the NDR indicator is related to a possible lack of cloudy conditions at some sites (e.g. in the Atacama desert in Chile with more than 320 clear nights per year), the reason simply being that the NDR indicator requires the presence of clouds to differentiate the bright and dark extensions into the NSB density histograms. This means that the NDR indicator can hardly be used for such astronomy-oriented sites which experience rare cloudy conditions. More

  • in

    Fitness costs associated with a GABA receptor mutation conferring dieldrin resistance in Aedes albopictus

    Agnew P, Berticat C, Bedhomme S, Sidobre C, Michalakis Y (2004) Parasitism increases and decreases the costs of insecticide resistance in mosquitoes. Evolution 58:579–586CAS 
    PubMed 
    Article 

    Google Scholar 
    Ahmad NA, Endersby-Harshman NM, Mohd Mazni NR, Mohd Zabari NZA, Amran SNS, Ridhuan Ghazali MK et al. (2020) Characterization of sodium channel mutations in the Dengue vector mosquitoes Aedes aegypti and Aedes albopictus within the context of ongoing Wolbachia releases in Kuala Lumpur, Malaysia. Insects 11:529PubMed Central 
    Article 

    Google Scholar 
    Alout H, Ndam NT, Sandeu MM, Djégbe I, Chandre F, Dabiré RK et al. (2013) Insecticide resistance alleles affect vector competence of Anopheles gambiae s.s. for Plasmodium falciparum field isolates. PLoS ONE 8:e63849CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Andreasen MH, ffrench-Constant RH (2002) In situ hybridization to the Rdl locus on polytene chromosome 3L of Anopheles stephensi. Med Vet Entomol 16:452–455CAS 
    PubMed 
    Article 

    Google Scholar 
    Assogba BS, Djogbénou LS, Milesi P, Berthomieu A, Perez J, Ayala D et al. (2015) An ace-1 gene duplication resorbs the fitness cost associated with resistance in Anopheles gambiae, the main malaria mosquito. Sci Rep. 5:14529CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Assogba BS, Milesi P, Djogbénou LS, Berthomieu A, Makoundou P, Baba-Moussa LS et al. (2016) The ace-1 locus is amplified in all resistant Anopheles gambiae mosquitoes: fitness consequences of homogeneous and heterogeneous duplications. PloS Biol 14:e2000618PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Atyame CM, Alout H, Mousson L, Vazeille M, Diallo M, Weill M et al. (2019) Insecticide resistance genes affect Culex quinquefasciatus vector competence for West Nile virus. Proc Biol Sci 286:20182273CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Auteri M, La Russa F, Blanda V, Torina A (2018) Insecticide resistance associated with kdr mutations in Aedes albopictus: an update on worldwide evidences. Biomed Res Int 2018:e3098575Article 

    Google Scholar 
    Berticat C, Boquien G, Raymond M, Chevillon C (2002) Insecticide resistance genes induce a mating competition cost in Culex pipiens mosquitoes. Genet Res 79:41–47Berticat C, Duron O, Heyse D, Raymond M (2004) Insecticide resistance genes confer a predation cost on mosquitoes, Culex pipiens. Genet Res 83:189–196CAS 
    PubMed 
    Article 

    Google Scholar 
    Bhatia SC, Deobhankar RB (1963) Reversion of dieldrin-resistance in the field population of A. culicifacies in Maharashtra State (erstwhile Bombay State), India. Indian J Malariol 17:339–351CAS 
    PubMed 

    Google Scholar 
    Bonizzoni M, Gasperi G, Chen X, James AA (2013) The invasive mosquito species Aedes albopictus: current knowledge and future perspectives. Trends Parasitol 29:460–468PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Bourguet D, Guillemaud T, Chevillon C, Raymond M (2004) Fitness costs of insecticide resistance in natural breeding sites of the mosquito Culex pipiens. Evolution 58:128–135PubMed 
    Article 

    Google Scholar 
    Brooke BD, Hunt RH, Coetzee M (2000) Resistance to dieldrin + fipronil assorts with chromosome inversion 2La in the malaria vector Anopheles gambiae. Med Vet Entomol 14:190–194CAS 
    PubMed 
    Article 

    Google Scholar 
    Buckingham SD, Biggin PC, Sattelle BM, Brown LA, Sattelle DB (2005) Insect GABA receptors: splicing, editing, and targeting by antiparasitics and insecticides. Mol Pharm 68:942–951CAS 
    Article 

    Google Scholar 
    Chen H, Li K, Wang X, Yang X, Lin Y, Cai F et al. (2016) First identification of kdr allele F1534S in VGSC gene and its association with resistance to pyrethroid insecticides in Aedes albopictus populations from Haikou City, Hainan Island, China. Infect Dis Poverty 5:31PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Davari B, Vatandoost H, Oshaghi MA, Ladonni H, Enayati AA, Shaeghi M et al. (2007) Selection of Anopheles stephensi with DDT and dieldrin and cross-resistance spectrum to pyrethroids and fipronil. Pestic Biochem Physiol 89:97–103CAS 
    Article 

    Google Scholar 
    Delatte H, Paupy C, Dehecq JS, Thiria J, Failloux AB, Fontenille D (2008) Aedes albopictus, vector of Chikungunya and Dengue viruses in Reunion Island: biology and control. Parasite 15:3–13CAS 
    PubMed 
    Article 

    Google Scholar 
    Deng J, Guo Y, Su X, Liu S, Yang W, Wu Y et al. (2021) Impact of deltamethrin-resistance in Aedes albopictus on its fitness cost and vector competence. PLoS Negl Trop Dis 15:e0009391CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Djogbénou L, Weill M, Hougard J-M, Raymond M, Akogbéto M, Chandre F (2007) Characterization of insensitive acetylcholinesterase (ace-1R) in Anopheles gambiae (Diptera: Culicidae): resistance levels and dominance. J Med Entomol 44:805–810PubMed 

    Google Scholar 
    Du W, Awolola TS, Howell P, Koekemoer LL, Brooke BD, Benedict MQ et al. (2005) Independent mutations in the Rdl locus confer dieldrin resistance to Anopheles gambiae and An. arabiensis. Insect Mol Biol 14:179–183CAS 
    PubMed 
    Article 

    Google Scholar 
    Duron O, Labbé P, Berticat C, Rousset F, Guillot S, Raymond M et al. (2006) High Wolbachia density correlates with cost of infection for insecticide resistant Culex pipiens mosquitoes. Evolution 60:303–314CAS 
    PubMed 
    Article 

    Google Scholar 
    ffrench-Constant RH, Rocheleau TA, Steichen JC, Chalmers AE (1993) A point mutation in a Drosophila GABA receptor confers insecticide resistance. Nature 363:449–451CAS 
    PubMed 
    Article 

    Google Scholar 
    ffrench-Constant RH, Anthony N, Aronstein K, Rocheleau T, Stilwell G (2000) Cyclodiene insecticide resistance: from molecular to population genetics. Annu Rev Entomol 45:449–466CAS 
    PubMed 
    Article 

    Google Scholar 
    Fox J, Weisberg S (2019) An R companion to applied regression, 3rd edn. SAGE, Thousand Oaks California, https://socialsciences.mcmaster.ca/jfox/Books/Companion/
    Google Scholar 
    Freeman JC, Smith LB, Silva JJ, Fan Y, Sun H, Scott JG (2021) Fitness studies of insecticide resistant strains: lessons learned and future directions. Pest Manag Sci 77:3847–3856CAS 
    PubMed 
    Article 

    Google Scholar 
    Gratz NG (2004) Critical review of the vector status of Aedes albopictus. Med Vet Entomol 18:215–227CAS 
    PubMed 
    Article 

    Google Scholar 
    Grau-Bové X, Tomlinson S, O’Reilly AO, Harding NJ, Miles A, Kwiatkowski D et al. (2020) Evolution of the insecticide target Rdl in African Anopheles is driven by interspecific and interkaryotypic introgression. Mol Biol Evol 37:2900–2917PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Grigoraki L, Lagnel J, Kioulos I, Kampouraki A, Morou E, Labbé P et al. (2015) Transcriptome profiling and genetic study reveal amplified carboxylesterase genes implicated in temephos resistance, in the Asian tiger mosquito Aedes albopictus. PLoS Negl Trop Dis 9:e0003771PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Hamon J, Garret-Jones C (1962) Insecticide-resistance in major vectors of malaria, and its operational importance. Bull World Health Organ, Geneva
    Google Scholar 
    Hartley CJ, Newcomb RD, Russell RJ, Yong CG, Stevens JR, Yeates DK et al. (2006) Amplification of DNA from preserved specimens shows blowflies were preadapted for the rapid evolution of insecticide resistance. Proc Natl Acad Sci USA 103:8757–8762CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Hemingway J, Ranson H (2000) Insecticide resistance in insect vectors of human disease. Annu Rev Entomol 45:371–391CAS 
    PubMed 
    Article 

    Google Scholar 
    Hemingway J, Hawkes NJ, McCarroll L, Ranson H (2004) The molecular basis of insecticide resistance in mosquitoes. Insect Biochem Mol Biol 34:653–665CAS 
    PubMed 
    Article 

    Google Scholar 
    Hosie AM, Baylis HA, Buckingham SD, Sattelle DB (1995) Actions of the insecticide fipronil, on dieldrin-sensitive and -resistant GABA receptors of Drosophila melanogaster. Br J Pharm 115:909–912CAS 
    Article 

    Google Scholar 
    Ishak IH, Riveron JM, Ibrahim SS, Stott R, Longbottom J, Irving H et al. (2016) The Cytochrome P450 gene CYP6P12 confers pyrethroid resistance in kdr-free Malaysian populations of the Dengue vector Aedes albopictus. Sci Rep. 6:24707CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Kasai S, Ng LC, Lam-Phua SG, Tang CS, Itokawa K, Komagata O et al. (2011) First detection of a putative knockdown resistance gene in major mosquito vector, Aedes albopictus. Jpn J Infect Dis 64:217–221CAS 
    PubMed 
    Article 

    Google Scholar 
    Kliot A, Ghanim M (2012) Fitness costs associated with insecticide resistance. Pest Manag Sci 68:1431–1437CAS 
    PubMed 
    Article 

    Google Scholar 
    Kolaczinski J, Curtis C (2001) Laboratory evaluation of fipronil, a phenylpyrazole insecticide, against adult Anopheles (Diptera: Culicidae) and investigation of its possible cross-resistance with dieldrin in Anopheles stephensi. Pest Manag Sci 57:41–45CAS 
    PubMed 
    Article 

    Google Scholar 
    Kraemer MU, Sinka ME, Duda KA, Mylne AQ, Shearer FM, Barker CM et al. (2015) The global distribution of the arbovirus vectors Aedes aegypti and Ae. albopictus. Elife 4:e08347PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Labbé P, David J-P, Alout H, Milesi P, Djogbénou L, Pasteur N et al. (2017) 14 – Evolution of resistance to insecticide in disease vectors. In: Tibayrenc M (ed) Genetics and Evolution of Infectious Diseases, Second Edition. Elsevier, London, p 313–339Chapter 

    Google Scholar 
    Latreille AC, Milesi P, Magalon H, Mavingui P, Atyame CM (2019) High genetic diversity but no geographical structure of Aedes albopictus populations in Réunion Island. Parasit Vectors 12:597PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Lebon C, Alout H, Zafihita S, Dehecq JS, Weill M, Tortosa P et al. (2022) Spatio-temporal dynamics of a dieldrin resistance gene in Aedes albopictus and Culex quinquefasciatus populations from Reunion Island. J Insect Sci 22:4PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Lebon C, Soupapoule K, Wilkinson DA, Goff GL, Damiens D, Gouagna LC (2018) Laboratory evaluation of the effects of sterilizing doses of γ-rays from Caesium-137 source on the daily flight activity and flight performance of Aedes albopictus males. PLoS ONE 13:e0202236PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Li Y, Xu J, Zhong D, Zhang H, Yang W, Zhou G et al. (2018) Evidence for multiple-insecticide resistance in urban Aedes albopictus populations in southern China. Parasit Vectors 11:4PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Low VL, Vinnie-Siow WY, Lim YAL, Tan TK, Leong CS, Chen CD et al. (2015) First molecular genotyping of A302S mutation in the gamma aminobutyric acid (GABA) receptor in Aedes albopictus from Malaysia. Trop Biomed 32:554–556CAS 
    PubMed 

    Google Scholar 
    McKenzie BA, Wilson AE, Zohdy S (2019) Aedes albopictus is a competent vector of Zika virus: a meta-analysis. PLoS ONE 14:e0216794CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Milesi P, Pocquet N, Labbé P (2013) BioRssay: A R script for bioassay analyses. http://www.isem.univ-montp2.fr/recherche/equipes/genomique-de-ladaptation/personnel/labbepierrick/Moyes CL, Vontas J, Martins AJ, Ng LC, Koou SY, Dusfour I et al. (2017) Contemporary status of insecticide resistance in the major Aedes vectors of arboviruses infecting humans. PLoS Negl Trop Dis 11:e0005625PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Ozoe Y, Kita T, Ozoe F, Nakao T, Sato K, Hirase K (2013) Insecticidal 3-benzamido-N-phenylbenzamides specifically bind with high affinity to a novel allosteric site in housefly GABA receptors. Pestic Biochem Physiol 107:285–292CAS 
    PubMed 
    Article 

    Google Scholar 
    Paupy C, Ollomo B, Kamgang B, Moutailler S, Rousset D, Demanou M et al. (2009) Comparative role of Aedes albopictus and Aedes aegypti in the emergence of Dengue and Chikungunya in central Africa. Vector Borne Zoonotic Dis 10:259–266Article 

    Google Scholar 
    Peakall R, Smouse PE (2012) GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research—an update. Bioinformatics 28:2537–2539CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Platt N, Kwiatkowska RM, Irving H, Diabaté A, Dabire R, Wondji CS (2015) Target-site resistance mutations (kdr and RDL), but not metabolic resistance, negatively impact male mating competiveness in the malaria vector Anopheles gambiae. Heredity 115:243–252CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    R Core Team (2019) R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria, https://www.R-project.org/
    Google Scholar 
    Ranson H, Burhani J, Lumjuan N, Black WCI (2010) Insecticide resistance in Dengue vectors. TropIKA.net [online] 1. http://journal.tropika.net/scielo.php?script=sci_arttext&pid=S2078-86062010000100003&lng=en&nrm=iso. Accessed 03 March 2022Raymond M, Berticat C, Weill M, Pasteur N, Chevillon C (2001) Insecticide resistance in the mosquito Culex pipiens: what have we learned about adaptation? Genetica 112–113:287–296PubMed 
    Article 

    Google Scholar 
    Renault P, Solet J-L, Sissoko D, Balleydier E, Larrieu S, Filleul L et al. (2007) A major epidemic of Chikungunya virus infection on Réunion Island, France, 2005–2006. Am J Trop Med Hy 77:727–731Article 

    Google Scholar 
    Rowland M (1991a) Behaviour and fitness of γHCH/dieldrin resistant and susceptible female Anopheles gambiae and An. stephensi mosquitoes in the absence of insecticide. Med Vet Entomol 5:193–206CAS 
    PubMed 
    Article 

    Google Scholar 
    Rowland M (1991b) Activity and mating competitiveness of γHCH/dieldrin resistant and susceptible male and virgin female Anopheles gambiae and An. stephensi mosquitoes, with assessment of an insecticide-rotation strategy. Med Vet Entomol 5:207–222CAS 
    PubMed 
    Article 

    Google Scholar 
    Russell VL (2021) Emmeans: Estimated Marginal Means, aka Least-Squares Means. R package version 1.5.5.1. https://CRAN.R-project.org/package=emmeansSu X, Guo Y, Deng J, Xu J, Zhou G, Zhou T et al. (2019) Fast emerging insecticide resistance in Aedes albopictus in Guangzhou, China: alarm to the Dengue epidemic. PLoS Negl Trop Dis 13:e0007665CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Tantely ML, Tortosa P, Alout H, Berticat C, Berthomieu A, Rutee A et al. (2010) Insecticide resistance in Culex pipiens quinquefasciatus and Aedes albopictus mosquitoes from La Réunion Island. Insect Biochem Mol Biol 40:317–324CAS 
    PubMed 
    Article 

    Google Scholar 
    Taskin BG, Dogaroglu T, Kilic S, Dogac E, Taskin V (2016) Seasonal dynamics of insecticide resistance, multiple resistance, and morphometric variation in field populations of Culex pipiens. Pestic Biochem Physiol 129:14–27CAS 
    PubMed 
    Article 

    Google Scholar 
    Taylor‐Wells J, Brooke BD, Bermudez I, Jones AK (2015) The neonicotinoid imidacloprid, and the pyrethroid deltamethrin, are antagonists of the insect Rdl GABA receptor. J Neurochem 135:705–713PubMed 
    Article 

    Google Scholar 
    Therneau T (2015) A Package for Survival Analysis in S. R package version 2.38. https://CRAN.R-project.org/package=survivalThompson M, Shotkoski F, ffrench-Constant R (1993) Cloning and sequencing of the cylodienne insecticide resistance from the yellow fewer Aedes aegypti. FEBS Lett 325:187–190CAS 
    PubMed 
    Article 

    Google Scholar 
    Tsetsarkin KA, Vanlandingham DL, McGee CE, Higgs S (2007) A single mutation in Chikungunya virus affects vector specificity and epidemic potential. PLoS Pathog 3:e201PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Vontas J, Kioulos E, Pavlidi N, Morou E, della Torre A, Ranson H (2012) Insecticide resistance in the major Dengue vectors Aedes albopictus and Aedes aegypti. Pestic Biochem Physiol 104:126–131CAS 
    Article 

    Google Scholar 
    Wondji CS, Dabire RK, Tukur Z, Irving H, Djouaka R, Morgan JC (2011) Identification and distribution of a GABA receptor mutation conferring dieldrin resistance in the malaria vector Anopheles funestus in Africa. Insect Biochem Mol Biol 41:484–491CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Xu J, Bonizzoni M, Zhong D, Zhou G, Cai S, Li Y et al. (2016) Multi-country survey revealed prevalent and novel F1534S mutation in voltage-gated sodium channel (VGSC) gene in Aedes albopictus. PLoS Negl Trop Dis 10:e0004696PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Yang C, Huang Z, Li M, Feng X, Qiu X (2017) RDL mutations predict multiple insecticide resistance in Anopheles sinensis in Guangxi, China. Malar J 16:482PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Zhou X, Yang C, Liu N, Li M, Tong Y, Zeng X et al. (2019) Knockdown resistance (kdr) mutations within seventeen field populations of Aedes albopictus from Beijing China: first report of a novel V1016G mutation and evolutionary origins of kdr haplotypes. Parasit Vectors 12:180PubMed 
    PubMed Central 
    Article 

    Google Scholar  More

  • in

    Author Correction: Widespread extinction debts and colonization credits in United States breeding bird communities

    In the version of this article initially published, there were errors in equations and notations in the Methods “Model development” subsection which arose during manuscript preparation; the errors affect presentation of the study but not the analysis, results, or code provided with the article. Clarifications to text and equations follow.In Equation (1), “N” replaces “Normal”; in Equations (2), (3), (7) and in text directly below Equations (3), (5) and (7), “ys,i,z” now replaces “Δxs,t1, t2.” In the two paragraphs below Equation (2), “t2 = 2016” and “t1 = 2001” now replace “2016” and “2001” in five instances. Further, Equations (5)–(7) have been revised as follows:$$begin{array}{ll}fleft( {x_{s,t}} right) = {{{mathrm{exp}}}} & left( {beta _0 + mathop {sum }limits_{i = 1}^{I = 5} beta _{1,i} x_{s,i,t} + mathop {sum }limits_{i = 1}^{I = 5} mathop {sum }limits_{k = i}^{K = 5} beta _{2,i,k}x_{s,i,t}x_{k,s,t}}right. \ & quad quad left. {+ mathop {sum }limits_{i = 1}^{I = 5} mathop {sum }limits_{k = 1, k neq i}^{K = 5} beta _{3,i,k}x_{s,i,t}x_{k,s,t}} right)end{array} {rm{Revised}} {rm{Eq}}. (5)$$$$begin{array}{ll}fleft( {x_{s,t}} right) \ = expleft( {beta _0 + mathop {sum }limits_{i = 1}^{I = 5} mathop {sum }limits_{j = 1}^{J = 2} beta _{0,i,j,}x_{i,s,t}^j + mathop {sum }limits_{i = 1}^{I = 5} mathop {sum }limits_{k = i + 1}^{K = 6} beta _{1,i,k}x_{i,s,t}x_{k,s,t}} right) {mathrm{Original}} {rm{Eq}}. (5)end{array}$$$$y_{s,i,z} = left{ {begin{array}{*{20}{l}} {y_{s,i,1} = left| {Delta x_{s,i}} right|,} hfill & {y_{s,i,2} = 0,} hfill & {{{{mathrm{if}}}},Delta x_{s,i} < 0} hfill \ {y_{s,i,1} = 0,} hfill & {y_{s,i,2} = Delta x_{s,i}} hfill & {{{{mathrm{otherwise}}}}} hfill end{array}} right. {rm{Revised}} {rm{Eq}}. (6)$$$$x_{i,s,} = left{ {begin{array}{*{20}{l}} {x_{1,i,s} = left| {Delta x_{i,s}} right|,} hfill & {x_{2,i,s} = 0,} hfill & {if,Delta x_{i,s} < 0} hfill \ {x_{1,i,s} = 0,} hfill & {x_{2,i,s} = Delta x_{i,s},} hfill & {otherwise} hfill end{array}} right. {rm{Original}} {rm{Eq}}. (6)$$$$omega left( {y_{s,i,z};gamma } right) = {{{mathrm{exp}}}}left( {mathop {sum }limits_{i = 1}^{I = 5} mathop {sum }limits_{z = 1}^{Z = 2} - gamma _{i,z} y_{s,i,z}} right) {rm{Revised}} {rm{Eq}}. (7)$$$$omega left( {Delta x_{s,t_1,t_2};gamma } right) = expleft( {mathop {sum }limits_{i = 1}^{I = 5} - gamma _{i,z}Delta x_{z,s,i}} right) {rm{Original}} {rm{Eq}}. (7)$$All changes have been made in the HTML and PDF versions of the article. More

  • in

    Resolving malaria’s dry-season dilemma

    Seasonal fluctuations in animal population dynamics are among the most fundamental attributes of life on Earth. A long recognized but poorly understood example is the dramatic seasonal fluctuation in the abundance of malaria vectors in the semi-arid savannah and Sahel regions of Africa. In these regions, the vector mosquitoes largely disappear during a prolonged 3- to 8-month dry season, when lack of rain causes the aquatic larval habitats to disappear. As a result, malaria transmission plummets. When the rains return, the mosquito vectors rapidly reappear, leading to a resurgence of malaria transmission. How the vector populations are able to persist through the prolonged dry season and rapidly rebound with the onset of rains is referred to as the ‘dry-season malaria paradox’, and has remained an enduring mystery of malariology for nearly 100 years. Writing in Nature Ecology & Evolution, Faiman et al.1 help to resolve this mystery by using an innovative isotopic labelling strategy: they demonstrate that at least approximately 20% of the local population of the malaria vector Anopheles coluzzi in the West African Sahel survive the dry season locally by undergoing summer dormancy, known as aestivation. More

  • in

    Isotopic evidence that aestivation allows malaria mosquitoes to persist through the dry season in the Sahel

    Adamou, A. et al. The contribution of aestivating mosquitoes to the persistence of Anopheles gambiae in the Sahel. Malar. J. 10, 151 (2011).PubMed 
    PubMed Central 

    Google Scholar 
    Huestis, D. L. et al. Seasonal variation in metabolic rate, flight activity and body size of Anopheles gambiae in the Sahel. J. Exp. Biol. 215, 2013–2021 (2012).PubMed 
    PubMed Central 

    Google Scholar 
    Huestis, D. L. et al. Variation in metabolic rate of Anopheles gambiae and A. arabiensis in a Sahelian village. J. Exp. Biol. 214, 2345–2353 (2011).PubMed 
    PubMed Central 

    Google Scholar 
    Lehmann, T. et al. Aestivation of the African malaria mosquito, Anopheles gambiae in the Sahel. Am. J. Trop. Med. Hyg. 83, 601–606 (2010).PubMed 
    PubMed Central 

    Google Scholar 
    Yaro, A. S. et al. Dry season reproductive depression of Anopheles gambiae in the Sahel. J. Insect Physiol. 58, 1050–1059 (2012).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Omer, S. M. & Cloudsley-Thompson, J. L. Survival of female Anopheles gambiae Giles through a 9-month dry season in Sudan. Bull. World Health Organ. 42, 319 (1970).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Omer, S. M. & Cloudsley-Thompson, J. L. Dry season biology of Anopheles gambiae Giles in the Sudan. Nature 217, 879–880 (1968).
    Google Scholar 
    Holstein, M. H. Biology of Anopheles gambiae (1954). World Health Organization.Andrade, C. M. et al. Increased circulation time of Plasmodium falciparum underlies persistent asymptomatic infection in the dry season. Nat. Med. 26, 1929–1940 (2020).CAS 
    PubMed 

    Google Scholar 
    Coulibaly, D. et al. Spatio-temporal dynamics of asymptomatic malaria: bridging the gap between annual malaria resurgences in a Sahelian environment. Am. J. Trop. Med. Hyg. 27, 1761–1769 (2017).
    Google Scholar 
    Gillies, M. & Wilkes, T. A study of the age-composition of populations of Anopheles gambiae Giles and A. funestus Giles in north-eastern Tanzania. Bull. Entomol. Res. 56, 237–262 (1965).CAS 
    PubMed 

    Google Scholar 
    Gillies, M. T. & De Meillon, B. The Anophelinae of Africa south of the Sahara (Ethiopian Zoogeographical Region) (Johannesburg: South African Institute for Medical Research, 1968).Dao, A. et al. Signatures of aestivation and migration in Sahelian malaria mosquito populations. Nature 516, 387–390 (2014).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Thomson, J. G. Malaria in Nyasaland. Proc. R. Soc. Med. 28, 391–404 (1934).
    Google Scholar 
    Huestis, D. L. et al. Windborne long-distance migration of malaria mosquitoes in the Sahel. Nature 574, 404–408 (2019).CAS 
    PubMed 

    Google Scholar 
    Lambert, B., North, A., Burt, A. & Godfray, H. C. J. The use of driving endonuclease genes to suppress mosquito vectors of malaria in temporally variable environments. Malar. J. 17, 154 (2018).PubMed 
    PubMed Central 

    Google Scholar 
    Verhulst, N. O., Loonen, J. A. C. M. & Takken, W. Advances in methods for colour marking of mosquitoes. Parasit. Vectors 6, 200 (2013).PubMed 
    PubMed Central 

    Google Scholar 
    Hagler, J. R. & Jackson, C. G. Methods for marking insects: current techniques and future prospects. Annu. Rev. Entomol. 46, 511–543 (2001).CAS 
    PubMed 

    Google Scholar 
    Hamer, G. L. et al. Dispersal of adult culex mosquitoes in an urban West Nile virus hotspot: a mark–capture study incorporating stable isotope enrichment of natural larval habitats. PLoS Negl. Trop. Dis. 8, e2768 (2014).PubMed 
    PubMed Central 

    Google Scholar 
    Hamer, G. L. et al. Evaluation of a stable isotope method to mark naturally-breeding larval mosquitoes for adult dispersal studies. J. Med. Entomol. 49, 61–70 (2012).CAS 
    PubMed 

    Google Scholar 
    Opiyo, M. A. et al. Using stable isotopes of carbon and nitrogen to mark wild populations of Anopheles and Aedes mosquitoes in south-eastern Tanzania. PLoS ONE 11, e0159067 (2016).PubMed 
    PubMed Central 

    Google Scholar 
    Hood-Nowotny, R., Mayr, L. & Knols, B. Use of carbon-13 as a population marker for Anopheles arabiensis in a sterile insect technique (SIT) context. Malar. J. 5, 6 (2006).PubMed 
    PubMed Central 

    Google Scholar 
    Hood-Nowotny, R. & Knols, B. G. J. Stable isotope methods in biological and ecological studies of arthropods. Entomol. Exp. Appl. 124, 3–16 (2007).CAS 

    Google Scholar 
    Hood-Nowotny, R. et al. Intrinsic and synthetic stable isotope marking of tsetse flies. J. Insect Sci. 11, 79 (2011).PubMed 
    PubMed Central 

    Google Scholar 
    Atzrodt, J., Derdau, V., Kerr, W. J. & Reid, M. Deuterium- and tritium-labelled compounds: applications in the life sciences. Angew. Chem. Int. Ed. 57, 1758–1784 (2018).CAS 

    Google Scholar 
    Copia, L., Wassenaar, L. I., Terzer-Wassmuth, S., Belachew, D. L. & Araguas-Araguas, L. J. Comparative evaluation of 2H- versus 3H-based enrichment factor determination on the uncertainty and accuracy of low-level tritium analyses of environmental waters. Appl. Radiat. Isot. 176, 109850 (2021).CAS 
    PubMed 

    Google Scholar 
    Begon, M., Harper, J. & Townsend, C. Ecology: Individuals, Populations and Communities (Blackwell Science, 1996).Faiman, R. et al. Marking mosquitoes in their natural larval sites using 2H-enriched water: a promising approach for tracking over extended temporal and spatial scales. Methods Ecol. Evol. 10, 1274–1285 (2019).PubMed 
    PubMed Central 

    Google Scholar 
    Florkin, M. Chemical Zoology: Arthropoda Part B (Academic Press, 2014).Hackman, R. H. & Goldberg, M. Studies on chitin VI. The nature of alpha-and beta-chitins. Aust. J. Biol. Sci. 18, 935–946 (1965).CAS 
    PubMed 

    Google Scholar 
    Faiman, R. et al. Quantifying flight aptitude variation in wild Anopheles gambiae in order to identify long-distance migrants. Malar. J. 19, 263 (2020).PubMed 
    PubMed Central 

    Google Scholar 
    Huestis, D. L. & Lehmann, T. Ecophysiology of Anopheles gambiae s.l.: persistence in the Sahel. Infect. Genet. Evol. 28, 648–661 (2014).PubMed 
    PubMed Central 

    Google Scholar 
    Lehmann, T. et al. Seasonal variation in spatial distributions of Anopheles gambiae in a Sahelian village: evidence for aestivation. J. Med. Entomol. 51, 27–38 (2014).PubMed 

    Google Scholar 
    Costantini, C. et al. Density, survival and dispersal of Anopheles gambiae complex mosquitoes in a West African Sudan savanna village. Med. Vet. Entomol. 10, 203–219 (1996).CAS 
    PubMed 

    Google Scholar 
    Toure, Y. T. et al. Mark–release–recapture experiments with Anopheles gambiae s.l. in Banambani Village, Mali, to determine population size and structure. Med. Vet. Entomol. 12, 74–83 (1998).CAS 
    PubMed 

    Google Scholar 
    Faiman, R. et al. A novel fluorescence and DNA combination for versatile, long-term marking of mosquitoes. Methods Ecol. Evol. https://doi.org/10.1111/2041-210X.13592 (2021).Brattström, O., Bensch, S., Wassenaar, L. I., Hobson, K. A. & Åkesson, S. Understanding the migration ecology of European red admirals Vanessa atalanta using stable hydrogen isotopes. Ecography 33, 720–729 (2010).
    Google Scholar 
    Hobson, K. A., Jinguji, H., Ichikawa, Y., Kusack, J. W. & Anderson, R. C. Long-distance migration of the globe skimmer dragonfly to Japan revealed using stable hydrogen (δ 2H) isotopes. Environ. Entomol. 50, 247–255 (2020).
    Google Scholar 
    Schilling, E. G. et al. Phenological and isotopic evidence for migration as a life history strategy in Aeshna canadensis (family: Aeshnidae) dragonflies. Ecol. Entomol. 46, 209–219 (2021).
    Google Scholar 
    Girard, P., Hillaire-Marcel, C. & Oga, M. S. Determining the recharge mode of Sahelian aquifers using water isotopes. J. Hydrol. 197, 189–202 (1997).CAS 

    Google Scholar 
    Gutiérrez-Expósito, C., Ramírez, F., Afán, I., Forero, M. & Hobson, K. A. Toward a deuterium feather isoscape for sub-Saharan Africa: progress, challenges and the path ahead. PLoS ONE https://doi.org/10.1371/journal.pone.0135938 (2015).Lutz, A., Thomas, J. M. & Panorska, A. Environmental controls on stable isotope precipitation values over Mali and Niger, West Africa. Environ. Earth Sci. 62, 1749–1759 (2011).CAS 

    Google Scholar 
    Risi, C. et al. Understanding the Sahelian water budget through the isotopic composition of water vapor and precipitation. J. Geophys. Res. Atmos. 115, 1–23 (2010).
    Google Scholar 
    Tremoy, G. et al. A 1-year long δ18O record of water vapor in Niamey (Niger) reveals insightful atmospheric processes at different timescales. Geophys. Res. Lett. 39, 1–5 (2012).
    Google Scholar 
    Terzer‐Wassmuth, S., Wassenaar, L. I., Welker, J. M., Araguás-Araguás, L. J. Improved high‐resolution global and regionalized isoscapes of δ18O, δ2H and d‐excess in precipitation. Hydrol. Process. 35 (2021).Hobson, K. A. et al. A multi-isotope (δ13C, δ15N, δ2H) feather isoscape to assign Afrotropical migrant birds to origins. Ecosphere 3, art44 (2012).
    Google Scholar 
    Diuk-Wasser, M. A. et al. Effect of rice cultivation patterns on malaria vector abundance in rice-growing villages in Mali. Am. J. Trop. Med. Hyg. 76, 869–874 (2007).PubMed 

    Google Scholar 
    Sogoba, N. et al. Malaria transmission dynamics in Niono, Mali: the effect of the irrigation systems. Acta Trop. 101, 232–240 (2007).PubMed 

    Google Scholar 
    Florio, J. et al. Diversity, dynamics, direction, and magnitude of high-altitude migrating insects in the Sahel. Sci. Rep. 10, 20523 (2020).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Wilkins, E. E., Howell, P. I. & Benedict, M. Q. IMP PCR primers detect single nucleotide polymorphisms for Anopheles gambiae species identification, Mopti and Savanna rDNA types, and resistance to dieldrin in Anopheles arabiensis. Malar. J. 5, 125 (2006).PubMed 
    PubMed Central 

    Google Scholar 
    Wassenaar, L. I. & Hobson, K. A. Comparative equilibration and online technique for determination of non-exchangeable hydrogen of keratins for use in animal migration studies. Isotopes Environ. Health Stud. 39, 211–217 (2003).CAS 
    PubMed 

    Google Scholar 
    Chesson, L. A., Podlesak, D. W., Cerling, T. E. & Ehleringer, J. R. Evaluating uncertainty in the calculation of non-exchangeable hydrogen fractions within organic materials. Rapid Commun. Mass Spectrom. 23, 1275–1280 (2009).CAS 
    PubMed 

    Google Scholar 
    Schimmelmann, A. Determination of the concentration and stable isotopic composition of nonexchangeable hydrogen in organic matter. Anal. Chem. 63, 2456–2459 (1991).CAS 

    Google Scholar 
    Speakman, J. Doubly Labelled Water: Theory and Practice (Chapman & Hall, 1997).Base SAS 9.4 Procedures Guide (SAS Institute, 2015).Cade, B. S. & N, B. R. A gentle introduction to quantile regression for ecologists. Front. Ecol. Environ. 1, 412–420 (2003).
    Google Scholar 
    SAS/STAT® 15.1 User’s Guide (SAS Institute, 2018).Mcclintock, B. T. et al. Uncovering ecological state dynamics with hidden Markov models. Ecol. Lett. 23, 1878–1903 (2020).PubMed 
    PubMed Central 

    Google Scholar 
    Issam, M., Naulet, N., Martin, M. L. & Martin, G. J. A site-specific and multielement approach to the determination of liquid–vapor isotope fractionation parameters: the case of alcohols. J. Phys. Chem. 94, 8303–8309 (1990).
    Google Scholar 
    Linderstrøm-Lang, C. U. & Vaslow, F. Isotope effect on the vapor pressures of water–ethanol and deuterium oxide–ethanol-d mixtures. J. Phys. Chem. 72, 2645–2650 (1968).
    Google Scholar 
    Ventura, M. & Jeppesen, E. Effects of fixation on freshwater invertebrate carbon and nitrogen isotope composition and its arithmetic correction. Hydrobiologia 632, 297–308 (2009).CAS 

    Google Scholar  More

  • in

    Unique thermal sensitivity imposes a cold-water energetic barrier for vertical migrators

    Robison, B. H. Conservation of deep pelagic biodiversity. Conserv. Biol. 23, 847–858 (2009).
    Google Scholar 
    Fernandez-Alamo, M. A. & Färber-Lorda, J. Zooplankton and the oceanography of the eastern tropical Pacific: a review. Prog. Oceanogr. 69, 318–359 (2006).
    Google Scholar 
    Bianchi, D., Galbraith, E. D., Carozza, D. A., Mislan, K. A. S. & Stock, C. A. Intensification of open-ocean oxygen depletion by vertically migrating animals. Nat. Geosci. 6, 545–548 (2013).CAS 

    Google Scholar 
    Steinberg, D. K. & Landry, M. R. Zooplankton and the ocean carbon cycle. Annu. Rev. Mar. Sci. 9, 413–444 (2017).
    Google Scholar 
    Kiko, R. & Hauss, H. On the estimation of zooplankton-mediated active fluxes in oxygen minimum zones regions. Front. Mar. Sci. https://doi.org/10.3389/fmars.2019.00741 (2019).Longhurst, A., Bedo, A., Harrison, W., Head, E. & Sameoto, D. Vertical flux of respiratory carbon by oceanic diel migrant biota. Deep Sea Res. Part I 37, 685–694 (1990).CAS 

    Google Scholar 
    Elder, L. E. & Seibel, B. A. The thermal stress response to diel vertical migration in the hyperiid amphipod, Phronima sedentaria. Comp. Biochem. Physiol. A 187, 20–26 (2015).CAS 

    Google Scholar 
    Tremblay, N., Gomez-Gutierrez, J., Zenteno-Savin, T., Robinson, C. J. & Sanchez-Velascoa, L. Role of oxidative stress in seasonal and daily vertical migration of three krill species in the Gulf of California. Limnol. Oceanogr. 55, 2570–2584 (2010).CAS 

    Google Scholar 
    Lopes, A. R. et al. Oxidative stress in deep scattering layers: heat shock response and antioxidant enzymes activities of myctophid fishes thriving in oxygen minimum zones. Deep Sea Res. Part I 82, 10–16 (2013).CAS 

    Google Scholar 
    Seibel, B. A., Schneider, J., Kaartvedt, S., Wishner, K. F. & Daly, K. L. Hypoxia tolerance and metabolic suppression in oxygen minimum zone euphausiids: implications for ocean deoxygenation and biogeochemical cycles. Integr. Comp. Biol. https://doi.org/10.1093/icb/icw091 (2016).Seibel, B. A. et al. Metabolic suppression during protracted exposure to hypoxia in the jumbo squid, Dosidicus gigas, living in an oxygen minimum zone. J. Exp. Biol. 217, 2710–2716 (2014).
    Google Scholar 
    Wishner, K. F. et al. Ocean deoxygenation and zooplankton: very small oxygen differences matter. Sci. Adv. 4, eaau5180 (2018).CAS 

    Google Scholar 
    Koslow, J. A., Goericke, R., Lara-Lopez, A. & Watson, W. Impact of declining intermediate-water oxygen on deepwater fishes in the California Current. Mar. Ecol. Prog. Ser. 436, 207–218 (2011).
    Google Scholar 
    Oschlies, A. A committed fourfold increase in ocean oxygen loss. Nat. Commun. 12, 2307 (2021).CAS 

    Google Scholar 
    Wishner, K. F., Seibel, B. A. & Outram, D. Ocean deoxygenation and copepods: coping with oxygen minimum zone variability. Biogeosciences 17, 2315–2339 (2020).
    Google Scholar 
    Stramma, L. et al. Expansion of oxygen minimum zones may reduce available habitat for tropical pelagic fishes. Nat. Clim. Change 2, 33–37 (2012).CAS 

    Google Scholar 
    Köhn, E. E., Münnich, M., Vogt, M., Desmmet, F. & Gruber, N. Strong habitat compression by extreme shoaling events of hypoxic waters in the Eastern Pacific. J. Geophys. Res. Oceans 127, e2022JC018429 (2022).
    Google Scholar 
    Poloczanska, E. S. et al. Global imprint of climate change on marine life. Nat. Clim. Change 3, 919–925 (2013).
    Google Scholar 
    Pinsky, M. L., Selden, R. L. & Kitchel, Z. J. Climate-driven shifts in marine species ranges: scaling from organisms to communities. Annu. Rev. Mar. Sci. 12, 153–179 (2020).
    Google Scholar 
    Cavole, L. M. et al. Biological impacts of the 2013–2015 warm-water anomaly in the northeast Pacific: winners, losers, and the future. Oceanography 29, 273–285 (2016).
    Google Scholar 
    Lavaniegosa, B. E., Jiménez-Herrera, M. A. & Ambriz-Arreola, I. Unusually low euphausiid biomass during the warm years of 2014–2016 in the transition zone of the California Current. Deep Sea Res. Part II 1, 69–170 (2019).
    Google Scholar 
    Lilly, L. E. & Ohman, M. D. Euphausiid spatial displacements and habitat shifts in the southern California Current system in response to El Niño variability. Prog. Oceanogr. 193, 102544 (2021).
    Google Scholar 
    Zeidberg, L. D. & Robison, B. H. Invasive range expansion by the Humboldt squid, Dosidicus gigas, in the eastern North Pacific. Proc. Natl Acad. Sci. USA 104, 12948–12950 (2007).CAS 

    Google Scholar 
    Szesciorka, A. R. et al. Timing is everything: drivers of interannual variability in blue whale migration. Sci. Rep. 10, 7710 (2020).CAS 

    Google Scholar 
    Hoving, H.-J. et al. Extreme plasticity in life‐history strategy allows a migratory predator (jumbo squid) to cope with a changing climate. Glob. Change Biol. 19, 2089–2103 (2013).
    Google Scholar 
    Boscolo-Galazzo, F. et al. Temperature controls carbon cycling and biological evolution in the ocean twilight zone. Science 371, 1148–1152 (2021).CAS 

    Google Scholar 
    Deutsch, C., Ferrel, A., Seibel, B. A., Pörtner, H.-O. & Huey, R. B. Climate change tightens a metabolic constraint on marine habitats. Science 348, 1132–1135 (2015).CAS 

    Google Scholar 
    Seibel, B. A. & Deutsch, C. Oxygen supply capacity in animals evolves to meet maximum demand at the current oxygen partial pressure regardless of size or temperature. J. Exp. Biol. 223, jeb210492 (2020).
    Google Scholar 
    Deutsch, C., Penn, J. L. & Seibel, B. A. Diverse hypoxia and thermal tolerances shape biogeography of marine animals. Nature 585, 557–562 (2020).CAS 

    Google Scholar 
    Childress, J. J. Are there physiological and biochemical adaptations of metabolism in deep-sea animals? Trends Ecol. Evol. 10, 30–36 (1995).CAS 

    Google Scholar 
    Seibel, B. A. & Drazen, J. C. The rate of metabolism in marine animals: environmental constraints, ecological demands and energetic opportunities. Philos. Trans. R. Soc. B. 362, 2061–2078 (2007).CAS 

    Google Scholar 
    Seibel, B. A. et al. Oxygen supply capacity breathes new life into the critical oxygen partial pressure (Pcrit). J. Exp. Biol. 224, jeb242210 (2021).
    Google Scholar 
    Childress, J. J. & Seibel, B. A. Life at stable low oxygen: adaptations of animals to oceanic oxygen minimum layers. J. Exp. Biol. 201, 1223–1232 (1998).CAS 

    Google Scholar 
    Garcia, H. E., et al. World Ocean Atlas 2018, Volume 3: Dissolved Oxygen, Apparent Oxygen Utilization, and Oxygen Saturation (NOAA/NESDIS, 2019).Locarnini, R. A., et. al. World Ocean Atlas 2018, Volume 1: Temperature (NOAA/NESDIS, 2019).Maas, A. E., Frazar, S., Outram, D., Seibel, B. A. & Wishner, K. F. Fine-scale vertical distribution of macroplankton and micronekton in an eastern tropical North Pacific in association with an oxygen minimum zone. J. Plankton Res. 36, 1557–1575 (2014).
    Google Scholar 
    Rosa, R. & Seibel, B. A. Synergistic effect of climate-related variables suggests future physiological impairment in a top oceanic predator. Proc. Natl Acad. Sci. USA 52, 20776–20780 (2008).
    Google Scholar 
    Halsey, L. G., Killen, S. S., Clark, T. D. & Norin, T. Exploring key issues of aerobic scope interpretation in ectotherms: absolute versus factorial. Rev. Fish. Biol. Fish. 28, 405–415 (2018).
    Google Scholar 
    Peterson, C. C., Nagy, K. A. & Diamond, J. Sustained metabolic scope. Proc. Natl Acad. Sci. USA 87, 2324–2328 (1990).CAS 

    Google Scholar 
    Seibel, B. A., Luu, B. E., Tessier, S. N., Towanda, T. & Storey, K. B. Metabolic suppression in the pelagic crab, Pleuroncodes planipes, in oxygen minimum zones. Comp. Biochem. Physiol. A 224, 88–97 (2018).CAS 

    Google Scholar 
    Hadj-Moussa, H., Logan, S. M., Seibel, B. A. & Storey, K. B. Potential role for microRNA in regulating hypoxia-induced metabolic suppression in the jumbo squid? BBA Gene Regul. Mech. 1861, 586–593 (2018).CAS 

    Google Scholar 
    Torres, J. J. & Childress, J. J. Relationship of oxygen consumption to swimming speed in Euphausia pacifica. Mar. Biol. 74, 79–86 (1983).
    Google Scholar 
    Cohen, J. H. & Forward, R. B. Jr. Zooplankton diel vertical migration—a review of proximate control. Oceanogr. Mar. Biol. Annu. Rev. 47, 77–110 (2009).
    Google Scholar 
    Gilly, W. F. et al. Locomotion and behavior of Humboldt squid, Dosidicus gigas, in relation to natural hypoxia in the Gulf of California, Mexico. J. Exp. Biol. 215, 3175–3190 (2012).
    Google Scholar 
    Jaffe, J. S., Ohman, M. D. & De Robertis, A. Sonar estimates of daytime activity levels of Euphausia pacifica in Saanich inlet. Can. J. Fish. Aquat. Sci. 56, 2000–2010 (1999).
    Google Scholar 
    Klevjer, T. A. & Kaartvedt, S. Krill (Meganyctiphanes norvegica) swim faster at night. Limnol. Oceanogr. 56, 765–774 (2011).
    Google Scholar 
    Backus, R. H. et al. Ceratoscopelus maderensis: pecuiiar sound-scattering layer identified with this myctophid fish. Science 160, 991–993 (1968).CAS 

    Google Scholar 
    Barham, E. G. in Proceedings of an International Symposium on Biological Sound Scattering in the Ocean (ed. Farquhar, G. B.) 100–118 (Superintendent of Documents, 1971).Sanders, N. K. & Childress, J. J. A comparison of the respiratory function of the haemocyanins of vertically migrating and non-migrating pelagic, deep-sea Oplophorid shrimps. J. Exp. Biol. 152, 167–187 (1990).
    Google Scholar 
    Seibel, B. A. Critical depth in the jumbo squid, Dosidicus gigas (Ommastrephidae), living in oxygen minimum zones II. Blood-oxygen binding. Deep Sea Res. Part II 95, 139–144 (2013).CAS 

    Google Scholar 
    Pörtner, H.-O., Bock, C. & Mark, F. C. Oxygen- and capacity-limited thermal tolerance: bridging ecology and physiology. J. Exp. Biol. 220, 2685–2696 (2017).
    Google Scholar 
    Laffoley, D. & Baxter, J. M. Ocean Deoxygenation: Everyone’s Problem—Causes, Impacts, Consequences and Solutions (IUCN, 2019).Birk, M. A. Respirometry: Tools for Conducting and Analyzing Respirometry Experiments. R version 1.4.0 http://cran.r-project.org/package=respirometry (2021).Huang, B. et al. Improvements of the daily optimum interpolation sea surface temperature (DOISST) Version 2.1. J. Clim. 34, 2923–2939 (2021).
    Google Scholar  More

  • in

    Oceanic vertical migrators in a warming world

    Seibel, B. A. & Birk, M. A. Nat. Clim. Change https://doi.org/10.1038/s41558-022-01491-6 (2022).Article 

    Google Scholar 
    Urban, M. C. et al. Science 353, aad8466 (2016).Article 

    Google Scholar 
    Pörtner, H. O. & Knust, R. Science 315, 95–97 (2007).Article 

    Google Scholar 
    Verberk, W. C. E. P., Bilton, D. T., Calosi, P. & Spicer, J. I. Ecology 92, 1565–1572 (2011).Article 

    Google Scholar 
    Rubalcaba, J. G., Verberk, W. C., Hendriks, A. J., Saris, B. & Woods, H. A. Proc. Natl Acad. Sci. USA 117, 31963–31968 (2020).CAS 
    Article 

    Google Scholar 
    Deutsch, C. et al. Proc. Natl Acad. Sci. USA 119, e2201345119 (2022).CAS 
    Article 

    Google Scholar 
    Stramma, L. et al. Nat. Clim. Change 2, 33–37 (2012).CAS 
    Article 

    Google Scholar 
    Vergés, A. et al. Proc. R. Soc. B 281, 20140846 (2014).Article 

    Google Scholar  More

  • in

    Plankton response to global warming is characterized by non-uniform shifts in assemblage composition since the last ice age

    Brett, C. E. Sequence stratigraphy, paleoecology, and evolution: biotic clues and responses to sea-level fluctuations. Palaios 13, 241–262 (1998).Article 

    Google Scholar 
    Brett, C. E., Hendy, A. J. W., Bartholomew, A. J., Bonelli, J. R. & McLaughlin, P. I. Response of shallow marine biotas to sea-level fluctuations: a review of faunal replacement and the process of habitat tracking. Palaios 22, 228–244 (2007).Article 

    Google Scholar 
    Parmesan, C. Ecological and evolutionary responses to recent climate change. Annu. Rev. Ecol. Evol. Syst. 37, 637–669 (2006).Article 

    Google Scholar 
    Root, T. L. et al. Fingerprints of global warming on wild animals and plants. Nature 421, 57–60 (2003).CAS 
    PubMed 
    Article 

    Google Scholar 
    Walther, G.-R. et al. Ecological responses to recent climate change. Nature 416, 389–395 (2002).CAS 
    PubMed 
    Article 

    Google Scholar 
    Lenoir, J. et al. Species better track climate warming in the oceans than on land. Nat. Ecol. Evol. 4, 1044–1059 (2020).PubMed 
    Article 

    Google Scholar 
    Poloczanska, E. S. et al. Global imprint of climate change on marine life. Nat. Clim. Change 3, 919–925 (2013).Article 

    Google Scholar 
    Rillo, M. C., Woolley, S. & Hillebrand, H. Drivers of global pre‐industrial patterns of species turnover in planktonic foraminifera. Ecography 2022, e05892 (2021).Article 

    Google Scholar 
    Van der Putten, W. H., Macel, M. & Visser, M. E. Predicting species distribution and abundance responses to climate change: why it is essential to include biotic interactions across trophic levels. Phil. Trans. R. Soc. B 365, 2025–2034 (2010).PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Antão, L. H. et al. Temperature-related biodiversity change across temperate marine and terrestrial systems. Nat. Ecol. Evol. 4, 927–933 (2020).PubMed 
    Article 

    Google Scholar 
    Chen, I. C. et al. Asymmetric boundary shifts of tropical montane Lepidoptera over four decades of climate warming. Glob. Ecol. Biogeogr. 20, 34–45 (2011).Article 

    Google Scholar 
    García Molinos, J. et al. Climate velocity and the future global redistribution of marine biodiversity. Nat. Clim. Change 6, 83–88 (2015).Article 

    Google Scholar 
    Beaugrand, G., Edwards, M., Raybaud, V., Goberville, E. & Kirby, R. R. Future vulnerability of marine biodiversity compared with contemporary and past changes. Nat. Clim. Change 5, 695–701 (2015).Article 

    Google Scholar 
    Benedetti, F. et al. Major restructuring of marine plankton assemblages under global warming. Nat. Commun. 12, 5226 (2021).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Occhipinti-Ambrogi, A. Global change and marine communities: alien species and climate change. Mar. Pollut. Bull. 55, 342–352 (2007).CAS 
    PubMed 
    Article 

    Google Scholar 
    Williams, J. W. & Jackson, S. T. Novel climates, no-analog communities, and ecological surprises. Front. Ecol. Environ. 5, 475–482 (2007).Article 

    Google Scholar 
    Burrows, M. T. et al. Ocean community warming responses explained by thermal affinities and temperature gradients. Nat. Clim. Change 9, 959–963 (2019).Article 

    Google Scholar 
    Dornelas, M. et al. BioTIME: a database of biodiversity time series for the Anthropocene. Glob. Ecol. Biogeogr. 27, 760–786 (2018).PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Jonkers, L. et al. Integrating palaeoclimate time series with rich metadata for uncertainty modelling: strategy and documentation of the PalMod 130k marine palaeoclimate data synthesis. Earth Syst. Sci. Data 12, 1053–1081 (2020).Article 

    Google Scholar 
    Buitenhuis, E. T. et al. MAREDAT: towards a world atlas of MARine Ecosystem DATa. Earth Syst. Sci. Data 5, 227–239 (2013).Article 

    Google Scholar 
    Yasuhara, M., Tittensor, D. P., Hillebrand, H. & Worm, B. Combining marine macroecology and palaeoecology in understanding biodiversity: microfossils as a model. Biol. Rev. 92, 199–215 (2017).PubMed 
    Article 

    Google Scholar 
    Aze, T. et al. A phylogeny of Cenozoic macroperforate planktonic foraminifera from fossil data. Biol. Rev. 86, 900–927 (2011).PubMed 
    Article 

    Google Scholar 
    Takagi, H. et al. Characterizing photosymbiosis in modern planktonic foraminifera. Biogeosciences 16, 3377–3396 (2019).CAS 
    Article 

    Google Scholar 
    Schiebel, R. & Hemleben, C. Planktic Foraminifers in the Modern Ocean (Springer, 2017).Morey, A. E., Mix, A. C. & Pisias, N. G. Planktonic foraminiferal assemblages preserved in surface sediments correspond to multiple environment variables. Quat. Sci. Rev. 24, 925–950 (2005).Article 

    Google Scholar 
    Fenton, I. S., Pearson, P. N., Dunkley Jones, T. & Purvis, A. Environmental predictors of diversity in recent planktonic foraminifera as recorded in marine sediments. PLoS ONE 11, e0165522 (2016).PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Rutherford, S., D’Hondt, S. & Prell, W. Environmental controls on the geographic distribution of zooplankton diversity. Nature 400, 749–753 (1999).CAS 
    Article 

    Google Scholar 
    Tittensor, D. P. et al. Global patterns and predictors of marine biodiversity across taxa. Nature 466, 1098–1101 (2010).CAS 
    PubMed 
    Article 

    Google Scholar 
    Yasuhara, M., Hunt, G., Dowsett, H. J., Robinson, M. M. & Stoll, D. K. Latitudinal species diversity gradient of marine zooplankton for the last three million years. Ecol. Lett. 15, 1174–1179 (2012).PubMed 
    Article 

    Google Scholar 
    Yasuhara, M. et al. Past and future decline of tropical pelagic biodiversity. Proc. Natl Acad. Sci. USA 117, 12891–12896 (2020).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Jonkers, L., Hillebrand, H. & Kucera, M. Global change drives modern plankton communities away from the pre-industrial state. Nature 570, 372–375 (2019).CAS 
    PubMed 
    Article 

    Google Scholar 
    Beaugrand, G., Reid, P. C., Ibañez, F., Lindley, J. A. & Edwards, M. Reorganization of North Atlantic marine copepod biodiversity and climate. Science 296, 1692–1694 (2002).CAS 
    PubMed 
    Article 

    Google Scholar 
    Hinder, S. L. et al. Changes in marine dinoflagellate and diatom abundance under climate change. Nat. Clim. Change 2, 271–275 (2012).Article 

    Google Scholar 
    Southward, A. J., Hawkins, S. J. & Burrows, M. T. Seventy years’ observations of changes in distribution and abundance of zooplankton and intertidal organisms in the western English Channel in relation to rising sea temperature. J. Therm. Biol. 20, 127–155 (1995).Article 

    Google Scholar 
    Fenton, I. S. et al. Triton, a new species-level database of Cenozoic planktonic foraminiferal occurrences. Sci. Data 8, 160 (2021).PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Kucera, M., Rosell-Melé, A., Schneider, R., Waelbroeck, C. & Weinelt, M. Multiproxy approach for the reconstruction of the glacial ocean surface (MARGO). Quat. Sci. Rev. 24, 813–819 (2005).Kucera, M. et al. Reconstruction of sea-surface temperatures from assemblages of planktonic foraminifera: multi-technique approach based on geographically constrained calibration data sets and its application to glacial Atlantic and Pacific Oceans. Quat. Sci. Rev. 24, 951–998 (2005).Article 

    Google Scholar 
    Siccha, M. & Kucera, M. ForCenS, a curated database of planktonic foraminifera census counts in marine surface sediment samples. Sci. Data 4, 170109 (2017).PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Ezard, T. H. G., Aze, T., Pearson, P. N. & Purvis, A. Interplay between changing climate and species’ ecology drives macroevolutionary dynamics. Science 332, 349–351 (2011).CAS 
    PubMed 
    Article 

    Google Scholar 
    Fenton, I. S. et al. The impact of Cenozoic cooling on assemblage diversity in planktonic foraminifera. Phil. Trans. R. Soc. B 371, 20150224 (2016).PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Lowery, C. M. & Fraass, A. J. Morphospace expansion paces taxonomic diversification after end Cretaceous mass extinction. Nat. Ecol. Evol. 3, 900–904 (2019).PubMed 
    Article 

    Google Scholar 
    Wade, B. S., Pearson, P. N., Berggren, W. A. & Pälike, H. Review and revision of Cenozoic tropical planktonic foraminiferal biostratigraphy and calibration to the geomagnetic polarity and astronomical time scale. Earth Sci. Rev. 104, 111–142 (2011).Article 

    Google Scholar 
    Antell, G. S., Fenton, I. S., Valdes, P. J. & Saupe, E. E. Thermal niches of planktonic foraminifera are static throughout glacial-interglacial climate change. Proc. Natl. Acad. Sci. USA https://doi.org/10.1073/pnas.2017105118 (2021).Fauth, J. E. et al. Simplifying the jargon of community ecology: a conceptual approach. Am. Nat. 147, 282–286 (1996).Article 

    Google Scholar 
    Jackson, S. T. & Overpeck, J. T. Responses of plant populations and communities to environmental changes of the late Quaternary. Paleobiology 26, 194–220 (2000).Article 

    Google Scholar 
    Bard, E., Rostek, F., Turon, J.-L. & Gendreau, S. Hydrological impact of Heinrich events in the subtropical Northeast Atlantic. Science 289, 1321–1324 (2000).CAS 
    PubMed 
    Article 

    Google Scholar 
    Broecker, W. S. Massive iceberg discharges as triggers for global climate change. Nature 372, 421–424 (1994).CAS 
    Article 

    Google Scholar 
    Ruddiman, W. F. Late Quaternary deposition of ice-rafted sand in the subpolar North Atlantic (lat 40° to 65°N). Geol. Soc. Am. Bull. 88, 1813–1827 (1977).Article 

    Google Scholar 
    Hubbell, S. P. The Unified Neutral Theory of Biodiversity and Biogeography (Princeton Univ. Press, 2001).Liow, L. H., Van Valen, L. & Stenseth, N. C. Red Queen: from populations to taxa and communities. Trends Ecol. Evol. 26, 349–358 (2011).PubMed 
    Article 

    Google Scholar 
    Hillebrand, H. et al. Biodiversity change is uncoupled from species richness trends: consequences for conservation and monitoring. J. Appl. Ecol. 55, 169–184 (2018).Article 

    Google Scholar 
    Jackson, S. T. & Sax, D. F. Balancing biodiversity in a changing environment: extinction debt, immigration credit and species turnover. Trends Ecol. Evol. 25, 153–160 (2010).PubMed 
    Article 

    Google Scholar 
    Williams, J. W., Ordonez, A. & Svenning, J. C. A unifying framework for studying and managing climate-driven rates of ecological change. Nat. Ecol. Evol. 5, 17–26 (2021).PubMed 
    Article 

    Google Scholar 
    Van Meerbeeck, C. J., Renssen, H. & Roche, D. M. How did Marine Isotope Stage 3 and Last Glacial Maximum climates differ? Perspectives from equilibrium simulations. Clim. Past 5, 33–51 (2009).Article 

    Google Scholar 
    Jonkers, L. & Kučera, M. Global analysis of seasonality in the shell flux of extant planktonic Foraminifera. Biogeosciences 12, 2207–2226 (2015).Article 

    Google Scholar 
    Ofstad, S. et al. Development, productivity, and seasonality of living planktonic foraminiferal faunas and Limacina helicina in an area of intense methane seepage in the Barents Sea. J. Geophys. Res. Biogeosci. 125, e2019JG005387 (2020).CAS 
    Article 

    Google Scholar 
    Bova, S., Rosenthal, Y., Liu, Z., Godad, S. P. & Yan, M. Seasonal origin of the thermal maxima at the Holocene and the last interglacial. Nature 589, 548–553 (2021).CAS 
    PubMed 
    Article 

    Google Scholar 
    Rillo, M. C. et al. On the mismatch in the strength of competition among fossil and modern species of planktonic Foraminifera. Glob. Ecol. Biogeogr. 28, 1866–1878 (2019).Article 

    Google Scholar 
    Lisiecki, L. E. & Stern, J. V. Regional and global benthic δ18O stacks for the last glacial cycle. Paleoceanography 31, 1368–1394 (2016).Article 

    Google Scholar 
    Reimer, P. J. et al. IntCal13 and Marine13 radiocarbon age calibration curves 0–50,000 years cal BP. Radiocarbon 55, 1869–1887 (2013).CAS 
    Article 

    Google Scholar 
    Butzin, M., Köhler, P. & Lohmann, G. Marine radiocarbon reservoir age simulations for the past 50,000 years. Geophys. Res. Lett. 44, 8473–8480 (2017).CAS 
    Article 

    Google Scholar 
    Langner, M. & Mulitza, S. Technical Note: PaleoDataView—A software toolbox for the collection, homogenization and visualization of marine proxy data. Clim 15, 2067–2072 (2019).
    Google Scholar 
    Mix, A. C., Bard, E. & Schneider, R. Environmental processes of the ice age: land, oceans, glaciers (EPILOG). Quat. Sci. Rev. 20, 627–657 (2001).Article 

    Google Scholar 
    Osman, M. B. et al. Globally resolved surface temperatures since the Last Glacial Maximum. Nature 599, 239–244 (2021).CAS 
    PubMed 
    Article 

    Google Scholar 
    Horn, H. S. Measurement of ‘overlap’ in comparative ecological studies. Am. Nat. 100, 419–424 (1966).Article 

    Google Scholar 
    Jost, L., Chao, A. & Chazdon, R. L. in Biological diversity: frontiers in measurement and assessment (eds Anne E. Magurran & Brian J. McGill) 66–84 (Oxford University Press, 2011).Ferrier, S., Manion, G., Elith, J. & Richardson, K. Using generalized dissimilarity modelling to analyse and predict patterns of beta diversity in regional biodiversity assessment. Divers. Distrib. 13, 252–264 (2007).Article 

    Google Scholar 
    Shannon, C. E. A mathematical theory of communication. Bell Syst. Tech. J. 27, 379–423 (1948).Article 

    Google Scholar 
    R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2022).Wickham, H. et al. Welcome to the Tidyverse. J. Open Source Softw. 4, 1686 (2019).Article 

    Google Scholar 
    Firke, S. janitor: Simple tools for examining and cleaning dirty data. R package version 2.1.0 https://CRAN.R-project.org/package=janitor (2021).Oksanen, J. et al. vegan: Community ecology package. R package version 2.5-7 https://CRAN.R-project.org/package=vegan (2020).Hallett, L. M. et al. codyn: an R package of community dynamics metrics. Methods Ecol. Evol. 7, 1146–1151 (2016).Article 

    Google Scholar 
    Juggins, S. rioja: Analysis of quaternary science data. R package version 0.9-26 https://cran.r-project.org/package=rioja (2020).Lê, S., Josse, J. & Husson, F. FactoMineR: an R package for multivariate analysis. J. Stat. Softw. 25, 1–18 (2008).Article 

    Google Scholar 
    Wickham, H. ggplot2: Elegant Graphics for Data Analysis (Springer, 2016).Hijmans, R. J. raster: Geographic data analysis and modeling. R package version 3.4-13 https://CRAN.R-project.org/package=raster (2021).Garnier, S. viridis: Default color maps from ‘matplotlib’. R package version 0.6.1 https://CRAN.R-project.org/package=viridis (2021.)Locarnini, R. A. et al. World Ocean Atlas 2018, Vol. 1: Temperature. NOAA Atlas NESDIS 81 (NOAA, 2019). More