Monitoring of Short-Term Heavy Metal Deposition by Accumulation in Epiphytic Lichens (Hypogymnia Physodes (L.) Nyl.)

Many lichens are very sensitive to air pollution due to their symbiotic nature. However, they are generally less sensitive to toxic effects of trace elements; therefore they can be used as accumulator organisms for estimating concentrations of these elements in the environment. Heavy metal accumulation in lichens is a commonly used bioindication method for assessing heavy metal ambient levels. An active biomonitoring method was used for the determination of short-term accumulation of zinc, lead, arsenic and cadmium by epiphytic lichens transplanted at different localities in Slovenia polluted by heavy metals. The atomic absorption spectrophotometry method was used for the determination of heavy metal content in lichen thalli. The content of heavy metals in lichens collected in the background area with clean air (Rogla, Pohorje Mountains) was used as the reference value. Lichens were transplanted from Rogla using the branch transplantation technique for a period of 6 months in the surroundings of Slovenian thermal power plants (Velenje, otanj, Zavodnje, Veliki Vrh, Vnajnarje, and Dobovec) and close to the lead and zinc ore smelter at erjav. The monthly accumulation of heavy metals was comparable within years at selected locations. Heavy metal pollution was the highest at erjav despite of remediation of lead and zinc ore smelter. The monthly accumulation of all four heavy metals was statistically significantly higher in lichens exposed at erjav than at other localities.     

Near-field dispersion from instantaneous sources in the surface layer

This paper considers the near-field dispersion of an ensemble of tracer particles released instantaneously from an elevated source into an adiabatic surface layer. By modelling the Lagrangian vertical velocity as a Markov process which obeys the Langevin equation, we show analytically that the mean vertical drift velocity w(t) is w()=bu _*(1–e ^–(1+)), where is time since release (nondimensionalized with the Lagrangian time scale at the source), b Batchelor's constant, and u _*, the friction velocity. Hence, the mean height and mean depth of the ensemble are calculated. Although the derivation is formally valid only when 1, the predictions for w, mean height and mean depth are consistent in the downstream limit ( 1) with surface-layer Lagrangian similarity theory and with the diffusion equation. By comparing the analytical predictions with numerical, randomflight solutions of the Langevin equation, the analytical predictions are shown to be good approximations at all times, both near-field and far-field.     

Reactions of alkoxy radicals under atmospheric conditions: The relative importance of decomposition versus reaction with O2

The reactions of alkoxy radicals determine to a large extent the products formed during the atmospheric degradations of emitted organic compounds. Experimental data concerning the decompositions, 1,5-H shift isomerizations and reactions with O₂ of several classes of alkoxy radicals are inconsistent with literature estimations of their absolute or relative rate constants. An alternative, although empirical, method for assessing the relative importance under atmospheric conditions of the reactions of alkoxy radicals with O₂ versus decomposition was derived. This estimation method utilizes the differences in the heats of reaction, (H)=(H_{decomposition}–H_{O 2 reaction}), between these two reactions pathways. For (H)[22–0.5(H_{O 2 reaction})], alkoxy radical decomposition dominates over the reaction with O₂ at room temperature and atmospheric pressure of air, while for (H)[25-0.5(H_{O 2 reaction})], the O₂ reaction dominates over decomposition (where the units of H are in kcal mol^–1). The utility and shortcomings of this approach are discussed. It is concluded that further studies concerning the reactions of alkoxy radicals are needed.     

Statistics of Surface-Layer Turbulence Over Terrain with Metre-Scale Heterogeneity

Refuge has patchy vegetation in sandy soil. During midday and at night, the surface sources and sinks for heat and moisture may thus be different. Although the Sevilleta is broad and level, its metre-scale heterogeneity could therefore violate an assumption on which Monin-Obukhov similarity theory (MOST) relies. To test the applicability of MOST in such a setting, we measured the standard deviations of vertical (_w) and longitudinal velocity (_u), temperature (_t), and humidity (_q), the temperature-humidity covariance (¯tq), and the temperature skewness (S_t). Dividing the former five quantities by the appropriate flux scales (u_*, _*, and q_*) yielded the nondimensional statistics _w/u_*, _u/u_*, _t/|t_*|, _q/|q_*|, and ¯tq/t_*q_*. _w/u_*, _t/|t_*|, and S_t have magnitudes and variations with stability similar to those reported in the literature and, thus, seem to obey MOST. Though _u/u_* is often presumed not to obey MOST, our _u/u_* data also agree with MOST scaling arguments. While _q/|q_*| has the same dependence on stability as _t/|t_*|, its magnitude is 28% larger. When we ignore ¯tq/t_*q_* values measured during sunrise and sunset transitions – when MOST is not expected to apply – this statistic has essentially the same magnitude and stability dependence as (_t/t_*)². In a flow that truly obeys MOST, (_t/t_*)², (_q/q_*)², and ¯tq/t_*q_* should all have the same functional form. That (_q/q_*)² differs from the other two suggests that the Sevilleta has an interesting surface not compatible with MOST. The sources of humidity reflect the patchiness while, despite the patchiness, the sources of heat seem uniformly distributed.     

Analytical solutions for the Ekman layer

The PBL equation that governs the transition from the constant-stress surface layer to the geostrophic wind in a neutrally stratified atmosphere for which the eddy viscosityK(z) is assumed to vary smoothly from the surface-layer value U _*z (0.4,U _*=friction velocity,z=elevation) to the geostrophic asymptoteK _GU _*d forzd is solved through an expansion in fd/U _*1 (f=Coriolis parameter). The resulting solution is separated into Ekman's constant-K solution an inner component that reduces to the classical logarithmic form forzd and isO() relative to the Ekman component forzd. The approximationKU _*d is supported by the solution of Nee and Kovasznay's phenomenological transport equation forK(z), which yieldsK–U _*d exp(–z/d), where is an empirical constant for which observation implies, 1. The parametersA andB in Kazanskii and Monin's similarity relation forG/U _* (G=geostrophic velocity) are determined as functions of . The predicted values ofG/U _* and the turning angle are in agreement with the observed values for the Leipzig wind profile. The predicted value ofB based on the assumption of asymptotically constantK is 4.5, while that based on the Nee-Kovasznay model is 5.1; these compare with the observed value of 4.7 for the Leipzig profile. A thermal wind correction, an asymptotic solution for arbitraryK(z) and 1, and an exact (unrestricted ) solution forK(z)=U _*d[1–exp(–z/d)] are developed in appendices.     

Transport-Dissipation Analytical Solutions to the E-∈Turbulence Model and their Role in Predictions of the Neutral ABL

E- turbulence model predictions of the neutralatmospheric boundary layer (NABL) are reinvestigated to determine thecause for turbulence overpredictions found in previous applications. Analytical solutions to the coupled E and equations for the case of steady balance between transport and dissipation terms, the dominant balance just below the NABL top, are derived. It is found that analytical turbulence profiles laminarizeat a finite height only for values of closure parameter ratio c ₂ /_e equal toor slightly greater than one, with laminarization as z for greater . The point = 2 is additionally foundthat where analytical turbulent length scale (l) profilesmade a transition from ones ofdecreasing ( < 2) to increasing ( > 2)values with height. Numerically predicted profiles near the NABL topare consistent with analytical findings. The height-increasingvalues of l predicted throughout the NABL with standard values ofclosure parameters thus appear a consequence of 2.5(> 2), implied by these values (c ₂ = 1.92, _{= 1.3}, _{e = 1}). Comparison of numericalpredictions with DNS data shows that turbulence overpredictions obtained with standard-valued parameters are rectifiedby resetting and _e to 1.1 and 1.6, respectively, giving, with c ₂ = 1.92, 1.3, and laminarization of the NABL's cappingtransport-dissipation region at a finite height.     

Observation of lake breeze penetration and subsequent development of the thermal internal boundary layer for the Nanticoke II Shoreline Diffusion Experiment

The turbulent structure of the lake breeze penetration and subsequent development of the thermal internal boundary layer (TIBL) was observed using a kytoon-mounted ultrasonic anemometer-thermometer. The lake breeze penetrated with an upward rolling motion associated with the upward flow near the lake breeze front. After the lake breeze front passed, the behaviors of the velocity and temperature at the top of the lake breeze layer were similar to those found in convective boundary layers (CBL). Comparing _gq/_*, _u /w _* and _w /w _* between the present observation of TIBL development after the passage of the lake breeze front and CBL data from the literature, the /_* values showed reasonable agreement; however, _u /w _* and _w /W_* had smaller values in the TIBL than in the CBL at higher altitudes. This is due to the differences in the mean velocity profiles. While the CBL has a uniform velocity profile, the TIBL has a peak at lower elevation due to the lake breeze penetration; the velocity then decreases with height.Present address: The Institute of Behavioral Science, 1-35-7 Yoyogi, Tokyo 151, Japan.     

