Nadim K. Copty

Stochastic analysis of pumping test drawdown data in heterogeneous geologic formations

Pumping tests are often used to estimate effective values of the transmissivity and the storativity of the perturbed aquifer portion surrounding the pumping well and observation points. The interpretation of such pumping tests is based on analytic or semi-analytic methods that are normally developed under the assumption of homogeneity in the vicinity of the well. The purpose of this paper is to examine the impact of the local-scale heterogeneity of the transmissivity on the transient drawdown due to pumping. The log-transmissivity field is modeled as a multi-variate Gaussian random spatial function. Using Monte Carlo simulations, the aquifer response to a pumping test is simulated for various parameter values of the log-transmissivity variance and integral scale. The transient drawdown of the heterogeneous system is compared to the drawdown of an equivalent homogeneous aquifer with transmissivity equal to the effective transmissivity. The results of this study indicate that the impact of local heterogeneity expressed by the log-transmissivity variance is most pronounced in the early-time drawdown rates at or near the pumping well. The integral scale of the log-transmissivity on the other hand influences the time needed for the drawdown rate of the heterogeneous system to approach that of the equivalent homogeneous system.

A New Method for the Interpretation of Pumping Tests in Leaky Aquifers

A novel methodology for the interpretation of pumping tests in leaky aquifer systems, referred to as the double inflection point (DIP) method, is presented. The method is based on the analysis of the first and second derivatives of the drawdown with respect to log time for the estimation of the flow parameters. Like commonly used analysis procedures, such as the type-curve approach developed by Walton (1962) and the inflection point method developed by Hantush (1956), the mathematical development of the DIP method is based on the assumption of homogeneity of the leaky aquifer layers. However, contrary to the two methods developed by Hantush and Walton, the new method does not need any fitting process. In homogeneous media, the two classic methods and the one proposed here provide exact results for transmissivity, storativity, and leakage factor when aquifer storage is neglected and the recharging aquifer is unperturbed. The real advantage of the DIP method comes when applying all methods independently to a test in a heterogeneous aquifer, where each method yields parameter values that are weighted differently, and thus each method provides different information about the heterogeneity distribution. Therefore, the methods are complementary and not competitive. In particular, the combination of the DIP method and Hantush method is shown to lead to the identification of contrasts between the local transmissivity in the vicinity of the well and the equivalent transmissivity of the perturbed aquifer volume.

Impact of dilution on the transport of poly(acrylic acid) supported magnetite nanoparticles in porous media

This paper investigates the impact of dilution on the mobility of magnetite nanoparticles surface coated with poly(acrylic acid) (PAA). Transport experiments were conducted in a water-saturated sand-packed column for input nanoparticle solutions with total Fe concentrations ranging from 100 to 600mg/L. Particle size analysis of the synthesized nanoparticle solutions showed that PAA provides good size stability for Fe concentrations as low as about 1mg/L. Time-moment analysis of the nanoparticle breakthrough curves, on the other hand, revealed that nanoparticle mass recovery from the column decreased consistently with dilution, with greater attenuation, sharper fronts and longer tails compared to that of the tracer. Particle size analysis of the eluted solutions shows that the nanoparticle size is negatively correlated with nanoparticle concentration. Modeling results suggest that the decrease in nanoparticle mobility with input concentration can be represented using a kinetic time-dependent deposition term with finite deposition capacity and a kinetic detachment term. For field applications, the increase in particle size and detachment resulting from dilution means reduced transport efficiency of nanoparticles and reaction potential with travel distance.

Leaching potential of nano-scale titanium dioxide in fresh municipal solid waste.

With the rapid development in nanotechnology in recent years, the number of commercially available products containing engineered nanomaterials (ENMs) has increased significantly. It is expected that large fractions of these ENMs will end up in landfills for final disposal. Despite the wide use of ENMs, little data is available on their fate within landfills. This study examined the leaching behavior of nanoscale titanium dioxide (nano-TiO2), one of the mostly used ENMs, in fresh municipal solid wastes (MSWs). Batch reactors containing municipal waste samples were spiked with a range of nano-TiO2 concentrations at different pH and ionic strength conditions. The Ti concentrations in leachate decreased rapidly and reached steady state after about 12-24 h. Results suggest that, for the environmental conditions considered, approximately 3-19% of the added nano-TiO2 remained in leachate. Batch tests conducted with individual synthetically-prepared solid waste components also showed low leaching potential (5.2% for organic waste, 3.3% for glass, 1.7% for both textile and paper and 0.6% for metal), indicating that all components of MSW contributed to the retention of the nano-TiO2 mass within the solid matrix.

Estimating transmissivity from single‐well pumping tests in heterogeneous aquifers

Although aquifers are naturally heterogeneous, the interpretation of pumping tests is commonly performed under the assumption of aquifer homogeneity. This yields interpreted hydraulic parameters averaged over a domain of uncertain extent which disguises their relation to the underlying heterogeneity. In this study, we numerically investigate the sensitivity of the transient drawdown at the pumping well, to nonuniform distributions of transmissivity in confined aquifers. Frechet kernels and their time derivative are used to estimate two spatially averaged transmissivities, denoted the equivalent and interpreted transmissivity, Teq and Tin, respectively, for the case of single-well pumping tests. Interrelating Teq and Tin is achieved by modeling Tin in terms of a distance dependent, radially heterogeneous field. In weakly heterogeneous aquifers, Teq approximates TPW, the local transmissivity at the pumped well. With increasing degree of heterogeneity, Teq deviates from TPW as pumping propagates. Tin starts at TPW, approaching the spatial geometric mean of transmissivity during late pumping times. Limits of the proposed spatial weighting functions are investigated by treating the interpreted storativity, Sest, as an indicator for flow connectivity. It is shown numerically that the spatial weights for Teq and Tin agree well to the underlying heterogeneity if Sest<1. Finally, implications for applying the concepts of Teq and Tin to heterogeneous domains, and, for real world applications are discussed. It is found that time-dependent spatial averages of Tin agree well with estimates of the interpreted transmissivity from the Continuous-Derivation method.

Laboratory‐scale experiments and numerical modeling of cosolvent flushing of multi-component NAPLs in saturated porous media

This study examines the mechanistic processes governing multiphase flow of a water-cosolvent-NAPL system in saturated porous media. Laboratory batch and column flushing experiments were conducted to determine the equilibrium properties of pure NAPL and synthetically prepared NAPL mixtures as well as NAPL recovery mechanisms for different water-ethanol contents. The effect of contact time was investigated by considering different steady and intermittent flow velocities. A modified version of multiphase flow simulator (UTCHEM) was used to compare the multiphase model simulations with the column experiment results. The effect of employing different grid geometries (1D, 2D, 3D), heterogeneity and different initial NAPL saturation configurations was also examined in the model. It is shown that the change in velocity affects the mass transfer rate between phases as well as the ultimate NAPL recovery percentage. The experiments with low flow rate flushing of pure NAPL and the 3D UTCHEM simulations gave similar effluent concentrations and NAPL cumulative recoveries. Model simulations over-estimated NAPL recovery for high specific discharges and rate-limited mass transfer, suggesting a constant mass transfer coefficient for the entire flushing experiment may not be valid. When multi-component NAPLs are present, the dissolution rate of individual organic compounds (namely, toluene and benzene) into the ethanol-water flushing solution is found not to correlate with their equilibrium solubility values.

Impact of deforestation on soil carbon stock and its spatial distribution in the Western Black Sea Region of Turkey

Land use management is one of the most critical factors influencing soil carbon storage and the global carbon cycle. This study evaluates the impact of land use change on the soil carbon stock in the Karasu region of Turkey which in the last two decades has undergone substantial deforestation to expand hazelnut plantations. Analysis of seasonal soil data indicated that the carbon content decreased rapidly with depth for both land uses. Statistical analyses indicated that the difference between the surface carbon stock (defined over 0-5 cm depth) in agricultural and forested areas is statistically significant (Agricultural = 1.74 kg/m(2), Forested = 2.09 kg/m(2), p = 0.014). On the other hand, the average carbon stocks estimated over the 0-1 m depth were 12.36 and 12.12 kg/m(2) in forested and agricultural soils, respectively. The carbon stock (defined over 1 m depth) in the two land uses were not significantly different which is attributed in part to the negative correlation between carbon stock and bulk density (-0.353, p < 0.01). The soil carbon stock over the entire study area was mapped using a conditional kriging approach which jointly uses the collected soil carbon data and satellite-based land use images. Based on the kriging map, the spatially soil carbon stock (0-1 m dept) ranged about 2 kg/m(2) in highly developed areas to more than 23 kg/m(2) in intensively cultivated areas as well as the averaged soil carbon stock (0-1 m depth) was estimated as 10.4 kg/m(2).

Leaching of nano-ZnO in municipal solid waste.

Despite widespread use of engineered nanomaterials (ENMs) in commercial products and their potential disposal in landfills, the fate of ENMs in solid waste environments are still not well understood. In this study, the leaching behavior of nano ZnO -one of the most used ENMs- in fresh municipal solid waste (MSW) was investigated. Batch reactors containing municipal solid waste samples were spiked with three different types of nano ZnO having different surface stabilization. The leaching of ZnO was examined under acidic, basic and elevated ionic strength (IS) conditions. The results of the 3-day batch tests showed that the percent of the added nano-ZnO mass retained within the solid waste matrix ranged between 80% and 93% on average for the three types of nano-ZnO tested. The pH and IS conditions did not significantly influence the leaching behavior of ZnO. To further analyze the behavior of ZnO in the MSW matrix, a kinetic particle deposition/detachment model was developed. The model was able to reproduce the main trends of the batch experiments. Reaction rate constants for the batch tests ranged from 0.01 to 0.4 1/hr, reflecting the rapid deposition of nano-ZnO within the MSW matrix.

Effect of temperature on cosolvent flooding for the enhanced solubilization and mobilization of NAPLs in porous media.

This paper examines the potential for enhanced NAPL recovery from the subsurface through the combined application of hot water and cosolvent flushing. Batch experiments were conducted to determine the effect of temperature on fluid properties and the multiphase behavior of the ethanol-water-toluene system and to assess the impact of temperature on the capillary, Bond and total trapping numbers and on flooding stability. Column flooding experiments were also conducted to evaluate toluene NAPL recovery efficiency for different ethanol contents and flushing solution temperatures. The ethanol content considered ranged from 20 to 100% by mass, while the flushing solution temperatures were varied from 10 to 40°C. It is shown that small variations in the system temperature can strongly influence the solubilization, mobilization and stability of the multiphase system, but that the impact of temperature on the enhanced NAPL recovery is also dependent on the ethanol content of the flushing solution. The impact of hot water on NAPL recovery was most pronounced at intermediate ethanol contents (40-60% by mass) where the increase in system temperature led to enhanced NAPL solubilization as well as NAPL mobilization. This study demonstrates that coupling of hot water with in situ cosolvent flooding is a potentially effective remedial alternative that can optimize NAPL recovery while reducing the amount of chemicals injected into the subsurface.

Biosorption of neodymium on Chlorella vulgaris in aqueous solution obtained from hard disk drive magnets

In recent years, biosorption is being considered as an environmental friendly technology for the recovery of rare earth metals (REE). This study investigates the optimal conditions for the biosorption of neodymium (Nd) from an aqueous solution derived from hard drive disk magnets using green microalgae (Chlorella vulgaris). The parameters considered include solution pH, temperature and biosorbent dosage. Best-fit equilibrium as well as kinetic biosorption models were also developed. At the optimal pH of 5, the maximum experimental Nd uptakes at 21, 35 and 50°C and an initial Nd concentration of 250 mg/L were 126.13, 157.40 and 77.10 mg/g, respectively. Analysis of the optimal equilibrium sorption data showed that the data fitted well (R2 = 0.98) to the Langmuir isotherm model, with maximum monolayer coverage capacity (qmax) of 188.68 mg/g, and Langmuir isotherm constant (KL) of 0.029 L/mg. The corresponding separation factor (RL) is 0.12 indicating that the equilibrium sorption was favorable. The sorption kinetics of Nd ion follows well a pseudo-second order model (R2>0.99), even at low initial concentrations. These results show that Chlorella vulgaris has greater biosorption affinity for Nd than activated carbon and other algae types such as: A. Gracilis, Sargassum sp. and A. Densus.