ORIGINAL RESEARCH article

Front. Ecol. Evol., 10 October 2022

Sec. Paleoecology

Volume 10 - 2022 | https://doi.org/10.3389/fevo.2022.957038

Reconstructing diets of hunted sika deer from Torihama Shell Midden site (ca. 6,000 years ago) by dental microwear texture analysis

  • 1. Shiga Prefectural Association for Cultural Heritage, Otsu, Shiga, Japan

  • 2. Department of Archaeology and Ethnology, Faculty of Letters, Keio University, Minato-ku, Tokyo, Japan

  • 3. Department of Natural Environmental Studies, Graduate School of Frontier Sciences, The University of Tokyo, Kashiwa, Chiba, Japan

Article metrics

View details

1

Citations

3,8k

Views

801

Downloads

Abstract

Sika deer (Cervus nippon) is the most abundant ruminant in the Japanese archipelago and has been the primary hunting target, including during the prehistoric ages. Abundant skeletal remains of sika deer have been excavated from archeological sites of the Jomon periods (ca. 15,000–2,400 BP). We reconstructed the feeding habits of sika deer from the Torihama Shell Midden site in Fukui Prefecture, western Japan. The Torihama site is one of the most well-preserved archeological sites of the Jomon period, and we investigated materials from the layers of the Early Jomon period (ca. 6,000 BP). In this study, we obtained three-dimensional tooth surface texture from the lower molars of excavated deer and measured microwear texture using international surface roughness parameters (ISO 25178-2) to infer their habitat use. Next, we estimated the percentage of grasses in the diet by using the reference dataset for extant sika deer with known diets. The results show that the Torihama deer overlapped with mixed feeding and grazing sika deer populations. Moreover, the proportion of grasses in the diet was estimated to be 50.7% on average but showed a wide range among the Torihama deer. This result implied that Torihama deer were mixed feeders of dicot leaves and grasses and had a flexible diet adapted to the vegetation of its habitat. Our results support the paleoenvironmental estimation that the Torihama Shell Midden site contained mixed vegetation of evergreen, deciduous, and coniferous trees around the Mikata Five Lakes in the Early Jomon period. These findings provide insights into the highly plastic diets of the extant sika deer in the Japanese archipelago.

Introduction

Sika deer (Cervus nippon), a medium-sized species of deer that inhabit a wide area of East Asia, from Amurland in the north to Vietnam in the south, are distributed in the Japanese archipelago from Hokkaido to the Kerama Islands, where there are five subspecies (Ohtaishi, 1986; Miura, 2008). Sika deer are thought to have migrated from the Asian continent to the Japanese archipelago during MIS12 (approximately 430,000 years ago) (Kawamura, 2009). Although the acidic soil in Japan is not conducive to the preservation organic remains, animal skeletal remains have been found, along with artifacts in prehistoric shell middens, caves, and low wetland sites. Skeletal remains of sika deer and wild boar (Sus scrofa) have often been excavated from those archeological sites, and these animals are thought to be among the major animals hunted during the Jomon period (ca. 15,000–2,400 BP) (Nishimoto, 1991).

From Torihama shell midden, abundant skeletal remains of sika deer have been excavated (Figure 1). This site was occupied by Paleolithic people from the Incipient to Early Jomon period. Because of the environmental condition of wetlands, this site has yielded various artifacts, for example, Jomon pottery and stone and bone tools; botanical artifacts, namely, fibers and wooden tools, represented by a well-preserved comb coated with red lacquer (Torihama Shell Midden Study Group, 1987); and many plant and animal remains. Therefore, Torihama Shell Midden has been one of the most intensively studied archeological sites in Japan. This research has been conducted by an interdisciplinary research team of archeologists and natural science researchers motivated by this slogan, “the reconstruction of the Jomon lifestyle focusing on their subsistence” (Kojima, 2015). Analyses of abundant sika deer and wild boar remains have been conducted based on various objectives and methodologies, for example, species identification and the quantitative analysis of each skeletal element (Morikawa, 1963; Nishida, 1979; Inami, 1983; Shigehara et al., 1991; Uchiyama, 2000; Anezaki et al., 2005; Sato, 2021a), the estimation of hunting seasons (Nishida, 1980; Ohtaishi, 1980; Uchiyama, 2005; Sato, 2021b), the examination of cut marks (Hongo, 1991), and the manufacturing of bone tools (Yamakawa, 1992a,b).

FIGURE 1

Studies have also been conducted to reconstruct the paleoenvironment and the use of plant resources by analyzing pollen (Yasuda, 1979), seed (Kasahara, 1983, 1984), phytolith (Sasaki, 1983), natural wood (Noshiro and Suzuki, 1990), and wood tools (Noshiro et al., 1996; Kudo et al., 2016a) in the sediments of the site. In addition, pollen analysis of lakes and lowland sediments has provided insights into long-term climatic changes around the site (Yasuda, 1982; Takahara and Takeoka, 1992; Nakagawa et al., 2002, 2005; Kitagawa et al., 2018).

The vegetation around the site in the Early Jomon period was estimated to be evergreen broadleaved forests, consisting of evergreen oak (Quercus subgen Cyclobalanopsis), chinquapin (Catanopisis sieboldii), and camellia (Camellia japonica Linn.) on the top of the hill; Japanese cedar (Cryptomeria japonica) on the hillsides; and deciduous broadleaved forests consisting of willow (Salix) and alder (Alnus sect. Gymnothyrsus) and Japanese ash (Fraxinus) in the waterside low areas (Noshiro and Suzuki, 1990). Subsequently, Yoshikawa et al. (2016) analyzed excavated pollen and seeds; their results suggest that the flora around the site was mixed vegetation: evergreen trees, broadleaved trees, and coniferous trees.

Therefore, clarifying the diets of sika deer in such a vegetation environment would help uncover the habitat use of the deer in the ecosystem and humans’ hunting activities. To estimate the diets of the excavated deer, we applied dental microwear texture analysis (DMTA) for the first time in the Jomon archeological sites. The DMTA measured the surface roughness of microscopic tooth enamel surfaces in three dimensions, obtained by optical profilometers such as confocal microscopes. The application of DMTA has been expanded to various types of vertebrates for dietary estimation (DeSantis, 2016) and has been most intensively applied to herbivores (e.g., Merceron et al., 2010; Berlioz et al., 2017), including extant, extinct, and excavated animals (Kubo and Fujita, 2021).

