Predicting organic acid concentration from UV/vis spectrometry measurements - A comparison of machine learning techniques

The concentration of organic acids in anaerobic digesters is one of the most critical parameters for monitoring and advanced control of anaerobic digestion processes. Thus, a reliable online-measurement system is absolutely necessary. A novel approach to obtaining these measurements indirectly and online using UV/vis spectroscopic probes, in conjunction with powerful pattern recognition methods, is presented in this paper. An UV/vis spectroscopic probe from S::CAN is used in combination with a custom-built dilution system to monitor the absorption of fully fermented sludge at a spectrum from 200 to 750 nm. Advanced pattern recognition methods are then used to map the non-linear relationship between measured absorption spectra to laboratory measurements of organic acid concentrations. Linear discriminant analysis, generalized discriminant analysis (GerDA), support vector machines (SVM), relevance vector machines, random forest and neural networks are investigated for this purpose and their performance compared. To validate the approach, online measurements have been taken at a full-scale 1.3-MW industrial biogas plant. Results show that whereas some of the methods considered do not yield satisfactory results, accurate prediction of organic acid concentration ranges can be obtained with both GerDA and SVM-based classifiers, with classification rates in excess of 87% achieved on test data.

Pure water vapor continuum measurements between 3100 and 4400 cm1: Evidence for water dimer absorption in near atmospheric conditions

Despite the potentially important role that water dimers may play in the Earth’s energy balance, there is still a lack of firm evidence for absorption of radiation by dimers in near-atmospheric conditions. We present results of the first high-resolution laboratory measurements of the water vapor continuum absorption within the 3100–4400 cm1 spectral region at a range of near-room temperatures. The analysis indicates a large contribution of dimer absorption to the water vapor continuum, significantly in excess of that predicted by other modern representations of the continuum. The temperature dependence agrees well with that expected for dimers.

Comparison of emission rate values for odour and odorous chemicals derived from two sampling devices

Field and laboratory measurements identified a complex relationship between odour emission rates provided by the US EPA dynamic emission chamber and the University of New South Wales wind tunnel. Using a range of model compounds in an aqueous odour source, we demonstrate that emission rates derived from the wind tunnel and flux chamber are a function of the solubility of the materials being emitted, the concentrations of the materials within the liquid; and the aerodynamic conditions within the device – either velocity in the wind tunnel, or flushing rate for the flux chamber. The ratio of wind tunnel to flux chamber odour emission rates (OU m-2 s) ranged from about 60:1 to 112:1. The emission rates of the model odorants varied from about 40:1 to over 600:1. These results may provide, for the first time, a basis for the development of a model allowing an odour emission rate derived from either device to be used for odour dispersion modelling.

Investigation of the effect of organics on the water uptake of marine aerosols

Water uptake refers to the ability of atmospheric particles to take up water vapour from the surrounding atmosphere. This is an important property that affects particle size and phase and therefore influences many characteristics of aerosols relevant to air quality and climate. However, the water uptake properties of many important atmospheric aerosol systems, including those related to the oceans, are still not fully understood. Therefore, the primary aim of this PhD research program was to investigate the water uptake properties of marine aerosols. In particular, the effect of organics on marine aerosol water uptake was investigated. Field campaigns were conducted at remote coastal sites on the east coast of Australia (Agnes Water; March-April 2007) and west coast of Ireland (Mace Head; June 2007), and laboratory measurements were performed on bubble-generated sea spray aerosols. A combined Volatility-Hygroscopicity-Tandem Differential Mobility Analyser (VH-TDMA) was employed in all experiments. This system probes the changes in the hygroscopic properties of nanoparticles as volatile organic components are progressively evaporated. It also allows particle composition to be inferred from combined volatility-hygroscopicity measurements. Frequent new particle formation and growth events were observed during the Agnes Water campaign. The VH-TDMA was used to investigate freshly nucleated particles (17-22.5 nm) and it was found that the condensation of sulphate and/or organic vapours was responsible for driving particle growth during the events. Aitken mode particles (~40 nm) were also measured with the VH-TDMA. In 3 out of 18 VH-TDMA scans evaporation of a volatile, organic component caused a very large increase in hygroscopicity that could only be explained by an increase in the absolute water uptake of the particle residuals, and not merely an increase in their relative hygroscopicity. This indicated the presence of organic components that were suppressing the hygroscopic growth of mixed particles on the timescale of humidification in the VH-TDMA (6.5 secs). It was suggested that the suppression of water uptake was caused by either a reduced rate of hygroscopic growth due to the presence of organic films, or organic-inorganic interactions in solution droplets that had a negative effect on hygroscopicity. Mixed organic-inorganic particles were rarely observed by the VH-TDMA during the summer campaign conducted at Mace Head. The majority of particles below 100 nm in clean, marine air appeared to be sulphates neutralised to varying degrees by ammonia. On one unique day, 26 June 2007, particularly large concentrations of sulphate aerosol were observed and identified as volcanic emissions from Iceland. The degree of neutralisation of the sulphate aerosol by ammonia was calculated by the VH-TDMA and found to compare well with the same quantity measured by an aerosol mass spectrometer. This was an important verification of the VH-TMDA‘s ability to identify ammoniated sulphate aerosols based on the simultaneous measurement of aerosol volatility and hygroscopicity. A series of measurements were also conducted on sea spray aerosols generated from Moreton Bay seawater samples in a laboratory-based bubble chamber. Accumulation mode sea spray particles (38-173 nm) were found to contain only a minor organic fraction (< 10%) that had little effect on particle hygroscopicity. These results are important because previous studies have observed that accumulation mode sea spray particles are predominantly organic (~80% organic mass fraction). The work presented here suggests that this is not always the case, and that there may be currently unknown factors that are controlling the transfer of organics to the aerosol phase during the bubble bursting process. Taken together, the results of this research program have significantly improved our understanding of organic-containing marine aerosols and the way they interact with water vapour in the atmosphere.

Characterization of frequency‐dependent magnetic susceptibility in UXO electromagnetic geophysics

Soils at many locations that have their origin in volcanic parent material and have undergone extensive weathering often exhibit strong frequency-dependent magnetic susceptibilities. The presence of such susceptibility has a profound effect on electromagnetic induction data acquired in such environments. Their transient electromagnetic response is characterized by a t-1 decay that is strong enough to mask UXO responses. In a field study and associated laboratory work on characterizing the frequency-dependent magnetic susceptibility and its influence on transient electromagnetic data, we collected soil samples on the surface and in soil pits from the Island of Kaho'olawe, Hawaii, and measured their frequency dependent magnetic susceptibilities. We present the details of the field investigation, confirm previous theoretical work with field and laboratory measurements, characterize the susceptibility with a Cole-Cole model, and investigate the response specific to the measured susceptibility.

Suomen maankuoren kivilajikoostumuksen tulkinta seismisistä aineistoista

Aim of this study is to investigate composition of the crust in Finland using seismic wide-angle velocity models and laboratory measurements on P- and S-wave velocities of different rock types. The velocities adopted from wide-angle velocity models were compared with laboratory velocities of different rock types corrected for the crustal PT conditions in the study area. The wide-angle velocity models indicate that the P-wave velocity does not only increase step-wise at boundaries of major crustal layers, but there is also gradual increase of velocity within the layers. On the other hand, the laboratory measurements of velocities indicate that no single rock type is able to provide the gradual downward increasing trends. Thus, there must be gradual vertical changes in rock composition. The downward increase of velocities indicates that the composition of the crust becomes gradually more mafic with increasing depth. Even though single rock types cannot simulate the wide-angle model velocities, it can be done with a mixture of rock types. There are a large number of rock type mixtures giving the correct P-wave velocities. Therefore, the inverse solution of rock types and their proportions from velocities is a non-unique problem if only P-wave velocities is available. Amount of the possible rock type mixtures can be limitted using S-wave velocities, reflection seismic results and other geological and geophysical results of the study area. Crustal model FINMIX-2 is presented in this study and it suggest that the crustal velocity profiles can be simulated with rock type mixtures, where the upper crust consists of felsic gneisses and granitic-granodioritic rocks with a minor contribution of quartzite, amphibolite and diabase. In the middle crust the amphibolite proportion increases. The lower crust consists of tonalitic gneiss, mafic garnet granulite, hornblendite, pyroxenite and minor mafic eclogite. This composition model is in agreement with deep crustal kimberlite-hosted xenolith data in eastern Finland and reflectivity of the FIRE (Finnish Reflection Experiment). According to FINMIX-2 model the Moho is deeper and the crustal composition is a more mafic than an average global continental model would suggest. Composition models of southern Finland are quite similar than FINMIX-2 model. However, there are minor differencies between the models, which indicates areal differences of composition. Models of northern Finland shows that the crustal thickness is smaller than southern Finland and composition of the upper crust is different. Density profiles calculated from the lithological models suggest that there is practically no density contrast at Moho in areas of the high-velocity lower crust. This implies that crustal thickness in the central Fennoscandian Shield may have been controlled by the densities of the lower crustal and upper mantle rocks.

Simulation blocks for TOSSIM-T2

We develop several hardware and software simulation blocks for the TinyOS-2 (TOSSIM-T2) simulator. The choice of simulated hardware platform is the popular MICA2 mote. While the hardware simulation elements comprise of radio and external flash memory, the software blocks include an environment noise model, packet delivery model and an energy estimator block for the complete system. The hardware radio block uses the software environment noise model to sample the noise floor. The packet delivery model is built by establishing the SNR-PRR curve for the MICA2 system. The energy estimator block models energy consumption by Micro Controller Unit(MCU), Radio, LEDs, and external flash memory. Using the manufacturerpsilas data sheets we provide an estimate of the energy consumed by the hardware during transmission, reception and also track several of the MCUs states with the associated energy consumption. To study the effectiveness of this work, we take a case study of a paper presented in [1]. We obtain three sets of results for energy consumption through mathematical analysis, simulation using the blocks built into PowerTossim-T2 and finally laboratory measurements. Since there is a significant match between these result sets, we propose our blocks for T2 community to effectively test their application energy requirements and node life times.

Monte Carlo Simulations of Molecular Clusters in Nucleation

A better understanding of the limiting step in a first order phase transition, the nucleation process, is of major importance to a variety of scientific fields ranging from atmospheric sciences to nanotechnology and even to cosmology. This is due to the fact that in most phase transitions the new phase is separated from the mother phase by a free energy barrier. This barrier is crossed in a process called nucleation. Nowadays it is considered that a significant fraction of all atmospheric particles is produced by vapor-to liquid nucleation. In atmospheric sciences, as well as in other scientific fields, the theoretical treatment of nucleation is mostly based on a theory known as the Classical Nucleation Theory. However, the Classical Nucleation Theory is known to have only a limited success in predicting the rate at which vapor-to-liquid nucleation takes place at given conditions. This thesis studies the unary homogeneous vapor-to-liquid nucleation from a statistical mechanics viewpoint. We apply Monte Carlo simulations of molecular clusters to calculate the free energy barrier separating the vapor and liquid phases and compare our results against the laboratory measurements and Classical Nucleation Theory predictions. According to our results, the work of adding a monomer to a cluster in equilibrium vapour is accurately described by the liquid drop model applied by the Classical Nucleation Theory, once the clusters are larger than some threshold size. The threshold cluster sizes contain only a few or some tens of molecules depending on the interaction potential and temperature. However, the error made in modeling the smallest of clusters as liquid drops results in an erroneous absolute value for the cluster work of formation throughout the size range, as predicted by the McGraw-Laaksonen scaling law. By calculating correction factors to Classical Nucleation Theory predictions for the nucleation barriers of argon and water, we show that the corrected predictions produce nucleation rates that are in good comparison with experiments. For the smallest clusters, the deviation between the simulation results and the liquid drop values are accurately modelled by the low order virial coefficients at modest temperatures and vapour densities, or in other words, in the validity range of the non-interacting cluster theory by Frenkel, Band and Bilj. Our results do not indicate a need for a size dependent replacement free energy correction. The results also indicate that Classical Nucleation Theory predicts the size of the critical cluster correctly. We also presents a new method for the calculation of the equilibrium vapour density, surface tension size dependence and planar surface tension directly from cluster simulations. We also show how the size dependence of the cluster surface tension in equimolar surface is a function of virial coefficients, a result confirmed by our cluster simulations.

Monte Carlo simulations of molecular clusters in nucleation

A better understanding of the limiting step in a first order phase transition, the nucleation process, is of major importance to a variety of scientific fields ranging from atmospheric sciences to nanotechnology and even to cosmology. This is due to the fact that in most phase transitions the new phase is separated from the mother phase by a free energy barrier. This barrier is crossed in a process called nucleation. Nowadays it is considered that a significant fraction of all atmospheric particles is produced by vapor-to liquid nucleation. In atmospheric sciences, as well as in other scientific fields, the theoretical treatment of nucleation is mostly based on a theory known as the Classical Nucleation Theory. However, the Classical Nucleation Theory is known to have only a limited success in predicting the rate at which vapor-to-liquid nucleation takes place at given conditions. This thesis studies the unary homogeneous vapor-to-liquid nucleation from a statistical mechanics viewpoint. We apply Monte Carlo simulations of molecular clusters to calculate the free energy barrier separating the vapor and liquid phases and compare our results against the laboratory measurements and Classical Nucleation Theory predictions. According to our results, the work of adding a monomer to a cluster in equilibrium vapour is accurately described by the liquid drop model applied by the Classical Nucleation Theory, once the clusters are larger than some threshold size. The threshold cluster sizes contain only a few or some tens of molecules depending on the interaction potential and temperature. However, the error made in modeling the smallest of clusters as liquid drops results in an erroneous absolute value for the cluster work of formation throughout the size range, as predicted by the McGraw-Laaksonen scaling law. By calculating correction factors to Classical Nucleation Theory predictions for the nucleation barriers of argon and water, we show that the corrected predictions produce nucleation rates that are in good comparison with experiments. For the smallest clusters, the deviation between the simulation results and the liquid drop values are accurately modelled by the low order virial coefficients at modest temperatures and vapour densities, or in other words, in the validity range of the non-interacting cluster theory by Frenkel, Band and Bilj. Our results do not indicate a need for a size dependent replacement free energy correction. The results also indicate that Classical Nucleation Theory predicts the size of the critical cluster correctly. We also presents a new method for the calculation of the equilibrium vapour density, surface tension size dependence and planar surface tension directly from cluster simulations. We also show how the size dependence of the cluster surface tension in equimolar surface is a function of virial coefficients, a result confirmed by our cluster simulations.

Svekofennialaisen orogenian ylityönnön terminen mallinnus Itä-Suomen arkeeisella kratonilla

This study reviews the thermophysical properties of rocks and the lithosphere and describes a one dimensional thermal numeric model of hypothetical 20 km thick overthrust plate obtruded on to the Archeaen craton in the Svecofennian orogeny (1.92-1.77 Ga). The objective is to find out if the overthrust plate and its radiogenic heat sources were able to produce the thermal effects observed on the current erosion level of the Archaean craton. Heat transfer in lithosphere is assumed conductive, and advective heat transfer due to melting and melt transfer is supposed negligible. The study area is located in the Eastern Finland, approximately current Kainuu and Northern Karelia regions, east from the most active orogenic belt (Raahe-Ladoga zone), so that orogenic magmatism can be neglected. Physical parameters and boundary conditions for the model are from different earlier published sources: deep seismic profiles (rock variation in depth), laboratory measurements (heat production and conductivity of rocks), field measurements (heat flow densities), and pT(t) estimations from the Finnish precambrian to estimate the size and thickness of the sheet. Comparison of the modelling results to previous K-Ar datings and other pTt estimations show, that the effect of the overthrust sheet has been adequate (max. T 450°C at 4 kbar) to produce the K-Ar resetting ages measured from the Archaean bedrock at current erosion level. No other kind of thermal activation in lithosphere is required. Results show possibly very minor partial melting in upper middle crust underneath the overthrust sheet.

Microbiological characterisation of soils : Evaluation of some critical steps in data collection and experimental design

The study of soil microbiota and their activities is central to the understanding of many ecosystem processes such as decomposition and nutrient cycling. The collection of microbiological data from soils generally involves several sequential steps of sampling, pretreatment and laboratory measurements. The reliability of results is dependent on reliable methods in every step. The aim of this thesis was to critically evaluate some central methods and procedures used in soil microbiological studies in order to increase our understanding of the factors that affect the measurement results and to provide guidance and new approaches for the design of experiments. The thesis focuses on four major themes: 1) soil microbiological heterogeneity and sampling, 2) storage of soil samples, 3) DNA extraction from soil, and 4) quantification of specific microbial groups by the most-probable-number (MPN) procedure. Soil heterogeneity and sampling are discussed as a single theme because knowledge on spatial (horizontal and vertical) and temporal variation is crucial when designing sampling procedures. Comparison of adjacent forest, meadow and cropped field plots showed that land use has a strong impact on the degree of horizontal variation of soil enzyme activities and bacterial community structure. However, regardless of the land use, the variation of microbiological characteristics appeared not to have predictable spatial structure at 0.5-10 m. Temporal and soil depth-related patterns were studied in relation to plant growth in cropped soil. The results showed that most enzyme activities and microbial biomass have a clear decreasing trend in the top 40 cm soil profile and a temporal pattern during the growing season. A new procedure for sampling of soil microbiological characteristics based on stratified sampling and pre-characterisation of samples was developed. A practical example demonstrated the potential of the new procedure to reduce the analysis efforts involved in laborious microbiological measurements without loss of precision. The investigation of storage of soil samples revealed that freezing (-20 °C) of small sample aliquots retains the activity of hydrolytic enzymes and the structure of the bacterial community in different soil matrices relatively well whereas air-drying cannot be recommended as a storage method for soil microbiological properties due to large reductions in activity. Freezing below -70 °C was the preferred method of storage for samples with high organic matter content. Comparison of different direct DNA extraction methods showed that the cell lysis treatment has a strong impact on the molecular size of DNA obtained and on the bacterial community structure detected. An improved MPN method for the enumeration of soil naphthalene degraders was introduced as an alternative to more complex MPN protocols or the DNA-based quantification approach. The main advantage of the new method is the simple protocol and the possibility to analyse a large number of samples and replicates simultaneously.

A Microscopic Theory of the f.c.c.-b.c.c. Interface .

We present a first-principles theory of the equilibrium b.c.c.-f.c.c. interface at coexistence using the density functional method. We assume that the interfacial region has local body-centred tetragonal (b.c.t.) symmetry and predict typical interfacial widths to be of order 2 to 3 lattice spacings with typical energies close to 0.05 J/m2. These quantities are in good agreement with laboratory measurements on coherent interfaces.

Simulation Blocks for TOSSIM-T2"

We develop several hardware and software simulation blocks for the TinyOS-2 (TOSSIM-T2) simulator. The choice of simulated hardware platform is the popular MICA2 mote. While the hardware simulation elements comprise of radio and external flash memory, the software blocks include an environment noise model, packet delivery model and an energy estimator block for the complete system. The hardware radio block uses the software environment noise model to sample the noise floor.The packet delivery model is built by establishing the SNR-PRR curve for the MICA2 system. The energy estimator block models energy consumption by Micro Controller Unit(MCU), Radio,LEDs, and external flash memory. Using the manufacturer’s data sheets we provide an estimate of the energy consumed by the hardware during transmission, reception and also track several of the MCUs states with the associated energy consumption. To study the effectiveness of this work, we take a case study of a paper presented in [1]. We obtain three sets of results for energy consumption through mathematical analysis, simulation using the blocks built into PowerTossim-T2 and finally laboratory measurements. Since there is a significant match between these result sets, we propose our blocks for T2 community to effectively test their application energy requirements and node life times.

Sarcopenia em pacientes idosos com doença renal crônica em hemodiálise

Os pacientes idosos em hemodiálise (HD) são altamente suscetíveis ao desenvolvimento de sarcopenia, devido ao processo natural de envelhecimento e ao catabolismo induzido pelo procedimento de HD. O objetivo deste estudo foi primeiro, avaliar a prevalência de sarcopenia, dinapenia e atrofia muscular em um grupo de pacientes idosos em HD; o segundo, avaliar se os critérios aplicados para o diagnóstico de sarcopenia, propostos por sociedades internacionais, são capazes de distinguir os pacientes com pior condição clínica, estado nutricional e qualidade de vida. Este estudo multicêntrico e transversal incluiu 94 pacientes idosos em HD (> 60 anos) de cinco centros de diálise. Todos os participantes foram submetidos à avaliação antropométrica, de composição corporal, força de preensão manual (FPM), laboratorial e avaliação da condição nutricional através da avaliação subjetiva global de 7 pontos (AGS-7p). Adicionalmente, os participantes responderam um questionário de qualidade de vida. Para o diagnóstico de sarcopenia, foram adotados os critérios propostos por sociedades internacionais, que englobam os parâmetros indicativos de baixa massa muscular e baixa função muscular. Para a massa muscular adotou-se o índice de massa magra (IMM) < percentil 20 para o gênero e faixa etária de uma população de referência, avaliado a partir da massa magra obtida pelo somatório de dobras cutâneas. Para o critério de baixa função muscular, adotou-se a FPM < percentil 10 para o gênero, faixa etária e o braço utilizado de uma população de referência. Os pacientes foram classificados como Sarcopênicos (baixo IMM associado à baixa FPM); Dinapênicos (baixa FPM) e Atrofia muscular (baixo IMM). A sarcopenia estava presente em 13.8% dos pacientes, enquanto a dinapenia foi observada em 37.2% e a atrofia muscular em 35.1%. A sarcopenia foi capaz de distinguir os pacientes que possuíam maior comprometimento do estado nutricional e da composição corporal. O critério de função muscular (isoladamente ou em combinação com a massa muscular) também identificou os pacientes com pior qualidade de vida. Em conclusão, a prevalência de sarcopenia foi observada em 13,8% do grupo. Entretanto, ao usar apenas critérios que indicam redução da força ou massa muscular, esta prevalência aumentou para 30%. A condição de sarcopenia distinguiu pacientes com pior estado nutricional e qualidade de vida.

Interaction between tsunami waves and isolated conical Islands

A newly developed computer model, which solves the horizontal two-dimensional Boussinesq equations using a total variation diminishing Lax-Wendroff scheme, has been used to study the runup of solitary waves, with various heights, on idealized conical islands consisting of side slopes of different angles. This numerical model has first been validated against high-quality laboratory measurements of solitary wave runups on a uniform plane slope and on an isoliated conical island, with satisfactory agreement being achieved. An extensive parametric study concerning the effects of the wave height and island slope on the solitary wave runup has subsequently been carried out. Strong wave shoaling and diffraction effects have been observed for all the cases investigated. The relationship between the runup height and wave height has been obtained and compared with that for the case on uniform plane slopes. It has been found that the runup on a conical island is generally lower than that on a uniform plane slope, as a result of the two-dimensional effect. The correlation between the runup with the side slope of an island has also been identified, with higher runups on milder slopes. This comprehensive study on the soliton runup on islands is relevant to the protection of coastal and inland regions from extreme wave attacks. © the Coastal Education & Research Foundation 2012.