JEQ Journal of Natural Resources and Life Sciences Education
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
QUICK SEARCH:   [advanced]


     

Published online 2 February 2006
Published in J Environ Qual 35:468-478 (2006)
DOI: 10.2134/jeq2005.0112
© 2006 American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA
This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via ISI Web of Science (5)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Whicker, J. J.
Right arrow Articles by Breshears, D. D.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Whicker, J. J.
Right arrow Articles by Breshears, D. D.
GeoRef
Right arrow GeoRef Citation
Agricola
Right arrow Articles by Whicker, J. J.
Right arrow Articles by Breshears, D. D.
Related Collections
Right arrow Watershed and Landscape Processes
Right arrow Ecosystem Management
Right arrow Soil Erosion
Right arrow Fire
Right arrow Radionuclides

TECHNICAL REPORTS

Landscape and Watershed Processes

Increased Wind Erosion from Forest Wildfire: Implications for Contaminant-Related Risks

Jeffrey J. Whickera,*, John E. Pinder, IIIb and David D. Breshearsc

a Los Alamos National Laboratory, Health Physics Measurements Group, Mail Stop J573, Los Alamos, NM 87545
b Colorado State University, Department of Environmental and Radiological Health Sciences, 1618 Campus Delivery, Ft. Collins, CO 80523
c University of Arizona, School of Natural Resources and Institute for the Study of Planet Earth and Department of Ecology and Evolutionary Biology, Biological Sciences East 325, P.O. Box 210043, Tucson, AZ 85721-0043

* Corresponding author (jjwhicker{at}lanl.gov)

Received for publication March 31, 2005.

   ABSTRACT
TOP
 ABSTRACT
INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Assessments of contaminant-related human and ecological risk require estimation of transport rates, but few data exist on wind-driven transport rates in nonagricultural systems, particularly in response to ecosystem disturbances such as forest wildfire and also relative to water-driven transport. The Cerro Grande wildfire in May of 2000 burned across ponderosa pine (Pinus ponderosa Douglas ex P.&C. Lawson var. scopulorum Englem.) forest within Los Alamos National Laboratory in northern New Mexico, where contaminant transport and associated post-fire inhalation risks are of concern. In response, the objectives of this study were to measure and compare wind-driven horizontal and vertical dust fluxes, metrics of transport related to wind erosion, for 3 yr for sites differentially affected by the Cerro Grande wildfire: unburned, moderately burned (fire mostly confined to ground vegetation), and severely burned (crown fire). Wind-driven dust flux was significantly greater in both types of burned areas relative to unburned areas, by more than one order of magnitude initially and by two to three times 1 yr after the fire. Unexpectedly, the elevated dust fluxes did not decrease during the second and third years in burned areas, apparently because ongoing drought delayed post-fire recovery. Our estimates enable assessment of amplification in contaminant-related risks following a major type of disturbance-wildfire, which is expected to increase in intensity and frequency due to climate change. More generally, our results highlight the importance of considering wind- as well as water-driven transport and erosion, particularly following disturbance, for ecosystem biogeochemistry in general and human and ecological risk assessment in particular.

Abbreviations: ANOVA, analysis of variance • Big Spring Number Eight (BSNE) samplers • HDF, horizontal dust flux • LCM, linear contrast of means • LANL, Los Alamos National Laboratory • VDF, vertical dust flux


   INTRODUCTION
TOP
 ABSTRACT
 INTRODUCTION
MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
PROCESSES driving erosion and transport of soil and associated nutrients and contaminants are of central concern for a diverse set of issues related to ecosystem management and risk assessment. Soil loss and redistribution through erosion can result in significant short- and long-term effects on human and ecosystem health (Whicker and Schultz, 1982; Saxton, 1995; Schlesinger and Pilmanis, 1998; Griffin et al., 2001; Toy et al., 2002; Whitford, 2002; Kaiser, 2004). Accelerated soil loss is considered to be one of the most pressing environmental problems, particularly in parts of Africa, South Asia, and Latin America (Toy et al., 2002; Kaiser, 2004), prompting some to refer to soil as an endangered ecosystem (Pimentel, 2000). Increased rates of soil erosion are largely the direct and indirect results of anthropogenic land use practices, especially those that cause intense, large-scale disturbances to vegetation communities (McNeill and Winiwarter, 2004). Accelerated soil erosion can be associated with high rates of nutrient loss from and redistribution within ecosystems and watersheds that lower soil quality (Baker and Jemison, 1991; Tongway and Ludwig, 1997; Schlesinger and Pilmanis, 1998; Whitford et al., 1998; Aguiar and Sala, 1999; Okin et al., 2001; Ludwig et al., 2002; Toy et al., 2002; Selmants et al., 2003; Miller, 2004; Wardle et al., 2004).

Soil erosion is driven by both water and wind, but most studies of soil erosion in nonagricultural settings have focused on water-driven rather than wind-driven processes. However, recent data estimating the relative rates of water- and wind-driven erosion and transport suggest that wind-driven processes have similar, and in many cases greater, impact on loss and local redistribution of soil in semiarid grassland, shrubland, and forest ecosystems (Breshears et al., 2003). Wind erosion is also important at regional and global scales, where dust storms impact human and environmental health and could provide important feedbacks to climate variability and change (Tegen et al., 1996; Griffin et al., 2001).

