APPLICATION OF AN IMPUTATION METHOD FOR GEOSPATIAL INVENTORY OF FOREST STRUCTURAL ATTRIBUTES ACROSS MULTIPLE SPATIAL SCALES IN THE LAKE STATES, U.S.A.

Credible spatial information characterizing the structure and site quality of forests is critical to sustainable forest management and planning, especially given the increasing demands and threats to forest products and services. Forest managers and planners are required to evaluate forest conditions over a broad range of scales, contingent on operational or reporting requirements. Traditionally, forest inventory estimates are generated via a design-based approach that involves generalizing sample plot measurements to characterize an unknown population across a larger area of interest. However, field plot measurements are costly and as a consequence spatial coverage is limited. Remote sensing technologies have shown remarkable success in augmenting limited sample plot data to generate stand- and landscape-level spatial predictions of forest inventory attributes. Further enhancement of forest inventory approaches that couple field measurements with cutting edge remotely sensed and geospatial datasets are essential to sustainable forest management. We evaluated a novel Random Forest based k Nearest Neighbors (RF-kNN) imputation approach to couple remote sensing and geospatial data with field inventory collected by different sampling methods to generate forest inventory information across large spatial extents. The forest inventory data collected by the FIA program of US Forest Service was integrated with optical remote sensing and other geospatial datasets to produce biomass distribution maps for a part of the Lake States and species-specific site index maps for the entire Lake State. Targeting small-area application of the state-of-art remote sensing, LiDAR (light detection and ranging) data was integrated with the field data collected by an inexpensive method, called variable plot sampling, in the Ford Forest of Michigan Tech to derive standing volume map in a cost-effective way. The outputs of the RF-kNN imputation were compared with independent validation datasets and extant map products based on different sampling and modeling strategies. The RF-kNN modeling approach was found to be very effective, especially for large-area estimation, and produced results statistically equivalent to the field observations or the estimates derived from secondary data sources. The models are useful to resource managers for operational and strategic purposes.

Tropospheric ozone and CO over North Atlantic for the past decade

We investigate how declines in US emissions of CO and O3 precursors have impacted the lower free troposphere over the North Atlantic. We use seasonal observations for O3 and CO from the PICO-NARE project for the period covering 2001 to 2010. Observations are used to verify model output generated by the GEOS-Chem 3-D global chemical transport model. Additional satellite data for CO from AIRS/Aqua and for O3 from TES/Aura were also used to provide additional comparisons; particularly for fall, winter, and spring when PICO-NARE coverage is sparse. We find GEOS-Chem captures the seasonal cycle for CO and O3 well compared to PICO-NARE data. For CO, GEOS-Chem is biased low, particularly in spring which is in agreement with findings from previous studies. GEOS-Chem is 24.7 +/- 5.2 ppbv (1-σ) low compared to PICO-NARE summer CO data while AIRS is 14.2 +/- 6.6 ppbv high. AIRS does not show nearly as much variation as seen with GEOS-Chem or the Pico data, and goes from being lower than PICO-NARE data in winter and spring, to higher in summer and fall. Both TES and GEOS-Chem match the seasonal ozone cycle well for all seasons when compared with observations. Model results for O3 show GEOS-Chem is 6.67 +/- 2.63 ppbv high compared to PICO-NARE summer measurements and TES was 3.91 +/- 4.2 ppbv higher. Pico data, model results, and AIRS all show declines in CO and O3 for the summer period from 2001 to 2010. Limited availability of TES data prevents us from using it in trend analysis. For summer CO Pico, GEOS-Chem, and AIRS results show declines of 1.32, 0.368, and 0.548 ppbv/year respectively. For summer O3, Pico and GEOS-Chem show declines of -0.726 and -0.583 ppbv/year respectively. In other seasons, both model and AIRS show declining CO, particularly in the fall. GEOS-Chem results show a fall decline of 0.798 ppbv/year and AIRS shows a decline of 0.8372 ppbv/year. Winter and spring CO declines are 0.393 and 0.307 for GEOS-Chem, and 0.455 and 0.566 for AIRS. GEOS-Chem shows declining O3 in other seasons as well; with fall being the season of greatest decrease and winter being the least. Model results for fall, winter, and spring are 0.856, 0.117, and 0.570 ppbv/year respectively. Given the availability of data we are most confident in summer results and thus find that summer CO and O3 have declined in lower free troposphere of the North Atlantic region of the Azores. Sensitivity studies for CO and O3 at Pico were conducted by turning off North American fossil fuel emissions in GEOS-Chem. Model results show that North America fossil fuel emissions contribute 8.57 ppbv CO and 4.03 ppbv O3 to Pico. The magnitude of modeled trends declines in all seasons without North American fossil fuel emissions except for summer CO. The increase in summer CO declines may be due to a decline of 5.24 ppbv/year trend in biomass burning emissions over the study period; this is higher than the 2.33 ppbv/year North American anthropogenic CO model decline. Winter O3 is the only season which goes from showing a negative trend to a positive trend.

