Vertebrate microfossil bonebeds (VMBs)—localized concentrations of small resilient vertebrate hard parts—are commonly studied to recover otherwise rarely found small-bodied taxa, and to document relative taxonomic abundance and species richness in ancient vertebrate communities. Analyses of taphonomic comparability among VMBs have often found significant differences in size and shape distributions, and thus considered them to be non-isotaphonomic. Such outcomes of “strict” statistical tests of isotaphonomy suggest discouraging limits on the potential for broad, comparative paleoecological reconstruction using VMBs. Yet it is not surprising that sensitive statistical tests highlight variations among VMB sites, especially given the general lack of clarity with regard to the definition of “strict” isotaphonomic comparability. We rigorously sampled and compared six VMB localities representing two distinct paleoenvironments (channel and pond/lake) of the Upper Cretaceous Judith River Formation to evaluate biases related to sampling strategies and depositional context. Few defining distinctions in bioclast size and shape are evident in surface collections, and most site-to-site comparisons of sieved collections are indistinguishable (p≤0.003). These results provide a strong case for taphonomic equivalence among the majority of Judith River VMBs, and bode well for future studies of paleoecology, particularly in relation to investigations of faunal membership and community structure in Late Cretaceous wetland ecosystems. The taphonomic comparability of pond/lake and channel-hosted VMBs in the Judith River Formation is also consistent with a formative model that contends that channel-hosted VMBs were reworked from pre-existing pond/lake assemblages, and thus share taphonomic history.

Although extinction risk has been found to have a consistent negative relationship with geographic range across wide temporal and taxonomic scales, the effect has been difficult to disentangle from factors such as sampling, ecological niche, or clade. In addition, studies of extinction risk have focused on benthic invertebrates with less work on planktic taxa. We employed a global set of 1114 planktic graptolite species from the Ordovician to lower Devonian to analyze the predictive power of species’ traits and abiotic factors on extinction risk, combining general linear models (GLMs), partial least-squares regression (PLSR), and permutation tests. Factors included measures of geographic range, sampling, and graptolite-specific factors such as clade, biofacies affiliation, shallow water tolerance, and age cohorts split at the base of the Katian and Rhuddanian stages.

The percent variance in durations explained varied substantially between taxon subsets from 12% to 45%. Overall commonness, the correlated effects of geographic range and sampling, was the strongest, most consistent factor (12–30% variance explained), with clade and age cohort adding up to 18% and other factors <10%. Surprisingly, geographic range alone contributed little explanatory power (<5%). It is likely that this is a consequence of a nonlinear relationship between geographic range and extinction risk, wherein the largest reductions in extinction risk are gained from moderate expansion of small geographic ranges. Thus, even large differences in range size between graptolite species did not lead to a proportionate difference in extinction risk because of the large average ranges of these species. Finally, we emphasize that the common practice of determining the geographic range of taxa from the union of all occurrences over their duration poses a substantial risk of overestimating the geographic scope of the realized ecological niche and, thus, of further conflating sampling effects on observed duration with the biological effects of range size on extinction risk.

Sauropodomorpha included the largest known terrestrial vertebrates and was the first dinosaur clade to achieve a global distribution. This success is associated with their early adoption of herbivory, and sauropod gigantism has been hypothesized to be a specialization for bulk feeding and obligate high-fiber herbivory. Here, we apply a combination of biomechanical character analysis and comparative phylogenetic methods with the aim of quantifying the evolutionary mechanics of the sauropodomorph feeding apparatus. We test for the role of convergence to common feeding function and divergence toward functional optima across sauropodomorph evolution, quantify the rate of evolution for functional characters, and test for coincident evolutionary rate shifts in craniodental functional characters and body mass. Results identify a functional shift toward increased cranial robustness, increased bite force, and the onset of static occlusion at the base of the Sauropoda, consistent with a shift toward bulk feeding. Trends toward similarity in functional characters are observed in Diplodocoidea and Titanosauriformes. However, diplodocids and titanosaurs retain significant craniodental functional differences, and evidence for convergent adoption of a common “adaptive zone” between them is weak. Modeling of craniodental character and body-mass evolution demonstrates that these functional shifts were not correlated with evolutionary rate shifts. Instead, a significant correlation between body mass and characters related to bite force and cranial robustness suggests a correlated-progression evolutionary mode, with positive-feedback loops between body mass and dietary specializations fueling sauropod gigantism.

