Bioturbation can increase time averaging by downward and upward movements of young and old shells within the entire mixed layer and by accelerating the burial of shells into a sequestration zone (SZ), allowing them to bypass the uppermost taphonomically active zone (TAZ). However, bioturbation can increase shell disintegration concurrently, neutralizing the positive effects of mixing on time averaging. Bioirrigation by oxygenated pore-water promotes carbonate dissolution in the TAZ, and biomixing itself can mill shells weakened by dissolution or microbial maceration, and/or expose them to damage at the sediment–water interface. Here, we fit transition rate matrices to bivalve age–frequency distributions from four sediment cores from the southern California middle shelf (50–75 m) to assess the competing effects of bioturbation on disintegration and time averaging, exploiting a strong gradient in rates of sediment accumulation and bioturbation created by historic wastewater pollution. We find that disintegration covaries positively with mixing at all four sites, in accord with the scenario where bioturbation ultimately fuels carbonate disintegration. Both mixing and disintegration rates decline abruptly at the base of the 20- to 40-cm-thick, age-homogenized surface mixed layer at the three well-bioturbated sites, despite different rates of sediment accumulation. In contrast, mixing and disintegration rates are very low in the upper 25 cm at an effluent site with legacy sediment toxicity, despite recolonization by bioirrigating lucinid bivalves. Assemblages that formed during maximum wastewater emissions vary strongly in time averaging, with millennial scales at the low-sediment accumulation non-effluent sites, a centennial scale at the effluent site where sediment accumulation was high but bioturbation recovered quickly, and a decadal scale at the second high-sedimentation effluent site where bioturbation remained low for decades. Thus, even though disintegration rates covary positively with mixing rates, reducing postmortem shell survival, bioturbation has the net effect of increasing the time averaging of skeletal remains on this warm-temperate siliciclastic shelf.

Understanding how time averaging changes during burial is essential for using Holocene and Anthropocene cores to analyze ecosystem change, given the many ways in which time averaging affects biodiversity measures. Here, we use transition-rate matrices to explore how the extent of time averaging changes downcore as shells transit through a taphonomically complex mixed layer into permanently buried historical layers: this is a null model, without any temporal changes in rates of sedimentation or bioturbation, to contrast with downcore patterns that might be produced by human activity. Assuming stochastic burial and exhumation movements of shells between increments within the mixed layer and stochastic disintegration within increments, we find that almost all combinations of net sedimentation, mixing, and disintegration produce a downcore increase in time averaging (interquartile range [IQR] of shell ages), this trend is typically associated with a decrease in kurtosis and skewness and by a shift from right-skewed to symmetrical age distributions. A downcore increase in time averaging is thus the null expectation wherever bioturbation generates an internally structured mixed layer (i.e., a surface, well-mixed layer is underlain by an incompletely mixed layer): under these conditions, shells are mixed throughout the entire mixed layer at a slower rate than they are buried below it by sedimentation. This downcore trend created by mixing is further amplified by the downcore decline in disintegration rate. We find that transition-rate matrices accurately reproduce the downcore changes in IQR, skewness, and kurtosis observed in bivalve assemblages from the southern California shelf. The right-skewed shell age-frequency distributions typical of surface death assemblages—the focus of most actualistic research—might be fossilized under exceptional conditions of episodic anoxia or sudden burial. However, such right-skewed assemblages will typically not survive transit through the surface mixed layer into subsurface historical layers: they are geologically transient. The deep-time fossil record will be dominated instead by the more time-averaged assemblages with weakly skewed age distributions that form in the lower parts of the mixed layer.