Because tooth surfaces retain dietary signals for a shorter period (days to weeks) than stable carbon and nitrogen isotopes of soft tissues and bones (Matsubayashi and Tayasu, 2019; Winkler et al., 2020), seasonal differences in a population can be detected, providing insights into hunting seasons and sites occupied by humans (Rivals et al., 2009, 2015). These results reveal human activity but are based on observations at low magnification; thus, they are not applicable to the method used in this study to quantitatively observe the three-dimensional (3D) shape of microwear. Kubo et al. (2017) and Kubo and Fujita (2021) have provided the reference dataset of extant sika deer sampled from various environments with quantitative dietary information; thus, the dental microwear texture (DMT) of the excavated deer can be compared with theirs.

In this study, we estimate the foraging behavior of the hunted sika deer from the Torihama Shell Midden site (hereafter “Torihama deer”) by using DMTA to infer the habitat use of the deer that lived during the Early Jomon period (ca. 6,000 BP).

Materials and methods

Materials

We used mandibles of the Torihama deer (n = 56) housed in the Wakasa History Museum and analyzed lower second molars (M2), which were not separated from the mandibles (Figure 2 and Supplementary Data 1). According to the correspondence between age in months and eruption of the molars (Niimi, 1997), the mandibles were 29 months old or older. Therefore, we observed the second molars, which had begun to wear down and had relatively few defects. They might have been derived from the same individual because the hemi-mandibles had been excavated in isolated conditions from the archeological site. However, in order to obtain the maximum number of samples, we included the left and right sides in the analysis. They were excavated during the ninth survey in 1984 (Torihama Shell Midden Study Group, 1985) and obtained from sediment layers yielding potteries of the Early Jomon period (Kitashirakawa Kaso II types). The layers, organic soils or shell deposits consisting primarily of freshwater mollusks, were dated to ca. 5,900–5,655 cal BP by radiocarbon dating (Kudo et al., 2016b). The deer remains showed evidence of human use, namely, scars from removal and dismantling and traces of extracting bone marrow and brain (Inami, 1983; Hongo, 1991; Yamakawa, 1992a; Uchiyama, 2000; Sato, 2021a). Additionally, many stone arrowheads, thought to have been used for hunting, were excavated (Torihama Shell Midden Study Group, 1987). These suggested that the bones probably originated from hunted individuals and were disposed of in the lowlands after being used as resources for food and tools. We expected those materials to include dead individuals from all seasons. The Torihama Shell Midden site yielded the evidence of hunting activity throughout the year, though it seemed to be intensively used from fall to winter based on the death-season of sika deer remains which were estimated from observations of antlers, the cementum layer, and mandibular tooth eruption (Nishida, 1980; Ohtaishi, 1980; Sato, 2021b).

FIGURE 2

We compared the data from the archeological remains with the published DMT data of extant sika deer. Among the 15 sika deer populations in Japan that Kubo and Fujita (2021) analyzed, we selected three populations (Figure 1 and Table 1): typical grazing (Kinkazan Island), browsing (Tsushima Island), and mixed-feeding (Shimane) populations. The Shimane population was selected because of its habitat’s geological proximity and environmental similarity to the Torihama Shell Midden location. We selected the Kinkazan and Tsushima populations because of their contrasting north–south diets in the Japanese archipelago, and many samples have been analyzed. Sika deer on Kinkazan Island inhabit a small island near the Oshika Peninsula in Miyagi Prefecture, facing the Pacific Ocean. The main vegetation is deciduous broadleaved forest with some open grasslands. Kinkazan deer feed mainly on grasses such as silver grass (Miscanthus sinensis) and lawn grass (Zoysia japonica) (Takatsuki, 1980). The Tsushima deer inhabited an evergreen broadleaved forest. Stomach content analysis revealed that they fed mainly on woody plants, followed by herbaceous plants, seeds, and fruits; graminoids were a minor component of their diet (Suda, 1997). Deer from the Shimane population inhabited the Misen Mountains of Shimane Prefecture, in the San-in region, at the same latitude as the Torihama site. Evergreen broadleaved forests are the primary vegetation in the area. Stomach content analysis of Shimane deer showed that they consumed both graminoids and leaves of broadleaved trees and forbs year-round, although there was seasonal variation (Shimane Prefectural Government, 2002). The dietary differences among the three populations were represented in the percentage of grasses in the diets: 65.8% in Kinkazan (grazers), 38.1% in Shimane (mixed feeders), and 3.4% in Tsushima (browsers) (Kubo and Yamada, 2014; Kubo and Fujita, 2021).

TABLE 1

Region of the materialsSpecimen no.Scan no.Annual mean temperature (°C)Annual precipitation (mm)Major vegetation% Grass
in diet*1
Kinkazan Island3212311.91,150Deciduous broad-leaves forest with some open grasslands65.8
Shimane2710014.61,685Evergreen broad-leaves forest38.1
Tsushima Island2812915.82,235Evergreen broad-leaves forest3.4
Torihama shell middenM2:56 (L:31, R:25)20515*22,000*2Evergreen broad-leaves forest with Japanease cedar and Deciduous broad-leaves

Habitat environment and food habits related to the materials in this study.

*1We referenced Kubo and Yamada (2014) and Kubo and Fujita (2021) which contained the quantitative dietary data showing the consumption rate of graminoids obtained from stomach content analysis and fecal analysis. *2Temperature and precipitation around the site are estimated from 5,900 to 5,600 cal.bp based on the results of Kitagawa et al. (2018).

Methods: Dental microwear texture analysis

Molding of tooth enamel surfaces of Torihama deer

We used the methodology of Kubo et al. (2017) and made molds of the M2. To remove dirt adhered during the sedimentation, we cleaned the molar occlusal surfaces with cotton swabs soaked with acetone. Next, the surface was molded using high-resolution A-silicone dental impression material (Dr. Silicon regular type, BSA Sakurai, Japan). We selected the occlusal enamel band as our target for scanning (Figure 2) because we expected it to be less affected than other parts by factors other than diet (Schulz et al., 2010).

Obtaining dental microwear texture (DMT) data and calculating DMT parameters

