Next Article in Journal
Hot Air and Microwave Combined Drying of Potato Monitored by Infrared Thermography
Next Article in Special Issue
Evaluation of the Radiological and Chemical Risk for Public Health from Flour Sample Investigation
Previous Article in Journal
An All-Batch Loss for Constructing Prediction Intervals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessment of Natural Radioactivity and Radiological Risks in River Sediments from Calabria (Southern Italy)

1
Department of Mathematics and Informatics, Physics and Earth Sciences (MIFT), University of Messina, Viale F. Stagno d’Alcontres 31, 98166 Messina, Italy
2
Istituto Nazionale di Geofisica e Vulcanologia (INGV), Sezione di Palermo, Via Ugo La Malfa 153, 90146 Palermo and Milazzo Office, Via dei Mille 46, 98057 Milazzo (ME), Italy
3
Stazione Zoologica Anton Dohrn (SZN), Villa Comunale, 80121 Napoli, Italy
4
Environmental Protection Agency of Calabria, Italy (ARPACal), Department of Reggio Calabria, Via Troncovito SNC, 89135 Reggio Calabria, Italy
5
Department of Biomedical and Dental Sciences and Morphofunctional Imaging SASTAS Section, University of Messina, Viale Annunziata, 98168 Messina, Italy
*
Author to whom correspondence should be addressed.
Submission received: 20 January 2021 / Revised: 5 February 2021 / Accepted: 11 February 2021 / Published: 15 February 2021
(This article belongs to the Special Issue Advances in Environmental Applied Physics)

Abstract

:
This study was developed to carry out a comprehensive radiological assessment of natural radioactivity for river sediment samples from Calabria, southern Italy, and to define a baseline background for the area on a radiation map. In the studied area, elevated levels of natural radionuclides are expected, due to the outcropping acidic intrusive and metamorphic rocks from which the radioactive elements derive. To identify and quantify the natural radioisotopes, ninety river sediment samples from nine selected coastal sampling points (ten samples for each point) were collected as representative of the Ionian and the Tyrrhenian coastline of Calabria. The samples were analyzed using a gamma ray spectrometer equipped with a high-purity germanium (HPGe) detector. The values of mean activity concentrations of 226Ra, 232Th and 40K measured for the studied samples are (21.3 ± 6.3) Bq kg−1, (30.3 ± 4.5) Bq kg−1 and (849 ± 79) Bq kg−1, respectively. The calculated radiological hazard indices showed average values of 63 nGy h−1 (absorbed dose rate), 0.078 mSv y−1 (effective dose outdoors), 0.111 mSv y−1 (effective dose indoors), 63 Bq kg−1 (radium equivalent), 0.35 (Hex), 0.41 (Hin), 0.50 (activity concentration index) and 458 µSv y−1 (Annual Gonadal Equivalent Dose, AGED). In order to delineate the spatial distribution of natural radionuclides on the radiological map and to identify the areas with low, medium and high radioactivity values, the Surfer 10 software was employed. Finally, the multivariate statistical analysis was performed to deduce the interdependency and any existing relationships between the radiological indices and the concentrations of the radionuclides. The results of this study, also compared with values of other locations of the Italian Peninsula characterized by similar local geological conditions, can be used as a baseline for future investigations about radioactivity background in the investigated area.

1. Introduction

Radioactivity naturally occurs in the environment. Aside from the cosmic radiation which bombards the earth from the outside of the atmosphere, the predominant natural source of human exposure to radiation doses is due to the terrestrial radiation which includes radionuclides of 238U, 232Th, 235U series and 40K [1].
It is established that the abundance of these primordial radionuclides depends on the local environment’s geology [2,3]. Moreover, natural and anthropogenic events such as flooding, radionuclide uptake by plants and the use of fertilizers, usually provoke an alteration in the distribution of the abundance levels of the radionuclides in the environment [4,5,6].
Studies of natural radiation are needed to establish reference levels, especially in areas where the hazard of radioactive exposure may be higher due to the lithological features of the geological context.
In this framework, central and southern Calabria (Southern Italy) can be considered as “risky” areas, since granitoids and metamorphic rocks that represent sources of radionuclides widely outcrop along the Serre and the Aspromonte Massifs. Moreover, recent studies highlighted the presence of heavy minerals, such as monazite and zircons, in the Calabrian rocks [7], which might contain high amounts of radioactive U and Th.
In this study, fluvial deposits from nine different selected sites from the Ionian and Tyrrhenian coastline of Calabria were analyzed. River sediments, consisting of mineral particles with different size and chemical composition, are considered reliable long-term indicators of river pollution by radionuclides, because water pollution components are deposited in the sediments [8]. Major sources of natural radionuclides in sediments have different possible origins, like weathering and recycling by erosion of terrestrial minerals and rocks (igneous or metamorphic) containing 40K and radionuclides of uranium and thorium radioactive series, and also rainfall and other depositional phenomena such as landslides and precipitations [9]. The main aim of this paper is to identify and quantify natural gamma-emitting radionuclides and to evaluate the radiological hazard effects [10]. For this purpose, High-Purity Germanium (HPGe) gamma spectrometry was used to measure the 226Ra (in secular equilibrium with 238U), 232Th and 40K activity concentration [11]. Statistical and radiological maps were produced in order to detect the relationship among values and to identify the most hazardous zone.
Moreover, since sandy natural materials were widely used in the field of civil construction, a series of radiological indices were calculated to evaluate the potential hazard connected with the building use of these materials.

