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.

DEVELOPMENT OF AN OPTIMIZATION MODEL FOR BIOFUEL FACILITY SIZE AND LOCATION AND A SIMULATION MODEL FOR DESIGN OF A BIOFUEL SUPPLY CHAIN

A range of societal issues have been caused by fossil fuel consumption in the transportation sector in the United States (U.S.), including health related air pollution, climate change, the dependence on imported oil, and other oil related national security concerns. Biofuels production from various lignocellulosic biomass types such as wood, forest residues, and agriculture residues have the potential to replace a substantial portion of the total fossil fuel consumption. This research focuses on locating biofuel facilities and designing the biofuel supply chain to minimize the overall cost. For this purpose an integrated methodology was proposed by combining the GIS technology with simulation and optimization modeling methods. The GIS based methodology was used as a precursor for selecting biofuel facility locations by employing a series of decision factors. The resulted candidate sites for biofuel production served as inputs for simulation and optimization modeling. As a precursor to simulation or optimization modeling, the GIS-based methodology was used to preselect potential biofuel facility locations for biofuel production from forest biomass. Candidate locations were selected based on a set of evaluation criteria, including: county boundaries, a railroad transportation network, a state/federal road transportation network, water body (rivers, lakes, etc.) dispersion, city and village dispersion, a population census, biomass production, and no co-location with co-fired power plants. The simulation and optimization models were built around key supply activities including biomass harvesting/forwarding, transportation and storage. The built onsite storage served for spring breakup period where road restrictions were in place and truck transportation on certain roads was limited. Both models were evaluated using multiple performance indicators, including cost (consisting of the delivered feedstock cost, and inventory holding cost), energy consumption, and GHG emissions. The impact of energy consumption and GHG emissions were expressed in monetary terms to keep consistent with cost. Compared with the optimization model, the simulation model represents a more dynamic look at a 20-year operation by considering the impacts associated with building inventory at the biorefinery to address the limited availability of biomass feedstock during the spring breakup period. The number of trucks required per day was estimated and the inventory level all year around was tracked. Through the exchange of information across different procedures (harvesting, transportation, and biomass feedstock processing procedures), a smooth flow of biomass from harvesting areas to a biofuel facility was implemented. The optimization model was developed to address issues related to locating multiple biofuel facilities simultaneously. The size of the potential biofuel facility is set up with an upper bound of 50 MGY and a lower bound of 30 MGY. The optimization model is a static, Mathematical Programming Language (MPL)-based application which allows for sensitivity analysis by changing inputs to evaluate different scenarios. It was found that annual biofuel demand and biomass availability impacts the optimal results of biofuel facility locations and sizes.

ANOVA kernels and RKHS of zero mean functions for model-based sensitivity analysis

Given a reproducing kernel Hilbert space (H,〈.,.〉)(H,〈.,.〉) of real-valued functions and a suitable measure μμ over the source space D⊂RD⊂R, we decompose HH as the sum of a subspace of centered functions for μμ and its orthogonal in HH. This decomposition leads to a special case of ANOVA kernels, for which the functional ANOVA representation of the best predictor can be elegantly derived, either in an interpolation or regularization framework. The proposed kernels appear to be particularly convenient for analyzing the effect of each (group of) variable(s) and computing sensitivity indices without recursivity.

Investigating the sensitivity of hurricane intensity and trajectory to sea surface temperatures using the regional model WRF

The influence of sea surface temperature (SST) anomalies on the hurricane characteristics are investigated in a set of sensitivity experiments employing the Weather Research and Forecasting (WRF) model. The idealised experiments are performed for the case of Hurricane Katrina in 2005. The first set of sensitivity experiments with basin-wide changes of the SST magnitude shows that the intensity goes along with changes in the SST, i.e., an increase in SST leads to an intensification of Katrina. Additionally, the trajectory is shifted to the west (east), with increasing (decreasing) SSTs. The main reason is a strengthening of the background flow. The second set of experiments investigates the influence of Loop Current eddies idealised by localised SST anomalies. The intensity of Hurricane Katrina is enhanced with increasing SSTs close to the core of a tropical cyclone. Negative nearby SST anomalies reduce the intensity. The trajectory only changes if positive SST anomalies are located west or north of the hurricane centre. In this case the hurricane is attracted by the SST anomaly which causes an additional moisture source and increased vertical winds.

Regional climate model simulations for Europe at 6 and 0.2 k BP: sensitivity to changes in anthropogenic deforestation

This study aims to evaluate the direct effects of anthropogenic deforestation on simulated climate at two contrasting periods in the Holocene, ~6 and ~0.2 k BP in Europe. We apply We apply the Rossby Centre regional climate model RCA3, a regional climate model with 50 km spatial resolution, for both time periods, considering three alternative descriptions of the past vegetation: (i) potential natural vegetation (V) simulated by the dynamic vegetation model LPJ-GUESS, (ii) potential vegetation with anthropogenic land use (deforestation) from the HYDE3.1 (History Database of the Global Environment) scenario (V + H3.1), and (iii) potential vegetation with anthropogenic land use from the KK10 scenario (V + KK10). The climate model results show that the simulated effects of deforestation depend on both local/regional climate and vegetation characteristics. At ~6 k BP the extent of simulated deforestation in Europe is generally small, but there are areas where deforestation is large enough to produce significant differences in summer temperatures of 0.5–1 °C. At ~0.2 k BP, extensive deforestation, particularly according to the KK10 model, leads to significant temperature differences in large parts of Europe in both winter and summer. In winter, deforestation leads to lower temperatures because of the differences in albedo between forested and unforested areas, particularly in the snow-covered regions. In summer, deforestation leads to higher temperatures in central and eastern Europe because evapotranspiration from unforested areas is lower than from forests. Summer evaporation is already limited in the southernmost parts of Europe under potential vegetation conditions and, therefore, cannot become much lower. Accordingly, the albedo effect dominates in southern Europe also in summer, which implies that deforestation causes a decrease in temperatures. Differences in summer temperature due to deforestation range from −1 °C in south-western Europe to +1 °C in eastern Europe. The choice of anthropogenic land-cover scenario has a significant influence on the simulated climate, but uncertainties in palaeoclimate proxy data for the two time periods do not allow for a definitive discrimination among climate model results.

On the reduced sensitivity of the Atlantic overturning to Greenland ice sheet melting in projections: a multi-model assessment

Large uncertainties exist concerning the impact of Greenland ice sheet melting on the Atlantic meridional overturning circulation (AMOC) in the future, partly due to different sensitivity of the AMOC to freshwater input in the North Atlantic among climate models. Here we analyse five projections from different coupled ocean–atmosphere models with an additional 0.1 Sv (1 Sv = 10 6 m3/s) of freshwater released around Greenland between 2050 and 2089. We find on average a further weakening of the AMOC at 26°N of 1.1 ± 0.6 Sv representing a 27 ± 14% supplementary weakening in 2080–2089, as compared to the weakening relative to 2006–2015 due to the effect of the external forcing only. This weakening is lower than what has been found with the same ensemble of models in an identical experimen - tal set-up but under recent historical climate conditions. This lower sensitivity in a warmer world is explained by two main factors. First, a tendency of decoupling is detected between the surface and the deep ocean caused by an increased thermal stratification in the North Atlantic under the effect of global warming. This induces a shoaling of ocean deep ventilation through convection hence ventilating only intermediate levels. The second important effect concerns the so-called Canary Current freshwater leakage; a process by which additionally released fresh water in the North Atlantic leaks along the Canary Current and escapes the convection zones towards the subtropical area. This leakage is increasing in a warming climate, which is a consequence of decreasing gyres asymmetry due to changes in Ekman rumping. We suggest that these modifications are related with the northward shift of the jet stream in a warmer world. For these two reasons the AMOC is less susceptible to freshwater perturbations (near the deep water formation sides) in the North Atlantic as compared to the recent historical climate conditions. Finally, we propose a bilinear model that accounts for the two former processes to give a conceptual explanation about the decreasing AMOC sensitivity due to freshwater input. Within the limit of this bilinear model, we find that 62 ± 8% of the reduction in sensitivity is related with the changes in gyre asymmetry and freshwater leakage and 38 ± 8% is due to the reduction in deep ocean ventilation associated with the increased stratification in the North Atlantic.