We scanned the dental impressions by using a confocal laser microscope (VK-9700, Keyence, Osaka, Japan) equipped with a 100× long-distance lens (N.A. = 0.95). Three-dimensional data was obtained with a lateral (x, y) sampling interval of 0.138 μm and a vertical resolution (z) of 0.001 μm. Occlusal enamel facets were scanned on areas of 140 × 105 μm, and grayscale images of the tooth surfaces were acquired for graphical comparison. The 3D surface data were processed by applying the methods in Yamada et al. (2018), Aiba et al. (2019), and Kubo and Fujita (2021), using MountainsMap Imaging Topography (7.4.8872, Digital Surf, Besançon, France). The methods provide an objective evaluation of quantitative data and scratch depth by measuring microwear three-dimensionally. First, we scanned and trimmed the 135 × 100 μm area of the microwear impressions, avoiding dust and large scratches. Because the scanned images were mirror images of real tooth surfaces, the coordinates were mirrored in the x- and z-axes. Second, through the least square plane by subtraction, the surfaces were leveled to remove the inclination of the molds. Next, we applied a robust Gaussian filter with a cutoff value of 0.8 μm to remove measurement noise (S-filter as defined in ISO25178) and the form removal function (polynomial of increasing power = 2) to remove the large-scale curvature of the enamel bands (F-operation in ISO25178). We used the automated outlier removal function of MountainsMap, which removes features with a slope > 80°, and a threshold to remove the upper and lower 0.1% of the data. Scans with a non-measured point above 0.5% of total data points were discarded from the following analyses. The non-measured points were filled by using the smoothing function of MountainsMap. Finally, the values of 31 parameters of ISO 25178-2 were calculated, of which the names and definitions are shown in Table 2. These 31 ISO parameters were used in Kubo and Fujita (2021) to investigate the relationship between DMTA parameters and the extant sika deer diet. Because we used the same analytical template of the same software (MountainsMap) that Kubo and Fujita (2021) used, their published parameters are directly comparable to ours.

TABLE 2

CategoryParameterDescriptionUnit
Height parametersSqStandard deviation of the height distributionμm
SskSkewness of the height distributionno unit
SkuKurtosis of the height distributionno unit
SpMaximum peak height, height between the highest peak and the mean planeμm
SvMaximum pit height, depth between the mean plane and the deepest pitμm
SzMaximum height, sum of the maximum peak height and the maximum pit height (Sp + Sv)μm
SaArithmetic mean heightμm

Spatial parametersSalAutocorrelation length (s = 0.2)μm
StrTexture aspect ratio (s = 0.2)no unit

Hybrid parametersSdqRoot mean square gradientno unit
SdrDeveloped interfacial area ratio%

Functional parametersSmrAreal material ratio, ratio of the area of the material at a specified height c (c = 1 μm under the highest peak)%
SmcInverse areal material ratio, height at which a given areal material ratio (p = 10%)μm
SxpPeak extreme height, difference in height between the p and q material ratio (p = 50%, q = 97.5%)μm
SkDistance between the highest and lowest level of the core surfaceμm
SpkAverage height of the protruding peaks above the core surfaceμm
SvkAverage height of the protruding dales below the core surfaceμm

Functional volume parametersVmMaterial volume at a given material ratio (p = 10%)μm3/μm2
VvVoid volume at a given material ratio (p = 10%)μm3/μm2
VmcMaterial volume of the core at a given material ratio (p = 10%, q = 80%)μm3/μm2
VvcVoid volume of the core (p = 10%, q = 80%)μm3/μm2
VvvVoid volume of the dale at a given material ratio (q = 80%)μm3/μm2

Feature parametersSpdDensity of peaks1 /μm2
SpcArithmetic mean peak curvature1 /μm
S10zTen-point heightμm
S5pFive-point peak heightμm
S5vFive-point pit heightμm
SdaClosed dale areaμm2
ShaClosed hill areaμm2
SdvClosed dale volumeμm3
ShvClosed hill volumeμm3

Names and definitions of 31 DMT parameters.

Statistical analyses

To clarify the feeding habits of the Torihama deer, we compared their DMT parameters with those of the three populations of sika deer (Kinkazan Island, Shimane, and Tsushima Island). We conducted principal component analysis (PCA) followed by a varimax rotation to interpret the 31 parameters and summarize them into a few principal components (PCs). We also conducted the Mann–Whitney U-test to determine differences in the parameters and the PC scores between Torihama and the three populations. The statistical analyses were conducted using SPSS statistics ver.28 (IBM, Chicago, USA).

Estimation of the percentage of grasses (% grass) in the Torihama deer diet

Kubo and Fujita (2021) presented an equation, as follows, to estimate the % grass in diets from an ISO 25178-2 parameter, Sk, which showed the highest correlation with the % grass in the diet in their analysis:

where logit (% grass) is transformed into the % grass by the equation

We used those two equations to estimate the % grass in the diet of 56 Torihama deer.

Results

Assessment of tooth microwear of Torihama deer by two-dimensional images

Representative microwear images and 3D models of the tooth surfaces are shown in Figure 3. Of the two-dimensional (2D) images, we excluded those in which the tooth surface was clearly damaged or the non-measured point was above 0.5% of the total data points. As a result, of the 224 scans (N = 224), we used data from 205 scans (Supplementary Data 1). Although some materials showed heavily worn surfaces, microwear features were well-preserved in the Torihama deer specimens, and the diagenetic alteration was considered negligible. Small pits were commonly observed among the Torihama deer samples. Some samples had flat surfaces with a visible enamel prism structure (Figure 3A: FTS84-086), and others showed abundant scratches (Figure 3B: FTS084-096 and Figure 3C: FTS84-054). Coarse scratches were observed on the surface of some individuals (Figure 3C).

FIGURE 3

FIGURE 4

Comparison between Torihama deer and extant deer from representative populations

The PCA of the 31 DMT parameters showed that the first, second, and third components explained approximately 52.6, 15.9, and 12.8% of the total variance, respectively. Table 3 indicates the factor loadings of the 31 parameters. The first component (PC1) can be interpreted as overall surface roughness because 22 of the 31 parameters had factor loadings higher than the critical threshold value of 0.4, and 20 parameters were larger than 0.8. The height (e.g., Sq, Sz, Sa) and volume (e.g., Vm, Vv, Vmc, Vvc, Vvv) parameters had positive loads; thus, the larger the PC1 value, the larger the surface relief. On the other hand, parameters related to surface segmentation (Sda, Sha, Sdv, Shv, and Spd) contributed significantly to PC2. The factor loadings of Sda, Sha, Sdv, and Shv, which indicate the areas and volumes of segmented hills and dales on the surface, were positive values, and that of Spd, which indicates the density of the peaks, was negative. Therefore, PC2 can be interpreted as the fineness of the microwear features, with the smaller values indicating a surface characterized by abundant fine features.