Sexual dimorphism is common in many extant animals, but it is difficult to demonstrate in fossil species. Working with material from the Late Cretaceous of the U.S. Coastal Plain, we herein analyze sexual dimorphism in ostracodes from the superfamily Cytheroidea, a group whose extant members have males that are relatively more elongate than females. We digitized outlines of more than 6000 individual ostracode valves or carapaces, extracted size (area) and shape (length-to-height ratio) information, and used finite mixture models to assess hypotheses of sexual dimorphism. Male and female clusters can be discerned in nearly all populations with sufficient data, resulting in estimates of size and shape dimorphism for 142 populations across 106 species; an additional nine samples are interpreted to consist only of females. Dimorphism patterns varied across taxa, especially for body size: males range from 30% larger to 20% smaller than females. Magnitudes of sexual dimorphism are generally stable within species across time and space; we can demonstrate substantial evolutionary changes in dimorphism in only one species, Haplocytheridea renfroensis. Several lines of evidence indicate that patterns of sexual dimorphism in these ostracodes reflect male investment in reproduction, suggesting that this study system has the potential to capture variation in sexual selection through the fossil record.

Among extant crinoids, the feather stars are the most diverse and occupy the greatest bathymetric range, being especially common in reef environments. Feather stars possess a variety of morphological, behavioral and physiological traits that have been hypothesized to be critical to their success, especially in their ability to cope with predation. However, knowledge of their predators is exceptionally scant, consisting primarily of circumstantial evidence of attacks by fishes. In this study the question whether regular echinoids, recently shown to consume stalked crinoids, also consume feather stars is explored. Aquarium observations indicate that regular echinoids find feather stars palatable, including feather stars known to be distasteful to fish, and that regular echinoids can capture and eat live feather stars, including those known to swim. Gut-content analyses of the echinoid Araeosoma fenestratum (Thomson, 1872), which is commonly observed with large populations of the feather star Koehlermetra porrecta (Carpenter, 1888) in video transects from marine canyons off the coast of France, revealed elements of feather stars in the guts of 6 of 13 individuals. The high proportion of crinoid material (up to 90%), and the presence of articulated crinoid skeletal elements in the gut of A. fenestratum, suggest that these echinoids consumed at least some live crinoids, although they may have also ingested some postmortem remains found in the sediment. Additionally, photographic evidence from the northeast Atlantic suggests that another regular echinoid, Cidaris cidaris (Linnaeus, 1758), preys on feather stars. Thus in spite of the broad suite of antipredatory adaptations, feather stars are today subject to predation by regular echinoids and may have been since the Mesozoic, when this group of crinoids first appeared.

Placodontia were a group of marine reptiles that lived in shallow nearshore environments during the Triassic. Based on tooth morphology it has been inferred that they were durophagous, but tooth morphology differs among species: placodontoid placodonts have teeth described as hemispherical, and the teeth of more highly nested taxa within the cyamodontoid placodonts have been described as flat. In contrast, the sister taxon to the placodonts, Palatodonta bleekeri, like many other marine reptiles, has tall pointed teeth for eating soft-bodied prey. The goals of this paper are to quantify these different tooth morphologies and compare tooth shape among taxa and with a functionally “optimal” tooth. To quantify tooth morphology we measured the radius of curvature (RoC) of the occlusal surface by fitting spheres to 3D surface scans or computed microtomographic scans. Large RoCs correspond to flatter teeth, while teeth with smaller RoCs are pointier; positive RoCs have convex occlusal surfaces, and a negative RoC indicates that the occlusal surface of the tooth is concave. We found the placodontoid taxa have teeth with smaller RoCs than more highly nested taxa, and palatine teeth tend to be flatter and closer to the optimal morphology than maxillary teeth. Within one well-nested clade, the placochelyids, the rearmost palatine teeth have a more complex morphology than the predicted optimal tooth, with an overall concave occlusal surface with a small, medial cusp. These findings are in keeping with the hypothesis that placodonts were specialized durophagous predators with teeth modified to break hard prey items while resisting tooth failure.

