Analysis of atmospheric transport and the effects of seasonal and interannual transport variability on measurements made at the Pico Mountain observatory

A non-hierarchical K-means algorithm is used to cluster 47 years (1960–2006) of 10-day HYSPLIT backward trajectories to the Pico Mountain (PM) observatory on a seasonal basis. The resulting cluster centers identify the major transport pathways and collectively comprise a long-term climatology of transport to the observatory. The transport climatology improves our ability to interpret the observations made there and our understanding of pollution source regions to the station and the central North Atlantic region. I determine which pathways dominate transport to the observatory and examine the impacts of these transport patterns on the O3, NOy, NOx, and CO measurements made there during 2001–2006. Transport from the U.S., Canada, and the Atlantic most frequently reaches the station, but Europe, east Africa, and the Pacific can also contribute significantly depending on the season. Transport from Canada was correlated with the North Atlantic Oscillation (NAO) in spring and winter, and transport from the Pacific was uncorrelated with the NAO. The highest CO and O3 are observed during spring. Summer is also characterized by high CO and O3 and the highest NOy and NOx of any season. Previous studies at the station attributed the summer time high CO and O3 to transport of boreal wildfire emissions (for 2002–2004), and boreal fires continued to affect the station during 2005 and 2006. The particle dispersion model FLEXPART was used to calculate anthropogenic and biomass-burning CO tracer values at the station in an attempt to identify the regions responsible for the high CO and O3 observations during spring and biomass-burning impacts in summer.

Cladophora, mass transport, and the nearshore phosphorus shunt

Excessive Cladophora growth in the Great Lakes has led to beach fouling and the temporary closure of nuclear power plants and has been associated with avian botulism and the persistence of human pathogens. As the growth-limiting nutrient for Cladophora, phosphorus is the appropriate target for management efforts. Dreissenids (zebra and quagga mussels) have the ability to capture particulate phase phosphorus (otherwise unavailable to Cladophora) and release it in a soluble, available form. The significance of this potential nutrient source is, in part, influenced by the interplay between phosphorus flux from the mussel bed and turbulent mixing in establishing the phosphorus levels to which Cladophora is exposed. It is hypothesized that under quiescent conditions phosphorus will accumulate near the sediment-water interface, setting up vertical phosphorus gradients and favorable conditions for resource delivery to Cladophora. These gradients would be eliminated under conditions of wind mixing, reducing the significance of the dreissenid-mediated nutrient contribution. Soluble reactive phosphorus (SRP) levels were monitored over dreissenid beds (densities on the order of 350•m-2 and 3000∙m-2) at a site 8 m deep in Lake Michigan. Monitoring was based on the deployment of Modified Hesslein Samplers which collected samples for SRP analysis over a distance of 34 cm above the bottom in 2.5 cm intervals. Deployment intervals were established to capture a wind regime (calm, windy) that persisted for an interval consistent with the sampler equilibration time of 7 hours. Results indicate that increased mussel density leads to an increased concentration boundary layer; increased wind speed leads to entrainment of the concentration boundary layer; and increased duration of quiescent periods leads to an increased concentration boundary layer. This concentration boundary layer is of ecological significance and forms in the region inhabited by Cladophora

IMPACT OF 2000-2050 CLIMATE CHANGE ON GLOBAL ATMOSPHERIC TRANSPORT AND DEPOSITION OF MERCURY

Since it is very toxic and accumulates in organisms, particularly in fish, mercury is a very important pollutant and one of the most studies. And this concern over the toxicity and human health risks of mercury has prompted efforts to regulate anthropogenic emissions. As mercury pollution problem is getting increasingly serious, we are curious about how serious this problem will be in the future. What is more, how the climate change in the future will affect the mercury concentration in the atmosphere. So we investigate the impact of climate change on mercury concentration in the atmosphere. We focus on the comparison between the mercury data for year 2000 and for year 2050. The GEOS-Chem model shows that the mercury concentrations for all tracers (1 to 3), elemental mercury (Hg(0)), divalent mercury (Hg(II)) and primary particulate mercury (Hg(P)) have differences between 2000 and 2050 in most regions over the world. From the model results, we can see the climate change from 2000 to 2050 would decrease Hg(0) surface concentration in most of the world. The driving factors of Hg(0) surface concentration changes are natural emissions(ocean and vegetation) and the transformation reactions between Hg(0) and Hg(II). The climate change from 2000 to 2050 would increase Hg(II) surface concentration in most of mid-latitude continental parts of the world while decreasing Hg(II) surface concentration in most of high-latitude part of the world. The driving factors of Hg(II) surface concentration changes is deposition amount change (majorly wet deposition) from 2000 to 2050 and the transformation reactions between Hg(0) and Hg(II). Climate change would increase Hg(P) concentration in most of mid-latitude area of the world and meanwhile decrease Hg(P) concentration in most of high-latitude regions of the world. For the Hg(P) concentration changes, the major driving factor is the deposition amount change (mainly wet deposition) from 2000 to 2050.

Study of water use and supply in the district of Independencia, Peru

Peru is a developing country with abundant fresh water resources, yet the lack of infrastructure leaves much of the population without access to safe water for domestic uses. The author of this report was a Peace Corps Volunteer in the sector of water & sanitation in the district of Independencia, Ica, Peru. Independencia is located in the arid coastal region of the country, receiving on average 15 mm of rain annually. The water source for this district comes from the Pisco River, originating in the Andean highlands and outflowing into the Pacific Ocean near the town of Pisco, Peru. The objectives of this report are to assess the water supply and sanitation practices, model the existing water distribution system, and make recommendations for future expansion of the distribution system in the district of Independencia, Peru. The assessment of water supply will be based on the results from community surveys done in the district of Independencia, water quality testing done by a detachment of the U.S. Navy, as well as on the results of a hydraulic model built in EPANET 2.0 to represent the distribution system. Sanitation practice assessments will be based on the surveys as well as observations from the author while living in Peru. Recommendations for system expansions will be made based on results from the EPANET model and the municipality’s technical report for the existing distribution system. Household water use and sanitation surveys were conducted with 84 families in the district revealing that upwards of 85% store their domestic water in regularly washed containers with lids. Over 80% of those surveyed are drinking water that is treated, mostly boiled. Of those surveyed, over 95% reported washing their hands and over 60% mentioned at least one critical time for hand washing when asked for specific instances. From the surveys, it was also discovered that over 80% of houses are properly disposing of excrement, in either latrines or septic tanks. There were 43 families interviewed with children five years of age or under, and just over 18% reported the child had a case of diarrhea within the last month at the time of the interview. Finally, from the surveys it was calculated that the average water use per person per day is about 22 liters. Water quality testing carried out by a detachment of the U.S. Navy revealed that the water intended for consumption in the houses surveyed was not suitable for consumption, with a median E. coli most probable number of 47/100 ml for the 61 houses sampled. The median total coliforms was 3,000 colony forming units per 100 ml. EPANET was used to simulate the water delivery system and evaluate its performance. EPANET is designed for continuous water delivery systems, assuming all pipes are always flowing full. To account for the intermittent nature of the system, multiple EPANET network models were created to simulate how water is routed to the different parts of the system throughout the day. The models were created from interviews with the water technicians and a map of the system created using handheld GPS units. The purpose is to analyze the performance of the water system that services approximately 13,276 people in the district of Independencia, Peru, as well as provide recommendations for future growth and improvement of the service level. Performance evaluation of the existing system is based on meeting 25 liters per person per day while maintaining positive pressure at all nodes in the network. The future performance is based on meeting a minimum pressure of 20 psi in the main line, as proposed by Chase (2000). The EPANET model results yield an average nodal pressure for all communities of 71 psi, with a range from 1.3 – 160 psi. Thus, if the current water delivery schedule obtained from the local municipality is followed, all communities should have sufficient pressure to deliver 25 l/p/d, with the exception of Los Rosales, which can only supply 3.25 l/p/d. However, if the line to Los Rosales were increased from one to four inches, the system could supply this community with 25 l/p/d. The district of Independencia could greatly benefit from increasing the service level to 24-hour water delivery and a minimum of 50 l/p/d, so that communities without reliable access due to insufficient pressure would become equal beneficiaries of this invaluable resource. To evaluate the feasibility of this, EPANET was used to model the system with a range of population growth rates, system lifetimes, and demands. In order to meet a minimum pressure of 20 psi in the main line, the 6-inch diameter main line must be increased and approximately two miles of trench must be excavated up to 30 feet deep. The sections of the main line that must be excavated are mile 0-1 and 1.5-2.5, and the first 3.4 miles of the main line must be increased from 6 to 16 inches, contracting to 10 inches for the remaining 5.8 miles. Doing this would allow 24-hour water delivery and provide 50 l/p/d for a range of population growth rates and system lifetimes. It is expected that improving the water delivery service would reduce the morbidity and mortality from diarrheal diseases by decreasing the recontamination of the water due to transport and household storage, as well as by maintaining continuous pressure in the system to prevent infiltration of contaminated groundwater. However, this expansion must be carefully planned so as not to affect aquatic ecosystems or other districts utilizing water from the Pisco River. It is recommended that stream gaging of the Pisco River and precipitation monitoring of the surrounding watershed is initiated in order to begin a hydrological study that would be integrated into the district’s water resource planning. It is also recommended that the district begin routine water quality testing, with the results available to the public.

Long-Range Pollution Transport: Trans-Atlantic Mechanisms and Lagrangian Modeling Methods

Over the past several decades, it has become apparent that anthropogenic activities have resulted in the large-scale enhancement of the levels of many trace gases throughout the troposphere. More recently, attention has been given to the transport pathway taken by these emissions as they are dispersed throughout the atmosphere. The transport pathway determines the physical characteristics of emissions plumes and therefore plays an important role in the chemical transformations that can occur downwind of source regions. For example, the production of ozone (O3) is strongly dependent upon the transport its precursors undergo. O3 can initially be formed within air masses while still over polluted source regions. These polluted air masses can experience continued O3 production or O3 destruction downwind, depending on the air mass's chemical and transport characteristics. At present, however, there are a number of uncertainties in the relationships between transport and O3 production in the North Atlantic lower free troposphere. The first phase of the study presented here used measurements made at the Pico Mountain observatory and model simulations to determine transport pathways for US emissions to the observatory. The Pico Mountain observatory was established in the summer of 2001 in order to address the need to understand the relationships between transport and O3 production. Measurements from the observatory were analyzed in conjunction with model simulations from the Lagrangian particle dispersion model (LPDM), FLEX-PART, in order to determine the transport pathway for events observed at the Pico Mountain observatory during July 2003. A total of 16 events were observed, 4 of which were analyzed in detail. The transport time for these 16 events varied from 4.5 to 7 days, while the transport altitudes over the ocean ranged from 2-8 km, but were typically less than 3 km. In three of the case studies, eastward advection and transport in a weak warm conveyor belt (WCB) airflow was responsible for the export of North American emissions into the FT, while transport in the FT was governed by easterly winds driven by the Azores/Bermuda High (ABH) and transient northerly lows. In the fourth case study, North American emissions were lofted to 6-8 km in a WCB before being entrained in the same cyclone's dry airstream and transported down to the observatory. The results of this study show that the lower marine FT may provide an important transport environment where O3 production may continue, in contrast to transport in the marine boundary layer, where O3 destruction is believed to dominate. The second phase of the study presented here focused on improving the analysis methods that are available with LPDMs. While LPDMs are popular and useful for the analysis of atmospheric trace gas measurements, identifying the transport pathway of emissions from their source to a receptor (the Pico Mountain observatory in our case) using the standard gridded model output, particularly during complex meteorological scenarios can be difficult can be difficult or impossible. The transport study in phase 1 was limited to only 1 month out of more than 3 years of available data and included only 4 case studies out of the 16 events specifically due to this confounding factor. The second phase of this study addressed this difficulty by presenting a method to clearly and easily identify the pathway taken by only those emissions that arrive at a receptor at a particular time, by combining the standard gridded output from forward (i.e., concentrations) and backward (i.e., residence time) LPDM simulations, greatly simplifying similar analyses. The ability of the method to successfully determine the source-to-receptor pathway, restoring this Lagrangian information that is lost when the data are gridded, is proven by comparing the pathway determined from this method with the particle trajectories from both the forward and backward models. A sample analysis is also presented, demonstrating that this method is more accurate and easier to use than existing methods using standard LPDM products. Finally, we discuss potential future work that would be possible by combining the backward LPDM simulation with gridded data from other sources (e.g., chemical transport models) to obtain a Lagrangian sampling of the air that will eventually arrive at a receptor.

Nitrogen removal and sustainability of vertical flow constructed wetlands for small scale wastewater treatment

The challenge for wastewater professionals is to design and operate treatment processes that support human well being and are environmentally sensitive throughout the life-cycle. This research focuses on one technology for small-scale wastewater treatment: the vertical flow constructed wetland (VFCW), which is herein investigated for the capacity to remove ammonium and nitrate nitrogen from wastewater. Hydraulic regime and presence/absence of vegetation are the basis for a three-phase bench scale experiment to determine oxygen transfer and nitrogen fate in VFCWs. Results show that 90% NH4+-N removal is achieved in aerobic downflow columns, 60% NO3--N removal occurs in anaerobic upflow columns, and 60% removal of total nitrogen can be achieved in downflow-upflow in-series. The experimental results are studied further using a variably saturated flow and reactive transport model, which allows a mechanistic explanation of the fate and transport of oxygen and nitrogen. The model clarifies the mechanisms of oxygen transport and nitrogen consumption, and clarifies the need for readily biodegradable COD for denitrification. A VFCW is then compared to a horizontal flow constructed wetland (HFCW) for life cycle environmental impacts. High areal emissions of greenhouse gases from VFCWs compared to HFCWs are the driver for the study. The assessment shows that because a VFCW is only 25% of the volume of an HFCW designed for the same treatment quality, the VFCW has only 25-30% of HFCW impacts over 12 impact categories and 3 damage categories. Results show that impacts could be reduced by design improvements. Design recommendations are downflow wetlands for nitrification, upflow wetlands for denitrification, series wetlands for total nitrogen removal, hydraulic load of 142 L/m2d, 30 cm downflow wetland depth, 1.0 m upflow wetland depth, recycle, vegetation and medium-grained sand. These improvements will optimize nitrogen removal, minimize gaseous emissions, and reduce wetland material requirements, thus reducing environmental impact without sacrificing wastewater treatment quality.

Electron transport in E x B devices

A Hall thruster, an E × B device used for in-space propulsion, utilizes an axial electric field to electrostatically accelerate plasma propellant from the spacecraft. The axial electric field is created by positively biasing the anode so that the positivelycharged ions may be accelerated (repelled) from the thruster, which produces thrust. However, plasma electrons are much smaller than ions and may be accelerated much more quickly toward the anode; if electrons were not impeded, a "short circuit" due to the electron flow would eliminate the thrust mechanism. Therefore, a magnetic field serves to "magnetize" plasma electrons internal to the thruster and confines them in gyro-orbits within the discharge channel. Without outside factors electrons would be confined indefinitely; however, electron-neutral collisions provide a mechanism to free electrons from their orbits allowing electrons to cross the magnetic field toward the anode, where this process is described by classical transport theory. To make matters worse, cross-field electron transport has been observed to be 100-1000 times that predicted by classical collisional theory, providing an efficiency loss mechanism and an obstacle for modeling and simulations in Hall thrusters. The main difficulty in studying electron transport in Hall thrusters is the coupling that exists between the plasma and the fields, where the plasma creates and yet is influenced by the electric field. A device has been constructed at MTU’s Isp Lab, the Hall Electron Mobility Gage, which was designed specifically to study electron transport in E × B devices, where the coupling between the plasma and electric field was virtually eliminated. In this device the two most cited contributors to electron transport in Hall thrusters, fluctuation-induced transport, and wall effects, were absent. Removing the dielectric walls and plasma fluctuations, while maintaining the field environment in vacuum, has allowed the study of electron dynamics in Hall thruster fields where the electrons behave as test particles in prescribed fields, greatly simplifying the environment. Therefore, it was possible to observe any effects on transport not linked to the cited mechanisms, and it was possible to observe trends of the enhanced mobility with control parameters of electric and magnetic fields and neutral density– parameters that are not independently variable in a Hall thruster. The result of the investigation was the observation of electron transport that was ~ 20-100 times the classical prediction. The cross-field electron transport in the Mobility Gage was generally lower than that found in a Hall thruster so these findings do not negate the possibility of fluctuations and/or wall collisions contributing to transport in a Hall thruster. However, this research led to the observation of enhanced cross-field transport that had not been previously isolated in Hall thruster fields, which is not reliant on momentum-transfer collisions, wall collisions or fluctuations.

An investigation of phenolic glycoside and condensed tannin homeostasis in Populus by salicyl alcohol feeding to cell cultures and by transgenic manipulation of the sucrose transporter, PTSUT4, in planta

Secondary metabolites play an important role in plant protection against biotic and abiotic stress. In Populus, phenolic glycosides (PGs) and condensed tannins (CTs) are two such groups of compounds derived from the common phenylpropanoid pathway. The basal levels and the inducibility of PGs and CTs depend on genetic as well as environmental factors, such as soil nitrogen (N) level. Carbohydrate allocation, transport and sink strength also affect PG and CT levels. A negative correlation between the levels of PGs and CTs was observed in several studies. However, the molecular mechanism underlying such relation is not known. We used a cell culture system to understand negative correlation of PGs and CTs. Under normal culture conditions, neither salicin nor higher-order PGs accumulated in cell cultures. Several factors, such as hormones, light, organelles and precursors were discussed in the context of aspen suspension cells’ inability to synthesize PGs. Salicin and its isomer, isosalicin, were detected in cell cultures fed with salicyl alcohol, salicylaldehyde and helicin. At higher levels (5 mM) of salicyl alcohol feeding, accumulation of salicins led to reduced CT production in the cells. Based on metabolic and gene expression data, the CT reduction in salicin-accumulating cells is partly a result of regulatory changes at the transcriptional level affecting carbon partitioning between growth processes, and phenylpropanoid CT biosynthesis. Based on molecular studies, the glycosyltransferases, GT1-2 and GT1-246, may function in glycosylation of simple phenolics, such as salicyl alcohol in cell cultures. The uptake of such glycosides into vacuole may be mediated to some extent by tonoplast localized multidrug-resistance associated protein transporters, PtMRP1 and PtMRP6. In Populus, sucrose is the common transported carbohydrate and its transport is possibly regulated by sucrose transporters (SUTs). SUTs are also capable of transporting simple PGs, such as salicin. Therefore, we characterized the SUT gene family in Populus and investigated, by transgenic analysis, the possible role of the most abundantly expressed member, PtSUT4, in PG-CT homeostasis using plants grown under varying nitrogen regimes. PtSUT4 transgenic plants were phenotypically similar to the wildtype plants except that the leaf area-to-stem volume ratio was higher for transgenic plants. In SUT4 transgenics, levels of non-structural carbohydrates, such as sucrose and starch, were altered in mature leaves. The levels of PGs and CTs were lower in green tissues of transgenic plants under N-replete, but were higher under N-depleted conditions, compared to the levels in wildtype plants. Based on our results, SUT4 partly regulates N-level dependent PG-CT homeostasis by differential carbohydrate allocation.

Theoretical estimation of optical absorption and photoluminescence in nanostructured silicon; an approach to improve efficiency in nanostructured solar cells

Renewable energy is growing in demand, and thus the the manufacture of solar cells and photovoltaic arrays has advanced dramatically in recent years. This is proved by the fact that the photovoltaic production has doubled every 2 years, increasing by an average of 48% each year since 2002. Covering the general overview of solar cell working, and its model, this thesis will start with the three generations of photovoltaic solar cell technology, and move to the motivation of dedicating research to nanostructured solar cell. For the current generation solar cells, among several factors, like photon capture, photon reflection, carrier generation by photons, carrier transport and collection, the efficiency also depends on the absorption of photons. The absorption coefficient,α, and its dependence on the wavelength, λ, is of major concern to improve the efficiency. Nano-silicon structures (quantum wells and quantum dots) have a unique advantage compared to bulk and thin film crystalline silicon that multiple direct and indirect band gaps can be realized by appropriate size control of the quantum wells. This enables multiple wavelength photons of the solar spectrum to be absorbed efficiently. There is limited research on the calculation of absorption coefficient in nano structures of silicon. We present a theoretical approach to calculate the absorption coefficient using quantum mechanical calculations on the interaction of photons with the electrons of the valence band. One model is that the oscillator strength of the direct optical transitions is enhanced by the quantumconfinement effect in Si nanocrystallites. These kinds of quantum wells can be realized in practice in porous silicon. The absorption coefficient shows a peak of 64638.2 cm-1 at = 343 nm at photon energy of ξ = 3.49 eV ( = 355.532 nm). I have shown that a large value of absorption coefficient α comparable to that of bulk silicon is possible in silicon QDs because of carrier confinement. Our results have shown that we can enhance the absorption coefficient by an order of 10, and at the same time a nearly constant absorption coefficient curve over the visible spectrum. The validity of plots is verified by the correlation with experimental photoluminescence plots. A very generic comparison for the efficiency of p-i-n junction solar cell is given for a cell incorporating QDs and sans QDs. The design and fabrication technique is discussed in brief. I have shown that by using QDs in the intrinsic region of a cell, we can improve the efficiency by a factor of 1.865 times. Thus for a solar cell of efficiency of 26% for first generation solar cell, we can improve the efficiency to nearly 48.5% on using QDs.

INTEGRATED MANAGEMENT OF THE INVASIVE SPECIES CIRSIUM ARVENSE, CANADA THISTLE, AND PHRAGMITES AUSTRALIS, COMMON REED: USING ECOLOGICALLY BASED CONTROL MEASURES AT EACH STAGE OF THE INVASION PROCESS

Invasive plant species threaten natural areas by reducing biodiversity and altering ecosystem functions. They also impact agriculture by reducing crop and livestock productivity. Millions of dollars are spent on invasive species control each year, and traditionally, herbicides are used to manage invasive species. Herbicides have human and environmental health risks associated with them; therefore, it is essential that land managers and stakeholders attempt to reduce these risks by utilizing the principles of integrated weed management. Integrated weed management is a practice that incorporates a variety of measures and focuses on the ecology of the invasive plant to manage it. Roadways are high risk areas that have high incidence of invasive species. Roadways act as conduits for invasive species spread and are ideal harborages for population growth; therefore, roadways should be a primary target for invasive species control. There are four stages in the invasion process which an invasive species must overcome: transport, establishment, spread, and impact. The aim of this dissertation was to focus on these four stages and examine the mechanisms underlying the progression from one stage to the next, while also developing integrated weed management strategies. The target species were Phragmites australis, common reed, and Cisrium arvense, Canada thistle. The transport and establishment risks of P. australis can be reduced by removing rhizome fragments from soil when roadside maintenance is performed. The establishment and spread of C. arvense can be reduced by planting particular resistant species, e.g. Heterotheca villosa, especially those that can reduce light transmittance to the soil. Finally, the spread and impact of C. arvense can be mitigated on roadsides through the use of the herbicide aminopyralid. The risks associated with herbicide drift produced by application equipment can be reduced by using the Wet-Blade herbicide application system.

LIFE CYCLE ASSESSMENTS (LCAs) OF PYROLYSIS-BASED GASOLINE AND DIESEL FROM DIFFERENT REGIONAL FEEDSTOCKS: CORN STOVER, SWITCHGRASS, SUGAR CANE BAGASSE, WASTE WOOD, GUINEA GRASS, ALGAE, AND ALBIZIA

Renewable hydrocarbon biofuels are being investigated as possible alternatives to conventional liquid transportation fossil fuels like gasoline, kerosene (aviation fuel), and diesel. A diverse range of biomass feedstocks such as corn stover, sugarcane bagasse, switchgrass, waste wood, and algae, are being evaluated as candidates for pyrolysis and catalytic upgrading to produce drop-in hydrocarbon fuels. This research has developed preliminary life cycle assessments (LCA) for each feedstock-specific pathway and compared the greenhouse gas (GHG) emissions of the hydrocarbon biofuels to current fossil fuels. As a comprehensive study, this analysis attempts to account for all of the GHG emissions associated with each feedstock pathway through the entire life cycle. Emissions from all stages including feedstock production, land use change, pyrolysis, stabilizing the pyrolysis oil for transport and storage, and upgrading the stabilized pyrolysis oil to a hydrocarbon fuel are included. In addition to GHG emissions, the energy requirements and water use have been evaluated over the entire life cycle. The goal of this research is to help understand the relative advantages and disadvantages of the feedstocks and the resultant hydrocarbon biofuels based on three environmental indicators; GHG emissions, energy demand, and water utilization. Results indicate that liquid hydrocarbon biofuels produced through this pyrolysis-based pathway can achieve greenhouse gas emission savings of greater than 50% compared to petroleum fuels, thus potentially qualifying these biofuels under the US EPA RFS2 program. GHG emissions from biofuels ranged from 10.7-74.3 g/MJ from biofuels derived from sugarcane bagasse and wild algae at the extremes of this range, respectively. The cumulative energy demand (CED) shows that energy in every biofuel process is primarily from renewable biomass and the remaining energy demand is mostly from fossil fuels. The CED for biofuel range from 1.25-3.25 MJ/MJ from biofuels derived from sugarcane bagasse to wild algae respectively, while the other feedstock-derived biofuels are around 2 MJ/MJ. Water utilization is primarily from cooling water use during the pyrolysis stage if irrigation is not used during the feedstock production stage. Water use ranges from 1.7 - 17.2 gallons of water per kg of biofuel from sugarcane bagasse to open pond algae, respectively.