The authors present here a new method to image reinforcement elements inside thick structures and the results of a demonstration measurement performed on a mock-up wall built at Los Alamos National Laboratory. The method, referred to as “multiple scattering muon radiography”, relies on the use of cosmic-ray muons as probes. The work described in this article was performed to prove the viability of the technique as a means to image the interior of the dome of Florence Cathedral Santa Maria del Fiore, one of the UNESCO World Heritage sites and among the highest profile buildings in existence. Its result shows the effectiveness of the technique as a tool to radiograph thick structures and image denser object inside them.

The cathedral of Santa Maria del Fiore (Saint Mary of the Flower), shown in Figure 1, was founded in 1296 and is the main church of Florence. The church of Florence is one of the UNESCO World Heritage sites. Its majestic dome (“Cupola” in Italian) was constructed between 1420 and 1436 under the direction of Filippo Brunelleschi and now dominates the Florence skyline.

FIG. 1.

Pictures of the Cathedral of Santa Maria del Fiore, images by authors. Left insert: the facade. Right insert: part of the main nave with Brunelleschi’s cupola above the transept.

FIG. 1.

Pictures of the Cathedral of Santa Maria del Fiore, images by authors. Left insert: the facade. Right insert: part of the main nave with Brunelleschi’s cupola above the transept.

Close modal

The Cupola consists of two shells, with the inner shell 2.25 m and the outer shell 0.8 cm thick at the base, separated by a 1.2 m and connected by spurs or buttresses.

It is built out of large blocks of sandstone up to an elevation of about 7 m above its base. From here to the base of the lantern, the masonry is made of clay bricks of different sizes and arranged according to different patterns with the purpose of ensuring the Cupola’s structural stability. At the time of its construction, the Cupola was the largest dome ever built and it is still not clear how Brunelleschi managed to build it without a temporary supporting structure (centering). Three pairs of macigno (large stones) chains are believed to reinforce the Cupola, although only the existence of the lowest pair has been proven. A wooden chain was also installed 7.75 m above the base of the Cupola. The base of the dome is octagonal in shape, the apothem of the polygon describing its inner base surface measures 20.46 m.1 

According to some scholars, there are iron chains inside the masonry, although investigations with metal detectors have failed to provide conclusive evidence. What is certain is that thousands of Florentine pounds of iron (1 Florentine pound = 0.34 kg) were purchased during the construction of the Dome though only a small fraction of it is currently accounted for.

The Dome has been affected by cracks for centuries; some of them are currently 6-8 cm wide. It is likely that they appeared shortly after the completion of the Cupola, although they were first mentioned in 1639. The location of the main cracks is shown in Figure 2 and Figure 3. The width of the cracks increases at a rate of 7.5 mm/century.1 

FIG. 2.

plan (left) and cross-section (right) of Santa Maria del Fiore showing the location of the main cracks. Reprinted from Ref. 1 with permission.

FIG. 2.

plan (left) and cross-section (right) of Santa Maria del Fiore showing the location of the main cracks. Reprinted from Ref. 1 with permission.

Close modal
FIG. 3.

Two dimensional representation of the dome showing the main cracks. Reprinted from Ref. 1 with permission.

FIG. 3.

Two dimensional representation of the dome showing the main cracks. Reprinted from Ref. 1 with permission.

Close modal

Detailed knowledge of the structure of the Cupola would benefit the finite element calculation models that are used to evaluate its behavior under static and dynamic (earthquake) conditions. Only the parts of the monument that are directly visible have been extensively studied because the techniques capable of surveilling the interior structure are limited. Brunelleschi was very secretive about his work and purposely did not leave technical drawings.

The authors realized that multiple scattering muon radiography could be used to image the inside of the Dome’s walls. Muons are charged leptons with a mass of 105.6 MeV/c2, about 200 times the mass of an electron. Cosmic-ray muons are naturally produced in the atmosphere by the interactions of primary cosmic rays, mostly protons, with the nuclei present in the atmosphere itself. They reach the surface of the Earth at a rate of about 1/cm2/min with a mean energy of ∼4 GeV.2 They have no hadronic interaction with nucleons, and their relatively large mass limits energy loss due to bremsstrahlung radiation. These properties allow energetic muons to penetrate large amounts of material that are inaccessible to other particles. Muons do possess an electric charge and undergo Coulomb scattering off nuclei as they pass through matter.

The first use of cosmic ray muons for radiography was in 1955, when George measured muon attenuation to determine the overburden of rock above a tunnel.3 This was followed by Alvarez et al., who used this method to confirm the Second Pyramid of Giza did not contain any undiscovered chambers,4 and more recently by several groups examining geologic features.5–9 

A different method, developed at Los Alamos National Laboratory, uses measurements of the multiple scattering angle of individual muons passing through an object to create tomographic images of the object’s interior structure.10 This technique was originally developed to inspect cargo containers for illicit trafficking of nuclear material,11,12 and has since been applied to detection of Special Nuclear Materials,13 industrial corrosion,14 nuclear reactors,15–17 and will be used to determine the condition of the damaged cores of the Fukushima Daiichi nuclear power plant.18,19

The method is based on the fact that muons, when travelling through materials undergo many single Coulomb scatterings with the charged atomic nuclei. As a result, they are deviated from their initial trajectory and exit the material with a different direction from the original one.

The theory of multiple coulomb scattering was developed by Moliere and Bethe.20–22 It predicts, for the muons that go through an object, a Gaussian polar angle distribution with tails. The dominant part of such distribution can be described by:

d N d θ = 1 2 π σ 2 e θ 2 2 θ 0 2
(1)

where θ is the polar angle and θ0 is given by

θ 0 = 14 . 1 MeV p β l X 0
(2)

with p being the muon momentum, β its velocity, and X0 the radiation length for the material. By measuring the distribution of cosmic muons scattering angles after they traverse an object of interest, an image of the object itself can be reconstructed as described in the paper by Borozdin et al.18 Since muons naturally and constantly shower the Earth, no additional radiation source is needed, making this technique completely safe; this avoids the complications related to the handling of artificial radiation sources or particle accelerators.

By using multiple scattering muon radiography on the inner wall of the Florence Dome, it is possible to:

  • Identify and image iron elements inside the masonry of the dome.

  • Determine a crack’s profile inside the wall, sheding light on how the wall itself is built. In fact, cracks propagate along weakness in the masonry, normally along the mortar as compared to the bricks. This information in particular could tell whether the inner wall is entirely made of bricks or if rubble masonry was used.

  • Image the transition between the stone and the brick masonry at the base of the Dome.

In order to prove the feasibility of this measurement, a demonstration measurement was performed at LANL during the summer of 2015.

If any reinforcement structure is present in the masonry of the Cupola’s wall, it most likely has to be located inside the inner shell. This shell is much thicker than the outer shell, and thus plays a more important role in sustaining the whole structure.

A concrete wall having the same thickness, in radiation lengths, as the inner wall of the Dome was built. The masonry of the Dome is mostly made of clay bricks. In order to calculate their radiation length a chemical composition and a density indicative of clay bricks from the same historical period and geographical area was obtained from an analysis of ancient bricks by Fernandes et al.23 The bricks were assumed to have the following chemical composition, the percentages being weight fractions: 55% SiO2, 25% Al2O3, 10% FeO2, 5% K2O, 5% H2O. The density was chosen to be 1.6 g/cm3, a typical value for the clay bricks manufactured in the same period of time and geographical area. The mock-up wall was built from standard radiation shielding concrete bricks whose radiation length was 26.7 g/cm2 based on the estimates provided in the Review of Particle Physics2 and whose density was measured to be (2.13±0.13) g/cm3. From these values, the thickness of a concrete wall equivalent to 2.25 m of clay was calculated to be (1.74 ±0.01) m. Since the size of the available bricks was not an integer fraction of 1.74 m, the thickness of the actual wall was actually (1.73±0.01) m.

Three iron bars were placed inside this mock-up wall, in a dedicated empty space purposely left in the middle of it. The mock-up wall is shown in Figure 4, the iron bars are shown in Figure 5.

FIG. 4.

The mock-up wall after its completion. Image by authors.

FIG. 4.

The mock-up wall after its completion. Image by authors.

Close modal
FIG. 5.

The iron bars inside the wall. Image by authors.

FIG. 5.

The iron bars inside the wall. Image by authors.

Close modal

The bars had rectangular or square cross-sections with the following dimensions:

  • 4.76 cm x 5 cm.

  • 2 cm x 3 cm, the bars in the Cupola wall are believed to be this size.24 

  • 10 cm x 10 cm.

The LANL Mini Muon Tracker (MMT),25 was used to perform the demonstration measurement. The detector consists of two modules made of 576 drift tubes arranged in planes. Each of the two modules consists of 6 planes of drift tubes, three of which are oriented along the horizontal X direction, and three of which are oriented along the horizontal Y direction, perpendicular to X. Each module can independently track cosmic rays and the trajectory information can be used to generate images of the object located between them.

The two modules of the MMT were deployed on opposite sides of the wall as shown in Figure 6.

FIG. 6.

The Mini Muon Tracker taking data around the mock up of the Cupola wall. The lower detector is still enclosed in a shed after a previous measurement campaign during which it needed to ne weather-proof. Image by authors.

FIG. 6.

