Reconstructing Long-Term Changes in Avian Populations Using Lake Sediments: Opening a Window Onto the Past

The lack of long-term monitoring data for many wildlife populations is a limiting factor in establishing meaningful and achievable conservation goals. Even for well-monitored species, time series are often very short relative to the timescales required to understand a population’s baseline conditions before the contemporary period of increased human impacts. To fill in this critical information gap, techniques have been developed to use sedimentary archives to provide insights into long-term population dynamics over timescales of decades to millennia. Lake and pond sediments receiving animal inputs (e.g., feces, feathers) typically preserve a record of ecological and environmental information that reflects past changes in population size and dynamics. With a focus on bird-related studies, we review the development and use of several paleolimnological proxies to reconstruct past colony sizes, including trace metals, isotopes, lipid biomolecules, diatoms, pollen and non-pollen palynomorphs, invertebrate sub-fossils, pigments, and others. We summarize how animal-influenced sediments, cored from around the world, have been successfully used in addressing some of the most challenging questions in conservation biology, namely: How dynamic are populations on long-term timescales? How may populations respond to climate change? How have populations responded to human intrusion? Finally, we conclude with an assessment of the current state of the field, challenges to overcome, and future potential for research.


INTRODUCTION
Traditionally, wildlife monitoring programs require continuous and sustained investment, and are thus costly, labor-intensive, and potentially distressing for the organism (Witmer, 2005). Furthermore, although modern monitoring methods provide estimates or indices of contemporary population size, most are relatively recent in implementation and thus provide limited insights into population trends that naturally change over many generations and varying ecological conditions. As a result, resources are typically allocated to populations that are currently of interest or concern, and thus most management programs are responsive to population changes as opposed to predictive of population change. This unfortunate lack of long-term data is well established. For example, Bonebrake et al. (2010) demonstrated that only 15% of studies examining population declines used data that were older than 100 years and that only 42% used monitoring data collected annually. In the same vein, several reviews have reiterated the importance of long-term perspectives to effectively address various environmental issues and to formulate effective decision-making in conservation biology (Willis and Birks, 2006;Froyd and Willis, 2008;Lotze and Worm, 2009;Dietl and Flessa, 2011;Seddon et al., 2014;Dietl et al., 2015).
As with most wildlife, many bird populations around the world are in decline (e.g., Şekercioǧlu et al., 2004;Inger et al., 2014;Paleczny et al., 2015;Beresford et al., 2019). Since 1970, there has been a staggering net loss of nearly three billion birds in North America (Rosenberg et al., 2019), and this estimate does not incorporate the known declines in seabird populations (Paleczny et al., 2015). Compounding this growing conservation issue is the reality that, for the majority of species, there is a lack of adequate long-term monitoring data to effectively consider baseline population size or natural variability when establishing management strategies and targets. As a result, it is very challenging to disentangle whether a species is declining as a result of anthropogenic activities and thus requires human intervention and/or whether the population is within the range of long-term natural variability and thus requires a different approach to maintain natural dynamics.
In response to the dearth of long-term monitoring data, a variety of methods have been developed to identify long-term trends in bird populations and diet, termed "paleornithology" (reviewed in Wood and De Pietri, 2015). To identify the demographic history on an evolutionary time scale (i.e., >10,000 years), ancient DNA (or aDNA) has been used. For example, genetic methods have been successfully employed to identify the population trends of the extinct passenger pigeon (Ectopistes migratorius) (Hung et al., 2014), the population recovery of the crested ibis (Nipponia nippon) (Feng et al., 2019), and population trends in Adélie penguins (Pygoscelis adeliae) and emperor penguins (Aptenodytes forsteri) in response to climate changes (Li et al., 2014). However, aDNA is less effective for more recent population trends (i.e., <10,000 years) because insufficient time has elapsed for the accumulation of evidence of evolutionary change. In the critical period from the present until ∼10,000 years BP (before present, where 0 is 1950 CE), we argue that paleolimnology, or the study of lake and pond sediments as historical archives, can be a highly effective approach to access long-term population demographics. In this review, we provide an accessible summary of current developments and future directions on how paleoecological approaches can be useful to reconstruct long-term animal trends, with a focus on avian colonies. We provide an overview of the paleolimnological approach and its methodologies, geared toward researchers who are unfamiliar with these techniques. Also, we present a variety of case studies that address critical questions that are broadly important in modern conservation biology. Ultimately, we believe that combining paleoecological techniques with modern wildlife monitoring approaches will provide new perspectives to address pressing conservation questions, help identify drivers of population declines, understand natural population dynamics, and provide a better understanding of how populations have changed in the past to establish predictions of future population trends.

What Is Paleolimnology?
Paleolimnology is the study of the physical, chemical, and biological information stored in lake and river sediments (Smol, 2008). Sediments provide archives of various indicators of environmental change within a given body of water, from the surrounding catchment, and from atmospheric deposition. As such, lake sediments can provide long-term perspectives that predate historical records, and thus can be used to help place recent monitoring data into an ecologically relevant temporal context.
The use of paleolimnology to reconstruct past population dynamics is built on the well-established role of certain animals as effective biological vectors of nutrients and other products that are capable of modifying their environment (e.g., Blais et al., 2007). These biological vectors, typically shortened to "biovectors, " can have various influences on the environment, broadly classified as trophic level effects and transport effects (Bauer and Hoye, 2014). Trophic level effects occur when large groups of organisms alter the food web through changing the availability of a resource, such as desert locusts (Schistocerca gregaria) consuming and devastating crops (Sánchez-Zapata et al., 2007) or zebra mussels (Dreissena polymorpha) selectively filtering desirable algae from the water, resulting in blooms of undesirable and toxic Microcystis (Vanderploeg et al., 2001). Transport effects are changes in an ecosystem resulting from the collection, concentration, and transportation of large amounts of nutrients, plant propagules, and parasites from one location to another. For example, anadromous fishes and colonial seabirds transport and concentrate nutrients, metals, and contaminants from their aquatic feeding grounds to their terrestrial breeding territories (e.g., Evenset et al., 2007a,b;Michelutti et al., 2009;Holtgrieve and Schindler, 2011). Because of the vast influences of biovectors on the environment, many are considered "ecosystem engineers, " or organisms that biotically or abiotically change the availability of resources for other species (Jones et al., 1994).
Many of the ecosystem changes linked to biovectors are preserved in nearby bodies of water and can be reconstructed to understand changes in past biovector abundance (Figure 1). Briefly, a site that receives substantial waste products (either washed or blown in or from direct inputs) to a body of water is selected. Lacustrine environments (i.e., standing water) are preferred, as the sediments tend to accumulate chronologically with minimal mixing, in contrast to those in fluvial environments (i.e., running water) (Smol, 2008). There is also a limited but growing body of literature using similar paleolimnological approaches on peat (Outridge et al., 2016;Davidson et al., 2018;Groff et al., 2020), directly from nests (Burnham et al., 2009), and on guano deposits themselves (Nocera et al., 2012;Gallant et al., 2020Gallant et al., , 2021. Next, the sample is collected using a sediment corer, which, in its simplest form, is a tube inserted into the sediment to retrieve an undisturbed vertical profile of sediment [see Glew et al. (2001) for a detailed overview]. The sediment core is then sectioned into discrete intervals and dated using a variety of physical [e.g., laminated sediment layers (varves) or volcanic ash layers (tephra)] and radioisotopic (e.g., 210 Pb, 137 Cs, 14 C) techniques [see Appleby (2001) and Björck and Wohlfarth (2001) for detailed overviews]. Finally, sediments are analyzed for a variety of proxies that directly (e.g., geochemistry, δ 15 N, lipid biomolecules) and indirectly (e.g., diatoms, invertebrate subfossils, pollen and non-pollen palynomorphs) track the presence and influences of the biovectors (Figure 2). Typically, multiple proxies are used for a "weight-of-evidence" approach (Smol, 2008) because increasing the number of corroborative proxies strengthens the accuracy of the reconstruction and minimizes the confounding impacts of extraneous influences.
Some of the earliest paleolimnological reconstructions of biovectors linked changes in sediment geochemistry to avian inputs to make inferences about past Adélie penguin rookery size (Zale, 1994). The technique was then expanded across the Antarctic in a series of studies aimed at understanding the long-term population changes in response to past climate (e.g., Sun et al., 2000Sun et al., , 2004Huang et al., 2009a). In the Northern Hemisphere, similar reconstructions conducted on Pacific sockeye salmon (Oncorhynchus nerka) identified striking fluctuations in numbers that were correlated to climate change and other environmental stressors (Finney et al., 2000(Finney et al., , 2002Gregory-Eaves et al., 2003;Selbie et al., 2007). Afterward, work from Cape Vera (Devon Island, Nunavut) in the Canadian High Arctic successfully expanded the field by demonstrating the utility of a variety of proxies that could be successfully used to reconstruct seabird numbers with sediments, including contaminants, stable isotopes (e.g., δ 15 N), invertebrate subfossils, and lipid biomolecules (Blais et al., 2007;Evenset et al., 2007a,b;Michelutti et al., 2009Michelutti et al., , 2011Cheng et al., 2016Cheng et al., , 2021. Since these studies, the use of sedimentary archives in biovector work has rapidly expanded to other organisms (e.g., Navarro et al., 2017;Gallant et al., 2020) and spatially to other regions (Supplementary Table 1).
Although a variety of animal groups can be tracked using sediment archives (e.g., butterflies and moths, fish, birds, bats), the focus of this review will be on aquatic birds, a globally widespread biovector with behaviors that are conducive to paleolimnological studies. Aquatic birds, especially seabirds, are an ideal study group for population reconstructions as 95% of seabirds are colonial (Danchin and Wagner, 1997), thus forming large congregations that concentrate their waste products (Otero et al., 2018), which can be archived in lakes. Seabirds are also highly visible sentinel species of marine ecosystem functioning, and therefore changes in observed seabird numbers may provide useful information about unobserved components of the marine environment (e.g., Parsons et al., 2008;Hazen et al., 2019;Velarde et al., 2019). As a result, the well-documented, global decline in seabird numbers highlights the pressing need for information that can lead to effective conservation and management. In Figure 2, we depict some of the key paleolimnological proxies that will be discussed in this review.

Bio-Elements
One of the most useful paleolimnological proxies for reconstructing bird colony size is sedimentary elements that are primarily of biogenic origins. As seabirds typically occupy upper trophic positions of the food web, they may carry high amounts of elements as a result of the concentrating effects of bioaccumulation and biomagnification (Mallory and Braune, 2012). Importantly, the cocktail of elements and their concentrations vary depending on a variety of factors, including age, species, location of feeding grounds, and, especially, diet and trophic position (Honda et al., 1990;Lock et al., 1992;Michelutti et al., 2010b). Based on these factors, various seabird species may introduce differing combinations of elements (often termed "bio-elements") to their breeding grounds that can be used to track changes in colony size with sediment records. This technique assumes that the diet and trophic position of the tracked species remains consistent and that the elements excreted by the study species remain similar throughout the record.
Elements accumulated in a bird's body can be released to the local environment via feces, dropped or regurgitated food and stomach oils, as well as eggshells, feathers, and carcasses of deceased adults or young (Liu et al., 2006a;De La Peña-Lastra, 2021). These products then ultimately drain into a nearby body of water, where they can be tracked in dated lake sediments. As such, changes in the concentrations of trace metals linked to avian deposition can inform changes in the colony size.
Element concentrations, however, do not have a one-to-one relationship with colony size, and a variety of factors must be considered when making inferences about numbers. For example, it is important to consider natural geological sources of metals in the sediments such as the soils of the catchment, hydrology, and sedimentation conditions, as well as potential post-depositional factors that may confound interpretations (Boës et al., 2011). Therefore, normalization is often used to distinguish which portion of the measured sedimentary metal load is most likely derived from the seabirds. A twotiered approach is recommended. First, sieving is suggested to minimize the influences of grain size. Next, elements-of-interest should be normalized against a conservative lithogenic element (i.e., one that is geochemically stable and therefore changes in its concentration would reflect post-depositional processes, such as erosion) (Kersten and Smedes, 2002). Normalization simply entails dividing the element of interest by a conservative lithogenic element. Effective elements for normalization include aluminum (Al), lithium (Li), rubidium (Rb), titanium (Ti), and zirconium (Zr) (Loring, 1991;Boës et al., 2011), with Al and Ti the most used in paleolimnological studies (Roberts et al., 2017;Gao et al., 2018c;Duda et al., 2020a,c).
Once normalized, it is important to distinguish which elements are representative of seabird inputs and which are FIGURE 1 | A schematic showing how biovectors can influence inland bodies of water and how sediments can be used to reconstruct past colony dynamics. (A) Biovectors, in this case storm-petrels, typically occupy a high-trophic level and are gregarious on land. When feeding in the open ocean, they collect and concentrate marine nutrients (i.e., N, P) and other products, including stable isotopes (e.g., δ 15 N), bio-elements (e.g., Cd, Zn, Hg), and contaminants. (B) Biovectors transport these nutrients and other products from their feeding grounds to terrestrial breeding colonies. (C) The nutrients and other products are deposited in high concentrations onto the landscape via wastes in the forms of feces, dropped or regurgitated food, eggshells, feathers, and carcasses. (D) In the presence of a receiving body of water, such as a lake or pond, the waste products alter the water chemistry and, subsequently, the biotic assemblages, which are then archived in a waterbody's sediments, thus creating a continuous long-term record of changes in the colony size. By retrieving a dated sediment core from the biovector-influenced pond and analyzing it for paleolimnological proxies, it is possible to reconstruct changes in past colony size.
from the catchment. This can be done with R-cluster analysis (Huang et al., 2009a), biogenic enrichment factors , or an index of geoaccumulation (Müller, 1979;Förstner et al., 1993). R-cluster analysis (sometimes called R mode cluster analysis") is a statistical technique used to partition and identify similar groups of elements and infer their origin. With an understanding of the elements introduced by the species of interest, it is possible to identify a group of elements of the same origin. As demonstrated in Figure 3, Liu et al. (2005) examined 24 elements from a sediment core collected in Ardley Island in the Antarctic that was influenced by a major colony of gentoo penguins (Pygoscelis papua), Adélie penguins, and chinstrap penguins (Pygoscelis antarctica). Using R-cluster analysis, followed by a Pearson correlation, the authors identified phosphorus (P), copper (Cu), strontium (Sr), zinc (Zn), selenium (Se), calcium (Ca), fluorine (F), and sulfur (S) as bio-elements with the same origins that were highly correlated to one another (Figure 3). By comparing this grouping to previously published chemical analyses of penguin droppings, the authors confirmed that the group of bio-elements is representative of the species of interest, and so could be used to understand changes in penguin numbers back in time. R-cluster analysis, which has been used worldwide, is often considered the most effective tool to identify ornithogenic elements if a nearby reference site (i.e., a nearby lake without avian inputs) is not available to distinguish the natural, background element concentrations (Zale, 1994;Sun et al., 2000Sun et al., , 2004Liu et al., 2005Liu et al., , 2006cHuang et al., 2009aHuang et al., , 2011Xu et al., 2011Xu et al., , 2016Xu et al., , 2019Xu et al., , 2020Roberts et al., 2017).
In situations with geochemical measurements from a reference site and a biovector fecal sample, it is possible to use biogenic enrichment factors (B) to more specifically identify the elements introduced by the target species ). B is calculated as the ratio of the concentration of an element in the feces to the concentration of the same element in the surface sediment of a reference site. Thus, any B value greater than one indicates that the concentration in feces is higher than in the reference conditions, and therefore the element concentration in the sediment is likely to be of ornithogenic origin, specific to the focal species. This technique provides more specific markers of avian inputs compared to R-cluster analysis, but it is less commonly used given that fecal samples are also needed. Biogenic enrichment factors have been useful in the Canadian High Arctic to track changes in northern fulmar (Fulmarus glacialis), where Brimble et al. (2009) identified arsenic (As), cadmium (Cd), potassium (K), P, and Zn as primarily fulmarderived elements. Similar approaches were used on islands in the Northwest Atlantic Ocean to link Cd and Zn to inputs from Leach's storm-petrel (Hydrobates leucorhous) (Duda et al., 2020a,c).
The final method to identify the bio-elements introduced by a seabird is with an index of geoaccumulation (I geo , sometimes "geoaccumulation index"). This method is similar to biogenic enrichment factors, although it uses bird-influenced sediments FIGURE 2 | A visual representation of the main paleolimnological proxies discussed in this review that can be used to reconstruct past biovector numbers. Increases and decreases in colony size, and thus the amount of biovector waste inputs, can result in changes in the relative abundance and prevalence of the proxies, including (A) bio-elements; (B) stable isotope ratios; (C) lipid biomolecules; (D) diatoms; (E) invertebrate sub-fossils; (F) terrestrial pollen and non-pollen palynomorphs; and (G) pigments.
instead of feces to compare with reference conditions (Müller, 1979;Förstner et al., 1993). The equation for I geo for each element of interest is given as where C n is the element concentration in the biovector influenced sediment, B n is the element concentration in local bedrock, and 1.5 is a factor to correct for possible variations in lithology. The result is then compared to a table to infer the level of contamination of the element by the seabirds (Table 1).
Using I geo , inputs from emperor penguins were linked to the enrichment of mercury (Hg), total phosphorus, Zn, Se, and Cd (Huang et al., , 2016). An important consideration when selecting the elements to track changes in avian populations is the potential interactions of the element with sediment composition and limnological conditions. For example, although P is highly elevated in most avian guano, it is also affected by desorption and resorption reactions in sediment (Koski-Vähälä et al., 2001) and can be highly mobile in anoxic sediments (Ginn et al., 2012). These events are triggered by the release of P bound to iron (Fe)-and Al-rich sediments, which can occur during periods of elevated pH, increased silicon, and anaerobic conditions (Koski-Vähälä et al., 2001). This P release, when not accounted for, can result in erroneous interpretations of the seabird colony size. Other elements, such as As and Zn, may also become labile in sediments in response to changes in pH, salinity, temperature, organic matter, and, in particular, redox conditions (Couture et al., 2010;Wang et al., 2015;Kouassi et al., 2019). Although bio-elements can provide an effective proxy for avian inputs, it is important to always consider any potential post-depositional effects when interpreting the results.

Isotope Ratios
Stable isotope analysis has a long history of use in paleolimnological reconstructions (Leng, 2006). The use of stable isotopes is predicated on the observation that the ratio of the heavier and lighter isotopes changes based on their source. In paleolimnological studies, the most common stable isotopes investigated are nitrogen (N) and carbon (C). Their isotopic ratios (R), 15 N/ 14 N and 13 C/ 12 C, respectively, are typically expressed as delta (δ) notation against an established standard as below.
The number is then multiplied by 1,000 to be expressed as "per mille ( )." The standard notation for N and C is, therefore, δ 15 N and δ 13 C.
FIGURE 3 | A workflow diagram illustrating how one may evaluate which elements are likely related to a biovector using R-cluster analysis. This example uses data from Liu et al. (2005) where authors determined that S, P 2 O 5 , F, Zn, Se, Sr, CaO, and Cu were related to penguin inputs. (A) R-cluster analysis creates a dendrogram of all measured elements. Major groups are determined based on their linkage distance and compared to known fecal data from the bird-of-interest to determine which grouping represents bio-elements. (B) A Pearson correlation is used to corroborate that the elements are of similar origins. (C) The resulting grouping of elements are positively correlated and represent input from the bird-of-interest. Higher contamination intensity equates to higher enrichment by the biovector.
δ 15 N values in sediments influenced by biovectors are enriched relative to reference conditions (i.e., atmospheric δ 15 N) because the heavier 15 N isotope is retained in tissues by organisms during amino acid decomposition, whereas the lighter 14 N isotope, because of its weaker bond with adjacent atoms, is preferentially excreted during transamination and deamination reactions (Adams and Sterner, 2000). This results in stepwise enrichment up the food web, averaging an increase of ∼3.4 per trophic level (Minagawa and Wada, 1984). In paleolimnology, the enriched δ 15 N is deposited onto the landscape by a nearby bird colony or roost and is archived in sediments via avian guano carcasses, feathers, and eggshells. Therefore, a higher δ 15 N value in sediments typically relates to greater bird numbers until a plateau is reached equal to the δ 15 N value of the refuse (Duda et al., 2020a), with minor changes from environmental factors such as ammonification, remineralization, nitrification, and denitrification (e.g., Cifuentes et al., 1989;Savage et al., 2004;Torres et al., 2012). Since δ 15 N is linked to the trophic status of the biovectors, the relative increase compared to background conditions is higher in upper-trophic-level organisms than in lower-trophic-level organisms (Michelutti et al., 2010b).
Although an indirect measure of avian inputs, changes in sediment δ 13 C values can also be used to track changes in seabird numbers. Instead of directly measuring inputs from guano, as is the case with δ 15 N, deviations in δ 13 C values may track changes in production linked with marine-derived changes in primary productivity. However, when considering recent changes in δ 13 C values, it is important to account for depleting atmospheric δ 13 C due to increased CO 2 production from fossil fuel burning (known as the Suess Effect; Keeling, 1979). A common standardization method is to normalize to pre-1840 levels of δ 13 C, based on the methods presented in Schelske and Hodell (1995). Once normalized, δ 13 C can have both a positive and a negative relationship with changes in the bird numbers, and therefore a thorough understanding of the system and the use of supporting proxies is recommended. Typically, δ 13 C values are positively correlated to colony size (i.e., δ 13 C is enriched as colony size increases). This relationship principally reflects an increase in primary production and organic matter introduction resulting from guano fertilization (Sun et al., 2000;Yuan et al., 2010;Liu et al., 2013;Conroy et al., 2015;Davidson et al., 2018;Gao et al., 2018c;Cheng et al., 2021). However, in some cases, δ 13 C values can also have a negative relationship with colony size (i.e., δ 13 C is depleted as colony size increases). For example, Liu et al. (2006b) linked the increase in penguin colony size to δ 13 C depletion, which was a result of vegetation trampling by the penguins and toxification from their guano, leading to a diminishment of the overall vegetation. Similarly, Duda et al. (2020a) noted that burrowing Leach's storm-petrels caused an accelerated rate of inorganic material delivery to the lake, which caused a depletion of δ 13 C as the colony increased in size.
Although less common, studies have also demonstrated the potential of HCl-soluble strontium isotopes to track seabird guano inputs, and therefore changes in long-term population size (as 87 Sr/ 86 Sr; Sun et al., 2005;Liu et al., 2007). Strontium isotopes are effective ornithological proxies that track the source of nutrients to an environment without transformation by biological or chemical processes (Miller et al., 1993). In penguin-influenced sediments from Antarctica, Sun et al. (2005) determined that the hydrochloric acid (HCL)-soluble 87 Sr/ 86 Sr reflected penguin guano ratios, whereas the HCl-insoluble 87 Sr/ 86 Sr reflected bedrock ratios. These results, supported by a follow-up study (Liu et al., 2007), suggest that using HClsoluble 87 Sr/ 86 Sr can be used to track marine-derived nutrient enrichment from seabirds.

Lipid Biomolecules
Lipid biomolecules are water-insoluble fats produced by all biological life. Of particular interest in paleolimnological reconstructions are sterols and stanols. Sterols are natural, unsaturated steroid alcohols that represent an important group of compounds for the functioning and structure of biological membranes (Bretscher, 1973). Stanols are the reduction products of sterols by the process of hydrogenation, which is typically mediated by microbes in animal gastrointestinal tracts or in the environment (Reeves and Patton, 2001;Bull et al., 2002). Sterols and stanols are discharged in varying concentrations and combinations in the feces of animals (Leeming et al., 1996;Harrault et al., 2019) and are emerging as important proxies to directly track animal inputs to lake sediments over decadal-and century-scale periods (D'Anjou et al., 2012;Hargan et al., 2019;Schroeter et al., 2020;Cheng et al., 2021). In general, sterols and stanols are hydrophobic compounds with a high affinity for sorption onto organic matter and suspended particles, and therefore they typically become rapidly incorporated into sediments of aquatic ecosystems (Korosi et al., 2015). Once in the sediments, the degradation of fecal sterols is likely minimal (Bartlett, 1987;Benfenati et al., 1994), which makes them ideal for tracking changes in the abundance of seabirds over time.
Recent proof-of-concept studies have demonstrated that the ratio of cholesterol (cholest-5-en-3β-ol) to the sum of cholesterol and sitosterol (β-sitosterol), termed the "seabird index, " tracked significant enrichment in pond sediments closest to a nesting colony of northern fulmars in the High Arctic (Cheng et al., 2016(Cheng et al., , 2021. This ratio is used to control for cholesterol in the sediment that is introduced from higher plants and algae as opposed to having ornithogenic origins (Cheng et al., 2016). Furthermore, enrichment in a similar seabird index was observed in temperate ponds receiving fecal inputs from waterbirds with higher trophic positions (Hargan et al., 2018). Unlike in the Arctic and sub-Arctic locations, where transformations of sterols in the environment are largely limited by the cold temperatures and short growing seasons, sterols in temperate locations may undergo hydrogenation to stanols in the environment. Thus, Hargan et al. (2018) expanded the seabird index for tracking waterbird (i.e., double-crested cormorant (Phalacrocorax auritus) and ring-billed gull (Larus delawarensis)) populations from the sediments of shallow, temperate ponds, by including cholesterol and sitosterol and their microbially reduced products, 5αcholestanol and stigmastanol, respectively. Overall, the use of these new seabird sterol indices has been validated through significant correlations with other traditional and independent seabird proxies (Cheng et al., 2016;Hargan et al., 2018;Yang et al., 2021). These studies indicate that both cholesterol concentrations and the "seabird sterol indices" can provide an indirect measure of seabird presence in sediment cores from lakes and ponds where seabird-derived enrichment is evident. These lipid proxies, especially when considered as part of a multi-proxy framework that also includes δ 15 N and biogenic elements, can be useful to infer avian population trends through time Duda et al., 2020b,c;Yang et al., 2021).

Diatoms
Bird colonies can dramatically alter the limnological conditions of receiving waterbodies (Otero et al., 2018). These limnological changes typically include increased concentrations of nutrients, metals, and dissolved organic carbon. Diatoms (Class: Bacillariophyceae) are microscopic, single-celled eukaryotic algae found in most wet habitats in the world. Changes in diatom assemblage composition are the most commonly used approach to track changes in water quality in paleolimnological assessments, as they are ideal indicators of ecosystem change (Smol and Stoermer, 2010) because their taxonomically identifiable siliceous cell walls are typically well-preserved in a broad spectrum of sediments (Battarbee et al., 2001). Importantly, diatom taxa have well-documented environmental optima and tolerances, which can be used to infer changes in nutrients, pH and other limnological variables potentially influenced by bird inputs.
Typically, with the introduction and/or growth of a seabird colony influence, the diatom assemblage inhabiting a body of water will shift from taxa indicative of low nutrients to those indicative of high nutrients (e.g., Brooks et al., 2012;Stewart et al., 2015Stewart et al., , 2019. In other instances, depending on the species producing the seabird guano and conditions of the surrounding ecosystem, the diatoms can also reflect changes in pH (Duda et al., 2020a,c) or conductivity (Duda et al., 2019). Given their ubiquity, diatoms have been successfully used to draw inferences about historical seabird populations around the world (Michelutti et al., 2009;Keatley et al., 2011;Brooks et al., 2012;MacDonald et al., 2015;Stewart et al., 2015Stewart et al., , 2019Roberts et al., 2017;Wasiłowska et al., 2017;Luoto et al., 2019;Duda et al., 2020a,c).

Invertebrate Sub-Fossils
The most common invertebrate indicators used in paleolimnology are chironomids (Order: Diptera), Cladocera, and Ostracods (Class: Crustacea). The remains of other invertebrates and larger organisms, such as fish, are often not abundant in sediments and typically do not have enough direct ecological link to seabird inputs to act as effective paleolimnological proxies.
Chironomid larva are the most abundant insect remains in lake sediments (Hall and Herhardt, 2002). Paleolimnologists use the distinct chitinous mouthparts and head capsule characteristics to identify remains. As ecological indicators for biovector reconstructions, chironomid communities are predominantly used to infer changes in bottom-water oxygen, acidity, and metal concentrations (Walker, 2001). As a bird colony grows in size, several limnological changes that affect the assemblage composition of chironomid larvae may occur. For example, increased nutrients typically result in increased algal biomass (i.e., increased organic matter), which in turn may increase benthic respiration/decomposition, leading to lower deepwater oxygen concentrations or even anoxia. Chironomids can also be affected by changes in pH and metal concentrations introduced by seabird guano. Changes in chironomid subfossil assemblages have been used to infer long-term changes in different seabird populations Luoto et al., 2019;Duda et al., 2020a,c;Lim et al., 2021), but often, in oligotrophic ecosystems, production is so low that head capsule remains are too few to adequately characterize assemblage composition (Quinlan and Smol, 2001). In such low productivity ecosystems, increases in head capsule concentration and diversity have been related to enhanced food availability from seabird nutrient introductions (Michelutti et al., 2009;Brooks et al., 2012;Luoto et al., 2014Luoto et al., , 2015Stewart et al., 2020).
Compared to chironomids, other invertebrate sub-fossils have been sparsely used to reconstruct changes in lake water characteristics linked to bird populations. Ostracod communities respond to changes in nutrient concentrations, but their identifying shells rarely preserve well in acidic conditions (Holmes, 2001). Similarly, cladoceran community structure can respond to habitat changes and lake trophic status. To date, however, higher abundances of Cladocera have mainly been linked to track increased food availability (i.e., higher productivity) resulting from avian nutrient enrichment (Griffiths et al., 2010;Luoto et al., 2014Luoto et al., , 2015.

Terrestrial Pollen and Non-pollen palynomorphs
Avian nutrient inputs typically also have a significant influence on terrestrial ecosystems; however, the nature and magnitude of these influences can be highly variable depending on the quantity and quality of the inputs, as well as the receiving environment. In some environments, avian-derived nutrients enhance productivity, resulting in a more fertilized, productive, and diverse ecosystem (Anderson and Polis, 1999;Croll et al., 2005). Alternatively, too much nutrient deposition can result in reductions in vegetation and seed bank richness, as observed on some islands inhabited by double-crested cormorants (Boutin et al., 2011).
Regardless of the type of influence, changes to the terrestrial environment preserved in nearby lake sediments can be used to make inferences about the seabird population size. There are two methods to understand past changes in vegetation: palynology and macrofossil analysis. Palynology is the study of microscopic pollen and non-pollen palynomorphs (NPP, or microfossils on pollen slides that are not pollen, taxonomically belonging to a wide variety of fungi, animals, algae, and plants), while macrofossil analysis examines a variety of vegetative remains, such as seeds, twigs, and leaves. Generally, palynology is preferred over macrofossil analysis as pollen and NPP are much more abundant and taxonomically identifiable in sediments. However, macrofossils are occasionally used in paleolimnological studies to understand changes in local vegetation (reviewed in Birks, 2001).
Various paleolimnological studies have used palynology to make inferences about long-term changes in avian populations. For example, Conroy et al. (2015) used a combination of changes in δ 15 N and δ 13 C values with pollen concentrations in lakes from the Galápagos Islands to make inferences about marine productivity and long-term changes in seabird population size. The authors described a positive relationship between increasing red-footed booby (Sula sula) numbers and red mangrove (Rhizophora mangle), the dominant plant species in the study region. Conroy et al. (2015) found that the seabirds provided the nutrients needed for an increase in the mangrove abundance. The mangrove pollen, alongside δ 15 N and δ 13 C measurements, provided a more robust understanding of the booby colony trends.
There are also more direct palynological proxies of seabird inputs. Spores of Sporormiella, an obligate coprophilous fungus, have been used to track mammal populations over long periods (Davis, 1987). However, these fungal spores could also be used to show long-term changes in the population size of endemic avian populations (Wood et al., 2011). Fungal hyphae are another example of a palynological proxy for seabird presence, as hyphae are nearly exclusively terrestrial in origin and are not typically observed in aquatic environments. However, in a study examining the effects of Leach's storm-petrel dynamics on the terrestrial environment, the authors linked the presence of fungal hyphae in the sediments to the presence of stormpetrels (Duda et al., 2020b). This relationship exists because storm-petrels are burrow-nesting seabirds that disturb soils and increase the availability of terrestrial components, which are then introduced into a body of water during runoff events.

Pigments
As birds introduce nutrients to an ecosystem, they typically also increase the overall primary production and alter the algal communities in any affected water column. However, unlike siliceous diatom valves and chrysophyte cysts, most soft-bodied algae do not preserve well in sediments. Therefore, a more complete picture of aquatic phototroph community composition can be inferred from directly measuring fossil pigments preserved in sediments (Leavitt and Hodgson, 2001) or using inferred pigment concentrations derived from visible and near-infrared spectroscopy (VNIRS). A commonly used pigment is chlorophyll-a, which is generally considered a robust indicator of whole-lake production (Michelutti et al., 2010a;Keatley et al., 2011;Chen et al., 2013;Stewart et al., 2015Stewart et al., , 2019Michelutti and Smol, 2016;Duda et al., 2020b,c). Other pigments related to diatoms (e.g., diatoxanthin) and cyanobacteria (e.g., echinenone, myxoxanthophyll) have also been linked to seabird enrichment in the Arctic (Keatley et al., 2011). Similarly, declines in the concentrations of aphanizophyll and canthaxanthin, pigments related to nitrogen-fixing cyanobacteria, have been correlated to a supplemented nitrogen supply related to burgeoning geese populations (MacDonald et al., 2015). In other cases, variations in the full range of reflectance as identified by VNIRS have been used to infer past seabird population sizes Xu et al., 2012Xu et al., , 2019Wu et al., 2017aWu et al., ,b, 2018.

Supporting and Novel Proxies
Numerous other proxies are generally not considered effective to reconstruct colony trends when used independently; however, when coupled with other proxies, they can provide compelling collaborative evidence for changing bird populations. For example, because bird inputs typically increase primary production and therefore the organic matter in the ecosystem, the loss-on-ignition (LOI, literally the amount of combustible material) has been used to identify changes in organic content of the sediments (Brooks et al., 2012;Luoto et al., 2014Luoto et al., , 2015MacDonald et al., 2015;Gao et al., 2018b;Yang et al., 2019). Similarly, magnetic susceptibility, or changes in the magnetism of sediments that are correlated to changes in particle size, can be associated with increased erosion caused by birds disturbing shorelines and catchments (Luoto et al., 2014(Luoto et al., , 2015Nie et al., 2015). Finally, n-alkanoic acids (often as n-alkyls) are lipid biomarkers that originate from the tissues of organic matter. As the carbon chain length differs depending on the source (Chen et al., 2020a), it is possible to associate the known influences of a bird colony on the environment to changes in n-alkyl variety and abundance. Often, short-chain n-alkyls (e.g., C 18 , C 24 ), which are proxies of aquatic moss and algae, have been used to make inferences on changes in colony size as bird droppings provide the nutrients for moss and algae proliferation Huang et al., 2010;Hu et al., 2013).
New proxies that strengthen overall colony size reconstructions and associated inferences are continually established. For example, Zheng et al. (2015) used differences in stable mercury isotope fractionation ( 199 Hg/ 201 Hg) to identify bioaccumulation pathways in sediments that accumulated inputs from penguins and seals. Chen et al. (2020b) identified that penguins increase the amount of organic and acid volatile S relative to sulfate and pyrite S. Lastly, Nie et al. (2014a) linked temporal trends in the concentration of rare earth elements with penguin populations.

CASE STUDIES: APPLICATIONS OF SEDIMENT ARCHIVES IN AVIAN CONSERVATION BIOLOGY
Paleoecological reconstructions of bird colonies provide population managers with new insights regarding changes in bird colonies over timescales vastly predating anything available using modern monitoring methods. These reconstructions offer an opportunity to address a range of questions that have the potential to inform conservation and management strategies for populations of concern.
How Dynamic Are Populations on Long-Term Timescales?
As many bird species are long-lived and highly mobile, they have the potential to form metapopulations (i.e., networks of local populations that are connected by dispersal of individuals; Esler, 2000). Metapopulations are vital in sustaining viable and effective population sizes when faced with local and regional stressors as well as taking advantage of favorable conditions when they emerge. Local stressors include changing habitat availability and influxes of predators, while regional stressors include a variety of forces, including changes in climate and ocean currents. As described below, paleolimnological approaches can be used to assess the effects of local and regional influences on bird colonies on the scale of decades to millennia; leading to an understanding of metapopulation dynamics.
Recent paleolimnological work on Antarctic penguin colonies has revealed the importance of natural phenomena relating to population size. At a local scale, Gao et al. (2018c) used ∼1,600-years of changes in δ 15 N and δ 13 C values, as well as geochemistry from lake sediment cores, to show that the Adélie penguin colonies in the Vestfold Hills of East Antarctica were largely controlled by krill abundance, which itself is correlated to trends in solar minima and maxima as inferred from a nearby ice core. This relationship occurs because increases in total solar irradiance (TSI) results in increases in energy flux to the oceanic ecosystem and thus increases phytoplankton biomass and primary productivity. This bottom-up increase changes the food web to support a greater number of penguins (Saba et al., 2014).
Regional influences can also drive changes in the bird population. Near the Antarctic Peninsula, Roberts et al. (2017) used a ∼8,500-year paleolimnological record from Ardley Island to identify maxima in penguin populations that were nearly extirpated by eruptions of the Deception Island volcano. Using dated ash (i.e., tephra) layers and a variety of proxies to reconstruct the penguin colony trends, the authors concluded that poisonous aerosols (e.g., SO 2 , F) and trace elements (e.g., Cd, As, and Pb) from the initial volcanic eruptions, followed by the smothering effects of ash deposition and prolonged exposure to abrasive glass particles while rearing and feeding, probably caused the site-wide population decline. Notably, the penguin population recovered at this site after ∼400 to 800 years (Roberts et al., 2017).
For some bird populations, a combination of local and regional drivers can cause dramatic changes in colony size. On islands in the Northwest Atlantic Ocean, paleolimnological reconstructions of Leach's storm-petrel colonies have suggested intercolonial movement as an important behavior favoring long-term persistence (Duda et al., 2020a,c). Using several sediment cores and multiple proxies (Supplementary Table 1) from lakes draining the two major storm-petrel colonies on Baccalieu Island (Newfoundland and Labrador, Canada) and Grand Colombier Island (St. Pierre and Miquelon, France), the authors reconstructed up to ∼5,800 years of population dynamics. Although these two sites are only ∼270 km apart, the two colonies peaked in storm-petrel numbers at different times (Figure 4). These paleolimnological data, combined with natural history observations (Pollet et al., 2020) and genetic studies (Paterson and Snyder, 1999;Bicknell et al., 2012), strongly indicate that Leach's storm-petrels require a network of interconnected colonies to maintain the metapopulation. Possible drivers of these fluctuations may include local factors, such as predators and suitable habitat availability, or regional changes, such as shifts in the Labrador Current. Regardless of the driver(s) of population size, long-term data, such as those provided by lake sediments near nesting colonies, are the only known means to explore natural population dynamics over long periods that predate modern monitoring programs. Understanding colony dynamics before human arrival can also address shifting baseline syndrome in conservation biology by providing longer temporal records for contemporary population trends (e.g., Soga and Gaston, 2018). For example, the available colony size estimates for Leach's storm-petrels on Baccalieu Island from 1984 (Sklepkovych and Montevecchi, 1989) and 2013 (Wilhelm et al., 2020) indicated that the colony has been in decline since the 1980s. The paleolimnological data corroborated the decline, but also showed that the colony size in the 1980s was at its highest in nearly two millennia. Therefore, had the original 1980s survey been conducted earlier, in the 1960s for instance, the storm-petrel colony would have appeared to be stable compared to 2013 (Figure 4). Moreover, had the original 1980s survey been conducted in the early 1800s, the colony would have appeared to be increasing in size (Figure 4). In contrast, the Grand Colombier colony appeared to be stable based on four surveys since the 1980s. However, when the timescale was pulled back, it was clear that the colony is now only a fraction (∼16% of the potential carrying capacity) of what was possible before permanent European settlement in the 19th century (Figure 4; Duda et al., 2020a).
These studies provide the opportunity to consider bird population changes before human interferences and in the context of a larger, dynamic population that may be declining as part of natural cycles. Whether a population is currently increasing or decreasing is no longer a sufficient indication of status; identifying when the natural dynamics are disrupted becomes the goal in population monitoring. For example, a colony may be increasing in numbers but only as a result of shortterm environmental conditions that are optimal for the species and thus conservation action may still be warranted. Similarly, a colony may be declining in numbers but may still remain within natural population variability over longer periods. Under these circumstances, a management approach that maintains these ordinary fluctuations and does not focus on the recent declines may be the most effective strategy.

How May Populations Respond to Climate Change?
Because of the multifaceted ecological changes related to recent climate warming, it is difficult to predict how bird populations may respond (Mason et al., 2019). In some regions, warming may result in increased abundance and range expansion (e.g., Frishkoff et al., 2016). By the same token, climate warming may result in changes in vegetation that may have adverse habitats for different taxa (e.g., Kuussaari et al., 2009). How can we predict changes in seabird populations with future warming and thereby provide the most effective conservation and management? Reconstructing how bird populations responded to past known climate events may act as a guide to future changes.
A series of paleoenvironmental reconstructions of Antarctic penguin colonies have identified the colonization dates and periods of "penguin optima" (originally coined by Baroni and Orombelli, 1994), during which penguin populations boomed in response to warmer and wetter climates during the Holocene (Figure 5). The timing of penguin arrival in the Antarctic varies by region but appears to be linked to increased availability of breeding habitat and open-water conditions coinciding with deglaciation (Emslie et al., 2007;Gao et al., 2018a). Post colonization, the period of penguin optimum has also varied depending on the study region (Huang et al., 2009a; Figure 5) but generally occurs during the late Holocene, from ∼2,000 to 4,000 years BP (Huang et al., 2009a;Sun et al., 2013). Based on several paleolimnological studies, Adélie penguins prefer warmer periods with deglaciation and larger ice-free areas for better prey conditions (Huang et al., 2009a;Gao et al., 2018c). However, the potential benefits of climate warming observed in the past may have already been exceeded. In the Ardley Peninsula, south of Tierra del Fuego, the colony of gentoo, Adélie, and chinstrap penguins that had been stable for the last ∼1,300 years appears now to be in decline .
Seabird populations in the tropics have recorded more variable responses to past climates than those occurring in polar regions. In the Xisha Islands (or Hoàng Sa in Viet Nam), the dominant red-footed booby colony reached maximum size during the cooling period of the Little Ice Age (1400 to  1850 CE) (Xu et al., 2016). In contrast to what was observed in the Antarctic, cooling may have had a positive influence on the seabird abundance in the tropics. For example, the expansion of ice cover in polar regions can increase marine primary productivity and food availability (Liu et al., 2002;Shiozaki and Chen, 2013;Xu et al., 2016). Similarly, on Dongdao Island (or Ð o Linh Côn in Viet Nam) from the same Xisha Islands archipelago, the bird population peaked at multiple times over the ∼1350-year record in correspondence to reconstructed sea surface temperature (SST) below 27 • C (Liu et al., 2006c). The authors indicated that a likely reason for the increased seabird abundance was better food supplies. Given the complex interactions of changing climates with entire marine ecosystems, hindcasting may provide some of the most effective predictions of the future of seabird dynamics.
A key benefit of hindcasting population changes is that it covers far longer periods of environmental variation and is largely free from human influences (depending on the location and length of the time series) when compared to modern monitoring programs, which typically incorporate the impacts by human interference as an additional factor, alongside current climate change.

How Have Populations Responded to Human Intrusion?
Before regular wildlife surveying began, the double-crested cormorant population in Lake Ontario was nearly extirpated as a result of dichlorodiphenyltrichloroethane (DDT) exposure (Wires and Cuthbert, 2006). However, following conservation efforts and the ban of DDT, the cormorant population grew from ∼136 breeding pairs in the 1970s to over 100,000 in the early 2000s (Weseloh et al., 2012). Because of this rapid growth and lack of pre-decline monitoring, cormorants were blamed for negatively impacting sport fishing and destroying natural habitats. Government programs aimed at reducing the cormorant populations (e.g., culling, extended hunting programs) were instigated in 2020 (Hobson, 2021). However, without long-term monitoring data, the question remained: What is "normal"? What were population numbers like before significant human interventions? Using a combination of paleolimnological proxies (Supplementary Table 1), it became evident that the modern cormorant population is currently larger than it has been in at least the last ∼150 years (Stewart et al., 2015(Stewart et al., , 2019. With these data, it is possible to estimate the historical size of cormorant populations and also consider the differences in the new environment driven by increasing anthropogenic influence, which together can be used to help set informed and achievable management goals. Although paleolimnological data suggested that the current double-crested cormorant population in Lake Ontario is now more numerous than in the past ∼150 years, the large colonies are principally supported in the new environment by introduced and nuisance warm-water fish species that are growing in numbers due to climate change, such as round goby (Neogobius melanostomus) and alewife (Alosa pseudoharengus) (Madura and Jones, 2016). By preying on these warm-water fish species, it could be argued that cormorants act as a biological control that contributes to maintaining natural biodiversity and supports commercially valuable native cold-water taxa, such as lake trout (Salvelinus namaycush) and lake whitefish (Coregonus clupeaformis) (Johnson et al., 2010(Johnson et al., , 2015. Common eider (Somateria mollissima borealis) is a culturally and economically important northern sea duck, but limited survey data exist at the more remote portions of its range, such as the Canadian Arctic and West Greenland. However, the available information suggested that the West Greenland breeding population had declined by ∼81% over 40 years (Merkel, 2004). Given the lack of regular surveying and numerous potential drivers of the decline (e.g., climate warming, increased hunting pressures), it was challenging to provide informed management advice for a species that supports an important Indigenous harvest. In response, Hargan et al. (2019) used a paleolimnological approach to supplement available records to elucidate the cause of decline (Figure 6). With a reconstructed population record, now up to ∼400 years long, it was clear that the eider declines coincided with a growing human population with higher hunting efficiency due to increased access to firearms and motorized vehicles. From this research, it was apparent that harvest was a strong driver of this population and that this would need to be considered when setting harvest regulations in the future.

CURRENT LIMITATIONS AND POTENTIAL SOLUTIONS
Long-term studies provide very high value in informing policy that can directly address modern ecological issues (Hughes et al., 2017). However, long-term studies are limited due to the substantial time and labor required to sustain the work. As a result, most surveying begins only after a population has experienced a noticeable change, and thus it becomes difficult to consider the natural, long-term dynamics of a species before human interference. Recent developments in paleoecological approaches have made substantial improvements to the application and interpretation of long-term population changes. However, there remain potential advancements for the development of proxies, ecologically important populations to explore, and other challenges to overcome. The following describes potential future directions of paleolimnological reconstructions of populations.
What Is the Relationship Between Paleolimnological Proxies and Absolute Colony Size?
Currently, paleolimnological inferences of colony size represent indices that generally reflect relative trends (i.e., timing and magnitudes of change) and cannot yet be used to definitively assign absolute numbers of birds at a site. This limitation is due to the complex and indirect relationship between bird inputs and the observed changes in sediments. For an estimation of past seabird numbers, several variables, such as the number of birds in the waterbody's catchment, the amount of refuse produced per bird, and post-depositional processes need to be estimated. This estimation is further complicated by the non-linear relationship between bird inputs and paleolimnological proxies [e.g., δ 15 N has a logarithmic relationship with bird numbers, reaching a plateau at the δ 15 N value of the species' guano (Duda et al., 2020a)].
To overcome these limitations, researchers have developed a variety of methods to summarize the paleolimnological data to estimate the trends in colony size. For instance, statistical techniques such as "Q-factor analysis" can provide information on changes in the relative abundance of bird-related bio-elements in the core, and this can be used to infer relative changes in colony size (e.g., Liu et al., 2006c;Huang et al., 2009a). Similarly, multiple proxies have been combined using an algorithmic mean of standardized data to provide an understanding of changes in colony size relative to the data available in the sediment core (Duda et al., 2020a,c). Building on this, Roberts et al. (2017) modelled colony sizes with wide confidence intervals for penguins breeding on Ardley Island in the Antarctic, and Duda et al. (2020a) used known habitat-specific bird densities to estimate the relative difference between modern and past colony size for a colony of Leach's storm-petrels breeding on Grand Colombier Island in the Northwest Atlantic Ocean.
Ground-truthing of paleolimnological reconstructions against well-understood, long-term surveying data is required to further develop the relationships between proxies and bird numbers. Although some studies have corroborated survey data with paleolimnological data (e.g., Duda et al., 2020c), monitoring data are typically very limited on timescales relevant for paleoecological comparisons. However, some bird colonies with a well-established, long-term monitoring program could provide the data required to calibrate the relationship between paleolimnological proxies and seabird numbers. For example, the wandering albatross (Diomedea exulans), breeding on Bird Island of South Georgia Island, has been monitored since 1961 (Wooller et al., 1992;BirdLife International, 2017) and has nearby freshwater bodies that may be viable for paleolimnological reconstruction.

Availability of Sampling Locations
Paleolimnological reconstructions are currently limited to locations where biovector inputs drain and are preserved in a nearby body of water. The body of water must be deep enough so that it does not experience significant mixing of the sediments by strong winds or storms. Further, the body of water cannot be too big, as this would result in the avian inputs being dwarfed compared to those from the environment. However, the paleolimnological techniques described in this review are not limited to bodies of water and can be expanded to peatlands and directly to guano deposits. For example, recent work has effectively used peat deposits, which may occur in locations without bodies of water, to track various seabirds using the North Water polynya in northwest Greenland (Outridge et al., 2016;Davidson et al., 2018) and seabirds nesting on the Falkland Islands (or Islas Malvinas in Argentina) (Groff et al., 2020). Similarly, some species produce geochronological accumulations of wastes based on their roosting behavior. For example, a study of accumulated guano and egested insect remains from a chimney swift (Chaetura pelagica) roost in Kingston, Canada, revealed that DDT application diminished Coleoptera (i.e., beetle) abundance, which altered insect communities and reduced the population size of the insectivorous swift and probably other aerial insectivores (Nocera et al., 2012). Beyond chimney swifts, other bird species that produce substantial accumulations of guano that can be explored using paleoecological techniques, such as gyrfalcon (Falco rusticolus) from Greenland (Burnham et al., 2009) and, potentially, oilbirds (locally known as guácharo; Steatornis caripensis) and Guanay cormorants (Leucocarbo bougainvilliorum).
Numerous regions remain unexplored. Based on our review of 69 studies that explored accessing avian dynamics using sediment cores (Supplementary Table 1), we found that the majority were from Antarctica (n = 36, or 52%) and North America (n = 16, or 23%). The remaining continents are less well studied, with relatively few studies in Asia (n = 10, or 14%), Europe (n = 5, or 7%), Oceania (n = 1, or 1%), South America (n = 1, or 1%), and Africa (n = 0). This geographical spread suggests that using lake sediments for population reconstruction is in its infancy and there remain many more colonies and sites to examine (Otero et al., 2018). The methodology, however, is entirely transferable to other regions, and the only requirement is a sediment archive that receives significant biovector inputs.

Sensitivity to Migration and Movement
A paleolimnological reconstruction from a single site is not sensitive to animal movements. In other words, it is impossible to determine whether a colony has truly declined in numbers or has simply relocated on the basis of a single sediment core. However, this issue can be remedied through the use of multiple sediment cores that span an entire island or multiple sites within a bird-inhabited region. Then, concurrent changes across multiple sites may suggest region-wide shifts in colony size, whereas asynchronous declines may suggest colony relocation.
For example, synchronous shifts in multiple sediments cores have been used to establish both increases (e.g., Stewart et al., 2015Stewart et al., , 2019 and decreases (e.g., Hargan et al., 2019;Duda et al., 2020c) in a bird population size. In contrast, asynchronous changes in sediment cores were used to establish colony redistribution (Duda et al., 2020a). Therefore, whenever possible, we recommend that multiple sediment cores from a study region influenced by the target species be examined.

Mixed Species Signals
Often, seabird colonies consist of several species and thus the paleolimnological signal is a mixture of all seabirds. Although not intrinsically a problem, conservation questions are often directed at a specific species rather than overall seabird numbers. This concern must be addressed when examining past trends, where the assumption is made that the inferred trends are representative of the same species throughout the record. To address this issue, we recommend using markers that may be more specific to the species of interest, such as geochemistry that relates to the organism in question or more broadly using sterol and stanol concentrations. For example, on islands in Lake Ontario, the inputs from double-crested cormorants and ringbilled gulls could be distinguished based on differences in their diets (Hargan et al., 2018;Stewart et al., 2020). Paleolimnological proxies that are more species-specific should continue to be developed to further distinguish the inputs from various species (see "Environmental DNA" and "Bile Acids").

Environmental DNA
Although several paleolimnological proxies have been established to effectively explore animal population dynamics, others are still being developed. For instance, environmental DNA (eDNA), sometimes called sedimentary DNA (sedDNA) or sedimentary ancient DNA (sedaDNA), has recently emerged as a potentially valuable proxy to understand changes with high taxonomic resolution (Domaizon et al., 2017;Capo et al., 2021). The eDNA archived in sediments can be from a variety of biological sources, such as animals, microbes, and vegetation, and understanding species changes in these various taxonomic groups can provide a holistic understanding of how the species affected the ecosystem. Although eDNA has not been used to track changes in avian population size, analogous studies have established the efficacy of the technique. For example, Giguet-Covex et al. (2014) used DNA metabarcoding of lake sediments to identify the presence and absence of mammal and plant DNA over the last ∼6,500 years, for the purpose of identifying periods of overgrazing by agricultural fauna and deforestation by the human population. This study established the ability to track changes in animals and vegetation using sedDNA, and similar approaches can be extended to avian inputs. However, methodological limitations and problems with sedDNA preservation remain important challenges to overcome (Capo et al., 2021). Alternatively, instead of directly tracking animal inputs, it may be possible to use eDNA to identify microbes introduced to the lake ecosystem or indirectly infer changes in avian colony size based on trends in the past microbial community. For instance, it may be possible to identify unique microbes of avian intestines and feces, whose presence in sediments may indicate avian inputs. One such bacterial genus is Catellicoccus, which is ubiquitous in the gut microbiome of numerous avian species, including the colonial thick-billed murre (Uria lomvia) (Góngora et al., 2021). The presence of Catellicoccus has been used to track gull fecal contamination in coastal environments (Sinigalliano et al., 2013) and, by extension, it may be possible to track avian trends if it were present in sediments.
Alternatively, changes in microbial communities identified with eDNA may also be used to track changes in sediments, similar to other biological proxies described earlier (e.g., diatoms, sub-fossil invertebrates). For example, Capo et al. (2019) examined annual variations in microbial community from a 40year sediment record to identify how the aquatic microbiota changed in response to climate. The authors found increases in the relative abundance of Chlorophyta, Dinophyceae, and Cyanobacteria (i.e., taxa correlated to increased phosphorus) during periods of warmer winter temperatures. By extension, similar methods may be used to track changes in avian numbers, which have a strong effect on nutrient concentrations of the aquatic and terrestrial environments.

Bile Acids
Bile acids are an additional emerging paleolimnological proxy, similar to lipid biomolecules, that can be used to track fecal inputs in sediments. Bile is originally produced within the liver of an organism and, when it reaches the intestines after a meal, is biotransformed by the metabolism of microorganisms (Ridlon et al., 2016). The bile acids are then released into the environment via excretion. As bile acid structures vary depending on the animal class (Bull et al., 2002;Hofmann et al., 2010), they may be used to track changes in avian colony size, with potentially a higher resolution for species identification than other lipid molecules or geochemistry.

CONCLUSION
Paleolimnological approaches can provide scientists with the long-term context needed to determine whether a population is in decline or is within the natural range of variability. As a result, resources can be effectively allocated to protect a species with realistic management goals. In our review of 69 studies, we demonstrated the efficacy of sediment archives to add the critical temporal context needed to understand baseline conditions and natural variability of avian populations for the last 10,000 years. Ultimately, by combining paleoecological techniques with modern wildlife monitoring approaches, it will be possible to access temporal perspectives that are otherwise lost to time.

DATA AVAILABILITY STATEMENT
A comprehensive table of the studies examined in this review is provided in the Supplementary Material.

AUTHOR CONTRIBUTIONS
MPD conceived the idea for this review. All authors contributed text and revisions to the manuscript, and approved the final version.

FUNDING
This review was primarily supported by the Natural Sciences and Engineering Research Council of Canada via Discovery grants to the authors as well as a graduate fellowship to MPD.

ACKNOWLEDGMENTS
The art for Figures 1, 2 was created by Kara Perilli (https:// karaperilli.myportfolio.com/) and is reproduced here with her permission. Thank you to Steven Emslie for his comments on