Mass Flux Rate Research Articles

Theoretical and numerical analysis of a spreading opposed-flow diffusion flame

This article describes the macroscopic and microscopic features of flames spreading over solid-fuel surfaces by examining and comparing three models. The first model examines ignition and flame spread over a solid-fuel surface using a two-dimensional numerical simulation code. This model employs variable density, variable thermophysical properties and one-step global finite-rate chemistry. The second model, a macroscopic ‘field’ model, is solved in terms of the mixture fraction ( Z ) and total enthalpy ( H ) functions. Comparisons are made with numerical predictions for primitive quantities: temperature, species distributions and velocity fields; and derived quantities: heat flux, mass flux, mixture fraction, enthalpy function and flame stretch rate. The third model yields a ‘localized’ flame structure description near the flame attachment point. Theoretical formulas are produced for the quenching distance, the leading edge heat flux, and the flame structure, as characterized by reactivity, temperature field and species distributions. The analytical predictions are compared with numerical simulations to derive flame microstructure scaling parameters.

ATOMIC JETS FROM CLASS 0 SOURCES DETECTED BYSPITZER: THE CASE OF L1448-C

We present Spitzer-IRS spectra obtained along the molecular jet from the Class 0 source L1448-C (or L1448-mm). Atomic lines from the fundamental transitions of [FeII], [SiII] and [SI] have been detected showing, for the first time, the presence of an embedded atomic jet at low excitation. Pure rotational H$_2$ lines are also detected, and a decrease of the atomic/molecular emission ratio is observed within 1 arcmin from the driving source. Additional ground based spectra (UKIRT/UIST) were obtained to further constrain the H$_2$ excitation along the jet axis and, combined with the 0--0 lines, have been compared with bow-shock models. From the different line ratios, we find that the atomic gas is characterized by an electron density n_e ~ 200-1000 cm^{-3}, a temperature T_e < 2500 K and an ionization fraction <~ 10^{-2}; the excitation conditions of the atomic jet are thus very different from those found in more evolved Class I and Class II jets. We also infer that only a fraction (0.05-0.2) of Fe and Si is in gaseous form, indicating that dust still plays a major role in the depletion of refractory elements. A comparison with the SiO abundance recently derived in the jet from an analysis of several SiO sub-mm transitions, shows that the Si/SiO abundance ratio is ~100, and thus that most of the silicon released from grains by sputtering and grain-grain collisions remains in atomic form. Finally, estimates of the atomic and molecular mass flux rates have been derived: values of the order of ~10$^{-6}$ and ~10$^{-7}$ M$_{\sun}$ yr$^{-1}$ are inferred from the [SI]25$\mu$m and H$_2$ line luminosities, respectively. A comparison with the momentum flux of the CO molecular outflow suggests that the detected atomic jet has the power to drive the large scale outflow.

Evaluating deep updraft formulation in NCAR CAM3 with high‐resolution WRF simulations during ARM TWP‐ICE

The updraft formulation used in NCAR CAM3 deep convection parameterization assumes that the mass flux for a single updraft increases exponentially with height to its top and detrainment is confined only to a thin layer at the updraft top. These assumptions are evaluated against three‐dimensional high‐resolution simulations from the Weather Research and Forecast (WRF) model during the monsoon period of the DOE Atmospheric Radiation Measurement (ARM) Program Tropical Warm Pool–International Cloud Experiment (TWP‐ICE). Analyses of the WRF‐generated updrafts suggest that the mass flux for a single updraft increases with height below the top of the conditionally unstable layer and decreases above. Detrainment may dominate above the conditionally unstable layer rather than only over a thin layer at the updraft top. It is argued that the assumed updraft mass flux profile in CAM3 might be unrealistic in many cases because the updraft acceleration is affected by other drag processes in addition to entrainment. Our analyses suggest that the CAM3‐parameterized convection could be too active and, as a result, excess moisture and heat could be transported to the upper troposphere by the parameterized convection. Future improvement is envisioned.

O I] sub-arcsecond study of a microjet from an intermediate mass young star: RY Tauri

