INTRODUCTION
Organic residues that are preserved on archaeological pottery can yield a lot of information on the dietary habits of past populations. Furthermore, since most of these remains are the results of cooking or related activities, there is a clear chronological association with the period during which the pottery was used and by extension with the occupation of the site. Pottery residue is therefore often sampled for radiocarbon (14C) dating. There are several ways in which these residues can survive up to the present day, either being absorbed inside a vessel wall or encrusted on the internal and external pottery surfaces.
Most residue analyses have focused on the inside of pottery, where the presence of proteins and lipids indicates a clear association with the processing of animal or plant resources. These types of residues form attractive samples for 14C dating. However, dates on “food crusts” do not always correspond to dates obtained on other kinds of evidence. This led to discussions about the reliability of accelerator mass spectrometry (AMS) dates on food residues, especially where the processing of aquatic resources is expected (Fischer and Heinemeier Reference Fischer and Heinemeier2003; Craig Reference Craig2004; Hart and Lovis Reference Hart and Lovis2007; Boudin et al. Reference Boudin, Van Strydonck and Crombé2009, Reference Boudin, Van Strydonck, Crombé, De Clercq, van Dierendonck, Jongepier, Ervynck and Lentacker2010; Philippsen Reference Philippsen2013; Roper Reference Roper2013; Philippsen and Meadows Reference Philippsen and Meadows2014; Roper Reference Roper2014; Hart and Lovis Reference Hart and Lovis2014; Heron and Craig Reference Heron and Craig2015). Residues sticking to the outer pottery surface, on the other hand, are almost never analyzed, not even when they are used for 14C dating.
Outer Surface Residues
Carbonized material on the outside of pottery is commonly considered to be soot, derived from the combustion of wood. Apart from carbon, soot may be composed of various other products such as organic tars and resins. When resin droplets in the rising smoke come in contact with the cooler ceramic surface, they can solidify, forming a hard and waterproof layer (Skibo Reference Skibo2015).
Most research on soot has been conducted within the context of vessel use and use-alteration traces (Beck et al. Reference Beck, Skibo and Hally2002; Hally Reference Hally1983; Skibo Reference Skibo1992, Reference Skibo2013, Reference Skibo2015). Experiments by Hally (Reference Hally1983) and Skibo (Reference Skibo1992) point out that temperature is a key variable in the formation and patterning of soot, as it cannot form on surfaces approaching 400°C. Soot is therefore unlikely to be found on a vessel’s base, where it is in close contact with the fire. It will rather adhere to higher parts of the vessel, in particular when water inside the vessel is cooling down the ceramic surface.
The determination of residue as soot is usually based on visual inspection. A distinction is often made between thick, charred deposits or “food crusts” on the inside and a rather thin layer of black, sometimes glossy coating on the outside of pottery. In some cases, however, traces of food residue are observed on the external surface, mostly near the rim, as a result of cooking contents boiling over or even leaking through pores from the inside out (Andersen and Malmros Reference Andersen and Malmros1985:85; Koch Reference Koch1998:117; Evershed Reference Evershed2008:904). Even though food crusts may be present on the outer pottery surface, external residues are never analyzed prior to 14C dating. In fact, with few exceptions (e.g. Oudemans and Boon Reference Oudemans and Boon1991; Oudemans 2006), hardly any research has been done on the composition of these residues. As is the case with food crusts, 14C dates obtained on outer surface residues often do not correspond with the expected pottery age. Moreover, in many cases these dates also do not agree with those on food crusts coming from the same type of pottery or potsherd. These dating discrepancies are usually explained as the result of a reservoir effect on inner food crusts or of an old wood effect on outer soot (Fischer and Heinemeier Reference Fischer and Heinemeier2003; Miyata et al. Reference Miyata, Minami, Onbe, Sakamoto, Matzuzaki, Nakamura and Imamura2011; Philippsen Reference Philippsen2012, Reference Philippsen2013; Piezonka et al. Reference Piezonka, Meadows, Hartz, Kostyleva, Nedomolkina, Ivanishcheva, Kosorukova and Terberger2016).
This study confronts 14C dates on internal and external surface residues preserved on Mesolithic/Neolithic pottery from Bazel Sluis (Belgium) with the results of bulk isotope analysis and thermally assisted hydrolysis and methylation pyrolysis-gas chromatography-mass spectrometry (THM-GC-MS).
MATERIALS AND METHODS
Site and Samples
The prehistoric wetland site of Bazel Sluis is situated in northwestern Belgium, on the left bank of the lower Scheldt River (51°08′09′′N, 4°19′23′′E; Figure 1). It is located on the top and gentle slope of a former sand ridge, most likely a point bar. Apart from numerous burnt and unburnt animal bones and ecofacts, this site yielded thousands of pottery fragments and lithic artifacts dating from the Middle Mesolithic to the Middle Neolithic (Crombé et al. Reference Crombé, Sergant, Perdaen, Meylemans and Deforce2015a). The selected pottery fragments come from the eastern slope of the sand ridge, where the occupation levels have been covered by peat and fluvial sediments (Crombé et al. Reference Crombé, Verhegge, Deforce, Meylemans and Robinson2015b), offering favorable preservation conditions for plant material, unburnt animal bone, and organic pottery residues.

