An Evaluation of Modeled Plume Injection Height with Satellite-Derived Observed Plume Height

Plume injection height influences plume transport characteristics, such as range and potential for dilution. We evaluated plume injection height from a predictive wildland fire smoke transport model over the contiguous United States (U.S.) from 2006 to 2008 using satellite-derived information, including plume top heights from the Multi-angle Imaging SpectroRadiometer (MISR) Plume Height Climatology Project and aerosol vertical profiles from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP). While significant geographic variability was found in the comparison between modeled plumes and satellite-detected plumes, modeled plume heights were lower overall. In the eastern U.S., satellite-detected and modeled plume heights were similar (median height 671 and 660 m respectively). Both satellite-derived and modeled plume injection heights were higher in the western U.S. (2345 and 1172 m, respectively). Comparisons of modeled plume injection height to satellite-derived plume height at the fire location (R = 0.1) were generally worse than comparisons done downwind of the fire (R = 0.22). This suggests that the exact injection height is not as important as placement of the plume in the correct transport layer for transport modeling. OPEN ACCESS Atmosphere 2012, 3 104


Introduction
During the past decade, there has been an increase in both the number and intensity of wildfires in the western U.S. and the wildfire season is expected to lengthen as temperatures warm [1].In addition, U.S. air quality regulations have been tightening.Due to these regulations, attention has been focused on modeling and forecasting smoke impacts from wildfires and efforts are underway to develop and improve smoke modeling and forecasting systems [2][3][4][5][6] that can accurately predict smoke impacts on visibility and concentrations of pollutants of concern to human health.
For this paper, we define plume injection height and plume top as the vertical zone in which a buoyant plume begins to transport horizontally away from its source, and the top of that zone respectively.Plume injection height is important for modeling wildfire smoke plume transport, footprints, and surface concentrations.The injection height of the plume determines the type of winds and turbulence the plume will experience.For example, extremely buoyant wildfire smoke columns produce plumes that break through the atmospheric boundary layer and are transported hundreds to thousands of kilometers downwind in the free-stream transport winds [7,8].In contrast, smoke plumes remaining within the boundary layer often become well-mixed in regions near a fire [9].A wildfire smoke column may have several plume injection heights and this vertical distribution of plume mass is a key input for dispersion and air quality models.Model results, such as surface concentrations, are sensitive to the amount of plume mass injected at various heights because the associated winds and turbulence directly influence the plume footprint and dilution.Accurately characterizing smoke plumes and their vertical structures is necessary to produce useful local-, regional-, and national-scale smoke predictions.
Buoyant plume injection height is a function of both the heat release rate from the source and the prevailing meteorology [10].The rate of heat release is directly related to the rate of biomass consumption.Unfortunately, estimates of biomass consumption for a given fire are highly uncertain [2,5,11,12].Moreover, the amount and rate of heat release from fires varies by combustion phase making hourly heat release difficult to estimate.Even if all of the physical parameters for a specific fire are well known, modeling plume injection height is difficult due to the difference between the spatial resolution of the model and the scale of the variation of the burn activity.Finally, wildfires often occur over a large area, and within that area differences in terrain, vegetation, and fuel moisture lead to different consumption characteristics.This heterogeneity often leads to the total heat released splitting into several plume cores [13].
Previous studies have identified a need to improve the spatial and temporal characterization of heat release from fires [4,14].Efforts are underway to quantify the uncertainties associated with smoke modeling and develop algorithms that are more representative of real-world smoke transport [9,10,[15][16][17].However, due to the difficulty of acquiring information about plume height from real fires, there has been a lack of data to compare observed plume heights with model estimates.
One way to characterize smoke plume behavior is through the use of information derived from satellite data.In recent years, satellite data products were used to evaluate the vertical structure and smoke plume characteristics of wildland fires.Several research efforts have characterized smoke plumes from observations with the goal of improving model performance [18][19][20][21], but most comparisons to model outputs have been limited to case studies, with one recent exception [22].This paper describes a study that compared plume height data derived from the MISR onboard Terra and the CALIOP onboard the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) satellite against modeled plume height estimates for a large number of fires in the U.S. from 2006 through 2008.The use of both MISR and CALIOP provides complementary information [23].To our knowledge, it is the first study that compares CALIPSO data to multiple modeled smoke plumes beyond individual case studies.