Sensitivity of the WRF model to PBL parametrisations and nesting techniques: evaluation of wind storms over complex terrain

Simulating surface wind over complex terrain is a challenge in regional climate modelling. Therefore, this study aims at identifying a set-up of the Weather Research and Forecasting Model (WRF) model that minimises system- atic errors of surface winds in hindcast simulations. Major factors of the model configuration are tested to find a suitable set-up: the horizontal resolution, the planetary boundary layer (PBL) parameterisation scheme and the way the WRF is nested to the driving data set. Hence, a number of sensitivity simulations at a spatial resolution of 2 km are carried out and compared to observations. Given the importance of wind storms, the analysis is based on case studies of 24 historical wind storms that caused great economic damage in Switzerland. Each of these events is downscaled using eight different model set-ups, but sharing the same driving data set. The results show that the lack of representation of the unresolved topography leads to a general overestimation of wind speed in WRF. However, this bias can be substantially reduced by using a PBL scheme that explicitly considers the effects of non-resolved topography, which also improves the spatial structure of wind speed over Switzerland. The wind direction, although generally well reproduced, is not very sensitive to the PBL scheme. Further sensitivity tests include four types of nesting methods: nesting only at the boundaries of the outermost domain, analysis nudging, spectral nudging, and the so-called re-forecast method, where the simulation is frequently restarted. These simulations show that restricting the freedom of the model to develop large-scale disturbances slightly increases the temporal agreement with the observations, at the same time that it further reduces the overestimation of wind speed, especially for maximum wind peaks. The model performance is also evaluated in the outermost domains, where the resolution is coarser. The results demonstrate the important role of horizontal resolution, where the step from 6 to 2 km significantly improves model performance. In summary, the combination of a grid size of 2 km, the non-local PBL scheme modified to explicitly account for non-resolved orography, as well as analysis or spectral nudging, is a superior combination when dynamical downscaling is aimed at reproducing real wind fields.

Sleep-Wake Cycle Dysfunction in the TgCRND8 Mouse Model of Alzheimer's Disease: From Early to Advanced Pathological Stages.

In addition to cognitive decline, individuals affected by Alzheimer's disease (AD) can experience important neuropsychiatric symptoms including sleep disturbances. We characterized the sleep-wake cycle in the TgCRND8 mouse model of AD, which overexpresses a mutant human form of amyloid precursor protein resulting in high levels of β-amyloid and plaque formation by 3 months of age. Polysomnographic recordings in freely-moving mice were conducted to study sleep-wake cycle architecture at 3, 7 and 11 months of age and corresponding levels of β-amyloid in brain regions regulating sleep-wake states were measured. At all ages, TgCRND8 mice showed increased wakefulness and reduced non-rapid eye movement (NREM) sleep during the resting and active phases. Increased wakefulness in TgCRND8 mice was accompanied by a shift in the waking power spectrum towards fast frequency oscillations in the beta (14-20 Hz) and low gamma range (20-50 Hz). Given the phenotype of hyperarousal observed in TgCRND8 mice, the role of noradrenergic transmission in the promotion of arousal, and previous work reporting an early disruption of the noradrenergic system in TgCRND8, we tested the effects of the alpha-1-adrenoreceptor antagonist, prazosin, on sleep-wake patterns in TgCRND8 and non-transgenic (NTg) mice. We found that a lower dose (2 mg/kg) of prazosin increased NREM sleep in NTg but not in TgCRND8 mice, whereas a higher dose (5 mg/kg) increased NREM sleep in both genotypes, suggesting altered sensitivity to noradrenergic blockade in TgCRND8 mice. Collectively our results demonstrate that amyloidosis in TgCRND8 mice is associated with sleep-wake cycle dysfunction, characterized by hyperarousal, validating this model as a tool towards understanding the relationship between β-amyloid overproduction and disrupted sleep-wake patterns in AD.