Wind erosion can also be important in determining risk levels associated with contaminant transport due to potential inhalation of airborne contaminated soils (Whicker and Schultz, 1982; Larney et al., 1999). Inhalation of radioactive aerosols from the wind-driven transport of soils contaminated with low levels of radionuclides, such as 137Cs and 239Pu, has been the focus of risk assessments at many Department of Energy and Department of Defense facilities across the USA (Anspaugh et al., 1975; National Academy of Science, 1989). These risk assessments are generally made under assumptions of undisturbed environments and often do not consider the potential for increased mobility of soil-bound contaminants following severe environmental disturbances, though these are of public and scientific concern (Whicker et al., 2004a).

Soil erosion can increase greatly following environmental disturbances that lower vegetation cover and disturb surface soils due to the nonlinear relationship between ground cover and erosion rate for both wind and water erosion (Fryrear, 1985; Johansen et al., 2001). There are limited studies, however, quantifying changes in erosion rates following environmental disturbance, particularly for wind erosion (Zobeck et al., 1989; Whicker et al., 2002). Yet those studies suggest that wind erosion rates can increase by more than an order of magnitude following disturbance. Disturbances in semiarid ecosystems are common and result from phenomena such as drought, insect breakouts, disease, fire, or through anthropogenic activities such as tree harvesting (Freedman, 1995; Rodgers, 1996). Among the many types of disturbance that can reduce ground cover, wildfire stands out as perhaps having the greatest potential for rapidly reducing ground cover and thereby increasing erosion rates (Johansen et al., 2001; Nyhan et al., 2001; Whicker et al., 2002). While there are numerous studies on the response of water erosion to wildfire (Nyhan et al., 2001; Johansen et al., 2001), there are few data on wind erosion rates in semiarid forests (Breshears et al., 2003) and even fewer data evaluating post-fire changes in wind-driven erosion and transport in this important ecosystem. In the case of wildfire, post-fire climate or treatments such as seeding and mulch cover should alter the rates of recovery in ground cover (Veenis, 2000), but studies evaluating such post-fire recovery effects on wind erosion are lacking.

The effects of wildfire on soil erosion are of increasing concern in the ponderosa pine forests of the semiarid west because a combination of overgrazing and fire suppression has led to overgrown tree stands with high fuel loads that are vulnerable to catastrophic wildfires (Friederici, 2003). The forests are particularly susceptible during periods of drought (Swetnam and Betancourt, 1998), which can occur over extended intervals (McCabe et al., 2004). In summary, there are few data that allow assessment of how wind erosion might increase in fire-vulnerable, semiarid forests following wildfire, or on how persistent possible post-fire increases in wind erosion might be.

The objective of this study was to compare the relative rates of wind-driven dust flux in burned and unburned ponderosa pine forest, as related to characteristics of local vegetation, soil, micrometeorology, and burn severity. Specifically, we (i) compared horizontal dust flux across a gradient of burn severities as a function of sampling height, (ii) showed that horizontal dust flux measures are correlates of direct wind erosion as measured by vertical dust flux, and (iii) analyzed temporal trends in horizontal dust flux for a period from 1 to 3 yr following the fire and related the trends to the recovery of forest floor vegetation and litter cover.


   MATERIALS AND METHODS
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Study Sites
Our study was conducted on the Pajarito Plateau located in the Jemez Mountains of northern New Mexico, USA, near Los Alamos. In May of 2000, the Cerro Grande fire burned more than 19 000 ha of mostly dense ponderosa pine forest in the Jemez Mountains. This fire also burned approximately 3000 ha of the LANL (Los Alamos National Laboratory; Los Alamos National Lab., 2000). Numerous locations on LANL property with the potential to be contaminated with either radioactive (e.g., depleted U and Pu) or chemical materials (Fresquez et al., 1998; Veenis, 2000) were burned (Kraig et al., 2001), and because the fire dramatically removed live tree density and ground cover (vegetation and litter), it exposed these sites to water and wind erosion (Fryrear, 1985; Johansen et al., 2001; Pinder et al., 2004). Wind transport in burned areas could be an especially important mechanism for redistribution, concentration, and off-site transport of these contaminated soils (Anspaugh et al., 1975; Whicker and Schultz, 1982; Gonzales et al., 2001; Breshears et al., 2003). While the initial levels of public intakes of these contaminants during the Cerro Grande fire has been evaluated and considered to have presented little health risk (Los Alamos National Lab., 2000; Kraig et al., 2001), there is still concern about the long-term, post-fire effects on continued dust emissions from LANL (Los Alamos National Lab., 2002).

The study was conducted at six sites along the western edge of LANL (35°52' N lat; 106°21' W long) at an elevation of about 2400 m (Fig. 1). Measurements of wind-driven dust were made from May 2001 through July 2003. The sites are on relatively flat mesa tops with slopes <10% (Breshears et al., 2003). Vegetation was dominated by ponderosa pine before the fire. To ensure that locally eroding material was captured, the sampling sites were centrally located within a relatively homogeneous area of vegetation with at least 100 m downwind distance along the dominant wind directions (from the west, southwest, or northwest) from any road or new disturbance. The criteria for the 100-m distance was based on data from Baldocchi (1997) that suggested that most of the eroding soil in a thick forest is collected at a distance of 10 m (analog for unburned forest), and the distance should be <100 m for vegetated but treeless areas, which is an analog for the severely burned area (Stout, 1990; Zobeck et al., 2003). Therefore, the selected distance of 100 m was considered to be conservative for all sampling sites. A line transect was established in the center of each sampling site along which two or three dust sampling stations were positioned every 30 m.


