Next Article in Journal
Sources Profiles of Volatile Organic Compounds (VOCs) Measured in a Typical Industrial Process in Wuhan, Central China
Next Article in Special Issue
The Generation and Forecast of Extreme Winds during the Origin and Progression of the 2017 Tubbs Fire
Previous Article in Journal
Lightning Characteristics and Electric Charge Structure of a Hail-Producing Thunderstorm on the Eastern Qinghai–Tibetan Plateau
Previous Article in Special Issue
Development and Application of a Hot-Dry-Windy Index (HDW) Climatology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental Design of a Prescribed Burn Instrumentation

1
Department of Atmospheric Sciences, University of Utah, Salt Lake City, UT 84112-0110, USA
2
Department of Mathematical and Statistical Sciences, University of Colorado Denver, Denver, CO 80217-3364, USA
*
Author to whom correspondence should be addressed.
Atmosphere 2018, 9(8), 296; https://doi.org/10.3390/atmos9080296
Submission received: 16 June 2018 / Revised: 22 July 2018 / Accepted: 26 July 2018 / Published: 29 July 2018
(This article belongs to the Special Issue Fire and the Atmosphere)

Abstract

:
Observational data collected during experiments, such as the planned Fire and Smoke Model Evaluation Experiment (FASMEE), are critical for evaluating and transitioning coupled fire-atmosphere models like WRF-SFIRE and WRF-SFIRE-CHEM into operational use. Historical meteorological data, representing typical weather conditions for the anticipated burn locations and times, have been processed to initialize and run a set of simulations representing the planned experimental burns. Based on an analysis of these numerical simulations, this paper provides recommendations on the experimental setup such as size and duration of the burns, and optimal sensor placement. New techniques are developed to initialize coupled fire-atmosphere simulations with weather conditions typical of the planned burn locations and times. The variation and sensitivity analysis of the simulation design to model parameters performed by repeated Latin Hypercube Sampling is used to assess the locations of the sensors. The simulations provide the locations for the measurements that maximize the expected variation of the sensor outputs with varying the model parameters.

1. Introduction

Recent advancements in fire-atmosphere numerical modeling have increased the number of physical processes integrated into these coupled models. This greater complexity allows for a more comprehensive representation of the coupled interactions and feedbacks between the fire and the atmosphere. However, as a consequence of these advancements, data requirements for model initialization and validation have increased as well. As coupled fire-atmosphere models utilize local flow properties to parameterize fire progression, emissions, and plume rise, e.g., [1], integrated in situ measurements are needed for model validation and future development.
Fire and smoke models are being increasingly relied upon for wildland fire decision making and planning. However, many models are used without adequate validation and evaluation due to the lack of suitable data. Accurate estimates of smoke emissions from wildland fires are highly dependent on reliable characterizations of the area burned, pre-burn biomass of vegetation, fuels, and fuel consumption [2,3,4,5]. Smoke dispersion, on the other hand, strongly depends on the large-scale meteorology and plume injection height [6,7], which is affected by the fire geometry and intensity [8] impacted by local (fire-affected) circulations. Therefore, characterizing fire-atmosphere interactions, including wildland fire behavior and plume dynamics, is fundamental for improving estimates of smoke production and dispersion. The proposed observational data to be collected during the Fire and Smoke Model Evaluation Experiment (FASMEE) [9,10] are critical for evaluating and transitioning coupled fire-atmosphere models like WRF-SFIRE and WRF-SFIRE-CHEM [1,11,12] into operational use.
As the resolution of weather forecasting models increases, the fire-atmosphere interactions and smoke plumes (that had to be treated as sub-grid-scale processes requiring a simplified treatment through external parameterizations) become explicitly res/olvable. Thanks to improved computational capabilities, operational applications of coupled-fire atmosphere models e.g., Israeli Matash system based on WRF-SFIRE, become increasingly feasible [1,13]. The rapid increase in the resolution of numerical weather prediction products in recent years opens new avenues for development in the near-to-moderate future of integrated systems e.g., WRF-SFIRE-CHEM [12], capable of resolving, in a fully-coupled way, fire progression, plume rise, smoke dispersion and chemical transformations.
To advance current operational modeling capabilities, rigorous testing, evaluation and model improvements are needed [14,15]. This would allow for assessing model performance, uncertainties and the key sources of these uncertainties that need improvements.
Currently, there are insufficient observational data available to facilitate this work, especially in the context of coupled fire-atmosphere models requiring integrated datasets to provide information about the fuel, energy release, local micrometeorology, plume dynamics, and chemistry. Without such observational data, the accuracy and validity of even the most sophisticated physics-based modeling approaches cannot be fully assessed. Further, known gaps in scientific understanding, such as relationships between heat release, wind fields, plume dynamics and feedbacks to wildland fire behavior, cannot be addressed. This paper represents an attempt to: (1) identify the most critical observational needs; (2) assist in the planning of actual experimental burns; and (3) help optimize sensor placement during experimental burns. Additional details for the methods employed in this work are explained in [16].

2. Methods

2.1. Coupled Atmosphere-Fire Simulation

This study utilizes a suite of numerical simulations, performed using the coupled fire-atmosphere model WRF-SFIRE, which was designed to assist with the planning of the experimental phase of FASMEE. These are simulations of experimental burns to be conducted at the Fishlake, North Kaibab and Fort Stewart sites. Since numerical weather prediction models are unable to forecast weather conditions for burns planned years in advance, we opted to run simulations of experimental burns driven by weather conditions typical of the burn sites during the FASMEE burning seasons. The simulations are driven by time-varying, 3-dimensional historical reanalysis data representing the most typical weather conditions that meet burn requirements and account for the possible local diurnal variability in weather conditions. The most typical days are defined by a statistical analysis of historical weather observations from remote automated weather stations closest to the burn sites for days that meet the burn criteria. The statistical methodology used to select the most typical days is described in Section 2.2. These days serve as proxies for weather conditions driving the simulations of the experimental burns.
To support simulations of prescribed burns, models need to accommodate for arbitrary ignition lines that are either continuous (drop-torch ignitions) or consist of a discrete succession of drops of incendiary devices, dispensed from an all-terrain vehicle or an aircraft. The incendiary drops are at arbitrary locations and times, and the distances between drops can be smaller or larger than the fire simulation mesh size, which presents additional modeling challenges. WRF-SFIRE simulates continuous drop-torch ignition at steady speeds along straight lines, developed for the simulation of the FireFlux experiment [17]. A generalization of this scheme to complex ignition patterns is unwieldy, so, in this project, we have developed a new ignition mechanism by building on the earlier work on perimeter ignitions [18,19]: instead of modeling the fire spread until the time of the perimeter is reached, perimeter ignition in WRF-SFIRE prescribes the fire arrival time so that an appropriate atmospheric circulation due to the fire forcing can develop. The new scheme allows the user to specify the latest fire arrival time at the nodes (grid points) of the fire simulation mesh as an additional input array. This enables an integration of the natural fire progression simulated by the model with the ignitions specified by the fire arrival time.
In WRF-SFIRE, the state of the fire is represented by a so-called level set function, defined in the simulation domain. The fire is burning where the level set function is negative, and an ignition at a given location is implemented by making the level set function take a negative value at that location. However, the level set function is only defined by its values on the nodes of the fire simulation mesh, which are typically spaced several meters to tens of meters apart, while the incendiary device drops are at arbitrary locations. To achieve a sub-fire-mesh representation of the ignition process, the code takes the minimum of the level set function and cones (pointing down) with vertices at the individual ignition points. This approach allows the ignition locations and times not to coincide with the mesh points or the time steps in the model. The aperture of the cones is given by the imposed spread rate in the early stages of the ignition before the fire reaches the mesh node.
Numerical simulations for the planned experimental burns were performed for the typical days (meeting burn requirements for the specific sites) obtained from the analysis described in Section 2.2, in particular, Figure 1. The numerical experiments were conducted for Fishlake, North Kaibab, and Fort Stewart burn sites, and all used a multiscale setup of telescopic domains of grid scales gradually decreasing from 12 km to 148 m. Each simulation was performed with 30 m-resolution topography, land use and fuel maps from LANDFIRE, and driven by time-varying large-scale weather forcing obtained from Northern American Regional Reanalysis (NARR, [20]). Model configurations used for these runs are presented in Table 1. The domain configurations are shown in Figure 2.
The Fishlake burn has been started as specified in the burn plan by a line ignition following the local mountain crest (see Figure 2a). For North Kaibab burn, a point ignition was used, as the ignition procedure wasn’t specified for this site at the time of our simulations. The Fort Stewart simulations were executed with point-, single-line, and multi-line ignitions in order to cover a range of possible ignition scenarios.
We have used the WRFx system [21] to conveniently set up the simulations. WRFx can download the selected atmospheric product data for the initial and boundary conditions, process the USGS topography data and LANFDFIRE fuel data for the simulation domain, set up, execute, and monitor the WRF-SFIRE simulation, and postprocess the output into geolocated images, and stream the results to a visualization server.
WRF-SFIRE [22,23] evolved from the Coupled Atmosphere-Wildland Fire Environment (CAWFE) [24,25] and it has been a part of WRF since release 3.2 as WRF-Fire [23,26]. After release 3.3, WRF-SFIRE has continued to develop outside of WRF release [1,12]. The version in WRF release was recently selected as a foundation of the operational Colorado Fire Prediction System (CO-FPS) [27] and the level set method was improved [28].

2.2. Statistical Analysis of Typical Burn Dates

