Quantification of air quality in space and time and its effects on health

The long-term adverse effects on health associated with air pollution exposure can be estimated using either cohort or spatio-temporal ecological designs. In a cohort study, the health status of a cohort of people are assessed periodically over a number of years, and then related to estimated ambient pollution concentrations in the cities in which they live. However, such cohort studies are expensive and time consuming to implement, due to the long-term follow up required for the cohort. Therefore, spatio-temporal ecological studies are also being used to estimate the long-term health effects of air pollution as they are easy to implement due to the routine availability of the required data. Spatio-temporal ecological studies estimate the health impact of air pollution by utilising geographical and temporal contrasts in air pollution and disease risk across $n$ contiguous small-areas, such as census tracts or electoral wards, for multiple time periods. The disease data are counts of the numbers of disease cases occurring in each areal unit and time period, and thus Poisson log-linear models are typically used for the analysis. The linear predictor includes pollutant concentrations and known confounders such as socio-economic deprivation. However, as the disease data typically contain residual spatial or spatio-temporal autocorrelation after the covariate effects have been accounted for, these known covariates are augmented by a set of random effects. One key problem in these studies is estimating spatially representative pollution concentrations in each areal which are typically estimated by applying Kriging to data from a sparse monitoring network, or by computing averages over modelled concentrations (grid level) from an atmospheric dispersion model. The aim of this thesis is to investigate the health effects of long-term exposure to Nitrogen Dioxide (NO2) and Particular matter (PM10) in mainland Scotland, UK. In order to have an initial impression about the air pollution health effects in mainland Scotland, chapter 3 presents a standard epidemiological study using a benchmark method. The remaining main chapters (4, 5, 6) cover the main methodological focus in this thesis which has been threefold: (i) how to better estimate pollution by developing a multivariate spatio-temporal fusion model that relates monitored and modelled pollution data over space, time and pollutant; (ii) how to simultaneously estimate the joint effects of multiple pollutants; and (iii) how to allow for the uncertainty in the estimated pollution concentrations when estimating their health effects. Specifically, chapters 4 and 5 are developed to achieve (i), while chapter 6 focuses on (ii) and (iii). In chapter 4, I propose an integrated model for estimating the long-term health effects of NO2, that fuses modelled and measured pollution data to provide improved predictions of areal level pollution concentrations and hence health effects. The air pollution fusion model proposed is a Bayesian space-time linear regression model for relating the measured concentrations to the modelled concentrations for a single pollutant, whilst allowing for additional covariate information such as site type (e.g. roadside, rural, etc) and temperature. However, it is known that some pollutants might be correlated because they may be generated by common processes or be driven by similar factors such as meteorology. The correlation between pollutants can help to predict one pollutant by borrowing strength from the others. Therefore, in chapter 5, I propose a multi-pollutant model which is a multivariate spatio-temporal fusion model that extends the single pollutant model in chapter 4, which relates monitored and modelled pollution data over space, time and pollutant to predict pollution across mainland Scotland. Considering that we are exposed to multiple pollutants simultaneously because the air we breathe contains a complex mixture of particle and gas phase pollutants, the health effects of exposure to multiple pollutants have been investigated in chapter 6. Therefore, this is a natural extension to the single pollutant health effects in chapter 4. Given NO2 and PM10 are highly correlated (multicollinearity issue) in my data, I first propose a temporally-varying linear model to regress one pollutant (e.g. NO2) against another (e.g. PM10) and then use the residuals in the disease model as well as PM10, thus investigating the health effects of exposure to both pollutants simultaneously. Another issue considered in chapter 6 is to allow for the uncertainty in the estimated pollution concentrations when estimating their health effects. There are in total four approaches being developed to adjust the exposure uncertainty. Finally, chapter 7 summarises the work contained within this thesis and discusses the implications for future research.

Abordagem rítmica probabilística aplicada em análises de fluxo e estabilidade de taludes

Tese (doutorado)—Universidade de Brasília, Faculdade de Tecnologia, Programa de Pós-Graduação em Geotecnia, 2015.

Caractérisation du fonctionnement d'une hydrolienne à membrane ondulante pour la récupération de l'énergie des courants marins

This manuscript presents three approaches : analytical, experimental and numerical, to study the behaviour of a flexible membrane tidal energy converter. This technology, developed by the EEL Energy company, is based on periodic deformations of a pre-stressed flexible structure. Energy converters, located on each side of the device, are set into motion by the wave-like motion. In the analytical model, the membrane is represented by a linear beam model at one dimension and the flow by a 3 dimensions potential fluid. The fluid forces are evaluated by the elongated body theory. Energy is dissipated all over the length of the membrane. A 20th scale experimental prototype has been designed with micro-dampers to simulate the power take-off. Trials have allowed to validate the undulating membrane energy converter concept. A numerical model has been developed. Each element of the device is represented and the energy dissipation is done by dampers element with a damping law linear to damper velocity. Comparison of the three approaches validates their ability to represent the membrane behaviour without damping. The energy dissipation applied with the analytical model is clearly different from the two other models because of the location (where the energy is dissipated) and damping law. The two others show a similar behaviour and the same order of power take off repartition but value of power take off are underestimated by the numerical model. This three approaches have allowed to put forward key-parameters on which depend the behaviour of the membrane and the parametric study highlights the complementarity and the advantage of developing three approaches in parallel to answer industrial optimization problems. To make the link between trials in flume tank and sea trials, a 1/6th prototype has been built. To do so, the change of scale was studied. The behaviour of both prototypes is compared and differences could be explained by differences of boundary conditions and confinement effects. To evaluated membrane long-term behaviour at sea, a method of ageing accelerated by temperature and fatigue tests have been carried out on prototype materials samples submerged in sea water.

Pore size distributions derived from adsorption isotherms, immersion calorimetry, and isosteric heats: A comparative study

We compare the pore size distribution of a well-characterized activated carbon derived from model-dependent, adsorption integral equation (AIE) methods with those from model-independent, immersion calorimetry and isosteric heat analyses. The AIE approach applied to nitrogen gave a mean pore width of 0.57 nm; the CO2 distribution exhibited wider dispersion. Spherical model application to CO2 and diffusion limitations for nitrogen and argon were proposed as primary reasons for inconsistency. Immersion enthalpy revealed a sharp decrease in available area equivalent to a cut-off due to molecular exclusion when the accessible surface was assessed against probe kinetic diameter. Mean pore width was identified as 0.58 ± 0.02 nm, endorsing the underlying assumptions for the nitrogen-based AIE approach. A comparison of the zero-coverage isosteric heat of adsorption for various non-polar adsorptives by the porous test sample was compared with the same adsorptives in contact with a non-porous reference adsorbent, leading to an energy ratio or adsorption enhancement factor. A linear relationship between the energy ratio and probe kinetic diameter indicated a primary pore size at 0.59 nm. The advantage of this enthalpy, model-independent methods over AIE were due to no assumptions regarding probe molecular shape, and no assumptions for pore shape and/or connectivity.

Estudo da perda de carga e calor no poço injetor no processo De drenagem gravitacional assistida por vapor (SAGD)

The oil companies in the area in general are looking for new technologies that can increase the recovery factor of oil contained in reservoirs. These investments are mainly aimed at reducing the costs of projects which are high. Steam injection is one of these special methods of recovery in which steam is injected into the reservoir in order to reduce the viscosity of the oil and make it more mobile. The process assisted gravity drainage steam (SAGD) using steam injection in its mechanism, as well as two parallel horizontal wells. In this process steam is injected through the horizontal injection well, then a vapor chamber is formed by heating the oil in the reservoir and, by the action of gravitational forces, this oil is drained down to where the production well. This study aims to analyze the influence of pressure drop and heat along the injection well in the SAGD process. Numerical simulations were performed using the thermal simulator STARS of CMG (Computer Modeling Group). The parameters studied were the thermal conductivity of the formation, the flow of steam injection, the inner diameter of the column, the steam quality and temperature. A factorial design was used to verify the influence of the parameters studied in the recovery factor. We also analyzed different injection flow rates for the model with pressure drop and no pressure drop, as well as different maximum flow rates of oil production. Finally, we performed an economic analysis of the two models in order to check the profitability of the projects studied. The results showed that the pressure drop in injection well have a significant influence on the SAGD process.

On the Stability of Supersonic Boundary Layers with Injection

The problem of supersonic flow over a 5 degree half-angle cone with injection of gas through a porous section on the body into the boundary layer is studied experimentally. Three injected gases are used: helium, nitrogen, and RC318 (octafluorocyclobutane). Experiments are performed in a Mach 4 Ludwieg tube with nitrogen as the free stream gas. Shaping of the injector section relative to the rest of the body is found to admit a "tuned" injection rate which minimizes the strength of shock waves formed by injection. A high-speed schlieren imaging system with a framing rate of 290 kHz is used to study the instability in the region of flow downstream of injection, referred to as the injection layer. This work provides the first experimental data on the wavelength, convective speed, and frequency of the instability in such a flow. The stability characteristics of the injection layer are found to be very similar to those of a free shear layer. The findings of this work present a new paradigm for future stability analyses of supersonic flow with injection.

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.

Breaking the curve with candels: a Bayesian approach to reveal the non-universality of the dust-attenuation law at high redshift

Dust attenuation affects nearly all observational aspects of galaxy evolution, yet very little is known about the form of the dust-attenuation law in the distant universe. Here, we model the spectral energy distributions of galaxies at z ~ 1.5–3 from CANDELS with rest-frame UV to near-IR imaging under different assumptions about the dust law, and compare the amount of inferred attenuated light with the observed infrared (IR) luminosities. Some individual galaxies show strong Bayesian evidence in preference of one dust law over another, and this preference agrees with their observed location on the plane of infrared excess (IRX, L_TIR/L_UV) and UV slope (β). We generalize the shape of the dust law with an empirical model, A_ λ,σ =E(B-V)k_ λ (λ / λ v)^ σ where k_λ is the dust law of Calzetti et al., and show that there exists a correlation between the color excess E(B-V) and tilt δ with δ =(0.62±0.05)log(E(B-V))+(0.26±0.02). Galaxies with high color excess have a shallower, starburst-like law, and those with low color excess have a steeper, SMC-like law. Surprisingly, the galaxies in our sample show no correlation between the shape of the dust law and stellar mass, star formation rate, or β. The change in the dust law with color excess is consistent with a model where attenuation is caused by scattering, a mixed star–dust geometry, and/or trends with stellar population age, metallicity, and dust grain size. This rest-frame UV-to-near-IR method shows potential to constrain the dust law at even higher redshifts (z>3).

Multifactorial analysis of changes in body mass index across the adult life course: a study with 65 years of follow-up

Background: Although the negative consequences on health of being obese are well known, most adults gain weight across the lifespan. The general increase in body mass index (BMI) is mainly considered to originate from behavioral and environmental changes; however, few studies have evaluated the influence of these factors on change in BMI in the presence of genetic risk. We aimed to study the influence of multifactorial causes of change in BMI, over 65 years. Methods and Findings: Totally, 6130 participants from TwinGene, who had up to five assessments, and 536 from the Swedish Adoption/Twin Study of Aging, who had up to 12 assessments, ranging over 65 years were included. The influence of lifestyle factors, birth cohort, cardiometabolic diseases and an individual obesity genetic risk score (OGRS) based on 32 single nucleotide polymorphisms on change in BMI was evaluated with a growth model. For both sexes, BMI increased from early adulthood to age of 65 years, after which the increase leveled off; BMI declined after age of 80 years. A higher OGRS, birth after 1925 and cardiometabolic diseases were associated with higher average BMI and a steeper increase in BMI prior to 65 years of age. Among men, few factors were identified that influence BMI trajectories in late life, whereas for women type 2 diabetes mellitus and dementia were associated with a steeper decrease in BMI after the age of 65 years. Conclusions: There are two turning points in BMI in late adulthood, one at the age of 65 years and one at the age 80 years. Factors associated with an increase in BMI in midlife were not associated with an increase in BMI after the age of 65 years. These findings indicate that the causes and consequences of change in BMI differ across the lifespan. Current health recommendations need to be adjusted accordingly.

Applying fuzzy logic to comparative distribution modelling: a case study with two sympatric amphibians.

We modelled the distributions of two toads (Bufo bufo and Epidalea calamita) in the Iberian Peninsula using the favourability function, which makes predictions directly comparable for different species and allows fuzzy logic operations to relate different models. The fuzzy intersection between individual models, representing favourability for the presence of both species simultaneously, was compared with another favourability model built on the presences shared by both species. The fuzzy union between individual models, representing favourability for the presence of any of the two species, was compared with another favourabilitymodel based on the presences of either or both of them. The fuzzy intersections between favourability for each species and the complementary of favourability for the other (corresponding to the logical operation “A and not B”) were compared with models of exclusive presence of one species versus the exclusive presence of the other. The results of modelling combined species data were highly similar to those of fuzzy logic operations between individual models, proving fuzzy logic and the favourability function valuable for comparative distribution modelling. We highlight several advantages of fuzzy logic over other forms of combining distribution models, including the possibility to combine multiple species models for management and conservation planning.

Above ground biomass functions with vegetation indices for multiple use systems of two evergreen oaks

Remote sensing is a promising approach for above ground biomass estimation, as forest parameters can be obtained indirectly. The analysis in space and time is quite straight forward due to the flexibility of the method to determine forest crown parameters with remote sensing. It can be used to evaluate and monitoring for example the development of a forest area in time and the impact of disturbances, such as silvicultural practices or deforestation. The vegetation indices, which condense data in a quantitative numeric manner, have been used to estimate several forest parameters, such as the volume, basal area and above ground biomass. The objective of this study was the development of allometric functions to estimate above ground biomass using vegetation indices as independent variables. The vegetation indices used were the Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI), Simple Ratio (SR) and Soil-Adjusted Vegetation Index (SAVI). QuickBird satellite data, with 0.70 m of spatial resolution, was orthorectified, geometrically and atmospheric corrected, and the digital number were converted to top of atmosphere reflectance (ToA). Forest inventory data and published allometric functions at tree level were used to estimate above ground biomass per plot. Linear functions were fitted for the monospecies and multispecies stands of two evergreen oaks (Quercus suber and Quercus rotundifolia) in multiple use systems, montados. The allometric above ground biomass functions were fitted considering the mean and the median of each vegetation index per grid as independent variable. Species composition as a dummy variable was also considered as an independent variable. The linear functions with better performance are those with mean NDVI or mean SR as independent variable. Noteworthy is that the two better functions for monospecies cork oak stands have median NDVI or median SR as independent variable. When species composition dummy variables are included in the function (with stepwise regression) the best model has median NDVI as independent variable. The vegetation indices with the worse model performance were EVI and SAVI.

A single nucleotide polymorphism in NEUROD1 is associated with production traits in nelore beef cattle.

Feed efficiency and carcass characteristics are late-measured traits. The detection of molecular markers associated with them can help breeding programs to select animals early in life, and to predict breeding values with high accuracy. The objective of this study was to identify polymorphisms in the functional and positional candidate gene NEUROD1 (neurogenic differentiation 1), and investigate their associations with production traits in reference families of Nelore cattle. A total of 585 steers were used, from 34 sires chosen to represent the variability of this breed. By sequencing 14 animals with extreme residual feed intake (RFI) values, seven single nucleotide polymorphisms (SNPs) in NEUROD1 were identified. The investigation of marker effects on the target traits RFI, backfat thickness (BFT), ribeye area (REA), average body weight (ABW), and metabolic body weight (MBW) was performed with a mixed model using the restricted maximum likelihood method. SNP1062, which changes cytosine for guanine, had no significant association with RFI or REA. However, we found an additive effect on ABW (P ≤ 0.05) and MBW (P ≤ 0.05), with an estimated allele substitution effect of -1.59 and -0.93 kg0.75, respectively. A dominant effect of this SNP for BFT was also found (P ≤ 0.010). Our results are the first that identify NEUROD1 as a candidate that affects BFT, ABW, and MBW. Once confirmed, the inclusion of this SNP in dense panels may improve the accuracy of genomic selection for these traits in Nelore beef cattle as this SNP is not currently represented on SNP chips.

Economic and Environmental Optimization of Ornamental Stone Block Cutting with a 3D Algorithm

This thesis focuses on finding the optimum block cutting dimensions in terms of the environmental and economic factors by using a 3D algorithm for a limestone quarry in Foggia, Italy. The environmental concerns of quarrying operations are mainly: energy consumption, material waste, and pollution. The main economic concerns are the block recovery, the selling prices, and the production costs. Fractures adversely affect the block recovery ratio. With a fracture model, block production can be optimized. In this research, the waste volume produced by quarrying was minimised to increase the recovery ratio and ensure economic benefits. SlabCutOpt is a software developed at DICAM–University of Bologna for block cutting optimization which tests different cutting angles on the x-y-z planes to offer up alternative cutting methods. The program tests several block sizes and outputs the optimal result for each entry. By using SlabCutOpt, ten different block dimensions were analysed, the results indicated the maximum number of non-intersecting blocks for each dimension. After analysing the outputs, the block named number 1 with the dimensions ‘1mx1mx1m’ had the highest recovery ratio as 43% and the total Relative Money Value (RMV) with a value of 22829. Dimension number 1, also had the lowest waste volume, with a value of 3953.25 m3, for the total bench. For cutting the total bench volume of 6932.25m3, the diamond wire cutter had the lowest dust emission values for the block with the dimension ‘2mx2mx2m’, with a value of 24m3. When compared with the Eco-Label standards, block dimensions having surface area values lower than 15m2, were found to fit the natural resource waste criteria of the label, as the threshold required 25% of minimum recovery [1]. Due to the relativity of production costs, together with the Eco-Label threshold, the research recommends the selection of the blocks with a surface area value between 6m2 and 14m2.

Urban air quality: new high-resolution modeling approaches and forecasting tools for citizens

Air pollution is one of the greatest health risks in the world. At the same time, the strong correlation with climate change, as well as with Urban Heat Island and Heat Waves, make more intense the effects of all these phenomena. A good air quality and high levels of thermal comfort are the big goals to be reached in urban areas in coming years. Air quality forecast help decision makers to improve air quality and public health strategies, mitigating the occurrence of acute air pollution episodes. Air quality forecasting approaches combine an ensemble of models to provide forecasts from global to regional air pollution and downscaling for selected countries and regions. The development of models dedicated to urban air quality issues requires a good set of data regarding the urban morphology and building material characteristics. Only few examples of air quality forecast system at urban scale exist in the literature and often they are limited to selected cities. This thesis develops by setting up a methodology for the development of a forecasting tool. The forecasting tool can be adapted to all cities and uses a new parametrization for vegetated areas. The parametrization method, based on aerodynamic parameters, produce the urban spatially varying roughness. At the core of the forecasting tool there is a dispersion model (urban scale) used in forecasting mode, and the meteorological and background concentration forecasts provided by two regional numerical weather forecasting models. The tool produces the 1-day spatial forecast of NO2, PM10, O3 concentration, the air temperature, the air humidity and BLQ-Air index values. The tool is automatized to run every day, the maps produced are displayed on the e-Globus platform, updated every day. The results obtained indicate that the forecasting output were in good agreement with the observed measurements.

Wind loading predictions using computational fluid dynamics

Using Computational Wind Engineering, CWE, for solving wind-related problems is still a challenging task today, mainly due to the high computational cost required to obtain trustworthy simulations. In particular, the Large Eddy Simulation, LES, has been widely used for evaluating wind loads on buildings. The present thesis assesses the capability of LES as a design tool for wind loading predictions through three cases. The first case is using LES for simulating the wind field around a ground-mounted rectangular prism in Atmospheric Boundary Layer (ABL) flow. The numerical results are validated with experimental results for seven wind attack angles, giving a global understanding of the model performance. The case with the worst model behaviour is investigated, including the spatial distribution of the pressure coefficients and their discrepancies with respect to experimental results. The effects of some numerical parameters are investigated for this case to understand their effectiveness in modifying the obtained numerical results. The second case is using LES for investigating the wind effects on a real high-rise building, aiming at validating the performance of LES as a design tool in practical applications. The numerical results are validated with the experimental results in terms of the distribution of the pressure statistics and the global forces. The mesh sensitivity and the computational cost are discussed. The third case is using LES for studying the wind effects on the new large-span roof over the Bologna stadium. The dynamic responses are analyzed and design envelopes for the structure are obtained. Although it is a numerical simulation before the traditional wind tunnel tests, i.e. the validation of the numerical results are not performed, the preliminary evaluations can effectively inform later investigations and provide the final design processes with deeper confidence regarding the absence of potentially unexpected behaviours.