Spatial Variability of Turbulent Fluxes in the Roughness Sublayer of an Even-Aged Pine Forest

The spatial variability of turbulent flow statistics in the roughness sublayer (RSL) of a uniform even-aged 14 m (= h) tall loblolly pine forest was investigated experimentally. Using seven existing walkup towers at this stand, high frequency velocity, temperature, water vapour and carbon dioxide concentrations were measured at 15.5 m above the ground surface from October 6 to 10 in 1997. These seven towers were separated by at least 100m from each other. The objective of this study was to examine whether single tower turbulence statistics measurements represent the flow properties of RSL turbulence above a uniform even-aged managed loblolly pine forest as a best-case scenario for natural forested ecosystems. From the intensive space-time series measurements, it was demonstrated that standard deviations of longitudinal and vertical velocities (_u, _w) and temperature (_T) are more planar homogeneous than their vertical flux of momentum (u_* ²) and sensible heat (H) counterparts. Also, the measured H is more horizontally homogeneous when compared to fluxes of other scalar entities such as CO₂ and water vapour. While the spatial variability in fluxes was significant (>15 %), this unique data set confirmed that single tower measurements represent the canonical structure of single-point RSL turbulence statistics, especially flux-variance relationships. Implications to extending the moving-equilibrium hypothesis for RSL flows are discussed. The spatial variability in all RSL flow variables was not constant in time and varied strongly with spatially averaged friction velocity u_*, especially when u_* was small. It is shown that flow properties derived from two-point temporal statistics such as correlation functions are more sensitive to local variability in leaf area density when compared to single point flow statistics. Specifically, that the local relationship between the reciprocal of the vertical velocity integral time scale (I_w) and the arrival frequency of organized structures (/h) predicted from a mixing-layer theory exhibited dependence on the local leaf area index. The broader implications of these findings to the measurement and modelling of RSL flows are also discussed.     

Determination of similarity functions of the resistance laws for the planetary boundary layer using surface-layer similarity functions

Functional forms of the universal similarity functions A, B (for wind components parallel and normal to the surface stress), and C (for potential temperature difference) are determined based on the generalized theory of the resistance laws for the Planetary Boundary Layer (PBL). The similarity-profile functions for the surface layer are matched with the velocity and temperature-defect profiles that are assumed to have shapes modified by certain powers of nondimensional height z/h, where h is the PBL height. The powers of the outer-layer profile functions are determined, so that the functions become negligible in the surface layer. To close the temperature defect law, an assumption that the temperature gradient across the top of the PBL is continuous with the stratification of the overlying atmosphere is used. The result of this assumption is that nondimensional momentum and temperature profiles in the PBL can be described in terms of four basic ratios: (1) roughness ratio = /h (2) scale-height ratio =|f|h/u_*, (3) ambient stratification parameter =h/_*, and (4) stability parameter =h/L, where L is the Monin-Obukhov length, z₀ is the surface roughness, is the upper-air stratification, u _* is the friction velocity, and _* is the temperature scale at the surface. For stable conditions, the scale-height ratio can be related to the atmospheric stability and the upperair stratification, and the generalized similarity and Rossby number similarity theories become identical. Under appropriate boundary conditions, function A is explicitly dependent on the stability parameter , while B is a function of scale-height ratio , which in turn depends on the stability. Function C is shown to be dependent on the stability and the upper-air stratification, due to the closure assumption used for the temperature profile.The suggested functional forms are compared with other empirical approximations by several authors. The general framework used to determine the functional forms needs to be tested against good boundary-layer measurements.     

Experiments On Buoyant Plume Dispersion In A Laboratory Convection Tank

Plume dispersion in the convective boundary layer (CBL) is investigated experimentally in a laboratory convection tank. The focusis on highly-buoyant plumes that loft near or become trapped in the CBL capping inversion and resistdownward mixing. Such plumes are defined by dimensionless buoyancy fluxes F_* 0.1, where F_* = F_b/(U w_* ² z_i), F_b is the stack buoyancy flux,U is the mean wind speed, w_* is the convective velocity scale, and z_i is the CBL depth. The aim is to obtain statistically-reliable mean (C) and root-mean-square (rms, _c) concentration fields as a function of F_* and the dimensionless distance X = w_*x/(U z_i), where x is the distance downstream of the source.The experiments reveal the following mainresults: (1) For 3 X 4and F_* 0.1, the crosswind-integrated concentration (CWIC) fields exhibit distinctly uniform profiles below z_i with a CWIC maximum aloft, in contrast to the nonuniform profiles obtained earlier by Willis and Deardorff. (2) The lateral dispersion (_y) variation with X is consistent with Taylor's theory for _* 0.1 and a buoyancy-enhanced dispersion, _y/z_i F_* ^1/3X^2/3, forF_* = 0.2 and 0.4. (3) The entrapment, the plume fraction above z_i, has a mean (E) that follows a systematic variationwith X and F_*, and a variability (_e/E) that is broad ( 0.3 to 2) near the source but subsides to 0.25 far downstream. (4) Vertical profiles of the concentration fluctuation intensity (c/C) are uniform for z < z_i and X > 1.5, but exhibit significant increases: (a) at the surface and close to the source (X 1.5), and(b) in the entrainment zone. (5) The cumulative distribution functions (CDFs) of the scaled concentration fluctuations (c/_c) separate into mixed-layer and entrainment-layer CDFs for X 2, with the mixed-layer group collapsing to a single distribution independent of z.These are the first experiments to obtain all components of the lateral and vertical dispersion parameters (rms meander, relative dispersion, total dispersion) for continuous buoyant releases in a convection tank. They also are the first tank experiments to demonstrate agreement with field observations of: (1) the scaled ground-level concentration along the plume centreline, and (2) the dimensionless lateral dispersion _y/z_i of buoyant plumes.     

Turbulence in a Katabatic Flow

Turbulence measurements performed in a stable boundary layer over the sloping ice surface of the Vatnajökull in Iceland are described. The boundary layer, in which katabatic forces are stronger than the large-scale forces, has a structure that closely resembles that of a stable boundary layer overlying a flat land surface, although there are some important differences. In order to compare the two situations the set-up of the instruments on an ice cap in Iceland was reproduced on a flat grass surface at Cabauw, the Netherlands. Wind speed and temperature gradients were calculated and combined with flux measurements made with a sonic anemometer in order to obtain the local stability functions _m and _h as a function of the local stability parameter z/L. Unlike the situation at Cabauw, where _m was linear as a function of z/L, in the katabatically forced boundary layer, the dependence of _m on stability was found to be non-linear and related to the height of the wind maximum. Thermal stratification and the depth of the stable boundary layer however seem to be rather similar under these two different forcing conditions.Furthermore, measurements on the ice were used to construct the energy balance. These showed good agreement between observed melt and components contributing to the energy balance: net radiation (supplying 55% of the energy), sensible heat flux (30%) and latent heat flux (15%).Local sources and sinks in the turbulent kinetic energy budget are summed and indicate a reasonable balance in near-neutral conditions but not in more stable situations. The standard deviation of the velocity fluctuations _u, _v, and _w, can be scaled satisfactorily with the local friction velocity u_* and the standard deviation of the temperature fluctuation with the local temperature scale _*.     

Numerical Simulation of Strongly Stratified Flow Over a Three-Dimensional Hill

A numerical study of stably stratified flow over a three-dimensional hill is presented. Large-eddy simulation is used here to examine in detail the laboratory experimental flows described in the landmark work of Hunt and Snyder about stratified flow over a hill. The flow is linearly stratified and U/Nh is varied from 0.2 to 1.0. Here N and U are the buoyancy frequency and freestream velocity respectively, and h is the height of the hill. The Reynolds number based on the hill height is varied from 365 to 2968. The characteristic flow patterns at various values of U/Nh have been obtained and they are in good agreement with earlier theoretical and experimental results. It is shown that the flow field cannot be predicted by Drazin's theory when recirculation exists at the leeside of the hill even at UNh 1. The wake structure agrees well with a two-dimensional wake assumption when U/Nh 1 but lee waves start to influence the wake structure as U/Nh increases. The dividing-streamline heights obtained in the simulation are in accordance with experimental results and Sheppard's formula. The energy loss along the dividing streamline due to friction/turbulence approximately offsets the energy gained from pressure field. When lee waves are present, linear theory always underestimates the amplitude and overestimates the wavelength of three-dimensional lee waves. The simulated variations of drag coefficients with the parameterK (=ND/ U) are qualitatively consistent with experimental data and linear theory. Here D is the depth of the tank.     

The use of σ T as a temperature scale in the atmospheric surface layer

The standard deviation of temperature _T is proposed as a temperature scale and as a velocity scale to describe the behaviour of turbulent flows in the Atmospheric Surface Layer (ASL), instead of _* andu _* of the Monin—Obukhov similarity theory, and ofT _f andU _f used for free convection stability conditions. On the basis of experimental evidence reported in the literature, it is shown that _T T _f andv _* U _f in the free convection region, and _{T *} andv _* U _* in nearneutral and stable conditions. This implies that the proposed scales can be applied for all stabilities. Furthermore, a new length scale is proposed and its relation with Obukhov length is given. Also, a simple semi-empirical expression is presented with which _T andv _* can be evaluated in a rather simple way. Some examples of practical applications are given, e.g., a stability classification for unstable conditions.     

Microwave vegetation optical depth and inverse modelling of soil emissivity using Nimbus/SMMR satellite observations

Summary A radiative transfer model has been used to determine the large scale effective 6.6 GHz and 37 GHz optical depths of the vegetation cover. Knowledge of the vegetation optical depth is important for satellite-based large scale soil moisture monitoring using microwave radiometry. The study is based on actual observed large scale surface soil moisture data and observed dual polarization 6.6 and 37 GHz Nimbus/SMMR brightness temperatures over a 3-year period. The derived optical depths have been compared with microwave polarization differences and polarization ratios in both frequencies and with Normalized Difference Vegetation Index (NDVI) values from NOAA/AVHRR. A synergistic approach to derive surface soil emissivity from satellite observed brightness temperatures by inverse modelling is described. This approach improves the relationship between satellite derived surface emissivity and large scale top soil moisture fromR ²=0.45 (no correction for vegetation) toR ²=0.72 (after correction for vegetation). This study also confirms the relationship between the microwave-based MPDI and NDVI earlier described and explained in the literature.List of Symbols f frequency [Hz] - f _i(p) fractional absorption at polarizationp - h surface roughness - h h cos² - H horizontal polarization - n _i complex index of refraction - p polarization (H orV) - R _s microwave surface reflectivity - T _B(p) brightness temperature at polarizationp - T ^* normalized brightness temperature - T polarization difference (T _v-T _H) - T _s temperature of soil surface - T _c temperature of canopy - T _max daily maximum air temperature - T _min daily minimum air temperature - V vertical polarization - soil moisture distribution factor; also used for the constant to partition the influence of bound and free water components to the dielectric constant of the mixture - empirical complex constant related to soil texture - microwave transmissivity of vegetation (=e ^–) - ^* effective transmissivity of vegetation (assuming =0) - microwave emissivity - _s emissivity of smooth soil surface - _rs emissivity of rough soil surface - _vs emissivity of vegetated surface - soil moisture content (% vol.) - K dielectric constant [F·m^–1] - K _fw dielectric constant of free water [F·m^–1] - K _ss dielectric constant of soil solids [F·m^–1] - K _m dielectric constant of mixture [F·m^–1] - K _o permittivity of free space [8.854·10^–12 F·m^–1] - high frequency limit ofK _wf [F·m^–1] - wavelength [m] - incidence angle [degrees from nadir] - polarization ratio (T _H/T _V) - b soil bulk density [gr·cm^–3] - s soil particle density [gr·cm^–3] - R surface reflectivity in red portion of spectrum - NIR surface reflectivity in near infrared portion of spectrum - _eff effective conductivity of soil extract [mS·cm^–1] - vegetation optical depth - _6.6 vegetation optical depth at 6.6 GHz - ₃₇ vegetation optical depth at 37 GHz - ^* effective vegetation optical depth (assuming =0) - single scattering albedo of vegetation With 12 Figures     

Equilibrium evaporation beneath a growing convective boundary layer

Expressions for the equilibrium surface Bowen ratio ( _s ) and equilibrium evaporation are derived for a growing convective boundary layer (CBL) in terms of the Bowen ratio at the top of the mixed layer _i and the entrainment parameter A _R . If A_R is put equal to zero, the solution for _s becomes-that previously obtained for the zero entrainment or closed box model. The Priestley-Taylor parameter is also calculated and plotted in terms ofA _R and _i . Realistic combinations of the atmospheric parameters give values of in the range 1.1 to 1.4.     

Final results of the CONDORS convective diffusion experiment

TheConvectiveDiffusionObserved byRemoteSensors (CONDORS) field experiment conducted at the Boulder Atmospheric Observatory used innovative techniques to obtain three-dimensional mappings of plume concentration fields, /Q, of oil fog detected by lidar and chaff detected by Doppler radar. It included extensive meteorological measurements and, in 1983, tracer gases measured at a single sampling arc. Final results from ten hours of elevated and surface release data are summarized here. Many intercomparisons were made. Oil fog /Q measured 40m above the arc are mostly in good agreement withSF ₆ values, except in a few instances with large spacial inhomogeneities over short distances. After a correction scheme was applied to compensate for the effect of its settling speed, chaff dy/Q agreed well with those of oil except in two cases of oil fog hot spots. Mass or frequency distribution vs. azimuth or elevation angle comparisons were made for chaff, oil, and wind, with mostly good agreements. Spacial standard deviations, _y and _z, of chaff and oil agree overall and are consistent at short range with velocity standard deviations _vand _w 0.6w^* (the convective scale velocity), as measured atz>100m. Surface release _y is enhanced up to 60% at smallx, consistent with the Prairie Grass measurements and with larger _v and reduced wind speed measured near the surface. Decreased _y at small dimensionless average times is also noted. Finally, convectively scaled dy, C _y, were plotted versus dimensionlessx andz for oil, chaff, and corrected chaff for each 30–60 min period. Aggregated CONDORSC _y fields compare well with laboratory tank and LES numerical simulations; surface-released oil fog compares expecially well with the tank experiments. However, large deviations from the norm occurred in individual averaging periods; these deviations correlated strongly with anomalies in measured distributions.On assignment to the US Environmental Protection Agency, Atmospheric Research and Exposure Assessment Laboratory, RTP, NC.     

The behaviour of a mixed-layer model of the convective boundary layer coupled to a big leaf model of surface energy partitioning

Daily mean values of the Priestley-Taylor coefficient, ¯, are derived from a simple model of the daily growth of a convective boundary layer. For a particular control set of driving environmental variables, ¯ is related to the prescribed bulk surface resistance, r_S by 1/¯ = 1/₀ + mr_S for parameters ₀ and m. The dependence of the parameters ₀ and m on weather is explored and a potential use of this linear relation to provide information about regional values of r_S is indicated.     

Local Similarity Relationships In The Urban Boundary Layer

To investigate turbulent structures in an urban boundary layer (UBL) with many tallbuildings, a number of non-dimensional variable groups based on turbulent observationsfrom a 325-m meteorological tower in the urban area of Beijing, China, are analyzedin the framework of local similarity. The extension of surface-layer similarity to localsimilarity in the stable and unstable boundary layer is also discussed. According to localsimilarity, dimensionless quantities of variables: e.g., velocity and temperature standarddeviations _i/u_*l (i=u,v,w) and_T/T_*l,correlation coefficients of uw and wT covariance, gradients of wind and temperature_m and _h, and dissipation rates of turbulent kinetic energy (TKE) andtemperature variance and _N can be represented as a functiononly of a local stability parameter z/, where is the local Obukhovlength and z is the height above ground. The average dissipation rates of TKE andtemperature variance are computed by using the u spectrum, and the uw and wTcospectra in the inertial subrange. The functions above were found to be in a goodagreement with observational behaviour of turbulence under unstable conditions, butthere were obvious differences in the stable air.     

Scavenging Efficiency of ‘Aerosol Carbon’ and Sulfate in Supercooled Clouds at Mt. Sonnblick (3106 m a.s.l., Austria)

Cloud water and interstitial aerosol samples collected at Mt. Sonnblick (SBO) were analyzed for sulfate and aerosol carbon to calculate in-cloud scavenging efficiencies. Scavenging efficiencies for sulfate (_SO) ranged from 0.52 to 0.99 with an average of 0.80. Aerosol carbon was scavenged less efficiently with an average value (_AC) of 0.45 and minimum and maximum values of 0.14 and 0.81, respectively. Both _SO and _AC showed a marked, but slightly different, dependence on the liquid water content (LWC) of the cloud. At low LWC, _SO increased with rising LWC until it reached a relatively constant value of 0.83 above an LWC of 0.3 g/m³. In the case of aerosol carbon, we obtained a more gradual increase of _AC up to an LWC of 0.5 g/m³. At higher LWCs, _{_} remained relatively constant at 0.60. As the differences between _SO and _A varied across the LWC range observed at SBO, we assume that part of the aerosol carbon was incorporated into the cloud droplets independently from sulfate. This hypothesis is supported by size classified aerosol measurements. The differences in the size distributions of sulfate and total carbon point to a partially external mixture. Thus, the different chemical nature and the differences in the size and mixing state of the aerosol particles are the most likely candidates for the differences in the scavenging behavior.