Coupled weather-fire-chemistry simulations (Section 2.1) provide important technological support for burn planning, operation, and analysis. To be relevant, these simulations are initialized with weather conditions typical of the burn location. In this analysis, “typical” is defined by a statistical analysis procedure applied to weather-station data records.
To start, consider a weather state defined by a single quantity e.g., air temperature. Then, one might define a “typical” day as one that has temperature closest to the historical sample mean of the temperature from all days considered.
More generally, a weather state is defined as a vector of temperature and multiple additional quantities, namely relative humidity, and the speed, direction, and gust of the horizontal wind vector. Such data present a challenge, because the quantities have different units, value ranges, and statistical character, which make them incommensurate. As an example, the input data (meteorological states comprising temperature, relative humidity, wind speed, direction and gust) for the “typical day” analysis are displayed in Figure 3 for one of the stations. Such figures are used as a “sanity check”, i.e., a minimal quality-control of the data (indeed, some spurious values were discovered this way), as well as to visualize the burn requirements for state values, day-of-year and time-of-day, and the relative amounts of state values that do and do not meet the burn requirements. They also indicate that the variability of the raw data is far from that of a multivariate Gaussian random process.
For the analysis, the bounds on the value ranges are first removed by transforming from relative humidity, speed, direction and gust, to obtain the weather state column vector
x = T , T d , u , v , T ,
where T is the temperature, T d is the dew point, u and v are eastward and westward components of the wind vector, and is the logarithm of the wind gust. The transformed data and their statistics are presented in Figure 4. It is evident from the scatter plots that the transformed state is still not very impressively Gaussian, but certainly is much more Gaussian than the original (temperature, relative humidity, wind speed, direction and gust) description. Also, the ellipse tilts and major-minor axis ratios confirm the general impression of the scatter of data points that meet the burn requirements.
Even after this transformation, the weather-state vector entries still have different units, and their collection for all days considered varies more in some directions and less in others, so that measuring closeness by a simple root of sum of squared differences is not appropriate. Instead, we choose as typical days the days with the least Mahalanobis distance from the sample mean. The sample mean vector is computed as
x = 1 N t = 1 N x t ,
where the sum is taken over all measurement times t when x t satisfies the burn criteria (Table 2) in the considered years regardless of burn seasons. The sample covariance matrix is then computed by
C = x x x x T
and the Mahalanobis distance of a vector x t from the sample mean is defined by
δ t = C 1 / 2 x t x = x t x T C 1 x t x ,
where the Cholesky factorization of C and the Mahalanobis deviation vector at time t are
C 1 / 2 C 1 / 2 T = C , δ t = C 1 / 2 x t x .
Please note that by definition, the δ t have zero sample mean and their sample covariance is identity, so the variance
var δ = δ T δ = δ 2 = 5 ,
which is useful to interpret the values of δ t .
The Mahalanobis distance puts more weight on the differences between weather state vectors along the directions where the collection of the weather state vectors varies less, on average. Another interpretation of this method is as a maximum likelihood: if the weather state vectors were a random sample from a multivariate Gaussian distribution, the days with smaller Mahalanobis distances from the sample mean would be the more likely ones. Finally, this analysis does not directly consider steadiness i.e., what happens before or after the time t that minimizes δ t . In theory, and to the extent the statistics are truly Gaussian, δ t 2 might be interpreted as a quadratic potential (like a harmonic oscillator has), so typical x t should be near equilibrium.
Software written in MATLAB (MathWorks, Natick, MA, USA) which implements the determination of typical days is provided in Supplementary Materials.

2.3. Statistical Optimization of Sensor Placement

It is most useful to place (the limited number of) sensors where the observations or parameters of interest change or show the largest variance. To estimate the sensitivity of measurements to model parameters, a global sensitivity analysis with repeated Latin Hypercube Sampling [29,30,31] was used. Consider a simple precursor method: run several simulations with the values of the parameters of interest chosen randomly and then compute the sample variance of the output over the simulations at various potential sensor locations. The method utilized here improves on such a conventional approach in two aspects: repeated Latin Hypercube Sampling (rLHS) allows for a drastic reduction in the number of simulations to be made by a systematic coverage of the parameter space, and the Sobol variance decomposition [32] can distinguish the effect of individual parameters. This method was originally developed to evaluate computer simulations at the Los Alamos National Laboratory [29,30] and was recently used for guiding sensor placement for oil reservoir engineering [33].

2.3.1. Repeated Latin Hypercube Sampling and Sobol Variance Decomposition

Consider a model with L parameters, divide the range of each parameter into N intervals of equal probability, and for each interval choose the probability midpoint as a sampling point 1 , , N . Then set up N parameter vectors by choosing the sampling points of each of the L parameter randomly, without repeating the same sampling point in any of the parameters. The simulation is run on each parameter vector. See Figure 5 for an example.
Repeated Latin Hypercube Sampling (rLHS) consists of r repeats of LHS, with new random permutations of the choices of the sampling points for each parameter. The complete rLHS scheme requires r N simulations, with each simulation receiving a vector of L parameters.
Each of the r N simulations delivers its output Y. Now fix the value of parameter X i at one of the sampling points x i j and take the average of the outputs Y over all simulations where parameter had that value, X i = x i j . The result is a function of the value of the parameter X i . This function is called the conditional expectation of Y given X i , and we denote it by Y i . Since the parameter X i attains N sampling values with equal probability 1 / N , the conditional mean Y i also attains N possible values with equal probability 1 / N each. From those values of Y i , we compute its sample variance var Y i in the usual way except that we use the denominator N instead of N + 1 in the expression for the variance. We also compute the sample variance var Y over all r N values, again using r N in the denominator instead of the usual r N + 1 .
The importance of this procedure lies in the following interpretation:
  • var Y is the total variability of the output Y over the range of the parameters X 1 , , X L . We take this as the strength of the signal from the sensor.
  • var Y i , called the Variance of the Conditional Expectation (VCE), is an estimate of the variability of the output Y due to the parameter X i .
  • The ratio
    eff Y i = var Y i / var Y
    is an estimate of the relative effect of parameter X i on the output Y. It is called the correlation ratio [30] or the sensitivity index [31] of the parameter X i .
This interpretation is motivated by the case when the model is linear,
Y = a 1 X 1 + + a L X L ,
where the coefficients a 1 , , a L are constant and the inputs X 1 , , X L are independent random variables with unit variance. Then, it can be shown that under suitable assumptions (which are certainly satisfied e.g., when X i have Gaussian distribution), the VCE from the rLHS procedure with Sobol’s variance decomposition provides asymptotically exact approximation of the VCE of Y given X i for large number of samples and intervals:
var Y i a i 2
as the number of repetitions r and the number of intervals N . That is, the method recovers asymptotically the variance of the output Y due to one parameter X i from the signal Y that combines the effect of many independent parameters, by probing the model. If the parameters are not independent, or the model is not linear, the recovery is only approximate.
MATLAB software which implements the rLHS and validates the implementation by observing the convergence in (1) computationally is provided in Supplementary Materials.
In the present application, Y is defined as one of measured values at some location e.g., the vertical wind velocity at a fixed time, longitude, latitude, and height. If the effect of one particular parameter is of interest, it makes sense to look for locations where the sensitivity index of that parameter is large. Then the signal, due to that parameter, stands out from the rest.
In the case where there is a small number of repetitions and intervals, the convergence of the sample variance of the measured value Y as well as the convergence in (1) may be far from the asymptotic range. However, since the measured values Y are taken from the same set of simulations and, in particular, when the measured values Y are the same measured quantity at nearby locations, their sampling errors are correlated and comparing their sample variances as in the figures in Section 3.2 is justified.
Finally, since the quantities of interest will be associated with a fixed location while the fire is moving, there is only a significant signal when the fire or smoke is advected over the sensor location. Hence, the sums of time series of variances associated with a fixed location, which represents the total variability captured by the sensor over time are also of an interest.

3. Results

3.1. Determination of Typical Burn Days

From Section 2.2, the Mahalanobis distance δ t of the weather vector at a station from the long-term mean is dimensionless, its mean value squared is δ 2 = 5 , and the most “typical” day is defined as one with the smallest δ t in the given burn time window. The resulting identification of typical burn days is in Figure 1 and Table 3. See Appendix A for further details.

3.2. Results from Numerical Experiments

3.2.1. Simulations of Planned Experimental Burns in Fishlake

The analysis of the expected plume characteristics has been performed for the Fishlake simulations. Three simulations were run for the day of 09/03/2014, (which has been identified based on the weather observations from the FSHU1 weather station as the most typical day meeting the burn criteria defined in Table 2):
  • reference simulation without fire,
  • standard simulation with fire, and
  • simulation with fire and doubled heat flux from the fire to the atmosphere.