2. Geological Setting

The investigated river sediments come from the Ionian and Tyrrhenian seaboard of the Calabrian coastline (Figure 1), and they derive from the dismantling of the Serre Massifs and Aspromonte due to erosion processes. These Massifs belong to the Calabria-Peloritani Orogen (CPO), which is the peri-Mediterranean orogenic Alpine nappe system that comprises the whole Calabria and the north-eastern sector of Sicily [12]. It includes the Sila and Catena Costiera Massifs in northern Calabria, the Serre and Aspromonte Massifs in central and southern Calabria and the Peloritani Mountains in Sicily [12].
The Serre Massif (Central Calabria) is characterized by outcrops of granitoids and Variscan metamorphic rocks. A thin early Mesozoic sedimentary cover is present only along the Ionian margin. According to recent studies [12,13,14], the Serre Massif represents a nearly complete continental crust section, which is divided into lower, middle and upper crustal portions. The lower one is composed of granulitic and migmatitic complexes, whereas the middle one is made up by late-Variscan acidic granitoid rocks, such as tonalites, granodiorites and granites. Conversely, the upper crustal portion consists of high- to low-grade metapelites (Mammola and Stilo-Pazzano complexes) [13,14].
The Aspromonte Massif (Calabria), together with the Peloritani Mountains (Sicily), constitutes the southern sector of the CPO (Calabria Peloritani Orogen). This Massif was generated by the juxtaposition of three different tectonic units. The Madonna di Polsi Unit is the lowermost unit and it is composed of greenschist to amphibolite facies metamorphic rocks showing only Alpine metamorphism. The higher-most unit is represented by the Stilo Unit, made up of phyllites, schists and paragneisses derived from the variscan metamorphism of Paleozoic sediments. The Aspromonte Unit, which is sandwiched between the Stilo and Madonna di Polsi Unit, exhibits high-grade metamorphites and late-variscan granitoids [12,13,14].
Paleozoic and early Mesozoic rocks of both Serre and Aspromonte Massif are covered by sedimentary deposits ranging in age from the Late Cretaceous to the Pleistocene.

3. Materials and Methods

3.1. Samples Collection and Preparation

Ninety samples of river sediments, around 1 kg each, were collected in nine selected sites of the Ionian and Tyrrhenian coast of Calabria for subsequent laboratory gamma spectrometric measurements. Locations, coordinates and labels of the collected samples are indicated in Figure 1 and in Table 1.
Samples were collected with a metal sampler and stored in labeled plastic boxes. Adequate care was taken to prevent their contamination, particularly during transportation to the laboratory. Samples were oven-dried for 24 h at 110 °C until moisture was completely removed and constant mass was attained. They were grounded and further homogenized by passing them through a 2 mm sieve using a sieve shaker. After, they were successfully inserted in Marinelli hermetically sealed containers of 1 L capacity. After 40 days, the secular radioactive equilibrium between 226Ra and their daughter products was attained and samples were ready for gamma spectrometry counting.

3.2. Gamma Spectrometry Measurements

Specific activity measurements of the investigated samples were performed by using a high-purity germanium (HPGe) detector (Ortec, Oak Ridge, TN, USA) for gamma spectrometry analysis.
In order to reduce statistical uncertainty, samples were counted for 70,000 seconds and spectra were analyzed in order to obtain the activity concentration of 226Ra, 232Th and 40K. The 226Ra-specific activity was calculated through a weighted average of 214Pb and 214Bi (in secular equilibrium) specific activities. The 232Th activity concentration was determined by using the 911.21 and 968.97 keV 228Ac γ-ray lines; for 40K, the evaluation was performed from its γ-line at 1460.8 keV. The experimental setup was composed by a positive biased Ortec HPGe detector (GEM) (Full Width at Half Maximum, FWHM) of 1:85 keV, peak to Compton ratio of 64:1, relative efficiency of 40% at the 1:33 MeV 60Co-line), placed inside lead wells to shield the background radiation environment. Efficiency and energy calibrations were performed using a multipeak Marinelli geometry gamma source (AK-5901) of 1 L capacity, covering the energy range, 59.54–1836 keV, customized to reproduce the exact geometries of samples in a water-equivalent epoxy resin matrix [15].
The Gamma Vision (Ortec) software was used for data acquisition and analysis [16].
The activity concentration (Bq kg−1) of the investigated radionuclides was calculated using the following formula [17]:
C = N E ε E t γ d M
where NE indicates the net area of a peak at energy E, εE and γd are the efficiency and yield of the photopeak at energy E respectively, M is the mass of the sample (kg) and t is the live time (s).
The measurement uncertainty is a combined standard one at coverage factor k = 2, taking into account the following components: uncertainty counting statistics, uncertainty in nuclear data library, uncertainty due to calibration efficiency and uncertainty about the quantity of the sample [18].
The quality of the gamma spectrometry experimental results was certified by the Italian Accreditation Body (ACCREDIA) [19]. This implies the continued verification (with annual periodicity) of the maintenance of the gamma spectrometry method performance characteristics.
In particular, with regard to the blanks, they are generally determined by acquiring an empty lead well spectrum or an ultrapure water sample in the sample holder of the desired geometry (for determination of natural-emitting gamma radionuclides) or, in the case of samples requiring a material of support (e.g., atmospheric particulate aspirated on filter), by acquiring the support of clean sampling, under the same geometric conditions as the sample. The measurement of the blanks is performed on a quarterly basis.
Regarding the precision and accuracy assessment, we proceed as follows:
For the precision, with the double test method, the repeatability of the method is verified over time. A certified reference material (also containing the radionuclides K-40, Ra-226 and Th-232) is analyzed twice, with the activity concentration of the radionuclide of interest defined as X1 (first measurement) and X2 (second measurement). The probability level p = 0.95 is considered. The following formula is applied:
| x 1 x 2 | 2 · s r · t
where t is the Student variable and sr is the standard deviation of repeatability obtained in the validation phase.
For the accuracy, it is evaluated by comparing the reference (certified sample) and the measured values of the radionuclide of interest, taking into account the uncertainty. This is done using the u-test, with the following acceptability criterion:
u t e s t = measured reference u m e a s 2 u r e f 2 2
where umeas2 and uref2 are the uncertainties of the measured and reference values, respectively.

