1. Introduction
Fulgurites (from the Latin fulgur, meaning thunderbolt) are natural hollow glass tubes formed in silica, quartzose sand, rocks or soil by lightning strikes. They are formed when temperatures above 1800 °C, originating from a lightning bolt, instantaneously melt silica and other soil minerals on a conductive surface and grains are fused together; the fulgurite tube is the cooled product. This process occurs over a short period (c. 1 s), and leaves evidence of the lightning path and its dispersion over the surface. Fulgurites are usually formed in the topsoil and do not typically extend deeper than 40 cm into the less conductive subsoil. Instead, they run parallel to the surface of the ground (Carter et al., Reference Carter, Hargreaves, Kee, Pasek and Edwards2010a). However, in some cases, fulgurites can be several metres deep (Martin Crespo et al., Reference Martin, Lozano and Gonzalez2009).
They are identified as a variety of the mineraloid lechatelierite (an amorphous form of silicon dioxide, SiO2), although their absolute chemical composition is dependent on the physical and chemical properties of the soil (Gailliot, Reference Gailliot2016; Elmi et al., Reference Elmi, Chen, Goldsby and Gieré2017). Their classification has generally been limited to varieties that are formed below the surface of the ground, such as soil, rock, sand or clay fulgurites (Pasek et al., Reference Pasek, Block and Pasek2012). Occurrences of rock fulgurites associated with steel pylons of overhead electric transmission lines have also been reported (Kassi et al., Reference Kassi, Kasi, Friis and Kakar2013).
In contrast, exogenic fulgurites, also referred to as ‘droplet fulgurites’ or as ‘type V’ fulgurites, are a fairly new concept. This class of fulgurites refers to liquefied materials resulting from a powerful lightning strike (>100 GW) that are thrown into the atmosphere above the lightning’s point of impact (i.e. they are ejected) and solidify in the air. Exogenic fulgurites feature an amorphous and often bubbly appearance due to the rapid manner in which the airborne, liquefied materials cool down. They are generally dark green in colour (as a result of a moderate iron oxide content), differing from traditional fulgurites, which are typically carrot, brown or tan. In addition, they have a smoother, glassier texture than conventional fulgurites, which have a gritty, sandy feel to them (Pasek et al., Reference Pasek, Block and Pasek2012).
Although there are descriptions of fulgurites dating from the 12th century (such as the quartzofeld–spathic fulgurite named mazintarincan in the Lapidary of King Alfonso X, from 1250), no exogenic fulgurites had been reported until modern times. To the best of the authors’ knowledge, Mohling in 2004 was the first to document the occurrence of an exogenic fulgurite south of the town of Elko, in Elko Hills, NE Nevada, USA, in 1997 (Fig. 1a) (Mohling, Reference Mohling2004). The second documented case, reported by Walter in 2008, occurred in Oswego, New York, USA (Fig. 1b) (Walter, Reference Walter2011). The third documented occurrence was in August 2013 in a suburb of Dallas, Texas, USA (Fig. 1c).
In ancient Iberia, Celtiberians and Lusitanians revered a deity named Neto or Neton, the god of lightning, but it was the Romans who showed more religious veneration for places that had been struck by lightning bolts. It was a Roman practice to surround such sacred places with an enclosure open at the top (puteal) or to build a shrine (Macrobius and Kaster, Reference Macrobius and Kaster2011).
Tiedra (Valladolid, Castile and León, Spain) was built on the site of a pre-Roman town known as Amallobriga, subsequently settled by the Romans. At the Cerro de la Ermita archaeological site there is a religious building erected in honour of Our Lady of Tiedra Vieja that retains the Roman custom of making votive offerings. Religious amulets and abundant Hispano-Roman pottery have been discovered around the hermitage (Martín-Gil and Martín-Gil, Reference Martín-Gil and Martin-Gil2001; Sanz-Mínguez and Sobrino-González, Reference Sanz-Mínguez and Sobrino-González2013). Along with these remains, intriguing droplets and vitrified masses were also found, which were initially misinterpreted as archaeologically anomalous material and discarded. In this study, their appropriate characterization as fulgurites, in accordance with the primitive sacred origin granted to the place, is reported. It is worth noting that the droplets investigated in this work showed no evidence of magnetic behaviour, thus ruling out the possibility of their being iron foundry slags.
This is the first report on exogenic fulgurites in Spain, given that a thorough search of the relevant literature yielded only two related articles in the Iberian peninsula, one on a natural fulgurite found in Bustarviejo, Madrid, Spain (Garcia-Guinea et al., Reference Garcia-Guinea, Furio, Fernandez-Hernan, Bustillo, Crespo-Feo, Correcher, Sanchez-Muñoz, Matesanz and Gucsik2009), and the other on a fulgurite generated by a small electricity pylon in Torre de Moncorvo, Bragança, Portugal (Martin Crespo et al., Reference Martin, Lozano and Gonzalez2009), neither of which were type-V.
2. Experimental
2.a. Location
The samples under study were found near the Roman road from Tiedra to Pobladura de Sotiedra (41.6521978° N, 5.274904° W), as shown in Figure 2, in September 1988.
The geology of the surrounding area comprises layers of sand and clays classified into the Neogene period, Miocene epoch, Tortonian stage (Vindovonian) (Fig. 2d). The dominant mineral species are quartz, kaolinite, muscovite and calcite. For a thorough mineralogical and petrological study of the site, the interested reader is referred to the work by Jiménez Fuentes and García Marcos (Reference Jiménez Fuentes and García Marcos1980).
2.b. Materials
The area in which the fulgurites were found was an agricultural plot with a topsoil layer rich in ashes. The samples were located at a depth of several centimetres before being removed by a mould plough. They were protected by finely divided material. Their morphology suggests that they had barely experienced erosion and transport since their fragile surface was well preserved. The only appreciable deterioration, compared to recent fulgurites, was the loss of the peduncles, which may have occurred when they were unearthed. In fact, another fulgurite found within the contents of the ash altar at the temple of Lykaian Zeus at Mount Lykaion in Greece (Romano and Voyatzis, Reference Romano and Voyatzis2010) also appeared as well conserved as the ones reported herein.
The droplet-like exogenic fulgurites (Fig. 3a) showed a grey-green hue and ranged in size from 0.7 to 1.5 cm. All featured smooth, lustrous surfaces on all sides, suggesting that they had solidified in the air. None of them floated in water.
Apropos of the masses of the molten glass pieces (composites, Fig. 3b) close to the effusive centre, they showed vesicular bodies to which the droplets must have been attached. These ‘exogenic fulgurite-matrix’ pieces showed the vitreous aspect of droplets on their outer surface but their interior appeared as a heterogeneous mixture with fully differentiated components (phases or regions), which could be visually distinguished from one another: one of greenish appearance, one of brown hue and another entirely white. A small sample of the outer surface of one of these pieces was subjected to a carbon, hydrogen and nitrogen (CHN) analysis. The rest of the analyses were conducted using non-destructive techniques, given the categorization of the specimens as archaeological objects. Inner and outer parts of these samples are indicated as i (interior) and e (exterior), respectively.
2.c. Methods
Given the impossibility of grinding the samples to a fine powder to prepare KBr pellets, the vibrational spectra were first characterized using a Thermo Scientific (Waltham, MA, USA) Nicolet iS50 Fourier-Transform Infrared (FTIR) spectrometer equipped with an inbuilt diamond attenuated total reflection (ATR) system. The spectra were collected in the 400–4000 cm−1 region at room temperature with a 1 cm−1 spectral resolution; 64 scans were co-added and the resulting interferogram was averaged. The ATR-FTIR spectra were corrected using the advanced ATR correction algorithm (Nunn and Nishikida, Reference Nunn and Nishikida2008) available in the OMNIC™ software suite.
FT-Raman spectra in the 4000–50 cm−1 range, at 3 cm−1 spectral resolution, were acquired also on non-ground samples with a micro-Raman spectroscopy system. A Horiba (Kyoto, Japan) T64000 apparatus, with focal length 0.64 m, f/7.5 aperture, diffraction gratings with 1800 gr mm−1, single channel detection (with a R943 photomultiplier), and a N2-cooled 1” charge-coupled device (CCD) detector with 1024 × 256 pixels, was used. The 1064 nm line of Nd:YAG laser was used as excitation source and laser power was set at 100 mW. Images were taken with an Olympus BH2 (Shinjuku, Tokyo, Japan) microscope.
The X-ray powder diffractograms of the samples were obtained using a Bruker (Billerica, MA, USA) D8 Advance Bragg–Brentano diffractometer, in reflection geometry, with CuKα (λ = 1.54 Å) radiation and using crystalline silicon as a standard. X-ray generator: 40 kV, 40 mA; 1D LynxEye detector; scan mode: continuous PSD fast; 2θ = 5–80°, in 0.02° steps, 1 s/step. Diffractograms were analysed in DIFFRAC.SUITE EVA software, using the most up-to-date ICDD (International Centre for Diffraction Data) and COD (Crystallography Open Database) databases for peak identification. It should be clarified that, since the archaeological specimens could not be ground to a fine powder, samples had to be mounted in such a way that the exposed surface would be as flat as possible (at the expense of reducing the intensity of low-angle peaks and adding noise to the diffractograms, which prevented subsequent Rietveld refinement).
Major (except Na and C) and minor elements were determined on non-ground samples by using a Niton XL3t XRF Analyzer (Thermo Fisher Scientific), using the TestAll™ Geo mode for a rapid analysis of major and trace elements (i.e. an overall composition of the sample). X-ray tube: Au anode, 50 kV, 200 µA. Data were processed using the Niton Data Transfer (NDT™) PC software suite.
C content in the outer layer of one of the ‘matrix’ samples was determined by dry combustion in a CHN-2000 Analyzer (LECO Corp., Saint Joseph, MI, USA). The sample was placed in a tared tin capsule and weighted using a microbalance. The combustion temperature was set to 950 °C and ethylenediaminetetraacetic acid (CAS 60-00-4, Sigma Aldrich, ≥99 %) was used as the analytical standard.
FTIR and CHN analyses were conducted at LTI-Palencia (Universidad de Valladolid, Spain), X-ray powder diffraction (XRPD) at the TAIL-UC facility (Universidade de Coimbra, Portugal), Raman spectroscopy at the Química–Física Molecular Research Unit lab (Universidade de Coimbra, Portugal), and X-ray fluorescence (XRF) at the GeoSciences Centre (Universidade de Coimbra, Portugal).
3. Results
3.a. Vibrational characterization
3.a.1. Infrared spectroscopy
In the molten glass-matrix composites (Fig. 4a), greenish and brown phases shared FTIR bands at 946, 873 and 425 cm−1. The first band, at 946 cm−1, is characteristic of silicate and aluminate glasses containing calcium. The second band, at 873 cm−1, can be attributed either to a Si–N bond, to calcite (CaCO3), or to C60 or C59N fullerenes, present in other fulgurites reported in the literature (Daly et al., Reference Daly, Buseck, Williams and Lewis1993; Heymann, Reference Heymann1998). The third band, at 425 cm−1, may be assigned to any magnesium carbonate (magnesite, MgCO3 or other minerals of the magnesite–siderite series, MgCO3 to FeCO3) (Brusentsova et al., Reference Brusentsova, Peale, Maukonen, Harlow, Boesenberg and Ebel2010). The greenish phase showed a weak band at 1506 cm−1 that can be attributed either to a ν 3 carbonate vibration (Saikia et al., Reference Saikia, Parthasarathy, Sarmah and Baruah2008) or to a C–C stretching in an aromatic ring. The brown phase is characterized by a band at 1457 cm−1 from either a C59N azafulleren or a pressure-induced phase pC60.
The white phase (also depicted in Fig. 4a) showed FTIR bands at 1084 cm−1 (cristobalite, SiO2), at 456 cm−1 (attributed to the shock compressed quartz, lechatelierite, SiO2) and the characteristic peaks for quartz (doublet at 799 and 780 cm−1 and singlet at 695 cm−1) (Foster and Walker, Reference Foster and Walker1984). The absorption peaks at 780–800 cm−1 are due to Si–O–Si symmetrical stretching vibration (tetrahedral–tetrahedral ion vibrations), and the absorption peak at 695 cm−1 (Si–O–Si symmetrical bending vibration) arises due to octahedral site symmetry (Parthasarathy et al., Reference Parthasarathy, Kunwar and Srinivasan2001; Saikia et al., Reference Saikia, Parthasarathy, Sarmah and Baruah2008). The two bands (not shown in Fig. 4a) that appeared at 2360 and 2342 cm−1 correspond to the absorption doublet of CO2.
The FTIR spectrum of the inner part of the droplet-fulgurites (Fig. 4b) showed bands at 1456 cm−1 (compatible with N-doped graphene or fullerene allotropes C60 and C59N), at 1081 cm−1 (Si–O–Si asymmetrical stretching vibration from quartz), at 796 cm−1 (Si–O–Si symmetrical stretching vibration for cristobalite), at 780 and 799 cm−1 (Si–O–Si symmetrical stretching vibration for quartz), at 695 cm−1 (quartz), at 621 cm−1 (cristobalite) (Foster and Walker, Reference Foster and Walker1984) and at 460 cm−1 (biotite, K(Mg,Fe)3(AlSi3O10)(F,OH)2; or phlogopite, KMg3(AlSi3)O10(OH)2). The doublet of CO2 also appeared at 2360 and 2342 cm−1 (not shown).
The outer crust of the droplets (Fig. 4b) shared several bands with the spectrum of the inner part (those arising from cristobalite and quartz), and showed additional weak bands at 878 cm−1 (attributable to CaCO3, vaterite or calcite), at 680 cm−1 (pyroxenes such as hedenbergite, FeCaSi2;O6; diopside, MgCaSi2O6; or augite (Ca,Na)(Mg,Fe,Al,Ti)(Si,Al)2O6), at 670 cm−1 (magnetite feature), and at 406 cm−1 (polaronic torsion in C–N–C bonds). The 2360 and 2342 cm−1 doublet of CO2 (not shown) was also present.
A probable feature during a lightning strike is that silicon in quartz is reduced to silicon metal, but in the exogenic fulgurites studied herein the typical phonon line of silicon at 520 cm−1 and the band at 970 cm−1 (related to disorder in a glass matrix) were not detected by FTIR. Other unique species to fulgurites, including silicides such as naquite (previously named fersilicite, FeSi) and silicon carbide (SiC), and phosphides in the form of schreibersite (Fe3P (Fe,Ni)3P) and TiP (Sheffer et al., Reference Sheffer, Melosh, Jarnot and Lauretta2003), were not detected either with this characterization technique, but will be discussed in subsequent sections.
3.a.2. Raman spectroscopy
For the outer part of the lightning-bolt induced vitrified composites, weak Raman lines were recorded at 320, 395 and 670 cm−1 (diopside and augite), at 1015 cm−1 (hedenbergite), and at 1120 cm−1 and 1606 cm−1 (hydrocarbons, C-C asymmetric stretching and C=C stretching, respectively (Carter et al., Reference Carter, Hargreaves, Kee, Pasek and Edwards2010a, b)) (Fig. 5a). The inner part (Fig. 5b) showed weak lines in the 400–2750 cm−1 region at 488 cm−1 (S–S disulphide), at 1295 cm−1 (polyenic chains), at 1442 cm−1 (compatible with C70 (Jeoung et al., Reference Jeoung, Kim, Kim and Kim1995)), at 1467 cm−1 (compatible with C60 (Lefrant et al., Reference Lefrant, Faulques, Godon, Buisson, Auban-Senzier, Jerome, Fabre, Rassat, Zahab, Lambert and Bernier1993) and amorphous diamond (Prawer et al., Reference Prawer, Nugent and Jamieson1998)), at 2444 cm−1 (compatible with graphite (Essene and Fisher, Reference Essene and Fisher1986)) and at 2730 cm−1 (compatible with graphite), in addition to a sharp multiple band attributed to unidentified hydrocarbon molecules with peaks at 2852 cm−1 (C–H2 symmetric and asymmetric, respectively), at 2884 cm−1 (CH3 symmetric stretching) and at 2935 cm−1 (–CH2 and –CH3) (Fig. 5b).
For the droplet-fulgurites (Fig. 5c) a series of sharp Raman lines were observed at 249 cm−1 (Ag mode, FeSi2), 229 cm−1 and 415 cm−1 (cristobalite), at 300 and 617 cm−1 (tentatively ascribed to silicon (Temple and Hathaway, Reference Temple and Hathaway1973)), at 505 cm−1 (albite), at 669 cm−1 (pyroxenes), at 830 cm−1 (graphitic structures), at c. 1065 cm–1 (carbonates), at 1116 cm–1 (muscovite) and at 1323 cm−1 (graphite (Bokobza et al., Reference Bokobza, Bruneel and Couzi2015)).
3.b. X-ray powder diffraction
In the diffractogram of the outer part of the matrix composites (grey-green glassy phase), and specifically in sample fulg5e, characteristic patterns for diopside (COD 9005340), augite (COD 1200006), hedenbergite (COD 9010328) and quartz (COD 1011176) could be observed (Fig. 6a). In the inner part of the same sample (white phase, sample fulg5i) cristobalite (COD 9008226), quartz (COD 9012600) and calcite (COD 1010962) were identified (Fig. 6b). For the outer surface of fulg6 matrix composite, peaks associated with magnesioferrite (Mg(Fe3+)2O4, COD 9001469), diopside and, in smaller quantities, maghemite (γ-Fe2O3, PDF 39-1346), hercynite (Fe2+Al2O4, COD 9013623), mcgovernite (Mn19Zn3(AsO4)3(AsO3)(SiO4)3(OH)21, PDF 25-0531) and silicon carbide (PDF 02-1462) were indexed.
Maghemite presence was evidenced in fulg1 droplet fulgurite by peaks at 2θ angles of 30° (220), 36° (311), 43° (400) and 57° (511) (Fig. 6d). The absence of peaks at 2θ = 10.5°, 17° and 20.5° questioned the presence of C60 postulated from FTIR, but magnesioferrite, naquite silicide or fersilicite (PDF 38-1397) and silicon carbide, not identified by FTIR, were found. The diffraction pattern for outer surface of fulg2 (Fig. 6e) resembled those of fulg5e and fulg6e, with diopside, magnesioferrite and quartz as the main mineral components.
3.c. Elemental analyses
Proportions of elements (in wt %), determined by XRF, are summarized in Table 1. The high Fe proportion for sample fulg1 was in good agreement with the composition based on magnesioferrite, naquite (fersilicite) and maghemite identified by XRPD. Combinations of Si/Al, Si/Fe and Si/Mg with values of around 5.0, 3.7 and 23.0 for sample fulg2 lead to a bulk assignation to diopside, augite and magnesioferrite. Composition of samples fulg3, fulg4 (not studied by other techniques) and fulg5e (studied by FTIR, Raman and XRPD) were also compatible with pyroxenes. Fulg6 exhibited Ca and Si percentages similar to those of fulg2. It is worth noting that Fulg4 exhibited a relatively high Zr content (0.17 %).
As regards C content, a CHN analysis was carried out for the outer crust of one of the ‘matrix’ samples (Fulg5e), obtaining C, H and N percentages of 1 %, 0.2 % and 0.2 %, respectively. This composition is compatible with carbonaceous materials based in elemental forms and allotropes of elemental carbon such as graphite, C60 and carbon black.
4. Discussion
The pressure and temperature conditions during lightning strikes are extreme and exceed the capabilities of the equipment available in synthesis laboratories, making it impossible to reconstruct how the components in fulgurites are formed (in spite of the efforts to reproduce the process by using, for instance, laser ablation (Martinez-Ramirez et al., Reference Martinez-Ramirez, Diaz, Camacho and Viehland2013)). What is known is that, in a rapid discharge, which occurs in a lightning bolt, the thermal energy is first consumed by vaporizing and melting the soil and subsequently diffuses radially outward from the spark where it melts and bakes the soil (Pasek and Pasek, Reference Pasek and Pasek2017).
The presence of amorphous quartz in the composition for the white phase in the exogenic fulgurite-matrix aggregates suggests that the abrupt electrical (Joule) heating of the soil surface yields high temperatures, which produces a thin melt layer on the surface, which then cools adiabatically. Cristobalite was found in the presumably hotter regions of the fulgurite, adjacent to the glassy melt.
The cooling of the melted glass would take place in the ozonized atmosphere thereby developed and would give rise to a baked glass that appears as a dark crust in the outer part of the pieces of the matrix and in the droplets that were spewed. The glass showed a composition of pyroxene, iron oxide and silicides. The low presence of carbonaceous materials in the glass suggests a high pyrolysis temperature induced in organic substances and a rapid quenching of the melt-trapped COx gases produced during heating (see Carter et al., (Reference Carter, Hargreaves, Kee, Pasek and Edwards2010a, b) and references therein).
The core, or inner part, of the exofulgurites was found to be richer in cristobalite, while the composition of the outer parts showed higher diopside, maghemite and albite contents and more zirconium. This variation could be due to a number of reasons, including the melting processes.
Although the exogenic fulgurites found in the archaeological site of Tiedra share some components (shocked quartz and cristobalite) with the quartzofeldspatic fulgurites found in Bustarviejo (Garcia-Guinea et al., Reference Garcia-Guinea, Furio, Fernandez-Hernan, Bustillo, Crespo-Feo, Correcher, Sanchez-Muñoz, Matesanz and Gucsik2009), they do not contain the Si–Fe–Al alloys present in the latter. Naquite, iron oxides, silicon carbide and piroxenes were identified instead. The comparison with the composition of exogenic fulgurites from Elko Hills and Oswego has not been performed due to lack of data (their discoverers have not yet published any chemical analysis results). With regard to the comparison with the Texas exogenic fulgurites, in spite of the similarities in terms of appearance, there were substantial differences in terms of composition: in our preliminary analyses of one of the Dallas fulgurites (their discoverer has not yet studied them in detail) we have found quartz, calcite, feldspars and iron oxides, but not unique species such as naquite, schreibersite, SiC or TiP.
The presence of buckminsterfulleren, C60 (a naturally occurring molecule produced when a lightning bolt strikes a hardwood tree) and other carbon compounds such as graphite, reported by other authors (Essene and Fisher, Reference Essene and Fisher1986; Daly et al., Reference Daly, Buseck, Williams and Lewis1993; Frank et al., Reference Frank, Jehlička and Hamplová2006), albeit suggested on the basis of the FTIR and Raman spectra of the samples reported herein, could not be confirmed by XRPD.
Correlation between results has been carried out through a concordance analysis based on the score of the number of matches (3, 2 or 1) for the mineral species that appeared in the identification profiles obtained through each of the techniques and taking into account their elucidative power (XRD > Raman = FTIR = XRF). Accordingly, the presence of cristobalite and pyroxenes (such as diopside, hedenbergite and augite) would be a substantiated and clear result (three occurrences); quartz/lechatelierite, maghemite, diopside, magnesioferrite, naquite and C60 or C59 may be regarded as components present with a high degree of certainty (two occurrences); and calcite, albite, moscovite, biotite, phlogopite, graphite and CO2 can be deemed as suggested components (one occurrence).
Apropos of the samples’ dating, it should be clarified that, although there is no certainty on the fact that the fulgurites under study date from Roman times, there is a high probability, provided that they were found along with many fragments of terra sigillata hispanica, and without contamination from modern remains. All the findings appeared in an area of 2 m × 2 m.
Finally, in regard to the association between sacred archaeological sites and vitrified materials, previously reported for vitrified stones in some dolmen (Abrunhosa et al., Reference Abrunhosa, Gonçalves and da Cruz1995) and for the circular stone enclosures of the Iron Age known as vitrified hill-forts (Cook et al., Reference Cook, Watson, Cook, Erskine, Jacobsson, Miller and Stetkiewicz2014), it would be further supported from the finding of glassy exogenic fulgurites in this Roman archaeological site.
5. Conclusions
The characterization of exogenic fulgurites found in the archaeological site of Tiedra represents the first study on Spanish exogenic fulgurites to date. The droplet fulgurites reported herein can be described as quartzo-pyroxene materials the composition of which includes naquite and CSi silicides; pyroxenes such as augite, diopside and hedenbergite; iron oxides such as maghemite and magnesioferrite; shocked quartz and neoformed cristobalite. Highly vesicular bodies and droplets are compatible with the siliciclastic sands (Tortonian) in the area.
The presented mineralogical and chemical data may help to shed light on high-temperature processes in geological systems and dynamic rock transformation in ultrafast events, and may be useful for gaining an insight into fulgurite glass formation processes.
Author ORCIDs
Pablo Martín-Ramos, 0000-0003-2713-2786; Francisco PSC Gil, 0000-0001-7546-0288; Francisco J Martín-Gil, 0000-0003-4989-6197; Jesús Martín-Gil, 0000-0001-9921-2465
Acknowledgements
Access to the TAIL-UC facility funded under the QREN-Mais Centro project ICT-2009-02-012-1980 and to the GeoSciences Centre of the Department of Earth Sciences XRF equipment at FCTUC (Prof. Lídia M. Gil Catarino) is acknowledged. P.M.R. acknowledges the financial support of Santander Universidades through the ‘Becas Iberoamérica Jóvenes Profesores e Investigadores, España’ scholarship programme. Dr Melissa Boyd is gratefully acknowledged for her diligent proofreading of this paper.
Declaration of interest
The authors declare no conflict of interest.