Evaluation of the effect of meteorological data resolution on Lagrangian particle dispersion simulations using the ETEX experiment

This paper presents results from a series of numerical experiments designed to evaluate operational long-range dispersion model simulations, and to investigate the effect of different temporal and spatial resolution of meteorological data from numerical weather prediction models on these simulations. Results of Lagrangian particle dispersion simulations of the first tracer release of the European Tracer Experiment (ETEX) are presented and compared with measured tracer concentrations. The use of analyzed data of higher resolution from the European Center for Medium-Range Weather Forecasts (ECMWF) model produced significantly better agreement between the concentrations predicted with the dispersion model and the ETEX measurements than the use of lower resolution Navy Operational Global Atmospheric Prediction System (NOGAPS) forecast data. Numerical experiments were performed in which the ECMWF model data with lower vertical resolution (4 instead of 7 levels below 500 mb), lower temporal resolution (12 h instead of 6 h intervals), and lower horizontal resolution (2.5° instead of 0.5°) were used. Degrading the horizontal or temporal resolution of the ECMWF data resulted in decreased accuracy of the dispersion simulations. These results indicate that flow features resolved by the numerical weather prediction model data at approximately 45 km horizontal grid spacing and 6 h time intervals, but not resolved at 225 km spacing and 12 h intervals, made an important contribution to the long-range dispersion.     

Determination of dibenzo-p-dioxin (DBD) and dibenzofuran (DBF) in paper pulp mill defoamer additives and in paper pulp samples

An analytical procedure has been developed for quantitatively measuring dibenzo-p-dioxin (DBD) and dibenzofuran (DBF) present at part-per-billion levels in paper mill defoamers and paper pulp using gas chromatography coupled with either low resolution mass spectrometry (GC-LRMS) or high resolution mass spectrometry (GC-HRMS).     

Influence of model grid resolution on NO2 vertical column densities over East Asia

NO₂ vertical column densities (VCDs) over East Asia in June and December 2007 were simulated by the Community Multi-scale Air Quality (CMAQ) version 4.7.1 using an updated and more elaborate version of the Regional Emission Inventory in Asia (REAS) version 2. The modeling system could reasonably capture observed spatiotemporal changes of NO₂ VCDs by satellite sensors, the Global Ozone Monitoring Experiment-2 (GOME-2), the Scanning Imaging Absorption Spectrometer for Atmospheric Cartography (SCIAMACHY), and the Ozone Monitoring Instrument (OMI), even at the coarsest horizontal resolution of 80 km. The CMAQ simulations were performed in a sequence of three horizontal resolutions (80 km, 40 km, and 20 km) for June and December 2007 to investigate the influence of changes of horizontal resolution on the obtained NO₂ VCDs. CMAQ-simulated NO₂ VCDs generally increased with improvements in resolution from 80 km to 40 km and then to 20 km. Increases in the CMAQ-simulated NO₂ VCDs were greater for the change from 80 km to 40 km than for those from 40 km and 20 km, in which the increases of NO₂ VCDs due to the improvement of horizontal resolution were approached convergence at the horizontal resolution of approximately 20 km. Conversely, no clear convergences in NO₂ VCDs changes were found at near Tokyo and over the East China Sea. The biases of the NO₂ VCDs simulated at a resolution of 20 km against the satellite retrievals were ～36% near Beijing (CHN1) and ～78% near Shanghai (CHN2) in summer; these errors were found to be comparable to the horizontal resolution-dependent errors, which were 18–25% at CHN1 and 44–58% at CHN2 from 80 km to 40 km. Conversely, the influence of changes of horizontal resolution in winter was relatively less compared to that in summer.

Implications: NO₂ VCDs over East Asia in June and December 2007 were simulated using CMAQ version 4.7.1 and REAS version 2. The modeling system could reasonably capture observed spatiotemporal changes of NO₂ VCDs by satellite sensors. The CMAQ simulations were performed in a sequence of three horizontal resolutions, 80, 40, and 20 km, to investigate the influence of changes of horizontal resolution on the obtained NO₂ VCDs. The results suggested that the influence of changes of horizontal resolution was larger in summer compared to that in winter. The magnitude of the influence was comparable to the biases of the NO₂ VCDs simulated at a resolution of 20 km against the satellite retrievals.     