3.3. Radiological Indices

3.3.1. Absorbed Gamma Dose Rate

The absorbed gamma dose rate calculation is the first major step to evaluate the health risk [20]. It was estimated as follows [21]:
D (nGy·h−1) = 0.462·CRa + 0.604·CTh + 0.0417·CK
where CRa, CTh and CK are the mean activity concentrations (Bq kg−1) of 226Ra, 232Th and 40K in the river sediments, respectively.

3.3.2. The Annual Effective Dose Equivalent

Annual estimated average effective dose equivalent (AEDE) received by an individual was calculated using the following equations with an outdoor occupancy of 20% and 80% for indoors [22]:
AEDE (outdoor) (mSv·y−1) = absorbed dose (nGy·h−1) × 8760·h × 0.7·Sv·Gy−1 × 0.2 × 10−6
AEDE (indoor) (mSv·y−1) = absorbed dose (nGy·h−1) × 8760·h × 0.7·Sv·Gy−1 × 0.8 × 10−6

3.3.3. Radium Equivalent Activity

The activity concentrations of the radionuclides measured in the analyzed samples do not give a direct measure of the γ-radiation dose incurred by the population. To estimate the gamma radiation dose due to activity concentrations of 226Ra, 232Th and 40K in the river samples, especially when used as components of building construction, the radium equivalent activity is regularly utilized.
It is an index which describes the activities of 226Ra, 232Th and 40Kin a single activity term. The gamma energy doses from these radionuclides in river sediments are different even if they are present in the same amount. Radium equivalent activity (Raeq) estimated in Bq kg−1 is evaluated with conditions that 1 Bq kg−1 of 226Ra or 1.43 Bq kg−1 of 232Th or 0.077 Bq kg−1 of 40K produce equal gamma dose rate.
The radium equivalent activity was calculated as [23]:
Raeq (Bq·kg−1) = CRa + 1.43 CTh + 0.077·CK

3.3.4. Hazard Indices

External and internal radiation hazard indices were defined to limit the radiation dose to 1 mSv y−1 [24].
The external hazard index (Hex) was calculated using the given equation:
Hex = (CRa/370 + CTh/259 + CK/4810) ≤ 1
It must be lower than the unity for the radiation hazard to be negligible.
The internal hazard index (Hin) gives the internal exposure to carcinogenic radon and its short-lived progeny and it is given by the following formula:
Hin = (CRa/185 + CTh/259 + CK/4810) ≤ 1

3.3.5. Activity Concentration Index

Another radiation hazard, called the activity concentration index (I), has been defined by the European Commission [25] and it is given below:
I = (CRa/300 + CTh/200 + CK/3000)
The activity concentration index is correlated with the annual dose rate due to the excess external gamma radiation caused by superficial material. Values of I ≤ 1 correspond to a criterion of 1 mSv y−1 [25]. Thus, this index should be used only as a screening tool to identify materials that might be of concern to be used as construction materials [26].

3.3.6. Annual Gonadal Equivalent Dose (AGED)

Some organs exhibit high sensitivity to radiation as they are highly susceptible to harmful effects of radiation. The gonads, the active bone marrow and the bone surface cells are considered to be the organs of importance due to their high radiosensitivity. The annual gonadal equivalent dose (AGED) due to the specific activities of 226Ra, 232Th and 40K is evaluated using the following formula [27]:
AGED (µSv·y−1) = 3.09·CRa + 4.18·CTh+ 0.314·CK
where CRa, CTh and CK are the mean activity concentrations of 226Ra, 232Th and 40K, respectively.

3.4. Delineation of Radiological Map

Surfer 10 software was employed to generate a radiological map for the study area [28]. The software uses the kriging interpolation method, used in geostatistics to predict the value of a random variable over a spatial region. Using the values of radiological measurements at all the sampling sites with different geographic coordinates, kriging allows the prediction of the radiological values throughout the regions where there are no experimental observations in the entire study area.
This method is widely employed in hydrogeology and environmental science, among others.

3.5. Statistical Analysis

The XLSTAT statistical software for Windows was used for all statistical calculations [29].
With the aim to individuate the presence of the relationships among the original variables (Pearson correlation analysis), an exploratory method (Principal Component Analysis, PCA) was performed. The PCA elaboration ensures the reduction of the data dimensionality, whereas the combinations of variables identified by the PCs provide the greatest contribution to sample variability. However, before the elaboration, the logarithmic transformation of the dataset was conducted and the appropriateness of the dataset was preliminarily checked by the Kaiser–Meyer–Olkin (KMO) test and the Bartlett test (sphericity test), as already reported in our previous work [30].

