The current state of knowledge of ecosystems and ecosystem services in Russia: A status report

This paper focusses on a conceptual overview of ways to address a comprehensive analysis of ecosystem services (ES) in a country as large and heterogeneous as Russia. As a first step, a methodology for assessing the services for the federal subjects of Russia was chosen, i.e., its constituent provinces and similar entities, in physical terms. Russia harbors a great diversity of natural conditions and ecosystems which are suppliers of ES, and likewise a variety of the socio-economic conditions that shape the demand for these services and their consumption. The methodological approach described permits several important tasks to be addressed: the evaluation of the degree of satisfaction of people’s needs for ES, the identification of ecological donor and acceptor regions, and zoning of the country’s territory for ES assessment. The next step is to prepare a prototype of a National Report on ES in Russia, for which we are presenting the planned structure.     

Development of cleaner-burning brick kilns in Ciudad Juarez, Chihuahua, Mexico

The following results provide a comparison between net airborne contamination produced by the traditional form of kiln used in Northern Mexico and by those modified according to a design by Dr. Robert Marquez. What has become known as the MK style kiln was intended to significantly reduce contaminant emissions. The concept involves covering the kiln with a dome and channeling the output of an active kiln through a second, identical loaded kiln for its additional filtration of the effluents. Kilns of a pair are connected via clay brick channels. The roles are reversed after the initial kiln is refilled. Significant reductions in the particulate and gaseous emissions were achieved in the prototype system, but a connectional problem with recent kiln pairs has also limited the degree of operational success. The problem did not mask the potential of the MK kiln, as will be shown. Additional anticipated benefits to the owners of MK kilns, such as reduced operating cycles and decreased quantities of fuel, also have been verified. Key measurements made during all of the burns were of aerosol densities and buoyancies in the flues, kiln temperatures, and, on a number of occasions, chemical analyses of both aerosol and gaseous effluents. Continuous time histories of aerosol densities for most burns (of a total of -40) provide a basis for examining features and the effects of differing styles of operation with respect to burn efficiency and net contaminant masses. Covering the active kiln with a dome produces a net reduction in dry aerosol effluent mass of a factor between 5 and 10, whereas the addition of a filter kiln produces a net reduction of about a factor of 2. The use of used motor oil as a fuel further reduced aerosol contamination by -1 order of magnitude.     

The Role of Viscosity in Fabric Filtration

Abstract

An approach for measuring point-source emissions of volatile organic compounds (VOCs), acidic vapors, and other species is presented. The amount emitted is determined by directly measuring the actual weight gain of an adsorbent bed over a period of time, which is a cumulative rather than a grabbed sample. As a result, wide fluctuations of concentration and erratic flow behavior during sampling are accommodated with no apparent effect on the accuracy of the measured emission rate. The emission rate is determined by a mass balance including the mass change of the sorbent, as well as the influent and effluent humidities.

Validation tests used a known mass flow rate of vapor in a carrier gas, which was compared with the amount measured. The vapor was a single VOC, a mixture of VOCs, or a mixture of a VOC with water. Conditions studied were the compound or mixture of compounds, concentration, carrier gas, flow rate, and adsorbent. In some tests the VOC was admitted intermittently. The VOCs included n-hexane, acetone, toluene, vinyl acetate, and 1,1,1 trichloroethane. For 105 tests, the average absolute discrepancy of the delivered and measured emission rates was 6.8% and the standard deviation was 3.4%.     

The integrated biomarker response revisited: optimization to avoid misuse

The growing need to evaluate the quality of aquatic ecosystems led to the development of numerous monitoring tools. Among them, the development of biomarker-based procedures, that combine precocity and relevance, is recommended. However, multi-biomarker approaches are often hard to interpret, and produce results that are not easy to integrate in the environmental policies framework. Integrative index have been developed, and one of the most used is the integrated biomarker response (IBR). However, an analysis of available literature demonstrated that the IBR suffers from a frequent misuse and a bias in its calculation. Then, we propose here a new calculation method based on both a more simple formula and a permutation procedure. Together, these improvements should rightly avoid the misuse and bias that were recorded. Additionally, a case study illustrates how the new procedure enabled to perform a reliable classification of site along a pollution gradient based on biomarker responses used in the IBR calculations.     

Effect of diffusion limitation and substrate inhibition on steady states of a biofilm reactor treating a single pollutant

The occurrence of multiple steady states in a toluene biodegrading, diffusion-limited biofilm under aerobic conditions was investigated by computer models: one steady-state, and one nonsteady-state. Two stable and one unstable intermediate steady-state were identified in a narrow set of combinations of parameters values. The nonsteady-state model predicts conditions that evolve to a steady state that is within 0.02–1% of the solution of the steady-state model, depending on the number of grid points used, confirming the algorithms are valid. Multiple steady states occur if, (1) a biofilm is exposed to a constant gas-phase pollution concentration, which exceeds or undershoots a certain threshold, (2) in a narrow range of parameter values and (3) provided that the pollutant degradation follows Haldane kinetics. Such a biofilm displays half-saturation (i.e., Michaelis-Menten)-like apparent (“falsified”) kinetics from a concentration range starting at zero up to the occurrence of a second steady state. Multiple steady states and falsified kinetics can negatively affect a biofilter and the experimental determination of kinetic parameters, respectively.

Implications: The occurrence of multiple steady states in a VOC treating biofilm, shows the significant impact of degradation kinetics and diffusion limitation on the biofilm behavior. Moreover, the implied possible sudden drop of removal efficiency of a biofilter, based on the occurrence of multiple steady states lead to possible bottle-necks in biofilter application and operation.     

Effect of Advanced Oxidants Generated Via Ultraviolet Light on a Sequentially Loaded and Regenerated Granular Activated Carbon Biofilter

Abstract

The objective of this research was to investigate a sequentially loaded and regenerated granular activated carbon (GAC) biofilter system and to determine whether regenerative ozonation/advanced oxidation could improve the removal and biodegradation of a volatile organic compound from a contaminated airstream. Bench-scale reactors were constructed to operate in a manner analogous to a commercially available system manufactured by Terr-Aqua Environmental Systems (only with longer contact time). The GAC system consisted of two GAC biofilter beds that operated in a cyclical manner. On a given day, the first GAC bed adsorbed methyl isobutyl ketone from a simulated waste airstream, while the second bed underwent regeneration; then on the next day, the second bed was in the adsorption mode while the first was regenerated.

Three bench-scale systems were used to compare the performance under three operating conditions: (1) ozone/ associated oxidant regeneration of a GAC biofilter system that was seeded with microorganisms from a field site, (2) a humid air regeneration of a seeded GAC biofilter, and (3) a humid air regeneration of an unseeded GAC biofilter. For the advanced oxidant regenerated GAC biofilter, a maximum removal efficiency of >95% was achieved with an empty bed contact time of 148 sec and an influent concentration of 125 ppm methyl isobutyl ketone, and 90–95% was achieved at 148-sec empty bed contact time and a 1150-ppm influent.     

Effect of advanced oxidants generated via ultraviolet light on a sequentially loaded and regenerated granular activated carbon biofilter

The objective of this research was to investigate a sequentially loaded and regenerated granular activated carbon (GAC) biofilter system and to determine whether regenerative ozonation/advanced oxidation could improve the removal and biodegradation of a volatile organic compound from a contaminated airstream. Bench-scale reactors were constructed to operate in a manner analogous to a commercially available system manufactured by Terr-Aqua Environmental Systems (only with longer contact time). The GAC system consisted of two GAC biofilter beds that operated in a cyclical manner. On a given day, the first GAC bed adsorbed methyl isobutyl ketone from a simulated waste airstream, while the second bed underwent regeneration; then on the next day, the second bed was in the adsorption mode while the first was regenerated. Three bench-scale systems were used to compare the performance under three operating conditions: (1) ozone/ associated oxidant regeneration of a GAC biofilter system that was seeded with microorganisms from a field site, (2) a humid air regeneration of a seeded GAC biofilter, and (3) a humid air regeneration of an unseeded GAC biofilter. For the advanced oxidant regenerated GAC biofilter, a maximum removal efficiency of >95% was achieved with an empty bed contact time of 148 sec and an influent concentration of 125 ppm methyl isobutyl ketone, and 90-95% was achieved at 148-sec empty bed contact time and a 1150-ppm influent.     

Integrated catalytic wet air oxidation and aerobic biological treatment in a municipal WWTP of a high-strength o-cresol wastewater

This study examines the feasibility of coupling a Catalytic Wet Air Oxidation (CWAO), with activated carbon (AC) as catalyst, and an aerobic biological treatment to treat a high-strength o-cresol wastewater. Two goals are pursued: (a) To determine the effect of the main AC/CWAO intermediates on the activated sludge of a municipal WasteWater Treatment Plant (WWTP) and (b) To demonstrate the feasibility of coupling the AC/CWAO effluent as a part of the influent of a municipal WWTP. In a previous study, a high-strength o-cresol wastewater was treated by AC/CWAO aiming to establish the distribution of intermediates and the biodegradability enhancement. In this work, the biodegradability, toxicity and inhibition of the most relevant intermediates detected in the AC/CWAO effluent were determined by respirometry. Also, the results of a pilot scale municipal WWTP study for an integrated AC/CWAO-aerobic biological treatment of this effluent are presented. The biodegradation parameters (i.e. maximum oxygen uptake rate and oxygen consumption) of main AC/CWAO intermediates allowed the classification of the intermediates into readily biodegradable, inert or toxic/inhibitory compounds. This detailed study, allowed to understand the biodegradability enhancement exhibited by an AC/CWAO effluent and to achieve a successful strategy for coupling the AC/CWAO step with an aerobic biological treatment for a high-strength o-cresol wastewater. Using 30%, as COD, of AC/CWAO effluent in the inlet to the pilot scale WWTP, the integrated AC/CWAO-biological treatment achieved a 98% of total COD removal and, particularly, a 91% of AC/CWAO effluent COD removal without any undesirable effect on the biomass.     

Formulation of a soil-pesticide transport model based on a compartmental approach

A semianalytical soil-pesticide transport model is formulated based on a compartmental approach to determine spatial and temporal variations of pesticide residues across a soil profile. The compartmental model is implemented by drawing an analogy between a series of continuous-flow stirred tank reactors and a soil horizon that consists of multiple perfectly mixed compartments. The analogy is strengthened by exploiting a relation between the compartment series and the conventional convective-dispersive equation (CDE) for vertical transport in the soil. Consequently, the number of compartments in the model formulation is not free, but dictated as a function of transport parameters. The model formulation allows consideration of arbitrary boundary value specifications and also, for some cases, spatially varying initial concentration profiles. Sorption kinetics is represented via a two-site model that involves a linear sorption isotherm and a first-order irreversible sorption or a radial diffusive penetrating model. For these three cases, analysis of the compartmental model allows the resultant concentration profiles to be expressed in terms of the Poisson distribution. When a nonlinear kinetic sorption model is used to simulate the sorption processes, an analytical solution is not found and a numerical approach is required.     

Stabilization of APC residues from waste incineration with ferrous sulfate on a semi-industrial scale

A stabilization method for air pollution control (APC) residues from municipal solid waste incineration (MSWI) involving mixing of the residue with water and FeSO4 has been demonstrated on a semi-industrial scale on three types of APC residues: a semidry (SD) APC residue, a fly ash (FA), and an FA mixed with sludge (FAS) from a wet flue gas cleaning system. The process was performed in batches of 165-175 kg residue. It generates a wastewater that is highly saline but has a low content of heavy metals such as Cd, Cr, and Pb. The stabilized and raw residues have been subject to a range of leaching tests: the batch leaching test, the pH-static leaching test, the availability test, and the column test. These tests showed that the stabilized residues have remarkably improved leaching properties, especially with respect to Pb but also with respect to Cd, Cu, and Zn. The release of Pb was reduced by a factor of 250-36,000.     

The effects of coal-smoke pollutants on the growth, yield and leaf epidermal features of Abelmoschus esculentus moench

In Abelmoschus esculentus var. Pusa Sawani, coal-smoke pollutants, experienced close to a thermal power plant, led to decreases of plant height and stem diameter, jeopardised the production of leaves and flower buds and stimulated leaf and flower fall, leading to a decrease in the number of fruits set per plant. The length of fruits was affected a little, and the circumference and fresh weight suffered significantly, causing a decrease in plant yield. The net productivity of stems, roots and leaves incurred losses of 22.5%, 24.0% and 37.5%, respectively, and resulted in a 28.3% loss in the total net productivity of polluted plants. A decrease in stomatal density, aperture and index, and a lower density of epidermal cells of greater dimensions were observed on both the leaf surfaces. Abnormalities of the stomatal complexes were also observed on the upper surface of polluted leaves.     

Natural phillipsite as a matrix for a slow-release formulation of oxamyl

The performance of phillipsite as a matrix for slow-release formulation of oxamyl [N,N-dimethyl-2-methylcarbamoyl-oxymino-2-(methylthio)acetamide] was tested. The adsorption kinetics followed a first-order law, and the adsorption isotherm fitted well in a two-surface Langmuir model, suggesting a double mechanism of interaction between oxamyl and the sorbent. The sorption mechanism, studied by FTIR, provided two fractions of oxamyl. The first one is sorbed on the mineral surface, linked by H-bonding, and the second one is constituted by a multilayer of oxamyl molecules linked by a water bridge between them. The release kinetics of oxamyl from a substratum zeolite-oxamyl also follows a first-order law, with two stages that correspond to both fractions of oxamyl previously detected.     

Boiler briquette coal versus raw coal: Part II--Energy, greenhouse gas, and air quality implications

The objective of this paper is to conduct an integrated analysis of the energy, greenhouse gas, and air quality impacts of a new type of boiler briquette coal (BB-coal) in contrast to those of the raw coal from which the BB-coal was formulated (R-coal). The analysis is based on the source emissions data and other relevant data collected in the present study and employs approaches including the construction of carbon, energy, and sulfur balances. The results show that replacing R-coal with BB-coal as the fuel for boilers such as the one tested would have multiple benefits, including a 37% increase in boiler thermal efficiency, a 25% reduction in fuel demand, a 26% reduction in CO2 emission, a 17% reduction in CO emission, a 63% reduction in SO2 emission, a 97% reduction in fly ash and fly ash carbon emission, a 22% reduction in PM2.5 mass emission, and a 30% reduction in total emission of five toxic hazardous air pollutant (HAP) metals contained in PM10. These benefits can be achieved with no changes in boiler hardware and with a relatively small amount of tradeoffs: a 30% increase in PM10 mass emission and a 9-16% increase in fuel cost.     

Optimization of the intermittent aeration in a full-scale wastewater treatment plant biological reactor for nitrogen removal.

The optimization of the intermittent aeration in a full-scale biological reactor treating municipal and industrial wastewater has been studied by means of an experimental design, where the main factors considered have been the on-off period of aeration and the dissolved-oxygen set point. The objective of the work has been to determine the optimal configuration of the aeration control to achieve a simultaneous removal of carbon (expressed as chemical oxygen demand [COD]) and nitrogen. The results were obtained in a full-scale bioreactor located at the wastewater treatment plant of Sant Celoni (Barcelona, Spain), where a nitrogen removal based on a nitrification-denitrification process is being tested. It has been shown that the crucial parameter to ensure a complete removal of COD and nitrogen is to give the system a long time of on-off aeration to complete both the nitrification and the denitrification processes, whereas the dissolved-oxygen set point has a minor influence on the overall performance. The results presented can be applied to a great number of similar systems, in which a nitrogen removal, based on a nitrification-denitrification scheme, is currently being tested or implemented.     

Rotary Kiln Incineration II. Laboratory-Scale Desorption and Kiln-Simulator Studies—Solids

With landfill costs increasing and regulations on landfilling becoming more stringent, alternatives to conventional hazardous waste treatment strategies are becoming more desirable. Incineration Is presently a permanent, proven solution for the disposal of most organic contaminants, but also a costly one, especially in the case of solids which require some auxiliary fuel. The goal of this research is to develop an understanding of the phenomena associated with the evolution of contaminants from solids In the primary combustor of an Incineration system. A four-fold approach is being used. First, a bench-scale particle characterization reactor was developed to study the transport phenomena on a particle basis, where the controlling processes are mainly intraparticle. Second, a bed-characterization reactor was built to examine the controlling transport phenomena within a bed of particles, where the processes are primarily interparticle. The results of these studies can be applied to any primary combustor. A pilot-scale rotary kiln was developed to study the evolution of contaminants from solids within a realistic temperature and rotation environment. Finally, in situ measurements are being obtained from a full-scale rotary-kiln.

This paper describes results obtained in a study using a commercial sorbent contaminated with toluene. The data are from the particle-characterization reactor and the rotary-kiln simulator. The results show that the method of contamination and charge size do not have a large effect on desorption, while temperature and contaminant concentration are important parameters In the evolution of contaminants in a rotary kiln.     

Calcium carbonate phosphate binding ion exchange filtration and accelerated denitrification improve public health standards and combat eutrophication in aquatic ecosystems.

Cultural eutrophication, the process by which a lake becomes rich in dissolved nutrients as a result of point and nonpoint pollutant sources, is a major cause of the loss of natural lake ecosystems throughout the world. The process occurs naturally in all lakes, but phosphate-rich nutrient runoff from sources such as storm drains and agricultural runoff is a major cause of excess phosphate-induced eutrophication. Especially in Madrona Marsh, one of the last remaining vernal marshes in the greater Los Angeles area, California, cultural eutrophication has become a major problem. In this study, calcium carbonate was found to be an excellent phosphate binder, reducing up to 70% of the phosphates in a given sample of water, and it posed relatively negligent ecological repercussions. This study involved the testing of this principle in both the laboratory and the real ecosystem. A calcium carbonate lacing procedure was first carried out to determine its efficacy in Madrona Marsh. Through this, ammonia was found to interfere with the solubility of calcium carbonate and therefore to be a hindrance to the reduction of phosphate. Therefore, various approaches for reduction of ammonia were tested, including aeration, use of fiber growth media, and plants, mainly Caulerpa verticellata, chosen for it hardiness, primarily in an attempt to increase population of Nitrobacter and Nitrosomonas. All were successful in moderately reducing ammonia levels. In addition, soil sampling, sediment analysis, microscopic plant analysis, microorganism and macroinvertebrate identification, and rate law formulations were conducted. The effect of phosphate and ammonia reduction on the populations of enterobacteria was also an important focus of this experiment. Varying concentrations of phosphate, ammonia, and calcium carbonate in conjunction with phosphate were tested in Madrona Marsh to determine their effects on the populations of enteropathogens on nonspecific blood agar, MacConkey agar, and Hektoen agar. Initial analyses suggest a strong correlation between phosphate concentrations and bacterial populations; a 66% decrease in phosphate resulted in a 35% reduction in bacterial populations and a 45% reduction in enteropathogenic populations. Likewise, a strong correlation was shown between calcium carbonate concentrations and bacterial reduction greater than that which can be attributed to the phosphate reduction alone. This was followed by the construction of various phosphate binding calcium carbonate filters, which used the ion exchange principle, including a spring loading filter, PVC pipe filter, and a galvanized filter. All were tested with the aid of Stoke's law formulation. The experiment was extremely successful in designing a working phosphate-binding and ammonia-reducing filter, and a large-scale agitator-clarifier filter system is currently being planned for construction in Madrona Marsh; this filter will reduce phosphate and ammonia levels substantially in the following years, bringing ecological, economical, and health-related improvements to the overall ecosystem and habitat.     

Sustainable development and the worldwide market economy: a necessary but difficult reconciliation

The only satisfactory concept of development is a global one, of which the usual notion of sustainable development is only a part. The only economic system providing efficiency is the market, with a high level of responsibility given to decentralised and autonomous decision-makers. However, the market economy does not automatically generate globally desirable development. It has to be complemented by public intervention, leading in a single country to a mixed economy, and in a world-set of countries to international agreements, creating two kinds of international order, at two levels: very restrictive between countries with common basic political and social values, and provisionally more flexible and concerning only few very fundamental topics at a world-wide level.     

Effects of pore volume-transmissivity correlation on transport phenomena

The relevant velocity that describes transport phenomena in a porous medium is the pore velocity. For this reason, one needs not only to describe the variability of transmissivity, which fully determines the Darcy velocity field for given source terms and boundary conditions, but also any variability of the pore volume. We demonstrate that hydraulically equivalent media with exactly the same transmissivity field can produce dramatic differences in the displacement of a solute if they have different pore volume distributions. In particular, we demonstrate that correlation between pore volume and transmissivity leads to a much smoother and more homogeneous solute distribution. This was observed in a laboratory experiment performed in artificial fractures made of two plexiglass plates into which a space-dependent aperture distribution was milled. Using visualization by a light transmission technique, we observe that the solute behaviour is much smoother and more regular after the fractures are filled with glass powder, which plays the role of a homogeneous fault gouge material. This is due to a perfect correlation between pore volume and transmissivity that causes pore velocity to be not directly dependent on the transmissivity, but only indirectly through the hydraulic gradient, which is a much smoother function due to the diffusive behaviour of the flow equation acting as a filter. This smoothing property of the pore volume-transmissivity correlation is also supported by numerical simulations of tracer tests in a dipole flow field. Three different conceptual models are used: an empty fracture, a rough-walled fracture filled with a homogeneous material and a parallel-plate fracture with a heterogeneous fault gouge. All three models are hydraulically equivalent, yet they have a different pore volume distribution. Even if piezometric heads and specific flow rates are exactly the same at any point of the domain, the transport process differs dramatically. These differences make it important to discriminate in situ among different conceptual models in order to simulate correctly the transport phenomena. For this reason, we study the solute breakthrough and recovery curves at the extraction wells. Our numerical case studies show that discrimination on the basis of such data might be impossible except under very favourable conditions, i.e. the integral scale of the transmissivity field has to be known and small compared to the dipole size. If the latter conditions are satisfied, discrimination between the rough-walled fracture filled with a homogeneous material and the other two models becomes possible, whereas the parallel-plate fracture with a heterogeneous fault gouge and the empty fracture still show identifiability problems. The latter may be solved by inspection of aperture and pressure testing.     

A multiple testing approach for hazard evaluation of complex mixtures in the aquatic environment: the use of diesel oil as a model

Traditional single species toxicity tests and multiple component laboratory-scaled microcosm assays were combined to assess the toxicological hazard of diesel oil, a model complex mixture, to a model aquatic environment. The immediate impact of diesel oil dosed on a freshwater community was studied in a model pond microcosm over 14 days: a 7-day dosage and a 7-day recovery period. A multicomponent laboratory microcosm was designed to monitor the biological effects of diesel oil (1.0 mg litre(-1)) on four components: water, sediment (soil + microbiota), plants (aquatic macrophytes and algae), and animals (zooplanktonic and zoobenthic invertebrates). To determine the sensitivity of each part of the community to diesel oil contamination and how this model community recovered when the oil dissipated, limnological, toxicological, and microbiological variables were considered. Our model revealed these significant occurrences during the spill period: first, a community production and respiration perturbation, characterized in the water column by a decrease in dissolved oxygen and redox potential and a concomitant increase in alkalinity and conductivity; second, marked changes in microbiota of sediments that included bacterial heterotrophic dominance and a high heterotrophic index (0.6), increased bacterial productivity, and the marked increases in numbers of saprophytic bacteria (10 x) and bacterial oil degraders (1000 x); and third, column water acutely toxic (100% mortality) to two model taxa: Selenastrum capricornutum and Daphnia magna. Following the simulated clean-up procedure to remove the oil slick, the recovery period of this freshwater microcosm was characterized by a return to control values. This experimental design emphasized monitoring toxicological responses in aquatic microcosm; hence, we proposed the term 'toxicosm' to describe this approach to aquatic toxicological hazard evaluation. The toxicosm as a valuable toxicological tool for screening aquatic contaminants was demonstrated using diesel oil as a model complex mixture.     

Predictive modelling of NAPL injection tests in variable aperture spatially correlated fractures.

In preparation for a field experiment where a NAPL will be injected into a fractured sandstone aquifer, a 2D invasion percolation model of DNAPL migration in a single horizontal fracture with varying aperture has been developed. This simulation investigated the effect of spatially correlated fracture aperture fields on pressure-saturation relationships as a function of variable aperture mean, standard deviation, and spatial correlation statistics under hydrostatic conditions. Results from spatially correlated variable aperture fields can be significantly different from random fields. Longer ranges decreased entry pressures and higher standard deviations decreased nonwetting phase saturations. Mean aperture is the major control on displacement pressure, entry pressure and the form of the pressure-saturation curve. Simulation results using statistical parameters for a variable aperture natural sandstone fracture as described by Yeo et al. [International Journal of Rock Mechanics and Mining Sciences 35 (1998) 1051] closely resemble a Brooks-Corey pressure-saturation function, and exhibit a distinct entry pressure followed by a rapid increase in nonwetting phase saturation. Graphical estimates of entry pressure provide a good approximation of the critical aperture controlling the formation of a continuous nonwetting phase pathway in a variable aperture fracture. Consequently, we show that multiphase flow concepts developed for porous media can successfully be applied to variable aperture fractures. Entry pressure correlates well to the mean aperture in these simulations when using a Gaussian distribution of fracture aperture. Interpreted aperture distributions from NAPL injection experiments do not fit the true distribution well at low nonwetting phase saturations, but do at higher saturations above the entry pressure. Consequently, true, mechanical aperture variation within a fracture plane cannot be determined from NAPL injection experiments either in the field or laboratory.