ASSESSMENT OF ASPEN AND NORTHERN HARDWOODS EXTENT IN THIRTY THREE COUNTIES OF UPPER/LOWER MICHIGAN

The importance of the United States' wood and wood byproducts as biomass feedstocks is increasing as the concern about security and sustainability of global energy production continues to rise. Thus, second generation woody feedstock sources in Michigan, e.g., hybrid poplar and hybrid willow (Populus spp.), are viewed as a potential source of biomass for the proposed biofuel ethanol production plant in Kinross, MI. It is important to gain an understanding of the spatial distribution of current feedstock sources, harvesting accessibility via the transportation infrastructure and land ownerships in order to ensure long-term feedstock extent. This research provides insights into the current extent of aspen and northern hardwoods, and an assessment of potential for expanding the area of these feedstock sources based on pre-European settlement conditions. A geographic information system (GIS) was developed to compile available geospatial data for 33 counties located within 150 miles of the Kinross facility. These include present day and pre-European settlement land use/cover, soils, road infrastructure, and land ownerships. The results suggest that a significant amount of northern hardwoods has been converted to other land use/cover types since European settlement, and the "scattering" of aspen stands has increased. Furthermore, a significant amount of woody biomass is available in close proximity to the existing road network, which can be effectively utilized as feedstock. Potential aspen and northern hardwoods restoration areas are identified in the vicinity of road networks which can be used for future woody feedstock production.

Modeling and classifying variable width riparian zones utilizing digital elevation models, flood height data, digital soil data and national wetlands inventory : a new approach for riparian zone delineation

Riparian zones are dynamic, transitional ecosystems between aquatic and terrestrial ecosystems with well defined vegetation and soil characteristics. Development of an all-encompassing definition for riparian ecotones, because of their high variability, is challenging. However, there are two primary factors that all riparian ecotones are dependent on: the watercourse and its associated floodplain. Previous approaches to riparian boundary delineation have utilized fixed width buffers, but this methodology has proven to be inadequate as it only takes the watercourse into consideration and ignores critical geomorphology, associated vegetation and soil characteristics. Our approach offers advantages over other previously used methods by utilizing: the geospatial modeling capabilities of ArcMap GIS; a better sampling technique along the water course that can distinguish the 50-year flood plain, which is the optimal hydrologic descriptor of riparian ecotones; the Soil Survey Database (SSURGO) and National Wetland Inventory (NWI) databases to distinguish contiguous areas beyond the 50-year plain; and land use/cover characteristics associated with the delineated riparian zones. The model utilizes spatial data readily available from Federal and State agencies and geospatial clearinghouses. An accuracy assessment was performed to assess the impact of varying the 50-year flood height, changing the DEM spatial resolution (1, 3, 5 and 10m), and positional inaccuracies with the National Hydrography Dataset (NHD) streams layer on the boundary placement of the delineated variable width riparian ecotones area. The result of this study is a robust and automated GIS based model attached to ESRI ArcMap software to delineate and classify variable-width riparian ecotones.

Parameter identification of a copper-zeolite SCR catalyst model using reactor data

The selective catalytic reduction system is a well established technology for NOx emissions control in diesel engines. A one dimensional, single channel selective catalytic reduction (SCR) model was previously developed using Oak Ridge National Laboratory (ORNL) generated reactor data for an iron-zeolite catalyst system. Calibration of this model to fit the experimental reactor data collected at ORNL for a copper-zeolite SCR catalyst is presented. Initially a test protocol was developed in order to investigate the different phenomena responsible for the SCR system response. A SCR model with two distinct types of storage sites was used. The calibration process was started with storage capacity calculations for the catalyst sample. Then the chemical kinetics occurring at each segment of the protocol was investigated. The reactions included in this model were adsorption, desorption, standard SCR, fast SCR, slow SCR, NH3 Oxidation, NO oxidation and N2O formation. The reaction rates were identified for each temperature using a time domain optimization approach. Assuming an Arrhenius form of the reaction rates, activation energies and pre-exponential parameters were fit to the reaction rates. The results indicate that the Arrhenius form is appropriate and the reaction scheme used allows the model to fit to the experimental data and also for use in real world engine studies.

Geospatial analysis of western juniper (Juniperus occidentalis) encroachment utilizing remotely sensed data

Utilizing remote sensing methods to assess landscape-scale ecological change are rapidly becoming a dominant force in the natural sciences. Powerful and robust non-parametric statistical methods are also actively being developed to compliment the unique characteristics of remotely sensed data. The focus of this research is to utilize these powerful, robust remote sensing and statistical approaches to shed light on woody plant encroachment into native grasslands--a troubling ecological phenomenon occurring throughout the world. Specifically, this research investigates western juniper encroachment within the sage-steppe ecosystem of the western USA. Western juniper trees are native to the intermountain west and are ecologically important by means of providing structural diversity and habitat for many species. However, after nearly 150 years of post-European settlement changes to this threatened ecosystem, natural ecological processes such as fire regimes no longer limit the range of western juniper to rocky refugia and other areas protected from short fire return intervals that are historically common to the region. Consequently, sage-steppe communities with high juniper densities exhibit negative impacts, such as reduced structural diversity, degraded wildlife habitat and ultimately the loss of biodiversity. Much of today's sage-steppe ecosystem is transitioning to juniper woodlands. Additionally, the majority of western juniper woodlands have not reached their full potential in both range and density. The first section of this research investigates the biophysical drivers responsible for juniper expansion patterns observed in the sage-steppe ecosystem. The second section is a comprehensive accuracy assessment of classification methods used to identify juniper tree cover from multispectral 1 m spatial resolution aerial imagery.

Development of a 1-D Catalyzed Diesel Particulate Filter Model for Simulation of the Performance and the Oxidation of Particulate Matter and Nitrogen Oxides Using Passive Oxidation and Active Regeneration Engine Experimental Data

Back-pressure on a diesel engine equipped with an aftertreatment system is a function of the pressure drop across the individual components of the aftertreatment system, typically, a diesel oxidation catalyst (DOC), catalyzed particulate filter (CPF) and selective catalytic reduction (SCR) catalyst. Pressure drop across the CPF is a function of the mass flow rate and the temperature of the exhaust flowing through it as well as the mass of particulate matter (PM) retained in the substrate wall and the cake layer that forms on the substrate wall. Therefore, in order to control the back-pressure on the engine at low levels and to minimize the fuel consumption, it is important to control the PM mass retained in the CPF. Chemical reactions involving the oxidation of PM under passive oxidation and active regeneration conditions can be utilized with computer numerical models in the engine control unit (ECU) to control the pressure drop across the CPF. Hence, understanding and predicting the filtration and oxidation of PM in the CPF and the effect of these processes on the pressure drop across the CPF are necessary for developing control strategies for the aftertreatment system to reduce back-pressure on the engine and in turn fuel consumption particularly from active regeneration. Numerical modeling of CPF's has been proven to reduce development time and the cost of aftertreatment systems used in production as well as to facilitate understanding of the internal processes occurring during different operating conditions that the particulate filter is subjected to. A numerical model of the CPF was developed in this research work which was calibrated to data from passive oxidation and active regeneration experiments in order to determine the kinetic parameters for oxidation of PM and nitrogen oxides along with the model filtration parameters. The research results include the comparison between the model and the experimental data for pressure drop, PM mass retained, filtration efficiencies, CPF outlet gas temperatures and species (NO2) concentrations out of the CPF. Comparisons of PM oxidation reaction rates obtained from the model calibration to the data from the experiments for ULSD, 10 and 20% biodiesel-blended fuels are presented.