Using stochastic volatility models to analyse weekly ozone averages in Mexico City

In this paper we make use of some stochastic volatility models to analyse the behaviour of a weekly ozone average measurements series. The models considered here have been used previously in problems related to financial time series. Two models are considered and their parameters are estimated using a Bayesian approach based on Markov chain Monte Carlo (MCMC) methods. Both models are applied to the data provided by the monitoring network of the Metropolitan Area of Mexico City. The selection of the best model for that specific data set is performed using the Deviance Information Criterion and the Conditional Predictive Ordinate method.

Meteomap: Generation Of Meteorological Variables For Distributed Physically-Based Hydrological Modeling

Distributed energy and water balance models require time-series surfaces of the meteorological variables involved in hydrological processes. Most of the hydrological GIS-based models apply simple interpolation techniques to extrapolate the point scale values registered at weather stations at a watershed scale. In mountainous areas, where the monitoring network ineffectively covers the complex terrain heterogeneity, simple geostatistical methods for spatial interpolation are not always representative enough, and algorithms that explicitly or implicitly account for the features creating strong local gradients in the meteorological variables must be applied. Originally developed as a meteorological pre-processing tool for a complete hydrological model (WiMMed), MeteoMap has become an independent software. The individual interpolation algorithms used to approximate the spatial distribution of each meteorological variable were carefully selected taking into account both, the specific variable being mapped, and the common lack of input data from Mediterranean mountainous areas. They include corrections with height for both rainfall and temperature (Herrero et al., 2007), and topographic corrections for solar radiation (Aguilar et al., 2010). MeteoMap is a GIS-based freeware upon registration. Input data include weather station records and topographic data and the output consists of tables and maps of the meteorological variables at hourly, daily, predefined rainfall event duration or annual scales. It offers its own pre and post-processing tools, including video outlook, map printing and the possibility of exporting the maps to images or ASCII ArcGIS formats. This study presents the friendly user interface of the software and shows some case studies with applications to hydrological modeling.

Implementação de um Módulo de Supervisão para um Sistema de Detecção de Vazamentos em Dutos de Petróleo

The transport of fluids through pipes is used in the oil industry, being the pipelines an important link in the logistics flow of fluids. However, the pipelines suffer deterioration in their walls caused by several factors which may cause loss of fluids to the environment, justifying the investment in techniques and methods of leak detection to minimize fluid loss and environmental damage. This work presents the development of a supervisory module in order to inform to the operator the leakage in the pipeline monitored in the shortest time possible, in order that the operator log procedure that entails the end of the leak. This module is a component of a system designed to detect leaks in oil pipelines using sonic technology, wavelets and neural networks. The plant used in the development and testing of the module presented here was the system of tanks of LAMP, and its LAN, as monitoring network. The proposal consists of, basically, two stages. Initially, assess the performance of the communication infrastructure of the supervisory module. Later, simulate leaks so that the DSP sends information to the supervisory performs the calculation of the location of leaks and indicate to which sensor the leak is closer, and using the system of tanks of LAMP, capture the pressure in the pipeline monitored by piezoresistive sensors, this information being processed by the DSP and sent to the supervisory to be presented to the user in real time

Stochastic simulation of time-series models combined with geostatistics to predict water-table scenarios in a Guarani Aquifer System outcrop area, Brazil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Genetic Structure of Populations of the Rice-Infecting Pathogen Rhizoctonia solani AG-1 IA from China

Sheath blight disease (SBD) on rice, caused by Rhizoctonia solani AG-1 IA, is one of the most devastating rice diseases on a global basis, including China (in Eastern Asia), the world's largest rice-growing country. We analyzed the population genetics of nine rice-infecting populations from China using nine microsatellite loci. One allopatric population from India (Southern Asia) was included in the analyses. In total, 300 different multilocus genotypes were found among 572 fungal isolates. Clonal fractions within rice fields were 16 to 95%, suggesting that sclerotia were a major source of primary inoculum in some fields. Global Phi(ST) statistics (Phi(ST) = 42.49; P <= 0.001) were consistent with a relatively high level of differentiation among populations overall; however, pairwise comparisons gave nonsignificant R(ST) values, consistent with contemporary gene flow among five of the populations. Four of these populations were located along the Yangtze River tributary network. Gene flow followed an isolation-by-distance model consistent with restricted long-distance migration. Historical migration rates were reconstructed and yielded values that explained the current levels of population subdivision. Except for one population which appeared to be strictly clonal, all populations showed evidence of a mixed reproductive mode, including both asexual and sexual reproduction. One population had a strictly recombining structure (all loci were in Hardy-Weinberg equilibrium) but the remaining populations from China and the one from India exhibited varying degrees of sexual reproduction. Six populations showed significant F(IS) values consistent with inbreeding.

Condições hidrológicas de uma bacia hidrográfica com diferentes usos do solo na região de taquaritinga, SP

In areas where human activities, as agriculture, are developed land use contributes to physical, chemical and biological characteristics of water. This study aimed to identify physical and chemical variations of the water monitoring network in 4 first magnitude watersheds with different land use/occupation in the hydrographic basin of the Glori Farm Creek in the city of Taquaritinga, State of Sao Paulo, from February to July, 2005. The methodology consisted of weekly sample water collections to naalyze turbidity, temperature, dissolved oxygen, pH and electric conductivity. The characterization of the hydrological conditions of the watersheds showed that agricultural activities, including the productive system of the sugarcane cultivation, have affected the quality of water resources.

Ações de educação ambiental para conservação e preservação do Córrego Bela Vista, Piacatu-SP

Pós-graduação em Geografia - FCT

Estudo para correção dos sensores da rede de monitoramento do campo elétrico em função do campo elétrico de tempo bom

A monitoring network of atmospheric electric field covering the Vale do Paraiba region was implemented. The sensors were located on different sites with different altitude and geographic topology. The present work reports the study carried on those sensors in order to verify the necessity of using some correction factor to the measured local electric field intensity due to effects of local environment. The measurements were done in continuous 24 hours per day with the data recorded on registers in each device accumulating information during a period of four months. The relation between the electric field values by each sensor was compared to the reference located on Sao Jose dos Campos city using the same period. In a graphical analysis using the local field intensity and the reference, the data were fitted to a straight line obtained by minimum square method. Variation up to 95% was observed between the field values in some sensors. Another method was proposed, comparing the mean values of the electric field in a function of time. The variation in some sensors reached up to 133%. We conclude that the variations are due to local atmospheric conditions and no correction factor is required on the electric field sensors

A non-homogeneous poisson model with spatial anisotropy applied to ozone data from Mexico City

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Stochastic simulation of time-series models combined with geostatistics to predict water-table scenarios in a Guarani Aquifer System outcrop area, Brazil

Stochastic methods based on time-series modeling combined with geostatistics can be useful tools to describe the variability of water-table levels in time and space and to account for uncertainty. Monitoring water-level networks can give information about the dynamic of the aquifer domain in both dimensions. Time-series modeling is an elegant way to treat monitoring data without the complexity of physical mechanistic models. Time-series model predictions can be interpolated spatially, with the spatial differences in water-table dynamics determined by the spatial variation in the system properties and the temporal variation driven by the dynamics of the inputs into the system. An integration of stochastic methods is presented, based on time-series modeling and geostatistics as a framework to predict water levels for decision making in groundwater management and land-use planning. The methodology is applied in a case study in a Guarani Aquifer System (GAS) outcrop area located in the southeastern part of Brazil. Communication of results in a clear and understandable form, via simulated scenarios, is discussed as an alternative, when translating scientific knowledge into applications of stochastic hydrogeology in large aquifers with limited monitoring network coverage like the GAS.

Numerical simulations of magma chamber dynamics at Campi Flegrei, and associated seismicity, deformation and gravity changes

Understanding the complex relationships between quantities measured by volcanic monitoring network and shallow magma processes is a crucial headway for the comprehension of volcanic processes and a more realistic evaluation of the associated hazard. This question is very relevant at Campi Flegrei, a volcanic quiescent caldera immediately north-west of Napoli (Italy). The system activity shows a high fumarole release and periodic ground slow movement (bradyseism) with high seismicity. This activity, with the high people density and the presence of military and industrial buildings, makes Campi Flegrei one of the areas with higher volcanic hazard in the world. In such a context my thesis has been focused on magma dynamics due to the refilling of shallow magma chambers, and on the geophysical signals detectable by seismic, deformative and gravimetric monitoring networks that are associated with this phenomenologies. Indeed, the refilling of magma chambers is a process frequently occurring just before a volcanic eruption; therefore, the faculty of identifying this dynamics by means of recorded signal analysis is important to evaluate the short term volcanic hazard. The space-time evolution of dynamics due to injection of new magma in the magma chamber has been studied performing numerical simulations with, and implementing additional features in, the code GALES (Longo et al., 2006), recently developed and still on the upgrade at the Istituto Nazionale di Geofisica e Vulcanologia in Pisa (Italy). GALES is a finite element code based on a physico-mathematical two dimensional, transient model able to treat fluids as multiphase homogeneous mixtures, compressible to incompressible. The fundamental equations of mass, momentum and energy balance are discretised both in time and space using the Galerkin Least-Squares and discontinuity-capturing stabilisation technique. The physical properties of the mixture are computed as a function of local conditions of magma composition, pressure and temperature.The model features enable to study a broad range of phenomenologies characterizing pre and sin-eruptive magma dynamics in a wide domain from the volcanic crater to deep magma feeding zones. The study of displacement field associated with the simulated fluid dynamics has been carried out with a numerical code developed by the Geophysical group at the University College Dublin (O’Brien and Bean, 2004b), with whom we started a very profitable collaboration. In this code, the seismic wave propagation in heterogeneous media with free surface (e.g. the Earth’s surface) is simulated using a discrete elastic lattice where particle interactions are controlled by the Hooke’s law. This method allows to consider medium heterogeneities and complex topography. The initial and boundary conditions for the simulations have been defined within a coordinate project (INGV-DPC 2004-06 V3_2 “Research on active volcanoes, precursors, scenarios, hazard and risk - Campi Flegrei”), to which this thesis contributes, and many researchers experienced on Campi Flegrei in volcanological, seismic, petrological, geochemical fields, etc. collaborate. Numerical simulations of magma and rock dynamis have been coupled as described in the thesis. The first part of the thesis consists of a parametric study aimed at understanding the eect of the presence in magma of carbon dioxide in magma in the convection dynamics. Indeed, the presence of this volatile was relevant in many Campi Flegrei eruptions, including some eruptions commonly considered as reference for a future activity of this volcano. A set of simulations considering an elliptical magma chamber, compositionally uniform, refilled from below by a magma with volatile content equal or dierent from that of the resident magma has been performed. To do this, a multicomponent non-ideal magma saturation model (Papale et al., 2006) that considers the simultaneous presence of CO2 and H2O, has been implemented in GALES. Results show that the presence of CO2 in the incoming magma increases its buoyancy force promoting convection ad mixing. The simulated dynamics produce pressure transients with frequency and amplitude in the sensitivity range of modern geophysical monitoring networks such as the one installed at Campi Flegrei . In the second part, simulations more related with the Campi Flegrei volcanic system have been performed. The simulated system has been defined on the basis of conditions consistent with the bulk of knowledge of Campi Flegrei and in particular of the Agnano-Monte Spina eruption (4100 B.P.), commonly considered as reference for a future high intensity eruption in this area. The magmatic system has been modelled as a long dyke refilling a small shallow magma chamber; magmas with trachytic and phonolitic composition and variable volatile content of H2O and CO2 have been considered. The simulations have been carried out changing the condition of magma injection, the system configuration (magma chamber geometry, dyke size) and the resident and refilling magma composition and volatile content, in order to study the influence of these factors on the simulated dynamics. Simulation results allow to follow each step of the gas-rich magma ascent in the denser magma, highlighting the details of magma convection and mixing. In particular, the presence of more CO2 in the deep magma results in more ecient and faster dynamics. Through this simulations the variation of the gravimetric field has been determined. Afterward, the space-time distribution of stress resulting from numerical simulations have been used as boundary conditions for the simulations of the displacement field imposed by the magmatic dynamics on rocks. The properties of the simulated domain (rock density, P and S wave velocities) have been based on data from literature on active and passive tomographic experiments, obtained through a collaboration with A. Zollo at the Dept. of Physics of the Federici II Univeristy in Napoli. The elasto-dynamics simulations allow to determine the variations of the space-time distribution of deformation and the seismic signal associated with the studied magmatic dynamics. In particular, results show that these dynamics induce deformations similar to those measured at Campi Flegrei and seismic signals with energies concentrated on the typical frequency bands observed in volcanic areas. The present work shows that an approach based on the solution of equations describing the physics of processes within a magmatic fluid and the surrounding rock system is able to recognise and describe the relationships between geophysical signals detectable on the surface and deep magma dynamics. Therefore, the results suggest that the combined study of geophysical data and informations from numerical simulations can allow in a near future a more ecient evaluation of the short term volcanic hazard.

Modeling, design and experimental characterization of Micro-Electro-Mechanical-Systems for gas chromatographic applications

Design parameters, process flows, electro-thermal-fluidic simulations and experimental characterizations of Micro-Electro-Mechanical-Systems (MEMS) suited for gas-chromatographic (GC) applications are presented and thoroughly described in this thesis, whose topic belongs to the research activities the Institute for Microelectronics and Microsystems (IMM)-Bologna is involved since several years, i.e. the development of micro-systems for chemical analysis, based on silicon micro-machining techniques and able to perform analysis of complex gaseous mixtures, especially in the field of environmental monitoring. In this regard, attention has been focused on the development of micro-fabricated devices to be employed in a portable mini-GC system for the analysis of aromatic Volatile Organic Compounds (VOC) like Benzene, Toluene, Ethyl-benzene and Xylene (BTEX), i.e. chemical compounds which can significantly affect environment and human health because of their demonstrated carcinogenicity (benzene) or toxicity (toluene, xylene) even at parts per billion (ppb) concentrations. The most significant results achieved through the laboratory functional characterization of the mini-GC system have been reported, together with in-field analysis results carried out in a station of the Bologna air monitoring network and compared with those provided by a commercial GC system. The development of more advanced prototypes of micro-fabricated devices specifically suited for FAST-GC have been also presented (silicon capillary columns, Ultra-Low-Power (ULP) Metal OXide (MOX) sensor, Thermal Conductivity Detector (TCD)), together with the technological processes for their fabrication. The experimentally demonstrated very high sensitivity of ULP-MOX sensors to VOCs, coupled with the extremely low power consumption, makes the developed ULP-MOX sensor the most performing metal oxide sensor reported up to now in literature, while preliminary test results proved that the developed silicon capillary columns are capable of performances comparable to those of the best fused silica capillary columns. Finally, the development and the validation of a coupled electro-thermal Finite Element Model suited for both steady-state and transient analysis of the micro-devices has been described, and subsequently implemented with a fluidic part to investigate devices behaviour in presence of a gas flowing with certain volumetric flow rates.

Sviluppo di tecniche per la progettazione delle reti di monitoraggio della qualita' dell'aria

L’obiettivo del lavoro consiste nell’implementare una metodologia operativa volta alla progettazione di reti di monitoraggio e di campagne di misura della qualità dell’aria con l’utilizzo del laboratorio mobile, ottimizzando le posizioni dei dispositivi di campionamento rispetto a differenti obiettivi e criteri di scelta. La revisione e l’analisi degli approcci e delle indicazioni fornite dalla normativa di riferimento e dai diversi autori di lavori scientifici ha permesso di proporre un approccio metodologico costituito da due fasi operative principali, che è stato applicato ad un caso studio rappresentato dal territorio della provincia di Ravenna. La metodologia implementata prevede l’integrazione di numerosi strumenti di supporto alla valutazione dello stato di qualità dell’aria e degli effetti che gli inquinanti atmosferici possono generare su specifici recettori sensibili (popolazione residente, vegetazione, beni materiali). In particolare, la metodologia integra approcci di disaggregazione degli inventari delle emissioni attraverso l’utilizzo di variabili proxy, strumenti modellistici per la simulazione della dispersione degli inquinanti in atmosfera ed algoritmi di allocazione degli strumenti di monitoraggio attraverso la massimizzazione (o minimizzazione) di specifiche funzioni obiettivo. La procedura di allocazione sviluppata è stata automatizzata attraverso lo sviluppo di un software che, mediante un’interfaccia grafica di interrogazione, consente di identificare delle aree ottimali per realizzare le diverse campagne di monitoraggio

INTEGRATED ASSESSMENT OF ANTHROPOGENIC, CLIMATE, AND POLICY INDUCED CHANGES OF PHOSPHORUS EXPORT IN THE UNITED STATES LAURENTIAN GREAT LAKES WATERSHEDS

Anthropogenic activities have increased phosphorus (P) loading in tributaries to the Laurentian Great Lakes resulting in eutrophication in small bays to most notably, Lake Erie. Changes to surface water quality from P loading have resulted in billions of dollars in damage and threaten the health of the world’s largest freshwater resource. To understand the factors affecting P delivery with projected increasing urban lands and biofuels expansion, two spatially explicit models were coupled. The coupled models predict that the majority of the basin will experience a significant increase in urban area P sources while the agriculture intensity and forest sources of P will decrease. Changes in P loading across the basin will be highly variable spatially. Additionally, the impacts of climate change on high precipitation events across the Great Lakes were examined. Using historical regression relationships on phosphorus concentrations, key Great Lakes tributaries were found to have future changes including decreasing total loads and increases to high-flow loading events. The urbanized Cuyahoga watersheds exhibits the most vulnerability to these climate-induced changes with increases in total loading and storm loading , while the forested Au Sable watershed exhibits greater resilience. Finally, the monitoring network currently in place for sampling the amount of phosphorus entering the U.S. Great Lakes was examined with a focus on the challenges to monitoring. Based on these interviews, the research identified three issues that policy makers interested in maintaining an effective phosphorus monitoring network in the Great Lakes should consider: first, that the policy objectives driving different monitoring programs vary, which results in different patterns of sampling design and frequency; second, that these differences complicate efforts to encourage collaboration; and third, that methods of funding sampling programs vary from agency to agency, further complicating efforts to generate sufficient long-term data to improve our understanding of phosphorus into the Great Lakes. The dissertation combines these three areas of research to present the potential future impacts of P loading in the Great Lakes as anthropogenic activities, climate and monitoring changes. These manuscripts report new experimental data for future sources, loading and climate impacts on phosphorus.

Table 2 Overview of the established BIOTA Observatories in Africa & Appendix

The international, interdisciplinary biodiversity research project BIOTA AFRICA initiated a standardized biodiversity monitoring network along climatic gradients across the African continent. Due to an identified lack of adequate monitoring designs, BIOTA AFRICA developed and implemented the standardized BIOTA Biodiversity Observatories, that meet the following criteria (a) enable long-term monitoring of biodiversity, potential driving factors, and relevant indicators with adequate spatial and temporal resolution, (b) facilitate comparability of data generated within different ecosystems, (c) allow integration of many disciplines, (d) allow spatial up-scaling, and (e) be applicable within a network approach. A BIOTA Observatory encompasses an area of 1 km2 and is subdivided into 100 1-ha plots. For meeting the needs of sampling of different organism groups, the hectare plot is again subdivided into standardized subplots, whose sizes follow a geometric series. To allow for different sampling intensities but at the same time to characterize the whole square kilometer, the number of hectare plots to be sampled depends on the requirements of the respective discipline. A hierarchical ranking of the hectare plots ensures that all disciplines monitor as many hectare plots jointly as possible. The BIOTA Observatory design assures repeated, multidisciplinary standardized inventories of biodiversity and its environmental drivers, including options for spatial up- and downscaling and different sampling intensities. BIOTA Observatories have been installed along climatic and landscape gradients in Morocco, West Africa, and southern Africa. In regions with varying land use, several BIOTA Observatories are situated close to each other to analyze management effects.