The first run serves as a reference point for assessing the impact of the fire on the boundary layer when compared to the second run. The third run is intended to account for possible underestimation of the fire heat fluxes under high fuel load conditions characteristic for the Fishlake burn site, as well as to assess sensitivity of the simulated plume height and vertical velocities to the fire heat flux.
As shown in Figure 6, for the standard and doubled heat flux fire simulations, maximum updraft velocities tended to spike during the ignition and then gradually rise over the next 4 h. The maximum updraft velocities reach 5 m/s in the standard simulation and twice as much in the simulation with the doubled heat flux. The elevations where maximum updrafts occur are significantly less sensitive to the fire heat flux enhancement. In both cases, they vary between 170 m and 1000 m. These simulations suggest that from the standpoint of the plume characterization it is important that the experimental burn lasts long enough so that the smoke column fully evolves (at least a couple of hours). The vertical in-plume velocity measurements (LIDAR vertical scans, UAV platforms) should allow for plume scanning up to at least 1.5 km and be able to withstand vertical velocities up to 10 m/s. To provide a guideline for the in-plume measurements of chemical species, the simulations were executed with fire smoke represented as a passive tracer. Figure 7a shows vertical profiles of the smoke concentration normalized by its maximum value for the baseline simulation and the case with doubled heat flux. Both cases show similar concentration profiles with a maximum near the ground and a single peak located around 800 m above the ground level. In the double-heat case smoke penetrates into higher elevations. The time series in Figure 7b show a rapid plume rise during the ignition, followed by a 2 h period of constant plume top height, followed by a sudden drop and a gradual increase with some further fluctuations. As the heat released during the ignition procedure is also doubled in the double heat flux case, the differences in the plume top heights between the cases are most evident during the first 2 h after ignition. Later, as the plume fully develops, the plume heights in these two simulations become similar (within 250 m). The concentration plots presented in Figure 7a suggest that plume sampling should be focused on the layer 3000 m AGL (Above Ground Level), with a particular focus on the 500 m and 1500 m AGL, where the simulations tend to produce the highest tracer concentrations. It has to be noted that as the plume evolves in response to the changing atmospheric conditions and the fire behavior the elevation of the plume layer may change rapidly.
To assess the fire’s impact on boundary layer flow in the baseline and double heat cases, the horizontal fire-induced winds are computed as a difference between the U and V wind speed components from the fire runs and the reference no-fire run. The time series of the levels where the maximum fire-induced winds occur and their magnitudes for both fire simulations are presented in Figure 8. The elevations of the fire-induced wind maxima are similar in both simulations, suggesting that the level at which the fire accelerates flow most rapidly is not particularly sensitive to the fire intensity (heat flux). However, the magnitudes of the fire-induced wind are. The maximum winds in the baseline simulation reach 7 m/s, while in the double heat flux case they reach 9 m/s. The simulations indicate that the maximum wind flow is either close to the ground (15–25 m AGL) or aloft, mostly between 800 m and 1200 m AGL. Therefore, from the measurement standpoint, the optimal sampling strategy would contain a combination of a meterological tower sampling the near-surface flow and Doppler LiDAR(s) providing wind speed measurements at elevations at least up to 1500 m AGL.

3.2.2. Simulations of Planned Experimental Burns in North Kaibab

To provide insight into the second western burn at the North Kaibab site, a simulation was performed for day 09/05/2008, identified based on the weather observations from the QLBA3 weather station as the most typical day meeting the burn criteria defined in Table 2. In the North Kaibab case, Figure 9a shows that the smoke layer is more uniformly distributed across the first 1500 m than in the Fishlake case, with no evident concentration peak. Figure 9b shows lower simulated plume top heights than in the Fishlake case. They increase gradually over the first 4 h and stabilize at around 2000 m, 200–700 m lower than in the Fishlake runs.
The evolution of the fire-induced flow is presented in Figure 10. The fire-induced winds, in this case, are similar in magnitude to the Fishlake baseline case, but their time evolution is different. Contrary to the winds from the Fishlake simulation showing a strong initial spike in fire-induced winds followed by a gradual increase, the North Kaibab winds ramp up to about 4 m/s, reaching a maximum of about 7 m/s at 2.5 h into the simulation, and then keep decreasing. The timing of the peak surface winds corresponds to the time when the maximum updraft is located relatively close to the surface (below 500 m AGL). The convective updraft evolves gradually in time, reaching maximum values close to 6 m/s after 5 h since ignition. Please note that as a result of a point ignition (less intense than the line ignition in the Fishlake burn), this evolution is not disturbed by the initial ignition-induced spike evident in the Fishlake run (Figure 8). The updraft analysis shown in Figure 10b indicates that over the first one-and-a-half h, the maximum updrafts are located at 2500–3000 m AGL (probably as an effect of fire-induced disturbance inducing vertical mixing at high elevations). As the fire evolves, updrafts in the plume core intensify and become dominant, which is indicated by the elevations of the maximum updrafts going back to 500–1000 m AGL.

3.2.3. Simulations of Planned Experimental Burns in Fort Stewart