Methods
Plume top height data collected by MISR and CALIOP were compared to the modeled vertical mass distribution for 227 fires and their plumes that occurred from 2006 through 2008 in the contiguous U.S.
The overall technical approach for this work involved four general steps: (1) data acquisition and processing; (2) identification of fires and corresponding CALIOP and MISR satellite observations from 2006 through 2008; (3) production of modeled smoke predictions for comparison to the satellite observations; and (4) comparison and assessment of the modeled plume height predictions to the satellite derived observations.

Data Acquisition and Processing
MISR and CALIOP measure plume heights in distinct and complementary ways.MISR relies on stereoscopic analysis of plumes based on a minimum of two angular views to extract plume heights [18].The CALIOP instrument uses active lidar to measure aerosol vertical profiles [24].
In addition to the satellite data, several other sources of data and software applications were used to acquire and prepare the data used for this study: • The Hazard Mapping System (HMS) Smoke Product [25] was used to identify when and where smoke plumes were present from 2006 through 2008 [26].• Daily fire locations from the Satellite Mapping Automated Reanalysis Tool for Fire Incident Reconciliation (SMARTFIRE) fire information system [27] were used to identify latitudinal and longitudinal locations of fires and burn area estimates for fires present from 2006 through 2008 [28].

• The Satellite Nadir Track Sensor Observation Service provided by the Vis Analysis Systems
Technologies team at the University of Alabama in Huntsville [29] was used to determine CALIPSO orbit paths and identify days when these paths intersected a smoke plume in the HMS data set.• CALIOP Lidar Level-2 5 km aerosol data acquired from the Atmospheric Science Data Center (ASDC) [30] were used to compare observed smoke plume heights to modeled plume heights [31].
• CALIOP Lidar Level-1B 532 nm wavelength attenuated backscatter data were acquired from the ASDC and used for data visualization and qualitative assessment of smoke plume characteristics and PM 2.5 concentrations.• Smoke plume summary data acquired from the MISR Plume Height Climatology Project [32] were used to compare observed smoke plume heights to modeled plume heights [33].Specifically, we used the median wind corrected plume heights.• The BlueSky framework smoke modeling system was used to develop smoke plume height and fire emission and heat estimates for comparison to the CALIOP and MISR satellite observations [2].• The Community Multiscale Air Quality (CMAQ) model was used to transport and diffuse injected smoke plumes to compare with CALIOP aerosol data [34].

