REVIEW article

Front. Environ. Sci., 09 May 2023

Sec. Soil Processes

Volume 11 - 2023 | https://doi.org/10.3389/fenvs.2023.1137496

Soil and groundwater contamination by crude oil spillage: A review and implications for remediation projects in Nigeria

  • 1. Pan African University Life and Earth Sciences Institute, Ibadan, Nigeria

  • 2. Department of Geology, University of Ibadan, Ibadan, Nigeria

  • 3. Department of Environmental Sciences, University of Toledo, Toledo, OH, United States

Abstract

Management of soil and groundwater resources has been recognized as essential to meeting the sustainable development goals of Agenda 2063 of the African Union. As Africa’s fastest growing population with over 200 million people, Nigeria is responsible for leading the continent’s environmental sustainability goal. Nigeria has seen a sizable number of crude oil spillages that have contaminated its soils and groundwater resources, and several of these contaminated sites are to be cleaned up yet. There needs to be more scientific data to design an effective cleanup and to manage the soil and groundwater resources effectively. So far, the only extensive crude oil-contaminated site remediation project documented is on Ogoniland in the Niger Delta region of Nigeria. However, this project resulted in less effective, albeit temporary, cleanup solutions. This review presents a state-of-the-art synthesis of research on soil and groundwater contamination by crude oil. It includes sections on processes, measurements, predictions, and management, as well as an analysis of the state and challenges in Nigeria. In-depth field, laboratory, and computer models for crude oil contamination investigation have been developed with over 60 years of significant research. However, studies and case projects in Nigeria have relied on point sampling to determine the concentration of crude oil contaminants in soil and groundwater. This method offers limited information on the solute concentration and hydraulic distribution, which regulates pollutant mobility within the subsurface. The absence of baseline and high-resolution subsurface characterization data has also resulted in a need for more process-based knowledge to direct the development of site-specific remediation strategies. As a result, it is challenging to design a conceptual model that is detailed enough to help with predictions of the flow dynamics of crude oil contaminants in the unsaturated and saturated zones. It is anticipated that this review will stimulate further multidisciplinary research through site evaluation and monitoring to provide reliable information that can be used to develop appropriate model-based remediation solutions.

1 Introduction

Since the discovery of Oil Creek in Northwestern Pennsylvania by Edwin Drake in 1859 (Wicks, 2009), the demand for crude oil exploration has increased over time. With new advances and technologies in conventional and unconventional exploration of oil and gas resources, crude oil has become the principal source of revenue growth in most oil-rich nations (). Following the discovery of crude oil, various exploration, production, and distribution activities have resulted in massive oil spills worldwide (). According to the International Tanker Owners Pollution Federation Limited (ITOPF) oil spill statistics, accidental large-scale oil spills contribute to a significant volume of contamination in the global environment (). The most notable oil spills include the Castillo de Bellver tanker carrying 1.5 million barrels of crude oil that capsized and sank in 420 m of water (); the release of over 6 million barrels of crude oil into the Arabian Gulf during the Gulf War (); and the release of 3.5 million barrels of crude oil from the collision of the Atlantic Empress and Aegan Captain in the Caribbean Sea (). Although oil spills have decreased in recent years, the long-term impact of the previously recorded oil spill still affects some productive and vulnerable compartments of the ecosystem ().

Crude oil spills affect different compartments of the environment, such as air, surface water, plants, soils, and groundwater systems. The near-surface soil and groundwater system (critical zone) is characterized by an elevated level of heterogeneity, which makes the introduction of crude oil with various physical and chemical properties more challenging to investigate (Theodorakis and Walker, 2014). Subsurface fluid flow dynamics are pretty heterogeneous due to variation in the soil’s physical, chemical, and biological properties. For instance, the organic matter content of the soil increases, so does its complexity, which is influenced by abiotic soil structure and biotic variation and results in diverse bacterial communities (; ). The presence of crude oil contaminants in soil and groundwater generates an additional fluid phase due to their physical qualities, resulting in a multiphase fluid flow system. The partitioning of contaminants between different phases affects the processes that govern the fate and transport of crude oil contaminants in the subsurface (; ). Investigating these processes requires a detailed site conceptualization that involves a multidisciplinary approach (Williams et al., 2006). This approach focuses on investigating the various physiochemical and biological processes and their interplay as they affect the fate and transport of crude oil contaminants (; ). A well-defined conceptual site model of contaminant flow and transport includes knowledge of the shallow subsurface physicochemical configuration and hydraulic and transport properties (Rolle et al., 2011). The conceptual model is expected to result in a hypothesis and predictions that can be used to understand the future implications of the movement of crude oil contaminants and their interactions with the sources, flow channels, and soil and groundwater systems ().

In Nigeria, thousands of records of crude oil spillage, whose causes include sabotage and theft, equipment failure, and corrosion, have been identified (). The threat of crude oil contamination to human health has resulted in multiple lawsuits, crises, and bloodshed (). The most recent historic verdict in favor of four farmers and environmental activists declared Shell Company accountable for the oil spill in Ibibio land, with an unknown amount to be paid to the farmers (). The National Oil Spill Detection and Response Agency (NOSDRA) estimated in 2021 that over 92 percent of crude oil-contaminated sites remained unremediated, with minimal scheduled cleanup. Most crude oil spills are more than 10 years old, creating a continuing source of contamination in the soil and groundwater systems (). Crude oil contamination remains a critical component of soil and groundwater contamination in Nigeria, with most oil-producing states recording considerable numbers of crude oil spillages, leading to several land and water systems remaining in deteriorated conditions, unmonitored, and remediation efforts that are currently ineffective. Until recently, there has been little research that summarized the state of knowledge on crude oil contamination in Nigeria, despite many reports of oil spillage as documented by NOSDRA. The current condition of many contaminated sites cannot be ascertained because comprehensive site investigations have not been carried out except for the recent Ogoniland cleanup, the federal government’s only significant remediation effort. Despite the United Nations Environment Programme (UNEP) outlining specific remediation strategies for contaminated soil, no particular strategy was found to be the most efficient for the Niger Delta region. (Zabbey et al., 2017). This can be attributed to a lack of detailed subsurface heterogeneity information that will boost the success rate of any remediation exercise, causing the initiative to fall short of its intended goal. A few contaminated site investigations recorded in the literature tend to focus more on soil and groundwater geochemical characterization by comparing concentration values with acceptable standard values, making a first-order assessment of contaminant sites in the Niger Delta a challenge, particularly on a regional scale (; ; ). The lack of interest from the government body is linked to the perceived economic value of groundwater resources in Nigeria. There is a need to review the practical approaches to contaminant site investigation by providing the current status of soil and groundwater processes that control the movement of crude oil contaminants.

In this study, we outline the physical, chemical, and biological processes that determine the fate and movement of crude oil contaminants in soil and groundwater systems using global examples but emphasizing the Nigerian context. This is intended to provide adequate knowledge of crude oil mass transfer, natural attenuation, and the appropriate management technique to stimulate more scientific discussion on Nigeria’s characterization and remediation of crude oil-contaminated sites. In addition, we identify possible knowledge gaps in Nigerian contaminant site characterization and remediation research. This review examines the existing conceptual model and research that explains the processes that govern the fate and transport of crude oil contaminants. We critically analyze this literature and the diverse, more specialized process-based literature that contributes to our understanding of mass transfer, crude oil transport, and attenuation in unsaturated and saturated zones.

2 Crude oil contaminant fate and transport processes in unconsolidated sediments

The evolution of crude oil contamination research begins with the assumption that crude oil introduced into the subsurface gravitationally flows down to the water table and either floats like a pancake on top of the water table, displacing nearly all the water and air in the aquifer’s pore space, or penetrates to the water table’s base before settling, depending on its density properties (; , ; ). However, studies have shown that crude oil coexists with water in the aquifer’s pore network, and the degree of crude oil contaminant saturation is highly dependent on the lithology and fluid properties (; ). Typically, this is determined from water and soil samples collected from wells and boreholes with BTEX (benzene, toluene, ethylbenzene, and xylene) concentrations significantly higher than the standard values (). Upon the release of crude oil into the subsurface, crude oil migrates downward under the influence of gravity and is laterally subjected to capillary forces (Figure 1). The capillary forces caused the contaminant mass to be immobilized within the soil pores, though lateral spreading within the unsaturated zone could be observed depending on the contaminant’s volume and the permeability of such a layer (Sethi and Di Molfetta, 2019). The downward migration is impeded upon getting to the capillary fringe as the buoyancy force and increased water content limit the extent of vertical migration, leading to lateral migration in the direction of the groundwater flow gradient ().

FIGURE 1

However, the accumulation of crude oil contaminants and radial flow can occur if the rate of crude oil contaminant movement from the surface is greater than that of the lateral migration (). The dissolved plume provides an extended source of contamination as soluble constituents are transformed into water and solid phases (). This transformation leads to changes in the mass and composition of the contaminant with time, which is directly related to the concentration variation along the flow path. These temporal and spatial transformations are governed by advection, dispersion, adsorption and desorption, volatilization, partitioning, evaporation, and biodegradation ().

2.1 Advection, dispersion, and diffusion

Advection is the most dominant solute transport mechanism in soils whose magnitude is primarily controlled by hydraulic gradient distribution along the flow path (de Vries et al., 2017). Advection describes the bulk movement of crude oil contaminants under a hydraulic head gradient (). In the unsaturated zone, the rate of the advection process has been identified to be controlled by the rate of infiltration, which is directly related to the quantity of contaminant transport in the aqueous phase toward the groundwater (; ; ). Its effect is significantly observed in layered permeable deposits with predominantly parallel bedding planes (). This process causes the interlayer solute transfer to decline along the permeable layers, thus increasing the concentration in the less permeable layers. This, in turn, leads to the development of a dispersed concentration gradient in the direction of flow (). Due to composition and temperature variation, advection transport in the unsaturated zone is greatly affected by barometric pressure water, movement through infiltration, water table fluctuation, and density gradient (). Though advection is not the most dominant contaminant transport process in the unsaturated zone, its relevance is observed in the highly permeable, relatively deep unsaturated zone (Šimůnek and Genuchten, 2016). ) revealed that a pressure-driven advection contributes up to 15% toward the net flux of crude oil contaminants in the unsaturated zone. Advection-induced flow tends to have more significant horizontal components as it is significantly driven by the lateral groundwater flow, though its effect is extensively seen at an early stage of the crude oil spill. Advection-driven crude oil contaminant migration occurs alongside dispersion and diffusion. Though the spatial influences of the latter two processes are controlled by pore-scale heterogeneity, their importance cannot be neglected (; ; ).

The irregularities of the soil pores and the flow velocity pattern within a pore influence the dispersion process (). Dispersion is commonly considered a second-order process characterized by variations in groundwater velocity around the mean advective velocity (). Under transient flow, dispersion at the pore scale tends to be visible at short intervals and extremely diffuses at long intervals. A more extensive pore size distribution generates more significant velocity differences in solute transport, resulting in a greater spread. This would subsequently contribute to the observed non-ideal behavior of crude oil contaminants. At short and long travel durations, the mixing and spreading of solute plumes traveling through saturated porous media exhibit unexpected dynamics and statistics. Recent research has demonstrated that a non-linear compound-specific parameterization of transverse dispersion is required to define transverse mixing at the Darcy scale (). In addition, transverse mixing in heterogeneous porous media lowers plume fingering caused by longitudinal dispersion (). Diffusion acting through small-scale stratigraphic heterogeneities in porous media can account for the longitudinal dispersion commonly observed in groundwater contaminant plumes and tracer pulse migration ().

Diffusion controls the mass transfer process in a non-equilibrium fluid phase. The diffusion process is the dominant mass transfer process in the unsaturated zone (You et al., 2020). Over time, diffusion tends to diminish acute disparities in concentration distributions. Transverse mixing of hydrocarbon contaminants along the plume margin can influence long-term diffusion-induced migration, posing difficulty in the upscaling of solute transport models (). Diffusion is, therefore, significant not just conceptually since it is the only mechanism that results in an efficient exchange of solute mass in fluids but also quantitatively and practically (). Under normal groundwater flow conditions, the pore diffusion term is quantitatively significant for the local hydrodynamic dispersion coefficients. This negates the assumption of ), who generalized advection and dispersion as the major driving forces in contaminant migration in the Niger Delta region. However, his work forms a base for further investigation to extend the knowledge base of the various processes that control contaminant migration within the area. In addition, when groundwater flow is primarily dominated by advection, diffusion plays a significant role in mechanical dispersion due to the development of poor mixing and compound-specific concentration gradients in the pore channels. Scheibe et al. (2015a) investigated the combined effect of pore-scale advection and intra-granular diffusion on solute transport. They concluded that pore-scale diffusive mass transfer processes between preferential flow paths and relatively immobile zones within the primary porosity might substantially impact transport, particularly in low-concentration tails. ) reported that if reaction rates or pore-size distribution variation is substantial, the magnitude of the rate of molecular diffusion across pores via the rock matrix can influence the geometry of the concentration gradients. As a result, it is expected that reaction rates and associated microbial growth will increase at pore throats in regions with shorter diffusion zones. More recently, Yoon et al. (2019) utilized calcite dissolution as an example to demonstrate how changes in pore-scale velocity reduce the reaction rate between dissolved and solid phases. They revealed that diffusion-induced contaminant transport is more visible in low-permeability zones. Therefore, diffusion explains poor results from many groundwater remediations because at contaminant sites where transmissive zones are in contact with low-permeability zones, the diffusion process is classified as a first-order process.

The advection and dispersion processes form the basis of crude oil contaminant modeling in groundwater contamination. Over time, there has been concern about the accuracy of the advection and dispersion equation used in contaminant site modeling because this approach depicts an ideal subsurface condition that rarely pictures the true subsurface field heterogeneity. ) suggested that solute transport time and pore-scale geometry control the magnitude of the effect of dispersion; thus, it is assumed that advection-induced flow is homogenous most times, leading to an underestimation of the dispersion coefficient. This was further confirmed by ), who revealed that at meter-scale heterogeneity, effective dispersion is more sensitive to pore-scale heterogeneity. Solute mixing by hydrodynamic dispersion reflects the heterogeneity of the aquifer at a scale smaller than that associated with the measurement of advection. Thus, the advection–dispersion process is believed to likely overestimate changes in concentration as a dispersion-based attenuation process which sometimes does not accurately simulate release from matrix diffusion storage. The diffusion-dominated process, compared to the dispersion-dominated process, has become essential due to the heterogeneity and anisotropy concept governing crude oil transport (). A new development in the groundwater contamination modeling approach is trending toward a heterogeneity-dominated advection–diffusion approach (). The interaction of the aforementioned processes results in an infinite number of scenarios that control how crude oil contaminants migrate in the subsurface. This poses a significant challenge for environmental managers in identifying primary crude oil controlling factors and considering them in evaluating and optimizing the associated risk (Swartjes, 2011). In addition, the interplay of soil, water, and contaminant poses a substantial challenge for researchers in forecasting important parameters that aid in understanding subsurface dynamics and contamination evolution (). Depending on the scale of the process, the physical, chemical, and biological parameters can be determined through laboratory and field research (; ).

Until recently, no studies summarized the state of contamination of freshwater and coastal aquifers in Nigeria. Crude oil contamination of soil and groundwater in Nigeria is widespread, with little information about the quantity trapped within the unsaturated zone and the groundwater system, making it difficult to understand the critical process that controls its migration. There is expected to be a paradigm shift in future contaminant research in Nigeria. The current understanding of the behavior of crude oil contaminants is largely deficient for properly optimizing remediation strategies (). This new shift should be data-driven () and involve detailed laboratory experiments simulating the subsurface site conditions from which different processes can be investigated before further upscaling to field-scale research. The initial approach of PINOT sampling investigation needs to be replaced by long-term monitoring of periodic changes in contaminant concentrations, water tables, and biological facies. These will serve as a base for further investigation of knowledge gaps relating to uncontrolled groundwater extraction, climate change, and sea level rise.

3 Approaches to investigate the fate and transport of crude oil contaminants

To better understand the processes in terms of their individual and combined influence on the fate and transport of crude oil pollutants, multiphase and reactive transport processes must be investigated at both the experimental and field scales. The combination of these two approaches has led to the development of predictive models that allow for subsurface modeling of contaminants’ evolution, monitoring their movement, and predicting their final fate (). These long-term estimates help evaluate natural attenuation or cleanup operations and determine the best reaction in the event of a spill mishap. In addition, data from quantitative and qualitative field research are regularly used in the calibration and validation of prediction models, serving as a guide for more effective and efficient remediation implementation. A summary of the selected literature on experimental and field-based approaches is summarized in Table 1 and discussed in the following sections.

TABLE 1

AuthorTopicGeology/source typeAimMethodologyResult(s)Conclusion(s)Contribution to knowledge
Van Geel and Sykes (1994)Laboratory and model simulations of a LNAPL spill in a variably-saturated sand, Laboratory experiment and image analysis techniquesGeology: well-sorted sandTo simulate LNAPL spill using LNAPL pressure–saturation distribution and image analysisEBD: 150/120/6.0 cm• The increase in water pressure did not correspond to a sharp increase in water saturation, as indicated by the static water–air capillary pressure–saturation relationship• The deviation from static conditions could also be partially attributed to the interfacial spreading of the LNAPL ahead of the LNAPL front• The image analysis techniques offer an excellent non-destructive tool for determining a NAPL saturation profile during a laboratory experiment
Source: n-heptaneEBM: stainless steel/plate glass• The capillary pressure–saturation relationship appeared to be dynamic prior to the arrival of the LNAPL front and the impending significant increase in water pressure
Hydrogeology: grain size distribution, pressure data, inflow, and outflow rates
Image visualization
Fingering of dense non-aqueous phase liquid in porous media: Experimental investigationGeology: sandTo investigate the propagation of dense non-aqueous phase liquids (DNAPLs) in water-saturated homogenous porous mediaEBD: 45/450/120 cm• The mobility of DNAPLs is a function of the capillary pressure of the pore spaces• Gravitational, viscous, and capillary forces were observed to have a varied influence on the flow behavior• Identify the sensitivity of fingered patterns and mobility of DNAPLs to the grain size variation, contaminant viscosity, and spill volume
Source: trichloroethylene, trichloroethane, and dibutyl phthalate (0.5–1 l each)EBM: stainless steel/plate glass• The density property of DNAPLs reduces the adequate pressure at the front end of the mobile contaminant• The ganglia of DNAPLs will be trapped whenever the displacement front passes and the water table fluctuation will not remobilize the entrapped DNAPLs
Hydrogeology: grain size distribution• At a microscopic scale, the heterogeneity of the porous media will create redistribution due to hydrodynamic dispersion
Image visualization• At a macroscopic scale, subsurface heterogeneity restricts the development of DNAPL fingering
Experimental study of movement and distribution of dense organic contaminants in heterogeneous aquifersGeology: sandTo observe the influence of aquifer heterogeneities on the movement and distribution of immiscible contaminantsEBD: 122/5/183 cm• The saturation measurements made with a dual-gamma system were susceptible to the path length of the gamma beam and hence to the internal thickness of the plume• Initial water saturation significantly impacts the transport of the fine soil’s organic capillary forces that entrap most of the organic contaminants• They demonstrated the influence of grain size distribution, water saturation, viscosity, and capillary forces at the pore scale on the migration of DNAPLs
Source: dibutyl phthalate and 1,1,1-trichloroethaneEBM: transparent stainless steel and plexiglass• The presence of coarse lenses did not influence the transport of the organic contaminant• In porous media, the amount of spill is controlled by capillary forces, and the viscosity of the organic material controls its movement
Hydro geochemical: grain size and saturation distribution• The fine lens altered the migration and entrapped the organic contaminant.
Image visualization• The impediment of the organic material in the fine lens leads to the horizontal accumulation of the organic material along the base of the coarse layer.
Measurement of fluid contents by light transmission in transient three-phase oil-water-air systems in sandGeology: sandTo develop a method by which fluid content can be measured rapidly in three-phase systemsEBD: 45/1/55 cm• Packing density differences produced minor variations in the intensity of the transmitted light• Visualizing the hue and intensity of the images allows the qualitative and quantitative characterization of water and liquid contents as they spread within the system• They presented the applicability of the light transmission method in estimating fluid contents in transient air–oil–water flow in sandy porous media
Source: oilEBM: thick poly-carbonate walls• The RGB color image could clearly distinguish the various oil, air, and water contents
Image visualization: Synchrotron X-rays light transmission method• Both two- and three-phase systems exhibit typical, unstable flow behavior
• Unstable finger behavior in three-phase systems exhibits many similarities to that in two-phase systems
Evidence for microbial enhanced electrical conductivity in crude oil-contaminated sedimentsGeology: fine to medium-grainedTo investigate changes in electrical conductivity during microbial mineralization of crude oils in a mesoscale column experiment that partially simulates field conditionsEBD: 80 cm• Higher numbers of oil-degrading microbial population shifts in microbial community structure, and elevated concentrations of cations occur within the contaminated column• The mechanism for the high bulk conductivity value observed in crude oil-contaminated sediments is related to enhance mineral weathering from metabolic by-products due to microbial degradation of the crude oil• They revealed the ability of geophysical methodologies to investigate microbial activity and geomicrobiological processes
Source: dieselEBM: polyvinyl chloride pipe geochemical analysis• The residual and diesel-saturated sediment layers were the most biologically active and the most bio-geochemically altered
Microbial analysis
Geophysics
Two-dimensional laboratory simulation of LNAPL infiltration and redistribution in the vadose zoneGeology: sandTo investigate the immiscible flow of LNAPLs in the unsaturated zone under homogeneous soil conditionsEBD: 120/8/180 cm• At the initial stage, LNAPL spread is majorly driven by gravity and, eventually, capillary pressure as it moves through the porous medium. Therefore, scaling two-phase P–S functions from their water relation and predicting three-phase P–S relations from these scaled two-phase functions, models should include the dependence of the residual wetting fluid saturation on the fluid pair under consideration• The knowledge of the distribution of residual NAPLs is essential as it can act as a constant long-term source of contamination with NAPL vapors migrating through the vadose zone and with NAPL components partitioning into precipitation water• They show the importance of residual and wetting fluid saturation in predicting three-phase P–S relations from scaled 2-phase P–S functions
Source: Soltrol 220 (isoparaffinic solvent)EBM: Perspex/glass
Hydro geochemical: grain size distribution, saturation distribution, water, and LNAPL pressure measurement
Image visualization
Estimation of contaminant transport parameters for a tropical sand in a sand tank modelGeology: poorly graded sandTo fabricate a 3D intermediate-scale sand tank model, characterize the porous medium material and the flowEBD: 180/30/80 cm• Longitudinal dispersion is dominant in the direction of flow• The characterized sand tank model is intended to serve as an economical and versatile tool for simulating and understanding field conditions of soil and groundwater contamination for similar tropical aquifer materials• They estimate the quantitative advection and dispersion parameters for further use in the contaminant transport model
Source: saltEBM: plexiglass• Estimated average seepage velocity values at outflow ports 1 and 2 are v = 4.13 × 10−6 m/s and v = 1.14 × 10−5 m/s, respectively.
Tracer test• The dispersion coefficients and dispersivity from the experiment represent the porous medium hydrodynamic parameters obtained for sandy aquifer materials
Non-aqueous phase liquids distribution in three-fluid phase systems in double-porosity soil media: Experimental investigation using image analysisGeology: sand and kaolin clayTo study the migration of both LNAPLs and DNAPLs in double-porosity soil in a 2D flow chamber in three fluid-phase systems (air–NAPLS–water) using the light transmission visualizationEBD: 45/30/1 cm• DNAPLs moved horizontally during the initial migration minutes in double-porosity soil before continuing its migration vertically downward to penetrate the water table• Pore size and geometry significantly influence the media’s porosity, thereby controlling the migration velocity and crude oil contamination• They revealed that light transmission visualization provides a non-intrusive and non-destructive technique for studying multiphase flow in double-porosity soil media where rapid changes in fluid distribution in the entire flow domain are not easy to measure using conventional tools
Source: toluene and tetra-chloroethyleneEBM: acrylic material• NAPL migration rates in the double-porosity medium were much faster than those of the single-porosity medium.• The chemical properties of NAPL have a significant influence on its migration in porous media
Light transmission visualization• Inter-aggregate pores in the double-porosity soil contribute to the increasing velocity of fluid migration through porous media• Toluene’s migration velocity was much faster than tetra-chloroethylene’s migration velocity
• The wettability of fluids and capillary pressure characteristics in the soil pores were influential factors in fluid migration within porous media
Sualla et al. (2018)Geo-electric Laboratory Simulation for Characterization of Crude oil Impacted Coastal SandsGeology: Coastal Sand and clayTo simulate Crude oil spills on coastal sands in a laboratory model box, a characterization of spills using the electrical methodEBD: 150/70/50 cm• The high resistive anomaly reflected the presence of the crude oil plume within the coastal sand medium.• The resistivity signatures obtained revealed that the migration of the Crude oil plume was more prominently associated with the hydraulic distribution of the porous media.• They showed the application of geophysics in Crude oil spill delineation in a coastal sand environment.
Source: Crude oilEBM: Plank• The magnitude of the resistive anomaly was slightly reduced due to the influence of the saline water.
• The crude oil response was imaged as a high apparent resistivity anomaly in the partially saturated soil.
Assessment of LNAPL in Subsurface under Fluctuating Groundwater Table using 2D Sand Tank ExperimentsGeology: clean sandTo investigate the fate and transport of toluene, a light non-aqueous phase liquids (LNAPLs), in the subsurface under dynamic groundwater table conditionsEBD: 150/10/120 cm• The large LNAPL pool area gets developed under rapidly fluctuating groundwater conditions, which significantly enhances the dissolution rate and contributes to a high concentration of dissolved LNAPL to the receiving receptors• The transport of dissolved LNAPL plumes is comparatively fast in the case of rapidly fluctuating groundwater, resulting in closely spaced concentration isolines of the toluene-containing plume• They revealed the importance of understanding the effect of varying subsurface water level conditions on the fate and transport of LNAPLs
Source: salt/tolueneEBM: stainless steel• Dissolution is highly affected by the groundwater table, which may cause loss of contaminant mass as a dissolved phase• The biodegradation rate is directly related to the volume of the dissolved plume
Tracer test• The microbial growth increased as the plume moved from the LNAPL pool toward the low gradient
Microbial analysis• Lowering the groundwater table increases air-filled porosity in the smear zone, influencing vapor phase LNAPL partitioning
Numerical modeling

Summary of selected experimental research on the fate and transport of crude oil.

3.1 Experimental approach

An experimental technique has been applied in a controlled environment under regulated conditions to improve the understanding of the dynamics of crude oil contaminants in both the saturated and unsaturated zones. Because of physical, chemical, and biological heterogeneity, field-scale research is typically complicated, but this complexity is decreased in an experimental context, allowing for high-resolution, microscopic studies. Multiphase studies have been carried out in sediment-packed one-dimensional (1D) columns, two-dimensional (2D) slabs, and three-dimensional (3D) tanks. carried out 2D multiphase flow experiments in variably saturated sand and unsaturated sand with varying porosities, with an emphasis on the effect of the capillary pressure–saturation relationship, the sensitivity of finger pattern development to grain size distribution, spill condition, and contaminant properties on the flow dynamics. At the capillary fringe, a reduction in vertical migration at the capillary fringe was observed as the presence of crude oil contaminants led to the depression of the water table. Also, the injection rate plays a vital role in the magnitude of vertical and horizontal migration under anisotropic hydraulic conductivity conditions; a lower injection rate promotes vertical migration, and a higher rate promotes horizontal migration. further carried out a quantitative 2D laboratory experiment to investigate the immiscible flow of crude oil in the vadose zone under transient flow. Their results demonstrate the adequacy of constructive pressure–saturation relations used in multiphase flow models. They proposed that hydraulic relations for three-fluid-phase models should account for residual NAPL saturation occurring after NAPL drainage as residual NAPL can act as a non-negligible constant source of contamination. ) identified the effect of using a single-rate mass transfer to describe a more extended timescale flow system because the effective dispersivity may increase at a longer scale system, resulting in the underestimation of contaminant spread within the pores.

Experimental characterization of crude oil contaminant entrapment and movement within unsaturated pores has shown that the dispersion process can cause the development of fingers due to pore-scale heterogeneity and fluid interactions (). Furthermore, residual waters are typically trapped as isolated macroscopic clusters as a result of crude oil displacement of water without buoyancy in water-wet heterogeneous sediments, with the degree of entrapment depending on the heterogeneity structure and crude oil penetration rate ().

3.2 Field-based approach

Field-based measurement is required to obtain in situ subsurface information at contaminated sites. These approaches include, but are not limited to, core sample measurement, partitioning tracer techniques, geophysics measurement, and geochemical analysis. Several studies have provided a quantitative appraisal of the advective property of hydrocarbon contaminants estimated from tracer tests, revealing additional information about the dominant process that controls the migration of contaminants in the subsurface. Although discovered that non-equilibrium partitioning could enhance results for high-saturation zones, partitioning tracer studies are mostly utilized for crude oil residual saturation estimation. Using a geochemical method, described the biodegradation process of crude oil in detail. Their research revealed that aerobic conditions result in a higher rate of degradation than anaerobic settings, and the presence of chlorine atoms enhances biodegradation. Changes in alkalinity, microbial populations, and preferential elimination of BTEX components along the transport path can be used to map subsurface attenuation. discovered that tidal oscillations in a hydraulically connected water body add complexity to subsurface contaminant models. Also, capillary fringe affects the fraction of benzene biodegraded in saturated and unsaturated zones (). observed a complex distribution of crude oil in the subsurface oil, whose saturation distribution has been influenced by sediment heterogeneities and water table fluctuation. Hydrogeologic studies indicated that the crude oil migrated laterally within the capillary fringe, and the center of the oil body depressed the water-saturated zone boundary.

Due to the spatial challenges involved with contaminant movement within the subsurface, in situ field hydrogeological and biogeochemical approaches are limited in capturing the spatial and temporal variation in properties and processes of interest (; ). Tomographic methods, including hydraulic tomography () and tracer tomography (), have been developed in hydrogeology to improve the spatial resolution of flow and transport properties. However, the resolved spatial scale is yet limited to a few meters-scale. Non-destructive geophysical techniques provide high-resolution spatial information about the subsurface condition and the fate of contaminants. Geophysical approaches have been demonstrated to be effective in contaminated site characterization and monitoring as they provide quantitative information on the contaminant flow process in both saturated and unsaturated media (). In addition, the geophysical technique is sensitive enough to detect biogeochemical changes caused by degradation and redox reactions, making it useful for environmental monitoring (). The use of electrical resistivity tomography to determine the spatial distribution of hydrocarbon pollution has received much attention (). Electrical resistivity can be utilized to discover active sources of a dissolved plume by targeting minimal abnormal resistivity levels, identifying changes to the soil matrix, and determining the presence of a dissolved plume (). This information can be utilized to create geoelectric models that can be compared with hydrocarbon concentration measurements for high-sensitivity interpretation (). However, due to the direct relationship between electrical resistivity, water saturation, porosity, and soil texture, the sole usage of the electrical resistivity method in measuring the degree of biodegradation creates either an overestimation or underestimation (; ). According to , for electrical resistivity to image the pollutant distribution, the apparent bulk electrical conductivity of the saturated zone plume must be three to five times higher than the background values. As biodegradation occurs, contaminated zones become more conductive and may be masked by conductive soil ().

Furthermore, microbial adaptation to petroleum as a carbon source can significantly impact the surrounding medium’s electrical signature (). Changes in the vertical organization of microbial communities may also result from increased geoelectric measurements. This could guide microbial ecology sampling during natural bioremediation procedures and monitoring (). Bacterial biodegradation of hydrocarbons alters groundwater redox conditions, forming zones with varying redox potentials. This results in the formation of redox gradients, which generates an electrochemical source for the SP signal via natural geo-battery processes (; ). Significant negative SP anomalies have been recorded at hydrocarbon-contaminated sites, and a linear association between SP signals and redox potential measurements has been established. identified a rapid shift in SP signatures from negative values in polluted zones to positive values in clean zones, allowing the contaminated and clean zones to be distinguished and reflecting the extent of alteration at the contaminated site. This method provided an indirect indicator of bacterial activity and the presence of biofilms and mineral precipitates, which are likely to have played a critical role in the electrical current conduction of the plume and the oxidizing conditions outside the plume. identified an IP response related to intrinsic biodegradation in contaminated areas. This response was caused by the interfacial properties of microbial and mineral fluids, which caused significant electrical impulses in the polluted zones. With IP, valuable information on the mineralization of hydrocarbon pollutants as a result of contaminant degradation and dissolution can be obtained. These data could show the buildup of metallic minerals in hydrocarbon-contaminated soil (). Changes in groundwater chemistry observed from spectral-induced polarization data at an ongoing remediation site after acetate input and withdrawal can be utilized to find ideal conditions for contaminant extraction (). Because IP is sensitive to geochemical changes, the phase response can be used to understand spatial heterogeneity caused by pore size variability and free phase contaminant distribution ().

3.3 Numerical modeling

Multiphase models have proven to help model the fate and transport of crude oil pollutants in soil and groundwater. Several scholars have utilized various numerical methods to solve the multiphase flow. These numerical models were developed to replicate specific environmental processes in the laboratory and field under specified environmental conditions. This approach provides a better understanding of the controlling factors influencing the fate and transport of crude oil contaminants (). Such factors include the influence of subsurface heterogeneity and capillary pressure–saturation–relative permeability (Pc-S-kr) relationships. The numerical models successfully reproduced the general features of experimental data sets when appropriate representations of heterogeneity and Pc-S-kr functions were used (). However, the dearth of detailed multiphase field data in different hydrogeological settings has hindered the ability to conduct comprehensive field-scale model testing. The unpredictability of the source-sink scenario in the unsaturated zone can be attributed to the combined actions of dissolution, biodegradation, and volatilization. The coupled system of moisture flow, airflow, and corresponding transport equations modeling demonstrated a gradual decrease in source intensity, corresponding to a decrease in the mass transfer process (). The simulation of a two-phase steady-state flow in the capillary force-dominated regime revealed a low sensitivity of the interfacial area, capillary pressure, and relative permeability to the initial fluid distribution, ruling out the use of capillary pressure measured in a displacement experiment for steady-state infiltration flow analysis (). ) confirmed this by observing specific differences in the relative permeability prediction and attributing them to inlet pore variation, which regulates the dispersion of pollutants. Multispecies reactive transport modeling results are influenced by variable degrees of heterogeneity and associated uncertainty in the aquifer’s physical parameters ().

Several numerical models have been developed to understand biodegradation and natural attenuation better. For example, simulated the movement of dissolved crude oils in the 2D domain under the influence of oxygen-limited biodegradation, reaeration, and anaerobic biodegradation. This study discovered that aerobic degradation controls most of the remediation process, with stripping and volatilization accounting for a minor portion of the process. ) analyzed crude oil pollutant attenuation in groundwater utilizing first-order solute decay due to advection, dispersion, sorption, and biodegradation processes under aerobic and anaerobic circumstances. ) used estimated flow across the vadose zone to demonstrate the effects of LNAPLs on water table aquifers. Weaver (1996) and Zheng et al. (2010) investigated groundwater seepage as a function of time, precipitation infiltration, and excess pore-water pressure in the saturated and unsaturated zones. Siracusa et al. (2007) investigated subsurface pollutant reactions, microbial metabolisms, and microbial transport kinetics. Because multiphase flow models are computationally expensive, alternatives to continuum-based modeling have been developed, including analytic solutions, invasion-percolation modeling, and dissipative particle dynamics modeling methodologies. This model has been used to simulate steady-state or transient flow and advective and dispersive transport of dissolved oxygen and crude oil plumes.

4 Crude oil contaminant fate and transport implications of soil and groundwater resources

Most crude oil contaminant studies focus on developing an effective remediation plan while considering subsurface heterogeneity and the pollutant’s physical and chemical properties. A good cleanup requires a comprehensive investigation of the polluted location. However, evaluating so many variables is complex and time-consuming. Understanding baseline subsurface conditions before contamination, subsurface hydraulic variation in relation to contaminant distribution, and the impact of external influences such as groundwater fluctuation, climate change, and groundwater extraction are all critical components of crude oil contaminant research (). In the control of crude oil environmental contaminants, future event mitigation and site rehabilitation are commonly used. Given the harmful effects of crude oil pollutants on the environment and human health, developing an efficient crude oil spill mitigation strategy and responding swiftly are vital. According to , crude oil contamination and soil/aquifer vulnerability are linked to bulk movement advection, and crude oil contamination spreads in saturated and unsaturated zones. In saturated and unsaturated zones, crude oil contamination is linked to groundwater level, net recharge, and hydraulic conductivity on climate change; changes in groundwater levels and increasing sea levels are highlighted since they are critical topics of contaminant research and vital for managing contaminants.

Groundwater accounts for the vast majority of the world’s fresh water. Depending on the recharge rate and underlying quality, groundwater can be found at depths ranging from 1 m in wetlands and marshy regions to several meters in other habitats (Winter, 1999). Groundwater, unlike soil, has a significant degree of mobility. As a result, the contaminant can be transported from one site to another via groundwater. Contaminants can travel from the upper soil layers to the lower soil layers, eventually reaching the groundwater zone and surface waters (). Because water travels considerably more slowly in aquifers than in surface water systems, it might take years, decades, or millennia for infiltrated water and pollutants to transit from aquifers to the point of discharge (Xia et al., 2010). This duration can vary according to the kind of pollution, aquifer characteristics, and flow path length. As a result, there is ample time for close interactions between the water and the soil and rock material, resulting in mineral composition changes in the water. Because of the significant variations in the geometry of hydrological systems and the soil and rock materials, the movement of contaminants from the source to sink has a varying impact on groundwater quality (natural background concentration) (). Naturally, disturbed groundwater will occasionally pose intolerable risks to human health or the environment (). This is demonstrated by the arsenic poisoning of many people in Bangladesh and West Bengal, India, caused by a natural phenomenon rather than industrial contamination (). In the event of immobile contaminants in the upper soil layers, leaching into groundwater could last for decades. For example, identified dieldrin in quantities exceeding the drinking water requirements in several wells at the old US military installation Clark Air Base in the Philippines, 16 years after dieldrin was last used there.

5 Scientific challenges in fate and transport of crude oil contaminant research

In the early field investigations of crude oil contaminant flow, it was discovered that the subsurface was heterogeneous due to a variety of variable factors: capillary pressures have been shown to reduce crude oil contaminant relative permeability and prevent infiltration into fine-grained layers, and hydraulic gradients cause down-dip migration, spreading, pooling, and trapping of crude oil contaminants by these fine-grained (frequently water-saturated) layers (). In addition, it was discovered through experiments with a multiphase flow that fingers occur in crude oil contamination infiltration due to pore-scale and macro-scale variability. All of these factors point to the difficulty in pinpointing crude oil pollution source zones (Swartjes et al., 2012). Many decades have passed without an answer to whether contaminant transport models can be trusted. The reliability of groundwater flow and pollutant transport models is a tough challenge, even in homogenous porous media. It is vital to separate reality from it in order to employ models.

5.1 Hydraulic heterogeneity

The diversity of hydrogeological conditions and changing crude oil chemistry are expected to make field deployment of these systems difficult. The subsurface geometry must be considered in terms of vast amounts of subsoil because plumes frequently traverse soil volumes ranging from several hundred to thousands of cubic meters in a year (Swartjes and Grima, 2011). Even though these enormous volumes are typically formed of numerous irregular layers of different materials and commonly include entities of distinct subsurface material, the employment of contaminant transport models requires the production of a more or less regular-looking profile (Vallero, 2021). The typical, irregular, geochemical input parameters for massive aquifer material must be established. As a result, the description of both hydraulic and geochemical input parameters contributes to uncertainty in expected groundwater quality (). It is generally recognized that the availability of sufficient input parameters is crucial to the accuracy of model computations in any other model. When attempting to provide a complete picture of heterogeneities, it is vital to consider the consequences of interactions with hydraulically related aquifers and surface water. Because of the high concentration of saltwater in coastal aquifers with a freshwater–saltwater boundary, seawater intrusion becomes a dynamic and three-dimensional phenomenon that causes water quality oscillations on horizontal and vertical scale interactions. Also, contact angle hysteresis, wettability, complex pore geometry, and field-scale heterogeneity may arise, thus adding in another degree of complexity (). Despite advances in image processing and micro-modeling approaches, turning experimental results into data suitable for evaluating macroscopic, continuum-based multiphase flow models remains challenging, especially for large-scale studies. Multiphase modeling at the pore scale has also revealed how pore-scale processes affect macroscale constitutive relations (). It is essential to assume a precise hydraulic gradient over a relatively broad area to be accurate in hydraulic situations.

5.2 Spatial and temporal heterogeneity

The concepts of homogeneity and heterogeneity in a system are inextricably linked to the scale at which the investigation is carried out. On a large scale, the same groundwater body can be termed homogeneous, whereas on a small scale, it can be deemed heterogeneous (). As a result, the study’s scope must be addressed from the start of any project. This scale contrasts between adjacent clean groundwater and public water supply groundwater sources that are only a few kilometers away when it comes to groundwater. For example, when building a risk management solution, the optimum size is determined by the type of solution and changes based on whether the solution is for source elimination or pathway interception (). Field sites must be characterized in various hydrogeological situations to evaluate models that predict the fate of crude oil contaminants. Steefel et al. (2005) described the impact of pore-scale processes on the macroscopic constitutive relationships driving crude oil contaminant saturation, relative porosity (permeability), entrapment, and mass transfer rate. In recent years, non-invasive approaches for determining the source zone architecture of crude oil pollutants have made significant progress.

Due to their complexity, it is still difficult to estimate field concentration distributions and temporal oscillations in crude oil contamination sources. A multi-constituent crude oil contaminant’s chemical composition changes over time as the constituents dissolve and volatilize at varying rates (). As a result, the effluent concentration histories for components dissolved in NAPLs differ. The breakdown of crude oil pollutants and changes in source architecture may result in a nonmonotonic drop in effluent concentration history. It has been proven in laboratory trials that heterogeneity, NAPL saturation, and interaction with moving water affect dissolution. In the presence of a solvent, the degradation of dissolved solutes can also affect the concentration gradient that promotes interphase mass transfer, thus creating a new complexity. However, recently developed hybrid models that dynamically combine pore-scale models for determining constitutive interactions with macroscale flow and transport show promise, but they face several computational and practical hurdles and are limited in their use. The key focus of future research should be on creating techniques for upscaling constitutive relations for field-scale studies.

5.3 Current state of crude oil contamination studies in Nigeria

Crude oil resources continue to be an important and growing source of revenue in Nigeria. A total of 13 multinational firms operate 606 oil fields (355 onshore and 251 offshore), 5,284 oil wells, and 7,000 km of oil and gas pipelines in the country, with the Niger Delta region serving as the epicenter of its production and processing operations (). The Niger Delta covers an area of roughly 20,000 square kilometers in southern Nigeria and accounts for approximately 8% of the country’s total land area. It is one of the world’s largest wetland and mangrove forest ecosystems. Sabotage, pipeline vandalism, well blowouts, and engineering failure have caused several crude oil spill accidents in Nigeria. Since the start of the Nigerian oil industry in 1960, 13 million tonnes of crude oil has been reported to have spilled in the Niger Delta due to terrorism, pipeline vandalism, and engineering failure (Figures 2A–C; Figure 3). The Nigerian Oil Spill Detection and Response Agency (NOSDRA) recently revealed statistics indicating that there have been a total of 12,566 confirmed oil spill occurrences from 2006 to the present, resulting in a total of 694,890 barrels of crude oil spilled. Crude oil contamination accounts for 81% of all-known oil leak occurrences (Figure 2D). Not only have agricultural resources in this region been poorly impacted due to crude oil contamination, but also surface water and groundwater resources have been rendered unfit for use in affected areas. The impact of the pipeline rupture that occurred between 1994 and 1996 in Baruwa, Lagos, is still noticeable, according to Uduebor and Ola (2016) and , as a hydrocarbon thickness of 0.72 m was still discovered. Despite the contamination level being known, the reliance on unmanaged natural attenuation as a remediation strategy is grossly insufficient. The contamination still poses a significant risk to the region’s soil and groundwater resources, 20 years after the reported crude oil spill. Several studies on the environmental impact of oil spills in the Niger Delta and other tropical places worldwide have repeatedly revealed that areas immediately exposed to large or recurring oil spills or leaks frequently face long-term environmental challenges. Though there has been a reduction in the occurrence of oil spillage (Figure 4), the long-term effect of the previous spillage has created a continuous effect on land, soil, and the ecosystem at large. Even though the region’s potential for long-term prosperity has to be fully realized yet, its future is jeopardized by a slew of environmental issues, the most serious of which is oil contamination.

FIGURE 2

FIGURE 3

FIGURE 4

).

The influence of crude oil contamination on surface water and aquatic ecosystems in Nigeria, as well as its agriculture, biological diversity, human health, and groundwater sources, has been studied (; Uko et al., 2020). With limited attention paid to the influence on groundwater contamination, the majority of this research concentrated on understanding the advection and dispersion processes of crude oil pollution on surface water (). Even yet, dissolved sources of contamination from surface water eventually have an impact on the groundwater system. Most of the investigation that has been carried out on subsurface crude oil contamination uses a point sample approach to assess the total crude oil content in soil and water samples. who measured the residual concentration values of soils and groundwater within the Ijegun oil spills in 2008 between 2009 and 2010 could only confirm the amount of contamination; no information was supplied as connected to the contaminant’s migratory pattern. The unregulated release of crude oil pollutants has increased the number of heavy metals in soils, plants, and groundwater, directly impacting human health and wellbeing (). Following a spill, the subsurface properties are likely to change due to the biodegradation process caused by microbial activity. As a result, knowledge of these alterations is crucial in developing a site-specific remediation strategy. The porosity and permeability of impacted soils have been reported to diminish due to the establishment of free-phase crude oil adhesion, generating a ganglion along the pore as a result of the replacement of water molecules by crude oils (). This has the potential to be very important in shaping migratory patterns throughout the system, among other things.

Furthermore, as demonstrated by various researchers on the high degree of pollution, comprehending the distribution paths of contaminants in the subsurface has not been well investigated. According to , the parameters influencing the hazard posed by hazardous substances include their concentrations and distribution, the channels through which they are exposed, and the length of time they are exposed. The flow path of crude oil contamination in the Niger Delta region remains unclear. As of 2020, two scholarly publications have been published exploring the advection and dispersion processes of pollutant migration within the Niger Delta region. According to and Ugbena et al. (2020), the groundwater flow velocity is 105 m per day and 121 m per day, respectively. This shows how easily toxins can spread across a large area, harming a more significant number of groundwater supplies. However, the scope of this investigation is constrained by the assumption that the salt tracer’s flow velocity is similar to the flow velocity of crude oil plumes. This non-conformity is thought to be caused by the density and viscosity of the crude oil pollution in the subsoil. Due to a lack of control information, it was difficult to ascertain the baseline subsurface condition and hydraulic distribution before the spill, making it impossible to determine how contamination affected the physical and flow parameters. This knowledge will help carry out a more effective remediation exercise in the real world. It is expected that future research focusing on flow path characterization of contaminant transport will be carried out to provide necessary information for futuristic prediction and identification of a suitable remediation approach.

Attempts have been made to clean up hazardous places. who conducted the first systematic investigation of oil contamination in Ogoniland discovered that, despite the natural conditions for crude oil degradation being favorable due to the region’s high temperature and rainfall, a low recovery rate had been recorded. The most often employed remediation approach, Remediation by Enhanced Natural Attenuation (RENA) (Figure 5), has been inefficient since soil pollution persists even after sites have been certified as remediated. detected residual contamination at a site that had undergone remediation with a 3.6-m smear zone that posed a significant danger to the environment. The applicability of the RENA technique has been demonstrated to be insufficient due to the level of pollution on both the vertical and lateral sides (Sam et al., 2017). The pollution level in the Nigeria Delta region has been estimated to be higher than 5 m, and leachates have been identified in groundwater, raising questions about the effectiveness of the RENA strategy’s tilling activity.

FIGURE 5

5.4 Management responses to crude oil contamination in Nigeria

The Niger Delta has been identified as one of the world’s most heavily oil-impacted places. The Niger Delta has been labeled “one of the world’s most highly oil-impacted areas” due to oil exploration activities in the region for more than 5 decades (Zabbey and Uyi, 2014). The Restoration by Enhanced Natural Attenuation (RENA) strategy is Nigeria’s popular contaminated land remediation technology. However, RENA is an unsuccessful technique for handling the problem due to the volume and extent of the spills, with oil entering the soil to depths of more than 5 m and leaking into groundwater aquifers (). The most prominent assessment of the level of contamination in the region is the United Nations Environment Programme (UNEP) report, (Environmental Assessment of Ogoniland) (UNEP, 2011). One of this study’s most significant findings was the need to develop and deploy oil pollution mitigation measures. Pollution policies are inadequate in Nigeria, which has been connected to several concerns, including a fragmented government structure, a lack of transparency in decision-making, and ineffective policy implementation, to name a few.

The Petroleum Act of 1969 established comprehensive petroleum regulations to prevent environmental pollution in diverse environmental media such as air, water, and soil. Other notable pieces of legislation were enacted, such as the Harmful Waste Act (1988) and the Environmental Impact Assessment Act (1992), but these were not designed to tackle the problem of polluted land management in the first place. Following an increase in oil production, oil pollution incidents (), and community protests, particularly in Ogoniland (), the Nigerian government responded by enacting the Environmental Guidelines and Standards for the Petroleum Industry (EGASPIN) (2002). The EGASPIN spans from 2002 to the present and acts as the regulatory underpinning for Nigeria’s current environmental mandate. Many of EGASPIN’s environmental quality standard criteria were adapted from other countries, especially the United States (mainly from the USA). This has been attributed to Nigeria’s lack of technical ability and expertise in designing such rules and procedures. Concerns have been raised concerning the lack of contextualization in the current guidance and the fact that land use, soil type, and total organic carbon circumstances in the United States differ from those in Nigeria. Due to this disparity, the suitability and effectiveness of the risk assessment and management suggestions may be jeopardized.

6 Conclusion

The current state of research and management strategies in Nigeria, as well as our understanding of the transport processes and modeling conditions that govern the fate of crude oil contaminants in unconsolidated sediments, has been reviewed. Individual processes were discovered to be reasonably understood, and the utility of field and experimental research in estimating model parameters was demonstrated. Despite more than 50 years of research, there are still questions about crude oil contamination. Crude oil contamination is still a hot topic, thanks to complicated field-scale subsurface heterogeneity fueling ongoing research. Most previous studies of processes controlling the fate of crude oil contaminants have relied on laboratory-scaled experiments or various forms of numerical and analytical simulation. Although this provides pore-scale information about individual processes and their interactions, its application to field conditions necessitates further consideration of scaling and field-scale heterogeneity. The accuracy of numerical models has thus become a problem, the solution to which is embedded in detailed subsurface monitoring involving the integration of hydraulics, geophysical, and hydrological information with geological controls and biochemical processes. Due to a lack of active contaminant research and remediation, Nigeria urgently requires a paradigm shift in crude oil contaminant management. In Nigeria, the current understanding of the behavior of crude oil contaminants is woefully inadequate for optimizing management and remediation approaches. According to this review, the most significant shortcoming of crude oil contaminant studies in Nigeria is the lack of intensive baseline and monitoring studies to understand the subsurface heterogeneity and hydraulic variation that control the fate and transport of crude oil contaminants in the subsurface. In addition, no detailed subsurface quantitative framework research was conducted across all known contaminant sites in the country. A long-term perspective on monitoring contaminant behavior is required for a successful remediation exercise. Given how contaminant management is carried out in Nigeria, a national guide to detailed contaminant site management must include the best remediation practices and improve monitoring and evaluation methods.

Statements

Author contributions

MA: conceptualization, journal review, framework of the manuscript, and first draft preparation. MO: conceptualization, supervision, and review and editing. KD: conceptualization, supervision, and review and editing.

Funding

This review paper was funded in part by the African Union.

Acknowledgments

The authors wish to thank the African Union for providing financial support (scholarship) through the Pan African University Initiative.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

  • 1

    AbbasM.JardaniA.MachourN.DupontJ. P. (2018). Geophysical and geochemical characterisation of a site impacted by hydrocarbon contamination undergoing biodegradation. Near Surf. Geophys.16 (2), 176192. 10.3997/1873-0604.2017061

  • 2

    Abdel AalG. Z. A.SlaterL. D.AtekwanaE. A. (2006). Induced-polarization measurements on unconsolidated sediments from a site of active hydrocarbon biodegradation. Geophysics71 (2), H13H24. 10.1190/1.2187760

  • 3

    AdipahS. (2018). Introduction of petroleum hydrocarbons contaminants and its human effects. J. Environ. Sci. Public Health3 (1). 10.26502/jesph.96120043

  • 4

    AhrenholzB.NiessnerJ.HelmigR.KrafczykM. (2011). Pore‐scale determination of parameters for macroscale modeling of evaporation processes in porous media. Water Resour. Res.47 (7), 125. 10.1029/2010wr009519

  • 5

    AlazaizaM. Y. D.NgienS. K.BobM. M.KamaruddinS. A.IshakW. M. F. (2018). Non-aqueous phase liquids distribution in three-fluid phase systems in double-porosity soil media: Experimental investigation using image analysis. Groundw. Sustain. Dev.7, 133142. 10.1016/j.gsd.2018.04.002

  • 6

    AlinnorI. J.NwachukwuM. A. (2013). Determination of total petroleum hydrocarbon in soil and groundwater samples in some communities in Rivers State, Nigeria. J. Environ. Chem. Ecotoxicol.5 (11), 292297.

  • 7

    AllenJ. P.AtekwanaE. A.AtekwanaE. A.DurisJ. W.WerkemaD. D.RossbachS. (2007). The microbial community structure in petroleum-contaminated sediments corresponds to geophysical signatures. Appl. Environ. Microbiol.73 (9), 28602870. 10.1128/AEM.01752-06

  • 8

    AlmutairiM. S. (2018). An assessment of remediation strategies for Kuwaiti oil lakes. Environ. Geotech.5 (6), 345355. 10.1680/jenge.16.00022

  • 9

    AminnajiM.RabbaniA.NiasarJ.BabaeiM. (2019). Effects of pore‐scale heterogeneity on macroscopic NAPL dissolution efficiency: A two‐scale numerical simulation study. Water Resour. Res.55 (11), 87798799. 10.1029/2019wr026035

  • 10

    AndersonM. P. (1989). Hydrogeologic facies models to delineate large-scale spatial trends in glacial and glaciofluvial sediments. Geol. Soc. Am. Bull.101 (4), 501511. 10.1130/0016-7606(1989)101<0501:hfmtdl>2.3.co;2

  • 11

    AndersonM. P. (1984). Movement of contaminants in groundwater: Groundwater transport-advection and dispersion. Groundw. Contam., 3745.

  • 12

    AndersonM. P.WoessnerW. W.HuntR. J. (2015). Applied groundwater modeling: Simulation of flow and advective transport. Massachusetts, United States: Academic Press.

  • 13

    AratoA.WehrerM.BiróB.GodioA. (2014). Integration of geophysical, geochemical and microbiological data for a comprehensive small-scale characterizationof an aged LNAPL-contaminated site. Environ. Sci. Pollut. Res.21 (15), 89488963. 10.1007/s11356-013-2171-2

  • 14

    Arrubarrena-MorenoM.Arango-GalvánC. (2013). Use of electrical resistivity tomography in the study of soil pollution caused by hydrocarbons: Case study in Puebla (México). Bol. La Soc. Geol. Mex.65 (2), 419426. 10.18268/BSGM2013v65n2a21

  • 15

    ArtiolaJ. F.BrusseauM. L.LiuN. (2019). “The role of environmental monitoring in pollution science,” in Environmental and pollution science. 3rd ed., Issue i (Amsterdam, Netherlands: Elsevier Inc). 10.1016/B978-0-12-814719-1.00010-0

  • 16

    AtekwanaE. A.AtekwanaE. A.WerkemaD. D.AllenJ. P.SmartL. A.DurisJ. W.et al (2004). Evidence for microbial enhanced electrical conductivity in hydrocarbon-contaminated sediments: Microbial enhanced conductivity. Geophys. Res. Lett.31 (23). 10.1029/2004GL021359

  • 17

    AzimiR.VaezihirA.LenhardR. J.HassanizadehS. M. (2020). Evaluation of LNAPL behavior in water table inter-fluctuate zone under groundwater drawdown condition. Water12 (9), 2337. 10.3390/w12092337

  • 18

    AzizR.Joekar-NiasarV.Martinez-FerrerP. (2018). Pore-scale insights into transport and mixing in steady-state two-phase flow in porous media. Int. J. Multiph. Flow109, 5162. 10.1016/j.ijmultiphaseflow.2018.07.006

  • 19

    BadejoO. T.NwiloP. C. (2004). Management of oil spill dispersal along the Nigerian coastal areas. Dept.of Surv. Geoinformatics, Univ. Lagos.

  • 20

    BasuS.YadavB. K.MathurS.GuptaP. K. (2020). In situ bioremediation of toluene‐polluted vadose zone: Integrated column and wetland study. CLEAN–Soil, Air, Water48 (5-6), 2000118. 10.1002/clen.202000118

  • 21

    BennettP. D.BrumbachB.FarmerT. W.IVP. L. F.HathewayA. W. (1999). Remedy selection for cleanup of uncontrolled waste sites. Pract. Periodical Hazard. Toxic, Radioact. Waste Manag.3 (1), 2334. 10.1061/(asce)1090-025x(1999)3:1(23)

  • 22

    BiswasB.QiF.BiswasJ. K.WijayawardenaA.KhanM. A. I.NaiduR. (2018). The fate of chemical pollutants with soil properties and processes in the climate change paradigm—a review. Soil Syst.2 (3), 51. 10.3390/soilsystems2030051

  • 23

    BodoT. (2019). Deep issues behind the crisis in the Niger delta region: The case of oil exploration in ogoniland, rivers state, Nigeria. Asian J. Geogr. Res.2 (1), 112. 10.9734/ajgr/2019/v2i130078

  • 24

    BrusseauM. L.ArtiolaJ. F. (2019). “Chapter 12 - chemical contaminants,” in Environmental and pollution science. 3rd ed. (Amsterdam, The Netherland: Elsevier Inc). 10.1016/B978-0-12-814719-1.00012-4

  • 25

    BrusseauM. L. (2019b). “Chapter 15 - subsurface pollution,” in Environmental and pollution science. 3rd ed. (Amsterdam, The Netherland: Elsevier Inc). 10.1016/B978-0-12-814719-1.00015-X

  • 26

    BrusseauM. L. (2019a). “Physical processes affecting contaminant transport and fate,” in Environmental and pollution science. 3rd ed. (Amsterdam, The Netherland: Elsevier Inc). 10.1016/B978-0-12-814719-1.00007-0

  • 27

    BrusseauM. L.WalkerD. B.FitzsimmonsK., (2019a). “Physical-chemical characteristics of water,” in Environmental and pollution science. 3rd ed. (Amsterdam, The Netherland: Elsevier Inc). 10.1016/B978-0-12-814719-1.00003-3

  • 28

    CerkirgeH. M.Al-RabehA. H.GunayN. (1992). Use of three generations of oil spill models during the Gulf War oil spills. Proc. 15th Arct. Mar. Oil Spill Program Tech. Seminar, 93105.

  • 29

    ChangS. W.ChungI. M.KimI. H.JooJ. C.MoonH. S. (2021). Importance of infiltration rates for fate and transport of benzene in high-tiered risk-based assessment considering Korean site-specific factors at contaminated sites. Water13 (24), 3646. 10.3390/w13243646

  • 30

    Che-AlotaV.AtekwanaE. A.AtekwanaE. A.SauckW. A.WerkemaD. D. (2009). Temporal geophysical signatures from contaminant-mass remediation. Geophysics74 (4), B113B123. 10.1190/1.3139769

  • 31

    ChenP.HubbardS. S.WilliamsK. H.Flores OrozcoA.KemnaA. (2012b). Estimating the spatiotemporal distribution of geochemical parameters associated with biostimulation using spectral induced polarization data and hierarchical Bayesian models. Water Resour. Res.48 (5). 10.1029/2011WR010992

  • 32

    ChenP.LiY.LanG.LiuB.ZhouH. (2012a). “Oil spills detection and monitoring using airborn thermal infrared remote sensing in dalian Xingang oil pipeline explosion,” in Proceedings of the 2012 2nd International Conference on Remote Sensing, Environment and Transportation Engineering, Nanjing, China, 01-03 June 2012. 10.1109/RSETE.2012.6260647

  • 33

    Chukwu OkeahG. O.Ushie OyegunL.UwadiaeC., (2019). Analysis of soil quality in a deltaic crude oil polluted environment Niger delta Nigeria. South Asian Res. J. Nat. Prod.,2(4).

  • 34

    CirpkaO. A.LevenC.SchwedeR.DoroK.BastianP.IppischO.et al (2014). “Tomographic methods in hydrogeology,” in Tomography of the earth’s crust: From geophysical sounding to real-time monitoring (Cham: Springer), 157176.

  • 35

    CocârţăD. M.StoianM. A.KarademirA. (2017). Crude oil contaminated sites: Evaluation by using risk assessment approach. Sustainability9 (8), 1365. 10.3390/su9081365

  • 36

    DagoO.JustinA.HarcourtP. (2018). Groundwater study for the identification of crude oil pollution sources within a hydrocarbon processing facility in Niger delta, Nigeria. Int. J. Appl. Geol. Geophys.6 (3). 10.9790/0990-0603035762

  • 37

    DarnaultC. J. G.DiCarloD. A.BautersT. W. J.JacobsonA. R.ThroopJ. A.MontemagnoC. D.et al (2001). Measurement of fluid contents by light transmission in transient three-phase oil-water-air systems in sand. Water Resour. Res.37 (7), 18591868. 10.1029/2000WR900380

  • 38

    Delgado-RodríguezO.ShevninV.Peinado-GuevaraH.de Guevara-TorresM. L. (2018). Characterization of Crude oil-contaminated sites based on geoelectrical methods of geophysical exploration. Geophysics, 8593.

  • 39

    DillardL. A.EssaidH. I.HerkelrathW. N. (1997). Multiphase flow modeling of a crude‐oil spill site with a bimodal permeability distribution. Water Resour. Res.33 (7), 16171632. 10.1029/97wr00857

  • 40

    DohertyV. F.OtitolojuA. A. (2013). Monitoring of soil and groundwater contamination following a pipeline explosion and petroleum product spillage in Ijegun, Lagos Nigeria. Environ. Monit. Assess.185, 41594170. 10.1007/s10661-012-2858-8

  • 41

    DongS.DaiZ.LiJ.ZhouW. (2018). The scale dependence of dispersivity in multi-facies heterogeneous formations. Carbonates Evaporites33 (1), 161165. 10.1007/s13146-018-0421-6

  • 42

    DornerS. M.AndersonW. B.SlawsonR. M.KouwenN.HuckP. M. (2006). Hydrologic modeling of pathogen fate and transport. Environ. Sci. Technol.40 (15), 47464753. 10.1021/es060426z

  • 43

    DoroK. O.CirpkaO. A.LevenC. (2015). Tracer tomography: Design concepts and field experiments using heat as a tracer. Groundwater53 (1), 139148. 10.1111/gwat.12299

  • 44

    DoroK. O.EhosiokeS.AizebeokhaiA. P. (2020). Sustainable soil and water resources management in Nigeria: The need for a data-driven policy approach. Sustainability; 12(10):4204. 10.3390/su12104204

  • 45

    DoroK. O.EmmanuelE. D.AdebayoM. B.BankC. G.WescottD. J.MickleburghH. L. (2022b). Time-lapse electrical resistivity tomography imaging of buried human remains in simulated mass and individual graves. Front. Environ. Sci.10, 501. 10.3389/fenvs.2022.882496

  • 46

    DoroK. O.KolapkarA. M.BeckerA. M. (2021). “Using shallow subsurface geophysical models to guide restoration of old agricultural fields in northwestern Ohio,” in SEG/AAPG/SEPM first international meeting for applied geoscience & energy (Texas, United States: AAPG).

  • 47

    DoroK. O.LevenC.CirpkaO. A. (2013). Delineating subsurface heterogeneity at a loop of River Steinlach using geophysical and hydrogeological methods." Environ. earth Sci.69 (2), 335348. 10.1007/s12665-013-2316-0

  • 48

    DoroK. O.StoikopoulosN. P.BankC. G.FerrisF. G. (2022a). Self-potential time series reveal emergent behavior in soil organic matter dynamics. Sci. Rep.12 (1), 13531. 10.1038/s41598-022-17914-5

  • 49

    DragunJ. (2021). “Recovery techniques and treatment technologies for petroleum and petroleum products in soil and groundwater,” in Petroleum contaminated soils (Florida, United States: CRC Press), 211217.

  • 50

    EarlN.CartwrightC. D.HorrocksS. J.WorboysM.SwiftS.KirtonJ. A.et al (2003). Review of the fate and transport of selected contaminants in the soil environment. Draft R and D Technical Report P5-079/TR1, 190.

  • 51

    EbuehiO.AbiboI.ShekwoloP.SigismundK.AdokiA.OkoroI. (2005). Remediation of crude oil contaminated soil by enhanced natural attenuation technique. J. Appl. Sci. Enviromental Manag.

  • 52

    EhirimC. N.AdizuaO. F.OkorieI. P. C. (2016). Geoelectrical characterization of matured PetroleumHydrocarbon impacted soil in port harcourt, Nigeria. Asian J. Earth Sci.9 (1), 915. 10.3923/ajes.2016.9.15

  • 53

    EjibaI. V.OnyaS. C.AdamsO. K. (2016). Impact of oil pollution on livelihood: Evidence from the Niger delta region of Nigeria. J. Sci. Res. Rep.12 (5), 112. 10.9734/jsrr/2016/26633

  • 54

    EspositoG.AllersmaH. G. B.SelvaduraiA. P. S. (1999). Centrifuge modeling of LNAPL transport in partially saturated sand. J. Geotechnical Geoenvironmental Eng.125 (12), 10661071. 10.1061/(asce)1090-0241(1999)125:12(1066

  • 55

    EssaidH. I.BekinsB. A.CozzarelliI. M. (2015). Organic contam-inant transport and fate in the subsurface: Evolution of knowledge and understanding. Water Resour. Researcg51 (7), 48614902. 10.1002/2015WR017121

  • 56

    EssaidH. I.HerkelrathW. N.HessK. M. (1993). Simulation of fluid distributions observed at a crude oil spill site incorporating hysteresis, oil entrapment, and spatial variability of hydraulic properties. Water Resour. Res.29 (6), 17531770. 10.1029/93wr00370

  • 57

    EtkinD. S.WelchJ. (2005). Oil spill intelligence report international oil spill database: Trends in oil spill volumes and frequency, Int. Oil Spill Conf., 2005. 10.7901/2169-3358-1997-1-949

  • 58

    FagerlundF.IllangasekareT. H.NiemiA. (2007). Nonaqueous-phase liquid infiltration and immobilization in heterogeneous media: 1. Experimental methods and two-layered reference case. Vadose Zone J.6 (3), 471482. 10.2136/vzj2006.0171

  • 59

    FagerlundF.NiemiA. (2007). A partially coupled, fraction-by-fraction modelling approach to the subsurface migration of gasoline spills. J. Contam. hydrology89 (3-4), 174198. 10.1016/j.jconhyd.2006.08.004

  • 60

    FangY. Z.MaZ. M.LuoY. Y. (2013). Research on vulnerability of aquifer to petroleum hydrocarbon basing on the ANN and GIS. Appl. Mech. Mater.319, 331336. 10.4028/www.scientific.net/amm.319.331

  • 61

    GaonkarO. D.Suresh KumarG.NambiI. M. (2016). Numerical modelling on fate and transport of coupled adsorption and biodegradation of pesticides in an unsaturated porous medium. ISH J. Hydraulic Eng.22 (3), 236246. 10.1080/09715010.2016.1166073

  • 62

    GengX.BoufadelM. C.CuiF. (2017). Numerical modeling of subsurface release and fate of benzene and toluene in coastal aquifers subjected to tides. J. Hydrology551, 793803. 10.1016/j.jhydrol.2016.10.039

  • 63

    GiadomF. D.AkpokodjeE. G.TseA. C. (2015). Determination of migration rates of contaminants in a hydrocarbon-polluted site using non-reactive tracer test in the Niger Delta, Nigeria. Environ. Earth Sci.74 (1), 879888. 10.1007/s12665-015-4094-3

  • 64

    GiadomF. D.TseA. C. (2014). Groundwater contamination and environmental risk assessment of a hydrocarbon contaminated site in Eastern Niger Delta, Nigeria. J. Environ. Earth Sci.5 (14), 123133.

  • 65

    GrathwohlP.EberhardtC.KlenkI.RügnerH.MaierU. (2020). “Transverse vertical dispersivity in aquifer materials: Implications on mass transfer across the capillary fringe and the length of steady state plumes,” in Groundwater (Florida, United States: CRC Press), 1920.

  • 66

    GuptaP. K.YadavB.YadavB. K. (2019). Assessment of LNAPL in subsurface under fluctuating groundwater table using 2D sand tank experiments. J. Environ. Eng.145 (9). 10.1061/(asce)ee.1943-7870.0001560

  • 67

    GüvenO.MolzF. J.MelvilleJ. G. (1986). Comment on “An Advection‐Diffusion Concept for Solute Transport in Heterogeneous Unconsolidated Geological Deposits” by Gillham et al. Water Resour. Res.22 (1), 8991. 10.1029/WR022i001p00089

  • 68

    HadleyP. W.NewellC. (2014). The new potential for understanding groundwater contaminant transport. Groundwater52 (2), 174186. 10.1111/gwat.12135

  • 69

    HadleyP. W.NewellC. (2014). The new potential for understanding groundwater contaminant transport. Groundwater52 (2), 656658. 10.1111/gwat.12246

  • 70

    HaggertyR.HarveyC. F.Freiherr von SchwerinC.MeigsL. C. (2004). What controls the apparent timescale of solute mass transfer in aquifers and soils? A comparison of experimental results. Water Resour. Res.40 (1). 10.1029/2002wr001716

  • 71

    HeibV. I.NeuweilerI.OchsS.FärberA. (2011). Experimental investigation on front morphology for two-phase flow in heterogeneous porous media. Water Resour. Res.47 (10). 10.1029/2011WR010612

  • 72

    HeldR. J.IllangasekareT. H. (1995a). Fingering of dense nonaqueous phase liquids in porous media: 1. Experimental investigation. Water Resour. Res.31 (5), 12131222. 10.1029/95WR00428

  • 73

    HochstetlerD. L.RolleM.ChiognaG.HabererC. M.GrathwohlP.KitanidisP. K. (2013). Effects of compound-specific transverse mixing on steady-state reactive plumes: Insights from pore-scale simulations and Darcy-scale experiments. Adv. water Resour.54, 110. 10.1016/j.advwatres.2012.12.007

  • 74

    HornS. A.NealC. P., (2005). The Atlantic empress sinking-a large spill without environmental disaster. 2005 International Oil Spill Conference, IOSC, 510

  • 75

    HuangX.ZhouW.DengD. (2020). Validation of pore network modeling for determination of two-phase transport in fibrous porous media. Sci. Rep.10 (1), 20852. 10.1038/s41598-020-74581-0

  • 76

    IllangasekareT. H.RamseyJ. L.JensenK. H.ButtsM. B. (1995). Experimental study of movement and distribution of dense organic contaminants in heterogeneous aquifers. J. Contam. Hydrology20 (1–2), 125. 10.1016/0169-7722(95)00045-W

  • 77

    IovenittiJ. L.NitaoJ. J.BishopD. J. (1992). Vadose zone investigations at the lawrence livermore national laboratory superfund site: An overview. Washington, D.C: Second USA/CIS Joint Conference Environmental Hydrology and Hydrogeology, 117.

  • 78

    ITOPF (2021). Oil tanker spill statistics. Available at: https://www.itopf.org/knowledge-resources/data-statistics/statistics/.

  • 79

    JarvisN.KoestelJ.LarsboM. (2016). Understanding preferential flow in the vadose zone: Recent advances and future prospects. Vadose Zone J.15 (12), 11. 10.2136/vzj2016.09.0075

  • 80

    KechavarziC.SogaK.IllangasekareT. H. (2005). Two-dimensional laboratory simulation of LNAPL infiltration and redistribution in the vadose zone. J. Contam. Hydrology76 (3–4), 211233. 10.1016/j.jconhyd.2004.09.001

  • 81

    KimH. S.SwansonJ. C. (2001). “Fate and transport modeling of contaminants,” in Salem sound. Report to the marine monitoring and research technical series. MMRTS-01-01.

  • 82

    KimM.HongS. H.WonJ.YimU. H.JungJ. H.HaS. Y.et al (2013). Petroleum hydrocarbon contaminations in the intertidal seawater after the Hebei Spirit oil spill – effect of tidal cycle on the TPH concentrations and the chromatographic characterization of seawater extracts. Water Res.47 (2), 758768. 10.1016/j.watres.2012.10.050

  • 83

    KoromaS.AratoA.GodioA. (2015). Analyzing geophysical signature of a hydrocarbon-contaminated soil using geoelectrical surveys. Environ. Earth Sci.74 (4), 29372948. 10.1007/s12665-015-4326-6

  • 84

    KovalickW. W.JrMontgomeryR. H., 2014. Developing a Program for contaminated site management in low and middle income countries.

  • 85

    LiZ.Galindo-TorresS.YanG.ScheuermannA.LiL. (2019). Pore-scale simulations of simultaneous steady-state two-phase flow dynamics using a lattice Boltzmann model: Interfacial area, capillary pressure and relative permeability. Transp. Porous Media129 (1), 295320. 10.1007/s11242-019-01288-w

  • 86

    LindénO.PålssonJ. (2013). Oil contamination in ogoniland, Niger delta. Ambio42 (6), 685701. 10.1007/s13280-013-0412-8

  • 87

    LittleD. I.HoltzmannK.GundlachE. R.GalperinY. (2018). Sediment hydrocarbons in former mangrove areas, southern ogoniland, eastern Niger delta, Nigeria. Coast. Res. Libr.25, 323342. 10.1007/978-3-319-73016-5_14

  • 88

    LiuM.ZhangY.HuL. (2010). Modeling the fate and transport of halohydrocarbons contaminants in groundwater aquifers of a northern city. 4th Int. Conf. Bioinforma. Biomed. Eng. ICBBE. 10.1109/ICBBE.2010.5515046

  • 89

    LuG.ZhengC.WolfsbergA., (2005). Effect of uncertain hydraulic conductivity on the simulated fate and transport of BTEX compounds at a field site. J. Environ. Eng., 131, 767, 776. 10.1061/(asce)0733-9372(2005)131:5(767)

  • 90

    MaierR. M. (2019). Biological processes affecting contaminants transport and fate. Environ. Pollut. Sci., 131146. 10.1016/b978-0-12-814719-1.00009-4

  • 91

    MajorW.VersteegR.Day-LewisF. D.JohnsonT. C.LaneJ. W. (2014). Optimized enhanced bioremediation through 4D geophysical monitoring and autonomous data collection. Process. Analysis ESTCP Proj.1–98.

  • 92

    MandocdocM.DavidC. P. (2006). Dieldrin contamination of the Groundwater in a former US military Base. Clean36, 1011.

  • 93

    MansiA. H.CastilloM. P.BernasconiG. (2017). Controlled laboratory test for the investigation of LNAPL contamination using a 2.0 GHz ground penetrating radar. Boll. Geofis. Teor. Appl.58 (3).

  • 94

    MaoD.LuL.RevilA.ZuoY.HintonJ.RenZ. J. (2016). Geophysical monitoring of hydrocarbon-contaminated soils remediated with a bioelectrochemical system. Environ. Sci. Technol.50 (15), 82058213. 10.1021/acs.est.6b00535

  • 95

    MichelJ.FingasM. (2016). “Oil Spills: Causes, consequences, prevention, and countermeasures,” in Fossil fuels: Current status and future directions (USA: University of South Carolina), 159201.

  • 96

    MoldanA. G. S.JacksonL. F.McGibbonS.Van Der WesthuizenJ. (1985). Some aspects of the Castillo de Bellver oilspill. Mar. Pollut. Bull.16 (3), 97102. 10.1016/0025-326X(85)90530-2

  • 97

    MolinsS.MayerK. U.AmosR. T.BekinsB. A. (2010). Vadose zone attenuation of organic compounds at a crude oil spill site - interactions between biogeochemical reactions and multicomponent gas transport. J. Contam. Hydrology112 (1–4), 1529. 10.1016/j.jconhyd.2009.09.002

  • 98

    MoreiraA.InnocentiC.HeleneL. P.CôrtesR. P. A. (2017). Integration of geoelectric and geochemical data in the evaluation of natural attenuation in a diesel contaminated site in São Manuel (Brazil). Geofísica Int.56 (3), 115. 10.22201/igeof.00167169p.2017.56.3.1814

  • 99

    Moreno-BarberoE.IllangasekareT. H. (2006). Influence of dense nonaqueous phase liquid pool morphology on the performance of partitioning tracer tests: Evaluation of the equilibrium assumption. Water Resour. Res.42 (4). 10.1029/2005WR004074

  • 100

    MoslenM.AigberuaA. (2018). Heavy metals and hydrocarbon contamination of surface water in azuabie Creek within bonny estuary, Nigeria. J. Appl. Sci. Environ. Manag.22 (7), 1083. 10.4314/jasem.v22i7.15

  • 101

    MugunthanP.MohanR.ThornburgT.MartinA.MahoneyM. (2016). Considerations for modeling contaminant transport at confined disposal sites. March1–31.

  • 102

    MyersT. (2012). Potential contaminant pathways from hydraulically fractured shale to aquifers. Groundwater50 (6), 872882. 10.1111/j.1745-6584.2012.00933.x

  • 103

    PrincewillN.ShvedaK., 2022, Shell escaped liability for oil spills in Nigeria for years.

  • 104

    NaiduR. (2013). Recent advances in contaminated site remediation. Water, Air, & Soil Pollut.224 (12), 17051711. 10.1007/s11270-013-1705-z

  • 105

    NaudetV.RevilA. (2005). A sandbox experiment to investigate bacteria‐mediated redox processes on self‐potential signals. Geophys. Res. Lett.32 (11), L11405. 10.1029/2005gl022735

  • 106

    NewellC. J.ConnorJ. A. (1998). Characteristics of dissolved petroleum Crude oil plumes. D.C., United States: American Petroleum Institute.

  • 107

    NoelC.GourryJ. C.DeparisJ.BlessingM.IgnatiadisI.GuimbaudC. (2016). Combining geoelectrical measurements and CO2Analyses to monitor the enhanced bioremediation of hydrocarbon-contaminated soils: A field implementation. Appl. Environ. Soil Sci.2016, 115. 10.1155/2016/1480976

  • 108

    NwachukwuA. N.OsuagwuJ. C. (2014). Effects of oil spillage on groundwater quality in Nigeria. Am. J. Eng. Res.3 (6), 271274.

  • 109

    NwankwoalaH.AmadiA.OmofuophuE.IbrahimH. (2020). Risk evaluation and modeling of soils contaminated with polycyclic aromatic hydrocarbons (PAHs) in parts of bonny island, Niger delta, Nigeria. Ann. Civ. Environ. Eng.4 (1), 015026. 10.29328/journal.acee.1001021

  • 110

    NwiloP. C.BadejoO. T. (2005). “Oil spill problems and management in the Niger Delta,” in International oil spill conference (D.C., United States: American Petroleum Institute), 567570.

  • 111

    OhanmuE. O.BakoS. P.OhanmuE.OhanmuO. O. (2018). Environmental implications, properties and attributes of crude oil in the oil-producing states of Nigeria. Ecologia9 (1), 19. 10.3923/ecologia.2019.1.9

  • 112

    OjuriO. O.OlaS. A.RudolphD. L.BarkerJ. F. (2010). Contamination potential of tar sand exploitation in the Western Niger-delta of Nigeria: Baseline studies. AU J. Technol.14 (2), 119128. 10.1007/s10064-009-0239-5

  • 113

    OlaS. A.FadugbaO. G.AdekoyaJ. A.BabatolaJ. O.AjayiO. O.OjuriO. O.et al (2019). Petroleum hydrocarbon contaminated groundwater remediation using 21st century technology: Baruwa community, Lagos State Nigeria as a case study. IOP Conf. Ser. Mater. Sci. Eng.640 (1), 012090. 10.1088/1757-899x/640/1/012090

  • 114

    OsaghaeE. E. (1995). The Ogoni uprising: Oil politics, minority agitation and the future of the Nigerian state. Afr. Aff.94 (376), 325344. 10.1093/oxfordjournals.afraf.a098833

  • 115

    OssaiI. C.AhmedA.HassanA.HamidF. S. (2020). Remediation of soil and water contaminated with petroleum hydrocarbon: A review. Environ. Technol. Innovation17, 100526. 10.1016/j.eti.2019.100526

  • 116

    OsujiL. C.EzebuiroP. E. (2006). Hydrocarbon contamination of a typical mangrove floor in Niger Delta, Nigeria. Int. J. Environ. Sci. Technol.3 (3), 313320. 10.1007/BF03325939

  • 117

    OvdatH.BerkowitzB. (2006). Pore-scale study of drainage displacement under combined capillary and gravity effects in index-matched porous media. Water Resour. Res.42 (6). 10.1029/2005WR004553

  • 118

    OyinkuroO.WariebiK. (2017). Hydrocarbon spill site characterization by electrical resistivity tomography and ground penetrating radar methods - a review. Asian J. Environ. Ecol.4 (3), 19. 10.9734/ajee/2017/36742

  • 119

    PanY.ZhangQ.YuY.TongY.WuW.ZhouY.et al (2021). Three-dimensional migration and resistivity characteristics of crude oil in heterogeneous soil layers. Environ. Pollut.268, 115309. 10.1016/j.envpol.2020.115309

  • 120

    ParkerJ. C. (1989). Multiphase flow and transport in porous media. Rev. Geophys.27 (3), 311328. 10.1029/rg027i003p00311

  • 121

    PepperI. L.GerbaC. P.GentryT. J. (2015). Environmental microbiology. San Diego, CA: Academic Press.

  • 122

    PetersG.SvanströmM. (2019). “Modelling environmental transport and fate of contaminants,” in Environmental sustainability for engineers and applied scientists (Cambridge, England: Cambridge University Press). 10.1017/9781316711408.005

  • 123

    PetriB. G.FučíkR.IllangasekareT. H.SmitsK. M.ChristJ. A.SakakiT.et al (2015). Effect of NAPL source morphology on mass transfer in the vadose zone. Groundwater53 (5), 685698. 10.1111/gwat.12284

  • 124

    QinX. S.HuangG. H.HeL. (2009). Simulation and optimization technologies for petroleum waste management and remediation process control. J. Environ. Manag.90 (1), 5476. 10.1016/j.jenvman.2008.07.002

  • 125

    RahbehM. E.MohtarR. H. (2007). Application of multiphase transport models to field remediation by air sparging and soil vapor extraction. J. Hazard. Mater.143 (1–2), 156170. 10.1016/j.jhazmat.2006.09.098

  • 126

    ReidenbachM. A.KoseffJ. R.MonismithS. G. (2007). Laboratory experiments of fine-scale mixing and mass transport within a coral canopy. Phys. Fluids19 (7), 075107. 10.1063/1.2752189

  • 127

    RenW.ErshadniaR.WallaceC. D.LaBolleE. M.DaiZ.de BarrosF. P.et al (2022). Evaluating the effects of multiscale heterogeneous sediments on solute mixing and effective dispersion. Water Resour. Res.58 (9), e2021WR031886. 10.1029/2021wr031886

  • 128

    RezaeiH.Bozorg-HaddadO.LoáicigaH. A. (2017). Effect of hydraulic conductivity uncertainty on in situ bioremediation of groundwater contaminated with dissolved petroleum Crude oils. J. Irrigation Drainage Eng.143 (12), 04017049.

  • 129

    RiccardiC.BerardiS.Di BasilioM.GariazzoC.GiardiP.VillariniM. (2001). Environmental assessment of a site contaminated by organic compounds. J. Environ. Sci. Health - Part A Toxic/Hazardous Subst. Environ. Eng.36 (6), 957970. 10.1081/ESE-100104124

  • 130

    RifaiH. S.BedientP. B. (1990). Comparison of biodegradation kinetics with an instantaneous reaction model for groundwater. Water Resour. Res.26 (4), 637645. 10.1029/WR026i004p00637

  • 131

    RolleM.KitanidisP. K. (2014). Effects of compound-specific dilution on transient transport and solute breakthrough: A pore-scale analysis. Adv. Water Resour.71, 186199. 10.1016/j.advwatres.2014.06.012

  • 132

    RolleM.MaierU.GrathwohlP. (2011). “Contaminant fate and reactive transport in groundwater,” in Dealing with contaminated sites) (Dordrecht: Springer), 851885.

  • 133

    SamK.CoulonF.PrpichG. (2017). Management of petroleum hydrocarbon contaminated sites in Nigeria: Current challenges and future direction. Land Use Policy64, 133144. 10.1016/j.landusepol.2017.01.051

  • 134

    ScheibeT. D.MurphyE. M.ChenX.RiceA. K.CarrollK. C.PalmerB. J.et al (2015a). An analysis platform for multiscale hydrogeologic modeling with emphasis on hybrid multiscale methods. Groundwater53, 3856. 10.1111/gwat.12179

  • 135

    SeigneurN.MayerK. U.SteefelC. I. (2019). Reactive transport in evolving porous media. Rev. Mineralogy Geochem.85 (1), 197238. 10.2138/rmg.2019.85.7

  • 136

    SethiR.Di MolfettaA. (2019). “Transport of immiscible fluids,” in Groundwater engineering (Cham: Springer), 249262.

  • 137

    ŠimůnekJ.van GenuchtenM. T. (2016). “Contaminant transport in the unsaturated zone: Theory and modeling,” in The handbook of groundwater engineering (Amsterdam The Nethherland: CRC Press), 221254.

  • 138

    SiracusaG.La RosaA. D.MusumeciL.MaiolinoG. (2007). Modelling of contaminant migration in unsaturated soils. WIT Trans. Ecol. Environ.102. 10.2495/SDP070441

  • 139

    SteefelC. I.DePaoloD. J.LichtnerP. C. (2005). Reactive transport modeling: An essential tool and a new research approach for the Earth sciences. Earth Planet. Sci. Lett.240 (3-4), 539558. 10.1016/j.epsl.2005.09.017

  • 140

    SwartjesF. A.GrimaJ. (2011). “Groundwater-related risk assessment,” in Dealing with contaminated sites (Dordrecht: Springer), 749785.

  • 141

    SwartjesF. A. (2011). “Introduction to contaminated site management,” in Dealing with contaminated sites (Dordrecht: Springer), 389. 10.1007/978-90-481-9757-6

  • 142

    SwartjesF. A.RutgersM.LijzenJ. P. A.JanssenP. J. C. M.OtteP. F.WintersenA.et al (2012). State of the art of contaminated site management in The Netherlands: Policy framework and risk assessment tools. Sci. Total Environ.427, 110. 10.1016/j.scitotenv.2012.02.078

  • 143

    TheodorakisC.WalkerT. D. (2014). “Environmental processes,” in Encyclopedia of toxicology. Third Edition3 (Amsterdam, Netherlands: Elsevier). 10.1016/B978-0-12-386454-3.01011-32

  • 144

    UdueborM. A.OlaS. A. (2016). Hydrocarbon remediation by natural attenuation at Baruwa, Lagos Nigeria. Electron. J. Geotechnical Eng.21, 501512.

  • 145

    UgbenaK. G.TseA. C.AkpokodjeE. G. (2020). Modeling soluble contaminant migration by advection process in subsurface water in the Eastern Niger Delta, Nigeria. World Sci. News139 (2), 102114.

  • 146

    UkoM.UdotongI.OfonU.UmanaS.AbrahamN. (2020). Effect of hydrocarbon contamination on the microbial diversity of freshwater sediments within akwa ibom state, Nigeria. J. Chem. Environ. Biol. Eng.8, 32. 10.11648/j.jcebe.20200402.11

  • 147

    ValleroD. (2021). Environmental systems science: Theory and practical applications. Amsterdam, Netherlands: Elsevier.

  • 148

    Van GeelP. J.SykesJ. F. (1994). Laboratory and model simulations of a LNAPL spill in a variably-saturated sand, 1. Laboratory experiment and image analysis techniques. J. Contam. Hydrology17 (1), 125. 10.1016/0169-7722(94)90075-2

  • 149

    WeaverJ. W. (1996). “November. Application of the Crude oil spill screening model to field sites,” in Non-aqueous phase liquids (NAPLs) in subsurface environment: Assessment and remediation, 788799.

  • 150

    WicksF. (2009). The oil age. Mech. Eng.131 (08), 4245. 10.1115/1.2009-aug-6

  • 151

    WilliamsS. D.LaddD. E.FarmerJ. J. (2006). Fate and transport of petroleum crude oils in soil and ground water at big south fork national river and recreation area. Washington: Tennessee and Kentucky, 20022003. U. S. Geological Survey Scientific Investigations Report 2005-5104.

  • 152

    WinterT. C. (1999). Ground water and surface water: A single resource. Collingdale, Pennsylvania: Diane Publishing, 1139.

  • 153

    XiaY.LiH.BoufadelM. C.SharifiY. (2010). Hydrodynamic factors affecting the persistence of the Exxon Valdez oil in a shallow bedrock beach. Water Resour. Res.46 (10). 10.1029/2010WR009179

  • 154

    YoonH.ChojnickiK. N.MartinezM. J. (2019). Pore-scale analysis of calcium carbonate precipitation and dissolution kinetics in a microfluidic device. Environ. Sci. Technol.53 (24), 1423314242. 10.1021/acs.est.9b01634

  • 155

    YouX.LiuS.DaiC.GuoY.ZhongG.DuanY. (2020). “Contaminant occurrence and migration between high- and low-permeability zones in groundwater systems: A review,” in Science of the total environment (Ammsterdam Netherland: Elsevier B.V), 743. 10.1016/j.scitotenv.2020.140703

  • 156

    ZabbeyN.SamK.OnyebuchiA. T. (2017). Remediation of contaminated lands in the Niger Delta, Nigeria: Prospects and challenges. Sci. Total Environ.586, 952965. 10.1016/j.scitotenv.2017.02.075

  • 157

    ZabbeyN.UyiH. (2014). Community responses of intertidal soft-bottom macrozoobenthos to oil pollution in a tropical mangrove ecosystem, Niger Delta, Nigeria. Mar. Pollut. Bull.82 (1), 167174. 10.1016/j.marpolbul.2014.03.002

  • 158

    ZhengC.WeaverJ.TonkinM. (2010). MT3DMS, A modular three-dimensional multispecies transport model – user guide to the crude oil spill source (HSS) package MT3DMS. D.C., United States: U.S. Environmental Protection Agency.

Summary

Keywords

crude oil, soil, groundwater, contaminant, fate and transport

Citation

Adeniran MA, Oladunjoye MA and Doro KO (2023) Soil and groundwater contamination by crude oil spillage: A review and implications for remediation projects in Nigeria. Front. Environ. Sci. 11:1137496. doi: 10.3389/fenvs.2023.1137496

Received

04 January 2023

Accepted

04 April 2023

Published

09 May 2023

Volume

11 - 2023

Edited by

Sören Thiele-Bruhn, University of Trier, Germany

Reviewed by

Bassam Tawabini, King Fahd University of Petroleum and Minerals, Saudi Arabia

Manuel Cánovas, Catholic University of the North, Chile

Updates

Copyright

*Correspondence: Margaret A. Adeniran,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics