Assessment of radiological hazards from soil samples in the Northeastern area of Burkina Faso

The activity concentrations of main naturally occurring radionuclides were determined with the aim of assessing the radiation hazards in the Northeastern part of Burkina Faso. Soil samples were taken and analyzed by the gamma-ray spectrometry method. The ranges of specific activity of 238U, 232Th and40K are 15.25±2.59Bq·kg-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$15.25 \pm 2.59 \, \text{Bq} \cdot {\text{kg}^{-1}}$$\end{document} to 109.57±3.21Bq·kg-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$109.57 \pm 3.21 \, \text{Bq} \cdot {\text{kg}^{-1}}$$\end{document}, 19.56±1.65Bq·kg-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$19.56 \pm 1.65 \, \text{Bq} \cdot {\text{kg}^{-1}}$$\end{document} to 44.88±2.49Bq·kg-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$44.88 \pm 2.49 \text{Bq} \cdot {\text{kg}^{-1}}$$\end{document} and 125.74±4.40Bq·kg-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$125.74 \pm 4.40 \, \text{Bq} \cdot {\text{kg}^{-1}}$$\end{document} to 705.85±10.79Bq·kg-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$705.85 \pm 10.79 \, \text{Bq} \cdot {\text{kg}^{-1}}$$\end{document} respectively. Three high background radioactivity areas have been identified for Uranium, Thorium and Potassium on different places in the study area. The radiological hazards indices vary from 62.87Bq·kg-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$62.87 \, \text{Bq} \cdot {\text{kg}^{-1}}$$\end{document} to 189.47Bq·kg-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$189.47 \, \text{Bq} \cdot {\text{kg}^{-1}}$$\end{document} for the Radium equivalent activity (Raeq\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Ra_{eq}$$\end{document}), 29.5ηGy·h-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$29.5 \, \eta \text{Gy} \cdot {\text{h}^{-1}}$$\end{document} to 86.65ηGy·h-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$86.65 \, \eta \text{Gy} \cdot {\text{h}^{-1}}$$\end{document} for the absorbed dose rate (DR\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D_R$$\end{document}), 0.18mSv·y-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.18 \, \text{mSv} \cdot {\text{y}^{-1}}$$\end{document} to 0.53mSv·y-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.53 \, \text{mSv} \cdot {\text{y}^{-1}}$$\end{document} for the annual effective dose rate (AEDE), 0.21 to 0.81 for the internal hazard index (Hin\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_{in}$$\end{document}), 0.17 to 0.51 for the external hazard index (Hex\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_{ex}$$\end{document}), 1.27×10-4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.27 \times 10^{-4}$$\end{document} to 3.73×10-4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3.73 \times 10^{-4}$$\end{document} for the Excess Lifetime Cancer Risk (ELCR) and 209.01μ0.18Sv·y-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$209.01{\mkern 1mu} \mu \, 0.18{\mkern 1mu} Sv \cdot y^{{ - 1}}$$\end{document} to 594.21μSv·y-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$594.21 \, \mu \text{Sv} \cdot {\text{y}^{-1}}$$\end{document} for the annual gonadal dose equivalent. The average values of indices are well below their permissible limit. However, at Niapsi the obtained values for absorbed dose rate, annual effective dose rate and excess lifetime cancer risk are little above their permissible limit.


Introduction
Radioactivity is a natural phenomenon manifested by the emission of ionizing radiations. These radiations are imperceptible by human senses and can only be detected by special devices. Radioactivity is omnipresent and it is found everywhere in different concentrations levels in the environment, drinking water, foodstuffs, air, soil and even inside the human body. When the radioactivity level is high enough, it can induce harmful effects on human health. The main elements responsible for radioactivity are essentially 238 U , 235 U , 232 Th and 40 K . These terrestrial radionuclides have been there since the formation of Earth, 4.5 billion years ago. They are unevenly distributed in space, depending on the geological nature of the soil and they are the main source of exposition of human body to ionizing radiation. There are two kinds of exposure to ionizing radiation: external exposure and internal exposure. External exposure is manifested by the interaction of gamma, alpha or beta radiation emitted by a radionuclide with the human body. The gammaray radiation is very dangerous and difficult to stop. For this, several centimeters of lead or several meters of concrete are needed to absorb it. Alpha radiation is also harmful to health but, fortunately, it has a very low penetrating power and can be stopped with a simple sheet of paper or a few centimeters of air. The beta radiation is quite penetrating and one need several meters of air or a few millimeters of aluminum to stop it. The internal exposure occurs when the subject incorporates by inhalation or by ingestion a radionuclide which decays inside the organism. Radioactivity could be harmful for human health and there are some places where the radioactivity is significantly higher than normal areas. In Burkina Faso, an airborne study revealed the presence of uranium anomalies in the Northeastern part of the country [1]. The airborne study was carried out at a given altitude above the ground and, it has not been possible to determine whether the anomaly presented radiological risks or not. In this part of the country, gold mining activities which lead to the release of soil dust from pits, are widely practiced. It would be important to perform studies in order to know the radiological impact of these activities, on the health of the populations. In order to bring some precisions on the anomaly, radioactivity has been measured on the ground through soil samples. In the present study, the concentration of main naturally occurring radionuclides is determined and the radiation hazards associated with the presence of these radionuclides is assessed.

Study area
The present study was performed in the NorthEastern part of Burkina Faso. Samples have been taken in the villages of Tangarce, Soula, Pougbonne and Koulfou located in the north part of Gnagna in Bogande ( N 13.00242 • , W 0.13644 • ) and the villages of Niapsi (located in Solhan) and Gourgnel located in Sebba ( N 13.43669 • , E 0.53598 • ). The climate is Sudano-Sahelian and composed of two seasons: A dry season which runs from November to May and a rainy season from June to October. The dry season is characterized by a cold period from December to February accompanied by a dry wind (Harmattan) which takes place in the Sahara. The dry season is dominated by heat which starts in March. During this season, extreme temperatures are reached (45 °C in April and 15 °C in January). In addition, the average falling rainwater during the rainy season is from 500 to 800 mm.
Previously, studies have been performed in the study area, in order to highlight the uranium anomalies first and then, to determine the gamma dose rate level at one meter above the ground for various depths by using Iden-tiFinder device [2,3]. These works allowed to identify some relatively high background radioactivity points where it has been recommended to perform radioanalytical investigations. For more precision and in order to know the nature and amount of radionuclides responsible of these anomalies, the present study is undertaken, and soil samples have been taken to the laboratory for analysing by the gamma-ray spectrometry method [4,5].

Sample collection and preparation
The identiFINDER R-400 (model idF2-ULK-NGH) provided with a room temperature detector of sodium iodide (NaI), is used to get an estimation of the radioactivity level in the study area. Then, sample has been collected from some geographic points where radioactivity is relatively high by digging the ground with a manual core drill (collector) fitted with three compartments which must be screwed together before using : A first compartment of the collector is an approximately twenty (20) centimeters high cylinder. The base of this cylinder has been cut to obtain interlocking hooks allowing it to easily sink into the ground following an applied rotation and to be able to carry the soil between the teeth of the hooks. A second compartment is a rod attached to the top of the first compartment. It connects the sample collector to the sleeve. A third compartment is a sleeve connected to the rod. A physical manual force applied to this part of the collector induce a rotation of the cylinder which lead to the extraction of soil sample.
In the study area, the identiFINDER R-400 is turned on to look for a relatively high background coordinate points where gamma dose rate is higher than 100 Sv ⋅ h −1 . Once a point is located, sampling points are randomly chosen around this point and its GPS (model Garmin eTrex 10) coordinates are noted first then, after having assembled the various compartments of the collector, a stratum of soil is taken at the depth of up to 20 cm of the ground [3]. Ten (10) soil samples were taken and analyzed by gammaray spectrometry respectively from the localities of Soula, Tangarce, Pougbonne, Koulfou, Niapsi and Gourgnel (see figure 1). Samples were put in a plastic bags which were stapled with a piece of coded paper and then, an associated sample collection form is filled. Collected samples have been prepared from April to June of 2019. The dried samples have been sieved, weighted and put inside 100 cm 3 bottle for at least four (04) weeks to reach the secular equilibrium of 238 U and 232 Th with their respective daughter progenies.

Radioactivity measurement
A High Purity Germanium (HPGe) detector model (GC4010) based gamma-ray spectrometer was employed for samples analysis. The resolution of the spectrometry system was 1.8 KeV at 1330 KeV gamma-ray line of 60 Co and 0.875 KeV at 122 KeV energy line. This model of detector is a coaxial one of 60.9 mm diameter and 63.3 mm height, equiped with a hybrid cooling system by a compressor that condenses nitrogen. In addition to the detector, the system consists of a preamplifier, an analog to digital converter (ADC), a multi-channel analyzer (MCA) consisting of 8192 channels and a display screen.
In order to provide accurate analysis of results, the spectrometer system was adequately calibrated in energy, resolution and efficiency. A spectrum is aquired for 216000 s each week to get the background which will be substracted by sample spectrum. The sample is placed on the detector for 24 hours for spectrum acquisition. There are basically two methods for measuring radioactivity by gamma spectrometry. The direct method, which exploits the gamma radiation emitted by the deexcitation of the radioactive nuclei one to analyze. However, not all nuclei emit gamma radiation. This is the case for example with 232 Th and 238 U . For this, we have recourse to the indirect method which consists in exploiting the secular balance between these nuclei and their descendants. Thus, to determine the activity of 238 U , we considered, in its progeny, radionuclides such as 234 Th The efficiency of the analysis procedure was tested with a quality control sample. This is the sample3 (IAEA-TEL-2019-03) provided by the International Atomic Energy Agency for the proficiency test. The results obtained are presented in Table 1. A good approximation is obtained between the measured values and the target values. The minimum detectable activity (MDA) of the radionuclides 238 U , 232 Th and 40 K is calculated and is 3.77 Bq ⋅ Kg −1 , 5.35 Bq ⋅ Kg −1 and 9.37 Bq ⋅ Kg −1 respectively. The Radium equivalent ( Ra eq ) activity is a radiological risk index that is used to assess the harmful effects of ionizing radiation on health. This index is associated to the internal and external exposure due to radon and its daughter progenies. The expression of Ra eq is obtained by assuming that 370 Bq ⋅ Kg −1 of 238 U , 259 Bq ⋅ Kg −1 of 232 Th and 4810 Bq ⋅ Kg −1 of 40 K produce equal gamma dose. The Radium equivalent activity is mathematically defined by [6]:

Fig. 1 Study area and sample locations
Where A U , A Th and A K are the specific activity concentration of 238 U , 232 Th and 40 K respectively. The recommended maximum limit value of Ra eq activity is 370 Bq ⋅ Kg −1 .

Absorbed dose ( D R )
The gamma dose rate is a measure of the gamma radiation energy released by radionuclides per unit time and unit mass of exposed material. It is expressed in Gray per hour ( Gy ⋅ h −1 ) and represents the level of a body exposure to an external source of radiation. The absorbed dose at one meter ( 1 m ) above the ground, due to a radionuclide, can be determined by the specific activity of this radionuclide. Theoretical calculations have allowed to establish relation between specific activity and gamma dose rate which expression is given by the following formula [6]:

Annual effective dose equivalent (AEDE)
The annual effective dose equivalent is the amount energy of gamma radiation released by a radioactive source which takes into account the radio-sensitivity of different organs of the human body. It is expressed by the following relation [6]:

Annual gonadal dose equivalent (AGDE)
The annual gonadal dose equivalent gives an estimation of the annual equivalent dose received in the gonads (reproductive organ). The AGDE for an exposition to natural source of radioactivity is given by the following formula [6]:

Radiation hazard index
External hazard index ( H ex ) The external hazard index is a parameter used to evaluate the radiological risks of materials due to an external exposure to gamma radiation. The recommended value of H ex must be less than 1 mSv ⋅ y −1 and its expression is given by the following relation [6]: The internal hazard index is a parameter used to assess the risks associated with the presence of natural radionuclides in building materials due to an internal exposure to radon ( 222 Rn ) and its daughter progenies by inhalation of alpha particles. Its expression is given by the following formula :

Excess lifetime cancer risk (ELCR)
The carcinogenic effects of gamma radiation due to ingestion, inhalation and external exposure to a radioactive sources are characterized by an estimation of cancer occurring probability during a specific lifetime. This probability is determined by the calculation of the excess lifetime cancer risk (ELCR). Its expression is given by the following relation [7]: In this expression, DL is the duration of life (70 years) and RF is the risk factor for cancer considering the stochastic effects. RF is equal to 0.05 Sv −1 for the public [7].

Spatial distribution
The activity concentration of 238 U , 232 Th et 40 K were measured by gamma-ray spectrometry method for various soil samples collected at the depth of up to 20 cm in the Northeastern area of Burkina Faso. The specific activity values and their respective standard deviations (SD) of the above natural radionuclides  . 2). This is due to the contribution of localities like Pougbonne  For each sample, the proportionality ratio between concentrations of 238 U and 232 Th (U/Th), 238 U and 40 K (U/K) was calculated. The mean value of the U/Th is clearly higher than one (01) while that of the U/K is clearly lower than one (01). This means that in most of the samples, 238 U concentration is higher than that of 232 Th but lower than that of 40 K . The activity concentration of 238 U is two times higher than that of the 232 Th at Niapsi village and the activity concentration of 40 K worth twenty (20) times that of the 238 U at Gourgnel village (see Table 2).
The results presented show that the maximum concentration of 238 U measured at Niapsi is three times higher than that of the world average. It is greater than seven times the minimum concentration measured at Tangarce. The maximum 232 Th concentration measured at Niapsi is twice that of Tangarce and the 40 K concentration at Niapsi is slightly lower than that obtained at Tangarce. In the study conducted by [2], the gamma dose rate measured at Tangarce was equal to the world average of 0.050 Sv ⋅ h −1 and its maximum value measured at Niapsi was 0.300 Sv ⋅ h −1 , which is six times that of Tangarce. This ratio is comparable to that of the uranium concentrations. We can deduce that the radiological anomaly of the Niapsi locality is an uranium one. The high dose of radioactivity measured in this locality is mainly due to the presence of 238 U . The maximum 232 Th concentration has been measured at Pougbonne and it is twice that of the minimum concentration obtained at Tangarce. The 238 U concentration of Pougbonne is worth five times that of Tangarce while the 40 K concentration of Pougbonne is slightly lower than that of Tangarce. The gamma dose rate measured at Pougbonne is three times the gamma dose rate measured at Tangarce. It is true that in this situation, the 238 U concentration of the locality is quite high compared to Tangarce but 232 Th concentration ratio is comparable to that of the gamma dose rates of these two localities. We can therefore say that the anomaly observed at Pougbonne is mainly due to 232 Th . The maximum 40 K concentration has been measured at Gourgnel and it is 2.76 times that of the concentration obtained at Tangarce. The gamma dose rate measured at Gourgnel is approximately 2.6 times that obtained at Tangarce. This similarity of proportions allows us to affirm that the anomaly observed at Gourgnel is mainly due to the presence of 40 K . The main natural radionuclides have been identified and quantified. The results show a predominance of 238 U , 232 Th and 40 K respectively in the localities of Niapsi, Pougbonne and Gourgnel. These results are important because they will make it possible not only to be able to orient research within the framework of a mining exploration activity but also to assess the radiological risks on human health. The village of Gourgnel is particularly rich in 40 K and this could reflect the soil fertility of this locality. A comparison of the results with values reported worldwide from various sample collection depths, presented in Table 3, shows that 238 U activities concentration in the study area is higher than those reported by [8][9][10][11][12][13][14][15][16] and lower than that of Kerala in India [17]. The 232 Th activity concentration is just below the value reported at Kerala and Karnataka and higher than those of the other localities. For the 40 K , the activity concentration in the study area is lower than the value of five reported activity concentrations. The activity concentration calculation of 238 U , 232 Th and 40 K from soil samples collected in the northeastern area of Burkina Faso reveals that except the 40 K , activity concentrations of above radionuclides are well above that of the world average. This confirms the existence of some radiological anomalies. To know the impact of these anomalies on people health, the radiological hazards will be assess in the following lines.