Identification of Fires and Corresponding MISR and CALIOP Satellite Observations
The HMS database was used to identify satellite detectable smoke plumes from fires from 2006 through 2008.The HMS data set includes daily electronic (GIS-based) maps of smoke plume extents as determined and digitized by trained satellite analysts.HMS analysts examine operational satellite imagery, particularly from MODIS and GOES, and qualitatively delineate visible areas of smoke.The SMARTFIRE information system was used to query fire information including fire location and area burned estimates for days and times corresponding to the HMS smoke plume data.SMARTFIRE integrates and reconciles human-recorded wildfire incident data from Incident Status Summary (ICS-209) reports with satellite-detected fire data from by HMS.The HMS fire data set is a combination of automated fire detection and human quality control.It includes automated fire detections from two geostationary satellites (GOES-East and GOES-West) and five polar orbiting satellites (Moderate Resolution Imaging Spectroradiometer [MODIS] Terra and Aqua and the Advanced High-Resolution Radiometer [AVHRR] on three NOAA spacecraft).
Next, the SMARTFIRE datas were used to query the MISR plume height satellite observation database for days and times that correspond to the SMARTFIRE data set.The MISR Plume Height Climatology Project data archive contains information about North American smoke plume characteristics for 2002 and from 2004 through 2007.The native data set consists of 61 smoke plume parameters corresponding to a single geographic coordinate location.The following MISR parameters were acquired for this study: plume location (latitude and longitude) and date, median and maximum plume height, median and maximum wind corrected plume height, and median and maximum wind corrected plume height using a planar correction.The MISR plume height data were matched to fires from the SMARTFIRE database by date and location.Matched MISR-SMARTFIRE plumes were determined by selecting the nearest same-day fire location in the SMARTFIRE database (within 10 km) for fires greater than 40.5 hectares (100 acres).A total of 163 records from the MISR plume height database were matched to SMARTFIRE fire locations for comparison.
CALIOP satellite data were acquired based on days and times when the CALIPSO satellite orbit path intersected an HMS smoke plume.Daily CALIPSO orbit path data were overlaid with daily HMS smoke plume data for 2006 through 2008 and the geographic intersection points were recorded as candidate times and locations when the satellite was likely to detect a smoke plume.For days with CALIPSO-smoke plume intersections, map images were created and inspected.Because this study is concerned with plume injection height rather than plume transport, we preferentially selected days when the CALIPSO orbit path was close to the origin of a smoke plume, as determined by the SMARTFIRE fire locations.The maximum allowed distance between the CALIPSO orbit path and the fire origin was 50 km.Figure 1 shows an example of the process for identifying days and times when the CALIPSO orbit paths intersected HMS smoke plumes.

Development of Modeled Smoke Predictions
Heat and emissions from wildland fires can be estimated as shown below in Equation (1) [35]: Where: M(x) is the heat or mass emitted of species x; A is the area burned; AFL is the available fuel load including live and dead vegetation and soils in mass per unit area; β is the burning efficiency or the fraction of fuel consumed, and Ef(x) is the emission factor for heat release (energy emitted per mass of fuel consumed) or for chemical species x (mass of chemical species x emitted per mass of fuel consumed).The BlueSky modeling framework [2] was used to construct the modeling pathway.The BlueSky framework links many different model applications together to produce emissions and pollutant transport estimates from wild and prescribed fires.Figure 2 illustrates the process pathway and model applications that were used to estimate emissions and generate the aerosol mass distribution predictions to compare to the satellite observations.Figure 2. Model pathway used to develop estimates of plume height and smoke predictions for comparison to MISR and CALIOP observations.FCCS, Consume3.0,and FEPS are used as implemented in the BlueSky Framework, version 3.0.0[36]. 1

BlueSky Framework
This study employed the most recent smoke modeling systems available in the BlueSky framework to calculate each term in the emission equation.Available fuel load was calculated using the Fuel Characterization and Classification System (FCCS) [43,44], which provides a map of over 100 fuelbed models at 1-km resolution.Biomass consumed was calculated using Consume 3.0 [39], which combines fuelbed information and weather data to produce estimates of fuel consumption by burn phase (i.e., flaming, smoldering, and residual) and fuel bed strata (e.g., forest canopy, soil, ladder fuels).Daily bulk consumption was then allocated hourly using a time profile developed by the Western Regional Air Partnership (WRAP) [45].Emissions factors, including heat, were modeled using the Fire Emission Production Simulator (FEPS) [40], which modulates emission factors based on combustion efficiencies.Results of this modeling included hourly emissions and heat release rates for each fire studied.Plume injection heights were calculated using Briggs plume rise methods for buoyancy-dominated plumes [46] as implemented in FEPS.Briggs plume rise is a function of buoyancy, ambient wind speed, and static stability.FEPS computes the buoyancy parameter as 2.58 × 10 −6 times the heat release rate.Heat release is calculated as shown in Equation 1and, therefore, depends on all previous modeling steps.The wind speed was defined by the BlueSky default value of 6 mph.Pasquill stability class [47] is set to unstable (stability class A) during daytime hours (defined as 0600 to 1800 local), and ranges from slightly unstable to slightly stable (stability class C through E) during nighttime hours depending on wind speed.The Briggs stability parameter is based on the vertical temperature gradient that FEPS defines based on Pasquill stability class.The FEPS injection height calculator produces hourly estimates of the plume top for each fire.
The modeled plume top data were used to compare to median wind corrected heights from the MISR plume height climatology project.The MISR heights are based on fitting a plane to the wind corrected heights derived at each pixel within an identified plume, removing all points more than 1.5 standard deviations from the plane, and taking the median of the remaining points [33].
Because the CALIOP footprint is so narrow, it rarely passes directly over a fire.CALIOP plume detections typically occur downwind of the source fire.Therefore, it was not possible to distinguish plume injection height from plume transport in the CALIOP data set.To compare the CALIOP plume data to modeled plume data, the CMAQ model (CMAQ version 4.5.1 was used for this study) [42] was used to generate downwind plume heights and PM 2.5 concentrations for fires corresponding to the CALIOP data set.
The Sparse Matrix Operating Kernel Emissions (SMOKE) system, version 2.3, [41] was used to post-process emissions estimates from FEPS to create CMAQ-ready emissions data files that contain spatially and temporally resolved emissions estimates [48][49][50][51].The output of the SMOKE system was a field of vertically distributed fire emissions, which were merged with day-specific biogenic and anthropogenic emissions, generated using SMOKE to produce the final merged three-dimensional emissions input files for the CMAQ model.Emissions estimates of heat release and several primary pollutants (e.g., PM 2.5 , CO, sulfur dioxide [SO 2 ], NO x ) were produced for input to the CMAQ model.
The CMAQ model was used to produce three-dimensional fields of PM 2.5 concentrations from fires during 2006 through 2008 on the national 36-km Regional Planning Organization grid, covering the contiguous U.S. The model simulations are driven by meteorological data from Version 3.7 of the Fifth-Generation Pennsylvania State University/National Center for Atmospheric Research mesoscale model (MM5).The MM5 modeling was performed on a 29-layer vertical grid, with vertical resolution of 50 to 75 m in the lower troposphere gradually stretching to coarser resolution at the upper layers of the grid.The CMAQ modeling was performed on a 17-layer vertical grid, which aligns with the MM5 vertical grid.
To compare the CMAQ model output with the CALIOP observations, vertical profiles of PM 2.5 concentrations from wildfires were extracted from the model output.These concentrations were bi-linearly interpolated from the four nearest grid cell centers to the satellite orbit path locations.The CALIOP orbit paths were matched to the CMAQ modeled data by day and time.Data were eliminated for times when a one-to-one match between the CALIOP and CMAQ data could not be made.A total of N = 64 CALIOP-CMAQ data pairs were included in the analysis and the modeled PM 2.5 concentration data were compared to the CALIOP aerosol data for the data pairs.Figure 3 shows an example of the comparison of the CALIOP aerosol data and the CMAQ modeled PM 2.5 concentrations.These comparisons were made to qualitatively examine the spatial distribution, extent, and concentrations of PM 2.5 plumes produced by the CMAQ model.Plume height information from the CALIOP data was compared with plume height information from the CMAQ model.For CALIOP, plume top height was calculated as the average of the maximum heights in the aerosol feature mask product.The CMAQ model output does not explicitly contain plume top data.To allow the comparison of the CMAQ modeled data and the CALIOP plume height data, a CMAQ plume top equivalent value was derived.The CMAQ model generates predictions of threedimensional pollutant concentrations as they undergo transport, chemical transformation, and dispersion.Concentrations of PM 2.5 from smoke are typically highest slightly above the location from which they are emitted.As a smoke plume rises vertically, atmospheric dispersion and mixing cause pollutant concentrations to decrease.The CMAQ plume height equivalent was derived by identifying the vertical height at which the modeled PM 2.5 concentration was 10% of the maximum PM 2.5 concentration within the vertical column.Use of a threshold is necessary because vertical diffusion in a coarse vertical grid structure mixes small amounts of PM 2.5 to altitudes well above the bulk of the modeled plume.The selection of this threshold is arbitrary, and changing its value would affect biases in predicted plume top height; however, other threshold values investigated had minimal effect on the overall correlation.Figure 4 is a conceptual illustration of the methodology used to derive the CMAQ plume top equivalent.

