NISP, Bone Fragmentation, and the Measurement of Taxonomic Abundance

Zooarchaeologists have long recognized that the number of identified specimens (NISP) is dependent on the degree to which bones are fragmented, but attempts are rarely made to control for the effects of fragmentation on NISP. This paper provides insight into those effects by presenting both a formal model of the relationship between NISP and fragmentation and experimental data on that relationship. The experimental data have practical implications regarding the effectiveness of potential measures of bone fragmentation, suggesting that specimen size—which can be determined easily through digital image analysis—is more useful than other variables that have been or might be used as fragmentation measures.     

Evaluating the Economic Significance of Sharks,Skates, and Rays (Elasmobranchs) in Prehistoric Economies

Establishing the significance of elasmobranchs in ancient economies is complicated by a variety of biological and taphonomic problems that hinder NISP, MNI, and weight measures and comparisons with other fish or faunal classes. To help address these problems, we present length and live weight estimates for 53 elasmobranch specimens, along with bone weight to edible meat weight ratios for 38 specimens. We argue that quantification is best accomplished using NISP and MNI measures in conjunction with meat weight estimates obtained from modern specimens. These data illustrate the importance of using multiple quantitative measures (i.e., NISP, MNI, weight, etc.) when evaluating the economic significance of various taxa within faunal assemblages.     

Critical review of the MNI (minimum number of individuals) as a zooarchaeological unit of quantification

The use of minimum number of individuals (MNI) in the analysis of mammal archaeofaunal assemblages has received intense criticism for its derived nature and its dependence of biasing variables. Some authors have argued that similar taxonomic variability as documented by MNI can be achieved with other less biased measuring units, such as number of identified specimens (NISP) and derivates thereof. The present study is the first experimentally controlled test that shows the degree of bias of estimates of MNI and NISP by different analysts. It shows that the margins of error of both measuring units are independent and that MNI can be more accurately estimated than NISP, despite its “derived” inferential nature.     

Problems of Identification and Quantification in Archaeozoological Analysis,Part I: Insights from a Blind Test

In archaeozoology, counts are generally considered as replicable data that accurately represent the initial abundances of elements, individuals, or taxa, although perhaps only at the ordinal scale. However, few studies have tested these assumptions with control data. To improve our knowledge of these issues, we conducted a blind test that involved the analysis of two large experimental samples composed of modern ungulate specimens of known element and taxon. Because the samples differed in level of fragmentation, the blind test provides substantial information on the impact of bone processing on faunal identification and quantification. Our results suggest that Number of Identified SPecimens (NISP) and Minimum Number of Elements (MNE) provide measures of abundance for whole assemblages and for samples limited to non-long bones that are both replicable and accurate at the ratio scale. However, the same metrics generally failed, even at the ordinal level, to predict abundances in analyses restricted to long bones and long bone portions. Given these mixed results, it seems judicious, in agreement with the current majority view among archaeozoologists, to treat faunal tallies as ordinal-level information. Despite issues of reproducibility and the difficulty of aggregating counts with MNE, the blind test also indicates that this measure is more robust at predicting skeletal abundances than NISP. Substantial variations in rates of long bone fragmentation and identification probably explain the poorer performance of NISP in the blind test.     

The goodness of fit and statistical significance of seriation solutions

This paper presents a method for determining the statistical significance and goodness of fit of seriation solutions based on correspondence analysis. The basic assumption is that seriation solutions corresponding to relative chronological sequences should have unimodal distributions of types across assemblages in a sequence. In order to evaluate significance of chronological signal in data, the observed number of modes of a seriation solution based on correspondence analysis is compared to a distribution of number of modes generated by randomization of the original data table. A quantitative measure of the goodness of fit is presented. The method is tested on two data sets: 1) archaeological data with a known chronological ordering 2) non-archaeological data without a chronological patterning. The method successfully detected a significant chronological signal in the former, and failed to detect it in the latter.     

Animal Resource Exploitation among the Virgin River Puebloans in the American Southwest

Abstract

