Multidimensional, multicomponent, subsurface reactive transport in nonuniform velocity fields: code verification using an advective reactive streamtube approach

High performance computing has made possible the development of high resolution, multidimensional, multicomponent reactive transport models that can be used to analyze complex geochemical environments. However, as increasingly complex processes are included in these models, the accuracy of the numerical formulation coupling the nonlinear processes becomes difficult to verify. Analytical solutions are not available for realistically complex problems and benchmark solutions are not generally available for specific problems. We present an advective reactive streamtube (ARS) transport technique that efficiently provides accurate solutions of nonlinear multicomponent reactive transport in nonuniform multidimensional velocity fields. These solutions can be compared with results from Eulerian-based advection-dispersion-reaction models to evaluate the accuracy of the numerical formulation used. The ARS technique includes mixed equilibrium and kinetic complexation and precipitation-dissolution reactions subject to the following assumptions: (1) transport is purely advective (i.e., no explicit diffusion or dispersion), and (2) chemistry is described by a canonical system of reactions that evolves with time and is unaffected by position in space. Results from the ARS technique are compared with results from the massively parallel, multicomponent reactive transport model MCTRACKER on a test problem involving irreversible oxidation of organic carbon and reaction of the oxidation products with two immobile mineral phases, gypsum and calcite, and fifteen aqueous complexes. Truncation error, operator splitting error, and the nonlinear transformation of these errors in the high-resolution reactive transport model are identified for this problem.     

Hydraulic anisotropy characterization of pneumatic-fractured sediments using azimuthal self potential gradient

The pneumatic fracturing technique is used to enhance the permeability and porosity of tight unconsolidated soils (e.g. clays), thereby improving the effectiveness of remediation treatments. Azimuthal self potential gradient (ASPG) surveys were performed on a compacted, unconsolidated clay block in order to evaluate their potential to delineate contaminant migration pathways in a mechanically-induced fracture network. Azimuthal resistivity (ARS) measurements were also made for comparative purposes. Following similar procedures to those used in the field, compressed kaolinite sediments were pneumatically fractured and the resulting fracture geometry characterized from strike analysis of visible fractures combined with strike data from optical borehole televiewer (BHTV) imaging. We subsequently injected a simulated treatment (electrolyte/dye) into the fractures. Both ASPG and ARS data exhibit anisotropic geoelectric signatures resulting from the fracturing. Self potentials observed during injection of electrolyte are consistent with electrokinetic theory and previous laboratory results on a fracture block model. Visual (polar plot) analysis and linear regression of cross plots show ASPG lobes are correlated with azimuths of high fracture strike density, evidence that the ASPG anisotropy is a proxy measure of hydraulic anisotropy created by the pneumatic fracturing. However, ARS data are uncorrelated with fracture strike maxima and resistivity anisotropy is probably dominated by enhanced surface conduction along azimuths of weak ‘starter paths’ formed from pulverization of the clay and increases in interfacial surface area. We find the magnitude of electrokinetic SP scales with the applied N₂ gas pressure gradient (ΔPN₂) for any particular hydraulically-active fracture set and that the positive lobe of the ASPG anomaly indicates the flow direction within the fracture network. These findings demonstrate the use of ASPG in characterizing the effectiveness of (1) pneumatic fracturing and (2) defining likely flow directions of remedial treatments in unconsolidated sediments and rock.     

Effect of sunn hemp (Crotalaria juncea L.) cutting date and planting density on weed suppression in Georgia,USA

A field study was conducted in 2008 and 2009 at the USDA, ARS, Plant Genetic Resources Conservation Unit in Griffin, GA, to investigate weed suppression by sunn hemp (Crotalaria juncea L). The objectives were to (1) evaluate the effects of apical meristem removal (AMR) at three dates [5, 6, and 7 wks after planting (WAP) on May 14, 2008 and May 21, 2009] and (2) assess the impact of seeding rates (11, 28, and 45 kg ha^?1) on weed biomass reduction. Weed species were identified at 4, 8, and 12 wks after sunn hemp planting. Sunn hemp cutting date had no significant effect on weed suppression in 2008 but significant differences for grass weeds at 4, 8, and 12 WAP and for yellow nutsedge at 8 and 12 WAP did occur when compared to the control in 2009. In comparison to the sunn hemp-free control plot in 2009, all three seeding rates had reduced grass weed dry weights at 4, 8, and 12 WAP. The total mass of yellow nutsedge when grown with sunn hemp was reduced compared to the total mass of yellow nutsedge grown in the weedy check for all seeding rates at 8 and 12 WAP. Lower grass weed biomass was observed by 12 WAP for cutting dates and seeding rates during 2008 and 2009. Sunn hemp cutting date and seeding rate reduced branch numbers in both years. The reduction in sunn hemp seeding rates revealed a decrease in weed populations.     

