Multivariate analysis and spatial distribution of suspended particulate metals of Abakaliki and Enugu in Southeastern Nigeria

Background: The atmosphere is usually contaminated by heavy metals from different natural and anthropogenic sources. The levels of heavy metals in the atmosphere often become harmful and cause environmental pollution with serious health consequences to humans when their concentrations exceed certain threshold limits. The determination of the atmospherically suspended particulate matter in general and toxic trace metallic concentrations in particular, is, therefore, necessary for the evaluation of the sources and fate of atmospherically deposited suspended particulate matter. The assessment of the concentration levels of suspended particulate trace metals like Zn and Pb in the cities of Enugu and Abakaliki southeastern Nigeria was performed to assess the concentrations of atmospherically suspended metals in the area, their sources, and their pollution status. Results: Results of the study showed that Zn concentration levels were not significant at all the sampling points, while the concentrations of Pb particulate matter which ranged from 0.11 to 0.91 μg/m exceeded the international revised EPA Pb standard of 0.15 μg/m but within the approved Federal Ministry of Environment (FMEnv) Nigeria, 24 h long duration of 1 μg/m, and National Ambient Air Quality Standard (NAAQS) 1 h duration (1 μg/m). The enrichment factor (EF) showed that Pb had a dominant crustal source, while Zn was determined to be of minor to moderate crust enrichment. The multivariate analysis carried out using principal component analysis (PCA) based on varimax rotation and Hierarchical Cluster Analysis (HCA) plotted in biplot and dendrogram diagrams revealed the levels of dependency and grouping, whereas Pearson’s Correlation Coefficient (PC) revealed a strong link between the two metallic pollutants Conclusion: This study therefore which applied statistical and pollution assessment of atmospherically suspended particulate matter has revealed the importance of monitoring atmospherically suspended materials resulting from anthropogenic sources like mining, quarrying, construction, etc. Similarly, the use of statistical techniques like PCA, HCA, and Pearson’s Correlation Coefficient to the study of ambient air quality has proved to be a useful tool in the assessment of atmospheric pollution. The findings of the present study therefore will be invaluable in the monitoring of particulate matter pollution of the atmosphere and their control measures, especially in mining districts.

to a sustainable environment and is, therefore, a public health concern (Ibe et al. 2017;Ukaogo et al. 2020;Fagorite et al. 2021). Although several environmental measures are usually put in place at preferred locations in some major cities especially in developed countries to reduce and control the problem of air pollution problem, however many other emerging cities and old towns in developing countries are still facing these air pollution challenges (Ukaogo et al. 2020;Fagorite et al. 2021). The air pollutants can be transported to become part of biochemical cycles in the environment and may even enter the human body via the food chain as part of less harmful reactive products or harmless metabolites (Cimboláková et al. 2020;Shah 2021). They also may form part of the toxic environmental chemicals, that are persistent, toxic, with bioaccumulative potentials which have therefore necessitated the general public concerns to reduce their impact on human health and the environment (Ali et al. 2019;Briffa et al. 2020).
Therefore several researchers have worked on various processes, techniques, and methods for the removal of heavy metals in the environment using a combination of adsorbents, and noted that adsorbent mixtures have high removal efficiencies for heavy metals (El-Shamy et al. 2017a, b, 2018Migahed et al. 2017;Shehata et al. 2019). Other research works also revealed that they improve corrosion-resistant rate for mild steel after adding sodium haxa meta phosphate, especially in steel pipelines. The proposition, therefore, is that these research areas are evolving; and plant biomass, meso-porous adsorbents, and clays are some of the key adsorbents for heavy metal removal, metal recovery in flow systems for both batch and continuous processes (El-Shamy et al. 2017a, b, 2018Migahed et al. 2017;Shehata et al. 2019). A similar challenge has been observed in the atmosphere as a result of atmospherically suspended and deposited particulate matter and metals, which are usually hazardous (Awan and Ahmed 2011). The Suspended Particulate Matter (SPM) present in the environment is well-known to cause harmful effects because the human tissue can absorb them via breathing. Additionally, when these suspended particulates containing some heavy metals in different forms are ingested or inhaled, they can result in very damaging health hazards (Obioh et al. 2005;Gray et al. 2015). Then again, climate and vegetation are also affected by SPM and this most often may lead to the occurrence of microbial pollution in the soil and the entire ecosystem (Sahu and Basti 2021). They can also cause photochemical smog, haze, and mists when the environment exceeds the natural capacity to dissipate, dilute or even absorb them, hence the need to study these atmospherically suspended pollutants (Chen et al. 2021).
The major sources of trace metallic elements in the atmosphere are often related to combustion, mining, and construction which result from particulate emission sources according to Mukherjee and Agrawal (2017). The main source of these suspended particulate matters and atmospherically suspended trace metals (Pb, Zn, Cd, etc.) are usually associated with anthropogenic sources rather than geogenic sources (Mukherjee and Agrawal 2017). This source has such components as traffic, transportation, mineral exploration, gas flaring, and power generation which have played dominant roles in air pollution studies (Maduna and Tomasic 2017).
Though the incidences of environmental pollution caused by heavy metals have been studied in different areas of the world (Abdul-Wahab 2006; Awan and Ahmed 2011;Adam 2013;Ikamaise et al. 2013;Seibokuro et al. 2018;Sintorini 2018;Rodriguez et al. 2019), there is however a paucity of data on the pollution and ecological risk assessment of atmospherically suspended trace metals. In Nigeria and especially within the study area, there is virtually no available literature on the challenges of SPM and heavy metals in the atmosphere. However, a lot of works have been done on ambient air quality study especially with respect to the gaseous criteria pollutants (Ibe et al. 2017). Similarly, the application of geostatistical techniques and pollution assessment models in evaluating the metallic concentrations, source, and fate of atmospherically suspended heavy metals in the atmosphere is quite new and has not been previously studied in the region. In addition, a thorough review of the studies revealed that pollutant concentration levels vary with the sources and distances to the pollutant source. Also, there is a general paradigm that pollutant concentration levels tend to escalate in areas witnessing emerging industrialization and urbanization owing to typical emission control problems in those cities. Hence, a regular comprehensive examination of particle geochemistry and ambient air quality change is pertinent (Douglas 2009;Basha et al. 2010;Maduna and Tomasic 2017).
Enugu and Abakaliki towns of southeastern Nigeria are located in the former east-central state and are associated with mining activities like quarrying. Minerals mined in the mining districts of Enugu and environs are largely associated with Coal while minerals mined from the several mining districts of Abakiliki and environs include lead (Pb), Zinc (Zn), salt, etc. Whereas Enugu is the Coal City and the former capital of defunct Southeastern Nigeria there is massive urbanization and associated activities in Abakaliki. Currently, both cities are capital cities of Enugu and Ebonyi states respectively with the attendant population increase, vehicular activities, industrial activities, etc. There have been several pollution sources resulting in high loadings of pollutants in the Page 3 of 13 Ichu et al. Bull Natl Res Cent (2021) 45:159 ambient air within these cities (Ngele and Onwu 2015;Awokola et al. 2020). Because of the potentially detrimental effects of these toxic suspended heavy metals when in high concentrations in the atmosphere, the main objective of this work is, therefore, to assess the levels of Zn and Pb suspended particulate trace metals in the study area, their primary sources and their ecological and pollution risks. Both areas are associated with high levels of air pollution owing to several anthropogenic activities like mining. Despite the huge risk of atmospherically suspended metallic pollutants in this region, little or no studies have been done in this region to assess the metallic concentrations and their sources. Also, there is a scarcity of the application of statistical methods and pollution assessment techniques in the analysis of SPM in the region. Hence, various statistical techniques were used in evaluating the SPM/trace metals for spatial and multivariate analysis of the pollutants for a significant data reduction and interpretation of datasets.

Methods
In the present study, samples trapped within the study area on glass fiber filters by a high-volume sampler for quantification of Zn and Pb metals were analyzed using Atomic Absorption Spectroscopy via wet HNO 3 acid digestion. Also. Arc GIS version 9.0 was used for mapping and geospatial analysis (See Additional file 1).

Study area
Abakaliki (6°20′N and 8°06′E)-lead mining and farming city, located in Ebonyi State and Enugu (6°30′N and 7°30′E)-the city of Coal, located in Enugu State, both in the Southeast geopolitical zone of Nigeria. Both cities are currently the capital cities of Enugu and Ebonyi states. Abakaliki and Enugu have rich cultural heritage and thriving civilizations, and ethnic settlements in both states. As of the 2006 Nigeria census, Abakaliki had a population of 141,438 of which constitutes 48.45% males and 51.55% females, while Enugu had a population of 772, 664 in 2006 in which 49.2% were females and 52.25% males (FRN 2009). Abakaliki has an average annual precipitation of 1250.7 mm per year while the Enugu region has an average annual rainfall of 2000 mm with May and October as the wettest months. The mean highest and lowest temperatures are 33-27 °C and 29-26 °C for Abakaliki and Enugu and Abakaliki respectively. Ebonyi state is a rich agricultural area and a leading producer of potatoes, yam, maize, cassava, rice, and beans within the region, with few commercial mines, and several solid mineral resources including zinc, lead metal, etc. On the contrary, the first European settlers in 1909 arrived in Enugu and made the first coal shipment by 1914 which largely fueled growth and development in the state. It is primarily a public sector-driven economy with few scattered industries food and beverage, brewery, cement, pharmaceutical and pottery, and machinery productions, and some mining activities (Ogbonna et al. 2014;Aloha et al. 2016).

Sampling and analysis
Ambient air monitoring was carried out at five different sites in the two locations namely ACM, ARS, AIN, AMN, and ARR for Abakaliki locations, and ECM, ERS, EIN, ECC, and ERR in Enugu locations and the references shown in Table 1. Note that from Table 1 the "A" and "E" stand for Abakaliki and Enugu, while the CM, RS, IN, MN, and RR represents the commercial area, residential area, industrial area, mining area, and rural area respectively. Also, Fig. 1 is a representation of the sampling area within the study group, and locations adapted from the 2021 Google Earth Map and the location map adapted from the Nigerian Geological Survey Agency (NGSA) at 1:5000 and 1:10,000 scale respectively.
The determination of Zinc and Lead in the particulate matter (dust) was collected on glass fiber filter paper by Whatman within a sample area profile of 50 × 60 cm 2 (0.2 mm pore-size and about 98% collection efficiency) using a particulate monitor manufactured by Environmental Devices Corp., USA, operated at a constant flow rate (1.3 m 3 /min) and programmed to collect 12 h. samples. Using (U.S EPA 1990) model, the volume sampler was placed at a height of 4 m above ground level, while the recorded air volume passing through each filter was 1500 m 3 and glass fiber filter papers were weighed before and after sampling. Hence, the mass of air passing through a sampler in 12 h is about 60-120 mg. It is important to note that the procedure for sampling was not followed on certain days due to excessive wind, heat, or equipment breakdown; however, it was rescheduled to different days and samples collected. A total of 20 samples were taken for each of the pollutants at the 10 different locations during August when the rains subsided for 1-2 weeks. An overview of the meteorological conditions is presented in Table 1.
After sampling, the filter paper was placed inside dark and air-tight plastic containers and taken to the laboratory for analysis of lead and zinc metals. Using the modified procedures described by Ercilla-Montserrat et al. (2018) and Basha et al. (2010), the glass filter papers were initially washed with deionized water, before drying in the oven at 105 °C to a constant weight. Furthermore, the filter papers were sliced into pieces and introduced into a plastic bottle containing 125 mL of 20% nitric acid solution placed inside a water bath shaker at 60 °C for 300 rpm. The obtained solution was poured into Page 4 of 13 Ichu et al. Bull Natl Res Cent (2021) 45:159 a 250-mL borosilicate glass beaker and the bottle was rinsed 3 times with the deionized water. The borosilicate beaker was then placed over a heated flame at 150 °C temperature to evaporate the nitric acid solution to 20 mL. Furthermore, another 20 mL of analytical grade nitric acid was added to the beaker and boiling continued till a 10 mL volume was left. The aliquot was now filtered and filtrate stored in a 100 mL amber bottle in the dark until further analysis was conducted (Basha et al. 2007). The blank values were calculated by analyzing unused filter papers via the procedure described above. The samples were analyzed using (AAS) Atomic Absorption Spectroscopy (Buck Model 210 VGP) equipped with a digital readout system (U.S EPA 1999). The flow rates of air and acetylene from the storage cylinder were adjusted and analysis was performed accordingly. The concentrations of Pb and Zn in the samples were calculated by comparison with the standard curves of the respective metals. The detection limits of the chemical analysis for Pb and Zn were 0.085 and 0.052 mg/L respectively. Moreover, for each measurement at a particular site within a location, the mean was taken after sampling three times, and the calculated difference between the three determinants was less than 4.5%. The calculated recovery efficiency for AAS was 94-97% of a known amount of the metal spiked unto clean filter papers before extraction and recovery. Specific hollow cathode lamps were used for each metal at their respective wavelengths (Gharaibeh et al. 2010). Table 1 shows the referenced location of the air monitoring sampling area. The particulate heavy metal concentrations in Abakaliki and a reference station ARR during August break are summarized in Table 2. Also, the particulate heavy metals concentrations in Enugu and a reference station ERR during August break are shown in Table 3. While the Pearsons Correlation Table is captured in Table 4, the mean of the three readings for each sampling site, median, and standard deviation as shown in the Tables. The spatial distribution of Pb within Abakaliki and Enugu is shown in Figs. 2 and 3 while the distribution of Zn metal in the sampling points is depicted in Fig. 4. The statistical biplot showing principal components and data subsets is shown in Fig. 5. Also, Fig. 6 represents the plot of the concentrations of Pb and Zn at Abakaliki and Enugu sampling sites respectively. Lastly, again the dendrogram showing the concentrations at the different sampling sites is shown in Fig. 7.

Heavy metal concentrations
The numerical summary of metal concentration data obtained from ten monitoring sites including the  Ichu et al. Bull Natl Res Cent (2021) 45:159 reference sites during the August break study period is presented in Tables 2 and 3. The mean metal concentrations (Zn) remained at 0.030 µg/m 3 at the reference sites ARR during the August break but were determined to be below the detection limit at the reference site ERR during the same period (Table 2) exhibiting no significant spatial variation. The mean concentration of Pb metal was 0.21 µg/m 3 and 0.18 µg/m 3 in ARR and ERR respectively, and therefore do not differ significantly.
The highest levels of Pb and Zn are found to be 0.91 µg/ m 3 (AMN) and 0.050 µg/m 3 (ECM) respectively. The highest Pb level at ARR and ERR is close to the approved FMEnv 24 h. long duration of 1 µg/m 3 , NAAQS for 1 h. duration (1 ug/m 3 ), and 1.4 µg/m 3 of (NESREA) National environmental standards and regulation agency, but it is a higher value when compared to the revised EPA lead standard of 0.15 µg/m 3 as shown in Fig. 1 (Obioh et al. 2005;Aloha et al. 2016;Frank et al. 2019). This toxic level of Pb puts the vulnerable human populace, especially children and the terrestrial ecosystem at risk. Research has shown that Pb exposures negatively affect the oxygen-carrying capacity of the blood and also result in decreased growth and reproduction among plants and animals in the ecosystems (Frank et al. 2019;Benibo et al. 2020;Ghosh et al. 2021). The sources of Pb within that area are mainly from lead-zinc mining and the heavyduty trucks that transport goods and raw materials from one town to the other. The heavy metals in SPM values for Page 6 of 13 Ichu et al. Bull Natl Res Cent (2021) 45:159 the measuring period ranged from 0.023 to 0.040 µg/m 3 for Zn and 0.11 to 0.91 µg/m 3 for Pb, which were about many times higher than EPA reference values but a little lower than NAAQS and NESREA recommendations as depicted in Fig. 2. Zinc metal is not listed as an occupational hazard neither is it carcinogenic; hence there is no established national ambient air quality standard for zinc except ingestion of zinc metal fumes for long durations. Also, the graph in Fig. 4 shows that zinc concentration is fairly constant and unaffected by relative humidity and wind speed in Enugu as similarly observed in Abakaliki. The mean concentrations of the metals derived from the values of each site are shown in Tables 2 and 3. In comparison, the reference sites ERR and ARR had lower mean concentrations for Pb at all sites (0.21 µg/m 3 ). The mean concentration of Zn particulate matter was 0.030 µg/m 3 at both reference sites (ARR and ERR) as seen in Fig. 4.
The mean concentration of Pb was found to be over two orders of magnitude than the EPA standard. Thus, in this study, the particulate metal concentration levels can be listed in the order of ranking as follows: Zn at Abakaliki (ARR < AIN < ACM < ARS < AMN) and Zn at Enugu (ERR < ERS < EIN < EMN < ECM); Pb at Abakaliki (ARS < ARR < AIN < ACM < AMN) and Pb at Enugu (ERS < ERR < EIN < ERS < ECM). Hence the magnitude of metal order can be divided into values ≤ 0.025 ≤ 0.05 µg/ m 3 for Zn and ≤ 0.05 ≤ 0.1 µg/m 3 for Pb. About 50% of zinc concentration levels were lower than 0.025 µg/m 3 while the remaining 50% were less than 0.05 µg/m 3 . The Pb grouping also shows that all sites exceeded 0.05 µg/m 3 concentration levels. Thus, the relative arrangement shows no specific ordering for the metals but the dominance of Pb within the study sites. This result shows that Pb is the highest contributor of particulate matter load in the environment within the measured locations (Offor et al. 2016;Aloha et al. 2017;Ojekunle et al. 2018;Ichu et al. 2021).

Pearson's correlation
The Pearson's correlation (PC) is generally used to describe the proportion of the total variance in the obtained data and is explained by a linear model of the variables under consideration. The Pearson's correlation ranges from − 1 to 1, and higher absolute values indicate better dependency among the variables. The formula for Pearson's correlation between two random variables x and y is given below in Eq. 1 as: where σ x is known as the standard deviation of the variable x, σ 2 x is known as the variance of x, σ y is known as (1) p xy = σ 2 xy σ 2 x σ 2 y   Page 7 of 13 Ichu et al. Bull Natl Res Cent (2021) 45:159 the standard deviation of the variable y, σ 2 y is known as the variance of y and σ 2 xy is known as the covariance of the variables of x and y.
When two variables are found to be independent, Pearson's correlation coefficient will be 0 but the inverse is sometimes not always true because there are instances   Ichu et al. Bull Natl Res Cent (2021) 45:159 that have dependent variables and Pearson's correlation is unable to detect the dependency. Although it can exist, because Pearson's correlation recognizes the linear correlations between variables that are jointly, and normally distributed. However, two variables can be correlated but unable to follow linear correlation and Pearson's correlation failing to detect the dependency. One very useful factor in Pearson's correlation is the measurement of dependency is the p value. It provides reliable information about the probability that a given dataset will be contrary to the hypothesis, in other words, the determined Pearson's correlation is not significant despite the value of P 2 xy . Thus, it is generally accepted that a p value ≤ of 0.05 shows a significant correlation and a p value > 0.05 is not a significant correlation. Also, the coefficient of determination ′, will give the percentage of PC between two variables. For instance, a P xy value of 0.5 and P 2 xy value of 0.25 means that 25% of y values are explained by values of x (Bermudez-Edo et al. 2018;Hůnová 2020).
The results from Pearson's correlation data are presented in Table 4 above. The p value is greater than 0.05 at Abakaliki and Abakaliki/Enugu (Zn/Zn) which means that the correlation is not significant, hence the data is inconsistent with the hypothesis. But at all the sites in Enugu and Enugu/Abakaliki (Pb/Pb), the p value was less than 0.05 which revealed that it is significant showing that the data is consistent with the hypothesis. The Pearson's correlation (r) values of Abakaliki and Abakaliki/ Enugu also reflect a statistical negative correlation values and indicate that the levels of particulate heavy metals in Abakaliki differ from each other and are reasonable to study each as an independent variable. They also indicate a weak link for air pollutants moving from Abakaliki to Enugu. However, at Enugu and Enugu/Abakaliki sites it appears to contribute significantly (positive correlation) to the recorded concentrations. This also confirms that a strong link exists between the two pollutants at Enugu and the potential movement of pollutants from Enugu to   Ichu et al. Bull Natl Res Cent (2021) 45:159 Abakaliki. Hence, selecting a single variable that incorporates two or more particulate matter is justifiable.

Principal component analysis (PCA)
The Principal Component Analysis (PCA) is known to be a powerful tool commonly used to reduce the dimension of data. The identified principal components are typically expected to account for most of the variability of the datasets acquired from the different samples. Hence, the PCA uses the identified principal component to differentiate samples or subsets of samples that are seen to be responsible for the variations between the different groups. Moreover, the contribution from each data set can be visualized in the scatter plot of any two principal components to reveal the relative variation of the unrelated data sets (Yu et al. 2016). In this study, PCA was conducted by dividing the datasets into two subsets one for Zn and the other for Pb. As can be seen in Fig. 5, it depicted the biplot of two principal components of datasets. Using the varimax rotation criteria, new variable vari-factors (VFs) were formed and reduced to PCs after rotating the obtained PCs (Al-Anzi et al. 2016). The Eigenvalues of the correlation matrix provided information on how many principal components should be considered because it is the principal component that accounts for most of the variance in the observed variable. Hence in the first obtained Eigenvalue as shown in Fig. 3 only 54.42% of data is covered while the second Eigenvalue covered the remaining 45.58% of data. A view of the screen plot shows how many principal components that are to be considered. It shows a visual indication of only 2 principal components to be considered which were adopted in the plot (Simeonov et al. 2004). From the extracted Eigenvalues, the coefficient principal component 1 is mainly the presence of Zn and Pb particulate metals while the coefficient of principal component 2 was mainly Pb particulate matter. The biplot displays both the loadings and the scores for the two selected components. The scatter plots show the different sample locations within Enugu and Abakaliki. Hence sampling site ACM 3 relies on Zn as the principal component which indicates that the primary pollutant at ACM 3 is Zn particulate matter while, ACM 4,EMN 2,EMN 4) are clustered around the Pb particulate matter which indicated that the primary source of those sampling sites is the Lead (Pb) particulate matter. But for locations like EIN 2,ERS 2,and ARR 2) their primary pollutant source cannot be established from the PCA. These observations from the data correlate with the simple plot shown Page 10 of 13 Ichu et al. Bull Natl Res Cent (2021) 45:159 below in Fig. 6 where the Zn particulate matter at ACM 1-4 was not statistically indicative and identified to be the primary pollutant (principal component 1) at ACM, while the Pb particulate matter at ACM 1-4 and ECM 1-4 were significantly obvious, hence the primary pollutants (principal component) at these sampling sites (Jassim et al. 2018).

Cluster analysis
Hierarchical Cluster Analysis (HCA) is a useful multivariate tool for finding patterns and groupings within a given data set; each one of them representing concentrations of the particulate matter (Zn and Pb) (Núñez-Alonso et al. 2019). Cluster analysis involves splitting a given set of data into several group observations with unique characteristics in terms of common values or attributes of the group. Hence hierarchical cluster analysis aims to maximize between-group variance and further minimize within variance in the same group. A major advantage is that any number of variables can be used to group members of the given sample (Saksena et al. 2003).
Hierarchical cluster analysis will show how these particulate heavy metals and sampling sites relate to each other and thus plotted in a dendrogram using origin 2.0 Software as shown in Fig. 6. A two-cluster solution was selected to perform analysis for each particulate matter. The dendrogram shows that there are five primary clusters (cluster 1, cluster 2, cluster 3, cluster 4, and cluster 5) built into two other clusters when an imaginary straight line is drawn across from 0.2 scalings. The cluster also shows that these particulate metals relate with each other beyond their sampling sites. These can be seen for instance in Fig. 7 that sample points 1, 3, and 7 (cluster 1) is closely related to sampling points 5 and 21 on the same cluster. Similarly, on the far right-hand side (cluster 4), it shows that sampling points 2 and 30 are more closely related than 27 and 31 in cluster 5. Also, from cluster 2, sampling points 11, 13, and 15 are commonly related to 10, 14, and 16 than any data in cluster 1 and cluster 3 and 5. Hence there are multiple levels of clustering showing similarity of measurement and relationship in terms of levels of concentration among the data sets regardless of the measured particulate matter. Thus cluster 2 is closer in similarity to any of cluster 1 and cluster 3 than 4 and 5 because the distance between them is almost twice.
However, it took half of the dendrogram to join cluster 1 and cluster 2 together which means that they are quite different but more closely related to each other in terms of the levels of concentration than cluster 3 which took 3/4 of the dendrogram. Finally, the distance needed for clusters 4 and 5 to join clusters 1, 2, and 3 took the other half of the dendrogram which means that the concentration levels from clusters 4 and 5 are entirely dissimilar from other sampling points.

Enrichment factor (EF)
The Enrichment factor (EF) is applied in the statistical analysis of data to identify the anthropogenic source of metallic pollutants. It is useful in the determination of the degree of enrichment of a particular element compared to the relative abundance of that element in crustal material (Kothai et al. 2011). From our study, the crustal EF is calculated using Iron (Fe) as the reference material. The elemental compositions of soil in Enugu and Abakaliki are used from Chibuike et al. (2019) and (Okolo et al. 2013), whereas the aerosol levels are also calculated from Offor et al. (2016). The formula used for enrichment calculation is EF i = (i/j) air /(i/j) crust . where EF i is the enrichment factor of Zn and Pb species of I, J. where I air is the content of the measured species in the examined environment, and J air is the content of the examined metal in the reference environment. I crust is the value of the reference element in the examined environment and J crust is the value of the reference element in the reference environment (Shaari et al. 2015). From our results, all the EF for Pb at all stations in Enugu was below 9 suggesting that the crustal source is dominant. While EF level at Abakaliki was between 10 and 11 meaning that the non-crustal source is the dominant source. For all calculated EF, Zn values were between 3 and 5 in both towns which indicate minor to moderate enrichment of crustal source.

Pollution indices
The overall suspended particulate matter pollution indices and their severity were determined using the Contamination Factor (CF) and Pollution Load Index (PLI). The CF and PLI were calculated using the straightforward mathematical expression by muller in 1969 and Tomlinson in 1980. Detailed explanations can be found in recent publications by Nagarajan et al. (2019) and Kowalska et al. (2018). The identification of contamination level of trace/heavy metal in the environmental matrix analysis is determined using Cf, while the PLI factors into consideration, the potential contribution of all elements to indicate the extent of pollution in a particular location. In this analysis, the soil metal concentrations and soil background concentrations were obtained using methods described by Basha et al. (2010) and Ercilla-Montserrat et al. (2018).
The results showed that CF for (Zn) at any of the Abakaliki sampling sites was 1 = CF < 3 which implies moderate contamination, nevertheless CF for (Zn) at Enugu sampling sites revealed CF < 1 at any sampling site and refers to low contamination levels. The calculated PLI (Pb) values for the entire Abakaliki and Enugu location were 1.318 and 1.603 (> 1), hence an indication of Pb pollution. On the other hand, PLI (Zn) values for the entire Abakaliki was > 1 (1.145) while PLI (Zn) values for the entire Enugu (0.0381) were less < 1 and indicates no pollution. The overall effect is that Abakaliki has Zn particulate matter enrichment, while the entire sampling location and site is at risk of Pb particulate matter pollution. The findings depict the relationship that exist between polluted soil and atmospheric suspended particulate matter in the study location.

Conclusions
The concentration levels of suspended Zn and Pb particulate matter were investigated at 10 sampling sites, during the August break period of the rainy season in Abakaliki and Enugu towns of Southeastern, Nigeria. Pb concentrations levels almost exhibited significant differences among the sites, owing to distance and source of emission while Zn levels showed no significant levels at all locations. The most dominant concentration levels were at AMN 1-4, ECM 1-4, and ECC 1-4 which represents Onu Nwafor/Okposhi Umughara quarry in Abakaliki and, Ogbete main market in Enugu and Coal camp in Enugu respectively. The results of PCA, HCA, and Pearson's correlation applied to ambient air quality monitoring study have proven to be useful towards the assessment of air quality not only showing primary pollutants but also the strong link that exists between Zn and Pb pollutants at Enugu and Abakaliki, which provided a deeper understanding into some of the mechanism involved. The enrichment factor showed that Pb had a crustal dominant source. The final results can be useful for local authorities in designing and planning particulate matter air pollution control measures for heavy metals during the August break. Hence biannual study is recommended in future work.