Modeling of vapor intrusion from hydrocarbon-contaminated sources accounting for aerobic and anaerobic biodegradation

Soil benzene emissions and inhalation health risks affected by various land covers

Accurate quantification of soil volatile organic compounds (VOCs) flux is crucial for assessing inhalation environmental health risks and developing region-specific remediation strategies. However, land cover significantly influences VOCs emissions from soil. This study investigated benzene, a representative VOCs, using a laboratory flux chamber and numerical simulations to evaluate its release patterns under different surface covers, including bare soil (no cover), clay brick, cement, and grass. In the experiment, gaseous benzene was collected using an adsorption tube filled with Tenax-TA adsorbent. The collected samples were subsequently analyzed using thermal desorption coupled with gas chromatography-mass spectrometry. By integrating these findings with environmental health risk assessment methodologies, we developed a tailored approach for assessing inhalation health risks at benzene-contaminated sites with varying land covers. Additionally, we conducted application studies of this method across various scenarios. The results indicate that soil benzene emissions could be reduced by using low-permeability coverings such as clay brick and cement, as well as by planting vegetation. The average fluxes of benzene through covering materials were of the order of 1.22 × 10-2, 4.37 × 10-3, 2.47 × 10-3, and 9.88 × 10-4 mg m-2·s-1 for bare soil, clay brick, grass, and cement, respectively. The application of clay brick and cement coverings on the soil surface results in more pollutants remaining in the soil in liquid and adsorbed states, making them less likely to volatilize. The inhalation carcinogenic risk (CR) values for soil benzene at an abandoned oil refinery site in Northwestern China under bare soil, brick, and cement cover are 1.3 × 10-6, 1.22 × 10-6, and 9.73 × 10-7, respectively. Low-permeability covers such as clay brick and cement reduces the inhalation CR of gaseous benzene from the surface soil, and delays the growth trend of cumulative inhalation CR.

Approximate analytical model for transient transport and oxygen-limited biodegradation of vapor-phase petroleum hydrocarbon compound in soil

Volatile organic compounds (VOCs) contamination may occur in subsurface soil due to various reasons and pose great threat to people. Petroleum hydrocarbon compound (PHC) is a typical kind of VOC, which can readily biodegrade in an aerobic environment. The biodegradation of vapor-phase PHC in the vadose zone consumes oxygen in the soil, which leads to the change in aerobic and anaerobic zones but has not been studied by the existing analytical models. In this study, a one-dimensional analytical model is developed to simulate the transient diffusion and oxygen-limited biodegradation of PHC vapor in homogeneous soil. Laplace transformation and Laplace inversion of the Talbot method are adopted to derive the solution. At any given time, the thickness of aerobic zone is determined by the dichotomy method. The analytical model is verified against numerical simulation and experimental results first and parametric study is then conducted. The transient migration of PHC vapor can be divided into three stages including the pure aerobic zone stage (Stage I), aerobic-anaerobic zones co-existence stage (Stage II), and steady-state stage (Stage III). The proposed analytical model should be adopted to accommodate scenarios where the transient effect is significant (Stage II), including high source concentration, deep contaminant source, high biodegradation capacity, and high water saturation. The applicability of this model to determine the breakthrough time for better vapor intrusion assessment is also evaluated. Lower first-order biodegradation rate, higher source concentration, and shallower source depth all lead to smaller breakthrough time.

Two-dimensional analytical solution for subsurface volatile organic compounds vapor diffusion from a point source in layered unsaturated zone

Although migration of subsurface volatile organic compounds (VOCs) from contaminant sources in unsaturated soil widely exists, the related analytical models are quite limited. A two-dimensional analytical solution is hence developed to simulate vapor diffusion from the subsurface contaminant source in the layered unsaturated zone. The contaminant source is simplified as a point source leaking at a constant rate. The influences of several important factors, including thickness of stagnant air layer, depth of groundwater table, source characteristics and soil layering characteristics, on vapor migration in subsurface soil are comprehensively investigated by the present model. Soil type does not affect the normalized vapor concentration profile for homogeneous soil, which is not valid for layered soil. The width and effective diffusivity of the upward diffusion pathway and downward diffusion pathway are favorable indexes to evaluate the intensity of subsurface vapor horizontal diffusion. The single-layer capillary fringe assumption overestimates the vapor plume, the assumption can give acceptable result for coarse soil while it is recommended to divide the soil into several layers based on the water-filled porosity profile for fine soil.

