Exploração e análise de dados coletados pelo sistema integrado de monitoramento ambiental utilizando técnicas de visual analytics

Pós-graduação em Ciências Cartográficas - FCT

Genome-Wide Analysis in Brazilian Xavante Indians Reveals Low Degree of Admixture

Characterization of population genetic variation and structure can be used as tools for research in human genetics and population isolates are of great interest. The aim of the present study was to characterize the genetic structure of Xavante Indians and compare it with other populations. The Xavante, an indigenous population living in Brazilian Central Plateau, is one of the largest native groups in Brazil. A subset of 53 unrelated subjects was selected from the initial sample of 300 Xavante Indians. Using 86,197 markers, Xavante were compared with all populations of HapMap Phase III and HGDP-CEPH projects and with a Southeast Brazilian population sample to establish its population structure. Principal Components Analysis showed that the Xavante Indians are concentrated in the Amerindian axis near other populations of known Amerindian ancestry such as Karitiana, Pima, Surui and Maya and a low degree of genetic admixture was observed. This is consistent with the historical records of bottlenecks experience and cultural isolation. By calculating pair-wise F-st statistics we characterized the genetic differentiation between Xavante Indians and representative populations of the HapMap and from HGDP-CEPH project. We found that the genetic differentiation between Xavante Indians and populations of Ameridian, Asian, European, and African ancestry increased progressively. Our results indicate that the Xavante is a population that remained genetically isolated over the past decades and can offer advantages for genome-wide mapping studies of inherited disorders.

Energía Geotérmica, ¿la gran desconocida? : Estudio prospectivo del potencial geotérmico de la isla de Gran Canaria

[ES]En la actualidad, las Islas Canarias se abastecen energéticamente de combustibles fósiles aproximadamente en un 90%. La Energía Geotérmica es un recurso energético renovable con muchas posibilidades de explotación en el Archipiélago Canario para la producción de energía; así lo demuestra su naturaleza volcánica, registros históricos y actividad sismotectónica. Se analizan las características necesarias desde varios puntos de vista para evaluar y cuantificar el potencial geotérmico de la isla de Gran Canaria

Mapping the increasing risk of human alveolar echinococcosis in Limburg, The Netherlands

The parasite Echinococcus multilocularis was first detected in The Netherlands in 1996 and repeated studies have shown that the parasite subsequently spread in the local population of foxes in the province of Limburg. It was not possible to quantify the human risk of alveolar echinococcosis because no relationship between the amount of parasite eggs in the environment and the probability of infection in humans was known. Here, we used the spread of the parasite in The Netherlands as a predictor, together with recently published historical records of the epidemiology of alveolar echinococcosis in Switzerland, to achieve a relative quantification of the risk. Based on these analyses, the human risk in Limburg was simulated and up to three human cases are predicted by 2018. We conclude that the epidemiology of alveolar echinococcosis in The Netherlands might have changed from a period of negligible risk in the past to a period of increasing risk in the forthcoming years.

Winning coal at 78° North : mining, contingency and the Chaîne Opératoire in old Longyear City

The purpose of this thesis is to analyze the evolution of an early 20th century mining system in Spitsbergen as applied by Boston-based Arctic Coal Company (ACC). This analysis will address the following questions: Did the system evolve in a linear, technological-based fashion? Or was the progression more a product of interactions and negotiations with the natural and human landscapes present during the time of occupation? Answers to these questions will be sought through review of historical records and material residues identified during the 2008 field examination on Spitsbergen. The Arctic Coal Company’s flagship mine, ACC Mine No. 1, will serve as the focus for this analysis. The mine was the company’s largest undertaking during its occupation of Longyear Valley and today exhibits a large collection of related features and artifacts. The study will emphasize on the material record within an analysis of technical, environmental and social influences that guided the course of the mining system. The intent of this thesis is a better understanding of how a particular resource extraction industry took root in the Arctic.

Impacts of changing climate and land use on vegetation dynamics in a Mediterranean ecosystem: insights from paleoecology and dynamic modeling

Forests near the Mediterranean coast have been shaped by millennia of human disturbance. Consequently, ecological studies relying on modern observations or historical records may have difficulty assessing natural vegetation dynamics under current and future climate. We combined a sedimentary pollen record from Lago di Massacciucoli, Tuscany, Italy with simulations from the LandClim dynamic vegetation model to determine what vegetation preceded intense human disturbance, how past changes in vegetation relate to fire and browsing, and the potential of an extinct vegetation type under present climate. We simulated vegetation dynamics near Lago di Massaciucoli for the last 7,000 years using a local chironomid-inferred temperature reconstruction with combinations of three fire regimes (small infrequent, large infrequent, small frequent) and three browsing intensities (no browsing, light browsing, and moderate browsing), and compared model output to pollen data. Simulations with low disturbance support pollen-inferred evidence for a mixed forest dominated by Quercus ilex (a Mediterranean species) and Abies alba (a montane species). Whereas pollen data record the collapse of A. alba after 6000 cal yr bp, simulated populations expanded with declining summer temperatures during the late Holocene. Simulations with increased fire and browsing are consistent with evidence for expansion by deciduous species after A. alba collapsed. According to our combined paleo-environmental and modeling evidence, mixed Q. ilex and A. alba forests remain possible with current climate and limited disturbance, and provide a viable management objective for ecosystems near the Mediterranean coast and in regions that are expected to experience a mediterranean-type climate in the future.

Long-term perspective on wildfires in the western USA

Understanding the causes and consequences of wildfires in forests of the western United States requires integrated information about fire, climate changes, and human activity on multiple temporal scales. We use sedimentary charcoal accumulation rates to construct long-term variations in fire during the past 3,000 y in the American West and compare this record to independent fire-history data from historical records and fire scars. There has been a slight decline in burning over the past 3,000 y, with the lowest levels attained during the 20th century and during the Little Ice Age (LIA, ca. 1400-1700 CE Common Era]). Prominent peaks in forest fires occurred during the Medieval Climate Anomaly (ca. 950-1250 CE) and during the 1800s. Analysis of climate reconstructions beginning from 500 CE and population data show that temperature and drought predict changes in biomass burning up to the late 1800s CE. Since the late 1800s, human activities and the ecological effects of recent high fire activity caused a large, abrupt decline in burning similar to the LIA fire decline. Consequently, there is now a forest ``fire deficit'' in the western United States attributable to the combined effects of human activities, ecological, and climate changes. Large fires in the late 20th and 21st century fires have begun to address the fire deficit, but it is continuing to grow.

Reconstructing 4000 years of mass movement and tsunami history in a deep peri-Alpine lake (Lake Geneva, France-Switzerland)

The study of mass movements in lake sediments provides insights into past natural hazards at historic and prehistoric timescales. Sediments from the deep basin of Lake Geneva reveal a succession of six large-scale (volumes of 22 × 106 to 250 × 106 m3) mass-transport deposits, associated with five mass-movement events within 2600 years (4000 cal bp to 563 ad). The mass-transport deposits result from: (i) lateral slope failures (mass-transport deposit B at 3895 ± 225 cal bp and mass-transport deposits A and C at 3683 ± 128 cal bp); and (ii) Rhône delta collapses (mass-transport deposits D to G dated at 2650 ± 150 cal bp, 2185 ± 85 cal bp, 1920 ± 120 cal bp and 563 ad, respectively). Mass-transport deposits A and C were most probably triggered by an earthquake, whereas the Rhône delta collapses were likely to be due to sediment overload with a rockfall as the external trigger (mass-transport deposit G, the Tauredunum event in 563 ad known from historical records), an earthquake (mass-transport deposit E) or unknown external triggers (mass-transport deposits D and F). Independent of their origin and trigger mechanisms, numerical simulations show that all of these recorded mass-transport deposits are large enough to have generated at least metre-scale tsunamis during mass movement initiation. Since the Tauredunum event in 563 ad, two small-scale (volumes of 1 to 2 × 106 m3) mass-transport deposits (H and I) are present in the seismic record, both of which are associated with small lateral slope failures. Mass-transport deposits H and I might be related to earthquakes in Lausanne/Geneva (possibly) 1322 ad and Aigle 1584 ad, respectively. The sedimentary record of the deep basin of Lake Geneva, in combination with the historical record, show that during the past 3695 years, at least six tsunamis were generated by mass movements, indicating that the tsunami hazard in the Lake Geneva region should not be neglected, although such events are not frequent with a recurrence time of 0·0016 yr−1.

Extreme River Floods in Western Switzerland and the Lake of Constance Region in the Period Prior to Instrumental Measurements

The floods that occurred on the Aare and Rhine rivers in May 2015 and the mostly successful handling of this event in terms of flood protection measures are a good reminder of how important it is to comprehend the causes and processes involved in such natural hazards. While the needed data series of gauge measurements and peak discharge calculations reach back to the 19th century, historical records dating further back in time can provide additional and useful information to help understanding extreme flood events and to evaluate prevention measures such as river dams and corrections undertaken prior to instrumental measurements. In my PhD project I will use a wide range of historical sources to assess and quantify past extreme flood events. It is part of the SNF-funded project “Reconstruction of the Genesis, Process and Impact of Major Pre-instrumental Flood Events of Major Swiss Rivers Including a Peak Discharge Quantification” and will cover the research locations Fribourg (Saane R.), Burgdorf (Emme R.), Thun, Bern (both Aare R.), and the Lake of Constance at the locations Lindau, Constance and Rorschach. My main goals are to provide a long time series of quantitative data for extreme flood events, to discuss the occurring changes in these data, and to evaluate the impact of the aforementioned human influences on the drainage system. Extracting information given in account books from the towns of Basel and Solothurn may also enable me to assess the frequency and seasonality of less severe river floods. Finally, historical information will be used for remodeling the historical hydrological regime to homogenize the historical data series to modern day conditions and thus make it comparable to the data provided by instrumental measurements. The method I will apply for processing all information provided by historical sources such as chronicles, newspapers, institutional records, as well as flood marks, paintings and archeological evidence has been developed and successfully applied to the site of Basel by Wetter et al. (2011). They have also shown that data homogenization is possible by reconstructing previous stream flow conditions using historical river profiles and by carefully observing and re-constructing human changes of the river bed and its surroundings. Taken all information into account, peak discharges for past extreme flood events will be calculated with a one-dimensional hydrological model.

Aedes aegypti, Aedes albopictus, and Dengue virus in Harris county: An estimate of risk

Recent outbreaks of dengue fever (DF) along the United States/Mexico border, coupled with the high number of reported cases in Mexico suggest that there is the possibility for DF emergence in Houston, Texas1,2. To determine the presence of DF, populations of Aedes aegypti and Aedes albopictus were identified and tested for dengue virus. Maps were created to identify "hot spots" (Figure 1) based on historical data on Ae. aegypti and Ae. albopictus, demographic information, and locations of human cases of dengue fever. BG Sentinel Traps®, in conjunction with BG Lure® attractant, octanol and dry ice, were used to collect mosquitoes, which were then tested for presence of dengue virus using ELISA techniques. All samples tested were negative for dengue virus (DV). Survival of DV ultimately comes down to whether or not it will be vectored by a mosquito to a susceptible human host. The presence of infected humans and contact with the mosquito vectors are two critical factors necessary in the establishment of DF. Historical records indicate the presence of Ae. aegypti and Ae. albopictus in Harris County, which would support localized dengue transmission if infected individuals are present.^ (1) Brunkard JM, Robles-Lopez JL, Ramirez J, Cifuentes E, Rothenberg SJ, Hunsperger EA, Moore CG, Brussolo RM, Villarreal NA, Haddad BM, 2007. Dengue fever seroprevalence and risk factors, Texas-Mexico border, 2004. Emerg Infect Dis 13: 1477-1483. (2) Ramos MM, Mohammed H, Zielinski-Gutierrez E, Hayden MH, Lopez JL, Fournier M, Trujillo AR, Burton R, Brunkard JM, Anaya-Lopez L, Banicki AA, Morales PK, Smith B, Munoz JL, Waterman SH, 2008. Epidemic dengue and dengue hemorrhagic fever at the Texas-Mexico Border: results of a household-based seroepidemiologic survey, December 2005. Am J Trop Med Hyg 78: 364-369.^

Impact of environmental regulations on lead concentrations in Texas ambient air and water: 1970--1990

In the last two decades, the significance of lead has been addressed in a number of environmental regulations at the national and state levels. This project investigated the environmental regulations (Clean Air Act and Amendments, 1970-1990 and Clean Water Act of 1977) and their cumulative effects on lead in ambient air and water in the state of Texas. For this purpose, historical records from the Texas Water Development Board, Texas Natural Resources Conservation Commission, and the United States Geological Survey have been assembled and analyzed for temporal and spatial trends. These trends might correspond to the phase out of lead in gasoline and other regulations.^ This study concluded that there is a significant correlation (p $\leq$.001) between environmental regulations of lead in gasoline and the concentration of lead in ambient air. Lead concentrations in ambient air have been reduced by over 90 percent in the past twenty years. An overall significant difference (p $\leq$.001) was found in mean (94, 15 respectively) lead concentrations in surface water between two time periods, one at the beginning of the twenty year period and one at the end of the study period. There has been an overall reduction of lead concentrations in surface water in Texas of approximately 84 percent. However, this reduction cannot be statistically associated with any one regulation. Groundwater data could not be analyzed for lead concentrations because of limitations of reporting data as "less than". Approximately two percent of the groundwater data was analyzed by Oneway ANOVA and no significant difference was found between the means (18, 19 respectively) of two time periods, 1977-1979 and 1988-1990. This data is consistent with the regulations having a contributory affect on declining concentrations, but other factors cannot be ruled out as having added to these declines. This study can also serve as a starting point for a more in-depth study of environmental regulations and their impact on the environment. ^

Timing of the most recent Neoglacial advance and retreat in Potter Cove, Antarctic Peninsula

The timing of the most recent Neoglacial advance in the Antarctic Peninsula is important for establishing global climate teleconnections and providing important post-glacial rebound corrections to gravity-based satellite measurements of ice loss. However, obtaining accurate ages from terrestrial geomorphic and sedimentary indicators of the most recent Neoglacial advance in Antarctica has been hampered by the lack of historical records and the difficulty of dating materials in Antarctica. Here we use a new approach to dating flights of raised beaches in the South Shetland Islands of the northern Antarctic Peninsula to bracket the age of a Neoglacial advance that occurred between 1500 and 1700 AD, broadly synchronous with compilations for the timing of the Little Ice Age in the northern hemisphere. Our approach is based on optically stimulated luminescence of the underside of buried cobbles to obtain the age of beaches previously shown to have been deposited immediately inside and outside the moraines of the most recent Neoglacial advance. In addition, these beaches mark the timing of an apparent change in the rate of isostatic rebound thought to be in response to the same glacial advance within the South Shetland Islands. We use a Maxwell viscoelastic model of glacial-isostatic adjustment (GIA) to determine whether the rates of uplift calculated from the raised beaches are realistic given the limited constraints on the ice advance during this most recent Neoglacial advance. Our rebound model suggests that the subsequent melting of an additional 16-22% increase in the volume of ice within the South Shetland Islands would result in a subsequent uplift rate of 12.5 mm/yr that lasted until 1840 AD resulting in a cumulative uplift of 2.5 m. This uplift rate and magnitude are in close agreement with observed rates and magnitudes calculated from the raised beaches since the most recent Neoglacial advance along the South Shetland Islands and falls within the range of uplift rates from similar settings such as Alaska.

Modelización de costas rocosas acantiladas

El retroceso de las costas acantiladas es un fenómeno muy extendido sobre los litorales rocosos expuestos a la incidencia combinada de los procesos marinos y meteorológicos que se dan en la franja costera. Este fenómeno se revela violentamente como movimientos gravitacionales del terreno esporádicos, pudiendo causar pérdidas materiales y/o humanas. Aunque el conocimiento de estos riesgos de erosión resulta de vital importancia para la correcta gestión de la costa, el desarrollo de modelos predictivos se encuentra limitado desde el punto de vista geomorfológico debido a la complejidad e interacción de los procesos de desarrollo espacio-temporal que tienen lugar en la zona costera. Los modelos de predicción publicados son escasos y con importantes inconvenientes: a) extrapolación, extienden la información de registros históricos; b) empíricos, sobre registros históricos estudian la respuesta al cambio de un parámetro; c) estocásticos, determinan la cadencia y magnitud de los eventos futuros extrapolando las distribuciones de probabilidad extraídas de catálogos históricos; d) proceso-respuesta, de estabilidad y propagación del error inexplorada; e) en Ecuaciones en Derivadas Parciales, computacionalmente costosos y poco exactos. La primera parte de esta tesis detalla las principales características de los modelos más recientes de cada tipo y, para los más habitualmente utilizados, se indican sus rangos de aplicación, ventajas e inconvenientes. Finalmente como síntesis de los procesos más relevantes que contemplan los modelos revisados, se presenta un diagrama conceptual de la recesión costera, donde se recogen los procesos más influyentes que deben ser tenidos en cuenta, a la hora de utilizar o crear un modelo de recesión costera con el objetivo de evaluar la peligrosidad (tiempo/frecuencia) del fenómeno a medio-corto plazo. En esta tesis se desarrolla un modelo de proceso-respuesta de retroceso de acantilados costeros que incorpora el comportamiento geomecánico de materiales cuya resistencia a compresión no supere los 5 MPa. El modelo simula la evolución espaciotemporal de un perfil-2D del acantilado que puede estar formado por materiales heterogéneos. Para ello, se acoplan la dinámica marina: nivel medio del mar, cambios en el nivel medio del lago, mareas y oleaje; con la evolución del terreno: erosión, desprendimiento rocoso y formación de talud de derrubios. El modelo en sus diferentes variantes es capaz de incluir el análisis de la estabilidad geomecánica de los materiales, el efecto de los derrubios presentes al pie del acantilado, el efecto del agua subterránea, la playa, el run-up, cambios en el nivel medio del mar o cambios (estacionales o interanuales) en el nivel medio de la masa de agua (lagos). Se ha estudiado el error de discretización del modelo y su propagación en el tiempo a partir de las soluciones exactas para los dos primeros periodos de marea para diferentes aproximaciones numéricas tanto en tiempo como en espacio. Los resultados obtenidos han permitido justificar las elecciones que minimizan el error y los métodos de aproximación más adecuados para su posterior uso en la modelización. El modelo ha sido validado frente a datos reales en la costa de Holderness, Yorkshire, Reino Unido; y en la costa norte del lago Erie, Ontario, Canadá. Los resultados obtenidos presentan un importante avance en los modelos de recesión costera, especialmente en su relación con las condiciones geomecánicas del medio, la influencia del agua subterránea, la verticalización de los perfiles rocosos y su respuesta ante condiciones variables producidas por el cambio climático (por ejemplo, nivel medio del mar, cambios en los niveles de lago, etc.). The recession of coastal cliffs is a widespread phenomenon on the rocky shores that are exposed to the combined incidence of marine and meteorological processes that occur in the shoreline. This phenomenon is revealed violently and occasionally, as gravitational movements of the ground and can cause material or human losses. Although knowledge of the risks of erosion is vital for the proper management of the coast, the development of cliff erosion predictive models is limited by the complex interactions between environmental processes and material properties over a range of temporal and spatial scales. Published prediction models are scarce and present important drawbacks: extrapolation, that extend historical records to the future; empirical, that based on historical records studies the system response against the change in one parameter; stochastic, that represent of cliff behaviour based on assumptions regarding the magnitude and frequency of events in a probabilistic framework based on historical records; process-response, stability and error propagation unexplored; PDE´s, highly computationally expensive and not very accurate. The first part of this thesis describes the main features of the latest models of each type and, for the most commonly used, their ranges of application, advantages and disadvantages are given. Finally as a synthesis of the most relevant processes that include the revised models, a conceptual diagram of coastal recession is presented. This conceptual model includes the most influential processes that must be taken into account when using or creating a model of coastal recession to evaluate the dangerousness (time/frequency) of the phenomenon to medium-short term. A new process-response coastal recession model developed in this thesis has been designed to incorporate the behavioural and mechanical characteristics of coastal cliffs which are composed of with materials whose compressive strength is less than 5 MPa. The model simulates the spatial and temporal evolution of a cliff-2D profile that can consist of heterogeneous materials. To do so, marine dynamics: mean sea level, waves, tides, lake seasonal changes; is coupled with the evolution of land recession: erosion, cliff face failure and associated protective colluvial wedge. The model in its different variants can include analysis of material geomechanical stability, the effect of debris present at the cliff foot, groundwater effects, beach and run-up effects, changes in the mean sea level or changes (seasonal or inter-annual) in the mean lake level. Computational implementation and study of different numerical resolution techniques, in both time and space approximations, and the produced errors are exposed and analysed for the first two tidal periods. The results obtained in the errors analysis allow us to operate the model with a configuration that minimizes the error of the approximation methods. The model is validated through profile evolution assessment at various locations of coastline retreat on the Holderness Coast, Yorkshire, UK and on the north coast of Lake Erie, Ontario, Canada. The results represent an important stepforward in linking material properties to the processes of cliff recession, in considering the effect of groundwater charge and the slope oversteeping and their response to changing conditions caused by climate change (i.e. sea level, changes in lakes levels, etc.).

Diseño de un modelo de riesgo integral de incendios forestales mediante técnicas multicriterio y su automatización en sistemas de información geográfica : una aplicación en la Comunidad Valenciana

Los incendios forestales son la principal causa de mortalidad de árboles en la Europa mediterránea y constituyen la amenaza más seria para los ecosistemas forestales españoles. En la Comunidad Valenciana, diariamente se despliega cerca de un centenar de vehículos de vigilancia, cuya distribución se apoya, fundamentalmente, en un índice de riesgo de incendios calculado en función de las condiciones meteorológicas. La tesis se centra en el diseño y validación de un nuevo índice de riesgo integrado de incendios, especialmente adaptado a la región mediterránea y que facilite el proceso de toma de decisiones en la distribución diaria de los medios de vigilancia contra incendios forestales. El índice adopta el enfoque de riesgo integrado introducido en la última década y que incluye dos componentes de riesgo: el peligro de ignición y la vulnerabilidad. El primero representa la probabilidad de que se inicie un fuego y el peligro potencial para que se propague, mientras que la vulnerabilidad tiene en cuenta las características del territorio y los efectos potenciales del fuego sobre el mismo. Para el cálculo del peligro potencial se han identificado indicadores relativos a los agentes naturales y humanos causantes de incendios, la ocurrencia histórica y el estado de los combustibles, extremo muy relacionado con la meteorología y las especies. En cuanto a la vulnerabilidad se han empleado indicadores representativos de los efectos potenciales del incendio (comportamiento del fuego, infraestructuras de defensa), como de las características del terreno (valor, capacidad de regeneración…). Todos estos indicadores constituyen una estructura jerárquica en la que, siguiendo las recomendaciones de la Comisión europea para índices de riesgo de incendios, se han incluido indicadores representativos del riesgo a corto plazo y a largo plazo. El cálculo del valor final del índice se ha llevado a cabo mediante la progresiva agregación de los componentes que forman cada uno de los niveles de la estructura jerárquica del índice y su integración final. Puesto que las técnicas de decisión multicriterio están especialmente orientadas a tratar con problemas basados en estructuras jerárquicas, se ha aplicado el método TOPSIS para obtener la integración final del modelo. Se ha introducido en el modelo la opinión de los expertos, mediante la ponderación de cada uno de los componentes del índice. Se ha utilizado el método AHP, para obtener las ponderaciones de cada experto y su integración en un único peso por cada indicador. Para la validación del índice se han empleado los modelos de Ecuaciones de Estimación Generalizadas, que tienen en cuenta posibles respuestas correlacionadas. Para llevarla a cabo se emplearon los datos de oficiales de incendios ocurridos durante el período 1994 al 2003, referenciados a una cuadrícula de 10x10 km empleando la ocurrencia de incendios y su superficie, como variables dependientes. Los resultados de la validación muestran un buen funcionamiento del subíndice de peligro de ocurrencia con un alto grado de correlación entre el subíndice y la ocurrencia, un buen ajuste del modelo logístico y un buen poder discriminante. Por su parte, el subíndice de vulnerabilidad no ha presentado una correlación significativa entre sus valores y la superficie de los incendios, lo que no descarta su validez, ya que algunos de sus componentes tienen un carácter subjetivo, independiente de la superficie incendiada. En general el índice presenta un buen funcionamiento para la distribución de los medios de vigilancia en función del peligro de inicio. No obstante, se identifican y discuten nuevas líneas de investigación que podrían conducir a una mejora del ajuste global del índice. En concreto se plantea la necesidad de estudiar más profundamente la aparente correlación que existe en la provincia de Valencia entre la superficie forestal que ocupa cada cuadrícula de 10 km del territorio y su riesgo de incendios y que parece que a menor superficie forestal, mayor riesgo de incendio. Otros aspectos a investigar son la sensibilidad de los pesos de cada componente o la introducción de factores relativos a los medios potenciales de extinción en el subíndice de vulnerabilidad. Summary Forest fires are the main cause of tree mortality in Mediterranean Europe and the most serious threat to the Spanisf forest. In the Spanish autonomous region of Valencia, forest administration deploys a mobile fleet of 100 surveillance vehicles in forest land whose allocation is based on meteorological index of wildlandfire risk. This thesis is focused on the design and validation of a new Integrated Wildland Fire Risk Index proposed to efficient allocation of vehicles and specially adapted to the Mediterranean conditions. Following the approaches of integrated risk developed last decade, the index includes two risk components: Wildland Fire Danger and Vulnerability. The former represents the probability a fire ignites and the potential hazard of fire propagation or spread danger, while vulnerability accounts for characteristics of the land and potential effects of fire. To calculate the Wildland Fire Danger, indicators of ignition and spread danger have been identified, including human and natural occurrence agents, fuel conditions, historical occurrence and spread rate. Regarding vulnerability se han empleado indicadores representativos de los efectos potenciales del incendio (comportamiento del fuego, infraestructurasd de defensa), como de las características del terreno (valor, capacidad de regeneración…). These indicators make up the hierarchical structure for the index, which, following the criteria of the European Commission both short and long-term indicators have been included. Integration consists of the progressive aggregation of the components that make up every level in risk the index and, after that, the integration of these levels to obtain a unique value for the index. As Munticriteria methods are oriented to deal with hierarchically structured problems and with situations in which conflicting goals prevail, TOPSIS method is used in the integration of components. Multicriteria methods were also used to incorporate expert opinion in weighting of indicators and to carry out the aggregation process into the final index. The Analytic Hierarchy Process method was used to aggregate experts' opinions on each component into a single value. Generalized Estimation Equations, which account for possible correlated responses, were used to validate the index. Historical records of daily occurrence for the period from 1994 to 2003, referred to a 10x10-km-grid cell, as well as the extent of the fires were the dependant variables. The results of validation showed good Wildland Fire Danger component performance, with high correlation degree between Danger and occurrence, a good fit of the logistic model used and a good discrimination power. The vulnerability component has not showed a significant correlation between their values and surface fires, which does not mean the index is not valid, because of the subjective character of some of its components, independent of the surface of the fires. Overall, the index could be used to optimize the preventing resources allocation. Nevertheless, new researching lines are identified and discussed to improve the overall performance of the index. More specifically the need of study the inverse relationship between the value of the wildfire Fire Danger component and the forested surface of each 10 - km cell is set out. Other points to be researched are the sensitivity of the index component´s weight and the possibility of taking into account indicators related to fire fighting resources to make up the vulnerability component.

Estimación de la peligrosidad sísmica mediante protocolos de juicios de expertos

Este trabajo estudia la aportación que los métodos de agregación de juicios de expertos pueden realizar en el cálculo de la peligrosidad sísmica de emplazamientos. Se han realizado cálculos en dos emplazamientos de la Península Ibérica: Mugardos (La Coruña) y Cofrentes (Valencia) que están sometidos a regímenes tectónicos distintos y que, además, alojan instalaciones industriales de gran responsabilidad. Las zonas de estudio, de 320 Km de radio, son independientes. Se ha aplicado un planteamiento probabilista a la estimación de la tasa anual de superación de valores de la aceleración horizontal de pico y se ha utilizado el Método de Montecarlo para incorporar a los resultados la incertidumbre presente en los datos relativos a la definición de cada fuente sismogenética y de su sismicidad. Los cálculos se han operado mediante un programa de ordenador, desarrollado para este trabajo, que utiliza la metodología propuesta por el Senior Seismic Hazard Analysis Commitee (1997) para la NRC. La primera conclusión de los resultados ha sido que la Atenuación es la fuente principal de incertidumbre en las estimaciones de peligrosidad en ambos casos. Dada la dificultad de completar los datos históricos disponibles de esta variable se ha estudiado el comportamiento de cuatro métodos matemáticos de agregación de juicios de expertos a la hora de estimar una ley de atenuación en un emplazamiento. Los datos de partida se han obtenido del Catálogo de Isosistas del IGN. Los sismos utilizados como variables raíz se han elegido con el criterio de cubrir uniformemente la serie histórica disponible y los valores de magnitud observados. Se ha asignado un panel de expertos particular a cada uno de los dos emplazamientos y se han aplicado a sus juicios los métodos de Cooke, equipesos, Apostolakis_Mosleh y Morris. Sus propuestas se han comparado con los datos reales para juzgar su eficacia y su facilidad de operación. A partir de los resultados se ha concluido que el método de Cooke ha mostrado el comportamiento más eficiente y robusto para ambos emplazamientos. Este método, además, ha permitido identificar, razonadamente, a aquellos expertos que no deberían haberse introducido en un panel. The present work analyses the possible contribution of the mathematical methods of aggregation in the assessment of Seismic Hazzard. Two sites, in the Iberian Peninsula, have been considered: Mugardos ( La Coruña) and Cofrentes (Valencia).Both of them are subjected to different tectonic regimes an both accommodate high value industrial plants. Their areas of concern, with radius of 320 Km, are not overlapping. A probabilistic approach has been applied in the assessment the annual probability of exceedence of the horizontal peak acceleration. The Montecarlo Method has allowed to transfer the uncertainty in the models and parameters to the final results. A computer program has been developed for this purpose. The methodology proposed by the Senior Seismic Analysis Committee (1997) for the NRC has been considered. Attenuation in Ground motion has been proved to be the main source of uncertainty in seismic hazard for both sites. Taking into account the difficulties to complete existing historical data in this subject the performance of four mathematical methods of aggregation has been studied. Original data have been obtained from the catalogs of the Spanish National Institute of Geography. The seismic events considered were chosen to cover evenly the historical records and the observed values of magnitude. A panel of experts have been applied to each site and four aggregation methods have been developed : equal weights, Cooke, Apostolakis-Mosleh and Morris The four proposals have been compaired with the actual data to judge their performance and ease of application. The results have shown that the Method of Cooke have proved the most efficient and robust for both sites. This method, besides, allow the reasoned identification of those experts who should be rejected from the panel