Metodología para la estimación de indicadores armonizados a partir de los inventarios forestales nacionales europeos con especial énfasis en la biodiversidad forestal

Durante las últimas décadas el objetivo principal de la silvicultura y la gestión forestal en Europa ha pasado de ser la producción de madera a ser la gestión sostenible de los ecosistemas, por lo que se deben considerar todos los bienes y servicios que proporcionan los bosques. En consecuencia, es necesario contar con información forestal periódica de diversos indicadores forestales a nivel europeo para apoyar el desarrollo y la implementación de políticas medioambientales y que se realice una gestión adecuada. Para ello, se requiere un seguimiento intensivo sobre el estado de los bosques, por lo que los Inventarios Forestales Nacionales (IFN), (principal fuente de información forestal a gran escala), han aumentado el número de variables muestreadas para cumplir con los crecientes requerimientos de información. Sin embargo, las estimaciones proporcionadas por los diferentes países no son fácilmente comparables debido a las diferencias en las definiciones, los diseños de muestreo, las variables medidas y los protocolos de medición. Por esto, la armonización de los datos que proporcionan los diferentes países es fundamental para la contar con una información forestal sólida y fiable en la Unión europea (UE). La presente tesis tiene dos objetivos principales: (i) establecer el diseño de una metodología para evaluar la biodiversidad forestal en el marco del Inventario forestal nacional de España teniendo en cuenta las diferentes iniciativas nacionales e internacionales, con el objetivo de producir estimaciones comparables con las de otros países de la UE y (ii) armonizar los indicadores más relevantes para satisfacer los requerimientos nacionales e internacionales. Como consecuencia del estudio realizado para alcanzar el primer objetivo, la metodología diseñada para estimar la biodiversidad fue adoptada por el Tercer Inventario forestal nacional. Ésta se componía de indicadores agrupados en: cobertura del suelo, composición de árboles y especies de arbustos, riqueza de especies herbáceas y helechos, especies amenazadas, estructura, madera muerta, y líquenes epífitos. Tras el análisis del diseño metodológico y de los datos proporcionados, se observó la conveniencia de modificarla con el fin de optimizar los costes, viabilidad, calidad y cantidad de los datos registrados. En consecuencia, en el Cuarto Inventario Forestal Nacional se aplica una metodología modificada, puesto que se eliminó el muestreo de especies herbáceas y helechos, de líquenes epífitos y de especies amenazadas, se modificaron los protocolos de la toma de datos de estructura y madera muerta y se añadió el muestreo de especies invasoras, edad, ramoneo y grado de naturalidad de la masa. En lo que se refiere al segundo objetivo, se ha avanzado en la armonización de tres grupos de variables considerados como relevantes en el marco de los IFN: los indicadores de vegetación no arbórea (que juegan un papel relevante en los ecosistemas, es donde existe la mayor diversidad de plantas y hasta ahora no se conocían los datos muestreados en los IFN), la determinación de los árboles añosos (que tienen un importante papel como nicho ecológico y su identificación es especialmente relevante para la evaluación de la biodiversidad forestal) y el bosque disponible para el suministro de madera (indicador básico de los requerimientos internacionales de información forestal). Se llevó a cabo un estudio completo de la posible armonización de los indicadores de la vegetación no arbórea en los IFN. Para ello, se identificaron y analizaron las diferentes definiciones y diseños de muestreo empleados por los IFN, se establecieron definiciones de referencia y se propusieron y analizaron dos indicadores que pudiesen ser armonizados: MSC (mean species cover) que corresponde a la media de la fracción de cabida cubierta de cada especie por tipo de bosque y MTC (mean total cover). Se estableció una nueva metodología que permite identificar los árboles añosos con los datos proporcionados por los inventarios forestales nacionales con el objetivo de proporcionar una herramienta eficaz para facilitar la gestión forestal considerando la diversidad de los sistemas forestales. Se analizó el concepto de "bosque disponible para el suministro de madera" (FAWS) estudiando la consistencia de la información internacional disponible con el fin de armonizar su estimación y de proporcionar recomendaciones para satisfacer los requerimientos europeos. Como resultado, se elaboró una nueva definición de referencia de FAWS (que será adoptada por el proceso paneuropeo) y se analiza el impacto de la adopción de esta nueva definición en siete países europeos. El trabajo realizado en esta tesis, puede facilitar el suministrar y/o armonizar parcial o totalmente casi la mitad de los indicadores de información forestal solicitados por los requerimientos internacionales (47%). De éstos, prácticamente un 85% tienen relación con los datos inventariados empleando la metodología propuesta para la estimación de la biodiversidad forestal, y el resto, con el establecimiento de la definición de bosque disponible para el suministro de madera. No obstante, y pese a que esta tesis supone un avance importante, queda patente que las necesidades de información forestal son cambiantes y es imprescindible continuar el proceso de armonización de los IFN europeos. ABSTRACT Over the last few decades, the objectives on forestry and forest management in Europe have shifted from being primarily focused on wood production to sustainable ecosystem management, which should consider all the goods and services provided by the forest. Therefore, there is a continued need for forest indicators and assessments at EU level to support the development and implementation of a number of European environmental policies and to conduct a proper forest management. To address these questions, intensive monitoring on the status of forests is required. Therefore, the scope of National Forest Inventories (NFIs), (primary source of data for national and large-area assessments), has been broadened to include new variables to meet these increasing information requirements. However, estimates produced by different countries are not easily comparable because of differences in NFI definitions, plot configurations, measured variables, and measurement protocols. As consequence, harmonizing data produced at national level is essential for the production of sound EU forest information. The present thesis has two main aims: (i) to establish a methodology design to assess forest biodiversity in the frame of the Spanish National Forest Inventory taking into account the different national and international initiatives with the intention to produce comparable estimates with other EU countries and (ii) to harmonize relevant indicators for national and international requirements. In consequence of the work done related to the first objective, the established methodology to estimate forest biodiversity was adopted and launched under the Third National Forest Inventory. It was composed of indicators grouped into: cover, woody species composition, richness of herbaceous species and ferns, endangered species, stand structure, dead wood, and epiphytic lichens. This methodology was analyzed considering the provided data, time costs, feasibility, and requirements. Consequently, in the ongoing Fourth National Forest Inventory a modified methodology is applied: sampling of herbaceous species and ferns, epiphytic lichens and endangered species were removed, protocols regarding structure and deadwood were modified, and sampling of invasive species, age, browsing impact and naturalness were added. As regards the second objective, progress has been made in harmonizing three groups of variables considered relevant in the context of IFN: Indicators of non-tree vegetation (which play an important role in forest ecosystems, it is where the highest diversity of plants occur and so far the related sampled data in NFIs were not known), the identification of old-growth trees (which have an important role as ecological niche and its identification is especially relevant for the assessment of forest biodiversity) and the available forest for wood supply (basic indicator of international forestry information requirements). A complete analysis of ground vegetation harmonization possibilities within NFIs frame was carried on by identifying and analyzing the different definitions and sampling techniques used by NFIs, providing reference definitions related to ground vegetation and proposing and analyzing two ground vegetation harmonized indicators: “Mean species cover” (MSC) and “Mean total cover” (MTC) for shrubs by European forest categories. A new methodology based on NFI data was established with the aim to provide an efficient tool for policy makers to estimate the number of old-growth trees and thus to be able to perform the analysis of the effect of forest management on the diversity associated to forest systems. The concept of “forest available for wood supply” (FAWS) was discussed and clarified, analyzing the consistency of the available international information on FAWS in order to provide recommendations for data harmonization at European level regarding National Forest Inventories (NFIs). As a result, a new reference definition of FAWS was provided (which will be adopted in the pan-European process) and the consequences of the use of this new definition in seven European countries are analyzed. The studies carried on in this thesis, can facilitate the supply and/or harmonization partially or fully of almost half of the forest indicators (47%) needed for international requirements. Of these, nearly 85% are related to inventoried data using the proposed methodology for the estimation of forest biodiversity, and the rest, with the establishment of the definition of forest available for wood supply. However, despite this thesis imply an important development, forest information needs are changing and it is imperative to continue the process of harmonization of European NFIs.

Análisis técnico y financiero de las máquinas de elevación

