Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Evaluating the Photoprotective Effects of Ochre on Human Skin by In Vivo SPF Assessment: Implications for Human Evolution, Adaptation and Dispersal

  • Riaan F. Rifkin ,

    riaanrifkin@gmail.com

    Affiliations Institute for Archaeology, History, Culture and Religion, University of Bergen, Bergen, Norway, Evolutionary Studies Institute, University of the Witwatersrand, Johannesburg, South Africa

  • Laure Dayet,

    Affiliation Centre National de la Recherche Scientifique, Unité Mixte de Recherche, University of Bordeaux, Pessac, France

  • Alain Queffelec,

    Affiliation Centre National de la Recherche Scientifique, Unité Mixte de Recherche, University of Bordeaux, Pessac, France

  • Beverley Summers,

    Affiliation Photobiology Laboratory, Department of Pharmacy, University of Limpopo, Medunsa, South Africa

  • Marlize Lategan,

    Affiliation Photobiology Laboratory, Department of Pharmacy, University of Limpopo, Medunsa, South Africa

  • Francesco d’Errico

    Affiliations Centre National de la Recherche Scientifique, Unité Mixte de Recherche, University of Bordeaux, Pessac, France, Evolutionary Studies Institute, University of the Witwatersrand, Johannesburg, South Africa

Abstract

Archaeological indicators of cognitively modern behaviour become increasingly prevalent during the African Middle Stone Age (MSA). Although the exploitation of ochre is viewed as a key feature of the emergence of modern human behaviour, the uses to which ochre and ochre-based mixtures were put remain ambiguous. Here we present the results of an experimental study exploring the efficacy of ochre as a topical photoprotective compound. This is achieved through the in vivo calculation of the sun protection factor (SPF) values of ochre samples obtained from Ovahimba women (Kunene Region, Northern Namibia) and the Palaeozoic Bokkeveld Group deposits of the Cape Supergroup (Western Cape Province, South Africa). We employ visible spectroscopy, energy-dispersive X-ray fluorescence (ED-XRF), X-ray diffraction (XRD) and granulometric analyses to characterise ochre samples. The capacity of ochre to inhibit the susceptibility of humans to the harmful effects of exposure to ultraviolet radiation (UVR) is confirmed and the mechanisms implicated in the efficacy of ochre as a sunscreen identified. It is posited that the habitual application of ochre may have represented a crucial innovation for MSA humans by limiting the adverse effects of ultraviolet exposure. This may have facilitated the colonisation of geographic regions largely unfavourable to the constitutive skin colour of newly arriving populations.

Introduction

Climate exerts a significant influence on ecosystems, communities and populations [14]. Of the climate-driven selective pressures that operated during the evolutionary history of the genus Homo, negotiating the risks and benefits of persistent exposure to sunlight presented an enduring challenge. Sunlight is an essential environmental factor in nearly all ecosystems [5, 6]. Most important in terms of human health are long-wave ultraviolet A (UVA) (315–400 nm) and medium-wave UVB (280–315 nm) radiation [7]. Whereas UVA causes DNA damage, skin aging and mild erythema in large doses, UVB is responsible for sunburn with subsequent DNA damage and skin cancers. Short-wave UVC (200–280 nm) is the most dangerous type of UVR but is largely absorbed by the ozone layer [6].

The beneficial effects of moderate sunlight exposure are numerous and several positive correlations exist between adequate UVR exposure, vitamin D synthesis and human fertility [8, 9, 10]. Vitamin D (the active form being 1.25 dihydroxyvitamin D3 or 1.25[OH]2D3) in turn reduces the incidence of rheumatoid arthritis, coronary heart disease, diabetes, multiple sclerosis, rickets, schizophrenia, autoimmune disease and several types of cancer [1114]. Conversely, excessive UVR exposure frequently causes malignant skin diseases [11, 15] and mitochondrial (mtDNA) and nuclear (nDNA) damage [1618]. Quantitative (mtDNA copy number) and qualitative (mutations) alterations in mtDNA in turn result in neurodegenerative and metabolic diseases [19] and amplifies immune suppression [20]. Although UVR exposure is the foremost cause of skin cancers, it is also a primary source of vitamin D in light-skinned individuals. As vitamin D protects against cancer and many other diseases, there is some disagreement concerning the health benefits and risks associated with UVR [8].

It has also been proposed that excessive UVR may have had influenced evolutionary processes, including extinction events such as the disappearance of megafauna during the late Pleistocene [2123] and the extinction of some hominin species [24, 25]. It has been demonstrated that there is a correlation between hominin extinction events and the 41 ka orbital cycle attributable to change in our planets obliquity [25]. The first peaked 1.5 million years ago (mya) and closely followed the Olduvai (1.67 to 1.95 mya) and Gilsa (1.57 mya) geomagnetic excursions [26]. This event appears to coincide with the extinction of Paranthropus robustus, P. bosei, Homo ergaster, H. rudolfensis and H. habilis [25]. A second pulse of extinctions, affecting H. neanderthalensis and H. floresiensis, occurred during the last 120 ka and corresponds loosely with the Blake (138 ka to 100 ka) and Laschamp excursions [2729]. The Laschamp excursion occurred at ~ 40 ka and has been cited as the most dramatic geomagnetic event to coincide with the extinction of H. neanderthalensis [24, 25, 30, 31]. Following the variable climate experienced during Marine Isotope Stage (MIS) 3 between 60 ka and 40 ka, the climate began to cool and by 30 ka Neanderthals were extinct. Proposed reasons for their demise include climate change, catastrophic natural events, disease epidemics, food shortages and the inability to compete with H. sapiens or other predatory species [3239]. Because Neanderthals have been shown to have had lightly pigmented skin [40, 41], it has been suggested that increasing exposure to UVR may have been a factor in their disappearance [24, 25, 42].

These hypotheses are difficult to test because of the uncertainties concerning the dating of extinctions events and in calculating the impact of UVR changes at different latitudes and among different populations. Moreover, the absence of the MC1R gene variant (R307G) in Neanderthal and Denisovan genomes suggest that considerable variability in skin colour existed among archaic hominins, and that some Neanderthal and Denisovan groups had darker and therefore more UVR resistant skin tones [43, 44]. Fluctuations in UVR intensity may nevertheless have generated a selective pressure on human populations that could have affected the life expectancy of portions of those populations.

During the Laschamp excursion, the dipole geomagnetic moment (the intensity of the magnetic field currently expressed as 7.78 x 1022 A.m2), decreased by a factor of ~4 (from 10.5 x 1022 A.m2 at 49 ka to 2.6 x 1022 A.m2 at 41 ka), resulting in a twofold increase in cosmogenic nuclide deposition rates [45, 46]. At the same time, the amount of solar-derived UVR reaching earth’s surface increased by >25% [47, 48]. To put this surge in UVR into perspective, UVA and UVB exposure in sub-Saharan Africa currently exceeds 250 watts (W) per active area (m2) and rarely falls below 175 W/m2 [14]. UVR exposure levels >250 W/m2 present extreme exposure categories and are rated 10 (defined as posing a ‘very high risk of harm from unprotected sun exposure’) to 11+ (‘extreme risk of harm from unprotected sun exposure’) [49]. Exposure levels raised by 25%, such as that brought on by the Laschamp event (at >300 W/m2), would therefore have posed a sudden and severe health risk.

Strategies for dealing with increasing UVR exposure may have comprised both biological adaptations and behavioural advances [6, 50]. The use of complex forms of clothing may have presented one option in terms of limiting the adverse effects of excessive UVR. The ability to transform animal hide into leather is one of the most important technological skills mastered by humans [34] and genetic research suggests that this may have occurred between 107 ka and 72 ka [51, 52]. In lower latitudes characterised by high temperatures, clothing may have presented a disadvantage, such as inducing increasing perspiration and diminishing thermal insulation due to evaporative cooling. In such regions, the topical application of ochre may have presented an important novel behavioural innovation. Resembling the positive impact of technological innovation in enhancing human subsistence strategies [5356], the topical application of ochre may therefore have served a significant function in terms of limiting the adverse effects of UVR exposure.

In this paper we recall the physiological mechanisms involved in the UVR protection and skin pigmentation, identify viable ethnographic candidates for the photoprotective use of red ochre, summarize the archaeological evidence attesting such a use in the past and, more importantly, experimentally test the efficacy of ochre as a topical photoprotective agent. As a result, we are able to demonstrate that certain types of iron-rich minerals, with or without binders, provide a perceptible degree of protection against UVR. We also evaluate the hypothesis that ochre may have been used as a sunscreen in the MSA and the Later Stone Age (LSA), discuss the advantages that such an innovation may have conveyed to our ancestors and propose a tentative scenario for its origin and geographic spread.

UVR Protective Mechanisms

Human skin derives most of its pigmentation from two types of melanin, a brownish-black eumelanin and a reddish-yellow pheomelanin. Whereas constitutive skin colour is the genetically determined cutaneous melanin pigmentation generated without solar influence, facultative skin colour is the short-lived tanning reaction elicited by exposure to UVR [57]. Many theories have been proposed to explain the variation in human skin pigmentation. For the development of skin darkening, examples include the shielding of sweat glands and blood vessels [58], protection from cancer and the overproduction of vitamin D [59], guarding against folate photodestruction [60] and defense against microorganisms [61] and infection [62]. Hypotheses for skin lightening comprise sexual selection, adaptation to cold UVR deficient climates and the enhancement of vitamin D photoproduction [12]. An additional evolutionarily significant effect of UVR involves the metabolism of folate, a B vitamin essential for the modulation of gene expression and the production and repair of DNA [63]. Acute folate deficiencies result in male infertility and the increased prevalence of birth defects [64, 65]. Maintaining adequate levels of folate is essential for human reproduction, and the evolution of dark skin pigmentation was an effective means to achieve this [57, 66].

The most commonly accepted factor explaining the geographic distribution of human skin pigmentation is UVR exposure [63, 67]. Evidence for natural selection operating at low latitudes (establishing and maintaining dark pigmentation under high UVR conditions), and high latitudes (the development of light pigmentation under low UVR conditions) confirms that human skin coloration is a Darwinian adaptation, with the gradient of skin colours observed from the equator to the poles being the product of two clines operating over a spatially varying optimum of UVR distribution [68, 69]. The correlation between pigmentation in indigenous populations and latitude is therefore traceable to the correlation between skin colour and the intensity of UVR exposure [60, 65, 7076] (Table 1).

