A possible point of contact between cosmic ray physics and archaeology: muon absorption radiography at the Tharros Phoenician-Roman site ANNALS OF GEOPHYSICS, 60, 1, 2017, S0107; doi:10.4401/ag-7382 S0107 A possible point of contact between cosmic ray physics and archaeology: muon absorption radiography at the Tharros Phoenician-Roman site Lorenzo Bonechi1,*, Fabio Ambrosino2,3, Luigi Cimmino2,3, Raffaello D’Alessandro1,4, Nicola Mori1,4, Pasquale Noli3, Giulio Saracino2,3, Paolo Strolin2,3, Lorenzo Viliani1,4 1 Istituto Nazionale di Fisica Nucleare (INFN), Sezione di Firenze, Sesto Fiorentino, Florence, Italy 2 Università di Napoli Federico II, Dipartimento di Fisica, Naples, Italy 3 Istituto Nazionale di Fisica Nucleare (INFN), Sezione di Napoli, Naples, Italy 4 Università di Firenze, Dipartimento di Fisica e Astronomia, Sesto Fiorentino, Florence, Italy ABSTRACT Several on-going activities exploiting the potential of the muon-absorption radiography technique are focusing on the study of large volcanic struc- tures located in different geographical regions. Nonetheless, the possible application of this non-destructive surveying method to other fields is also under consideration by different groups. Looking backward to the his- tory of muon radiography, we can learn how the first important success of this technique was achieved in the archaeological field by the physicist Luis W. Alvarez in the 1960s. Further examples of possible applications can be found concerning geological and mining applications. In this paper a possible application of muon absorption radiography in the context of the archaeological research is proposed. Results of a simulation of a sim- plified case study are presented. 1. Introduction Many activities exploiting the muon absorption ra- diography (or muon transmission radiography) tech- nique are currently under consideration or development in fields other than volcanology, as partially reported in other contributions to this workshop. This technique, which can be properly contextualized in the generic framework of monitoring or imaging of material struc- tures, has the potential to be applied in many different areas. Each area of application, and often each applica- tion itself, requires an independent feasibility study for this technique to be validated. This is due to several rea- sons, which can be summarized in the following main points: 1) the different materials and material thicknesses encountered by the incident muons, 2) the muon energy range relevant for the selected case study, 3) the relative detector-target geome-try, 4) the type of detector which can be used for the measurement and 5) the purpose of the measurement itself. The first practical cases where muon absorption ra- diography was successfully applied were in the fields of geology and archaeology. In the 1950s E.P. George meas- ured the thickness of the ice burden over the Guthega- Munyang tunnel in Australia [Agostinelli et al. 2003]. A decade later L.W. Alvarez studied the interior of the second pyramid of Chephren in Egypt, looking for a hypothetic hidden chamber [Bonechi et al. 2015]. In the latter case, Alvarez was able to demonstrate the ab- sence of any large chamber hidden over the already dis- covered one. Drawing on the idea of Alvarez, in this work we propose a possible application of muon absorption ra- diography in the framework of Archaeology, based on the use of a robust, low cost and low power consump- tion particle detector. In Section 2 a description of the case study is presented; in Section 3 the preliminary re- sults of a simulation are shortly discussed; finally, in Section 4 the first activities for the production of a pro- totype detector are presented. 2. Detecting hidden tombs at Tharros The Tharros site (see Figure 1) is a beautiful Punic- Roman site located at the south end of the Sinis Penin- sula near the town of Oristano in the Sardinia Island (Italy). In 2014 we had contacts with the managers of the site who illustrated some relevant details about the an- Article history Received July 5, 2016; accepted October 5, 2016. Subject classification: Muon radiography, Muon absorption, Archaeology, Tharros necropolis. cient town. In summer, we had the possibility to visit the site to evaluate a possible usefulness of exploiting muon radiography for a survey of unknown parts of the site. Apart from the main ruins of Tharros, related to the old town itself, a large area of the peninsula hosts an important necropolis, which is still hidden under- ground for the greatest part. Only a few points have been brought to light by excavations and later covered again to preserve the status of the site. Part of these ex- cavations are still visible (see Figure 2). Surveys by means of a GPR device have been only partially conducted due to the difficulty of moving this kind of device over large areas covered by the Mediterranean scrub, which is protected by the Italian legislation. Therefore, the presence of tombs in several parts of the site (e.g. Fig- ure 3) is currently only a guess. Thanks to the muon radiography technique, the placement of a muon detector in a fixed measuring po- sition, possibly ten or fifteen metres underground or under a cliff, would give the possibility to make a sur- vey of some hundred square metre surface searching for empty cavities. The minimum volume of a hidden tomb that can be reasonably detected with this method varies from one to a few cubic metres, assuming the tomb is located a few metres underground. To assess the possibility of identifying such hypothetic hidden tombs by exploiting the information carried by the de- tected muons, a data analysis method [Anastasio et al. 2013] has been proposed, based on the comparison of the measured flux with a simulated expected flux, cal- culated for the same data acquisition time. The required simulations must take into account all the relevant aspects of the problem, from the functional shape of the muon energy/angular spectra at ground level, to its modification while traversing the material along the detector’s field of view. This method is based on the back-projection of the detected and simulated muon tracks to parallel planes at increasing distances from the detector, or to spherical surfaces centred on the detector’s centre for increasing radius. Two hit map distributions are therefore produced for each surface and then subtracted from each other. As a result of the different muon ab- sorption in the real and in the simulated cases, the re- sulting “difference map” should present regions with an accumulation or deficit of hits beyond the statistical fluctuations in correspondence of a sufficiently large unknown empty cavity or a place containing a material much heavier than the surrounding earth or rock. 3. Preliminary study with the Geant4 high energy physics simulation package Preliminary simulations of a simplified case study have been implemented based on the Geant4 simula- tion package [Agostinelli et al. 2003], a toolkit that is commonly used in high energy physics for simulating the passage of particles through matter. Details of this study are published in Bonechi et al. [2015]. A homo- geneous terrain with a flat surface has been considered, with a generic “reasonable” chemical composition. A realistic muon detector, based on the technology used for the MURAY experiment [Anastasio et al. 2013], has been positioned 7 m underground and a one cubic metre empty cavity has been placed 4 m undersoil on its field BONECHI ET AL. 2 Figure 1. View of the Tharros remains. Figure 2. Tombs found in the Tharros necropolis. Figure 3. The region of the excavation shown in Figure 2. A necrop- olis is thought to occupy large parts of this place. 3 of view, 8 m far from it. A sketch of the com-plete geom- etry implemented in Geant4 is shown in Figure 4. A realistic muon flux generator has been in- tegrated in the Geant4 code. This generator is based on experimental measurements of cosmic rays at ground level by the ADAMO experiment [Bonechi et al. 2005] in the momentum range between 100 MeV/c and 130 GeV/c, which is relevant for this application. Two different simulated data sets, with and with- out the presence of the empty cavity, have been gener- ated and the results compared (see Bonechi et al. [2015] for details). For each of these simulations a number of muon events, corresponding to a data acquisition time ranging from two weeks up to 130 days, entering the detector’s field of view and passing through a region surrounding the cavity has been generated. Muon tracks are reconstructed and the final trajectories saved to file only if crossing all tracking planes. The sets of recon- structed tracks are back projected to vertical planes or spherical surfaces and a two dimensional hit density map is produced for both simulations in correspondence of each surface, as shortly described at the end of the pre- vious section. From the comparison of the two simulations, a hit map difference is obtained for all the selected surfaces. Due to geometrical reasons, a focused image of the cav- ity is expected as the back projection surface ap- proaches the position of the cavity. In Figures 5, 6 and 7 the difference maps on spherical surfaces of radii 1 m, 4 m and 7 m are shown in the azimuth-zenith plane. These figures show how the image of the cavity gets more detailed while the radius of the back projec- tion sphere increase from 1 m to 7 m (see the RMS val- ues of the distributions). As shown in Bonechi et al. [2015], the best focusing of the image of the cavity can be obtained when the back projection sphere (or plane) crosses a region in proximity of the real position of the cavity itself. In Fig- ure 8 the dependence of the width of the recon- structed image along the azimuth axis, as function of the radius of the back projection sphere (equivalent to the distance from the centre of the detector) is shown. The minimum of this function is found for a dis- tance d = 6.1 ± 0.7 m, approximately compatible with the position of the face of the cavity closer to the detector. The angle subtended in this example by the cavity with respect to the detector’s centre is approximately 125 mrad. A detector with an angular resolution of 10÷50 mrad is therefore required for this measurement. 4. Development of a prototype borehole tracking detector Based on the experience derived from the MURAY experiment, a prototype borehole muon tracker is under development for a demonstration of the technique pro- MUON RADIOGRAPHY IN ARCHAEOLOGY Figure 4. A sketch of the geometry implemented in the GEANT4 simulation. Figure 5. Projection on a 1 m radius spherical surface. Figure 6. Projection on a 4 m radius spherical surface. Figure 7. Projection on a 7m radius spherical sur-face. posed in Bonechi et al. [2015]. This detector is a 40 cm × 40 cm × 40 cm cubic tracker made of two x and two y planes. Each plane is made of 22 scintillator bars with isosceles triangular section (Figure 9). The scintillator material is UPS-923A by De-tecRad, which has a good light yield and a fast rising time. The bar section has a 40 mm base and 20 mm height and the total length of each bar is 40 cm. Each bar is enveloped in an aluminized mylar film to guarantee the optical insulation of the scintillator from the external environment (Figure 10). The triangular bars are assembled one with the base facing upwards and the adjacent one with the base facing downwards, in such a way to fill a 20 mm height parallelepiped volume (with oblique faces at two ends). Signals from the two adjacent bars crossed by a muon can be used to determine the crossing point of the muon tracks with a better resolution with respect to the simple digital algorithm. For the single tracking plane, made of 22 bars, we expect a spatial resolution around 3÷4 mm, according to previous tests with similar bars [Carlà 2009/2010]. Considering two planes at 40 cm distance, it means therefore an angular resolution of 10 mrad or better. Each bar is readout by two silicon photomultiplier (SiPM), one on each side. The selected model for this de- vice is ASD-NUV4S-P by AdvanSiD, which is a 4×4 mm2 device optimized for near ultra violet light detection. A transparent epoxy layer protects the surface of each SiPM, therefore the optical contact between SiPM and scintilla- tor bar can be achieved by means of a thin layer of opti- cal grease directly applied to the surface of the SiPM. These sensors are kept faced to the triangular faces of the bars thorough dedicated plastic supports produced at a 3D printer. A photo of such a plastic support is shown in Figure 11. These parts are fixed head on to the bars thanks to standard double sided tape. Some pieces of kapton tape is finally applied for safety to fix the support laterally. The readout electronics is currently under devel- opment. It is based on the EASIROC chip by Omega, a 32 channels chip dedicated to readout SiPM detectors. BONECHI ET AL. 4 Figure 8. Dependence of the RMS value associ-ated to the image of the cavity as a function of the radius of the back projection sphere. A minimum can be identified for a radius around 6-7 m. Figure 9. One scintillator bar used for the con-struction of the pro- totype muon tracker. Figure 10. One scintillator bar enveloped in the aluminized mylar film. Figure 11. Plastic support use to fix the SiPM sensors to the scintil- lator bars. At the centre a slot used to host the sensor is visible, with a small hole to drive the signal to the readout electronics. 5 The same chip is used for the MURAY and MURAVES projects. 5. Conclusions The possibility to exploit the potential of muon ab- sorption radiography for surveys in the archaeological field is under consideration. In particular, the application of this technique for the study of the necropolis of the Tharros site in the Sardinia island is under considera- tion. For this purpose a prototype borehole muon tracker based on scintillators and Silicon Photomulti- pliers in under development. The readout electronics will be similar the MURAY and MURAVES projects. The first test of the technique is foreseen for 2016 in a controlled underground environment at the INFN Unit of Florence (Italy). Acknowledgements. We would like to thank prof. Anna Chiara Fariselli, who has the ministerial concession for archaeological research at Tharros, dr. Federica Boschi from University of Bologna and dr. Roberto Carrus, president of Cooperativa “Penisola del Sinis”, for their availability for discussions and for the invitation to visit the site of Tharros. We would like also to thank Mr. Massimo Falorsi, head of the workshop of the Department of Physics and Astronomy of University of Florence, and Mr. Marco Manetti of INFN Unit of Florence for the production of scintillators and mechanics. References Agostinelli, S., et al. (2003). Geant4 - a simulation toolkit, Nucl. Instrum. Meth. A, 506, 250-303. Anastasio, A., et al. (2013). The MURAY experiment. An application of SiPM technology to the under- standing of volcanic phenomena, Nucl. Instrum. Meth. A, 718, 134-137. Bonechi, L., M. Bongi, D. Fedele, M. Grandi, S. Riccia- rini and E. Vannuccini (2005). Development of the ADAMO detector: test with cosmic rays at different zenith angles, Proceedings of the 29th International Cosmic Ray Conference, 9, 283-286. Bonechi, L., R. D’Alessandro, N. Mori and L. Viliani (2015). A projective reconstruction method of un- derground or hidden structures using atmospheric muon absorption data, J. Instrum., 10, P02003. Carlà, M. (2009/2010). Misura della risoluzione spaziale di un odoscopio realizzato con scintillatori a sezione triangolare, Bachelor’s Degree in Physics, Academic Year 2009/2010, University of Florence, Italy (in Italian). *Corresponding author: Lorenzo Bonechi, Istituto Nazionale di Fisica Nucleare (INFN), Sezione di Firenze, Sesto Fiorentino, Florence, Italy; email: lorenzo.bonechi@fi.infn.it. © 2017 by the Istituto Nazionale di Geofisica e Vulcanologia. All rights reserved. MUON RADIOGRAPHY IN ARCHAEOLOGY << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles false /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Warning /CompatibilityLevel 1.3 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.1000 /ColorConversionStrategy /LeaveColorUnchanged /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams true /MaxSubsetPct 100 /Optimize false /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments false /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile (None) /AlwaysEmbed [ true /AndaleMono /Apple-Chancery /Arial-Black /Arial-BoldItalicMT /Arial-BoldMT /Arial-ItalicMT /ArialMT /CapitalsRegular /Charcoal /Chicago /ComicSansMS /ComicSansMS-Bold /Courier /Courier-Bold /CourierNewPS-BoldItalicMT /CourierNewPS-BoldMT /CourierNewPS-ItalicMT /CourierNewPSMT /GadgetRegular /Geneva /Georgia /Georgia-Bold /Georgia-BoldItalic /Georgia-Italic /Helvetica /Helvetica-Bold /HelveticaInserat-Roman /HoeflerText-Black /HoeflerText-BlackItalic /HoeflerText-Italic /HoeflerText-Ornaments /HoeflerText-Regular /Impact /Monaco /NewYork /Palatino-Bold /Palatino-BoldItalic /Palatino-Italic /Palatino-Roman /SandRegular /Skia-Regular /Symbol /TechnoRegular /TextileRegular /Times-Bold /Times-BoldItalic /Times-Italic /Times-Roman /TimesNewRomanPS-BoldItalicMT /TimesNewRomanPS-BoldMT /TimesNewRomanPS-ItalicMT /TimesNewRomanPSMT /Trebuchet-BoldItalic /TrebuchetMS /TrebuchetMS-Bold /TrebuchetMS-Italic /Verdana /Verdana-Bold /Verdana-BoldItalic /Verdana-Italic /Webdings ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 150 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.10000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 150 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.10000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.08250 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile (None) /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName (http://www.color.org) /PDFXTrapped /Unknown /CreateJDFFile false /SyntheticBoldness 1.000000 /Description << /ENU (Use these settings to create PDF documents with higher image resolution for high quality pre-press printing. The PDF documents can be opened with Acrobat and Reader 5.0 and later. These settings require font embedding.) /JPN /FRA /DEU /PTB /DAN /NLD /ESP /SUO /NOR /SVE /KOR /CHS /CHT /ITA >> >> setdistillerparams << /HWResolution [2400 2400] /PageSize [595.000 842.000] >> setpagedevice