A review of the literature documenting faunal remains recovered from 16 archaeological sites attributed to the Virjin River Puebloan (VRP) tradition (ca. 300 B.C.–A.D. 1200) of the American Southwest was conducted to identify the importance of wild animal resources to the subsistence strategy. The Number of Individual Specimens (NISP) was counted (NISP = 16,386) and sorted by taxa and animal size categories to explore differences in the distributions of faunal remains by archaeological recovery strategy, site type, physiographic location, and chronology. The most significant factor in the differential use of animal resources among the VRPs was physiographic location. Desert lowland sites had greater species diversity. Sites in the intermediate St. George Basin had more large mammals such as deer, and sites located on the higher plateau had more small mammals such as lagomorphs. These results indicate that the subsistence strategy varied across diverse environmental conditions, but animal resources were always an integral part of the VRP adaptation until their abandonment of the region in the early 13th century A.D.     

A new method for the quantitative analysis of cutmark micromorphology

Cutmarks found on the fossilised bones of butchered animals provide direct evidence for the procurement of meat through technological means. As such, they hold some of the oldest available information on cognitive ability and behaviour in human evolution. Here we present a new method that allows a three-dimensional reconstruction of cutmark morphology and the quantification of profile parameters. We have tested this new technique on cutmarks that were experimentally inflicted on a pig rib using a steel knife and an un-retouched flint flake at different angles. The method allows for the cross-sectional shape, the sharpness and depths of the resulting cutmarks to be quantified. The data show that knife mark sections are characterised by a V-shape or √-shape depending on the inclination of the knife. Cutmarks produced with the flint flake were less clearly defined and generally less sharp than those produced by the knife. We discuss the method's potential to provide new information on butchery technique and cognitive abilities developed by the human lineage, from the earliest tool-using hominins through to modern Homo sapiens.     

Controlling the False Discovery Rate: A New Application to Account for Multiple and Dependent Tests in Local Statistics of Spatial Association

Assessing the significance of multiple and dependent comparisons is an important, and often ignored, issue that becomes more critical as the size of data sets increases. If not accounted for, false-positive differences are very likely to be identified. The need to address this issue has led to the development of a myriad of procedures to account for multiple testing. The simplest and most widely used technique is the Bonferroni method, which controls the probability that a true null hypothesis is incorrectly rejected. However, it is a very conservative procedure. As a result, the larger the data set the greater the chances that truly significant differences will be missed. In 1995, a new criterion, the false discovery rate (FDR), was proposed to control the proportion of false declarations of significance among those individual deviations from null hypotheses considered to be significant. It is more powerful than all previously proposed methods. Multiple and dependent comparisons are also fundamental in spatial analysis. As the number of locations increases, assessing the significance of local statistics of spatial association becomes a complex matter. In this article we use empirical and simulated data to evaluate the use of the FDR approach in appraising the occurrence of clusters detected by local indicators of spatial association. Results show a significant gain in identification of meaningful clusters when controlling the FDR, in comparison to more conservative approaches. When no control is adopted, false clusters are likely to be identified. If a conservative approach is used, clusters are only partially identified and true clusters are largely missed. In contrast, when the FDR approach is adopted, clusters are fully identified. Incorporating a correction for spatial dependence to conservative methods improves the results, but not enough to match those obtained by the FDR approach.     

试论十二桥文化的生业方式——以动物考古学研究为中心

通过对十二桥文化成都平原类型和渝东类型典型遗址所出动物骨骼的可鉴定标本、最小个体数和肉量估算的统计来考察其生业方式,结果显示前者的肉食资源以家畜为主,而后者则以狩猎为主。与香炉石文化的生业方式进行比较,并认为在早期文明阶段,同一考古学文化的不同类型其生业方式不尽相同,而不同的考古学文化间的生业方式则可能趋同。     

Cutmark systematics: Analyzing morphometrics and spatial patterning at Palangana

Cutmarks have played an important role in addressing whether our hominid ancestors were hunters or scavengers, describing ritual modification of human bone, defining the origin of metallurgy, and highlighting the diversity of prehistoric butchering behavior. The widespread occurrence of cutmarks on animal bones and their variability allows archaeologists to use this kind of evidence to address a broad range of questions. One goal in examining cutmark diversity is to identify diagnostic cutmarks of prehistoric butchering, processing, and consuming behaviors. Linking cutmarks to specific activities allows us to test fine-grained hypotheses concerning the nature of an archaeological site, and to do this a systematic method for addressing variability in both the appearance and anatomical position of cutmarks is essential. An analysis of caribou bones collected by Lewis Binford from the Palangana site in Alaska is used to build and test a diagnostic cutmark classification using their morphometric and spatial properties. This case study demonstrates that cutmarks can be used to simultaneously address large-scale anthropological questions and reveal intra-site behavioral variability in the archaeological record.     

The fractal geometry of ancient Maya settlement

Ancient Maya settlement patterns exhibit fractal geometry both within communities and across regions. Fractals are self-similar sets of fractional dimension. In this paper, we show how Maya settlement patterns are logically and statistically self-similar. We demonstrate how to measure the fractal dimensions (or Hausdorff–Besicovitch dimensions) of several data sets. We describe nonlinear dynamical processes, such as chaotic and self-organized critical systems, that generate fractal patterns. As an illustration, we show that the fractal dimensions calculated for some Maya settlement patterns are similar to those produced by warfare, supporting recent claims that warfare is a significant factor in Maya settlement patterning.     

Archaeological Quantification of Pottery: The Rims Count Adjusted using the Modulus of Rupture (MR)

Archaeological quantification is a recurrent issue in research about pottery, its typologies and its distribution. We accept the validity of other methods of quantification—sherd count, minimum number of individuals (MNI) or sherd weight—but the methodology that we have proposed for quantification of assemblages of archaeological contexts is the rims count, which has to be transformed into coefficients of reference through a correction using the modulus of rupture (MR). Such correctors are obtained through measuring the percentage of preserved rim of a significant number of sherds of each type and establishing the average of that percentage. This quantification method is easily applicable to all pottery types and it is also statistically reliable. Besides, it can be used in any study in which the gross number of rims is published. Finally, in the case of ceramic transport containers, a second correction can be applied by multiplying the corrected coefficient (number of rims × MR) by its average capacity (AC), another corrector that will allow us to gather statistics according to the litres of transported product. We believe that the rims count (the easiest part to classify) is a fast, relatively easy and very reliable method that needs to be corrected using the MR.     

Mammalian Microfaunal Remains from Khonkho Wankane (Late Formative Period), Mollo Kontu (Middle Horizon Period) and Pukara de Khonkho (Late Intermediate Period) in the Bolivian Altiplano

Paleoecological and taphonomic analysis of the mammalian microfauna from three Bolivian altiplano sites is presented: Khonkho Wankane [Late Formative period; number of identified specimens (NISP) = 1019], Mollo Kontu (Middle Horizon period; NISP = 610) and Pukara de Khonkho (Late Intermediate period; NISP = 165). The following species, all of which are present in this region today, were identified: Akodon boliviensis (Bolivian grass mouse), Andinomys edax (Andean mouse), Auliscomys sublimis (Andean leaf‐eared mouse), Cavia porcellus (guinea pig/cuy), Chinchilla chinchilla (chinchilla), Ctenomys leucodon (white‐toothed tuco‐tuco), Eligmodontia puerulus (Andean gerbil mouse), Galea musteloides (common yellow‐toothed cavy), Lagidium viscacia (vizcacha), Neotomys ebriosus (Andean swamp rat),and Oligoryzomys destructor (destructive pygmy rice rat). Some direct human involvement with the deposition of some of the larger species (Cavia porcellus, Chinchilla chinchilla, Ctenomys leucodon and Lagidium viscacia) was noted, although C. porcellus was identified only from a single element from Mollo Kontu. Most deposition was likely due to natural in‐burrow deaths and other natural attrition of these species. No significant shift in environment during the occupation of these sites was detected. Copyright © 2012 John Wiley & Sons, Ltd.     

Relative frequency of butchering cutmarks produced by obsidian and flint: an experimental approach

An experimental approach was used to determine how stone tool materials affect the frequency of observed cutmarks. Five whitetail deer were butchered by an expert butcher using flint and obsidian tools and the skeletons were prepared and analyzed. The authors found that the deer butchered with obsidian tools generally displayed fewer incidents of cutmarks than their flint-butchered counterparts. Consideration of this variable may prove to be important to archaeologists attempting to answer questions of specialization or when making comparisons of cutmark data between sites.     

Problems of Identification and Quantification in Archaeozoological Analysis,Part II: Presentation of an Alternative Counting Method

Archaeozoologists commonly use Number of Identified SPecimens (NISP) and Minimum Number of Elements (MNE) as measures of anatomical abundances. According to a blind test examining the reproducibility and accuracy of identifications of ungulate remains (Morin et al., Part I, Journal of Archaeological Method and Theory, doi: 10.1007/s10816-016-9300-4), NISP provides estimates of skeletal abundances that are less robust than those based on MNE. However, although results were improved with the latter method, MNE is not free of problems. Here, we show through an analysis of paired NISP-MNE data for 24 classes of elements that MNE is prone to inflate the representation of rare parts (as measured by NISP), a phenomenon more strongly expressed in certain elements than in others. Moreover, some elements show a wide scatter of points, which raises issues of data reproducibility. MNE is also known for being seriously affected by aggregation methods. These fundamental problems severely undermine the value of MNE as a measure of abundance. This article introduces an alternative counting method that avoids many of the weaknesses of MNE. This counting method, called the Number of Distinct Elements (NDE), focuses on the occurrence of pre-determined, invariant landmarks counted on mutually exclusive specimens. Preliminary experimental results suggest that NDE counts are robust predictors of skeletal, and perhaps taxonomic, abundances. Moreover, the NDE approach eliminates the complex and time-consuming task of spreading or drawing specimens to identify fragment overlap. Furthermore, NDE values are additive and easy to calculate. Given these features, the NDE approach represents a compelling alternative to MNE in archaeozoological analysis.     

Exploring the Use of Computer Vision Metrics for Spatial Pattern Comparison

Detection of changes in spatial processes has long been of interest to quantitative geographers seeking to test models, validate theories, and anticipate change. Given the current “data-rich” environment of today, it may be time to reconsider the methodological approaches used for quantifying change in spatial processes. New tools emerging from computer vision research may hold particular potential to make significant advances in quantifying changes in spatial processes. In this article, two comparative indices from computer vision, the structural similarity (SSIM) index, and the complex wavelet structural similarity (CWSSIM) index were examined for their utility in the comparison of real and simulated spatial data sets. Gaussian Markov random fields were simulated and compared with both metrics. A case study into comparison of snow water equivalent spatial patterns over northern Canada was used to explore the properties of these indices on real-world data. CWSSIM was found to be less sensitive than SSIM to changing window dimension. The CWSSIM appears to have significant potential in characterizing change and/or similarity; distinguishing between map pairs that possess subtle structural differences. Further research is required to explore the utility of these approaches for empirical comparison cases of different forms of landscape change and in comparison to human judgments of spatial pattern differences.     

New light on feasting and deposition: exploring accumulation history through taphonomic analysis at later prehistoric middens in Britain

This paper explores the potential of using taphonomic analysis to reconstruct broad-scale variation in patterns of consumption and deposition at six later prehistoric midden sites in the UK. These sites comprise large accumulations of material culture, dominated by faunal and ceramic fragments, presumed to result from feasting events during the Bronze Age–Iron Age transition. New artefact and landscape studies have refined the characterisation of these sites (Tubb 2011a, b; Waddington 2009), but little research has focussed on accumulation history. This research uses simple statistical analyses on a large dataset (NISP >20,000) to compare the prevalence of bone modifications between midden sites. Crucially, significant differences in modification are not interpreted directly. Additional tests are undertaken to assess whether variation in assemblage composition could account for differences in modification. Previous research has demonstrated that certain elements and the remains of certain taxa are more likely to exhibit modification, and consequently, a prevalence of these specimens could account for differences, rather than their resulting from varied depositional treatment (Madgwick and Mulville 2012). Therefore, patterns of modification are only interpreted once compositional differences can be discounted from responsibility. The study is intentionally broad in its focus and assesses whether large-scale inter-site differences in depositional practice can be reconstructed. Clear patterns are observed with some middens accumulating predominantly through rapid, large-scale deposits and others building up through smaller, more gradual deposits and being subject to greater disturbance and bone movement. These findings have implications for our understanding of ritualised consumption and deposition at the Bronze Age–Iron Age transition in Britain.     

Assessing Protocols for Identifying Pacific Island Archaeological Fish Remains: The Contribution of Vertebrae

Three fish bone identification protocols used for determining taxa composition for Pacific island archaeofaunal assemblages are evaluated. The protocols include using the following: (1) the most commonly identified five paired cranial bones and ‘specials’ or unique elements; (2) an expanded number of cranial bones; and (3) the less common inclusion of all vertebrae. Explicit identification and quantification protocols are outlined for systematically incorporating all vertebrae which, predictably, increases the number of identified specimens for an assemblage, thus providing more bones useful for reconstructing live fish biomass (weight and length). Significantly, a range of unique archaeological vertebrae are useful for calculating minimum number of individuals. Using a well‐preserved assemblage from Henderson Island, Pitcairn Group, southeast Polynesia, numbering 6480 fish bones (concentration index = 21 580 m³), we demonstrate differences in rank‐order abundance from three taxon identification protocols. For example, when using all vertebrae grouper (Serranidae) and surgeonfishes (Acanthuridae) are more numerically equivalent than when relying mostly on cranial bones for identification for minimum number of individuals and number of identified specimens. This has important implications for making comparisons between sites or across regions where different identification protocols were used. This pilot study demonstrates that using all vertebrae for taxon identification and quantification, not just unique hypurals (terminal vertebrae) or those from sharks and rays (Elasmobranchii), should be standard practice for identifying a greater number of bones to taxon and thereby providing better reconstructions of prehistoric fishing and subsistence practices in the Pacific. Copyright © 2013 John Wiley & Sons, Ltd.     

Pascal programs for computing taxonomic abundance in samples of fossil mammals

The number of identified specimens (NISP) and the minimum number of individuals (MNI) are the two most widely used measures of taxonomic abundance in samples of fossil bones. Each has disadvantages which can be largely offset when the two are used together and when both are calculated and presented for each skeletal part of each taxon in a sample. However, joint use does not reduce the tedium and likelihood of error involved in MNI calculation, or affect the possibility that MNIs may not be comparable among samples, because different analysts employ different and often vague calculation assumptions. The integrated Pascal computer programs described in this article dramatically cut the time and effort involved in MNI calculation, while simultaneously reducing the likelihood of calculation error. They also ensure that calculation is based on explicit, though flexible assumptions. Finally, they provide an inexpensive and compact way to store and retrieve bone data for manipulation by other programs.     

Quantifying flake scar patterning on cores using 3D recording techniques

A key feature of stone artefact morphology is the arrangement and patterning of negative flake scars left on flakes and cores. Scar patterning is often treated as a rough guide to identifying methods of core preparation and reduction and usually forms a key component of lithic typologies and other systems of analysis. However, scar patterns are often complex and difficult to capture using traditional measurement or classificatory techniques, particularly where flakes are thick and irregular, or where cores are flaked on many sides from a number of platforms. Three-dimensional analysis of flake scars is now more feasible using digital 3D surface scanning technologies, or using three-dimensional measurement tools such as the Microscribe-3DX which is now widely used in biological geometric morphometric studies more generally [D.C. Adams, F.J. Rohlf, D.E. Slice, Ital. J. Zool., 71 (2004) 5–16; F.J. Rohlf, L.F. Marcus, Trends Ecol. Evol. 8 (1993) 129–132. [1 and 34]]. This paper develops a mathematical formula for describing scar patterning using vectors calculated from the start and end points of flake scars recorded in three dimensions.