thumbnail
Table 1. The classification of human skin types pertaining to responses to UVR exposure [7075].

Constitutive values in brackets refer to the Von Luschan chromatic skin colour classification scale [76]. ‘Exposure’ refers to the amount of time during which UVR exposure would not induce perceptible damage to the skin. RSF is a factor characterizing the protective effect of a sunscreen against the generation of free radicals.

https://doi.org/10.1371/journal.pone.0136090.t001

Differences in skin pigmentation do not result from differences in the number of melanocytes in the skin, but from differences in the melanogenic activity, the type of melanin produced in melanosomes and the size, number and packaging of melanosomes, with melanin content of melanosomes ranging from 17.9% to 72.3%. While melanin content in the epidermis differs by two-fold in Asian, and white skin, dark skin contains higher levels of melanin (3 to 6 fold), has more eumelanin and larger melanosomes than white skin. As a result, and whereas a dark epidermis allows only 7.4% of UVR-B and 17.5% of UVR-A to penetrate the skin, 24% UVR-B and 55% UVR-A passes through white skin [77].

That darker pigmentation conferred significant adaptive benefits is shown by the conservation of the MC1R, SLC45A2 and 70 other genetic loci involved in pigment production [7880]. This indicates that there was strong evolutionary pressure to retain pigmentation among sub-Saharan Africa populations [62, 65]. Unlike SLC45A2, which occurs only amongst Europeans, SLC24A5, TYRP1 and KITLG are present at low frequencies in some sub-Saharan populations [80]. While these alleles may have arisen in and spread out of Africa with the early migrations of modern humans, many may derive from a series of admixture events with groups containing Eurasian genetic ancestry [81].

Ethnographic Evidence

Ethnographic accounts illustrating the use of red ochre as a cosmetic are abundant and have been reported for San hunter-gatherers [8289], Khoe pastoralists [90], and Tswana [91] and Xhosa [92] agro-pastoralists. Foremost examples of the habitual use of red ochre as a body cosmetic comprise the Cushitic-speaking Hamar in Southern Ethiopia [93, 94] and the Ovahimba of Angola and Namibia [95100]. The Ovahimba form part of a complex socio-economic assemblage of hunter-pastoralists who inhabit this region. They presumably represent a former Herero community who remained in the Kunene Region on a southward migration from the Central Lakes region into Namibia, or the final settlers in the terminal destination on a westward migration into the area after AD 1500 [99] (Fig 1). Recent archaeological excavations indicate that humans have inhabited the Kunene Region from at least 220 ka [101].