Radiological parameters
The study area is a predominantly rural locality where farming and gold mining are the most practiced activities.
Mining boom phenomenon has caused the displacement of populations to these localities and sometimes children and women leave school and agriculture activities in favor of gold mining. These vulnerable people are permanently found in the pits and in contact with extracted soil dust which could contains a significant concentrations of radionuclides. The present study has shown that the locality has a relatively high concentration of 238 U and 232 Th . Thus radiological risks were assessed by the calculation of radiological parameters which values were compared to authorized international standards and results are shown in Tables 4 and 5.
Radium equivalent Ra eq calculated using the activity concentrations of 238 U , 232 Th and 40 K of the study area varies from 62.87 Bq ⋅ Kg −1 to 189.47 Bq ⋅ Kg −1 . The mean value of Ra eq is 129.55 Bq ⋅ Kg −1 , which is much less than threshold value of 370 Bq ⋅ Kg −1 . Table 4 shows that the minimum value of the absorbed gamma dose rate ( D R ) is obtained equal to 29.5 Gy ⋅ h −1 at Tangarce, while the   Figure 3 shows that for Niapsi village, D R is a little higher than this maximum limit value. The values of annual effective dose equivalent (AEDE) ranged from 0.18 mSv ⋅ y −1 to 0.53 mSv ⋅ y −1 with a mean value of 0.37 mSv ⋅ y −1 , which is less than the recommended value of 0.8 mSv ⋅ y −1 while its value at Niapsi is higher than the threshold one as it is shown on figure 3. The annual gonadal dose equivalent (AGDE) ranged from 209.01 Sv ⋅ y −1 to 594 Sv ⋅ y −1 and the mean value in the study area is 413.36 Sv ⋅ y −1 ,   Fig. 3 shows that the ELCR level is higher than the recommended value of 2.90 × 10 −4 in Poubonne and Niapsi villages. The radiological risk assessment revealed that in the localities of Niapsi and Pougbonne, the value of the excess lifetime cancer risk (ELCR) is significantly higher than the recommended value. Ramasamy et al. [17] showed that in the locality of Kerala in India, the average ELCR value is six times that of the recommended value. In this locality there are at least 35, 000 new cases of cancer each year and between 70% and 90% of cancers are caused by environmental effects. It is therefore wise to give particular importance to this parameter and to pay attention on the naturally occurring radioactivity in Niapsi and Pougbonne.

Conclusion
The specific activity concentrations of 238 U and 232 Th obtained in the villages of Pougbonne and Niapsi are clearly higher than the world average. The 238 U concentration at Niapsi is three times higher than that of the world average, which could confirm the presence of uranium anomaly in the locality. Gourgnel is the only one locality in the study area where 40 K concentration is much higher than the world average. The assessment of radiological hazard shows that, at Niapsi and sometimes at Pougbonne, the values of gamma dose rate ( D R ), annual effective dose equivalent (AEDE) and excess lifetime cancer risk (ELCR) are higher than the recommended maximum limit values. In the study area, it would be important to perform a study to see the evolution of the radiological parameters with the depth.
Funding The authors have no relevant financial or non-financial interests to disclose. The authors have no competing interests to declare that are relevant to the content of this article. All authors certify that they have no affiliations with or involvement in any organization or entity with any financial interest or non-financial interest in the subject matter or materials discussed in this manuscript. The authors have no financial or proprietary interests in any material discussed in this article. The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

Conflict of interest
On behalf of all authors, the corresponding author states that there is no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.