Determination of Sulfur Oxides in the Flue Gases of the Pulping Processes

The body of information presented in this paper is directed to the operating personnel and process engineers employed in the power and recovery departments of a chemical pulping operation. The proper evaluation of the total analytical and sampling system (TASS), to be used in the determination of sulfur oxides is as important as a proper analytical and recording system (ARS). The presence of other sulfur gaseous compounds and particulates could greatly influence the results of the determination.

The analytical method employed determines sulfur dioxide and trioxide from an aliquot of the trapping solution, 3% hydrogen peroxide and 8 0% isopropyl alcohol respectively. The aliquot is titrated with barium perchlorate in the presence of Thorin indicator. The results of evaluating the method indicated negligible interference from the presence of hydrogen sulfide, mercaptans and nitrogen oxides. A blank correction of 15 parts per million (ppm) is recommended whenever 100 ppm of hydrogen sulfide or more are simultaneously present in the gas stream. Particulaies are shown to interfere either by addition or subtraction. Sulfate particulates that will add to the determination must be removed, but in doing so, care must be exerted to avoid surface-contacting conditions that promote reaction between carbonates and the sulfur oxides. The integrated method of sampling and analysis will permit determinations from a flue gas with sulfur oxides concentrations of 30 ppm and above. The relative standard deviation improves from 10% at 100 ppm SO₂ to 2.6% at 1000 ppm SO₂. In both cases, sulfides were present.     

U.S. Environmental Protection Agency’s Super-fund Innovative Technology Evaluation of Pneumatic Fracturing ExtractionSM

The U.S. Environmental Protection Agency (EPA), in cooperation with Accutech Remedial Systems (ARS) and the New Jersey Institute of Technology (NJIT), performed a field demonstration of Pneumatic Fracturing Extraction (PFE)^SM for the removal of chlorinated-volatile organic compounds (VOCs) from vadose zones of low permeability. The demonstration was conducted in the fall of 1992 at an industrial park in Somerville, New Jersey, where removal of VOC contamination in shale bedrock was required to comply with New Jersey’s Environmental Cleanup Responsibility Act (ECRA). During the demonstration, airflow and contaminant concentrations were monitored to establish a database against which the developer’s claims about the technology were evaluated. The developer contended that PFE would increase extracted airfkiw rates from the subsurface formation by at least 100 percent and would increase the mass removal rate for the key contaminant, trichloroethene (TCE), by at least 50 percent. Also, during the demonstration hot-gas injection was evaluated. Based on comparisons of four-hour test results before and after fracturing, airflow rates increased more than 600 percent, and TCE mass-removal rates increased about 675 percent. The increase in TCE mass-removal rates appeared to be a result, primarily, of the increased airflow. In addition, the extracted air contained significantly higher concentrations of other VOCs after fracturing. Using data developed in the four-hour postfracture test, the estimated cost for a hypothetical one-year clean-up is $140 per pound of TCE removed, or $140 per ton of soil contaminated with one pound of TCE. Experiments to evaluate the effects of injecting heated air, at 200 to 250°F, into the vadose zone gave inconclusive results.     

Interactive effects of ozone and elevated carbon dioxide on the growth and physiology of black cherry, green ash, and yellow-poplar seedlings

