Synthesis of Carbon Nanomaterials through Up-Cycling Agricultural and Municipal Solid Wastes

This work addresses the synthesis of carbon nanomaterials (CNMs) by up-cycling common solid wastes. These feedstocks could supersede the use of costly and often toxic or highly flammable chemicals, such as hydrocarbon gases, carbon monoxide, and hydrogen, which are commonly used as feedstocks in current nanomanufacturing processes for CNMs. Agricultural sugar cane bagasse and corn residues, scrap tire chips, and postconsumer polyethylene (PE) and polyethylene terephthalate (PET) bottle shreddings were either thermally treated by sole pyrolysis or by sequential pyrolysis and partial oxidation. The resulting gaseous carbon-bearing effluents were then channeled into a heated reactor. CNMs, including carbon nanotubes, were catalytically synthesized therein on stainless steel meshes. This work revealed that the structure of the resulting CNMs is determined by the feedstock type, through the disparate mixtures of carbon-bearing gases generated when different feedstocks are pyrolyzed. CNM characterization was conducted by scanning and transmission electron microscopy as well as by Raman spectroscopy and by thermogravimetric analysis. Gas chromatography was used to characterize the gases in the synthesis chamber. This work demonstrated an alternative method for efficient manufacturing of CNMs using both biodegradable and nonbiodegradable agricultural and municipal carbonaceous wastes.

RegCM3 nested in HadAM3 scenarios A2 and B2: projected changes in extratropical cyclogenesis, temperature and precipitation over the South Atlantic Ocean

The RegCM3 (Regional Climate Model-version 3) was nested in HadAM3 model to simulate present (1975-1989, referred hereafter as RegHad) and two future climate scenarios (A2 and B2 from 2071 to 2085, referred as RegA2 and RegB2) over the South America (SA) and South Atlantic Ocean (SAO). Projected changes in the air temperature, precipitation, low level circulation and cyclogenesis climatology were investigated. The cyclogenesis were identified using an automatic scheme for tracking based on the minimum of relative vorticity (zeta) from 10-m height wind. During summer, a general decrease (increase) in the precipitation is projected by RegA2 and RegB2 over the northeastern SA (center-west and south Brazil, north Argentina and Uruguay). For winter, an anomalous low level anticyclonic circulation is associated with the reduction in the rainfall over the central part of southern Brazil in RegA2 and RegB2 scenarios. Similar to HadAM3, RegCM3 projects larger warming in A2 scenario. For the present climate, when compared to HadAM3, RegHad defines better both the location of the main cyclogenetic areas and its annual cycle near southwestern SAO. The projections indicate a reduction in the total number of cyclones of -7.2% and -4.7% for RegA2 and RegB2, respectively, while HadAM3 reduction is -4.5% for both scenarios. The decrease is larger for initially intense cyclones (zeta <=-<= 2.5 x 10(-5) s-(1)): -20.9% (RegA2) and -11.3% (RegB2). For the lifetime, distance traveled and mean velocity of the cyclones, the A2 and B2 scenarios present mean values close to the present climate ( 3 days, 1900 km, and 9 m s(-1), respectively). Regarding the initial mean vorticity of the systems, RegB2 simulates values similar to the present climate, but they are initially weaker in RegA2. In general, RegA2 and RegB2 show a large decrease in the number of cyclones over the southern SAO due to an anticyclonic anomaly covering SAO between 30-55A degrees S. The reduction is larger in the scenario with higher concentrations of greenhouse gases (RegA2).

Aerosols in Amazonia: Natural Biogenic Particles and Large Scale Biomass Burning Impacts.

