Search

Modern humans have smaller faces relative to Middle and Late Pleistocene members of the genus Homo. While facial reduction and differences in shape have been shown to increase biting efficiency in Homo sapiens relative to these hominins, facial size reduction has also been said to decrease our ability to resist masticatory loads. This study compares crania of Homo heidelbergensis and H. sapiens with respect to mechanical advantages of masticatory muscles, force production efficiency, strains experienced by the cranium and modes of deformation during simulated biting. Analyses utilize X-ray computed tomography (CT) scan-based 3D models of a recent modern human and two H. heidelbergensis. While having muscles of similar cross-sectional area to H. heidelbergensis, our results confirm that the modern human masticatory system is more efficient at converting muscle forces into bite forces. Thus, it can produce higher bite forces than Broken Hill for equal muscle input forces. This difference is the result of alterations in relative in and out-lever arm lengths associated with well-known differences in midfacial prognathism. Apparently at odds with this increased efficiency is the finding that the modern human cranium deforms more, resulting in greater strain magnitudes than Broken Hill when biting at the equivalent tooth. Hence, the facial reduction that characterizes modern humans may not have evolved as a result of selection for force production efficiency. These findings provide further evidence for a degree of uncoupling between form and function in the masticatory system of modern humans. This may reflect the impact of food preparation technologies. These data also support previous suggestions that differences in bite force production efficiency can be considered a spandrel, primarily driven by the midfacial reduction in H. sapiens that occurred for other reasons. Midfacial reduction plausibly resulted in a number of other significant changes in morphology, such as the development of a chin, which has itself been the subject of debate as to whether or not it represents a mechanical adaptation or a spandrel.

Alveolar bone, together with the underlying trabecular bone, fulfils an important role in providing structural support against masticatory forces. Diseases such as osteoporosis or periodontitis cause alveolar bone resorption which weakens this structural support and is a major cause of tooth loss. However, the functional relationship between alveolar bone remodelling within the molar region and masticatory forces is not well understood. This study investigated this relationship by comparing mammalian species with different diets and functional loading (Felis catus, Cercocebus atys, Homo sapiens, Sus scrofa, Oryctolagus cuniculus, Ovis aries). We performed histomorphometric analyses of trabecular bone morphology (bone volume fraction, trabecular thickness and trabecular spacing) and quantified the variation of bone and tooth root volumes along the tooth row. A principal component analysis and non-parametric MANOVA showed statistically significant differences in trabecular bone morphology between species with contrasting functional loading, but these differences were not seen in sub-adult specimens. Our results support a strong, but complex link between masticatory function and trabecular bone morphology. Further understanding of a potential functional relationship could aid the diagnosis and treatment of mandibular diseases causing alveolar bone resorption, and guide the design and evaluation of dental implants.

Finite element analysis (FEA) is increasingly applied in skeletal biomechanical research in general, and in fossil studies in particular. Underlying such studies is the principle that FEA provides results that approximate reality. This paper provides further understanding of the reliability of FEA by presenting a validation study in which the deformations experienced by a real cadaveric human cranium are compared to those of an FE model of that cranium under equivalent simulated loading. Furthermore, model sensitivity to simplifications in segmentation and material properties is also assessed. Our results show that absolute deformations are not accurately predicted, but the distribution of the regions of relatively high and low strains, and so the modes of global deformation, are reasonably approximated.

Island subspecies of the Arctic fox Vulpes lagopus differ morphologically from the mainland subspecies. In particular, differences in cranial form may reflect varied biomechanical adaptations associated with hunting and feeding behaviours. We tested the hypothesis that the observed cranial differences between two island foxes (living on two North Pacific islands) and those living on the mainland have no impact on biomechanical performance during simulated biting. 3D cranial models of three Arctic fox subspecies were compared based on biomechanical parameters (e.g. local strain and large-scale deformation). Finite elements (FE) analyses were used to simulate equivalent biting loads, and geometric morphometrics was used to compare the modes of deformation among the models. The results showed differences in local strains and modes of global deformation among the three subspecies; the mainland subspecies was particularly distinct from the island subspecies. The representative cranium of the mainland subspecies experienced higher strain than that of both island subspecies during all bites. However, the findings highlight issues that arise when relating biomechanical performance, measured via FE analyses, to the foods consumed rather than to the mechanical properties of the individual’s diet. Additional physical properties data for each prey type are necessary to determine the extent to which the present findings reflect biomechanical adaptations to diet and prey acquisition.