Potted seedlings of black cherry (Prunus serotina Ehrh.) (BC), green ash (Fraxinus pennsylvanica Marsh.) (GA), and yellow-poplar (Liriodendron tulipifera L.) (YP) were exposed to one of the four treatments: (1) charcoal-filtered air (CF) at ambient CO(2) (control); (2) twice ambient O(3) (2 x O(3)); (3) twice ambient CO(2) (650 microl l(-1)) plus CF air (2 x CO(2)); or (4) twice ambient CO(2) (650 microl l(-1)) plus twice ambient O(3) (2 x CO(2) + 2 x O(3)). The treatments were duplicated in eight continuously stirred tank reactors for 10 weeks. Gas exchange was measured during the last 3 weeks of treatment and all seedlings were destructively harvested after 10 weeks. Significant interactive effects of O(3) and CO(2) on the gas exchange of all three species were limited. The effects of elevated CO(2) and O(3), singly and combined, on light-saturated net photosynthesis (A(max)) and stomatal conductance (g(s)) were inconsistent across species. In all three species, elevated O(3) had no effect on g(s). Elevated CO(2) significantly increased A(max) in GA and YP foliage, and decreased g(s) in YP foliage. Maximum carbon exchange rates and quantum efficiencies derived from light-response curves increased, while compensation irradiance and dark respiration decreased in all three species when exposed to 2 x CO(2). Elevated O(3) affected few of these parameters but any change that was observed was opposite to that from exposure to 2 x CO(2)-air. Interactive effects of CO(2) and O(3) on light-response parameters were limited. Carboxylation efficiencies, derived from CO(2)-response curves (A/C(i) curves) decreased only in YP foliage exposed to 2 x CO(2)-air. In general, growth was significantly stimulated by 2 x CO(2) in all three species; though there were few significant growth responses following exposure to 2 x O(3) or the combination of 2 x CO(2) plus 2 x O(3). Results indicate that responses to interacting stressors such as O(3) and CO(2) are species specific.     

Simultaneous conversion of CHClF(2) and CH(3)Br to CH(2)CF(2)

Gas phase reaction of CHClF(2) with CH(3)Br in an alumina tube reactor at 773-1123 K as a function of various input ratios of CH(3)Br to CHClF(2) is presented. The major products detected include C(2)F(4), CH(2)CF(2), and CH(4). Minor products include CH(3)Cl, CHF(3), C(2)H(4), C(2)H(2), CH(2)CF-CF(3), and C(2)H(3)F. The reaction produces a high yield of CH(2)CF(2) (53% based on CHClF(2) feed) at 1123 K and an input molar ratio of CH(3)Br to CHClF(2) of 1.8, suggesting that the reaction potentially can be developed as a process to convert two ozone depleting substances (CHClF(2) and CH(3)Br) to a highly valuable chemical, CH(2)CF(2). The reaction of CHClF(2) with CH(3)Cl and CH(3)I was also investigated under similar reaction conditions, to assist in understanding the reaction chemistry involved in the reaction of CHClF(2) with CH(3)Br.     

Occurrence of alternative flame retardants in indoor dust from New Zealand: indoor sources and human exposure assessment

Due to worldwide restrictions on polybrominated diphenyl ethers (PBDEs), the demand for alternative flame retardants (AFRs), such as organophosphate flame retardants (OPFRs), novel brominated FRs (NBFRs) and hexabromocyclododecanes (HBCDs), has recently increased. Little is known about human exposure to NBFRs and OPFRs and that their levels in dust have been scarcely evaluated worldwide. To increase the knowledge regarding these chemicals, we measured concentrations of five major NBFRs, ten OPFRs and three HBCD isomers in indoor dust from New Zealand homes. Dust samples were taken from living room floors (n=34) and from mattresses of the same houses (n=16). Concentrations (ngg(-1)) of NBFRs were: 1,2-bis(2,4,6-tribromophenoxy)ethane (BTBPE) (<2-175), decabromodiphenyl ethane (DBDPE) (<5-1430), 2-ethylhexyl-2,3,4,5-tetrabromobenzoate (TBB) (<2-2285) and bis(2-ethylhexyl)-3,4,5,6-tetrabromophthalate (TBPH) (<2-640). For OPFRs, concentrations (ngg(-1)) ranged between: tri-ethyl-phosphate (TEP) (<10-235), tri-n-butyl-phosphate (TnBP) (<20-7545), tris-(2-chloroethyl)-phosphate (TCEP) (<20-7605), tris-(1-chloro-2-propyl) phosphate (TCPP) (20-7615), tri-(2-butoxyethyl)-phosphate (TBEP) (50-27325), tris-(2,3-dichloropropyl)-phosphate (TDCPP) (20-16560), tri-phenyl-phosphate (TPhP) (20-35190), and tri-cresyl-phosphate (TCP) (<50-3760). HBCD concentrations fell in the range <2-4100ngg(-1). BTBPE, DBDPE, TBPH, TBEP, and TnBP showed significant positive correlation (p<0.05) between their concentrations in mattresses and the corresponding floor dust (n=16). These data were used to derive a range of plausible exposure scenarios. Although the estimated exposure is well below the corresponding reference doses (RfDs), caution is needed given the likely future increase in use of these FRs and the currently unknown contribution to human exposure by other pathways such as inhalation and diet.     

