United States Environmental Protection Agency Environmental Sciences Research Laboratory Research Triangle Park NC 27711 Research and Development EPA-600/S3-84-001 Jan. 1984 &EPA Project Summary Natural Hydrocarbon Emission Rate Measurements from Selected Forest Sites Brian K. Lamb, Hal H. Westberg, Timothy Quarles, and Donald L. Flyckt Studies of biogenic hydrocarbon emis- sions were conducted in a hardwood forest in Pennsylvania during 1979 and in a coniferous forest in Washington during 1980. The principal objective of the studies was to compare a branch enclosure method with a micrometeor- ological gradient technique for measur- ing biogenic hydrocarbon fluxes for forested areas. A second important objective for the Pennsylvania work was to develop a regional natural hydro- carbon emission inventory for use in the Northeast Regional Oxidant Study. Isoprene emission fluxes determined by the gradient profile procedure in the deciduous forest agreed reasonably well with those measured using the enclosure technique. The isoprene flux from the gradient profile data was 10% higher than the enclosure flux at 30°C, but was approximately three times greater than the enclosure flux at 20° C. The differences at the low temperatures possibly were caused by the lack of profile data at the lower temperatures. The combined enclosure and gradient profile data were correlated with ambi- ent temperature to the same degree as the correlation of each data set alone. In the Washington study, the alpha-pinene flux as measured by the gradient profile method ranged from 76 to 1,320 //g/mz-hr whereas, the range deter- mined using the enclosure method was 9 to 700 /ug/m2-hr. The mean fluxes from the two methods were within the estimated limits of uncertainty. Alpha- pinene fluxes determined with the gradi- ent profile method increased exponen- tially with increasing relative humidity. Emission fluxes calculated from the branch enclosure samples were not correlated with ambient relative humid- ity, but were strongly correlated with temperature when wet and dry branches were considered separately. Total biogenic hydrocarbon emissions from the state of Pennsylvania were calculated to be 3400 tons/day during August. Approximately 75% of these emissions were from forested lands, and the remainder were from agricul- tural crops, primarily corn. The forest emissions were approximately evenly divided among isoprene-emitting and non-isoprene-emitting hardwoods and softwood trees. This Project Summary was developed by £PA's Environmental Sciences Re- search Laboratory, Research Triangle Park, NC, to announce key findings of the research project that is fully docu- mented in a separate report of the same title (see Project Report ordering in- formation at back). Introduction The significance of naturally emitted hydrocarbons in rural atmospheric chem- istry remains uncertain. This uncertainty arises in part from a lack of information concerning the magnitude of sources and the distribution of hydrocarbon species in the atmosphere. Current estimates of natural hydrocarbon emission rates have been obtained using a vegetation enclo- sure technique, a micrometeorological gradient profile method, and an energy balance/Bowen ratio approach. Confi- dence in these methods has been limited by uncertainties about the effects of enclosing vegetation and in measure- ments of small vertical gradients of ------- temperature, wind speed, and hydro- carbon concentrations above a forest. These experimental problems have pre- vented widespread agreement upon actu- al hydrocarbon emission rates. In this report, field measurements of biogenic hydrocarbon emission rates are reported for the predominant species of trees and some agricultural crops in the northeastern U.S. and for a coniferous forest in the Pacific Northwest. These data were obtained using a vegetation enclosure method and, for the forests, a micrometeorological gradient profile tech- nique. The major objectives of the investi- gation were (1) to obtain natural hydro- carbon emission rates from vegetation in the northeastern U.S. as a basis for modeling natural hydrocarbon fluxes in the Northeast Regional Oxidant Study (NEROS), (2) to compare emission fluxes determined by a vegetation enclosure technique and a micrometeorological gradient profile method, and (3) to com- pare the chemical composition of essen- tial plant oils with the composition of gas emissions from selected vegetation. Procedure In the forest studies, the flux of a particular hydrocarbon was calculated from surface layer theory, based upon measurements of vertical wind speed, temperature, and hydrocarbon concentra- tion profiles collected along a tower rising above the forest canopy. For the hard- wood forest, 30-min average air samples were collected at six levels on the tower using Teflon sampling pumps and 6-L Tedlar bags. For the coniferous forest, 30- min average air samples were collected at five levels on a tower in stainless steel tubes packed with Tenax-GC adsorbent. Vegetation enclosure samples in each of the study areas were obtained using a branch enclosure method. This involved enclosing a branch in a 100-LTedlar bag, collecting a background sample from the bag, filling the bag with hydrocarbon-free air, and, after a measured length of time, collecting a second sample of the bag air.. When used to develop an area emission inventory, the individual branch emission rates (yug/g-hr) were multiplied by an appropriate biomass factor (g/m2) to give the area emission flux (//g/m2-hr). All hydrocarbon samples were analyzed with Perkin-Elmer and Hewlett-Packard gas chromatographs equipped with subambi- ent temperature programming capabil- ities. The hardwood forest was located in gently sloping terrain approximately 10 km northeast of York, Pennsylvania. A vegetation inventory of the forest indicat- ed that the woodland was a second growth oak and chestnut forest. The biomass factor for isoprene-emitting spe- cies was determined to equal 379 g/m2 from the vegetation inventory and an empirical relationship between the diam- eter of the tree at breast height (DBH) and the biomass. The biomass factor for non- isoprene emitting species was 206 g/m2. The average canopy height was 20 m. The coniferous forest was located 36 km southeast of Seattle, Washington, in the A. E. Thompson Research Forest (University of Washington). The forest consisted almost entirely of Douglas fir, with red alder, western hemlock, and western red cedar present in small amounts. The age of the stand was 52 years and the average canopy height was 31 m. The biomass factor was estimated to equal 830 g/m2 from a relationship between biomass and the age of the stand developed specifically for the Thomp- son forest. Results Northeastern Biogenic Hydrocarbon Emissions Results from the Pennsylvania study indicated that volatile hydrocarbon emis- sions can be divided into four categories based on similarities m composition and emission source. The first two categories are represented by hardwood trees that (1) emit isoprene and (2) do not emit isoprene. The two remaining groups are (3) monoterpene emitting softwood trees and (4) agricultural crops. Included in the group of isoprene emitting hardwood trees measured were oak, black locust, and sycamore. Isoprene, which is only emitted during daylight hours, accounted for 78% of the total volatile hydrocarbon emissions from this group. As shown in Figure 1, the emission rates were directly related to temperature in an exponential manner, with 90% of the variation in emission rates associated with temperature variations. The regres- sion relationship between temperature and emission rate predicts a total hydro- carbon emission rate of 6.1 /ug/g-hr at 25°C. The non-isoprene-emitting hardwood trees measured included black gum, sassafras, tulip tree, red maple, dogwood, red hickory, black cherry, beech, silver maple, and birch. The hydrocarbon emis- sions from these species were also related to temperature in a exponential manner, with 91% of the variation in emission rates associated with temper- . ature fluctuations. The regression rela- tionship between temperature and emis- sion rates shown in Figure 1 predicts an emission rate of 3.4 /jg/g-hr at 25°C. Major identified components of the vola- tile emissions were the monoterpenes a- pinene, sabinene, /3-pinene, myrcene, camphene, /3-phellandrene, A3-carene, and linalool. The monoterpene-emitting softwood species measured included eastern white pine, Virginia pine, eastern hemlock, and pitch pine. The volatile emissions from this group consisted of a-pinene, sabi- nene, camphene, /3-pinene, myrcene, fl- phellandrene, limonene, and A3-carene. These compounds comprised 50-95% of the volatile emissions. Correlation be- tween emission rate and temperature was not observed for this group; thus, the release of the monoterpenes must be controlled by some temperature-inde- pendent process. The agricultural crop category meas- ured included tobacco, corn, clover, alfal- fa, and pasture. This group is poorly defined at present and may require further divisions. The members emitted compounds eluting near the monoter- penes and very small amounts of iso- prene. There appeared to be no relation- ship between emisison rate and temper- ature. The emission rates established in this study can be combined with vegetation surveys and biomass relationships to provide estimates of the magnitude of biogenic emissions. Total biogenic hydro- carbon emissions from the state of Penn- sylvania were calculated to be 3400 tons/day during August. Approximately 74% of these emissions were from forest- ed lands. The balance was estimated to result primarily from corn (25%) and other agricultural crops and pasture. The forest emissions were approximately equally divided among isoprene-emitting hard- woods (890 tons/day), non-isoprene-emit- ting hardwoods (710 tons/day), and conifer species (910 tons/day). In order to gain insight into the pro- cesses that control hydrocarbon emis- sions from vegetation, extracts of essen- tial oils from leaf and branch samples were obtained at the Pennsylvania study site. For the monoterpene-emitting soft- woods, the composition of the oils in the needles and branches match the range of compounds identified in the emission samples quite closely. This agreement suggests that volatilization of oils is the major source of emissions from these * species. For isoprene-emitting hard- ------- 100.0 60.0 20.0 10.0 6.0 o ^ 2.0 1.0 - 0.6 ~ 0.4 10 15 30 35 Figure 1. 20 25 Temperature °C Total nonmethane hydrocarbon emission rate as a function of temperature for isoprene-emittinghardwoods (PA A 1-—)andnon-isoprene-emitting hardwoods (PA • 2—/ woods, the extraction results showed little similarity to the emission samples and no isoprene was measured in the extraction samples. Comparison of Vegetation Enclosure and Micrometeorological Gradient Techniques In Pennsylvania, isoprene fluxes meas- ured by both techniques were observed to vary exponentially with temperature. The isoprene flux estimated from the enclo- sure data at 20°C ambient was 890 /j- g/m2-hr and the flux calculated from the gradient profile data was 2510 Ai-g/m2- hr. At 30°C ambient, the difference between estimates obtained using the two methods was approximately 10% (7300/^-g/m2-hrfrom the enclosure data and 8000 /u-g/m2-hr from the gradient profile method). The differences at 30°C were within the estimated experimental uncertainty of the two methods of meas- urement, but the results at the lower temperature were not within the error limits. However, very few gradient profile experiments were conducted at low tem- peratures. Furthermore, as shown in Figure 2, the combined data from the enclosure and gradient profile measure- ments were correlated with ambient temperature to the same degree as either data set alone. Thus, the data obtained in Pennsylvania indicate that the two inde- pendent methodsfor measuring isoprene flux yield results that are in reasonable agreement. In Washington, alpha-pmene flux meas- ured via the gradient profile technique increased exponentially with increasing relative humidity at relative humidities greater than approximately 60%. The calculated flux was also positively corre- lated with the ratio of relative humidity to wind speed as shown in Figure 3. Alpha- pinene fluxes were only weakly corre- lated with ambient temperature. The mean flux of alpha-pinene measured in seven gradient profile experiments was 440 /jg/m2-hr (median = 230 A/g/m2-hr), with a range of 76 to 1320 Aig/m2-hr. The mean alpha-pinene flux estimated from 13 enclosure samples was 150 mg/m2-hr (median = 46 //g/m2-hr) with a range of 9 to 700 //g/m2-hr. No relation between ambient relative humidity and alpha-pinene flux was observed in the enclosure data. However, when the sam- ples were grouped according to wet and dry branches, a distinct exponential rela- tionship with increasing ambient temper- ature was evident, as shown in Figure 4. No single environmental parameter was correlated with the fluxes determined by both methods. As a result, it was not possible to compare directly fluxes pre- dicted for a common reference point, as was done for isoprene as a function of temperature. The mean flux of alpha- pinene for the gradient profile experi- ments and the mean flux predicted by the enclosure samples were within experi- mental uncertainties estimated for the two methods. The wide range of fluxes observed was not a result of random experimental errors. Rather, changes in environmental conditions such as temper- ature and humidity were observed to affect the emission rates dramatically. In view of the strong dependence of flux upon humidity, the lowfluxesdetermined with the enclosure method suggested that humidity inside the bag was generally lower than ambient. The error limits estimated for the two techniques were based on standard error analysis procedures. Data from typical experiments were combined with esti- mated or measured uncertainties in the independent variables to calculate the overall error and degree of reproducibility in each flux method. In Pennsylvania, the error in the gradient profile method was estimated to equal ± 35%; in Washington, the error was approximately ± 55%. The major sources of error were uncertainties in the zero plane displacement height and the hydrocarbon concentration gradient. The gradient profile method is estimated to be reproducible to within less than ± 17% for test conducted at a specific site. The overall error estimated for a typical enclosure sample was ± 40%, and the technique was estimated to be reproduc- ible to within ± 20% for samples collected in a particular forest. Summary and Conclusions For a typical northeastern deciduous forest, isoprene emissions accounted for approximately 78% of the total non- ------- 10.000 1.000 0) 100 Figure 2. Gradient r = 0.65 Gradient + Enclosure r = 0.63 Enclosure (T*mKJ r = 0.60 Enclosure (T r = 0.85 10 20 Temperature °C 30 Isoprene flux as a function of ambient temperature measured at 6m above the canopy by the gradient profile method (line 1). the gradient profile and enclosure method (line 2). and only the enclosure method (line 3) for Pennsylvania. Line 4 is from the Pennsylvania enclosure data as a function of enclosure temperature. ed here show that the enclosure method yields emission fluxes quite similar to those obtained with a totally independent technique. The enclosure system is portable, easily operated by one person and not limited to idealized sites. However, care is required to use the enclosure method. Enclosure conditions must be closely monitored in order to relate the data to ambient conditions. Biomass factors must be developed from site-specific inventories and representative biomass relationships. In comparison the gradient profile method involves considerable effort and instru- mentation, the site requirements are restrictive to the point of being impracti- cal, and the uncertainties in specifying the zero plane displacement height and property gradients reduce the applicability of the technique to the kind of compar- ative studies described in this report. methane hydrocarbon emission rate. Isoprene concentrations and isoprene emis- sion fluxes increased exponentially with ambient temperature. At 30°C, the emission flux of isoprene was approxi- mately 8000 jug/m2-hr. Results from the vegetation enclosure method and the micrometeorological technique generally agreed well. The emission rates of alpha-pinene measured by the gradient profile method in a northwestern Douglas fir forest were closely correlated with relative humidity. Emission rates measured by the enclo- sure method were not correlated with ambient relative humidity. However, emis- sion rates from wet branches were an order of magnitude higher than those for dry branches. Alpha-pinene emission rates from both wet and dry branches increased exponentially with increasing ambient temperature. During the autumn sampling period, the emission flux of alpha-pinene was less than 1800/ug/m2- hr and ambient concentrations above the canopy were less than 1 jug/m3. In this case, the range of fluxes observed via the enclosure and gradient methods over- lapped, but the median values did not agree within the estimated uncertainty of the methods. The comparisons between the two measurement techniques in Pennsylvania and in Washington have better defined the reliability of emission flux estimates. Although the usefulness of data collected with the enclosure method has been disparaged previously, the results report- ------- 1,000 = 0.77 •c \ <0 i 700 40 20 Figure 3. 40 60 RH/uz (%/m/s) 80 100 120 Alpha-pinene flux as a function of relative humidity/wind speed measured by the gradient profile method. ------- 1000 500 200 100 (0 I. 50 20 10 Dry Branches = 092 10 12 14 Temperature (°Cj 16 18 Figure 4. Alpha-pinene flux as a function of ambient temperature measured by the enclosure method. ------- |