The study of morphological variation in the hominin fossil record has been transformed in recent years by the advent of high resolution 3D imaging combined with improved geometric morphometric (GM) toolkits. In parallel, increasing numbers of studies have applied finite elements analysis (FEA) to the study of skeletal mechanics in fossil and extant hominoid material. While FEA studies of fossils are becoming ever more popular they are constrained by the difficulties of reconstruction and by the uncertainty that inevitably attaches to the estimation of forces and material properties. Adding to these modelling difficulties it is still unclear how FEA analyses should best deal with species variation. Comparative studies of skeletal form and function can be further advanced by applying tools from the GM toolkit to the inputs and outputs of FEA studies. First they facilitate virtual reconstruction of damaged material and can be used to rapidly create 3D models of skeletal structures. Second, GM methods allow variation to be accounted for in FEA by warping models to represent mean and extreme forms of interest. Third, GM methods can be applied to compare FEA outputs the ways in which skeletal elements deform when loaded. Model comparisons are hampered by differences in material properties, forces and size among models but how deformations from FEA are impacted by these parameters is increasingly well understood, allowing them to be taken into account in comparing FEA outputs. In this paper we review recent advances in the application of GM in relation to FEA studies of craniofacial form in hominins, providing examples from our recent work and a critical appraisal of the state of the art. info:eu-repo/semantics/publishedVersion

In recent years finite element analysis (FEA) has emerged as a useful tool for the analysis of skeletal form-function relationships. While this approach has obvious appeal for the study of fossil specimens, such material is often fragmentary with disrupted internal architecture and can contain matrix that leads to errors in accurate segmentation. Here we examine the effects of varying the detail of segmentation and material properties of teeth on the performance of a finite element model of a Macaca fascicularis cranium within a comparative functional framework. Cranial deformations were compared using strain maps to assess differences in strain contours and Procrustes size and shape analyses, from geometric morphometrics, were employed to compare large scale deformations. We show that a macaque model subjected to biting can be made solid, and teeth altered in material properties, with minimal impact on large scale modes of deformation. The models clustered tightly by bite point rather than by modeling simplification approach, and fell out as being distinct from another species. However localized fluctuations in predicted strain magnitudes were recorded with different modeling approaches, particularly over the alveolar region. This study indicates that, while any model simplification should be undertaken with care and attention to its effects, future applications of FEA to fossils with unknown internal architecture may produce reliable results with regard to general modes of deformation, even when detail of internal bone architecture cannot be reliably modeled.

Objectives Changes to soft- and hard-tissue components of the masticatory complex during development can impact functional performance by altering muscle excursion potential, maximum muscle forces, and the efficiency of force transfer to specific bitepoints. Within Macaca fascicularis, older individuals exploit larger, more mechanically resistant food items and more frequently utilize wide-gape jaw postures. We therefore predict that key architectural and biomechanical variables will scale during ontogeny to maximize bite force and gape potential within older, larger-bodied individuals. Materials and methods We analyzed 26 specimens of M. fascicularis, representing a full developmental spectrum. The temporalis, superficial masseter, and deep masseter were dissected to determine muscle mass, fiber length, and physiologic cross-sectional area (PCSA). Lever-arm lengths were also measured for each muscle, alongside the height of the temporomandibular joint (TMJ) and basicranial length. These variables were scaled against two biomechanical variables (jaw length and condyle-molar length) to determine relative developmental changes within these parameters. Results During ontogeny, muscle mass, fiber length, and PCSA scaled with positive allometry relative to jaw length and condyle-molar length within all muscles. TMJ height also scaled with positive allometry, while muscle lever arms scaled with isometry relative to jaw length and with positive allometry (temporalis) or isometry (superficial and deep masseter) relative to condyle-molar length. Conclusion Larger individuals demonstrate adaptations during development towards maximizing gape potential and bite force potential at both an anterior and posterior bitepoint. These data provide anatomical evidence to support field observations of dietary and behavioral differences between juvenile and adult M. fascicularis.

Dielectric measurements provide valuable information about the properties of materials, and could be used to classify and identify the source of objects, in fields such as archaeology. Current methods of identification are all partly destructive, so an innovative electromagnetic method developed by the authors, based on resonant cavity perturbation (RCP), provides an attractive, non-destructive alternative. A problem with traditional RCP is that the changes in frequency and Q-factor vary with the object’s shape; however we overcome this by creating a replica of the object, from a material whose dielectric properties are known. Then, by combining three separate perturbations with orthogonal field directions, due firstly to the object and then to its replica, we eliminate the shape dependency, and thus determine the object’s dielectric constant and loss factor. After developing the theory of this novel DRM technique, we demonstrate the principle using a set of geometric shapes made in both polytetrafluoroethylene (PTFE) and a 3D printed material. Further measurements then enable second-order terms to be included in the model, improving its accuracy. Finally, DRM is shown to be capable of distinguishing two irregularly shaped objects of different materials. Potential applications of DRM include determining the provenance of pottery, glasses and flints, and distinguishing ivory from bone. These would be of interest to customs and environmental agencies, as well as museum curators and archaeologists.

In 2015 an engraved shale pendant was found during excavations at the Early Mesolithic site of Star Carr, UK. Engraved motifs on Mesolithic pendants are extremely rare, with the exception of amber pendants from southern Scandinavia. The artwork on the pendant is the earliest known Mesolithic art in Britain; the 'barbed line' motif is comparable to styles on the Continent, particularly in Denmark. When it was first uncovered the lines were barely visible but using a range of digital imaging techniques it has been possible to examine them in detail and determine the style of engraving as well as the order in which the lines might have been made. In addition, microwear and residue analyses were applied to examine whether the pendant showed signs that it had been strung or worn, and whether the lines had been made more visible through the application of pigments, as has been suggested for some Danish amber pendants. This approach of using multiple scientific and analytical techniques has not been used previously and provides a methodology for the examination of similar artefacts in the future.

Shamanic belief systems represent the first form of religious practice visible within the global archaeological record. Here we report on the earliest known evidence of shamanic costume: modified red deer crania headdresses from the Early Holocene site of Star Carr (c. 11 kya). More than 90% of the examples from prehistoric Europe come from this one site, establishing it as a place of outstanding shamanistic/cosmological significance. Our work, involving a programme of experimental replication, analysis of macroscopic traces, organic residue analysis and 3D image acquisition, metrology and visualisation, represents the first attempt to understand the manufacturing processes used to create these artefacts. The results produced were unexpected-rather than being carefully crafted objects, elements of their production can only be described as expedient.

An accurate, dynamic, functional model of the skull that can be used to predict muscle forces, bite forces, and joint reaction forces would have many uses across a broad range of disciplines. One major issue however with musculoskeletal analyses is that of muscle activation pattern indeterminacy. A very large number of possible muscle force combinations will satisfy a particular functional task. This makes predicting physiological muscle recruitment patterns difficult. Here we describe in detail the process of development of a complex multibody computer model of a primate skull (Macaca fascicularis), that aims to predict muscle recruitment patterns during biting. Using optimisation criteria based on minimisation of muscle stress we predict working to balancing side muscle force ratios, peak bite forces, and joint reaction forces during unilateral biting. Validation of such models is problematic; however we have shown comparable working to balancing muscle activity and TMJ reaction ratios during biting to those observed in vivo and that peak predicted bite forces compare well to published experimental data. To our knowledge the complexity of the musculoskeletal model is greater than any previously reported for a primate. This complexity, when compared to more simple representations provides more nuanced insights into the functioning of masticatory muscles. Thus, we have shown muscle activity to vary throughout individual muscle groups, which enables them to function optimally during specific masticatory tasks. This model will be utilised in future studies into the functioning of the masticatory apparatus.

Biomechanical analyses are commonly conducted to investigate how craniofacial form relates to function, particularly in relation to dietary adaptations. However, in the absence of corresponding muscle activation patterns, incomplete muscle data recorded experimentally for different individuals during different feeding tasks are frequently substituted. This study uses finite element analysis (FEA) to examine the sensitivity of the mechanical response of a Macaca fascicularis cranium to varying muscle activation patterns predicted via multibody dynamic analysis. Relative to the effects of varying bite location, the consequences of simulated variations in muscle activation patterns and of the inclusion/exclusion of whole muscle groups were investigated. The resulting cranial deformations were compared using two approaches; strain maps and geometric morphometric analyses. The results indicate that, with bite force magnitude controlled, the variations among the mechanical responses of the cranium to bite location far outweigh those observed as a consequence of varying muscle activations. However, zygomatic deformation was an exception, with the activation levels of superficial masseter being most influential in this regard. The anterior portion of temporalis deforms the cranial vault, but the remaining muscles have less profound effects. This study for the first time systematically quantifies the sensitivity of an FEA model of a primate skull to widely varying masticatory muscle activations and finds that, with the exception of the zygomatic arch, reasonable variants of muscle loading for a second molar bite have considerably less effect on cranial deformation and the resulting strain map than does varying molar bite point. The implication is that FEA models of biting crania will generally produce acceptable estimates of deformation under load as long as muscle activations and forces are reasonably approximated. In any one FEA study, the biological significance of the error in applied muscle forces is best judged against the magnitude of the effect that is being investigated.