Context. High-resolution studies of microjets in T Tauri stars (cTTs) reveal key information on the jet collimation and launching mechanism, but only a handful of systems have been mapped so far. Aims. We wish to perform a detailed study of the microjet from the 2 Myoung star RY Tau, to investigate the influence of its higher stellar mass and claimed close binarity on jet properties. Methods. Spectro-imaging observations of RY Tau were obtained in (O i)λ6300 with resolutions of 0. 4 and 135 km s −1 ,u sing the integral field spectrograph OASIS at the Canada-France-Hawaii Telescope. Deconvolved images reach a resolution of 0. 2. Results. The blueshifted jet is detected within 2 �� of the central star. We determine its PA, collimation, 2D kinematics, mass-flux rate, ejection to accretion ratio, and transverse velocity shifts taking accurately into account errors due to finite signal to noise ratio. The RY Tau system is shown to provide important constraints to several models of steady MHD ejection. Conclusions. The remarkably similar properties of the RY Tau microjet compared to jets from lower mass cTTs gives support to the common belief that the jet launching mechanism is universal over a broad range of stellar masses. The proximity between the jet PA and the PA of the photocenter variations observed by Hipparcos calls into question the interpretation of the latter in terms of binarity of RY Tau. Partial occultation events of the photosphere may offer an alternative explanation.

DENSITY OSCILLATION FROM NANOSCALE TO MACROSCALE

The effect of channel width on the density structure of confined fluid in the nano-/micro-channels is examined by equilibrium molecular dynamics (EMD) simulation. It was found that the density oscillation occurs near the wall in both cases of the macroscale or nanoscale confined flow. There exists a threshold channel width L threshold , when channel width H&lt;L threshold , density oscillates throughout the channel. When H&gt;L threshold , L threshold is constant and about 5–6 molecular diameter long, and the density becomes uniform beyond this threshold layer. A newly defined ch number may serve to be the parameter to compare similarity in the micro-/nano-scale channel flow. Moreover, the effect of the density oscillation on fluid mass flux rate is examined quantitatively. The result shows that the effect should be considered when the channel width is below 5 molecular diameter.

Mass estimations of ejecta from Strombolian explosions by inversion of Doppler radar measurements

We present a new method for estimating particle loading parameters (mass, number, volume) of eruptive jets by inversion of echo power data measured using a volcano Doppler radar (VOLDORAD) during typical Strombolian activity from the southeast (SE) crater of Mount Etna on 4 July 2001. Derived parameters such as mass flux, particle kinetic and thermal energy, and particle concentration are also estimated. The inversion algorithm uses the complete Mie (1908) formulation of electromagnetic scattering by spherical particles to generate synthetic backscattered power values. In a first data inversion model (termed the polydisperse model), the particle size distribution (PSD) is characterized by a scaled Weibull function. The mode of the distribution is inferred from particle terminal velocities measured by Doppler radar for each explosion. The distribution shape factor is found to be 2.3 from Chouet et al.'s (1974) data for typical Strombolian activity, corresponding to the lognormal PSDs commonly characteristic of other Strombolian deposits. The polydisperse model inversion converges toward the Weibull scale factor producing the best fit between synthetic and measured backscattered power. A cruder, alternative monodisperse model is evaluated on the basis of a single size distribution assumption, the accuracy of which lies within 25% of that of the polydisperse model. Although less accurate, the monodisperse model, being much faster, may be useful for rapid estimation of physical parameters during real‐time volcano monitoring. Results are illustrated for two explosions at Mount Etna with contrasted particle loads. Estimates from the polydisperse model give 58,000 and 206,000 kg as maxima for the total mass of pyroclasts, 26,400 and 73,600 kg s−1 for mass flux rates, 38 and 135 m3 (22 and 76 m3 equivalent magma volume) for the pyroclast volumes, and 0.02–0.4 and 0.06–0.12 kg m−3 for particle concentrations, respectively. The time‐averaged kinetic energy released is found to be equal to 4.2 × 107 and 3.9 × 108 J, and thermal energy is estimated at 8.4 × 1010 and 3 × 1011 J.

Molecular jets driven by high-mass protostars: a detailed study of the IRAS 20126+4104 jet

We present here an extensive analysis of the protostellar jet driven by IRAS 20126+4104, deriving the kinematical, dynamical, and physical conditions of the H2 gas along the flow. The jet has been investigated by means of near-IR H2 and [FeII] narrow-band imaging, high resolution spectroscopy of the 1-0S(1) line (2.12 um), NIR (0.9-2.5 um) low resolution spectroscopy, along with ISO-SWS and LWS spectra (from 2.4 to 200 um). The flow shows a complex morphology. In addition to the large-scale jet precession presented in previous studies, we detect a small-scale wiggling close to the source, that may indicate the presence of a multiple system. The peak radial velocities of the H2 knots range from -42 to -14 km s^-1 in the blue lobe, and from -8 to 47 km s^-1 in the red lobe. The low resolution spectra are rich in H_2 emission, and relatively faint [FeII] (NIR), [OI] and [CII] (FIR) emission is observed in the region close to the source. A warm H2 gas component has an average excitation temperature that ranges between 2000 K and 2500 K. Additionally, the ISO-SWS spectrum reveals the presence of a cold component (520 K), that strongly contributes to the radiative cooling of the flow and plays a major role in the dynamics of the flow. The estimated L(H2) of the jet is 8.2+/-0.7 L_sun, suggesting that IRAS20126+4104 has an accretion rate significantly increased compared to low-mass YSOs. This is also supported by the derived mass flux rate from the H2 lines (Mflux(H2)~7.5x10^-4 M_sun yr^-1). The comparison between the H2 and the outflow parameters strongly indicates that the jet is driving, at least partially, the outflow. As already found for low-mass protostellar jets, the measured H2 outflow luminosity is tightly related to the source bolometric luminosity.

Evaluation and adaptation of a regional climate model for the Horn of Africa: rainfall climatology and interannual variability

AbstractThis study evaluates the ability of the Abdus Salam International Center for Theoretical Physics (ICTP) version 3 Regional Climate Model (RegCM3) to reproduce the observed rainfall amounts and distribution over the topographically varied region of the Horn of Africa. Simulations are performed for the widespread very dry 1984 and locally very wet 1996 years using all of the alternative convective schemes available in the ICTP RegCM3. Extensive comparisons of simulations for 1984 and 1996 reveal that the Emanuel scheme best captures the rainfall patterns and interannual variations for the two extreme years over the Horn of Africa. This scheme therefore is selected for specific adaptation for the region, despite its excessive rainfall estimation there and its demonstrated deficiencies over other regions of Africa not of interest in this study.To reduce rainfall overestimation by the Emanuel scheme, several sensitivity experiments are performed by varying the key parameters that control the rate of convective mass flux, the fraction of condensed water converted to precipitation, and the heating and moistening characteristics of the environment. Results show that the amount of condensed water that ultimately falls out as rain crucially affects simulated rainfall amounts. With an appropriate adjustment to this key parameter, the excessive rainfall amount produced by the Emanuel scheme is reduced substantially. Furthermore, evaluation of RegCM3 simulations for 1982–1999 shows that the modified Emanuel convective scheme not only reproduces the 18‐year average rainfall realistically but also captures the interannual variability adequately over the Horn of Africa. The correlation between the modified Emanuel‐simulated and Ethiopian station rainfall is quite strong (+0.66). This customized ICTP RegCM3 model now can be used with confidence for the Horn of Africa to study regional rainfall processes and variability and to dynamically downscale seasonal rainfall forecasts. Copyright © 2008 Royal Meteorological Society

Boiling Heat Transfer Rates for Small Precisely Placed Water Droplets on a Heated Horizontal Plate

Two small horizontal surfaces, heated to temperatures up to 220°C, were cooled by small (50–300 mm diameter) room-temperature droplets at 1 atmosphere pressure. One surface was a 10×10 mm thin-film nichrome heater that was used to measure heat fluxes below 100 W/cm2. The other surface, used for fluxes in excess of 100 W/cm2, was a solid copper heater with an 8×8 mm exposed surface. A continuous jet droplet generator coupled with two mutually perpendicular deflection plates was used to manipulate the path of constant diameter water droplets so that the impact of the drops could be precisely located on the heated surfaces. The droplet generator and the deflection plates were employed so that the effect of the impact frequency, droplet diameter, droplet velocity and spacing on the resulting heat transfer rates could be studied under controlled conditions. Optimal droplet spacing between 0.75 and 1.5 times the droplet diameter increased the critical heat flux approximately 30 percent above the value that was achieved when the drops were deposited in one location. For area-averaged mass flow rates less than about 0.08 g/(cm2s), there was no trend in the critical heat flux with the Weber number. However, for larger mass flux rates, the critical heat flux increased with an increasing Weber number. The measured critical heat flux values were roughly twice the heat flux of traditional pool boiling for identical superheat temperatures. Two droplet cooling dimensionless critical heat flux correlations are proposed as a function of Weber and Strouhal numbers; one for a single stream of drops and the other for drops that are spaced across the heated surface. The correlation for the spaced droplets is a function of a dimensionless droplet spacing on the heater.

Criteria for piloted ignition of combustible solids

The onset of piloted ignition (flash point) of combustible polymers is predicted by a gas phase combustion energy density of 1.9 MJ/m 3 that describes the lower flammability limit of fuel vapor–air mixtures. The mass flux (1 g/m 2 s) and (virtual) heat release rate (24 kW/m 2) of the solid at incipient ignition derived from this observation are in general agreement with experimental data for piloted ignition of plastics. Sustained ignition marks the transition from lower-limit to stoichiometric burning and the critical heat release rate increases proportionally to 66 kW/m 2. A critical heat release rate, unique among ignition criteria, is independent of material composition because of the implicit coupling between the gas and condensed phase processes.

Effect of P/M value on the preforms and microstructures of spray formed 70Si30Al alloy

A novel 70Si30Al alloy was prepared by the spray forming process for electronic packaging materials. The effect of the ratio of atomization pressure to metal melt mass flux rate ( PIM) on the preforms and microstructures of the spray-deposited 70Si30Al alloy was studied. The results indicate that the PIM value has a considerable influence on the formation of the preforms and the optimal value is in the range of 0.209–0.231 MPa/(kg-min -1). The microstructure of the spray formed 70Si30Al alloy is fine and homogenous. The primary silicon phases distributing in the aluminum matrix evenly are fine and irregular. The aluminum matrix is divided into wo groups: supersaturated α-Al phase or α-Al phase and Al-Si pseudoeutectic phase or Al-Si eutectic phase.

Recipes for stellar jets: results of combined optical/infrared diagnostics

We examine the conditions of the plasma along a sample of 'classical' Herbig-Haro jets located in the Orion and Vela star forming regions, through combined optical-infrared spectral diagnostics. Our sample includes HH 111, HH 34, HH 83, HH 73, HH 24 C/E, HH 24 J, observed at moderate spatial/spectral resolution. The obtained spectra cover a wide wavelength range from 0.6-2.5 um, including many transitions from regions of different excitation conditions. This allows us to probe the density and temperature stratification which characterises the cooling zones behind the shock fronts along the jet. The derived physical parameters (such as the extinction, the electron density and temperature, the ionisation fraction, and the total density) are used to estimate the depletion onto dust grains of Calcium and Iron with respect to solar abundances. This turns out to be between 70% and 0% for Ca and ~90% for Fe, suggesting that the weak shocks present in the beams are not capable of completely destroying the dust grains. We then derive the mass flux rates (Mdot_jet is on average 5 10^-8 M_solar yr^-1) and the associated linear momentum fluxes. The latter are higher than, or of the same order as, those measured in the coaxial molecular flows, suggesting that the flows are jet driven. Finally, we discuss differences between jets in our sample.

LIQUID FILM SPRAYED ON A STRETCHING CYLINDER

ABSTRACT The fluid dynamics of a film sprayed on a stretching cylinder is solved using similarity transform. The problem is governed by a stretching Reynolds number R and the mass flux rate, which is related to the film thickness. Perturbation solutions for small and large R compare well with exact numerical solutions of the universal curves. The heat transfer problem is discussed.

Lidar characterization of crystalline silica generation and transport from a sand and gravel plant

Light detection and ranging (Lidar) remote sensing two-dimensional vertical and horizontal scans collected downwind of a sand and gravel plant were used to evaluate the generation and transport of geologic fugitive dust emitted by quarry operations. The lidar data give unsurpassed spatial resolution of the emitted dust, but lack quantitative particulate matter (PM) mass concentration data. Estimates of the airborne PM 10 and crystalline silica concentrations were determined using linear relationships between point monitor PM 10 and quartz content data with the lidar backscatter signal collected from the point monitor location. Lidar vertical profiles at different distances downwind from the plant were used to quantify the PM 10 and quartz horizontal fluxes at 2-m vertical resolution as well as off-site emission factors. Emission factors on the order of 65–110 kg of PM 10 (10–30 kg quartz) per daily truck activity or 2–4 kg/t product shipped (0.5–1 kg quartz/t) were quantified for this facility. The lidar results identify numerous elevated plumes at heights >30 m and maximum plume heights of 100 m that cannot be practically sampled by conventional point sampler arrays. The PM 10 and quartz mass flux was greatest at 10–25 m height and decreased with distance from the main operation. Measures of facility activity were useful for explaining differences in mass flux and emission rates between days. The study results highlight the capabilities of lidar remote sensing for determining the spatial distribution of fugitive dust emitted by area sources with intermittent and spatially diverse dust generation rates.