Los nombres que se usan actualmente para las grúas los ponen las empresas fabricantes y muy frecuentemente no tienen relación con su tipología ni con su capacidad. Por otra parte, es de uso común en la construcción, llamar a las grúas usando su tonelaje nominal que coincide normalmente con su capacidad máxima que se obtiene a radio mínimo. Existe una controversia por el uso de este valor ya que no suele definir bien la capacidad de las maquinas. En cuanto el radio de trabajo se aleja de sus valores mínimos, las grúas están limitadas por el momento de vuelco que no tiene porque comportarse de manera proporcional o ni siquiera relacionada con el valor de la capacidad nominal. Esto hace que comparar grúas mediante sus capacidades nominales (que son sus denominaciones) pueda inducir a errores importantes. Como alternativa, se pretende estudiar el uso de momento máximo de vuelco MLM por sus siglas en ingles (Maximum Load Moment) para intentar definir la capacidad real de las grúas. Se procede a realizar un análisis técnico y financiero de grúas con respecto a ambos valores mencionados con objeto de poder determinar cual de los dos parámetros es más fiable a la hora de definir la capacidad real de estas maquinas. Para ello, se seleccionan dentro de las tres tipologías más importantes por su presencia e importancia en la construcción (grúas de celosía sobre cadenas, grúas telescópicas sobre camión y grúas torre) nueve grúas de distintos tamaños y capacidades con objeto de analizar una serie de parámetros técnicos y sus costes. Se realizan de este modo diversas comparativas analizando los resultados en función de las tipologías y de los tamaños de las distintas maquinas. Para cada máquina se obtienen las capacidades y los momentos de vuelco correspondientes a distintos radios de trabajo. Asimismo, se obtiene el MLM y el coste hora de cada grúa, este último como suma de la amortización de la máquina, intereses del capital invertido, consumos, mantenimiento y coste del operador. Los resultados muestran las claras deficiencias del tonelaje nominal como valor de referencia para definir la capacidad de las grúas ya que grúas con el mismo tonelaje nominal pueden dar valores de capacidad de tres a uno (e incluso mayores) cuando los radios de trabajo son importantes. A raiz de este análisis se propone el uso del MLM en lugar del tonelaje nominal para la denominación de las grúas ya que es un parámetro mucho más fiable. Siendo conscientes de la dificultad que supone un cambio de esta entidad al tratarse de un uso común a nivel mundial, se indican posibles actuaciones concretas que puedan ir avanzando en esa dirección como seria por ejemplo la nomenclatura oficial de los fabricantes usando el MLM dentro del nombre de la grúa que también podría incluir la tipología o al menos alguna actuación legislativa sencilla como obligar al fabricante a indicar este valor en las tablas y características de cada máquina. El ratio analizado Coste horario de la grúa / MLM resulta ser de gran interés y permite llegar a la conclusión que en todas las tipologías de grúas, la eficiencia del coste por hora y por la capacidad (dada por el MLM) aumenta al aumentar la capacidad de la grúa. Cuando los tamaños de cada tipología se reducen, esta eficiencia disminuye y en algunos casos incluso drasticamente. La tendencia del mundo de la construcción de prefabricación y modularización que conlleva pesos y dimensiones de cargas cada vez más grandes, demandan cada vez grúas de mayor capacidad y se podría pensar en un primer momento que ante un crecimiento de capacidades tan significativo, el coste de las grúas se podría disparar y por ello disminuir la eficiencia de estas máquinas. A la vista de los resultados obtenidos con este análisis, no solo no ocurre este problema sino que se observa que dicho aumento de tamaños y capacidades de grúas redunda en un aumento de su eficiencia en cualquiera de las tipologías de estas maquinas que han sido estudiadas. The crane names that are actually used are given by crane manufacturers and, very frequently, they do not have any relationship with the crane type nor with its capacity. On the other hand, it is common in construction to use the nominal capacity (which corresponds in general to the capacity at minimum radius) as crane name. The use of this figure is controversial since it does not really reflect the real crane capacity. When the working radius increases a certain amount from its minimum values, the crane capacity starts to be limited by the crane tipping load and the moment is not really related to the value of the nominal capacity. Therefore, comparing cranes by their nominal capacity (their names) can lead to important mistakes. As an alternative, the use of the maximum load moment (MLM) can be studied for a better definition of real crane capacity. A technical and financial analysis of cranes is conducted using both parameters to determine which one is more reliable in order to define crane’s real capacity. For this purpose, nine cranes with different sizes and capacities will be selected within the most relevant crane types (crawler lattice boom cranes, telescopic truck mounted cranes and tower cranes) in order to analyze several parameters. The technical and economic results will be compared according to the crane types and sizes of the machines. For each machine, capacities and load moments are obtained for several working radius as well as MLM and hourly costs of cranes. Hourly cost is calculated adding up depreciation, interests of invested capital, consumables, maintenance and operator’s cost. The results show clear limitations for the use of nominal capacity as a reference value for crane definition since cranes with the same nominal capacity can have capacity differences of 3 to 1 (or even bigger) when working on important radius. From this analysis, the use of MLM as crane name instead of nominal capacity is proposed since it is a much more reliable figure. Being aware of the difficulty of such change since nominal capacity is commonly used worldwide; specific actions are suggested to progress in that direction. One good example would be that manufacturers would include MLM in their official crane names which could also include the type as well. Even legal action can be taken by simply requiring to state this figure in the crane charts and characteristics of every machine. The analyzed ratio: hourly cost / MLM is really interesting since it leads to the conclusion that for all crane types, the efficiency of the hourly cost divided by capacity (given by MLM) increases when the crane capacity is higher. When crane sizes are smaller, this efficiency is lower and can fall dramatically in certain cases. The developments in the construction world regarding prefabrication and modularization mean bigger weights and dimensions, which create a demand for bigger crane capacities. On a first approach, it could be thought that crane costs could rise significantly because of this capacity hugh increase reducing in this way crane efficiency. From the results obtained here, it is clear that it is definitely not the case but the capacity increase of cranes will end up in higher efficiency levels for all crane types that have been studied.

Contribución a la localización en interiores basada en comunicaciones inalámbricas

Determinar con buena precisión la posición en la que se encuentra un terminal móvil, cuando éste se halla inmerso en un entorno de interior (centros comerciales, edificios de oficinas, aeropuertos, estaciones, túneles, etc), es el pilar básico sobre el que se sustentan un gran número de aplicaciones y servicios. Muchos de esos servicios se encuentran ya disponibles en entornos de exterior, aunque los entornos de interior se prestan a otros servicios específicos para ellos. Ese número, sin embargo, podría ser significativamente mayor de lo que actualmente es, si no fuera necesaria una costosa infraestructura para llevar a cabo el posicionamiento con la precisión adecuada a cada uno de los hipotéticos servicios. O, igualmente, si la citada infraestructura pudiera tener otros usos distintos, además del relacionado con el posicionamiento. La usabilidad de la misma infraestructura para otros fines distintos ofrecería la oportunidad de que la misma estuviera ya presente en las diferentes localizaciones, porque ha sido previamente desplegada para esos otros usos; o bien facilitaría su despliegue, porque el coste de esa operación ofreciera un mayor retorno de usabilidad para quien lo realiza. Las tecnologías inalámbricas de comunicaciones basadas en radiofrecuencia, ya en uso para las comunicaciones de voz y datos (móviles, WLAN, etc), cumplen el requisito anteriormente indicado y, por tanto, facilitarían el crecimiento de las aplicaciones y servicios basados en el posicionamiento, en el caso de poderse emplear para ello. Sin embargo, determinar la posición con el nivel de precisión adecuado mediante el uso de estas tecnologías, es un importante reto hoy en día. El presente trabajo pretende aportar avances significativos en este campo. A lo largo del mismo se llevará a cabo, en primer lugar, un estudio de los principales algoritmos y técnicas auxiliares de posicionamiento aplicables en entornos de interior. La revisión se centrará en aquellos que sean aptos tanto para tecnologías móviles de última generación como para entornos WLAN. Con ello, se pretende poner de relieve las ventajas e inconvenientes de cada uno de estos algoritmos, teniendo como motivación final su aplicabilidad tanto al mundo de las redes móviles 3G y 4G (en especial a las femtoceldas y small-cells LTE) como al indicado entorno WLAN; y teniendo siempre presente que el objetivo último es que vayan a ser usados en interiores. La principal conclusión de esa revisión es que las técnicas de triangulación, comúnmente empleadas para realizar la localización en entornos de exterior, se muestran inútiles en los entornos de interior, debido a efectos adversos propios de este tipo de entornos como la pérdida de visión directa o los caminos múltiples en el recorrido de la señal. Los métodos de huella radioeléctrica, más conocidos bajo el término inglés “fingerprinting”, que se basan en la comparación de los valores de potencia de señal que se están recibiendo en el momento de llevar a cabo el posicionamiento por un terminal móvil, frente a los valores registrados en un mapa radio de potencias, elaborado durante una fase inicial de calibración, aparecen como los mejores de entre los posibles para los escenarios de interior. Sin embargo, estos sistemas se ven también afectados por otros problemas, como por ejemplo los importantes trabajos a realizar para ponerlos en marcha, y la variabilidad del canal. Frente a ellos, en el presente trabajo se presentan dos contribuciones originales para mejorar los sistemas basados en los métodos fingerprinting. La primera de esas contribuciones describe un método para determinar, de manera sencilla, las características básicas del sistema a nivel del número de muestras necesarias para crear el mapa radio de la huella radioeléctrica de referencia, junto al número mínimo de emisores de radiofrecuencia que habrá que desplegar; todo ello, a partir de unos requerimientos iniciales relacionados con el error y la precisión buscados en el posicionamiento a realizar, a los que uniremos los datos correspondientes a las dimensiones y realidad física del entorno. De esa forma, se establecen unas pautas iniciales a la hora de dimensionar el sistema, y se combaten los efectos negativos que, sobre el coste o el rendimiento del sistema en su conjunto, son debidos a un despliegue ineficiente de los emisores de radiofrecuencia y de los puntos de captura de su huella. La segunda contribución incrementa la precisión resultante del sistema en tiempo real, gracias a una técnica de recalibración automática del mapa radio de potencias. Esta técnica tiene en cuenta las medidas reportadas continuamente por unos pocos puntos de referencia estáticos, estratégicamente distribuidos en el entorno, para recalcular y actualizar las potencias registradas en el mapa radio. Un beneficio adicional a nivel operativo de la citada técnica, es la prolongación del tiempo de usabilidad fiable del sistema, bajando la frecuencia en la que se requiere volver a capturar el mapa radio de potencias completo. Las mejoras anteriormente citadas serán de aplicación directa en la mejora de los mecanismos de posicionamiento en interiores basados en la infraestructura inalámbrica de comunicaciones de voz y datos. A partir de ahí, esa mejora será extensible y de aplicabilidad sobre los servicios de localización (conocimiento personal del lugar donde uno mismo se encuentra), monitorización (conocimiento por terceros del citado lugar) y seguimiento (monitorización prolongada en el tiempo), ya que todos ellas toman como base un correcto posicionamiento para un adecuado desempeño. ABSTRACT To find the position where a mobile is located with good accuracy, when it is immersed in an indoor environment (shopping centers, office buildings, airports, stations, tunnels, etc.), is the cornerstone on which a large number of applications and services are supported. Many of these services are already available in outdoor environments, although the indoor environments are suitable for other services that are specific for it. That number, however, could be significantly higher than now, if an expensive infrastructure were not required to perform the positioning service with adequate precision, for each one of the hypothetical services. Or, equally, whether that infrastructure may have other different uses beyond the ones associated with positioning. The usability of the same infrastructure for purposes other than positioning could give the opportunity of having it already available in the different locations, because it was previously deployed for these other uses; or facilitate its deployment, because the cost of that operation would offer a higher return on usability for the deployer. Wireless technologies based on radio communications, already in use for voice and data communications (mobile, WLAN, etc), meet the requirement of additional usability and, therefore, could facilitate the growth of applications and services based on positioning, in the case of being able to use it. However, determining the position with the appropriate degree of accuracy using these technologies is a major challenge today. This paper provides significant advances in this field. Along this work, a study about the main algorithms and auxiliar techniques related with indoor positioning will be initially carried out. The review will be focused in those that are suitable to be used with both last generation mobile technologies and WLAN environments. By doing this, it is tried to highlight the advantages and disadvantages of each one of these algorithms, having as final motivation their applicability both in the world of 3G and 4G mobile networks (especially in femtocells and small-cells of LTE) and in the WLAN world; and having always in mind that the final aim is to use it in indoor environments. The main conclusion of that review is that triangulation techniques, commonly used for localization in outdoor environments, are useless in indoor environments due to adverse effects of such environments as loss of sight or multipaths. Triangulation techniques used for external locations are useless due to adverse effects like the lack of line of sight or multipath. Fingerprinting methods, based on the comparison of Received Signal Strength values measured by the mobile phone with a radio map of RSSI Recorded during the calibration phase, arise as the best methods for indoor scenarios. However, these systems are also affected by other problems, for example the important load of tasks to be done to have the system ready to work, and the variability of the channel. In front of them, in this paper we present two original contributions to improve the fingerprinting methods based systems. The first one of these contributions describes a method for find, in a simple way, the basic characteristics of the system at the level of the number of samples needed to create the radio map inside the referenced fingerprint, and also by the minimum number of radio frequency emitters that are needed to be deployed; and both of them coming from some initial requirements for the system related to the error and accuracy in positioning wanted to have, which it will be joined the data corresponding to the dimensions and physical reality of the environment. Thus, some initial guidelines when dimensioning the system will be in place, and the negative effects into the cost or into the performance of the whole system, due to an inefficient deployment of the radio frequency emitters and of the radio map capture points, will be minimized. The second contribution increases the resulting accuracy of the system when working in real time, thanks to a technique of automatic recalibration of the power measurements stored in the radio map. This technique takes into account the continuous measures reported by a few static reference points, strategically distributed in the environment, to recalculate and update the measurements stored into the map radio. An additional benefit at operational level of such technique, is the extension of the reliable time of the system, decreasing the periodicity required to recapture the radio map within full measurements. The above mentioned improvements are directly applicable to improve indoor positioning mechanisms based on voice and data wireless communications infrastructure. From there, that improvement will be also extensible and applicable to location services (personal knowledge of the location where oneself is), monitoring (knowledge by other people of your location) and monitoring (prolonged monitoring over time) as all of them are based in a correct positioning for proper performance.