4. Results and Discussion

4.1. 226. Ra, 232Th and 40K Activity Concentration in the River Sediments, Spatial Distribution and Classification of Hazard Regions

The activity concentrations of the three naturally occurring 226Ra, 232Th and 40K radionuclides in river sediments across the study area were determined for the sampling points reported in Table 1.
The obtained results are reported in Table 2. As shown in the Table, the highest average concentration, (31.3 ± 8.3) Bq kg−1, of 226Ra was obtained for samples from the site ID8 (Siderno–Novito river), whereas the lowest mean value, (14.5 ± 4.9) Bq kg−1, was detected for sediments from the site ID1 (Gioia Tauro–Budello river). Concerning 232Th, samples from the site ID8 are characterized by the highest average concentration, (43.5 ± 4.9) Bq kg−1, while those from the site ID2 (Gioia Tauro–Petrace river) have the lowest mean value, (20.3 ± 2.9) Bq kg−1. Sediments from the site ID1 show the highest average concentration of 40K, with (1088 ± 98) Bq kg−1, while the lowest mean concentration of (643 ± 53) Bq kg−1 was obtained for samples from the site ID 6 (Africo–Laverde river).
Across all sampling locations, the activity concentration values of 226Ra ranged between 10.2 and 33.9 Bq kg−1, with overall mean value of 21.3 ± 6.3 Bq kg−1, that of 232Th varied from 16.5 to 48.6 Bq kg−1 and the overall average value is (30.3 ± 4.5) Bq kg−1. The 40K ranged between 528 and 1450 Bq kg−1, with an overall average value of (849 ± 79) Bq kg−1.
The average values of 226Ra, 232Th and 40K obtained in this study, for each sampling site, were compared with the world average value reported in Reference [1]. All the average values of 226Ra are lower than that of the world average (35 Bq kg−1). For 232Th, the mean values are lower than that of the world average (30 Bq kg−1) except for samples from the sites ID3 (39.5 ± 4.8 Bq kg−1), ID6 (32.5 ± 9.2 Bq kg−1), ID7 (34.5 ± 5.6 Bq kg−1) and ID8 (43.5 ± 4.9 Bq kg−1). The average concentration of 40K is higher than the 420 Bq/kg average concentration reported in Reference [1] for all the analyzed samples.
Contour maps, showing the spatial distributions of 226Ra, 232Th and 40K, are presented in Figure 2.
As can be seen in the 226Ra map (Figure 2a), there are four distinct regions. The activity concentration of 226Ra in these regions was classified as low (CRa < 15 Bq kg−1), moderate (15 Bq kg−1 < CRa < 20 Bq kg−1), high (20 Bq kg−1 < CRa < 25 Bq kg−1) and very high (CRa > 25 Bq kg−1), which are deep blue, light blue, green and yellow (and red) areas on the map. In the case of 232Th, two regions of very low and moderately low activity concentrations were identified. As seen from Figure 2b, the southern part of the study area as well as the mid-west and the northern part exhibit low activity concentrations of 232Th with typical activities concentrations less than 33 Bq kg−1. Otherwise, the middle-east zone of the study area exhibits high activity, higher than 35 Bq kg−1. The delineated map of 40K (Figure 2c) shows two distinct zones of concentration. The blue zone in the southern part belongs to the low-concentration region (lower than 700 Bq kg−1). Whereas, other areas are classified as high-concentration zones (higher than 700 Bq kg−1).

4.2. Dose Assessment and Hazard Indices

Table 3 reports the calculated average values of adsorbed and effective doses together with the hazard indices. The absorbed gamma dose rate was estimated using Equation (4). Measured values range from 52 to 79 nGy h−1, with an average value of 63 nGy h−1. This average value is near the world mean value of 60 nGy/h [1], although 66% of the analyzed samples showed values that are in excess than the world mean. The variability of the adsorbed and effective doses is attributable to the different lithologic components of the considered areas.
Equation (5) was used to evaluate the annual outdoor effective dose due to the activities of 226Ra, 232Th and 40K in the analyzed samples. The obtained values range from 0.064 to 0.097 mSv y−1, with an average value of 0.078 mSv y−1. In a similar way, the annual indoor effective dose (Equation (6)) ranges from 0.070 to 0.323 mSv y−1, with a mean value of 0.111 mSv y−1. The total average annual effective dose is therefore of 0.189 mSv y−1. This value is lower than the worldwide average of annual effective dose (0.48 mSv y−1 [1]). The mean value is also lower than 1 mSv y−1, which is set as the maximum limit by Reference [31].
The radium equivalent activity (Raeq) was calculated using Equation (7) in order to determine the suitability of the samples when they are used as components in building materials. The obtained Raeq values vary from 52 to 79 Bq kg−1 with an average value of 63 Bq kg−1, which is lower than 370 Bq kg−1 set as the upper limit for building materials [29]. This indicates that the investigated samples may not be hazardous if used in the field of civil construction. The minimum value of radium equivalent activity was obtained for the river sediments come from the site ID4 (Reggio Cal.–Calopinace river), while the maximum one characterizes the samples collected at the site ID8 (Siderno–Novito river).
The external and internal hazard indices, given by Equations (8) and (9), are lower than the unity for all investigated samples, thus the radiological risks can be considered negligible.
The activity concentration index, calculated with Equation (10), is again lower than the unity for all analyzed samples. Thus, these river sediments are not be of concern to be used as construction materials.
The obtained values of annual gonadal equivalent dose (AGED), calculated by Equation (11), range from 378 (site ID4) to 566 μSv y−1 (site ID8), with an average of 458 μSv y−1, for the investigated samples. These values are higher than the world average value for AGED, for sediments, which in the literature was found to be 300 μSv y−1 [32].

4.3. Statistical Features

The radiological parameters obtained from the measured activity concentrations of 40K, 226Ra and 232Th were subjected to Pearson’s correlation analysis, with the aim to deduce the interdependency and identify any existing relationships between the radiological indices and the concentrations of the radionuclides. This will help in making valid assessments concerning the nature and distribution of these radionuclides in the study area. The results of the Pearson’s correlation coefficients are presented in Table 4. Among the negative correlations, the most interesting are those between 40K and 226Ra (−0.287) and between 40K and 232Th (−0.305), while a perfect positive correlation is found between absorbed dose rate and effective dose outdoors (1.00), and between Raeq and absorbed dose rate (1.00) and effective dose outdoors (1.00). The correlations of the three radionuclides are positive and in agreement with all the radiological parameters considered in this study at p < 0.01, except for correlations between 226Ra and Hin (0.739) and between 232Th and Hin (0.719), obtained at p < 0.05.
Moreover, 226Ra and 232Th are negatively correlated with the annual effective dose equivalent indoors (−0.378 and −0.281, respectively). The radiological parameters are positively and significantly correlated with each other, showing direct relationships.
In order to better define the relationships among the specific activities of 226Ra, 232Th and 40K, the PCA multivariate statistical method was carried out. The geological feature of sediments is the most important variable, neglecting the others (erosive power of rivers, solid transport, selective erosion, weathering, etc.), which is useful to evaluate the degree of relationship between the isotopic specific activities and the samples’ geolocation. Results of the PCA allow to group the sampling sites in four clusters that seem to be comparable to each other (Figure 3) in terms of activity concentration of 226Ra, 232Th and 40K. Specifically, sites ID1 and ID2 show positive correlations with 40K-specific activity but negative linear relationships with 226Ra and 232Th ones, whereas sites IDs 3, 4 and 5 are negatively correlated to the activity concentrations of 226Ra, 232Th and 40K. The sites ID6 and ID7 show instead a positive correlation with values of 226Ra- and 232Th-specific activities, and exhibit a negative correlation with 40K. Finally, sites ID8 and ID9 show positive correlations of the three investigated radionuclides’ activity concentrations. We believe that the differences of 226Ra-, 232Th- and 40K-specific activities among the identified clusters mainly depend on the compositional/mineralogical features of rocks from which the analyzed river sediments derive. The sediments collected along the rivers located in the Tyrrhenian side (ID1 and ID2) record the chemical fingerprint of the granitoid rocks from the Serre Massif. On the contrary, the Ionian river sediments (ID6, ID7, ID8 and ID9) reflect the chemical and mineralogical composition of the outcropping Late Cretaceous and Cenozoic terrigeous sediments. Moreover, the ID6 and ID7 rivers are characterized by a component provided by the erosion of the Aspromonte Unit, whereas clasts from the upper crustal complexes (Stilo-Pazzano and Mammola) and from granitoids supply the ID8 and ID9 rivers. Regarding the rivers on the western and southern margin of the Aspromonte Massif (ID3, ID4 and ID5), the main lithologies occurring within the drainage basins come from the Aspromonte and the Stilo Units.

5. Conclusions

The activity concentration and the distribution of natural terrestrial radionuclides 226Ra, 232Th and 40K for river sediments from selected locations of the Calabrian coastline (Southern Italy) were measured using HPGe gamma spectrometry. The obtained average values of the absorbed dose rate in air, AEDE (outdoor and indoor), Raeq, Hex and Hin, I and AGDE are comparable with average worldwide ranges and they fell within the maximum recommended values.
The background radiation maps for the study area evidence that the Tyrrhenian side shows, in general, the lowest activity concentrations of 226Ra and 232Th, whereas the Ionian zone is characterized by the highest ones. Concerning the 40K, an increasing of the specific activity was defined from south to north. Multivariate statistical analysis allowed to assess the nature and distribution of natural radionuclides in the investigated area and highlighted the relationships between chemical/mineralogical features of the local outcropping lithologies.
Data provided in this study will serve to contribute as a reference database for future researches in the area, and also as a pilot scheme to provide a comprehensive background radiation map for Italy, which is not yet available.

Author Contributions