Insights into vapour intrusion phenomena: Current outlook and preferential pathway scenario

Vapour intrusion (VI) is the phenomenon by which volatile organic compounds (VOCs) migrate from the subsurface source through the soil and enter into the overlying buildings, affecting the indoor air quality and ultimately causing health hazards to the occupants. Health risk assessments associated with hydrocarbon contaminated sites and recommendations of site closure are often made by quantifying the VI risks using mathematical models known as 'vapour intrusion models' (VIM). In order to predict the health risk, various factors such as the lithological and geochemical conditions of the subsurface, environmental conditions, building operational conditions etc. are commonly evaluated using VIMs. Use of these models can overlook the role of preferential pathways like highly permeable subsurface layers and utility lines which act as the path of least resistance for vapour transport, which can increase the VI risks. The extensive networks of utility lines and sanitary sewer systems in urban areas can significantly exacerbate the uncertainty of VI investigations. The backfill materials like sand and gravel surrounding the utility lines can allow the vapours to easily pass through due to their high porosity as compared to natural formations. Hence, failure to understand the role of preferential pathways on the fate and transport of VOC in the vadose zone can result in more conservative predictions of indoor air vapour concentrations and wrong clean up approaches. This comprehensive review outlines the vapour transport mechanisms, factors influencing VI, VIMs and the role of preferential pathways in predicting indoor air vapour concentrations.

Refinement of the gradient method for the estimation of natural source zone depletion at petroleum contaminated sites

Rates of natural source zone depletion (NSZD) are increasingly being used to aid remedial decision making and light non-aqueous phase liquid (LNAPL) longevity estimates at petroleum release sites. Current NSZD estimate methods, based on analyses of carbon dioxide (CO₂) and oxygen (O₂) soil-gas concentration gradients ("gradient method") assume linear concentration profiles with depth. This assumption can underestimate the concentration gradients especially above LNAPL sources that are typically characterized by curvilinear or semi-curvilinear O₂ and CO₂ concentration profiles. In this work, we proposed a new method that relies on calculating the O₂ and CO₂ concentration gradient using a first-order reaction model. The method requires an estimate of the diffusive reaction length that can be easily derived from soil-gas concentration data. A simple step-by-step guide for applying the new method is provided. Nomographs were also developed to facilitate method application. Application of the nomographs using field data from published literature showed that NSZD rates could be underestimated by nearly an order of magnitude if reactivity in the vadose zone is not accounted for. The new method helps refine NSZD rates estimation and improve risk-based decision making at certain petroleum contaminated sites.

Vapor Intrusion Investigations and Decision-Making: A Critical Review

At sites impacted by volatile organic compounds (VOCs), vapor intrusion (VI) is the pathway with the greatest potential to result in actual human exposure. Since sites with VI were first widely publicized in late 1990s, the scientific understanding of VI has evolved considerably. The VI conceptual model has been extended beyond relatively simple scenarios to include nuances, such as biological and hydrogeological factors that may limit the potential for VI and alternative pathways, such as preferential pathways and direct building contact/infiltration that may enhance VI in some cases. Regulatory guidance documents typically recommend initial concentration- or distance-based screening to evaluate whether VI may be a concern, followed by a multiple-lines-of-evidence (MLE) investigation approach for sites that do not screen out. These recommendations for detailed evaluation of VI currently focus on monitoring of VOC concentrations in groundwater, soil gas, and indoor air and can be supplemented by other lines of evidence. In this Critical Review, we summarize key elements important to VI site characterization, provide the status and current understanding, and highlight data interpretation challenges, as well as innovative tools developed to help overcome the challenges. Although there have been significant advances in the understanding of VI in the past 20 years, limitations and knowledge gaps in screening, investigation methods, and modeling approaches still exist. Potential areas for further research include improved initial screening methods that account for the site-specific role of barriers, improved understanding of preferential pathways, and systematic study of buildings and infrastructure other than single-family residences.

Development of a modular vapor intrusion model with variably saturated and non-isothermal vadose zone

Human health risk assessment at hydrocarbon-contaminated sites requires a critical evaluation of the exposure pathways of volatile organic compounds including assessments of vapor exposure in indoor air. Although there are a number of vapor intrusion models (VIM) currently available, they rarely reproduce actual properties of soils in the vadose zone. At best, users of such models assume averaged parameters for the vadose zone based on information generated elsewhere. The objective of this study was to develop a one-dimensional steady-state VIM, indoorCARE™ model, that considers vertical spatial variations of the degree of saturation (or effective air-filled porosity) and temperature of the vadose zone. The indoorCARE™ model was developed using a quasi-analytical equation that (1) solves the coupled equations governing soil-water movement driven by pressure head and a soil heat transport module describing conduction of heat and (2) provides a VIM that accommodates various types of conceptual site model (CSM) scenarios. The indoorCARE™ model is applicable to both chlorinated hydrocarbon and petroleum hydrocarbon (PHC) contaminated sites. The model incorporates biodegradations of PHCs at a range of CSM scenarios. The results demonstrate that predictions of indoor vapor concentrations made with the indoorCARE™ model are close to those of the J&E and BioVapor models under homogeneous vadose zone conditions. The newly developed model under heterogeneous vadose zone conditions demonstrated improved predictions of indoor vapor concentrations. The research study presented a more accurate and more realistic way to evaluate potential human health risks associated with the soil-vapor-to-indoor-air pathways.

Impact of potent bioremediation enhancing plant extracts on Raoultella ornithinolytica properties

Long-term contact of microorganisms with different compounds in the environment can cause significant changes in cell metabolism. Surfactants adsorption on cell surface or incorporation in the cell membrane, lead to their modification, which helps microorganisms adopt to the conditions of metabolic stress. The main objective of this study was to investigate the effects of three saponin-reach plant extracts from Hedera helix, Saponaria officinalis and Sapindus mucorossi on growth and adaptation of Raoultella ornithinolytica to high concentrations of these substances. For this purpose we investigated cell surface properties, membrane fatty acids and genetic changes of the microorganisms. The results revealed that prolonged exposure of the microorganisms to high concentrations of these surfactants can induce genetic changes of their genes. Moreover, the adaptation to contact with high concentrations of saponins was also associated with changes in composition of fatty acids responsible for the stabilisation of membrane structure and the increase in membrane permeability. The changes affected also the outer layer of cells. A significant increase (p < 0.05) in the cell surface hydrophobicity of tested strain was also observed. The cells after long-term contact with S. officinalis and S. mucorossi acquire properties that may be favourable in hydrophobic substances bioremediation.

Estimating the oxygenated zone beneath building foundations for petroleum vapor intrusion assessment

Previous studies show that aerobic biodegradation can effectively reduce hydrocarbon soil gas concentrations by orders of magnitude. Increasingly, oxygen limited biodegradation is being included in petroleum vapor intrusion (PVI) guidance for risk assessment at leaking underground storage tank sites. The application of PVI risk screening tools is aided by the knowledge of subslab oxygen conditions, which, however, are not commonly measured during site investigations. Here we introduce an algebraically explicit analytical method that can estimate oxygen conditions beneath the building slab, for PVI scenarios with impervious or pervious building foundations. Simulation results by this new model are then used to illustrate the role of site-specific conditions in determining the oxygen replenishment below the building for both scenarios. Furthermore, critical slab-width-to-source-depth ratios and critical source depths for the establishment of a subslab "oxygen shadow" (i.e. anoxic zone below the building) are provided as a function of key parameters such as vapor source concentration, effective diffusion coefficients of concrete and building depth. For impervious slab scenarios the obtained results are shown in good agreement with findings by previous studies and further support the recommendation by U.S. EPA about the inapplicability of vertical exclusion distances for scenarios involving large buildings and high source concentrations. For pervious slabs, results by this new model indicate that even relatively low effective diffusion coefficients of concrete can facilitate the oxygen transport into the subsurface below the building and create oxygenated conditions below the whole slab foundation favorable for petroleum vapor biodegradation.

Role of the source to building lateral separation distance in petroleum vapor intrusion

The adoption of source to building separation distances to screen sites that need further field investigation is becoming a common practice for the evaluation of the vapor intrusion pathway at sites contaminated by petroleum hydrocarbons. Namely, for the source to building vertical distance, the screening criteria for petroleum vapor intrusion have been deeply investigated in the recent literature and fully addressed in the recent guidelines issued by ITRC and U.S.EPA. Conversely, due to the lack of field and modeling studies, the source to building lateral distance received relatively low attention. To address this issue, in this work we present a steady-state vapor intrusion analytical model incorporating a piecewise first-order aerobic biodegradation limited by oxygen availability that accounts for lateral source to building separation. The developed model can be used to evaluate the role and relevance of lateral vapor attenuation as well as to provide a site-specific assessment of the lateral screening distances needed to attenuate vapor concentrations to risk-based values. The simulation outcomes showed to be consistent with field data and 3-D numerical modeling results reported in previous studies and, for shallow sources, with the screening criteria recommended by U.S.EPA for the vertical separation distance. Indeed, although petroleum vapors can cover maximum lateral distances up to 25-30m, as highlighted by the comparison of model outputs with field evidences of vapor migration in the subsurface, simulation results by this new model indicated that, regardless of the source concentration and depth, 6m and 7m lateral distances are sufficient to attenuate petroleum vapors below risk-based values for groundwater and soil sources, respectively. However, for deep sources (>5m) and for low to moderate source concentrations (benzene concentrations lower than 5mg/L in groundwater and 0.5mg/kg in soil) the above criteria were found extremely conservative as the model results indicated that for such scenarios the lateral screening distance may be set equal to zero.

Sensitivity and uncertainty analysis for Abreu & Johnson numerical vapor intrusion model

This study conducted one-at-a-time (OAT) sensitivity and uncertainty analysis for a numerical vapor intrusion model for nine input parameters, including soil porosity, soil moisture, soil air permeability, aerobic biodegradation rate, building depressurization, crack width, floor thickness, building volume, and indoor air exchange rate. Simulations were performed for three soil types (clay, silt, and sand), two source depths (3 and 8m), and two source concentrations (1 and 400 g/m(3)). Model sensitivity and uncertainty for shallow and high-concentration vapor sources (3m and 400 g/m(3)) are much smaller than for deep and low-concentration sources (8m and 1g/m(3)). For high-concentration sources, soil air permeability, indoor air exchange rate, and building depressurization (for high permeable soil like sand) are key contributors to model output uncertainty. For low-concentration sources, soil porosity, soil moisture, aerobic biodegradation rate and soil gas permeability are key contributors to model output uncertainty. Another important finding is that impacts of aerobic biodegradation on vapor intrusion potential of petroleum hydrocarbons are negligible when vapor source concentration is high, because of insufficient oxygen supply that limits aerobic biodegradation activities.

A two-dimensional analytical model of petroleum vapor intrusion

In this study we present an analytical solution of a two-dimensional petroleum vapor intrusion model, which incorporates a steady-state diffusion-dominated vapor transport in a homogeneous soil and piecewise first-order aerobic biodegradation limited by oxygen availability. This new model can help practitioners to easily generate two-dimensional soil gas concentration profiles for both hydrocarbons and oxygen and estimate hydrocarbon indoor air concentrations as a function of site-specific conditions such as source strength and depth, reaction rate constant, soil characteristics and building features. The soil gas concentration profiles generated by this new model are shown in good agreement with three-dimensional numerical simulations and two-dimensional measured soil gas data from a field study. This implies that for cases involving diffusion dominated soil gas transport, steady state conditions and homogenous source and soil, this analytical model can be used as a fast and easy-to-use risk screening tool by replicating the results of 3-D numerical simulations but with much less computational effort.

Evaluation of site-specific lateral inclusion zone for vapor intrusion based on an analytical approach

In 2002, U.S. EPA proposed a general buffer zone of approximately 100 feet (30 m) laterally to determine which buildings to include in vapor intrusion (VI) investigations. However, this screening distance can be threatened by factors such as extensive surface pavements. Under such circumstances, EPA recommended investigating soil vapor migration distance on a site-specific basis. To serve this purpose, we present an analytical model (AAMLPH) as an alternative to estimate lateral VI screening distances at chlorinated compound-contaminated sites. Based on a previously introduced model (AAML), AAMLPH is developed by considering the effects of impervious surface cover and soil geology heterogeneities, providing predictions consistent with the three-dimensional (3-D) numerical simulated results. By employing risk-based and contribution-based screening levels of subslab concentrations (50 and 500 μg/m(3), respectively) and source-to-subslab attenuation factor (0.001 and 0.01, respectively), AAMLPH suggests that buildings greater than 30 m from a plume boundary can still be affected by VI in the presence of any two of the three factors, which are high source vapor concentration, shallow source and significant surface cover. This finding justifies the concern that EPA has expressed about the application of the 30 m lateral separation distance in the presence of physical barriers (e.g., asphalt covers or ice) at the ground surface.

A Petroleum Vapor Intrusion Model Involving Upward Advective Soil Gas Flow Due to Methane Generation

At petroleum vapor intrusion (PVI) sites at which there is significant methane generation, upward advective soil gas transport may be observed. To evaluate the health and explosion risks that may exist under such scenarios, a one-dimensional analytical model describing these processes is introduced in this study. This new model accounts for both advective and diffusive transport in soil gas and couples this with a piecewise first-order aerobic biodegradation model, limited by oxygen availability. The predicted results from the new model are shown to be in good agreement with the simulation results obtained from a three-dimensional numerical model. These results suggest that this analytical model is suitable for describing cases involving open ground surface beyond the foundation edge, serving as the primary oxygen source. This new analytical model indicates that the major contribution of upward advection to indoor air concentration could be limited to the increase of soil gas entry rate, since the oxygen in soil might already be depleted owing to the associated high methane source vapor concentration.

Vapor intrusion screening model for the evaluation of risk-based vertical exclusion distances at petroleum contaminated sites

The key role of biodegradation in attenuating the migration of petroleum hydrocarbon vapors into the indoor environments has been deeply investigated in the last decades. Very recently, empirical screening levels for the separation distance from the source, above which the potential for vapor intrusion can be considered negligible, were defined. In this paper, an analytical solution that allows one to predict risk-based vertical screening distances for hydrocarbons compounds is presented. The proposed solution relies on a 1-D vapor intrusion model that incorporates a piecewise first-order aerobic biodegradation limited by oxygen availability and accounts also for the effect of the building footprint. The model predictions are shown to be consistent with the results obtained using a 3-D numerical model and with the empirical screening criteria defined by U.S.EPA and CRC care. However, the different simulations carried out show that in some specific cases (e.g., large building footprint, high methane concentration, and low attenuation in the capillary fringe), the respect of these empirical screening criteria could be insufficient to guarantee soil-gas concentrations below acceptable risk-based levels.

Estimation of contaminant subslab concentration in petroleum vapor intrusion

In this study, the development and partial validation are presented for an analytical approximation method for prediction of subslab contaminant concentrations in PVI. The method involves combining an analytic approximation to soil vapor transport with a piecewise first-order biodegradation model (together called the Analytic Approximation Method, including Biodegradation, AAMB), the result of which calculation provides an estimate of contaminant subslab concentrations, independent of building operation conditions. Comparisons with three-dimensional (3-D) simulations and another PVI screening tool, BioVapor, show that the AAMB is suitable for application in a scenario involving a building with an impermeable foundation surrounded by open ground surface, where the atmosphere is regarded as the primary oxygen source. Predictions from the AAMB can be used to determine the required vertical source-building separation, given a subslab screening concentration, allowing identification of buildings at risk for PVI. This equation shows that the "vertical screening distance" suggested by U.S. EPA is sufficient in most cases, as long as the total petroleum hydrocarbon (TPH) soil gas concentration at the vapor source does not exceed 50-100mg/L. When the TPH soil gas concentration of the vapor source approaches a typical limit, i.e. 400mg/L, the "vertical screening distance" required would be much greater.

Numerical model investigation for potential methane explosion and benzene vapor intrusion associated with high-ethanol blend releases

Ethanol-blended fuel releases usually stimulate methanogenesis in the subsurface, which could pose an explosion risk if methane accumulates in a confined space above the ground where ignitable conditions exist. Ethanol-derived methane may also increase the vapor intrusion potential of toxic fuel hydrocarbons by stimulating the depletion of oxygen by methanotrophs, and thus inhibiting aerobic biodegradation of hydrocarbon vapors. To assess these processes, a three-dimensional numerical vapor intrusion model was used to simulate the degradation, migration, and intrusion pathway of methane and benzene under different site conditions. Simulations show that methane is unlikely to build up to pose an explosion hazard (5% v/v) if diffusion is the only mass transport mechanism through the deeper vadose zone. However, if methanogenic activity near the source zone is sufficiently high to cause advective gas transport, then the methane indoor concentration may exceed the flammable threshold under simulated conditions. During subsurface migration, methane biodegradation could consume soil oxygen that would otherwise be available to support hydrocarbon degradation, and increase the vapor intrusion potential for benzene. Vapor intrusion would also be exacerbated if methanogenic activity results in sufficiently high pressure to cause advective gas transport in the unsaturated zone. Overall, our simulations show that current approaches to manage the vapor intrusion risk for conventional fuel released might need to be modified when dealing with some high ethanol blend fuel (i.e., E20 up to E95) releases.