Recent developments in simulating musculoskeletal functioning in the craniofacial complex using multibody dynamic analysis and finite elements analysis enable comprehensive virtual investigations into musculoskeletal form and function. Because the growth of the craniofacial skeleton is strongly influenced by mechanical functioning, these methods have potential in investigating the normal and abnormal development of the skull: loading history during development can be predicted and bony adaptations to these loads simulated. Thus these methods can be used to predict the impact of altered loading or modifications of skull form early in ontogeny on the subsequent development of structures. Combining functional models with geometric morphometric methods (GMM), which are principally concerned with the study of variations of form, offers the opportunity to examine variations in form during development and the covariations between form and factors such as functional performance. Such a combination of functional models and GMM can potentially be applied in many useful ways, for example: to build and modify functional models, to assess the outcomes of remodelling studies by comparing the results with morphological changes during ontogeny, and to compare the outcomes of finite element analyses within a multivariate framework. Studies using these tools can not only investigate the development of the skull but also the mechanical processes and thus to some degree, behaviours underlying the development of variation among extant and fossil skeletal elements. By bringing together these tools from quite different comparative traditions, a novel and potentially powerful framework for simulation and statistical biomechanical analyses of form and function emerges. This paper reviews these recent developments in the context of the evolutionary and functional influences on skull development.

The development of virtual methods for anatomical reconstruction and functional simulation of skeletal structures offers great promise in evolutionary and ontogenetic investigations of form–function relationships. Key developments reviewed here include geometric morphometric methods for the analysis and visualization of variations in form (size and shape), finite element methods for the prediction of mechanical performance of skeletal structures under load and multibody dynamics methods for the simulation and prediction of musculoskeletal function. These techniques are all used in studies of form and function in biology, but only recently have they been combined in novel ways to facilitate biomechanical modelling that takes account of variations in form, can statistically compare performance, and relate performance to form and its covariates. Here we provide several examples that illustrate how these approaches can be combined and we highlight areas that require further investigation and development before we can claim a mature theory and toolkit for a statistical biomechanical framework that unites these methods.

Finite element analysis (FEA) is a modelling technique increasingly used in anatomical studies investigating skeletal form and function. In the case of the cranium this approach has been applied to both living and fossil taxa to (for example) investigate how form relates to function or infer diet or behaviour. However, FE models of complex musculoskeletal structures always rely on simplified representations because it is impossible completely to image and represent every detail of skeletal morphology, variations in material properties and the complexities of loading at all spatial and temporal scales. The effects of necessary simplifications merit investigation. To this end, this study focuses on one aspect, model geometry, which is particularly pertinent to fossil material where taphonomic processes often destroy the finer details of anatomy or in models built from clinical CTs where the resolution is limited and anatomical details are lost. We manipulated the details of a finite element (FE) model of an adult human male cranium and examined the impact on model performance. First, using digital speckle interferometry, we directly measured strains from the infraorbital region and frontal process of the maxilla of the physical cranium under simplified loading conditions, simulating incisor biting. These measured strains were then compared with predicted values from FE models with simplified geometries that included modifications to model resolution, and how cancellous bone and the thin bones of the circum-nasal and maxillary regions were represented. Distributions of regions of relatively high and low principal strains and principal strain vector magnitudes and directions, predicted by the most detailed FE model, are generally similar to those achieved in vitro. Representing cancellous bone as solid cortical bone lowers strain magnitudes substantially but the mode of deformation of the FE model is relatively constant. In contrast, omitting thin plates of bone in the circum-nasal region affects both mode and magnitude of deformation. Our findings provide a useful frame of reference with regard to the effects of simplifications on the performance of FE models of the cranium and call for caution in the interpretation and comparison of FEA results.

During development the masticatory apparatus undergoes significant changes; including changes in muscle size and orientation, and in dental form due to replacement and wear. As body size increases during development, absolute masticatory forces are predicted to increase. However the effect of simultaneous changes to muscle lever mechanics and dental topography upon masticatory performance is currently unknown. Cercocebus atys (the sooty mangabey) provides an ideal opportunity to examine this. Over 50% of the diet among all post-weaning individuals centres upon exploiting a mechanically stress resistant seed. Given this dietary consistency, it may be expected that the masticatory system will be integrated across development; such that individuals of all ages are capable of processing this food item. Using virtual modelling techniques, changes in bite pressure along the postcanine dental row were predicted at various stages of development. Results indicate that while maximum muscle force production increase...

In recent years finite element analysis (FEA) has emerged as a useful tool for the analysis of skeletal form-function relationships. While this approach has obvious appeal for the study of fossil specimens, such material is often fragmentary with disrupted internal architecture and can contain matrix that leads to errors in accurate segmentation. Here we examine the effects of varying the detail of segmentation and material properties of teeth on the performance of a finite element model of a Macaca fascicularis cranium within a comparative functional framework. Cranial deformations were compared using strain maps to assess differences in strain contours and Procrustes size and shape analyses, from geometric morphometrics, were employed to compare large scale deformations. We show that a macaque model subjected to biting can be made solid, and teeth altered in material properties, with minimal impact on large scale modes of deformation. The models clustered tightly by bite point rather than by modeling simplification approach, and fell out as being distinct from another species. However localized fluctuations in predicted strain magnitudes were recorded with different modeling approaches, particularly over the alveolar region. This study indicates that, while any model simplification should be undertaken with care and attention to its effects, future applications of FEA to fossils with unknown internal architecture may produce reliable results with regard to general modes of deformation, even when detail of internal bone architecture cannot be reliably modeled.

Multi-body dynamics is a powerful engineering tool which is becoming increasingly popular for the simulation and analysis of skull biomechanics. This paper presents the first application of multi-body dynamics to analyse the biomechanics of the rabbit skull. A model has been constructed through the combination of manual dissection and three-dimensional imaging techniques (magnetic resonance imaging and micro-computed tomography). Individual muscles are represented with multiple layers, thus more accurately modelling muscle fibres with complex lines of action. Model validity was sought through comparing experimentally measured maximum incisor bite forces with those predicted by the model. Simulations of molar biting highlighted the ability of the masticatory system to alter recruitment of two muscle groups, in order to generate shearing or crushing movements. Molar shearing is capable of processing a food bolus in all three orthogonal directions, whereas molar crushing and incisor biting are predominately directed vertically. Simulations also show that the masticatory system is adapted to process foods through several cycles with low muscle activations, presumably in order to prevent rapidly fatiguing fast fibres during repeated chewing cycles. Our study demonstrates the usefulness of a validated multi-body dynamics model for investigating feeding biomechanics in the rabbit, and shows the potential for complementing and eventually reducing in vivo experiments.

The development of virtual methods for anatomical reconstruction and functional simulation of skeletal structures offers great promise in evolutionary and ontogenetic investigations of form–function relationships. Key developments reviewed here include geometric morphometric methods for the analysis and visualization of variations in form (size and shape), finite element methods for the prediction of mechanical performance of skeletal structures under load and multibody dynamics methods for the simulation and prediction of musculoskeletal function. These techniques are all used in studies of form and function in biology, but only recently have they been combined in novel ways to facilitate biomechanical modelling that takes account of variations in form, can statistically compare performance, and relate performance to form and its covariates. Here we provide several examples that illustrate how these approaches can be combined and we highlight areas that require further investigation and development before we can claim a mature theory and toolkit for a statistical biomechanical framework that unites these methods.

Finite elements analysis (FEA) is now used routinely to interpret skeletal form in terms of function in both medical and biological applications. To produce accurate predictions from FEA models, it is essential that the loading due to muscle action is applied in a physiologically reasonable manner. However, it is common for muscle forces to be represented as simple force vectors applied at a few nodes on the model's surface. It is certainly rare for any wrapping of the muscles to be considered, and yet wrapping not only alters the directions of muscle forces but also applies an additional compressive load from the muscle belly directly to the underlying bone surface. This paper presents a method of applying muscle wrapping to high-resolution voxel-based finite element (FE) models. Such voxel-based models have a number of advantages over standard (geometry-based) FE models, but the increased resolution with which the load can be distributed over a model's surface is particularly advantageous, reflecting more closely how muscle fibre attachments are distributed. In this paper, the development, application and validation of a muscle wrapping method is illustrated using a simple cylinder. The algorithm: (1) calculates the shortest path over the surface of a bone given the points of origin and ultimate attachment of the muscle fibres; (2) fits a Non-Uniform Rational B-Spline (NURBS) curve from the shortest path and calculates its tangent, normal vectors and curvatures so that normal and tangential components of the muscle force can be calculated and applied along the fibre; and (3) automatically distributes the loads between adjacent fibres to cover the bone surface with a fully distributed muscle force, as is observed in vivo. Finally, we present a practical application of this approach to the wrapping of the temporalis muscle around the cranium of a macaque skull.