The Large Scale Biosphere Atmosphere Experiment in Amazonia (LBA) is a long term (20 years) research effort aimed at the understanding of the functioning of the Amazonian ecosystem. In particular, the strong biosphere-atmosphere interaction is a key component looking at the exchange processes between vegetation and the atmosphere, focusing on aerosol particles. Two aerosol components are the most visible: The natural biogenic emissions of aerosols and VOCs, and the biomass burning emissions. A large effort was done to characterize natural biogenic aerosols that showed detailed organic characterization and optical properties. The biomass burning component in Amazonia is important in term of aerosol and trace gases emissions, with deforestation rates decreasing, from 27,000 Km2 in 2004 to about 5,000 Km2 in 2011. Biomass burning emissions in Amazonia increases concentrations of aerosol particles, CO, ozone and other species, and also change the surface radiation balance in a significant way. Long term monitoring of aerosols and trace gases were performed in two sites: a background site in Central Amazonia, 55 Km North of Manaus (called ZF2 ecological reservation) and a monitoring station in Porto Velho, Rondonia state, a site heavily impacted by biomass burning smoke. Several instruments were operated to measured aerosol size distribution, optical properties (absorption and scattering at several wavelengths), composition of organic (OC/EC) and inorganic components among other measurements. AERONET and MODIS measurements from 5 long term sites show a large year-to year variability due to climatic and socio-economic issues. Aerosol optical depths of more than 4 at 550nm was observed frequently over biomass burning areas. In the pristine Amazonian atmosphere, aerosol scattering coefficients ranged between 1 and 200 Mm-1 at 450 nm, while absorption ranged between 1 and 20 Mm-1 at 637 nm. A strong seasonal behavior was observed, with greater aerosol loadings during the dry season (Jul-Nov) as compared to the wet season (Dec-Jun). During the wet season in Manaus, aerosol scattering (450 nm) and absorption (637 nm) coefficients averaged, respectively, 14 and 0.9 Mm-1. Angstrom exponents for scattering were lower during the wet season (1.6) in comparison to the dry season (1.9), which is consistent with the shift from biomass burning aerosols, predominant in the fine mode, to biogenic aerosols, predominant in the coarse mode. Single scattering albedo, calculated at 637 nm, did not show a significant seasonal variation, averaging 0.86. In Porto Velho, even in the wet season it was possible to observe an impact from anthropogenic aerosol. Black Carbon was measured at a high 20 ug/m³ in the dry season, showing strong aerosol absorption. This work presents a general description of the aerosol optical properties in Amazonia, both during the Amazonian wet and dry seasons.

Modelos de circulación global y procesos de mezcla en el océano

La difusividad diapicna en el océano es uno de los parámetros más desconocidos en los modelos climáticos actuales. Su importancia radica en que es uno de los principales factores de transporte de calor hacia capas más profundas del océano. Las medidas de esta difusividad son variables e insuficientes para confeccionar un mapa global con estos valores. A través de una amplia revisión bibliográfica hasta el año 2009 del tema se encontró que el sistema climático es extremadamente sensible a la difusividad diapicna, donde el escalado del Océano Pacífico Sur, con una potencia de su coeficiente de difusividad o kv de 0.63, resultó ser más sensible a los cambios en el coeficiente de difusividad diapicna que el Océano Atlántico con una potencia de kv de 0.44 , se pone de manifiesto así la necesidad de esclarecer los esquemas de mezcla, esquemas de clausura y sus parametrizaciones a través de Modelos de Circulación Global (GCMs) y Modelos de Complejidad Intermedia del Sistema Terrestre (EMICs), dentro del marco de un posible cambio climático y un calentamiento global debido al aumento de las emisiones de gases de efecto invernadero. Así, el objetivo principal de este trabajo es comprender la sensibilidad del sistema climático a la difusividad diapicna en el océano a través de los GCMs y los EMICs. Para esto es necesario el análisis de los posibles esquemas de mezcla diapicna con el objetivo final de encontrar el modelo óptimo que permita predecir la evolución del sistema climático, el estudio de todas las variables que influyen en el mismo, y la correcta simulación en largos periodos de tiempo. The diapycnal diffusivity in the ocean is one of the least known parameters in current climate models. Measurements of this diffusivity are sparse and insufficient for compiling a global map. Through a lengthy review of the literature through 2009 found that the climate system is extremely sensitive to the diapycnal diffusivity, where in the South Pacific scales with the 0.63 power of the diapycnal diffusion, in contrasts to the scales with the 0.44 power of the diapycnal diffusion of North Atlantic. Therefore, the South Pacific is more sensitive than the North Atlantic. All this evidenced the need to clarify the schemes of mixing and its parameterisations through Global Circulation Models (GCMs) and Earth Models of Intermediate Complexity (EMICs) within a context of possible climate change and global warming due to increased of emissions of greenhouse gases. Thus, the main objective of this work understands the sensitivity of the climate system to diapycnal diffusivity in the ocean through the GCMs and EMICs. This requires the analysis of possible schemes of diapycnal mixing with the ultimate goal of finding the optimal model to predict the evolution of the climate system, the study of all variables that affect it and the correct simulation over long periods of time.

Applicazione di strumenti e metodi per la gestione ambientale di un territorio vasto - Sviluppo di un sistema di gestione ambientale in un ente pubblico secondo lo schema europeo EMAS e analisi di un sistema energetico territoriale

Environmental Management includes many components, among which we can include Environmental Management Systems (EMS), Environmental Reporting and Analysis, Environmental Information Systems and Environmental Communication. In this work two applications are presented: the developement and implementation of an Environmental Management System in local administrations, according to the European scheme "EMAS", and the analysis of a territorial energy system through scenario building and environmental sustainability assessment. Both applications are linked by the same objective, which is the quest for more scientifically sound elements; in fact, both EMS and energy planning are oftec carachterized by localism and poor comparability. Emergy synthesis, proposed by ecologist H.T. Odum and described in his book "Environmental Accounting: Emergy and Environmental Decision Making" (1996) has been chosen and applied as an environmental evaluation tool, in order complete the analysis with an assessment of the "global value" of goods and processes. In particular, eMergy syntesis has been applied in order to improve the evaluation of the significance of environmental aspects in an EMS, and in order to evaluate the environmental performance of three scenarios of future evolution of the energy system. Regarding EMS, in this work an application of an EMS together with the CLEAR methodology for environmental accounting is discussed, in order to improve the identification of the environmental aspects; data regarding environmental aspects and significant ones for 4 local authorities are also presented, together with a preliminary proposal for the integration of the assessment of the significance of environmental aspects with eMergy synthesis. Regarding the analysis of an energy system, in this work the carachterization of the current situation is presented together with the overall energy balance and the evaluation of the emissions of greenhouse gases; moreover, three scenarios of future evolution are described and discussed. The scenarios have been realized with the support of the LEAP software ("Long Term Energy Alternatives Planning System" by SEI - "Stockholm Environment Institute"). Finally, the eMergy synthesis of the current situation and of the three scenarios is shown.

NMR spectroscopy and imaging of hyperpolarized gases

Since the discovery of the nuclear magnetic resonance (NMR) phenomenon, countless NMR techniques have been developed that are today indispensable tools in physics, chemistry, biology, and medicine. As one of the main obstacles in NMR is its notorious lack of sensitivity, different hyperpolarization (HP) methods have been established to increase signals up to several orders of magnitude. In this work, different aspects of magnetic resonance, using HP noble gases, are studied, hereby combining different disciplines of research. The first part examines new fundamental effects in NMR of HP gases, in theory and experiment. The spin echo phenomenon, which provides the basis of numerous modern experiments, is studied in detail in the gas phase. The changes of the echo signal in terms of amplitude, shape, and position, due to the fast translational motion, are described by an extension of the existing theory and computer simulations. With this knowledge as a prerequisite, the detection of intermolecular double-quantum coherences was accomplished for the first time in the gas phase. The second part of this thesis focuses on the development of a practical method to enhance the dissolution process of HP 129Xe, without loss of polarization or shortening of T1. Two different setups for application in NMR spectroscopy and magnetic resonance imaging (MRI) are presented. The continuous operation allows biological and multidimensional spectroscopy in solutions. Also, first in vitro MRI images with dissolved HP 129Xe as contrast agent were obtained at a clinical scanner.

Interactions between the hydrological cycle and atmospheric chemistry

The land-atmosphere exchange of atmospheric trace gases is sensitive to meteorological conditions and climate change. It contributes in turn to the atmospheric radiative forcing through its effects on tropospheric chemistry. The interactions between the hydrological cycle and atmospheric processes are intricate and often involve different levels of feedbacks. The Earth system model EMAC is used in this thesis to assess the direct role of the land surface components of the terrestrial hydrological cycle in the emissions, deposition and transport of key trace gases that control tropospheric chemistry. It is also used to examine its indirect role in changing the tropospheric chemical composition through the feedbacks between the atmospheric and the terrestrial branches of the hydrological cycle. Selected features of the hydrological cycle in EMAC are evaluated using observations from different data sources. The interactions between precipitation and the water vapor column, from the atmospheric branch of the hydrological cycle, and evapotranspiration, from its terrestrial branch, are assessed specially for tropical regions. The impacts of changes in the land surface hydrology on surface exchanges and the oxidizing chemistry of the atmosphere are assessed through two sensitivity simulations. In the first, a new parametrization for rainfall interception in the densely vegetated areas in the tropics is implemented, and its effects are assessed. The second study involves the application of a soil moisture forcing that replaces the model calculated soil moisture. Both experiments have a large impact on the local hydrological cycle, dry deposition of soluble and insoluble gases, emissions of isoprene through changes in surface temperature and the Planetary Boundary Layer height. Additionally the soil moisture forcing causes changes in local vertical transport and large-scale circulation. The changes in trace gas exchanges affect the oxidation capacity of the atmosphere through changes in OH, O$_3$, NO$_x$ concentrations.

A farm-level programming model to compare the atmospheric impact of conventional and organic farming

A model is developed to represent the activity of a farm using the method of linear programming. Two are the main components of the model, the balance of soil fertility and the livestock nutrition. According to the first, the farm is supposed to have a total requirement of nitrogen, which is to be accomplished either through internal sources (manure) or through external sources (fertilisers). The second component describes the animal husbandry as having a nutritional requirement which must be satisfied through the internal production of arable crops or the acquisition of feed from the market. The farmer is supposed to maximise total net income from the agricultural and the zoo-technical activities by choosing one rotation among those available for climate and acclivity. The perspective of the analysis is one of a short period: the structure of the farm is supposed to be fixed without possibility to change the allocation of permanent crops and the amount of animal husbandry. The model is integrated with an environmental module that describes the role of the farm within the carbon-nitrogen cycle. On the one hand the farm allows storing carbon through the photosynthesis of the plants and the accumulation of carbon in the soil; on the other some activities of the farm emit greenhouse gases into the atmosphere. The model is tested for some representative farms of the Emilia-Romagna region, showing to be capable to give different results for conventional and organic farming and providing first results concerning the different atmospheric impact. Relevant data about the representative farms and the feasible rotations are extracted from the FADN database, with an integration of the coefficients from the literature.

Teoria e pratica dell'architettura solare - Morfologia, rendimento, strategia progettuale

La recente Direttiva 31/2010 dell’Unione Europea impone agli stati membri di riorganizzare il quadro legislativo nazionale in materia di prestazione energetica degli edifici, affinchè tutte le nuove costruzioni presentino dal 1° gennaio 2021 un bilancio energetico tendente allo zero; termine peraltro anticipato al 1° gennaio 2019 per gli edifici pubblici. La concezione di edifici a energia “quasi” zero (nZEB) parte dal presupposto di un involucro energeticamente di standard passivo per arrivare a compensare, attraverso la produzione preferibilmente in sito di energia da fonti rinnovabili, gli esigui consumi richiesti su base annuale. In quest’ottica la riconsiderazione delle potenzialità dell’architettura solare individua degli strumenti concreti e delle valide metodologie per supportare la progettazione di involucri sempre più performanti che sfruttino pienamente una risorsa inesauribile, diffusa e alla portata di tutti come quella solare. Tutto ciò in considerazione anche della non più procrastinabile necessità di ridurre il carico energetico imputabile agli edifici, responsabili come noto di oltre il 40% dei consumi mondiali e del 24% delle emissioni di gas climalteranti. Secondo queste premesse la ricerca pone come centrale il tema dell’integrazione dei sistemi di guadagno termico, cosiddetti passivi, e di produzione energetica, cosiddetti attivi, da fonte solare nell’involucro architettonico. Il percorso sia analitico che operativo effettuato si è posto la finalità di fornire degli strumenti metodologici e pratici al progetto dell’architettura, bisognoso di un nuovo approccio integrato mirato al raggiungimento degli obiettivi di risparmio energetico. Attraverso una ricognizione generale del concetto di architettura solare e dei presupposti teorici e terminologici che stanno alla base della stessa, la ricerca ha prefigurato tre tipologie di esito finale: una codificazione delle morfologie ricorrenti nelle realizzazioni solari, un’analisi comparata del rendimento solare nelle principali aggregazioni tipologiche edilizie e una parte importante di verifica progettuale dove sono stati applicati gli assunti delle categorie precedenti

Quality assessment of a landslide inventory map and its application to land‐use planning. A case study in the Northern Apennines (Emilia‐Romagna region, Italy)

Landslide hazard and risk are growing as a consequence of climate change and demographic pressure. Land‐use planning represents a powerful tool to manage this socio‐economic problem and build sustainable and landslide resilient communities. Landslide inventory maps are a cornerstone of land‐use planning and, consequently, their quality assessment represents a burning issue. This work aimed to define the quality parameters of a landslide inventory and assess its spatial and temporal accuracy with regard to its possible applications to land‐use planning. In this sense, I proceeded according to a two‐steps approach. An overall assessment of the accuracy of data geographic positioning was performed on four case study sites located in the Italian Northern Apennines. The quantification of the overall spatial and temporal accuracy, instead, focused on the Dorgola Valley (Province of Reggio Emilia). The assessment of spatial accuracy involved a comparison between remotely sensed and field survey data, as well as an innovative fuzzylike analysis of a multi‐temporal landslide inventory map. Conversely, long‐ and short‐term landslide temporal persistence was appraised over a period of 60 years with the aid of 18 remotely sensed image sets. These results were eventually compared with the current Territorial Plan for Provincial Coordination (PTCP) of the Province of Reggio Emilia. The outcome of this work suggested that geomorphologically detected and mapped landslides are a significant approximation of a more complex reality. In order to convey to the end‐users this intrinsic uncertainty, a new form of cartographic representation is needed. In this sense, a fuzzy raster landslide map may be an option. With regard to land‐use planning, landslide inventory maps, if appropriately updated, confirmed to be essential decision‐support tools. This research, however, proved that their spatial and temporal uncertainty discourages any direct use as zoning maps, especially when zoning itself is associated to statutory or advisory regulations.

Design of an indirectly fired gas turbine integrated with an organic rankine cycle unit for combined heat and power production

In the last years, the European countries have paid increasing attention to renewable sources and greenhouse emissions. The Council of the European Union and the European Parliament have established ambitious targets for the next years. In this scenario, biomass plays a prominent role since its life cycle produces a zero net carbon dioxide emission. Additionally, biomass can ensure plant operation continuity thanks to its availability and storage ability. Several conventional systems running on biomass are available at the moment. Most of them are performant either in the large-scale or in the small power range. The absence of an efficient system on the small-middle scale inspired this thesis project. The object is an innovative plant based on a wet indirectly fired gas turbine (WIFGT) integrated with an organic Rankine cycle (ORC) unit for combined heat and power production. The WIFGT is a performant system in the small-middle power range; the ORC cycle is capable of giving value to low-temperature heat sources. Their integration is investigated in this thesis with the aim of carrying out a preliminary design of the components. The targeted plant output is around 200 kW in order not to need a wide cultivation area and to avoid biomass shipping. Existing in-house simulation tools are used: They are adapted to this purpose. Firstly the WIFGT + ORC model is built; Zero-dimensional models of heat exchangers, compressor, turbines, furnace, dryer and pump are used. Different fluids are selected but toluene and benzene turn out to be the most suitable. In the indirectly fired gas turbine a pressure ratio around 4 leads to the highest efficiency. From the thermodynamic analysis the system shows an electric efficiency of 38%, outdoing other conventional plants in the same power range. The combined plant is designed to recover thermal energy: Water is used as coolant in the condenser. It is heated from 60°C up to 90°C, ensuring the possibility of space heating. Mono-dimensional models are used to design the heat exchange equipment. Different types of heat exchangers are chosen depending on the working temperature. A finned-plate heat exchanger is selected for the WIFGT heat transfer equipment due to the high temperature, oxidizing and corrosive environment. A once-through boiler with finned tubes is chosen to vaporize the organic fluid in the ORC. A plate heat exchanger is chosen for the condenser and recuperator. A quasi-monodimensional model for single-stage axial turbine is implemented to design both the WIFGT and the ORC turbine. The system simulation after the components design shows an electric efficiency around 34% with a decrease by 10% compared to the zero-dimensional analysis. The work exhibits the system potentiality compared to the existing plants from both technical and economic point of view.

Simultaneous stable isotope analysis of methane and nitrous oxide on ice core samples

Methane and nitrous oxide are important greenhouse gases which show a strong increase in atmospheric mixing ratios since pre-industrial time as well as large variations during past climate changes. The understanding of their biogeochemical cycles can be improved using stable isotope analysis. However, high-precision isotope measurements on air trapped in ice cores are challenging because of the high susceptibility to contamination and fractionation. Here, we present a dry extraction system for combined CH4 and N2O stable isotope analysis from ice core air, using an ice grating device. The system allows simultaneous analysis of δD(CH4) or δ13C(CH4), together with δ15N(N2O), δ18O(N2O) and δ15N(NO+ fragment) on a single ice core sample, using two isotope mass spectrometry systems. The optimum quantity of ice for analysis is about 600 g with typical "Holocene" mixing ratios for CH4 and N2O. In this case, the reproducibility (1σ ) is 2.1‰ for δD(CH4), 0.18‰ for δ13C(CH4), 0.51‰ for δ15N(N2O), 0.69‰ for δ18O(N2O) and 1.12‰ for δ15N(NO+ fragment). For smaller amounts of ice the standard deviation increases, particularly for N2O isotopologues. For both gases, small-scale intercalibrations using air and/or ice samples have been carried out in collaboration with other institutes that are currently involved in isotope measurements of ice core air. Significant differences are shown between the calibration scales, but those offsets are consistent and can therefore be corrected for.

Extended frequency range hearing thresholds and otoacoustic emissions in acute acoustic trauma

We sought to evaluate the relative value of pure tone audiometry (PTA), extended high-frequency audiometry (EFA) and transiently evoked otoacoustic emissions (OAE) and distortion products when monitoring acute acoustic trauma (AAT).

Modeling of diesel particulate filter filtration and regeneration for transient driving schedules

