Volcanic processes and human exposure as elements to build a risk model for Volcan de Fuego, Guatemala

The activity of Fuego volcano during the 1999 - 2013 eruptive episode is studied through field, remote sensing and observatory records. Mapping of the deposits allows quantifying the erupted volumes and areas affected by the largest eruptions during this period. A wide range of volcanic processes results in a diversity of products and associated deposits, including minor airfall tephra, rockfall avalanches, lava flows, and pyroclastic flows. The activity can be characterized by long term, low level background activity, and sporadic larger explosive eruptions. Although the background activity erupts lava and ash at a low rate (~ 0.1 m3/s), the persistence of such activity over time results in a significant contribution (~ 30%) to the eruption budget during the studied period. Larger eruptions produced the majority of the volume of products during the studied period, mainly during three large events (May 21, 1999, June 29, 2003, and September 13, 2012), mostly in the form of pyroclastic flows. A total volume of ~ 1.4 x 108 m3 was estimated from the mapped deposits and the estimated background eruption rate. Posterior remobilization of pyroclastic flow material by stream erosion in the highly confined Barranca channels leads to lahar generation, either by normal rainfall, or by extreme rainfall events. A reassessment of the types of products and volumes erupted during the decade of 1970's allows comparing the activity happening since 1999 with the older activity, and suggests that many of the eruptive phenomena at Fuego may have similar mechanisms, despite the differences in scale between. The deposits of large pyroclastic flows erupted during the 1970's are remarkably similar in appearance to the deposit of pyroclastic flows from the 1999 - 2013 period, despite their much larger volume; this is also the case for prehistoric eruptions. Radiocarbon dating of pyroclastic flow deposits suggests that Fuego has produced large eruptions many times during the last ~ 2 ka, including larger eruptions during the last 500 years, which has important hazard implications. A survey was conducted among the local residents living near to the volcano, about their expectations of possible future crises. The results show that people are aware of the risk they could face in case of a large eruption, and therefore they are willing to evacuate in such case. However, their decision to evacuate may also be influenced by the conditions in which the evacuation could take place. If the evacuation represents a potential loss of their livelihood or property they will be more hesitant to leave their villages during a large eruption. The prospect of facing hardship conditions during the evacuation and in the shelters may further cause reluctance to evacuate. A short discussion on some of the issues regarding risk assessment and management through an early warning system is presented in the last chapter.

Energy harvesting from body motion using rotational micro-generation

Autonomous system applications are typically limited by the power supply operational lifetime when battery replacement is difficult or costly. A trade-off between battery size and battery life is usually calculated to determine the device capability and lifespan. As a result, energy harvesting research has gained importance as society searches for alternative energy sources for power generation. For instance, energy harvesting has been a proven alternative for powering solar-based calculators and self-winding wristwatches. Thus, the use of energy harvesting technology can make it possible to assist or replace batteries for portable, wearable, or surgically-implantable autonomous systems. Applications such as cardiac pacemakers or electrical stimulation applications can benefit from this approach since the number of surgeries for battery replacement can be reduced or eliminated. Research on energy scavenging from body motion has been investigated to evaluate the feasibility of powering wearable or implantable systems. Energy from walking has been previously extracted using generators placed on shoes, backpacks, and knee braces while producing power levels ranging from milliwatts to watts. The research presented in this paper examines the available power from walking and running at several body locations. The ankle, knee, hip, chest, wrist, elbow, upper arm, side of the head, and back of the head were the chosen target localizations. Joints were preferred since they experience the most drastic acceleration changes. For this, a motor-driven treadmill test was performed on 11 healthy individuals at several walking (1-4 mph) and running (2-5 mph) speeds. The treadmill test provided the acceleration magnitudes from the listed body locations. Power can be estimated from the treadmill evaluation since it is proportional to the acceleration and frequency of occurrence. Available power output from walking was determined to be greater than 1mW/cm³ for most body locations while being over 10mW/cm³ at the foot and ankle locations. Available power from running was found to be almost 10 times higher than that from walking. Most energy harvester topologies use linear generator approaches that are well suited to fixed-frequency vibrations with sub-millimeter amplitude oscillations. In contrast, body motion is characterized with a wide frequency spectrum and larger amplitudes. A generator prototype based on self-winding wristwatches is deemed to be appropriate for harvesting body motion since it is not limited to operate at fixed-frequencies or restricted displacements. Electromagnetic generation is typically favored because of its slightly higher power output per unit volume. Then, a nonharmonic oscillating rotational energy scavenger prototype is proposed to harness body motion. The electromagnetic generator follows the approach from small wind turbine designs that overcome the lack of a gearbox by using a larger number of coil and magnets arrangements. The device presented here is composed of a rotor with multiple-pole permanent magnets having an eccentric weight and a stator composed of stacked planar coils. The rotor oscillations induce a voltage on the planar coil due to the eccentric mass unbalance produced by body motion. A meso-scale prototype device was then built and evaluated for energy generation. The meso-scale casing and rotor were constructed on PMMA with the help of a CNC mill machine. Commercially available discrete magnets were encased in a 25mm rotor. Commercial copper-coated polyimide film was employed to manufacture the planar coils using MEMS fabrication processes. Jewel bearings were used to finalize the arrangement. The prototypes were also tested at the listed body locations. A meso-scale generator with a 2-layer coil was capable to extract up to 234 µW of power at the ankle while walking at 3mph with a 2cm³ prototype for a power density of 117 µW/cm³. This dissertation presents the analysis of available power from walking and running at different speeds and the development of an unobtrusive miniature energy harvesting generator for body motion. Power generation indicates the possibility of powering devices by extracting energy from body motion.

IDENTIFICATION OF GENES CONTROLLING BIOLOGICAL PROCESSES AND PATHWAYS THROUGH STATISTICAL ANALYSIS AND NETWORK RECONSTRUCTION

The developmental processes and functions of an organism are controlled by the genes and the proteins that are derived from these genes. The identification of key genes and the reconstruction of gene networks can provide a model to help us understand the regulatory mechanisms for the initiation and progression of biological processes or functional abnormalities (e.g. diseases) in living organisms. In this dissertation, I have developed statistical methods to identify the genes and transcription factors (TFs) involved in biological processes, constructed their regulatory networks, and also evaluated some existing association methods to find robust methods for coexpression analyses. Two kinds of data sets were used for this work: genotype data and gene expression microarray data. On the basis of these data sets, this dissertation has two major parts, together forming six chapters. The first part deals with developing association methods for rare variants using genotype data (chapter 4 and 5). The second part deals with developing and/or evaluating statistical methods to identify genes and TFs involved in biological processes, and construction of their regulatory networks using gene expression data (chapter 2, 3, and 6). For the first part, I have developed two methods to find the groupwise association of rare variants with given diseases or traits. The first method is based on kernel machine learning and can be applied to both quantitative as well as qualitative traits. Simulation results showed that the proposed method has improved power over the existing weighted sum method (WS) in most settings. The second method uses multiple phenotypes to select a few top significant genes. It then finds the association of each gene with each phenotype while controlling the population stratification by adjusting the data for ancestry using principal components. This method was applied to GAW 17 data and was able to find several disease risk genes. For the second part, I have worked on three problems. First problem involved evaluation of eight gene association methods. A very comprehensive comparison of these methods with further analysis clearly demonstrates the distinct and common performance of these eight gene association methods. For the second problem, an algorithm named the bottom-up graphical Gaussian model was developed to identify the TFs that regulate pathway genes and reconstruct their hierarchical regulatory networks. This algorithm has produced very significant results and it is the first report to produce such hierarchical networks for these pathways. The third problem dealt with developing another algorithm called the top-down graphical Gaussian model that identifies the network governed by a specific TF. The network produced by the algorithm is proven to be of very high accuracy.

Characterizing the effects of turbulence on sediment oxygen and sediment nitrate demand using flow-through reactors

Onondaga Lake has received the municipal effluent and industrial waste from the city of Syracuse for more than a century. Historically, 75 metric tons of mercury were discharged to the lake by chlor-alkali facilities. These legacy deposits of mercury now exist primarily in the lake sediments. Under anoxic conditions, methylmercury is produced in the sediments and can be released to the overlying water. Natural sedimentation processes are continuously burying the mercury deeper into the sediments. Eventually, the mercury will be buried to a depth where it no longer has an impact on the overlying water. In the interim, electron acceptor amendment systems can be installed to retard these chemical releases while the lake naturally recovers. Electron acceptor amendment systems are designed to meet the sediment oxygen demand in the sediment and maintain manageable hypolimnion oxygen concentrations. Historically, designs of these systems have been under designed resulting in failure. This stems from a mischaracterization of the sediment oxygen demand. Turbulence at the sediment water interface has been shown to impact sediment oxygen demand. The turbulence introduced by the electron amendment system can thus increase the sediment oxygen demand, resulting in system failure if turbulence is not factored into the design. Sediment cores were gathered and operated to steady state under several well characterized turbulence conditions. The relationship between sediment oxygen/nitrate demand and turbulence was then quantified and plotted. A maximum demand was exhibited at or above a fluid velocity of 2.0 mm•s-1. Below this velocity, demand decreased rapidly with fluid velocity as zero velocity was approached. Similar relationships were displayed by both oxygen and nitrate cores.

Dynamics of settling charged particles in turbulence : theory and experiments

It has been proposed that inertial clustering may lead to an increased collision rate of water droplets in clouds. Atmospheric clouds and electrosprays contain electrically charged particles embedded in turbulent flows, often under the influence of an externally imposed, approximately uniform gravitational or electric force. In this thesis, we present the investigation of charged inertial particles embedded in turbulence. We have developed a theoretical description for the dynamics of such systems of charged, sedimenting particles in turbulence, allowing radial distribution functions to be predicted for both monodisperse and bidisperse particle size distributions. The governing parameters are the particle Stokes number (particle inertial time scale relative to turbulence dissipation time scale), the Coulomb-turbulence parameter (ratio of Coulomb ’terminalar speed to turbulence dissipation velocity scale), and the settling parameter (the ratio of the gravitational terminal speed to turbulence dissipation velocity scale). For the monodispersion particles, The peak in the radial distribution function is well predicted by the balance between the particle terminal velocity under Coulomb repulsion and a time-averaged ’drift’ velocity obtained from the nonuniform sampling of fluid strain and rotation due to finite particle inertia. The theory is compared to measured radial distribution functions for water particles in homogeneous, isotropic air turbulence. The radial distribution functions are obtained from particle positions measured in three dimensions using digital holography. The measurements support the general theoretical expression, consisting of a power law increase in particle clustering due to particle response to dissipative turbulent eddies, modulated by an exponential electrostatic interaction term. Both terms are modified as a result of the gravitational diffusion-like term, and the role of ’gravity’ is explored by imposing a macroscopic uniform electric field to create an enhanced, effective gravity. The relation between the radial distribution functions and inward mean radial relative velocity is established for charged particles.

DIFFUSE-INTERFACE FIELD APPROACH TO MODELING SELF-ASSEMBLY OF HETEROGENEOUS COLLOIDAL SYSTEMS AND RELATED DIPOLE-DIPOLE INTERACTION PHENOMENA

Colloid self-assembly under external control is a new route to fabrication of advanced materials with novel microstructures and appealing functionalities. The kinetic processes of colloidal self-assembly have attracted great interests also because they are similar to many atomic level kinetic processes of materials. In the past decades, rapid technological progresses have been achieved on producing shape-anisotropic, patchy, core-shell structured particles and particles with electric/magnetic charges/dipoles, which greatly enriched the self-assembled structures. Multi-phase carrier liquids offer new route to controlling colloidal self-assembly. Therefore, heterogeneity is the essential characteristics of colloid system, while so far there still lacks a model that is able to efficiently incorporate these possible heterogeneities. This thesis is mainly devoted to development of a model and computational study on the complex colloid system through a diffuse-interface field approach (DIFA), recently developed by Wang et al. This meso-scale model is able to describe arbitrary particle shape and arbitrary charge/dipole distribution on the surface or body of particles. Within the framework of DIFA, a Gibbs-Duhem-type formula is introduced to treat Laplace pressure in multi-liquid-phase colloidal system and it obeys Young-Laplace equation. The model is thus capable to quantitatively study important capillarity related phenomena. Extensive computer simulations are performed to study the fundamental behavior of heterogeneous colloidal system. The role of Laplace pressure is revealed in determining the mechanical equilibrium of shape-anisotropic particles at fluid interfaces. In particular, it is found that the Laplace pressure plays a critical role in maintaining the stability of capillary bridges between close particles, which sheds light on a novel route to in situ firming compact but fragile colloidal microstructures via capillary bridges. Simulation results also show that competition between like-charge repulsion, dipole-dipole interaction and Brownian motion dictates the degree of aggregation of heterogeneously charged particles. Assembly and alignment of particles with magnetic dipoles under external field is studied. Finally, extended studies on the role of dipole-dipole interaction are performed for ferromagnetic and ferroelectric domain phenomena. The results reveal that the internal field generated by dipoles competes with external field to determine the dipole-domain evolution in ferroic materials.

Spatial and temporal patterns in ungulate-ecosystem interactions

Ungulates are important components of a variety of ecosystems worldwide. This dissertation integrates aspects of ungulate and forest ecology to increase our understanding of how they work together in ways that are of interest to natural resource managers, educators, and those who are simply curious about nature. Although animal ecology and ecosystem ecology are often studied separately, one of the general goals of this dissertation is to examine how they interact across spatial and temporal scales. Forest ecosystems are heterogeneous across a range of scales. Spatial and temporal habitat use patterns of forest ungulates tend to be congregated in patches where food and/or cover are readily available. Ungulates interact with ecosystem processes by selectively foraging on plants and excreting waste products in concentrated patches. Positive feedbacks may develop where these activities increase the value of habitat through soil fertilization or the alteration of plant chemistry and architecture. Heterogeneity in ecosystem processes and plant community structure, observed at both stand and local scales, may be the integrated outcome of feedbacks between ungulate behavior and abiotic resource gradients. The first chapter of this dissertation briefly discusses pertinent background information on ungulate ecology, with a focus on white-tailed deer (Odocoileus virginianus) in the Upper Great Lakes region and moose (Alces acles) in Isle Royale National Park, Michigan, USA. The second chapter demonstrates why ecological context is important for studying ungulate ecology in forest ecosystems. Excluding deer from eastern hemlock (Tsuga canadensis) stands, which deer use primarily as winter cover, resulted in less spatial complexity in soil reactive nitrogen and greater complexity in diffuse light compared to unfenced stands. The spatial patterning of herbaceous-layer cover was more similar to nitrogen where deer were present, and was a combination of nitrogen and light within deer exclosures. This relationship depends on the seasonal timing of deer habitat use because deer fertilize the soil during winter, but leave during the growing season. The third chapter draws upon an eight-year, 39-stand data set of deer fecal pellet counts in hemlock stands to estimate the amount of nitrogen that deer are depositing in hemlock stands each winter. In stands of high winter deer use, deer-excreted nitrogen inputs consistently exceeded those of atmospheric deposition at the stand scale. At the neighborhood scale, deer-excreted nitrogen was often in excess of atmospheric deposition due to the patchy distribution of deer habitat use. Spatial patterns in habitat use were consistent over the eight-year study at both stand and neighborhood scales. The fourth chapter explores how foraging selectivity by moose interacts with an abiotic resource gradient to influence forest structure and composition. Soil depth on Isle Royale varies from east to west according to glacial history. Fir saplings growing in deeper soils on the west side are generally more palatable forage for moose (lower foliar C:N) than those growing in shallower soils on the east side. Therefore, saplings growing in better conditions are less likely to reach the canopy due to moose browsing, and fir is a smaller overstory component on the west side. Lastly, chapter five focuses on issues surrounding eastern hemlock regeneration failure, which is a habitat type that is important to many wildlife species. Increasing hemlock on the landscape is complicated by several factors including disturbance regime and climate change, in addition to the influence of deer.

Stilling Basin Scour Remediation Using Air Injection and Flat Plate Extension

The South Florida Water Management District (SFWMD) is responsible for managing over 2500 miles of waterways and hundreds of water control structures. Many of these control structures are experiencing erosion, known as scour, of the sediment downstream of the structure. Laboratory experiments were conducted in order to investigate the effectiveness of two-dimensional air diffusers and plate extensions (without air injection) on a 1/30 scale model of one of SFWMD gated spillway structures, the S65E gated spillway. A literature review examining the results of similar studies was conducted. The experimental design for this research was based off of previous work done on the same model. Scour of the riverbed downstream of gated spillway structures has the potential to cause serious damage, as it can expose the foundation of the structure, which can lead to collapse. This type of scour has been studied previously, but it continues to pose a risk to water control structures and needs to be studied further. The hydraulic scour channel used to conduct experiments contains a head tank, flow straighteners, gated spillway, stilling basin, scour chamber, sediment trap, and tailwater tank. Experiments were performed with two types of air diffusers. The first was a hollow, acrylic, triangular end sill with air injection holes on the upstream face, allowing for air injection upstream. The second diffuser was a hollow, acrylic rectangle that extended from the triangular end sill with air injection holes in the top face, allowing for vertical air injection, perpendicular to flow. Detailed flow and bed measurements were taken for six trials for each diffuser ranging from no air injection to 5 rows of 70 holes of 0.04" diameter. It was found that with both diffusers, the maximum amount of air injection reduced scour the most. Detailed velocity measurements were taken for each case and turbulence statistics were analyzed to determine why air injection reduces scour. It was determined that air injection reduces streamwise velocity and turbulence. Another set of experiments was performed using an acrylic extension plate with no air injection to minimize energy costs. Ten different plate lengths were tested. It was found that the location of deepest scour moved further downstream with each plate length. The 32-cm plate is recommended here. Detailed velocity measurements were taken after the cases with the 32-cm plate and no plate had reached equilibrium. This was done to better understand the flow patterns in order to determine what causes the scour reduction with the extension plates. The extension plate reduces the volume of scour, but more importantly translates the deepest point of scour downstream from the structure, lessening the risk of damage.