Seasonal variation in dissolved gaseous mercury and total mercury concentrations in Juam Reservoir, Korea

Dissolved gaseous mercury (DGM) and total mercury (TM) concentrations were measured in Juam Reservoir, Korea. DGM concentrations were higher in spring (64+/-13pgL(-1)) and summer (109+/-15pgL(-1)), and lower in fall (20+/-2pgL(-1)) and winter (23+/-6pgL(-1)). In contrast, TM concentrations were higher in fall (3.2+/-0.1ngL(-1)) and winter (3.3+/-0.1ngL(-1)) than in spring (2.3+/-0.1ngL(-1)) and summer (2.2+/-0.4ngL(-1)). DGM concentrations were correlated with water temperature (p<0.0001), ORP (p<0.0001), UV intensity (UV-A: p=0.008; UV-B: p=0.003), and DOC concentration (p=0.0107). DGM concentrations varied diurnally with UV intensity. The average summer DGM (109+/-15pgL(-1)) and TM (2.2+/-0.4ngL(-1)) concentrations in Juam Reservoir were higher than the averages for North American lakes (DGM=38+/-16pgL(-1); TM=1.0+/-1.2ngL(-1)), but lower than levels reported for Baihua Reservoir in China.     

Pentachlorophenol association with fulvic acids from recycled wastes

The apparent water solubility of pentachlorophenol was measured at pH=6 and at 25 degrees C in pure water, aqueous solutions of three salts (NaCl, KNO(3) and CaCl(2) at 0.010, 0.10 and 1.0M) and in aqueous solutions of three fulvic acids samples extracted from a natural soil (sFA), composted sewage sludge (csFA) and composted livestock's material (lsFA). A solubility enhancement method was developed for the measurement of partition coefficients (K(oc), L/kg organic carbon). Pentachlorophenol associates strongly with the fulvic acid samples and the calculated K(oc) were the following (averages and standard deviations): (sFA) (211+/-22) x 10(2), (csFA) (253+/-26) x 10(2), (lsFA) (235+/-10) x 10(2). For comparison purposes the K(oc) for pyrene were also calculated for the three FA samples and were the following: (sFA) (119+/-10) x 10(2), (csFA) (239+/-21) x 10(2), (lsFA) (92+/-10) x 10(2). The analysis of variance (one-way ANOVA) of the effect of the type of FA sample on the solubilization of pentachlorophenol and pyrene shows that this factor causes significant differences on the aqueous solubilization of these two organic substances.     

Identification of highly brominated analogues of Q1 in marine mammals

Three novel halogenated organic compounds (HOCs) have been identified in the blubber of marine mammals from coastal New England with the molecular formulae C(9)H(3)N(2)Br(6)Cl, C(9)H(3)N(2)Br(7), and C(9)H(4)N(2)Br(5)Cl. They were identified using high and low resolution gas chromatography mass spectrometry (GCMS) and appear to be highly brominated analogues of Q1, a heptachlorinated HOC suspected to be naturally produced. These compounds were found in Atlantic white sided dolphin (Lagenorhynchus acutus), bottlenose dolphin (Tursiops truncatus), common dolphin (Delphinus delphis), Risso's dolphin (Grampus griseus), harbor porpoise (Phocoena phocoena), beluga whale (Delphinapterus leucas), fin whale (Balaenoptera physalus), grey seal (Halichoerus grypus), harp seal (Phoca groenlandica) and a potential food source (Loligo pealei) with concentrations as high as 2.7 microg/g (lipid weight). The regiospecificity of C(9)H(3)N(2)Br(6)Cl is suggestive of a biogenic origin. Debromination of C(9)H(3)N(2)Br(6)Cl may be significant in the formation of C(9)H(4)N(2)Br(5)Cl.     

Photocatalytic oxidation of NOx using composite sheets containing TiO2 and a metal compound