Figure 1
View larger version (17K):
[in this window]
[in a new window]
  		Fig. 1. Location of sampling plots relative to New Mexico, USA, and New Mexico State Highway 501, which borders the west side of Los Alamos National Laboratory. The burn gradient is represented as variously shaded boxes. Within each burn category, the primary sampling sites are represented as larger rectangles and secondary sites represented as smaller rectangles. The map also shows the locations of the TA-6 meteorological tower and the sites used to establish relationships between HDF (horizontal dust flux) and VDF (vertical dust flux).

 
For the purposes of this study, the range of forest fire damage was classified into three broad categories including unburned, moderately burned, and severely burned (Pinder et al., 2004). Unburned sites were those where there was no evidence of fire damage. Moderately burned sites were those where the fire had consumed the ground vegetation and litter and had scorched but not consumed the pine needles. Some tree mortality occurred at moderate sites. At severely burned sites, the fire had consumed the ground vegetation and litter and had extended up into the canopy, consuming all the needles and killing the vast majority of the trees. Remediation practices of seeding with annual and biennial forbs and mulching with straw were applied to these severely burned sites within a few months following the fire (Veenis, 2000). Within each burn category, two sites were selected to assess variability within burn categories. Because of equipment limitations, one site was termed the primary site and had greater sampling intensities than the other site, which was termed the secondary site.

Site Characterization
To document the magnitude of fire effects and measure the rate of system recovery for important variables impacting dust flux, measurements of vegetation cover, litter cover, tree canopy cover, and soil texture were made at each site (Pinder et al., 2004). Fire effects on individual trees were measured as mortality and canopy scorching (i.e., the percentage of the tree's canopy with scorched needles). Fire recovery was measured as the development of litter and vegetation cover on the soil. Micrometeorological measurements were made in the primary severely burned and unburned plots, and also at the TA-6 (Technical Area-6) meteorological station operated by LANL, which was within 3 km of all sampling areas.

Fire Effects on Forest Structure
Tree structure in burned and unburned sections of forest was characterized by measuring tree density and tree heights in 100-m2 plots. There were five randomly selected plots within each primary site, and four plots within each secondary site. Tree heights were measured using a forester's altimeter. The percentage of live tree canopy cover was measured in June 2001 using a concave mirror spherical densiometer that reflects a vertical view of the tree canopy (Lemmon, 1956). The spherical densiometer measurements were made at 15 and 12 locations in primary and secondary sampling sites, respectively.

Recovery of ground vegetation and litter cover was assessed using digital photographs of 1-m2 plots taken at early spring, midsummer, and late summer for the length of the project (June 2001–August 2003; Pinder et al., 2004). These pictures were taken from a height of approximately 2 m and at an angle from the vertical of about 25°. The 3.2 megapixel images were magnified (≥2x), and the cover percentage for each of the vegetation–cover plots determined by visual estimation of the percentage of cover of live vegetation and litter within 25 subplots. Vegetation and litter cover were separately evaluated within each subplot. A separate measure of the total cover was estimated from 81 points dispersed throughout the 1-m2 plot (Pinder et al., 2004).

Soil Sampling
Separate soil samples were collected in July 2003 from five random locations within each primary site and three random locations in each secondary site. Twenty 20-mm-deep and 38-mm-diameter cores were collected around a 1-m2 frame at each soil sampling location and composited. Litter cover was carefully scraped aside to collect only the upper horizon of soil. The soil was then air dried and submitted to the Colorado State University's Soil, Water, and Plant Testing Laboratory for measures of sand, silt, and clay percentage by mass.

Meteorological Conditions
Meteorological data collection stations (Davis Instruments Corp., Hayward, CA), which measured wind velocities at 1 m, precipitation, and temperature in 15-min averages, were installed at the primary severely burned and the primary unburned sites. Additional meteorological data were taken from a LANL-operated 46-m-tall meteorological tower located within 3 km of all the sampling sites. The LANL TA-6 meteorological station provided 15-min-averaged data on a number of meteorological conditions such as precipitation amounts and wind velocities at 11.5, 23, 46, and 92 m (Rishel et al., 2003).

Horizontal Dust Flux Measurements
Horizontal dust flux (HDF) is a measure of material transport associated with wind erosion (Breshears et al., 2003) and was measured using BSNE (Big Spring Number Eight) samplers (Fryrear, 1986). A tail fin attached to each BSNE sampler orients their 10-cm2 opening into the wind where airborne dust enters and is deposited onto a collection pan as winds decelerate through wider portions of the sampler. The BSNE field dust collectors have been extensively tested and have good sampling efficiency for soils with high fractions of sand and silt (Fryrear, 1986; Goossens and Offer, 2000) such as those that are abundant at LANL (Nyhan et al., 1978). The HDF was calculated by using the mass of the dust collected (dried at 50°C to a constant mass) divided by the area of the sampler opening (10 cm2) and the length of the sampling period. The units of HDF were then expressed as grams per square meter per day. This value represents a measure of the wind-driven mass of dust flowing horizontally along the earth surface at a particular height as measured by a single sampler. It should be noted that some other researchers have limited the term dust to describe the fine soil fraction but not the larger sand and silt fraction (Goossens and Riksen, 2004). However, others have used "dust" in a more generic sense to mean a suspension of nonspecific solid particles in air that have been formed by the mechanical disintegration of a parent material and whose size generally ranges from 1 to 10 000 µm (Hinds, 1982; Pye, 1987), although the particles collected in this study were suspended a significant period of time and have diameters generally <500 µm (Hinds, 1982; Pye, 1987; Stout and Zobeck, 1996). Therefore, we use the term dust in a more generic sense, without reference to particle size or composition (e.g., mineral or organic).

Multiple BSNE samplers may be mounted at different heights to measure gradients in HDF with increasing distance above the soil surface. Lack of highly time-resolved measurements of particle size and wind velocity in the sampling areas precluded corrections for collection efficiency, so a collection efficiency of 100% was assumed and will bias the measured HDF estimates low.

Three dust sampling stations were used in the primary sampling sites and two were used in the secondary sampling sites. A typical BSNE sampling station consisted of three BSNE samplers at heights of 0.25, 0.5, and 1 m. Additionally, recognizing that others have found sharp gradients of mass flux with sampling height (Zobeck et al., 2003), a limited number of BSNE measurements were made at heights of 0.07 and 0.18 m during the second and third year of the study. This allowed better assessment of HDF profiles with sampling height because most wind-driven dust mass occurs near soil surfaces, with about 99% occurring below a height of 1 m (Stout and Zobeck, 1996). A nonlinear least-squares method was used to determine the parameters for an exponential relationship between height and the height-specific average HDF. This exponential relationship was then integrated to determine the total horizontal mass passing through a 1-m-wide length of ground surface (Stout and Zobeck, 1996; Gillette et al., 1997). These self-orienting samplers also integrate horizontal flux for all wind directions, so the integrated horizontal flux represents the horizontal flux passing through a rotating 1-m-wide gate up to a height of 1 m (Breshears et al., 2003).

Measurements of HDF were collected beginning in May 2001 and extended through July 2003. Samples were collected generally at weekly to biweekly intervals in 2001 from June through November, in 2002 from January through November, and from March through July in 2003. Many of these collections were obtained from the same or overlapping intervals in different years. Constraints of resources imposed the midyear to midyear sampling program, and the complicating effects that seasonal weather patterns of precipitation and wind speeds would have on HDF data were not immediately appreciated. The effects of these seasonal influences on HDF and the complications imposed on comparing data among years will be discussed below.

Relationships between the HDF and meteorological conditions were established using 1-m wind and precipitation data from the meteorological station and HDF measurements from three nearby BSNE sampling stations. These HDF measurements were made from May 2001 to September 2002 in an open, unforested area.

Vertical Dust Flux Measurements
Contrary to HDF, vertical dust flux (VDF) measurements made at heights >1 m indicate dust that is available for long-distance transport and generally consists of smaller, respirable particles that are important for health risk assessments (Stout and Zobeck, 1996; Breshears et al., 2003). Therefore, VDF was measured at the TA-6 meteorological tower, which is located in the middle of an ~200-m opening in the forest that contains mostly low-lying (~0.5 m in height) vegetation (Fig. 1). The VDF was calculated as the product of the eddy diffusivity coefficient (Kz) and the concentration gradient in the vertical direction (Blackadar, 1997), with Kz being equal to the product of the von Karmon constant (0.4), the sampling height of the wind velocity measurements, and the friction velocity (u*). Friction velocities were measured at a height of 11.5 m by a sonic anemometer at the TA-6 tower (Rishel et al., 2003) and the dust mass concentrations were measured by drawing air through glass-fiber filters placed at 1- and 3-m sampling heights as described in Whicker et al. (2002). Friction velocity was measured between August 2000 and August 2001, and a subset of friction velocity measurements obtained at night, when neutral atmospheric stability was more likely, was used in these calculations.

The VDF measurements were statistically correlated to HDF measurements made from the BSNE samplers also located near the TA-6 meteorological tower using data from overlapping time periods. The VDF–HDF sites, the Davis meteorological station, and TA-6 meteorological station were arranged in a triangular pattern with approximately 200 m per side (Fig. 1).

Statistical Analyses
Statistical comparisons of HDF among burned and unburned areas were performed using a mixed-model ANOVA (analysis of variance) where burn severity, sampling periods, and the interaction of burn severity and sampling dates were fixed effects, and sites and specific BSNE locations within sites were random effects (Milliken and Johnson, 1984). Statistical computations were performed using the Type 3 option of PROC MIXED of the SAS System (Littell et al., 1996). Hypotheses of specific effects, such as (i) the effects of specific burn severities and (ii) changes in the effects of burn severities among years, were tested using F ratios computed for LCMs (linear contrasts of means; Littell et al., 1996). Due to the differential effects of precipitation and wind speeds on HDF at different elevations, separate analyses were performed for the 0.25-, 0.5-, and 1.0-m elevations, and separately for all and for dry sampling periods. The relationship between the HDF measurements made during dry sampling periods and sampling height was determined using a simple exponential function and integrated to obtain estimates of total horizontal flux through a 1-m-wide vertical plane. This analysis also allowed relative comparisons with other wind and water erosion measurements reported in similar studies.


   RESULTS
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Impact of Fire on Tree and Ground Cover
Measurements made in 2001 (details in Pinder et al., 2004) show large decreases in tree canopy cover along the burn severity gradient, with almost no live tree cover in the severely burned area, and tree cover in the moderately burned about 50% of that found in the unburned location (Table 1). Our sampling locations also captured some of the heterogeneity of tree density and tree cover typically found in ponderosa pine forests, as we found significant differences in tree cover within burned and unburned types.


View this table:
[in this window]
[in a new window]
  		Table 1. Site characterization data for the burned and unburned plots. Values provided are means ± 1 standard deviation.

 
Statistical comparisons for vegetation cover showed significantly more vegetation cover in the severely burned areas than the moderately burned and unburned areas, which was due to both viable biological residual and the extensive remediation efforts in the severely burned areas. While the remediation effort probably increased the initial surge in vegetation cover following the fire, we found that the vegetation cover in the severely burned sites decreased each year from 2001 to 2003, which is probably the result of ongoing drought persisting throughout the study period (Pinder et al., 2004).

Statistical comparisons for litter cover showed that there was significantly greater litter cover in the unburned sites (mostly dead pine needles) than that found in either the moderately or severely burned sites. Litter cover in both the moderately and severely burned sites increased from 2001 to 2002, with no further increase from 2002 to 2003 (Pinder et al., 2004). Litter cover on moderately burned sites was composed of recently shed, scorched pine needles over bare earth that was exposed by the fire.

Trends in total cover (vegetation plus litter) showed greatest total cover in the unburned sites, next was the moderately burned sites, and finally, the lowest total cover was found in the severely burned sites (Table 1). Total cover increased from 2001 to 2002 in both the severely and moderately burned locations mostly from additional litter accumulation. There was no change in total cover with time in the unburned location (Pinder et al., 2004).

Soils are dominated by sand and silt; textural classes represented are loams, silt loams, and a sandy loam (Table 1). A statistical comparison of soil texture showed no significant differences in soil composition between the different burn locations and there was more variation within burn types than among burn types (Pinder et al., 2004).

Impact of Fire on Horizontal Dust Flux
The HDF varied across time and among burn categories (Fig. 2), with a generally increasing trend in HDF with burn severity (Table 2). The magnitude of the differences in HDF among burned sites declined with sampling height (Fig. 3a).


Figure 2
View larger version (19K):
[in this window]
[in a new window]
  		Fig. 2. Mean HDF (horizontal dust flux) at 1-m sampling heights categorized by burn type.

 

View this table:
[in this window]
[in a new window]
  		Table 2. The HDF (horizontal dust flux) on burned and unburned areas in 2001, 2002, and 2003. Data are shown as means ± 1 standard deviation and are computed for data from both primary and secondary sites from all sampling periods. Data in parentheses are derived from dry sampling periods (<1.3 mm precipitation).

 

Figure 3
View larger version (12K):
[in this window]
[in a new window]
  		Fig. 3. Mean and 1 standard error for (a) HDF (horizontal dust flux) categorized by sampling heights and averaged across all sampling periods and (b) HDFs that were averaged across dry sampling periods with precipitation <1.3 mm.

 
The ANOVA showed that the HDF measured on severe burn areas were significantly (all LCM F values > 4.5; df = 1, >465; P < 0.05) greater than those on unburned areas in all years and at all sampling heights. The profile of HDF with elevation showed greater increases at the 0.25-m height than at the 1.0-m height for all 3 yr. For example, mean (±standard error of mean) increases in HDF in 2003 were 7.79 ± 0.74, 2.52 ± 0.19, and 2.04 ± 0.19 at 0.25-, 0.50-, and 1.0-m heights, respectively.

The response of HDF to the moderate burn was less distinct than that for the severe burn, but the mean HDF on the moderate burns was always greater than that on the unburned areas. Mixed-model ANOVA and LCM analyses based on these data did not indicate any significant difference between moderate and unburned areas; however, these comparisons were compromised by skewness in the data (Table 2), as indicated by standard deviation to mean ratios being >0.3. The LCM analyses on logarithmic transformations of HDF, which moderate the effects of skewness, indicated significantly (P < 0.05) greater HDF on moderate burn sites at all time periods. The HDF measurements <0 were arbitrarily redefined as 0.03 before logarithmic transformations. Mean HDF on severe burn areas were significantly greater than those on moderate burn areas at all heights (all LCM F values > 19; df = 1, 4; P < 0.05). The major difference between the response of HDF on severe and moderate areas was the greater increase on severe areas at the 0.25-m elevation. For example, the mean (±standard error) differences between severe and moderate burn areas in 2003 were 6.56 ± 0.68, 1.58 ± 0.18, and 0.94 ± 0.19 at 0.25, 0.50, and 1.0 m, respectively.

Analysis of Fire Effects on Horizontal Dust Flux during Dry Periods
Field observations of mud splash on the outside of the lowest BSNE samplers suggested the potential for additional collection of soil in the sampler from rain splash, thus prompting an analysis of HDF for both all sampling periods and periods with little or no rain (Table 2, Fig. 3a and 3b). The potential impact of precipitation and rain splash on HDF measurements was further investigated by correlating HDF measurements and total precipitation. The HDF measurements for this investigation were made adjacent to the TA-6 meteorological tower and the correlation showed that, as total precipitation in the sampling interval increased, there were corresponding increases in HDF at 0.25 m (F = 17.3; df = 1, 46; P < 0.01) and 0.5 m (F = 15.9; df = 1, 45; P < 0.01) as shown in Fig. 4, Column a. The increasing HDF with increasing precipitation probably results from rain splash during rain events. Although no direct measurements of rain splash were made, it is a plausible explanation based on (i) observations of mud splatter on sampler surfaces, (ii) rain splash being an effective mechanism for soil particle transport (Finkle, 1986), (iii) the possibilities that splashed particles could enter samplers through either the BSNE front opening or open top vent, and (iv) the fact that there was no statistically significant (F = 0.32; df = 1, 46; P ≥ 0.05) increase in HDF at 1.0 m (a height that is generally greater than the height of rain splash; Finkle, 1986) with increasing precipitation. There was, however, a statistically significant (F = 6.1; df = 1, 46; P 0.05) increase in HDF at 1.0 m with increasing mean daily maximum wind speed (Fig. 4, Column b). Wind speeds had no statistically significant (P > 0.10) effects on HDF at 0.25 and 0.5 m, possibly due to the overriding effects of rain splash.


Figure 4
View larger version (24K):
[in this window]
[in a new window]
  		Fig. 4. The means of the horizontal dust flux as a function of precipitation (Column a) and wind velocity (Column b). The relationships are plotted individually for sampling heights of 25, 50, and 100 cm.

 
To separate the effects of wind and rain splash, the HDFs were additionally compared among burned and unburned areas during dry periods. For this analysis, we selected sampling periods where precipitation during the sampling interval (generally 1–2 wk) was <0.13 cm. This precipitation threshold was considered to be low enough to have very little effect on the HDF (Fig. 4, Column a), yet provides enough sampling periods for robust statistical analysis. Although fewer sampling periods are available for statistical comparisons, the data indicate greater HDF on burned areas during dry periods than on unburned areas (Fig. 3b), consistent with the analysis using all data as shown in Fig. 3a. At 0.25 m, the HDF measures across all 3 yr were significantly greater on severe burn (LCM F = 368.2; df = 1, 3; P ≤ 0.01) and moderate burn areas (LCM F = 56.5; df = 1,3; P ≤ 0.01). At 0.5 m, the HDF was significantly greater on severe burn areas (LCM F = 24.8; df = 1, 3; P ≤ 0.05) but not significantly greater on moderate burn areas (P > 0.10). At 1.0 m, the HDF were statistically significantly greater on both severe areas (LCM F = 703.6; df = 1, 3; P ≤ 0.01) and moderate areas (LCM F = 84.5; df = 1, 3; P ≤ 0.01).

The mean HDF for dry periods did not dramatically decrease between sampling heights of 25 to 100 cm as seen in the means from all sampling periods (Fig. 3b). We used the limited set of HDF measurements at heights of 7 and 18 cm in each of the burned and unburned primary sites to determine a best-fit exponential relationship between sampling height (h) and mean HDF using data from the dry periods (to avoid potential rain-splash effects). Integrating these functions over the five sampling heights provided the total mass flux per unit width. A simple exponential model of the form Y = a exp(bh) was used to describe this relationship, with a and b being estimated using an iterative least-squares method (StatSoft, 1994). The values for a and b and the associated standard errors were 2.64 ± 1.19 and –0.87 ± 0.29, respectively, for the severely burned area, and 1.36 ± 1.21 and –0.65 ± 0.32, respectively, for the moderately burned location. The relationships were statistically significant (P < 0.05), but the coefficients of determination (R2) for these models were generally low and were 0.14 and 0.07 for the severely and moderately burned areas, respectively. Integrating these functions to a 1-m height gave an estimated total mass passing through a 1-m-wide length of 1.76 and 1.00 g m–1 d–1 for the severely burned and moderately burned areas, respectively. The height profile of HDF in the unburned area did not vary significantly with sampling height, so a comparable estimate of HDF is simply the flux averaged across all three heights multiplied by the range of sampling heights (1 m). This gives an estimated mean of 0.74 g m–1 d–1 passing through the similar 1-m-wide length in the unburned areas.

Relationship between Horizontal and Vertical Dust Flux
A significant linear relationship was found between the HDF measurements (dependent variable) and the VDF measurements (P < 0.05, R2 = 0.37). The intercept and slope with standard errors were 0.02 ± 0.02 and 0.05 ± 0.02, respectively. A plot of this relationship is shown in Fig. 5. This correlation suggests a link between HDF and VDF and the erosion and long-distance transport of soil and soil-bound contaminants that is implied by measures of VDF.


Figure 5
View larger version (10K):
[in this window]
[in a new window]
  		Fig. 5. Relationship between the VDF (vertical dust flux) measurements and the HDF (horizontal dust flux) measured at the TA-6 meteorological station. The regression equation is VDF = 0.02 + 0.05(HDF), R2 = 0.37.

 
Forest Recovery and Horizontal Dust Flux
There was little indication of declining HDF with time since the fire. Measures of HDF were similar among years or showed trends of increasing HDF. The LCM comparing overlapping sampling intervals in severely burned areas indicated no increase in HDF from 2001 to 2002 at any elevation (P > 0.05) but significant increases in HDF between 2002 and 2003 at 0.25 m (LCM F = 7.06; df = 1, 470; P < 0.01) and 0.50 m (LCM F = 18.1; df = 1, 470; P < 0.01). There was no significant increase from 2002 to 2003 at 1.0 m (P ≥ 0.05). Linear contrasts of means using logarithmically transformed data from the moderately burned plots indicated no increase in HDF from 2001 to 2002 at any elevation (P ≥ 0.05) but significant increases in HDF between 2002 and 2003 at 0.50 and 1.0 m (LCM Fs > 7; df = 1, 470; P < 0.01).


   DISCUSSION
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
CONCLUSIONS
 REFERENCES
 
There have been numerous studies on the effects of forest fires on water erosion and runoff, but this study is among the first to document wildfire effects on wind erosion in a forest ecosystem and one that is particularly vulnerable to wildfire (Covington, 2003; Friederici, 2003). Our results suggest that wildfire increases wind erosion. Specifically, HDF was significantly increased in burned areas, especially in the severely burned areas. The likely causes of increased HDF in the burned areas are decreased vegetation, litter, and tree cover that resulted in higher surface wind velocities (Whicker et al., 2002).

The use of BSNE samplers in this study differs somewhat from their use in previous studies. We were interested in obtaining and comparing measures of average dust flux occurring in differentially burned areas. The samplers allow collection during extended periods of time from representative areas. In some past studies, the samplers were used for short durations under well-monitored specific conditions to obtain event-based measurements to quantify relationships between factors that affect the entrainment and suspension of particles. Our finding that the impact of the fire on HDF was most pronounced at the lower sampling heights was probably affected by rain splash at the 0.25- and 0.5-m sampling heights, and rain splash should be considered in future studies. Despite the potential effects of rain splash on HDF measurement, the statistical conclusion of increased HDF in burned areas held during dry periods, but the effect of forest burning was less. Increased rain splash in the burned areas may be explained by the facts that the burned areas had relatively more bare soil and the trees in the unburned areas moderate the effect by reducing raindrop sizes and velocities in the unburned areas. Regarding contaminated soils, increased rain splash and post-fire vegetation in burned areas suggests higher transfer rates of the contaminant to plant surfaces and subsequent ingestion by grazing animals (Dreicer et al., 1984). The lack of a pronounced decrease in HDF with sampling height during the dry periods was unexpected (Stout and Zobeck, 1996) and may be explained if a large fraction of the collected dust is not made up of saltating particles, but rather "background" atmospheric dust from distant sources filtered through the forest.

The possibility that much of the collected dust could be from distant sources is also supported by the percentage of total ground cover (Table 2), which is especially large relative to many wind erosion studies in agricultural settings. Results from Siddoway et al. (1965) and Gregory (1982) suggest that wind erosion is essentially controlled in agricultural fields with total ground cover >50%, as was measured in this study. There are few studies, however, in nonagricultural systems to predict how general relationships developed in agricultural systems translate into other environments. For example, the fact that total cover was <50% may limit (not eliminate) wind erosion in more natural settings due to the patchy nature of vegetation and cover in these landscapes.

We were not able to begin measurements for this study until ~1 yr after the fire, and we may have missed important information on the full impact of the fire on dust flux. Specifically, one would expect that the effects of the fire on soil erosion would be greatest immediately following the fire because of significant loss of vegetation and litter cover and lack of recovery time. For the time period studied, however, we found that increased HDF in the burned areas persisted through July of 2003, 3 yr post-fire, despite extensive remediation (Veenis, 2000). Comparisons of burn severities on tree and ground cover indicated that the effects of the Cerro Grande fire was similar to that for other fires in ponderosa pine forests, but the rates of recovery of ground vegetation and litter were slower than those for other fires (Pinder et al., 2004). The slower recoveries, which occurred even on those areas that were seeded and mulched, are probably due to the effects of ongoing drought.

Ecological Implications
Airborne dust is constantly being transported through environments (Pye, 1987), and this study showed that disturbances such as wildfire can significantly increase dust flux, and the increase can persist for years. This study reports estimates of 1-m integrated dust flux values of slightly more than 1 kg of dust passing through each meter in a burned forest per year, about 2.4 times that found in unburned forests and about five times that found for water erosion on shallow slopes (Johansen et al., 2001; Breshears et al., 2003). This could be a significant amount of increased dust mass when integrated over large tracts of forested land and over multiple years. The soil particle sizes and chemical constituents of the dust were not measured in this study, but the dust may be expected to be enriched in soil fines containing soil organic materials and associated soil nutrients, both important components for soil productivity. Future work should consider the chemical makeup of the dust to more fully understand the ecological impact of this enhanced redistribution of soil.

Soil is a key integrating factor in forest ecosystems because of the strong coupling between soil quality and the biotic components in the ecosystem that it supports. Loss of soil and associated nutrients is considered an important indicator of ecosystem function (Tongway and Ludwig, 1997; Whitford et al., 1998). The response of wind erosion to wildfire suggests that it could also respond to other disturbances; thus, wind erosion may be an important metric to be included in the assessment of impacts of ecosystem disturbance, forest management, and the trajectory of the recovery (Miller, 2004; Whicker et al., 2004b; Neff et al., 2005).

Post-Fire Transport and Erosion by Water and Wind
Wildfire increases both wind and water erosion, but there are few, if any, comparisons of the relative effects. Johansen et al. (2001) studied water erosion on mild slopes (4–8%) using rainfall simulation on unburned plots and plots burned over by the Cerro Grande fire (in the same primary severely burned plot used in this study). Using the data from the Johansen et al. (2001) water erosion study and this wind erosion study, the effects of forest fire are compared qualitatively for severely burned areas.

The amount of increased wind erosion is estimated using the mean HDF values measured at 1 m during the study (averages of 0.73, 1.18, and 1.76 g m–2 d–1 at the unburned, moderate burn, and severe burn areas, respectively) and the empirically determined relationship between measured HDF and the VDF (Fig. 5). Estimated annual VDFs were 21 g m–2 yr–1 [e.g., 365 d yr–1(0.020 + 0.73[0.053])] for unburned forests, 30 g m–2 yr–1 for moderately burned forests, and 41 g m–2 yr–1 for severely burned forests. These estimates of upward net VDFs are most applicable to the moderately and severely burned locations because of similar ground and tree cover to the TA-6 location where they were made. The vertical flux estimates suggest an approximate 95% relative increase in vertical flux and an absolute increase of 20 g m–2 yr–1 in the severely burned area. Though a significant increase, the measured erosion rate in the burned forest is still quite low relative to agricultural fields, which can also be considered a disturbed environment (Lal et al., 2003).

The 95% increase in wind erosion appears small when compared with the >20x post-fire increase in sediment yield (approximate means of 3–70 kg ha–1 mm–1) through water erosion observed at LANL by Johansen et al. (2001). Breshears et al. (2003) used (i) the frequency and magnitude of precipitation events for the LANL site and (ii) the rainfall simulation plot data of Johansen et al. (2001) on an unburned forest plot to estimate water-driven erosion rates of approximately 1 g m–2 yr–1 before the fire. Assuming a 20-fold increase in water erosion in severely burned areas gives an estimated water erosion rate of 20 g m–2 yr–1 after the fire. This would suggest an increase of 19 g m–2 yr–1 due to water erosion following a fire; however, this water erosion study was conducted within a few months of the Cerro Grande fire and a significant fraction of the collected mass was ash. Our wind erosion measurements started the following year, when much of the ash had eroded away and some vegetation recovery had occurred. Thus, this comparison is likely to be biased toward greater water erosion.

The previous assessment of the impacts of fire on wind erosion rates is based on two assumptions concerning the HDF measures in the unburned forest. The first is that the relationship between HDF and VDF observed for the TA-6 location is also valid for the forest. It is reasonable to challenge this assumption. The levels of canopy coverage and the tree heights (Table 1) would suggest less vertical air movement and mixing. Thus, the ratio of VDF to HDF may be less for unburned forests. The application of the relationship shown in Fig. 5 to severe burn areas is more plausible because of the loss of the effects of canopy coverage and foliage mass on wind mixing following the fire. The second assumption is that all the material measured as HDF in the forest is local erosion. This assumption may not be valid. The collected material may be partly or mostly deposition of airborne particulates from distant sources filtering through the canopy. Failure of either of these assumptions results in overestimation of vertical flux for the forest and a corresponding underestimation of the increased rate of wind erosion. Thus, the comparison of the increases in water and wind erosion rates following fire as presented above may be biased and underestimate the relative importance of wind erosion. This underestimation may be as much as 200% [100% (41 – 0)/(41 – 21)].

Implications for Transport of Soil-Bound Contaminants
Rates of wind-driven transport and erosion underlie any assessment for human or ecological risk where the wind-driven pathway is of concern. Our results provide data on such rates in a major ecosystem type for the western USA-ponderosa pine forest-and provide some of the first estimates of how much those rates increase in response to a major type of ecosystem disturbance-wildfire. The specific risks will be dependent on other additional factors such as size distribution of dust particles, the associated size distribution for contaminant particles, and the contaminant chemical properties (Whicker and Schultz, 1982), the later two of which will vary with contaminant type. Nonetheless, our findings provide underlying rates relevant for assessing risk for a diverse suite of contaminants. Post-fire water erosion at our site was related to 137Cs transport (Johansen et al., 2003); we expect that post-fire wind erosion will be similarly related to transport of 137Cs and other soil-bound contaminants, although this will require future study.

Regarding contaminant transport, wildfires probably increase local redistribution of contaminated soils through wind, and the effects can last for years depending on the recovery of the vegetation and litter cover. This local redistribution may concentrate contaminants in vegetation and in areas of aeolian deposition such as in areas of higher canopy cover or leeward sides of topographic objects such as hills and canyons where the wind velocity is lower (Whicker and Schultz, 1982; Ritchie and McHenry, 1990; Coppinger et al., 1991; Sutherland et al., 1991; Gonzales et al., 2001; Toy et al., 2002). Importantly, forest fires may also increase long-range (several kilometers) transport of contaminants and increase exposures to the public, as shown in the relationship between the VDF and the HDF found in this study. As an example, LANL (2002) shows that depleted U air concentrations at the LANL boundaries have significantly increased following the Cerro Grande fire, though the levels remain low and within regulated safety standards. This study, however, confirms the value of data-based risk assessments following environmental disturbance.


   CONCLUSIONS
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
REFERENCES
 
Wildfire has been shown to increase HDF, especially in severely burned locations and, to a lesser degree, in moderately burned locations. The increased dust fluxes continued for several years following the fire, probably because of continuing drought, which has slowed the recovery of ground vegetation. The finding of increased HDF and VDF implies that wind-driven mobility of soil-bound contaminants may have also increased. This potential increase in contaminant mobility at LANL following the Cerro Grande fire is supported by increased air concentrations of depleted U, which has been distributed in the local environment through explosive testing (Los Alamos National Lab., 2002). On a broader scale, the results of this study have important implications for soil quality, biogeochemistry, ecosystem health, and risks associated with contaminant transport.


   ACKNOWLEDGMENTS
 
This work was funded, in part, by the Technology, Development, Enhancement, and Application program at Los Alamos National Laboratory and the remainder was funded by the Department of Energy under contract W7405 ENG-36. The support of Shawna Eisele, Jeffrey Hoffman, Dr. Robert Devine, and Dr. Robert Murphy is greatly appreciated. We would like to thank Kristine N. Baker for data collection and the following individuals for their valuable contributions to this study: Dr. Shawki Ibrahim, Dr. Jeff Collett, Dr. John Zimbrick, and Dr. Ted Zobeck.


   REFERENCES
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 





This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in Pu