Evaluation of an upscaled model for DNAPL dissolution kinetics in heterogeneous aquifers

Estimates of contaminant fluxes from DNAPL sources as a function of time and DNAPL mass reduction are important to assess the long-term sustainability and costs of monitored natural attenuation and to determine the benefits of partial source removal. We investigate the accuracy of the upscaled mass transfer function (MTF) proposed by Parker and Park [Parker JC, Park E. Modeling field-scale dense nonaqueous phase liquid dissolution kinetics in heterogeneous aquifers. WRR 2004;40:W05109] to describe field-scale dissolved phase fluxes from DNAPL sources for a range of scenarios generated using high-resolution 3-D numerical simulations of DNAPL infiltration and long-term dissolved phase transport. The results indicate the upscaled MTF is capable of accurately describing field-scale DNAPL dissolution rates as a function of time. For finger-dominated source regions, an empirical mass depletion exponent in the MTF takes on values greater than one which results in predicted mass flux rates that decrease continuously with diminishing DNAPL mass over time. Lens-dominated regions exhibit depletion exponents less than one, which results in more step-function like mass flux versus time behavior. Mass fluxes from DNAPL sources exhibiting both lens- and finger-dominated subregions were less accurately described by the simple MTF, but were well described by a dual-continuum model of the same form for each subregion. The practicality of calibrating a dual-continuum model will likely depend on the feasibility of obtaining spatially resolved field measurements of contaminant fluxes or concentrations associated with the subregions using multilevel sampling or some other means.

Shear-induced material transfer across the core-mantle boundary aided by the post-perovskite phase transition

Abstract We present a novel mechanical model for the extraction of outer core material upwards across the CMB into the mantle side region of D” and subsequent interaction with the post-perovskite (ppv) phase transition. A strong requirement of the model is that the D” region behaves as a poro-viscoelastic granular material with dilatant properties. Using new ab-initio estimates of the ppv shear modulus, we show how shear-enhanced dilation promoted by downwelling mantle sets up an instability that drives local fluid flow. If loading rates locally exceed c. 10−12 s−1, calculated core metal upwelling rates are &gt;10−4 m/s, far in excess of previous estimates based on static percolation or capillary flow. Associated mass flux rates are sufficient to deliver 0.5% outer core mass to D” in &lt; 106 yr, provided the minimum required loading rate is maintained. Core metal transported upwards into D” may cause local rapid changes in electrical and thermal conductivity and rheology that if preserved, may account for some of the observed small wavelength heterogeneties (e.g. PKP scattering) there.

Experimental Studies on the Thermal and Moisture Properties of Rammed Earth Used in Adobe Buildings in China

Among the rich diversity of Chinese vernacular buildings, adobe building is a special kind of building using rammed earth for its external walls. The primary aim of this study is to examine quantitatively the thermal and moisture properties of the rammed earth. This paper presents the experimental results of these properties, which include the density, thermal conductivity, specific heat capacity, moisture sorption isotherms and the surface absorption/desorption mass flux rates.

Modeling field-scale multiple tracer injection at a low-level waste disposal site in fractured rocks: Effect of multiscale heterogeneity and source term uncertainty on conceptual understanding of mass transfer processes

Multiple factors may affect the scale-up of laboratory multi-tracer injection into structured porous media to the field. Under transient flow conditions and with multiscale heterogeneities in the field, previous attempts to scale-up laboratory experiments have not answered definitely the questions about the governing mechanisms and the spatial extent of the influence of small-scale mass transfer processes such as matrix diffusion. The objective of this research is to investigate the effects of multiscale heterogeneity, mechanistic and site model conceptualization, and source term density effect on elucidating and interpreting tracer movement in the field. Tracer release and monitoring information previously obtained in a field campaign of multiple, conservative tracer injection under natural hydraulic gradients at a low-level waste disposal site in eastern Tennessee, United States, is used for the research. A suite of two-pore-domain, or fracture-matrix, groundwater flow and transport models are calibrated and used to conduct model parameter and prediction uncertainty analyses. These efforts are facilitated by a novel nested Latin-hypercube sampling technique. Our results verify, at field scale, a multiple-pore-domain, multiscale mechanistic conceptual model that was used previously to interpret only laboratory observations. The results also suggest that, integrated over the entire field site, mass flux rates attributable to small-scale mass transfer are comparable to that of field-scale solute transport. The uncertainty analyses show that fracture spacing is the most important model parameter and model prediction uncertainty is relatively higher at the interface between the preferred flow path and its parent bedrock. The comparisons of site conceptual models indicate that the effect of matrix diffusion may be confined to the immediate neighborhood of the preferential flow path. Finally, because the relatively large amount of tracer needed for field studies, it is likely that source term density effect may exaggerate or obscure the effect of matrix diffusion on the movement of tracers from the preferred flow path into the bedrock.

Internal dynamics of a paleoaccretionary wedge: insights from combined isotope tectonochronology and sandbox modelling of the South-Central Chilean forearc

Forearc accretionary wedges are cyclic systems in which material is frontally and/or basally accreted. Material cycling involves underthrusting, subduction, underplating, exhumation, erosion, transfer to the trench and underthrusting again. In this study we present a novel, tectonochronologic approach to constrain long-term exhumation rates of basally accreted wedge complexes, based on isotopic dating of structural features, on petrological data and sandbox analogue simulations. Congruence between the structural inventory in nature and structures generated in scaled sandbox experiments allows detailed insights into wedge dynamics. For the present-day surface material of the paleoaccretionary wedge of South-Central Chile (Valdivia area, 40°S), published U–Pb ages of detrital zircon place a maximum age of ∼278 Ma for subduction. Prograde metamorphism at transitional greenschist to blueschist facies conditions (420 °C, 8–9 kbar) was immediately followed by progressive penetrative deformation associated with basal accretion, dated at ∼250–245 Ma using Rb/Sr internal mineral isochrons. The accretion process involved duplex tectonics and antiformal stacking, with formation of near-horizontal mylonitic shear zones at around 241 Ma. Continuous basal accretion at depth gave rise to an extensional tectonic regime at higher structural levels. Both semi-ductile, small-scale extensional shear zones and post-kinematic vein mineralizations yield Rb/Sr ages of ∼235 Ma. Tension gashes, representing the latest isotopically dateable stage of structural evolution, were formed at ∼210 Ma, at conditions of ∼230 °C at 1.5–3 kbar, as constrained by fluid inclusion data. Zircon fission track data indicate final cooling to below ∼200 °C at 186±24 Ma. The results suggest continuous basal accretion for at least 50 Ma, with long-term average exhumation rates of 0.6±0.2 mm/a, most probably outbalanced by similar long-term average erosion rates. Changing plate boundary conditions at about 210–200 Ma terminated the accretion process, as evident from a dramatic decrease of exhumation rates at that time. Since then, the paleoaccretionary wedge remained stably in place despite its delicate geotectonic position within the Andean active margin. The tectonochronologic approach complements thermochronologic and geomorphologic methods of exhumation research as it provides direct constraints on mass flux rates even for high-temperature increments of P– T trajectories.

Sorbent Wicking Device for Sampling Hydrophobic Organic Compounds in Unsaturated Soil Pore Water. II: Chemical Capture, Recovery, and Analysis

Experiments investigated the chemical capture and quantitative recovery of chlorinated hydrophobic organic compounds (HOCs) sampled from the pore-water of unsaturated soils with a new in situ sorbent wicking device. The sampling device is comprised of an annular porous stainless-steel interface encasing granular-activated carbon and a fiberglass wick. Laboratory column studies using Ottawa sand spiked with tri- and tetrachlorobenzenes evaluated the accuracy of the sampling device in assessing average pore-water concentrations. The findings from these experiments showed that the sorbent-wick sampler provided reliable and repeatable estimates of the pore-water concentrations of the chlorobenzenes in the unsaturated soil columns at pore-water concentrations as low as several parts per billion. Field trials were conducted in a land biotreatment unit containing industrial soil/sludge materials to assess the availability of trace levels of aqueous-phase polychlorinated biphenyls (PCBs). The land treatment unit had undergone 5 years of intrinsic remediation and the test results indicated that the samplers were able to replicate the aqueous-phase PCB homolog pattern measured in ex situ aqueous extraction tests. Orthosubstituted PCBs were dominant in sorbent-wick samples and aqueous extracts; these PCBs appear relatively recalcitrant in land biotreatment units. The sorbent-wicking device offers the ability to obtain time-averaged mass and volumetric flux rates of HOCs in land treatment units.