Conceptualization, F.C.; methodology, F.C. and A.F.M.; validation, F.C. and G.S.; formal analysis, F.C.; investigation, F.C., A.F.M., G.S., M.D.B., M.R.F., D.R. and F.I.; resources, G.S. and G.B.; data curation, F.C. and A.F.M.; writing—original draft preparation, F.C.; supervision, G.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. United Nations Scientific Committee on the Effects of Atomic Radiation (UNSCEAR). Sources and Effects of Ionizing Radiation. Report of the United Nations Scientific Committee on the Effects of Atomic Radiation to the General Assembly; United Nations: New York, NY, USA, 2000. [Google Scholar]
  2. Baeza, A.; Del Rio, M.; Jimenez, A.; Miro, C.; Paniagua, J. Influence of geology and soil particle size on the surface-area/volume activity ratio for natural radionuclides. J. Radioanal. Nucl. Chem. 1995, 189, 289–299. [Google Scholar] [CrossRef]
  3. Tzortzis, M.; Tsertos, H.; Christofides, S.; Christodoulides, G. Gamma-ray measurements of naturally occurring radioactive samples from Cyprus characteristic geological rocks. Radiat. Meas. 2003, 37, 221–229. [Google Scholar] [CrossRef] [Green Version]
  4. Navas, A.; Soto, J.; Machin, J. Edaphic and physiographic factors affecting the distribution of natural gamma-emitting radionuclides in the soils of the Arnas catchment in the Central Spanish Pyrenees. Eur. J. Soil Sci. 2002, 53, 629–638. [Google Scholar] [CrossRef] [Green Version]
  5. Pulhani, V.A.; Dafauti, S.; Heqde, A.G.; Sharma, R.M.; Mishra, U.C. Uptake and distribution of natural radioactivity in wheat plants from soil. J. Environ. Radioact. 2005, 79, 331–346. [Google Scholar] [CrossRef]
  6. Chandrajith, R.; Seneviratna, S.; Wickramaarachchi, K.; Attanayake, T.; Aturaliya, T.N.C.; Dissanayake, C.B. Natural radionuclides and trace elements in rice field soils in relation to fertilizer application: Study of a chronic kidney disease area in Sri Lanka. Environ. Earth Sci. 2010, 60, 193–201. [Google Scholar] [CrossRef]
  7. Fornelli, A.; Piccarreta, G.; Del Moro, A.; Acquafredda, P. Multi-stage melting in the lower crust of the Serre (southern Italy). J. Petr. 2002, 43, 2191. [Google Scholar] [CrossRef]
  8. Sabatino, G.; Di Bella, M.; Caridi, F.; Italiano, F.; Romano, D.; Magazù, S.; Gnisci, A.; Faggio, G.; Messina, G.; Santangelo, S.; et al. Radiological assessment, mineralogy and geochemistry of the heavy-mineral placers from the Calabrian coast (South Italy). J. Instrum. 2019, 14, P05015. [Google Scholar] [CrossRef]
  9. El-Gamal, A.; Nasr, S.; El-Taher, A. Study of the spatial distribution of natural radioactivity in the upper Egypt Nile river sediments. Radiat. Meas. 2007, 42, 457–465. [Google Scholar] [CrossRef]
  10. Caridi, F.; Marguccio, S.; Belvedere, A.; D’Agostino, M.; Belmusto, G. A methodological approach to a radioactive sample analysis with low-level γ-ray spectrometry. J. Instrum. 2018, 13, P09022. [Google Scholar] [CrossRef]
  11. Caridi, F.; Marguccio, S.; Belvedere, A.; D’Agostino, M.; Belmusto, G.; Durante, G.; Trozzo, R.; Fullone, F. Natural radioactivity measurements and dosimetric evaluations in soil samples with a high content of NORM. Eur. Phys. J. Plus 2017, 132, 56. [Google Scholar] [CrossRef]
  12. Cirrincione, R.; Fazio, E.; Fiannacca, P.; Ortolano, G.; Pezzino, A.; Punturo, R. The Calabria-Peloritani Orogen, a composite terrane in Central Mediterranean; Its overall architecture and geodynamic significance for a pre-Alpine scenario around the Tethyan basin. Periodico di Mineralogia 2015, 84, 701–749. [Google Scholar] [CrossRef]
  13. Caggianelli, A.; Liotta, D.; Prosser, G.; Ranalli, G. Pressure-temperature evolution of the late Hercyian Calabria continental crust: Compatibility with post-collisional extensional tectonics. Terra Nova 2007, 19, 502. [Google Scholar] [CrossRef]
  14. Angì, G.; Cirrincione, R.; Fazio, E.; Fiannacca, P.; Ortolano, G.; Pezzino, A. Metamorphic evolution of preserved Hercynian crustal section in the Serre massif (Calabria-Peloritani Orogen, southern Italy). Lithos 2010, 115, 237. [Google Scholar] [CrossRef]
  15. Caridi, F.; D’Agostino, M.; Belvedere, A.; Marguccio, S.; Belmusto, G.; Gatto, M.F. Diagnostics techniques and dosimetric evaluations for environmental radioactivity investigations. J. Instrum. 2016, 11, C10012. [Google Scholar] [CrossRef]
  16. Caridi, F.; D’Agostino, M.; Messina, M.; Belvedere, A.; Marguccio, S.; Belmusto, G.; Marcianò, G.; Grioli, L. Lichens as environmental risk detectors. Eur. Phys. J. Plus 2017, 132, 189. [Google Scholar] [CrossRef]
  17. Caridi, F.; Belvedere, A.; D’Agostino, M.; Marguccio, S.; Marino, G.; Messina, M.; Belmusto, G. An investigation on airborne particulate radioactivity, heavy metals and polycyclic aromatic hydrocarbons composition in Calabrian selected sites, southern Italy. Ind. J. Environ. Protect. 2019, 39, 321–326. [Google Scholar]
  18. Caridi, F.; Messina, M.; Belvedere, A.; D’Agostino, M.; Marguccio, S.; Settineri, L.; Belmusto, G. Food salt characterization in terms of radioactivity and metals contamination. Appl. Sci. 2019, 9, 2882. [Google Scholar] [CrossRef] [Green Version]
  19. Caridi, F.; D’Agostino, M.; Belvedere, A. Radioactivity in Calabrian (Southern Italy) wild boar meat. Appl. Sci. 2020, 10, 3580. [Google Scholar] [CrossRef]
  20. Ramasamy, V.; Suresh, G.; Meenakshisundaram, V.; Ponnusamy, V. Horizontal and vertical characterization of radionuclides and minerals in river sediments. Appl. Rad. Isot. 2011, 69, 184–195. [Google Scholar] [CrossRef] [PubMed]
  21. Caridi, F.; Pappaterra, D.; Belmusto, G.; Messina, M.; Belvedere, A.; D’Agostino, M.; Settineri, L. Radioactivity and heavy metals concentration in Italian (Calabrian) DOC wines. Appl. Sci. 2019, 9, 4584. [Google Scholar] [CrossRef] [Green Version]
  22. Caridi, F.; Messina, M.; Faggio, G.; Santangelo, S.; Messina, G.; Belmusto, G. Radioactivity, radiological risk and metal pollution assessment in marine sediments from Calabrian selected areas, southern Italy. Eur. Phys. J. Plus 2018, 133, 65. [Google Scholar] [CrossRef]
  23. Torrisi, L.; Visco, A.M.; Campo, N.; Caridi, F. Pulsed laser treatments of polyethylene films. Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. Atoms 2010, 268, 3117–3121. [Google Scholar] [CrossRef]
  24. Beretka, J.; Matthew, P.J. Natural radioactivity of Australian building materials, waste and by-products. Health Phys. 1985, 48, 87–95. [Google Scholar] [CrossRef]
  25. Official Journal of the European Union. Council Directive 2013/59/Euratom; Official Journal of the European Union: Luxembourg, 2015. [Google Scholar]
  26. Ravisankar, R.; Vanasundari, K.; Chandrasekaran, A.; Rajalakshmi, A.; Suganya, M.; Vijayagopal, P.; Meenakshisundaram, V. Measurement of natural radioactivity in building materials of Namakkal, Tamil Nadu, India using gamma-ray spectrometry. Appl. Rad. And Isot. 2012, 70, 699–704. [Google Scholar] [CrossRef] [PubMed]
  27. Darwish, D.A.E.; Abul-Nasr, K.T.M.; El-Khayatt, A.M. The assessment of natural radioactivity and its associated radiological hazards and dose parameters in granite samples from South Sinai, Egypt. J. Radiat. Res. Appl. Sci. 2015, 8, 17–25. [Google Scholar] [CrossRef] [Green Version]
  28. Surfer 10 Software. Available online: https://www.goldensoftware.com/products/surfer (accessed on 17 November 2020).
  29. XLSTAT: Statistical Software & Data Analysis Add-On for Excel. Available online: https://www.xlstat.com/en/ (accessed on 16 November 2020).
  30. Mottese, A.F.; Sabatino, G.; Di Bella, M.; Fede, M.R.; Caridi, F.; Parisi, F.; Marcianò, G.; Caccamo, M.T.; Italiano, F.; Yuce, G.; et al. Environmental screening for the assessment of potentially toxic elements content in PGI soils from the Mediterranean region (Italy and Turkey). Environ. Earth Sci. 2020, 79, 499. [Google Scholar] [CrossRef]
  31. International Commission on Radiological Protection (ICRP). 1990 Recommendations of the International Commission on Radiological Protection. ICRP Publication 60. Ann. ICRP 1991, 21, 1–3. [Google Scholar]
  32. Chandrasekaran, A.; Ravisankar, R.; Senthilkumar, G.; Thillaivelavan, K.; Dhinakaran, B.; Vijayagpal, P.; Bramha, S.N.; Venkatraman, B. Spatial distribution and lifetime cancer risk due to gamma radioactivity in Yelagiri Hills, Tamilnadu, India. Egypt J. Basic Appl. Sci. 2014, 1, 38–48. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Geological sketch map of central-southern Calabria; yellow stars: sampling sites.
Figure 1. Geological sketch map of central-southern Calabria; yellow stars: sampling sites.
Applsci 11 01729 g001
Figure 2. Contour maps showing the spatial distributions of (a) 226Ra, (b) 232Th and (c) 40K in the investigated region.
Figure 2. Contour maps showing the spatial distributions of (a) 226Ra, (b) 232Th and (c) 40K in the investigated region.
Applsci 11 01729 g002
Figure 3. Two-dimensional (2D) plots of the first two Principal Components (PCs) obtained through Principal Component Analysis (PCA) elaboration starting from 40K, 232Th and 226Ra activity concentrations (a). Mean concentration three-dimensional (3D) map of 40K, 232Th and 226Ra into the investigated sites (b).
Figure 3. Two-dimensional (2D) plots of the first two Principal Components (PCs) obtained through Principal Component Analysis (PCA) elaboration starting from 40K, 232Th and 226Ra activity concentrations (a). Mean concentration three-dimensional (3D) map of 40K, 232Th and 226Ra into the investigated sites (b).
Applsci 11 01729 g003
Table 1. Labels and Global Positioning System (GPS) locations of the sampling sites.
Table 1. Labels and Global Positioning System (GPS) locations of the sampling sites.
Site IDSampling SiteGPS Position
LatitudeLongitude
1Gioia Tauro–Budello river38.43 15.9075
2Gioia Tauro–Petrace river38.420277815.8839
3Reggio Cal.–Gallico river38.173055615.6506
4Reggio Cal.–Calopinace river38.100555615.6897
5Bova Marina–Amendolea river37.92833315.8878
6Africo–Laverde river38.929166716.8338
7Bovalivo–Bonamico river38.127516.1591
8Siderno–Novito river38.250555616.2794
9Caulonia–Allaro river38.348055616.4727
Table 2. Activity concentrations of 226Ra, 232Th and 40K in the investigated sampling sites.
Table 2. Activity concentrations of 226Ra, 232Th and 40K in the investigated sampling sites.
Scheme 226.No. of Samples226Ra
(Bq kg1)
232Th
(Bq kg1)
40K
(Bq kg1)
RangeMeanRangeMeanRangeMean
11010.2–23.614.5 ± 4.916.5–48.622.7 ± 3.2961–14501088 ± 98
21013.6–22.918.2 ± 4.417.5–23.820.3 ± 2.9949–10281009 ± 88
31015.1–19.616.2 ± 4.337.2–41.539.5 ± 4.8712–918815 ± 77
41015.1–18.816.5 ± 6.822.5–26.324.1 ± 2.9591–811722 ± 84
51019.1–21.720.6 ± 2.320.3–46.827.8 ± 3.9678–762702 ± 66
61019.3–29.426.2 ± 7.525.1–36.532.5 ± 9.2528–691643 ± 53
71023.6–31.826.5 ± 9.528.2–43.534.5 ± 5.6608–749715 ± 60
81026.6–33.931.3 ± 8.338.2–48.643.5 ± 4.9871–942917 ± 86
91017.8–26.821.9 ± 8.525.3–33.827.7 ± 3.41002–10791030 ± 98
Overall9010.2–33.921.3 ± 6.316.5–48.630.3 ± 4.5528–1450849 ± 79
Table 3. Radiological indices in the investigated sampling sites.
Table 3. Radiological indices in the investigated sampling sites.
Scheme 1.Absorbed Dose Rate (nGyh1)Effective Dose Outdoors (mSv y1)Effective Dose Indoors (mSv y1)Raeq
(Bq kg−1)
HexHinIAnnual Gonadal Equivalent Dose (AGED) (µSvy1)
1660.0810.323660.350.390.52481
2630.0770.084630.330.390.50457
3660.0800.088650.360.410.52471
4520.0640.070520.280.330.42378
5560.0680.075550.300.360.44400
6580.0720.079580.330.400.46418
7630.0770.085630.350.430.50450
8790.0970.106790.440.530.63566
9700.0860.094690.380.440.55506
Mean630.0780.111630.350.410.50458
Table 4. Pearson correlation matrix among the considered variables.
Table 4. Pearson correlation matrix among the considered variables.
Pearson Correlation Matrix (n).
Variables226Ra232Th40KAbsorbed Dose Rate (nGy h−1)Effective Dose Outdoors (mSv y−1)Effective Dose Indoors (mSv y−1)Raeq (Bq kg−1)HexHinGamma index IAGED µSv y−1
226Ra10.649−0.2870.4680.465−0.3780.4680.5590.7390.4640.417
232Th0.6491−0.3050.5460.536−0.2810.5460.6660.7190.5530.494
40K−0.287−0.30510.5950.6030.6070.5960.4650.3010.5910.647
Absorbed Dose Rate (nGy/h)0.4680.5460.59511.0000.2361.0000.9880.9360.9990.998
Effective Dose outdoors (mSv/y)0.4650.5360.6031.00010.2451.0000.9860.9330.9980.998
Effective Dose indoors (mSv/y)−0.378−0.2810.6070.2360.24510.2370.1480.0020.2230.275
Raeq (Bq/kg)0.4680.5460.5961.0001.0000.23710.9880.9360.9990.998
Hex0.5590.6660.4650.9880.9860.1480.98810.9700.9880.976
Hin0.7390.7190.3010.9360.9330.0020.9360.97010.9350.912
Gamma index I0.4640.5530.5910.9990.9980.2230.9990.9880.93510.996
AGDE (µSv/y)0.4170.4940.6470.9980.9980.2750.9980.9760.9120.9961
Bold values refer to p < 0.05.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Caridi, F.; Di Bella, M.; Sabatino, G.; Belmusto, G.; Fede, M.R.; Romano, D.; Italiano, F.; Mottese, A.F. Assessment of Natural Radioactivity and Radiological Risks in River Sediments from Calabria (Southern Italy). Appl. Sci. 2021, 11, 1729. https://0-doi-org.brum.beds.ac.uk/10.3390/app11041729

AMA Style

Caridi F, Di Bella M, Sabatino G, Belmusto G, Fede MR, Romano D, Italiano F, Mottese AF. Assessment of Natural Radioactivity and Radiological Risks in River Sediments from Calabria (Southern Italy). Applied Sciences. 2021; 11(4):1729. https://0-doi-org.brum.beds.ac.uk/10.3390/app11041729

Chicago/Turabian Style

Caridi, Francesco, Marcella Di Bella, Giuseppe Sabatino, Giovanna Belmusto, Maria Rita Fede, Davide Romano, Francesco Italiano, and Antonio Francesco Mottese. 2021. "Assessment of Natural Radioactivity and Radiological Risks in River Sediments from Calabria (Southern Italy)" Applied Sciences 11, no. 4: 1729. https://0-doi-org.brum.beds.ac.uk/10.3390/app11041729

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop