Time-dependent multivariate and spectroscopic characterisation of oil residue in Niger Delta soil

In this paper, we present a detailed evaluation of changes in the oil residue in soil following a spill using weathering indices obtained from analytical instruments such as UV, IR, GC, and 1H NMR, and chemometrics based on the time of spill in the Niger Delta region of Nigeria. UV, IR and 1H NMR spectra of eight (8) oil residue samples were analyzed. The PCA of the UV and IR spectrometric index showed that the first two PCs accounted for 87 and 71% of the variance of the index, respectively. The detailed results suggested that the absorption ratios A225/256 and A248/278 from UV were good estimators for petroleum of different weathering profiles and the presence of different types of di- and poly-aromatics, nitrogen, sulphur, and oxygen (NSO) containing compounds. Similarly, sulphoxide, aromatic, and carbonyl index obtained from IR would be more valuable in evaluating changes in oil residue over time. An 84% PC obtained for NMR indicators described for weathered crude oil was the best at explaining structural changes compared to the region defined for fresh heavy crude oil. These models showed good predictive ability for the crude-oil residue composition and could be used to provide a rapid assessment of compositional differences in crude-oil residue following a spill.


Introduction
The realisation that crude oil contamination of the environment is extensive and signicant has emerged largely as a result of increased awareness with improved evaluation methods. Statistics showed that there is a gradual decrease in the number of occurrence of oil spills, however, medium and large-scale spills persist globally. 1 The oil-fed economic and industrial growth that has occurred in Nigeria has caused environmental damage. 2 Studies have shown that the quantity of oil spilled over 50 years was put at 9-13 million barrels, an equivalent of 50 Exxon Valdez spills. 3 Knowledge of the distribution of the major structural classes of hydrocarbons present in crude oil residue is therefore needed for effective implementation of remediation strategies for crude oil spill polluted sites. Once entering the environment, an oil spill is subject to a variety of weathering processes such as evaporation, dissolution, dispersion, ushing owing to wave energy, emulsication, photochemical oxidation, microbial biodegradation, adsorption to suspended matter, and deposition onto surfaces. [4][5][6][7] Thus, characterizations of the individual molecular components of the degraded crude oil are quite challenging due to the complex composition of crude oils.
Researchers have used analytical and physical chemistry techniques to monitor the degradation and fate of crude oil in the environment. Techniques such as gas chromatography (GC), gas chromatography-mass spectrometry (GC-MS) 8 highperformance liquid chromatography (HPLC), thin layer chromatography (TLC), 9 and ultraviolet-visible spectroscopy (UV) 10 have been employed in the analysis of spilled oil. Other studies have combined these analytical techniques with pattern recognition methods such as principal component analysis (PCA) and hierarchical cluster analysis (HCA) to classify weathered oil. [11][12][13] This has driven studies of aged polluted systems to evaluate the compositional changes of the oil residue following a spill. 4,14 As a spectroscopic tool, UV-visible measurement is useful for many pure chemicals. However, complex mixtures like crude oils with high absorption coefficients are also amenable to UV application. 15 This technique is increasingly being employed for in-eld applications for laboratory studies of crude oils and petroleum asphaltenes. [16][17][18][19] Direct quantication of functional groups in asphaltenes in heavy crude oil and the effect of dispersants and additives on their blends have also been studied using UV spectroscopy. [20][21][22][23] UV studies have shown that the aromatic moiety through p-p* and dipole interactions, is one of the dominant contributors to asphaltene self-association 24 and have been used to match weathered oil with source of spill 25 as well as determine compositional changes in persistent and photochemically transformed water-soluble constituents of industrial crude oil and natural seep oil in seawater. 26 Absorption maxima and absorption ratio derived from UV measurements have been used [27][28][29][30] to differentiate between SARA fractions, fresh oil, and degraded samples as well as determine the probable type of hydrocarbon in seawater. These ratios have been used successfully to differentiate the classes of aromatic (low molecular hydrocarbon with nearly identical aromatic values, mono-aromatic hydrocarbon, and diand polyaromatic hydrocarbons) found in the samples. Fourier Transform Infrared spectroscopy (FTIR) presents another unique approach to determining the changes in composition and distribution of certain chemical functionalities or classes of compounds in crude oil. 31,32 It requires minimum sample treatment and the acquisition of spectra takes a few minutes and offers acceptable repeatability in liquid samples. 33,34 IR spectroscopy has been applied in agriculture, 35 fuels and energy, 36,37 and petroleum geochemistry. 13,32,[38][39][40][41][42] Weathering indices calculated from aliphatic and aromatic biomarkers derived from the oil residue were used to ascertain the sources of hydrocarbon contamination and the extent of degradation of the spilled oils. 43 These biomarkers are highly reliable because of their specicity, hydrophobic nature, and long residence time in the environment. 44,45 Gas chromatography (GC) or gas chromatography-mass spectroscopy (GC-MS) sometimes gives poor resolution due to severe overlaps of the numerous peaks. 46 Unaltered constituents of crude oil cannot be characterized by direct GC analysis, hence the derivatization or other treatments are usually needed. The need for high accuracy, precision, quality assurance/quality control measures required in the use of chromatographic technique oen makes the use of GC tasking. NMR spectroscopy has been proposed as a powerful and highly useful technique for the structural elucidation of various crude oil residues and their fractions following spillages. 47 Although NMR spectroscopy is well known in organic chemistry, it has only been applied to the study of crude oil fractions in the last few decades. [47][48][49][50] It is explicit in describing structural parameters of crude oil and its fractions. [51][52][53][54][55] Multidimensional NMR experiments such as correlation spectroscopy (COSY), heteronuclear correlation (HETCOR), and heteronuclear single quantum coherence (HSQC) spectroscopy have yielded important information about 1 H-1 H and 1 H-13 C connectivity in petroleum molecules as well as the chemistry of crude oil with respect to functional groups present. 55,56 As a result of the variations in environmental conditions and time, true composition of oil residue varies as it ages in soil. Therefore, there is a need to measure eld samples as it is impossible to reproduce similar conditions in the laboratory. The Niger Delta region provides an excellent background for such study especially for clean-up and remediation work. [57][58][59] Thus, this research aims to (a) determine the level of degradation using weathering index for oil residue obtained from UV, IR, NMR spectroscopic tools, and gas chromatography, and (b) develop a model that will give an instantaneous view of the structural composition and changes in crude oil residues following a spill by using principal component analysis.

Study area and sample collection
Niger Delta region in Nigeria is the delta of the Niger River that sits directly on the Gulf of Guinea on the Atlantic Ocean. This study was carried out in Edo, Delta, and Rivers States of the Niger Delta region of Nigeria. Fig. 1 shows the map of the Niger Delta region and sampling points. The sampling was done in May 2015. The study involved eight (8) notable oil spilled sampling sites for the survey while fresh crude oil was used as the control. The choice of sampling was based on a site with a history and time of oil spill. Soil samples (S1-S6) were obtained from the oil spilled sites, while samples coded B1 and B2 were obtained by biostimulation; a simulated process involving spiking of unpolluted soil samples with a combination of fresh crude oil and cow dung acting as the natural fertilizer. The cow dung contained petroleum utilizing bacteria to improve the degradation of the oil residue overtime. The fresh crude oil (FC) used as the control sample was also subjected to the same procedure as described for the oil residue samples.

Extraction and fractionation of oil residue
The extraction and fractionation of the oil residue from soil have been described elsewhere. 60 The extracted oil and fractions obtained were used for subsequent gas chromatographic (GC), ultraviolet-visible spectroscopy (UV), infrared spectroscopy (IR), and nuclear magnetic resonance (NMR) analyses.

Ultraviolet analysis
The absorbance of the oil residue and fresh crude oil in hexane (20 mg L À1 ) was measured on a UV-2600 spectrophotometer (Shimadzu, Japan). Selected absorption maxima and absorption ratios that dened hydrocarbon content and type were used to determine the variation of oil residue composition with the age of spill. The selected absorption ratios used were A 205 /A 215 , A 228 / A 256 , A 248 /A 267 , and A 248 /A 278 . 61,62 Standard solutions of 10 mg L À1 crude oil in hexane were prepared and successively analysed in duplicate on a UV-vis spectrophotometer in the range of 200-450 nm in order to assess the correlation between compositional changes and time of spills. The different absorbance signals (A) depicting naphthenic and aromatic compounds in the resulting sample spectra were plotted against the corresponding wavelengths (l) of 230 nm and 260 nm, respectively. 63 A linear regression curve of variations in absorbance at specic wavelength versus time was obtained; and the line equation and the coefficient of determination (R 2 ) were evaluated.

Infrared spectroscopic analysis
Attenuated Total Reectance Infra-red (ATR-IR) analyses were carried out for oil residue and fresh crude oil using a Bruker-Alpha FTIR Spectrophotometer (Bruker Optics GmbH Ettlinger, Germany). Data were recorded in the MIR spectral range of 400-4000 cm À1 at a spectral resolution of 2 cm À1 . Twentyfour (24) scans were averaged for both the backgrounds and the sample spectra. All experiments were done at ambient conditions (T ¼ 302 K). The functional and structural band areas were used to generate the IR spectrometric index based on the methods of Pieri et al. 64 and Ganz and Kalkreuth. 65 The average area in the dened region was used to calculate the spectrometric index, and the standard deviation of all measurements was found to be within the acceptable limit. The main areas used have been described by Permanyer et al. 66 and are given in Table S1. †

Gas chromatographic-ame ionization detector (GC-FID) analysis
The saturate fraction of the oil residues was analysed for total petroleum hydrocarbon (TPH) using the method previously described by Udoetok and Osuji. 67 Individual n-alkanes were identied based on the retention times of the standards of nalkane n-C8 to n-C40 (Sigma Aldrich). The concentrations of each n-alkane in the samples were subsequently determined from the calibration curves of individual n-alkanes.

1 H-NMR measurement
The 1 H-spectra were obtained using a 5% (w/w) oil residue in CDCl 3 (99.8%, Aldrich) using a 5 mm probe at a spinning rate of 10 Hz and a temperature of 301.15 K. The spectrum was recorded on a 500 MHz Bruker Avance NMR spectrometer while 30 pulses (Bruker zg30 pulse sequence) were used with a delay time of 2 s (sweep width 3900.25 Hz). Sixteen scans were averaged for each spectrum. The automated phase and baseline Nudelmann et al. 68 Rios et al. 69 Mohammad and Azeredo 70 Poveda and Molina 71 corrections were applied to get reproducible integral values. The spectra were integrated three times and the mean values were used in the calculations. To assess the repeatability of 1 H chemical shi at 500 MHz, some samples were independently determined in triplicate and a standard deviation of 0.003 ppm was obtained. The dened regions in the NMR spectra obtained were then integrated as presented in Table 1 and were used to obtain the percentage areas of each region. Table 2 describes the assignments of all the regions integrated.

Principal component analysis
Owing to the variations in the chemical compositions of the oil residue over time, the predened data set consisting of variables obtained from UV absorbance ratio, infrared spectrometric index, GC weathering index, and NMR region areas integrated for the samples were utilized for the PCA analysis.
Principal component analysis (PCA) is a combination of variables describing the major relationship in a set of data. It helps to give a projection of the most important information contained in a set of data that do not appear ordinarily. Data from the study was analyzed using OriginPro2019b Statistics soware and visualized using SIMCA 17.0 (Umetrics). Principal component analysis (PCA) was computed to identify signicant principal components in the data as well as possible loadings and similarities between the oil residue and the time of the spill.
The PCA was carried out by the varimax normalized rotation method. A scree plot and Kaiser criterion with an eigenvalue greater than 1 was performed to determine the number of important components present in the data, to provide a visual interrelationship, detect and interpret sample patterns, similarities, and disparities. 72,73 3 Results and discussions 3.1 Spectra features and absorbance ratio determined from UV-vis spectroscopy . This would suggest that a complex mixture of compounds such as saturates, aromatics, resins, and asphaltene were responsible for these absorptions. 74 Since it is difficult to assign absorption maxima to a single compound in an oil mixture, we, therefore, took advantage of the similarity in the chemical nature of fractions. For example, an absorption band owing to the differences in the aromatic composition of the sample would arise as a result of the ability of substituents present in the aromatic moiety to go into electronic interactions such as conjugation or hyperconjugation (in the case of alkyl substituent) with the p-bond system of the ring leading to Aliphatic hydrogens in methyl or methylene b or further from an aromatic ring H 3 Alicyclic hydrogens in position b to two aromatic rings (naphthenic rings and methylenes) H 4 Aliphatic hydrogens in methyl or methylene g to an aromatic ring, which can be attached in a g position to another or the same ring H 5 Olenic hydrogen H 6 H f to aromatic rings CH 2   The UV-visible spectrum of oil residues from samples S1-B2 and FC.
a shi to a longer or shorter wavelength. Also, degradations oen lead to the formation of auxochromes or groups (polar/ resin in nature) in the molecules that create an additional 2 to 5 nm to a wavelength of major absorption bands. 75 On a sample basis, younger residues B1, B2, S5, and S6 (Table S1 †) show absorption maxima between 226-230 nm and 254-262 nm. This suggests the presence of unsubstituted polyaromatic compounds such as pyrene (245 nm), anthracene (255 nm), and phenanthrene (255 nm) as the major constituents. 76 Older residues such as samples S1, S2, and S3 have absorption maxima at 278, 274, and 264 nm, respectively compared to fresh crude oil with absorption maxima at 256 nm. It could be that these samples contained alkylated polyaromatic compounds with contributions from three and four aromatic systems that oen give rise to absorption at a longer wavelength. 77 In this study, four absorption ratios (Table S2 †) calculated from the UV visible spectra were used to evaluate the compositional changes in time as a result of weathering processes. The absorption ratio A 205/215 was between 1.213 and 1.525; older samples (S1, S2, and S3) had relatively constant values while the values for younger residue (<5 years) were slightly higher. These may be attributed to the presence of soluble low molecular weight hydrocarbon and aromatic compounds such as benzene, xylene, and toluene. 30 The higher values of A 205/215 obtained for fresh crude oil supported this evidence. The values calculated for A 228/256 ratio were higher than A 205/215 .
The ratio A 228/256 has been described as a good estimator for petroleum input and may be constant for oil from the same source and different for oil of other origin or weathering prole. 78 In our case, this ratio increased with the time of the spill and may be as a result of an increase in long-chain carbon or branched chain hydrocarbon attached to aromatic moieties. The third and fourth ratios dened as A 248/267 and A 248/278 supported by the R 2 value of 0.9117 in Fig. 3b were indicators of the presence of different types of di-and polyaromatics 62 and compounds containing nitrogen, sulphur, and oxygen (NSO) atoms that could cause the bathochromic shi. Fig. 3a and b show the linear relationships of absorbance of oil residue with the time of spill at 230 nm (benzenic compounds) and 262 nm (naphthenic compounds), respectively. An optimum concentration of 20 mg L À1 of the oil residue was used for the absorbance measurement because above this concentration, the oil residues began to aggregate and the maximum absorbance was disrupted as a result of noise 74 and low concentrations below 10 mg L À1 could not be used because the technique could not detect species below this concentration. 79 The absorbance values of samples S2, S3, S4, S5, B1, and B2 gave fairly acceptable and positive regression values (R 2 ) of 0.8463 and 0.9117 at 230 nm and 262 nm, respectively. In assessing the validity of the linear regression obtained, the absorbance value of the fresh crude oil gave a time scale of less than one year at 230 nm. This unexpected value gives credence to the linear relationship between the absorbance at the selected wavelength and the time of spill. The absorbance values of samples S1 and S6 (omitted in the plot) which were higher than expected might have been strongly inuenced by their compositions that caused the overlapping of the absorbances of the different constituents in them.

Spectra features and spectrometric index determined from infra-red spectroscopy
The use of ATR to simplify and evaluate mid-infrared spectroscopy (MIR) is well documented. 80 The ATR is noted for accuracy because, in MIR, the fundamental bands are specic, sharp, and sensitive to include the absorption band of aliphatic C-H bands and additional bands originating from groups containing oxygen, sulphur, and nitrogen. Table S3 † shows the assignment of the bands observed in a typical crude oil IR assay while the results of the IR spectrometric index calculated are shown in Table S4. † Differences in the structural composition were not obvious by mere visual inspection as observed in the IR spectra (Fig. 4a) of whole S1-FC oil residue samples. In order to illustrate the differences in the fractions obtained, a representative spectrum obtained from saturates, aromatics, and polar fractions of sample S1 is shown in Fig. 4b. Appearance and absence of bands at the region that denes each fraction showed the presence of functional groups in the oil residues. Generally, saturates show characteristic bands at 2962 to 2872 cm À1 as a result of the C-H stretching vibration from CH 3 groups while bending vibrations of CH 3 appear between 1475 and 1375 cm À1 . 81 The presence of straight-chain alkane is also conrmed by the band at 721 cm À1 . Evidence of the aromaticity of the sample is supported by the presence of bands at 900 to 675 cm À1 and 1605 cm À1 . A representative spectrum (Fig. 4c) of the resin fraction of sample B1 shows broad and strong bands between 1330 and 1100 cm À1 assigned to the S]O functional group 82 and 1030 cm À1 which is typical of the aromatic ethers and sulphoxide groups. 32 The resin fraction of sample B1 that contains phenols and organic acids shows bands for O-H and C]O at 3100-3600 cm À1 and 1710-1680 cm À1 , respectively. The application of the spectrometric index obtained from the IR region to PCA, yielded useful information on the pattern of distribution of NSO species in samples with the time of spill.

Gas chromatography (n-alkane distribution)
Normal-alkanes give useful information about the extent of weathering or freshness of crude oil following a spill. Samples with different spill histories can be compared by using the distribution of the n-alkanes biomarkers (pristane and phytane) and the total concentrations of the individual alkane. 83 Table  S5 † shows the GC weathering indices of the saturate fractions (F1) of the oil residue obtained from different locations. Evaporative losses removed the low molecular weight alkanes for all the samples and then biodegradation and photooxidation started to affect the nature and distribution of the classes of biomarkers, complex alkanes, and naphthenic compounds present in the oil residues. 83 Compared to the oil residues, the fresh crude oil is characterized by alkanes from C 11 -C 40 (S1) while there are no n-alkanes below C 14 in samples S1-B2. Pristane (Pr) and phytane (Ph) are ubiquitous in crude oil and are resistant to weathering 84 and have been used as indicators of petroleum contamination. 85 Pristine/phytane (Pr/ Ph) ratio ranged from 0.74-1.56, indicating petroleum origin source. 86 According to Didyk et al., 87 a high concentration of phytane could also be explained by low microbial activities in a reducing condition as is the case of older residues with Pr/Py values less than 1 except for S1F1 and values greater than 1, for younger residue. The values of Pr/n-C 17 and Ph/n-C 18 ranged from 2.7-6.6 and 2.46-6.76, respectively. This result is in agreement with the report by Faboya et al. 88 and is an indication of the loss of n-C 17 and n-C 18 relative to pristane and phytane. The carbon preference index (CPI) values calculated from the nalkane distribution in the present study ranged from 1.29-1.42. CPI values of n-alkanes close to 1 are considered to emanate from petroleum sources. 89 There was no observable or drastic change in the ratio (Pr + n-C 17 )/(Ph + n-C 18 ) for the oil residue and as such may not be useful in spotting the obvious differences in the structural compositions of the oil residue.

Principal component analysis of UV, IR, and GC index
Weathering Indices obtained from UV (Table S2 †), IR (Table  S4 †), and GC (Table S5 †) measurements relevant to the elucidation of structural and composition changes in oil residue following a spill were used in the PCA analysis. Table 3 shows the rst two PCs obtained from the spectroscopic and GC weathering index. The selection of two PC scores was appropriate considering the percent variance in each of the PCs. For the absorbance ratio, PC 1 and PC 2 accounted for 62 and 25% of the total variance, respectively. Its PC 1 was inuenced by positive values of A 248 /A 278 . (62%), A 248 /A 267 (59%) and A 228 /A 256 (45%) in that order, and represented the younger residue and biostimulated samples. PC 2 has a high correlation (87%) with A 205 /A 215 suggesting a high concentration of aliphatic compounds. 68 This shows that the oil residue and the fresh crude oil are characterized by the saturate fraction and gives reliable information on the extent of weathering, irrespective of the degradation process involved. The pattern in the biplot (Fig. 5a) clearly shows that oil residues S1, S2, and S3 are very similar in composition and are different from the younger residues such as S4, S6, and S5. These results agree with our hypothesis since they cluster together in the PCA biplot. Similarly, Fig. 5b indicates a good evolution of the sample with time in an anti-clockwise manner.
Other studies have reported the property of crude oil either by using the SARA fractions or the whole crude. 41,80 The relationship between the IR spectrometric index and the time of spill was made obvious by the PCA biplot. From our results, biostimulated samples B1 and B2 and older oil residues (S1, S2, and S3) were more of carboxylic and aliphatic acids. Samples B1 and B2 had lower P A values compared to other non-assisted oil residues (S1-S6). This decrease is consistent with the loss of more degradable hydrocarbon as microorganisms prefer small molecules. 68 In the PCA determination, no pre-classication of the oil residue was performed. The rst and second PCs represented 49.5 and 22% of the total variance, respectively. The Euclidean distance between the samples in the plot represented real differences in the IR structural index. The rst PC was largely inuenced by aromatic index (52%), aliphatic index (50%), sulphoxide index (46%), and long-chain index (37%) in that order. The second PC was inuenced by the branched index (51%) and carbonyl index (44%). Fig. 5c also shows that the most positive values were represented by older oil residues (S2, S3, and S4) with the sulphoxide content while B1 and B2 had a good correlation with carbonyl index suggesting biodegradation. The younger residues correlated with aromatic index for sample S6 suggesting the transformation of aromatic compounds. However, the unlikely correlation of sample S1 (10 years) may suggest an increase in polyaromatic hydrocarbon. This distribution into groups based on the time of spill (Fig. 5d) and a 71% obtained for the total variance suggests that IR  spectrometric index (carbonyl, sulphoxide, and aromatic index) would be a better indicator to account for structural changes in crude oil following a spill. It is usually assumed that the major transformation occurs in the rst year of the environmental exposure but the biplot shows that recent spills may have similar characteristics as the aged oil residue. PCA analysis of the GC weathering index gave an 85% total variance explained by the rst two PCs ( Table 3). The rst PC (71%) represented the younger residues and fresh crude oil. The biplot (Fig. 5a) shows how the parameters are inuenced by the samples with high Pr/ Ph, low Pr/C 17 , and Phy/C 18 ) suggesting high concentration of aliphatic compound. 68 The second PC (14%) was associated with older oil residues inuenced by CPI and (Pr + n-C 17 )/(Ph + n-C 18 ). This suggests the transformation of aliphatic compounds into a polar compound and a possible increase in the aromatic component. Fig. 6f shows the distribution of the time of spill. Transformation occurs in the rst few years of environmental exposure. 68

Nuclear magnetic resonance
Marshall and Rodgers 90 have described crude oil as the most complex compositional organic mixture. In this work, we exploited the limitation of NMR as a tool for the structural elucidation of a single-component system and applied it to the determination of the composition of a multi-component system such as oil residue. As expected, the NMR spectrum showed an overlapping prole ( Fig. 6a and b) making it difficult to differentiate between samples. Fig. 7 shows a representative integrated NMR spectrum (sample S1) using the areas described by Poveda   This clearly shows that NMR indicators described for weathered crude could best explain compositional differences than the region dened for fresh heavy crude oil.  Fig. 8a shows the PCA biplot of the sample distribution with NMR parameters described by Nudelmann et al. 68 In this case, there was a clear pattern recognition between the oil residue and the NMR structural features. Older residues (S1, S2, S3, and possibly S4) were characterized by H 2 , H 4 , and H COOH . Apart from the H COOH , there are aliphatic hydrogens in methyl or methylene in b or g position to an aromatic ring and are fairly resistant to degradation. 91 Free fatty acids with an alkane chain lengths of C 14 -C 26 range have been reported in petroleum sand samples and the typical acids found were alkanoic, mono-and diunsaturated C 18 -alkenoic acids. 92 Acid components present in   the polar fractions could be carboxylic acids for species with two or more oxygen atoms, phenols, and also NSO species, such as pyrroles, indoles, and thiols, 93 as well as esters. 94 The formation of carboxylic acids and hydroxyl functionality is also as a result of microbial degradation of petroleum hydrocarbons following several oxidative processes. 95,96 Samples B1 and B2 hold interesting features in that they are characterized mainly by H 1 and H 3 (Fig. 8a). This indicates that the already identied microbial consortium in the cow dung such as Acinetobacter spp., Bacillus spp., Pseudomonas spp., Alcaligenes spp., and Serratia spp. 97 that are heterotrophic and petroleum-utilizing in nature might have degraded the aromatics (benzenoid) and phenolic component into nutrient, biomass, and CO 2 via biodegradation and through their intermediates. 98,99 The gas chromatography-mass spectrometry and IR analyses of the polar fraction of the oil residues show that samples S1, S2, and S3, B1, and B2 contained carboxylic acids and alcohols. 100 Younger oil residues (S5 and S6) are characterized by H A , H 2 , and H 4 . These samples may contain long side-chain n-alkyl aromatic with side carbon atoms ranging from C 7 -C 27 . 49 Fig. 8b shows the time evolution of the oil residues with NMR parameters. The older residues (>5 years) clustered and younger residues (<5 years) evolved in a clockwise direction, suggesting that time evolution and environmental conditions could affect degradation in a predictable pattern. Reports by Zhou et al. 101 shows that data from both chromatographic and NMR parameters agree well with the level of biodegradation of crude oil sample from different sources. Fig. 8c shows the biplot from Rios et al. 69 NMR parameters. It appeared that older residues (S1, S2, and S3), as well as residues that might have undergone enhanced degradation (B1 and B2), were characterized by H A1 and H 3 . It is clear that samples S6 and FC described by H 2 only did not show any proximity to samples and so might not be a structural feature to consider when evaluating older residues. Fig. 8d shows the time evolution of the sample with NMR parameters from Rios et al. 69 In this spatial distribution of time, the degraded oil residues with a time of spill greater than 5 years and biostimulated samples (B1 and B2) were grouped. The maximum variance obtained from this biplot might probably be as a result of fewer NMR parameter terms used. Biplots (Fig. 8e) 71 NMR parameters (Fig. 8g) did not show any clear pattern. This was the same for the time evolution (Fig. 8h) and so was not discussed further. It has been suggested that the environmental evolution of the oil spill samples when reported may not reect the true nature of samples in the eld, but would give an instantaneous view on the level of weathering in the area under study. 6

Conclusions
This work has shown that compositional changes in crude oil following a spill at different times can be modelled using GC weathering index, UV absorption ratios, IR spectrometric indices, NMR spectroscopy and Chemometrics to give a good prediction on the functional groups in each of the oil residue irrespective of the environmental exposure or conditions. Pr/Ph ratio and CPI values close to unity from the GC index showed that the oil residues were all derived from petroleum sources. Results also suggested that the absorption ratios A 225/256 and A 248/278 from UV were good estimators for petroleum of different weathering proles while sulphoxide, aromatic and carbonyl index obtained from IR would be more valuable in evaluating changes in oil residues over time. For both the UV and IR spectroscopic measurements, two PCs were captured that explained 87% and 71% of the total variance, respectively. NMR indicated that b and g alkyl-substituted aromatic, aliphatic, and -OH and COOH type compounds provided the structural difference in oil residues with the time of the spill. The Niger Delta region provided an excellent point for this study as a result of its notable and continual oil spill. The NMR region described for weathered oil residue could be applied to obtain structural differences. PCA indicated that the time evolution of the sample followed a predictable pattern and might be useful when the age of a spill is known.

Conflicts of interest
There are no conicts to declare.