Figure 1 Palaeogeographical map of the Scheldt estuary around 6300 cal BP (modified after Vos and Van Heeringen Reference Vos and Van Heeringen1997) with indication of the prehistoric sites mentioned in the text: Bazel Sluis (B) and Doel Deurganckdok (D). Key: 1. Pleistocene coversand area; 2. peatland; 3. tidal area (mudflats and saltmarshes); 4. beaches and dunes; 5. North Sea, tidal inlets and tidal channels.
Eight sherds with residues on both the internal and external surfaces were selected for analysis and 14C dating. They represent rim, wall, and base fragments from vessels belonging to the same pottery group. Based on decorative elements (rim decoration, knobs, perforations) this pottery could be related to that of the Final Mesolithic-Neolithic Swifterbant culture, known from nearby sites at Doel “Deurganckdok” (Crombé et al. Reference Crombé, Boudin and Van Strydonk2008; Crombé Reference Crombé2010). At Doel, the 14C dates on food crusts and on organic material seem to be incompatible, the former clustering in the first half of the 6th millennium BP and the latter around the middle and in the second half of the 6th millennium BP (Boudin et al. Reference Boudin, Van Strydonck and Crombé2009, Reference Boudin, Van Strydonck, Crombé, De Clercq, van Dierendonck, Jongepier, Ervynck and Lentacker2010). The older 14C dates probably result from a freshwater reservoir effect on the food crusts, since these sites are located in the Scheldt River valley and cremated bones of freshwater fish were found during the excavations (Boudin 2009:697; Crombé et al. Reference Crombé, Boudin and Van Strydonk2008).
For the selected Bazel pottery, indirect 14C dates are available on unburnt animal bone originating from the same stratigraphic layer (Table 1) (Meylemans et al. Reference Meylemans, Perdaen, Sergant, Bastiaens, Crombé, Debruyne, Deforce, Du Rang, Ervynck, Lentacker, Storme and Van Neer2016). Based on these three dates, the pottery is expected to date between 5350 and 5000 BP or 4350 and 3700 cal BC. The calibration method is described in the section below.
Table 1 AMS 14C determinations on collagen from animal bones excavated at Bazel Sluis (Meylemans et al. Reference Meylemans, Perdaen, Sergant, Bastiaens, Crombé, Debruyne, Deforce, Du Rang, Ervynck, Lentacker, Storme and Van Neer2016).

The inner surface residues consist of thick layers of charred crusts, while the outer surfaces are patched with dusty black particles and/or a glossy black film or coating (Figure 2), in line with the descriptions of pottery residues from other archaeological sites. Visually, these residue types seem to be of different origins.

Figure 2 Organic residues preserved on the outer (left) and inner (right) surfaces of sherd 2/12/17/6 from Bazel Sluis.
AMS 14C Dating
All residues were pretreated with the acid-alkali-acid (AAA) method. 14C determinations were measured on the AMS at the Royal Institute for Cultural Heritage, Brussels (Lab code RICH-) (Boudin et al. Reference Boudin, Van Strydonck, van den Brande, Synal and Wacker2016). CO2 was obtained by sample combustion in the presence of CuO and Ag. Graphitization was performed with H2 over a Fe catalyst. Targets were prepared at the Royal Institute for Cultural Heritage in Brussels (Belgium) (Van Strydonck et al. 1990–Reference Van Strydonck and Van der Borg1991). 14C calibrations were performed using OxCal version 4.2 (Bronk Ramsey Reference Bronk Ramsey2009) and the IntCal13 calibration curve date (Reimer et al. Reference Reimer, Bard, Bayliss, Beck, Blackwell, Bronk Ramsey, Buck, Cheng, Edwards, Friedrich, Grootes, Guilderson, Haflidason, Hajdas, Hatté, Heaton, Hoffmann, Hogg, Hughen, Kaiser, Kromer, Manning, Niu, Reimer, Richards, Scott, Southon, Staff, Turney and van der Plicht2013).
Stable Carbon and Nitrogen Isotope Analysis
The δ13C and δ15N analyses were performed in duplicate on a Thermo Flash EA/HT elemental analyzer, coupled to a Thermo Delta V Advantage Isotope Ratio Mass Spectrometer via ConFlo IV interface (Thermo Fischer Scientific, Bremen, Germany). Standards used were IAEA-N1, IAEA-C6, and internally calibrated acetanilide. Analytical precision was 0.25‰ for both δ13C and δ15N based on multiple measurements of the standard acetanilide.
THM-GC-MS
In total 7 residue samples from both the inside and outside of different pottery fragments have been analyzed using THM-GC-MS. Prior to analysis, 20 μL 2.5wt% tetramethylammonium hydroxide (TMAH) in methanol solution was added to a tiny aliquot of sample (not weighed, but in the order of a few tens of micrograms). This solution also contained 100 ng/μL heptadecanoic acid (C17:0) as an internal standard. The content of this vial was well-mixed to homogenize, and 2 μL was transferred to the stainless-steel pyrolysis cup (Frontier Lab Eco-cup LF) with auto-Rx glass fiber disc. The cup was pyrolyzed at 480°C in a He atmosphere for 12 s.
Analyses were carried out with a Frontier Lab Multi-Shot Pyrolyzer (3030D) in conjunction with an autoshot sampler AS-1020ET (both Frontier Lab, Japan). The interface and the injector of the chromatographic system were kept at 300°C. The analytical column was directly coupled to the pyrolyzer via a custom made split device (split ratio 20), minimizing dead volume and improving the signal. For the chromatographic separations, a TraceGC gas chromatograph (Thermo), hyphenated with a PolarisQ Ion Trap mass spectrometer (Thermo), was used. Separations were accomplished on a SLB-5ms capillary column (Supelco, 20 m × 0.18 mm i.d.×0.18 μm film thickness) applying the following temperature program: initially the oven temperature was maintained at 35°C for 1 min after pyrolysis. Next, a 10°C/min gradient was applied up to 240°C; finally, the column was heated to a temperature of 315°C at a rate of 6°C/min; this temperature was maintained for 5 min. The carrier gas was helium 6.0 grade at a constant flow of 0.9 mL/min. The MS transfer line temperature was kept at 290°C. Ionization was carried out in the ion volume of the ion trap mass spectrometer under the standard EI positive mode at 70 eV. The mass spectrometer was scanned in the 35–650 atomic mass unit (amu) range, with a cycle time of 0.59 s.
RESULTS AND DISCUSSION
AMS 14C Dating
The 14C ages for charred residues on the inner and outer surfaces of the Bazel pottery are given in Table 2. Sample 8 did not contain enough residue on the outside to date.
Table 2 AMS 14C determinations on inner and outer surface residues preserved on pottery fragments from the prehistoric site of Bazel Sluis. When multiple 14C dates from the same sample are combined (prior to calibration) a χ2 test (Shennan Reference Shennan1988) is performed to establish their internal consistency. If the χ2 value (T) is higher than the 5% confidence limit of 3.8, the combined 14C dates are not coherent and the test fails. For the χ2 tests, the uncalibrated 14C dates on both types of residue from each sample were combined using OxCal v. 4.2 (Bronk Ramsey Reference Bronk Ramsey2009).

The outer surface residues systematically date younger than the food crusts. Chi-squared (χ2) tests indicate a statistically significant difference in most samples for which both residues have been dated. Only samples 4 and 7 pass the test. The extent to which the 14C dates on food crusts and outer surface residues differ seems to vary quite a lot between samples, with a discrepancy of up to approximately 400 14C yr for sample 3.
The 14C dates from organic residues preserved on pottery from other archaeological sites show similar deviations between inner and outer surface residues (Table 3). In all cases but one—a recently published date on the outer crust of Narva pottery at Veksa 3 (Piezonka et al. Reference Piezonka, Meadows, Hartz, Kostyleva, Nedomolkina, Ivanishcheva, Kosorukova and Terberger2016)—residues on the outer surfaces yield younger ages. Again there seems to be a considerable variation as to how much these ages differ from those obtained on food crusts. Sample SLAS-1713, a potsherd from the Ertebølle site of Schlamersdorf in northern Germany, probably is not representative here. This sample presents an exceptionally large dating discrepancy that is difficult to explain. According to Philippsen (Reference Philippsen2013), a freshwater reservoir effect is however likely for the food crusts on the inside of this pottery.
Table 3 AMS 14C determinations on paired organic residue preserved on pottery from archaeological sites in Denmark, Germany, Russia, and Japan.

Fischer and Heinemeier (Reference Fischer and Heinemeier2003) and Miyata et al. (Reference Miyata, Minami, Onbe, Sakamoto, Matzuzaki, Nakamura and Imamura2011) find reason for the dating discrepancies in a (freshwater) reservoir effect on food crusts—in the light of younger 14C dates on associated organic material—and/or in an old wood effect on the outer surface soot. Miyata et al. (Reference Miyata, Minami, Onbe, Sakamoto, Matzuzaki, Nakamura and Imamura2011), in a recent article about this same dating problem, mention diagenesis as a third possible cause. According to these authors, the less carbonized food crusts on the inner surfaces could be more vulnerable to the influence of humic acids in the soil. However, a harsh AAA treatment prior to AMS 14C dating, as was done for this study, should remove all humic acids. Piezonka et al. (Reference Piezonka, Meadows, Hartz, Kostyleva, Nedomolkina, Ivanishcheva, Kosorukova and Terberger2016) emphasize that there is currently no direct evidence for a significant freshwater reservoir effect at Veksa 3, but that high δ15N values for both inner and outer surface residues on sample VE-2007/118 suggest the use of aquatic resources. Regardless, on this site there seems to be no difference between the 14C dates on both types of residue.
At Bazel Sluis, not only do the outer surface residues date systematically younger than the food crusts, but both types of residue yield 14C dates that are too old in comparison to dates obtained on associated animal bone. An old wood effect on the outer residues in combination with a (larger) freshwater reservoir effect on the inner residues is a possible explanation for this discrepancy. However, in order to gain insight into the reasons for this age difference we need to look at the biochemical composition of both residue types.
Stable Carbon and Nitrogen Isotope Analysis
The δ13C, δ15N values, and atomic C:N determined for charred residues on the inner and outer surfaces of the Bazel pottery are given in Table 4 and plotted in Figure 3. Samples 1 and 8 did not contain enough residue on the outside to analyze.

Figure 3 Plot of the δ13C and δ15N values determined for charred residues on the inner and outer surfaces of the Bazel pottery. The analytical precision is 0.25% for both δ13C and δ15N.
Table 4 δ13C, δ15N values and atomic C:N ratios determined for charred residues on the inner and outer surfaces of the Bazel pottery. The right-most column presents the analysis numbers for residues analyzed with THM-GC-MS at the Royal Institute for Cultural Heritage (Brussels, Belgium).

The isotopic values of inner and outer surface residues cluster within the same range, between approximately –26.5 and –28.5% δ13C and between approximately +6 and +10% δ15N. The only exceptions are the inner residues of sample 2 in nitrogen values and of sample 1 in carbon values. Unfortunately, for the latter sample we do not have comparative data for the external surface. Sample 3 presents the biggest difference in carbon values between the inner and outer surface residues.
Looking at these carbon and nitrogen values, it could well be that aquatic resources have been processed inside these vessels. Craig et al. (Reference Craig, Forster, Andersen, Koch, Crombé, Milner, Stern, Bailey and Heron2007) defined a bulk isotopic theory to distinguish aquatic products in charred food residue, according to which relatively light δ13C values (<–25%) and δ15N values between +6 and +10% indicate that freshwater fish is a likely component of these residues, and probably the major component for samples with high δ15N values (+8 to +10%). At Bazel, both burnt and unburnt freshwater fish bones were collected during the excavation, but it is not clear which of these remains are related to the same occupation phase as the pottery fragments discussed in this article. However, the presence of burnt fish bones and the location of this site near the Scheldt River, in combination with 14C dates on food crusts that are too old, can be considered as indicators for the processing of freshwater fish.
If we apply Craig’s bulk isotopic theory on the Bazel pottery, freshwater fish is likely to be present in all residues, including those from the external surfaces. Regardless of the presence of aquatic resources, the clustering of isotopic values from both types of residue indicates that they at least have some similar components.
Finally, it is noted that the outer surface residues systematically have higher C:N ratios than the inner food crusts. This has also been noticed by Miyata et al. (Reference Miyata, Minami, Onbe, Sakamoto, Matzuzaki, Nakamura and Imamura2011) for the Jomon pottery from the Irienaiko site in Japan and by Philippsen (Reference Philippsen2012) for residues that had formed on copies of prehistoric vessels during cooking experiments. Based on the isotopic ratios and the C:N values, the latter author suspected a possible presence of food residue in the outer crusts, as a result of food boiling over. Philippsen (Reference Philippsen2012:125) considered the higher C:N ratios for the outer crusts as an indication that not only food residue but also soot was present on the exterior.
THM-GC-MS
To determine their biochemical composition, 4 samples of inner surface residue and 3 samples of outer surface residue have been analyzed using THM-GC-MS (Table 4). For both types of residue similar chromatograms were obtained (Figure 4), stressing the results of the bulk stable isotope analyses, i.e. that many of the same organic compounds can be found on the inside and outside of the Bazel pottery.

Figure 4 Gas chromatograms (THM-GC-MS) of the inner and outer surface residues of sample 4 from Bazel Sluis. “B” stands for compounds also present in blank samples. “Cx:y” indicates a fatty acid with (x-axis) the number of carbon atoms and (y-axis) the number of unsaturated bonds. A solution containing C17:0 as an internal standard was added to each sample.
Both pottery surfaces show clear indicators for the presence of food residue. Many fatty acids—mainly palmitic (C16:0) and stearic (C18:0)—are measured, with oleic acid (C18:1) as the most common unsaturated fatty acid. It is noted that the ratio between palmitic (P) and stearic (S) acids is more variable within the outer surface residues (P:S 0.47–1.2) as compared to the inner surface residues (P:S 0.52–0.57). With the exception of the outer crust on sample 3, all residues seem to contain cholesterol-related compounds. In addition, all residues show at least some protein fragments. On the inside of sample 5 and the outside of sample 7 palmitone, the ketone of palmitic acid, is detected. Although some of these compounds may derive from plant material, several elements point toward the presence of animal fats and the processing of animal resources inside this pottery. It is not yet clear whether these are related to aquatic or terrestrial animal species.
In all residues one or several polycyclic aromatic hydrocarbons (PAHs) are detected. It concerns benzo[e]pyrene or benzo[k]fluorantheneFootnote 1 , possibly in combination with indeno[1,2,3-cd]pyrene for the outer crust of sample 3. These PAHs derive from the incomplete combustion of organic compounds, probably animal fats or fuel (wood?). A study by Oudemans and Boon (Reference Oudemans and Boon1991) also revealed the presence of PAHs in residues on the outside of pottery, using Py-GC-MS, which was interpreted as being indicative for smoke condensates or “soot” from cooking over wood fires (Oudemans 2006:38). In order to test this premise, a reference sample of modern soot was analyzed by using the same technique. This led to the detection of several markers (in the form of PAHs) for soot mentioned in the literature (Wei et al. Reference Wei, Fang, Cao and Schreiner2011). However, inside the pottery residues from Bazel only one of these markers, benzo[e]pyrene or benzo[k]fluoranthene, was detected. It is possible that the analyzed samples were too small to detect more PAHs, or maybe other techniques, such as Raman spectroscopy (Coccato et al. Reference Coccato, Jehlicka, Moens and Vandenabeele2015), would be more successful in identifying soot. It seems unlikely that soot is totally absent from the ceramic surface, especially on the outside of the vessels.
It is clear from the analyses that food remains are not only present on the inner pottery surface but also on the outer pottery surface. Most of these remains are located near the rim, but they also seem to occur on wall and even base fragments. During cooking inside these vessels, foodstuff probably boiled over, or organic compounds migrated through the vessel wall and burnt on the exterior (Evershed Reference Evershed2008). The possible presence of food residue on the outside of archaeological pottery is an important factor to take into account when sampling material for 14C dating, as it may be subject to the same dating problems as inner food crusts, i.e. a reservoir effect.
However, if only food remains were present on either side of the vessel then there would not be a divergence between their 14C dates. We suggest that this discrepancy results from the mixed presence of food residue and soot on the external pottery surfaces, as opposed to a more homogeneous composition of food residues encrusted on the inside of this pottery. Because the (freshwater) reservoir effect is probably much larger than a possible old wood effect, sampling a mixture of soot and food crusts may lead to younger dates on the external pottery surfaces. Because the residues on the external surfaces are mixed, variations in the age differences of the residues between different pottery fragments could be related to sampling location and quantity, affecting the relative proportions of food and soot present within the dated residue. For example, on the external surface of sample 7 it is theoretically possible that mostly food residue was preserved or sampled, resulting in little or no age difference between the inside and outside. Sample 3, by contrast, presents the largest discrepancy, possibly indicating that the residue sampled on the external surface mostly consisted of soot.
CONCLUSION
This study demonstrates the presence of food residue on the external surfaces of pottery from Bazel Sluis. Using THM-GC-MS it was possible to detect many fatty acids, cholesterol-related compounds and protein fragments, proving that the outer surface crust does not entirely consist of soot. This means that outer surface residues might also be subject to a reservoir age. Chemical analyses (e.g. THM-GC-MS) should be conducted prior to 14C dating to determine the composition of these residues.
THM-GC-MS has also led to the detection of benzo[e]pyrene/benzo[k]fluoranthene, a possible soot marker. However, not enough markers (PAHs) were detected to prove soot, although its presence can at least be presumed for the outer pottery surfaces. Other techniques, such as Raman spectroscopy, might be better suited to identify soot.
14C dates on outer surface residues at Bazel and also on other sites are often younger than those obtained on inner food crusts, but still are not compatible with dates obtained on associated plant and/or bone samples. We suggest that (1) these food crusts are subject to a reservoir effect and (2) the younger dates on the outer surface residues are due to a mixed presence of soot and food residue. Because the reservoir effect is larger than a possible old wood effect on soot, the mixture of both types of residue on the outside leads to younger dates than those obtained on more homogeneous food crusts sampled from the internal pottery surface.
ACKNOWLEDGMENTS
This research is a collaboration between the department of Archaeology of Ghent University and the Royal Institute for Cultural Heritage. The authors wish to thank the Flanders Heritage Agency for providing the sampled material. Part of this research was funded by the “Special Research Fund” (BOF) of Ghent University, in the context of the research project “Agropastoral Neolithic Ways of Life in the Sandy Lowland of Belgium.”
 
 







