Local and distant residence times of contaminants in multi-compartment models. Part II: application to assessing environmental mobility and long-range atmospheric transport

In Part I, the concepts of inherent, local and distant residence times (DRTs) were reviewed as metrics of the extent to which chemical discharges or emissions in one region or box are transported to distant regions. In this second part, the concepts are applied to geographically relevant systems to illustrate their applicability to the assessment of chemicals for long-range transport potential (LRTP). It is shown that the relative ranking of chemicals as characterized by the DRT method is similar to that of the characteristic travel distance concept. A DRT source-receptor matrix is developed that can express the chemical-specific potential of source regions to contaminate a specific receptor region of concern such as the Arctic. The matrix can be modified to identify for a specific source region the likely destinations of emissions as well as to assess the relative vulnerability of regions in the global environment to contaminants of concern.     

Validation of trajectory statistical methods

A number of works have been trying to validate various trajectory statistical methods (TSMs), mostly through subjective comparison with known sources. Here in a more comprehensive and quantitative approach three trajectory statistical methods (potential source contribution function (PSCF), concentration field method (CF), and redistributed concentration field method (RCF)) were subjected to two validation approaches: validation with virtual and real sources under idealised conditions, where the effects of dispersion and removal of the trace substance are excluded, and comparisons with the EMEP SO₂ emission inventory under realistic conditions.The best performance was achieved in an idealised situation with about 78% common spatial variance between the EMEP emission inventory and the trajectory statistical reconstruction of the EMEP emission inventory with the RCF method, whereas the real world experiments for SO₂ on an European scale resulted in a much lower performance with 33% common spatial variance between the EMEP SO₂ emission inventory and the trajectory statistical reconstruction with the PSCF method.The experiments suggest that the limitation of the accuracy and spatial range of TSMs are rooted in the simplified transport process described just by trajectory paths. If one links these limitations with the concept of the mean residence time of the considered trace substance, a temporal and spatial scope can be deduced, within which the effect of the simplification of the transport process is restricted and useful information can be expected from TSMs. The lower values of the mean residence time for SO₂ range from 9 to 17 h, which were deduced from the decay approach, where an exponential decay, respectively, removal of SO₂ was built into the trajectory statistical procedure. The values derived from the optimum real world validation experiment place the upper range of the mean residence time to about 60 h or 2.5 days. Both figures are within the range of mean residence times for SO₂ cited in literature. Through the validation experiments of this work the rule of thumb, not to trust TSMs beyond the mean residence time of the substance, has become palpable. Nevertheless TSMs and related methods are computationally fast procedures, which deliver first hints on potential source areas, if applied within the frame of the mean residence time of the considered substance.     

Volatile organic compound emissions from municipal solid waste disposal sites: A case study of Mumbai,India

Improper solid waste management leads to aesthetic and environmental problems. Emission of volatile organic compounds (VOCs) is one of the problems from uncontrolled dumpsite. VOCs are well known to be hazardous to human health and many of them are known or potential carcinogens. They also contribute to ozone formation at ground level and climate change as well. The qualitative and quantitative analysis of VOCs emitting from two municipal waste (MSW) disposal sites in Mumbai, India, namely Deonar and Malad, are presented in this paper. Air at dumpsites was sampled and analyzed on gas chromatography–mass spectrometry (GC-MS) in accordance with U.S. Environmental Protection Agency (EPA) TO-17 compendium method for analysis of toxic compounds. As many as 64 VOCs were qualitatively identified, among which 13 are listed under Hazardous Air Pollutants (HAPs). Study of environmental distribution of a few major VOCs indicates that although air is the principal compartment of residence, they also get considerably partitioned in soil and vegetation. The CO₂ equivalent of target VOCs from the landfills in Malad and Deonar shows that the total yearly emissions are 7.89E+03 and 8.08E+02 kg, respectively. The total per hour ozone production from major VOCs was found to be 5.34E-01 ppb in Deonar and 9.55E-02 ppb in Malad. The total carcinogenic risk for the workers in the dumpsite considering all target HAPs are calculated to be 275 persons in 1 million in Deonar and 139 persons in 1 million in Malad.