A Model Study of the Circulation in the Pearl River Estuary (PRE) and Its Adjacent Coastal Waters: 2. Sensitivity Experiments

The Princeton Ocean Model is used to study the circulation in the Pear River Estuary (PRE) and the adjacent coastal waters in the winter and summer seasons. Wong et al. [2003] compares the simulation results with the in situ measurements collected during the Pearl River Estuary Pollution Project (PREPP). In this paper, sensitivity experiments are carried out to examine the plume and the associated frontal dynamics in response to seasonal discharges and monsoon winds. During the winter, convergence between the seaward spreading plume water and the saline coastal water sets up a salinity front that aligns from the northeast to the southwest inside the PRE. During the summer the plume water fills the PRE at the surface and spreads eastward in the coastal waters in response to the prevailing southwesterly monsoon. The overall alignment of the plume is from the northwest to the southeast. The subsurface front is similar to that in the winter and summer except that the summer front is closer to the mouth and the winter front closer to the head of the estuary. Inside the PRE, bottom flows are always toward the head of the estuary, attributed to the density gradient associated with the plume front. In contrast, bottom flows in the shelf change from offshore in winter to onshore in summer, reflecting respectively the wintertime downwelling and summertime upwelling. Wind also plays an essential role in controlling the plume at the surface. An easterly wind drives the plume westward regardless winter or summer. The eastward spreading of the plume during the summer can be attributed to the southerly component of the wind. On the other hand, the surface area of the plume is positively proportional to the amount of discharge.

Assessment of the effect on statistical power of regression model misspecification by using techniques of mathematical statistics and simulation study

Objectives. This paper seeks to assess the effect on statistical power of regression model misspecification in a variety of situations. ^ Methods and results. The effect of misspecification in regression can be approximated by evaluating the correlation between the correct specification and the misspecification of the outcome variable (Harris 2010).In this paper, three misspecified models (linear, categorical and fractional polynomial) were considered. In the first section, the mathematical method of calculating the correlation between correct and misspecified models with simple mathematical forms was derived and demonstrated. In the second section, data from the National Health and Nutrition Examination Survey (NHANES 2007-2008) were used to examine such correlations. Our study shows that comparing to linear or categorical models, the fractional polynomial models, with the higher correlations, provided a better approximation of the true relationship, which was illustrated by LOESS regression. In the third section, we present the results of simulation studies that demonstrate overall misspecification in regression can produce marked decreases in power with small sample sizes. However, the categorical model had greatest power, ranging from 0.877 to 0.936 depending on sample size and outcome variable used. The power of fractional polynomial model was close to that of linear model, which ranged from 0.69 to 0.83, and appeared to be affected by the increased degrees of freedom of this model.^ Conclusion. Correlations between alternative model specifications can be used to provide a good approximation of the effect on statistical power of misspecification when the sample size is large. When model specifications have known simple mathematical forms, such correlations can be calculated mathematically. Actual public health data from NHANES 2007-2008 were used as examples to demonstrate the situations with unknown or complex correct model specification. Simulation of power for misspecified models confirmed the results based on correlation methods but also illustrated the effect of model degrees of freedom on power.^

Case studies of the effects of a total solar eclipse on weather with varying degrees of cloud cover using a local numerical weather prediction model