Modeling atmospheric mercury deposition in the vicinity of power plants

Two mathematical models of the atmospheric fate and transport of mercury (Hg), an Eulerian grid-based model and a Gaussian plume model, are used to calculate the atmospheric deposition of Hg in the vicinity (i.e., within 50 km) of five coal-fired power plants. The former is applied using two different horizontal resolutions: coarse (84 km) and fine (16.7 km). More than 96% of the power plant Hg emissions are calculated with the plume model to be transported beyond 50 km from the plants. The grid-based model predicts a lower fraction to be transported beyond 50 km: >91% with a coarse resolution and >95% with a fine resolution. The contribution of the power plant emissions to total Hg deposition within a radius of 50 km from the plants is calculated to be <8% with the plume model, <14% with the Eulerian model with a coarse resolution, and <10% with the Eulerian model with a fine resolution. The Eulerian grid-based model predicts greater local impacts than the plume model because of artificially enhanced vertical dispersion; the former predicts about twice as much Hg deposition as the latter when the area considered is commensurate with the resolution of the grid-based model. If one compares the local impacts for an area that is significantly less than the grid-based model resolution, then the grid-based model may predict lower local deposition than the plume model, because two compensating errors affect the results obtained with the grid-based model: initial dilution of the power plant emissions within one or more grid cells and enhanced vertical mixing to the ground.     

Method verification for determination of tetrachlorodibenzodioxin in soil

Methods for the determination of tetrachlorodibenzodioxin (TCDD) in surface and subsurface soil have been verified by spiking experiments conducted at a nominal level of twice the limit of detection. The methods involve Soxhlet extraction, silica and alumina clean-up and analysis by high resolution gas chromatography/ high resolution mass spectrometry. Mass Spectrometry done under profile data control is described. The precision (accuracy) found for surface and subsurface soils were 10.8% (85%) and 10.8% (95%) respectively. Application of the surface soil procedure is illustrated by analysis of soils from the Northeast United States.     

Modeling Atmospheric Mercury Deposition in the Vicinity of Power Plants

Abstract

Two mathematical models of the atmospheric fate and transport of mercury (Hg), an Eulerian grid–based model and a Gaussian plume model, are used to calculate the atmospheric deposition of Hg in the vicinity (i.e., within 50 km) of five coal–fired power plants. The former is applied using two different horizontal resolutions: coarse (84 km) and fine (16.7 km). More than 96% of the power plant Hg emissions are calculated with the plume model to be transported beyond 50 km from the plants. The grid–based model predicts a lower fraction to be transported beyond 50 km: >91% with a coarse resolution and >95% with a fine resolution. The contribution of the power plant emissions to total Hg deposition within a radius of 50 km from the plants is calculated to be <8% with the plume model, <14% with the Eulerian model with a coarse resolution, and <10% with the Eulerian model with a fine resolution. The Eulerian grid–based model predicts greater local impacts than the plume model because of artificially enhanced vertical dispersion; the former predicts about twice as much Hg deposition as the latter when the area considered is commensurate with the resolution of the grid–based model. If one compares the local impacts for an area that is significantly less than the grid–based model resolution, then the grid–based model may predict lower local deposition than the plume model, because two compensating errors affect the results obtained with the grid–based model: initial dilution of the power plant emissions within one or more grid cells and enhanced vertical mixing to the ground.     

Experimental verification for real-time environmental emergency response system: WSPEEDI by European tracer experiment

Japan Atomic Energy Research Institute has developed an emergency response system WSPEEDI to forecast long-range atmospheric dispersions of radionuclides discharged into the atmosphere. The latest version of WSPEEDI consists of an atmospheric dynamic model MM5 for calculating meteorological fields and a particle random-walk model for atmospheric dispersion. The performance of WSPEEDI was evaluated by data obtained from a field tracer experiment over Europe (ETEX) in this paper. The model validation was done with respect to the following points: (1) the dependence of model accuracy on the temporal and spatial resolutions of the meteorological fields and (2) the superiority of an atmospheric dynamic model over a mass-consistent wind model. Regarding (1), it was shown that the calculation accuracy of the new version with high temporal resolution was improved, especially at the edge of the plume. Moreover, although the increase in horizontal spatial resolution of the old version had no substantial effect on the model performance, increase in horizontal resolution of the new version contributed to the significant improvement of the calculation accuracy. These results showed that the dynamically calculated meteorological field with the spatial resolution of the meso-β–γ scale greatly improved calculation accuracy.     

Estimating PM2.5 over southern Sweden using space-borne optical measurements

In the present study Bremen aerosol retrieval (BAER) columnar aerosol optical thickness (AOT) data, according to moderate resolution imaging spectroradiometer (MODIS) and medium resolution imaging sensor (MERIS) level 1 calibrated satellite data, have been compared with AOT data obtained with the MODIS and MERIS retrieval algorithms (NASA and ESA, respectively) and by AErosol RObotic NETwork (AERONET). Relatively good agreement is found between these different instruments and algorithms. The R² and relative RMSD were 0.86 and 31% for MODIS when comparing with AERONET and 0.92 and 21% for MERIS. The aerosols investigated were influenced by low relative humidity. During this period, a relatively large range of aerosol loadings were detected; from continental background aerosol to particles emitted from agricultural fires. In this study, empirical relationships between BAER columnar AOT and ground-measured PM_2.5 have been estimated. Linear relationships, with R² values of 0.58 and 0.59, were obtained according to MERIS and MODIS data, respectively. The slopes of the regression of AOT versus PM_2.5 are lower than previous studies, but this could easily be explained by considering the effect of hygroscopic growth. The present AOT–PM_2.5 relationship has been applied on MERIS full resolution data over the urban area of Stockholm and the results have been compared with particle mass concentrations from dispersion model calculations. It seems that the satellite data with the 300 m resolution can resolve the expected increased concentrations due to emissions along the main highways close to the city. Significant uncertainties in the spatial distribution of PM_2.5 across land/ocean boundaries were particularly evident when analyzing the high resolution satellite data.     

Polychlorinated dibenzo-p-dioxins/furans and polychlorinated biphenyls in fresh fishes from Qiantangjiang River, China

Muscles of two species of fish collected in Qiangtangjiang River were analyzed for polychlorinated dibenzo-p-dioxins/dibenzofurans (PCDD/Fs) and polychlorinated biphenyls (PCBs) by high resolution gas chromatography/high resolution mass spectrometry (HRGC/HRMS). Average concentrations of total PCDD/Fs in Crucian carp and White Amur bream were 1.14-7.88 pg g(-1)ww and 1.02-8.18 pg g(-1)ww respectively. Corresponding values for PCBs were 469.36-10972 pg g(-1)ww and 364.79-4948.0 pg g(-1)ww. Average total WHO-TEQ of Crucian carp and White Amur bream from five areas ranged from 0.25 to 2.33 pg g(-1)ww and 0.30 to 1.70 pg g(-1)ww. The contamination level was compared with other studies and risk assessment was discussed.     

Levels of dioxins and dioxin-like PCBs in food samples on the Greek market

Food intake is the main source of exposure to dioxin-like compounds for humans. The results of a surveillance programme on polychlorinated dibenzo-p-dioxins (PCDDs), dibenzofurans (PCDFs) and co-planar polychlorinated biphenyls (co-PCBs) in 77 food samples from the Greek market and producers are presented. The study included the analyses of milk and dairy products, meat and meat products, fish, vegetable oil, eggs, fruit, vegetable and rice collected between August and December 2002. After extraction, extracts were cleaned up on a series of carbon column chromatography, silica gel, alumina chromatography, and then analysed by high resolution gas chromatography/high resolution mass spectrometry. All samples had a dioxin content far below the EC Regulation (2375/2001/EC) limits.     