The Mini Muon Tracker taking data around the mock up of the Cupola wall. The lower detector is still enclosed in a shed after a previous measurement campaign during which it needed to ne weather-proof. Image by authors.

Close modal

The rate of cosmic-ray muons increases as their direction approaches the vertical. For this reason, one of the two tracking detectors was raised by 121.92 cm (4 ft) to increase the counting rate and reduce the data taking time. We collected data for 35 days, during which we collected a total of 3,112,277 muon tracks through both detectors.

An image of the vertical plane passing through the iron bars was reconstructed using the algorithm described in the paper by Morris et al.26 and the resulting image is shown in the left insert in Figure 7.

FIG. 7.

Left: Image of the three iron bars described in the text inside the concrete mock-up wall, obtained using 35 days of data collected with the MMT. Right: Image of the same three bars obtained from the output of a Monte Carlo simulation.

FIG. 7.

Left: Image of the three iron bars described in the text inside the concrete mock-up wall, obtained using 35 days of data collected with the MMT. Right: Image of the same three bars obtained from the output of a Monte Carlo simulation.

Close modal

The pixels in the image have 1 cm x 1 cm size and their intensity is proportional to the areal density parameter extracted from the fit.10 

All the iron bars are visible, including the thinnest one. For comparison, a Monte Carlo simulation of the same experimental set-up was performed with GEANT427 and its data was analyzed with the same analysis technique. The number of simulated tracks going through both detectors used in the analysis was the same as the number of through going tracks collected in the real experiment. The areal density image obtained from the Monte Carlo simulation is also shown in Figure 7. Both the three bars and the concrete wall are uniform and homogeneous along the vertical direction. The 2D images shown in Figure 7 were therefore integrated along the vertical direction, and the 1D plot shown in Figure 8 was obtained.

FIG. 8.

The integral of the areal density obtained from the data (blue) and from the Monte Carlo simulation described in the text (red).

FIG. 8.

The integral of the areal density obtained from the data (blue) and from the Monte Carlo simulation described in the text (red).

Close modal

In order to better reproduce the experimental data, the Monte Carlo image had to be scaled by a factor of 1/1.29. The authors believe this to be due to differences between the muon momentum spectrum used in the simulation and the actual muon momentum at the measurement site, as already noticed in one of their previous works.28 

The contrast and the resolution seen in the two images are very similar. A slight difference between the two histograms in Figure 8 can be noticed at the extremes of the horizontal scale (x < 360 mm and x > 960 mm). This is probably due to bias in the experimental data in proximity with the edges of the detector’s field of view.

The good agreement between the data and the simulation demonstrates that the technique is well understood and that the performances of future measurements can be confidently predicted.

The demonstration measurement discussed here was performed in Los Alamos, NM, USA, at an altitude of about 2,230 m above sea level, to be compared with the 50 m elevation of Florence, Italy.

Since the near horizontal muon flux at sea level is about 5 % lower2 than it is at ∼2,200 m elevation a measurement at sea level, in general, would take 5% longer than the same measurement performed at 2,200 m altitude.

The authors focused here on the problem of locating iron reinforcement structures inside a thick wall. Showing the potential of muon radiography to image the profile of a crack inside a thick wall requires a different mock-up wall, with a 6-8 cm wide crack in it. Multiple measurements need to be taken at different angles around the crack to obtain a 3D tomographic image, and further development on software to reconstruct the actual tomographic image is also necessary. The authors are currently working on this problem, which will be the subject of a dedicated publication.

The authors have shown that multiple scattering muon radiography can effectively be used to locate and image iron bars inside of a 173 cm thick concrete wall, making the method ideal for investigating the interior of Brunelleschi’s dome.

This method can, in general, be used to radiograph thick structures, imaging denser objects inside them (reinforcement bars being an outstanding example). Applications include, but are not limited to, the monitoring of building and civil constructions (e.g. dams, bridges) to assess their structural stability and the imaging of the interior of thick archaeological artifacts.

This work was funded by the laboratory Directed Research and Development program at Los Alamos National Laboratory.

The authors would like to thank Prof. Carlo Blasi, Dr. Michele Fanelli and Dr. Federica Ottoni for the useful discussion and advice.

1.
M.
Fanelli
and
G.
Fanelli
,
Brunelleschi’s Cupola. Past and Present of an Architectural Masterpiece
(
Mandragora
,
Florence
,
2004
).
2.
K. A.
Olive
 et al, “
The review of Particle Physics
,”
Chinese Phys. C
38
, no. 090001 (
2014
).
3.
E. P.
George
, “
Cosmi Rays Measure Overburden of Tunnel
,”
Commonw. Eng.
455
457
(
1955
).
4.
L.
Alvarez
, “
Search for Hidden Chambers in the Pyramids
,”
Science (80-. )
167
,
832
839
(
1970
).
5.
K.
Nagamine
,
M.
Iwasaki
,
K.
Shimomura
, and
K.
Ishida
, “
Method of probing inner-structure of geophysical substance with the horizontal cosmic-ray muons and possible application to volcanic eruption prediction.
,”
Nucl. Instruments Methods Phys. Res. Sect. A Accel. Spectrometers, Detect. Assoc. Equip.
356
(
2
),
585
595
(
1995
).
6.
N.
Lesparre
 et al, “
Geophysical muon imaging: feasibility and limits
,”
Geophys. J. Int.
183
,
1348
1361
(
2010
).
7.
L.
Olah
 et al, “
Cosmic Muon Detection for Geophysical Applications
,”
Adv. High Energy Phys.
vol. Article ID (
2013
).
8.
C.
Carloganu
 et al, “
Towards a Muon Radiography of the Puy de Dome
,”
Geosci. Methods, Instrum. Data Syst.
2
(
1
),
55
60
(
2013
).
9.
F.
Ambrosino
 et al, “
The MU-RAY project: detector technology and first data from Mt. Vesuvius.
,”
J. Instrum.
9
(
2014
).
10.
L.
Schultz
 et al, “
Image reconstruction and material Z discrimination via cosmic ray muon radiography.
,”
Nucl. Instruments Methods Phys. Res. Sect. A Accel. Spectrometers, Detect. Assoc. Equip.
519
(
3
),
687
694
(
2004
).
11.
K. N.
Borozdin
 et al, “
Surveillance: Radiographic imaging with cosmic-ray muons
,”
Nature
422
(
6929
), (
2003
).
12.
W.
Priedhorsky
 et al, “
Detection of high-Z objects using multiple scattering of cosmic ray muons.
,”
Rev. Sci. Instrum.
74
(
10
), (
2003
).
13.
E.
Guardincerri
 et al, “
Detecting special nuclear material using muon-induced neutron emission.
,”
Nucl. Instruments Methods Phys. Res. Sect. A Accel. Spectrometers, Detect. Assoc. Equip.
789
,
109
113
(
2015
).
14.
J. M.
Durham
 et al, “
Tests of cosmic ray radiography for power industry applications.
,”
AIP Adv.
5
,
184909
(
2015
).
15.
T.
Sugita
 et al, “
Cosmic-ray muon radiography of UO2 fuel assembly.
,”
J. Nucl. Sci. Technol.
51
(
7–8
),
1024
1031
(
2014
).
16.
C. L.
Morris
 et al, “
Analysis of muon radiography of the Toshiba nuclear critical assembly reactor.
,”
Appl. Phys. Lett.
104
,
024110
(
2014
).
17.
J. O.
Perry
 et al, “
Imaging a nuclear reactor using cosmic ray muons
,”
J. Appl. Phys.
113
,
184909
(
2013
).
18.
K.
Borozdin
 et al, “
Cosmic Ray Radiography of the Damaged Cores of the Fukushima Reactors.
,”
AIP Adv.
2
,
042128
(
2012
).
19.
H.
Miyadera
 et al, “
Imaging Fukushima Daiichi reactors with muon.
,”
AIP Adv.
3
,
052133
(
2013
).
20.
H. A.
Bethe
, “
Molière’s Theory of Multiple Scattering
,”
Phys. Rev.
89
,
1256
(
1953
).
21.
G.
Moliere
, “
No Title
,”
J. Phys. Sci.
2
,
133
(
1947
).
22.
G.
Moliere
, “
No
,”
J. Phys. Sci.
2
,
78
(
1948
).
23.
F. M.
Fernandes
 et al, “
Ancient Clay Bricks: Manufacture and Properties.
,” in
Materials, Technologies and Practice in Historic Heritage
, edited by
M. B. D.
(
Springer Science+ Business Media
,
2009
).
24.
C.
Blasi
, Private communication,
Florence
(
2015
).
25.
J. O.
Perry
, “
Advanced Applications Of Cosmic-Ray Muon Radiography
,”
The University of New Mexico
(
2013
).
26.
C. L.
Morris
 et al, “
A new method for imaging nuclear threats using cosmic ray muons.
,”
AIP Adv.
3
,
082128
(
2013
).
27.
S.
Agostinelli
 et al, “
GEANT4 - a simulation toolkit.
,”
Nucl. Instruments Methods Phys. Res. Sect. A Accel. Spectrometers, Detect. Assoc. Equip.
506
(
3
),
250
303
(
2003
).
28.
J. M.
Durham
 et al, “
Cosmic Ray Muon Imaging of Spent Nuclear Fuel in Dry Storage Casks.
,”
J. Nucl. Mater. Manag.
(
2015
).