Numerical and taxonomic resolution of compositional data sets affects investigators’ abilities to detect and measure relationships between communities and environmental factors. We test whether varying numerical (untransformed, square-root- and fourth-root-transformed relative abundance and presence–absence data) and taxonomic (species, genera, families) resolutions reveals different insights into early to middle Miocene molluscan communities along bathymetric and salinity gradients. The marine subtidal has a more even species-abundance distribution, a higher number of rare species, and higher species:family and species:genus ratios than the three habitats—marine and estuarine intertidal, estuarine subtidal—with higher fluctuations in salinity and other physical parameters. Taxonomic aggregation and numerical transformation of data result in very different ordinations, although all habitats differ significantly from one another at all taxonomic and numerical levels. Rank correlations between species-level and higher-taxon, among-sample dissimilarities are very high for proportional abundance and decrease strongly with increasing numerical transformation, most notably in the two intertidal habitats. The proportion of variation explained by depth is highest for family-level data, decreases gradually with numerical transformation, and is higher in marine than in estuarine habitats. The proportion of variation explained by salinity is highest for species-level data, increases gradually with numerical transformation, and is higher in subtidal than in intertidal habitats. Therefore, there is no single best numerical and taxonomic resolution for the discrimination of communities along environmental gradients: the “best” resolution depends on the environmental factor considered and the nature of community response to it. Different numerical and taxonomic transformations capture unique aspects of metacommunity assembly along environmental gradients that are not detectable at a single level of resolution. We suggest that simultaneous analyses of community gradients at multiple taxonomic and numerical resolutions provide novel insights into processes responsible for spatial and temporal community stability.

The late Pleistocene megafaunal extinctions may have been the first extinctions directly related to human activity, but in North America the close temporal proximity of human arrival and the Younger Dryas climate event has hindered efforts to identify the ultimate extinction cause. Previous work evaluating the roles of climate change and human activity in the North American megafaunal extinction has been stymied by a reliance on geographic binning, yielding contradictory results among researchers. We used a fine-scale geospatial approach in combination with 95 megafaunal last-appearance and 75 human first-appearance radiocarbon dates to evaluate the North American megafaunal extinction. We used kriging to create interpolated first- and last-appearance surfaces from calibrated radiocarbon dates in combination with their geographic autocorrelation. We found substantial evidence for overlap between megafaunal and human populations in many but not all areas, in some cases exceeding 3000 years of predicted overlap. We also found that overlap was highly regional: megafauna had last appearances in Alaska before humans first appeared, but did not have last appearances in the Great Lakes region until several thousand years after the first recorded human appearances. Overlap in the Great Lakes region exceeds uncertainty in radiocarbon measurements or methodological uncertainty and would be even greater with sampling-derived confidence intervals. The kriged maps of last megafaunal occurrence are consistent with climate as a primary driver in some areas, but we cannot eliminate human influence from all regions. The late Pleistocene megafaunal extinction was highly variable in timing and duration of human overlap across the continent, and future analyses should take these regional trends into account.

Patterns of turnover in the mammalian fossil record have been interpreted as reflecting “pulses” of originations and extinctions hypothesized to be driven by climate change. However, criteria for determining what constitutes a meaningful pulse have been idiosyncratic, and investigations of turnover patterns in mammals have yielded mixed results.

This study presents simple simulations of fossil records in which origination and extinction probabilities for each lineage are held constant. Nonetheless, the total number of turnover events per time bin varies stochastically, producing statistical “noise.” Various simulation and analytical assumptions are examined to determine their impact on the type I error rate (i.e., how often “pulses” are detected in a purely stochastic process).

Results suggest that simple analytical parameters (length of time bins and turnover-pulse criterion) have predictable and straightforward effects on false-positive rates. Furthermore, “pulses” of turnover of a magnitude similar to that observed in the terrestrial mammalian fossil record may be quite common under realistic analytical conditions.

The null turnover model offers a practical way to evaluate the significance of observed turnover events in future empirical studies of the fossil record. In evaluating the significance of a “pulse” of fossil origination or extinction events, analytical parameters can be explored using this null model to determine the approximate type I error rate for a set of parameters. Because false-positive rates are shown to be quite high, functional trait-based approaches may offer more reliable indicators of the impact of climate change on turnover dynamics.

Studying the deep-time origins of macroecological phenomena can help us to understand their long-term drivers. Given the considerable spatiotemporal bias of the terrestrial fossil record, it behooves us to understand how much biological information is lost. The aim of this study is to establish whether latitudinal diversity gradients are detectable in a biased terrestrial fossil record. I develop a simulated fossilization approach, weighting the probability of terrestrial mammal species appearing in the fossil record based on body size and geographic-range size; larger species with larger range sizes are more likely to enter the fossil record. I create simulated fossil localities from the modern North American mammal record. I vary the percentage of species successfully fossilized and estimate the magnitude of the latitudinal diversity gradient (slope of the richness gradient and degree of species turnover). I find that estimates of the latitudinal diversity gradient are sensitive to the loss of species with small body size and geographic-range sizes. In some cases, simulated fossil-record bias completely obliterates evidence of declining richness with latitude, a phenomenon that is not ameliorated by the application of nonparametric richness estimation. However, if the rate of preservation is medium (50% of species) to high (75% of species), the magnitude of the latitudinal diversity gradient can be reliably estimated. Similarly, changes in the diversity gradient estimates are largely explained by differences in the diversity–climate relationship among iterations, suggesting that these relationships may be measurable in the fossil record.

A central premise of bioclimatic envelope modeling is the assumption of niche conservatism. Whereas such assumptions are testable in modern populations, it is unclear whether niche conservatism holds over deeper time spans and over very large geographic ranges. Hyaenids occupied a diversity of ecological niches over time and space, and until the end-Pleistocene they occurred in Europe and most of Asia, with Asian populations of Crocuta suggested as being genetically distinct from their closest living relatives. Further, little is known regarding whether and how the dietary ecology of extinct populations of Crocuta differed from those of their extant African counterparts. Here, we use a multiproxy approach to assess an assumption of conserved dietary ecology in late Pleistocene extant spotted hyenas via finite element analysis, dental microwear texture analysis, and a novel dental macrowear method (i.e., whether teeth are minimally, moderately, or extremely worn, as defined by degree of dentin exposure) proposed here. Results from finite element simulations of the masticatory apparatus of Chinese and African Crocuta demonstrate lower skull stiffness and higher stress in the orbital region of the former when biting with carnassial teeth, suggesting that Chinese Crocuta could not process prey with the same degree of efficiency as extant Crocuta crocuta. Dental microwear texture data further support this interpretation, as Chinese Crocuta have intermediate and indistinguishable complexity values (indicative of hard-object feeding) between the extant African lion (Panthera leo) and extant hyenas (C. crocuta, Hyaena hyaena, and Parahyaena brunnea), being most similar to the omnivorous P. brunnea. The use of dental macrowear to infer dietary behavior may also be possible in extinct taxa, as evinced by dietary correlations between extant African feliforms and dental macrowear assignments. Collectively, this multiproxy analysis suggests that Chinese Crocuta may have exhibited dietary behavior distinct from that of living C. crocuta, and assumptions of niche conservatism may mask significant dietary variation in species broadly distributed in time and space.

The robusticity of the weight-bearing limbs of large terrestrial animals is expected to increase at a more rapid rate than in their smaller relatives. This scaling has been hypothesized to allow large species to maintain stresses in the limb bones that are similar to those seen in smaller ones. Curvilinear scaling has previously been found in mammals and nonavian theropods but has not been demonstrated in birds. In this study, polynomial regressions of leg-bone length and circumference in terrestrial flightless birds were carried out to test for a relationship similar to that seen in nonavian theropods. Flightless birds exhibit curvilinear scaling, with the femora of large taxa becoming thicker relative to length at a greater rate than in smaller taxa. Evidence was found for nonlinear scaling in the leg bones of nonavian theropods. However, unlike in avians, there is also phylogenetic variation between taxonomic groups, with tyrannosaur leg bones in particular scaling differently than other groups. Phylogenetically corrected quadratic regressions and separate analyses of taxonomic groupings found little phylogenetic variation in flightless birds. It is suggested here that the nonlinear scaling seen in avian femora is due to the need to maintain the position of the knee under a more anterior center of mass, thereby restricting femoral length. The femur of nonavian theropods is not so constrained, with greater variability of the linear scaling relationships between clades. Phylogenetic variation in limb-bone scaling may broaden the errors for mass-predictive scaling equations based on limb-bone measurements of nonavian theropods.

The demonstration of sexual dimorphism in the fossil record can provide vital information about the role that sexual selection has played in the evolution of life. However, statistically robust inferences of sexual dimorphism in fossil organisms are exceedingly difficult to establish, owing to issues of sample size, experimental control, and methodology. This is particularly so in the case of dinosaurs, for which sexual dimorphism has been posited in many species, yet quantifiable data are often lacking. This study presents the first statistical investigation of sexual dimorphism across Dinosauria. It revisits prior analyses that purport to find quantitative evidence for sexual dimorphism in nine dinosaur species. After the available morphological data were subjected to a suite of statistical tests (normality and unimodality tests and mixture modeling), no evidence for sexual dimorphism was found in any of the examined taxa, contrary to conventional wisdom. This is not to say that dinosaurs were not sexually dimorphic (phylogenetic inference suggests they may well have been), only that the available evidence precludes its detection. A priori knowledge of the sexes would greatly facilitate the assessment of sexual dimorphism in the fossil record, and it is suggested that unambiguous indicators of sex (e.g., presence of eggs, embryos, medullary bone) be used to this end.

Chitinous arm hooks (onychites) of belemnoid coleoid cephalopods are widely distributed in Mesozoic sediments. Due to their relative abundance and variable morphology compared with the single, bullet-shaped, belemnite rostrum, arm hooks came into the focus of micropaleontologists as a promising index fossil group for the Jurassic–Cretaceous rock record and have been the target of functional, ecological, and phylogenetic interpretations in the past. Based on three well-preserved arm crowns of the Toarcian diplobelid Chondroteuthis wunnenbergi, we analyzed the shape of a total of 87 micro-hooks. The arm crown of Chondroteuthis is unique in having uniserial rather than biserial hooks. The first application of elliptic Fourier shape analysis to the arm weapons of belemnoid coleoids allows for the distinction of four micro-hook morphotypes and the quantification of shape variation within these morphotypes. Based on the best-preserved arm crown, we reconstructed the distribution of morphotypes within the arm crown and along a single arm. Our quantitative data support former observations that smaller hooks were found close to the mouth and at the most distal arm parts, while the largest hooks were found in the central part of the arm crown. Furthermore, we found a distinct arm differentiation, as not every arm was equipped with the same hook morphotype. Here, we report the functional specialization of the belemnoid arm crown for the first time and speculate about the potential function of the four morphotypes based on comparisons with modern cephalopods. Our analyses suggest a highly adapted functional morphology and intra-individual distribution of belemnoid hooks serving distinct purposes mainly during prey capture.

Our knowledge of the diversity, ecology, and phylogeny of Mesozoic birds has increased significantly during recent decades, yet our understanding of their flight competence remains poor. Wing loading (WL) and aspect ratio (AR) are two aerodynamically relevant parameters, as they relate to energy costs of aerial locomotion and flight maneuverability. They can be calculated in living birds (i.e., Neornithes) from body mass (BM), wingspan (B), and lift surface (S L). However, the estimates for extinct birds can be subject to biases from statistical issues, phylogeny, locomotor adaptations, and diagenetic compaction. Here we develop a sequential approach for generating reliable multivariate models that allow estimation of measurements necessary to determine WL and AR in the main clades of non-neornithine Mesozoic birds. The strength of our predictions is supported by the use of those variables that show similar scaling patterns in modern and stem taxa (i.e., non-neornithine birds) and the similarity of our predictions with measurements obtained from fossils preserving wing outlines. In addition, although our WL and AR values are based on estimates (BM, B, and S L) that have an associated error, there is no cumulative error in their calculation, and both parameters show low prediction errors. Therefore, we present the first taxonomically broad, error-calibrated estimation of these two important aerodynamic parameters in non-neornithine birds. Such estimates show that the WL and AR of the non-neornithine birds here analyzed fall within the range of variation of modern birds (i.e., Neornithes). Our results indicate that most modern flight modes (e.g., continuous flapping, flap and gliding, flap and bounding, thermal soaring) were possible for the wide range of non-neornithine avian taxa; we found no evidence for the presence of dynamic soaring among these early birds.

Estimating biodiversity and its variations through geologic time is a notoriously difficult task, due to several taphonomic and methodological effects that make the reconstructed signal potentially distinct from the unknown, original one. Through a simulation approach, we examine the effect of a major, surprisingly still understudied, source of potential disturbance: the effect of time discretization through biochronological construction, which generates spurious coexistences of taxa within discrete time intervals (i.e., biozones), and thus potentially makes continuous- and discrete-time biodiversity curves very different. Focusing on the taxonomic-richness dimension of biodiversity (including estimates of origination and extinction rates), our approach relies on generation of random continuous-time richness curves, which are then time-discretized to estimate the noise generated by this manipulation. A broad spectrum of data-set parameters (including average taxon longevity and biozone duration, total number of taxa, and simulated time interval) is evaluated through sensitivity analysis. We show that the deteriorating effect of time discretization on the richness signal depends highly on such parameters, most particularly on average biozone duration and taxonomic longevity because of their direct relationship with the number of false coexistences generated by time discretization. With several worst-case but realistic parameter combinations (e.g., when relatively short-lived taxa are analyzed in a long-ranging biozone framework), the original and time-discretized richness curves can ultimately show a very weak to zero correlation, making these two time series independent. Based on these simulation results, we propose a simple algorithm allowing the back-transformation of a discrete-time taxonomic-richness data set, as customarily constructed by paleontologists, into a continuous-time data set. We show that the reconstructed richness curve obtained this way fits the original signal much more closely, even when the parameter combination of the original data set is particularly adverse to an effective time-discretized reconstruction.

Sauropod bone histology has provided a great deal of insight into the life history of these enormous animals. However, because of high growth rates, annual growth rings are not common in sauropod long bones, so directly measuring growth rates and determining sexual maturity require alternative measures. Histological ontogenetic stages (HOS) have been established to describe the changes in bone histology through development for basal Macronaria and Diplodocoidea, and subsequently for Titanosauria. Despite this, the current HOS model is not able to discriminate bone tissues in late ontogeny, when sauropods had reached asymptotic size and continued to live into senescence but their long bones became extensively remodeled by secondary osteons and all primary bone was destroyed. Here we establish remodeling stages (RS) to characterize the Haversian bone development through ontogeny in eight sauropod taxa (Apatosaurinae, Giraffatitan brancai, Camarasaurus spp., Dicraeosaurus spp., Ampelosaurus atacis, Phuwiangosaurus sirindhornae, Magyarosaurus dacus, and Alamosaurus sanjuanensis) and find significant correlation of RS with corresponding femur length (CFL) for the studied taxa, with the exception of Dicraeosaurus and Magyarosaurus. Remodeling stages are based on the maximum number of observable generations of crosscutting osteons from the innermost, mid-, and outermost part of the cortex. The correlation with CFL indicates that secondary osteons present an ontogenetic signal that could extend the histological ontogenetic stages. Remodeling stages also provide additional insight into the changes in histology through ontogeny for Sauropoda. This method has the potential to be used in other taxa, such as thyreophorans and many ornithischians, that develop Haversian tissue through development.

Despite more than a century of interest, body-mass estimation in the fossil record remains contentious, particularly when estimating the body mass of taxa outside the size scope of living animals. One estimation approach uses humeral and femoral (stylopodial) circumferences collected from extant (living) terrestrial vertebrates to infer the body masses of extinct tetrapods through scaling models. When applied to very large extinct taxa, extant-based scaling approaches incur obvious methodological extrapolations leading some to suggest that they may overestimate the body masses of large terrestrial vertebrates. Here, I test the implicit assumption of such assertions: that a quadratic model provides a better fit to the combined humeral and femoral circumferences-to-body mass relationship. I then examine the extrapolation potential of these models through a series of subsetting exercises in which lower body-mass sets are used to estimate larger sets. Model fitting recovered greater support for the original linear model, and a nonsignificant second-degree term indicates that the quadratic relationship is statistically linear. Nevertheless, some statistical support was obtained for the quadratic model, and application of the quadratic model to a series of dinosaurs provides lower mass estimates at larger sizes that are more consistent with recent estimates using a minimum convex-hull (MCH) approach. Given this consistency, a quadratic model may be preferred at this time. Still, caution is advised; extrapolations of quadratic functions are unpredictable compared with linear functions. Further research testing the MCH approach (e.g., the use of a universal upscaling factor) may shed light on the linear versus quadratic nature of the relationship between the combined femoral and humeral circumferences and body mass.

Most geometric morphometric studies are underpinned by sets of photographs of specimens. The camera lens distorts the images it takes, and the extent of the distortion will depend on factors such as the make and model of the lens and camera and user-controlled variation such as the zoom of the lens. Any study that uses populations of geometric data digitized from photographs will have shape variation introduced into the data set simply by the photographic process. We illustrate the nature and magnitude of this error using a 30-specimen data set of Recent New Zealand Mactridae (Mollusca: Bivalvia), using only a single camera and camera lens with four different photographic setups. We then illustrate the use of retrodeformation in Adobe Photoshop and test the magnitude of the variation in the data set using multivariate Procrustes analysis of variance. The effect of photographic method on the variance in the data set is significant, systematic, and predictable and, if not accounted for, could lead to misleading results, suggest clustering of specimens in ordinations that has no biological basis, or induce artificial oversplitting of taxa. Recommendations to minimize and quantify distortion include: (1) that studies avoid mixing data sets from different cameras, lenses, or photographic setups; (2) that studies avoid placing specimens or scale bars near the edges of the photographs; (3) that the same camera settings are maintained (as much as practical) for every image in a data set; (4) that care is taken when using full-frame cameras; and (5) that a reference grid is used to correct for or quantify distortion.