thumbnail
Fig 1. The location of the Kunene Region (red shaded area) in north-western Namibia (sourced from https://en.wikipedia.org/wiki/Template:Location_map_Namibia/doc).

https://doi.org/10.1371/journal.pone.0136090.g001

Ovahimba women are renowned for covering their bodies, hair and personal attire with a red ochre-based substance, ‘otjise’, which consists of equal parts of milk-derived clarified butter and red ochre powder [100]. Alexander (1838) observed that both women and men are ‘fond of greasing the skin’ with otjise [102], and personal observations and ethnographic interviews conducted by two of us (RFR, FDE) confirm this statement. The habitual application of otjise by men only declined during the 1960s, largely as a result of the presence of the South African Defense Force in the region and the subsequent employment of some Ovahimba as soldiers (Hepuite Venjakera, Nakara Tjimbosi and Tharhirwa Shawari pers. comm. 2012 [103, 104].

Clarified butter is produced by women and red ochre is mined or bought from Zemba land-owners or travelling merchants. While ochre processing occurs arbitrarily and forms part of the general sociable settings of daily life, the application of otjise occurs within the confines of women’s huts. Post menarche women mix roughly equal amounts of ochre powder and clarified butter together between the palms and apply 4.2mg of otjise per 1cm2 of skin surface. Women therefore require 60g of red ochre powder to cover their bodies (every 2 to 3 days) and 320g to treat their hair (every 2 to 4 weeks) (Fig 2) [100]. Otjise features in initiation ceremonies, is applied by men when they are to be wed or when they undertake extensive journeys and used to cover human corpses prior to interment [95, 97]. Personal belongings such as leather headdresses and aprons, wooden containers and jewellery are also covered with the compound, apparently for both symbolic and functional (preservative) purposes. Whereas infants are intentionally covered with the mixture, regular contact between women and children result in the substantial transfer of the red compound to toddlers. Ethnographic interviews conducted by two of us (RFR, FDE) (refer to http://www.canal-u.tv/video/cerimes/the_color_of_the_ovahimba.15970) revealed that, besides the intrinsic symbolic significance of red ochre, it also fulfills several secondary functional roles, including that as a topical sun-protection element.

thumbnail
Fig 2. The production and application of otjise.

(A) Red ochre powder is obtained by grinding chunks (inset—scale bar is 1cm) between a round upper and a flat lower grindstone, (B) after which it is mixed (at a 1:1 ratio between the palms) with milk-derived clarified butter and (C) applied to the hair, body and ornaments.

https://doi.org/10.1371/journal.pone.0136090.g002

Archaeological Evidence

The exploitation of ochre of variable colour and geological origin is not restricted to Homo sapiens [105107] and it does not reflect a sudden change signifying a shift away from archaic to modern and symbolic patterns of behaviour. Red ochre residues have been found inside perforated shell beads from African and Near Eastern MSA and Levantine Mousterian sites dated from 92 ka to 60 ka [108112]. Since strung beads were in all probability worn around the neck or wrists [113], these residues most likely transferred from human bodies onto the shells. But as direct evidence on how and why ochre was utilised during the MSA and LSA is not available, diverse interpretations for its usage have been proposed. Although the symbolic and functional significance of ochre is a hotly debated topic [114123], these uses are not necessarily mutually exclusive.

Cultural innovations seeking to reduce the impact of UVR, such as the habitual topical application of powdered ochre, could also have reduced the selective pressure induced by increasing exposure to UVR. Changes in climate and ecology have been shown to provoke speciation and extinction [25, 124127] and it has been proposed that these may have acted as selective pressures for enhanced human cognition [128131]. MIS 5e to 5d represent the most dynamic periods of innovation and may have generated the conditions for the expansion of populations within and out of Africa [4]. Examples of geographically confined cultural traditions that may represent the first instances of socially mediated symbolic responses to changing ecological conditions, include abstract engraved designs on ochre [132135] and perforated marine shell beads [106, 108113]. It has been shown that cultural niche construction [136] played a key role in positively influencing human evolutionary outcomes [137, 138]. Given that dark-skinned individuals are also subject to the harmful effects of UVR [68, 69], the development of effective sun-protection strategies was essential. One solution may have entailed the topical application of organic or inorganic agents perceived to reduce the detrimental effects of UVR.

Materials and Methods

Permission to carry out the Namibian ethnographic research was provided by the Namibian Ministry of Culture and the Namibian Filming Authority. All individuals referred to and represented in this manuscript have given written informed consent to publish case details, including data and photographs. Permission to collect ochre samples from the Palaeozoic Bokkeveld Group deposits was provided by private landowners. Informed consent for in vivo SPF experiments was obtained in written form and approved by the Photobiology Laboratory of Medunsa Campus of the University of Limpopo (Sefako Makgatho University of Health Sciences) Research and Ethics Committee (Project Number MREC/H/158/2014). One of the authors (RFR) participated in the in vivo SPF trials.

Ochre

The in vivo experimental results reported here follows from an initial in vitro SPF assessment of 24 ochre samples derived from ethnographic sources and geological resources likely to have been exploited during the MSA [100, 139]. Samples one to six derive from Okamanga, Ovinjange and Otjongoro villages in the Kunene Region of Northern Namibia. These were acquired by Ovahimba women in Opuwo and derive from Ombumbuu red ochre mines near Ruacana. Red ochre derived this geological source appears to be particularly sought after by Ovahimba women across the Kunene Region. The samples were processed into powder by Ovahimba women by crushing and grinding ochre chunks between a round upper and a flat lower grindstone. Samples 7 to 12 were acquired in Opuwo and were experimentally processed by direct grinding onto a coarse quartzite surface. This method resembles the techniques used to produce ochre powder during the MSA [121]. An additional twelve samples, comprising yellow, grey and red specimens, were collected from the Palaeozoic Bokkeveld Group deposits of the Cape Supergroup [140]. These are mostly fissile shales and comprise the most likely source of ochre at Blombos Cave [133]. Six samples were processed into powder by direct grinding onto a coarse quartzitic sandstone surface, and six by way of the technique employed by Ovahimba women. Table 2 presents information concerning the samples subjected to UVR during the experiments.

thumbnail
Table 2. Experimental ochre samples subjected to SPF, ED-XRF, XRD, granulometric, colorimetric and UVR analyses.

https://doi.org/10.1371/journal.pone.0136090.t002

Binders

For the SPF assessments, ochre samples were mixed (at a 1:1 ratio) with two types of organic binders. Clarified butter was procured from Otjongoro village and subcutaneous fat from a kudu (Tragelaphus strepsiceros) antelope. As per the SANS 1557/ISO 24444 in vivo SPF testing protocol [6], we applied 2mg of experimental materials per 1cm2 of skin surface during the experiments. Although Ovahimba women typically apply 4.2mg otjise per cm2, we adhered to the standard 2mg/cm2 ratio to acquire comparable SPF values. Given that Ovahimba women apply an additional 2.2mg per cm2, the actual SPF values are likely to be significantly higher than the results presented here.

Analytical methods

Ultraviolet radiation (UVR).

The assessment of the efficiency of a sunscreen is based on the value of its ‘sun protection factor’ (SPF), which reflects the degree of protection against UVR-induced erythema. At a certain exposure dose of UVB, minimally detectable redness can be observed. This is referred to as the minimal erythema dose (MED). SPF is the ratio of MED after sunscreen application and baseline (untreated) MED. Whereas SPF values generally denote the efficiency of a sunscreen to protect the skin from medium-wave UVB (280–315 nm) radiation, UVA is also a component of the solar simulated radiation which forms part of the in vivo SPF test. UVAPF denotes its UVA or long-wave ultraviolet (315–400 nm) protection factor. A sunscreen with an SPF of 2 filters out 50% of UV, an SPF of 15 filters out 93% UV and an SPF of 50 filters out 98% UV [6, 141]. In theory, the application of a product with SPF of 5 provides sun protection for five times longer than unprotected skin. Modern sunscreens typically comprise UVR absorbing chemicals that attenuate the amount of UVR reaching viable cells in the skin. Sunscreens contain chemical filters that absorb UVB radiation and physical filters that attenuate UVA and UVB radiation through reflection and scattering [6, 50]. Whereas organic molecules are selected for their UVR absorption capacity (cinnamates, benzophenone and butylmethoxy dibenzoyl methane), inorganic particles (zinc, titanium and iron oxides) are selected for their capacity to absorb, reflect or scatter UVR [142]. Sunscreens do not completely prevent photodamage but they do alter the UVR spectrum that reaches skin cells [7].

The UVR protection capacities (SPF and UVAPF values) of ochre samples were established by means of a series of initial in vitro and in vivo experiments performed under controlled conditions at the Photobiology Laboratory, Medunsa Campus (Sefako Makgatho University), South Africa [143]. As an initial screening method, experimental samples were subjected to in vitro protection tests [100]. The samples were applied to Transpore tape at a ratio of 2mg/cm2 and analysed according to the SANS 1557:1992 procedure. Critical Wavelength (CW) was determined using the Optometrics SPF 290 method. This method was chosen as there is unlikely to be photodegradation of the samples, hence the more recent ISO/SANS 24443 were not considered necessary. The CW is the wavelength below which 90% of the UV protection is situated. Accordingly, the higher the CW, the higher the UVA protection ranges of the products. Following the initial in vitro experiments, the six ochre samples with the highest SPF values (red ochre samples 1, 5, 7, 10, 16 and 18) were selected for in vivo SPF assessment. Informed consent was obtained from all research participants and the principles of the Declaration of Helsinki were strictly adhered to. The in vivo SPF test method followed the SANS 1557/ISO 24444 SPF test protocol. Tests were performed on three test subjects, designated A, B and C, with skin phototype II, and conducted using an application quantity of 2mg/cm2. The UV irradiation sources consisted of a calibrated Multiport Solar Simulator. An SPF 15 reference standard was included in the test and erythema was assessed visually by two experienced assessors 24 hours after irradiation.

Visible spectroscopy.

To determine whether colour played a role in the UVR reflectance and absorption capabilities of ochre, visible spectroscopy was employed to obtain L*a*b* values for each sample. We used an Avantes AvaSpec2048 fibre optic spectrometer equipped with a 2048 pixel CCD detector set to operate in the retrodiffusion mode. This instrument is equipped with an optical fibre probe, which is set in contact with the powder samples stored in transparent plastic bags. An AvaLight-HAL was used as an illumination source. The equipment is calibrated with a Halon D65 white reference sample in the same lighting conditions as for the archaeological samples. The colour parameters were obtained by Avasoft 7.5 software.

Particle size analyses.

Since particle size distribution may play an important role in determining the UVR protection capacities of clay minerals [144], ochre samples were subjected to particle size analyses. We used a Horiba LA950 (0.01μm—3000μm) laser micro-granulometer and analysed the data with NextSpec 7.10 software. The Mie solution to Maxwell’s equations provides the basis for measuring the size of particles through the scattering of electromagnetic radiation [145]. We employed this solution to calculate the particle size distribution in aqueous solution (refractive index 1,333). For sediments not predominantly comprised of iron oxides and hydroxides, refractive indices for the Mie calculation were based on the standard sediment refractive index of 1.55i–0.01i. For the samples consisting primarily of hematite, calculations were made with a refractive index of 2.94i–0.01i. The pre-treatment of samples included suspension in sodium hexametaphosphate (5 g/L) for 12 hours and 60 seconds of ultrasonification in the micro-granulometer to achieve optimal dispersion. The results are expressed through the use of D10, D25 and D50 values, the size of the particles that represent the 10th, 25th and the 50th percentile of the cumulative granulometric distribution. A sample with a D10 value of 15μm is a sample in which 10% of the particles are smaller or equal to 15μm. For the ternary plots, standard geological scales and terminologies were employed. Clay, silt and sand categories comprise the amounts of particles smaller than 7μm, between 7μm and 63μm and between 63μm and 2mm, respectively. The limit of 7μm is based on the demonstration by Konert and Vandenberghe (1997) [146] that this is the value that correlates better, when laser granulometry is used, to the particle size value of the clay fraction (<2μm) of previous sedimentological methods.

Energy dispersive X-ray fluorescence (ED-XRF).

The ED-XRF measurements were carried out using a portable SPECTRO xSORT X-ray fluorescence spectrometer from AMETEK equipped with a silicon drift detector (SDD) and a low power W X-ray tube with an excitation source of 40 kV. Samples were positioned above a 7mm diameter spectrometer aperture and analysed from below for 60 seconds. For data treatment, we used the peak count rates of all detected elements and quantitative data for a selection of elements. The quantitative data were calculated according to a calibration operated with AMETEK X-LabPro software. This calibration was constructed by using 12 certified and local references and allows for the semi-quantification of 7 major and trace elements among the most abundant in ferruginous rocks (normalised concentrations are presented in oxide weight) [147].

X-ray diffraction (XRD).

Ochre samples were analysed by way of X-ray diffraction using a silicon calibrated Panalytical diffractometer with a Bragg-Brentano (theta-theta) configuration. Data acquisition times of 30 minutes allowed an angle range of 8 to 80° with a resolution of 0.02°. The diffraction spectra were analysed using EVA application software coupled with the JCPDS-ICDD database. The mineral species were identified by their respective combinations of two or three characteristic peaks. This first qualitative determination was validated and quantified using TOPAS software. The best fit for each measure was determined using Rietveld’s refinement method [148]. This procedure ensured that the determinations made with EVA did not overlook information masked by overlapping peaks and allowed precise quantification of the relative proportions of each mineral.

Results

Characterisation of ochre samples

Colourimetry.

The colorimetric values of samples are presented in Table 3. Pigment powder produced by Ovahimba women appear colorimetrically homogeneous and display differences in L*a*b* coordinates (ΔE) of less than 2.4. This is well below 5, the value generally considered as the limit above which a significant colour difference can be perceived by humans [149]. The South African Bokkeveld Group red ochre samples display the same range of hue and chroma (similar a* and b* values) as the other red ochre samples. Only lightness (L*) values are dissimilar, with Ovahimba red ochres being darker.

thumbnail
Table 3. L*a*b* colourimetric properties of ethnographic and experimental ochre samples with L* representing lightness (100 white / 0 black), a* indicating variations between red (+) and green (-), and b* representing variations between yellow (+) and blue components (-).

The first sample was arbitrarily selected as the reference sample and the Delta E (ΔE) values calculated accordingly.

https://doi.org/10.1371/journal.pone.0136090.t003

It has been shown that different ochre processing techniques result in differences in pigment powder consistency and colour [121]. In this study, ochre powder produced by direct grinding presents higher chroma (a* and b* values) than powder obtained by grinding with lower and upper grindstones. Noticeable differences in colour resulting from different processing techniques are exhibited by the Napier 15 (ground directly onto a coarse quartzite surface as evidenced by examples derived from MSA and LSA contexts) and 21 (ground conventionally between an upper and lower grindstone such as practiced by the Ovahimba) and the Napier 16 (direct) and 22 (conventional) samples.

Granulometry.

Grain size analysis identifies three broad groups of samples including 1) very fine powders, 2) fine powders and 3) coarse powders (Table 4). The first group has a D10 value close to 1μm and a clay-size proportion approximating or exceeding 60%. This is the case for the samples from Opuwo which were ground directly onto quartzite surfaces. The fine powders display D10 values ranging between 1.5μm and 4μm and a clay content <50%. Most of the samples processed by Ovahimba women fall into this category. The third group is characterised by a D10 close to 5 or higher and a clay-sized content <30%. The samples experimentally processed by way of the technique employed by the Ovahimba belong to this group. Differences in grain size between samples ground by the Ovahimba and experimentally by the authors are significant. Results suggest that, in order to be effective, the Ovahimba technique must be applied to fine-grained iron-rich sources. The application of this technique to process ochre with a high proportion of quartz, such as those from Cape Point, produces coarse-grained pigment powder (Table 4). Every sample shows multimodal grain size distribution (Fig 3) which is classical in grinding processes and especially in these materials comprised of multiple mineral phases of dissimilar hardness.

thumbnail
Fig 3. Representative grain size distributions for each group.

https://doi.org/10.1371/journal.pone.0136090.g003

thumbnail
Table 4. Results for particle size analyses of ochre samples.

https://doi.org/10.1371/journal.pone.0136090.t004

Mineralogical composition.

Analysed samples comprise variable amounts of iron oxides, clays and quartz (Table 5). Goethite and hematite are present in almost all samples, with hematite being the primary constituent. No hematite is detected in the yellow (13 and 19) and grey ochre samples (14 and 20), the former containing only goethite and the latter being the only sample in which iron oxide could not be detected. In comparison with other red samples, Ovahimba red ochre typically contains less quartz and clay minerals (illite and kaolinite) and more iron oxides (Fig 4). Most of the Ovahimba samples are composed of more than 60% iron oxides/oxy-hydroxides. In the samples derived from Cape Point, the concentration in iron oxide/oxy-hydroxide does not exceed 20%.

thumbnail
Fig 4. XRD spectra recorded for ochre samples 5 (Ovinjange fine ground ochre powder), 15 (Napier medium hard red shale chunk) and 19 (Napier soft yellow limonite chunk).

https://doi.org/10.1371/journal.pone.0136090.g004

thumbnail
Table 5. Results of XRD analyses of ethnographic and experimental ochre samples.

https://doi.org/10.1371/journal.pone.0136090.t005

Elemental composition.