Implications: This paper describes the hazards of VOC emission from open dumpsites, a common practice, in an Indian metro city. The subsequent partitioning of the emitted VOCs in other environmental compartment from air is presented. The global warming potential and the health hazards to the dumpsite workers from the emitted VOCs have also been estimated.     

Kinematic analysis of solute mass flows in rock fractures with spatially random parameters

Field data of physical properties in heterogeneous crystalline bedrock, like porosity and fracture aperture, is associated with uncertainty that can have a significant impact on the analysis of solute transport in rock fractures. Solutions to the central temporal moments of the residence time probability density function (PDF) are derived in a closed form for a solute Dirac pulse. The solutions are based on a model that takes into account advection along the fracture plane, diffusion into the rock matrix and sorption kinetics in the rock matrix. The most relevant rock properties including fracture aperture and several matrix properties as well as flow velocity are assumed to be spatially random along transport pathways. The mass transport is first solved in a general form along one-dimensional pathways, but the results can be extended to multi-dimensional flows simply by substituting the expected travel time for inert water parcels. Based on data obtained with rock samples taken at Asp? Hard Rock Laboratory in Sweden, the solutions indicate that the heterogeneity of the rock properties contributes to increasing significantly both the variance and the skewness of the residence time probability density function for a pulse travelling in a fracture. The Asp? data suggests that the bias introduced in the variance of the residence time PDF by neglecting the effect of heterogeneity of the rock properties on the radionuclide migration is very large for fractures thinner than a few tenths of a millimetre.     

Indicators for persistence and long-range transport potential as derived from multicompartment chemistry-transport modelling

Total environmental and compartmental residence times as a measure for persistence as well as indicators for long-range transport potential (LRTP) have been derived from global geo-referenced modelling and LRTP is characterized in two geographic directions for the first time. A dynamic multicompartment chemistry-transport model (MCTM) was used to study the fate of the insecticides DDT and alpha-hexachlorocyclohexane (alpha-HCH) during the first 2 years upon entry. The indicators for LRTP were defined such as to address the tendencies of substance distributions to migrate ('plume displacement', PD) and to spread into remote areas ('spatial spreading', SS). The indicators deliver values as function of time upon entry. With the aim to address the effect of location of entry on environmental fate, scenarios of emission from the territories of seven countries were studied. It was found that the effect of location of entry on the spatial scale of countries (400-4000 km) is significant for the compartmental distribution and the inter-compartmental mass exchange fluxes (e.g., number of atmospheric cycles, 'hops'). Location of entry introduces uncertainties in the order of a factor of 5 for the total environmental residence time, tau(overall), and a factor of 5-20 for PD and SS. For the 2nd year upon entry into the environment, tau(overall)=317-1527 days are predicted for DDT and 101-463 days for alpha-HCH. The influence of location of entry does affect the substance ranking, i.e. we cannot simply state that DDT is more persistent than alpha-HCH, but for one scenario studied, application in China, the opposite is predicted. Precipitation patterns proved to be significant, besides other climate parameters, for atmospheric residence time. Integration of the location of entry in chemicals risk assessments is therefore recommended. In general, persistence and some indicators for LRTP, pertinent to their definition, refer to the fate of a large fraction, e.g., 63% (=1-1/e) or 90%, but not the total substance burden. The choice of this fraction may have the consequence of a normative step which defines the spatial and temporal extensions of a related chemicals risk assessment and may affect substance ranking.     

Estimating the relevance of engineered carbonaceous nanoparticle facilitated transport of hydrophobic organic contaminants in porous media

Naturally occurring nanoparticles (NP) enhance the transport of hydrophobic organic contaminants (HOCs) in porous media. In addition, the debate on the environmental impact of engineered nanoparticles (ENP) has become increasingly important. HOC bind strongly to carbonaceous ENP. Thus, carbonaceous ENP may also act as carriers for contaminant transport and might be important when compared to existing transport processes. ENP bound transport is strongly linked to the sorption behavior, and other carbonaceous ENP-specific properties. In our analysis the HOC-ENP sorption mechanism, as well as ENP size and ENP residence time, was of major importance. Our results show that depending on ENP size, sorption kinetics and residence time in the system, the ENP bound transport can be estimated either as (1) negligible, (2) enhancing contaminant transport, or (3) should be assessed by reactive transport modeling. One major challenge to this field is the current lack of data for HOC-ENP desorption kinetics.     

Generic scenario for the impact of cooling lubricants into the hydrosphere

According to directive 93/67/EEC of the European Commission, this paper deals with the generic exposition evaluation of the environmental concentration of cooling Lubricant chemicals from the metal working industry into the hydrosphere. After considering the relevant life-cycle steps and the selection of a representative point source for the 'reasonable worst case', the emission per day Elocal(water) is determined. It leads to the predicted environmental concentration (PEClocal(water)) for the local stage in the compartment water. In order to gain the PEClocal(water) for an example--a corrosion inhibitor as additive in a cooling lubricant--, the relevant emission paths and the corresponding representative point source are described for the reasonable worst case. For non-water-miscible cooling lubricants, none of the operations in the life cycle leads to a release into the compartment water. To evaluate the hazard potential for cooling lubricant chemicals, the complete risk assessment has to be done. Also, the assessment has to be done for all high production volume chemicals, new substances and existing hazard chemicals. This means that even industrial categories like chemicals used in the textile industry or biocides and others have to be evaluated.     

Human population intake fractions and environmental fate factors of toxic pollutants in life cycle impact assessment

The present paper outlines an update of the fate and exposure part of the fate, exposure and effects model USES-LCA. The new fate and exposure module of USES-LCA was applied to calculate human population intake fractions and fate factors of the freshwater, marine and terrestrial environment for 3393 substances, including neutral organics, dissociating organics and inorganics, emitted to 7 different emission compartments. The human population intake fraction is on average 10(-5)-10(-8) for organics and 10(-3)-10(-4) for inorganics, depending on the emission compartment considered. Chemical-specific human population intake fractions can be 1-2.7 orders of magnitude higher or lower compared to the typical estimates. For inorganics, the human population intake fractions highly depend on the assumption that exposure via food products can be modelled with constant bioconcentration factors. The environmental fate factor is on average 10(-11)-10(-18) days m(-3) for organics and 10(-10)-10(-12) days m(-3) for inorganics, depending on the receiving environment and the emission compartment considered. Chemical-specific environmental fate factors can be 1-8 orders of magnitude higher or lower compared to the typical estimates. The largest differences between the new and old version of USES-LCA are found for emissions to air and soil. This is caused by a significant change in the structure of the air and soil compartments in the new version of USES-LCA, i.e. the distinction between rural and urban air, including rain-no rain conditions and including soil depth dependent intermedia transport.     

Combustion sources of particles. 1. Health relevance and source signatures

Combustion processes result in generation of a large number of particle and gaseous products that create health and environmental risks. Of particular importance are the very small particles that are emitted in large quantities from all the combustion sources, and that have been shown to be potentially more significant in terms of their impact on health than larger particles. To control and mitigate the particles with a view of health and environmental risk reduction, a good understanding is necessary of the relative and absolute contribution from the emission sources to the airborne concentrations. This understanding could only be achieved by developing source signature libraries through direct emission measurements from the sources on one hand, and by measuring particle concentrations in the air, and apportioning them to the specific local and distant sources using the signatures, on the other hand. This paper is a review of particle characteristics that are used as source signatures as well as their general advantages and limitations. The second part of the paper reviews source signatures of the most common combustion pollution sources.     

Multicompartmental fate of persistent substances