A numerical investigation of oxygen concentration dependence on biodegradation rate laws in vapor intrusion

In subsurface vapor intrusion, aerobic biodegradation has been considered as a major environmental factor that determines the soil gas concentration attenuation factors for contaminants such as petroleum hydrocarbons. The site investigation has shown that oxygen can play an important role in this biodegradation rate, and this paper explores the influence of oxygen concentration on biodegradation reactions included in vapor intrusion (VI) models. Two different three dimensional (3-D) numerical models of vapor intrusion were explored for their sensitivity to the form of the biodegradation rate law. A second order biodegradation rate law, explicitly including oxygen concentration dependence, was introduced into one model. The results indicate that the aerobic/anoxic interface depth is determined by the ratio of contaminant source vapor to atmospheric oxygen concentration, and that the contaminant concentration profile in the aerobic zone was significantly influenced by the choice of rate law.

A review of vapor intrusion models

A complete vapor intrusion (VI) model, describing vapor entry of volatile organic chemicals (VOCs) into buildings located on contaminated sites, generally consists of two main parts: one part describing vapor transport in the soil and the other describing its entry into the building. Modeling the soil vapor transport part involves either analytically or numerically solving the equations of vapor advection and diffusion in the subsurface. Contaminant biodegradation must often also be included in this simulation, and can increase the difficulty of obtaining a solution, especially when explicitly considering coupled oxygen transport and consumption. The models of contaminant building entry pathway are often coupled to calculations of indoor air contaminant concentration, and both are influenced by building construction and operational features. The description of entry pathway involves consideration of building foundation characteristics, while calculation of indoor air contaminant levels requires characterization of building enclosed space and air exchange within this. This review summarizes existing VI models, and discusses the limits of current screening tools commonly used in this field.

Role of natural attenuation in modeling the leaching of contaminants in the risk analysis framework

Natural attenuation (NA) processes occurring in the subsurface can significantly affect the impact on groundwater from contamination sources located in the vadose zone, especially when mobile and readily biodegradable compounds, such as BTEX, are present. Besides, in the last decades several studies have shown natural attenuation to take place also for more persistent compounds, such as Polycyclic Aromatic Hydrocarbons (PAHs). Nevertheless, common risk analysis frameworks, based on the ASTM RBCA (Risk Based Corrective Action) approach, do not include NA pathways in the fate and transport models, thus possibly leading to an overestimation of the calculated risk. The aim of this study was to provide an insight on the relevance of the different key natural attenuation processes usually taking place in the subsurface and to highlight for which contamination scenarios their inclusion in the risk-analysis framework could provide a more realistic risk assessment. To this end, an analytical model accounting for source depletion and biodegradation, dispersion and diffusion during leaching was developed and applied to several contamination scenarios. These scenarios included contamination by BTEX, characterized by relatively high mobility and biodegradation rate, and PAHs, i.e. a more persistent class of compounds. The obtained results showed that BTEX are likely to be attenuated in the source zone due to their mobility and ready biodegradation (assuming biodegradation constant rates in the order of 0.01-1 d(-1)). Instead, attenuation along transport through the vadose zone was found to be less important, as the residence time of the contaminant in the unsaturated zone is often too low with respect to the time required to get a relevant biodegradation of BTEX. On the other hand, heavier compounds such as PAHs, were found to be attenuated during leaching since the residence time in the vadose zone can reach values up to thousands of years. In these cases, even with the relatively slow biodegradation rate of PAHs, in the order of 0.0001-0.001 d(-1), attenuation can result significant. These conclusions were also confirmed by comparing the model results with experimental data collected at an hydrocarbon-contaminated site. The proposed model, that neglects the transport of NAPLs, could be easily included in the risk-analysis framework, allowing to get a more realistic assessment of risks, while keeping the intrinsic simplicity of the ASTM-RBCA approach.