Chlorinated fluorenes and alkylfluorenes in bleached kraft pulp and pulp mill discharges

The occurrence and levels of polychlorinated fluorenes (PCFLs), methylfluorenes (C₁-PCFLs) and dimethylfluorenes (C₂-PCFLs) were studied by high resolution gas chromatography/high resolution mass spectrometry (HRGC/HRMS) in hardwood and softwood pulp, effluent and biosludge from a bleached kraft mill. Semiquantitative concentrations of these compounds were estimated using selected ion monitoring and a mass spectral resolution of 10,000. The identification of PCFLs and C₁-PCFLs was confirmed by the chlorination mixtures of fluorene and 2-methylfluorene.

Several congeners were found for each compound groups studied and isomer patterns were similar in bleached pulp and biosludge. The approximated sum levels of di- and trichloro isomers of PCFLs, C₁-PCFLs and C₂-PCFLs in the biosludge were about 2, 5 and 4 ng/g dry weight (dw), respectively. The levels of the corresponding compounds in the softwood pulp were 0.5, 0.8 and 0.2 ng/g dw, respectively. Of the studied compounds, only C₁-diCFLs (0.01 ng/g dw) were detected in the hardwood pulp. In the effluent, PCFLs were non-detectable, but traces of C₁-diCFLs (0.5 ng/L) and C₂-pentaCFLs (0.06 ng/L) were found. These were observed mainly as dissolved in the filtrate and not as particle-bound.     

Dioxin like compounds from municipal waste incinerator emissions: assessment of the presence of polychlorinated naphthalenes

Polychloronaphthalenes (PCN) were identified and quantified in emission samples collected from five different municipal waste incinerators (MWI). Polychlorodibenzo-p-dioxins (PCDD) and polychlorodibenzofurans (PCDF) were also determined to find a possible relationship between these classes of organochlorinated compounds. The analyses of PCDD/PCDF and PCN were carried out by high resolution gas chromatography coupled with high resolution mass spectrometry using a positive electron ionization source and operating in the selected ion monitoring analyzer mode (HRGC-HRMS/EI(+)-SIM). The total levels of PCN varied from 1.08 up to 21.36 ng/Nm3 (mono- to octachlorinated) and 0.33 to 5.72 ng/Nm3 (tetra- to octachlorinated), whereas the levels of PCDD/PCDF ranged between 1.14 and 276.26 ng/Nm3 (0.01 and 5 ng I-TEQ/Nm3), depending on the type of the MWI. These findings do not corroborate a PCN and PCDD/PCDF correlation.     

The determination of toxaphene in environmental samples by negative ion electron capture high resolution mass spectrometry

A method is described for the analysis of toxaphene in environmental samples by high resolution gas chromatography high resolution mass spectrometry in selected ion electron capture negative ion mode. Toxaphene is detected by monitoring [M-Cl]- fragments ions of chlorobornanes (Cl6-Cl10) and chlorobornenes (Cl6-Cl7) using formal ion ratio criteria. The method shows high sensitivity and excellent specificity for the detection of chlorobornanes with negligible interference from of a wide range of common environmental organochlorines including chlordanes, PCBs and PCB-oxygen adducts. The analytical scheme is applicable to sediment and biological tissue samples and typically provide sample detection limits of less than 0.2 ng/g technical toxaphene for a 10 g sample.     

Compilation of a database on the composition of anthropogenic VOC emissions for atmospheric modeling in Europe

To analyse and generate air pollution control strategies and policies, e.g. efficient abatement strategies or action plans that lead to a fulfilment of air quality aims, atmospheric dispersion models (CTMs) have to be used. These models include a chemical model, where the numerous volatile organic compounds (VOCs) species are lumped together in classes. On the other hand, emission inventories usually report only total non-methane VOC (NMVOC), but not a subdivision into these classes. Thus, VOC species profiles are needed that resolve total NMVOC emission data. The objective of this publication is to present the results of a compilation of VOC species profiles that dissolve total VOC into single-species profiles for all relevant anthropogenic emission source categories and the European situation. As in atmospheric dispersion models usually modules for generating biogenic emissions are directly included, only anthropogenic emissions are addressed. VOC species profiles for 87 emission source categories have been developed. The underlying data base can be used to generate the data for all chemical mechanisms. The species profiles have been generated using recent measurements and studies on VOC species resolution and thus represent the current state of knowledge in this area. The results can be used to create input data for atmospheric dispersion models in Europe.The profiles, especially those for solvent use, still show large uncertainties. There is still an enormous need for further measurements to achieve an improved species resolution. In addition, the solvent use directive and the DECOPAINT directive of the European Commission will result in a change of the composition of paints; more water-based and high-solid paints will be used; thus the species resolution will change drastically in the next years. Of course, the species resolution for combustion and production processes also requires further improvement.     

Isomeric analysis of BTEXs in the atmosphere using beta-cyclodextrin capillary chromatography coupled with thermal desorption and mass spectrometry

An analytical method capable of determining trace levels of BTEX-aromatics (benzene, toluene, ethylbenzene, m-, p- and o-xylenes) in the atmosphere with as high resolution as possible has been developed. The method is based on the preconcentration of air samples using a multibed tube (Carbopack C, Carbograph 1) at ambient temperature, followed by thermal desorption, and analysis of aromatic species by a beta-cyclodextrin capillary chromatography coupled with mass spectrometry. The resolution achieved was sufficient for individual separation of BTEXs as well as m- and p-xylenes. The BTEX-ratios have been determined in an air tunnel and in on-road, suburban and rural forest atmosphere. The ethylbenzene/m-xylene ratios could provide a deep insight into anthropogenic related NMHC patterns at different locations and under different meteorological conditions and may reflect photochemical processes in the best way.     

Evaluation of a long-range lagrangian dispersion model with ETEX

Performance of a Lagrangian dispersion model was examined in connection with its dependency on the boundary layer modelling and the input data resolution. The European Tracer Experiment (ETEX) data were used as reference. According to the sensitivity analysis of the model performance, the long-range dispersion model with the sparse input data was not noticeably different from that with the finer resolution data. The assumption of the prescribed constant mixing depth did not largely degrade the prediction results as compared with the simulation results with the temporally changing boundary layer. It is, therefore, concluded that the model is practical, considering the limited input data in the operational mode. However, it was also pointed out that the parameterization for the horizontal and vertical diffusion processes used in the present model enhanced the growth of plume. The improvement of input data resolution in time and space caused further dispersion of tracer deterministically. These resulted in the underestimation of the maximum concentration and the unfocussed concentration distribution map although the mean concentration was predicted fairly well.     

Concentration of polychlorinated dibenzo-p-dioxins, polychlorinated dibenzofurans and dioxin-like PCBs in human adipose tissue from Turkish men

There is no previous report from Turkey on chemically determined polychlorinated dibenzo-para-dioxins (PCDDs), polychlorinated dibenzofurans (PCDFs), and dioxin-like polychlorinated biphenyls (PCBs) in human tissues expressed as World Health Organization (WHO) toxic equivalents (TEQs). The objective of this study was to determine the occurrence of PCDDs/Fs, and dioxin-like PCBs in the general adult Turkish population. For this reason we measured adipose tissue concentrations of PCDDs/Fs and dioxin-like PCBs in 23 Turkish men living in Ankara,Turkey in 2004. PCDD/F concentrations ranged between 3.2 and 19.7 pg WHO-TEQ/g fat (5.34 and 42.7 WHO-TEQ/g fat, respectively, including dioxin-like PCBs). The mean concentrations of WHO(PCDD/F)-TEQ and WHO(PCB)-TEQ were 9.2 and 6.67 pg/g on a lipid basis , respectively. Samples were analyzed for PCDD/F and twelve dioxin-like PCB congeners using high resolution gas chromatography/high resolution mass spectrometry (HRGC/HRMS). This study is very important since it is the first report on PCDDs/Fs and dioxin-like PCB contamination in human adipose tissue from Turkey.     

Analysis of polychlorinated dibenzo-p-dioxins and furans in various aqueous samples in Taiwan

Twenty-one fresh water samples were taken from reservoirs, wells and drinking water treatment plants in Taiwan for measurement of polychlorinated dibenzo-p-dioxins and dibenzofurans (PCDD/F) concentrations. Due to low concentration of PCDD/Fs, water samples (400-1000 L) were collected by using the on-site large volume pre-concentration system and analyzed by high resolution gas chromatography (HRGC)/high resolution mass spectrometer (HRMS). Different characteristics of dioxin profiles in surface water and groundwater are observed. The PCDD/F concentrations in fresh water samples varied from 0.001 to 0.265 pg WHO₉₈-TEQ/L which are lower than the standard (12 pg WHO₉₈-TEQ/L) set by the Taiwan Environmental Protection Administration (TEPA) for drinking water. Only a small fraction (<5%) of the total PCDD/F concentration in fresh water samples was water-soluble. As a result of their hydrophobicity, PCDD/Fs are effectively removed along with suspended solid via the water treatment process (average 98% removal efficiency).     

Assessment of human exposure to PCDDs,PCDFs and Co-PCBs using hair as a human pollution indicator sample I: Development of analytical method for human hair and evaluation for exposure assessment

Dioxins including polychlorinated dibenzo-p-dioxins (PCDDs), polychlorinated dibenzofurans (PCDFs) and co-planar polychlorinated biphenyls (Co-PCBs) are highly toxic. Even at environmental pollution levels, they cause hormonal damage in women, and they have been shown to induce immunosuppression and genital function damage in humans. In this study, a new method using isotope dilution was established to detect PCDDs, PCDFs and Co-PCBs in human hair. This method, comprised of washing and cutting of hair, alkaline decomposition, hexane extraction, multilayer silica gel column chromatography, high performance liquid chromatography with a porous graphite carbon column and analysis by high resolution gas chromatography/high resolution mass spectrometry, enabled us to analyze PCDDs, PCDFs and Co-PCBs at trace levels of less than pg/g with good reproducibility. In addition, there was a correlation between some isomers in human hair and blood collected from identical donors. Human hair analysis is useful to evaluate human risk assessment including that due to environmental pollution.     

Evaluating the ability of a numerical weather prediction model to forecast tracer concentrations during ETEX 2

In this paper the meteorological processes responsible for transporting tracer during the second ETEX (European Tracer EXperiment) release are determined using the UK Met Office Unified Model (UM). The UM predicted distribution of tracer is also compared with observations from the ETEX campaign. The dominant meteorological process is a warm conveyor belt which transports large amounts of tracer away from the surface up to a height of 4 km over a 36 h period. Convection is also an important process, transporting tracer to heights of up to 8 km. Potential sources of error when using an operational numerical weather prediction model to forecast air quality are also investigated. These potential sources of error include model dynamics, model resolution and model physics. In the UM a semi-Lagrangian monotonic advection scheme is used with cubic polynomial interpolation. This can predict unrealistic negative values of tracer which are subsequently set to zero, and hence results in an overprediction of tracer concentrations. In order to conserve mass in the UM tracer simulations it was necessary to include a flux corrected transport method. Model resolution can also affect the accuracy of predicted tracer distributions. Low resolution simulations (50 km grid length) were unable to resolve a change in wind direction observed during ETEX 2, this led to an error in the transport direction and hence an error in tracer distribution. High resolution simulations (12 km grid length) captured the change in wind direction and hence produced a tracer distribution that compared better with the observations. The representation of convective mixing was found to have a large effect on the vertical transport of tracer. Turning off the convective mixing parameterisation in the UM significantly reduced the vertical transport of tracer. Finally, air quality forecasts were found to be sensitive to the timing of synoptic scale features. Errors in the position of the cold front relative to the tracer release location of only 1 h resulted in changes in the predicted tracer concentrations that were of the same order of magnitude as the absolute tracer concentrations.