Background, Aim and Scope Modelling of the fate of environmental chemicals can be done by relatively simple multi-media box models or using complex atmospheric transport models. It was the aim of this work to compare the results obtained for both types of models using a small set of non-ionic and non-polar or moderately polar organic chemicals, known to be distributed over long distances. Materials and Methods Predictions of multimedia exposure models of different types, namely three multimedia mass-balance box models (MBMs), two in the steady state and one in the non-steady state mode, and one non-steady state multicompartment chemistry-atmospheric transport model (MCTM), are compared for the first time. The models used are SimpleBox, Chemrange, the MPI-MBM and the MPI-MCTM. The target parameters addressed are compartmental distributions (i.e. mass fractions in the compartments), overall environmental residence time (i.e. overall persistence and eventually including other final sinks, such as loss to the deep sea) and a measure for the long-range transport potential. These are derived for atrazine, benz-[a]-pyrene, DDT, α and γ-hexachlorocyclohexane, methyl parathion and various modes of substance entry into the model world. Results and Discussion Compartmental distributions in steady state were compared. Steady state needed 2–10 years to be established in the MCTM. The highest fraction of the substances in air is predicted by the MCTM. Accordingly, the other models predict longer substance persistence in most cases. The results suggest that temperature affects the compartmental distribution more in the box models, while it is only one among many climate factors acting in the transport model. The representation of final sinks in the models, e.g. burial in the sediment, is key for model-based compartmental distribution and persistence predictions. There is a tendency of MBMs to overestimate substance sinks in air and to underestimate atmospheric transport velocity as a consequence of the neglection of the temporal and spatial variabilities of these parameters. Therefore, the long-range transport potential in air derived from MCTM simulations exceeds the one from Chemrange in most cases and least for substances which undergo slow degradation in air. Conclusions and Perspectives MBMs should be improved such as to ascertain that the significance of the atmosphere for the multicompartmental cycling is not systematically underestimated. Both types of models should be improved such as to cover degradation in air in the particle-bound state and transport via ocean currents. A detailed understanding of the deviations observed in this work and elsewhere should be gained and multimedia fate box models could then be ‘tuned in’ to match better the results of comprehensive multicompartmental transport models. ESS-Submission Editor: Prof. Dr. Michael Matthies (matthies@uos.de)     

Emission of respirable dust from stone quarrying,potential health effects,and its management

Mining of minerals exerts adverse pressure on different compartments of environment directly or indirectly. Air is the worst affected environmental matrix, and it can carry the harmful effect of pollutants generated from mining activity even to distant places. The present study was undertaken to estimate the emission of particulate matter (PM_2.5 and PM₁₀) from different activities undertaken in stone quarrying in Mahendragarh, Haryana. The results obtained from the present study indicated that drilling, blasting, crushing, and transport of mined material are chiefly responsible for the generation of dust. Whereas drilling, blasting, and loading were responsible for emission of higher fraction of PM₁₀, crushing and re-suspension of roadside dust from movement of vehicles resulted in generation of relatively higher fraction of finer dust (PM_2.5). Modelling the transport of dust over the Hybrid Single-Particle Lagrangian Integrated Trajectory model revealed that the emitted particle may move up to the distance of about 40 km within 4 h of emission under average meteorological conditions. Fourier transform infrared (FTIR) spectroscopy analysis of dust confirmed the presence of calcite and gypsum, thus confirming the source as mining. The study concluded that generation of PM_2.5-sized particles may impose serious respiratory health effects over the workers engaged in mining, crushing, and transportation of sandstone. Apart from it, population residing downwind of the mining area is particularly vulnerable to the pulmonary effects due to inhalation of dust.

     

Surface ozone at the Swiss Alpine site Arosa: the hemispheric background and the influence of large-scale anthropogenic emissions

An innovative and effective method using isentropic trajectory analysis based on the residence time of air masses over the polluted region of Europe was successfully applied to categorize surface ozone amounts at Arosa, Switzerland during 1996–1997. The “European representative” background ozone seasonal cycle at Arosa is associated with long-range transport of North Atlantic air masses, and displays the spring maximum–summer minimum with an annual average of 35 ppb. The photochemical ozone production due to the intense large-scale anthropogenic emission over Europe is estimated as high as 20 ppb in summer, whereas it is insignificant in winter. European sources contribute an annual net ozone production of 9–12 ppb at Arosa. Comparison with the selected regional representative site in Western Europe shows similar results indicating that the categorized ozone data at Arosa by this technique could be regarded as a representative for northern hemispheric mid-latitudes.     

Quantifying the global fractionation of polychlorinated biphenyls

Due to the wide range of their physical-chemical properties, polychlorinated biphenyls (PCBs) have played an important role in the derivation of the global fractionation hypothesis, which predicts changes in the composition of persistent organic pollutant mixtures with latitude. Recent historical emission estimates, the derivation of an internally consistent property data set, in combination with a zonally averaged global fate and transport model, allow a quantitative investigation of the compositional shifts PCBs experience as a function of environmental compartment, latitude and time. Model simulations reproduce the higher relative abundance of lighter PCB congeners with increasing latitude, observed in air and soil, and quantify the relative importance of partitioning, persistence and emissions in establishing PCB patterns. Compositional variations consistent with global fractionation, as well as inverted concentration profiles with higher levels in the Arctic than at lower latitudes, are consistent with only minor fractions of the global PCB inventory being transferred northward.     

The reactive transport of trichloroethene is influenced by residence time and microbial numbers

The dechlorination rate in a flow-through porous matrix can be described by the species specific dechlorination rate observed in a liquid batch unless mass transport limitations prevail. This hypothesis was examined by comparing dechlorination rates in liquid batch with that in column experiments at various flow rates (3-9-12 cm day(-1)). Columns were loaded with an inoculated sand and eluted with a medium containing 1mM trichloroethene (TCE) for 247 days. Dechlorination in the column treatments increased with decreasing flow rate, illustrating the effect of the longer residence time. Zeroth order TCE or cis-DCE degradation rates were 4-7 folds larger in columns than in corresponding batch systems which could be explained by the higher measured Geobacter and Dehalococcoides numbers per unit pore volume in the columns. The microbial numbers also explained the variability in dechlorination rate among flow rate treatments marked by a large elution of the dechlorinating species' yield as flow increased. Stop flow events did not reveal mass transport limitations for dechlorination. We conclude that flow rate effects on reactive transport of TCE in this coarse sand are explained by residence time and by microbial transport and that mass transport limitations in this porous matrix are limited.     

Atmospheric lead and bromine in Germany: post-abatement levels,variabilities and trends

INTENTION, GOAL, SCOPE, BACKGROUND: Abatement measures since the 1970s have depleted lead and bromine levels in the atmosphere over large parts of Europe. Our knowledge of the atmospheric cycling of these elements while a several decade-long period of intensive mobilization reaches its end is incomplete. OBJECTIVE: We have characterized the trends in the atmospheric levels of Pb and Br and present-day temporal and spatial variabilities. METHODS: This was achieved by short-term (weeks) and long-term (years) measurements of particulate Pb and Br at various sites in Germany. Samples of atmospheric particulate matter were collected on filter membranes and analyzed by x-ray fluorescence. RESULTS AND DISCUSSION: Average Pb levels at rural and urban inland sites did not exceed concentrations in background aerosols, sampled at a Baltic Sea coastal site, by more than a factor of 3. Due to sea salt, bromide inland particulate Br levels are below those at a coastal site. There, non-sea salt Br, however, is significant as well. Urban Pb and Br levels are not necessarily higher than rural levels. The concentration levels have decreased in such a way that the previously common source, local vehicular traffic emissions, is no longer predominant. Regional rather than local sources have been increasingly determining the concentrations since the 1990s. This is more pronounced for Br than for Pb. We found indications for coal burning and long-range transport as significant Pb sources. For particulate Pb species, a range of ages (elapsed time since Pb emission) has been found. This range shows two maxima corresponding to characteristic times of 72 and 24 h. CONCLUSIONS: The (mean) atmospheric residence time of particulate Pb is longer than the residence time of Br, in particular in the wintertime. The chemical species contributing to atmospheric Pb should be addressed in future studies. RECOMMENDATION AND OUTLOOK: Clearly, despite effective abatement measures, atmospheric Pb will continue to be dominated by anthropogenic mobilization. The influence from long-range transport can be expected to decrease with the effectiveness of abatement programmes in neighbouring countries of the region.     

Modeling water exchange and contaminant transport through a Baltic coastal region

The water exchange of the Baltic coastal zone is characterized by its seasonally varying regimes. In the safety assessment of a potential repository for spent nuclear fuel, it is important to assess the consequences of a hypothetical leak of radionuclides through the seabed into a waterborne transport phase. In particular, estimates of the associated residence times in the near-shore coastal zone are of interest. There are several methods to quantify such measures, of which three are presented here. Using the coastal location of Forsmark (Sweden) as an example, methods based on passive tracers, particle trajectories, and the average age distribution of exogeneous water parcels are compared for a representative one-year cycle. Tracer-based methods can simulate diffusivity more realistically than the other methods. Trajectory-based methods can handle Lagrangian dispersion processes due to advection but neglect diffusion on the sub-grid scale. The method based on the concept of average age (AvA) of exogeneous water can include all such sources simultaneously not only boundary water bodies but also various (fresh)-water discharges. Due to the inclusion of sub-grid diffusion this method gives a smoother measure of the water renewal. It is shown that backward in time trajectories and AvA-times are basically equipollent methods, yielding correlated results within the limits set by the diffusivity.     

Hazardous air pollutants in industrial area of Mumbai - India

Hazardous Air Pollutants (HAPs) have a potential to be distributed into different component of environment with varying persistence. In the current study fourteen HAPs have been quantified in the air using TO-17 method in an industrial area of Mumbai. The distribution of these HAPs in different environmental compartments have been calculated using multi media mass balance model, TaPL3, along with long range transport potential and persistence. Results show that most of the target compounds partition mostly in air. Phenol and trifluralin, partition predominantly into soil while ethyl benzene and xylene partition predominantly into vegetation compartment. Naphthalene has the highest persistence followed by ethyl benzene, xylene and 1,1,1 trihloro ethane. Long range transport potential is maximum for 1,1,1 trichloroethane. Assessment of human health risk in terms of non-carcinogenic hazard and carcinogenic risk due to exposure to HAPs. have been estimated for industrial workers and residents in the study area considering all possible exposure routes using the output from TaPL3 model. The overall carcinogenic risk for residents and workers are estimated as high as unity along with very high hazard potential.     

Global environmental cycling of γ-HCH and DDT in the 1980s – A study using a coupled atmosphere and ocean general circulation model

A coupled atmosphere–ocean general circulation model, ECHAM5-MPIOM, was used to study the multicompartmental cycling and long-range transport of persistent and semivolatile organics. Multiphase systems in air and ocean are covered by submodels for atmospheric aerosols, HAM, and marine biogeochemistry, HAMOCC5, respectively. The model, furthermore, encompasses 2D surface compartments, i.e. top soil, vegetation surfaces and sea-ice. The total environmental fate of γ-hexachlorocyclohexane (γ-HCH, lindane) and dichlorophenyltrichloroethane (DDT) in agriculture were studied.DDT is mostly present in the soils, the water-soluble γ-HCH in soils and ocean. DDT has the longest residence time in almost all compartments. Quasi-steady state with regard to substance accumulation is reached within a few years in air and vegetation surfaces. In seawater the partitioning to suspended and sinking particles contributes to the vertical transport of substances. On the global scale deep water formation is, however, found to be more efficient. Up to 30% of DDT but only less than 0.2% of γ-HCH in seawater are stored in particulate matter.On the time scale studied (1 decade) and on global scale substance transport in the environment is determined by the fast atmospheric circulation. The meridional transport mechanism, for both compounds, is significantly enhanced by multi-hopping. Net meridional transport in the ocean is effective only regionally, mostly by currents along the western boundaries of Africa and the Americas. The total environmental burdens of the substances experience a net northward migration from their source regions, which is more pronounced for DDT than for γ-HCH. Due to the application distribution, however, after 10 years of simulation 21% of the global environmental burden of γ-HCH and 12% of DDT have accumulated in the Arctic.