Numerical simulations for Fort Stewart were performed with the main intent of providing an insight into the impact of the ignition procedure on plume evolution. For prescribed burns that require complex aerial ignitions, multiple ignition lines are formed in a relatively short period of time, with significant heat fluxes coming directly from the igniting agents like the gasoline mixture in Helitorch or ping-pong balls (DAID) systems. As the aerial ignition procedure is generally fast and difficult to precisely capture by scanning IR systems which do not provide a continuous static field of view, it is important to know how important the ignition itself is for further plume evolution. To assess that, five different ignition procedures were simulated: a single ignition point, a single ignition line of two different thicknesses, a set of three parallel lines, and a set of five parallel lines. All ignition lines were oriented approximately from northwest to southeast (perpendicular to the mean wind direction). They were 2 m thick, except for the single bold ignition line which was 30 m thick. The lateral rate of spread within the ignition lines was set to 0.5 m/s in order to assure gradual initial heat release. Each line was ignited to its full length in 200 s. In the multiple-line cases, ignitions were performed sequentially, and the spacing between the lines was 250 m. The results from these simulations are presented in Figure 11. The time series of the maximum vertical velocity and the plume height indicate that the ignition itself plays a significant role in updraft evolution, especially during the initial phase of the burn. The time series from different ignition procedures only start to converge 7 h after ignition. These results suggest that for relatively short experimental fires (expected to last a couple of hours), accurate documentation of the ignition procedure is critical for a realistic representation of the initial plume evolution in subsequent numerical simulations.
It should be also noted that the initial interactions between ignition lines are hard to account for in numerical experiments performed using models with parameterized fire spread. In these models the ignition process itself is crudely simplified, limiting their capabilities in terms of rendering these effects. Therefore, simple ignition patterns are desired in experimental burns so that the fire can be realistically initialized in the model.
Aside from the simulations targeting the ignition effects, additional simulations were performed for 04/22/2014 (a typical day based on statistical analysis of observations from the KCWV weather station (Figure 2c). Results from this numerical experiment are presented in Figure 12. The Fort Stewart site is significantly different from south-western Fishlake and North Kaibab sites in terms of both fuel and topography. The latter ones are located in areas of complex terrain where steep elevation strongly impacts both vertical and horizontal flows, making fire-induced flows harder to detect. It seems that in complex terrain the destabilization induced by the fire can impact the atmospheric state at significant distances downwind from the plume core. Strong upslope winds can induce updrafts of a similar magnitude to the fire itself, and the interaction between the flow and the mountain ridges may lead to significant variations in the horizontal winds that become visible only when differences between fire and no-fire cases are computed. As a consequence, the heights of maximum fire-induced winds and vertical velocities fluctuate strongly in the North Kaibab simulations but are relatively steady in the Fort Stewart burn. In Fort Stewart for instance, the elevations of strongest fire-induced winds are confined to a shallow layer between 0 and 140 m AGL, with the majority of strongest fire winds occurring within the 20 m layer above the ground (Figure 12a). These results suggest that a combination of meteorological towers and LIDARs located close to fire line should be able to sample most of the maxima in the fire-induced flow. In the Fort Stewart simulation, both fire-induced winds and updraft velocities evolve very quickly, reaching a quasi-steady state just after one hour from the ignition. These simulations suggest, that despite the fuel load and expected fire intensity are lower for Fort Stewart than for Fishlake burns, in the context of characterization of the fire-induced circulation, the Fort Stewart site has a potential to provide better and easier to interpret results relative to the more complex southwestern sites.

3.3. Sensitivity Study of Sensor Placement

We have performed a sensitivity study following the rLHS with variance decomposition methodology, described in Section 2.3.1. We have chosen the Fishlake simulation standard case from Section 3.2.1 as the base case.

3.3.1. Sampling Setup

We have chosen L = 7 simulation parameters, which were:
  • 10 h fuel moisture content, varying from 0.04 to 0.14 water mass/dry fuel mass. 1 h fuel moisture was 10 h fuel moisture minus 0.01, 100 h fuel moisture was 10 h fuel moisture plus 0.01, and live fuel moisture was 0.78. These values were entered as initial moisture values and did not change with time. See [1,19] for a further description of the fuel moisture in WRF-SFIRE.
  • Heat extinction depth, varying from 6 m to 50 m. The fire heat flux is entered into WRF boundary layer with exponential decay, rather than all into the bottom layer of cells. The heat is apportioned depending on the height of the cell center above the terrain, with weight 1 at the ground and at the heat extinction depth. This gradual heat insertion is a parameterization for unresolved mixing and radiative heat transfer.
  • Heat flux multiplier, varying from 0.5 to 2. The heat flux multiplier was chosen as a measure of the fire effect on the atmosphere; unlike the fuel load, it does not influence the rate of spread.
  • Multiplier for the omnidirectional component R in the approximate fire rate of spread (ROS) formula following [34],
    R O S = R + R W + R S
  • Multiplier for the wind-induced ROS component R W
  • Multiplier for the slope-induced ROS component R S
  • The simulation day, selected from the “typical” burn days
The first 6 parameters were constant over the whole fire simulation domain. Since the 6 parameters are positive, first they were transformed by taking the logarithm and then a Gaussian distribution centered within the transformed interval was determined so that 90% probability was between the lowest and the highest value. The transformed interval was then divided into subintervals of equal probability, sampling points chosen as the middle probability values, and transformed back into the original physical space. The resulting sampling values are in Table 4.
Total r = 200 repetition was done, resulting in r N = 1000 simulations. The selection of the sample points in the first replicant is in Figure 5. For the purpose of the sensitivity analysis, the output of the model is the value of a quantity of interest at a location in the physical space (a node of the simulation mesh on the Earth’s surface, or interpolated to a specified height). We then visualize the values associated with the quantity of interest, such as the VCE (Variance of the Conditional Expectation, an estimate of the variability of the output due to the parameter change, cf., Section 2.3.1) as a function of location on the Earth surface. Please note that in WRF, the vertical position of mesh nodes is not fixed, rather it is derived from the geopotential height, which is a part of the solution, and it changes with time.
We consider the following quantities of interest:
  • The vertical velocity vector component W, interpolated to a given height above the terrain, or to a given altitude above the sea level.
  • The smoke intensity (the concentration of WRF tracer tr17_1), interpolated to a given height above the terrain, or to a given altitude above the sea level.
  • Plume-top height, derived from the smoke concentration.
Because some extreme values of the sampling parameters caused instability in the WRF surface scheme during and immediately after ignition (when the updraft carrying the heat away is not yet established), we have simulated the initial 30 min when the ignition occurs using the standard case and only then switched to the modified parameter values. We show the results from simulations running for an additional 7.5 h, during which the model state was captured and analyzed at hourly intervals. The resultant 7 time frames were averaged to provide a picture corresponding to the 7 h of the burn.

3.4. Computational Results for the Fishlake Burn Simulation

In the following section selected results of the repeated Latin Hypercube Sampling (rLHS) and the analysis of variance described above are shown. In all computations reported here, the average variance is computed as variance for each simulation day separately, and then averaged over the simulation days.
The vertical level of highest variances in the updraft velocities (1200 m above the ground) have been selected for visualization based on the vertical velocity variance time series plotted in the left panel of Figure 13. Based on a similar time series of the smoke concentrations and the plume top height estimates form the baseline case (Figure 6), the analysis level for smoke has been set to 1400 m. Figure 13b shows the variance of the vertical wind velocity at 1200 m above the terrain, computed from 25 rLHS-sampled simulations. Two well-defined maxima are located outside of the northwestern plot boundary, indicating optimal locations for upper-level vertical velocity measurements. Local variance reaches up to 3 m 2 / s 2 , which corresponds to local deviations in the vertical velocity of about 1.7 m/s between the sampled simulations, which is high enough to be sampled by a Doppler LiDAR, Radar or Sodar. However, the fact that the updraft variances are concentrated over small regions indicate that optimal placement of the measurement devices may play a critical role in constraining model parameters based on the vertical velocity observations.
Figure 14 shows maps of the mean and the normalized variance of the plume top height. The former informs about the expected vertical plume extent which should be taken into account when the plume sampling strategies are considered. The latter shows regions of most pronounced impact of the model parameters on the plume rise. Interestingly, the region directly downwind from the fire, where the plume reaches highest elevations, is not the same as the region of highest plume height average variance. The plume height average variance exhibits high values over two symmetrical regions located on both sides of the plume axis and are pushed much more downwind from the fire compared to the mean plume top height. These results indicate that while the localized vertical smoke sampling (corkscrew path) may be optimal for assessing the maximum plume top height, the sampling region must be significantly extended downwind in order to adequately sample the plume top height variance.
Figure 15a shows the map of the mean smoke concentration at 1400 m, and Figure 15b is the map of the smoke concentration variance. The mean smoke concentration is intended to inform the measurements where the probability of successful smoke sampling is the highest. The variance on the other hand, shows where the tested model parameters have strongest impact on the smoke concentrations. It is noteworthy that high values of smoke concentrations variances are found over much larger region than the vertical velocity variance shown in Figure 13b, which suggest that vertical velocity measurements are more suitable for local platforms like ground LiDARs, while the smoke measurements could be performed from airborne platforms covering larger areas.
The rLHS analysis not only informs where the measurement should be taken in order to constrain model parameters, but also allows to find relative contribution of the tested parameters to variances in the variables of interest. We present results of the first-order variance decomposition of vertical velocity at 1200 m, smoke concentration at 1400 m, and the plume top height, attributed to the most important simulation parameters. These spatial plots of the sensitivity indices inform about the relative importance of the tested parameters (such as fuel moisture, heat extinction depth, fire heat flux, etc.) as well as indicate regions of highest impact of these parameters on the variables of interest.
The most critical parameters controlling the vertical velocity, smoke concentration and the plume top height are shown in Figure 16. The heat flux multiplier contributes to the variance of vertical velocity, smoke concentration and the plume top height up to 50%. Unsurprisingly, the fire heat flux plays the major role in driving the plume rise, which indicates the importance of a detailed fire heat characterization during experimental burns. The second most significant contribution (up to 40%) comes from the heat extinction depth, defining the depth over which the fire heat flux is distributed vertically in the model. The extinction depth seems to have an impact of similar magnitude on the variables of interest like the heat flux multiplier, but its character is different. Figure 16 suggests that the heat flux multiplier has the opposite effect on the extinction depth. For instance, the vertical velocities are impacted by the heat flux multiplier mostly downwind (northeast) from the fire (see left top panel). The extinction depth, on the other hand, affects mostly vertical velocities upwind from the fire (see right top panel). Similarly, the smoke concentration is impacted by the heat flux in the north-eastern part of the domain, while the extinction depth seems to impact mostly the south-western part. A similar situation can be observed in the bottom panels of Figure 16, showing the impact of these parameters on the plume top height. Here, the heat flux impact is confined to the very core of the plume, while the extinction depth has the most pronounced impact on the sides of the plume. The differences in these patterns indicate that the placement of the sensor optimal form the perspective of constraining the heat flux may be not optimal for constraining the extinction depth and vice-versa. Also, the patterns of the heat flux contribution seem generally more organized than the ones of the extinction depth. The impact of the extinction depth seems to have non-local character especially when vertical velocity and smoke concentration is considered, spreading across large areas upwind from the fire. That suggests that estimating the extinction depth through the observations of vertical velocities and smoke concentration may be difficult, and may require direct tower-based observations of the vertical attenuation of the fire heat flux.
The analysis presented above focuses on the impact of model parameters on measurable variables like the vertical velocity or the plume height. In that sense, it informs where the signal coming from the changes in the parameters is expected to be highest and consequently, shows locations where measurements should be taken to detect it. As the weather conditions on the burn day cannot be predicated at the moment of experimental planning, the analysis presented in this study has been employed to find most typical historical days for all the burn stations. Five most typical days have been included in this analysis. On the top of the baseline simulation, performed for 09/03/2014, analogous runs have been executed for other 4 typical days listed in Table 3. These 5 days are treated as an ensemble of most typical scenarios for the Fishlake experimental burn. Shown below are the means and average variances of the variables of interest computed across the simulations performed for these days. Figure 17 shows maps of mean and variance of the smoke concentration at 1400 m. A pattern can be noticed there, indicating general north-northwest smoke dispersion. The fact that regions of maxima in both mean concentration and the variance can be easily identified, indicates that the typical days are similar to each other in terms of a general flow pattern in the Fishlake region. It should be noted that the wind direction was not a part of burn requirements for this burn site which defined only the wind speed, temperature, and relative humidity. Since typical days show relatively consistent flow pattern, it can be concluded that typical days in terms of the air humidity, wind speed, and moisture are generally associated with south-southwesterly flow.
The maps of the vertical velocity variance and the mean plume top height are shown in Figure 18. A clear maximum in the updraft fluctuations is evident north-west from the end of the ignition line, as well as a well-defined north-northeasterly plume reaching 3200 m. The statistical analysis of simulations performed for the most typical day such as the one shown in Figure 17 and Figure 18, are intended to be used as a guideline in the planning stages of experimental burns when optimal locations of smoke and vertical velocity measurements are being considered.

4. Conclusions

4.1. Estimation of Typical Days

If the meteorological state vector x had a multivariate Gaussian distribution, then the most typical day as defined in this work would also be the day with the most probable state. As Figure 1 shows, a very typical day need not occur very often, or within an interval of similarly typical days; rather, it is a day when the state is in a rigorous sense very close to its sample mean. Indeed, for a non-Gaussian distribution (the case here), a day with state close to its mean need not occur very often, even approximately. From this standpoint, the analysis performed here assures that statistically, any historical day meeting the burn requirements will be similar to the identified typical day, which fits the purpose of this work. In future work, an alternative to consider would be to pick a state (and thus, the day it occurs) randomly from a multivariate bin containing states that do occur simultaneously very often. However, before following this alternative, one would have to choose among many competing theories about how to choose the sizes of such a bin robustly.

4.2. Results from Numerical Experiments

The numerical experiments presented in this study were performed and analyzed to aid with the preparation of the FASMEE experimental plan. Since we do not know how similar the typical burn days are to actual burn days, the results presented here cannot provide deterministic forecasts. Thus, we can only provide examples of possible realizations of the planned experimental burns. A timely simulation (possibly with a truncated sensitivity study) can be run immediately before the burn to inform the deployment of observation resources, as it was successfully done for the FireFlux II experiment [35].
It has to be emphasized, that presented results of numerical simulations of anticipated burns should be treated as guidelines only. Even though WRF-SFIRE has been validated against vertical velocities measured during the FireFlux burn [17], and the simulated plume top heights compared favorably with MISR observations [12], its ability to realistically resolve the wildfire plume dynamics, and the near surface flow induced by a wildland fire remain to be assessed. Especially, in the context of simplified representation of the canopy flow in WRF, which may limit its ability to render near-surface winds in forested regions.
The results showing model sensitivity to grid resolution (not shown here) indicates that local updraft velocities may be significantly underestimated in the model simulations when the model mesh is not fine enough to resolve small-scale plume features. Therefore the simulated vertical velocities should be treated as lower limits for the expected values. The simulated updraft strengths suggest that measurement platforms used during the experiments should be robust enough to sustain updrafts exceeding 10 m/s. Our simulations suggest that the strongest updrafts will be observed in the Fishlake burn. However, due to the influence of complex local topography, exact identification of the fire impact on the local circulation may be difficult there. Therefore, the Fort Stewart burn performed in fat terrain would serve as a great complementary burn that could provide a clearer picture of fire-induced circulations.
One of the most important aspects of experimental burns is their durations. As indicated by the western U.S. simulations, the burn experiments should last at least a couple of hours in order to fully capture the fire plume evolution, and consequently, the burn plots should be big enough to enable fire progression over this period of time. For small experimental burns, the fire may reach the end of the burn plot before the plume fully evolves.
Unlike the fire-induced horizontal winds, expected to reach maximum speeds relatively close to the ground, maximum vertical velocities are expected at significantly higher levels. Simulations of the experimental burns indicate strong updrafts located up to a couple of kilometers aloft, so Doppler LiDARs and other platforms taking vertical wind measurements should provide adequate vertical coverage. Regarding smoke measurements, the simulations suggest that the vertical smoke concentration profiles and the plume top heights may vary significantly between the experimental burns. The Fishlake burn of expected highest intensity is associated with an elevated smoke concentration peak and higher plume tops than other lower-intensity burns. Therefore, it seems beneficial to adjust the measurement heights accordingly to assure optimal plume sampling.
Numerical simulations of the Fort Stewart burn indicate a strong sensitivity of the initial plume evolution to the fire ignition pattern. Therefore, from the standpoint of model validations, a detailed characterization of the ignition procedure is necessary. Also, as the fire initialization is crudely implemented in the model, interactions between simultaneously ignited fire lines are difficult to capture. Therefore, in experimental burns, simple ignitions are preferable, especially in cases when continuous IR observations at sufficient temporal and horizontal resolutions are not possible. Similarly, the sensitivity of the simulated plume top heights and vertical velocities to the fire heat flux (demonstrated based on the Fishlake simulations), indicates that detailed characterization of fire heat fluxes (radiative and convective) is critical to model evaluation. Convective heat fluxes from the fire must be observed during the experimental burns, in addition to radiative fluxes.

4.3. Sensitivity Study of Sensor Placement

The results presented in Section 3.3.1 illustrate how numerical modeling combined with advanced statistical methods can objectively guide planning of experimental burns. There are multiple possible applications of this method. If an experiment is performed with the intention to constrain certain model parameters, such analysis can help assess the feasibility of this approach, and design an optimal sampling strategy. In the analyzed Fishlake burn, for instance, the results indicate that assessing the extinction depth through the vertical and smoke concentration measurements may not be feasible, which leads to a recommendation of direct measurements of the fire heat flux at various heights. This method also allows identification of the most critical model parameters that control plume dynamics. Here the heat flux has been identified as the one of highest importance for the plume development, which emphasizes the value of accurate characterization of fire heat fluxes during the experimental burns. As this method utilizes the results of statistical analysis of typical burn days it also informs about expected local flow patterns, plume dispersion and fire behavior in a probabilistic sense. If the decisions about sensor placements have to be made in advance, such analysis may guide them. In case of the Fishlake burn this analysis indicates that generally the area north-northeast from the burn plot should be considered for placing the measurement platforms.
We have performed only a pilot analysis of variations for several sampled variants of the Fishlake experimental burn. Listed below are some preliminary observations and implications for future research:
  • The statistical analysis of runs with carefully sampled parameter sets was shown to provide clear guidelines on placing the measurements in space and time.
  • Also, some measurements were found more affected by certain parameters, which can inform what parameters can be indirectly constrained by observations.
  • Analyses of the variance in smoke concentration vertical velocity and plume top height from runs executed for different days and with different parameter sets show similar patterns indicating that typical day statistics managed to identify statistically similar days from the standpoint of plume rise and dispersion.
  • The variability of the plume at higher altitudes (here, 1000–1400 m above the ground) is concentrated up to couple kilometers downwind from the plot, which suggests that sampling right above the fire, optimal for the fire heat flux measurements, may not be optimal for sampling most active parts of the plume.
  • Additional parameters can be considered. The cost of one repetition does not increase with the number of parameters, but more repetitions need to be done for statistical convergence. The variability of the outputs will increase and the added parameters will help model additional uncertainty, which is always present in reality.

4.4. Computing Resources

The rLHS computations reported here consisted of 1000 simulations executed in about a week to reduce the sampling error. The wall clock time of one Fishlake simulation is about 2 h on 3600 Intel cores, and it produces about 15 GB of data. Counting repetition e.g., for runs failed for various reasons, for development, and for experimentations, about 1,000,000 core hours were used and 50 TB of disk space was needed to store the outputs. After cleaning and elimination of redundant files and keeping only the final version, the final data set was 12 TB. The statistical processing was done and the visualizations were generated on an auxiliary high-memory cluster in MATLAB. The analysis of the outputs of the simulations was performed in serial mode, which significantly contributed to the total analysis time. Overall, several thousands of hours on the auxiliary cluster were used.
As multiple simulations can be run concurrently, we suggest performing a similar analysis prior to the burn using available forecast data. With sufficient supercomputer allocation, this is feasible.

Supplementary Materials

The WRF-SFIRE and WRFx software (wrfxpy, wrfxctrl, wrxfweb) are available at https://github.com/openwfm. Selected visualizations can be accessed at http://demo.openwfm.org/fasmee. MATLAB software developed for the experimental design is available at https://github.com/openwfm/design.

Author Contributions

A.K. set up and analyzed the computational experiments and wrote the introduction, the computational experiments section, and the conclusion. A.F. devised the statistical analysis of typical burn days, wrote the statistical analysis section, and suggested the use of rLHS and Sobol variance decomposition for the sensitivity analysis. J.M. made the custom developments within the WRF-SFIRE and WRFx softwares to support this work, devised and performed the sensitivity analysis, wrote the sensitivity analysis section, and finalized the paper for submission.

Funding

This research was funded in part by the Joint Fire Science Program under grant JFSP 16-4-5-03 as a part of the Fire and Smoke Model Evaluation Experiment (FASMEE) Phase I, the National Science Foundation under grants DMS-1216481, ICER-1664175, and NASA under grant NNX13AH59G to Colorado State University, Sher Schranz, PI, CIRA in affiliation with NOAA/ESRL.

Acknowledgments

The authors would like to acknowledge high-performance computing support from Yellowstone (ark:/85065/d7wd3xhc) and Cheyenne (doi:10.5065/D6RX99HX) provided by NCAR’s Computational and Information Systems Laboratory, sponsored by the National Science Foundation. The computing support from the University of Utah Center for High Performance Computing and the Center for Computational Mathematics, University of Colorado Denver, is greatly appreciated. The authors would like to thank Derek Mallia for reading the paper and suggesting improvements, and anonymous reviewers for helpful comments, which resulted in improving the paper.

Conflicts of Interest

The authors declare no conflict of interest. The funding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

Appendix A. Statistics of Weather Data

This appendix contains a visualization of the remaining weather station data to complement Figure 3 and Figure 4 in the description of the method in Section 2.2.
Figure A1. Station data (vertical) vs. month number (horizontal) for station KCWV (Claxton Evans County Airport, 32.19510° N, 81.86960° W 111 ft in Georgia) from 2008 to 2016 with 2 burn windows per year. Points that do and do not meet the burn-requirements (Table 2) are dark and light gray, respectively, and enumerated in the legends and titles, respectively. The abscissa is t ( mod 12 ) for all the years indicated in the title.
Figure A1. Station data (vertical) vs. month number (horizontal) for station KCWV (Claxton Evans County Airport, 32.19510° N, 81.86960° W 111 ft in Georgia) from 2008 to 2016 with 2 burn windows per year. Points that do and do not meet the burn-requirements (Table 2) are dark and light gray, respectively, and enumerated in the legends and titles, respectively. The abscissa is t ( mod 12 ) for all the years indicated in the title.
Atmosphere 09 00296 g0a1
Figure A2. Station data (vertical) vs. month number (horizontal) for station KLHW (Ft. Stewart, 31.88333° N, 81.56667° W, 46 ft in Georgia) from 2002 to 2016 with 2 burn windows per year. Points that do and do not meet the burn-requirements (Table 2) are dark and light gray, respectively, and enumerated in the legends and titles, respectively. The abscissa is t ( mod 12 ) for all the years indicated in the title.
Figure A2. Station data (vertical) vs. month number (horizontal) for station KLHW (Ft. Stewart, 31.88333° N, 81.56667° W, 46 ft in Georgia) from 2002 to 2016 with 2 burn windows per year. Points that do and do not meet the burn-requirements (Table 2) are dark and light gray, respectively, and enumerated in the legends and titles, respectively. The abscissa is t ( mod 12 ) for all the years indicated in the title.
Atmosphere 09 00296 g0a2
Figure A3. Station data (vertical) vs. month number (horizontal) for station LCSS1 (Savriv, 33.33305° N, 81.591667° W, 275 ft in South Carolina) from 2004 to 2016. Points that do and do not meet the burn-requirements (Table 2) are dark and light gray, respectively, and enumerated in the legends and titles, respectively. The abscissa is t ( mod 12 ) for all the years indicated in the title.
Figure A3. Station data (vertical) vs. month number (horizontal) for station LCSS1 (Savriv, 33.33305° N, 81.591667° W, 275 ft in South Carolina) from 2004 to 2016. Points that do and do not meet the burn-requirements (Table 2) are dark and light gray, respectively, and enumerated in the legends and titles, respectively. The abscissa is t ( mod 12 ) for all the years indicated in the title.
Atmosphere 09 00296 g0a3
Figure A4. Station data (vertical) vs. month number (horizontal) for station QLBA3 (Lindbergh Hill, 36.285556° N, 112.078611° W, 8800 ft in Arizona) for 1999–2016. Points that do and do not meet the burn-requirements (Table 2) are dark and light gray, respectively, and enumerated in the legends and titles, respectively. The abscissa is t ( mod 12 ) for all the years indicated in the title.
Figure A4. Station data (vertical) vs. month number (horizontal) for station QLBA3 (Lindbergh Hill, 36.285556° N, 112.078611° W, 8800 ft in Arizona) for 1999–2016. Points that do and do not meet the burn-requirements (Table 2) are dark and light gray, respectively, and enumerated in the legends and titles, respectively. The abscissa is t ( mod 12 ) for all the years indicated in the title.
Atmosphere 09 00296 g0a4
Figure A5. Station data (vertical) vs. month number (horizontal) for station TT084 (Fishlake D4 Pt #4, 38.960319° N, 111.405983° W, 8523 ft in Utah) for 2012–2016. Points that do and do not meet the burn-requirements (Table 2) are dark and light gray, respectively, and enumerated in the legends and titles, respectively. The abscissa is t ( mod 12 ) for all the years indicated in the title.
Figure A5. Station data (vertical) vs. month number (horizontal) for station TT084 (Fishlake D4 Pt #4, 38.960319° N, 111.405983° W, 8523 ft in Utah) for 2012–2016. Points that do and do not meet the burn-requirements (Table 2) are dark and light gray, respectively, and enumerated in the legends and titles, respectively. The abscissa is t ( mod 12 ) for all the years indicated in the title.
Atmosphere 09 00296 g0a5
Figure A6. As in Figure 4 but for station KCVW.
Figure A6. As in Figure 4 but for station KCVW.
Atmosphere 09 00296 g0a6
Figure A7. As in Figure 4 but for station KLHW.
Figure A7. As in Figure 4 but for station KLHW.
Atmosphere 09 00296 g0a7
Figure A8. As in Figure 4 but for station LCSS1.
Figure A8. As in Figure 4 but for station LCSS1.
Atmosphere 09 00296 g0a8
Figure A9. As in Figure 4 but for station QLBA3.
Figure A9. As in Figure 4 but for station QLBA3.
Atmosphere 09 00296 g0a9
Figure A10. As in Figure 4 but for station TT084.
Figure A10. As in Figure 4 but for station TT084.
Atmosphere 09 00296 g0a10

References

  1. Mandel, J.; Amram, S.; Beezley, J.D.; Kelman, G.; Kochanski, A.K.; Kondratenko, V.Y.; Lynn, B.H.; Regev, B.; Vejmelka, M. Recent advances and applications of WRF-SFIRE. Nat. Hazards Earth Syst. Sci. 2014, 14, 2829–2845. [Google Scholar] [CrossRef]
  2. Larkin, N.K.; O’Neill, S.M.; Solomon, R.; Raffuse, S.; Strand, T.; Sullivan, D.C.; Krull, C.; Rorig, M.; Peterson, J.; Ferguson, S.A. The BlueSky smoke modeling framework. Int. J. Wildland Fire 2009, 18, 906–920. [Google Scholar] [CrossRef]
  3. Seiler, W.; Crutzen, P.J. Estimates of gross and net fluxes of carbon between the biosphere and the atmosphere from biomass burning. Clim. Chang. 1980, 2, 207–247. [Google Scholar] [CrossRef]
  4. Urbanski, S.P.; Hao, W.M.; Nordgren, B. The wildland fire emission inventory: Western United States emission estimates and an evaluation of uncertainty. Atmos. Chem. Phys. 2011, 11, 12973–13000. [Google Scholar] [CrossRef] [Green Version]
  5. Liu, Y.; Qu, J.J.; Wang, W.; Hao, X. Estimates of Wildland Fire Emissions. In Remote Sensing and Modeling Applications to Wildland Fires; Qu, J.J., Sommers, W.T., Yang, R., Riebau, A.R., Eds.; Springer: Berlin/Heidelberg, Germany, 2013; pp. 117–133. [Google Scholar] [CrossRef]
  6. Colarco, P.R.; Schoeberl, M.R.; Doddridge, B.G.; Marufu, L.T.; Torres, O.; Welton, E.J. Transport of smoke from Canadian forest fires to the surface near Washington, D.C.: Injection height, entrainment, and optical properties. J. Geophys. Res. Atmos. 2004, 109, D06203. [Google Scholar] [CrossRef]
  7. Liu, Y.; Achtemeier, G.; Goodrick, S. Sensitivity of air quality simulation to smoke plume rise. J. Appl. Remote Sens. 2008, 2, 021503. [Google Scholar] [CrossRef]
  8. Goodrick, S.L.; Achtemeier, G.L.; Larkin, N.K.; Liu, Y.; Strand, T.M. Modelling smoke transport from wildland fires: A review. Int. J. Wildland Fire 2013, 22, 83–94. [Google Scholar] [CrossRef]
  9. Ottmar, R.; Brown, T.J.; French, N.H.F.; Larkin, N.K. Fire and Smoke Model Evaluation Experiment (FASMEE) Study Plan. 2017. Available online: https://www.fasmee.net/wp-content/uploads/2017/07/FASMEE_StudyPlan_Final_07-11-17.pdf (accessed on 30 June 2018).
  10. Liu, Y.; Kochanski, A.; Baker, K.; Mell, R.; Linn, R.; Paugam, R.; Mandel, J.; Fournier, A.; Jenkins, M.A.; Goodrick, S.; et al. Fire and Smoke Model Evaluation Experiment (FASMEE): Modeling Gaps and Data Needs. In Proceedings of the 2nd International Smoke Symposium, Long Beach, CA, USA, 14–17 November 2016; Available online: https://www.fs.fed.us/rm/pubs_journals/2017/rmrs_2017_liu_y001.pdf (accessed on 30 June 2018).
  11. Mandel, J.; Beezley, J.D.; Kochanski, A.K. Coupled atmosphere-wildland fire modeling with WRF 3.3 and SFIRE 2011. Geosci. Model Dev. 2011, 4, 591–610. [Google Scholar] [CrossRef] [Green Version]
  12. Kochanski, A.K.; Jenkins, M.A.; Yedinak, K.; Mandel, J.; Beezley, J.; Lamb, B. Toward an integrated system for fire, smoke, and air quality simulations. Int. J. Wildland Fire 2016, 25, 534–546. [Google Scholar] [CrossRef]
  13. The “Matash” Fire-Prediction System. 2017. Available online: http://mops.gov.il/English/HomelandSecurityENG/NFServices/Pages/FirePredictionSystem.aspx (accessed on 14 May 2017).
  14. Val Martin, M.; Logan, J.A.; Kahn, R.A.; Leung, F.Y.; Nelson, D.L.; Diner, D.J. Smoke injection heights from fires in North America: Analysis of 5 years of satellite observations. Atmos. Chem. Phys. 2010, 10, 1491–1510. [Google Scholar] [CrossRef]
  15. Mallia, D.V.; Kochanski, A.K.; Urbanski, S.P.; Lin, J.C. Optimizing Smoke and Plume Rise Modeling Approaches at Local Scales. Atmosphere 2018, 9, 166. [Google Scholar] [CrossRef]
  16. Kochanski, A.; Mandel, J.; Fournier, A.; Jenkins, M.A. Modeling Support for FASMEE Experimental Design Using WRF-SFIRE-CHEM; Final Report for Joint Fire Science Program Project JFSP 16-4-05-3, 2017. Available online: https://www.firescience.gov/projects/16-4-05-3/project/16-4-05-3_final_report.pdf (accessed on 15 July 2018).
  17. Kochanski, A.K.; Jenkins, M.A.; Mandel, J.; Beezley, J.D.; Clements, C.B.; Krueger, S. Evaluation of WRF-SFIRE Performance with Field Observations from the FireFlux experiment. Geosci. Model Dev. 2013, 6, 1109–1126. [Google Scholar] [CrossRef]
  18. Kochanski, A.K.; Jenkins, M.A.; Kondratenko, V.; Schranz, S.; Vejmelka, M.; Clements, C.B.; Davis, B. Ignition from fire perimeter and assimilation into a coupled fire-atmosphere model. Proceedings for the 5th International Fire Behavior and Fuels Conference, Portland, OR, USA, 11–15 April 2016; International Association of Wildland Fire: Missoula, MT, USA, 2016; pp. 80–85. Available online: https://www.iawfonline.org/conference-proceedings (accessed on 27 July 2018).
  19. Mandel, J.; Beezley, J.D.; Kochanski, A.K.; Kondratenko, V.Y.; Kim, M. Assimilation of Perimeter Data and Coupling with Fuel Moisture in a Wildland Fire—Atmosphere DDDAS. Procedia Comput. Sci. 2012, 9, 1100–1109. [Google Scholar] [CrossRef]
  20. Mesinger, F.; DiMego, G.; Kalnay, E.; Mitchell, K.; Shafran, P.C.; Ebisuzaki, W.; Jović, D.; Woollen, J.; Rogers, E.; Berbery, E.H.; et al. North American Regional Reanalysis. Bull. Am. Meteorol. Soc. 2006, 87, 343–360. [Google Scholar] [CrossRef] [Green Version]
  21. Vejmelka, M. WRFxpy. 2016. Available online: https://github.com/openwfm/wrfxpy (accessed on 15 June 2018).
  22. Mandel, J.; Beezley, J.D.; Coen, J.L.; Kim, M. Data Assimilation for Wildland Fires: Ensemble Kalman filters in coupled atmosphere-surface models. IEEE Control Syst. Mag. 2009, 29, 47–65. [Google Scholar] [CrossRef]
  23. Mandel, J.; Beezley, J.D.; Kochanski, A.K.; Kondratenko, V.Y.; Zhang, L.; Anderson, E.; Daniels, J., II; Silva, C.T.; Johnson, C.R. A Wildland Fire Modeling and Visualization Environment; Paper 6.4, Ninth Symposium on Fire and Forest Meteorology, Palm Springs, October 2011. Available online: http://ams.confex.com/ams/9FIRE/webprogram/Paper192277.html (accessed on 27 July 2018).
  24. Clark, T.L.; Jenkins, M.A.; Coen, J.; Packham, D. A Coupled Atmospheric-Fire Model: Convective Feedback on Fire Line Dynamics. J. Appl. Meteorol. 1996, 35, 875–901. [Google Scholar] [CrossRef]
  25. Clark, T.L.; Coen, J.; Latham, D. Description of a Coupled Atmosphere-Fire Model. Int. J. Wildland Fire 2004, 13, 49–64. [Google Scholar] [CrossRef]
  26. Coen, J.L.; Cameron, M.; Michalakes, J.; Patton, E.G.; Riggan, P.J.; Yedinak, K. WRF-Fire: Coupled Weather-Wildland Fire Modeling with the Weather Research and Forecasting Model. J. Appl. Meteorol. Climatol. 2013, 52, 16–38. [Google Scholar] [CrossRef]
  27. Jiménez, P.A.; Muñoz-Esparza, D.; Kosović, B. A High Resolution Coupled Fire-Atmosphere Forecasting System to Minimize the Impacts of Wildland Fires: Applications to the Chimney Tops II Wildland Event. Atmosphere 2018, 9, 197. [Google Scholar] [CrossRef]
  28. Muñoz-Esparza, D.; Kosović, B.; Jiménez, P.A.; Coen, J.L. An Accurate Fire-Spread Algorithm in the Weather Research and Forecasting Model Using the Level-Set Method. J. Adv. Model. Earth Syst. 2018, 10, 908–926. [Google Scholar] [CrossRef] [Green Version]
  29. McKay, M.D.; Beckman, R.J.; Conover, W.J. A comparison of three methods for selecting values of input variables in the analysis of output from a computer code. Technometrics 1979, 21, 239–245. [Google Scholar] [CrossRef]
  30. McKay, M.D. Evaluating Prediction Uncertainty; Technical Report NUREG/CR-6311 LA-12915-MS; Los Alamos National Laboratory: Los Alamos, NM, USA, 1995; Available online: http://www.iaea.org/inis/collection/NCLCollectionStore/_Public/26/051/26051087.pdf (accessed on 16 October 2016).
  31. Saltelli, A.; Tarantola, S.; Campolongo, F.; Ratto, M. Sensitivity Analysis in Practice; John Wiley & Sons, Ltd.: Chichester, UK, 2004; p. xii+219. [Google Scholar]
  32. Sobol’, I.M. Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates. Math. Comput. Simul. 2001, 55, 271–280. [Google Scholar] [CrossRef]
  33. Chugunov, N.; Ramakrishnan, T.S.; Moldoveanu, N.; Fournier, A.; Osypov, K.S.; Altundas, Y.B. Targeted Survey Design Under Uncertainty. US Patent Application No. 0278110 A1, 18 September 2014. [Google Scholar]
  34. Rothermel, R.C. A Mathematical Model for Predicting Fire Spread in Wildland Fires; USDA Forest Service Research Paper INT-115, 1972. Available online: https://www.fs.fed.us/rm/pubs_int/int_rp115.pdf (accessed on 15 March 2018).
  35. Clements, C.B.; Davis, B.; Seto, D.; Contezac, J.; Kochanski, A.; Filippi, J.B.; Lareau, N.; Barboni, T.; Butler, B.; Krueger, S.; et al. Overview of the 2013 FireFlux-II Grass Fire Field Experiment. In Advances in Forest Fire Research; Viegas, D.X., Ed.; Coimbra University Press: Coimbra, Portugal, 2014; pp. 392–400. [Google Scholar] [CrossRef]
Figure 1. Time series of Mahalanobis (standardized) deviation norm δ t for each station (6 rows). Points for states that do and do not meet burn requirements are dark and light gray, respectively. For each station, the most typical day in the burn window has the smallest δ t and is listed in the title and indicated with a ∇ marker. The abscissa is t ( mod 12 ) for all the years indicated in the plot titles.
Figure 1. Time series of Mahalanobis (standardized) deviation norm δ t for each station (6 rows). Points for states that do and do not meet burn requirements are dark and light gray, respectively. For each station, the most typical day in the burn window has the smallest δ t and is listed in the title and indicated with a ∇ marker. The abscissa is t ( mod 12 ) for all the years indicated in the plot titles.
Atmosphere 09 00296 g001
Figure 2. WRF domain setups for burn simulation: (a) Fishlake; (b) North Kaibab; (c) Fort Stewart, with indication of nearby meteorological stations used to define typical days.
Figure 2. WRF domain setups for burn simulation: (a) Fishlake; (b) North Kaibab; (c) Fort Stewart, with indication of nearby meteorological stations used to define typical days.
Atmosphere 09 00296 g002
Figure 3. Station data (vertical) vs. month number (horizontal) for station FSHU1 (Fish Lake ranger station near Koosharem 10ENE, 38.55167° N, 111.72278° W, 8880 ft in Utah) in years 2001–2016. Points that do and do not meet the burn-requirements (Table 2) are dark and light gray, respectively, and enumerated in the legends and titles, respectively. The abscissa is t ( mod 12 ) for all the years indicated in the title.
Figure 3. Station data (vertical) vs. month number (horizontal) for station FSHU1 (Fish Lake ranger station near Koosharem 10ENE, 38.55167° N, 111.72278° W, 8880 ft in Utah) in years 2001–2016. Points that do and do not meet the burn-requirements (Table 2) are dark and light gray, respectively, and enumerated in the legends and titles, respectively. The abscissa is t ( mod 12 ) for all the years indicated in the title.
Atmosphere 09 00296 g003
Figure 4. Scatter plots of FSHU1 pairs of state variables. Black and gray indicate values that are and are not constrained by the burn requirements. The marginal covariance of a pair of variables is depicted by the green ellipse centered at the mean and with axes equal to the deviations in the principal directions. The ellipse projections on the coordinate axes have lengths equaling twice the corresponding marginal standard deviations. Ellipse area decreases with increasing correlation (if tilting up) or decreasing anti-correlation (if tilting down); a circle would indicate zero correlation. Only the values that satisfy the minimum and maximum values for burn requirements in Table 2 are used to compute the covariances, in the considered years regardless of burn season.
Figure 4. Scatter plots of FSHU1 pairs of state variables. Black and gray indicate values that are and are not constrained by the burn requirements. The marginal covariance of a pair of variables is depicted by the green ellipse centered at the mean and with axes equal to the deviations in the principal directions. The ellipse projections on the coordinate axes have lengths equaling twice the corresponding marginal standard deviations. Ellipse area decreases with increasing correlation (if tilting up) or decreasing anti-correlation (if tilting down); a circle would indicate zero correlation. Only the values that satisfy the minimum and maximum values for burn requirements in Table 2 are used to compute the covariances, in the considered years regardless of burn season.
Atmosphere 09 00296 g004
Figure 5. Latin Hypercube Sampling (LHS). Each of L = 7 parameters takes the values from its N = 5 sampling points. The rows are random permutations of the numbers 1 to 5, while the columns are the parameter vectors on which the simulations will be run.
Figure 5. Latin Hypercube Sampling (LHS). Each of L = 7 parameters takes the values from its N = 5 sampling points. The rows are random permutations of the numbers 1 to 5, while the columns are the parameter vectors on which the simulations will be run.
Atmosphere 09 00296 g005
Figure 6. Time series of the maximum updraft strength and location from the Fishlake simulation for 09/03/2014. (a) Standard simulation; (b) Simulation with doubled fire heat flux released to the atmosphere.
Figure 6. Time series of the maximum updraft strength and location from the Fishlake simulation for 09/03/2014. (a) Standard simulation; (b) Simulation with doubled fire heat flux released to the atmosphere.
Atmosphere 09 00296 g006
Figure 7. Smoke analysis for the Fishlake burn executed for 09/03/2014. (a) Vertical profile of normalized smoke concentration; (b) Time series of the plume top height.
Figure 7. Smoke analysis for the Fishlake burn executed for 09/03/2014. (a) Vertical profile of normalized smoke concentration; (b) Time series of the plume top height.
Atmosphere 09 00296 g007
Figure 8. Analysis of the fire-induced winds for the Fishlake simulation executed for 09/03/2014. (a) Time series of the maximum fire-induced horizontal winds and their heights in the baseline case; (b) Same in the double heat case.
Figure 8. Analysis of the fire-induced winds for the Fishlake simulation executed for 09/03/2014. (a) Time series of the maximum fire-induced horizontal winds and their heights in the baseline case; (b) Same in the double heat case.
Atmosphere 09 00296 g008
Figure 9. Smoke analysis for the North Kaibab burn simulation executed for 09/05/2008. (a) Vertical profile of normalized smoke concentrations; (b) Time series of the simulated plume top height.
Figure 9. Smoke analysis for the North Kaibab burn simulation executed for 09/05/2008. (a) Vertical profile of normalized smoke concentrations; (b) Time series of the simulated plume top height.
Atmosphere 09 00296 g009
Figure 10. Analysis of the fire-induced wind and updraft for the North Kaibab simulation executed for 09/05/2008. (a) Time series of the maximum fire-induced horizontal winds; (b) TTime series of the maximum fire-induced updraft.
Figure 10. Analysis of the fire-induced wind and updraft for the North Kaibab simulation executed for 09/05/2008. (a) Time series of the maximum fire-induced horizontal winds; (b) TTime series of the maximum fire-induced updraft.
Atmosphere 09 00296 g010
Figure 11. Results from numerical simulations for Fort Stewart (02/14/2013) with different ignition procedures. (a) Time series of the simulated maximum fire-induced updraft; (b) Same for plume top height.
Figure 11. Results from numerical simulations for Fort Stewart (02/14/2013) with different ignition procedures. (a) Time series of the simulated maximum fire-induced updraft; (b) Same for plume top height.
Atmosphere 09 00296 g011
Figure 12. Analysis of the fire-induced wind and updraft for the Fort Stewart simulation executed for 04/22/2014. (a) Time series of the fire-induced maximum horizontal wind speed; (b) Same for the updraft.
Figure 12. Analysis of the fire-induced wind and updraft for the Fort Stewart simulation executed for 04/22/2014. (a) Time series of the fire-induced maximum horizontal wind speed; (b) Same for the updraft.
Atmosphere 09 00296 g012
Figure 13. (a) The time series of the variance of vertical velocity at different heights above the ground from simulations executed with parameters from the rLHS sampling; (b) Map of the variance of the vertical velocity at 1200 m with the indication of the burn plot (green contour) and the ignition (red line).
Figure 13. (a) The time series of the variance of vertical velocity at different heights above the ground from simulations executed with parameters from the rLHS sampling; (b) Map of the variance of the vertical velocity at 1200 m with the indication of the burn plot (green contour) and the ignition (red line).
Atmosphere 09 00296 g013
Figure 14. (a) The mean plume top height; (b) The normalized variance in the plume height. The green contour represents the burning plot boundary and the red line shows the ignition.
Figure 14. (a) The mean plume top height; (b) The normalized variance in the plume height. The green contour represents the burning plot boundary and the red line shows the ignition.
Atmosphere 09 00296 g014
Figure 15. (a) The mean normalized smoke concentration at 1400 m above the terrain; (b) The normalized smoke concentration average variance at 1400 m above the terrain. The green contour represents the burning plot boundary, the red line shows the ignition.
Figure 15. (a) The mean normalized smoke concentration at 1400 m above the terrain; (b) The normalized smoke concentration average variance at 1400 m above the terrain. The green contour represents the burning plot boundary, the red line shows the ignition.
Atmosphere 09 00296 g015
Figure 16. Sensitivity of the vertical velocity at 1200 m (top), smoke concentration at 1400 m (middle), and plume top height (bottom) to the fire heat flux multiplier (left column) and the heat extinction depth (right column).
Figure 16. Sensitivity of the vertical velocity at 1200 m (top), smoke concentration at 1400 m (middle), and plume top height (bottom) to the fire heat flux multiplier (left column) and the heat extinction depth (right column).
Atmosphere 09 00296 g016
Figure 17. (a) The mean normalized smoke concentration at 1400 m above the ground computed from 5 most typical days; (b) The normalized average variance of the smoke concentration at 1400 m computed from 5 most typical days. Green contour represents burn unit boundaries, red line shows the ignition.
Figure 17. (a) The mean normalized smoke concentration at 1400 m above the ground computed from 5 most typical days; (b) The normalized average variance of the smoke concentration at 1400 m computed from 5 most typical days. Green contour represents burn unit boundaries, red line shows the ignition.
Atmosphere 09 00296 g017
Figure 18. (a) The vertical velocity average variance at 1000 m computed from 5 most typical days; (b) The mean plume top height computed from 5 most typical days. Green contour represents burn unit boundaries, white line shows the ignition.
Figure 18. (a) The vertical velocity average variance at 1000 m computed from 5 most typical days; (b) The mean plume top height computed from 5 most typical days. Green contour represents burn unit boundaries, white line shows the ignition.
Atmosphere 09 00296 g018
Table 1. Parameters of the experimental burn simulations.
Table 1. Parameters of the experimental burn simulations.
Burn SiteFishlakeNorth KaibabFort Stewart
Meteo forcingNARRNARRNARR
Number of domains5 (realistic)5 (realistic)5 (realistic)
Domain sizes(XYZ) 97 × 97 × 41 97 × 97 × 41 97 × 97 × 41
Model top13.2 km13.2 km13.2 km
Horizontal resolution12 km/4 km/1.33 km/444 m/148 m12 km/4 km/1.33 km/444 m/148 m12 km/4 km/1.33 km/444 m/148 m
Vertical resolution5.3 m–2233 m5.3 m–2233 m5.3 m–2233 m
Fire mesh resolution29.6 m29.6 m29.6 m
Total number of runs7410
IgnitionHelicopterPointStraight Line/Point
Simulation start09/03/2014 00:00 UTC09/05/2008 00:00 UTC04/22/2014 00:00 UTC
09/11/2016 00:00 UTC09/19/2015 00:00 UTC04/27/2009 00:00 UTC
09/22/2012 00:00 UTC09/01/2001 00:00 UTC02/27/2013 00:00 UTC
09/26/2015 00:00 UTC09/02/2011 00:00 UTC02/26/2008 00:00 UTC
Ignition Time15:00 UTC (9:00 local)15:00 UTC (9:00 local)15:00 UTC (11:00 local)
Simulation Length48 h48 h48 h
Fire Output Interval5 min5 min5 min
Time step (d05)0.5 s0.5 s0.74 s
1 h fuel moisture6.0%6.0%15.0%
10 h fuel moisture8.0%8.0%13.0%
100 h fuel moisture9.0%9.0%13.0%
1000 h fuel moisture12.0%12.0%18.0%
Table 2. Stations, time windows (mm/dd hh:mm UTC) and value limits for potential burns (burn windows). T—temperature, ϕ —relative humidity, s—wind speed.
Table 2. Stations, time windows (mm/dd hh:mm UTC) and value limits for potential burns (burn windows). T—temperature, ϕ —relative humidity, s—wind speed.
StationBurn SeasonMin TMax TMin ϕ Max ϕ Min sMax s
FSHU109/01 16:00 UTC–10/31 18:00 UTC61 F85 F16%22%015 mph
KCWV02/20 16:00 UTC–02/28 18:00 UTC60903055620
04/15 13:00UTC–04/30 14:00 UTC60903055620
KLHW02/20 16:00 UTC–02/28 18:00 UTC60903055620
04/15 13:00 UTC–04/30 14:00 UTC60903055620
LCSS101/01 15:00 UTC–02/01 19:00 UTC60903055620
12/01 15:00 UTC–12/31 19:00 UTC60903055620
QLBA309/01 17:00 UTC–10/31 19:00 UTC61851622015
TT08409/01 16:00 UTC–10/31 18:00 UTC61851622015
Table 3. Summary of most typical days identified for the burn sites.
Table 3. Summary of most typical days identified for the burn sites.
Burn SiteFishlakeNorth KaibabFort Stewart
Weather StationFSHU1TT084QLBA3KCWVKLHW
Typical days09/03/201409/05/201209/05/200804/22/201404/27/2009
09/26/201509/04/201209/19/200502/24/201302/26/2013
Table 4. Sampling points for Latin Hypercube Sampling.
Table 4. Sampling points for Latin Hypercube Sampling.
Sampling Point12345
10-fuel moisture (kg/kg)0.04590.06130.07480.09140.1219
Heat extinction depth (m)7.583012.353117.320524.285439.5620
Heat flux multiplier (1)0.58270.80171.00001.24731.7161
Multiplier for R (1)0.58270.80171.00001.24731.7161
Multiplier for R W (1)0.58270.80171.00001.24731.7161
Multiplier for R S (1)0.58270.80171.00001.24731.7161
Simulation day (date)09/03/201409/11/201609/22/201209/26/201509/27/2015

Share and Cite

MDPI and ACS Style

Kochanski, A.K.; Fournier, A.; Mandel, J. Experimental Design of a Prescribed Burn Instrumentation. Atmosphere 2018, 9, 296. https://doi.org/10.3390/atmos9080296

AMA Style

Kochanski AK, Fournier A, Mandel J. Experimental Design of a Prescribed Burn Instrumentation. Atmosphere. 2018; 9(8):296. https://doi.org/10.3390/atmos9080296

Chicago/Turabian Style

Kochanski, Adam K., Aimé Fournier, and Jan Mandel. 2018. "Experimental Design of a Prescribed Burn Instrumentation" Atmosphere 9, no. 8: 296. https://doi.org/10.3390/atmos9080296

APA Style

Kochanski, A. K., Fournier, A., & Mandel, J. (2018). Experimental Design of a Prescribed Burn Instrumentation. Atmosphere, 9(8), 296. https://doi.org/10.3390/atmos9080296

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop