19 research outputs found
methodological approach for the estimation of recoverable lnapl volumes in aquifer and for the analysis of residual phase behaviour in laboratory
Free and residual phases are the most abundant phases of LNAPL (light non aqueous phase liquid) in the subsoil; therefore the free phase volume estimation and the understanding of residual phase behaviour are fundamental for an effective remediation. The volume estimation of free phase present in a contaminated site was carried out through the application of the Pancake Model and the Vertical Equilibrium Model. This estimation shows a remarkable difference between the two models and between the different delimitation area methods employed. In particular, the estimated volumes are lower for the Vertical Equilibrium Model and the higher difference is observed for a 200x200 mesh. The results underline also some critical points as the amount of the effective porosity; indeed the sensitivity analysis shows that a reduction of this parameter produce a variation of the estimated volume until to 20%. The behaviour of residual phase was analysed through lab-scale column tests carried out using three different porous media and toluene as contaminant. The higher residual saturation and lower dissolution in finer materials reveal the influence of porous media on toluene behaviour. The lower dissolution is confirmed also by the modelling realized using a traditional approach and an approach based on experimental results. The risk analysis, carried out applying a traditional approach and a new experimental approach, shows that the last one allow to have a more reliable hazard index which contemplates the site specific conditions
Approccio metodologico per la stima dei volumi di LNAPL recuperabili in falda e per l’analisi del comportamento della fase residuale in laboratorio
NAPLs (non aqueous phase liquids), for their chemical-physical properties, are the most common and harmful contaminants in groundwater as they represent potential long-term sources of contamination and are carcinogenic or toxic for human health (Baciocchi et al., 2010). Based on their relative density with respect to water, NAPLs are classified as DNAPLs (denser than water) and LNAPLs (lighter than water). The former ones are the subject of this study. In the vadose zone, immediately following the release, LNAPLs typically migrate downwards under the influence of gravity. This vertical migration continues only if the capillary forces exceed the residual soil retention capacity (CL:AIRE, 2014). If this happens, LNAPLs will continue to migrate downwards until the water table is reached. Once the LNAPL is in contact with the capillary fringe, the contaminant starts spreading laterally (Baldi and Pacciani, 1997; Brost and DeVaull, 2000), unless sufficient LNAPL potential energy exists for it to displace water and penetrate the water table (CL:AIRE, 2014). LNAPL can enter into the pores and dislocate water only if it reaches and exceeds the entry pressure, which is inversely proportional to the pore throat radius. The entry pressure will be greater, for smaller radius values, due to the inverse relationship between this pressure and pore throat radius. Hence, LNAPL will enter more easily in the large pore. The entry pressure and capillary forces are two of the key factors that control the behaviour of LNAPL and its partitioning in the residual phase (adsorption to solids particles due to capillary forces) and free phase (a liquid separated phase immiscible in water). Conversely, the partitioning in vapour (in the vadose zone) and dissolved phase (which forms plumes of contaminants in the phreatic zone) is controlled by Henry’s constant and solubility (Dippenaar et al., 2005; Jeong and Charbeneau, 2014; Pankow and Cherry, 1996). The free and residual phase represent the 99% of LNAPL in the subsoil and are difficult to detect and quantify, so they are subject of this study in order also to understand the fate and impact of these contaminants on drinkable water quality.
In order to reach this focus, a contaminated site, located in Sicily (Southern Italy), which is characterized by a diffused contamination by LNAPLs in free and dissolved phase, chlorinated solvents, agricultural fertilizers and heavy metals, was used as a tool to estimate the volume of the free phase. In particular, the volume of the free phase present in the site has been quantified using two different conceptual models: the Pancake Model and the Vertical Equilibrium Model.
According to the Pancake Model, the migration of LNAPL to the water table and its lateral spreading through the capillary fringe creates a buoyant pool with uniform and constant saturation (Baldi and Pacciani, 1997; CL:AIRE, 2014). Since hydrocarbons are assumed immiscible in water, the free phase is suspended on the capillary fringe above the water table and so the thickness measured in the monitoring well is an apparent thickness (Dippenaar et al., 2005; Gruszczenski, 1987; Testa and Paczkowski, 1989). The difference between apparent and real thickness is due to the absence of the capillary fringe in the monitoring well. In fact, this absence conducts to a lower level of water table in the well and so free product flows more easily in the monitoring well creating an exaggerated thickness of supernatant (Hughes et al., 1988). In addition, the weight of free phase depresses above the water table in the monitoring well facilitating the flow of free product in the well. Since according to the Pancake Model there is a difference between the apparent and the real thickness, it is necessary to correct the thickness measured in the well. There are different ways to correct this measure such as empirical factors (e.g. four according to De Pastrovich et al., 1979), or factors derived by field test (recharge test and baildown test). In the studied contaminated site, baildown tests have been used to calculate the exaggeration factor in order to estimate the real thickness of free phase in the aquifer. The baildown test consists of pumping only the supernatant (Hughes et al, 1988) or supernatant and water (Gruszczenski, 1987) from the well. Initially, before the pumping begins, it is necessary to measure the apparent thickness of the product in the well. During the test, the level of air/LNAPL interface and LNAPL/water interface are measured and their difference provides the free product thickness. The real thickness is considered as the measured thickness corresponding to the inflection point in the graph LNAPL/water interface vs time.
The volume of the free product existing in the site at June 2013, has been calculated using the Pancake Model and Thiessen polygons and grid at regular square mesh (100 m x 100 m and 200 m x 200 m). First, we have identified the monitoring wells with supernatant and the measured apparent thickness, then we have subdivided the monitoring wells based on the type of free product detected (diesel with more of 70% of C10-C30, gasoline with more of 70% of C6-C9 and mixtures of gasoline and diesel). Based on this first evaluation, different exaggeration factors, obtained by baildown tests, have been used to calculate the real thickness. Once obtained the real thickness for every monitoring well, the specific volume has been calculated as the product of the real thickness and the effective porosity (that is 0.25 for sands and 0.20 for silty sands). The specific volume has been then multiplied for the area associated at each monitoring point using the Thiessen polygons and the grid at regular mesh. The total volume of the free product obtained using the Thiessen polygons was about 9000 m3, while using the grid at 200x200 m mesh and 100x100 m mesh the volumes were about 8700 m3 and 4800 m3, respectively.
The Vertical Equilibrium Model instead, can been used only when there is not persistence of release of contaminant (Lundegard and Mudford, 1998), i.e. the scenario considered in the study site. In such case, it is assumed that there is not a discrete layer of LNAPL floating on the water table, but that LNAPL can penetrate below the water table. In addition, the pore fraction occupied by product is less than 100% due to the presence of other fluids such as air and water, and the LNAPL saturation varies with the depth (ITRC, 2009; Lundegard and Mudford, 1998). The shape of LNAPL saturation curves can be regular (shark fin) if there is a homogeneous aquifer or irregular if the aquifer is heterogeneous. Hence, according to the Vertical Equilibrium Model, the relationship between the LNAPL thickness in the monitoring well and its specific volume in the aquifer is linked by the capillary properties of the soil and the LNAPL characteristics. Usually, given a certain thickness in the well, the specific volume will be lower if a finer material constitutes the aquifer (Lundegard and Mudford, 1998). Therefore, the Vertical Equilibrium Model requires the knowledge of the saturation profiles to calculate the specific volume. These curves can be obtained in different ways. For instance, they can estimated based on empirical approaches (measuring saturation in core samples), by analytical modelling or with software as LDRM (LNAPL Distribution and Recovery Model, distributed by American Petroleum Institute) which requires information about the characteristics of the LNAPL and the affected aquifer.
For the calculation of the LNAPL volume existing in the site at June 2013, we used the Vertical Equilibrium Model by employing the LDRM approach which, based on the input information, provides the specific volume (Dn) and the recoverable specific volume (Rn). Then, also in this case we have used the Thiessen polygons and the grid at regular mesh to estimate the total volume. The results of the calculation are about 5700 m3 for the Thiessen polygons, and about 4300 m3 and 3000 m3 for 200x200 and 100x100 grid at regular mesh, respectively. The recoverable volume calculation shows that about 70% of the product present in the site can be removed by pumping.
The comparisons between the results obtained using the two models and different methods to calculate areas, show that there is a difference of thousands of cubic meter between the estimated volumes. In particular, the Vertical Equilibrium Model estimates lower volumes than the Pancake Model; this is probably due to the LNAPL saturation considered in the two models (100% in Pancake Model and less of 100% in Vertical Equilibrium Model). In addition, it was found that for both models the estimated volume was higher for 200x200 regular mesh than for 100x100 regular mesh. The difference between volumes estimated with the Pancake Model and the Vertical Equilibrium Model was about 51% using 200x200 mesh and about 37% using 100x100 mesh. This result allows to suppose that the use of smaller mesh can lead at a reduction of the differences, but at this stage the available data do not permit to use smaller mesh and confirm this theory. The use of the Thiessen polygons, that show a difference of 37% between volumes estimated with the two different models, reveals another problem related to the construction method of these polygons; in fact, since they are created as a function of the distance and spatial distribution of the monitoring wells, they have different shape and dimension and this influences the volume calculation. Monitoring wells with similar specific volume can be associated to polygons with areas much different and so the volume calculated for one well can be very different from that of another well. The comparison of volumes obtained by the Pancake Model considering only the piezometers and those obtained considering the piezometers and the wells confirms the influence of the Thiessen polygons on the volume estimation. In fact, in the first calculation, the volume of supernatant was about 9000 m3, while in the second one instead it was about 6800 m3.
In addition to the uncertainty due to the area delimitation, other critical points in the calculation of the LNAPL volume have been encountered, such as 1) measure of the product thickness in the wells, because it can be influenced by the removal of the product that can be done the days before the measurement; 2) presence of different types of product; 3) lack of some site-specific data as α and N parameters of Van Genuchten (API, 2001), the porosity, the irreducible saturation of water and the interfacial tensions.
A sensitivity analysis has been done about the porosity and it has revealed that a reduction of this parameter (from 0.25 to 0.20 for sands and from 0.20 to 0.15 for silty sands) provides a reduction of about 20-22% of specific volume for sands and of about 25% for silty sands.
The recovery of LNAPL volumes present in the site is ongoing through the existing and new recovery wells. In fact, to date, the realization of 8 new wells and 3 piezometers for the LNAPL recovery have been foreseen. These new wells will accelerate the recovery of free LNAPL until the residual saturation will be reached, i.e. when the remaining volume of LNAPL will not be recoverable and it will remain as residual phase.
The residual phase has been the focus of the second part of this work. Indeed, in order to accurately plan and realize an effective remediation coupled with a good cost/benefits ratio, it is necessary to understand its behaviour and interaction with the impacted soil and groundwater. To this end, we have carried out lab-scale column tests using different porous media (glass spheres, sandy soil A with 9% of silt and 1% of clay, sandy soil B with 14% of silt and 2% of clay) and toluene as contaminant. The cylindrical column (12.7 cm x 2.9 cm) was packed, from bottom to top, with 2 cm of large glass spheres (φ=5 mm), 1 cm of small glass spheres (φ=3 mm), 2 glass-microfibre filters (0.7 μm), 6.7 cm of soil (dried in an oven at 110°C) corresponding to approximately 74 g of soil (glass spheres φ=6 mm in the control test), 2 glass-microfibre filters (0.7 μm), 1 cm of small glass spheres and 2 cm of large glass spheres. The column was first saturated, using a peristaltic pump, with demineralized water from the bottom to remove gas bubbles trapped in the porous media and to estimate (based on the water trapped in the column) the pore volume (PV) of the column. Then, toluene was fed into the column to simulate the movement of organic liquid into the saturated zone (U.S. EPA, 1990) until the water present in the column was completely displaced by toluene, indicating complete saturation with the organic liquid (Powers et al.,1992; US. EPA, 1990). The column was then flushed from the bottom with demineralized water at relatively flow rates (from 0.5 to 1 ml/min) for approximately 25 PV to displace the free phase of toluene. After this step, it was possible to calculate the residual saturation, i.e. the saturation at which the NAPL becomes discontinuous and it is immobilized by capillary forces under ambient groundwater flow conditions (Mercer and Cohen, 1990). The residual saturation has been calculated as the ratio between volume of toluene remained in the soil and volume of toluene injected in the soil. Successively, to assess the dissolution kinetics of toluene from the residual phase entrapped in the column, demineralized water was fed into the column and effluent samples were analysed for toluene by static headspace gas chromatography-mass spectrometry (HSS-GC-MS) using fluorobenzene as internal standard. The results of these experiments showed that the residual saturation is function of the grain size distribution; in fact, residual saturations for glass spheres, soil A and soil B were respectively 36%, 70% and 77%. These data showed that residual saturation is inversely proportional to grain size dimension, indeed it increases with decreasing grain size. Dissolution profiles showed that the solubility of toluene in water was never reached during the test, except in some cases in the first samples collected by glass sphere test and soil A test. Conversely, in the soil B test, the solubility was never reached. In addition, it has been observed that the concentration of toluene decreases more quickly with the increase of particle size. This is due to capillary forces that in the finer soil are higher and so hold more toluene conducting to higher residual saturation and lower dissolution of residual phase. Experimental results are also been modelled using a first order kinetic model according to which the cumulative mass released by the soil is:
Mout (cum, measured) = Σ Cmeasured,iΔVi (1)
where Cmeasured (mg/l) is the toluene concentration measured at the column outlet at the i-th interval and ΔV(l) is the volume of water fed into the column during the i-th interval.
This mass was compared with the mass calculated following the traditional approach that considers the concentration of toluene in water as the solubility:
Mout (cum, equilibrium) = Σ Csol * ΔV= Csol * npv * PV (2)
where Csol (mg/l) is the solubility of the contaminant in the water, ΔV (l) is the volume of water fed into the column during the i-th interval, npv (-) is the number of pore volumes fed into the column and PV (l) is the pore volume of the column.
By combining eq. (1) and (2), it is possible to estimate the ratio (CF) between theoretical cumulative mass assuming saturation condition (eq. 2) and the experimental cumulative mass (eq. 1), both calculated for a specific number of pore volumes delivered (npv*):
CF (mpv*) =Mout (cum, equilibrium; npv*)/ Mout (cum; npv*) (3)
The higher is the CF values obtained through eq. (3), the greater will be the overestimation of the mass release calculated assuming equilibrium conditions between the eluate and the NAPL. The results of modelling show that there is a good correlation between experimental results and the first order model. The assumption according to which the eluate from the residual LNAPL is saturated with toluene may lead to an overestimation of the true dissolved toluene concentration in water.
Finally, the experimental results were used also to carry out a risk analysis. In particular, an effective hazard index is directly proportional to hazard index (HI) estimated with traditional approach (ASTM-Risk-based Corrective Action) and inversely proportional to the correction factor (CF). The results of risk analysis showed that the traditional ASTM-RBCA approach can lead to an overestimation of the hazard index for human health and this is particularly significant for sources characterized by lower lengths along the groundwater direction and in scenarios characterized by higher groundwater velocity. Based on these findings, it was also found that the approach proposed in some countries (Carlon, 2007) to estimate the risk to groundwater resources based on conservative concentration thresholds may provide an inaccurate perception of the effective impact on groundwater quality. Indeed, the use of steady-state transport models combined with linear equilibrium partitioning model results in a situation that is representative only of the very early stages after the contamination event. The results indicated that after a relatively short time, the contaminant concentrations in groundwater are expected to decrease by orders of magnitude. This behaviour cannot be properly evaluated by a simple comparison with threshold values, whereas it can be easily accounted for by adopting a risk assessment procedure (especially if coupled with the approach proposed in this work) that calculates the risk based on the effective cumulative dose over the entire exposure duration
Monitoring of composition of LNAPL: essential tool for the estimation of free LNAPL specific volumes
Light Non-Aqueous Phase Liquid, Pancake Model, Vertical Equilibrium Mode
Arsenic-fluoride co-contamination in groundwater: Background and anomalies in a volcanic-sedimentary aquifer in central Italy
Abstract In a volcanic-sedimentary aquifer in central Italy, we investigate the co-existence of arsenic and fluoride in groundwater, aiming at identifying the most probable processes deductible at regional/groundwater body scale leading to the observed co-contamination in groundwater. Further, the areas at risk for human health where high concentrations can produce a significant risk to human health have been investigated. The study area is located in Latium (Central Italy) where silica-undersaturated alkali-potassic formations of Plio-Pleistocene age largely outcrop above marine and continental sand and clay deposits (Neogene) and continental alluvial deposits (Lower Pleistocene–Middle Pleistocene). Geochemical data from groundwater at 322 wells and 76 springs have been analyzed through statistical methods including clustering/PCA and geostatistical analysis. The results show exceedances of the drinking water standards for F (1.5 mg/L) and As (10.0 μg/L) in 29% and 55% of the sampled groundwater, respectively. Multivariate statistics suggest a widespread process of water-rock interaction with the K-alkaline volcanic formations releasing As, F, K, Si, V, Rb and PO4 to the groundwater. As and F show a good correlation (Pearson's r = 0.61, Spearman's rs = 0.59) and define a separate PCA component, confirming that their background in groundwater might be governed by a common process. Kriging interpolations have been used to study the spatial distribution of the two parameters, identifying areas with the highest concentrations and highest probability of exceeding the standards for human consumption. Moreover, by resampling the As-F data with the jackknife technique it was possible to identify the variations of their correlation index in the study area, due to specific As or F anomalies. While in the peripheral areas of the volcanic districts, dominated by sedimentary deposits, the As-F correlation index does not present important fluctuations, Indicator Kriging shows specific As or F correlation anomalies within the volcanic groundwater bodies and along the Tyrrhenian coastline. These anomalies seem to correspond to the zones with the highest thermal flux and/or are located near important structural lineaments. Fluoride correlation anomalies close to mining sites (fluorite) have also been observed. We hypothesize that, unlike the regional co-contamination, these local anomalies are related to the upwelling of geothermal fluids along fracture/fault systems that mix with cold groundwater, or to the interaction with mineral deposits particularly enriched of these elements
Proposal for an integrated method of natural background levels assessment in groundwater
Natural Background Level (NBL), geochemical approach, statistical methods
A proposal for groundwater sampling guidelines: application to a case study in Southern Latium
Groundwater monitoring , body status, best sampling techniques, quaternary alluvial, pyroclastic deposit
Physiological Profiling and Functional Diversity of Groundwater Microbial Communities in a Municipal Solid Waste Landfill Area
The disposal of municipal solid wastes in landfills represents a major threat for aquifer environments at the global scale. The aim of this study was to explore how groundwater geochemical characteristics can influence the microbial community functioning and the potential degradation patterns of selected organic substrates in response to different levels of landfill-induced alterations. Groundwaters collected from a landfill area were monitored by assessing major physical-chemical parameters and the microbiological contamination levels (total coliforms and fecal indicators—Colilert-18). The aquatic microbial community was further characterized by flow cytometry and Biolog EcoPlatesTM assay. Three groundwater conditions (i.e., pristine, mixed, and altered) were identified according to their distinct geochemical profiles. The altered groundwaters showed relatively higher values of organic matter concentration and total cell counts, along with the presence of fecal indicator bacteria, in comparison to samples from pristine and mixed conditions. The kinetic profiles of the Biolog substrate degradation showed that the microbial community thriving in altered conditions was relatively more efficient in metabolizing a larger number of organic substrates, including those with complex molecular structures. We concluded that the assessment of physiological profiling and functional diversity at the microbial community level could represent a supportive tool to understand the potential consequences of the organic contamination of impacted aquifers, thus complementing the current strategies for groundwater management
Methodological approach for the estimation of recoverable LNAPL volumes in aquifer and for the analysis of residual phase behaviour in laboratory
Free and residual phases are the most abundant phases of LNAPL (light non aqueous phase liquid) in the subsoil; therefore the free phase volume estimation and the understanding of residual phase behaviour are fundamental for an effective remediation. The volume estimation of free phase present in a contaminated site was carried out through the application of the Pancake Model and the Vertical Equilibrium Model. This estimation shows a remarkable difference between the two models and between the different delimitation area methods employed. In particular, the estimated volumes are lower for the Vertical Equilibrium Model and the higher difference is observed for a 200x200 mesh. The results underline also some critical points as the amount of the effective porosity; indeed the sensitivity analysis shows that a reduction of this parameter produce a variation of the estimated volume until to 20%. The behaviour of residual phase was analysed through lab-scale column tests carried out using three different porous media and toluene as contaminant. The higher residual saturation and lower dissolution in finer materials reveal the influence of porous media on toluene behaviour. The lower dissolution is confirmed also by the modelling realized using a traditional approach and an approach based on experimental results. The risk analysis, carried out applying a traditional approach and a new experimental approach, shows that the last one allow to have a more reliable hazard index which contemplates the site specific conditions
Free LNAPL Volume Estimation by Pancake Model and Vertical Equilibrium Model: Comparison of Results, Limitations, and Critical Points
Light nonaqueous phase liquids (LNAPLs), due to their low solubility, dissolve slowly, acting as a long-term source of water contamination, and consequently they represent an important environmental issue. In the subsoil, more than 99% of spilled LNAPL remains as adsorbed and free phase; therefore, the volume estimation of free phase, obtained in this case through two different conceptual models (Pancake Model and Vertical Equilibrium Model), is considered a fundamental step for a correct site remediation. According to the first model, the LNAPL floating on the water table and its saturation is up to 100%; instead, according to the second one, the LNAPL can penetrate below the water table and the coexistence of LNAPL, water, and air in the pore fraction, leads to a lower LNAPL saturation, variable with the depth. Actually, in subsoil LNAPL and water saturations vary with depth due to the influence of capillarity, leading to the inaccuracy of Pancake Model assumption. Despite the evident limitation of Pancake Model, both models were applied, coupled with area calculations with Thiessen polygons and grid at regular mesh, to roughly estimate the free LNAPL volume existing in a contaminated site. The volume estimation carried out, considering the LNAPL type and its features, the soil type, and relative effective porosity, provides estimates of volumes having differences up to thousands of cubic meters. The results analysis shows that this estimation has several critical points such as area definition and the lack of site-specific data (e.g., porosity). Indeed, the sensitivity analysis for porosity shows that a reduction of this parameter provides a 20% reduction of estimated volume