Fe and Si are the two major elements detected by ED-XRF (Table 6). The relative contents of Fe (Fe2O3), Si (SiO2) and other minor elements were calculated following calibration and normalization to 100%. CaO and K2O are only present as minor components, except in the sample from Blombos (sample 14) which contain higher proportions of K. TiO2 is also present as a minor constituent. Compared to samples collected from the Bokkeveld shales, the ratio of Fe2O3 to SiO2 is significantly higher in most of the Ovahimba samples. Two samples from Opuwo (11 and 12) present lower Fe contents and lower ratio of Fe2O3 to SiO2 than the other Ovahimba samples. This is consistent with the XRD results. Only the samples from Cape Point (18 and 24) present Fe2O3 to SiO2 ratios as high as those recorded for the Opuwo samples.

thumbnail
Table 6. Results of ED-XRF analyses of ethnographic and experimental ochre samples.

https://doi.org/10.1371/journal.pone.0136090.t006

XRD and ED-XRF data are in general agreement with regards the iron oxide content of the samples. Significant discrepancies between the calculated Fe2O3 and hematite and goethite concentrations concern the samples from Napier (15 and 16) and Cape Point (18 and 24), where the former is twice as high as the latter. This could be explained by a low degree of goethite and hematite crystallinity, leading to an underestimation of their concentration by the Rietveld refinement technique. XRD and ED-XRF are also concordant for K2O and illite/muscovite or microcline content respectively for Cape Point and Ovahimba samples.

SPF values.

The in vitro SPF of the samples range from 1.9 (± 0.1) to 13.1 (± 1.1) (Table 7). Pigment powder produced by Ovahimba women consistently exhibit SPF values higher than those recorded for specimens processed experimentally by the authors. Given their high in vitro SPF and UVAPF values, samples 1, 5, 7, 10, 16 and 18 were subjected to in vivo SPF analysis. The in vivo SPF of the samples, applied to human skin after being mixed with clarified butter, range from 3.5 to 10.7 (Table 8). As indicated by the results obtained for subjects A, B and C, actual and average SPF values are not markedly divergent and experimental samples do display a degree of consistency in terms of their respective SPF values. Sample 18 (Cape Point—hard shale-derived red ochre) provides the highest SPF values, followed by sample 16 (Napier—medium hard shale-derived red ochre) and sample 10 (Opuwo—fine-grained red ochre).

thumbnail
Table 7. In vivo SPF and in vitro UVAPF values and critical wavelengths (CW) of ethnographic and experimental ochre powder samples.

https://doi.org/10.1371/journal.pone.0136090.t007

thumbnail
Table 8. SPF ratings of ochre samples 1, 5, 7, 10, 16 and 18, mixed with clarified butter, and obtained by in vivo assessment on three test subjects (A, B and C).

https://doi.org/10.1371/journal.pone.0136090.t008

SPF values obtained by in vivo assessment for red ochre samples mixed with both clarified butter and with animal fat on test subject A range from 4.5 to 10.0. When mixed with animal fat, SPF values range from 5.4 to 9.8. The type of binder used, whether it is dairy derived clarified butter or subcutaneous animal fat, appears to have no marked effects on the SPF capacity of ochre. A comparison of average SPF values obtained for samples applied in dry form and with binders (clarified butter and animal fat), indicates that dry ochre powder provides a degree of protection against UVR similar to ochre powder mixed with clarified butter and with animal fat. This suggests that it is not necessary to mix ochre powder with a binder to achieve satisfactory SPF values.

Correlations between SPF and chemical and structural characteristics

A positive correlation is observed between in vitro SPF values and the proportion of particles <2μm, hematite and goethite content and Fe2O3 and V2O3 normalised percentages (p = <0.05) (Table 9). A negative correlation is observed between SPF values and the granulometric data D10 to D50, clay mineral content and SiO2. Similar trends are observed for the in vitro UVAPF values.

thumbnail
Table 9. Pearson’s correlation table indicating the relationship between SPF values on the one hand and the granulometric, mineralogical and elemental composition of the pigment samples on the other hand.

Data in bold emphasize prominent positive or negative correlations.

https://doi.org/10.1371/journal.pone.0136090.t009

Figs 5 and 6 explore the relationship between SPF and grain size composition (Fig 5A and 5B), mineral content (Fig 5C and 5D), elemental composition (Fig 5E and 5F) and colour (Fig 6A and 6B). The ternary diagram correlating grain size and SPF (Fig 5A) reveals that the finer the ochre powder, the higher the photoprotective effect (SPF values). This effect becomes particularly intense when the mean size of the particles, at the first 25 centiles (D25) of the pigment grain size distribution, ranges between 1μm and 5μm (Fig 5B). SPF values derived from samples that were processed by way of different techniques are also of interest. SPF values tend to be higher when the grinding process resulted in the incidence of smaller particle sizes. Results also indicate that, when abraded directly on coarse grindstone surfaces, less homogeneous specimens, such as those from the Bokkeveld Group, can generate SPF values comparable to those obtained from homogeneous specimens, such as the Kunene Region (Ombumbuu) samples.

thumbnail
Fig 5. Ternary diagrams and biplots illustrating the relationship between SPF and grain size composition (A, B), mineral content (C, D) and elemental composition (E, F).

https://doi.org/10.1371/journal.pone.0136090.g005

thumbnail
Fig 6. Ternary diagram (A) and biplot (B) illustrating the relationship between SPF values and ochre colour.

https://doi.org/10.1371/journal.pone.0136090.g006

The ternary diagram and biplot (Fig 5C and 5D) illustrate the relationship between SPF and mineralogical composition and confirm the positive correlation between iron-rich minerals and higher SPF values. This is further supported by the results obtained for samples 11 and 12. Despite deriving from the Kunene Region, both exhibit low SPF values. This is largely due to their low hematite and goethite content. These diagrams furthermore show that samples in which iron-rich minerals are represented solely by goethite possess comparable SPF values to those in which hematite or hematite and goethite are present in the same quantities. Elemental analyses reveal a similar pattern. The correlation between the bulk normalised Fe2O3 content and in vitro SPF is shown to be significant, but not linear (Fig 5E and 5F). SPF is positively correlated with ochre samples in which Fe2O3 contents range between 5% and 50%, but is not further influenced once this element constitutes >60% of the sample. Conversely, TiO2 does not appear to influence the in vitro SPF values, although this element is known for its photoprotective capacity [142, 144]. This may be due to its low concentration within the analysed samples (<2%). Finally, the colorimetric properties of the samples are variably correlated with SPF values (Fig 6A and 6B). SPF is negatively correlated with lightness (L*) and positively correlated with increasing a* (red) and b* values. Samples with the brightest red hues, reflecting substantial hematite content, consistently display the highest SPF values. Goethite-rich yellow samples, and samples that have been ground directly onto coarse stone surfaces, appear to provide increased protection from UVR. This also seems to correlate with intensified a* and b* values and greater brightness, which results from this processing technique.

Iron oxide concentration and particle sizes elucidate most, but not all, of the variations observed in SPF and UVAPF values obtained by in vitro experiments. In many instances, other minerals appear to also influence the photoprotective efficacy of powdered ochre. A case in point involves the sample from Napier (13 and 19) which present higher SPF values than samples with similar particle size distribution (15 and 16 and 21 and 22). They differ from the latter examples in terms of higher clay mineral content, lower quartz concentration and the predominance of goethite as opposed to hematite. This suggests that variations in mineral composition do influence the overall photoprotective capacity of ochre. Sample 16 presents a similar particle size distribution and fall within the same range of Fe2O3 content than most Bokkeveld ochre samples processed experimentally using the direct method (sample 15). Nevertheless, this sample exhibits a significantly higher SPF value, which could be explained in terms of the higher concentration of kaolinite present. Whereas directly ground Ovahimba samples (7 to 10) exhibit similar particle size distribution (Fig 5A) and minor variations in mineral composition, they exhibit a surprisingly broad range of SPF values. Other variables, such as iron oxide/oxy-hydroxide crystal size and shape may also play a role in SPF value. It has been shown that crystal shape influences the optical properties of Fe oxides [150], and it is therefore predictable that it also affects the photoprotective properties of Fe-rich ochres.

Discussion and Conclusion

Our experiments demonstrate that red ochre holds significant photoprotective capacities under both in vitro and in vivo conditions. High iron oxide content and smaller grain sizes correlate with greater photoprotective capacity (SPF). Results also confirm that it is not essential to combine powdered ochre with an organic binder to achieve optimal UVA and UVB protection. The use of binders may nevertheless be advantageous in terms of fixing ochre particles to the skin surface and therefore extending the duration of protection. The high photoprotective effect of the single yellow ochre sample was unexpected. Samples from the same location (Napier) present variable photoprotective (SPF) values which does not fit the general pattern we observe for the other samples [100, 142, 144]. Because this represents a minor proportion of the analysed instances, our discussion will focus on the broader results which suggest that the finer the ochre powder and the higher its hematite content the higher the photoprotective effect.

Although the influence of regional availability on ochre procurement strategies has been highlighted by Watts (2002) [114] and Dayet et al. (2013) [123], it is clear that MSA and LSA people deliberately selected ochre depending on its properties, and not only on its ease of acquisition. Colourimetric properties appear to have also been significant in terms of LSA and MSA ochre selection strategies. Several MSA specimens exhibit what may represent traces of streak-testing [133, 151], suggesting that ochre was consciously selected according to specific criteria. At Klasies River in South Africa, dark red ochre with unique colourimetric and mineralogical properties and therefore different from samples that were ground to extract powder, was selected for engraving purposes at 100 ka [135]. At 90 ka at Qafzeh Cave in Israel, local ochre sources were ignored in favour of more remote, but finer-grained, red sedimentary shales some 60 km from the site [152]. At 65 ka at Diepkloof Rock Shelter, and even though red ochre occurs inside the shelter, ochre from sources more than 20 km away from the site were preferentially processed [123].

Because the colour of ochre is generally indicative of elemental content and mineral composition [150], colour may have been a key criterion for ochre selection, depending on the intended purpose or applications in mind. At 100 ka, the Mousterian and Middle Palaeolithic inhabitants of the Es-Skhul (Israel) collected yellow goethite-rich ochre from 80 km away. These specimens were subsequently heated to induce colour transformations to red [153]. Based on our experimental results, one should, in addition to yellow examples, therefore expect to find substantial amounts of dark red shale-derived ochres in southern African MSA contexts. Ochre with a strong red hue does in fact appear to predominate in MSA contexts from 130 ka [114, 115, 133]. But colour-based selection does not necessarily imply a preference based on aesthetic or symbolic considerations. In addition to colour, ochre may have been selected for other physical and chemical characteristics such as grain size, processability and surface covering capacity. Cultural reasons, such as symbolic connotations to particular geological resources, may have also played a role. The widespread preference for and use of Ombumbuu red ochre, amongst the Ovahimba, suggest that colour-independent properties may also play an important role in the selection and preference of red ochre from this source.

Although ochre is attested in southern and eastern African and Levantian sites dated to at least 300 ka [151158], its use only becomes ubiquitous during the last interglacial. During the early part of MIS 5, from 127 ka to 116 ka when temperature, rainfall and vegetation patterns are likely to have resembled those of today, human populations became archaeologically increasingly visible. The essentially only detectable behavioral signature of this transformation, the widespread exploitation of ochre, designates an increasingly complex degree of cultural complexity emerging in southern Africa at the end of MIS 5. We are cautious to argue in favour of a direct linear relationship between ochre frequency and geomagnetic excursion events and acknowledge that further research is required to substantiate this likelihood. The timing of the Blake and post-Blake geomagnetic excursions, estimated to have spanned the period between 138 ka to 100 ka [28, 29], and, more specifically, to have occurred at 120 ka to 115 ka and 100 ka to 95 ka, respectively [159], appears to both precede and coincide with the increasing prevalence of red ochre exploitation at Blombos Cave and at Klasies River Cave 1 after 120 ka [117, 133, 135]. Increased cosmogenic nuclide deposition rates, indicated by long lasting relative paleointensity records associated with full geomagnetic polarity reversals, are confirmed for both the Blake and the post-Blake events [159]. There is currently limited evidence for increasing or declining ochre exploitation following the Laschamp excursion after 40 ka. This is because we have very little information concerning the amounts and types of ochre exploited in southern Africa between 40 ka and 20 ka [160]. It is therefore difficult to assess the relationship, if any, between the suite of adaptations recorded before 60 ka and those that emerged with the LSA. Although a linear relationship between ochre frequency and insolation cannot be demonstrated at this point, the initial large-scale exploitation of red ochre does seem to coincide with the Blake excursion event after 120 ka. However, and in order to determine the exact time period in which the use of ochre was uniquely due to a cultural adaptation against insolation, one should probably consider much earlier periods, perhaps from 300 ka or more. The precise motivation for this increasing focus on red ochre is unclear. It might represent both an expansion and an exaptation of various former functions fulfilled by this material. Significantly, recent evidence derived from Sibudu Cave [161] indicate that a mixture of red ochre and casein from milk, possibly obtained from hunted lactating wild bovids, was produced by the inhabitants at 49 ka. This liquid mixture, consisting of powdered red ochre and milk, was in all probability used as a paint-like medium that could have been applied to human skin.

Our results provide a viable foundation for the emergence of ochre use, suggesting that the topical application of red ochre may have been a key innovation that functioned primarily to restrict the adverse effects of UVR exposure [100]. The skin colour of southern African MSA Homo sapiens is likely to have resembled the olive-brown ‘Capoid’ skin type IV (Von Luschan 16 to 21) typical of some indigenous southern African groups [76, 139]. This skin type is characteristic of more than 90% of southern San individuals, has a natural SPF of about 3.5 that of Phototype II Caucasian skin, a high tanning capacity and lower susceptibility to UVR-induced damage. In terms of avoiding vitamin D deficiencies and negating the detrimental effects of UVR, and given the correlation between adequate UVR exposure and human health, an olive-brown skin tone represents the optimal skin type for inhabitants of southern and eastern sub-Saharan Africa. It therefore stands to reason that, given the natural SPF of 3.5 for skin type IV, ochres exhibiting SPF and UVAPF values approximating or exceeding 10, represented here by the dark red ochre specimens 10, 16 and 18, are likely to have provided sufficient protection from UVR without initiating vitamin D deficiencies [139]. The topical application of red ochre may consequently have facilitated the colonisation of areas unfavourable to the constitutive skin colour of foreign immigrants, thereby increasing the ability of newcomers to compete with better adapted local populations. Subsequently, such functional innovations may have reinforced existing cultural mechanisms related to individual, intra- and inter-group identification, which in turn resulted in the emergence of ever more complex and diversified symbolic material cultures. At an individual level, ochre could have fulfilled a central role in processes of information exchange by visually enhancing the bodies, faces, hair and ornaments and attire of individuals. At the group level, ochre may have served to stimulate the formation of social and economic relationships, reinforce group solidarity and ensure cooperation or discern different statuses, achievements or affiliations amongst group members. Few Darwinian models have tried to explain the emergence of ochre use in human history. The Female Cosmetic Coalitions (FCC) model [162, 163] proposes an interesting explanatory framework for the increased exploitation of red ochre based on evolutionary constrains and the symbolic embodiment of relationships between men and women. The sunscreen model, proposed here, is not in contradiction with the FCC model in the sense that it may well explain the emergence of the phenomenon at a stage when the various symbolic functions of ochre were not yet in place. While it is difficult to establish a precise scenario for the emergence of such an innovation, one could envision a situation in which the habitual use of red ochre as a sunscreen may have originally arisen locally, perhaps in response to changing UVR exposure rates produced by orbital climatic cycles, and subsequently provided an adaptive advantage for migrating populations. It is likely that the use of ochre as a sunscreen emerged repeatedly and that this innovation periodically vanished and reoccurred, a hypothesis consistent with the chronological and geographic discontinuity in ochre occurrences during this time span. We do not know, at present, whether such applications may have facilitated the latitudinal displacement of populations and whether it coincided with symbolic connotations. It can, however, be argued that the increase in ochre use corresponding to the beginning of the last interglacial in all probability reflects a gradual, albeit significant, change of magnitude in the role of ochre in human culture, which has not come to an end since it emerged [139].

Combined with experimental replication and multi-proxy analyses, ethnographic examples of red ochre exploitation enabled the identification of noteworthy ancestral functions, while also revealing the extraordinary degree of social complexity that it has contributed to since its initial exploitation. Although the transition in archaeology from data to social interpretation is challenging, the importance of ochre exploitation during the MSA can only be expressed in terms of complex multi-functional explanations. This is so because the advantages offered by the habitual exploitation of ochre during the evolution of Homo sapiens were numerous. Although exclusively functional uses of ochre are rare, it can be argued that they may have been more common in the past and that the use of ochre for symbolic purposes in ethnographic contexts does not imply similar use by early MSA humans. Ochre could also have functioned symbolically before the emergence of anatomically modern Homo sapiens after 195 ka [164]), but it may in fact have had several functional uses early on which subsequently evolved to comprise both symbolic significance and functional value. The hypothesized and confirmed functional applications of red ochre should not therefore be viewed as in contradiction with a symbolic use of ochre, but instead as indicative of the possible exploitation and extensive application of a raw material with remarkable colourimetric, structural and chemical properties during the MSA.

Acknowledgments

We thank the Ovahimba participants, in particular Nakara, Mokadera, Bendokarampi and Tharhirwa (Otjongoro, Ovinjange and Oukongo villages) for sharing their knowledge with us. We thank the Namibian Ministry of Culture, the director of the National Museum of Namibia (Ester Moombolah-Goagoses), Eugene Marais (National Museum of Namibia) and Luc Fabre (French Embassy in Namibia) for supporting our research. Anzel Veldman (National Museum of Namibia) and Sennobia Katjiuongua (National Heritage Council of Namibia) are thanked for their support in the field. Renata Garcia Moreno Mazel (CNRS UMR PACEA, University of Bordeaux) is thanked for preliminary analyses.

Author Contributions

Conceived and designed the experiments: RFR FDE LD AQ BS. Performed the experiments: RFR ML BS LD AQ. Analyzed the data: AQ LD FDE RFR ML BS. Contributed reagents/materials/analysis tools: AQ LD BS ML. Wrote the paper: RFR FDE AQ LD.

References

  1. 1. Blome MW, Cohen AS, Tryon CA, Brooks AS, Russell J (2012) The environmental context for the origins of modern human diversity: A synthesis of regional variability in African climate 150,000–30,000 years ago. J Hum Evol 62 (5): 563–592. pmid:22513381
  2. 2. Nyman T, Linder HP, Pen C, Malm T, Wahlberg N (2012) Climate-driven diversity dynamics in plants and plant-feeding Insects. Ecol Lett 14: 1–10.
  3. 3. Condamine FL, Rolland J, Morlon H (2013) Macroevolutionary perspectives to environmental change. Ecol Lett 16: 72–85.
  4. 4. Ziegler M, Simon MH, Hall IR, Barker S, Stringer C, Zahn R (2013) Development of Middle Stone Age innovation linked to rapid climate change. Nat Commun 4 (1905):
  5. 5. Paul ND, Gwynn-Jones D (2003) Ecological roles of solar UV radiation: Towards an integrated approach. Trends Ecol Evol 18 (1): 48–55.
  6. 6. Lautenschlager S, Wulf HC, Pittelkow MR (2007) Photoprotection. Lancet 370: 528–537. pmid:17693182
  7. 7. Gasparro FP (2000) Sunscreens, skin photobiology and skin cancer: The need for UVA protection and evaluation of efficacy. Environ Health Perspect 108 (1): 71–78.
  8. 8. Moan J, Porojnicu AC, Dahlback A, Setlow RB (2007) Addressing the health benefits and risks, involving vitamin D or skin cancer, of increased sun exposure. Proc Natl Acad Sci U S A 105 (2): 668–673.
  9. 9. Reichrath J (2007) Vitamin D and the skin: An ancient friend, revisited. Exp Dermatol 16: 618–625. pmid:17576242
  10. 10. Yuen AWC, Jablonski NG (2010) Vitamin D: In the evolution of human skin colour. Med Hypotheses 74: 39–44. pmid:19717244
  11. 11. Reichrath J (2006) The challenge resulting from positive and negative effects of sunlight: How much solar UV exposure is appropriate to balance between risks of vitamin D deficiency and skin cancer? Prog Biophys Mol Biol 92: 9–16. pmid:16603232
  12. 12. Juzeniene A, Setlow R, Porojnicu A, Hykkerud Steindal A, Moan J (2009) Development of different human skin colours: A review highlighting photobiological and photobiophysical aspects. J Photochem Photobiol B B: Biology 96: 93–100. pmid:19481954
  13. 13. Shoenfeld N, Amital H, Shoenfield Y (2009) The effect of melanism and vitamin D synthesis on the incidence of autoimmune disease. Nat Clin Pract Rheumatol 5 (2): 99–105. pmid:19182816
  14. 14. Wright CY, Norval M, Summers B, Davids L, Coetzee G, Oriowo MO (2012) The impact of solar ultraviolet radiation on human health in sub-Saharan Africa. South Afr J Sci 108 (11–12): 1–6.
  15. 15. Biniek K, Levi K, Dauskardt RH (2012) Solar UV radiation reduces the barrier function of human skin. Proc Natl Acad Sci U S A 109 (42): 17111–17116. pmid:23027968
  16. 16. Krishnan KJ, Harbottle A, Birch-Machin MA (2004) The use of a 3895 bp mitochondrial DNA deletion as a marker for sunlight exposure in human skin. J Investig Dermatol 123: 1020–1024. pmid:15610508
  17. 17. Kadekaro AL, Kavanagh R, Kanto H, Terzieva S, Hauser J, Kobayashi N, et al. (2005) α-Melanocortin and endothelin-1 activate antiapoptotic pathways and reduce DNA damage in human melanocytes. Cancer Res 65 (10): 4292–4299. pmid:15899821
  18. 18. Gendron SP, Bastien N, Mallet JD, Rochette PJ (2013) The 3895-bp mitochondrial DNA deletion in the human eye: a potential involvement in corneal ageing and macular degeneration. Mutagen 28 (2): 197–204.
  19. 19. Chen TC, Lu Z, Holick MF (2010) Photobiology of vitamin D. Holick , Springer NY.
  20. 20. Ullrich SE, Kripke ML, Ananthaswamy HN (2002) Mechanisms underlying UV-induced immune suppression: Implications for sunscreen design. Exp Dermatol 11 (1): 1–4.
  21. 21. Knudsen MF, Riisager P, Donadini F, Snowball I, Muscheler R, Korhonen K, et al. (2008) Variations in the geomagnetic dipole moment during the Holocene and the past 50 kyr. Earth Planet Sci Lett 272 (1–2): 319–329.
  22. 22. Faith JT, Surovell TA (2009) Synchronous extinction of North America’s Pleistocene mammals. Proc Natl Acad Sci U S A 106 (49): 20641–20645. pmid:19934040
  23. 23. LaViolette PA (2011) Evidence for a solar flare cause of the Pleistocene mass extinction Radiocarbon 53 (2): 303–323.
  24. 24. Valet J-P, Valladas H (2010) The Laschamp-Mono lake geomagnetic events and the extinction of Neanderthals: A causal link or a coincidence? Quat Sci Rev 29: 3887–3893.
  25. 25. Grove M (2012) Amplitudes of orbitally induced climatic cycles and patterns of hominin speciation. J Archaeol Sci 39: 3085–3094.
  26. 26. Laj C, Channell JET (2008) Geomagnetic excursions. Kono, Elsevier, Amsterdam.
  27. 27. Thouveny N, Carcaillet J, Moreno E, Leduc G, Nérini D (2004) Geomagnetic moment variation and palaeomagnetic excursions since 400 kyr BP: A stacked record from sedimentary sequences of the Portuguese margin. Earth Planet Sci Lett 219: 377–396.
  28. 28. Bourne M, Niocaill CM, Thomas AL, Knudsen MF, Henderson GM (2012) Rapid directional changes associated with a 6.5 kyr-long Blake geomagnetic excursion at the Blake-Bahama Outer Ridge. Earth Planet Sci Lett 333–334: 21–34.
  29. 29. Osete M, Martın-Chivelet J, Rossi C, Edwards RL, Egli R, Muñoz-García MB, et al. (2012) The Blake geomagnetic excursion recorded in a radiometrically dated speleothem. Earth Planet Sci Lett 353–354: 173–181.
  30. 30. Guillou H, Singer BS, Laj C, Kissel C, Scaillet S, Jicha BR (2004) On the age of the Laschamp geomagnetic excursion. Earth Planet Sci Lett 227: 331–343.
  31. 31. Plenier G, Valet J-P, Guérin G, Lefèvre JC, Carter-Stiglitz B (2007) Origin and age of the directions recorded during the Laschamp event in the Chaîne des Puys (France). Earth Planet Sci Lett 259: 414–431.
  32. 32. Finlayson C, Carrión JS (2007) Rapid ecological turnover and its impact on Neanderthal and other human populations. Trends Ecol Evol 22 (4): 213–222. pmid:17300854
  33. 33. Banks WE, d'Errico F, Peterson AT, Kageyama M, Sima A, Sánchez-Goñi M-F (2008) Neanderthal extinction by competitive exclusion. PLOS ONE
  34. 34. Gilligan I (2008) Neanderthal extinction and modern human behaviour: The role of climate change and clothing. World Archaeology 39 (4): 499–514.
  35. 35. Dennell RW, Martinón-Torres M, Bermúdez de Castro JM (2011) Hominin variability, climatic instability and population demography in Middle Pleistocene Europe Quaternary Science Reviews 30: 1511–1524.
  36. 36. Mellars P, French JC (2011) Tenfold population increase in Western Europe at the Neanderthal-to-modern human transition. Science 333: 623–627. pmid:21798948
  37. 37. Müller UC, Pross J, Tzedakis PC, Gamble C, Kotthoff U, Schmiedl G, et al. (2011) The role of climate in the spread of modern humans into Europe. Quat Sci Rev 30: 273–279.
  38. 38. Sørensen B (2011) Demography and the extinction of European Neanderthals. J Anthropol Archaeol 30: 17–29.
  39. 39. Discamps E (2013) Ungulate biomass fluctuations endured by Middle and Early Upper Paleolithic societies (SW France, MIS 5–3): The contributions of modern analogs and cave hyena paleodemography. Quat Int 337: 64–79.
  40. 40. Lalueza-Fox C, Römpler H, Caramelli D, Stäubert C, Catalano G, Hughes D, et al. (2007) A melanocortin 1 receptor allele suggests varying pigmentation among Neanderthals. Science 318: 1453–1455. pmid:17962522
  41. 41. Robinson KC, Fisher DE (2009) Specification and loss of melanocyte stem cells. Semin Cell Dev Biol 20: 111–116. pmid:19124082
  42. 42. Mitra D, Luo X, Morgan A, Wang J, Hoang MP, Lo J, et al. (2012) An ultraviolet-radiation-independent pathway to melanoma carcinogenesis in the red hair/fair skin background. Nature 15 (7424):449–453.
  43. 43. Cerqueira CC, Pãixa-Côrtes VR, Zambra FM, Salzano FM, Hünemeier T, Bortolini C (2012) Predicting Homo pigmentation phenotype through genomic data: From Neanderthal to James Watson. Am J Hum Biol 24: 705–709. pmid:22411106
  44. 44. Meyer M, Kircher M, Gansauge M, Li H, Racimo F, Mallick S, et al. (2012) A high coverage genome sequence from an archaic Denisovan individual. Science 338 (6104): 222–226. pmid:22936568
  45. 45. Ménabréaz L, Thouveny N, Bourlès DL, Deschamps P, Hamelin B, Demory F (2011) The Laschamp geomagnetic dipole low expressed as a cosmogenic 10Be atmospheric overproduction at ~ 41 ka. Earth Planet Sci Lett 31 (4–5): 305–317.
  46. 46. Nowaczyk NR, Arz HW, Frank U, Kind J, Plessen B (2012) Dynamics of the Laschamp geomagnetic excursion from Black Sea sediments. Earth Planet Sci Lett 351–352: 54–69.
  47. 47. Hodgson DA, Sabbe K, Vyverman W, Squier AH, Verleyen E, Leavitt PR, et al. (2005) Late Pleistocene record of elevated UV radiation in an Antarctic lake. Earth Planet Sci Lett 236: 765–772.
  48. 48. Winkler H, Sinnhuber M, Notholt J, Kallenrode M, Steinhilber F, Voght J, et al. (2008) Modelling impacts of geomagnetic field variations on middle atmospheric ozone responses to solar proton events on long timescales. J Geophys Res 113: 1–11.
  49. 49. World Health Organisation 2002. Global solar UV index: A practical guide. WHO/SDE/OEH/02.2.
  50. 50. Gilaberte Y, González S (2010) Update on photoprotection. Actas Dermosifiliogr 101 (8): 659–672. pmid:20965009
  51. 51. Light JE, Reed DL (2009) Multigene analysis of phylogenetic relationships and divergence times of primate sucking lice (Phthiraptera: Anoplura). Mol Phylogenet Evol 50: 376–390. pmid:19027083
  52. 52. Toups MA, Kitchen A, Light JE, Reed DL (2010) Origin of clothing lice indicates early clothing use by anatomically modern humans in Africa. Mol Biol Evol 28 (1): 29–32. pmid:20823373
  53. 53. Brown KS, Marean CW, Herries AIR, Jacobs Z, Tribolo C, Braun D, et al. (2009) Fire as an engineering tool of early modern humans. Science 325: 859–862. pmid:19679810
  54. 54. Mourre V, Villa P, Henshilwood CS (2010) Early use of pressure flaking on lithic artefacts at Blombos Cave, South Africa. Science 330: 659–662. pmid:21030655
  55. 55. Dusseldorp GL (2011) Tracking the influence of technological change on Middle Stone Age hunting strategies in South Africa. Quat Int 270: 70.79.
  56. 56. Lombard M (2011) Quartz-tipped arrows older than 60 ka: Further use-trace evidence from Sibudu, KwaZulu-Natal, South Africa. J Archaeol Sci 38: 1918–1930.
  57. 57. Jablonski NG (2004) The evolution of human skin and skin colour. Ann Rev Anthropol 33: 585–623.
  58. 58. Zihlman AL, Cohn BA (1988) The adaptive response of human skin to the savannah. J Hum Evol 33: 397–409.
  59. 59. Gilchrest BA (2008) Sun exposure and vitamin D sufficiency. American Journal of Clinical Nutrition 88: 570–577.
  60. 60. Jablonski NG, Chaplin G (2000) The evolution of human skin coloration. J Hum Evol 39 (1): 57–106. pmid:10896812
  61. 61. MacKintosh JA (2001) The antimicrobial properties of melanocytes, melanosomes and melanin and the evolution of black skin. J Theor Biol 211: 101–113. pmid:11419954
  62. 62. Elias M, Williams ML (2013) Re-appraisal of current theories for the development and loss of epidermal pigmentation in hominins and modern humans. J Hum Evol 64 (6): 687–692. pmid:23478159
  63. 63. Jablonski NG (2012) The evolution of human skin colouration and its relevance to health in the modern world. J R Coll Physicians Edinb 42 (1): 58–63. pmid:22441067
  64. 64. Tam TT, Juzeniene A, Steindal AH, Iani V, Moan J (2009) Photodegradation of 5-methyltetrahydrofolate in the presence of uroporphyrin. Photochem Photobiol 94: 201–204.
  65. 65. Jablonski NG, Chaplin G (2010) Human skin pigmentation as an adaptation to UV radiation. Proc Natl Acad Sci U S A 107 (2): 8962–8968.
  66. 66. Sturm RA (2009) Molecular genetics of human pigmentation diversity. Hum Mol Genet
  67. 67. Parra EJ (2007) Human pigmentation variation: Evolution, genetic basis and implications for public health. Yearb Phys Anthropol 50: 85–105.
  68. 68. Wagner JK, Parra EJ, Norton HL, Jovel C, Schriver MD (2002) Skin responses to ultraviolet radiation: Effects of constitutive pigmentation, sex and ancestry. Pigment Cell Melanoma Res 15: 385–390.
  69. 69. Tadokoro T, Kobayashi N, Zmudzka BZ, Ito S, Wakamatsu K, Yamaguchi Y, et al. (2003) UV-induced DNA damage and melanin content in human skin differing in racial/ethnic origin. FASEB J 17 (9): 1177–1179. pmid:12692083
  70. 70. Fluhr JW, Darlenski R, Berardesca E (2008) Ethnic groups and sensitive skin: Two examples of special populations in dermatology. Drug Discov Today Dis Mech 5 (2): 249–263.
  71. 71. Tishkoff SA, Kidd KK (2004) Implications of biogeography of human populations for race and medicine. Nat Genet 36 (11): 21–27.
  72. 72. Fitzpatrick TB (1988) The validity and practicality of sun-reactive skin types I through VI. Arch Dermatol 124 (6): 869–871. pmid:3377516
  73. 73. Cripps DJ (1981) Natural and artificial photoprotection. J Investig Dermatol 76: 154–157.
  74. 74. Herrling T, Jung K, Fuchs J (2007) The role of melanin as protector against free radicals in skin and its role as free radical indicator in hair. Spectrochim Acta A Mol Biomol Spectrosc 69 (5): 1429–1435. pmid:17988942
  75. 75. Young AR (2006) Acute effects of UVR on human eyes and skin. Prog Biophys Mol Biol 92: 80–85. pmid:16600340
  76. 76. Swiatoniowski AK, Quillen EE, Shriver MD, Jablonski NG (2013) Technical note: Comparing von Luschan skin color tiles and modern spectrophotometry for measuring human skin pigmentation. Am J Phys Anthropol 151 (2): 325–330. pmid:23633083
  77. 77. Brenner M, Hearing BJ (2008) The protective role of melanin against UV damage in human skin. Photochem Photobiol 84 (3): 539–549. pmid:18435612
  78. 78. Lao O, de Gruijter JM, van Duijn K, Navarro A, Kayser M (2006) Signatures of positive selection in genes associated with human skin pigmentation as revealed from analyses of single nucleotide polymorphisms. Ann Hum Genet 71: 354–369.
  79. 79. Barbujani G, Colonna V (2010) Human genome diversity: Frequently asked questions. Trends Genet 26: 285–295. pmid:20471132
  80. 80. Beleza S, Johnson NA, Candille SI, Absher DM, Lopes J, Campos J, et al. (2012) Genetic architecture of skin and eye color in an African European admixed population. PLOS Genet 9 (3): 1–15.
  81. 81. Pickrell JK, Patterson N, Loh P, Lipson M, Berger B, Stoneking M, et al. (2013) Ancient west Eurasian ancestry in southern and eastern Africa. Proc Natl Acad Sci U S A 11 (7): 2632–2637.
  82. 82. Burchell WJ (1822) Travels in the interior of southern Africa. Struik, Cape Town.
  83. 83. Methuen HH (1846) Life in the wilderness. Richard Bentley, London.
  84. 84. Stow GW (1905) The native races of South Africa. Struik, Cape Town.
  85. 85. Bleek WHI, Lloyd L (1911) Specimens of Bushman folklore. George Allen, London.
  86. 86. Marshall L (1961) Sharing, talking, and giving: Relief of social tensions among the! Kung Bushmen. Africa 31 (3): 231–249.
  87. 87. Marshall L (1976) The! Kung of Nyae Nyae. Harvard University Press, Cambridge.
  88. 88. Silberbauer GB (1965) Report to the Government of Bechuanaland on the Bushman survey. Bechuanaland Government Press, Gaberone.
  89. 89. Viegas Guerreiro M (1968) Bochimanes! khu de Angola; estudo ethnografico. Instituto de Investigacao Cientifica de Angola. Hunta de Investigacoes do Ultramar, Lisbon.
  90. 90. Thunberg CP (1775) Travels at the Cape of Good Hope. Van Riebeek Society, Cape Town.
  91. 91. Campbell J (1816) Travels in South Africa: Undertaken at the request of the London Missionary Society. Flagg and Gould, London.
  92. 92. Gordon JR (1786) Cape Travels. Brenthurst Press, Johannesburg.
  93. 93. Strecker IA (1976) Traditional life and prospects for socio-economic development in the Hamar administrative district of southern Gamu Gofa. A report to the relief and rehabilitation commission of the provisional military government of Ethiopia, Addis Ababa.
  94. 94. Lydell J, Strecker I (1979) The Hamar of southern Ethiopia. Renner Verlag, Hohenschäftlarn.
  95. 95. Galton F (1853) The narrative of an explorer in tropical South Africa. John Murray, London.
  96. 96. Johnston HH (1884) On the races of the Congo and the Portuguese colonies in West Africa. J Roy Anthropol Inst Great Brit Ireland 13 (1884): 461–479.
  97. 97. Tönjes H (1911) Ovamboland: Country people mission with particular reference to the largest tribe, the Kwanyama. Namibia Scientific Society, Windhoek.
  98. 98. Kuntz J (1912) Ovatschimba im nordlichen Kaokofeld (Deutsch-Sadwestafrika). Petermanns Geogr Mitt 25: 206.
  99. 99. Vedder H (1928) The Herero. Hahn, Vedder, Fourie, Cape Times, Cape Town.
  100. 100. Rifkin RF, d’Errico F, Summers B, Dayet L (2015) Assessing the photoprotective effects of red ochre on human skin by in vitro laboratory experiments. South Afr J Sci 111 (3/4) 1–8.
  101. 101. Nicoll K (2010) Geomorphic development and Middle Stone Age archaeology of the Lower Cunene River, Namibia-Angola Border. Quat Sci Rev 29: 1419–1431.
  102. 102. Alexander JE (1838) An expedition of discovery into the interior of southern Africa. Carey and Hart, Philadelphia, PA.
  103. 103. Gear S (2002) Wishing us away: Challenges facing ex-combatants in the 'new' South Africa. Viol Transit 8: 1–33.
  104. 104. Gewald J-B (2010) Remote but in contact with history and the world. Proc Natl Acad Sci U S A 107 (18): 75.
  105. 105. Soressi M, d’Errico F (2007) Pigments, gravures, parures: Les comportements symboliques controversés des Néandertaliens. Vandermeersch, Maureille, Éditions du CTHS, Paris.
  106. 106. Zilhão J, Angelucci DE, Badal-García E, d’Errico F, Daniel F, Dayet L, et al. (2010) Symbolic use of marine shells and mineral pigments by Iberian Neanderthals. Proc Natl Acad Sci U S A 107 (3): 1023–1028. pmid:20080653
  107. 107. Roebroeks W, Sier MJ, Kellberg Nielsen T, De Loecker D, Maria Parés J, Arps CES, et al. (2012) Use of red ochre by early Neanderthals. Proc Am Philos Soc 109 (6): 1889–1894.
  108. 108. d’Errico F, Henshilwood CS, Vanhaeren M, van Niekerk K (2005) Nassarius kraussianus shell beads from Blombos Cave: Evidence for symbolic behaviour in the Middle Stone Age. J Hum Evol 48: 3–24. pmid:15656934
  109. 109. d’Errico F, Vanhaeren M, Henshilwood CS, Lawson G, Maureille B, Gambier D, et al. (2009) From the origin of language to the diversification of languages: What can archaeology and palaeoanthropology say? d’Errico Hombert, John Benjamins, Amsterdam.
  110. 110. Vanhaeren M, d’Errico F, Stringer C, James SL, Todd JA, Mienis HK (2006) Middle Palaeolithic shell beads in Israel and Algeria. Science 312: 1785–1788. pmid:16794076
  111. 111. Vanhaeren M, d’Errico F, van Niekerk KL, Henshilwood CS, Erasmus RM (2013) Thinking strings: Additional evidence for personal ornament use in the Middle Stone Age at Blombos Cave, South Africa. J Hum Evol 64 (6): 500–517. pmid:23498114
  112. 112. Bouzouggar A, Barton N, Vanhaeren M, d’Errico F, Collcutt S, Higham T, et al. (2007) 82 000-year-old shell beads from North Africa and implications for the origins of modern human behaviour. Proc Natl Acad Sci U S A 104: 9964–9969. pmid:17548808
  113. 113. Bar-Yosef Mayer DE, Vandermeersch B, Bar-Yosef O (2009) Shells and ochre in Middle Palaeolithic Qafzeh Cave, Israel: Indications for modern behaviour. J Hum Evol 56: 307–314. pmid:19285591
  114. 114. Watts I (2002) Ochre in the Middle Stone Age of southern Africa: Ritualized display or hide preservative? S Afr Archaeol Bull 57: 1–14.
  115. 115. Watts I (2009) Red ochre, body-painting, and language: Interpreting the Blombos ochre. Botha , Knight , Oxford University Press, Oxford.
  116. 116. Henshilwood CS, d’Errico F, Vanhaeren M, van Niekerk K, Jacobs Z (2004) Middle Stone Age shell beads from South Africa. Science 304: 404. pmid:15087540
  117. 117. Henshilwood CS, d’Errico F, van Niekerk KL, Coquinot Y, Jacobs Z, Lauritzen S-E, et al. (2011) A 100,000-year-old ochre-processing workshop at Blombos Cave, South Africa. Science 334 (6053): 219–222. pmid:21998386
  118. 118. Wadley L (2005) Ochre crayons or waste products? Replications compared with MSA crayons from Sibudu Cave, South Africa. Before Farming 2005 (3): 1–12.
  119. 119. Wadley L (2006) Revisiting cultural modernity and the role of ochre in the Middle Stone Age. Soodyall , Jonathan Ball, Johannesburg.
  120. 120. Rifkin RF (2011) Assessing the efficacy of red ochre as a prehistoric hide-tanning ingredient. J Afr Archaeol 9 (2): 1–28.
  121. 121. Rifkin RF (2012) Processing ochre in the Middle Stone Age: Testing the inference of prehistoric behaviours from actualistically derived experimental data. J Anthropol Archaeol 31: 174–195.
  122. 122. Hodgskiss T (2012) An investigation into the properties of the ochre from Sibudu, KwaZulu-Natal, South Africa. South Afr Humanit 24: 99–120.
  123. 123. Dayet L, Texier P-J, Daniel F, Porraz G (2013) Ochre resources from the Middle Stone Age sequence of Diepkloof Rock Shelter, Western Cape, South Africa. J Archaeol Sci 40 (9): 3492–3505.
  124. 124. Dynesius M, Jansson R (2000) Evolutionary consequences of changes in species’ geographical distributions driven by Milankovitch climate oscillations. Proc Natl Acad Sci U S A 97 (16): 9115–9120. pmid:10922067
  125. 125. Hancock AM, Witonsky DB, Gordon AS, Eshel G, Pritchard JK, Coop G, et al. (2008) Adaptations to climate in candidate genes for common metabolic disorders. PLOS Genet 4 (2): 1–13.
  126. 126. Hancock AM, Alkorta-Aranburu G, Witonsky DB, Di Rienzo A (2010) Adaptations to new environments in humans: The role of subtle allele frequency shifts. Philos Trans R Soc Lond B Biol Sci 365 (1552): 2459–2468. pmid:20643735
  127. 127. Ezard TH, Aze T, Pearson PN, Purvis A (2011) Interplay between changing climate and species’ ecology drives macroevolutionary dynamics. Science 332 (6027): 349–351. pmid:21493859
  128. 128. Kanazawa S (2008) Temperature and evolutionary novelty as forces behind the evolution of general intelligence. Intell 36 (2): 99–108.
  129. 129. Kanazawa S (2011) Intelligence and physical attractiveness. Intell 39 (1): 7–14.
  130. 130. Kummu M, Varis O (2011) The world by latitudes: A global analysis of human population, development level and environment across the north-south axis over the past half century. Appl Geogr 31: 495–507.
  131. 131. Rindermann H, Woodley MA, Stratford J (2012) Haplogroups as evolutionary markers of cognitive ability. Intell 40 (4): 362–375.
  132. 132. Mackay A, Welz A (2008) Engraved ochre from a Middle Stone Age context at Klein Kliphuis in the Western Cape of South Africa. J Archaeol Sci 35: 1521–1532.
  133. 133. Henshilwood CS, d’Errico F, Watts I (2009) Engraved ochres from Middle Stone Age levels at Blombos Cave, South Africa. J Hum Evol 57: 27–47. pmid:19487016
  134. 134. Watts I (2010) The pigments from Pinnacle Point Cave 13B, Western Cape, South Africa. J Hum Evol 59: 392–411. pmid:20934093
  135. 135. d’Errico F, Garcia Moreno R, Rifkin RF (2012) Technological, elemental and colorimetric analysis of an engraved ochre fragment from the Middle Stone Age levels of Klasies River Cave 1, South Africa. J Archaeol Sci 39: 942–952.
  136. 136. Laland KN, Odling-Smee J, Myles S (2010) How culture shaped the human genome: Bringing genetics and the human sciences together. Nature 11: 137–148.
  137. 137. Odling-Smee FJ, Laland KN, Feldman MW (2003) Niche construction: The neglected process in evolution. Princeton University Press, Princeton, NJ.
  138. 138. Nettle D (2009) Ecological influences on human behavioural diversity: A review of recent findings. Trends Ecol Evol 24 (11): 618–624. pmid:19683831
  139. 139. Rifkin RF, d’Errico F, Summers B, Dayet L (2015) Assessing the photoprotective effects of red ochre on human skin by in vitro laboratory experiments. S Afr J Sci 111 (3/4) 1–8.
  140. 140. Vorster CJ (2002) Simplified geology and iron deposits: South Africa, Lesotho and Swaziland, Council for Geoscience, Pretoria.
  141. 141. Bernerd F, Vioux C, LeJeune F, Asselineau D (2003) The sun protection factor (SPF) inadequately defines broad spectrum photoprotection: Demonstration using skin reconstructed in vitro exposed to UVA, UVB or UV-solar simulated radiation. Eur J Dermatol 13: 242–249. pmid:12804982
  142. 142. Hoang-Minh T, Le TL, Kasbohm J, Gieré R (2011) Substituting non-natural agents in UV-protection cream by a mixture of clay with Ganoderma pfeifferi extract. Appl Clay Sci 53: 66–72.
  143. 143. Summers B (2012) Sun protection factor of South African-tested sunscreens. S Afr Med J 102 (12): 897.
  144. 144. Hoang-Minh T, Le TL, Kasbohm J, Gieré R (2010) UV-protection characteristics of some clays. Appl Clay Sci 48: 349–357.
  145. 145. Mie G (1908) Beiträge zur Optik trüber Medien speziell kolloidaler Goldlösungen. Ann Phys 25: 377–445.
  146. 146. Konert M, Vandenberghe J (1997) Comparison of laser grain size analysis with pipette and sieve analysis: a solution for the underestimation of the clay fraction. Sedimentology 44 (3): 523–535.
  147. 147. Dayet L, d’Errico F, Garcia-Moreno R (2014) Searching for consistencies in Châtelperronian pigment use. J Archaeol Sci 44: 180–193.
  148. 148. Rietveld H (1969) A profile refinement method for nuclear and magnetic structures. J Appl Crystallogr 2: 65–71.
  149. 149. Witzel RF, Burnham RW, Onley JW (1973) Threshold and suprathreshold perceptual color differences. J Opt Soc Am 63 (5): 615–625. pmid:4731340
  150. 150. Elias M, Chartier C, Prévot G, Garay H, Vignaud C (2006) The colour of ochres explained by their composition. Materials Science and Engineering 127: 70–80.
  151. 151. Barham LS (1998) Possible early pigment use in south-central Africa. Curr Anthropol 39: 703–710.
  152. 152. Hovers E, Ilani S, Bar-Yosef O, Vandermeersch B (2003) An early case of colour symbolism: Ochre use by modern humans in Qafzeh Cave. Curr Anthropol 44 (4): 491–522.
  153. 153. Salomon H, Vignaud D, Coquinot Y, Beck L, Stringer C, Strivay D, et al. (2012) Selection and heating of colouring materials in the Mousterian levels of Es-Skhul (c. 100 000 years BP, Mount Carmel, Israel). Archaeometry 54 (4): 698–722.
  154. 154. Barham L (2002) Backed tools in Middle Pleistocene central Africa and their evolutionary significance. J Hum Evol 43: 585–603. pmid:12457850
  155. 155. Tryon CA, McBrearty S (2002) Tephrostratigraphy and the Acheulean to Middle Stone Age transition in the Kapthurin Formation, Kenya. J Hum Evol 42: 211–235. pmid:11795975
  156. 156. Cruz-Uribe K, Klein RG, Avery G, Avery M, Halkett D, Hart T, et al. (2003) Excavation of buried Late Acheulean (Mid-Quaternary) land surfaces at Duinefontein 2, Western Cape Province, South Africa. J of Archaeol Sci 30: 559–575.
  157. 157. van Peer P, Fullagar R, Stokes S, Bailey RM, Moeyersons J, Steenhoudt F, et al. (2003) The Early to Middle Stone Age transition and the emergence of modern human behaviour at site 8-B-11, Sai Island, Sudan. J Hum Evol 45 (2): 187–193. pmid:14529653
  158. 158. Beaumont PB, Vogel J (2006) On a timescale for the past million years of human history in central South Africa. S Afr J Sci 102: 217–228.
  159. 159. Thouveny N, Bourlés DL, Saracco G, Carcaillet JT, Bassinot F (2008) Paleoclimatic context of geomagnetic dipole lows and excursions in the Brunhes, clue for an orbital influence on the geodynamo? Earth Planet Sci Lett 275: 269–284.
  160. 160. d’Errico F, Backwell L, Villa P, Degano I, Lucejko JK, Manford MK, et al. (2012). Early evidence of San material culture represented by organic artifacts from Border Cave, South Africa. Proc Natl Acad Sci U S A 109 (33): 13214–13219. pmid:22847420
  161. 161. Villa P, Pollarolo L, Degano I, Birolo L, Pasero M, Biagioni C, et al. (2015) A milk and ochre paint mixture used 49,000 years ago at Sibudu, South Africa. PLOS ONE
  162. 162. Knight C, Power C, Watts I (1995) The human symbolic revolution: A Darwinian account. Camb Archaeol J 5 (1): 75–114.
  163. 163. Power C, Aiello L (1997) Female proto-symbolic strategies. Hager , Routledge, London.
  164. 164. McDougall I, Brown FH, Fleagle JG (2005) Stratigraphic placement and age of modern humans from Kibish, Ethiopia. Nature 433: 733–736. pmid:15716951