LiDAR remote sensing applied to forest resources assessment

Disponer de información precisa y actualizada de inventario forestal es una pieza clave para mejorar la gestión forestal sostenible y para proponer y evaluar políticas de conservación de bosques que permitan la reducción de emisiones de carbono debidas a la deforestación y degradación forestal (REDD). En este sentido, la tecnología LiDAR ha demostrado ser una herramienta perfecta para caracterizar y estimar de forma continua y en áreas extensas la estructura del bosque y las principales variables de inventario forestal. Variables como la biomasa, el número de pies, el volumen de madera, la altura dominante, el diámetro o la altura media son estimadas con una calidad comparable a los inventarios tradicionales de campo. La presente tesis se centra en analizar la aplicación de los denominados métodos de masa de inventario forestal con datos LIDAR bajo diferentes condiciones y características de masa forestal (bosque templados puros y mixtos) y utilizando diferentes bases de datos LiDAR (información proveniente de vuelo nacionales e información capturada de forma específica). Como consecuencia de lo anterior, se profundiza en la generación de inventarios forestales continuos con LiDAR en grandes áreas. Los métodos de masa se basan en la búsqueda de relaciones estadísticas entre variables predictoras derivadas de la nube de puntos LiDAR y las variables de inventario forestal medidas en campo con el objeto de generar una cartografía continua de inventario forestal. El rápido desarrollo de esta tecnología en los últimos años ha llevado a muchos países a implantar programas nacionales de captura de información LiDAR aerotransportada. Estos vuelos nacionales no están pensados ni diseñados para fines forestales por lo que es necesaria la evaluación de la validez de esta información LiDAR para la descripción de la estructura del bosque y la medición de variables forestales. Esta información podría suponer una drástica reducción de costes en la generación de información continua de alta resolución de inventario forestal. En el capítulo 2 se evalúa la estimación de variables forestales a partir de la información LiDAR capturada en el marco del Plan Nacional de Ortofotografía Aérea (PNOA-LiDAR) en España. Para ello se compara un vuelo específico diseñado para inventario forestal con la información de la misma zona capturada dentro del PNOA-LiDAR. El caso de estudio muestra cómo el ángulo de escaneo, la pendiente y orientación del terreno afectan de forma estadísticamente significativa, aunque con pequeñas diferencias, a la estimación de biomasa y variables de estructura forestal derivadas del LiDAR. La cobertura de copas resultó más afectada por estos factores que los percentiles de alturas. Considerando toda la zona de estudio, la estimación de la biomasa con ambas bases de datos no presentó diferencias estadísticamente significativas. Las simulaciones realizadas muestran que las diferencias medias en la estimación de biomasa entre un vuelo específico y el vuelo nacional podrán superar el 4% en áreas abruptas, con ángulos de escaneo altos y cuando la pendiente de la ladera no esté orientada hacia la línea de escaneo. En el capítulo 3 se desarrolla un estudio en masas mixtas y puras de pino silvestre y haya, con un enfoque multi-fuente empleando toda la información disponible (vuelos LiDAR nacionales de baja densidad de puntos, imágenes satelitales Landsat y parcelas permanentes del inventario forestal nacional español). Se concluye que este enfoque multi-fuente es adecuado para realizar inventarios forestales continuos de alta resolución en grandes superficies. Los errores obtenidos en la fase de ajuste y de validación de los modelos de área basimétrica y volumen son similares a los registrados por otros autores (usando un vuelo específico y parcelas de campo específicas). Se observan errores mayores en la variable número de pies que los encontrados en la literatura, que pueden ser explicados por la influencia de la metodología de parcelas de radio variable en esta variable. En los capítulos 4 y 5 se evalúan los métodos de masa para estimar biomasa y densidad de carbono en bosques tropicales. Para ello se trabaja con datos del Parque Nacional Volcán Poás (Costa Rica) en dos situaciones diferentes: i) se dispone de una cobertura completa LiDAR del área de estudio (capitulo 4) y ii) la cobertura LiDAR completa no es técnica o económicamente posible y se combina una cobertura incompleta de LiDAR con imágenes Landsat e información auxiliar para la estimación de biomasa y carbono (capitulo 5). En el capítulo 4 se valida un modelo LiDAR general de estimación de biomasa aérea en bosques tropicales y se compara con los resultados obtenidos con un modelo ajustado de forma específica para el área de estudio. Ambos modelos están basados en la variable altura media de copas (TCH por sus siglas en inglés) derivada del modelo digital LiDAR de altura de la vegetación. Los resultados en el área de estudio muestran que el modelo general es una alternativa fiable al ajuste de modelos específicos y que la biomasa aérea puede ser estimada en una nueva zona midiendo en campo únicamente la variable área basimétrica (BA). Para mejorar la aplicación de esta metodología es necesario definir en futuros trabajos procedimientos adecuados de medición de la variable área basimétrica en campo (localización, tamaño y forma de las parcelas de campo). La relación entre la altura media de copas del LiDAR y el área basimétrica (Coeficiente de Stock) obtenida en el área de estudio varía localmente. Por tanto es necesario contar con más información de campo para caracterizar la variabilidad del Coeficiente de Stock entre zonas de vida y si estrategias como la estratificación pueden reducir los errores en la estimación de biomasa y carbono en bosques tropicales. En el capítulo 5 se concluye que la combinación de una muestra sistemática de información LiDAR con una cobertura completa de imagen satelital de moderada resolución (e información auxiliar) es una alternativa efectiva para la realización de inventarios continuos en bosques tropicales. Esta metodología permite estimar altura de la vegetación, biomasa y carbono en grandes zonas donde la captura de una cobertura completa de LiDAR y la realización de un gran volumen de trabajo de campo es económica o/y técnicamente inviable. Las alternativas examinadas para la predicción de biomasa a partir de imágenes Landsat muestran una ligera disminución del coeficiente de determinación y un pequeño aumento del RMSE cuando la cobertura de LiDAR es reducida de forma considerable. Los resultados indican que la altura de la vegetación, la biomasa y la densidad de carbono pueden ser estimadas en bosques tropicales de forma adecuada usando coberturas de LIDAR bajas (entre el 5% y el 20% del área de estudio). ABSTRACT The availability of accurate and updated forest data is essential for improving sustainable forest management, promoting forest conservation policies and reducing carbon emissions from deforestation and forest degradation (REDD). In this sense, LiDAR technology proves to be a clear-cut tool for characterizing forest structure in large areas and assessing main forest-stand variables. Forest variables such as biomass, stem volume, basal area, mean diameter, mean height, dominant height, and stem number can be thus predicted with better or comparable quality than with costly traditional field inventories. In this thesis, it is analysed the potential of LiDAR technology for the estimation of plot-level forest variables under a range of conditions (conifer & broadleaf temperate forests and tropical forests) and different LiDAR capture characteristics (nationwide LiDAR information vs. specific forest LiDAR data). This study evaluates the application of LiDAR-based plot-level methods in large areas. These methods are based on statistical relationships between predictor variables (derived from airborne data) and field-measured variables to generate wall to wall forest inventories. The fast development of this technology in recent years has led to an increasing availability of national LiDAR datasets, usually developed for multiple purposes throughout an expanding number of countries and regions. The evaluation of the validity of nationwide LiDAR databases (not designed specifically for forest purposes) is needed and presents a great opportunity for substantially reducing the costs of forest inventories. In chapter 2, the suitability of Spanish nationwide LiDAR flight (PNOA) to estimate forest variables is analyzed and compared to a specifically forest designed LiDAR flight. This study case shows that scan angle, terrain slope and aspect significantly affect the assessment of most of the LiDAR-derived forest variables and biomass estimation. Especially, the estimation of canopy cover is more affected than height percentiles. Considering the entire study area, biomass estimations from both databases do not show significant differences. Simulations show that differences in biomass could be larger (more than 4%) only in particular situations, such as steep areas when the slopes are non-oriented towards the scan lines and the scan angles are larger than 15º. In chapter 3, a multi-source approach is developed, integrating available databases such as nationwide LiDAR flights, Landsat imagery and permanent field plots from SNFI, with good resultos in the generation of wall to wall forest inventories. Volume and basal area errors are similar to those obtained by other authors (using specific LiDAR flights and field plots) for the same species. Errors in the estimation of stem number are larger than literature values as a consequence of the great influence that variable-radius plots, as used in SNFI, have on this variable. In chapters 4 and 5 wall to wall plot-level methodologies to estimate aboveground biomass and carbon density in tropical forest are evaluated. The study area is located in the Poas Volcano National Park (Costa Rica) and two different situations are analyzed: i) available complete LiDAR coverage (chapter 4) and ii) a complete LiDAR coverage is not available and wall to wall estimation is carried out combining LiDAR, Landsat and ancillary data (chapter 5). In chapter 4, a general aboveground biomass plot-level LiDAR model for tropical forest (Asner & Mascaro, 2014) is validated and a specific model for the study area is fitted. Both LiDAR plot-level models are based on the top-of-canopy height (TCH) variable that is derived from the LiDAR digital canopy model. Results show that the pantropical plot-level LiDAR methodology is a reliable alternative to the development of specific models for tropical forests and thus, aboveground biomass in a new study area could be estimated by only measuring basal area (BA). Applying this methodology, the definition of precise BA field measurement procedures (e.g. location, size and shape of the field plots) is decisive to achieve reliable results in future studies. The relation between BA and TCH (Stocking Coefficient) obtained in our study area in Costa Rica varied locally. Therefore, more field work is needed for assessing Stocking Coefficient variations between different life zones and the influence of the stratification of the study areas in tropical forests on the reduction of uncertainty. In chapter 5, the combination of systematic LiDAR information sampling and full coverage Landsat imagery (and ancillary data) prove to be an effective alternative for forest inventories in tropical areas. This methodology allows estimating wall to wall vegetation height, biomass and carbon density in large areas where full LiDAR coverage and traditional field work are technically and/or economically unfeasible. Carbon density prediction using Landsat imaginery shows a slight decrease in the determination coefficient and an increase in RMSE when harshly decreasing LiDAR coverage area. Results indicate that feasible estimates of vegetation height, biomass and carbon density can be accomplished using low LiDAR coverage areas (between 5% and 20% of the total area) in tropical locations.

Análisis, modelado e igualación no lineal de arrays paramétricos acústicos

La presente Tesis analiza y desarrolla metodología específica que permite la caracterización de sistemas de transmisión acústicos basados en el fenómeno del array paramétrico. Este tipo de estructuras es considerado como uno de los sistemas más representativos de la acústica no lineal con amplias posibilidades tecnológicas. Los arrays paramétricos aprovechan la no linealidad del medio aéreo para obtener en recepción señales en el margen sónico a partir de señales ultrasónicas en emisión. Por desgracia, este procedimiento implica que la señal transmitida y la recibida guardan una relación compleja, que incluye una fuerte ecualización así como una distorsión apreciable por el oyente. Este hecho reduce claramente la posibilidad de obtener sistemas acústicos de gran fidelidad. Hasta ahora, los esfuerzos tecnológicos dirigidos al diseño de sistemas comerciales han tratado de paliar esta falta de fidelidad mediante técnicas de preprocesado fuertemente dependientes de los modelos físicos teóricos. Estos están basados en la ecuación de propagación de onda no lineal. En esta Tesis se propone un nuevo enfoque: la obtención de una representación completa del sistema mediante series de Volterra que permita inferir un sistema de compensación computacionalmente ligero y fiable. La dificultad que entraña la correcta extracción de esta representación obliga a desarrollar una metodología completa de identificación adaptada a este tipo de estructuras. Así, a la hora de aplicar métodos de identificación se hace indispensable la determinación de ciertas características iniciales que favorezcan la parametrización del sistema. En esta Tesis se propone una metodología propia que extrae estas condiciones iniciales. Con estos datos, nos encontramos en disposición de plantear un sistema completo de identificación no lineal basado en señales pseudoaleatorias, que aumenta la fiabilidad de la descripción del sistema, posibilitando tanto la inferencia de la estructura basada en bloques subyacente, como el diseño de mecanismos de compensación adecuados. A su vez, en este escenario concreto en el que intervienen procesos de modulación, factores como el punto de trabajo o las características físicas del transductor, hacen inviables los algoritmos de caracterización habituales. Incluyendo el método de identificación propuesto. Con el fin de eliminar esta problemática se propone una serie de nuevos algoritmos de corrección que permiten la aplicación de la caracterización. Las capacidades de estos nuevos algoritmos se pondrán a prueba sobre un prototipo físico, diseñado a tal efecto. Para ello, se propondrán la metodología y los mecanismos de instrumentación necesarios para llevar a cabo el diseño, la identificación del sistema y su posible corrección, todo ello mediante técnicas de procesado digital previas al sistema de transducción. Los algoritmos se evaluarán en términos de error de modelado a partir de la señal de salida del sistema real frente a la salida sintetizada a partir del modelo estimado. Esta estrategia asegura la posibilidad de aplicar técnicas de compensación ya que éstas son sensibles a errores de estima en módulo y fase. La calidad del sistema final se evaluará en términos de fase, coloración y distorsión no lineal mediante un test propuesto a lo largo de este discurso, como paso previo a una futura evaluación subjetiva. ABSTRACT This Thesis presents a specific methodology for the characterization of acoustic transmission systems based on the parametric array phenomenon. These structures are well-known representatives of the nonlinear acoustics field and display large technological opportunities. Parametric arrays exploit the nonlinear behavior of air to obtain sonic signals at the receptors’side, which were generated within the ultrasonic range. The underlying physical process redunds in a complex relationship between the transmitted and received signals. This includes both a strong equalization and an appreciable distortion for a human listener. High fidelity, acoustic equipment based on this phenomenon is therefore difficult to design. Until recently, efforts devoted to this enterprise have focused in fidelity enhancement based on physically-informed, pre-processing schemes. These derive directly from the nonlinear form of the wave equation. However, online limited enhancement has been achieved. In this Thesis we propose a novel approach: the evaluation of a complete representation of the system through its projection onto the Volterra series, which allows the posterior inference of a computationally light and reliable compensation scheme. The main difficulty in the derivation of such representation strives from the need of a complete identification methodology, suitable for this particular type of structures. As an example, whenever identification techniques are involved, we require preliminary estimates on certain parameters that contribute to the correct parameterization of the system. In this Thesis we propose a methodology to derive such initial values from simple measures. Once these information is made available, a complete identification scheme is required for nonlinear systems based on pseudorandom signals. These contribute to the robustness and fidelity of the resulting model, and facilitate both the inference of the underlying structure, which we subdivide into a simple block-oriented construction, and the design of the corresponding compensation structure. In a scenario such as this where frequency modulations occur, one must control exogenous factors such as devices’ operation point and the physical properties of the transducer. These may conflict with the principia behind the standard identification procedures, as it is the case. With this idea in mind, the Thesis includes a series of novel correction algorithms that facilitate the application of the characterization results onto the system compensation. The proposed algorithms are tested on a prototype that was designed and built for this purpose. The methodology and instrumentation required for its design, the identification of the overall acoustic system and its correction are all based on signal processing techniques, focusing on the system front-end, i.e. prior to transduction. Results are evaluated in terms of input-output modelling error, considering a synthetic construction of the system. This criterion ensures that compensation techniques may actually be introduced, since these are highly sensible to estimation errors both on the envelope and the phase of the signals involved. Finally, the quality of the overall system will be evaluated in terms of phase, spectral color and nonlinear distortion; by means of a test protocol specifically devised for this Thesis, as a prior step for a future, subjective quality evaluation.

Validación de HP UV-LED para la detección en continuo de contaminantes en aguas continentales mediante fluorescencia

La preservación del medio ambiente, el avance en las técnicas para que el impacto de la actividad humana sobre la fauna y flora sea lo menor posible, hacen que se deban monitorizar los diversos indicadores de calidad. El presente estudio viene motivado debido a que actualmente existen sistemas de medida y control en continuo de la calidad de las aguas, al margen de los estudios de laboratorio por toma de muestras, a través de los cuales se obtienen indicadores de calidad. El desarrollo tecnológico en analizadores en continuo para la medida de fósforo, amonio, DBO y otros, hacen que cada vez se consiga un control más exhaustivo de la calidad en tiempo real. Sin embargo, la detección temprana de contaminantes que no deben encontrarse presentes en el agua, hacen que el desarrollo de sensores de detección de estos contaminantes sea de gran utilidad. A este respecto, las técnicas mediante fluorescencia presentan enormes ventajas, ya que no existe contacto directo con la muestra, reduciéndose el desgaste y alargando el tiempo entre mantenimientos, como se ha comprobado en numerosos desarrollos con tecnología láser. Para la producir fluorescencia, tradicionalmente se vienen utilizando en el laboratorio principalmente lámparas de gas y monocromadores. Los nuevos LED de alta potencia en el espectro ultravioleta son una alternativa muy interesante que además puede ser aplicada en los mencionados sistemas de medición en continuo. En este trabajo se realiza un estudio de viabilidad de estos dispositivos como fuentes de excitación para la producción de fluorescencia tomando como contaminantes los hidrocarburos. El funcionamiento en estaciones en continuo hace que se tenga que realizar además ensayos de vida acelerados, así como estudios de modos de trabajo. Al respecto de la fluorescencia producida, se estudia la influencia de factores que pueden afectar a las medidas, tales como la temperatura. El estudio del espectro y su análisis para la identificación del contaminante es otro de los puntos desarrollados en este trabajo. Por último, y dado que la monitorización se realiza en modo continuo, es necesario un sistema de comunicaciones compacto y fiable: en este apartado se analizan los metamateriales como solución tecnológica, ya que se adapta perfectamente a la filosofía de estas estaciones de medición. ABSTRACT Currently the monitoring of quality indicators is a need to preserve the environment and minimize the impact of human activity on the fauna and flora. Currently there are measuring systems and continuous monitoring of water quality, regardless of sampling laboratory studies, through which quality indicators are obtained. Technological development in continuous analyzers for the measurement of phosphorus, ammonia, BOD and others increasingly make a more comprehensive real-time quality control is achieved. However, early detection of contaminants that should not be present in the water, make the development of sensors for detecting these contaminants is very useful. In this regard, fluorescence techniques have huge advantages, since there is no direct contact with the sample, reducing wear and extending the time between maintenance, as has been demonstrated in numerous developments in laser technology. To produce fluoresce, traditionally are being used mainly gas lamps and monochromators at the laboratory. The new high-power LEDs in the ultraviolet spectrum are a very interesting alternative that can also be applied in the above continuous measurement systems. In this paper a viability study of these devices as excitation sources to produce fluorescence using hydrocarbon as contaminants is performed. The stations in continuous operation makes it necessary to also perform accelerated life tests and studies operating modes. In regard to the fluorescence produced, the influence of factors that may affect the measurements, such as temperature is studied. The study of the spectrum and analysis to identify the contaminant is another of the points developed in this work. Finally, since the monitoring is carried out in continuous mode, a compact and reliable communication is necessary: in this section metamaterials as a technological solution is analyzed since it fits perfectly with the philosophy of these measuring stations.

Sistemas de ayuda al diagnóstico y a la terapia funcional en enfermedades neurodegenerativas

El incremento de la esperanza de vida en los países desarrollados (más de 80 años en 2013), está suponiendo un crecimiento considerable en la incidencia y prevalencia de enfermedades discapacitantes, que si bien pueden aparecer a edades tempranas, son más frecuentes en la tercera edad, o en sus inmediaciones. Enfermedades neuro-degenerativas que suponen un gran hándicap funcional, pues algunas de ellas están asociadas a movimientos involuntarios de determinadas partes del cuerpo, sobre todo de las extremidades. Tareas cotidianas como la ingesta de alimento, vestirse, escribir, interactuar con el ordenador, etc… pueden llegar a ser grandes retos para las personas que las padecen. El diagnóstico precoz y certero resulta fundamental para la prescripción de la terapia o tratamiento óptimo. Teniendo en cuenta incluso que en muchos casos, por desgracia la mayoría, sólo se puede actuar para mitigar los síntomas, y no para sanarlos, al menos de momento. Aun así, acertar de manera temprana en el diagnóstico supone proporcionar al enfermo una mayor calidad de vida durante mucho más tiempo, por lo cual el esfuerzo merece, y mucho, la pena. Los enfermos de Párkinson y de temblor esencial suponen un porcentaje importante de la casuística clínica en los trastornos del movimiento que impiden llevar una vida normal, que producen una discapacidad física y una no menos importante exclusión social. Las vías de tratamiento son dispares de ahí que sea crítico acertar en el diagnóstico lo antes posible. Hasta la actualidad, los profesionales y expertos en medicina, utilizan unas escalas cualitativas para diferenciar la patología y su grado de afectación. Dichas escalas también se utilizan para efectuar un seguimiento clínico y registrar la historia del paciente. En esta tesis se propone una serie de métodos de análisis y de identificación/clasificación de los tipos de temblor asociados a la enfermedad de Párkinson y el temblor esencial. Empleando técnicas de inteligencia artificial basadas en clasificadores inteligentes: redes neuronales (MLP y LVQ) y máquinas de soporte vectorial (SVM), a partir del desarrollo e implantación de un sistema para la medida y análisis objetiva del temblor: DIMETER. Dicho sistema además de ser una herramienta eficaz para la ayuda al diagnóstico, presenta también las capacidades necesarias para proporcionar un seguimiento riguroso y fiable de la evolución de cada paciente. ABSTRACT The increase in life expectancy in developed countries in more than 80 years (data belongs to 2013), is assuming considerable growth in the incidence and prevalence of disabling diseases. Although they may appear at an early age, they are more common in the elderly ages or in its vicinity. Nuero-degenerative diseases that are a major functional handicap, as some of them are associated with involuntary movements of certain body parts, especially of the limbs. Everyday tasks such as food intake, dressing, writing, interact with the computer, etc ... can become large debris for people who suffer. Early and accurate diagnosis is crucial for prescribing optimal therapy or treatment. Even taking into account that in many cases, unfortunately the majority, can only act to mitigate the symptoms, not to cure them, at least for now. Nevertheless, early diagnosis may provide the patient a better quality of life for much longer time, so the effort is worth, and much, grief. Sufferers of Parkinson's and essential tremor represent a significant percentage of clinical casuistry in movement disorders that prevent a normal life, leading to physical disability and not least social exclusion. There are various treatment methods, which makes it necessary the immediate diagnosis. Up to date, professionals and medical experts, use a qualitative scale to differentiate the disease and degree of involvement. Therefore, those scales are used in clinical follow-up. In this thesis, several methods of analysis and identification / classification of types of tremor associated with Parkinson's disease and essential tremor are proposed. Using artificial intelligence techniques based on intelligent classification: neural networks (MLP and LVQ) and support vector machines (SVM), starting from the development and implementation of a system for measuring and objective analysis of the tremor: DIMETER. This system besides being an effective tool to aid diagnosis, it also has the necessary capabilities to provide a rigorous and reliable monitoring of the evolution of each patient.

Taxonomía y caracterización del arbusto forrajero Atriplex halimus L.

El objetivo de esta Tesis Doctoral ha sido la revisión de la taxonomía infraespecífica de la especie Atriplex halimus L., la delimitación del área geográfica que ocupan sus taxones y su caracterización morfológica. En base a las citas bibliográficas y material de herbario consultado se ha trazado un mapa de distribución de la especie, que delimita su distribución como planta espontánea a la Cuenca Mediterránea, desierto de Siria, Macizo de Hoggar y Sahara Occidental. El resto debe considerarse poblaciones naturalizadas o cultivadas. La taxonomía actual, basada en relaciones filogenéticas, nos indica que lo adecuado es separar las poblaciones de A. halimus en dos clados que se caracterizan por tener diferente nivel de ploidía. Para definir el área de distribución de los clados diploide y tetraploide se ha determinado la ploidía de 19 poblaciones en el espacio comprendido entre las dos áreas de distribución conocidas. Y como resultado el área de distribución de ambos clados puede definirse trazando una recta que une el estrecho de Gibraltar con Estambul, las poblaciones al sur de esta línea son tetraploides a las que habría que añadir las poblaciones del este y sur de la isla de Cerdeña que son también tetraploides. Se han determinado cuales son los parámetros morfológicos que permiten caracterizar ambos clados, para lo cual se han estudiado 52 poblaciones que comprenden toda el área de distribución de la especie. Se han evaluado caracteres cualitativos y cuantitativos de valvas fructíferas, plántulas cultivadas en cámara y arbustos de 2 a 4 años situados en la misma parcela. El único elemento fiable para determinar la pertenencia al clado diploide o tetraploide es la medición del nivel de ploidía. Ninguno de los caracteres morfológicos estudiados es plenamente fiable para determinar si los individuos de una población pertenecen al clado diploide o tetraploide a pesar de realizarse las mediciones sobre poblaciones cultivadas en igualdad de condiciones de suelo y clima, y en condiciones muy favorables para la expresión de caracteres cualitativos y cuantitativos. Los caracteres cualitativos observados que guardan una mayor correlación con el nivel de ploidía son: la base de las hojas de las plántulas, el porte del arbusto, la ramificación, el ángulo de inserción de las ramas y la longitud de la inflorescencia. Respecto a los datos cuantitativos los más relevantes han sido, la longitud de los cotiledones y la longitud del limbo de la 1ª, 3ª o 5ª hoja de la plántula o la razón entre este valor y la anchura del limbo. ABSTRACT The aim of this Thesis has been in the revision of the infraspecies taxa of Atriplex halimus L., the delimitation of their natural growing area and their morphological characterization. Using references to relevant literature and consultation of herbariums, it has been possible to map the distribution of the species as spontaneous plant in the Mediterranean Basin, Desert of Syria, Massif of Hoggar and Western Sahara. The remaining populations must be considered as naturalized or cultivated. The current taxonomy, based on phylogenetic relationships, indicates that A.halimus should be divided in two clades. These clades are characterized by their different ploidy level. Samples of nineteen populations collected from localizations between the defined distribution area of both clades have been characterized for their ploidy level. As a result, it can be seen that the line of clade separation can be trazed from Istambul tothe Strait of Gibraltar, where the populations north to this line are classified as diploid and those south of the line are clsassified as tetraploid. The only exception of this rule is in Sardinia, hwre populations to the West and South of the Island are also classified as tetraploid. The morphological parameters, which include the characterization of both clades, have been defined, and the study of fifty-two populations from the entire distribution area according to these parameters have been analysed. The qualitative and quantitative characteristics of the fruit valves have been evaluated, with seedlings grown in a culture chamber, and 2-4 years old shrubs grown in the same orchard. None of the morphological characters under study are fully reliable to determine whether individuals in a population belong to clade diploid or tetraploid. This is despite measurements being made on populations grown under the same conditions of soil and climate, and in very favorable conditions for the expression of qualitative and quantitative traits. The qualitative characteristics that show a higher correlation with the ploidy level are: the base of the leaves of the seedlings, the bearing bush, the branch arrangement, the angle of insertion of the branches and the inflorescence length. With regards to the quantitative data the most relevant were the length of the cotyledons and the length of the blade of the 1st, 3rd or 5th leaves of the seedling, and the ratio of any of this values and the leaf blade width. The only reliable way to determine the diploid or tetraploid clade membership is the measurement of ploidy level.

Estudio acústico de la sala 8 de los cines Kinépolis de Ciudad de la Imagen

Este proyecto consistira en la realization de un estudio aciistico sobre la sala 8 de los cines Kinepolis de Ciudad de la Imagen, que dispone de 408 butacas. Los cines Kinepolis es uno de los mayores complejos multisala de Europa. Cuenta con mas de 9.200 butacas en total distribuidas en 25 salas a las que se accede mediante dos grandes pasillos conectados por el hall. En 1998, ano de su apertura, el complejo recibio el Record Guinness a la sala cinematografica mas grande del mundo, que dispone de 996 butacas. El objetivo de este proyecto es conseguir caracterizar acusticamente una sala de cine a traves de la medicion de parametros acusticos de la sala y de un modelo virtual de la misma. Para llevar a cabo el proyecto, primero se van a realizar tanto una medicion geometrica como acustica de la sala mediante el sistema de medida DIRAC. Los resultados de estas mediciones nos serviran para construir y validar un modelo virtual de la sala real con el software de simulacion EASE. La medicion acustica se va a realizar con el sistema de medicion DIRAC. Este sistema nos dara information sobre una amplia variedad de parametros acusticos. En este proyecto no se va a trabajar con todos ellos, solo con los mas significativos. Estos se describen a continuacion en la introduccion teorica. La medicion geometrica nos va a servir para construir un modelo virtual que tenga las mismas dimensiones que la sala original. Esta medicion la realizaremos mediante un medidor laser y una cinta metrica. Una vez construido el modelo virtual, se procedera a su validacion. Este proceso se realiza ajustando el tiempo de reverberacion del modelo mediante la introduccion de distintos materiales acusticos en las superficies del mismo, de manera que, variando la absorcion de la sala, el tiempo de reverberacion promedio del modelo se asemeje lo mas posible al medido en la sala real. Este proceso tiene como objetivo comprobar que el modelo virtual tiene un comportamiento acustico similar al de la sala real. Es necesario validar adecuadamente el modelo para que las comparaciones y conclusiones sean fiables. Por ultimo, tras la simulacion acustica del modelo, se compararan los resultados simulados con los medidos en la sala. En este proceso se contrastaran algunos de los parametros que guardan relation con el tiempo de reverberacion. De esta manera se verificara si el tiempo de reverberacion es o no un parametro acustico fiable para la validacion de un modelo virtual de una sala de cine. Anteriormente se han realizado proyectos iguales de otras salas de diferente tamano de Kinepolis. El objetivo de realizar el mismo estudio en distintas salas, es comprobar si el tamano de la sala influye en la validacion de los modelos virtuales mediante el tiempo de reverberacion. ABSTRACT. This Project consists on the development of an acoustic research of the movie theater 8 of the Kinepolis complex in Ciudad de la Imagen, Madrid. This room has 408 spots. Kinepolis is one of the biggest multiplex complexes in Europe. It has 9,200 locations disposed in 25 rooms. There are two large corridors which give access to all of theaters. In the middle of the structure, there is the main hall that connects these corridors. In 1998, at the time when the complex was open, it was awarded with the Record Guinness for the biggest theater in the world, which has 996 locations. The target of this project is to successfully characterize the acoustics of a movie theater through reverberation time and a virtual model. In order to reach this goal, in the first place, we are going to perform both, an acoustic and a geometric measurement of the room using DIRAC measurement system. The results of these measures will allow us to build and validate a virtual model of the room, using the simulation software EASE. We are going to use the DIRAC system in order to accomplish the acoustic measure. This operation gives us a huge variety of acoustic parameters. Not all of these are going to be used for this research, only the most significant ones. These are described in the theoretical introduction. The geometric measure is essential to help us to build the virtual model, because the model has to be exactly equal as the real room. This measurement will be performed with an electronic distance meter and a measuring tape. Once the virtual model is finished, it will be proved. This validation process will be realized by adjusting the reverberation time in the model. We will change the walls materials, therefore, the overall absorption of the room will change. We want the model reverberation time resemble to the real one. This practice is going to ensure that the model acoustic performance is close to the real one. In addition, it has to be successfully validate of we want the future comparisons to be reliable. Finally, after the model virtual simulation, we will compare the simulated results with the measure in the room. In this process, we will compare not only the reverberation time, but others parameters that keep relation with the reverberation time. We will verify this way, if the reverberation time is or is not an appropriate acoustic parameter to validate a virtual model of a movie theater. There have been done others similar acoustic researches in different theaters with different sizes. The aim of performing similar researches in different rooms is to determine if the size of the room infers in the validation process.

Crecimiento forestal en el bosque tropical de montaña: efectos de la diversidad florística y de la manipulación de nutrientes

El bosque tropical de montaña, es considerado zona de mega diversidad y de alto grado de endemismo, por las diferentes zonas ecológicas que presenta. Durante las últimas décadas estos bosques han recibido mayor atención por parte de investigadores, sin embargo, aún existe poca investigación en determinar cuáles son las respuestas de los bosques a los cambios ambientales a los que son sometidos. Estos bosques están sufriendo serias amenazas como pérdida de cobertura vegetal y cambios en los ciclos de nutrientes. El trabajo se dividió en cuatro objetivos específicos. i) Caracterización y análisis de patrones altitudinales de la riqueza de especies en el bosque tropical de montaña en el sur del Ecuador; con el fin de conocer cómo varía la diversidad de especies riqueza lo largo de un gradiente altitudinal. ii) Conocer los patrones espaciales del crecimiento en tres remanentes boscosos de un bosque tropical de montaña para determinar cómo la vecindad y la semejanza funcional de ésta influyen en el crecimiento forestal. iii) Conocer los efectos de la fertilización en el crecimiento diamétrico de especies arbóreas, en el bosque tropical de montaña; se analizó cómo reaccionan los árboles a la adición de nutrientes N y P en tres tipos de bosque. iv) Saber la respuesta de la comunidad de árboles a la adición de nutrientes en el bosque montano andino; este objetivo se basó con el supuesto de la deficiencia de tres tipos de nutrientes N, P y Ca, en esta formación boscosa y cómo reaccionan los árboles a la adición de nutrientes. El presente trabajo se llevó a cabo, en el bosque tropical de montaña que se encuentra localizada en la parte adyacente del Parque Nacional Podocarpus (PNP) en la cordillera del Consuelo, forma parte de la cadena oriental de los Andes del sur del Ecuador El trabajo de desarrollo entre los años 2008 y 2014. Para abordar el primer objetivo se establecieron 54 parcelas ubicadas aleatoriamente a lo largo de un gradiente altitudinal (3 niveles de altitud) y se e midieron e identificaron todos los individuos mayores a 5cm de DAP. Se construyó una filogenia con Phylocom y se calcularon diferentes componentes de diversidad para cada parcela ( riqueza taxonómica, diversidad filogenética y edad media de las especies). Ajustando modelos lineares se contrastó el efecto de la altitud sobre dichos componentes y se vio que la riqueza taxonómica y la edad media de las especies aumentaron con la altitud, en sentido contrario a las predicciones de la "hipótesis del conservadurismo tropical" (Tropical Conservatism Hypothesis). Para abordar el segundo objetivo se realizó una remedición de todos los árboles cartografiados en tres parcelas permanentes de alrededor de 5000 m2 cada una, representativas de tres estados diferentes de la sucesión del bosque montano. A partir de las coordenadas y de los datos de registrados, y empleando diferentes funciones de correlación de marca se analizó la distribución espacial del tamaño y del crecimiento relativo y del tamaño. Se constató que mientras que el tamaño de los árboles presentó una correlación espacial negativa, el crecimiento presentó correlación espacial positiva, en ambos casos a distancias cortas. El rango y la magnitud de ambas correlaciones aumentaron al avanzar la sucesión. La distribución espacial del crecimiento mostró una correlación negativa con la distribución espacial de tamaños. Por otro lado, la distribución espacial del crecimiento mostró una correlación negativa para árboles semejantes funcionalmente y positiva cuando se calculó entre árboles con diferente estrategia funcional. En conjunto, los resultados obtenidos señalan un aumento de la importancia de procesos competitivos y una mayor estructuración espacial del crecimiento y de la distribución de tamaños al avanzar la sucesión. Para el tercer y cuarto objetivo se instalaron 52 parcelas distribuidas en bloques donde se fertilizaron dos veces al año durante 6,4 años, se identificaron todos los individuos mayores a 10 cm de DAP, y se midió el crecimiento diamétrico durante estos años Con la adición de nutrientes realizada a los diferentes tipos de bosque en la gradiente altitudinal, encontramos que el efecto sobre el crecimiento diamétrico en la comunidad varia con el rango altitudinal, y el tipo de nutriente, analizando a nivel de las especies, en la mayoría de los casos las especies comunes no tuvieron cambios significativos a la adición de nutrientes. Los resultados de este estudio aportan nuevas evidencias para el entendimiento de la diversidad, estructura y dinámica de los bosques tropicales de montaña. ABSTRACT The montane tropical forest is considered a megadiverse habitat that harbor an enormous degree of endemism. This is mainly due to the high degree of environmental heterogeneity found and the presence of different well defined ecological areas. These forests have received more attention during the last decades, however, the information regarding the responses of these forests to environmental change, is still scarce. These forests are seriously endangered and are suffering serious threats, such as loss of vegetative cover, changes in the nutrient cycles. The work was divided in four specific objectives: i) Characterization and analysis of the species richness altitudinal patterns in the montane tropical forest of south Ecuador. Specifically, how species diversity changes along altitudinal gradients. ii) Exploring the spatial patterns of tree growth in three remnants of a montane tropical forest, and analyze how tree neighborhood and functional similarity among trees influence tree growth. Tropical Conservatism Hypothesis iii) Understanding the effects of fertilization in arboreal species growth (increase in diameter) of the montane tropical forest. Specifically we studied the effects of P and N addition on three different forests across an altitudinal gradient. iv) Know the response of the community of trees to the addition of nutrients in the Andean montane forest; this objective was based on the supposition of deficiency of three types of nutrients: P, N and Ca in this forest all formation and how the trees react to the addition of these nutrients. The present work was carried out in the montane tropical forest located in Bombuscaro, San Francisco and Cajanuma close to Podocarpus National Park (PNP) on Consuelo mountain range (Andean oriental range) at South of Ecuador. Field work was carried out during 2008 and 2014. To address the first objective, we randomly placed 54 plots along an altitudinal gradient. In these plots, every individual larger than 5 cm of DBH was measured and identified. A phylogeny was build with Phylocom and different diversity components (taxonomic richness, phylogenetic diversity and average species age) were computed for each plot. Linear models were used to test the effects of altitude on the diversity components. Our results showed that, contrary to the Tropical Conservatism Hypothesis, both taxonomic richness and average species age increased with altitude. To address our second objective, all mapped trees in three successional permanent plots (around ~5000 m2 each) were re-measured. Using different mark correlation functions, we analyzed the spatial distribution of tree-size and tree relative growth rate. Whereas tree size showed negative spatial correlation at fine spatial scales, relative growth rate showed positive correlation at the same scales. The range and magnitude of those correlations increased along successional stage. The spatial distribution of the relative growth rate was negatively correlated with the spatial distribution of tree sizes. Additionally, we found that the spatial correlation of the relative growth rate was negative for functionally similar trees and positive when computed for functionally different trees. In synthesis, our results point to an increase of competitive processes and strong spatial structure of relative growth rate and tree size along succession. For the third and fourth objectives, 52 plots were placed in a block design and were fertilized twice a year for 6,4 years. In these plots all the individuals with DBH > 10 cm were identified, and the diametrical growth was measured during these years. The nutrient addition at the three different altitude forests, revealed that the effect on the diametrical growth in the community varied with the altitudinal range. When analyzed at species level, the addition of nutrients was no significant in most cases. These results represent new evidences that will improved our understanding of diversity patterns and structure, and the dynamics of tropical montane forests.

Relationship between residual stresses and mechanical behavior in ceramic composites

En este trabajo, materiales de tipo alúmina/Y-TZP (ZrO2 tetragonal, estabilizada con 3 mol. % Y2O3), como sistema cerámico popular por sus mejoradas propiedades mecánicas en comparación con las cerámicas de alúmina puras, han sido estudiados en términos de propiedades mecánicas y tensiones residuales. El novedoso método de colado en cinta, consistente en el apilamiento de cintas de cerámica verde a temperatura ambiente y el uso de bajas presiones, se ha escogido para la presente investigación con el fin de poder aprovechar al máximo el futuro desarrollo de materiales laminados de alúmina-óxido de circonio. Se han determinado las propiedades de los materiales obtenidos por este nuevo método de procesamiento comparándolas con las de los materiales obtenidos por “slip casting”, con el fin de analizar si el método propuesto afecta a la microestructura y, por tanto, a las propiedades mecánicas y tensiones residuales propias de estos materiales. Para analizar la idoneidad del proceso de fabricación, utilizado para evitar la presencia de discontinuidades en las intercaras entre las láminas así como otros fenómenos que puedan interferir con las propiedades mecánicas, se estudiaron materiales cerámicos con la misma composición en cintas. Por otra parte también se analizó el efecto de la adición de óxido de circonio sobre la aparición de tensiónes residuales en cerámicas Al2O3/Y-TZP, teniendo en cuenta su notable influencia sobre las propiedades microestructurales y mecánicas de los materiales, así como el requisito de co-sinterización de capas con diferentes materiales compuestos en materiales laminados. La caracterización del material incluye la determinación de la densidad, el análisis de la microestructura, la obtención de las propiedades mecánicas (módulo de elasticidad, dureza, resistencia a la flexión y tenacidad de fractura) así como de las tensiones residuales. En combinación con otros métodos de medida tradicionales, la nanoindentación también se empleó como una técnica adicional para la medida del módulo de elasticidad y de la dureza. Por otro lado, diferentes técnicas de difracción con neutrones, tanto las basadas en longitud de onda constante (CW) como en tiempo de vuelo (TOF), han sido empleadas para la medición fiable de la deformación residual a través del grosor en muestras a granel. Las tensiones residuales fueron determinadas con elevada precisión, aplicando además métodos de análisis apropiados, como por ejemplo el refinamiento de Rietveld. Las diferentes fases en cerámicas sinterizadas, especialmente las de zirconia, se examinaron con detalle mediante el análisis de Rietveld, teniendo en cuenta el complicado polimorfismo del Óxido de Zirconio (ZrO2) así como las posibles transformaciones de fase durante el proceso de fabricación. Los efectos del contenido de Y-TZP en combinación con el nuevo método de procesamiento sobre la microestructura, el rendimiento mecánico y las tensiones residuales de los materiales estudiados (Al2O3/Y-TZP) se resumen en el presente trabajo. Finalmente, los mecanismos de endurecimiento, especialmente los relacionados con las tensiones residuales, son igualmente discutidos. In present work, Alumina/Y-TZP (tetragonal ZrO2 stabilized with 3 mol% Y2O3) materials, as an popular ceramic system with improved mechanical properties compared with the pure alumina ceramics, have been studied in terms of mechanical properties and residual stresses. The novel tape casting method, which involved the stacking of green ceramics tapes at room temperature and using low pressures, is selected for manufacturing and investigation, in order to take full advantage of the future development of alumina-zirconia laminated materials. Features of materials obtained by the new processing method are determined and compared with those of materials obtained by conventional slip casting in a plaster mold, in order to study whether the proposed method of processing affects microstructure and thereby the mechanical properties and residual stresses characteristics of materials. To analyse the adequacy of the manufacturing process used to avoid the presence of discontinuities at the interfaces between the sheets and other phenomena that interfere with the mechanical properties, ceramic materials with the same composition in tapes were investigated. Moreover, the effect of addition of zirconia on residual stress development of Al2O3/Y-TZP ceramics were taken into investigations, considering its significantly influence on the microstructure and mechanical properties of materials as well as the requirement of co-sintering of layers with different composites in laminated materials. The characterization includes density, microstructure, mechanical properties (elastic modulus, hardness, flexure strength and fracture toughness) and residual stresses. Except of the traditional measurement methods, nanoindentation technique was also used as an additional measurement of the elastic modulus and hardness. Neutron diffraction, both the constant-wavelength (CW) and time-of-flight (TOF) neutron diffraction techniques, has been used for reliable through-thickness residual strain measurement in bulk samples. Residual stresses were precisely determined combined with appropriate analysis methods, e.g. the Rietveld refinement. The phase compositions in sintered ceramics especially the ones of zirconia were accurately examined by Rietveld analysis, considering the complex polymorph of ZrO2 and the possible phase transformation during manufacturing process. Effects of Y-TZP content and the new processing method on the microstructure, mechanical performance and residual stresses were finally summarized in present studied Al2O3/Y-TZP materials. The toughening mechanisms, especially the residual stresses related toughening, were theoretically discussed.

Calibración de resistencias termométricas para la transferencia de custodia de hidrocarburos líquidos refinados

El principal objetivo del presente proyecto consiste en la creación de un laboratorio de metrología de hidrocarburos líquidos para la calibración de resistencias termométricas durante la transferencia de custodia. En estas operaciones de transferencia es muy importante llevar el control de la cantidad de hidrocarburo transferido, por ello tiene una gran relevancia la medida de la temperatura. El método más fiable de medición es el realizado por resistencias termométricas tipo Pt-100, la cual nos proporciona datos de incertidumbre aproximados a cero. Estos sistemas de medición tenderán a descalibrarse, aumentando la incertidumbre, pudiendo provocar pérdidas económicas. Para su calibración este proyecto presenta la realización de un laboratorio de metrología en el marco ISO-IEC 17025. Dicha calibración se realizará bajo el procedimiento operativo para la calibración de termómetros de resistencia de platino. Por otra parte, el equipamiento del laboratorio también tiene que ser calibrado, en este caso se hará externamente.

Modelización de la gasificación de lodos de depuradoras urbanas

La gasificación de lodos de depuración consiste en la conversión termoquímica del lodo por oxidación parcial a alta temperatura mediante un agente gasificante, que generalmente es aire, oxígeno o vapor de agua. Se trata de una tecnología de gran interés, ya que consigue reducir la masa de estos residuos y permite el aprovechamiento de los gases formados, tanto en la generación de energía térmica y/o eléctrica como en la síntesis de productos químicos orgánicos y combustibles líquidos. Debido a la complejidad de este proceso, es útil el uso de modelos que faciliten su estudio de forma fiable y a bajo coste. El presente Proyecto Fin de Carrera se centra en el diseño de un modelo adimensional de equilibrio en estado estacionario basado en la minimización de la energía libre de Gibbs. Para ello, se ha empleado el software de simulación de procesos Aspen Plus, que posee una amplia base de datos de propiedades físicas y permite gran flexibilidad en el manejo de sólidos. Para la elaboración del modelo se han asumido las hipótesis de mezcla perfecta dentro del reactor y operación isoterma. El gasificador se ha considerado de lecho fluidizado burbujeante, al permitir un buen control de la temperatura y una alta transferencia de materia y energía entre el sólido y el agente gasificante. El modelo desarrollado consta de cuatro etapas. La primera reproduce el proceso de pirólisis o descomposición térmica de los componentes del lodo en ausencia de agente gasificante. En la segunda etapa se simula que todo el nitrógeno y el azufre contenidos en el lodo se transforman en amoniaco y ácido sulfhídrico, respectivamente. En la tercera etapa se produce la gasificación en dos reactores. El primer gasificador alcanza el equilibrio químico mediante la minimización de la energía libre de Gibbs del sistema. En el segundo reactor se establece un equilibrio restringido por medio de la especificación de una aproximación de temperatura para cada reacción. Este método permite validar los resultados del modelo con datos reales. En la última etapa se separa el residuo carbonoso o char (compuesto por carbono y cenizas) del gas de salida, formado por N2, H2, CO, CO2, CH4 (supuesto como único hidrocarburo presente), NH3, H2S y H2O. Este gas debe ser depurado mediante equipos de limpieza aguas abajo. Los resultados de la simulación del modelo han sido validados frente a los valores obtenidos en ensayos previos llevados a cabo en la planta de gasificación a escala de laboratorio ubicada en el Departamento de Ingeniería Química Industrial y del Medio Ambiente de la Escuela Técnica Superior de Ingenieros Industriales de la Universidad Politécnica de Madrid. Estos resultados han mostrado muy buena concordancia con los obtenidos experimentalmente, con un error inferior al 7% en todos los parámetros analizados en el caso de gasificación con aire y menor al 13% cuando se utiliza una mezcla aire/vapor de agua como agente gasificante. Se ha realizado un análisis de sensibilidad con el fin de estudiar la influencia de las condiciones de operación (temperatura, ratio equivalente y ratio vapor/biomasa) sobre los resultados del proceso modelado (composición, producción y poder calorífico inferior de los gases, conversión de carbono y eficiencia de la gasificación). Para ello, se han llevado a cabo diferentes simulaciones modificando la temperatura de gasificación entre 750ºC y 850ºC, el ratio equivalente (ER) entre 0,2 y 0,4 y el ratio vapor/biomasa (S/B) entre 0 y 1. Como ya ocurriera con la validación del modelo, los resultados de las simulaciones bajo las distintas condiciones de gasificación se ajustan de forma satisfactoria a los valores experimentales. Se ha encontrado que un aumento en la temperatura mejora la cantidad y la calidad del gas producido y, por tanto, la eficiencia del proceso. Un incremento del ratio equivalente reduce la concentración de CO y H2 en el gas y, en consecuencia, también su poder calorífico. Sin embargo, valores bajos del ratio equivalente disminuyen la producción de gases y la conversión de carbono. La alimentación de vapor de agua en el sistema mejora todos los parámetros analizados. Por tanto, dentro del rango estudiado, las condiciones de operación que optimizan el proceso de gasificación de lodos consisten en el empleo de mezclas aire/vapor de agua como agente gasificante, una temperatura de 850ºC y un ER de 0,3.

Factores anatómicos, dendrométricos y climáticos implicados en la producción de resina de Pinus pinaster AIT. : aplicación a la mejora de los métodos de resinación

Pinus pinaster Ait. es la conífera con mayor área de distribución en la Península Ibérica y es, a día de hoy, la única especie resinada en nuestro país. La inducción del flujo de resina al exterior para su recolección a través de distintos tipos de heridas ha sido practicada desde hace miles de años por distintas culturas. En todos los casos, las técnicas desarrolladas se basan en la estimulación del característico sistema de defensa de las pináceas. En los últimos siete años se viene observando una tendencia de incremento sustancial de la superficie resinada en España, acompañada por avances tecnológicos dirigidos a la mecanización y mejora de estimulantes. El aprovechamiento resinero se perfila como un sector estratégico en la generación de empleo rural y la conservación de ecosistemas. La industria resinera demanda métodos de extracción más eficaces, una selvicultura adecuada y actualizada, y condiciones laborales de los resineros más dignas con objeto de llegar a ser competitiva en el mercado internacional. Este trabajo se centra en ampliar el conocimiento sobre el sistema de defensa de P. pinaster, concretamente sobre las estructuras y procesos que pueden afectar a la producción de resina. Se analizan las relaciones entre las características anatómicas del xilema, destacando las relacionadas con los canales resiníferos, las variables dendrométricas y dasométricas de la masa y el flujo de resina (objetivo 1). Se estudia cómo estas relaciones son moduladas por las heridas de resinación dependiendo de la técnica de resinación aplicada (objetivo 2), el clima y el balance hídrico del suelo (objetivo 3). El material vegetal, las muestras de suelo y los datos de producción de resina y climáticos usados en esta tesis han sido recogidos en tres montes de utilidad pública; MUP 101 en Armuña, MUP 108 en Melque de Cercos y MUP 117 en Nieva (en esta última solo se recogieron los datos de producciones), todos ellos pinares monoespecíficos de P. pinaster localizados en la denominada Tierra de Pinares Segoviana. En los árboles de nuestro estudio se han aplicado cuatro métodos de resinación: método de pica de corteza con estimulante y método mecanizado con estimulante, ambos en sentido ascendente y descendente. En los trabajos realizados para el análisis de la influencia de la anatomía constitutiva en la producción de resina (objetivo 1) y el efecto del clima (objetivo 3), se obtuvieron muestras del xilema de 26 árboles resinados en Melque de Cercos y Armuña y 12 árboles control sin resinar. Para caracterizar los pies estudiados, se midió la altura, diámetro normal y porcentaje de copa viva. Las muestras de tejido fueron recogidas en una zona del tronco a una distancia del límite de la herida considerada en la bibliografía como no afectada (anatomía constitutiva). Para el análisis de las alteraciones anatómicas inducidas por la herida (objetivo 2), se recogieron muestras en ocho de los individuos en los que se habían realizado los distintos métodos de resinación descritos y en cinco árboles control. Se obtuvieron ocho muestras de tejido distribuidas en la parte superior, inferior, lateral y centro de la herida de cada uno de los árboles resinados. Para establecer las diferencias en la producción de resina según el método de resinación, se analizaron las producciones de 561 árboles resinados en 2012 con estos cuatro métodos en Nieva. Los principales resultados de estos trabajos muestran que la producción de resina está ligada al volumen de canales (axiales y radiales) y a la frecuencia de canales radiales existentes en el árbol antes de efectuar ninguna herida (sistema constitutivo). De esta manera, los árboles grandes productores de resina mostraron una red de canales más densa que aquellos con producciones medias. Una vez realizada la herida de resinación, observamos una disminución del ancho del anillo de crecimiento y del tamaño medio de los canales axiales a la vez que se incrementaba la frecuencia y área ocupada por mm2 de anillo de estos canales. Estos cambios perduraron en el árbol durante al menos tres años y fueron distintos dependiendo de la localización en el entorno de la herida y del método de resinación. Las respuestas más intensas a la herida se observaron el año siguiente a la realización de la misma, en dirección axial, para las distancias más próximas al límite de la herida y para los métodos de resinación en sentido ascendente. Además, se ha constatado que como consecuencia de las heridas de resinación se produjeron cambios en la anatomía del xilema en zonas alejadas de la herida, tanto en el año de la herida como años posteriores. Es decir, se observó una respuesta sistémica del árbol. Respecto al papel del clima como regulador de la respuesta del árbol, se ha evidenciado que la temperatura, la radiación y la ETP influyeron en la producción de resina, no solo durante la campaña de resinación, sino también durante los meses anteriores. El déficit hídrico favoreció la producción y la formación de canales axiales pero, a partir de un determinado umbral, esa relación se invirtió y las lluvias estivales incrementaron la producción. Algunas de estas variables climáticas se asociaron a cambios en el tamaño y frecuencia de las estructuras secretoras, las cuales posiblemente modulan la respuesta defensiva de la planta. La dendrometría del árbol (evaluada a través del diámetro normal, altura y porcentaje de copa viva), la densidad de la masa y el tipo de suelo influyeron en el potencial de producción de resina de P. pinaster. Árboles más vigorosos, parcelas con menores densidades y suelos con más capacidad para la retención de agua y nutrientes presentaron producciones mayores. Estos trabajos se complementan en anexos con una caracterización del sistema socio-ecológico del pinar en resinación. En ese trabajo se identifican sus potenciales servicios ecosistémicos y se evalúa su grado de vinculación con el aprovechamiento resinero con objeto de valorar su funcionalidad y aproximar una valoración económica de modo que sea posible apreciar la importancia económica de los mismos. Para concluir, podemos resaltar que son necesarios más trabajos de carácter científico para avanzar en la comprensión de los procesos anatómicos y fisiológicos que regulan la secreción de resina en P. pinaster y sus interacciones con el medio. Esto permitiría avances certeros hacia el desarrollo de métodos de extracción más eficaces, una selvicultura óptima, el reconocimiento de los beneficios socio-ecológicos y económicos del aprovechamiento y, de manera general, una bibliografía amplia y fiable para la consulta y desarrollo de futuras mejoras que posibiliten la reactivación y conservación de la resinación como aprovechamiento rentable. ABSTRACT Pinus pinaster Ait. is the most widespread conifer in Spain and is now the only species tapped for its oleoresin. External induction of resin secretion, based on the defense system of Pinus trees, has been performed by humans since Classical times through various methods. The socio-economic implication of this practice in Spain justifies a new approach to improve tapping methodology and understand the effects of this activity on the tree. In the last five years, sharp increases in the price of natural resins, accompanied by technological advances directed toward mechanization, have made resin tapping a strategic activity for rural development and forest conservation. The resin industry demands more efficient tapping methods and forest management plans as a way to increase competitiveness in a global market. In this way, this work focuses on the study of the defense system of P. pinaster, with the aim to understand the effects of anatomical and physiological characteristics and environmental conditions on resin yield. The relationships between anatomical variables -with special focus on resin canals-, dendrometric and dasometric variables, and resin yield will be evaluated (objective 1). The tapping wound effects (objective 2) and the intra- and inter-annual variability of climate conditions and soil water availability influence (objective 3) on resin yield will be also studied. The plant and soil material and the resin yield and climatic data used in this thesis have been collected in stands of three public forests of P. pinaster; Armuña, Melque de Cercos and Nieva, located in Segovia (Central Spain). Trees were tapped using two different methods: mechanized or traditional tool, in both upwards and downwards direction. Twenty-six tapped trees of contrasting resin yield classes and twelve non-tapped (control) trees, growing in two locations (Armuña y Melque de Cercos) with the same climate but different stand density and soil characteristics, were selected for studying the role of tree size, xylem anatomy at distal parts aside from the tapping wound (objective 1) and climate influence (objective 3) on resin yield. Concerning the tree defenses induced by the tapping wound (objective 2), the xylem of eight trees, tapped with the two described methods in both upwards and downwards direction, were analyzed. From each tapped tree, eight cores were collected at different locations and varying distances from the tapping wound. In each core, a histological analysis was made. Growth ring width, earlywood and latewood width, and axial canal frequency, area, mean size and location were measured. The effect of the tapping method on resin yield was assessed in 561 P. pinaster tapped trees in a stand in Nieva. In tissues not affected by the tapping wound, the frequency of radial resin canals and the total volume of resin canals were related to resin yield. The frequency of radial canals and the resin yield were strongly related to tree diameter and percentage of live crown. High area of axial resin canals per mm2 was related to high yielding trees, but only in the location with higher plant density and poorer soil quality. In tapped trees, an increase in axial canal frequency and area was found during the three years following the start of tapping activity, suggesting that canal formation is a systemic induced response to wounding. The highest mean annual resin yield was found using the traditional tool in upwards direction, which also induced the highest increase in axial canal frequency and area. The lowest yield was found for mechanized tapping, which showed no differences between the upwards and downwards directions. The strongest induction of systemic induced responses in terms of resin canal frequency and area was detected one year after tapping for upwards tapping. This suggests the involvement of signaling processes that spread mainly upwards, and the importance of adaptive processes as a defense against periodic insect attacks. Intra-annual variation in resin yield was strongly correlated with temperature, solar radiation, potential evapotranspiration and soil water deficit. Inter-annual variation in resin yield and resin canal abundance were correlated with temperature and water deficit in spring, but above a certain threshold of cumulated water deficit in summer rainfall favored resin yield. Under adverse climate scenarios where resource optimization is desirable, a reduced tapping season during the warmest months (June–September) would be advisable, assuming a very small production loss relative to traditional tapping season. Similarly, in years with a rainy summer and/or dry spring, a slightly longer tapping season could be suggested, as resin yield increases after these events. Tree diameter and percentage of live crown, and radial resin canal frequency could be useful criteria for estimating resin yields in P. pinaster. Vigorous trees in lower density stands and growing up in good quality soils will be the most productive. These conclusions could be applied to improve tapping management and breeding programs. These works are complemented with socio-ecological characterization, the identification of the main ecosystem services and an assessment of the possible economic impact derived from the tapping practice. To conclude, more scientific studies are necessary for understanding the anatomical and physiological processes behind resin synthesis and their interactions with the environment. This would afford further progresses towards an extensive and reliable bibliography and improved tapping methods and optimal selvicultural guide lines.

Eficacia de un programa de tratamiento cognitivo-conductal para trastornos por estrés traumático, depresivos y de ansiedad a largo plazo en víctimas de atentados terroristas

OBJETIVO Se analiza la eficacia de la terapia cognitivo-conductual centrada en el trauma en víctimas de atentados terroristas con trastorno por estrés postraumático (TEPT) o trastornos de ansiedad o depresivos a muy largo plazo (> 5 años después del atentado). Se pretende, así, engrosar el limitado conocimiento científico sobre la eficacia de los tratamientos en esta población concreta. MÉTODO Un grupo de 120 víctimas directas e indirectas de atentados terroristas cometidos en España una media de 18,43 años atrás, que presentaban aislada o concurrentemente trastorno por estrés postraumático (TEPT; 49,2%), trastorno depresivo mayor (TDM; 40,8%) y otros trastornos de ansiedad (75,8%), fueron aleatorizados en 16 sesiones de Terapia Cognitivo-Conductual (TCC) centrada en el trauma (grupo experimental) o 16 semanas de espera (grupo control). La TCC centrada en el trauma estaba basada en el protocolo de exposición prolongada de Foa y Rothbaum (1998) pero se le añadieron técnicas específicas para el manejo de los trastornos del estado de ánimo y trastornos de ansiedad, se hizo un mayor énfasis en la psicoeducación sobre las reacciones post traumáticas y en la reestructuración cognitiva y se añadieron estrategias de motivación, regulación emocional y terapia narrativa. De las 60 personas que fueron asignadas al grupo experimental, 25 (41,6%) rechazaron la terapia y 10 (28,6%) la abandonaron una vez iniciada. RESULTADOS Las personas del grupo experimental que completaron el tratamiento (n = 25) tenían puntuaciones medias en sintomatología postraumática (PCL), depresiva (BDIL-II) y de ansiedad (BAI) significativamente inferiores a las obtenidas en el pretratamiento. Ninguna de ellas presentó en el postratamiento TEPT ni TDM y se redujo significativamente el porcentaje de otros trastornos de ansiedad. Estos resultados fueron significativamente mejores que los del grupo de control. Se obtuvieron tamaños del efecto pre/post grandes (d PCL = 1,39; d BDI-II = 1,19; d BAI = 1,20) y grandes y medios entre grupos (d PCL = 0,84; d BDI-II = 0,69; d BAI = 0,88). Un 78,3%, 91,7% y 91,7% obtuvieron una puntuación media inferior a la puntuación C para la PCL, el BDI-II y el BAI, y un 56,5%, 45,8% y 54,2% mejoraron según el índice de cambio fiable. Los beneficios se mantuvieron a los 6 meses de seguimiento y el análisis de intención de tratar (n = 35) y el análisis de todas las personas que habían completado el tratamiento (incluyendo al grupo control; n = 40) confirmaron la significación de los datos. La terapia tenía efectos mejores sobre la comorbilidad triple (TEPT, trastornos depresivos y ansiosos), pero no se encontró un efecto positivo de la adición de psicofármacos, en todo caso, peores resultados (por ejemplo, en la puntuación del BAI y la muestra de intención de tratar). CONCLUSIONES Los resultados de este estudio, comparables a los de otros estudios específicos previos, sugieren que la TCC centrada en el trauma y adaptada a víctimas del terrorismo que sufren TEPT o trastornos de ansiedad o depresivos muchos años después de los atentados terroristas, es eficaz. Aunque se recomienda mejorar las técnicas de motivación, hacer nuevos estudios con tamaños muestrales más amplios, donde se realicen análisis de componentes, análisis de la fiabilidad y adherencia al protocolo y atendiendo a diferentes grupos de víctimas y factores de riesgo, los resultados de esta tesis doctoral sustentan firmemente la idea de que la terapia cognitivo-conductual centrada en el trauma debe ser el tratamiento de elección para víctimas del terrorismo con TEPT, TDM u otros trastornos de ansiedad, incluso cuando éstos se presentan de manera comórbida y crónica.