Statistics and Analysis Used in Comparing Plume Heights
Statistical analyses were performed to compare the extent to which observed plume heights from the MISR and CALIOP data agreed with the FEPS and CMAQ model results, respectively.Basic descriptive statistics (e.g., mean, median, standard deviation) were calculated for each data set.For each comparison, a Kolmogorov-Smirnof test was performed to determine the statistical significance between the satellite-derived and modeled plume height data.The data were also analyzed by region to determine if there was a regional influence on the statistical results.Plume height estimates were further investigated by examining parameters that drive the plume injection height calculations, such as area burned and fuel loading.

Comparison of MISR Observed Plume Heights to FEPS Plume Heights
Results of the comparisons between the MISR observed plume heights (median wind corrected) and the FEPS modeled plume heights show that the satellite-derived plume-top height has statistically significantly (p < 0.01) higher median value than the modeled plume-top height (Table 1).In addition, the overall range of the modeled plume heights was larger than the range of observed plume heights.Examining the paired data for all the individual plumes shows that there is only a weak correlation between the data sets (r 2 = 0.1). Figure 5 shows the results of the regression analysis for the observed plume heights (x-axis) versus the modeled plume heights (y-axis).Other studies have found better correlation using MISR maximum wind corrected heights instead of median wind corrected heights [52]; however, our results using the maximum wind corrected height were slightly worse (r 2 = 0.09).While there was little correlation between MISR observed and FEPS modeled plume heights when all data values were considered, the spatial distribution of modeled and observed plume heights across the U.S. agree qualitatively.The data points in Figure 5 are broken out by three regions, delineated in Figure 6.There are several differences between the fire regimes in the southeastern and western U.S., including the type, density, and moisture of the vegetation.In addition, most of the area burned in the Southeast is the result of prescribed fires, while there is a much higher occurrence of wildfires in the western U.S. [53].Wildland fires in the southeastern U.S. generally produce lower plume heights compared to those in the western U.S. The difference between observed and modeled plume heights is larger in the mountain and Pacific Northwest regions, where the median observed plume heights are statistically significantly greater than the modeled plume heights but some modeled plume heights are several times higher than those observed by MISR.In the Southeast, there is no statistically significant difference between observed and modeled plume heights.Distribution statistics by region are given in Table 2.
In our modeling, plume injection height is controlled by heat release, which is in turn driven by a combination of area burned, fuel loading, and combustion efficiency.The assumed diurnal time profile also affects the heat released, but this effect is constant among all fires.Because area burned and fuel loadings are highly variable, we examined the FEPS modeled plume injection height compared to the MISR observations as a function of modeled area burned and fuel loading.Estimates of area burned have a significant impact on the FEPS modeled plume injection height (Figure 7).Modeled plume heights for fires with an area burned of fewer than 800 hectares (1976 acres) are significantly lower than the MISR observed plume heights for the same fires.Likewise, modeled fires with an area burned greater than 2000 hectares (4940 acres) produced plume rises much greater than MISR observed data.Fuel loading is the total biomass available to burn per unit area of land.Fuel loading values vary greatly across the U.S. and sometimes within individual fires.The plume height model was examined to determine how it performs under different fuel loading scenarios.Figure 8 shows a notched box-whisker plot of the modeled (red) and observed (blue) plume height as a function of fuel loading.For fires with fuel loading values of about 2.25 kg of fuel/square meter (10 tons/acre) or less, the median observed plume height was statistically significantly higher than the median modeled plume height.For fires with fuel loading values greater than 2.25 kg/square meter, there was no statistically significant difference between the modeled and observed median plume height.

Comparison of CALIOP Observed Plume Heights to CMAQ Modeled Plume Heights
Table 3 shows the results of the statistical comparison between the CALIOP observed plume height and the CMAQ modeled plume height equivalent.The median CMAQ plume height equivalent was statistically significantly lower than the median CALIOP observed plume height.The range of plume height values was similar; however, the distribution of plume height values was significantly different (p = 0.005).A weak but positive correlation was observed between the CALIOP observed plume height and the CMAQ plume height equivalent (R 2 = 0.22).Figure 9 shows the results of the regression analysis of the CALIOP observations (x-axis) and the modeled plume height data (y-axis).While there was significant scatter in the results, most of the modeled results (>90%) are within a factor of two difference of the CALIOP observations.The spatial pattern of plume heights in the CALIOP-CMAQ data set is shown by different symbols in Figure 9 and on a map in Figure 10.There was no significant difference between the CALIOP and CMAQ data by region.Figure 10 shows a map of the fire locations and corresponding plume height data from the CALIOP observations (green bars) and the modeled plume heights (orange bars).Qualitatively, the CMAQ model compares favorably with the CALIOP measurements in some regions, such as central Idaho.Distribution statistics are not provided by region for the CALIOP analysis due to the small sample size.

Discussion and Conclusions
We compared plume injection height estimates from the default model path within the BlueSky Framework predictive modeling system to a MISR-derived database of plume top heights and CALIOP aerosol vertical profiles.MISR and CALIOP are unique in their ability to vertically resolve atmospheric aerosol signals globally and are thus ideally suited to observing wildland fire plume injection height and evaluating modeled plume injection heights.
On average, our modeled smoke plume tops were lower than the observed plume tops derived from both MISR and CALIOP.The modeled plume tops exhibited much greater variability than the MISR heights and there was a weak correlation between the two data sets.We also input the modeled plume tops into a chemical transport model and compared the resulting vertical profiles with CALIOP aerosol data.These data were weakly correlated and the modeled plume heights were lower than the CALIOP heights.
When examined regionally, the difference between observed and modeled plume tops of wild and prescribed fires that occurred in the southeastern U.S. was not significant; however, the observed and modeled plume heights in the Northwest were significantly different.In addition, under certain circumstances, modeled plume heights of fires in the Northwest were extremely high.These higher plume height differences tended to occur when the size of a fire, or the area burned, was very large.The regional differences in observed and modeled plume heights do not appear to be solely the result of vegetation and fuel loading.In the Southeast, fires tend to be relatively small, many of which are controlled or prescribed burns.Small and/or controlled burns tend to be less variable than large, uncontrolled fires.Because the factors that affect plume height tend to be less variable and more bounded, the plume height algorithms produce consistent results for fires in the Southeast.
Fire burn areas in the Northwest tend to be much larger, with lower fuel moisture and more diverse vegetation, or fuel loading.Fires with large burn areas and overall low fuel moistures produce significant heat, which greatly affects plume height in the plume injection height algorithm.Also, as large fires burn heterogeneously in multiple locations, the total heat released is divided into multiple plume cores, resulting in lower heat released per plume and lower plume injection heights [54].A key factor for large fires that is not considered in the assessed approach is the interaction of the plume with the ambient atmosphere.More sophisticated approaches are being applied in other modeling systems [55] and should be explored.However, the key variables necessary to drive plume injection height models, such as heat flux and the number of independent plume cores within a wild or prescribed fire, remain unavailable for predictive models.
The Briggs approach to plume height calculations depends on heat release, atmospheric stability, and wind speed.The use of default diurnal profiles of atmospheric stabilities and static wind speed in the calculation of plume injection heights is a weakness in the current modeling system.The plume height equivalent data from the model post-transport (CMAQ) showed less variability than the plume height data from the model pre-transport (FEPS) when compared to remotely sensed plume heights.In many cases plume heights were modeled well above the atmospheric boundary layer.While some smoke plumes are injected into the free troposphere in reality [56], most smoke within our modeling domain remains within the boundary layer.val Martin, et al. [20] estimated that 4-12% of plumes from fires are injected above the boundary layer in North America, although that includes the more energetic boreal fires of Canada and Alaska.
There is a lack of available data for the parameters that drive heat-release and plume-rise estimates.Existing plume height algorithms and dispersion models attempt to calculate plume injection height based on a series of physical processes to estimate heat emissions estimates, and some include atmospheric processes, but they are hampered by both a lack of available inputs and limited availability of observed data to test process-based algorithms.The Briggs equations underlying the modeled heights used here are commonly used in many different dispersion and air quality modeling systems.To improve this method for wildland fires, the inclusion of dynamic wind speed and atmospheric stability should be explored.In the absence of dynamic meteorology, average, region-specific plume injection heights may be a more representative technique, particularly in the case of large wildfires that burn across large areas of heterogeneous terrain.As new methods and approaches are developed, MISR and CALIPSO can both provide useful comparison data sets.

Figure 1 .
Figure 1.Example of the method for identifying days and times when the CALIPSO orbit path intersected an HMS smoke plume.The HMS smoke plume data are indicated by blue shading, the CALIPSO orbit paths are shown as green lines, and red circles are the SMARTFIRE fire locations.Red lines indicate times when the CALIPSO orbit path intersected the HMS smoke plume.

Figure 4 .
Figure 4. Conceptual illustration of the methodology used to derive the CMAQ plume top equivalent.

Figure 5 .
Figure 5. Regression analysis results for modeled plume heights (y-axis) versus the MISR observed plume heights (x-axis).Data points are color coded by region as shown on Figure 6.Also shown on the plot are the linear regression (solid line) and 1:2, 1:1, and 2:1 lines (dotted black lines).Note that a single value is off scale and not shown (MISR height = 2312 m; Modeled height = 18,699 m).

Figure 6 .
Figure 6.Map of the fire locations and corresponding plume height data from the modeled (red bars) and satellite derived (blue bars) data.Regions are color coded to match data points shown in Figure 5.

Figure 7 .
Figure 7. Notched box-whisker plot of modeled (red) and MISR observed (blue) plume heights as a function of modeled area burned.The notch is centered on the median concentration, widening to the size of the box to illustrate a 95% confidence interval in the median concentration value.The edges of the box illustrate 25th and 75th percentile concentrations.The whiskers indicate the lowest and highest values that are within 1.5 times the interquartile range (IQR).Outliers (denoted by plus signs) fall between 1.5 and 3 times the IQR.

Figure 8 .
Figure 8. Notched box-whisker plot of the modeled (red) and MISR observed (blue) plume heights as a function of modeled fuel loading.

Figure 9 .
Figure 9. Results of the regression analysis of the CALIOP observations (x-axis) and the modeled plume height data (y-axis).Data points are color coded by region as shown on Figure 10.Also shown on the plot are the linear regression (solid line) and 1:2, 1:1, and 2:1 lines (dotted black lines).

Figure 10 .
Figure 10.Map of the fire locations used for the comparison of CALIOP observed plume heights (blue bars) to modeled plume heights (red bars).Regions are color coded to match data points shown in Figure 9.

Table 1 .
Comparison of the overall distributions of observed median wind corrected plume top heights from the MISR instrument and modeled plume top heights from FEPS.Plume height units are meters (m) above ground level.

Table 2 .
Comparison of the overall distributions of observed median wind corrected plume top heights from the MISR instrument and modeled plume top heights from FEPS by region.Plume height units are meters (m) above ground level.Region definitions are shown in Figure6.

Table 3 .
Comparison of the overall distributions of observed plume top heights from the CALIOP instrument and modeled plume heights from CMAQ as in Table1.Plume height units are meters above ground level.