Search

Assessing the impact of a land-use change (LUC) or change in land-use management on nonpoint source-driven groundwater quality in heterogeneous aquifers requires complex analysis. Stochastic methods have been used to account for prediction uncertainty but at high computational cost, which significantly limits the application of these approaches. As an efficient alternative, this study evaluates the application of a meta-analytical solution for evaluating the change in contaminant breakthrough curves at extraction wells in response to LUC. The solution uses the concentration percentiles from a reference stochastic simulation of water flow and solute transport in a groundwater system, assuming a reference land-use distribution pattern. Reference land-use controls the spatially variable rates of both, recharge and contaminant mass loading. The effect of a LUC is evaluated by scaling the ratio between the reference and the new (post-LUC) average input concentrations. The validity of the proposed meta-analysis tool is tested by comparing the results of the meta-analytical solution with those from a full stochastic simulation of the post-LUC scenario. Simulation results show that the accuracy of the meta-analytical solution is best when the regional average recharge rates for both pre- and post-LUC remain approximately unchanged, for any change in contaminant mass loading. Results also indicate that changes in spatial variability and pattern of the recharge rate do not significantly impact the flow field, travel times, and resulting concentrations, if the magnitude of local recharge remains about the same. Lastly, the results show large variability among wells of (and—for an individual well—uncertainty about) the time lag between the time of LUC and the time of consequential effective change in concentrations across wells in the affected region, captured here using statistical metrics.

Non-point source (NPS) pollution has degraded groundwater quality of unconsolidated sedimentary basins over many decades. Properly conceptualizing NPS pollution from the well scale to the regional scale leads to complex and expensive numerical models: key controlling factors of NPS pollution – recharge rate, leakage of pollutants, and soil and aquifer hydraulic properties – are spatially and, for recharge and pollutant leakage, temporally variable. This leads to high uncertainty in predicting well pollution. On the other hand, concentration levels of some key NPS contaminants (salinity, nitrate) vary within a limited range (

Lagrangian methods, such as the random-walk particle-tracking (RWPT), are often qualified as a potentially valuable alternative to error-prone Eulerian methods for simulating solute transport in unsaturated porous media. Yet, the RWPT method has not yet been validated against - and compared to - currently used Eulerian solutions for simulating solute transport under a range of typical unsaturated conditions. This paper presents a new implementation of the RWPT approach for advective - dispersive transport problems under variably saturated conditions. We first show that, as previously demonstrated for a heterogeneous dispersion tensor, using an interpolation scheme in the RWPT algorithm performs well for problems with abrupt changes in the water content. The new model is then compared against a simple 1D uniform transport problem, for which an analytical solution exist, and against a variety of 1D and 3D numerical solutions using the different Eulerian schemes implemented in Hydrus software suite. Results show that, while the Eulerian solutions significantly suffer from numerical dispersion in case of a coarse spatial discretization of the simulation domain, the new Lagrangian model provides accurate solutions for all problems. Furthermore, RWPT reproduces accurately solute transport for typical unsaturated flow conditions (infiltration, evaporation). Moreover, the Lagrangian model appears to be orders of magnitude faster than its Eulerian alternative to solve a 3D heterogeneous problem. Thus, RWPT should be seen as an attractive, stable and efficient alternative for simulating solute transport in the vadose zone, especially in case of complex and large problems.

Author(s): Guo, Z; Fogg, GE; Henri, CV | Abstract: Regional scale transport models are needed to support the long-term evaluation of groundwater quality and to develop management strategies aiming to prevent serious groundwater degradation. The purpose of this study is to evaluate the capacity of a previously developed upscaling approach to adequately describe the main solute transport processes, including the capture of late-time tails under changing boundary conditions. Potential factors that impact the performance of upscaling methods, including temporal variations in mass transfer rates and mass distributions, were investigated. Advective-dispersive contaminant transport in a 3-D heterogeneous domain was simulated and used as a reference solution. The equivalent transport under homogeneous flow conditions was then evaluated by applying the multirate mass transfer (MRMT) model. The random walk particle tracking method was used to solve the solute transport for heterogeneous and homogeneous MRMT scenarios under steady state and transient conditions. The results indicate that the MRMT model can capture the tails satisfactorily for plumes transported with ambient steady state flow fields at all studied scales using the same parameters. However, when the boundary conditions change in either local, plume, or regional scale, the mass transfer model calibrated for transport under steady state conditions cannot accurately reproduce the tailings observed for the heterogeneous scenario. The deteriorating impacts of transient boundary conditions on the upscaled model are more significant for regions where the flow fields are dramatically affected, which highlights the poor applicability of the MRMT approach for complex field settings. This finding also has implications for the suitability of other potential upscaling approaches.

Heterogeneities in natural porous formations, mainly manifested through the hydraulic conductivity (K) and, to a lesser degree, the porosity (ϕ), largely control subsurface flow and solute transport. The influence of the heterogeneous structure of K on transport processes has been widely studied, whereas less attention is dedicated to the joint heterogeneity of conductivity and porosity fields. Our study employs Monte Carlo simulations to investigate the coupled effect of K − ϕ spatial variability on the transport behavior (and uncertainty) of conservative and reactive plumes within a 3D aquifer domain. We explore multiple scenarios, characterized by different levels of heterogeneity of the geological properties, and compare the computational results from the joint K − ϕ heterogeneous system with the results originating from generally adopted constant ϕ conditions. In our study, the spatially variable K − ϕ fields are positively correlated. We statistically analyze key Environmental Performance Metrics: first arrival times and peak mass fluxes for non-reactive species and increased lifetime cancer risk for reactive chlorinated solvents. The conservative transport simulations show that considering coupled K − ϕ fields decreases the plume dispersion, increases both the first arrival times of solutes and the peak mass fluxes at the observation planes. A positive correlation between aquifer connectivity and peak mass fluxes is identified for both homogeneous and heterogeneous ϕ. Our conservative transport results indicate that the relevance of ϕ variability can depend on the metric of interest, the control plane-source distance as well as the level of heterogeneity of the conductivity field. The analysis on reactive transport shows that ϕ variability only slightly affects the mean increased lifetime cancer risk at the control planes but leads to a considerable reduction of the cancer risk uncertainty. We also see that the sensitivity of cancer risk towards ϕ heterogeneity can be influenced by the level of variability of the conductivity field, the source-to-control plane distance, but is not affected by the manner in which the contaminant concentration is computed.

Author(s): Pauloo, RA; Fogg, GE; Guo, Z; Henri, CV | Abstract: Regional-scale groundwater quality degradation from nonpoint source pollution threatens the long-term sustainability of major alluvial aquifer-aquitard systems worldwide. Models are needed to simulate regionally upscaled transport processes in heterogeneous geology that produces non-Fickian (anomalous) transport. In this investigation, we show the impact of nonhorizontal and nonvertical mean flow direction on non-Fickian transport within a 3D, kilometer-scale, typical alluvial aquifer-aquitard system. Notably, vertical anisotropy in K and seasonal pumping and recharge can fundamentally change hydraulic gradients and shift the mean flow direction between mostly horizontal and mostly vertical flow. Detailed 3D flow and transport simulations under varying mean flow directions indicate that alterations to hydraulic gradients which control the mean flow direction can lead to increasingly non-Fickian transport. Under mostly horizontal flow, diffusion and slow advection dominant low-K facies slow mass transfer rates from low-K material, and preferential flow along connected high-K networks causes increased spatial spreading along the mean flow direction. Conversely, predominantly vertical flow caused by spatially distributed pumping and recharge shifts mass transfer processes in low-K material from diffusion and slow advection dominant to advection dominant, resulting in vertically oriented particle trajectories that compactly migrate through high-K and low-K facies alike, and increasingly Fickian transport. Thus, mean flow direction transience driven by vertical anisotropy in K and seasonal pumping and recharge can create oscillating transport patterns, that range from persistently non-Fickian to more Fickian. Results illustrate the regional-scale hydrogeologic factors that explain why upscaled transport models fail to capture non-Fickian effects resulting from mean flow direction transience.

Author(s): Guo, Zhilin; Pauloo, Rich; Fogg, Graham E; Henri, Christopher; Zheng, Chunmiao | Abstract: Regional scale transport models are needed to support the long-term evaluation of groundwater quality and to develop management strategies aiming to prevent serious groundwater degradation. The transport dominant process, advection or diffusion, was identified for flow fields with different primary flow directions. The capacities of Multi-Rate Mass Transfer (MRMT) and adaptive Multi-rate Mass Transfer (aMMT), modified from MRMT by updating mass transfer rates with changing velocities, to adequately describe the main solute transport processes, including the capture of late-time tails under changing boundary conditions were evaluated. Advective-dispersive contaminant transport simulated in a 3D heterogeneous medium was used as a reference solution. Equivalent transport under homogeneous flow conditions was then evaluated by applying the MRMT or aMMT models for upscaling. Results indicated that for advection-dominated transport, both the MRMT and aMMT methods can upscale the anomalous transport dynamics affected by sub-grid heterogeneity under transient flow conditions. Whereas, for diffusion-dominated systems, the MRMT model failed to capture the tails of tracer breakthrough curves (BTCs) after the boundary condition changed, but the results from the aMMT model were significantly improved. However, if the overall flow direction changed, both MRMT and aMMT failed to represent the BTC tail generated by the heterogeneous system. In this study, an indicator that describe the primary flow direction in anisotropic heterogeneous domain was developed, and the relationship between the flow direction and the dominant transport process was investigated. The ranges of the indicator, within which the advection or diffusion is dominant, are determined. Therefore, this study not only show the capability of upscaling methods on describing the transport that dominated by different processes, but provides a guide on choosing upscaling methods in field site, which supports long-term management of groundwater.

Author(s): Guo, Z; Henri, CV; Fogg, GE; Zhang, Y; Zheng, C | Abstract: The long-term evaluation of regional-scale groundwater quality needs efficient upscaling methods for transient flow. Upscaling techniques, such as the Multirate Mass Transfer (MRMT) method with constant upscaling parameters, have been used for transport with steady-state flow, yet the upscaling parameters (i.e., rate coefficients) may be time dependent. This study proposed and validated an adaptive MRMT (aMMT) method by allowing the mass transfer coefficients in MRMT to change with the flow field. Advective-dispersive contaminant transport simulated in a 3-D heterogeneous medium was used as a reference solution. Equivalent transport under homogeneous flow conditions was evaluated by applying the MRMT and aMMT models for upscaling. The relationship between mass transfer coefficients and flow rates was fitted under steady-state flow driven by various hydraulic gradients. A power law relationship was obtained, which was then used to update the mass transfer coefficients in each stress period under transient flow conditions in the aMMT method. Results indicated that for advection-dominated transport, both the MRMT and aMMT methods can upscale the anomalous transport dynamics affected by subgrid heterogeneity under transient flow conditions. Whereas for diffusion-dominated systems, the MRMT model failed to capture the tails of tracer breakthrough curves after the boundary condition changed, but the results from the aMMT model were significantly improved. However, if the overall flow direction changed, both MRMT and aMMT failed to represent the breakthrough curve tail generated by the heterogeneous system. The results point toward a promising path for upscaling transport in complex aquifers with transient flow.

Upscaled transport models are needed to address regional-scale groundwater quality degradation in major aquifers worldwide, and these models must accurately characterize anomalous (non-Fickian) tra...

Summary Facies delineation is defined as the separation of geological units with distinct intrinsic characteristics (grain size, hydraulic conductivity, mineralogical composition). A major challenge in this area stems from the fact that only a few scattered pieces of hydrogeological information are available to delineate geological facies. Several methods to delineate facies are available in the literature, ranging from those based only on existing hard data, to those including secondary data or external knowledge about sedimentological patterns. This paper describes a methodology to use kernel regression methods as an effective tool for facies delineation. The method uses both the spatial and the actual sampled values to produce, for each individual hard data point, a locally adaptive steering kernel function, self-adjusting the principal directions of the local anisotropic kernels to the direction of highest local spatial correlation. The method is shown to outperform the nearest neighbor classification method in a number of synthetic aquifers whenever the available number of hard data is small and randomly distributed in space. In the case of exhaustive sampling, the steering kernel regression method converges to the true solution. Simulations ran in a suite of synthetic examples are used to explore the selection of kernel parameters in typical field settings. It is shown that, in practice, a rule of thumb can be used to obtain suboptimal results. The performance of the method is demonstrated to significantly improve when external information regarding facies proportions is incorporated. Remarkably, the method allows for a reasonable reconstruction of the facies connectivity patterns, shown in terms of breakthrough curves performance.

An edited version of this paper was published by AGU. Copyright (2015) American Geophysical Union.

This is the peer reviewed version of the following article: [Ding, D., Benson, D. A., Fernàndez‐Garcia, D., Henri, C. V., Hyndman, D. W., Phanikumar, M. S., & Bolster, D. (2017). Elimination of the reaction rate “scale effect”: Application of the Lagrangian reactive particle‐tracking method to simulate mixing‐limited, field‐scale biodegradation at the Schoolcraft (MI, USA) site. Water Resources Research, 53, 10,411–10,432. https://doi.org/10.1002/2017WR021103], which has been published in final form at https://doi.org/10.1002/2017WR021103. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Self-Archiving. Measured (or empirically fitted) reaction rates at groundwater remediation sites are typically much lower than those found in the same material at the batch or laboratory scale. The reduced rates are commonly attributed to poorer mixing at the larger scales. A variety of methods have been proposed to account for this scaling effect in reactive transport. In this study, we use the Lagrangian particle-tracking and reaction (PTR) method to simulate a field bioremediation experiment at the Schoolcraft, MI site. A denitrifying bacterium, Pseudomonas Stutzeri strain KC (KC), was injected to the aquifer, along with sufficient substrate, to degrade the contaminant, carbon tetrachloride (CT), under anaerobic conditions. The PTR method simulates chemical reactions through probabilistic rules of particle collisions, interactions, and transformations to address the scale effect (lower apparent reaction rates for each level of upscaling, from batch to column to field scale). In contrast to a prior Eulerian reaction model, the PTR method is able to match the field-scale experiment using the rate coefficients obtained from batch experiments.

When laboratory-measured chemical reaction rates are used in simulations at the field-scale, the models typically overpredict the apparent reaction rates. The discrepancy is primarily due to poorer mixing of chemically distinct waters at the larger scale. As a result, realistic field-scale predictions require accurate simulation of the degree of mixing between fluids. The Lagrangian particle-tracking (PT) method is a now-standard way to simulate the transport of conservative or sorbing solutes. The method’s main advantage is the absence of numerical dispersion (and its artificial mixing) when simulating advection. New algorithms allow particles of different species to interact in nonlinear (e.g., bimolecular) reactions. Therefore, the PT methods hold a promise of more accurate field-scale simulation of reactive transport because they eliminate the masking effects of spurious mixing due to advection errors inherent in grid-based methods. A hypothetical field-scale reaction scenario is constructed and run in PT and Eulerian (finite-volume/finite-difference) simulators. Grid-based advection schemes considered here include 1st- to 3rd-order spatially accurate total-variation-diminishing flux-limiting schemes, both of which are widely used in current transport/reaction codes. A homogeneous velocity field in which the Courant number is everywhere unity, so that the chosen Eulerian methods incur no error when simulating advection, shows that both the Eulerian and PT methods can achieve convergence in the L1 (integrated concentration) norm, but neither shows stricter pointwise convergence. In this specific case with a constant dispersion coefficient and bimolecular reaction A + B → P , the correct total amount of product is 0.221MA0, where MA0 is the original mass of reactant A. When the Courant number drops, the grid-based simulations can show remarkable errors due to spurious over- and under-mixing. In a heterogeneous velocity field (keeping the same constant and isotropic dispersion), the PT simulations show an increased reaction total from 0.221MA0 to 0.372MA0 due to fluid deformation, while the 1st-order Eulerian simulations using ≈ 106 cells (with a classical grid Peclet number Δx/αL of 10) have total product of 0.53MA0, or approximately twice as much additional reaction due to advection error. The 3rd-order TVD algorithm fares better, with total product of 0.394MA0, or about 1.14 times the increased reaction total. A very strict requirement on grid Peclet numbers for Eulerian simulations will be required for realistic reactions because of their nonlinear nature. We analytically estimate the magnitude of the effect for the end-member cases of very fast and very slow reactions and show that in either case, the mass produced is proportional to 1 / P e , where Pe is the Peclet number. Therefore, extra mass is produced according to D , where the dispersion includes any numerical dispersion error. We test two PT methods, one that kills particles upon reaction and another that decrements a particle’s mass. For the bimolecular reaction studied here, the computational demands of the particle-killing methods are much smaller than, and the particle-number-preserving algorithm are on par with, the fastest Eulerian methods.

[1] The acquisition of reliable data sets representative of hydrological regimes and their variations is a critical concern for water resource assessment. For the subsurface, traditional approaches based on probe measurements, core analysis, and well data can be laborious, expensive, and highly intrusive, while only yielding sparse data sets. For this study, an innovative field survey, merging relative microgravimetry, magnetic resonance soundings, and hydrological measurements, was conducted to evaluate both surface and subsurface water storage variations in a semiarid Sahelian area. The instrumental setup was implemented in the lower part of a typical hillslope feeding to a temporary pond. Weekly measurements were carried out using relative spring gravimeters during 3 months of the rainy season in 2009 over a 350 × 500 m2 network of 12 microgravity stations. Gravity variations of small to medium amplitude (≤220 nm s−2) were measured with accuracies better than 50 nm s−2, revealing significant variations of the water storage at small time (from 1 week up to 3 months) and space (from a couple of meters up to a few hundred meters) scales. Consistent spatial organization of the water storage variations were detected, suggesting high infiltration at the outlet of a small gully. The comparison with hydrological measurements and magnetic resonance soundings involved that most of the microgravity variations came from the heterogeneity in the vadose zone. The results highlight the potential of time lapse microgravity surveys for detecting intraseasonal water storage variations and providing rich space-time data sets for process investigation or hydrological model calibration/evaluation.

The interplay between the spatial variability of the aquifer hydraulic properties, mass transfer due to sub-grid heterogeneity and chemical reactions often complicates reactive transport simulations. It is well documented that hydro-biochemical properties are ubiquitously heterogeneous and that diffusion and slow advection at the sub-grid scale typically leads to the conceptualization of an aquifer as a multi-porosity system. Within this context, chemical reactions taking place in mobile/immobile water regions can be substantially different between each other. This paper presents a particle-based method that can efficiently simulate heterogeneity, network reactions and multi-rate mass transfer. The approach is based on the development of transition probabilities that describe the likelihood that particles belonging to a given species and mobile/immobile domain at a given time will be transformed into another species and mobile/immobile domain afterwards. The joint effect of mass transfer and sequential degradation is shown to be non-trivial. A characteristic rebound of degradation products can be observed. This late rebound of concentrations is not driven by any change in the flow regime (e.g., pumping ceases in the pump-and-treat remediation strategy) but due to the natural interplay between mass transfer and chemical reactions. To illustrate that the method can simultaneously represent mass transfer, spatially varying properties and network reactions without numerical problems, we have simulated the degradation of tetrachloroethylene (PCE) in a three-dimensional fully heterogeneous aquifer subjected to rate-limited mass transfer. Two types of degradation modes were considered to compare the effect of an active biofilm with that of clay pods present in the aquifer. Results of the two scenarios display significantly differences. Biofilms that promote the degradation of compounds in an immobile region are shown to significantly enhance degradation, rapidly producing daughter products and less tailing.

Modeling multispecies reactive transport in natural systems with strong heterogeneities and complex biochemical reactions is a major challenge for assessing groundwater polluted sites with organic and inorganic contaminants. A large variety of these contaminants react according to serial-parallel reaction networks commonly simplified by a combination of first-order kinetic reactions. In this context, a random-walk particle tracking method is presented. This method is capable of efficiently simulating the motion of particles affected by first-order network reactions in three-dimensional systems, which are represented by spatially variable physical and biochemical coefficients described at high resolution. The approach is based on the development of transition probabilities that describe the likelihood that particles belonging to a given species and location at a given time will be transformed into and moved to another species and location afterward. These probabilities are derived from the solution matrix of the spatial moments governing equations. The method is fully coupled with reactions, free of numerical dispersion and overcomes the inherent numerical problems stemming from the incorporation of heterogeneities to reactive transport codes. In doing this, we demonstrate that the motion of particles follows a standard random walk with time-dependent effective retardation and dispersion parameters that depend on the initial and final chemical state of the particle. The behavior of effective parameters develops as a result of differential retardation effects among species. Moreover, explicit analytic solutions of the transition probability matrix and related particle motions are provided for serial reactions. An example of the effect of heterogeneity on the dechlorination of organic solvents in a threedimensional random porous media shows that the power-law behavior typically observed in conservative tracers breakthrough curves can be largely compromised by the effect of biochemical reactions.

We have studied the reactionspp→ppπ+π-,K+p→K+pπ+π−π, π+p→ π+,pπ+π− and π−p →π+π− at 147 GeV/c using the 30-inch Fermilab hybrid system. All four reactions were detected with the same apparatus and analyzed in the same way. The energy dependence of the channel cross section was found to beAp−0.6+B for thepp reaction andAp−1+B for the other three. About 90% of the cross section at 147 GeV/c can be accounted for by either beam or target diffraction. Some of the remaining cross section may come from double Pomeron exchange reactions which we tried to isolate. We have tested the hypothesis of a factorizable Pomeron and our data indicates a violation of this hypothesis. We show that the 3π mass enhancement in the mass region 1.2–1.4 GeV is diffractively produced in the π± beam reactions. Fourprong, four-constraint and six-prong, four-constraint cross sections are reported.

We have studied the inclusive production of K/sup asteriskplus-or-minus/(890) and Y/sup asteriskplus-or-minus/(1385) in pp, ..pi../sup +/p, and K/sup +/p interactions at 147 GeV/c. The experiment used the Fermilab 30-inch hydrogen bubble chamber with the hybrid spectrometer system. Results are based on a sample of 1916 observed K/sub S/ and 932 observed ..lambda... Inclusive cross sections are given for K/sup asteriskplus-or-minus/ and Y/sup asteriskplus-or-minus/ production from the three beams, and comparisons are made with experiments at other energies. Feynman-x and transverse-momentum-squared distributions are also calculated. The results suggest that the K/sup asterisk-/ is entirely produced in the central region, while the K/sup asterisk+/ includes a component from beam fragmentation. Comparisons are made with the additive quark model.

The gross features of the data are reviewed and the momentum tensor components, 〈S〉, and 〈T〉 are studied. Evidence for coplanarity effects is observed. A probabilistic method to search for jets is proposed. The jets found by this method are seen to be similar to highp T jets frompp interactionse + e − annihilations andvp interactions. The effect of subtracting the leading particles and the effect of lowp T diffractive and non diffractive components are investigated.

The Fermilab hybrid 30-in. bubble-chamber spectrometer was exposed to a tagged 147-GeV/c positive beam containing ..pi../sup +/, K/sup +/, and p. A sample of 3003 K/sup +/p, 19 410 pp, and 20 745 ..pi../sup +/p interactions is used to derive sigma/sub n/, , f/sub 2//sup tsc/c, and /D for each beam particle. These values are compared to values obtained at other, mostly lower, beam momenta. The overall dependence of on E/sub a/, the available center-of-mass energy, for these three reactions as well as ..pi../sup -/p and pp interactions has been determined.

In an experiment with the 30-inch Hybrid Spectrometer at Fermilab we have obtained the inclusive and semi-inclusive production cross sections of the ϱ0 meson using a conventional background subtraction technique. Production cross sections for the ϱ0 are derived as a function of the Feynman scaling variablex, and the transverse variablespt2 andEt=(pt2+M2)1/2. The longitudinal distributions are compared with the (1−x) dependence of the proton and meson valence quark structure functions, using various forms of recombination and fragmentation models. The transverse distributions are compared with thermodynamic models. We give density matrix elements for the ϱ0 production from pions in the extreme forward region.

The ratio of hadronic π+/π− production has been studied in the target fragmentation region of π+ p, π− p,K + p andp p collisions at 147 GeV/c beam momentum. The difference between this ratio in the π+ p and π− p data is less than at lower energies but still present. No difference exists anymore for the π+/π− ratio at this energy betweenK + p andpp collisions. Exclusion of diffraction dissociation and resonance production does not change the main conclusions. The energy dependence of the ratio for π+ p and π− p collisions gives evidence for a universal asymptotic ratio for all beams. Thep t -dependence at low-p t is in agreement with the expectation from geometrical correlation of the valence quarks.

Results are presented from a study of inclusive neutral strange particle production by a 147 GeV/c tagged π+/K+/p beam in the Fermilab 30-inch hydrogen bubble chamber. The experiment made use of the proportional hybrid spectrometer system. Results are based on 995 KS0, 485 Λ, and 83 Λ found in a sample of 132 000 pictures. Cross sections are given for inclusive production of these particles by each of the three beam particles, and comparisons are made with measurements at other energies. Topological cross sections are also calculated, and KNO multiplicity scaling is investigated. Distributions are presented of invariant cross sections as functions of the Feynman scaling variable x and c.m. rapidity y. The transverse momentum-squared distributions with their fitted slopes are also given. Comparisons are made of the production characteristics for the three beam types.

The energy dependence of the average of the charged multiplicity and its dispersion in π + /K + /p interaction on protons at 147 GeV/ c is found to be the same as in e + e − annihilations if an “effective energy” variable is used instead of the total energy. The effective energy S eff is defined as the invariant mass of all secondaries left after the two leading particles have been removed. Fitting the expression aS eff b to the average charge multiplicity 〈 n ch 〉, we find the power b to be in good agreement with the value of 0.25 predicted by Fermi's statistical model and by Landau's hydrodynamical model.

An analysis is presented of the production of jets in low-ptπ±p,k+p,pp collisions at 147 GeV/c collected at Fermilab in the rapid cycling hydrogen bubble chamber hybrid spectrometer. It is shown that selected events are more planar than predicted by a cylindrical phase space, but that 3-jet and 4-jet like structures seen in the momentum flows are reasonably reproduced by the cylindrical phase space events subjected to the same mathematical manipulation required to build momentum flows.

We have searched for the production of long-lived R-..delta../sup + +/ hadrons, which contain a gluino in addition to three up quarks, in pp, ..pi../sup +/p, and K/sup +/p interactions at 147 GeV/c. For R-..delta../sup + +/'s of lifetime greater than 10/sup -9/ sec, and laboratory momentum greater than 2 GeV/c, we place 90%-confidence-level upper limits of 6.1 ..mu..b in pp collisions, 4.4 ..mu..b in ..pi../sup +/p, and 29 ..mu..b in K/sup +/p.

The measurement of 12D 0 and 10D ± exclusive branching ratios are presented. The analysis is based on 608 spatially resolved charm particle decays produced in 360 GeV/c pgr– p and 400 GeV/cpp interactions.

A comparative analysis of the global properties of 400 GeV/c proton-proton interactions with and without charm production is presented. Multiplicities, momentum distributions, sphericity and thrust of these interactions are compared.

Charm data from 360 GeV/c π−p and 400 GeV/cpp interactions are used to give results onD-meson branching ratios. The analysis is based on 438 charm events containing 608 spatially resolved charm particle decays. We present topological branching ratios, as well as the following inclusive branching ratios ofD-mesons into kaons: $$\begin{gathered} B(D^ \pm \to K^ \mp + anything) = 0.17 \pm 0.07, \hfill \\ B(D^ \pm \to K^ \pm + anything) = 0.08_{ - 0.05}^{ + 0.06} , \hfill \\ B(D^0 \to K^ \pm + anything) = 0.42 \pm 0.08, \hfill \\ B(D^0 \to K^ + + anything) = 0.03_{ - 0.02}^{ + 0.05} , \hfill \\ \end{gathered} $$ and the following semi-electronic branching ratios ofD-mesons: $$\begin{gathered} B(D^ \pm \to e^ \pm + anything) = 0.20_{ - 0.07}^{ + 0.09} , \hfill \\ B(D^0 \to e^ \pm + anything) = 0.15 \pm 0.05. \hfill \\ \end{gathered} $$

We report on the Λc properties observed with the LEBC-EHS set-up at CERN SPS with 360 GeV c π − and 400 GeV c p beams. The lifetime is (1.2−0.3+0.5)×10−13 s and the mass (2284.7 ± 2.3 ± 0.5)MeV c 2 .

The release of industrial contaminants into the subsurface has led to a rapid degradation of groundwater resources. Contamination caused by Dense Non-Aqueous Phase Liquids (DNAPLs) is particularly severe owing to their limited solubility, slow dissolution and in many cases high toxicity. A greater insight into how the DNAPL source zone behavior and the contaminant release towards the aquifer impact human health risk is crucial for an appropriate risk management. Risk analysis is further complicated by the uncertainty in aquifer properties and contaminant conditions. This study focuses on the impact of the DNAPL release mode on the human health risk propagation along the aquifer under uncertain conditions. Contaminant concentrations released from the source zone are described using a screening approach with a set of parameters representing several scenarios of DNAPL architecture. The uncertainty in the hydraulic properties is systematically accounted for by high-resolution Monte Carlo simulations. We simulate the release and the transport of the chlorinated solvent perchloroethylene and its carcinogenic degradation products in randomly heterogeneous porous media. The human health risk posed by the chemical mixture of these contaminants is characterized by the low-order statistics and the probability density function of common risk metrics. We show that the zone of high risk (hot spot) is independent of the DNAPL mass release mode, and that the risk amplitude is mostly controlled by heterogeneities and by the source zone architecture. The risk is lower and less uncertain when the source zone is formed mostly by ganglia than by pools. We also illustrate how the source zone efficiency (intensity of the water flux crossing the source zone) affects the risk posed by an exposure to the chemical mixture. Results display that high source zone efficiencies are counter-intuitively beneficial, decreasing the risk because of a reduction in the time available for the production of the highly toxic subspecies.

A study of the properties of charm particles produced in 400 GeV/cpp interactions is reported. The experiment was performed using the high resolution hydrogen bubble chamber LEBC in association with the European Hybrid Spectrometer at the CERN SPS. Details of the experimental set-up and operational procedures are given and the methods to extract samples of charm decays are discussed. Results are presented on the intrinsic properties of charm particles (masses, lifetimes, decay modes and branching ratios), adding, whenever appropriate, the relevant information obtained in a similar study made with 360 GeV/cpgr – p interactions. The hadroproduction properties of charm states (total and differential cross sections, correlations) are presented and discussed in the context of current QCD inspired phenomenology.