TABLE 3

PC1PC2PC3
Variance explained (%)52.65915.93312.832
Sq0.9030.2300.329
Ssk0.0800.0560.940
Sku–0.037–0.081−−0.894
Sp0.8870.1090.348
Sv0.9100.204–0.191
Sz0.9600.1700.068
Sa0.8750.2280.392
Smr−−0.822–0.107–0.321
Smc0.8560.2130.439
Sxp0.9190.2450.131
Sal0.2240.2890.593
Str–0.006–0.053–0.184
Sdq0.895–0.393–0.101
Sdr0.866−−0.405–0.076
Vm0.8170.1790.384
Vv0.8580.2130.439
Vmc0.8410.2260.443
Vvc0.8370.2070.476
Vvv0.912–0.231–0.012
Spd0.026−−0.752–0.302
Spc0.5060.095–0.007
S10z0.8810.081–0.087
S5p0.756–0.1090.068
S5v0.8240.231–0.206
Sda0.0870.9310.123
Sha–0.0790.9460.030
Sdv0.3380.8330.021
Shv0.1840.878–0.028
Sk0.892–0.1490.224
Spk0.898–0.303–0.092
Svk0.812–0.093−−0.413

Factor loadings obtained from PCA conducted using 31DMT parameters.

Parameters with factor loadings higher than the critical threshold of 0.4 are shown in bold.

We found significant differences between the Torihama deer and the three populations of extant deer (P < 0.05; Table 4). The comparison with the Kinkazan deer showed statistically significant differences in 23 parameters, except for some hybrid, functional, and feature parameters. Statistically significant differences between the Torihama deer and the Shimane deer were found in 21 parameters, excluding the spatial, functional volume, and feature parameters. The largest number of statistical differences (28 parameters) was found in the comparison between the Torihama deer and the Tsushima deer, indicating that the difference between the Torihama deer and the Tsushima deer is larger than that between the deer in Kinkazan and Shimane.

TABLE 4

Torihama vs. Kinkazan
Torihama vs. Shimane
Torihama vs. Tsushima
ParameterzPzPzP
Sq−−4.910< 0.0012.1000.0362.6470.008
Ssk−−4.944< 0.001−−3.723< 0.001−−2.9600.003
Sku5.465< 0.0013.703< 0.0013.1120.002
Sp−−3.817< 0.0012.7510.0063.596< 0.001
Sv–0.9980.3194.782< 0.0014.735< 0.001
Sz−−2.5590.0114.219< 0.0014.678< 0.001
Sa−−5.152< 0.0011.5940.1112.2110.027
Smr3.1310.002−−2.6630.008−−3.359< 0.001
Smc−−5.152< 0.0011.3030.1931.9930.046
Sxp−−3.869< 0.0012.8090.0053.2550.001
Sal0.2950.768–0.7780.4370.2470.805
Str4.779< 0.0011.5260.1272.7520.006
Sdq−−2.3770.0176.250< 0.0016.841< 0.001
Sdr–1.6310.1036.610< 0.0017.003< 0.001
Vm−−4.294< 0.0011.3800.1682.0590.039
Vv−−5.126< 0.0011.3610.1742.0210.043
Vmc−−5.447< 0.0011.0980.2721.8410.066
Vvc−−5.248< 0.0010.9530.3411.7840.074
Vvv−−3.340< 0.0013.0420.0023.615< 0.001
Spd3.0880.0023.791< 0.0015.845< 0.001
Spc−−2.8970.0042.9260.0032.7520.006
S10z–0.2520.8014.481< 0.0014.583< 0.001
S5p0.7890.4303.528< 0.0014.346< 0.001
S5v–1.1020.2714.219< 0.0013.568< 0.001
Sda−−3.652< 0.001−−3.392< 0.001−−5.238< 0.001
Sha−−2.9150.004−−3.626< 0.001−−5.987< 0.001
Sdv−−2.6280.0090.2330.816−−2.6850.007
Shv−−3.478< 0.001–0.3010.763−−3.909< 0.001
Sk−−5.456< 0.0014.335< 0.0014.479< 0.001
Spk–1.9340.0535.434< 0.0016.433< 0.001
Svk–1.0150.3105.550< 0.0014.706< 0.001
PC1−−3.2180.0014.831< 0.001−−5.371< 0.001
PC2–1.7170.086−−2.6540.0085.105< 0.001
PC3−−4.589< 0.001−−4.121< 0.0012.3340.020

Comparison of DMT parameters and principal component scores between Torihama deer and extant deer.

P-values of Mann–Whitney U-tests are shown. Significant (P < 0.05) differences are shown in bold.

The results of the statistical comparison of the PC scores between the Torihama and extant deer populations are also shown in Table 4. A scatter plot of PC1 and PC2 scores are shown in Figure 4 for the populations of Torihama, Shimane, Kinkazan, and Tsushima. Deer on Kinkazan Island, typical grazers, were characterized by higher PC1 scores (i.e., large surface roughness), whereas those deer on Tsushima Island, typical browsers, were located at the lower bottom. The Shimane population, mixed feeders, was located between them. Thus, we found a separation between the grazing, mixed feeding, and browsing populations of the extant deer, although the Shimane and Tsushima populations overlapped. The Torihama deer overlapped with the distribution of the grazing and mixed-feeding sika deer but not with the browsing population.

Estimation of % grass in the diet of Torihama deer and its variation within the assemblage

We applied the estimation equation of Kubo and Fujita (2021) to estimate the % grass in the diet of 56 Torihama deer (Figure 5). The Torihama deer were estimated to have consumed a diet of, on average, 50.7% grass, which supported a general mixed-feeding habit of the Torihama deer. Browsing and grazing individuals occurred with a similar frequency, and there was a wide range in % grass within the Torihama deer population.

FIGURE 5

Discussion

In this study, we used DMTA to estimate the diets of Torihama deer inhabited approximately 6,000 years ago. Deer mandibles excavated from the site were in satisfactory preservation conditions; thus, observing the antemortem microwear was possible. The results of the comparison with the extant sika deer showed that the Torihama deer overlapped with the populations that had mixed-feeding and grazing-feeding habits and showed a wide range in % grass in the diet.

Reconstruction of the diet of the Torihama deer based on extant sika deer references

The results of PCA on 31 parameters of surface roughness showed that the Torihama deer overlapped with the grazing and mixed-feeding deer. However, they did not overlap with the browsing Tsushima deer. The Tsushima deer did not have well-developed microwear. A reason for this observation is that their diets were less abrasive than those of the Torihama deer. For example, the former ate dicot leaves, which polished tooth surfaces rather than abrading them. The absence of overlap between the Tsushima and the Torihama deer indicates that the diet of the Torihama deer was not primarily dicots. The Torihama deer overlapped with the Kinkazan deer, a typical grazer type that developed microwear by feeding on grass, and with the Shimane deer, a mixed feeder type, which demonstrates the possibility that the Torihama deer were feeding on both grass plants and the leaves and fruits of dicotyledonous plants.

Furthermore, some Torihama deer developed abundant scratches on their tooth surfaces, which were finer than those of grazing deer on Kinkazan Island (Figure 3B). These individuals did not overlap with the Kinkazan deer in the scatter plot of PC1 and PC2 (Figure 4). These results may indicate that some Torihama deer foraged on grass species, which differ from those foraged by the Kinkazan deer. The Kinkazan deer feed primarily on short-statured lawn grass (Z. japonica) (Takatsuki, 1980), which may force them to ingest the soil secondarily when feeding on lawn grass at the ground level (Rivals et al., 2014; Rivals and Takatsuki, 2015). Microwear can also be affected by sand and dust attached to foods (Maas, 1994; Gügel et al., 2001; Lucas et al., 2013). However, experimental feeding has clarified that the effects of external abrasives are mediated by the rumination process to some extent (Ackermans et al., 2020; Schulz et al., 2020). The comparison of 15 populations of sika deer demonstrated that the Kinkazan deer had rougher surfaces than other grazing deer (Kubo and Fujita, 2021). Therefore, the Torihama deer might not have fed on short grasses contaminated by soil. Around the Torihama Shell Midden site grew graminoid (Poaceae and Cyperaceae) plants, indicated by the analyses of pollen and macrobotanical remains (Kasahara, 1983, 1984; Kitagawa et al., 2018). The analysis of phytoliths has also identified grasses such as the reed Phragmites and Bambusoideae species (Sasaki, 1983). Therefore, the Torihama deer could also have fed on these grasses when they were available.

The flexible diet of Torihama deer

The proportion of grass in the diets of the Torihama deer was estimated to be 50.7% on average, with a wide range of 0 to 100% among these individuals (Figure 5). This estimate indicates that Torihama deer are mixed feeders and may reflect the dietary flexibility of sika deer during the Jomon period. It is also assumed that the Torihama deer would change diets according to seasonal availability. The Torihama deer in this study were over 2 years old, representing the diets of weaned individuals. Furthermore, the average age in the Torihama population was estimated to be 6.96 years (Koike and Ohtaishi, 1985); thus, the tendencies of sexually mature individuals can be observed. We could not estimate the mortality seasons for each material used in this study. However, the literature has suggested that the people who lived at the Torihama Shell Midden site may have hunted sika deer year-round, mainly during fall and winter (Nishida, 1980; Ohtaishi, 1980; Sato, 2021b).

The result of stomach content analysis of the Shimane populations, which we used as the reference in this study, showed a seasonal change in diets of deer with a high consumption rate of grasses in summer, seeds and nuts in fall, and tree leaves in winter (Shimane Prefectural Government, 2002). Such seasonal change in diets was also identified for sika deer populations in Hyogo, the closest area to the archeological site; the consumption of seeds and nuts was high in fall; and that of evergreen leaves was high in winter (Yokoyama, 2009). Therefore, although we did not estimate the mortality seasons, the individuals that died in each season were analyzed by the DMTA. Based on the seasonal changes in the diets of the extant sika deer with geographical proximity and the background information that the Torihama deer may have been hunted year-round, the DMT of the Torihama deer may reflect their seasonal dietary changes. As mentioned in the introduction, the flora around the site was mixed vegetation of evergreen trees, broadleaved trees, and coniferous trees, which allowed the Torihama deer to change their forage plants seasonally.

Additionally, the extant sika deer have dietary flexibility and can inhabit a wide range of habitat environments in Japan (Takatsuki, 2009). Takatsuki (2006) suggested that the ecological plasticity of the sika deer is associated with the extinction of large herbivores in the Japanese archipelago around the Pleistocene–Holocene boundary, through the process that vacant ecological niches after those extinctions were occupied by a survived species, sika deer. The argument of Takatsuki (2006) is vague regarding the cause–effect relationship, namely, “sika deer could survive the extinction because they had the ability to adapt to changing environments” or “since there were vacant ecological niches after the extinction, the sika deer became ecologically diverse according to the respective habitats.” To clarify this issue, information on the ecological variability of Pleistocene and Holocene sika deer, as well as other extinct ruminants, is necessary. This study is the first step toward this clarification. Our study indicates that sika deer in the Early Jomon period (approximately 6,000 years ago) had a flexible feeding habit of adapting to the vegetation of their habitats. Further investigation of excavated sika deer from older archeological sites than the site we examined would provide insights into when and how they attained high ecological plasticity.

Conclusion

We investigated the feeding habits of the sika deer in the Early Jomon period by analyzing the DMT of abundant skeletal remains excavated from an archeological site. The Torihama deer were mixed feeders, with a wide range in the estimated percentage of foraged grasses. These results are consistent with the estimate of paleo-vegetation around the Torihama Shell Midden site. Similar to the extant deer populations inhabiting the evergreen broadleaved forest in Honshu, the Torihama deer might have changed their foraging plants seasonally, depending on the growth of plants around their habitat. Their flexible diets are important for understanding how they survived the extinction of large mammals at the end of the Pleistocene and are widely present around the Japanese archipelago today. Most of the Japanese extant samples were collected during the opening of hunting seasons, and in the unhunted population of Kinkazan deer, they were individuals that died of natural causes. Therefore, because the comparative data are mainly from autumn to early spring, further research on seasonal variation in microwear should be conducted by obtaining samples that can be compared among seasons.

For further research, DMTA results will be able to further elucidate the paleoecology of deer by including studies of the mortality season and range of behavior in individuals.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material. The original scan data of all individuals in common 3D surface format (“.sur” format) and 3D surface analytical files including raw scans and procedures of analyses in DMTA standard format (“.mnt”) were uploaded to an online repository (Zenodo): https://doi.org/10.5281/zenodo.7058460.

Ethics statement

Ethical review and approval was not required for the animal study because it is focused on archeological materials.

Author contributions

KS: conceptualization, investigation, formal analysis, data curation, visualization, and writing—original draft preparation. TS: conceptualization, resources, supervision, and writing—review and editing. MOK: resources, supervision, and writing—review and editing. All authors contributed to the article and approved the submitted version.

Funding

This study was partially funded by JSPS KAKENHI (Grant No. 16H03106) to TS and the Keio University Doctorate Student Grant-in-Aid Program from Ushioda Memorial Fund to KS.

Acknowledgments

We thank the following institutions and people for their support. The Wakasa History Museum permitted us to use the collections. Mayumi Ajimoto (Wakasa History Museum) and Hideaki Kojima (Wakasa Town Section of History and Culture) provided us with suggestions and support. In addition, we received valuable comments from professors of the Department of Archaeology and Ethnology, Keio University, and the Mammal Society of Japan. We would like to thank Enago (www.enago.jp) for the English language review. Finally, we appreciate the two reviewers for their constructive comments on the earlier version of this manuscript.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fevo.2022.957038/full#supplementary-material

References

  • 1

    AckermansN. L.WinklerD. E.MartinL. F.KaiserT. M.ClaussM.HattJ.-M. (2020). Dust and grit matter: Abrasives of different size lead to opposing dental microwear textures in experimentally fed sheep (Ovis aries).J. Exp. Biol.223:jeb220442. 10.1242/jeb.220442

  • 2

    AibaK.MiuraS.KuboM. O. (2019). Dental microwear texture analysis in two ruminants, Japanese serow (Capricornis crispus) and sika deer (Cervus nippon), from central Japan.Mamm. Study44183192. 10.3106/ms2018-0081

  • 3

    AnezakiT.NishimotoT.NiimiM. (2005). Faunal analysis of the Torihama shell midden, materials from the area 1 of the 1985 excavation (In Japanese).Torihama Shell Midden Pap.4-51930.

  • 4

    BerliozE.AzoritC.BlondelC.RuizM. S. T.MerceronG. (2017). Deer in an arid habitat: Dental microwear textures track feeding adaptability.Hystrix28222230. 10.4404/hystrix-28.2-12048

  • 5

    DeSantisL. (2016). Dental microwear textures: Reconstructing diets of fossil mammals.Surf. Topogr.4:e023002. 10.1088/2051-672X/4/2/023002

  • 6

    GügelI. L.GrupeG.KunzelmannK. H. (2001). Simulation of dental microwear: Characteristic traces by opal phytoliths give clues to ancient human dietary behavior.Am. J. Phys. Anthropol.114124138. 10.1002/1096-8644(200102)114:2<124::AID-AJPA1012<3.0.CO;2-S

  • 7

    HongoH. (1991). Analysis of cut marks on mammal remains from section 85 of the Torihama shell midden (In Japanese).Bull. Natl. Mus. Jpn. Hist.29149195.

  • 8

    InamiM. (1983). “Remains of sika deer and wild boar from the Torihama shell midden,” in The summary report of excavation in 1981 and 1982 and research results of Torihama shell midden – research 3 on archaeological site located in lowland area from the early Jomon period (In Japanese), ed.Torihama Shell Midden Study Group (Wakasa Town: Board of Fukui Prefecture and Wakasa Museum of History and Folklore), 6581.

  • 9

    KasaharaY. (1983). “The detection and identification of plant seeds from Torihama shell midden (the 6th excavation),” in The summary report of excavation in 1981 and 1982 and research results of Torihama shell midden – research 3 on archaeological site located in lowland area from the early Jomon period (In Japanese), ed.Torihama Shell Midden Study Group (Wakasa Town: Board of Fukui Prefecture and Wakasa Museum of History and Folklore), 4764.

  • 10

    KasaharaY. (1984). “The detection and identification of plant seeds from the Torihama shell midden (7th excavation),” in The summary report of excavation in 1983 and research results of Torihama shell midden – research 4 on archaeological site located in lowland area from the early Jomon period (In Japanese), ed.Torihama Shell Midden Study Group (Wakasa Town: Board of Fukui Prefecture and Wakasa Museum of History and Folklore), 4990.

  • 11

    KawamuraY. (2009). “Fossil record of sika deer in Japan,” in Sika deer, edsMcCulloughD. R.TakatsukiS.KajiK. (Tokyo: Springer), 1125. 10.1007/978-4-431-09429-6_2

  • 12

    KitagawaJ.KojimaH.YoshidaT.YasudaY. (2018). Adaptations of the early Jomon people in their settlement relocation to climate change around Lake Mikata, Central Japan.Archaeol. Res. Asia166677. 10.1016/j.ara.2018.03.002

  • 13

    KoikeH.OhtaishiN. (1985). Prehistoric hunting pressure estimated by the age composition of excavated sika deer (Cervus nippon) using the annual layer of tooth cement.Archaeol. Sci.12443456. 10.1016/0305-4403(85)90004-4

  • 14

    KojimaH. (2015). The consideration of the excavation and study history of Torihama shell midden and its contribution to the zooarchaeological study (In Japanese with English abstract).Zooarchaeology321124.

  • 15

    KuboM. O.FujitaM. (2021). Diets of Pleistocene insular dwarf deer revealed by dental microwear texture analysis.Paleogeogr. Palaeoclimatol. Palaeoecol.562:e110098. 10.1016/j.palaeo.2020.110098

  • 16

    KuboM. O.YamadaE. (2014). The inter-relationship between dietary and environmental properties and tooth wear: Comparisons of mesowear, molar wear rate, and hypsodonty index of extant sika deer populations.PLoS One9:e90745. 10.1371/journal.pone.0090745

  • 17

    KuboM. O.YamadaE.KuboT.KohnoN. (2017). Dental microwear texture analysis of extant sika deer with considerations on inter-microscope variability and surface preparation protocols.Biosurf. Biotribol.3155165. 10.1016/j.bsbt.2017.11.006

  • 18

    KudoY.SuzukiM.NoshiroS.AjimotoM.AmitaniK. (2016a). Radiocarbon dating of Castanea crenata wood from the incipient and initial Jomon periods excavated from the Torihama shell midden site, Fukui Prefecture, Japan (In Japanese with English abstract).Jpn. Assoc. Hist. Bot.245968. 10.34596/hisbot.24.2_59

  • 19

    KudoY.AmitaniK.YoshikawaJ.SasakiY.AjimotoM.NoshiroS. (2016b). Radiocarbon dating of plant macrofossils excavated from the Torihama shell midden site, Fukui Prefecture, Japan—a chronological reconsideration of sediment layers and types of pottery from the incipient to early Jomon periods (In Japanese with English abstract).Jpn. Assoc. Hist. Bot.244358. 10.34596/hisbot.24.2_43

  • 20

    LucasP. W.OmarR.Al-FadhalahK.AlmusallamA. S.HenryA. G.MichaelS.et al (2013). Mechanisms and causes of wear in tooth enamel: Implications for hominin diets.J. R. Soc. Interface10:e20120923. 10.1098/rsif.2012.0923

  • 21

    MaasM. C. (1994). A scanning electron microscopic study of in vitro abrasion of mammalian tooth enamel under compressive loads.Arch. Oral Biol.39111. 10.1016/0003-9969(94)90028-0

  • 22

    MatsubayashiJ.TayasuI. (2019). Collagen turnover and isotopic records in cortical bone.J. Archaeol. Sci.1063744. 10.1016/j.jas.2019.03.010

  • 23

    MerceronG.EscarguelG.AngibaultJ.-M.Verheyden-TixierH. (2010). Can dental microwear textures record inter-individual dietary variations?PLoS One5:e9542. 10.1371/journal.pone.0009542

  • 24

    MiuraS. (2008). “Sika deer,” in A guide to the mammals of Japan, ed.Japan Wildlife Research Center (Tokyo: Tokai university press), 110111.

  • 25

    MorikawaM. (1963). Some problems concerning the Torihama shell midden site in Fukui Prefecture (In Japanese).Mater. Cult.11934.

  • 26

    MorikawaM. (2002). The Torihama shell midden- time capsule of the Jomon people.Tokyo: Miraisya.

  • 27

    NakagawaT.KitagawaH.YasudaY.TarasovP. E.GotandaK.SawaiY. (2005). Pollen/event stratigraphy of the varved sediment of Lake Suigetsu, Central Japan from 15,701 to 10,217 SG vyr BP (Suigetsu varve years before present): Description, interpretation, and correlation with other regions.Quat. Sci. Rev.2416911701. 10.1016/j.quascirev.2004.06.022

  • 28

    NakagawaT.TarasovP. E.NishidaK.GotandaK.YasudaY. (2002). Quantitative pollen-based climate reconstruction in Central Japan: Application to surface and Late Quaternary spectra.Quat. Sci. Rev.2120992113. 10.1016/S0277-3791(02)00014-8

  • 29

    NiimiM. (1997). The preliminary report of the determination of deer hunting seasons by the teeth eruption and wearing (In Japanese with English abstract).Zooarchaeology92132.

  • 30

    NishidaM. (1979). “Animal remains,” in Torihama shell midden – research 1 on archaeological site located in lowland area from the early Jomon period (In Japanese), ed.Torihama Shell Midden Study Group (Wakasa Town: Board of Fukui Prefecture), 164166.

  • 31

    NishidaM. (1980). Food resources and subsistence activities in the Jomon period on the organic remains of the Torihama shell midden (In Japanese).Anthropol. Q.11341.

  • 32

    NishimotoT. (1991). Hunted deer and wild boar in the Jomon period (In Japanese).Kodai91114132.

  • 33

    NoshiroS.SuzukiM. (1990). A palaeobotanical study of natural woods from the Torihama shell midden and reconstruction of the forest vegetation (In Japanese with English abstract).Bull. Jpn. Sea Res. Inst.2263152.

  • 34

    NoshiroS.SuzukiM.AmitaniK. (1996). Species of wooden products from the Torihama shell midden (In Japanese).Torihama Shell Midden Pap.12379.

  • 35

    OhtaishiN. (1980). Estimation of sex, age, and season of death using mandibles of Cervus nippon excavated from an archaeological site (In Japanese).Archaeol. Nat. Sci.135174.

  • 36

    OhtaishiN. (1986). Preliminary memorandum of classification, distribution and geographic variation on sika deer (In Japanese).Honyurui Kagaku531317.

  • 37

    RivalsF.SchulzE.KaiserT. M. (2009). A new application of dental wear analyses: Estimation of duration of hominid occupations in archaeological localities.J. Hum. Evol.56329339. 10.1016/j.jhevol.2008.11.005

  • 38

    RivalsF.TakatsukiS.AlbertR. M.MaciàL. (2014). Bamboo feeding and tooth wear of three sika deer (Cervus nippon) populations from northern Japan.J. Mammal.9510431053. 10.1644/14-MAMM-A-097

  • 39

    RivalsF.TakatsukiS. (2015). Within-island local variations in tooth wear of sika deer (Cervus nippon centralis) in northern Japan.Mammal. Biol.80333339. 10.1016/j.mambio.2015.02.001

  • 40

    RivalsF.PrignanoL.SemprebonG. M.LozanoS. (2015). A tool for determining duration of mortality events in archaeological assemblanges using extant ungulate microwear.Sci. Rep.5:e17330. 10.1038/srep17330

  • 41

    SasakiS. (1983). “Plant opal analysis of soil from Torihama shell midden,” in The summary report of excavation in 1981 and 1982 and research results of Torihama shell midden – research 3 on archaeological site located in lowland area from the early Jomon period (In Japanese), ed.Torihama Shell Midden Study Group (Wakasa Town: Board of Fukui Prefecture and Wakasa Museum of History and Folklore), 3741.

  • 42

    SatoK. (2021a). The resource utilization of sika deer from Torihama shell midden (In Japanese).Annu. Bull. Shiga Prefect. Assoc. Cult. Herit.34112.

  • 43

    SatoK. (2021b). Death-season of sika deer remains from the Torihama shell midden: Nondestructive analyses on mandibular molariform tooth development (In Japanese with English abstract).Zooarchaeology381122.

  • 44

    SchulzE.CalandraI.KaiserT. M. (2010). Applying tribology to teeth of hoofed mammals.Scanning32162182. 10.1002/sca.20181

  • 45

    SchulzE.WinklerD. E.ClaussM.CarlssonJ.AckermansN. L.MartinL. F.et al (2020). Everything matters: Molar microwear texture in goats (Capra aegagrus hircus) fed diets of different abrasiveness.Palaeogeogr. Palaeoclimatol. Palaeoecol.552:109783. 10.1016/j.palaeo.2020.109783

  • 46

    ShigeharaN.HongoH.AmitaniK. (1991). Mammalian remains from the Torihama shell midden excavation in 1985 (In Japanese).Bull. Natl. Mus. Jpn. Hist.29329341.

  • 47

    Shimane Prefectural Government (2002). A survey on the sika deer (Cervus nippon) in Misen mountainous region on the Shimane Peninsula, Japan (In Japanese).Shimane: Shimane Prefectural Government.

  • 48

    SudaK. (1997). Rumen contents and food selectivity of sika deer (Cervus nippon) on Tsushima Islands (In Japanese with English abstract).Wildl. Conserv. Jpn.2125134.

  • 49

    TakaharaH.TakeokaM. (1992). Vegetation history since the last glacial period in the Mikata lowland, the sea of Japan area, western Japan.Ecol. Res.7371386. 10.1007/BF02347104

  • 50

    TakatsukiS. (1980). Food habits of sika deer on Kinkazan Island.Sci. Rep. Tohoku Univ. Ser. IV Biol.38731.

  • 51

    TakatsukiS. (2006). Ecological history of sika deer.Tokyo: University of Tokyo Press.

  • 52

    TakatsukiS. (2009). “Geographical variations in food habits of sika deer: The Northern Grazer vs. the Southern Browser,” in Sika deer, edsMcCulloughD. R.TakatsukiS.KajiK. (Tokyo: Springer), 231237. 10.1007/978-4-431-09429-6_17

  • 53

    Torihama Shell Midden Study Group (1985). Summary report of excavation in 1984 and research results of Torihama shell midden – research 5 on archaeological site located in lowland area from the early Jomon period (In Japanese).Wakasa Town: Board of Fukui Prefecture and Wakasa Museum of History and Folklore.

  • 54

    Torihama Shell Midden Study Group (1987). Torihama shell midden – summary of research from 1980-1985 (In Japanese).Wakasa Town: Board of Fukui Prefecture and Wakasa Museum of History and Folklore.

  • 55

    UchiyamaJ. (2000). The problem of deer and wild boar in Torihama shell midden: The function of the site revealed by the remains of sika deer and wild boar excavated in 1984 (In Japanese).Torihama Shell Midden Pap.2122.

  • 56

    UchiyamaJ. (2005). Hunting seasonality of sika deer (Cervus nippon) and wild boar (Sus scrofa) at Torihama shell midden (Early Jomon): A zooarchaeological analysis of the faunal remains of the excavation in 1984 (In Japanese).Torihama Shell Midden Pap.4-579113.

  • 57

    WinklerD. E.Schulz-KornasE.KaiserT. M.CodronD.LeichliterJ.HummelJ.et al (2020). The turnover of dental microwear texture: Testing the “last supper” effect in small mammals in a controlled feeding experiment.Palaeogeogr. Palaeoclimatol. Palaeoecol.557:e1099930. 10.1016/j.palaeo.2020.109930

  • 58

    YamadaE.KuboO. M.KuboT.KohnoN. (2018). Three-dimensional tooth surface texture analysis on stall-fed and wild boars (Sus scrofa).PLoS One13:e0204719. 10.1371/journal.pone.0204719

  • 59

    YamakawaF. (1992a). Method of making bone piercing tools in the Jomon: Analysis of animal remains excavated from Torihama shell midden (In Japanese).J. Archaeol. Soc. Nippon7861106.

  • 60

    YamakawaF. (1992b). Bone piercing tools from Torihama shell midden (In Japanese).Bull. Wakasa Mus. Hist. Folk.4142.

  • 61

    YasudaY. (1979). “Pollen analysis,” in Torihama shell midden – research 1 on archaeological site located in lowland area from the early Jomon period (In Japanese), ed.Torihama Shell Midden Study Group (Wakasa Town: Board of Fukui Prefecture), 176196.

  • 62

    YasudaY. (1982). Pollen analytical study of the sediment from the Lake Mikata in Fukui Prefecture, Central Japan –especially on the fluctuation of precipitation since the Last Glacial age on the side of Sea of Japan (In Japanese with English abstract).Q. Res.21255271. 10.4116/jaqua.21.255

  • 63

    YokoyamaM. (2009). “Biology of sika deer in Hyogo: Characteristics of reproduction, food habits, growth, and condition,” in Sika deer, edsMcCulloughD. R.TakatsukiS.KajiK. (Tokyo: Springer), 193205. 10.1007/978-4-431-09429-6_14

  • 64

    YoshikawaM.YoshikawaJ.NoshiroS.KudoY.SasakiY.SuzukiM.et al (2016). Vegetation history and use of plant resources from the incipient to early Jomon periods at the Torihama shell midden site, Fukui Prefecture, central Japan (In Japanese with English abstract).Jpn. Assoc. Hist. Bot.246982. 10.34596/hisbot.24.2_69

  • 65

    YoshiokaK. (1973). Plant-geography (In Japanese).Tokyo: Kyouritsu-Shuppan.

Summary

Keywords

Archaeology, Jomon period, Torihama Shell Midden site, Sika deer, Microwear, ISO 25178

Citation

Sato K, Sato T and Kubo MO (2022) Reconstructing diets of hunted sika deer from Torihama Shell Midden site (ca. 6,000 years ago) by dental microwear texture analysis. Front. Ecol. Evol. 10:957038. doi: 10.3389/fevo.2022.957038

Received

30 May 2022

Accepted

16 September 2022

Published

10 October 2022

Volume

10 - 2022

Edited by

Ferran Estebaranz-Sánchez, Milá y Fontanals Institution for Research in Humanities (CSIC), Spain

Reviewed by

Florent Rivals, Catalan Institution for Research and Advanced Studies (ICREA), Spain; Lua Valenzuela-Suau, Universitat de les Illes Balears, Spain

Updates

Copyright

*Correspondence: Koyo Sato,

This article was submitted to Paleoecology, a section of the journal Frontiers in Ecology and Evolution

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics