Biogeochemical processes in a clay formation in situ experiment: Part F - Reactive transport modelling

Reactive transport modelling was used to simulate solute transport, thermodynamic reactions, ion exchange and biodegradation in the Porewater Chemistry (PC) experiment at the Mont Terri Rock Laboratory. Simulations show that the most important chemical processes controlling the fluid composition within the borehole and the surrounding formation during the experiment are ion exchange, biodegradation and dissolution/precipitation reactions involving pyrite and carbonate minerals. In contrast, thermodynamic mineral dissolution/precipitation reactions involving alumo-silicate minerals have little impact on the fluid composition on the time-scale of the experiment. With the accurate description of the initial chemical condition in the formation in combination with kinetic formulations describing the different stages of bacterial activities, it has been possible to reproduce the evolution of important system parameters, such as the pH, redox potential, total organic C. dissolved inorganic C and SO(4) concentration. Leaching of glycerol from the pH-electrode may be the primary source of organic material that initiated bacterial growth, which caused the chemical perturbation in the borehole. Results from these simulations are consistent with data from the over-coring and demonstrate that the Opalinus Clay has a high buffering capacity in terms of chemical perturbations caused by bacterial activity. This buffering capacity can be attributed to the carbonate system as well as to the reactivity of clay surfaces.

Fine structure of hippocampal mossy fiber synapses following rapid high-pressure freezing

Synapses of hippocampal neurons play important roles in learning and memory processes and are involved in aberrant hippocampal function in temporal lobe epilepsy. Major neuronal types in the hippocampus as well as their input and output synapses are well known, but it has remained an open question to what extent conventional electron microscopy (EM) has provided us with the real appearance of synaptic fine structure under in vivo conditions. There is reason to assume that conventional aldehyde fixation and dehydration lead to protein denaturation and tissue shrinkage, likely associated with the occurrence of artifacts. However, realistic fine-structural data of synapses are required for our understanding of the transmission process and for its simulation. Here, we used high-pressure freezing and cryosubstitution of hippocampal tissue that was not subjected to aldehyde fixation and dehydration in ethanol to monitor the fine structure of an identified synapse in the hippocampal CA3 region, that is, the synapse between granule cell axons, the mossy fibers, and the proximal dendrites of CA3 pyramidal neurons. Our results showed that high-pressure freezing nicely preserved ultrastructural detail of this particular synapse and allowed us to study rapid structural changes associated with synaptic plasticity.

Sudden Increase and Regrowth of Fecal Coliforms and Escherichia Coli in Wastewater Biosolids After High Solids Centrifugation

Treatment plants that operate either thermophilic or mesophilic anaerobic digesters with centrifugal dewatering processes have consistently observed densities of fecal coliform and Escherichia coli, both indicator bacteria, that decrease during digestion but then increase after dewatering and storage. The increases have been characterized as two separate phenomena to explain this observation: 1) “Sudden Increase,” or SI, which is defined as the increase that occurs immediately after dewatering and 2) “regrowth,” which is defined as an increase during storage of cake samples over a period of hours or days. The SI observation appears to be more prevalent with biosolids that are generated with thermophilic processes and dewatered by centrifugation. Both thermophilic and mesophilic digesters with centrifuge dewatering processes have observed the regrowth phenomena. This research hypothesizes that the SI phenomenon is due to the presence of viable nonculturable (VNC) bacteria that are reactivated during dewatering. In other words, the bacteria were always present but were not enumerated by standard culturing methods (SCM). Analysis of the E. coli density in thermally treated solids by SCMs and quantitative real-time polymerase chain reaction (qPCR) indicated that E. coli densities are often underestimated by SCM. When analyzed with qPCR, the E. coli density after digestion can be 4-5 orders of magnitude greater than the non-detect levels identified by SCMs, which supports the non-culturable hypothesis. The VNC state describes a condition where bacteria are alive but unable to sustain the metabolic process needed for cellular division. Supplements added to culturing media were investigated to determine if the resuscitation of VNC bacteria could be enhanced. The autoinducer molecules Nhexanoyl- L-Homoserine lactone (C6-HSL), 3-oxo-N-octanoyl-L-Homoserine lactone (3-oxo- C8-HSL), and norepinephrine were unable to induce the resuscitation of VNC E. coli. Additional sampling was performed to determine if autoinducer molecules, peroxides, or other as of yet unknown inhibitory agents and toxins could be removed from biosolids during SCM. Culture media supplemented with the peroxide degrading compounds catalase, α-ketoglutaric acid, and sodium pyruvate was unable to resuscitate non-culturable E. coli. The additions of bentonite and exponential growth phase E. coli cell-free supernatant to culturing media were also unable to increase the culturability of E. coli. To remove inhibitory agents and toxins, a cell washing technique was employed prior to performing SCM; however, this cell washing technique may have increased cellular stresses that inhibited resuscitation since cell densities decreased. A novel laboratory-scale dewatering process was also investigated to determine if the SI and regrowth phenomena observed in full-scale centrifugal dewatering could be mimicked in the laboratory using a lab shearing device. Fecal coliform and E. coli densities in laboratory prepared cake samples were observed to be an order of magnitude higher than full-scale dewatered cakes. Additionally, the laboratory-scale dewatering process was able to resuscitate fecal coliforms and E. coli in stored sludge such that the density increased by 4-5 orders of magnitude from nondetect values. Lastly, the addition of aluminum sulfate during centrifuge dewatering at a full-scale utility produced an increased regrowth of fecal coliforms and E. coli that was sustained for 5 days.

Monitoring and Modeling of the Hydraulic and Sediment Processes at In-Stream Restoration Structures in the Vicinity of a Bridge Crossing

The long-term performance of infrastructure depends on reliable and sustainable designs. Many of Pennsylvania’s streams experience sediment transport problems that increase maintenance costs and lower structural integrity of bridge crossings. A stream restoration project is one common mitigation measure used to correct such problems at bridge crossings. Specifically, in an attempt to alleviate aggradation problems with the Old Route 15 Bridge crossing on White Deer Creek, in White Deer, PA, two in-stream structures (rock cross vanes) and several bank stabilization features were installed along with a complete channel redevelopment. The objectives of this research were to characterize the hydraulic and sediment transport processes occurring at the White Deer Creek site, and to investigate, through physical and mathematical modeling, the use of instream restoration structures. The goal is to be able to use the results of this study to prevent aggradation or other sediment related problems in the vicinity of bridges through improved design considerations. Monitoring and modeling indicate that the study site on White Deer Creek is currently unstable, experiencing general channel down-cutting, bank erosion, and several local areas of increased aggradation and degradation of the channel bed. An in-stream structure installed upstream of the Old Route 15 Bridge failed by sediment burial caused by the high sediment load that White Deer Creek is transporting as well as the backwater effects caused by the bridge crossing. The in-stream structure installed downstream of the Old Route 15 Bridge is beginning to fail because of the alignment of the structure with the approach direction of flow from upstream of the restoration structure.

Paleoenvironment and Paleoecology of a Late Paleocene High-Latitude Terrestrial Succession, Arkose Ridge Formation at Box Canyon, Southern Talkeetna Mountains, Alaska

Paleogene sedimentary rocks of the Arkose Ridge Formation (Talkeetna Mountains, Alaska) preserve a record of a fluvial-lacustrine depositional environment and its forested ecosystem in an active basin among the convergent margin tectonic processes that shaped southern Alaska. An -800 m measured succession at Box Canyon indicates braid-plain deposition with predominantly gravelly deposits low in the exposure to sandy and muddy facies associations below an overlying lava flow sequence. U-Pb geochronology on zircons from a tuff and a sandstone within the measured section, as well as an Ar/Ar date from the overlying lava constrain the age of the sedimentary succession to between similar to 59 Ma and 48 Ma Fossil plant remains occur throughout the Arkose Ridge Formation as poorly-preserved coalified woody debris and fragmentary leaf impressions. At Box Canyon, however, a thin la-custrine depositional lens of rhythmically laminated mudrocks yielded fish fossils and a well-preserved floral assemblage including foliage and reproductive organs representing conifers, sphenopsids, monocots, and dicots. Leaf physiognomic methods to estimate paleoclimate were applied to the dicot leaf collection and indicate warm temperate paleotemperatures (-11-15 +/- -4 degrees C MAT) and elevated paleoprecipitation (-120 cm/yr MAP) estimates as compared to modem conditions; results that are parallel with previously published estimates from the partly coeval Chickaloon Formation deposited in more distal depositional environments in the same basin. The low abundance of leaf herbivory in the Box Canyon dicot assemblage (-9% of leaves damaged) is also similar to the results from assemblages in the meander-plain depositional systems of the Chickaloon. This new suite of data informs models of the tectonostratigraphic evolution of southern Alaska and the developing understanding of terrestrial paleoecology and paleoclimate at high latitudes during the Late Paleocene-Early Eocene greenhouse climate phase. (c) 2014 Elsevier B.V. All rights reserved.

Prion removal capacity of plasma protein manufacturing processes

BACKGROUND The variant Creutzfeldt-Jakob disease incidence peaked a decade ago and has since declined. Based on epidemiologic evidence, the causative agent, pathogenic prion, has not constituted a tangible contamination threat to large-scale manufacturing of human plasma-derived proteins. Nonetheless, manufacturers have studied the prion removal capabilities of various manufacturing steps to better understand product safety. Collectively analyzing the results could reveal experimental reproducibility and detect trends and mechanisms driving prion removal. STUDY DESIGN AND METHODS Plasma Protein Therapeutics Association member companies collected more than 200 prion removal studies on plasma protein manufacturing steps, including precipitation, adsorption, chromatography, and filtration, as well as combined steps. The studies used a range of model spiking agents and bench-scale process replicas. The results were grouped based on key manufacturing variables to identify factors impacting removal. The log reduction values of a group are presented for comparison. RESULTS Overall prion removal capacities evaluated by independent groups were in good agreement. The removal capacity evaluated using biochemical assays was consistent with prion infectivity removal measured by animal bioassays. Similar reduction values were observed for a given step using various spiking agents, except highly purified prion protein in some circumstances. Comparison between combined and single-step studies revealed complementary or overlapping removal mechanisms. Steps with high removal capacities represent the conditions where the physiochemical differences between prions and therapeutic proteins are most significant. CONCLUSION The results support the intrinsic ability of certain plasma protein manufacturing steps to remove prions in case of an unlikely contamination, providing a safeguard to products.

A HIDDEN MARKOV MODEL FOR JOINT ESTIMATION OF GENOTYPE AND COPY NUMBER IN HIGH-THROUGHPUT SNP CHIPS

Amplifications and deletions of chromosomal DNA, as well as copy-neutral loss of heterozygosity have been associated with diseases processes. High-throughput single nucleotide polymorphism (SNP) arrays are useful for making genome-wide estimates of copy number and genotype calls. Because neighboring SNPs in high throughput SNP arrays are likely to have dependent copy number and genotype due to the underlying haplotype structure and linkage disequilibrium, hidden Markov models (HMM) may be useful for improving genotype calls and copy number estimates that do not incorporate information from nearby SNPs. We improve previous approaches that utilize a HMM framework for inference in high throughput SNP arrays by integrating copy number, genotype calls, and the corresponding confidence scores when available. Using simulated data, we demonstrate how confidence scores control smoothing in a probabilistic framework. Software for fitting HMMs to SNP array data is available in the R package ICE.

Thrombospondin-1 mediates platelet adhesion at high shear via glycoprotein Ib (GPIb): an alternative/backup mechanism to von Willebrand factor

Acute thrombotic arterial occlusion is the leading cause of morbidity and mortality in the Western world. Von Willebrand factor is thought to be the only indispensable adhesive substrate to promote thrombus formation in high shear environments. We found that thrombospondin-1, a glycoprotein enriched in arteriosclerotic plaques, might function as an alternative substrate for thrombus formation. Platelets adhered to thrombospondin-1 in a shear dependent manner with an optimum shear as found in stenosed arteries. Adhesion is extremely firm, with no detachment of platelets up to a shear rate of 4000 s(-1). Experiments using platelets from a patient completely lacking von Willebrand factor showed that von Willebrand factor is not involved in platelet binding to thrombospondin-1. Platelet adhesion to thrombospondin-1 is not mediated via beta3-integrins or GPIa. CD36 partially mediates the adhesion of pre-activated platelets. We identified GPIb as high shear adhesion-receptor for thrombospondin-1. Soluble GPIb, as well as antibodies against the GPIb, blocked platelet adhesion almost completely. The new discovered thrombospondin-1-GPIb adhesion axis under arterial shear conditions might be important, not only during thrombus formation but also for pathological processes where other cells bind to the endothelium or subendothelium, including arteriosclerosis, inflammation and tumor metastasis, and a promising therapeutic target.

Effect of mesh distortion on the accuracy of high order vorticity-velocity CFD approaches

The numerical solution of the incompressible Navier-Stokes equations offers an alternative to experimental analysis of fluid-structure interaction (FSI). We would save a lot of time and effort and help cut back on costs, if we are able to accurately model systems by these numerical solutions. These advantages are even more obvious when considering huge structures like bridges, high rise buildings or even wind turbine blades with diameters as large as 200 meters. The modeling of such processes, however, involves complex multiphysics problems along with complex geometries. This thesis focuses on a novel vorticity-velocity formulation called the Kinematic Laplacian Equation (KLE) to solve the incompressible Navier-stokes equations for such FSI problems. This scheme allows for the implementation of robust adaptive ordinary differential equations (ODE) time integration schemes, allowing us to tackle each problem as a separate module. The current algortihm for the KLE uses an unstructured quadrilateral mesh, formed by dividing each triangle of an unstructured triangular mesh into three quadrilaterals for spatial discretization. This research deals with determining a suitable measure of mesh quality based on the physics of the problems being tackled. This is followed by exploring methods to improve the quality of quadrilateral elements obtained from the triangles and thereby improving the overall mesh quality. A series of numerical experiments were designed and conducted for this purpose and the results obtained were tested on different geometries with varying degrees of mesh density.

Effect of high order interpolation in the stability and efficiency of the time-integration process in vorticity-velocity CFD algorithms

The numerical solution of the incompressible Navier-Stokes Equations offers an effective alternative to the experimental analysis of Fluid-Structure interaction i.e. dynamical coupling between a fluid and a solid which otherwise is very complex, time consuming and very expensive. To have a method which can accurately model these types of mechanical systems by numerical solutions becomes a great option, since these advantages are even more obvious when considering huge structures like bridges, high rise buildings, or even wind turbine blades with diameters as large as 200 meters. The modeling of such processes, however, involves complex multiphysics problems along with complex geometries. This thesis focuses on a novel vorticity-velocity formulation called the KLE to solve the incompressible Navier-stokes equations for such FSI problems. This scheme allows for the implementation of robust adaptive ODE time integration schemes and thus allows us to tackle the various multiphysics problems as separate modules. The current algorithm for KLE employs a structured or unstructured mesh for spatial discretization and it allows the use of a self-adaptive or fixed time step ODE solver while dealing with unsteady problems. This research deals with the analysis of the effects of the Courant-Friedrichs-Lewy (CFL) condition for KLE when applied to unsteady Stoke’s problem. The objective is to conduct a numerical analysis for stability and, hence, for convergence. Our results confirmthat the time step ∆t is constrained by the CFL-like condition ∆t ≤ const. hα, where h denotes the variable that represents spatial discretization.

Ultraviolet digital imaging of volcanic plumes : implementation and application to magmatic processes at basaltic volcanoes

Magmatic volatiles play a crucial role in volcanism, from magma production at depth to generation of seismic phenomena to control of eruption style. Accordingly, many models of volcano dynamics rely heavily on behavior of such volatiles. Yet measurements of emission rates of volcanic gases have historically been limited, which has restricted model verification to processes on the order of days or longer. UV cameras are a recent advancement in the field of remote sensing of volcanic SO2 emissions. They offer enhanced temporal and spatial resolution over previous measurement techniques, but need development before they can be widely adopted and achieve the promise of integration with other geophysical datasets. Large datasets require a means by which to quickly and efficiently use imagery to calculate emission rates. We present a suite of programs designed to semi-automatically determine emission rates of SO2 from series of UV images. Extraction of high temporal resolution SO2 emission rates via this software facilitates comparison of gas data to geophysical data for the purposes of evaluating models of volcanic activity and has already proven useful at several volcanoes. Integrated UV camera and seismic measurements recorded in January 2009 at Fuego volcano, Guatemala, provide new insight into the system’s shallow conduit processes. High temporal resolution SO2 data reveal patterns of SO2 emission rate relative to explosions and seismic tremor that indicate tremor and degassing share a common source process. Progressive decreases in emission rate appear to represent inhibition of gas loss from magma as a result of rheological stiffening in the upper conduit. Measurements of emission rate from two closely-spaced vents, made possible by the high spatial resolution of the camera, help constrain this model. UV camera measurements at Kilauea volcano, Hawaii, in May of 2010 captured two occurrences of lava filling and draining within the summit vent. Accompanying high lava stands were diminished SO2 emission rates, decreased seismic and infrasonic tremor, minor deflation, and slowed lava lake surface velocity. Incorporation of UV camera data into the multi-parameter dataset gives credence to the likelihood of shallow gas accumulation as the cause of such events.

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.

Understanding volcanic processes using UV camera measurements of sulfur dioxide and coincident infrasound and seismicity

The exsolution of volatiles from magma maintains an important control on volcanic eruption styles. The nucleation, growth, and connectivity of bubbles during magma ascent provide the driving force behind eruptions, and the rate, volume, and ease of gas exsolution can affect eruptive activity. Volcanic plumes are the observable consequence of this magmatic degassing, and remote sensing techniques allow us to quantify changes in gas exsolution. However, until recently the methods used to measure volcanic plumes did not have the capability of detecting rapid changes in degassing on the scale of standard geophysical observations. The advent of the UV camera now makes high sample rate gas measurements possible. This type of dataset can then be compared to other volcanic observations to provide an in depth picture of degassing mechanisms in the shallow conduit. The goals of this research are to develop a robust methodology for UV camera field measurements of volcanic plumes, and utilize this data in conjunction with seismoacoustic records to illuminate degassing processes. Field and laboratory experiments were conducted to determine the effects of imaging conditions, vignetting, exposure time, calibration technique, and filter usage on the UV camera sulfur dioxide measurements. Using the best practices determined from these studies, a field campaign was undertaken at Volcán de Pacaya, Guatemala. Coincident plume sulfur dioxide measurements, acoustic recordings, and seismic observations were collected and analyzed jointly. The results provide insight into the small explosive features, variations in degassing rate, and plumbing system of this complex volcanic system. This research provides useful information for determining volcanic hazard at Pacaya, and demonstrates the potential of the UV camera in multiparameter studies.

High-resolution computer simulation of electrophoretic mobilization in isoelectric focusing

Cationic and anionic electrophoretic mobilization for focusing of hemoglobins (Hb's) in the presence of 100 carrier ampholytes covering a pI range of 6.00-7.98 was studied by computer simulation at a constant current density of 300 A/m(2). Electropherograms that would be produced by whole column imaging and by single detectors placed at different locations along the focusing column are presented. Upon mobilization, peak heights of the Hb zones decrease, but the zones retain a relatively sharp constant profile and are migrating at a constant velocity. A further peak decrease occurs during readjustment at the locations of the original buffer/column interfaces, indicating that detection sensitivity is the lowest at these locations. An anionic carrier ampholyte mobility smaller than that of its cationic species produces a cathodic drift which is smaller than the transport rate used for electrophoretic mobilization. Compared to the case with equal mobilities of carrier ampholyte species, a small increase (decrease) is predicted for the cationic (anionic) mobilization rate within the focusing column. Simulation data suggest that electrophoretic mobilization after focusing and focusing with concurrent electrophoretic mobilization are comparable isotachophoretic processes that occur when there is an uninterrupted flux of an ion through the focusing column. Cathodic drift caused by unequal mobilities of the species of carrier ampholytes, electrophoretic mobilization, and decomposition occurring at the pH gradient edges are related electrophoretic processes.