Due to their high thermal efficiency, diesel engines have excellent fuel economy and have been widely used as a power source for many vehicles. Diesel engines emit less greenhouse gases (carbon dioxide) compared with gasoline engines. However, diesel engines emit large amounts of particulate matter (PM) which can imperil human health. The best way to reduce the particulate matter is by using the Diesel Particulate Filter (DPF) system which consists of a wall-flow monolith which can trap particulates, and the DPF can be periodically regenerated to remove the collected particulates. The estimation of the PM mass accumulated in the DPF and total pressure drop across the filter are very important in order to determine when to carry out the active regeneration for the DPF. In this project, by developing a filtration model and a pressure drop model, we can estimate the PM mass and the total pressure drop, then, these two models can be linked with a regeneration model which has been developed previously to predict when to regenerate the filter. There results of this project were: 1 Reproduce a filtration model and simulate the processes of filtration. By studying the deep bed filtration and cake filtration, stages and quantity of mass accumulated in the DPF can be estimated. It was found that the filtration efficiency increases faster during the deep-bed filtration than that during the cake filtration. A “unit collector” theory was used in our filtration model which can explain the mechanism of the filtration very well. 2 Perform a parametric study on the pressure drop model for changes in engine exhaust flow rate, deposit layer thickness, and inlet temperature. It was found that there are five primary variables impacting the pressure drop in the DPF which are temperature gradient along the channel, deposit layer thickness, deposit layer permeability, wall thickness, and wall permeability. 3 Link the filtration model and the pressure drop model with the regeneration model to determine the time to carry out the regeneration of the DPF. It was found that the regeneration should be initiated when the cake layer is at a certain thickness, since a cake layer with either too big or too small an amount of particulates will need more thermal energy to reach a higher regeneration efficiency. 4 Formulate diesel particulate trap regeneration strategies for real world driving conditions to find out the best desirable conditions for DPF regeneration. It was found that the regeneration should be initiated when the vehicle’s speed is high and during which there should not be any stops from the vehicle. Moreover, the regeneration duration is about 120 seconds and the inlet temperature for the regeneration is 710K.

Development of an optimization model for biofuel facility size and location and a simulation model for design of a biofuel supply chain

To mitigate greenhouse gas (GHG) emissions and reduce U.S. dependence on imported oil, the United States (U.S.) is pursuing several options to create biofuels from renewable woody biomass (hereafter referred to as “biomass”). Because of the distributed nature of biomass feedstock, the cost and complexity of biomass recovery operations has significant challenges that hinder increased biomass utilization for energy production. To facilitate the exploration of a wide variety of conditions that promise profitable biomass utilization and tapping unused forest residues, it is proposed to develop biofuel supply chain models based on optimization and simulation approaches. The biofuel supply chain is structured around four components: biofuel facility locations and sizes, biomass harvesting/forwarding, transportation, and storage. A Geographic Information System (GIS) based approach is proposed as a first step for selecting potential facility locations for biofuel production from forest biomass based on a set of evaluation criteria, such as accessibility to biomass, railway/road transportation network, water body and workforce. The development of optimization and simulation models is also proposed. The results of the models will be used to determine (1) the number, location, and size of the biofuel facilities, and (2) the amounts of biomass to be transported between the harvesting areas and the biofuel facilities over a 20-year timeframe. The multi-criteria objective is to minimize the weighted sum of the delivered feedstock cost, energy consumption, and GHG emissions simultaneously. Finally, a series of sensitivity analyses will be conducted to identify the sensitivity of the decisions, such as the optimal site selected for the biofuel facility, to changes in influential parameters, such as biomass availability and transportation fuel price. Intellectual Merit The proposed research will facilitate the exploration of a wide variety of conditions that promise profitable biomass utilization in the renewable biofuel industry. The GIS-based facility location analysis considers a series of factors which have not been considered simultaneously in previous research. Location analysis is critical to the financial success of producing biofuel. The modeling of woody biomass supply chains using both optimization and simulation, combing with the GIS-based approach as a precursor, have not been done to date. The optimization and simulation models can help to ensure the economic and environmental viability and sustainability of the entire biofuel supply chain at both the strategic design level and the operational planning level. Broader Impacts The proposed models for biorefineries can be applied to other types of manufacturing or processing operations using biomass. This is because the biomass feedstock supply chain is similar, if not the same, for biorefineries, biomass fired or co-fired power plants, or torrefaction/pelletization operations. Additionally, the research results of this research will continue to be disseminated internationally through publications in journals, such as Biomass and Bioenergy, and Renewable Energy, and presentations at conferences, such as the 2011 Industrial Engineering Research Conference. For example, part of the research work related to biofuel facility identification has been published: Zhang, Johnson and Sutherland [2011] (see Appendix A). There will also be opportunities for the Michigan Tech campus community to learn about the research through the Sustainable Future Institute.