Spurious and functional correlates of the isotopic composition of a generalist across a tropical rainforest landscape
© McGlynn et al; licensee BioMed Central Ltd. 2009
Received: 29 July 2009
Accepted: 24 November 2009
Published: 24 November 2009
The isotopic composition of generalist consumers may be expected to vary in space as a consequence of spatial heterogeneity in isotope ratios, the abundance of resources, and competition. We aim to account for the spatial variation in the carbon and nitrogen isotopic composition of a generalized predatory species across a 500 ha. tropical rain forest landscape. We test competing models to account for relative influence of resources and competitors to the carbon and nitrogen isotopic enrichment of gypsy ants (Aphaenogaster araneoides), taking into account site-specific differences in baseline isotope ratios.
We found that 75% of the variance in the fraction of 15N in the tissue of A. araneoides was accounted by one environmental parameter, the concentration of soil phosphorus. After taking into account landscape-scale variation in baseline resources, the most parsimonious model indicated that colony growth and leaf litter biomass accounted for nearly all of the variance in the δ15N discrimination factor, whereas the δ13C discrimination factor was most parsimoniously associated with colony size and the rate of leaf litter decomposition. There was no indication that competitor density or diversity accounted for spatial differences in the isotopic composition of gypsy ants.
Across a 500 ha. landscape, soil phosphorus accounted for spatial variation in baseline nitrogen isotope ratios. The δ15N discrimination factor of a higher order consumer in this food web was structured by bottom-up influences - the quantity and decomposition rate of leaf litter. Stable isotope studies on the trophic biology of consumers may benefit from explicit spatial design to account for edaphic properties that alter the baseline at fine spatial grains.
It is generally thought that stable isotope ratios vary across terrestrial landscapes. The extent of this spatial variation is enigmatic, and not explicitly built into the design of experiments wherein stable isotopes of carbon and nitrogen are employed to infer the trophic biology of consumers . Typically, a landscape may be sampled broadly, with samples pooled together based on the assumption that baseline isotopic ratios will not vary appreciably within one kilometer (e.g., [2, 3]). Nevertheless, to our knowledge there has been no study to test whether baseline stable isotope ratios are consistent across short distances. Moreover, no study has heretofore elucidated the contribution of spatial variation in isotope ratios to the isotopic variation in the tissue of consumers, relative to the effects of biotic factors shaping spatial patterns in isotopic fractionation.
In a tropical rain forest, nearly half of the aboveground primary production falls to the ground as litter, soon to release its carbon and nutrients into the diverse detrital food web [4, 5]. In tropical rain forest leaf litter, a warehouse of biological diversity, little is known about the functional biology of even the most abundant organisms . The processes that govern detrital food web structure, particularly in tropical forests, are poorly known. Caution is required in the use of isotopes to draw inferences about the trophic position of animals across space, as there is great variance introduced by spatial heterogeneity in the baseline isotope ratios and unknown differences in rates of isotopic discrimination .
The leaf litter food web of tropical rainforests is replete with ants, the most ecologically dominant guild of predators in this environment. Most ants in the leaf litter are believed to be generalists in diet, consuming prey from a variety of trophic levels , and little is known about the environmental factors that may shape diet and isotopic discrimination. Investigations into the variation in the trophic ecology of generalists have focused on effects of heterogeneity in competition  as well as variation in the resource base itself .
Tropical forests demonstrate great heterogeneity in the density of arthropods within leaf litter. This heterogeneity tracks limiting nutrients, especially phosphorus [11, 12]. Because animals are known to shift their diet in response to changes in resource availability, we may expect that spatial variation in the resource base may account for variation in prey selection by ants in leaf litter. Consequently, parameters of the resource base - including phosphorus, the biomass of leaf litter and the rate of its decomposition - may affect isotopic composition of higher-level consumers.
On the other hand, competitor density and diversity are more often considered to influence the trophic position of generalist animals (e.g., [13, 14]). Animals with broad diets may shift their feeding to underexploited resources when faced with competition . The effect of competitors on diet should covary with density and diversity along environmental gradients [16, 17]. Tropical rain forest leaf litter demonstrates great heterogeneity in the density and diversity of ants, which compete with one another for energy and nutrients derived from the resource base, the litter itself [18, 19]. Leaf litter ant densities correspond to characteristics of the resource base as well as microclimate and microhabitat [20, 21]. Thus we might expect the density or diversity of competing ants to account for any variation in the isotopic composition of a litter-foraging ant species.
Stable isotopes are often used to infer the relative trophic position of animals and how trophic position varies across a landscape [1, 22]. To establish relative trophic position of a single species across the landscape of a single forest, stable isotope approaches integrate the trophic position of the entire consumer chain, provided that the baseline isotopic signature is measured for each locality . Shifts in the concentration of stable isotope of nitrogen, 15N, have been used to infer that generalist species have variable diets that respond to proximate environmental conditions , though it is unclear how environmental conditions altered the isotope ratio of baseline resources. Interspecific studies using 15N have confirmed that species of ants differ in trophic level, with some ant species occupying trophic positions similar to herbivores [23, 24]. Tillberg et al.  tested for intraspecific δ15N variation in five species of ants. They found remarkably broad variation within a species, even within a single locality, suggesting that life history and small-scale environmental factors may interact to result in differences in the isotopic composition of ants.
In the present study, we evaluate the carbon and nitrogen isotopic composition of a generalist species, the gyspy ant Aphaenogaster araneoides, across the landscape of a single tropical rain forest. We demonstrate environmental parameters that robustly predict the isotope ratio of A.araneoides, and then account for landscape scale variation to test competing models to account for isotope discrimination factors.
Correlation (r) between soil parameters and isotopic composition of A. araneoides.
Correlation (r) with δ15N
Correlation (r) with δ13C
Regression models testing for the most parsimonious associations with the relative trophic position of A. araneoides (δ15N A. araneoides - δ15 Nbaseline).
model (independent factors)
Leaf litter decomposition rate
Log (mass, standing litter)
Total soil P (0-10 cm depth)
Colony growth, leaf litter decomposition rate
Colony growth, log (mass, standing litter)
Colony growth, total soil P
Colony growth, litter-ant density
Colony growth, litter-ant richness
Regression models to test for the most parsimonious associations with (δ13C A. araneoides - δ13Cbaseline).
Model (independent factors)
Leaf litter decomposition rate
Log (mass, standing litter)
Total soil P (0-10 cm depth)
Total soil P, colony size
Total soil P, leaf litter decomposition rate
Leaf litter decomposition rate, colony size
Total soil P, leaf litter decomposition rate, colony size
The variation in δ15N and δ13C in A. araneoides within this 500-ha forest landscape was of equivalent magnitude to that reported from other ant species that were sampled over a much greater geographic extent by Tillberg et al . It was surprising to find that soil phosphorus was such a tight predictor of the δ15N of A. araneoides, considering the proximity of sampling sites and multiple unknown trophic links between soil and this common predator. We found this relationship occurred because the variance in phosphorus throughout the 500 ha. forest was predictive of δ15N of leaf litter, which is the base of the food web for A. araneoides.
Across this forest landscape, the standing stocks and rates of decomposition of the leaf litter were principal determinants of the isotopic composition of gypsy ants. If one assumes that δ15N discrimination indicates relative trophic position and that the trophic shift of 15N is about 3, then gypsy ants range from one to two trophic steps above leaf litter, with a higher trophic position positively associated with the biomass of standing leaf litter. Such assumptions about δ15N discrimination may not withstand continued scrutiny . Regardless, the rate of decomposition of leaf litter at each sample site was the environmental parameter that best accounted for the source of C in the diet of the local gypsy ants. These findings suggest that the basic properties of leaf litter, and the abiotic and biotic factors that govern these properties in tropical rain forests, are essentially regulating the isotopic composition of a generalist predator in the complex environment of the leaf litter food web.
The degree of δ15C discrimination from litter to ants was slight, though the variation that occurred was well accounted by environmental parameters. We found that ants with greater δ13C discrimination factors came from sites with faster leaf litter decomposition and greater mean colony sizes. We can infer that the source of C varies across the fertility gradient, though we will have less solid conclusions about the steps in the food web that can account for this variation, as the fractionation of C in litter food webs is open to speculation. Carbon experiences variable rates of fractionation depending on the nature of the trophic step . As the rate of decomposition is expected to be closely tied to rates of microbial respiration, we may expect that the weighted source of C in the litter food web - from air, litter and soil - may vary along the decomposition gradient. Complicating matters further is the possibility that δ13C varies along a vertical gradient, because at the ground level much of the CO2 assimilated by leaves is of respiratory origin, but in the canopy much of the CO2 fixed in photosynthesis is of atmospheric origin. The C in forest-floor leaf litter is from a mélange of leaves from all vertical strata, and we are not able to create a defensible mixing model to evaluate C source. Nevertheless, because δ15N and δ13C tightly covaried among species of collembola along another decomposition gradient , we might surmise that the differences in δ13C among sites may be attributed to differences in trophic structure above the level of primary consumers. This might explain the association with colony size, as nutritional demands of the colony may change with age .
The present results reinforce the view that stable isotope ratios for a given species - regardless of its diet or other habits - are apt to vary across space. We have found that much of this variance may be attributed to ambient variation in the isotopic composition of the base of the food web. In this circumstance, the discrimination factor of δ15N is governed by P abundance. The mechanism behind the spatial coupling of P and δ15N remains open to speculation. Two non-competing possibilities include differential rates of N fixation and differential physiological rates of N fractionation under conditions of P limitation .
We found that the N stable isotope ratio within gypsy ant tissue varied across the range of 4 per mil, a magnitude on par with similar results for other ant species, from collections spanning much larger geographic areas . After accounting for the baseline isotopic composition at each locality within the forest, we found that dietary breadth was not as broad as suggested by the ant-tissue δ15N values alone. This finding suggests that comparisons of δ15N for sessile organisms (or organisms that do not move large distances) must be constrained to very short distances unless baseline data are collected from each locality.
Phosphorus limitation closely tracked the δ15N of gypsy ants because the δ15N of the base of the food web varied predictably with phosphorus at a fine spatial scale. The δ15N discrimination factor of gypsy ants was most parsimoniously explained by the mass of leaf litter on the forest floor. Spatial variation in δ13C in ants, relative to the resource base, was almost wholly accounted by colony life history and the rate of leaf litter decomposition. In terrestrial environments, the failure to incorporate fine-scale information on the isotopic composition of baseline resources may erroneously identify environmental factors predicting trophic position.
The study was carried out along a well-documented soil-fertility gradient within the old-growth lowland tropical forest at La Selva Biological Station, in N.E. Costa Rica. This forest (50-250 m elevation) receives ca. 4 m of rain annually . More information about La Selva is available at http://www.ots.ac.cr. We sited our studies at seven forest plots of the CARBONO Project plot network [29, 30], chosen to represent the broadest range of ant density, ant nest sizes and soil nutrient concentrations in this landscape, based on previous studies [11, 21]. The plots were 0.5 ha and were located across ca. 500 ha of the old-growth forest landscape.
At each of the seven sampling sites, three colonies of the gypsy ant, Aphaenogaster araneoides, were located by feeding a forager and following her back to her nest. Each colony was selected to be at least 10 m from all other colonies used in the study, to prevent overlap or adjacency of home ranges to ensure independence of sampling, based on known measures of home range area . All colonies were collected in their entirety by excavation. We considered a nest to be fully excavated when we reached a discrete terminal chamber, typically containing early instar brood unattended by workers that had typically already been collected. Any nest that we judged to be not completely collected was excluded and another nest was collected in its stead. Freshly collected colonies were kept alive for at least eight hours inside re-sealable plastic bags before being frozen at -20°C for at least 24 h. We then counted all individuals in each colony to determine colony size (the number of adult workers) and an index of colony growth (the ratio of adult workers and worker pupae). After drying each colony at 60°C to constant mass, we took three replicate subsamples to characterize the colony with stable isotope analysis. Each subsample consisted of a thorax and six legs of a single individual ant selected to be ca. 1 mg to provide sufficient mass for accurate determination of isotope ratios. The gasters were intentionally excluded to prevent inclusion of residual food. For several other ant species, comparisons of stable isotope values between heads and thoraxes have indicated no significant differences between these tissue types .
The isotopic composition of the standing leaf litter at the time and site of each colony collection was used as the baseline of the food web for A. araneoides. At each ant sampling site, three forest-floor fine-litter samples (including leaf, twig, and reproductive materials), each ca. 10 L in volume, were collected from three haphazardly selected nearby locations. These samples were washed to remove soil and then dried at 60°C to constant mass. Each dried sample was ground to powder in a Wiley Mill and homogenized, three subsamples of ca. 2 mg were analyzed for isotopic composition, and the mean of these values was used to characterize each sampling site.
Stable isotope ratios are represented as δ15N and δ13C, representing per mil (‰) proportionality of heavy:light isotopes, relative to a universal standard for each. Samples were weighed into tin capsules in a microbalance to the nearest milligram. Ratios of stable isotopes in the ant tissue were measured with a PDZ Europa 20/-20 isotope ratio mass spectrometer at the UC Davis Stable Isotope Facility. For each sample, values of δ15N and δ13C were calibrated using values from established laboratory standards, run every 12 samples, calibrated against NIST Standard Reference Materials IAEA-N1, IAEA-N2, IAEA-N3, IAEA-CH7, and NBS-22.
We use (δ15NA.araneoides- δ15Nbaseline) as the δ15N discrimination factor of A. araneoides, in which the baseline is local (collection-site) standing leaf litter. The δ13C discrimination factor of A. araneoides was similarly evaluated, using (δ13CA.araneoides- δ13Cbaseline).
For each of the seven sampling sites, litter-dwelling ant density and richness were estimated with an "intensive sampling" protocol . Ten 1 m2 quadrats were established at 10 m intervals along a linear transect. Within each quadrat, all litter-ant nests were collected, sorted to species, and measured for size and growth estimates using the same method applied to A. araneoides. The species richness of the litter-ants was estimated using the mean number of species per quadrat; the density of litter ants was estimated using the mean number of adult litter ants nesting in quadrat.
We estimated the rate of leaf litter decomposition at each of the seven sampling sites with 2 successive 1-yr deployments of sets of litterbags containing a common litter, 10 g of freshly fallen and dried leaves of the pioneer tree Cecropia obtusifolia. Each litterbag was sized ca. 100 cm2 and was constructed with 55 μm vinyl mesh. Three sets of five litterbags each were placed at each site in June 2006 and June 2007. Bags were sampled destructively in the following time sequence: 2 wk, 4 wk, 8 wk, 26 wk, 52 wk, for a total of six time points including t = 0. Using the mean values for each time step per plot, decomposition rate, k, was calculated based on the exponential rate of decay, e-kt, and for each site we averaged the resulting k values from the two measurement series (2006-2007 and 2007-2008).
Total nutrients in the surface soil were determined for each of the seven sampling sites (plots) based on compositing six regularly-spaced soil cores (0-10 cm depth) from each 0.5 ha plot. After being air-dried, sieved (2 mm) and ground, the samples were analyzed at the Institute of Soil Science and Forest Nutrition, University of Göttingen, Germany by HNO3-pressure extraction and ICP-AES (see for details). While indices of available P have been found to be poorly related to ecosystem function at this site (e.g., , both stand structure and ecosystem processes significantly vary across the 3-fold within-forest range in total soil P at La Selva . Total soil P was also previously shown to be a robust predictor of the variation in density of litter arthropods across the La Selva landscape .
To test the hypotheses that resource base, competition and life history are associated with trophic level and C source, we created multiple regression models for all variables with δ15N discrimination factor (δ15NA.araneoides- δ15Nbaseline) and δ13C discrimination factor (δ13CA.araneoides- δ13Cbaseline) as dependent variables. These multiple regression models should be interpreted cautiously given the number of comparisons and degrees of freedom. We used AIC  to select the factor(s) most closely associated with the dependent variable and to evaluate competing models for within-forest variance in the δ15N and δ13C of A. araneoides.
Support came from a Sally Casanova Memorial RSCAAP award, NSF (EAR-0421178 and OISE-0749047), Occidental College, Howard Hughes Medical Institute and a LS-AMP award to CSU Dominguez Hills. Support for the CARBONO plots was provided by the National Science Foundation (DEB-9629245), the Andrew W. Mellon Foundation, and Conservation International's TEAM Initiative. Field and laboratory work was assisted by B. Dixon, G. Lam, and J. Wegner. The manuscript was improved with comments from D.A. Clark, C. Huang, S. Menke, J. Thomlinson, C. Tillberg and two anonymous reviewers.
- Post D: Using stable isotopes to estimate trophic position: models, methods and assumptions. Ecology. 2002, 83 (3): 703-718.View Article
- Pasquaud S, Lobry J, Elie P: Facing the necessity of describing estuarine ecosystems: a review of food web ecology study techniques. Hydrobiologia. 2007, 588 (1): 159-172. 10.1007/s10750-007-0660-3.View Article
- McNabb D, Halaj J, Wise D: Inferring trophic positions of generalist predators and their linkage to the detrital food web in agroecosystems: a stable isotope analysis. Pedobiologia-International Journal of Soil Biology. 2001, 45 (4): 289-297. 10.1078/0031-4056-00087.View Article
- Seastedt T: The Role of Microarthropods in Decomposition and Mineralization Processes. Annual Reviews in Entomology. 1984, 29 (1): 25-46. 10.1146/annurev.en.29.010184.000325.View Article
- Kaspari M, Garcia M, Harms K, Santana M, Wright S, Yavitt J: Multiple nutrients limit litterfall and decomposition in a tropical forest. Ecology Letters. 2008, 11 (1): 35-43.PubMed
- Moore J, Berlow E, Coleman D, Ruiter P, Dong Q, Hastings A, Johnson N, McCann K, Melville K, Morin P: Detritus, trophic dynamics and biodiversity. Ecology Letters. 2004, 7 (7): 584-600. 10.1111/j.1461-0248.2004.00606.x.View Article
- Martinez del Rio C, Wolf N, Carleton SA, Gannes LZ: Isotopic ecology ten years after a call for more laboratory experiments. Biological Reviews. 2009, 84: 91-111. 10.1111/j.1469-185X.2008.00064.x.View Article
- Hölldobler B, Wilson EO: The ants. 1990, Cambridge, Mass.: Harvard University PressView Article
- Polis G, Anderson W, Holt R: Toward an integration of landscape and food web ecology. Annual Review of Ecology and Systematics. 1997, 28: 289-316. 10.1146/annurev.ecolsys.28.1.289.View Article
- Hishi T, Hyodo F, Saitoh S, Takeda H: The feeding habits of collembola along decomposition gradients using stable carbon and nitrogen isotope analyses. Soil Biology & Biochemistry. 2007, 39 (7): 1820-1823. 10.1016/j.soilbio.2007.01.028.View Article
- McGlynn T, Salinas D, Dunn R, Wood T, Lawrence D, Clark D: Phosphorus Limits Tropical Rain Forest Litter Fauna. Biotropica. 2007, 39 (1): 50-53. 10.1111/j.1744-7429.2006.00241.x.View Article
- Townsend A, Asner G, Cleveland C: The biogeochemical heterogeneity of tropical forests. Trends in Ecology & Evolution. 2008, 23 (8): 424-431. 10.1016/j.tree.2008.04.009.View Article
- Polis G, Strong D: Food web complexity and community dynamics. American Naturalist. 1996, 147 (5): 813-846. 10.1086/285880.View Article
- van Veen F, Morris R, Godfray H: Apparent Competition, Quantitative Food Webs, and the Structure of Phytophagous Insect Communities. Annual Review of Entomology. 2006, 51: 187-208. 10.1146/annurev.ento.51.110104.151120.View ArticlePubMed
- Cruz-Rivera E, Hay ME: Can quantity replace quality? Food choice, compensatory feeding, and fitness of marine mesograzers. Ecology. 2000, 81 (1): 201-219.View Article
- Duffy J, Cardinale B, France K, McIntyre P, Thebault E, Loreau M: The functional role of biodiversity in ecosystems: incorporating trophic complexity. Ecology Letters. 2007, 10 (6): 522-538. 10.1111/j.1461-0248.2007.01037.x.View ArticlePubMed
- Kneitel JM: Intermediate-consumer identity and resources alter a food web with omnivory. Journal of Animal Ecology. 2007, 76 (4): 651-659. 10.1111/j.1365-2656.2007.01250.x.View ArticlePubMed
- Levings SC, Franks NR: Patterns of nest dispersion in a tropical ground ant community. Ecology. 1982, 63: 338-344. 10.2307/1938951.View Article
- McGlynn TP: Ants on the move: Resource limitation of a litter-nesting ant community in Costa Rica. Biotropica. 2006, 38: 419-427. 10.1111/j.1744-7429.2006.00153.x.View Article
- Kaspari M: Litter ant patchiness at the 1-m2 scale: disturbance dynamics in three Neotropical forests. Oecologia. 1996, 107: 265-273. 10.1007/BF00327911.View Article
- McGlynn TP, Fawcett RM, Clark DA: Litter biomass and nutrient determinants of ant density, nest size and growth in a Costa Rican tropical wet forest. Biotropica. 2009, 41: 234-240. 10.1111/j.1744-7429.2008.00465.x.View Article
- Thompson D, Bury S, Hobson K, Wassenaar L, Shannon J: Stable isotopes in ecological studies. Oecologia. 2005, 144 (4): 517-519. 10.1007/s00442-005-0171-8.View ArticlePubMed
- Blüthgen N, Gebauer G, Fiedler K: Disentangling a rainforest food web using stable isotopes: dietary diversity in a species-rich ant community. Oecologia. 2003, 137: 426-435. 10.1007/s00442-003-1347-8.View ArticlePubMed
- Davidson DW, Cook SC, Snelling RR, Chua TH: Explaining the abundance of ants in lowland tropical rainforest canopies. Science. 2003, 300: 969-972. 10.1126/science.1082074.View ArticlePubMed
- Tillberg CV, McCarthy DP, Dolezal AG, Suarez AV: Measuring the tropic ecology of ants using stable isotopes. Insect Soc. 2006, 53: 65-69. 10.1007/s00040-005-0836-7.View Article
- McCutchan J, Lewis W, Kendall C, McGrath C: Variation in trophic shift for stable isotope ratios of carbon, nitrogen, and sulfur. Oikos. 2003, 102 (2): 378-390. 10.1034/j.1600-0706.2003.12098.x.View Article
- Clarkson B, Schipper L, Moyersoen B, Silvester W: Foliar 15 N natural abundance indicates phosphorus limitation of bog species. Oecologia. 2005, 144 (4): 550-557. 10.1007/s00442-005-0033-4.View ArticlePubMed
- McDade L, Bawa KS, Hespenheide HA, Hartshorn GS: La Selva: Ecology and Natural History of a Neotropical Rain Forest. 1994, University Of Chicago Press
- Clark D, Clark D: Landscape-scale variation in forest structure and biomass in a tropical rain forest. Forest Ecology and Management. 2000, 137 (1-3): 185-198. 10.1016/S0378-1127(99)00327-8.View Article
- Espeleta J, Clark D: Multi-scale variation in fine-root biomass in a tropical rain forest: a seven-year study. Ecological Monographs. 2007, 77 (3): 377-404. 10.1890/06-1257.1.View Article
- McGlynn TP, Carr RA, Carson JH, Buma J: Frequent nest relocation in the ant Aphaenogaster araneoides: resources, competition, and natural enemies. Oikos. 2004, 106 (3): 611-611. 10.1111/j.0030-1299.2004.13135.x.View Article
- Bestelmeyer BT, Agosti D, Alonso LE, Brandão CRF, Brown WL, Delabie JHC, Silvestre R: Field techniques for the study of ground-dwelling ants: an overview, description and evaluation. Ants: standard methods for measuring and monitoring biodiversity. Edited by: Agosti D, Majer J, Alonso LE, Schultz T. 2000, Washington D.C.: Smithsonian Institution Press, 122-144.
- Veldkamp E, Becker A, Schwendenmann L, Clark D, Schulte-Bisping H: Substantial labile carbon stocks and microbial activity in deeply weathered soils below a tropical wet forest. Global Change Biology. 2003, 9 (8): 1171-1184. 10.1046/j.1365-2486.2003.00656.x.View Article
- Powers J, Treseder K, Lerdau M: Fine roots, arbuscular mycorrhizal hyphae and soil nutrients in four neotropical rain forests: patterns across large geographic distances. New Phytologist. 2005, 165 (3): 913-921. 10.1111/j.1469-8137.2004.01279.x.View ArticlePubMed
- Johnson J, Omland K: Model selection in ecology and evolution. Trends in Ecology & Evolution. 2004, 19 (2): 101-108. 10.1016/j.tree.2003.10.013.View Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.