The importance of renewable energies for the European electricity market is growing rapidly. This presents transmission grids and the power market in general with new challenges which stem from the higher spatiotemporal variability of power generation. This uncertainty is due to the fact that renewable power production results from weather phenomena, thus making it difficult to plan and control. We present a sensitivity study of a total solar eclipse in central Europe in March. The weather in Germany and Europe was modeled using the German Weather Service's local area models COSMO-DE and COSMO-EU, respectively (http://www.cosmo-model.org/). The simulations were performed with and without considering a solar eclipse for the following 3 situations: 1. An idealized, clear-sky situation for the entire model area (Europe, COSMO-EU) 2. A real weather situation with mostly cloudy skies (Germany, COSMO-DE) 3. A real weather situation with mostly clear skies (Germany, COSMO-DE) The data should help to evaluate the effects of a total solar eclipse on the weather in the planetary boundary layer. The results show that a total solar eclipse has significant effects particularly on the main variables for renewable energy production, such as solar irradiation and temperature near the ground.

Modeling the sensitivity of precipitation oxygen isotopes to the Last Glacial Maximum boundary conditions: A study using Community Atmosphere Model (CAM3.0)

To understand the validity of d18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM). A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM) values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST), and orbital parameters) were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (d18Oprecip) in response to individual climate factors. The change in topography (due to the change in land ice cover) played a significant role in reducing the surface temperature and d18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and d18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3). Large reductions in d18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the d18Oprecip distribution among the simulations.

Model results of sensitivity experiments for marine nitrogen cycle in four NetCDF files

The marine nitrogen (N) inventory is thought to be stabilized by negative feedback mechanisms that reduce N inventory excursions relative to the more slowly overturning phosphorus inventory. Using a global biogeochemical ocean circulation model we show that negative feedbacks stabilizing the N inventory cannot persist if a close spatial association of N2 fixation and denitrification occurs. In our idealized model experiments, nitrogen deficient waters, generated by denitrification, stimulate local N2 fixation activity. But, because of stoichiometric constraints, the denitrification of newly fixed nitrogen leads to a net loss of N. This can enhance the N deficit, thereby triggering additional fixation in a vicious cycle, ultimately leading to a runaway N loss. To break this vicious cycle, and allow for stabilizing negative feedbacks to occur, inputs of new N need to be spatially decoupled from denitrification. Our idealized model experiments suggest that factors such as iron limitation or dissolved organic matter cycling can promote such decoupling and allow for negative feedbacks that stabilize the N inventory. Conversely, close spatial co-location of N2 fixation and denitrification could lead to net N loss.

Power output fluctuations in large scale PV plants: one year observations with one second resolution and a derived analytic model

The variable nature of the irradiance can produce significant fluctuations in the power generated by large grid-connected photovoltaic (PV) plants. Experimental 1 s data were collected throughout a year from six PV plants, 18 MWp in total. Then, the dependence of short (below 10 min) power fluctuation on PV plant size has been investigated. The analysis focuses on the study of fluctuation frequency as well as the maximum fluctuation value registered. An analytic model able to describe the frequency of a given fluctuation for a certain day is proposed

Sensitivity of a distributed temperature-radiation index melt model based on AWS observations and surface energy balance fluxes, Hurd Peninsula glaciers, Livingston Island, Antarctica

We use an automatic weather station and surface mass balance dataset spanning four melt seasons collected on Hurd Peninsula Glaciers, South Shetland Islands, to investigate the point surface energy balance, to determine the absolute and relative contribution of the various energy fluxes acting on the glacier surface and to estimate the sensitivity of melt to ambient temperature changes. Long-wave incoming radiation is the main energy source for melt, while short-wave radiation is the most important flux controlling the variation of both seasonal and daily mean surface energy balance. Short-wave and long-wave radiation fluxes do, in general, balance each other, resulting in a high correspondence between daily mean net radiation flux and available melt energy flux. We calibrate a distributed melt model driven by air temperature and an expression for the incoming short-wave radiation. The model is calibrated with the data from one of the melt seasons and validated with the data of the three remaining seasons. The model results deviate at most 140 mm w.e. from the corresponding observations using the glaciological method. The model is very sensitive to changes in ambient temperature: a 0.5 ◦ C increase results in 56 % higher melt rates.