The photocatalytic oxidation of nitrogen oxides (NO(x)) over titanium dioxide (TiO(2)) sheets containing metal compounds (MCs) had been studied. Calcium oxide (CaO), magnesium oxide (MgO), calcium carbonate (CaCO(3)), aluminium oxide (Al(2)O(3)) and ferric oxide (Fe(2)O(3)) were used as MCs. Al(2)O(3) and Fe(2)O(3) added to the TiO(2) sheet did not affect the photooxidation of nitrogen oxides (NO(x)). The CaO sheet treated with TiO(2) sol had the greatest efficiency as a NO(x) remover under UV irradiation. It is believed that CaO has a high adsorptivity for nitrogen dioxide (NO(2)) and nitric acid (HNO(3)). The amount of NO(x) removed by a TiO(2) sheet including MC showed a tendency to increase with increasing pH of the MC suspension, i.e. there is a good correlation between the alkalinity of the MC and the retention of NO(2) and HNO(3).     

Investigation of a new approach to decompose two potent greenhouse gases: photoreduction of SF(6) and SF(5)CF(3) in the presence of acetone

In this paper, we addressed the utilization of photochemical method as an innovative technology for the destruction and removal of two potent greenhouse gases, SF(6) and SF(5)CF(3). The destruction and removal efficiency (DRE) of the process was determined as a function of excitation wavelength, irradiation time, initial ratio of acetone to SF(5)X (X represented F or CF(3)), initial SF(5)X concentration, additive oxygen and water vapor concentration. A complete removal was achieved by a radiation period of 55min and 120min for SF(6)-CH(3)COCH(3) system and SF(5)CF(3)-CH(3)COCH(3) system respectively under 184.9nm irradiation. Extra addition of water vapor can enhance DRE by approximately 6% points in both systems. Further studies with GC/MS and FT-IR proved that no hazardous products such as S(2)F(10), SO(2)F(2), SOF(2), SOF(4) were generated in this process.     

Fluxes of oxidised and reduced nitrogen above a mixed coniferous forest exposed to various nitrogen emission sources

Concentrations of nitrogen gases (NH(3), NO(2), NO, HONO and HNO(3)) and particles (pNH(4) and pNO(3)) were measured over a mixed coniferous forest impacted by high nitrogen loads. Nitrogen dioxide (NO(2)) represented the main nitrogen form, followed by nitric oxide (NO) and ammonia (NH(3)). A combination of gradient method (NH(3) and NO(x)) and resistance modelling techniques (HNO(3), HONO, pNH(4) and pNO(3)) was used to calculate dry deposition of nitrogen compounds. Net flux of NH(3) amounted to -64 ng N m(-2) s(-1) over the measuring period. Net fluxes of NO(x) were upward (8.5 ng N m(-2) s(-1)) with highest emission in the morning. Fluxes of other gases or aerosols substantially contributed to dry deposition. Total nitrogen deposition was estimated at -48 kg N ha(-1) yr(-1) and consisted for almost 80% of NH(x). Comparison of throughfall nitrogen with total deposition suggested substantial uptake of reduced N (+/-15 kg N ha(-1) yr(-1)) within the canopy.     

A two-generational reproductive toxicity study of zinc in rats

A two-generation reproductive toxicity study of zinc chloride (ZnCl(2)) was conducted in rats. F(o) male and female rats were administered 0.00 (control), 7.50 (low), 15.00 (mid) and 30.00 (high) mg/kg/day of ZnCl(2). Selected F(1) male and female rats were exposed to the same doses received by their parents (F(o)). Exposure of F(0) parental rats to ZnCl(2) showed significant reduction in fertility, viability (days 0 and 4), and the body weight of F(1) pups from the high-dose group but caused no effects on litter size, weaning index, and sex ratio. Similarly, the continued exposure of F(1) parental rats to ZnCl(2) also reduced fertility, liter size, viability (day 0), and the body weight of F(2) pups within the high-dose group but caused no effects on weaning index and sex ratio. Exposure of ZnCl(2) to F(0) and F(1) parental males resulted in a significant reduction in their body weights, and the F(0) and F(1) parental females did not show any significant difference in their body weights compared to their control groups. However, the postpartum dam weights of both F(0) and F(1) female rats were significantly reduced compared to their controls. Exposure of ZnCl(2) to F(o) and F(1) generation parental rats did not produce any significant change of their clinical signs as well as their clinical pathology parameters, except the alkaline phosphotase (ALK) level, which showed an upward trend in both sexes of both generations. Exposure of ZnCl(2) to F(0) rats resulted in a reduction of brain, liver, kidney, spleen and seminal vesicles weights of males and in the spleen and uterus of females. Similarly, exposure of F(1) rats to ZnCl(2) also resulted in reduction of brain, liver, kidney, adrenal, spleen, prostate and seminal vesicles weights of males and in spleen and uterus of females. ZnCl(2) exposure resulted in grossly observed gastro-intestianla (GI) tract, lymphoreticular/hematopoietic, and reproductive tract lesions in parental rats in both generations. Reduced body fat was also recorded in F(1) parental rats.     

Biodegradability and toxicity assessment of trans-chlordane photochemical treatment

The removal of trans-chlordane (C(10)H(6)Cl(8)) from aqueous solutions was studied using UV, UV/H(2)O(2), UV/H(2)O(2)/Fe(2+), UV/TiO(2), or UV/TiO(2)/H(2)O(2) treatment using either UV/Vis blue lamps or UVC lamps (254 nm). H(2)O(2), FeSO(4) and TiO(2) were added at 1700, 456, and 2500 mgL(-1), respectively. trans-Chlordane was not significantly removed in non-irradiated controls and in samples irradiated with UV/Vis. It was also not removed in the absence of surfactant Triton X-114 added at 250 mgL(-1). In the presence of the surfactant, trans-chlordane concentration was reduced by 95-100% after 48 h of UVC and UVC/H(2)O(2) treatments and 70-80% after UVC/H(2)O(2)/Fe(2+), UVC/TiO(2) and UVC/H(2)O(2)/TiO(2) treatments. Based on these results, UVC, UVC/H(2)O(2) and UVC/TiO(2) treatments were further investigated. UVC treatment supported the highest pollutant removal (100% in 48 h), dechlorination efficiency (81% in 48 h), and detoxification to Lepidium sativum seed germination and activated sludge respiration although irradiated samples remained toxic to Chlorella vulgaris. Biodegradation of the UVC irradiated samples removed the source of algae toxicity but this could not be clearly attributed to the removal of trans-chlordane photoproducts because the surfactant interfered with the chemical and biological assays. Evidence was found that trans-chlordane was photodegraded through photolysis causing its successive dechlorination. trans-Chlordane removal was well described by a first order kinetic model at a rate of 0.21±0.01h(-1) at the 95% confidence interval.     

Investigation of heavy metal accumulation in Polygonum thunbergii for phytoextraction

In this study, cadmium (II), lead (II), copper (II) and zinc (II) were determined in Polygonum thunbergii and soil from the Mankyung River watershed, Korea. Soil samples contained detectable lead (<17.5 g g(-1)), copper (<8.4 g g(-1)) and zinc (<24.5 g g(-1)), whereas cadmium was undetectable. Whole plants of P. thunbergii contained detectable lead (<320.8 g g(-1)), copper (<863.2 g g(-1)) and zinc (<2427.3 g g(-1)), whereas cadmium was detectable only in the stem (<7.4 g g(-1)) and root (<10.1 g g(-1)). Whole plant concentrations were very different for each metal, particularly in the case of zinc. The mean content of heavy metal in the whole plants increased in the order of cadmium (8.5 g g(-1))相似文献   

Development of the high spatial resolution EMICAT2000 emission model for air pollutants from the north-eastern Iberian Peninsula (Catalonia, Spain)

Ozone (O(3)) pollution episodes take place in Catalonia (NE of the Iberian Peninsula), mainly during summertime. The complex O(3) behaviour could be understood by using a Chemical Transport Model (CTM). Emission inventories provide the spatial and temporal emissions distribution of the O(3) precursors and other pollutants required by this approach. We developed the EMICAT2000 model with high spatial (cells of 1 km(2)) and temporal (1h) resolutions, to estimate the emissions during the year 2000 from Catalonia. Total annual emissions were 107 kt yr(-1) of NO(x), 137 kt yr(-1) of NMVOC, 267 kt yr(-1) of CO, 65 kt yr(-1) of SO(2), 24 kt yr(-1) of TSP and 32,175 kt yr(-1) of equivalent CO(2). Main NO(x) sources are on-road traffic (58%) and industries (38%). Main NMVOC sources are on-road traffic (36%), vegetation (34%) and use of solvents (13%). Speciation was established according to the Carbon Bond IV mechanism. EMICAT2000 generates directly the data files required for the third generation CTM Models-3/CMAQ.