Taphonomy plays diverse roles in paleobiology. These include assessing sample quality relevant to ecologic, biogeographic, and evolutionary questions, diagnosing the roles of various taphonomic agents, processes and circumstances in generating the sedimentary and fossil records, and reconstructing the dynamics of organic recycling over time as a part of Earth history. Major advances over the past 15 years have occurred in understanding (1) the controls on preservation, especially the ecology and biogeochemistry of soft-tissue preservation, and the dominance of biological versus physical agents in the destruction of remains from all major taxonomic groups (plants, invertebrates, vertebrates); (2) scales of spatial and temporal resolution, particularly the relatively minor role of out-of-habitat transport contrasted with the major effects of time-averaging; (3) quantitative compositional fidelity; that is, the degree to which different types of assemblages reflect the species composition and abundance of source faunas and floras; and (4) large-scale variations through time in preservational regimes (megabiases), caused by the evolution of new bodyplans and behavioral capabilities, and by broad-scale changes in climate, tectonics, and geochemistry of Earth surface systems. Paleobiological questions regarding major trends in biodiversity, major extinctions and recoveries, timing of cladogenesis and rates of evolution, and the role of environmental forcing in evolution all entail issues appropriate for taphonomic analysis, and a wide range of strategies are being developed to minimize the impact of sample incompleteness and bias. These include taphonomically robust metrics of paleontologic patterns, gap analysis, equalizing samples via rarefaction, inferences about preservation probability, isotaphonomic comparisons, taphonomic control taxa, and modeling of artificial fossil assemblages based on modern analogues. All of this work is yielding a more quantitative assessment of both the positive and negative aspects of paleobiological samples. Comparisons and syntheses of patterns across major groups and over a wider range of temporal and spatial scales present a challenging and exciting agenda for taphonomy in the coming decades.

Naturally time-averaged accumulations of skeletal remains—death assemblages—provide reliable, albeit temporally coarse, information on the species composition and structure of communities in diverse settings, and their mismatch with local living communities usually signals recent human-driven ecological change. Here, we present the first test of live–dead mismatch as an indicator of human stress using ostracodes. On three islands along a gradient of human population density in the Bahamas, we compared the similarity of living and death assemblages in 10 lakes with relatively low levels of human stress to live–dead similarity in 11 physically comparable lakes subject to industrial, agricultural, or other human activities currently or in the past. We find that live–dead agreement in pristine lakes is consistently excellent, boding well for using death assemblages in modern-day and paleolimnological biodiversity assessments. In most comparison of physically similar paired lakes, sample-level live–dead mismatch in both taxonomic composition and species’ rank abundance is on average significantly greater in the stressed lakes; live–dead agreement is not lower in all samples from stressed lakes, but is more variable. When samples are pooled for lake-level and island-level comparisons, stressed lakes still yield lower live–dead agreement, but the significance of the difference with pristine lakes decreases—species that occur dead-only (or alive-only) in one sample are likely to occur alive (or dead) in other samples. Interisland differences in live–dead agreement are congruent with, but not significantly correlated with, differences in human population density. This situation arises from heterogeneity in the timing and magnitudes of stresses and in the extent of poststress recovery. Live–dead mismatch in ostracode assemblages thus may be a reliable indicator of human impact at the sample level with the potential to be a widely applicable tool for identifying impacted habitats and, perhaps, monitoring the progress of their recovery.

Laboratory experiments on regular echinoids indicate that low water temperatures retard organic decomposition far more effectively than anoxia, and that the primary role of anoxia in the preservation of articulated multi-element calcareous skeletons may be in excluding scavenging organisms. When tumbled at 20 rpm, specimens that were first allowed to decay for two days in warm seawater (30°C) disintegrated more than six times faster than specimens treated at room temperature (23°C) and more than an order of magnitude faster than specimens treated in cool water (11°C). In contrast, the effects of aerobic versus anerobic decay on disintegration rates were insignificant. The longer the period that specimens were allowed to decay before tumbling, the greater the rate at which specimens disintegrated, until a threshold time that appears to mark the decomposition of collagenous ligaments. This required a few days at 30°C, about two weeks at 23°C, and more than 4 weeks at 11°C for Strongylocentrotus. Up until this threshold, coronas disintegrate by a combination of cross-plate fractures and separation along plate sutures; cross-plate fractures thus can be taphonomic in origin and are not necessarily related to predation. Specimens decayed for longer-than-threshold periods of time disintegrate virtually instantaneously upon tumbling by sutural separation only. Undisturbed coronas can remain intact for months, sufficient time for epibiont occupation. Rates of disintegration were documented semi-quantitatively by recognizing seven stages of test disarticulation, and quantitatively by tensometer measures of test strength and toughness. The effects of temperature and oxygen on decay and the existence of a decay threshold in disintegration should apply at least in a qualitative sense to many other animals whose skeletons consist of multiple, collagen-bound elements.

Regular echinoids should still be perceived as taphonomically fragile organisms, but our results suggest the potential for latitudinal as well as bathymetric gradients in the preservation of fossil echinoid faunas. Echinoid preservation under any given set of conditions should also be a function of taxonomic differences in test construction (particularly stereom interlocking along plate sutures) as suggested by previous workers, although our experiments indicate that these effects should only be significant among post-threshold specimens. A survey of regular echinoids from Upper Cretaceous white chalk facies of Britain substantiates the basic experimental patterns, yielding examples of all disarticulation stages and significant taxonomic differences in quality of preservation. A diverse array of borers and encrusters on fossil coronas also corroborates the post-mortem persistence of some tests on mid-latitude seafloors.

The richness (number of species) and evenness (uniformity of species abundances) of death assemblages can differ from corresponding living communities due to processes such as between-habitat transport, environmental condensation, and differential taphonomic destruction. Analysis of 132 single-census live-dead comparisons of benthic molluscs from a variety of soft-bottom marine settings indicates that on average evenness does not differ greatly between live and dead assemblages, regardless of the particular depositional setting or grain size of associated sediment. However, individual death assemblages can deviate quite substantially from their corresponding living assemblages, especially if processed using a fine mesh. In addition, death assemblages collected using sieves with 2 mm mesh or coarser showed consistently and significantly greater evenness than corresponding living assemblages. These results are encouraging for broad-scale assessments of evenness in the fossil record based on the comparison of average values (rather than for individual assemblages) and where trends in evenness are the aim of the study.

Our live-dead comparisons of richness sample-size corrected by rarefaction revealed that death assemblages were on average ~1.45 times richer than the corresponding living assemblages regardless of rarefied size. In 63.6% of death assemblages both dead richness and dead evenness were greater than live, suggesting sufficient time-averaging to catch significant random or directional changes in the living community and/or introduction of individuals from outside the sampled habitat. In 12.9% of collections both dead richness and dead evenness were less than live, suggesting either rapid loss of dead shells so that dead diversity is depressed below the local living community or selective loss of taphonomically vulnerable taxa. In 18.2% of data sets dead richness was elevated but dead evenness was depressed relative to live: these are interpreted to reflect the addition of low-evenness allochthonous material. The remaining 4.5% of data sets had elevated dead evenness but depressed dead richness, suggesting that live and dead in this case may not be closely related.

In seven available time series, temporal volatility in living communities over 6–24 months was considerable but could not account for observed (mostly higher) evenness values in corresponding death assemblages, whose evenness and composition were quite stable in the few examined studies. A densely sampled spatial transect shows that changes in living-assemblage evenness along an environmental gradient were preserved in the corresponding death assemblages, although dead evenness at any location on the gradient was substantially higher than living evenness.

Quantifying the effects of taphonomic processes on species abundances in time-averaged death assemblages (DAs) is pivotal for paleoecological inference. However, fidelity estimates based on conventional “live-dead” comparisons are fundamentally ambiguous: (1) data on living assemblages (LAs) are based on a very short period of sampling and thus do not account for biological variability in the LA, (2) LAs are sampled at the same time as the DA and thus do not necessarily reflect past LAs that contributed to the DA, (3) compositions of LAs and DAs can be autocorrelated owing to shared cohorts, and (4) fidelity estimates are cross-scale estimates because DAs are time-averaged and LAs are not. Some portion of raw (total) live-dead (LD) variation in species composition thus arises from incomplete sampling of LAs and from biological temporal variation among LAs (together = premortem component of LD variation), as contrasted with new variation created by interspecific variation in population turnover and preservation rates and by the time-averaging of skeletal input (together = postmortem component of LD variation). To tackle these problems, we introduce a modified test for homogeneity of multivariate dispersions (HMD) in order to (1) account for temporal autocorrelation in composition between LAs and DAs and (2) decompose total LD compositional variation into premortem and postmortem components, and we use simulations to evaluate the contribution of within-habitat time-averaging on the postmortem component. Applying this approach to 31 marine molluscan data sets, each consisting of spatial replicates of LAs and DAs in a single habitat, we find that total LD variation is driven largely by variation among LAs. However, genuinely postmortem processes have significant effects on composition in 25–65% of data sets (depending on the metric) when the effects of temporal autocorrelation are taken into account using HMD. Had we ignored the effects of autocorrelation, the effects of postmortem processes would have been negligible, inflating the similarity between LAs and DAs. Simulations show that within-habitat time-averaging does not increase total LD variation to a large degree—it increases total LD variation mainly via increasing species richness, and decreases total LD variation by reducing dispersion among DAs. The postmortem component of LD variation thus arises from differential turnover and preservation and multi-habitat time-averaging. Moreover, postmortem processes have less effect on the compositions of DAs in habitats characterized by high variability among LAs than they have on DAs in temporally stable habitats, a previously unrecognized first-order factor in estimating postmortem sources of compositional variation in DAs.

Four basic types of skeletal concentrations are modeled in terms of changes in sedimentation rate alone. The model categorizes fossil concentrations on the relatively objective basis of their bed contacts, and uses this criterion to infer directional shifts in net sedimentation. This radical simplification of accumulation histories, in which hardpart input is held constant, yields a surprisingly powerful model capable of predicting a broad spectrum of taphonomic and paleobiologic phenomena. Type I concentrations grade from less fossiliferous sediments and terminate in omission surfaces; if hardpart supply is held constant, they record a slowdown from positive to zero net sedimentation. Type II concentrations are the same as Type I but terminate in erosion surfaces (slowdown to negative net sedimentation), and Type III and IV concentrations are characterized by basal erosion or omission surfaces, respectively, grade upward into less fossiliferous sediments, and record increases in net sedimentation from negative or zero rates to positive rates. According to the model, samples collected from successive horizons within any of these shell beds will differ in the degree and type of post-mortem bias owing to differing histories of hardpart exposure at the depositional interface. Moreover, because rates of sediment accumulation govern the abundance of hardparts at the depositional interface and thus many of the physical characteristics of the benthic habitat, the dynamics of fossil accumulation have direct consequences for the structure of benthic communities (taphonomic feedback) and for ecologically controlled species morphometry.

The model is highly robust to fluctuations in hardpart input, as judged by its ability to correctly infer modes of formation of concentrations in synthetic stratigraphic sections. In addition, field examples of Type I–IV concentrations show independent evidence of formation during intervals of reduced net sedimentation, and many exhibit trends in taphonomic and paleobiologic features expected from the postulated changes in net sedimentation. The model thus provides a testable working hypothesis for the accumulation of fossil material in a wide range of environments, and should be applicable to concentrations of any taxonomic composition, state of preservation, or geologic age. The power and robustness of this heuristic model in fact argue that fossil-rich and fossil-poor strata provide fundamentally different records of past conditions, and that sedimentation rather than hardpart input is the primary control on the nature of the fossil record.

All else being equal, species with short life spans are expected to be overrepresented in time-averaged death assemblages relative to their standing abundance in the living community, but the magnitude of the distortion of proportional abundance and assemblage evenness has received little attention. Here, information from 30 data sets on the living and dead abundances of marine bivalves in local habitats is combined with a global compilation of bivalve life spans to determine whether bias from mortality rate can explain observed differences in species proportional abundances. Although bivalve maximum life spans range from one to 75 years in these data sets, indicating annual mortality rates of 0.97 to 0.09, the “life span bias” (LB) of a species–the difference between its proportional abundance expected dead and that observed alive–is consistently small in magnitude (average change <2%, maximum about 20%) and random in sign relative to observed discordance (OD = difference between that species' proportional abundance observed dead and that observed alive). The aggregate result for 413 living species occurrences is a significantly positive but weak correlation of OD to LB, with only 10% of variation in OD explained. The model performs better among longer-lived species than among shorter-lived species, probably because longer-lived species conform better to the model assumption that species maintain a constant proportional abundance in the living assemblage over time. Among individual data sets, only seven exhibit significant positive correlations between OD and LB. The model also under-predicts the cases where a death assemblage is dominated by a species that is shorter lived than the dominant species in the living assemblage, indicating that some factor(s) other than or in addition to mortality rate is responsible for OD. We can find no evidence of preservational bias linked to life span, for example through body size. This negative outcome reflects a weak biological relationship between life span and living abundance among bivalves in local habitats, contrary to the terrestrial paradigm, and points toward a simpler model of time-averaged death assemblage formation where higher abundances reflect (under-sampled) past populations. Contrary to long-held expectations, variation in population turnover among species is not a major source of taphonomic bias in time-averaged death assemblages among bivalves and perhaps among other marine groups: bias must arise largely from other factors.

Although only a few studies have explicitly evaluated live-dead agreement of species and community responses to environmental and spatial gradients, paleoecological analyses implicitly assume that death assemblages capture these gradients accurately. We use nine data sets from modern, relatively undisturbed coastal study areas to evaluate how the response of living molluscan assemblages to environmental gradients (water depth and seafloor type; “environmental component” of a gradient) and geographic separation (“spatial component”) is captured by their death assemblages. We find that:

1. Living assemblages vary in composition either in response to environmental gradients alone (consistent with a species-sorting model) or in response to a combination of environmental and spatial gradients (mass-effect model). None of the living assemblages support the neutral model (or the patch-dynamic model), in which variation in species abundance is related to the spatial configuration of stations alone. These findings also support assumptions that mollusk species consistently differ in responses to environmental gradients, and suggest that in the absence of postmortem bias, environmental gradients might be accurately captured by variation in species composition among death assemblages. Death assemblages do in fact respond uniquely to environmental gradients, and show a stronger response when abundances are square-root transformed to downplay the impact of numerically abundant species and increase the effect of rare species.

2. Species' niche positions (position of maximum abundance) along bathymetric and sedimentary gradients in death assemblages show significantly positive rank correlations to species positions in living assemblages in seven of nine data sets (both square-root-transformed and presence-absence data).

3. The proportion of compositional variation explained by environmental gradients in death assemblages is similar to that of counterpart living assemblages. Death assemblages thus show the same ability to capture environmental gradients as do living assemblages. In some instances compositional dissimilarities in death assemblages show higher rank correlation with spatial distances than with environmental gradients, but spatial structure in community composition is mainly driven by spatially structured environmental gradients.

4. Death assemblages correctly identify the dominance of niche metacommunity models in mollusk communities, as revealed by counterpart living assemblages. This analysis of the environmental resolution of death assemblages thus supports fine-scale niche and paleoenvironmental analyses using molluscan fossil records. In spite of taphonomic processes and time-averaging effects that modify community composition, death assemblages largely capture the response of living communities to environmental gradients, partly because of redundancy in community structure that is inherently associated with multispecies assemblages. The molluscan data sets show some degree of redundancy as evidenced by the presence of at least two mutually exclusive subsets of species that replicate the community structure, and simple simulations show that between-sample relationships can be preserved and remain significant even when a large proportion of species is randomly removed from data sets.

This paper tests whether the most common fossil brachiopod, gastropod, and bivalve genera also have intrinsically more durable shells. Commonness was quantified using occurrence frequency of the 450 most frequently occurring genera of these groups in the Paleobiology Database (PBDB). Durability was scored for each taxon on the basis of shell size, thickness, reinforcement (ribs, folds, spines), mineralogy, and microstructural organic content. Contrary to taphonomic expectation, common genera in the PBDB are as likely to be small, thin-shelled, and unreinforced as large, thick-shelled, ribbed, folded, or spiny. In fact, only six of the 30 tests we performed showed a statistically significant relationship between durability and occurrence frequency, and these six tests were equally divided in supporting or contradicting the taphonomic expectation. Thus, for the most commonly occurring genera in these three important groups, taphonomic effects are either neutral with respect to durability or compensated for by other factors (e.g., less durable taxa were more common in the original communities). These results suggest that biological information is retained in the occurrence frequency patterns of our target groups.

Phanerozoic trends in shell and life habit traits linked to postmortem durability were evaluated for the most common fossil brachiopod, gastropod, and bivalve genera in order to test for changes in taphonomic bias. Using the Paleobiology Database, we tabulated occurrence frequencies of genera for 48 intervals of ∼11 Myr duration. The most frequently occurring genera, cumulatively representing 40% of occurrences in each time bin, were scored for intrinsic durability on the basis of shell size, reinforcement (ribs, folds, and spines), life habit, and mineralogy.

Shell durability is positively correlated with the number of genera in a time bin, but durability traits exhibit different temporal patterns across higher taxa, with notable offsets in the timing of changes in these traits. We find no evidence for temporal decreases in durability that would indicate taphonomic bias at the Phanerozoic scale among commonly occurring genera. Also, all three groups show a remarkable stability in mean shell size through the Phanerozoic, an unlikely pattern if strong size-filtering taphonomic megabiases were affecting the fossil record of shelly faunas. Moreover, small shell sizes are attained in the early Paleozoic in brachiopods and in the latest Paleozoic in gastropods but are steady in bivalves; unreinforced shells are common to all groups across the entire Phanerozoic; organophosphatic and aragonitic shells dominate only the oldest and youngest time bins; and microstructures having high organic content are most common in the oldest time bins.

In most cases, the timing of changes in durability-related traits is inconsistent with a late Mesozoic Marine Revolution. The post-Paleozoic increase in mean gastropod reinforcement occurs in the early Triassic, suggesting either an earlier appearance and expansion of durophagous predators or other drivers. Increases in shell durability hypothesized to be the result of increased predation in the late Mesozoic are not evident in the common genera examined here. Infaunal life habit does increase in the late Mesozoic, but it does not become more common than levels already attained during the Paleozoic, and only among bivalves does the elevated late Mesozoic level persist through the Holocene.

These temporal patterns suggest control on the occurrence of durability-related traits by individual evolutionary histories rather than taphonomic megabiases. Our findings do not mean taphonomic biases are absent from the fossil record, but rather that their effects apparently have had little net effect on the relative occurrence of shell traits generally thought to confer higher preservation potential over long time scales.

Paleoecological analyses that test for spatial or temporal variation in diversity must consider not only sampling and preservation bias, but also the effects of temporal scale (i.e., time-averaging). The species-time relationship (STR) describes how species diversity increases with the elapsed time of observation, but its consequences for assessing the effects of time-averaging on diversity of fossil assemblages remain poorly explored. Here, we use a neutral, dispersal-limited model of metacommunity dynamics, with parameters estimated from living assemblages of 31 molluscan data sets, to model the effects of within-habitat time-averaging on the mean composition and multivariate dispersion of assemblages, on diversity at point (single station) and habitat scales (pooled multiple stations), and on beta diversity. We hold sample size constant in STRs to isolate the effects of time-averaging from sampling effects. With increasing within-habitat time-averaging, stochastic switching in the identity of species in living (dispersal-limited) assemblages (1) decreases the proportional abundance of abundant species, reducing the steepness of the rank-abundance distribution, and (2) increases the proportional richness of rare, temporally short-lived species that immigrate from the neutral metacommunity with many rare species. These two effects together (1) can shift the mean composition away from the non-averaged (dispersal-limited) assemblages toward averaged assemblages that are less limited by dispersal, resembling that of the metacommunity; (2) allow the point and habitat diversity to increase toward metacommunity diversity under a given sample size (i.e., the diversity in averaged assemblages is inflated relative to non-averaged assemblages); and (3) reduce beta diversity because species unique to individual stations become shared by other stations when limited by a larger but static species pool. Surprisingly, these scale-dependent changes occur at fixed sample sizes and can become significant after only a few decades or centuries of time-averaging, and are accomplished without invoking ecological succession, environmental changes, or selective postmortem preservation. Time-averaging results in less inflation of diversity at habitat than at point scales; paleoecological studies should thus analyze data at multiple spatial scales, including that of the habitat where multiple bulk samples have been pooled in order to minimize time-averaging effects. The diversity of assemblages that have accumulated over 1000 years at point and habitat scales is expected to be inflated by an average of 2.1 and 1.6, respectively. This degree of inflation is slightly higher than that observed in molluscan death assemblages at these same spatial scales (1.8 and 1.3). Thus, neutral metacommunity models provide useful quantitative constraints on directional but predictable effects of time-averaging. They provide minimal estimates for the rate of increase in diversity with time-averaging because they assume no change in environmental conditions and in the composition of the metacommunity within the window of averaging.

Despite extensive paleoecological analyses of spatial and temporal turnover in species composition, the fidelity with which time-averaged death assemblages capture variation in species composition and diversity partitioning of living communities remains unexplored. Do death assemblages vary in composition between sites to a lesser degree than do living assemblages, as would be predicted from time-averaging? And is the higher number of species observed in death relative to living assemblages reduced with increasing spatial scale? We quantify the preservation of spatial and temporal variation in species composition using 11 regional data sets based on samples of living molluscan communities and their co-occurring time-averaged death assemblages. (1) Compositional dissimilarities among living assemblages (LA) within data sets are significantly positively rank-correlated to dissimilarities among counterpart pairs of death assemblages (DA), demonstrating that pairwise dissimilarity within a study area has a good preservation potential in the fossil record. Dissimilarity indices that downplay the abundance of dominant species return the highest live-dead agreement of variation in species composition. (2) The average variation in species composition (average dissimilarity) is consistently smaller in DAs than in LAs (9 of 11 data sets). This damping of variation might arise from DAs generally having a larger sample size, but the reduction by ∼10–20% mostly persists even in size-standardized analyses (4 to 7 of 11 data sets, depending on metric). Beta diversity expressed by the number of compositionally distinct communities is also significantly reduced in death assemblages in size-standardized analyses (by ∼25%). This damping of variation and reduction in beta diversity is in accord with the loss of temporal resolution expected from time-averaging, without invoking taphonomic bias (from differential preservation or postmortem transportation) or sample-size effects. The loss of temporal resolution should directly reduce temporal variation, and assuming time-for-space substitution owing to random walk within one habitat and/or temporal habitat shifting, it also decreases spatial variation in species composition. (3) DAs are more diverse than LAs at the alpha scale, but the difference is reduced at gamma scales because partitioning of alpha and beta components differs significantly between LAs and DAs. This indicates that the effects of time-averaging are reduced with increasing spatial scale. Thus, overall, time-averaged molluscan DAs do capture variation among samples of the living assemblage, but they tend to damp the magnitude of variation, making them a conservative means of inferring change over time or variation among regions in species composition and diversity. Rates of temporal and spatial species turnover documented in the fossil record are thus expected to be depressed relative to the turnover rates that are predicted by models of community dynamics, which assume higher temporal resolution. Finally, the capture by DAs of underlying variation in the LA implies little variation in the net preservation potential of death assemblages across environments, despite the different taphonomic pathways suggested by taphofacies studies.

Contrary to the geological stereotype of pure-carbonate reef platforms, approximately 50% of shallow shelf area in the Tropics is accumulating siliciclastic and mixed siliciclastic-carbonate sediments. Taphonomic characterization of these settings is thus essential for assessing variation among major facies types within the Tropics, as well as for eventual comparison with higher-latitude settings. Our grab samples and dredge samples of bivalve death assemblages from nine stations in five subtidal habitats in a large marine embayment of Caribbean Panama (Bocas del Toro) provide the first actualistic information on the taphonomic condition of shells in Recent tropical siliciclastic sediments. Focusing on unambiguous damage to bivalve shell interiors, we found that the quality of shell preservation in fine-grained siliciclastics is superb: commonly «10% of specimens are affected by encrustation, boring, edge-rounding fine-scale surface alteration via dissolution, microbioerosion maceration. Pure-carbonate and mixed siliciclastic-carbonate environments containing hard substrata (patch reefs, Halimeda gravelly sand, mud among patch reefs) contain higher numbers of more severely damaged shells (generally >25%) and also higher diversities of fossilizable encrusters and borers. Disarticulation and fragmentation are pervasive across all environments and are probably related to predation rather than to postmortem processes. As in other shallow subtidal study areas, the taxonomic compositions of death assemblages have not been homogenized by postmortem transport but show high spatial fidelity to the distribution of living species. Assemblages from the five sedimentary environments have distinct taphonomic signatures, but the strongest differences are between the two fine-grained, exclusively soft-sediment siliciclastic environments on the one hand and the three environments containing hard substrata on the other. Experimental tests for rates and agents of damage, still in progress, indicate that the most critical environmental variables are exhumation cycles and burial rate. Bivalve death assemblages from Bocas del Toro demonstrate that damage levels in tropical fine-grained siliciclastic environments are much lower than in closely associated reefs and algal sands suggest a less filtered record of biological information.

Bivalve death assemblages from subtidal environments within the tropical Bocas del Toro embayment of Caribbean Panama permit a test of the extent to which levels of damage are determined by the intrinsic nature of shell supply (proportion of epifaunal species, thick shells, calcitic shells, low-organic microstructures), as opposed to the extrinsic postmortem environment that shells experience. Only damage to interior surfaces of shells was used, to ensure that damage was unambiguously postmortem in origin. We find that facies-level differences in patterns of damage (the rank order importance of postmortem encrustation, boring, edge-rounding, fine-scale surface degradation) are overwhelmingly controlled by environmental conditions: in each environment, all subsets of the death assemblage present the same damage profile. The composition of shell supply affects only the intensity of the taphonomic signature (i.e., percentage of shells affected) only in environments containing hard substrata (patch reefs, Halimeda gravelly sand, mud among patch reefs). In these environments, epifauna, whether aragonitic or calcitic and whether thin or thick, exhibit significantly higher damage than co-occurring infauna, probably due to the initial period of seafloor exposure they typically experience after death. Thick shells (>0.5 mm), regardless of life habit or mineralogy, are damaged more frequently than thin shells, probably because of selective colonization by fouling organisms. Calcitic shells show no consistently greater frequency of damage than aragonitic shells high-organic microstructures yield mixed patterns. Taphofacies surveys in such depositional systems could thus be confidently based on any subset of the fauna, including diagenetically residual assemblages of calcitic shells and thick-shelled molds. Further tests are needed to determine whether the higher levels of damage observed on some subsets of shells are a consequence of greater time-averaging (thus lower temporal resolution), greater exposure time, preferential attack (potential bias in relative abundance), or some combination of these. Paleobiologically, however, the implication is that ecological subsets of bivalve assemblages are not isotaphonomic, either in tangible damage or in probable bias, within hard-substrate environments, although they may be within soft-sediment environments. In actualistic studies, targeting broad classes of taxa for comparison across environments maximizes our ability to extrapolate taphonomic guidelines into the fossil record, where life habits, skeletal types shallow subtidal habitats have dramatically different patterns of abundance and deployment.

Age-frequency distributions of dead skeletal material on the landscape or seabed—information on the time that has elapsed since the death of individuals—provide decadal- to millennial-scale perspectives both on the history of production and on the processes that lead to skeletal disintegration and burial. So far, however, models quantifying the dynamics of skeletal loss have assumed that skeletal production is constant during time-averaged accumulation. Here, to improve inferences in conservation paleobiology and historical ecology, we evaluate the joint effects of temporally variable production and skeletal loss on postmortem age-frequency distributions (AFDs) to determine how to detect fluctuations in production over the recent past from AFDs. We show that, relative to the true timing of past production pulses, the modes of AFDs will be shifted to younger age cohorts, causing the true age of past pulses to be underestimated. This shift in the apparent timing of a past pulse in production will be stronger where loss rates are high and/or the rate of decline in production is slow; also, a single pulse coupled with a declining loss rate can, under some circumstances, generate a bimodal distribution. We apply these models to death assemblages of the bivalve Nuculana taphria from the Southern California continental shelf, finding that: (1) an onshore-offshore gradient in time averaging is dominated by a gradient in the timing of production, reflecting the tracking of shallow-water habitats under a sea-level rise, rather than by a gradient in disintegration and sequestration rates, which remain constant with water depth; and (2) loss-corrected model-based estimates of the timing of past production are in good agreement with likely past changes in local production based on an independent sea-level curve.