Groundwater vulnerability and risk mapping using GIS, modeling and a fuzzy logic tool

A groundwater vulnerability and risk mapping assessment, based on a source-pathway-receptor approach, is presented for an urban coastal aquifer in northeastern Brazil. A modified version of the DRASTIC methodology was used to map the intrinsic and specific groundwater vulnerability of a 292 km(2) study area. A fuzzy hierarchy methodology was adopted to evaluate the potential contaminant source index, including diffuse and point sources. Numerical modeling was performed for delineation of well capture zones, using MODFLOW and MODPATH. The integration of these elements provided the mechanism to assess groundwater pollution risks and identify areas that must be prioritized in terms of groundwater monitoring and restriction on use. A groundwater quality index based on nitrate and chloride concentrations was calculated, which had a positive correlation with the specific vulnerability index.     

Ecohydrological Responses to Diversion of Groundwater: Case Study of a Deep-Rock Repository for Spent Nuclear Fuel in Sweden

Planning and license applications concerning groundwater diversion in areas containing water-dependent or water-favored habitats must take into account both hydrological effects and associated ecological consequences. There is at present no established methodology to assess such ecohydrological responses. Thus, this paper describes a new stepwise methodology to assess ecohydrological responses to groundwater diversion from, e.g., water-drained pits, shafts, tunnels, and caverns in rock below the groundwater table. The methodology is illustrated using the planned deep-rock repository for spent nuclear fuel at Forsmark in central Sweden as a case study, offering access to a unique hydrological and ecological dataset. The case study demonstrates that results of ecohydrological assessments can provide useful inputs to planning of monitoring programs and mitigation measures in infrastructure projects. As a result of the assessment, artificial water supply to wetlands is planned in order to preserve biological diversity, nature values, and vulnerable species.     

Predictions of long-term performance of granular iron permeable reactive barriers: field-scale evaluation

Long-term performance is a key consideration for the granular iron permeable reactive barrier (PRB) technology because the economic benefit relies on sustainable operation for substantial periods of time. However, predictions on the long-term performance have been limited mainly because of the lack of reliable modeling tools. This study evaluated the predictive capability of a recently-developed reactive transport model at two field-scale PRBs, both having relatively high concentrations of dissolved carbonate in the native groundwater. The first site, with 8 years of available monitoring data, was a funnel-and-gate installation, with a low groundwater velocity through the gate (about 0.12 m d(-1)). The loss in iron reactivity caused by secondary mineral precipitation was small, maintaining relatively high removal rates for chlorinated organics. The simulated concentrations for most constituents in the groundwater were within the range of the monitoring data. The second site, with monitoring data available for 5 years, was a continuous wall PRB, designed for a groundwater velocity of 0.9 m d(-1). A comparison of measured and simulated aqueous concentrations suggested that the average groundwater velocity through the PRB could be lower than the design value by a factor of two or more. The distribution and amounts of carbonate minerals measured in core samples supported the decreased groundwater velocity used in the simulation. The generally good agreement between the simulated and measured aqueous and solid-phase data suggest that the model could be an effective tool for predicting long-term performance of granular iron PRBs, particularly in groundwater with high concentrations of carbonate.     

Impacts of swine manure pits on groundwater quality

Manure deep-pits are commonly used to store manure at confined animal feeding operations. However, previous to this study little information had been collected on the impacts of deep-pits on groundwater quality to provide science-based guidance in formulating regulations and waste management strategies that address risks to human health and the environment. Groundwater quality has been monitored since January 1999 at two hog finishing facilities in Illinois that use deep-pit systems for manure storage. Groundwater samples were collected on a monthly basis and analyzed for inorganic and bacteriological constituent concentrations. The two sites are located in areas with geologic environments representing different vulnerabilities for local groundwater contamination. One site is underlain by more than 6 m of clayey silt, and 7-36 m of shale. Concentrations of chloride, ammonium, phosphate, and potassium indicated that local groundwater quality had not been significantly impacted by pit leakage from this facility. Nitrate concentrations were elevated near the pit, often exceeding the 10 mg N/l drinking water standard. Isotopic nitrate signatures suggested that the nitrate was likely derived from soil organic matter and fertilizer applied to adjacent crop fields. At the other site, sandstone is located 4.6-6.1 m below land surface. Chloride concentrations and delta15N and delta15O values of dissolved nitrate indicated that this facility may have limited and localized impacts on groundwater. Other constituents, including ammonia, potassium, phosphate, and sodium were generally at or less than background concentrations. Trace- and heavy-metal concentrations in groundwater samples collected from both facilities were at concentrations less than drinking water standards. The concentration of inorganic constituents in the groundwater would not likely impact human health. Fecal streptococcus bacteria were detected at least once in groundwater from all monitoring wells at both sites. Fecal streptococcus was more common and at greater concentrations than fecal coliform. The microbiological data suggest that filtration of bacteria by soils may not be as effective as commonly assumed. The presence of fecal bacteria in the shallow groundwater may pose a significant threat to human health if the ground water is used for drinking. Both facilities are less than 4 years old and the short-term impacts of these manure storage facilities on groundwater quality have been limited. Continued monitoring of these facilities will determine if they have a long-term impact on groundwater resources.     

Persistence of LNAPL sources: relationship between risk reduction and LNAPL recovery

Light nonaqueous phase liquids (LNAPLs), such as fuels, are the source of much soil and groundwater contamination. Though the mobility of LNAPLs is limited in many environments, dissolved-phase components, such as benzene, can produce groundwater plumes that are more mobile than the LNAPL source. In such a setting, it is commonly assumed that recovery of the LNAPL will result in a reduction in risk associated with the dissolved phase. This paper synthesizes several existing multiphase and chemical transport solutions into a single linked methodology that predicts concentrations of soluble constituents within and downgradient of LNAPL source zones from dissolution of those constituents into groundwater flowing through and below LNAPL sources. This approach has been applied to a variety of LNAPL spill conditions. For biodegradable compounds, these analyses show that the period of time where the dissolved-phase plume is expanding is very small compared to the duration of most LNAPL sources, and that the downgradient extent is generally less than about 100 m for BTEX compounds. Therefore, the risk to receptors, as measured by the maximum downgradient extent of dissolved-phase plume or the maximum concentration of these compounds at a downgradient receptor, is generally unrelated to the longevity of the LNAPL sources. The maximum downgradient extent of the dissolved-phase plume is determined almost entirely by the groundwater velocity and the biodegradation rate. These analyses further demonstrate that recovery of LNAPL by hydraulic methods is often ineffective at reducing risk. Except in coarse-grained soils or intermediate soils with significant LNAPL saturations, free-product recovery approaches do not result in significant reductions in the longevity of downgradient dissolved-phase contamination. Further, for biodegradable constituents, remediation does not result in a near-term decrease in the downgradient extent of contamination. Cleanup methods that act to change the composition of the LNAPL source are more effective at reducing the downgradient concentrations, particularly for fine-grained soils or when LNAPL saturations are low.     

Temporal variability of groundwater chemistry in shallow and deep aquifers of Araihazar, Bangladesh

Samples were collected every 2-4 weeks from a set of 37 monitoring wells over a period of 2-3 years in Araihazar, Bangladesh, to evaluate the temporal variability of groundwater composition for As and other constituents. The monitoring wells are grouped in 6 nests and span the 5-91 m depth range. Concentrations of As, Ca, Fe, K, Mg, Mn, Na, P, and S were measured by high-resolution ICPMS with a precision of 5% or better; concentrations of Cl were measured by ion chromatography. In shallow wells <30 m deep, As and P concentrations generally varied by <30%, whereas concentrations of the major ions (Na, K, Mg, Ca and Cl) and the redox-sensitive elements (Fe, Mn, and S) varied over time by up to +/-90%. In wells tapping the deeper aquifers >30 m often below clay layers concentrations of groundwater As were much lower and varied by <10%. The concentrations of major cations also varied by <10% in these deep aquifers. In contrast, the concentration of redox-sensitive constituents Fe, S, and Mn in deep aquifers varied by up to 97% over time. Thus, strong decoupling between variations in As and Fe concentrations is evident in groundwaters from shallow and deep aquifers. Comparison of the time series data with groundwater ages determined by (3)H/(3)He and (14)C dating shows that large seasonal or inter-annual variations in major cation and chloride concentrations are restricted to shallow aquifers and groundwater recharged <5 years ago. There is no corresponding change in As concentrations despite having significant variations of redox sensitive constituents in these very young waters. This is attributed to chemical buffering due to rapid equilibrium between solute and solid As. At two sites where the As content of groundwater in existing shallow wells averages 102 microg/L (range: <5 to 648 microg/L; n=118) and 272 microg/L (range: 10 to 485 microg/L; n=65), respectively, a systematic long-term decline in As concentrations lends support to the notion that flushing may slowly deplete an aquifer of As. Shallow aquifer water with >5 years (3)H/(3)He age show a constant As:P molar ratio of 9.6 over time, suggesting common mechanisms of mobilization.     

Integrated methodology for assessing the HCH groundwater pollution at the multi-source contaminated mega-site Bitterfeld/Wolfen

A large-scale groundwater contamination characterises the Pleistocene groundwater system of the former industrial and abandoned mining region Bitterfeld/Wolfen, Eastern Germany. For more than a century, local chemical production and extensive lignite mining caused a complex contaminant release from local production areas and related dump sites. Today, organic pollutants (mainly organochlorines) are present in all compartments of the environment at high concentration levels. An integrated methodology for characterising the current situation of pollution as well as the future fate development of hazardous substances is highly required to decide on further management and remediation strategies. Data analyses have been performed on regional groundwater monitoring data from about 10 years, containing approximately 3,500 samples, and up to 180 individual organic parameters from almost 250 observation wells. Run-off measurements as well as water samples were taken biweekly from local creeks during a period of 18 months. A kriging interpolation procedure was applied on groundwater analytics to generate continuous distribution patterns of the nodal contaminant samples. High-resolution geological 3-D modelling serves as a database for a regional 3-D groundwater flow model. Simulation results support the future fate assessment of contaminants. A first conceptual model of the contamination has been developed to characterise the contamination in regional surface waters and groundwater. A reliable explanation of the variant hexachlorocyclohexane (HCH) occurrence within the two local aquifer systems has been derived from the regionalised distribution patterns. Simulation results from groundwater flow modelling provide a better understanding of the future pollutant migration paths and support the overall site characterisation. The presented case study indicates that an integrated assessment of large-scale groundwater contaminations often needs more data than only from local groundwater monitoring. The developed methodology is appropriate to assess POP-contaminated mega-sites including, e.g. HCH deposits. Although HCH isomers are relevant groundwater pollutants at this site, further organochlorine pollutants are present at considerably higher levels. The study demonstrates that an effective evaluation of the current situation of contamination as well as of the related future fate development requires detailed information of the entire observed system.     

Analysis of munitions constituents in groundwater using a field-portable GC-MS

The use of munitions constituents (MCs) at military installations can produce soil and groundwater contamination that requires periodic monitoring even after training or manufacturing activities have ceased. Traditional groundwater monitoring methods require large volumes of aqueous samples (e.g., 2-4 L) to be shipped under chain of custody, to fixed laboratories for analysis. The samples must also be packed on ice and shielded from light to minimize degradation that may occur during transport and storage. The laboratory’s turn-around time for sample analysis and reporting can be as long as 45 d. This process hinders the reporting of data to customers in a timely manner; yields data that are not necessarily representative of current site conditions owing to the lag time between sample collection and reporting; and incurs significant shipping costs for samples.The current work compares a field portable Gas Chromatograph-Mass Spectrometer (GC-MS) for analysis of MCs on-site with traditional laboratory-based analysis using High Performance Liquid Chromatography with UV absorption detection. The field method provides near real-time (within ∼1 h of sampling) concentrations of MCs in groundwater samples. Mass spectrometry provides reliable confirmation of MCs and a means to identify unknown compounds that are potential false positives for methods with UV and other non-selective detectors.     

Application of Unknown Groundwater Pollution Source Release History Estimation Methodology to Distributed Sources Incorporating Surface-Groundwater Interactions

Sources of contamination of groundwater are often difficult to characterize. However, it is essential for effective remediation of polluted groundwater resources. This study demonstrates an application of the linked simulation-optimization based methodology to estimate the release history from spatially distributed sources of pollution at an illustrative abandoned mine-site. In linked simulation-optimization approaches a numerical groundwater flow and transport simulation model is linked to the optimization model. In this study, topographic and geologic characteristics of the abandoned mine-site were simulated using a three-dimensional (3D) numerical groundwater flow model. Transport of contaminant in the groundwater was simulated using a 3D transient advective-dispersive contaminant transport model. Adsorption or chemical reaction of the contaminant was not considered in the contaminant transport model. Adaptive simulated annealing (ASA) was employed for solving the optimization problem. An optimization algorithm generates the candidate solutions corresponding to various unknown groundwater source characteristics. The candidate solutions are used as input in the numerical groundwater transport simulation model to generate the concentration of pollutant in the study area. This information is used to calculate the objective function value, which is utilized by the optimization algorithm to improve the candidate solution. This process continues until an optimal solution is obtained. Optimal solutions obtained in this study show that the linked simulation-optimization based methodology is potentially applicable for the characterization of spatially distributed pollutant sources, typically present at abandoned mine-sites.     

Arsenic forensic modeling leads to decommission of downgradient/off-site monitoring wells

Arsenic was detected at concentrations exceeding the regulatory limit of 0.010 milligrams per liter (mg/L) in an off-site bedrock monitoring well downgradient of a former electroplating facility in Merrimack, New Hampshire. The bedrock underlying the site is associated with naturally occurring high concentrations of groundwater arsenic. Geochemical modeling was used to evaluate whether the arsenic in bedrock groundwater at the off-site monitoring location was site-related or naturally occurring. The hydrogeochemical signature of the off-site bedrock well did not resemble signatures of site-impacted bedrock wells. Multiple lines of evidence support that the arsenic observed in off-site bedrock groundwater was not a result of adverse impacts from site-related groundwater contamination.     

Cost-effective sampling network design for contaminant plume monitoring under general hydrogeological conditions

A new simulation-optimization methodology is developed for cost-effective sampling network design associated with long-term monitoring of large-scale contaminant plumes. The new methodology is similar in concept to the one presented by Reed et al. (Reed, P.M., Minsker, B.S., Valocchi, A.J., 2000a. Cost-effective long-term groundwater monitoring design using a genetic algorithm and global mass interpolation. Water Resour. Res. 36 (12), 3731-3741) in that an optimization model based on a genetic algorithm is coupled with a flow and transport simulator and a global mass estimator to search for optimal sampling strategies. However, this study introduces the first and second moments of a three-dimensional contaminant plume as new constraints in the optimization formulation, and demonstrates the proposed methodology through a real-world application. The new moment constraints significantly increase the accuracy of the plume interpolated from the sampled data relative to the plume simulated by the transport model. The plume interpolation approaches employed in this study are ordinary kriging (OK) and inverse distance weighting (IDW). The proposed methodology is applied to the monitoring of plume evolution during a pump-and-treat operation at a large field site. It is shown that potential cost savings up to 65.6% may be achieved without any significant loss of accuracy in mass and moment estimations. The IDW-based interpolation method is computationally more efficient than the OK-based method and results in more potential cost savings. However, the OK-based method leads to more accurate mass and moment estimations. A comparison of the sampling designs obtained with and without the moment constraints points to their importance in ensuring a robust long-term monitoring design that is both cost-effective and accurate in mass and moment estimations. Additional analysis demonstrates the sensitivity of the optimal sampling design to the various coefficients included in the objective function of the optimization model.     

Development of a risk assessment methodology for evaluating potential impacts associated with contaminated mud disposal in the marine environment

In order to assess impacts associated with disposal of contaminated mud arising from Hong Kong's dredging and reclamation projects, a methodology has been formulated to determine the level of risk posed by consumption of seafood/marine prey species to humans and to the Chinese White Dolphin (Sousa chinensis). This methodology improves on previously used techniques by incorporating risks for organic contaminants, accounting for doses from sources other than seafood, and incorporating additional local knowledge on Sousa chinensis behaviour. It thus represents an advance in risk assessment techniques and a new integration of risk assessment and monitoring in environmental management.     

Analysis of groundwater contamination using concentration-time series recorded during an integral pumping test: bias introduced by strong concentration gradients within the plume

When only few monitoring wells are available to assess the extent and level of groundwater contamination, inversion of concentration breakthrough curves acquired during an integral pumping test can be used as an alternative quantification method. The idea is to use concentration-time series recorded during integral pumping tests through an inversion technique to estimate contaminant mass fluxes crossing a control plane. In this paper, we examine how a longitudinal concentration gradient along a contaminant plume length scale affects the estimated inversed-concentration distribution and its associated mass flux. The analytically inversed-concentration distribution at the imaginary control plane (ICP) is compared to a numerically generated concentration distribution, treating the latter one as a "real contaminant plume" characterized by the presence of a longitudinal concentration gradient. It is found that the analytically inversed-concentration can lead to overestimation or underestimation of concentration distribution values depending on the transport time period and dispersivity values. At lower dispersivity values, with shorter transport time periods, the analytically inversed-concentration distribution overestimates the "real" concentration distribution. A better fit of the estimated concentration distribution to the "real" one is observed when the transport time period increases, i.e. when the advective front has already crossed the ICP. However, for higher dispersivity values, underestimation of the real concentration distribution is observed. Deviation of the inversed-concentration distribution from the "real" one is assessed for a site-specific concentration gradient term. A concentration gradient adjusted contaminant mass flux is thus formulated to evaluate groundwater contamination levels at a given time period through an ICP. This concentration gradient ratio can indicate whether the ICP is well positioned to evaluate accurately contaminant mass fluxes which are representative of groundwater contamination levels.     

Analysis of groundwater contamination using concentration–time series recorded during an integral pumping test: Bias introduced by strong concentration gradients within the plume

When only few monitoring wells are available to assess the extent and level of groundwater contamination, inversion of concentration breakthrough curves acquired during an integral pumping test can be used as an alternative quantification method. The idea is to use concentration–time series recorded during integral pumping tests through an inversion technique to estimate contaminant mass fluxes crossing a control plane. In this paper, we examine how a longitudinal concentration gradient along a contaminant plume length scale affects the estimated inversed-concentration distribution and its associated mass flux. The analytically inversed-concentration distribution at the imaginary control plane (ICP) is compared to a numerically generated concentration distribution, treating the latter one as a “real contaminant plume” characterized by the presence of a longitudinal concentration gradient. It is found that the analytically inversed-concentration can lead to overestimation or underestimation of concentration distribution values depending on the transport time period and dispersivity values. At lower dispersivity values, with shorter transport time periods, the analytically inversed-concentration distribution overestimates the “real” concentration distribution.A better fit of the estimated concentration distribution to the “real” one is observed when the transport time period increases, i.e. when the advective front has already crossed the ICP. However, for higher dispersivity values, underestimation of the real concentration distribution is observed. Deviation of the inversed-concentration distribution from the “real” one is assessed for a site-specific concentration gradient term. A concentration gradient adjusted contaminant mass flux is thus formulated to evaluate groundwater contamination levels at a given time period through an ICP. This concentration gradient ratio can indicate whether the ICP is well positioned to evaluate accurately contaminant mass fluxes which are representative of groundwater contamination levels.     

Heuristic space–time design of monitoring wells for contaminant plume characterization in stochastic flow fields

An optimization methodology for designing groundwater quality monitoring networks applicable to stochastic flow fields is presented and evaluated. The approach sets itself apart from previous techniques by incorporating the time dimension directly into the objective function. This function is extremized using a directed partial enumeration strategy guided by physical considerations related to transport processes. The result is a set of monitoring well locations and a sampling schedule that minimizes plume characterization error while satisfying constraints on the maximum number of wells and allowable number of active wells. The method is evaluated using hypothetical plumes with varying degrees of heterogeneity. Results indicate that the proposed approach is successful in generating near-optimal sampling networks that satisfy all imposed constraints. Monitoring networks with as little as three active wells and a total of 12 wells are found to provide adequate plume characterization for low toxicity contaminants.     

Treatment of groundwater contaminated with gasoline components by an ozone/UV process

In this paper, the treatment of real groundwater samples contaminated with gasoline components, such as benzene, toluene, ethylbenzene, and xylene (BTEX), methyl tert-butyl ether (MTBE), tert-butyl alcohol (TBA), and other gasoline constituents in terms of total petroleum hydrocarbons as gasoline (TPHg) by an ozone/UV process was investigated. The treatment was conducted in a semi-batch reactor under different experimental conditions by varying ozone gas dosage and incident UV light intensity. The groundwater samples contained BTEX compounds, MTBE, TBA, and TPHg in the ranges of 5-10000, 3000-5500, 80-1400, and 2400-20000mugl(-1), respectively. The ozone/UV process was very effective compared to ozonation in the removal of the gasoline components from the groundwater samples. For the various gasoline constituents, more than 99% removal efficiency was achieved for the ozone/UV process and the removal efficiency for ozonation was as low as 27%. The net ozone consumed per mol of organic carbon (from BTEX, MTBE, and TBA) oxidized varied in the range of 5-60 for different types of groundwater samples treated by the ozone/UV process. In ozonation experiments, it was observed that the presence of sufficient amount of iron in groundwater samples improved the removal of BTEX, MTBE, TBA, and TPHg.     

Molarity provides better clarity – using molar-based data when evaluating chlorinated volatile organic compound impacted sites

Abstract

Chlorinated volatile organic compounds are common constituents observed at many contaminated groundwater sites. Common industry practice has been to measure these constituents in groundwater on a weight concentration basis (e.g. µg/l, mg/l). This paper highlights the use of molar-based concentrations, especially in the case of parent-daughter degradation sequences commonly observed with the chloroethene, chloroethane, and chloromethane families of compounds. Converting to molar-based concentrations provides the practitioner greater insight into groundwater plume behavior including better evaluation of degradation processes, remedial progress, possible commingling, and/or sourcing. For example, this paper provides a tank analogy to evaluate whether the commonly misinterpreted observation of “DCE-stall” may be occurring at a site. Multiple examples of the benefits of using molar-based concentration data are also summarized in a project case study presented herein. As demonstrated in this paper, molarity does provide better clarity and can be a powerful evaluation tool in the groundwater practitioner’s toolbox.     

The fate of hydrogen peroxide as an oxygen source for bioremediation activities within saturated aquifer systems

In situ bioremediation is an innovative technique for the remediation of contaminated aquifers that involves the use of microorganisms to remediate soils and groundwaters polluted by hazardous substances. During its application, this process may require the addition of nutrients and/or electron acceptors to stimulate appropriate biological activity. Hydrogen peroxide has been commonly used as an oxygen source because of the limited concentrations of oxygen that can be transferred into the groundwater using above-ground aeration followed by reinjection of the oxygenated groundwater into the aquifer or subsurface air sparging of the aquifer. Because of several potential interactions of H2O2 with various aquifer material constituents, its decomposition may be too rapid, making effective introduction of the H2O2 into targeted treatment zones extremely difficult and costly. Therefore, a bench-scale study was conducted to determine the fate of H2O2 within subsurface aquifer environments. The purpose of this investigation was to identify those aquifer constituents, both biotic and abiotic, that are most active in controlling the fate of H2O2. The decomposition rates of H2O2 were determined using both equilibrated water samples and soil slurries. Results showed H2O2 decomposition to be effected by several commonly found inorganic soil components; however, biologically mediated catalytic reactions were determined to be the most substantial.     

A Monitoring Network Design Procedure for Three-Dimensional (3D) Groundwater Contaminant Source Identification

Finding the location and concentration of contaminant sources is an important step in groundwater remediation and management. This discovery typically requires the solution of an inverse problem. This inverse problem can be formulated as an optimization problem where the objective function is the sum of the square of the errors between the observed and predicted values of contaminant concentration at the observation wells. Studies show that the source identification accuracy is dependent on the observation locations (i.e., network geometry) and frequency of sampling; thus, finding a set of optimal monitoring well locations is very important for characterizing the source. The objective of this study is to propose a sensitivity-based method for optimal placement of monitoring wells by incorporating two uncertainties: the source location and hydraulic conductivity. An optimality metric called D-optimality in combination with a distance metric, which tends to make monitoring locations as far apart from each other as possible, is developed for finding optimal monitoring well locations for source identification. To address uncertainty in hydraulic conductivity, an integration method of multiple well designs is proposed based on multiple hydraulic conductivity realizations. Genetic algorithm is used as a search technique for this discrete combinatorial optimization problem. This procedure was applied to a hypothetical problem based on the well-known Borden Site data in Canada. The results show that the criterion-based selection proposed in this paper provides improved source identification performance when compared to uniformly distributed placement of wells.     

A partially coupled, fraction-by-fraction modelling approach to the subsurface migration of gasoline spills

The subsurface spreading behaviour of gasoline, as well as several other common soil- and groundwater pollutants (e.g. diesel, creosote), is complicated by the fact that it is a mixture of hundreds of different constituents, behaving differently with respect to e.g. dissolution, volatilisation, adsorption and biodegradation. Especially for scenarios where the non-aqueous phase liquid (NAPL) phase is highly mobile, such as for sudden spills in connection with accidents, it is necessary to simultaneously analyse the migration of the NAPL and its individual components in order to assess risks and environmental impacts. Although a few fully coupled, multi-phase, multi-constituent models exist, such models are highly complex and may be time consuming to use. A new, somewhat simplified methodology for modelling the subsurface migration of gasoline while taking its multi-constituent nature into account is therefore introduced here. Constituents with similar properties are grouped together into eight fractions. The migration of each fraction in the aqueous and gaseous phases as well as adsorption is modelled separately using a single-constituent multi-phase flow model, while the movement of the free-phase gasoline is essentially the same for all fractions. The modelling is done stepwise to allow updating of the free-phase gasoline composition at certain time intervals. The output is the concentration of the eight different fractions in the aqueous, gaseous, free gasoline and solid phases with time. The approach is evaluated by comparing it to a fully coupled multi-phase, multi-constituent numerical simulator in the modelling of a typical accident-type spill scenario, based on a tanker accident in northern Sweden. Here the PCFF method produces results similar to those of the more sophisticated, fully coupled model. The benefit of the method is that it is easy to use and can be applied to any single-constituent multi-phase numerical simulator, which in turn may have different strengths in incorporating various processes. The results demonstrate that the different fractions have significantly different migration behaviours and although the methodology involves some simplifications, it is a considerable improvement compared to modelling the gasoline constituents completely individually or as one single mixture.